The Influence of the Earth’s Oblateness on the Energy Integral and Some Characteristics of a Spacecraft’s Orbit

Taking into account the oblateness of the Earth with the help of the second zonal harmonic of the potential of the Earth’s gravitational field modifies the integral of the energy of near-Earth motion in comparison with the Keplerian model of analysis. On the basis of this integral, an explicit analytical structure of the change in the main orbital parameter—the Keplerian energy constant—is obtained. This method of taking into account the oblateness of the Earth also makes it possible to estimate the correction, in comparison with the standard unperturbed analysis, in the departure speed of a spacecraft from the near-Earth waiting orbit when flying to the Moon or another planet.


FORMULATION OF THE PROBLEM
During a spacecraft's flight or, in general, for a material point in the central Newtonian gravitational field of the Earth, the motion of a spacecraft in a nonrotating geocentric geoequatorial coordinate system satisfies the equation (1) where r (x, y, z) is the radius vector of a point (SC; r = |r|), * is the transposition sign, and μ is the gravitational parameter of the Earth; (2) is potential of the Earth's gravity without taking into account the disturbances. System (1) corresponds to Keplerian motions in the gravitational field of a spherical, homogeneous Earth. One of the most important first integrals in this case is the energy integral [1][2][3][4][5][6]: (3) where V = |V|, V being the geocentric speed of the point; h k is the Keplerian energy constant; and element a in the case of elliptical motion when h k < 0 is the semimajor axis of the orbit. In hyperbolic motion, when h k > 0, a < 0, modulus |a| = α [3] has a geometric sense. Let us consider how integral (3) changes when taking into account the perturbation from the oblateness of the Earth as a body of revolution, as well as its use in the analysis of the orbital motion of a material point.

GENERALIZED INTEGRAL OF ENERGY TAKING INTO ACCOUNT
THE OBLATENESS OF THE EARTH In the simplest case of analyzing motion near the Earth to take into account its oblateness in potential of the Earth U to main member U 0 , second zonal harmonic U 2 is added [2,3,5,6]: where ϕ is the geocentric latitude of the spacecraft; J 2 = -C 20 is the zonal harmonic coefficient of the second order; R E is the average equatorial radius of the Earth, R E ≈ 6378.137 km; J 2 ≈ 1082.63 × 10 -6 [6]; and ε ≈ 2.63328 × 10 10 km 5 /s 2 . Equation of motion (1) changes: (6) where U(r) is potential (4), (2), (5). To the main acceleration (1), a disturbing a P , which is the gradient of the function U 2 and in rectangular coordinates it will be written in the form [2]: Theorem. The total mechanical energy of motion of a material point in model (6), where potential U(r) corresponds to (4), (2), and (5), is constant on the trajectory of the point.
Indeed, in this case, the acceleration in (6) is determined by the single-valued scalar force function U = U(r), the motion occurs in a potential field, force function U = U(r) is a potential that does not depend on time. It then follows from the general theorem of mechanics [1,4] that, on the trajectory of point dh/dt = 0, by virtue of equation of motion (6), h = const. It is possible to show this directly using (6), (1), (7). Therefore, on any trajectory of a point the energy (8) is constant, we have the first integral: We will call this integral the "generalized energy integral," bearing in mind that it generalizes integral (3) to the case of taking into account the Earth's oblateness when calculating the trajectory of a material point the SC. Remark 1. It follows from (8), (9) that, in this model of the potential, for a given energy constant h, the movement of a point occurs in the area of space = The difference from the Keplerian case is visible.
Remark 2. This approach can be applied to a more complete model of zonal harmonics. Using the model of the Earth as a body of revolution symmetric about the equatorial plane, for example, and adding a zonal harmonic of the fourth order to potential U in addition to the second one, the generalized energy integral takes the form One can also use the entire expansion of the potential in zonal harmonics, which allows one to improve the accuracy and take into account the asymmetry about the equator: Here, are nth-order Legendre polynomials. For the Earth, the following coefficients hold after the second order: J 3 = -2.53 × 10 -6 and J 4 = -1.61 × 10 -6 [6], i.e., around three orders of magnitude less than J 2 .
Remark 3. In this case, the potential axisymmetric force field is also the integral of the axial angular momentum of the point [4]: M z = (e z , [r, V]) = m = const, where e z is the ort along the axis of rotation of the Earth.

CHANGE IN THE KEPLERIAN ENERGY CONSTANT ON ORBITS OF DEPARTURE
TO THE MOON AND PLANETS Having written out integral (9) for initial point of the trajectory x 0 (r 0 , V 0 , t 0 ) and for some other point x f (r f , V f , t f ), we obtain the ratio (10) Let us apply it to the analysis of the spacecraft departure trajectories from the Earth to the Moon and planets.
It follows from (10) that the change in the Keplerian energy constants in accurate analysis satisfies the relationship (11) (11a) where i and ω are the inclination and argument of the orbit perigee and u and θ are the latitude argument and true point anomaly.
Remark 4. Energy integral (9) and exact relation (11) for the change in Δh k are valid for any orbit-in particular, for any values of inclination and eccentricity.
If we proceed from Keplerian constant h k to semiaxis a (3) and linearize according to a, then we obtain variation Δa as a first approximation from (11): During the spacecraft's flight to the planet, the departure from the Earth will occur in a hyperbolic orbit (in an osculating approximation). For purposes of assessment, we take the velocity "at infinity" for this orbit as V ∞ = 3-4 km/s. In this case, distance r f increases without bound and the penultimate term in (11), 2U 2f ,  we assume that a highly elongated near-parabolic elliptical orbit has a perigee at starting point r π = r 0 ≈ 6578 km; initial distance at apogee r α corresponds to distance to the Moon r M when the spacecraft approaches the Moon; r α ≥ (r M -Δr α ); r M ≈ 360000-405000 km; and Δr α < 0 is the correction for decreasing r α due to the oblateness of the Earth, Δr α ≈ 2Δa. In (11), distance r f increases in the process of spacecraft motion from r 0 to r M < ∞, while the penultimate term in (11) decreases to a small value (~10 -6 km 2 /s 2 ) and the limiting change in Keplerian constant h Kl is For numerical estimates, for the initial point of departure to the Moon and planets from the near-Earth reference orbit, we take, for definiteness, u 0 = 0, which is close to the characteristics of interplanetary and lunar flights. In formulas (13), (14), then, ϕ 0 = 0 and z 0 = 0, and we obtain (14a) In this case, for departures from the Earth both to the planet in a hyperbolic orbit and to the Moon in an elongated elliptical orbit, Δh Kl ≈ -0.0617 km 2 /s 2 . Numerical calculations have confirmed these estimates [7,8]. When h K < 0, this change in the Keplerian energy constant (14) corresponds to a change in the semimajor axis of the orbit, in accordance with (12): Δa ≈ -6200 km at a 0 = 200000 km, Δa ≈ -8900 km at a 0 = 220000 km, and Δa ≈ -56000 km at a 0 = 600000 km. In practice, this change occurs quickly, during the first 3 h or so of the spacecraft's initial flight from the Earth, with an increase in distance r f up to ~70000 km. aged. To simplify this procedure, using elementary identical trigonometric transformations, we lead 2U 2f to the sum of constant c 0θ and several terms of the form (2ε/p 3 ) where n j and m j are whole numbers: This constant gives the "average" shift of Keplerian energy constant h K due to the oblateness of the Earth according to (17). In this problem, parameter θ is convenient, but it is customary to do averaging over mean anomaly M. Then, after some transformations, we obtain the function value 2U 2f that is "average" over M: Remark 5. Let the celestial body not fly away from the Earth, but approach it, entering its atmosphere. The change in the Keplerian energy constant with time then occurs in the opposite direction. Moreover, distance from the Earth r decreases with time and the value of potential U 2 increases to the final value corresponding to the entry into the atmosphere, r ~ 6478 km. Remark 6. If a celestial body (spacecraft, asteroid, comet) approaches the Earth, moving at a great distance (r ~ 300000-400000 km) from the Earth in an elliptical near-parabolic orbit for which -Δh K < h K < 0, then, due to the oblateness of the Earth, the Keplerian energy constant can increase to a positive value and this body will approach the Earth in a hyperbolic orbit. Another case is possible in which, under the influence of the Earth's oblateness, the orbit changes its structure from hyperbolic to elliptical. CONCLUSIONS A simplified analysis of the effect of the oblateness of the Earth as a body of revolution-based on the zonal harmonics of the gravitational potential-makes it possible to use the energy integral, which generalizes the energy integral in the Keplerian case. This makes it possible to consider some of the qualitative features of a spacecraft's motion during flight to the Moon and plan-ets and during the return to the Earth, as well as celestial bodies, asteroids, and comets closely approaching the Earth. Oblateness of the Earth can cause a change in the structure of the orbits of these bodies-from elliptical to hyperbolic and vice versa.