High-fidelity trajectory optimization for aeroassisted vehicles using variable order pseudospectral method

In this study, the problem of time-optimal reconnaissance trajectory design for the aeroassisted vehicle is considered. Different from most works reported previously, we explore the feasibility of applying a high-order aeroassisted vehicle dynamic model to plan the optimal flight trajectory such that the gap between the simulated model and the real system can be narrowed. A highly-constrained optimal control model containing six-degree-of-freedom vehicle dynamics is established. To solve the formulated high-order trajectory planning model, a pipelined optimization strategy is illustrated. This approach is based on the variable order Radau pseudospectral method, indicating that the mesh grid used for discretizing the continuous system experiences several adaption iterations. Utilization of such a strategy can potentially smooth the flight trajectory and improve the algorithm convergence ability. Numerical simulations are reported to demonstrate some key features of the optimized flight trajectory. A number of comparative studies are also provided to verify the effectiveness of the applied method as well as the high-order trajectory planning model.     

绳系卫星轨道转移的最优控制

考虑系绳的弹性以及复杂状态和控制约束的作用,研究了绳系卫星面内轨道转移的最优控制问题。借助Gauss伪谱算法,将绳系卫星轨道转移的连续时间最优控制问题离散为大规模动态规划问题,进而利用非线性规划方法进行求解。通过数值模拟计算了子星最优转移轨道及最优控制力。结果表明:在满足相关约束的条件下,通过调节系绳张力可将子星从主星下方转移到上方的平衡位置,精确地实现子星轨道转移,并使得轨道转移过程呈现出良好的光滑性和对称性。最后基于协态映射定理对解的最优性进行了验证。     

航天器轨迹优化的通用数值方法

给出了航天器轨迹优化的一种通用数值仿真方法，该方法是由静态参数优化和动态参数优化构成，其中，静态参数优化采用可变误差多面体算法，动态参数优化采用基于最优控制理论的共轭梯度方法。     

Optimization method for mission analysis of aeroassisted orbital transfer vehicles

This paper is about an optimization method which has been developed to deal with trajectory optimization and mission analysis of an aeroassisted orbital transfer vehicle (OTV), in the context of preliminary design studies. Although this kind of trajectory can already be computed with existing trajectory optimization tools, we need a faster and robust tool which can be integrated as a “black box” in a multidisciplinary design process, in order to study rapidly many different OTV concepts and missions. In this context, our objective is not to get a very precise “optimal trajectory”, as existing “heavy” optimization tools do, but a solution precise enough to give a good insight of the performance (namely, the apogee altitude variation) and the mechanical and thermal loads. Incidentally, the solution obtained may also be used as an initial guess for a more precise trajectory optimization tool. To achieve this goal, we have studied parametric formulations of the control law, with optimization of the switching times. This development has been done considering a low lift-to-drag ratio vehicle (controlled only with the bank angle), like the aerocapture-designed version of the Mars Sample Return Orbiter. The cost function to be minimized is the heat flux, which is a key parameter for the multidisciplinary design of this kind of vehicle. The parametric formulation eventually chosen yields a good level of precision and robustness. Also, the study has been pushed further with the optimization of some mission parameters in the same process, in order to get directly preliminary answers to some trade-off issues in the mission analysis, like the choice of the initial perigee altitude.     

基于NFTET的高超声速飞行器再入容错制导

针对以X-33为对象的三自由度高超声速飞行器,采用相邻可行轨迹存在定理(NFTET)设计了容错制导律以解决再入段执行器发生故障的轨迹重构问题。在标称情况下采用预测校正算法生成满足再入过程约束和终端约束要求的再入轨迹;当执行器发生故障时,飞行器气动参数、结构和舵面力矩都可能发生不可预测的变化,原先的轨迹不再满足制导要求,因此需要设计新型容错制导律。针对实际再入制导模型,基于NFTET设计容错制导算法对轨迹进行重构,得到满足故障情况下制导任务的可行轨迹。从仿真结果中可以看出,容错制导算法生成的新轨迹重新回到了约束范围之内,轨迹呈收敛趋势,使得高超声速飞行器从故障恢复到正常飞行状态,提高了飞行器的自主容错能力。     

空天飞行器弹道/轨道一体化设计

弹道/轨道一体化设计是解决空天飞行器发射入轨和轨道转移问题的一种全新思路。针对目前存在的空天飞行器弹道/轨道一体化设计问题,通过改进非开普勒轨道方程的方法建立飞行器在连续推力、气动力、引力以及摄动力等多种力作用下的弹道/轨道一体化设计动力学模型;提出基于轨道设计反方法的弹道/轨道一体化设计方法。其创新点主要体现在：通过整合连续推力、气动力、引力以及摄动力等多种作用力达到了统一弹道/轨道模型的目的;提出了基于傅里叶级数形状方法的轨道设计方法,该方法相比于之前的逆多项式法,可以处理带推力约束的轨道设计问题;由于在弹道段采用类似于轨道设计反方法的设计思想设计弹道,使得弹道和轨道两段轨迹的设计方法也达到了统一,致使从模型和设计方法的角度都体现了弹道/轨道设计的统一性,解决了传统分段设计方法是在不同段采用不同的模型和方法,很难体现出一体化设计思想的问题。仿真分析表明本文提出的弹道/轨道一体化设计方法是可行和有效的。     

近地小推力转移轨道的加权组合制导策略

 研究了近地小推力转移轨道的制导问题，给出了一种基于局部最优控制律的自主制导算法。推导出了各改进春分点根数对应的局部最优控制律；通过最优推力分配和目标偏差两个策略，对各局部最优控制律进行动态加权组合，从而有效减少了制导律的设计参数。在此基础上，针对燃料最省转移轨道，定义了一种新的发动机开关函数。采用遗传/逐次二次规划混合优化算法计算了最优制导参数。与传统算法相比，该制导算法是一种闭环制导算法，能够实现飞行器的自主制导，并且制导过程中无需对制导参数进行更新。以地球低轨到高轨的小推力转移为例，采用该方法分别求解了时间和燃料最省转移问题，并与传统算法进行了比较分析。数值结果验证了该算法的有效性。     

航天飞机自动着陆轨迹优化设计

采用沿约束轨迹的推演计算方法及优化原理对航天飞机自动着陆飞行轨迹进行优化设计。首先根据过载、垂直下降率、连续性和平稳性要求确定飞行轨迹的基本形状，然后采用两点边值优化方法对轨迹参数进行优化设计以获得理想的接地动压。与传统的和制导系统联系在一起的轨迹设计方法相比，该方法具有简便、快速、可重复使用、鲁棒性强等特点。     

Optimal pure proportional navigation for maneuvering targets

The optimal pure proportional navigation (PPN) guidance law with time-varying navigation gains is considered. Unlike the conventional optimal PPN approach where linearized model was assumed in the optimization process, this work exploits the exact nonlinear formulation of PPN to derive analytically the optimal time trajectory of the navigation gain to minimize a performance index which is a weighted sum of the final time and the integral of the squared acceleration. It is verified that the PPN scheme with constant navigation gain is not only optimal in the vicinity of the interception point, but also optimal for the whole trajectory, if the navigation constant is designed by the methodology proposed here. Based on the optimization results for nonmaneuvering targets, a recursive optimal PPN scheme is proposed for maneuvering targets, wherein the optimal navigation gain and time-to-go are predicted recursively during the interception, and trajectory and performance of the interceptor guided by optimal recursive PPN scheme are evaluated analytically.     

基于HP-RPM的高超声速滑翔飞行器轨迹优化

针对含有航路点、禁飞区约束的再入突防轨迹优化问题,提出了基于HP自适应Radau伪谱法(HP-RPM)的分段轨迹优化策略,给出了在含有热流密度、过载、动压、航路点和禁飞区等约束条件下的再入轨迹优化模型;利用HP-RPM对含有再入的最优控制问题进行离散化,将其转为非线性规划问题,并根据航路点所在的位置,对再入轨迹进行分段,以再入滑翔飞行的时间最短为仿真目标函数进行仿真计算。仿真结果表明,此方法可以生成一条满足多种约束条件的高精度优化轨迹,并且用时较短。     

Observer-type Kalman innovation filter for uncertain linear systems

An observer-type of Kalman innovation filtering algorithm to find a practically implementable "best" Kalman filter, and such an algorithm based on the evolutionary programming (EP) optima-search technique, are proposed, for linear discrete-time systems with time-invariant unknown-but-hounded plant and noise uncertainties. The worst-case parameter set from the stochastic uncertain system represented by the interval form with respect to the implemented "best" filter is also found in this work for demonstrating the effectiveness of the proposed filtering scheme. The new EP-based algorithm utilizes the global optima-searching capability of EP to find the optimal Kalman filter and state estimates at every iteration, which include both the best possible worst case Interval and the optimal nominal trajectory of the Kalman filtering estimates of the system state vectors. Simulation results are included to show that the new algorithm yields more accurate estimates and is less conservative as compared with other related robust filtering schemes     

基于多模型预测的再入飞行器制导方法

在具有控制和轨迹约束的再入制导中,提出了一种新的基于多模型预测再入飞行器纵向制导方法。首先对再入飞行器的纵向运动方程沿着标准轨道线性化,基于分段仿射约束系统建模理论逼近原非线性系统,并采用线性矩阵不等式(LMIs)优化技术离线将保证闭环系统稳定的终端二次代价项求出;然后,在每个制导周期内,基于多模型预测在线求解构造好的有限时域优化目标,从而获得控制量的增量,并把第一个控制增量分量与标准轨道的控制量叠加后形成全量控制,用于实际再入轨道的制导,而在下一个制导周期基于新的状态参数重复上面的优化过程。数字仿真结果证实了所提方法的有效性。     

面内轨道转移过程中的绳系系统摆振特性研究

轨道转移过程中的绳系系统处于非开普勒轨道,导致系统呈现复杂的动力学行为并影响着星体的飞行安全,因此研究系统摆振特性具有重要的理论和实际意义。针对复杂的非线性和强耦合问题,利用动量矩定理建立绳系系统姿态动力学方程,以切向常值加速度轨道转移为任务背景,给出了系统质心运动轨迹;通过分析面内摆角的静态分岔现象,推导了面内、面外摆角的一阶摆动解析解;引入经典的珠点模型,研究系绳纵向和横向的振动特性,并分析了系绳摆动与系绳振动之间的耦合关系。仿真结果表明:面内轨道转移过程中,面内、面外摆角以固定的频率绕平衡位置做往返摆动,摆动频率大小以及平衡位置的变化均与系绳长度、推力加速度和所处轨道密切相关,面内摆角摆动频率接近轨道角频率时会引起共振现象,系绳在轨道转移过程中会出现大幅度横向振动等现象。     

基于单参数迭代的TAEM在线轨迹生成方法

针对大升力体轨道再入飞行器末端能量管理（TAEM）段制导控制能力强、末端约束不惟一的问题，将TAEM段分为动压跟踪和着陆预备2个阶段，设计了不同的纵向轨迹剖面，从而将TAEM段在线轨迹生成问题转化为单参数搜索问题。第1阶段设计标称动压剖面为纵向参考轨迹，使得飞行器过程约束得到保证。第2阶段纵向剖面设计为标称高度剖面，从而使得末端点高度和倾角约束得到保证。根据末端动压误差设计修正律，迭代修正第一阶段动压剖面，从而使得最终的纵向轨迹满足所有的状态约束。在线轨迹递推采用以时间为自变量的数值积分，递推过程引入闭环制导律，通过实时修正攻角跟踪纵向剖面，修正倾侧角跟踪地面轨迹，从而保证在线生成的轨迹符合物理特性，降低闭环制导难度。在考虑初期再入末端大范围状态散布情况下，数值仿真显示了所提算法的鲁棒性。     

高超声速飞行器周期巡航轨迹设计与优化

针对高超声速飞行器,给出了参数建模方法并分析了其气动特性,建立了乘波体外形高超声速飞行器的纵向平面运动学模型;提出了高超声速飞行器巡航段周期轨迹的方案及相应的轨迹优化方法.采用高斯伪谱法解决周期性轨道优化问题,将原有的连续周期非线性优化问题转化为多段离散优化问题,并采用SQP算法求得最优解.仿真结果表明,在巡航航程相同的情况下,周期轨迹比稳态轨迹更省燃料.     

垂直起降运载火箭返回轨迹不确定性优化

针对垂直起降运载火箭一子级在返回着陆的过程中存在的参数不确定性,提出了一种基于非侵入式多项式混沌展开的序列优化和可靠度评估的返回轨迹不确定性优化方法。首先,建立了返回多飞行段轨迹在确定性条件下的优化模型。然后,为同时兼顾轨迹的鲁棒性和可靠性,建立了由鲁棒最优目标函数、基于可靠度的路径约束和鲁棒等式约束组成的不确定性返回轨迹优化模型。最后,基于非侵入式多项式混沌展开方法对鲁棒目标函数和等式约束进行量化处理,将原随机鲁棒优化问题转化为高维状态空间中的等价确定性优化问题;为提高路径约束的可靠度评估效率,基于非侵入式多项式混沌展开方法对最可能点法进行改进,进一步发展了序列优化和可靠度评估策略。数值仿真结果表明,所提出的不确定性优化方法具有较好的鲁棒性,可以满足工程可靠性指标要求,同时还具有较高的精度和计算效率。     

On Optimal Allocation and Guidance Laws for Linear lnterception and Rendezvous Problems

This paper considers the problem of guiding N interceptors so that they intercept or rendezvous with N targets at the same time T. It is assumed that the interceptors and targets are described by linear time-varying differential equations (deterministic or stochastic). Under certain assumptions both the optimal (static or dynamic) interceptor-target allocation and the guidance of each interceptor for the optimal allocation can be obtained with modest real-time and storage computer requirements.     

基于周期平均的固定翼双旋弹弹道修正方法

针对一种滚转稳定的固定翼双旋弹,提出了一种平均控制力大小可调节的弹道修正方法,并据此设计了该弹的制导与控制方法。根据固定翼双旋弹的高速旋转特性,在对其角运动方程进行简化后进行角运动分析,得到了弹体的合法向力与固定鸭舵偏转角度之间的关系,然后基于周期平均的概念提出了一种弹道修正组件以不同转速、不同振幅旋转以产生大小可控平均法向力的弹道修正方法,并通过六自由度弹道仿真验证了此方法的可行性。分析结果表明,相比于弹道修正组件固定的传统弹道修正方法,这种新型的弹道修正方法可以通过连续地控制平均法向力的大小和方向,实现固定翼双旋弹的制导与控制,消除弹道偏差,提高命中精度。     

Direct Optimization of Low-thrust Many-revolution Earth-orbit Transfers

Low-thrust Earth-orbit transfers with 10^?5-order thrust-to-weight ratios involve a large number of orbital revolutions which poses a real challenge to trajectory optimization. This article develops a direct method to optimize minimum-time low-thrust many-revolution Earth-orbit transfers. A parameterized control law in each orbit, in the form of the true optimal control, is proposed, and the time history of the parameters governing the control law is interpolated through a finite number of nodal values. The orbital averaging method is used to significantly reduce the computational workload and the trajectory optimization is conducted based on the orbital averaging dynamics expressed by nonsingular equinoctial elements. Furthermore, Earth's shadowing and perturbation effects are taken into account. The optimal transfer problem is thus converted to the parameter optimization problem that can be solved by nonlinear programming. Taking advantage of the mapping between the parameterized control law and the Lyapunov control law, a technique is proposed to acquire good initial guesses for optimization variables, which results in enlarged convergence domain of the direct optimization method. Numerical examples of optimal Earth-orbit transfers are presented.     

基于轨道方程线性化的中段飞行轨道修正方法

对于深空探测、大范围远程交会、抵近观察等中段飞行时间较长的航天任务,中段轨道修正是完成任务的重要保证.基于摄动制导的思路,通过对轨道方程进行线性化,研究了二体条件下的中段轨道修正方法.在近地航天任务中,各种摄动力特别是J2项对终端位置精度的影响不可忽略,因此研究了考虑J2摄动影响时线性化方程的补偿方法.最后通过数值仿真...