Evolution of Motion of a Binary Planet

A model of a binary planet, consisting of a material point of small mass and a deformable viscoelastic sphere, is suggested. The center of mass of the binary planet moves in the gravitational field of a central body in the plane, which contains planets forming the binary planet. A deformable spherical planet rotates around the axis orthogonal to the plane of planetary motion. Planet deformations are described by the linear theory of viscoelasticity. It is shown that with an appropriate approximation of the gravitational potential, there is a class of quasicircular orbits, when the eccentricities of an orbit of the center of mass of a binary planet and an orbit, describing mutual planet motion, are equal to zero. The further evolution of motion is investigated in this class of orbits with the use of the canonical Poincare–Andoyer variables. Corresponding averaged equations are found, and phase pictures are constructed; the stability of stationary solutions is investigated on the basis of equations in variations. For the Solar system planets with their satellites, forming binary planets, the application of the presented model allows us to conclude that satellites sooner or later will fall on the corresponding planets.     

Secular evolution of rotary motion of a charged satellite in a decaying orbit

An electrostatically charged Earth satellite whose orbit is decaying due to the Earths oblateness is considered. Secular perturbations of the orbit are taken into account: they are caused by the second zonal harmonic of the geopotential. These perturbations represent deviations of the longitude of the ascending node and perigee argument, the orbit form being invariable and the orbit inclination to the equatorial plane being constant. The attitude rotary motion of the satellite under the action of perturbing moments of the gravitational and Lorentz forces is studied. The magnetic field of the Earth is taken in a quadrupole approximation. The evolution of the satellites rotary motion is investigated on the basis of new differential equations in s-parameters specially constructed for this purpose. Using the method of averaging, basic regularities of the secular evolution of rotary motion of a screened satellite are revealed. It is found that the rotary motion of a charged satellite essentially depends on the quadrupole component of the geomagnetic potential.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 2, 2005, pp. 111–125.Original Russian Text Copyright © 2005 by Tikhonov.     

Translational–Rotational Motion of a Viscoelastic Sphere in Gravitational Field of an Attracting Center and a Satellite

We study the translational–rotational motion of a planet modeled by a viscoelastic sphere in the gravitational fields of an immovable attracting center and a satellite modeled as material points. The satellite and the planet move with respect to their common center of mass that, in turn, moves with respect to the attracting center. The exact system of equations of motion of the considered mechanical system is deduced from the D'Alembert–Lagrange variational principle. The method of separation of motions is applied to the obtained system of equations and an approximate system of ordinary differential equations is deduced which describes the translational–rotational motion of the planet and its satellite, taking into account the perturbations caused by elasticity and dissipation. An analysis of the deformed state of the viscoelastic planet under the action of gravitational forces and forces of inertia is carried out. It is demonstrated that in the steady-state motion, when energy dissipation vanishes, the planet's center of mass and the satellite move along circular orbits with respect to the attracting center, being located on a single line with it. The viscoelastic planet in its steady-state motion is immovable in the orbital frame of reference. It is demonstrated that this steady-state motion is unstable.     

Evolution of Rotational Motion of a Symmetrical Satellite with Viscoelastic Rods

We study the motion of a symmetrical satellite with a pair of flexible viscoelastic rods in a central Newtonian gravitational field. A restricted problem formulation is considered, when the satellite's center of mass moves along a fixed circular orbit. A small parameter is introduced which is inversely proportional to the stiffness of flexible elements. Another small parameter is equal to the ratio of the squared orbital angular velocity and the squared magnitude of the initial angular velocity of the satellite. In order to describe the satellite rotational motion relative to the center of mass, we use the canonical Andoyer variables. In the undisturbed formulation of the problem, i.e., at = 0 and = 0, these variables are the action–angle variables. Equations describing the evolution of motion are derived by an asymptotic method which combines the method of separating motions for systems with an infinite number of degrees of freedom and the Krylov–Bogolyubov method for systems with fast and slow variables. The manifolds of stationary motions are found, and their stability is investigated on the basis of equations in variations. Phase portraits are constructed which describe the rotational motion of a satellite at the stage of slow dissipative evolution.     

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.     

Reconstruction of spacecraft rotational motion using a Kalman filter

Quasi-static microaccelerations of four satellites of the Foton series (nos. 11, 12, M-2, M-3) were monitored as follows. First, according to measurements of onboard sensors obtained in a certain time interval, spacecraft rotational motion was reconstructed in this interval. Then, along the found motion, microacceleration at a given onboard point was calculated according to the known formula as a function of time. The motion was reconstructed by the least squares method using the solutions to the equations of satellite rotational motion. The time intervals in which these equations make reconstruction possible were from one to five orbital revolutions. This length is increased with the modulus of the satellite angular velocity. To get an idea on microaccelerations and satellite motion during an entire flight, the motion was reconstructed in several tens of such intervals. This paper proposes a method for motion reconstruction suitable for an interval of arbitrary length. The method is based on the Kalman filter. We preliminary describe a new version of the method for reconstructing uncontrolled satellite rotational motion from magnetic measurements using the least squares method, which is essentially used to construct the Kalman filter. The results of comparison of both methods are presented using the data obtained on a flight of the Foton M-3.     

The Dynamics of a Satellite-Gyrostat with a Single Nonzero Component of the Vector of Gyrostatic Moment

The dynamics of a satellite-gyrostat moving in the central Newtonian force field along a circular orbit is studied. In the particular case when the vector of gyrostatic moment is parallel to one of the satellite’s principal central axes of inertia, all the equilibrium states are determined. For each equilibrium orientation, sufficient conditions of stability are obtained as a result of the analysis of the generalized energy integral, and necessary conditions of stability are determined as a result of analysis of the linearized equations of motion. The evolution of regions of validity for the conditions of stability of equilibrium positions are studied in detail depending on the parameters of the problem. All the bifurcation values of the parameters at which qualitative changes of the regions of stability take place are determined.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 4, 2005, pp. 283–294.Original Russian Text Copyright © 2005 by Sarychev, Mirer, Degtyarev.     

Estimating the accuracy of the technique of reconstructing the rotational motion of a satellite based on the measurements of its angular velocity and the magnetic field of the Earth

The paper has studied the accuracy of the technique that allows the rotational motion of the Earth artificial satellites (AES) to be reconstructed based on the data of onboard measurements of angular velocity vectors and the strength of the Earth magnetic field (EMF). The technique is based on kinematic equations of the rotational motion of a rigid body. Both types of measurement data collected over some time interval have been processed jointly. The angular velocity measurements have been approximated using convenient formulas, which are substituted into the kinematic differential equations for the quaternion that specifies the transition from the body-fixed coordinate system of a satellite to the inertial coordinate system. Thus obtained equations represent a kinematic model of the rotational motion of a satellite. The solution of these equations, which approximate real motion, has been found by the least-square method from the condition of best fitting between the data of measurements of the EMF strength vector and its calculated values. The accuracy of the technique has been estimated by processing the data obtained from the board of the service module of the International Space Station (ISS). The reconstruction of station motion using the aforementioned technique has been compared with the telemetry data on the actual motion of the station. The technique has allowed us to reconstruct the station motion in the orbital orientation mode with a maximum error less than 0.6° and the turns with a maximal error of less than 1.2°.     

Evolution of asteroid orbits at the 3 : 1 their mean motion resonance with Jupiter (planar problem)

The 3 : 1 mean motion resonance is studied in the planar elliptic restricted three-body problem (Sun-Jupiter-asteroid). Using double numerical averaging, the equations are constructed that describe the secular evolution of eccentricity and perihelion longitude of the asteroid orbit. The region of adiabatic chaos is isolated in the phase space of the system under study.     

Determination of the Attitude Motion of a Satellite Based on a Kinematical Model of Motion and Using Data of Measurements of Its Angular Velocity and the Earth’s Magnetic Field Strength

An integral statistical procedure of determination of the attitude motion of a satellite using the data of onboard measurements of angular velocity vectors and the strength of the Earth’s magnetic field (EMF) is suggested. The procedure uses only the equations of kinematics of a solid body and is applicable to determining both controlled and uncontrollable motions of a satellite at any external mechanical moments acting upon it. When applying this procedure, the data of measurements of both types, accumulated during a certain interval of time, are processed jointly. The data of measuring the angular velocity are smoothed by discrete Fourier series, and these series are substituted into kinematical Poisson equations for elements of the matrix of transition from a satellite-fixed coordinate system to the orbital coordinate system. The equations thus obtained represent a kinematical model of the satellite motion. The solution to these equations (which approximate the actual motion of a satellite) is found from the condition of the best (in the sense of the least squares method) fit of the data of measuring the EMF strength vector to its calculated values. The results of testing the suggested procedure using the data of measurements of the angular velocity vectors onboard the Foton-12 satellite and measurements of EMF strengths are presented.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 4, 2005, pp. 295–305.Original Russian Text Copyright © 2005 by Abrashkin, Volkov, Voronov, Egorov, Kazakova, Pankratov, Sazonov, Semkin.     

The Use of Quaternions in the Optimal Control Problems of Motion of the Center of Mass of a Spacecraft in a Newtonian Gravitational Field: I

The problem of optimal control is considered for the motion of the center of mass of a spacecraft in a central Newtonian gravitational field. For solving the problem, two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems. Both the variants have a quaternion variable among the phase variables. In the first variant this variable characterizes the orientation of an instantaneous orbit of the spacecraft and (simultaneously) the spacecraft location in this orbit, while in the second variant only the instantaneous orbit orientation is specified by it. The suggested equations are convenient in the respect that they allow the general three-dimensional problem of optimal control by the motion of the spacecraft center of mass to be considered as a composition of two interrelated problems. In the first variant these problems are (1) the problem of control of the shape and size of the spacecraft orbit and (2) the problem of control of the orientation of a spacecraft orbit and the spacecraft location in this orbit. The second variant treats (1) the problem of control of the shape and size of the spacecraft orbit and the orbit location of the spacecraft and (2) the problem of control of the orientation of the spacecraft orbit. The use of quaternion variables makes this consideration most efficient. The problem of optimal control is solved on the basis of the maximum principle. Several first integrals of the systems of equations of the boundary value problems of the maximum principle are found. Transformations are suggested that reduce the dimensions of the systems of differential equations of boundary value problems (without complicating them). Geometrical interpretations are given to the transformations and first integrals. The relation of the vectorial first integral of one of the derived systems of equations (which is an analog of the well-known vectorial first integral of the studied problem of optimal control) with the found quaternion first integral is considered. In this paper, which is the first part of the work, we consider the models of motion of the spacecraft center of mass that employ quaternion variables. The problem of optimal control by the motion of the spacecraft center of mass is investigated on the basis of the first variant of equations of motion. An example of a numerical solution of the problem is given.     

The Use of Quaternions in the Optimal Control Problems of Motion of the Center of Mass of a Spacecraft in a Newtonian Gravitational Field: II

The problem of a rendezvous in the central Newtonian gravitational field is considered for a controlled spacecraft and an uncontrollable spacecraft moving along an elliptic Keplerian orbit. For solving the problem, two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems and using quaternion variables to describe the orientations of these coordinate systems. In the first variant of the equations of motion a quaternion variable characterizes the orientation of an instantaneous orbit of the spacecraft and the spacecraft location in the orbit, while in the second variant it characterizes the orientation of the plane of the spacecraft instantaneous orbit and the location of a generalized pericenter in the orbit. The quaternion variable used in the second variant of the equations of motion is a quaternion osculating element of the spacecraft orbit. The problem of a rendezvous of two spacecraft is formulated as a problem of optimal control by the motion of the center of mass of a controlled spacecraft with a movable right end of the trajectory, and it is solved on the basis of Pontryagin's maximum principle.     

Nonlinear Evolution of a Balloon Satellite Orbit

The motion of a spherically symmetric balloon satellite near the equatorial plane is considered. Taking the Earth's oblateness and solar light pressure into account, the integral of motion can be obtained under certain simplifications. The eccentricity is related to the solar angle which represents an angle between pericenter and the Sun. This analytical approximation describes a large and complicated evolution of the eccentricity in corresponding areas of the phase space and the space of parameters. Phase portraits contain fixed saddle points and separatrices that divide different types of oscillations of the eccentricity. In the unsimplified problem, separatrices break down, and specific stochastic motions arise. The aims of the present study are (1) evaluation of the accuracy of analytical approximation with the help of numerical integration using a sufficiently complete model of motion and (2) numerical investigation of stochastic motions and dimensions of stochastic zones in the region of broken separatrices for an adequate model of motion. For a balloon satellite with a semimajor axis of 2.15 Earth's radii and a windage of 30 cm²/g the dimensions of a stochastic zone in eccentricity and solar angle are 10^–5and 0.1°, respectively. The analytical approximation describes the orbit evolution in the right way, except for the cases of large eccentricities, e> 0.4, which corresponds to a pericenter height of less than 1400 km, where the atmospheric drag is already significant.     

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.     

The use of a movable telescoping end mass system for the time-optimal control of spinning spacecraft

The time-optimal control of a spin-stabilized spacecraft with a movable telescoping appendage (boom), is considered analytically and numerically. The motion of a control mass at the end of the boom is determined such that the terminal time will be minimized for two-axis control of a symmetric spacecraft. The equations of rotational motion are linearized about the desired state of spin about the symmetry axis. The equations for the transverse angular velocity components have the form of a coupled two dimensional harmonic oscillator with boom motion as a control force. The control function which brings the system to the desired state is known to be a series of positive and negative pulses. If the initial state is such that the system can be driven to rest in a single switch, the responses, switching and final times, and required boom motion may be determined analytically. Some typical numerical results based on these solutions are discussed.     

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.     

Study of the mode of angular velocity damping for a spacecraft at non-standard situation

Non-standard situation on a spacecraft (Earth’s satellite) is considered, when there are no measurements of the spacecraft’s angular velocity component relative to one of its body axes. Angular velocity measurements are used in controlling spacecraft’s attitude motion by means of flywheels. The arising problem is to study the operation of standard control algorithms in the absence of some necessary measurements. In this work this problem is solved for the algorithm ensuring the damping of spacecraft’s angular velocity. Such a damping is shown to be possible not for all initial conditions of motion. In the general case one of two possible final modes is realized, each described by stable steady-state solutions of the equations of motion. In one of them, the spacecraft’s angular velocity component relative to the axis, for which the measurements are absent, is nonzero. The estimates of the regions of attraction are obtained for these steady-state solutions by numerical calculations. A simple technique is suggested that allows one to eliminate the initial conditions of the angular velocity damping mode from the attraction region of an undesirable solution. Several realizations of this mode that have taken place are reconstructed. This reconstruction was carried out using approximations of telemetry values of the angular velocity components and the total angular momentum of flywheels, obtained at the non-standard situation, by solutions of the equations of spacecraft’s rotational motion.     

Asymptotic study of a complete magnetic attitude control cycle providing a single-axis orientation

The angular motion of an axisymmetrical satellite equipped with the active magnetic attitude control system is examined. Attitude control system has to ensure necessary orientation of the axis of symmetry in the inertial space. It implements the following strategy: coarse reorientation of the axis of symmetry with nutation damping or “-Bdot” without initial detumbling; spinning-up about the axis of symmetry to achieve the property of a gyro; fine reorientation of the axis in the inertial space. Dynamics of the satellite is analytically studied using averaging technique on the complete control loop consisting of five algorithms. Solutions of the equations of motion are obtained in terms of quadratures for most cases or even in closed-form. The latter allowed to study the dependence of motion parameters including time-response with respect to the orbit inclination and other parameters for all algorithms.     

Rotational motion of <Emphasis Type="Italic">Foton M-4</Emphasis>

The actual controlled rotational motion of the Foton M-4 satellite is reconstructed for the mode of single-axis solar orientation. The reconstruction was carried out using data of onboard measurements of vectors of angular velocity and the strength of the Earth’s magnetic field. The reconstruction method is based on the reconstruction of the kinematic equations of the rotational motion of a solid body. According to the method, measurement data of both types collected at a certain time interval are processed together. Measurements of the angular velocity are interpolated by piecewise-linear functions, which are substituted in kinematic differential equations for a quaternion that defines the transition from the satellite instrument coordinate system to the inertial coordinate system. The obtained equations represent the kinematic model of the satellite rotational motion. A solution of these equations that approximates the actual motion is derived from the condition of the best (in the sense of the least squares method) match between the measurement data of the strength vector of the Earth’s magnetic field and its calculated values. The described method makes it possible to reconstruct the actual rotational satellite motion using one solution of kinematic equations over time intervals longer than 10 h. The found reconstructions have been used to calculate the residual microaccelerations.     

Motion of three viscoelastic planets in the field of their mutual attraction forces

Translational-rotational motion of three planets modeled by viscoelastic balls in the gravitational field of mutual attraction is studied in this paper. The system of equations of motion for the mechanical system under consideration is deduced from the d’Alembert-Lagrange variational principle. Using the method of separation of motions, an approximate system of ordinary differential equations, describing the translational-rotational motion of the planets, is obtained with taking into account perturbations caused by elasticity and dissipation. The found steady-state motion of the system is an analog to triangular libration points in the classical three-body problem.