胡冬生，王俊峰，張普卓

（1. 中國(guó)運(yùn)載火箭技術(shù)研究院研究發(fā)展中心，北京，100076；2. 北京宇航系統(tǒng)工程研究所，北京，100076）

基于分解策略的SSO任務(wù)發(fā)射軌道優(yōu)化研究

胡冬生1，王俊峰2，張普卓2

（1. 中國(guó)運(yùn)載火箭技術(shù)研究院研究發(fā)展中心，北京，100076；2. 北京宇航系統(tǒng)工程研究所，北京，100076）

700 km太陽(yáng)同步軌道為運(yùn)載火箭彈道設(shè)計(jì)中典型的發(fā)射任務(wù)。采用軌跡分解優(yōu)化策略對(duì)此類任務(wù)的發(fā)射軌道進(jìn)行研究，分析了中間軌道近地點(diǎn)高度、入軌點(diǎn)真近點(diǎn)角兩個(gè)約束變量對(duì)火箭運(yùn)載能力的影響，并重點(diǎn)結(jié)合氣動(dòng)熱流限制對(duì)約束變量的取值進(jìn)行了計(jì)算分析。采用的研究思路和得出的結(jié)論對(duì)于運(yùn)載火箭上升段彈道的優(yōu)化設(shè)計(jì)以及提升運(yùn)算收斂性和設(shè)計(jì)效率均有一定的意義。

分解優(yōu)化策略；發(fā)射軌道；約束變量；氣動(dòng)熱流

0 引 言

700 km太陽(yáng)同步軌道（Sun Synchronous Orbit，SSO），廣泛應(yīng)用于地球資源勘察、氣象/災(zāi)害監(jiān)測(cè)、軍事偵察等對(duì)地觀測(cè)領(lǐng)域，很多情況下運(yùn)載火箭的設(shè)計(jì)基于700 km SSO運(yùn)載能力作為技術(shù)指標(biāo)。對(duì)于該類發(fā)射任務(wù)，運(yùn)載火箭上升段彈道往往采用有中間軌道的發(fā)射彈道，通過(guò)無(wú)動(dòng)力滑行，利用重力來(lái)改變速度方向，進(jìn)而提升火箭性能，解決火箭動(dòng)力飛行時(shí)間短的問(wèn)題[1,2]；如此帶來(lái)設(shè)計(jì)變量的增多，彈道設(shè)計(jì)及優(yōu)化復(fù)雜?；鸺w行過(guò)程包含上升段、中間無(wú)動(dòng)力轉(zhuǎn)移段和軌道圓化段，各個(gè)飛行段不同的彈道特性增強(qiáng)了整個(gè)問(wèn)題的非線性，算法收斂變得更加困難[3]。在這類多軌跡段的彈道優(yōu)化問(wèn)題的求解上，軌跡分解優(yōu)化策略是一種有效的手段。本文應(yīng)用分解優(yōu)化策略對(duì)SSO任務(wù)的發(fā)射軌道進(jìn)行分析研究。

1 分解優(yōu)化策略設(shè)計(jì)思路

1.1 分解優(yōu)化策略概述

分解優(yōu)化策略是將整個(gè)任務(wù)分為不同的飛行段，如上升段、再入段、滑行段等，每個(gè)飛行段通過(guò)設(shè)置不同的約束目標(biāo)值產(chǎn)生一個(gè)滿秩的子問(wèn)題，可以應(yīng)用迭代法求解。同時(shí)，這些子問(wèn)題中的約束又作為主問(wèn)題中的獨(dú)立變量，結(jié)合整個(gè)問(wèn)題的優(yōu)化目標(biāo)應(yīng)用梯度投影法等算法對(duì)這些獨(dú)立變量進(jìn)行優(yōu)化。在整個(gè)問(wèn)題的求解中，第1級(jí)是各個(gè)子問(wèn)題的迭代求解，第2級(jí)是主問(wèn)題的優(yōu)化。采用兩級(jí)分解優(yōu)化策略后，問(wèn)題的收斂大大加快，提高了優(yōu)化效率。

Petersen等應(yīng)用兩級(jí)軌跡分解算法對(duì)執(zhí)行地球同步轉(zhuǎn)移軌道（Geostationary transfer oribt，GTO）任務(wù)的火箭發(fā)射彈道進(jìn)行了研究，將整個(gè)彈道分解為上升段、軌道轉(zhuǎn)移段及軌道圓化段，并引入衛(wèi)星質(zhì)量、上升段關(guān)機(jī)速度及轉(zhuǎn)移軌道傾角作為子問(wèn)題的目標(biāo)約束及主問(wèn)題的獨(dú)立變量，進(jìn)而以運(yùn)載能力為目標(biāo)進(jìn)行優(yōu)化得到約束變量的值[3]；Beltracchi等應(yīng)用分解策略對(duì)地球靜止軌道（Geostationary Earth Orbit，GEO）任務(wù)開(kāi)展分析，將整個(gè)彈道優(yōu)化問(wèn)題分解成助推級(jí)及上面級(jí)兩個(gè)飛行段，以中間停泊軌道參數(shù)作為子問(wèn)題的目標(biāo)約束（近地點(diǎn)高度、遠(yuǎn)地點(diǎn)高度和近地點(diǎn)輻角），對(duì)子問(wèn)題和主問(wèn)題均開(kāi)展優(yōu)化運(yùn)算，通過(guò)優(yōu)化停泊軌道參數(shù)提升GEO任務(wù)的運(yùn)載能力，同時(shí)針對(duì)停泊軌道參數(shù)對(duì)運(yùn)載能力的靈敏度進(jìn)行計(jì)算分析，以獲得不同參數(shù)取值對(duì)運(yùn)載能力的影響[4]。

1.2 SSO任務(wù)發(fā)射軌道分解策略設(shè)計(jì)模型

對(duì)于三級(jí)運(yùn)載火箭（三級(jí)可二次啟動(dòng)）執(zhí)行700 km SSO任務(wù)的彈道設(shè)計(jì)而言，其飛行過(guò)程為：火箭起飛至三級(jí)一次關(guān)機(jī)進(jìn)入中間轉(zhuǎn)移軌道，滑行一段時(shí)間后，在700 km高度左右二次點(diǎn)火進(jìn)行軌道圓化，最終進(jìn)入目標(biāo)軌道。圖1為運(yùn)載火箭執(zhí)行700 km SSO任務(wù)發(fā)射軌道示意。

圖1 運(yùn)載火箭執(zhí)行700km SSO任務(wù)發(fā)射軌道示意

通常除了發(fā)射方位角可通過(guò)與軌道傾角的關(guān)系進(jìn)行獨(dú)立設(shè)計(jì)外，還有一級(jí)最大攻角αm、二級(jí)程序角斜率dφcx2、三級(jí)一次飛行段程序角斜率dφcx31、三級(jí)一次關(guān)機(jī)時(shí)間tk31、三級(jí)滑行時(shí)間tc、三級(jí)二次點(diǎn)火時(shí)刻程序角 φcx32、三級(jí)二次飛行段程序角斜率 dφcx32、三級(jí)二次工作時(shí)間tk32共8個(gè)設(shè)計(jì)變量，而目標(biāo)約束僅有3個(gè)：700 km SSO入軌點(diǎn)的v，r，θ。設(shè)計(jì)變量個(gè)數(shù)遠(yuǎn)遠(yuǎn)多于約束個(gè)數(shù)。

一般情況下，三級(jí)二次工作時(shí)間較短，該段程序角斜率對(duì)結(jié)果影響很小，可認(rèn)為該段按常值程序角飛行，即dφcx32=0。對(duì)于其余的設(shè)計(jì)變量，可以根據(jù)分解策略引入若干中間約束，將全飛行段分為3個(gè)小段，每個(gè)小段的設(shè)計(jì)變量與約束個(gè)數(shù)相等，從而應(yīng)用迭代法求解方程組的方式來(lái)求出設(shè)計(jì)變量值。引入的3個(gè)約束分別為一級(jí)關(guān)機(jī)點(diǎn)當(dāng)?shù)貜椀纼A角 θk、中間軌道近地點(diǎn)高度hp及中間軌道入軌點(diǎn)位置f（即真近點(diǎn)角），運(yùn)載能力僅與這3個(gè)參數(shù)有關(guān)，可表示成三者的函數(shù)。同時(shí)，添加了 1個(gè)不等式約束：末級(jí)剩余推進(jìn)劑量mpr≥0 kg。因此，彈道優(yōu)化的目的即為尋找一級(jí)關(guān)機(jī)點(diǎn)θk、中間軌道hp及入軌點(diǎn)f，使得火箭運(yùn)載能力最大。尤其需要開(kāi)展中間軌道 hp及入軌點(diǎn) f對(duì)火箭彈道設(shè)計(jì)的影響，給出合理取值范圍，以有效指導(dǎo)發(fā)射軌道的設(shè)計(jì)和優(yōu)化。由此可將原問(wèn)題轉(zhuǎn)化為 1個(gè)主問(wèn)題和 3個(gè)子問(wèn)題，如圖2所示。

圖2 三級(jí)運(yùn)載火箭執(zhí)行SSO任務(wù)分解策略優(yōu)化方法

1.3 仿真計(jì)算效果

以某三級(jí)液體運(yùn)載火箭為例，進(jìn)行彈道優(yōu)化計(jì)算和分析。彈道計(jì)算中采用橢球形地球模型，考慮J2項(xiàng)引力攝動(dòng)，所用彈道計(jì)算模型見(jiàn)文獻(xiàn)[5]。

采用傳統(tǒng)方法求解時(shí)，算法對(duì)初值的設(shè)置較為敏感，常出現(xiàn)運(yùn)算發(fā)散需要修改初值并重新啟動(dòng)運(yùn)行的情況，且在數(shù)十次迭代以內(nèi)難以收斂。分解優(yōu)化策略雖然增加了求解問(wèn)題的數(shù)量，但每個(gè)子問(wèn)題的求解均比較簡(jiǎn)單，且僅需迭代4次左右即可完成收斂，子問(wèn)題的順利求解加快了主問(wèn)題的求解速度。相對(duì)而言，分解優(yōu)化策略大大增加了求解的可靠性和魯棒性[6]。

2 優(yōu)化變量對(duì)火箭性能的影響分析

2.1 不同近地點(diǎn)高度及真近點(diǎn)角情況下的運(yùn)載能力

為計(jì)算分析不同近地點(diǎn)高度及真近點(diǎn)角對(duì)運(yùn)載能力的影響，分別計(jì)算近地點(diǎn)高度取值區(qū)間為20～250 km情況下，不同真近點(diǎn)角時(shí)（對(duì)應(yīng)中間軌道不同入軌點(diǎn)）對(duì)應(yīng)的運(yùn)載能力。在每次計(jì)算中，暫不考慮落區(qū)約束，設(shè)定近地點(diǎn)高度及真近點(diǎn)角為固定值，優(yōu)化一級(jí)關(guān)機(jī)點(diǎn)當(dāng)?shù)貜椀纼A角，使火箭運(yùn)載能力最大。當(dāng)近地點(diǎn)高度小于80 km時(shí)（含亞軌道），人為設(shè)定火箭進(jìn)入中間軌道時(shí)高度應(yīng)大于80 km，且在近地點(diǎn)以后入軌。

近地點(diǎn)高度為80～250 km情況下運(yùn)載能力隨真近點(diǎn)角變化曲線見(jiàn)圖3。為便于比較，將計(jì)算出的運(yùn)載能力統(tǒng)一除以某固定值，形成無(wú)量綱運(yùn)載能力。

圖3 不同近地點(diǎn)高度運(yùn)載能力曲線

由圖3可以看出：

a）運(yùn)載能力隨著近地點(diǎn)高度的降低而增加。近地點(diǎn)高度降低時(shí)，火箭飛行過(guò)程中的重力損失減小，有助于運(yùn)載能力的提升。

b）對(duì)于相同近地點(diǎn)高度，真近點(diǎn)角的取值存在一個(gè)最優(yōu)值，使得運(yùn)載能力最大；該最優(yōu)值位于近地點(diǎn)附近，且隨著近地點(diǎn)高度的升高而增大（真近點(diǎn)角為負(fù)值表示火箭在近地點(diǎn)之前進(jìn)入中間軌道）；近地點(diǎn)附近入軌可基本保證火箭按照霍曼變軌的方式實(shí)現(xiàn)軌道轉(zhuǎn)移，在充分利用重力轉(zhuǎn)彎改變速度方向的同時(shí)抬高飛行高度，從而達(dá)到性能優(yōu)化。

c）對(duì)于圖3的曲線簇，同一運(yùn)載能力可能與不同曲線或同一曲線存在多個(gè)交點(diǎn)，因此可以通過(guò)設(shè)置不同近地點(diǎn)高度和真近點(diǎn)角組合來(lái)達(dá)到相同運(yùn)載能力。

近地點(diǎn)高度為20～80 km情況下運(yùn)載能力隨真近點(diǎn)角變化曲線見(jiàn)圖 4。此時(shí)，運(yùn)載能力隨真近點(diǎn)角的變化曲線與大于80 km時(shí)的情況不同，運(yùn)載能力最優(yōu)值出現(xiàn)在近地點(diǎn)以后入軌情況（對(duì)應(yīng)真近點(diǎn)角40～60°）?？梢?jiàn)，當(dāng)近地點(diǎn)高度過(guò)小時(shí)，受到在近地點(diǎn)以后入軌及入軌高度的影響，火箭不能充分利用霍曼變軌來(lái)達(dá)到性能優(yōu)化，運(yùn)載能力反而大幅降低。因此，若對(duì)火箭在中間軌道的滑行時(shí)間以及測(cè)控條件沒(méi)有約束，則應(yīng)盡量避免采用近地點(diǎn)高度過(guò)低的中間軌道，尤其是亞軌道。

圖4 不同近地點(diǎn)高度運(yùn)載能力曲線（20～80km）

2.2 最低回落高度及氣動(dòng)熱流變化情況

在彈道計(jì)算中，當(dāng)中間軌道入軌點(diǎn)高度過(guò)低或入軌點(diǎn)彈道傾角過(guò)小時(shí)，火箭飛行彈道會(huì)出現(xiàn)一個(gè)“回落”的階段。一方面，入軌點(diǎn)高度過(guò)低會(huì)造成入軌點(diǎn)速度變大，在各級(jí)推重比及發(fā)動(dòng)機(jī)工作時(shí)間一定的情況下，火箭需通過(guò)拉低彈道來(lái)增大飛行速度；另一方面，入軌點(diǎn)彈道傾角為負(fù)值時(shí)，則必然需要通過(guò)飛行高度的壓低來(lái)滿足入軌要求。在彈道“回落”的階段飛行高度有所降低，之后高度持續(xù)升高，因而存在一個(gè)“最低回落高度”。中間軌道近地點(diǎn)高度80 km、入軌點(diǎn)真近點(diǎn)角20°時(shí)火箭飛行高度隨時(shí)間變化曲線見(jiàn)圖5。

圖5 飛行高度隨時(shí)間變化曲線

由圖5可以看出，最低回落高度出現(xiàn)在567.468 s、高度約90.7 km處。

高度回落帶來(lái)的影響是：由于此時(shí)衛(wèi)星整流罩已分離（整流罩一般在二級(jí)飛行段分離，分離高度約100～120 km之間），衛(wèi)星裸露在外界空間，仍將受到稀薄大氣的熱流作用。高度降低后，大氣密度增加，而同時(shí)火箭飛行速度相比整流罩分離時(shí)增大很多，這樣就有可能使熱流快速回升，對(duì)衛(wèi)星產(chǎn)生不利影響。計(jì)算表明，在相同近地點(diǎn)高度情況下，真近點(diǎn)角越小，熱流再次回升的值越大，衛(wèi)星的飛行環(huán)境越惡劣。

在最低回落高度時(shí)刻，大氣還具有一定密度，且火箭飛行速度始終在持續(xù)增大。根據(jù)熱流計(jì)算公式，熱流再次回升至最大點(diǎn)（可稱為“最大回升熱流”）往往在最低回落高度對(duì)應(yīng)時(shí)間點(diǎn)之后。如在近地點(diǎn)高度80 km、真近點(diǎn)角20°，最大回升熱流點(diǎn)在最低回落高度點(diǎn)之后8 s，如圖6所示。

圖6 火箭飛行熱流曲線

為了將最大回升熱流限制在衛(wèi)星可承受的范圍內(nèi)（如不大于1 135 W/m2），需要將最低回落高度抬高至一定值，可通過(guò)增大真近點(diǎn)角來(lái)實(shí)現(xiàn)（其對(duì)應(yīng)的真近點(diǎn)角可稱為“臨界真近點(diǎn)角”），或者抬高中間軌道近地點(diǎn)高度，計(jì)算結(jié)果見(jiàn)表1。表1表明，最低回落高度應(yīng)不小于132 km，才能滿足衛(wèi)星熱流要求，且運(yùn)載能力隨著近地點(diǎn)高度的增大而有所提升。對(duì)于近地點(diǎn)高度150 km以上的情況，最優(yōu)運(yùn)載能力在0°真近點(diǎn)角附近，此時(shí)高度回落及熱流回升均不明顯，均能滿足衛(wèi)星熱流約束。

表1 不同近地點(diǎn)高度對(duì)應(yīng)的臨界真近點(diǎn)角及相關(guān)參數(shù)

2.3 近地點(diǎn)高度及真近點(diǎn)角的選擇

為了滿足熱流限制，最低回落高度應(yīng)不小于一定值（如132 km），則應(yīng)對(duì)近地點(diǎn)高度及真近點(diǎn)角的取值做一定限制。結(jié)合前面計(jì)算結(jié)果，可根據(jù)中間軌道近地點(diǎn)高度分兩種情況進(jìn)行討論：

a）近地點(diǎn)高度低于150 km（含亞軌道）。

若考慮熱流約束，則運(yùn)載能力相比2.1節(jié)中的最優(yōu)值將有不同程度的降低。以近地點(diǎn)高度80 km為例，運(yùn)載能力將由最優(yōu)值0.941降至0.892以下，甚至小于近地點(diǎn)高度 150 km、真近點(diǎn)角為 0°時(shí)的運(yùn)載能力0.917。

b）近地點(diǎn)高度高于150 km。

在近地點(diǎn)入軌基本可保證運(yùn)載能力最大，且最低回落高度不至于過(guò)低、回升熱流不至于過(guò)大。當(dāng)近地點(diǎn)高度為150 km時(shí)，回升熱流在800 W/m2以下。

根據(jù)上述計(jì)算及分析，考慮熱流因素后，較低近地點(diǎn)高度對(duì)應(yīng)的運(yùn)載能力可以由較高近地點(diǎn)高度、近地點(diǎn)附近入軌（即真近點(diǎn)角約為0°）來(lái)實(shí)現(xiàn)，甚至獲得更大的運(yùn)載能力。針對(duì)本文所提及的火箭構(gòu)型，考慮一定余量后，近地點(diǎn)高度選取200 km左右為宜，過(guò)高或過(guò)低，都將帶來(lái)運(yùn)載能力的損失。

可見(jiàn)，近地點(diǎn)高度與真近點(diǎn)角的選擇是一個(gè)綜合權(quán)衡的過(guò)程，不能僅僅局限于運(yùn)載能力的考慮給出一個(gè)最優(yōu)值，因此該類彈道優(yōu)化問(wèn)題不存在所謂絕對(duì)最優(yōu)值。工程設(shè)計(jì)中只能依據(jù)物理規(guī)律及各項(xiàng)約束條件選取一個(gè)大概的范圍，該取值能使問(wèn)題的求解得到簡(jiǎn)化，提高設(shè)計(jì)效率、滿足設(shè)計(jì)要求，同時(shí)又不至于改變問(wèn)題的本質(zhì)，火箭性能不至于大幅降低。

3 結(jié) 論

應(yīng)用分解優(yōu)化策略對(duì)700 km SSO任務(wù)的發(fā)射軌道進(jìn)行研究，提出一種求解彈道優(yōu)化問(wèn)題的思路，可較大程度上加快計(jì)算的收斂，提高計(jì)算魯棒性和設(shè)計(jì)效率，并對(duì)類似的其它彈道設(shè)計(jì)問(wèn)題均有借鑒價(jià)值。

重點(diǎn)分析了近地點(diǎn)高度、入軌點(diǎn)真近點(diǎn)角兩個(gè)中間約束變量對(duì)運(yùn)載能力的影響，并結(jié)合氣動(dòng)熱流限制得出了一些定性及定量的結(jié)論，對(duì)于SSO任務(wù)發(fā)射軌道設(shè)計(jì)具有一定的意義：

a）設(shè)置不同的中間軌道近地點(diǎn)高度及入軌點(diǎn)真近點(diǎn)角，可得到相同運(yùn)載能力，在彈道設(shè)計(jì)時(shí)應(yīng)綜合考慮氣動(dòng)熱流等其它約束條件來(lái)確定設(shè)計(jì)變量取值。

b）中間軌道采用霍曼變軌方式使火箭運(yùn)載能力得到較大提升，且近地點(diǎn)高度越低，對(duì)運(yùn)載能力的提升越有利。但當(dāng)近地點(diǎn)高度降低至一定值時(shí)，將使衛(wèi)星受到嚴(yán)重的熱流作用，對(duì)衛(wèi)星產(chǎn)生不利的影響。

c）中間軌道近地點(diǎn)高度設(shè)置在200 km左右，并在近地點(diǎn)入軌，可基本兼顧運(yùn)載能力優(yōu)化及熱流等約束。

應(yīng)注意的是，文中給出的優(yōu)化結(jié)果會(huì)隨著運(yùn)載火箭總體參數(shù)的不同而有所差異。在工程應(yīng)用中，應(yīng)根據(jù)測(cè)控、末級(jí)滑行時(shí)間、落區(qū)等實(shí)際情況綜合權(quán)衡，在滿足各項(xiàng)約束的前提下來(lái)考慮目標(biāo)的優(yōu)化，以選擇合適的中間軌道設(shè)計(jì)變量值。

[1]Arianspace. VEGA User’s Manual[M]. Paris: Programmes Division Technical Service, 2002.

[2]The Boeing Company. DELTA Ⅳ payload planners guide[M]. California: Space and Communications Group, 2000.

[3]Petersen F M, Cornick D E, and Brauer G L. A two-level trajectory decomposition algorithm featuring optimal intermedia target selection[J]. Journal of Spacecraft and Rockets, 1977，14(11): 676-682.

[4]Beltracchi T J. Decomposition approach to solving the all-up trajectory optimization problem[J]. Journal of Guidance, Control, and Dynamics, 1992, 15(3): 707-716.

[5]賈沛然, 陳克俊, 何力. 遠(yuǎn)程火箭彈道學(xué)[M]. 長(zhǎng)沙：國(guó)防科技大學(xué)出版社, 2004.

[6]Brauer G L, Cornick D E, Stevenson R. Capabilities and applications of the Program to optimize simulated trajectories program (POST) program summary document[R]. NASA CR-2770, 1977.

Launch Trajectory Optimization for SSO Mission Based on Decomposition Algorithm

Hu Dong-sheng1, Wang Jun-feng2, Zhang Pu-zhuo2
(1. R&D Center, China Academy of Launch Vehicle Technology, Beijing, 100076; 2. Beijing Institute of Astronautical Systems Engineering, Beijing, 100076)

In this paper the decomposition algorithm is proposed for the launch trajectory optimization of the SSO mission. The simulation is conducted to analyze the influence of two constraint variables to launch capability, such as perigee altitude and injection true anomaly of the intermediate orbit. Lastly, the value of the constraint variables is studied according to aerodynamic heating limitation. The research method and analysis results are important to trajectory optimization of launch vehicle, as well as the convergence of calculation and design efficiency.

Decomposition algorithm; Launch trajectory; Constraint variables; Aerodynamic heating

V412.4

A

1004-7182(2017)01-0010-05

10.7654/j.issn.1004-7182.20170103

2015-04-16；

2015-08-20

胡冬生（1982-），男，高級(jí)工程師，主要研究方向?yàn)檫\(yùn)載火箭彈道設(shè)計(jì)