## Orbit Due To Stronger Central Force

3 a lal EFFECT OF NONPOINT MASS FORCES

### Ibl PERIGEE ADVANCE FROM AN EQUATORIAL ORBIT

Fig. 3-19. Motion of the Lines of Apsides Due to the Earth's Oblateness less, which causes the line of apsides to rotate opposite the direction of motion while the satellite is over the pole. In a polar orbit (i'=90 deg), a satellite spends part of the time over the equator where the gravitational force is larger than the mean and part of the time over the poles where it is smaller. It can be shown that for polar orbits the net «ffect is rotation of the line of apsides opposite the direction of motion, although the rate of rotation is less than for equatorial orbits.

As described in Section 5.2 and Appendix G, the oblateness of the Earth is treated analytically by expanding the gravitational potential in a series of spherical harmonics. The first term in the expansion provides the force resulting from a point mass. The second term, called J2 perturbations (from the second or J2 coefficient in the expansion), represents the modification to the mean force due to the oblateness of the Earth. Higher order terms represent effects due to the deviation of the shape of the Earth from a simple ellipsoid. Values of the principal coefficients are given in Appendix L; additional coefficients are given by Allen [1973].

The rate of change of the orbital elements arising from the J2 perturbation is conveniently expressed in terms of the mean angular motion, n, equal to the rate of change of mean anomaly (see, for example, Roy [1965]):

where the nominal mean angular motion, Rq, is given from Table 3-2 by

The anomalistic period (i.e., perigee to perigee) is

and the mean anomaly at time t is

where M0 is the mean anomaly at the epoch /„. Note that the mean anomaly is measured relative to the moving perigee. Each of these equations reduces to the nonperturbed form when J2~0.

When only the secular variations due to the J2 term are considered, the right ascension of the ascending node at time t is given by:

= a0- 2.06474 x 1014a~7/2( 1 - e2)~2(cos i)(f - i0) + 6 (J22) (3-41d)

where (/ — tg) is in days, Q is in degrees, a is the semimajor axis in kilometres, e is the eccentricity, i is the inclination, and is £2 at the epoch time /„. In the final two forms above, n is approximated by n^ which may produce an error as large as 0.1% in the dSl/dt term.

If the product of the three terms in a, e, and / in Eq. (3-4Id) has the value - 4.7737 X10"15 km~7/2, the rotation rate of the node will be =0.9856 deg/day or 1 rotation per year. Such an orbit is called Sun synchronous because the orientation of the orbit plane will remain nearly fixed relative to the Sun as the Earth moves in its orbit. Thus, the spacecraft will continuously view the surface of the Earth at the same local time at any given latitude. For any satellite in a circular orbit, the local mean time*, T, at which the satellite is over latitude A is given by (see Appendix J and Eq. (A-6)):