CN103063345B - 一种检测轴销所承受径向力的方法 - Google Patents
一种检测轴销所承受径向力的方法 Download PDFInfo
- Publication number
- CN103063345B CN103063345B CN201210554820.6A CN201210554820A CN103063345B CN 103063345 B CN103063345 B CN 103063345B CN 201210554820 A CN201210554820 A CN 201210554820A CN 103063345 B CN103063345 B CN 103063345B
- Authority
- CN
- China
- Prior art keywords
- radial force
- pivot pin
- epsiv
- radial
- foil
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Force Measurement Appropriate To Specific Purposes (AREA)
- Measurement Of Force In General (AREA)
Abstract
本发明公开了一种轴销式力传感器,包括轴销和至少四个应变片,每个应变片均匀布置在轴销的同一圆周径向上。当轴销承受径向力时,应变片产生应变信号输出,通过计算获得轴销所承受的径向力大小和方向。本发明还提供了检测轴销所承受径向力的方法。本发明解决检测轴销所承受多方向复杂的随机径向力的问题,提高检测的精度和减少测量误差。
Description
技术领域
本发明涉及一种传感器,更具体地说涉及一种轴销式力传感器及检测轴销所承受径向力的方法。
背景技术
传统的测试轴销受力的方法是把机械装置中的轴和销改造成为应变式测力传感器(即轴销式传感器)。根据不同用途,方便地安装在两结构联接处,既能起到替代原有轴销的功能,又能起到测力传感器的作用,从而使整个测力系统的机械部件大大简化。这种方法能把该装置系统许多困难的测力问题变得极为简单,而又不需要改变结构和增加零件。轴销式传感器体现了电子衡器小型化、模块化、集成化、智能化的发展方向,是将承载器和传感器合二为一的一种集成结构,是测量轴承、滑轮等构件的径向载荷的专用传感器。因此,轴销式传感器在机械行业有着广泛的应用。
传统的轴销式传感器,双剪型电阻应变计粘贴在轴销的一个方向位置上,通过成惠斯通电桥进行测量,通过简单的力学公式计算获得径向力大小。因此,传统的轴销式传感器,只能检测单方向径向力,无法检测变方向径向力。也就是说,当轴销式传感器设计后,轴销式传感器只能检测一个方向的径向力,其它方向的径向力检测结果误差很大。因此,传统的轴销式传感器,对于变方向载荷的测量却无能为力,对于径向力的方向测量也无能为力。在工程机械工作过程中总是承受多方向复杂的随机载荷。因此,传统的轴销式传感器在工程机械测力方面的应用受到了极大的限制。
发明内容
有鉴于此,本发明提出一种轴销式力传感器及检测轴销所承受径向力的方法,以解决检测轴销所承受多方向复杂的随机径向力的问题,提高检测的精度和减少测量误差。
一方面,本发明提供了一种轴销式力传感器,包括轴销和至少四个应变片,每个应变片均匀布置在轴销的同一圆周径向上。
进一步地,位置相对的两个应变片为一组,每组应变片构成一个电桥电路输出,每组应变片分别代表不同方向的应变值。
进一步地,沿轴销外表面上开设有多个凹槽,应变片位于凹槽内。
进一步地,沿轴销的轴线开设有通孔,应变片位于通孔的内表面上。
另一方面还提供了一种检测轴销所承受径向力的方法,包括:
步骤1:在轴销的不同径向上布置有至少四个应变片,位置相对的两个应变片为一组,每组应变片构成一个电桥电路输出,每组应变片分别代表不同方向的应变值;
步骤2:在标定试验机上对轴销的圆周上进行加载径向力,采集与每个加载径向力对应的应变片的应变信号,
步骤3;建立预测模型;通过预测模型预测轴销所承受的径向力大小和方向。
进一步地,在步骤2与步骤3之间,还包括建立神径网络模型或者支持向量机模型;将步骤2采集的每个加载径向力对应的应变片的应变信号数据输入到神径网络模型或者支持向量机模型中进行训练;得到径向力大小及方位的预测模型。
进一步地,在步骤1中,应变片均匀布置在轴销的同一圆周径向上。
另一方面还提供另一种检测轴销所承受径向力的方法,包括:
在轴销的不同径向上布置有四个应变片,相邻两个应变片径向相隔90度,位置相对的两个应变片为一组,每组应变片分别接入电桥电路,每组应变片产生的应变信号单独输出,分别代表X、Y方向的应变值εx+εxy、εy+εyx;
在X、Y方向对轴销加载径向力,获取x、y方向加载径向力与应变信号之间的系数k1、k2;
在X、Y其它方向对轴销加载径向力,获取正交效应系数k12、k21;
根据下面公式计算x、y方向的真实分力:
由p1、p2合成总的载荷并计算力的方位角θ=arctan(p2/p1)。
进一步地,在X、Y其它方向对轴销加载径向力,该径向力在X、Y上的真实分力为P1和P2表示,实测载荷在X、Y方向的分力为F1和F2表示;根据下面公式计算正交效应系数k12、k21;
F1=p1+k12p2=k1εx+k12k2εy=k1(εx+εxy);
F2=p2+k21p1=k2εy+k21k1εx=k2(εy+εyx)。
进一步地,四个应变片均匀布置在轴销的同一圆周径向上。
本发明提出一种轴销式力传感器及检测轴销所承受径向力的方法有益效果在于:其一,可以检测轴销所承受多方向复杂的随机径向力的问题;其二,检测的精度高;其三,可以测量轴销所承受径向力的大小和径向力作用方位。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明一种轴销式传感器结构示意图;
图2为图1的A—A向示意图;
图3为图1的应变片电桥电路接线示意图;
图4为本发明另一种轴销式传感器结构示意图;
图5为图4的A—A向示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
下面结合图1至图3,对本发明的优选实施例作进一步详细说明。如图1和图2所示,一种轴销式传感器,包括轴销1和四个应变片3a、3b、3c、3d,四个应变片均匀分布在轴销1同一圆周径向上,即位于轴销1同一圆周面、同一横截面上。
四个应变片3a、3b、3c、3d通过粘贴工艺粘贴在的四个凹槽2中。凹槽2为四方形,两端采用半圆形过渡结构,该凹槽2结构可以防止轴销1弹性变形时,在凹槽2处应力集中,造成检测轴销1所承受径向力测量不准,提高了检测径向力精度。同时,四个应变片位于同一圆周径向上,可以提高检测轴销1所承受径向力精度,防止测量误差的产生。
如图3所示,四个应变片3a、3b、3c、3d分成两组电桥电路输出。位于Y方向,相对的两个应变片3a、3d和两个电阻R1、R2组成一组电桥电路输出,电源一端连接于R1和应变片3a之间,另一端连接于R2和应变片3d之间;电桥电路输出电路U0一端连接于两个应变片3a、3d之间,另一端连接于两个电阻R1和R2之间。U0代表Y方向上,两个应变片3a、3d产生的应变信号输出。
同理,位于X方向,相对的两个应变片3b、3c组成一组电桥电路输出。U0代表X方向上,两个应变片3b、3c产生的应变信号输出。
上述的电阻R1和R2元件可以集成在两个应变片3a、3d中。
如图4和图5所示,另一种轴销式传感器与图1和图2所示的轴销式传感器不同在于,轴销1沿轴线方向开设有通孔,在通孔的内表面上粘贴有四个应变片3a、3b、3c、3d,该四个应变片3a、3b、3c、3d直接粘贴在通孔的内表面,轴销1对称设置,轴销1承受径向力后,不会产生应力集中,防止实测数值与实际径向力产生偏差,提高了检测轴销1所承受径向力的精度,同时结构简单,安装方便,而且可以保护四个应变片3a、3b、3c、3d不受到破坏。
轴销式传感器通过采集四个应变片3a、3b、3c、3d产生的应变信号,经过计算获得轴销所承受径向力的大小及力的作用方位角度,从而达到测量轴销圆周360度任何方向的径向力。
上述的轴销式传感器,应变片数量也可以均匀布置在同一圆周上的3个、5个、6个、7个、8个、9个、10个等。
本发明为了检测轴销圆周360度任何方向的径向力和作用力方位角度,还提供一种检测轴销所承受径向力的方法。
(1)首先将四个应变片3a、3b、3c、3d均匀布置在轴销同一圆周径向上,采用图3的电桥电路连接方式。
(2)采集训练样本:在标定试验机上对轴销进行加载,分别采集与加载力F(方位角为θ)对应的四个应变片的应变信号,采集的应变信号序列用4×n矩阵X表示,Xij(i=1,2,3,4;j=1,…,n)表示在载荷Fj(方位角为θj)作用下第i个应变片的输出,n表示采样序列的长度,X即为标定试验中采集到的训练样本,F,θ为目标输出。
(3)建立神经网络模型或支持向量机模型:运用Matlab神经网络工具箱,建立神经网络模型或支持向量机模型。
(4)学习训练:将(2)中所述训练样本X和目标输出F,θ输入到(3)中所建模型,对模型进行训练。神经网络和SVM采用不同的训练方法,前者最小化经验风险;后者则最小化期望风险。两者的通用做法都是通过将样本集的输入/输出模式反复作用于模型,使模型按照一定的学习算法自动调节自身的参数,当模型的实际输出满足期望要求时,则认为学习圆满结束。使用训练集中的很多样本重复网络的训练,使得在所有样本上的经验风险或期望风险达到最小时,得到理想的轴销式力传感器的模型。具体来说,首先将训练集中每一个样本(Xj,Fj)、(Xj,θj),j=1,…,n,输入模型,计算出模型实际输出的力的大小和力的方位角然后求模型对该模型的预测误差然后计算在所有样本上的累计误差最后根据误差dF,dθ(支持向量机中还考虑期望误差),调整权矩阵W,直到对整个样本集来说,误差不超过规定的范围,即得到理想的轴销式力传感器的预测模型。
(5)实测应变信号:用轴销式力传感器代替原有结构中的轴销,安装在结构连接处,在结构工作过程中,用应变仪分别采集四个应变片的应变信号,记为Y。
(6)预测阶段:将(5)中所述现场实测应变信号Y作为(4)中所述轴销式力传感器预测模型的输入,预测轴销实际受力大小及力的方位角。也就是说,对实际采集的四个应变值的每一个组合Yj,将其输入训练好的轴销式力传感器模型中,模型都会预测出与Yj相对应的径向力Fj以及力的作用方位角θj,对于研究对象在一段连续工作时间内采集到的应变组合向量Y,将其输入训练好的轴销式力传感器模型中,模型就会预测出研究对象在相应工作时间内承受的径向力F以及力的作用方位角θ。
神经网络技术或支持向量机技术在内的智能识别技术是一种对多源不确定性信息进行综合处理的理论和方法,克服了传统的数据处理方法大多针对单一信息源数据,缺乏对多源信息数据的协同利用、综合处理的缺点。基于数据的机器学习是现代智能识别技术中的重要方面,机器学习的目的是根据给定的训练样本求对某系统输入输出之间依赖关系的估计,使它能够对未知输出作出尽可能准确的预测,即研究从观测数据(样本)出发寻找规律,利用这些规律对未来数据或无法观测的数据进行预测。具体来说,本发明所采用的智能识别技术包括神经网络技术和支持向量机技术。所述的神经网络技术和生物神经网络一样,信息处理是通过神经元来实现的,其实质是一个多输入、单输出的非线性系统,其功能为:对每个输入信息加权;对各加权后的信息求和;通过转移函数求输出。所述支持向量机技术是建立在统计学习理论的VC维理论和结构风险最小原理基础上的,根据有限的样本信息在模型的复杂性(即对特定训练样本的学习精度)和学习能力(即无错误地识别任意样本的能力)之间寻求最佳折衷,以期获得最好的推广能力,它在解决小样本、非线性及高维模式识别中表现出许多特有的优势。
上述方法的应变片数量也可以采用位于不同径向上的6个、8个、10个等。应变片在布置上,也可以将Y方向应变片3a、3d与X方向应变片3b、3c错开,相隔一定距离,不布置在轴销的同一圆周面上。
为了检测轴销圆周360度任何方向的径向力和作用力方位角度,还提供另一种检测轴销所承受径向力的方法。该方法考虑了力的正交效应。
(1)首先四个应变片3a、3b、3c、3d采用图3的电桥电路连接方式,四个应变片被分成两组,位置相对的两个应变片为一组,每组应变片分别接入半桥电路,每组应变片产生的应变信号单独输出。相对的两个应变片3b、3c组成一组电桥电路输出代表X的应变值εx+εxy,εx代表径向力在X方向的分力产生应变值,εxy代表径向力在Y方向的分力引起的X方向横向应变值;相对的两个应变片3a、3d组成一组电桥电路输出代表Y的应变值εy+εyx,εy代表径向力在Y方向的分力产生应变值,εyx代表径向力在X方向的分力引起的Y方向横向应变值。
(2)单向载荷作用下,对传感器进行标定,分别求出x、y方向载荷与应变之间的系数k1、k2。也就是说,只在传感器的x方向加载一组力p1[i],得到x方向的一组相应的总应变εx[i](此处取x方向的总应变εx[i]=|εb[i]-εc[i]|,在X方向上应变片3b的应变值εb[i]与应变片3c的应变值εc[i]之差值),用Matlab对εx[i]与p1[i]进行拟合,得到系数k1。同理,对y方向进行相同的处理,得到系数k2。
(3)变方向载荷作用下,对传感器进行标定,求出正交效应系数k12、k21。即当载荷方向变化时,也就是说除X和Y方向径向载荷外,其它方向径向载荷时,由于传感器同时受两个沿正交方向的载荷的作用,实际测出的两个正交分力中都分别混杂有正交方向分力所引起的正交效应力。也就是说,X方向的分力引起的Y方向横向应变值εyx;Y方向的分力引起的X方向横向应变值εxy。
由于测到的应变中已经包含了正交效应的影响,所以实测载荷和真实载荷并不相等。
如下式:
F1=p1+k12p2=k1εx+k12k2εy=k1(εx+εxy),
F2=p2+k21p1=k2εy+k21k1εx=k2(εy+εyx),
根据实测载荷F1、F2和真实载荷p1、p2可求出正交效应系数k12、k21。
式中,F1、F2为实测得出的两个分力;p1、p2为两个互为正交方向的真实分力;k12为p2在其正交方向引起的正交效应系数;k21为p1在其正交方向引起的正交效应系数。也就是说,由于正交效应的影响,两个实测的分力F1、F2与两个真实的分力p1、p2并不相等。
(4)实测应变信号:用全周型轴销式力传感器代替原有结构中的轴销,安装在结构连接处,在结构工作过程中,采集应变信号,记x、y方向的实测应变信号分别为εx+εxy、εy+εyx。
(5)根据下式计算x、y方向的真实分力:
(6)由p1、p2合成总的载荷并计算力的方位角θ=arctan(p2/p1)。
在实际应用当中,根据条件和需要,选择上述任何一种数据处理方法均可方便地实现变方向径向载荷的测量。
上述方法的应变片在布置上,也可以将Y方向应变片3a、3d与X方向应变片3b、3c错开,相隔一定距离,不布置在轴销的同一圆周面上。
本发明克服了传统轴销式传感器不能测量变方向径向载荷的问题,变方向径向载荷是指轴销圆周360度任何方向径向载荷。使得测量工程机械用滑动轴承、滑轮等装置的变方向径向载荷成为可能,扩大了轴销式传感器的应用范围。结构简单、使用方便,并且通过合理地选择传感器的参数,可以达到很高的测量精度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种检测轴销所承受径向力的方法,其特征在于,包括:
在轴销的不同径向上布置有四个应变片,相邻两个应变片径向相隔90度,位置相对的两个应变片为一组,每组应变片分别接入电桥电路,每组应变片产生的应变信号单独输出,分别代表X、Y方向的应变值εx+εxy、εy+εyx;
在X、Y方向对轴销加载径向力,获取x、y方向加载径向力与应变信号之间的系数k1、k2;
在X、Y其它方向对轴销加载径向力,获取正交效应系数k12、k21;
根据下面公式计算x、y方向的真实分力:
由p1、p2合成总的载荷并计算力的方位角θ=arctan(p2/p1)。
2.根据权利要求1所述的方法,其特征在于,
在X、Y其它方向对轴销加载径向力,该径向力在X、Y上的真实分力为P1和P2表示,实测载荷在X、Y方向的分力为F1和F2表示;根据下面公式计算正交效应系数k12、k21;
F1=p1+k12p2=k1εx+k12k2εy=k1(εx+εxy);
F2=p2+k21p1=k2εy+k21k1εx=k2(εy+εyx)。
3.根据权利要求1或2所述的方法,其特征在于,四个应变片均匀布置在轴销的同一圆周径向上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210554820.6A CN103063345B (zh) | 2012-12-19 | 2012-12-19 | 一种检测轴销所承受径向力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210554820.6A CN103063345B (zh) | 2012-12-19 | 2012-12-19 | 一种检测轴销所承受径向力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103063345A CN103063345A (zh) | 2013-04-24 |
CN103063345B true CN103063345B (zh) | 2015-06-10 |
Family
ID=48106064
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210554820.6A Expired - Fee Related CN103063345B (zh) | 2012-12-19 | 2012-12-19 | 一种检测轴销所承受径向力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103063345B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105092121B (zh) * | 2015-08-11 | 2017-11-03 | 中国航空工业集团公司西安飞机设计研究所 | 用于测量刚性管的径向力的方法 |
CN105157890B (zh) * | 2015-09-14 | 2018-04-03 | 重庆长安汽车股份有限公司 | 一种消除摆臂球头受力测量误差的方法 |
CN106438727A (zh) * | 2016-09-08 | 2017-02-22 | 南京工业大学 | 具有超载报警功能的回转支承 |
CN107330218A (zh) * | 2017-07-13 | 2017-11-07 | 徐工集团工程机械有限公司 | 轴销式传感器及其径向力标定与计算方法、装置和系统 |
CN111811714A (zh) * | 2020-07-15 | 2020-10-23 | 中国船舶工业集团公司第七0八研究所 | 一种喷水推进操舵倒航机构销轴径向力测量装置 |
CN112729651B (zh) * | 2021-04-02 | 2021-07-30 | 博鼎精工智能科技(山东)有限公司 | 轴销式力传感器及农机悬挂装置及土壤阻力测量方法 |
CN114370960B (zh) * | 2021-12-29 | 2024-01-26 | 浙江清华柔性电子技术研究院 | 拉杆载荷测量方法、装置、系统及存储介质 |
CN114577379B (zh) * | 2022-02-28 | 2023-05-12 | 太原理工大学 | 销轴式力传感器径向力大小及方向的测量方法 |
CN115371858A (zh) * | 2022-08-03 | 2022-11-22 | 中联重科股份有限公司 | 销轴传感组件的信号处理方法、系统、装置及测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2450647Y (zh) * | 2000-11-09 | 2001-09-26 | 张家港保税区中诚国际工贸有限公司 | 动态轴销式测力传感器 |
CN101103259A (zh) * | 2004-12-24 | 2008-01-09 | 空中客车英国有限公司 | 用于测量轴承销所承受的负载的设备和方法 |
CN102269618A (zh) * | 2010-05-21 | 2011-12-07 | 杰富意先进技术株式会社 | 称重装置 |
CN202177476U (zh) * | 2011-07-19 | 2012-03-28 | 宁波金旭电子有限公司 | 轴销式传感器 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003334729A (ja) * | 2002-05-14 | 2003-11-25 | Mitsubishi Heavy Ind Ltd | 工作機械の主軸装置、該装置に用いる荷重検出器とその測定方法及び予圧調整方法 |
-
2012
- 2012-12-19 CN CN201210554820.6A patent/CN103063345B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2450647Y (zh) * | 2000-11-09 | 2001-09-26 | 张家港保税区中诚国际工贸有限公司 | 动态轴销式测力传感器 |
CN101103259A (zh) * | 2004-12-24 | 2008-01-09 | 空中客车英国有限公司 | 用于测量轴承销所承受的负载的设备和方法 |
CN102269618A (zh) * | 2010-05-21 | 2011-12-07 | 杰富意先进技术株式会社 | 称重装置 |
CN202177476U (zh) * | 2011-07-19 | 2012-03-28 | 宁波金旭电子有限公司 | 轴销式传感器 |
Also Published As
Publication number | Publication date |
---|---|
CN103063345A (zh) | 2013-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103063345B (zh) | 一种检测轴销所承受径向力的方法 | |
CN103063341B (zh) | 一种检测轴销所承受径向力的方法 | |
CN103245928B (zh) | 方向可调的均匀磁场和均匀一阶梯度磁场的方法及装置 | |
CN101907448B (zh) | 一种基于双目三维视觉的深度测量方法 | |
CN111900731B (zh) | 一种基于pmu的电力系统状态估计性能评价方法 | |
CN105203055A (zh) | 一种关节式坐标测量机动态误差补偿方法 | |
CN110285781B (zh) | 一种相对于基准面的平面平行度快速评定方法 | |
CN105841857B (zh) | 一种并联式五维力传感器 | |
CN104677257A (zh) | 基于球面电容原理的精密球铰链关节间隙检测方法 | |
Zhang et al. | Health factor extraction of lithium-ion batteries based on discrete wavelet transform and soh prediction based on catboost | |
Antonucci et al. | Li-ion battery modeling and state of charge estimation method including the hysteresis effect | |
CN110263826A (zh) | 含噪声非线性过程故障检测模型的构建方法及其检测方法 | |
CN103940453B (zh) | 一种提高传感器测量精度方法 | |
CN205719350U (zh) | 一种并联式五维力传感器 | |
CN107587955B (zh) | 基于深信度网络的火箭发动机推力偏移量的标定方法 | |
CN109101759A (zh) | 一种基于正逆响应面法的参数识别方法 | |
CN113973403A (zh) | 基于结构离散测点拓扑的温致应变场重分布智能感知方法 | |
Wang et al. | Decoupling research of a three-dimensional force tactile sensor based on radical basis function neural network | |
CN107478170B (zh) | 一种光纤应变花和智能涂层集成传感器的实现方法 | |
CN206056844U (zh) | 一种轴动态承载力测量试验装置 | |
CN105806943B (zh) | 基于离散小波的复合材料梁结构损伤定位方法及系统 | |
Lu et al. | Application of Mechanical Electronic Engineering Technology in Sensor Measurement System | |
Zeng et al. | An Intelligent Correlation Real-Time Analysis Method for the Mechanical Properties of Members in Super-Span Valve Hall Grid Structure Hoisting Process | |
Wang et al. | Application of stir tool force measuring dynamometer for friction stir welding of aluminum alloys | |
Senthilkumar et al. | Non destructive fabric weight measurement using capacitance principle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150610 Termination date: 20181219 |