CN112464429B - 一种小推力航天器轨道根数长期演化的极大值估计方法 - Google Patents
一种小推力航天器轨道根数长期演化的极大值估计方法 Download PDFInfo
- Publication number
- CN112464429B CN112464429B CN202011072998.8A CN202011072998A CN112464429B CN 112464429 B CN112464429 B CN 112464429B CN 202011072998 A CN202011072998 A CN 202011072998A CN 112464429 B CN112464429 B CN 112464429B
- Authority
- CN
- China
- Prior art keywords
- maximum value
- period
- variation
- track
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明提出一种小推力航天器轨道根数长期演化的极大值估计方法,属于航天器动力学与控制领域。本发明首先分别建立关于半长轴、偏心率、轨道倾角、升交点赤经、近拱点幅角、真近点角六个轨道根数的摄动方程,然后将前五个方程分别除以第六个方程,得到前五个轨道根数分别关于真近点角的微分方程;基于微分代数原理,利用轨道动力学特性,分析推导出小推力航天器轨道根数关于轨道周期的长期演化变化量极大值估计方法。本发明能够快速计算出航天器轨道根数变化量极大值,计算消耗低,适用范围广,可提高航天器轨道根数变化量极大值估计的计算效率与适用范围。
Description
技术领域
本发明属于航天器动力学与控制领域,特别涉及一种小推力航天器轨道根数长期演化的极大值估计方法。
背景技术
离子电推进等推进技术的发展带动了航天任务的设计应用的拓展。连续小推力推进以其高推进效率在航天器在轨服务、航天器轨道保持和编队交会等航天任务中有着广泛应用。
对于航天器长期在轨运行的任务周期,以及小推力轨迹优化等问题,通常需要进行轨道递推演化来确定航天器轨道变化范围以及确定变轨能力。如何以简便快速的方式进行小推力航天器长期轨道根数演化能够为轨迹优化及航天任务设计提供便利和参考。
当前小推力航天器轨道根数递推通常采用两种方式:一种是数值积分的方式,利用积分器对航天器从初始轨道进行积分,对于长轨道积分时间的情况以及轨迹优化等情况需要很大的计算量以及计算时间消耗;另一种半解析估计方式,一般是利用简化的小推力形式,只考虑单一确定方向的小推力推进,进行单一小推力方向的小推力航天器轨道根数估计,只能用于确定在此种小推力作用下,轨道根数的变化情况,无法确定小推力航天器轨道根数的变化范围边界。
发明内容
本发明旨在解决当前小推力航天器长期轨道演化轨道根数变化极大值估计时所面临的计算量大,适用范围窄的不足之处,提出一种小推力航天器轨道根数长期演化的极大值估计方法。本发明计算消耗低,适用范围广,可提高航天器轨道根数变化量极大值估计的计算效率与适用范围。
本发明提出一种小推力航天器轨道根数长期演化的极大值估计方法,其特征在于,包括以下步骤:
1)建立航天器摄动方程:
对于高斯形式的拉格朗日行星运动方程,航天器轨道根数a,e,i,Ω,ω,f分别代表半长轴、偏心率、轨道倾角、升交点赤经、近拱点幅角、真近点角;
令T=[av,an,ah]为小推力航天器的摄动加速度,其中av代表摄动加速度在密切轨道平面内沿速度方向的分量,an代表摄动加速度在垂直速度方向的分量,ah代表摄动加速度在垂直于轨道平面的分量;
建立航天器摄动方程表达式如下:
其中,θ=ω+f表示真近点角与近地点幅角之和;r,v,h,μ分别为卫星地心距、速率、角动量标量和中心天体引力常数;t表示时间变量;
2)摄动方程转换;
对于步骤1)得到的摄动方程式(1)-(6),将式(1)-(5)分别除以式(6),得到轨道根数a,e,i,Ω,ω关于真近点角f的微分方程,从而得到以真近点角为自变量的高斯摄动方程式(7)-(11),表达式如下:
3)计算航天器轨道根数长期变化量的极大值;具体步骤如下:
3-1)设定半长轴、偏心率、轨道倾角、升交点赤经、近拱点幅角和真近点角的初值分别为a0,e0,i0,Ω0,ω0,f0,设定表示使得轨道根数在一个周期变化量取极大值的最优摄动加速度,其中三个分量分别代表最优摄动加速度在密切轨道平面内沿速度方向、垂直速度方向以及垂直于轨道平面的分量;
3-2)计算半长轴、偏心率、近拱点幅角长期变化量的极大值;具体步骤如下:
3-2-1)设定初始周期数N=1,则对于半长轴、偏心率、近拱点幅角存在:aN-1=a0,eN-1=e0,ωN-1=ω0;其中,aN-1为第N-1个周期的半长轴,eN-1为第N-1个周期的偏心率,ωN-1为第N-1个周期的近拱点幅角;
3-2-2)从半长轴、偏心率、近拱点幅角中任意选取一个轨道根数;
3-2-3)计算步骤3-2-2)选取的轨道根数在第N-1个周期的变化量极大值;方法如下:
3-2-3-1)对于半长轴;
3-2-3-2)对于偏心率;
其中,||T||表示加速度矢量的模;AN-1,BN-1为用于简化加速度表示形式所引入的中间量;
3-2-3-3)对于近拱点幅角;
使得近拱点幅角在第N-1个周期变化率最大的最优摄动加速度的三个分量分别为:
其中,Aω,N-1,Bω,N-1,Cω,N-1为中间量,计算表达式分别如下:
3-2-4)计算步骤3-2-2)选取的轨道根数在第N个周期的值为:
3-2-5)令N=N+1,然后重新返回步骤3-2-3),直至N的取值到达设定的数值时,得到该轨道根数在第N-1个周期变化的极大值;
3-3)计算轨道倾角、升交点赤经长期变化量的极大值;
从轨道倾角和升交点赤经中任意选取一个轨道根数,计算该轨道根数在第N个周期的变化量极大值,方法如下:
3-3-1)对于轨道倾角:
使得轨道倾角一个周期变化量取得极大值的最优摄动加速度在密切轨道平面内垂直速度方向分量为:
将式(18)的结果代入式(9)里积分,得到轨道倾角在一个周期的变化量极大值为:
其中H,G是积分常数;
则在设定的第N个周期时,轨道倾角变化极大值iN为:
3-3-2)对于升交点赤经:
使得升交点赤经一个周期变化量取得极大值的最优摄动加速度在密切轨道平面内垂直速度方向分量为:
将式(21)的结果代入式(10)里积分,得到升交点赤经在一个周期的变化量极大值为:
则在设定第N个周期时,升交点赤经变化极大值ΩN表示为:
将式(20)代入式(23)中积分,得到在第N个周期时,升交点赤经变化的极大值为:
本发明的特点及有益效果在于:
1.本发明基于微分代数原理,利用轨道动力学特性,分析推导出小推力航天器轨道根数关于轨道周期的长期演化变化量极大值估计方法,能够快速计算出航天器轨道根数变化量极大值。相对于数值递推计算方式和单小推力方向的轨道根数变化估计方法,本发明方法属于一种半解析方式,具有计算效率高,适用范围广等优点。
2.本发明中,对于部分轨道根数的长期演化,轨道根数变化量估计采用迭代算法,能够快速计算出相应的轨道根数变化量极大值,相比于数值积分算法能节省大量计算量。另外本方法可以针对航天器轨道周期量级上的轨道根数变化量估计,能够针对短期或者长期的轨道根数变化估计预报。
3.在航天器动力学与控制领域,本发明可以应用于小推力轨迹优化、航天器轨道根数变化预报、航天任务可行性分析以及轨道长期演化分析。能够适用于不同类型的轨道情况,对不同轨道根数情况均能较好地进行估计,与数值准确方法相比,具有估计精度高,计算效率高等优点。
附图说明
图1为本发明方法的整体流程图;
图2为本发明实施例中半长轴变化量极大值曲线示意图;
图3为本发明实施例中偏心率变化量极大值曲线示意图;
图4为本发明实施例中近拱点幅角变化量极大值曲线示意图;
图5为本发明实施例中轨道倾角变化量极大值曲线示意图;
图6为本发明实施例中升交点赤经化量极大值曲线示意图。
具体实施方式
本发明提出一种小推力航天器轨道根数长期演化的极大值估计方法,下面结合附图和具体实施例进一步详细说明如下。
本发明提出一种小推力航天器轨道根数长期演化的极大值估计方法,整体流程如图1所示,包括以下步骤:
1)建立航天器摄动方程:
对于高斯形式的拉格朗日行星运动方程,航天器轨道根数a,e,i,Ω,ω,f分别代表半长轴、偏心率、轨道倾角、升交点赤经、近拱点幅角、真近点角。令T=[av,an,ah]为小推力航天器的摄动加速度,其中三个分量分别代表摄动加速度在密切轨道平面内沿速度方向、垂直速度方向以及垂直于轨道平面的分量。
建立航天器摄动方程表达式如下:
其中,θ=ω+f表示真近点角与近地点幅角之和。r,v,h,μ分别为卫星地心距、速率、角动量标量和中心天体引力常数。t表示时间变量。
2)摄动方程转换;
对于步骤1)得到的高斯型摄动方程式(1)-(6),摄动加速度对真近点角的改变量为相对小量,仅考虑二体情况下真近点角的变化,用前五个轨道根数的微分方程式(1)—(5)分别除以第六个方程式(6),可以得到前五个轨道根数a,e,i,Ω,ω关于真近点角f的微分方程,从而得到以真近点角为自变量的高斯摄动方程式(7)-(11),表达式如下:
3)计算航天器轨道根数长期变化量的极大值;具体步骤如下:
3-1)设定任意满足开普勒椭圆轨道的轨道根数初值并分别记作a0,e0,i0,Ω0,ω0,f0。同时设定用带上标*的加速度表示使得轨道根数在一个周期变化量取极大值的最优加速度,其中三个分量分别代表最优摄动加速度在密切轨道平面内沿速度方向、垂直速度方向以及垂直于轨道平面的分量。
3-2)计算半长轴、偏心率、近拱点幅角长期变化量的极大值;具体步骤如下:
3-2-1)设定初始周期数N=1,则对于半长轴、偏心率、近拱点幅角存在:aN-1=a0,eN-1=e0,ωN-1=ω0;
3-2-2)从半长轴、偏心率、近拱点幅角中任意选取一个轨道根数;
3-2-3)计算步骤3-2-2)选取的轨道根数在第N-1个周期的变化量极大值;方法如下:
3-2-3-1)对于半长轴;
3-2-3-2)对于偏心率;
其中,||T||表示加速度矢量的模。AN-1,BN-1为用于简化加速度表示形式所引入的中间量,eN-1为第N-1个周期的偏心率。
3-2-3-3)对于近拱点幅角;
利用最优性原理可以得到使得近拱点幅角在第N-1个周期变化率最大的最优摄动加速度的三个分量分别为:
其中为了简化上述加速度分量计算中所引入的中间量Aω,N-1,Bω,N-1,Cω,N-1可以表示为:
3-2-4)按照下式更新轨道根数值,得到步骤3-2-2)选取的轨道根数在第N个周期的值为:
3-2-5)令N=N+1,然后重新返回步骤3-2-3),直至N的取值到达设定的数值,即可得到该轨道根数在第N-1个周期变化的极大值;
3-3)计算轨道倾角、升交点赤经长期变化量的极大值;
从轨道倾角和升交点赤经中任意选取一个轨道根数,计算该选取的轨道根数在第N个周期的变化量极大值,方法如下:
3-3-1)对于轨道倾角:
利用最优性原理得到使得轨道倾角一个周期变化量取得极大值的最优摄动加速度在密切轨道平面内垂直速度方向分量为:
将式(18)的结果代入轨道倾角的微分方程式(9)里积分可以得到轨道倾角在一个周期的变化量极大值(为常数)为:
其中H,G是积分常数;
则在设定的第N个周期时,轨道倾角变化极大值iN为:
3-3-2)对于升交点赤经:
利用最优性原理得到使得升交点赤经一个周期变化量取得极大值的最优摄动加速度在密切轨道平面内垂直速度方向分量为:
将式(21)的结果代入升交点赤经微分方程式(10)里积分可以得到升交点赤经在一个周期的变化量极大值为:
则在设定第N个周期时,升交点赤经变化极大值ΩN可以表示为:
将轨道倾角变化的极大值方程式(20)代入上式(23)中积分,可以得到在第N个周期时,升交点赤经变化的极大值为:
至此,即可以得到上述五个轨道根数的长期演化变化量极大值估计值。
实例验证
本实施例中,航天器初始轨道根数及长期演化目标轨道圈数等仿真参数如下表所示:
表1实例验证轨道根数初值设置表
设定小推力加速度大小为||T||=10-4m/s2。
本实施例中步骤3)的详细步骤包括:
3)计算航天器轨道根数长期变化量的极大值;具体步骤如下:
3-2)计算半长轴、偏心率、近拱点幅角长期变化量的极大值;具体步骤如下:
3-2-1)设定初始周期数N=1,对于给定的初始轨道的轨道根数有aN-1=7500km,eN-1=0.001,ωN-1=10;
3-2-2)从半长轴、偏心率、近拱点幅角中任意选取一个轨道根数;
3-2-3)计算步骤3-2-2)选取的轨道根数在第N-1个周期的变化量极大值;方法如下:
3-2-3-1)对于半长轴;
3-2-3-2)对于偏心率;
3-2-3-3)对于近拱点幅角;
其中:
3-2-4)按照下式更新轨道根数值,得到步骤3-2-2)选取的轨道根数在第N个周期的值为:
3-2-5)令N=N+1,然后重新返回步骤3-2-3),直至取值N=[2500,3000,300],即可到该轨道根数在设定周期变化的极大值;
图2、图3和图4分别为半长轴、偏心率、近拱点幅角变化极大值与数值解对比图,本发明所计算的对于不同的周期数N下的轨道根数变化曲线与精确的数值计算的轨道根数变化曲线相比两者基本保持一致,说明本实例采用的计算半长轴、偏心率、近拱点幅角变化极大值的方法能够很好的近似精确的轨道根数变化极大值,验证了本方法的正确性和有效性。
3-3)计算轨道倾角、升交点赤经长期变化的极大值;具体方法如下:
对于轨道倾角和升交点赤经,在进行这两个轨道根数的长期变化量近似时,忽略中其他轨道根数的变化以及中除轨道倾角i之外其他轨道根数的变化,即计算轨道倾角和升交点赤经变化的极大值时a0=7500km,e0=0.001,ω0=10。
从轨道倾角和升交点赤经中任意选取一个轨道根数,计算选取的轨道根数在第N个周期的变化量极大值;方法如下
3-3-1)对于轨道倾角:
利用最优性原理得到使得轨道倾角一个周期变化量取得极大值的最优摄动加速度在密切轨道平面内垂直速度方向分量为:
将式(18)的结果代入轨道倾角的微分方程式(9)里积分可以得到轨道在一个周期的变化量极大值(为常数):
其中H,G是积分常数;
则在设定的周期数N=3000时,轨道倾角变化极大值iN为:
其中i0=10。
3-3-2)对于升交点赤经:
利用最优性原理得到使得升交点赤经一个周期变化量取得极大值的最优摄动加速度在密切轨道平面内垂直速度方向分量为:
将式(21)的结果代入升交点赤经微分方程式(10)里积分可以得到轨道在一个周期的变化量极大值:
则在设定的周期数N时,升交点赤经变化极大值ΩN可以表示为:
将轨道倾角变化的极大值方程式(20)代入上式(23)中积分,可以得到在周期数为N=3000时,升交点赤经变化的极大值为:
其中,Ω0=30,i0=10,代入上式即可得升交点赤经在周期数为3000时变化的极大值。
图5、图6分别为轨道倾角、升交点赤经变化的极大值对比所示,实线和虚线分别表示本发明所示方法计算的轨道根数变化极大值和精确的数值极大值解。通过图5、图6可以发现本发明方法计算的轨道倾角与升交点赤经变化的极大值与精确数值解计算结果基本一致,说明本方法在可以精确估计轨道倾角和升交点赤经变化的极大值,表明本方法的计算的正确性和有效性。
综上所述,本实施例采用一种小推力轨道根数长期演化变化的极大值估计方法,应用于估计航天器轨道根数变化的极大值范围估计,与精确的数值仿真计算结果对比,能够精确有效地估计航天器在不同周期下的轨道根数变化的极大值。本实施例针对任意形式的航天器开普勒椭圆轨道,基于微分代数和轨道动力学特性,能够方便地应用于各种形式的轨道以及不同轨道周期下的轨道根数变化极大值估计。
Claims (1)
1.一种小推力航天器轨道根数长期演化的极大值估计方法,其特征在于,包括以下步骤:
1)建立航天器摄动方程:
对于高斯形式的拉格朗日行星运动方程,航天器轨道根数a,e,i,Ω,ω,f分别代表半长轴、偏心率、轨道倾角、升交点赤经、近拱点幅角、真近点角;
令T=[av,an,ah]为小推力航天器的摄动加速度,其中av代表摄动加速度在密切轨道平面内沿速度方向的分量,an代表摄动加速度在垂直速度方向的分量,ah代表摄动加速度在垂直于轨道平面的分量;
建立航天器摄动方程表达式如下:
其中,θ=ω+f表示真近点角与近地点幅角之和;r,v,h,μ分别为卫星地心距、速率、角动量标量和中心天体引力常数;t表示时间变量;
2)摄动方程转换;
对于步骤1)得到的摄动方程式(1)-(6),将式(1)-(5)分别除以式(6),得到轨道根数a,e,i,Ω,ω关于真近点角f的微分方程,从而得到以真近点角为自变量的高斯摄动方程式(7)-(11),表达式如下:
3)计算航天器轨道根数长期变化量的极大值;具体步骤如下:
3-1)设定半长轴、偏心率、轨道倾角、升交点赤经、近拱点幅角和真近点角的初值分别为a0,e0,i0,Ω0,ω0,f0,设定表示使得轨道根数在一个周期变化量取极大值的最优摄动加速度,其中三个分量分别代表最优摄动加速度在密切轨道平面内沿速度方向、垂直速度方向以及垂直于轨道平面的分量;
3-2)计算半长轴、偏心率、近拱点幅角长期变化量的极大值;具体步骤如下:
3-2-1)设定初始周期数N=1,则对于半长轴、偏心率、近拱点幅角存在:aN-1=a0,eN-1=e0,ωN-1=ω0;其中,aN-1为第N-1个周期的半长轴,eN-1为第N-1个周期的偏心率,ωN-1为第N-1个周期的近拱点幅角;
3-2-2)从半长轴、偏心率、近拱点幅角中任意选取一个轨道根数;
3-2-3)计算步骤3-2-2)选取的轨道根数在第N-1个周期的变化量极大值;方法如下:
3-2-3-1)对于半长轴;
3-2-3-2)对于偏心率;
AN-1=2(eN-1+cosf)
其中,||T||表示加速度矢量的模;AN-1,BN-1为用于简化加速度表示形式所引入的中间量;
3-2-3-3)对于近拱点幅角;
使得近拱点幅角在第N-1个周期变化率最大的最优摄动加速度的三个分量分别为:
其中,Aω,N-1,Bω,N-1,Cω,N-1为中间量,计算表达式分别如下:
3-2-4)计算步骤3-2-2)选取的轨道根数在第N个周期的值为:
3-2-5)令N=N+1,然后重新返回步骤3-2-3),直至N的取值到达设定的数值时,得到该轨道根数在第N-1个周期变化的极大值;
3-3)计算轨道倾角、升交点赤经长期变化量的极大值;
从轨道倾角和升交点赤经中任意选取一个轨道根数,计算该轨道根数在第N个周期的变化量极大值,方法如下:
3-3-1)对于轨道倾角:
使得轨道倾角一个周期变化量取得极大值的最优摄动加速度在密切轨道平面内垂直速度方向分量为:
将式(18)的结果代入式(9)里积分,得到轨道倾角在一个周期的变化量极大值为:
其中H,G是积分常数;
则在设定的第N个周期时,轨道倾角变化极大值iN为:
3-3-2)对于升交点赤经:
使得升交点赤经一个周期变化量取得极大值的最优摄动加速度在密切轨道平面内垂直速度方向分量为:
将式(21)的结果代入式(10)里积分,得到升交点赤经在一个周期的变化量极大值为:
则在设定第N个周期时,升交点赤经变化极大值ΩN表示为:
将式(20)代入式(23)中积分,得到在第N个周期时,升交点赤经变化的极大值为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011072998.8A CN112464429B (zh) | 2020-10-09 | 2020-10-09 | 一种小推力航天器轨道根数长期演化的极大值估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011072998.8A CN112464429B (zh) | 2020-10-09 | 2020-10-09 | 一种小推力航天器轨道根数长期演化的极大值估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464429A CN112464429A (zh) | 2021-03-09 |
CN112464429B true CN112464429B (zh) | 2022-09-23 |
Family
ID=74833165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011072998.8A Active CN112464429B (zh) | 2020-10-09 | 2020-10-09 | 一种小推力航天器轨道根数长期演化的极大值估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464429B (zh) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110789739B (zh) * | 2019-11-08 | 2020-12-11 | 中国人民解放军国防科技大学 | 一种j2摄动下长时间轨道交会最优速度增量快速估计方法 |
CN111268176B (zh) * | 2020-01-17 | 2021-06-11 | 中国人民解放军国防科技大学 | 一种摄动轨道四脉冲交会快速优化方法 |
CN111547274B (zh) * | 2020-03-19 | 2023-08-29 | 上海航天控制技术研究所 | 一种航天器高精度自主目标预报方法 |
-
2020
- 2020-10-09 CN CN202011072998.8A patent/CN112464429B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112464429A (zh) | 2021-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xu et al. | Quasi-equilibrium glide adaptive guidance for hypersonic vehicles | |
CN104309822B (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
CN105865459B (zh) | 一种考虑视线角约束的小天体接近段制导方法 | |
CN109613822B (zh) | 基于非线性模型预测控制的空间系绳系统展开控制方法 | |
CN110007681B (zh) | 一种利用连续推进器实现绳系编队自旋稳定展开优化方法 | |
CN109269504B (zh) | 一种具有末端约束的姿态机动路径规划方法 | |
CN110466805B (zh) | 基于优化制导参数的小行星着陆制导方法 | |
CN112572835B (zh) | 一种具有姿态切换的卫星在轨角动量管理及控制方法 | |
CN111731513B (zh) | 一种基于单脉冲轨控的高精度引力场中回归轨道维持方法 | |
CN109657256A (zh) | 一种高精度弹道式再入标称返回轨道仿真方法 | |
CN115202380B (zh) | 基于分段多项式的地外天体着陆轨迹规划方法 | |
CN115373423A (zh) | 一种用于商业卫星的编队捕获方法 | |
CN113343442B (zh) | 一种求解固定时间有限燃料多脉冲转移轨道的方法及系统 | |
CN112464429B (zh) | 一种小推力航天器轨道根数长期演化的极大值估计方法 | |
CN108082538B (zh) | 一种考虑始末约束的多体系统低能量捕获轨道方法 | |
Doupe et al. | Optimal attitude control of agile spacecraft using combined reaction wheel and control moment gyroscope arrays | |
CN111814313A (zh) | 一种高精度引力场中回归轨道设计方法 | |
CN115892519A (zh) | 一种用于近距离航天器轨道脉冲博弈的航天器控制方法 | |
CN113128032B (zh) | 一种基于轨道解析摄动解的交点时刻及位置求解算法 | |
CN115848648A (zh) | 一种动能撞击小行星的伪谱轨道优化方法和装置 | |
Song et al. | A lunar cargo mission design strategy using variable low thrust | |
CN116520719A (zh) | 适用于挠性航天器姿态动力学全物理仿真的方法及系统 | |
Cianciolo et al. | Autonomous aerobraking development software: Phase 2 summary | |
CN112455725A (zh) | 一种脉冲轨道转移向有限推力轨道转移转换的方法 | |
Kluever et al. | Trajectory-tracking guidance law for low-thrust earth-orbit transfers |
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 |