CN102935898A - 空间纯引力轨道万有引力摄动在轨飞行验证方法 - Google Patents
空间纯引力轨道万有引力摄动在轨飞行验证方法 Download PDFInfo
- Publication number
- CN102935898A CN102935898A CN2012105060387A CN201210506038A CN102935898A CN 102935898 A CN102935898 A CN 102935898A CN 2012105060387 A CN2012105060387 A CN 2012105060387A CN 201210506038 A CN201210506038 A CN 201210506038A CN 102935898 A CN102935898 A CN 102935898A
- Authority
- CN
- China
- Prior art keywords
- orbit
- perturbation
- universal gravitation
- satellite
- flight
- 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
- 238000012795 verification Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 title abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 24
- 238000010200 validation analysis Methods 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 230000005484 gravity Effects 0.000 abstract description 4
- 238000007689 inspection Methods 0.000 abstract 1
- 230000005656 quantum gravity Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000005658 nuclear physics Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种空间纯引力轨道万有引力摄动在轨飞行验证方法,1,将纯引力飞行卫星在模式一下运行一定周期,记录轨道观测数据,2,将卫星绕径向主轴做180°旋转,作为模式二,再次运行与步骤1中相同的周期并记录轨道观测数据,3,分别得出两种模式下的轨道半长轴长期变化,做差得到两种模式下轨道半长轴变化量之差的实际观测值;4,通过轨道摄动方程,得出两种模式下万有引力摄动导致的轨道半长轴变化量之差的理论预测值;5,对比理论预测值与实际观测值,验证理论预测值的准确性。本发明在纯引力飞行任务初始化阶段,利用卫星的不同飞行状态,对纯引力轨道万有引力摄动这一最主要的摄动因素,进行在轨验证,为万有引力摄动计算方法提供检验。
Description
技术领域
本发明涉及航天动力学技术领域,尤其是涉及一种空间纯引力轨道万有引力摄动在轨飞行验证方法。
背景技术
一些以基础物理实验为任务目标的空间任务,如探测引力波和检验广义相对论的LISA和ASTROD任务,需要参考质量沿着纯引力轨道飞行(参见期刊《经典与量子引力》(Classical andQuantum Gravity)2003年第20卷的文章“LISA的集成模型(The LISA integrated model)”和期刊《原子核物理B》(Nuclear Physics B)2007年第166卷153-158页文章“ASTROD(激光天文动力学)and ASTROD I”)。同样,利用参考质量沿着近地纯引力轨道飞行,并获取参考质量的纯引力轨道,能够用于精确测量地球重力场(参见期刊《国际宇航科学院院刊》(ActaAstronautica)2012年特刊文章“采用精密编队飞行技术获取纯引力轨道(Acquirement of puregravity orbit using precision formation flying technology)”)。对于这些任务的科学目标而言,卫星作用在内部参考质量上的万有引力是一个主要的干扰力,影响纯引力轨道的性能水平(参见《经典与量子引力》(Classical and Quantum Gravity)2004年第21卷第5期S653-S660页的文章“当前的LISA残余加速度误差估计(Current error estimates for LISA spurious accelerations)”)。精确掌握纯引力轨道的万有引力摄动,是提高纯引力轨道性能的基础。
现有技术中,LISA模型团队建立了万有引力干扰的数值计算方法,采用卫星有限单元模型提供的结点质量和位置,并将每个单元近似为质点计算其对参考质量的引力、力矩和梯度作用,然后对所有单元求和得到整体量(参见《经典与量子引力》(Classical and Quantum Gravity)2005年第22卷第10期S395-S402页的文章“LISA自引力分析模型(Self-gravity modeling forLISA)”)。谷振丰等提出了利用专门的补偿质量抑制万有引力摄动的方法,即在卫星上安装补偿质量块,调整卫星的质量分布,降低卫星对参考质量块的万有引力(参见《国际宇航科学院院刊》(Acta Astronautica)2012年第81卷第2期635-644页的文章“内编队系统纯引力轨道万有引力摄动补偿抑制(Compensation of gravitational attraction disturbance to pure gravity orbit forInner Formation Flying System)”)。
但是,这些都是基于理论方法、数值计算开展的,并且采用了一些近似方法,需要进行实验验证。只有对万有引力摄动计算和抑制结果进行检验,才能更好地推动纯引力飞行系统的工程实施。
发明内容
本发明的目的在于设计一种新型的空间纯引力轨道万有引力摄动在轨飞行验证方法,解决上述问题。
为了实现上述目的,本发明采用的技术方案如下:
一种空间纯引力轨道万有引力摄动在轨飞行验证方法,包括步骤如下:
步骤1、将纯引力轨道飞行卫星在模式一下运行一定周期,观测卫星轨道,并记录轨道观测数据;
步骤2、将步骤1中所述的卫星绕着沿轨道面径向的主轴做180°旋转,从而将飞行方向的万有引力摄动反向,形成模式二;所述的卫星在所述模式二下再次运行与步骤1中相同的周期,同样精确观测卫星轨道,并记录轨道观测数据;
步骤3、根据步骤1和步骤2的轨道观测数据,分别得出两种模式下的轨道半长轴长期变化,做差,得到两种模式下的卫星轨道半长轴变化量之差的实际观测值;
步骤4、根据两种模式下的万有引力摄动计算结果,通过卫星轨道摄动方程得出两种模式下万有引力摄动导致的轨道半长轴变化量计算值,进而得到两种模式下轨道半长轴变化量之差的理论预测值;
步骤5、将步骤4所得的理论预测值与步骤3所得的实际观测值进行对比,如果两者一致,则验证了理论预测值的准确性,从而间接验证了万有引力摄动计算结果的准确性,说明了万有引力摄动计算方法的有效性;否则,说明万有引力摄动计算方法无效。
优选的,所述的模式一,是指纯引力轨道飞行卫星飞行方向的万有引力摄动沿着前进方向,又称为标称模式,该模式下万有引力摄动会使得参考质量轨道半长轴逐渐增大;
所述的模式二,是指纯引力轨道飞行卫星绕着沿轨道面径向的主轴做180°旋转,从而将飞行方向的万有引力摄动反向,该模式下万有引力摄动会使得参考质量轨道半长轴逐渐减小。
优选的,步骤1中,所述的一定周期,是指15日至1个月。
优选的,所述卫星轨道摄动方程是指:
其中ainner为参考质量轨道半长轴,einner为参考质量轨道偏心率,finner为参考质量轨道真近点角,ninner为参考质量轨道角速率;fx为万有引力摄动在参考质量轨道径向上的分量,fy为万有引力摄动在参考质量轨道飞行方向上的分量。
本发明的技术方案中中步骤1、2、3作为整体,与步骤4之间是顺序不分先后的;既可以先进行步骤1、2、3,再进行步骤4;也可以先进行步骤4,再进行步骤1、2、3;也可以在进行步骤1、2、3的同时进行步骤4。
本发明的目的是在纯引力轨道飞行卫星中,根据轨道动力学原理对纯引力轨道万有引力摄动进行在轨飞行验证。
本发明空间纯引力轨道万有引力摄动在轨飞行验证方法包括:将纯引力轨道飞行卫星在两种不同模式下运行、利用轨道观测确定两种模式下卫星轨道半长轴变化量之差、利用万有引力摄动计算结果预测两种模式下卫星轨道半长轴变化量之差、比较两种模式下卫星轨道半长轴变化量之差的理论预测值与实际观测值。
所述的纯引力轨道飞行卫星在两种不同模式下运行,是先将卫星在标称模式下运行一定周期,然后将卫星绕着沿轨道面径向的主轴做180°旋转,从而将飞行方向的万有引力摄动反向,在这第二种模式下再次运行相同的周期。所述的利用轨道观测确定两种模式下卫星轨道半长轴变化量之差,是对两种模式下运行的卫星进行轨道观测,然后利用轨道观测数据分析得出两种模式下的轨道半长轴长期变化,做差得到两种模式下的卫星轨道半长轴变化量之差的实际观测值。所述的利用万有引力摄动计算结果预测两种模式下卫星轨道半长轴变化量之差,是根据两种模式下的万有引力摄动计算结果,通过卫星轨道摄动方程得出两种模式下万有引力摄动导致的轨道半长轴变化量计算值,进而得到两种模式下轨道半长轴变化量之差的理论预测值。所述的比较两种模式下卫星轨道半长轴变化量之差的理论预测值与实际观测值,是将两种模式下卫星轨道半长轴变化量之差的理论预测值与实际观测值进行对比,如果两者一致,则验证了理论预测值的准确性,从而间接验证了万有引力摄动计算结果的准确性,也说明了万有引力摄动计算方法的有效性。
本发明的有益效果可以总结如下:
1,本发明在纯引力轨道飞行卫星入轨之后的初始化阶段,可以利用卫星的不同飞行状态,对纯引力轨道万有引力摄动这一最主要的摄动因素,进行在轨验证,为万有引力摄动计算方法提供检验。
2,本发明方法简单、实施成本低廉,测量结果精确。
具体实施方式
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
一种空间纯引力轨道万有引力摄动在轨飞行验证方法,包括步骤如下:
步骤1、将纯引力轨道飞行卫星在模式一下运行一定周期,观测卫星轨道,并记录轨道观测数据;
步骤2、将步骤1中所述的卫星绕着沿轨道面径向的主轴做180°旋转,从而将飞行方向的万有引力摄动反向,形成模式二;所述的卫星在所述模式二下再次运行与步骤1中相同的周期,同样精确观测卫星轨道,并记录轨道观测数据;
步骤3、根据步骤1和步骤2的轨道观测数据,分别得出两种模式下的轨道半长轴长期变化,做差,得到两种模式下的卫星轨道半长轴变化量之差的实际观测值;
步骤4、根据两种模式下的万有引力摄动计算结果,通过卫星轨道摄动方程得出两种模式下万有引力摄动导致的轨道半长轴变化量计算值,进而得到两种模式下轨道半长轴变化量之差的理论预测值;
步骤5、将步骤4所得的理论预测值与步骤3所得的实际观测值进行对比,如果两者一致,则验证了理论预测值的准确性,从而间接验证了万有引力摄动计算结果的准确性,说明了万有引力摄动计算方法的有效性;否则,说明万有引力摄动计算方法无效。
其中,所述的模式一,是指纯引力轨道飞行卫星飞行方向的万有引力摄动沿着前进方向,又称为标称模式,该模式下万有引力摄动会使得参考质量轨道半长轴逐渐增大;
所述的模式二,是指纯引力轨道飞行卫星绕着沿轨道面径向的主轴做180°旋转,从而将飞行方向的万有引力摄动反向,该模式下万有引力摄动会使得参考质量轨道半长轴逐渐减小。
步骤1中,所述的一定周期,是指15日至1个月。
所述卫星轨道摄动方程是指:
其中ainner为参考质量轨道半长轴,einner为参考质量轨道偏心率,finner为参考质量轨道真近点角,ninner为参考质量轨道角速率;fx为万有引力摄动在参考质量轨道径向上的分量,fy为万有引力摄动在参考质量轨道飞行方向上的分量。
以下举例说明具体的验证过程:
纯引力轨道飞行卫星包含一个参考质量,卫星具有相对位置测量、精密轨道跟踪、三轴姿态稳定控制能力。
纯引力轨道飞行卫星中万有引力摄动对参考质量的作用是连续的,导致参质量逐渐偏离理想的纯引力轨道,则同样使得参考质量轨道半长轴发生长期变化。根据卫星轨道摄动方程可知,在万有引力摄动下参考质量轨道半长轴的变化为
其中ainner为参考质量轨道半长轴,einner为参考质量轨道偏心率,finner为参考质量轨道真近点角,ninner为参考质量轨道角速率;fx为万有引力摄动在参考质量轨道径向上的分量,fy为万有引力摄动在参考质量轨道飞行方向上的分量。
设计纯引力轨道飞行卫星的偏心率einner为0,则轨道半长轴在万有引力摄动作用下的变化率为
纯引力轨道飞行卫星具有两种不同的工作模式:
模式一,纯引力轨道飞行卫星飞行方向的万有引力摄动沿着前进方向。由式(2)可知,该模式下万有引力摄动会使得参考质量轨道半长轴逐渐增大。
模式二,纯引力轨道飞行卫星绕着沿轨道面径向的主轴做180°旋转,从而将飞行方向的万有引力摄动反向。由式(2)可知,该模式下万有引力摄动会使得参考质量轨道半长轴逐渐减小。
在完成任务初始化之后,纯引力轨道飞行卫星先在模式一下运行一个月,并精确测量其轨道。然后纯引力轨道飞行卫星切换到模式二下运行一个月,同样精确测量其轨道。
在完成任务飞行之后,通过对轨道测量数据的事后处理,分别得出纯引力轨道飞行卫星在模式一、模式二下的轨道半长轴变化量,从而得到两种模式下的轨道半长轴变化量之差的实际观测值。
以纯引力轨道飞行卫星万有引力摄动计算结果为基础,利用式(2)计算轨道半长轴变化率,通过数值积分得出两种模式下万有引力摄动导致的轨道半长轴变化量计算值,进而得到两种模式下轨道半长轴变化量之差的理论预测值。
模式一与模式二中,参考质量轨道运动摄动因素的主要差别是飞行方向的万有引力摄动不同,这是两种模式下轨道半长轴变化量之差的主导因素。对比其理论预测值与实际观测值,如果两者一致,则验证了理论预测值的准确性,从而间接验证了万有引力摄动计算结果的准确性,也进一步说明了万有引力摄动计算方法的有效性。
综上可见,本发明在纯引力轨道飞行卫星入轨之后的初始化阶段,可以利用卫星的不同飞行状态,对纯引力轨道万有引力摄动这一最主要的摄动因素,进行在轨验证,为万有引力摄动计算方法提供检验。
以上通过具体的和优选的实施例详细的描述了本发明,但本领域技术人员应该明白,本发明并不局限于以上所述实施例,凡在本发明的精神和原则之内,所作的任何修改、等同替换等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种空间纯引力轨道万有引力摄动在轨飞行验证方法,其特征在于,包括步骤如下:
步骤1、将纯引力轨道飞行卫星在模式一下运行一定周期,观测卫星轨道,并记录轨道观测数据;
步骤2、将步骤1中所述的卫星绕着沿轨道面径向的主轴做180°旋转,从而将飞行方向的万有引力摄动反向,形成模式二;所述的卫星在所述模式二下再次运行与步骤1中相同的周期,同样精确观测卫星轨道,并记录轨道观测数据;
步骤3、根据步骤1和步骤2的轨道观测数据,分别得出两种模式下的轨道半长轴长期变化,做差,得到两种模式下的卫星轨道半长轴变化量之差的实际观测值;
步骤4、根据两种模式下的万有引力摄动计算结果,通过卫星轨道摄动方程得出两种模式下万有引力摄动导致的轨道半长轴变化量计算值,进而得到两种模式下轨道半长轴变化量之差的理论预测值;
步骤5、将步骤4所得的理论预测值与步骤3所得的实际观测值进行对比,如果两者一致,则验证了理论预测值的准确性,从而间接验证了万有引力摄动计算结果的准确性,说明了万有引力摄动计算方法的有效性;否则,说明万有引力摄动计算方法无效。
2.根据权利要求1所述的空间纯引力轨道万有引力摄动在轨飞行验证方法,其特征在于:所述的模式一,是指纯引力轨道飞行卫星飞行方向的万有引力摄动沿着前进方向,又称为标称模式,该模式下万有引力摄动会使得参考质量轨道半长轴逐渐增大;
所述的模式二,是指纯引力轨道飞行卫星绕着沿轨道面径向的主轴做180°旋转,从而将飞行方向的万有引力摄动反向,该模式下万有引力摄动会使得参考质量轨道半长轴逐渐减小。
3.根据权利要求2所述的空间纯引力轨道万有引力摄动在轨飞行验证方法,其特征在于:步骤1中,所述的一定周期,是指15日至1个月。
4.根据权利要求1述的空间纯引力轨道万有引力摄动在轨飞行验证方法,其特征在于:所述卫星轨道摄动方程是指:
其中ainner为参考质量轨道半长轴,einner为参考质量轨道偏心率,finner为参考质量轨道真近点角,ninner为参考质量轨道角速率;fx为万有引力摄动在参考质量轨道径向上的分量,fy为万有引力摄动在参考质量轨道飞行方向上的分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210506038.7A CN102935898B (zh) | 2012-11-30 | 2012-11-30 | 空间纯引力轨道万有引力摄动在轨飞行验证方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210506038.7A CN102935898B (zh) | 2012-11-30 | 2012-11-30 | 空间纯引力轨道万有引力摄动在轨飞行验证方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102935898A true CN102935898A (zh) | 2013-02-20 |
CN102935898B CN102935898B (zh) | 2015-05-06 |
Family
ID=47694834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210506038.7A Expired - Fee Related CN102935898B (zh) | 2012-11-30 | 2012-11-30 | 空间纯引力轨道万有引力摄动在轨飞行验证方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102935898B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104596520A (zh) * | 2015-01-20 | 2015-05-06 | 中国人民解放军63920部队 | 一种基于半长轴滤波值连续性的航天器轨道机动检测方法 |
CN104608940A (zh) * | 2014-12-22 | 2015-05-13 | 汪宇 | 引力场曲率驱动引擎的方法 |
CN113987407A (zh) * | 2021-10-28 | 2022-01-28 | 中国西安卫星测控中心 | 一种地球静止卫星非球形引力切向漂移加速度计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2691129A1 (fr) * | 1992-05-18 | 1993-11-19 | Aerospatiale | Satellite géodésique de faible masse à rétroréflecteur à correction d'aberration de vitesse. |
CN102360084A (zh) * | 2011-07-26 | 2012-02-22 | 清华大学 | 基于四面体质元划分的纯引力轨道万有引力干扰计算方法 |
CN102785785A (zh) * | 2012-09-10 | 2012-11-21 | 清华大学 | 利用外航天器自旋抑制纯引力轨道万有引力干扰的方法 |
-
2012
- 2012-11-30 CN CN201210506038.7A patent/CN102935898B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2691129A1 (fr) * | 1992-05-18 | 1993-11-19 | Aerospatiale | Satellite géodésique de faible masse à rétroréflecteur à correction d'aberration de vitesse. |
CN102360084A (zh) * | 2011-07-26 | 2012-02-22 | 清华大学 | 基于四面体质元划分的纯引力轨道万有引力干扰计算方法 |
CN102785785A (zh) * | 2012-09-10 | 2012-11-21 | 清华大学 | 利用外航天器自旋抑制纯引力轨道万有引力干扰的方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104608940A (zh) * | 2014-12-22 | 2015-05-13 | 汪宇 | 引力场曲率驱动引擎的方法 |
CN104596520A (zh) * | 2015-01-20 | 2015-05-06 | 中国人民解放军63920部队 | 一种基于半长轴滤波值连续性的航天器轨道机动检测方法 |
CN104596520B (zh) * | 2015-01-20 | 2017-07-28 | 中国人民解放军63920部队 | 一种基于半长轴滤波值连续性的航天器轨道机动检测方法 |
CN113987407A (zh) * | 2021-10-28 | 2022-01-28 | 中国西安卫星测控中心 | 一种地球静止卫星非球形引力切向漂移加速度计算方法 |
CN113987407B (zh) * | 2021-10-28 | 2022-07-08 | 中国西安卫星测控中心 | 一种地球静止卫星非球形引力切向漂移加速度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102935898B (zh) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Williams et al. | Tides on the Moon: Theory and determination of dissipation | |
Searcy et al. | Magnetometer-only attitude determination using novel two-step Kalman filter approach | |
CN103017774B (zh) | 单探测器脉冲星导航方法 | |
Ciufolini et al. | Testing gravitational physics with satellite laser ranging | |
WO2019191294A1 (en) | Methods and apparatus for in-situ measurements of atmospheric density | |
Renzetti | First results from LARES: An analysis | |
Santoli et al. | ISA, a High Sensitivity Accelerometer in the Interplanetary Space: Updates after the Near-Earth Commissioning Phase of Italian Spring Accelerometer–ISA | |
CN105300597A (zh) | 三轴气浮台质心调平衡方法及装置 | |
CN108959734A (zh) | 一种基于实时递推太阳光压力矩辨识方法及系统 | |
CN102935898B (zh) | 空间纯引力轨道万有引力摄动在轨飞行验证方法 | |
Navabi et al. | A comparative study of dynamics models for satellite formation flying-Cartesian ordinary differential equations description | |
CN103091723B (zh) | 降低重力卫星质心调整误差对地球重力场精度影响的方法 | |
CN102735265B (zh) | 一种基于陀螺漂移估计值的星敏感器周期性故障检测方法 | |
Sun et al. | Researching on the compensation technology of rotating mechanism error in single-axis rotation strapdown inertial navigation system | |
Congedo et al. | Testing the strong equivalence principle with spacecraft ranging towards the nearby Lagrangian points | |
CN102785785B (zh) | 利用外航天器自旋抑制纯引力轨道万有引力干扰的方法 | |
Gutsche et al. | Podcast: Precise orbit determination software for leo satellites | |
Adnane et al. | Reliable Kalman filtering for satellite attitude estimation under gyroscope partial failure | |
Ho et al. | Extended two-body problem for rotating rigid bodies | |
Yu et al. | Real-time onboard orbit determination using GPS navigation solutions | |
Cheon et al. | Disturbance reduction on the small satellite actuator | |
Cha et al. | Coarse alignment of lunar exploration rover using accelerometer and sun sensor | |
De Marchi et al. | Space tests of the strong equivalence principle: BepiColombo and the Sun–Earth Lagrangian points opportunity | |
Aleshin et al. | Horizontal stabilization of the two-degree-of-freedom platform of a uniaxial wheeled module tracking a given trajectory over an underlying surface | |
CN102636817A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150506 |