葉 嚴(yán)，姚志敏，李金明，劉 波

?

大口徑EFP侵徹典型裝甲鋼板過程的數(shù)值仿真與試驗(yàn)研究

葉 嚴(yán)1，姚志敏1，李金明2，劉 波1

(1.軍械工程學(xué)院導(dǎo)彈工程系，河北 石家莊，050003；2.軍械工程學(xué)院彈藥工程系，河北 石家莊，050003）

通過大口徑EFP高速侵徹裝甲靶板試驗(yàn)，以及利用AUTODYN-3D有限元仿真軟件對(duì)整個(gè)侵徹過程進(jìn)行數(shù)值模擬，研究了EFP開坑、穩(wěn)定侵徹、尾翼侵徹和沖塞貫穿形成二次破片的物理過程，模擬結(jié)果與試驗(yàn)現(xiàn)象和理論分析均吻合較好，并從原理上分析了試驗(yàn)中各宏觀現(xiàn)象產(chǎn)生的原因。研究結(jié)果不但認(rèn)識(shí)了EFP侵徹裝甲靶板的機(jī)理，也可以為聚能裝藥對(duì)典型裝甲毀傷評(píng)估提供參考。

EFP；侵徹機(jī)理；裝甲靶板；數(shù)值仿真

爆炸成型彈丸(EFP)利用聚能原理，通過裝藥的爆轟作用，使高溫高壓的爆轟產(chǎn)物作用于金屬藥型罩上，使藥型罩材料發(fā)生極大的塑性變形而閉合形成具有較高質(zhì)心速度(1 500~3 000m/s)和一定結(jié)構(gòu)形狀的彈丸，利用動(dòng)能侵徹目標(biāo)[1]。目前，國(guó)內(nèi)外對(duì)EFP的成型及其影響因素和侵徹過程研究比較多[2-4]，但從原理上對(duì)整個(gè)侵徹過程進(jìn)行細(xì)致分析的還比較少[5]。由于EFP在大炸高下能保持完整的彈丸特性來攻擊裝甲目標(biāo)，其威力可達(dá)到在1 000倍裝藥口徑的距離上穿透厚度相當(dāng)于1倍裝藥口徑的均質(zhì)裝甲[6]，所以抗EFP沖擊是判斷坦克裝甲和混凝土工事防護(hù)能力的一個(gè)重要依據(jù)，研究EFP對(duì)裝甲靶板的高速侵徹效應(yīng)很有意義。

由于EFP高速侵徹裝甲靶板過程十分復(fù)雜，僅靠理論分析和實(shí)驗(yàn)很難認(rèn)識(shí)整個(gè)過程中所產(chǎn)生的侵徹效應(yīng)。為此，本文在試驗(yàn)呈現(xiàn)侵徹現(xiàn)象的同時(shí)，結(jié)合AUTODYN有限元仿真軟件[7]，采用光滑粒子流體動(dòng)力（SPH）方法計(jì)算，對(duì)EFP從形成到侵徹完成的整個(gè)過程進(jìn)行了數(shù)值模擬，并分析了侵徹過程中產(chǎn)生各現(xiàn)象的原因，從機(jī)理上研究了EFP對(duì)裝甲靶板的高速侵徹效應(yīng)。

1 試驗(yàn)布置與試驗(yàn)結(jié)果

1.1 EFP戰(zhàn)斗部結(jié)構(gòu)

圖1為EFP戰(zhàn)斗部結(jié)構(gòu)圖，EFP裝藥主要由藥型罩、裝藥、殼體和起爆裝置等組成，其中裝藥直徑D為140mm，藥型罩的內(nèi)曲率半徑1為144mm，外曲率半徑2為137mm，罩頂壁厚為6mm。

圖1 EFP戰(zhàn)斗部結(jié)構(gòu)圖

1.2 試驗(yàn)原理及布置

本文選用裝甲靶板作為侵徹對(duì)象，主要是因?yàn)閲?guó)內(nèi)外普遍采用此材料進(jìn)行EFP的威力考核，具有普遍意義。試驗(yàn)中，聚能裝藥選用大口徑140mm EFP試驗(yàn)彈，炸高為1 000mm，裝藥為8701炸藥，裝藥高度102mm，藥型罩材料為紫銅，變壁厚結(jié)構(gòu)。裝甲靶板的尺寸有兩種，一種為300mm×300mm×40mm；另一種為400mm×400mm×45mm。靶板由專用鋼制靶板支架固定，分別有0°、30°和45°的靶板支架（由于后效作用比較大，每做一次實(shí)驗(yàn)換一個(gè)靶板支架），置于回收水箱上，后效靶板置于水箱上方，距靶板1 000mm。在靶板上垂直放置一個(gè)高為1 000mm的炸高筒，并將EFP戰(zhàn)斗部放在炸高筒上，藥型罩口部朝下，正對(duì)靶板。校準(zhǔn)完成后，采用制式電雷管起爆。試驗(yàn)布置情況見圖2。

圖2 試驗(yàn)布置情況示意圖

1.3 試驗(yàn)結(jié)果

試驗(yàn)時(shí)對(duì)5發(fā)試驗(yàn)彈進(jìn)行編號(hào)，并記錄每發(fā)彈丸侵徹靶板的孔徑、出入口直徑、前后靶塑性變形區(qū)等數(shù)據(jù)。典型試驗(yàn)結(jié)果照片如圖3所示，試驗(yàn)數(shù)據(jù)見表1 。

圖3 典型試驗(yàn)結(jié)果照片

表1 試驗(yàn)方案與試驗(yàn)結(jié)果

如圖3和表1所示，5次試驗(yàn)均出現(xiàn)以下現(xiàn)象：靶板入口卷邊花瓣?duì)钇茐?；出口具有明顯拉伸斷裂特征的外翻花瓣形穿孔；入口直徑大于出口直徑；前靶出現(xiàn)塑性變形區(qū)；從垂直侵徹到斜侵徹，穿孔的形狀由圓形變成橢圓形，并且傾斜角度越大橢圓長(zhǎng)短半軸比越大。

2 EFP侵徹過程數(shù)值模擬

2.1 基本假設(shè)

通常情況下，對(duì)某一問題的仿真分析必須先做一些假設(shè)，將復(fù)雜目標(biāo)簡(jiǎn)化成可以建立數(shù)值仿真模型并保證目標(biāo)真實(shí)特性的簡(jiǎn)單目標(biāo)，這就要求對(duì)一特定的環(huán)境參數(shù)、幾何條件進(jìn)行簡(jiǎn)化，然后再運(yùn)用數(shù)學(xué)計(jì)算的方法來對(duì)目標(biāo)進(jìn)行量化計(jì)算模擬。采用Backman和Goldsmith對(duì)于侵徹的假設(shè)如下：（1）靶板的響應(yīng)集中在幾倍于彈體直徑的區(qū)域中，可忽略遠(yuǎn)離侵徹點(diǎn)的復(fù)雜邊界條件；（2）忽略包括摩擦在內(nèi)的熱效應(yīng)；（3）假定撞擊面局部為平面。

2.2 計(jì)算模型

由于聚能裝藥結(jié)構(gòu)比較復(fù)雜，在AUTODYN中無法對(duì)其進(jìn)行網(wǎng)格剖分，故采用LS-DYNA建立炸藥、藥型罩以及靶板的有限元網(wǎng)格；為降低計(jì)算量，充分考慮模型的對(duì)稱性，只建立其1/4模型，并生成文件；將LS-DYNA中生成的文件通過接口導(dǎo)入AUTODYN中，再將上述有限元網(wǎng)格通過SPH粒子替換，并賦予對(duì)應(yīng)材料，炸藥、藥型罩和靶板分別設(shè)置粒子間距為0.6mm、0.4mm和0.6mm，變光滑長(zhǎng)度，總粒子數(shù)為51 600個(gè)；EFP戰(zhàn)斗部及靶板SPH模型建立過程如圖4所示。

圖4 EFP戰(zhàn)斗部和靶板有限元模型

設(shè)置侵蝕準(zhǔn)則，并保留侵蝕單元節(jié)點(diǎn)的初始慣性運(yùn)動(dòng)特性，以此來模擬靶后破片的飛散效果。根據(jù)模擬方案，考慮到SPH方法計(jì)算效率偏低，且對(duì)計(jì)算機(jī)要求較高，數(shù)值計(jì)算采用配置為Intel Xeon E5-2690 CPU，64 GB內(nèi)存的工作站進(jìn)行。

2.3 材料模型與參數(shù)

該仿真采用AUTODYN-3D，采用的單位是cm-g-μs，各部分材料參數(shù)大部分來自AUTODYN標(biāo)準(zhǔn)材料庫，數(shù)值模擬中藥型罩為紫銅、靶板材料為裝甲鋼，材料模型如表2所示。

表2 材料模型

藥型罩及靶板材料模型均選用為Johnson-Cook 材料模型，并通過Gruneisen 狀態(tài)方程描述其變形過程；8701 炸藥采用JWL狀態(tài)方程描述。首先對(duì)EFP垂直侵徹40mm厚均質(zhì)裝甲靶板進(jìn)行了數(shù)值模擬，為驗(yàn)證材料模型、對(duì)應(yīng)參數(shù)以及數(shù)值模型的可靠性，將所得結(jié)果與試驗(yàn)所得結(jié)果進(jìn)行對(duì)比驗(yàn)證。

2.4 仿真試驗(yàn)結(jié)果分析

為直觀反映EFP侵徹靶板的整個(gè)過程，選取具有代表性的各階段的剖面圖，如圖5所示。

圖5 侵徹過程

由圖5可知，EFP對(duì)靶板的侵徹可分為4個(gè)階段：

（1）開坑階段(侵徹深度小于等于侵徹體直徑)：此時(shí)的碰撞速度最高，產(chǎn)生的碰撞應(yīng)力也最大，遠(yuǎn)遠(yuǎn)超過了EFP和靶板金屬材料的動(dòng)態(tài)強(qiáng)度，使靶板材料在碰撞的局部區(qū)域內(nèi)發(fā)生破壞、變形，并向抗力最小的方向飛濺排出，形成了靶前破片。這時(shí)EFP上只作用有慣性力和壓縮力，并不斷地破壞、飛濺。同時(shí)在靶板表面產(chǎn)生了不斷擴(kuò)大的彈坑，在靶內(nèi)建立起了相對(duì)穩(wěn)定的高壓、高應(yīng)變和高變形率狀態(tài)，使靶板周圍快速形成塑性變形區(qū)和高溫區(qū)，提供了有利于侵徹正常進(jìn)行的條件，并自碰撞點(diǎn)向靶板中傳入較強(qiáng)的球形沖擊波，同時(shí)在靶板表面產(chǎn)生反射波和拉伸波。由于頭部相對(duì)較尖，接觸界面形成小面積凹陷，該階段持續(xù)時(shí)間非常短暫，約4μs左右。圖5中69~72μs時(shí)均屬于開坑階段。

（2）穩(wěn)定侵徹階段（鼓包形成）：在開坑階段完成之后，彈靶材料繼續(xù)不斷地破壞、飛濺，EFP繼續(xù)沖擊靶板，沖擊速度明顯低于EFP速度，在此過程中EFP不斷縮短，質(zhì)點(diǎn)向側(cè)向流動(dòng)擴(kuò)孔，使侵徹孔徑明顯大于EFP杵體直徑。同時(shí)，杵體由于受到高溫和磨蝕的作用，發(fā)生塑性變形，使其質(zhì)量不斷減少，動(dòng)能降低，截面積增大，靶板內(nèi)部孔徑也隨之增大，出現(xiàn)新表面。破碎的侵徹體碎片又不斷地撞擊彈坑的新表面，如此重復(fù)使彈坑不斷加深。當(dāng)EFP侵徹到靶板一定深度時(shí)，由于穩(wěn)定侵徹及入射波和反射波的共同作用，在靶板的背面產(chǎn)生金屬材料的塑性流動(dòng)，進(jìn)而形成鼓包區(qū)，鼓包近似為球缺形。如圖5中72~84μs為穩(wěn)定侵徹階段。

（3）尾翼侵徹階段：在84μs時(shí)尾翼開始撞擊靶板，此時(shí)，EFP加速度增至最大，速度衰減最快，靶板孔徑也增至最大，由于尾翼對(duì)入口的撞擊作用，造成入口直徑的進(jìn)一步增大。約96μs時(shí)，靶板背面出現(xiàn)明顯隆起，該現(xiàn)象對(duì)沖擊過程將產(chǎn)生一定的影響，但由于材料沒有出現(xiàn)破壞，又在一定程度上限制了該影響。因此，整個(gè)沖擊過程中，靶板背面對(duì)侵徹的影響都很小，近似分析中可以不予考慮。

（4）沖塞貫穿階段（鼓包的破裂直至靶后破片的形成）：隨著侵徹的繼續(xù)進(jìn)行，EFP的磨蝕量增大，其截面積不斷減小，孔徑也隨之減小，在EFP的侵徹和破片反擠的不斷作用下，鼓包的高度不斷增加，并因拉應(yīng)力的持續(xù)增大而開始自外表面破裂。隨后，隨著靶板的抗剪切性能的下降，鼓包的高度繼續(xù)增加，最終彈坑周圍的鼓包區(qū)沿鼓包周邊的應(yīng)力集中及初始裂紋區(qū)域產(chǎn)生拉伸斷裂，鼓包完全破裂。這時(shí)剩余EFP從鼓包中沖出，其后面跟隨著侵徹體碎片和靶板材料的崩落碎片，形成具有殺傷力的靶后破片。如圖5中114~168μs為沖塞貫穿階段。

侵徹結(jié)束時(shí)裝甲靶板的破壞特征主要有以下3種：一是入口呈現(xiàn)卷邊花瓣?duì)钇茐?，造成該現(xiàn)象主要是由于侵徹過程中靶板和EFP均呈明顯的流體動(dòng)力特性，在高速侵徹時(shí)質(zhì)點(diǎn)流動(dòng)產(chǎn)生的離心力和應(yīng)力波反射產(chǎn)生的強(qiáng)拉伸應(yīng)力共同作用，推動(dòng)入口邊緣外翻，使入口呈現(xiàn)花瓣?duì)罹磉吰茐?，這與試驗(yàn)現(xiàn)象基本一致；二是出口呈現(xiàn)外翻花瓣形穿孔且具有明顯的拉伸斷裂特征，主要是因?yàn)镋FP在靶板內(nèi)運(yùn)動(dòng)推動(dòng)靶板背面出現(xiàn)隆起，隆起部分在侵徹過程后期受到EFP的推動(dòng)進(jìn)一步變形，最后拉伸應(yīng)力超過材料的拉伸強(qiáng)度，在EFP四周產(chǎn)生星形裂縫，EFP貫穿靶板后，靶板的拉伸應(yīng)力會(huì)把已產(chǎn)生裂縫的邊緣拉住，并在強(qiáng)度較弱的地方斷裂，形成出口處的外翻花瓣形破壞；三是靶板入口直徑明顯大于出口直徑，且出現(xiàn)兩頭直徑大中間孔徑小的現(xiàn)象。分析侵徹過程可知，EFP尾翼對(duì)入口的擴(kuò)孔作用是造成入口直徑較大的主要原因，而侵徹過程中因磨蝕使得EFP直徑和質(zhì)量減少造成了靶板孔徑縮小，隨著EFP速度的降低，材料的主要變形破壞機(jī)制雖然仍是流體動(dòng)力特征，但其塑性逐漸呈現(xiàn)，EFP出現(xiàn)一定的蘑菇頭形狀，造成侵徹孔徑又有一定的擴(kuò)大。

3 結(jié)論

本文通過數(shù)值模擬、試驗(yàn)研究和理論分析等方法得出以下結(jié)論：（1）從原理上給出了靶板出入口呈現(xiàn)不同花瓣?duì)钇茐?、沖塞和入口直徑明顯大于出口直徑等現(xiàn)象的物理原因。（2）從垂直侵徹到傾斜侵徹，穿孔的形狀由圓形變成橢圓形，并且傾斜角度越大橢圓長(zhǎng)短半軸比越大。（3）分析整個(gè)侵徹過程，可以得到一些提高裝甲防護(hù)能力的措施，即采用爆炸反應(yīng)裝甲(ERA)來破壞或防止EFP開坑等。同時(shí)也可以發(fā)現(xiàn)EFP侵徹靶板除了需要較高的速度外，還需有足夠的長(zhǎng)度消耗在侵徹過程中，為此可通過適當(dāng)方法(如改變藥型罩結(jié)構(gòu))來增加EFP的長(zhǎng)徑比，從而提高其侵徹能力。

[1] 李向東,錢建平,曹兵等.彈藥概論[M].北京:國(guó)防工業(yè)出版社,2004.

[2] 周翔,龍?jiān)?岳小兵.76mm口徑EFP成形過程數(shù)值模擬及影響因素研究[J]. 彈道學(xué)報(bào), 2003,15( 2): 59-63.

[3] 魏濤,黃德武,呂國(guó)皓,等.藥型罩厚度和炸藥性能參數(shù)對(duì)爆炸成型彈丸的影響[J]. 沈陽理工大學(xué)學(xué)報(bào), 2007, 26(6): 54-57.

[4] 高永宏,劉天生,王旭.小口徑雙金屬藥型罩EFP的侵徹效應(yīng)研究[J].彈箭與制導(dǎo)學(xué)報(bào), 2013 (04):83-86.

[5] 李傳增,王樹山,榮竹.爆炸成型彈丸對(duì)裝甲靶板的高速?zèng)_擊效應(yīng)研究[J].振動(dòng)與沖擊,2011(04):91-94.

[6] 王儒策.彈藥工程[M].北京:北京理上大學(xué)出版社,2002.

[7] Century Dynamics Inc. Autodyn user manual[M].Century Dynamics Inc, San Ramon, 2005

Numerical Simulation and Experimental Research on the Process of Large Caliber EFP Penetrating the Typical Armor Plate

YE Yan1, YAO Zhi-min1, LI Jin-ming2, LIU Bo1

（1. Missile Engineering Department, Ordnance Engineering College, Shijiazhuang, 050003；2.Ammunition Engineering Department, Ordnance Engineering College, Shijiazhuang, 050003）

Through the test of large caliber EFP penetrating into the armor target plate speedy, as well as the numerical simulation of the whole penetration process using the finite element simulation software AUTODYN 3D, the EFP open pit, steady penetration, tail penetration and plugging formed through a secondary fragment of the physical process, were studied. Simulation results with the experimental phenomenon and theoretical analysis are in good agreement, and the reasons of the macroscopic phenomena generation were analyzed in principle. The research results not only show the mechanism of EFP penetrating armour target board, but also provide a reference for shaped charge of typical armored damage assessment.

EFP；Penetration mechanism；Armor target；Numerical simulation

1003-1480（2016）02-0017-04

TJ410.3

A

2015-12-22

葉嚴(yán)（1991 -），男，在讀碩士研究生，主要從事導(dǎo)彈裝備仿真與控制研究。