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.     

Periodic motions of a nearly dynamically symmetric satellite in the neighborhood of hyperboloidal precession

The motion of a satellite close to a dynamically symmetric solid body in a Newtonian gravitational field over a circular orbit is studied. The system of differential equations describing the body’s motion is close to a system with cyclic coordinate. New classes of periodic motions are constructed in the neighborhood of a known partial solution to an unperturbed problem, hyperboloidal precession of a dynamically symmetric satellite. In the resonance case, when the ratio of one frequency of small oscillations of a reduced system with two degrees of freedom in the neighborhood of a stable equilibrium position to the frequency of cyclic coordinate variation is close to an integer number, there exist one or three families of periodic motions that are analytic in terms of fractional powers of a small parameter. A study of stability of these motions was performed with the help of KAM (Kolmogorov-Arnold-Moser) theoty. Faling the described resonance there exists a unique family of periodic motions that is analytic in terms of integer powers of a small parameter. The check-up of stability of these motrons was carried out. We distinguished the cases of parametric resonance, resonances of the third and fourth orders, and a non-resonant case. In the resonance cases our study relies on well-known results on stability of Hamiltonian systems during resonances [1]. In the non-resonant case we use the KAM theory [2].     

The family of potentials admitting integration of the perturbed two-body problem in regular coordinates

An infinite system of potentials is presented that admits separation of regular variables in the perturbed two-body problem. The regular coordinates are constructed using a specially selected L-matrix. An explicit solution to the problem in the elliptical case is constructed. In the general case the solution is reduced to inversion of hyper-elliptic integrals. The cases of motion with and without constraints are considered. The results of numerical experiments are presented.     

Dynamics and chaos control of asymmetric gyrostat satellites

The motion of a free gyrostat consisting of a platform with a triaxial ellipsoid of inertia and a rotor with a slight asymmetry with respect to the axis of rotation is considered. Dimensionless equations of motion for a system with perturbations caused by the small asymmetries of the rotor are written in Andoyer-Deprit variables. These perturbations result in a chaotic layer in the separatrix vicinity. Heteroclinic and homoclinic trajectories are written in analytical form for gyrostats with different ratios of their moments of inertia. These trajectories are used to construct a modified Melnikov function, and to produce control that eliminates separatrix chaos. The Poincare sections and Melnikov function are constructed via numerical modeling that demonstrates the effectiveness of control.     

Through optimization of branching injection trajectories by the Pontryagin maximum principle using stochasticmodels

A problem of the through optimization of aerospace vehicle branching ascent trajectories in view of atmospheric disturbances (wind and thermodynamic parameters variations) and constraints on admissible areas of separated parts fall points is considered. The problem solution is based on the Pontryagin maximum principle. The problem with a continuum of branches corresponding random disturbances was reduced to the more standard one with a finite set of branches being the most critical (at the given event probability) with regard to the functional considered. The analytical synthesis of such critical disturbances profiles is obtained on basis of the Bliss formula, connecting a functional change with variations of right parts of motion equations. The developed technique and its efficiency are demonstrated by an example for a typical launcher.     

Optimal control of reorientation of a spacecraft using free trajectory method

The problem of optimal control over spatial reorientation of a spacecraft is considered. The functional having a sense of propellant consumption is minimized. The analytical solution to the formulated problem is presented. It is shown that the optimal solution can be found in the class of two-impulse control at which the spacecraft’s turn is performed along a free motion trajectory. In order to improve the accuracy of spacecraft guidance into a specified angular position, methods of control are suggested that realize the method of free trajectories. The synthesized controls are invariant with respect to both external perturbations and parametric errors. The results of mathematical modeling are presented that demonstrate high efficiency of developed control algorithms. Propellant consumption for realizing a programmed turn is numerically estimated taking into account considerable gravitational and aerodynamic moments acting upon the spacecraft.     

Some issues of time-optimal control over a spacecraft programmed turn

The problem of optimal turn of a spacecraft from an arbitrary initial position to a final specified angular position in a minimum time is considered and solved. A case is investigated, when the constraint on spacecraft’s angular momentum during the turn is essential. Based on the quaternion method a solution to the posed problem has been found, and an optimal control program is constructed taking the constraints on controlling moment into account. The optimal control is found in the class of regular motions. A condition (calculation expression) is presented for determining the moment to begin braking with the use of measurements of current motion parameters, which considerably improves the accuracy of putting the spacecraft into a preset position. For a dynamically symmetrical spacecraft the solution to the problem of optimal control by the spacecraft spatial turn is presented in analytical form (expressions in elementary functions). An example of mathematical modeling of the spacecraft motion dynamics under optimal control over reorientation is given.     

基于SQP方法的常推力月球软着陆轨道优化方法

月球软着陆是未来月球探测中的一项关键技术。针对这项技术，本文给出了一种基于SQP方法的常推力月球软着陆轨道优化方法。该方法通过将常推力月球软着陆轨道离散化，利用离散点处状态连续作为约束条件，把常推力月球软着陆轨道优化问题归结为一个非线性规划问题，对于此问题的求解，其初值均为有物理意义的状态和控制量，从而避免了采用传统优化方法在解决此优化问题时对没有物理意义变量初值的猜测。最后，利用SQP方法求解了此轨道优化问题。仿真计算结果表明这种离散化的方法应用于此轨道优化问题可以避免传统轨道优化方法对初值敏感的问题。     

Optimization of Multi-Orbit Transfers between Noncoplanar Elliptic Orbits

Using the maximum principle formalism, the problem of optimizing interorbital transfer between two noncoplanar elliptic orbits is reduced to solution of a boundary value problem for a system of ordinary differential equations. In order to solve the resulting boundary value problem numerically, the numerical homotopic method or modified Newton's method is used. When solving the boundary value problem, the right-hand sides of differential equations of motion are averaged numerically. Efficient software is developed, and a large number of optimal trajectories are calculated using it. As a result of analysis of these numerical data, new high-quality results are obtained. Specifically, a bifurcation of optimal solutions is found, the existence of critical inclination is demonstrated, and a partial classification of the structure of optimal control is performed.     

Optimization of interplanetary trajectories for spacecraft with ideally regulated engines using the continuation method

The problem of optimization of interplanetary trajectories is considered for spacecraft with a small-thrust ideally regulated engine. When the maximum principle is used, determination of the optimal trajectory is reduced to solution of a two-point boundary value problem for a system of ordinary differential equations. In order to solve this boundary value problem, the method of continuation in parameter is used, and with the help of it the formal reduction of the boundary value problem to a Cauchy problem is performed. Different variants of the continuation method are considered, including the method of continuation in the gravitational parameter which allows one to find extreme trajectories with a preset angular distance. The issues of numerical realization of the continuation method are discussed, and numerical examples of its use for solving the problems of optimization of interplanetary trajectories are presented.     

Theory of physical libration of the Moon with the liquid core: Forced librations

For a two-layer model of the Moon that consists of a solid nonspherical mantle and an ellipsoidal homogeneous liquid core, a theory of forced librations under the effect of gravitational Earth’s moments has been developed. The motion of the Moon over its orbit has been described by the high-accuracy theory of DE/LE-4 orbital motion. Tables have been constructed that present forced librations of the Moon caused by the second harmonic of its force function, in the neighborhood of its motion according to the generalized Cassini laws. Disturbances of the first-order with respect to dynamic compressions of the Moon and its core are obtained in analytical form for Andoyer variables and Poincare variables and for the projection of the angular velocity vector of Moon’s mantle rotation and the Poincare coordinate system (relative to which core’s liquid accomplishes simple motion) on its major central axes of inertia, as well as for the classical variables in the Moon libration theory, etc. Constructed tables of the forced librations theory give the amplitudes and periods of librations and combinations of arguments of the orbital motion theory that correspond to libration parameters. The interpretation of basic variations has been given and a comparison with the previous theories has been carried out, in particular with the modern empirical theory constructed based on the laser observation data.     

不确定条件下采用协方差描述函数法的再入轨迹鲁棒优化

Uncertainties are taken into account in reentry trajectory optimization and a robust trajectory optimization model is constructed based on the robust optimization theory. There are computational difficulties in solving this robust trajectory optimization problem due to the random variables associated with the uncertainties. To overcome these difficulties, the covariance analysis describing function technique (CADET) is employed to convert the robust trajectory optimization model into an equivalent deterministic trajectory optimization formulation, which is then solvable by using the existing pseudo spectral method. In the case study, a robust reentry trajectory optimization problem considering the uncertainties of aerodynamic parameters is solved to obtain the robust optimal trajectories. By comparing with the traditional deterministic optimal trajectories, the robust optimal trajectories are significantly less sensitive to the uncertainties of aerodynamic parameters, showing the effectiveness of the presented method.      

伴随卫星轨道保持

文中利用伴随轨道方程推导出伴随卫星相对运动状态转移方程。基于这一状态转移方程，研究了双脉冲校正中的伴随运动状态的变化，并得出需要的脉冲控制量的解析式。利用遗传算法对脉冲控制量进行了优化设计，求得使总脉冲最小的最优变轨时间，脉冲控制量和轨迹。     

The Motion of a Material Point in a Newtonian Field under the Action of Equiangular Thrust

The problems of investigation and optimization of the motion of spacecraft are extensively discussed in the literature. Nevertheless, in many cases a large variety of qualitative characteristics of their motion and of the form of their trajectories are still unclear. In this paper we consider a plane equiangular acceleration of a spacecraft both in a Newtonian field and in its absence (at a large distance from the center of attraction). The general equation of a trajectory of plane acceleration is presented with the introduction of a new variable, an index of an exponent, which allows one to obtain convenient solutions at different values of the time-independent angle of inclination of the vector of thrust to the spacecraft's radius vector (i.e., when equiangular acceleration takes place). Asymptotic solutions are constructed and an interesting fact is revealed. Namely, it is shown that when the center of attraction exists or is absent, for all initial conditions the trajectories appearing at the above equiangular acceleration of a material point tend to the standard logarithmic spirals at a large distance from the center. Specifically, when the value of transverse (perpendicular to the radius vector) thrust is constant, there appears a logarithmic spiral with an angle of inclination to the radius vector equal to 35.264°. Different forms of the trajectory of equiangular acceleration of spacecraft at a low thrust are also studied. The results obtained can be useful for the investigation and choice of optimum space trajectories.     

6-DOF robust adaptive terminal sliding mode control for spacecraft formation flying

This paper addresses the tracking control problem of the leader–follower spacecraft formation, by which we mean that the relative motion between the leader and the follower is required to track a desired time-varying trajectory given in advance. Using dual number, the six-degree-of-freedom motion of the follower spacecraft relative to the leader spacecraft is modeled, where the coupling effect between the translational motion and the rotational one is accounted. A robust adaptive terminal sliding mode control law, including the adaptive algorithms, is proposed to ensure the finite time convergence of the relative motion tracking errors despite the presence of model uncertainties and external disturbances, based on which a modified controller is furthermore developed to solve the dual-equilibrium problem caused by dual quaternion representation. In addition, to alleviate the chattering, hyperbolic tangent function is adopted to substitute for the sign function. And by theoretical analysis, it is proved that the tracking error in such case will converge to a neighborhood of the origin in finite time. Finally, numerical simulations are performed to demonstrate the validity of the proposed approaches.     

间接Legendre伪谱法的欠驱动航天器姿态运动轨迹跟踪

针对仅带有两组喷气推力器的非轴对称欠驱动刚性航天器，提出一种基于间接Legendre伪谱法的姿态运动轨迹跟踪控制算法。首先采用Legendre伪谱法(LPM)离线规划出系统的最短时间姿态机动参考轨迹。接着将实际运行轨迹与参考轨迹之间的偏差作为变量，根据Pontryagin极小值原理必要条件把系统姿态运动跟踪问题转化为一个两点边值问题(TPBVP)。最后采用 Legendre-Gauss-Lobatto(LGL)点将此两点边值问题离散转化为一个线性方程组来求解，避免了对传统Riccati微分方程的积分运算。数值仿真校验了本文基于间接Legendre伪谱法的姿态运动轨迹跟踪控制算法的有效性。     

Analysis of evolution of multi-orbit perturbed trajectories of a spacecraft in a central field of attraction

Properties of differential equations of multi-orbit trajectory motion of a spacecraft are investigated analytically. The spacecraft moves under the action of small perturbations (in particular, low thrust) in the plane of a central Newtonian field of attraction. The conditions are specified for existence of a partial singular aperiodic solution, in the neighborhood of which the behavior of osculating elements changes sharply. In this case, phase variables (the angular position of the pericenter and the true anomaly) are found to undergo the sharpest changes. The exact superposition of solutions is suggested for the equations of motion transformed to the form of a quasi-linear, weakly non-stationary system: a partial singular aperiodic solution and fast solutions oscillating around it. Asymptotic representations are obtained for both components of the superposition. They are fairly exact in the region of smallness of perturbing terms at a long variation of the argument.     

自由返回轨道的P平面参数快速设计方法

针对自由返回轨道求解过程中地心轨道类型变化造成的B平面参数方法计算失败问题,提出一种基于P平面参数的自由返回轨道快速设计方法。首先,基于轨道半通径参数的普适性,给出了不同轨道类型的P平面参数定义,建立了以P平面参数为求解目标量的自由返回轨道求解模型。其次,给出了基于P平面参数的自由返回轨道快速设计方法,在构建的瞬时地月惯性系下,以平面双二体自由返回轨道作为初值,实现了高精度力模型下的自由返回轨道快速求解。对8种构型自由返回轨道的设计结果表明,P平面参数具有类似于B平面参数的大收敛域,且有效解决了轨道类型变化对计算的影响,可直接应用于中国后续月球探测任务轨道设计。     

Optimal control of the deployment process of solar wings on spacecraft

The optimal attitude control problem of spacecraft during the stretching process of solar wings is investigated in this paper. The dynamical equations of the nonholonomic system are derived from the conservation principle of the angular momentum of the multibody system. Attitude control of the spacecraft with internal motion is reduced to a nonholonomic motion planning problem. The spacecraft attitude control is transformed into the steering problem for a drift free control system. The optimal solution for steering a spacecraft with solar wings is presented. The controlled motion of spacecraft is simulated for two cases. The numerical results demonstrate the effectiveness of the optimal control approach.