CN102424116B - 一种静止轨道卫星变轨策略优化方法 - Google Patents
一种静止轨道卫星变轨策略优化方法 Download PDFInfo
- Publication number
- CN102424116B CN102424116B CN2011104096283A CN201110409628A CN102424116B CN 102424116 B CN102424116 B CN 102424116B CN 2011104096283 A CN2011104096283 A CN 2011104096283A CN 201110409628 A CN201110409628 A CN 201110409628A CN 102424116 B CN102424116 B CN 102424116B
- Authority
- CN
- China
- Prior art keywords
- rail
- orbit
- satellite
- turns
- transfer orbit
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000001360 synchronised effect Effects 0.000 claims description 27
- 239000003380 propellant Substances 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000005484 gravity Effects 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 239000002699 waste material Substances 0.000 claims description 3
- 239000007788 liquid Substances 0.000 abstract 2
- 238000005070 sampling Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
Images
Landscapes
- Machines For Laying And Maintaining Railways (AREA)
Abstract
一种静止轨道卫星变轨策略优化方法,一是确定变轨次数、变轨圈次及各次变轨的控制量;二是确定各次变轨的时间和推力方向。目前静止轨道卫星的发射过程通常是,首先由运载火箭将卫星发射到一个带有倾角的大椭圆转移轨道,然后通过卫星自身所带的液体发动机进行若干次远/近地点变轨,转移至地球同步轨道,最后通过轨道倾角修正和圆化实现静止轨道。从转移轨道经过若干次远/近地点变轨实现静止轨道对卫星而言操作比较复杂,因此变轨次数一般不宜设置过多,以免增加变轨的复杂性和风险;另一方面,考虑到卫星的液体发动机能力以及变轨期间的弧段损失等因素,变轨次数也不宜太少。
Description
技术领域
本发明涉及一种静止轨道卫星变轨策略优化方法。
背景技术
静止轨道卫星变轨策略设计属于轨道机动类问题的一种,在数学上可采用极大值极小值方法来求解,或者划为非线性规划类问题,采取相应的算法求解。根据实际情况,非线性规划问题有多种算法,一般可以采用穷举法,它的特点是简单、但是比较费时、计算量大的一种算法,它的缺点是需要采样点足够密集才能找到比较精确的全局最优解,而且当函数的性质不够清楚的时候难以定出合理的采样范围和采样密度,另外当需要优化的参数增多时,计算量会迅速增大。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种静止轨道卫星变轨策略优化方法,合理确定变轨策略设计的各个约束条件,以减少变轨策略设计过程中的人工干预以及计算时间和计算量。
本发明的技术解决方案是:
一种静止轨道卫星变轨策略优化方法,步骤如下:
(1)确定卫星从当前轨道出发到目标轨道所需要的变轨次数n,具体通过如下步骤进行:
(1.1)通过公式计算从当前轨道以脉冲方式机动至目标轨道所需总速度增量dv,其中,v0和v1变轨前后卫星在远地点的速度,且 其中,μ为地心引力常数,ra为卫星处于远地点的地心距离;a0,a1表示卫星所处当前轨道的半长轴;δ表示两轨道面夹角;
(1.2)通过公式dm=m0(1-exp(-dv/Isp/g))计算从当前轨道机动至目标轨道所需推进剂质量dm,其中,m0为起飞重量,Isp为发动机比冲,g为重力常数;
(1.3)通过公式dt=dm/(F/Isp/g)计算从当前轨道机动至目标轨道发动机所需工作时间dt,其中,F表示发动机推力;
(1.4)根据步骤(1.3)中得到的从当前轨道机动至目标轨道发动机所需工作时间dt估算从目标轨道机动至目标轨道所需变轨次数n,在满足变轨次数尽可能少的前提下,使发动机平均工作时间控制在指标Tengine之内,指标Tengine可由发动机单次点火时间长度限制值除以1.5得到;
(2)确定各次变轨的发动机熄火条件ai及相邻两次变轨之间的漂移圈数qi,所述熄火条件ai即为轨道半长轴的大小,具体通过如下步骤进行:
(2.1)根据预定的入轨参数确定转移轨道类型为同步转移轨道或者超同步转移轨道;
(2.2)对于同步转移轨道,确定首次变轨和末次变轨的熄火条件ai及漂移圈数qi,同步转移轨道变轨包括若干次远地点变轨,超同步转移轨道变轨包括若干次远地点变轨和一次近地点变轨;对于超同步转移轨道,确定首次变轨和最后两次变轨的熄火条件及漂移圈数;
对于同步转移轨道或者超同步转移轨道,根据漂移圈数确定原则,由λ0,可确定首次远地点变轨的漂移圈数q1,λ0,分别为同步转移轨道和超同步转移轨道在星箭分离时刻卫星的星下点地理经度和经度漂移率,λ0,均根据转移轨道参数计算,是通过转移轨道类型可直接确定的数据;
所述漂移圈数确定原则为:两次变轨之间卫星漂移圈数最少、每次变轨均在地面测控范围内进行以及最终实现目标轨道;
同步转移轨道和超同步转移轨道卫星变轨首次变轨的熄火条件的确定原则是:在发动机单次工作时间范围内使半长轴尽可能大,即使得熄火条件尽可能大;
对于超同步转移轨道末次远地点变轨,即倒数第二次变轨,有:λ′obj=λobj+20°,发动机熄火条件an′由超同步转移轨道参数计算得到,根据漂移圈数确定原则qn′为1圈;
对于超同步转移轨道,末次变轨为近地点变轨:变轨的星下点经度为卫星定点轨位,经度漂移率为0度/圈,末次变轨的发动机熄火条件an为42164km,根据漂移圈数确定原则qn为0.5圈;
(2.3)确定其余的变轨熄火条件及漂移圈数;
在约束条件下,根据λ0、q1、λn′和确定余下的变轨熄火条件及漂移圈数;λn′是指第n′次变轨的星下点经度,是指第n′次变轨结束后的经度漂移率,所述约束条件是指是两次变轨之间圈数最少、变轨点在地面测控范围内以及发动机单次工作时间不超自身限制;
(3)确定每次变轨的点火方向δi和点火时刻ti,具体步骤如下:
(3.1)以点火方向δi和点火时刻ti为优化参数,定义优化目标函数F(X)为F(X)=Δmi(ti,δi)+ΔmR,其中,Δmi表示第i次变轨所需推进剂,ΔmR表示从第i次变轨后的轨道出发到达目标轨道所需推进剂,
(3.2)以公式F(X)=Δmi(ti,δi)+ΔmR所计算的推进剂消耗量F(X)最小为原则确定各点火方向和点火时刻;
(4)根据步骤(1)中得到的变轨次数、步骤(2)中得到的各次变轨的发动机熄火条件、相邻两次变轨之间的漂移圈数、步骤(3)中得到的每次变轨的点火方向和点火时刻计算n次变轨后所实现的目标轨道,并计算其轨位,根据其与所需定点轨位的差值Δλ,对第2次变轨熄火条件进行调整。
所述从第i次变轨后的轨道出发到达目标轨道所需推进剂ΔmR通过如下方式计算:
(2.1)通过公式计算从当前轨道以脉冲方式机动至目标轨道所需速度增量dv,其中,v0和v1变轨前后卫星在远地点的速度,且 μ为地心引力常数,ra,rp分别为卫星处于远地点、近地点时与地心之间的距离;a0,a1表示卫星所处当前轨道的半长轴;δ表示两轨道面夹角;
(2.2)通过公式ΔmR=m0(1-exp(-dv/Isp/g))计算从当前轨道机动至目标轨道所需推进剂质量ΔmR,其中,m0为起飞重量,Isp为发动机比冲,g为重力常数。
所述步骤(2.2)中对于超同步转移轨道末次远地点变轨,发动机熄火条件an′由超同步转移轨道参数计算得到,具体过程为: 其中为地球赤道半径,h0为同步轨道高度,ha为转移轨道远地点高度,a,e表示转移轨道半长轴和偏心率。
所述步骤(2.3)中在约束条件下,根据λ0、a1、q1、λn′和确定余下的变轨熄火条件及漂移圈数,具体过程为:由λ0、q1可确定λ1,由λ1、a1可确定λ2,根据a1与an′的差值大小并结合漂移圈数确定原则来确定a2,由λ2、a2可确定q2,依此类推,确定第一次变轨和最后一次远地点变轨之间的各次变轨熄火条件和漂移圈数。
所述步骤(4)中对第2次变轨熄火条件进行调整通过如下方式进行:
本发明与现有技术相比的有益效果是:
本发明方法在现有静止轨道卫星变轨策略设计方法的基础上,重新分析变轨设计的各项限制因素,合理确定优化过程中需要考虑的约束条件,计算量小,计算时间短,自动化程度高,人工干预少。
附图说明
图1为本发明方法流程图;
具体实施方式
如图1所示,本发明提供了一种静止轨道卫星变轨策略优化方法,步骤如下:
(1)确定卫星从当前轨道出发到目标轨道所需要的变轨次数n,具体通过如下步骤进行:
(1.1)通过公式计算从当前轨道以脉冲方式机动至目标轨道所需总速度增量dv,其中,v0和v1变轨前后卫星在远地点的速度,且 其中,μ为地心引力常数,ra为卫星处于远地点的地心距离;a0,a1表示卫星所处当前轨道的半长轴;δ表示两轨道面夹角;
(1.2)通过公式dm=m0(1-exp(-dv/Isp/g))计算从当前轨道机动至目标轨道所需推进剂质量dm,其中,m0为起飞重量,Isp为发动机比冲,g为重力常数;
(1.3)通过公式dt=dm/(F/Isp/g)计算从当前轨道机动至目标轨道发动机所需工作时间dt,其中,F表示发动机推力;
(1.4)根据步骤(1.3)中得到的从当前轨道机动至目标轨道发动机所需工作时间dt估算从目标轨道机动至目标轨道所需变轨次数n,在满足变轨次数尽可能少的前提下,使发动机平均工作时间控制在指标Tengine之内,指标Tengine可由发动机单次点火时间长度限制值除以1.5得到;
(2)确定各次变轨的发动机熄火条件ai及相邻两次变轨之间的漂移圈数qi,所述熄火条件ai即为轨道半长轴的大小,具体通过如下步骤进行:
(2.1)根据预定的入轨参数确定转移轨道类型为同步转移轨道或者超同步转移轨道;
(2.2)对于同步转移轨道,确定首次变轨和末次变轨的熄火条件ai及漂移圈数qi,同步转移轨道变轨包括若干次远地点变轨,超同步转移轨道变轨包括若干次远地点变轨和一次近地点变轨;对于超同步转移轨道,确定首次变轨和最后两次变轨的熄火条件ai及漂移圈数qi;
对于同步转移轨道或者超同步转移轨道,根据漂移圈数确定原则,由λ0,可确定首次远地点变轨的漂移圈数q0,λ0,分别为同步转移轨道和超同步转移轨道在星箭分离时刻卫星的星下点地理经度和经度漂移率,λ0,均根据转移轨道参数计算,是通过转移轨道类型可直接确定的数据;
所述漂移圈数确定原则为:两次变轨之间卫星漂移圈数最少、每次变轨均在地面测控范围内进行以及最终实现目标轨道;
同步转移轨道和超同步转移轨道卫星变轨首次变轨的熄火条件的确定原则是:在发动机单次工作时间范围内使半长轴尽可能大,即使熄火条件尽可能大;
对于超同步转移轨道末次远地点变轨,即倒数第二次变轨,有:λ′obj=λobj+20°,发动机熄火条件an′由超同步转移轨道参数计算得到,根据漂移圈数确定原则qn′为1圈;
对于超同步转移轨道,末次变轨为近地点变轨:变轨的星下点经度为卫星定点轨位,经度漂移率为0度/圈,末次变轨的发动机熄火条件an为42164km,根据漂移圈数确定原则qn为0.5圈;
(2.3)确定其余的变轨熄火条件及漂移圈数;
在约束条件下,根据λ0、q0、λn′和确定余下的变轨熄火条件及漂移圈数;λn′是指第n′次变轨的星下点经度,是指第n′次变轨结束后的经度漂移率,所述约束条件是指是两次变轨之间圈数最少、变轨点在地面测控范围内以及发动机单次工作时间不超自身限制;
(3)确定每次变轨的点火方向δi和点火时刻ti,具体步骤如下:
(3.1)以点火方向δi和点火时刻ti为优化参数,定义优化目标函数F(X)为F(X)=Δmi(ti,δi)+ΔmR,其中,Δmi表示第i次变轨所需推进剂,ΔmR表示从第i次变轨后的轨道出发到达目标轨道所需推进剂,
(3.2)以公式F(X)=Δmi(ti,δi)+ΔmR所计算的推进剂消耗量F(X)最小为原则确定各点火方向和点火时刻。其中Δmi通过公式计算,Δt表示本次点火时间长度,Δt由各次变轨的熄火条件ai和备选的发动机点火时刻ti、点火方向δi决定,本次点火时间长度Δt通过求解卫星动力学方程得到;ΔmR的计算方法同(1)部分的①②;
所述从第i次变轨后的轨道出发到达目标轨道所需推进剂ΔmR通过如下方式计算:
(2.1)通过公式计算从当前轨道以脉冲方式机动至目标轨道所需速度增量dv,其中,v0和v1变轨前后卫星在远地点的速度,且 μ为地心引力常数,ra,rp分别为卫星处于远地点、近地点时与地心之间的距离;a0,a1表示卫星所处当前轨道的半长轴;δ表示两轨道面夹角;
(2.2)通过公式ΔmR=m0(1-exp(-dv/Isp/g))计算从当前轨道机动至目标轨道所需推进剂质量ΔmR,其中,m0为起飞重量,Isp为发动机比冲,g为重力常数。
所述步骤(2.2)中对于超同步转移轨道末次远地点变轨,发动机熄火条件an′由超同步转移轨道参数计算得到,具体为: 其中为地球赤道半径,h0为同步轨道高度,ha为转移轨道远地点高度,a,e表示转移轨道半长轴和偏心率。
所述步骤(2.3)中在约束条件下,根据λ0、q0、λn′和确定余下的变轨熄火条件及漂移圈数具体通过如下方式进行:根据λ0、a1、q1、λn′和确定余下的变轨熄火条件及漂移圈数,具体过程为:由λ0、q1可确定λ1,由λ1、a1可确定λ2,根据a1与an′的差值大小并结合a1的确定原则来确定a2(a1的确定原则即为漂移圈数确定原则),由λ2、a2可确定q2。依此类推,确定第一次变轨和最后一次远地点变轨之间的各次变轨熄火条件和漂移圈数。
(4)根据步骤(1)中得到的变轨次数、步骤(2)中得到的各次变轨的发动机熄火条件、相邻两次变轨之间的漂移圈数、步骤(3)中得到的每次变轨的点火方向和点火时刻计算n次变轨后所实现的目标轨道,并计算其轨位,根据其与所需定点轨位的差值Δλ,对第2次变轨熄火条件进行调整,使得Δλ等于0,即实现卫星最终定点轨位。
对第2次变轨熄火条件进行调整通过如下方式进行:
最后将T代入中计算调整后第二次变轨的熄火条件a2。
Claims (2)
1.一种静止轨道卫星变轨策略优化方法,其特征在于步骤如下:
(1)确定卫星从当前轨道出发到目标轨道所需要的变轨次数n,具体通过如下步骤进行:
(1.1)通过公式计算从当前轨道以脉冲方式机动至目标轨道所需总速度增量dv,其中,v0和v1变轨前后卫星在远地点的速度,且其中,μ为地心引力常数,ra为卫星处于远地点的地心距离;a0,a1表示卫星所处当前轨道的半长轴;δ表示两轨道面夹角;
(1.2)通过公式dm=m0(1-exp(-dv/Isp/g))计算从当前轨道机动至目标轨道所需推进剂质量dm,其中,m0为起飞重量,Isp为发动机比冲,g为重力常数;
(1.3)通过公式dt=dm/(F/Isp/g)计算从当前轨道机动至目标轨道发动机所需工作时间dt,其中,F表示发动机推力;
(1.4)根据步骤(1.3)中得到的从当前轨道机动至目标轨道发动机所需工作时间dt估算从当前轨道机动至目标轨道所需变轨次数n,在满足变轨次数尽可能少的前提下,使发动机平均工作时间控制在指标Tengine之内,指标Tengine可由发动机单次点火时间长度限制值除以1.5得到;
(2)确定各次变轨的发动机熄火条件ai及相邻两次变轨之间的漂移圈数qi,所述熄火条件ai即为轨道半长轴的大小,具体通过如下步骤进行:
(2.1)根据预定的入轨参数确定转移轨道类型为同步转移轨道或者超同步转移轨道;
(2.2)对于同步转移轨道,确定首次变轨和末次变轨的熄火条件ai及漂移圈数qi,同步转移轨道变轨包括若干次远地点变轨,超同步转移轨道变轨包括若干次远地点变轨和一次近地点变轨;对于超同步转移轨道,确定首次变轨和最后两次变轨的熄火条件及漂移圈数;
对于同步转移轨道或者超同步转移轨道,根据漂移圈数确定原则,由λ0,可确定首次远地点变轨的漂移圈数q1,λ0,分别为同步转移轨道和超同步转移轨道在星箭分离时刻卫星的星下点地理经度和经度漂移率,λ0,均根据转移轨道参数计算,是通过转移轨道类型可直接确定的数据;
所述漂移圈数确定原则为:两次变轨之间卫星漂移圈数最少、每次变轨均在地面测控范围内进行以及最终实现目标轨道;
同步转移轨道和超同步转移轨道卫星变轨首次变轨的熄火条件的确定原则是:在发动机单次工作时间范围内使半长轴尽可能大,即使得熄火条件尽可能大;
对于超同步转移轨道末次远地点变轨,即倒数第二次变轨,有:λ′obj=λobj+20°,发动机熄火条件an′由超同步转移轨道参数计算得到,根据漂移圈数确定原则qn′为1圈;
对于超同步转移轨道,末次变轨为近地点变轨:变轨的星下点经度为卫星定点轨位,经度漂移率为0度/圈,末次变轨的发动机熄火条件an为42164km,根据漂移圈数确定原则qn为0.5圈;
对于超同步转移轨道末次远地点变轨,发动机熄火条件an′由超同步转移轨道参数计算得到,具体过程为: 其中为地球赤道半径,h0为同步轨道高度,ha为转移轨道远地点高度,a,e表示转移轨道半长轴和偏心率;
(2.3)确定其余的变轨熄火条件及漂移圈数;
在约束条件下,根据λ0、q1、λn′和确定余下的变轨熄火条件及漂移圈数;λn′是指第n′次变轨的星下点经度,是指第n′次变轨结束后的经度漂移率,所述约束条件是指是两次变轨之间圈数最少、变轨点在地面测控范围内以及发动机单次工作时间不超自身限制;
所述在约束条件下,根据λ0、a1、q1、λn′和确定余下的变轨熄火条件及漂移圈数,具体过程为:由λ0、q1可确定λ1,由λ1、a1可确定λ2,根据a1与an′的差值大小并结合漂移圈数确定原则来确定a2,由λ2、a2可确定q2,依此类推,确定第一次变轨和最后一次远地点变轨之间的各次变轨熄火条件和漂移圈数;
(3)确定每次变轨的点火方向δi和点火时刻ti,具体步骤如下:.
(3.1)以点火方向δi和点火时刻ti为优化参数,定义优化目标函数F(X)为F(X)=Δmi(ti,δi)+ΔmR,其中,Δmi表示第i次变轨所需推进剂,ΔmR表示从第i次变轨后的轨道出发到达目标轨道所需推进剂,
(3.2)以公式F(X)=Δmi(ti,δi)+ΔmR所计算的推进剂消耗量F(X)最小为原则确定各点火方向和点火时刻;
(4)根据步骤(1)中得到的变轨次数、步骤(2)中得到的各次变轨的发动机熄火条件、相邻两次变轨之间的漂移圈数、步骤(3)中得到的每次变轨的点火方向和点火时刻计算n次变轨后所实现的目标轨道,并计算其轨位,根据其与所需定点轨位的差值Δλ,对第2次变轨熄火条件进行调整;
所述对第2次变轨熄火条件进行调整通过如下方式进行:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011104096283A CN102424116B (zh) | 2011-12-08 | 2011-12-08 | 一种静止轨道卫星变轨策略优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011104096283A CN102424116B (zh) | 2011-12-08 | 2011-12-08 | 一种静止轨道卫星变轨策略优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102424116A CN102424116A (zh) | 2012-04-25 |
CN102424116B true CN102424116B (zh) | 2013-11-20 |
Family
ID=45958180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011104096283A Expired - Fee Related CN102424116B (zh) | 2011-12-08 | 2011-12-08 | 一种静止轨道卫星变轨策略优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102424116B (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880184B (zh) * | 2012-10-24 | 2015-05-27 | 北京控制工程研究所 | 一种静止轨道卫星自主轨道控制方法 |
CN103144784B (zh) * | 2013-02-26 | 2015-12-23 | 北京空间飞行器总体设计部 | 一种适用于导航geo卫星的离轨方法 |
CN103303496B (zh) * | 2013-06-28 | 2015-09-09 | 上海新跃仪表厂 | 一种微幅调整轨道高度的方法 |
CN103412563B (zh) * | 2013-07-25 | 2015-09-23 | 北京控制工程研究所 | 一种预估轨控综合效率和推进剂消耗量的方法 |
CN103593505B (zh) * | 2013-10-16 | 2016-06-01 | 中国空间技术研究院 | 一种卫星轨控发动机安装参数确定方法 |
CN103885456B (zh) * | 2014-03-31 | 2019-02-15 | 航天东方红卫星有限公司 | 一箭多星异轨面分布转移方法 |
CN104015938B (zh) * | 2014-05-07 | 2016-05-25 | 北京理工大学 | 一种电推进静止轨道卫星的位置保持方法 |
CN104443432B (zh) * | 2014-11-25 | 2016-06-15 | 哈尔滨工业大学 | 一种卫星有限推力共面圆轨道自主轨道转移制导方法 |
CN104714554B (zh) * | 2015-03-17 | 2017-07-28 | 北京航空航天大学 | 一种共位式卫星编队的构造方法 |
CN105607478B (zh) * | 2016-01-21 | 2017-07-07 | 北京理工大学 | 地球静止轨道航天器电推进转移轨道控制方法 |
CN106156414B (zh) * | 2016-06-30 | 2019-04-26 | 北京润科通用技术有限公司 | 一种卫星轨道仿真方法及装置 |
CN106570316B (zh) * | 2016-10-20 | 2018-12-21 | 北京空间飞行器总体设计部 | 基于推进剂预算的低轨椭圆轨道卫星成功入轨判定方法 |
CN107554820B (zh) * | 2017-09-04 | 2019-08-20 | 上海航天控制技术研究所 | 一种基于增益在线估计的卫星自主轨道保持方法 |
CN107826269B (zh) * | 2017-09-18 | 2019-10-22 | 北京控制工程研究所 | 一种适用于地球静止轨道卫星平台的近地点变轨方法 |
CN107885917A (zh) * | 2017-10-27 | 2018-04-06 | 中国地质大学(武汉) | 基于异面变轨策略的卫星星座重构方法、设备及存储设备 |
CN108216687B (zh) * | 2017-12-25 | 2019-12-20 | 中国空间技术研究院 | 基于粒子群算法的geo卫星变轨策略计算方法、***及介质 |
CN108408083B (zh) * | 2018-03-02 | 2019-02-26 | 北京空间技术研制试验中心 | 载人航天器在轨运行风险防控方法 |
CN109625323B (zh) * | 2018-11-09 | 2021-07-20 | 中国科学院空间应用工程与技术中心 | 一种卫星化学推进变轨方法及*** |
CN109460049B (zh) * | 2018-11-14 | 2021-11-16 | 北京控制工程研究所 | 基于惯性指向模式的地球同步轨道卫星远地点变轨方法 |
CN111338367B (zh) * | 2020-02-27 | 2022-10-04 | 中国西安卫星测控中心 | 一种偏心率冻结同轨双脉冲控制的中间轨道确定方法 |
CN113602531B (zh) * | 2021-07-05 | 2023-05-12 | 上海卫星工程研究所 | 异常分离情况下组合体geo变轨策略生成方法及*** |
CN113525721B (zh) * | 2021-07-29 | 2023-10-03 | 陕西星邑空间技术有限公司 | 卫星轨道变换方法、装置、电子设备及存储介质 |
CN113968360B (zh) * | 2021-08-09 | 2024-03-15 | 中国空间技术研究院 | 一种静止轨道卫星星上自主电推进变轨方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2747102B1 (fr) * | 1996-04-05 | 1998-06-26 | Europ Propulsion | Procede et systeme de mise en orbite d'un vehicule spatial avec des propulseurs a forte impulsion specifique |
FR2768399B1 (fr) * | 1997-09-12 | 1999-11-12 | Aerospatiale | Procede de mise sur orbite geostationnaire d'un satellite artificiel |
CN101214860A (zh) * | 2007-12-26 | 2008-07-09 | 北京控制工程研究所 | 轨控过程中自主选取定姿方式的方法 |
CN100493993C (zh) * | 2007-12-26 | 2009-06-03 | 北京控制工程研究所 | 一种卫星的自主变轨方法 |
-
2011
- 2011-12-08 CN CN2011104096283A patent/CN102424116B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102424116A (zh) | 2012-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102424116B (zh) | 一种静止轨道卫星变轨策略优化方法 | |
Bandikova et al. | GRACE accelerometer data transplant | |
CN109625323B (zh) | 一种卫星化学推进变轨方法及*** | |
US8489258B2 (en) | Propulsive guidance for atmospheric skip entry trajectories | |
CA2948860C (en) | Orbit transfer method for a spacecraft using a continuous or quasi-continuous thrust and embedded driving system for implementing such a method | |
RU2737644C2 (ru) | Энергоэффективное маневрирование спутника | |
US11753188B2 (en) | Optimized power balanced low thrust transfer orbits utilizing split thruster execution | |
CN103257653A (zh) | 一种基于燃料消耗优化的卫星编队构形控制方法 | |
Mazal et al. | Closed-loop distance-keeping for long-term satellite cluster flight | |
US11396388B2 (en) | Optimized power balanced variable thrust transfer orbits to minimize an electric orbit raising duration | |
US6845950B1 (en) | System for high efficiency spacecraft orbit transfer | |
CN110032768B (zh) | 一种使用精确动力学模型的四脉冲轨道交会优化方法 | |
CN115639580A (zh) | 卫星自主导航方法 | |
CN103678814B (zh) | 临界倾角近圆轨道的偏心率预偏置设计方法 | |
CN101186236A (zh) | 一种减少航天器重力损失的变轨方法 | |
DK2586711T3 (en) | Method and system for control of a unit of at least two satellites adapted to provide a service | |
CN107589756A (zh) | 一种奔月卫星编队初始化方法 | |
Wertz et al. | Autonomous constellation maintenance | |
O'Shaughnessy et al. | MESSENGER's use of solar sailing for cost and risk reduction | |
Shang et al. | Design and optimization of low-thrust orbital phasing maneuver | |
Gill et al. | First results from a hardware-in-the-loop demonstration of closed-loop autonomous formation flying | |
Leomanni et al. | Precise attitude control of all-electric GEO spacecraft using Xenon microthrusters | |
Herman | Balancing, turning, saving special AOCS operations to extend the GRACE mission | |
RU2596004C2 (ru) | Способ управления движением космического аппарата на активном участке его выведения на орбиту искусственного спутника планеты | |
Nugnes et al. | Multi-revolution low-thrust trajectory optimisation using differential dynamic programming in orbital element formulation |
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: 20131120 |