深空探测器多次引力辅助转移轨道全局搜索

介绍行星引力辅助近似模型,包括近行星点无动力和有动力两种方式;采用开普勒轨道拼接法与双脉冲Lambert算法,建立多次引力辅助转移轨道的参数优化模型;通过广度优先搜索算法对引力辅助行星序列进行穷举,对发射窗口和天体间转移时间进行离散化网格搜索,每一步网格搜索后均采用合理的定界剪枝方法,减少后续计算量。这种全局搜索方法不需要提前指定引力辅助行星序列,并可得到对应不同飞行时间以及不同引力辅助次数的搜索结果,获得若干多天体引力辅助转移轨道初步结果,为进一步局部优化设计奠定基础。文章给出几组全局搜索算例,得到从地球到火星、木星和土星的多次引力辅助转移轨道,验证了全局搜索方法的有效性。     

Design of transfer trajectories between resonant orbits in the Earth–Moon restricted problem

The application of dynamical systems techniques to mission design has demonstrated that employing invariant manifolds and resonant flybys enables previously unknown trajectory options and potentially reduces the ΔV$\mathrm{\Delta }V$ requirements. In this investigation, planar and three-dimensional resonant orbits are analyzed and cataloged in the Earth–Moon system and the associated invariant manifold structures are computed and visualized with the aid of higher-dimensional Poincaré maps. The relationship between the manifold trajectories associated with multiple resonant orbits is explored through the maps with the objective of constructing resonant transfer arcs. As a result, planar and three-dimensional homoclinic- and heteroclinic-type trajectories between unstable periodic resonant orbits are identified in the Earth–Moon system. To further illustrate the applicability of 2D and 3D resonant orbits in preliminary trajectory design, planar transfers to the vicinity of L₅ and an out-of-plane transfer to a 3D periodic orbit, one that tours the entire Earth–Moon system, are constructed. The design process exploits the invariant manifolds associated with orbits in resonance with the Moon as transfer mechanisms.     

Low-energy Earth–Moon transfers involving manifolds through isomorphic mapping

Analysis and design of low-energy transfers to the Moon has been a subject of great interest for decades. Exterior and interior transfers, based on the transit through the regions where the collinear libration points are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is concerned with a geometrical approach for low-energy Earth-to-Moon mission analysis, based on isomorphic mapping. The isomorphic mapping of trajectories allows a visual, intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Two types of Earth-to-Moon missions are considered. The first mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a periodic orbit around the Moon. The second mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a capture (non-periodic) orbit around the Moon. In both cases three velocity impulses are needed to perform the transfer: the first at an unknown initial point along the low Earth orbit, the second at injection on the stable manifold, the third at injection in the final (periodic or capture) orbit. The final goal is in finding the optimization parameters, which are represented by the locations, directions, and magnitudes of the velocity impulses such that the overall delta-v of the transfer is minimized. This work proves how isomorphic mapping (in two distinct forms) can be profitably employed to optimize such transfers, by determining in a geometrical fashion the desired optimization parameters that minimize the delta-v budget required to perform the transfer.     

Manifold dynamics in the Earth–Moon system via isomorphic mapping with application to spacecraft end-of-life strategies

Recently, manifold dynamics has assumed an increasing relevance for analysis and design of low-energy missions, both in the Earth–Moon system and in alternative multibody environments. With regard to lunar missions, exterior and interior transfers, based on the transit through the regions where the collinear libration points L₁ and L₂ are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is focused on the definition and use of a special isomorphic mapping for low-energy mission analysis. A convenient set of cylindrical coordinates is employed to describe the spacecraft dynamics (i.e. position and velocity), in the context of the circular restricted three-body problem, used to model the spacecraft motion in the Earth–Moon system. This isomorphic mapping of trajectories allows the identification and intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Heteroclinic connections, i.e. the trajectories that belong to both the stable and the unstable manifolds of two distinct periodic orbits, can be easily detected by means of this representation. This paper illustrates the use of isomorphic mapping for finding (a) periodic orbits, (b) heteroclinic connections between trajectories emanating from two Lyapunov orbits, the first at L₁, and the second at L₂, and (c) heteroclinic connections between trajectories emanating from the Lyapunov orbit at L₁ and from a particular unstable lunar orbit. Heteroclinic trajectories are asymptotic trajectories that travels at zero-propellant cost. In practical situations, a modest delta-v budget is required to perform transfers along the manifolds. This circumstance implies the possibility of performing complex missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining suitable, convenient end-of-life strategies for spacecraft orbiting the Earth. Seven distinct options are identified, and lead to placing the spacecraft into the final disposal orbit, which is either (a) a lunar capture orbit, (b) a lunar impact trajectory, (c) a stable lunar periodic orbit, or (d) an outer orbit, never approaching the Earth or the Moon. Two remarkable properties that relate the velocity variations with the spacecraft energy are employed for the purpose of identifying the optimal locations, magnitudes, and directions of the velocity impulses needed to perform the seven transfer trajectories. The overall performance of each end-of-life strategy is evaluated in terms of time of flight and propellant budget.     

Existence of nonimpact motions along a wire rope fixed to an extended spacecraft

The possibility of nonimpact tension of a cable after its weakening when a small load moves along the cable whose ends are fixed on a massive dumbbell-like spacecraft which is moving in a steady-state manner along a circular orbit is considered. The conditions of existence and classification of the trajectories of nonimpact motion, including periodic ones, are presented.     

地月循环轨道动力学建模与计算研究

根据地月循环轨道的概念,按照生成第二类周期轨道的弧段进行分类,并讨论了其共振性与对称性在轨道设计与应用中的作用。然后归纳了三种循环轨道的动力学建模与计算方法及其轨道延拓策略,最后总结了三种方法的利弊和应用轨道类型。对地月系统循环轨道的研究和分析,能够为我国未来载人登月工程提供一种新的思路与理论支持。     

Optimization of low-energy resonant hopping transfers between planetary moons

In response to the scientific interest in Jupiter's Galilean moons, NASA and ESA have plans to send orbiting missions to Europa and Ganymede, respectively. The inter-moon transfers of the Jovian system offer obvious advantages in terms of scientific return, but are also challenging to design and optimize due in part to the large, often chaotic, sensitivities associated with repeated close encounters of the planetary moons. The approach outlined in this paper confronts this shortcoming by exploiting the multi-body dynamics with a patched three-body model to enable multiple “resonant-hopping” gravity assists. Initial conditions of unstable resonant orbits are pre-computed and provide starting points for the elusive initial guess associated with the highly nonlinear optimization problem. The core of the optimization algorithm relies on a fast and robust multiple-shooting technique to provide better controllability and reduce the sensitivities associated with the close approach trajectories. The complexity of the optimization problem is also reduced with the help of the Tisserand–Poincaré (T–P) graph that provides a simple way to target trajectories in the patched three-body problem. Preliminary numerical results of inter-moon transfers in the Jovian system are presented. For example, using only 59 m/s and 158 days, a spacecraft can transfer between a close resonant orbit of Ganymede and a close resonant orbit of Europa.     

Earth–Moon libration point orbit stationkeeping: Theory,modeling, and operations

Collinear Earth–Moon libration points have emerged as locations with immediate applications. These libration point orbits are inherently unstable and must be maintained regularly which constrains operations and maneuver locations. Stationkeeping is challenging due to relatively short time scales for divergence, effects of large orbital eccentricity of the secondary body, and third-body perturbations. Using the Acceleration Reconnection and Turbulence and Electrodynamics of the Moon's Interaction with the Sun (ARTEMIS) mission orbit as a platform, the fundamental behavior of the trajectories is explored using Poincaré maps in the circular restricted three-body problem. Operational stationkeeping results obtained using the Optimal Continuation Strategy are presented and compared to orbit stability information generated from mode analysis based in dynamical systems theory.     

地月高能共振循环轨道的快速计算及最优选择

针对现有高能共振循环轨道计算方法存在计算量大、有可能改变轨道共振特性和不能构造共振比大于2.3的地月循环轨道等缺点,本文提出了一种地月圆型限制性三体问题下高能共振循环轨道的快速计算方法。首先根据轨道在月球附近的组成弧段对高能共振循环轨道进行分类;然后根据轨道类型构建二体开普勒椭圆轨道;再进一步计算圆型限制性三体问题下的地月高能共振循环轨道;最后根据能量、稳定性、时间周期、近地点高度和近月点高度对所计算出的地月高能共振循环轨道进行最优选择。仿真结果表明,本文所提出的方法简单有效,能够计算出共振比为5:2的地月高能共振循环轨道。     

飞月轨道引力捕获设计方法研究

利用太阳引力摄动与月球绕飞设计的地月转移轨道（飞月轨道），与霍曼转移相比，虽然飞行时间较长（约三、四个月），但可显著节省速度增量（可达１５０米／秒），对无人月球探测器尤为适合。应用平面圆型限制性四体问题动力学模型，选择从月球出发的初始条件。借助“地心距－时间曲线”，从平面圆型限制性四体问题转换为一般的限制性四体问题。通过典型模拟计算，分析负向积分（从月球轨道出发）初始轨道参数及太阳方位对月球探测器     

基于轨迹成型法星际小推力转移轨道快速设计

给出了一种基于轨迹成型法星际小推力转移轨道的快速设计方法.首先介绍了基于轨迹成型法的有关内容,并通过与Hohmann变轨相对比验证了该方法描述小推力变轨的可用性.鉴于这种方法数学模型收敛性强,计算速度快的特点,以地火转移轨道为例,提出了不同推力条件下发射窗口的搜索方法.最后,利用基于轨迹成型法和真实数值仿真之间误差小、线性度好的特点,通过几次简单的线性迭代设计出小推力转移轨道.     

从地球到日-火系平动点轨道的转移

平动点轨道特殊的空间位置及动力学特征,使其在深空探测中具有重要的应用。以日-火系平动点轨道(Lissajous与Halo轨道)任务为目标,结合平动点轨道的不变流形理论,研究了小推力转移问题。首先给出了圆型限制性三体动力学模型下平动点附近不变流形(稳定和不稳定流形)高阶分析解以及相应的计算实例。接着以流形分析解为基础,建立了初始小推力轨道优化模型,并利用改进的协作进化算法求解初始小推力轨道。最后将初始轨道离散,采用多点打靶法将最优控制问题转化为参数优化问题,并用序列二次规划方法(SQP)求解。仿真结果证明轨道设计方法的有效性。     

A low-thrust transfer between the Earth–Moon and Sun–Earth systems based on invariant manifolds

A low-energy, low-thrust transfer between two halo orbits associated with two coupled three-body systems is studied in this paper. The transfer is composed of a ballistic departure, a ballistic insertion and a powered phase using low-thrust propulsion to connect these two trajectories. The ballistic departure and insertion are computed by constructing the unstable and stable invariant manifolds of the corresponding halo orbits, and a complete low-energy transfer based on the patched invariant manifolds is optimized using the particle swarm optimization (PSO) algorithm on the criterion of smallest velocity discontinuity and limited position discontinuity (less than 1 km). Then, the result is expropriated as the boundary conditions for the subsequent low-thrust trajectory design. The fuel-optimal problem is formulated using the calculus of variations and Pontryagin's Maximum Principle in a complete four-body dynamical environment. Then, a typical bang–bang control is derived and solved using the indirect method combined with a homotopic technique. The contributions of the present work mainly consist of two points. Firstly, the global search method proposed in this paper is simply handled using the PSO algorithm, a number of feasible solutions in a fairly wide range can be delivered without a priori or perfect knowledge of the transfers. Secondly, the indirect optimization method is used in the low-thrust trajectory design and the derivations of the first-order necessary conditions are simplified with a modified controlled, restricted four-body model.     

滚转对旋成体的菱形花纹的影响

这篇文章给出了滚转对旋成体的菱形花纹的影响。给出了任意形状旋成体滚转时的菱形花纹歪扭轨迹的计算公式。给出了三种物面形状的物体的ｎ（ｘ）计算公式：（１）锥（２）母线由直线和圆弧组成的尖头旋成体（３）球头锥。给出了参数ｋ＝０．０２５－０．５滚转圆锥的菱形花纹轨迹的空间图。计算了４种来流条件下的滚转尖头旋成体的菱形花纹轨迹。计算了３组模型和不同来流条件下滚转球头锥的菱形花纹轨迹，并对计算结果进行了分析比较     

载人登月绕月自由返回轨道与窗口精确快速设计

针对载人登月任务背景及工程约束,提出一种轨道与窗口一体化设计方法。通过两次坐标转化,将自由返回轨道设计参数解耦为近月点独立变量。在双二体假设下,通过4段二体轨道拼接完成自由返回轨道初值快速搜索及匹配近地停泊轨道(LEO)面的月窗口,其结果作为下一步采用序列二次规划（SQP）迭代求解高精度动力学模型轨道参数的初值,在该条精确轨道近月点时刻90 min邻域内产生可以匹配LEO地月转移入轨相位的零窗口轨道。算例表明,该流程能够精确快速地完成具有复杂任务背景及苛刻工程要求的载人登月绕月自由返回轨道与窗口设计问题。     

结合行星借力飞行技术的小推力转移轨道初始设计

针对结合行星借力和小推力技术的行星际转移轨道设计问题,提出一种基于形状逼近策略的初始设计方法。采用改进的逆六次多项式策略计算小推力弧段,通过引入B平面参数和推进器开关点时间系数实现行星借力和推滑混合轨道的拼接,将初始设计问题转化为求解混合整数非线性规划问题。为降低规划模型求解难度,通过参数变换对模型进行简化处理,并采用具有全局大范围搜索能力的改进微分进化算法求解最优设计参数。数值结果表明：相比正弦指数曲线设计方法,本文方法可以有效对交会型转移轨道进行设计,并且可以提供更少燃料消耗的探测机会。
     

PVM环境下提高并行计算加速比的数值实验研究

通过采用区域分割技术和拼接网络的接法，在ＰＶＭ并行环境下，对三高超音速绕流场场进行了多机并行计算，得到了较加速比，通过数值实验研究，讨论了负载平衡、不同的网格规模对加速比的影响。     

圆转方塞式喷管型面设计和试验研究

提出了圆转方塞式喷管的内喷管和塞锥型面的设计方法，内喷管用圆弧和抛物线近似，塞锥型面用抛物线和三次曲线近似，设计了一单元圆转方塞式喷管试验发动机。并采用气氧作氧化剂，气氢作燃料，进行了点火热试研究。介绍了试验发动机的结构与设计参数，以及试验系统组成和点火方式，给出了试验发动机照片、试验结果照片、测量参数曲线和性能数据处理。试验结果表明，试验发动机具有较高的热试效率：在三个不同工作高度下，喷管推力系数效率在93％-98％之间，说明圆转方塞式喷管的型面设计和试验方法是可行的。     

Spatial approaches to moons from resonance relative to invariant manifolds

In this study, the final approach to a moon or other body from resonance is explored and compared to the invariant manifolds of unstable periodic orbits. It is shown that the stable manifolds of planar Lyapunov orbits can act as a guide for the periods or resonances that are required for the final approach in both the planar and spatial problems. Previously developed techniques for the planar problem are expanded for use with resonances and used for comparison with trajectories approaching a moon from these resonances. A similar technique is then used for exploring the relationship of invariant manifolds to approach trajectories in the spatial problem. It is shown that the invariant manifolds of unstable periodic orbits provide insight into the trajectory design, and they can be used as a guide to the more direct approach trajectories.     

星历模型下地月自由返回全飞行过程轨道设计

面向载人登月任务需要,针对星历模型下具备自由返回能力的地月转移轨道设计问题进行了研究。在三体模型下对地月三维自由返回轨道进行了求解,得到了地月空间内的自由返回轨道分布情况。在二体模型假设下对近月段的三脉冲变轨进行了求解,给出了变平面机动的计算方法。进一步提出了两轮逐次优化修正策略,分别以高度和再入走廊为主要约束,采用内点法和SQP算法在高精度星历模型下对自由返回轨道初值进行逐次优化修正。之后,采用SQP算法在星历模型下对近月三脉冲变轨进行优化修正,得到了星历模型下的自由返回+近月三脉冲变轨地月转移策略。仿真校验结果表明本文提出的方法能够在给定约束下有效求解星历模型下具备自由返回能力的地月转移轨道,为载人登月任务的转移轨道设计提供参考。