范俊余，方秦，張亞棟，陳力

(解放軍理工大學(xué) 工程兵工程學(xué)院，江蘇 南京210007)

0 引言

乳化炸藥是一類應(yīng)用廣泛的工業(yè)炸藥，其性能隨品種、配方乃至生產(chǎn)廠家的不同而有差異。TNT是一種常用的軍用炸藥，各種炸藥爆炸威力等參數(shù)計(jì)算以TNT 裝藥為爆源基礎(chǔ)。研究確定乳化炸藥的TNT 當(dāng)量系數(shù)，對(duì)于抗爆結(jié)構(gòu)設(shè)計(jì)、確定安全距離、爆炸事故分析等方面都有重要意義［1-3］。傳統(tǒng)上，常用炸藥的TNT 當(dāng)量系數(shù)都是根據(jù)能量相似原理按爆熱換算得到的，在常用炸藥手冊(cè)中均可查到。按爆熱換算主要基于做功能力的考慮，當(dāng)用于爆炸沖擊波計(jì)算時(shí)往往產(chǎn)生較大的誤差。為此，提出了由實(shí)際爆炸輸出來確定炸藥TNT 當(dāng)量的方法，即將同樣距離處得到相同爆炸峰值壓力(或正壓沖量)時(shí)的TNT 藥量與炸藥藥量的比值作為該炸藥的TNT 當(dāng)量系數(shù)［3］。文獻(xiàn)［4］采用這種方法對(duì)一種巖石型乳化炸藥的TNT 當(dāng)量系數(shù)進(jìn)行了試驗(yàn)研究，給出10～1 000 kPa 壓力范圍時(shí)該炸藥的等壓力TNT當(dāng)量系數(shù)平均值為0.708.一般情況下，用這種方法確定等壓力當(dāng)量系數(shù)與等沖量當(dāng)量系數(shù)時(shí)，二者均隨比例距離而變化，實(shí)際應(yīng)用時(shí)需加以注意。

本文采用模型試驗(yàn)和數(shù)值模擬相結(jié)合的方法，研究確定了2 號(hào)巖石乳化炸藥的等壓力TNT 當(dāng)量系數(shù)，為相關(guān)的研究和應(yīng)用提供了依據(jù)。

1 試驗(yàn)方案

試驗(yàn)分為空氣中自由場(chǎng)爆炸和模型通道外爆炸兩種情況。

1.1 空氣中自由場(chǎng)爆炸試驗(yàn)

采用標(biāo)準(zhǔn)TNT 藥塊和某2 號(hào)巖石乳化炸藥爆炸進(jìn)行對(duì)比試驗(yàn)，各炮藥量見表1.裝藥離地面1.5 m，在距裝藥中心50 cm、70 cm、90 cm 處各安裝1 只自由場(chǎng)壓力傳感器來測(cè)量爆炸沖擊波。

表1 空中爆炸試驗(yàn)藥量Tab.1 Explosive quantities for the explosion test in free air

1.2 模型通道外爆炸試驗(yàn)

圖1 炸藥位置示意圖Fig.1 Setup of explosion

如圖1所示，用1 cm 厚鋼板制作一個(gè)斷面寬×高為370 mm×360 mm、深為500 mm 的通道模型，置于大塊鋼板上，炸藥安裝在模型前面50 cm 處，各炮藥量見表2.分別用5 只壓電晶體壓力計(jì)和2 只壓阻式壓力計(jì)來量測(cè)模型底部封堵鋼板上的沖擊波壓力。壓力計(jì)安裝位置見圖2.

圖2 壓力傳感器位置示意圖Fig.2 Setup of pressure sensors

表2 通道外爆炸試驗(yàn)藥量Tab.2 Explosive quantities for the explosion test outside of corridor

2 試驗(yàn)結(jié)果及其分析

2.1 空氣中自由場(chǎng)爆炸試驗(yàn)

為建立巖石乳化炸藥與TNT 的等壓力藥量換算系數(shù)，將兩種炸藥的壓力試驗(yàn)數(shù)據(jù)分別按下式進(jìn)行擬合:

根據(jù)最小二乘法原理，對(duì)巖石乳化炸藥與TNT分別擬合出如下公式:

式中ΔpTNT、ΔpREE分別為TNT 和巖石乳化炸藥的沖擊波超壓峰值(MPa);分別為TNT 和巖石乳化炸藥的比例爆距

圖3是TNT 裝藥的試驗(yàn)數(shù)據(jù)及其擬合曲線與相關(guān)的經(jīng)驗(yàn)公式對(duì)比。從圖3可以看到，本文的試驗(yàn)數(shù)據(jù)與文獻(xiàn)［5］和文獻(xiàn)［6］吻合較好。巖石乳化炸藥的試驗(yàn)值和擬合曲線見圖4.

圖3 TNT 炸藥空氣中自由場(chǎng)爆炸的曲線Fig.3 Relation between peak pressure and scaling distance for the TNT explosions in free air

圖4 巖石乳化炸藥空氣中自由場(chǎng)爆炸的曲線Fig.4 Relation between peak pressure and scaling distance for the rock emulsion explosives explosions in free air

巖石乳化炸藥的TNT 當(dāng)量系數(shù)可表示為

式中mREE、mTNT分別為巖石乳化炸藥和TNT 的質(zhì)量(kg)。

文獻(xiàn)［4］給出了某巖石乳化炸藥的TNT 當(dāng)量系數(shù)與其比例爆距的關(guān)系圖。本文計(jì)算獲得曲線與文獻(xiàn)［5］給出的曲線的對(duì)比，見圖5，可以看出二者相差較大，這表明不同乳化炸藥的等壓力TNT 當(dāng)量系數(shù)差異顯著。

實(shí)際使用時(shí)，通常將某個(gè)比例爆距范圍內(nèi)的當(dāng)量系數(shù)平均值作為TNT 壓力當(dāng)量系數(shù)。對(duì)于本文研究情況，這一平均值為0.609.

圖5 巖石乳化炸藥等壓力TNT 當(dāng)量系數(shù)與比例距離關(guān)系(空氣中自由場(chǎng)爆炸)Fig.5 Relation between TNT equivalent coefficient of rock emulsion explosives and scaling distance (in free air)

2.2 模型通道外爆炸試驗(yàn)

由模型通道外爆炸試驗(yàn)結(jié)果擬合得到的兩種炸藥的爆炸沖擊波峰值計(jì)算公式分別為

將(5)式和各經(jīng)驗(yàn)曲線畫于圖6.從圖6看到，TNT 試驗(yàn)值比經(jīng)驗(yàn)值相差不大。巖石乳化炸藥的試驗(yàn)值和擬合曲線見圖7.

圖6 TNT 通道外爆炸時(shí)通道內(nèi)曲線Fig.6 Relation between peak pressure and scaling distance for TNT explosion outside of corridor

與前文類似，可以計(jì)算得到模型通道外爆炸情況下，比例爆距為1.698～3.24 m/之間時(shí)，巖石乳化炸藥的等壓力TNT 當(dāng)量系數(shù)的范圍為0.512～0.618，平均值為0.582.等壓力當(dāng)量系數(shù)隨比例距離的變化見圖8.

對(duì)比圖8與圖5看出，巖石乳化炸藥的當(dāng)量系數(shù)還與爆炸環(huán)境有很大關(guān)系。

圖7 巖石乳化炸藥通道外爆炸時(shí)通道內(nèi)曲線Fig.7 Relation between peak pressure and scaling distance for rock emulsion explosive explosions outside of corridor

圖8 巖石乳化炸藥等壓力TNT 當(dāng)量系數(shù)與比例距離關(guān)系(通道外爆炸)Fig.8 Relation between TNT equivalent coefficient of rock emulsion explosives and scaling distance(outside of corridor)

3 數(shù)值模擬

為評(píng)價(jià)試驗(yàn)結(jié)果的適用性，使用AUTODYN 有限元軟件進(jìn)行模擬分析。

3.1 材料模型

計(jì)算涉及到的材料模型主要有空氣、炸藥兩種材料模型?？諝獠牧夏Ｐ筒捎肁UTODYN 中的理想氣體狀態(tài)方程，模型的材料參數(shù)列于表3.

表3 空氣材料模型計(jì)算參數(shù)Tab.3 Calculating parameters of air material model

TNT 材料模型采用如下JWL 狀態(tài)方程:

式中:A、B、R1、R2、ω 是材料常數(shù);p 是壓力;Vr是相對(duì)體積;E0TNT是TNT 的初始內(nèi)能;D 是TNT 起爆速度。TNT 炸藥模型的參數(shù)取值如表4所示。

表4 TNT 材料模型計(jì)算參數(shù)Tab.4 Calculating parameters of TNT material model

圖9是49 g TNT 裝藥時(shí)距裝藥中心70 cm 處測(cè)點(diǎn)A1 的壓力時(shí)程曲線對(duì)比圖。從圖中可以看出，數(shù)值模擬曲線與試驗(yàn)曲線基本吻合。

3.2 空氣中自由場(chǎng)爆炸試驗(yàn)的數(shù)值模擬

空氣和炸藥采用Euler 單元，采用多物質(zhì)單元算法，采用2 mm×2 mm×2 mm 的網(wǎng)格大小。

圖9 A1 測(cè)點(diǎn)壓力時(shí)程曲線的數(shù)值與試驗(yàn)對(duì)比(mTNT=49 g)Fig.9 Comparison of pressure-time histories at measured point A1 between the numerical results and test data(mTNT=49 g)

圖10和圖11分別給出自由場(chǎng)中TNT 爆炸的壓力峰值和比沖量與比例爆距的關(guān)系圖。從圖中可以看出，比例距離小于2.0 m/時(shí)數(shù)值模擬得出的壓力峰值與試驗(yàn)測(cè)得數(shù)據(jù)和文獻(xiàn)［7］計(jì)算值基本吻合，但當(dāng)比例爆距大于2 m時(shí)超壓的數(shù)值模擬值比試驗(yàn)值和文獻(xiàn)［7］的經(jīng)驗(yàn)值略大。

由于缺少巖石乳化炸藥的狀態(tài)方程及其參數(shù)，本文取當(dāng)量系數(shù)0.609 將巖石乳化炸藥換算成TNT 炸藥，對(duì)半無限自由場(chǎng)中空爆試驗(yàn)進(jìn)行計(jì)算模擬。

圖10 曲線的數(shù)值與試驗(yàn)對(duì)比(空氣中自由場(chǎng)TNT 爆炸)Fig.10 Comparison ofcurves between the numerical results and test data (TNT explosion in free air)

圖11 比沖量與比例距離關(guān)系的數(shù)值與試驗(yàn)對(duì)比(空氣中自由場(chǎng)TNT 爆炸)Fig.11 Comparison of specific impulse versus scaling curves between the numerical results and test data (TNT explosion in free air)

圖12是70 g 巖石乳化炸藥爆炸時(shí)距裝藥中心70 cm 處測(cè)點(diǎn)A1 的壓力時(shí)程曲線圖。從圖中可以看出，數(shù)值模擬曲線與試驗(yàn)曲線基本吻合，數(shù)值模擬的壓力峰值比試驗(yàn)測(cè)得的壓力峰值略高。

圖13為數(shù)值模擬得出的巖石乳化炸藥壓力峰值與比例爆距的關(guān)系圖。從圖中可以看出，巖石乳化炸藥壓力峰值與比例爆距的變化規(guī)律的數(shù)值模擬結(jié)果與試驗(yàn)數(shù)據(jù)吻合較好。

3.3 模型通道外爆炸試驗(yàn)的數(shù)值模擬

通道外爆炸試驗(yàn)的數(shù)值分析模型見圖14所示。同樣，由于缺少巖石乳化炸藥的狀態(tài)方程及其參數(shù)，本文取當(dāng)量系數(shù)0.582 將巖石乳化炸藥換算成TNT炸藥，對(duì)通道外爆炸試驗(yàn)進(jìn)行計(jì)算模擬。模型通道采用剛體材料模型，炸藥和空氣采用歐拉單元模擬，所有單元均采用三維6 面體8 節(jié)點(diǎn)單元模擬，空氣與剛體之間采用流體—固體耦合的算法，取1/2 模型進(jìn)行分析。

圖12 A1 測(cè)點(diǎn)壓力時(shí)程曲線的數(shù)值與試驗(yàn)對(duì)比(mREE=70 g)Fig.12 Comparison of pressure-time histories at measured point A1 between the numerical results and test data (mREE=70 g)

圖13 的數(shù)值與試驗(yàn)對(duì)比(空氣中自由場(chǎng)巖石乳化炸藥爆炸)Fig.13 Comparison of curves between the numerical results and test data (rock emulsion explosives explosion in free air)

圖14 通道外爆炸的有限元計(jì)算模型Fig.14 Finite element model for explosions outside of corridor

圖15給出不同質(zhì)量的TNT 裝藥在通道爆炸時(shí)P6 測(cè)點(diǎn)壓力的試驗(yàn)與數(shù)值模擬結(jié)果的對(duì)比圖。

圖16給出了不同裝藥量巖石乳化炸藥通道外爆炸時(shí)P6 測(cè)點(diǎn)的試驗(yàn)壓力曲線與數(shù)值模擬結(jié)果的對(duì)比。

圖15 P6 測(cè)點(diǎn)壓力時(shí)程的數(shù)值與試驗(yàn)對(duì)比(TNT 在通道外爆炸)Fig.15 Comparison of pressure-time histories of P6 between the numerical results and test data (TNT explosion outside of corridor)

從圖15和圖16可以看出，不同藥量情況下，TNT 和巖石乳化炸藥通道爆炸條件下的數(shù)值模擬曲線與試驗(yàn)曲線基本吻合。這表明，在本文討論的范圍內(nèi)，將巖石乳化炸藥按當(dāng)量系數(shù)0.582 換算成TNT 對(duì)通道外爆炸情況進(jìn)行模擬計(jì)算，得到的壓力與比例距離關(guān)系曲線是可以接受的。

圖16 P6 測(cè)點(diǎn)壓力時(shí)程曲線的數(shù)值與試驗(yàn)對(duì)比(巖石乳化炸藥在通道外爆炸)Fig.16 Comparison of pressure-time histories of P6 between the numerical results and test data (rock emulsion explosives explosion outside of corridor)

4 結(jié)論

本文對(duì)TNT 和某2 號(hào)巖石乳化炸藥分別進(jìn)行了自由場(chǎng)和模型通道外爆炸試驗(yàn)研究，根據(jù)試驗(yàn)數(shù)據(jù)確定了巖石乳化炸藥的TNT 當(dāng)量系數(shù)，并進(jìn)一步對(duì)試驗(yàn)情況采用AUTODYN 有限元軟件進(jìn)行了數(shù)值模擬，主要結(jié)論有:

1)影響巖石乳化炸藥TNT 當(dāng)量系數(shù)的因素比較多，除了品種、生產(chǎn)廠家等因素外，爆炸試驗(yàn)條件也有明顯影響，實(shí)際應(yīng)用時(shí)應(yīng)注意加以區(qū)分。

2)在空氣中自由場(chǎng)爆炸條件下，當(dāng)巖石乳化炸藥的比例爆距為0.849～2.111 m/時(shí)，其等壓力TNT 當(dāng)量系數(shù)為0.522～0.656，平均值為0.609.

3)在文中通道外爆炸條件下，當(dāng)巖石乳化炸藥的比例爆距為1.70～3.24 m/時(shí)，其等壓力TNT當(dāng)量系數(shù)范圍為0.512～0.618，平均值為0.582.由于通道外爆炸與通道及距離條件有關(guān)，上述結(jié)論僅適用于本文的試驗(yàn)情況。

4)對(duì)于空氣中自由場(chǎng)爆炸和通道外爆炸兩種情況，將巖石乳化炸藥按所求得的當(dāng)量系數(shù)換算成TNT 裝藥進(jìn)行爆炸波參數(shù)的數(shù)值模擬計(jì)算是可行的。

References)

［1］ 東兆星，邵鵬.爆破工程［M］.北京:中國建筑工業(yè)出版社，2005.DONG Zhao-xing，SHAO Peng.Blasting engineering［M］.Beijing:China Building Industry Press，2005.(in Chinese)

［2］ Napadesky H S.TNT equivalency investigations，ADA007557［R］.US:DTIC，279-300.

［3］ TANG Ming-jun，PENG Jin-hua.A new conception on TNT equivalence［C］∥Proceedings of China-Japan Seminar on Energetic Materials.Nanjing:Nanjng University of Science and Technology，The University of Tokyo，1996:139-148.

［4］ 喬小玲，胡毅亭，彭金華，等.巖石型乳化炸藥的TNT 當(dāng)量［J］.爆破器材，1998，27(6):5-8.QIAO Xiao-ling，HU Yi-ting，PENG Jin-hua，et al.TNT equivalence of a rock emulsion explosive［J］.Explosive Materials，1998，27(6):5-8.(in Chinese)

［5］ 美國陸軍工程兵水道試驗(yàn)站.常規(guī)武器防護(hù)設(shè)計(jì)原理［M］.方秦，譯.南京:解放軍工程兵工程學(xué)院，1997:40-66.Waterway Experimental Station，Corps of Engineers.Fundamentals of protective design for conventional weapons［M］.FANG Qin，translated.Nanjing:The Engineering Institute of Engineer Corps，PLA，1997:40-66.(in Chinese)

［6］ Henrych J.The dynamics of explosion and its use［M］.Amsterdam:Elsevier，1979.

［7］ 李曉軍，張殿臣，李清獻(xiàn)，等.常規(guī)武器破壞效應(yīng)與工程防護(hù)技術(shù)［R］.洛陽:總參工程兵科研三所，2001.LI Xiao-jun，ZHANG Dian-chen，LI Qing-xian，et al.Damage effects of conventional weapons and engineering protection techniques［R］.Luoyang:The Third Research Institute of Engineering Corps，General Staff of PLA，2001.(in Chinese)