劉 輝，梁其誠(chéng)

(海南省軍區(qū)某軍械修理所,海南 海口 570100)

炮口制退器是火炮的關(guān)鍵部件之一，其主要作用是發(fā)射時(shí)消耗一部分后坐能量，減小火炮后坐對(duì)炮架的作用力，對(duì)炮口制退器的性能影響因素進(jìn)行分析和改進(jìn)，有助于提高火炮的機(jī)動(dòng)性，增強(qiáng)戰(zhàn)斗力。國(guó)內(nèi)對(duì)炮口制退器的研究起步稍晚，但進(jìn)展很快，成果顯著[1-3]。目前炮口制退器作用機(jī)理和相關(guān)理論已相對(duì)成熟，但隨著新的火炮武器裝備的研制，炮口制退器依然是國(guó)內(nèi)外火炮研究領(lǐng)域的研究熱點(diǎn)。

衡量炮口制退器的性能的首要指標(biāo)是制退器的制退效率和炮口沖擊波超壓值，在提升制退效率的同時(shí)往往會(huì)增大炮口沖擊波超壓值，如何提高制退器效率并降低炮口沖擊波超壓值是筆者著力解決的主要問(wèn)題。通過(guò)對(duì)已有研究成果的總結(jié)和對(duì)本制退器內(nèi)部結(jié)構(gòu)的理論計(jì)算，已得知對(duì)制退器性能影響的靈敏因子，即側(cè)孔傾角的大小影響靈敏度最大，其次為側(cè)孔面積，最后為中央彈孔面積。筆者通過(guò)對(duì)制退器的側(cè)孔傾角、側(cè)孔面積和中央彈孔面積的參數(shù)進(jìn)行改進(jìn)后，應(yīng)用計(jì)算流體動(dòng)力學(xué)仿真軟件FLUENT建立原制退器和改進(jìn)制退器的二維仿真模型，仿真對(duì)比分析其效率和炮口沖擊波超壓值，證明改進(jìn)后的制退器具有更高的制退效率和較低的炮口沖擊波超壓值。

1 理論計(jì)算

應(yīng)用改進(jìn)的奧爾洛夫法對(duì)制退器進(jìn)行理論計(jì)算，各側(cè)孔角度、側(cè)孔面積和中央彈孔面積對(duì)炮口制退器效率的靈敏度的影響如圖1～圖3所示。

根據(jù)上述計(jì)算結(jié)果取制退器效率最大時(shí)的側(cè)孔傾角、側(cè)孔面積和中央彈孔面積相應(yīng)的結(jié)構(gòu)參數(shù)，其原結(jié)構(gòu)和優(yōu)化結(jié)構(gòu)參數(shù)如表1和表2所示。

表1 原炮口制退器主要參數(shù)

表2 改進(jìn)后炮口制退器主要參數(shù)

2 仿真對(duì)比分析

炮口制退器在炮彈發(fā)射瞬間，氣流為沖擊式，流動(dòng)狀態(tài)呈現(xiàn)無(wú)序的混亂形式，制退器各部位壓力、各質(zhì)點(diǎn)的速度隨時(shí)間和空間的變化發(fā)生不規(guī)則的脈動(dòng)，而且炮口氣流為多組分、氣固兩相并伴隨有化學(xué)反應(yīng)，想要建立一個(gè)準(zhǔn)確且全面的數(shù)學(xué)模型很困難[4-6]。把模型進(jìn)行適當(dāng)?shù)暮?jiǎn)化處理，炮彈發(fā)射瞬間有沖擊波，忽略炮彈對(duì)流場(chǎng)的影響，簡(jiǎn)化為炮口制退器的入口定義為速度，該研究對(duì)象為非定常流問(wèn)題，通過(guò)設(shè)計(jì)時(shí)間步長(zhǎng)和時(shí)間步數(shù)來(lái)研究火藥后效期的流場(chǎng)分布，把彈丸離開(kāi)炮口時(shí)定義為零時(shí)刻。

2.1 控制方程數(shù)學(xué)模型

選擇歐拉方程，不考慮化學(xué)反應(yīng)、重力以及溫度等因素。

(1)

(2)

方程的積分形式表示為：

(3)

非定常流的控制方程的離散格式，選為二階迎風(fēng)格式，流場(chǎng)迭代選擇Simple算法。

2.2 建立模型和定義相關(guān)參數(shù)

根據(jù)內(nèi)彈道學(xué)理論，火藥燃?xì)庾饔脧膿舭l(fā)到后效期劃分為三個(gè)階段：?jiǎn)?dòng)階段、彈丸在膛內(nèi)運(yùn)動(dòng)階段和火藥氣體后效階段。這里主要研究火藥后效期階段流場(chǎng)的分布，火藥后效期的氣體平均速度根據(jù)裝藥號(hào)得出為900 m/s，因此入口定義其速度為900 m/s；時(shí)間步長(zhǎng)定義為0.001 s；時(shí)間步數(shù)設(shè)置為10，即迭代時(shí)間為0.01 s。環(huán)境設(shè)置為常溫常壓下，流體為空氣，初始化后設(shè)置參數(shù)的收斂精度為默認(rèn)值0.001[7]。

圖4為炮口制退器的實(shí)物圖和三維仿真流體模型，圖5為改進(jìn)后制退器的三維仿真模型，在炮口制退器周圍劃分2 m×2 m的邊界區(qū)域。

2.3 仿真結(jié)果分析

圖6為原制退器和改進(jìn)后制退器在2 ms時(shí)的速度等值線圖，對(duì)比可知，改進(jìn)后的制退器流量的分配比更加均勻，各側(cè)孔出口處的氣體速度差值減少，這樣使得制退器外部的紊流強(qiáng)度減小，同理外部的壓強(qiáng)也會(huì)趨于平衡，從而使得火炮周圍的超壓值降低。

圖7為原制退器和改進(jìn)后制退器在2 ms時(shí)的動(dòng)壓等值線圖，對(duì)比可知，原制退器的動(dòng)壓主要分布在制退器第2對(duì)側(cè)孔的側(cè)后方，改進(jìn)后的制退器各側(cè)孔外圍的動(dòng)壓分布較均衡，往后側(cè)噴射的氣流被分散往制退器兩側(cè)噴射，相當(dāng)于改進(jìn)后壓力的沖擊面積增大，從而使得對(duì)單點(diǎn)的壓力值降低。

圖8為側(cè)孔擋板在0.3 ms時(shí)刻所受的動(dòng)壓力曲線圖，圖9為后壁在1 ms時(shí)刻所受的總壓曲線圖，通過(guò)數(shù)據(jù)整理得到原制退器和改進(jìn)后制退器側(cè)孔所受動(dòng)壓的最大壓力值和后壁所受的最大總壓值如表3所示。

表3 原制退器和改進(jìn)后制退器擋板最大動(dòng)壓和后壁最大總壓

制退器的制退效果主要通過(guò)制退器側(cè)孔的擋板受力實(shí)現(xiàn)，由表3可知，經(jīng)過(guò)改進(jìn)的制退器各側(cè)孔在0.3 ms和1 ms時(shí)刻的受力比原制退器受力更加均勻平衡，并且改進(jìn)的制退器所受的軸向合力更大，因此其制退效率也將更高。在0.1 ms和1 ms時(shí)刻原制退器后壁總壓分別195 kPa和78 kPa，而改進(jìn)后的制退器為39 kPa和34 kPa，所以改進(jìn)后制退器對(duì)后壁的沖擊超壓值要小得多。

3 結(jié) 論

制退器的制退力主要為氣流對(duì)側(cè)孔的軸向作用力之和，由表3得知改進(jìn)制退器對(duì)側(cè)孔的軸向合力比原制退器對(duì)側(cè)孔的軸向合力大得多，說(shuō)明改進(jìn)制退器的制退效率更高。以上仿真分析在制退器的四周建立了2 m×2 m的邊界區(qū)域，通過(guò)仿真得知改進(jìn)制退器對(duì)后壁沖擊力要比原制退器小得多，說(shuō)明改進(jìn)制退器降低了炮口沖擊波超壓值。綜上所述，改進(jìn)制退器性能比原制退器優(yōu)越，具有較高的實(shí)用價(jià)值。

[1] 李鴻志,尤國(guó)釗.炮口沖擊波的形成和分布規(guī)律以及對(duì)減小沖擊波途徑的分析[J].華東工程學(xué)院學(xué)報(bào),1977(1):26-48.

LI Hong-zhi,YOU Guo-zhao.Analysis of the formation and distribution of the muzzle shock wave and the approach to reduce the shock wave[J]. Journal of East China Institute of Engineering, 1977(1):26-48.(in Chinese)

[2] 江坤.炮口制退器優(yōu)化設(shè)計(jì)理論和方法研究[D].南京：南京理工大學(xué),2007.

JIANG Kun.Study on the optimization design theory and the method of muzzle brake[D].Nanjing:Nanjing University of Science and Technology,2007.(in Chinese)

[3] 樂(lè)貴高,馬大為,馮勇,等.某火炮膛口流場(chǎng)數(shù)值仿真[J].兵工學(xué)報(bào),2004,25(1):19-21.

LE Gui-gao,MA Da-wei,FENG Yong,et al.Numerical simulation of a gun muzzle flow field[J].Acta Armamentarii,2004,25(1):19-21.(in Chinese)

[4] 王仕松,鄭堅(jiān),賈長(zhǎng)治,等.帶制退器的膛口流場(chǎng)數(shù)值模擬[J].火力與指揮控制,2011,36(2):148-151.

WANG Shi-song,ZHENG Jian,JIA Chang-zhi,et al. Numerical simulation of the muzzle flow field with muzzle brake[J].Fire and Command Control,2011,36(2):148-151.(in Chinese)

[5] 李爍.炮口制退器效率的計(jì)算與研究[D].太原：中北大學(xué),2012.

LI Shuo. Calculation and research on muzzle brake efficiency[D].Taiyuan：North University of China,2012. (in Chinese)

[6] 張煥好,陳志華,姜孝海,等.在炮架上測(cè)試炮口制退器效率的方法[J].兵工學(xué)報(bào),2012,33 (5):623-629.

ZHANG Huan-hao,CHEN Zhi-hua,JIANG Xiao-hai, et al. The method of testing muzzle brake efficiency in the gun carriage[J]. Acta Armamentarii,2012,33 (5):623-629. (in Chinese)

[7] 朱紅鈞,林元華,謝龍漢.FLUENT流體分析及仿真實(shí)用教程[M].北京：人民郵電出版社,2010.

ZHU Hong-jun,LIN Yuan-hua,XIE Long-han.FLUENT analysis and simulation of fluid and practical tutorial[M].Beijing: People’s Posts and Telecommunications Press,2010. (in Chinese)