月地转移轨道快速设计与特性分析

对采用直接大气再入方式的月地转移轨道,考虑大气再入界面参数和地面落点位置约束,提出了一种基于双二体模型的快速设计方法。该方法分为内外两层迭代循环,内层循环使月心段轨道和地心段轨道在月球影响球边界处连续,并采用Lambert问题与Newton-Raphson法相集合的方法求解满足再入角约束的地心段轨道参数;外层循环通过调整地心段轨道倾角和轨道置入时间使月地转移轨道满足地面落点位置约束。分析表明,存在四种类型的月地转移轨道满足大气再入界面约束,分别为降 降型、降 升型、升 降型和升 升型。在此基础上,对四种类型月地转移轨道的近地点地心距、置入分布点、再入点分布等特性进行了分析。仿真结果验证了所提出方法的有效性。     

地球同步轨道卫星混合推进转移轨道特性分析

以地球同步轨道卫星转移轨道设计为背景，针对全化学推进燃料消耗大和全电推进转移时间长的问题，开展了化学 电混合推进转移轨道优化设计与特性分析。首先，讨论了轨道倾角和近地点幅角变化对混合推进转移轨道的影响。研究表明，在混合推进优化设计中需要将轨道倾角作为优化变量之一。然后，以近地点半径、远地点半径、轨道倾角为优化变量生成搜索网格，得到过渡轨道集。针对每条过渡轨道，构建化学推进转移段和电推进转移段。其中化学推进段采用单圈兰伯特转移解算，电推进段采用混合法优化。最后，以燃料消耗和转移时间为指标，在搜索域内开展解算分析，研究了混合推进轨道在整个搜索域内的变化趋势。该方法可以提供具有不同燃料消耗和转移时间的混合推进转移解集，拓宽了解空间，可供轨道设计人员根据任务约束灵活选用。     

圆及椭圆轨道释放绳系子星的方法研究

研究绳系卫星系统质心沿椭圆轨道运动时，子星释放的最佳初始条件，运用变量替代简化动力学控制方程；对圆轨道运动的母星，采用轨道转移方法释放子星，计算弹射变轨参数及过渡轨道终了真近点角。本文取消母星质量远大于子星质量的假设，考虑子星释放过程中绳系卫星系统质心随子星运动的变化，所得结果更具有普遍意义。数值模拟结果表明，所述方法的有效性。     

航天器交会中的Lambert问题

应用Lagrange转移时间方程研究空间交会中的Lambert问题,包括经典Lambert问题(飞行弧段不足一圈的椭圆型轨道转移)与多圈Lambert问题(飞行圈数超过一圈的轨道转移),阐述转移轨道的几何特性与转移轨道类型,分析转移时间与转移轨道参数及变轨速度增量之间的关系。对航天器交会中常用的圆轨道之间的双冲量转移,给定转移角与转移时间,阐述最小变轨速度增量所对应的转移圈数与轨道参数的求解方法,提出满足最小变轨速度增量要求的轨道转移的图解法。对给定的初始分离角与交会时间,按最小变轨速度增量要求,确定航天器交会的初始漂移时间、双冲量轨道转移时间与终端停泊时间。     

一种小推力借力飞行转移轨道初始设计方法

提出一种基于正弦指数函数的小推力借力飞行转移轨道初始设计方法。建立极坐标形式的正弦指数函数表达式,用于模拟小推力转移轨道,并用绕圈参数和飞行路径角对转移轨道的参数进行表征;在约束条件下对转移轨道参数进行离散化处理,求解转移轨道与目标星投射轨道在参考面内的交点,计算到达目标轨道时刻的极角,进而得出小推力转移轨道的初始设计参数;设计了地球—木星的火星借力小推力转移轨道,仿真结果验证了该方法在小推力转移轨道初始设计中的快速性与准确性。     

首次整星出口水平全面提升尼日利亚通信卫星-1上天并定点

□□2007年5月14日00:01,我国承制的尼日利亚通信卫星-1(NIGCOMSAT-1)由长征-3B火箭发射至近地点为209km、远地点为41951km、轨道倾角为25°的超地球同步转移轨道.这次发射是根据我国与尼日利亚签订的尼日利亚通信卫星-1轨道交付合同实施的,也是我国首次以火箭、卫星及发射支持的整体方式,为国际用户提供商业卫星服务.     

月球探测器发射机会分析

发射月球探测器实际上是使探测器与月球交会的问题,由于月球位置的变化,选择不同的交会日期相应的地月转移轨道是不相同的。文章分析了月球在一个恒星月内位置变化与相应的地月转移轨道升交点赤经和近地点幅角变化的关系,在这个基础上讨论了发射机会和发射窗口问题。     

中星6B通信卫星成功发射：飞吧!让太空纪录作古

2007年7月5日傍晚的西昌，晚霞消退，暮色渐浓。20点08分，长征三号乙运载火箭在巨大的轰鸣声中踏上了长征系列运载火箭的第101次征程。火箭拖着长长的桔红色火焰，托举起中星6B通信卫星飞向浩渺太空。约26分钟后，从西安卫星测控中心传来消息：星箭分离，卫星进入近地点高度206千米、远地点高度50030千米、轨道倾角24．3度的地球同步转移轨道。     

火卫一周期准卫星轨道及入轨分析

围绕火卫一的准卫星轨道（QSOs）因其具有良好的稳定性，是火卫一探测任务最为实用的轨道。在平面圆型限制性三体问题模型下，利用庞加莱截面和KAM环迭代方法探究了准卫星轨道的周期轨道族，并给出不同能量准卫星周期轨道的初始条件。针对火卫一周期准卫星轨道入轨，提出一种转移轨道设计方法:对准卫星周期轨道调整速度后进行反向积分，直至离开火卫一邻近区域，从而得到由火星环绕轨道向火卫一周期准卫星轨道的转移轨道，并调整转移轨道参数对燃料与时间消耗进行优化。研究结果表明，当周期准卫星轨道能量处于特定区间时，存在特定速度脉冲区间，可利用火卫一引力实现较少燃料消耗的轨道转移；在该速度脉冲区间中，通过选取较小的速度脉冲，可缩短转移时间。      

飞向晕轨道的探测器轨道优化

晕轨道的稳定流形为从地球到晕轨道的转移轨道设计提供了便利．以往都采用在晕轨道上的目标点施加脉冲,这样,稳定流形只是为转移轨道的设计提供一个初始猜想,探测器并没有运行在稳定流形上,因而并未真正利用稳定流形节省燃料的优势．利用基于序优化理论的微分修正法,研究从晕轨道近地点稳定流形上不同点进入稳定流形所需要的燃料消耗,寻找燃耗最少的转移轨道．仿真表明,对于晕轨道近地点入轨,找到的稳定流形射入点机动比以往的晕轨道入轨点机动节省约33％的燃料消耗．此外,还对晕轨道上不同入轨点的入轨代价进行了研究,得到了晕轨道近地点入轨的最小燃耗解．     

月地转移轨道快速设计方法

月地转移轨道设计一般分为初步轨道设计和精确轨道设计.其中,初步轨道设计的准确性是确保后续精确轨道设计收敛的关键.提出了一种基于Lambert算法的月地转移轨道快速设计方法.以出月球影响球的时刻、位置和速度为中间变量,将轨道分为地心段和月心段分别进行计算.将探测器飞出月球影响球至指定再入点的地心段轨道简化为一个Lambert问题进行求解,提出了通过牛顿迭代法求解月地转移轨道Lambert问题的方法,避免了Lambert问题求解时大量的超几何函数和级数计算,提高了计算效率.在月心段轨道的快速计算中,提出了根据探测器出影响球速度矢量、月球停泊轨道倾角和近月点高度计算月心双曲线轨道根数的新方法.通过迭代计算,使得两段轨道在月球影响球处的位置和速度连续,从而获得一条完整的满足两端约束的双二体月地转移轨道.该方法计算速度快,精度相对较高.计算结果可以作为后续精确轨道设计的初值.      

地月自由返回轨道设计

载人月球探测任务需要着重考虑任务的安全性,能否“自由返回”是很重要的一项指标。在轨道设计上存在两项要求,首先是地月转移,然后是经月球近旁转向并判断能否返回地球,即实现飞越月球以后的无动力返回。文章结合精确的行星历表(DE405),首先比较不同力模型中同一条轨道的返回情况,以确定初始轨道设计过程中需要考虑的摄动力;然后分析了返回时近地距的各种影响因素,包括出发时刻、转移时间、倾角组合等,提出了一种可使返回近地距逐渐向设定的目标值靠拢的单向搜索方法,该方法简单易行,可信度高,能够满足轨道设计时要求的各项约束条件。该工作对中国将要进行的载人月球探测任务中需要用到的自由返回轨道的设计具有一定的参考价值。     

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

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

月面上升的轨迹优化与参数分析

面向月球采样返回任务分析需求,对月面上升段的轨迹优化及燃料消耗影响因素进行了研究。基于上升器运动模型,建立以燃料消耗最优为目标考虑入轨约束的轨迹优化模型,通过Gauss伪谱法和序列二次规划求解上升过程最优推力方向。改变运动模型中的初始推重比、入轨约束中的目标轨道参数,根据轨迹优化结果得到对应的燃料消耗,分析了这些因素对上升器燃料消耗的影响。针对上升器非共面起飞的问题,提出了上升偏航、升交点调整、倾角调整3种方案,从燃料消耗的角度分析了各方案的适用情况,为未来工程应用提供参考。     

部署retro-GEO巡视器的月球借力飞行轨道分析

retro-GEO是指逆行（retrograde）地球静止轨道（geostationary Earth orbit, GEO），该轨道与GEO轨道高度相同或相近，但倾角为180°，安装在retro-GEO卫星上的巡视器可每12h对GEO带空间资产附近碎片环境安全预警。直接西向发射retro-GEO卫星存在地面测控和发射能耗较大的困难。基于平面四体模型，为降低设计变量敏感性，以近月点参数为设计变量，建立了部署retro-GEO巡视器的月球借力飞行轨道设计模型，利用轨道动力学模型延拓策略，得到该类轨道绕月后返回地球飞行时长只能约为114.79h，该结论可用于求解该类轨道高精度轨道动力学模型解。     

气动辅助异面变轨在多碎片清除中的应用分析

碎片清除飞行器异面变轨需要消耗大量燃料.从气动辅助异面变轨优化设计及被清除碎片轨道高度差值、倾角差值等参数对变轨性能的影响出发,比较分析了优化气动辅助异面变轨与双脉冲霍曼轨道转移的燃料节约量,研究了不同轨道高度差对于实施气动辅助变轨燃料节约量的影响.当地球静止轨道（GEO）与低地轨道（LEO）间气动辅助变轨优化速度增量约为1.55km·s-1、质量面积比172kg·m-2、比冲310s、轨道倾角变化16°时,燃料节约率约为45%.对比研究了不同轨道高度差LEO轨道间实施气动辅助变轨的燃料节约情况.结果表明:随着轨道高度的增加,气动辅助优化效率逐渐降低;在相同高度轨道间实施异面变轨,随着轨道倾角的增加,气动辅助变轨燃料节约率先增大后减小,倾角改变量约为20°时,燃料节约率最大;当轨道倾角为5°时,采用气动辅助变轨和双脉冲变轨的燃料消耗量相同.      

LEO-LEO掩星事件持续时间及卫星间相对角速度仿真分析

通过对LEO-LEO掩星事件仿真,分析卫星轨道参数对掩星事件持续时间和卫星间相对角速度的影响.结果表明:对掩星事件持续时间影响较大的LEO-LEO卫星轨道参数依次为轨道倾角、升交点赤经、轨道高度,近地点角距对掩星事件持续时间影响不大;LEO-LEO卫星轨道倾角互补时,掩星事件持续时间最短,为98s;两颗LEO卫星轨道倾角之和或者升交点赤经之差在120°-240°范围内,掩星持续时间约100~150s;两颗LEO卫星轨道越高,掩星持续时间越短;掩星发生时,LEO卫星间水平方向相对角速度最大值约为0.14（°）·-1,俯仰方向相对角速度最大值约为0.078（°）·s-1.      

横向连续推力小偏心率人工冻结轨道设计

非临界倾角自然冻结轨道要求近地点幅角为90°, 偏心率在1× 10-3量级, 并且偏心率与轨道倾角和半长轴有特定的对应关系, 苛刻的要求极大地限制了轨道根数选择的多样性. 研究表明, 施加常幅值、方向半轨道周期切换的横向连续小推力可以形成任意轨道根数的小偏心率人工冻结轨道, 放宽自然冻结轨道对轨道根数的严苛限制. 维持小偏心率人工冻结轨道的横向连续推力a_T与半长轴的平方根成正比, 与偏心率大小成反比; 在临界倾角i_c附近, 推力a_T与倾角偏置量Δi近似呈正比, 且长期维持的燃料消耗较小. 此外, 在进行仿真计算时还需考虑拟平均轨道根数与瞬时轨道根数的转换.      

Integrated design of solar thermal propulsion system for microsatellite with liquid ammonia as propellant

This paper presents a design of solar thermal propulsion (STP) system for microsatellite with liquid ammonia as propellant. The system was equipped with two concentrators, which were respectively placed in the tank and thrust chamber for propellant supply and heating. A platelet heat exchanger was adopted to heat the propellant in the chamber, and the fluid–solid coupling effect between the wall and the gas was considered. Meanwhile, the effects of satellite mass, initial orbit, nozzle size and target temperature on the performance of STP system were analyzed. The results show that for microsatellites with a total mass of 100 kg, the STP system can fully heat the propellant to more than 2050 K, generate an intermittent thrust of about 26 N, and enable the satellite to obtain a velocity increment of more than 1470 m/s within 19 days, consuming only 42 kg of propellant, which can directly meet the transfer mission from the geostationary transfer orbit (GTO) to the geostationary orbit (GEO). The maximum velocity increment could reach more than 1950 m/s when the propellant was completely consumed; Changing the mass and initial orbit of the satellite will not affect the thrust and specific impulse. Satellites with smaller mass will spend less time and propellant during orbit transfer. The lower is the perigee height of the initial orbit, the greater is the propellant consumption, while the shorter is the time of orbit transfer; The reduction of nozzle throat size and target temperature will lead to the increase of specific impulse and the decrease of orbital transfer time, but the reduction of thrust.