CN106406085B - 基于跨尺度模型的空间机械臂轨迹跟踪控制方法 - Google Patents

基于跨尺度模型的空间机械臂轨迹跟踪控制方法 Download PDF

Info

Publication number
CN106406085B
CN106406085B CN201610145908.0A CN201610145908A CN106406085B CN 106406085 B CN106406085 B CN 106406085B CN 201610145908 A CN201610145908 A CN 201610145908A CN 106406085 B CN106406085 B CN 106406085B
Authority
CN
China
Prior art keywords
space manipulator
space
joint
neural network
control method
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
CN201610145908.0A
Other languages
English (en)
Other versions
CN106406085A (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN201610145908.0A priority Critical patent/CN106406085B/zh
Publication of CN106406085A publication Critical patent/CN106406085A/zh
Application granted granted Critical
Publication of CN106406085B publication Critical patent/CN106406085B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)
  • Manipulator (AREA)

Abstract

一种基于跨尺度模型的空间机械臂轨迹跟踪控制方法,在分析漂浮基空间机械臂***建模过程中存在参数与非参数跨尺度特征的情况下,对机械臂关节空间进行实时在线跟踪控制。该控制方法引入径向基神经网络对空间机械臂动力学模型中存在跨尺度特征的变化项进行逼近,利用神经网络的学***面2连杆空间机械臂为例进行仿真验证,实现了空间机械臂关节空间内对期望轨迹的快速精确跟踪。

Description

基于跨尺度模型的空间机械臂轨迹跟踪控制方法
技术领域
本发明属于智能控制与***仿真技术领域,尤其是涉及一种基于跨尺度模型的空间机械臂轨迹跟踪控制方法。
背景技术
随着空间技术的不断发展,太空探索活动进一步延伸。但太空环境具有微重力、高真空、强辐射、大温差等特点,在这样危险的环境中,采用空间机械臂协助或代替宇航员来完成大量艰巨危险的任务成为世界各空间大国的一致目标。
与地面机械臂的一个显著差异是空间机械臂的基座是运动的,是一种十分复杂的多输入-多输出的强耦合非线性时变***,使得空间机械臂的控制问题与地面机械臂相比具有许多新的特点。空间机械臂***数学模型中的跨尺度特征主要表现在参数与非参数变化的跨尺度。参数的跨尺度特征主要表现在难以精确得到的动力学、运动学参数,如各部分质心位置、转动惯量、负载质量等,同时,还有很多低阶时变参数,如随着燃料消耗基座质量的变化等。非参数的跨尺度特征无法用定常参数描述,如机械臂低速运动时的非线性摩擦力、结构共振和一些高阶未建模误差等。一些针对地面机械臂能够达到较好效果的控制方法,未必适用于空间机械臂。针对轨迹跟踪控制问题,目前常用的控制方法主要有PID控制、自适应控制、鲁棒控制和智能控制。
Parlaktuna和Ozkan将漂浮基空间机械臂的控制问题由惯性空间转化到关节空间内,在得到空间机械臂***参数线性化的动力学方程后,设计了一种关节空间内漂浮基空间机械臂轨迹跟踪的PD控制方法。但是PID控制属于线性控制方法,忽略了空间机械臂***中的非线性因素和外界干扰,同时,PID控制往往需要较大的控制能量,不适用于跟踪精度要求较高的情况。所以当空间机械臂***存在参数或非参数跨尺度特征时,控制效果并不理想。
Wang H和Xie Y针对漂浮基空间机械臂设计了一种递推自适应控制方法,利用参数自适应律实时估计控制参数;梁捷、陈力设计了一种标称计算力矩控制附加自适应模糊补偿控制的复合控制方法,能够有效地克服空间机械臂***未知参数的影响;张福海等设计了一种笛卡儿空间内的自适应轨迹跟踪控制方法,在保证惯量矩阵可逆的同时,又可以实时估计控制参数。然而上述自适应控制方法仅有效地克服了参数变化对空间机械臂***的影响,当漂浮基空间机械臂***存在外部扰动等非参数的跨尺度特征时,单纯地采用自适应控制方法难以保证空间机械臂***的稳定性,需要与其它先进的控制策略相结合提高空间机械臂***鲁棒性。
谢立敏等针对空间机械臂关节控制输入力矩幅值受限且空间机械臂***存在不确定参数的复杂情况,设计了一种鲁棒自适应混合控制方法,对参数进行鲁棒自适应调节;Pazelli等针对存在参数变化影响和外部干扰的漂浮基空间机械臂***,对各类非线性H控制方法进行研究分析。然而上述鲁棒控制方法是以先验知识上界为基础设计的,是一种比较保守的控制策略,因此不是最佳控制。
郭益深和陈力利用径向基神经网络,提出了一种无需机械臂动力学模型的自适应神经网络控制方法,但并没有讨论模型存在跨尺度特征时的解决方法;谢箭等提出了一种针对漂浮基空间机械臂的神经网络自适应控制方法,通过径向基神经网络逼近模型的非线性函数和不确定性上界,提出的自适应控制律保证了权值的有界性,但是所设计的自适应律较为复杂,影响计算速度;张文辉等设计了一种径向基神经网络鲁棒自适应控制方法,应用于漂浮基空间机械臂***,雷霆针对控制力矩受限的漂浮基空间机械臂***,设计了一种神经网络自适应控制方法,然而这两种方法针对参数变化所设计的补偿律包含了动力学模型的全部信息,其中模型的标称部分为已知信息,在补偿律中属于冗余部分。
发明内容
本发明的目的是提供一种神经网络自适应控制方法,针对动力学模型中具有参数及非参数跨尺度的空间机械臂***,实现关节空间对期望轨迹的快速精确跟踪。
为实现上述目的,本发明提供一种基于跨尺度模型的空间机械臂轨迹跟踪控制方法,其特征在于:设计神经网络自适应控制律,将空间机械臂***动力学模型中存在跨尺度特征的变化项表示为利用神经网络对变化项f进行逼近,从而实现对变化项f的补偿,神经网络控制律为v为用于克服神经网络逼近误差的鲁棒项,其值为v=Kvsgn(r);误差函数为神经网络的形式为径向基神经网络,径向基神经网络的输入取理想的逼近算法为则网络的输出为径向基神经网络权值调整自适应律为D0、C0为对象的名义模型,D0=D-ΔD,C0=C-ΔC,ΔD、ΔC为建模误差矩阵,d是总和扰动,e=qd-q和分别为关节角跟踪误差和角速度跟踪误差,qd和q分别为期望和实际的关节矢量,KP、KI分别是正定比例和积分增益矩阵,Kv为鲁棒项系数,为神经网络的权值向量,为高斯基函数的输出向量,ci为网络第i个节点的中心矢量,bi为节点i的基宽度参数。
与现有技术相比本发明的有益效果是:
考虑到空间机械臂***中的建模误差和外界干扰,用径向基神经网络对空间机械臂***动力学模型中存在跨尺度特征的参数及非参数项f进行在线逼近,f仅包括建模误差ΔD(q)、以及未知干扰无需考虑已知的标称模型。利用神经网络的学习能力,有效抑制了参数与非参数跨尺度变化对空间机械臂***的影响,自适应律可在线调整神经网络权值,保证了权值的有界性,解决了未知上界有界的问题。
本发明仿真结果得出,在2s内,关节1和关节2的角位移、角速度迅速跟踪上期望轨迹。关节角跟踪误差稳定在±5×10-3rad以内,关节角速度误差稳定在±5×10-3rad/s以内,实现空间机械臂关节空间内对期望轨迹的快速精确跟踪。
附图说明
图1是平面2连杆空间机械臂模型图;
图2是神经网络自适应控制方法结构框图;
图3是关节1角度跟踪随时间变化的曲线;
图4是关节2角度跟踪随时间变化的曲线;
图5是关节1角速度跟踪随时间变化的曲线;
图6是关节2角速度跟踪随时间变化的曲线;
图7是关节1关节角跟踪误差随时间变化的曲线;
图8是关节2关节角跟踪误差随时间变化的曲线;
图9是关节1关节角速度跟踪误差随时间变化的曲线;
图10是关节2关节角速度跟踪误差随时间变化的曲线;
图11是关节1的控制力矩随时间变化的曲线;
图12是关节2的控制力矩随时间变化的曲线。
其中:∑I,惯性坐标系;∑0,运动基座坐标系;O,惯性坐标系原点;CM,空间机械臂***总质心;Bi,刚体i,机械臂的第i个连杆;B0,运动基座,连杆0;Ci,连杆i的质心;C0,基座的质心;ri∈R2,连杆i质心的位置矢量;r0,基座质心的位置矢量;rc∈R2,空间机械臂***质心CM的位置矢量;pi∈R2,连杆i的位置矢量;ai,从关节Ji到连杆i质心的矢量;bi,从连杆i质心到关节Ji+1的矢量;b0,从基座质心到关节J1的矢量。
具体实施方式
平面2连杆空间机械臂模型如图1,由可自由漂浮的运动基座B0和两个臂杆B1、B2组成。
空间机械臂***各项动力学参数如表1所描述,由基座的初始位置和姿态角和连杆1、连杆2的初始姿态角组成的矢量为[qb,qs]T=[x,y,q0,q1,q2]T,基座和连杆1、连杆2的初始速度矢量为各项参数及期望轨迹的初始数值如表2所示。
表1平面2连杆空间机械臂***参数表
表2空间机械臂神经网络自适应控制仿真初始值
设置控制参数为KP=diag{100,100,100,100,100},KI=diag{250,250,250,250,250},Kv=0.2,FW=diag{0.0005,0.0005,0.0005,0.0005,0.0005},根据图2的控制方法结构框图进行仿真验证,本发明以平面2连杆漂浮基空间机械臂***为研究对象,动力学方程为
其中,q=[q1 q2]T为关节角位移量,D(q)5×5为空间机械臂的惯性矩阵,表示包括非线性离心力和哥氏力的矩阵,τ为控制力矩。
空间机械臂的动力学方程满足如下性质:
性质1惯性矩阵D(q)是对称、正定、有界矩阵。
性质2选择适当的可使D(q)和满足
性质3存在kc>0及正定函数使得
性质4给定误差矩阵满足ΔDl≤||ΔD||≤ΔDh,ΔCl≤||ΔC||≤ΔCh,其中h和l分别为上下界值。
为了实现关节空间内对期望轨迹的快速精确跟踪,建模时还需考虑空间机械臂***存在参数与非参数跨尺度特征。引入外部扰动,可以将动力学方程式改写为如下形式
其中,是总和扰动,包括摩擦力矩扰动和其他外部扰动。
在实际工程中,对象的实际模型很难得到,即无法得到精确的D(q)、只能建立理想的名义模型。将动力学方程写成理想模型和存在跨尺度特征的变化项之和的形式,则可以表示为
其中,D0、C0为对象的名义模型,D0(q)=D(q)-ΔD(q),ΔD(q)、为误差矩阵,为包括建模误差和外部干扰力矩等存在跨尺度特征的参数与非参数项,为一未知非线性时变函数,具体形式为
由于D(q)是可逆的,可得
设计误差函数为
其中,e=qd-q和分别为关节角跟踪误差和加速度跟踪误差,qd和q分别为期望和实际的关节矢量,KP、KI分别是正定比例和积分增益矩阵。
求导可得
导出等效控制律
因此得到稳定的闭环***为
针对名义模型,控制律设计为
考虑到未知干扰,可以得出
由此可见,模型中存在跨尺度特征的变化项为
采用径向基神经网络逼近f,网络输入取则网络输出为
设计控制律为
其中,v=Kvsgn(r)为鲁棒项,用于克服神经网络逼近误差造成的影响。
定义Lyapunov函数为
其中,D和FW为正定阵,即则V是正定的。
求导,并结合性质2可得
整理可得
根据条件考虑到v=Kvsgn(r),可得
则可设计如下自适应律以调整径向基神经网络的权值
于是可得
为了方便描述,定义
于是
由性质3和性质4可取KvIn>|Q|,其中In=[1,1,…,1]T∈Rn,则
控制结构的输入为关节角的期望轨迹,输出关节角的实际值作为负反馈与关节角的实际值作比较,根据关节角的跟踪误差、误差函数和鲁棒项设计神经网络自适应控制***,仿真结果如图3-图12所示。
关节1和关节2的角度随时间变化曲线如图3、图4,角速度随时间变化曲线如图5、图6。其中的红色虚线表示轨迹跟踪的期望值,蓝色实线表示实际的关节矢量。可以看出,关节1和关节2的角位移、角速度在2s内迅速跟踪上期望轨迹。
关节1和关节2的角度跟踪误差随时间变化曲线如图7、图8,可以看出,关节角跟踪误差保持在±5×10-3rad范围内。
关节1和关节2的角速度跟踪误差随时间变化曲线如图9、图10,可以看出,关节角速度跟踪误差保持在±5×10-3rad/s范围内。
关节1和关节2的控制力矩随时间变化曲线如图11、图12,可以看出,各关节控制力矩保持在可实现的范围内。
运动基座的位置及姿态角随时间变化的曲线,可以看出,连杆1和连杆2的运动造成基座位置和姿态的变化较为平缓,适用于漂浮基空间机械臂的轨迹跟踪控制。
实验结果验证了神经网络自适应控制算法的有效性,对于动力学模型中存在参数与非参数跨尺度特征的情况,能够快速地在线跟踪期望轨迹,具有一定的鲁棒性和抗干扰性。

Claims (1)

1.一种基于跨尺度模型的空间机械臂轨迹跟踪控制方法,其特征在于:设计神经网络自适应控制律,将空间机械臂***动力学模型中存在跨尺度特征的变化项表示为利用神经网络对变化项f进行逼近,从而实现对变化项f的补偿,神经网络控制律为v为用于克服神经网络逼近误差的鲁棒项,其值为v=Kvsgn(r);误差函数为神经网络的形式为径向基神经网络,径向基神经网络的输入取理想的逼近算法为则网络的输出为径向基神经网络权值调整自适应律为D0、C0为对象的名义模型,D0=D-ΔD,C0=C-ΔC,ΔD、ΔC为建模误差矩阵,d是总和扰动,e=qd-q和分别为关节角跟踪误差和角速度跟踪误差,qd和q分别为期望和实际的关节矢量,KP、KI分别是正定比例和积分增益矩阵,Kv为鲁棒项系数,为神经网络的权值向量,为高斯基函数的输出向量,ci为网络第i个节点的中心矢量,bi为节点i的基宽度参数,D为空间机械臂的惯性矩阵,C表示包括非线性离心力和哥氏力的矩阵,FW为正定阵。
CN201610145908.0A 2016-03-15 2016-03-15 基于跨尺度模型的空间机械臂轨迹跟踪控制方法 Expired - Fee Related CN106406085B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610145908.0A CN106406085B (zh) 2016-03-15 2016-03-15 基于跨尺度模型的空间机械臂轨迹跟踪控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610145908.0A CN106406085B (zh) 2016-03-15 2016-03-15 基于跨尺度模型的空间机械臂轨迹跟踪控制方法

Publications (2)

Publication Number Publication Date
CN106406085A CN106406085A (zh) 2017-02-15
CN106406085B true CN106406085B (zh) 2019-02-01

Family

ID=58007139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610145908.0A Expired - Fee Related CN106406085B (zh) 2016-03-15 2016-03-15 基于跨尺度模型的空间机械臂轨迹跟踪控制方法

Country Status (1)

Country Link
CN (1) CN106406085B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107203141A (zh) * 2017-08-02 2017-09-26 合肥工业大学 一种机械臂分散化神经鲁棒控制的轨迹跟踪算法
CN107505835A (zh) * 2017-08-11 2017-12-22 广东工业大学 一种控制机械手运动的方法、装置及***
CN108445768B (zh) * 2018-05-29 2020-12-25 福州大学 空间机器人操作空间轨迹跟踪的增广自适应模糊控制方法
CN109176529B (zh) * 2018-10-19 2021-11-30 福州大学 一种空间机器人协调运动的自适应模糊控制方法
CN109227550A (zh) * 2018-11-12 2019-01-18 吉林大学 一种基于rbf神经网络的机械臂控制方法
CN111061216B (zh) * 2019-12-28 2022-11-15 哈尔滨工业大学 基于二元样条尺度函数的贴片机运动***智能控制方法
CN111399397B (zh) * 2020-04-01 2022-03-04 合肥工业大学 机器人的控制方法、控制器及控制***
CN112152539A (zh) * 2020-09-29 2020-12-29 中国船舶重工集团公司第七二四研究所 一种神经网络补偿电机负载力矩观测器实现方法
CN113219825B (zh) * 2021-03-26 2023-04-25 齐鲁工业大学 一种四足机器人单腿轨迹跟踪控制方法及***
CN113253610B (zh) * 2021-04-20 2021-10-01 中国科学院自动化研究所 飞行器控制方法及装置
CN113370205B (zh) * 2021-05-08 2022-06-17 浙江工业大学 一种基于机器学习的Baxter机械臂轨迹跟踪控制方法
CN113296393A (zh) * 2021-05-27 2021-08-24 安徽工业大学 基于自调整模糊迭代学习的二连杆机械臂控制方法、装置及介质
CN114516047B (zh) * 2022-02-14 2024-07-09 安徽大学 基于径向基神经网络终端滑模控制机械臂轨迹方法及***
CN114700938B (zh) * 2022-03-04 2023-06-16 华南理工大学 一种基于跳增益积分神经网络的冗余机械臂运动规划方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5428559A (en) * 1991-08-14 1995-06-27 Kabushiki Kaisha Toshiba Predictive control method and apparatus
JPH07223181A (ja) * 1994-02-09 1995-08-22 Yaskawa Electric Corp マニピュレータの干渉チェック装置
CN101630144A (zh) * 2009-08-18 2010-01-20 湖南大学 电子节气门的自学习逆模型控制方法
CN102289204A (zh) * 2011-06-03 2011-12-21 华南理工大学 基于确定学习理论的机械臂通用控制方法
CN102501251A (zh) * 2011-11-08 2012-06-20 北京邮电大学 带有动态摩擦补偿的机械臂关节定位控制方法
CN104589349A (zh) * 2015-01-16 2015-05-06 西北工业大学 一种混合悬浮微重力环境下带有单关节机械臂的组合体自主控制方法
CN105223808A (zh) * 2015-06-24 2016-01-06 浙江工业大学 基于神经网络动态面滑模控制的机械臂***饱和补偿控制方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5428559A (en) * 1991-08-14 1995-06-27 Kabushiki Kaisha Toshiba Predictive control method and apparatus
JPH07223181A (ja) * 1994-02-09 1995-08-22 Yaskawa Electric Corp マニピュレータの干渉チェック装置
CN101630144A (zh) * 2009-08-18 2010-01-20 湖南大学 电子节气门的自学习逆模型控制方法
CN102289204A (zh) * 2011-06-03 2011-12-21 华南理工大学 基于确定学习理论的机械臂通用控制方法
CN102501251A (zh) * 2011-11-08 2012-06-20 北京邮电大学 带有动态摩擦补偿的机械臂关节定位控制方法
CN104589349A (zh) * 2015-01-16 2015-05-06 西北工业大学 一种混合悬浮微重力环境下带有单关节机械臂的组合体自主控制方法
CN105223808A (zh) * 2015-06-24 2016-01-06 浙江工业大学 基于神经网络动态面滑模控制的机械臂***饱和补偿控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Research on Control Method for Time Delay of Internet-Based Tele-Operation Manipulators;Wei Gao 等;《International Journal of Future Generation Communication and Networking》;20130831;第6卷(第4期);第203-216页
Trajectory planning of a manipulator of a space robot using a neural network;T. Mineta 等;《IEEE》;20020806;第34-39页

Also Published As

Publication number Publication date
CN106406085A (zh) 2017-02-15

Similar Documents

Publication Publication Date Title
CN106406085B (zh) 基于跨尺度模型的空间机械臂轨迹跟踪控制方法
CN108942924B (zh) 基于多层神经网络的模型不确定性机械臂运动控制方法
CN105549598B (zh) 一种二维运动移动机器人的迭代学习轨迹跟踪控制及其鲁棒优化方法
CN109465825A (zh) 机械臂柔性关节的rbf神经网络自适应动态面控制方法
CN102411304B (zh) 一种航天器小角度姿态机动控制参数优化方法
CN107160398B (zh) 基于确定学习的全状态受限刚性机械臂安全可靠控制方法
CN109176525A (zh) 一种基于rbf的移动机械手自适应控制方法
Zhou et al. Adaptive NN control for nonlinear systems with uncertainty based on dynamic surface control
CN105785762A (zh) 一种基于自适应反步滑模的两轴惯性稳定平台高精度控制方法
CN110181510A (zh) 一种基于时延估计与模糊逻辑的机械臂轨迹跟踪控制方法
CN109683624A (zh) 用于小型无人直升机姿态控制的非线性鲁棒控制方法
CN108445768A (zh) 空间机器人操作空间轨迹跟踪的增广自适应模糊控制方法
Alqaisi et al. Adaptive sliding mode control based on RBF neural network approximation for quadrotor
CN105652667A (zh) 一种模型不确定双关节机械手的高精度轨迹跟踪控制方法
Qiu et al. Model predictive control for uncalibrated and constrained image-based visual servoing without joint velocity measurements
Sainzaya et al. LQR control with refined PID to balance rotary inverted pendulum with time-varying uncertainty
Kizir et al. Fuzzy control of a real time inverted pendulum system
CN108693776A (zh) 一种三自由度Delta并联机器人的鲁棒控制方法
Guo et al. The robot arm control based on rbf with incremental pid and sliding mode robustness
Niu et al. The adaptive control based on BP neural network identification for two-wheeled robot
CN112363398B (zh) 一种控制输入受限下的桥吊***有限时间滑模控制***及方法
Martins et al. Trajectory tracking of a nonholonomic mobile robot with parametric and nonparametric uncertainties: A proposed neural control
Benhellal et al. Decoupled adaptive neuro-fuzzy sliding mode control applied in a 3D crane system
Yan-dong et al. Adaptive RBFNN formation control of multi-mobile robots with actuator dynamics
Hu et al. An efficient neural controller for a nonholonomic mobile robot

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190201

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