月球卫星轨道力学综述

月球探测器的运动通常可分为3个阶段，这3个阶段分别对应3种不同类型的轨道：近地停泊轨道、向月飞行的过渡轨道与环月飞行的月球卫星轨道。近地停泊轨道实为一种地球卫星轨道；过渡轨道则涉及不同的过渡方式(大推力或小推力等)；环月飞行的月球卫星轨道则与地球卫星轨道有很多不同之处，它决不是地球卫星轨道的简单克隆。针对这一点，全面阐述月球卫星的轨道力学问题，特别是环月飞行中的一些热点问题，如轨道摄动解的构造、近月点高度的下降及其涉及的卫星轨道寿命、各种特殊卫星(如太阳同步卫星和冻结轨道卫星等)的轨道特征、月球卫星定轨等。     

考虑地球扁率摄动影响的初轨计算方法

在二体问题意义下的短弧定轨,Laplace型方法是最主要最典型的一种初轨计算方法。若测角资料达到10^-4-10^-5精度(相当于2″—20″之间),那么要使定轨精度达到与其相应的程度,地球非球形引力位中的扁率项摄动应该考虑,在此前提下,同样可以采用相应的Laplace型定轨方法。即给出这种严格包含扁率摄动的初轨计算方法的原理和具体计算过程以及计算实例,除采用多资料定轨方法外,这种方法也是提高初轨计算精度的一种途径,它同样可用于多资料的情况,这种方法对于大扁率主天体(即中心天体)的卫星定轨将更有实用价值。     

关于星座小卫星的编队飞行问题

从轨道力学角度来看星座小卫星编队飞行和星星跟踪中的伴飞，遵循着如下动力学机制：(1)在各小卫星绕地球运动过程中轨道摄动变化的主要特征决定了星-星之间的空间构形，(2)当星星之间相互距离较近时，在退化的限制性三体问题(实为限制性二体问题)中，共线秤动点附近的条件周期运动亦可在一定时间内制约星-星之间的空间构形．将具体阐明这两种动力学机制的原理和相应的星星之间的相对构形，并用仿真计算来证实这两种动力学机制的适用范围，为星座小卫星编队飞行和在伴飞运动过程中进行轨控提供理论依据和具体的轨控条件．     

月球物理天平动对环月轨道器运动的影响

月球物理天平动是月球赤道在空间真实的摆动，会导致月球引力场在空间坐标系中的变化，从而引起环月轨道器(以下称为月球卫星)的轨道变化，这与地球的岁差章动现象对地球卫星轨道的影响类似．采用类似对地球岁差章动的处理方法，讨论月球物理天平动对月球卫星轨道的影响，给出相应的引力位的变化及卫星轨道的摄动解，清楚地表明了月球卫星轨道的变化规律，并和数值解进行了比对，从定性和定量方面作一讨论．     

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

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

同步卫星短弧定轨

同步卫星受到摄动力的影响,它的实际轨道有一点漂移．卫星需要不断的调轨调姿,以保证其正常运行．为了研究卫星在几小时,甚至更短的时间内的轨迹情况,采用短弧段定轨法．用动力学方法进行短弧定轨,分别研究1小时和15分钟定轨并进行比较,目的是为了在同步轨道卫星变轨后,能尽快地为卫星提供精密的预报轨道．此外,在系列短弧定轨后,得到精密轨道系列,为研究轨道变化的力学因素及研究短弧中卫星转发器时延变化规律等提供依据．     

地-月系平动点轨道的特征及其相关问题

当探测器定点在地-月系共线平动点L_1、L_2附近的halo轨道或Lissajous轨道时,由于其固有的动力学特征,通常是被人们置于地-月系质心旋转坐标系中展现其几何特征.其实,它们同样是环绕地球运行的Kepler轨道,这类探测器实为地球的远地卫星.但由于其自身所具有的不稳定性特征,在轨道外推中,初值误差的传播程度远比一般的环绕型探测器轨道外推显著.这在轨道设计、运行控制和地面测控等领域都是需要重视的问题.尽管如此,除在构造这类轨道变化的受摄分析解时遇到困难外,对其定轨等问题,与一般远地卫星类似,并无其他特殊要求.将具体给出该类轨道由于不稳定特性引起误差快速传播的定量状态和相应的理论分析,以及实际应用中的短弧定轨和相应的高精度轨道预报方法,并附有实测资料进行定轨结果的检验.     

关于共线平动点的特征及其在深空探测中的应用

系统阐述了小天体运动对应的圆型限制性三体问题共线平动点的强不稳定性特征,以及其附近的条件周期轨道——晕轨道(Halo Orbit)的存在、相应解的构造。这种特殊的轨道形式和共线平动点附近的弱稳定走廊,可分别用于在深空特殊位置附近定点有各种科学探测目标的探测器和向节能轨道过渡的通道。     

一种解析定轨方法

本文给出了人造地球卫星轨道计算的一种解析方法,定轨方案中摄动计算考虑了地球引力场非球形摄动的J2,J3,J4的长期项,长周期项,J2短周期项,大气阻力,太阳光压及日月引力摄动的长期项。初始根数改正量估计采用微分轨道改进算法。在定轨迭代收敛后,残差的均方根误差在5″左右,资料使用率超过80%。     

地球非球形引力位中田谐项摄动的有关问题

在人造卫星绕地球运动中,地球非球形引力摄动是最重要的讨论了容易被忽视的田谐项摄动,尽管它对低轨卫星的影响,只相当于J2项的二阶量,又是短周期效应,但它却包含了大10多倍的地球自转项,必须给以重视.还导出包含全部阶次田谐项的摄动解,并分离出地球自转项,对轨道半长径a还增加了(J2×Jlm)联合摄动的地球自转项,既为理论分析提供依据又可用于分析法定轨和预报.     

满足一定约束条件的登月飞行轨道的设计

讨论满足约束条件的登月飞行轨道的设计问题，将约束条件分类为只与太阳，月球，地球，飞行器和观测站之间的相对位置有关的运动学约束条件以及小及到飞行器轨道云动的动力学约束条件，在考虑登月飞行轨道的特征后，给出一种设计满足约束条件的标准飞行轨道的方法，并将方法应用于不同约束条件下的我国登月飞行以及月球卫星的轨道预测计。     

Some orbital characteristics of lunar artificial satellites

In this paper we present an analytical theory with numerical simulations to study the orbital motion of lunar artificial satellites. We consider the problem of an artificial satellite perturbed by the non-uniform distribution of mass of the Moon and by a third-body in elliptical orbit (Earth is considered). Legendre polynomials are expanded in powers of the eccentricity up to the degree four and are used for the disturbing potential due to the third-body. We show a new approximated equation to compute the critical semi-major axis for the orbit of the satellite. Lie-Hori perturbation method up to the second-order is applied to eliminate the terms of short-period of the disturbing potential. Coupling terms are analyzed. Emphasis is given to the case of frozen orbits and critical inclination. Numerical simulations for hypothetical lunar artificial satellites are performed, considering that the perturbations are acting together or one at a time.     

Spacecraft Orbit Determination with The B-spline Approximation Method

It is known that the dynamical orbit determination is the most common way to get the precise orbits of spacecraft. However, it is hard to build up the precise dynamical model of spacecraft sometimes. In order to solve this problem, the technique of the orbit determination with the B-spline approximation method based on the theory of function approximation is presented in this article. In order to verify the effectiveness of this method, simulative orbit determinations in the cases of LEO (Low Earth Orbit), MEO (Medium Earth Orbit), and HEO (Highly Eccentric Orbit) satellites are performed, and it is shown that this method has a reliable accuracy and stable solution. The approach can be performed in both the conventional celestial coordinate system and the conventional terrestrial coordinate system. The spacecraft's position and velocity can be calculated directly with the B-spline approximation method, it needs not to integrate the dynamical equations, nor to calculate the state transfer matrix, thus the burden of calculations in the orbit determination is reduced substantially relative to the dynamical orbit determination method. The technique not only has a certain theoretical significance, but also can serve as a conventional algorithm in the spacecraft orbit determination.     

On analytic modeling of lunar perturbations of artificial satellites of the earth

Two different procedures for analytically modeling the effects of the Moon's direct gravitational force on artificial Earth satellites are discussed from theoretical and numerical viewpoints. One is developed using classical series expansions of inclination and eccentricity for both the satellite and the Moon, and the other employs a method of averaging. Both solutions are seen to have advantages, but it is shown that while the former can be more accurate in special situations, the latter is quicker and more practical for the general orbit determination problem where observed data is used to correct the orbit in near real time.This work was sponsored with the support of the Department of the Air Force under contract F19628-85-C-0002. The views expressed are those of the author and do not reflect the official policy or position of the US Government.     

基于双频星载GPS数据的LEO卫星运动学定轨研究

运动学定轨是星载GPS特有的定轨方法,该方法不依赖于任何力学模型(地球重力场、大气阻力及太阳辐射压等),尤其适用于受大气阻力影响严重的低轨卫星定轨.基于双频星载GPS数据,研究了运动学定轨原理,讨论了数据预处理方法,建立了一套非差运动学定轨算法.并以GRACE (Gravity Recovery And Climate Experiment)-A、B卫星2008年2月实测数据作为试算验证了本研究方法的有效性和可靠性.GRACE 卫星实测数据计算结果表明:运动学定轨能达到5 cm精度(相对于SLR (Satellite Laser Ranging)),与动力学和简化动力学定轨精度相当.     

Analytical methods for the orbits of artificial satellites of the Moon

The motion of a close artificial satellite of the Moon is considered. The principal perturbations taken into account are caused by the nonsphericity of the Moon and the attraction of the Earth and the Sun. To begin with, the expansions of the disturbing functions due to the nonsphericity of the primary body and the action of the disturbing mass-point body have been derived. The second expansion is produced in terms of the Keplerian elements of a satellite and the spherical coordinates of the disturbing body. Both expansions are valid for an arbitrary reference plane. The motion of a satellite of the Moon is studied in the selenocentric coordinate system referred to the Lunar equator and rotating with respect to the fixed ecliptic system. However, the coordinate exes in the equatorial plane are chosen so that the angular speed of rotation of the system is small. The motion of the satellite is described by means of the contact elements which enable one to utilize the conventional Lagrange's planetary equations and may be regarded as the generalization of the notion of the osculating elements to the case of the disturbing function depending not only o the coordinates and the time but on the velocities as well. Two methods are proposed to represent the motion of Lunar satellites over long intervals of time: the von Zeipel method and the Euler method of analytical integration with application of the variation-of-elements technique at every step of integration. The second method is exposed in great detail.Presented at the Meeting of Commission 7 of the IAU on Analytical Methods for the Orbits of Artificial Celestial Objects 14-th General Assembly of the IAU, Brighton, 1970.     

Approximate analytic method for high-apogee twelve-hour orbits of artificial Earth’s satellites

We propose an approach to the study of the evolution of high-apogee twelve-hour orbits of artificial Earth’s satellites. We describe parameters of the motion model used for the artificial Earth’s satellite such that the principal gravitational perturbations of the Moon and Sun, nonsphericity of the Earth, and perturbations from the light pressure force are approximately taken into account. To solve the system of averaged equations describing the evolution of the orbit parameters of an artificial satellite, we use both numeric and analytic methods. To select initial parameters of the twelve-hour orbit, we assume that the path of the satellite along the surface of the Earth is stable. Results obtained by the analytic method and by the numerical integration of the evolving system are compared. For intervals of several years, we obtain estimates of oscillation periods and amplitudes for orbital elements. To verify the results and estimate the precision of the method, we use the numerical integration of rigorous (not averaged) equations of motion of the artificial satellite: they take into account forces acting on the satellite substantially more completely and precisely. The described method can be applied not only to the investigation of orbit evolutions of artificial satellites of the Earth; it can be applied to the investigation of the orbit evolution for other planets of the Solar system provided that the corresponding research problem will arise in the future and the considered special class of resonance orbits of satellites will be used for that purpose.     

Lunar perturbations of artificial satellites of the earth

Lunisolar perturbations of an artificial satellite for general terms of the disturbing function were derived by Kaula (1962). However, his formulas use equatorial elements for the Moon and do not give a definite algorithm for computational procedures. As Kozai (1966, 1973) noted, both inclination and node of the Moon's orbit with respect to the equator of the Earth are not simple functions of time, while the same elements with respect to the ecliptic are well approximated by a constant and a linear function of time, respectively. In the present work, we obtain the disturbing function for the Lunar perturbations using ecliptic elements for the Moon and equatorial elements for the satellite. Secular, long-period, and short-period perturbations are then computed, with the expressions kept in closed form in both inclination and eccentricity of the satellite. Alternative expressions for short-period perturbations of high satellites are also given, assuming small values of the eccentricity. The Moon's position is specified by the inclination, node, argument of perigee, true (or mean) longitude, and its radius vector from the center of the Earth. We can then apply the results to numerical integration by using coordinates of the Moon from ephemeris tapes or to analytical representation by using results from lunar theory, with the Moon's motion represented by a precessing and rotating elliptical orbit.     

Numerical models for the study of motion of lunar satellites

The present study deals with numerical modeling of the elliptic restricted three-body problem as well as of the perturbed elliptic restricted three-body (Earth-Moon-Satellite) problem by a fourth body (Sun). Two numerical algorithms are established and investigated. The first is based on the method of the series solution of the differential equations and the second is based on a 5th-order Runge-Kutta method. The applications concern the solution of the equations and integrals of motion of the circular and elliptical restricted three-body problem as well as the search for periodic orbits of the natural satellites of the Moon in the Earth-Moon system in both cases in which the Moon describes circular or elliptical orbit around the Earth before the perturbations induced by the Sun. After the introduction of the perturbations in the Earth-Moon-Satellite system the motions of the Moon and the Satellite are studied with the same initial conditions which give periodic orbits for the unperturbed elliptic problem.