CN105574261A - 一种月球借力约束的地月平动点转移轨道设计方法 - Google Patents
一种月球借力约束的地月平动点转移轨道设计方法 Download PDFInfo
- Publication number
- CN105574261A CN105574261A CN201510937574.6A CN201510937574A CN105574261A CN 105574261 A CN105574261 A CN 105574261A CN 201510937574 A CN201510937574 A CN 201510937574A CN 105574261 A CN105574261 A CN 105574261A
- Authority
- CN
- China
- Prior art keywords
- moon
- constraint
- spacecraft
- centerdot
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种月球借力约束的地月平动点转移轨道设计方法,属于航天器轨道设计与优化领域。选择目标Halo轨道入轨点位置与期望的月球借力约束大小,利用方程(1)对入轨点的状态量进行积分,至航天器到达预期的月球借力位置。通过相关的优化算法对入轨点状态量进行修正调整,使得借力位置的状态满足方程(3)的约束条件,进而确定航天器进入Halo轨道所需要的机动速度增量大小。通过合理选择月球借力约束条件及Halo轨道入轨点等参数,能够有效地设计出满足任务要求的低耗能转移轨道。
Description
技术领域
本发明涉及一种月球借力约束的地月平动点转移轨道设计方法,属于航天器轨道设计与优化领域。
背景技术
近年来,地月***L2平动点(EML2)附近Halo轨道成为航天工程及应用领域的研究重点。通过EML2点Halo轨道,探测器能够观测月球背面并为月球着陆任务提供导航策略。同时,围绕EML2点能够设计出到达月球、火星与其它行星低耗能的转移轨道,可以为探测深空环境提供最佳的落脚点(徐明,徐世杰.地-月系平动点及Halo轨道的应用研究[J].宇航学报.2006,27(4):695-699)。因此,为了降低平动点转移轨道设计过程中的燃料消耗,有必要对探测器采用的飞行方式进行研究。
对于EML2点转移轨道设计问题,轨道转移方式通常采用霍曼直接转移,间接转移,弱稳定转移与月球借力转移。例如Zazzera采用遗传算法与序列二次规划的混合优化方法,对不考虑月球借力的多条间接转移轨道进行分析(ZazzeraFB,TopputoF,MassariM.Assessmentofmissiondesignincludingutilizationoflibrationpointsandweakstabilityboundaries[R].ESTECContractNo.18147/01/NL/MV,2004);Parker基于弱稳定边界与不变流形转移策略,通过连接日地***与地月***的不变流形,设计了从低地球轨道到达EML2点Halo轨道的长时间飞行转移轨道(ParkerJS.Low-EnergyBallisticLunarTransfers[D].Colorado:UniversityofColorado,2007);Gordon利用稳定流形理论与月球借力技术,构造了两脉冲转移轨道,但降低燃耗的能力有限;Li则在Gordon的基础上,设计了低燃耗的三脉冲转移轨道方案,但并未给出Halo轨道入轨点(HOI)的选择策略与分析所考虑的月球借力约束及其影响,在工程应用中存在着一定的局限性(GordonDP.Transferstoearth-moonL2haloorbitsusinglunarproximityandinvariantmanifolds[D].Indiana:PurdueUniversity,2008;LiMT,ZhengJH.Impulsivelunarhalotransfersusingthestablemanifoldsandlunarflybys[J].ActaAstronautica,2010,66:1481-1492)。
因此,对以EML2点Halo轨道为目标轨道的深空探测器而言,在低燃耗条件下设计转移轨道的技术成为一个关键性问题,如何选择月球借力约束条件,并综合考虑Halo轨道入轨点的影响成为突出问题。
发明内容
本发明的目的是为了降低平动点转移轨道设计过程中的燃耗问题。本发明公开的一种月球借力约束的地月平动点转移轨道设计方法,通过合理选择月球借力约束条件及Halo轨道入轨点等参数,能够有效地设计出满足任务要求的低耗能转移轨道。
本发明的目的是通过下述技术方案实现的:
一种月球借力约束的地月平动点转移轨道设计方法,具体步骤如下:
步骤一、航天器施加第一次脉冲,由目标Halo轨道反推至月球借力位置;
选择目标Halo轨道入轨点位置与期望的月球借力约束大小,利用方程(1)对入轨点的状态量进行积分,至航天器到达预期的月球借力位置。通过相关的优化算法对入轨点状态量进行修正调整,使得借力位置的状态满足方程(3)的约束条件,进而确定航天器进入Halo轨道所需要的机动速度增量大小。
步骤二、针对不同的任务轨道,航天器在借力点施加不同的脉冲次数,使得航天器由月球借力位置进入地-月转移轨道段,最终反推至地球停泊轨道。
当所述任务轨道包含绕月飞行的地-月转移轨道时,(a)选择绕月飞行轨迹末端满足的约束大小,利用方程(1)对步骤一中确定的月球借力点状态量进行积分,通过相关的优化算法对状态量进行修正调整,使得借力位置的状态满足方程(3)的约束条件,进而确定航天器绕月飞行所需的速度增量大小。(b)选择期望的地球停泊轨道约束大小,利用方程(1)对(a)中确定的绕月飞行轨迹末端的状态进行积分,至航天器到达预期的地球停泊轨道位置,通过相关的优化算法对绕月飞行轨迹末端的状态量进行修正调整,进而确定航天器进入地-月转移轨道段所需的速度增量大小、地球逃逸点速度增量大小与总飞行时间。
当所述任务轨道为地-月转移轨道时,选择期望的地球停泊轨道约束大小,利用方程(1)对步骤一中确定的月球借力点状态量进行积分,至航天器到达预期的地球停泊轨道位置。通过相关的优化算法对月球借力点状态量进行修正调整,进而确定航天器在月球借力点与地球逃逸点的速度增量大小与总飞行时间。
以上两种任务轨道,在地球停泊轨道上都满足高度约束与航迹角约束,并且逃逸点机动速度增量方向与航天器速度方向共线。
步骤三、整个平动点转移轨道在设计过程中采用的是逆向积分的策略,因此,步骤一和步骤二中施加的脉冲与实际工程中航天器转移轨道中机动顺序相反。
步骤一所述方程(1)为在地球-月球质心旋转系下建立的航天器动力学方程。
其中坐标系的原点为***的质心,x轴由地球指向月球,z轴与地-月***的角速度方向一致,y轴与x轴,z轴垂直,构成右手坐标系。
在质心旋转系下,航天器的动力学模型表示为
其中表征航天器在质心旋转系中的位置、速度与加速度状态量,μ表示***的质量系数,为航天器与地球之间的距离,为航天器与月球之间的距离。
所述月球借力位置的约束条件为:在转移轨道设计过程中,当月球借力位置考虑不同的约束条件时,将对整个轨道转移过程的速度增量与飞行时间等参数产生较大的影响。为了达到降低燃耗的目的,考虑月球借力位置的约束条件分别为:近月高度hm,相对月球飞行航迹角γm与借力方位角δm,可表示为:
其中,rfm=[xfm,yfm,zfm]T与为旋转系中借力点位置与速度矢量,Pfm=[xfm-1+μ,yfm,zfm]T为航天器相对于月球质心的位置矢量,Rm为月球半径,矢量
同时,为了充分利用目标Halo轨道的特性,采用逆向积分的策略设计地月平动点转移轨道。在保证高度hm与航迹角γm约束的条件下,改变借力方位角δm的取值大小,能够得到不同约束条件时燃料消耗变化情况,进而确定设计方案的各项参数。
所述采用逆向积分的策略,即为将步骤一中的公式(1)逆向积分;
步骤一、步骤二中的优化算法为微分修正算法,其作为一种局部优化算法,具有快速收敛特性,基于此算法对由公式(1)积分确定的轨迹进行调整,可以保证轨迹末端的状态量满足给定的约束条件。首先,选取自由变量C与约束矢量F(C),满足
其中,[Vx,Vy,Vz]T为Halo轨道入轨点三轴方向速度矢量,T为Halo轨道至月球借力位置的飞行时间;(hm,γm,δm)由公式(2)确定,为给定的期望约束值。
然后,推导约束矢量关于控制变量的偏导数关系
其中,符号δ表示变分关系,为约束矢量关于速度变量Vx的偏导数。最后利用最小二乘法对自由变量C进行微分修正,迭代求解能够保证最终的转移轨道满足给定的约束条件。
有益效果
1、本发明公开的一种月球借力约束的地月平动点转移轨道设计方法,由于同时考虑了多种月球借力约束条件,相比于不考虑月球借力或少约束情况的速度增量小,进而降低了任务燃料消耗。
2、本发明公开的一种月球借力约束的地月平动点转移轨道设计方法,根据不同的飞行时间与燃料消耗等任务要求,选择不同的目标Halo轨道、不同的Halo轨道入轨位置、航天器在月球借力点满足的约束条件,可以完成相应的转移轨道设计任务。
3、本发明公开的一种月球借力约束的地月平动点转移轨道设计方法,在月球借力点多施加一次脉冲机动,能够形成绕月飞行轨道,从而近距离的对月球背面进行观测,为未来月球背面搭建无线电望远镜时进行物质运输投放任务提供参考。
附图说明
图1为月球借力约束下的三脉冲与绕月飞行四脉冲轨道设计示意图;
图2为一种月球借力约束的地月平动点转移轨道方法流程图;
图3为公式(1)中地球-月球质心旋转坐标系的示意图;
图4为步骤一中月球借力点约束条件示意图;
图5为本发明实施例1不同借力约束条件的三脉冲转移轨道;
图6为本发明实施例1不同参数下转移轨道所需速度增量对比分析图;
图7为本发明实施例2不同借力约束条件的绕月飞行四脉冲转移轨道。
具体实施方式
为了更好地说明本发明的目的和优点,下面将结合附图对本发明作进一步的详细说明。
实施例1
如图1、图2所示,针对三脉冲设计方案,一种月球借力约束的地月平动点转移轨道设计方法,包括如下步骤:
步骤一、航天器施加第一次脉冲,由目标Halo轨道反推至月球借力位置,月球借力点约束如图4所示;
选择法向幅值Az=10000km的Halo轨道为目标轨道,期望的月球借力约束条件满足:近月高度飞行航迹角借力方位角选取Halo轨道入轨点状态为[1.12507,-0.08292,-0.00461,-0.03415,0.06426,-0.04819]T,利用公式(1)对Halo轨道入轨点状态量进行积分,在满足航迹角约束的位置中,距离月球最近点的高度为hm=1.3978×104km,借力方位角δm=16.434°,不满足期望的借力高度与方位角约束,采用相关的优化算法对入轨点状态量进行修正调整,使得借力位置的状态量满足公式(3)的约束矢量,进而确定航天器进入Halo轨道所需的机动速度增量大小为39.62m/s。
步骤二、航天器在月球借力位置施加第二次脉冲,使航天器由月球借力位置进入地-月转移轨道段,最终反推至地球停泊轨道。
选取期望的地球停泊轨道约束条件满足:轨道高度飞行航迹角利用公式(1)对月球借力位置状态进行积分,通过相关的优化算法进一步对月球借力位置状态量修正调整,直至满足给定的停泊轨道约束条件,则计算出航天器在月球借力点的速度增量大小为190.24m/s,总飞行时间为13.63天,从地球停泊轨道逃逸的速度增量为3120.33m/s。
当月球借力位置考虑少约束条件,即仅包含高度约束与航迹角约束时,设计整个平动点转移轨道所需的燃耗为:地球逃逸点为3133.41m/s,月球借力点为685.01m/s,Halo轨道入轨点为9.67m/s,而总飞行时间为13.25天。比较以上两种情况可知,在飞行时间相近条件下,考虑方位角约束后,总速度增量显著降低,两者相差477.9m/s。
针对以上两种情况设计了三脉冲转移轨道,如图5所示,其中虚线为仅考虑借力高度与航迹角约束的转移轨道,实线为添加借力方位角后的转移轨道。
为了确定Halo轨道入轨点的选择依据,选取不同目标Halo轨道、不同Halo轨道入轨位置、航天器在月球借力点满足的约束条件进行设计与分析,不同参数下的轨道所需速度增量的对比分析图如图6所示。由图6可知:对于相同编号入轨点,Halo轨道入轨点速度增量随着法向幅值的增大而增加,借力点机动速度增量不断下降,总速度增量逐步增大,范围满足ΔVtotal=(3320,3340)m/s。同时,能够满足约束的入轨点数目随着借力方位角的不断减小而减少,综合分析数据,能够较好地满足月球借力高度、航迹角与借力方位角约束的Halo轨道入轨点选取区间为(90,160)与(240,315)。以上参数能够为平动点轨道的设计与应用、参数的选择等方面提供依据。
步骤三、整个平动点转移轨道在设计过程中采用的是逆向积分策略,因此,步骤一和步骤二中施加的脉冲顺序与实际工程中航天器转移轨道中机动顺序相反。即航天器首先施加一次脉冲从地球停泊轨道逃逸,进入地-月转移轨道段,到达满足任务要求的月球借力位置,施加第二次脉冲,进入月球-Halo转移轨道,最终在Halo轨道入轨点施加第三次脉冲,实现目标轨道捕获。
实施例2
如图1、图2所示,针对绕月四脉冲设计方案,一种月球借力约束的地月平动点转移轨道设计方法,包括如下步骤:
步骤一、航天器施加第一次脉冲,由目标Halo轨道反推至月球借力位置,月球借力点约束如图4所示;
选择法向幅值Az=5000km的Halo轨道为目标轨道,期望的月球借力约束条件满足:近月高度飞行航迹角借力方位角选取Halo轨道入轨点状态为[1.12135,-0.05816,-0.00740,-0.01368,0.12341,-0.00535]T,利用公式(1)对Halo轨道入轨点状态量进行积分,在满足航迹角约束的位置中,距离月球最近点的高度为hm=2.18664×104km,借力方位角δm=7.54890°,不满足期望的借力不满足期望的借力高度与方位角约束,采用相关的优化算法对入轨点状态量进行修正调整,使得借力位置的状态量满足公式(3)的约束矢量,进而确定航天器进入Halo轨道所需的机动速度增量大小为15.61m/s。
步骤二、航天器在月球借力位置施加第二次脉冲,形成绕月飞行轨迹,在施加第三次脉冲,使航天器由月球借力位置进入地-月转移轨道段,最终反推至地球停泊轨道。
(a)选择绕月飞行轨迹末端满足的约束条件为:近月高度飞行航迹角利用公式(1)对步骤一中确定的月球借力点状态量[0.98420,-0.00307,-0.00032,1.42238,-1.69059,-0.03850]T进行积分,并结合相关的优化算法对状态量进行修正调整,直至满足高度与航迹角约束条件,进而确定航天器绕月飞行所需的速度增量大小为51.38m/s,飞行时间1.33天。(b)选择期望的地球停泊轨道约束条件:停泊轨道高度航迹角利用公式(1)对(a)中确定的绕月飞行轨迹末端的状态量进行积分,状态量为[0.98485,-0.00467,-0.00038,1.67400,-1.07063,-0.04063]T。通过相关优化算法对状态量进行修正调整,直至满足给定的地球停泊轨道约束条件。最终确定航天器进入地-月转移轨道段所需的速度增量大小为250.92m/s,地球逃逸点为3118.35m/s与总飞行时间为14.72天。
当月球借力位置仅考虑高度约束与航迹角约束时,总速度增量大小为3490.896m/s,总飞行时间为14.52天,,比较以上两种情况,在飞行时间相近条件下,考虑方位角约束后,总速度增量降低,两者相差54.636m/s。
针对以上两种情况设计了绕月飞行四脉冲转移轨道,如图7所示,其中虚线为仅考虑借力高度与航迹角约束的转移轨道,实线为添加借力方位角后的转移轨道。
实施例1与实施例2所述公式(1)为在地球-月球质心旋转系下建立航天器的动力学方程。
如图3所示,其中坐标系的原点为***的质心,x轴由地球指向月球,z轴与地-月***的角速度方向一致,y轴与x轴,z轴垂直构成右手坐标系。
在质心旋转系下,航天器的动力学模型表示为
其中表征航天器在质心旋转系中的位置、速度与加速度状态量,μ表示***的质量系数,地-月***中,取值为航天器与地球之间的距离,为航天器与月球间的距离。
所述月球借力位置的约束条件为:在转移轨道设计过程中,当月球借力位置考虑不同的约束条件时,将对整个轨道转移过程的速度增量与飞行时间等参数产生较大的影响。为了达到降低燃耗的目的,考虑月球借力位置的约束条件分别为:近月高度hm,相对月球飞行航迹角γm与借力方位角δm,可表示为:
其中,rfm=[xfm,yfm,zfm]T与为旋转系中借力点位置与速度矢量,Pfm=[xfm-1+μ,yfm,zfm]T为航天器相对于月球质心的位置矢量,Rm为月球半径,矢量
同时,为了充分利用目标Halo轨道的特性,采用逆向积分的策略设计地月平动点转移轨道。在保证高度hm与航迹角γm约束的条件下,改变借力方位角δm的取值大小,能够得到不同约束条件时燃料消耗变化情况,进而确定设计方案的各项参数。
所述采用逆向积分的策略,即为将步骤一中的公式(1)逆向积分;
步骤一、步骤二中的优化算法为微分修正算法,其作为一种局部优化算法,具有快速收敛特性,基于此算法对由公式(1)积分确定的轨迹进行调整,可以保证轨迹末端的状态量满足给定的约束条件。首先,选取自由变量C与约束矢量F(C),满足
其中,[Vx,Vy,Vz]T为Halo轨道入轨点三轴方向速度矢量,T为Halo轨道至月球借力位置的飞行时间;(hm,γm,δm)由公式(2)确定,为给定的期望约束值。
然后,推导约束矢量关于控制变量的偏导数关系
其中,符号δ表示变分关系,为约束矢量关于速度变量Vx的偏导数。最后利用最小二乘法对自由变量C进行微分修正,迭代求解能够保证最终的转移轨道满足给定的约束条件。
尽管已经给出并描述了本发明的实施例,本领域的普通技术人员都能够理解,因此,在不脱离本发明设计原理与方法的情况下,可以对给定的实施例进行修改,参数替换。
Claims (8)
1.一种月球借力约束的地月平动点转移轨道设计方法,其特征在于:具体步骤如下:
步骤一、航天器施加第一次脉冲,由目标Halo轨道反推至月球借力位置;
选择目标Halo轨道入轨点位置与期望的月球借力约束大小,利用方程(1)对入轨点的状态量进行积分,至航天器到达预期的月球借力位置;通过相关的优化算法对入轨点状态量进行修正调整,使得借力位置的状态满足方程(3)的约束条件,进而确定航天器进入Halo轨道所需要的机动速度增量大小;
步骤二、针对不同的任务轨道,航天器在借力点施加不同的脉冲次数,使得航天器由月球借力位置进入地-月转移轨道段,最终反推至地球停泊轨道;
步骤三、整个平动点转移轨道在设计过程中采用的是逆向积分的策略,因此,步骤一和步骤二中施加的脉冲与实际工程中航天器转移轨道中机动顺序相反。
2.如权利要求1所述的一种月球借力约束的地月平动点转移轨道设计方法,其特征在于:当所述任务轨道包含绕月飞行的地-月转移轨道时,(a)选择绕月飞行轨迹末端满足的约束大小,利用方程(1)对步骤一中确定的月球借力点状态量进行积分,通过相关的优化算法对状态量进行修正调整,使得借力位置的状态满足方程(3)的约束条件,进而确定航天器绕月飞行所需的速度增量大小;(b)选择期望的地球停泊轨道约束大小,利用方程(1)对(a)中确定的绕月飞行轨迹末端的状态进行积分,至航天器到达预期的地球停泊轨道位置,通过相关的优化算法对绕月飞行轨迹末端的状态量进行修正调整,进而确定航天器进入地-月转移轨道段所需的速度增量大小、地球逃逸点速度增量大小与总飞行时间。
3.如权利要求1所述的一种月球借力约束的地月平动点转移轨道设计方法,其特征在于:当所述任务轨道为地-月转移轨道时,选择期望的地球停泊轨道约束大小,利用方程(1)对步骤一种确定的月球借力点状态量进行积分,至航天器到达预期的地球停泊轨道位置;通过相关的优化算法对月球借力点状态量进行修正调整,进而确定航天器在月球借力点与地球逃逸点的速度增量大小与总飞行时间。
4.如权利要求1或2或3所述的一种月球借力约束的地月平动点转移轨道设计方法,其特征在于:所述两种任务轨道,在地球停泊轨道上都满足高度约束与航迹角约束,并且逃逸点机动速度增量方向与航天器速度方向共线。
5.如权利要求1所述的一种月球借力约束的地月平动点转移轨道设计方法,其特征在于:步骤一所述公式(1)为在地球-月球质心旋转系下建立的航天器动力学方程;
其中坐标系的原点为***的质心,x轴由地球指向月球,z轴与地-月***的角速度方向一致,y轴与x轴,z轴垂直,构成右手坐标系;
在质心旋转系下,航天器的动力学模型表示为
其中表征航天器在质心旋转系中的位置、速度与加速度状态量,μ表示***的质量系数,为航天器与地球之间的距离,为航天器与月球之间的距离。
6.如权利要求1或2或3所述的一种月球借力约束的地月平动点转移轨道设计方法,其特征在于:所述月球借力位置的约束条件为:在转移轨道设计过程中,当月球借力位置考虑不同的约束条件时,将对整个轨道转移过程的速度增量与飞行时间等参数产生较大的影响;为了达到降低燃耗的目的,考虑月球借力位置的约束条件分别为:近月高度hm,相对月球飞行航迹角γm与借力方位角δm,可表示为:
其中,rfm=[xfm,yfm,zfm]T与为旋转系中借力点位置与速度矢量,Pfm=[xfm-1+μ,yfm,zfm]T为航天器相对于月球质心的位置矢量,Rm为月球半径,矢量
同时,为了充分利用目标Halo轨道的特性,采用逆向积分的策略设计地月平动点转移轨道;在保证高度hm与航迹角γm约束的条件下,改变借力方位角δm的取值大小,能够得到不同约束条件时燃料消耗变化情况,进而确定设计方案的各项参数。
7.如权利要求6所述的一种月球借力约束的地月平动点转移轨道设计方法,其特征在于:所述采用逆向积分的策略,即为将步骤一中的公式(1)逆向积分。
8.如权利要求1所述的一种月球借力约束的地月平动点转移轨道设计方法,其特征在于:步骤一、步骤二中的所述的优化算法为微分修正算法,基于此算法对由公式(1)积分确定的轨迹进行调整,可以保证轨迹末端的状态量满足给定的约束条件;首先,选取自由变量C与约束矢量F(C),满足
其中,[Vx,Vy,Vz]T为Halo轨道入轨点三轴方向速度矢量,T为Halo轨道至月球借力位置的飞行时间;(hm,γm,δm)由公式(2)确定,为给定的期望约束值;
然后,推导约束矢量关于控制变量的偏导数关系
其中,符号δ表示变分关系,为约束矢量关于速度变量Vx的偏导数;最后利用最小二乘法对自由变量C进行微分修正,迭代求解能够保证最终的转移轨道满足给定的约束条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510937574.6A CN105574261B (zh) | 2015-12-15 | 2015-12-15 | 一种月球借力约束的地月平动点转移轨道设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510937574.6A CN105574261B (zh) | 2015-12-15 | 2015-12-15 | 一种月球借力约束的地月平动点转移轨道设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105574261A true CN105574261A (zh) | 2016-05-11 |
CN105574261B CN105574261B (zh) | 2018-09-21 |
Family
ID=55884391
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510937574.6A Active CN105574261B (zh) | 2015-12-15 | 2015-12-15 | 一种月球借力约束的地月平动点转移轨道设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105574261B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105912020A (zh) * | 2016-05-26 | 2016-08-31 | 北京航空航天大学 | 一种弱稳定边界区域内航天器编队飞行方法 |
CN106682274A (zh) * | 2016-12-05 | 2017-05-17 | 北京理工大学 | 考虑振幅约束的一种Halo轨道在轨保持方法 |
CN107506505A (zh) * | 2016-06-14 | 2017-12-22 | 北京空间技术研制试验中心 | 高精度地月自由返回轨道设计方法 |
CN108519958A (zh) * | 2018-02-05 | 2018-09-11 | 中国人民解放军国防科技大学 | 一种解析构造航天器追逃界栅和判断捕获逃逸区域的方法 |
CN110096726A (zh) * | 2019-02-21 | 2019-08-06 | 上海卫星工程研究所 | 基于月球借力的geo卫星应急转移轨道快速优化设计方法 |
CN110733667A (zh) * | 2019-09-29 | 2020-01-31 | 北京空间技术研制试验中心 | 地月平动点轨道间转移设计方法 |
CN111737814A (zh) * | 2020-05-25 | 2020-10-02 | 北京理工大学 | 一种利用地月旁近的平衡点周期轨道逃逸方法 |
CN111753398A (zh) * | 2020-05-25 | 2020-10-09 | 北京理工大学 | 一种基于状态转移矩阵的平衡点转移轨道重构方法 |
CN112009727A (zh) * | 2020-08-21 | 2020-12-01 | 北京空间技术研制试验中心 | 平动点轨道最优小推力转移分段设计方法 |
CN115343960A (zh) * | 2022-10-19 | 2022-11-15 | 北京理工大学 | 一种地月***中航天器光照阴影规避控制方法 |
CN116541968A (zh) * | 2023-06-28 | 2023-08-04 | 南京航空航天大学 | 一种地月dro最优转移轨道的确定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1746025A1 (en) * | 2005-07-20 | 2007-01-24 | Japan Aerospace Exploration Agency | Method for designing an orbit of a spacecraft |
CN103488830A (zh) * | 2013-09-16 | 2014-01-01 | 北京航空航天大学 | 一种基于Cycler轨道的地月往返的任务仿真*** |
CN104657559A (zh) * | 2015-03-04 | 2015-05-27 | 中国人民解放军国防科学技术大学 | 基于圆柱型速度参数截面的地月自由返回轨道设计方法 |
-
2015
- 2015-12-15 CN CN201510937574.6A patent/CN105574261B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1746025A1 (en) * | 2005-07-20 | 2007-01-24 | Japan Aerospace Exploration Agency | Method for designing an orbit of a spacecraft |
CN103488830A (zh) * | 2013-09-16 | 2014-01-01 | 北京航空航天大学 | 一种基于Cycler轨道的地月往返的任务仿真*** |
CN104657559A (zh) * | 2015-03-04 | 2015-05-27 | 中国人民解放军国防科学技术大学 | 基于圆柱型速度参数截面的地月自由返回轨道设计方法 |
Non-Patent Citations (2)
Title |
---|
JEFFREY S PARKER: "TARGETING LOW-ENERGY BALLISTIC LUNAR TRANSFERS", 《JOURNAL OF THE ASTRONAUTICAL SCIENCES》 * |
李明涛 等: "基于月球借力的LLZ点Hlao轨道转移轨道研究", 《中国宇航学会深空探测技术专业委员会第五届学术年会论文集》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105912020A (zh) * | 2016-05-26 | 2016-08-31 | 北京航空航天大学 | 一种弱稳定边界区域内航天器编队飞行方法 |
CN105912020B (zh) * | 2016-05-26 | 2018-10-09 | 北京航空航天大学 | 一种弱稳定边界区域内航天器编队飞行方法 |
CN107506505A (zh) * | 2016-06-14 | 2017-12-22 | 北京空间技术研制试验中心 | 高精度地月自由返回轨道设计方法 |
CN107506505B (zh) * | 2016-06-14 | 2020-10-23 | 北京空间技术研制试验中心 | 高精度地月自由返回轨道设计方法 |
CN106682274B (zh) * | 2016-12-05 | 2019-12-13 | 北京理工大学 | 考虑振幅约束的一种Halo轨道在轨保持方法 |
CN106682274A (zh) * | 2016-12-05 | 2017-05-17 | 北京理工大学 | 考虑振幅约束的一种Halo轨道在轨保持方法 |
CN108519958B (zh) * | 2018-02-05 | 2022-02-08 | 中国人民解放军国防科技大学 | 一种解析构造航天器追逃界栅和判断捕获逃逸区域的方法 |
CN108519958A (zh) * | 2018-02-05 | 2018-09-11 | 中国人民解放军国防科技大学 | 一种解析构造航天器追逃界栅和判断捕获逃逸区域的方法 |
CN110096726A (zh) * | 2019-02-21 | 2019-08-06 | 上海卫星工程研究所 | 基于月球借力的geo卫星应急转移轨道快速优化设计方法 |
CN110096726B (zh) * | 2019-02-21 | 2023-08-11 | 上海卫星工程研究所 | 基于月球借力的geo卫星应急转移轨道快速优化设计方法 |
CN110733667A (zh) * | 2019-09-29 | 2020-01-31 | 北京空间技术研制试验中心 | 地月平动点轨道间转移设计方法 |
CN111737814A (zh) * | 2020-05-25 | 2020-10-02 | 北京理工大学 | 一种利用地月旁近的平衡点周期轨道逃逸方法 |
CN111753398A (zh) * | 2020-05-25 | 2020-10-09 | 北京理工大学 | 一种基于状态转移矩阵的平衡点转移轨道重构方法 |
CN111753398B (zh) * | 2020-05-25 | 2023-01-06 | 北京理工大学 | 一种基于状态转移矩阵的平衡点转移轨道重构方法 |
CN112009727A (zh) * | 2020-08-21 | 2020-12-01 | 北京空间技术研制试验中心 | 平动点轨道最优小推力转移分段设计方法 |
CN115343960A (zh) * | 2022-10-19 | 2022-11-15 | 北京理工大学 | 一种地月***中航天器光照阴影规避控制方法 |
CN116541968A (zh) * | 2023-06-28 | 2023-08-04 | 南京航空航天大学 | 一种地月dro最优转移轨道的确定方法 |
CN116541968B (zh) * | 2023-06-28 | 2023-09-22 | 南京航空航天大学 | 一种地月dro最优转移轨道的确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105574261B (zh) | 2018-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105574261A (zh) | 一种月球借力约束的地月平动点转移轨道设计方法 | |
CN106697333B (zh) | 一种航天器轨道控制策略的鲁棒性分析方法 | |
CN104309822B (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
CN106679674B (zh) | 基于星历模型的地月L2点Halo轨道阴影分析方法 | |
CN105912819B (zh) | 一种地月l1拉格朗日点转移轨道的快速设计方法 | |
CN109911249B (zh) | 低推重比飞行器的星际转移有限推力入轨迭代制导方法 | |
CN109470252A (zh) | 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 | |
CN101794154A (zh) | 一种编队卫星相对轨道与姿态解耦控制方法 | |
CN110053788B (zh) | 一种考虑复杂摄动的星座长期保持控制频次估计方法 | |
CN106021784B (zh) | 一种基于两层优化策略的全轨迹优化设计方法 | |
CN108100307A (zh) | 一种用于复杂约束下低能量小天体精确探测轨道转移方法 | |
CN112016187B (zh) | 一种基于混合动力的近地小行星交会任务轨道优化方法 | |
CN113343442B (zh) | 一种求解固定时间有限燃料多脉冲转移轨道的方法及*** | |
CN108128484A (zh) | 一种基于线性二次型调节器的双星***轨道保持方法 | |
CN109839116A (zh) | 一种速度增量最小的交会接近方法 | |
CN109977543B (zh) | 基于侧向优先的三维剖面覆盖区域计算方法、***及介质 | |
Kos et al. | Altair descent and ascent reference trajectory design and initial dispersion analyses | |
Mao et al. | REAL-TIME ONLINE TRAJECTORY PLANNING AND GUIDANCE FOR TERMINAL AREA ENERGY MANAGEMENT OF UNMANNED AERIAL VEHICLE. | |
CN108082538B (zh) | 一种考虑始末约束的多体***低能量捕获轨道方法 | |
CN107506505A (zh) | 高精度地月自由返回轨道设计方法 | |
Kang et al. | Nanosat formation flying design for SNIPE mission | |
CN115892519A (zh) | 一种用于近距离航天器轨道脉冲博弈的航天器控制方法 | |
CN114117631B (zh) | 一种带有最优终端时间估计的火箭回收轨迹优化方法 | |
Duan et al. | Trajectory tracking and online replanning for Mars entry | |
CN112009727B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |