彈道導(dǎo)彈燃?xì)鈴椛浼夹g(shù)

程洪杰, 陳 力, 林 睿

(火箭軍工程大學(xué)兵器發(fā)射理論與技術(shù)國(guó)家重點(diǎn)學(xué)科實(shí)驗(yàn)室, 西安 710025)

0 引言

隨著現(xiàn)代空間偵察技術(shù)的快速發(fā)展,對(duì)導(dǎo)彈隱蔽生存和機(jī)動(dòng)性也提出了更高的要求,彈射相對(duì)于自力發(fā)射能有效提高導(dǎo)彈速度和射程,簡(jiǎn)化陣地設(shè)備以提高地面機(jī)動(dòng)能力,減少燒蝕和沖擊,降低陣地紅外特征[1],因而對(duì)滿足導(dǎo)彈特定作戰(zhàn)需求有巨大優(yōu)勢(shì)。

早期的導(dǎo)彈彈射技術(shù)主要采用液壓式彈射和壓縮空氣式彈射,20世紀(jì)50年代末期以來,由于固體火箭發(fā)動(dòng)機(jī)相關(guān)技術(shù)的迅速發(fā)展以及內(nèi)彈道學(xué)理論的更新和完善,使彈射裝置小型化成為可能,燃?xì)鈴椛浼夹g(shù)和燃?xì)庹羝麖椛浼夹g(shù)應(yīng)運(yùn)而生[2]。

1 基于零維內(nèi)彈道模型的研究方法

零維內(nèi)彈道模型在定量分析內(nèi)彈道性能變化趨勢(shì)及結(jié)構(gòu)設(shè)計(jì)優(yōu)化上具有明顯優(yōu)勢(shì)。文獻(xiàn)[2]和文獻(xiàn)[3]給出了經(jīng)典內(nèi)彈道學(xué)的整套計(jì)算和分析方法,為內(nèi)彈道設(shè)計(jì)及優(yōu)化奠定了理論基礎(chǔ)。

1.1 內(nèi)彈道性能影響因素研究

高、低壓室壓強(qiáng)曲線及最大壓強(qiáng)、導(dǎo)彈的最大加速度及出筒速度是衡量?jī)?nèi)彈道性能的重要特征量。影響內(nèi)彈道性能的主要因素有:

1)高壓室裝藥。火藥的燃速直接決定了高壓室壓強(qiáng)曲線的形狀及其最大值,影響火藥燃速的主要因素如表1所示。

2)低壓室初始容積。當(dāng)?shù)蛪菏业某跏既莘e減小時(shí),低壓室壓強(qiáng)上升速度加快,沖擊振動(dòng)效應(yīng)變強(qiáng),彈底燃?xì)鉁囟戎饾u增大,初始容積的殘余空氣具有阻隔高溫燃?xì)鉁p少?gòu)椀谉g的作用[7];倘若初始容積過大,低壓室壓力上升緩慢,將嚴(yán)重影響導(dǎo)彈出筒速度。此外,初容較小還將會(huì)導(dǎo)致噴管流動(dòng)狀態(tài)為非臨界狀態(tài),對(duì)高壓室形成反壓作用,導(dǎo)致高壓室壓力出現(xiàn)波紋狀波動(dòng)[4]。

3)噴喉截面積。喉部直徑是彈射器內(nèi)彈道設(shè)計(jì)中的重要參數(shù),喉部通氣面積直接決定燃?xì)獾牧髁看笮?從而影響對(duì)導(dǎo)彈的做功能力。噴口越小,離筒速度越大,但噴管過小則后噴不足,導(dǎo)彈過載和后座加大,嚴(yán)重影響發(fā)射精度[8]。

4)環(huán)境溫度。文獻(xiàn)[9]通過燃?xì)獍l(fā)生器在高溫(+60 ℃)、常溫(+20 ℃)和低溫(-40 ℃)的熱試車實(shí)驗(yàn)和高、常、低溫的彈射實(shí)驗(yàn),發(fā)現(xiàn)燃?xì)獍l(fā)生器的壓強(qiáng)峰隨溫度升高而變大,低溫彈射后勁不足,影響導(dǎo)彈加速度和出筒速度。

表1 影響高壓室壓強(qiáng)變化的因素

1.2 內(nèi)彈道結(jié)構(gòu)優(yōu)化設(shè)計(jì)研究

彈射裝置的結(jié)構(gòu)尺寸對(duì)工質(zhì)的流場(chǎng)形態(tài)有重要影響,所以結(jié)構(gòu)優(yōu)化成為研究?jī)?nèi)彈道性能的重要方向。合理選擇裝藥和噴管尺寸,可以消除發(fā)射過程中導(dǎo)彈的加速度高頻脈沖過載[10],而推遲后噴、增加發(fā)射筒端面積和減小尾噴管喉部面積可以提高導(dǎo)彈出筒速度[8]。

文獻(xiàn)[11]針對(duì)低壓室峰值壓強(qiáng)過大的問題,采用增大低壓室初始容積、減小主裝藥初始燃面和增大噴管喉徑的方法,有效的解決了低壓室壓強(qiáng)過大的問題。文獻(xiàn)[12]改進(jìn)單活塞氣缸式彈射裝置,設(shè)計(jì)了雙級(jí)提拉式彈射裝置,增加導(dǎo)彈有效加速行程,在保證低過載的情況下,提高了導(dǎo)彈出筒速度。

1.3 零維內(nèi)彈道的“修正”方法

1)試湊法

文獻(xiàn)[9]通過實(shí)驗(yàn)獲取了彈射器在高溫(+60 ℃)和低溫(-40 ℃)下的“壓強(qiáng)-時(shí)間”和“加速度-時(shí)間”曲線,并用經(jīng)典內(nèi)彈道學(xué)作出理論曲線,用試湊法修正理論曲線,使其逐步逼近試驗(yàn)曲線。為經(jīng)典內(nèi)彈道學(xué)的完善提供了新的思路。

2)“半實(shí)驗(yàn)半理論”的方法

文獻(xiàn)[13]采用燃?xì)獍l(fā)生器實(shí)驗(yàn)數(shù)據(jù)繪成的高壓室的壓強(qiáng)-時(shí)間曲線,以此替代經(jīng)典內(nèi)彈道學(xué)從燃面理論推導(dǎo)的曲線,從而建立低壓室內(nèi)彈道及運(yùn)動(dòng)學(xué)求解方程組,經(jīng)仿真,與純理論計(jì)算相比,誤差小于2%,對(duì)于燃?xì)獍l(fā)生器的設(shè)計(jì)具有指導(dǎo)意義。

2 多維模型的數(shù)值仿真研究

隨著計(jì)算機(jī)性能的提高以及計(jì)算流體力學(xué)(CFD)理論尤其是動(dòng)網(wǎng)格技術(shù)的發(fā)展,使得多維數(shù)值計(jì)算的方法在內(nèi)彈道研究中得到廣泛應(yīng)用[14-15]。如圖1所示,彈射內(nèi)彈道模型對(duì)稱特性顯著,故通常采用二維軸對(duì)稱模型。

圖1 燃?xì)鈴椛溲b置結(jié)構(gòu)示意圖

2.1 零維模型與二維模型優(yōu)缺點(diǎn)

零維理論不考慮燃燒物之間化學(xué)反應(yīng),認(rèn)為氣體在空間上是均勻的,所以并不能精確描述內(nèi)彈道變化過程和流場(chǎng)特性[7],但零維模型具有計(jì)算量小、速度快的優(yōu)勢(shì)。二維軸對(duì)稱模型能準(zhǔn)確地給出高、低壓室的流動(dòng)細(xì)節(jié)和壓力分布情況[16],二者計(jì)算結(jié)果相差不大[17]。

2.2 二次燃燒現(xiàn)象研究

文獻(xiàn)[16]和文獻(xiàn)[17]是單相燃?xì)鈴椛鋬?nèi)彈道模型,但燃?xì)獠⒉皇菃我唤M分,常帶有大量富燃?xì)怏w[18],與低壓室氧氣發(fā)生二次燃燒現(xiàn)象[19],形成了初始?jí)簭?qiáng)峰,造成流場(chǎng)溫度、壓力和荷載上升,出筒時(shí)間縮短,出筒加速度降低,出筒速度提高[20]。文獻(xiàn)[21]進(jìn)一步研究了噴管入口燃燒產(chǎn)物壓力和組分濃度比值對(duì)內(nèi)彈道和荷載的影響,得出結(jié)論:隨著噴管入口壓力增大,低壓室氧氣消耗時(shí)間縮短,導(dǎo)彈出筒時(shí)間縮短,出筒速度增加,加速度峰值增大;隨著噴管入口CO與H2濃度比值變大,低壓室氧氣消耗時(shí)間變長(zhǎng),導(dǎo)彈出筒速度減小,加速度峰值減小。

2.3 結(jié)構(gòu)對(duì)二次燃燒的影響

含有環(huán)形腔結(jié)構(gòu)會(huì)改變?nèi)細(xì)饬鲾U(kuò)散方向,減小燃?xì)馀c空氣的接觸面積,有效降低二次燃燒的沖擊;環(huán)形腔開口向上時(shí),壓強(qiáng)和加速度變化趨于緩慢,導(dǎo)彈出筒速度變小,導(dǎo)彈出筒時(shí)間延長(zhǎng)[22]。文獻(xiàn)[23]表明壁面障礙物會(huì)降低二次壓力峰的沖擊,延遲導(dǎo)彈出筒時(shí)間,減小出筒速度,當(dāng)障礙物高度為300 mm時(shí),荷載變化平穩(wěn),出筒時(shí)間較短且出筒速度較高。

3 縮比實(shí)驗(yàn)

對(duì)于大型彈射器,由于經(jīng)濟(jì)、技術(shù)等條件的限制,使得用縮比模型試驗(yàn)代替實(shí)物原型試驗(yàn)的方法得到廣泛應(yīng)用[24-25]。縮比有如圖2的兩種方案,文獻(xiàn)[26]運(yùn)用“化四維因次為三維因次”的方法,從理論上推導(dǎo)了模型Ⅰ(奇異相似)和模型Ⅱ(幾何相似)與原型的內(nèi)彈道參量對(duì)應(yīng)相等的條件。奇異相似方案在直觀考察裝藥、高壓室工作壓強(qiáng)等方面具有優(yōu)勢(shì),宜優(yōu)先選用[2]。

圖2 發(fā)射筒系統(tǒng)模型

4 目前仍存在的問題及發(fā)展趨勢(shì)

4.1 多變量的耦合研究

在導(dǎo)彈彈射過程中,藥柱質(zhì)量誤差、侵蝕燃燒、燃面燃燒規(guī)律、噴喉尺寸、低壓室端面積及初始容積均對(duì)內(nèi)彈道性能有直接影響,高壓室的長(zhǎng)徑比也對(duì)內(nèi)彈道壓力曲線有重要影響[27],且上述因素之間互相具有耦合關(guān)系。現(xiàn)有文獻(xiàn)雖然有研究少量耦合現(xiàn)象,但變量數(shù)較少,多變量的耦合尚待進(jìn)一步研究,如:低壓室初始容積過大,二次燃燒現(xiàn)象嚴(yán)重,彈底溫度過高,且壓力上升緩慢;低壓室初始容積過小,則對(duì)彈體沖擊效應(yīng)明顯,過載過大,且低壓室易出現(xiàn)較大壓強(qiáng)峰,對(duì)高壓室造成反壓,同時(shí)也給發(fā)射筒增加強(qiáng)度負(fù)擔(dān)。如何調(diào)整高壓室尺寸、裝藥設(shè)計(jì)、噴管尺寸和低壓室尺寸,使內(nèi)彈道性能達(dá)到最優(yōu),這具有重要意義。

4.2 現(xiàn)代智能算法在線優(yōu)選結(jié)構(gòu)參數(shù)

由于內(nèi)彈道結(jié)構(gòu)參數(shù)耦合及設(shè)計(jì)方程組欠約束,傳統(tǒng)的內(nèi)彈道設(shè)計(jì)需多次試算逐步逼近可行方案,然而并不能保證為最優(yōu)方案[2]。所以修正零維模型使其正向計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合,再利用強(qiáng)大的現(xiàn)代智能算法對(duì)結(jié)構(gòu)參數(shù)進(jìn)行在線優(yōu)選匹配,對(duì)使內(nèi)彈道性能、能量利用效率及空間利用效率耦合達(dá)到最優(yōu)具有重要意義。

4.3 一體化內(nèi)彈道模型

在傳統(tǒng)的CFD研究中,低壓室流場(chǎng)入口邊界條件通常來自燃?xì)獍l(fā)生器的實(shí)驗(yàn)數(shù)據(jù),在超臨界狀態(tài)下,這種方法是可行的,但實(shí)際情況中可能出現(xiàn)亞臨界狀態(tài),此時(shí)方法就存在誤差,所以考慮推進(jìn)劑的燃燒過程,并建立一體化內(nèi)彈道模型顯得非常必要。

4.4 多相流的研究

鋁粉具有高燃燒熱,加入固體推進(jìn)劑中,可以提高火箭發(fā)動(dòng)機(jī)的比沖,減少不穩(wěn)定燃燒[28]。添加鋁粉后的推進(jìn)劑燃燒除了產(chǎn)生氣固兩相燃?xì)?還帶有熔融狀態(tài)的Al2O3液滴,之后Al2O3液滴放出熱量相變?yōu)楣腆w顆粒并聚集結(jié)塊[29-30],撞擊在燃燒室、發(fā)射筒上、導(dǎo)彈底部,形成侵蝕,鋁粉的燃燒過程還會(huì)引起一氧化碳、二氧化碳和氫氣等組分以及溫度場(chǎng)的變化[31]。多相流中顆粒、液滴、氣泡在壁面上的碰撞和聚集對(duì)內(nèi)彈道的溫度場(chǎng)、速度場(chǎng)及壓力場(chǎng)均有重要影響[32],對(duì)于多相流如何影響彈射內(nèi)彈道性能及荷載特性,尚待進(jìn)一步研究。

4.5 低燃溫推進(jìn)劑對(duì)內(nèi)彈道影響研究

燃?xì)鉁囟冗^高導(dǎo)致對(duì)導(dǎo)彈底部有燒蝕是燃?xì)鈴椛浼夹g(shù)最大的弊端,隨著低燃溫火藥技術(shù)的成熟,燃?xì)獍l(fā)生劑燃速在15～35 mm/s可調(diào),燃燒溫度在1 000～1 500 ℃[33]。為了降低燃溫,推進(jìn)劑里面添加了降速降溫劑以及多種催化劑,或者采用改變金屬(催化劑)粒徑的方法[34-35],這些物質(zhì)的添加會(huì)增加燃?xì)饨M分的復(fù)雜性,目前針對(duì)低燃溫推進(jìn)劑對(duì)內(nèi)彈道性能的影響研究相對(duì)較少。

5 結(jié)束語

現(xiàn)行燃?xì)鈴椛涞难芯慷嗖捎美碚摵蛯?shí)驗(yàn)相結(jié)合的方式,盡管已取得諸多成果,但有關(guān)文獻(xiàn)的研究仍然存在諸多不足。文中從多變量耦合、現(xiàn)代智能算法在線優(yōu)選結(jié)構(gòu)參數(shù)、建立一體化內(nèi)彈道模型、多相流以及低燃溫推進(jìn)劑5個(gè)方面指出了未來需要解決的問題及發(fā)展的趨勢(shì)。

參考文獻(xiàn):

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

[2] 譚大成. 彈射內(nèi)彈道學(xué) [M]. 北京: 北京理工大學(xué)出版社, 2015: 2-3.

[3] 袁曾鳳. 火箭導(dǎo)彈彈射內(nèi)彈道學(xué) [M]. 北京: 北京工業(yè)學(xué)院出版社, 1987: 112-132.

[4] 謝偉, 王漢平. 提拉式彈射內(nèi)彈道特性的影響因素分析 [J]. 固體火箭技術(shù), 2016, 39(1): 146-150.

[5] 王志健, 杜佳佳. 動(dòng)網(wǎng)格在固體火箭發(fā)動(dòng)機(jī)非穩(wěn)態(tài)工作過程中的應(yīng)用 [J]. 固體火箭技術(shù), 2008, 31(4): 350-353.

[6] 王天輝, 陳慶貴, 何超. 燃?xì)獍l(fā)生器內(nèi)彈道設(shè)計(jì)計(jì)算 [J]. 現(xiàn)代防御技術(shù), 2014, 42(2): 56-60.

[7] 牛鈺森. 自彈式發(fā)射內(nèi)彈道流場(chǎng)特性研究 [D]. 北京: 北京理工大學(xué), 2016: 111-114.

[8] 郝曉琴, 畢世華. 高速導(dǎo)彈彈射內(nèi)彈道優(yōu)化設(shè)計(jì)研究 [J]. 彈箭與制導(dǎo)學(xué)報(bào), 2005, 25(3): 520-523.

[9] 劉可平. 某型彈射器高/低溫工作特性研究 [D]. 上海: 上海交通大學(xué), 2008: 61-71.

[10] EDQUIST C T. Prediction of the launch pulse for gas-generator-launched missiles [J]. Journal of Propulsion and Power, 1990, 6(6): 705-712.

[11] 楊文, 屠小昌, 陳靜, 等. 某小型彈體發(fā)射系統(tǒng)低壓發(fā)射室降壓方案研究 [J]. 火工品, 2016(1): 1-5.

[12] 白鵬英, 喬軍. 雙級(jí)氣缸式彈射裝置內(nèi)彈道分析 [J]. 現(xiàn)代防御技術(shù), 2007, 35(4): 44-49.

[13] 惠衛(wèi)華, 鮑福廷, 劉旸. 考慮低燃溫燃?xì)獍l(fā)生器試驗(yàn)的彈射器內(nèi)彈道性能預(yù)示 [J]. 固體火箭技術(shù), 2013, 36(6): 715-719.

[14] LIU Yongquan, XI Anmin, LIU Hongfei, et al. An interior ballistic simulation of the gas-steam missile ejection [C]∥IEEE. 2010 Wase International Conference on Information Engineering. [S.l.]:IEEE, 2010: 235-237.

[15] LIU Y Q, XI A M, LIU H F. Numerical study of interior trajectory of missile ejection based on mechanics [J]. Advanced Materials Research, 2012, 485: 616-619.

[16] 趙雁鵬, 王學(xué)智, 劉少偉, 等. 垂直彈射系統(tǒng)彈射初始階段內(nèi)流場(chǎng)數(shù)值分析 [J]. 空軍工程大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 15(4): 5-8.

[17] 譚大成, 苗佩云. 彈射器低壓室二維內(nèi)彈道模型及數(shù)值研究 [J]. 彈箭與制導(dǎo)學(xué)報(bào), 2006, 26(4): 224-226.

[18] 張光喜, 周為民, 張鋼錘, 等. 固體火箭發(fā)動(dòng)機(jī)尾焰流場(chǎng)特性研究 [J]. 固體火箭技術(shù), 2008, 31(1): 19-23.

[19] 胡曉磊, 樂貴高, 李仁鳳, 等. 燃?xì)鈴椛浒l(fā)射筒內(nèi)燃?xì)?空氣二次燃燒現(xiàn)象研究 [J]. 彈道學(xué)報(bào), 2014, 26(4): 76-81.

[20] 胡曉磊, 王輝, 樂貴高, 等. 二次燃燒對(duì)燃?xì)鈴椛漭d荷和內(nèi)彈道影響數(shù)值研究 [J]. 固體火箭技術(shù), 2015, 38(6): 776-781.

[21] 李仁鳳, 樂貴高, 馬大為. 燃燒產(chǎn)物特性對(duì)燃?xì)鈴椛鋬?nèi)彈道與載荷的影響研究 [J]. 兵工學(xué)報(bào), 2016, 37(2): 245-252.

[22] 胡曉磊, 樂貴高, 馬大為, 等. 環(huán)形腔對(duì)燃?xì)鈴椛涑跞菔叶稳紵绊憯?shù)值研究 [J]. 兵工學(xué)報(bào), 2015, 36(6): 1024-1032.

[23] 李仁鳳, 樂貴高, 馬大為, 等. 壁面障礙物對(duì)燃?xì)鈴椛淞鲌?chǎng)和內(nèi)彈道的影響 [J]. 固體火箭技術(shù), 2016, 39(3): 444-450.

[24] 郝海鵬. 同心筒內(nèi)氣流特性的實(shí)驗(yàn)與數(shù)值研究 [D]. 北京: 北京理工大學(xué), 2016.

[25] 呂翔, 李江, 陳劍, 等. 變深度水下發(fā)射系統(tǒng)內(nèi)彈道實(shí)驗(yàn)研究 [J]. 固體火箭技術(shù), 2012, 35(1): 24-28.

[26] 邵徉. 發(fā)射筒內(nèi)彈道的相似準(zhǔn)則 [J]. 導(dǎo)彈與航天運(yùn)載技術(shù), 1995(1): 58-68.

[27] CHENG C, ZHANG X. Numerical modeling and investigation of two-phase reactive flow in a high-low pressure chambers system [J]. Applied Thermal Engineering, 2016, 99: 244-252.

[28] SHAFQAT W, XIE K, LIU Y. Numerical simulation of multi-phase combustion flow in solid rocket motors with metalized propellant [J]. Journal of Aerospace Power, 2009, 24(7): 1654-1660.

[29] PRICE E W, SIGMAN R K. Combusion of aluminized solid propellants [M]. [S.l.]:AIAA, 2000: 663-687.

[30] DRAZIN P G, REID W H. Hydrodynamic stability [M]. Cambridge: The Press Syndicate of the University Press, 2004.

[31] 劉叢林, 郜冶, 賀征. 鋁顆粒在多相流中燃燒的數(shù)值模擬 [J]. 彈箭與制導(dǎo)學(xué)報(bào), 2012, 32(3): 118-120.

[32] RASHAD M, ZHANG X B, SADEK H E. Interior ballistics two-phase reactive flow model applied to large caliber guided projectile-gun system [J]. Advances in Mechanical Engineering, 2014, 2014(2): 677-681.

[33] 喬應(yīng)克, 魯國(guó)林. 導(dǎo)彈彈射用低溫燃?xì)獍l(fā)生劑技術(shù)研究 [C]∥中國(guó)宇航學(xué)會(huì): 中國(guó)宇航學(xué)會(huì)固體火箭推進(jìn)第22屆年會(huì)論文集(推進(jìn)劑分冊(cè)). [出版地不詳]:[出版者不詳], 2005: 194-198.

[34] 秦能, 賈延斌. 過渡金屬催化劑對(duì)低燃溫雙基推進(jìn)劑性能的影響 [J]. 火炸藥學(xué)報(bào), 2014, 37(2): 73-77.

[35] 吳雄崗, 李笑江, 宋桂賢, 等. 鋁粉粒徑對(duì)改性雙基推進(jìn)劑燃燒性能的影響 [J]. 火炸藥學(xué)報(bào), 2010, 33(3): 80-83.