卫星轨道摄动频谱分析

本文基于Ｌｉｅ级数的方法对卫星轨道摄动进行了频谱分析。这种摄动是由地球引力场的所有带谱调和项和田谐调和项引起的。为了适用小偏心率轨道，采用了Ｈｉｌｌ变量来描述卫星轨道，并在所获得的谱分析式中保留了所有与偏心率ｅ无关的项。在这个基础上又经过简单的变换给出了卫星位置摄动三分量的谱分析式。     

天问一号环绕器中继轨道受摄共振特性分析

“天问一号”环绕器在环火中继轨道为“祝融”火星车提供中继服务时,星下点出现了过着陆点纬度时的经度先增大再减小的特殊漂移现象,需要对该现象产生的原因进行深入分析。基于摄动理论设计了一系列数值实验,提出了一种星下点过指定纬度时经度的简化计算方法;并结合对田谐项引力场位函数的解析推导发现,由于中继轨道周期与火星自转周期比接近1:3,导致3阶3次田谐项产生了轨道共振,从而使得半长轴长周期运动有着较为显著的振荡,进而与带谐项摄动联合导致了该特殊现象。研究揭示了“天问一号”环绕器的轨道共振现象,发现了高阶田谐项对非同步轨道上的环火卫星仍能产生较为明显的轨道共振现象。研究结果可为后续火星环绕探测任务轨道设计提供重要参考。     

Artificial satellites orbits in 2:1 resonance: GPS constellation

In this work, the resonance problem in the artificial satellites motion is studied. The development of the geopotential includes the zonal harmonics J₂₀ and J₄₀ and the tesseral harmonics J₂₂ and J₄₂. Through an averaging procedure and successive Mathieu transformations, the order of dynamical system is reduced and the final system is solved by numerical integration. In the simplified dynamical model, three critical angles are studied. The half-width of the separatrix is calculated through a linearized model which describes the behavior of the dynamical system in a neighborhood of each critical angle. Through the resonance overlap criterion the possible regular and irregular motions are investigated by the time behavior of the semi-major axis, argument of perigee and eccentricity. The largest Lyapunov exponent is used as tool to verify the chaotic motion.     

Artificial satellites orbits in 2:1 resonance: GPS constellation

In this work, the resonance problem in the artificial satellites motion is studied. The development of the geopotential includes the zonal harmonics J₂₀ and J₄₀ and the tesseral harmonics J₂₂ and J₄₂. Through an averaging procedure and successive Mathieu transformations, the order of dynamical system is reduced and the final system is solved by numerical integration. In the simplified dynamical model, three critical angles are studied. The half-width of the separatrix is calculated through a linearized model which describes the behavior of the dynamical system in a neighborhood of each critical angle. Through the resonance overlap criterion the possible regular and irregular motions are investigated by the time behavior of the semi-major axis, argument of perigee and eccentricity. The largest Lyapunov exponent is used as tool to verify the chaotic motion.     

Collisions of stars by oscillating orbits of the second kind

Ordinary estimations of the number of star collisions in our galaxy—by simple kinematic considerations—lead to a very small number of such collisions: about one or even less every millions of years. However star collisions can occur through the following indirect way which has a much higher probability. (a) Binary stars are very common in our galaxy, about 30–50% of the stars. (b) If two binary stars meet a triple system can be formed by an ordinary exchange type motion. (c) A triple system is generally decomposed into the “inner orbit” (i.e. the relative orbit of the two nearest stars) and the “outer orbit” (i.e. the relative orbit of the third star with respect to the center of mass of the two nearest stars). The major axes of these two orbits have generally small perturbations and it is the same for the eccentricity of the outer orbit. On the contrary, if the relative inclination of the two orbits is large, the perturbations of the eccentricity of the inner orbit are important and can even in some cases lead to an eccentricity equal to one, that is to a collision of the two stars of the inner orbit.Such orbits can be called “oscillating orbits of the second kind”, indeed the first oscillating orbits—conceived by Khilmi and described for the first time in an example by Sitnikov—have unbounded mutual distances r_ij, but the system always come back to small sizes, it has an infinite number of very large expansions followed by strong contractions and, in the three-body case, an upper bound of lim inf (r_1.2 + r_1.3 + r_2.3) can be given in terms of the three masses and the integrals of motion. For the oscillating orbits of the second kind the mutual distances r_ij are bounded, but the velocities are unbounded (i.e. lim inf r_ij = 0 for at least one r_ij) and the system goes to a collision if the bodies have non-zero radius even small. The analytical study of the oscillating orbits of the second kind is a part of the general analytical study of the three-body problem, a part which must be valid for large eccentricities and large inclinations. The use of Delaunay's variables and of a Von Zeipel transformation lead to a first order integrable approximation, valid for any eccentricities and any inclinations, and giving the following results: (a) The oscillating orbits of the second kind occur when the angular momentum of the outer orbit has a modulus sufficiently close to the modulus of the total angular momentum of the three-body system. Hence these orbits occur for inclinations in the vicinity of 90°. (b) The oscillating orbits represent a set of positive measure of phase space and the first order study allows to give a rough estimation of the probability of collisions—even for stars of infinitely small radius. This probability, for given initial major axes and eccentricities and for isotropic arbitrary initial orientations, is generally of the order of m₃$\text{R}\text{M}$ (m₃ being the mass of the outer star, M the total mass and R the ratio of the period of the inner orbit to that of outer orbit).One question remains to be solved: how many collisions of stars are due to that phenomenon? That question is difficult because the probability of formation of a triple system by a random meeting of two binaries is very uneasy to estimate. However it seems that, compared to the usual evaluations based on pure kinematic considerations without gravitational effects, the number of collisions must be multiplied by a factor between one thousand and one million.     

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.     

Mapping orbits with low station keeping costs for constellations of satellites based on the integral over the time of the perturbing forces

The present paper is concerned with the search for orbits that have potential to require low fuel consumption for station-keeping maneuvers for constellations of satellites. The method used to study this problem is based on the integral over the time of the undesired perturbing forces. This integral measures the change of velocity caused by the perturbation forces acting on the satellite, so mapping orbits that are less perturbed, which generates good candidates for orbits that requires low fuel consumption for station-keeping maneuvers. The integral over the time depends only on the orbit of the spacecraft and the dynamical system considered. The type of engine and the control technique applied to the spacecraft are not considered to search for those orbits. It can be a good strategy to be applied for a first mapping of orbits. For this search, it is analyzed the integral of orbits with different values of the Keplerian elements in order to find the best ones with respect to this criterion. The perturbations considered are the ones caused by the third body, which includes the Sun and the Moon, and the J₂ term of the geopotential. The results presented here show numerical simulations to obtain the integral of those perturbing forces for different orbits. The GPS and the Molniya constellations are used as examples for those calculations.     

Applications of multi-body dynamical environments: The ARTEMIS transfer trajectory design

The application of forces in multi-body dynamical environments to permit the transfer of spacecraft from Earth orbit to Sun–Earth weak stability regions and then return to the Earth–Moon libration (L₁ and L₂) orbits has been successfully accomplished for the first time. This demonstrated that transfer is a positive step in the realization of a design process that can be used to transfer spacecraft with minimal Delta-V expenditures. Initialized using gravity assists to overcome fuel constraints; the ARTEMIS trajectory design has successfully placed two spacecrafts into Earth–Moon libration orbits by means of these applications.     

太阳帆航天器三维人工平动点变化特性研究

本文研究了太阳帆航天器在光压因子、锥角和钟角共同作用下的三维人工平动点的变化特性。在光压因子较小的时候，五个平动点会拓展成五个不相连的“人工”平动点面。随着光压因子增大，平动点面SL 3 ，SL 4 和SL5将逐渐扩大并互相融合，最终延展至SL 1 与之融合。而平动点面SL2则始终保持独立的球面，只随着光压因子的增大而扩大但不与其它平动点面发生融合。平动点位置的改变，意味着对应的周期轨道也随之改变，这为平动点周期轨道的转移等任务提供了有效参考。     

卫星热设计中β角在不同轨道下的变化规律分析

阳光矢量与卫星轨道平面之间夹角是卫星热控设计中至关重要的一个设计参数,其不仅与到达卫星表面的阳光热流密度密切相关,也与到达卫星表面的地球反射热流密度密切相关。文章根据β角的变化规律将传统意义上的轨道分为了同步轨道和非同步轨道两种,并结合升黄赤经夹角的变化规律分析了不同轨道下β角的变化规律,以及相应的热设计特点。     

Decay of high eccentricity satellite orbits with air drag using uniformly regular KS canonical elements

A non-singular analytical theory for the motion of high eccentricity satellite orbits under the influence of air drag is developed in terms of the Uniformly Regular Kustaanheimo and Stiefel (URKS) canonical elements, by assuming the atmosphere to be oblate with variation of density scale height with altitude. The series expansions include up to fourth-order terms of an independent variable Δ=λ² (function of eccentric anomaly) and c (a small parameter dependent on the flattening of the atmosphere). Only two of the nine equations are solved analytically due to symmetry in the equations of motion. Comparison of the important orbital parameters semi-major axis and eccentricity up to 1000 revolutions, obtained with the present analytical solution and the KS theory, shows the superiority of the present solution over the KS elements analytical solution. The theory can be used effectively for the orbital decay of aero-assisted orbital transfer orbits during mission planning.     

Control of satellite imaging formations in multi-body regimes

Libration point orbits may be ideal locations for satellite imaging formations. Therefore, control of these arrays in multi-body regimes is critical. A continuous feedback control algorithm is developed that maintains a formation of satellites in motion that is bounded relative to a halo orbit. This algorithm is derived based on the dynamic characteristics of the phase space near periodic orbits in the circular restricted three-body problem (CR3BP). By adjusting parameters of the control algorithm appropriately, satellites in the formation follow trajectories that are particularly advantageous to imaging arrays. Image reconstruction and coverage of the (u, v) plane are simulated for interferometric satellite configurations, demonstrating potential applications of the algorithm and the resulting motion.     

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.     

State propagation in an uncertain asteroid gravity field

Various spacecraft have been and will be sent to asteroids to characterize them. Generally, an asteroid's gravity field is very irregular and not accurately known when compared to the gravity field of a major planet, Earth in particular. It has been well studied that the irregularity significantly affects the trajectory of an orbiting spacecraft, and causes it to impact or to escape from the asteroid. Complementary to that, this paper focuses on the influence of the limited knowledge of this gravity field on the evolution of the spacecraft's orbit. It develops a general method by which this influence can be quantified. This method comprises specific Monte Carlo simulations with a discrete set of low-altitude orbits, taking into account the uncertainties in the gravity-field parameters. For illustration purposes, it is applied to two different asteroids. Already after three revolutions, the gravity-field uncertainties propagate to significant position uncertainties; this specifically holds for prograde orbits, and around the smaller asteroid. Applying this robust and accurate method helps mission designers and planners to assess the risk posed by gravity uncertainties, and take appropriate measures such as choosing the most favorable orbital geometries and/or lowering the orbit more slowly.     

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

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

On the features of long-term evolution of a high-apogee orbit of the SPECTR-R spacecraft

The practical tasks related to qualitative investigation of long-term evolution of high-apogee orbits of artificial Earth satellites (AES), for which the main perturbing factors are gravitational perturbations from the Moon and the Sun, are considered. Attention is given to the problem of the ballistic lifetime of similar orbits, and the issues associated with possibilities of the correction of orbits for ensuring the required duration of their ballistic lifetime are considered. The orbit of the SPECTR-R spacecraft launched in July of 2011 is considered as an example.     

关于旋进椭圆的讨论

根据Bertrand定理关于轨道闭合的条件,我们知道当中心引力的形式为径向距离的幂次型函数时,并不总能导致闭合轨道,即运动的质点不一定能返回到它自己的轨迹上去。这个定理指出:一切运动质点其初始条件稍微偏离圆轨道的要求条件时,只有当引力满足平方反比律(即牛顿万有引力定律)或线性定律(即虎克定律)时,轨道才是闭合的。 本文试图采用分析力学中作用角变量的方法,对一般情形下的轨道闭合条件进行讨论。这里我们看到使初始轨道接近于圆轨道的假设是不必要的,而且力的类型可推广到包括平方反比和其它幂次型的力。我们证明了轨道闭合的判别公式,并得到在上述情形下轨道非闭合的结论。 最后,如果中心引力包括有r~(-2)项和r~(-3)项时,象相对论性改正后所得到的那样,轨道就不再是闭合的,而是一旋进椭圆形轨道,椭圆的旋进角速度可以容易地计算出来。     

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.     

太阳帆绕地球周期轨道研究

  地球同步和太阳同步卫星在各个领域有着广泛的应用。静止轨道是一种特殊的地球同步轨道,轨道资源有限。利用化学推进或电推进可以实现轨道高度不同的同步轨道,如悬挂轨道,但需要消耗较多的燃料,工程上无法承受。本文考虑利用太阳帆实现地球同步和太阳同步轨道。太阳光压力在轨道平面内沿拱线方向,选择光压力与平面的夹角使得轨道平面的旋转速率与太阳光同步。研究表明,设计合适的半长轴和偏心率可以使得轨道旋转速率与地球自转速率一致。假设太阳光与赤道平面平行,可以得到准静止轨道,太阳帆将在传统静止轨道的附近运动,星下点的经度将在一个固定值附近振动。实际上太阳光是与黄道面平行,黄道面与赤道面之间存在夹角。考虑黄赤交角的情况下,太阳帆将在一定纬度和经度范围内运动。适合于对某个区域进行长期观测任务。     

A Concept of Guidance for an Air-Launch Vehicle with Compensation of Initial Downrange and Launch Time Errors at Direct Insertion into a Rendezvous Point in the Orbit

A flying launcher (airplane carrier) can generate initial errors in position and time of launch. In order to compensate for these errors, one should have two control parameters in addition to those that provide for a spacecraft's injection into a preset orbit. We suggest the concept of controlling the trajectory of injection by choosing thrust values (within allowable regions of control) of second-stage engines or/and of a space booster of the Polyot carrier launcher. As an example, a rendezvous of the spacecraft at the end of its boost phase with the International Space Station (ISS) is considered. The methodology of the suggested approach can be extended to other mobile systems of launch to rendezvous orbits.