梁 斌,馮高鵬,周婕群,魏雪婷,朱永清

(中國(guó)工程物理研究院 總體工程研究所,四川 綿陽(yáng) 621900)

大質(zhì)量鎢破片沖擊多層間隙蓋板炸藥數(shù)值模擬與試驗(yàn)研究

梁 斌,馮高鵬,周婕群,魏雪婷,朱永清

(中國(guó)工程物理研究院 總體工程研究所,四川 綿陽(yáng) 621900)

為分析大質(zhì)量破片對(duì)多層蓋板炸藥沖擊起爆閾值的影響因素,采用SPH和FEM相結(jié)合的方法開展了不同蓋板厚度和不同蓋板間隔情況對(duì)起爆閾值影響的數(shù)值模擬。在此基礎(chǔ)上,采用二級(jí)輕氣炮作為加載手段,進(jìn)行了大質(zhì)量鎢合金破片沖擊不同蓋板及間隙情況下炸藥沖擊起爆試驗(yàn)研究,獲得了起爆閾值。數(shù)值模擬與試驗(yàn)結(jié)果相一致,表明采用SPH與FEM相結(jié)合的方法可以較好地模擬此類問(wèn)題。通過(guò)擬合得到了大質(zhì)量破片沖擊起爆閾值速度與沖擊角度之間的經(jīng)驗(yàn)關(guān)系式。分析結(jié)果表明,大質(zhì)量破片沖擊產(chǎn)生的大量二次破片對(duì)起爆閾值速度和起爆延遲時(shí)間具有較大的影響。

爆炸力學(xué);沖擊起爆;數(shù)值模擬;光滑粒子法;起爆閾值

Abstract:In order to investigate the influence factor of large-mass fragments on shock initiation,the smoothed particle hydrodynamics(SPH)and finite element method(FEM)was combined,and the influence of different thinckness and different spacing of cover plates on initiation threshold was numerically simulated.The two-stage gas gun was applied,and the experiment of large-mass tungsten fragments impacting explosive was carried out.The critical initiation velocity was obtained.The simulation results are consistent with experimental results,and the SPH/FEM coupling method can be applied to analyze the velocity shock initiation problem.The correlations between critical initiation velocity and oblique impact angle were obtained by data fitting.Results show that the secondary fragments deriving from tungsten fragments have an important effect on the initiation velocity and the initiation delay-time.

Keywords:mechanics of explosion;shock initiation;numerical simulation;smoothed particle hydrodynamics;critical initiation velocity

裝藥沖擊安全性、反導(dǎo)戰(zhàn)斗部設(shè)計(jì)以及炸藥存儲(chǔ)安全性等工程問(wèn)題中,裝藥的沖擊起爆特性備受關(guān)注。反導(dǎo)戰(zhàn)斗部通常采用裝藥爆炸產(chǎn)生的高速預(yù)制/半預(yù)制小質(zhì)量破片(幾克至十幾克)引爆來(lái)襲導(dǎo)彈裝藥,實(shí)現(xiàn)對(duì)導(dǎo)彈的有效攔截。但隨著制導(dǎo)彈藥防護(hù)性能的增強(qiáng),傳統(tǒng)的反導(dǎo)彈藥產(chǎn)生的小質(zhì)量破片難以引爆來(lái)襲導(dǎo)彈,其反導(dǎo)性能受到嚴(yán)峻挑戰(zhàn)。據(jù)此,國(guó)外研制了新型的大質(zhì)量重金屬破片戰(zhàn)斗部,如愛國(guó)者-3(PAC-3)動(dòng)能欄截彈,其采用24枚200余克的破片作為毀傷元[1],成功攔截“飛毛腿”等戰(zhàn)術(shù)導(dǎo)彈。

對(duì)大多數(shù)戰(zhàn)術(shù)導(dǎo)彈戰(zhàn)斗部裝藥結(jié)構(gòu),可將此問(wèn)題簡(jiǎn)化為破片對(duì)具有一定厚度和間隙的多層蓋板裝藥的沖擊起爆問(wèn)題,涉及破片穿甲過(guò)程,殼體與炸藥在高溫、高壓、高應(yīng)變率下非定常的多相復(fù)雜反應(yīng)過(guò)程[2-6],從理論上研究該問(wèn)題存在較大的難度和挑戰(zhàn)性。目前研究者廣泛采用數(shù)值模擬并由少量試驗(yàn)對(duì)其進(jìn)行驗(yàn)證的方法進(jìn)行研究。當(dāng)前,國(guó)內(nèi)外學(xué)者對(duì)小質(zhì)量破片沖擊單層無(wú)間隙蓋板炸藥這類問(wèn)題進(jìn)行了大量研究[7-10],但是由于高速加載以及防護(hù)等困難,對(duì)于超過(guò)100 g的大質(zhì)量破片沖擊起爆研究較少[11],對(duì)于大質(zhì)量破片沖擊帶間隙的多層蓋板炸藥起爆問(wèn)題方面的研究鮮見公開報(bào)道。鑒于此,本文采用SPH(smoothed particle hydrodynamics)和FEM(finite element method)相結(jié)合的數(shù)值建模方法,運(yùn)用“升-降”法對(duì)大質(zhì)量鎢破片沖擊多層間隙蓋板炸藥起爆性能進(jìn)行了模擬。在此基礎(chǔ)上,采用二級(jí)輕氣炮作為加載手段,開展了大質(zhì)量鎢破片沖擊不同間隙蓋板炸藥起爆試驗(yàn),驗(yàn)證了數(shù)值模擬結(jié)果;獲得了影響大質(zhì)量破片沖擊多層間隔蓋板炸藥起爆閾值的主要因素。

1 數(shù)值建模方法與材料模型驗(yàn)證

對(duì)于破片沖擊起爆問(wèn)題,當(dāng)前廣泛采用的是Lee和Tarver[4-5]在Cochran方程的基礎(chǔ)上提出的反應(yīng)速率方程及點(diǎn)火—燃燒—快反應(yīng)三項(xiàng)形式的反應(yīng)速率方程[6-7],其計(jì)算結(jié)果與其PoP(progresses of a propagating)和壓力歷程曲線都符合較好,因此本文采用該模型研究沖擊起爆過(guò)程。高速?zèng)_擊過(guò)程中,破片和靶板將出現(xiàn)較大變形和破壞,單獨(dú)采用Lagrange、ALE(arbitrary langrange eluer)方法均難以重現(xiàn)整個(gè)物理過(guò)程,為此,本文采用有限元與光滑粒子相結(jié)合的方法對(duì)該問(wèn)題進(jìn)行建模與分析。

1.1 數(shù)值建模

本文主要基于文獻(xiàn)[9]試驗(yàn),并參考文獻(xiàn)[10-12]中數(shù)值模擬材料模型及參數(shù),建立破片沖擊帶蓋板炸藥的沖擊起爆模型,通過(guò)數(shù)值模擬與試驗(yàn)結(jié)果比對(duì),驗(yàn)證數(shù)值模型。依據(jù)試驗(yàn),建立不同直徑鎢珠沖擊3.65 mm厚的鉭蓋板及與之接觸50 mm厚的B炸藥。

建模中,鎢珠采用SPH粒子,粒子直徑為0.5 mm。蓋板、炸藥采用Lagrange網(wǎng)格。根據(jù)炸藥爆轟波寬度,網(wǎng)格劃分尺寸控制在0.5 mm以內(nèi)。鎢珠與蓋板、炸藥以及蓋板與炸藥之間采用接觸算法,根據(jù)對(duì)稱性,建立四分之一模型,在炸藥不同位置設(shè)置測(cè)量點(diǎn),如圖1所示。

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

炸藥的沖擊響應(yīng)采用Lee-Tarver三項(xiàng)式點(diǎn)火增長(zhǎng)模型模擬[4,13],該方程可很好地模擬非均勻炸藥的沖擊起爆特性:

(1)

式中:F是反應(yīng)率,即反應(yīng)炸藥與未反應(yīng)炸藥之比;a,b,c,d,G1,G2,I,u,s,w,e,f為炸藥參數(shù)。

采用JWL狀態(tài)方程描述炸藥的反應(yīng)物和未反應(yīng)物,該方程描述了壓力與相對(duì)體積V、單位體積以及初始能量E的關(guān)系:

(2)

式中:ω,R1,R2,A,B為炸藥常數(shù);E為CJ爆轟釋放的內(nèi)能。

本文采用高能B炸藥,主要參數(shù)如表1、表2[10,13-14]所示。

表1 B炸藥Lee-Tarver模型參數(shù)

表2 B炸藥JWL狀態(tài)方程參數(shù)

應(yīng)用Johnson-Cook強(qiáng)度模型和Gruneison狀態(tài)方程[14-15]描述鎢和坦材料。其主要參數(shù)取自AUTODYN材料庫(kù)[13]。

1.3 數(shù)值模擬模型驗(yàn)證

為模擬不同直徑D鎢珠對(duì)蓋板炸藥的沖擊起爆閾值,采用了“升降法”。鎢破片沖擊起爆過(guò)程中炸藥壓力與反應(yīng)率云紋如圖2、圖3所示,從圖中可以看出,破片高速?zèng)_擊情況下,在幾μs時(shí)炸藥局部反應(yīng)率就達(dá)到1,大約10 μs時(shí),炸藥峰值壓力就上升到29 GPa,隨后剩余炸藥很快反應(yīng)。單個(gè)破片沖擊帶蓋板B炸藥的起爆壓力閾值約為7 GPa,與文獻(xiàn)[2]中給出的起爆閾值一致。

單個(gè)和多個(gè)鎢球沖擊接觸蓋板炸藥的影響因素及炸藥反應(yīng)過(guò)程,具體見文獻(xiàn)[10,15]。

表3給出了不同直徑鎢珠沖擊蓋板炸藥起爆閾值。數(shù)值模擬計(jì)算得到的起爆閾值和文獻(xiàn)中沖擊起爆試驗(yàn)結(jié)果一致,驗(yàn)證了本文數(shù)值建模方法及模型參數(shù)的合理性。

表3 不同直徑鎢珠沖擊蓋板炸藥起爆閾值對(duì)比

2 大質(zhì)量鎢破片沖擊蓋板炸藥數(shù)值建模與分析

2.1 SPH與FEM相結(jié)合方法沖擊起爆數(shù)值建模

破片與蓋板在高速?zèng)_擊情況下存在較大的變形和破壞,其網(wǎng)格存在畸變,導(dǎo)致計(jì)算精度降低和計(jì)算時(shí)間步長(zhǎng)變小。為此,破片和玻璃鋼、鋁板和鋼板采用SPH粒子,炸藥采用Lagrange單元,SPH粒子直徑和炸藥單元尺寸均為0.5 mm。破片與靶體以及靶體之間采用接觸算法,為減小邊界的影響,炸藥側(cè)面和下表面采用無(wú)反射邊界條件。

對(duì)2種破片、4種靶體進(jìn)行建模,如表4、表5所示。表中,m1,l1,b1,h2分別為破片的質(zhì)量、長(zhǎng)、寬和高。由于破片著角θ超過(guò)80°后出現(xiàn)跳飛,因此著角θ計(jì)算范圍取值從0°到80°,間隔5°。

表4 破片結(jié)構(gòu)參數(shù)

表5 靶體結(jié)構(gòu)參數(shù)

圖4給出了鎢破片1(記為F1)和鎢破片2(記為F2)以著角θ=70°沖擊靶體4(記為T4)的數(shù)值計(jì)算模型,為了測(cè)取炸藥在沖擊過(guò)程中的響應(yīng),在炸藥中設(shè)置了不同測(cè)取點(diǎn),記為G1,…,G7。

2.2 多層間隔蓋板炸藥沖擊起爆數(shù)值模擬結(jié)果

采用升降法對(duì)不同情況沖擊過(guò)程進(jìn)行了數(shù)值模擬,根據(jù)起爆結(jié)果調(diào)整破片著靶速度,其調(diào)整速度最小間隔為25 m/s,直至獲得起爆和未起爆的速度值。圖5、圖6分別給出了破片F(xiàn)1以60°著角沖擊不同蓋板炸藥壓力與炸藥反應(yīng)率云圖。由圖可見,破片F(xiàn)1以1 475 m/s速度沖擊無(wú)間隙蓋板炸藥,在50 μs時(shí)刻炸藥起爆,而破片1以同樣角度沖擊3層間隔蓋板炸藥,需要2 025 m/s的速度才能起爆,從炸藥中監(jiān)測(cè)點(diǎn)可看出,其起爆所需時(shí)間約為70 μs,起爆時(shí)間較沖擊無(wú)間隙蓋板炸藥靶體情況長(zhǎng)。

圖7、圖8給出了破片F(xiàn)1以3 250 m/s速度、80°著角沖擊靶板T4的過(guò)程及監(jiān)測(cè)點(diǎn)壓力曲線,從圖7中可看出,高速?zèng)_擊過(guò)程中靶板破碎產(chǎn)生大量的二次破片,并與沖擊方向呈一定角度向四周飛散,其中一部分破片沖擊下層蓋板。盡管鎢破片在100.0 μs時(shí)已從第一層靶面跳飛,但大量二次破片對(duì)下層蓋板炸藥持續(xù)沖擊,致使炸藥在約110 μs時(shí)刻被引爆,見圖8。

圖9給出了鎢破片對(duì)不同蓋板炸藥沖擊起爆的著角與速度閾值。從圖中可看出,破片著角小于30°時(shí),著角對(duì)起爆閾值影響較小,當(dāng)著角超過(guò)60°時(shí),起爆閾值隨著角增加迅速增大,當(dāng)著角超過(guò)70°時(shí),起爆閾值的增加尤為突出,主要是因?yàn)榇藭r(shí)破片出現(xiàn)跳飛,起爆主要靠二次破片的累積和持續(xù)的沖擊效應(yīng)引爆炸藥。

靶板間隙對(duì)引爆閾值也有較大的影響,對(duì)于本文中計(jì)算條件,多層間隙蓋板炸藥的沖擊起爆閾值較蓋板炸藥接觸情況起爆閾值高20%～40%,其差距隨著角的增加而增大。在有間隙情況下,鎢破片開始階段沖擊上層靶板產(chǎn)生的沖擊波難以直接傳入炸藥中,只有在破片穿透復(fù)合材料與鋁靶后沖擊與炸藥接觸的蓋板產(chǎn)生的沖擊波才能對(duì)炸藥起爆產(chǎn)生直接作用。如果破片著角大于臨界跳飛角,則鎢破片將不能直接沖擊與炸藥接觸的蓋板,此時(shí)炸藥的起爆僅靠二次破片沖擊作用[11],因此,此時(shí)起爆閾值與無(wú)間隙情況差距進(jìn)一步增大。

3 大質(zhì)量破片沖擊多層間隔蓋板炸藥試驗(yàn)

3.1 試驗(yàn)布局

為驗(yàn)證上述計(jì)算結(jié)果以及分析彈靶交匯條件對(duì)起爆閾值的影響,設(shè)計(jì)了不同角度鎢破片沖擊間隔蓋板炸藥試驗(yàn)。破片及試驗(yàn)布局如圖10所示。破片通過(guò)二級(jí)輕氣炮發(fā)射并加速到預(yù)定的速度。蓋板炸藥試驗(yàn)件、見證板及壓力傳感器均通過(guò)安裝板安裝在密閉防爆容器中。壓力傳感器安裝在見證板前端,測(cè)量破片沖擊后爆炸容器中的壓力。

3.2 試驗(yàn)狀態(tài)判定

破片沖擊蓋板炸藥可能出現(xiàn)未反應(yīng)、燃燒、局部爆轟、完全爆轟等狀態(tài)。判斷依據(jù)如下:

①未反應(yīng)。壓力信號(hào)與典型應(yīng)力波振蕩衰減波形相似,試驗(yàn)后蓋板變形不明顯,靶室有大量殘余顆粒狀炸藥。

②爆燃。壓力信號(hào)無(wú)爆轟特征波形,試驗(yàn)后蓋板向前變形不明顯,側(cè)板有明顯燃燒痕跡,且側(cè)板與上掛板之間間隙有較多粉末狀殘余炸藥。

③局部爆轟。壓力波形在破片撞擊蓋板后緊隨出現(xiàn)峰值,試驗(yàn)后蓋板向前、側(cè)板向外變形較明顯,側(cè)板與頂板之間的間隙、側(cè)板上方有較多蠟狀殘余炸藥以及明顯的燃燒痕跡。

④完全爆轟。壓力波形和典型的爆轟波波形相似,試驗(yàn)后蓋板炸藥部分形成比炸藥直徑稍大的圓形破孔,側(cè)板向外嚴(yán)重變形,頂板及其支撐結(jié)構(gòu)變形隆起,靶室內(nèi)壁、頂板、側(cè)板上方無(wú)殘余炸藥。

具體判定過(guò)程中,壓力波形作為首要判斷依據(jù),然后根據(jù)試驗(yàn)后蓋板、側(cè)板等情況綜合分析進(jìn)行判定。

3.3 試驗(yàn)引爆閾值確定方法

破片撞擊蓋板炸藥的反應(yīng)十分復(fù)雜,由于諸如碰撞壓力,環(huán)境溫度、濕度,炸藥與蓋板結(jié)合緊密情況等多種因素的作用,使其爆炸程度在引爆閾值速度附近與撞擊速度并不是嚴(yán)格的單調(diào)線性關(guān)系。

引爆閾值速度計(jì)算原則:某一碰撞角下引爆閾值速度取局部爆炸狀態(tài)的最低速度與未爆狀態(tài)最高速度的平均值。若某碰撞角度下無(wú)未爆狀態(tài),則閾值取爆燃與局部爆炸相近速度的平均值。若某碰撞角度下既無(wú)未爆狀態(tài)又無(wú)局部爆炸狀態(tài),則其閾值取最高爆燃狀態(tài)速度。

3.4 沖擊起爆試驗(yàn)結(jié)果與分析

試驗(yàn)獲得16發(fā)有效數(shù)據(jù)。不同狀態(tài)下測(cè)試獲得的典型壓力歷程曲線如圖11所示。圖中,YB-14-ch4表示第14發(fā)試驗(yàn)中第4個(gè)傳感器所測(cè)數(shù)據(jù),其余依次類推。破片以70°著角沖擊間隙蓋板炸藥試驗(yàn)后見證板破壞以及蓋板炸藥如圖12和圖13所示。

70°碰撞角狀態(tài)共進(jìn)行6次試驗(yàn),數(shù)據(jù)采集成功5次,具體如表6所示。其中YB-14和YB-15試驗(yàn)壓力信號(hào)于2 000 μs附近與典型應(yīng)力波振蕩衰減波形相似(圖11(a)),試驗(yàn)后蓋板向前變形不明顯,靶室內(nèi)有殘余顆粒狀炸藥(圖12),狀態(tài)判斷為未爆。YB-10試驗(yàn)壓力波形在2 000 μs附近出現(xiàn)峰值(圖11(c)),試驗(yàn)后蓋板向前、側(cè)板向外變形較明顯,側(cè)板與頂板之間的間隙、側(cè)板上方有較多蠟狀殘余炸藥以及明顯燃燒痕跡(圖13),狀態(tài)判為局部爆炸。該狀態(tài)下引爆閾值速度按上述計(jì)算方法,取YB-10,YB-14發(fā)射速度平均值為2 360 m/s,綜合上述分析認(rèn)為該值為70°著角沖擊情況下的起爆閾值。

分析試驗(yàn)結(jié)果,著角θ分別為30°,60°,70°,80°下引爆閾值分別為1 716 m/s,2 085 m/s,2 420 m/s,3 150 m/s。對(duì)試驗(yàn)引爆閾值進(jìn)行指數(shù)方程擬合,獲得引爆閾值經(jīng)驗(yàn)曲線如圖14(a)。圖14(b)給出了數(shù)值模擬與試驗(yàn)引爆閾值的對(duì)比,并對(duì)數(shù)值模擬閾值進(jìn)行指數(shù)方程擬合,其中,擬合方程為y=A1exp(x/t1)+y0,R2=0.996,y0=1.510±0.034,A1=0.036±0.011,t1=21.199±1.704。

表6 破片以70°著角沖擊時(shí)試驗(yàn)狀態(tài)

從圖14可知,破片沖擊引爆閾值vp隨沖擊角度增大呈指數(shù)增大,當(dāng)沖擊角度超過(guò)70°后,引爆閾值急劇增加,這主要是由于超過(guò)一定角度后,破片沖擊出現(xiàn)跳飛情況,主要依靠沖擊產(chǎn)生的二次破片來(lái)引爆炸藥,其趨勢(shì)與前述數(shù)值模擬結(jié)果一致。從圖14(b)可看出,數(shù)值模擬獲得的起爆閾值與試驗(yàn)獲得的起爆閾值吻合較好,說(shuō)明本文采用的數(shù)值建模與分析方法具有一定的合理性,可用于此類問(wèn)題的模擬分析。

4 結(jié)束語(yǔ)

采用SPH與FEM相結(jié)合的方法建立了破片沖擊多層間隔蓋板炸藥數(shù)值模型,開展了破片在不同沖擊角度下的沖擊引爆試驗(yàn)研究,獲得以下初步結(jié)論:

①數(shù)值模擬獲得的引爆閾值與試驗(yàn)結(jié)果一致,說(shuō)明采用SPH與FEM相結(jié)合的方法可有效模擬大質(zhì)量破片對(duì)多層間隔蓋板炸藥沖擊起爆問(wèn)題;

②大質(zhì)量破片對(duì)多層間隔蓋板炸藥沖擊產(chǎn)生的二次破片對(duì)起爆閾值具有重要影響,其起爆時(shí)間相對(duì)較長(zhǎng),二次破片與大質(zhì)量破片共同引起炸藥的反應(yīng);

③大質(zhì)量破片對(duì)多層間隙蓋板炸藥的引爆閾值隨破片沖擊角度的增加呈指數(shù)增長(zhǎng)。當(dāng)著角超過(guò)70°時(shí),引爆閾值急劇增大,著角大于80°時(shí),因破片出現(xiàn)跳飛,炸藥難以引爆。

由于影響大質(zhì)量破片對(duì)多層間隙蓋板炸藥沖擊起爆閾值的因素較多,如炸藥配方、裝填工藝、內(nèi)部缺陷等,同時(shí)還受某些隨機(jī)因素的影響,很難對(duì)某類炸藥給出明確的起爆閾值,工程中往往給出特定條件下的起爆概率,其起爆機(jī)理及影響因素還需開展大量的研究工作。

致謝中國(guó)工程物理研究院化工材料研究所對(duì)試驗(yàn)裝藥、空氣動(dòng)力研究中心第四研究所對(duì)試驗(yàn)實(shí)施給予了大力支持,在此表示感謝。

[1] 孫亞力.不斷改進(jìn)的“愛國(guó)者”地空導(dǎo)彈系統(tǒng)[J].現(xiàn)代軍事,2000,30(2):27-29. SUN Ya-li.Improvement of “PAC” air-to-ground missile system[J].Modern Armament,2000,30(2):27-29.(in Chinese)

[2] CHEN Wei-dong,ZHANG Zhong,LIU Jia-liang.Numerical simulation and analysis of shock initation of shielded explosive impacted by fragments[J].Acta Armamentii,2009,30(9):1 187-1 191.

[3] 張寶平,張慶明,黃風(fēng)雷.爆轟物理學(xué)[M].北京:兵器工業(yè)出版社,2001:340-344. ZHANG Bao-ping,ZHANG Qing-ming,HUANG Feng-lei.Detonation physics[M].Beijing:Ordnance Industry Press,2001:340-344.(in Chinese)

[4] HELD M.Initiation phenomena of uncovered or covered high explosive charges by jets or projectile[C]//Proceedings of the 3rd International Autumn Seminar on Propellants,Explosive,and Pyrotechnics.Chengdu:Sichuan Publishing House of Science & Technology,1999:273-289.

[5] LEE E L,TARVER C M.Phenomenological model of shock initiation in heterogeneous explosives[J].Journal of Phys,Fluids,1980,23(12):2 362-2 372.

[6] CHEN J K,CHING H K,ALLAHDADI F A.Shock-induced detonation of high explosives by high velocity impact[J].Journal of Mechanics of Materials and Structures,2007,2(9):1 701-1 721.

[7] URTIEW P A,VANDER K S,TARVER C M,et al.Shock initiation experiments and modeling of composition B and C-4[C]//Proceedings of the 13th International Detonation Symposium.Norfolk,VA:Office of Naval Research,2006:55-61.

[8] 王樹山,李朝君,馬曉飛,等.鎢合金破片對(duì)屏蔽裝藥撞擊起爆的實(shí)驗(yàn)研究[J].兵工學(xué)報(bào),2001,22(2):189-191. WANG Shu-shan,LI Chao-jun,MA Xiao-fei,et al.An experimental study on the initiation of covered charge impacted by tungsten alloy fragments[J].Acta Armamentii,2001,22(2):189-191.(in Chinese)

[9] MURPHY M J,LEE E L.Modeling shock initiation in composition B[C]//The 10th International Detonation Symposium.Los Alamos,NM:Los Alamos National Laboratory,1993:73-84.

[10] 梁斌,馮高鵬,魏雪婷.多枚破片沖擊引爆帶蓋板炸藥數(shù)值模擬分析[J].彈箭與制導(dǎo)學(xué)報(bào),2013,33(6):43-47. LIANG Bin,FENG Gao-peng,WEI Xue-ting.Numerical simulation on shock initiation of composition explosive of cover board subjected to multi-fragment[J].Journal of Projectiles,Rockets,Missiles and Guidance,2013,33(6):43-47.(in Chinese)

[11] 趙海軍,盧永剛,梁斌.多層蓋板炸藥在鎢破片撞擊下起爆數(shù)值模擬研究[J].四川兵工學(xué)報(bào),2015(5):75-78. ZHAO Hai-jun,LU Yong-gang,LIANG Bin.Study on numerical simulation of multilayer plate explosive detonated by tungsten fragment impact[J].Journal of Sichuan Ordnance,2015(5):75-78.(in Chinese)

[12] LAWRENCE W,STARKENBERG J,KRZEWINSKI B.Secondary effects on projectile impact shock initiation[C]//Proceedings of the 13th International Detonation Symposium.Norfolk,VA:Office of Naval Research,2006:839-846.

[13] 崔凱華,洪滔,曹結(jié)東.射彈沖擊帶蓋板Comp B裝藥起爆過(guò)程數(shù)值模擬[J].含能材料,2010,18(3):286-289. CUI Kai-hua,HONG Tao,CAO Jie-dong.Numerical simulation of shock initiation in covered composition B by projectile impact[J].Chinese Journal of Energetic Materials,2010,18(3):286-289.(in Chinese)

[14] Century Dynamics Corporation.AUTODYN users manual[M].California:Century Dynamics Corporation,2005:234-249.

[15] 鄭平剛,趙鋒,張振宇.圓頭鋼彈對(duì)PBX-9404炸藥撞擊起爆的數(shù)值模擬[J].高壓物理學(xué)報(bào),2002,16(1):29-34. ZHENG Ping-gang,ZHAO Feng,ZHANG Zhen-yu.Numerical simulation of initiation of PBX-9404 impacted by round-nosed steel projectiles[J].Chinese Journal of High Pressure Physics,2002,16(1):29-34.(in Chinese)

NumericalSimulationandExperimentalStudyontheExplosiveWithMulti-layerIntervalTargetsImpactedbyLarge-massTungstenFragments

LIANG Bin,FENG Gao-peng,ZHOU Jie-qun,WEI Xue-ting,ZHU Yong-qing

(Institute of Systems of Engineering,China Academy of Engineering Physics,Mianyang 621900,China)

2016-12-10

國(guó)家自然科學(xué)基金項(xiàng)目(11672278)

梁斌(1976- ),男,高級(jí)工程師,博士,研究方向?yàn)槌Ｒ?guī)武器研制及毀傷評(píng)估。E-mail:lb_110119@163.com。

TJ038;TH213.3

A

1004-499X(2017)03-0056-07