CN113220007B - 执行机构故障的挠性航天器有限时间姿态协同控制方法 - Google Patents
执行机构故障的挠性航天器有限时间姿态协同控制方法 Download PDFInfo
- Publication number
- CN113220007B CN113220007B CN202110528568.0A CN202110528568A CN113220007B CN 113220007 B CN113220007 B CN 113220007B CN 202110528568 A CN202110528568 A CN 202110528568A CN 113220007 B CN113220007 B CN 113220007B
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- attitude
- flexible
- matrix
- formula
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 230000007246 mechanism Effects 0.000 title claims abstract description 14
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 42
- 239000011159 matrix material Substances 0.000 claims abstract description 40
- 238000004891 communication Methods 0.000 claims abstract description 9
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 14
- 238000013461 design Methods 0.000 claims description 11
- 150000001875 compounds Chemical class 0.000 claims description 8
- 238000013528 artificial neural network Methods 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000011217 control strategy Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 13
- 238000004088 simulation Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
- G05D1/0816—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
- G05D1/0833—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using limited authority control
Landscapes
- Engineering & Computer Science (AREA)
- Computer Security & Cryptography (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供一种执行机构故障的挠性航天器有限时间姿态协同控制方法,步骤一:基于旋转矩阵建立无退绕的挠性航天器姿态运动数学模型;步骤二:获取挠性航天器编队中各成员的姿态信息及参考轨迹,建立相对姿态动力学方程;步骤三:基于无向图描述编队系统各成员之间的通信拓扑结构,并结合集中误差设计快速终端滑模变量;步骤四:考虑到未知的外界扰动、难以观测挠性动态和时变的惯性参数,设计自适应容错控制器。本发明基于自适应控制对其不确定的动态进行实时估计和补偿,显著增强了航天器姿态协同容错控制策略的实用性。本发明尤其是最小学习参数算法的应用降低了控制系统对计算资源的需求,从而达到减少硬件负载的目的,具有很强的工程意义。
Description
技术领域
本发明涉及航天器姿态控制技术领域,具体为执行机构故障的挠性航天器有限时间姿态协同控制方法。
背景技术
随着太空事业的不断发展,人类对航天器的数量、质量和多样性提出了更高的要求。与传统的大型单体卫星系统相比,小型航天器编队在灵活性、容错性和成本方面拥有无可比拟的优势,可以更好地完成包括地面监测、在轨服务及合成孔径成像在内的各种太空任务。作为其中的关键技术之一,研究航天器编队姿态协同控制问题无疑具有极高的现实意义。
值得注意的是,考虑到外界环境的多变性和太空任务的复杂性,高性能航天器姿态协同控制器的设计仍然面临着诸多挑战。首先是由于环境扰动及参数摄动引起的系统不确定性问题,主要包括未知的干扰力矩和时变的惯性参数;随后,航天器自身附带的挠性结构也会对其姿态控制产生不容忽视的影响,例如太阳能帆板、太空机械臂、通讯天线等;最后,在轨航天器系统在漫长的服役周期中不可避免会出现各种故障现象,其中最常见的类型即为执行机构故障,包括部分失效故障及随机漂移故障等。因此,目前迫切地需要一种针对挠性航天器编队的姿态协同控制算法,以保证编队系统在外界扰动、参数摄动及执行机构故障条件下的鲁棒性和可靠性。
发明内容
本发明的目的在于提供解决挠性航天器姿态协同控制问题的执行机构故障的挠性航天器有限时间姿态协同控制方法。
本发明的目的是这样实现的:
执行机构故障的挠性航天器有限时间姿态协同控制方法,包括如下步骤:
步骤一:基于旋转矩阵建立无退绕的挠性航天器姿态运动数学模型;
步骤二:获取挠性航天器编队中各成员的姿态信息及参考轨迹,建立相对姿态动力学方程;
步骤三:基于无向图描述编队系统各成员之间的通信拓扑结构,并结合集中误差设计快速终端滑模变量;
步骤四:考虑到未知的外界扰动、难以观测的挠性动态和时变的惯性参数,设计自适应容错控制器。
本发明还包括这样一些结构特征:
1.所述步骤一具体为:
考虑航天器编队系统中存在n个成员,则针对第i个航天器可建立如下姿态动力学模型:
式中,Ri∈SO(3)表示该航天器的旋转矩阵,为其姿态旋转角速率,正定矩阵为惯量矩阵,是三阶模态下挠性附件与中心刚体的耦合矩阵,为挠性附件的三阶模态坐标,和分别为控制输入力矩和外界干扰力矩,和分别为阻尼矩阵及刚度矩阵。另外,需要引入运算符号×和执行机构故障函数fi(ui):
式中,Ki=diag(Ki1,Ki2,Ki3)表示执行机构乘性故障系数,且Kim(m=1,2,3)满足0<θi≤Kim≤1,θi>0表示能够维持执行机构正常运行的最低限度;表示执行机构的加性故障,且满足其中ψi为大于零的未知常数。
2.所述步骤二具体为:
进一步引入姿态误差向量ei:
基于上述阐述,挠性航天器编队姿态协同控制系统的相对姿态动力学模型可以表达为:
式中,矩阵函数Ei,Fi以及Gi分别为:
为了简化该编队系统模型,引入以下向量和矩阵:
E=diag(E1,…,En),J=diag(J1,…,Jn)
则挠性航天器编队姿态协同控制系统的数学模型可以被表述为:
3.所述步骤三具体为:
设计快速终端滑模变量:
式中,ki(i=1,2,3)均为大于零的设计参数,以A表示航天器编队通信拓扑无向图的加权邻接矩阵,aij为A的元素,若航天器i可以获取航天器j的状态信息,有aij>0;否则aij=0。在无向图中恒有aij=aji。
L=D-A
4.所述步骤四具体为:
航天器挠性附件的动态可以做如下处理:
另一方面,考虑到难以获知航天器精确的惯性参数,采用径向基神经网络对相关的不确定性进行逼近和补偿,其具体方法为:
式中,Wi=[wi1,wi2,...,wip]T为权值矩阵,基函数hi(x)=[hi1(x),...,hip(x)]T选择高斯函数,x=[x1,x2,...,xm]为输入向量,即有:
式中cij∈Rm为中心向量,bij∈R为函数宽度。
式中ki(i=4,5,6,7,8,9),μi(i=1,2,3),α,β,ξ,σ均为大于零的设计参数,且参数γ满足0<γ<1。
与现有技术相比,本发明的有益效果是:
1.与大多数发明中将航天器进行纯刚体动力学建模的方法不同,本发明考虑到了航天器挠性附件对姿态控制的不利影响,并基于自适应控制对其不确定的动态进行实时估计和补偿,显著增强了航天器姿态协同容错控制策略的实用性。
2.本发明针对太空任务中航天器的惯性参数会不断发生微小变化的现实问题,利用径向基神经网络对相关的未知动态实现了精确逼近,尤其是最小学习参数算法的应用有效降低了控制系统对计算资源的需求,从而达到减少硬件负载的目的,具有很强的工程意义。
附图说明
图1为本发明中挠性航天器有限时间姿态协同控制方法的流程图;
图2为挠性航天器编队中各成员的相对姿态轨迹;
图3为挠性航天器编队中各成员的相对角速度轨迹;
图4为挠性航天器编队中各成员的滑模变量曲线;
图5为挠性航天器编队中各成员的控制力矩曲线。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
如图1所示,本发明设计的执行机构故障的挠性航天器有限时间姿态协同控制方案,包括以下步骤:
第一步,基于旋转矩阵建立无退绕的挠性航天器姿态运动数学模型。考虑航天器编队系统中存在n个成员,则针对第i个航天器可建立如下姿态动力学模型:
式中,Ri∈SO(3)表示该航天器的旋转矩阵,为其姿态旋转角速率,正定矩阵为惯量矩阵,是三阶模态下挠性附件与中心刚体的耦合矩阵,为挠性附件的三阶模态坐标,和分别为控制输入力矩和外界干扰力矩,和分别为阻尼矩阵及刚度矩阵。另外,需要引入运算符号×和执行机构故障函数fi(ui):
式中,Ki=diag(Ki1,Ki2,Ki3)表示执行机构乘性故障系数,且Kim(m=1,2,3)满足0<θi≤Kim≤1,θi>0表示能够维持执行机构正常运行的最低限度;表示执行机构的加性故障,且满足其中ψi为大于零的未知常数。
进一步引入姿态误差向量ei:
基于上述阐述,挠性航天器编队姿态协同控制系统的相对姿态动力学模型可以表达为:
式中,矩阵函数Ei,Fi以及Gi分别为:
为了简化该编队系统模型,引入以下向量和矩阵:
E=diag(E1,…,En),J=diag(J1,…,Jn)
则挠性航天器编队姿态协同控制系统的数学模型可以被表述为:
为方便后续的控制器设计,针对本例中的挠性航天器编队系统,给出如下假设:
假设二:作用于挠性航天器的外界干扰di是有界的,即满足||di||≤Di,其中,Di为大于零的未知常数。
第三步,基于无向图描述编队系统各成员之间的通信拓扑结构,并结合集中误差设计快速终端滑模变量。其中,滑模变量被设计为:
式中,ki(i=1,2,3)均为大于零的设计参数,以A表示航天器编队通信拓扑无向图的加权邻接矩阵,aij为A的元素,若航天器i可以获取航天器j的状态信息,有aij>0;否则aij=0。在无向图中恒有aij=aji。
L=D-A
第四步,考虑到未知的外界扰动、难以观测的挠性动态和时变的惯性参数,设计自适应容错控制器。
首先,航天器挠性附件的动态可以做如下处理:
另一方面,考虑到难以获知航天器精确的惯性参数,采用径向基神经网络对相关的不确定性进行逼近和补偿,其具体方法为:
式中,Wi=[wi1,wi2,...,wip]T为权值矩阵,基函数hi(x)=[hi1(x),...,hip(x)]T选择高斯函数,x=[x1,x2,...,xm]为输入向量,即有:
式中cij∈Rm为中心向量,bij∈R为函数宽度。
式中ki(i=4,5,6,7,8,9),μi(i=1,2,3),α,β,ξ,σ均为大于零的设计参数,且参数γ满足0<γ<1。
接下来将通过仿真实例对上述控制器的性能进行展示和验证。针对由四个挠性航天器组成的编队系统,利用无向图描述其通信拓扑结构,首先给出加权邻接矩阵为:
四个航天器的惯性参数为:
其初始的姿态参数及外界扰动被设定为:
ω1(0)=ω2(0)=ω3(0)=ω4(0)=0,di=0.1×[sin(0.1t),cos(0.2t),sin(0.2t)]TN·m.
ξi,1=0.056,ξi,2=0.086,ξi,3=0.08,ξi,4=0.025
ιi,1=1.0973,ιi,2=1.2761,ιi,3=1.6538,ιi,4=2.2893
航天器编队的参考姿态轨迹设定为:
Rd(0)=I,ωd=0.15×[sin(t/50),-cos(t/50),-sin(t/60)]Trad/s
假设执行机构的故障参数为:
本例中将控制器的各参数取值为:
α=0.05;β=0.05,ξ=0.05,k1=0.5,k2=0.5,k3=0.02,k4=20,k5=20,k6=0.01
k7=0.01,k8=0.01,k9=0.01,γ=0.9,σ=0.05.
仿真结果如图2-5所示,由此可知本发明提出的执行机构故障的挠性航天器姿态协同控制算法具备有限时间稳定性能。
Claims (2)
1.执行机构故障的挠性航天器有限时间姿态协同控制方法,其特征是,包括如下步骤:
步骤一:基于旋转矩阵建立无退绕的挠性航天器姿态运动数学模型;
考虑航天器编队系统中存在n个成员,则针对第i个航天器可建立如下姿态动力学模型:
式中,ωi为第i个航天器的姿态旋转角速度,Ji为第i个航天器的惯量矩阵,fi(ui)为第i个航天器的执行机构故障函数,n为挠性航天器编队中的成员数量,ηi为第i个航天器的挠性附件的模态坐标,δi为第i个航天器挠性附件与中心刚体的耦合矩阵,di为第i个航天器的外界干扰力矩,ui为第i个航天器的控制输入力矩,Ri∈SO(3)表示该航天器的旋转矩阵,为其姿态旋转角速率,正定矩阵为惯量矩阵,是三阶模态下挠性附件与中心刚体的耦合矩阵,为挠性附件的三阶模态坐标,和分别为控制输入力矩和外界干扰力矩,和分别为阻尼矩阵及刚度矩阵;另外,需要引入运算符号×和执行机构故障函数fi(ui):
式中,Ki=diag(Ki1,Ki2,Ki3)表示执行机构乘性故障系数,且Kim(m=1,2,3)满足0<θi≤Kim≤1,θi>0表示能够维持执行机构正常运行的最低限度;表示执行机构的加性故障,且满足其中ψi为大于零的未知常数;
步骤二:获取挠性航天器编队中各成员的姿态信息及参考轨迹,建立相对姿态动力学方程;
进一步引入姿态误差向量ei:
基于上述阐述,挠性航天器编队姿态协同控制系统的相对姿态动力学模型可以表达为:
式中,矩阵函数Ei,Fi以及Gi分别为:
为了简化该编队系统模型,引入以下向量和矩阵:
E=diag(E1,…,En),J=diag(J1,…,Jn)
则挠性航天器编队姿态协同控制系统的数学模型可以被表述为:
步骤三:基于无向图描述编队系统各成员之间的通信拓扑结构,并结合集中误差设计快速终端滑模变量;
式中,ej为第j(j=1,2,...,n)个航天器的姿态误差向量,ki(i=1,2,3)均为大于零的设计参数,以A表示航天器编队通信拓扑无向图的加权邻接矩阵,aij为A的元素,若航天器i可以获取航天器j的状态信息,有aij>0;否则aij=0;在无向图中恒有aij=aji;
为方便后续的设计和推导,给出如下定义:
L=D-A
步骤四:考虑到未知的外界扰动、难以观测的挠性动态和时变的惯性参数,设计自适应容错控制器。
2.根据权利要求1所述的执行机构故障的挠性航天器有限时间姿态协同控制方法,其特征是,所述步骤四具体为:
航天器挠性附件的动态可以做如下处理:
另一方面,考虑到难以获知航天器精确的惯性参数,采用径向基神经网络对相关的不确定性进行逼近和补偿,其具体方法为:
式中,Wi=[wi1,wi2,...,wip]T为权值矩阵,基函数hi(x)=[hi1(x),...,hip(x)]T选择高斯函数,x=[x1,x2,...,xm]为输入向量,即有:
式中,cij∈Rm为中心向量,bij∈R为函数宽度;
式中,ki(i=4,5,6,7,8,9),μi(i=1,2,3),α,β,ξ,σ均为大于零的设计参数,且参数γ满足0<γ<1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110528568.0A CN113220007B (zh) | 2021-05-14 | 2021-05-14 | 执行机构故障的挠性航天器有限时间姿态协同控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110528568.0A CN113220007B (zh) | 2021-05-14 | 2021-05-14 | 执行机构故障的挠性航天器有限时间姿态协同控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113220007A CN113220007A (zh) | 2021-08-06 |
CN113220007B true CN113220007B (zh) | 2022-11-18 |
Family
ID=77092006
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110528568.0A Active CN113220007B (zh) | 2021-05-14 | 2021-05-14 | 执行机构故障的挠性航天器有限时间姿态协同控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113220007B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114489096B (zh) * | 2021-12-14 | 2023-11-10 | 中国运载火箭技术研究院 | 基于多机械臂驱动航天器位姿一体化鲁棒动力学控制方法 |
CN115014363B (zh) * | 2022-06-08 | 2023-01-24 | 哈尔滨工业大学(深圳) | 一种基于旋转矩阵拓扑结构的姿态路径规划方法、系统、设备和介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109459931A (zh) * | 2018-05-09 | 2019-03-12 | 南京理工大学 | 一种航天器编队有限时间姿态容错控制方法 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2447792A1 (en) * | 2005-09-19 | 2012-05-02 | Cleveland State University | Controllers, observer, and applications thereof |
CN105204495A (zh) * | 2015-09-24 | 2015-12-30 | 哈尔滨工程大学 | 一种基于虚拟执行器的气垫船执行器故障处理方法 |
CN105843240B (zh) * | 2016-04-08 | 2019-01-22 | 北京航空航天大学 | 一种考虑执行器故障的航天器姿态积分滑模容错控制方法 |
CN106406086B (zh) * | 2016-05-26 | 2019-05-07 | 北京航空航天大学 | 一种基于滑模干扰观测器的大挠性航天器干扰补偿方法 |
CN107577145B (zh) * | 2017-08-25 | 2020-06-09 | 湘潭大学 | 编队飞行航天器反步滑模控制方法 |
CN108153322B (zh) * | 2017-12-06 | 2019-03-29 | 北京航空航天大学 | 一种考虑时变的转动惯量的航天器姿态跟踪自适应容错控制方法 |
CN108333949B (zh) * | 2018-03-09 | 2020-10-02 | 南京航空航天大学 | 针对多智能体系统执行器失效故障的滑模容错控制方法 |
CN108919651B (zh) * | 2018-09-26 | 2020-09-01 | 北京航空航天大学 | 一种考虑执行器输入约束的高超声速飞行器自适应容错控制方法 |
CN109213184B (zh) * | 2018-11-06 | 2021-06-08 | 哈尔滨工业大学(深圳) | 挠性航天器的有限时间多模态滑模姿态控制算法 |
CN109683628B (zh) * | 2018-12-26 | 2022-01-25 | 哈尔滨工程大学 | 一种基于有限时间分布式速度观测器的航天器相对位置控制方法 |
CN111258221B (zh) * | 2020-01-21 | 2022-09-02 | 中国西安卫星测控中心 | 一种基于自适应滑模理论的航天器容错控制方法 |
CN112357119A (zh) * | 2020-07-31 | 2021-02-12 | 盐城工学院 | 一种输入受限的有限时间姿态协同跟踪容错控制方法 |
-
2021
- 2021-05-14 CN CN202110528568.0A patent/CN113220007B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109459931A (zh) * | 2018-05-09 | 2019-03-12 | 南京理工大学 | 一种航天器编队有限时间姿态容错控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113220007A (zh) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Van | An enhanced tracking control of marine surface vessels based on adaptive integral sliding mode control and disturbance observer | |
CN107703742B (zh) | 一种柔性航天器传感器故障调节方法 | |
Chen et al. | Consensus control for multiple AUVs under imperfect information caused by communication faults | |
Chen et al. | Guaranteed transient performance based control with input saturation for near space vehicles | |
Zhang et al. | Hybrid fuzzy adaptive fault-tolerant control for a class of uncertain nonlinear systems with unmeasured states | |
CN113220007B (zh) | 执行机构故障的挠性航天器有限时间姿态协同控制方法 | |
Lee | Geometric controls for a tethered quadrotor UAV | |
CN110488603B (zh) | 一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法 | |
CN113589689B (zh) | 一种基于多参数自适应神经网络的滑模控制器设计方法 | |
CN110543183B (zh) | 一种考虑执行器受限问题的刚体飞行器固定时间姿态跟踪控制方法 | |
Wang et al. | Fault‐tolerant spacecraft attitude control under actuator saturation and without angular velocity | |
Erdong et al. | Robust attitude synchronisation controllers design for spacecraft formation | |
CN109188910B (zh) | 一种刚性飞行器的自适应神经网络容错跟踪控制方法 | |
CN111857181A (zh) | 分布式航天器编队自适应有限时间姿态协同跟踪控制方法 | |
Chen et al. | Sliding-mode-based attitude tracking control of spacecraft under reaction wheel uncertainties | |
Liu et al. | Adaptive fault‐tolerant formation control for quadrotors with actuator faults | |
Sun | Constrained adaptive fault-tolerant attitude tracking control of rigid spacecraft | |
Huang et al. | Distributed flocking control of quad-rotor UAVs with obstacle avoidance under the parallel-triggered scheme | |
CN114879515A (zh) | 基于学习神经网络的航天器姿态重构容错控制方法 | |
CN108958043B (zh) | 一种刚性飞行器的固定时间自适应姿态容错控制方法 | |
Mustafa et al. | Adaptive backstepping sliding mode control based on nonlinear disturbance observer for trajectory tracking of robotic manipulator | |
CN110488855B (zh) | 一种基于神经网络估计的刚性飞行器自适应固定时间姿态容错控制方法 | |
Ye et al. | Adaptive switching control for hypersonic vehicle with uncertain control direction | |
Guiming et al. | Coordinated multiple spacecraft attitude control with communication time delays and uncertainties | |
Yin et al. | Fixed-time attitude cooperative fault-tolerant control with prescribed performance for heterogeneous multiple satellite |
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 |