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

Equilibrium positions of a weight on a cable fixed to a dumbbell-shaped space station moving along a circular geocentric orbit

The motion of a space object in the gravitational field of the Earth is considered. The object consists of an extended space station and a weight, which is free to move along the cable fixed to the ends of the station. It is assumed that the station is composed of two masses coupled by a weightless rod, while the cable is weightless and non-stretched. The equations of motion of such a system are derived for the case when the motion proceeds in a single plane, while the center of mass of the system moves along a circular geocentric orbit. The conditions of the cable tension (conditions of being on tie) are derived. The phase portrait of the weight motion along the cable is constructed when the station is oriented to the attracting center or is perpendicular to this position. The possibility to leave the tie in this case is analyzed. Equilibrium configurations of the system are found, i.e., such motions of the object under consideration at which the weight does not change its position relative to the station. Lyapunov stability of such configurations is analyzed for two situations: when the station is composed of equal masses and when masses at the ends of the station are different. In particular, for the case of different masses it is established that there exist such positions of equilibrium at which the dumbbell is located at an angle to the direction to the attracting center. In some cases these positions can be stabilized (if the weight is fixed on the cable).     

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.     

Evolution of the Motion of a Viscoelastic Sphere in a Central Newtonian Field

The evolution of translational-rotational motion of a viscoelastic sphere in a central Newtonian field is studied. By the method of separating motions and averaging in the generalized Andoyer–Delaunay variables, equations are derived that describe the evolution of motion. The analysis of the approximate equations obtained is performed. The condition of existence of a steady-state motion is found, and its stability is investigated using the equations in variations.     

Analysis of the orbital motion of the asteroid Apophis’ satellite

We have analyzed the orbital disturbed spacecraft motion near an asteroid. The equations of the asteroidocentric spacecraft motion have been used with regard to three perturbations from celestial bodies, the asteroid’s nonsphericity, and solar radiation pressure. It has been shown that the orbital parameters of the main spacecraft and a small satellite with a radio beacon can be selected such that the orbits are rather stable for a fairly long period of time, i.e., a few weeks for the main spacecraft with an orbit initial radius of ~0.5 km and a few years before approaching Apophis with the Earth in 2029, for a small satellite at an orbit initial radius of ~1.5 km. The initial orientation of the spacecraft orbital plane perpendicular to the sunward direction is optimal from the point of view of the stability of the spacecraft flight near an asteroid.     

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.     

Equilibrium attitude and stability of a spacecraft on a stationary orbit around an asteroid

The stationary orbits around an asteroid, if exist, can be used for communication and navigation purposes just as around the Earth. The equilibrium attitude and stability of a rigid spacecraft on a stationary orbit around a uniformly-rotating asteroid are studied. The linearized equations of attitude motion are obtained under the small motion assumption. Then, the equilibrium attitude is determined in both cases of a general and a symmetrical spacecraft. Due to the higher-order inertia integrals of the spacecraft, the equilibrium attitude is slightly away from zero Euler angles. Then necessary conditions of stability of this conservative system are analyzed based on the linearized equations of motion. The effects of different parameters, including the harmonic coefficients C₂₀ and C₂₂ of the asteroid and higher-order inertia integrals of the spacecraft, on the stability are assessed and compared. Due to the significantly non-spherical shape and rapid rotation of the asteroid, the effects of the harmonic coefficients C₂₀ and C₂₂ are very significant, while effects of the third- and fourth-order inertia integrals of the spacecraft can be neglected. Considering a spacecraft on a stationary orbit around an example asteroid, we show that the classical stability domain predicted by the Beletskii–DeBra–Delp method on a circular orbit in a central gravity field is modified due to the non-spherical mass distribution of the asteroid. Our results are confirmed by a numerical simulation.     

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.     

Calculation and Study of Limited Orbits around the <Emphasis Type="Italic">L</Emphasis><Subscript>2</Subscript> Libration Point of the Sun–Earth System

A procedure has been proposed for calculating limited orbits around the L₂ libration points of the Sun–Earth system. The motion of a spacecraft in the vicinity of the libration point has been considered a superposition of three components, i.e., decreasing (stable), increasing (unstable), and limited. The proposed procedure makes it possible to correct the state vector of the spacecraft so as to neutralize the unstable component of the motion. Using this procedure, the calculation of orbits around various types of libration points has been carried out and the dependence on the orbit type on the initial conditions has been studied.     

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.     

圆型限制性三体问题相对运动解析研究

针对圆型限制性三体问题共线平动点附近周期/拟周期轨道下的相对运动问题，提出一种新的、通用的解析研究方法。在周期/拟周期轨道近似解析解的基础上，结合微分修正方法，获得了精确的周期/拟周期轨道。对周期/拟周期轨道的单值矩阵进行分析，同时借鉴Floquet理论核心思想，建立了六个相对运动模态，并将相对运动表示为六个相对运动模态的线性组合，获得了相对运动的近似解析解。最后在地-月系统圆型限制性三体问题下，以L1点作为研究对象，分别以Halo轨道、Lissajous轨道和Lyapunov轨道为参考轨道，对相对运动模态和相对运动进行仿真分析，说明了相对运动模态的正确性以及相对运动近似解析解的有效性。     

Generalized problem of two fixed centers

In the problem under consideration, the libration points are determined, and their stability is investigated. The surfaces of zero velocity are constructed, and the regions where motion is possible are found. The limiting case of the problem is considered.     

系统参数对电动力绳系动力学的影响

将电动力绳系(EDT)的主星质量、子星质量、绳系质量以及绳系中的电流视为系统参数，研究这些参数对系统的摆动动力学和轨道动力学的影响。哑铃模型下的电动力绳系摆动动力学方程存在不稳定的周期解，通过Floquet理论来衡量周期解的不稳定程度，从而研究各系统参数对摆动动力学的影响。建立了用春分点轨道元素的形式描述的电动力绳系轨道动力学方程，并以降轨时间来衡量电动力绳系的降轨效率，从而研究系统参数对轨道动力学的影响。运用算例对周期解迁移矩阵的特征值、降轨时间随各系统参数的变化关系进行了仿真，分别得出了各系统参数对系统摆动动力学和轨道动力学的影响。综合本文的仿真结果，并考虑实际发射及空间运行中的其它因素，对电动力绳系的设计和降轨策略提出了建议。     

Solving Optimization Problems for the Flight Trajectories of a Spacecraft with a High-Thrust Jet Engine in Pulse Formulation for an Arbitrary Gravitational Field in a Vacuum

A mathematically well-posed technique is suggested to obtain first-order necessary conditions of local optimality for the problems of optimization to be solved in a pulse formulation for flight trajectories of a spacecraft with a high-thrust jet engine (HTJE) in an arbitrary gravitational field in vacuum. The technique is based on the Lagrange principle of derestriction for conditional extremum problems in a function space. It allows one to formalize an algorithm of change from the problems of optimization to a boundary-value problem for a system of ordinary differential equations in the case of any optimization problem for which the pulse formulation makes sense. In this work, such a change is made for the case of optimizing the flight trajectories of a spacecraft with a HTJE when terminal and intermediate conditions (like equalities, inequalities, and the terminal functional of minimization) are taken in a general form. As an example of the application of the suggested technique, we consider in this work, within the framework of a bounded circular three-point problem in pulse formulation, the problem of constructing the flight trajectories of a spacecraft with a HTJE through one or several libration points (including the case of going through all libration points) of the Earth–Moon system. The spacecraft is launched from a circular orbit of an Earth's artificial satellite and, upon passing through a point (or points) of libration, returns to the initial orbit. The expenditure of mass (characteristic velocity) is minimized at a restricted time of transfer.     

单通道控制旋转弹系统动力学建模与Hopf分岔

为准确预测旋转弹系统的锥形运动形态并判断其稳定性，提出一对鸭舵引起的气动不对称性以及可能引起复杂非线性动力学特性的非线性因素，建立能准确描述单通道控制旋转弹系统姿态运动的复数形式的动力学模型并分析讨论其动态稳定性条件。利用分岔理论对单通道控制旋转弹系统开展Hopf分岔研究，给出了Hopf分岔发生的判断准则；推导了用于判定极限环稳定的第一Lyapunov系数。数值仿真结果验证了条件的正确性与有效性并发现了拟周期运动及混沌运动。研究结果为旋转弹的控制参数设计及结构参数设计提供理论参考。     

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

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

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

同步双小行星系统共振轨道设计

研究了同步双小行星系统中共振轨道的设计方法及演化规律。首先，基于双椭球模型建立探测器运动方程，并给出共振轨道初值选取方法。然后，利用改进并行打靶法，提出一种双小行星系统平面共振轨道两步修正方法。同时结合稳定性理论及分岔理论，给出双小行星系统三维共振轨道生成和延拓方法；最后，以双小行星系统1999KW4为例，设计了共振比为1∶1,1∶2,1∶3,1∶4,2∶3的平面和空间共振轨道族，并分析了共振轨道的特性及轨道周期和轨道能量的变化规律。给出的双小行星系统中共振轨道的设计方法具有普适性，对未来双小行星系统探测任务中的轨道设计具有一定的参考意义与借鉴价值。     

Theory of type III radio bursts in the interplanetary space. II. Study of instability of the electron beam-solar wind plasma system in a linear approximation

Frequency characteristics of disturbances of a one-dimensionally inhomogeneous electron beam-solar wind plasma system are studied in the geometrical optics approximation on the basis of the Maxwell equations closed by the material equation obtained earlier. The corresponding dispersion equation is derived and solved. It is found that resonance interaction of a wave with an electron beam can occur only at two spatial points. Perhaps, such a short-time (point-like) mechanism of the resonance clarifies one of the main problems of physics of electron beams generated by solar flares: their time of existence is much longer than the time following from the previous theoretical estimates of the beam energy loss rate due to radiation.     

On Saturn’s rotation relative to a center of mass under the action of the gravitational moments of the Sun and Jupiter

Saturn’s rotation relative to a center of mass is considered within an elliptic restricted three-body problem. It is assumed that Saturn is a solid under the action of gravity of the Sun and Jupiter. The motions of Saturn and Jupiter are considered elliptic with small eccentricities e_S and e_J, respectively; the mean motion of Jupiter n_J is also small. We obtain the averaged Hamiltonian function for a small parameter of ε = n_J and integrals of evolution equations. The main effects of the influence of Jupiter on Saturn’s rotation are described: (α) the evolution of the constant parameters of regular precession for the angular momentum vector I₂; (β) the occurrence of new libration zones of oscillations I₂ near the plane of the celestial equator parallel to the plane of the Jupiter’s orbit; (γ) the occurrence of additional unstable equilibria of vector I₂ at the points of the north and south poles of the celestial sphere and, as a result, the existence of homoclinic trajectories; and (δ) the existence of periodic trajectories with arbitrarily large periods near the homoclinic trajectory. It is shown that the effects of (β), (γ), and (δ) are caused by the eccentricity e of the Jupiter’s orbit and are practically independent of Jupiter’s mass (within satellite approximation).