趙 晟，王躍峰

(1.國(guó)防科技大學(xué)，湖南長(zhǎng)沙 410073;2.中國(guó)空間技術(shù)研究院，北京 100094)

深空探測(cè)返回是空間探索的一個(gè)重要發(fā)展方向，深空取樣返回是進(jìn)行深空探索，對(duì)外星球生命、物質(zhì)進(jìn)行研究的重要手段.中國(guó)已提出月球取樣返回計(jì)劃，隨著火星、金星等更遠(yuǎn)星體探測(cè)任務(wù)的執(zhí)行，外星體采樣返回將成為必然的發(fā)展趨勢(shì)［1］.深空探測(cè)返回具有以下特點(diǎn):再入速度高、不同任務(wù)的航程需求變化大、再入點(diǎn)散布大.

返回技術(shù)是集航天動(dòng)力學(xué)、空氣動(dòng)力學(xué)、防熱結(jié)構(gòu)學(xué)、控制學(xué)及導(dǎo)航學(xué)、火箭發(fā)動(dòng)機(jī)技術(shù)、測(cè)控技術(shù)、回收和著陸技術(shù)等學(xué)科的一門(mén)綜合性學(xué)科［2］.國(guó)內(nèi)目前掌握的返回技術(shù)為彈道式返回和標(biāo)準(zhǔn)軌道半升力返回.彈道式返回在返回過(guò)程中不產(chǎn)生升力或者產(chǎn)生較小的升力，但不對(duì)升力進(jìn)行控制的返回方式.標(biāo)準(zhǔn)軌道返回的返回器通過(guò)質(zhì)心偏置獲得升力，其標(biāo)準(zhǔn)軌道在地面發(fā)射前已經(jīng)注入，在再入大氣過(guò)程中進(jìn)行攝動(dòng)制導(dǎo)，保證著陸精度.此外，國(guó)外深空返回中也有半升力預(yù)測(cè)制導(dǎo)返回，它采用與標(biāo)準(zhǔn)軌道制導(dǎo)同樣的升力獲取方式，在返回過(guò)程中實(shí)時(shí)規(guī)劃飛行航跡，從而提高高速、長(zhǎng)航程再入飛行的著陸精度.

從國(guó)外的深空探測(cè)返回經(jīng)驗(yàn)看，均采用簡(jiǎn)單可靠的彈道式返回方式或有較高著陸精度的預(yù)測(cè)制導(dǎo)返回:美國(guó)的深空取樣返回器“星塵號(hào)”和“起源號(hào)”均采用彈道式再入;美國(guó)早期的探月飛行器“阿波羅系列”和蘇聯(lián)的探月飛行器“探測(cè)器6號(hào)”均采用預(yù)測(cè)制導(dǎo)返回［1］.

本文首先建立球冠倒錐外形飛行器在大氣層內(nèi)飛行的動(dòng)力學(xué)模型，然后對(duì)彈道式返回和預(yù)測(cè)制導(dǎo)返回兩種返回方式進(jìn)行飛行仿真，最后從返回飛行性能及其工程應(yīng)用技術(shù)要求兩個(gè)方面出發(fā)進(jìn)行工程應(yīng)用分析.為深空再入返回飛行的返回方案選擇提供指導(dǎo).

1 大氣層內(nèi)再入飛行動(dòng)力學(xué)建模

圖1 球冠倒錐飛行器

彈道式返回在再入過(guò)程中進(jìn)行自旋控制，保證以零升力返回;半升力返回通過(guò)調(diào)整傾側(cè)角，改變升力在縱、橫平面的分量達(dá)到控制飛行航跡的效果.

根據(jù)飛行器的外形及其在再入過(guò)程的飛行情況建立動(dòng)力學(xué)模型，存在以下假設(shè)條件:

1)飛行器始終以配平攻角飛行;

2)不考慮地球扁率;

3)地球萬(wàn)有引力滿(mǎn)足反平方關(guān)系.

根據(jù)動(dòng)力學(xué)理論，飛行器在大氣層再入飛行的動(dòng)力學(xué)模型為［3］

式中:ψ為方向角，指當(dāng)時(shí)飛行速度矢量與當(dāng)?shù)卣龞|方向的夾角;γ為飛行航跡傾角，指當(dāng)時(shí)飛行速度矢量與當(dāng)?shù)厮骄€(xiàn)的夾角;φ為滾轉(zhuǎn)角，由制導(dǎo)律確定;ω為地球自轉(zhuǎn)角速度，ω=7.2722×10－5rad/s;L為氣動(dòng)升力為氣動(dòng)阻力，D=6378.14km .

CM教科書(shū)有理數(shù)內(nèi)容中涉及學(xué)科相關(guān)聯(lián)背景素材的例題有5.08%，而RJ版教科書(shū)例題中未涉及與其它學(xué)科相關(guān)的背景素材，表明CM教科書(shū)更加注重引入與其它學(xué)科相關(guān)聯(lián)的背景知識(shí)，如數(shù)學(xué)在物理、化學(xué)和生物乃至體育等學(xué)科的應(yīng)用，促使學(xué)生意識(shí)到數(shù)學(xué)的廣泛用途．

2 彈道式和半升力預(yù)測(cè)制導(dǎo)返回飛行仿真

本文對(duì)高速再入問(wèn)題分別采用彈道式返回和預(yù)測(cè)制導(dǎo)返回兩種返回方式的飛行過(guò)程進(jìn)行仿真分析，飛行器以美國(guó)獵戶(hù)座外形為基礎(chǔ).采用彈道式返回時(shí)，假設(shè)質(zhì)心在縱向?qū)ΨQ(chēng)軸上，飛行過(guò)程升力系數(shù)為零.表1給出了獵戶(hù)座的外形尺寸和飛行器質(zhì)量［4］.

表1 航天器參數(shù)

2.1 彈道式返回仿真

根據(jù)彈道式返回飛行航程相對(duì)于再入角的靈敏度，以及實(shí)際再入姿態(tài)精度控制能力，選擇對(duì)－8°再入飛行過(guò)程仿真，以分析氣動(dòng)過(guò)載和熱流密度等性能參數(shù)特點(diǎn).

彈道式返回飛行持續(xù)200s左右，飛行航程為900km;從仿真曲線(xiàn)a)可以看出，彈道式返回由于再入角大(絕對(duì)值)，而且飛行過(guò)程沒(méi)有升力作用，所以飛行航跡較短，最終直接再入大氣;從速度變化曲線(xiàn)b)中可以看出，再入飛行過(guò)程隨著飛行高度下降，飛行密度增加，速度變化率進(jìn)一步增大，從而導(dǎo)致較大的氣動(dòng)過(guò)載和熱流密度峰值.曲線(xiàn)c)和d)分別為該返回過(guò)程對(duì)應(yīng)的氣動(dòng)過(guò)載荷熱流密度隨時(shí)間的變化曲線(xiàn).從圖中可知，本文仿真的返回器以－8°彈道式再入時(shí)，峰值氣動(dòng)過(guò)載可達(dá)25g，駐點(diǎn)熱流密度峰值達(dá)3.4×106W/m2.

圖2 彈道式返回仿真曲線(xiàn)

2.2 半升力預(yù)測(cè)制導(dǎo)返回仿真

預(yù)測(cè)制導(dǎo)返回在返回過(guò)程中需根據(jù)實(shí)際飛行情況、落點(diǎn)要求持續(xù)實(shí)時(shí)進(jìn)行傾側(cè)角制導(dǎo)計(jì)算，本文采用數(shù)值預(yù)測(cè)制導(dǎo)算法進(jìn)行制導(dǎo)計(jì)算.

半升力預(yù)測(cè)制導(dǎo)能夠?qū)崟r(shí)進(jìn)行制導(dǎo)計(jì)算，再入角變化率小，適合采用較淺再入，本文選擇－6°再入，4000km目標(biāo)航程問(wèn)題進(jìn)行仿真研究［5-7］.

從仿真曲線(xiàn)a)可以看出，航天器再入下降至一定高度(約60km)時(shí)，在升力控制下跳躍飛行來(lái)保證飛行航程需要，對(duì)應(yīng)于仿真曲線(xiàn)b)，在跳躍飛行段氣動(dòng)力較小，飛行器的速度變化率比較平緩，即實(shí)現(xiàn)平緩減速，避免了較大的氣動(dòng)過(guò)載.

仿真曲線(xiàn)b)～d)為對(duì)應(yīng)的航程、氣動(dòng)過(guò)載和熱流密度的變化曲線(xiàn).從圖中可以看出，通過(guò)航跡規(guī)劃，航天器的飛行航程達(dá)到約3960km，與目標(biāo)航程相差40km，是因?yàn)榭紤]到開(kāi)傘減速，仿真計(jì)算終止于10km.

與彈道式返回相比，該種長(zhǎng)距離再入飛行任務(wù)，采用預(yù)測(cè)制導(dǎo)返回具有較高的精度.從飛行過(guò)程的氣動(dòng)過(guò)載和熱流密度可以看出，通過(guò)航跡規(guī)劃，氣動(dòng)峰值分散在前后兩部分，最大氣動(dòng)過(guò)載僅為4.2g左右;熱流密度峰值出現(xiàn)在第一次下降飛行段，最大峰值為2×106W/m2.

3 彈道式和半升力預(yù)測(cè)制導(dǎo)返回工程應(yīng)用分析

通過(guò)上述對(duì)彈道式返回和預(yù)測(cè)制導(dǎo)返回的仿真分析比較，結(jié)合工程實(shí)際，可進(jìn)行深空彈道式返回和預(yù)測(cè)制導(dǎo)返回的對(duì)比.

3.1 彈道式返回飛行優(yōu)點(diǎn)及缺點(diǎn)分析:

彈道式返回在工程實(shí)際應(yīng)用中具有以下優(yōu)點(diǎn):

1)彈道式返回采用旋成體，外形設(shè)計(jì)簡(jiǎn)單;

2)彈道式返回飛行過(guò)程中，始終以恒定的角速度自旋，返回器外部受熱均勻，熱防護(hù)設(shè)計(jì)簡(jiǎn)單;

3)彈道式返回過(guò)程不需要制導(dǎo)，所以返回器質(zhì)量較輕，在滿(mǎn)足再入條件之后，飛行安全性高;

4)技術(shù)成熟度水平高，工程應(yīng)用實(shí)現(xiàn)的成本低.

彈道式返回在任務(wù)適應(yīng)性方面存在一些不足，可以歸納為以下幾點(diǎn):

1)彈道式返回航程不可調(diào)，飛行航跡相對(duì)再入角非常敏感，所以要想獲得較小的落點(diǎn)散布需要再入點(diǎn)飛行狀態(tài)誤差很小，并且為減小再入飛行過(guò)程中氣動(dòng)不確定性和風(fēng)等影響，需選擇較大的再入角，來(lái)保證飛行任務(wù)成功;

2)彈道式返回的再入角大，而且再入飛行過(guò)程中沒(méi)有升力可利用，所以飛行航跡較陡，在飛行過(guò)程中會(huì)產(chǎn)生較大的氣動(dòng)過(guò)載和熱流密度，這對(duì)防熱材料要求提高.并且該種氣動(dòng)過(guò)載環(huán)境不適用于載人返回等低過(guò)載深空返回任務(wù).

3.2 半升力預(yù)測(cè)制導(dǎo)返回飛行優(yōu)點(diǎn)及缺點(diǎn)分析:

半升力預(yù)測(cè)制導(dǎo)返回作為一種新的返回方式，有許多新的優(yōu)點(diǎn):

1)預(yù)測(cè)制導(dǎo)返回可以在軌規(guī)劃航跡，所以具有較大的任務(wù)適應(yīng)性;

2)返回飛行過(guò)程，利用升力可以調(diào)整航跡，拓寬再入走廊，從而降低對(duì)再入點(diǎn)精度的控制要求;

3)可以在較淺的再入條件下，進(jìn)行精確制導(dǎo)返回，從而大大降低氣動(dòng)過(guò)載，有利于深空載人等任務(wù)的展開(kāi).同時(shí)，峰值熱力密度較小，對(duì)防熱材料選用要求低;

4)通過(guò)跳躍飛行，可以有效增加飛行航程，有利于更廣泛的選擇落區(qū).

半升力預(yù)測(cè)制導(dǎo)返回在進(jìn)行工程應(yīng)用時(shí)相比與彈道式返回有一些缺點(diǎn):

1)預(yù)測(cè)制導(dǎo)返回的制導(dǎo)方案設(shè)計(jì)復(fù)雜，在工程實(shí)際應(yīng)用中，需要與控制系統(tǒng)相結(jié)合，增加了返回器設(shè)計(jì)的難度;

2)返回飛行持續(xù)時(shí)間長(zhǎng)，加熱總量大，所以需要增加防熱材料的使用，可能增加了返回器的質(zhì)量.

4 結(jié)論

從以上分析可知，深空彈道式返回可以簡(jiǎn)化設(shè)計(jì)，縮短研制周期，降低成本，適合于進(jìn)行深空返回技術(shù)的前期驗(yàn)證和無(wú)人深空探測(cè)任務(wù)的完成;半升力預(yù)測(cè)制導(dǎo)返回是有效解決深空高速再入返回問(wèn)題的途徑，對(duì)于載人或者其攜帶有對(duì)再入飛行過(guò)載有較高要求的有效載荷時(shí)，適宜采用半升力預(yù)測(cè)制導(dǎo)返回實(shí)現(xiàn).

［1］ 韓鴻碩，陳杰.21世紀(jì)國(guó)外深空探測(cè)發(fā)展計(jì)劃及進(jìn)展［J］.航天器工程，2008，17(3):2-22

［2］ 王希季.航天器進(jìn)入與返回技術(shù)［M］.北京:中國(guó)宇航出版社，1991

［3］ Frank J R.Re-entry vehicle dynamics［M］.New York:America Institute of Aeronautics and Astronatics Inc.，1984

［4］ Bairstow S H，Barton G H.Orion reentry guidance with extended range capability using PredGuid［C］.AIAA Guidance，Navigation and Control Conference and Exhibit，Hilton Head，South Carolina，August 20-23，2007

［5］ Graves C A，Harpold J C.Apollo experience report mission planning for Apollo entry［R］.NASA TN D-6725，19720013191，1972

［6］ 陸平，朱亮，敬忠良，等.探月返回跳躍式再入制導(dǎo)［C］.全國(guó)第十三屆空間及運(yùn)動(dòng)體控制技術(shù)學(xué)術(shù)會(huì)議，宜昌，2008

［7］ Chapman D R.An approximate analytical method for studying entry into planetary atmospheres［D］.Ames Aeronautical Laboratory Moffett Field，Calif.，1958