A specific case of stability of cylindrical precession of a satellite

In a central Newtonian gravitational field, the motion of a dynamically symmetrical satellite along an elliptical orbit of arbitrary eccentricity is considered. The particular motion of the satellite is known when its axis of symmetry is perpendicular to the orbit plane, and the satellite rotates about this axis with a constant angular velocity (cylindrical precession). A nonlinear analysis of stability of this motion has been performed under the assumption that the geometry of the satellite mass corresponds to a thin plate. At small values of orbit eccentricity e the analysis is analytical, while numerical analysis is used for arbitrary values of e.     

Orbital and relative motion of a tethered satellite system

A simplified model for the orbital and relative motion of a tethered satellite system is presented. The tether acts as a light elastic string with small structural damping but without bending stiffness. Its mass is taken into account in the calculation of the total kinetic and potential energies of the tethered system. This formulation allows the inclusion of the complete gravity gradient influence on the dynamics of the system. The resulting three-dimensional motion is separated in the centre of mass orbital motion on the one hand and the relative motion of the end-bodies on the other. No restrictions on length of the tether or on mass ratio of the end-masses are imposed. It is found that the frequencies and amplitudes of the longitudinal tether oscillations are realistic as long as the tether mass is less than that of the subsatellite.     

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.     

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.     

Variation of the Satellite Orbit Altitude by the Light Pressure Force

The possibility of the uncontrolled increase of the altitude of an almost circular satellite orbit by the force of the light pressure is investigated. The satellite is equipped with a damper and a system of mirrors (solar batteries can serve as such a system). The flight of the satellite takes place in the mode of a single-axis gravitational orientation, the axis of its minimum principal central moment of inertia makes a small angle with the local vertical and the motion of the satellite around this axis constitutes forced oscillations under the impact of the moment of force of the light pressure. The form of the oscillations and the initial orbit are chosen so that the transverse component of the force of the light pressure acting upon the satellite be positive and the semimajor axis of the orbit would continuously increase. As this takes place, the orbit remains almost circular. We investigate the evolution of the orbit over an extended time interval by the method which employs separate integration of the equations of the orbital and rotational motions of the satellite. The method includes outer and inner cycles. The outer cycle involves the numerical integration of the averaged equations of motion of the satellite center of mass. The inner cycle serves to calculate the right-hand sides of these equations. It amounts to constructing an asymptotically stable periodic motion of the satellite in the mode of a single-axis gravitational orientation for current values of the orbit elements and to averaging the equations of the orbital motion along it. It is demonstrated that the monotone increase of the semimajor axis takes place during the first 15 years of motion. In actuality, the semimajor axis oscillates with a period of about 60 years. The eccentricity and inclination of the orbit remain close to their initial values.     

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.     

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

Averaging of the equation of motion of a satellite with respect to its center of mass in degenerate cases: 1. Oscillations in the absolute frame of reference

The problem of attitude oscillations of a satellite with a small dynamic asymmetry in the plane of the orbit leads to a system degenerate to the fifth order from the point of view of the method of averaging. An explicit expression for the dominant term is obtained by integration in the complex plane. The recurrence procedure of calculating the higher approximations of the method of averaging is considered, as well as an approach to the analysis of the structure of derived expressions.     

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.     

J2摄动作用下近地轨道卫星编队构形长期演化机理分析

研究近地轨道卫星编队构形在引力摄动作用下的长期演化机理。把用转移矩阵法得出的相对运动分析解表达成一般轨道根数形式，为研究摄动影响而进一步表达成无奇点轨道根数形式；根据各轨道根数的J2摄动解分别研究编队构形在轨道坐标系三个方向上的振幅与相位的长期变化情况，揭示了构形引力摄动的长期演化机理；给出编队构形的摄动表达式，并以仿真实例验证其可信性；最后给出构形长期引力摄动的相关结论，并对编队构形设计提出建议，这将为研究编队构形的规划技术带来方便。     

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.     

Dynamics of a gyrostat satellite with the vector of gyrostatic moment in the principal plane of inertia

Dynamics of a gyrostat satellite moving along a circular orbit in a central Newtonian field of force is investigated. In a particular case, when the gyrostatic moment vector lies in one of the satellite’s principal central planes of inertia, all positions of equilibrium are determined, and the conditions of their existence are analyzed. Also determined are bifurcation values of dimensionless parameters, at which the number of equilibrium positions changes. As a result of analysis of the generalized energy integral, for each equilibrium orientation the sufficient conditions of stability are derived. Evolution of the regions where the sufficient conditions of stability are valid is investigated under variation of the system’s dimensionless parameters.     

基于在轨运行状态的卫星稳定性与可靠性研究

为解决卫星在轨运行故障率较高的问题,文章分析了卫星在轨运行的特点,定义了卫星在轨工作状态及其与系统稳定性变化的关系。研究表明:对在轨运行安全稳定状态进行控制,能使状态保持或向高级跃迁,可提升卫星稳态及可靠性。文章最后提出了提高卫星在轨运行稳定性与可靠性的在轨运行策略与技术方向。     

Families of Periodic Solutions to the Beletsky Equation

A review of achievements in the investigation of the planar periodic oscillations and rotations of a satellite around its center of mass, which moves along an elliptic orbit in the central gravitational field, is presented. These oscillations and rotations of the satellite are described by an ordinary differential equation (Beletsky equation) of the second order with periodic coefficients and two parameters. The equation is equivalent to a periodic Hamiltonian system with a single degree of freedom and has a singularity. It turned out that two-parameter families of the generalized periodic solutions to this simple equation form complicated structures of a new type. A comparison of numerous separate results made it possible to outline a sufficiently unified and complete picture of the location and structure of the families of generalized 2-periodic solutions with an integer number of rotation. These families are compared with known data related to the resonance rotations of the natural satellites in the solar system.     

低轨星座对静止轨道卫星同频干扰指标研究

随着全球低轨星座系统的快速部署及应用,非静止轨道卫星对静止轨道卫星系统的同频干扰问题日益凸显。为保护静止轨道卫星系统不受有害干扰,国际电联现行无线电规则中规定了非静止轨道卫星系统应满足的等效功率通量密度限值或干扰噪声比限值,但所适用频段及限值标准仍在不断修订完善中。为了精确定量评估低轨星座对静止轨道卫星的同频干扰,通过对等效功率通量密度限值和干扰噪声比限值的确定方法及演进历程进行研究,分析了上述限值与静止轨道卫星实际被干扰的对应关系,并以典型倾斜轨道星座和极地轨道星座为例,分析了不同干扰指标对卫星系统间干扰判定结果的影响,当静止轨道卫星系统的实际被干扰门限低于国际电联制定等效功率通量密度限值选用的参考门限时,建议低轨星座选择干扰噪声比限值作为开展频率干扰分析与判定的依据。     

Relationship between the fluxes of relativistic electrons at geosynchronous orbit and the level of ULF activity on the Earth’s surface and in the solar wind during the 23rd solar activity cycle

We present the results of a cross-correlation analysis made on the basis of Spearman’s rank correlation method. The quantities to correlate are daily values of the fluence of energetic electrons at a geosynchronous orbit, intensities of ground and interplanetary ultra-low-frequency (ULF) oscillations in the Pc5 range, and parameters of the solar wind. The period under analysis is the 23rd cycle of solar activity, 1996–2006. Daily (from 6 h to 18 h of LT) magnetic data at two diametrically opposite observatories of the Intermagnet network are taken as ground-based measurements. The fluxes of electrons with energies higher than 2 MeV were measured by the geosynchronous GOES satellites. The data of magnetometers and plasma instruments installed on ACE and WIND spacecraft were used for analysis of the solar wind parameters and of the oscillations of the interplanetary magnetic field (IMF). Some results elucidating the role played by interplanetary ULF waves in the processes of generation of magneospheric oscillations and acceleration of energetic electrons are obtained. Among them are (i) high and stable correlation of ground ULF oscillations with waves in the solar wind; (ii) closer link of mean daily amplitudes of both interplanetary and ground oscillations with ‘tomorrow’ values of the solar wind velocity than with current values; and (iii) correlation of the intensity of ULF waves in the solar wind, normalized to the IMF magnitude, with fluxes of relativistic electrons in the magnetosphere.     

Stability of resonance rotation of a satellite with respect to its center of mass in the orbit plane

The stability of resonance oscillations and rotations of a satellite in the plane of its orbit in the case when the difference of the moments of inertia with respect to the principal axes lying in the orbit plane is small is determined at a given rotation number m by the sign of function Φ_m(e), introduced by F.L. Chernous’ko in 1963. In this paper, convenient analytical representations of functions Φ_m(e) are described in the form of integrals and series of Bessel functions regular at e → 1^?. Values of Φ_m(1) are calculated in explicit form. A theorem about the double asymptotic form of functions Φ_m(e) at m → ∞ and e → 1^? is proved by the saddlepoint method.     

近赤道卫星绕飞编队的运动学分析及轨道设计研究

首先对一般运动学方法在进行近赤道卫星编队相对运动分析和轨道设计时存在的问题进行了分析 ;其次介绍了一种基于零倾角轨道变换的运动学新方法 ,用其对近赤道卫星编队中参考卫星轨道倾角对环绕卫星轨道根数的影响进行了研究 ;最后对零倾角卫星编队相对运动方程的线性化误差进行了理论分析 ,并利用数值仿真对两种运动学方法的相对运动方程线性化误差进行了比较分析。数值仿真的结果表明 :由基于零倾角轨道变换运动学方法得到的相对运动方程的线性化误差不随参考卫星轨道倾角改变 ,而由一般运动学得到的相对运动方程的线性化误差随着卫星编队接近赤道而呈非线性增大。     

北斗在轨卫星广播星历精度分析

为分析北斗广播星历的精度,采用激光测距资料和精密星历作为参考进行分析。由于没有提供BDS-3试验卫星精密星历,本文设计了3天弧段的DBS-3试验卫星精密定轨试验,解算的轨道精度约为50 cm,钟差精度约为2 ns。之后联合国际GNSS监测评估系统iGMAS发布的事后BDS-2精密星历,对BDS-3/BDS-2广播星历轨道、钟差精度进行分析。为准确评估广播星历轨道精度,还采用SLR观测数据作为外部检核手段分析北斗在轨卫星广播星历轨道精度。统计了广播星历轨道误差的1D RMS和钟差误差的STD,结果表明：多数钟差精度优于 10 ns ,轨道精度优于5 m。另外为分析广播星历卫星钟的时频性能,计算了在轨卫星钟差的频率稳定度、漂移率和准确度,试验结果表明：BDS-3试验卫星各项性能指标普遍优于BDS-2,频率准确度、频率漂移率和频率稳定性指标计算结果显示BDS-3卫星相对于BDS-2卫星分别提升了42.8%,22.5%,9.5%。