CN1039301C - 卫星姿态控制的测试设备及其应用 - Google Patents

卫星姿态控制的测试设备及其应用 Download PDF

Info

Publication number
CN1039301C
CN1039301C CN92111379A CN92111379A CN1039301C CN 1039301 C CN1039301 C CN 1039301C CN 92111379 A CN92111379 A CN 92111379A CN 92111379 A CN92111379 A CN 92111379A CN 1039301 C CN1039301 C CN 1039301C
Authority
CN
China
Prior art keywords
overbar
vector
satellite
sun
centerdot
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.)
Expired - Fee Related
Application number
CN92111379A
Other languages
English (en)
Other versions
CN1076419A (zh
Inventor
M·苏劳尔
H·比特纳
W·菲希特
H·D·费希尔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Airbus Defence and Space GmbH
Original Assignee
Messerschmitt Bolkow Blohm AG
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Messerschmitt Bolkow Blohm AG filed Critical Messerschmitt Bolkow Blohm AG
Publication of CN1076419A publication Critical patent/CN1076419A/zh
Application granted granted Critical
Publication of CN1039301C publication Critical patent/CN1039301C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/36Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
    • B64G1/363Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using sun sensors
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/242Orbits and trajectories
    • B64G1/2427Transfer orbits
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems
    • B64G1/245Attitude control algorithms for spacecraft attitude control
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/242Orbits and trajectories
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/36Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
    • B64G1/365Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using horizon or Earth sensors
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/36Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
    • B64G1/369Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using gyroscopes as attitude sensors

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种三轴稳定卫星姿态控制,其特征在于:在一预定平面太阳敏感器范围,在与垂直方向测量范围局限于平面两侧夹角范围,并仅一个单轴测量、积分速度陀螺仪,其测量轴与太阳敏感器设置平面成一角度为±|β|,其角度为:|β|≥(π/2)-α2max。本发明还提供了使用上述有关测试设备,以确定卫星旋转速度的和捕获太阳的相对于由太阳光矢量给定的参考方向的位置误差的方法。

Description

卫星姿态控制的测试设备及其应用
本发明涉及一种三轴稳定卫星对太阳和地球定向的测试设备(具有太阳敏感器,以测定在卫星本体座标系XYZ中的太阳方向,太阳矢量S=(SX,SY,SZ)T,|S|=1,具有速度陀螺仪,以测定卫星的旋转速度矢量ω=(ωX,ωY,ωZ)T及其测定卫星旋转速度和卫星偏移的测试方法,以及使用这种测试设备时,执行姿态控制机动飞行的控制系统及控制方法。
有一种最早命名为这种类型的测试设备是出自H.Bittner在1977.11 ESA SP-128,1977.10.3~6,在Nordwijk举行的AOCS国际会议上公开的“The Attitude Determination and control Subsytemof the Inteelsat V Spacecraft”。这种公知类型的卫星具有一附加设置的、在卫星本体座标系三轴X、Y、Z测量的速度陀螺仪包以及同样附加设置的双视太阳敏感器,其中一个以负Y轴为中心,其覆盖XZ平面的一半,因而其垂直方向,即Y轴方向具有一定的宽度;另一个则以正X轴为中心,覆盖XY平面的
Figure C9211137900061
,同样在垂直方向上具有一定宽度,这种敏感器设备用于:卫星体的周转,也就是说,对卫星旋转速度矢量的分量进行跟踪测量,一旦太阳位于太阳敏感器的视场中,通过确定从卫星指向太阳中心点的太阳矢量即可独立进行卫星定向。这种由敏感器提供的测量结果将提供给卫星的姿态控制系统使用,从而达到改变其姿态的目的,也就是说发出旋转指令,使得其在转移轨道以及最终卫星轨道进行必要的机动飞行即能及时消除已经存在的姿态偏差。
用于测量、控制和衰减卫星周转运行使用的速度陀螺仪是一种结构非常复杂的电机精密仪器。因而相应造价昂贵并易受干扰,为了在宇宙空间特有的环境条件下满足其可靠性要求,尤其是负有长期使命时,这种仪器必须装设许多附加设备,从而使其成本和附加重量增加。如果在一个包含多个速度陀螺仪的陀螺包中,一根测量轴失效,一般而言,残留的其它陀螺仪轴也不能再使用。带有许多功能部件的系统换效的概率相应增加,带有大量易受干扰的电机部件的测量系统也存有效大的失效风险。
本发明的发明任务是:提从一种简单、造价低并且可靠的上述测试设备。这种测试设备能够使得:所有在转移轨道以及最终卫星轨道,尤其是在同步轨道中的卫星能自动地作必要的机动飞行此外,易受干扰的部件数量将尽可能减至最小。
必要的机动飞行包括有:从起始姿捕获太阳,从太阳定向搜地球,远地点机动飞行校准,当在转移轨道飞行期间在原规定方向定向和运动衰减中稳定姿态以及在同步轨道中对地球和太阳再搜索。
本发明提出的发明任务是通过以下方式来实现,即在一预定的平面内,太阳敏感器的测量范围负括周角α1,其中0≤α1≤2π,以及包括在与之垂直的方向上、在平面两侧的被限定的角度范围α2max≥α2≥-α2max2max<π/2),并且仅有一个单轴测量的积分速度陀螺仪,其测量轴G=(GX、GY、GZ)T,|G|=1,与设在太阳敏感器上的平面成一角度±|β|,其角度为:|β|≤(2/π)-α2max
该测试设备包括有多个并排列好的太阳敏感器,使得太阳矢量S的分量能在一平面在周角范围(0≤α1≤2π)测量,且在该平面的垂直方向能在一限定的角度范围|α2|≤α2max2max<90°)来进行测量。在此,以一个卫星体座标系为基础,特别是一种带Roll轴X,Nick轴Y以及Gier轴Z的笛卡儿XYZ系统,其中在同步轨道中一般而言Roll轴在轨道方向,Gier轴指向地球中心点和Nick轴垂直于轨道平面。下面提及的平面为测量平量。该平面应由太阳敏感器测量范捕捉住,例如XZ-平面。垂直方向,即Y轴方向,太阳敏感器测量范围则受到限制(±α2max)。太阳敏感器这样限定的测量范围能通过多个单一带有不同指向光轴的太阳敏感器相互毗连的、即重叠搭接的视场而构成。
另一实质性的特征是:以一个单轴测量、积分速度陀螺仪替代已公知的卫星三轴测量,速度陀螺仪包。其测量轴,通过一个单位矢量G来代表,并且与太阳敏感器测量平面成一角度β,角度大小为(π/2)-α2max。速度陀螺仪的测量轴也应位于从坐标系原点O为顶点的双锥体之内,双锥体的对称轴垂直于太阳敏感器的测量平面,如果以XZ平面作为太阳敏感器的测量平面,则以Y轴作为其对称轴。
关于本发明的进一步构成是:对速度陀螺仪测量轴允许的空间角度范围仍有限制。也即,带有测量轴 G T· S Ri≥ε>0。此处 S Ri为可预见的运转相位(1)的太阳参考矢量,在这些相位中,各相应的规定方向长时间地保持。
稍后,关于速度陀螺仪测量轴的可靠定向的原因变得更加清楚。后面给出的是测量方法,通过这些方法将能在一预先给定的参考方向测定卫星旋转速度和卫星偏移,当然是以依据本发明的测试设备为前提。后面又给出使用依据本发明的测试设备进行姿态控制机动飞行的一个控制系统以及其提到的相关使用方法、控制方法,尤其是给出了在所述控制系统中使用的控制方程式。
下面将根据附图更详细地描述本发明。附图以简图的方式给出:
图1:一颗地球卫星的转移轨道和同步轨道与地球和太阳之间的几何关系,
图2a、2b:太阳敏感器的视场,
图3:速度陀螺测量轴的允许空间角度范围,
图4:在使用依据本发明提出的测试设备时,其控制系统的方框图。
本发明可用于同步卫星,图1简要地表示了一颗这种卫星的转移轨道和同步卫星轨道的一部分。这颗卫星首先在合近地点(轨道上距地球最近的点)的椭圆形轨道上运行,其远地点(轨道上距地球最远的点)已经由基本为环形的确定的地球同步卫星轨道的半径所确定。同步卫星轨道的平面几乎与地球赤道平面一致。转移轨道平面则相对同步轨道的平面倾斜。
在这颗卫星最后进入地球同步轨道之前,它必须在与运载工具脱离之后在转移轨道上作一些机动飞行。在图1中表示了在转移轨道上的三个不同的位置P1~P3,卫星在这些位置上总是定向于不同的方向。在P1位置卫星的负Z轴指出太阳的方向,并以指令旋转速ωc(旋转速度矢量ωc)绕卫星——太阳中点连线转动。在P2位置,已经完成地球捕获(地球搜索),因此卫星的正Z轴向着地球中心定向。借助监测Z轴方向的红外线地球敏感器(IRS)形成上述定向,图中简略地表示了红外线地球敏感器相互垂直的测量方向。利用地球敏感器和太阳敏感器,上述定向使得在进入P3位置的远地点机动飞行之前建立三光轴参考系以完成陀螺定向和形成卫星正确定向成为可能。由于远地点发动机推力在负Z轴上的作用,对于这种机动飞行,卫星的正Z轴必须指向地球同步轨道与转移轨道的交点。
上述对地球和/或太阳的定向一结束,上述用于转移轨道的机动飞行(也就是太阳捕获和地球捕获)在地球同步轨道上将成为必不可少。很明显,本发明也可能用于其它的转移轨道构型和卫星构型,尤其是远地点发动机设备和地球敏感器。
在图2a中示出了如何校准太阳敏感器对称于XZ平面以及旋转对称于卫星本体坐标系中Y轴的视场也就是其测量范围,该测量范围包括一个相对KZ平面的角度范围α2max≥α2≥-α2max,这个角度范围必须相对Y轴旋转对称并因此形成一整个环视(0≤α1≤2π)。一个围绕Y轴旋转对称的空间角度范围不包括在太阳敏感器的测量范围内,所述空间角度范围的形状是一个张角为2×[(π/2-α2max]的在坐标系原点上对中的双锥。
在图2b中示出了如何在XZ平面中利用三个具有光轴OA1至OA3且相互相间120°设置的太阳敏感器来实现环视。在该实施例中,三个太阳敏感器中的任一个均在XZ平面内具有一个绕自身光轴对中的测量范围,该范围分别为2α11、2α12和2α13,例如各α为±67°。因此,在该具体情况下,在相邻的测量范围之间形成14°的重叠。在Y轴方向上,任一测量范围均为±α2max
很明显,各太阳敏感器处于另外的状况也是可能的,但在这些状况中总的来说必须在卫星本体坐标系的一个平面内(如XZ平面)以垂直于此平面的一个确定宽度提供上述环视。
也可以利用两个单轴测量太阳敏感器替代一个双轴测量太阳敏感器,各单轴太阳敏感器的测量方向绕各自的光轴转动90°即可。
图3a表示速度陀螺测量轴的允许空间角度范围。其中表示的是一个向卫星本体坐标系XYZ原点O对中的单位圆(半径为1)的侧视图。单位矢量G代表速度陀螺测量轴。单位圆上的单位矢量G的终点允许进入阴影线所示空间角度区域。该区域相对Y轴旋转对称并包含由XZ平面所测角β2的区域|β2|≥(π/2)-α2max。很显然,矢量G的位置并非如图3a所示限制于YZ平面中,而是首先允许其进入0≤α1≤2π的整个区域。
图3b以透视的方式示出了已经从图3a得知的实际情形,然而在此对速度陀螺(矢量G)测量轴的允许空间角度范围有另一个限制。通过在XZ平面处于X轴两测的角度范围±γx,并通过该角度绕Z轴的旋转来形成不许可的空间角度范围。图3b中单位圆上画阴影线所表示的四个区域保留下来作为速度陀螺测量轴的允许空间角度范围。
本发明的测试设备能够利用太阳敏感器确定太阳矢量S=(SX,SY,SZ)T,利用速度陀螺确定卫星旋转速度矢量ω的与其测量轴相平行的分量ωm,这些数值都是用于姿态控制的重要测量。值根据单位矢量(|S|=1)定义,太阳矢量S表示太阳位于卫星本体坐标系XZY的哪个方向上,速度陀螺测量轴由一个单位矢量G=(GX,GY,GZ)T表示(|G|=1)。卫星旋转速度矢量ω的分量ωm为:
       ωm=GT,ω
上标T表示一个矢量或一个矩阵的转置,因此矢量G被定义为列矢量,在上述标积中使用的矢量GT被定义为行矢量。两上矢量的标积中一个作为列矢量,另一个作为行矢量,这是众所周知的。
利用两个测量值S和ωm可确定用于姿态控制的两个重要的矢量,即卫星旋转速度矢量ω和相对一预先给定的参考姿态的卫星偏移矢量d。
卫星旋转速度矢量ω按以下公式计算得出:
        ω=ωpq    (1a) ω ‾ q = - S ‾ ~ · S ‾ . - - - ( 1 a )
        ωp=S(ST· ω) (1c) S ‾ T · ω ‾ = ω m - G ‾ T · ω ‾ q G ‾ T · S ‾ · - - - ( 1 d )
在这里,ωq、ωp分别表示与太阳矢量S垂直或平行的卫星转速矢量的分量,
Figure C9211137900133
表示一个变换表示的矢量积
Figure C9211137900134
的斜对称矩阵。 S ‾ ~ = 0 - S Z S Y S Z 0 - S X - S Y S X 0
ST·ω,GT·ωc,GT·S为相应矢量的标积,
Figure C9211137900136
为太阳矢量的时间导数,这个由太阳敏感器设备的测量信号(S)得出的导数可通过数值微分或通过(最好两极的)高通滤波器的传递函数得出: S ‾ ^ . = ( T 1 + T 2 ) S 2 + S ( 1 + T 1 S ) ( 1 + T 2 S ) · S ‾ - - - ( 1 e )
实际上,在按照方程(1a)至(1d)实施运算校验标量积GT·S数值时,当其数值超出了规定下限值(如:0.1),为了省去方程式(1d)和(1c)计算时估计数字困难,在方程式(1a)中矢量 ω p可简单地设为0。
卫星旋转速度矢量ω可以根据测定的太阳矢量S、其时间导数 和陀螺信号ωm及矩阵M按以下公式计算得出: ω ‾ = ( M ‾ T · M ) - 1 · M T S ‾ ^ . ω m , M ‾ = 0 - S Z S Y S Z 0 - S X - S Y S X 0 G X G Y G Z - - - ( 1 f ) 其也可由以下附加公式得出: S ‾ ~ G ‾ T · ω ‾ ≡ M ‾ · ω ‾ = ^ S ‾ . ω m - - - ( 1 g )
根据线性理论这样测定的ω估计值对旋转速度矢量而言会至少产生二次幂的误差,在实践中不好的条件矩阵M也会出现数值问题,例如,当太阳矢量S垂直位于速度陀螺仪的测量轴G时。如以上所述,开始提到的确定旋转速度矢量ω的方案是一种简单的方案,可以避免这类数值问题。
卫星偏移矢量d=(dX,dY,dZ)T表明了测定的太阳矢量S和一个太阳参数矢量SR之间的相互关系,矢量SR给出太阳矢量S的规定方向,当太阳矢量S不与太阳参考矢量 S R叠合时,即产生一个偏移量。通过卫星从参考方向旋转而产生一个这样
的偏移,在此分别用、θ和ψ表示卫星绕坐标轴X、Y、Z旋转的角度大小(、θ、ψ是欧拉(Euler)角)。由这种混合旋转描述的卫星偏移矢量d在小偏移时其分量可直接采用欧拉(Euler)角来表示:d=(、θ、ψ)T,依据测量值S和ωm,矢量d按下列公式计算确定:
d=dq+dp        (2a) d ‾ q = - S ‾ ~ R · ( S ‾ - S ‾ R ) ≡ - S ‾ ~ R · S - - - ( 2 b )
dp=SR(SR T·d)  (2c) S ‾ R T · d ‾ = ξ - G ‾ T · d ‾ q G ‾ T · S ‾ R - - - ( 2 d ) ξ ≡ G ‾ T · d ‾ = ∫ o t ω m dt - - - ( 2 e )
在这:dq、dp分别表示与太阳参考矢量SR(|SR|=1)垂直及平行的卫星偏移矢量d的分量,SR(类似于S)表示由SR矢量分量构成的斜对称矢量积矩阵,ξ表示速度陀螺仪的积分测量信号。
在这种情况下,按照方程式(2a)~(2e)实行运算时同样也应通过控制标量积GT·SR这样简便的方法来避免数值运算困难,如在确定旋转速度矢量ω时已描述过的一样。此外卫星偏移矢量d也可采用与计算旋转速度矢量ω类似的方法借助矩阵MR来确定,即按以下公式计算得出: d ‾ = ( M ‾ R T · M ‾ R ) · M ‾ R T S ‾ - S ‾ R ξ - - - ( 2 f )
以上计算公式是以小偏移的线性理论为前提,在这,还可采用一种近似的计算方法(E是单位矩阵): ( E ‾ - d ‾ ~ ) · S ‾ R = S ‾ S ‾ ~ R · d ‾ = S ‾ - S ‾ R ( = - d ‾ ~ · S ‾ R l )
   GT·d=ξ S ‾ ~ R G ‾ T · d ‾ ≡ M ‾ R · d ‾ = S ‾ - S ‾ R ξ - - - ( 2 g )
图4示出了一个三轴稳定地球卫星执行转移轨道和搜索机动飞行时依据本发明测试设备的控制系统方框图示例。卫星1的特征参数包括:旋转速度矢量ω和相对预先给定方向的太阳参考矢量SR的偏移量(Euler角、θ、ψ)(偏移矢量d)。只要太阳位于其可测范围,借助于太阳敏感器2即可测出与太阳参考矢量偏离的实际太阳矢量S,从而可以测出其偏移程度,并以一相应信号(Sunpresence=Sp)表示。一个单轴测量,积分速度陀螺仪3测出与测量轴平行的卫星旋转速度矢量ω的分量ωm。将测量值ωm和S按照公式(1a)至(1d)或(2a)至(2e)进行测量数据处理单元4,从而得出输出矢量值S,d和ω。一个逻辑单元5通过外来指令和内在的逻辑信号(如Sp)负责:在与之相连的控制器8中采用相应的控制原理;通过开关6、7将必要的参考值,即太阳参考矢量SR和指令旋转速度矢量ωc接合于系统中;在测量数据处理单元4中确定必须的矢量值。在控制器8的输出端给出一个推动调节执行装置9的矢量调节参数u,从而给出对应的卫星旋转力矩Tc
矢量值信号线路如图4双箭头所示,标量值或逻辑信号线路通过简单的实线示出。逻辑运算和数据处理公式和控制原理均能使用逻辑块和电子网络线路或卫星上的计算器来实现。
对于地球搜索和地球定向很重要的地球敏感器在图4中没有绘出。
执行转移轨道机动运行的控制原理以矢量描述的一般形式如下给出:
u=-KDLD(ω-ωC)+KpLP(SR·S)  (3a)
在此KP=diag(KPX,KPY,KPZ)和KD=diag(KDX,KDY,KDZ)为姿态放大系数对角矩阵和速度放大系数对角阵;Lp和LD为信号限制值;ωc为指令旋转速度矢量(如搜索运动时和太阳定向时)。姿态和速度放大系数考虑了调节执行装置(例如姿态控制喷管)的性质和配置,在绕不同的卫星轴旋转时大致会有不同的杠杆臂或不同的惯性力矩,但其目的是绕所述的轴产生相同的旋转力矩分量。信号限制值带有运算符号并使在方程式(3a)中给出的调节量矢量u的两个矢量分量按预先给定的限制不超过额度值,从而保证不产生围绕想要对准的规定方向的强烈摆动。
对绕飞行器各轴旋转的纯旋转速度进行调节时,例如在太阳定向时绕Z轴旋转(图1中P1位置),姿态放大系数矩阵(KPX,KPY,KPZ)的相应分量可设为0。
一种最重要的机动飞行随时需要进行,即当太阳不位于太阳敏感器的测量范围内并且因此完全确定旋转速度矢量ω已不可能是需要这种机动飞行。这种情况在卫星任何起始姿态或从规定方向作大的改变时容易出现。因而首先要采用控制原理
u=-KD·G·GT·ω  (3b)
来进行控制,即将卫星旋转速度ω的用单轴速度陀螺仪可测分量ωm=GT·ω调整为0。卫星仅能绕一根垂直于速度陀螺仪测量轴的轴进行旋转。
在此要描述本发明最重要的技术特征,这些特征确定涉及太阳敏感器测量范围(|α2|≤α2max)的速度陀螺仪测量轴的允许定向,即测量轴应位于以下范围|β2|≥(π/2)-α2max,也就是保证与测量轴垂直的方向位于太阳敏感器测量范围内。当卫星绕与测量轴垂直的方向旋转时,因为其旋转时卫星覆盖整个球面,太阳根据一定的旋转速度或早或迟进入太阳敏感器的测量范围。如果相反,速度陀螺仪测量轴不位于上面确定的允许范围内,则总是存在一绕与测量轴垂直的方向旋转对称的双球形空间范围,其不能落入太阳敏感器的测量范围。
如果使用方程式(3b)的控制原理引起的、不过不是公知的绕与测量轴垂直的轴的旋转速度太慢甚至根本不转,那么太阳在一给定的时间间隔内也照耀不到太阳敏感器的测量范围内,所以可采取相应控制措施。如图4通过一个逻辑信号简图所表示的一样,通过直接控制调节执行装置9给卫星施加一给定数信的旋转力矩脉冲从而获得一相应旋转速度(例如绕位于XZ平面的一旋转轴转动),太阳也因此进入其视场即太阳敏感器的测量范围中。
以上提及的和如图3b所描述的对速度陀螺仪测量轴定向允许范围的进一步限制,根据如下:图3b基于一种假定:在太阳矢量S与正Z轴互相重叠时,已预设了一个运转相位或该相位是可预见的。当速度陀螺仪的测量轴位于XY平面时,测量与XY平面垂直的旋转速度分量是不可能的。因为在这种情况下,速度陀螺仪只能在XY平面上测量旋转速度分量,并且太阳敏感器也不能测量与该平面垂直的旋转速度分量,因为此分量准确位于太阳矢量S方向即Z轴方向。尽管卫星旋转,太阳矢量仍保持其在坐标系中的准确定向。所以如图3所示给出了一个速度陀螺仪测量轴的禁区范围±γx。如果测量轴不位于禁区范围只要γx值不要选择得过小,速度陀螺仪就能测出Z方向的旋转速度分量。一般条件是:速度陀螺仪在太阳矢量S的方向的旋转速度分量应能够测量自不需说,测量轴代表矢量GT和太阳矢量S的标量积也不应为0,即必须有一确定最小值:
        GT·s≥ε>0
根据禁区角度范围±γx的不同的具体情形,给定值ε也相应不同,即ε=Sin(γx)。如果给定了一个工作状态,在此状态太阳矢量S在X轴方向定向。则速度陀螺仪测量轴则应禁止进入Z轴两侧的和旋转对称X轴的角度范围±γx。上述条件也适应于太阻矢量S的其它参考定向。
附图说明:
图1:转移轨道几何关系和飞行器定向
Sonne                                太阳
Neigungzekilptik                     黄道倾斜
Zielorbit                            目标轨道
quatorebene                        赤道平面
Transperbahn                         转移轨道
Erde                                 地球
Figure C9211137900192
quator                             赤道
Inklination der Transferbahn         转移轨道偏移
图4:执行转移轨道机动飞行的控制系统
Kdo’s                               指令
5-Logik                     逻辑单元
4-Messdatenaufbereitumg     测量数据处理单元
8-Regler                    控制器
9-Stellglieder              调节执行装置
1-Satellit                  卫星
3-Geschw.kreisel            速度陀螺仪
2-Sonnensensoren            太阳敏感器

Claims (9)

1.用于三轴稳定卫星姿态控制的测试设备,借助太阳放感器测定在一卫星本体座标系XYZ中太阳方向太阳矢量S=(SX,SY,SZ)T,借助速度陀螺仪测定卫星的旋转速度矢量ω=(ωX,ωY,ωZ)T,其特征在于:在一预定平面内包括周角α1,其中0≤α1≤2π,以及包括在与之垂直方向上在平面两侧的被限定的夹角范围α2max≥α2≥-α2max2max<π/2),并且仅有一个单轴测量的积分速度陀螺仪,其测量轴G=(GX、GY、GZ)T,|G|=1,与设在太阳敏感器上的平面成一角度±|β|,其角度为:
       |β|≤(2/π)-α2max
2.根据权利要求1所述的测试设备,其特征在于:带有测量轴G的速度陀螺仪如此确定方向,即总是满足以下附加条件:
       GT·SRi≥ε>0
在此SRi为可预见的以指数i表示运转相位的太阳参考矢量,在这些相位相应的规定方向长时间地保持。
3.使用如权利要求1所述的测试设备以确定卫星旋转速度的测定方法,其特征在于:借助太阳敏感器测定的太阳矢量S及由此推导出的其时间导数
Figure C9211137900021
,加之获取的速度陀螺仪的测量信号ωm,卫星旋转速度矢量ω即按以下公式计算: ω ‾ = ω ‾ q + S ‾ · ω ‾ m · G ‾ T · ω ‾ q G ‾ T · S ‾ ω ‾ q = · S ‾ ~ · S ‾ ^ . 其矩阵 S ‾ ~ ≡ 0 - S Z S Y S Z 0 - S X - S Y S X 0
4.使用据权利要求1所述的测试用以确定卫星旋转速度的测定方法,其特征在于:借助太阳敏感器测定的太阳矢量S和由此推导出的其时间导数 ,以及速度陀螺仪的测量信号ωm,卫星旋转速度矢量ω即按以下公式计算出: ω ‾ ≡ ( M ‾ T · M ) - 1 · M T · ^ S ‾ . ω m 其矩阵 M ≡ S ‾ ~ G ‾ T -
5.使用据权利要求1所述的测试设备用以相对太阳矢量S的一预先给定的参考方向,确定卫星偏移的测定方法,其特征在于:借助太阳敏感器测定的太阳矢量S和代表参考方向的太阳参考矢量SR以及速度陀螺仪的测量信号ωm,卫星偏移矢量d按以下公式算出: d ‾ = d ‾ q + S ‾ R ξ - G ‾ T · d ‾ q G ‾ T · S ‾ R d ‾ q = - S ‾ ~ R · S ‾ 其中 ξ ≡ ∫ o t ω m dt 其矩阵 S ‾ ~ R ≡ 0 - S RZ S RY S RZ 0 - S RX - S RY S RX 0
6.使用据权利要求1所述的测试设备用以在小偏移情况下,相对一个太阳矢量S的预先给出的参考方向确定卫星偏移的测定方法,其特征在于:从借助太阳敏感器测定的太阳矢量S和代表参考方向的太阳参考矢量SR以及速度陀螺仪的测量信号ωm,卫星偏移矢量d按以下公式计算出: d ‾ = ( M ‾ R T · M ‾ R ) · M ‾ R T S ‾ - S ‾ R ξ 其矩阵 M R = S ‾ ~ R G ‾ T
7.使用据权利要求1所述的测试设备用于机动控制三轴稳定卫星位置的控制系统,带有若干绕卫星本体座标轴(XYZ)施加调节力矩的调节执行装置和一个与这些执行装置前置联接的控制器,其特征在于一个测量数据处理装置(4)和一个开关装置(5,6,7),装置(4)用于按照由权利要求3或4之一给出的一种测定方法,确定卫星旋转速度矢量ω、按照由权利要求5或6给出的一种测定方法,确定卫星偏移矢量d以及确定太阳矢量S,装置(4)还给控制器提供相应的输入信号,开关装置(5,6,7)用于将太阳参考矢量SR和/或一个指令卫星放置速度矢量ωc送到控制器(8)的输入端。
8.据权利要求7所述在控制系统中采用的控制方法,其特征在于:控制器(8)以调节量矢量u的形式构成一个控制调节执行装置(9)的输出信号: u ‾ = - K ‾ D L D ( ω ‾ - ω ‾ C ) + K ‾ P L P ( S ‾ ~ R - S ‾ )
其中Kp和KD表示姿态及速度放大系数的对角矩阵,Lp和LD则表示信号限制值。
9.在一根据权利要求7所述的控制系统中使用的从卫星任一起始姿态开始捕获太阳的控制方法,在起始姿态太阳不位于太阳敏感器的可视范围,其特征在于:控制器(8)首先对调节执行装置(9)给出一个输出信号(矢量u)u=-KD·G·GT·ω,从而将卫星放置速度矢量ω的可由速度陀螺仪测量的分量GT·ω调节到0,此外,太阳没在一个预先给定期限内照射在太阳敏感器的视场时,调节执行装置(9)将会得到一个控制信号,通过此信号装置,(9)给出一个具有给定值的瞬时脉冲,以便卫星绕一根位于测量平面(XZ平面)内的并通过坐标系(XYZ)原点O的轴旋转,使太阳进入太阳敏感器的视场。
CN92111379A 1991-09-06 1992-09-05 卫星姿态控制的测试设备及其应用 Expired - Fee Related CN1039301C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DEP4129630.3 1991-09-06
DE4129630A DE4129630A1 (de) 1991-09-06 1991-09-06 Messanordnung und regelungssystem zur lageregelung eines dreiachsenstabilisierten satelliten sowie zugehoerige mess- und regelverfahren

Publications (2)

Publication Number Publication Date
CN1076419A CN1076419A (zh) 1993-09-22
CN1039301C true CN1039301C (zh) 1998-07-29

Family

ID=6439991

Family Applications (2)

Application Number Title Priority Date Filing Date
CN92111379A Expired - Fee Related CN1039301C (zh) 1991-09-06 1992-09-05 卫星姿态控制的测试设备及其应用
CN92111380A Expired - Fee Related CN1039302C (zh) 1991-09-06 1992-09-05 三轴稳定对地定向卫星的姿态控制系统和其搜索太阳和地球的方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN92111380A Expired - Fee Related CN1039302C (zh) 1991-09-06 1992-09-05 三轴稳定对地定向卫星的姿态控制系统和其搜索太阳和地球的方法

Country Status (7)

Country Link
US (2) US5535965A (zh)
EP (2) EP0601061B1 (zh)
JP (2) JP2637288B2 (zh)
CN (2) CN1039301C (zh)
CA (2) CA2117192C (zh)
DE (3) DE4129630A1 (zh)
WO (2) WO1993004922A1 (zh)

Families Citing this family (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19703629A1 (de) * 1997-01-31 1998-08-06 Daimler Benz Aerospace Ag Verfahren zur bordautonomen Bestimmung der Position eines Satelliten
US5597142A (en) * 1995-03-06 1997-01-28 Space Systems/Loral, Inc. Spacecraft acquisition of orientation by scan of earth sensor field of view
DE19510371C1 (de) * 1995-03-22 1996-10-31 Daimler Benz Aerospace Ag Verfahren zur Sonnensuche für einen dreiachsenstabilisierten Satelliten und dreiachsenstabilisierter Satellit
EP0795806B1 (de) * 1995-08-11 2001-12-05 Astrium GmbH Vorrichtung zur Lageregelung und Stabilisierung eines Raumfahrzeuges
DE19816978C1 (de) * 1998-04-17 1999-11-04 Daimler Chrysler Ag Verfahren zur Identifizierung eines fehlerhaft messenden Sensors in einem Raumfahrzeug
US6076774A (en) * 1998-08-12 2000-06-20 Hughes Electronics Corporation Fuel and thermal optimal spiral earth acquisition
US6145790A (en) * 1998-09-22 2000-11-14 Hughes Electronics Corporation Attitude determination system and method
US6152403A (en) * 1998-11-11 2000-11-28 Hughes Electronics Corporation Gyroscopic calibration methods for spacecraft
US6253125B1 (en) * 2000-03-01 2001-06-26 Space Systems/Loral, Inc. Method and apparatus for generating orbital data
US6317660B1 (en) * 2000-05-25 2001-11-13 Space Systems/Loral, Inc. Method for using satellite state vector prediction to provide satellite sensor automatic scan inhibit and/or sensor switching
CN100451898C (zh) * 2005-12-14 2009-01-14 上海微小卫星工程中心 微小卫星的姿态控制方法及系统
US7874519B2 (en) 2006-02-25 2011-01-25 Space Systems/Loral, Inc. Spacecraft three-axis attitude acquisition from sun direction measurement
ATE485220T1 (de) * 2007-07-30 2010-11-15 Astrium Gmbh Vorrichtung und verfahren zur bahnbestimmung und -vorhersage von satelliten, die signale an benutzer aussenden
FR2932163B1 (fr) * 2008-06-09 2010-06-11 Astrium Sas Procede de commande d'attitude de satellite et satellite commande en attitude
CN101509775B (zh) * 2009-03-18 2010-05-12 清华大学 一种多孔阵列式太阳敏感器的图像快速识别方法及装置
CN101712381B (zh) * 2009-11-13 2013-01-02 北京航空航天大学 一种基于多敏感器的定姿系统
CN101858969A (zh) * 2010-03-26 2010-10-13 航天东方红卫星有限公司 一种基于优化虚拟转动的卫星目标姿态预先确定方法
CN102176163B (zh) * 2010-12-10 2012-11-14 航天东方红卫星有限公司 一种任务观测持续时间的确定方法
CN102114918B (zh) * 2010-12-31 2013-03-20 北京航空航天大学 一种基于多速率敏感器组合定姿的姿控反馈回路
FR2980176A1 (fr) * 2011-09-19 2013-03-22 Astrium Sas Procede de controle d'attitude d'un satellite et satellite commande en attitude
CN102514734B (zh) * 2011-10-27 2013-11-27 北京航空航天大学 基于日地系统Halo轨道探测器构型与姿态指向的姿态递推方法
CN102681440A (zh) * 2012-05-02 2012-09-19 中国西安卫星测控中心 一种ns位保模式下脉冲修正方法
CN103264776B (zh) * 2013-05-30 2015-04-22 中国空间技术研究院 一种基于信息融合的控制系统工作模式设置及切换方法
CN103472846B (zh) * 2013-08-23 2015-10-21 北京控制工程研究所 一种基于多个太阳电池片的太阳捕获方法
FR3013685B1 (fr) * 2013-11-25 2017-05-19 Astrium Sas Procede et dispositif de commande d'une phase d'acquisition du soleil par un engin spatial
CN104071355B (zh) * 2014-06-12 2016-03-30 上海微小卫星工程中心 一种卫星姿态稳定控制方法及装置
CN107463191B (zh) * 2017-06-18 2020-01-07 珠海磐磊智能科技有限公司 控制力矩陀螺仪系统及行驶装置
CN108639385B (zh) * 2018-05-15 2021-04-13 浙江大学 一种无基准最简姿控系统快速稳定控制的实现方法
FR3093998B1 (fr) * 2019-03-20 2021-09-10 Airbus Defence & Space Sas Procédé de contrôle d’attitude d’un satellite en mode survie sans connaissance a priori de l’heure locale de l’orbite du satellite
CN111007865B (zh) * 2019-12-18 2023-07-14 中国人民解放军国防科技大学 以对日指向偏差为约束的卫星平稳对地定向方法
CN111268178B (zh) * 2020-03-16 2021-10-15 北京微动航科技术有限公司 卫星姿态控制飞轮对、卫星姿态控制系统和卫星
CN111483618B (zh) * 2020-04-09 2021-10-01 上海航天控制技术研究所 一种适用于火星探测的多指向驱动任务自主管理控制方法
CN112208797B (zh) * 2020-10-22 2022-07-01 上海卫星工程研究所 深空探测器二维天线电轴方向标定期间的姿态控制方法及系统
CN112977889B (zh) * 2021-03-19 2022-08-30 航天科工空间工程发展有限公司 一种基于太阳敏感器和地球敏感器的卫星姿态捕获方法
CN113386979B (zh) * 2021-06-03 2022-12-13 长光卫星技术股份有限公司 一种自适应规避太阳的数传姿态规划方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3937423A (en) * 1974-01-25 1976-02-10 Hughes Aircraft Company Nutation and roll error angle correction means
US4288051A (en) * 1976-09-18 1981-09-08 Messerschmitt-Bolkow-Blohm Gesellschaft Mit Beschrankter Haftung Method and apparatus for controlling a satellite
EP0394897A1 (fr) * 1989-04-24 1990-10-31 Alcatel Espace Procédé de mise à poste d'un satellite de télécommunications géostationnaire
US5042752A (en) * 1987-09-16 1991-08-27 Messerschmitt-Bolkow-Blohm Gmbh Apparatus for controlling the attitude of and stabilizing an elastic body

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE1039771B (de) * 1957-07-20 1958-09-25 Boelkow Entwicklungen Kg Mess- und Regeleinrichtung fuer sehr kleine Geschwindigkeiten
US3717413A (en) * 1969-11-20 1973-02-20 Tokyo Shibaura Electric Co Sun sensing system for a flying body
US3722297A (en) * 1970-11-19 1973-03-27 Sperry Rand Corp Fluid bearing gyroscope
US4071211A (en) * 1976-09-23 1978-01-31 Rca Corporation Momentum biased active three-axis satellite attitude control system
FR2397625A1 (fr) * 1977-07-11 1979-02-09 Radiotechnique Compelec Senseur d'orientation muni de deux cellules solaires
DE2749868C3 (de) * 1977-11-08 1980-05-22 Messerschmitt-Boelkow-Blohm Gmbh, 8000 Muenchen Sonnen- und Erderfassungsverfahren für Satelliten
DE2951125C2 (de) * 1979-12-19 1982-12-02 Messerschmitt-Bölkow-Blohm GmbH, 8000 München Verfahren zur Ermittlung der wahren Lotrichtung von Luft- und Raumfahrzeugen
FR2583873B1 (fr) * 1985-06-20 1987-09-11 Matra Procede et dispositif d'injection de satellite sur orbite geostationnaire avec stabilisation suivant les trois axes
DE3606636C1 (de) * 1986-02-28 1987-11-05 Messerschmitt Boelkow Blohm Verfahren zur Bestimmung von Erdmagnetfeldkomponenten bezueglich eines satellitenfesten Koordinatensystems
DE3634192A1 (de) * 1986-10-08 1988-04-21 Bodenseewerk Geraetetech Vorrichtung zur messung der rollrate oder rollage eines flugkoerpers
DE3734941A1 (de) * 1987-10-15 1989-04-27 Messerschmitt Boelkow Blohm Verfahren zur kalibrierung der kreisel eines dreiachsenstabilisierten satelliten
FR2637564B1 (fr) * 1988-10-06 1994-10-14 Aerospatiale Procede et systeme de controle autonome d'orbite d'un satellite geostationnaire
DE3842937A1 (de) * 1988-12-21 1990-07-05 Telefunken Systemtechnik Interferometer zur messtechnischen ermittlung von winkelgeschwindigkeiten im dreidimensionalen raum
JPH02274697A (ja) * 1989-04-14 1990-11-08 Toshiba Corp 宇宙航行体の姿勢制御装置
DE3922761C2 (de) * 1989-07-11 1994-08-25 Deutsche Aerospace Verfahren und Vorrichtung zum Ausrichten eines geostationären Satelliten in einer Referenzrichtung
US5020744A (en) * 1990-01-12 1991-06-04 General Electric Company Method for acquiring three-axis earth pointing attitude for an initially spinning spacecraft
US5080307A (en) * 1990-05-14 1992-01-14 Hughes Aircraft Company Spacecraft earth-pointing attitude acquisition method
US5130931A (en) * 1990-07-13 1992-07-14 General Electric Company Spacecraft attitude and velocity control system
DE4029215A1 (de) * 1990-09-14 1992-04-23 Deutsche Forsch Luft Raumfahrt Verfahren zur genauen messung raeumlicher winkel, trajektorien, konturen und bewegungsvorgaengen sowie schwereanomalien mit kreiseln und inertialsystemen
US5255879A (en) * 1991-11-27 1993-10-26 Hughes Aircraft Company Three axes stabilized spacecraft and method of sun acquisition

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3937423A (en) * 1974-01-25 1976-02-10 Hughes Aircraft Company Nutation and roll error angle correction means
US4288051A (en) * 1976-09-18 1981-09-08 Messerschmitt-Bolkow-Blohm Gesellschaft Mit Beschrankter Haftung Method and apparatus for controlling a satellite
US5042752A (en) * 1987-09-16 1991-08-27 Messerschmitt-Bolkow-Blohm Gmbh Apparatus for controlling the attitude of and stabilizing an elastic body
EP0394897A1 (fr) * 1989-04-24 1990-10-31 Alcatel Espace Procédé de mise à poste d'un satellite de télécommunications géostationnaire

Also Published As

Publication number Publication date
CA2117090C (en) 2002-02-19
CA2117192A1 (en) 1993-03-18
EP0601061A1 (de) 1994-06-15
JP2635821B2 (ja) 1997-07-30
CA2117192C (en) 2002-11-12
WO1993004923A1 (de) 1993-03-18
DE59202171D1 (de) 1995-06-14
WO1993004922A1 (de) 1993-03-18
DE59202172D1 (de) 1995-06-14
CA2117090A1 (en) 1993-03-18
US5535965A (en) 1996-07-16
US5558305A (en) 1996-09-24
CN1076419A (zh) 1993-09-22
CN1039302C (zh) 1998-07-29
EP0601061B1 (de) 1995-05-10
JP2637288B2 (ja) 1997-08-06
EP0601051A1 (de) 1994-06-15
JPH06510499A (ja) 1994-11-24
EP0601051B1 (de) 1995-05-10
JPH06510502A (ja) 1994-11-24
CN1071887A (zh) 1993-05-12
DE4129630A1 (de) 1993-05-06

Similar Documents

Publication Publication Date Title
CN1039301C (zh) 卫星姿态控制的测试设备及其应用
DE69630989T2 (de) Universelles Lagerregelungssystem für Raumfahrzeug
US7823836B2 (en) Optimal sun safe attitude for satellite ground tracking
Takada et al. Control moment gyro singularity-avoidance steering control based on singular-surface cost function
CA2040463C (en) Spacecraft earth-pointing attitude acquisition method
CN105905317A (zh) 一种卫星对日定向控制系统及其控制方法
DE69501728T2 (de) Künstlicher satellit mit aerodynamischen und magnetischen drehmomentgebern und steuerungsverfahren für einen derartigen satelliten
CN103389099B (zh) 基于x射线脉冲星的航天器姿态和位置测量系统及方法
CN111897357A (zh) 一种卫星对地扫描的姿态跟踪控制方法
CA1265225A (en) Pointing compensation system for spacecraft instruments
JPH0820269B2 (ja) 3軸安定宇宙船の天体探索および捕捉方法
US6347262B1 (en) Minimum fuel attitude and nutation controller for spinning spacecraft
CN110466803B (zh) 基于等倾角姿态控制的自旋稳定卫星姿态预测方法
DE3886714T2 (de) Nutationsensor und nutationsteuerungssystem für einen stabilisierten doppelspinsatellit.
US3171612A (en) Satellite attitude control mechanism and method
Karataş Leo satellites: dynamic modelling, simulations and some nonlinear attitude control techniques
RU2021173C1 (ru) Способ ориентации космического аппарата на планету
Hasegawa A study on rotation and its application for attitude reference system
DE3587199T2 (de) Kollimationskompensiersystem fuer raumfahrzeuginstrumente.
Wan et al. High precision satellite attitude control based on feedforward compensation
Celani Reduced-attitude stabilization for spacecraft boresight pointing using magnetorquers
Bolotskikh Analysis of Attitude Control System Flight Results of the Earth-Remote Sensing Nanosatellite OrbiCraft-Zorkiy Aleksey Bolotskikh а, Danil Ivanovb*, Stepan Tkachevc, Yaroslav Mashtakovd, Roman Zharkikhe
Mobley et al. MSX Attitude determination and control hardware
RU2021174C1 (ru) Способ ориентации космического аппарата по местной вертикали планеты
Davis Coordinate systems for the space shuttle program

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
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee