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.     

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

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

Quaternion regularization in celestial mechanics and astrodynamics and trajectory motion control. I

Regularization problems in celestial mechanics and astrodynamics are considered. The fundamental regular quaternion models of celestial mechanics and astrodynamics are presented. It is shown that the efficiency of analytical investigation and numerical solution of boundary problems of optimal trajectory motion control of spacecraft may be increased using quaternion astrodynamics models. The regularization problem of celestial mechanics and astrodynamics that implies eliminating the feature, which arises in the equations of the two-body problem in case of impact of the second body with the central body, is considered in the first section of the paper. The quaternion method for regularizing the equations of the perturbed spatial two-body problem suggested by the author is presented; the method is compared with Kustaanheimo-Stiefel (KS) regularization. Demonstrative geometric and kinematic interpretations of regularizing transformations are provided. Regular quaternion equations for the two-body problem, which generalize the regular Kustaanheimo-Stiefel equations, as well as regular equations in quaternion osculating elements and quaternion regular equations for perturbed central motion of a material point, are considered. The papers on quaternion regularization in celestial mechanics and astrodynamics are briefly analyzed.     

解析求解非线性相对运动的最优重构问题

椭圆轨道相对运动模型的线性化导致其在大尺度相对运动应用中精度不能满足需求。针对任意椭圆轨道上的大尺度航天器编队最优重构问题,提出一种基于椭圆轨道非线性相对运动模型的近似解析求解方法。首先通过变分法建立了非线性最优重构问题的数学模型;然后采用摄动法,以偏近点角为积分变量求得了不含特殊积分的解析开环最优控制,有效地避免了真近点角域下最优控制解所含有的特殊积分。仿真验证了所求最优控制的有效性和优越性,结果表明在相对运动尺度较大时,相比基于椭圆轨道线性化模型的最优控制,在燃耗保持相近的情况下,所求非线性控制有效地降低了重构误差。     

The Initial Stage of “Space Debris” Cloud Formation during Planetary Satellite Breakdown

The process of formation of a particles cloud (space debris) during breakdown of the Earth’s Artificial Satellite (EAS) is considered. A system of equations of motion is derived, and its analytical solution in the linear approximation for circular orbits is presented. The distinction between this solution and an exact numerical solution is indicated. The cases of satellite breakdown at the apogee and perigee on an elliptic orbit are also considered, and some characteristics of these processes are compared.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 4, 2005, pp. 269–273.Original Russian Text Copyright © 2005 by Gorelov, Zarubkin.     

非线性条件下编队卫星周期性相对运动条件

利用编队卫星机械能守恒原理，提出了非线性条件下求解编队卫星周期性相对运动条件的新方法，给出了非线性周期相对运动的初始条件。编队卫星相对距离较近时，利用非线性周期运动条件，可修正Hill方程的初始条件，抑制编队卫星的长期漂移。编队卫星相对距离较大，非线性因素不可忽略时，利用非线性周期运动条件，可找到不需消耗任何燃料的周期性相对运动轨道。最后的数值仿真结果验证了该方法的正确性。     

Evolution of Motion of a Binary Planet

A model of a binary planet, consisting of a material point of small mass and a deformable viscoelastic sphere, is suggested. The center of mass of the binary planet moves in the gravitational field of a central body in the plane, which contains planets forming the binary planet. A deformable spherical planet rotates around the axis orthogonal to the plane of planetary motion. Planet deformations are described by the linear theory of viscoelasticity. It is shown that with an appropriate approximation of the gravitational potential, there is a class of quasicircular orbits, when the eccentricities of an orbit of the center of mass of a binary planet and an orbit, describing mutual planet motion, are equal to zero. The further evolution of motion is investigated in this class of orbits with the use of the canonical Poincare–Andoyer variables. Corresponding averaged equations are found, and phase pictures are constructed; the stability of stationary solutions is investigated on the basis of equations in variations. For the Solar system planets with their satellites, forming binary planets, the application of the presented model allows us to conclude that satellites sooner or later will fall on the corresponding planets.     

访问多个特定相对位置的航天器轨道设计

在航天器轨道设计问题中,将惯性空间中经典的吉布斯三矢量定轨方法拓展到相对运动空间中,给出了一种相对运动条件下的三矢量定轨方法。针对已知轨道的目标航天器,以及二个或三个给定的空间相对位置,基于相对运动方程,提出了设计跟随航天器飞行轨道的数值方法。以轨道面共面或异面,以及目标航天器轨道形状为椭圆或圆,将问题分为四种情况进行约束条件和自由变量个数的分析讨论。对于自由变量个数多于约束方程的情况,额外给定周期重访约束,将各种情况下的特定相对位置访问问题转化为一至二维的非线性方程(组)求解问题。对一维方程求解采用分段黄金分割+割线法进行快速求解;对二维方程组通过网格法搜索迭代初值并通过牛顿迭代快速求解。进一步基于线性模型的解,采用微分修正方法求解了各情况下J2摄动模型下的结果。数值算例验证了提出方法的正确性及有效性。     

高阶带谐摄动下卫星相对运动精确动力学模型

在考虑J 2 项摄动的卫星相对运动精确动力学模型基础上,推导了可以精确到任意阶带谐非球形引力摄动的卫星相对运动精确动力学模型,还给出了一组简洁的基于参考卫星轨道要素（Reference Satellite Variables）的微分方程,能够描述精确到任意阶带谐非球形引力的卫星地球轨道运动。由于在模型推导过程中没有引入任何简化,新推导的结果适用于任意偏心率的地球轨道。     

Cubature Kalman filtering for relative spacecraft attitude and position estimation

A novel relative spacecraft attitude and position estimation approach based on cubature Kalman filter is derived. The integrated sensor suit comprises the gyro sensors on each spacecraft and a vision-based navigation system on the deputy spacecraft. In the traditional algorithm, an assumption that the chief?s body frame coincides with its Local Vertical Local Horizontal (LVLH) frame is made to construct the line-of-sight observations for convenience. To solve this problem, two relative quaternions that map the chief?s LVLH frame to the deputy and chief body frames are involved. The general relative equations of motion for eccentric orbits are used to describe the positional dynamics. The implementation equations for the cubature Kalman filter are derived. Simulation results indicate that the proposed filter provides more accurate estimates of relative attitude and position over than the extended Kalman filter.     

基于三阶解析解的小行星平衡点附近halo轨道确定方法研究

尽管周期解的存在性已经被证明,但要在给定的动力学系统中寻找到满足一定精度要求的周期解依然是一件极富挑战性的工作.提出如下方法确定小行星平衡点附近精确的周期轨道(halo轨道).首先扩展运动方程:将小行星平衡点附近轨道运动方程的右端项在平衡点处展成三阶幂级数.从而将非线性运动学方程扩展为拟线性微分方程.然后求近似解析解:应用Lindstedt-Poincaré方法求解扩展后的运动方程组,将周期解和其运动频率展开成三阶幂级数,并将二者代人扩展后的拟线性微分方程中.这样就可以得到三个不同阶的线性运动方程,逐次求解三个微分方程并消除解中的永年项即可得到hal.轨道的三阶解析解.最后微分校正:将周期轨道在三阶解析解附近线性化,得到状态转移矩阵,并使用状态转移矩阵和轨道终端状态的偏差修正轨道初值,从而得到满足精度要求的精确引力场中的halo轨道.     

The motion of a particle in perturbed field of an attracting center

The integrable case of the perturbed two-body problem is considered. The perturbation is determined by the potential of a special form. The L-matrix is chosen in such a way that partial separation of variables should take place in regular coordinates. Integration of the equations of motion of the problem under consideration is made. The solutions are expressed through elliptic functions. The orbits for various cases are constructed. The results of numerical calculations are given.     

椭圆轨道航天器自主接近的制导律研究

椭圆轨道自主交会接近制导律问题是当今空间技术领域中的一个研究不足。基于T\|H方程 ,设计了两种仅利用相对运动信息、控制力定常的椭圆轨道自主接近制导律：滑模制导律和 势函数制导律。从两种制导律的特性分析,滑模制导律对系统的运行轨迹更加苛刻,而势函 数制导律更加充分地利用了系统的自由运动,具有更好的控制性能。数值仿真结果也表明： 对 于相同的制导目标、制导时间和制导精度,势函数制导律所需的喷气开关次数和速度增量均较少。     

Biaxial Mode of Rotation of a Satellite in the Orbit Plane

A mode of motion of a satellite with respect to its center of mass is studied, which is called the biaxial rotation in the orbit plane. In this mode of rotation, an elongated and nearly dynamically symmetric satellite rotates around the longitudinal axis, which, in turn, rotates around the normal to the plane of an orbit; the angular velocity of rotation around the longitudinal axis is several times larger than the orbital angular velocity, deviations of this axis from the orbit plane are small. Such a rotation is convenient in the case when it is required to secure a sufficiently uniform illumination of the satellite's surface by the Sun at a comparatively small angular velocity of the satellite. The investigation consists of the numerical integration of equations of the satellite's motion, which take into account gravitational and restoring aerodynamic moments, as well as the evolution of the orbit. At high orbits, the mode of the biaxial rotation is conserved for an appreciable length of time, and at low orbits it is destroyed due to the impact of the aerodynamic moment. The orbit altitudes and the method of constructing the initial conditions of motion that guarantee a sufficiently prolonged period of existence of this mode are specified.     

Translational–Rotational Motion of a Viscoelastic Sphere in Gravitational Field of an Attracting Center and a Satellite

We study the translational–rotational motion of a planet modeled by a viscoelastic sphere in the gravitational fields of an immovable attracting center and a satellite modeled as material points. The satellite and the planet move with respect to their common center of mass that, in turn, moves with respect to the attracting center. The exact system of equations of motion of the considered mechanical system is deduced from the D'Alembert–Lagrange variational principle. The method of separation of motions is applied to the obtained system of equations and an approximate system of ordinary differential equations is deduced which describes the translational–rotational motion of the planet and its satellite, taking into account the perturbations caused by elasticity and dissipation. An analysis of the deformed state of the viscoelastic planet under the action of gravitational forces and forces of inertia is carried out. It is demonstrated that in the steady-state motion, when energy dissipation vanishes, the planet's center of mass and the satellite move along circular orbits with respect to the attracting center, being located on a single line with it. The viscoelastic planet in its steady-state motion is immovable in the orbital frame of reference. It is demonstrated that this steady-state motion is unstable.     

追踪星跟踪空间非合作目标的相对轨道设计

对空间非合作目标跟踪飞行可以执行观测或监视等任务。首先从一般性出发对追踪星与非合作目标之间的椭圆轨道相对运动方程进行分析,给出具有任意初始条件的相对运动方程解析表达式。其次,对追踪星沿航向跟踪目标并考虑约束条件时的相对轨道设计进行分析后,给出设计追踪星轨道的方法,该方法使得追踪星在保持对地定向的同时也满足测量敏感器的约束条件。最后通过数学仿真进行了验证。     

Application of the stroboscopic method to equinoctial elements

In this paper the stroboscopic method is applied to the equinoctial elements which avoid the singularities of circular and/or equatorial orbits. The Lagrange equations for the variation of parameters are formulated using respectively one of the three longitudes as fast angular variable. It is shown how the first-order theory of the stroboscopic method can be developed. The perturbation by the gravity potential of the central body and the third-body perturbation are considered in detail. The paper concludes with a few analytical results.     

邻近近圆轨道两个飞行器的最优交会

邻近近圆轨道上两个飞行器的相对运动可以用线性化Ｃ—Ｗ方程描述。本文利用线性化Ｃ—Ｗ方程研究了常推力作用下两个飞行器的固定时间最省燃料交会问题，给出了最优推力作用的开关函数和最优推力方向的计算方法。     

Formation design in elliptical orbit using relative orbit elements

A new set of relative orbit elements is strictly defined through spherical geometry. The exact transformation equations between the new relative orbit elements and classical-orbital elements are derived. A new relative motion model with no singularity problem is derived based on the relative orbit elements, which are suitable for both elliptical and circular reference orbits. The in-plane and out-of-plane relative motion can be completely decoupled based on the new model. The inverse transformation of state transfer matrix is obtained to analyze perturbation effects and control strategy. The geometric characteristics of relative motion can be easily described using the relative eccentricity/inclination vector method. The proposed method and conclusions are validated by simulation through some typical examples. This paper improves the basic theory of relative orbit elements and unifies the expressions of the elliptical and near-circular close relative motion.