席運(yùn)志， 李軍偉， 陳雪莉， 韓磊， 王寧飛

(1.北京理工大學(xué) 宇航學(xué)院, 北京 100081； 2.西安近代化學(xué)研究所 燃燒與爆炸技術(shù)重點(diǎn)實(shí)驗(yàn)室, 陜西 西安 710065)

0 引言

隨著高能推進(jìn)劑的廣泛應(yīng)用，近年來多型固體發(fā)動(dòng)機(jī)頻繁出現(xiàn)燃燒不穩(wěn)定現(xiàn)象，嚴(yán)重影響裝備安全性和可靠性[1-2]。因此，固體發(fā)動(dòng)機(jī)燃燒不穩(wěn)定問題備受關(guān)注，亟需研究和解決。由于理論模型的不完備，當(dāng)前主要通過實(shí)驗(yàn)方法獲取推進(jìn)劑壓強(qiáng)耦合響應(yīng)函數(shù)，對發(fā)動(dòng)機(jī)的燃燒穩(wěn)定性進(jìn)行預(yù)估[3-4]。在已發(fā)展的測試方案中[5]，T型燃燒器及其改進(jìn)方案[6-9]具有結(jié)構(gòu)簡單、操作便捷等優(yōu)點(diǎn)而被廣泛應(yīng)用，但其也存在費(fèi)用高、測試誤差高達(dá)30%～50%、不易開展低頻實(shí)驗(yàn)等缺點(diǎn)[3]。為克服T型燃燒器的上述不足，Brown等[10-12]提出了旋轉(zhuǎn)閥法，該方法具有測試頻域廣、經(jīng)濟(jì)適用性好、測試結(jié)果最接近實(shí)際發(fā)動(dòng)機(jī)及可開展含鋁推進(jìn)劑實(shí)驗(yàn)等優(yōu)點(diǎn)。因此，旋轉(zhuǎn)閥法是T型燃燒器的一種良好替代或補(bǔ)充方案。

旋轉(zhuǎn)閥法對實(shí)驗(yàn)測控系統(tǒng)的測試精度要求極高，需要精確測量燃燒室振蕩壓強(qiáng)相對于旋轉(zhuǎn)閥次級(jí)排氣通道排氣面積(以下簡稱次級(jí)排氣面積)變化的相位延遲角，該相位延遲角的測試精度是影響旋轉(zhuǎn)閥法測試結(jié)果的關(guān)鍵因素[3,13-14]。為此，Brown等[10]采用外加輔助冷流裝置和接觸式探針相結(jié)合的方案間接獲取次級(jí)排氣面積，但存在探針磨損、測試誤差大、數(shù)據(jù)量多、裝置及操作流程復(fù)雜等缺點(diǎn)。上述不足限制了旋轉(zhuǎn)閥法的廣泛應(yīng)用，文獻(xiàn)[15]基于旋轉(zhuǎn)閥裝置開展了固體發(fā)動(dòng)機(jī)冷氣環(huán)境下的壓力振蕩特性研究，但并給出相位延遲角測試方案，因此，如何準(zhǔn)確測量相位延遲角是當(dāng)前亟需研究并解決的課題。

文獻(xiàn)[16-17]表明，現(xiàn)階段固體發(fā)動(dòng)機(jī)中出現(xiàn)的燃燒不穩(wěn)定多發(fā)生在200 Hz以下的低頻域。鑒于T型燃燒器對低頻域測試的不足，本文建立了一套基于旋轉(zhuǎn)閥的固體推進(jìn)劑壓強(qiáng)耦合響應(yīng)測試裝置，提出了一種非接觸式的圓光柵組件法用于測量相位延遲角。通過理論計(jì)算和冷流實(shí)驗(yàn)相互對比，驗(yàn)證了相位延遲角測試方案的有效性。同時(shí)，開展多頻率點(diǎn)的固體推進(jìn)劑點(diǎn)火實(shí)驗(yàn)，驗(yàn)證了測試系統(tǒng)的可靠性。

1 實(shí)驗(yàn)方法和理論模型

1.1 實(shí)驗(yàn)裝置

基于旋轉(zhuǎn)閥的固體推進(jìn)劑壓強(qiáng)耦合響應(yīng)函數(shù)實(shí)驗(yàn)裝置如圖 1所示，主要由小型實(shí)驗(yàn)發(fā)動(dòng)機(jī)、旋轉(zhuǎn)閥、圓光柵定位組件三部分組成。對于小型實(shí)驗(yàn)發(fā)動(dòng)機(jī)(見圖 1(a))，為保證推進(jìn)劑恒定的燃燒特性[18]，避免速度耦合或侵蝕燃燒的影響，實(shí)驗(yàn)采用端燃裝藥。發(fā)動(dòng)機(jī)左端主噴管(喉部面積12.56 mm2)用于維持燃燒室內(nèi)一定的平均壓強(qiáng)，燃燒室中部開孔并與旋轉(zhuǎn)閥相連通。在主噴管一側(cè)接入點(diǎn)火線，引燃黑火藥包產(chǎn)生的高溫高壓燃?xì)饪裳杆冱c(diǎn)燃整個(gè)推進(jìn)劑裝藥燃面。在主噴管段和燃燒室各安裝1個(gè)西安杰誠傳感器測控技術(shù)有限公司生產(chǎn)的CYG4100型高頻壓強(qiáng)傳感器(量程0～10 MPa，響應(yīng)頻率0～100 kHz)，編號(hào)為1號(hào)和2號(hào)，用于測量燃燒室內(nèi)動(dòng)態(tài)壓強(qiáng)；壓強(qiáng)數(shù)據(jù)由江蘇東華測試技術(shù)股份有限公司生產(chǎn)的DH5922D型16通道數(shù)據(jù)采集儀采集，采樣頻率為20 kHz. 其中，燃燒室內(nèi)徑及體積分別為76 mm和450 cm3，其他尺寸等比例設(shè)計(jì)。

圖1 旋轉(zhuǎn)閥實(shí)驗(yàn)裝置(點(diǎn)火)Fig.1 Ignition experimental setup with rotary valve assemble (ignition)

旋轉(zhuǎn)閥采用臺(tái)達(dá)電子企業(yè)管理有限公司生產(chǎn)的ECMA- E21320RS型伺服電機(jī)驅(qū)動(dòng)，轉(zhuǎn)速范圍為0～3 000 r/min，其內(nèi)部剖面示意如圖 1(b)所示，主要由轉(zhuǎn)子、轉(zhuǎn)子軸、定子、聯(lián)軸器等組成。其中，轉(zhuǎn)子與定子均為石墨材質(zhì)，轉(zhuǎn)子軸為中空軸，材質(zhì)為30CrMnSiA，轉(zhuǎn)子和轉(zhuǎn)子軸裝配在一起形成轉(zhuǎn)動(dòng)部件同步轉(zhuǎn)動(dòng)。同時(shí)，轉(zhuǎn)子軸與伺服電機(jī)由聯(lián)軸器連接，故轉(zhuǎn)子轉(zhuǎn)速由伺服電機(jī)控制。轉(zhuǎn)子外徑為76 mm，其周向等間距開有23個(gè)半徑為2.5 mm的轉(zhuǎn)子排氣通道(RED)。定子與發(fā)動(dòng)機(jī)燃燒室裝配在一起，燃燒室內(nèi)高壓氣體可通過定子中心位置的定子排氣通道(SED)排氣。其中SED與RED中心對齊。為消除諧波組分[9]及安裝精度因素的影響，SED截面形狀設(shè)計(jì)為矩形，其長為7.5 mm、寬為5 mm. 在轉(zhuǎn)子轉(zhuǎn)動(dòng)過程，高壓氣體可從燃燒室內(nèi)經(jīng)23個(gè)RED與SED形成的周期性“開- 閉”次級(jí)排氣通道排出，引發(fā)燃燒室內(nèi)的壓強(qiáng)振蕩。壓強(qiáng)振蕩頻率由伺服電機(jī)轉(zhuǎn)速?zèng)Q定，本裝置當(dāng)前階段適用于300 Hz以內(nèi)頻率點(diǎn)。

對于次級(jí)排氣面積，實(shí)驗(yàn)采用由英國雷尼紹有限公司生產(chǎn)的RESM20USA057型圓光柵、T2021- 30A型讀數(shù)頭、Ti0004型細(xì)分盒及美國泰道有限公司生產(chǎn)的PMAC2A- ETH型運(yùn)動(dòng)控制卡，組成的圓光柵定位組件間接獲取，如圖 1(c)所示。結(jié)合圖 1(b)，圓光柵安裝于轉(zhuǎn)子軸并隨軸同步轉(zhuǎn)動(dòng)，其外徑為57 mm，表面刻有9 000條相鄰間距為20 μm柵格線和1條光柵零點(diǎn)線。在轉(zhuǎn)動(dòng)過程中，該柵格線變化可由固定于端蓋的非接觸式讀數(shù)頭讀取。同時(shí)，細(xì)分盒與讀數(shù)頭配合使用并將柵格信號(hào)進(jìn)行4等分處理，即細(xì)分盒每輸出36 000個(gè)柵格信號(hào)則表示轉(zhuǎn)子轉(zhuǎn)動(dòng)1周，對應(yīng)排氣23次，每次排氣占用柵格信號(hào)點(diǎn)1 565個(gè)。此外，細(xì)分盒將柵格信號(hào)實(shí)時(shí)反饋給運(yùn)動(dòng)控制卡，用于修正伺服電機(jī)轉(zhuǎn)速與運(yùn)動(dòng)控制卡指令之間的誤差，即圓光柵、運(yùn)動(dòng)控制卡和伺服電機(jī)構(gòu)成閉環(huán)控制系統(tǒng)。其中，圓光柵表面的光柵零點(diǎn)線為零位校準(zhǔn)點(diǎn)，柵格信號(hào)同步傳輸給數(shù)據(jù)采集儀，為進(jìn)一步定位RED位置和計(jì)算次級(jí)排氣面積提供數(shù)據(jù)支持。

1.2 理論模型

1.2.1 壓強(qiáng)耦合響應(yīng)函數(shù)求解模型

旋轉(zhuǎn)閥實(shí)驗(yàn)裝置采用小型發(fā)動(dòng)機(jī)，其燃燒室長度遠(yuǎn)小于壓強(qiáng)振蕩波長。因此，由次級(jí)排氣面積周期性變化引發(fā)的燃燒室壓強(qiáng)振蕩在空間分布均勻，振幅處處相同，即燃燒室壓強(qiáng)處于整體振蕩狀態(tài)。同時(shí)，由于采用小容積發(fā)動(dòng)機(jī)開展實(shí)驗(yàn)，可進(jìn)一步忽略壁面熱傳導(dǎo)損失[3,10,19]。綜上，對于小振幅簡諧振蕩，文獻(xiàn)[10]及文獻(xiàn)[20]給出了固體推進(jìn)劑壓強(qiáng)耦合響應(yīng)函數(shù)實(shí)部和虛部的定義式：

(1)

式中：Rp，r和Rp，im分別為壓強(qiáng)耦合響應(yīng)函數(shù)實(shí)部與虛部；γ為燃?xì)獗葻岜?；θp,s為相位延遲角；Ω為無量綱頻率；ψa、αa分別表示排氣面積及壓強(qiáng)波動(dòng)的無量綱幅值，

(2)

(3)

固體推進(jìn)劑壓強(qiáng)耦合響應(yīng)函數(shù)的實(shí)部用于表征發(fā)動(dòng)機(jī)出現(xiàn)燃燒不穩(wěn)定的可能性，其數(shù)值越大則發(fā)動(dòng)機(jī)出現(xiàn)不穩(wěn)定燃燒的可能性越高[13-14]。由(1)式可知，響應(yīng)函數(shù)實(shí)部大小由γ、αa、ψa及θp,s共同決定，其中，γ可由美國宇航局化學(xué)平衡與應(yīng)用軟件(NASA-CEA)獲取，αa可由壓強(qiáng)傳感器獲取，ψa可由(2)式計(jì)算獲取，而θp,s由燃燒室氣體動(dòng)力學(xué)和推進(jìn)劑瞬態(tài)燃燒響應(yīng)共同決定，只能通過實(shí)驗(yàn)測量獲取。

根據(jù)(1)式及旋轉(zhuǎn)閥的工作特點(diǎn)，可使用不同喉徑的主噴管改變?nèi)紵移骄鶋簭?qiáng)，通過控制轉(zhuǎn)子轉(zhuǎn)速來調(diào)節(jié)被測頻率及排面積變化規(guī)律，從而獲得不同工作壓強(qiáng)、振蕩頻率條件下固體推進(jìn)劑壓強(qiáng)耦合響應(yīng)函數(shù)。測試方案框架如圖 2所示，其中p-t為工作壓強(qiáng)隨時(shí)間變化曲線，S-t為旋轉(zhuǎn)閥排氣面積隨時(shí)間變化曲線。

圖2 旋轉(zhuǎn)閥測試方案框架圖Fig.2 Framework of rotary valve test program

1.2.2 次級(jí)排氣面積求解模型

旋轉(zhuǎn)閥排氣面積S由主噴管面積Sn和次級(jí)排氣面積Sr組成，前者恒定，后者處于動(dòng)態(tài)變化狀態(tài)。在旋轉(zhuǎn)閥轉(zhuǎn)動(dòng)過程，RED與SED形成周期性開- 閉的次級(jí)排氣通道(見圖 3)。

圖3 旋轉(zhuǎn)閥轉(zhuǎn)動(dòng)示意圖Fig.3 Rotary motion schematic of rotary valve

圖4 次級(jí)排氣面積求解示意圖Fig.4 Schematic diagram of the solution to the secondary exhaust area

f=nrn，

(4)

式中：nr為RED個(gè)數(shù)。

旋轉(zhuǎn)閥轉(zhuǎn)動(dòng)過程中的轉(zhuǎn)子轉(zhuǎn)動(dòng)角速度為

ω=2πn.

(5)

結(jié)合圖 1(b)中的型面示意及圖 3可知，次級(jí)排氣面積Sr的大小由RED位置決定，其求解示意如圖4所示。

當(dāng)0≤Δx<2R0時(shí)，次級(jí)排氣面積Sr由0 mm2逐漸增加至最大值：

(6)

式中：

θ1=arccos(1-Δx/R0).

(7)

當(dāng)RED轉(zhuǎn)至A3位置時(shí)，次級(jí)排氣面積Sr最大，而后逐漸減小至A5位置的零值。因此，當(dāng)2R0<Δx≤4R0時(shí)，

(8)

式中：

θ2=arccos(3-Δx/R0).

(9)

圓光柵定位組件的光柵零點(diǎn)線與一個(gè)RED的中心位置對正。當(dāng)該RED處于A3位置(見圖 3(c))時(shí)，讀數(shù)頭獲取的柵格信號(hào)數(shù)Nc=0，此時(shí)Sr最大。一個(gè)排氣周期對應(yīng)的柵格信號(hào)總數(shù)Nc0=1 565，引入中間變量Nc1=(Nc-783)，則當(dāng)Nc1=0時(shí)RED處于A1位置，此時(shí)Sr=0. 因此，Δx與Nc1的關(guān)系滿足(10)式：

(10)

式中：Dr為轉(zhuǎn)子外直徑。

在旋轉(zhuǎn)閥轉(zhuǎn)動(dòng)期間，轉(zhuǎn)子上的23個(gè)RED與SED形成周期性開- 閉的排氣通道，結(jié)合(6)式～(10)式可推導(dǎo)出柵格信號(hào)數(shù)Nc1與次級(jí)排氣面積Sr之間的關(guān)系為

(11)

(12)

式中：j為已完成的排氣周期個(gè)數(shù)。

實(shí)驗(yàn)中，將圓光柵組件獲取的柵格信號(hào)數(shù)據(jù)代入(12)式即可獲取次級(jí)排氣面積變化規(guī)律，為進(jìn)一步獲取相位延遲角奠定基礎(chǔ)。

2 圓光柵定位方法有效性驗(yàn)證

相位延遲角的測試精度決定了壓強(qiáng)耦合響應(yīng)函數(shù)的準(zhǔn)確性，故點(diǎn)火實(shí)驗(yàn)前需要驗(yàn)證圓光柵定位組件測試次級(jí)排氣面積并獲取相位延遲角的有效性。為此，本文建立了次級(jí)排氣面積Sr與燃燒壓強(qiáng)振蕩之間的理論計(jì)算模型，并對應(yīng)開展了流體介質(zhì)為氮?dú)獾睦錃鈱?shí)驗(yàn)。

2.1 冷氣條件下的壓強(qiáng)振蕩計(jì)算模型

對于冷氣實(shí)驗(yàn)，燃燒室內(nèi)氣體的動(dòng)態(tài)質(zhì)量滿足(13)式：

(13)

供氣質(zhì)量流率由(14)式[20]定義為

(14)

式中：A為供氣入口面積；ξ為阻力系數(shù)[21]；pi為供氣壓強(qiáng)；Rg為氣體常數(shù)；Tg為氣體溫度。

次級(jí)排氣通道內(nèi)的氣體流動(dòng)劇烈變化且處于雍塞狀態(tài)，排氣質(zhì)量流率由(15)式[22]定義為

(15)

式中：Sr由(11)式定義；CD為排氣流量系數(shù)[22-23]，

(16)

結(jié)合氣體狀態(tài)方程，(14)式～(16)式代入(13)式，得

(17)

式中：V為燃燒室自由容積。

基于(17)式可以計(jì)算不同振蕩頻率及供氣壓強(qiáng)對燃燒室壓強(qiáng)振蕩的影響規(guī)律，理論計(jì)算參數(shù)來源于冷氣實(shí)驗(yàn)，具體如表1所示。

表1 冷氣實(shí)驗(yàn)中主要參數(shù)(298 K)Tab.1 Main parameters for cold-gas test(298 K)

2.2 冷氣實(shí)驗(yàn)

基于圖1，在1號(hào)傳感器的位置接入進(jìn)氣管路、電磁閥、減壓閥和高壓氮?dú)鈨?chǔ)氣瓶等部件，設(shè)計(jì)如圖5所示的冷氣實(shí)驗(yàn)裝置[15]。實(shí)驗(yàn)中，高壓氣體經(jīng)減壓閥降低至實(shí)驗(yàn)需求壓強(qiáng)后，進(jìn)入燃燒室并從旋轉(zhuǎn)閥排出，引發(fā)燃燒室壓強(qiáng)振蕩。其中，主噴管安裝聲阻塞不排氣，壓強(qiáng)數(shù)據(jù)及柵格信號(hào)數(shù)據(jù)分別通過線路1、線路2傳輸給數(shù)據(jù)采集儀。

圖5 旋轉(zhuǎn)閥冷氣實(shí)驗(yàn)裝置Fig.5 Schematic diagram of a cold-gas flow experimental apparatus

對于固體發(fā)動(dòng)機(jī)領(lǐng)域的燃燒不穩(wěn)定，其燃燒室壓強(qiáng)振蕩主要以近似正弦波的形式出現(xiàn)[4]。以23 Hz為例(見圖 6(a))，對其壓強(qiáng)振蕩曲線進(jìn)行非線性(正弦)擬合。結(jié)果表明：理論計(jì)算的校正決定系數(shù)為0.999 65，略大于實(shí)驗(yàn)的0.996 06. 其中，文獻(xiàn)[24]表明：校正決定系數(shù)越接近于1，擬合曲線和原始數(shù)據(jù)點(diǎn)匹配度越高。因此，理論計(jì)算和實(shí)驗(yàn)獲取的燃燒室壓強(qiáng)均按正弦波的形式變化。

為量化實(shí)驗(yàn)與理論計(jì)算之間的差異，定義壓強(qiáng)延遲時(shí)間tp,s為

tp,s=ts-tp，

(18)

對于正弦波振蕩，相位延遲角可由(19)式獲?。?/p>

θp,s=2πftp,s.

(19)

對于給定的振蕩頻率，相位延遲角θp,s由壓強(qiáng)延遲時(shí)間tp,s決定。故基于圓光柵定位組件獲取次級(jí)排氣面積和相位延遲角的方法是有效的。

圖6 冷氣實(shí)驗(yàn)與理論計(jì)算結(jié)果對比Fig.6 Comparison of experimental and theoretically calculated results of cold gas

3 固體推進(jìn)劑點(diǎn)火實(shí)驗(yàn)

3.1 固體推進(jìn)劑參數(shù)

點(diǎn)火實(shí)驗(yàn)采用密度1.62 g/cm3的不含鋁且穩(wěn)定性高的雙基推進(jìn)劑，其配方及燃速參數(shù)分別如表3和表4所示。

實(shí)驗(yàn)推進(jìn)劑裝藥結(jié)構(gòu)如圖7所示，裝藥燃面直徑72 mm，裝藥長度35 mm，包覆層厚度為2 mm. 其中，點(diǎn)火藥為2號(hào)小粒黑，藥量5 g.

表2 壓強(qiáng)延遲時(shí)間和相位延遲角(冷氣)Tab.2 Pressure delay times and phase delayangles (cold-gas)

表3 推進(jìn)劑配方參數(shù)Tab.3 Parameters of solid propellant formulation

表4 推進(jìn)劑燃速(20 ℃)Tab.4 Burning rates of solid propellant(20 ℃)

圖7 推進(jìn)劑裝藥Fig.7 Schematic diagram of solid propellant charge

基于圖1裝置和圖7裝藥開展振蕩頻率分別為23 Hz、46 Hz、69 Hz、115 Hz、138 Hz的推進(jìn)劑點(diǎn)火實(shí)驗(yàn)，獲取燃燒室內(nèi)壓強(qiáng)振蕩數(shù)據(jù)，以23 Hz、69 Hz、138 Hz為例，其燃燒室壓強(qiáng)隨時(shí)間變化如圖8所示。

圖8 燃燒室壓強(qiáng)隨時(shí)間變化曲線(點(diǎn)火)Fig.8 Variation in pressure in conbustion chamber with time (ignition test)

3.2 實(shí)驗(yàn)結(jié)果及分析

由圖8(a)可知，在t=0.16 s時(shí)，瞬時(shí)點(diǎn)火壓強(qiáng)峰值達(dá)到10 MPa，隨后在旋轉(zhuǎn)閥間歇排氣作用下，燃燒室壓強(qiáng)在2.9 MPa附近呈周期性振蕩。由局部視圖可知，噴管段(1號(hào),pn)和燃燒室(2號(hào),pc)壓強(qiáng)均按正弦波形振蕩且變化規(guī)律一致，進(jìn)一步印證燃燒室壓強(qiáng)處于整體振蕩狀態(tài)。由于1號(hào)處傳感器在實(shí)驗(yàn)時(shí)未進(jìn)行二次標(biāo)定，導(dǎo)致前者測量壓強(qiáng)數(shù)據(jù)比后者同步高出約0.16 MPa，故后續(xù)將采用2號(hào)處的壓強(qiáng)數(shù)據(jù)開展分析研究。

由于高溫、高壓燃?xì)獬掷m(xù)周期性的經(jīng)次級(jí)排氣通道從中空轉(zhuǎn)子軸排出，使得氣流所經(jīng)定子、轉(zhuǎn)子及轉(zhuǎn)子軸受熱膨脹，導(dǎo)致轉(zhuǎn)子與定子、轉(zhuǎn)子與旋轉(zhuǎn)閥體內(nèi)壁之間的摩擦阻力增大。當(dāng)旋轉(zhuǎn)閥轉(zhuǎn)動(dòng)所需的扭矩持續(xù)大于電機(jī)額定扭矩時(shí)，電機(jī)將發(fā)生停轉(zhuǎn)，燃燒室壓強(qiáng)不再呈周期性振蕩。如23 Hz工況的t=2.47 s(見圖8(a))，69 Hz的t=2.94 s(見圖8(b))，138 Hz的t=1.48 s(見圖8(c))均出現(xiàn)電機(jī)停轉(zhuǎn)現(xiàn)象。以23 Hz為例，其點(diǎn)火時(shí)刻和電機(jī)停轉(zhuǎn)時(shí)刻對應(yīng)的旋轉(zhuǎn)閥工作狀態(tài)，如圖9所示。在點(diǎn)火瞬間(t=0.16 s)，高溫高壓氣體攜帶固體顆粒從主噴管和中空轉(zhuǎn)子軸排出，安全閥(設(shè)計(jì)閾值壓強(qiáng)15 MPa)芯極(ABS材質(zhì))被點(diǎn)亮。當(dāng)發(fā)動(dòng)機(jī)持續(xù)工作至t=2.47 s時(shí)電機(jī)停轉(zhuǎn)，轉(zhuǎn)子被燒紅。

圖9 旋轉(zhuǎn)閥工作過程的瞬態(tài)畫面(23 Hz)Fig.9 Transient pictures of rotary valve in working process(23 Hz)

對于固體推進(jìn)劑壓強(qiáng)耦合響應(yīng)函數(shù)，實(shí)際所需的有效數(shù)據(jù)段為3～5個(gè)振蕩周期，如圖 8中各頻率點(diǎn)的局部視圖，因此電機(jī)停轉(zhuǎn)對實(shí)驗(yàn)測試沒有影響。

參考圖6數(shù)據(jù)處理方法，對上述5組頻率點(diǎn)壓強(qiáng)及排氣面積數(shù)據(jù)進(jìn)行處理，獲取燃燒室壓強(qiáng)波動(dòng)無量綱幅值αa、相位延遲角θp,s等壓強(qiáng)振蕩表征參數(shù)如表5所示，進(jìn)一步繪制αa和θp,s隨頻率變化曲線如圖 10所示。

表5 壓強(qiáng)振蕩表征參數(shù)Tab.5 Characteristic parameters of pressure oscillation

圖10 壓強(qiáng)波動(dòng)無量綱幅值和相位延遲角隨頻率變化曲線Fig.10 Variation in pressure oscillation dimensionless amplitude and phase delay angle with oscillation frequency

由表 5可知，23 Hz和46 Hz 2個(gè)頻率點(diǎn)的燃燒室平均壓強(qiáng)維持在2.9 MPa附近，余下3組隨著頻率的增加，平均壓強(qiáng)由69 Hz的2.764 MPa降低至138 Hz的2.523 MPa. 而壓強(qiáng)波動(dòng)幅值p′則隨著頻率的增加逐漸由23 Hz的0.278 MPa降低至138 Hz的0.036 MPa. 結(jié)合冷流實(shí)驗(yàn)(見圖6)，與p′同樣隨著頻率的增加而逐漸減小，故頻率影響燃燒室平均壓強(qiáng)和壓強(qiáng)波動(dòng)幅值。同樣地，壓強(qiáng)波動(dòng)無量綱幅值αa和相位延遲角θp,s也隨著頻率的增加而逐漸減小(見圖10)，分別由23 Hz的0.096 3和0.531π rad降低至138 Hz的0.014 3和0.501π rad.

3.3 壓強(qiáng)耦合響應(yīng)函數(shù)

圖11 壓強(qiáng)耦合響應(yīng)函數(shù)隨頻率變化Fig.11 Variation in pressure-coupled response function with oscillation frequency

4 結(jié)論

本文建立了一套基于旋轉(zhuǎn)閥的固體推進(jìn)劑壓強(qiáng)耦合響應(yīng)測試裝置，提出了一種非接觸式的圓光柵組件法用于測量相位延遲角。通過理論計(jì)算和冷流實(shí)驗(yàn)相互對比，驗(yàn)證了圓光柵組件法的有效性。同時(shí)，開展多頻率點(diǎn)的固體推進(jìn)劑點(diǎn)火實(shí)驗(yàn)，驗(yàn)證了測試系統(tǒng)的可靠性。得到如下主要結(jié)論：

1) 冷氣實(shí)驗(yàn)與理論計(jì)算結(jié)果對比誤差小于4.35%，驗(yàn)證了圓光柵定位組件方案測量相位延遲時(shí)間和相位延遲角的有效性與可行性。

2) 5次推進(jìn)劑點(diǎn)火實(shí)驗(yàn)，驗(yàn)證了旋轉(zhuǎn)閥法測試系統(tǒng)在高溫、高壓環(huán)境下工作的可靠性。

3) 點(diǎn)火實(shí)驗(yàn)獲取了5個(gè)頻率點(diǎn)下的固體推進(jìn)劑的壓強(qiáng)耦合響應(yīng)函數(shù)，驗(yàn)證了理論分析方法與數(shù)據(jù)處理流程的科學(xué)性。