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.     

On dynamics of a plasma ring rotating in the magnetic field of a central body: Magneto-gyroscopic waves. Problems of stability and quantization

Based on a mathematical model described in [1], some new aspects of the dynamics of a thin planar plasma ring rotating in the magnetic field of a central body are considered. The dipole field is considered assuming that the dipole has a small eccentricity, and the dipole axis is inclined at a small angle to the central body’s axis of rotation. Emphasis is placed on the problem of stability of the ring’s stationary rotation. Unlike [1], the disturbed motion is considered which has a character of eddy magneto-gyroscopic waves. The original mathematical model is reduced to a system of finite-difference equations whose asymptotic analytical solution is obtained. It is demonstrated that some “elite” rings characterized by integral quantum numbers are long-living, while “lethal” or unstable rings (antirings) are associated with half-integer quantum numbers. As a result, an evolutionally rife rotating ring of magnetized plasma turns out to be stratified into a large number of narrow elite rings separated by gaps whose positions correspond to antirings. The regions of possible existence of elite rings in near-central body space are considered. Quantum numbers determining elite eigenvalues of the mean sector velocity (normalized in a certain manner) of a ring coincide with the quantum numbers appearing in the solution to the Schrödinger equation for a hydrogen atom. Perturbations of elite orbits corresponding to these quantum numbers satisfy the de Brogli quantum-mechanical condition. This is one more illustration of the isomorphism of quantization in microcosm and macrocosm.     

二维复合目标电磁散射的仿真研究

为了对E极化波和H极化波入射情形下二维任意形状金属和两种不同电磁参数损耗介质复合目标的电磁散射进行仿真,建立了复合结构目标的电磁积分方程,采用矩量法(MoM)和共轭梯度迭代法(CGM)求解此积分方程。依据所建数学模型编程对一种金属-蜂窝吸波介质复合结构翼面的雷达散射截面(RCS)进行了仿真分析。     

访问多个特定相对位置的航天器轨道设计

在航天器轨道设计问题中,将惯性空间中经典的吉布斯三矢量定轨方法拓展到相对运动空间中,给出了一种相对运动条件下的三矢量定轨方法。针对已知轨道的目标航天器,以及二个或三个给定的空间相对位置,基于相对运动方程,提出了设计跟随航天器飞行轨道的数值方法。以轨道面共面或异面,以及目标航天器轨道形状为椭圆或圆,将问题分为四种情况进行约束条件和自由变量个数的分析讨论。对于自由变量个数多于约束方程的情况,额外给定周期重访约束,将各种情况下的特定相对位置访问问题转化为一至二维的非线性方程(组)求解问题。对一维方程求解采用分段黄金分割+割线法进行快速求解;对二维方程组通过网格法搜索迭代初值并通过牛顿迭代快速求解。进一步基于线性模型的解,采用微分修正方法求解了各情况下J2摄动模型下的结果。数值算例验证了提出方法的正确性及有效性。     

Hori method for generalized canonical systems

In this paper, some special features on the canonical version of Hori method, when it is applied to generalized canonical systems (systems of differential equations described by a Hamiltonian function linear in the momenta), are presented. Two different procedures, based on a new approach for the integration theory recently presented for the canonical version, are proposed for determining the new Hamiltonian and the generating function for systems whose differential equations for the coordinates describe a periodic system with one fast phase. These procedures are equivalent and they are directly related to the canonical transformations defined by the general solution of the integrable kernel of the Hamiltonian. They provide the same near-identity transformation for the coordinates obtained through the non-canonical version of Hori method. It is also shown that these procedures are connected to the classic averaging principle through a canonical transformation. As examples, asymptotic solutions of a non-linear oscillations problem and of the elliptic perturbed problem are discussed.     

考虑状态约束的二体旋转库仑卫星系统重构控制

研究考虑控制输入饱和与状态约束的深空旋转二体库仑卫星构型控制问题,只采用卫星之间的库仑力作为控制力,提出一种基于反步法的非线性控制方法。首先推导了二体库仑卫星在地-月系平动点附近的相对运动方程,利用旋转二体库仑卫星的特性,对方程进行了简化。为了完成禁止相对运动区域的回避,设计了新的状态限制辅助函数,结合抗饱和方法与反步法得到了二体库仑卫星的构型控制器。接着证明了由于状态限制辅助函数的加入,所设计的控制器可以保证卫星相对运动不超出限制范围。进一步应用Lyapunov稳定性定理证明了其闭环系统的一致最终有界性。最后在Matlab/Simulink平台上进行了仿真校验,结果表明了方法的有效性。     

弯曲平板的电磁散射特性分析

根据弯曲导体平板的特点,将电流在水平与竖直向分别进行离散,由电流连续性方程获得电荷离散公式,得到电场积分方程（EFIE）的离散形式。用泰勒展开消除格林函数积分中的奇异值,并通过积分近似求解,给出了一种改进的双共轭梯度迭代法（CGS）快速求解大型线性方程组。仿真结果表明：与原双共轭梯度迭代法相比,该法的收敛速度更快、精度更高。     

空间最优交会

在相对体制下,本文讨论了混合指标下的最优交会问题。除给出了最优交会体制导律的结构形式外,还提供了求解最优交会制导律的参数方程组。该方程具有迭代求解形式,给数字求解带来极大方便。     

A new treatment of Lambertian mechanics

Lambert's formulas for the time of orbital flight between two points in space are rederived by first establishing a universal differential equation governing the time function, independent of the conic type of the trajectory, the focal characteristics of the trajectory sector, and the range angle. A unified form of Lambert's formulas is then obtained as the general solution of the differential equation, and the various forms of the classical Lambert's formulas are obtained as its particular solutions under different boundary conditions. Following this basic treatment, various hypergeometric expansions for Lambert's time function and its derivatives are developed, and the behavior of the function and its implications in the solution of Lambert's problem and the isochronous trajectories are briefly reviewed. Finally, a short comparison of the present treatment with those found in current literature on Lambertian Mechanics is made and briefly discussed.     

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.     

Determination of the Quasistatic Component of Microaccelerations at the MirStation

We describe the determination of the quasistatic component of microaccelerations, as it was done during the space experiments with the instruments DACON and ALICE-2. This component was calculated using the telemetric information related to the motion of the station with respect to its center of mass. The information consists of the values of the station angular velocity vector and the quaternion which specifies its orientation, determined at discrete instants of time. The quaternion is determined with a step of about 1 min, the angular velocity with a step of about 10 s. The information is used in the following manner. At first, the quaternion components corresponding to some time interval are smoothed by splines. Then, employing the obtained splines and kinematic equations, the angular velocity and acceleration of the station are calculated on this interval. Finally, the microacceleration is calculated as a function of time at the point of the location of the instrument. The data of measurements of the angular velocity are used for the purpose of control. As a rule, the available telemetric information allows one to find the quasi-static component for the entire time interval of carrying out the experiment. Examples of determination of this component for some experiments are presented. A comparison with the results of calculating the microacceleration quasistatic component by other methods is made.     

粘弹塑性界面的断裂特性

研究张开型粘弹塑性界面断裂。用傅立叶正、余弦变换及逐段定积分变换方法将边值问题的控制方程化为奇异积分方程组。解方程后计算了裂纹尖端塑性区尺寸及裂纹尖端张开位移(COD：crack—tip opening displacement)，并给出了能量释放率算式。结果表明，裂纹尖端塑性区尺寸和COD均随两种材料的最小屈服极限的增加而减小；随时间的增大，COD先增长后衰减，最后渐近地逼近于定值。     

模型修正中一种特殊广义逆的性质与功能

本文对一种特殊的广义逆矩阵在模型修正中使用时可能具有的性质与功能给予了间接的论证,从而为这种广义逆的进一步应用给予了提示。该论证表明,这种广义逆不仅具有增补方程个数的作用,而且还具有对修正模型的参考基增补高阶理论模态信息的功能。通过本文论证,文[2][4]之解的唯一性问题明确解决了,即为一个实用完备模态空间中的确定解。     

Orbital rendezvous mission planning using mixed integer nonlinear programming

The rendezvous and docking mission is usually divided into several phases, and the mission planning is performed phase by phase. A new planning method using mixed integer nonlinear programming, which investigates single phase parameters and phase connecting parameters simultaneously, is proposed to improve the rendezvous mission’s overall performance. The design variables are composed of integers and continuous-valued numbers. The integer part consists of the parameters for station-keeping and sensor-switching, the number of maneuvers in each rendezvous phase and the number of repeating periods to start the rendezvous mission. The continuous part consists of the orbital transfer time and the station-keeping duration. The objective function is a combination of the propellant consumed, the sun angle which represents the power available, and the terminal precision of each rendezvous phase. The operational requirements for the spacecraft–ground communication, sun illumination and the sensor transition are considered. The simple genetic algorithm, which is a combination of the integer-coded and real-coded genetic algorithm, is chosen to obtain the optimal solution. A practical rendezvous mission planning problem is solved by the proposed method. The results show that the method proposed can solve the integral rendezvous mission planning problem effectively, and the solution obtained can satisfy the operational constraints and has a good overall performance.     

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.     

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.     

正交各向异性层合圆板的非轴对称弯曲

利用层合板的宏观应力应变关系，从三维板壳理论入手，建立了对称铺设正交各向异性层合圆板在外载荷作用下弯曲问题的微分方程。考虑一个复合材料整星隔振适配器与卫星所组成系统的横向运动等效模型，计算出特殊正交各向异性层合圆环板在集中弯矩作用下的解析解，并与有限元计算结果进行对比，证明了解的正确性。     

Theory of type III radio bursts in the interplanetary space: I. Derivation of permittivity tensor for the system consisting of electron beam and solar wind plasma

The tensor of permittivity for the system “electron beam - plasma of the interplanetary space” is derived in the approximation of geometrical optics. The problem is one-dimensional; all parameters such as density of the beam and of the solar wind plasma, and the induction of the interplanetary magnetic field are assumed to be dependent only on the distance to the Sun. The beam is generated by an active region during a solar flare, and it is a source of radio bursts of type III in the interplanetary space. The tensor of permittivity was obtained to close field equations by a material equation. On the basis of these equations it becomes possible to study theoretically the amplitude-frequency characteristics of the radio bursts as disturbances of the above-described beam-plasma system.     

天基空间目标监视的短弧段定轨技术

利用两个短弧段的天基测角资料实现对空间目标的轨道确定是天基空间目标监视系统需解决的重要问题之一。提出利用二体轨道的角动量和能量守恒方程计算空间目标初轨的新方法,该方法将短弧段的主要测量信息表示为弧段属性,由两个短弧段的弧段属性建立守恒方程,采用变量替换法求解守恒方程获得多个初轨,通过方差分析从中选择合适的初轨。针对轨道改进中的迭代发散和收敛于局部极小点的问题,提出了选取多个迭代初值进行轨道改进的采样方法。仿真结果表明,初轨确定算法是可行的,轨道改进能成功解算最小二乘轨道。     

航天员出舱牵引构型设计与结构响应性能评估

针对航天员出舱过程所用安全绳锁挂频繁的问题,进行了安全高效的航天员出舱牵引系统构型设计,可作为太空行走的辅助导向装置。基于材料结构弹塑性力学方程与含应变耦合项的热传导方程,推导结构热力响应可计算数学模型,利用Newmark方法以及Crank-Nicolson数值格式分别对热弹性方程和热传导方程进行时间域上的离散,构造金属结构响应无条件稳定隐式有限元格式,搭建了温度场与位移场相互影响热力耦合有限元算法(FEM)性能评估软件系统,以此对牵引系统的导轨进行热力耦合数值计算,得到在太阳辐射外热流载荷作用下具有较强抗变形能力的结论。计算了航天员采用牵引系统出舱时对导轨产生的拉力,并借助ANSYS Workbench对导轨进行静力学求解,验证了该牵引系统具备结构强度可靠性,可为空间站建设舱外作业平台提供重要理论分析依据与设计参考。