苗潤(rùn)源，程淑杰，梁安定，梁爭(zhēng)峰

(西安近代化學(xué)研究所， 西安 710065)

爆炸成型彈丸(explosively formed penetrators,EFP)具有更高的打擊比動(dòng)能以及良好的氣動(dòng)外形，大大提高了穿甲能力，可直接“擊爆”來(lái)襲目標(biāo)的戰(zhàn)斗部，對(duì)于反厚壁類精確制導(dǎo)彈藥具有重要的意義[1-2]。周向MLEFP具有傳統(tǒng)EFP的優(yōu)點(diǎn)外，同時(shí)也具有彈丸空間分布密度高的特點(diǎn)，成為近年來(lái)研究的熱點(diǎn)[3-5]。目前關(guān)于周向MLEFP的研究主要集中在總結(jié)藥型罩曲率半徑、壁厚、殼體等宏觀結(jié)構(gòu)參數(shù)以及起爆方式等對(duì)周向MLEFP成型的影響規(guī)律，并且藥型罩結(jié)構(gòu)大多采用球缺型藥型罩，而對(duì)于大錐角藥型罩微元壓垮速度的變化對(duì)周向MLEFP的成型影響規(guī)律研究較少[6-9]?；谏鲜鲈?，本文采用ANSYS/LS-DYNA顯式動(dòng)力分析軟件對(duì)多個(gè)具有不同大錐角藥型罩結(jié)構(gòu)參數(shù)的周向MLEFP的成型結(jié)果進(jìn)行了數(shù)值模擬，并通過(guò)將彈丸成型結(jié)果與藥型罩微元壓垮速度進(jìn)行結(jié)合對(duì)比分析，得到了藥型罩微元壓垮速度對(duì)周向MLEFP成型性能的影響規(guī)律[10-11]，并通過(guò)戰(zhàn)斗部試驗(yàn)對(duì)所得到的結(jié)論進(jìn)行驗(yàn)證。該研究對(duì)于以后周向MLEFP戰(zhàn)斗部的設(shè)計(jì)可提供一定的指導(dǎo)。

1 數(shù)值模擬模型

1.1 模型組成及結(jié)構(gòu)參數(shù)

建立如圖1所示的周向MLEFP數(shù)值模擬計(jì)算模型，在殼體圓周方向上通過(guò)預(yù)制刻槽的方式均勻分布18個(gè)藥型罩，每個(gè)藥型罩的結(jié)構(gòu)參數(shù)均相同[12-13]，建模所采用的單位制為cm-g-μs。為減小數(shù)值模擬計(jì)算量，建立了1/4模型。

整個(gè)計(jì)算模型共分為3個(gè)PART，PART1為炸藥、PART2為空氣、PART3為藥型罩，其中炸藥材料選用高能混合炸藥，藥型罩材料選用純鐵。

該戰(zhàn)斗部的計(jì)算模型網(wǎng)絡(luò)總數(shù)：part1的網(wǎng)格總數(shù)為406 620個(gè)，part2的網(wǎng)格總數(shù)為302 670個(gè)，part3的網(wǎng)格總數(shù)為12 096個(gè)。

圖1 周向MLEFP數(shù)值模擬計(jì)算模型示意圖

圖2為藥型罩結(jié)構(gòu)參數(shù)示意圖。其中D為戰(zhàn)斗部外徑，H為藥型罩高度，K為藥型罩頂點(diǎn)裝藥高度，δ1、δ2、R1、R2分別為藥型罩頂點(diǎn)厚度、邊緣厚度、內(nèi)徑和外徑。

圖2 藥型罩結(jié)構(gòu)示意圖

為研究藥型罩壓垮速度的變化對(duì)彈丸成型性能的影響，依次改變藥型罩曲率半徑、頂點(diǎn)及邊緣厚度比值、頂點(diǎn)裝填比，設(shè)計(jì)表1所示的4個(gè)方案藥型罩參數(shù)。所有藥型罩參數(shù)方案保持戰(zhàn)斗部外徑D、藥型罩高度H大小不變，并且H/D=0.5。

1.2 材料參數(shù)

使用顯式動(dòng)力分析有限元程序ANSYS/LS-DYNA進(jìn)行數(shù)值計(jì)算，為避免單元發(fā)生畸變影響計(jì)算過(guò)程，炸藥及空氣采用ALE算法，藥型罩采用拉格朗日算法，并通過(guò)關(guān)鍵字CONSTRINED_LAGRANGE_IN-SOLID進(jìn)行流固耦合[14-15]。

表1 藥型罩參數(shù)

藥型罩本構(gòu)模型采用John-Cook本構(gòu)模型，狀態(tài)方程為Gruneisen狀態(tài)方程，Gruneisen狀態(tài)方程表達(dá)式為：

(γ0+aμ)E

(1)

式(1)中:P為壓力;ρ0為初始密度；C為VS-VP曲線截距；S1、S2、S3為VS-VP曲線斜率系數(shù)；γ0為Gruneisenγ值；a為一階修正量。 對(duì)于純鐵材料，初始密度為7.85 g/cm3，剪切模量為0.77。

主裝藥本構(gòu)模型為HIGH_EXPLOSIVE_BURN，狀態(tài)方程為JWL狀態(tài)方程，并且JWL狀態(tài)方程表達(dá)式為：

(2)

表2為JWL狀態(tài)方程相關(guān)參數(shù)值。

表2 JWL狀態(tài)方程參數(shù)

主裝藥采用的本構(gòu)模型中，裝藥密度為1.78 g/cm3，爆速為8 200 m/s，CJ壓力為0.372 Mbar。

空氣采用EOS-LINEAR-POLYNOMIAL線性多項(xiàng)式狀態(tài)方程和MAT-NULL空材料模型，并且空氣密度為1.293×10-3g/cm3。

為保證彈丸成型過(guò)程中爆轟波的對(duì)稱性，起爆方式為戰(zhàn)斗部中心點(diǎn)起爆。

2 數(shù)值模擬結(jié)果及分析

2.1 數(shù)值模擬結(jié)果

對(duì)表1中4種采用不同藥型罩結(jié)構(gòu)參數(shù)的周向MLEFP成型結(jié)果進(jìn)行模擬，圖3(a)～圖3(d)分別為表1中4種彈丸最終成型結(jié)果。

由圖3彈丸成型結(jié)果可以看出，方案1中彈丸并未成型；方案2相比于方案1減小了曲率半徑，中彈丸頭部壓實(shí)，但長(zhǎng)徑比較??；方案3相比于方案2增加了藥型罩邊緣厚度后，彈丸長(zhǎng)徑比增加，成型效果較理想；方案4減小了藥型罩頂點(diǎn)裝填比，彈丸速度相比于方案3彈丸速度下降了15%，彈丸未成型。

圖3 周向MLEFP計(jì)算模型示意圖

2.2 數(shù)值模擬結(jié)果分析

分別選取每個(gè)方案中從藥型罩頂點(diǎn)到邊緣的10個(gè)微元，并得到每個(gè)微元壓垮速度的軸向及徑向分量。圖4、圖5分別為每個(gè)方案中藥型罩微元壓垮速度的軸向及徑向分量隨微元位置變化的分布曲線。

圖4 藥型罩微元壓垮速度軸向分量分布曲線

圖5 藥型罩微元壓垮速度徑向分量分布曲線

由圖4可以看出，方案1和2中藥型罩頂點(diǎn)及邊緣壓垮速度軸向差值均為180 m/s，彈丸成型結(jié)果均不理想；方案3中藥型罩頂點(diǎn)與邊緣的壓垮速度軸向差值達(dá)到了220 m/s，大于方案1和方案2中的藥型罩頂點(diǎn)與邊緣的軸向速度差，形成的彈丸長(zhǎng)徑比最大，彈丸成型最理想；方案4中藥型罩頂點(diǎn)與邊緣壓垮速度軸向差值達(dá)到210 m/s，但由于藥型罩頂點(diǎn)壓垮速度軸向值相比于方案3，由1 900 m/s下降到1 650 m/s彈丸未成型。

由圖5可以看出，方案2中藥型罩微元壓垮速度徑向差值大于方案1中的差值，形成的彈丸頭部壓實(shí)，但長(zhǎng)徑比太小。方案4藥型罩微元壓垮速度徑向差值最大，但由于藥型罩頂點(diǎn)壓垮速度軸向值太小，彈丸未成型。

以上數(shù)值模擬結(jié)果分析表明，藥型罩微元壓垮速度梯度與表1方案3中藥型罩微元壓垮速度梯度相同時(shí)，彈丸成型性能最好。

3 戰(zhàn)斗部試驗(yàn)驗(yàn)證

依據(jù)數(shù)值模擬所得到的結(jié)論，選取表1方案3中的參數(shù)作為周向MLEFP戰(zhàn)斗部藥型罩設(shè)計(jì)參數(shù)，設(shè)計(jì)如圖6所示的三層周向MLEFP戰(zhàn)斗部進(jìn)行戰(zhàn)斗部毀傷性能試驗(yàn)，藥型罩的材料為純鐵。

圖6 戰(zhàn)斗部結(jié)構(gòu)示意圖

圖6中戰(zhàn)斗部藥型罩的材料為純鐵，主裝藥為高能熔鑄炸藥。

圖7為戰(zhàn)斗部靜爆試驗(yàn)場(chǎng)布局照片，戰(zhàn)斗部豎直放置，距離戰(zhàn)斗部中心 6 m半徑處布置20 mm厚的Q235鋼板。引爆戰(zhàn)斗部后，通過(guò)統(tǒng)計(jì)彈丸對(duì)20 mm厚的Q235鋼板的穿甲率，對(duì)表1中F3彈丸成型性能進(jìn)行驗(yàn)證。

圖7 戰(zhàn)斗部靜爆試驗(yàn)場(chǎng)布局照片

圖8為引爆戰(zhàn)斗部后，周向MLEFP在目標(biāo)距離處的鋼板的穿孔分布圖。

圖8 周向MLEFP鋼板穿孔照片

依據(jù)試驗(yàn)結(jié)果統(tǒng)計(jì)，周向MLEFP對(duì)6 m半徑處20 mm厚Q235鋼板的穿甲率為100%；靶板上彈丸穿孔的形狀大小基本與方案3中通過(guò)數(shù)值模擬得到的彈丸形狀一致，均為扁平狀。

以上試驗(yàn)結(jié)果表明：當(dāng)周向MLEFP初速不低于1 800 m/s且藥型罩頂點(diǎn)與邊緣壓跨速度軸向速度差值達(dá)到220 m/s時(shí)，周向MLEFP能夠成型性能良好，并能穩(wěn)定飛行，具有很好的毀傷效果。

4 結(jié)論

1) 藥型罩頂點(diǎn)與邊緣壓垮速度軸向差值是彈丸成型的主要影響因素；增加藥型罩微元壓垮速度徑向差值有利于彈丸的壓合，但不利于增加彈丸的長(zhǎng)徑比。

2) 對(duì)于藥型罩高度與戰(zhàn)斗部最大直徑比值為0.5的周向MLEFP，初速不低于1800 m/s時(shí)并且藥型罩頂點(diǎn)與邊緣的壓垮速度軸向差值達(dá)到220 m/s時(shí)，形成的彈丸長(zhǎng)徑比大，能夠穩(wěn)定飛行，并具有很好的毀傷效應(yīng)。

3) 藥型罩頂點(diǎn)壓垮速度軸向值下降過(guò)多時(shí)，增加周向MLEFP成型時(shí)頂點(diǎn)與邊緣壓垮速度軸向差值有利于提高彈丸的成型性能。