液體火箭注氣式蓄壓器工作過程仿真及控制規(guī)律研究

張青松，范瑞祥，張 兵，胡久輝

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

0 引 言

液體運(yùn)載火箭通常采用蓄壓器來抑制POGO振動(dòng)，通過設(shè)計(jì)合理的蓄壓器性能參數(shù)，將動(dòng)力系統(tǒng)的頻率調(diào)整到合適的范圍內(nèi)，以降低或消除動(dòng)力系統(tǒng)推力脈動(dòng)與結(jié)構(gòu)縱向振動(dòng)之間耦合，進(jìn)而達(dá)到抑制火箭縱向耦合振動(dòng)的目的[1,2]。根據(jù)蓄壓器的結(jié)構(gòu)特點(diǎn)，可以將蓄壓器分為貯氣式和注氣式兩大類。

目前國外低溫液體火箭一般采用注氣式蓄壓器來抑制POGO振動(dòng)，中國液體運(yùn)載火箭目前均采用金屬膜盒式蓄壓器抑制POGO振動(dòng)[3]。注氣式蓄壓器中的工作氣體與推進(jìn)劑直接接觸，沒有彈性結(jié)構(gòu)隔離，更適應(yīng)于低溫工作環(huán)境和大能量值需求。注氣式蓄壓器在中國運(yùn)載火箭上還沒有應(yīng)用的先例，因此對此類蓄壓器的工作特性研究還很少。文獻(xiàn)[4]對注氣式蓄壓器的工作特性及其對推進(jìn)系統(tǒng)的影響開展了仿真分析，由于在分析時(shí)設(shè)定了“恒溫、恒流量、恒出口壓力”等理想化假設(shè)條件，難以反映實(shí)際的動(dòng)態(tài)過程。文獻(xiàn)[5]介紹了Ares I火箭二級注氣式蓄壓器的方案論證情況，給出了用于預(yù)估蓄壓器單次充氣所需氦氣用量的計(jì)算模型，并對比分析了不同方案對蓄壓器慣性系數(shù)的影響。文獻(xiàn)[6]和[7]研究了宇宙神5火箭一級發(fā)動(dòng)機(jī)關(guān)機(jī)前注氣式蓄壓器排氣所引起的箭體結(jié)構(gòu)低頻振動(dòng)問題，分析了注氣式蓄壓器在發(fā)動(dòng)機(jī)關(guān)機(jī)前的排氣需求及其對結(jié)構(gòu)系統(tǒng)的影響。

本文以注氣式蓄壓器為研究對象，基于實(shí)際物理過程建立描述飛行過程中蓄壓器內(nèi)壓力、溫度以及氣腔容積變化規(guī)律的狀態(tài)變量模型。針對一個(gè)多余氣體外排的注氣式蓄壓器系統(tǒng)，開展不同入口壓力條件下的工作過程仿真分析和充放氣控制參數(shù)研究，目的是建立定容積注氣式蓄壓器的控制規(guī)律和參數(shù)優(yōu)化設(shè)計(jì)方法。

1 注氣式蓄壓器的系統(tǒng)方案

液體運(yùn)載火箭POGO抑制用注氣式蓄壓器，自帶充氣和放氣系統(tǒng)，蓄壓器內(nèi)的氣體與推進(jìn)劑直接接觸，在火箭飛行過程中不斷有氣體注入蓄壓器的氣腔中，通常在蓄壓器內(nèi)設(shè)置限位管來控制蓄壓器的氣腔容積。根據(jù)蓄壓器內(nèi)多余氣體排放方式的不同，可以將注氣式蓄壓器進(jìn)一步分為兩大類：一類是多余氣體排入動(dòng)力系統(tǒng)，另一類是氣體排出箭體外部，對動(dòng)力系統(tǒng)的正常工作不產(chǎn)生任何影響。

圖1是一個(gè)多余氣體外排式的注氣式蓄壓器系統(tǒng)方案，該蓄壓器由環(huán)形氣腔、環(huán)形液腔、充氣系統(tǒng)和排放系統(tǒng)組成。氣液界面的位置由排放系統(tǒng)的泄出口位置確定，蓄壓器的環(huán)形液腔與液氧輸送管之間通過一圈連通孔相連，為降低蓄壓器的慣性系數(shù)和阻力系數(shù)，連通孔的面積應(yīng)該盡量大一些。蓄壓器的充氣系統(tǒng)主要由氣瓶、過濾器、充氣控制電磁閥、限流孔板、氣體擴(kuò)散器和充氣管路組成，火箭點(diǎn)火時(shí)打開充氣控制電磁閥，在火箭飛行過程中持續(xù)向蓄壓器內(nèi)充入氦氣。蓄壓器的排放系統(tǒng)由泄出控制元件、泄出控制電液閥和相應(yīng)管路組成，火箭點(diǎn)火前打開泄出電液閥，向火箭外部排放泄出口附近的氣體或液體，保證蓄壓器的氣腔容積在設(shè)計(jì)指標(biāo)要求的范圍之內(nèi)，提供全箭POGO抑制所需的蓄壓器工作容積。

圖1 注氣式蓄壓器系統(tǒng)原理 Fig.1 Schematic Diagram of Gas Filled Accumulator System

2 注氣式蓄壓器的動(dòng)力學(xué)模型

注氣式蓄壓器在火箭飛行過程中不斷有氣體注入到蓄壓器的氣腔中，多余氣體通過溢出管路排入輸送管中或排出箭體外部，氣液界面在溢出口附近維持動(dòng)態(tài)平衡，蓄壓器充氣、放氣參數(shù)設(shè)計(jì)的重要目標(biāo)是盡量控制氣腔容積保持相對恒定。為保證注氣式蓄壓器正常工作，應(yīng)使蓄壓器的充氣流量、泄出能力、蓄壓器安裝位置處的壓力變化規(guī)律三者相適應(yīng)，確保蓄壓器的工作容積在火箭飛行過程中滿足設(shè)計(jì)指標(biāo)要求。

注氣式蓄壓器在工作過程中不斷有氣體進(jìn)入和流出，系統(tǒng)的連續(xù)方程和能量守恒方程如式（1）、式（2）所示[8]：

式中um、U分別為注氣式蓄壓器氣腔內(nèi)的氣體質(zhì)量、內(nèi)能；im為進(jìn)出蓄壓器氣腔的氣體流量；ih為進(jìn)出蓄壓器的氣體比焓；Q為蓄壓器內(nèi)氣體與外界的換熱；W為蓄壓器內(nèi)氣體對外所作的功。

上世紀(jì)60年代起，氯化鐵逐漸成為電路板及印刷業(yè)中廣泛應(yīng)用的蝕刻劑。90年代開始用于銅版畫制版。與氯化銅比較，它不產(chǎn)生有害氣體，安全性更高。一般的銅版畫工作室條件下，基于氯化鐵的蝕刻劑毒性低，更為安全，也更為便捷。

在式（1）、式（2）的基礎(chǔ)上，結(jié)合低溫下實(shí)際氣體的狀態(tài)方程、內(nèi)能和焓的一般關(guān)系式，建立如下模型描述注氣式蓄壓器工作過程中的壓力和溫度變化：

式中uρ、uP、uT、uV、uh分別為注氣式蓄壓器氣腔內(nèi)的氣體密度、壓力、溫度、容積、比焓；glQ·、gwQ·分別為蓄壓器內(nèi)氣體與液面之間、氣體與結(jié)構(gòu)壁面之間的換熱熱流。公式中的偏導(dǎo)數(shù)系數(shù)項(xiàng)根據(jù)實(shí)際氣體的狀態(tài)方程以及比熱容確定。

蓄壓器的氣腔容積變化與進(jìn)出蓄壓器的推進(jìn)劑流量相關(guān)，受蓄壓器氣腔壓力和蓄壓器安裝位置處輸送管路中的壓力變化的影響。在臨界工況下，蓄壓器的排放流量按如下公式進(jìn)行計(jì)算：

式中mlq為蓄壓器的液體排放流量；lC為蓄壓器泄出孔排液狀態(tài)的流量系數(shù)；lρ為推進(jìn)劑的密度；sP為推進(jìn)劑的飽和蒸汽壓；mgq為蓄壓器的氣體排放流量；gC為蓄壓器泄出孔排氣狀態(tài)的流量系數(shù)；A為泄出孔面積；gA為氣體排出部分的面積，是蓄壓器氣枕容積的函數(shù)，與泄出控制元件的具體結(jié)構(gòu)形式相關(guān)；R、k分別為氣體常數(shù)和絕熱指數(shù)。

計(jì)算蓄壓器內(nèi)氣體與液面之間的換熱以及氣體與結(jié)構(gòu)壁面之間的換熱可以參考如下方程[9]：

式中g(shù)lh為氣液之間的換熱系數(shù)；glA為氣液換熱面積；LsT為蓄壓器內(nèi)推進(jìn)劑的飽和溫度；k為導(dǎo)熱系數(shù)；L為特征長度；Gr為格拉曉夫數(shù)；Pr為普朗特?cái)?shù)。

式中g(shù)wh為氣枕與結(jié)構(gòu)壁面之間的換熱系數(shù)；gwA為氣枕與結(jié)構(gòu)壁面之間的換熱面積；wT為與氣枕接觸的殼體結(jié)構(gòu)溫度；μ為粘度系數(shù)；xN為軸向過載系數(shù)；β為容積膨脹系數(shù)；pC為定壓比熱；m·為進(jìn)入蓄壓器的氣體流量；dA為氣體擴(kuò)散器出口面積；ξ為與消能器結(jié)構(gòu)型式有關(guān)的修正因子。

3 注氣式蓄壓器工作過程仿真及充放氣控制研究

3.1 注氣式蓄壓器工作過程仿真

對圖1所示的注氣式蓄壓器系統(tǒng)，由存放在液氧貯箱內(nèi)的單獨(dú)冷氦氣瓶提供氣源，蓄壓器的充氣、放氣過程通過電磁閥、電液閥實(shí)施控制，為提高系統(tǒng)可靠性，盡量簡化系統(tǒng)方案，在火箭飛行過程中在能滿足設(shè)計(jì)指標(biāo)要求的前提下不再對電磁閥、電液閥進(jìn)行調(diào)節(jié)控制。針對一個(gè)設(shè)計(jì)氣枕容積為20 L、容積控制偏差在±1 L之內(nèi)的注氣式蓄壓器方案，本文對液氧輸送系統(tǒng)、注氣式蓄壓器系統(tǒng)開展工作過程的仿真分析。 圖2是一個(gè)給定的輸送管出口壓力變化規(guī)律所對應(yīng)的蓄壓器氣腔容積仿真結(jié)果，從圖2可以看出，在大部分時(shí)間內(nèi)蓄壓器容積均能控制在設(shè)計(jì)要求的范圍之內(nèi)，對系統(tǒng)的工作狀態(tài)和外部環(huán)境變化有較強(qiáng)的適應(yīng)能力，初步驗(yàn)證了關(guān)鍵設(shè)計(jì)參數(shù)的正確性。在飛行后期助推器分離導(dǎo)致芯一級液氧輸送管出口壓力突然下降，由此引起注氣式蓄壓器的氣腔容積在此分離時(shí)刻突然升高至29 L附近，超出了蓄壓器的容積控制范圍，隨后蓄壓器內(nèi)的氣體通過泄出孔排出，蓄壓器容積降至設(shè)計(jì)指標(biāo)范圍內(nèi)。由于輸送管路中壓力突然下降導(dǎo)致蓄壓器氣腔容積增大、氣液界面高度降低，蓄壓器內(nèi)的氣體有可能通過連通孔進(jìn)入輸送管路中，影響發(fā)動(dòng)機(jī)的正常工作，所以蓄壓器液腔與輸送管之間的連通孔位置設(shè)計(jì)需要充分考慮此類工況，應(yīng)在連通孔以上預(yù)留出足夠的液腔容積，以防止蓄壓器內(nèi)的氣體進(jìn)入輸送管中。

圖2 飛行過程中蓄壓器氣腔容積及輸送管壓力變化 Fig.2 Variation of Accumulator Gas Volume and Pressure in Pipeline Through Flight

圖3 為飛行過程中蓄壓器排出氣體流量及排出推進(jìn)劑流量的仿真結(jié)果曲線，隨著飛行過程中氣瓶壓力的下降，注入蓄壓器中的氣體流量降低，蓄壓器的排氣流量也逐漸下降。隨著輸送管出口壓力的升高，溢出口中液體排放部分面積增大，推進(jìn)劑的排出流量逐漸增加，蓄壓器在飛行全程的液氧排放總量約為53 kg。

圖3 飛行過程中蓄壓器的排氣流量和排液流量變化 Fig.3 Variation of Helium and Liquid Oxygen Mass Flow Rate Released from Accumulator Through Flight

3.2 充放氣控制規(guī)律研究

3.2.1 不同充氣孔板面積的影響

在泄出路等效孔板直徑確定的情況下，對充氣路不同孔板面積下注氣式蓄壓器系統(tǒng)的工作過程進(jìn)行仿真計(jì)算。圖4、圖5分別給出了蓄壓器氣腔容積、氦氣排放流量的仿真計(jì)算結(jié)果。當(dāng)充氣孔板面積較小時(shí)，蓄壓器在飛行過程中維持在較低容積狀態(tài)，氦氣排放流量較低，液氧排放流量較大。增大充氣孔板面積，蓄壓器氣腔容積在控制范圍內(nèi)略有增大，同時(shí)氦氣排放流量增大、液氧排放流量降低。在較大的充氣孔板面積變化范圍內(nèi)，該注氣式蓄壓器均能在預(yù)期的氣腔容積控制范圍之內(nèi)正常工作，對供氣流量的變化表現(xiàn)出較強(qiáng)的適應(yīng)性。當(dāng)充氣孔板面積過大時(shí)，充氣流量將大于泄出路的最大排氣能力，蓄壓器在工作初期氣腔容積超出了設(shè)計(jì)指標(biāo)范圍，氣液界面低于溢出口位置較多。

圖4 不同充氣孔板面積下的蓄壓器氣腔容積 Fig.4 Variation of Accumulator Gas Volume with Different Diameter of Charging Orifice

圖5 不同充氣孔板面積下的氦氣排放流量 Fig.5 Variation of Released Helium Mass Flow Rate with Different Diameter of Charging Orifice

3.2.2 不同等效泄出孔直徑的影響

在充氣路孔板面積已經(jīng)確定的情況下，對泄出路不同等效孔板直徑下注氣式蓄壓器系統(tǒng)的工作過程進(jìn)行仿真計(jì)算。圖6、圖7分別給出了蓄壓器氣腔容積、氦氣排放流量的仿真計(jì)算結(jié)果。當(dāng)?shù)刃钩隹字睆捷^小時(shí)，泄出路的最大排氣能力小于充氣流量，蓄壓器在工作前期以及助推器分離后的一段時(shí)間內(nèi)氣腔容積超出了設(shè)計(jì)指標(biāo)范圍，氣液界面低于溢出口位置較多。增大泄出路等效孔板直徑，蓄壓器氣腔容積在控制范圍內(nèi)有所減小，同時(shí)在飛行過程中液氧的排放流量也將會增加，而氣體的排放流量將趨于一致；等效泄出孔直徑增大后，增強(qiáng)了蓄壓器的排氣能力，助推器分離后蓄壓器氣腔容積也會以更快的速度回到設(shè)計(jì)指標(biāo)范圍內(nèi)。仿真計(jì)算表明，在較大的等效泄出孔直徑變化范圍內(nèi)，該注氣式蓄壓器均能在預(yù)期的氣腔容積控制范圍之內(nèi)正常工作，對泄出流量的變化表現(xiàn)出了較強(qiáng)的適應(yīng)性。

圖6 不同等效泄出孔板直徑下的蓄壓器氣腔容積 Fig.6 Variation of Accumulator Gas Volume with Different Equivalent Diameter of Venting Orifice

圖7 不同等效泄出孔板直徑下的氦氣排放流量 Fig.7 Variation of Released Helium Mass Flow Rate with Different Equivalent Diameter of Venting Orifice

3.2.3 輸送管出口不同壓力曲線的影響

輸送管出口壓力受貯箱增壓控制以及火箭飛行過載的影響，對注氣式蓄壓器的容積變化以及充放氣控制規(guī)律影響較大。圖8給出了3種具有代表性的液氧輸送管出口壓力變化及蓄壓器氣腔容積變化規(guī)律的仿真結(jié)果：分別代表了在火箭飛行過程中蓄壓器入口壓力“瞬態(tài)快速下降”、“相對平穩(wěn)”、“逐漸升高”3種典型變化工況。用一組相同的充氣控制元件、放氣控制元件以及氣瓶系統(tǒng)，針對圖8所示的3種蓄壓器安裝位置處輸送管中的壓力變化，對本文所研究的注氣式蓄壓器開展仿真分析。

圖8 液氧輸送管不同出口壓力工況下蓄壓器氣腔容積變化 Fig.8 Simulation of Accumulator Gas Volume Variation with Different Pressure in Pipeline

3種不同的輸送管出口壓力工況下，在火箭飛行的大部分時(shí)間內(nèi)，蓄壓器氣腔容積都能在設(shè)定的控制范圍內(nèi)正常變化，滿足設(shè)計(jì)指標(biāo)要求。工況3由于輸送管出口壓力突然大幅下降，導(dǎo)致蓄壓器氣腔容積突然增大，短時(shí)超過了設(shè)計(jì)指標(biāo)所要求的范圍，隨后又降到了正常的容積控制范圍之內(nèi)。針對此類情況，注氣式蓄壓器設(shè)計(jì)時(shí)需要保證在蓄壓器氣腔容積快速增大的過程中，多余氣體不會進(jìn)入到輸送系統(tǒng)中，以避免影響發(fā)動(dòng)機(jī)的正常工作。

4 結(jié) 論

本文建立了注氣式蓄壓器的工作過程仿真模型，針對多余氣體外排的注氣式蓄壓器系統(tǒng)方案，開展了蓄壓器的工作過程仿真和充放氣控制規(guī)律研究。主要結(jié)論如下：

a）為保證注氣式蓄壓器能在設(shè)計(jì)指標(biāo)要求的范圍內(nèi)正常工作，蓄壓器的充氣控制元件和放氣控制元件應(yīng)該協(xié)調(diào)匹配：設(shè)計(jì)充氣流量需要小于各飛行工況下泄出路的最大放氣能力。

b）對于已經(jīng)確定的蓄壓器充氣流量而言，在確保泄出路最大放氣能力大于充氣流量的前提下，應(yīng)盡量減小泄出路的有效流通面積，以降低火箭飛行過程中蓄壓器的液氧消耗量。

c）針對輸送管中壓力存在突然下降的飛行工況，蓄壓器氣腔容積會突然增大，泄出路的最大放氣能力應(yīng)滿足輸送管中壓力突降時(shí)的蓄壓器泄壓要求；在氣瓶儲氣量充足的情況下，適當(dāng)增大注氣式蓄壓器充氣控制元件的流通面積，有助于降低火箭飛行過程中注氣式蓄壓器的液氧消耗量。