CN108802427B - 预紧式并联六维加速度传感器及其测量和灵敏度分析方法 - Google Patents

预紧式并联六维加速度传感器及其测量和灵敏度分析方法 Download PDF

Info

Publication number
CN108802427B
CN108802427B CN201810599723.6A CN201810599723A CN108802427B CN 108802427 B CN108802427 B CN 108802427B CN 201810599723 A CN201810599723 A CN 201810599723A CN 108802427 B CN108802427 B CN 108802427B
Authority
CN
China
Prior art keywords
tightening
sensor
branched chain
sensitivity
acceleration
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
Application number
CN201810599723.6A
Other languages
English (en)
Other versions
CN108802427A (zh
Inventor
张文妍
王艳
李成刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Forestry University
Original Assignee
Nanjing Forestry University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Forestry University filed Critical Nanjing Forestry University
Priority to CN201810599723.6A priority Critical patent/CN108802427B/zh
Publication of CN108802427A publication Critical patent/CN108802427A/zh
Application granted granted Critical
Publication of CN108802427B publication Critical patent/CN108802427B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P15/00Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
    • G01P15/18Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration in two or more dimensions
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P15/00Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
    • G01P15/02Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration by making use of inertia forces using solid seismic masses
    • G01P15/08Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration by making use of inertia forces using solid seismic masses with conversion into electric or magnetic values
    • G01P15/09Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration by making use of inertia forces using solid seismic masses with conversion into electric or magnetic values by piezoelectric pick-up
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P21/00Testing or calibrating of apparatus or devices covered by the preceding groups

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

本发明公开了一种预紧式并联六维加速度传感器及其测量方法和灵敏度分析方法。该传感器采用八支链并联机构作为传感器的弹性体结构,具有更高的刚度和更大的承载能力;可通过预紧装置对压电陶瓷施加预紧力,在预紧之后,该传感器能够有效避免压电陶瓷输出信号失真的问题,使得传感器具有更好的测量精度。该测量方法只需预先在计算机中编写程序,就能够实现六维加速度的实时测量,而其灵敏度分析方法通过挖掘传感器并联机构中固有约束关系,建立力协调方程,将传感器的灵敏度模型超静定问题求解转化为静定问题求解,从而推导得到传感器的灵敏度的解析表达式,为后续传感器的结构优化设计提供理论指导。

Description

预紧式并联六维加速度传感器及其测量和灵敏度分析方法
技术领域
本发明属于惯性测量技术领域,具体涉及一种预紧式并联六维加速度传感器及其测量方法和灵敏度分析方法。
背景技术
加速度传感器作为获取信息的关键部件,在工程测振、工业控制、惯性导航与制导等领域发挥着重要作用。传统的低维加速度传感器的理论与技术已日臻完善,市场上的各种加速度传感器产品也屡见不鲜,但多以单轴的传感器为主。随着信息的多元化发展,信息测量要求的提高,人们不再仅仅局限于纯线加速度和纯角加速度参量的测量,对加速度传感器提出了更高的要求。
对于一个在三维空间运动的物体,若想获得信噪比良好的数据,还原物体的真实运动状态,我们需要知道其沿三个坐标轴的线加速度和绕三个坐标轴的角加速度。例如,机械臂末端执行器在捕获、抓取物体的过程中会产生较大的振动,另外,机械臂的柔性关节处也会产生残余振动,会对机器人控制性能造成影响,而振动具有不确定性,因此需要获取机械臂末端的六维加速度信息以实现机械臂惯性力的补偿和残余振动的抑制。现有惯性导航系统中的测量装置都是采用线加速度计和陀螺仪以获得载体相对于惯性坐标系的实时速度、姿态和航向等信息。目前线加速度计的技术已经非常成熟,但传统的机械式陀螺仪会存在一定的交叉耦合,对外界环境较为敏感,且不能承受较大的振动冲击,测量精度较低,光纤陀螺、环形激光陀螺等光学陀螺因其制造成本较高、技术难度大、不易微型化等特点还未能广泛应用。因此,研制一种不采用陀螺仪也能测量载体完整的六维加速度信息成为惯性领域重点关注的问题。
目前国内外对于六维加速度传感器的研究还不够系统全面,仍处于理论探索阶段。传统的六维加速度测量方法一般是利用多个多轴加速度计组合作为惯性元件来代替陀螺测量载体的角速度参量,但多个加速度计的组合安装方式使其质心位置并不一致,使得系统对加速度计的加工工艺及安装精度要求太高,结构难以微型化,且在实际应用中还存在可观的交叉耦合,角速度解算误差积累较快,导致其传感精度较低。
并联机构以结构紧凑且稳定、承载能力高、无误差累积等固有特点,成为一类六维力(加速度)传感器的结构实现形式。申请号201010249900.1,名称为“冗余预紧式并联六维加速度传感器及其测量方法”的中国专利公开了一种基于9—SPS并联机构构型的压电式六维加速度传感器,其充分结合了并联机构和压电陶瓷的优点,解耦效率高,并且克服了基于多个单轴加速度传感器组合测量系统的缺点。但其不足之处在于并联机构中引入了三重复合铰,结构复杂、加工难度大,实际工程中很难保证三个柔性球铰的转动中心完全重合于一点;且在工作过程中,压电陶瓷不能承受较大的拉力,在使用前需要施进行预紧,但由于结构的限制,该构型无法沿支链的轴向施加预紧载荷。
尽管六维加速度传感器的研究尚处于理论探索阶段,但其应用趋势却日益加强,应用前景也日益广阔,因此,有必要设计一种新型可预紧的并联机构构型来充当六维加速度传感器的弹性体结构。同时,为了促进六维加速度传感器产品化,有必要针对传感器的性能指标进行探索,为后续传感器结构优化设计提供依据。
发明内容
发明目的:为了克服现有技术的缺陷,本发明提供了一种预紧式并联六维加速度传感器该六维加速度传感器具有更高的刚度和更大的承载能力,同时也具有更好的测量精度。
同时,本发明还提供了上述预紧式并联六维加速度传感器的测量方法和灵敏度分析方法。
技术方案:本发明所述的预紧式并联六维加速度传感器,包括外壳、位于外壳内的立方体质量块、八条预紧支链;所述预紧支链包括第一球铰连接杆、第二球铰连接杆、压电陶瓷和用以对压电陶瓷施加预紧力的预紧装置,所述第一球铰连接杆和第二球铰连接杆分别连接在所述压电陶瓷的两端;预紧支链的一端通过第一球铰连接杆连接于所述立方体质量块上,预紧支链的另一端通过第二球铰连接杆连接于所述外壳的壁面上;
其中,所述立方体质量块的中心为质心,立方体质量块具有六个表面,设相邻且相互垂直的三个表面分别为第一表面、第二表面和第三表面,平行于所述第三表面的为第四表面;所述八条预紧支链分别为第一预紧支链、第二预紧支链、第三预紧支链、第四预紧支链、与第三预紧支链关于所述质心对称的第五预紧支链、与第二预紧支链关于质心对称的第六预紧支链、与第一预紧支链关于质心对称的第七预紧支链、与第四预紧支链关于质心对称的第八预紧支链;
所述第一预紧支链的延伸方向垂直于所述第一表面,第一预紧支链的一端连接于第一表面、第二表面、第四表面相交的顶点;所述第二预紧支链的延伸方向垂直于所述第二表面,第二预紧支链的一端连接于第一表面、第二表面、第三表面相交的顶点;所述第三预紧支链的延伸方向垂直于所述第三表面,第三预紧支链的一端连接于第二表面和第三表面相交棱边的中点;所述第四预紧支链的延伸方向垂直于所述第一表面,第四预紧支链的一端连接于第一表面和第三表面相交棱边的中点。
有益效果:本发明通过设置八条预紧支链并联机构作为传感器的弹性体拓扑结构,相较于传统的Stewart并联机构增加了两条冗余支链,一方面,该传感器具有更高的刚度和更大的承载能力;另一方面,可以通过预紧装置对压电陶瓷施加预紧力,在预紧之后,该传感器能够有效避免压电陶瓷输出信号失真的问题,使得传感器具有更好的测量精度。同时,该传感器结构简单,加工难度低,有利于六维加速度传感器的微型化。
具体的,所述预紧装置包括具有内螺纹的套筒、具有外螺纹的两个锁紧螺母、两个弹性件;所述两个锁紧螺母一一对应分别套设在所述第一球铰连接杆和所述第二球铰连接杆上且分别螺纹连接于套筒的内部两端;第一球铰连接杆和第二球铰连接杆在与压电陶瓷连接处均具有台肩;所述两个弹性件的两端分别抵靠在对应的锁紧螺母和台肩上。
进一步的,所述外壳的各个壁面一一对应平行于所述立方体质量块的各个表面,外壳在至少一个壁面上设置有连通外壳内外的预紧装置调节口,从而便于对各个预紧支链进行预紧调节。
本发明所述的使用预紧式并联六维加速度传感器的测量方法采用的技术方案是,包括如下步骤:
(a)对八条预紧支链的压电陶瓷施加初始预紧力,压电陶瓷始终保持受压状态;
(b)将预紧式并联六维加速度传感器刚性固定于待测载体上;
(c)将预紧式并联六维加速度传感器上八个压电陶瓷的两极通过电荷输出导线分别连接电荷放大器八通道的输入端,借助电荷放大器、数据采集卡实现电荷信号的调理和采集,将其转换为数字信号供计算机分析处理,最终得到待测载体的实时加速度。
所述步骤(b)中,首先分别在外壳和地面上固结坐标系{W}和{I},在运动初始时刻,两个坐标系重合且坐标系原点位于立方体质量块的质心,三个坐标轴方向与立方体质量块的三条正交棱边相平行;
待测载体实时加速度的计算步骤如下:
(1)通过公式Qi=Ci·Vi,i=1~8,将电荷放大器八个通道输出的电压值转换为输入的电荷量,式中,Qi为第i通道的实际电荷量,也即第i组压电陶瓷两极产生的电荷量;Ci为电荷转换系数,其表达式为Ci=Sq/Su,其中,Sq为电荷放大器的输入电荷灵敏度,Su为电荷放大器的输出电压灵敏度,Sq、Su均可调;Vi为第i通道的输出电压,含正负号;
(2)分别通过公式ΔLi=Qi/(Ki·d33)和fi=Qi/d33,将八个压电陶瓷两极产生的电荷量转换为各预紧支链的压电陶瓷的动态变形量和各预紧支链的压电陶瓷所受到的压力;式中,ΔLi为压电陶瓷的动态变形量,fi为各预紧支链的压电陶瓷所受到的压力,即预紧支链所受的轴向力,ΔLi和fi的正负号由Qi决定;Ki为第i个压电陶瓷的等效刚度,d33为第i组压电陶瓷的等效压电系数;
(3)假设立方体质量块的质量和半边长分别为m和n,根据牛顿-欧拉法构建传感器系统的动力学方程组,
其中,立方体质量块的质心的绝对线加速度aI表示为
Figure GDA0002383049790000041
质量块的绝对角加速度εI表示为
Figure GDA0002383049790000042
式中,
Figure GDA0002383049790000043
为{W}系相对于{I}系的旋转矩阵,通过引入辅助角速度ω*的方式求解得出;
(4)将八条预紧支链产生的轴向力fi
Figure GDA0002383049790000051
代入系统的动力学方程组,即求得立方体质量块的绝对线加速度aI和绝对角加速度εI,也即待测载体的实时加速度。
有益效果:通过将传感器的外壳刚性固定于待测载体上并随载体一起运动,那么传感器的外壳的绝对加速度即为待测载体的六维加速度。传感器外壳的绝对加速度可由立方体质量块相对于惯性系的绝对运动参量和立方体质量块相对于外壳的运动参量推算得到,但由于压电陶瓷的刚度较大,质量块与外壳之间的相对运动可忽略不计,因此,立方体质量块的质心的绝对线加速度和绝对角加速度即为待测载体的六维加速度。通过上述测量方法,预先在计算机中编写程序,利用预紧式并联六维加速度传感器就能够实现六维加速度的实时测量,且在工作过程中,能够有效避免压电陶瓷在承受一定拉力时所造成的输出信号失真现象。
对应于上述预紧式并联六维加速度传感器,本发明提供的该传感器的灵敏度分析方法可采用如下技术方案:
灵敏度包括线加速度灵敏度和角加速度灵敏度,依次推导传感器沿x、y和z三个坐标轴方向的线加速度灵敏度和角加速度灵敏度;将沿三个坐标轴方向的线加速度灵敏度的最低值作为传感器的线加速度灵敏度,将沿三个坐标轴方向的角加速度灵敏度的最低值作为传感器的角加速度灵敏度。
定义所述传感器的灵敏度为:加速度分量每增加一个单位值时对应的八条预紧支链所受轴向力的变化量的绝对值之和。
构建传感器系统的六维加速度参量与八条预紧支链的轴向力之间的映射关系,得到预紧式并联六维加速度传感器所有预紧支链在纯线加速度或者纯角加速度驱动情况下的轴向力情况;推导预紧式并联六维加速度传感器灵敏度数学模型,包括如下步骤:
S1:选取外壳或与外壳刚性连接的载体为研究对象,对预紧式并联六维加速度传感器模型进行简化,构建传感器在纯线加速度运动或者纯角加速度运动情况下的动力学方程组;根据牛顿第二定律列写传感器沿x轴方向的纯线加速度ax和ax+Δax驱动下立方体质量块的两组动力学方程:
Figure GDA0002383049790000061
Figure GDA0002383049790000062
式中,ax和ax+Δax为沿x轴方向输入传感器的已知加速度参量;f1~f8为八条预紧支链的轴向力;m为立方体质量块的质量,n为立方体质量块的半边长;
S2:基于并联机构的正向运动学理论挖掘预紧式并联六维加速度传感器各个支链的变形协调条件,推导出该传感器的两个力协调方程:
Figure GDA0002383049790000063
Figure GDA0002383049790000064
式中,fi
Figure GDA0002383049790000065
i=1~8分别为第i条预紧支链所受的轴向力和初始预紧力;ki,i=1~8为第i条预紧支链的轴向刚度,由于八条预紧支链结构完全相同,则ki=k,则两个协调方程分别整理为:
Figure GDA0002383049790000066
Figure GDA0002383049790000067
S3:结合式(1)、式(5)和式(6),得纯线加速度ax驱动对应的预紧支链的轴向力f,写成矩阵的形式:
f=C-1D1+C-1D2 (7)
式中,
Figure GDA0002383049790000071
f=[f1 f2 f3 f4 f5 f6 f7 f8]T
D1=[max 0 0 0 0 0 0 0]T
D2=[0 0 0 0 0 0 b1 b2]T
Figure GDA0002383049790000072
Figure GDA0002383049790000073
S4:结合式(2)、式(5)和式(6),得纯线加速度ax+Δax驱动对应的预紧支链的轴向力f*,写成矩阵的形式:
Figure GDA0002383049790000074
式中,
Figure GDA0002383049790000075
S5:将式(7)和式(8)的等号两边分别相减,得到:
Figure GDA0002383049790000076
S6:根据所述传感器的灵敏度的定义,得到所述传感器沿x轴方向的线加速度灵敏度Sax的定义式
Figure GDA0002383049790000077
并结合式(9),推得传感器沿x轴方向的线加速度灵敏度
Figure GDA0002383049790000078
S7:按照步骤S1~S6推导得到所述传感器沿y轴方向的线加速度灵敏度
Figure GDA0002383049790000079
传感器沿z轴方向的线加速度灵敏度
Figure GDA0002383049790000081
传感器沿x轴方向的角加速度灵敏度
Figure GDA0002383049790000082
传感器沿y轴方向的角加速度灵敏度
Figure GDA0002383049790000083
和传感器沿z轴方向的角加速度灵敏度
Figure GDA0002383049790000084
S8:最终推导出该传感器的线加速度灵敏度为Sa=min(Sax,Say,Saz)=m,角加速度灵敏度为Sε=min(Sεx,Sεy,Sεz)=mn。
有益效果:通过上述灵敏度分析方法,实现对预紧式并联六维加速度传感器的灵敏度的推导,挖掘预紧式并联六维加速度传感器各个预紧支链的变形协调条件,得到传感器的力协调条件,并与传感器的动力学方程结合从而推导得到传感器的灵敏度数学模型,有利于后续对传感器进行结构优化,同时有利于将预紧式并联六维加速度传感器产品化。
附图说明
图1是本发明预紧式并联六维加速度传感器的结构示意图;
图2是本发明的预紧支链的剖视结构图;
图3是本发明的三维模型示意图;
图4是本发明的灵敏度数学模型推导过程示意图。
具体实施方式
以下将结合附图及具体实施例对本发明的结构及工作过程进行详细说明。
首先,参考图1所示,本发明提供一种预紧式并联六维加速度传感器,包括一个外壳1、一个立方体质量块2和八条预紧支链3,下面分别介绍。
外壳1为刚性结构,本实施例中设计为T型结构,具有呈板状且带支脚的水平部分和大致呈立方体结构的竖直部分,并在其竖直部分开设有大致呈立方体结构的容置腔。
立方体质量块2设置于外壳1的容置腔内,外壳1的各个壁面一一对应平行于所述立方体质量块2的各个表面。
请一并结合图2,预紧支链3包括第一球铰连接杆4、第二球铰连接杆5、压电陶瓷6和预紧装置。每条预紧支链3上均具有一组压电陶瓷6,用于敏感各支链的应变与应力,压电陶瓷6为若干压电陶瓷片堆叠而成,压电陶瓷的两端分别设有正电极和负电极,自所述正电极和负电极分别延伸出电荷传输导线并焊接在两侧。所述第一球铰连接杆5和第二球铰连接杆6分别连接在所述压电陶瓷6的两端,所述预紧装置连接所述第一球铰连接杆4和第二球铰连接杆5用以对压电陶瓷6施加预紧力,预紧支链3的一端通过第一球铰连接杆4连接于所述立方体质量块上,预紧支链的另一端通过第二球铰连接杆5连接于所述外壳1的壁面上。
具体的,预紧装置包括具有内螺纹的套筒71、具有外螺纹的两个锁紧螺母72、两个弹性件73;所述两个锁紧螺母72一一对应分别套设在所述第一球铰连接杆4和所述第二球铰连接杆5上,且分别螺纹连接于套筒71的内部两端;同时,第一球铰连接杆4和第二球铰连接杆5均通过环氧树脂胶与压电陶瓷6粘接,第一球铰连接杆4和第二球铰连接杆5在与压电陶瓷6的连接处均具有台肩8;所述两个弹性件73一一对应分别围绕于第一球铰连接杆4和第二球铰连接杆5上,弹性件73的两端分别抵靠在对应的锁紧螺母72和台肩8上。通过旋动锁紧螺母72,可以改变弹性件73的形变量,进而对压电陶瓷6施加预紧力。
进一步的,为了便于对预紧装置进行调节,外壳1在至少一个壁面上设置有连通外壳1内外的预紧装置调节口11。本实施例中,顶部的壁面为可拆卸式连接在左右两侧壁面之间,顶部壁面的前后两侧分别预留有预紧装置调节口11。
请参阅图3,其中A1~A8分别示意八条预紧支链3的第一球铰连接杆4与立方体质量块2的连接节点,B1~B8分别示意八条预紧支链3的第二球铰连接杆5与立方体质量块2的连接节点,P1~P8分别示意八条预紧支链3的压电陶瓷6对应的信号。为了便于描述,这里以P1~P8分别作为八条预紧支链3的标识。
所述立方体质量块2的中心为质心,立方体质量块2具有六个表面,为了便于描述,设相邻且相互垂直的三个表面分别为第一表面、第二表面和第三表面,平行于第三表面的为第四表面;本实施例中第一表面为图1和图3中示意的立方体质量块2的下表面,观察图时直接面对的在竖直方向上的两个表面中,右侧为第二表面,左侧为第四表面。
所述八条预紧支链3分别为第一预紧支链P1、第二预紧支链P2、第三预紧支链P3、第四预紧支链P4、与第三预紧支链P3关于所述质心对称的第五预紧支链P5、与第二预紧支链P2关于质心对称的第六预紧支链P6、与第一预紧支链P1关于质心对称的第七预紧支链P7、与第四预紧支链P4关于质心对称的第八预紧支链P8
所述第一预紧支链P1的延伸方向垂直于所述第一表面,第一预紧支链P1的一端连接于第一表面、第二表面、第四表面相交的顶点;所述第二预紧支链P2的延伸方向垂直于所述第二表面,第二预紧支链P2的一端连接于第一表面、第二表面、第三表面相交的顶点;所述第三预紧支链P3的延伸方向垂直于所述第三表面,第三预紧支链P3的一端连接于第二表面和第三表面相交棱边的中点;所述第四预紧支链P4的延伸方向垂直于所述第一表面,第四预紧支链P4的一端连接于第一表面和第二表面相交棱边的中点。
实际测量时,将该传感器的外壳1刚性固定于待测载体上并随载体一起运动,外壳1的绝对加速度即为待测载体的六维加速度。而传感器的外壳1的绝对加速度可由立方体质量块2相对于惯性系的绝对运动参量和立方体质量块2相对于外壳1的运动参量推算得到,但由于压电陶瓷6的刚度较大,立方体质量块2与外壳1之间的相对运动可忽略不计,因此,立方体质量块2质心的绝对线加速度和绝对角加速度即为待测载体的六维加速度。
本发明也提供了一种预紧式并联六维加速度传感器的测量方法,该方法的步骤如下:
在工作之前,对八条预紧支链的压电陶瓷施加初始预紧力,通过旋动锁紧螺母调节碟簧的压缩量以改变预紧力的大小,保证在工作过程中压电陶瓷始终保持受压状态。
然后将预紧式并联六维加速度传感器刚性固定于待测载体上。并将预紧式并联六维加速度传感器上八个压电陶瓷的两极通过电荷输出导线分别连接电荷放大器八通道的输入端,借助电荷放大器、数据采集卡实现电荷信号的调理和采集,将其转换为数字信号供计算机分析处理,最终得到待测载体的实时加速度。
具体的,首先分别在外壳和地面上固结坐标系{W}和{I},在运动初始时刻,两个坐标系重合且坐标系原点位于立方体质量块的质心,三个坐标轴方向与立方体质量块的三条正交棱边相平行;利用传感器上八组压电陶瓷的输出电荷量计算待测载体实时加速度,待测载体实时加速度的计算步骤如下:
(1)通过公式Qi=Ci·Vi,i=1~8,将电荷放大器八个通道输出的电压值转换为输入的电荷量。式中,Qi为第i通道的实际电荷量,也即第i组压电陶瓷两极产生的电荷量;Ci为电荷转换系数,其表达式为Ci=Sq/Su,其中,Sq为电荷放大器的输入电荷灵敏度,Su为电荷放大器的输出电压灵敏度,Sq、Su均可调;Vi为第i通道的输出电压,含正负号;
(2)分别通过公式ΔLi=Qi/(Ki·d33)和fi=Qi/d33,将八个压电陶瓷两极产生的电荷量转换为各预紧支链的压电陶瓷的动态变形量和各预紧支链的压电陶瓷所受到的压力;式中,ΔLi为压电陶瓷的动态变形量,fi为各预紧支链的压电陶瓷所受到的压力,即预紧支链所受的轴向力,ΔLi和fi的正负号由Qi决定;Ki为第i个压电陶瓷的等效刚度,d33为第i组压电陶瓷的等效压电系数,与压电陶瓷的型号、作用力方向等因素有关;
(3)假设立方体质量块的质量和半边长分别为m和n,根据牛顿-欧拉法构建传感器系统的动力学方程组,
其中,立方体质量块的质心的绝对线加速度aI可表示为
Figure GDA0002383049790000111
质量块的绝对角加速度εI可表示为
Figure GDA0002383049790000112
式中,
Figure GDA0002383049790000113
为{W}系相对于{I}系的旋转矩阵,可通过引入辅助角速度ω*的方式求解得出;参考文献《基于并联机构的的六维加速度传感器的方案设计及建模研究》中通过引入辅助角速度ω*的方式,进而求解出
Figure GDA0002383049790000114
(4)将八条预紧支链产生的轴向力fi
Figure GDA0002383049790000115
代入系统的动力学方程组,即可求得立方体质量块的绝对线加速度aI和绝对角加速度εI,也即待测载体的实时加速度。
可以看出:只要按照以上步骤预先在计算机中编写程序,利用本发明所提供的预紧式并联六维加速度传感器就能够实现六维加速度的实时测量,且在工作过程中,能够有效避免压电陶瓷在承受一定拉力时所造成的输出信号失真现象。
如图4所示,本发明还提供了该预紧式并联六维加速度传感器的灵敏度分析方法,在对上述传感器进行灵敏度建模时,首先构建传感器系统的两组动力学方程,然后通过挖掘预紧式并联六维加速度传感器各个支链的变形协调条件,推导出该传感器的两个力协调方程,将传感器的灵敏度模型超静定问题求解转化为静定问题求解,从而得到传感器加速度灵敏度的解析表达式。
具体的,灵敏度包括线加速度灵敏度和角加速度灵敏度,依次推导传感器沿x、y和z三个坐标轴方向的线加速度灵敏度和角加速度灵敏度;将沿三个坐标轴方向的线加速度灵敏度的最低值作为传感器的线加速度灵敏度,将沿三个坐标轴方向的角加速度灵敏度的最低值作为传感器的角加速度灵敏度。
参考文献《传感器原理及工程应用》中对传感器灵敏度的定义:输出量增量Δy与引起输出量增量Δy的响应输入量增量Δx之比,将所述预紧式并联六维加速度传感器的灵敏度定义为:加速度分量每增加一个单位值时对应的八条预紧支链所受轴向力的变化量的绝对值之和。
进一步,构建传感器系统的六维加速度参量与八条预紧支链的轴向力之间的映射关系,得到传感器所有预紧支链在纯线加速度或者纯角加速度驱动情况下的轴向力情况;推导传感器灵敏度数学模型,包括如下步骤:
S1:选取外壳或与外壳刚性连接的载体为研究对象,对预紧式并联六维加速度传感器模型进行简化,构建传感器在纯线加速度运动或者纯角加速度运动情况下的动力学方程组;根据牛顿第二定律列写传感器沿x轴方向的纯线加速度ax和ax+Δax驱动下立方体质量块的两组动力学方程:
Figure GDA0002383049790000121
Figure GDA0002383049790000131
式中,ax和ax+Δax为沿x轴方向输入传感器的已知加速度参量;f1~f8为八条预紧支链的轴向力;m为立方体质量块的质量,n为立方体质量块的半边长;
S2:基于并联机构的正向运动学理论挖掘预紧式并联六维加速度传感器各个支链的变形协调条件,推导出该传感器的两个力协调方程:
Figure GDA0002383049790000132
Figure GDA0002383049790000133
式中,fi
Figure GDA0002383049790000134
i=1~8分别为第i条预紧支链所受的轴向力和初始预紧力;ki,i=1~8为第i条预紧支链的轴向刚度,由于八条预紧支链结构完全相同,则ki=k,则两个协调方程可分别整理为:
Figure GDA0002383049790000135
Figure GDA0002383049790000136
S3:结合式(1)、式(5)和式(6),可得纯线加速度ax驱动对应的预紧支链的轴向力f,写成矩阵的形式:
f=C-1D1+C-1D2 (7)
式中,
Figure GDA0002383049790000141
f=[f1 f2 f3 f4 f5 f6 f7 f8]T
D1=[max 0 0 0 0 0 0 0]T
D2=[0 0 0 0 0 0 b1 b2]T
Figure GDA0002383049790000142
Figure GDA0002383049790000143
S4:结合式(2)、式(5)和式(6),可得纯线加速度ax+Δax驱动对应的预紧支链的轴向力f*,写成矩阵的形式:
Figure GDA0002383049790000144
式中,
Figure GDA0002383049790000145
S5:将式(7)和式(8)的等号两边分别相减,得到:
Figure GDA0002383049790000146
S6:根据所述传感器的灵敏度的定义,得到所述传感器沿x轴方向的线加速度灵敏度Sax的定义式
Figure GDA0002383049790000147
并结合式(9),推得传感器沿x轴方向的线加速度灵敏度
Figure GDA0002383049790000148
S7:按照步骤S1~S6推导得到所述传感器沿y轴方向的线加速度灵敏度
Figure GDA0002383049790000149
传感器沿z轴方向的线加速度灵敏度
Figure GDA0002383049790000151
传感器沿x轴方向的角加速度灵敏度
Figure GDA0002383049790000152
传感器沿y轴方向的角加速度灵敏度
Figure GDA0002383049790000153
和传感器沿z轴方向的角加速度灵敏度
Figure GDA0002383049790000154
S8:最终推导出该传感器的线加速度灵敏度为Sa=min(Sax,Say,Saz)=m,角加速度灵敏度为Sε=min(Sεx,Sεy,Sεz)=mn。
至此,传感器的线加速度灵敏度和角加速度灵敏度的数学模型已推导完毕。
可以看出,在前述的推导过程中,通过挖掘支链的变形协调条件建立补充方程,将传感器的灵敏度数学模型求解问题转化为静定问题求解,得到传感器的六维加速度参量与支链轴向力之间的简单映射关系,从而推导出传感器的加速度灵敏度的解析表达式。
以上所述仅是本发明的优选实施方式,应当指出在不脱离本发明的构思的前提下,还可以做出若干推演或替代,这些推演或替代都应视为本发明的保护范围。

Claims (9)

1.一种预紧式并联六维加速度传感器,其特征在于,包括外壳(1)、位于外壳内的立方体质量块(2)、八条预紧支链(3);所述预紧支链包括第一球铰连接杆(4)、第二球铰连接杆(5)、压电陶瓷(6)和用以对压电陶瓷施加预紧力的预紧装置,所述第一球铰连接杆和第二球铰连接杆分别连接在所述压电陶瓷的两端;预紧支链的一端通过第一球铰连接杆连接于所述立方体质量块上,预紧支链的另一端通过第二球铰连接杆连接于所述外壳的壁面上;
其中,所述立方体质量块的中心为质心,立方体质量块具有六个表面,设相邻且相互垂直的三个表面分别为第一表面、第二表面和第三表面,平行于所述第三表面的为第四表面;所述八条预紧支链分别为第一预紧支链、第二预紧支链、第三预紧支链、第四预紧支链、与第三预紧支链关于所述质心对称的第五预紧支链、与第二预紧支链关于质心对称的第六预紧支链、与第一预紧支链关于质心对称的第七预紧支链、与第四预紧支链关于质心对称的第八预紧支链;
所述第一预紧支链的延伸方向垂直于所述第一表面,第一预紧支链的一端连接于第一表面、第二表面、第四表面相交的顶点;所述第二预紧支链的延伸方向垂直于所述第二表面,第二预紧支链的一端连接于第一表面、第二表面、第三表面相交的顶点;所述第三预紧支链的延伸方向垂直于所述第三表面,第三预紧支链的一端连接于第二表面和第三表面相交棱边的中点;所述第四预紧支链的延伸方向垂直于所述第一表面,第四预紧支链的一端连接于第一表面和第三表面相交棱边的中点。
2.根据权利要求1所述的预紧式并联六维加速度传感器,其特征在于,所述预紧装置包括具有内螺纹的套筒(71)、具有外螺纹的两个锁紧螺母(72)、两个弹性件(73);所述两个锁紧螺母一一对应分别套设在所述第一球铰连接杆和所述第二球铰连接杆上且分别螺纹连接于套筒的内部两端;第一球铰连接杆和第二球铰连接杆在与压电陶瓷连接处均具有台肩(8);所述两个弹性件的两端分别抵靠在对应的锁紧螺母和台肩上。
3.根据权利要求1所述的预紧式并联六维加速度传感器,其特征在于,所述外壳的各个壁面一一对应平行于所述立方体质量块的各个表面,外壳在至少一个壁面上设置有连通外壳内外的预紧装置调节口(11)。
4.根据权利要求1所述的预紧式并联六维加速度传感器,其特征在于,所述压电陶瓷为若干压电陶瓷片堆叠而成,压电陶瓷的两端分别设有正电极和负电极,自所述正电极和负电极分别延伸出电荷传输导线。
5.一种使用如权利要求1~4任一所述预紧式并联六维加速度传感器的测量方法,其特征在于,包括如下步骤:
(a)对八条预紧支链的压电陶瓷施加初始预紧力,压电陶瓷始终保持受压状态;
(b)将预紧式并联六维加速度传感器刚性固定于待测载体上;
(c)将预紧式并联六维加速度传感器上八个压电陶瓷的两极通过电荷输出导线分别连接电荷放大器八通道的输入端,借助电荷放大器、数据采集卡实现电荷信号的调理和采集,将其转换为数字信号供计算机分析处理,最终得到待测载体的实时加速度。
6.根据权利要求5所述的测量方法,其特征在于,首先分别在外壳和地面上固结坐标系{W}和{I},在运动初始时刻,两个坐标系重合且坐标系原点位于立方体质量块的质心,三个坐标轴方向与立方体质量块的三条正交棱边相平行;
待测载体实时加速度的计算步骤如下:
(1)通过公式Qi=Ci·Vi,i=1~8,将电荷放大器八个通道输出的电压值转换为输入的电荷量,式中,Qi为第i通道的实际电荷量,也即第i组压电陶瓷两极产生的电荷量;Ci为电荷转换系数,其表达式为Ci=Sq/Su,其中,Sq为电荷放大器的输入电荷灵敏度,Su为电荷放大器的输出电压灵敏度,Sq、Su均可调;Vi为第i通道的输出电压,含正负号;
(2)分别通过公式ΔLi=Qi/(Ki·d33)和fi=Qi/d33,将八个压电陶瓷两极产生的电荷量转换为各预紧支链的压电陶瓷的动态变形量和各预紧支链的压电陶瓷所受到的压力;式中,ΔLi为压电陶瓷的动态变形量,fi为各预紧支链的压电陶瓷所受到的压力,即预紧支链所受的轴向力,ΔLi和fi的正负号由Qi决定;Ki为第i个压电陶瓷的等效刚度,d33为第i组压电陶瓷的等效压电系数;
(3)假设立方体质量块的质量和半边长分别为m和n,根据牛顿-欧拉法构建传感器系统的动力学方程组,
其中,立方体质量块的质心的绝对线加速度aI表示为
Figure FDA0002383049780000031
质量块的绝对角加速度εI表示为
Figure FDA0002383049780000032
式中,
Figure FDA0002383049780000033
为{W}系相对于{I}系的旋转矩阵,通过引入辅助角速度ω*的方式求解得出;
(4)将八条预紧支链产生的轴向力fi
Figure FDA0002383049780000034
代入系统的动力学方程组,即求得立方体质量块的绝对线加速度aI和绝对角加速度εI,也即待测载体的实时加速度。
7.一种如权利要求1~4任一所述预紧式并联六维加速度传感器的灵敏度分析方法,其特征在于,灵敏度包括线加速度灵敏度和角加速度灵敏度,依次推导传感器沿x、y和z三个坐标轴方向的线加速度灵敏度和角加速度灵敏度;将沿三个坐标轴方向的线加速度灵敏度的最低值作为传感器的线加速度灵敏度,将沿三个坐标轴方向的角加速度灵敏度的最低值作为传感器的角加速度灵敏度。
8.根据权利要求7所述的灵敏度分析方法,其特征在于,将所述传感器的灵敏度定义为:加速度分量每增加一个单位值时对应的八条预紧支链所受轴向力的变化量的绝对值之和。
9.根据权利要求8所述的灵敏度分析方法,其特征在于,构建传感器系统的六维加速度参量与八条预紧支链的轴向力之间的映射关系,得到预紧式并联六维加速度传感器所有预紧支链在纯线加速度或者纯角加速度驱动情况下的轴向力情况;推导预紧式并联六维加速度传感器灵敏度数学模型,包括如下步骤:
S1:选取外壳或与外壳刚性连接的载体为研究对象,对预紧式并联六维加速度传感器模型进行简化,构建传感器在纯线加速度运动或者纯角加速度运动情况下的动力学方程组;根据牛顿第二定律列写传感器沿x轴方向的纯线加速度ax和ax+Δax驱动下立方体质量块的两组动力学方程:
Figure FDA0002383049780000041
Figure FDA0002383049780000042
式中,ax和ax+Δax为沿x轴方向输入传感器的已知加速度参量;f1~f8为八条预紧支链的轴向力;m为立方体质量块的质量,n为立方体质量块的半边长;
S2:基于并联机构的正向运动学理论挖掘预紧式并联六维加速度传感器各个支链的变形协调条件,推导出该传感器的两个力协调方程:
Figure FDA0002383049780000043
Figure FDA0002383049780000044
式中,fi
Figure FDA0002383049780000045
分别为第i条预紧支链所受的轴向力和初始预紧力;ki,i=1~8为第i条预紧支链的轴向刚度,由于八条预紧支链结构完全相同,则ki=k,则两个协调方程分别整理为:
Figure FDA0002383049780000046
Figure FDA0002383049780000047
S3:结合式(1)、式(5)和式(6),得纯线加速度ax驱动对应的预紧支链的轴向力f,写成矩阵的形式:
f=C-1D1+C-1D2 (7)
式中,
Figure FDA0002383049780000051
f=[f1 f2 f3 f4 f5 f6 f7 f8]T
D1=[max 0 0 0 0 0 0 0]T
D2=[0 0 0 0 0 0 b1 b2]T
Figure FDA0002383049780000052
Figure FDA0002383049780000053
S4:结合式(2)、式(5)和式(6),得纯线加速度ax+Δax驱动对应的预紧支链的轴向力f*,写成矩阵的形式:
f*=C-1D1 *+C-1D2 (8)
式中,
Figure FDA0002383049780000054
S5:将式(7)和式(8)的等号两边分别相减,得到:
Figure FDA0002383049780000055
S6:根据所述传感器的灵敏度的定义,得到所述传感器沿x轴方向的线加速度灵敏度Sax的定义式
Figure FDA0002383049780000056
并结合式(9),推得传感器沿x轴方向的线加速度灵敏度
Figure FDA0002383049780000057
S7:按照步骤S1~S6推导得到所述传感器沿y轴方向的线加速度灵敏度
Figure FDA0002383049780000058
传感器沿z轴方向的线加速度灵敏度
Figure FDA0002383049780000061
传感器沿x轴方向的角加速度灵敏度
Figure FDA0002383049780000062
传感器沿y轴方向的角加速度灵敏度
Figure FDA0002383049780000063
和传感器沿z轴方向的角加速度灵敏度
Figure FDA0002383049780000064
S8:最终推导出该传感器的线加速度灵敏度为Sa=min(Sax,Say,Saz)=m,角加速度灵敏度为Sε=min(Sεx,Sεy,Sεz)=mn。
CN201810599723.6A 2018-06-11 2018-06-11 预紧式并联六维加速度传感器及其测量和灵敏度分析方法 Expired - Fee Related CN108802427B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810599723.6A CN108802427B (zh) 2018-06-11 2018-06-11 预紧式并联六维加速度传感器及其测量和灵敏度分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810599723.6A CN108802427B (zh) 2018-06-11 2018-06-11 预紧式并联六维加速度传感器及其测量和灵敏度分析方法

Publications (2)

Publication Number Publication Date
CN108802427A CN108802427A (zh) 2018-11-13
CN108802427B true CN108802427B (zh) 2020-06-02

Family

ID=64085254

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810599723.6A Expired - Fee Related CN108802427B (zh) 2018-06-11 2018-06-11 预紧式并联六维加速度传感器及其测量和灵敏度分析方法

Country Status (1)

Country Link
CN (1) CN108802427B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109765842A (zh) * 2019-01-17 2019-05-17 湖北工程学院 一种微定位工作台
CN111750902B (zh) * 2020-06-24 2021-12-24 兰州理工大学 多灵敏度压电式螺栓松动监测装置及其使用与识别方法
CN113466491B (zh) * 2021-07-01 2023-12-05 兰州空间技术物理研究所 一种卫星加速度计敏感结构
CN114624468B (zh) * 2022-05-17 2022-09-02 山东利恩斯智能科技有限公司 一种防水型六维振动传感器及其测量方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5224380A (en) * 1990-05-21 1993-07-06 The University Of Maryland Superconducting six-axis accelerometer
CN101034095A (zh) * 2007-04-17 2007-09-12 燕山大学 六维正交解耦rss加速度传感器
CN101336377A (zh) * 2005-12-29 2008-12-31 舍塞尔公司 用于获取六分量地震数据的系统
CN101493473A (zh) * 2009-03-04 2009-07-29 南京航空航天大学 基于6-ss并联机构的压电式六维加速度传感器
CN101949954A (zh) * 2010-08-10 2011-01-19 南京航空航天大学 冗余并联式六维加速度传感器及其测量方法
CN104048790A (zh) * 2014-06-06 2014-09-17 南京航空航天大学 对偶正交六维力传感器及测量方法
CN105352647A (zh) * 2015-10-22 2016-02-24 哈尔滨工业大学 一种机器人末端十二维传感器及其设计方法
CN105353168A (zh) * 2015-12-07 2016-02-24 南京林业大学 一种预紧型六维加速度传感器及其故障自诊断方法
CN107643425A (zh) * 2017-08-18 2018-01-30 南京林业大学 六维加速度传感器及故障诊断方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7178401B2 (en) * 2005-04-14 2007-02-20 General Electric Company Three axis accelerometer with variable axis sensitivity

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5224380A (en) * 1990-05-21 1993-07-06 The University Of Maryland Superconducting six-axis accelerometer
CN101336377A (zh) * 2005-12-29 2008-12-31 舍塞尔公司 用于获取六分量地震数据的系统
CN101034095A (zh) * 2007-04-17 2007-09-12 燕山大学 六维正交解耦rss加速度传感器
CN101493473A (zh) * 2009-03-04 2009-07-29 南京航空航天大学 基于6-ss并联机构的压电式六维加速度传感器
CN101949954A (zh) * 2010-08-10 2011-01-19 南京航空航天大学 冗余并联式六维加速度传感器及其测量方法
CN104048790A (zh) * 2014-06-06 2014-09-17 南京航空航天大学 对偶正交六维力传感器及测量方法
CN105352647A (zh) * 2015-10-22 2016-02-24 哈尔滨工业大学 一种机器人末端十二维传感器及其设计方法
CN105353168A (zh) * 2015-12-07 2016-02-24 南京林业大学 一种预紧型六维加速度传感器及其故障自诊断方法
CN107643425A (zh) * 2017-08-18 2018-01-30 南京林业大学 六维加速度传感器及故障诊断方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
"Single-mass 6-DOF isotropic accelerometer with segmented PSD sensors";Chapsky V 等;《sensors & actuators》;20071231;第135卷(第2期);558-569页 *
"六维加速度传感器的研究现状及发展趋势";尤晶晶 等;《振动与冲击》;20151231;第34卷(第11期);150-172页 *
"基于冗余并联机构的压电式六维加速度传感器研究";尤晶晶;《中国优秀博士学位论文全文数据库》;20161231;I140-71 *
"并联式六维加速度传感器数学建模及精度特性研究";夏玉辉;《中国优秀硕士学位论文全文数据库》;20151231;I140-302 *
"预紧式并联六维加速度传感器的解耦算法研究";尤晶晶 等;《仪器仪表学报》;20170531;第38卷(第5期);1216-1225页 *

Also Published As

Publication number Publication date
CN108802427A (zh) 2018-11-13

Similar Documents

Publication Publication Date Title
CN108802427B (zh) 预紧式并联六维加速度传感器及其测量和灵敏度分析方法
CN101949954B (zh) 冗余并联式六维加速度传感器及其测量方法
Sun et al. Design and optimization of a novel six-axis force/torque sensor for space robot
CN112630813B (zh) 基于捷联惯导和北斗卫星导航系统的无人机姿态测量方法
CN103575932B (zh) 一种mems压阻式加速度计
CN110259879B (zh) 用于电动Stewart结构的无力反馈隔振控制方法和系统
You et al. Decoupling algorithm and maximum operation frequency of a novel parallel type six-axis accelerometer
Wang et al. A six-degree-of-freedom acceleration sensing method based on six coplanar single-axis accelerometers
CN105352647A (zh) 一种机器人末端十二维传感器及其设计方法
CN101493473A (zh) 基于6-ss并联机构的压电式六维加速度传感器
CN106767777A (zh) 一种新型水下机器人冗余捷联惯导装置
CN101464472B (zh) 基于九加速度敏感单元的六轴加速度传感器的布局方法
Liu et al. Ocean wave buoy based on parallel six-dimensional accelerometer
CN105353168A (zh) 一种预紧型六维加速度传感器及其故障自诊断方法
Vihonen et al. Geometry-aided angular acceleration sensing of rigid multi-body manipulator using mems rate gyros and linear accelerometers
Sun et al. Design and optimization of a novel six-axis force/torque sensor with good isotropy and high sensitivity
CN112857366B (zh) 一种基于压缩结构的光纤捷联惯导系统姿态解算方法
CN102721831B (zh) 基于折叠梁结构的一体差动式石英振梁加速度计
CN203606386U (zh) 一种mems压阻式加速度计
KR100327602B1 (ko) 3차원적인 개체의 움직임을 감지/측정하는 방법 및 그 장치
CN109668550B (zh) 一种全解耦的三自由度微机械陀螺
Portilla et al. A sensor based on a spherical parallel mechanism for the measurement of fluid velocity: Experimental development
Li et al. A Comparative Study on Two Methods of Decoupling a Six-Axis Accelerometer without and with a Gyroscope
You et al. Research on Hamiltonian dynamics of parallel type six-axis accelerometer
CN113848780A (zh) 一种多核异构处理器架构下的高机动平台姿态解算装置及方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20181113

Assignee: Jiangsu Sinosunflower Energy Technology Co.,Ltd.

Assignor: NANJING FORESTRY University

Contract record no.: X2020980008524

Denomination of invention: Preloaded parallel six dimensional acceleration sensor and its measurement and sensitivity analysis method

Granted publication date: 20200602

License type: Common License

Record date: 20201126

Application publication date: 20181113

Assignee: Jiangsu Baihe Machinery Co.,Ltd.

Assignor: NANJING FORESTRY University

Contract record no.: X2020980008577

Denomination of invention: Preloaded parallel six dimensional acceleration sensor and its measurement and sensitivity analysis method

Granted publication date: 20200602

License type: Common License

Record date: 20201127

EE01 Entry into force of recordation of patent licensing contract
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200602

Termination date: 20210611

CF01 Termination of patent right due to non-payment of annual fee