最优轨迹规划方法及其平动点航天器编队均衡耗能重构应用

由于均衡耗能航天器编队能够提高整体航天器编队服役时间,针对平动点航天器编队重构的均衡耗能最优轨迹规划问题,提出一种以状态、协态和控制三类变量插值为核心的求解非线性最优控制问题的新方法。基于连续时间表达的非线性最优控制问题通过变分原理转化为非线性方程组的求解,并进一步推导非线性方程组显示格式的Jacobi矩阵提高非线性方程组的计算效率。本文方法既满足最优控制理论的一阶必要条件又具有较大的收敛域;同时,不需要对协态初值准确猜测,避免了大规模非线性规划问题的求解。通过对中心航天器位置固定和无中心航天器两种情况的数值模拟,结果表明,本文方法对航天器编队重构轨迹规划问题能够达到均衡耗能的目标,具有一定的应用价值。     

非线性轨迹优化问题的保辛自适应求解方法

非线性轨迹优化问题一般是一个非线性最优控制问题。将非线性系统的最优控制问题导入到哈密顿体系的辛几何空间当中,基于对偶变量变分原理提出了求解非线性最优控制问题的一种保辛自适应方法。以时间区段两端协态作为独立变量,在时间区段内采用拉格朗日插值近似状态和协态变量,并利用对偶变量变分原理将非线性最优控制问题转化为非线性方程组的求解,保持了哈密顿系统的辛几何结构。并进一步,提出了基于多层次迭代的自适应算法,提高了非线性最优控制问题的求解效率。数值实验验证了该算法在求解非线性轨迹优化问题中的有效性。     

含J2摄动的有限时间航天器追逃博弈问题研究

航天器追逃博弈是航天器在轨捕获任务的一个重要问题,具有极高的军民两用双重价值．针对有限时间且考虑J2摄动的航天器追逃博弈问题,本文提出了一种精确的求解方法．该方法的核心思想是将有限时间的航天器追逃博弈问题建模为有限时间二人零和对策,则博弈中两航天器的最优控制策略可以转化为有限时间二人零和对策的鞍点解．在鞍点解的求解过程中,本文首先基于考虑J2摄动的非线性动力学方程,将两航天器动力学方程和始末边值条件与鞍点解必要条件结合得到两点边值问题,然后提出一种结合遗传算法和配点法的混合算法求解该两点边值问题以得到精确的鞍点解．本文利用数值仿真对所提方法的有效性进行了验证．结果表明：(i) 在航天器追逃博弈过程中,J2摄动对两航天器的最优控制策略具有较大影响;(ii) 所提方法能够精确求解出两航天器在有限时间的追逃博弈过程中的最优控制策略．     

多星编队构型重构全局优化策略1）--第七届全国空间轨道设计竞赛乙题解法

第七届全国空间轨道设计竞赛乙组题目以近地轨道卫星编队的燃料最优构型重构问题为背景,要求合理设计从星的相对飞行轨迹,构建尽可能多的目标构型,并考虑燃料消耗均衡分配问题.本文介绍了该题目的解法,包括问题分析、求解方案以及相关计算结果.由于燃料最优控制问题可用能量最优控制问题近似,文中采用能量最优连续推力优化模型近似计算燃料最优控制律.首先通过优化分析两两构型之间的燃料消耗,确定了构建4种目标构型的次序;然后综合考虑构型保持时间约束以及相位约束,采用混合PSO（particle swarm optimization）算法进行全局优化;最后将影响燃料消耗指标的从星相对飞行轨迹的连续推力段转化为bang-bang控制.     

多星编队构型重构全局优化策略——第七届全国空间轨道设计竞赛乙题解法

第七届全国空间轨道设计竞赛乙组题目以近地轨道卫星编队的燃料最优构型重构问题为背景,要求合理设计从星的相对飞行轨迹,构建尽可能多的目标构型,并考虑燃料消耗均衡分配问题.本文介绍了该题目的解法,包括问题分析、求解方案以及相关计算结果.由于燃料最优控制问题可用能量最优控制问题近似,文中采用能量最优连续推力优化模型近似计算燃料最优控制律.首先通过优化分析两两构型之间的燃料消耗,确定了构建4种目标构型的次序;然后综合考虑构型保持时间约束以及相位约束,采用混合PSO(particle swarm optimization)算法进行全局优化;最后将影响燃料消耗指标的从星相对飞行轨迹的连续推力段转化为bang-bang控制.     

非线性最小二乘跟踪的对偶变量变分方法

最小二乘跟踪方法是近几年提出的一种计算动力系统跟踪轨迹的方法.基于最小二乘跟踪的灵敏度分析算法可以有效避免传统的非线性系统灵敏度分析方法中的病态初值问题,因此其在混沌系统灵敏度分析方面有着重要的应用.针对非线性的最小二乘跟踪问题,首先将其重新描述为带有约束的非线性最优控制问题,引入协态变量并将系统的哈密顿函数表示为关于状态变量和协态变量的函数.然后将目标函数的积分时间离散化,根据对偶变量变分原理,以离散区间两端的状态变量作为独立变量,用Lagrange插值多项式近似离散区间内的状态变量和协态变量,进而将非线性最优控制问题转化为求解非线性方程组问题.这种算法无需对原问题做线性化处理,避免了复杂的线性化过程以及可能因此造成的误差,同时为求解非线性最小二乘跟踪问题提供了新的思路.根据最小二乘方法可以得到两条设计参数有微小变化的状态轨迹,基于这两条状态轨迹可进一步计算出系统关于设计参数的灵敏度,范德波振子作为数值算例验证了该方法在求解最小二乘跟踪问题以及计算非线性系统灵敏度时的有效性.      

航天器太阳阵伸展过程最优控制的遗传算法

本文讨论航天器太阳阵伸展中航天器姿态的最优控制问题。利用动量矩守恒原理导出带太阳阵航天器的动力方程,指出了系统的非」完整约束性质。将航天器太阳阵展开过程中主体姿态控制问题转化为非线性系统最优控制问题。厚优控制中引入遗传算法,替代传统的牛顿迭代方法。提出基于遗传算法的非完整运动规划最优控制算法。通过数值仿真表明,该方法对太阳阵伸展过程航天器姿态控制是有效的。     

带空间机械臂航天器姿态运动规划的数值算法研究

自由漂浮空间机械臂系统在无外力矩作用时,系统的动量矩守恒而成为非完整系统。利用这一特性本文研究了自由漂浮空间机械臂系统的三维姿态运动规划问题。导出带空间机械臂的航天器三维姿态运动数学模型,将系统的非完整运动规划问题转化为非线性系统最优控制问题,在最优控制中利用小波逼近控制输入规律,提出基于遗传算法的最优控制数值算法。通过数值仿真,表明该方法对带空间机械臂航天器系统的非完整姿态运动规划是有效的。     

航天器最优受控绕飞轨迹推力幅值延拓设计方法

针对航天器燃料最优可变周期绕飞轨迹的求解问题, 提出了一种以推力幅值为延拓参数的延拓方法. 问题求解从最为简单的双脉冲绕飞轨迹出发, 首先利用有限推力替代脉冲推力, 设定推力序列为“开— 关— 开”,然后逐步减小推力幅值, 最终得到最小推力绕飞轨迹; 此后, 再逐步增加推力幅值, 结合主矢量曲线判断最优推力开关序列, 将最小推力解延拓至有限推力以及脉冲推力燃料最优解. 该方法通过对推力幅值的延拓, 实现了有限推力bang-bang 控制与脉冲推力燃料最优绕飞轨迹优化问题的一并求解, 同时避免了最优控制问题中协态变量的随机猜测. 慢速与快速绕飞算例的优化结果验证了所提出方法的有效性.      

THRUST-AMPLITUDE CONTINUATION DESIGN APPROACH FOR SOLVING SPACECRAFT OPTIMAL CONTROLLED FLY-AROUND TRAJECTORY

针对航天器燃料最优可变周期绕飞轨迹的求解问题, 提出了一种以推力幅值为延拓参数的延拓方法. 问题求解从最为简单的双脉冲绕飞轨迹出发, 首先利用有限推力替代脉冲推力, 设定推力序列为“开— 关— 开”,然后逐步减小推力幅值, 最终得到最小推力绕飞轨迹; 此后, 再逐步增加推力幅值, 结合主矢量曲线判断最优推力开关序列, 将最小推力解延拓至有限推力以及脉冲推力燃料最优解. 该方法通过对推力幅值的延拓, 实现了有限推力bang-bang 控制与脉冲推力燃料最优绕飞轨迹优化问题的一并求解, 同时避免了最优控制问题中协态变量的随机猜测. 慢速与快速绕飞算例的优化结果验证了所提出方法的有效性.     

Optimal nonlinear feedback control of spacecraft rendezvous with finite low thrust between libration orbits

This paper presents the nonlinear closed-loop feedback control strategy for the spacecraft rendezvous problem with finite low thrust between libration orbits in the Sun–Earth system. The model of spacecraft rendezvous takes the perturbations in initial states, the actuator saturation limits, the measurement errors, and the external disturbance forces into consideration from an engineering point of view. The proposed nonlinear closed-loop feedback control strategy is not analytically explicit; rather, it is implemented by a rapid re-computation of the open-loop optimal control at each update instant. To guarantee the computational efficiency, a novel numerical algorithm for solving the open-loop optimal control is given. With the aid of the quasilinearization method, the open-loop optimal control problem is replaced successfully by a series of sparse symmetrical linear equations coupled with linear complementary problem, and the computational efficiency can be significantly increased. The numerical simulations of spacecraft rendezvous problems in the paper well demonstrate the robustness, high precision, and dominant real-time merits of the proposed closed-loop feedback control strategy.     

OPTIMAL MOTION PLANNING FOR A RIGID SPACECRAFT WITH TWO MOMENTUM WHEELS USING QUASI-NEWTON METHOD

An optimal motion planning scheme based on the quasi-Newton method is proposedfor a rigid spacecraft with two momentum wheels. A cost functional is introduced to incorporatethe control energy, the final state errors and the constraints on states. The motion planning fordetermining control inputs to minimize the cost functional is formulated as a nonlinear optimalcontrol problem. Using the control parametrization, one can transform the infinite dimensionaloptimal control problem to a finite dimensional one that is solved via the quasi-Newton methodsfor a feasible trajectory which satisfies the nonholonomic constraint. The optimal motion planningscheme was applied to a rigid spacecraft with two momentum wheels. The simulation results showthe effectiveness of the proposed optimal motion planning scheme.     

Optimal control of stretching process of flexible solar arrays on spacecraft based on a hybrid optimization strategy

The optimal control of multibody spacecraft during the stretching process of solar arrays is investigated,and a hybrid optimization strategy based on Gauss pseudospectral method(GPM) and direct shooting method(DSM) is presented. First, the elastic deformation of flexible solar arrays was described approximately by the assumed mode method, and a dynamic model was established by the second Lagrangian equation. Then, the nonholonomic motion planning problem is transformed into a nonlinear programming problem by using GPM. By giving fewer LG points, initial values of the state variables and control variables were obtained. A serial optimization framework was adopted to obtain the approximate optimal solution from a feasible solution. Finally, the control variables were discretized at LG points, and the precise optimal control inputs were obtained by DSM. The optimal trajectory of the system can be obtained through numerical integration. Through numerical simulation, the stretching process of solar arrays is stable with no detours, and the control inputs match the various constraints of actual conditions.The results indicate that the method is effective with good robustness.     

Optimal control of attitude for coupled-rigid-body spacecraft via Chebyshev-Gauss pseudospectral method

The attitude optimal control problem(OCP) of a two-rigid-body spacecraft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a discrete nonlinear programming(NLP) problem using the Chebyshev-Gauss pseudospectral method(CGPM). Solutions of the NLP problem can be obtained using the sequential quadratic programming(SQP) algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss(CG) points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform(FFT). To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and control variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.     

Symplectic multi-level method for solving nonlinear optimal control problem

By converting an optimal control problem for nonlinear systems to a Hamiltonian system,a symplecitc-preserving method is proposed.The state and costate variables are approximated by the Lagrange polynomial.The state variables at two ends of the time interval are taken as independent variables.Based on the dual variable principle,nonlinear optimal control problems are replaced with nonlinear equations.Furthermore,in the implementation of the symplectic algorithm,based on the 2N algorithm,a multilevel method is proposed.When the time grid is refined from low level to high level,the initial state and costate variables of the nonlinear equations can be obtained from the Lagrange interpolation at the low level grid to improve efficiency.Numerical simulations show the precision and the efficiency of the proposed algorithm in this paper.     

航天器集群编队最优单脉冲机动

对航天器集群编队最优单脉冲机动问题进行了研究. 针对不同的任务约束,基于非线性相对运动的周期性条件,以解析的思路分别研究了机动时刻给定和机动时刻未定情况下集群编队的最优单脉冲机动问题. 对于机动时刻给定的情况,从高斯变分方程和基于能量匹配条件的拉格朗日乘子法两个角度分别进行了探讨,将问题转化为对一元二次方程求极值或对一个单零点非线性方程求根;对于机动时刻未定的情况,将问题转化为对一个多零点非线性方程求根,通过傅里叶-贝塞尔级数展开可以得到任意高阶近似解. 对于每种情况,推导得到二范数意义下能量最省对应的最优参考长半轴,以及所施加的最优速度脉冲. 数值仿真验证了本文方法的正确性,并对仿真结果进行了解释和分析.      

HYBRID OPTIMIZATION STRATEGY FOR MOTION PLANNING OF SPACE ROBOT SYSTEM WITH PRISMATIC JOINT

研究了自由漂浮带滑移铰空间机器人非完整运动规划的最优控制问题,提出一种由高斯伪谱法求解可行解与直接打靶法求解最优解相结合的混合优化策略.首先,根据多体系统动力学理论建立空间机器人的动力学模型,给定系统的初始和目标位形,将空间机器人运动规划问题描述成博尔察（Bolza）型最优控制问题;然后,利用高斯伪谱法将最优控制问题离散为非线性规划问题,求解在较少勒让德-高斯（Legendre-Gauss,LG）点时状态变量和控制变量对应的可行解;最后,在LG点处离散控制变量,作为直接打靶法的初值,利用序列二次规划算法求解空间机器人系统的优化运动轨迹和最优控制输入.通过数值仿真,系统优化运动轨迹光滑平稳,最优控制输入也能很好地满足各种约束条件,仿真结果验证了该混合优化策略的鲁棒性和有效性.     

Control strategy of optimal deployment for spacecraft solar array system with initial state uncertainty

A control strategy combining feedforward control and feedback control is presented for the optimal deployment of a spacecraft solar array system with the initial state uncertainty. A dynamic equation of the spacecraft solar array system is established under the assumption that the initial linear momentum and angular momentum of the system are zero. In the design of feedforward control, the dissipation energy of each revolute joint is selected as the performance index of the system. A Legendre pseudospectral method (LPM) is used to transform the optimal control problem into a nonlinear programming problem. Then, a sequential quadratic programming algorithm is used to solve the nonlinear programming problem and offline generate the optimal reference trajectory of the system. In the design of feedback control, the dynamic equation is linearized along the reference trajectory in the presence of initial state errors. A trajectory tracking problem is converted to a two-point boundary value problem based on Pontryagin’s minimum principle. The LPM is used to discretize the two-point boundary value problem and transform it into a set of linear algebraic equations which can be easily calculated. Then, the closed-loop state feedback control law is designed based on the resulting optimal feedback control and achieves good performance in real time. Numerical simulations demonstrate the feasibility and effectiveness of the proposed control strategy.     

基于微分包含的绳系卫星时间最优释放控制

考虑系绳弹性的影响,建立了绳系卫星系统三维动力学模型,研究了在状态和控制约束下的绳系卫星非线性时间最优控制问题. 为缩减系统变量,控制律设计没有采用通常的状态空间模型,而是基于二阶微分包含,将连续时间最优控制问题离散为大规模动态规化问题,最后通过数值模拟验证了该方法的有效性.