梁振剛，蔣建偉

(1 北京理工大學爆炸科學與技術國家重點實驗室，北京 100081;2 沈陽理工大學裝備工程學院，沈陽 110159)

0 引言

周向MEFP(multiple explosively formed projectile)戰(zhàn)斗部爆炸后同時形成多個EFP彈丸，實現(xiàn)了EFP單點毀傷向多點毀傷的飛躍，可用于應對高空目標及輕型裝甲目標。

目前國內對 MEFP 的研究很多［1－3］，但對周向MEFP成型的研究較少。Richard Fong等設計了新型的周向聚焦MEFP戰(zhàn)斗部，可根據(jù)目標的不同，設計不同的裝藥結構，形成不同速度和形狀的MEFP［4］;張利等設計了一種新型多枚爆炸成型彈丸戰(zhàn)斗部，為小口徑MEFP戰(zhàn)斗部設計提供了參考［5］。文中設計了周向MEFP戰(zhàn)斗部的仿真模型，仿真分析了在不同藥型罩結構參數(shù)和戰(zhàn)斗部結構參數(shù)下EFP成型的特點，總結了藥型罩參數(shù)和戰(zhàn)斗部結構參數(shù)對周向MEFP成型規(guī)律的影響。

1 周向MEFP戰(zhàn)斗部及數(shù)值模擬

1.1 周向MEFP戰(zhàn)斗部模型及材料參數(shù)的選取

1)周向MEFP戰(zhàn)斗部模型建立

本周向MEFP戰(zhàn)斗部是外表面為棱柱形殼體并分層均勻排列若干個藥型罩，共設置了m×n個藥型罩，其中m為戰(zhàn)斗部軸向藥型罩分層層數(shù)，n為每層藥型罩周向分布數(shù)量。文中選用球缺型、等罩厚藥型罩進行研究，整個戰(zhàn)斗部的結構如圖1所示。

2)主要材料的選取

周向 MEFP選取 B炸藥或 OCTOL炸藥，采用JWL狀態(tài)方程，材料參數(shù)如表1。藥型罩的材料選用銅，殼體選用鋼4340，兩個材料均選用Linear狀態(tài)方程和Johnson-Cook強度模型。

圖1 周向MEFP結構圖

表1 炸藥參數(shù)

1.2 仿真模型的建立

圖2 全模型的1/32仿真模型

使用顯式動力學分析軟件AUTODYN-3D建立的模型，其中內部裝藥、端蓋和藥型罩均采用流固耦合的算法來滿足爆炸成型過程中大變形的要求。為了獲得較好的EFP成型速度與形狀，文中選用沿軸向中心線起爆方式［6］。周向MEFP模型較大，根據(jù)周向MEFP裝藥結構上下對稱的原則，可將全模型分成1/2模型，再根據(jù)藥型罩沿周向呈均勻分布的特點，設計出了全模型的1/(4×8)仿真模型(見圖2)，該裝藥側端面施加了位移約束。同時為了進一步降低計算量，提高計算速度，在不考慮上下端面對罩的稀疏波影響的前提下，根據(jù)各層EFP成型速度和形狀接近的規(guī)律及裝藥結構的對稱性，選用單個藥型罩裝藥仿真模型代替多層藥型罩仿真模型的方法，設計了簡化后的藥型罩仿真模型，并在裝藥側端面和上下端面均施加了位移約束(見圖3)。為了檢驗簡化后的仿真模型成型效果，同時數(shù)值模擬了上述兩種仿真模型的成型過程，模擬方案為:中心線起爆，炸藥選用B炸藥，裝藥外接圓柱直徑92 mm，殼體厚度3 mm，藥型罩軸向間距3 mm，藥型罩分為10層，每層周向分布8個罩，藥型罩口徑32 mm，藥型罩曲率半徑44 mm，藥型罩厚度3 mm。

圖3 簡化后的藥型罩仿真模型

1.3 仿真計算結果

兩種仿真模型的成型圖如圖4、圖5所示。

圖4 全模型的1/32仿真模型成型圖

圖5 簡化后的EFP成型圖

從圖4可看出，中心線起爆過程中，炸藥的爆轟波的傳播比較均衡，藥型罩沿中心軸線向外飛散同時會受到徑向擠壓和軸向拉伸，形成密實的壓垮式EFP［7］。從圖4和圖5可見，兩種仿真模型的成型形狀基本相同。

以靠近上下對稱面的藥型罩為1號，依次向上為2、3、4、5 號罩，各罩成型速度見表2。

表2 各層EFP成型速度統(tǒng)計

簡化仿真模型的EFP成型速度是1 970.1 m/s，由表2比較可知，1/32全模型的EFP成型速度與簡化模型的成型速度接近。仿真結果說明，簡化后的藥型罩仿真模型的成型速度、周向MEFP形狀與全模型的接近，所以在不考慮上下端蓋對藥型罩的稀疏波影響時，可以用簡化后的藥型罩仿真模型代替全模型進行周向MEFP成型規(guī)律的仿真研究。

2 藥型罩參數(shù)對周向MEFP成型影響

藥型罩各項參數(shù)中直接影響EFP成型性能的重要因素有藥型罩的罩口徑、曲率半徑、藥型罩的厚度，所以在研究周向MEFP成型規(guī)律時，研究這三個參數(shù)對MEFP成型影響是很必要的。本研究中選取簡化后的單個藥型罩仿真模型進行計算，周向MEFP戰(zhàn)斗部的殼體厚度C是3 mm，炸藥選用B炸藥，裝藥外接圓柱直徑D為92 mm，藥型罩軸向間距為3 mm，周向分布藥型罩數(shù)量為8個，罩口徑為d，EFP成型速度為V，長徑比為E，罩厚度為h，罩曲率半徑為r。

2.1 罩口徑與裝藥外接圓柱直徑比值對周向MEFP成型影響

在研究罩口徑d與裝藥外接圓柱直徑D比值對周向MEFP成型影響時，藥型罩的厚度設為3 mm，曲率半徑為40 mm，藥型罩的罩口徑與裝藥外接圓柱直徑比值分別選取 0.3、0.312、0.325、0.337 和 0.35，周向MEFP的成型速度V及長徑比隨d/D變化曲線結果如圖6和圖7所示。

圖6 EFP成型速度隨d/D比值變化曲線

圖7 EFP長徑比隨d/D比值變化曲線

由圖6可知隨著罩口徑與裝藥外接圓柱直徑比值的增大，EFP成型速度雖然略有降低，但下降幅度不大。由圖7可知，隨著罩口徑與裝藥外接圓柱直徑比值的加大，長徑比略增加，但幅度也不明顯，說明飛行中藥型罩受到徑向擠壓和軸向拉伸幅度變化不大。所以調整藥型罩口徑參數(shù)對MEFP成型影響不大。

2.2 罩厚度與裝藥外接圓柱直徑比值對周向MEFP裝藥成型影響

在研究罩厚度h與裝藥外接圓柱直徑D比值對周向MEFP成型影響時，設計的仿真方案數(shù)據(jù)如下:藥型罩的曲率半徑為40 mm，藥型罩的口徑為29.4 mm，罩厚度與裝藥外接圓柱直徑比值分別為0.023 8、0.026 8、0.029 8、0.032 7 和 0.035 7，周向MEFP成型速度V及長徑比隨h/D變化曲線如圖8和圖9所示。

圖8 EFP成型速度隨h/D比值變化曲線

圖9 長徑比隨h/D比值變化曲線

由圖8可看出:h/D值越大，EFP成型速度越慢，而且變化明顯。這主要因為隨著厚度的增大，藥型罩承受爆轟波的能力增強，藥型罩的質量也隨之增大所致。圖9的曲線顯示隨著罩厚度與裝藥外接圓柱直徑比值的增大，成型過程中徑向擠壓和軸向拉伸的幅度也降低，最終EFP的長徑比也明顯下降。

2.3 藥型罩曲率半徑與裝藥外接圓柱直徑比值對周向MEFP裝藥成型影響

在研究罩曲率半徑與裝藥外接圓柱直徑比值對周向MEFP成型影響時，設計的仿真方案數(shù)據(jù)如下:藥型罩的厚度為3 mm，藥型罩的口徑為32 mm，藥型罩曲率半徑r與裝藥外接圓柱直徑D比值分別為0.347 8、0.391 3、0.434 8、0.478 3 和 0.521 7，周向MEFP的成型速度V及長徑比隨罩r/D比值變化曲線如圖10和圖11所示。

圖10 EFP成型速度隨r/D比值變化曲線

圖11 EFP長徑比隨r/D比值變化曲線

由圖10可知，r/D比值的增大，對EFP成型速度影響不大。圖11顯示隨著r/D比值的增大，EFP的長徑比明顯減小，說明藥型罩受到徑向擠壓力增大，而軸向拉伸幅度變化不大，從而使長徑比降低。

3 戰(zhàn)斗部結構參數(shù)對周向MEFP成型影響

3.1 殼體壁厚對周向MEFP成型影響

周向MEFP戰(zhàn)斗部殼體壁厚對周向MEFP成型影響的仿真方案如下，裝藥外接圓柱直徑92 mm，藥型罩的口徑32 mm，藥型罩厚度3 mm，藥型罩曲率半徑36 mm，殼體壁厚分別選取2 mm、3 mm、4 mm和5 mm，仿真結果如圖12和圖13所示。

圖12 EFP成型速度隨殼體壁厚變化曲線

由圖12和圖13可看出殼體越厚，EFP成型速度越快，EFP的長徑比也隨之增加且增幅明顯。這是由于殼體厚度增大，將引起藥型罩變形過程中各微元向軸線并攏的速度加大。

圖13 EFP長徑比隨殼體厚度變化曲線

3.2 周向藥型罩分布數(shù)量對周向MEFP成型影響

在研究周向分布數(shù)量對周向MEFP成型影響時，炸藥選取B炸藥，裝藥外接圓柱直徑92 mm，藥型罩的口徑27.6 mm，藥型罩厚度3 mm，藥型罩曲率半徑40 mm，殼體壁厚3 mm，周向MEFP戰(zhàn)斗部側面每層藥型罩分布數(shù)量依次為8個、9個和10個，仿真結果如圖14和15所示。

圖14 EFP成型速度隨周向數(shù)量變化曲線

圖15 EFP長徑比隨周向數(shù)量變化曲線

由圖14和圖15可知，在藥量一定的情況下，隨著周向藥型罩分布數(shù)量的增加，EFP的成型速度及長徑比都降低，但幅度很小。

3.3 兩種不同炸藥對周向MEFP成型影響

仿真方案為裝藥外接圓柱直徑92 mm，藥型罩的口徑32 mm，藥型罩厚度2.5 mm，藥型罩曲率半徑36 mm，殼體壁厚3 mm，徑向4層藥型罩，每層分布8個藥型罩，炸藥分別選取B炸藥和OCTOL炸藥，仿真結果如表3。

表3 選用兩種炸藥的EFP成型數(shù)據(jù)

因為OCTOL炸藥的能量密度大于B炸藥，所以形成EFP成型速度也會更快一些，長徑比也相應增大。

4 結論

1)在不考慮端面稀疏波影響的前提下，文中設計的周向MEFP簡化仿真模型與全模型仿真模型的EFP成型形狀一致，速度接近，可以代替全模型仿真模型進行周向MEFP成型規(guī)律的研究。

2)在裝藥結構一定的情況下，調整藥型罩的口徑和厚度可調整EFP的動能，調整藥型罩的曲率半徑可調整EFP長徑比。

3)隨著殼體厚度的增加，EFP成型速度小幅加快，長徑比增幅比較明顯;周向藥型罩數(shù)量的分布對EFP成型速度和長徑比的影響均很小;高能炸藥會使EFP成型速度和長徑比均增加，其中長徑比增幅較明顯。

以上結果可為周向MEFP戰(zhàn)斗部結構設計提供參考。

［1］孔毓琦，郭銳.MEFP式戰(zhàn)斗部正交試驗優(yōu)化［J］.機械制造與自動化，2013(1):22－26.

［2］陳宇，曹紅松，王智軍，等.一種新型防空MEFP的設計與仿真［J］.現(xiàn)代防御技術，2013，41(2):201－206.

［3］吳峰，吳峰，孔毓琦，等.MEFP式戰(zhàn)斗部侵徹鋼靶的毀傷效應研究［J］.機械制造與自動化，2013(2):87－90.

［4］Cox M，Campbell A W.Corner-turning in TATB［C］∥Seventh Symposium(International)on Detonation.Dahlgren:Office of Naval Research，1981:624－633.

［5］張利，彭雪琴.周向多爆炸成型彈丸戰(zhàn)斗部仿真［J］.四川兵工學報，2014(9):32－35.

［6］王樹有，蔣建偉，門建兵.準球形EFP成型影響因素的數(shù)值模擬［J］.含能材料，2008，16(6):731－733.

［7］Joseph Carleone.Tactical Missle Warheads［M］.1993.