小天体着陆动力学参数不确定性影响分析

针对小天体不规则程度高、引力场复杂,且物理参数存在较高不确定性的问题,基于小天体着陆动力学方程线性化近似解析解,对各动力学参数不确定性的影响进行了分析。考虑动力学方程线性化带来的误差,引入线性化误差补偿校正方法,建立了探测器轨迹对动力学参数不确定性的敏感度方程。以小行星Eros 433为例,重点分析了目标小天体质量、自转角速度、引力势函数系数,以及探测器初始状态、推力加速度等动力学参数不确定性对探测器着陆轨迹的影响。数学仿真分析表明,针对本文选取的目标小天体,推力加速度扰动为主要影响因素,探测器初始状态的不确定性为次要影响因素,其他参数扰动的影响较小。     

小行星远距离抵近轨道

讨论了小行星引力一阶项可被忽略情况下的小行星远距离轨道设计及动力学。此时,航天器的运动受太阳引力和太阳光压的影响。航天器和小行星的加速度之差在这两者之间形成的独特的相对动力学,为航天器在小行星附近停驻与观测提供特定轨道。完整解决了小行星处于圆形日心轨道这一较简单情况,也考虑和阐述了椭圆轨道情况,并取得了一些初步结果。     

A fast Chebyshev polynomial method for calculating asteroid gravitational fields using space partitioning and cosine sampling

This paper presents a computationally fast method for solving gravitational accelerations near irregularly-shaped asteroids. This method is based on analytical three-dimensional Chebyshev polynomial approximation of the polyhedral gravity. For the purpose of improving the approximation accuracy, space partitioning schemes based on practical flight zones is used to avoid interpolation the whole space around the target asteroid. Specifically, a minimum ellipsoid close to the asteroid surface is defined to select the space for surrounding trajectories with safe distance and a cone connected to the surface is defined to select the space for descent trajectories. Moreover, interpolation points are sampled in a cosine sampling fashion according to the Chebyshev-Gauss-Lobatto nodes and a radial adaption technique. The performance of different space partitioning schemes is analyzed. The effectiveness of the proposed method is validated through simulations of solving gravitational accelerations at the test points near different shaped asteroids 1996 HW1, 433 Eros, 25143 Itokawa and 101955 Bennu.     

On tethered sample and mooring systems near irregular asteroids

A tethered asteroid sample and mooring system is investigated in this paper. In this system the spacecraft is moored to the surface of an irregular asteroid such as 216 Kleopatra by using a rocket-propelled anchor with a cable. The rocket-propelled anchor is a kind of space penetrator, which can inject into asteroids at high speeds generated by its own rocket engine. It can be used to explore the interior structure of asteroids, and it can also be used as a sample collector. When the sampling mission is done, the sample can be pulled back to the spacecraft with the anchor. Using this method, the spacecraft can be kept in a safe region in which it cannot be trapped by the gravitational field of the asteroid. This work is concerned with the dynamics of the tethered system near irregular asteroids. First, a shape model and gravitational field model of irregular asteroids are built. Then, the configuration and the stability of the tethered system are investigated, and the quasi-periodic motion near the equilibrium point of the tethered system is analyzed. Finally, the non-uniform density distribution of the asteroids is considered. The deployment process and the oscillation of the tethered system in the uncertain asteroid gravity field are simulated using the Monte Carlo method. The feasibility of the tethered asteroid sample and mooring system is proved.     

强不规则天体引力场中的动力学研究进展

小行星探测与彗星探测是深空探测的重要方面。一般来说,小行星和彗星因质量都不足以使得万有引力克服应力达到流体静力学平衡,而具有强不规则的外形。研究强不规则天体引力场中的动力学行为及其内在机制,是探测器被不规则天体捕获并对其形成近距离探测轨道的基础。从引力场模型和动力学行为两个方面综述了强不规则天体引力场中动力学的研究进展,在引力场模型的研究方面介绍了强不规则天体引力场建模的球谐函数摄动展开模型、简单特殊体模型及多面体模型的研究现状,在动力学机制的研究方面介绍了强不规则天体引力场中的周期轨道和拟周期轨道、平衡点、流形、分岔与共振以及混沌运动的研究现状,指出了这些方面研究的重点与难点。分析了强不规则体引力场中动力学的研究趋势。     

Orbital stability close to asteroid 624 Hektor using the polyhedral model

We investigate the orbital stability close to the unique L4-point Jupiter binary Trojan asteroid 624 Hektor. The gravitational potential of 624 Hektor is calculated using the polyhedron model with observational data of 2038 faces and 1021 vertexes. Previous studies have presented three different density values for 624 Hektor. The equilibrium points in the gravitational potential of 624 Hektor with different density values have been studied in detail. There are five equilibrium points in the gravitational potential of 624 Hektor no matter the density value. The positions, Jacobian, eigenvalues, topological cases, stability, as well as the Hessian matrix of the equilibrium points are investigated. For the three different density values the number, topological cases, and the stability of the equilibrium points with different density values are the same. However, the positions of the equilibrium points vary with the density value of the asteroid 624 Hektor. The outer equilibrium points move away from the asteroid’s mass center when the density increases, and the inner equilibrium point moves close to the asteroid’s mass center when the density increases. There exist unstable periodic orbits near the surface of 624 Hektor. We calculated an orbit near the primary’s equatorial plane of this binary Trojan asteroid; the results indicate that the orbit remains stable after 28.8375?d.     

Phobos-Grunt: Russian sample return mission

As an important milestone in the exploration of Mars and small bodies, a new generation space vehicle “Phobos-Grunt” is planned to be launched by the Russian Aviation and Space Agency. The project is optimized around a Phobos sample return mission and follow up missions targeted to study some main asteroid belt bodies, NEOs and short period comets. The principal constraint is use of the “Soyuz-Fregat” rather than the “Proton” launcher to accomplish these challenging goals. The vehicle design incorporates innovative SEP technology involving electrojet engines that allowed us to increase significantly the mission's energetic capabilities, as well as highly autonomous on-board systems. Basic criteria underlining the “Phobos-Grunt” mission scenario, scientific objectives and rationale including Mars observations during the vehicle's insertion into Mars orbit and Phobos approach maneuvers, are discussed and an opportunity for international cooperation is suggested.     

双小行星系统引力场建模方法研究

针对弱引力双小行星系统的引力场建模问题,本文采用复杂度和精度依次递增的球体–球体模型、椭球体–球体模型和改进的限制性椭球体–椭球体模型来进行引力场建模,并分别采用椭圆积分以及无积分环节、计算效率高的二阶二次球谐函数来表征引力势,从而比较精确地刻画双小行星系统和探测器构成的限制性全三体问题的动力学模型;针对双小行星系统1999KW4,对其不同的引力场模型进行了仿真研究,分别给出了不同模型下的等效势能函数曲面及零速度曲线,比较了不同模型下的平动点位置坐标偏差。结果表明,二阶二次球谐函数计算引力势的椭球体-椭球体模型计算精度高,复杂程度低,计算量更少,计算速度更快,能够较精确的对双小行星系统进行引力场建模。     

Finite-time asteroid hovering via multiple-overlapping-horizon multiple-model predictive control

This paper investigates the asteroid hovering problem using the Multiple-Overlapping-Horizon Multiple-Model Predictive Control method. The effectiveness of the predictive controllers in satisfying control constraints and minimizing the required control effort is making Model Predictive Control a desirable control method for asteroid exploration missions which consist of the asteroid hovering phase. However, the computational burden of Model Predictive Control is an obstacle to employing the asteroid’s complex gravitational field model. As an alternative option, the Multiple Horizon Multiple-Model Predictive Control method has been introduced previously, which could provide a solution with the less computational burden with respect to the nonlinear Model Predictive Control. It was shown that it is not necessary to deduce the exact dynamics model to predict the system’s behavior during a long period using this approach. However, the calculated control acceleration was not smooth enough because of the crisp borders of consecutive horizons, which may cause an image motion and degrades the geometric accuracy of high-resolution images in asteroid hovering missions. In this paper, the Multiple-Overlapping-Horizon Multiple-Model Predictive Control method is introduced instead to solve the problem of controlling acceleration fluctuations by overlapping consecutive horizons. Numerical simulation results are presented to validate the effectiveness of the proposed control method, and its advantage is demonstrated accordingly for the asteroid hovering problem in achieving the hovering position and velocity.     

Landing of hopping rovers on Irregularly-shaped small bodies using attitude control

A hopping rover that is driven only by internal or external attitude actuators is an ideal mobility approach for surface exploration of small solar system bodies. Without thrust control and grasping mechanisms, a hopping rover is mechanically simple to design and less prone to mechanical failures, but faces challenges during soft landing. It may rebound from the surface, causing deviations from its original landing site. In this paper, landing of a hopping rover on the surface of an asteroid is investigated, and a strategy using only attitude control to shorten the landing distance is proposed. Based on rigid body impact dynamics, the edge impact configuration is investigated in detail. The factors that affect the impact states of a cube-shaped hopping rover are studied. Then, controlled edge landing is analyzed, in which the post-impact velocity of the hopping rover is changed by controlling its attitude prior to impact. Three guidance schemes are developed, followed by attitude profile generation and finite time stable attitude control. Finally, simulations are performed on an ideal flat surface and uneven terrain. The results show that controlled edge landing can effectively reduce the landing distance and settling time, compared with uncontrolled landing. This study on hopping motion on the surface of an irregular-shaped asteroid with attitude control, can provide a reference for hopper trajectory plan in future asteroid surface explorations.     

Adaptive double-saturated control for hovering over an asteroid

Hovering over an irregular-shaped asteroid is particularly challenging due to the large gravitational uncertainties and various external disturbances. An adaptive control scheme considering commanded acceleration and its change-rate saturation for hovering is developed in this paper. Taking full advantage of terminal sliding-mode control theory, first, we convert the double-saturated control problem to a new equivalent system by introducing a special bounded function, in which just control input saturation needs to be considered. Then, a continuous finite-time saturated controller is designed for the new system with the assistance of an constructed auxiliary subsystem. Additionally, an adaptive law is devised for the controller to avoid the requirement of the unknown upper bounds of the disturbances, rendering the control scheme especially suitable to asteroid hovering missions. The finite-time stability of the whole closed-loop system is proved via Lyapunov analysis. Numerical simulation studies are carried out, and the results demonstrate the design features and the desired performance.     

Study of highly perturbed spacecraft formation dynamics via approximation

This paper explores methods for approximating and analyzing the dynamics of highly perturbed spacecraft formations with an emphasis on computationally efficient approaches. This facilitates on-board computation or rapid preliminary mission design analysis. Perturbed formation dynamics are often approximated as linear time-varying (LTV) systems, for which Floquet theory can be used to analyze the degree of system instability. Furthermore, the angular momentum of the relative orbital state can be computed with the approximate dynamics to provide additional insight. A general methodology is developed first and then applied to the problem of unstable formation dynamics in asteroid orbits. Here the dominant perturbative effects due to low-order gravitational harmonics and solar radiation pressure are modeled. Numerical simulations validate the approach and illustrate the approximation accuracy achieved.     

附着小天体的动态面鲁棒制导与控制方法

以精确附着小天体表面的任务为背景,提出一种基于扰动观测器(DOB)和动态面控制的附着小天体的制导与控制方法。根据探测器的初始条件与终端着陆条件规划了标称轨迹,并将引力场建模误差、参数摄动和外部干扰等视为总扰动,结合动态面控制和DOB设计了标称轨迹跟踪控制器。分析总扰动估计误差的渐进收敛性以及闭环标称轨迹跟踪控制系统的稳定性,并确定控制器参数选取条件。数值仿真结果表明,所设计的DOB可以有效地估计并抑制总扰动且闭环标称轨迹跟踪控制系统具有良好的稳定性和控制精度。     

附着小行星的视线制导规律

提出了一种航天器附着（软着陆）小行星的制导与控制规律。为了保证垂直软着陆，事先规划了满足约束的理想视线与视线角轨迹，通过设计滑模变结构控制器跟踪理想轨迹，实现了在小行星表面垂直软着陆；证明了制导与控制规律的鲁棒性；最后通过数学仿真验证了提出的制导与控制规律的可行性。     

小天体自主附着多滑模面鲁棒制导方法研究

小天体形状不规则及缺乏观测信息的特点使得小天体附近的动力学环境较为复杂,附着动力学模型存在较大不确定性。通过引入多滑模面鲁棒制导方法,分别设计2个滑模面,使探测器状态先后到达这2个滑模面,可实现指定时刻精确附着小天体的目标。通过选取参数的分析总结了制导律中相关参数的选取对燃料消耗的影响,给出了制导律相关参数选取原则。在存在外界环境扰动、初始状态误差和导航误差条件下,蒙特卡洛仿真结果表明:多滑面制导方法能够在小天体的不确知环境中实现高精度附着,且具有很好的鲁棒性。多滑模面制导方法精度高、鲁棒性好,且无需设计参考轨迹,实时性好,适合小天体自主精确附着的任务需求。     

小天体安全着陆与表面探测控制方法研究

针对探测器在小天体表面的弹跳式运动特点,提出一种表面运动控制方法。首先基于简化的运动模型,计算出探测器进行弹跳式转移所需的初始速度;然后通过对轨道方程的线性化,基于线性时变系统的求解方程,采用脉冲控制对每次弹跳的运动轨迹进行修正,增加控制精度;并分析转移距离与弹跳次数的关系;最后,以小行星1620 Geographos作为目标小天体进行仿真分析,验证了控制算法的有效性。结果表明,该控制方法可以应用于探测器在附着到小天体表面后的弹跳转移,也适用于弱引力环境下探测器的弹跳式探测。     

三小行星系统216 Kleopatra引力场中的动力学

主带三小行星系统216 Kleopatra是由主星216 Kleopatra及两个小月亮(moonlet)Alexhelios[S/2008(216)1]和Cleoselene[S/2008(216)2]组成。其中主星216Kleopatra是一个具有强不规则形状如哑铃的连接双星,大小为217km×94km×81km,外小月亮Alexhelios大小约为8.9km,内小月亮大小约为6.9km。其动力学行为具有非常丰富的科学内涵。研究了三小行星系统216Kleopatra自身的动力学机制及其引力场中探测器的运动规律,分析了主星质心固连系中探测器的动力学方程,给出了三小行星引力全多体问题的动力学方程及Jacobi积分,方程考虑了三个小行星的不规则外形、轨道与姿态。发现三小行星系统216Kleopatra主星引力场中一种新的周期轨道族的倍周期分岔。考虑主星的不规则精确外形与引力、两个小月亮的相互作用,研究了该三小行星系统的动力学构形。发现Kleopatra的强不规则几何外形及两个小月亮Alexhelios和Cleoselene的相互作用引起两个小月亮的轨道参数的显著变化。     

小行星探测发展综述

小天体上保存着太阳系形成初期的原始成分,同时可能蕴含着地球生命与水起源的重要线索,是研究太阳系起源和演化历史的活化石。近年来,小行星探测已成为主要航天国家深空探索领域的重点发展目标之一。简要总结了国际上小天体探测历程,对小行星探测的研究和发展趋势进行了综述,重点探讨了未来小行星探测任务面临的主要关键技术,并对中国后续开展小行星探测活动提出了相关建议。     

双星系统L1点悬停探测控制器设计与仿真

双小行星系统探测具有重要的科学意义,受其复杂动力学环境影响,探测任务极具挑战。利用球谐函数法对双星系统进行引力场建模,求解双星系统平动点,并选取其内部共线平动点L_1点作为双星系统悬停探测目标位置。采用航天探测实际任务中常用的脉冲推力式发动机,设计了一种原理简单、便于工程实现的常值切换bang-bang控制器。以69230Hermes双星系统为例,将Hermes近似为双椭球系统,仿真分析航天器在Hermes双星系统L_1点悬停飞行的控制效果,验证所提控制策略有效性。     

VLBI for better gravimetry in SELENE

The Japanese lunar explorer SELENE (SElenological and Engineering Explorer), to be launched in 2007, will for the first time utilize VLBI observations in lunar gravimetry investigations. This will particularly improve the accuracy to which the low degree gravitational harmonics and the gravity field near the limb can be measured, and when combined with Doppler measurements will enable three-dimensional information to be extracted. Differential VLBI Radio sources called VRAD experiment involves two on-board sub-satellites, Rstar and Vstar. These will be observed using differential VLBI to measure the trajectories of the satellites with the Japanese network named VERA (VLBI Exploration of Radio Astrometry) and an international VLBI network.