CN107844124B - 一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法 - Google Patents
一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法 Download PDFInfo
- Publication number
- CN107844124B CN107844124B CN201711245553.3A CN201711245553A CN107844124B CN 107844124 B CN107844124 B CN 107844124B CN 201711245553 A CN201711245553 A CN 201711245553A CN 107844124 B CN107844124 B CN 107844124B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- quadrotor
- axis
- gravity
- center
- 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 19
- 230000005484 gravity Effects 0.000 claims abstract description 28
- 238000013528 artificial neural network Methods 0.000 claims abstract description 12
- 230000001133 acceleration Effects 0.000 claims description 7
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 206010034719 Personality change Diseases 0.000 claims 2
- 230000007704 transition Effects 0.000 claims 1
- 238000013459 approach Methods 0.000 abstract description 3
- 230000006978 adaptation Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 5
- 230000002146 bilateral effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 238000007634 remodeling Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification 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/0825—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Feedback Control In General (AREA)
Abstract
本发明所述的四旋翼飞行器携带不平衡负载稳定飞行的控制方法,主要对携带不平衡负载的飞行器进行重新建模,串级PID控制其稳定飞行,利用神经网络自适应来逼近负载晃动,当负载的重心位置严重偏离出机体之外时仍然能够稳定飞行。
Description
技术领域:
本发明属于飞行控制技术领域。其突出特点是四旋翼飞行器在携带负载飞行的过程中重物突然偏离,使四旋翼飞行器的重心偏离出机身之外,破环其飞行的稳定性,但是在本文的改进的模型和算法下四旋翼飞行器依然能够稳定飞行。
背景技术:
随着微电子技术、传感器技术以及各种控制算法的进步,四旋翼飞行器成为人们研究的热点。它常被用来搜救、航拍、物流服务、军事侦察以及目标打击等诸多领域。2016年12月7日亚马逊在英国成功完成了首次无人机送货。2017年二月中旬顺丰快递公司也在开发无人机运送快递。国内外的很多学者研究了四旋翼飞行器稳定控制。多数针对挂在四旋翼质心位置下的摇摆不定的负载问题进行建模仿真研究。都未考虑到负载严重失衡的情况下飞行器的整体模型,针对这一问题,本发明将针对四旋翼飞行器在不平衡负载下稳定飞行进行研究,具有很大的实用价值。
发明内容:
本发明所要解决的技术问题是对携带不平衡负载的飞行器进行重新建模,串级PID控制其稳定飞行,利用神经网络自适应来逼近负载晃动,当负载的重心位置严重偏离出机体之外时仍然能够稳定飞行。
(1)根据起飞姿态信息,重新计算不平衡负载下飞行器重心位置
在不平衡负载下,以机体坐标系下,飞行器的重心位置由之前的O(0,0,0)移至G(xg,yg,0),重心位置通过从开机到离地这段时间的飞行参数来估计,因此,起飞时根据MPU6050传感器采集到的姿态信息推导出飞行器的姿态角和重心坐标之间存在如下的关系:
(2)力学平衡方程的建立
设四旋翼的总质量为M=m+m′,m′为负载的质量,四旋翼在空中悬停时四个螺旋桨的升力和与四旋翼的重力相等,F1~F4分别为四只螺旋桨产生的升力,可以写为公式(2),F1+F2+F3+F4=Mg (2)
每个电机电机到机身中心的距离均为l转动过程中每个螺旋桨产生的绕x轴和y轴的力矩平衡方程为:
悬停时还需要保证反扭矩为0,式中M1,M3,M2,M4为1-4号电机角速度为ω1,ω2,ω3,ω4产生的反转矩
M1+M3=M2+M4 (4)
结合(1)(2)(3)得到在不平衡负载下,飞行器在空中悬停的受力平衡方程组(5);
(3)角速度平衡方程的建立
ω1+ω3=ω2+ω4=ω (6)
为了同时保证式(5)和式(6)同时成立,可以对式(6)做如下修改:
(ω0+Δω13)+(ω0-Δω13)=(ω0+Δω24)+(ω0-Δω24) (7)
其中,ω0是飞行器在平衡负载下每只螺旋桨的平均转速,由式(5)和式(7),可以得出Δω13和Δω24分别为:
那么在不平衡负载下,四只螺旋桨的转速做下面的修正,达到预期位置和姿态;
4)不平衡负载下飞行器的运动学分析
将四旋翼飞行器看作一个在空间中运动的刚体在地面坐标系下的受力分析为
Fx、Fy和Fz分别是飞行器的总升力沿地面坐标系x轴、y轴和z轴方向上的分力Sx、Sy和Sz分别为飞行器沿x、y和z轴方向上的位移,其一阶导数表示速度,二阶导数表示加速度;
在不平衡负载下飞行器的动力学模型为:
θ,ψ的一阶导数表示角速度,二阶导数表示角加速度,Jx、Jy和Jz分别是飞行器三个轴的转动惯量;可以发现,当重心位置G(xg,yg,0)在机体坐标原点的时候,即当xg=0,yg=0时,该系统模型就和一般的平衡负载模型一致。
5)姿态解算
四旋翼飞行器飞行过程中期望的姿态角为0,应用欧拉角方法将四旋翼飞行器在空中飞行的姿态信息反馈到主控制器中,根据期望输出与实际输出的误差设计串级PID控制器控制,计算出控制量,转化为相应的控制信号经驱动电路后驱动四个电机工作,保持四旋翼飞行器稳定飞行。
2、一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法,其特征在于,所述串级 PID控制器考虑到四旋翼飞行器在工作过程中受到非线性的因素的影响,在串级PID的基础设计了神经网络自适应控制算法逼近系统的非线性部分,算法设计如下:
将公式(13)所描述的四旋翼飞行器运动方程写成如下形式:
具体推导步骤及证明过程如下
设计自适应律为:
设Lyapunov函数为
V=V1+V2
Υ为正常数,矩阵P为对称正定的且满足如下Lyapunov方程:
ATP+PA=-Q
附图说明
图1是不平衡负载下飞行器的受力图;
图2是挂重物重心位置变化图;
图3整个发明的流程图;
图4翻滚角仿真实验结果;
图5俯仰角仿真实验结果;
图6偏航角仿真实验结果;
图中曲线说明:本发明进行了仿真和实验验证,做了两组对比实验:
(1)在已有的串级PID为控制算法,分别控制常规的四旋翼模型(曲线1)和本发明推导的四旋翼模型(曲线2);
(2)以本发明所推导的模型为控制对象,分别用本发明的算法(曲线3)和已有发明的串级PID算法(曲线2)进行对比,发现本发明控制的四旋翼飞行器更稳定;
图7为实际飞行过程中负载在距中心10cm时翻滚角比较;
图8为实际飞行过程中负载在距中心10cm时俯仰角角比较;
图9为实际飞行过程中负载在距中心20cm时翻滚角比较;
图10为实际飞行过程中负载在距中心20cm时俯仰角比较;
图11为实际飞行过程中双边负载在距中心20cm时翻滚角比较;
图12为实际飞行过程中双边负载在距中心20cm时俯仰角比较;
具体实施方式
下面结合附图对本发明做详细的描述
本发明所述的四旋翼飞行器携带不平衡负载稳定飞行的控制方法,主要对携带不平衡负载的飞行器进行重新建模,串级PID控制其稳定飞行,利用神经网络自适应来逼近负载晃动,当负载的重心位置严重偏离出机体之外时仍然能够稳定飞行。
参阅图1,图2,所述的四旋翼飞行器携带不平衡负载稳定飞行的方法包括四旋翼飞行器的重新建模:
(1)根据起飞姿态信息,重新计算不平衡负载下飞行器重心位置
在不平衡负载下,以机体坐标系下,飞行器的重心位置由之前的O(0,0,0)移至G(xg,yg,0),重心位置通过从开机到离地这段时间的飞行参数来估计,因此,起飞时根据MPU6050传感器采集到的姿态信息推导出飞行器的姿态角和重心坐标之间存在如下的关系:
(2)力学平衡方程的建立
根据四旋翼飞行器在空间的悬停状态时,升力与重力大小相等,反扭矩和为0以及绕x 轴和y轴的力矩平衡方程,得到在不平衡负载下,飞行器在空中悬停的受力平衡方程组;
其中,M1、M2、M3和M4分别为四只螺旋桨对飞行器产生的反扭矩,四只电机到机身中心的距离均为l,F1~F4分别为四只螺旋桨产生的升力;
(3)角速度平衡方程的建立
ω1+ω3=ω2+ω4=ω (3)
为了同时保证式(2)和式(3)同时成立,可以对式(3)做如下修改:
(ω0+Δω13)+(ω0-Δω13)=(ω0+Δω24)+(ω0-Δω24) (4)
其中,ω0是飞行器在平衡负载下每只螺旋桨的平均转速,由式(2)和式(4),可以得出Δω13和Δω24分别为:
那么在不平衡负载下,四只螺旋桨的转速做下面的修正,达到预期位置和姿态;
4)不平衡负载下飞行器的运动学分析
在不平衡负载下飞行器的动力学模型为:
θ,ψ的一阶导数表示角速度,二阶导数表示角加速度,Jx、Jy和Jz分别是飞行器三个轴的转动惯量;Sx、Sy和Sz分别为飞行器沿x、y和z轴方向上的位移,其一阶导数表示速度,二阶导数表示加速度;可以发现,当重心位置G(xg,yg,0)在机体坐标原点的时候,即当xg=0,yg=0,时,该系统模型就和一般的平衡负载模型一致。
5)姿态解算
四旋翼飞行器飞行过程中期望的姿态角为0,应用欧拉角方法将四旋翼飞行器在空中飞行的姿态信息反馈到主控制器中,根据期望输出与实际输出的误差设计串级PID控制器控制,计算出控制量,转化为相应的控制信号经驱动电路后驱动四个电机工作,保持四旋翼飞行器稳定飞行。
2、参阅图3,所述的四旋翼飞行器携带不平衡负载稳定飞行的控制方法,其特征在于,所述串级PID控制器考虑到四旋翼飞行器在工作过程中受到非线性的因素的影响,在串级 PID的基础设计了神经网络自适应控制算法逼近系统的非线性部分,算法设计如下:
将公式(8)所描述的四旋翼飞行器方写成如下形式:
y=x
设计自适应律为:
设Lyapunov函数为
V=V1+V2
γ为正常数,矩阵P为对称正定的且满足如下Lyapunov方程:
ATP+PA=-Q
3.为了验证以上算法的正确性,选择恒拓科技的小型四旋翼飞行器为实验平台,其规格参数如下
表1四旋翼飞行器参数表
符号 | 数值 | 单位 |
M | 1.5 | kg |
G | 9.81 | m/s<sup>2</sup> |
l | 0.25 | m |
I<sub>xx</sub>,I<sub>yy</sub> | 0.0081 | Kg/m<sup>2</sup> |
I<sub>zz</sub> | 0.0142 | Kg/m<sup>2</sup> |
并做了几组对比实验负载质量分别为0.285kg和0.130kg。单边负载质量为0.285kg(占机身质量的1/5)。在飞行器的2号桨位置,分别将重物添加到距离飞行器中心单边10、20cm 处,以及两个重物同时悬挂在相邻双边20cm的位置,总重量为0.415kg(占机身质量的 28/100)。实际飞行结果曲线见附图7-12。
Claims (2)
1.一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法,其特征在于,包含以下步骤:
(1)根据起飞姿态信息,重新计算不平衡负载下飞行器重心位置
在不平衡负载下,以机体坐标系下,飞行器的重心位置由之前的0(0,0,0)移至G(xg,yg,0),其中,xg,yg分别为重心在x轴和y轴上的位移;重心位置通过从开机到离地这段时间的飞行参数来估计,因此,起飞时根据MPU6050传感器采集到的姿态信息推导出飞行器的姿态角和重心坐标之间存在如下的关系:
(2)力学平衡方程的建立
根据四旋翼飞行器在空间的悬停状态时,升力与重力大小相等,反扭矩和为0以及绕x轴和y轴的力矩平衡方程,得到在不平衡负载下,飞行器在空中悬停的受力平衡方程组;
其中,M1、M2、M3和M4分别为四只螺旋桨对飞行器产生的反扭矩,四只电机到机身中心的距离均为l,F1~F4分别为四只螺旋桨产生的升力;
(3)角速度平衡方程的建立
ω1+ω3=ω2+ω4=ω (3)
其中,ω1,ω2,ω3,ω4分别表示四只螺旋桨在悬停状态下的转速,ω表示1号和3号或者2号和4号螺旋桨的转速和;
为了同时保证式(2)和式(3)同时成立,可以对式(3)做如下修改:
(ω0+Δω13)+(ω0-Δω13)=(ω0+Δω24)+(ω0-Δω24) (4)
其中,ω0是飞行器在平衡负载下每只螺旋桨的平均转速,Δω13和Δω24分别表示1、3号螺旋桨和2、4号螺旋桨的角速度变化量;由式(2)和式(4),可以得出Δω13和Δω24分别为:
那么在不平衡负载下,四只螺旋桨的转速做下面的修正,达到预期位置和姿态;
4)不平衡负载下飞行器的运动学分析
在不平衡负载下飞行器的动力学模型为:
式(7)中,Jx,Jy和Jz分别表示飞行器沿x轴、y轴和z轴方向上的转动惯量;根据以上论述,结合欧拉公式其中M=(M1,M2,M3,M4),J=(Jx,Jy,Jz),ω=(ω1,ω2,ω3,ω4)分别表示扭矩向量、转动惯量向量和转速向量, 得到不平衡负载下的四旋翼飞行器的整体动力学模型:
θ,ψ分别表示飞行器沿x轴、y轴和z轴方向的姿态变化角度,其一阶导数均表示对应的角速度,二阶导数表示角加速度,Jx、Jy和Jz分别是飞行器三个轴的转动惯量;Sx、Sy和Sz分别为飞行器沿x、y和z轴方向上的位移,其一阶导数表示速度,二阶导数表示加速度;可以发现,当重心位置G(xg,yg,0)在机体坐标原点的时候,即当xg=0,yg=0,时,该系统模型就和一般的平衡负载模型一致;
5)姿态解算
四旋翼飞行器飞行过程中期望的姿态角为0,应用欧拉角方法将四旋翼飞行器在空中飞行的姿态信息反馈到主控制器中,根据期望输出与实际输出的误差设计串级PID控制器控制,计算出控制量,转化为相应的控制信号经驱动电路后驱动四个电机工作,保持四旋翼飞行器稳定飞行。
2.根据权利要求1所述的一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法,其特征在于,所述串级PID控制器考虑到四旋翼飞行器在工作过程中受到非线性的因素的影响,在串级PID的基础设计了神经网络自适应控制算法逼近系统的非线性部分,算法设计如下:
将公式(8)所描述的四旋翼飞行器方写成如下形式:
y=x
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711245553.3A CN107844124B (zh) | 2017-12-01 | 2017-12-01 | 一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711245553.3A CN107844124B (zh) | 2017-12-01 | 2017-12-01 | 一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107844124A CN107844124A (zh) | 2018-03-27 |
CN107844124B true CN107844124B (zh) | 2020-07-14 |
Family
ID=61664324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711245553.3A Active CN107844124B (zh) | 2017-12-01 | 2017-12-01 | 一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107844124B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375988A (zh) * | 2018-05-25 | 2018-08-07 | 哈尔滨工业大学 | 一种带有不平衡负载的四旋翼无人机位姿控制方法 |
CN110231822B (zh) * | 2019-06-04 | 2022-04-05 | 哈尔滨工程大学 | 一种舰船用变输出约束式无模型自适应航向控制方法 |
CN111781820B (zh) * | 2020-07-24 | 2023-06-02 | 滨州学院 | 基于重心动力学的航模pid算法控制方法 |
CN112882379B (zh) * | 2021-01-13 | 2022-09-30 | 哈尔滨工业大学 | 一种飞机纵向重心调配控制方法 |
CN113741502B (zh) * | 2021-09-15 | 2023-08-25 | 重庆大学 | 输出约束下带负载四旋翼无人机的神经网络控制方法 |
CN117694841B (zh) * | 2024-02-06 | 2024-04-30 | 成都中医药大学 | 一种重心轨迹数据获取方法 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4232977B2 (ja) * | 2004-12-06 | 2009-03-04 | 日本輸送機株式会社 | 攪拌搬送方法、及び搬送装置 |
US20160272310A1 (en) * | 2014-12-04 | 2016-09-22 | Elwha Llc | Reconfigurable unmanned aircraft system |
US9919797B2 (en) * | 2014-12-04 | 2018-03-20 | Elwha Llc | System and method for operation and management of reconfigurable unmanned aircraft |
US9469394B2 (en) * | 2015-03-10 | 2016-10-18 | Qualcomm Incorporated | Adjustable weight distribution for drone |
CN106444368A (zh) * | 2015-11-18 | 2017-02-22 | 南京航空航天大学 | 具有输入非线性的近空间飞行器预设性能姿态跟踪控制方法 |
CN105447269B (zh) * | 2015-12-15 | 2019-01-11 | 西北工业大学 | 一种飞机在定常过载状态下的非结构质量配平方法 |
CN105759832A (zh) * | 2016-05-20 | 2016-07-13 | 武汉科技大学 | 一种基于反演法的四旋翼飞行器滑模变结构控制方法 |
CN106200665A (zh) * | 2016-08-25 | 2016-12-07 | 东北大学 | 携带不确定负载的四轴飞行器的建模与自适应控制方法 |
CN106647792B (zh) * | 2017-01-25 | 2019-09-03 | 天津大学 | 用于无人机悬挂负载系统的抗扰控制方法 |
CN107144401A (zh) * | 2017-04-18 | 2017-09-08 | 南京信息工程大学 | 一种四旋翼飞行器转动惯量测量方法 |
CN106933104B (zh) * | 2017-04-21 | 2020-05-19 | 苏州工业职业技术学院 | 一种基于dic-pid的四旋翼飞行器姿态与位置的混合控制方法 |
CN106972798A (zh) * | 2017-05-27 | 2017-07-21 | 湖南大学 | 一种抑制电机不平衡负载下转矩脉动的控制方法及系统 |
-
2017
- 2017-12-01 CN CN201711245553.3A patent/CN107844124B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107844124A (zh) | 2018-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107844124B (zh) | 一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法 | |
CN105159305B (zh) | 一种基于滑模变结构的四旋翼飞行控制方法 | |
Ke et al. | Design and implementation of a hybrid UAV with model-based flight capabilities | |
CN106707754B (zh) | 一种基于切换系统的货物搬运旋翼无人机建模及自适应控制方法 | |
CN107765553B (zh) | 针对旋翼无人机吊挂运输系统的非线性控制方法 | |
Govdeli et al. | Unsteady aerodynamic modeling and control of pusher and tilt-rotor quadplane configurations | |
Kui et al. | Sliding mode control for a quadrotor slung load system | |
Zhang et al. | Three-dimensional fuzzy control of mini quadrotor UAV trajectory tracking under impact of wind disturbance | |
CN105759832A (zh) | 一种基于反演法的四旋翼飞行器滑模变结构控制方法 | |
CN104765272A (zh) | 一种基于pid神经元网络控制(pidnn)的四旋翼飞行器控制方法 | |
CN111026160B (zh) | 一种四旋翼无人机轨迹跟踪控制方法 | |
CN108803648A (zh) | 无人飞行器鲁棒姿态控制方法、装置及电子设备 | |
Panish et al. | Transition trajectory optimization for a tiltwing VTOL aircraft with leading-edge fluid injection active flow control | |
CN111061282A (zh) | 基于能量法的四旋翼无人机吊挂飞行系统控制方法 | |
Pedro et al. | Direct adaptive neural control of a quadrotor unmanned aerial vehicle | |
Hegde et al. | Transition flight modeling and robust control of a VTOL unmanned quad tilt-rotor aerial vehicle | |
Wang et al. | Stability analysis of tailsitters in vertical takeoff and landing flights | |
CN107678442A (zh) | 一种基于双模型下的四旋翼自主着船控制方法 | |
Kiani et al. | A new bio-inspired flying concept: The quad-flapper | |
Carter | Adaptive control of the transition from vertical to horizontal flight regime of a quad-tailsitter uav | |
Xu et al. | Modelling and hovering control of a novel multi-tandem ducted fan vehicle | |
CN111650954A (zh) | 基于深度学习的四旋翼无人机地面效应补偿降落控制方法 | |
CN116088549B (zh) | 一种尾座式垂直起降无人机姿态控制方法 | |
Leylek et al. | Use of compliant hinges to tailor flight dynamics of unmanned aircraft | |
Yao et al. | Trajectory tracking controller based on PID-NLADRC |
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 |