推進系統參數變化對液體火箭縱向振動響應的影響分析

唐 冶，方 勃，張業偉，譚立軍，于子文

（1．哈爾濱工業大學航天學院，哈爾濱 150001；2．上海大學上海市應用數學和力學研究所，上海 200072；3．沈陽航空航天大學航空航天工程學部，沈陽 110136；4．北京宇航系統工程研究所，北京 100076）

推進系統參數變化對液體火箭縱向振動響應的影響分析

唐 冶1，方 勃1，張業偉2，3，譚立軍1，于子文4

（1．哈爾濱工業大學航天學院，哈爾濱 150001；2．上海大學上海市應用數學和力學研究所，上海 200072；3．沈陽航空航天大學航空航天工程學部，沈陽 110136；4．北京宇航系統工程研究所，北京 100076）

為了優化液體火箭推進系統，避免POGO振動引起液體火箭低頻振動環境惡化對火箭飛行過程的不利影響，建立了液體火箭POGO振動系統的動力學模型，利用動態靈敏度技術，提出了液體火箭縱向振動響應對推進系統參數的靈敏度時域分析模型。通過數值仿真，得到了推進系統流體慣性、阻力和剛度參數以及泵的動態增益變化對液體火箭縱向振動響應的影響規律。研究結果表明，液體火箭縱向振動響應對流體慣性和阻力參數的敏感程度比流體剛度參數明顯大，泵動態增益的變化對液體火箭縱向振動響應的影響最大，泵前短管的流體剛度變化對液體火箭縱向振動響應的影響最小。為減小液體火箭縱向振動，進一步研究POGO振動特性提供參考。

液體火箭；POGO振動；縱向振動響應；推進系統參數；靈敏度

液體火箭系統非常復雜，在發射過程中伴隨著眾多的動力學現象，其中POGO振動是火箭結構系統典型的流固耦合振動之一。這種現象是指液體火箭推進系統的壓力脈動、流量脈動與箭體結構系統發生耦合而產生的自激振動［1］。國內外各種火箭發射經驗表明［2－4］，這種振動不僅給火箭帶來不利的負載，而且危害宇航員的安全。因此學者們研究了很多抑制方法來降低或者消除POGO振動。文獻［5－6］主要采用在推進系統適當位置上安蓄壓器裝置的被動抑制方式來改變推進系統的脈動頻率，使其遠離箭體結構系統的縱向振動頻率，避免發生這種不穩定的POGO振動。Marius［7］通過改變管路系統的壓力和流量，研究了液體火箭POGO振動的主動控制技術。目前，已經開展過大量的POGO的穩定性分析［8－9］和時域仿真［10－11］，但是這些工作不具備優化處理能力，需要進一步的分析。

結構的動態靈敏度技術在工程中用于研究系統振動特性對參數的敏感程度，估計參數變化的效果，指出系統隨參數變化的改變量，確定振動系統的優化搜索方向，達到系統優化的目的［12］。黃益民［13］分析了充液管道模態的參數靈敏度，對充液管道的共振可靠性進行了研究。文獻［14］利用動態靈敏度技術，提出了振動傳遞路徑系統對路徑參數的靈敏度方法，分析了振動接受結構系統響應對各種非粘性阻尼的靈敏度。

在液體火箭POGO振動系統參數處于穩定性的條件下，建立了POGO振動系統的動力學模型，利用動態靈敏度技術，提出了液體火箭縱向振動響應對推進系統參數的靈敏度時域分析模型。通過數值仿真分析，得到了推進系統參數變化對液體火箭縱向振動響應的影響規律。為進一步研究POGO振動特性提供參考。

1 液體火箭POGO振動系統的動力學模型

1.1 液體火箭推進系統的動力學模型

液體火箭推進系統主要由貯箱、不可壓縮管路、可壓縮管路、多通接頭、波紋管、蓄壓器、泵和燃燒室等單元部件組成。根據文獻［15］和文獻［16］，利用牛頓力學原理，可建立各單元的動力學方程組。

貯箱單元的動力學方程為

式中：βnj為第n階箭體結構振型在單位加速度下導致的貯箱出口處的模態壓力，·q·為第n階結構縱向振動

n系統的模態加速度。

不可壓縮管路單元動力學方程組為

式中：ρ為液體火箭推進系統中推進劑密度參數。

可壓縮管路單元動力學方程組為

式中：H2i和H2j分別為可壓縮管路單元入口和出口處的壓頭矢量。

多通接頭單元的動力學方程組為

式中：U為多通單元的出口管路數。波紋管單元的動力學方程組為

式中：g為重力加速度，r4i和r4j分別為波紋管單元入口和出口處的管壁振動的位移矢量。

蓄壓器單元的動力學方程組為

式中：α為泵的質量流增益因子，(m+)1為泵的動態增益。

燃燒室單元的動力學方程組為

式中：τ為燃燒時滯。

火箭發射過程中，液體推進劑流入推進系統的各個單元依次是貯箱、主管路、支管路，補償管路、蓄壓器單元、泵前短管、泵、泵后管路、燃燒室（見圖1），第一段管路為主管路，第二段管路為支管路，第三段管路為補償管路，第五段管路為泵前短管。

圖1 液體火箭推進系統的構造圖Fig．1 The structuralmap of propulsion system in liquid rockets

由于泵前各管路的容積特性較大，而泵后管路的容積特性較小，所以建立液體火箭推進系統的動力學模型時，泵前各管路作為可壓縮管路，而泵后管路作為不可壓縮管路。根據各管路單元、蓄壓器單元與泵單元的動力學方程，可列出這些單元入口與出口脈動壓力用脈動流量及其導數表示的動力學方程，以相鄰單元的共同節點脈動壓力相等為耦合條件，消去脈動壓力，可依序列出各個節點脈動流量為狀態變量的每兩個相鄰單元的耦合動力學方程。以推進系統中的所有節點脈動流量為狀態變量，將各個節點的動力學方程組合，進而得到液體火箭推進系統的動力學模型為

式中：MT、CT和KT分別為液體火箭推進系統的慣性矩陣，阻尼矩陣和剛度矩陣。FT為液體火箭推進系統所受到的外力列陣。x、x·和x··分別為液體火箭推進系統的狀態位移變量列陣、狀態速度變量列陣和狀態加速度變量列陣。

1.2 液體火箭結構系統的動力學模型

根據文獻［15］，采用振動模態理論，則液體火箭結構縱向振動系統的動力學模型為：

式中：Mn、ξn和ωn分別為第n階箭體結構縱向振動系統的模態質量、模態阻尼比和模態頻率。qn、q·n和Qn分別為第n階結構縱向系統的模態位移、模態速度和所受到的廣義力。φnk為結構縱向系統的第n階在第k個推進劑單元上的振型，Fk為第k個推進劑單元對結構縱向系統的作用力，由于多通接頭單元和波紋管單元對結構的作用非常小，因此根據文獻［15］。對結構縱向系統的作用力的推進劑單元主要是：

貯箱單元對第n階箭體結構振型的作用力為：

式中：φnj為貯箱出口處的第n階模態位移。

不可壓縮管路單元對結構的作用力為：

可壓縮管路單元對結構的作用力為：

式中：N2z為可壓縮管路單元中點處流動方向上的單位矢量。

蓄壓器單元對結構的作用力為：

泵單元對結構的作用力為

燃燒室單元對結構的作用力為：

式中：N7為燃燒室處的穩態流動方向的單位矢量。Cf為推力系數。

液體火箭縱向振動系統的位移可表示為各階模態位移之和，根據振動系統的模態截斷理論以及液體火箭POGO振動的低頻振動的特點，因此取箭體結構縱向系統的前兩階模態位移進行分析。

1.3 液體火箭POGO振動系統的動力學模型

根據液體火箭推進系統的動力學模型和箭體結構縱向振動系統的動力學模型以及它們之間的耦合關系，可建立液體火箭POGO振動系統的動力學模型為：

2 液體火箭縱向振動響應靈敏度時域模型

基于液體火箭POGO振動系統的動力學模型，利用Kronecker代數和矩陣微分理論，可導出液體火箭POGO振動系統響應對推進系統參數的靈敏度方程組為：

設δ為V中的任意一個變量（液體火箭推進系統的某個參數），根據式（17）和式（18），可導出液體火箭POGO振動系統響應對參數δ的靈敏度時域模型為：

3 數值結果與分析

由于蓄壓器的流體慣性和阻力相對管路小的多，并且蓄壓器對推進系統的作用是通過柔度表現的，所以數值分析中不考慮蓄壓器的流體慣性和阻力。采用Nemark積分方法對式（19）進行數值計算。

3.1 流體慣性參數變化對響應的影響

以液體火箭推進系統的流體慣性參數作為分析參數，其他參數固定不變，得到不同時刻液體火箭縱向振動響應對各個流體慣性參數的靈敏度（見圖2）。

圖2 液體火箭縱向振動響應對流體慣性參數的靈敏度隨時間變化圖Fig．2 The sensitivity of longitudinal vibration response in liquid rockets to the inertial parameters of fluid changes with time

從圖2可知，液體火箭縱向振動響應對主管路流體慣性、支管路流體慣性、補償管路流體慣性、泵前短管流體慣性、泵流體慣性和泵后管路流體慣性的靈敏度隨時間的變化。在流體慣性參數中，主管路流體慣性的變化對液體火箭縱向振動響應的影響最大，泵流體慣性的變化對液體火箭縱向振動響應的影響最小。

3.2 流體阻力參數變化對響應的影響

以液體火箭推進系統的流體阻力參數作為分析參數，其他參數固定不變，得到不同時刻液體火箭縱向振動響應對各個流體阻力參數的靈敏度（見圖3）。

圖3 液體火箭縱向振動響應對流體阻力參數的靈敏度隨時間變化圖Fig．3 The sensitivity of longitudinal vibration response in liquid rockets to the drag parameters of fluid changeswith time

從圖3可知，液體火箭縱向振動響應對主管路流體阻力、支管路流體阻力、補償管路流體阻力、泵前短管流體阻力、泵流體阻力和泵后管路流體阻力的靈敏度隨時間的變化。在流體阻力參數中，主管路流體阻力的變化對液體火箭縱向振動響應的影響最敏感，泵流體阻力的變化對液體火箭縱向振動響應的影響最不敏感。

在實際液體火箭發動機系統中，泵與貯箱都是固定在箭體結構上，而在它們之間的管路具有相當大的柔性。所以當火箭發射過程中，火箭整體處于振動狀態時，管路在同類參數變化中比泵的參數對推進系統的流量脈動影響大的多，根據POGO振動系統的機理可知，推進系統的流量脈動的改變能進一步的影響液體火箭縱向振動。在本文的POGO模型中，主管路是所有管路中流體慣性和阻力最大的管路，所以在同類管路參數中，它的變化對液體火箭縱向振動響應的影響應該是最大的，而泵應該是最小的。這與本文在流體慣性、阻力參數變化對液體火箭縱向振動響應影響分析中所得到的結果是一致。

3.3 流體剛度參數變化對響應的影響

以液體火箭推進系統的流體剛度參數作為分析參數，其他參數固定不變，得到不同時刻液體火箭縱向振動響應對各個流體剛度參數的靈敏度（見圖4）。

圖4 液體火箭縱向振動響應對流體剛度參數的靈敏度隨時間變化圖Fig．4 The sensitivity of longitudinal vibration response in liquid rockets to the stiffness parameters of fluid changes with time

從圖4可知，液體火箭縱向振動響應對主管路流體剛度、支管路流體剛度、補償管路流體剛度、蓄壓器流體剛度、泵前短管流體剛度和泵流體剛度的靈敏度隨時間的變化。在流體剛度參數中，蓄壓器流體剛度變化對液體火箭縱向振動響應的影響最強烈，泵前短管流體剛度變化對液體火箭縱向振動響應的影響最弱小。

在液體火箭推進系統中引入蓄壓器主要是通過它的柔度特性來吸收推進系統中的脈動能量，因為蓄壓器的柔度特性要比管路大的多，所以蓄壓器的流體剛度變化比管路對推進系統流量脈動的影響要大，從而在流體剛度參數中它的改變對液體火箭縱向振動的影響應該是最敏感的。這與流體剛度參數對液體火箭縱向振動的影響分析所得到的結果相一致。

3.4 泵的動態增益變化對響應的影響

以液體火箭推進系統的泵的動態增益作為分析參數，其他參數固定不變，得到不同時刻液體火箭縱向振動響應對泵的動態增益的靈敏度（見圖5）。

從圖5可知，液體火箭縱向振動響應對泵動態增益的靈敏度隨時間的變化。

圖5 液體火箭縱向振動響應對泵的動態增益的靈敏度隨時間變化圖Fig．5 The sensitivity of longitudinal vibration response in liquid rockets to dynamic gain of pump changeswith time

推進系統中泵的動態增益在液體火箭發動機系統中起著增加管路中流體壓力的作用，因此它的改變會直接導致火箭發動機燃燒室處的脈動壓力的變化，而此脈動壓力是直接作用于液體火箭縱向振動系統上，所以泵的動態增益改變應該是所有參數中對液體火箭縱向振動響應的影響最大。

比較圖2～圖5，泵動態增益的變化對液體火箭縱向振動響應的影響最敏感，泵前短管的流體剛度變化對液體火箭縱向振動響應的影響最不敏感，另外，從圖2～圖5還可以看出，液體火箭縱向振動響應對流體慣性和阻力參數的敏感程度比流體剛度參數明顯大，這是符合系統振動的常識。所以在設計推進系統的動力學參數時，要充分依序考慮泵的動態增益、流體慣性參數、流體阻力參數、流體剛度參數的情況，以便能減小液體火箭的縱向振動響應，改善火箭發射的振動環境。

4 結 論

基于液體火箭POGO振動系統的動力學模型，提出了液體火箭縱向振動響應對推進系統參數的靈敏度時域分析模型，分析了推進系統參數變化對液體火箭縱向振動響應的影響，得到的主要結論為：

（1）液體火箭縱向振動響應對流體慣性和阻力參數的敏感程度比流體剛度參數明顯大。

（2）泵動態增益的變化對液體火箭縱向振動響應的影響最大，泵前短管的流體剛度變化對液體火箭縱向振動響應的影響最小。

［1］Rubin S．Longitudinal instability of liquid rockets due to propulsion feedback（POGO）［J］．Journal of Spacecraft and Rockets，1966，3（8）：1188－1195．

［2］Nagai H，Noda K，Yamazaki I，et al．Status of H-Ⅱrocket first stage propulsion system［J］．Journal of Propulsion and Power，1992，8（2）：313－319．

［3］Meyers J F．DeltaⅡ-A new generation begins［J］．AIAA Journal，1989，89－2740：1－10．

［4］黃懷德．液體火箭的POGO振動研究［J］．振動工程學報，1987，1（1）：5－13．

HUANG Huai-de．Research into POGO vibration of liquid rockets［J］．Journal of Vibration Engineering，1987，1（1）：5－13．

［5］Ujino T．POGO Prevention of H-2 launch vehicle［J］．AIAA Journal，1994，94－1624－CP：2858－2867．

［6］Lack MH，Rubin S．Passive suppression of POGO on the space shuttle［R］．NASA CR－132452，1974．

［7］Marius T J．Investigation on active suppression of POGO［J］．AIAA Journal，1995，1995－3311－CP：1239－1249．

［8］Zhao Z H，Ren G X，Yu ZW，et al．Parameter study on POGO stability of liquid rockets［J］．Journal of Spacecraftand Rockets，2011，48（3）：537－541．

［9］Dotson K W，Rubin S，Sako B H．Mission-specific POGO stability analysis with correlated pump parameters［J］．Journal of Propulsion and Power，2005，21：619－626．

［10］Muller S，Breviere F，Kemilis A，et al．Influence pump cavitation process on POGO diagnosis for the A5E／CA upper stage［J］．AIAA Journal，2010，2010－6892：1－18．

［11］Ryan R S，Kiefling A K，Buchanan H J，etal．Simulation of Saturn V S－Ⅱstage propellant feedline dynamics［J］．Journal of Spacecraft and Rockets，1970，7（12）：1407－1412．

［12］Choi K M，Jo H K，Kim W H，et al．Sensitivity analysis of non-conservative eigensystems［J］．Journal of Sound and Vibration，2004，271：997－1011．

［13］黃益民，劉偉，劉永壽，等．充液管道模態的參數靈敏度及其共振可靠性分析［J］．振動與沖擊，2010，29（1）：193－196．

HUANG Yi-min，LIU Wei，LIU Yong-shou，et al．Parameter sensitivity and resonance reliability of a fluid-filled pipeline［J］．Journal of Vibration and Shock，2010，29（1）：193－196．

［14］Zhang Yi-min，Huang Xian-zhen，Zhao Qun．Sensitivity analysis for vibration transfer path systems with non-viscous damping［J］．Journal of Vibration and Control，2011，17（7）：1042－1048．

［15］Oppenheinm B W，Rubin S．Advanced POGO stability analysis for liquid rockets［J］．Journal of Spacecraft and Rockets，1993．30（3）：360－373．

［16］王其政，黃懷德，姚德源．結構耦合動力學［M］．北京：宇航出版社，1999．

Influences of parameters variation of propulsion system on longitudinal vibration response of liquid rockets

TANG Ye1，FANG Bo1，ZHANG Ye-wei2，3，TAN Li-jun1，YU Zi-wen4
（1．School of Astronautics，Harbin Institute of Technology，Harbin 150001，China；2．Shanghai Institute of Applied Mathematics and Mechanics，Shanghai200072，China；3．Department of Astronautics，Shenyang Aerospace University，Shenyang 110136，China；4．Beijing Institute of Astronautical Systems Engineering，Beijing 100076，China）

To optimize the propulsion system of liquid rockets and avoid the low frequency vibration environment deterioration of liquid rockets caused by POGO vibration，a dynamicmodel of POGO vibration system in liquid rocketswas established，and a time domain model for analysing the sensitivity of longitudinal vibration response of liquid rockets（LVRLR）with respect to parameters of propulsion system was proposed based on the dynamic sensitivity technique．Numerical calculations，show that the variations of inertial parameters，drag parameters and stiffness parameters of the fluid and the variation of dynamic gain of the pump for propulsion system have important effects on LVRLR and their influenceswere analysed．The results show the sensitive extents of LVRLR with respect to the inertial parameters and the drag parameters of fluid are sharply greater than the stiffness parameters of fluid，the influence of the dynamic gain variation of pump on LVRLR is themaximum，and the influence of the variation of fluid stiffness in the short pipe at pump front on LVRLR is theminimum．The results provide references to the reduction of longitudinal vibration of liquid rockets and to the further study of the characteristics of POGO vibration．

liquid rockets；POGO vibration；longitudinal vibration response；parameters of propulsion system；sensitivity

V434；O322

A

10．13465／j．cnki．jvs．2015．12．032

國家自然科學基金（11202140）

2014－03－17 修改稿收到日期：2014－06－06

唐冶男，博士生，1984年3月生

方勃男，教授，博士生導師，1964年8月生