基于DM-DSC的舰载机着舰自动复飞控制算法

针对理想复飞轨迹已知条件下的舰载机自动复飞控制问题,提出一种基于偏差模型的动态面控制（DM-DSC）算法。基于Radau伪谱法给出了舰载机着舰的最优复飞轨迹;根据得到的最优复飞轨迹及其所对应的控制方案,分别给出了速度子系统和高度子系统的偏差控制模型和反演（Backstepping）控制器,并通过引入动态面结构来获得虚拟控制量的微分信号,避免了Backstepping控制律求解过程中的“微分膨胀”问题;考虑到气动参数的不确定性及舰尾流场的干扰,采用线性扩张状态观测器（LESO）对控制模型中的干扰项进行估计和补偿,并设计抗饱和辅助系统来抑制控制饱和的不利影响;最后,基于Lyapunov方法证明闭环系统信号的有界性。仿真结果表明：所提算法具有良好的控制性能。     

电磁航天器编队动力学建模与运动规划方法

基于电磁航天器作用原理,利用拉格朗日方法,建立电磁航天器“绳系”动力学模型.基于偏差线性化动力学模型,以切向、径向编队为例,分析编队构型保持稳定性,设计构型保持控制律;将双星电磁航天器编队(EMFF,Electromagnetic Formation Flight)重构运动规划问题转化为标准优化问题,利用高斯伪谱优化方法进行求解.并提出序列控制策略,将该双星模型扩展应用于多星电磁编队构型重构问题,转化为多阶段运动规划问题利用多阶段优化方法进行求解.仿真结果表明本文的动力学建模方法和控制方法是可行的.      

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

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

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

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

基于自适应伪谱法的升力式飞行器火星进入段快速轨迹优化

针对升力式火星飞行器定点着陆任务的轨迹优化问题,给出了基于自适应伪谱法的快速优化算法.综合考虑探测器火星大气进入过程中的动力学约束、边界约束、路径约束以及控制约束条件,利用自适应伪谱法将轨迹优化问题转换为离散的非线性规划问题,采用序列二次规划算法进行求解,得到性能指标最优的进入轨迹.通过仿真验证,给出了实现火星进入过程燃料消耗最优的状态量和控制量轨迹.仿真结果表明,在Matlab中采用自适应伪谱法,能够在800s内采用267个配点,给出近似精度为10-6的火星进入过程中消耗能量最优的参考轨迹.      

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.     

基于预测校正法的低升阻比飞船再入轨迹设计

对低升阻比飞船高速再入地球大气层的轨迹规划问题,提出了基于数值预测校正NPC(Numeric Predictor-Corrector)法的预测制导方法,并用时间替代能量作为自变量从而规避奇异问题;由飞船再入的3自由度运动方程推导纵向平面内的运动方程,并使用预测-校正方法获得倾侧角数值;同时通过设计横向漏斗确定倾侧角的符号,最终得到当前需要的倾侧角.仿真结果表明,该方法能根据不同的再入条件得出对应的再入轨迹和倾侧角变化规律,其动压、热流率、过载也符合再入要求.     

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

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

平动点周期轨道间小推力转移的Gauss伪谱法

针对三体问题共线平动点附近周期轨道间的小推力转移问题,构造了一种新的形状函数,在此基础上提出了一种基于Gauss伪谱法的优化设计方法。首先,建立小推力轨道转移动力学模型,参考初始轨道和目标轨道的类型,构造一种新的形状函数以近似小推力转移轨道。为满足不同的约束要求,提出了振幅和相位按多项式变化的假设,推导了小推力转移轨道的近似解析解;然后利用Gauss伪谱法将小推力轨道转移的最优控制问题转化为非线性规划问题,并对推导的近似解析解进行解算和处理,为Gauss伪谱法求解非线性规划问题提供较为有效的控制变量的初始猜测值;最后以地月系统L1点附近Halo轨道间的小推力转移问题为例进行了仿真分析。仿真结果表明,小推力转移轨道近似解析解具备有效性和普适性,使得Gauss伪谱法的迭代效率提高55%以上,同时也表明Gauss伪谱法可有效解决平动点周期轨道间的小推力转移轨道优化设计问题。     

基于backstepping/RHO的变体飞机控制器设计

为了确保变体飞机在变体过程中的飞行稳定性,提出了一种并行修正控制器设计方案.采用backstepping方法进行标称控制器设计,提供基本的飞行稳定性和跟踪性能.使用雅可比线性化方法对变体飞机非线性方程进行线性化得到线性化模型,采用基于指令滤波器的滚动时域优化（RHO,Receding Horizon Optimal）方法进行修正控制器设计.在有限时域区间内实时解算修正控制量,对标称控制器进行补偿.从航迹角控制系统仿真结果可以看出,航迹角能够较好地跟踪指令信号,基本不受变体过程的影响,飞行控制系统满足实时性和鲁棒性的要求.     

基于GPU的快速有限元法求解密度场

为了快速计算分析利用视频测量方法测得的高速风洞试验密度场在扰动流场作用下的实验数据,针对密度场的数值求解问题,经过光线偏折理论分析密度场得到的二阶偏微分方程,对其研究实现了CPU串行有限元法求解。在此基础上提出了基于GPU的快速有限元求解密度场的方法,该方法经过对串行有限元法求解过程效率分析后,将耗时的神经网络拟合、总刚度矩阵和总载荷向量的求解进行了基于GPU的并行加速。实验结果表明:在精度满足实际工程要求的前提下,相对于CPU串行求解方法,所提方法可大大提高求解效率,且随着网格剖分成倍加密,其加速比成倍增加。      

叉流式热交换器动态特性的数值计算

建立了叉流式热交换器在热流体和冷流体都随时间变化的条件下的动态二维数学模型。使用一阶精度显式格式建立了热交换器的有限差分模型。该模型能够满足任意变化的温度及流量入口条件。经过实例计算表明,所得的理论热交换器出口平均温度变化曲线与相应的空中试飞测试曲线有较好的跟随性、吻合程度良好。     

直升机-吊挂耦合系统平衡特性和稳定性分析

基于单质点吊挂假设,建立了直升机-吊挂耦合系统的非线性动力学模型。在该假设下,吊挂将引入额外的自由度和约束,使运动方程增加4阶,且为隐式的微分代数方程组。通过将惯性力中的广义加速度项与广义速度的二次项分离,可将运动方程转化为显式的微分方程组。针对无吊挂的直升机本体和直升机-吊挂耦合系统这两种模型,采用直接数值方法,计算了它们以不同速度做零侧滑定直平飞的配平状态,并结合飞行试验数据进行对比。进一步对两种模型进行了小扰动线化处理,分析其运动模态并进行对比。结果表明,单质点吊挂会给系统引入两个新的运动模态,使得吊挂两个自由度上的运动与直升机高度响应产生耦合,同时改变直升机本体各个模态的特性,会使部分模态响应品质变坏。     

基于高斯伪谱法的化-电混合推进系统转移轨道优化设计

提出利用化-电混合模式推进系统完成地球同步卫星轨道转移任务,该推进系统极具应用前景,能够满足高有效载荷率、高入轨精度的工程实践需求。并针对基于该混合模式推进系统的转移轨道的优化方法展开研究,提出一种多阶段最优控制问题（OCP）的高斯伪谱法求解方法。该方法通过分段点的关联设置,将多个经高斯伪谱法转化而来的非线性规划问题转化为一个连贯的非线性规划问题(NLP)。两个不同算例的仿真过程与分析结果表明,该方法能够有效地解决多阶段非光滑连接轨道的优化问题,具有运算效率高、收敛性半径大、求解精度高等优点,可便捷地处理化 电混合模式推进系统的转移轨道优化设计问题。     

柔性杆子杆模型参数的遗传算法优化求解

研究了应用实数型的遗传算法解决柔性机器人子杆法建模中,子杆模型参数的优化求解问题.子杆模型参数形成高度非线性的方程组,它的求解是子杆法能否实现的关键.普通的牛顿迭代法及其变形,要求初始向量足够靠近解,一些所谓大范围收敛的算法又有方程组的Jacobian矩阵非奇异等要求,在此都难以应用.本文基于遗传算法的解法可求得模型参数的最优解,实例仿真显示了这一方法的可行性和高效率.      

升力体飞行器能量近似最优倾斜转弯飞行策略

升力体飞行器返回地球大气层内时受到热流、动压及过载等约束条件,为使得飞行器在倾斜转弯飞行过程中能量损失最小,需要研究一种能量最优的倾斜转弯机动飞行策略。本文的主要研究内容包括:从再入动力学模型出发,分析了升力体飞行器倾斜转弯的弹道特性,推导了终端速度与倾斜转弯幅度相关的解析解,提出了一种多约束条件下能量近似最优的倾斜转弯飞行策略;为进一步验证飞行策略的能量近似最优性,建立了能量最优的非线性轨迹优化模型,通过高斯伪谱法进行求解,获得能量最优的飞行轨迹。仿真结果表明,该飞行策略与优化方法获得的结果高度一致,并且该方法求解效率更高,工程应用性更强。      

基于伪光谱方法月球软着陆轨道快速优化设计

根据软着陆轨道的特征和优化算法的特点,对探月器软着陆轨道状态方程进行合理的转化处理,使其更适合优化数值解法求解。在此基础上,使用伪光谱方法将软着陆轨道优化问题转为一个约束参数优化问题,然后采用乘子法处理约束条件,采用变尺度法求解处理后的参数优化问题,仿真结果表明此方法具有收敛速度快、对初始控制量不敏感、鲁棒性强的优点。     

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

Influence of the 2:1 resonance in the orbits of GPS satellites

Due to the characteristics of their orbits the GPS satellites are submitted to the following main perturbations: terrestrial gravitational field, luni-solar gravitational attraction and solar radiation pressure (including the effects of the Earth's shadow). An additional perturbation arises due to the 2:1 commensurability of the orbital period of the satellite with the period of the Earth's rotation. An analytical theory is briefly presented to solve the equations of motion including the previously mentioned effects. The analytical solution, based on the Lie-Hori method, is compared with a numerical integration of the equations.