2. 深空自主导航与控制工信部重点实验室, 北京 100081;
3. 飞行器动力学与控制教育部重点实验室, 北京 100081
2. Key Laboratory of Autonomous Navigation and Control for Deep Space Exploration, Ministry of Industry and Information Technology, Beijing 100081, China;
3. Key Laboratory of Dynamics and Control of Flight Vehicle, Ministry of Education, Beijing 100081, China
小天体是太空环境中具有极高探测价值的一类天体。附着探测是提升人类对小天体认知的重要手段。相 比于大的行星天体,小天体具有质量轻、形状不规则的特点,这导致其附近呈现不规则弱引力场[1]。在这种引力场环境下,探测器在小天体附近的运动形式可以概括为:贴近小天体表面的轨道段、轨道运动与表面作用耦合的弹跳段、完全位于表面上的滑移滚动段。
针对小天体附近独特的运动形式,相关学者对其动力学特性和稳定性展开了研究。Guibout等构建了一个均匀旋转的均质三轴椭球模型,采用解析算法得到了一组平衡点,并对其稳定性进行了全面的分析[2]。Liu等采用一个旋转的均质立方体作为小天体模型,给出了忽略摩擦情况下立方体表面的平衡点和围绕平衡点的周期轨道,以及有摩擦情况下质点不从表面逃逸的条件[3]。Yu等基于多面体模型对粒子在Geographos小行星表面的运动行为进行了归纳和分析,并修正了多面体模型边界处的引力势函数奇异问题[4]。Tardivel等构建了球形探测器在不规则小行星表面的受力模型,并运用概率模型描述了小行星表面的障碍物情况[5]。
在小天体附近的弱引力环境下,任何微小的扰动都可能导致探测器离开小天体表面。在完成小天体附着时,为了避免反弹现象的发生,相关任务对着陆装置的构型都做了特殊的设计。JAXA的Hayabusa1任务中,着陆信标被设计成一个由若干小球填充的袋子,在信标与目标小行星发生碰撞后,通过小球的相互碰撞,将能量消耗掉,使得信标有较小的恢复系数[6]。ESA罗塞塔任务的着陆器菲莱通过起落架上的阻尼装置消耗掉一部分碰撞能量,冷喷气系统提供向下的推力,另外,菲莱会释放两个鱼叉进入地表,使自身固定[7-8]。TAG(touch-and-go) 是弱引力环境下,探测器完成对小天体表面采样收集的有效手段。在TAG过程中,探测器通常先在目标小天体上方进行悬停。待确定采样点以后,探测器下降到小天体表面,在只做几秒钟的停留完成采样任务后迅速上升离开小天体表面。TAG采样已经成功应用于Hayabusa任务[9],OSIRIS- Rex[10]任务也将采用这种方式完成采样返回。
由于小天体附近引力微弱,表面不能提供足够的摩擦力,传统的轮动式巡视器无法有效运动,因此小天体表面巡视器通常也采用弹跳式方法。针对探测器在小天体表面的弹跳移动问题,Ulamec对小天体探测任务中弹跳式探测器的应用价值进行了深入研究,并且在基于旋转质量块加速弹跳原理的基础上提出了撑杆跳弹跳概念[11]。前苏联的Phobos 2任务中携带了一个计划着陆于火卫一的弹跳探测器PROP-F,虽然任务失败,但弹跳探测器的概念对小天体探测具有很高的参考价值[13]。Hayabusa任务中,MINERVA巡视器的驱动装置为扭矩电机,它可以提供小于10 cm/s的跳跃速度。Yoshimitsu基于MINERVA巡视器建立了小天体表面弹跳过程动力学,并分析了不同环境参数下弹跳性能与应用前景[15]。Bellerose针对小天体弹跳移动方法,基于椭球体引力模型,提出了探测器移动距离和所需时间的解析计算方法,并基于滑模控制思想研究了单个和多个探测器控制方法[14]。Mège等针对小天体复杂地形设计了对附着环境具备鲁棒性的弹跳运动机构,弹跳距离可以从几米到几十米[12]。
弱引力环境下进行附着任务极易发生反弹,如果没有后续操作,探测器将长时间的在小天体表面弹跳运动,造成附着任务的失败。为了避免这一现象的发生,提高小天体附着探测任务的稳定性,本文提出一种表面运动控制方法,使得探测器在附着到小天体表面后能够以弹跳的方式进行运动,最终稳定停留在原着陆点或转移到新的目标点。同时,这一控制方法也适用于弱引力环境下探测器的弹跳式探测。
1 动力学模型 1.1 多面体引力场模型本文采用小行星1620 Geographos作为目标小天体,其多面体模型如图 1所示。定义小天体固连坐标系如下:原点位于小天体质心,z轴与小天体最大惯量轴即自转轴重合,x与y轴分别与最小和中间惯量轴重合,x,y,z三轴满足右手法则。
![]() |
图 1小行星1620 Geographos多面体模型Fig. 1Polyhedral model of asteroid Geographos |
把小天体近似为内部密度恒定的多面体以后,可以将引力势函数写成如下的形式:
$V = - \frac{1}{2}G\rho \sum\limits_{e = 1}^{{N_e}} {{r_e} \cdot {E_e} \cdot {r_e} \cdot {L_e}} + \frac{1}{2}{\rm{G}}\rho \sum\limits_{f = 1}^{{N_f}} {{r_f} \cdot {F_f} \cdot {r_f} \cdot {\omega _f}} $ | (1) |
其中:G表示万有引力常数;ρ表示小天体密度;re和rf分别表示场点到多面体边和面上任意点的距离;Le的值定义如下
$ {L_e} = \ln \frac{{{a_e} + {b_e} + {e_e}}}{{{a_e} + {b_e} - {e_e}}} $ | (2) |
其中:ae和be分别表示场点到任意边两顶点的距离;ee为该边的长度;Ff为构成边e的两个面的外法线矢量的点积,即
$ {{{F} }_f} = {{{n}}_f} \cdot {{{n}}_f} $ | (3) |
Ee定义如下
$ {{{E}}_e} = {{{n}}_{f1}}{{n}}_e^{f_1} + {{{n}}_{f2}}{{n}}_e^{f_2}\;\;\;\;\;1\; \leqslant f_1,\;f_2 \leqslant {N_f} $ | (4) |
其中:${{n}}_e^f$表示面f和边e的外法线方向矢量。
ωf的值定义如下
$ {\omega _f} = 2\arctan \frac{{{{{r}}_1} \cdot {{{r}}_2} \times {{{r}}_3}}}{{1 + {{{r}}_1} \cdot {{{r}}_2} + {{{r}}_2} \cdot {{{r}}_3} + {{{r}}_3} \cdot {{{r}}_1}}} $ | (5) |
根据以上多面体引力场模型,能够计算场点处的引力势能,通过求引力势函数的一阶和二阶导数,可以得到场点处的引力加速度和Hessian矩阵:
$\left\{ {\begin{array}{*{20}{l}} {\nabla V{\rm{ = }}G\rho \sum\limits_e {{E_e} \cdot {r_e}} \cdot {L_e} - G\rho \sum\limits_f {{F_f} \cdot {r_f}} \cdot {\omega _f}}\\ {{\nabla ^2}V{\rm{ = }} - G\rho \sum\limits_e {{E_e}} \cdot {L_e}{\rm{ + }}G\rho \sum\limits_f {{F_f}} \cdot {\omega _f}} \end{array}} \right.$ | (6) |
探测器在小天体附近的运动包含轨道段和表面段,小天体固连坐标系下的轨道段运动方程为
$\left\{ {\begin{array}{*{20}{l}} {\dot x = {v_x}}\\ {\dot y = {v_y}}\\ {\dot z = {v_z}}\\ {{v_x} = 2\omega {v_y} + {\omega ^2}x + {V_x} + {a_x}}\\ {{v_y} = - 2\omega {v_x} + {\omega ^2}y + {V_y} + {a_y}}\\ {{v_z} = {V_z} + {a_z}} \end{array}} \right.$ | (7) |
其中:ω表示目标小天体的自旋角速度;Vx、Vy和Vz表示沿坐标轴方向的引力加速度分量,可以由公式(6) 计算得到;ax、ay和az表示沿坐标轴方向的控制加速度。
1.3 表面段动力学方程表面段运动过程可以分解为法向和切向,分解过程如图 2所示。由于表面段运动时间很短,这里忽略掉探测器的滑移过程,仅考虑碰撞。假设探测器与小天体表面接触前具有入射速度vc,碰撞后的法向速度由恢复系数cr(0≤cr≤ 1) 来确定。
![]() |
图 2表面段运动分解Fig. 2Concept of surface dynamics |
$ {v_{n1}} = - {c_r}{v_{n0}} $ | (8) |
其中:vn0和vn1分别表示碰撞前和碰撞后的法向速度。切向运动满足如下的关系:
$ {v_{t1}} = {v_{t0}} - {\mu _r}\left( {{v_{n0}} + {v_{n1}}} \right) = {v_{t0}} - {\mu _r}\left( {1 + {c_r}} \right){v_{n0}} $ | (9) |
其中:vt0和vt1分别表示碰撞前和碰撞后的切向速度;μr表示当地的摩擦系数。
1.4 简化模型忽略小天体自旋和引力不规则变化的影响,可以将探测器的弹跳运动近似为抛物运动。通过对抛物运动的分析,可以近似得到通过一次弹跳将探测器从起始点转移到目标点所需要的速度vi。
定义矢量R为起始点指向单次弹跳目标点的矢量,g为起始点处的引力加速度矢量。构建R-g平面如图 3所示。R方向的单位矢量为
![]() |
图 3抛物运动模型Fig. 3Parabolic motion |
$ \hat {{R}} = \frac{{{R}}}{R} $ | (10) |
在R-g平面内,定义矢量l与R垂直,并指向小天体外侧。gl和gR分别为g沿l和R方向的分量,假设引力加速度g不变,给定期望弹跳时间∆t,那么vi沿l方向的速度分量为
$ {v_l} = \frac{{{g_l}\Delta t}}{2} $ | (11) |
vi沿R方向的速度分量为
$ {v_R} = \frac{{2R - {g_R}\Delta {t^2}}}{{2\Delta t}} $ | (12) |
将l方向和R方向的速度合成,可以得到一次弹跳将探测器从起始点转移到目标点所需要的速度vi为
$ {{{V}}_i} = {v_l}\hat {{l}} + {v_R}\hat {{R}} $ | (13) |
其中:$\hat {{l}}$为沿l方向的单位矢量。
由式(13) 可以近似得到通过一次弹跳将探测器从起始点转移到目标点所需要的速度vi,但是由于忽略了小天体自旋和引力不规则变化的影响,探测器的实际附着位置将会偏离预定的目标点。因此需要在探测器飞行过程中施加轨迹修正指令,对探测器的运动轨迹进行修正。
2 控制方法本文所提出的控制方法,可以概括如下:1) 在起始点与目标点之间设置多个中间路径点;2) 在起始点处,通过式(13) 求出探测器到第一个中间路径点所需的初始速度vi,并由脉冲发动机实现这一速度;3) 当探测器飞行一段时间后,施加轨迹修正指令,使探测器到达第一个中间路径点;4) 以当前路径点作为起始点,重复前三步,直到探测器抵达目标点;5) 施加反向速度脉冲,使探测器最终停留在目标点。表面运动控制方法示意图如图 4所示。本节将重点介绍轨迹修正算法的实现。
![]() |
图 4控制方法示意图Fig. 4Scenario of hopping exploration on the small body |
选取探测器在小天体固连坐标系下的位置和速度作为状态变量${{X}} = {\left[ {{{{r}}^{\rm{T}}},{{{V}}^{\rm{T}}}} \right]^{\rm{T}}} = {\left[ {x,y,z,{v_x},{v_y},{v_z}} \right]^{\rm{T}}}$。将轨迹修正指令施加处探测器所处的状态${{{X}}_0} = {\left[ {{{r}}_0^{\rm{T}},\;{{V}}_0^{\rm{T}}} \right]^{\rm{T}}}$ 选取为初始状态。此时,探测器满足轨道运动方程(7) , 控制加速度ax=ay=az=0,对方程(7) 进行泰勒展开可以得到
$\left[ {\begin{array}{*{20}{c}} {\dot x}\\ {\dot y}\\ {\dot z}\\ {{{\dot v}_x}}\\ {{{\dot v}_y}}\\ {{{\dot v}_z}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{v_x}}\\ {{v_y}}\\ {{v_z}}\\ {\;2\omega {v_y} + {\omega ^2}x + {V_{x0}} + {V_{xx0}}\left( {x - {x_0}} \right) + {V_{xy0}}\left( {y - {y_0}} \right) + {V_{xz0}}\left( {z - {z_0}} \right) + \;...\;}\\ {\; - 2\omega {v_x} + {\omega ^2}y + {V_{y0}} + {V_{xy0}}\left( {x - {x_0}} \right) + {V_{yy0}}\left( {y - {y_0}} \right) + {V_{yz0}}\left( {z - {z_0}} \right) + \;...\;}\\ {{V_{z0}} + {V_{xz0}}\left( {x - {x_0}} \right) + {V_{yz0}}\left( {y - {y_0}} \right) + {V_{zz0}}\left( {z - {z_0}} \right) + \;...\;} \end{array}} \right]$ | (14) |
下标0表示相关参数在初始状态处的取值。仅保留0阶项和1阶项,可以得到轨道运动的线性化方程为
$ \dot {{X}} = \left( {{{A}} + {{δ}} } \right){{X}} + {{u}} $ | (15) |
其中$\left( {{{A}},\;{{δ}} ,{{u}}} \right) \in \left( {{{\rm IR}^{{\text{6}} \times {\text{6}}}},\;{{\rm IR}^{{\text{6}} \times {\text{6}}}},{{\rm IR}^{{\text{6}} \times {\text{1}}}}} \right)$定义如下
$\begin{array}{*{20}{l}} {A = \left[ {\begin{array}{*{20}{c}} 0&0&0&1&0&0\\ 0&0&0&0&1&0\\ 0&0&0&0&0&1\\ {{\omega ^2}}&0&0&0&{2\omega }&0\\ 0&{{\omega ^2}}&0&{ - 2\omega }&0&0\\ 0&0&0&0&0&0 \end{array}} \right],}\\ {\delta = \left[ {\begin{array}{*{20}{c}} {{0_{3 \times 3}}}&{{0_{3 \times 3}}}\\ {{{\left. {\nabla {V^2}} \right|}_0}}&{{0_{3 \times 3}}} \end{array}} \right],\;u = {{\left[ {{0_{1 \times 3}}\;,\;{{\left( {{{\left. {\nabla V} \right|}_0}, - {{\left. {\nabla {V^2}} \right|}_0},{r_0}} \right)}^{\rm{T}}}} \right]}^{\rm{T}}}.} \end{array}$ |
${\left. {\nabla {{V}}} \right|_0}$和${\left. {\nabla {{{V}}^2}} \right|_0}$分别为引力势函数的一阶导数矩阵和二阶导数矩阵,定义如下
${\left. {\nabla V} \right|_0} = \left[ {\begin{array}{*{20}{c}} {{V_{x0}}}&0&0\\ 0&{{V_{y0}}}&0\\ 0&0&{{V_{z0}}} \end{array}} \right]\;\;{\left. {\nabla {V^2}} \right|_0} = \left[ {\begin{array}{*{20}{c}} {{V_{xx0}}}&{{V_{xy0}}}&{{V_{xz}}_0}\\ {{V_{xy0}}}&{{V_{yy0}}}&{{V_{yz0}}}\\ {{V_{xz0}}}&{{V_{yz0}}}&{{V_{zz0}}} \end{array}} \right]$ |
方程(15) 为线性定常系统,该系统的解为
$ {{X}}\left( t \right) = {{\rm e}^{\left( {{{A}} + {{δ}} } \right)t}}{{{X}}_0}{\rm{ + }} \rm{\int}_0^t {{e^{\left( {{{A}} + {{δ}}} \right)\left( {t - \tau } \right)}}{{u}}{\rm{d}}\tau } $ | (16) |
u为一个常值矩阵,${\left. {\frac{{{\partial ^2}{{V}}}}{{\partial {{{r}}^2}}}} \right|_0} + {\rm diag}\left( {{\omega ^2},{\omega ^2},0} \right)$是可逆的,所以
$ \rm{\int}_0^t {{e^{\left( {{{A}} + {{δ}} } \right)\left( {t - \tau } \right)}}{{u}}{\rm{d}}\tau } = {\left( {{{A}} + {{δ}} } \right)^{ - 1}}\left[ {{e^{\left( {{{A}} + {{δ}} } \right)t}} - {{{I}}_{6 \times 6}}} \right]{{u}} $ | (17) |
将方程(17) 代入方程(16) 可以得到
$ {{X}}\left( t \right) = {{\rm e}^{\left( {{{A}} + {{δ}} } \right)t}}{{{X}}_0}{\rm{ + }}{\left( {{{A}} + {{δ}} } \right)^{ - 1}}\left[ {{{\rm e}^{\left( {{{A}} + {{δ}} } \right)t}} - {{{I}}_{6 \times 6}}} \right]{{u}} $ | (18) |
已知探测器的当前位置r0和目标位置rt,并给定转移时间tf,通过对方程(18) 进行性逆运算,可以求出轨迹修正后的速度为
$ {{{v}}_e} = {{Φ}} {\left( {1\!\!:\!\!3,4\!\!:\!\!6} \right)^{ - 1}} \cdot \left( {{{{r}}_t} - {{Φ}} \left( {1\!\!:\!\!3,1\!\!:\!\!3} \right){{{r}}_0} - {{G}}\left( {1\!\!:\!\!3} \right)} \right) $ | (19) |
矩阵G和Φ的取值如下:
$ \left\{ \begin{array}{l} {{Φ}} = {{\rm e}^{\left( {{{A}} + {{δ}} } \right){t_f}}}\\ {{G}} = {\left( {{{A}} + {{δ}} } \right)^{ - 1}}\left[ {{{\rm e}^{\left( {{{A}} + {{δ}} } \right){t_f}}} - {{{I}}_{6 \times 6}}} \right]{{u}} \end{array} \right. $ | (20) |
由式(19) 可以求出轨迹修正点的速度增量为
$\Delta v = {v_e} - {v_0}$ | (21) |
以上分析给出了单次弹跳的控制过程,在探测器在小天体表面进行弹跳运动时,弹跳次数的设计应该考虑转移距离的影响。图 4给出了初始速度vi大小固定的情况下,即$v_l^2 + v_R^2 = C$(其中C为一常数) ,小行星Geographos表面某点的初始法向速度与转移距离的关系图。由图 5可以发现,当转移距离较远的时候,法向速度可能会超越当地的法向逃逸速度,这会使探测器在飞行一段时间后脱离引力束缚,离开小天体表面。另外,当转移距离较远时,误差积累也较大,这使得后续的轨迹修正过程燃耗增加。基于以上考虑,对于距离远的目标点,应该使探测器通过多次弹跳到达。
![]() |
图 5转移距离与法向初始速度Fig. 5Relationship between transfer distance and normal initial velocity |
假设探测器停留在小天体表面,位置矢量为${{{r}}_0} = {\left[ {554.7,5{\rm{7}}{\rm{.6}},9{\rm{60}}{\rm{.0}}} \right]^{\rm{T}}}$,目标位置矢量为${{{r}}_t} = \left[ {577.3,} \right.$ ${\left. {159.2,933.7} \right]^{\rm{T}}}$。由式(13) 可以求出初始速度vi,仿真运动轨迹如图 6所示。图中实线为真实运动轨迹,虚线为忽略小天体自旋和引力不规则变化的情况下得到的运动轨迹。可以看出,初始速度vi可以使探测器以弹跳的形式向目标点运动,但由于做了简化处理,实际着陆点会远远偏离目标点。所以,在探测器运动的过程,有必要进行轨迹修正。
![]() |
图 6弹跳运动轨迹Fig. 6Bouncing trajectory in 3D |
本组仿真采用与前一次仿真相同的初末位置,在这组仿真中,探测器受初始速度vi的驱动,以弹跳的形式向目标点运动,当探测器飞行一段时间并到达一定高度后,由式(21) 求出轨迹修正指令∆v,并最终完成一次弹跳运动。仿真结果如图 7所示。
![]() |
图 7固连坐标系下运动轨迹Fig. 7Trajectory histories in body-fixed coordinate |
由图 7可以看出,运动轨迹可以分为两段,前一段轨迹为轨迹修正指令作用前的结果,与前一次的仿真轨迹重合,后一段为轨迹修正指令作用后的运动轨迹。最终,探测器能够准确地抵达目标点。
3.3 远距离弹跳转移前两组仿真验证了文中所提的表面运动控制方法在单次弹跳运动里的控制效果。本组仿真用来验证该控制方法在远距离转移的控制效果。初始点的位置矢量为${\left[ {604.8,\; - 106.2,\;1005} \right]^{\rm{T}}}$,目标点的位置矢量为${\left[ {680.2,\;370.4,\;926.8} \right]^{\rm{T}}}$。分别设置3个中间路径点和4个中间路径点,使探测器以4次弹跳运动和5次弹跳运动的方式到达目标点,在目标点处,施加与当前速度等大小的反向速度脉冲,使探测器最终停留在目标点。运动轨迹如图 8所示。由仿真结果可以发现,探测器准确抵达目标点,并且在整个运动过程中没有发生逃逸。
![]() |
图 8固连坐标系下运动轨迹Fig. 8Trajectory histories in body-fixed coordinate |
本文针对探测器在小天体表面的弹跳式运动特点,提出了一种表面运动控制方法。该控制方法可以应用于探测器在附着到小天体表面后的弹跳转移,也适用于弱引力环境下探测器的弹跳式探测。
本文首先对探测器在小天体附近运动的引力场模型和动力学模型进行了介绍,并基于简化的运动模型,计算出探测器进行弹跳式转移所需的初始速度;通过对轨道方程的线性化,基于线性时变系统的求解方程,采用脉冲控制的方法对每次弹跳的运动轨迹进行修正,增加了控制精度;同时分析了转移距离与弹跳次数的关系;最后,以小行星1620 Geographos作为目标小天体进行仿真分析,验证了控制方法的有效性。
[1] | 崔平远, 乔栋. 小天体附近轨道动力学与控制研究现状与展望[J].力学进展, 2013, 43(5): 526-539. |
[2] | Guibout V, Scheeres D J. Stability of surface motion on a rotating ellipsoid[J]. Celestial Mechanics and Dynamical Astronomy, 2003, 87(3): 263-290.doi: 10.1023/B:CELE.0000005720.09027.ee |
[3] | Liu X, Baoyin H, Ma X. Dynamics of surface motion on a rotating massive homogeneous body[J]. Science China Physics, Mechanics and Astronomy, 2013, 56(4): 818-829.doi: 10.1007/s11433-013-5044-2 |
[4] | Yu Y, Baoyin H X. Modeling of migrating grains on asteroid's surface[J].Astrophysics and Space Science, 2015, 335(1): 43-56. |
[5] | Tardivel S, Scheeres D J, Michel P, et al. Contact motion on surface of asteroid[J]. Journal of Spacecraft and Rockets, 2014, 51(6): 1857-1871.doi: 10.2514/1.A32939 |
[6] | Sawai S, Kawaguchi J, Scheeres D, et al. Development of a target marker for landing on asteroids,[J]. Journal of Spacecraft and Rockets, 2001, 38(4): 601-608.doi: 10.2514/2.3723 |
[7] | Kömle N I, Ball A J, Kargl G, et al. Using the anchoring device of a comet lander to determine surface mechanical properties[J].Planetary and Space Science, 1997, 48(12): 1515-1538. |
[8] | Ulamec S, Biele J. Surface elements and landing strategies for small bodies missions-philae and beyond[J]. Advances in Space Research, 2009, 44(7): 847-858.doi: 10.1016/j.asr.2009.06.009 |
[9] | Kubota T,Otsuki M,Hashimoto T,et al. Touchdown dynamics for sampling in hayabusa mission[C]//In Proceedings of the AIAA/AAS Astrodynamics Specialist Conference and Exhibit. Keystone:AIAA,2006:2006-6539. |
[10] | Berry K,Sutter B,May A,et al. Osiris-rex touch-and-go (TAG) mission design and analysis[C]//36th Annual AAS Guidance and Control Conference. Breckenridge:AAS,2013:13-95. |
[11] | Ulamec S, Kucherenko V, Biele J, et al. Hopper concepts for small body landers[J]. Advances in Space Research, 2011, 47(3): 428-439.doi: 10.1016/j.asr.2010.09.013 |
[12] | Mège D, Gurgurewicz J, Grygorczuk J, et al. The Highland terrain hopper (hopter):concept and use cases of a new locomotion system for the exploration of low gravity solar system bodies[J].Acta Astronautica, 2016(121): 200-220. |
[13] | Sagdeev R Z, Zakharov A V. Brief history of the Phobos mission[J]. Nature, 1989, 341(6243): 581-585.doi: 10.1038/341581a0 |
[14] | Bellerose J,Scheeres D J. Dynamics and control for surface exploration of small bodies[C]//Proceedings of AIAA/AAS 2008 Astrodynamics Specialist Conference. Honolulu:AIAA,2008:2008-6251.http://cn.bing.com/academic/profile?id=496685cbd54a20a876a35b3987189fb3&encoded=0&v=paper_preview&mkt=zh-cn |
[15] | Yoshimitsu T, Kubota T, Nakatani I, et al. Micro-hopping robot for asteroid exploration[J].Acta Astronautica, 2003, 52(2): 441-446. |