CN106402089A - 一种基于耦合干扰观测器的级联电液伺服系统控制方法及系统 - Google Patents

一种基于耦合干扰观测器的级联电液伺服系统控制方法及系统 Download PDF

Info

Publication number
CN106402089A
CN106402089A CN201610924448.1A CN201610924448A CN106402089A CN 106402089 A CN106402089 A CN 106402089A CN 201610924448 A CN201610924448 A CN 201610924448A CN 106402089 A CN106402089 A CN 106402089A
Authority
CN
China
Prior art keywords
cascade
centerdot
hydraulic
alpha
servo
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.)
Granted
Application number
CN201610924448.1A
Other languages
English (en)
Other versions
CN106402089B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201610924448.1A priority Critical patent/CN106402089B/zh
Publication of CN106402089A publication Critical patent/CN106402089A/zh
Application granted granted Critical
Publication of CN106402089B publication Critical patent/CN106402089B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F15FLUID-PRESSURE ACTUATORS; HYDRAULICS OR PNEUMATICS IN GENERAL
    • F15BSYSTEMS ACTING BY MEANS OF FLUIDS IN GENERAL; FLUID-PRESSURE ACTUATORS, e.g. SERVOMOTORS; DETAILS OF FLUID-PRESSURE SYSTEMS, NOT OTHERWISE PROVIDED FOR
    • F15B21/00Common features of fluid actuator systems; Fluid-pressure actuator systems or details thereof, not covered by any other group of this subclass
    • F15B21/02Servomotor systems with programme control derived from a store or timing device; Control devices therefor
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]

Abstract

本发明公开了一种基于耦合干扰观测器的级联电液伺服系统控制方法,液压控制系统中高动态的外负载干扰,且外负载干扰与级联机电对象的多个机构运动相关,本发明采用基于耦合干扰观测器的反步控制方法来提高存在外负载干扰条件下电液伺服控制驱动级联机械臂的动态控制性能。电液伺服系统中的负载干扰采用耦合干扰观测器进行在线估计,然后设计反步控制律实现级联机械臂的位置跟踪。由于反步控制律中含有干扰估计值,对干扰进行实时补偿,从一定程度上抑制控制量恶化,并提高了级联机械臂关节运动的动态控制性能。

Description

一种基于耦合干扰观测器的级联电液伺服系统控制方法及 系统
技术领域
本发明涉及一种可用于级联电液伺服系统的控制方法,特别是级联电液伺服驱动多自由度机械臂、起重装置等机电对象。
背景技术
目前应用在多自由度机械臂、多关节机器人等机电对象上的控制执行机构一般为电液伺服执行器,由于其能输出有很高的控制力和力矩。机电对象最常见的驱动模式是全驱动模式。而在全驱动模式中,一般将考虑关节负载作为干扰量处理,不需要对其进行测量或估计,关节跟踪控制算法也不包含关节负载的估计值,这样将导致大动态负载干扰对关节运动跟踪的动态响应性能的影响很大,控制精度降低。因此,本专利设计一种耦合干扰观测器,估计多自由度机械臂各关节负载,并将估计值应用于级联电液伺服控制算法中,从而抑制控制变量剧烈恶化,提高大动态负载干扰情况下级联电液伺服控制系统的输出响应性能,并保证系统具有全局稳定性和有界收敛性。
发明内容
本发明的目的是克服目前级联电液伺服全驱动控制方法的不足,同时适用于高动态负载干扰的在线估计,可以防止控制变量剧烈恶化,提高级联电液伺服控制系统的跟踪动态性能。
本发明的技术方案是一种基于耦合干扰观测器的级联电液伺服系统控制方法,该方法包括:
步骤1:建立级联电液伺服执行器模型;
建立的级联电液伺服执行器模型为:
yi=xi1
i=1,…n
其中为第i个级联电液伺服执行器的4阶模型状态变量,yi为第i个液压缸输出位移,为第i个输出位移变化率,pLi为第i个液压缸输出的液压压力,xvi为第i个伺服阀阀芯位移,mi为第i个负载质量,ps为供油压力,Ap为对称缸横截面积,Ctl为液压缸总泄漏系数,Vt为液压缸容积,βe为液压油有效体积弹性模量,Cd为伺服阀流量系数,w为伺服阀面积梯度,ρ为液压油密度,K为负载刚度系数,b为液压油阻尼系数,FLi为第i个外负载压力,Ksv为伺服阀放大系数,Tsv为伺服阀一阶响应时间常数,k为反正切函数tanh(·)中指数项系数,ui为第i个伺服阀控制电压,Tsvi为第i个伺服阀一阶响应时间常数,Ksvi为第i个伺服阀放大系数;
步骤2:驱动电液伺服,实时获取电液伺服的反馈数据,包括:第i个液压缸输出位移、第i个液压缸输出位移变化率、第i个液压缸负载压力、第i个伺服阀阀芯位移,i=1,…n;
步骤3:利用反馈数据结合系统的状态误差计算控制变量;
步骤4:采用耦合干扰观测器对机械臂各关节负载进行在线估计;
步骤5:结合反馈数据、系统误差和负载干扰估计量计算反步控制律;
步骤6:根据反步控制律对级联电液伺服机构实时进行驱动。
进一步的,所述步骤3中系统误差zij(i=1,…,n,j=1,…,4)表示为
其中yid表示第i个液压缸期望位移指令,αij为反步控制律设计中虚拟控制变量,表示为:
其中
kij表示反步控制律中的控制参数,表示第i个伺服阀的耦合干扰观测器。
进一步的,所述步骤4中耦合干扰观测器表示如下:
其中状态变化率表示为
耦合干扰观测器参数为Kdij=-Kdji≠0,(i≠j),即第i个关节负载干扰估计值不仅与本执行器的状态变量估计误差xi2有关,而且也与其他n-1个执行器的状态变量估计误差xj2(j=1,…,n,j≠i)有关。
进一步的,所述步骤5中反步控制律为:
其中i=1,…,n。
一种采用基于耦合干扰观测器的级联电液伺服系统,该系统为多自由度机械臂对象,包括n个机械连杆,n个电液伺服阀,n个双作用液压缸,1个伺服电机,1个定量柱塞泵,1个油箱;其中第i个连杆与第i+1个连杆之间铰接,称为机械臂第i个关节,在关节处设置光电编码器,用于测量第i个关节的运动角度和角速度,在第i个液压缸进油口和出油口设置第i个压力传感器,测量第i个液压缸的负载力,在定量柱塞泵出口安装1个压力表,监测系统的供油压力。
本发明的目的之三是提出耦合干扰观测器与非线性反步控制算法相结合的反步控制设计方法,既能对机械臂多个关节处的动态负载进行实时估计,同时也可以利用伺服阀控制电压对负载进行有效补偿,并约束液压缸输出位置误差,提高级联电液伺服控制系统的跟踪动态性能。
附图说明
图1为本发明的采用基于耦合干扰观测器和反步控制方法的电液伺服驱动级联机械臂机构示意图;
图2为本发明虚拟变量计算顺序示意图;
图3为本发明一种基于耦合干扰观测器的级联电液伺服控制方法流程图。
具体实施方式
以下提供本发明一种基于耦合干扰观测器和反步控制方法的电液伺服驱动级联机械臂机构的具体实时方式。
电液伺服执行器的模型为4阶模型,不考虑机械臂机构运动的模型,机械臂运动所需要的关节力矩作为电液伺服执行器的负载干扰考虑,简述如下:
1)电液伺服执行器建模
采用四阶模型描述伺服阀驱动液压缸回路的电液伺服执行器模型如下:
i=1,…n
其中为第i个级联电液伺服执行器的4阶模型状态变量,yi为第i个液压缸输出位移,为第i个输出位移变化率,pLi为第i个液压缸输出的液压压力,xvi为第i个伺服阀阀芯位移,mi为第i个负载质量,ps为供油压力,Ap为对称缸横截面积,Ctl为液压缸总泄漏系数,Vt为液压缸容积,βe为液压油有效体积弹性模量,Cd为伺服阀流量系数,w为伺服阀面积梯度,ρ为液压油密度,K为负载刚度系数,b为液压油阻尼系数,FLi为第i个外负载压力,Ksv为伺服阀放大系数,Tsv为伺服阀一阶响应时间常数,k为反正切函数tanh(·)中指数项系数,ui为第i个伺服阀控制电压。
2)级联机械臂关节处产生的负载干扰观测器表示如下:
其中为第i个干扰估计值,为状态变化率的估计值,Kdij=-Kdji≠0,(i≠j)为耦合干扰观测器参数,实际状态变化率由测量值xi2微分获得。
3)基于耦合干扰观测器的级联反步控制律
其中
i=1,…,n
kij表示反步控制律中的控制参数,系统误差zij(i=1,…,n,j=1,…,4)表示为
yid表示第i个液压缸期望位移指令,αij为反步控制律设计中虚拟控制变量。
一般的反步控制律设计是一个迭代过程,如图2所示。每个变量的计算顺序为: 根据公式(4)可以对α2、α3求导得到变化率然后联合公式(1),(2),(3),(4),(5),(6)得到最终的反步迭代控制律表示为
其中

Claims (5)

1.一种基于耦合干扰观测器的级联电液伺服系统控制方法,该方法包括:
步骤1:建立级联电液伺服执行器模型;
建立的级联电液伺服执行器模型为:
x · i 1 = x 2 x · i 2 = 1 m i ( - Kx i 1 - bx i 2 + A p x i 3 - F L i ) x · i 3 = - 4 β e A p V t x i 2 - 4 β e C t l V t x i 3 + 4 β e C d wK s v V t ρ p s - tanh ( kx i 4 ) x i 3 x i 4 x · i 4 = - 1 T s v i x i 4 + K s v i T s v i u i
yi=xi1
i=1,…n
其中为第i个级联电液伺服执行器的4阶模型状态变量,yi为第i个液压缸输出位移,为第i个输出位移变化率,pLi为第i个液压缸输出的液压压力,xvi为第i个伺服阀阀芯位移,mi为第i个负载质量,ps为供油压力,Ap为对称缸横截面积,Ctl为液压缸总泄漏系数,Vt为液压缸容积,βe为液压油有效体积弹性模量,Cd为伺服阀流量系数,w为伺服阀面积梯度,ρ为液压油密度,K为负载刚度系数,b为液压油阻尼系数,FLi为第i个外负载压力,Ksv为伺服阀放大系数,Tsv为伺服阀一阶响应时间常数,k为反正切函数tanh(·)中指数项系数,ui为第i个伺服阀控制电压,Tsvi为第i个伺服阀一阶响应时间常数,Ksvi为第i个伺服阀放大系数;
步骤2:驱动电液伺服,实时获取电液伺服的反馈数据,包括:第i个液压缸输出位移、第i个液压缸输出位移变化率、第i个液压缸负载压力、第i个伺服阀阀芯位移,i=1,…n;
步骤3:利用反馈数据结合系统的状态误差计算控制变量;
步骤4:采用耦合干扰观测器对机械臂各关节负载进行在线估计;
步骤5:结合反馈数据、系统误差和负载干扰估计量计算反步控制律;
步骤6:根据反步控制律对级联电液伺服机构实时进行驱动。
2.如权利要求1所述的一种基于耦合干扰观测器的级联电液伺服系统控制方法,其特征在于所述步骤3中系统误差zij(i=1,…,n,j=1,…,4)表示为
z i 1 = x i 1 - y i d z i j = x i j - α i ( j - 1 ) , ( i = 1 , ... , n , j = 2 , 3 , 4 )
其中yid表示第i个液压缸期望位移指令,αij为反步控制律设计中虚拟控制变量,表示为:
α i 1 = - k i 1 z i 1 + y · i d α i 2 = 1 g i 1 ( - k i 2 z i 2 - f i 1 - z i 1 + d ^ i + α · i 1 ) α i 3 = 1 g i 2 ( - k i 3 z i 3 - g i 1 z i 2 - f i 2 + α · i 2 ) , i = 1 , ... , n
其中
f i 1 = 1 m i ( - Kx i 1 - bx i 2 ) , g i 1 = A p m i
f i 2 = - 4 β e A p V t x i 2 - 4 β e C t l V t x i 3 , g i 2 = 4 β e C d wK s v V t ρ p s - tanh ( k u ) x i 3 ,
i=1,…,n
kij表示反步控制律中的控制参数,表示第i个伺服阀的耦合干扰观测器。
3.如权利要求1所述的一种基于耦合干扰观测器的级联电液伺服系统控制方法,其特征在于所述步骤4中耦合干扰观测器表示如下:
d ^ 1 = - K d 11 ( x 12 - x ^ 12 ) - ... - K d 1 i ( x i 2 - x ^ i 2 ) ... - K d 1 n ( x n 2 - x ^ n 2 ) . . . d ^ i = - K d i 1 ( x 12 - x ^ 12 ) - ... - K d i i ( x i 2 - x ^ i 2 ) ... - K d i n ( x n 2 - x ^ n 2 ) . . . d ^ n = - K d n 1 ( x 12 - x ^ 12 ) - ... - K d n i ( x i 2 - x ^ i 2 ) ... - K d n n ( x n 2 - x ^ n 2 ) ,
其中状态变化率表示为
x ^ · 12 = 1 m 1 ( - Kx 11 - bx 12 + A p x 13 ) - d ^ 1 . . . x ^ · i 2 = 1 m i ( - Kx i 1 - bx i 2 + A p x i 3 ) - d ^ i . . . x ^ · n 2 = 1 m n ( - Kx n 1 - bx n 2 + A p x n 3 ) - d ^ n ,
耦合干扰观测器参数为Kdij=-Kdji≠0,(i≠j),即第i个关节负载干扰估计值不仅与本执行器的状态变量估计误差xi2有关,而且也与其他n-1个执行器的状态变量估计误差xj2(j=1,…,n,j≠i)有关。
4.如权利要求1所述的一种基于耦合干扰观测器的级联电液伺服系统控制方法,其特征在于所述步骤5中反步控制律为:
α i 1 = - k i 1 z i 1 + y · i d α i 2 = 1 g i 1 ( - k i 2 z i 2 - f i 1 - z i 1 + d ^ i + α · i 1 ) α i 3 = 1 g i 2 ( - k i 3 z i 3 - g i 1 z i 2 - f i 2 + α · i 2 ) , i = 1 , ... , n
u i = 1 g i 3 ( - k i 4 z i 4 - g i 2 z i 3 - f i 3 + α · i 3 )
其中
5.一种采用基于耦合干扰观测器的级联电液伺服系统,该系统为多自由度机械臂对象,包括n个机械连杆,n个电液伺服阀,n个双作用液压缸,1个伺服电机,1个定量柱塞泵,1个油箱;其中第i个连杆与第i+1个连杆之间铰接,称为机械臂第i个关节,在关节处设置光电编码器,用于测量第i个关节的运动角度和角速度,在第i个液压缸进油口和出油口设置第i个压力传感器,测量第i个液压缸的负载力,在定量柱塞泵出口安装1个压力表,监测系统的供油压力。
CN201610924448.1A 2016-10-24 2016-10-24 一种基于耦合干扰观测器的级联电液伺服系统控制方法及系统 Expired - Fee Related CN106402089B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610924448.1A CN106402089B (zh) 2016-10-24 2016-10-24 一种基于耦合干扰观测器的级联电液伺服系统控制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610924448.1A CN106402089B (zh) 2016-10-24 2016-10-24 一种基于耦合干扰观测器的级联电液伺服系统控制方法及系统

Publications (2)

Publication Number Publication Date
CN106402089A true CN106402089A (zh) 2017-02-15
CN106402089B CN106402089B (zh) 2017-11-17

Family

ID=58012622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610924448.1A Expired - Fee Related CN106402089B (zh) 2016-10-24 2016-10-24 一种基于耦合干扰观测器的级联电液伺服系统控制方法及系统

Country Status (1)

Country Link
CN (1) CN106402089B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110081046A (zh) * 2019-05-27 2019-08-02 电子科技大学 一种基于反步控制的多电液伺服执行器跟踪同步控制方法
CN110107563A (zh) * 2019-05-27 2019-08-09 电子科技大学 存在负载干扰情况下多电液伺服执行器分布协同控制方法
CN110173487A (zh) * 2019-05-27 2019-08-27 电子科技大学 一种多电液伺服执行器在切换网络下的领导同步控制方法
CN110262226A (zh) * 2019-03-29 2019-09-20 浙江理工大学 一种欠驱动系统控制方法、装置、设备及介质
CN111258220A (zh) * 2020-01-21 2020-06-09 重庆邮电大学 一种基于干扰观测器的柔性机械臂级联控制方法及系统
CN111775142A (zh) * 2020-08-12 2020-10-16 电子科技大学 一种液压机械臂的模型辨识与自适应控制方法
CN112049837A (zh) * 2020-08-05 2020-12-08 东北大学 基于参数化建模技术的液压装置数字化样机设计方法
CN112068435A (zh) * 2020-09-09 2020-12-11 北京航空航天大学 基于干扰观测器的康复机械装置迭代学习控制方法及系统
CN112180742A (zh) * 2020-11-09 2021-01-05 山东大学 基于分数阶扰动观测器的海上起重机稳定控制方法及系统
CN114801301A (zh) * 2022-04-28 2022-07-29 重庆智能机器人研究院 伺服电缸压力机的控制方法、装置、电子设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103753526A (zh) * 2013-12-30 2014-04-30 重庆交通大学 可定位补偿的精密重载机械手
CN104317198A (zh) * 2014-10-21 2015-01-28 南京理工大学 具有时变输出约束的电液伺服系统非线性鲁棒位置控制方法
CN104635490A (zh) * 2014-12-15 2015-05-20 南京理工大学 一种单出杆液压缸位置伺服系统的输出反馈控制方法
CN105159076A (zh) * 2015-08-24 2015-12-16 南京理工大学 基于融合型自适应鲁棒的电液负载模拟器力控制方法
CN105182742A (zh) * 2015-07-23 2015-12-23 北京航空航天大学 一种弹性体飞行器自适应受限跟踪控制间接法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103753526A (zh) * 2013-12-30 2014-04-30 重庆交通大学 可定位补偿的精密重载机械手
CN104317198A (zh) * 2014-10-21 2015-01-28 南京理工大学 具有时变输出约束的电液伺服系统非线性鲁棒位置控制方法
CN104635490A (zh) * 2014-12-15 2015-05-20 南京理工大学 一种单出杆液压缸位置伺服系统的输出反馈控制方法
CN105182742A (zh) * 2015-07-23 2015-12-23 北京航空航天大学 一种弹性体飞行器自适应受限跟踪控制间接法
CN105159076A (zh) * 2015-08-24 2015-12-16 南京理工大学 基于融合型自适应鲁棒的电液负载模拟器力控制方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110262226B (zh) * 2019-03-29 2022-05-06 浙江理工大学 一种欠驱动系统控制方法、装置、设备及介质
CN110262226A (zh) * 2019-03-29 2019-09-20 浙江理工大学 一种欠驱动系统控制方法、装置、设备及介质
CN110107563A (zh) * 2019-05-27 2019-08-09 电子科技大学 存在负载干扰情况下多电液伺服执行器分布协同控制方法
CN110173487A (zh) * 2019-05-27 2019-08-27 电子科技大学 一种多电液伺服执行器在切换网络下的领导同步控制方法
CN110081046A (zh) * 2019-05-27 2019-08-02 电子科技大学 一种基于反步控制的多电液伺服执行器跟踪同步控制方法
CN111258220A (zh) * 2020-01-21 2020-06-09 重庆邮电大学 一种基于干扰观测器的柔性机械臂级联控制方法及系统
CN111258220B (zh) * 2020-01-21 2022-12-06 重庆邮电大学 一种基于干扰观测器的柔性机械臂级联控制方法及系统
CN112049837A (zh) * 2020-08-05 2020-12-08 东北大学 基于参数化建模技术的液压装置数字化样机设计方法
CN111775142A (zh) * 2020-08-12 2020-10-16 电子科技大学 一种液压机械臂的模型辨识与自适应控制方法
CN112068435A (zh) * 2020-09-09 2020-12-11 北京航空航天大学 基于干扰观测器的康复机械装置迭代学习控制方法及系统
CN112180742A (zh) * 2020-11-09 2021-01-05 山东大学 基于分数阶扰动观测器的海上起重机稳定控制方法及系统
CN112180742B (zh) * 2020-11-09 2021-08-03 山东大学 基于分数阶扰动观测器的海上起重机稳定控制方法及系统
CN114801301A (zh) * 2022-04-28 2022-07-29 重庆智能机器人研究院 伺服电缸压力机的控制方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
CN106402089B (zh) 2017-11-17

Similar Documents

Publication Publication Date Title
CN106402089B (zh) 一种基于耦合干扰观测器的级联电液伺服系统控制方法及系统
CN106438593B (zh) 一种存在参数不确定性和负载干扰的电液伺服控制方法及机械臂
Ye et al. Position control of nonlinear hydraulic system using an improved PSO based PID controller
Feng et al. Identification and compensation of non-linear friction for a electro-hydraulic system
CN103983453B (zh) 一种航空发动机的执行机构和传感器故障诊断的区分方法
CN110928182B (zh) 基于状态估计的液压伺服系统鲁棒自适应重复控制方法
Zhou et al. Adaptive robust control design for underwater multi-dof hydraulic manipulator
CN103406909B (zh) 一种机械臂系统的跟踪控制设备及方法
CN103674552A (zh) 航空发动机试车台推力加载/校准系统及自动控制方法
CN105538310B (zh) 一种基于衰减记忆滤波的电液伺服控制方法及2‑dof机械臂
CN110456646B (zh) 基于ekf的锚杆钻机钻臂摆角自适应滑模控制系统设计方法
CN110081046B (zh) 一种基于反步控制的多电液伺服执行器跟踪同步控制方法
CN105563489A (zh) 基于非线性自抗扰控制技术的柔性机械臂控制方法
CN105068564A (zh) 一种压电陶瓷作动器的位移控制方法
CN102501251A (zh) 带有动态摩擦补偿的机械臂关节定位控制方法
CN112000009B (zh) 一种基于状态与扰动估计的物料转送器强化学习控制方法
CN104698844A (zh) 液压位置伺服系统的不确定性补偿的滑模控制方法
Guenther et al. Cascade controlled pneumatic positioning system with LuGre model based friction compensation
CN111648758B (zh) 一种水井钻机推进装置无模型自适应控制方法及系统
CN110107563B (zh) 存在负载干扰情况下多电液伺服执行器分布协同控制方法
CN105425610A (zh) 一种作业型水下机器人液压推进器的推力控制仿真系统
CN111796511A (zh) 液压执行器的小波神经网络pid在线控制方法及系统
CN111775142B (zh) 一种液压机械臂的模型辨识与自适应控制方法
CN111086646B (zh) 一种飞机襟翼操纵系统的状态监测方法及系统
CN109814393B (zh) 一种悬吊漂浮物随动系统的降维观测器和控制器设计方法

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: 20171117

Termination date: 20201024

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