石軍磊，劉迎彬，胡曉艷，張旭光

(中北大學化工與環(huán)境學院，山西 太原 030051)

輔助藥型罩材料對超聚能射流性能影響的數(shù)值模擬

石軍磊，劉迎彬，胡曉艷，張旭光

(中北大學化工與環(huán)境學院，山西 太原 030051)

為了研究超聚能射流形成過程中輔助藥型罩材料對射流性能的影響，利用AUTODYN-2D軟件，采用高精度多物質求解器Euler-2D Multi-material對輔助藥型罩材料分別為Wu、Cu、Fe，錐形藥型罩材料為Al的截頂輔助藥型罩超聚能裝藥進行模擬。結果表明，與傳統(tǒng)聚能裝藥相比，超聚能裝藥結構形成的射流性能更優(yōu)。輔助藥型罩材料密度越大，超聚能射流形態(tài)和連續(xù)性越好；超聚能射流具有更高的頭部速度和能量。用Wu、Cu、Fe作輔助藥型罩材料時，超聚能射流的最大速度分別為14400、13300和13100m/s，最大能量分別為10.2×107、8.5×107、7.5×107J/kg，并且材料密度越大，起爆后相同時刻的射流形態(tài)越細越長。

AUTODYN；超聚能射流；輔助藥型罩；數(shù)值模擬

引 言

在反裝甲聚能裝藥戰(zhàn)斗部中，聚能金屬射流得到了廣泛應用。隨著裝甲防護系統(tǒng)的不斷改進，裝甲防護能力逐漸提升，聚能射流的侵徹能力也需要進一步提高。近年來，研究人員分別從藥型罩材料成分配比、含能材料藥型罩和裝藥結構等方面進行了研究，研制出更具侵徹能力的聚能裝藥戰(zhàn)斗部[1-2]。

俄羅斯V.F.Minin等[3]對截頂加輔助藥型罩裝藥結構進行了研究，結果表明，超聚能裝藥形成的射流與傳統(tǒng)射流相比具有更高的速度和穿深能力。王淦龍[4]以截頂輔助藥型罩為研究對象，結合數(shù)值模擬方法對裝藥結構在不同裝藥長徑比、藥型罩錐角和不同輔助藥型罩厚度下形成的超聚能射流的速度和能量進行了分析。王成等[5]對截頂加錐形鈦合金輔助藥型罩超聚能裝藥進行了研究，得出該裝藥結構可以明顯提高聚能射流的整體速度和有效質量分數(shù)的結論。李慶鑫等[6]對“蘑菇形”超聚能裝藥結構進行了研究，結果表明選用銅、鉭、鎢作藥型罩材料可以獲得良好成型效果的高速射流。

綜上所述，國內學者在超聚能裝藥相關領域進行了一定的研究，但對于輔助藥型罩材料對超聚能射流性能影響的報道較少。本研究以截頂輔助藥型罩為對象，通過改變輔助藥型罩的材料，對不同情況下超聚能射流的形成情況進行了數(shù)值模擬，探討了輔助藥型罩材料對超聚能射流性能的影響規(guī)律。

1 模擬計算

1.1 計算模型

在超聚能裝藥結構中，輔助藥型罩材料分別選用Wu、Cu、Fe，錐形罩材料為Al。本研究所涉及的模型具有軸對稱性質，采用二維建?？梢詽M足要求，因此選用AUTODYN-2D模擬軟件進行模擬。模型涉及多物質大變形運動，為提高計算精度，采用AUTODYN高精度多物質求解器Euler-2D Multi-material求解器進行數(shù)值模擬。

超聚能裝藥結構示意圖如圖1所示，錐形藥型罩為等壁厚，藥型罩的錐角為60°，厚度為1.73mm，材料為Al；輔助藥型罩厚度為4mm，直徑為20mm，材料分別為Wu、Cu、Fe，截頂處直徑為14mm，不同材料裝藥結構下炸藥均為HMX，殼體材料為鋼，厚度為2mm。模擬軟件建立模型如圖2所示。

圖1 超聚能裝藥結構示意圖Fig.1 Schematic diagram of hyper shaped charge structure

圖2 超聚能裝藥模型Fig.2 Model of hyper shaped charge

1.2 計算方法

炸藥爆炸后，在爆轟波的作用下，藥型罩被壓垮，形成高溫、高速、高能量密度的聚能射流，其過程十分復雜，為簡化理論模型，作出如下假設[7]：

(1)炸藥、藥型罩、殼體和空氣為連續(xù)均勻介質；(2)超聚能射流速度隨射流長度呈線性分布；(3)忽略超聚能射流的重力作用；(4)整個射流的形成過程為絕熱過程。

采用AUTDYN-2D Euler算法進行數(shù)值模擬，建立1/2模型。初始條件設置為Still Air，邊界條件為Flow-out。

炸藥采用奧克托金(HMX)，狀態(tài)方程為JWL，該方程能精確地描述爆轟氣體產(chǎn)物的體積、壓力、能量等特性。輔助藥型罩材料均從Autodyn-2D標準材料庫中選取，鎢選用材料TUNGSTEN，采用Shock狀態(tài)方程和Steinberg強度模型；銅選用COPPER，采用Shock狀態(tài)方程和PiecewiseJC強度模型；鐵選用IRON，采用Shock狀態(tài)方程和PiecewiseJC強度模型，鋁選用Al2024T351，采用Shock狀態(tài)方程和Johnso-Cook強度模型。殼體材料選用Steel 1006，采用Shock狀態(tài)方程和Johnso-Cook強度模型。各材料具體參數(shù)如表1和表2所示。

表1 HMX材料模型及其狀態(tài)方程參數(shù)

表2 不同材料的材料模型及其狀態(tài)方程參數(shù)

2 結果與討論

2.1 輔助藥型罩材料對超聚能裝藥射流形成過程的影響

采用控制變量法進行對比分析，裝藥高度為40mm，裝藥直徑為44mm，錐形藥型罩為等壁厚，材料為Al，錐角2α為60°，通過改變輔助藥型罩的材料，在不同結構下觀察不同輔助藥型罩材料時超聚能裝藥結構形成射流的頭部最大速度，以及起爆后10、30、50μs形成的射流速度及其形態(tài)，并對射流長度、連續(xù)性、能量進行分析，其中，輔助藥型罩材料為Wu、Cu、Fe。將形成的超聚能射流與不添加輔助藥型罩情況下的傳統(tǒng)射流進行對比分析，研究不同輔助藥型罩材料對聚能射流的作用。數(shù)值模擬過程中，以罩口軸線為坐標原點，起爆方式采用平面起爆，起爆面模型直徑位置坐標為(-66，10)和(-66，-10)，起爆位置示意圖如圖3所示。

圖3 線起爆位置示意圖Fig.3 Schematic diagram of line detonation

圖4～圖6分別為無輔助藥型罩的傳統(tǒng)聚能裝藥(以下簡稱無)和添加不同材料輔助藥型罩的超聚能裝藥起爆10、30和50μs后的射流形態(tài)。

圖4 10μs時的射流形態(tài)Fig.4 Jet shape corresponding in 10μs

超聚能裝藥起爆后，起爆處的爆轟波以極快的速度在炸藥內部傳播，當爆轟波傳到藥型罩后，藥型罩將承受極大的爆轟壓力(最大能達到200GPa，平均壓力大于20GPa)，在巨大的壓力作用下，藥型罩被壓垮，以2000～3000m/s的速度向藥型罩對稱軸閉合。被壓垮的藥型罩在極短的時間間隔內變形非常大，應變率達到104～107s-1，最大應變超過10。

從圖4可以看出，起爆10μs后，藥型罩還未完全壓垮，普通聚能裝藥射流形成過程較慢，從形態(tài)圖中可以看出，未加輔助藥型罩時，藥型罩壓垮形成的射流較粗，射流頭部呈圓形狀，并且杵體部分較大；當加入輔助藥型罩后，形成的射流頭部成細錐形，杵體也相對越少。隨著起爆時間的延續(xù)，射流被拉長，可以看出，輔助藥型罩能夠使形成的射流更加細長，杵體部分減少，射流有效含量較大，并且輔助藥型罩密度越大，射流越細越長。

圖5 30μs時射流形態(tài)Fig.5 Jet shape corresponding in 30μs

圖6 50μs時射流形態(tài)Fig.6 Jet shape corresponding in 50μs

由圖6可知，聚能裝藥起爆50μs后，無輔助藥型罩時，射流在前期的連續(xù)性較好，隨著時間的增加，斷裂開始明顯，并且斷裂間隙較大。Wu作輔助藥型罩材料時，形成的射流有5處斷裂間隙，由于Wu的密度和硬度較大，并沒有形成射流。Cu作輔助藥型罩材料時，在相同時刻形成的射流長度比Wu的短，并且有6處斷裂間隙。Fe作輔助藥型罩材料時，射流斷裂間隙較大，射流連續(xù)性較差。表4給出30μs和50μs時射流的長度。

表3 30 μs和50μs時超聚能射流長度

由以上結果可以得出，輔助藥型罩沒有充分形成射流，而是與射流尾部混合在一起。起爆30μs后，輔助藥型罩與錐形罩出現(xiàn)分離的趨勢，50μs時，兩者已完全分離。加上輔助藥型罩后，射流斷裂間隙減少，材料為Wu時的超聚能射流主體部分連續(xù)性最好。

與普通聚能裝藥相比，超聚能裝藥形成的超聚能射流連續(xù)性相對較好，增加了射流的有效長度，射流比較細長，杵體較少，輔助藥型罩密度越大，此特征越明顯。

2.2 輔助藥型罩材料對超聚能射流速度影響

在模型中心軸線上設置固定觀測點，以裝藥口徑的倍數(shù)設置高斯點11個，以觀測射流在軸線上的速度分布和能量變化。觀測點1的坐標為(0，0)，觀測點2的坐標為(20，0)，依次排列，如圖7所示。不同裝藥結構形成的射流速度如圖8所示。

圖7 高斯觀測點分布示意圖Fig.7 Schematic diagram of gauge observation points distribution

圖8 不同裝藥結構射流速度分布Fig.8 Velocity distribution of jet corresponding to different shaped charge

從圖8可以看出，普通聚能裝藥結構形成的聚能射流頭部最大速度只有 9200m/s，速度下降梯度較小。使用Wu作輔助藥型罩材料，其超聚能裝藥結構射流的頭部最大速度為14400m/s，隨后降到13000m/s，然后以較小的梯度逐漸遞減。表5給出了不同輔助藥型罩材料對超聚能裝藥結構射流的頭部速度的影響。

表5 在不同時間所測不同裝藥結構射流頭部速度

表5表明，用Wu作輔助藥型罩材料時，射流頭部最大速度較高，是普通聚能射流的1.54倍；使用Cu材料時，射流頭部最大速度是普通聚能射流的1.43倍；使用Fe材料時，射流頭部最大速度與Cu相差不大，為1.40倍。

使用Wu、Cu、Fe作輔助藥型罩材料時，10μs開始聚能射流頭部速度明顯高于普通射流頭部速度，并且材料密度越大，射流頭部速度也越大。同時可以看出，隨著時間的增加，速度梯度都是從大到小逐漸變化。

對于普通聚能裝藥結構，藥型罩兩側受到爆轟波的作用，從而沿軸線方向被壓垮。圖9為普通聚能裝藥和超聚能裝藥在起爆6μs后壓垮初始受爆轟波作用壓力云圖。由圖9可知，加上輔助藥型罩形成超聚能裝藥結構后，爆轟波對藥型罩的作用被改變，錐形罩兩側壓力明顯增大，同時，爆轟壓力通過輔助藥型罩作用于錐形罩截頂處，錐形罩壓垮情況也隨之改變。此時，錐形罩受到更大爆轟壓力后，產(chǎn)生更高的壓垮速度。因此，超聚能裝藥能夠產(chǎn)生速度更高的射流。

圖9 6μs時藥型罩受爆轟波壓力云圖 Fig.9 The pressure nephogram subjected to detonation wave for liner in 6μs

2.3 超聚能射流能量分析

在模型軸線上設置11個高斯點，如圖7所示，每個高斯點可以觀測出通過射流的頭尾速度和能量。藥型罩在爆轟波的作用下，瞬間被壓垮而形成射流，此過程會釋放很大的能量，聚能射流正是以高速、高能的特征對目標形成毀傷。圖10為不同時刻在高斯觀測點處聚能射流的能量分布。

圖10 不同時刻聚能射流能量分布Fig.10 Energy distribution of shaped charge jet at different times

由圖10可知，與普通聚能裝藥相比，超聚能裝藥在起爆后10μs前均有很高的能量，但之后大幅度降低。用Wu作輔助藥型罩材料時，最大能量超過108J/kg，有3個觀測點處的能量超過8×107J/kg；用Cu作藥型罩材料時，最大能量超過8×107J/kg；用Fe作藥型罩材料時，最大能量為7.5×107J/kg。普通射流最大能量不是在起爆之初生成，而是在50μs左右，約為2.8×107J/kg，比超聚能射流能量小的多。綜上可知，超聚能射流比普通射流具有更高的能量，并且最大能量與輔助藥型罩的材料有關，輔助藥型罩密度越大，射流能量越高。

從模擬結果可以看出，如果選擇合適的輔助藥型罩材料，超聚能射流具有很高的速度和內能；從射流形態(tài)上看，超聚能射流更細更長，連續(xù)性較好，射流的有效質量分數(shù)更大。這對聚能裝藥破甲技術的進一步研究和發(fā)展具有一定的意義。

3 結 論

(1)運用AUTODYN-2D軟件，采用高精度多物質求解器模擬計算了輔助藥型罩材料對超聚能射流形成過程的影響。結果表明，與普通聚能裝藥結構相比，超聚能裝藥結構形成的射流形態(tài)比較細長，并且在相同錐角情況下，超聚能射流形態(tài)與輔助藥型罩材料密度有關。密度越大，射流越細越長，并且射流連續(xù)性好，杵體越少，藥型罩利用率越高。

(2)超聚能裝藥形成的超聚能射流頭部最大速度較普通聚能裝藥大很多，材料的密度越大，射流頭部最大速度也越大。用Wu作輔助藥型罩材料時，超聚能射流頭部最大速度是普通聚能射流的1.54倍；用Cu作材料時，超聚能射流頭部最大速度是普通聚能射流的1.43倍；用Fe作材料時，超聚能射流頭部最大速度是普通射流的1.40倍。

(3)超聚能射流具有更高的內能，并且隨著輔助藥型罩材料密度的增加，射流能量也隨之增大。用Wu、Cu、Fe作輔助藥型罩材料時，對應的聚能射流最大內能分別為10.2×107、8.5×107、7.5×107J/kg，而普通聚能射流最大能量只有2.8×107J/kg。

[1] 高永宏,顧曉輝,王鳳英.粉末藥型罩的材料密度對聚能射流性能的影響[J].高壓物理學報,2013(4):556-560.GAO Yong-hong,GU Xiao-hui,WANG Feng-ying.Powder liner lensity effect on the performationce of shaped charge jet[J].Chinese Journal of High Pressure Physics,2013(4):556-560.

[2] 萬文乾,余道強,彭飛,等.含能材料藥型罩的爆炸成型及毀傷作用[J].爆炸與沖擊,2014(2):235-240.WAN Wen-qian,YU Dao-qiang,PENG Fei,et al.Formation and terminal effect of explosively-formed penetratormade by energetic materials[J].Explosion and Shock Waves,2014(2):235-240.

[3] Minin V F,Minin O V,Minin I V,et al.Physics hypercumulation and comdined shaped charges[C]∥2012 11th International Conference on Actual Problems of Electronics Instrument Engineering.Novosibirsk:[s.n.],2012:95-113.

[4] 王淦龍.超聚能射流數(shù)值模擬的初步研究[D].北京:北京理工大學,2015.WANG Jin-long.Elementary study of numerical simulation on hypercumulation[D].Beijing: Beijing Institute of Technology,2015.

[5] 王成,錢俊松,王萬軍.超聚能射流形成的數(shù)值模擬[C]∥北京力學會第20屆學術年會論文集.北京:北京力學會,2014:2.WANG Cheng,QIAN Jun-song,WANG Wan-jun.The numerical simulation on hypercumulation[C]∥ Proceedings of the 20th Annual Meetings of Insitute of Beijing.Beijing:Beijing Society of Theoretical Applied Mechanics,2014:2.

[6] 李慶鑫,王志軍,陳莉,等.一種超聚能裝藥結構的仿真[J].兵器裝備工程學報,2016(6):35-38.LI Qing-xin,WANG Zhi-jun,CHEN Li,et al.Simulation of research of super shaped charge structure[J].Journal of Ordnance Equipment Engineering,2016(6):35-38.

[7] 王偉力,李永勝,田傳勇.串聯(lián)戰(zhàn)斗部前級環(huán)形切割器的設計與試驗[J].火炸藥學報,2011,34(2):39-43.WANG Wei-li,LI Yong-sheng,TIAN Chuan-yong.Optimazition and test on front annular cutter of tandem warhead[J].Chinese Journal of Explosives & Propellants(Huozhayao Xuebao),2011,34(2):39-43.

Numerical Simulation of Effect of Material of Additional Liner on the Performances of Hypercumulation

SHI Jun-lei,LIU Ying-bin,HU Xiao-yan,ZHANG Xu-guang

(School of Chemical Engineering and Environment,North University of China,Taiyuan 030051,China)

To study the effect of material of additinal liner on the jet performances in the hypercumulation formation process,the numerical simulation of hyper shaped charge of a trunconical additinal liner,with material of additional liner is Wu,Cu and Fe,respectively and material of conical shaped charge liner is Al,was carried out using AUTODYN-2D simulation software and high precision multi material solver Euler-2D Multi-material.The results show that the maximum velocity of hypercumulation of additional liner material as Wu,Cu and Fe are 14400,13300 and 13100m/s,respectively,the maximum energy are 10.2×107,8.5×107and 7.5×107J/kg respectively,and the hypercumulation after detonation becomes slender and longger after the density of the material increases.

AUTODYN; hypercumulation; additional liner;numerical simulation

10.14077/j.issn.1007-7812.2017.01.014

2016-08-29；

2016-11-21

國家自然科學基金(No.11572292)

石軍磊(1991-)，男，碩士研究生，從事聚能裝藥和裝甲防護研究。E-mail:sh_junlei@163.com

劉迎彬(1985-)，男，博士，講師，從事彈藥毀傷防護與安全工程方面研究。E-mail:liuyb85@mail.ustc.edu.cn

TJ55;TJ410.3+33

A

1007-7812(2017)01-0069-06