基于LS-DYNA的大口徑EFP數(shù)值仿真

孫憲福，任 邯，姜 濤

(1. 92815 部隊(duì)，浙江 寧波 315717;2.海軍潛艇學(xué)院，山東 青島 266071)

爆炸成型彈丸(EFP)的形狀和速度與很多因素有關(guān)，如藥型罩的材料、藥型罩的曲率半徑、藥型罩的厚度、裝藥高度、裝藥有無(wú)殼體以及殼體的材料與厚度等［1］。本文僅對(duì)藥型罩曲率半徑、藥型罩厚度和裝藥高度3 個(gè)參數(shù)進(jìn)行數(shù)值計(jì)算，以得出其對(duì)EFP 性能的影響。與已有關(guān)于EFP 的實(shí)驗(yàn)與數(shù)值計(jì)算所不同的是，本文所模擬的均是大口徑EFP(500mm 口徑)的形成與性能參數(shù)。LS-DYNA(nonlinear dynamic analysis of structures in three dimension)，即三維動(dòng)態(tài)非線性結(jié)構(gòu)分析，是世界上最著名的通用顯式動(dòng)力分析程序，能模擬真實(shí)世界的各種復(fù)雜問題，特別適于求解各種二維、三維非線性結(jié)構(gòu)高速碰撞、爆炸和金屬成型等非線性動(dòng)力沖擊問題，同時(shí)也可求解傳熱、流體及流固耦合問題。

1 計(jì)算模型

1.1 控制方程

LS-DYNA 程序的主要算法采用 Lagrange 描述增量法［2］。

取初始時(shí)刻質(zhì)點(diǎn)坐標(biāo)為Xi(i=1，2，3)，在任意時(shí)刻t，該質(zhì)點(diǎn)的坐標(biāo)為xi(i=1，2，3)，則該質(zhì)點(diǎn)的運(yùn)動(dòng)方程為:

在t=0 時(shí)，初始條件為:

式中Vi為初始速度。

動(dòng)量方程:

式中:σij為柯西應(yīng)力;fi為單位質(zhì)量體積力;為加速度。

質(zhì)量守恒方程:

式中:ρ 為當(dāng)前質(zhì)量密度;ρ0為初始質(zhì)量密度;J 為雅克比行列式。

能量方程:

式中:V 為現(xiàn)時(shí)構(gòu)形的體積;˙εij為應(yīng)變率張量;q 為體積黏性阻力。

偏應(yīng)力:

壓力:

1.2 模型的參數(shù)

模型由2 部分組成:高能炸藥和藥型罩，二者均采用拉格朗日算法?？紤]到計(jì)算時(shí)間問題，建立1/4 模型，如圖1所示。其中高能炸藥劃分38 880個(gè)單元，藥型罩劃分3 880個(gè)單元。

圖1 EFP 的有限元模型

炸藥單元材料模型采用JWL 狀態(tài)方程描述:

紫銅材料模型采用GRUNEISEN 狀態(tài)方程描述［4］，材料沖擊壓縮時(shí)其表達(dá)式為:

材料膨脹時(shí)其表達(dá)式為:

2 數(shù)值計(jì)算結(jié)果

2.1 藥型罩曲率半徑的影響

藥型罩曲率半徑對(duì)EFP 形狀的影響如圖2 所示。

圖2 各曲率半徑對(duì)EFP 形狀的影響

從計(jì)算結(jié)果可以看出，藥型罩的曲率半徑越大，所產(chǎn)生彈丸越短粗，其質(zhì)心越靠后;藥型罩曲率半徑越小，所產(chǎn)生彈丸越細(xì)長(zhǎng)，其質(zhì)心越靠前。但是，當(dāng)藥型罩曲率半徑為30 cm時(shí)，彈丸形成過程中發(fā)生了斷裂。

藥型罩的曲率半徑對(duì)彈丸的速度有著同樣重要的影響［3］。表1 為各曲率半徑藥型罩結(jié)構(gòu)下彈丸的頭部速度。

表1 各曲率半徑藥型罩結(jié)構(gòu)下彈丸頭部速度

從速度計(jì)算結(jié)果可知:藥型罩曲率半徑越小，彈丸頭部獲得的最大速度越大，但是彈丸尾部獲得的速度差別卻不明顯，導(dǎo)致彈丸形狀有差異。即藥型罩曲率半徑小，由于軸向速度的差別大，將彈丸拉長(zhǎng)，彈丸就細(xì)長(zhǎng);藥型罩曲率半徑大，由于軸向速度差別小，彈丸被拉長(zhǎng)的幅度小，故彈丸就短粗［4］。

2.2 藥型罩厚度對(duì)彈丸速度的影響

藥型罩厚度對(duì)彈丸的影響主要體現(xiàn)在對(duì)彈丸速度的影響上。分別對(duì)裝藥高度20 cm，藥型罩曲率半徑50 cm，藥型罩口徑50 cm 條件下，藥型罩壁厚分別為0.5，1，2 cm 的情況作了數(shù)值計(jì)算。具體結(jié)果見圖3。

圖3 藥型罩厚度對(duì)彈丸頭部速度影響

從計(jì)算結(jié)果可以看出:藥型罩壁厚越厚，其獲得的速度越小。但是壁厚薄引起的彈丸形狀的改變、彈丸動(dòng)能的減少也不容忽視。當(dāng)罩的壁厚過小時(shí)，難以形成彈丸。

2.3 裝藥高度對(duì)彈丸速度的影響

裝藥高度的改變可影響爆轟波作用藥型罩姿態(tài)，也影響藥型罩獲得的能量。為了分析裝藥高度對(duì)彈丸的影響，在藥型罩曲率半徑50 cm，藥型罩壁厚1 cm，藥型罩口徑50 cm 條件下，分別對(duì)裝藥高度100、200、300 mm 的情況作了數(shù)值計(jì)算，結(jié)果見表2。

表2 各裝藥高度條件下彈丸頭部速度

數(shù)值計(jì)算結(jié)果表明:裝藥高度越大，藥型罩獲得的能量也越多，其速度也越大。但是當(dāng)裝藥高度增長(zhǎng)到一定值時(shí)，裝藥繼續(xù)增加對(duì)彈丸的速度將不會(huì)有影響或僅有很小影響。

3 結(jié)束語(yǔ)

由數(shù)值計(jì)算結(jié)果可以看出大口徑EFP 在各曲率半徑下的彈丸形狀。在曲率半徑與EFP 口徑比為1 左右時(shí)，形成的彈丸較密實(shí)，彈丸的長(zhǎng)度適中;不同于小口徑的EFP，大口徑EFP 在藥型罩厚度為5 ～10 mm(藥型罩厚度與EFP 口徑比為0.0167 ～0.033)時(shí)，形成的彈丸速度較大且形狀較好;在裝藥高度不小于200 mm 時(shí)，有利于彈丸獲得較高速度。

［1］李裕春，楊萬(wàn)江，沈蔚.藥型罩曲率半徑對(duì)爆炸成型彈丸的影響［J］.火工品，2003(1): 45-48.

［2］時(shí)黨勇，李裕春，張勝民.基于ANSYS/LS-DYNA8.1 進(jìn)行顯示動(dòng)力分析［M］.北京:清華大學(xué)出版社，2005.

［3］王成，琿壽榕，黃風(fēng)雷.同口徑破-破型串聯(lián)裝藥戰(zhàn)斗部的實(shí)驗(yàn)研究［J］.彈箭與制導(dǎo)學(xué)報(bào)，2003，22(2):31-34.

［4］謝文，龍淵，岳小兵，等.模擬爆炸成形彈丸對(duì)大間隔靶的侵徹試驗(yàn)［J］.兵工學(xué)報(bào)，2003，24(3):392-395.

［5］唐平江，顧文彬，羅智偉. 變壁厚殼體對(duì)形成帶尾翼EFP 的影響［J］.四川兵工學(xué)報(bào)，2009(1):38-40.

［6］朱峰，朱衛(wèi)華，王怡舒. 彈丸初速度對(duì)侵徹雙層鋼板影響的數(shù)值模擬［J］.四川兵工學(xué)報(bào)，2010(9):28-30.