The fifth-order analytical solution of the equations of motion of a satellite in orbit around a non-spherical planet

A high-precise analytical theory of a satellite in orbit around a non-spherical planet has been developed. The Poisson's small parameter method has been used. All secular and short-periodic perturbations proportional up to and including a product of five arbitrary harmonic coefficients of the planetary potential expansion are calculated. Long-periodic perturbations are derived with the accuracy of up to the fourth-order, inclusive. The influence of the high-order perturbations on the motion of ETALON-1 satellite has been investigated. The results of comparison of the numerical and analytical integration of the equations of its motion over a five year interval are as follows:

The equations of motion of an artificial satellite in nonsingular variables

The equations of motion of an artificial satellite are given in nonsingular variables. Any term in the geopotential is considered as well as luni-solar perturbations up to an arbitrary power ofr/r, r being the geocentric distance of the disturbing body. Resonances with tesseral harmonics and with the Moon or Sun are also considered. By neglecting the shadow effect, the disturbing function for solar radiation is also developed in nonsingular variables for the long periodic perturbations. Formulas are developed for implementation of the theory in actual computations.     

Quadrature solution for the general relativistic motion of a satellite or a planet

The Schwarzschild field of a central massM is used to derive the general relativistic motion of a particle in a bounded orbit aroundM. A quadrature gives the central angle as a quasi-periodic function (f) of an effective true anomalyf. The linear term is an infinite series, whose second term yields the usual rate of advance of pericenter. For an artificial satellite this may be as large as 17 of arc per year. The periodic part is a sine series, with coefficients containing the small parameter 2GM/c ² p, wherep is closely approximated by the classical semi-latus rectum. The radius vectorr is a Kepler-like function off.The essentially new features of the calculation are the appropriate factoring of a certain cubic polynomialF(p/r), the use of the above effective true anomalyf, and the introduction of an effective eccentric anomalyE. The latter serves to reduce the differential equation forf as a function of the timet, obtained by combining the solution for (f) with the relativistic integrals of motion, to a Kepler equation forE.Knowing the constants of the motion, one can then solve successively forE(t), f(t), r(t), and (t). This is best done as a variational calculation, comparing the relativistic orbiter with a classical orbiter having the same initial conditions. The resulting variations agree with those of Lass and Solloway, but the present method is quite different from theirs and results in a simpler algorithm. The results show that the radial and transverse corrections, r andr , arising from the Schwarzschild field, may be of the order of a kilometer for 1000 revolutions of an Earth satellite of orbital eccentricitye ₀0.6.For bounded motion, the cubic polynomialF(p/r) has three positive real zeros, the two smaller ones corresponding to apocenter and pericenter. The third and apparently non-physical one occurs forrSchwarzschild radius. It may thus correspond to the incipient fall of the orbiter into the central body, if the latter is a black hole.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 27–September 2, 1972.Research sponsored by NASA Goddard Space Flight Center under Contract No. NAS 5-11909.     

An alternate transition from the Lagrangian of a satellite to equations of motion

A procedure is outlined for a conservative system which makes it possible to go from a Lagrangian of a librating system to the corresponding equations of motion in the Eulerian form. The transition does not require a choice of rotational coordinates and makes use of angular velocities and direction cosines directly. The procedure thus synthesizes attractive features of two classical approaches and has far reaching consequences: it is particularly useful in formulating equations of motion for complex flexible systems of contemporary interest. For the case of a satellite with two flexible plate-type appendages, for example, the approach reduced the formulation time to one-third. The basic mathematical concepts are briefly touched upon in the beginning which help explain the subsequent development.     

On an alternative method for approximate solution of linear stochastic differential equations

An approximate method for solving formal linear stochastic differential equations of first order is proposed. On the basis of the Reynold averaging technique, the stochastic differential equation is transformed into an infinite hierarchical system. This infinite system is cut off in such a manner that uncorrelated (totally randomly) and totally correlated stochastic processes are exactly included. By this the applicability of the Reynold method is extended.     

On the roto-translatory motion of a satellite of an oblate primary

We deal with the problem of the motion of a triaxial satellite of an oblate primary of larger mass. We show that the treatment is simplified by using a canonical set of variables previously introduced by the authors, that allows a drastic reduction in the expansions of the potential. A general method to avoid the appearance of virtual singularities when the angles between certain planes are small is designed. Our approach is applicable either to natural or artificial satellites.     

Computerized series solution of relativistic equations of motion

A method of solution of the equations of planetary motion is described. It consists of the use of numerical general perturbations in orbital elements and in rectangular coordinates. The solution is expanded in Fourier series in the mean anomaly with the aid of harmonic analysis and computerized series manipulation techniques. A detailed application to the relativistic motion of the planet Mercury is described both for Schwarzschild and isotropic coordinates.Receipt delayed by the postal strike in Great Britain.     

A semi-analytic theory for the motion of a lunar satellite

A semi-analytical solution to the problem of the motion of a satellite of the moon is presented. Perturbative effects which are considered include those due to the attraction of the moon, earth, and sun, the non-sphericity of the moon's gravitational field, coupling of lower-order terms, solar radiation pressure, and physical libration. Short-period terms and intermediate-period terms, terms with the period of the moon's longitude, are produced by means of von Zeipel's method; it is proposed to obtain the secular perturbations, and those depending only on the argument of perilune, by numerical integration of the equations of motions. The short-period terms and intermediate-period terms are developed up to second order, where first order is 10^–2. The secular perturbations and perturbations dependent on the argument of perilune are obtained to third order.     

The motion of a satellite in an axi-symmetric gravitational field

The equations for the variation of the osculating elements of a satellite moving in an axi-symmetric gravitational field are integrated to yield the complete first-order perturbations for the elements of the orbit. The expressions obtained include the effects produced by the second to eighth spherical harmonics. The orbital elements are presented in the most general form of summations by means of Hansen coefficients. Due to their general forms it is a simple matter to estimate the perturbations of any higher harmonic by simply increasing the index of summation. Finally, this paper gives the respective general expressions for the secular perturbations of the orbital elements. The formulae presented should be useful for the reductions of Earth-satellite observations and geopotential studies based on them.List of Symbols semi-major axis - C _jk ⁿ (, ) cosine functions of and - e eccentricity of the orbit - f acceleration vector of perturbing force - f sin²t - i inclination of the orbit - J _n coefficients in the potential expansion - M mean anomaly - n mean motion - p semi-latus rectum of the orbit - R, S, andW components of the perturbing acceleration - r radius-vector of satellite - r magnitude ofr - S _jk ⁿ (, ) sine functions of and - T time of perigee passage - u argument of latitude - U gravitational potential - true anomaly - V perturbing potential - G(M₊+m) (gravitational constant times the sum of the masses of Earth and satellite) - _n,k coefficients ofR component of disturbing acceleration (funtions off) - _n,k coefficients ofS andW components of disturbing acceleration (functions off) - mean anomaly at timet=0 - X ₀ ^n,m zero-order Hansen coefficients - argument of perigee - right ascension of the ascending node     

Theory of the motion of an artificial Earth satellite

An improved analytical solution is obtained for the motion of an artificial Earth satellite under the combined influences of gravity and atmospheric drag. The gravitational model includes zonal harmonics throughJ ₄, and the atmospheric model assumes a nonrotating spherical power density function. The differential equations are developed through second order under the assumption that the second zonal harmonic and the drag coefficient are both first-order terms, while the remaining zonal harmonics are of second order.Canonical transformations and the method of averaging are used to obtain transformations of variables which significantly simplify the transformed differential equations. A solution for these transformed equations is found; and this solution, in conjunction with the transformations cited above, gives equations for computing the six osculating orbital elements which describe the orbital motion of the satellite. The solution is valid for all eccentricities greater than 0 and less than 0.1 and all inclinations not near 0^o or the critical inclination. Approximately ninety percent of the satellites currently in orbit satisfy all these restrictions.     

Air drag effect on the motion of an artificial Earth satellite

Two methods have been used to compute and compare the perturbations in perigee distance for an artificial Earth satellite. The two methods have used different air density models. The first (Helali, 1987) used the TD model, formulated by Sehnel (1986a), which contains terms that describe all the principal changes of the thermospheric density due to solar activity, geomagnetic activity, and the height. The second method (Davis, 1963) used a model of the density which takes into account the rotation of the atmosphere, the bulging atmosphere and the height. For different values of eccentricities from 0.001 to 0.05 we computed the perturbations P _r in the perigee distance at different heights from 200 to 350 km for both methods. The results show a good agreement for the computed values of P _r for different values of e (0 < e 0.02) in both methods at perigee heights from 250 to 350 km. Meanwhile, for perigee heights smaller than about 250 km we found a maximum difference in P _r amounting to 20 metres/revolution for e = 0.005 and 0.01.     

Third-order secular Solution of the variational equations of motion of a satellite in orbit around a non-spherical planet

We constructed an analytical theory of satellite motion up to the third order relative to the oblateness parameter of the Earth (J ₂). Equations of secular variations was developed for the first three orbital elements (a, e, i) of an artificial satellite. The secular variations are solved in a closed form.     

The relativistic equations of motion for a satellite in orbit about a finite-size,rotating Earth

The relativistic equations of motion are derived for N self-gravitating, rotating finite bodies. These equations are then applied to the near-Earth satellite orbit determination problem. The apparent change of the shape of the Earth from the Earth centered frame to the Solar System barycentric frame changes the value of the Newtonian potential term in the metric. This in turn leads to a simplification of the equations of motion in the barycentric frame.     

The geomagnetic effects on the motion of an electrically charged artificial satellite

The orbital effects of the Lorentz force on the motion of an electrically charged artificial satellite moving in the Earth's magnetic field are determined. The geomagnetic field is considered as a multipole potential field and the satellite electrical charge is supposed to be constant. The relativistic perturbations of the main geomagnetic field are discussed briefly. The results are concentrated on the determination of the secular changes, and numerical values are computed for the case of the LAGEOS satellite. The results are discussed in the context of a possible detection of the Lense-Thirring effect analyzing the orbital perturbations of the LAGEOS and LAGEOS X satellites.     

The motion of a lunar satellite

Presented in this theory is a semianalytical solution for the problem of the motion of a satellite in orbit around the moon. The principal perturbations on such a body are due to the nonspherical gravity field of the moon, the attraction of the earth, and, to a lesser degree, the attraction of the sun. The major part of the problem is solved by means of the celebrated von Zeipel Method, first successfully applied to the motion of an artificial earth satellite by Brouwer in 1959. After eliminating from the Hamiltonian all terms with the period of the satellite and those with the period of the moon, it is suggested to solve the remaining problem with the aid of numerical integration of the modified equations of motion.This theory was written in 1964 and presented as a dissertation to Yale University in 1965. Since then a great deal has been learned about the gravity field of the moon. It seems that quite a number of recently determined gravity coefficients would qualify as small quantities of order two. Hence, according to the truncation criteria employed, they should be considered in the present theory. However, the author has not endeavored to update the work accordingly. The final results, therefore, are incomplete in the lunar gravitational perturbations. Nevertheless, the theory does give the largest such variations and it does present the methods by which perturbations may be derived for any gravity terms not actually developed.     

On Pfaff's equations of motion in dynamics; Applications to satellite theory

In this article we study a form of equations of motion which is different from Lagrange's and Hamilton's equations: Pfaff's equations of motion. Pfaff's equations of motion were published in 1815 and are remarkably elegant as well as general, but still they are much less well known. Pfaff's equations can also be considered as the Euler-Lagrange equations derived from the linear Lagrangian rather than the usual Lagrangian which is quadratic in the velocity components. The article first treats the theory of changes of variables in Pfaff's equations and the connections with canonical equations as well as canonical transformations. Then the applications to the perturbed two-body problem are treated in detail. Finally, the Pfaffians are given in Hill variables and Scheifele variables. With these two sets of variables, the use of the true anomaly as independent variable is also considered.     

Tidal effects and the motion of a satellite

The effect of resonant planetary perturbations on the evolution of the orbit of a satellite driven by tidal forces is studied in this paper. The basic equations that govern it are similar to the equations found in orbit-orbit and in spin-orbit couplings. The general form of these equations is: A general treatment of such equations, proposed earlier (J. Kovalevsky, in Dynamical Trapping and Evolution of the Solar system, IAU Colloquium n^o74, V. V. Markellos and Y. Kozai, eds., 1983) is sketched.In particular, the effects of the large long periodic variations of the excentricity e' of the planet are analysed on an example taken from the lunar theory and the Earth's general theory due to Bretagnon.The argument of the well known planetary term =18 V-16T due to the tidal friction and quasi-periodic variations due to the presence of e' in the expression of the mean motion of the Moon. Their joint effect, has been to produce in the past resonant situations for this argument that repeated more than 100 times. Every such situation can be treated by equation (1).Numerical integration, using conditions that might have occurred while or similar other arguments were quasi resonant, have produced the following results: (a) In some cases, the argument becomes temporarily resonant. Between the capture to and the escape from the resonance, the semi-major axis undergoes oscillations, but the tidal secular evolution is stopped. (b) In other cases, the argument is not trapped into a resonant conditions, but the semi-major axis undergoes a quick change while d/dt is close to zero.A number of arguments that have been quasi resonant in the past history of the Earth-Moon system has been identified from the Chapront and Chapront-Touzé Lunar Theory. It appears that the phenomena described are frequent features in the evolution of the Lunar orbit.     

Poisson equations of rotational motion for a rigid triaxial body with application to a tumbling artificial satellite

Differential equations are derived for studying the effects of either conservative or nonconservative torques on the attitude motion of a tumbling triaxial rigid satellite. These equations, which are analogous to the Lagrange planetary equations for osculating elements, are then used to study the attitude motions of a rapidly spinning, triaxial, rigid satellite about its center of mass, which, in turn, is constrained to move in an elliptic orbit about an attracting point mass. The only torques considered are the gravity-gradient torques associated with an inverse-square field. The effects of oblateness of the central body on the orbit are included, in that, the apsidal line of the orbit is permitted to rotate at a constant rate while the orbital plane is permitted to precess (either posigrade or retrograde) at a constant rate with constant inclination.A method of averaging is used to obtain an intermediate set of averaged differential equations for the nonresonant, secular behavior of the osculating elements which describe the complete rotational motions of the body about its center of mass. The averaged differential equations are then integrated to obtain long-term secular solutions for the osculating elements. These solutions may be used to predict both the orientation of the body with respect to a nonrotating coordinate system and the motion of the rotational angular momentum about the center of mass. The complete development is valid to first order in (n/w ₀)², wheren is the satellite's orbital mean motion andw ₀ its initial rotational angular speed.     

The regular motion of an axisymmetrical satellite in the field of a triaxial body

This paper investigates the regular motions of an axisymmetrical satellite in the field of Newton's attraction of a triaxial body. Both the orbital and the self rotational motions of the two bodies are taken into consideration. The exact solutions are discussed using Poincaré's method of small parameter. In the decomposition of the force function all the harmonic terms up to the third order are taken into account.The results show the existence of eight solutions. The stability of the new group of solutions is discussed using two methods to get the necessary and sufficient conditions required for the stability of these motions.