Spin-stabilized satellite magnetic attitude control scheme without initial detumbling

The angular motion of an axisymmetrical satellite equipped with an active magnetic attitude control system is considered. The dynamics of the satellite are analytically studied on the whole control loop. The control loop is as follows: preliminary reorientation along with nutation damping, spinning about the axis of symmetry, then precise reorientation of the axis of symmetry in inertial space. Reorientation starts right after separation from the launch vehicle. Active magnetic attitude control system time-response with respect to its parameters is analyzed. It is proven that low-inclined orbit forces low control system time-response. Comparison with the common control scheme shows the time-response gain. Numerical analysis of the disturbances effect is carried out and good pointing accuracy is proved.     

Study of a bunch of three algorithms for magnetic control of attitude and spin rate of a spin-stabilized satellite

We consider the angular motion of an axi-symmetrical satellite equipped with an active magnetic attitude control system. Dynamics of the satellite is studied on the entire control loop, consisting of a bunch of three successively used algorithms. The control cycle includes the stages of nutation damping, spinning up the satellite about its symmetry axis, and reorienting the symmetry axis into a preset direction in the inertial space. The results are confirmed by numerical simulations.     

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.     

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.     

Design of satellite constellations with continuous coverage on elliptic orbits of <Emphasis Type="Italic">Molniya</Emphasis> type

New methods of choosing the structures of satellite constellations (SC) on elliptical orbits of the Molniya type are presented. The methods, using critical inclination and putting the orbit apogee in the Earth’s hemisphere with an area of continuous coverage, are based on geometrical analysis of two-dimensional representation of the coverage conditions and SC motion in the space of inertial longitude of the orbit ascending node and time. The coverage conditions are represented in the form of a certain region. Dynamics of all satellites in this space is represented by uniform motion along a straight line approximately parallel to the ordinate axis, while the satellite system forms a grid. The problem of choosing a minimal (as far as the number of satellites is concerned) SC configuration can be formulated as a search for the most sparse grid. The contemporary advanced methods of computational geometry serve as an algorithmic basis for the problem solution. Design of SC for continuous coverage of latitude belts with the use of kinematically regular systems is considered. A method of analyzing single-track systems for continuous coverage of arbitrary geographic regions is described, which makes a region at any time instant observable by at least one satellite of the system. As an example, SC on elliptical orbits are considered with periods of ～4, 12, and 24 hours.     

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.     

Stabilization of a Reentry Vehicle by a Partial Spin-up during Uncontrolled Descent

Stabilization of a reentry vehicle (RV) by a partial spin-up of it is considered for the case of uncontrolled descent into the atmosphere. In this case, the vehicle is a composite construction consisting of two rigid bodies, a return capsule and a stabilizing block, which is put in rotation. A model is developed for the spatial motion of the reentry vehicle considered as a system of coaxial rigid bodies rotating about a common axis of symmetry. The free motion is studied, and the stability of steady-state regimes is analyzed. The spatial motion of the system is considered for the case of a small asymmetry due to displacement of the axes of dynamic symmetry of the bodies with respect to the spin axis, and approximate solutions for the motion parameters of the free system are found.     

基于机器学习的卫星姿态控制律设计

传统的优化算法只能针对给定控制律的指定参数进行优化.机器学习的方法,不仅进化参数,而且可以按给定准则进化出卫星姿态控制律表达式.建立卫星姿态动力学模型及敏感器和执行机构的原理与误差模型,构成含噪声的单自由度的闭环数字仿真系统.统计一段时间内仿真结果的姿态精度和稳定度作为该控制律的适应度函数.针对姿态控制律选定合适的函数...     

Decay of high eccentricity satellite orbits with air drag using uniformly regular KS canonical elements

A non-singular analytical theory for the motion of high eccentricity satellite orbits under the influence of air drag is developed in terms of the Uniformly Regular Kustaanheimo and Stiefel (URKS) canonical elements, by assuming the atmosphere to be oblate with variation of density scale height with altitude. The series expansions include up to fourth-order terms of an independent variable Δ=λ² (function of eccentric anomaly) and c (a small parameter dependent on the flattening of the atmosphere). Only two of the nine equations are solved analytically due to symmetry in the equations of motion. Comparison of the important orbital parameters semi-major axis and eccentricity up to 1000 revolutions, obtained with the present analytical solution and the KS theory, shows the superiority of the present solution over the KS elements analytical solution. The theory can be used effectively for the orbital decay of aero-assisted orbital transfer orbits during mission planning.     

Effects of uncertainties and flexible dynamic contributions on the control of a spacecraft full-coupled model

One of the most important problems for performing a good design of the spacecraft attitude control law is connected to its robustness when some uncertainty parameters are present on the inertial and/or on the elastic characteristics of a satellite. These uncertainties are generally intrinsic on the modeling of complex structures and in the case of large flexible structures they can be also attributed to secondary effects associated to the elasticity. One of the most interesting issues in modeling large flexible space structures is associated to the evaluation of the inertia tensor which in general depends not only on the geometric ‘fixed’ characteristic of the satellite but also on its elastic displacements which of course in turn modify the ‘shape’ of the satellite. Usually these terms can be considered of a second order of magnitude if compared with the ones associated to the rigid part of a structure. However the increasing demand on the dimension of satellites due to the presence for instance of very large solar arrays (necessary to generate power) and/or large antennas has the necessity to investigate their effects on their global dynamic behavior in more details as a consequence. In the present paper a methodology based on classical Lagrangian approach coupled with a standard Finite Element tool has been used to derive the full dynamic equations of an orbiting flexible satellite under the actions of gravity, gravity gradient forces and attitude control. A particular attention has been paid to the study of the effects of flexibility on the inertial terms of the spacecraft which, as well known, influence its attitude dynamic behavior. Furthermore the effects of the attitude control authority and its robustness to the uncertainties on inertial and elastic parameters has been investigated and discussed.     

Control of satellite imaging formations in multi-body regimes

Libration point orbits may be ideal locations for satellite imaging formations. Therefore, control of these arrays in multi-body regimes is critical. A continuous feedback control algorithm is developed that maintains a formation of satellites in motion that is bounded relative to a halo orbit. This algorithm is derived based on the dynamic characteristics of the phase space near periodic orbits in the circular restricted three-body problem (CR3BP). By adjusting parameters of the control algorithm appropriately, satellites in the formation follow trajectories that are particularly advantageous to imaging arrays. Image reconstruction and coverage of the (u, v) plane are simulated for interferometric satellite configurations, demonstrating potential applications of the algorithm and the resulting motion.     

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.     

一种空间推扫型光学成像系统几何映射方法

卫星光学遥感器在轨成像会受到颤振的影响,因此图像品质受姿态稳定度的影响很大,而高分辨率图像受到平台误差的影响更加明显。通过研究空间相机的几何模型,提出了一种从地球坐标系到空间光学遥感器坐标系之间的转换关系;分析了在轨卫星的姿态误差和运动源,并在几何模型中加入了内外方位元素特征;然后进行了空间TDICCD相机的成像仿真实验。为了在像面上模拟颤振,分别进行不同模态的颤振仿真,并且对视线范围内多模态综合作用进行仿真。仿真实验是基于三个轴系方向的,并且计算了TDICCD仿真图像的几何畸变量,模型中几何畸变量的测量尺度达到亚像素级,结果有利于指导卫星平台和遥感器的参数设计。     

Determination of attitude motion of the Foton M-2 satellite using the data of onboard measurements of its angular velocity

We present the resutls of a prompt determination of the uncontrolled attitude motion of the Foton M-2 satellite, which was in orbit from May 31 to June 16, 2005. The data of onboard measurements of the angular velocity vector were used for this determination. The measurement sessions were carried out once a day, each lasting 83 min. Upon terminating a session, the data were transmitted to the ground to be processed using the least squares method and integrating the equations of motion of the satellite with respect to its center of mass. As a result of processing, the initial conditions of motion during a session were estimated, as well as parameters of the mathematical model used. The satellite’s actual motion is determined for 12 such sessions. The results obtained in flight completely described the satellite’s motion. This motion, having begun with a small angular velocity, gradually became faster, and in two days became close to the regular Euler precession of an axisymmetric solid body. On June 14, 2005 the angular velocity of the satellite with respect to its longitudinal axis was approximately 1.3 degrees per second, and the angular velocity projection onto a plane perpendicular to this axis had a magnitude of about 0.11 degrees per second. The results obtained are consistent with more precise results obtained later by processing the data on the Earth’s magnetic field measured on the same satellite, and they complement the latter in determination of the motion in the concluding segment of the flight, when no magnetic measurements were performed.     

On the planar inertial motion of a system of three coupled bodies

Numerous papers are devoted to studying the motion of a system (coupling) of two bodies in the Earth’s satellite orbit ([1–4] and others). The problem on the planar inertial motion of three bodies, coupled by a non-extensible weightless string having the form of an unfastened chain, is considered in the paper. Such a configuration can be represented, for example, by a system of two coupled spacecraft rotating around their common center of mass (in order to simulate the gravity force) in long-term space missions, when the third body (the lift) is located on a connecting cable. The bodies are considered to be the material points (particles).     

船载惯性测量单元姿态信息优化方法

随着海洋卫星通信技术的广泛应用，船载通信系统的需求量与日俱增。基于微电子系统的惯性测量单元主要用于测量船体姿态信息，实现姿态反馈与控制补偿，是船载卫星通信设备的核心器件。实际应用中通过卡尔曼滤波处理的惯性测量单元的姿态信息并不完善，经过分析惯性测量单元在运动过程中存在干扰性的有害加速度。分析了有害加速度的产生原因和对测量姿态的影响并提出姿态信息的优化方法。通过此优化方法可提升惯性测量单元的测量精度，增强船载通信系统的控制精度。     

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.     

基于T-H方程的卫星轨迹模型参考输出跟踪控制方法

当目标卫星沿椭圆轨道运行时,描述追踪星与目标星相对运动的线性化方程为T-H方程。将描述航天器相对运动的T-H方程变换为周期系统的状态空间形式,并给出卫星轨迹跟踪控制问题的数学描述。基于周期系统的参量Lyapunov方法和模型参考跟踪控制理论,提出了卫星轨迹跟踪控制器的设计方法。利用该方法设计了带有收敛速率保障的反馈镇定控制器和具有自由参数的前馈控制器。对追踪星相对目标星悬停任务进行了数值仿真,仿真结果表明提出的控制方案是有效的。     

Model development and adaptive imbalance vibration control of magnetic suspended system

A system model is developed to describe the translational and rotational motion of an active-magnetic-bearing-suspended rigid rotor of a single-gimbal control moment gyro (SCMG) onboard a rigid satellite. This model closely reflects the motion characteristics of the rotor by considering the dynamic and static imbalance as well as the coupling between the gimbal's and the rotor's motion on a satellite platform. Adaptive autocentering control is strictly constructed for the preceding rotor with unknown dynamic imbalance. The rotor achieves its rotation about the principal axis of inertia by identifying the little rotational angles from the geometric axis to the principal axis and then using the results to tune a stabilizing controller, which is composed of a decentralized PD controller with cross-axis proportional gains and high-pass and low-pass filters. The main disturbance in the wheel spinning can thereby be completely removed and the vibration acting on the satellite can be attenuated.     

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.