不同月球借力约束下的地月Halo轨道转移轨道设计

针对地月系L2点不同任务需求下的低耗能转移轨道设计问题,基于不变流形理论与混合优化技术,深入研究了不同月球借力约束与不同幅值Halo轨道的入轨点（简称HOI点）对转移轨道飞行时间与燃料消耗的影响,给出了HOI点选择策略。首先结合任务要求并考虑月球引力影响,在月球借力点施加不同约束条件,通过微分修正算法调整Halo轨道的稳定流形,设计月球到Halo轨道的转移轨道。采用遗传算法与微分修正算法相结合的混合优化策略,在同时考虑地球停泊轨道高度、倾角、升交点赤经与航迹角等多约束条件下,对燃料最优的地月转移轨道进行研究。最后,分析月球借力高度、借力方位角和不同HOI点对平动点转移轨道飞行时间与燃耗变化量的影响,对于考虑月球借力的地月平动点转移轨道设计与应用具有重要的参考价值。     

“小飞”的杰作：地月合影

据国防科工局消息,11月28日,探月工程三期再入返回飞行试验器服务舱飞抵地月系统拉格朗日L2点,实施了地月L2点绕飞期间第一次轨道维持控制。在完成地月L2点探测活动后,服务舱将返回月球并进入环月轨道开展后续拓展试验项目。     

双三体系统不变流形拼接成的低成本探月轨道

传统的探月轨道设计原理为二体模型框架下的Hohmann变轨理论,但1991年日本的Hiten探月器利用太阳的摄动,用比传统的方法更少的燃料完成了探月任务。利用三体问题非线性系统的不变流形设计了节省燃料的探月轨道。沿用JPL研究组的思路,将太阳-地球-月亮-航天器四体问题分解成太阳-地球-航天器和地球-月亮-航天器两个共面的圆形限制性三体问题,对Hiten类的探月轨道给出了更深刻的数学、力学解释;给出了流形的结构以及更合理的拼接方式;找到了发射位置、发射速度和拼接点;设计出了类似Hiten探月器的探月轨道,可比传统方法节省速度增量12%左右。结果证明了三体系统不变流形在登月轨道设计研究中的可行性和优越性。     

约束条件下的Halo轨道转移轨道设计

平动点任务的转移轨道往往存在约束条件,以往的研究集中于无约束条件下的Halo轨道转移轨道设计,研究约束条件下的Halo轨道的转移轨道设计问题.首先分析了平动点任务转移轨道的约束条件,然后给出了一种约束条件下Halo轨道转移轨道设计的一般方法,重点推导了考虑轨道高度、航迹角、轨道倾角、升交点约束的微分修正公式.然后以日地L,点附近的Halo轨道为目标轨道,在约束条件下设计了其转移轨道,仿真结果验证了本文方法的有效性.     

日地Halo轨道的多约束转移轨道分层微分修正设计

针对日-地系统L1点（简称SEL1点）Halo轨道转移轨道设计中存在的多约束与初值敏感性问题,提出一种基于分层微分修正与初值多项式的设计方法。首先定义平动点转移轨道设计过程中存在的约束条件,然后根据不同的终端约束条件,重点给出了同时考虑轨道高度、轨道倾角、升交点赤经与航迹角等多约束条件下的分层微分修正方法。通过分析约束变量与控制变量之间的关系,得到能够解决微分修正初值问题的初值表达式。最后在多约束条件下设计了从轨道高度为200km的地球停泊轨道到SEL1点Halo轨道的转移轨道。仿真结果表明,分层微分修正方法能够处理多约束问题,且初值表达式可以为微分修正提供良好的初始条件,从而保证算法收敛,方法具有较好的实用性。     

绕月飞行的大幅值逆行轨道研究

以地-月系为背景,研究了绕月飞行的大幅值逆行轨道(DRO)的轨道稳定性,及轨道转移等问题。数值结果表明：太阳引力对DRO稳定性有破坏作用,但仍能保持较长时间的绕飞。随后,利用与DRO相切的Lyapunov轨道研究了DRO的低能轨道转移：利用地月系LL1点Lyapunov轨道的不变流形,实现DRO的快速转移;利用
LL2点Lyapunov轨道作为弱稳定边界（WSB）转移的入口,实现DRO的低能转移。显然,得到的两类转移方式完全不同于以往的研究,且数值仿真表明了设计方法的可行性。

     

地月低能转移轨道设计方法研究

研究地月低能转移轨道的设计方法。这种低能转移轨道利用了弱稳定边界理论，通过太阳的引力摄动，使得探测器能够不经过减速就被月球俘获。与经典的霍曼转移相比，低能转移轨道呵节省约140m／s的速度脉冲。由于设计是基于叫体问题模型进行的，问题具有很强的非线性特性，寻找满足约束条件的转移轨道变得非常困难。常用的两点边值问题的解法在这里都失效。本文在研究地月低能转移轨道特忡的基础上，对一般地月转移轨道搜索的变步长爬山法进行改进，用来设计地月低能转移轨道。仿真结果验证了该方法的有效性。     

发射短讯

嫦娥二号卫星进入环月长期运行轨道 2010年10月1日18点59分,长征-3C运载火箭从西昌卫星发射中心将中国探月工程二期的先导星——嫦娥二号卫星准确送入地月转移轨道。嫦娥二号卫星踏上约112h的奔月之旅。     

圆型限制性三体问题相对运动解析研究

针对圆型限制性三体问题共线平动点附近周期/拟周期轨道下的相对运动问题，提出一种新的、通用的解析研究方法。在周期/拟周期轨道近似解析解的基础上，结合微分修正方法，获得了精确的周期/拟周期轨道。对周期/拟周期轨道的单值矩阵进行分析，同时借鉴Floquet理论核心思想，建立了六个相对运动模态，并将相对运动表示为六个相对运动模态的线性组合，获得了相对运动的近似解析解。最后在地-月系统圆型限制性三体问题下，以L1点作为研究对象，分别以Halo轨道、Lissajous轨道和Lyapunov轨道为参考轨道，对相对运动模态和相对运动进行仿真分析，说明了相对运动模态的正确性以及相对运动近似解析解的有效性。     

Halo轨道维持的线性周期控制策略

共线型平动点附近的Halo轨道具有指数不稳定性,轨道维持是必不可少的。推导了基于Halo轨道的误差动力学方程,并证明其一阶近似即为线性周期系统。以一次维持的作用时间为离散步长,并通过定常变换,将所得误差动力学化为线性离散定常系统;则仅需通过极点配置,即可实现Halo轨道镇定。研究结果表明,利用Halo轨道周期性设计的线性周期控制策略,可以满足轨道维持任务的需要。     

The past and future of the earth-moon system

The lunar orbit is presently expanding due, we believe, to tidal friction, i.e. the attraction of the moon for the tides it raises on the rotating Earth. The Moon may therefore have been significantly closer to the Earth in the distant past, a point of great interest to those studying the lunar origin. This work presents the results of the integration of the equations which govern the rates of change of the lunar orbit elements and the angular momentum of the Earth. Results are presented for both the past and future of the Earth-Moon system.     

利用太阳帆定点探测器在地一月系共线平动点附近的探讨

共线平动点附近的运动仅仅是条件稳定的,探测器的轨道需要经过控制才能维持在其附近.以地-月系11点和12点附近大振幅晕轨道的控制为例,探讨了太阳帆在定点这类探测器中的应用.首先,考虑了月球轨道的偏心率和太阳辐射的影响,给出了太阳帆对日定向的探测器轨道的低阶分析解,并在此基础上构造了在太阳系真实引力模型下一段时间内维持在共线平动点附近的拟周期轨道.然后,给出了两种利用太阳帆的控制方案,一是固定面质比而改变太阳帆法线的方向,另一是固定太阳帆对日定向而改变面质比,并对两种方案分别作了数值模拟.最后,文章探讨了测控误差及地、月影对轨道控制的影响.     

Halo orbits in the vicinity of the L 2 libration point in the Sun-Earth system

Methods are proposed for constructing the orbits of spacecraft remaining for long periods of time in the vicinity of the L ₂ libration point in the Sun-Earth system (so-called halo orbits), and the trajectories of uncontrolled flights from low near-Earth orbits to halo orbits. Halo orbits and flight trajectories are constructed in two stages: A suitable solution to a circular restricted three-body problem is first constructed and then transformed into the solution for a restricted four-body problem in view of the real motions of the Sun, Earth, and Moon. For a halo orbit, its prototype in the first stage is a combination of a periodic Lyapunov solution in the vicinity of the L ₂ point and lying in the plane of large-body motion, with the solution for the linear second-order system describing small deviations of the spacecraft from this plane along the periodic solution. The desired orbit is found as the solution to the three-body problem best approximating the prototype in the mean square. The constructed orbit serves as a similar prototype in the second stage. In both stages, the approximating solution is constructed by continuation along a parameter that is the length of the approximation interval. Flight trajectories are constructed in a similar manner. The prototype orbit in the first stage is a combination of a solution lying in the plane of large-body motion and a solution for a linear second-order system describing small deviations of the spacecraft from this plane. The planar solution begins near the Earth and over time tends toward the Lyapunov solution existing in the vicinity of the L ₂ point. The initial conditions of both prototypes and the approximating solutions correspond to the spacecraft’s departure from a low near-Earth orbit at a given distance, perigee, and inclination.     

Autonomous optical navigation for orbits around Earth–Moon collinear libration points

The analysis of optical navigation in an Earth–Moon libration point orbit is examined. Missions to libration points have been winning momentum during the last decades. Its unique characteristics make it suitable for a number of operational and scientific goals. Literature aimed to study dynamics, guidance and control of unstable orbits around collinear libration points is vast. In particular, several papers deal with the optimisation of the Δv budget associated to the station-keeping of these orbits. One of the results obtained in literature establishes the critical character of the Moon–Earth system in this aspect. The reason for this behaviour is twofold: high Δv cost and short optimal manoeuvre spacing. Optical autonomous navigation can address the issue of allowing a more flexible manoeuvre design. This technology has been selected to overcome similar difficulties in other critical scenarios. This paper analyses in detail this solution. A whole GNC system is defined to meet the requirements imposed by the unstable dynamic environment. Finally, a real simulation of a spacecraft following a halo orbit of the L2 Moon–Earth system is carried out to assess the actual capabilities of the optical navigation in this scenario.     

地月空间NRHO与DRO在月球探测中的应用研究

针对地月系统三体问题低能往返轨道转移在月球探测中的应用研究，结合天体借力飞行技术与混合优化技术，系统分析了不同目标轨道与借力方位对任务飞行时间与燃料消耗等关键参数的影响，给出了往返轨道设计初值的选择策略。针对轨道设计初值猜想问题，首先采用遗传算法与二体Lambert转移快速确定轨迹拼接点初值。在同时考虑近月点与近地点多约束条件下，基于序列二次规划算法与多重打靶法进一步对燃料最优的地月往返轨道进行研究，并推导了约束方程解析梯度提高设计效率。最后分析近月点高度、不同目标轨道的转移时间与燃耗变化特性，对于考虑月球借力的地月空间轨道往返转移设计及参数选取具有重要的参考价值。     

地-月低能耗转移轨道中途修正问题研究

采用地-月低能耗转移轨道的探测器从地球停泊轨道转移到极月轨道一般需要3~4个月时间,这类转移轨道对入轨精度有较高的要求。本文对地月转移轨道中途修正问题进行了研究。文中结合地-月低能耗转移轨道的特点,给出一种分段式多目标多次中途修正方案。利用显式制导结合牛顿迭代,分别以地球和月球作为中心天体求解兰伯特问题,在假设探测器各种轨道误差的基础上进行了蒙特卡罗仿真。采用该方法一般需要3~5次中途修正能够满足月球探测器环月轨道入轨精度要求,整个转移过程燃料消耗小于传统地月转移轨道。文中给出的仿真结果验证了该方案的可行性。     

解析计算在月球中继卫星Halo轨道设计中的应用

面向月球中继卫星工程轨道设计需求,研究解析计算方法在地月系L2点halo轨道设计中的应用问题。在讨论圆型限制性三体问题三阶解析近似计算方法的基础上,分析了解析计算与数值计算的差异,给出了解析近似计算在工程约束下的适用范围,进而提出了基于解析计算的轨道设计和特征筛选方法。分别采用解析初值和数值初值进行halo轨道外推,比对验证采用解析计算设计轨道的可行性。研究结果表明,解析计算方法适用于月球中继卫星轨道的初步设计、特征分析和构型筛选。     

嫦娥五号平动点拓展任务轨道方案研究

中国嫦娥五号探测器成功实现月球采样返回任务,为最大限度利用任务资源,研究了利用嫦娥五号轨道器的平动点拓展任务轨道方案,设计了平动点轨道及其转移轨道.首先,给出了任务轨道设计的轨道动力学模型,包括圆型限制性三体问题模型和精确力模型.其次,基于嫦娥二号和嫦娥5T1平动点拓展任务设计经验,介绍了平动点轨道直接转移与入轨等轨道...     

定时定点月面着陆全程轨道控制设计

针对定时定点月面着陆的目标要求，提出了全程轨道控制设计方法。进行了包括地月转移、近月制动、环月降轨和动力下降的全程轨道控制的分段设计和联合规划，实现在入轨轨道偏差条件下的定时定点月面着陆。分别构建了中途修正、近月制动、环月降轨三段轨道控制的规划变量和目标参数；根据轨道倾角建立了动力下降点与着陆点的匹配转换关系。设计了中途修正、近月制动、环月降轨、动力下降的全程轨道控制策略的联合规划。建立了着陆位置偏差与轨道倾角偏差、着陆时间偏差与轨道半长轴偏差的修正关系，修正设计了中途修正目标倾角和近月制动目标半长轴。仿真算例表明，在入轨偏差轨道条件下，保证了中途修正后的飞行轨道与标称轨道基本一致，实现了与标称状态基本一致的定时定点月面着陆。可应用于月球着陆、月球采样返回以及载人登月等实施月面定时定点着陆任务的轨道设计和控制实施。     

Manifold dynamics in the Earth–Moon system via isomorphic mapping with application to spacecraft end-of-life strategies

Recently, manifold dynamics has assumed an increasing relevance for analysis and design of low-energy missions, both in the Earth–Moon system and in alternative multibody environments. With regard to lunar missions, exterior and interior transfers, based on the transit through the regions where the collinear libration points L₁ and L₂ are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is focused on the definition and use of a special isomorphic mapping for low-energy mission analysis. A convenient set of cylindrical coordinates is employed to describe the spacecraft dynamics (i.e. position and velocity), in the context of the circular restricted three-body problem, used to model the spacecraft motion in the Earth–Moon system. This isomorphic mapping of trajectories allows the identification and intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Heteroclinic connections, i.e. the trajectories that belong to both the stable and the unstable manifolds of two distinct periodic orbits, can be easily detected by means of this representation. This paper illustrates the use of isomorphic mapping for finding (a) periodic orbits, (b) heteroclinic connections between trajectories emanating from two Lyapunov orbits, the first at L₁, and the second at L₂, and (c) heteroclinic connections between trajectories emanating from the Lyapunov orbit at L₁ and from a particular unstable lunar orbit. Heteroclinic trajectories are asymptotic trajectories that travels at zero-propellant cost. In practical situations, a modest delta-v budget is required to perform transfers along the manifolds. This circumstance implies the possibility of performing complex missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining suitable, convenient end-of-life strategies for spacecraft orbiting the Earth. Seven distinct options are identified, and lead to placing the spacecraft into the final disposal orbit, which is either (a) a lunar capture orbit, (b) a lunar impact trajectory, (c) a stable lunar periodic orbit, or (d) an outer orbit, never approaching the Earth or the Moon. Two remarkable properties that relate the velocity variations with the spacecraft energy are employed for the purpose of identifying the optimal locations, magnitudes, and directions of the velocity impulses needed to perform the seven transfer trajectories. The overall performance of each end-of-life strategy is evaluated in terms of time of flight and propellant budget.