太阳帆日心悬浮轨道附近的相对运动

研究悬浮轨道附近的编队飞行问题，推导了悬浮轨道附近的相 对运动方程. 由于编队太阳帆之间的距离与帆到太阳的距离的比值为小量，将相对运动方程 在悬浮轨道附近线性化，得到了线性化方程. 基于该线性化方程，考虑了悬浮轨道附近的几 种编队控制方法，只需通过调节太阳帆的姿态来进行简单的控制. 其中包括一种被动编队控 制策略，该控制策略具有实现简单、稳定区域大的特点，具有很好的工程应用前景. 最后基 于非线性方程对每种编队策略进行了数值仿真验证，数值结果表明该控制方法能实现编 队.     

Motion of the moonlet in the binary system 243 Ida

The motion of the moonlet Dactyl in the binary system 243 Ida is investigated in this paper. First, periodic orbits in the vicinity of the primary are calculated, including the orbits around the equilibrium points and large-scale orbits. The Floquet multipliers' topological cases of periodic orbits are calculated to study the orbits' stabilities. During the continuation of the retrograde near-circular orbits near the equatorial plane, two period-doubling bifurcations and one Neimark–Sacker bifurcation occur one by one, leading to two stable regions and two unstable regions. Bifurcations occur at the boundaries of these regions. Periodic orbits in the stable regions are all stable, but in the unstable regions are all unstable. Moreover, many quasi-periodic orbits exist near the equatorial plane. Long-term integration indicates that a particle in a quasi-periodic orbit runs in a space like a tire. Quasi-periodic orbits in different regions have different styles of motion indicated by the Poincare sections. There is the possibility that moonlet Dactyl is in a quasi-periodic orbit near the stable region I, which is enlightening for the stability of the binary system.     

Optimal Control of Deployment of a Tethered Subsatellite

One of the most important operations during a tethered satellite system mission is the deployment of a subsatellite from a space ship. We restrict tothe simple but practically important case that the system ismoving on a circular orbit around the Earth. The main problem duringdeployment due to gravity gradient is that the two satellites do not move along the straight radial relative equilibrium position which is stable for a tether of constant length. Instead, deploymentleads to an unstable motion with respect to the radial relativeequilibrium configuration. Therefore we introduce an optimal control strategy using theMaximum Principle to achieve a force controlled deployment of the tethered subsatellite from the radial relative equilibrium position close to the space ship to the radial relative equilibrium position far away from the space ship.     

Adaptive control for autonomous rendezvous of spacecraft on elliptical orbit

A strategy for spacecraft autonomous rendezvous on an elliptical orbit in situation of no orbit information is developed. Lawden equation is used to describe relative motion of two spacecraft. Then an adaptive gain factor is introduced, and an adaptive control law for auton- omous rendezvous on the elliptical orbit is designed using Lyapunov approach. The relative motion is proved to be ultimately bounded under this control law, and the final relative position error can achieve the expected magnitude. Simulation results indicate that the adaptive control law can realize autonomous rendezvous on the elliptical orbit with relative state information only.     

GENERAL FORM RESONANCE ORBITS

为了利用较小的推力使航天器的轨道产生较大的变化,可以利用共振原理来研究航天器的运动,称这样一类非开普勒轨道为共振轨道。将圆频率作为变量,通过合理地选择轨道描述参数、时间尺度和推力描述方式建立了一般形式的共振轨道模型,并基于仿真分析研究了共振轨道圆频率对共振轨道的影响。通过对地球——火星共振转移轨道的算例进行仿真分析,初步研究了共振轨道在星际探测轨道设计中应用的效果。研究结果表明：圆频率改变将对推力峰值产生影响;共振轨道在星际探测中的应用是可行的,并且在能量消耗方面优于Lambert 轨道。     

Long-term evolution of mid-altitude quasi-satellite orbits

Quasi-satellite orbits are of great interest for the exploration of planetary moons because of their dynamical features and close proximity with respect to the surface of scientifically relevant objects like Phobos and Deimos. This paper explores the equations of the elliptical Hill problem, offering a new analytical insight into the long-term evolution of mid-altitude quasi-satellite orbits. Our developments are based on the Yamanaka–Ankersen solution of the Tschauner–Hempel equations and capture the effects of the secondary’s gravity and orbital eccentricity on the shape and orientation of near-equatorial retrograde relative trajectories. The analytical solution of the in-plane and out-of-plane components of the secular motion is achieved by averaging over the relative longitude of a spacecraft as seen from the co-rotating frame of the two primaries. Developments are validated against the numerical integration of quasi-periodic trajectories that densely cover the surface of three-dimensional invariant tori. This analysis confirms the stable nature of quasi-satellite orbits and provides new tools for future spacecraft missions such as the Martian Moons eXploration envisaged by JAXA.

     

A new control law based on characteristic exponent assignment for libration point orbit maintenance

A new method is developed for stabilizing motion on collinear libration point orbits using the formalism of the circular restricted three body problem. Linearization about the collinear libration point orbits yields an unstable linear parameter-varying system with periodic coefficients. Given the variational equations, an innovative control law based on characteristic exponent assignment is introduced for libration point orbit maintenance. A numerical simulation choosing the Richardson＇s third order approximation for a halo orbit as a nominal orbit is conducted, and the results demonstrate the effectiveness of this control law.     

Averaging analyses for spacecraft orbital motions around asteroids

This paper summarizes a few cases of spacecraft orbital motion around asteroid for which averaging method can be applied, i.e., when central body rotates slowly, fast, and when a spacecraft is near to the resonant orbits between the spacecraft mean motion and the central body＇s rotation. Averaging conditions for these cases are given. As a major extension, a few classes of near resonant orbits are analyzed by the averaging method. Then some resulted conclusions of these averaging analyses are applied to understand the stabil- ity regions in a numerical experiment. Some stability conclu- sions are obtained. As a typical example, it is shown in detail that near circular 1 ： 2 resonant orbit is always unstable.     

Displaced orbits generated by solar sails for the hyperbolic and degenerated cases

Displaced non-Keplerian orbits above planetary bodies can be achieved by orientating the solar sail normal to the sun line. The dynamical systems techniques are employed to analyze the nonlinear dynamics of a displaced orbit and different topologies of equilibria are yielded from the basic configurations of Hill’s region, which have a saddlenode bifurcation point at the degenerated case. The solar sail near hyperbolic or degenerated equilibrium is quite unstable. Therefore, a controller preserving Hamiltonian structure is presented to stabilize the solar sail near hyperbolic or degenerated equilibrium, and to generate the stable Lissajous orbits that stay stable inside the stabilizing region of the controller. The main contribution of this paper is that the controller preserving Hamiltonian structure not only changes the instability of the equilibrium, but also makes the modified elliptic equilibrium become unique for the controlled system. The allocation law of the controller on the sail’s attitude and lightness number is obtained, which verifies that the controller is realizable.     

非自旋航天器混沌姿态运动及其参数开闭环控制

研究万有引力场中受大气阻力且存在结构内阻尼的非自旋航天器在椭圆轨道上平面天平动的混沌及其参数开闭环控制问题．在建立数学模型的基础上确定出现混沌的必要条件并数值验证混沌的存在性，提出非线性振动系统混沌运动的参数开闭环控制并应用于控制航天器的混沌姿态运动．     

<Emphasis Type="Italic">J</Emphasis><Subscript>2</Subscript> invariant relative orbits via differential correction algorithm

This paper describes a practical method for finding the invariant orbits in J ₂ relative dynamics. Working with the Hamiltonian model of the relative motion including the J ₂ perturbation, the effective differential correction algorithm for finding periodic orbits in three-body problem is extended to formation flying of Earth’s orbiters. Rather than using orbital elements, the analysis is done directly in physical space, which makes a direct connection with physical requirements. The asymptotic behavior of the invariant orbit is indicated by its stable and unstable manifolds. The period of the relative orbits is proved numerically to be slightly different from the ascending node period of the leader satellite, and a preliminary explanation for this phenomenon is presented. Then the compatibility between J ₂ invariant orbit and desired relative geometry is considered, and the design procedure for the initial values of the compatible configuration is proposed. The influences of measure errors on the invariant orbit are also investigated by the Monte–Carlo simulation. The project supported by the Innovation Foundation of Beihang University for Ph.D. Graduates, and the National Natural Science Foundation of China (60535010).     

Fast estimation of stable regions in real models

The combined use of periodic orbit computation and Fast Lyapunov Indicator Maps of systems with three degrees of freedom is proposed as a tool for estimating qualitative stability in corresponding real models. Application to spacecraft mission analysis reveals very useful in finding real orbits with very long term stability properties.     

纯引力轨道检验质量的相对测量技术

纯引力轨道是物体在太空仅受引力作用的运行轨道, 通过构造纯引 力轨道, 可实现超高精度的空间引力探测, 也可为科学实验提供超稳定卫星 平台. 作为纯引力轨道构造的核心, 检验质量的相对测量不仅提供了部分任 务科学数据, 还为航天器平台的跟踪控制提供输入. 首先, 描述了纯引力轨道 的概念内涵, 总结了它在卫星重力测量、引力波探测等方面的应用情况. 其 次, 综述了不同任务对相对测量的需求, 给出了电容式测量、磁感应测量和 光学测量的原理, 总结了各自的优缺点. 根据检验质量的姿态运动, 将检验质 量质心相对状态解算问题分为3 类, 给出了基于检验质量姿态动力学与表面 建模的典型解算模型和质心速度估计方法. 最后分析了非引力干扰的理论计 算、地面实验验证和在轨实验验证问题.      

陀螺进动与强迫进动轨道

陀螺进动与轨道进动现象的相似性归因于二者的动力学相似性. 通过类比二者的动力学模型,提出了一类强迫进动轨道. 若以圆轨道为初始轨道,通过施加常值法向力可以实现一种特殊的悬浮型强迫进动轨道. 采用四元数建模方法求解了这种强迫进动轨道的进动规律,给出了解析表达式,据此分析了这种轨道的性质. 分析结果表明这种强迫进动轨道与初始圆轨道在同一球面上,且与初始位置相切. 其角速度为进动角速度与初始轨道角速度的合成,是一种悬浮轨道,即属于非开普勒轨道. 悬浮轨道在地球观测、行星际科学、天文观测、无线电通讯以及地球工程等领域具有潜在应用前景. 从强迫进动的角度出发所作的分析为悬浮轨道的实现提供了一种新途径.     

The geometry of the orbits of artificial satellites

A novel approach to the study of the orbits of artificial satellites is presented. Emphasis is placed upon the basic geometry and other aspects of satellite motion which are of first importance to satellite engineering. The motion of the orbital plane as a rigid body is introduced and a non-elliptical orbit motion in this plane is defined. The plane orbit so defined possesses the very desirable feature of representing a succession of satellite positions and hence reveals the true motion of the satellite. An analytical treatment yields a completely general second order theory of earth satellite motion which is suitable for engineering purposes. In the latter development, particular attention is paid to the apsidal motion of the orbit and the concomitant resonance effects at the critical orbit inclination. The basic nonlinear features of the apsidal motion, which have not been recognized in earlier theories, are incorporated in the analytical development so as to produce a theory valid at all angles of inclination of the orbit.     

Coupled trajectory and attitude stability of displaced orbits

Coupled trajectory and attitude stability of displaced solar orbits is studied by using sailcraft with a kind of two-folding construction with two unequal rectangular plates forming a right angle. Three-dimensional coupled trajectory and attitude equations are developed for the coupled dynamical system, and the results show that all three types of displaced solar orbits widely referenced can be achieved through selecting an appropriate size of the two-folding sail. An anal- ysis of the corresponding linear stability of the trajectory and attitude coupled system is carried out, and both trajectory and attitude linearly stable orbits are found to exist in a small range of parameters, whose non-linear stability is then examined via numerical simulations. Finally, passively stable orbits are found to have weak stability, and such passive means of station-keeping are attractive and useful in practice because of its simplicity.     

THE EVOLUTION CHARACTERIZATION OF LIBRATIN POINT ORBITS AND APPLICATION RESEARCH

三体问题中, 轨道的受力和运动规律非常复杂. 对于特定的任务, 如何选择轨道的初始解是一大难题.针对平面三体问题, 利 用近拱点庞加莱映射, 对平动点顺行轨道和逆行轨道的长期和短期演化性质进行分析.根据轨道的初始状态将其分为逃逸轨道和捕获轨道.对于逃逸轨道, 给出了同宿轨道和异宿轨道的设计方法, 并利用两级微分修正法消除了拼接点处的位置不连续问题.对于捕获轨道, 得到了几类典型的周期和准周期轨道.对逆行轨道的演化性质进行分析时发现, 逆行轨道通常为准周期轨道, 比顺行轨道更加稳定.利用近拱点庞加莱映射可以快速确定不同类型轨道对应的初始状态, 为特定任务需求下的轨道设计提供了一种快速而有效的选择方案.     

Trajectory optimization and applications using high performance solar sails

The high performance solar sail can enable fast missions to the outer solar system and produce exotic non-Keplerian orbits. As there is no fuel consumption, mission trajectories for solar sail spacecraft are typically optimized with respect to flight time. Several investigations focused on interstellar probe missions have been made, including optimal methods and new objective functions. Two modes of interstellar mission trajectories, namely “direct flyby” and “angular momentum reversal trajectory”, are compared and discussed. As a foundation, a 3D non-dimensional dynamic model for an ideal plane solar sail is introduced as well as an optimal control framework. A newly found periodic double angular momentum reversal trajectory is presented, and some properties and potential applications of this kind of inverse orbits are illustrated. The method how to achieve the minimum periodic inverse orbit is also briefly elucidated.     

The Fast Norm Vector Indicator (FNVI) method: a new dynamical parameter for detecting order and chaos in Hamiltonian systems

In the present article, we introduce and also deploy a new, simple, very fast, and efficient method, the Fast Norm Vector Indicator (FNVI) in order to distinguish rapidly and with certainty between ordered and chaotic motion in Hamiltonian systems. This distinction is based on the different behavior of the FNVI for the two cases: the indicator after a very short transient period of fluctuation displays a nearly constant value for regular orbits, while it continues to fluctuate significantly for chaotic orbits. In order to quantify the results obtained by the FNVI method, we establish the dFNVI, which is the quantified numerical version of the FNVI. A thorough study of the method??s ability to achieve an early and clear detection of an orbit??s behavior is presented both in two and three degrees of freedom (2D and 3D) Hamiltonians. Exploiting the advantages of the dFNVI method, we demonstrate how one can rapidly identify even tiny regions of order or chaos in the phase space of Hamiltonian systems. The new method can also be applied in order to follow the time evolution of sticky orbits. A detailed comparison between the new FNVI method and some other well-known dynamical methods of chaos detection reveals the great efficiency and the leading role of this new dynamical indicator.     

Control of a Hyperchaotic Discrete System

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓｍｕｃｈａｔｔｅｎｔｉｏｎｈａｓｂｅｅｎｐａｉｄｔｏｃｏｎｔｒｏｌｌｉｎｇｃｈａｏｓｂｅｃａｕｓｅｏｆｉｔｓｔｈｅｏｒｅｔｉｃｉｍｐｏｒｔａｎｃｅａｎｄｐｏｓｓｉｂｌｅａｐｐｌｉｃａｔｉｏｎｓ[1,2 ].Ｃｏｎｔｒｏｌｏｆｈｙｐｅｒｃｈａｏｓｈａｓｂｅｅｎｓｔｕｄｉｅｄａｓｗｅｌｌ[3].Ａｍｏｎｇｓｅｖｅｒａｌａｐｐｒｏａｃｈｅｓ,ｔｈｅＯＧＹｍｅｔｈｏｄ[4 ]ｐｒｏｐｏｓｅｄｂｙＯｔｔｅｔａｌｉｓｖｅｒｙａｔｔｒａｃｔｉｖｅ .Ｉｎｔｈｅｍｅｔｈｏｄ ,ｔｈ…