龔宇鵬，謝 天，林曉輝

(哈爾濱工業(yè)大學(xué) 衛(wèi)星技術(shù)研究所，黑龍江 哈爾濱 150080)

0 引言

隨著航天技術(shù)的飛速發(fā)展，各航天大國開展了一系列深空探測(cè)活動(dòng)，并將目光投向了太陽系中的小行星、彗星等小型天體。這些天體數(shù)量眾多，對(duì)人類研究太陽系和恒星的成因及演化歷史，探索生命起源，獲取太空資源具有重要作用[1]。在深空探測(cè)中，小型天體本身的特殊形狀導(dǎo)致了其引力場(chǎng)的特殊性，不能將它們單純視作中心引力場(chǎng)，因此在航天器對(duì)其接近和著陸的過程中，可以將小型天體視作多個(gè)引力中心的組合，從而精確描述其引力特性[2]。

目前，只有美國“近地小行星交會(huì)”(NEAR)探測(cè)器、日本“隼鳥一號(hào)”探測(cè)器、歐洲太空局“羅塞塔號(hào)”彗星探測(cè)器與小行星表面有過直接接觸[3-5]。2000年2月，NEAR探測(cè)器在設(shè)計(jì)任務(wù)之外完成了在“愛神”小行星上的軟著陸，獲得了寶貴的探測(cè)資料和成果。“隼鳥一號(hào)”探測(cè)器于2005年11月在“絲川”小行星表面軟著陸并采集了土壤樣品。2014年11月，由“羅塞塔號(hào)”彗星探測(cè)器釋放的“菲萊”探測(cè)器成功登陸“67P/楚留莫夫-格拉希門克”彗星，但由于探測(cè)器著陸位置不良，導(dǎo)致其在很長(zhǎng)時(shí)間內(nèi)處于休眠狀態(tài)。另外，美國“黎明號(hào)”探測(cè)器、日本“隼鳥二號(hào)”、中國“嫦娥二號(hào)”探測(cè)器也已成功開展了對(duì)小行星的探測(cè)任務(wù)。其中“嫦娥二號(hào)”探測(cè)器于2012年12月成功近距離飛躍“戰(zhàn)神-圖塔蒂斯”小行星，實(shí)現(xiàn)了中國首次對(duì)小行星的近距離觀測(cè)[6]，走出了中國深空探測(cè)極其重要的一步。

研究小行星探測(cè)器運(yùn)動(dòng)及控制方式能為小行星探測(cè)方案設(shè)計(jì)提供依據(jù)，對(duì)深空探測(cè)具有重要意義。目前國內(nèi)外成功使用的行星表面探測(cè)器均為輪式車型探測(cè)器。雖然輪式車型探測(cè)器在月球、火星等大質(zhì)量天體上的使用獲得了成功，但在小行星的微重力場(chǎng)下，由于探測(cè)器車輪與地面之間的接觸力非常小，使車輪打滑而移動(dòng)緩慢，且在小行星粗糙的表面，一個(gè)小擾動(dòng)就會(huì)使探測(cè)器脫離地面，因此，輪式車型探測(cè)器并不適用于小行星表面探測(cè)，需要尋求其他可行的方案。

美國、日本、歐洲各國和地區(qū)提出了一種新的行走方案——跳躍行走，并對(duì)其開展了相關(guān)研究[7]。跳躍行走能很好地適應(yīng)小行星表面的微重力環(huán)境。探測(cè)器從小行星表面起跳后，在空中進(jìn)行彈道運(yùn)動(dòng)并落在另一個(gè)位置，從而實(shí)現(xiàn)了行走的目的。探測(cè)器在跳躍行走過程中能輕易越過幾倍甚至幾十倍于自身尺寸的障礙物，能適應(yīng)小行星復(fù)雜和不可預(yù)測(cè)的地形地貌，且探測(cè)器的運(yùn)動(dòng)范圍大，能耗低。因此，跳躍行走是小行星表面巡游探測(cè)比較可行的方案。

國內(nèi)外針對(duì)探測(cè)器在小行星表面跳躍行走的實(shí)現(xiàn)開展了大量研究。如HANAZAWA等[6-7]設(shè)計(jì)了一種輻條輪機(jī)構(gòu)，可在松軟的行星表面進(jìn)行移動(dòng)；柏龍等[9]設(shè)計(jì)了用于小行星探測(cè)的跳躍機(jī)構(gòu)，通過折疊機(jī)構(gòu)將彈性勢(shì)能轉(zhuǎn)變?yōu)樘綔y(cè)器彈跳的動(dòng)能，具有極強(qiáng)的越障和行走能力；GAJAMOHAN等[10-11]提出了利用反作用飛輪進(jìn)行跳躍的立方體模型，可通過飛輪制動(dòng)，短時(shí)間內(nèi)獲得可使探測(cè)器起跳的較大力矩，但其對(duì)起跳后探測(cè)器的姿態(tài)、跳躍過程的控制問題沒有深入研究；PAVONE等[12-13]設(shè)計(jì)了一款仿刺猬機(jī)器人，其通過剛性桿與地面接觸，并通過反作用飛輪實(shí)現(xiàn)對(duì)起跳和移動(dòng)的控制，但未對(duì)機(jī)器人的移動(dòng)策略進(jìn)行進(jìn)一步研究；SHEN等[14-15]針對(duì)探測(cè)器在小行星表面的跳躍運(yùn)動(dòng)進(jìn)行了研究，但未考慮探測(cè)器在運(yùn)動(dòng)過程中的控制問題。

采用跳躍機(jī)構(gòu)會(huì)增加小行星探測(cè)器質(zhì)量及成本，容易降低探測(cè)器可靠性，在不增加額外跳躍執(zhí)行機(jī)構(gòu)的情況下，基于反作用飛輪的小行星跳躍探測(cè)器具有結(jié)構(gòu)簡(jiǎn)單、可靠性強(qiáng)的優(yōu)勢(shì)。為此，本文針對(duì)輪控探測(cè)器建模及跳躍過程的控制進(jìn)行研究。

1 跳躍移動(dòng)原理

目前，對(duì)小行星跳躍探測(cè)器跳躍機(jī)理的研究主要分為兩類：一是通過跳躍機(jī)構(gòu)將擠壓力作用于行星表面，使探測(cè)器獲得水平速度的彈跳跳躍方式；二是通過反作用飛輪的加速改變探測(cè)器與地面的接觸力，使探測(cè)器起跳。

1.1 彈跳跳躍原理

彈跳跳躍的靈感主要來源于將跳躍移動(dòng)作為主要移動(dòng)方式的跳蚤、青蛙、袋鼠等動(dòng)物。它們一般具有強(qiáng)壯的下肢，在跳躍時(shí)依靠下肢折疊產(chǎn)生彈力，并利用上肢調(diào)整跳躍姿態(tài)，以實(shí)現(xiàn)跳躍移動(dòng)[8-9]。

彈跳探測(cè)器的跳躍移動(dòng)需要通過跳躍機(jī)構(gòu)來執(zhí)行，其原理為：跳躍機(jī)構(gòu)通過自身折疊儲(chǔ)存彈性勢(shì)能，積累后迅速釋放，使地面支持力增加，同時(shí)產(chǎn)生前向的地面摩擦力，使探測(cè)器躍起并向前移動(dòng)。

圖1 彈跳跳躍原理示意Fig.1 Principle of bounce hopping motion

根據(jù)不同的工程應(yīng)用，彈跳探測(cè)器所采用的跳躍機(jī)構(gòu)或運(yùn)動(dòng)方式各有不同，但總結(jié)其運(yùn)動(dòng)特性，可以得知彈跳跳躍方式有以下幾個(gè)特點(diǎn)：

1)具有極強(qiáng)的越障能力，可以越過相較于自身尺寸數(shù)十倍的障礙物，落腳面小，能夠適應(yīng)復(fù)雜地形與未知環(huán)境。

2)能更好地適應(yīng)小行星的微重力環(huán)境，彈跳行走過程能耗較小，并能夠解決地面摩擦力不足的問題。

3)彈跳探測(cè)器一般具有較輕的質(zhì)量和較大的底面，可有效減小壓強(qiáng)，在松軟的行星表面跳躍效果較好。

盡管彈跳跳躍能很好解決探測(cè)器在微重力環(huán)境中的移動(dòng)問題，但由于引入了跳躍機(jī)構(gòu)，增加了探測(cè)器的結(jié)構(gòu)質(zhì)量，降低了可靠性，且難以在跳躍過程中對(duì)探測(cè)器姿態(tài)、運(yùn)動(dòng)方向進(jìn)行控制，無法保證其穩(wěn)定性，容易發(fā)生傾翻，造成儀器及設(shè)備的損壞。

1.2 基于反作用飛輪的跳躍原理

基于反作用飛輪跳躍的起跳過程如圖2所示，主要包括：

1)反作用飛輪從初始狀態(tài)開始加速，隨著飛輪力矩的逐漸增大，為平衡飛輪力矩，地面等效支持力的作用點(diǎn)將逐漸向探測(cè)器起跳支點(diǎn)D處移動(dòng)。當(dāng)?shù)孛嬷С至ν耆饔糜谄鹛c(diǎn)且重力力矩?zé)o法與飛輪力矩相平衡時(shí)，與倒立擺起擺時(shí)的原理相似，探測(cè)器將繞起跳支點(diǎn)發(fā)生轉(zhuǎn)動(dòng)。

2)當(dāng)探測(cè)器發(fā)生傾斜后，重力相對(duì)于支點(diǎn)處的力矩將隨傾角的增加而減小，而飛輪仍處于加速狀態(tài)，因此探測(cè)器將繞其支點(diǎn)加速旋轉(zhuǎn)。此時(shí)，探測(cè)器質(zhì)心的水平速度與豎直速度都將逐漸增加，直到探測(cè)器重力不足以作為以支點(diǎn)D為圓心的加速圓周運(yùn)動(dòng)的向心力時(shí)，探測(cè)器躍起。

3)在探測(cè)器躍起后，由于其在旋轉(zhuǎn)過程中已獲得了一定的水平及豎直速度，探測(cè)器將脫離小行星表面進(jìn)行拋體運(yùn)動(dòng)。在探測(cè)器起跳后，飛輪可以逐漸減速，提供反作用力矩，使探測(cè)器姿態(tài)恢復(fù)平穩(wěn)。

圖2 探測(cè)器起跳過程的二維模型Fig.2 2D model of probe takeoff process

探測(cè)器的成功起跳是在小行星表面進(jìn)行跳躍探測(cè)的基礎(chǔ)。為了研究起跳過程，需要探討能夠使探測(cè)器產(chǎn)生翻滾或跳躍的飛輪力矩最小值，并且需要知道施加力矩與探測(cè)器跳躍瞬間運(yùn)動(dòng)參數(shù)的關(guān)系。對(duì)探測(cè)器起跳的二維動(dòng)力學(xué)模型進(jìn)行分析。根據(jù)牛頓第二定律和動(dòng)量矩原理，得到飛輪起跳過程中相對(duì)于起跳支點(diǎn)D的動(dòng)力學(xué)方程，即

(1)

假設(shè)探測(cè)器兩頂點(diǎn)之間的夾角為2β，則由式(1)可得探測(cè)器所需最小起跳控制力矩的公式，即

Tmin>mgLsinβ

(2)

式(2)也是跳躍探測(cè)器能夠起跳的條件。假設(shè)探測(cè)器模型為質(zhì)量分布對(duì)稱的均勻立方體，則L為立方體半對(duì)角線的長(zhǎng)度。由式(2)可知：L越大，則飛輪起跳所需的控制力矩就越大；同時(shí)，探測(cè)器起跳所需的控制力矩也與探測(cè)器2個(gè)頂點(diǎn)之間的夾角有關(guān)，當(dāng)L不變時(shí)，該夾角越大，探測(cè)器起跳所需要的飛輪輸出力矩也就越大。

2 探測(cè)器動(dòng)力學(xué)建模

2.1 探測(cè)器姿態(tài)動(dòng)力學(xué)模型

將小行星跳躍探測(cè)器假設(shè)為1個(gè)剛性立方體(見圖3)，探測(cè)器的質(zhì)量均勻分布，在質(zhì)心處有三軸正交安裝的反作用飛輪，通過飛輪加速及制動(dòng)改變探測(cè)器的姿態(tài)，實(shí)現(xiàn)跳躍移動(dòng)。

圖3 小行星跳躍探測(cè)器三維模型Fig.3 3D model of asteroid hopping probe

為明確反作用飛輪對(duì)于探測(cè)器姿態(tài)的控制作用，需要推導(dǎo)飛輪與探測(cè)器的姿態(tài)動(dòng)力學(xué)方程。根據(jù)歐拉方程可得

(3)

式中：下標(biāo)“I”代表探測(cè)器在慣性坐標(biāo)系下的時(shí)間導(dǎo)數(shù)；下標(biāo)“B”代表探測(cè)器在本體坐標(biāo)系下的時(shí)間導(dǎo)數(shù)；H為探測(cè)器的角動(dòng)量。

由歐拉方程可以推導(dǎo)得到在系統(tǒng)質(zhì)心坐標(biāo)系Os-xsyszs下帶慣性輪的多剛體姿態(tài)動(dòng)力學(xué)方程，即

(4)

式中：I為系統(tǒng)關(guān)于其質(zhì)心的慣性張量在系統(tǒng)質(zhì)心坐標(biāo)系下的矩陣；ω為探測(cè)器的角速度矢量在系統(tǒng)質(zhì)心坐標(biāo)系下的矩陣；U為各動(dòng)量輪組成的安裝矩陣；Jw為各動(dòng)量輪轉(zhuǎn)動(dòng)慣量組成的對(duì)角陣；Ωw為各動(dòng)量輪的轉(zhuǎn)速組成的列陣；Tb為外力矩在系統(tǒng)質(zhì)心坐標(biāo)系下的分量列陣，Tb=[TxTyTz]T=rOsX×F，rOsX為小行星坐標(biāo)系下探測(cè)器質(zhì)心到與地面碰撞點(diǎn)的向量,F為探測(cè)器所受的合外力。

探測(cè)器內(nèi)部的3個(gè)飛輪是正交安裝的，安裝矩陣U為單位陣，則

(5)

(6)

式中：Mk為飛輪k的控制輸入力矩，k=B，C，D。通過控制Mk可執(zhí)行不同探測(cè)器控制策略。

2.2 反作用飛輪模型

反作用飛輪的工作原理是動(dòng)量矩定理，即星體的總動(dòng)量矩矢量對(duì)時(shí)間的導(dǎo)數(shù)等于作用于星體上外力矩矢量的總和。反作用飛輪通常作為執(zhí)行機(jī)構(gòu)對(duì)航天器進(jìn)行姿態(tài)控制，本文中反作用飛輪是探測(cè)器姿態(tài)控制及跳躍運(yùn)動(dòng)的唯一執(zhí)行機(jī)構(gòu)，即采用飛輪加速既能使探測(cè)器具有運(yùn)動(dòng)速度，還可利用其進(jìn)行姿態(tài)調(diào)整和控制。

反作用飛輪系統(tǒng)通?？珊?jiǎn)化為一階系統(tǒng)進(jìn)行控制，即

(7)

式中：Tw=RJ/(KeKm)為飛輪電動(dòng)機(jī)的機(jī)電時(shí)間常數(shù)；R為電動(dòng)機(jī)阻抗；J為反作用飛輪轉(zhuǎn)動(dòng)慣量；Ke為電勢(shì)常數(shù)；Km為力矩常數(shù)；Kw=1/Ke為飛輪電動(dòng)機(jī)的增益系數(shù)。

輪控系統(tǒng)的實(shí)際應(yīng)用一般更關(guān)心反作用飛輪產(chǎn)生的實(shí)際力矩M對(duì)期望控制信號(hào)Tc的響應(yīng)，因此希望將控制器環(huán)境生成的期望控制力矩，作為飛輪的輸入(控制)信號(hào)進(jìn)行輸入。在飛輪電機(jī)的模型中加入電樞電流反饋，調(diào)節(jié)電機(jī)電樞電流，產(chǎn)生所需要的電磁驅(qū)動(dòng)力矩，它們之間的關(guān)系為

(8)

(9)

式中：Tc(t)為飛輪的期望控制力矩函數(shù)；i(t)為電樞電流；T′w=R/KP；KP為比例系數(shù)。式(9)為力矩模式飛輪的基本數(shù)學(xué)模型，可看成1個(gè)慣性環(huán)節(jié)。對(duì)于性能較好的慣性輪，慣性環(huán)節(jié)的時(shí)間常數(shù)可達(dá)到毫秒級(jí)。本文選取T′w為0.01 s。

本文選用的反作用飛輪性能指標(biāo)見表1。

表1 反作用飛輪的性能指標(biāo)

2.3 探測(cè)器與地面的接觸力模型

由于在飛輪加速使探測(cè)器起跳的過程中需要計(jì)算碰撞時(shí)間和碰撞力等信息，而傳統(tǒng)的瞬間碰撞模型無法揭示這些信息，因此采用等效彈簧阻尼模型，將探測(cè)器與地面之間的接觸等效為1個(gè)彈簧阻尼模型。在等效彈簧阻尼模型中，法向受力由基于Hertz定律和非線性阻尼函數(shù)的非線性彈簧-阻尼模型計(jì)算得到[16-17]。切向力由簡(jiǎn)化的Karnopp摩擦力模型計(jì)算得到[18]。

探測(cè)器與地面的法向碰撞力用彈簧接觸力來描述，由Hertz定律確定力的大小，用阻尼系數(shù)表征接觸過程的能量損失?；贖ertz定律和非線性阻尼函數(shù)的非線性彈簧-阻尼模型表達(dá)式為

(10)

(11)

(12)

為描述探測(cè)器的切向碰撞力，需要考慮到探測(cè)器在跳躍過程中可能發(fā)生的粘滯與滑動(dòng)兩種運(yùn)動(dòng)狀態(tài)，因此引入簡(jiǎn)化的Karnopp模型對(duì)探測(cè)器切向摩擦力進(jìn)行描述。Karnopp首先定義了一個(gè)速度轉(zhuǎn)換容限VD，當(dāng)切向速度小于VD時(shí)，將探測(cè)器狀態(tài)視為粘滯，否則為滑動(dòng)狀態(tài)。因此，針對(duì)簡(jiǎn)化后的切向碰撞力為

(13)

式中：Fs=μ0Fn為靜摩擦力。

3 探測(cè)器控制過程及控制律設(shè)計(jì)

跳躍探測(cè)器在小行星表面的跳躍運(yùn)動(dòng)主要依靠其內(nèi)部的飛輪來驅(qū)動(dòng)，飛輪的控制力矩是整個(gè)跳躍系統(tǒng)的控制輸入。在開環(huán)條件下，探測(cè)器能在一定程度上適應(yīng)地形變化，但探測(cè)器的路徑會(huì)與預(yù)定的路徑有較大偏差，無法準(zhǔn)確到達(dá)預(yù)定地點(diǎn)進(jìn)行探測(cè)。通過合理的控制算法控制飛輪的輸入力矩，就能夠控制飛輪的轉(zhuǎn)速和角加速度，實(shí)現(xiàn)對(duì)探測(cè)器跳躍的速度、方向等的控制，實(shí)現(xiàn)探測(cè)器對(duì)行星表面多個(gè)預(yù)定地點(diǎn)的準(zhǔn)確到達(dá)。在整個(gè)探測(cè)過程中，需要控制探測(cè)器的速度不能過大，防止由于速度超過小行星逃逸速度而脫離小行星的引力場(chǎng)，也要注意飛輪的飽和及卸載問題。另外，跳躍探測(cè)器從地面躍起后，由于其具有初始旋轉(zhuǎn)速度，會(huì)在空中進(jìn)行旋轉(zhuǎn)和翻滾。如果不對(duì)其進(jìn)行控制，則探測(cè)器落地時(shí)隨機(jī)的姿態(tài)會(huì)使下次跳躍變得十分復(fù)雜。因此，在實(shí)際應(yīng)用中期望每次都以固定的姿態(tài)及角速度落回地面，這就需要在每次探測(cè)器跳躍的過程中通過飛輪的控制來調(diào)整其姿態(tài)，以控制探測(cè)器落地姿態(tài)。

探測(cè)器起跳或與地面發(fā)生碰撞后，均會(huì)產(chǎn)生一定的角速度。若探測(cè)器內(nèi)部角速度具有耦合時(shí)，探測(cè)器還可能具有角加速度，這會(huì)使探測(cè)器的姿態(tài)運(yùn)動(dòng)變得復(fù)雜，不利于探測(cè)器的控制，因此主要在探測(cè)器滯空時(shí)對(duì)探測(cè)器姿態(tài)進(jìn)行調(diào)整。

為了對(duì)探測(cè)器進(jìn)行姿態(tài)穩(wěn)定控制，需要引入探測(cè)器本體坐標(biāo)系與目標(biāo)姿態(tài)坐標(biāo)系間的誤差四元數(shù)Qe=[qe0qe]T∈R4×1，假設(shè)小行星本體系下的探測(cè)器當(dāng)前姿態(tài)四元數(shù)為Q，目標(biāo)姿態(tài)四元數(shù)Qf=[qf0qf]T，則誤差四元數(shù)

Qe=Q-1Qf

(14)

式中：Qe滿足的約束條件為

qe02+qe12+qe22+qe32=1

(15)

(16)

目標(biāo)姿態(tài)四元數(shù)可根據(jù)不同的姿態(tài)任務(wù)確定，當(dāng)前姿態(tài)四元數(shù)可由姿態(tài)測(cè)量系統(tǒng)測(cè)量得到。假設(shè)探測(cè)器當(dāng)前的姿態(tài)角速度為ω，目標(biāo)角速度為ωf，均定義在航天器本體坐標(biāo)系中，則定義誤差角速度

ωe=ω-ωf

(17)

將式(17)代入式(4)，可得到誤差動(dòng)力學(xué)與運(yùn)動(dòng)學(xué)方程，即

(ωe+ωf)×[I(ωe+ωf)+JwΩw]+

(18)

(19)

本節(jié)中探測(cè)器的姿態(tài)機(jī)動(dòng)控制問題就是在給定目標(biāo)姿態(tài)四元數(shù)下，設(shè)計(jì)控制律使得當(dāng)前姿態(tài)四元數(shù)在系統(tǒng)存在參數(shù)攝動(dòng)及外界干擾的情況下，能在有限時(shí)間內(nèi)Qe→[±1 0 0 0]T，ωe→0。

為了達(dá)到更好的控制效果，考慮到飛輪慣性環(huán)節(jié)及克服摩擦力矩的干擾，本文在控制環(huán)節(jié)采取的控制率為

(20)

式中：KL為控制器的回路增益矩陣；k1，k2，k3分別為比例、微分、積分控制項(xiàng)。經(jīng)調(diào)試，k1=0.02，k2=0.09，k3=0.000 1。

4 跳躍移動(dòng)仿真校驗(yàn)

為驗(yàn)證探測(cè)器基于反作用飛輪進(jìn)行跳躍移動(dòng)方案的可行性，需要對(duì)其跳躍過程進(jìn)行仿真。假設(shè)探測(cè)器飛輪通過施加控制使探測(cè)器起跳后保持勻速前進(jìn)，仿真條件設(shè)置如下。

1)小行星表面環(huán)境及物理參數(shù)。根據(jù)現(xiàn)有的小行星環(huán)境，假設(shè)小行星表面重力加速度g=0.001 m/s2，逃逸速度vescape=2.5 m/s，探測(cè)器與小行星表面摩擦系數(shù)μ=0.45，靜摩擦系數(shù)μ0=0.47，恢復(fù)系數(shù)e=0.95，剛度系數(shù)K=0.5 GPa。

2)探測(cè)器的基本參數(shù)。質(zhì)量為5 kg，轉(zhuǎn)動(dòng)慣量I=diag([0.055 0.055 0.055])kg·m2，探測(cè)器邊長(zhǎng)2l=0.4 m。沿探測(cè)器3個(gè)飛輪轉(zhuǎn)軸施加的最大飛輪力矩大小Tc=0.01 N·m。

3)仿真初始條件。探測(cè)器初始角速度ω0=[0 0 0]Trad/s，反作用飛輪初始轉(zhuǎn)速Ω=[0 0 0]Trad/s，探測(cè)器初始姿態(tài)參數(shù)Q=[1 0 0 0]T。

圖4為通過飛輪加速起跳后的探測(cè)器質(zhì)心運(yùn)動(dòng)軌跡。由圖可見，對(duì)反作用飛輪施加控制，可實(shí)現(xiàn)探測(cè)器在小行星表面跳躍行走。行走過程仿真時(shí)間為800 s，得到探測(cè)器的起跳時(shí)間為3.16 s，跳躍的最大高度為1.33 m。

圖4 探測(cè)器質(zhì)心運(yùn)動(dòng)軌跡Fig.4 Trajectory of probe centroid

探測(cè)器質(zhì)心高度變化如圖5所示。從圖中可以清晰看出探測(cè)器每一次跳躍的過程。探測(cè)器在起跳后達(dá)到了最高跳躍高度，繼續(xù)跳躍移動(dòng)后，其滯空時(shí)間與跳躍高度逐漸減小，直至達(dá)到穩(wěn)定的前向跳躍。

圖5 探測(cè)器質(zhì)心高度變化Fig.5 Variation of probe’s centroid height

圖6 探測(cè)器質(zhì)心速度變化Fig.6 Variation of probe’s centroid velocity

探測(cè)器各方向速度分量隨時(shí)間的變化情況如圖6所示。圖中：水平方向速度分量vx與vy的變化曲線形狀為階梯形。由于跳躍探測(cè)器在空中飛行時(shí)只受到法向重力作用，切向受力為0，因此速度大小和方向都未發(fā)生改變。但當(dāng)探測(cè)器與地面發(fā)生碰撞時(shí)，在水平方向上受到了摩擦力的作用，使質(zhì)心水平速度在短時(shí)間內(nèi)發(fā)生了較大變化；豎直方向速度分量vz的曲線為鋸齒形，這是因?yàn)樵诳罩酗w行時(shí)探測(cè)器受到的重力作用、與地面碰撞時(shí)受到豎直方向接觸力的作用導(dǎo)致的。由于與地面碰撞的時(shí)間很短，因此每次碰撞時(shí)質(zhì)心速度和角速度均會(huì)出現(xiàn)突變。

由以上結(jié)果可看出，小行星探測(cè)器可以僅通過飛輪在小行星表面微重力的環(huán)境下實(shí)現(xiàn)跳躍移動(dòng)。對(duì)跳躍時(shí)的探測(cè)器姿態(tài)施加控制，可維持或改變跳躍方向，達(dá)到向不同方向移動(dòng)的目的。由于每次跳躍的距離及跳躍高度相對(duì)穩(wěn)定，使用反作用飛輪與使用彈跳機(jī)構(gòu)的跳躍移動(dòng)相比，可以更好地對(duì)探測(cè)器的飛行姿態(tài)及跳躍過程施加控制，更適合于微重力環(huán)境下的小行星探測(cè)任務(wù)的執(zhí)行。

5 結(jié)束語

本文針對(duì)小行星探測(cè)器跳躍行走方案的動(dòng)力學(xué)和控制策略開展研究。分析了國內(nèi)外小行星跳躍探測(cè)器的發(fā)展現(xiàn)狀，建立了跳躍探測(cè)器的動(dòng)力學(xué)模型，給出了反作用飛輪在探測(cè)器跳躍移動(dòng)中所需提供的最小力矩，設(shè)計(jì)了探測(cè)器姿態(tài)控制律，進(jìn)行了仿真檢驗(yàn)。結(jié)果表明：探測(cè)器能進(jìn)行完整的姿態(tài)控制過程，可以執(zhí)行小行星表面的探測(cè)任務(wù)。對(duì)跳躍探測(cè)器跳躍行走過程進(jìn)行了仿真，探測(cè)器通過反作用飛輪加速提供的力矩產(chǎn)生跳躍運(yùn)動(dòng)，驗(yàn)證了方案的可行性?？紤]到小行星表面的地形變化及非均勻的重力場(chǎng)環(huán)境，未來將研究跳躍探測(cè)器在不同地形環(huán)境及非均勻重力場(chǎng)下的跳躍運(yùn)動(dòng)，針對(duì)不同地形下所采取的跳躍移動(dòng)策略及相應(yīng)的控制方法進(jìn)行設(shè)計(jì)，以滿足不同條件下小行星著陸探測(cè)任務(wù)的需求。