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.     

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.     

Earth–Moon low energy trajectory optimization in the real system

The problem of the Earth–Moon low energy trajectory optimization in the real system (the model defined by the JPL ephemeris DE405) is considered in this paper. First, this problem is investigated in the model of circular restricted three-body problem, since the fuel consumption is closely related to the Jacobi integral of the transfer trajectory, a method based on Jacobi integral is proposed and eight optimal trajectories are obtained. These optimal trajectories provide initial information (the flight time and the braking velocity impulse) to search the optimal low energy trajectories in the real system through optimization techniques. Considering the merit and drawback of particle swarm optimization and differential evolution algorithm in solving the space trajectory problem, an improved cooperative evolutionary algorithm is put forward. Result shows that the low energy trajectories in the real system are more fuel-efficient than the corresponding ones under the circular restricted three-body problem.     

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.     

Multiobjective genetic optimization of Earth–Moon trajectories in the restricted four-body problem

In this paper, optimal trajectories of a spacecraft traveling from Earth to Moon using impulsive maneuvers (ΔV maneuvers) are investigated. The total flight time and the summation of impulsive maneuvers ΔV are the objective functions to be minimized. The main celestial bodies influencing the motion of the spacecraft in this journey are Sun, Earth and Moon. Therefore, a three-dimensional restricted four-body problem (R4BP) model is utilized to represent the motion of the spacecraft in the gravitational field of these celestial bodies. The total ΔV of the maneuvers is minimized by eliminating the ΔV required for capturing the spacecraft by Moon. In this regard, only a mid-course impulsive maneuver is utilized for Moon ballistic capture. To achieve such trajectories, the optimization problem is parameterized with respect to the orbital elements of the ballistic capture orbits around Moon, the arrival date and a mid-course maneuver time. The equations of motion are solved backward in time with three impulsive maneuvers up to a specified low Earth parking orbit. The results show high potential and capability of this type of parameterization in finding several Pareto-optimal trajectories. Using the non-dominated sorting genetic algorithm with crowding distance sorting (NSGA-II) for the resulting multiobjective optimization problem, several trajectories are discovered. The resulting trajectories of the presented scheme permit alternative trade-off studies by designers incorporating higher level information and mission priorities.     

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.     

考虑月球升交点黄经的地月低能转移轨道

地月转移轨道研究的基本模型主要是限制性二体模型和限制性三体模型。在限制性三体建模下,基于轨道分类理论的地月脉冲变轨比传统的霍曼变轨节省了更多的燃料(至少13%以上)。单独的"地球-月球-飞行器"的限制性三体问题无法应用轨道分类的理论进行轨道设计,必须加入太阳。文章考虑了月球升交点赤经、黄白交角,以及正确的地月系统稳定流形的方向,整个"太阳-地球-月球-飞行器"系统拆分为"地球-月球-飞行器"的平面圆型限制性三体问题,和"太阳-地球-飞行器"三维圆型限制性三体问题进行分析,并利用Matlab仿真验证了方案的可行性。     

Optimal station-keeping near Earth–Moon collinear libration points using continuous and impulsive maneuvers

In this study the gravitational perturbations of the Sun and other planets are modeled on the dynamics near the Earth–Moon Lagrange points and optimal continuous and discrete station-keeping maneuvers are found to maintain spacecraft about these points. The most critical perturbation effect near the L₁ and L₂ Lagrange points of the Earth–Moon is the ellipticity of the Moon’s orbit and the Sun’s gravity, respectively. These perturbations deviate the spacecraft from its nominal orbit and have been modeled through a restricted five-body problem (R5BP) formulation compatible with circular restricted three-body problem (CR3BP). The continuous control or impulsive maneuvers can compensate the deviation and keep the spacecraft on the closed orbit about the Lagrange point. The continuous control has been computed using linear quadratic regulator (LQR) and is compared with nonlinear programming (NP). The multiple shooting (MS) has been used for the computation of impulsive maneuvers to keep the trajectory closed and subsequently an optimized MS (OMS) method and multiple impulses optimization (MIO) method have been introduced, which minimize the summation of multiple impulses. In these two methods the spacecraft is allowed to deviate from the nominal orbit; however, the spacecraft trajectory should close itself. In this manner, some closed or nearly closed trajectories around the Earth–Moon Lagrange points are found that need almost zero station-keeping maneuver.     

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.     

Simulation software for the spiral trajectory of our Moon

At 4.56 Ga, the accretion of the slowly rotating Solar Nebula led to the formation of Sun and its Planets in the plane of disc of accretion. Moon was formed by accretion from a circumterrestrial disk of debris generated by the glancing angle impact of the young Earth by a Mars size planetary embryo at about 4.5 Ga at a distance of 15,000 km. The Moon since then has migrated to the present position of 384,400 km from the center of the Earth. In course of this outward migration it has slowed down the spin rate of Earth and caused the lengthening of diurnal day length from 5 h initially to 24 h presently. The basic mechanics of Earth–Moon System has been worked out and theoretical determination of lengthening of day curve is carried out. This theoretical lengthening of day curve is compared with the observed lengthening of day curve based on paleobotanical evidences, ancient tidalites and Australian Banded Iron Formation. There is a remarkable correspondence between the two curves except for intermittent deviations due to geographical and geophysical factors. Based on the theoretical curve of lengthening of day, an empirical formula for the lunar orbital radius expansion is determined. Based on this empirical formula, simulation software is developed that gives the correct evolution of the semi-major axis (a) of our Moon for any time span from the inception to the time chosen under study. For mathematical simplicity the system is considered to be a two body rotating system throughout its evolutionary history of 4.5 Gyrs. This simulation draws the Moon’s spiral trajectory from its inception to any subsequent epoch. The terminal epoch is an input to the simulation software to arrive at the spiral trajectory of the Moon from the inception to the given epoch. The basic mechanics of Earth–Moon System and this simulation can be generalized to lay the foundation of simulation software for any Planet–Satellite pair or any Sun-Planet pair in our Solar System or Star-Planet pair in any Extra-Solar System. The basic dynamics has been found to be valid for Star–Planet pair also. So this Simulation Methodology can as well be applied to study the migratory evolution of Gas Giants also.     

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.     

Resonant motion, ballistic escape, and their applications in astrodynamics

A special set of solutions governing the motion of a particle, subject to the gravitational attractions of the Earth, the Moon, and, eventually, the Sun, is discussed in this paper. These solutions, called resonant orbits, correspond to a special motion where the particle is in resonance with the Moon. For a restricted set of initial conditions the particle performs a resonance transition in the vicinity of the Moon. In this paper, the nature of the resonance transition is investigated under the perspective of the dynamical system theory and the energy approach. In particular, using a new definition of weak stability boundary, we show that the resonance transition mechanism is strictly related to the concept of weak capture. This is shown through a carefully computed set of Poincaré surfaces, at different energy levels, on which both the weak stability boundary and the resonant orbits are represented. It is numerically demonstrated that resonance transitioning orbits pass through the weak stability boundaries. In the second part of the paper the solar perturbation is taken into account, and the motion of the resonant orbits is studied within a four-body dynamics. We show that, for a wide class of initial conditions, the particle escapes from the Earth–Moon system and targets an heliocentric orbit. This is a free ejection called a ballistic escape. Astrodynamical applications are discussed.     

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.     

Low-thrust transfer to Backflip orbits

The aim of the work is to design a low-thrust transfer from a Low Earth Orbit to a “useful” periodic orbit in the Earth–Moon Circular Restricted Three Body Model (CR3BP). A useful periodic orbit is here intended as one that moves both in the Earth–Moon plane and out of this plane without any requirements of propellant mass. This is achieved by exploiting a particular class of periodic orbits named Backflip orbits, enabled by the CR3BP. The unique characteristics of this class of periodic solutions allow the design of an almost planar transfer from a geocentric orbit and the use of the Backflip intrinsic characteristics to explore the geospace out of the Earth–Moon plane. The main advantage of this approach is that periodic plane changes can be obtained by performing an almost planar transfer. In order to save propellant mass, so as to increase the scientific payload of the mission, a low-powered transfer is considered. This foresees a thrusting phase to gain energy from a departing circular geocentric orbit and a second thrusting phase to match the state of the target Backflip orbit, separated by an intermediate ballistic phase. This results in a combined application of a low-thrust manoeuvre and of a periodical solution in the CR3BP to realize a new class of missions to explore the Earth–Moon neighbourhoods in a quite inexpensive way. In addition, a low-thrust transit between two different Backflip orbits is analyzed and considered as a possible extension of the proposed mission. Thus, also a Backflip-to-Backflip transfer is addressed where a low-powered probe is able to experience periodic excursions above and below the Earth–Moon plane only performing almost planar and very short transfers.     

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

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

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.     

Libration point orbit rendezvous using PWPF modulated terminal sliding mode control

The near-range rendezvous problem of two libration point orbit spacecraft in the Earth–Moon system is studied using the terminal sliding mode control which enables a time-fixed process with the flight time prescribed a priori. The underlying dynamics are the full nonlinear equations of motion for a complete Solar System model. For practical purposes, two means of pulse-width pulse-frequency (PWPF) modulation are employed to realize the theoretical continuous control with a series of thrust pulses. Extensive simulations with major errors taken into account show that the sliding mode controller can successfully guide the chaser to a given staging node with the final position and velocity errors, on average, lower than 20 m and 1 mm/s, respectively. Compared with the glideslope guidance previously studied, the proposed approach outperforms the former by saving approximately 50–60% of total delta-v.     

Optimal guidance based on receding horizon control for low-thrust transfer to libration point orbits

This paper addresses the design and computation of a guidance law for a transfer mission from an orbit near the Earth to a halo orbit around the libration point L₂ in the Sun–Earth system. The guidance law, which is designed based on receding horizon control and compensates for launch velocity errors that are introduced by inaccuracies of the launch vehicle, is solved using the generating function method. During the design of the closed-loop guidance law, the entire transfer mission, which is considered a nonlinear optimal control problem, is evaluated to obtain a nominal reference trajectory. Using the launch velocity errors and the uncertainty of the model, a spacecraft controlled by the proposed guidance law tracks the reference trajectory. Furthermore, the original Riccati differential equation in the receding horizon control algorithm is replaced by an equivalent convenient form of the Riccati differential equation that is based on the generating function. The high-efficiency solution of the equivalent equation avoids the online direct integration of the original Riccati differential equation, which significantly increases the computational efficiency for the receding horizon control problem. Numerical simulations using a nonlinear bicircular four-body model demonstrate the capabilities of the proposed receding horizon guidance law for the transfer mission. In addition, the generating function method improves the computational efficiency by at least one order of magnitude over the backward sweep method in solving the receding horizon control problem.     

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

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