Stability of L4 for radiated rigid primaries in resonant cases

The non-linear stability of the triangular libration point L₄ of the restricted three-body problem is studied under the presence of third- and fourth-order resonances, when the more massive primary is a triaxial rigid body and source of radiation. In this study, Markeev's theorems are applied with the help of Moser's theorem. It is found that the stability of the triangular libration point is unstable in the third-order resonance case and in the fourth-order resonance case, this is stable or unstable depending on A₁ and A₂, and a source of radiation parameter α, where A₁, A₂ depend upon the lengths of the semi-axes of the triaxial rigid body.     

Attitude dynamics of a rigid body on a Keplerian orbit: A simplification

An infinitestimal contact transformation is proposed to simplify at first order the Hamiltonian representing the attitude of a triaxial rigid body on a Keplerian orbit around a mass point. The simplified problem reduces to the Euler-Poinsot model, but with moments of inertia depending on time through the longitude in orbit. Should the orbit be circular, the moments of inertia would be constant.     

Dynamics of a triaxial gyrostat at a Lagrangian equilibrium of a binary asteroid

The non-canonical Hamiltonian dynamics of a triaxial gyrostat in Newtonian interaction with two punctual masses is considered. This serves as a model for the study of the attitude dynamics of a spacecraft located at a Lagrangian equilibrium point of the system formed by a binary asteroid and a spacecraft. Using geometric-mechanics methods, the approximated dynamics that arises when developing the potential in series of Legendre functions and truncating the series to the second harmonics is studied. Working in the reduced problem, the existence of equilibria in Lagrangian form are studied, in analogy with classic results on the topic. In this way, the classical results on equilibria of the three-body problem, as well as other results by different authors that use more conventional techniques for the case of rigid bodies, are generalized. The rotational Poisson dynamics of a spacecraft located at a Lagrangian equilibrium and the study of the nonlinear stability of some important equilibria are considered. The analysis is done in vectorial form avoiding the use of canonical variables and the tedious expressions associated with them.     

Out-of-plane equilibrium points and their stability in the Robe’s problem with oblateness and triaxiality

This paper examines the existence and stability of the out-of-plane equilibrium points of a third body of infinitesimal mass when the equations of motion are written in the three dimensional form under the set up of the Robe’s circular restricted three-body problem, in which the hydrostatic equilibrium figure of the first primary is an oblate spheroid and the second one is a triaxial rigid body under the full buoyancy force of the fluid. The existence of the out of orbital plane equilibrium points lying on the xz-plane is noticed. These points are however unstable in the linear sense.     

Fast computation of Jacobian elliptic functions and incomplete elliptic integrals for constant values of elliptic parameter and elliptic characteristic

In order to accelerate the numerical evaluation of torque-free rotation of triaxial rigid bodies, we present a fast method to compute various kinds of elliptic functions for a series of the elliptic argument when the elliptic parameter and the elliptic characteristic are fixed. The functions we evaluate are the Jacobian elliptic functions and the incomplete elliptic integral of the second and third kinds regarded as a function of that of the first kind. The key technique is the utilization of the Maclaurin series expansion and the addition theorems with respect to the elliptic argument. The new method is around 25 times faster than the method using the incomplete elliptic integral of general kind and around 70 times faster than the method using mathematical libraries given in the latest version of Numerical Recipes.     

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.     

A study of non-collinear libration points under the combined effects of radiation pressure and Stokes drag

We study the non-collinear libration points in the frame work of photo-gravitational circular restricted three-body problem with Stokes drag acting as a dissipative force and considering the more massive primary as a radiating body and the less massive primary as a triaxial rigid body. The combined effects of radiation pressure and Stokes drag on the existence and stability of non-collinear libration points is analyzed. It is found that there exist two non-collinear libration points and are asymptotically stable in the interval 0.6149 ≤ q ≤ 1 for μ = 0.01, where q and μ are the radiation factor and mass ratio, respectively.     

On the free rotation of a rigid body

It is well known that the equations governing the motion of a freely-rotating rigid body possess an exact analytical solution, involving Jacobi's elliptic functions. Andoyer (1923) and Deprit (1967) have shown that the problem may be very usefully reduced to a one-degree-of-freedom Hamiltonian system. When two of the body's principal moments of inertia are very nearly equal, the Hamiltonian system has the same form as the Ideal Resonance Problem. In earlier publications (Jupp, 1969, 1972, 1973), the author has constructed formal power-series solutions of the latter problem.In this article, the general solution of the Ideal Resonance Problem is employed to formulate a second-order formal series solution of the problem of a freely-rotating rigid body which has two of its principal moments of inertia differing by a small quantity. This solution is firstly expressed in terms of the mean elements, and then in terms of the initial conditions. The latter solution is global in nature being applicable over the whole phase plane. It is demonstrated that the exact solution and the second-order formal series solution, written in terms of the initial conditions, differ by terms of at most third order in the small parameter, over the whole domain of possible motions. This serves as an important check on the general results published in the earlier articles.     

The Non-Linear Stability of L4 in the Restricted Three-Body Problem when the Bigger Primary is a Triaxial Rigid Body

The non-linear stability of L ₄ in the restricted three-body problem has been studied when the bigger primary is a triaxial rigid body with its equatorial plane coincident with the plane of motion. It is found that L ₄ is stable in the range of linear stability except for three mass ratios:

Stability of triangular points in the photogravitational CR3BP with Poynting-Robertson drag and a smaller triaxial primary

The stability of triangular equilibrium points in the framework of the circular restricted three-body problem (CR3BP) is investigated for a test particle of infinitesimal mass in the vicinity of two massive bodies (primaries), when the bigger primary is a source of radiation and the smaller one is a triaxial rigid body with one of the axes as the axis of symmetry and its equatorial plane coinciding with the plane of motion, under the Poynting-Robertson (P-R) drag effect as a result of the radiating primary. It is found that the involved parameters influence the position of triangular points and their linear stability. It is noted that these points are unstable in the presence of Poynting-Robertson drag effect and conditionally stable in the absence of it.     

Effects of radiation and triaxiality of primaries on triangular equilibrium points in elliptic restricted three body problem

This paper studies the motion of an infinitesimal mass around triangular equilibrium points in the elliptic restricted three body problem assuming bigger primary as a source of radiation and the smaller one a triaxial rigid body. A practical application of this case could be the study of motion of a satellite under the effect of Sun and Earth. We have exploited the method of averaging used by Grebnikov (Nauka, Moscow, revised 1986) throughout the analysis of stability of the system. The critical mass ratio depends on the radiation pressure, oblateness, eccentricity and semi major axis of the elliptic orbits and the range of stability decreases as the radiation parameter increases.     

On the restricted three-body problem with generalized forces

By introducing general functions which depend on distance, a general scheme which determines the equilibrium solutions for the generalized restricted three-body problem is given. Applications to problems such as primaries considered as rigid bodies, influence of the radiation pressure of the primaries, and a combination of radiation pressure and rigid body are presented.     

Green's functions of the induction equation on regions with boundary. III. Half-space

In two earlier papers (BRÄUER and RÄDLER 1986, 1987) the evolution of a magnetic field was considered which pervades an electrically conducting fluid and its non-conducting surroundings. A construction principle for Green's functions of the corresponding initial value problem was proposed, and worked out for the case in which the fluid fills a spherical region. Now the principle is applied to the case of a fluid body occupying a half-space. Green's functions are constructed for arbitrary motions of the fluid. More concrete results are derived for shear flow, and explicit expressions of Green's functions are given for rigid body motion.     

Effects of the Triaxiality on the Rotation of Celestial Bodies: Application to the Earth,Mars and Eros

In this paper we discuss the influence of the triaxiality of a celestialbody on its free rotation, i.e. in absence of any external gravitationalperturbation. We compare the results obtained through two different analytical formalisms, one established from Andoyer variables by usingHamiltonian theory, the other one from Euler's variables by usingLagrangian equations. We also give a very accurate formulation of thepolar motion (polhody) in the case of a small amplitude of this motion.Then, we carry out a numerical integration of the problem, with aRunge–Kutta–Felberg algorithm, and for the two kinds of methods above, that we apply to three different celestial bodies considered as rigid : the Earth, Mars, and Eros. The reason of this choice is that each of this body corresponds to a more or less triaxial shape.In the case of the Earth and Mars we show the good agreement betweenanalytical and numerical determinations of the polar motion, and theamplitude of the effect related to the triaxial shape of the body, whichis far from being negligible, with some influence on the polhody of theorder of 10 cm for the Earth, and 1 m for Mars. In the case of Eros, weuse recent output data given by the NEAR probe, to determine in detailthe nature of its free rotational motion, characterized by the presence ofimportant oscillations for the Euler angles due to the particularly largetriaxial shape of the asteroid.     

Robe’s Circular Restricted Three-Body Problem Under Oblate and Triaxial Primaries

This paper analyzes Robe??s circular restricted three-body problem when the hydrostatic equilibrium figure of the first primary is assumed to be an oblate spheroid, the shape of the second primary is considered as a triaxial rigid body, and the full buoyancy force of the fluid is taken into account. It is found that there is an equilibrium point near the center of the first primary, another equilibrium point exists on the line joining the centers of the primaries and there exist infinite number of equilibrium points on an ellipse in the orbital plane of the second primary. It is also observed that under certain conditions, all these equilibrium points can be stable. The most interesting and distinguishable results of this study are the existence of elliptical points and their stability.     

Stability of the classical type of relative equilibria of a rigid body in the J 2 problem

The motion of a point mass in the J ₂ problem is generalized to that of a rigid body in a J ₂ gravity field. The linear and nonlinear stability of the classical type of relative equilibria of the rigid body, which have been obtained in our previous paper, are studied in the framework of geometric mechanics with the second-order gravitational potential. Non-canonical Hamiltonian structure of the problem, i.e., Poisson tensor, Casimir functions and equations of motion, are obtained through a Poisson reduction process by means of the symmetry of the problem. The linear system matrix at the relative equilibria is given through the multiplication of the Poisson tensor and Hessian matrix of the variational Lagrangian. Based on the characteristic equation of the linear system matrix, the conditions of linear stability of the relative equilibria are obtained. The conditions of nonlinear stability of the relative equilibria are derived with the energy-Casimir method through the projected Hessian matrix of the variational Lagrangian. With the stability conditions obtained, both the linear and nonlinear stability of the relative equilibria are investigated in details in a wide range of the parameters of the gravity field and the rigid body. We find that both the zonal harmonic J ₂ and the characteristic dimension of the rigid body have significant effects on the linear and nonlinear stability. Similar to the classical attitude stability in a central gravity field, the linear stability region is also consisted of two regions that are analogues of the Lagrange region and the DeBra-Delp region respectively. The nonlinear stability region is the subset of the linear stability region in the first quadrant that is the analogue of the Lagrange region. Our results are very useful for the studies on the motion of natural satellites in our solar system.     

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 libration points of axisymmetric satellite in the gravitation field of two triaxial rigid bodies

In this paper we consider the restricted problem of three rigid bodies (an axisymmetric satellite in the gravitation field of two triaxial primaries). The collinear and triangular equilibrium solutions are obtained. The effect of the primaries on the location of the libration points of a spherical satellite has been studied numerically.     

Exact solution of a triaxial gyrostat with one rotor

The problem of the attitude dynamics of a triaxial gyrostat under no external torques and one constant internal rotor, is a three degrees-of-freedom system, although thanks to the existence of integrals of motion it can be reduced to only one degree-of-freedom problem. We introduce coordinates to represent the orbits of constant angular momentum as a flow on a sphere. This representation shows that the problem is equivalent to a quadratic Hamiltonian depending on two parameters. We find the exact solution of the orbits in terms of elliptic functions. By making use of properties of elliptic functions we find the solution at each region of the parametric partition from the solution of one region. We also prove that heteroclinic orbits are planar curves.