线性伪谱模型预测能量最优姿态机动控制方法

针对大气层外飞行器大角度姿态机动控制问题,提出了一种能量最优的线性伪谱模型预测大角度姿态机动控制方法。首先,通过离线弹道规划获得满足初始、终端约束且能量最优的姿态机动控制轨迹;然后,以离线弹道为基准对姿态动力学方程进行小扰动线性化处理,获得以状态偏差为自变量的线性误差传播方程;最后,以能量最优作为性能指标,通过高斯伪谱法对原问题进行离散,推导获得满足终端偏差修正的控制解析表达式。数值计算和蒙特卡罗仿真表明,该方法不仅计算精度高、求解速度快,满足实时计算要求,而且具有较强的鲁棒性,能够实时消除各种干扰。此外,在同等控制精度条件下,该方法相对传统线性二次型调节器（LQR）跟踪方法,能量消耗减小10%。      

Indirect optimization for finite thrust orbit transfer and cooperative rendezvous using an initial guess generator

In this paper, the finite thrust nonplanar transfer and cooperative rendezvous are investigated. The indirect method is adopted in the optimization process. The energy-optimal to fuel-optimal continuation is performed by using the homotopy method. An initial guess generator (IGG) is designed based on a simplified dynamical model to overcome the difficulty in providing the initial guess. Some reasonable simplifications on the dynamical model are applied based on the near-circular near-coplanar transfer. The initial guess can be obtained by analytically solving a system of linear equations. Therefore, the IGG has high computational efficiency. Moreover, the terminal time and state are estimated by IGG to deal with the cooperative rendezvous problem. Numerical examples are presented and the validity of the “IGG-Homotopy” combined algorithm is verified.     

Spacecraft fuel-optimal and balancing maneuvers for a class of formation reconfiguration problems

This paper presents fuel optimal and balancing methodologies for reconfiguring multiple spacecraft in formation subject to a Newtonian gravity field. For a kind of continuous-thrust propulsion system, a fuel-optimal control problem is formulated to minimize the integral squared control subject to the linearized Hill or Clohessy–Wiltshire dynamics of relative motion with respect to a circular reference orbit. Palmer’s analytical solution for general reconfiguration is adapted to maneuvers between projected circular orbits, resulting in the optimal fuel consumption index as a function of configuration parameters such as orbit radius, phase angle, and transfer time. Parametric analyses reveal unique characteristics of individual fuel optimality and gross fuel consumption: for an arbitrary selection of initial/terminal orbit radii, (i) there exist special transfer times such that individual fuel consumption is optimally uniform for all phase angles, and (ii) the total fuel expenditure for a group of three or more spacecraft is invariant for the relatively same configuration with respect to the departure phase. These results serve to effectively design fuel balancing strategies for formation reconfiguration of multiple spacecraft.     

Solar sails for perturbation relief: Application to asteroids

A major cause of spacecraft orbital variation comes from natural perturbations, which, in close proximity of a body, are dominated by its non-spherical nature. For small bodies, such as asteroids, these effects can be considerable, given their uneven (and uncertain) mass distribution. Solar sail technology is proposed to reduce or eliminate the net secular effects of the irregular gravity field on the orbit. Initially, a sensitivity analysis will be carried out on the system which will show high sensitivity to changes in initial conditions. This presents a challenge for optimisation methods which require an initial guess of the solution. As such, the Genetic Algorithm (GA) is proposed as the preferred optimisation method as this requires no initial guess from the user. A multi-objective optimisation is performed which aims to achieve a periodic orbit whilst also minimising the effort required by the sail to do so. Given the system sensitivity, the control law for one orbit is not necessarily applicable for any subsequent orbit. Therefore, a new method of updating the control law for subsequent orbits is presented, based on linearisation and use of a Control Transition Matrix (CTM). The techniques will later find application in a multiple asteroid rendezvous mission with a solar sail as the primary propulsion system.     

冗余飞轮姿控系统控制分配与重构研究

针对具有多冗余执行机构的航天器姿态控制系统,给出了处理冗余的控制分配方案,该方案考虑了幅值和速度约束条件并实现了二次最优．仿真验证了分配环节的有效性,并比较了不同求解算法下的动态分配效果．针对执行机构可能出现的故障,提出了结合故障诊断与隔离系统的执行环节控制重构方法,实现了机构故障下的容错控制,通过卫星姿态跟踪控制仿真,验证了分配和重构环节与控制回路的相容性．     

Hybrid optimization for multiple-impulse reconfiguration trajectories of satellite formation flying

A hybrid optimization method is developed for fuel-optimal reconfigurations of a group of satellites flying in formation. The genetic algorithm performs a global search to find two-impulse trajectories, and primer vector analysis finds multiple-impulsive local optimal trajectories with the two-impulse trajectories as initial guesses. Hybrid optimization finds globally optimal trajectories for formation reconfigurations, including formation resizing, reassignment and reorientation maneuvers. Multiple-impulse trajectories reduce the fuel consumption from the two-impulse trajectories by up to 4.4% for those maneuvers. In real missions, satellites can follow two-impulse trajectories to gain the advantage of a smaller number of impulses, with the cost of slightly more propellant. The qualitative characteristics of the optimal trajectories are analyzed from the number of optimal trajectories found by hybrid optimization.     

全电推进卫星轨道优化的推力同伦解法

全电推进卫星的入轨过程是一个典型的多圈小推力轨道优化问题，由于其推力器加速度小，变轨圈数多，造成其最优理论解的求解较困难。为解决该问题，利用最优控制理论建立了全电推进卫星变轨优化的间接法模型，将变轨优化问题转化为协态变量初值猜测的两点边值问题。从大推力问题开始，通过遗传算法获得大范围猜测值并结合系列二次规划方法获得大推力的精确解。采用推力同伦思想，使用逐渐缩小推力的方式完成小推力问题的求解。仿真算例表明，采用推力同伦的方法，通过数十次的推力缩减即可有效解决多达上百圈变轨的静止轨道全电推进卫星入轨优化问题。     

Energy-optimal reconfiguration of satellite formation flying in the presence of uncertainties

In this study, a two-step control methodology is developed for energy-optimal reconfiguration of satellites in formation in the presence of uncertainties or external disturbances. First, based on a linear deterministic system model, an optimal control law is analytically determined such that a satellite maneuvers from an initial state to a final state relative to another satellite. The structure of this optimal solution is predetermined and simply given by a linear combination of the fundamental matrix solutions associated with the original equations of relative motion. Only the coefficients are to be determined to satisfy given initial and final conditions. In the second step, an uncertain nonlinear formation system is considered and a robust adaptive controller is designed to compensate for the effects of uncertainties or disturbances that the formation system may encounter. Although the control strategy is inspired by sliding mode control, it produces smooth control signals, thereby avoiding chattering. Also, an adaptation law is added such that the uncertainty or disturbance effects are effectively and quickly eliminated without a priori information about them. The combination of these two controllers guarantees that the satellite accurately tracks the optimal path in the unknown environment. Numerical simulations demonstrate the effectiveness and accuracy of the proposed two-step control methodology, in which a satellite formation is optimally reconfigured under unknown environmental disturbances.     

再入飞行器标称攻角优化设计

再入飞行器的标称攻角在弹道规划以及飞行器覆盖能力分析中起到重要作用,由于再入飞行中气动加热严重,过载和动压约束严格,给标称攻角的设计带来很大困难.针对弹道射面内最大纵程和最小总热载荷问题,在考虑热流、动压和过载约束下分别进行标称飞行攻角的优化设计.首先将过程约束转化为对控制量攻角的约束,将需要优化的标称攻角通过分段线性函数参数化,把最优控制问题转化为4个参数的寻优问题,然后利用遗传算法获得参数的初始猜想,并设计序列二次规划(SQP,Sequential Quadratic Programming)算法求解.仿真结果显示该方法能够快速获取再入标称飞行攻角,为再入轨迹优化和制导总体设计提供参考.     

Fuel-optimal and Energy-optimal guidance schemes for lunar soft landing at a desired location

Two guidance schemes (i) fuel-optimal (ii) energy-optimal to realize soft landing at a desired location on the moon are developed using the optimal control laws. The optimal control laws are obtained by solving a two-point boundary value problem formulated based on Pontryagin’s principle. The guidance laws, adapted from the optimal control laws, are obtained as a function of unknown co-state variables. Differential Transformation (DT) technique is employed to determine the unknown co-states at each time instant of landing trajectory using the information on the current vehicle state, target landing site (loaded on-board apriori) and the time-to-go. The numerical integration of co-state dynamics is avoided due to the DT based approach. The time-to-go, a critical parameter for any guidance scheme, is computed and updated real time using a simple strategy which uses the current and end states. The simple strategy for time-to-go works well even when the shape of the trajectory is nonlinear. Extensive analysis is carried out to evaluate and compare the proposed guidance schemes. Further, the proposed schemes are compared with other popular guidance schemes. The DT based proposed schemes help quantify the landing masses for fuel-optimal and energy-optimal objectives. Other features of the proposed schemes are that they do not assume constant gravity field and independent of reference trajectory.     

常值推力下面内轨道优化的一种改进间接法

针对常值推力下航天器面内轨道转移燃耗最省的轨道优化问题,利用极大值原理导出了最优轨迹下推力方向角应满足的控制方程,结合动力学方程建立了一种求解航天器面内最优转移轨道的改进间接法,及其在推力方向角调节能力受限条件下的应用方法。由于避免了协态变量微分方程组的求解,改进间接法相对于传统间接法降低了初值猜测的难度和计算量;与采用Gauss伪谱法求解相比,所建立的改进间接法求解结果精度更高,数值光滑性更好。仿真算例表明:推力方向角调节能力受限会改善推力方向角变化规律,降低推力方向角变化范围;就燃耗而言,推力越大燃耗越多,优化轨道节省燃耗更加显著。     

Optimal solar sail trajectory approximation with finite Fourier series

The heliocentric transfer of a solar sail-based spacecraft is usually studied from an optimal perspective, by looking for the control law that minimizes the total flight time. The optimal control problem can be solved either with an indirect approach, whose solution is difficult to obtain due to its sensitivity to an initial guess of the costates, or with a direct method, which requires a good estimate of a feasible (guess) trajectory. This work presents a procedure to generate an approximate optimal trajectory through a finite Fourier series. The minimum time problem is solved using a nonlinear programming solver, in which the optimization parameters are the coefficients of the Fourier series and the positions of the spacecraft along the initial and target orbits. Suitable constraints are enforced on the direction and magnitude of the sail propulsive acceleration vector in order to obtain feasible solutions. A comparison with the numerical results from an indirect approach shows that the proposed method provides a good approximation of the optimal trajectory with a small computational effort.     

航天器与空间碎片的混合编队重构控制与应用

针对空间碎片清理问题,提出了一种利用航天器与空间碎片混合编队队形重构控制技术捕获碎片的方法。首先,分析了地/月—日系L2拉格朗日平动点附近的限制性三体环境,并建立了编队卫星相对运动动力学模型;其次,提出了以太阳光压力作为航天器与空间碎片编队队形重构的控制力,实现各从星接近空间碎片的目的;最后,设计了基于线性二次型的最优控制器,并在Matlab/Simulink环境下进行仿真实验。仿真结果表明该方法可控制从星到达期望的位置(空间碎片的位置),且太阳帆板的姿态变化在可控范围内,进而证明了该方案可以应用于复杂空间环境下的碎片清理任务。     

高超末段机动突防/精确打击弹道建模与优化

针对高超声速飞行器再入末段机动突防、精确打击问题,从最优控制角度出发,提出了一种考虑拦截弹动力学特性的最优机动突防弹道优化方法,获得了高超声速飞行器的最大机动能力。该方法将拦截弹运动模型引入突防弹道优化的模型中,通过施加约束限制拦截弹按照比例导引律飞行。根据飞行任务和交战双方的弹道特点分段,结合各段的任务和特性,分别提出了突防性能指标和精确打击性能指标等,并通过加权函数将各个独立、矛盾的性能指标统一,建立了多对象、多段、多约束机动突防弹道优化模型,采用Radau多段伪谱法（MRPM）进行求解。针对该问题求解的初值敏感、可行域窄等问题,提出了一系列弹道优化策略,提高了收敛速度和求解精度,最终获得了最优机动弹道,并通过协态映射原理对其最优性进行了验证。结果表明,该方法能充分发挥高超声速飞行器的机动能力,获得满足落点精度要求的突防弹道,相对已有方法,将脱靶量提高了1~2个量级。灵敏度分析表明,该弹道对拦截弹的发射时间不敏感。      

Optimal spacecraft formation establishment and reconfiguration propelled by the geomagnetic Lorentz force

The Lorentz force acting on an electrostatically charged spacecraft as it moves through the planetary magnetic field could be utilized as propellantless electromagnetic propulsion for orbital maneuvering, such as spacecraft formation establishment and formation reconfiguration. By assuming that the Earth’s magnetic field could be modeled as a tilted dipole located at the center of Earth that corotates with Earth, a dynamical model that describes the relative orbital motion of Lorentz spacecraft is developed. Based on the proposed dynamical model, the energy-optimal open-loop trajectories of control inputs, namely, the required specific charges of Lorentz spacecraft, for Lorentz-propelled spacecraft formation establishment or reconfiguration problems with both fixed and free final conditions constraints are derived via Gauss pseudospectral method. The effect of the magnetic dipole tilt angle on the optimal control inputs and the relative transfer trajectories for formation establishment or reconfiguration is also investigated by comparisons with the results derived from a nontilted dipole model. Furthermore, a closed-loop integral sliding mode controller is designed to guarantee the trajectory tracking in the presence of external disturbances and modeling errors. The stability of the closed-loop system is proved by a Lyapunov-based approach. Numerical simulations are presented to verify the validity of the proposed open-loop control methods and demonstrate the performance of the closed-loop controller. Also, the results indicate the dipole tilt angle should be considered when designing control strategies for Lorentz-propelled spacecraft formation establishment or reconfiguration.     

Rapid design and optimization of low-thrust rendezvous/interception trajectory for asteroid deflection missions

Asteroid deflection techniques are essential in order to protect the Earth from catastrophic impacts by hazardous asteroids. Rapid design and optimization of low-thrust rendezvous/interception trajectories is considered as one of the key technologies to successfully deflect potentially hazardous asteroids. In this paper, we address a general framework for the rapid design and optimization of low-thrust rendezvous/interception trajectories for future asteroid deflection missions. The design and optimization process includes three closely associated steps. Firstly, shape-based approaches and genetic algorithm (GA) are adopted to perform preliminary design, which provides a reasonable initial guess for subsequent accurate optimization. Secondly, Radau pseudospectral method is utilized to transcribe the low-thrust trajectory optimization problem into a discrete nonlinear programming (NLP) problem. Finally, sequential quadratic programming (SQP) is used to efficiently solve the nonlinear programming problem and obtain the optimal low-thrust rendezvous/interception trajectories. The rapid design and optimization algorithms developed in this paper are validated by three simulation cases with different performance indexes and boundary constraints.     

Robust optimal sliding mode control for the deployment of Coulomb spacecraft formation flying

Inter-spacecraft electrostatic force (Coulomb force) is desirable for close formation flying control because of its propellant-less and free contaminate characteristics attributed to the propellant exhaust emission. This paper presents robust optimal sliding mode control to deal with the problem of thruster saturation in tracking the formation trajectory for Coulomb spacecraft formation flying. The robust controller design is based on optimal control theory as a linear quadratic system, and it is augmented with an integral sliding mode control technique. The stability of the closed-loop system is guaranteed using the second Lyapunov method. The developed controller outperforms the existing ones, because it has a higher degree of fine-tuning to cope with the uncertainty. Numerical simulations are employed to confirm the efficiency of the developed controller.     

Multi-agent Q-Learning control of spacecraft formation flying reconfiguration trajectories

This paper presents a novel approach based on multi-agent reinforcement learning for spacecraft formation flying reconfiguration tracking problems. In this scheme, spacecrafts learn the control strategy via transfer learning. For this matter, a new generalized discounted value function is introduced for the tracking problems. Due to the digital nature of spacecraft computer systems, local optimal controllers are developed for the spacecrafts in discrete-time. The stability of the controller is proven. Two Q-learning algorithms are proposed, in each of which the optimal control solution is learned on-line without knowledge about the system dynamics. In the first algorithm, each agent learns the optimal control independently. In the second one, each agent shares the learned information with other agents. Next, the collision avoidance capability is provided. The effectiveness of the presented schemes is verified through simulations and compared with each other.     

Convex optimisation approach to constrained fuel optimal control of spacecraft in close relative motion

This paper describes an interesting and powerful approach to the constrained fuel-optimal control of spacecraft in close relative motion. The proposed approach is well suited for problems under linear dynamic equations, therefore perfectly fitting to the case of spacecraft flying in close relative motion. If the solution of the optimisation is approximated as a polynomial with respect to the time variable, then the problem can be approached with a technique developed in the control engineering community, known as “Sum Of Squares” (SOS), and the constraints can be reduced to bounds on the polynomials. Such a technique allows rewriting polynomial bounding problems in the form of convex optimisation problems, at the cost of a certain amount of conservatism. The principles of the techniques are explained and some application related to spacecraft flying in close relative motion are shown.     

多操纵面飞机综合重构飞行控制方法

为提高飞机飞行安全性能,将被动重构控制与主动重构控制相结合,提出了一种综合重构飞行控制系统设计方法.采用轨迹线性化控制方法进行飞行控制律设计,使系统具有较强的鲁棒性.针对多操纵面飞机,提出面向故障的操纵面管理概念,通过基排序最优控制分配方法以基底的形式对操纵面进行重新组合,并按照期望控制目标进行排序,实现在操纵面故障下的控制重构.仿真结果表明,该控制系统能有效处理各种典型操纵面故障,在操纵面发生故障时仍能快速跟踪控制指令,保证较好的飞行性能.