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.     

Method of continuation for optimization of interplanetary low-thrust trajectories

The problem of local optimization of interplanetary low-thrust trajectories is considered with the use of the maximum principle and continuation numerical methods. Two types of problems are analyzed: problems with limited power and problems with limited thrust. The latter problem is generalized by introducing the dependence of thrust and specific impulse on available electric power. In order to reduce the problem of optimal control to a boundary value problem, the Pontryagin maximum principle is used, and then, using the continuation method, this boundary value problem is reduced to the Cauchy problem. Variants of the continuation method for optimizing low-thrust trajectories are presented in the paper, including a new method of continuation for the limited thrust problem, which does not require any choice of the initial approximation for boundary values of conjugate variables.     

Optimization of transfer trajectories to the Apophis asteroid for spacecraft with high and low thrust

The problem of optimization of a spacecraft transfer to the Apophis asteroid is investigated. The scheme of transfer under analysis includes a geocentric stage of boosting the spacecraft with high thrust, a heliocentric stage of control by a low thrust engine, and a stage of deceleration with injection to an orbit of the asteroid’s satellite. In doing this, the problem of optimal control is solved for cases of ideal and piecewise-constant low thrust, and the optimal magnitude and direction of spacecraft’s hyperbolic velocity “at infinity” during departure from the Earth are determined. The spacecraft trajectories are found based on a specially developed comprehensive method of optimization. This method combines the method of dynamic programming at the first stage of analysis and the Pontryagin maximum principle at the concluding stage, together with the parameter continuation method. The estimates are obtained for the spacecraft’s final mass and for the payload mass that can be delivered to the asteroid using the Soyuz-Fregat carrier launcher.     

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.     

Planar Problem of an Optimal Transfer of a Low-Thrust Spacecraft from High-Elliptic to Geosynchronous Orbit

Low-thrust flights from high-elliptic orbits are of considerable interest, since they allow one to decrease (compared to high-thrust flights) the propulsion consumption and to reduce the flight duration. At the same time, in comparison with the spiral unwinding flights from low near-circular orbits, this scheme minimizes the harmful effect of the radiation belts. Based on the maximum principle, the problem of optimization is reduced to a two-point boundary value problem, which is solved numerically using the modified Newton method. A method is suggested to obtain the initial approximation for solving the boundary value problem. The method takes advantage of the idea of transition from an approximately optimal trajectory to the optimal one. Two problems, which have different low-thrust models, are considered: one with permanently acting low thrust and the other with the possibility of turning it on/off. In both cases no restrictions are imposed on the thrust direction. A comparison of these problems is made. We investigated (i) what gain in the final mass can be attained when passing from the first to the second problem, (ii) at the cost of what loss in flight duration this can be achieved, and (iii) what changes in the optimal program of control must be done in this case.     

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.     

Optimal flights to near-Earth asteroid

The spacecraft flights to the Near-Earth asteroid in order to give an impact influence on the asteroid, correct its orbit and prevent the asteroid’s collision with the Earth are analyzed.In the first part, the impulse flights are analyzed in the Lambert approach. There are determined the optimal trajectories maximizing the asteroid deviation from the Earth.In the second part, the flights with the chemical and electric-jet engines are analyzed. The high thrust is used to launch the spacecraft from the geocentric orbit, and the low thrust is applied for the heliocentric motion. On the base of optimal impulse transfer, the optimal low thrust trajectories are determined using Pontryagin maximum principle.The numerical results are given for the flight to the asteroid Toutatis. Parameters of the spacecraft impact on the asteroid are determined. The asteroid deviation from the Earth caused by the spacecraft influence is presented.     

Solving Optimization Problems for the Flight Trajectories of a Spacecraft with a High-Thrust Jet Engine in Pulse Formulation for an Arbitrary Gravitational Field in a Vacuum

A mathematically well-posed technique is suggested to obtain first-order necessary conditions of local optimality for the problems of optimization to be solved in a pulse formulation for flight trajectories of a spacecraft with a high-thrust jet engine (HTJE) in an arbitrary gravitational field in vacuum. The technique is based on the Lagrange principle of derestriction for conditional extremum problems in a function space. It allows one to formalize an algorithm of change from the problems of optimization to a boundary-value problem for a system of ordinary differential equations in the case of any optimization problem for which the pulse formulation makes sense. In this work, such a change is made for the case of optimizing the flight trajectories of a spacecraft with a HTJE when terminal and intermediate conditions (like equalities, inequalities, and the terminal functional of minimization) are taken in a general form. As an example of the application of the suggested technique, we consider in this work, within the framework of a bounded circular three-point problem in pulse formulation, the problem of constructing the flight trajectories of a spacecraft with a HTJE through one or several libration points (including the case of going through all libration points) of the Earth–Moon system. The spacecraft is launched from a circular orbit of an Earth's artificial satellite and, upon passing through a point (or points) of libration, returns to the initial orbit. The expenditure of mass (characteristic velocity) is minimized at a restricted time of transfer.     

Design feasibility via ascent optimality for next-generation spacecraft

This paper deals with the optimization of the ascent trajectories for single-stage-sub-orbit (SSSO), single-stage-to-orbit (SSTO), and two-stage-to-orbit (TSTO) rocket-powered spacecraft. The maximum payload weight problem is studied for different values of the engine specific impulse and spacecraft structural factor.The main conclusions are that: feasibility of SSSO spacecraft is guaranteed for all the parameter combinations considered; feasibility of SSTO spacecraft depends strongly on the parameter combination chosen; not only feasibility of TSTO spacecraft is guaranteed for all the parameter combinations considered, but the TSTO payload is several times the SSTO payload.Improvements in engine specific impulse and spacecraft structural factor are desirable and crucial for SSTO feasibility; indeed, aerodynamic improvements do not yield significant improvements in payload.For SSSO, SSTO, and TSTO spacecraft, simple engineering approximations are developed connecting the maximum payload weight to the engine specific impulse and spacecraft structural factor. With reference to the specific impulse/structural factor domain, these engineering approximations lead to the construction of zero-payload lines separating the feasibility region (positive payload) from the unfeasibility region (negative payload).     

Online generation of quasi-optimal spacecraft rendezvous trajectories

A direct method for rapid generation of combined time-propellant near-optimal trajectories of proximity maneuvers of a chaser spacecraft required to dock a target one, with predetermined thrust history along a master direction, is presented. The predetermined thrust history is generated by applying the Pontryagin maximum principle. The new direct method, already implemented and tested on board real aircraft, is based on three concepts: high-order polynomials as reference functions, preset on–off sequence of a master control, and reduction of the optimization problem to the determination of a small set of parameters. Presetting the master control, the remaining controls act as slaves, guarantying the chaser to move along the desired path. Seeking of the optimum strategy is transformed into a nonlinear programming problem, and then numerically solved through an ad hoc algorithm in accelerated time scale. Examples are reported to prove the rapidness of the approach to generate a sub-optimal docking trajectory.     

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.     

UKF参数估计在航天器气动辅助变轨问题中的应用

为快速精确地求解气动辅助变轨问题,提出一种基于无损卡尔曼滤波(UKF)参数估计的数值求解方法。首先,针对气动辅助变轨问题,利用极大值原理将其转化为对应的两点边值问题;然后,以协态变量的初值作为估计参数,以末端条件为期望观测值,将该两点边值问题转化为参数估计问题,并应用UKF滤波算法求解。该算法基于估计理论,避免了计算传统数值方法所需要的梯度矩阵,同时克服了猜测协态变量初值的困难,降低了求解气动辅助变轨问题的难度。数值仿真表明,该算法结构简单,求解效率高,具有良好的鲁棒性。     

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.     

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.     

基于凸优化理论的含约束月球定点着陆轨道优化

针对月球精确定点软着陆问题,考虑导航及障碍检测敏感器视场约束及制动发动机推力大小约束,对月球动力下降段轨道优化方法进行了研究。首先建立了含约束条件的三维定点软着陆轨道优化问题模型,根据庞德亚金极小值原理推导了最优推力开关方程,并给出了推力奇异区间不存在的证明。针对优化模型中的复杂非线性约束,引入凸优化理论将问题转化为二阶锥优化问题,并采用内点法求解了最优标称轨迹。最后给出了月球软着陆制动段、接近段的仿真结果,验证了该着陆轨道优化方法的有效性。     

Optimal impulse point-to-elliptical orbit transfers of a spacecraft

The paper deals with energetically optimal multi-impulse transfer of a spacecraft in the central Newtonian gravity field near a planet. At the initial state of the transfer the distance from the spacecraft to the center of attraction, its radial and transversal velocity projections are known. At the end of the transfer the spacecraft must be located in the elliptical orbit with the given area and energy constants. The distance from the spacecraft to the center of attraction is bounded above and below, the transfer time being unspecified. The initial orbit intersects the inner boundary of the given ring.All the optimal solutions have been obtained by analytical way. A number of new solutions has been found for the given problem in comparison with the case of the transfer from the orbit at the free initial point.Up to five impulses can be applied on the optimal trajectories. The numerical simulation of the problem is carried out. It shows that all obtained solutions give not only local but global optimal energetic input on the corresponding conditions.     

应用非线性规划求解异面最优轨道转移问题

研究了一种应用非线性规化求解有限推力作用下异面最优轨道转移问题的方法。采用改进春分点根素形式的高斯行星方程，从庞德里亚金极小值原理出发，将有限推力作用下异面最优轨道转移问题转化为两点边值问题；在考虑边界条件、横截条件及开关函数的前提下，将两点边值问题转化为针对协状态初值等的参数优化问题；最后应用非线性规划方法求解形成的参数优化问题。本方法特点是能得到开关函数从而得到最优发动机开关机逻辑。文章最后通过一个仿真计算，演示了完整的求解过程，验证了方法的有效性。     

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

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

月球最优软着陆两点边值问题的数值解法

对于月球软着陆,燃耗最优是制导过程的基本要求.文中首先应用极大值原理设计了最优着陆制导控制律,此时求解最优轨迹变成一个两点边值问题(TPBVP).本文利用一种基于初值猜测技术的打靶法求解这个两点边值问题,得到软着陆最优轨迹.结果表明该方法可有效改善迭代计算,具有一定的优越性.