劉玉磊 楊碩 崔金平 丁永強

摘要： 針對氣控喉栓式可調噴管方案的結構設計問題， 分析了影響方案的關鍵尺寸以及方案優(yōu)選的目標， 通過試驗設計和數(shù)值仿真的方法實現(xiàn)了對方案關鍵尺寸的優(yōu)化設計， 并依據該設計完成冷試原理樣機的加工以及試驗。 試驗結果表明所設計的氣控喉栓式可調噴管方案原理可行， 所采用的試驗設計和數(shù)值仿真方法能夠滿足方案設計與優(yōu)化的要求。

關鍵詞： 可調噴管； 喉栓； 試驗設計； 數(shù)值仿真； 優(yōu)化

中圖分類號： TJ763； V435文獻標識碼： A文章編號： 1673-5048（2018）04-0062-05

0引言

固體火箭發(fā)動機推力調節(jié)技術能夠根據彈道優(yōu)化的需求， 實時調節(jié)發(fā)動機推力大小， 合理分配發(fā)動機能量， 提高導彈的機動性與靈活性， 使導彈具備多任務空戰(zhàn)的能力， 滿足未來空戰(zhàn)目標多樣性的需求， 具有強烈的軍事需求背景。 在諸多發(fā)動機推力調節(jié)方案中， 喉栓式可調噴管方案以其結構簡單、 工作可靠和易于工程使用等優(yōu)點， 逐步成為該領域研究的熱點。 該方案只需在噴管中嵌入一根喉栓， 就能實現(xiàn)推力大小的調節(jié)[1-3]， 推力調節(jié)比大， 理論上具有無級調節(jié)的能力。

目前， 國內外對喉栓式可調噴管方案的研究主要集中在理論研究以及數(shù)值模擬等方面[4-10]， 對詳細的方案設計以及試驗研究方面則鮮有報道。 本文針對氣控喉栓式可調噴管方案的結構設計問題， 分析了影響方案的關鍵尺寸以及方案優(yōu)選的目標， 通過試驗設計和數(shù)值仿真的方法實現(xiàn)了對方案關鍵尺寸的優(yōu)化設計， 并依據該設計完成冷試原理樣機的加工及試驗。

1問題描述

氣控喉栓式可調噴管方案結構如圖1所示， 主要包括噴管殼體、 喉栓支撐體以及喉栓三部分， 喉栓支撐體與喉栓共同組成控制氣室， 控制氣室通過氣路與高壓氣瓶及調壓閥設備相連。

當噴管內存在氣流流動時， 喉栓外表面會受到氣流的作用， 若此時控制氣室內壓力較低， 則喉栓做打開運動， 噴管喉部面積增大。 若增加控制氣室內壓力， 則喉栓內表面合力增加， 當內表面合力大于外表面合力時， 喉栓開始做關閉運動， 噴管喉部面積減小。 在恒流量或恒燃面的情況下， 噴管喉部面積的減小會導致來流壓強的增大， 作用在喉栓外表面的壓力增大， 最終使喉栓在一個新的位置達到平衡， 喉部面積達到新的穩(wěn)定值。

從分析來看， 該方案設計的關鍵在于使喉栓在關閉過程中， 控制氣室內對應的平衡壓強呈上升趨勢， 即喉栓外表面受到的氣流作用的合力呈增加趨勢。 另外， 由于喉栓的嵌入會擾亂噴管中的氣流流動， 使噴管產生復雜的激波系以及流動分離等現(xiàn)象， 造成總壓損失[11]。 因此， 氣控喉栓式可調噴管的優(yōu)化設計需考慮喉栓的型面與噴管型面的相互作用。 減小因喉栓插入而引起的總壓損失， 也是方案設計需要考慮的問題。

2分析建模

2.1典型工況

本研究中對于可調噴管裝置方案設計的可行性驗證采用冷流試驗的辦法。 在冷流試驗中， 以空氣模擬燃氣的作用。 目前， 受冷流試驗設備條件的限制， 無法實現(xiàn)在喉栓運動過程中流量的變化（達到符合裝藥燃燒的效果）， 因此在數(shù)值仿真及原理驗證試驗中均設定空氣流量為恒定值。

2.2仿真模型及目標

本研究的主要目的在于探索可行的氣控喉栓式可調噴管的型面設計， 對于氣室支撐結構的繞流及損失不做研究， 因此問題可簡化為二維軸對稱模型， 結構示意圖如圖2所示。 其中A， Rt， L， Lk為可設計尺寸， 其他尺寸為約束尺寸。 圖2中的仿真設計參數(shù)及取值如表1所示。

其他參數(shù)說明：

喉部半徑=喉部半徑系數(shù)（cR）×喉栓半徑（R0）； 喉部直段長度=喉部直段長度系數(shù)（cL）×喉部半徑； 喉栓行程=喉栓行程系數(shù)（cK）×喉栓總行程。

針對以上仿真設計參數(shù)， 采用全因子實驗設計方法， 生成72個仿真計算樣本點并進行計算， 依據計算結果選取喉栓關閉過程中， 控制氣室內的平衡壓強呈上升趨勢， 且出口總壓恢復系數(shù)高的設計方案作為優(yōu)選方案。

3結果分析

在入口質量流量不變的情況下， 各設計參數(shù)對入口總壓的主效應如圖3所示。

航空兵器2018年第4期劉玉磊， 等： 喉栓式可調噴管方案設計與試驗研究可以看出， 入口總壓隨喉栓行程系數(shù)的增大呈加速上升的趨勢， 隨喉部半徑系數(shù)的增大呈近似線性減少的趨勢。 另外可以看出， 喉部直段長度系數(shù)和噴管收斂半角對入口總壓的影響相對較小。

各設計參數(shù)對總壓恢復系數(shù)的主效應如圖4所示， 文中總壓恢復系數(shù)定義為噴管出口總壓與來流總壓之比。

可以看出， 喉栓行程系數(shù)對總壓恢復系數(shù)的影響最大， 且隨著開度的增加， 總壓恢復系數(shù)呈加速下降的趨勢； 總壓恢復系數(shù)受噴管收斂半角和喉部半徑系數(shù)影響相對較小， 且隨著噴管收斂半角的增加呈線性下降趨勢， 隨喉部半徑系數(shù)的增加呈線性上升趨勢； 喉部直段長度系數(shù)對總壓恢復系數(shù)的影響最小， 總壓恢復系數(shù)隨喉部直段長度系數(shù)增加略有下降。

各設計參數(shù)對喉栓軸向力的主效應如圖5所示。

可以看出， 喉栓軸向力的變化主要受喉栓行程和喉部半徑系數(shù)的影響。 隨著喉栓行程的增加， 喉栓軸向力呈加速上升的趨勢， 隨喉部半徑系數(shù)的增加接近線性下降。

喉栓式可調噴管是通過調節(jié)喉栓的開度實現(xiàn)噴管喉部面積調節(jié)， 因此喉栓的開度（即喉栓行程系數(shù)）只是喉栓工作狀態(tài)的反映， 不能作為設計優(yōu)化變量。 綜合以上分析， 除了喉栓行程系數(shù)， 對總壓恢復系數(shù)和喉栓軸向力影響最大的參數(shù)是喉部半徑系數(shù)， 另外考慮到要使總壓恢復系數(shù)較高， 應選取較小的噴管收斂半角和喉部直段長度系數(shù)， 因此選取A=45°， cL=0.4時， 喉栓軸向力以及總壓恢復系數(shù)隨喉栓行程和喉部半徑系數(shù)的變化情況做進一步分析。

喉栓軸向力隨喉栓行程以及喉部半徑系數(shù)的變化情況如圖6所示。

可以看出， 在相同開度情況下， 喉栓軸向力與喉部半徑系數(shù)呈明顯的負相關； 在喉部半徑系數(shù)一定的情況下， 喉栓軸向力隨著喉栓行程的增加均呈加速上升趨勢， 但上升的幅度則隨著喉部半徑系數(shù)的增加而減小。

總壓恢復系數(shù)隨喉栓行程以及喉部半徑系數(shù)的變化情況如圖7所示。

圖中可以看出， 在相同開度情況下， 總壓恢復系數(shù)與喉部半徑系數(shù)呈明顯的負相關； 在喉部半徑系數(shù)一定的情況下， 隨著喉栓行程的增加均有下降趨勢， 但下降的幅度則隨著喉部半徑系數(shù)的增加而減小。

由上面分析可以看出， 當cR取較大值（0.7）時， 雖然在不同喉栓行程下， 總壓恢復系數(shù)都能保持最高值， 但喉栓上軸向力隨開度的上升趨勢較弱， 對氣室平衡壓力的控制提出了較高的要求， 鑒于目前試驗中采用的是手動調壓的方式， 所以很可能出現(xiàn)過調的情況， 無法達到原理驗證的目的。 而當選取較小的cR值（0.5）時， 可以看到喉栓上軸向力隨開度的上升趨勢較強， 但是較總壓恢復系數(shù)要低得多。 因此， 綜合考慮總壓損失以及保證方案可行， 最終選取喉部半徑系數(shù)為0.6。

4試驗驗證及分析

氣控喉栓式可調噴管冷試原理驗證試驗連接如圖8所示。

控制氣通過調壓閥連接高壓氣瓶， 并利用調壓閥手動控制其輸出壓力。 試驗過程中先開啟來流系統(tǒng)并調節(jié)至預先設定的壓力值， 隨后手動調節(jié)控制氣室壓力做先上升后下降的變化過程， 試驗過程中控制氣室壓力（P2）和來流壓力（P1）的變化情況如圖9所示。

圖中可以看出， 來流壓力在32 s開始趨于穩(wěn)定， 此時來流靜壓約0.83 MPa。 控制氣室調壓閥在43 s開始動作， 之后控制氣室壓力迅速上升， 到約57.5 s時二者基本平衡。 進一步提高控制氣室壓力， 則喉栓開始關閉運動， 喉部面積減小， 來流壓力也隨之增大。 在96 s時控制氣室壓力達到最大值1.9 MPa， 此時來流壓力也達到最大值2.3 MPa。 隨后調節(jié)控制氣室壓力逐步下降， 同樣的來流壓力也隨之呈下降趨勢， 到132 s時二者重新達到一致， 此時喉栓達到最大開度， 隨后進一步降低控制氣室壓力， 來流壓力亦不隨之變化。

取喉栓動作過程中來流壓力（P1）與控制氣室壓力（P2）的比值如圖10所示。

從圖中可以看出， 在喉栓動作段內來流壓力與控制氣室壓力之比呈先增加后減少的變化趨勢。 在喉栓關閉段內， 喉部開度逐漸減小， 壓比則不斷增加， 反之， 在喉栓打開段內， 喉部開度逐漸變大， 壓比則減小。 在打開和關閉的過程中， 壓比對同一開度的重復性較好。

5結論

噴管總壓恢復系數(shù)受喉栓行程的影響最大， 且隨著喉栓行程的增加， 總壓恢復系數(shù)呈加速下降趨勢； 受噴管收斂半角和喉部半徑系數(shù)影響相對較?。?喉部直段長度系數(shù)對總壓恢復系數(shù)的影響很小， 可忽略。

喉栓表面力的變化主要受喉栓行程和喉部半徑系數(shù)的影響， 且隨著喉栓行程的增加， 喉栓表面力呈加速上升的趨勢， 隨喉部半徑系數(shù)的增加接近線性下降。

試驗結果表明， 本文提出的氣控喉栓式可調噴管方案， 能夠通過調節(jié)控制氣室壓強實現(xiàn)對喉栓位置的調節(jié)， 最終實現(xiàn)噴管喉部面積的調節(jié)。 所采用的試驗設計和基于仿真的優(yōu)化方法能夠滿足方案設計的要求。

參考文獻：

[1] 張淑慧， 胡波， 孟雅桃. 推力可控固體火箭發(fā)動機應用及發(fā)展[J]. 固體火箭技術， 2002， 25（4）： 12-15.

Zhang Shuhui， Hu Bo， Meng Yatao. Application and Development of Controllable Solid Rocket Motor[J]. Journal of Solid Rocket Technology， 2002， 25（4） ： 12-15.（in Chinese）

[2] Smithkent R， Loh HT， Chwalowski P. Analytical Contouring of Pintle Nozzle Exit Cone Using Computational Fluid Dynamics， AIAA 1995-2877[R]. 1995.

[3] 李娟， 李江. 喉栓式固體火箭發(fā)動機噴管性能影響研究[J]. 彈箭與制導學報， 2007， 27（3）： 154-160.

Li Juan， Li Jiang. Analysis for Nozzle Performance of Pintle Controllable Thrust Solid Rocket Motor[J]. Journal of Projectiles， Rockets， Missiles and Guidance， 2007， 27（3）： 154-160.（in Chinese）

[4] 唐金蘭， 宋慧敏， 李進賢， 等. 基于動網格的喉栓式推力可調噴管內流場數(shù)值模擬[J].固體火箭技術， 2014， 37（5）： 634-639.

Tang Jinlan， Song Huimin， Li Jinxian， et al. Numerical Simulation of PintleControlled Nozzle Flow Field Based on Dynamic Grid[J]. Journal of Solid Rocket Technology， 2014， 37（5）： 634-639.（in Chinese）

[5] Moureau V， Lartigue G， Sommerer Y， et al.Numerical Methods for Unsteady Compressible MultiComponent Reacting Flows on Fixed and Moving Grids[J].Journal of Computational Physics， 2005， 202（2）： 710-736.

[6] 黃河峽， 譚慧俊， 周唯陽， 等.一種帶中心體的TBCC可調噴管的設計與仿真[J].南京航空航天大學學報， 2013， 45（5）： 658-664.

Huang Hexia， Tan Huijun， Zhou Weiyang， et al. Design and Computational Study of TBCC Variable Nozzle with Central Body[J]. Journal of Nanjing University of Aeronautics & Astronautics， 2013， 45（5）： 658-664.（in Chinese）

[7] 張留歡， 徐驚雷， 莫建偉. 二元非對稱噴管可調方案試驗研究[J].航空學報， 2013， 34（4）： 772-778.

Zhang Liuhuan， Xu Jinglei， Mo Jianwei. Experimental Study of 2D Adjustable Asymmetric Nozzles[J]. Acta Aeronautica et Astronautica Sinica， 2013， 34（4）： 772-778.（in Chinese）

[8] 范志鵬， 徐驚雷， 呂鄭， 等.型面旋轉變馬赫數(shù)風洞噴管的優(yōu)化設計[J].航空學報， 2014， 35（5）： 1216-1225.

Fan Zhipeng， Xu Jinglei， Lü Zheng， et al. Optimization Design of Variable Mach Number Wind Tunnel Nozzle by Rotating Profile[J]. Acta Aeronautica et Astronautica Sinica， 2014， 35（5）： 1216-1225.（in Chinese）

[9] 李娟， 李江， 王毅林， 等.喉栓式變推力發(fā)動機性能研究[J].固體火箭技術， 2007， 30（6）： 505-509.

Li Juan， Li Jiang， Wang Yilin， et al. Study on Performance of Pintle Controlled Thrust Solid Rocket Motor[J]. Journal of Solid Rocket Technology， 2007， 30（6）： 505-509.（in Chinese）

[10] 滑利輝， 田維平， 甘曉松， 等.喉栓式推力可調發(fā)動機噴管流場數(shù)值模擬[J].固體火箭技術， 2008， 31（4）： 344-349.

Hua Lihui， Tian Weiping， Gan Xiaosong， et al. Numerical Simulation on Steady Flow Field of Variable Thrust Motor Nozzle with Pintle[J]. Journal of Solid Rocket Technology， 2008， 31（4）： 344-349.（in Chinese）

[11] 成沉， 鮑福廷， 劉旸， 等.基于響應面法的喉栓式噴管型面優(yōu)化設計[J].固體火箭技術， 2014， 37（1）： 47-51.

Cheng Chen， Bao Futing， Liu Yang， et al. Optimization Design for Contour of Pintle Nozzle Based on Response Surface Method[J]. Journal of Solid Rocket Technology， 2014， 37（1）： 47-51.（in Chinese）

Abstract： The structural design of the pintlecontrolled regulated nozzle is discussed. The key dimensions affecting the design and the objective of the optimization are analyzed. Through the experimental design and numerical simulation method， the optimal design of the key dimensions is realized. According to the design， the principle prototype and the cold experimental verification are finished. The experimental results show that the design principle of the pintlecontrolled regulated nozzle is feasible， and the simulation method can meet the requirements of the design and optimization of the program.

Key words： regulated nozzle； pintlecontrolled； experimental design； numerical simulation； optimization