The lunar orbit revisited,III

In this paper we present an investigation on the tidal evolution of a system of three bodies: the Earth, the Moon and the Sun. Equations are derived including dissipation in the planet caused by the tidal interaction between the planet and the satellite and between the planet and the sun. Dissipation within the Moon is included as well. The set of differential equations obtained is valid as long as the solar disturbances dominate the perturbations on the satellite's motion due to the oblateness of the planet, namelya/R _e greater than 15, and closer than that point equations derived in a preceding paper are used.The result shows the Moon was closer to the Earth in the past than now with an inclination to the ecliptic greater than today, whereas the obliquity was smaller. Toward the past, the inclination to the Earth's equator begins decreasing to 12° fora/R _e=12 and suddenly grows. During the first stage the results are weakly dependant on the magnitude of the dissipation within the satellite, whereas the distance of the closest approach and the prior history are strongly dependent on that dissipation. In particular, the crossing of the Roche limit can be avoided.     

Fast Earth–Moon transfers with ballistic capture

This contribution deals with fast Earth–Moon transfers with ballistic capture in the patched three-body model. We compute ensembles of preliminary solutions using a model that takes into account the relative inclination of the orbital planes of the primaries. The ballistic capture orbits around the Moon are obtained relying on the hyperbolic invariant structures associated to the collinear Lagrangian points of the Earth–Moon system, and the Sun–Earth system portion of the transfers are quasi-periodic orbits obtained by a genetic algorithm. The trajectories are designed to be good initial guesses to search optimal cost-efficient short-time Earth–Moon transfers with ballistic capture in more realistic models.     

The effect of the Earth’s oblateness on the Moon’s physical libration in latitude

The Moon’s physical libration in latitude generated by gravitational forces caused by the Earth’s oblateness has been examined by a vector analytical method. Libration oscillations are described by a close set of five linear inhomogeneous differential equations, the dispersion equation has five roots, one of which is zero. A complete solution is obtained. It is revealed that the Earth’s oblateness: a) has little effect on the instantaneous axis of Moon’s rotation, but causes an oscillatory rotation of the body of the Moon with an amplitude of 0.072″ and pulsation period of 16.88 Julian years; b) causes small nutations of poles of the orbit and of the ecliptic along tight spirals, which occupy a disk with a cut in a center and with radius of 0.072″. Perturbations caused by the spherical Earth generate: a) physical librations in latitude with an amplitude of 34.275″; b) nutational motion for centers of small spiral nutations of orbit (ecliptic) pole over ellipses with semi-major axes of 113.850″ (85.158″) and the first pole rotates round the second one along a circle with radius of 28.691″; c) nutation of the Moon’s celestial pole over an ellipse with a semi-major axis of 45.04″ and with an axes ratio of about 0.004 with a period of T = 27.212 days. The principal ellipse’s axis is directed tangentially with respect to the precession circumference, along which the celestial pole moves nonuniformly nearly in one dimension. In contrast to the accepted concept, the latitude does not change while the Moon’s poles of rotation move. The dynamical reason for the inclination of the Moon’s mean equator with respect to the ecliptic is oblateness of the body of the Moon.     

The evolution of the lunar orbit revisited,II

We present here a model for the tidal evolution of an isolated two-body system. Equations are derived, including the dissipation in the planet as in the satellite, in a frequency dependent lag model. The set of differential equations obtained is still valid for large eccentricity, as well as for all inclinations. The reference plane chosen enables us to study the evolution for both the orbital plane and the equatorial plane.The results obtained show the Moon, after having approached the Earth with small variations for the inclination and the eccentricity, exhibits strong increase for the two parameters in the vicinity of the closest approach. In every case the eccentricity tends towards the value 1, whereas the variations of the in clinations are dependent on the magnitude of the dissipation in the satellite.Some qualitative results are also investigated for the final behaviour of satellites such as Triton and the Galilean satellites.     

Earth satellites in resonance with the Moon and the Sun as objects of laser ranging. Analytical solution for their motion

Distant Earth satellites repeating nearly periodically their configurations both with the Moon and the Sun may appear to be even more convenient carriers of laser retroreflectors than the Moon, when geodetic applications are of primary interest. The analytical solution for the motion of a satellite in resonance both with the Moon and the Sun has been outlined in this paper, the periodic orbit of the planar restricted four-body problem being taken as an intermediary. The Von Zeipel transformation gives the Hamiltonian not depending on the fast variables. The stationary solution for this Hamiltonian has been found. Then the non-homogeneous variation equations have been formed, taking into account the orbital eccentricities of the Moon and the Sun. The solution of these equations has been obtained and its accuracy has been tested by numerical integration.Presented at the XXII Congress of the International Astronautical Federation, Brussels, Belgium September 20–25, 1971.     

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.     

Fuel optimization for low-thrust Earth–Moon transfer via indirect optimal control

The problem of designing low-energy transfers between the Earth and the Moon has attracted recently a major interest from the scientific community. In this paper, an indirect optimal control approach is used to determine minimum-fuel low-thrust transfers between a low Earth orbit and a Lunar orbit in the Sun–Earth–Moon Bicircular Restricted Four-Body Problem. First, the optimal control problem is formulated and its necessary optimality conditions are derived from Pontryagin’s Maximum Principle. Then, two different solution methods are proposed to overcome the numerical difficulties arising from the huge sensitivity of the problem’s state and costate equations. The first one consists in the use of continuation techniques. The second one is based on a massive exploration of the set of unknown variables appearing in the optimality conditions. The dimension of the search space is reduced by considering adapted variables leading to a reduction of the computational time. The trajectories found are classified in several families according to their shape, transfer duration and fuel expenditure. Finally, an analysis based on the dynamical structure provided by the invariant manifolds of the two underlying Circular Restricted Three-Body Problems, Earth–Moon and Sun–Earth is presented leading to a physical interpretation of the different families of trajectories.     

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.     

Computation of weak stability boundaries: Sun–Jupiter system

We compute the weak stability boundary in the planar circular restricted three-body problem starting from the algorithmic definition, and its generalization by García and Gómez. In addition, we consider a new set of primaries, Sun–Jupiter, to replace the case of Earth–Moon considered in previous studies. Numerical enhancements are described and compared to previous methods. This includes defining the equations of motion in polar coordinates and a modified numerical scheme for the derivation of both stable sets and their boundaries. These enhancements decrease the computational time. New results are obtained by considering the Sun–Jupiter case which we compare to the Earth–Moon case.     

Toward the construction of a new precession-nutation theory of nonrigid Earth

The accuracy of the rigid Earth solution SMART97 is 2?μ as over the time interval (1968, 2023), accuracy showed by the comparison with a numerical integration using the positions of the Moon, the Sun, and the planets given by DE403. To obtain a nonrigid Earth solution, we use the transfer function of Mathews et al. (2000) and , to keep the precision of our rigid Earth solution in the computation of the geophysical effects, we apply this transfer function to the Earth's angular velocity vector in order to avoid the inherent approximations of the classical methods. Moreover the perturbations of the third component of the angular velocity vector are taken into account. Lastly, we take into account, in an iterative process, the second order perturbations due to the geophysical effects. The results are compared with the Herring solution (1996) published in the IERS Conventions.     

Effect of Moon perturbation on the energy curves and equilibrium points in the Sun–Earth–Moon system

In this paper, we have considered that the Moon motion around the Earth is a source of a perturbation for the infinitesimal body motion in the Sun–Earth system. The perturbation effect is analyzed by using the Sun–Earth–Moon bi–circular model (BCM). We have determined the effect of this perturbation on the Lagrangian points and zero velocity curves. We have obtained the motion of infinitesimal body in the neighborhood of the equivalent equilibria of the triangular equilibrium points. Moreover, to know the nature of the trajectory, we have estimated the first order Lyapunov characteristic exponents of the trajectory emanating from the vicinity of the triangular equilibrium point in the proposed system. It is noticed that due to the generated perturbation by the Moon motion, the results are affected significantly, and the Jacobian constant is fluctuated periodically as the Moon is moving around the Earth. Finally, we emphasize that this model could be applicable to send either satellite or telescope for deep space exploration.     

Geostationary secular dynamics revisited: application to high area-to-mass ratio objects

The long-term dynamics of the geostationary Earth orbits (GEO) is revisited through the application of canonical perturbation theory. We consider a Hamiltonian model accounting for all major perturbations: geopotential at order and degree two, lunisolar perturbations with a realistic model for the Sun and Moon orbits, and solar radiation pressure. The long-term dynamics of the GEO region has been studied both numerically and analytically, in view of the relevance of such studies to the issue of space debris or to the disposal of GEO satellites. Past studies focused on the orbital evolution of objects around a nominal solution, hereafter called the forced equilibrium solution, which shows a particularly strong dependence on the area-to-mass ratio. Here, we (i) give theoretical estimates for the long-term behavior of such orbits, and (ii) we examine the nature of the forced equilibrium itself. In the lowest approximation, the forced equilibrium implies motion with a constant non-zero average ‘forced eccentricity’, as well as a constant non-zero average inclination, otherwise known in satellite dynamics as the inclination of the invariant ‘Laplace plane’. Using a higher order normal form, we demonstrate that this equilibrium actually represents not a point in phase space, but a trajectory taking place on a lower-dimensional torus. We give analytical expressions for this special trajectory, and we compare our results to those found by numerical orbit propagation. We finally discuss the use of proper elements, i.e., approximate integrals of motion for the GEO orbits.     

月球卫星轨道力学综述

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

Using space manifold dynamics to deploy a small satellite constellation around the Moon

The possibility of communicating with the far side of the Moon is essential for keeping a continuous radio link with lunar orbiting spacecraft, as well as with manned or unmanned surface facilities in locations characterized by poor coverage from Earth. If the exploration and the exploitation of the Moon must be sustainable in the medium/long term, we need to develop the capability to realize and service such facilities at an affordable cost. Minimizing the spacecraft mass and the number of launches is a driving parameter to this end. The aim of this study is to show how Space Manifold Dynamics can be profitably applied in order to launch three small spacecraft onboard the same launch vehicle and send them to different orbits around the Moon with no significant difference in the Delta-V budgets. Internal manifold transfers are considered to minimize also the transfer time. The approach used is the following: we used the linearized solution of the equations of motion in the Circular Restricted Three Body Problem to determine a first–guess state vector associated with the Weak Stability Boundary regions, either around the collinear Lagrangian point L1 or around the Moon. The resulting vector is then used as initial state in a numerical backward-integration sequence that outputs a trajectory on a manifold. The dynamical model used in the numerical integration is four-body and non-circular, i.e. the perturbations of the Sun and the lunar orbital eccentricity are accounted for. The trajectory found in this way is used as the principal segment of the lunar transfer. After separation, with minor maneuvers each satellite is injected into different orbits that lead to ballistic capture around the Moon. Finally, one or more circularization maneuvers are needed in order to achieve the final circular orbits. The whole mission profile, from launch to insertion into the final lunar orbits, is modeled numerically with the commercial software STK.     

Sun-synchronous orbits near critical inclination

The long period dynamics of Sun-synchronous orbits near the critical inclination 116.6° are investigated. It is known that, at the critical inclination, the average perigee location is unchanged by Earth oblateness. For certain values of semimajor axis and eccentricity, orbit plane precession caused by Earth oblateness is synchronous with the mean orbital motion of the apparent Sun (a Sun-synchronism). Sun-synchronous orbits have been used extensively in meteorological and remote sensing satellite missions. Gravitational perturbations arising from an aspherical Earth, the Moon, and the Sun cause long period fluctuations in the mean argument of perigee, eccentricity, inclination, and ascending node. Double resonance occurs because slow oscillations in the perigee and Sun-referenced ascending node are coupled through the solar gravity gradient. It is shown that the total number and infinitesimal stability of equilibrium solutions can change abruptly over the Sun-synchronous range of semimajor axis values (1.54 to 1.70 Earth radii). The effect of direct solar radiation pressure upon certain stable equilibria is investigated.     

月球探测器过渡轨道的短弧定轨方法

论述的短弧定轨,是指在无先验信息情况下又避开多变元迭代的初轨计算方法,它需要相应的动力学问题有一能反映短弧内达到一定精度的近似分析解．探测器进入月球引力作用范围后接近月球时可以处理成相对月球的受摄二体问题,而在地球附近,则可处理成相对地球的受摄二体问题,但在整个过渡段的力模型只能处理成一个受摄的限制性三体问题．而限制性三体问题无分析解,即使在月球引力作用范围外,对于大推力脉冲式的过渡方式,相对地球的变化椭圆轨道的偏心率很大(超过Laplace极限),在考虑月球引力摄动时亦无法构造摄动分析解．就此问题,考虑在地球非球形引力(只包含J2项)和月球引力共同作用下,构造了探测器飞抵月球过渡轨道段的时间幂级数解,在此基础上给出一种受摄二体问题意义下的初轨计算方法,经数值验证,定轨方法有效,可供地面测控系统参考．     

Analytical theory of a lunar artificial satellite with third body perturbations

We present here the first numerical results of our analytical theory of an artificial satellite of the Moon. The perturbation method used is the Lie Transform for averaging the Hamiltonian of the problem, in canonical variables: short-period terms (linked to l, the mean anomaly) are eliminated first. We achieved a quite complete averaged model with the main four perturbations, which are: the synchronous rotation of the Moon (rate ), the oblateness J ₂ of the Moon, the triaxiality C ₂₂ of the Moon ( ) and the major third body effect of the Earth (ELP2000). The solution is developed in powers of small factors linked to these perturbations up to second-order; the initial perturbations being sorted ( is first-order while the others are second-order). The results are obtained in a closed form, without any series developments in eccentricity nor inclination, so the solution apply for a wide range of values. Numerical integrations are performed in order to validate our analytical theory. The effect of each perturbation is presented progressively and separately as far as possible, in order to achieve a better understanding of the underlying mechanisms. We also highlight the important fact that it is necessary to adapt the initial conditions from averaged to osculating values in order to validate our averaged model dedicated to mission analysis purposes.     

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.     

A motivating exploration on lunar craters and low-energy dynamics in the Earth–Moon system

It is known that most of the craters on the surface of the Moon were created by the collision of minor bodies of the Solar System. Main Belt Asteroids, which can approach the terrestrial planets as a consequence of different types of resonance, are actually the main responsible for this phenomenon. Our aim is to investigate the impact distributions on the lunar surface that low-energy dynamics can provide. As a first approximation, we exploit the hyberbolic invariant manifolds associated with the central invariant manifold around the equilibrium point L ₂ of the Earth–Moon system within the framework of the Circular Restricted Three-Body Problem. Taking transit trajectories at several energy levels, we look for orbits intersecting the surface of the Moon and we attempt to define a relationship between longitude and latitude of arrival and lunar craters density. Then, we add the gravitational effect of the Sun by considering the Bicircular Restricted Four-Body Problem. In the former case, as main outcome, we observe a more relevant bombardment at the apex of the lunar surface, and a percentage of impact which is almost constant and whose value depends on the assumed Earth–Moon distance d_EM. In the latter, it seems that the Earth–Moon and Earth–Moon–Sun relative distances and the initial phase of the Sun θ ₀ play a crucial role on the impact distribution. The leading side focusing becomes more and more evident as d_EM decreases and there seems to exist values of θ ₀ more favorable to produce impacts with the Moon. Moreover, the presence of the Sun makes some trajectories to collide with the Earth. The corresponding quantity floats between 1 and 5 percent. As further exploration, we assume an uniform density of impact on the lunar surface, looking for the regions in the Earth–Moon neighbourhood these colliding trajectories have to come from. It turns out that low-energy ejecta originated from high-energy impacts are also responsible of the phenomenon we are considering.