開口向上的環(huán)形腔寬度對彈射載荷的影響

郭佳肄，胡曉磊，孫船斌，張學鋒，馬大為

(1.安徽工業(yè)大學 a.工程實踐與創(chuàng)新教育中心； b.機械工程學院; c.計算機學院， 安徽 馬鞍山 243002；2.南京理工大學 機械工程學院， 南京 210094)

【裝備理論與裝備技術(shù)】

開口向上的環(huán)形腔寬度對彈射載荷的影響

郭佳肄1a，胡曉磊1b，孫船斌1b，張學鋒1c，馬大為2

(1.安徽工業(yè)大學 a.工程實踐與創(chuàng)新教育中心； b.機械工程學院; c.計算機學院， 安徽 馬鞍山 243002；2.南京理工大學 機械工程學院， 南京 210094)

為了研究開口向上的環(huán)形腔寬度對低溫燃氣彈射初容室內(nèi)載荷的影響，采用計算流體動力學方法，結(jié)合SST湍流模型、有限速率/渦耗散燃燒模型和域動分層動網(wǎng)格技術(shù)，建立密閉空間導彈彈射過程二次燃燒模型。在與實驗對比驗證的基礎(chǔ)上，結(jié)合數(shù)值模擬研究了開口向上的環(huán)形腔寬度對低溫燃氣彈射載荷的影響。結(jié)果表明：環(huán)形腔寬度為110 mm時，彈射載荷變化最平穩(wěn)，且無明顯雙峰現(xiàn)象。研究結(jié)果為導彈燃氣彈射的初容室結(jié)構(gòu)設(shè)計提供參考。

計算流體動力學；低溫燃氣彈射；二次燃燒；載荷

低溫燃氣彈射具有操作簡單和無需熱防護等優(yōu)勢，在國內(nèi)外導彈發(fā)射中廣泛應(yīng)用[1]。由于低溫推進劑燃燒后產(chǎn)生大量未完全燃燒的氣體組分，在空氣中產(chǎn)生二次燃燒現(xiàn)象，造成初容室內(nèi)載荷突然增加，影響導彈出筒參數(shù)。針對二次燃燒問題，國內(nèi)外主要采用有限速率/渦耗散模型進行研究。Luan[2]采用該方法結(jié)合氣固兩相流理論研究了煤的燃燒過程。Dharavath[3]采用有限速率/渦耗散模型建立了氫燃料超燃沖壓發(fā)動機三維流動模型，獲得的燃燒波速度在近尾跡區(qū)域與實驗吻合較好。遲宏偉[4]采用熱解氣體有限速率/渦耗散模型研究了沖壓發(fā)動機燃燒室中PMMA自點火性能。李仁鳳[5]運用該方法結(jié)合動網(wǎng)格技術(shù)研究了低溫推進劑燃燒產(chǎn)物對低溫彈射內(nèi)彈道與載荷的影響。胡曉磊、 樂貴高和馬大為[6]研究了環(huán)形腔的開口方向?qū)椛漭d荷和內(nèi)彈道的影響，結(jié)果表明環(huán)形腔開口向上時，燃氣彈射載荷低于無環(huán)形腔和開口向下的環(huán)形腔結(jié)構(gòu)。本文在文獻[6]的基礎(chǔ)上，研究環(huán)形腔的寬度對低溫彈射載荷的影響。

1 物理模型和計算方法

1.1 彈射裝置結(jié)構(gòu)

低溫燃氣彈射初容室結(jié)構(gòu)如圖1所示。從燃氣發(fā)生器噴出的低溫燃氣在初容室內(nèi)二次燃燒，推動尾罩向上運動。

圖1 低溫燃氣彈射初容室結(jié)構(gòu)示意圖

1.2 控制方程

由于低溫燃氣彈射初容室具有軸對稱結(jié)構(gòu)特性，采用二維軸對稱多組分N-S方程

(1)

式中符號意義參見文獻[6]。

湍流模型采用RNGk-ε湍流模型[7]，二次燃燒采用有限速率/渦耗散模型[8]，采用域動分層動網(wǎng)格技術(shù)[9]模擬導彈尾罩的運動過程。

1.3 網(wǎng)格模型和邊界條件

本文建立了三種網(wǎng)格密度流動模型，分別為10.2萬網(wǎng)格、6.7萬網(wǎng)格和2.2萬網(wǎng)格。提取P點壓力和溫度載荷作為比較對象。通過分析圖2壓力曲線和圖3溫度曲線可見，2.2萬網(wǎng)格的P點溫度與6.7萬和10.2萬網(wǎng)格相差較大。綜合數(shù)值計算經(jīng)濟性，采用6.7萬初始計算網(wǎng)格進行數(shù)值計算。圖4為包含開口向上的環(huán)形腔結(jié)構(gòu)網(wǎng)格模型。燃氣室壓力隨時間變化曲線如圖5所示。燃氣發(fā)生器噴噴管入口燃燒產(chǎn)物的組分由CEA軟件計算，噴管入口和初容室內(nèi)初始氣體質(zhì)量分數(shù)如圖6所示。

圖2 不同網(wǎng)格下觀測點壓力曲線

圖3 不同網(wǎng)格下觀測點溫度曲線

圖4 網(wǎng)格模型

圖5 入口壓力曲線

圖6 噴管入口氣體質(zhì)量分數(shù)

采用有限體積法對控制方程進行離散，壓力梯度、動量方程和湍流方程采用二階迎風格式。

1.4 數(shù)值方法驗證.

為了驗證數(shù)值方法的有效性，對比了無環(huán)形腔下觀測點P的壓力數(shù)值計算結(jié)果和實驗結(jié)果，如圖7所示。初容室內(nèi)最大壓力直接影響導彈發(fā)射時的對地載荷，其精度直接影響數(shù)值仿真的參考價值。實驗結(jié)果與仿真結(jié)果均顯示了第一個壓力峰值最大壓力所在位置，實驗結(jié)果與數(shù)值模擬結(jié)果誤差為6.02%。表明本文的數(shù)值計算方法對燃氣彈射的工程設(shè)計具有指導價值。

圖7 觀測點壓力數(shù)值計算曲線與實驗曲線

2 結(jié)果與分析

圖8和圖9分別為5種工況的發(fā)射筒內(nèi)觀測點壓力和溫度隨曲線，表1為五種環(huán)形腔寬度工況的觀測點壓力數(shù)據(jù)。

從圖8和表1可見，環(huán)形腔寬度為60 mm時，在0.261 2t0時刻，存在第一個壓力峰值，為0.739 5p0。隨著環(huán)形腔寬度的增加，第一個壓力峰值依次減小。環(huán)形腔寬度為130 mm時，觀測點壓力曲線僅存在一個壓力峰值，發(fā)生在0.726 4t0時刻，壓力為0.712 5P0。環(huán)形腔寬度在60～110 mm之間時，隨著環(huán)形腔寬度的減小，第二個壓力峰值依次減小。其中，當環(huán)形腔寬度為110 mm時，觀測點第二個壓力峰值最小，為0.680 8t0。可見環(huán)形腔寬度達到130 mm時，抑制了二次反應(yīng)過程，造成筒內(nèi)壓力出現(xiàn)“單峰”現(xiàn)象。從表1還可以看出，環(huán)形腔寬度為110 mm時，筒內(nèi)第一個壓力峰值最小，為0.704 1P0。由此可見，環(huán)形腔寬度為110 mm時，發(fā)射筒內(nèi)觀測點壓力變化最平穩(wěn)。同時，初容室內(nèi)的最大壓力峰值直接影響導彈的對地載荷，因此開口向上的環(huán)形腔寬度為110 mm時對發(fā)射地面的損傷小于其他四種工況，在一定程度上提高了發(fā)射安全性。從圖9可見，五種環(huán)形腔寬度下的觀測點溫度變化趨勢相同，都是先升高，后降低。觀測點溫度在0.87T0左右。這說明環(huán)形腔的寬度對初容室內(nèi)溫度影響較小。

綜合5種算例壓力云圖、溫度云圖以及觀測點載荷變化規(guī)律可以看出，環(huán)形腔寬度為110 mm時，發(fā)射筒內(nèi)壓力載荷變化最平穩(wěn)。

圖8 觀測點壓力曲線

圖9 觀測點溫度曲線

表1 觀測點壓力數(shù)據(jù)

3 結(jié)論

建立了5種寬度的環(huán)形腔結(jié)構(gòu)模型，通過與實驗對比，驗證了模型的有效性。研究了五種不同寬度的開口向上的環(huán)形腔結(jié)構(gòu)，數(shù)值仿真結(jié)果表明，當環(huán)形腔寬度為110 mm 時，導彈出筒過程中，初容室內(nèi)壓力無明顯“雙峰”現(xiàn)象，且壓力峰值最小，對地載荷也最小。說明對初容室環(huán)形腔寬度進行合理改進后，可以提高發(fā)射安全性。

[1] 李廣裕.戰(zhàn)略導彈彈射技術(shù)的發(fā)展[J].國外導彈與航天運載器,1990(7):38-49.

[2] LUAN YAN TSAN,CHYOU YAU-PIN,WANG TING.Numerical analysis of gasification perform via finite-rate model in a cross-type two-stage gasifier[J].International Journal of Heat and Transfer,2013,57:558-566.

[3] DHARAVATH M,MANNA P.CHAKRABORTY D.Thermo chemical exploration of hydrogen combustion in generic scramjet combustion[J].Aerospace Science and Technology,2013,24(1):264-274.

[4] 遲宏偉,魏志軍,王利和,等.固體燃料超燃沖壓發(fā)動機燃燒室中PMMA自點火性能數(shù)值研究[J].推進技術(shù),2014,35(6):799-808.

[5] 李仁鳳,樂貴高,馬大為.燃燒產(chǎn)物特性對燃氣彈射內(nèi)彈道與載荷的影響研究[J].兵工學報,2016,37(2):245-252.

[6] 胡曉磊,樂貴高,馬大為,等.環(huán)形腔對燃氣彈射發(fā)射筒二次反應(yīng)影響數(shù)值研究[J].兵工學報,2015,36(6):1024-1032.

[7] SHI T H.A newk-εEddy-viscosity model for high Reynolds number turbulent flows-model development and validation[J].Computers & Fluids，1995，24(3)：227-238.

[8] GUESSAB A,ARIS A,BOUNIF A.Simulation of turbulent piloted methane non-premixed flame based on combination of finite-rate/eddy-dissipation model[J].Mechanika,2013,19(6):657-664.

[9] 姜毅，郝繼光，傅德彬.導彈發(fā)射過程三維非定常數(shù)值模擬[J].兵工學報，2008，29(8)：911-915.

[10]胡曉磊，樂貴高，馬大為.燃氣彈射發(fā)射筒內(nèi)燃氣-空氣二次燃燒現(xiàn)象研究[J].彈道學報，2014,26(4)：76-81.

(責任編輯周江川)

Chamber Width Influence on Gas-Ejection Load

GUO Jia-yi1a, HU Xiao-lei1b, SUN Chuan-bin1b, ZHANG Xue-feng1c, MA Da-wei2

(1.a.Center of Engineering Practice and Innovation; b.School of Mechanical Engineering;c.School of Computer Science and Technology, Anhui University of Technology, Ma’anshan 243002, China;2.School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China)

To study the upward annular chamber width influence on lower-temperature gas-ejection load, the computational fluid dynamics methods, SST turbulent model, Finite-rate/Dissipation model and dynamic mesh method were adopted to establish the ejection secondary-combustion model compared with experimental results and the obstacle width influence on ejection load. Results show that while the annular chamber width is 110 mm, the load change most gently and without double peak value. Results can provide reference for designing initial chamber structure.

computational fluid dynamics; lower-temperature ejection; secondary combustion; load

2016-09-10；

郭佳肄(1986—)，女，碩士，主要從事計算流體動力學研究。

10.11809/scbgxb2016.12.010

郭佳肄，胡曉磊，孫船斌，等.開口向上的環(huán)形腔寬度對彈射載荷的影響[J].兵器裝備工程學報，2016(12):42-44.

format:GUO Jia-yi, HU Xiao-lei, SUN Chuan-bin,et al.Chamber Width Influence on Gas-Ejection Load[J].Journal of Ordnance Equipment Engineering,2016(12):42-44.

TJ768

A

2096-2304(2016)12-0042-04

修回日期：2016-09-30