全太阳系深空探测轨迹全局优化(英文)

针对太阳系中全部的248997颗行星的探测问题,给出了一种关于探测飞行器的深空探测全局四维轨迹(t,x,y,z)优化方案,即飞行器从地球发射进入太阳系并采用小推力控制,优化方案的性能指标为飞行器与太阳系中全部行星中相遇和交会的星的数量最多并且燃料消耗最少。本方案给出了四维飞行轨迹进行全局优化的一套算法,该算法由搜索算法和四维轨迹优化算法组成。此搜索算法从太阳系的248997颗行星中寻找获得尽可能多的经过近地球3维走廊内的行星;而四维轨迹优化算法由改进的动态规划算法、基于最优控制理论的共轭梯度算法和静态参数优化算法组成,其中静态参数优化算法用于搜索最优发射时间窗口。基于该组合算法,通过长时间的大规模的飞行数字仿真,最终计算出探测器的四维最优飞行轨迹,在一年内路过了太阳系中全部行星中的12颗行星。     

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

The problem of a rendezvous in the central Newtonian gravitational field is considered for a controlled spacecraft and an uncontrollable spacecraft moving along an elliptic Keplerian orbit. For solving the problem, two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems and using quaternion variables to describe the orientations of these coordinate systems. In the first variant of the equations of motion a quaternion variable characterizes the orientation of an instantaneous orbit of the spacecraft and the spacecraft location in the orbit, while in the second variant it characterizes the orientation of the plane of the spacecraft instantaneous orbit and the location of a generalized pericenter in the orbit. The quaternion variable used in the second variant of the equations of motion is a quaternion osculating element of the spacecraft orbit. The problem of a rendezvous of two spacecraft is formulated as a problem of optimal control by the motion of the center of mass of a controlled spacecraft with a movable right end of the trajectory, and it is solved on the basis of Pontryagin's maximum principle.     

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

The results of numerical solution of the problem of a rendezvous in the central Newtonian gravitational field of a controlled spacecraft with an uncontrollable spacecraft moving along an elliptic Keplerian orbit are presented. Two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems and using quaternion variables to describe the orientations of these coordinate systems. The problem of a rendezvous of two spacecraft is formulated [1, 2] as a problem of optimal control by the motion of the center of mass of a controlled spacecraft with a movable right end of the trajectory, and it is solved on the basis of Pontryagin's maximum principle. The paper is a continuation of papers [1, 2], where the problem of a rendezvous of two spacecraft has been considered theoretically using the two above variants of the equations of motion for the center of mass of the controlled spacecraft.     

Optimal Control of a Programmed Turn of a Spacecraft

A problem of optimal turn of a spacecraft is considered. The time of turn is minimized, as well as the functional having a meaning of the propellant consumption. An analytical solution to the problem stated is derived. It is demonstrated that the solution optimal in this sense belongs to a class of two-impulse controls, under which a spacecraft executes the turn along the trajectory of its free motion. The solution obtained in this paper differs from earlier available solutions considerably. The estimations of the propellant consumption for a realization of the programmed turn are made.     

Quaternion regularization and trajectory motion control in celestial mechanics and astrodynamics: II

Problems of regularization in celestial mechanics and astrodynamics are considered, and basic regular quaternion models for celestial mechanics and astrodynamics are presented. It is shown that the effectiveness of analytical studies and numerical solutions to boundary value problems of controlling the trajectory motion of spacecraft can be improved by using quaternion models of astrodynamics. In this second part of the paper, specific singularity-type features (division by zero) are considered. They result from using classical equations in angular variables (particularly in Euler variables) in celestial mechanics and astrodynamics and can be eliminated by using Euler (Rodrigues-Hamilton) parameters and Hamilton quaternions. Basic regular (in the above sense) quaternion models of celestial mechanics and astrodynamics are considered; these include equations of trajectory motion written in nonholonomic, orbital, and ideal moving trihedrals whose rotational motions are described by Euler parameters and quaternions of turn; and quaternion equations of instantaneous orbit orientation of a celestial body (spacecraft). New quaternion regular equations are derived for the perturbed three-dimensional two-body problem (spacecraft trajectory motion). These equations are constructed using ideal rectangular Hansen coordinates and quaternion variables, and they have additional advantages over those known for regular Kustaanheimo-Stiefel equations.     

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.     

One optimization problem for trajectories of spacecraft rendezvous mission to a group of asteroids

The optimization problem is considered for the trajectory of a spacecraft mission to a group of asteroids. The ratio of the final spacecraft mass to the flight time is maximized. The spacecraft is controlled by changing the value and direction of the jet engine thrust (small thrust). The motion of the Earth, asteroids, and the spacecraft proceeds in the central Newtonian gravitational field of the Sun. The Earth and asteroids are considered as point objects moving in preset elliptical orbits. The spacecraft departure from the Earth is considered in the context of the method of a point-like sphere of action, and the excess of hyperbolic velocity is limited. It is required sequentially to have a rendezvous with asteroids from four various groups, one from each group; it is necessary to be on the first three asteroids for no less than 90 days. The trajectory is finished by arrival at the last asteroid. Constraints on the time of departure from the Earth, flight duration, and final mass are taken into account in this problem.     

一种航天器空间机动轨道的改进形状设计方法

针对带推力约束的航天器三维空间机动轨道初始设计问题,提出了一种基于傅立叶级数展开的改进形状设计方法。首先,在柱坐标系下建立了航天器运动模型,并将基于傅立叶级数展开的形状方法推广到三维空间的机动轨道初始设计;然后,基于所得到的空间初始机动轨道,采用直接配点法进行了完整三维空间机动轨道的优化设计。仿真结果表明,提出的方法可以为带有任务约束的航天器三维空间机动轨道的优化设计提供更优的初始参数及其解析解,为空间机动轨道设计提供了新的方法。     

Modes of uncontrolled rotational motion of the <Emphasis Type="Italic">Progress M-29M</Emphasis> spacecraft

We have reconstructed the uncontrolled rotational motion of the Progress M-29M transport cargo spacecraft in the single-axis solar orientation mode (the so-called sunward spin) and in the mode of the gravitational orientation of a rotating satellite. The modes were implemented on April 3–7, 2016 as a part of preparation for experiments with the DAKON convection sensor onboard the Progress spacecraft. The reconstruction was performed by integral statistical techniques using the measurements of the spacecraft’s angular velocity and electric current from its solar arrays. The measurement data obtained in a certain time interval have been jointly processed using the least-squares method by integrating the equations of the spacecraft’s motion relative to the center of mass. As a result of processing, the initial conditions of motion and parameters of the mathematical model have been estimated. The motion in the sunward spin mode is the rotation of the spacecraft with an angular velocity of 2.2 deg/s about the normal to the plane of solar arrays; the normal is oriented toward the Sun or forms a small angle with this direction. The duration of the mode is several orbit passes. The reconstruction has been performed over time intervals of up to 1 h. As a result, the actual rotational motion of the spacecraft relative to the Earth–Sun direction was obtained. In the gravitational orientation mode, the spacecraft was rotated about its longitudinal axis with an angular velocity of 0.1–0.2 deg/s; the longitudinal axis executed small oscillated relative to the local vertical. The reconstruction of motion relative to the orbital coordinate system was performed in time intervals of up to 7 h using only the angularvelocity measurements. The measurements of the electric current from solar arrays were used for verification.     

欠驱动刚性航天器姿态的非完整运动规划粒子群算法

研究欠驱动刚性航天器姿态的非完整运动规划问题。众所周知航天器利用三个动量飞轮可以控制其姿态和任意定位，当其中一轮失效，航天器动力学方程表现为不可控。在系统角动量为零的情况下，系统的姿态控制问题可转化为无漂移系统的运动规划问题。基于粒子群优化技术设计了欠驱动刚性航天器姿态的非完整运动规划算法。通过数值仿真，并和遗传算法进行了比较，结果表明该方法对欠驱动航天器姿态运动规划是有效的。     

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.     

一种基于EPSO的航天器交会轨迹优化方法

研究了航天器在固定时间内燃料最省的多脉冲交会问题,提出了一种基于种群熵粒子群优化 (Population Entropy based Particle Swarm Optimization,EPSO)算法的交会轨迹优化设计方法。采用线性化C\|W方程描述航天器的相对运动,以能耗最优为控制目标,得到了基于连续推力的最优转移轨迹,用于确定脉冲点的位置。考虑工程实用性,采用多脉冲控制方法,利用脉冲点的位置参数建立了以脉冲点时间间隔为决策变量的优化目标函数,并用EPSO算法进行求解。在EPSO中,种群熵描述粒子在搜索空间中位置分布的混乱程度,并通过上一代的种群熵确定下一代的搜索空间,从而减少搜索空间的浪费,提高了算法的搜索速度和收敛精度。仿真结果表明,算法本身具有良好的优化性能,适用于航天器轨迹优化。     

Multi-structural estimation of an integrated navigation system of a reusable spacecraft

An approach to the synthesis of an integrated navigation system is considered for a reusable space-craft that performs an arbitrary spatial maneuver under the conditions of internal and external disturbances. The offered approach provides for a noise-suppressing solution of the navigation problem, both in a regular mode of spacecraft motion, and during its descent along the unplanned trajectory.     

Quaternion regularization in celestial mechanics and astrodynamics and trajectory motion control. I

Regularization problems in celestial mechanics and astrodynamics are considered. The fundamental regular quaternion models of celestial mechanics and astrodynamics are presented. It is shown that the efficiency of analytical investigation and numerical solution of boundary problems of optimal trajectory motion control of spacecraft may be increased using quaternion astrodynamics models. The regularization problem of celestial mechanics and astrodynamics that implies eliminating the feature, which arises in the equations of the two-body problem in case of impact of the second body with the central body, is considered in the first section of the paper. The quaternion method for regularizing the equations of the perturbed spatial two-body problem suggested by the author is presented; the method is compared with Kustaanheimo-Stiefel (KS) regularization. Demonstrative geometric and kinematic interpretations of regularizing transformations are provided. Regular quaternion equations for the two-body problem, which generalize the regular Kustaanheimo-Stiefel equations, as well as regular equations in quaternion osculating elements and quaternion regular equations for perturbed central motion of a material point, are considered. The papers on quaternion regularization in celestial mechanics and astrodynamics are briefly analyzed.     

万有引力场中带挠性太阳帆板航天器的姿态稳定性

本文讨论带双侧挠性太阳帆板航天器在万有引力场中的姿态运动，建立带挠性帆板航天器的欧拉方程和帆板强迫振动方程。利用Ｇａｌｅｒｋｉｎ方法对动力学方程离散化，利用Ｋｅｌｖｉｎ－Ｔａｉｔ－Ｃｈｅｔａｙｅｖ定量判断航天器在轨道坐标系内相对平衡的稳定性。导出适用于任意阶模态的解析形式稳定性充分条件。     

月球探测器直接软着陆最优轨道设计

研究月球探测器直接软着陆最优轨道的设计问题。首先根据探测器直接软着陆的特点，提出了有限推力最省燃料的最优轨道设计问题；然后利用有限推力月面软着陆的最优推力控制方向的计算公式，研究了边值条件和计算方法；最后通过直接软着陆最优轨道的算例及结果分析，发现开始制动高度越低越省能量；推力方向可变时比不可变时节省能量；推力大小可变相当于采用了多级制动，对安全定点着陆非常有利。     

Ballistics, navigation, and motion control for a spacecraft during its landing on the surface of Phobos

Basic concepts and algorithms laid as foundations of the scheme of landing on the Martian moon Phobos (developed for the Phobos-Grunt project) are presented. The conditions ensuring the landing are discussed. Algorithms of onboard navigation and control are described. The equations of spacecraft motion with respect to Phobos are considered, as well as their use for correction of the spacecraft motion. The algorithm of estimation of the spacecraft’s state vector using measurements with a laser altimeter and Doppler meter of velocity and distance is presented. A system for modeling the landing with a firmware complex including a prototype of the onboard computer is described.     

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.     

间接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.