Alternative transfer to the Earth–Moon Lagrangian points L4 and L5 using lunar gravity assist

Lagrangian points L4 and L5 lie at 60° ahead of and behind the Moon in its orbit with respect to the Earth. Each one of them is a third point of an equilateral triangle with the base of the line defined by those two bodies. These Lagrangian points are stable for the Earth–Moon mass ratio. As so, these Lagrangian points represent remarkable positions to host astronomical observatories or space stations. However, this same distance characteristic may be a challenge for periodic servicing mission. This paper studies elliptic trajectories from an Earth circular parking orbit to reach the Moon’s sphere of influence and apply a swing-by maneuver in order to re-direct the path of a spacecraft to a vicinity of the Lagrangian points L4 and L5. Once the geocentric transfer orbit and the initial impulsive thrust have been determined, the goal is to establish the angle at which the geocentric trajectory crosses the lunar sphere of influence in such a way that when the spacecraft leaves the Moon’s gravitational field, its trajectory and velocity with respect to the Earth change in order to the spacecraft arrives at L4 and L5. In this work, the planar Circular Restricted Three Body Problem approximation is used and in order to avoid solving a two boundary problem, the patched-conic approximation is considered.     

Design and analysis of an extended mission of CE-2: From lunar orbit to Sun–Earth L2 region

Chang’E-2 (CE-2) has firstly successfully achieved the exploring mission from lunar orbit to Sun–Earth L2 region. In this paper, we discuss the design problem of transfer trajectory and at the same time analyze the visible segment of Tracking, Telemetry & Control (TT&C) system for this mission. Firstly, the four-body problem of Sun–Earth–Moon and Spacecraft can be decoupled in two different three-body problems (Sun–Earth + Moon Restricted Three-Body Problems (RTBPs) and Earth–Moon ephemeris model). Then, the transfer trajectory segments in different model are computed, respectively, and patched by Poincaré sections. The full-flight trajectory including transfer trajectory from lunar orbit to Sun–Earth L2 region and target Lissajous orbit is obtained by the differential correction method. Finally, the visibility of TT&C system at the key time is analyzed. Actual execution of CE-2 extended mission shows that the trajectory design of CE-2 mission is feasible.     

Two-manoeuvres transfers between LEOs and Lissajous orbits in the Earth–Moon system

The purpose of this work is to compute transfer trajectories from a given Low Earth Orbit (LEO) to a nominal Lissajous quasi-periodic orbit either around the point L₁ or the point L₂ in the Earth–Moon system. This is achieved by adopting the Circular Restricted Three-Body Problem (CR3BP) as force model and applying the tools of Dynamical Systems Theory.     

地月空间的远距离逆行轨道族及其分岔研究

地月系统中存在着一类绕月逆行、高度稳定的轨道族,称为远距离逆行轨道族(DRO)。以圆型限制性三体问题(CR3BP)为动力学模型研究了DRO轨道族周边的动力系统结构。利用Broucke稳定性图寻找分叉点,判断分叉类型,基于数值延拓计算分岔后产生的一系列新轨道分支。分叉类型主要有切分叉与多倍周期分叉(从3倍周期开始),轨道维度包含平面轨道族与三维轨道族。计算新轨道族的特征,包括形状、周期、能量、稳定性、双曲流形结构等。探讨周期轨道的轨道周期与能量的关系,以几何化的方式展现分叉结构、多周期轨道的双曲流形结构等。该动力结构将为基于DRO轨道族的地月空间任务提供重要的理论支持。      

Evolution of periodic orbits near the Lagrangian point L2

A study of the evolution of the periodic and the quasi-periodic orbits near the Lagrangian point L₂, which is located to the right of the smaller primary on the line joining the primaries and whose distance from the more massive primary is greater than the distance between the primaries, in the framework of restricted three-body problem for the Sun–Jupiter, Earth–Moon (relatively large mass ratio) and Saturn–Titan (relatively small mass ratio) systems is made. Two families of periodic orbits around the smaller primary are identified using the Poincaré surface of section method – family I (initially elliptical, gradually becomes egg-shaped with the increase in the Jacobi constant C and elongated towards the more massive primary) and family II (initially egg-shaped orbits elongated towards L₂ and gradually becomes elliptical with the increase in C). The family I in the Sun–Jupiter and Saturn–Titan systems contains two separatrix caused by third-order and fourth-order resonances, while the Earth–Moon system has only one separatrix which is caused by third-order resonances. Also in the Sun–Jupiter and the Saturn–Titan systems, family I merge with family II, around Jacobian constant 3.0393 and 3.0163, respectively, while in the Earth–Moon system, family II evolves separately from two different branches. The two branches merge at C = 3.184515. In the Earth–Moon system, the family II contains a separatrix due to third-order resonances which is absent in the other two systems.     

平动点周期轨道间小推力转移的Gauss伪谱法

针对三体问题共线平动点附近周期轨道间的小推力转移问题,构造了一种新的形状函数,在此基础上提出了一种基于Gauss伪谱法的优化设计方法。首先,建立小推力轨道转移动力学模型,参考初始轨道和目标轨道的类型,构造一种新的形状函数以近似小推力转移轨道。为满足不同的约束要求,提出了振幅和相位按多项式变化的假设,推导了小推力转移轨道的近似解析解;然后利用Gauss伪谱法将小推力轨道转移的最优控制问题转化为非线性规划问题,并对推导的近似解析解进行解算和处理,为Gauss伪谱法求解非线性规划问题提供较为有效的控制变量的初始猜测值;最后以地月系统L1点附近Halo轨道间的小推力转移问题为例进行了仿真分析。仿真结果表明,小推力转移轨道近似解析解具备有效性和普适性,使得Gauss伪谱法的迭代效率提高55%以上,同时也表明Gauss伪谱法可有效解决平动点周期轨道间的小推力转移轨道优化设计问题。     

从日地系统L2出发借力月球飞越近地小行星

对于停留在日地系统L2的“嫦娥2号”探测器,其后续飞行方案有多个选项,例如主动撞月或重返月球轨道、返回地球轨道或再入大气、飞往地月系统L1/L2或日地系统L1、进入深空飞越近地小行星(最终,“嫦娥2号”于2012年12月13日成功地实现了对Toutatis小行星的近距离飞越)。探讨上述的飞行方案需要对飞行轨道进行初步设计,总的速度脉冲限制在100 m/s以内并且需要考虑探测器同时受到太阳、地球、月球的引力作用。本研究设计了探测器从日地系统L2出发借力月球实现Toutatis小行星飞越的飞行方案,与直接飞越方案相比,借力月球可以进一步节省探测器的燃料消耗,其等效速度脉冲设计值为58.47 m/s。     

高精度模型下Halo轨道设计研究

针对未来地月L2点Halo轨道空间站长期停泊任务,研究了高精度模型下Halo轨道设计方法。首先,详细推导了圆型限制性三体问题(Circular Restricted Three Body Problem,CR3BP)质心会合坐标系与高精度模型地心J2000坐标系之间的转换关系,并在此基础上,将CR3BP下的闭合Halo轨道转换到地心J2000坐标系得到了高精度模型下Halo轨道迭代初值。其次,采用序列二次规划(Sequence Quadratic Program,SQP)构造多层迭代格式,在高精度模型下对初值进行逐层修正。最后,通过仿真测试,验证了该方法的可行性与有效性。该研究结果可为未来平动点任务标称轨道设计方案的制定提供参考。     

Libration point orbits for lunar global positioning systems

With the development of lunar exploration, a lunar global positioning system (LGPS) is demanded for both on-ground and in-flight lunar exploration missions. The traditional configuration of constellation requires at least eighteen satellites to cover the whole lunar surface continuously. In this paper, the configurations of the libration point orbits (LPOs) constellations are investigated. By using the constellations on the Earth–Moon _L1$L_1$ and _L2$L_2$ LPOs, the basic functions of the LGPS can be realized by using eight to fourteen satellites. First, the LPO and the combinations of LPOs, which can be used in the constellations of the LGPS, are investigated. The criteria and procedures of the configuration design are introduced. Second, the configurations of LPOs constellations are investigated in the Earth–Moon circular-restricted three-body problem (CR3BP). The size of the LPOs and the distribution of the satellites on these LPOs are determined by using an exhaustive algorithm and a global optimization method, respectively. The key performance parameters of these constellations are computed. Third, the constellations with good performance in the CR3BP are redesigned in the more accurate Earth–Moon based Sun-perturbed bicircular four-body problem (B4BP). Moreover, in order to avoid the ground coverage problem caused by the perturbation of the Sun, some modifications are implemented, and the configuration of the no blind area LGPS in the B4BP is obtained.     

Periodic attitude motions along planar orbits in the elliptic restricted three-body problem

A way to improve the accuracy of the three-body problem model is taking into account the eccentricity of primary attractors. Elliptic Restricted Three-Body Problem (ER3BP) is a model for studying spacecraft trajectory within the three-body problem such that the orbital eccentricity of primaries is reflected in it. As the principal cause of perturbation in the employed dynamical model, the primaries eccentricity changes the structure of orbits compared to the ideal Circular Restricted Three-Body Problem (CR3BP). It also changes the attitude behavior of a spacecraft revolving along periodic orbits in this regime. In this paper, the coupled orbit-attitude dynamics of a spacecraft in the ER3BP are exploited to find precise periodic solutions as the spacecraft is considered to be in planar orbits around Lagrangian points and Distant Retrograde Orbits (DRO). Periodic solutions are repetitious behaviors in which spacecraft whole dynamics are repeated periodically, these periodic behaviors are the main interest of this study because they are beneficial for future mission designs and allow delineation of the system’s governing dynamics. Previous studies laid the foundation for spacecraft stability analysis or studying pitch motion of spacecraft in the ER3BP regime. While in this paper, at first, initial guesses for correction algorithms were derived through verified search methods, then correction algorithms were used to refine calculated orbit-attitude periodic behaviors. Periodic orbits and full periodic solutions are portrayed and compared to previous studies and simpler models. Natural periodic solutions are valuable information eventuate in the longer functional lifetime of spacecraft. Since the problem assumption considered in this paper is much closer to real mission conditions, these results may be the means to use natural bounded motions in the actual operational environment.     

Multiple Sun–Earth Saddle Point flybys for LISA Pathfinder

LISA Pathfinder is an ESA mission due to be launched in the next two years. The gravity gradiometer onboard has the sensitivity required to test predictions by gravitational theories proposed as alternatives to Dark Matter such as TeVeS. Within the Solar System measurable effects are predicted only in the vicinity of gravitational saddle points (SP). For this reason it has been proposed to fly LPF by the Earth–Sun SP, at some 259,000 km from Earth. This could be done in an extension to the nominal mission which uses a Lissajous orbit about the Earth–Sun L1 point. The responsibility for LPF mission design lies with ESA/ESOC, who have designed the transfer trajectories, orbits about L1, and station keeping strategies. This article describes an analysis performed by Astrium to support a suggestion for a possible mission extension to a saddle point crossing. With only very limited fuel availability, reaching the saddle point is a significant challenge. In this article, we present recent advances in the work on trajectory design. It is demonstrated that reaching the SP is feasible once the LPF mission is completed. Furthermore, in a significant enhancement, it is demonstrated that trajectories including more than one SP flyby are possible, thus improving the science return for this proposed mission extension.     

Families of halo orbits in the elliptic restricted three-body problem for a solar sail with reflectivity control devices

Solar sail halo orbits designed in the Sun-Earth circular restricted three-body problem (CR3BP) provide inefficient reference orbits for station-keeping since the disturbance due to the eccentricity of the Earth’s orbit has to be compensated for. This paper presents a strategy to compute families of halo orbits around the collinear artificial equilibrium points in the Sun-Earth elliptic restricted three-body problem (ER3BP) for a solar sail with reflectivity control devices (RCDs). In this non-autonomous model, periodic halo orbits only exist when their periods are equal to integer multiples of one year. Here multi-revolution halo orbits with periods equal to integer multiples of one year are constructed in the CR3BP and then used as seeds to numerically continue the halo orbits in the ER3BP. The linear stability of the orbits is analyzed which shows that the in-plane motion is unstable while the out-of-plane motion is neutrally stable and a bifurcation is identified. Finally, station-keeping is performed which shows that a reference orbit designed in the ER3BP is significantly more efficient than that designed in the CR3BP, while the addition of RCDs improve station-keeping performance and robustness to uncertainty in the sail lightness number.     

A lunar cargo mission design strategy using variable low thrust

A design technique for a near optimal, Earth–Moon transfer trajectory using continuous variable low thrust is proposed. For the Earth–Moon transfer trajectory, analytical and numerical methods are combined to formulate the trajectory optimization problem. The basic concept of the proposed technique is to utilize analytically optimized solutions when the spacecraft is flying near a central body where the transfer trajectories are nearly circular shaped, and to use a numerical optimization method to match the spacecraft’s states to establish a final near optimal trajectory. The plasma thruster is considered as the main propulsion system which is currently being developed for crewed/cargo missions for interplanetary flight. The gravitational effects of the 3rd body and geopotential effects are included during the trajectory optimization process. With the proposed design technique, Earth–Moon transfer trajectory is successfully designed with the plasma thruster having a thrust direction sequence of “fixed-varied-fixed” and a thrust acceleration sequence of “constant-variable-constant”. As this strategy has the characteristics of a lesser computational load, little sensitivity to initial conditions, and obtaining solutions quickly, this method can be utilized in the initial scoping studies for mission design and analysis. Additionally, derived near optimal trajectory solution can be used as for initial trajectory solution for further detailed optimization problem. The demonstrated results will give various insights into future lunar cargo trajectories using plasma thrusters with continuous variable low thrust, establishing approximate costs as well as trajectory characteristics.     

Hierarchical global search for fuel-optimal impulsive transfers between lunar libration orbits

The attention to the periodic orbit in the Earth-Moon restricted three-body system continues to grow due to its special environment and locations. This research investigates the feasibility of constructing fuel-optimal single and multiple impulse transfers between unstable periodic orbits at L₁ and L₂ points. Invariant manifolds, which could provide the appropriate initial trajectories for optimization, are analyzed deeply to enable previously unknown orbit options and potentially to reduce mission cost. A global search strategy based on comparing the orbital state of the unstable and stable manifolds, incorporated with low-thrust techniques, is performed to seek a suitable matching point for maneuver application. Then the sequential quadratic programming (SQP) is adopted to further optimize the velocity increment and obtain the single/multiple impulse optimal transfers. The associated constraint gradients are derived to achieve higher accuracy and rapidity of the algorithm. To highlight the effectivity of the transfer scheme, three-dimensional low-energy transfers between different types and spatial regions of performing single and multiple impulses are explored. The total Delta-V required varies between a few meters per second and tens of meters per second, and the related flight time is about several weeks, mainly depending on the energy of periodic orbits and the invariant manifold structure. The results obtained in this paper can provide a useful reference for the selection of escape and capture site along the manifolds, maneuver magnitude and transfer time.     

Resonant quasi-periodic near-rectilinear Halo orbits in the Elliptic-Circular Earth-Moon-Sun Problem

Motivated by the near-future re-exploration of the cislunar space, this paper investigates dynamical substitutes of the Earth-Moon’s resonant Near-Rectilinear Halo Orbits (NRHOs) under the Elliptic-Circular Restricted Four-Body Problem formulation of the Earth-Moon-Sun system. This model considers that the Earth and Moon move in elliptical orbits about each other and that a third body, the Sun, moves in a circular orbit about the Earth-Moon barycenter. By making use of this higher-fidelity dynamical model, we are able to incorporate the Sun’s influence and the Moon’s eccentricity, two of the most significant perturbations of the cislunar environment. As a result of these perturbations, resonant periodic NRHOs of the Earth-Moon Circular Restricted Three-Body Problem (CR3BP) are hereby replaced by two-dimensional quasi-periodic tori that better represent the dynamical evolution of satellites near the vicinity of the Moon. We present the steps and algorithms needed to compute these dynamical structures in the Elliptic-Circular model and subsequently assess their utility for spacecraft missions. We focus on the planned orbit for the NASA-led Lunar Gateway mission, a 9:2 synodic resonant L₂ southern NRHO, as well as on the 4:1 synodic and 4:1 sidereal resonances, due to the proximity to the nominal orbit and their advantageous dynamical properties. We verify that the dynamical equivalents of these orbits preserve key dynamical attributes such as eclipse avoidance and near-linear stability. Furthermore, we find that the higher dimensionality of quasi-periodic solutions offers interesting alternatives to mission designers in terms of phasing maneuvers and low-altitude scientific observations.     

Properties of transit trajectory in the restricted three and four-body problem

In the restricted three-body problem if the Jacobi constant is just below the value corresponding to Lagrangian point only a little neck exists around the equilibrium point and capture trajectories are indicated as low-energy. Capture properties depend on the dynamics around these critical points and qualitative results can be obtained using linearized systems. In this paper, to study transit trajectory properties in the restricted three and four-body problem, the Earth–Moon–Sun–Satellite system is considered as example and studied using different models. In the circular restricted three-body problem (Earth–Moon–Satellite), transit, non transit and asymptotic trajectories, are easily identified by using the principal reference frame. Dynamics around Lagrangian point are then studied introducing the Moon eccentricity into the elliptical restricted three-body model. A preferential region for transit orbit is individuated and studied as a function of eigenvalue properties. To introduce the Sun effect, the bi-circular four-body model is considered and dynamics around Lagrangian points studied as a function of angular distance between Earth–Sun and Earth–Moon line. Finally, results obtained in the elliptical three-body model and bi-circular four-body model, are compared with numerical simulations using real Sun–Moon–Earth ephemeris.     

火卫一周期准卫星轨道及入轨分析

围绕火卫一的准卫星轨道（QSOs）因其具有良好的稳定性,是火卫一探测任务最为实用的轨道。在平面圆型限制性三体问题模型下,利用庞加莱截面和KAM环迭代方法探究了准卫星轨道的周期轨道族,并给出不同能量准卫星周期轨道的初始条件。针对火卫一周期准卫星轨道入轨,提出一种转移轨道设计方法:对准卫星周期轨道调整速度后进行反向积分,直至离开火卫一邻近区域,从而得到由火星环绕轨道向火卫一周期准卫星轨道的转移轨道,并调整转移轨道参数对燃料与时间消耗进行优化。研究结果表明,当周期准卫星轨道能量处于特定区间时,存在特定速度脉冲区间,可利用火卫一引力实现较少燃料消耗的轨道转移;在该速度脉冲区间中,通过选取较小的速度脉冲,可缩短转移时间。      

基于三体系统平动点轨道的环日全景任务轨道设计

在空间开展太阳观测是研究太阳活动周、太阳爆发、极端天气等事件起源的重要手段。环日全景探测计划是为实现从黄道面360°全方位观察太阳行星际空间而提出的。本文针对环日全景探测计划,构建了基于三体系统平动点低能量轨道的环日全景轨道部署方法。该方法以日–地L1/L2点Halo轨道幅值及Halo轨道离轨点为变量,以转移轨道飞行时间、入轨机动大小为评价指标,基于三体系统不变流形构建环日全景的转移轨道,并开展轨道优化设计。采用等高线图对设计变量及任务成本进行全局分析。仿真计算表明,轨道部署无法同时满足飞行时间最短与入轨机动最小的要求。设计了轨道机动约束条件下的最优飞行时间解,并给出了基于长三甲运载火箭的一箭双星发射及入轨方案。      

Mapping orbits around the asteroid 2001SN263

The present paper has the goal of mapping orbits, with respect to the perturbations, for a spacecraft traveling around the asteroid 2001SN₂₆₃. This asteroid is a triple system, which center of mass is in an elliptic orbit around the Sun. The perturbations considered in the present model are the ones due to the oblateness of the central body, the gravity field of the two satellite bodies (Beta and Gamma), the Sun, the Moon, the asteroids Vesta, Pallas and Ceres and all the planets of the Solar System. This mapping is important, because it shows the relative importance of each force for a given orbit for the spacecraft, helping to make a decision about which forces need to be included in the model for a given accuracy and nominal orbit. Another important application of this type of mapping is to find orbits that are less perturbed, since it is expected that those orbits have good potential to require a smaller number of station-keeping maneuvers. Simulations under different conditions are made to find those orbits. The main reason to study those trajectories is that, currently, there are several institutions in Brazil studying the possibility to make a mission to send a spacecraft to this asteroid (the so-called ASTER mission), because there are many important scientific studies that can be performed in that system. The results showed that Gamma is the main perturbing body, followed by Beta (10 times smaller) and the group Sun–Mars-oblateness of Alpha, with perturbations 1000 times weaker than the effects of Gamma. The other bodies have perturbations 10⁷ times smaller. The results also showed that circular and polar orbits are less perturbed, when compared to elliptical and equatorial orbits. Regarding the semi-major axis, an internal orbit is the best choice, followed by a larger external orbit. The inclination of the orbit plays an important role, and there are values for the inclination where the perturbations show minimum and maximum values, so it is important to make a good decision on those values.     

Dynamics of a spacecraft and normalization around Lagrangian points in the Neptune–Triton system

The problem of a spacecraft orbiting the Neptune–Triton system is presented. The new ingredients in this restricted three body problem are the Neptune oblateness and the high inclined and retrograde motion of Triton. First we present some interesting simulations showing the role played by the oblateness on a Neptune’s satellite, disturbed by Triton. We also give an extensive numerical exploration in the case when the spacecraft orbits Triton, considering Sun, Neptune and its planetary oblateness as disturbers. In the plane a × I (a = semi-major axis, I = inclination), we give a plot of the stable regions where the massless body can survive for thousand of years. Retrograde and direct orbits were considered and as usual, the region of stability is much more significant for the case of direct orbit of the spacecraft (Triton’s orbit is retrograde). Next we explore the dynamics in a vicinity of the Lagrangian points. The Birkhoff normalization is constructed around L₂, followed by its reduction to the center manifold. In this reduced dynamics, a convenient Poincaré section shows the interplay of the Lyapunov and halo periodic orbits, Lissajous and quasi-halo tori as well as the stable and unstable manifolds of the planar Lyapunov orbit. To show the effect of the oblateness, the planar Lyapunov family emanating from the Lagrangian points and three-dimensional halo orbits are obtained by the numerical continuation method.