Optimization of transfers under constraints on the thrust direction: II

This paper completes the study of optimal transfers with constraints imposed on the thrust vector direction that was opened by paper [1]. The linear inhomogeneous and homogeneous constraints on the thrust direction are considered (specified either by equalities or inequalities), as well as mixed constraints. Some examples of the constraints are presented. A modified method of the transporting trajectory is applied in order to find the optimal transfer under the linear constraints on the thrust direction. This method also gives the necessary condition for a transfer possibility at a given constraint on the thrust direction. A numerical example is considered, in which the propellant consumption is analyzed for the cases of transfers with and without constraints.     

A Modification of the Method of Transporting Trajectory

Further development of an approximate method for optimizing a flight with an ideally controlled small thrust is proposed. The method is based on the employment of the transporting trajectory and considered in [1–3]. A detailed analysis of the means of improving the accuracy of this method suggested in [2, 3] is carried out, and the solution is presented in finite form. The proposed approach is applied to the flights making flybys of many celestial bodies. In the case of small bodies the solution is also obtained in finite form. A numerical example is considered confirming the high efficiency of this method.     

小推力最优轨迹协态估计的高效机器学习方法

针对变比冲小推力轨迹间接优化中的协态变量初值猜测问题,提出了一种基于机器学习的协态变量初值高精度高效估计方法。首先,基于标称最优轨迹延拓,建立了状态量边值高扰动上限情形下的数据集生成方法,并分析了扰动上限对求解效率的影响。然后,构建了基于位置速度、轨道根数和改进春分点轨道根数多形式状态量组合输入的人工神经网络(ANN)映射关系,分析并优化了神经网络结构。将提出的方法应用于深空探测小推力转移场景,仿真结果表明该方法相对于标称轨迹直接扰动的数据集生成方法及单一形式状态量输入的人工神经网络映射方法,均有效地提升了求解收敛率,能够高效高精度地估计协态变量初值,实现轨迹快速优化。     

Fuel optimum low-thrust elliptic transfer using numerical averaging

Low-thrust electric propulsion is increasingly being used for spacecraft missions primarily due to its high propellant efficiency. As a result, a simple and fast method for low-thrust trajectory optimization is of great value for preliminary mission planning. However, few low-thrust trajectory tools are appropriate for preliminary mission design studies. The method presented in this paper provides quick and accurate solutions for a wide range of transfers by using numerical orbital averaging to improve solution convergence and include orbital perturbations. Thus, preliminary trajectories can be obtained for transfers which involve many revolutions about the primary body. This method considers minimum fuel transfers using first-order averaging to obtain the fuel optimum rates of change of the equinoctial orbital elements in terms of each other and the Lagrange multipliers. Constraints on thrust and power, as well as minimum periapsis, are implemented and the equations are averaged numerically using a Gausian quadrature. The use of numerical averaging allows for more complex orbital perturbations to be added in the future without great difficulty. The effects of zonal gravity harmonics, solar radiation pressure, and thrust limitations due to shadowing are included in this study. The solution to a transfer which minimizes the square of the thrust magnitude is used as a preliminary guess for the minimum fuel problem, thus allowing for faster convergence to a wider range of problems. Results from this model are shown to provide a reduction in propellant mass required over previous minimum fuel solutions.     

小行星探测最优两脉冲交会轨道设计与分析

小行星探测已经成为新世纪深空探测的一个新热点和未来世界航天发展的一个新方向。转移轨道的设计和探测目标可接近性的分析是小行星探测的关键技术之一。现利用了任意两个非共面非共轴椭圆轨道之间的最优两脉冲转移方法，对我国提出的探测Ivar小行星的交会转移轨道进行了设计与分析，给出了全局最优两脉冲交会轨道的设计参数，并利用此方法对近地小行星的可接近性进行了分析和排序，给出了可接近性较好的40颗近地小行星的转移轨道设计参数。这些研究结果对于近地小行星探测任务的目标选择和发射机会的预测都有重要的参考价值。     

Optimal ballistically captured Earth–Moon transfers

The optimality of a low-energy Earth–Moon transfer terminating in ballistic capture is examined for the first time using primer vector theory. An optimal control problem is formed with the following free variables: the location, time, and magnitude of the transfer insertion burn, and the transfer time. A constraint is placed on the initial state of the spacecraft to bind it to a given initial orbit around a first body, and on the final state of the spacecraft to limit its Keplerian energy with respect to a second body. Optimal transfers in the system are shown to meet certain conditions placed on the primer vector and its time derivative. A two point boundary value problem containing these necessary conditions is created for use in targeting optimal transfers. The two point boundary value problem is then applied to the ballistic lunar capture problem, and an optimal trajectory is shown. Additionally, the problem is then modified to fix the time of transfer, allowing for optimal multi-impulse transfers. The tradeoff between transfer time and fuel cost is shown for Earth–Moon ballistic lunar capture transfers.     

行星际Halo轨道转移的多目标优化设计

研究了星历约束下不同太阳—行星系统Halo轨道间转移的多目标优化设计问题。分析了直接和间接两种转移方式的特点,并引入伪流形技术加快了Halo轨道的逃逸和捕获速度。构建了两种转移机制的多目标优化模型。对于直接转移方式,采用伪流形双向拼接策略实现了转移轨道的构建;对于间接转移方式,通过近拱点庞加莱映射与双曲超速匹配完成了转移轨道的拼接。进一步,采用多项式样条函数对伪流形进行逼近,提高了伪流形的计算效率。两种转移机制的轨道优化设计都可以归结为简单的多变量无约束优化问题,采用非支配快速排序遗传算法NSGA-II求解。对地球—火星Halo轨道间的转移进行了多目标优化设计,校验了本文方法的有效性。     

Inter-orbital low-thrust transfers in an arbitrary field of forces

Low-thrust transfers between preset orbits are considered in the presence of perturbations of different origin. A simple method of finding the transfer trajectory is suggested, based on linearization of motion near reference orbits. The required accuracy of calculations is achieved by way of increasing the number of reference orbits. The method can also be used in the case of a large number of revolutions around the attracting center: no averaging of motion is required in this case. The suggested method is applicable as well, when the final orbit is specified partially and when there are constraints on the thrust direction. The optimal solution to the linearized problem is not optimal for the original problem; closeness of solutions to these two problems is estimated using a numerical example. Capabilities of the method are also illustrated by examples.     

Low-energy Earth–Moon transfers involving manifolds through isomorphic mapping

Analysis and design of low-energy transfers to the Moon has been a subject of great interest for decades. Exterior and interior transfers, based on the transit through the regions where the collinear libration points 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 concerned with a geometrical approach for low-energy Earth-to-Moon mission analysis, based on isomorphic mapping. The isomorphic mapping of trajectories allows a visual, 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. Two types of Earth-to-Moon missions are considered. The first mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a periodic orbit around the Moon. The second mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a capture (non-periodic) orbit around the Moon. In both cases three velocity impulses are needed to perform the transfer: the first at an unknown initial point along the low Earth orbit, the second at injection on the stable manifold, the third at injection in the final (periodic or capture) orbit. The final goal is in finding the optimization parameters, which are represented by the locations, directions, and magnitudes of the velocity impulses such that the overall delta-v of the transfer is minimized. This work proves how isomorphic mapping (in two distinct forms) can be profitably employed to optimize such transfers, by determining in a geometrical fashion the desired optimization parameters that minimize the delta-v budget required to perform the transfer.     

Optimization of passive orbit with the use of gravity maneuver

A high-precision method of calculating gravitational interactions is applied in order to determine optimal trajectories. A number of problems, necessary for determination of optimal parameters at a launch of a spacecraft and during its flyby near celestial bodies, are considered. The spacecraft trajectory was determined by numerical integration of the equations of passive motion of the spacecraft and of the equations of motion for planets, the Sun, and the Moon. The optimal trajectory of the spacecraft approaching the Sun is determined by fitting its initial conditions.     

载人月面着陆与上升飞行器辐射环境适应性研究

载人月面着陆与上升飞行器是在载人月球探测任务中用于月球轨道与月面之间人员、货物往返的运输工具。针对载人月面着陆与上升飞行器全寿命周期内的辐射环境进行了分析,对近地轨道、地月转移及环月轨道的辐射环境中的高能带电粒子能谱进行了对比,对其辐射环境适应性进行了评估,提出了针对性的辐射防护建议,供月面着陆与上升飞行器设计参考。     

A low-thrust transfer between the Earth–Moon and Sun–Earth systems based on invariant manifolds

A low-energy, low-thrust transfer between two halo orbits associated with two coupled three-body systems is studied in this paper. The transfer is composed of a ballistic departure, a ballistic insertion and a powered phase using low-thrust propulsion to connect these two trajectories. The ballistic departure and insertion are computed by constructing the unstable and stable invariant manifolds of the corresponding halo orbits, and a complete low-energy transfer based on the patched invariant manifolds is optimized using the particle swarm optimization (PSO) algorithm on the criterion of smallest velocity discontinuity and limited position discontinuity (less than 1 km). Then, the result is expropriated as the boundary conditions for the subsequent low-thrust trajectory design. The fuel-optimal problem is formulated using the calculus of variations and Pontryagin's Maximum Principle in a complete four-body dynamical environment. Then, a typical bang–bang control is derived and solved using the indirect method combined with a homotopic technique. The contributions of the present work mainly consist of two points. Firstly, the global search method proposed in this paper is simply handled using the PSO algorithm, a number of feasible solutions in a fairly wide range can be delivered without a priori or perfect knowledge of the transfers. Secondly, the indirect optimization method is used in the low-thrust trajectory design and the derivations of the first-order necessary conditions are simplified with a modified controlled, restricted four-body model.     

Optimization of low-thrust orbit transfers with constraints

Approximate numerical methods of optimization are presented for multi-orbit noncoplanar orbit transfers of low-thrust spacecraft. The linear representation of derivatives of boundary parameters is used in the vicinity of a reference trajectory with its discretization into small segments. For each segment a set of pseudo-impulses is introduced, representing possible directions of the thrust vector. In order to solve essentially nonlinear problems, the iterative process of upgrading the reference trajectory is used. At each iteration the linear programming problem of high dimensionality is solved, its boundary conditions being represented in the form of a linear matrix equation. Interval constraints are considered in the form of blocking the propulsion system operation on shadow segments of the orbit, as well as cycling constraints, and constraints on total duration of maneuvers at certain trajectory segments. The results of comparison with solutions obtained by other methods are presented together with examples illustrating the convergence of iterative processes. Optimizations of the trajectories for launching geosynchronous satellites to their orbits and of the trajectories of a noncoplanar transfer from low to high-elliptic Molniya orbit are considered under these constraints.     

考虑机动系数的标准轨迹再入制导方法

针对再入过程中标准轨迹的实际长度需要通过多次更新或迭代求解的问题，提出了一种基于机动系数的通用再入轨迹设计方法。该方法将机动系数定义为真实轨迹长度与初始纵程之比，将再入过程中的机动性进行了量化，可以快速获得到达指定目标可行的轨迹长度；采用真实轨迹长度作为设计参考阻力加速度剖面的依据，避免了轨迹长度的迭代，简化了再入轨迹的生成流程；轨迹曲率问题采用动态航向偏差走廊的方法，控制终端航向偏差、剩余航程满足设计需求。设计轨迹跟踪控制器进行参考轨迹跟踪，完成再入制导。在机动系数区间内指定机动系数进行了数值仿真。仿真结果表明，所提出的标准轨迹制导方法能够快速生成满足路径及终端约束的标准轨迹，且轨迹跟踪效果良好，有较好的应用潜力。     

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

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

制导工具误差折合的遗传主成分方法

制导工具误差折合是导弹精度评定中的重要问题,传统主成分方法中主要成分选择根据试验弹道环境函数矩阵特征值选择,仅考虑了试验弹道的影响,而在试验弹道中作为主要成分,在全程弹道中并不一定是主要成分。提出了充分考虑试验弹道与全程弹道特点的主成分确定方法,并利用遗传算法确定最佳主成分子集。仿真计算表明,与最小二乘和经典主成分方法相比,遗传主成分算法获得的全程弹道遥外差与落点偏差更接近于真实值。     

吸气式超声速导弹爬升段多约束轨迹优化

针对吸气式超声速导弹飞行过程多约束及强耦合的特性,研究了超声速导弹爬升段的轨迹优化设计问题。考虑吸气式推进系统与气动力、飞行轨迹的耦合,对超声速导弹冲压发动机的性能进行分析,揭示了吸气式发动机推力、静压裕度以及余气系数随飞行状态的变化规律;在考虑过载、动压、终端弹道参数及发动机参数等约束的条件下,建立多约束条件下的轨迹优化模型,提出一种适用于此类飞行器飞行轨迹与推力规律的优化设计方法,并对最小油耗的爬升弹道进行优化设计分析。仿真结果表明,该方法能有效解决吸气式超声速导弹多约束轨迹优化问题,可为吸气式超声速导弹的弹道规划与制导律设计提供参考。     

Definition and analysis of properties of optimal three-impulse point-to-orbit transfers under time constraint

This paper deals with energetically optimal multi-impulse transfers of a spacecraft in the central Newtonian gravitational field near a planet. The transfer from a point on initial orbit to the final orbit with the given angular momentum and energy constants is considered. The transfer time is bounded above.With the distance from spacecraft to planet limited and the time free, such parameters of given orbits are chosen that the 3-impulse apsidal transfer T_r is optimal with an intermediate impulse at the maximum distance. On the basis of necessary optimality conditions an algorithm is developed to numerically determine the desired optimal transfer trajectory T_t under time constraint, the apsidal trajectory T_r being taken as initial approach. From the geometry and energy viewpoints, both trajectories T_t and T_r are close to each other. The trajectory T_t is also 3-impulsive, all impulses on it are nonapsidal. The distance from the planet is larger and the sum of impulses is less for this trajectory than for the initial trajectory T_r with the same transfer time.The simplified solution of the problem is constructed producing good approximation to the exact numerical optimization results. The solution asymptotics is found when the transfer time tends to infinity.     

助推-滑翔导弹射程管理技术研究

针对助推-滑翔导弹全程弹道设计问题,提出了基于能量管理的射程管理技术,采用助推段能量管理机动和滑翔段阻力加速度能量管理方法,研究并提出了具体的射程管理方案,分析了不同射程管理方案对射程的影响及其射程区间,验证了通过能量管理实现射程管理的可行性,确定了助推-滑翔导弹射程的可覆盖范围,并给出了一组特定射程下的飞行参数。研究结果表明,通过能量管理技术可实现大范围的射程调节,最小射程可到最大射程的49.3%,采用该射程管理技术可实现助推-滑翔导弹弹道快速、灵活设计,为其发射参数的装订提供了一种新的途径。