胡 齊，李 永，姚 燦，劉錦濤

(1.北京控制工程研究所，北京100190; 2.北京市高效能及綠色宇航推進(jìn)工程技術(shù)研究中心，北京100190)

大容量推進(jìn)劑貯箱液體晃動(dòng)性能試驗(yàn)*

胡 齊1，2，李 永1，2，姚 燦1，2，劉錦濤1，2

(1.北京控制工程研究所，北京100190; 2.北京市高效能及綠色宇航推進(jìn)工程技術(shù)研究中心，北京100190)

對(duì)某種內(nèi)帶推進(jìn)劑管理裝置(PMD)的大容量推進(jìn)劑貯箱內(nèi)液體晃動(dòng)性能進(jìn)行試驗(yàn).開(kāi)展不同充液比工況下空殼貯箱液體晃動(dòng)試驗(yàn)，分析比較晃動(dòng)試驗(yàn)結(jié)果與采用等效動(dòng)力學(xué)模型的液體晃動(dòng)理論計(jì)算結(jié)果，兩者結(jié)果一致吻合，試驗(yàn)系統(tǒng)可靠性和理論模型的正確性得到良好驗(yàn)證.然后開(kāi)展了不同充液比工況下內(nèi)帶PMD貯箱液體晃動(dòng)試驗(yàn)，試驗(yàn)結(jié)果表明貯箱內(nèi)液體晃動(dòng)性能受內(nèi)部PMD結(jié)構(gòu)影響明顯.該試驗(yàn)研究結(jié)果為運(yùn)載火箭和衛(wèi)星的姿態(tài)和軌道控制系統(tǒng)的設(shè)計(jì)優(yōu)化提供重要參考和數(shù)據(jù)支撐.

液體晃動(dòng);晃動(dòng)頻率;阻尼;大容量推進(jìn)劑貯箱;推進(jìn)劑管理裝置(PMD)

0 引言

本文針對(duì)某一種大容量推進(jìn)劑貯箱，該貯箱內(nèi)部推進(jìn)劑管理裝置(PMD)帶有防晃錐以及中間隔板，采用共振自由衰減晃動(dòng)試驗(yàn)方法和強(qiáng)迫晃動(dòng)試驗(yàn)方法開(kāi)展貯箱液體晃動(dòng)性能試驗(yàn)研究.

1 晃動(dòng)試驗(yàn)研究

1.1 試驗(yàn)研究對(duì)象

推進(jìn)劑貯箱液體晃動(dòng)試驗(yàn)的試驗(yàn)對(duì)象為模型貯箱，其容積為1170L，屬于大容量貯箱.該模型貯箱由有機(jī)玻璃內(nèi)殼和鋁合金框架組成，其內(nèi)殼輪廓尺寸與實(shí)際貯箱內(nèi)腔尺寸相同，而且內(nèi)部安放PMD.模型貯箱的PMD為實(shí)際貯箱PMD的簡(jiǎn)化結(jié)構(gòu)，包括中間隔板和防晃錐，模型貯箱由中間隔板隔開(kāi)分成上艙和下艙兩個(gè)部分，試驗(yàn)時(shí)其上艙與下艙液體相對(duì)獨(dú)立.模型貯箱坐標(biāo)如圖1所示，安裝時(shí)其z軸為豎直方向.

圖1 模型貯箱結(jié)構(gòu)示意圖Fig.1 Structural sketch of model tank

1.2 試驗(yàn)內(nèi)容

通過(guò)晃動(dòng)試驗(yàn)及參數(shù)測(cè)定來(lái)獲取相關(guān)參數(shù)，包括晃動(dòng)頻率、晃動(dòng)質(zhì)量、晃動(dòng)阻尼等.上述參數(shù)可通過(guò)自由衰減試驗(yàn)和強(qiáng)迫晃動(dòng)試驗(yàn)擬合得到.其中自由衰減試驗(yàn)是為了獲取晃動(dòng)頻率、晃動(dòng)質(zhì)心高度、晃動(dòng)阻尼比.強(qiáng)迫晃動(dòng)可以獲取得到晃動(dòng)質(zhì)量、總質(zhì)量等參數(shù).試驗(yàn)內(nèi)容及工況詳見(jiàn)表1.

1.3 試驗(yàn)方法

初中英語(yǔ)老師應(yīng)全面深入的了解每一位學(xué)生英語(yǔ)學(xué)習(xí)的實(shí)際情況，知道每一位學(xué)生的英語(yǔ)學(xué)習(xí)習(xí)慣以及英語(yǔ)學(xué)習(xí)態(tài)度，根據(jù)班級(jí)內(nèi)學(xué)生總數(shù)和學(xué)生英語(yǔ)成績(jī)的差距來(lái)確定英語(yǔ)分層數(shù)，堅(jiān)持學(xué)生自由選擇和老師微調(diào)整相結(jié)合的原則來(lái)將班級(jí)內(nèi)的學(xué)生靈活劃分到相應(yīng)的學(xué)習(xí)小組之中。英語(yǔ)老師要根據(jù)不同層次學(xué)生整體情況制定合理的學(xué)習(xí)計(jì)劃和學(xué)習(xí)目標(biāo)。

進(jìn)行強(qiáng)迫晃動(dòng)試驗(yàn)時(shí)采用定位移激勵(lì)，啟動(dòng)振動(dòng)臺(tái)，調(diào)節(jié)振動(dòng)頻率，待運(yùn)動(dòng)平穩(wěn)后開(kāi)始記錄試件的晃動(dòng)位移、力和力矩的時(shí)間歷程，記錄50個(gè)振動(dòng)周期后停止;根據(jù)記錄的時(shí)間歷程數(shù)據(jù)計(jì)算位移、力和力矩的幅值、相位和頻率，并估計(jì)頻響函數(shù)值;改變振動(dòng)頻率后重復(fù)進(jìn)行記錄，直至完成所有頻率點(diǎn)的試驗(yàn).

進(jìn)行自由衰減晃動(dòng)試驗(yàn)時(shí)，啟動(dòng)振動(dòng)臺(tái)，調(diào)整振動(dòng)頻率至液體共振頻率，待晃幅達(dá)到一定幅度后將方框與推拉臺(tái)脫離后鎖定，記錄晃動(dòng)力的時(shí)間歷程數(shù)據(jù)，計(jì)算晃動(dòng)阻尼系數(shù)和晃動(dòng)頻率，重復(fù)5～7次.

表1 液體晃動(dòng)試驗(yàn)內(nèi)容及工況Tab.1 Contents and conditions of liquid sloshing experiment

1.4 液體晃動(dòng)等效動(dòng)力學(xué)模型理論

對(duì)于平底圓筒內(nèi)的液體晃動(dòng)問(wèn)題是存在解析解，但由于模型貯箱明顯不屬于平底圓筒的構(gòu)型，本貯箱液體晃動(dòng)問(wèn)題是不存在解析解.對(duì)于模型貯箱豎直安裝且z軸向上進(jìn)行橫向(x方向)晃動(dòng)試驗(yàn)的情況，當(dāng)液面高度位于試驗(yàn)貯箱柱段位置時(shí)，按液面等效將貯箱等效為平底圓筒，或采用文獻(xiàn)［13］中球底圓筒的實(shí)驗(yàn)結(jié)果.當(dāng)液面高度位于下半球或上半球高度內(nèi)可采用文獻(xiàn)［13］中球形貯箱的近似解;也可按液面等效將貯箱等效為平底圓筒，位于下半球時(shí)按淺箱處理，液面高度取液體全部上移和全部下移的平均值，位于上半球時(shí)按深箱處理，液面高度取液體上表面，由此可推算出貯箱的等效動(dòng)力學(xué)模型參數(shù).

這種等效動(dòng)力學(xué)模型有兩方面的優(yōu)點(diǎn)，一是便于控制系統(tǒng)建立穩(wěn)定性分析模型，二是對(duì)于非光滑貯箱(如本模型貯箱內(nèi)帶PMD)，晃動(dòng)問(wèn)題沒(méi)有理論解，或理論解誤差較大，此時(shí)等效動(dòng)力學(xué)模型仍然適用，但是模型參數(shù)無(wú)法得到理論解，只能通過(guò)晃動(dòng)試驗(yàn)來(lái)獲取.

2 試驗(yàn)結(jié)果與分析

2.1 空殼模型貯箱晃動(dòng)試驗(yàn)

空殼模型貯箱在豎直安裝方式下沿橫向(x方向)進(jìn)行晃動(dòng)試驗(yàn)，試驗(yàn)介質(zhì)為水，如圖2所示.

圖2 空殼貯箱豎直安裝x向晃動(dòng)示意圖Fig.2 Sketch of slosh when upright setting about the tank without PMD in x direction

分析比較空殼模型貯箱的液體晃動(dòng)等效動(dòng)力學(xué)模型理論計(jì)算結(jié)果與試驗(yàn)結(jié)果，晃動(dòng)頻率、晃動(dòng)質(zhì)心高度、擬合的總質(zhì)量均隨液體充液比增加而增加，如圖3～5所示.但是晃動(dòng)質(zhì)量均隨充液比的增加呈現(xiàn)先增大后減小的趨勢(shì)，在充液比0.65附近時(shí)晃動(dòng)質(zhì)量達(dá)到最大值，如圖6所示.

根據(jù)空殼模型貯箱晃動(dòng)試驗(yàn)結(jié)果可知，試驗(yàn)結(jié)果與理論計(jì)算結(jié)果一致吻合，試驗(yàn)系統(tǒng)誤差在允許范圍之內(nèi)，試驗(yàn)系統(tǒng)可靠性和理論模型的正確性得到驗(yàn)證.

圖3 空殼貯箱晃動(dòng)頻率理論與試驗(yàn)結(jié)果對(duì)比Fig.3 Comparison between theoretic and experimental slosh frequencies in the tank without PMD

圖4 空殼貯箱晃動(dòng)質(zhì)心高度與充液比的關(guān)系Fig.4 Relation between slosh centroid height and filling ratio in the tank without PMD

圖5 空殼貯箱總質(zhì)量的擬合值、理論值與實(shí)際稱重結(jié)果的對(duì)比Fig.5 Comparison of overall weight between simulant theoretic and experimental results in the tank without PMD

圖6 空殼貯箱晃動(dòng)質(zhì)量理論與試驗(yàn)結(jié)果對(duì)比Fig.6 Comparison between theoretic and experimental slosh weight in the tank without PMD

2.2 帶PMD模型貯箱上艙晃動(dòng)試驗(yàn)

帶PMD模型貯箱在豎直安裝方式下沿橫向(x方向)進(jìn)行晃動(dòng)試驗(yàn)，如圖7所示.

圖7 帶PMD模型貯箱豎直安裝x向晃動(dòng)示意圖Fig.7 Sketch of slosh when upright setting about the tank with PMD in x direction

從晃動(dòng)頻率與充液比關(guān)系曲線可知，帶PMD貯箱的晃動(dòng)頻率范圍覆蓋0.5～1.1 Hz，空殼貯箱的晃動(dòng)頻率隨著充液比的增大單調(diào)增加，帶PMD貯箱的晃動(dòng)頻率隨充液比總體呈非單調(diào)變化趨勢(shì)，如圖8所示.由此可知，貯箱液體晃動(dòng)頻率受貯箱內(nèi)PMD結(jié)構(gòu)影響明顯.

充液比從0.35到0.85工況下，防晃錐產(chǎn)生箱底效應(yīng)，浸沒(méi)防晃錐區(qū)域液體的一階固有頻率較空殼貯箱的頻率值低.在0.4充液比時(shí)，在此液面以下的防晃錐椎體部分無(wú)開(kāi)孔，防晃錐起著類似縮小貯箱直徑的作用，使得液體頻率較0.35充液比時(shí)顯著增加;在0.45充液比時(shí)，液位處于貯箱中間柱段與PMD中間隔板連接處，液體可自由通過(guò)防晃錐上的開(kāi)孔;在0.85充液時(shí)，液體剛浸沒(méi)整個(gè)防晃錐，即0.85充液比以上液體位置高過(guò)防晃錐的錐頂;在0.95充液時(shí)，帶 PMD貯箱充混合液的試驗(yàn)值介于等效平底圓筒理論值和球形貯箱理論值之間.

圖8 晃動(dòng)頻率與充液比關(guān)系Fig.8 Relation between slosh frequency and filling ratio

液體晃動(dòng)擺點(diǎn)高度隨充液比總體呈減少趨勢(shì)，由于擺點(diǎn)高度與一階晃動(dòng)頻率有關(guān)，在0.45充液比附近擺點(diǎn)高度值較大，如圖9所示.

圖9 擺點(diǎn)高度與充液比關(guān)系Fig.9 Relation between sway height and filling ratio

由帶PMD貯箱自由衰減試驗(yàn)數(shù)據(jù)結(jié)果可知，阻尼呈現(xiàn)非單調(diào)變化趨勢(shì).在0.35充液比時(shí)，液體位于防晃錐根部位置，貯箱晃動(dòng)時(shí)液體受到防晃錐的影響相對(duì)較小，因此其阻尼值出現(xiàn)了減小的現(xiàn)象;在0.40低充液比時(shí)，液面之下的防晃錐上無(wú)開(kāi)孔，防晃錐對(duì)液體的阻礙作用十分明顯，因此試驗(yàn)阻尼值明顯增加.在0.45～0.85充液比時(shí)，液體受到防晃錐的錐角和開(kāi)孔的共同影響，阻尼隨著充液比的增加而減少.在0.89和0.95充液比時(shí)，液面位置已經(jīng)超過(guò)防晃錐的錐頂，防晃錐對(duì)液體晃動(dòng)的阻尼影響減少，0.95充液比時(shí)阻尼值最小，如圖10所示.

圖10 阻尼與充液比關(guān)系Fig.10 Relation between slosh damping and filling ratio

如圖11所示，帶PMD貯箱的晃動(dòng)質(zhì)量隨充液比的增加呈現(xiàn)先增大后減小的趨勢(shì)，空殼貯箱的晃動(dòng)質(zhì)量在0.65附近最大，帶PMD貯箱在0.76充液比附近時(shí)晃動(dòng)質(zhì)量最大.

圖11 空殼貯箱與帶PMD貯箱的晃動(dòng)質(zhì)量對(duì)比Fig.11 Comparison of slosh weight in the tank between without PMD and with PMD

如圖12所示，帶PMD貯箱的晃動(dòng)質(zhì)量占總質(zhì)量的比例隨著充液比的增加總體呈現(xiàn)先增加后減少的趨勢(shì)，晃動(dòng)質(zhì)量所占總質(zhì)量的比例在10%～40%之間，在充液比0.55～0.75之間晃動(dòng)質(zhì)量所占比例最大.

圖12 不同充液比條件下晃動(dòng)質(zhì)量占總質(zhì)量的比例關(guān)系Fig.12 Proportion of slosh weight in overall under different filling ratios

2.3 小結(jié)

1)試驗(yàn)結(jié)果與理論計(jì)算結(jié)果一致吻合，試驗(yàn)系統(tǒng)的可靠性和理論模型的正確性得到驗(yàn)證;

2)帶PMD貯箱的晃動(dòng)頻率受貯箱內(nèi)部結(jié)構(gòu)影響明顯，隨充液比總體頻率呈非單調(diào)變化趨勢(shì)，其覆蓋范圍達(dá)到0.5～1.1 Hz;

3)由于受到防晃錐的影響，阻尼呈現(xiàn)非單調(diào)變化趨勢(shì);

4)帶PMD貯箱的晃動(dòng)質(zhì)量隨充液比的增加先增大后減小，在0.76充液比附近時(shí)晃動(dòng)質(zhì)量最大;

5)帶PMD貯箱的液體晃動(dòng)擺點(diǎn)高度隨充液比總體呈減少趨勢(shì)，由于擺點(diǎn)高度與一階晃動(dòng)頻率有關(guān)，在0.45充液比附近擺點(diǎn)高度值較大.

3 結(jié)論

通過(guò)對(duì)大容量推進(jìn)劑貯箱液體晃動(dòng)性能的試驗(yàn)研究工作，得到以下結(jié)論:

1)空殼貯箱的晃動(dòng)試驗(yàn)驗(yàn)證了等效動(dòng)力學(xué)模型理論的正確性及晃動(dòng)試驗(yàn)結(jié)果的可靠性與準(zhǔn)確性;

2)大容量推進(jìn)劑貯箱內(nèi)部PMD上的中間隔板及防晃錐對(duì)液體晃動(dòng)特性影響明顯，為后續(xù)貯箱PMD結(jié)構(gòu)設(shè)計(jì)提供參考依據(jù);

3)通過(guò)晃動(dòng)試驗(yàn)獲取的大量試驗(yàn)數(shù)據(jù)及晃動(dòng)參數(shù)能夠較為準(zhǔn)確地體現(xiàn)出大容量推進(jìn)劑貯箱內(nèi)液體晃動(dòng)特性，同時(shí)該試驗(yàn)研究結(jié)果為運(yùn)載火箭和航天器的總體及姿態(tài)軌道控制系統(tǒng)的設(shè)計(jì)優(yōu)化提供重要參考和數(shù)據(jù)支撐.

［1］曲廣吉.衛(wèi)星動(dòng)力學(xué)工程［M］.北京:中國(guó)科學(xué)技術(shù)出版社，2000.

［2］楊旦旦，岳寶增.低重環(huán)境下旋轉(zhuǎn)軸對(duì)稱貯箱內(nèi)液體晃動(dòng)研究［J］.宇航學(xué)報(bào)，2013，34(7):917-925.YANG D D，YUE B Z.Research on sloshing in axsymmetrical containers under low gravity［J］.Journal of Astronautics，2013，34(7):917-925.

［3］HYUNG R S.Unstructured grid based Reynolds-averaged navier-stokes method for liquid tank sloshing［J］.Journal of Fluid Engineering，Transaction of the ASME，2005，127(3):572-582.

［4］FALTINSEN O M，ROGNEBAKLE O F，TIMOKHA A N.Classification of three-dimensional nonlinear sloshing in a square-base tank with finite depth［J］.Journal of Fluids and Structures，2005(20):81-103.

［5］LUPPES R，HELDER J A，VELDMAN A E P.The numerical simulation of liquid sloshing in microgravity［C］//European Conference on Computational Fluid Dynamics.TU Delft，The Netherlands，2006.

［6］ADLER S，WARSHAVSHY A，PERETZ A.Low-cost cold-gas reaction control system for sloshsat FLEVO small satellite［J］.Journal of Spacecraft and Rockets，2005，24(2):345-351.

［7］VELDMAN A E P，GERRITS J，LUPPES R，et al.The numerical simulation of liquid sloshing on board spacecraft［J］.Journal of Computational Physics，2007，224:82-99.

［8］修觀，李鑫.火箭三級(jí)無(wú)動(dòng)力飛行段晃動(dòng)的一種建模方法［J］.上海航天，2015，32(1):22-25.XIU G，LI X.A modeling method for propellant-sloshing of launch vehicle in third stage flight without power［J］.Aerospace Shanghai，2015，32(1):22-25.

［9］陸國(guó)平，李道奎，曾擎.大型平鋪貯箱衛(wèi)星的液體晃動(dòng)模型與晃動(dòng)力矩分析［J］.上海航天，2011，28(1): 18-22.LU G P，LI D K，ZENG Q.Liquid sloshing model and sloshing torque analysis of large satellite with tiled tank［J］.Aerospace Shanghai，2011，28(1):18-22.

［10］黃華，楊雷，張熇，等.航天器貯箱大幅液體晃動(dòng)三維質(zhì)心面等效模型研究［J］.宇航學(xué)報(bào)，2010，31(1): 55-59.HUANG H，YANG L，ZHANG H，et al.Research on 3D constraint surface model for large amplitude liquid sloshing on spacecraft tank［J］.Journal of Astronautics，2010，31(1):55-59.

［11］林宏，彭慧蓮，董鍇.推進(jìn)劑貯箱液體晃動(dòng)的仿真研究與驗(yàn)證［J］.強(qiáng)度與環(huán)境，2011，38(5):25-30.LING H，PENG H L，DONG K.Simulation and test verification for liquid sloshing of propellant tank［J］.Structure＆Environment Engineering，2011，38(5): 25-30.

［12］李松，高芳清，楊翊仁，等.水平圓柱腔內(nèi)液體晃動(dòng)特性試驗(yàn)研究［C］.第14屆全國(guó)反應(yīng)堆結(jié)構(gòu)力學(xué)會(huì)議.北京:中國(guó)力學(xué)學(xué)會(huì)，2006.

［13］ABRAMSON H N.The dynamics behavior of liquids in moving containers［R］.NASA SP-106，1996.

Experiment of Liquid Sloshing Performance in Bulky Propellant Tank

HU Qi，LI Yong，YAO Can，LIU Jintao
(1.Beijing Institute of Control Engineering，Beijing 100190，China; 2.Beijing Engineering Research Center of Propulsion Technology，Beijing 100190，China)

Aiming at the bulky propellant tank with propellant management device(PMD)，experiment research on liquid sloshing performance is carried out.The liquid sloshing tests in a tank without PMD are developed under different liquid filling ratios.Based on the equivalent mechanical model，the theoretical result of liquid sloshing is calculated.The results of tests and calculation are analyzed and compared.The results are consistent，both the reliability of test system and the validity of theoretical model are well verified.Then the liquid sloshing tests in the tank with PMD are developed under different liquid filling ratios，and the test results indicate that the structure of inner PMD significantly affects the performance of liquid sloshing in the tank.All the experiment research results provide an important reference and data support for designing and optimizing the attitude and orbit control system of rocket and spacecraft.

liquid sloshing;slosh frequency;damping bulky propellant tank;PMD

V43

A

1674-1579(2016)03-0044-05

10.3969/j.issn.1674-1579.2016.03.008

胡 齊(1985—)，男，工程師，研究方向?yàn)楹教焱七M(jìn)技術(shù)與微重力下流體機(jī)理;李 永(1977—)，男，研究員，研究方向?yàn)楹教焱七M(jìn)技術(shù)、微重力下流體機(jī)理與先進(jìn)流動(dòng)測(cè)量技術(shù);姚 燦(1983—)，男，高級(jí)工程師，研究方向?yàn)楹教焱七M(jìn)技術(shù);劉錦濤(1986—)，男，工程師，研究方向?yàn)楹教炱魍七M(jìn)技術(shù)與微重力流動(dòng)理論.

*國(guó)家自然科學(xué)基金資助項(xiàng)目(51406010).

2016-01-26