CN109213184A - 挠性航天器的有限时间多模态滑模姿态控制算法 - Google Patents
挠性航天器的有限时间多模态滑模姿态控制算法 Download PDFInfo
- Publication number
- CN109213184A CN109213184A CN201811312299.9A CN201811312299A CN109213184A CN 109213184 A CN109213184 A CN 109213184A CN 201811312299 A CN201811312299 A CN 201811312299A CN 109213184 A CN109213184 A CN 109213184A
- Authority
- CN
- China
- Prior art keywords
- flexible
- finite time
- mode
- modal
- sliding
- 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
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive 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
- G05B13/042—Adaptive 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 in which a parameter or coefficient is automatically adjusted to optimise the performance
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)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种挠性航天器的有限时间多模态滑模姿态控制算法,利用姿态四元数和欧拉轴/角表示方法建立挠性航天器运动学方程,采用混合坐标法对中心刚体带有挠性附件、外部干扰、转动惯量不确定性的挠性航天器建立动力学方程;采用多模态控制思想,基于Lyapunov有限时间稳定定理设计以下两种滑模控制律:针对挠性模态可测量的情况,设计一种多模态有限时间滑模控制律;针对挠性模态不可测量的情况设计一种基于动态观测器的多模态有限时间滑模控制律。本发明的有益效果是:采用本发明设计的姿态控制算法使挠性航天器较快实现姿态稳定,对外界干扰和转动惯量不确定性具有较鲁棒性。
Description
技术领域
本发明涉及挠性航天器控制技术领域,尤其涉及一种挠性航天器的有限时间多模态滑模姿态控制算法。
背景技术
传统的挠性航天器姿态滑模控制算法中,不考虑挠性航天器的转动惯量存在不确定性和外界干扰,并且传统的滑模控制算法只保证系统状态在单一滑模面滑动,并且无法保证系统有限时间稳定。
发明内容
为了解决现有技术中的问题,本发明提供了一种挠性航天器的有限时间多模态滑模姿态控制算法。本发明具体通过如下技术方案实现:
一种挠性航天器的有限时间多模态滑模姿态控制算法,利用姿态四元数和欧拉轴/角表示方法建立挠性航天器运动学方程,采用混合坐标法对中心刚体带有挠性附件、外部干扰、转动惯量不确定性的挠性航天器建立动力学方程;采用多模态控制思想,基于Lyapunov有限时间稳定定理设计以下两种滑模控制律:针对挠性模态可测量的情况,设计一种多模态有限时间滑模控制律;针对挠性模态不可测量的情况设计一种基于动态观测器的多模态有限时间滑模控制律。
作为本发明的进一步改进,以姿态四元数和欧拉轴/角为基础的挠性航天器运动学方程和动力学方程如下:
其中,q0,qv分别为姿态四元数的标量部分与向量部分;ω是航天器的姿态角;δ为挠性航天器的挠性部分与刚体主体之间的耦合矩阵;C,K分别为阻尼矩阵和刚度矩阵,
C=diag{2ξ1ωn1,2ξ2ωn2,...,2ξNωnN}
考虑N个弹性模态,其对应的自然角频率为ωni,i=1,2,...,N,对应的阻尼为ξi,i=1,2,...,N;Jmb为刚体部分的转动惯量,是Jmb的期望值,为转动惯量不确定系,u表示控制力矩,d表示有界外部干扰力矩,卫星带挠性附件,并且含有有界外部干扰力矩和转动惯量不确定性。
作为本发明的进一步改进,针对挠性模态可测的情形下,设计如下多模态滑模面:
设计的多模态有限时间挠性卫星姿态滑模控制律如下:
作为本发明的进一步改进,针对挠性模态不可测的情形下,设计如下多模态滑模面:
设计的基于动态观测器的有限时间滑模控制律如下:
正定对称矩阵P满足下述Lyapunov方程:
本发明的有益效果是:本发明解决了挠性航天器在飞行过程中存在外界干扰和转动惯量不确定性时的姿态控制问题,具有有限时间稳定优势,保证了挠性航天器姿态稳定的较好鲁棒性控制;采用本发明设计的姿态控制器保证航天器姿态的良好稳定性和模态振动的衰减,当卫星系统存在外界干扰和转动惯量不确定性时,航天器的姿态能在有限时间内较快达到稳定,能够有效消除外界干扰和转动惯量不确定性带来的影响。。
附图说明
图1是挠性模态可测的有限时间多模态滑模姿态控制系统的框图;
图2是挠性模态不可测的有限时间多模态滑模姿态控制系统的框图。。
具体实施方式
下面结合附图说明及具体实施方式对本发明进一步说明。
如图1至图2所示,本发明的挠性航天器的有限时间多模态滑模姿态控制算法,包括:
首先,以姿态四元数和欧拉轴/角为基础对存在有界外部干扰力矩和转动惯量不确定性的挠性航天器建立如下运动学方程和动力学方程:
其中,Jmb是主体惯量矩阵,是Jmb的期望值,为Jmb的逆,q0,qv分别为姿态四元数的标量部分与向量部分,e为欧拉轴,为旋转角,和分别为q0,qv,e和的一阶导数,I3表示3×3的单位阵,矩阵qv的乘积矩阵定义为表示为q1,q2,q3分别为qv的三个分量,即同理,e×为e的乘积矩阵,eT是e的转置操作,ω是航天器的姿态角,δ为挠性航天器的挠性部分与刚体主体之间的耦合矩阵,ψ为挠性附件的总速度,η为振动模态坐标向量,u为作用在航天器上的外部控制力矩,d是有界的外部干扰力矩,可为常值干扰或者正弦干扰,是转动惯量不确定项。其特征在于:航天器带有挠性附件,并且在运行过程中受各种干扰力矩的影响且存在惯性不确定性。
考虑以下两种情况:
(1)假设卫星的状态η和ψ完全可测量,针对挠性航天器(1)-(3)设计多模态有限时间滑模控制律。
S101、选取多模态滑模面
选取如下的多模态滑模切换函数:
其中,k1,k2,k3,α,β,γ均为正标量,且γ满足1/2<γ<1,e为欧拉角,||qv||为qv的范数。
由于角速度在切换点处应该是连续的,控制参数之间应满足如下关系:
k1=αk2,k2=βγ-1k3
并证明系统状态在所选各段滑模面进行有限时间滑动,闭环系统在最后一阶段滑模面上稳定运行,保证该系统s=0时有限时间稳定性得到保证。
证明:选取Lyapunov函数Vq为:
S102、设计控制律
设计以下多模态滑模控制律:
其中,k为常量,p为正标量且满足1>p>0;se是向量s的单位方向向量,满足se=s/||s||;λ是正标量,满足为矩阵的最大特征值,为外界干扰力矩d的上界。
在滑模控制律(5)作用下,从任意初始状态出发挠性航天器系统能在有限时间内到达滑模面s=0,即ω=-k3||qv||γe,并在滑模面上保持稳定运行。
以下根据多模态控制思想设计多模态控制律u,选取以下的Lyapunov函数:
令其中,c,a为常数,满足c>0,1>a>0,得到控制律u。
(2)在实际应用中模态η和ψ难以测量时,针对挠性航天器系统(1)-(3)设计基于动态观测器的滑模控制律。
采用下式的动态观测器:
定义如下模态观测误差:
其中,分别为η和ψ的估计值,I为单位矩阵,p-1是矩阵P的逆,正定对称矩阵P满足下述Lyapunov方程
其中,Q为正定矩阵。
S201、选取多模态滑模面
选取如下的多模态滑模切换函数:
并证明系统状态在所选各段滑模面进行有限时间滑动,闭环系统在最后一阶段滑模面上稳定运行,保证该系统s=0时有限时间稳定性得到保证。
证明:选取Lyapunov函数为:
S202、设计控制律
设计以下多模态滑模控制律:
在滑模控制律作用下,从任意初始状态出发挠性航天器系统能在有限时间内到达滑模面s=0,即ω=-k3||qv||γe,并在滑模面上保持稳定运行。
以下根据多模态控制思想设计多模态控制律u,选取以下的Lyapunov函数:
最后,通过计算证明其中,c>0,1>a>0,根据有限时间稳定定理可得在该控制律下可保证系统的有限时间稳定性。
本发明提供的一种挠性航天器的有限时间多模态滑模姿态控制算法,针对存在外界干扰和转动惯量不确定性的挠性航天器姿态控制问题,设计了一种有限时间多模态滑模姿态控制算法。该发明算法的目的在于解决挠性航天器在执行任务过程中存在有界干扰和转动惯量不确定性时的姿态控制及挠性附件的振动抑制问题。该发明利用姿态四元数和欧拉轴/角表示方法建立挠性航天器运动学方程,采用混合坐标法对中心刚体带有挠性附件、外部干扰、转动惯量不确定性的挠性航天器建立动力学方程;采用多模态控制思想,基于Lyapunov有限时间稳定定理设计一种多模态有限时间滑模控制律。并构造挠性模态观测器测量挠性状态变量,设计带挠性模态观测器的多模态有限时间滑模控制律。最后,运用MATLAB中的Simulink模块验证设计的控制算法的有效性。
以上所述仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (4)
1.一种挠性航天器的有限时间多模态滑模姿态控制算法,其特征在于:所述算法首先利用姿态四元数和欧拉轴/角表示方法建立挠性航天器运动学方程,采用混合坐标法对中心刚体带有挠性附件、外部干扰、转动惯量不确定性的挠性航天器建立动力学方程;采用多模态控制思想,基于Lyapunov有限时间稳定定理设计以下两种滑模控制律:针对挠性模态可测量的情况,设计基于状态反馈的有限时间滑模控制律;针对挠性模态不可测量的情况设计基于动态观测器的有限时间滑模控制律。
2.根据权利要求1所述的挠性航天器的有限时间多模态滑模姿态控制算法,其特征在于,以姿态四元数和欧拉轴/角为基础的挠性卫星的运动学方程和动力学方程如下:
其中,q0,qv分别为姿态四元数的标量部分与向量部分;ω是航天器的姿态角;δ为挠性航天器的挠性部分与刚体主体之间的耦合矩阵;C,K分别为阻尼矩阵和刚度矩阵,
C=diag{2ξ1ωn1,2ξ2ωn2,...,2ξNωnN}
考虑N个弹性模态,其对应的自然角频率为ωni,i=1,2,...,N,对应的阻尼为ξi,i=1,2,...,N;Jmb为刚体部分的转动惯量,是Jmb的期望值,为转动惯量不确定项,u表示控制力矩,d表示有界外部干扰力矩,卫星带挠性附件,并且含有有界外部干扰力矩和转动惯量不确定性。
3.根据权利要求2所述的挠性航天器的有限时间多模态滑模姿态控制算法,其特征在于,针对挠性模态可测的情形下,设计如下多模态滑模面:
设计的所述多模态有限时间挠性卫星姿态滑模控制律如下:
4.根据权利要求2所述的挠性航天器的有限时间多模态滑模姿态控制算法,其特征在于,针对挠性模态不可测的情形下,设计如下多模态滑模面:
设计的基于动态观测器的有限时间滑模控制律如下:
其中,正定对称矩阵p满足下述Lyapunov方程:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811312299.9A CN109213184B (zh) | 2018-11-06 | 2018-11-06 | 挠性航天器的有限时间多模态滑模姿态控制算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811312299.9A CN109213184B (zh) | 2018-11-06 | 2018-11-06 | 挠性航天器的有限时间多模态滑模姿态控制算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109213184A true CN109213184A (zh) | 2019-01-15 |
CN109213184B CN109213184B (zh) | 2021-06-08 |
Family
ID=64995461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811312299.9A Expired - Fee Related CN109213184B (zh) | 2018-11-06 | 2018-11-06 | 挠性航天器的有限时间多模态滑模姿态控制算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109213184B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507892A (zh) * | 2019-01-22 | 2019-03-22 | 哈尔滨工业大学(深圳) | 挠性航天器的自适应滑模姿态稳定控制方法 |
CN110083171A (zh) * | 2019-04-30 | 2019-08-02 | 哈尔滨工业大学(深圳) | 挠性航天器的动态滑模姿态跟踪控制的方法及系统 |
CN111498147A (zh) * | 2020-04-03 | 2020-08-07 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 挠性航天器的有限时间分段滑模姿态跟踪控制算法 |
CN113220007A (zh) * | 2021-05-14 | 2021-08-06 | 哈尔滨工程大学 | 执行机构故障的挠性航天器有限时间姿态协同控制方法 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6904353B1 (en) * | 2003-12-18 | 2005-06-07 | Honeywell International, Inc. | Method and system for sliding mode control of a turbocharger |
US7035725B1 (en) * | 2001-12-17 | 2006-04-25 | Hydro-Aire, Inc. | Sliding integral proportional (SIP) controller for aircraft skid control |
CN103412491A (zh) * | 2013-08-27 | 2013-11-27 | 北京理工大学 | 一种挠性航天器特征轴姿态机动指数时变滑模控制方法 |
CN104614995A (zh) * | 2015-02-16 | 2015-05-13 | 天津大学 | 一种二阶系统有限时间滑模控制器的通用设计方法 |
CN104898683A (zh) * | 2015-05-20 | 2015-09-09 | 哈尔滨工业大学 | 一种挠性卫星神经网络反步滑模姿态控制方法 |
CN104950899A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种固定时间收敛的飞行器姿态控制方法 |
WO2016073252A1 (en) * | 2014-11-03 | 2016-05-12 | Echogen Power Systems, L.L.C. | Active thrust management of a turbopump within a supercritical working fluid circuit in a heat engine system |
CN106194903A (zh) * | 2016-09-28 | 2016-12-07 | 西安理工大学 | 一种气动位置伺服系统的分数阶滑模变结构控制方法 |
CN106406086A (zh) * | 2016-05-26 | 2017-02-15 | 北京航空航天大学 | 一种基于滑模干扰观测器的大挠性航天器干扰补偿方法 |
US20170153650A1 (en) * | 2015-11-30 | 2017-06-01 | Metal Industries Research & Development Centre | Multiple rotors aircraft and control method |
CN107065913A (zh) * | 2017-05-10 | 2017-08-18 | 哈尔滨工业大学深圳研究生院 | 挠性航天器的滑模变结构姿态控制算法 |
CN107943062A (zh) * | 2017-09-13 | 2018-04-20 | 哈尔滨工业大学深圳研究生院 | 带外部干扰力矩的挠性卫星姿态滑模控制方法 |
CA3040564A1 (en) * | 2016-10-18 | 2018-04-26 | Deakin University | Thrust vectored multicopters |
CN108490786A (zh) * | 2018-04-27 | 2018-09-04 | 哈尔滨工业大学 | 一种基于终端滑模的高超声速飞行器鲁棒跟踪控制方法 |
-
2018
- 2018-11-06 CN CN201811312299.9A patent/CN109213184B/zh not_active Expired - Fee Related
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7035725B1 (en) * | 2001-12-17 | 2006-04-25 | Hydro-Aire, Inc. | Sliding integral proportional (SIP) controller for aircraft skid control |
US6904353B1 (en) * | 2003-12-18 | 2005-06-07 | Honeywell International, Inc. | Method and system for sliding mode control of a turbocharger |
CN103412491A (zh) * | 2013-08-27 | 2013-11-27 | 北京理工大学 | 一种挠性航天器特征轴姿态机动指数时变滑模控制方法 |
WO2016073252A1 (en) * | 2014-11-03 | 2016-05-12 | Echogen Power Systems, L.L.C. | Active thrust management of a turbopump within a supercritical working fluid circuit in a heat engine system |
CN104614995A (zh) * | 2015-02-16 | 2015-05-13 | 天津大学 | 一种二阶系统有限时间滑模控制器的通用设计方法 |
CN104898683A (zh) * | 2015-05-20 | 2015-09-09 | 哈尔滨工业大学 | 一种挠性卫星神经网络反步滑模姿态控制方法 |
CN104950899A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种固定时间收敛的飞行器姿态控制方法 |
US20170153650A1 (en) * | 2015-11-30 | 2017-06-01 | Metal Industries Research & Development Centre | Multiple rotors aircraft and control method |
CN106406086A (zh) * | 2016-05-26 | 2017-02-15 | 北京航空航天大学 | 一种基于滑模干扰观测器的大挠性航天器干扰补偿方法 |
CN106194903A (zh) * | 2016-09-28 | 2016-12-07 | 西安理工大学 | 一种气动位置伺服系统的分数阶滑模变结构控制方法 |
CA3040564A1 (en) * | 2016-10-18 | 2018-04-26 | Deakin University | Thrust vectored multicopters |
CN107065913A (zh) * | 2017-05-10 | 2017-08-18 | 哈尔滨工业大学深圳研究生院 | 挠性航天器的滑模变结构姿态控制算法 |
CN107943062A (zh) * | 2017-09-13 | 2018-04-20 | 哈尔滨工业大学深圳研究生院 | 带外部干扰力矩的挠性卫星姿态滑模控制方法 |
CN108490786A (zh) * | 2018-04-27 | 2018-09-04 | 哈尔滨工业大学 | 一种基于终端滑模的高超声速飞行器鲁棒跟踪控制方法 |
Non-Patent Citations (5)
Title |
---|
ABRAHAM VILLANUEVA,等: "Multi-mode flight sliding mode control system for a quadrotor", 《2015 INTERNATIONAL CONFERENCE ON UNMANNED AIRCRAFT SYSTEMS (ICUAS)》 * |
ZHU QING HUA,等: "Attitude control without angular velocity measurement for flexible satellites", 《CHINESE JOURNAL OF AERONAUTICS》 * |
刘赛娜,等: "空间自动对接多模态滑模控制", 《宇航学报》 * |
朱庆华,等: "基于动态滑模控制的挠性航天器姿态控制", 《控制理论与应用》 * |
杨新岩,等: "基于分段趋近律的航天器对地凝视姿态滑模控制", 《科学技术与工程》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507892A (zh) * | 2019-01-22 | 2019-03-22 | 哈尔滨工业大学(深圳) | 挠性航天器的自适应滑模姿态稳定控制方法 |
CN110083171A (zh) * | 2019-04-30 | 2019-08-02 | 哈尔滨工业大学(深圳) | 挠性航天器的动态滑模姿态跟踪控制的方法及系统 |
CN111498147A (zh) * | 2020-04-03 | 2020-08-07 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 挠性航天器的有限时间分段滑模姿态跟踪控制算法 |
CN111498147B (zh) * | 2020-04-03 | 2021-09-21 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 挠性航天器的有限时间分段滑模姿态跟踪控制算法 |
CN113220007A (zh) * | 2021-05-14 | 2021-08-06 | 哈尔滨工程大学 | 执行机构故障的挠性航天器有限时间姿态协同控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109213184B (zh) | 2021-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109213184A (zh) | 挠性航天器的有限时间多模态滑模姿态控制算法 | |
CN104723340B (zh) | 基于连接和阻尼配置的柔性关节机械臂的阻抗控制方法 | |
CN111498147B (zh) | 挠性航天器的有限时间分段滑模姿态跟踪控制算法 | |
Hu et al. | Coordinated adaptive robust contouring control of an industrial biaxial precision gantry with cogging force compensations | |
CN108241292B (zh) | 一种基于扩张状态观测器的水下机器人滑模控制方法 | |
CN110543184B (zh) | 一种刚性飞行器的固定时间神经网络控制方法 | |
CN108663936B (zh) | 模型不确定航天器无退绕姿态跟踪有限时间控制方法 | |
CN104656447A (zh) | 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法 | |
CN110543183B (zh) | 一种考虑执行器受限问题的刚体飞行器固定时间姿态跟踪控制方法 | |
CN109634291B (zh) | 一种基于改进型障碍李雅普诺夫函数的刚性飞行器姿态约束跟踪控制方法 | |
CN113306747A (zh) | 基于so(3)群的挠性航天器姿态稳定控制方法和系统 | |
CN111258216A (zh) | 一种适用于四旋翼飞行器的滑模重复控制器 | |
CN107678281A (zh) | 基于修正型罗德里格参数的挠性航天器自适应姿态控制律 | |
Dirksz et al. | Adaptive tracking control of fully actuated port-Hamiltonian mechanical systems | |
CN105912007A (zh) | 空间机械臂抗干扰姿态稳定的微分几何非线性控制方法 | |
CN111506095A (zh) | 一种双刚体特征点间饱和固定时间相对位姿跟踪控制方法 | |
CN109683480A (zh) | 考虑执行器故障的非线性机械系统类固定时间控制方法 | |
CN110488854B (zh) | 一种基于神经网络估计的刚性飞行器固定时间姿态跟踪控制方法 | |
CN115933725B (zh) | 一种刚柔液耦合航天器高精度姿态控制方法 | |
Hu et al. | Gradient projection of weighted Jacobian matrix method for inverse kinematics of a space robot with a controlled-floating base | |
CN107943062A (zh) | 带外部干扰力矩的挠性卫星姿态滑模控制方法 | |
CN110488855B (zh) | 一种基于神经网络估计的刚性飞行器自适应固定时间姿态容错控制方法 | |
Zhang et al. | Finite-time attitude optimization maneuver control for coupled spacecraft under attitude measurement errors and actuator faults | |
Hu | Semi-globally input-to-state stable controller design for flexible spacecraft attitude stabilization under bounded disturbances | |
CN113485397B (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 | ||
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: 20210608 Termination date: 20211106 |