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

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

On solving the problems of optimization of trajectories of many-revolution orbit transfers of spacecraft

The problem of optimal control over many-revolution spacecraft orbit transfers between circular coplanar orbits of satellites is considered. The spacecraft flight is controlled by a thrust vector of a jet engine with restricted thrust (JERT). The mass expenditure is minimized at a limited time of flight. The optimal control problem is solved based on the maximum principle. The boundary value problem of the maximum principle is solved numerically using the shooting method. A modified computation scheme of the shooting method is suggested (multi-point shooting), as well as a method (correlated with the scheme) of choosing the initial approximation with the use of a solution to the optimization problem in the impulse formulation. The scheme and method allow one to construct many-revolution spacecraft orbit transfers.     

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.     

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.     

Conditions of the Maximum Principle in the Problem of Optimal Control over an Aggregate of Dynamic Systems and Their Application to Solution of the Problems of Optimal Control of Spacecraft Motion

The necessary first-order conditions of strong local optimality (conditions of maximum principle) are considered for the problems of optimal control over a set of dynamic systems. To derive them a method is suggested based on the Lagrange principle of removing constraints in the problems on a conditional extremum in a functional space. An algorithm of conversion from the problem of optimal control of an aggregate of dynamic systems to a multipoint boundary value problem is suggested for a set of systems of ordinary differential equations with the complete set of conditions necessary for its solution. An example of application of the methods and algorithm proposed is considered: the solution of the problem of constructing the trajectories of a spacecraft flight at a constant altitude above a preset area (or above a preset point) of a planet's surface in a vacuum (for a planet with atmosphere beyond the atmosphere). The spacecraft is launched from a certain circular orbit of a planet's satellite. This orbit is to be determined (optimized). Then the satellite is injected to the desired trajectory segment (or desired point) of a flyby above the planet's surface at a specified altitude. After the flyby the satellite is returned to the initial circular orbit. A method is proposed of correct accounting for constraints imposed on overload (mixed restrictions of inequality type) and on the distance from the planet center: extended (nonpointlike) intermediate (phase) restrictions of the equality type.     

The Use of a Solar Electrojet Propulsion System for Jupiter Satellite Injection

The problem of optimization of the interplanetary trajectory of flight for a multistage spacecraft with high- and low-thrust engines into the Jupiter satellite orbit is considered. Low-thrust engines (stationary plasma engines) are used on a heliocentric flight segment. Their operation is maintained with electric power supply from solar batteries. The principal feasibility of the realization of such a project is shown, and estimations of the mass of a spacecraft placed into Jupiter's satellite orbit are presented.     

Earth–Moon libration point orbit stationkeeping: Theory,modeling, and operations

Collinear Earth–Moon libration points have emerged as locations with immediate applications. These libration point orbits are inherently unstable and must be maintained regularly which constrains operations and maneuver locations. Stationkeeping is challenging due to relatively short time scales for divergence, effects of large orbital eccentricity of the secondary body, and third-body perturbations. Using the Acceleration Reconnection and Turbulence and Electrodynamics of the Moon's Interaction with the Sun (ARTEMIS) mission orbit as a platform, the fundamental behavior of the trajectories is explored using Poincaré maps in the circular restricted three-body problem. Operational stationkeeping results obtained using the Optimal Continuation Strategy are presented and compared to orbit stability information generated from mode analysis based in dynamical systems theory.     

一种新的共线平动点拟周期轨道稳定保持策略

限制性三体问题下共线平动点附近的拟周期轨道在深空探测中具有重要的实际应用价值,得到了各航天大国的广泛重视。通过将动力学中心流形结构引入轨道控制方法的设计之中,得到了基于投影到中心流形的共线平动点拟周期轨道稳定保持策略。首先推导了会合坐标到中心流形坐标的正则变换方法,在此基础上设法通过引入轨道机动,将偏差状态点投影到中心流形上,从而达到消除不稳定分量的目的。该方法充分整合了平动点的动力学特性,并且也适用于周期轨道的稳定保持。通过对Lissajous轨道和晕轨道的数值仿真表明,该方法较以往方法具有更强的稳定性,能在显著降低轨控燃料消耗的基础上达到较好的稳定保持效果。     

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

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

Two-craft tether formation relative equilibria about circular orbits and libration points

The relative equilibria of a two spacecraft tether formation connected by line-of-sight elastic forces moving in the context of a restricted two-body system and a circularly restricted three-body system are investigated. For a two spacecraft formation moving in a central gravitational field, a common assumption is that the center of the circular orbit is located at the primary mass and the center of mass of the formation orbits around the primary in a great-circle orbit. The relative equilibrium is called great-circle if the center of mass of the formation moves on the plane with the center of the gravitational field residing on it; otherwise, it is called a nongreat-circle orbit. Previous research shows that nongreat-circle equilibria in low Earth orbits exhibit a deflection of about a degree from the great-circle equilibria when spacecraft with unequal masses are separated by 350 km. This paper studies these equilibria (radial, along-track and orbit-normal in circular Earth orbit and Earth–Moon Libration points) for a range of inter-craft distances and semi-major axes of the formation center of mass. In the context of a two-spacecraft Coulomb formation with separation distances on the order of dozens of meters, this paper shows that the equilibria deflections are negligible (less than 10^?6°) even for very heterogeneous mass distributions. Furthermore, the nongreat-circle equilibria conditions for a two spacecraft tether structure at the Lagrangian libration points are developed.     

Optimization of multi-mode rendezvous trajectories with constraints

Approximate numerical methods of optimization of spacecraft rendezvous trajectories are presented that make use of interior point algorithms for problems of linear programming of high dimensionality (tens to hundreds of thousands of variables). The basis of the methods is discretization of a trajectory into small segments in which maneuvers are allowed to be executed; for all segments sets of pseudo-impulses are introduced that determine the possible directions of the spacecraft thrust vector. The terminal conditions are presented in the form of a linear matrix equation. A matrix inequality for the sums of characteristic velocities of pseudo-impulses on each segment is used to make a transformation to the linear programming form. Spacecraft rendezvous trajectories are considered in the neighborhood of circular orbits with the use of multi-mode propulsion systems (including those with low thrust) and existence of boundary conditions at interior points and constraints on the time of operation of the propulsion system at separate segments of the trajectory.     

The use of solutions to problems of spacecraft trajectory optimization in impulse formulation when solving the problems of optimal control of trajectories of a spacecraft with limited thrust engine: I

In this first part of our paper, it is suggested to use solutions to boundary value problems in the optimization problems (in impulse formulation) for spacecraft trajectories in order to obtain the initial approximation, when boundary value problems of the maximum principle are solved numerically by the shooting method. The technique suggested is applied to the problems of optimal control over motion of the center of mass of a spacecraft controlled by the thrust vector of jet engine with limited thrust in an arbitrary gravitational field in a vacuum. The method is based on a modified (in comparison to the classic scheme) shooting method computation together with the method of continuation along a parameter (maximum reactive acceleration, initial thrust-to-weight ratio, or any other parameter equivalent to them). This technique allows one to obtain the initial approximation with a high precision, and it is applicable to a wide range of optimal control problems solved using the maximum principle, if the impulse formulation makes sense for these problems.     

从地球到日-火系平动点轨道的转移

平动点轨道特殊的空间位置及动力学特征,使其在深空探测中具有重要的应用。以日-火系平动点轨道(Lissajous与Halo轨道)任务为目标,结合平动点轨道的不变流形理论,研究了小推力转移问题。首先给出了圆型限制性三体动力学模型下平动点附近不变流形(稳定和不稳定流形)高阶分析解以及相应的计算实例。接着以流形分析解为基础,建立了初始小推力轨道优化模型,并利用改进的协作进化算法求解初始小推力轨道。最后将初始轨道离散,采用多点打靶法将最优控制问题转化为参数优化问题,并用序列二次规划方法(SQP)求解。仿真结果证明轨道设计方法的有效性。     

Optimization of low-energy resonant hopping transfers between planetary moons

In response to the scientific interest in Jupiter's Galilean moons, NASA and ESA have plans to send orbiting missions to Europa and Ganymede, respectively. The inter-moon transfers of the Jovian system offer obvious advantages in terms of scientific return, but are also challenging to design and optimize due in part to the large, often chaotic, sensitivities associated with repeated close encounters of the planetary moons. The approach outlined in this paper confronts this shortcoming by exploiting the multi-body dynamics with a patched three-body model to enable multiple “resonant-hopping” gravity assists. Initial conditions of unstable resonant orbits are pre-computed and provide starting points for the elusive initial guess associated with the highly nonlinear optimization problem. The core of the optimization algorithm relies on a fast and robust multiple-shooting technique to provide better controllability and reduce the sensitivities associated with the close approach trajectories. The complexity of the optimization problem is also reduced with the help of the Tisserand–Poincaré (T–P) graph that provides a simple way to target trajectories in the patched three-body problem. Preliminary numerical results of inter-moon transfers in the Jovian system are presented. For example, using only 59 m/s and 158 days, a spacecraft can transfer between a close resonant orbit of Ganymede and a close resonant orbit of Europa.     

Possibilities of Combining High- and Low-Thrust Engines in Flights to Mars

The use of combinations of chemical and electric jet engines in the spacecraft designs results in a multistage vehicle configuration and in related problems of the optimum distribution of masses between the stages, the problems of flight trajectory optimization, and the problems of choosing the design parameters of a spacecraft. The appropriate issues are considered using flights to Mars as an example. The conditions for the optimum matching of high and low thrust trajectory segments are presented. The model of the simultaneous optimization of the geocentric and heliocentric legs of the trajectory is proposed. One- and two-orbit optimum trajectories of flight are investigated and analyzed.     

Guidance of a spacecraft with low lift-to-drag ratio to a landing point

The problem of terminal control over a deorbiting spacecraft at the stage of its flight after leaving plasma (altitude of ∼40 km) is considered, the aim being to guide it to a preset landing point. The algorithm is based on a modification of the well-known method of proportional navigation, when a fixed point is the target. It is suggested to use satellite navigation systems (of the GLONASS or GPS types) and/or radio beacons, which should allow one to determine the spacecraft trajectory parameters with high precision. Single-channel control is performed by changing the roll angle according to current parameters of the trajectory, which ensures adaptability of the method. Examples of three-dimensional trajectories of flight are presented for a manned spacecraft with low lift-to-drag ratio (∼0.5), currently under design in Russia. The results of statistical modeling taking into account initial deviations of the trajectory parameters and wind disturbances are presented. A method of statistical choice of a reference trajectory for the guidance stage is suggested. A theoretical possibility of using the algorithm of spacecraft guidance (in case of in-light accident with a carrier launcher) to preset regions in the vicinity of launching route is demonstrated. A qualitative analysis of proportional navigation with a fixed target is presented.     

The use of solutions to problems of spacecraft trajectory optimization in impulse formulation when solving the problems of optimal control of trajectories of a spacecraft with limited thrust engine: II

As examples of application of the technique suggested in the first part of this work, the problems of optimizing the trajectories of spacecraft transfers between circular coplanar orbits are considered in this second part. During the transfer the spacecraft is controlled by the vector of thrust of a limited-thrust jet engine. The mass consumption is minimized for a limited time of transfer. Extreme trajectories with two and three powered sections (Homan-type and bi-elliptic transfer trajectories) are numerically determined. The solution of these well-studied problems allows one to compare the results of applying the suggested technique with the results of application of other previously used techniques.     

Analysis of spacecraft motion under constant circumferential propulsive acceleration

This paper reassesses the classical circumferential-thrust problem, in which a spacecraft orbiting around a primary body is subjected to a propulsive acceleration of constant modulus, whose direction is in the plane of the parking orbit and orthogonal to the spacecraft-primary line. In particular, a new formulation is proposed to obtain a reduction in the number of differential equations required for the study of the spacecraft propelled trajectory. The mathematical complexity of the problem may be further reduced assuming that both the propulsive acceleration modulus and the spacecraft distance from the primary body are sufficiently small. In that case, an approximate model is able to accurately describe the characteristics of the propelled trajectory when the parking orbit is circular. Finally, using the data obtained by numerical simulations, the approximate model is extended to generate a set of semi-analytical equations for the analysis of a classical escape mission scenario.     

About one problem of trajectory optimization

The optimization problem for trajectories of spacecraft flight from the Earth to an asteroid is considered in this paper. The flight is realized in the central Newtonian gravitational field of the Sun with a possibility of gravitational maneuvers near planets. Perturbation maneuvers are taken into account using the method of point area of action with a limitation on the flyby altitude. The spacecraft is controlled by changing the value and direction of the engine thrust. The problem is solved taking into account constraints on the launch time, flight duration, and minimum distance to the Sun.     

The Use of Quaternions in the Optimal Control Problems of Motion of the Center of Mass of a Spacecraft in a Newtonian Gravitational Field: I

The problem of optimal control is considered for the motion of the center of mass of a spacecraft in a central Newtonian gravitational field. For solving the problem, two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems. Both the variants have a quaternion variable among the phase variables. In the first variant this variable characterizes the orientation of an instantaneous orbit of the spacecraft and (simultaneously) the spacecraft location in this orbit, while in the second variant only the instantaneous orbit orientation is specified by it. The suggested equations are convenient in the respect that they allow the general three-dimensional problem of optimal control by the motion of the spacecraft center of mass to be considered as a composition of two interrelated problems. In the first variant these problems are (1) the problem of control of the shape and size of the spacecraft orbit and (2) the problem of control of the orientation of a spacecraft orbit and the spacecraft location in this orbit. The second variant treats (1) the problem of control of the shape and size of the spacecraft orbit and the orbit location of the spacecraft and (2) the problem of control of the orientation of the spacecraft orbit. The use of quaternion variables makes this consideration most efficient. The problem of optimal control is solved on the basis of the maximum principle. Several first integrals of the systems of equations of the boundary value problems of the maximum principle are found. Transformations are suggested that reduce the dimensions of the systems of differential equations of boundary value problems (without complicating them). Geometrical interpretations are given to the transformations and first integrals. The relation of the vectorial first integral of one of the derived systems of equations (which is an analog of the well-known vectorial first integral of the studied problem of optimal control) with the found quaternion first integral is considered. In this paper, which is the first part of the work, we consider the models of motion of the spacecraft center of mass that employ quaternion variables. The problem of optimal control by the motion of the spacecraft center of mass is investigated on the basis of the first variant of equations of motion. An example of a numerical solution of the problem is given.