CN114394260B - 地球同步卫星小推力南北保持控制解析计算方法 - Google Patents
地球同步卫星小推力南北保持控制解析计算方法 Download PDFInfo
- Publication number
- CN114394260B CN114394260B CN202111473573.2A CN202111473573A CN114394260B CN 114394260 B CN114394260 B CN 114394260B CN 202111473573 A CN202111473573 A CN 202111473573A CN 114394260 B CN114394260 B CN 114394260B
- Authority
- CN
- China
- Prior art keywords
- long
- ignition
- control
- north
- south
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 26
- 238000012423 maintenance Methods 0.000 title claims abstract description 20
- 230000007774 longterm Effects 0.000 claims abstract description 74
- 238000000034 method Methods 0.000 claims abstract description 12
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 16
- 230000001174 ascending effect Effects 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000003203 everyday effect Effects 0.000 claims description 3
- 239000000446 fuel Substances 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000013213 extrapolation Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开的地球同步卫星小推力南北保持控制解析计算方法,具体计算过程包括如下步骤:步骤1、解析计算地球同步卫星倾角矢量的长期摄动变化量;步骤2、计算电推力器点火序列,解析计算出整个规划周期内每一次电推点火的开关机时刻。本发明实现地球同步卫星小推力南北保持点火序列的计算,通过解析方程精确计算出同步卫星倾角矢量的长期项摄动变化量,点火序列计算简便精准,且点火时长接近理论最优,具有计算速度快、适用性强的优点。
Description
技术领域
本发明属于航天导航技术领域,具体涉及地球同步卫星小推力南北保持控制解析计算方法。
背景技术
通常来说,为满足星上的通信、光学等载荷的地面指向性或者覆盖性需求,地球同步卫星的倾角矢量要始终处于指定的较小范围内。在地球非球形引力和日月引力摄动的影响下,卫星的倾角矢量会逐渐漂移,为了克服这种自然摄动产生的倾角矢量漂移,而定期进行的倾角矢量控制,就是工程上俗称的“南北保持”。
在长期管理阶段,地球同步卫星的南北保持是最耗燃料的卫星控制工作,仅抵消一年的轨道倾角长期项摄动就需要给卫星提供约43m/s的速度增量,消耗数十千克燃料质量。随着小推力技术逐渐成熟,越来越多的地球静止卫星采用小推力推力器进行南北控制,小推力推力器最重要的优点是比冲大,提供相同速度增量,仅需要消耗化学推进器10%左右的燃料质量。然而,目前的小推力推力器推力一般只有几十到几百毫牛,实现南北保持需求的倾角控制量需要长时间开机(每天数小时)。此外,小推力推力器工质的动能来源是星上电能(化学推力器是燃料化学能),过长的开机时间对卫星能源供应是个挑战。因此,根据卫星应用约束合理设计电推力器南北控制策略,使卫星以尽可能少的开机时间,对于完成南北保持控制是很有意义的。
传统的化学推进南北控制是在倾角即将超出限制区域时,选取合适的时刻利用法,向速度增量将倾角矢量向年平均摄动方向的反向控制,只在倾角矢量ix方向到达边界限制时,对控制方向做一些调整。此策略的优点是,基本控制方向一年才做一次调整,便于操作,且控制方向也近似是燃料最优控制方向。年平均摄动方向的计算通常有两种方法——解析公式或轨道外推,解析公式计算会不可避免的存在一定误差,而轨道外推计算不能提出长周期项的影响。因此,如果在小推力条件下长期固定倾角矢量控制方向,方向误差造成的累积效应会使倾角矢量与控制目标显著偏差。
发明内容
本发明的目的是提供地球同步卫星小推力南北保持控制解析计算方法,解决了现有装配了电推力器的地球同步卫星南北保持任务,解析公式计算存在误差的问题。
本发明所采用的技术方案是,地球同步卫星小推力南北保持控制解析计算方法,具体计算过程包括如下步骤:
步骤1、解析计算地球同步卫星倾角矢量的长期摄动变化量;
步骤2、计算电推力器点火序列,解析计算出整个规划周期内每一次电推点火的开关机时刻。
本发明的其他特点还在于,
步骤1的具体步骤为:确定倾角矢量为(ixiy),其中ix=icosΩ,iy=isinΩ,i和Ω分别为地球同步卫星轨道要素的倾角和升交点赤经;计算规划时段t0到tf内,倾角矢量的长期摄动变化量(Δix)长期和(Δiy)长期。
步骤1具体按照以下步骤:
步骤1.1、计算地球非球形引力对倾角矢量的摄动量;
步骤1.2、计算太阳引力对倾角矢量的摄动量;
步骤1.3、计算月球引力对倾角矢量的摄动量;
步骤1.4、将步骤1.1至步骤1.3分别求得的地球非球形引力、太阳引力、月球引力对倾角矢量的摄动量求和,即为地球同步卫星倾角矢量的长期项摄动变化量。
步骤1.1中地球非球形引力对倾角矢量的摄动量计算公式为:
其中,n为卫星轨道角速度,a为卫星轨道半长轴,μe为地球引力常数,Re为地球半径,J2、J22、J31、J33、J4、J43、λ22、λ31、λ33、λ42、λ43是地球非球形四阶主项摄动系数,ix、iy为倾角矢量初值。
步骤1.2中太阳引力对倾角矢量的摄动量计算公式为:
其中,μs为太阳引力常数,rs为日地距离,Ωs为地心惯性系下的太阳轨道升交点赤经,is为太阳轨道倾角。
步骤1.3中月球引力对倾角矢量的摄动量计算公式为:
其中,μl为月球引力常数,rl为地月距离,Ωl为月球轨道升交点赤经,il为月球轨道倾角。
步骤1.4中地球同步卫星倾角矢量的长期项摄动变化量由上述三个部分组成,即
(Δix)长期=(Δix)非球形,长期+(Δix)太阳,长期+(Δix)月球,长期 (7)
(Δiy)长期=(Δiy)非球形,长期+(Δiy)太阳,长期+(Δiy)月球,长期 (8)。
步骤2的具体过程如下:
步骤2.1、根据步骤1中地球同步卫星倾角矢量的长期项摄动变化量得到一个周期南北保持的倾角矢量目标控制量;一个周期南北保持的倾角矢量目标控制量能够抵消周期内的长期项摄动量,即
(Δiy)control=-(Δiy)长期 (9)
(Δix)control=-(Δix)长期 (10);
步骤2.2、计算每次小推力点火时长,得到规划周期内的倾角控制量;
根据步骤2.1中一个周期南北保持的倾角矢量目标控制量得到控制方向角,控制方向角αi=arctan2[(Δiy)control,(Δix)control],南侧电推力器的点火中心赤经为InnerTwoPi(αi),北侧电推力器的点火中心赤经为InnerTwoPi(αi+π);
点火序列规划周期为N天,其中的M天进行电推点火(M≤N,为轨道确定、载荷控制预留空间),每天有2次南北控制的时机,电推点火总次数为2M,计算每次小推力点火时长为
Tfire=2arcsin[(Δi)controlVsωe/acc/M/4]/ωe (11)
其中,(Δi)control为规划周期内的倾角控制量,acc为电推力器法向合加速度,Vs为地球同步卫星标称轨道速度,ωe为地球自转角速度;
电推力器一天内的两次点火时刻t1,on、t1,off和t2,on、t2,off分别由南侧和北侧点火中心赤经和点火时长计算出来,设当天00:00:00.00时刻的格林尼治恒星时角为Θ0,卫星平经度为λ;
t1,on=InnerTwoPi(αi-λ-Θ0)/ωe-Tfire/2 (13)
t1,off=t1,on+Tfire (14)
t2,on=InnerTwoPi(αi+π-λ-Θ0)/ωe-Tfire/2 (15)
t2,off=t2,on+Tfire (16)
计算出来的t1,on、t1,off和t2,on、t2,off是从当天00:00:00.00时刻开始的积秒,以此类推,解析计算出整个规划周期内每一次电推点火的开关机时刻。
本发明的有益效果是,本发明地球同步卫星小推力南北保持控制计算方法,适应于小推力条件下的地球同步卫星南北位置保持控制机动规划,针对装配了电推力器的地球同步卫星南北保持任务,提出一种解析计算电推力器点火序列的方法,控制卫星抵消长期项倾角矢量摄动,使卫星倾角矢量在小范围内做周期性运动,从而实现地球同步卫星小推力南北保持点火序列的计算,通过解析方程精确计算出同步卫星倾角矢量的长期项摄动变化量,点火序列计算简便精准,且点火时长接近理论最优,提高了计算速度,满足快速规划计算的需求,具备广泛的适用性。
附图说明
图1是本发明一种实施例连续进行一年小推力南北控制的卫星倾角矢量变化情况图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的地球同步卫星小推力南北保持控制解析计算方法,具体包括如下步骤:
步骤1、解析计算地球同步卫星倾角矢量的长期摄动变化量;
步骤1的具体步骤为:确定倾角矢量为(ixiy),其中ix=icosΩ,iy=isinΩ,i和Ω分别为地球同步卫星轨道要素的倾角和升交点赤经;计算规划时段t0到tf内,倾角矢量的长期摄动变化量(Δix)长期和(Δiy)长期;
步骤1.1、计算地球非球形引力对倾角矢量的摄动量;
地球非球形引力对倾角矢量的摄动量计算公式为:
其中,n为卫星轨道角速度,a为卫星轨道半长轴,μe为地球引力常数,Re为地球半径,J2、J22、J31、J33、J4、J43、λ22、λ31、λ33、λ42、λ43是地球非球形四阶主项摄动系数,ix、iy为倾角矢量初值。
步骤1.2、计算太阳引力对倾角矢量的摄动量;
太阳引力对倾角矢量的摄动量计算公式为:
其中,μs为太阳引力常数,rs为日地距离,Ωs为地心惯性系下的太阳轨道升交点赤经,is为太阳轨道倾角。
步骤1.3、计算月球引力对倾角矢量的摄动量;
月球引力对倾角矢量的摄动量计算公式为:
其中,μl为月球引力常数,rl为地月距离,Ωl为月球轨道升交点赤经,il为月球轨道倾角。
步骤1.4、将步骤1.1至步骤1.3分别求得的地球非球形引力、太阳引力、月球引力对倾角矢量的摄动量求和,即为地球同步卫星倾角矢量的长期项摄动变化量。
地球同步卫星倾角矢量的长期项摄动变化量由上述三个部分组成,即
(Δix)长期=(Δix)非球形,长期+(Δix)太阳,长期+(Δix)月球,长期 (7)
(Δiy)长期=(Δiy)非球形,长期+(Δiy)太阳,长期+(Δiy)月球,长期 (8)
步骤2、计算电推力器点火序列,解析计算出整个规划周期内每一次电推点火的开关机时刻。
电推力器能够提供的推力只有几十/几百毫牛,完成南北保持任务需要平均每天开机数小时,所以工程上通常采取一次性计算包含多天多次控制的电推力器点火序列的方式进行。
步骤2.1、根据步骤1中地球同步卫星倾角矢量的长期项摄动变化量得到一个周期南北保持的倾角矢量目标控制量,一个周期南北保持的倾角矢量目标控制量能够抵消周期内的长期项摄动量,即
(Δiy)control=-(Δiy)长期 (9)
(Δix)control=-(Δix)长期 (10)步骤2.2、计算每次小推力点火时长,得到规划周期内的倾角控制量;
根据步骤2.1中一个周期南北保持的倾角矢量目标控制量得到控制方向角,控制方向角αi=arctan2[(Δiy)control,(Δix)control],南侧电推力器的点火中心赤经为InnerTwoPi(αi),北侧电推力器的点火中心赤经为InnerTwoPi(αi+π)。
点火序列规划周期为N天,其中的M天进行电推点火(M≤N,为轨道确定、载荷控制预留空间),每天有2次南北控制的时机,电推点火总次数为2M,计算每次小推力点火时长为
Tfire=2arcsin[(Δi)controlVsωe/acc/M/4]/ωe (11)
其中,(Δi)control为规划周期内的倾角控制量,acc为电推力器法向合加速度,Vs为地球同步卫星标称轨道速度,ωe为地球自转角速度。
电推力器一天内的两次点火时刻t1,on、t1,off和t2,on、t2,off分别由南侧和北侧点火中心赤经和点火时长计算出来。设当天00:00:00.00时刻的格林尼治恒星时角为Θ0,卫星平经度为λ。
t1,on=InnerTwoPi(αi-λ-Θ0)/ωe-Tfire/2 (13)
t1,off=t1,on+Tfire (14)
t2,on=InnerTwoPi(αi+π-λ-Θ0)/ωe-Tfire/2 (15)
t2,off=t2,on+Tfire (16)
计算出来的t1,on、t1,off和t2,on、t2,off是从当天00:00:00.00时刻开始的积秒,以此类推,解析计算出整个规划周期内每一次电推点火的开关机时刻。
通过上述方法,即可实现地球同步卫星小推力南北保持点火序列的计算,其特点是通过解析方程精确计算出同步卫星倾角矢量的长期项摄动变化量,点火序列计算简便精准,且点火时长接近理论最优。
实施例
本实施例以亚太6D卫星一个周期的小推力南北控制序列计算为例,如图1所示,南北控制周期为7天(2021年6月7日至2021年6月13日),其中前5天进行电推点火。
利用式(1)至式(8)计算7天的倾角矢量长期摄动变化量。
带入常数λ22=-14.929°、λ31=6.968°、λ33=20.994°、λ44=30.280°、J22=1.811528e-6、J31=2.2091169e-6、J33=0.2213602e-6、J44=7.63937899e-9、Re=6378.1363、μ=398600.4415和a=42165.700、n=7.292116e-5、ix=0.02°、iy=0.03°、tf-t0=86400.0×7。
带入常数μs=1.3271222e11、rs=1.496e8、Ωs=0.0°、is=23.4437°;
带入常数μl=4.90280107e3、rl=3.85e5(月球倾角il和升交点赤经Ωl是时变量,由当前日期计算);
(Δix)长期=(Δix)非球形,长期+(Δix)太阳,长期+(Δix)月球,长期 (7)
(Δiy)长期=(Δiy)非球形,长期+(Δiy)太阳,长期+(Δiy)月球,长期 (8)代数式(1)至式(6)计算结果
(Δix)长期=-0.0023°
(Δiy)长期=0.0169°
由式(9)和(10)计算倾角矢量目标控制量
(Δiy)control=-(Δiy)长期 (9)
(Δix)control=-(Δix)长期 (10)
(Δiy)control=0.0023°
(Δix)control=-0.0169°
控制方向角αi=arctan2[(Δiy)control,(Δix)control]=277.8198°,南侧电推力器的点火中心赤经为InnerTwoPi(αi),北侧电推力器的点火中心赤经为
InnerTwoPi(αi+π)。
N=7,M=5
带入卫星电推加速度acc,计算每次小推力点火时长为
Tfire=2arcsin[(Δi)controlVsωe/acc/M/4]/ωe (11)
Tfire=10785(s)
规划周期内的倾角控制量
(Δi)control=0.171°
由式(13)至式(16)计算周期内每次小推力点火的开关机时刻
t1,on=InnerTwoPi(αi-λ-Θ0)/ωe-Tfire/2 (13)
t1,off=t1,on+Tfire (14)
t2,on=InnerTwoPi(αi+π-λ-Θ0)/ωe-Tfire/2 (15)
t2,off=t2,on+Tfire (16)
例如,对2021年6月7日的两次点火有:
Θ0=261.5293°,亚太6D卫星的定点经度为134.0°E,带入式(13)至式(16)得到t1,on=26787.5(s),t1,off=37572.5(s),t2,on=69987.5(s),t2,off=80772.5(s),对应时刻如下表1所示。
表1亚太6D卫星对应时刻数据表
开机时刻t1,on | 开机时刻t1,off | 开机时刻t2,on | 开机时刻t2,off |
07:26:26.5 | 09:57:01.5 | 19:26:26.5 | 21:57:01.5 |
Claims (2)
1.地球同步卫星小推力南北保持控制解析计算方法,其特征在于,具体计算过程包括如下步骤:
步骤1、解析计算地球同步卫星倾角矢量的长期摄动变化量;
步骤2、计算电推力器点火序列,解析计算出整个规划周期内每一次电推点火的开关机时刻;
步骤1的具体步骤为:确定倾角矢量为(ix iy),其中ix=icosΩ,iy=isinΩ,i和Ω分别为地球同步卫星轨道要素的倾角和升交点赤经;计算规划时段t0到tf内,倾角矢量的长期摄动变化量(Δix)长期和(Δiy)长期;
所述步骤1具体按照以下步骤:
步骤1.1、计算地球非球形引力对倾角矢量的摄动量;
步骤1.2、计算太阳引力对倾角矢量的摄动量;
步骤1.3、计算月球引力对倾角矢量的摄动量;
步骤1.4、将步骤1.1至步骤1.3分别求得的地球非球形引力、太阳引力、月球引力对倾角矢量的摄动量求和,即为地球同步卫星倾角矢量的长期项摄动变化量;
所述步骤1.1中地球非球形引力对倾角矢量的摄动量计算公式为:
其中,n为卫星轨道角速度,a为卫星轨道半长轴,μe为地球引力常数,Re为地球半径,J2、J22、J31、J33、J4、J43、λ22、λ31、λ33、λ43是地球非球形四阶主项摄动系数,ix、iy为倾角矢量初值;
所述步骤1.2中太阳引力对倾角矢量的摄动量计算公式为:
其中,μs为太阳引力常数,rs为日地距离,Ωs为地心惯性系下的太阳轨道升交点赤经,is为太阳轨道倾角;
所述步骤1.3中月球引力对倾角矢量的摄动量计算公式为:
其中,μl为月球引力常数,rl为地月距离,Ωl为月球轨道升交点赤经,il为月球轨道倾角;
所述步骤1.4中地球同步卫星倾角矢量的长期项摄动变化量由上述三个部分组成,即
(Δix)长期=(Δix)非球形,长期+(Δix)太阳,长期+(Δix)月球,长期 (7)
(Δiy)长期=(Δiy)非球形,长期+(Δiy)太阳,长期+(Δiy)月球,长期 (8)。
2.如权利要求1所述的地球同步卫星小推力南北保持控制解析计算方法,其特征在于,所述步骤2的具体过程如下:
步骤2.1、根据步骤1中地球同步卫星倾角矢量的长期项摄动变化量得到一个周期南北保持的倾角矢量目标控制量;一个周期南北保持的倾角矢量目标控制量能够抵消周期内的长期项摄动量,即
(Δiy)control=-(Δiy)长期 (9)
(Δix)control=-(Δix)长期 (10);
步骤2.2、计算每次小推力点火时长,得到规划周期内的倾角控制量;
根据步骤2.1中一个周期南北保持的倾角矢量目标控制量得到控制方向角,控制方向角αi=arctan2[(Δiy)control,(Δix)control],南侧电推力器的点火中心赤经为InnerTwoPi(αi),北侧电推力器的点火中心赤经为InnerTwoPi(αi+π);
点火序列规划周期为N天,其中的M天进行电推点火,M≤N,为轨道确定、载荷控制预留空间,每天有2次南北控制的时机,电推点火总次数为2M,计算每次小推力点火时长为
Tfire=2arcsin[(Δi)controlVsωe/acc/M/4]/ωe (11)
其中,(Δi)control为规划周期内的倾角控制量,acc为电推力器法向合加速度,Vs为地球同步卫星标称轨道速度,ωe为地球自转角速度;
电推力器一天内的两次点火时刻t1,on、t1,off和t2,on、t2,off分别由南侧和北侧点火中心赤经和点火时长计算出来,设当天00:00:00.00时刻的格林尼治恒星时角为Θ0,卫星平经度为λ;
t1,on=InnerTwoPi(αi-λ-Θ0)/ωe-Tfire/2 (13)
t1,off=t1,on+Tfire (14)
t2,on=InnerTwoPi(αi+π-λ-Θ0)/ωe-Tfire/2 (15)
t2,off=t2,on+Tfire (16)
计算出来的t1,on、t1,off和t2,on、t2,off是从当天00:00:00.00时刻开始的积秒,以此类推,解析计算出整个规划周期内每一次电推点火的开关机时刻。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111473573.2A CN114394260B (zh) | 2021-11-29 | 2021-11-29 | 地球同步卫星小推力南北保持控制解析计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111473573.2A CN114394260B (zh) | 2021-11-29 | 2021-11-29 | 地球同步卫星小推力南北保持控制解析计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114394260A CN114394260A (zh) | 2022-04-26 |
CN114394260B true CN114394260B (zh) | 2023-10-27 |
Family
ID=81225229
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111473573.2A Active CN114394260B (zh) | 2021-11-29 | 2021-11-29 | 地球同步卫星小推力南北保持控制解析计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114394260B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115626303A (zh) * | 2022-09-27 | 2023-01-20 | 上海航天控制技术研究所 | 一种采用电推的地球同步轨道倾角保持方法 |
CN116946392B (zh) * | 2023-07-12 | 2024-03-05 | 中国西安卫星测控中心 | 基于多维姿态偏置的地球同步卫星电推进倾角控制方法 |
CN117163325B (zh) * | 2023-11-02 | 2024-01-02 | 北京控制工程研究所 | 考虑故障容错的多冗余姿控发动机分配方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104015938A (zh) * | 2014-05-07 | 2014-09-03 | 北京理工大学 | 一种电推进静止轨道卫星的位置保持方法 |
CN105197257A (zh) * | 2015-08-31 | 2015-12-30 | 中国空间技术研究院 | 一种分舱优化设计的桁架式geo卫星推力器布局方法 |
CN110254753A (zh) * | 2019-06-04 | 2019-09-20 | 北京理工大学 | 一种地球静止轨道卫星电推力器及其布局优化方法 |
CN113148236A (zh) * | 2021-04-25 | 2021-07-23 | 上海宇航系统工程研究所 | 一种基于霍尔推力器的南北位保方法 |
-
2021
- 2021-11-29 CN CN202111473573.2A patent/CN114394260B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104015938A (zh) * | 2014-05-07 | 2014-09-03 | 北京理工大学 | 一种电推进静止轨道卫星的位置保持方法 |
CN105197257A (zh) * | 2015-08-31 | 2015-12-30 | 中国空间技术研究院 | 一种分舱优化设计的桁架式geo卫星推力器布局方法 |
CN110254753A (zh) * | 2019-06-04 | 2019-09-20 | 北京理工大学 | 一种地球静止轨道卫星电推力器及其布局优化方法 |
CN113148236A (zh) * | 2021-04-25 | 2021-07-23 | 上海宇航系统工程研究所 | 一种基于霍尔推力器的南北位保方法 |
Non-Patent Citations (1)
Title |
---|
蒯政中.脉冲和电推进下的地球静止卫星位置保持和轨道转移策略研究.《中国优秀硕士学位论文全文数据库(电子期刊)工程科技II辑》.2018,(第1期),第C031-449页. * |
Also Published As
Publication number | Publication date |
---|---|
CN114394260A (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114394260B (zh) | 地球同步卫星小推力南北保持控制解析计算方法 | |
Foster et al. | Constellation phasing with differential drag on planet labs satellites | |
CN108490963B (zh) | 全电推进卫星电推力器故障模式下的位置保持方法及系统 | |
CN111268177B (zh) | 一种静止轨道卫星分散式闭环自主位置保持控制方法 | |
CN109656133B (zh) | 一种针对空间走廊跟踪观测的分布式卫星群优化设计方法 | |
CN110068845B (zh) | 一种基于平根数理论确定卫星理论轨道的方法 | |
CN110429974A (zh) | 基于回归轨道星座的快速对准方法和装置 | |
CN115015983B (zh) | 光学遥感星座的长期在轨维持方法 | |
Kos et al. | Altair descent and ascent reference trajectory design and initial dispersion analyses | |
Park et al. | Nanosatellite constellation deployment using on-board magnetic torquer interaction with space plasma | |
Jesick et al. | Navigation overview for the Mars atmosphere and volatile evolution mission | |
Golikov | THEONA—a numerical-analytical theory of motion of artificial satellites of celestial bodies | |
Shah | Automated station-keeping for satellite constellations | |
Oh et al. | Analytic Optimization of Mixed Chemical-Electric Orbit Raising Missions | |
Emma et al. | Algorithm for autonomous longitude and eccentricity control for geostationary spacecraft | |
Stanbridge et al. | Achievable force model accuracies for messenger in mercury orbit | |
RU2535353C2 (ru) | Способ удержания космического аппарата на геосинхронной 24-часовой орбите | |
Boone et al. | Resolution of Orbit Determination prediction instabilities at Titan during Cassini’s Solstice mission | |
Yazdi et al. | Analysis of parking orbits and transfer trajectories for mission design of cis-lunar space stations | |
COOK et al. | Return to the moon-The Lunar Observer mission | |
Chao | Semiautonomous stationkeeping of geosynchronous satellites | |
Ryne et al. | GRAIL orbit determination for the science phase and extended mission | |
McAdams et al. | MESSENGER–Six Primary Maneuvers, Six Planetary Flybys, and 6.6 Years to Mercury Orbit | |
Zhao et al. | Precise Intelligent Orbit Control Strategy on Geostationary Orbit Satellite | |
Manohar et al. | Trajectory reconstruction during thrusting phase of rockets using differential corrections |
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 |