MEFP式戰斗部正交試驗優化

孔毓琦，郭銳

（南京理工大學機械工程學院，江蘇南京 210094）

0 引言

多彈頭爆炸成型彈丸（multiple explosively formed projectile，MEFP），與一般的預制破片戰斗部技術相比，它具有對目標打擊毀傷率高、遠距離打擊能力強等特點?，F階段對于MEFP的研究多偏向于對其裝藥結構和起爆方式等方面，對于MEFP藥形罩的結構設計也多偏向于傳統的軸向布置，而圓周方向布置的MEFP尚缺乏系統的研究。在此，設計了一種新型周向式MEFP戰斗部，并選取其裝藥材料、曲率半徑、材料和壁厚等結構參數進行正交試驗分析，從而求得其在這種結構上的優化方案。

1 MEFP戰斗部優化設計方案（圖1）（圖2）

將藥形罩的裝藥類型和結構尺寸作為優化設計的對象。因此，可以把藥形罩的罩頂厚δ，裝藥類型N，內表面曲率半徑r1以及外表面曲率半徑r2設定為正交優化設計的4個因素，各個水平因素的最終取值見表1。

表1 正交設計各因素水平取值表

這里，用L9（34）正交表將罩頂厚、裝藥類型、內表面曲率半徑以及外表面曲率半徑4個因素所對應的9個水平列成9個數值計算方案如表2所示。

表2 各仿真計算方案因素取值

仿真完成后將各方案仿真得到的質量、剩余速度和剩余比動能等仿真計算結果作為各方案的評定指標。

2 MEFP成型過程仿真

在此，按照表2所列戰斗部結構方案，分別建立有限元仿真計算模型，計算得到不同方案的EFP質量m，剩余速度vf和剩余比動能Es等模擬結果，根據參數指標對計算結果進行極差分析，就可以得到優化方案的各組成因素。下面以方案3為例來進行仿真模擬。在此方案中采用多列藥形罩進行疊放，藥型罩截面為變壁厚球缺形式。為簡化建模過程，對列與列之間的EFP連接進行了斷開處理。

由圖3可以看出，在爆轟產物和爆轟壓力的作用下，大約在10 μs時罩頂微元開始被壓垮變形并流向藥形罩的對稱中心，這個時候藥形罩也同時被壓垮和變形，于是整個藥形罩就開始向前運動（圖3（b））。在大約20 μs時，藥形罩微元在對稱中心發生堆積現象，各微元之間相互碰撞、擠壓，就使得罩壁厚不斷增加，于是藥形罩微元中速度高的部分就會產生“射流”，而速度低的部分就會產生“杵體”（圖3（c））。在大約43 μs時，MEFP基本成型，所形成的破片具有比較均勻的質量分布和較高的速度。同時破片場的分布也很均勻（圖3（d））。

由圖4可知方案3形成的破片質量為15.511 g，遠大于一般的預制破片，同時其速度也達到1 658.9 m/s。下面把各方案EFP的最終成型圖展示如下（圖5）。

圖5 各仿真方案的最終成型圖

3 結果分析及優化方案確定

由仿真最終成型圖可以看出，破片近似與一個橢圓柱，其平均迎風面積近似為其表面積的四分之一。得到各方案破片尺寸及平均迎風面積見表3。

表3 各方案破片尺寸及平均迎風面積

根據指標要求有效殺傷半徑為20 m，由速度衰減公式可得各方案破片飛行20 m后的剩余速度vf，剩余動能Ef，剩余比動能Es的值如表4所示。

表4 各方案仿真計算結果

對以上各值進行極差分析，所得的結果如表5、表6和表7所示。

表5 質量m的極差分析結果

表6 剩余速度vf的極差分析結果

表7 剩余比動能Es的極差分析結果

根據表5、表6和表7，可以得到各指標隨各因素的變化情況如圖6所示，由此可以得到藥形罩結構參數對戰斗部指標的影響趨勢，從而得到一組最優的結構參數。

圖6 各方案指標隨各因素的變化圖

由計算結果可以看出，對于各指標，其重要性依次為:剩余比動能Es最重要，剩余速度vf次之，質量m最小。各指標隨因素的變化情況如圖6所示，根據極差分析可知:

罩頂厚δ對剩余比動能Es和質量m影響最大，δ越小Es越大，δ宜取2.2 mm;δ越大 m 越大，δ宜取 2.4 mm;罩頂厚δ對剩余速度vf是第二要指標，δ越小vf越大，δ宜取2.2 mm;故罩頂厚δ取2.2 mm最好。

同理可得，內表面曲率半徑r1取28 mm時各項指標最好，外表面曲率半徑r2取22 mm時各項指標最好，裝藥類型N為8701炸藥時各項指標最好。

綜上分析，MEFP戰斗部優化方案是;罩頂厚δ=2.2 mm，內表面曲率半徑r1=28 mm，外表面曲率半徑r2=22 mm，裝藥類型N為8701炸藥。

4 小結

在MEFP成型的過程中，藥形罩的裝藥材料、曲率半徑、材料和壁厚等結構參數都對MEFP的成型都產生著影響，所以有必要對以上因素對MEFP成型的影響進行正交試驗分析，把藥形罩的罩頂厚、內表面曲率半徑、外表面曲率半徑和裝藥類型確定為正交優化設計的4個因素，用L9（34）正交表來安排4個因素所對應的9個水平，得到9個仿真方案，將剩余比動能、剩余速度和質量作為優化設計評定指標。通過極差分析得到一種相對最優方案是:罩頂厚δ=2.2 mm，內表面曲率半徑r1=28 mm，外表面曲率半徑r2=22 mm，裝藥類型N為8701炸藥。

［1］曹兵.EFP成型機理及關鍵技術研究［D］.南京:南京航空航天大學博士后士學位論文，2001，2.

［2］李潤蔚.藥型罩及殼體參數影響EFP成型的數值計算研究［J］.彈箭與制導學報，1996，（3）:42-46.

［3］曹兵.多點起爆下裝藥爆轟及 EFP成型研究［D］.南京:南京理工大學博士論文，1998，6.

［4］龐勇，于川，桂毓林.球缺藥型罩爆炸成型彈丸數值模擬［J］.高壓物理學報，2005，（19）1:86-92.

［5］秦友花，周聽清.爆炸成型彈丸的試驗研究［J］.實驗力學，2002，17（2）:160-163.