Numerical Computation of Differential-Algebraic Equations for the Approximation of Artificial Satellite Trajectories and Planetary Ephemerides

摘要

The principle of virtual work and Lagrange's equations of motion are used to construct a system of differential equations for constrained spatial multibody system modeling. The differential equations are augmented with algebraic constraints representing the system being modeled. The resulting system is a high index differential-algebraic equation (DAE) which is cast as an ordinary differential equation (ODE) by differentiating the constraint equations twice. The initial conditions are the heliocentric rectangular equatorial generalized coordinates and their first time derivatives of the planets of the solar system and an artificial satellite. The ODE is computed using the integration subroutine LSODAR to generate the body generalized coordinates and time derivatives and hence produce the planetary ephemerides and satellite trajectories for a time interval. Computer simulation and graphical output indicate the satellite and planetary positions and the latter may be compared with those provided in the Astronomical Almanac. Constraint compliance is investigated to establish the accuracy of the computation.