Evolution of the Orbit of an Earth Satellite with a Solar Sail

The efficiency of using the light pressure of solar radiation for increasing the semimajor axis of the orbit of an Earth Satellite carrying a solar sail is estimated. The orbit is nearly circular and has an altitude of about 900 km. The satellite is in the mode of single-axis solar orientation: it rotates at an angular velocity of 1 deg/s around the axis of symmetry, which traces the direction to the Sun. This mode is maintained by the solar sail, which serves in this case as a solar stabilizer. The following method of increasing the semimajor axis of the orbit (which is equivalent to increasing the total energy of the satellite's orbital motion) is considered. On those sections of the orbit, where the angle between the light pressure force acting upon the sail and the vector of geocentric velocity of the satellite does not exceed a specified limit, the sail is functioning as a solar stabilizer. On those sections of the orbit, where the above-indicated angle exceeds this limit, the sail is furled by way of turning the edges of the petals towards the Sun. Such a control increases the semimajor axis by more than 150 km for three months of flight. In this case, the accuracy of solar orientation decreases insignificantly.     

Single-Axis Solar Orientation of a Satellite of the Earth

The possibility of using the mode of single-axis solar orientation is considered for a satellite placed into a nearly circular orbit with an altitude of 900 km and bearing a solar sail. The satellite (together with the sail) has an axisymmetric structure, its symmetry axis being the principal central axis of the maximum moment of inertia. The center of the sail pressure lies on this axis and is displaced with respect to the satellite's center of mass. The symmetry axis of the satellite is set to the Sun so that its center of mass would be located between the Sun and the pressure center and would rotate around this axis with an angular velocity of a few degrees per second. The satellite's axis of symmetry makes a slow precession under the action of the gravitational moment and the moment of light pressure forces. Though the maximum magnitudes of these moments are comparable, the moment of the light pressure forces dominates and controls the precession in such a way that the symmetry axis orientation to the Sun remains unchanged.     

Orbital stability of planar oscillations of a satellite in a circular orbit

We consider the attitude motion of a satellite with a circular orbit in a central Newtonian gravitational field. The satellite is a solid body whose mass geometry is that of a plate. A nonlinear analysis is made of orbital stability of planar oscillations of the satellite at which its middle or major axis of inertia is perpendicular to the orbit plane. At small amplitudes of oscillations the analysis of stability was made analytically, while for arbitrary amplitudes the numerical analysis was performed.     

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.     

Realization of modes of satellite attitude motion with a small level of microaccelerations using electromechanical executive devices

Quasi-static microaccelerations are estimated for a satellite specially designed to perform space experiments in the field of microgravity. Three modes of attitude motion of the spacecraft are considered: passive gravitational orientation, orbital orientation, and semi-passive gravitational orientation. In these modes the lengthwise axis of the satellite is directed along the local vertical, while solar arrays lie in the orbit plane. The second and third modes are maintained using electromechanical executive devices: flywheel engines or gyrodynes. Estimations of residual microaccelerations are performed with the help of mathematical modeling of satellite’s attitude motion under the action of gravitational and aerodynamic moments, as well as the moment produced by the gyro system. It is demonstrated that all modes ensure rather low level of quasi-static microaccelerations on the satellite and provide for a fairly narrow region of variation for the vector of residual microacceleration. The semi-passive gravitational orientation ensures also a limited proper angular momentum of the gyro system.     

Bifurcation Regime of a Synchronous Resonance in the Translational–Rotational Motion of Nonspherical Natural Satellites of Planets

Period-doubling bifurcations of the synchronous spin-orbit resonance in the motion of a nonspherical natural planetary satellite along the elliptic orbit are studied. The satellite spin axis is assumed to coincide with the axis of its largest principal moment of inertia and is perpendicular to the orbital plane. The period-doubling bifurcations take place when the value of satellite's dynamical asymmetry parameter falls in the parametric resonance domain. Theoretical dependences of the amplitude of the bifurcation oscillations of a satellite at the pericenter of its orbit upon the eccentricity and dynamical asymmetry parameter are investigated. Three different methods of calculating the amplitude of bifurcation oscillations are presented and compared. These theoretical estimates can be used to predict the opportunity to observe the bifurcation regime. The possibility of the occurrence of the bifurcation regime in the motion of natural planetary satellites is studied. It is concluded that the bifurcation regime is possible in the motion of Deimos, Epimetheus, Helen, Pandora, and Phobos. Phobos is the most probable candidate for finding the bifurcation regime of a synchronous rotation. The identification of such a regime would allow one to impose stringent constraints on the values of the inertial parameters of the satellite observed.     

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.     

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.     

The Mode of Gravitational Orientation for the International Space Station

The results of the determination of the uncontrolled attitude motion of the International Space Station during its unmanned flight in 1999 are presented. The data of onboard measurements of three components of the angular velocity are used for this determination. These data covering an interval of slightly less than one orbit were jointly processed by the least squares method, by integrating the equations of motion of the station relative to its center of mass. As a result of this processing, the initial conditions of the motion and the parameters of the mathematical model used were estimated. The actual motion of the station has been determined for 20 such intervals during April–November. Throughout these intervals, the station rotated about the axis of the minimum moment of inertia, the latter executing small oscillations relative to the local vertical. Such a mode, known as the mode of gravitational orientation of a rotating satellite or the mode of generalized gravitational orientation, was planned before the flight. The measurements were made to verify it. The quasistatic component of the microaccelerations aboard the station is estimated for this mode.     

Compensation of gravitational attraction disturbance to pure gravity orbit for Inner Formation Flying System

The Inner Formation Flying System (IFFS) consisting of an outer satellite and an inner satellite which is a solid sphere proof mass freely flying in the shield cavity can construct a pure gravity orbit to precisely measure the earth gravity field. The gravitational attraction on the inner satellite due to the outer satellite is a significant disturbance source to the pure gravity orbit and is required to be limited to 10⁻¹¹ m s⁻² order. However, the gravitational disturbance force was on 10⁻⁹ m s⁻² order actually and must be reduced by dedicated compensation mass blocks. The region of relative motion of the inner satellite about its nominal position is within 1 cm in dimension, which raises the complexity of the compensation blocks design. The iterative design strategy of the compensation blocks based on reducing the gravitational attraction at the nominal position of the inner satellite is presented, aiming to guarantee the gravitational force in the relative motion region within requirements after the compensation. The compensation blocks are designed according to the current status of IFFS, and the gravitational disturbance force in the region is reduced to 10⁻¹¹ ms⁻² order with minimized adding mass.     

Modeling of the aerodynamic moment acting upon a satellite

We consider the issues of modeling the moments of aerodynamic forces acting upon a satellite with gravitational system of stabilization. It is assumed that satellite orbits are almost circular with heights in the range 550–750 km. Simplified analytical expressions are suggested for the aerodynamic moment in the case when a satellite moves in the regime of gravitational orientation. Accuracy of the obtained expressions is estimated to be compared with that of expressions derived under the assumption of constant coefficient of frontal resistance. An analysis is made of short-periodic variations of the atmosphere density occurring due to orbital motion of a satellite. It is demonstrated that these variations can result in a substantial change of the aerodynamic moment, and their approximation by a truncated Fourier series is suggested. Estimates of the accuracy of the suggested approximation are given.     

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.     

On orbit control using gravity gradient effects

Possibility of orbit control using gravity gradient (GG) effects without any mass expulsion is discussed. For simplicity, a dumb-bell type satellite and circular orbits are mainly considered. It is shown that the GG effects can be applied to convert attitude torques into orbital torques and vice versa. In central gravitational force fields, maximum orbital torques or thrusts are available from the GG force when roll or pitch angle is $\text{±}\text{π}\text{4}$ provided that the attitude angle is null when the dumb-bell axis coincides with the local vertical. Such external torques as geomagnetic or solar wind pressure can be utilized to maintain the $\text{±}\text{π}\text{4}$ attitude, then the orbital torques are available forever. In non-central gravitational fields, without any external torque, the orbital radii of circular orbits can be increased by controlling the satellite attitude using electric energy. The use of the Earth's oblateness effects and the exterior Lunar potential is discussed.     

高精度重力场下回归轨道半解析优化设计

为解决太阳同步回归轨道的标称设计问题,提出一种基于高精度重力场的半解析优化方法。建立地球非球形引力摄动阶数为J15 的高精度重力场解析模型,并分离出引力摄动的长期项和长周期项。构建回归轨道从半长轴到平交点周期的对应关系,平交点周期变化随引力摄动阶数的提高而逐渐收敛。通过微分修正迭代算法所确定的半长轴相对于传统J2摄动模型的半长轴确定值具有更高的精度和更好的稳定性。考察摄动短周期项影响下的密切交点周期,结果表明其受初始位置(平近点角)影响较大,变化范围为0.015s,并由此给出精确回归轨道优化设计的基准：不同的初始位置上满足星下点轨迹严格回归的半长轴期望值。     

Gravity-oriented satellite dynamics subject to gravitational and active damping torques

The dynamics of the rotational motion of a satellite moving in the central Newtonian field of force over a circular orbit under the effect of gravitational and active damping torques, which depend on the satellite angular velocity projections, has been investigated. The paper proposes a method of determining all equilibrium positions (equilibrium orientations) of a satellite in the orbital coordinate system for specified values of damping coefficients and principal central moments of inertia. The conditions of their existence have been obtained. For a zero equilibrium position where the axes of the satellite-centered coordinate system coincide with the axes of the orbital coordinate system, the necessary and sufficient conditions for asymptotic stability are obtained using the Routh–Hurwitz criterion. A detailed analysis of the regions where the conditions of the asymptotic stability of a zero equilibrium position are fulfilled have been obtained depending on three dimensionless parameters of the problem, and the numerical study of the process of attenuation of satellite’s spatial oscillations for various damping coefficients has been carried out. It has been shown that there is a wide range of damping parameters from which, by choosing the necessary values, one can provide the asymptotic stability of satellite’s zero equilibrium position in the orbital coordinate system.     

The use of refined model of aerodynamic moment in the problem of reconstruction of the <Emphasis Type="Italic">Foton</Emphasis> satellite rotational motion

A new mathematical model of the uncontrolled rotational motion of the Foton satellite is presented. The model is based on the Euler dynamic equations of rigid body motion and takes into account the action upon the satellite of four external mechanical moments: gravitational, restoring aerodynamic, moment with constant components in the satellite-fixed coordinate system, and moment arising due to interaction of the Earth’s magnetic field with the satellite’s proper magnetic moment. To calculate the aerodynamic moment a special geometrical model of the outer satellite shell is used. Detailed form of the formulas giving above-mentioned moments in the equations of satellite motion is agreed with the form of the considered motion. Model testing is performed by determining with its help the rotational motion of the Foton M-2 satellite (it was in orbit from May 31, 2005 to June 16, 2005) using the data of the onboard measurements of the Earth’s magnetic field strength. The use of the new model has led to a relatively small improvement in the accuracy of the motion determination, but allowed us to obtain physically real estimates of some parameters.     

Modes of uncontrolled rotational motion of the <Emphasis Type="Italic">Progress M-29M</Emphasis> spacecraft

We have reconstructed the uncontrolled rotational motion of the Progress M-29M transport cargo spacecraft in the single-axis solar orientation mode (the so-called sunward spin) and in the mode of the gravitational orientation of a rotating satellite. The modes were implemented on April 3–7, 2016 as a part of preparation for experiments with the DAKON convection sensor onboard the Progress spacecraft. The reconstruction was performed by integral statistical techniques using the measurements of the spacecraft’s angular velocity and electric current from its solar arrays. The measurement data obtained in a certain time interval have been jointly processed using the least-squares method by integrating the equations of the spacecraft’s motion relative to the center of mass. As a result of processing, the initial conditions of motion and parameters of the mathematical model have been estimated. The motion in the sunward spin mode is the rotation of the spacecraft with an angular velocity of 2.2 deg/s about the normal to the plane of solar arrays; the normal is oriented toward the Sun or forms a small angle with this direction. The duration of the mode is several orbit passes. The reconstruction has been performed over time intervals of up to 1 h. As a result, the actual rotational motion of the spacecraft relative to the Earth–Sun direction was obtained. In the gravitational orientation mode, the spacecraft was rotated about its longitudinal axis with an angular velocity of 0.1–0.2 deg/s; the longitudinal axis executed small oscillated relative to the local vertical. The reconstruction of motion relative to the orbital coordinate system was performed in time intervals of up to 7 h using only the angularvelocity measurements. The measurements of the electric current from solar arrays were used for verification.     

地球磁场对带电赤道卫星的轨道半长轴和偏心率的摄动影响

本文利用摄动理论研究了带电的赤道卫星在地球磁场中做椭圆轨道运动时地磁摄动力对轨道半长径和偏心率的摄动影响。研究结果表明：地磁场对带电赤道卫星的轨道半长径没有摄动影响，既无周期摄动，也无长期摄动，但对轨道偏心率有摄动影响，且只有周期性摄动，而无长期摄动。当卫星自身带电量较大时，这种摄动影响必须予以考虑     

Satellite orbit perturbations in a dusty Martian atmosphere

In this paper we calculate the effect of atmospheric dust on the orbital elements of a satellite. Dust storms that originate in the Martian surface may evolve into global storms in the atmosphere that can last for months can affect low orbiter and lander missions. We model the dust as a velocity-square depended drag force acting on a satellite and we derive an appropriate disturbing function that accounts for the effect of dust on the orbit, using a Lagrangean formulation. A first-order perturbation solution of Lagrange's planetary equations of motion indicates that for a local dust storm cloud that has a possible density of 8.323×10⁻¹⁰ kg m⁻³ at an altitude of 100 km affects the orbital semimajor axis of a 1000 kg satellite up −0.142 m day⁻¹. Regional dust storms of the same density may affect the semimajor axis up to of −0.418 m day⁻¹. Other orbital elements are also affected but to a lesser extent.     

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.