王昕， 蔣建偉， 王樹(shù)有， 李梅

(北京理工大學(xué) 爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 北京 100081)

鎢球?qū)χ鎺ぱb藥的沖擊起爆數(shù)值模擬研究

王昕， 蔣建偉， 王樹(shù)有， 李梅

(北京理工大學(xué) 爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 北京 100081)

采用AUTODYN-3D數(shù)值軟件，開(kāi)展了鎢球與不同曲率半徑柱面帶鋼殼Comp B炸藥作用過(guò)程的數(shù)值模擬。分析了鎢球撞擊位置對(duì)炸藥沖擊起爆特性的影響，采用升降法獲得柱面帶殼裝藥的臨界沖擊起爆速度。結(jié)果表明：柱面帶殼裝藥沖擊起爆過(guò)程與平面帶殼裝藥相一致，炸藥起爆點(diǎn)發(fā)生在離炸藥和殼體界面一定距離處，且隨著速增加而越靠近交界面；相同條件下柱殼裝藥更易于起爆，其臨界起爆速度隨裝藥曲率半徑r增加呈現(xiàn)非線性增大，碰撞點(diǎn)偏移量δ=0時(shí)，r=∞的臨界起爆速度較r=40 mm時(shí)增加3.2%；隨偏移量δ的增大呈現(xiàn)指數(shù)增加，r=40 mm時(shí)，δ=0.94r的臨界起爆速度較δ=0時(shí)增加35.6%，較平面帶殼裝藥增加31.5%.

兵器科學(xué)與技術(shù)； 柱面帶殼裝藥； 起爆點(diǎn)； 臨界起爆速度

0 引言

對(duì)來(lái)襲彈藥進(jìn)行有效攔截毀傷是當(dāng)今彈藥技術(shù)研究的熱點(diǎn)問(wèn)題[1]。采用破片對(duì)其戰(zhàn)斗部高速侵徹實(shí)施引爆是最有效的手段，而問(wèn)題的實(shí)質(zhì)則是破片對(duì)帶殼裝藥的沖擊起爆。

國(guó)內(nèi)外學(xué)者針對(duì)這一問(wèn)題已進(jìn)行了廣泛的研究。從公布的文獻(xiàn)上來(lái)看，大多以破片沖擊平面帶殼裝藥模型為基礎(chǔ)，主要采用理論、數(shù)值模擬和試驗(yàn)方法開(kāi)展研究。包括不同材質(zhì)、形狀破片[2-5]對(duì)不同厚度殼體、炸藥[6]的沖擊起爆及起爆判據(jù)的建立[7-13]。然而，實(shí)際戰(zhàn)斗部多為圓柱形結(jié)構(gòu)，平面帶殼裝藥的沖擊起爆判據(jù)能否適用尚待驗(yàn)證。一些學(xué)者也曾開(kāi)展了相關(guān)試驗(yàn)和數(shù)值模擬分析，其中宋浦等[14]開(kāi)展某型破片式戰(zhàn)斗部對(duì)柱殼裝藥的撞擊毀傷試驗(yàn)研究，得到柱殼裝藥引燃、引爆的閾值范圍；辛建國(guó)等[15]等通過(guò)破片沖擊柱面薄殼裝藥試驗(yàn)，得到殼體和破片的機(jī)械作用是引爆薄殼裝藥的主要原因；江增榮等[1]則開(kāi)展破片材料對(duì)沖擊起爆影響的數(shù)值模擬研究，提出重金屬破片可提高反導(dǎo)戰(zhàn)斗部毀傷能力。但這些研究并未給出柱面帶殼裝藥和平面帶殼裝藥在起爆機(jī)理及起爆判據(jù)上的差異。

本文采用AUTODYN-3D數(shù)值軟件針對(duì)鎢球沖擊起爆柱面帶殼裝藥(簡(jiǎn)稱柱殼裝藥)問(wèn)題進(jìn)行數(shù)值模擬。通過(guò)分析鎢球以不同著速命中裝藥不同位置的沖擊起爆過(guò)程，揭示柱殼裝藥沖擊起爆機(jī)理，得到各參數(shù)對(duì)炸藥沖擊起爆閾值的影響規(guī)律。研究結(jié)果可為工程應(yīng)用提供參考。

1 數(shù)值模擬模型

1.1 物理模型

圖1為鎢球與柱殼裝藥作用的物理模型。其中鎢球質(zhì)量m；帶殼裝藥中炸藥曲率半徑r、長(zhǎng)度l，材料為Comp B炸藥；殼體厚度h，材料為鋼。為了簡(jiǎn)化計(jì)算，假設(shè)鎢球垂直于帶殼裝藥軸線入射。

圖1 鎢球與柱殼裝藥作用物理模型圖Fig.1 Physical model of tungsten fragment impacting on cylindrical charge with steel casing

圖2為鎢球碰撞點(diǎn)剖面示意及參數(shù)定義圖。其中O1、O2分別為帶殼裝藥軸心和鎢球中心；鎢球以著速v撞擊柱殼裝藥不同位置，定義碰撞點(diǎn)相對(duì)于O1的垂直距離為碰撞點(diǎn)偏移量，記為δ；O1O2與v夾角為撞擊角，記為θ，由幾何關(guān)系可得sinθ=δ/(r+h)；碰撞點(diǎn)處的殼體水平厚度記為hθ，即hθ=h/cosθ;沿碰撞點(diǎn)法線與切線方向?qū)⒅俜纸鉃榇怪彼俣葀v和滑移速度vs;在θ為0～75°范圍內(nèi)沿裝藥半徑方向每間隔7.5°設(shè)置一列觀測(cè)點(diǎn)，列號(hào)設(shè)為θ_1～θ_11，每列觀測(cè)點(diǎn)間間距4 mm.

圖2 鎢球碰撞點(diǎn)剖面示意及參數(shù)定義圖Fig.2 Impact point and parameter definition

1.2 數(shù)值模擬模型

采用AUTODYN-3D數(shù)值軟件建立離散化模型?？紤]模型對(duì)稱性，建立1/2模型以簡(jiǎn)化計(jì)算。圖3為建立的鎢球、殼體和炸藥的離散化網(wǎng)格模型，網(wǎng)格數(shù)量分別為2 048、2 880、192 000，采用Lagrange方法開(kāi)展計(jì)算。

圖3 離散化網(wǎng)格模型Fig.3 Discrete model

為了描述炸藥在沖擊作用下的起爆過(guò)程，其狀態(tài)方程采用Lee-Tarver狀態(tài)方程[16]：

dF/dt=I(1-F)b(μ-a)x+
G1(1-F)cFdpy+G2(1-F)eFgpz,

(1)

式中：F為燃燒質(zhì)量分?jǐn)?shù)，它在模擬爆轟過(guò)程中控制炸藥化學(xué)能的釋放；參數(shù)I和x為點(diǎn)火量沖擊強(qiáng)度及持續(xù)函數(shù)；μ為炸藥壓縮比；a為臨界壓縮度參數(shù)；b為點(diǎn)火項(xiàng)燃耗階數(shù)；參數(shù)G1和d為控制點(diǎn)火后早期增長(zhǎng)函數(shù)；c、y為燃燒項(xiàng)的燃耗階數(shù)和壓力冪數(shù)；參數(shù)G2和z為高壓反映率相關(guān)函數(shù)；e、g為常數(shù)；p為爆炸氣體壓力。炸藥未反應(yīng)物和反應(yīng)物均采用JWL狀態(tài)方程。Comp B炸藥模型參數(shù)取自AUTODYN材料庫(kù)，具體參數(shù)見(jiàn)表1.

表1 Comp B炸藥模型參數(shù)[17]

鎢球和殼體選用能較好描述材料大應(yīng)變、高應(yīng)變率及高溫度狀態(tài)的Johnson-Cook強(qiáng)度模型，材料強(qiáng)度模型、狀態(tài)方程和侵蝕準(zhǔn)則列于表2中。參數(shù)取自AUTODYN標(biāo)準(zhǔn)數(shù)據(jù)庫(kù)。

表2 鎢球、殼體材料模型

1.3 計(jì)算工況

為了研究鎢球?qū)χ鶜ぱb藥的沖擊起爆過(guò)程，確定命中位置、裝藥曲率半徑等參數(shù)對(duì)柱殼裝藥沖擊起爆的影響，設(shè)計(jì)了表3所示的計(jì)算工況，其中m為參考典型預(yù)制破片戰(zhàn)斗部中鎢球的質(zhì)量設(shè)定，依據(jù)常規(guī)彈藥戰(zhàn)斗部選取柱殼裝藥殼體厚度h和裝藥曲率半徑r，碰撞點(diǎn)偏移量δ則根據(jù)撞擊角θ分別為0°、15°、30°、45°、55°時(shí)的偏移距離制定。

采用“升降法”得到炸藥臨界起爆速度閾值，對(duì)比各工況炸藥起爆點(diǎn)位置、觀測(cè)點(diǎn)壓力變化，闡明柱殼裝藥與平面帶殼裝藥沖擊起爆的差異。

1.4 模型校驗(yàn)

材料參數(shù)的正確與否與計(jì)算結(jié)果正確性直接相關(guān)。采用經(jīng)典的皮克汀尼工程判據(jù)[9]對(duì)已建立的仿真模型進(jìn)行驗(yàn)證，表4為針對(duì)工況2～工況5選取兩種質(zhì)量鎢球撞擊起爆的臨界速度仿真值與理論值的對(duì)比。從表4可看出，仿真值與理論值誤差約5%，即認(rèn)定仿真算法的正確性。

表3 計(jì)算工況列表

表4 臨界起爆速度的仿真值與理論值對(duì)比

2 柱殼裝藥沖擊起爆過(guò)程

選取鎢球以碰撞點(diǎn)偏移量δ=0 mm對(duì)柱殼裝藥的沖擊起爆過(guò)程進(jìn)行詳細(xì)說(shuō)明。圖4(a)與圖4(b)、圖4(c)與圖4(d)分別為鎢球在兩種速度條件下撞擊柱殼裝藥和平面帶殼裝藥典型時(shí)刻應(yīng)力圖，圖5為觀測(cè)點(diǎn)壓力隨時(shí)間變化曲線。

圖4 鎢球撞擊柱形和平面帶殼裝藥典型時(shí)刻應(yīng)力圖Fig.4 Typical time-stress diagram of tungsten fragments impacting cylindrical and plate charges

圖5 炸藥內(nèi)觀測(cè)點(diǎn)壓力- 時(shí)間曲線Fig.5 Pressure-time curves of guages in charge

對(duì)于平面帶殼裝藥的沖擊起爆，在破片質(zhì)量和速度一定時(shí)，這類沖擊引爆機(jī)制與殼體厚薄程度有關(guān)。通常認(rèn)為當(dāng)殼體較薄時(shí)是沖擊起爆機(jī)制，而當(dāng)殼體較厚時(shí)是沖擊起爆和剪切起爆兩種機(jī)制并存[18]?？紤]到本文研究的殼體厚度，以沖擊起爆機(jī)制作為引爆機(jī)理開(kāi)展研究。

從圖4(a)、圖4(b)可以看出，在4 μs時(shí)刻內(nèi)，鎢球作用過(guò)程表現(xiàn)為對(duì)鋼殼的穿透和對(duì)裝藥的侵徹。v=2 826 m/s時(shí)，4 μs時(shí)刻距離殼體一定位置的炸藥處出現(xiàn)應(yīng)力階躍點(diǎn)，并迅速向外傳播，8 μs后形成穩(wěn)定傳播的爆轟波。而v=2 795 m/s時(shí)，沖擊波隨傳播距離逐漸衰減，炸藥未發(fā)生爆轟。對(duì)比圖4(b)、圖4(c)可以得到，柱殼和平面帶殼裝藥沖擊起爆過(guò)程基本一致，可確定柱殼裝藥的引爆機(jī)理同樣為沖擊起爆機(jī)制。

從圖5的觀測(cè)點(diǎn)壓力- 時(shí)間曲線可得到相同結(jié)論。由圖5(a)、圖5(b)可看出：在v=2 826 m/s條件下0～3 μs時(shí)間內(nèi)，第1列觀測(cè)點(diǎn)θ_1_1、觀測(cè)點(diǎn)θ_1_2變化趨勢(shì)及壓力值基本相同，在4 μs時(shí)刻觀測(cè)點(diǎn)θ_1_3壓力突然升高，隨后壓力繼續(xù)上升至一定峰值并穩(wěn)定傳播。此壓力的傳播即為爆轟的傳播，可確定起爆點(diǎn)在觀測(cè)點(diǎn)θ_1_3和觀測(cè)點(diǎn)θ_1_4附近。由于著速降低，圖5(a)中觀測(cè)點(diǎn)壓力逐漸降低，表明炸藥未發(fā)生穩(wěn)定爆轟。對(duì)比圖5(b)、圖5(d)可看出，v=2 826 m/s時(shí)平面帶殼裝藥內(nèi)觀測(cè)點(diǎn)壓力逐漸減小，炸藥未起爆。

圖6為不同著速下柱殼裝藥第1列觀測(cè)點(diǎn)壓力- 時(shí)間曲線。從圖6(a)中可看出，炸藥起爆點(diǎn)在觀測(cè)點(diǎn)θ_1_2和觀測(cè)點(diǎn)θ_1_3附近，圖6(b)顯示起爆點(diǎn)在觀測(cè)點(diǎn)θ_1_1和觀測(cè)點(diǎn)θ_1_2附近。與圖5(b)進(jìn)行對(duì)比，可得炸藥完全起爆時(shí)，起爆點(diǎn)隨著速增加逐漸向炸藥與殼體交界面靠近。

圖6 不同著速下炸藥內(nèi)觀測(cè)點(diǎn)壓力- 時(shí)間曲線Fig.6 Pressure-time curves of guages at different impact velocities

上述仿真結(jié)果可得：

1)柱殼裝藥和平面帶殼裝藥的沖擊起爆過(guò)程基本一致；

2)起爆點(diǎn)離裝藥和殼體交界面一定距離，隨著速增加，起爆點(diǎn)位置越靠近交界面。

3 柱殼裝藥的沖擊起爆影響因素分析

3.1 鎢球命中位置對(duì)柱殼裝藥沖擊起爆影響

3.1.1 命中位置對(duì)沖擊起爆現(xiàn)象的影響

對(duì)于平面帶殼裝藥，速度方向一定時(shí)，命中位置對(duì)撞擊角度并無(wú)影響；而對(duì)于柱殼裝藥，命中點(diǎn)位置不同意味著撞擊角的改變，沖擊起爆現(xiàn)象也將發(fā)生變化。

圖7為著速相同、命中殼體不同位置時(shí)典型時(shí)刻應(yīng)力圖。與圖4(a)相比，初始沖擊波形狀隨命中位置不同呈現(xiàn)非對(duì)稱性，壓力階躍點(diǎn)偏離碰撞點(diǎn)法線。碰撞點(diǎn)偏移量δ越大，炸藥起爆將越困難。

圖7 v=2 920 m/s鎢球命中不同位置時(shí)典型時(shí)刻炸藥沖擊起爆應(yīng)力圖Fig.7 Typical time-stress diagram of fragments impacting the different locations of cylindrical charge for v=2 920 m/s

3.1.2 命中位置對(duì)起爆點(diǎn)位置影響

提取計(jì)算文件中的觀測(cè)點(diǎn)壓力- 時(shí)間曲線以確定各工況炸藥起爆點(diǎn)位置變化。圖8為鎢球以偏移量δ=11.9 mm(θ=15°)撞擊柱殼裝藥時(shí)第2列、第3列、第4列觀測(cè)點(diǎn)壓力- 時(shí)間曲線。若起爆點(diǎn)位置位于第3列觀測(cè)點(diǎn)上，則沖擊波到達(dá)觀測(cè)點(diǎn)θ_4_2和觀測(cè)點(diǎn)θ_2_2的時(shí)間應(yīng)相同。而圖8顯示沖擊波觀測(cè)點(diǎn)到達(dá)順序?yàn)橛^測(cè)點(diǎn)θ_3_2、觀測(cè)點(diǎn)θ_4_2、觀測(cè)點(diǎn)θ_2_2，且觀測(cè)點(diǎn)θ_4_2壓力峰值高于觀測(cè)點(diǎn)θ_2_2，由此可確定起爆點(diǎn)位于第3列至第4列觀測(cè)點(diǎn)之間，即15°～22.5°之間。因此，撞擊角存在時(shí)起爆點(diǎn)偏離O1O2連線，位于O1O2和v夾角并距離殼體與炸藥界面一定距離。不同撞擊角下的起爆點(diǎn)位置也具有類似特征。

圖8 v=2 920 m/s和δ=11.9 mm時(shí)炸藥內(nèi)觀測(cè)點(diǎn)壓力- 時(shí)間曲線Fig.8 Pressure-time curves of guages for v=2 920 m/s and δ=11.9 mm

3.1.3 命中位置對(duì)臨界起爆速度的影響

通過(guò)“升降法”確定裝藥曲率半徑r=40 mm、殼體厚度h=6 mm條件下炸藥的臨界起爆速度vc，以偏移量δ=0時(shí)的臨界起爆速度vc0為基準(zhǔn)作歸一化處理。定義臨界起爆速度增量Δvc見(jiàn)(2)式，圖9為增量Δvc隨δ的變化關(guān)系。

(2)

由圖9可得：曲率半徑r=40 mm條件下，Δvc隨δ增大呈指數(shù)增加，δ=0.94r時(shí)的起爆速度較δ=0時(shí)增加35.6%，較平面帶殼裝藥增加31.5%，偏移量繼續(xù)增大將發(fā)生跳飛。因鎢球?qū)嶋H水平穿殼厚度hθ隨偏移量δ的增大而增大，穿透殼體消耗能量隨之增多，起爆炸藥也就變得困難。

圖9 臨界起爆速度增量隨偏移距離變化曲線Fig.9 Δvc vs. δ

為了消除碰撞點(diǎn)處殼體水平厚度hθ對(duì)炸藥沖擊起爆的影響，對(duì)hθ相同、偏移量δ不同時(shí)的炸藥臨界起爆速度進(jìn)行計(jì)算，并與相同殼體厚度、不同偏移量δ時(shí)的垂直速度vv增量進(jìn)行對(duì)比。圖10為hθ=6 mm時(shí)臨界起爆速度增量Δvc和h=6 mm時(shí)垂直速度增量Δvv隨δ的變化關(guān)系。從圖10可以看出：hθ=6 mm時(shí)Δvc隨δ增大而緩慢增加，但差距低于5%，證明即使保證hθ相同，命中位置依然對(duì)柱殼裝藥沖擊起爆產(chǎn)生影響；h=6 mm時(shí)垂直速度增量Δvv隨δ變化非線性遞減，差距至21%，說(shuō)明滑移速度vs對(duì)炸藥的沖擊起爆起積極作用。

圖10 臨界起爆速度和垂直速度增量隨偏移距離變化曲線Fig.10 Δvc and Δvv vs. δ

對(duì)鎢球命中殼體不同位置時(shí)炸藥起爆特性分析可得如下結(jié)果：

1)相同結(jié)構(gòu)條件下，柱殼裝藥的起爆隨碰撞點(diǎn)偏移量的增加而越困難；相同水平靶厚條件下，命中位置依然會(huì)影響柱殼裝藥沖擊起爆。

2)偏移量δ≠0時(shí)起爆點(diǎn)并不在O1O2連線上，而是位于O1O2和v的夾角并距離殼體與炸藥界面一定距離，起爆點(diǎn)位置不對(duì)稱。

3)垂直速度vv、滑移速度vs均會(huì)對(duì)炸藥起爆產(chǎn)生影響。

3.2 裝藥曲率半徑對(duì)柱殼裝藥沖擊起爆影響

取平面和柱殼裝藥(r=∞，r=40 mm)內(nèi)觀測(cè)點(diǎn)θ_1_2、觀測(cè)點(diǎn)θ_1_3、觀測(cè)點(diǎn)θ_1_4的壓力進(jìn)行比較。圖11為v=3 000 m/s時(shí)觀測(cè)點(diǎn)的壓力時(shí)程曲線。沖擊波到達(dá)平面帶殼裝藥的時(shí)間早于柱殼裝藥，但相同位置處柱殼裝藥內(nèi)的壓力值高于平面裝藥。

圖11 柱形和平面帶殼裝藥內(nèi)觀測(cè)點(diǎn)θ_1_2、觀測(cè)點(diǎn)θ_1_3、觀測(cè)點(diǎn)θ_1_4壓力時(shí)程曲線(v=3 000 m/s)Fig.11 Pressure-time curves of Gauges θ_1_2，θ_1_3，θ_1_4 in cylindrical and plate charges (v=3 000 m/s)

利用應(yīng)力波理論作定性解釋，鎢球撞擊柱殼裝藥時(shí)，會(huì)在殼體和炸藥界面處發(fā)生波系的反射和透射，在炸藥中傳入透射波。由于柱殼裝藥曲率半徑小于平面帶殼裝藥，傳入的波系會(huì)因曲率的存在而發(fā)生匯聚，相同撞擊條件下、炸藥內(nèi)同一位置處的壓力值便高于平面帶殼裝藥，因此，起爆柱殼裝藥易于起爆平面帶殼裝藥。

為了得到裝藥曲率半徑對(duì)柱殼裝藥臨界起爆速度的影響關(guān)系，選取鎢球撞擊r分別為40 mm、60 mm、75 mm、100 mm、∞(平板)的Comp B炸藥進(jìn)行仿真計(jì)算。以r=∞時(shí)的起爆速度為基準(zhǔn)做歸一化處理，獲得圖12臨界起爆速度增量隨裝藥曲率半徑倒數(shù)的關(guān)系。

圖12 臨界起爆速度增量隨曲率半徑倒數(shù)的變化關(guān)系(δ=0)Fig.12 Δvc vs. 1/r (δ=0)

由圖12可以看出：在偏移量δ=0時(shí)，臨界起爆速度增量隨1/r的增大非線性減??；r=∞(平板裝藥)時(shí)，裝藥臨界起爆速度最大，較r=40 mm時(shí)提高3.2%.

4 結(jié)論

采用AUTODYN-3D仿真軟件，對(duì)柱殼裝藥的沖擊起爆特性開(kāi)展研究，初步認(rèn)清了柱殼裝藥的沖擊起爆機(jī)理與規(guī)律，主要結(jié)論如下：

1)柱殼裝藥與平面帶殼裝藥起爆過(guò)程相一致：起爆點(diǎn)位于炸藥和殼體交界面一定距離。在臨界起爆速度以上時(shí)，起爆位置隨著速增加逐漸靠近交界面。相同條件下，起爆柱殼裝藥較易于起爆平面帶殼裝藥。

2)命中位置對(duì)柱殼裝藥沖擊起爆影響較大。臨界起爆速度隨偏移距離δ的增大呈指數(shù)增加，裝藥曲率半徑r=40 mm條件下，δ=0.94r時(shí)的起爆速度較δ=0時(shí)增加35.6%，較平面帶殼裝藥增加31.5%.

3)裝藥曲率半徑對(duì)柱殼裝藥的沖擊起爆有一定影響。臨界起爆速度隨曲率半徑增加呈現(xiàn)非線性增大，偏移距離δ=0條件下，r=∞時(shí)起爆速度較r=40 mm時(shí)提高3.2%.

References)

[1] 江增榮,李向榮,李世才，等.預(yù)制破片對(duì)戰(zhàn)斗部沖擊起爆數(shù)值模擬[J].彈道學(xué)報(bào),2009,21(1):9-13. JIANG Zeng-rong, LI Xiang-rong, LI Shi-cai, et al. Numerical simulation on shock initiation of performed fragment to warhead[J]. Journal of Ballistics,2009,21(1): 9-13.(in Chinese)

[2] 陳海利,蔣建偉,門建兵.破片對(duì)帶鋁殼炸藥的沖擊起爆數(shù)值模擬研究[J].高壓物理學(xué)報(bào),2006,20(1):109-112. CHEN Hai-li，JIANG Jian-wei，MEN Jian-bing. Numerical simulation of fragment impacting on charge with aluminum shell [J].Chinese Journal of High Pressure Physics,2006,20(1):109-112.(in Chinese)

[3] 李小笠,屈明,路中華, 等.三種破片對(duì)帶殼炸藥沖擊起爆能力的數(shù)值分析[J].彈道學(xué)報(bào),2009,21(4):72-75. LI Xiao-li, QU Ming, LU Zhong-hua, et al. Numerical analysis of impact initiation ability of three kinds of fragments on shell explosive [J]. Journal of Ballistics,2009, 21(4):72-75. (in Chinese)

[4] 童宗保,王金相,彭楚才，等.預(yù)制破片對(duì)屏蔽炸藥沖擊引爆研究[J].科學(xué)技術(shù)與工程,2013,14(7):173-177. TONG Zong-bao, WANG Jin-xiang, PENG Chu-cai, et al. Study onshock initiation of shielded explosive by prefabricated fragment [J]. Science Technology and Engineering, 2013,14(7):173-177. (in Chinese)

[5] 盧錦釗,智小琦,李娟娟，等.六棱柱破片沖擊起爆帶殼B炸藥速度閾值研究[J].彈箭與制導(dǎo)學(xué)報(bào),2016,36(1):82-84. LU Jin-zhao, ZHI Xiao-qi, LI Juan-juan, et al. Threshold velocity of hexagonal prism fragment impacting on Comp B with shell[J].Journal of Projectiles, Rockets, Missiles and Guidance, 2016,36(1):82-84. (in Chinese)

[6] 陳衛(wèi)東,張忠,劉家良，等.破片對(duì)屏蔽炸藥沖擊起爆的數(shù)值模擬與分析[J].兵工學(xué)報(bào).2009,30(9):1187-1191. CHEN Wei-dong, ZHANG Zhong, LIU Jia-liang, et al. Numerical simulation and analysis of shock initiation of shielded explosive by fragment [J].Acta Armamentarii, 2009,30(9):1187-1191. (in Chinese)

[7] Held M. Initiation criteria of high explosive at different projectiles or jet densities [J] .Propellants, Explosives, Pyrotechnics, 1996, 21(5):235-237.

[8] Roslund L A . Initiation of warhead fragments I：normal impacts, NOLTR73-124 [R]. Dahlgren, VA， US：Naval Surface Weapons Center,1973.

[9] Lloyd R M. Conventional warhead systems physics and engineering design[M].Reston, VA,USA: Progress in Astronautics and Aeronautics,1998:502-504.

[10] 方青,衛(wèi)玉章,張克明.射彈傾斜撞擊帶蓋板炸藥引發(fā)爆轟的條件[J].爆炸與沖擊,1997,17(2):154-158. FANG Qing，WEI Yu-zhang，ZHANG Ke-ming，et al．On the projectile oblique-impact initiation conditions for explosive covered with a plate [J]. Explosion and Shock Waves, 1997, 17 (2):154-158. (in Chinese)

[11] 張先鋒,趙有守,陳惠武.射彈沖擊引爆帶殼炸藥臨界條件[J].彈道學(xué)報(bào),2006,18(4)57-59. ZHANG Xian-feng, ZHAO You-shou, CHEN Hui-wu. Thecritical condition of shelled explosive initiated by projectile [J]. Journal of Ballistics, 2006,18(4)57-59. (in Chinese)

[12] 明路遙,力昌坤,程波，等.爆炸沖擊波作用下屏蔽裝藥模型起爆判據(jù)研究[J].四川兵工學(xué)報(bào),2006,37(4):95-97. MING Lu-yao, LI Chang-kun, CHENG Bo, et al. Study oninitiation criterion of shielding charge model under explosion shock wave [J]. Journal of Sichuan Ordnance, 2006,37(4):95-97. (in Chinese)

[13] 傅華,譚多望,李濤，等.鎢射彈引爆帶殼板炸藥閾值工程計(jì)算方法[J].含能材料,2008,16(1):100-102. FU Hua，TAN Duo-wang，LI Tao，et al．Engineering calculation on threshold velocity of covered explosives impacted by tungsten projectile [J].Chinese Journal of Energetic Materials, 2008, 16(1): 100- 102.(in Chinese)

[14] 宋浦,梁安定.破片對(duì)柱殼裝藥的撞擊毀傷試驗(yàn)研究[J].彈箭與制導(dǎo)學(xué)報(bào),2006,26(1)87-88. SONG Pu, LIANG An-ding. Experimental study on impact damage of fragment to cylinder charge [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2006,26(1):87-88. (in Chinese)

[15] 辛建國(guó),徐豫新,李超，等.破片沖擊柱面薄殼裝藥實(shí)驗(yàn)[J].兵工學(xué)報(bào),2014,35(2):222-227. XIN Jian-guo, XU Yu-xin, LI Chao, et al．Experiment of fragment impact on cylinder charge covered with thin shell [J]. Acta Armamentarii, 2014,35(2):222-227. (in Chinese)

[16] Lee E L, Tarver C M. Phenomenological model of shock initiation in heterogeneous explosive[J] Physics of Fluids, 1980, 23(12):2362.

[17] Murphy M J，Lee E L，Weston A M，et al. Modeling shock initiation in composition B[C]∥Proceedings of the 10th Detonation Symposium. Boston, MA, US: LLNL,1993.

[18] 賈憲振,陳松,楊建,等.雙破片同時(shí)撞擊對(duì)B炸藥沖擊起爆的數(shù)值模擬研究[J].高壓物理學(xué)報(bào), 2011,25(5):469-474. JIA Xian-zhen, CHEN Song, YANG Jian, et al. Numerical study of explosives initiation by simultaneous impact from two fragment [J]. Chinese Journal of High Pressure Physics,2011,25(5):469-474.(in Chinese)

Numerical Simulation on the Initiation of Cylindrical Covered Charge Impacted by Tungsten Sphere Fragment

WANG Xin, JIANG Jian-wei, WANG Shu-you, LI Mei

(State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China)

The interaction of tungsten fragment and cylindrical covered charge with different curvature radius is simulated by using AUTODYN-3D software. The influence of impacting position of fragment on the detonation characteristics of explosive is analyzed, and the “up-down” method is used to obtain the critical detonation velocity. The results show that the impact initiation progress of cylindrical covered charge coincides with that of plane covered charge. The detonation point is at a certain distance from the interface between explosive and casing, and gradually approaches to the interface as the fragment impact velocity increases. The cylindrical covered charge is easy to detonate under the same condition, and its critical detonation velocity increases nonlinearly with the increase in curvature radius. For collision point offsetδ=0, the critical detonation velocity is increased by 3.2% for infinite curvature radiusr=∞ compared with that forr=40 mm; the critical detonation velocity increases exponentially with the increase in offsetδ. Forr=40 mm, the critical detonation velocity forδ=0.94ris increased by 35.6% and 31.5% compared to that forδ=0 and plane covered charge, respectively.

ordnance science and technology; cylindrical covered charge; detonation point ;critical detonation velocity

2017-01-10

王昕(1990—)，女，博士研究生。E-mail:bitwangxin@bit.edu.cn

蔣建偉(1962—)，男，教授，博士生導(dǎo)師。E-mail:bitjjw@bit.edu.cn

TJ410.3+41

A

1000-1093(2017)08-1498-08

10.3969/j.issn.1000-1093.2017.08.006