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.     

General dynamics of a large class of flexible satellite systems

The paper presents a general formulation for librational dynamics of satellites with an arbitrary number, type and orientation of deploying flexible appendages. In particular, the case of beam-type flexible appendages deploying from a satellite in an arbitrary orbit is considered. The governing nonlinear, nonautonomous and coupled equations for vibration of the appendages and libration of the satellite are integrated numerically. Several cases of practical importance are considered making the system progressively more general and hence complex: (i) planar case representing pitch and appendage oscillations in the orbital plane; (ii) general attitude motion with planar vibrations of flexible members; (iii) above two cases together with the out-of-plane component of vibrations. Results show that under critical combinations of the system parameters the combined effect of flexibility and deployment can be substantial.     

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.     

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.     

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.     

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.     

Stability of triangle libration points in generalized restricted circular three-body problem

Equilibrium positions of a small-mass body are considered with respect to a precessing dumbbell. The dumbbell represents two rigidly fixed spherical gravitating bodies. Such a system can be considered as a model of a binary asteroid. Stability of relative equilibrium positions with equal distances from the small mass to the attracting centers is studied. By analogy with the classical restricted three-body problem, these positions are referred to as triangle libration points. It is shown that the character of stability of these libration points is determined by three constant parameters: nutation angle and angular velocity of precession, as well as the ratio of masses at the ends of the dumbbell. Stability conditions are derived in the linear approximation, and the regions of stability and instability in the space of problem parameters are constructed. The paper is a continuation of [1].     

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.     

Hydrodynamic Analogy for the Solution to the Doubly Averaged Hill Problem Disturbed by Precession of the Orbit Plane of a Disturbing Body

A hydrodynamic analogy for the solution to a satellite version of the restricted three-body problem is considered with allowance made for precession of the disturbing body orbit. Applying, in the first approximation, the analytical solution to a doubly averaged equation system obtained in 1961 by M.L. Lidov, we consider the precession of the disturbing body orbit as some disturbing factor. The appropriate model enables us to describe the general pattern of integral curves in the entire area of their existence as a picture of flow lines of the potential motion of some fictitious homogeneous incompressible fluid. For the synthesis of an adequate mathematical model, a method similar to the Pade approximation is used. The obtained model and its discrete continual analog allow us to present, in a pictorial form, the disturbance of orbits under the effect of precession of the disturbing body orbit.     

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].     

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

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

Transient attitude dynamics of satellites with deploying flexible appendages

The paper presents a general formulation for librational dynamics of satellites with an arbitrary number, types, and orientation of deploying flexible appendages. The generalized force term is incorporated making the formulation applicable to a wide variety of situations where aerodynamic forces, solar radiation, earth's magnetic field, etc. become significant. In particular, the case of a beam-type flexible appendage deploying from a satellite in an arbitrary orbit is considered. The corresponding nonlinear, non-autonomous equations for in-plane and out-of-plane vibrations are derived, allowing for the variation of mass density and flexural rigidity along the length with time dependent deployment velocity and spin rate. Next, the attention is focused on the linearized analysis of the in-plane vibrational equation using the assumed-mode method and its substantiation through numerical integration. Finally, the paper presents results for both steady-state and transient attitude behaviour for a representative gravity gradient configuration for a range of initial conditions and system parameters. Results show the combined effect of flexibility and deployment on the dynamics of the system to be substantial. Disturbance of the appendage can excite large amplitude librations. On the other hand, the converse situation is not necessarily true. Furthermore, Coriolis loading, induced by the extending appendages, can become a limiting factor in arriving at a deployment strategy; an effect not pointed out in the literature.     

卫星太阳能帆板的撞击问题

本文研究卫星帆板展开过程中的撞击问题，用子结构方法和单向递推组集方法行到经铰坐标和柔性体模态坐标为系统广义坐标的撞击阶段柔性多体系统的动力学方程，由与撞击有关约束方程的Ｌａｇｒａｎｇｅ乘子导出撞击力的计算公式，计算了卫星帆板在撞击过程中撞击力的最大值和卫星姿态角速度的变化规律。计算结果表明，对柔性体之间的碰撞，卫星姿态角速度在变化过程中呈现上下波动的趋势。     

电磁卫星编队位置跟踪滑模变结构控制

电磁卫星编队位置跟踪是存在输入延时和扰动不确定的控制问题,以卫星相对运动CW方程作为动力学模型,采用前馈加反馈的方式控制编队运动跟踪预定轨迹,前馈控制量由预定轨迹结合CW方程给出,反馈控制采用状态线性变换后的滑模变结构控制律.证明了滑模切换面可达以及滑模面上的运动与不确定扰动无关.数值仿真表明考虑时滞的控制相比未考虑时滞的控制,位置跟踪精度有所提高.     

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.     

簇状卫星姿态控制方程的通用建模方法

当前卫星通常都由一个中心刚体和其他挠性附件构成，每个挠性附件都与中心刚体直接相连，组成了一个簇状的刚柔耦合的多体系统。随着挠性附件的数目的不同，卫星系统的拓扑构型也将有所不同。即使对于同一颗卫星而言，在展开过程的不同阶段，卫星系统具有不同的拓扑构型。本文利用柔性多体系统动力学的单向递推组集方法，提出了一套通用的解决方案，以推导具有不同拓扑构型的簇状卫星的姿态控制方程，得出了姿态控制方程的系数矩阵。     

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.     

Determination of aerodynamic forces on satellites by theory and wind tunnel experiments

Theoretical and experimental aspects of satellite aerodynamics in free molecular flow have been studied. A FORTRAN-program package has been developed to calculate free molecular aerodynamic force and moment coefficients of general and special convex shapes, using an approximation of the surfaces by plane triangular elements.As a special case for concave bodies, the drag of a satellite of TD1A-similar shape has been investigated by an approximate theory and by wind tunnel experiments. Comparisons are made between experimental results, the approximate theory and exact theory.     

Nutation damping of spinning satellite with viscous annular dampers

A new assumption, that the fluid in an annular damper is supposed to be a variable slug, is taken to discuss the nutation damping of a spinning satellite with a partially filled viscous annular damper. The dynamical equations for the spinning satellite with the damper are established by application of the equations of the first order in normal form for the rigid body. Approximate analytical solutions of nutation angle attenuation are obtained by the successive approximation and the small parameter methods. The theoretical results from the variable slug can explain the experimental phenomena, which are not explainable by the results from the conventional rigid assumption, taken by most the former references. Thus the presented analytical solutions provide a more reliable basis for designing the annular damper.