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

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

Tesseral harmonic perturbations for high order and degree harmonics

A new formula has been derived for geopotential expressed in terms of orbital elements. The summation sequence was changed so that the terms of the same frequencies would be grouped and the generalized lumped coefficients were derived. The proposed formula has the same form for both odd and evenl-m.Applying Hori's perturbation method, new formulae were derived for tesseral harmonic perturbations in nonsingular orbital elements:l+g, h, e cosg,e sing, L, andH. We show the possibility of effective application of the derived formulae to the calculation of orbits of very low satellites taking into account the coefficients of tesseral harmonics of the Earth's gravitational field up to high orders and degrees. As an example the perturbations up to the order and degree of 90 for the orbit of GRM satellites were calculated. The calculations were carried out on an IBM AT personal computer.     

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.     

Theory of the motion of an artificial Earth satellite

An improved analytical solution is obtained for the motion of an artificial Earth satellite under the combined influences of gravity and atmospheric drag. The gravitational model includes zonal harmonics throughJ ₄, and the atmospheric model assumes a nonrotating spherical power density function. The differential equations are developed through second order under the assumption that the second zonal harmonic and the drag coefficient are both first-order terms, while the remaining zonal harmonics are of second order.Canonical transformations and the method of averaging are used to obtain transformations of variables which significantly simplify the transformed differential equations. A solution for these transformed equations is found; and this solution, in conjunction with the transformations cited above, gives equations for computing the six osculating orbital elements which describe the orbital motion of the satellite. The solution is valid for all eccentricities greater than 0 and less than 0.1 and all inclinations not near 0^o or the critical inclination. Approximately ninety percent of the satellites currently in orbit satisfy all these restrictions.     

Artificial Martian frozen orbits and Sun-Synchronous orbits using continuous low-thrust control

This paper analyses three types of artificial orbits around Mars pushed by continuous low-thrust control: artificial frozen orbits, artificial Sun-Synchronous orbits and artificial Sun-Synchronous frozen orbits. These artificial orbits have similar characteristics to natural frozen orbits and Sun-Synchronous orbits, and their orbital parameters can be selected arbitrarily by using continuous low-thrust control. One control strategy to achieve the artificial frozen orbit is using both the transverse and radial continuous low-thrust control, and another to achieve the artificial Sun-Synchronous orbit is using the normal continuous low-thrust control. These continuous low-thrust control strategies consider J ₂, J ₃, and J ₄ perturbations of Mars. It is proved that both control strategies can minimize characteristic velocity. Relevant formulas are derived, and numerical results are presented. Given the same initial orbital parameters, the control acceleration and characteristic velocity taking into account J ₂, J ₃, and J ₄ perturbations are similar to those taking into account J ₂ perturbations for both Mars and the Earth. The control thrust of the orbit around Mars is smaller than that around the Earth. The magnitude of the control acceleration of ASFOM-4 (named as Artificial Sun-Synchronous Frozen Orbit Method 4) is the lowest among these strategies and the characteristic velocity within one orbital period is only 0.5219 m/s for the artificial Sun-Synchronous frozen orbit around Mars. It is evident that the relationship among the control thrusts and the primary orbital parameters of Martian artificial orbits is always similar to that of the Earth. Simulation shows that the control scheme extends the orbital parameters’ selection range of three types of orbits around Mars, compared with the natural frozen orbit and Sun-Synchronous orbit.     

The motion of artificial satellites in the set of Eulerian redundant parameters

In this paper, the connections between orbit dynamics and rigid body dynamics are established throughout the Eulerian redundant parameters, the perturbation equations for any conic motion of artificial satellites are derived in terms of these parameters. A general recursive and stable computational algorithm is also established for the initial-value problem of the Eulerian parameters for satellites prediction in the Earth's gravitational field with axial symmetry. Applications of the algorithm are considered for the two cases of short and long term predictions. For the short-term prediction, we consider the problem of the final state prediction of some typical ballistic missiles in the geopotential model with zonal harmonic terms up to J ₃₆, while for the long-term prediction, we consider the perturbed J ₂ motion of Explorer 28 over 100 revolutions.     

Long-term motion of resonant satellites with arbitrary eccentricity and inclination

A first-order, semi-analytical method for the long-term motion of resonant satellites is introduced. The method provides long-term solutions, valid for nearly all eccentricities and inclinations, and for all commensurability ratios. The method allows the inclusion of all zonal and tesseral harmonics of a nonspherical planet.We present here an application of the method to a synchronous satellite includingonly theJ ₂ andJ ₂₂ harmonics. Global, long-term solutions for this problem are given for arbitrary values of eccentricity, argument of perigee and inclination.     

Effect of the Earth’s compression on the physical libration of the Moon

The effect of the Earth??s compression on the physical libration of the Moon is studied using a new vector method. The moment of gravitational forces exerted on the Moon by the oblate Earth is derived considering second order harmonics. The terms in the expression for this moment are arranged according to their order of magnitude. The contribution due to a spherically symmetric Earth proves to be greater by a factor of 1.34 × 10⁶ than a typical term allowing for the oblateness. A linearized Euler system of equations to describe the Moon??s rotation with allowance for external gravitational forces is given. A full solution of the differential equation describing the Moon??s libration in longitude is derived. This solution includes both arbitrary and forced oscillation harmonics that we studied earlier (perturbations due to a spherically symmetric Earth and the Sun) and new harmonics due to the Earth??s compression. We posed and solved the problem of spinorbital motion considering the orientation of the Earth??s rotation axis with regard to the axes of inertia of the Moon when it is at a random point in its orbit. The rotation axes of the Earth and the Moon are shown to become coplanar with each other when the orbiting Moon has an ecliptic longitude of L _? = 90° or L _? = 270°. The famous Cassini??s laws describing the motion of the Moon are supplemented by the rule for coplanarity when proper rotations in the Earth-Moon system are taken into account. When we consider the effect of the Earth??s compression on the Moon??s libration in longitude, a harmonic with an amplitude of 0.03?? and period of T ₈ = 9.300 Julian years appears. This amplitude exceeds the most noticeable harmonic due to the Sun by a factor of nearly 2.7. The effect of the Earth??s compression on the variation in spin angular velocity of the Moon proves to be negligible.     

Multi-SunSynchronous Orbits in the Solar System

Performances of a planetary observation system are strongly related to the choice of the orbit used. Trajectories with characteristics of periodicity are very useful for the assessment of time-varying phenomena and thus Periodic SunSynchronous and Periodic Multi-SunSynchronous Orbits are particularly suitable to this end. In this paper, the research into these kinds of orbits, previously proposed for the Earth and Mars, has been extended to planets of the Solar System and to their principal moons. In general, these trajectories are typically obtained under the hypothesis that the J₂ harmonic is predominant with respect to the other orbital perturbations, since this allows an analytical solution. However, the hypothesis of J₂ predominant is not always verified in the Solar System and so analytical techniques must be replaced by numerical simulations. Interesting results have been obtained for the planets Mars and Jupiter and for the moons Europa, Callisto and Titan, where periodic trajectories with reduced revisit times and low altitudes have been found. These solutions allow the observation of time-varying phenomena with high spatial and temporal resolution.     

Analytical solutions for J 2-perturbed unbounded equatorial orbits

While solutions for bounded orbits about oblate spheroidal planets have been presented before, similar solutions for unbounded motion are scarce. This paper develops solutions for unbounded motion in the equatorial plane of an oblate spheroidal planet, while taking into account only the J ₂ harmonic in the gravitational potential. Two cases are distinguished: A pseudo-parabolic motion, obtained for zero total specific energy, and a pseudo-hyperbolic motion, characterized by positive total specific energy. The solutions to the equations of motion are expressed using elliptic integrals. The pseudo-parabolic motion unveils a new orbit, termed herein the fish orbit, which has not been observed thus far in the perturbed two-body problem. The pseudo-hyperbolic solutions show that significant differences exist between the Keplerian flyby and the flyby performed under the the J ₂ zonal harmonic. Numerical simulations are used to quantify these differences.     

A New Reference Equipotential Surface, and Reference Ellipsoid for the Planet Mars

Using the shape model of Mars GTM090AA in terms of spherical harmonics complete to degree and order 90 and gravitational field model of Mars GGM2BC80 in terms of spherical harmonics complete to degree and order 80, both from Mars Global Surveyor (MGS) mission, the geometry (shape) and gravity potential value of reference equipotential surface of Mars (Areoid) are computed based on a constrained optimization problem. In this paper, the Areoid is defined as a reference equipotential surface, which best fits to the shape of Mars in least squares sense. The estimated gravity potential value of the Areoid from this study, i.e. W ₀ = (12,654,875 ± 69) (m²/s²), is used as one of the four fundamental gravity parameters of Mars namely, {W ₀, GM, ω, J ₂₀}, i.e. {Areoid’s gravity potential, gravitational constant of Mars, angular velocity of Mars, second zonal spherical harmonic of gravitational field expansion of Mars}, to compute a bi-axial reference ellipsoid of Somigliana-Pizzetti type as the hydrostatic approximate figure of Mars. The estimated values of semi-major and semi-minor axis of the computed reference ellipsoid of Mars are (3,395,428 ± 19) (m), and (3,377,678 ± 19) (m), respectively. Finally the computed Areoid is presented with respect to the computed reference ellipsoid.     

Long period variations of the motion of a satellite due to non-resonant tesseral harmonics of a gravity potential

The main effects of tesseral harmonics of a gravity potential expansion on the motion of a satellite, are short period variations as well as long period variations due to resonances. However, other smaller long period and secular variations can arise from interactions between tesseral terms of the same order. The analytical integration of these effects is developed, using numerical evaluation of Kaula eccentricity and inclination functions. Examples for some Earth's geodetic satellites show that secular effects can reach a few decameters per year. The secular variations can even reach several hundred of meters per year for the Mars natural satellite Phobos.     

Analytical approach using KS elements to short-term orbit predictions including J 2

Analytical solutions using KS elements are derived. The perturbation considered is the Earth's zonal harmonic J ₂. The series expansions include terms of fourth power in the eccentricity. Only two of the nine KS element equations are integrated analytically due to the reasons of symmetry. The analytical solution is suitable for short-term orbit computations. Numerical studies show that reasonably good estimates of the orbital elements can be obtained in one step of 10 to 30 degrees of eccentric anomaly for near-Earth orbits of moderate eccentricity. For application purposes, the analytical solution can be effectively used for onboard computation in the navigation and guidance packages, where the modelling of J ₂ effect becomes necessary.     

Regional positioning using a low Earth orbit satellite constellation

Global and regional satellite navigation systems are constellations orbiting the Earth and transmitting radio signals for determining position and velocity of users around the globe. The state-of-the-art navigation satellite systems are located in medium Earth orbits and geosynchronous Earth orbits and are characterized by high launching, building and maintenance costs. For applications that require only regional coverage, the continuous and global coverage that existing systems provide may be unnecessary. Thus, a nano-satellites-based regional navigation satellite system in Low Earth Orbit (LEO), with significantly reduced launching, building and maintenance costs, can be considered. Thus, this paper is aimed at developing a LEO constellation optimization and design method, using genetic algorithms and gradient-based optimization. The preliminary results of this study include 268 LEO constellations, aimed at regional navigation in an approximately 1000 km \(\times \) 1000 km area centered at the geographic coordinates [30, 30] degrees. The constellations performance is examined using simulations, and the figures of merit include total coverage time, revisit time, and geometric dilution of precision (GDOP) percentiles. The GDOP is a quantity that determines the positioning solution accuracy and solely depends on the spatial geometry of the satellites. Whereas the optimization method takes into account only the Earth’s second zonal harmonic coefficient, the simulations include the Earth’s gravitational field with zonal and tesseral harmonics up to degree 10 and order 10, Solar radiation pressure, drag, and the lunisolar gravitational perturbation.     

Analytical solutions to the four post-Newtonian effects in a near-Earth satellite orbit

On the basis of the results by Huang et al. (1990), this paper further discusses and analyses the four post-Newtonian effects in a near-Earth satellite orbit: the Schwarzschild solution, the post-Newtonian effects of the geodesic precession, the Lense-Thirring precession and the oblateness of the Earth. A full analytical solution to the effects including their direct perturbations and mixed perturbations due to the Newtonian oblateness (J ₂) perturbation and the Schwarzschild solution is obtained using the quasi-mean orbital element method analogous to the Kozai's mean orbital element one. Some perturbation properties of the post-Newtonian effects are revealed. The results obtained not only can provide a sound scientific basis for the precise determination of a man-made satellite orbit but also is suitable for similar mechanics systems, such as the motions of planets, asteroids and natural satellites.     

Models, figures and gravitational moments of Jupiter’s satellite Io: Effects of the second order approximation

The effects of the equilibrium figure theory to within terms of the second order in a small parameter α on figure parameters and gravitational moments of the Galilean satellite Io have been considered. Integro-differential equations of the theory of figure to second order have been first solved numerically. Relations between the low-order coefficients of the gravitational field for satellites in hydrostatic equilibrium are generalized according to the second order theory. To show the effects of the second approximation, two three-layer trial models of Io are used. The considered models of the Io’s interiors differ by the size and density of the core, while having the same thickness and density of the crust, and the mantle density difference is only 20 kg/m³. The corrections of second order in smallness to the gravitational moments J₂ and C₂₂ decrease the third decimal digit of model gravitational moments by two units. As the effects of third and forth harmonics are determined mostly by outer layers of Io, to distinguish between model mantle density, the gravitational moments J₄, C₄₂ and C₄₄ should be determined to accuracy with three or four decimal digits. The second order corrections mostly effect the semi-axis a, and less the semi-axes b and c.     

The effect of the dynamical parameters in the motion of the Galilean satellites of Jupiter

The construction of an analytical theory of the motion of the Galilean satellites of Jupiter requires that we keep track of the dynamical parameters, that is, the masses of the satellites, and the harmonic coefficients of the potential of the planet J₂ and J₄. This is realized here. But as in other theories the solution becomes partly numerical from the resolution of an autonomous system. The aim of this paper is to present a method to obtain developped solutions of this autonomous system. In these solutions the proper motions of the pericenters and nodes are obtained as short series developped in the neighbourhood of a numerical solution. We have used these results to obtain complementary terms in the general solution which give a complete representation of the motions with respect to the dynamical parameters.     

On the computation of the spherical harmonic terms needed during the numerical integration of the orbital motion of an artificial satellite

A method is presented for the accurate and efficient computation of the forces and their first derivatives arising from any number of zonal and tesseral terms in the Earth's gravitational potential. The basic formulae are recurrence relations between some solid spherical harmonics,V _n,m, associated with the standard polynomial ones.     

On error bounds and initialization in satellite orbit theories

The order of magnitude of the error is investigated for a first-order von Zeipel theory of satellite orbits in an axisymmetric force field, i.e., first-order long period and short-period effects are included along with second order secular rates. The treatment is valid for zero eccentricity and/or inclination. In the case where initial position and velocity vectors are known, the in-track position error over time intervals of order 1/J ₂ is kept at 0(J ₂ ²), like the other position errors and velocity errors, by calibration of the mean motion with the aid of the energy integral. The results are specifically applicable to accuracy comparisons of the Brouwer orbit prediction method with numerical integration. A modified calibration is presented for the general asymmetric force field which includes tesseral harmonics.     

Analytical Approach to the Motion of a Lunar Artificial Satellite

The canonical equations of motion of an artificial lunar satellite are formulated including the effects of the asphericity of the Moon comprising the harmonics J ₂, J ₂₂, J ₃, J ₃₁, J ₄ andJ ₅, the oblateness of the Earth up to the second zonal harmonic, as well as the disturbing function due to the attractions of the Earth and of the Sun (terms are retained up to order 10^-6 for the higher orbits and 10^-8 for the lower orbits). This revised version was published online in July 2006 with corrections to the Cover Date.