考慮液體晃動效應(yīng)的運載火箭等效轉(zhuǎn)動慣量對姿控系統(tǒng)穩(wěn)定裕度的影響

司 成，余延生，姚春燕

（1.上海宇航系統(tǒng)工程研究所，上海 201109；2.上海航天控制技術(shù)研究所，上海 200233）

0 引言

大型運載火箭目前多以液體推進劑作為動力，推進劑質(zhì)量占全箭總質(zhì)量約90%，液體晃動對火箭控制系統(tǒng)影響顯著，是動力學(xué)建模不可忽視的因素。工程上將液體運動等效為一定的力學(xué)模型，其中包含隨貯箱運動和晃動兩部分。前者用固定質(zhì)量塊表示，后者用彈簧－質(zhì)點表示，通過理論計算分析得到彈簧－質(zhì)點模型的特性參數(shù)，用于穩(wěn)定系統(tǒng)的設(shè)計。

總體在設(shè)計運載火箭總體參數(shù)時，須給出全箭繞縱軸的轉(zhuǎn)動慣量Jz和質(zhì)心位置Xz。傳統(tǒng)方法是將全部推進劑假設(shè)為固化物體以簡化設(shè)計，因簡化引起的誤差通過加大參數(shù)的偏差范圍加以考慮，同時將此偏差傳遞給穩(wěn)定系統(tǒng)設(shè)計，但較大的參數(shù)變化范圍增加了穩(wěn)定系統(tǒng)設(shè)計的難度。

現(xiàn)有資料表明，采用這種簡化的液態(tài)晃動處理方法后，全箭質(zhì)心位置與實際情形相差不大，但火箭液體貯箱的有效轉(zhuǎn)動慣量與液體固化時的轉(zhuǎn)動慣量之比可達15%～65%，誤差較明顯［1］。因此，穩(wěn)定系統(tǒng)參數(shù)設(shè)計時應(yīng)考慮推進劑晃動帶來的有效轉(zhuǎn)動慣量變化的影響。

針對某兩級火箭，本文定量分析了推進劑晃動引起的有效轉(zhuǎn)動慣量對箭體動力學(xué)參數(shù)的影響，同時根據(jù)火箭特點，在選取的特征秒點，仿真分析了對火箭穩(wěn)定系統(tǒng)裕度的影響。

1 液體晃動引起的等效轉(zhuǎn)動慣量變化

鑒于液體晃動對繞箭體縱軸（X軸）的轉(zhuǎn)動慣量Jx影響很小，本文討論僅限于對俯仰、偏航運動影響較大的繞橫軸（Y、Z軸）的轉(zhuǎn)動慣量Jz的變化。由于俯仰和偏航通道的對稱性，分析以俯仰通道為例。

可從運動機理分析液體晃動引起的箭體轉(zhuǎn)動慣量變化：淺箱轉(zhuǎn)動時，貯箱上下底對運動的約束作用較強，較多液體隨貯箱轉(zhuǎn)動，液體轉(zhuǎn)動慣量相對固化處理偏差不大；深箱轉(zhuǎn)動時，側(cè)壁對液體運動影響較大，較多液體隨貯箱轉(zhuǎn)動，液體轉(zhuǎn)動慣量下降不明顯；介于兩者之間的中等深度液體，上下底及側(cè)壁對液體運動的約束不強，轉(zhuǎn)動慣量下降較明顯［1］。

根據(jù)NASA的分析，在大部分液位高度，有效轉(zhuǎn)動慣量與固化液體轉(zhuǎn)動慣量相當(dāng)；但在少數(shù)液位狀態(tài)，兩者差別較大［2－3］。考慮液體晃動的貯箱等效轉(zhuǎn)動慣量與經(jīng)固化處理得到的轉(zhuǎn)動慣量隨液位高度的變化規(guī)律如圖1所示。圖中：h為液位高度／箱體；d為直徑。國內(nèi)某研究所針對液體有效轉(zhuǎn)動慣量的試驗數(shù)據(jù)支持了圖1的分析結(jié)論［4］。

1.1 液體晃動對一級轉(zhuǎn)動慣量Jz的影響

某火箭一級氧化劑箱、燃燒劑箱液位高度與直徑比為0～3。根據(jù)圖1，當(dāng)液體晃動高度h與貯箱直徑d相當(dāng)時，其等效轉(zhuǎn)動慣量下降最明顯，等效值約為固化計算值的19%。根據(jù)火箭相關(guān)數(shù)據(jù)，典型液位高度下氧化劑箱、燃燒劑箱等效轉(zhuǎn)動慣量變化分別見表1、2。

圖1 不同液位高度／箱體直徑比下等效轉(zhuǎn)動慣量與固化液體轉(zhuǎn)動慣量之比Fig.1 Relation between equivalent inertial moment and ratio of liquid depth／tank diameter

由圖1可知：在某些液位高度，液體晃動時貯箱等效轉(zhuǎn)動慣量與固化處理的比值較小，最小0.15左右，因此有必要分析充液貯箱對全箭轉(zhuǎn)動慣量的影響。

計算結(jié)果表明：在選取的典型液位高度上，考慮推進劑晃動，帶來的轉(zhuǎn)動慣量總計最大減小1.46%，小于參數(shù)偏差10%的要求。因此，用固化法計算獲得的全箭箭體轉(zhuǎn)動慣量在工程上有足夠的精確性。

究其原因，推進劑貯箱對全箭的慣量貢獻中包括繞自身縱軸項慣量及移軸項部分，后者由于距離全箭質(zhì)心距離較大而成為主要項，使液體晃動導(dǎo)致的轉(zhuǎn)動慣量減小的影響并不明顯。

1.2 液體晃動對二級轉(zhuǎn)動慣量Jz的影響

某運載火箭二級貯箱液位高度與直徑的比值為0～0.93。考慮火箭二級飛行過程中，加滿狀態(tài)下箱頂對液體晃動有較強的抑制作用，此時圖1理論分析結(jié)果適用性較差。因此二級分析考慮液位高度h／d＝0.77的35s時刻及h／d＝0.53的90s。此時液體晃動效應(yīng)較顯著，穩(wěn)定系統(tǒng)設(shè)計裕度較低。

表1 不同液位高度的氧箱等效轉(zhuǎn)動慣量Tab.1 Equivalent inertial moment under various liquid depth／tank diameter ratio（Y1－tank）

表2 燃箱等效轉(zhuǎn)動慣量Jz的變化Tab.2 Relation between equivalent inertial moment and the ratio of liquid depth／tank diameter（R1－tank）

計算所得的氧箱和燃箱等效轉(zhuǎn)動慣量變化分別見表3、4。

二級計算結(jié)果表明，液體晃動帶來的轉(zhuǎn)動慣量減小總計達全箭慣量的13.9%，超出了參數(shù)偏差10%的指標要求，對穩(wěn)定系統(tǒng)設(shè)計裕度會產(chǎn)生不利影響。

2 液體晃動對穩(wěn)定系統(tǒng)裕度的影響

由上述分析可知，液體晃動導(dǎo)致的一級箭體等效轉(zhuǎn)動慣量下降不明顯，二級箭體等效轉(zhuǎn)動慣量變化較大。

2.1 箭體穩(wěn)定系統(tǒng)動力學(xué)模型

全箭繞橫軸的轉(zhuǎn)動慣量變化，直接影響繞橫軸的轉(zhuǎn)動運動。箭體繞橫軸的轉(zhuǎn)動運動動力學(xué)模型為

式（1）～（7）為箭體運動動力學(xué)線化模型，穩(wěn)定系統(tǒng)以此為基礎(chǔ)，設(shè)計合適的校正網(wǎng)絡(luò)，使系統(tǒng)具有各項裕度。模型中，式（2）為力矩方程，各項力矩系數(shù)bi均為轉(zhuǎn)動慣量Jz的函數(shù)，當(dāng)轉(zhuǎn)動慣量變化時，力矩方程系數(shù)隨之而變，從而對系統(tǒng)設(shè)計的穩(wěn)定裕度產(chǎn)生影響。

表3 氧箱等效轉(zhuǎn)動慣量Jz的變化Tab.3 Equivalent inertial moment under various liquid depth／tank diameter ratio（Y2－tank）

表4 燃箱等效轉(zhuǎn)動慣量Jz的變化Tab.4 Equivalent inertial moment under various liquid depth／tank diameter ratio（R2－tank）

2.2 對穩(wěn)定系統(tǒng)設(shè)計裕度的影響

根據(jù)上述分析結(jié)果，選取二級35，90s液體晃動等效效應(yīng)影響最強的特征點，對系統(tǒng)剛體、晃動、彈性裕度的變化情況進行仿真。其中35s高頻段彈性裕度變化見表5。考慮晃動效應(yīng)后，轉(zhuǎn)動慣量的變化對系統(tǒng)裕度的影響主要集中在中、低頻段，而對高頻段的彈性裕度影響有限。

二級35s剛晃穩(wěn)定裕度仿真結(jié)果對見表6；二級90s剛晃穩(wěn)定裕度仿真結(jié)果見表7。

表5 考慮晃動效應(yīng)后二級35s高頻裕度Tab.5 Marge of system on high frequency in second stage at 35s

表6 考慮晃動效應(yīng)后二級35s剛晃穩(wěn)定裕度比較Tab.6 Marge of system on low－mid frequency in second stage at 35s

表7 考慮晃動效應(yīng)后二級90s剛晃穩(wěn)定裕度比較Tab.7 Marge of system on low－mid frequency in second stage at 90s

仿真結(jié)果表明：在液體晃動等效效應(yīng)較明顯的時刻，液體晃動對系統(tǒng)穩(wěn)定裕度有一定的影響，具體而言，對剛體低頻幅裕度（圓頻率0.6～0.7）影響較小；剛體高頻幅裕度（圓頻率約20）下降最大約2.2dB；晃動相裕度下降約10.7°（圓頻率5～6）。液體晃動的等效效應(yīng)對系統(tǒng)中、低頻段的幅、相裕度均產(chǎn)生一定的不利影響。

上述分析是在系統(tǒng)額定狀態(tài)下進行的。當(dāng)考慮其他參數(shù)偏差組合進行系統(tǒng)上、下限分析時，由于系統(tǒng)參數(shù)組合情況較惡劣，再加上液體晃動等效效應(yīng)對系統(tǒng)裕度的影響，35s時仿真結(jié)果，系統(tǒng)出現(xiàn)了晃動不穩(wěn)定的現(xiàn)象。

3 結(jié)束語

本文根據(jù)圓柱貯箱液體晃動試驗結(jié)果，分析了液體晃動對某型火箭繞橫軸的轉(zhuǎn)動慣量的影響，給出了考慮液體晃動和固化處理兩種處理方法全箭轉(zhuǎn)動慣量的變化情況，并在此基礎(chǔ)上，仿真分析了對姿控系統(tǒng)設(shè)計裕度的影響。結(jié)果表明：液體晃動的等效效應(yīng)對一級轉(zhuǎn)動慣量Jz的影響較小，對二級轉(zhuǎn)動慣量Jz的影響較大；進一步分析表明，等效轉(zhuǎn)動慣量對二級高頻穩(wěn)定裕度基本不產(chǎn)生影響，而系統(tǒng)的中、低頻穩(wěn)定裕度均有一定程度的下降。當(dāng)系統(tǒng)設(shè)計的剛體、晃動裕度較低時，應(yīng)當(dāng)關(guān)注液體晃動等效效應(yīng)對系統(tǒng)裕度的影響，防止姿控系統(tǒng)失穩(wěn)現(xiàn)象的發(fā)生。

［1］ 馬斌捷，劉 楨，林 宏.航天器貯箱液體有效轉(zhuǎn)動慣量［J］.導(dǎo)彈與航天運載技術(shù)，2011（2）：22－25.

［2］ ABRAMSON H N.The dynamic behavior of liquids in moving containers［R］.NASA SP－106，1966.

［3］ STOFAN A J，SUMMER I E.Experimental investigation of liquid sloshing in a scale－model centaur liquid－h(huán)ydrogen tank［R］.Lewis Research Center，1966.

［4］ 方良玉.圓柱形貯箱中的液體有效轉(zhuǎn)動慣量及其實驗研究［J］.強度與環(huán)境，1989，16（4）：27－32.

［5］ 導(dǎo)彈與航天叢書系列編輯委員會.控制系統(tǒng)（上）［M］.北京：宇航出版社，1989.