## X

Fig. 5-12. Denning the Orientation of an Orbit in Space. (See also Fig. 3-7.)

where Q, u, and < are the longitude of the ascending node, the argument of perigee, and the inclination; and x,y, and z are resolved in the coordinate system in which the elements are defined. It is possible to use Eqs. (5-33) and (5-34) to remove the explicit true anomaly dependence in Eqs. (5-36) through (5-38) and thus compute the position directly.

The velocity at time t may be found by applying the chain rule to Eqs. (5-36) through (5-38). Thus, for example, dx dx _3r_ 9£ 3x jip_ 3E ,, «q.

where u, 0, and < are assumed constant. The quantity dE/dt may be found by differentiating Kepler's equation to obtain lir-f ~«°^)d37 <«»>

The equations for the velocity are as follows:

I, = cos Q cos « — sin i2 sinco cos i m. =sini2cosa) + cosi2sintocos/

.mil n^sinwsin/

l2 — — cos fl sin w—sin fleos w cosí m2 — - sin fl sin « + cos fleos to cos / n2 = coswsin»

The above procedure is not coordinate-system-dependent; that is, the position and velocity will be in whatever coordinate system the elements of the orbit are defined. For Earth-orbiting spacecraft, the elements are usually given in geocentric inertia! coordinates, whereas for interplanetary orbits, the elements are usually given in heliocentric coordinates (see Section 2.2).

The simple two-body orbit generator described above may be modified to take into account noncentral forces (or forces from a third body). The procedure, which is known as the method of general perturbations (Section 5.5), is to obtain series solutions to the equations of motion in the form of perturbations to the orbit elements. These elements then become functions of time, and the method outlined above for solving Kepler's equation is applied with different elements each time the position is to be calculated. (Some simple results of general perturbations, in showing the effect of the oblateness of the Earth on the orbit of an Earth satellite, are given in Section 3.4, Eqs. (3-37), (3-38), (3-42), and (3-43).) Increasingly accurate descriptions of the actual motion of a spacecraft can be obtained by including an increasing number of perturbation terms (both periodic and secular). This forms the basis of the Brouwer method, which is a detailed application of the theory of general perturbations to the motion of artificial Earth satellites.

A simpler, more direct approach to the detailed calculation of orbits is to integrate the equations of motion directly, given the initial conditions. This is known as the method of special perturbations (Section 5.5). Models of all forces which are expected to be significant are included. Two commonly used integration schemes are those of Cowell (integration in rectangular coordinates) and Encke (calculation of an osculating conic section for which integration gives the differences between the real coordinates and the coordinates given by the conic section). Methods of this type are used in GTDS.

Utility and Accuracy of Two-Body Orbit Generators. The question arises as to the utility and accuracy of simple two-body orbit generators. Some comparisons, using as a reference an integration model from GTDS incorporating Sun, Moon, Earth harmonics, atmospheric drag, and solar radiation pressure are shown in Table 5-11. The references in Table 5-11 to "error with and without perturbations" refer to the inclusion in the two-body generator of the secular .perturbations in the right ascension of the ascending node, the argument of perigee, and the mean motion due to the second-order gravitational harmonics of the Earth. (These are the J2 perturbations discussed in Chapter 3.)

The near-circular orbits with moderate inclination and altitudes (cases I and 2) show improvement by including the perturbation. The highly elliptical orbit (case 4) and the low-altitude, high-inclination orbit (case 3) show no pronounced effect from including the perturbation (slight degradation in case 3 and slight improvement in case 4) and the absolute errors are large. The study also indicated

Table 5-11. Comparison of Simple Two-Body Orbit Generators. (From Shear ; see text for discussion.)

0 0