圆型限制性三体问题相对运动解析研究

针对圆型限制性三体问题共线平动点附近周期/拟周期轨道下的相对运动问题，提出一种新的、通用的解析研究方法。在周期/拟周期轨道近似解析解的基础上，结合微分修正方法，获得了精确的周期/拟周期轨道。对周期/拟周期轨道的单值矩阵进行分析，同时借鉴Floquet理论核心思想，建立了六个相对运动模态，并将相对运动表示为六个相对运动模态的线性组合，获得了相对运动的近似解析解。最后在地-月系统圆型限制性三体问题下，以L1点作为研究对象，分别以Halo轨道、Lissajous轨道和Lyapunov轨道为参考轨道，对相对运动模态和相对运动进行仿真分析，说明了相对运动模态的正确性以及相对运动近似解析解的有效性。     

Gravity field models from kinematic orbits of CHAMP,GRACE and GOCE satellites

The aim of our work is to generate Earth’s gravity field models from GPS positions of low Earth orbiters. Our inversion method is based on Newton’s second law, which relates the observed acceleration of the satellite with forces acting on it. The observed acceleration is obtained as numerical second derivative of kinematic positions. Observation equations are formulated using the gradient of the spherical harmonic expansion of the geopotential. Other forces are either modelled (lunisolar perturbations, tides) or provided by onboard measurements (nongravitational perturbations). From this linear regression model the geopotential harmonic coefficients are obtained.     

Space debris mitigation in geosynchronous orbit

In order to preserve the geosynchronous region, the Inter-Agency Space Debris Coordination Committee (IADC) proposed and endorsed a re-orbiting strategy for spacecraft at the end-of-life: they should be disposed above the synchronous altitude and passivated, to reduce the risk of inadvertent explosions. The recommended perigee altitude of the disposal orbit took into account all relevant perturbations and was a function of the expected perturbing acceleration induced by solar radiation pressure. It was intended to prevent any further interference with a properly defined geostationary protected region.     

The Investigation of Possible Approaches of Cataloged Space Objects to Manned Spacecraft

Initial orbital parameter errors are used to examine the miss distance between a spacecraft and an ensemble of tracked objects by a Monte Carlo-type analysis. The radial separation between orbits is evaluated and a keep-out zone is determined, which reduces the risk of collision to an acceptable level.An operational prediction methodology is suggested based on a catalog database, which identifies potentially hazardous approaches and computes the probability of collision for selected spacecraft. An example for the Mir Space Station is presented, which estimates the collision probability and the cross-sectional flux of cataloged objects for the time frame of interest. The results appear to be in good agreement with those of other space debris models.     

Resonant quasi-periodic near-rectilinear Halo orbits in the Elliptic-Circular Earth-Moon-Sun Problem

Motivated by the near-future re-exploration of the cislunar space, this paper investigates dynamical substitutes of the Earth-Moon’s resonant Near-Rectilinear Halo Orbits (NRHOs) under the Elliptic-Circular Restricted Four-Body Problem formulation of the Earth-Moon-Sun system. This model considers that the Earth and Moon move in elliptical orbits about each other and that a third body, the Sun, moves in a circular orbit about the Earth-Moon barycenter. By making use of this higher-fidelity dynamical model, we are able to incorporate the Sun’s influence and the Moon’s eccentricity, two of the most significant perturbations of the cislunar environment. As a result of these perturbations, resonant periodic NRHOs of the Earth-Moon Circular Restricted Three-Body Problem (CR3BP) are hereby replaced by two-dimensional quasi-periodic tori that better represent the dynamical evolution of satellites near the vicinity of the Moon. We present the steps and algorithms needed to compute these dynamical structures in the Elliptic-Circular model and subsequently assess their utility for spacecraft missions. We focus on the planned orbit for the NASA-led Lunar Gateway mission, a 9:2 synodic resonant L₂ southern NRHO, as well as on the 4:1 synodic and 4:1 sidereal resonances, due to the proximity to the nominal orbit and their advantageous dynamical properties. We verify that the dynamical equivalents of these orbits preserve key dynamical attributes such as eclipse avoidance and near-linear stability. Furthermore, we find that the higher dimensionality of quasi-periodic solutions offers interesting alternatives to mission designers in terms of phasing maneuvers and low-altitude scientific observations.     

地-月系平动点及Halo轨道的应用研究

地-月系统的平动点L1点及L2点的Halo轨道在探月工程中有重要的应用价值，可分别用于地月连续通信覆盖和月球背面的探测。由于在地-月系统中太阳的引力不可忽略，特别是在长时间作用以后，其动力学行为与摄动力较小的日-地系统有明显的不同。本文分析了如何利用太阳引力进入地-月系统的L1点及L2点的Halo轨道、以及由Halo轨道进入近月轨道的问题，两者综合起来构成了一条完整的地月低能转移轨道。研究结果对探月轨道设计有一定的参考价值。     

Effects of eccentricity of the reference orbit on multi-tethered satellite formations

The paper discusses the relevance of eccentric reference orbits on the dynamics of a tethered formation, when a massive cable model is included in the analysis of a multi-tethered satellite formation. The formations examined in this study are hub-and-spoke (HAS) and closed-hub-and-spoke (CHAS) configurations for in-plane and Earth-facing spin planes. Stability of the formations is studied by means of numerical simulation, together with the evaluation of the effects of eccentricity on tether elongation, agents relative position, and formation orientation and shape.     

Dynamical cartography of Earth satellite orbits

We have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term ($～$120?years) integrations of the non-averaged equations of motion. Hence, we characterize the long-term dynamics and the phase-space structure of the Earth-orbiter environment, starting from low altitudes ($～$400?km) and going up to the GEO region and beyond. This investigation was done in the framework of the EC-funded “ReDSHIFT” project, with the purpose of enabling the definition of passive debris removal strategies, based on the use of physical mechanisms inherent in the complex dynamics of the problem (i.e., resonances). Accordingly, the complicated interactions among resonances, generated by different perturbing forces (i.e., lunisolar gravity, solar radiation pressure, tesseral harmonics in the geopotential) are accurately depicted in our results, where we can identify the regions of phase space where the motion is regular and long-term stable and regions for which eccentricity growth and even instability due to chaotic behavior can emerge. The results are presented in an “atlas” of dynamical stability maps for different orbital zones, with a particular focus on the (drag-free) range of semimajor axes, where the perturbing effects of the Earth’s oblateness and lunisolar gravity are of comparable order. In some regions, the overlapping of the predominant lunisolar secular and semi-secular resonances furnish a number of interesting disposal hatches at moderate to low eccentricity orbits. All computations were repeated for an increased area-to-mass ratio, simulating the case of a satellite equipped with an on-board, area-augmenting device. We find that this would generally promote the deorbiting process, particularly at the transition region between LEO and MEO. Although direct reentry from very low eccentricities is very unlikely in most cases of interest, we find that a modest “delta-v” ($\mathrm{\Delta }V$) budget would be enough for satellites to be steered into a relatively short-lived resonance and achieve reentry into the Earth’s atmosphere within reasonable timescales ($～$50?years).     

基于非开普勒轨道的高超声速临近空间飞行器自主天文导航研究

针对航天器自主导航方法不适合高超声速临近空间飞行器的问题, 研究了基于非开普勒轨道的高超声速临近空间飞行器自主天文导航方案. 论述了基于非开普勒轨道的自主天文导航机理, 通过对高超声速临近空间飞行器受力分析, 建立了动力学方程; 利用矢量倒数法则推导出空间运动方程; 设计了基于非开普勒轨道的状态模型和基于星光折射间接敏感地平的观测模型, 采用卡尔曼滤波进行了仿真验证. 仿真结果表明, 基于非开普勒轨道的高超声速临近空间飞行器自主天文导航可达到较高的位置和速度精度.      

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.