王　超，吳國東，王志軍，許俊祥，孫加肖

(中北大學(xué) 機(jī)電工程學(xué)院，太原　030051)

?

中心孔護(hù)管對(duì)橢圓形藥型罩的數(shù)值模擬

王超，吳國東，王志軍，許俊祥，孫加肖

(中北大學(xué) 機(jī)電工程學(xué)院，太原030051)

摘要：為了獲得中心孔對(duì)橢圓形藥型罩的影響規(guī)律，采用動(dòng)力學(xué)有限元軟件AUTODYN-2D對(duì)不同中心孔直徑、護(hù)管厚度及起爆方式的橢圓形藥型罩成型裝藥進(jìn)行數(shù)值模擬研究；結(jié)果表明：當(dāng)中心孔直徑與橢圓形藥型罩直徑之比λ=0.25時(shí)，形成的侵徹體頭部速度最高、成型效果最佳；護(hù)管可以改變侵徹體的形態(tài)，但護(hù)管厚度d的變化對(duì)侵徹體的成型效果影響不大；中心孔橢圓形藥型罩的最佳起爆方式為正向環(huán)起爆。

關(guān)鍵詞：中心孔；護(hù)管；橢圓形藥型罩；侵徹體

本文引用格式：王超，吳國東，王志軍，等.中心孔護(hù)管對(duì)橢圓形藥型罩的數(shù)值模擬[J].兵器裝備工程學(xué)報(bào)，2016(7):60-62.

Citation format:WANG Chao, WU Guo-dong, WANG Zhi-jun, et al.Numerical Simulation of Oval Shaped Liners with A Center-Hole Protective Pipe[J].Journal of Ordnance Equipment Engineering,2016(7):60-62.

中心帶孔聚能裝藥[1]是一種為了滿足在藥型罩開孔工程上的需要同時(shí)能保證戰(zhàn)斗部威力的技術(shù)。中心開孔并帶有護(hù)管的聚能裝藥通常作為逆序起爆串聯(lián)戰(zhàn)斗部二級(jí)裝藥，作為單獨(dú)使用的戰(zhàn)斗部的研究國內(nèi)外還很少見。國內(nèi)雖有關(guān)于中心孔的文獻(xiàn)，但僅僅是關(guān)于中心孔直徑[2]對(duì)圓錐形藥型罩及球缺罩的研究，對(duì)護(hù)管厚度研究不多，同時(shí)沒有在相關(guān)文獻(xiàn)中查到有關(guān)起爆方式對(duì)中心孔藥型罩的研究。

對(duì)橢圓形藥型罩，研究其線性聚能裝藥[3]的研究比較多。但橢圓形藥型罩形成桿式侵徹體，頭部最大速度有限，通過改變中心孔的大小及護(hù)管的厚度，探討中心孔對(duì)其形成的侵徹體的頭部速度以及侵徹體形態(tài)的影響，為將來組合藥型罩以及多模戰(zhàn)斗部的研究提供理論依據(jù)。

1建立模型

1.1幾何模型和有限元模型

本文所研究的聚能裝藥結(jié)構(gòu)模型如圖1所示，藥型罩直徑為D=80 mm，裝藥高度為H=96 mm，殼體壁厚為2 mm，中心孔直徑為h，護(hù)管厚度為d，中心孔直徑與藥型罩直徑之比為λ。

圖1　聚能裝藥結(jié)構(gòu)模型

由于炸藥爆炸及藥型罩在壓垮過程中，材料變形比較大，因此主裝藥和藥型罩均采用Euler網(wǎng)格，計(jì)算模型中不考慮殼體對(duì)侵徹體形成的影響。由于該聚能裝藥結(jié)構(gòu)具有軸對(duì)稱性，本文采用有限元軟件AUTODYN-2D建立1/2有限元模型來模擬中心孔對(duì)藥型罩的影響，鏡像后的有限元模型如圖2所示。通過在空氣邊界上添加“Flow-out”邊界條件，消除邊界應(yīng)力影響，數(shù)值模擬中采用的單位制為mm-mg-ms。

圖2　鏡像后的有限元模型

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

本文藥型罩和護(hù)管材料均選用紫銅，均采用Shock狀態(tài)方程和Piecewise JC強(qiáng)度模型描述本構(gòu)關(guān)系。炸藥選用B炸藥，用JWL狀態(tài)方程描述，所有材料參數(shù)均來自AUTODYN軟件材料庫的內(nèi)置參數(shù)。

2中心孔橢圓形藥型罩侵徹體的形成過程

以λ=0.1為例，采用中心點(diǎn)起爆，分析侵徹體形成過程，如圖3所示。

首先護(hù)管被壓垮、變形，在t=10 μs時(shí)可以看出，護(hù)管形成的侵徹體被軸向擠出，隨著爆轟波的持續(xù)，橢圓形藥型罩頂部開始被擠壓變形，開始形成侵徹體，隨著時(shí)間的不斷增加，由于存在速度梯度，侵徹體頭部發(fā)生頸縮、拉斷，逐漸形成頭部速度較大的侵徹體，并且最終在侵徹體頭部形成少量的速度較高的粒子流，在t=60 μs時(shí)形成效果比較好，而在t=70 μs時(shí)整個(gè)侵徹體出現(xiàn)不同程度的頸縮。

3數(shù)值模擬結(jié)果及分析

3.1中心孔直徑對(duì)藥型罩的影響

分別取λ為0、0.1、0.15、0.2、0.25、0.3、0.35，來研究中心孔直徑對(duì)橢圓形藥型罩所形成的侵徹體頭部速度、有效射流長度等的影響。

圖3　侵徹體形成過程

由圖4可以看出，在裝藥結(jié)構(gòu)一定的情況下，中心孔直徑對(duì)形成的侵徹體的影響較大。當(dāng)λ≤0.25時(shí)，射流的頭部速度隨λ的增大不斷增加，在λ=0.25時(shí)，侵徹體頭部速度提升最大，最高可達(dá)4 382 m/s，提升幅度達(dá)到28.3%。當(dāng)λ>0.25時(shí)，隨著中心孔的不斷增大，藥型罩的質(zhì)量不斷減少，裝藥量也越來越少，造成射流的頭部速度不斷減小，同時(shí)也使得射流頭部形成少量的高速粒子流，速度最高可達(dá)6 000 m/s以上，但由于粒子的數(shù)量太少，忽略不計(jì)。侵徹體尾部速度變化不大，大致在400 m/s左右。

圖4　60 μs時(shí)不同λ值速度變化曲線

圖5　60 μs時(shí)不同λ值的射流長度

中心孔也會(huì)影響所形成的侵徹體射流長度。以45鋼為靶板，在侵徹過程中，侵徹截止速度大約為2 000 m/s，在此速度以上的射流為有效射流，從圖5可以看出，侵徹體的總長度和有效長度均在λ=0.25時(shí)最長，分別為188.5 mm和117.5 mm，漲幅可達(dá)34.2%和57.5%。由以上兩點(diǎn)可以看出，中心孔直徑對(duì)橢圓形藥型罩形成的侵徹體頭部速度、有效射流長度等都有明顯的影響。

3.2護(hù)管的厚度d對(duì)藥型罩的影響

改變護(hù)管厚度d的大小，分別取d為1，1.5，2，2.5，3，在60 μs時(shí)速度云圖如圖6所示。

隨著d的不斷增大，侵徹體的有效長度增加，但是速度梯度相差不大，除去高速粒子流，有效射流頭部速度相近，侵徹體前端的高速粒子流隨著d的增加，數(shù)量越來越多，并且頭部侵徹能力越來越強(qiáng)。生成的侵徹體的形態(tài)由橢圓形藥型罩時(shí)的桿式侵徹體變?yōu)榫勰苌淞鳌?/p>

3.3起爆方式對(duì)侵徹體形成的影響

由圖7可知，藥型罩起爆方式采用面起爆與正向環(huán)起爆時(shí)所形成的侵徹體的有效長度最長，由于速度梯度大，侵徹體頭部有少量的高速粒子流，而逆向環(huán)起爆，不形成高速粒子流，頭部形成微小的蘑菇頭，但起爆方式對(duì)侵徹體的形態(tài)沒有影響。

圖7　60 μs時(shí)不同起爆方式侵徹體形成圖

由于形成的高速粒子流太少，侵徹性能有限，因此頭部速度不考慮粒子流部分，侵徹體頭部與尾部速度如表1所示，正向環(huán)起爆下，藥型罩形成侵徹體頭部速度最高，相比點(diǎn)起爆，速度提升34.1%，尾部速度變化不大。綜合射流有效長度及其頭部速度可知，藥型罩聚能裝藥起爆方式采用正向環(huán)起爆效果最佳。

表1　不同起爆方式下侵徹體速度　m·s-1

4正交設(shè)計(jì)分析

通過正交設(shè)計(jì)方法對(duì)侵徹體速度在不同水平組合下進(jìn)行數(shù)值模擬，然后通過數(shù)據(jù)處理得出各因素的主次關(guān)系。

4.1因素水平標(biāo)的確定

以侵徹體的速度v為考察指標(biāo)，選取因素為起爆方式，中心孔直徑與藥型罩之比λ，護(hù)管厚度d，水平為3水平，經(jīng)前面討論，起爆方式選面起爆，正向環(huán)起爆，分別用1,2代替；λ選取 0.25,0.3；d選取2.5,3。因素水平如表2所示。

表2　因素水平

4.2標(biāo)頭設(shè)計(jì)與分析

選用正交表L4(23)安排計(jì)算，將3個(gè)因素依次放在L4(23) 的1、2、3列，根據(jù)正交表安排的4種方案依次進(jìn)行計(jì)算，并將結(jié)果填入右側(cè)，如表3所示。

表3　方案與結(jié)果

從極差排列出各因素對(duì)侵徹體速度影響的大小順序是：起爆方式、中心孔直徑h、護(hù)管厚度d，其中起爆方式是最主要的影響因素，不同的起爆方式對(duì)侵徹體速度的影響很大。

5結(jié)論

對(duì)橢圓形藥型罩而言，當(dāng)中心孔直徑與藥型罩直徑之比λ=0.25時(shí)，藥型罩形成的侵徹體頭部速度最高、有效射流長度最長、成型效果最好。

在λ=0.25時(shí)，隨著護(hù)管厚度d的變化，侵徹體頭部速度(除高速粒子流外)變化不大，但侵徹體的形態(tài)發(fā)生變化，由桿式侵徹體變?yōu)榫勰苌淞鳌?/p>

起爆方式不僅影響侵徹體的有效長度，還影響侵徹體頭部速度，采用正向環(huán)起爆時(shí)，侵徹體成型效果最佳。起爆方式是影響侵徹體速度的最主要因素。

參考文獻(xiàn)：

[1]高靖.一種新型聚能裝藥數(shù)值模擬與實(shí)驗(yàn)研究 [D].太原:中北大學(xué),2011.

[2]高靖.中心孔對(duì)聚能裝藥成形的影響[J].彈箭與制導(dǎo)學(xué)報(bào),2011,31(1)：110-112.

[3]武雙章，顧文彬.橢圓形罩線型聚能裝藥射流成型過程數(shù)值模擬研究[J].工程爆破,2013,19(5)：5-9.

[4]王成、付曉磊等.起爆方式對(duì)聚能射流性能影響的數(shù)值分析[J].北京理工大學(xué)學(xué)報(bào),2006,26(5)：401-404.

[5]王志軍,尹建平.彈藥學(xué)[M].2版.北京:北京理工大學(xué)出版社,2005.

(責(zé)任編輯周江川)

收稿日期：2016-01-29；修回日期：2016-02-29

作者簡介：王超(1988—)，男，碩士研究生，主要從事裝備理論與裝備技術(shù)研究。

doi:10.11809/scbgxb2016.07.014

中圖分類號(hào)：TJ413

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：2096-2304(2016)07-0060-04

Numerical Simulation of Oval Shaped Liners with A Center-Hole Protective Pipe

WANG Chao, WU Guo-dong, WANG Zhi-jun, XU Jun-xiang, SUN Jia-xiao

(College of Mechatronic Engineering, North University of China, Taiyuan 030051, China)

Abstract:In order to obtain the influence of the oval shaped liners with a center-hole, adopting the dynamic finite element software AUTODYN-2D, numerical simulation research with different central-hole diameter, protecting pipe thickness and the initiation of elliptic type medicine cover molding charge were processed. The results show that when lambda which the ratio of the diameter between the center-hole and the oval shaped liners is 0.25, the head speed of the penetration is highest, forming effect of the penetration is best; But the change of the thickness of the enlarge tube d had a little effect for shaping effect on penetration; The best way to detonation of oval shaped liners with a center-hole is positive ring initiation.

Key words:center-hole; protective pipe; oval shaped liner; penetration

【裝備理論與裝備技術(shù)】