基于弹道逃逸和小推力捕获的地月转移轨道设计

针对地月空间货运任务和环月轨道空间设施建设任务,提出一种弹道逃逸和小推力捕获相结合的新型地月轨道转移模式,并建立了一整套该类型轨道设计方法。首先,在三体模型假设下分别建立地心弹道逃逸轨道和月心小推力捕获轨道的二维极坐标动力学模型。对于弹道逃逸轨道,将地心旋转系对准角和地月转移加速速度增量作为控制变量,提出初值估计解析公式,并应用序列二次规划算法进行快速求解。对于小推力捕获轨道,以月心距为参考量设置与弹道逃逸轨道的拼接点约束,提出能量匹配方法预估飞行时间,采用最优螺旋轨道的初始伴随状态解析式预估近月点伴随变量初值。基于混合法和轨道逆推思想,采用人工免疫算法进行小推力捕获轨道求解。仿真结果表明,基于弹道逃逸和小推力捕获的地月轨道转移方式大幅降低了近月制动燃料消耗,能快速穿越地球辐射带,且飞行时间适中;同时,提出的轨道设计方法能快速搜索到基于弹道逃逸和小推力捕获的地月转移轨道,验证了该方法的有效性。     

基于大幅值Lyapunov轨道的地月转移轨道设计研究

基于平面圆形限制性三体问题模型,利用与绕月轨道相切的大幅值Lyapunov周期轨道,提出了一种新的地月转移轨道设计方法。根据Poincaré截面与限制性三体问题动力学系统对称性计算得到的大幅值Lyapunov轨道,通过与绕月轨道拼接,将地月转移问题转化为地球到大幅值Lyapunov轨道的转移问题。为保证探测器能够从近地轨道（LEO）切向逃逸到达大幅值Lyapunov轨道,通过计算其稳定流形,采用最近点作为Poincaré截面的终止条件求解探测器的初始状态,并根据初始状态完成地月轨道的设计。仿真结果表明,该地月转移策略相比于Hohmann转移,在同样只需要两次速度增量的前提下,约节约100 m/s的速度增量,该研究为地月转移轨道的设计提供了一种新思路。     

航天器燃耗最优轨道直接/间接混合法延拓求解

针对转移时间和始末状态固定的航天器燃耗最优轨道的求解,给出了一种延拓方法:以双脉冲轨道为初值,首先求解全程推进轨道,然后逐步增加推力幅值,应用直接/间接混合法依次求解所有推力幅值下的、满足包括开关函数在内的所有必要条件的转移轨道,包括连续和脉冲推力轨道。通过基于开关函数曲线变化趋势的开关序列预设方法,以及基于已有优化结果的延拓步长自适应方案,实现了延拓方法的自动运行。为实现该延拓方法,给出了适用于改进春分点根数模型的脉冲最优转移轨道主矢量必要条件,推导了无推力轨道段改进春分点根数协态变量状态转移矩阵。通过3个算例对延拓求解会遇到的不同情况进行了具体说明。延拓方法可以看作现有直接/间接混合法的进一步完善与拓展,延拓过程和结果有助于对燃耗最优轨道与系统参数之间的关联获得更为深刻的认识。     

行星际小推力轨道Lambert解及应用

 为提高行星际小推力转移轨道初始设计精度,提出了基于N次逆多项式逼近的半解析Lambert算法,并基于该算法发展了一种转移轨道初始设计方法。首先,采用N次逆多项式近似小推力轨道形状,应用推力方向假设和位置速度边界条件推导出部分系数及推力大小解析式。接着,分析了飞行时间约束和轨道动力学约束下解的存在性,并给出了关键系数的可行域。然后,利用探测器质量消耗方程建立了Lambert问题求解模型并加以解决。最后,基于所提Lambert算法,通过对连续推力约束进行降维,提出一种求解多圈非固定时间的行星际小推力转移轨道初始设计方法。分别以固定和非固定时间转移任务为例对所提Lambert算法和初始轨道设计方法进行了数学仿真,数值结果表明：相比传统6阶方法,所提Lambert算法在目标轨道半长轴为5 AU时可减少速度增量需求36.63%;所提初始设计方法与最优化方法设计结果接近,可为转移轨道的精确设计提供可行的设计初值。     

连续地月转移系统动力学研究与能量分析

为了研究新型连续地月转移系统的动力学及能量需求,采用Lagrange方法,在系绳为刚性杆假设的前提下,同时忽略第三体引力、地球扁率和系绳轴向变形等扰动因素的影响,建立了驱动型动量交换绳系卫星(MMET)系统的三维刚性动力学模型。对所建立的动力学模型进行了数值仿真及对比分析,仿真结果验证了所建模型的正确性。研究表明,外力矩对系统轨道运动参数影响甚小,对姿态运动参数影响明显。采用MMET方式进行载荷转移,推导出了实现载荷地月轨道转移所需的入口速度条件以及时间周期条件,并求解出了载荷在2次任务之间的时间间隔。给定初始条件下,当MMET系统以0.231 6 rad/s的旋转角速度绕其质心旋转1 448.5圈,其绕地心刚好运行5圈时,载荷可顺利进入地月转移轨道。最后,对连续地月转移系统实现载荷的地月转移进行了能量对比分析,结果表明,相同条件下,MMET载荷转移方式相比于传统脉冲变轨方式在载荷转移过程中消耗更少的能量。     

Optimization of multiple-impulse minimum-time rendezvous with impulse constraints using a hybrid genetic algorithm

The minimum-time multiple-impulse rendezvous with impulse constraints is investigated in this paper. Based on the Clohessy–Wiltshire (C–W) equations, an optimization model including several different kinds of impulse constraints such as the maximum impulse magnitude, the total velocity change magnitude and the time of imposing impulse for multiple-impulse minimum-time rendezvous is established. A generalized inverse matrix solution for linear equation is applied to avoid handling the terminal equality constraints. In order to obtain the global solution efficiently, a hybrid optimizer combining the advantages of a floating-coded genetic algorithm and simplex method is employed. A low-earth orbit multiple-impulse rendezvous problem is used as an example. The influence of the number of impulses, the optimization variables and the constraints on the solution is analyzed, and the different optimization algorithms are compared. Results indicate our proposed model and approach is effective in designing linearized minimum-time rendezvous trajectory with impulse constraints.     

基于解析梯度的微推进转移轨道优化方法

 研究基于离散脉冲策略的行星际微推进转移轨道优化问题,为改善收敛性和计算效率,提出了一种解析的梯度矩阵计算方法。首先,基于离散脉冲思想建立了行星际微推进转移轨道优化模型,通过对速度脉冲矢量进行参数变换,避免了传统优化模型初值猜测和搜索域难以界定的缺点;然后,基于变分原理分析了离散脉冲模型中轨道状态转移矩阵与状态变分之间的关系,并给出了多种轨道特征条件下状态变分的计算方法,该方法可扩展到多次状态不连续情况;在此基础上,推导了离散脉冲轨道优化模型中性能指标与轨道约束对寻优参数梯度矩阵的解析表达式。以地球-火星的燃料最省微推进转移为例对所提解析梯度方法进行了仿真验证。数值结果表明：本文推导的解析梯度矩阵是正确的,与传统数值梯度计算方法相比,解析梯度矩阵可有效改善寻优算法的收敛特性,并能够提高计算效率约47%。     

星际小推力转移轨道快速设计方法

 针对星际探测任务中燃料最省小推力转移轨道问题，提出一种基于标称轨道的快速设计方法。首先以具有相同端点时刻的无摄动标称轨道为参考，对传统的非线性轨道优化模型进行合理变换，将复杂的优化问题简化为可解的两点边值问题；然后基于标称轨道3个独立积分推导出解析的状态转移矩阵，并以此为基础导出了两点边值问题的最优解析解。该方法无需数值迭代，有效地克服了数值优化方法收敛性差、计算效率低的缺点。最后，以探测火星的小推力转移轨道为例对该方法进行了验证，与精确的数值结果相比，该方法计算的燃料消耗误差小于1%。     

小卫星星座的最优脉冲变轨研究

针对小卫星星座,进行星座发射中的最优脉冲式变轨研究,给出了形成星座的脉冲式变轨的基本原理;基于卫星相对运动状态转移方程,推导出了星座参脉冲式变轨的理论解,即所要施加的脉冲控制量的解析式,利用遗传算法,对双脉冲式变轨的脉冲控制量进行了优化计算,求得了使总变轨脉冲最小的最优变轨时间,最后,探讨了星座脉冲式变轨的工程实现现途径,为工程应用和研究提供参考。     

一种基于轨道根数约束的最优制导方法

针对航天器空间变轨的制导问题,研究了一种基于轨道根数约束的最优制导方法。在地心惯性坐标系下直接建立航天器的最优控制模型,给出了位置速度表达式和协态变量初值之间的关系;进一步,在不约束真近点角的前提下,推导了5个轨道根数的约束方程,并通过对轨道根数和最优控制理论中协态方程的特性分析,获得了另外两个约束方程。协态变量初值可直接通过求解7个完整约束方程组获得,进而得到最优推力方向。仿真验证了所提制导方法的有效性。     

Fuel-optimal trajectories of aeroassisted orbital transfer withplane change

The fuel-optimal control problem arising in noncoplanar orbital transfer employing aeroassist technology is addressed. The mission involves the transfer from high Earth orbit to low Earth orbit with plane change. The complete maneuver consists of a deorbit impulse to inject a vehicle from a circular orbit to an elliptic orbit for atmospheric entry, a boost impulse at the exit from the atmosphere for the vehicle to attain a desired orbital altitude, and a reorbit impulse to circularize the path of the vehicle. In order to minimize the total fuel consumption, a performance index is chosen as the sum of the deorbit, boost, and reorbit impulses. The application of optimization principles leads to a nonlinear, two-point, boundary value problem, which is solved by a multiple shooting method     

月球软着陆制动段飞行轨迹与制导律研究

对自环月轨道开始的月球软着陆制动段飞行轨迹和制导律进行了研究。建立了统一的均匀球体三维软着陆模型,采用燃料次优解析制导方法,通过单步优化获得了局部最优的推力角解析式。同时,在推力角中引入前馈项,用于消除初始位置和速度偏差对软着陆的影响。对制动段飞行轨迹的仿真分析验证了制导律设计的有效性。     

电推进在静止轨道空间碎片减缓中的应用

研究了电推进在静止轨道空间碎片减缓中应用的可行性和效果。通过分析电推进中离子推进的技术特点、发展和应用,阐明了用电推进将静止轨道卫星在寿命末期移高到“垃圾轨道”是可行的。提出了电推进采用沿飞行方向的常值推力,以多次双脉冲变轨方式工作的方案,并通过仿真计算,得出了该方案中推力大小、工作时间、特征速度、燃料消耗量等参数的确定方法和量化性指标,验证了电推进应用在空间碎片减缓中的作用和效果。探讨了电推进以脉冲方式工作的工程实现方案．     

Earth-moon Trajectory Optimization Using Solar Electric Propulsion

The optimization of the Earth-moon trajectory using solar electric propulsion is presented. A feasible method is proposed to optimize the transfer trajectory starting from a low Earth circular orbit (500 km altitude) to a low lunar circular orbit (200 km altitude). Due to the use of low-thrust solar electric propulsion, the entire transfer trajectory consists of hundreds or even thousands of orbital revolutions around the Earth and the moon. The Earth-orbit ascending (from low Earth orbit to high Earth orbit) and lunar descending (from high lunar orbit to low lunar orbit) trajectories in the presence of J2 perturbations and shadowing effect are computed by an analytic orbital averaging technique. A direct/indirect method is used to optimize the control steering for the trans-lunar trajectory segment, a segment from a high Earth orbit to a high lunar orbit, with a fixed thrust-coast-thrust engine sequence. For the trans-lunar trajectory segment, the equations of motion are expressed in the inertial coordinates about the Earth and the moon using a set of nonsingular equinoctial elements inclusive of the gravitational forces of the sun, the Earth, and the moon. By way of the analytic orbital averaging technique and the direct/indirect method, the Earth-moon transfer problem is converted to a parameter optimization problem, and the entire transfer trajectory is formulated and optimized in the form of a single nonlinear optimization problem with a small number of variables and constraints. Finally, an example of an Earth-moon transfer trajectory using solar electric propulsion is demonstrated.     

基于TH方程的椭圆轨道交会偏差分析

研究了考虑闭环控制偏差的椭圆轨道交会问题。基于推导的状态转移矩阵,采用线性协方差分析方法对交会椭圆轨道目标过程的导航偏差和控制偏差进行了分析,建立了闭环控制状态偏差和协方差的递推模型,给出了一种改善交会制导精度的修正算法。通过蒙特卡洛打靶仿真对线性协方差分析结果进行了验证,结果表明,本文所提出的偏差分析方法和修正算法精度较高,可用于椭圆轨道交会问题的研究。     

An epoch state filter for use with analytical orbit models of low earth satellites

Kalman filters provide a well established means for satellite orbit determination. In combination with space based sensors like GPS, DORIS or PRARE, accurate estimates of the spacecraft position and velocity can be obtained in real-time on-board the space vehicle. Traditionally, numerical methods of varying complexity are applied for propagating the state vector between measurements and updates of the state vector are referred to the epoch of the latest sensor output. In the present study a different approach is followed, which offers increased on-board autonomy and is particularly promising for small satellites with moderate accuracy requirements. An analytical orbit model is used to describe the spacecraft trajectory and mean elements at epoch are estimated instead of the instantaneous, osculating state vector. This adds the capability of performing on-board orbit prediction over time scales of up to one week, which is required, for instance, for the autonomous forecast of eclipse times or station contacts. Making use of the SGP4 orbit model that is compatible with NORAD twoline elements, an epoch state Kalman filter has been implemented and tested with GPS flight data of GPS/MET (MicroLab-1) and MOMSNAV (MIR). It is demonstrated that the proposed method provides an accuracy compatible with that of the analytical model and is robust enough to handle large data gaps in case of limited on-board resources for GPS operations. By adjusting the ballistic coefficient along with the mean elements, a considerable improvement of mid-term orbit predictions is achieved over methods that are restricted to the estimation of the state vector alone.     

基于连续同伦算法的天基仅测角初定轨方法

要进一步提高天基短弧初定轨的精度,在观测资料精度较高的情况下,仅考虑二体问题是不够的,还应考虑轨道摄动的影响。因此,基于无摄初轨的单位矢量法原理和矢量斜分解方法,给出了考虑摄动的天基仅测角初定轨单位矢量法。针对天基仅测角观测条件方程组求解过程中易出现迭代不收敛或收敛到平凡解的问题,引入连续同伦算法求解观测条件方程组,提出了单星观测方式下的空间目标天基仅测角初定轨方法,并通过数值仿真算例验证了该算法在较大范围的收敛性和数值稳定性。     

Fuel-optimal deorbit scheme of space debris using tethered space-tug based on pseudospectral method

This paper proposes a fuel-optimal deorbit scheme for space debris deorbit using tethered space tug. The scheme contains three stages named respectively as dragging, maintenance and swinging. In the first stage, the tug, propelled by continuous thrust, tows deorbit to a transfer orbit with a tether. Then in the second stage, the combination of the tug and the debris flies unpowered and uncontrolled to a swing point on the transfer orbit. Finally, in the third stage, the tug is propelled at the swing point and the rotation speed of the tethered system increases such that the debris obtains enough velocity increment. The trajectory optimization of the first stage is established considering the total fuel consumption of the three stages, whereas the dynamic model is simplified for computation efficiency. The solution to the optimal problem is obtained using a direct method based on Gauss pesudospectral discretization. Then a model predictive controller is designed to track the open-loop optimal reference trajectories, reducing the states’ deviations caused by model simplification and ignorance of perturbations. Furthermore, it is proved that the fuel-optimal swing point is the apogee of the transfer orbit. The paper analyzes the fuel consumption of a typical scenario and demonstrates effectiveness of the proposed deorbit scheme numerically.     

单点弹道精度与轨道半长轴远地点高度计算精度的映射关系

在仅使用单点位置、速度信息计算轨道的奈件下,针对轨道半长轴、远地点高度的精度问题,在轨道面内,应用活力公式和二体运动学理论推导得出了轨道计算精度与弹道测量精度间映射关系的解析表达式,并采用数值分析方法给出了不同的位置、速度误差与半长轴、远地点高度最大误差之间的数值关系.仿真结果表明,对于位置误差和速度误差大小分别为100 m和1 m/s的算例,半长轴最大误差和远地点高度最大误差分别约为2 km和4 km.基于此方法,可以将弹道误差传递至轨道参数误差,进一步分析故障误判和漏判概率;也可根据轨道参数精度要求反算弹道测量精度要求,以作为地面测量系统建设的技术依据.     

基于速度增量的空间绳系机器人中距离逼近过程最优轨迹规划

针对空间绳系机器人中距离逼近过程最优轨迹规划问题,提出了基于速度增量的多目标逼近轨迹优化方法,优化指标为总速度增量及逼近时间。首先建立逼近过程相对动力学模型及最优逼近轨迹优化模型,然后利用改进型非劣分类遗传算法得到相对逼近距离1.5 km内逼近轨迹的Pareto最优解。仿真结果表明,该方法可以揭示空间绳系机器人逼近距离1.5 km内逼近时间、燃料消耗、相对目标的面内视界角及速度增量次数之间的相互关系,能满足针对不同任务需求提供相应最优轨迹的要求。