航天器的起旋和消旋运动

本文利用轴对称刚体在自身力矩作用下,绕定点运动的一阶正规型运动微分方程代替欧拉方程,讨论航天器的起旋和消旋运动。给出两个卡尔丹角为小量时运动方程的解析积分,由此直接导出航天器相对质心的动量矩矢量及自旋轴在起旋和消旋运动中的变化规律。并以伽里略航天器作为具体对象进行了数值计算。     

The Use of Quaternions in the Optimal Control Problems of Motion of the Center of Mass of a Spacecraft in a Newtonian Gravitational Field: I

The problem of optimal control is considered for the motion of the center of mass of a spacecraft in a central Newtonian gravitational field. For solving the problem, two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems. Both the variants have a quaternion variable among the phase variables. In the first variant this variable characterizes the orientation of an instantaneous orbit of the spacecraft and (simultaneously) the spacecraft location in this orbit, while in the second variant only the instantaneous orbit orientation is specified by it. The suggested equations are convenient in the respect that they allow the general three-dimensional problem of optimal control by the motion of the spacecraft center of mass to be considered as a composition of two interrelated problems. In the first variant these problems are (1) the problem of control of the shape and size of the spacecraft orbit and (2) the problem of control of the orientation of a spacecraft orbit and the spacecraft location in this orbit. The second variant treats (1) the problem of control of the shape and size of the spacecraft orbit and the orbit location of the spacecraft and (2) the problem of control of the orientation of the spacecraft orbit. The use of quaternion variables makes this consideration most efficient. The problem of optimal control is solved on the basis of the maximum principle. Several first integrals of the systems of equations of the boundary value problems of the maximum principle are found. Transformations are suggested that reduce the dimensions of the systems of differential equations of boundary value problems (without complicating them). Geometrical interpretations are given to the transformations and first integrals. The relation of the vectorial first integral of one of the derived systems of equations (which is an analog of the well-known vectorial first integral of the studied problem of optimal control) with the found quaternion first integral is considered. In this paper, which is the first part of the work, we consider the models of motion of the spacecraft center of mass that employ quaternion variables. The problem of optimal control by the motion of the spacecraft center of mass is investigated on the basis of the first variant of equations of motion. An example of a numerical solution of the problem is given.     

The Use of Quaternions in the Optimal Control Problems of Motion of the Center of Mass of a Spacecraft in a Newtonian Gravitational Field: III

The results of numerical solution of the problem of a rendezvous in the central Newtonian gravitational field of a controlled spacecraft with an uncontrollable spacecraft moving along an elliptic Keplerian orbit are presented. Two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems and using quaternion variables to describe the orientations of these coordinate systems. The problem of a rendezvous of two spacecraft is formulated [1, 2] as a problem of optimal control by the motion of the center of mass of a controlled spacecraft with a movable right end of the trajectory, and it is solved on the basis of Pontryagin's maximum principle. The paper is a continuation of papers [1, 2], where the problem of a rendezvous of two spacecraft has been considered theoretically using the two above variants of the equations of motion for the center of mass of the controlled spacecraft.     

微重力环境中重力梯度加速度引起的流体晃动运动

本文研究了关于旋转轴在贮箱的非对称轴上且远离贮箱的几何中心情况下，流体在微重力环境中由重力梯度加速度诱发的晃动特性，建立了问题的数学模型并对模型进行了数值模拟。以高级Ｘ射线天文物理实验卫星（简称ＡＸＡＦ─Ｓ）作为研究对象，获得了由旋转运动引起的重力梯度加速度的数学表达式。晃动问题的数值计算以与卫星固连的非惯性坐标系为基础，目的是寻求一种较为易处理且适合于流体力学方程的边界和初始条件。通过数值计算获得了流体作用于卫星贮箱上的力和力矩。     

Exact analytic solution for the spin-up maneuver of an axially symmetric spacecraft

The problem of spinning-up an axially symmetric spacecraft subjected to an external torque constant in magnitude and parallel to the symmetry axis is considered. The existing exact analytic solution for an axially symmetric body is applied for the first time to this problem. The proposed solution is valid for any initial conditions of attitude and angular velocity and for any length of time and rotation amplitude. Furthermore, the proposed solution can be numerically evaluated up to any desired level of accuracy. Numerical experiments and comparison with an existing approximated solution and with the integration of the equations of motion are reported in the paper. Finally, a new approximated solution obtained from the exact one is introduced in this paper.     

The use of a movable telescoping end mass system for the time-optimal control of spinning spacecraft

The time-optimal control of a spin-stabilized spacecraft with a movable telescoping appendage (boom), is considered analytically and numerically. The motion of a control mass at the end of the boom is determined such that the terminal time will be minimized for two-axis control of a symmetric spacecraft. The equations of rotational motion are linearized about the desired state of spin about the symmetry axis. The equations for the transverse angular velocity components have the form of a coupled two dimensional harmonic oscillator with boom motion as a control force. The control function which brings the system to the desired state is known to be a series of positive and negative pulses. If the initial state is such that the system can be driven to rest in a single switch, the responses, switching and final times, and required boom motion may be determined analytically. Some typical numerical results based on these solutions are discussed.     

Quaternion regularization and trajectory motion control in celestial mechanics and astrodynamics: II

Problems of regularization in celestial mechanics and astrodynamics are considered, and basic regular quaternion models for celestial mechanics and astrodynamics are presented. It is shown that the effectiveness of analytical studies and numerical solutions to boundary value problems of controlling the trajectory motion of spacecraft can be improved by using quaternion models of astrodynamics. In this second part of the paper, specific singularity-type features (division by zero) are considered. They result from using classical equations in angular variables (particularly in Euler variables) in celestial mechanics and astrodynamics and can be eliminated by using Euler (Rodrigues-Hamilton) parameters and Hamilton quaternions. Basic regular (in the above sense) quaternion models of celestial mechanics and astrodynamics are considered; these include equations of trajectory motion written in nonholonomic, orbital, and ideal moving trihedrals whose rotational motions are described by Euler parameters and quaternions of turn; and quaternion equations of instantaneous orbit orientation of a celestial body (spacecraft). New quaternion regular equations are derived for the perturbed three-dimensional two-body problem (spacecraft trajectory motion). These equations are constructed using ideal rectangular Hansen coordinates and quaternion variables, and they have additional advantages over those known for regular Kustaanheimo-Stiefel equations.     

Optimal Control of a Programmed Turn of a Spacecraft

A problem of optimal turn of a spacecraft is considered. The time of turn is minimized, as well as the functional having a meaning of the propellant consumption. An analytical solution to the problem stated is derived. It is demonstrated that the solution optimal in this sense belongs to a class of two-impulse controls, under which a spacecraft executes the turn along the trajectory of its free motion. The solution obtained in this paper differs from earlier available solutions considerably. The estimations of the propellant consumption for a realization of the programmed turn are made.     

The Motion of Coaxial Bodies of Varying Composition on the Active Leg of Descent

The motion of a spacecraft (SC) with double rotation and variable mass on the active leg of its descent is considered. The SC consists of two coaxial bodies. The coaxial scheme is used for gyroscopic stabilization of the SC longitudinal axis by the method of partial spin-up. The equations of spatial motion of coaxial bodies of varying composition are derived and approximate solutions for the angles of spatial orientation are found. The condition of decreasing amplitude of nutation oscillations is obtained, which allows the estimation of efficiency of the stabilization by partial spin-up. The errors in the magnitude and direction of the vector of braking thrust are also determined.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 3, 2005, pp. 224–232.Original Russian Text Copyright © 2005 by Aslanov, Doroshin, Kruglov.     

Study of the mode of angular velocity damping for a spacecraft at non-standard situation

Non-standard situation on a spacecraft (Earth’s satellite) is considered, when there are no measurements of the spacecraft’s angular velocity component relative to one of its body axes. Angular velocity measurements are used in controlling spacecraft’s attitude motion by means of flywheels. The arising problem is to study the operation of standard control algorithms in the absence of some necessary measurements. In this work this problem is solved for the algorithm ensuring the damping of spacecraft’s angular velocity. Such a damping is shown to be possible not for all initial conditions of motion. In the general case one of two possible final modes is realized, each described by stable steady-state solutions of the equations of motion. In one of them, the spacecraft’s angular velocity component relative to the axis, for which the measurements are absent, is nonzero. The estimates of the regions of attraction are obtained for these steady-state solutions by numerical calculations. A simple technique is suggested that allows one to eliminate the initial conditions of the angular velocity damping mode from the attraction region of an undesirable solution. Several realizations of this mode that have taken place are reconstructed. This reconstruction was carried out using approximations of telemetry values of the angular velocity components and the total angular momentum of flywheels, obtained at the non-standard situation, by solutions of the equations of spacecraft’s rotational motion.     

Optimal control of the deployment process of solar wings on spacecraft

The optimal attitude control problem of spacecraft during the stretching process of solar wings is investigated in this paper. The dynamical equations of the nonholonomic system are derived from the conservation principle of the angular momentum of the multibody system. Attitude control of the spacecraft with internal motion is reduced to a nonholonomic motion planning problem. The spacecraft attitude control is transformed into the steering problem for a drift free control system. The optimal solution for steering a spacecraft with solar wings is presented. The controlled motion of spacecraft is simulated for two cases. The numerical results demonstrate the effectiveness of the optimal control approach.     

追踪星跟踪空间非合作目标的相对轨道设计

对空间非合作目标跟踪飞行可以执行观测或监视等任务。首先从一般性出发对追踪星与非合作目标之间的椭圆轨道相对运动方程进行分析,给出具有任意初始条件的相对运动方程解析表达式。其次,对追踪星沿航向跟踪目标并考虑约束条件时的相对轨道设计进行分析后,给出设计追踪星轨道的方法,该方法使得追踪星在保持对地定向的同时也满足测量敏感器的约束条件。最后通过数学仿真进行了验证。     

The Use of Quaternions in the Optimal Control Problems of Motion of the Center of Mass of a Spacecraft in a Newtonian Gravitational Field: II

The problem of a rendezvous in the central Newtonian gravitational field is considered for a controlled spacecraft and an uncontrollable spacecraft moving along an elliptic Keplerian orbit. For solving the problem, two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems and using quaternion variables to describe the orientations of these coordinate systems. In the first variant of the equations of motion a quaternion variable characterizes the orientation of an instantaneous orbit of the spacecraft and the spacecraft location in the orbit, while in the second variant it characterizes the orientation of the plane of the spacecraft instantaneous orbit and the location of a generalized pericenter in the orbit. The quaternion variable used in the second variant of the equations of motion is a quaternion osculating element of the spacecraft orbit. The problem of a rendezvous of two spacecraft is formulated as a problem of optimal control by the motion of the center of mass of a controlled spacecraft with a movable right end of the trajectory, and it is solved on the basis of Pontryagin's maximum principle.     

Stability of large damped flexible spacecraft with stored angular momentum

Vibrational stability of large flexible structurally damped spacecraft carrying internal angular momentum and undergoing large rigid body rotations is analysed modeling the systems as elastic continua. Initially, analytical solutions to the motion of rigid gyrostats under torque-free conditions are developed. The solutions to the gyrostats modeled as axisymmetric and triaxial spacecraft carrying three and two constant speed momentum wheels, respectively, with spin axes aligned with body principal axes are shown to be complicated. These represent extensions of solutions for simpler cases existing in the literature. Using these solutions and modal analysis, the vibrational equations are reduced to linear ordinary differential equations. Equations with periodically varying coefficients are analysed applying Floquet theory. Study of a few typical beam- and plate-like spacecraft configurations indicate that the introduction of a single reaction wheel into an axisymmetric satellite does not alter the stability criterion. However, introduction of constant speed rotors deteriorates vibrational stability. Effects of structural damping and vehicle inertia ratio are also studied.     

Influence of disturbances on the angular motion of a spacecraft in the powered section of its descent

The motion of a variable-mass spacecraft is considered in the powered section of a descending trajectory. Approximate analytical solutions are obtained for the angles of spatial orientation of the spacecraft, which allows one to analyze the nutation motion and to develop recommendations on the spacecraft’s mass configuration, providing the smallest possible deviations of the longitudinal axis and thrust vector from specified directions. The errors of stabilization of the spacecraft’s longitudinal axis are calculated by means of numerical integration of complete models and using the obtained analytical solutions, the results being in good agreement.     

Features of rotational motion of a spacecraft descending in the Martian atmosphere

Angular motion at atmospheric entry is studied in the paper for a spacecraft with a bi-harmonic moment characteristic. Special attention is given to the case when the spacecraft possesses two stable balanced positions, and, hence, it can oscillate in dense atmospheric layers in the ranges of small or large angles of attack. The averaged equations of spacecraft motion are derived, which allow one to increase the speed of calculations by several orders of magnitude. A real example is presented, which concerns a spacecraft specially designed for descending in the Martian atmosphere.     

Estimation of the Probability of Capture into a Resonance Mode of Motion for a Spacecraft during Its Descent in the Atmosphere

The estimation of the probability of capture into a resonance mode of motion is considered for a spacecraft with a small asymmetry during its entry into the atmosphere. It is assumed that the initial conditions of spacecraft motion are distributed uniformly in some sufficiently small domain. The problem is solved for the equations of spacecraft motion linear with respect to the angle of attack. An analytical estimate of the probability of the spacecraft capture into the resonance corresponding to an ascending branch of the velocity head is obtained. The emphasis in the analysis of the estimate is made on the effect of the spacecraft asymmetry type on the probability of capture. A comparison of the estimate with the results of numerical computation is carried out. A model problem concerning the construction of the domain of the spacecraft center of mass locations, most dangerous from the point of view of the realization of the stable resonant modes of motion, is solved.     

Determination of parameters of attitude motion of a small communication satellite using the data of measurements of current of solar panels

A communication satellite (small spacecraft) injected into a geosynchronous orbit is considered. Flywheel engines are used to control the rotational spacecraft motion. The spacecraft after the emergency situation has passed into a state of uncontrolled rotation. In this case, no direct telemetric information about parameters of its rotational motion was accessible. As a result, the problem arose to determine the rotational satellite motion according to the available indirect information: current taken from the solar panels. Telemetric measurements of solar panel current obtained on the time interval of a few hours were simultaneously processed by the least squares method integrating the equations of rotational satellite motion. We present the results of processing 10 intervals of the measurement data allowing one to determine the real rotational spacecraft motion and to estimate the total angular momentum of flywheel engines.     

Ballistics, navigation, and motion control for a spacecraft during its landing on the surface of Phobos

Basic concepts and algorithms laid as foundations of the scheme of landing on the Martian moon Phobos (developed for the Phobos-Grunt project) are presented. The conditions ensuring the landing are discussed. Algorithms of onboard navigation and control are described. The equations of spacecraft motion with respect to Phobos are considered, as well as their use for correction of the spacecraft motion. The algorithm of estimation of the spacecraft’s state vector using measurements with a laser altimeter and Doppler meter of velocity and distance is presented. A system for modeling the landing with a firmware complex including a prototype of the onboard computer is described.     

Optimization of passive orbit with the use of gravity maneuver

A high-precision method of calculating gravitational interactions is applied in order to determine optimal trajectories. A number of problems, necessary for determination of optimal parameters at a launch of a spacecraft and during its flyby near celestial bodies, are considered. The spacecraft trajectory was determined by numerical integration of the equations of passive motion of the spacecraft and of the equations of motion for planets, the Sun, and the Moon. The optimal trajectory of the spacecraft approaching the Sun is determined by fitting its initial conditions.