Trajectory optimization of multiple asteroids exploration with asteroid 2010TK7 as main target

In this study, the Earth’s Trojan asteroid 2010 TK₇ is selected as the rendezvous target. The multiple flyby sequence of asteroid exploration was proposed by optimizing the probe’s orbit. Impulsive maneuvers and low-thrust propulsion were used respectively to design the trajectories of the multiple asteroids exploration mission. Under impulsive maneuvers, gravity assist technique was adopted to reduce fuel consumption. First a reference orbit with only 2010 TK₇ as the rendezvous target was designed. Then five asteroids near the reference orbit were selected as candidates. Finally, we obtained a multiple asteroids exploration sequence of three asteroids based on gravity assist technique and genetic algorithm, and an additional velocity impulse of 0.4?km/s was required. In the subsequent section, a sixth-degree inverse polynomial shape-based method is applied to the low-thrust trajectory design of 2010 TK₇, and the exploration sequence under the action of low-thrust propulsion was provided.     

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

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

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

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

双体小行星系统平衡态与稳定性研究

由于双体小行星独特的运动形式可为行星的演化提供重要线索,因而成为小行星探测中的热点目标。基于双椭球体模型研究了双小行星系统的相对运动、平衡态及稳定性。首先基于双椭球的全二体模型建立了双星系统相对运动的动力学方程;其次利用拉格朗日方程,通过求解系统角动量和能量,建立了双星系统平衡态对应的状态约束;最后给出了通过零速度状态曲面判断双星系统平衡态稳定性的一般性方法,在此基础上分析了小行星物理参数变化对系统平衡态稳定性的影响。研究可为未来双体小行星系统探测任务中的轨道设计与控制提供重要的理论参考。     

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

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

近地小行星采矿与防御计划发展现状

近地小行星交会、绕飞、着陆与采样返回技术经过数10年的发展日趋成熟。美国的OSIRIS-Rex对C类小行星进行特征分析与采样,日本宇宙航空研究开发机构的“隼鸟-2号”任务目的是小行星深层采样。美国国家航空航天局（NASA）和欧洲空间局（ESA）的小行星探测任务开始转向行星防御领域。NASA的ARM（Asteroid Redirect Mission）计划是开展小行星抓捕与轨道重定向,ESA联合NASA提出了小行星撞击与偏转评估计划,拟对双星系统开展撞击实验,为行星防御提供技术积累。此外,行星资源公司和深空工业公司分别规划了小行星商业采矿的蓝图,并已开展相关的在轨技术验证。对近地小行星的探测历程进行了回顾,重点介绍了OSIRIS-Rex、“隼鸟-2号”、NASA和ESA的行星防御计划及小行星采矿公司的商业采矿战略规划,总结了未来开展行星防御与采矿的关键技术。     

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.     

小行星远距离抵近轨道

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

A mission concept for the low-cost large-scale exploration and characterisation of near earth objects

This paper presents the preliminary mission and science analysis of a new mission concept for the large scale, low-cost exploration of Near Earth Objects (NEOs). The concept is to enable close range observations of NEOs by performing close flybys of a series of NEOs at one of their nodal points, with pairs of small spacecraft flying in formation. The paper presents a preliminary assessment of accessible asteroids and multi-target tour trajectories from data available in the JPL small-body database.The main instruments on board each spacecraft are a camera and a LIDAR which together can be used for orbit determination, surface imaging, direct asteroid ranging and asteroid mass estimation via intersatellite ranging. The paper provides a qualitative and quantitative assessment of the measurable quantities during each flyby. In particular, the feasibility of a novel method of NEO mass estimation is assessed.     

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.     

小行星探测电推进系统方案研究

小行星探测以及资源开发与利用对国家抢占深空探测主动权和制高点有着不可估量的战略意义。电推进具有高比冲、长寿命和高度自主巡航等特点,小行星探测器采用电推进执行巡航阶段轨道机动任务,将大幅减少推进剂重量和提高载重比。调研了国外小行星探测的电推进系统方案,针对我国小行星探测对电推进系统的任务需求及现有电推力器的技术基础,提出了5种电推进系统方案,并进行多维度对比,对最优方案进行了设计和关键技术梳理。     

从日地系统L2出发借力月球飞越近地小行星

对于停留在日地系统L2的“嫦娥2号”探测器,其后续飞行方案有多个选项,例如主动撞月或重返月球轨道、返回地球轨道或再入大气、飞往地月系统L1/L2或日地系统L1、进入深空飞越近地小行星(最终,“嫦娥2号”于2012年12月13日成功地实现了对Toutatis小行星的近距离飞越)。探讨上述的飞行方案需要对飞行轨道进行初步设计,总的速度脉冲限制在100 m/s以内并且需要考虑探测器同时受到太阳、地球、月球的引力作用。本研究设计了探测器从日地系统L2出发借力月球实现Toutatis小行星飞越的飞行方案,与直接飞越方案相比,借力月球可以进一步节省探测器的燃料消耗,其等效速度脉冲设计值为58.47 m/s。     

Attitude stability and periodic attitudes of rigid spacecrafts on the stationary orbits around asteroid 216 Kleopatra

In this work, equilibrium attitude configurations, attitude stability and periodic attitude families are investigated for rigid spacecrafts moving on stationary orbits around asteroid 216 Kleopatra. The polyhedral approach is adopted to formulate the equations of rotational motion. In this dynamical model, six equilibrium attitude configurations with non-zero Euler angles are identified for a spacecraft moving on each stationary orbit. Then the linearized equations of attitude motion at equilibrium attitudes are derived. Based on the linear system, the necessary conditions of stability of equilibrium attitudes are provided, and stability domains on the spacecraft’s characteristic plane are obtained. It is found that the stability domains are distributed in the first and third quadrants of the characteristic plane and the stability domain in the third quadrant is separated into two regions by an unstable belt. Subsequently, we present the linear solution around a stable equilibrium attitude point, indicating that there are three types of elemental periodic attitudes. By means of numerical approaches, three fundamental families of periodic solutions are determined in the full attitude model.     

Optimal station-keeping near Earth–Moon collinear libration points using continuous and impulsive maneuvers

In this study the gravitational perturbations of the Sun and other planets are modeled on the dynamics near the Earth–Moon Lagrange points and optimal continuous and discrete station-keeping maneuvers are found to maintain spacecraft about these points. The most critical perturbation effect near the L₁ and L₂ Lagrange points of the Earth–Moon is the ellipticity of the Moon’s orbit and the Sun’s gravity, respectively. These perturbations deviate the spacecraft from its nominal orbit and have been modeled through a restricted five-body problem (R5BP) formulation compatible with circular restricted three-body problem (CR3BP). The continuous control or impulsive maneuvers can compensate the deviation and keep the spacecraft on the closed orbit about the Lagrange point. The continuous control has been computed using linear quadratic regulator (LQR) and is compared with nonlinear programming (NP). The multiple shooting (MS) has been used for the computation of impulsive maneuvers to keep the trajectory closed and subsequently an optimized MS (OMS) method and multiple impulses optimization (MIO) method have been introduced, which minimize the summation of multiple impulses. In these two methods the spacecraft is allowed to deviate from the nominal orbit; however, the spacecraft trajectory should close itself. In this manner, some closed or nearly closed trajectories around the Earth–Moon Lagrange points are found that need almost zero station-keeping maneuver.     

潜在威胁小行星碰撞防御的计算与分析

研究了采用碰撞的方式进行小行星防御的动力学问题。采用多面体模型来建立小行星的外形模型,以碎石堆模型来建立小行星的结构模型,计算了小行星受到与其密度和材质相同的球体高速碰撞过程和碰撞后的碎石分布。计算过程在考虑了小行星与碰撞球体的接触形变以及小行星内部组成碎石堆的接触形变条件下,计算了碎石堆内部的相互引力、法向接触力、切向静摩擦力、切向动摩擦力和滚动摩擦力矩。以小行星101955 Bennu(中文名贝努)为对象计算了潜在威胁小行星的碰撞防御过程的动力学行为。结果显示:采用高速碰撞的方法进行小行星防御可以有效地将小行星撞成大量碎小的石块,且该方法具有核爆的方法不可比拟的优势,即对空间环境无污染。     

利用绳系太阳帆减缓小行星自转的技术研究

提出一种利用太阳帆绳系系统逐渐减小小行星自转速率的方法。在系绳长度不变的情况下,利用太阳帆受到的太阳光压力使其始终保持与小行星同步,避免了由于小行星自转而引起的系绳缠绕问题。通过控制太阳帆使系绳始终拉紧,系绳中的拉力便可以持续提供一个与小行星自转方向相反的力矩,从而减小其自转速率。仿真结果表明,面积10⁶ m²的太阳帆,经过约86天可将小行星的自转消除,验证了该方法的有效性。     

Collinear equilibrium points in the elliptic R3BP with oblateness and radiation

We have studied the positions and stability of the collinear equilibrium points, L_1,2,3 of an infinitesimal body in the elliptic restricted three-body problem (ER3BP) when both primaries of the system are luminous and oblate spheroids moving in elliptic orbits around their common center of mass. We observe that their positions are affected by the radiation pressure forces, oblateness and the eccentricity of the orbits, but the stability character remains unchanged and are unstable. The effects of the parameters involved on the collinear points, in particular for the binary systems Achird, Luyten 726-8, Kruger 60, Alpha Centauri AB and Xi Bootis, and their stability in general have been investigated numerically using the analytical results obtained.     

小行星附近制导与控制研究综述

由于小行星具有引力场不规则、物理参数不确定性大、表面逃逸速度小等因素,使得小行星附近制导和控制极具挑战性。回顾了小行星探测任务的历史、现状和意义;针对小行星附近制导与控制研究的基础即小行星附近动力学,分析了研究现状;针对悬停、绕飞、转移、着陆等任务形式,详细介绍了轨道制导和控制方法方面的热点问题和研究现状;基于研究现状,列举了部分未来可进一步研究的方向和问题。     

Analytic and numerical treatment of motion of dust grain particle around triangular equilibrium points with post-AGB binary star and disc

This paper investigates the motion around the triangular equilibrium points, of a passively gravitating dust particle in the gravitational field of a low-mass post-AGB binary system, surrounded by circumbinary disc. The two bodies of the binary are modeled as a triaxial star and a radiating-oblate star. Due to small deviation of disc stars on circular orbits, we have assumed that the Coriolis and centrifugal forces of the stars are slightly perturbed. The triangular equilibrium points of the particle are found. These points are defined by, triaxiality of the primary star, oblateness and radiation of the secondary one and the gravitational potential from the disc mass. Further, when the disc mass increases, the particle moves nearer to the stars and farther away from the disc. In general, these equilibrium points are linearly stable when μ < μ_C; where μ is the mass ratio and μ_C is the critical mass function, defined by the parameters of the system. The effects of each of these parameters on the size of the stability region are stated, and the periodic motion around the stable points is examined. It is seen that the orbits are ellipses, and the orientation, eccentricities, lengths of the semi-major and semi-minor axes are influenced by the parameters of the problem. In particular, for our numerical linear stability analysis, we have taken an extremely depleted pulsating star, IRAS 11472-0800 as the post-AGB triaxial star, with a weakly-radiating young white dwarf star; G29-38 as the secondary. For this system, the stability result of the triangular points comes out different. Here, μ_C < μ throughout the entire range of the mass ratio and the critical mass function. Hence, the triangular equilibrium points are unstable. The stability of the orbits is tested using the Poincaré surfaces of section (Pss). The region of stability is controlled by the introduced parameters and the Jacobi constant.     

Modeling irregular small bodies gravity field via extreme learning machines and Bayesian optimization

Close proximity operations around small bodies are extremely challenging due to their uncertain dynamical environment. Autonomous guidance and navigation around small bodies require fast and accurate modeling of the gravitational field for potential on-board computation. In this paper, we investigate a model-based, data-driven approach to compute and predict the gravitational acceleration around irregular small bodies. More specifically, we employ Extreme Learning Machine (ELM) theories to design, train and validate Single-Layer Feedforward Networks (SLFN) capable of learning the relationship between the spacecraft position and the gravitational acceleration. ELM-base neural networks are trained without iterative tuning therefore dramatically reducing the training time. Analysis of performance in constant density models for asteroid 25143 Itokawa and comet 67/P Churyumov-Gerasimenko show that ELM-based SLFN are able learn the desired functional relationship both globally and in selected localized areas near the surface. The latter results in a robust neural algorithm for on-board, real-time calculation of the gravity field needed for guidance and control in close-proximity operations near the asteroid surface.