CN109470241B - 一种具备重力扰动自主补偿功能的惯性导航系统及方法 - Google Patents
一种具备重力扰动自主补偿功能的惯性导航系统及方法 Download PDFInfo
- Publication number
- CN109470241B CN109470241B CN201811407237.6A CN201811407237A CN109470241B CN 109470241 B CN109470241 B CN 109470241B CN 201811407237 A CN201811407237 A CN 201811407237A CN 109470241 B CN109470241 B CN 109470241B
- Authority
- CN
- China
- Prior art keywords
- gravity
- information
- inertial
- current position
- carrier
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明涉及一种具备重力扰动自主补偿功能的惯性导航系统及方法,其技术特点在于:包括平台式惯性导航系统、惯性信息测量模块和重力梯度测量模块;所述平台式惯性导航系统包括三轴平台框架、平台台体和平台底座;所述惯性信息测量模块包括惯性测量单元和重力测量模块;所述惯性测量单元输出线运动和角运动信息进行导航解算;所述三轴平台框架将平台台体稳定在地理坐标下,同时输出载体位置、速度、航向、姿态导航信息;所述重力测量模块和重力梯度测量模块同步测量载体当前位置的重力异常信息和全张量重力梯度信息,计算出载体当前位置的重力矢量信息,并将重力矢量信息用于惯性测量单元的导航解算。本发明能够提高惯性导航系统精度。
Description
技术领域
本发明属于惯性导航技术领域,涉及惯性导航系统,尤其是一种具备重力扰动自主补偿功能的惯性导航系统及方法。
背景技术
惯性导航是根据牛顿惯性原理,利用陀螺仪建立导航坐标系,利用加速度计测量载体运动加速度,经过一次积分得到运载体的速度,再经过一次积分得到运载体的地理位置的导航方法。传统的固定方位半解析式惯性导航系统基本原理如图3所示,通过三只单自由度陀螺仪和框架系统构成三轴稳定平台,通过陀螺仪敏感的角运动信息施加指令驱动力矩电机以达到隔离载体角运动的作用。
在系统定位方面,在稳定平台始终保持水平指北的前提下,在平台上水平安装三只相互垂直的加速度计,三只加速度计敏感轴分布沿平台北向轴、东向轴和垂向轴安装,用来分别测量载体北向、东向和垂向加速度。这些加速度信号除含有运载体相对地球的运动加速度以外,还含有哥氏加速度、离心加速度及重力矢量等有害加速度。因此在解算过程中进行有害加速度ABx、ABy和ABz的补偿,然后经过一次积分可得运载体速度分量,即
式中Vx(t)、Vy(t)和Vz(t)表示载体在地理坐标系下的北向、东向和垂向实时速度,Vx0、Vy0和Vz0表示载体在地理坐标系下的初始时刻北向和东向和垂向速度,Ax、Ay和Az表示北向、东向和垂向加速度计的输出加速度,ABx、ABy和ABz表示需要补偿的北向、东向和垂向有害加速度。传统惯导系统的有害加速度项为:
将所得到速度分量在经过一次积分运算,即可得到运载体相对地球的经度、纬度和高度的变化量。如果输入起始点的经度、纬度和高度,就可得到运载体所在地的经度、纬度和高度。
有害加速度的补偿是惯性导航系统解算过程中重要环节之一,而由于加速度计无法区分运动加速度和重力加速度,一般采用正常重力模型来代替实际重力值在解算中予以补偿。但由于实际地球的形状不规则、内部质量分布不均匀、参考椭球与大地水准面不完全重合等原因,重力矢量和正常重力模型的计算值有所偏差,这种偏差称为重力扰动。重力扰动分为重力异常和垂向偏差两部分,其中重力异常是实际重力矢量模值与模型重力模值之差,垂向偏差为实际重力方向与模型重力方向的夹角。重力扰动矢量可表示为:
δg=|g|-g0
g0=[0 0 -g0]T
g=[gx gy -gz]T
式中δg表示重力异常值,|g|表示当前位置实际重力矢量模值,g0表示正常重力模型下的计算值,gx表示重力矢量东向分量,gy表示重力矢量北向分量,gz表示重力矢量垂向分量。
在考虑重力扰动的条件下,则惯导系统中有害加速度可表示为:
随着惯性元件精度的不断提高和惯性导航系统技术的不断发展,重力扰动已成为高精度惯性导航系统的一项主要误差来源。传统的应对方法是采用专用的重力仪或重力梯度仪实时测量重力矢量,对惯性导航系统进行补偿。惯导系统与重力信息测量系统是分离的,存在系统复杂、可靠性差、成本高、数据信息传输延迟的问题。并且由于二者空间不一致,存在一定的测量误差。
发明内容
本发明的目的在于克服现有技术的不足,提出一种具备重力扰动自主补偿功能的惯性导航系统及方法,具备重力实时测量功能和重力梯度实时测量功能,利用重力、重力梯度的实时测量值计算有害加速度中的重力扰动项并在解算中予以补偿,以实现提高惯性导航系统精度的目的。
本发明解决其现实问题是采取以下技术方案实现的:
一种具备重力扰动自主补偿功能的惯性导航系统,包括平台式惯性导航系统、惯性信息测量模块和重力梯度测量模块;所述平台式惯性导航系统包括三轴平台框架、平台台体和平台底座;所述惯性信息测量模块包括惯性测量单元和重力测量模块;所述平台底座上安装有平台框架,在该平台框架上安装有平台台体;在该平台台体上固装有惯性测量单元和重力测量模块、重力梯度测量模块;所述惯性测量单元输出线运动和角运动信息进行导航解算;所述三轴平台框架将平台台体稳定在地理坐标下,同时输出载体位置、速度、航向、姿态导航信息;所述重力测量模块和重力梯度测量模块同步测量载体当前位置的重力异常信息和全张量重力梯度信息,计算出载体当前位置的重力矢量信息,并将重力矢量信息用于惯性测量单元的导航解算。
一种具备重力扰动自主补偿功能的惯性导航方法,包括以下步骤:
步骤1、惯性导航系统在执行导航功能的同时,同步完成重力和重力梯度实时测量;
步骤2、根据步骤1获得的重力梯度张量的测量值结合惯导系统解算出当前载体的位置信息,实时计算载体当前位置的重力矢量水平双分量;
步骤3、通过当前位置的重力异常值,结合步骤2中得到的当前位置的重力矢量水平分量,计算载体当前位置的重力矢量垂向分量,从而得到地理坐标系下的重力矢量;
步骤4、利用地理坐标系下的重力矢量三分量,在惯导系统解算中更高精度的补偿有害加速度中的重力矢量。
而且,所述步骤2的具体步骤包括:
(1)在系统正常工作时,平台上的重力测量组件输出当前位置的重力异常值δg,重力梯度测量组件输出当前位置的重力梯度全张量信息Γ,其中Γ在数学上可表示为:
式中,Γij(i,j=x,y,z)为梯度张量的分量,表示重力分量gi在j方向上的变化率;
(2)将所得到重力梯度信息中的水平分量对位置经过一次积分运算,即可得到当前位置重力矢量水平分量的变化量;如果输入起始点的重力矢量水平分量gx0及gy0,就可得到运载体所在地重力矢量的水平分量Δgx及Δgy,即:
式中gx(p)、gy(p)表示当前位置重力矢量的北向、东向分量,gx0、gy0表示初始位置重力矢量的北向、东向分量。
而且,所述步骤3的具体步骤包括:
(1)重力仪输出的载体当前位置的重力异常值δg,即实际重力矢量与模型重力矢量的模之差;模型重力矢量模g0的计算式为:
式中,ge=9.78049m/s2表示赤道海平面上的重力加速度;
(2)将重力仪输出的载体当前位置的重力异常值δg与模型重力矢量模g0相加,可得到载体当前位置的重力矢量模值|g|,即:
|g(p)|=g0(p)+δg(p)
(3)再将得到的载体当前位置的重力矢量模值|g(p)|进行坐标分解,即可得到重力矢量在垂向上的投影,即:
本发明的优点和有益效果:
1、本发明在具备惯性导航功能的同时还拥有实时进行重力异常测量、重力梯度测量能力,利用重力异常、重力梯度信息自主补偿惯性导航系统解算中的重力扰动误差,以达到提高惯性导航系统精度的目的。
2、本发明通过此惯性导航系统内部的重力测量组件和重力梯度测量组件计算运载体当前位置的重力矢量信息,补偿惯导解算过程中的重力扰动项,以提高系统的导航精度。
3、本发明的重力测量组件、重力梯度测量组件和加速度计测量组件共用惯性导航系统的平台,可达到减小系统体积、提高系统可靠性、降低系统成本的作用。
4、本发明的重力测量组件、重力梯度测量组件和惯性测量单元在空间上紧密安装在一起,可在最大程度上降低由于重力信息测量组件与惯性测量单元空间不一致带来的测量误差。
5、本发明的惯导系统还可利用重力和重力梯度的实时测量值与事先建立的重力数据库或重力梯度数据库进行匹配比对,利用重力信息与地理坐标信息的相关特性,估算运载体当前地理位置,以提高系统的位置精度。
附图说明
图1是本发明的新型惯性导航系统组成图;
图2是本发明的具备补偿重力扰动能力的新型惯性导航系统原理图;
图3是本发明的固定方位半解析式惯导系统原理图;
附图标记说明:
1-1-惯性测量单元;1-2-重力测量模块;2-重力梯度测量模块;3-平台框架;4-平台台体;5-平台底座。
具体实施方式
以下结合附图对本发明实施例作进一步详述:
一种具备重力扰动自主补偿功能的惯性导航系统,如图1所示,包括平台式惯性导航系统、惯性信息测量模块1和重力梯度测量模块2;所述平台式惯性导航系统包括隔离载体运动的三轴平台框架3、稳定在地理坐标系下的平台台体4和平台底座5;惯性信息测量模块1包括惯性测量单元1-1和重力测量模块1-2;
所述平台底座上安装有平台框架,在该平台框架上安装有平台台体;在该平台台体上固装有惯性测量单元和重力测量模块、重力梯度测量模块;所述惯性测量单元输出线运动和角运动信息进行导航解算;所述三轴平台框架将平台台体稳定在地理坐标下,同时输出载体位置、速度、航向、姿态等导航信息;所述重力测量模块和重力梯度测量模块同步测量载体当前位置的重力异常信息和全张量重力梯度信息,计算出载体当前位置的重力矢量信息,并将重力矢量信息用于惯性测量单元的导航解算,以提高惯性导航系统性能;平台底座则用于承载平台框架、惯性测量单元、重力测量模块和重力梯度测量模块。
本发明的新型惯性导航系统的工作原理为:
在惯性导航系统中增加重力测量组件和重力梯度测量模块,二者与系统中加速度计组件紧密安装,在稳定平台的作用下共同稳定在地理坐标系下,惯导系统正常工作时,稳定平台采用三环固定指北半解析式稳定平台控制方案,可以充分隔离载体的角运动影响,同时将导航用加速度计测量组件、重力测量模块和重力梯度测量模块稳定在地理坐标系下。
一种具备重力扰动自主补偿功能的惯性导航方法,如图2所示,包括以下步骤:
步骤1、惯性导航系统在执行导航功能的同时,同步完成重力和重力梯度实时测量;
步骤2、根据步骤1获得的重力梯度张量的测量值结合惯导系统解算出当前载体的位置信息,实时计算载体当前位置的重力矢量水平双分量;
所述步骤2的具体步骤包括:
(1)在系统正常工作时,平台上的重力测量组件输出当前位置的重力异常值δg,重力梯度测量组件输出当前位置的重力梯度全张量信息Γ,其中Γ在数学上可表示为:
式中,Γij(i,j=x,y,z)为梯度张量的分量,表示重力分量gi在j方向上的变化率;
(2)将所得到重力梯度信息中的水平分量对位置经过一次积分运算,即可得到当前位置重力矢量水平分量的变化量;如果输入起始点的重力矢量水平分量gx0及gy0,就可得到运载体所在地重力矢量的水平分量Δgx及Δgy,即:
式中gx(p)、gy(p)表示当前位置重力矢量的北向、东向分量,gx0、gy0表示初始位置重力矢量的北向、东向分量。
步骤3、通过当前位置的重力异常值,结合步骤2中得到的当前位置的重力矢量水平分量,计算载体当前位置的重力矢量垂向分量,从而得到地理坐标系下的重力矢量;
所述步骤3的具体步骤包括:
(1)重力仪输出的载体当前位置的重力异常值δg,即实际重力矢量与模型重力矢量的模之差;模型重力矢量模g0的计算式为:
式中,ge=9.78049m/s2表示赤道海平面上的重力加速度;
(2)将重力仪输出的载体当前位置的重力异常值δg与模型重力矢量模g0相加,可得到载体当前位置的重力矢量模值|g|,即:
|g(p)|=g0(p)+δg(p)
(3)再将得到的载体当前位置的重力矢量模值|g(p)|进行坐标分解,即可得到重力矢量在垂向上的投影,即:
步骤4、利用地理坐标系下的重力矢量三分量,在惯导系统解算中更高精度的补偿有害加速度中的重力矢量。
通过上述方法,可通过惯性系统上的重力测量组件和重力梯度测量组件得到载体当前位置的重力矢量三分量信息,以此补偿惯导解算有害加速度补偿中的重力扰动项,可达到提高惯性导航系统精度的目的。
此外,通过重力测量组件和重力梯度测量组件提供的重力、重力梯度实时信息与事先建立的重力数据库或重力梯度数据库进行匹配导航,可有效降低惯性导航系统定位误差的发散速度,延长重调周期。相比于现有利用分立的重力仪、重力梯度仪提供的重力信息与惯导信息进行组合导航,本惯导系统的重力、重力梯度测量组件与惯性测量单元空间距离更近,可有效降低二者测量点位置不一致引起的误差,使得重力测量信息更准确,匹配精度更高。
需要强调的是,本发明所述实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (3)
1.一种具备重力扰动自主补偿功能的惯性导航系统,其特征在于:包括平台式惯性导航系统、惯性信息测量模块和重力梯度测量模块;所述平台式惯性导航系统包括三轴平台框架、平台台体和平台底座;所述惯性信息测量模块包括惯性测量单元和重力测量模块;所述平台底座上安装有平台框架,在该平台框架上安装有平台台体;在该平台台体上固装有惯性测量单元和重力测量模块、重力梯度测量模块;所述惯性测量单元输出线运动和角运动信息进行导航解算;所述三轴平台框架将平台台体稳定在地理坐标下,同时输出载体位置、速度、航向、姿态导航信息;所述重力测量模块和重力梯度测量模块同步测量载体当前位置的重力异常信息和全张量重力梯度信息,计算出载体当前位置的重力矢量信息,并将重力矢量信息用于惯性测量单元的导航解算;
所述的一种具备重力扰动自主补偿功能的惯性导航系统的导航方法,包括以下步骤:
步骤1、惯性导航系统在执行导航功能的同时,同步完成重力和重力梯度实时测量;
步骤2、根据步骤1获得的重力梯度张量的测量值结合惯导系统解算出当前载体的位置信息,实时计算载体当前位置的重力矢量水平双分量;
步骤3、通过当前位置的重力异常值,结合步骤2中得到的当前位置的重力矢量水平分量,计算载体当前位置的重力矢量垂向分量,从而得到地理坐标系下的重力矢量;
步骤4、利用地理坐标系下的重力矢量三分量,在惯导系统解算中更高精度的补偿有害加速度中的重力矢量。
2.根据权利要求1所述的一种具备重力扰动自主补偿功能的惯性导航系统的导航方法,其特征在于:所述步骤2的具体步骤包括:
(1)在系统正常工作时,平台上的重力测量组件输出当前位置的重力异常值δg,重力梯度测量组件输出当前位置的重力梯度全张量信息Γ,其中Γ在数学上可表示为:
式中,Γij(i,j=x,y,z)为梯度张量的分量,表示重力分量gi在j方向上的变化率;
(2)将所得到重力梯度信息中的水平分量对位置经过一次积分运算,即可得到当前位置重力矢量水平分量的变化量;如果输入起始点的重力矢量水平分量gx0及gy0,就可得到运载体所在地重力矢量的水平分量Δgx及Δgy,即:
式中gx(p)、gy(p)表示当前位置重力矢量的北向、东向分量,gx0、gy0表示初始位置重力矢量的北向、东向分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811407237.6A CN109470241B (zh) | 2018-11-23 | 2018-11-23 | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811407237.6A CN109470241B (zh) | 2018-11-23 | 2018-11-23 | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109470241A CN109470241A (zh) | 2019-03-15 |
CN109470241B true CN109470241B (zh) | 2020-11-10 |
Family
ID=65674042
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811407237.6A Active CN109470241B (zh) | 2018-11-23 | 2018-11-23 | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109470241B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109883541B (zh) * | 2019-03-26 | 2020-04-17 | 北京交通大学 | 一种密封三轴调姿机构 |
CN110068876B (zh) * | 2019-05-30 | 2021-01-26 | 中国船舶重工集团公司第七0七研究所 | 基于载体自振动航空重力梯度仪运动误差补偿方法 |
CN112304310A (zh) * | 2019-07-23 | 2021-02-02 | 南京航空航天大学 | 一种基于陀螺仪信息的惯性导航方法 |
CN111044041B (zh) * | 2019-12-31 | 2021-04-27 | 北京理工大学 | 基于重力场三维特征的重力辅助惯性导航适配区选取方法 |
CN111141285B (zh) * | 2020-01-06 | 2021-01-08 | 中国自然资源航空物探遥感中心 | 一种航空重力测量装置 |
CN112762924B (zh) * | 2020-12-23 | 2023-07-14 | 北京机电工程研究所 | 基于重力梯度-地形异源数据匹配的导航定位方法 |
CN112729288B (zh) * | 2020-12-23 | 2023-07-14 | 北京机电工程研究所 | 一种重力梯度-地形异源数据匹配的导航定位系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105606093A (zh) * | 2016-01-29 | 2016-05-25 | 北京航空航天大学 | 基于重力实时补偿的惯性导航方法及装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102052922A (zh) * | 2010-11-19 | 2011-05-11 | 中国人民解放军海军工程大学 | 实际重力场影响惯性导航系统扰动重力补偿方法 |
CN103076640B (zh) * | 2013-01-17 | 2015-03-18 | 中国科学院测量与地球物理研究所 | 利用方差-协方差对角张量原理反演地球重力场的方法 |
CN103163562B (zh) * | 2013-02-01 | 2015-05-13 | 中国科学院测量与地球物理研究所 | 基于滤波原理的卫星重力梯度反演方法 |
CN104061932B (zh) * | 2014-06-10 | 2017-04-19 | 中国空间技术研究院 | 一种利用引力矢量和梯度张量进行导航定位的方法 |
CN104374401A (zh) * | 2014-10-15 | 2015-02-25 | 哈尔滨工程大学 | 一种捷联惯导初始对准中重力扰动的补偿方法 |
CN106052694B (zh) * | 2016-07-11 | 2017-03-15 | 中南大学 | 基于重力矢量及其梯度张量对单个运动物体进行定位跟踪的方法 |
-
2018
- 2018-11-23 CN CN201811407237.6A patent/CN109470241B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105606093A (zh) * | 2016-01-29 | 2016-05-25 | 北京航空航天大学 | 基于重力实时补偿的惯性导航方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109470241A (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109470241B (zh) | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 | |
CN108051866B (zh) | 基于捷联惯性/gps组合辅助水平角运动隔离的重力测量方法 | |
JP5068531B2 (ja) | 測定及び記憶された重力傾度を用いて慣性航法測定値の精度を改善する方法及びシステム | |
CN104880189B (zh) | 一种动中通天线低成本跟踪抗干扰方法 | |
CN103925930B (zh) | 一种重力仪双轴陀螺稳定平台航向误差效应的补偿方法 | |
CN111024074B (zh) | 一种基于递推最小二乘参数辨识的惯导速度误差确定方法 | |
CN112432642B (zh) | 一种重力灯塔与惯性导航融合定位方法及系统 | |
CN104374388A (zh) | 一种基于偏振光传感器的航姿测定方法 | |
CN109916395A (zh) | 一种姿态自主冗余组合导航算法 | |
CN111722295B (zh) | 一种水下捷联式重力测量数据处理方法 | |
CN105091907A (zh) | Sins/dvl组合中dvl方位安装误差估计方法 | |
CN103175528A (zh) | 基于捷联惯导系统的捷联罗经姿态测量方法 | |
CN113340298A (zh) | 一种惯导和双天线gnss外参标定方法 | |
CN117053782A (zh) | 一种水陆两栖机器人组合导航方法 | |
Zhang et al. | Research on accuracy enhancement of low-cost MEMS INS/GNSS integration for land vehicle navigation | |
CN110221331B (zh) | 基于状态变换的惯性/卫星组合导航动态滤波方法 | |
CN109084755B (zh) | 一种基于重力视速度与参数辨识的加速度计零偏估计方法 | |
RU2382988C1 (ru) | Бесплатформенная инерциальная система ориентации на "грубых" чувствительных элементах | |
CN113008229A (zh) | 一种基于低成本车载传感器的分布式自主组合导航方法 | |
CN111141285B (zh) | 一种航空重力测量装置 | |
CN112292578B (zh) | 大地水准面测量方法、测量装置、估计装置、计算用数据采集装置 | |
CN112798014A (zh) | 一种基于重力场球谐模型补偿垂线偏差的惯导自对准方法 | |
CN112229421A (zh) | 基于李群最优估计的捷联惯性导航晃动基座粗对准方法 | |
CN114061574B (zh) | 一种基于位置不变约束及零速校正的采煤机定姿定向方法 | |
CN106323226B (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 |