劉 輝，梁其誠

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

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

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

1 理論計算

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

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

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

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

2 仿真對比分析

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

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

選擇歐拉方程，不考慮化學反應、重力以及溫度等因素。

(1)

(2)

方程的積分形式表示為：

(3)

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

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

根據(jù)內(nèi)彈道學理論，火藥燃氣作用從擊發(fā)到后效期劃分為三個階段：啟動階段、彈丸在膛內(nèi)運動階段和火藥氣體后效階段。這里主要研究火藥后效期階段流場的分布，火藥后效期的氣體平均速度根據(jù)裝藥號得出為900 m/s，因此入口定義其速度為900 m/s；時間步長定義為0.001 s；時間步數(shù)設置為10，即迭代時間為0.01 s。環(huán)境設置為常溫常壓下，流體為空氣，初始化后設置參數(shù)的收斂精度為默認值0.001[7]。

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

2.3 仿真結(jié)果分析

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

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

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

表3 原制退器和改進后制退器擋板最大動壓和后壁最大總壓

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

3 結(jié) 論

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

[1] 李鴻志,尤國釗.炮口沖擊波的形成和分布規(guī)律以及對減小沖擊波途徑的分析[J].華東工程學院學報,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)化設計理論和方法研究[D].南京：南京理工大學,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] 樂貴高,馬大為,馮勇,等.某火炮膛口流場數(shù)值仿真[J].兵工學報,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] 王仕松,鄭堅,賈長治,等.帶制退器的膛口流場數(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] 李爍.炮口制退器效率的計算與研究[D].太原：中北大學,2012.

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

[6] 張煥好,陳志華,姜孝海,等.在炮架上測試炮口制退器效率的方法[J].兵工學報,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流體分析及仿真實用教程[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)