藥型罩參數(shù)對聚能射流水下運動的影響

歐靖，楊易釗，郁文杰，唐奎

武器裝備

藥型罩參數(shù)對聚能射流水下運動的影響

歐靖a，楊易釗a，郁文杰a，唐奎b

（南京理工大學 a. 錢學森學院 b. 瞬態(tài)物理國家重點實驗室，南京 210094）

研究藥型罩結(jié)構(gòu)參數(shù)對所形成的聚能射流在水中運動的影響，改進水中聚能射流的運動特性。采用多物質(zhì)單元 ALE 法就錐形罩射流對水介質(zhì)的侵徹進行數(shù)值模擬, 分析錐形裝藥結(jié)構(gòu)中藥型罩錐角和厚度對所形成的聚能射流侵徹水時運動參數(shù)的影響。錐形罩錐角大小及藥型罩厚度對聚能射流在水中的形狀、射流速度、加速度等有著明顯的影響。侵徹體進入水中10 cm后，藥型罩的錐角從30°增加到150°的過程中，剩余速度先增大、后減小，在90°時達到最高。藥型罩厚度為1.5～4 mm時，剩余速度變化起伏??；厚度為4～6 mm時，剩余速度開始大幅下降。當錐角為90°時，罩厚為4 mm的藥型罩所形成的射流在水中表現(xiàn)最好，形成的射流侵徹深度最長，侵徹水介質(zhì)10 cm后的剩余速度最大，存速能力最強。

聚能射流；水介質(zhì)；藥型罩；運動特性；ALE法；侵徹體

傳統(tǒng)水中爆破性戰(zhàn)斗部主要通過沖擊波和氣泡脈動對目標進行毀傷，其能量利用率低，毀傷效率低。聚能戰(zhàn)斗部的優(yōu)點是能量比較集中、裝藥利用率較高，故世界上的許多國家相繼開展了在水中運用聚能戰(zhàn)斗部的研究[1-10]。

在實驗和數(shù)值模擬研究方面，國內(nèi)許多學者進行了大量的研究工作，其中大多數(shù)的研究集中在爆炸成型彈丸[11-15]、普通射流[16-20]和桿式射流[21-25]等方面。曹兵[11]通過侵徹實驗研究了EFP戰(zhàn)斗部在空中和水下對靶板的破壞效果，發(fā)現(xiàn)水下戰(zhàn)斗部在水中爆炸形成的EFP對水下目標靶的破壞威力更大。王長利等[12]發(fā)現(xiàn)，聚能裝藥在水下爆炸形成的聚能侵徹體在與水的作用過程中，會形成彈前激波，并且該激波會追趕爆炸沖擊波，之后兩者疊加，使得沖擊波強度增強，從而解釋了曹兵的研究。張向榮等[19]采用Autodyn軟件對材料為鎢銅合金的射流在空氣和水中對鋼板的侵徹深度進行了數(shù)值模擬，結(jié)果表明，無論是在空氣中還是在水中侵徹，都存在一個最合適炸高。當炸高小于４倍的裝藥直徑時，在水中的侵徹深度會大于在空氣中的情形?？梢钥闯觯橘|(zhì)在模型中起著增大炸高的作用，并且水介質(zhì)也會在一定程度上侵蝕射流頭部。張會鎖等[20]考慮在相同裝藥條件下，采用ANSYS中集成的LS-DYNA模塊，進行有限元數(shù)值模擬，在深度考慮并模擬了各種藥型罩壁厚的工況下3種127型石油射孔彈的金屬射流，通過軟件進行數(shù)值仿真來模擬實際的形成過程，同時對模擬結(jié)果進行比較分析。結(jié)果表明，1.5 mm壁厚的藥型罩所產(chǎn)生的射流頭部速度高，動能較大。李明星等[21]采用Autodyn軟件研究了魚雷在水下發(fā)生接觸與非接觸爆炸時對靶板穿孔孔徑的影響，以及不同的聚能裝藥在考慮水下壓力時的爆炸成形運動特性。結(jié)果表明，爆炸成形彈丸頭部速度衰減最慢，桿式射流稍快，普通射流的速度衰減最快。王海福等[22]采用Autodyn軟件通過數(shù)值模擬研究了聚能裝藥在水下的作用行為，初步明確了藥型罩的半徑曲率、中心、形狀等平面數(shù)值對水下發(fā)射的聚能裝藥具有較大的影響作用。

針對藥型罩參數(shù)對射流水中運動特性的影響，國內(nèi)學者也進行了一些研究。周方毅等[26]針對圓錐–球缺組合藥型罩的不同參數(shù)，采用正交試驗設(shè)計的方法設(shè)計了16種戰(zhàn)斗部結(jié)構(gòu)在水中侵徹的模型。數(shù)值計算結(jié)果表明，圓錐–球缺組合藥型罩的結(jié)構(gòu)優(yōu)化方向為，圓錐罩角度可選60°，圓錐罩壁厚在0.2～0.5 cm選取，球缺罩圓心角應(yīng)大于150°，球缺罩的最佳壁厚為0.8 cm。廖莎莎等[27]運用實驗方法對比研究了2種不同金屬材料藥型罩形成的射流在水中的侵徹特性，結(jié)果表明，由鎢銅合金藥型罩形成的射流可以在水中保持較高的剩余速度。王海福等[22]通過研究藥型罩形狀對聚能裝藥水下作用行為，發(fā)現(xiàn)錐形裝藥結(jié)構(gòu)所形成聚能射流在水中運動速度的衰減最快，偏心亞半球形裝藥結(jié)構(gòu)所形成聚能射流在水中運動速度的衰減最慢。同時，他們針對偏心亞半球形銅罩的不同罩厚進行了數(shù)值模擬和實驗驗證，確定偏心亞半球形銅罩壁厚為1.8 mm時，侵徹體水下穿甲能力最優(yōu)。

本文針對錐形裝藥結(jié)構(gòu)所形成聚能射流在水中運動速度衰減過快的問題，采用多物質(zhì)單元 ALE 法進行數(shù)值模擬，研究了錐形罩射流對水介質(zhì)的侵徹，重點探討錐形裝藥結(jié)構(gòu)中藥型罩錐角和厚度對所形成的聚能射流侵徹水時運動參數(shù)的影響。

1 數(shù)值模擬

1.1 模型設(shè)計

圖1 錐形聚能裝藥戰(zhàn)斗部

根據(jù)實際戰(zhàn)斗部結(jié)構(gòu)的對稱性，并結(jié)合實際，考慮數(shù)值模擬計算所需的時間成本與長周期，按照要求建立了1/2戰(zhàn)斗部有限元模型，由炸藥、藥型罩、空氣、水、隔板5部分組成，如圖2所示。將ALE單元定義成流體，它們之間的連接為流固耦合。計算中所有單元類型都是選用8節(jié)點實體單元Solid164，如圖3所示。將聚能裝藥爆炸后模擬計算得到的射流用耦合的方式將壓力和能量沖擊到金屬隔板上，完成了對隔板的侵徹。擊穿隔板后射流進入水中，開始對水的侵徹。采用Euler網(wǎng)格對炸藥、水、空氣和藥型罩進行建模，采用Lagrange網(wǎng)格對隔板進行建模，而單元使用多物質(zhì)ALE算法。計算中，模型的單位制為g-cm-μs，網(wǎng)格邊長為1 mm，并采用0.1 mm厚度實體網(wǎng)格。這種簡化既可以充分利用多物質(zhì)ALE算法，又可以將模型尺寸大大減小。

圖2 數(shù)值模擬模型

圖3 1/2有限元模型網(wǎng)格劃分

1.2 基于ALE網(wǎng)格的流固耦合算法

根據(jù)網(wǎng)格劃分的不同，仿真計算方程的計算方法大體上可以分為歐拉法和拉格朗日法。歐拉法的計算原理是，在網(wǎng)格劃分上使介質(zhì)空間始終保持固定不動，物體運動，在歐拉網(wǎng)格空間中填滿介質(zhì)，其被歐拉網(wǎng)格劃分為諸多離散單元，這種方法適合用來處理大畸變的問題，常用于解決超高速撞擊的問題。但是物質(zhì)交界面必須要作專門的處理，當出現(xiàn)較復(fù)雜的問題時，計算網(wǎng)格的數(shù)目會大幅度提升，花費較長的計算時間。拉格朗日算法的計算原理是，在網(wǎng)格劃分上使物體始終保持固定不動，網(wǎng)格固定在物體內(nèi)部，隨物體運動而運動，且隨物體的變形而變形。在一定程度上，拉格朗日算法能更精確地追蹤到研究物體的邊界，所需的計算網(wǎng)格數(shù)也較少。但是它存在一定的缺陷，就是無法用于計算大畸變和大變形問題，因為當網(wǎng)格發(fā)生大變形時，計算不穩(wěn)定，可能會產(chǎn)生負體積、負質(zhì)量或者其他病態(tài)現(xiàn)象，使計算無法繼續(xù)進行。另一方面，由于計算網(wǎng)格需固定在材料內(nèi)，計算網(wǎng)格隨材料不同而畸變，一旦畸變嚴重，就必須重新對網(wǎng)格進行劃分，否則計算無法繼續(xù)下去[1]。

文中使用ANSYS中集成的LS-DYNA模塊中的ALE（Arbitrary Lagrangian-Eulerian）算法，ALE法是同時包括拉格朗日法與歐拉法的一種算法。ALE算法可以分為3步[13]：第一步是顯示Lagrange計算，即只將壓力梯度分布對速率和能量的影響作為參考指標，取動量方程中前一時刻的壓力數(shù)值，此時式子是顯式格式；第二步是用隱式格式的方法來解開動量方程，與此同時，將第一步所求的速度分量作為接下來一個迭代算法的初始值；第三步，在前面方程的基礎(chǔ)上重新劃分網(wǎng)格，并計算網(wǎng)格之間的運輸量?；贏LE算法的網(wǎng)格結(jié)構(gòu)，在LS-DYNA模塊中可以輕易地運行歐拉–拉格朗日耦合算法，從而實現(xiàn)流固耦合計算。ALE算法能夠很好地解決涉及網(wǎng)格大變形、材料流動的問題。

1.3 材料模型和狀態(tài)方程

1）模型中藥型罩材料選擇紫銅，材料模型選擇Steinberg模型進行描述，其運動的狀態(tài)使用Gruneisen狀態(tài)方程（1）進行近似模擬，狀態(tài)方程（2）為藥型罩膨脹狀態(tài)方程[28]。材料參數(shù)見表1。

式中：為聲速；1、2和3為-曲線斜率的系數(shù)；0為Gruneisen系數(shù)；是對0的一階修正[28]。

表1 紫銅材料參數(shù)

Tab.1 Material parameters of red copper

2）炸藥采用的是TNT，材料模型為High_ Explosive_Burn，狀態(tài)方程為JWl狀態(tài)方程[28]，見式（3）。材料以及式（3）中的輸入?yún)?shù)見表2。

3）隔板材料選擇45鋼，采用隨動硬化模型進行描述，材料參數(shù)見表3[28]。

式中，RI—沉積物中重金屬的潛在生態(tài)風險指數(shù)；重金屬的潛在生態(tài)風險系數(shù)；沉積物中重金屬i相對參比值的污染系數(shù)；Ci—重金屬i的實測含量；重金屬i的評價參比值，本研究以洞庭湖元素的背景值[4,19]作為參比值，Ba、V、Cr、Mn、Co、Ni、Cu、Zn、Pb的評價參比值依次為污染物i的毒性響應(yīng)系數(shù)，反映重金屬的毒性、污染水平和環(huán)境對重金屬污染的敏感程度[22-23]．參考林麗欽[22]和徐爭啟等[23]的研究結(jié)果，9種重金屬的毒性響應(yīng)系數(shù)取值分別為Co = Ni= Cu =Pb = 5 > Ba = V = Cr = 2> Mn = Zn = 1．沉積物中重金屬的潛在生態(tài)風險等級劃分如表1．

4）靶板材料為水，模型中水和空氣均采用MAT_NULL模型以及*EOS_Gruneisen狀態(tài)方程，材料參數(shù)見表4[28]。

表2 TNT材料參數(shù)

Tab.2 Material parameters of TNT

表3 45鋼材料參數(shù)

Tab.3 Material parameters of 45 # steel

2 不同錐角射流對水介質(zhì)侵徹的仿真結(jié)果

2.1 仿真工況

隔板厚度為0.5 mm，藥型罩厚度為2 mm，聚能戰(zhàn)斗部的侵徹炸高保持在50 mm。在確保整體裝藥直徑30 mm和頂部裝藥高度25 mm不變的前提下，通過改變模型的藥型罩錐角來建立不同的有限元模型。起爆點位于主裝藥頂部前后邊界中心位置處，取錐角為30°～150°，每隔30°確定不同仿真研究方案，對不同錐角方案進行數(shù)值模擬計算。

考慮到大型水面艦艇大部分采用多層的板架舷側(cè)防護結(jié)構(gòu)[29]，防護結(jié)構(gòu)一般有多道防雷隔艙，隔艙可以是空艙或者液艙，根據(jù)參考文獻[30]，水艙的寬度一般為10 cm，因此本文針對聚能射流侵入水深10 cm的情況進行討論。

2.2 仿真結(jié)果

不同錐角的藥型罩形成聚能射流擊穿隔板并侵入水深10 cm時的形態(tài)如圖4所示。

表4 空氣和水相關(guān)參數(shù)[23]

Tab.4 Relevant parameters of air and water[23]

圖4 不同錐角形成聚能射流侵入水深10 cm的形態(tài)

將仿真數(shù)據(jù)進行曲線擬合，可以得出聚能侵徹體侵入水深10 cm所需的時間，以及此時的剩余速度隨錐角變化曲線，如圖5所示。

2.3 分析與討論

從圖4和圖5可知，由于小錐角的藥型罩形成射流的杵體質(zhì)量大，射流質(zhì)量小，而射流頭部速度高，射流在擊穿隔板并入水后，頭部與杵體分離，頭部質(zhì)量較小，射流形狀分散，受到水的阻力的影響較為明顯，入水后速度衰減很快，存速能力較差。大錐角的藥型罩在TNT炸藥所產(chǎn)生的爆轟下也會產(chǎn)生射流和對應(yīng)的杵體，但由于杵體的質(zhì)量小，壓合角比較大，藥型罩大部分會轉(zhuǎn)換為頭部，形成直徑較大、速度相對較低的EFP，在擊穿隔板并入水后，整體質(zhì)量集中，因此其入水后速度衰減較慢，受到水的阻力的影響較小。

從圖4可以看出，藥型罩錐角過大時，由于形成的EFP直徑過大，且速度低，在入水前擊穿隔板時，能量消耗嚴重，入水后侵徹體質(zhì)量分散，實際侵徹部分質(zhì)量較小，不利于后續(xù)水中侵徹運動。由圖5可知，隨著錐角的增大，聚能侵徹體侵入水深10 cm的剩余速度先增大、后減小。當錐角從30°增加到90°時，剩余速度增加；當錐角從90°增加到150°時，剩余速度減小。

厚度為2 mm、錐角為90°的藥型罩形成射流侵徹過程如圖6所示。

圖5 射流水中侵徹時間和侵徹速度隨侵水深度的變化

圖6 厚度為2 mm、錐角為90°的藥型罩形成射流的侵徹過程

由圖6中射流的侵徹過程可以看出，圖5a中拐點的形成是由于射流進入水中，頭部和杵體分離，頭部分割水域形成了一個空腔，杵體在空腔中運動阻力較小，隨后頭部速度減弱，杵體在后部推動射流前進，使得頭部速度加快，產(chǎn)生了拐點。綜合數(shù)值模擬的結(jié)果，由圖4和圖5可以看出，錐角為90°的聚能裝藥所形成的聚能射流在水中表現(xiàn)較好，到達水深10 cm所需的時間較短，且射流的剩余速度最大。

3 不同罩厚射流對水介質(zhì)侵徹的仿真結(jié)果

3.1 仿真工況

從第2節(jié)的仿真結(jié)果分析可知，錐角為90°的聚能裝藥所形成的射流在水中表現(xiàn)較好。因此，保持藥型罩的錐角為90°，隔板厚度為0.5 mm，頂部裝藥的高度保持25 mm，直徑保持30 mm，起爆點仍位于主裝藥頂部前后邊界中心位置處，炸高保持在50 mm，只改變藥型罩的厚度，建立有限元模型。

3.2 仿真結(jié)果

不同厚度的藥型罩形成聚能射流擊穿隔板并侵入水深10cm時的形態(tài)見圖7。

將仿真數(shù)據(jù)進行曲線擬合，可以得出射流侵入水深10 cm所需的時間，以及此時的剩余速度隨厚度變化曲線，如圖8和圖9所示。

將計算所得加速度數(shù)據(jù)進行曲線擬合，可以得出射流侵入水深10 cm后，加速度隨厚度的變化曲線，如圖10所示。

圖7 不同厚度形成聚能射流侵入水深10 cm的形態(tài)

圖8 射流侵入水深10 cm所需的時間隨罩厚的變化曲線

圖9 射流侵入水深10 cm后剩余速度隨罩厚的變化曲線

圖10 射流侵入水深10 cm后加速度隨罩厚的變化曲線

3.3 分析與討論

分析圖8—10中的曲線可以得到如下結(jié)論：

1）隨著藥型罩厚度的增大，射流侵入水深10 cm所需的時間是逐漸增大的。當藥型罩厚度為1.5～4 mm時，所需時間平緩增加；當藥型罩厚度為4～7 mm時，所需時間快速增加。

2）隨著藥型罩厚度的增大，射流侵入水深10 cm后，剩余速度具有減小的趨勢。當藥型罩厚度為1.5～4 mm時，射流剩余速度變化起伏小；當藥型罩厚度為4～6 mm時，射流剩余速度開始大幅下降；當藥型罩厚度大于6 mm以后，射流剩余速度開始平緩減小。

3）隨著藥型罩厚度的增大，射流侵入水深10 cm后，速度的衰減具有減小的趨勢。

不同藥型罩厚度條件下，射流在水中的速度衰減程度不同：藥型罩厚度越小，射流初始的動能越大，射流質(zhì)量小，速度衰減越快；藥型罩厚度越大，射流初始的動能越小，射流質(zhì)量較大，速度衰減越緩慢。

這些現(xiàn)象是由于入水初速隨罩厚增加而降低，射流質(zhì)量逐漸增大，水中侵徹過程的速度衰減減緩。這說明罩厚增大是有利于射流水中侵徹能力的提高，但如果罩厚過大，會導(dǎo)致入水速度過低，對后期侵徹效果不利。

分析圖7中的射流形狀可知，厚度較小的藥型罩形成的射流在侵入水深10 cm后形狀分散，實際侵徹體質(zhì)量??；厚度較大的藥型罩形成的射流在侵入水深10 cm后，整體質(zhì)量較為集中，有利于提高后期侵徹效果。因此，在保持裝藥結(jié)構(gòu)不變的前提下，適當增加藥型罩的厚度，可以有效地提高射流在水中的速度抗衰減能力。綜合數(shù)值模擬的結(jié)果，由圖8—10可以看出，當錐角為90°時，罩厚為3～4 mm的藥型罩所形成的射流在水中表現(xiàn)較好，侵入水深10 cm后，剩余速度較大，且速度衰減較慢，具有較好的存速能力。

針對藥型罩厚度4 mm、錐角為90°的較優(yōu)方案，得到射流水中運動的形貌。圖11中展示了0～800 μs的時間段中，每200 μs射流在水中的運動狀態(tài)，其中不同顏色表示不同的流體密度?？梢钥闯觯朔桨杆幮驼之a(chǎn)生的射流呈現(xiàn)出梭子型，水下受阻較小，運動穩(wěn)定且頭部和杵體分離，具有二次毀傷的效果。

圖11 厚度為4 mm、錐角為90°形成射流的侵徹過程

4 結(jié)論

基于 Ls-dyna 有限元軟件，對水下聚能侵徹體的聚能裝藥設(shè)計開展了優(yōu)化，重點針對藥型罩結(jié)構(gòu)進行了仿真研究，分析了藥型罩厚度和錐角對聚能射流侵徹體水下成形和運動的影響，得到了聚能射流侵徹體成形及侵徹水介質(zhì)的全部過程。模擬結(jié)果表明，藥型罩的厚度與錐角會明顯影響其形成的聚能射流在水中的存速能力。

1）增加藥型罩的壁厚會增加射流侵徹體的質(zhì)量，從而提高其存速能力，但是會降低射流的初始動能。適當增加藥型罩的厚度能夠提高射流在水中的速度抗衰減能力，仿真結(jié)果表明，藥型罩厚度為4 mm時，射流的剩余速度最大，侵徹距離最深。

2）增加藥型罩錐角會使得形成的杵體質(zhì)量較小，轉(zhuǎn)換為頭部的藥型罩更多, 形成直徑較大、速度相對較低的EFP，但其受到水中阻力的影響較小。當藥型罩錐角為90°時，形成的射流侵入水深10 cm后，剩余速度最大。

3）根據(jù)數(shù)值模擬的結(jié)果可以看出，當錐角為90°時，罩厚為3～4 mm的藥型罩所形成的射流在水中的運動速度最快，表現(xiàn)較好，侵入水深10 cm后，剩余速度較大，且速度衰減較慢，具有較好的存速能力。

[1] 李裕春, 楊萬江, 沈蔚. 藥型罩曲率半徑對爆炸成型彈丸參數(shù)的影響[J]. 火工品, 2003(1): 45-48.

LI Yu-chun, YANG Wan-jiang, SHEN Wei. Influence of Liner Curvature Radius on Parameters of Explosively Formed Penetrator (EFP)[J]. Initiators & Pyrotechnics, 2003(1): 45-48.

[2] RAWLING W. The Challenge of Modernization: The Royal Canadian Navy and Antisubmarine Weapons, 1944-1945[J]. The Journal of Military History, 1999, 63(2): 355.

[3] HU F, WU H, FANG Q, et al. Impact Performance of Explosively Formed Projectile (EFP) into Concrete Targets[J]. International Journal of Impact Engineering, 2017, 109: 150-166.

[4] LI R, LI W B, WANG X M, et al. Effects of Control Parameters of Three-Point Initiation on the Formation of an Explosively Formed Projectile with Fins[J]. Shock Waves, 2018, 28(2): 191-204.

[5] CHAN S K. An Improvement in the Modified Finite Element Procedure for Underwater Shock Anaslysis[C]//Proceeding of 63nd Shock and Vibration Symposium. Las Cruces: The Shock and Vibration Exchange, 1992.

[6] 張志春, 孫新利, 孟會林, 等. EFP成型影響因素的數(shù)值模擬[J]. 彈箭與制導(dǎo)學報, 2006, 26(S1): 227-228.

ZHANG Zhi-chun, SUN Xin-li, MENG Hui-lin, et al. Numerical Simulation for the Parameter Study of Explosively Formed Projectile[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2006, 26(S1): 227-228.

[7] 王海福, 江增榮, 俞為民, 等. 桿式射流裝藥水下作用行為研究[J]. 北京理工大學學報, 2006, 26(3): 189-192.

WANG Hai-fu, JIANG Zeng-rong, YU Wei-min, et al. Behavior of Jetting Penetrator Charge Operating Underwater[J]. Transactions of Beijing Institute of Technology, 2006, 26(3): 189-192.

[8] 金鍵, 朱錫, 侯海量, 等. 水下爆炸載荷下艦船響應(yīng)與毀傷研究綜述[J]. 水下無人系統(tǒng)學報, 2017, 25(6): 396-409.

JIN Jian, ZHU Xi, HOU Hai-liang, et al. Review of Dynamic Response and Damage Mechanism of Ship Structure Subjected to Underwater Explosion Load[J]. Journal of Unmanned Undersea Systems, 2017, 25(6): 396-409.

[9] 尹群, 陳永念, 胡海巖. 水下爆炸研究的現(xiàn)狀和趨勢[J]. 造船技術(shù), 2003(6): 6-12.

YIN Qun, CHEN Yong-nian, HU Hai-yan. Present Situation and Trend of Underwater Explosion Research[J]. Journal of Marine Technology, 2003(6): 6-12.

[10] 宗智, 鄒麗, 劉謀斌, 等. 模擬二維水下爆炸問題的光滑粒子(SPH)方法[J]. 水動力學研究與進展A輯, 2007, 22(1): 61-67.

ZONG Zhi, ZOU Li, LIU Mou-bin, et al. SPH Simulation of Two-Dimensional Underwater Explosion[J]. Journal of Hydrodynamics (SerA), 2007, 22(1): 61-67.

[11] 曹兵. EFP戰(zhàn)斗部水下作用特性研究[J]. 火工品, 2007(3): 1-5.

CAO Bing. Study on the Performance of EFP Warhead Operating Underwater[J]. Initiators & Pyrotechnics, 2007(3): 1-5.

[12] 王長利, 周剛, 馬坤, 等. 聚能裝藥水下爆炸沖擊波載荷規(guī)律[J]. 高壓物理學報, 2017, 31(4): 453-461.

WANG Chang-li, ZHOU Gang, MA Kun, et al. Shockwave Characteristics of Shaped Charge Exploded Underwater[J]. Chinese Journal of High Pressure Physics, 2017, 31(4): 453-461.

[13] 凌榮輝, 錢立新, 唐平, 等. 聚能型魚雷戰(zhàn)斗部對潛艇目標毀傷研究[J]. 彈道學報, 2001, 13(2): 23-27.

LING Rong-hui, QIAN Li-xin, TANG Ping, et al. Target Damage Study of shaped-Charge Warhead of Antisubmarine Torpedo[J]. Journal of Ballistics, 2001, 13(2): 23-27.

[14] 胡功笠, 劉榮忠, 李斌, 等. 復(fù)合式魚雷戰(zhàn)斗部威力試驗研究[J]. 南京理工大學學報(自然科學版), 2005, 29(1): 6-8.

HU Gong-li, LIU Rong-zhong, LI Bin, et al. Power Experiments on Compound Torpedo Warhead Assembly[J]. Journal of Nanjing University of Science and Technology, 2005, 29(1): 6-8.

[15] 步相東, 王團盟. 魚雷聚能戰(zhàn)斗部自鍛彈丸水中運動特性仿真研究[J]. 魚雷技術(shù), 2006, 14(3): 44-47.

BU Xiang-dong, WANG Tuan-meng. Simulation Study on Kinematic Characteristic of Explosively Formed Projectile(EFP) in the Water for Torpedo Shaped Charge Warhead[J]. Torpedo Technology, 2006, 14(3): 44-47.

[16] LEE M, LONGORIA R G, WILSON D E. Cavity Dynamics in High-Speed Water Entry[J]. Physics of Fluids, 1997, 9(3): 540-550.

[17] 張世澤, 曹雄, 胡雙啟, 等. 錐形罩錐角對射流形成和侵徹影響的數(shù)值模擬[J]. 四川兵工學報, 2005, 26(6): 17-20.

ZHANG Shi-ze, CAO Xiong, HU Shuang-qi, et al. Numerical Simulation of the Influence of Cone Angle of Conical Cover on Jet Formation and Penetration[J]. Journal of Ordnance Equipment Engineering, 2005, 26(6): 17-20.

[18] 侯秀成, 楊金龍, 陳智剛, 等. 罩頂藥高影響藥型罩的射流轉(zhuǎn)化特性[J]. 彈箭與制導(dǎo)學報, 2017, 37(2): 67-70.

HOU Xiu-cheng, YANG Jin-long, CHEN Zhi-gang, et al. Study on Conversion Characteristics of Jet of Liner Influenced by Explosive Height above Liner[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2017, 37(2): 67-70.

[19] 張向榮, 黃風雷. 炸高對鎢銅射流空氣及水中侵徹的影響[J]. 北京理工大學學報, 2011, 31(3): 262-264.

ZHANG Xiang-rong, HUANG Feng-lei. Effects of Standoff on the Penetration of WCu-Pseudo-Alloy Shaped Charges in the Air and Water[J]. Transactions of Beijing Institute of Technology, 2011, 31(3): 262-264.

[20] 張會鎖, 趙捍東, 王芳, 等. 藥型罩壁厚對聚能射流影響的數(shù)值模擬[J]. 測井技術(shù), 2006, 30(1): 47-49.

ZHANG Hui-suo, ZHAO Han-dong, WANG Fang, et al. Numerical Simulation of Effect of Shaped Charge Jet on the Different Wall Thicknesses of Liners[J]. Well Logging Technology, 2006, 30(1): 47-49.

[21] 李明星, 王志軍, 伊建亞, 等. 不同聚能裝藥水下作用效果的對比分析[J]. 兵器材料科學與工程, 2017, 40(5): 95-99.

LI Ming-xing, WANG Zhi-jun, YI Jian-ya, et al. Comparative Analysis of Underwater Effect for Different Shaped Charge[J]. Ordnance Material Science and Engineering, 2017, 40(5): 95-99.

[22] 王海福, 江增榮, 李向榮. 藥型罩參數(shù)對聚能裝藥水下作用效應(yīng)的影響[J]. 北京理工大學學報, 2006, 26(5): 405-409.

WANG Hai-fu, JIANG Zeng-rong, LI Xiang-rong. Influences of Liner Parameters on the Effects of Shaped Charge Operating Underwater[J]. Transactions of Beijing Institute of Technology, 2006, 26(5): 405-409.

[23] 何洋揚, 龍源, 張朋軍, 等. 水介質(zhì)中組合式戰(zhàn)斗部成型技術(shù)研究[J]. 火工品, 2009(2): 52-56.

HE Yang-yang, LONG Yuan, ZHANG Peng-jun, et al. Study on Combined Shaped Charge Warhead Shaping in the Water[J]. Initiators & Pyrotechnics, 2009(2): 52-56.

[24] 李惠明, 陳智剛, 趙東華, 等. 典型戰(zhàn)斗部水下爆炸侵徹仿真[J]. 工程爆破, 2016, 22(4): 39-42.

LI Hui-ming, CHEN Zhi-gang, ZHAO Dong-hua, et al. Penetration Underwater Explosion by Typical Warhead Simulation[J]. Engineering Blasting, 2016, 22(4): 39-42.

[25] 吳晗玲, 段卓平, 汪永慶. 桿式射流形成的數(shù)值模擬研究[J]. 爆炸與沖擊, 2006, 26(4): 328-332.

WU Han-ling, DUAN Zhuo-ping, WANG Yong-qing. Simulation Investigation of Rod-Like Jets[J]. Explosion and Shock Waves, 2006, 26(4): 328-332.

[26] 周方毅, 詹發(fā)民, 吳曉鴻, 等. 圓錐-球缺藥型罩聚能戰(zhàn)斗部結(jié)構(gòu)優(yōu)化設(shè)計[J]. 爆破器材, 2014, 43(6): 43-47.

ZHOU Fang-yi, ZHAN Fa-min, WU Xiao-hong, et al. Optimal Design of Structure for Tapered and Spherical Combined Liner Shaped Charge Warhead[J]. Explosive Materials, 2014, 43(6): 43-47.

[27] 廖莎莎, 吳成, 畢世華. 鎢銅合金藥型罩水介質(zhì)中侵徹規(guī)律的實驗研究[J]. 彈箭與制導(dǎo)學報, 2011, 31(2): 99-101.

LIAO Sha-sha, WU Cheng, BI Shi-hua. A Penetrating Feature Investigation of Shaped Charge with Tungsten Copper Alloy Liner in Water[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2011, 31(2): 99-101.

[28] Livermore Software Technology Corporation. LS-DYNA Keywords User's Manual[K]. California: Livermore Software Technology Corporation, 2007.

[29] QIU Gang, HENKE S, GRABE J. Application of a Coupled Eulerian–Lagrangian Approach on Geomechanical Problems Involving Large Deformations[J]. Computers and Geotechnics, 2011, 38(1): 30-39.

[30] 朱錫, 張振華, 劉潤泉, 等. 水面艦艇舷側(cè)防雷艙結(jié)構(gòu)模型抗爆試驗研究[J]. 爆炸與沖擊, 2004, 24(2): 133-139.

ZHU Xi, ZHANG Zhen-hua, LIU Run-quan, et al. Experimental Study on the Explosion Resistance of Cabin near Shipboard of Surface Warship Subjected to Underwater Contact Explosion[J]. Explosion and Shock Waves, 2004, 24(2): 133-139.

Effects of Parameters of Charge Cover on Underwater Movement of Shaped Charge Jet

OU Jinga, YANG Yi-zhaoa, YU Wen-jiea, TANG Kuib

(a. Qian Xuesen College, b. National Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing 210094, China)

The work aims to study the effects of the structural parameters of the charge cover on the movement of shaped charge jet in water and improve the movement characteristics of shaped charge jet in water. Numerical simulation of water penetration by conical casing jet was carried out with the multi-material element ALE method, and the effects of cone angle and thickness on the movement parameters of shaped charge jet were analyzed. The results showed that the cone angle of conical cover and the thickness of charge cover had obvious effects on the shape, velocity and acceleration of shaped charge jet in water. The remaining velocity of the invading body after entering the water of 10 cm was analyzed, and it was found that in the process of increasing the cone angle of the charge cover from 30° to 150°, the remaining speed increased first and then decreased, and reached the fastest at 90°; when the thickness of the charge cover was between 1.5 mm and 4 mm, the remaining speed changed little, and when the thickness was between 4 mm and 6 mm, the remaining speed began to decrease significantly. According to the results of the numerical simulation in this work, when the cone angle is 90°, the jet formed by the charge cover with a thickness of 4 mm performs better in the water, the jet formed has the longest penetration depth, and the remaining speed after invading the aqueous medium of 10 cm is the largest, and the storage capacity is the strongest.

shaped charge jet; water medium; charge cover; motion characteristic; ALE method; penetration

2022-05-21；

2022-06-10

OU Jing (1999-), Male.

唐奎（1990—），男，博士。

Corresponding author：TANG Kui (1990-), Male, Doctor.

歐靖, 楊易釗, 郁文杰, 等. 藥型罩參數(shù)對聚能射流水下運動的影響[J]. 裝備環(huán)境工程, 2023, 20(2): 000-000.

O358

A

1672-9242(2023)02-0001-09

10.7643/ issn.1672-9242.2023.02.001

2022–05–21；

2022–06–10

國家自然科學基金（12102202）；南京理工大學本科生科研訓(xùn)練“百千萬”計劃（202010288157Z）

Fund：The National Natural Science Foundation of China (12102202); Nanjing University of Science and Technology Undergraduate Research training Project Named “Bai Qian Wan” (202010288157Z)

歐靖（1999—），男。

OU Jing, YANG Yi-zhao, YU Wen-jie, et al.Effects of Parameters of Charge Cover on Underwater Movement of Shaped Charge Jet[J]. Equipment Environmental Engineering, 2023, 20(2): 000-000.

責任編輯：劉世忠