The effect of the third-body perturbation is one of the most important mechanisms in delivering the major Earth orbiting objects into the regions where the atmosphere starts to decay. Several papers have studied the third-body perturbation of the Sun and the Moon on a satellite orbiting the Earth and major related works are listed as follows. In a series of papers (Kozai 1959, 1962, 1965), Kozai developed the secular and long period terms of the disturbing function due to lunisolar perturbations in terms of the orbital elements of a satellite, the Sun and the Moon, and applied those secular perturbations to asteroids with high inclination and eccentricity assuming that Jupiter is the perturbing body. Blitzer (1959) obtained the estimates of secular terms due to lunisolar disturbances. In the sequence, Cook (1962) used the Lagrange’s planetary equations to obtain expressions for the variations of elements during the revolution of a satellite and for the rates of variation of the elements. Giacaglia (1974) obtained the disturbing function for the disturbance of the Moon using the ecliptic elements for the Moon and the equatorial elements for a satellite. Hough (1981) studied the effects of lunisolar disturbances in orbits close to critical inclinations from the point of view of the geopotential of the Earth and concluded that the effects are significant in high altitudes. Broucke (2003) and Prado & Costa (1998), truncated the disturbing function of the third body after the second and the fourth order expansion in Legendre polynomials, respectively. After that, Costa (1998) expanded the order of this model to the eighth order. Solόrzano & Prado (2007) developed a semi-analytical study of the perturbation caused on a spacecraft by a third body with a single averaged model to eliminate the terms due to the short-lived periodic motion of the spacecraft. Prado (2003) studied the evolution of retrograde orbits around the Earth. The existence of circular, equatorial and frozen orbits was also considered. Costa & Prado (2002) discussed the critical angle of a third-body perturbation that is a value for the inclination such that any near-circular orbit with inclination below this value remains near- circular. Domingos et al. (2008) showed an analytical expansion to study the third-body perturbation for a case where the perturbing body is in an elliptical orbit, based on the expansion of the perturbing function in polynomials of Legendre.
In this paper, the study done by Broucke (2003), Prado (2003) and Domingos et al. (2008) is expanded, including the eccentricity and inclination of the perturbing body. A three-dimensional model is used to study the effects of a third-body gravitational potential on the spacecraft that orbits near a celestial body including the eccentricity and the inclination of the third-body in the formulation. There is no assumption in our study on the perturbed and perturbing orbits except the usual assumption that there are only three bodies involved in the system: the first body with mass,
This paper is organized as follows. In Section 2, we formulate the disturbing function (using the traditional expansion in Legendre polynomials) due to a third body truncated up to the second-order term. In Section 3, The Lagrange’s planetary equations for the variations of spacecraft’s orbital elements are shown. In Section 4, we simulate the study assuming a spacecraft in an elliptical inclined orbit around the Moon with its motion perturbed by the Earth.
About notations through this paper, we use the wellknown Keplerian elements: the semi-major axis,
For the models used in this research, it is necessary to calculate the parts of the disturbing function due to
The magnitude of the effect of the disturbing body on the orbit of a spacecraft depends on the position of the disturbing body in its orbit. The next step is to express cos
The last relation can be rewritten in a compact form as below:
where the non-vanishing coefficients are
The next step is to average expression (2), using Eq. (3), over the short period of the spacecraft as well as with respect to the distant perturbing body. The standard definition for the average can be written as:
Rewriting above equation gives
where the non-vanishing coefficients are
Using Lagrange planetary equations, the variation of the orbital elements due to the averaging disturbing potential «
In this section, we suppose a spacecraft in an elliptical inclined orbit around the Moon and its motion perturbed by the Earth for a duration time of 1,000 days. The initial conditions found in Prado (2003) and Domingos et al. (2008) ; a lunar satellite with
[Fig. 2.] Evolutions of orbital elements for spacecraft’s over 1,000 days with initial conditions a0=0.01 (3,844km), e0=0.01, i0=15°, and ω0=Ω0=0°. It is supposed that the initial conditions of the third body are a’= (384,400km), e’=0.1, and ω’=Ω’=0°, and 10° ≤ i ’ ≤ 80°, inc= i ’. (a) Eccentricity, (b) Inclination (rad), (c) Argument of periapsis (rad), (d) Right ascension of ascending node (rad)
[Fig. 3.] Evolutions of orbital elements for spacecraft’s over 1,000 days with initial conditions a0=0.01 (3,844km), e0=0.01, i0=80°, and ω0=Ω0=0°. It is supposed that the initial conditions of the third body are a’= (384,400km), e’=0.1, and ω’=Ω’=0°, and 10° ≤ i ’ ≤ 80°, inc= i ’. (a) Eccentricity, (b) Inclination (rad), (c) Argument of periapsis (rad), (d) Right ascension of ascending node (rad)
Beginning this study, the basic question is; what is the effect of inclined third body on the spacecraft orbital elements? For the answer, the effects of wide range of the inclinations of the third body (10°-80°) are presented. Using the initial conditions found in Prado (2003) and Domingos et al. (2008), the results of changing the third body’s inclinations over (10°-80°) on the model indicated that the effect of the third body’s inclinations cannot be neglected to be consistent with many previous studies.
In this paper, we have used three-dimensional model to analyze the perturbations of the third body potential on a spacecraft that orbits close to a heavenly body including the eccentricity and the inclination of the third-body of the perturbing body within the formulation.
All curves turned out to be smooth through the well-known double-average model based on the mean motion of the spacecraft and the perturbing body which plays an important role in the analysis of spacecrafts long term dynamics enabling the short-lived periodic terms to be eliminated.