基于SPH法的不同壁厚球缺型藥形罩桿式射流成型數(shù)值模擬

趙玉志 崔 瀚*

(沈陽工學(xué)院能源與水利學(xué)院,遼寧 撫順113122)

1 概述

桿狀射流是近幾年間在聚能裝藥方向上的熱點研究問題,它是通過采用某種藥型罩的結(jié)構(gòu)以及適當起爆方式,使壓垮藥型罩產(chǎn)生的聚能裝藥侵徹體在一定距離的飛行過程中被逐漸的拉長,但是桿式射流的斷裂時間相對滯后,這種聚能裝藥侵徹體形式具有較高的速度和較低速度梯度,能保持較大侵徹深度且有較大的侵徹孔徑[1]。桿狀射流由彈丸內(nèi)藥形罩形成,藥形罩的結(jié)構(gòu)直接影響影響桿狀射流狀態(tài)。

本文運用AUTODYN軟件中采用SPH法對不同壁厚結(jié)構(gòu)的球缺型藥型罩形成桿式射流進行仿真計算,最后進行結(jié)果輸出,得出不同壁厚的球缺型藥型罩形成桿式射流的過程數(shù)據(jù),并對其進行分析對比。

2 材料模型及參數(shù)

炸藥材料選為B炸藥,材料模型為None,狀態(tài)方程為JWL方程；藥形罩材料為CU-OFHC無氧銅,材料模型為Johnson Cook,狀態(tài)方程為Linear。無氧銅密度為8.9 g·cm-3,炸藥密度為1.619 g·cm-3。

球缺型藥型罩模型的幾何尺寸為:炸藥藥柱直徑取D1=55mm；炸藥藥柱高度取H1=54mm；球缺罩口部直徑取D2=52 mm；球缺罩高度取H2=17.2 mm；藥型罩中心壁厚為δ1,外側(cè)壁厚為δ2,聚能裝藥結(jié)構(gòu)如圖1所示,其中壁厚參數(shù)為控制變量。

圖1 藥型罩及炸藥結(jié)構(gòu)

3 仿真結(jié)果

3.1 等壁厚球缺型藥型罩桿式射流

此時取δ1=δ2=2mm,即藥形罩中心與外側(cè)厚度一致,經(jīng)過仿真不同時刻形成桿式射流形態(tài)如圖2所示,從圖中可以看出產(chǎn)生射流形狀較為均勻。提取的仿真結(jié)果如圖3所示,從圖3可知射流頭部速度在0.02 ms達到最大值,為3245.29 m/s,之后穩(wěn)定在2960±5m/s區(qū)間；尾部速度穩(wěn)定在1385±5m/s區(qū)間。中部速度在0.03 ms后穩(wěn)定在1740±5m/s區(qū)間。

圖2 等壁厚藥型罩桿式射流形成過程不同時刻狀態(tài)

圖3 等壁厚藥型罩桿式射流形成過程不同時刻速度變化曲線圖

3.2 中間厚外側(cè)薄球缺型藥型罩桿式射流

中間厚外側(cè)薄藥形罩取參數(shù)δ1=2.5 mm,δ2=1.5 mm,經(jīng)過仿真不同時刻形成桿式射流形態(tài)如圖4所示,從圖中可以發(fā)現(xiàn)產(chǎn)生的射流形狀頭部細,尾部粗,向軸線集中,射流頭部直徑約為4.1 mm,射流尾部直徑約為9.5 mm,射流中部直徑約為7.6 mm,射流長度約為235mm。提取的仿真結(jié)果如圖5所示,從圖中可知頭部速度在0.038 ms后達到最大值,為2136.30 m/s,之后穩(wěn)定在2050±5m/s區(qū)間；尾部速度在0.01 ms到0.03 ms之間速度波動較大,之后穩(wěn)定在800±20m/s區(qū)間；中部速度在0.03 ms后穩(wěn)定在630±5m/s區(qū)間。

圖4 中間厚外側(cè)薄藥型罩桿式射流形成過程不同時刻狀態(tài)

圖5 中間厚外側(cè)薄藥型罩桿式射流形成過程不同時刻速度變化曲線圖

3.3 中間薄外側(cè)厚球缺型藥型罩桿式射流

中間薄外側(cè)厚藥形罩厚度取值為δ1=1.5 mm,δ2=2.5 mm。經(jīng)過仿真不同時刻形成桿式射流形態(tài)如圖6所示,從圖中可以看出射流形狀分散,在0.05 ms后出現(xiàn)空心的情況,射流直徑約為6.5 mm,長度約為258mm。提取的仿真結(jié)果如圖7所示,從圖中可知,頭部速度在0.03 ms達到最大值,為2020.15 m/s,之后穩(wěn)定在1975±5m/s區(qū)間；尾部速度最大值1104m/s,隨后下降到250m/s左右；中部速度在0.03 ms后穩(wěn)定在975±5m/s區(qū)間。

圖6 中間薄外側(cè)厚藥型罩桿式射流形成過程不同時刻狀態(tài)

圖7 中間薄外側(cè)厚藥型罩桿式射流形成過程不同時刻速度變化曲線圖

3.4 綜合分析

從上述仿真結(jié)果可以看出,速度梯度對桿式射流的斷裂有著明顯的影響,具有的速度梯度越大,桿式射流產(chǎn)生斷裂的時間就越早,中間薄外側(cè)厚球缺型藥型罩形成的桿式射流頭部和尾部速度差以及速度梯度最大,在形成射流過程中最早產(chǎn)生斷裂現(xiàn)象,形成的射流最長,但是射流的密實性較差；中間厚外側(cè)薄球缺型藥型罩的頭部和尾部速度差同其速度梯度最小,產(chǎn)生斷裂的時間最晚,形成的射流長度適中；等壁厚球缺型藥型罩的頭部和尾部速度差和其產(chǎn)生的速度梯度適中,其產(chǎn)生斷裂的時間早于中間厚外側(cè)薄球缺型藥型罩,其具有的頭部速度和尾部速度是三種設(shè)計方案中最高的,而且其射流長度是最短的,三種方案的仿真對比數(shù)據(jù)如表1所示。

表1 球缺型藥型罩形成的桿式射流在0.1 ms時的主要參數(shù)表

4 結(jié)論

通過SPH方法對球缺型藥形罩形成桿式射流的過程進行數(shù)值模擬對比分析,等壁厚球缺型藥型罩形成的桿式射流速度最高,頭尾兩端速度差值最小,速度隨時間變化穩(wěn)定；能形成一定長度桿式射流,射流形狀為頭尾兩端大小較為均勻的柱狀。綜上所述,等壁厚球缺型藥型罩形成的桿式射流最為穩(wěn)定。