尹文娟，魏延明

(北京控制工程研究所，北京100190)

0 引言

雙組元統(tǒng)一推進氣路系統(tǒng)作為氣體增壓系統(tǒng)，其功能是通過對推進劑貯箱的增壓使推進劑輸送系統(tǒng)輸出相應(yīng)流量的推進劑為發(fā)動機提供氧化劑和燃料，因此在工作過程中其穩(wěn)定性必須符合一定要求才能保證發(fā)動機的正常工作.然而實際工作過程中，氣路系統(tǒng)總會受到外界和內(nèi)部一些因素的擾動，例如結(jié)構(gòu)參數(shù)和工作參數(shù)的變化、環(huán)境條件的改變等.如果氣路系統(tǒng)不穩(wěn)定，就會在任何微小的擾動作用下偏離原來的平衡狀態(tài)，并隨時間的推移而發(fā)散.目前，一般通過數(shù)學(xué)仿真和試驗的方法來分析氣路系統(tǒng)的動態(tài)特性［1-4］，而很少對氣路系統(tǒng)的穩(wěn)定性(抗干擾性)進行理論研究.傳統(tǒng)的穩(wěn)定性研究模型通常是基于對理想氣體狀態(tài)方程進行求導(dǎo)并采用等熵過程假設(shè)推導(dǎo)得到;文獻［5］針對某膜片式減壓閥的穩(wěn)定性進行了研究，但沒有對阻尼腔和卸荷腔建模;文獻［6-7］對單獨的減壓閥穩(wěn)定性進行了詳細研究，但并沒有涉及氣路系統(tǒng)各個組件間的耦合作用.

本文在雙組元推進氣路系統(tǒng)的數(shù)學(xué)模型以及減壓閥狀態(tài)空間模型的基礎(chǔ)上［6］進行拓展，采用線性化分析方法，分別建立了氣路系統(tǒng)組件的狀態(tài)空間模型，得到考慮各個組件間耦合的氣路系統(tǒng)以及減壓閥傳遞函數(shù)模型，探討了主要結(jié)構(gòu)參數(shù)對于氣路系統(tǒng)以及減壓閥穩(wěn)定性的影響規(guī)律.最后，提出了改善減壓閥和氣路系統(tǒng)穩(wěn)定性的主要措施.

1 氣路系統(tǒng)簡介

氣路系統(tǒng)的組成見圖1.雙組元推進系統(tǒng)氣路結(jié)構(gòu)復(fù)雜，包含減壓器、單向閥、氣瓶、貯箱、管路連接件等部件.增壓系統(tǒng)工作過程為:手動閥門打開后，氣瓶內(nèi)的高壓氣體經(jīng)過減壓閥減壓，之后分為兩路經(jīng)過單向閥，分別為氧化劑貯箱和燃料貯箱增壓.

圖1 氣路系統(tǒng)原理圖Fig.1 Gas system schematic

單向閥主要用于防止推進劑蒸汽逆向進入氣路系統(tǒng)而造成故障，其結(jié)構(gòu)原理見圖2.當(dāng)單向閥上、下游壓力差大于開啟壓差時閥芯打開，氣體從大、小閥芯與相應(yīng)閥座之間的縫隙及止回孔流向下游;當(dāng)上、下游壓力差小于開啟壓差時，閥芯運動件上所受的開啟力小于彈簧對閥芯的閉合力，閥芯關(guān)閉.

圖2 單向閥原理圖Fig.2 Check valve(CV)schematic

2 數(shù)學(xué)模型

文獻［6］給出了減壓閥的動態(tài)數(shù)學(xué)模型，單向閥的動態(tài)數(shù)學(xué)模型與減壓閥類似.同時，考慮到氣瓶輸出壓力為恒定值的工況，貯箱可以簡化成一個氣腔及氣體入口端口和一個液腔及液體出口端口組成的組合件，兩個容腔之間的隔膜看成是理想的幾何隔離面，兩側(cè)壓強相等.具體的數(shù)學(xué)公式不再贅述.

3 狀態(tài)空間模型

倘若考慮理想氣體，且不考慮環(huán)境溫度的影響，小擾動條件下，穩(wěn)態(tài)點附近的各個腔室溫度變化很小，因此這里忽略溫度變化帶來的影響.采用小偏差法對動態(tài)方程進行線性化處理，得到氣路系統(tǒng)在穩(wěn)定點附近的線性化狀態(tài)空間方程.

3.1 減壓閥狀態(tài)空間模型

文獻［6］給出考慮減壓閥下游流量為恒值的狀態(tài)空間模型，而在實際工作過程中，減壓閥下游流量和出口壓力輸出存在一種調(diào)節(jié)關(guān)系.考慮減壓閥下游流量對減壓閥出口壓力輸出的影響，假設(shè)減壓閥出口壓力等于單向閥入口壓力，下游流量為兩路單向閥入口流量之和，得到減壓閥狀態(tài)空間模型為:輸入變量:U1=［p1p7］，輸出變量:Y1=［p2］，狀態(tài)變量:X1=［hhdp2p3p4］T，

式中，系統(tǒng)矩陣A1，輸入與輸出矩陣B1C1分別為

式中，p1、p2、p3、p4分別為減壓閥高壓腔、低壓腔、反饋腔及卸荷腔壓力，p7為單向閥開啟腔壓力，T2、T3、T4分別為減壓閥低壓腔、反饋腔及卸荷腔溫度，V2、V3、V4分別為減壓閥低壓腔、反饋腔及卸荷腔容積，Qm1、Qm3、Qm4分別為減壓閥閥芯節(jié)流處、反饋孔處、卸荷孔處流量，Qm5為單向閥入口流入流量，R為氣體常數(shù)，h、hd、m分別為減壓閥閥芯運動組件位移、速度、質(zhì)量，Am為膜片有效面積，A2、A4、A1-A0分別為減壓閥低壓氣體、卸荷腔氣體以及高壓氣體作用在閥芯的有效面積，k、c分別為減壓閥彈性元件剛度，阻尼系數(shù)，dc3為反饋腔直徑，dc4為頂桿直徑，dc5為卸荷腔直徑.(*)0代表參數(shù)*在氣路系統(tǒng)穩(wěn)定工作狀態(tài)時的值.

3.2 傳遞函數(shù)模型

根據(jù)氣路系統(tǒng)各個組件的關(guān)系，得到氣路傳遞函數(shù)框圖如圖3.

圖3 氣路系統(tǒng)傳遞函數(shù)框圖Fig.3 Transfer function of gas system

圖中，P5、P8分別為單向閥出口壓力和貯箱氣體部分壓力，G1(s)為減壓閥的執(zhí)行元件傳遞函數(shù)，G2(s)為單向閥的執(zhí)行元件傳遞函數(shù)，G3(s)為貯箱的執(zhí)行元件傳遞函數(shù)，H2(s)為單向閥對于減壓閥出口的反饋回路傳遞函數(shù)，H3(s)為貯箱對于單向閥出口的反饋回路傳遞函數(shù).

根據(jù)圖3，可以分別得到從減壓閥入口到貯箱出口的氣路系統(tǒng)閉環(huán)傳遞函數(shù)以及帶有單向閥壓力、流量反饋回路的減壓閥閉環(huán)傳遞函數(shù).

4 動態(tài)響應(yīng)過程仿真與分析

考慮減壓閥入口壓力為恒定值的工況，對氣路系統(tǒng)從啟動到額定工況的工作過程進行動態(tài)響應(yīng)特性仿真.其中，工作介質(zhì)為氦氣，減壓閥入口壓力為12 MPa，系統(tǒng)溫度為 293.15 K.仿真中，在第 6 s，兩路貯箱下游閥門打開，在第12 s，下游閥門關(guān)閉.

從圖4的減壓閥、單向閥、貯箱氣體部分的壓力仿真曲線可以看出，在理想氣體絕熱流動的假設(shè)下，氣路系統(tǒng)的仿真曲線雖然在關(guān)閉階段和試驗曲線有所差異，沒有由傳熱引起的壓力升高再下降現(xiàn)象，但是對氣路系統(tǒng)在穩(wěn)態(tài)時刻的動態(tài)特性影響并不大.

圖4 氣路系統(tǒng)動態(tài)特性曲線Fig.4 Characteristics of gas system

由圖4所示的溫度變化可以看出，各腔室溫度雖然在下游閥門開啟之后一段時間里變化劇烈，但變化幅度并不大.因此，在穩(wěn)態(tài)工作過程，各個腔室的氣流狀態(tài)變化可近似為等溫過程.

5 動態(tài)穩(wěn)定性分析

一般，減壓閥的反饋孔直徑、彈簧阻尼系數(shù)、出口腔容積、閥芯運動組件質(zhì)量以及單向閥止回孔直徑對氣路穩(wěn)定性影響很大.因此，分別分析這幾個參數(shù)對氣路穩(wěn)定性的影響，同時，以主導(dǎo)閉環(huán)極點距離虛軸的距離作為穩(wěn)定性程度的判斷依據(jù).

5.1 減壓閥結(jié)構(gòu)參數(shù)對氣路穩(wěn)定性的影響

圖5為各個減壓閥相關(guān)結(jié)構(gòu)參數(shù)對于氣路穩(wěn)定性的影響規(guī)律示意圖.可以看出，使得氣路系統(tǒng)和減壓閥穩(wěn)定的減壓閥反饋孔直徑參數(shù)范圍是一致的.但在相同結(jié)構(gòu)參數(shù)下，氣路系統(tǒng)的穩(wěn)定性程度卻優(yōu)于減壓閥.事實上，在工作參數(shù)恒定，氣體流動均勻、絕熱的情況下，系統(tǒng)穩(wěn)定性范圍主要由閥門的閥芯運動規(guī)律決定.由于單向閥在穩(wěn)態(tài)工作點附近，大小閥芯的開度基本不會發(fā)生變化，因此，即使增加了單向閥和貯箱兩個組件，各個結(jié)構(gòu)參數(shù)對于氣路系統(tǒng)和減壓閥的穩(wěn)定性影響規(guī)律在范圍上是一致的.雖然單向閥和貯箱對氣路系統(tǒng)的穩(wěn)定性范圍基本不造成影響，卻使得氣路系統(tǒng)在各個結(jié)構(gòu)參數(shù)下的閉環(huán)主導(dǎo)極點距虛軸的距離存在一個極值，因此，在存在擾動的情況下，收斂速度大范圍慢于減壓閥.

改善減壓閥和氣路系統(tǒng)穩(wěn)定性主要有兩個途徑:適當(dāng)降低閥芯運動的靈敏度以及減緩反饋腔壓力對于閥芯運動的反饋.圖5(a)和圖5(b)分別表示不同反饋孔直徑和不同減壓閥下游管路容積對減壓閥以及氣路系統(tǒng)穩(wěn)定性的影響.可以看出，減小減壓閥反饋孔直徑和增大減壓閥下游管路容積都可以減緩減壓閥反饋腔壓力對于閥芯運動的反饋，提高減壓閥和氣路系統(tǒng)的穩(wěn)定性.同時，若減壓閥反饋孔直徑過小或者減壓閥下游管路容積過大則會使得減壓閥反饋腔的壓力變化延遲過大，穩(wěn)定性變差.如圖5(c)和圖5(d)所示，減小減壓閥閥芯運動組件質(zhì)量和增大減壓閥彈簧阻尼系數(shù)則是通過降低閥芯運動的靈敏度來提高減壓閥以及氣路系統(tǒng)的穩(wěn)定性.

5.2 單向閥結(jié)構(gòu)參數(shù)對氣路穩(wěn)定性的影響

圖6以減壓閥反饋孔直徑為例，研究了不同單向閥止回孔直徑對于減壓閥和氣路系統(tǒng)穩(wěn)定性的影響規(guī)律.

圖5 減壓閥結(jié)構(gòu)參數(shù)對穩(wěn)定性的影響Fig.5 Influences of PRV's structure parameters on the stability

由于不同的單向閥止回孔直徑不僅會影響下游的穩(wěn)態(tài)輸出，而且也會影響到其上游減壓閥的穩(wěn)態(tài)工作點，因此，從圖6看到，不同單向閥止回孔直徑對應(yīng)著不同的氣路穩(wěn)定性影響規(guī)律曲線.單向閥止回孔直徑越小，額定流量下的單向閥開啟腔和止回腔之前的壓差越大，施加給粘連在一起的大、小閥芯的力也就越大，單向閥閥芯開度越穩(wěn)定，氣路系統(tǒng)似乎越不容易產(chǎn)生振蕩.然而，從抗干擾性的角度來看，止回孔直徑的減小并不會影響使得減壓閥和氣路系統(tǒng)穩(wěn)定的結(jié)構(gòu)參數(shù)范圍.并且止回孔直徑越小，減壓閥和氣路系統(tǒng)的最大閉環(huán)主導(dǎo)極點距虛軸的距離反而大范圍減小，收斂速度變慢，穩(wěn)定性變差.

圖6 單向閥結(jié)構(gòu)參數(shù)對穩(wěn)定性的影響Fig.6 Influences of CV's structure parameters on the stability

6 結(jié)論

本文針對以往文獻對雙組元統(tǒng)一推進氣路系統(tǒng)抗干擾性理論研究方面的不足，在文獻［6］基礎(chǔ)上進行拓展，建立了考慮組件耦合的氣路系統(tǒng)閉環(huán)傳遞函數(shù)模型，探討了結(jié)構(gòu)參數(shù)對于氣路系統(tǒng)穩(wěn)定性的影響規(guī)律，從中得到以下結(jié)論:

(1)在工作參數(shù)恒定，氣體流動均勻、絕熱的情況下，系統(tǒng)的穩(wěn)定性范圍主要由閥門的閥芯運動規(guī)律決定.因此，即使增加了單向閥和貯箱兩個組件，各個結(jié)構(gòu)參數(shù)對于氣路系統(tǒng)和減壓閥的穩(wěn)定性影響規(guī)律在范圍上是一致的.

(2)雖然單向閥和貯箱對氣路系統(tǒng)的穩(wěn)定性范圍基本上不造成影響，卻使得氣路系統(tǒng)在各個結(jié)構(gòu)參數(shù)下的閉環(huán)主導(dǎo)極點距虛軸的距離存在極限值，使得其收斂速度大范圍慢于減壓閥.

(3)適當(dāng)減小反饋孔直徑、增大減壓閥下游管路容積、減小減壓閥閥芯運動組件質(zhì)量以及增大減壓閥彈簧阻尼系數(shù)都可以同時提高減壓閥和氣路系統(tǒng)的穩(wěn)定性.

(4)單向閥止回孔直徑越小，單向閥閥芯開度越穩(wěn)定，但是，減壓閥和氣路系統(tǒng)的最大閉環(huán)主導(dǎo)極點距虛軸的距離反而大范圍減小，收斂速度變慢，穩(wěn)定性變差，因而對單向閥止回孔直徑的選擇要進行綜合考慮.

［1］HURLBERT E，ABE J T.Methods used to investigate and resolve the space shuttle helium pressure regulator instability［C］//The 26thAIAA /ASME/SAE/ASEE Joint Propulsion Conference.Orlando:AIAA，1990:11-17.

［2］SEKITA R，MATSUDA M，NAKAMURA R.Pressure oscillation analyses of the pressure regulator for the HIIA propulsion system ［C］.The 26thAIAA/ASME/SAE/ASEE Joint Propulsion Conference.Alabama:AIAA，2003:1-9.

［3］陳曉琴.減壓閥充填過程動態(tài)特性仿真［J］.導(dǎo)彈與航天運載技術(shù)，2006(5):48-52.Chen X Q.Dynamic simulation of the pressure reducing valve in filling conditions［J］.Journal of Missile and Space Vehicle，2006(5):48-52.

［4］張雪梅，張黎輝，金廣明，等.減壓器動態(tài)過程的數(shù)學(xué)仿真［J］.航空動力學(xué)報，2004，19(4):110-114.ZHANG X M，ZHANG L H，JIN G M.Numerical simulation of the dynamic process for pressure regulator［J］.Journal of Aerospace Power，2004，19(4):110-114.

［5］尤裕榮，曾維亮.氣體減壓閥的穩(wěn)定性分析［J］.火箭推進，2009(10):34-38.YOU Y R，ZENG W L.Analysis on pneumatic pressure reducing valve stablity［J］.Journal of Rocket Propulsion，2009(10):34-38.

［6］尹文娟，魏延明.減壓閥動態(tài)響應(yīng)特性與穩(wěn)定性研究［J］.空間控制技術(shù)與應(yīng)用，2013，39(3):24-29.YIN W J，WEI Y M.Analysis on responding characteristics and stability of pneumatic pressure reducing valve［J］.Aerospace Control and Application，2013，39(3):24-29.

［7］NACI Z，GREG R L.Stability of gas pressure regulators［J］.Applied Mathematical Modeling，2008(32):61-82.