A study of evolution of spin-up mode of a satellite in the plane of its orbit

We investigate the mode of spinning up a low-orbit satellite in the plane of its orbit. In this mode the satellite rotates around its principal central axis of the minimum moment of inertia which executes small oscillations with respect to the normal to the orbit plane; the angular velocity of the rotation around this axis several times exceeds the mean orbital motion. Gravitational and restoring aerodynamic moments are taken into account in the satellite’s equations of motion. A small parameter characterizing deviation of the satellite from a dynamically symmetric shape is introduced into the equations. A two-dimensional integral surface of the equations of motion, describing quasi-steady-state rotations of the satellite close to cylindrical precession of the corresponding symmetrical satellite in a gravitational field, has been studied by the method of small parameter and numerically. Such quasi-steady-state rotations are suggested to be considered as unperturbed motions of the satellite in the spin-up mode. Investigation of the integral surface is reduced to numerical solution of a periodic boundary value problem of a certain auxiliary system of differential equations and to calculation of quasi-steady-state rotations by the two-cycle method. A possibility is demonstrated to construct quasi-steady rotations by way of minimization of a special quadratic functional.     

Rapid rotations of a satellite with a cavity filled with viscous fluid under the action of moments of gravity and light pressure forces

Rapid rotational motion of a dynamically asymmetric satellite relative to the center of mass is studied. The satellite has a cavity filled with viscous fluid at low Reynolds numbers, and it moves under the action of moments of gravity and light pressure forces. Orbital motions with an arbitrary eccentricity are supposed to be specified. The system, obtained after averaging over the Euler-Poinsot motion and applying the modified averaging method, is analyzed. The numerical analysis in the general case is performed, and the analytical study in the axial rotation vicinity is carried out. The motion in the specific case of a dynamically symmetric satellite is considered.     

One method of gravitational attitude control of a rotating satellite

The mode of spinning up a low-orbit satellite in the plane of its orbit is studied. In this mode, the satellite rotates around its longitudinal axis (principal central axis of the minimum moment of inertia), which executes small oscillations with respect to the normal to the orbit plane; the angular velocity of the rotation around the longitudinal axis is several tenths of a degree per second. Gravitational and restoring aerodynamic moments were taken into account in the equations of satellite’s motion, as well as a dissipative moment from eddy currents induced in the shell of the satellite by the Earth’s magnetic field. A small parameter characterizing deviation of the satellite from a dynamically symmetric shape and nongravitational external moments are introduced into the equations. A two-dimensional integral surface of the equations of motion, describing quasistationary rotations of the satellite close to cylindrical precession of the corresponding symmetrical satellite in a gravitational field, has been studied by the method of small parameter and numerically. We propose to consider such quasistationary rotations as unperturbed motions of the satellite in the spin-up mode.     

Periodic motions of a satellite-gyrostat relative to its center of mass under the action of gravitational torque

We investigated periodic motions of the axis of symmetry of a model satellite of the Earth, which are similar to the motions of the longitudinal axes of the Mir orbital station in 1999–2001 and the Foton-M3 satellite in 2007. The motions of these spacecraft represented weakly disturbed regular Euler precession with the angular momentum vector of motion relative to the center of mass close to the orbital plane. The direction of this vector during the motion was not practically changed. The model satellite represents an axisymmetric gyrostat with gyrostatic moment directed along the axis of symmetry. The satellite moves in a circular orbit and undergoes the action of the gravitational torque. The motion of the axis of symmetry of this satellite relative to the absolute space is described by fourth-order differential equations with periodic coefficients. The periodic solutions to this system with special symmetry properties are constructed using analytical and numerical methods.     

The motion of a synchronous satellite in a system similar to the earth-moon system

The stationary motions of a synchronous axisymmetric satellite are studied in the field of attraction by the Earth and a third body whose parameters are close to those of the Moon. Equations of motion are written in canonical variables that take into account the resonance character of the problem. The plots characterizing the dependence of the rotation parameters of the satellite relative to the center of mass on the elements of satellite’s translational motion are presented. A picture is given that represents the initial configuration of the system for implementing stationary motions.     

Three-frequency resonances in the problem of fast rotation of an asymmetric body in elliptical orbit

The paper presents a classification of resonance problems generated by the highest-degree term of the perturbing function in the problem of fast rotation of an asymmetric solid body in elliptic orbit in a central gravitational field. Explicit formulas are obtained for Hamiltonians of all Hamiltonian systems that deter-mine motions in the neighborhood of resonances. Basic resonance effects are described.     

Biaxial Mode of Rotation of a Satellite in the Orbit Plane

A mode of motion of a satellite with respect to its center of mass is studied, which is called the biaxial rotation in the orbit plane. In this mode of rotation, an elongated and nearly dynamically symmetric satellite rotates around the longitudinal axis, which, in turn, rotates around the normal to the plane of an orbit; the angular velocity of rotation around the longitudinal axis is several times larger than the orbital angular velocity, deviations of this axis from the orbit plane are small. Such a rotation is convenient in the case when it is required to secure a sufficiently uniform illumination of the satellite's surface by the Sun at a comparatively small angular velocity of the satellite. The investigation consists of the numerical integration of equations of the satellite's motion, which take into account gravitational and restoring aerodynamic moments, as well as the evolution of the orbit. At high orbits, the mode of the biaxial rotation is conserved for an appreciable length of time, and at low orbits it is destroyed due to the impact of the aerodynamic moment. The orbit altitudes and the method of constructing the initial conditions of motion that guarantee a sufficiently prolonged period of existence of this mode are specified.     

A specific case of stability of cylindrical precession of a satellite

In a central Newtonian gravitational field, the motion of a dynamically symmetrical satellite along an elliptical orbit of arbitrary eccentricity is considered. The particular motion of the satellite is known when its axis of symmetry is perpendicular to the orbit plane, and the satellite rotates about this axis with a constant angular velocity (cylindrical precession). A nonlinear analysis of stability of this motion has been performed under the assumption that the geometry of the satellite mass corresponds to a thin plate. At small values of orbit eccentricity e the analysis is analytical, while numerical analysis is used for arbitrary values of e.     

On plane oscillations of a pendulum with variable length suspended on the surface of a planet’s satellite

The problem of planar oscillations of a pendulum with variable length suspended on the Moon’s surface is considered. It is assumed that the Earth and Moon (or, in the general case, a planet and its satellite, or an asteroid and a spacecraft) revolve around the common center of mass in unperturbed elliptical Keplerian orbits. We discuss how the change in length of a pendulum can be used to compensate its oscillations. We wrote equations of motion, indicated a rule for the change in length of a pendulum, at which it has equilibrium positions relative to the coordinate system rotating together with the Moon and Earth. We study the necessary conditions for the stability of these motions. Chaotic dynamics of the pendulum is studied numerically and analytically.     

Variation of the Satellite Orbit Altitude by the Light Pressure Force

The possibility of the uncontrolled increase of the altitude of an almost circular satellite orbit by the force of the light pressure is investigated. The satellite is equipped with a damper and a system of mirrors (solar batteries can serve as such a system). The flight of the satellite takes place in the mode of a single-axis gravitational orientation, the axis of its minimum principal central moment of inertia makes a small angle with the local vertical and the motion of the satellite around this axis constitutes forced oscillations under the impact of the moment of force of the light pressure. The form of the oscillations and the initial orbit are chosen so that the transverse component of the force of the light pressure acting upon the satellite be positive and the semimajor axis of the orbit would continuously increase. As this takes place, the orbit remains almost circular. We investigate the evolution of the orbit over an extended time interval by the method which employs separate integration of the equations of the orbital and rotational motions of the satellite. The method includes outer and inner cycles. The outer cycle involves the numerical integration of the averaged equations of motion of the satellite center of mass. The inner cycle serves to calculate the right-hand sides of these equations. It amounts to constructing an asymptotically stable periodic motion of the satellite in the mode of a single-axis gravitational orientation for current values of the orbit elements and to averaging the equations of the orbital motion along it. It is demonstrated that the monotone increase of the semimajor axis takes place during the first 15 years of motion. In actuality, the semimajor axis oscillates with a period of about 60 years. The eccentricity and inclination of the orbit remain close to their initial values.     

Generalized restricted circular three-body problem as a model for dynamics of binary asteroids

Exploration of the Solar System has recently revealed the existence of a large number of asteroids with satellites, which has stimulated interest in studying the dynamics of such systems. This paper is dedicated to the analysis of the relative motion of a binary asteroid. The conditions of existence of such a system (i.e., when its components do not run away) are derived in the Introduction. Then it is assumed that the satellite has no significant effect on the motion of the main asteroid, the latter being modeled as a dumbbell-like precessing solid body. The equations of motion of this system are a two-parameter generalization of the corresponding equations of the restricted circular three-body problem. It is demonstrated that in the system under consideration there exist steady-state motions in which the small asteroid is equidistant from attracting centers at the ends of the dumbbell (an analog to triangle libration points). The conditions of existence of such motions are derived, and the positions with respect to the dumbbell are analyzed in detail. Examination of the stability of the triangle libration points is reduced to investigation of a characteristic equation of the sixth degree. The stability conditions are derived in the case when the main asteroid executes near-planar motion.     

Families of periodic solutions to Hamiltonian systems: Nonsymmetrical periodic solutions for a planar restricted three-body problem

A monodromy matrix calculated at a single arbitrary point of the periodic solution to a Hamiltonian system allows one to obtain both the direction of continuation for the family of solutions of the first (in Poincarés sense) kind and the multiplicity and direction of branching for periodic solutions of the second kind. In case of resonances 1 : 1 and 1 : 2 one needs to take into account the structure of elementary divisors of the monodromy matrix. Using the planar circular restricted three-body problem as an example, the infiniteness of the process of branching for a nonintegrable system and its finiteness for an integrable system are demonstrated. It is proved that periodic solutions of both first and second kinds which are obtained by continuation of symmetric periodic solutions of a restricted problem are also symmetric. The only exception is the case of resonance 1 : 1 and two second-order cells of the monodromy matrix in the Jordanian form. In this case, all periodic solutions of the second kind turned out to be nonsymmetrical. Examples of the families of nonsymmetrical periodic solutions are given.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 2, 2005, pp. 88–110.Original Russian Text Copyright © 2005 by Kreisman.     

One Class of Motions for a Satellite Carrying a Strong Magnet

The attitude motion of an artificial satellite carrying a strong magnet is studied. The approximate first integrals of the problem, i.e., adiabatic invariants, are indicated. The basic properties of the satellite motions close to the regular precessions with slowly varying parameters are established via the analysis of the adiabatic invariants.     

Application of the Method of Integral Manifolds for Construction of Resonant Curves for the Problem of Spacecraft Entry into the Atmosphere

The separation of motions into slow (precession) and fast (nutation) components in the problem of the entry of a spacecraft (SC) with a small asymmetry into the atmosphere is considered. For the separation of the slow and fast motions the method of integral manifolds is used together with the asymptotic method for singularly perturbed systems. The separation of motions allows one to isolate the frequencies that are functions of the slow variables of a system, and further on, after determining the integer relations between them, to construct the resonant curves (surfaces). This method gives the possibility to analyze the conditions of the emergence of resonances for a SC at angles of attack that are not small and when aerodynamic characteristics are nonlinear. Examples of the construction of resonant curves for a SC with typical aerodynamic characteristics are considered.     

On the one-body problem in the field of a tall cylinder in orbit

We describe the motion of a material point that is attracted by a tall, homogeneous cylinder. The cylinder moves on a circular orbit and points in the radial direction. The finite solution can be expanded in powers of two small parameters: the initial eccentricity of the orbital motion about the cylinder and the coupling parameter of the Coriolis force. Such expansions are given for the case of the indifferent longitudinal location. Solutions periodic in the polar angle are constructed. We do not find resonance effects or other large cumulative amplitude variations that might be due to the gravity-gradient and Coriolis force.     

Dynamics of an axisymmetric satellite under the action of gravitational and aerodynamic torques

Dynamics of attitude motion of an axisymmetric satellite moving in a circular orbit under the action of gravitational and aerodynamic torques is investigated. All equilibrium positions of the satellite in the orbital coordinate system are determined numerically, and sufficient conditions of stability of the equilibrium positions are derived.     

Maximum Lyapunov Exponents for Chaotic Rotation of Natural Planetary Satellites

Based on the Chirikov approach [1, 2] within the context of the theory of separatrix mappings, we suggest and substantiate a simple method for estimating the maximum Lyapunov characteristic exponent (MLCE) of motion in a chaotic layer in the neighborhood of nonlinear resonance separatrix of a Hamiltonian system under an asymmetric periodic perturbation. For a number of natural planetary satellites, using this method, the estimates are made of the MLCE of chaotic rotation (relative to the center of mass) in the main chaotic layer, i.e., in the chaotic layer in the neighborhood of a synchronous resonance separatrix. The value of the MLCE is determined by an orbital eccentricity and by the parameter of the satellite's dynamic asymmetry. The quantity inverse to the MLCE gives a typical time of predictable rotational motion.     

Planar oscillations of a vibrating dumbbell-like body in a central field of forces

The problem of planar motions of a dumbbell-like body of variable length in a central field of Newtonian attraction is considered both in the exact formulation and in satellite approximation. In the satellite approximation the true anomaly of the center of mass is used as an independent variable, which has allowed us to represent the equation of planar oscillations of the dumbbell in the form similar to the Beletskii equation. Some exact solutions to the inverse problem are given both in the complete formulation and in satellite approximation. Under an assumption of small orbit eccentricity and amplitude of the dumbbell vibrations the conditions of existence are found for families of almost periodic motions and splitting separatrices. The phenomena of alternation of regular and chaotic motions are established numerically, as well as chaos suppression with increasing frequency of vibrations. Using the method of averaging the stabilization of tangent equilibria is proved to be impossible.     

Satellite dynamics under the influence of gravitational and aerodynamic torques. A study of stability of equilibrium positions

The dynamics of the rotational motion of a satellite, moving in the central Newtonian force field under the influence of gravitational and aerodynamic torques, is investigated. The paper proposes a method for determining all equilibrium positions (equilibrium orientations) of a satellite in the orbital coordinate system for specified values of aerodynamic torque and the major central moments of inertia; the sufficient conditions for their existence are obtained. For each equilibrium orientation the sufficient stability conditions are obtained using the generalized energy integral as the Lyapunov function. The detailed numerical analysis of the regions where the stability conditions of the equilibrium positions are satisfied is carried out depending on four dimensionless parameters of the problem. It is shown that, in the general case, the number of satellite’s equilibrium positions, for which the sufficient stability conditions are satisfied, varies from 4 to 2 with an increase in the value of the aerodynamic torque magnitude.