满足多终端约束的二次曲线迭代制导方法研究

某些空间载荷会对入轨精度和入轨姿态同时提出很高的要求,应用于运载火箭控制系统的摄动制导方法的入轨精度无法满足要求,而传统的迭代制导方法无法约束终端的入轨姿态。为此提出了一种满足多终端约束的二次曲线迭代制导方法,该方法通过二次曲线形式规划整个真空段的飞行制导程序角,实时满足位置、速度与终端姿态约束,从而可以使得火箭以期望姿态角实现高精度入轨。算例结果表明,该方法能同时保证高精度的入轨指标和入轨姿态,并能适应偏差状态、约束姿态角变化与轨道根数小幅变化,具有较高的工程应用价值。     

基于T-H方程的多脉冲最优交会方法

针对椭圆参考轨道交会问题,采用T-H方程描述相对运动,提出了一种时间固定燃料最省的多脉冲最优交会方法,优化参数为交会脉冲及其施加时刻.当考虑到J2摄动或航天器初始相对距离较大时,用T-H方程进行状态预测其线性化误差一般不容忽略,而若用轨道积分预测则耗时较多,进而导致优化时间过长.针对此问题,提出一种采用前一优化脉冲节点的状态导出的轨道根数预测当前节点状态的预测方法.该方法简单实用,有效地加快了优化收敛速度.最后基于多脉冲优化解进行了数值轨道积分以验证交会精度.仿真结果表明,即使加上J2摄动,在初始相对距离为1 000 km时,该方法的终端位置精度仍能达到75 m.     

微小推力下小卫星的轨道递推与推力标定

针对采用微小推力进行轨道机动的小卫星,考虑复杂摄动力的基础设计了一种高精度轨道外推和推力在轨标定算法.首先,建立了考虑地球复杂摄动力和微小推力的小卫星轨道动力学模型;然后基于动力学模型,利用变步长龙格库塔算法,设计了对微小推力小卫星进行高精度轨道外推的方法.随后通过无迹卡尔曼滤波器(UKF),设计在轨标定算法,对存在误...     

深空探测转移轨道自主中途修正方法研究

针对深空转移轨道,提出以B平面参数为终端参数,采用脉冲控制和线性制导的自主中途修正方法.由自主导航系统定期确定探测器的当前位置和速度,之后利用精确动力学积分递推状态至B平面,得到打靶误差,若误差超出目标精度范围又在自主修正系统修正能力范围内,则立即利用以B平面参数为终端参数的线性制导公式并结合牛顿迭代计算出修正轨道的速度增量.利用中心差分公式计算终端参数对控制参数的敏感矩阵.蒙特卡罗仿真表明,在小偏差前提下该方法能够达到制导目标.     

基于夏氏最小二乘的轨道控制力系数辨识

在航天器轨道捕获、轨道维持和空间目标碰撞规避中都需要进行航天器轨道机动。针对航天器轨道机动过程中推力器的推力系数为装订常数,没有根据在轨工作实际进行优化而导致出现较大误差的情况,对控制力拟合系数进行辨识,作为修正控制参数以补偿轨道控制误差的依据,提高轨道控制精度。统计分析在轨管理的典型航天器平台及其发动机的轨道控制历史数据,分析轨道控制理论和在轨控制数据拟合建立轨道控制经验模型,用当前可测量的系统输入和输出预测系统输出的未来演变,得到不同工作情况下实际轨道控制误差与控制参数及其他主要影响因素之间关系的经验公式,为轨道控制策略决策提供参考。选取轨道半长轴控制量300m以上和300m以下的两类近地卫星,对其轨道控制历史数据进行分析,经实际数据测试,采用夏氏法进行推力系数拟合后预测的速度变化量精度较高。该种计算方法利用了轨道控制历史数据,计算方法简单,提高了轨道控制速度增量的预测精度,对轨道控制实施具有参考意义。     

长征运载火箭发射地火直接转移轨道研究

针对发射场射向范围和低温入轨级在停泊轨道上最长允许滑行时间严格受限下，工程设计火星探测任务发射轨道时，优化保障2~3周发射日期窗口的困难，以及精确设计拼接探测器分离点时刻轨道6根数耗时周期长等问题，提出了加入滑行时间限制再精确双向微分修正的设计算法（有别于传统B平面矢量法）。由地球影响球边界速度直接解出长征运载火箭入轨点轨道根数，将火箭飞行星下点弧段表示成滑行时间的分析式，快速找到满足诸多设计约束的初始发射轨道；并解决了考虑探测器小幅深空机动后火箭发射轨道的深入优化问题。精确力学模型下数值微分改进后，获得高精度计算结果，同等可比条件下运算精度不低于STK软件。实现了从地面起飞到抵达近火点的飞行轨道整体优化，确保了首次探火工程的顺利实施。      

空间交会目标航天器相位调整策略

目标航天器相位调整机动的目的是为交会对接任务做准备,是一个非固定时间交会问题.从调相的墓本原理出发,分析了大气阻力、升交点赤经漂移及轨道误差对相位角的影响,并将这些参数用于改进冲童对相位角影响的估算公式.将佑算公式与数值积分相结合,从摄动轨道机动参数计算、调相方向选取与机动圈次选择三个方面改进相位调整策略.仿真结果表明:提出的相位调整策略可以稳定高效地获得摄动轨道机动参数,能够满足终端精度要求,相对于直接根据相位角偏差正负来选择调相方向的策略可以节省推进剂消耗.     

一种地月转移轨道中途修正的制导算法

针对地月转移轨道中途修正问题，提出了一种求解修正速度增量的制导算法。该算法由初值设计和精确解求解两部分组成。首先，利用伪状态理论，通过简单迭代设计中途修正的初值，并通过Vinti预报方法修正了地球扁率的影响。然后，在求解精确解时，提出了一种基于伪状态理论的状态转移矩阵解析算法。该算法通过设计高精度的初值，降低了求解地月转移轨道中途修正问题的难度，而且避免了传统数值方法计算状态转移矩阵的复杂性。数值仿真结果表明，该算法可有效求解中途修正问题。     

中轨道Walker导航星座在轨备份方案优化设计

针对中轨道Walker导航星座在轨备份方案的优化问题，首先提出了在轨备份星轨位优化设计方法，考虑星座运行期间在轨备份星与工作卫星存在共同提供服务的情况，选取PDOP值和可见卫星数作为轨位优化指标，建立了轨位优化模型，并基于NSGA-Ⅱ算法对不同轨位下在轨备份星对导航星座服务性能的提升效果进行仿真分析；然后，基于在轨备份星轨位的优化结果，建立在轨备份星替换的轨道机动模型，并综合考虑速度增量和替换时间，确定了以替换时间最少为优化目标的在轨备份星替换方案。结果表明，提出的在轨备份方案能有效满足备份星的设计需求，增强导航星座的服务性能，实现故障卫星的快速替换，可为导航星座备份星的建设提供借鉴。     

严格回归轨道的管道导航方法研究

分析了作为参考轨道的严格回归轨道与卫星在轨运行状态的相对运动关系,提出近地遥感卫星的管道导航方法。由于参考轨道的设计只考虑高精度的地球非球形摄动,与在轨卫星的动力学环境存在差别,这导致两者之间存在切航向漂移。基于高精度的轨道动力学模型和位置确定方法,设计了卫星与参考轨道采样点的沿航向对齐算法,从而获取了卫星相对参考轨道采样点的相位时间偏差和卫星在参考轨道编队坐标系切航向平面内的相对运动轨迹,进而引入椭圆的“最小二乘适配法”获取相对运动轨迹的特征量。所研究的管道导航方法可应用于基于GNSS测量数据的卫星自主轨迹保持。     

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

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

Two-level optimization approach for Mars orbital long-duration,large non-coplanar rendezvous phasing maneuvers

A relative dynamics equation-set based on orbital element differences with J₂ effects is derived, based on which a two-level approach is proposed to optimize the Mars orbital rendezvous phasing with a large difference in the initial ascending node. The up-level problem uses the revolution deviation between the target spacecraft and the chaser as the design variable, and employs a linear search to find the optimum. The low-level problem uses the maneuver revolutions, locations, and impulses as the design variables, and is solved using a hybrid genetic algorithm combined with sequential quadratic programming. To improve the solution accuracy, an iteration method is developed to satisfy the terminal constraints of the absolute numerical integration trajectory. Test cases involving Mars sample return missions with large initial node differences are presented, which show that the relative dynamics, two-level optimization model, and hybrid optimization algorithm are efficient and robust. Compared with previously published results, the total velocity increment has been further reduced by utilizing this proposed approach. It is found that a five-impulse plan requires the least quantity of propellant, and a propellant-optimal minimum rendezvous duration exists for this long-duration, large non-coplanar rendezvous problem.     

载人月球极地探测定点返回轨道设计

针对载人月球极地探测任务，对定点返回轨道优化设计问题进行了研究。根据月球极地轨道的特性，介绍了三种返回轨道机动方案。结合三脉冲变轨方案，采用了从初步计算到精确计算的串行求解策略，对定点返回轨道进行优化设计。初步计算阶段，建立了基于近月点伪参数的三段二体拼接模型，将三脉冲机动段与月球逃逸段解耦，求解轨道初值；精确计算阶段，提出了两段拼接方法，分别进行逆向和正向高精度数值积分。经过仿真测试，验证了该策略求解的有效性和准确性。最后，通过大量的仿真计算，分析了定点返回轨道的特性。研究结论对未来载人月球极地探测定点返回轨道方案的设计具有重要的参考价值。     

Design and optimization of stable initial heliocentric formation on the example of LISA

Space-based gravitational wave (GW) detector such as the LISA (Laser Interferometer Space Antenna) mission requires high-precision and stability of the triangular formation. The dynamic environment of the detectors is complex, the science requirements for the formation are tight, and consequently, design and optimization of this high-standard formation with essentially many decision variables are very challenging. This paper studies the design and optimization of the stable initial formation of the heliocentric GW detector by taking the LISA as an example. The linearization method based on relative orbital elements is used for formation design in the two-body system. Three constraints are presented to reduce the number of decision variables to fourteen. The geometric features of the arm length and breathing angle of the triangular formation and the relative position of LISA to the Earth are analyzed and numerically verified in a high-fidelity dynamic model, from which the relationships of multiple requirements of LISA are studied, and a single index is summarized to simplify the optimization. Sobol sensitivity analysis is used to quantitatively evaluate the sensitivities of the decision variables to the cost function, with which a self-adaptive adjustment algorithm of the region of the variables is presented to improve the computational efficiency. The availability of the method to quickly and precisely find a stable initial formation in an extensive neighborhood of the nominal formation is verified by numerical simulation, where the best solution decreases about 47.54% of the arm length change from the requirement. This study shows that the initial formation should be deployed appropriately away from the Earth, and the gravitations of Venus and Jupiter should be utilized to maintain the formation stability.     

遗传算法在星座多星接近轨道优化中的应用

以三颗非共轨的Walker星座卫星为研究对象, 对航天器无需变轨与其接近的可能性进行研究. 将Lambert方法得到的航天器轨道作为初始轨道, 利用遗传算法对初始轨道进行优化. 对初始轨道在参考时刻位置和速度的改变量进行编码,形成对应的种群. 以航天器与星座卫星之间的最近距离为适应度函数, 通过种群的繁殖得到优化结果. 结合仿真算例, 分析了最小二乘算法和遗传算法在轨道优化中的优劣以及接近过程中轨道摄动的影响. 结果表明, 遗传算法适用于所提出的轨道改进问题. 研究结果可为单航天器无需变轨对星座多星接近问题提供理论依据.      

基于多目标优化与强化学习的空战机动决策

为了解决无人机自主空战中的机动决策问题,提出了一种将优化思想与机器学习相结合的机动决策模型。采用多目标优化方法作为决策模型核心,既解决了传统优化方法需要为多个优化目标设置权重的困难,又提高了决策模型的可拓展性;同时在多目标优化的基础上通过强化学习方法训练评价网络进行辅助决策,解决了决策模型在对抗时博弈性不足的缺点。为了测试决策模型的性能,以近距空战为背景,设计了3组仿真实验分别验证多目标优化方法的可行性、辅助决策网络的有效性以及决策模型的总体性能,仿真结果表明,决策模型可以对有机动的敌机进行有效的实时机动对抗。      

嫦娥一号卫星地月转移轨道中途修正分析

 中国第一颗月球探测卫星“嫦娥一号”于2007年10月24日成功发射并于11月7日顺利进入了距月面200km的科学探测的使命轨道,原计划在整个飞行过程中比较关键的轨道段是地月转移轨道段。要进行2至3次中途轨道修正,而实际的飞行结果只在第41h作了一次很小的修正,所用的速度增量是4.8m/s。基于有关的实测数据对此进行详细的分析,以期获得一些规律性的认识。     

基于线性协方差方法的交会对接误差分析

将线性协方差分析方法和蒙特卡罗仿真相结合,按交会任务和飞行特征把交会过程分为变轨飞行、自由飞行和中途速度修正三种特征段,研究了状态误差的传播规律和交会过程中各种误差对交会对接精度的影响。在变轨飞行段,分析了追踪航天器的姿态误差、控制系统性能状态估计误差,以及目标航天器轨道摄动对状态误差传播的影响。在自由飞行段,分析了追踪航天器估计状态误差的先验值和测轨误差对状态误差传播的影响。在中途速度修正段,分析了追踪航天器姿态误差和控制系统性能误差对状态误差传播的影响。仿真结果表明,误差分析方法设计合理,可以指导交会对接的轨道设计工作,能对已经设计好的交会策略进行误差分析和设计验证。     

再入飞行器末端能量管理段纵向剖面优化方法

摘要： 对于轨道再入飞行器,根据各阶段飞行特性和任务的不同被分为初期再入、末端能量管理以及进场着陆几个阶段.本文提出一种末端能量管理段航程及纵向剖面优选方法,通过对不同的初始航程和纵向剖面进行递推,以阻力板控制裕度最大为优化目标,选出最适应飞行器升阻特性的航程以及相应的高度动压剖面.考虑初始状态误差及气动特性偏差条件下的六自由度仿真验证所设计的纵向剖面的鲁棒性.