高 鵬,徐豫新,趙春龍,王樹山

(1 北京理工大學(xué)爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,北京 100081;2 中北大學(xué)機(jī)電工程學(xué)院,太原 030051)

動(dòng)能塊對(duì)屏蔽B炸藥沖擊引爆效應(yīng)研究*

高 鵬1,徐豫新1,趙春龍2,王樹山1

(1 北京理工大學(xué)爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,北京 100081;2 中北大學(xué)機(jī)電工程學(xué)院,太原 030051)

研究動(dòng)能塊對(duì)屏蔽B炸藥沖擊起爆機(jī)制與臨界起爆速度。采用AutoDyn-3D仿真軟件,對(duì)鎢合金動(dòng)能塊撞擊不同蓋板厚度屏蔽裝藥仿真計(jì)算,獲得著角0°～80°動(dòng)能塊引爆屏蔽裝藥臨界起爆速度,擬合獲得速度與著角函數(shù)關(guān)系。結(jié)果表明,臨界起爆速度在蓋板不同厚度條件下隨著角增大非線性變化。著角大于40°起爆閾值速度隨角度增加而增大,但增速不同;著角小于40°蓋板厚度大于10 mm速度隨著角增大而減小,蓋板厚度10 mm屏蔽裝藥臨界起爆速度具有隨機(jī)性。

爆炸力學(xué);沖擊起爆;動(dòng)能塊;臨界起爆速度

0 引言

對(duì)于大質(zhì)量動(dòng)能塊對(duì)屏蔽裝藥的沖擊起爆試驗(yàn),因加載條件多為二級(jí)輕氣炮,工作量極大且對(duì)于各物理量變化規(guī)律的直接測試是有難度的,因此目前大量采用數(shù)值仿真方法研究該問題。在仿真計(jì)算方法選取上,大部分采用拉格朗日算法和歐拉算法。近年來,隨著光滑粒子流體動(dòng)力法(SPH)的日益成熟,更多的爆炸問題使用了該種算法。時(shí)立國[5]采用SPH方法,模擬一維炸藥一端沖擊起爆的主要特性,驗(yàn)證了SPH方法在模擬爆轟問題方面的可行性。研究表明[6-7],SPH算法較傳統(tǒng)的拉格朗日(Lagrange)算法、歐拉(Euler)算法、ALE耦合算法等能更好的描述高速撞擊下材料的破碎行為,更適宜于高速條件下仿真分析。已有研究雖多,但針對(duì)動(dòng)能塊作用于屏蔽炸藥沖擊起爆的臨界起爆速度與蓋板厚度、著角關(guān)系的研究較少。

文中采用AUTODYN-3D軟件,使用SPH算法,建立了動(dòng)能塊對(duì)蓋板不同厚度下的屏蔽裝藥沖擊起爆數(shù)值仿真模型;據(jù)此,進(jìn)行了不同著角下,動(dòng)能塊對(duì)不同蓋板厚度下屏蔽裝藥沖擊起爆仿真,獲得了不同著角及蓋板厚度B炸藥起爆閾值,并分析了不同蓋板厚度起爆閾值隨碰撞角度變化規(guī)律。研究結(jié)果可對(duì)反導(dǎo)戰(zhàn)斗部設(shè)計(jì)提供技術(shù)基礎(chǔ)支持。

1 仿真計(jì)算模型

1.1 仿真模型離散化

建立仿真計(jì)算模型如圖1所示。單位制為cm-g-μs-GPa。動(dòng)能塊形狀選用常見的圓柱體動(dòng)能塊,大小為Φ16 mm×60 mm,屏蔽裝藥為不同厚度10#鋼板+Comp.B炸藥,炸藥周圍為鋁殼。

由于動(dòng)能塊及屏蔽裝藥模型均為對(duì)稱結(jié)構(gòu),為縮短計(jì)算時(shí)間,建立了動(dòng)能塊和屏蔽裝藥的1/2模型。采用Truegraid網(wǎng)格劃分軟件對(duì)兩模型劃分網(wǎng)格后,導(dǎo)入到AutoDyn軟件中仿真計(jì)算,網(wǎng)格尺寸1 mm。為觀察炸藥內(nèi)部壓力變化及判斷是否發(fā)生爆炸反應(yīng),在炸藥內(nèi)部設(shè)置了24個(gè)觀測點(diǎn),均位于1/2模型對(duì)稱平面內(nèi),每行觀測點(diǎn)距離0.5 cm,每行設(shè)置6個(gè)觀測點(diǎn),每個(gè)觀測點(diǎn)間隔3 cm,如圖1所示。

圖1 仿真計(jì)算模型

1.2 材料本構(gòu)模型及相關(guān)參數(shù)

根據(jù)已有小質(zhì)量破片對(duì)帶殼B炸藥沖擊起爆試驗(yàn)結(jié)果[8-10],驗(yàn)證獲得的數(shù)值仿真用材料模型后,進(jìn)行材料模型參數(shù)賦值。根據(jù)美國PAC-3[11]愛國者導(dǎo)彈,選擇圓柱體動(dòng)能塊,質(zhì)量為211 g,材料為鎢合金(W232),采用Shock狀態(tài)方程和Steinberg Guinan強(qiáng)度模型。

作為歐·亨利代表作品之一的《警察和贊美詩》早已進(jìn)入了研究者的視野。無論是小說的人物形象還是小說的作品主題，就連對(duì)其語言特色也早都有研究者進(jìn)行過分析?？上У氖牵谄湔Z言及相關(guān)問題的研究中，卻存在著以下諸多的局限性：

對(duì)于沖擊問題傳統(tǒng)采用拉格朗日和歐拉算法進(jìn)行求解,拉格朗日算法可更好的觀察到物質(zhì)點(diǎn)的運(yùn)動(dòng)過程,但對(duì)于高速?zèng)_擊時(shí)的大變形,會(huì)因網(wǎng)格畸變存在難以在不刪除畸變網(wǎng)格的情況下計(jì)算下去,歐拉算法雖能很好的解決這一問題,但無法很好的描述多物質(zhì)間的界面以及材料的固體失效斷裂行為。作為近20多年來逐步發(fā)展起來的SPH方法實(shí)質(zhì)是一種無網(wǎng)格方法,基本思想是將連續(xù)的流體(或固體)用相互作用的質(zhì)點(diǎn)組來描述,各個(gè)物質(zhì)點(diǎn)上承載各種物理量,包括質(zhì)量、速度等,通過求解質(zhì)點(diǎn)組的動(dòng)力學(xué)方程和跟蹤每個(gè)質(zhì)點(diǎn)的運(yùn)動(dòng)軌道,求得整個(gè)系統(tǒng)的力學(xué)行為;因此,它可避免極度大變形時(shí)網(wǎng)格扭曲而造成的精度破壞等問題,并且也能較為方便的處理不同介質(zhì)的交界面,特別適合于求解高速碰撞等動(dòng)態(tài)大變形問題。文中針對(duì)高速侵徹大變形特點(diǎn),對(duì)其使用光滑粒子流體動(dòng)力法(SPH)進(jìn)行計(jì)算。鋼板為STEEL 1006,使用Shock狀態(tài)方程和Johnson Cook強(qiáng)度模型,鋁殼為Al 7075-T6,采用Shock狀態(tài)方程和Steinberg Guinan強(qiáng)度模型。裝藥為B炸藥,采用Lee-Traver三項(xiàng)點(diǎn)火狀態(tài)方程[12]:

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

(1)

式中:F是燃燒質(zhì)量分?jǐn)?shù),它在模擬爆轟過程中控制著炸藥化學(xué)能的釋放;I、b、a、x、G1、c、d、y、G2、e、g、z為參數(shù)。參數(shù)a是阻止點(diǎn)火的臨界壓縮度;參數(shù)I和x控制著點(diǎn)火量,作為沖擊強(qiáng)度及持續(xù)的函數(shù);參數(shù)G1和d控制著點(diǎn)火后反應(yīng)的早期增長;參數(shù)G2和z確定高壓反應(yīng)率。各材料本構(gòu)模型及相關(guān)參數(shù)如表1所示。

表1 各材料本構(gòu)模型及相關(guān)參數(shù)

2 仿真計(jì)算結(jié)果及分析

2.1 計(jì)算結(jié)果

沖擊速度按“低-高”次序,共計(jì)進(jìn)行了36種工況沖擊起爆仿真計(jì)算,速度間隔50 m/s。因B炸藥的爆轟壓力為29.5 GPa,因此當(dāng)炸藥內(nèi)部壓力云圖或反應(yīng)度云圖出現(xiàn)爆轟波傳播,且炸藥內(nèi)觀測點(diǎn)壓力曲線形成脈沖峰,峰值超過20 GPa[13]時(shí),認(rèn)為炸藥發(fā)生了爆炸反應(yīng)。

如圖2、圖3所示,圖2中蓋板厚度為10 mm,動(dòng)能塊撞擊速度為1 300 m/s,壓力云圖中并未發(fā)生爆炸反應(yīng),觀測點(diǎn)處壓力峰值無明顯變化,最大峰值為2 GPa,可判斷裝藥在動(dòng)能塊撞擊下僅發(fā)生了爆燃反應(yīng),并未達(dá)到穩(wěn)定爆炸狀態(tài);圖3中動(dòng)能塊撞擊速度為1 350 m/s,壓力云圖中炸藥內(nèi)部產(chǎn)生了爆轟波陣面并向遠(yuǎn)處傳播,壓力曲線明顯形成壓力脈沖峰且峰值壓力超過20 GPa,表明裝藥發(fā)生了穩(wěn)定爆炸反應(yīng)。其他工況參考上述判斷方法進(jìn)行了仿真計(jì)算,結(jié)果如表2示,臨界起爆速度取未起爆速度與已起爆速度均值。

圖2 蓋板厚10 mm、碰撞角度70°、裝藥未起爆(速度1 300 m/s)的壓力云圖及觀測點(diǎn)壓力曲線

圖3 蓋板厚10 mm碰撞角度70°裝藥已起爆(速度1 350 m/s)

蓋板厚度/mm速度動(dòng)能塊碰撞角/(°)01020304050607080未起爆速度/(m/s)5507506509008508509501300150010已起爆速度/(m/s)600800700950900900100013501550臨界起爆速度/(m/s)57577567592587587597513251525未起爆速度/(m/s)950950950900900100010501300210020已起爆速度/(m/s)1000100010009509501050110013502150臨界起爆速度/(m/s)9759759759259251025107513252125未起爆速度/(m/s)13001050100010001050115012001550295030已起爆速度/(m/s)135011001050105011001200125016003000臨界起爆速度/(m/s)132510751025102510751175122515752975未起爆速度/(m/s)15001150110011001150130014002400390040已起爆速度/(m/s)155012001150115012001350145024503950臨界起爆速度/(m/s)152511751125112511751325142524253925

2.2 結(jié)果分析

由表2可見,當(dāng)動(dòng)能塊以大于40°撞擊屏蔽裝藥時(shí),裝藥臨界起爆速度隨蓋板厚度和著角的增加逐漸增大;但以小角度撞擊時(shí),蓋板厚度較小(10 mm)情況下臨界起爆速度大體上隨著角增大而增大,厚度較大(≥20 mm)情況下著角增大時(shí),起爆速度變化不大或逐漸減小。分析原因如下:

1)著角大于40°,隨蓋板厚度增大,動(dòng)能塊因高速侵徹過程中的侵蝕作用,質(zhì)量損失較大,欲引爆屏蔽裝藥,需提高速度,沖擊波傳播至炸藥形成熱點(diǎn),使裝藥發(fā)生穩(wěn)定爆炸反應(yīng)。速度提高產(chǎn)生了更高動(dòng)能,但大多損失在侵蝕蓋板上,因此在該種情況下動(dòng)能塊主要通過沖擊波引爆屏蔽裝藥。故大著角下臨界起爆速度隨著角穩(wěn)定單調(diào)增加。

2)著角小于40°,蓋板厚度大于等于20 mm時(shí),臨界起爆速度隨著角變化不大。臨界起爆速度與動(dòng)能塊作用炸藥的等效接觸面積有關(guān)。小于40°著角碰撞時(shí)動(dòng)能塊碎裂程度小,角度越大,接觸面積越大,比動(dòng)能越大。同時(shí),動(dòng)能塊侵徹厚蓋板產(chǎn)生崩落物降低了臨界起爆速度。此外,仿真所用圓柱動(dòng)能塊長徑比(3.75)與質(zhì)量(211 g)均較大,對(duì)厚蓋板具有更好的侵徹效果。

3)著角小于40°,蓋板厚度10 mm時(shí),臨界起爆速度隨角度增加而增大。分析該4種工況起爆過程發(fā)現(xiàn),由于蓋板較薄,動(dòng)能塊侵蝕蓋板后形狀質(zhì)量幾乎無損失,剩余速度和動(dòng)能很大,幾乎能夠以全部能量作用于炸藥,因此動(dòng)能塊主要依靠自身動(dòng)能侵徹蓋板后產(chǎn)生的應(yīng)力波持續(xù)壓縮炸藥產(chǎn)生熱點(diǎn),從而引爆裝藥,與蓋板較厚時(shí)規(guī)律不同。

4)當(dāng)著角為10°和30°,蓋板厚度10 mm工況下,出現(xiàn)了這兩種工況的臨界起爆速度高于同蓋板厚度相鄰角度計(jì)算結(jié)果的情況。分析兩種工況下的引爆過程發(fā)現(xiàn),動(dòng)能塊在侵蝕蓋板進(jìn)入裝藥后,對(duì)裝藥持續(xù)壓縮后才可引爆裝藥。可見,動(dòng)能塊撞擊屏蔽裝藥的反應(yīng)發(fā)展過程復(fù)雜,對(duì)于薄蓋板小撞擊角侵徹,前導(dǎo)沖擊波作用還是后續(xù)機(jī)械作用是難以區(qū)分的,引爆裝藥情況存在隨機(jī)性。

由于各工況間起爆機(jī)理不同,且動(dòng)能塊著角與蓋板厚度變化并非能夠直接等效成蓋板厚度的線性增加,所以難以采用傳統(tǒng)的Jacobs-Roslund(J-R)經(jīng)驗(yàn)公式進(jìn)行擬合。在此,根據(jù)表2中臨界起爆速度,采用較為簡單的三次多項(xiàng)式形式,擬合得到臨界起爆速度隨著角和蓋板厚度變化曲線,如圖4、圖5所示。

圖4 臨界起爆速度隨著角變化曲線

由圖4可見,蓋板厚度較薄時(shí),臨界起爆速度隨著角增大平緩增加,增速起初較緩慢,角度大于40°后加快;蓋板較厚時(shí),速度隨著角存在先增大后減小過程,且減小(增大)速度隨蓋板厚度增大而加快,在碰撞角度40°附近存在極小值點(diǎn),與上述原因分析一致。根據(jù)上述曲線,擬合得到臨界起爆速度隨碰撞角度在蓋板不同厚度條件下的函數(shù)表達(dá)式為:

(2)

式中:VCr為臨界起爆速度閾值(m/s);θ為撞擊角((°));h為蓋板厚度(mm)。

圖5 臨界起爆速度隨蓋板厚度變化曲線

由圖5可見,大碰撞角時(shí),起爆速度隨蓋板厚度增大穩(wěn)定增加,碰撞角為80°近似線性增加;小碰撞角時(shí),起爆速度隨蓋板厚度增大存在先增大后減小過程,在大厚度蓋板情況下較明顯,與上述原因分析一致。

3 結(jié)論

通過分析211 g動(dòng)能塊撞擊屏蔽B炸藥沖擊引爆效應(yīng)仿真計(jì)算,得到如下結(jié)論:

1)屏蔽B炸藥的蓋板厚度不同,起爆機(jī)理并不完全相同。著角小于40°時(shí)主要依靠動(dòng)能塊撞擊產(chǎn)生的機(jī)械作用引爆裝藥;大于40°時(shí)主要依靠透射沖擊波引爆裝藥。

2)臨界起爆閾值速度隨著角非線性變化,著角大于40°臨界起爆速度隨角度呈穩(wěn)定單調(diào)增加,小于40°臨界起爆速度變化規(guī)律與蓋板厚度有關(guān)。

3)蓋板大于20 mm情況下,引爆帶殼CompB炸藥的臨界起爆速度存在極小值點(diǎn),出現(xiàn)在著角40°附近。

[1] 南海陽, 韓曉明, 陳俊杰, 等. 一種新型的反TBM殺傷技術(shù)——定向含能動(dòng)能桿 [J]. 飛航導(dǎo)彈, 2013(9): 86-90.

[2] 江增榮, 李向榮, 李世才, 等. 預(yù)制破片對(duì)戰(zhàn)斗部沖擊起爆數(shù)值模擬 [J]. 彈道學(xué)報(bào), 2009, 21(1): 9-13.

[3] 梁斌, 馮高鵬, 魏雪婷. 多枚破片沖擊引爆帶蓋板炸藥數(shù)值模擬分析 [J]. 彈箭與制導(dǎo)學(xué)報(bào), 2013, 33(6): 62-69.

[4] 李俊龍, 李本勝, 王晶禹, 等. 反TBM戰(zhàn)斗部技術(shù)研究現(xiàn)狀及趨勢 [J]. 含能材料, 2015, 23(2): 189-198.

[5] 時(shí)立國, 湯文輝, 李楠. 一維炸藥沖擊起爆的SPH數(shù)值模擬研究 [C]∥第八屆全國沖擊動(dòng)力學(xué)學(xué)術(shù)討論會(huì), 2007: 222-228.

[6] 張忠, 陳衛(wèi)東, 陳浩. FEM/SPH耦合算法的射彈沖擊起爆仿真 [J]. 系統(tǒng)仿真學(xué)報(bào), 2009, 21(22): 7080-7083.

[7] 樂莉, 閆軍, 鐘秋海. 超高速撞擊仿真算法分析 [J]. 系統(tǒng)仿真學(xué)報(bào), 2004, 16(9): 1941-1943.

[8] 徐豫新. 破片毀傷效應(yīng)若干問題研究 [D]. 北京: 北京理工大學(xué), 2012.

[9] 馬曉飛, 李園, 徐豫新, 等. 破片對(duì)薄殼裝藥的沖擊起爆研究 [J]. 彈箭與制導(dǎo)學(xué)報(bào), 2009, 29(5): 133-138.

[10] 馬曉飛, 徐豫新, 李園, 等. 破片對(duì)模擬聚能戰(zhàn)斗部裝藥沖擊起爆研究 [J]. 戰(zhàn)術(shù)導(dǎo)彈技術(shù), 2010, 21(2): 29-32.

[11] 郭錫監(jiān), 高桂清. PAC-3反導(dǎo)系統(tǒng)作戰(zhàn)效能研究 [J]. 現(xiàn)代防御技術(shù), 2007, 35(2): 19-22.

[12] 陳海利, 蔣建偉, 門建兵. 破片對(duì)帶鋁殼炸藥的沖擊起爆數(shù)值模擬研究 [J]. 高壓物理學(xué)報(bào), 2006, 20(1): 110-111.

[13] 張寶坪, 張慶明, 黃風(fēng)雷. 爆轟物理學(xué) [M]. 北京: 兵器工業(yè)出版社, 2001: 160-161.

Study on Shock Detonation Effect of Kinetic Energy Block on Shielded Explosive B

GAO Peng1,XU Yuxin1,ZHAO Chunlong2,WANG Shushan1

(1 State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China; 2 School of Mechatronics Engineering, North University of China, Taiyuan 030051, China)

The critical velocity for ignition and shock initiation mechanism of kinetic energy block to shielded explosive of Comp. B was described. The tungsten-alloy kinetic energy block impacting different covered explosive was simulated by using Ansys AutoDyn-3D simulation software, and critical velocity for ignition of kinetic energy block detonating covered explosive at the angle of 0 to 80 was obtained. Through fitting, the functional relationship between velocity and collision angle was also obtained. The results showed that critical velocity for ignition changed nonlinearly with the increase of impact angle under conditions of different cover thickness. When impact angle was more than 40°, critical velocity for ignition of covered explosive increased with the increase of impact angle, but the growth rate was different. When impact angle was 40° or less and cover thickness was greater than 10 mm, critical velocity for ignition decreased with the increase of impact angle, and when the cover thickness was 10 mm, the critical velocity for ignition had characteristics of randomness.

explosion mechanics; impact initiation; kinetic energy block; critical initiation velocity

2015-11-29

高鵬(1990-),男,遼寧沈陽人,碩士研究生,研究方向:武器系統(tǒng)總體技術(shù)研究，通訊作者：徐豫新，講師，博士，E-mail:xuyuxin@bit.edu.cn。

TJ410.2

A