On the existence of the relative equilibria of a rigid body in the J 2 problem

The motion of a point mass in the J ₂ problem has been generalized to that of a rigid body in a J ₂ gravity field for new high-precision applications in the celestial mechanics and astrodynamics. Unlike the original J ₂ problem, the gravitational orbit-rotation coupling of the rigid body is considered in the generalized problem. The existence and properties of both the classical and non-classical relative equilibria of the rigid body are investigated in more details in the present paper based on our previous results. We nondimensionalize the system by the characteristic time and length to make the study more general. Through the study, it is found that the classical relative equilibria can always exist in the real physical situation. Numerical results suggest that the non-classical relative equilibria only can exist in the case of a negative J ₂, i.e., the central body is elongated; they cannot exist in the case of a positive J ₂ when the central body is oblate. In the case of a negative J ₂, the effect of the orbit-rotation coupling of the rigid body on the existence of the non-classical relative equilibria can be positive or negative, which depends on the values of J ₂ and the angular velocity Ω _e . The bifurcation from the classical relative equilibria, at which the non-classical relative equilibria appear, has been shown with different parameters of the system. Our results here have given more details of the relative equilibria than our previous paper, in which the existence conditions of the relative equilibria are derived and primarily studied. Our results have also extended the previous results on the relative equilibria of a rigid body in a central gravity field by taking into account the oblateness of the central body.     

Bifurcations of relative equilibria for one spheroidal and two spherical bodies

We discuss existence and bifurcations of non-collinear (Lagrangian) relative equilibria in a generalized three body problem. Specifically, one of the bodies is a spheroid (oblate or prolate) with its equatorial plane coincident with the plane of motion where only the “J ₂” term from its potential expansion is retained. We describe the bifurcations of relative equilibria as function of two parameters: J ₂ and the angular velocity of the system formed by the mass centers. We offer the values of the parameters where bifurcations in shape occur and discuss their physical meaning. We conclude with a general theorem on the number and the shape of relative equilibria.     

Equilibrium points in the photogravitational restricted four-body problem

We study numerically the photogravitational version of the problem of four bodies, where an infinitesimal particle is moving under the Newtonian gravitational attraction of three bodies which are finite, moving in circles around their center of mass fixed at the origin of the coordinate system, according to the solution of Lagrange where they are always at the vertices of an equilateral triangle. The fourth body does not affect the motion of the three bodies (primaries). We consider that the primary body m ₁ is dominant and is a source of radiation while the other two small primaries m ₂ and m ₃ are equal. In this case (photogravitational) we examine the linear stability of the Lagrange triangle solution. The allowed regions of motion as determined by the zero-velocity surface and corresponding equipotential curves, as well as the positions of the equilibrium points on the orbital plane are given. The existence and the number of the collinear and the non-collinear equilibrium points of the problem depends on the mass parameters of the primaries and the radiation factor q ₁. Critical masses m ₃ and radiation q ₁ associated with the existence and the number of the equilibrium points are given. The stability of the relative equilibrium solutions in all cases are also studied. In the last section we investigate the existence and location of the out of orbital plane equilibrium points of the problem. We found that such critical points exist. These points lie in the (x,z) plane in symmetrical positions with respect to (x,y) plane. The stability of these points are also examined.     

Nonlinear analytic solution for perturbed relative motion using differential equinoctial elements

This paper develops a nonlinear analytic solution for satellite relative motion in J₂-perturbed elliptic orbits by using the geometric method that can avoid directly solving the complex differential equations. The differential equinoctial elements (DEEs) are used to remove any singularities for zero-eccentricity or zero-inclination orbits. Based on the relationship between the relative states and the DEEs, state transition tensors (STTs) for transforming the osculating DEEs and propagating the mean DEEs have been derived. The formulation of these STTs has been split into a set of vector and matrix operations, which avoids directly expanding the complex second-order terms, and thus, the obtained STTs could be easy-to-understand and easy-to-code. Numerical results show that the proposed nonlinear solution is valid for zero-eccentricity and zero-inclination reference orbit and is more accurate than the previous linear or nonlinear methods for the long-term prediction of satellite relative motion.     

Effects of axis-symmetric body,zonal harmonic and potential from a belt on the stability of triangular libration points in the CR3BP

Within the frame work of the circular restricted three-body problem (CR3BP) we have examined the effect of axis-symmetric of the bigger primary, oblateness up to the zonal harmonic J ₄ of the smaller primary and gravitational potential from a belt (circular cluster of material points) on the linear stability of the triangular libration points. It is found that the positions of triangular libration points and their linear stability are affected by axis-symmetric of the bigger primary, oblateness up to J ₄ of the smaller primary and the potential created by the belt. The axis-symmetric of the bigger primary and the coefficient J ₂ of the smaller primary have destabilizing tendency, while the coefficient J ₄ of the smaller primary and the potential from the belt have stabilizing tendency. The overall effect of these perturbations has destabilizing tendency. This study can be useful in the investigation of motion of a particle near axis-symmetric—oblate bodies surrounded by a belt.     

Hamiltonian dynamics of a rigid body in a central gravitational field

This paper concerns the dynamics of a rigid body of finite extent moving under the influence of a central gravitational field. A principal motivation behind this paper is to reveal the hamiltonian structure of the n-body problem for masses of finite extent and to understand the approximation inherent to modeling the system as the motion of point masses. To this end, explicit account is taken of effects arising because of the finite extent of the moving body. In the spirit of Arnold and Smale, exact models of spin-orbit coupling are formulated, with particular attention given to the underlying Lie group framework. Hamiltonian structures associated with such models are carefully constructed and shown to benon-canonical. Special motions, namely relative equilibria, are investigated in detail and the notion of anon-great circle relative equilibrium is introduced. Non-great circle motions cannot arise in the point mass model. In our analysis, a variational characterization of relative equilibria is found to be very useful. Thereduced hamiltonian formulation introduced in this paper suggests a systematic approach to approximation of the underlying dynamics based on series expansion of the reduced hamiltonian. The latter part of the paper is concerned with rigorous derivations of nonlinear stability results for certain families of relative equilibria. Here Arnold's energy-Casimir method and Lagrange multiplier methods prove useful. This work was supported in part by the AFOSR University Research Initiative Program under grant AFOSR-87-0073, by AFOSR grant 89-0376, and by the National Science Foundation's Engineering Research Centers Program: NSFD CDR 8803012. The work of P.S. Krishnaprasad was also supported by the Army Research Office through the Mathematical Sciences Institute of Cornell University.     

New formulation of the two body problem using a continued fractional potential

In the present work, the two body problem using a potential of a continued fractions procedure is reformulated. The equations of motion for two bodies moving under their mutual gravity is constructed. The integrals of motion, angular momentum integral, center of mass integral, total mechanical energy integral are obtained. New orbit equation is obtained. Some special cases are followed directly. Some graphical illustrations are shown. The only included constant of the continued fraction procedure is adjusted so as to represent the so called J ₂ perturbation term of the Earth’s potential.     

The stability of periodic orbits in the three-body problem

A method is developed to study the stability of periodic motions of the three-body problem in a rotating frame of reference, based on the notion of surface of section. The method is linear and involves the computation of a 4×4 variational matrix by integrating numerically the differential equations for time intervals of the order of a period. Several properties of this matrix are proved and also it is shown that for a symmetric periodic motion it can be computed by integrating for half the period only.This linear stability analysis is used to study the stability of a family of periodic motions of three bodies with equal masses, in a rotating frame of reference. This family represents motion such that two bodies revolve around each other and the third body revolves around this binary system in the same direction to a distance which varies along the members of the family. It was found that a large part of the family, corresponding to the case where the distance of the third body from the binary system is larger than the dimensions of the binary system, represents stable motion. The nonlinear effects to the linear stability analysis are studied by computing the intersections of several perturbed orbits with the surface of sectiony ₃=0. In some cases more than 1000 intersections are computed. These numerical results indicate that linear stability implies stability to all orders, and this is true for quite large perturbations.     

Hamiltonian Dynamics of a Gyrostat in the N-Body Problem: Relative Equilibria

We consider the non-canonical Hamiltonian dynamics of a gyrostat in Newtonian interaction with n spherical rigid bodies. Using the symmetries of the system we carry out two reductions. Then, working in the reduced problem, we obtain the equations of motion, a Casimir function of the system and the equations that determine the relative equilibria. Global conditions for existence of relative equilibria are given. Besides, we give the variational characterization of these equilibria and three invariant manifolds of the problem; being calculated the equations of motion in these manifolds, which are described by means of a canonical Hamiltonian system. We give some Eulerian and Lagrangian equilibria for the four body problem with a gyrostat. Finally, certain classical problems of Celestial Mechanics are generalized.     

The equilibrium configurations of the restricted problem of 2+2 triaxial rigid bodies

The restricted 2+2 body problem is considered. The infinitesimal masses are replaced by triaxial rigid bodies and the equations of motion are derived in Lagrange form. Subsequently, the equilibrium solutions for the rotational and translational motion of the bodies are detected. These solutions are conveniently classified in groups according to the several combinations which are possible between the translational equilibria and the constant orientations of the bodies.     

The motion of artificial satellites in the set of Eulerian redundant parameters

In this paper, the connections between orbit dynamics and rigid body dynamics are established throughout the Eulerian redundant parameters, the perturbation equations for any conic motion of artificial satellites are derived in terms of these parameters. A general recursive and stable computational algorithm is also established for the initial-value problem of the Eulerian parameters for satellites prediction in the Earth's gravitational field with axial symmetry. Applications of the algorithm are considered for the two cases of short and long term predictions. For the short-term prediction, we consider the problem of the final state prediction of some typical ballistic missiles in the geopotential model with zonal harmonic terms up to J ₃₆, while for the long-term prediction, we consider the perturbed J ₂ motion of Explorer 28 over 100 revolutions.     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

Satellite orbital precessions caused by the first odd zonal J 3 multipole of a non-spherical body arbitrarily oriented in space

An astronomical body of mass M and radius R which is non-spherically symmetric generates a free space potential U which can be expanded in multipoles. As such, the trajectory of a test particle orbiting it is not a Keplerian ellipse fixed in the inertial space. The zonal harmonic coefficients J ₂,J ₃,… of the multipolar expansion of the potential cause cumulative orbital perturbations which can be either harmonic or secular over time scales larger than the unperturbed Keplerian orbital period T. Here, I calculate the averaged rates of change of the osculating Keplerian orbital elements due to the odd zonal harmonic J ₃ by assuming an arbitrary orientation of the body’s spin axis \(\hat{\boldsymbol{k}}\) . I use the Lagrange planetary equations, and I make a first-order calculation in J ₃. I do not make a-priori assumptions concerning the eccentricity e and the inclination i of the satellite’s orbit.     

Frozen orbits at high eccentricity and inclination: application to Mercury orbiter

We hereby study the stability of a massless probe orbiting around an oblate central body (planet or planetary satellite) perturbed by a third body, assumed to lay in the equatorial plane (Sun or Jupiter for example) using a Hamiltonian formalism. We are able to determine, in the parameters space, the location of the frozen orbits, namely orbits whose orbital elements remain constant on average, to characterize their stability/unstability and to compute the periods of the equilibria. The proposed theory is general enough, to be applied to a wide range of probes around planet or natural planetary satellites. The BepiColombo mission is used to motivate our analysis and to provide specific numerical data to check our analytical results. Finally, we also bring to the light that the coefficient J ₂ is able to protect against the increasing of the eccentricity due to the Kozai-Lidov effect and the coefficient J ₃ determines a shift of the equilibria.     

Effects of zonal harmonics and a circular cluster of material points on the stability of triangular equilibrium points in the R3BP

We have studied a modified version of the classical restricted three-body problem (CR3BP) where both primaries are considered as oblate spheroids and are surrounded by a homogeneous circular planar cluster of material points centered at the mass center of the system. In this dynamical model we have examined the effects of oblateness of both primaries up to zonal harmonic J ₄; together with gravitational potential from the circular cluster of material points on the existence and linear stability of the triangular equilibrium points. It is found that, the triangular points are stable for 0<μ<μ _c and unstable for $\mu_{c} \le \mu \le \frac{1}{2}$ , where μ _c is the critical mass ratio affected by the oblateness up to J ₄ of the primaries and potential from the circular cluster of material points. The coefficient J ₄ has stabilizing tendency, while J ₂ and the potential from the circular cluster of material points have destabilizing tendency. A practical application of this model could be the study of the motion of a dust particle near oblate bodies surrounded by a circular cluster of material points.     

Combined effect of oblateness,radiation and a circular cluster of material points on the stability of triangular liberation points in the R3BP

This paper studies the motion of an infinitesimal mass in the framework of the restricted three-body problem (R3BP) under the assumption that the primaries of the system are radiating-oblate spheroids, enclosed by a circular cluster of material points. It examines the effects of radiation and oblateness up to J ₄ of the primaries and the potential created by the circular cluster, on the linear stability of the liberation locations of the infinitesimal mass. The liberation points are found to be stable for 0<μ<μ _c and unstable for $\mu_{c}\le\mu\le\frac{1}{2}$ , where μ _c is the critical mass value depending on terms which involve parameters that characterize the oblateness, radiation forces and the circular cluster of material points. The oblateness up to J ₄ of the primaries and the gravitational potential from the circular cluster of material points have stabilizing propensities, while the radiation of the primaries and the oblateness up to J ₂ of the primaries have destabilizing tendencies. The combined effect of these perturbations on the stability of the triangular liberation points is that, it has stabilizing propensity.     

Integrable approximation of J 2-perturbed relative orbits

Most existing satellite relative motion theories utilize mean elements, and therefore cannot be used for calculating long-term bounded perturbed relative orbits. The goal of the current paper is to find an integrable approximation for the relative motion problem under the J ₂ perturbation, which is adequate for long-term prediction of bounded relative orbits with arbitrary inclinations. To that end, a radial intermediary Hamiltonian is utilized. The intermediary Hamiltonian retains the original structure of the full J ₂ Hamiltonian, excluding the latitude dependence. This formalism provides integrability via separation, a fact that is utilized for finding periodic relative orbits in a local-vertical local-horizontal frame and determine an initialization scheme that yields long-term boundedness of the relative distance. Numerical experiments show that the intermediary-based computation of orbits provides long-term bounded orbits in the full J ₂ problem for various inclinations. In addition, a test case is shown in which the radial intermediary-based initial conditions of the chief and deputy satellites yield bounded relative distance in a high-precision orbit propagator.     

Stability of Triangular Equilibrium Points in the Photogravitational Restricted Three-Body Problem with Oblateness and Potential from a Belt

We have examined the effects of oblateness up to J ₄ of the less massive primary and gravitational potential from a circum-binary belt on the linear stability of triangular equilibrium points in the circular restricted three-body problem, when the more massive primary emits electromagnetic radiation impinging on the other bodies of the system. Using analytical and numerical methods, we have found the triangular equilibrium points and examined their linear stability. The triangular equilibrium points move towards the line joining the primaries in the presence of any of these perturbations, except in the presence of oblateness up to J ₄ where the points move away from the line joining the primaries. It is observed that the triangular points are stable for 0 < μ < μ _c and unstable for \(\mu_{\mathrm{c}} \le \mu \le \frac {1}{2},\) where μ _c is the critical mass ratio affected by the oblateness up to J ₄ of the less massive primary, electromagnetic radiation of the more massive primary and potential from the belt, all of which have destabilizing tendencies, except the coefficient J₄ and the potential from the belt. A practical application of this model could be the study of motion of a dust particle near a radiating star and an oblate body surrounded by a belt.     

A highly symmetric restricted nine-body problem and the linear stability of its relative equilibria

We study a highly symmetric nine-body problem in which eight positive masses, called the primaries, move four by four, in two concentric circular motions such that their configuration is always a square for each group of four masses. The ninth body being of negligible mass and not influencing the motion of the eight primaries. We assume all the nine masses are in the same plane and that the masses of the primaries are \(m_{1}=m_{2}=m_{3}=m_{4}=\tilde{m}\) and m ₅=m ₆=m ₇=m ₈=m and the radii associated to the circular motion of the bodies with mass \(\tilde{m}\) is λ∈[λ ₀,1] and for the bodies with mass m is 1. We prove the existence of central configurations which characterize such arrangement of the primaries and we study the influence of the parameter λ, the ratio of the radii of the two circles, on the masses m and \(\tilde{m}\) . We use a synodical system of coordinates to eliminate the time dependence on the equations of motion. We show the existence of equilibria solutions symmetrically distributed on the four quadrants and their dependence on the parameter λ. Finally, we show that there can be 13, 17 or 25 equilibria solutions depending on the size of λ and we investigate their linear stability.