姿态角幅值约束下的太阳帆Lissajous轨道保持控制

太阳帆航天器以两姿态角作为轨道控制输入时, 其轨道动力学方程具有非仿射非线性特性. 通过人工平动点处线性化获得的线性系统可完成太阳帆航天器轨道保持控制器的分析与设计. 由于线性近似模型为有误差模型, 存在近似有效范围约束, 表现为轨道高度约束和姿态角幅值约束. 本文研究了姿态角幅值约束对线性近似模型有效性的影响, 通过计算给出满足近似误差要求的姿态角幅值约束. 当控制输入存在幅值约束时, 控制器轨道修正能力受到束缚. 通过研究姿态角幅值约束下的最大允许入轨误差, 设计了最大允许入轨误差下线性二次型调节器(LQR)用于轨道保持控制, 并将控制器应用于太阳帆日地三体系统非线性模型中, 实现了日地人工L₁点Lissajous轨道最大允许入轨误差的控制收敛和良好精度下的轨道保持控制.      

Halo轨道族延拓方法及特性研究

对Halo轨道周期和运动范围等特性的研究是平动点任务设计的首要前提。针对大幅值Halo轨道和完整Halo轨道族的应用需求及其数值计算问题,面向当前应用广泛的地月系和日-地月系共线平动点,基于延拓法研究了圆型限制性三体问题下的Halo轨道族数值计算和运动学特性,给出了Halo轨道族延拓计算方法。数值仿真了族参数选择对轨道族计算的影响,得到了地月系和日-地月系共线平动点的大范围南北Halo轨道族,同时给出了轨道族的轨道周期变化和空间位置变化特性。研究结果表明,固定延拓步长下,L1点Halo轨道族应选择会合坐标系x坐标作为族参数,L2点Halo轨道族应选择y方向速度或者周期T作为族参数。方法适用于任意三体系统平动点的周期轨道族计算,特别是对其中的状态转移矩阵简单修改后可用于完整力模型下的Halo轨道（族）的数值设计。     

Integrated guidance and control for solar sail station-keeping with optical degradation

Current control approaches for solar sail station-keeping on libration point orbits have not considered the degradation of the sail’s optical properties. However, significant optical degradation could lead to poor station-keeping performance or even complete failure. This paper presents an integrated guidance and control strategy to address this problem by updating the reference orbit based on in situ estimation. An exponential optical degradation model is incorporated into the solar radiation acceleration model, and an on-line reference orbit update approach is incorporated into the station-keeping, coupled with an active disturbance rejection controller. The reflection coefficient is estimated on-line and the reference orbit is updated discretely when the optical properties have degraded by a prescribed amount. This strategy provides discrete updates to the reference orbits such that the perturbation due to the optical degradation is maintained within a small range. These smaller perturbations can be dealt with by the controller’s robustness and station-keeping can be sustained for long durations even in the presence of large optical degradation.     

Periodic attitude motions along planar orbits in the elliptic restricted three-body problem

A way to improve the accuracy of the three-body problem model is taking into account the eccentricity of primary attractors. Elliptic Restricted Three-Body Problem (ER3BP) is a model for studying spacecraft trajectory within the three-body problem such that the orbital eccentricity of primaries is reflected in it. As the principal cause of perturbation in the employed dynamical model, the primaries eccentricity changes the structure of orbits compared to the ideal Circular Restricted Three-Body Problem (CR3BP). It also changes the attitude behavior of a spacecraft revolving along periodic orbits in this regime. In this paper, the coupled orbit-attitude dynamics of a spacecraft in the ER3BP are exploited to find precise periodic solutions as the spacecraft is considered to be in planar orbits around Lagrangian points and Distant Retrograde Orbits (DRO). Periodic solutions are repetitious behaviors in which spacecraft whole dynamics are repeated periodically, these periodic behaviors are the main interest of this study because they are beneficial for future mission designs and allow delineation of the system’s governing dynamics. Previous studies laid the foundation for spacecraft stability analysis or studying pitch motion of spacecraft in the ER3BP regime. While in this paper, at first, initial guesses for correction algorithms were derived through verified search methods, then correction algorithms were used to refine calculated orbit-attitude periodic behaviors. Periodic orbits and full periodic solutions are portrayed and compared to previous studies and simpler models. Natural periodic solutions are valuable information eventuate in the longer functional lifetime of spacecraft. Since the problem assumption considered in this paper is much closer to real mission conditions, these results may be the means to use natural bounded motions in the actual operational environment.     

地月空间的远距离逆行轨道族及其分岔研究

地月系统中存在着一类绕月逆行、高度稳定的轨道族,称为远距离逆行轨道族(DRO)。以圆型限制性三体问题(CR3BP)为动力学模型研究了DRO轨道族周边的动力系统结构。利用Broucke稳定性图寻找分叉点,判断分叉类型,基于数值延拓计算分岔后产生的一系列新轨道分支。分叉类型主要有切分叉与多倍周期分叉(从3倍周期开始),轨道维度包含平面轨道族与三维轨道族。计算新轨道族的特征,包括形状、周期、能量、稳定性、双曲流形结构等。探讨周期轨道的轨道周期与能量的关系,以几何化的方式展现分叉结构、多周期轨道的双曲流形结构等。该动力结构将为基于DRO轨道族的地月空间任务提供重要的理论支持。      

基于三体系统平动点轨道的环日全景任务轨道设计

在空间开展太阳观测是研究太阳活动周、太阳爆发、极端天气等事件起源的重要手段。环日全景探测计划是为实现从黄道面360°全方位观察太阳行星际空间而提出的。本文针对环日全景探测计划,构建了基于三体系统平动点低能量轨道的环日全景轨道部署方法。该方法以日–地L1/L2点Halo轨道幅值及Halo轨道离轨点为变量,以转移轨道飞行时间、入轨机动大小为评价指标,基于三体系统不变流形构建环日全景的转移轨道,并开展轨道优化设计。采用等高线图对设计变量及任务成本进行全局分析。仿真计算表明,轨道部署无法同时满足飞行时间最短与入轨机动最小的要求。设计了轨道机动约束条件下的最优飞行时间解,并给出了基于长三甲运载火箭的一箭双星发射及入轨方案。      

On two kinds of intermediate orbits for asteroid explorations

In the paper, two kinds of intermediate orbits for asteroid explorations are proposed. One is around the collinear libration points of the Sun-asteroid restricted three-body problem. The other is around the asteroid itself. The first kind of intermediate orbit is applicable to asteroids with known masses, while the second is suitable for asteroids with unknown or negligible masses. Analytical solutions of these two intermediate orbits in the simplified models are introduced first, and then numerical algorithms are used to refine them to obtain the true orbits in the real force model. At last, the problem of station-keeping is addressed. The linear optimal feedback control law is used, and numerical simulations are made to both kinds of intermediate orbits. The results show that both kinds of orbits are feasible. The cost is reasonable and mainly depends on the initial insertion error.     

日地平动点卫星两脉冲转移轨道设计

研究基于最小二乘微分修正方法的平动点卫星两脉冲转移轨道设计,推导了考虑高度和航迹角约束的微分修正公式,讨论了该方法的收敛性.以日地L₁点附近的Halo轨道为目标轨道,在圆型限制性三体问题模型下设计了其转移轨道,系统地研究了HOI(Halo Orbit Insertion)点和Halo轨道幅值对转移轨道的影响,给出了HOI点的选择策略,并讨论了应急情况下快速转移轨道设计.数值仿真验证了方法的有效性,选择Halo轨道靠近地球侧的点作HOI点可以获得飞行时间适中的转移轨道.     

地月L2点周期轨道的月球背面覆盖分析

地月L2点附近轨道具备独特的动力学和运动学特性,是月球背面探测任务的中继卫星首选布设位置。面向未来月球背面探测任务的中继通信需求,分析并研究了地月L2点周期轨道(halo轨道)对月球背面的覆盖。在圆型限制性三体问题模型下,研究并给出了halo轨道族延拓计算方法,基于延拓法设计了地月系大范围南北halo轨道族;给出了中继卫星的月球背面覆盖计算模型,定义了相应的时间覆盖因子;数值仿真了地月系南北halo轨道族的月球背面覆盖情况。研究结果表明:地月L2点周期轨道幅值和类型决定其对月面的覆盖性,幅值较小的轨道的月面整体覆盖性较好,幅值较大的轨道对月球南北极覆盖较好,南北族轨道分别有利于月球南北半球的覆盖。文章研究可为我国"嫦娥4号"月球背面探测任务的中继星轨道设计提供有益参考和借鉴。     

Configuration space and stability analysis of solar sail near-vertical earth-trailing orbits

Regions outside the reach of traditional propulsion systems or the ones that require significant propellant, may be reached by harnessing the solar radiation pressure and leveraging coupled dynamics to maneuver a sail-based spacecraft. Earth-trailing orbits have recently been investigated for getting a unique perspective of the Sun while maintaining the spacecraft in close proximity to Earth. Vertical orbits trailing the Earth exhibit the additional capability to view the Sun from above and below the ecliptic plane. In this work, families of sail-based orbits are explored for varying Earth-trailing angles and Z amplitudes in the Sun-Earth circular restricted three-body problem. Optimization is carried out to ensure that the non-traditional vertical orbits exhibit a constant pitch angle control history, as well as symmetry across the X-Y plane. The stability of the resulting orbit families is assessed using an extension of Flouquet theory to Differential Algebraic Equations. Results indicate that sail-based Earth-trailing vertical orbits can be more stable than traditional sub-L₁ sail-based vertical orbits.     

椭圆三体问题中的时间周期不变流形

借助有限时间Lyapunov指数（FTLE）定义拉格朗日拟序结构（LCS）,并以单摆系统为例阐述LCS与动力系统中不变流形之间的联系.利用LCS研究椭圆限制性三体问题（ER3BP）中的时间周期不变流形的性质.采用数值方法验证得到了两点结论：时间周期不变流形的内部是穿越轨道集,外部是非穿越轨道集;时间周期不变流形是轨道的不变集.     

从日地系统L2出发借力月球飞越近地小行星

对于停留在日地系统L2的“嫦娥2号”探测器,其后续飞行方案有多个选项,例如主动撞月或重返月球轨道、返回地球轨道或再入大气、飞往地月系统L1/L2或日地系统L1、进入深空飞越近地小行星(最终,“嫦娥2号”于2012年12月13日成功地实现了对Toutatis小行星的近距离飞越)。探讨上述的飞行方案需要对飞行轨道进行初步设计,总的速度脉冲限制在100 m/s以内并且需要考虑探测器同时受到太阳、地球、月球的引力作用。本研究设计了探测器从日地系统L2出发借力月球实现Toutatis小行星飞越的飞行方案,与直接飞越方案相比,借力月球可以进一步节省探测器的燃料消耗,其等效速度脉冲设计值为58.47 m/s。     

地月平动点中继应用轨道维持

地月平动点中继应用轨道对于月球背面探测具有十分重要的应用价值,由于地月平动点的不稳定性,必须进行轨道维持。文章研究了真实力模型下月球平动点中继应用轨道的维持。首先,基于限制性三体问题下平动点轨道的运动特性,研究了平动点轨道维持的数学模型与维持策略,提出了平动点轨道维持的连续环绕控制方法,并给出了轨道维持的Halo和Lissajous两种控制方式;其次,充分考虑各天体和光压摄动下,采用数值手段研究了不同幅值的地月平动点周期中继应用轨道的维持间隔与速度增量等。研究结果表明:Lissajous控制方式适用于月球平动点中继应用轨道的维持,在给定测控精度条件下,维持间隔约7.4d,速度增量优于20m/s/a。该方法已经成功应用于我国"嫦娥2号"日地平动点任务和"嫦娥5T1"地月平动点任务并获得了良好的控制效果,还可直接应用于我国未来"嫦娥4号"等月球背面探测任务。     

Collinear equilibrium points in the elliptic R3BP with oblateness and radiation

We have studied the positions and stability of the collinear equilibrium points, L_1,2,3 of an infinitesimal body in the elliptic restricted three-body problem (ER3BP) when both primaries of the system are luminous and oblate spheroids moving in elliptic orbits around their common center of mass. We observe that their positions are affected by the radiation pressure forces, oblateness and the eccentricity of the orbits, but the stability character remains unchanged and are unstable. The effects of the parameters involved on the collinear points, in particular for the binary systems Achird, Luyten 726-8, Kruger 60, Alpha Centauri AB and Xi Bootis, and their stability in general have been investigated numerically using the analytical results obtained.     

Low-energy near Earth asteroid capture using Earth flybys and aerobraking

Since the Sun-Earth libration points L₁ and L₂ are regarded as ideal locations for space science missions and candidate gateways for future crewed interplanetary missions, capturing near-Earth asteroids (NEAs) around the Sun-Earth L₁/L₂ points has generated significant interest. Therefore, this paper proposes the concept of coupling together a flyby of the Earth and then capturing small NEAs onto Sun–Earth L₁/L₂ periodic orbits. In this capture strategy, the Sun-Earth circular restricted three-body problem (CRTBP) is used to calculate target Lypaunov orbits and their invariant manifolds. A periapsis map is then employed to determine the required perigee of the Earth flyby. Moreover, depending on the perigee distance of the flyby, Earth flybys with and without aerobraking are investigated to design a transfer trajectory capturing a small NEA from its initial orbit to the stable manifolds associated with Sun-Earth L₁/L₂ periodic orbits. Finally, a global optimization is carried out, based on a detailed design procedure for NEA capture using an Earth flyby. Results show that the NEA capture strategies using an Earth flyby with and without aerobraking both have the potential to be of lower cost in terms of energy requirements than a direct NEA capture strategy without the Earth flyby. Moreover, NEA capture with an Earth flyby also has the potential for a shorter flight time compared to the NEA capture strategy without the Earth flyby.     

Optimal station-keeping near Earth–Moon collinear libration points using continuous and impulsive maneuvers

In this study the gravitational perturbations of the Sun and other planets are modeled on the dynamics near the Earth–Moon Lagrange points and optimal continuous and discrete station-keeping maneuvers are found to maintain spacecraft about these points. The most critical perturbation effect near the L₁ and L₂ Lagrange points of the Earth–Moon is the ellipticity of the Moon’s orbit and the Sun’s gravity, respectively. These perturbations deviate the spacecraft from its nominal orbit and have been modeled through a restricted five-body problem (R5BP) formulation compatible with circular restricted three-body problem (CR3BP). The continuous control or impulsive maneuvers can compensate the deviation and keep the spacecraft on the closed orbit about the Lagrange point. The continuous control has been computed using linear quadratic regulator (LQR) and is compared with nonlinear programming (NP). The multiple shooting (MS) has been used for the computation of impulsive maneuvers to keep the trajectory closed and subsequently an optimized MS (OMS) method and multiple impulses optimization (MIO) method have been introduced, which minimize the summation of multiple impulses. In these two methods the spacecraft is allowed to deviate from the nominal orbit; however, the spacecraft trajectory should close itself. In this manner, some closed or nearly closed trajectories around the Earth–Moon Lagrange points are found that need almost zero station-keeping maneuver.     

基于电推进卫星飞往Halo轨道的转移轨道研究

利用电推进及轨道力学的特性实现节能优化,将限制性三体问题中的稳定不变流形与小推力轨道优化相结合,研究全电推进卫星从地球停泊轨道飞向日地拉格朗日L2点Halo轨道的低消耗转移轨道.航天器的转移轨道分为逃逸段、拼接段与无动力滑行段.在逃逸段卫星沿速度方向加速脱离地球引力,拼接段采用Radau伪谱法进行优化,使航天器以最短时间到达目标Halo轨道的稳定不变流形上,随后航天器电推进系统关机,沿稳定不变流形无动力滑行至目标轨道.基于雅克比积分常数给出拼接段轨道初始猜测值,以先提高切向方向航天器能量避免了全程优化离散点过多难以求解的问题.仿真结果表明,该方法收敛速度较快,对平动点工程任务的初期轨道特性计算具有实际意义.      

连续小推力非开普勒悬浮轨道研究综述

连续小推力非开普勒悬浮轨道在深空探测与地球极地观测任务中有着重要的应用前景。归纳了电推进、太阳帆推进等连续小推力技术的发展历程与现状；阐述了日心、行星悬浮轨道的动力学特性、稳定性、轨道保持策略；分析了三体问题下人工拉格朗日点的优势及其在深空探测方面的应用；讨论了悬浮轨道编队飞行的研究方法与控制策略。最后针对小推力悬浮轨道研究发展面临的难题，提出了研究新思路和应用新方向。     

Improving magnetosphere in situ observations using solar sails

Past and current magnetosphere missions employ conventional spacecraft formations for in situ observations of the geomagnetic tail. Conventional spacecraft flying in inertially fixed Keplerian orbits are only aligned with the geomagnetic tail once per year, since the geomagnetic tail is always aligned with the Earth-Sun line, and therefore, rotates annually. Solar sails are able to artificially create sun-synchronous orbits such that the orbit apse line remains aligned with the geomagnetic tail line throughout the entire year. This continuous presence in the geomagnetic tail can significantly increase the science phase for magnetosphere missions. In this paper, the problem of solar sail formation design is explored using nonlinear programming to design optimal two-craft, triangle, and tetrahedron solar sail formations, in terms of formation quality and formation stability. The designed formations are directly compared to the formations used in NASA’s Magnetospheric Multi-Scale mission.