Dynamics of the orbital cable system

The Orbital Cable System (OCS) consists of a vehicle with an engine, a long connecting cable (string) and a satellite probe for low altitude geophysical explorations. The vehicle orbits the Earth at an altitude of 200–250 km and tows a probe at an altitude of 100–150 km. The air drag is compensated by the vehicle engine thrust. In the OCS dynamics study the cable weight, elasticity and air drag are taken into account. The OCS stationary motions are found, and stability of these motions is investigated. The OCS periodical motions that are close to the stationary motions are found. Simple ways of the OCS deployment in orbit are described.     

On the stability of stationary motions of a system of coaxial bodies

We consider the stability of stationary motions of a model of a spacecraft as a system of coaxial bodies with small asymmetry caused by the shift of the axes of dynamic symmetry of bodies relative to the axis of rotation. We determine the stationary motions of the system; their stability is studied with respect to both the projections of angular velocity and the position of the axis of rotation. The sufficient conditions for the stability of these stationary motions are obtained by constructing a Lyapunov function, and the necessary conditions are obtained by analyzing the corresponding linearized equations of perturbed motion.     

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.     

星载合成孔径雷达的运动效应

对星载合成孔径雷达的运动效应进行了详细的分析，给出了雷达回波的多普勒参数和距离徒动特性的一般关系式。详细论证了卫星运动对多普勒参数和成像的影响，导出了描述速度误差、加速度误差和姿态误差的一系列定量计算公式，并就运动误差的补偿方法进行了阐述。     

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.     

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.     

空间绳系机器人抓捕非合作目标的质量特性参数辨识

针对空间绳系机器人抓捕非合作目标/空间垃圾后需要对其进行回收/拖曳的精确控制问题,提出了一种利用抓捕后保持阶段的振动特性辨识目标参数的方法。首先,根据质量特性参数辨识的需要,推导了系统的动力学模型。不同于以往将本体卫星和被抓捕目标简化为质点的动力学模型,本文针对任意的目标抓捕位置,在考虑重力梯度影响的基础上,利用拉格朗日法获得系统各广义坐标的动力学公式。然后,分析非合作目标和系绳在后抓捕保持阶段的姿态运动。最后,在非合作目标与本体卫星没有任何信息交互的情况下,利用后抓捕阶段目标卫星和系绳特有的振动,并使用具有鲁棒性可遗忘因子的递推最小二乘法,提出了包括转动惯量和质心到任意抓捕点距离在内的质量特性参数辨识算法。     

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.     

IGSO和MEO卫星用于TWSTFT的卫星运动误差研究

针对倾斜地球同步轨道（IGSO）和中高地球轨道（MEO）卫星用于卫星双向时间传递（TWSTFT）时的卫星运动所引入的误差进行了研究,并对IGSO和MEO卫星应用于TWSTFT的可行性进行了分析。首先推导了TWSTFT中的卫星运动误差公式,并利用两行根数（TLE）对目前在轨的IGSO和MEO卫星的运动进行了计算,在此基础...     

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.     

Analysis of a double gimbaled reaction wheel spacecraft attitude stabilization system

Three axis attitude stabilization of a satellite using a single spinning reaction wheel mounted on a two degree-of-freedom passively and actively torqued gimbal system is investigated. The passive control is assumed to be provided by a spring-loaded damper mounted on each of the gimbal axes, while active control results from both the wheel acceleration and the torque applied about the gimbal axes. The stability of the uncontrolled and passively controlled systems is investigated analytically. For constant wheel speed the pitch motion is decoupled from the roll-yaw and gimbal motions. Control laws for the roll-yaw motion are developed based on pole clustering and linear optimal control theory. For the pitch motion control laws are obtained based on classical second order system theory. Estimation techniques are applied to the roll-yaw system for the case when the complete state may not be directly observable (in the absence of a fine yaw position sensor).     

Periodic motions of a nearly dynamically symmetric satellite in the neighborhood of hyperboloidal precession

The motion of a satellite close to a dynamically symmetric solid body in a Newtonian gravitational field over a circular orbit is studied. The system of differential equations describing the body’s motion is close to a system with cyclic coordinate. New classes of periodic motions are constructed in the neighborhood of a known partial solution to an unperturbed problem, hyperboloidal precession of a dynamically symmetric satellite. In the resonance case, when the ratio of one frequency of small oscillations of a reduced system with two degrees of freedom in the neighborhood of a stable equilibrium position to the frequency of cyclic coordinate variation is close to an integer number, there exist one or three families of periodic motions that are analytic in terms of fractional powers of a small parameter. A study of stability of these motions was performed with the help of KAM (Kolmogorov-Arnold-Moser) theoty. Faling the described resonance there exists a unique family of periodic motions that is analytic in terms of integer powers of a small parameter. The check-up of stability of these motrons was carried out. We distinguished the cases of parametric resonance, resonances of the third and fourth orders, and a non-resonant case. In the resonance cases our study relies on well-known results on stability of Hamiltonian systems during resonances [1]. In the non-resonant case we use the KAM theory [2].     

卫星双向时间频率传递中的Sagnac效应

 结合卫星双向时间频率传递技术,提出了在卫星运动情况下Sagnac效应解决方案。卫星双向时间频率传递技术是基于地球同步卫星进行的。地球同步卫星在各种摄动力的影响下,相对于地面上的某点不是绝对静止的,而是作小幅度日周期性运动。 Sagnac 效应与卫星和地面观测站的位置密切相关,卫星的运动直接导致了 Sagnac 效应也具有日周期变化的特征。使用高精度实测卫星轨道数据对 Sagnac 效应进行计算分析,结果表明：由卫星运动引起的Sagnac效应值具有与卫星运动一致的周日变化规律特征,大小达到几百皮秒的量级,与传统的将同步卫星作为静止点处理相比,提高了Sagnac效应误差的修正精度。该方案对于各种高精度卫星时间比对技术和卫星导航等领域具有重要的应用价值。     

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.     

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.     

The Use of Television Signals of Stationary Television Centers for Navigation Determination of Low-Orbit Satellites

For systems of autonomous satellite navigation, it is proposed to use television signals of stationary television centers positioned along the flight path in the band of satellite radio coverage as sources of primary navigation information. Television signal processing allows one to determine the time instants of transit through traverses of the television centers. An algorithm of calculating the parameters of satellite motion by times of transit through traverses is developed.     

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.     

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.     

绳系卫星系统中的周期运动

在椭圆轨道上飞行的绳系卫星系统的重要特征是系绳的方向角具有周期性运动，称为驻环。不稳定的驻环导致全系统的不稳定性。本文给出系统运动的数学模型、距离速率控制方法、求解驻环的数值迭代算法、驻环的稳定性和吸引域以及控制参数平面上驻环存在的区界和计算机仿真结果     

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.