卫星轨道偏心率的变化特征及其对轨道寿命的影响

轨道偏心率的变化极其重要,它是制约各类(不同高度)空间飞行体轨道寿命的关键因素之一.对于地球低轨卫星,主要受大气耗散作用的影响,而对环月(或环火星)低轨卫星,主要受非球形引力位中奇次带谐项的影响,会出现变幅较大的长周期变化,从而导致近星点高度hp在一段时间内有明显的下降趋势.对大偏心率轨道和高轨道,第三体的引力作用也会使e出现变幅较大的长周期变化,近星点高度hp也会有明显下降的现象,这都会影响卫星的轨道寿命,但这一动力学机制与大气耗散机制和非球形引力机制都不相同.即对轨道偏心率的变化特征及其对轨道寿命的影响作一综述.     

制约卫星轨道寿命的另一种机制

近点共振会导致太阳系小天体(小行星,自然卫星以及大行星和月球的人造卫星)的轨道偏心率出现变幅较大的长周期变化,特别是以月球和大行星为中心天体的大倾角轨道(确切地说是倾角接近90°的极轨道)卫星,由于类似的原因,偏心率的增大而导致近星距rp=a(1-e)≤ae(ae是中心天体的赤道半径),使其落到中心天体上,结束轨道寿命,这与耗散机制大不相同,因此将对其作理论分析,并以计算实例加以证实.     

第三体摄动分析解的一种表达式

在太阳系中,大行星、小行星和卫星（包括自然卫星和人造卫星）等对应的运动问题,都可以处理成受摄二体问题,而摄动源又多为第三体,作为第三体的摄动天体,有的比运动天体离中心天体近,有的则相反,前者称为内摄内体,全者则称为外摄天体,对一个具体的运动天体,可以同时出现这两个摄动天体,但是,只要运动天体与摄动天体的轨道都建立在以中心天体（质心）为坐标原点的同一坐标系内,那么在一定条件下（即除运动天体与摄动天体     

卫星径向位置摄动计算中的几个问题

为了对卫星轨道径向位置误差进行分析，本文将给出由地球非球形引力位（包括潮汐形变）引起的卫星径向位置摄动表达式，它将同时包含完整的卫星轨道偏心率的0阶和1阶项，并给出径向位置摄动空间分布的一种简单计算方法，它可明显地节省计算机时．     

中心体自转对天体轨道要素变化的后牛顿效应

本文给出了在三种引力理论为中心自转对天体轨道要素变化产生的后牛顿摄动效应的研究结果。研究结果表明：六个轨道要素除长钾不受摄动影响外其它五个要素均有周期摄动，特别升交点经度和近星点经度还有长期摄动效应。最后将文中的理论结论同前人的工作做了比较还应用于行星自转对卫星轨道要素变化的摄动效应计算上。作者在文[1]中研究了天体轨道要素变化的后牛顿效应，但在该文中并没有考虑中心体自转的影响。本文研究了三种引力理论（Einstein,Brans-Dick和Nordtvedt)中的这方面效应，并给出理论和数值的研究结果。     

星径向位置摄动计算中的几个问题

为了对卫星轨道径向位置误差进行分析，本文将给出由地球北球形引力位引起的卫星径向摄动表达式。它将同时包含完整的卫星轨道道偏心率的０阶和１阶项，并给出径向位置摄动空间分布的一种简单计算方法，它可明显地节省计算机时。     

关于月球低轨卫星运动的两个问题

对月球低轨卫星的轨道寿命特征和冻结轨道晶状态作了详尽的理论分析,给出它们与轨道倾角之间的关系以及它们相互之间的某种联系,并考虑低轨卫星的主要摄动源,在完整力模型下作了相应的模拟计算,不仅证实了理论分析的正确性,而且为环月运行探测器的轨道设计提供了极有参考价值的数值结果．     

火星非球形引力位田谐项联合摄动分析解

火星非球形引力场模型与地球有明显差别,其非球形引力位中的田谐项系数基本都要比地球的相应值大一个量级,尤其是J_2,2项（赤道椭率项）的大小接近它的动力学扁率项J₂.对于低轨探测器,若要使轨道外推1 d弧段的精度达到500 m(相当于标准单位10^-4量级),在构造环火探测器的轨道分析解时,田谐项与J₂项以及田谐项与田谐项之间的联合摄动不容忽视.根据摄动量级分析和构造的摄动分析解证实,上述联合摄动对轨道沿迹方向的影响可超过10^-4,并给出了数值验证.结果表明,与地球低轨卫星不同,在类似的问题中,构造环火卫星摄动分析解时,必须考虑这些联合摄动项的影响.     

关于数值求解天体运动方程的几个问题

本文讨论三个问题：1.在采用各种非辛（Symplectic）的数值积分器积分天体运动方程时，截断误差将引起人为的能量耗散，这一问题是不能用简单地在相应的力模型中加进一个人为的阻力因子而得以解决的，被歪曲的能量（或数值轨道）必须在积分过程的每一步用能量关系来进行校正，此即能量控制方法．2.当摄动加速度涉及到坐标轴的旋转时，如何在各种积分器中采用能量控制方法．3.对于大偏心率轨道，用数值方法求解相应运动方程时，积分步长必须随运动天体与中心天体之间的距离变化而改变，显然，这对所有积分器都是不方便的，特别是多步积分器．本义给出了一种步长均匀化的处理，可以使上述大偏心率轨道积分问题按定步长计算．     

利用太阳光压的大偏心率伴飞卫星轨道控制

提出了利用太阳帆进行大偏心率伴飞卫星轨道控制的方法.伴飞卫星围绕其惯量主轴做角速度恒定的自转,其惯量主轴在惯性系内指向保持不变.对伴飞卫星的控制分为轨道面的控制和轨道面内控制两部分.在控制过程中,优先考虑轨道面内的控制,在轨道面内控制不能进行(或者因为几何原因不能进行轨道面内控制)时,进行轨道面的控制.通过滑膜控制方法(Sliding Mode Control)计算轨道面内控制需要的控制力的方向和大小.得到需求的控制力要求后,推算出在控制过程中太阳帆相对于伴飞卫星主体的角度解析表达式.通过控制太阳帆的方向得到所需的不同的控制力.整个控制过程只针对伴飞卫星,主星处于自然飞行状态.最后对于这种控制方法进行数值验证.在无摄运动状态下通过控制系统进行伴飞轨道的轨道调整和误差消除,在考虑4阶非球形引力和第三体引力摄动情况下进行伴飞轨道的轨道维持.数值结果表明通过这种控制方法伴飞轨道能够保持轨道误差小于5 m.     

Two problems about the motion of low-moon-orbit satellites

A detailed theoretical analysis on the orbital lifetime and frozen orbit of low-moon-orbit satellites (LMOS) is carried out, and their relationships with the orbital inclination, as well as some mutual relationships are presented. Taking account of the main perturbing sources of low-orbit satellites, we carried out numerical simulations under a comprehensive force model, and the results not only confirm the correctness of the theoretical analysis, but also provide some valuable insights on the orbital design of LMOS.     

Orbital evolution of a planet with tidal dissipation in a restricted three-body system

The angle between planetary spin and the normal direction of an orbital plane is supposed to reveal a range of information about the associated planetary formation and evolution. Since the orbit's eccentricity and inclination oscillate periodically in a hierarchical triple body and tidal friction makes the spin parallel to the normal orientation of the orbital plane with a short timescale in an isolated binary system, we focus on the comprehensive effect of third body perturbation and tidal mechanism on the angle. Firstly, we extend the Hut tidal model(1981) to the general spatial case, adopting the equilibrium tide and weak friction hypothesis with constant delay time, which is suitable for arbitrary eccentricity and any angle ? between the planetary spin and normal orientation of the orbital plane. Furthermore, under the constraint of angular momentum conservation, the equations of orbital and ratational motion are given. Secondly, considering the coupled effects of tidal dissipation and third body perturbation, and adopting the quadrupole approximation as the third body perturbation effect, a comprehensive model is established by this work. Finally, we find that the ultimate evolution depends on the timescales of the third body and tidal friction. When the timescale of the third body is much shorter than that of tidal friction, the angle ? will oscillate for a long time,even over the whole evolution; when the timescale of the third body is observably larger than that of the tidal friction, the system may enter stable states, with the angle ? decaying to zero ultimately, and some cases may have a stable inclination beyond the critical value of Lidov-Kozai resonance. In addition, these dynamical evolutions depend on the initial values of the orbital elements and may aid in understanding the characteristics of the orbits of exoplanets.     

Averaged model to study long-term dynamics of a probe about Mercury

This paper provides a method for finding initial conditions of frozen orbits for a probe around Mercury. Frozen orbits are those whose orbital elements remain constant on average. Thus, at the same point in each orbit, the satellite always passes at the same altitude. This is very interesting for scientific missions that require close inspection of any celestial body. The orbital dynamics of an artificial satellite about Mercury is governed by the potential attraction of the main body. Besides the Keplerian attraction, we consider the inhomogeneities of the potential of the central body. We include secondary terms of Mercury gravity field from \(J_2\) up to \(J_6\), and the tesseral harmonics \(\overline{C}_{22}\) that is of the same magnitude than zonal \(J_2\). In the case of science missions about Mercury, it is also important to consider third-body perturbation (Sun). Circular restricted three body problem can not be applied to Mercury–Sun system due to its non-negligible orbital eccentricity. Besides the harmonics coefficients of Mercury’s gravitational potential, and the Sun gravitational perturbation, our average model also includes Solar acceleration pressure. This simplified model captures the majority of the dynamics of low and high orbits about Mercury. In order to capture the dominant characteristics of the dynamics, short-period terms of the system are removed applying a double-averaging technique. This algorithm is a two-fold process which firstly averages over the period of the satellite, and secondly averages with respect to the period of the third body. This simplified Hamiltonian model is introduced in the Lagrange Planetary equations. Thus, frozen orbits are characterized by a surface depending on three variables: the orbital semimajor axis, eccentricity and inclination. We find frozen orbits for an average altitude of 400 and 1000 km, which are the predicted values for the BepiColombo mission. Finally, the paper delves into the orbital stability of frozen orbits and the temporal evolution of the eccentricity of these orbits.     

Coordinate Additional Perturbations to Mars Orbiters and Choice of Corresponding Coordinate System

Similar to the study of the related problems of Earth satellites, in the research of the motion of Mars orbiter especially for low-orbit satellites, it is more appropriate to choose an epoch Mars-centered and Mars-equator reference system, which indeed is called the Mars-centered celestial coordinate system. In this system, the xy-plane and the direction of the x-axis correspond to the mean equator and mean equinox. Similar to the precession and nutation of the Earth, the wiggling of instantaneous Mars equator causes the coordinate additional perturbations in this Mars coordinate system. The paper quotes a method which is similar to the one used in dealing with the coordinate additional perturbations of Earth. According to this method, based on the IAU2000 Mars orientation model and under the precondition of a certain accuracy, we are able to figure out the precession part of the change of Mars gravitation. This lays the foundation for further study of its influence on the Mars orbiter's orbit of precession and the solution of the corresponding coordinate additional perturbations. The obtained analytical solution is easy to use. Compared with the numerical solution with higher accuracy, the result shows that the accuracy of this analytical solution could satisfy the general requirements in use. Therefore, our result verifies that a unified coordinate system, the Mars-centered celestial system in which J2000.0 is chosen as its current initial epoch, could be applied to deal with the relative problems of Mars orbiters, especially for low-orbit satellites. It is different from the method we previously used in dealing with the corresponding problems of Earth satellites, where we adopted the instantaneous equator and epoch (J1950.0) mean equinox as xy-plane and the direction of x -axis. In contrast, the coordinate transformation brings heavy workload and certain inconvenience in relative former works in which the prior system is used. If adopting the unified coordinate system, the transformation could be simply avoided and the computation load could be decreased significantly.     

Possible Effect of the Earth’s Inertial Induction on the Orbital Decay of LAGEOS

The theory of velocity dependent inertial induction, based upon extended Mach’s principle, has been able to generate many interesting results related to celestial mechanics and cosmological problems. Because of the extremely minute magnitude of the effect its presence can be detected through the motion of accurately observed bodies like Earth satellites. LAGEOS I and II are medium altitude satellites with nearly circular orbits. The motions of these satellites are accurately recorded and the past data of a few decades help to test many theories including the general theory of relativity. Therefore, it is hoped that the effect of the Earth’s inertial induction can have any detectable effect on the motion of these satellites. It is established that the semi-major axis of LAGEOS I is decreasing at the rate of 1.3 mm/d. As the atmospheric drag is negligible at that altitude, a proper explanation of the secular change has been wanting, and, therefore, this paper examines the effect of the Earth’s inertial induction effect on LAGEOS I. Past researches have established that Yarkovsky thermal drag, charged and neutral particle drag might be the possible mechanisms for this orbital decay. Inertial induction is found to generate a perturbing force that results in 0.33 mm/d decay of the semi major axis. Some other changes are also predicted and the phenomenon also helps to explain the observed changes in the orbits of a few other satellites. The results indicate the feasibility of the theory of inertial induction i.e. the dynamic gravitation phenomenon of the Earth on its satellites as a possible partial cause for orbital decay.     

Effects of the physical libration of the moon on lunar orbiters

Lunar physical libration, which is true oscillation of lunar equator in the space, alters the lunar gravitational field in the space coordinate system and affects the orbiting motion of lunar orbiters (hereafter called as lunar satellites) correspondingly. The effect is very similar to that of the precession and nutation on the earth satellites, and a similar treatment can be used. The variations in the gravitational force and in the orbit perturbation solution are clearly given in this paper together with numerical illustrations.     

Variations in the atmospheric neutral density at 145km

A least-squares multiple linear regression is performed on orbital decay density data obtained from precise orbital analysis of 22 low-perigee (130–160 km) Air Force satellites. Variations related to solar activity, the semi-annual effect, geomagnetic activity, and the zenith angle of the Sun are in agreement with the model of Jacchia (1971). Density variations in longitude and latitude are also deduced and compared with recent results from other investigations within this altitude regime.     

《大气一号》气球卫星轨道倾角变化分析

引起《大气一号》两颗气球卫星（ＤＱ－１Ａ和ＤＱ－１Ｂ）轨道倾角变化的摄动因素主要是太阳光压摄动、大气旋转和日月引力摄动。太阳光压摄动引起气球卫星轨道倾角增大，平均每天变化约０．００１７，大气旋转引起轨道倾角减小，平均每天变化不到０．０００１，但随着高度下降，变化量亦增大，陨落前达０．００２。本文根据卫星轨道摄动理论，给出气球卫星轨道倾角变化的一种定量分析方法，得到的分析结果为：（１）由太阳光压摄动