一種M型藥型罩形成環形射流的數值模擬研究

童宗保，王小春，程　劍，劉　燕，胡善寶，

劉智華，雷勇濤，胡潤娟，馮　梅

(湖南湘潭江南工業集團有限公司，湖南 湘潭　411207)

?

一種M型藥型罩形成環形射流的數值模擬研究

童宗保，王小春，程劍，劉燕，胡善寶，

劉智華，雷勇濤，胡潤娟，馮梅

(湖南湘潭江南工業集團有限公司，湖南 湘潭411207)

摘要：對串聯戰斗部前級聚能裝藥開孔不大的問題開展了研究，提出了一種M型藥型罩，通過分析該藥型罩結構設計的原理，運用ANSYS/LS-DYNA仿真軟件對該藥型罩形成環形聚能射流進行了數值模擬研究。結果表明：M型藥型罩能夠形成直徑較大的環形射流，射流持續穩定；軍用領域，能夠用作串聯戰斗部的前級藥型罩，對靶板切割成一定孔深且孔徑大的環形孔，方便后級隨進戰斗部對靶板形成充塞破壞；民用領域，可以對一定厚度的靶板進行快速打孔。

關鍵詞：M型藥型罩；環形切割；大孔徑;串聯戰斗部

Citation format:TONG Zong-bao, WANG Xiao-chun, CHEN Jian, et al.Study on Numerical Simulation of M Shaped Charged Liner form Annular Jet[J].Journal of Ordnance Equipment Engineering,2016(3):143-145.

聚能裝藥結構主要用來侵徹和破壞某些特殊的目標，如裝甲車、機場跑道、地下指揮所等一些典型的結構，在軍用領域，其多用在兩級串聯戰斗部的設計當中[1]，前級戰斗部首先作用，利用聚能效應對硬目標開孔為后級隨機戰斗部的侵徹開辟通道，問題的關鍵是設計前級聚能裝藥對目標切割出一個直徑足夠大的孔；對民用領域，可以用來石油開采，爆炸切割等。

目前，串聯戰斗部中前級使用的藥型罩大部分為錐形、雙錐形、半球形、喇叭形，弧錐形等[2-6]，這些藥型罩均不能在毀傷的目標上形成直徑足夠大的孔，不利于后級隨進戰斗部的進一步毀傷。盡管已經有學者對環形切割進行了研究，但藥型罩剖面還主要集中在錐形結構上[7]，針對剖面是M型藥型罩結構形成環形聚能射流過程、速度和孔徑大小的研究還不很多。

本研究提出了一種M型藥型罩，能夠在炸藥的爆炸作用下形成直徑足夠大的環形射流，在目標上形成大的開孔且具有一定的孔深，用作串聯戰斗部前級聚能裝藥對目標靶進行環形切割，方便后級隨進戰斗部對目標靶形成充塞破壞，而且還可以用來對一定厚度的鋼板進行快速打孔。本文作者以理論分析為基礎，通過數值模擬的方法對其形成直徑足夠大的射流進行仿真分析。

1理論分析

圖1為M型藥型罩以及其裝藥結構設計幾何模型，藥型罩的壁厚的為δ，兩邊對稱的“倒V”型罩錐角為α，兩對稱軸之間夾角為θ，兩對稱軸的交點O為炸藥的起爆點。

炸藥在O點起爆之后，就會形成球形爆轟波，當球形波到達藥型罩的罩頂A和A′時，就會對藥型罩進行壓垮，罩頂位置形成杵體，由于兩對稱軸在球形爆轟波的徑向，所以形成的聚能射流不會彎曲，而是兩道成直線的聚能射流。

圖1　M型藥型罩及裝藥結構幾何模型

2數值模擬研究

2.1計算模型和算法介紹

采用LS-DYNA仿真分析軟件對M型藥型罩形成環形聚能射流的過程進行數值模擬仿真，計算模型如圖2所示。有限元計算模型包括炸藥和藥型罩，藥型罩壁厚δ為3 mm，兩邊對稱的“倒V”型罩錐角α為61°，兩對稱軸之間夾角θ為16°，計算模型采用solid162二維實體單元進行網格劃分，炸藥和藥型罩之間采用CONTACT_2D_AUTOMATIC_SURFACE_TO_SURFACE接觸算法，由于藥型罩在形成射流的過程中，單元會產生大變形，因此要使用自適應網格技術。當計算時間達到60 μs時，炸藥基本爆轟完畢，對射流的形成影響很小，因此可以刪除炸藥的part以及炸藥與藥型罩之間的接觸，炸藥和藥型罩均采用Lagrange算法，計算單位為cm-g-μs。

圖2　計算模型

2.2材料模型和狀態方程

炸藥選用PBX9010炸藥，采用高能炸藥模型(HIGH_EXPLOSIVE_BURN)和JWL狀態方程進行描述，相關參數[8]見表1所示，JWL狀態方程精確的描述了在爆炸驅動過程中，爆轟氣體產物的壓力、體積、能量的特征，JWL狀態方程的表達式為

(1)

式(1)中：P是等熵壓力；V是爆轟產物的相對體積；A、B、R1、R2、ω為輸入的參數；E為體積內能。

表1　炸藥材料模型和狀態方程參數

藥型罩材料為銅，選用適用于大變形的Johnson-Cook材料本構模型，采用Gruneisen狀態方程來描述。其本構模型為

(2)

式(2)中：σeq是等效應力；A是屈服應力；B是應變硬化系數；N是應變硬化指數；M是溫度相關系數；ε*是無量綱塑性比，ε*=εeq/ε0；εeq為等效塑性應變；ε0為參考塑性應變，一般取ε0為1s-1；T*為無量綱溫度，T*=(T-Tr)/(TM-Tr)；T為當前溫度；Tr為室溫；TM為熔點溫度，當溫度T到達融化溫度TM時，屈服應力為零。

Gruneisen狀態方程的具體表達式在壓縮狀態下，即μ>0：

(3)

在膨脹狀態下，即μ<0:

(4)

表2　藥型罩材料模型和狀態方程參數

3數值模擬結果分析

3.1聚能射流形成分析

通過數值模擬的方法可以看到M型藥型罩形成聚能射流的全過程，一些典型時刻的結果如圖3所示。

圖3(a)為M型藥型罩形成環形聚能射流的初始狀態；圖3(b)為炸藥爆轟46 μs后，爆轟波剛傳播到M型藥型罩的頂端，與藥型罩剛接觸，接下來爆轟波開始壓垮藥型罩；圖3(c)為炸藥爆轟60 μs后，炸藥爆轟基本完成，此時炸藥對射流的影響很小，可以刪除炸藥的part以及炸藥與藥型罩之間的接觸；圖3(d)為M型藥型罩在爆轟波壓垮的作用下，杵體開始形成；圖3(e)為藥型罩閉合為射流的自由運動，由于速度梯度的存在，聚能射流在運動中逐漸拉伸；圖3(f)表示M型藥型罩在爆轟波的驅動下能夠形成細長而且穩定的射流。

圖3　聚能射流形成過程

3.2聚能射流速度分析

在聚能射流頭部選取編號為6581的節點，通過后處理軟件觀看其速度隨時間變化曲線，結果如圖4所示。

圖4　時間-速度曲線

根據圖4分析可知：總體來說，M型藥型罩在炸藥爆轟的驅動作用下形成的射流速度不高，達不到錐形罩的射流速度。約125 μs時刻，射流速度達到最大值，約為3 300 m/s，之后逐漸衰減。射流速度不高的原因可能跟裝藥結構和藥型罩的材料有關，接下來的工作是對裝藥結構和藥型罩材料進行進一步的優化，提高射流的速度。

3.3形成環形射流孔徑分析

圖5為環形聚能射流頭部的局部放大圖，通過后處理軟件可以測量射流頭部兩點之間的距離，即為環形射流的孔徑大小。

圖5　環形射流孔徑

由圖5分析可知，該M型藥型罩形成的環形聚能射流的孔徑為184 mm，孔徑很大，說明該藥型罩可以用作串聯戰斗部的前級藥型罩，碰靶瞬間前級裝藥起爆，爆轟波壓垮M型藥型罩，形成的環形射流可以對目標靶進行環形切割，在目標靶上形成一個具有一定深度且孔徑很大的環形孔，該孔徑大于現在很多反坦克導彈和反艦導彈后級隨進戰斗部的口徑[10]，可以方便后級隨進戰斗部對其形成充塞破壞，進而對內部目標毀傷。

4結論

本研究針對一種新型的M型藥型罩形成聚能射流的過程進行了數值模擬仿真分析，得出了以下結論：

1) 該M型的藥型罩在炸藥爆轟的驅動下被壓垮，能夠形成細長而且穩定的射流。

2) 可能由于裝藥結構和藥型罩材料的影響，導致該M型藥型罩形成的聚能射流頭部速度不高。接下來的工作將對裝藥結構和藥型罩的材料進行進一步的優化，提高聚能射流速度并且開展試驗研究，通過試驗驗證該藥型罩設計的合理性以及侵徹毀傷金屬靶板的威力。

3) 該藥型罩形成的環形射流孔徑為184 mm，大于目前大部分的反坦克導彈和反艦導彈的后級隨進戰斗部的外徑，可以用作串聯戰斗部的前級藥型罩，在目標靶上形成一定的外徑深度且直徑較大環形孔，方便后級隨進戰斗部對目標靶形成充塞破壞。

參考文獻：

[1]常敬臻,盧永剛,孫傳杰.聚能裝藥對磚墻結構靶體開孔效能試驗研究[J].兵器材料科學與工程,2009(6):39-41.

[2]張毅,茍瑞君,陳亞紅,等.變壁厚藥型罩對爆轟波對撞EFP成型影響的數值模擬研究[J].兵器材料科學與過程,2015(2):40-43.

[3]李永勝,王偉力,田傳勇.環形切割器藥型罩方案優化設計[J].彈箭與制導學報,2011(4):118-120.

[4]周方毅,詹發民,吳曉鴻,等.圓錐-球缺藥型罩聚能戰斗部結構優化設計[J].爆破器材,2014(6):43-47.

[5]張利，李剛.基于分層序列法的EFP戰斗部藥型罩結構優化與仿真[J].四川兵工學報，2013(8)：23-25.

[6]陳忠勇，李文彬，高旭東，等.藥型罩結構參數對侵徹體成型的影響[J].四川兵工學報，2011(2)：12-15.

[7]傅 磊，王偉力，黃雪峰，等.環型聚能裝藥侵徹靶板能力影響因素分析[J].海 軍 航 空 工 程 學 院 學 報，2015(2):151-155.

[8]呂愿宏,王續躍,李曉杰,等.粉末藥型罩石油射孔彈侵徹性能的數值模擬[J].工程爆破,2015(2):1-4.

[9]時黨勇,李裕春,張勝明.基于ANSYS/LS-DYNA8.1進行顯示動力分析[M].北京:清華大學出版社,2005.

[10]王東生.國外反坦克導彈的發展分析[J].飛航導彈,2010(5):52-57.

(責任編輯唐定國)

Study on Numerical Simulation of M Shaped Charged Liner form Annular Jet

TONG Zong-bao, WANG Xiao-chun, CHEN Jian, LIU Yan, HU Shan-bao,LIU Zhi-hua, LEI Yong-tao, HU Run-juan, FENG Mei

(Hunan Xiangtan Jiangnan Industry Group Co., Ltd., Xiangtan 411207, China)

Abstract:The problem that the precursory shaped charge in the tandem warhead cannot form big hole was studied. We presented a new shaped charge liner of M and the design principle of its structure was analyzed, and numerial simulated the process of this shaped charge liner form annular jet with the ANSYS/LS-DYNA software. The result shows that M shaped charge liner can form annular jet with big hole and the jet is sustained and steady. In military fields, it can be used as a tandem warhead of pre-shaped charge liner and can be cut into the target with a depth and large circular hole which can convenient make the stage warhead to destroy the target with plug. In civil field, it can fastly drill in the target with a certain thickness.

Key words:M shaped charger liner；annular cut；big hole；tandem warhead

文章編號：1006-0707(2016)03-0143-04

中圖分類號：TJ41

文獻標識碼：A

doi:10.11809/scbgxb2016.03.034

作者簡介：童宗保(1988—)，男，助理工程師，主要從事戰斗部設計研究。

收稿日期：2015-08-31；修回日期：2015-09-15

本文引用格式：童宗保，王小春，程劍，等.一種M型藥型罩形成環形射流的數值模擬研究[J].兵器裝備工程學報，2016(3):143-145.

【機械制造與檢測技術】