CL-20基含鋁炸藥爆轟波陣面法向速度與曲率的關(guān)系

沈　飛，王　輝，袁建飛，任新聯(lián)

(西安近代化學(xué)研究所，陜西西安710065)

?

CL-20基含鋁炸藥爆轟波陣面法向速度與曲率的關(guān)系

沈飛，王輝，袁建飛，任新聯(lián)

(西安近代化學(xué)研究所，陜西西安710065)

摘要：采用高速掃描相機及電探針測速法測量了具有相同鋁含量的CL-20基和RDX基含鋁炸藥的擬定態(tài)爆轟波形及爆速，分析了炸藥波陣面法向速度Dn與曲率κ之間的函數(shù)關(guān)系。結(jié)果表明，CL-20基含鋁炸藥的爆轟波陣面較RDX基含鋁炸藥的平坦，其法向爆速受曲率效應(yīng)的影響也較RDX基含鋁炸藥的小。當(dāng)κ>0.005mm-1時，其法向爆速的下降速率明顯小于RDX基含鋁炸藥；當(dāng)κ<0.005mm-1時，其法向爆速的下降速率略高于RDX基含鋁炸藥。

關(guān)鍵詞：爆炸力學(xué)；含鋁炸藥；非理想爆轟；Dspan(κ)關(guān)系；CL-20；高速掃描；擬定態(tài)波形

引言

含鋁炸藥是一類高密度、高爆熱和高威力的混合炸藥，由于其爆轟反應(yīng)區(qū)較寬，爆轟波的傳播過程中具有明顯的非理想特性，使得在研究爆轟波繞射、爆轟波與介質(zhì)相互作用、拐角效應(yīng)等方面遇到了較多困難[1]。對于這類問題，目前主要采用爆轟沖擊動力學(xué)(DSD)法進(jìn)行研究。DSD法可以對爆轟反應(yīng)流動方程組進(jìn)行解耦處理，即將其分解為爆轟沖擊波陣面的發(fā)展方程和一維擬定態(tài)反應(yīng)區(qū)動力學(xué)的常微分方程組[2-3]。在該方法的分析過程中需要了解爆轟波陣面曲率κ對法向爆速Dn的影響規(guī)律，即Dn(κ)關(guān)系式，由于該關(guān)系式只與炸藥的本構(gòu)性質(zhì)有關(guān)，所以一般需要通過曲面爆轟波傳播實驗進(jìn)行標(biāo)定。然而，這方面目前已有的成果主要是針對PBX-9502[4-6]、JB-9014[7]、JBO-9021[8]、RDX/TNT[9]、HMX/TNT[9]，對于含鋁炸藥Dn(κ)關(guān)系的研究還較少。

本研究采用高速掃描相機及電探針測速系統(tǒng)分別測量擬定態(tài)條件下CL-20基含鋁炸藥的爆轟波形及爆速，根據(jù)實驗結(jié)果分析了爆轟波陣面法向速度Dn與曲率κ之間的函數(shù)關(guān)系，并與常用的RDX基含鋁炸藥進(jìn)行對比，以期為該炸藥爆轟特性的深入研究提供參考。

1實驗

1.1樣品制備

CL-20基及RDX基含鋁炸藥配方見表1。采用壓裝工藝將兩種含鋁炸藥壓制成Φ50mm×50mm的藥柱，藥柱密度分別為1.96、1.75g/cm3。

表1　兩種含鋁炸藥配方

待測樣品由6節(jié)藥柱粘接而成，用JH-14作傳爆藥，藥柱尺寸為Φ25mm×25mm，密度為1.67g/cm3，實驗前將傳爆藥柱粘接在主裝藥一端的中心位置處。

1.2實驗裝置

采用SJZ-15型轉(zhuǎn)鏡式高速掃描相機掃描主裝藥柱端面的擬定態(tài)爆轟波形，實驗裝置如圖1所示。其中，主裝藥柱固定于V形槽內(nèi)，以減小邊界約束條件對爆轟波形的影響；掃描相機的光學(xué)狹縫通過反射鏡對準(zhǔn)主裝藥柱端面的直徑，并通過一組電探針測定炸藥的爆速。為了保證底片中波形的分辨率，相機的掃描速度設(shè)定為6mm/μs。

圖1　實驗裝置示意圖Fig.1　Schematic diagram of the test apparatus

2結(jié)果與分析

2.1擬定態(tài)波形

高速相機記錄的爆轟波形如圖2所示，橫向表示裝藥直徑，縱向表示掃描時間，從圖2可以清晰地看出波形的前沿。由于主裝藥的長徑比較大，可以使爆轟波演化為擬定態(tài)波形，此外，第一根電探針距離起爆端較遠(yuǎn)，此時爆速已基本穩(wěn)定，探針?biāo)鶞y得的爆速可認(rèn)為是擬定態(tài)爆速。

通過對底片進(jìn)行數(shù)字化判讀，并結(jié)合相機的掃描速度及底片的放大比，可得到爆轟波到達(dá)藥柱端面的時間，再將獲得的時間點與探針測量的擬定態(tài)爆速相乘，即可獲得波陣面曲線。目前，對于波陣面曲線擬合的函數(shù)一般選擇擬合精度較高的ln[cos(r)]級數(shù)，同時該函數(shù)也具有較高的數(shù)值計算效率，其具體表達(dá)式為[10]：

(1)

式中：r為波陣面上截面圓的半徑，mm；R為主裝藥的半徑，mm；z(r)為波陣面曲線，mm；ai及b均為擬合參數(shù)。

圖2　爆轟波形照片F(xiàn)ig.2　Photographs of detonation wave shape

在擬合過程中，正常n=2時可很好地模擬波陣面曲線。兩種炸藥爆轟波陣面曲線的擬合效果如圖3所示。

圖3　爆轟波波陣面擬合曲線Fig.3　Fitting curves of wave front

由圖3可以看出，CA-1炸藥擬定態(tài)爆轟波形的彎曲程度明顯小于RA-1炸藥，在裝藥邊界處，CA-1炸藥波陣面的z值約為0.8mm，而RA-1炸藥波陣面的z值接近1.2mm。所得擬合參數(shù)值及探針測得的擬定態(tài)爆速值(D0)均列于表2。

表2　波形擬合參數(shù)

2.2爆速與曲率的關(guān)系

當(dāng)爆轟波達(dá)到擬定態(tài)波形時，炸藥爆轟波法向傳播速度Dn與擬定態(tài)爆速D0之間的關(guān)系如圖4所示。

圖4　波陣面Dn與D0關(guān)系圖Fig.4　Relation of Dnand D0at wave front

根據(jù)圖4中的幾何關(guān)系，可得出：

(2)

式中：θ為波陣面法向與藥柱軸線方向的夾角。爆轟波陣面上的曲率κ用公式(3)計算：

(3)

當(dāng)r=0時，公式右端的兩項相等。由公式(1)~(3)可計算出爆轟波陣面的Dn(κ)關(guān)系曲線，結(jié)果如圖5所示。

圖5　兩種含鋁炸藥的Dn(κ)關(guān)系曲線Fig.5　Relation curves of Dn(κ) for the two explosives

由圖5可以看出，曲率κ從藥柱軸線至邊界兩側(cè)逐漸增大，且法向爆速Dn逐漸減小，CA-1炸藥波陣面的κ從0.004mm-1增加至0.01mm-1時，Dn從8.47mm/μs降至8.45mm/μs左右，κ繼續(xù)增大時，Dn的下降速度明顯提高；RA-1炸藥波陣面的κ從0.006mm-1增加至0.018mm-1時，Dn從7.91mm/μs降至7.86mm/μs左右，κ繼續(xù)增大時，Dn的下降速度明顯提高。此外，CA-1炸藥波陣面的最大曲率約為0.012mm-1，對應(yīng)的Dn較D0約降低0.03mm/μs，而RA-1炸藥波陣面的最大曲率約為0.023mm-1，對應(yīng)的Dn較D0約降低了0.07mm/μs。這是因為爆轟波陣面的彎曲程度反映了反應(yīng)區(qū)能量的損耗，從Dn的下降幅度可以看出爆轟傳播過程中沿著波陣面從軸線向兩側(cè)流動的能量大小[2,11]，因此，圖5中的數(shù)據(jù)表明，CA-1炸藥在爆轟傳播過程中沿波陣面從軸線向兩側(cè)流動的能量明顯小于RA-1炸藥。

由于難以通過理論分析獲得一般反應(yīng)速率形式下的Dn(κ)關(guān)系，通常采用經(jīng)驗公式對實驗數(shù)據(jù)進(jìn)行擬合，而簡單的線性函數(shù)又明顯不適用于描述圖5中的曲線，本研究采用一種效果較好的非線性函數(shù)形式，其具體形式為[5,10]：

(4)

式中：DCJ為炸藥的CJ爆速；C1、C2、C3、f、α、β、γ均為擬合參數(shù)。通過對圖5中的數(shù)據(jù)點進(jìn)行擬合，得到兩種炸藥Dn(κ)關(guān)系的擬合參數(shù)計算值，結(jié)果見表3。

表3　兩種含鋁炸藥Dn(κ)關(guān)系擬合參數(shù)計算值

為了進(jìn)一步比較Dn相對于DCJ的衰減狀況，設(shè)η=(DCJ-Dn)/DCJ，則公式(4)可變?yōu)椋?/p>

(5)

結(jié)合表3中參數(shù)，可得到η-κ關(guān)系曲線，如圖6所示。圖6可看出，當(dāng)κ小于0.005mm-1時，CA-1的η值略高于RA-1；當(dāng)κ大于0.005mm-1時，CA-1的η值明顯低于RA-1，且在κ達(dá)到最大值時，RA-1的η值接近1.2%，而此時CA-1的η值小于0.5%，這說明CA-1炸藥爆轟波在傳播過程中受曲率效應(yīng)的影響比RA-1小，這可能是由于CL-20的爆速和爆壓較高，使得含鋁炸藥爆轟波傳播過程中的非理想程度得以降低，但具體的影響規(guī)律還需要進(jìn)一步研究。

圖6　兩種含鋁炸藥的η-κ關(guān)系曲線Fig.6　Relation curves of η-κ for the two explosives

3結(jié)論

(1)與RDX基含鋁炸藥相比，CL-20基含鋁炸藥的爆轟波陣面明顯平坦，其波陣面的最大曲率約為0.012mm-1，RDX基含鋁炸藥波陣面的最大曲率約為0.023mm-1。

(2)CL-20基含鋁炸藥波陣面的法向爆速受到曲率效應(yīng)的影響較RDX基含鋁炸藥小，其法向爆速的最低值較擬定態(tài)爆速降低約0.03mm/μs，而RDX基含鋁炸藥法向爆速的最低值較擬定態(tài)爆速降低約0.07mm/μs。

(3)κ<0.005mm時，CL-20基含鋁炸藥法向爆速的衰減系數(shù)大于RDX基含鋁炸藥；而κ>0.005mm時，CL-20基含鋁炸藥法向爆速的衰減系數(shù)則小于RDX基含鋁炸藥。

參考文獻(xiàn)：

[1]裴明敬，田朝陽，胡華權(quán)，等. 鋁粉在溫壓炸藥爆炸過程中的響應(yīng)分析[J]. 火炸藥學(xué)報, 2013, 36(4): 7-12.

PEI Ming-jing, TIAN Zhao-yang, HU Hua-quan, et al. Response analysis of aluminum in the process of thermobaric explosive detonation[J]. Chinese Journal of Explosives and propellants, 2013, 36(4): 7-12.

[2]孫承緯, 衛(wèi)玉章, 周之奎. 應(yīng)用爆轟物理[M]. 北京: 國防工業(yè)出版社, 2000.

[3]Lambert D E, Stewart D S, Yoo S, et al. Experimental validation of detonation shock dynamics in condensed explosives[J]. Journal of Fluid Mechanics, 2006, 546: 227-253.

[4]Bdzil J B, Davis W C. Detonation shock dynamics(DSD) calibration for PBX 9502[C]∥Proceedings of Tenth Symposium (International) on Detonation. Boston: Office of Naval Research, 1993: 146-149.

[5]Aslam T D, Bdzil J B, Hill L G. Extensions to DSD theory: Analysis of PBX 9502 rate stick date[C]∥Proceedings of 11th International Detonation Symposium. Colorado: Office of Naval Research, 1998: 21-29.

[6]向梅, 饒國寧, 彭金華. 復(fù)合結(jié)構(gòu)裝藥爆轟波爆速與曲率關(guān)系的數(shù)值模擬[J]. 火炸藥學(xué)報, 2010, 33(4): 53-55.

XIANG Mei, RAO Guo-ning, PENG Jin-hua. Numerical simulation on the relationship between detonation velocity and curvature of composite charge structure ammunition[J]. Chinese Journal of Explosives and Propellants, 2010, 33(4): 53-55.

[7]譚多望, 鄒立勇, 張光升, 等. 低溫下JB-9014鈍感炸藥DSD參數(shù)研究[J]. 高壓物理學(xué)報, 2012, 26(4): 475-480.

TAN Duo-wang, ZOU Li-yong, ZHANG Guang-sheng, et al. Detonation shock dynamics calibration of JB-9014 explosive at low temperature[J]. Chinese Journal of High Pressure Physics, 2012, 26(4): 475-480.

[8]汪斌, 譚多望, 趙繼波, 等. 常溫下JBO-9021高能鈍感炸藥直徑效應(yīng)實驗[J]. 爆炸與沖擊, 2012, 32(5): 490-494.

WANG Bin, TAN Duo-wang, ZHAO Ji-bo, et al. Diameter effect of JBO-9021 rate sticks at room temperature[J]. Explosion and Shock Waves, 2012, 32(5): 490-494.

[9]張宏亮, 黃風(fēng)雷. RDX/TNT、HMX/TNT炸藥Dn(κ)關(guān)系研究[J]. 爆炸與沖擊, 2012, 32(5): 495-500.

ZHANG Hong-liang, HUANG Feng-lei. Study on theDn(κ) relation for the RDX/TNT and HMX/TNT rate sticks[J]. Explosion and Shock Waves, 2012, 32(5): 495-500.

[10]Bdzil J B. Steady-state two-dimension detonation[J]. Journal of Fluid Mechanics, 1981, 108: 195-226.

[11]趙繼波, 趙鋒, 譚多望, 等. 對鈍感炸藥Dn(κ)關(guān)系式一種改進(jìn)的探討[J]. 高壓物理學(xué)報, 2006, 20(3): 301-307.

ZHAO Ji-bo, ZHAO Feng, TAN Duo-wang, et al. Discussion on the improvement ofDn(κ) relation of IHE[J]. Chinese Journal of High Pressure Physics, 2006, 20(3): 301-307.

Relationship Between Normal Velocity and Curvature of Detonation Wave

Front for CL-20-based Aluminized Explosive

SHEN Fei, WANG Hui, YUAN Jian-fei, REN Xin-lian

(Xi′an Modern Chemistry Research Institute, Xi′an 710065, China)

Abstract:The steady state wave shapes and detonation velocities of CL-20-based and RDX-based aluminized explosives with same Al content were measured by a high speed scanning camera and electrical probe measuring velocity method. The function relationship between the normal velocity (Dn) and curvature (κ) of the detonation wave front of explosives was analyzed. Results show that the detonation wave front of CL-20-based aluminized explosive is flat than that of RDX-based aluminized explosive. The effect of curvature effect on normal detonation velocity of CL-20-based aluminized explosive is smaller than that of RDX-based aluminized explosive. Whenκ>0.005mm-1, the reduced rate ofDnof CL-20-based aluminized explosive is obviously less than that of RDX-based aluminized explosive. Whenκ<0.005mm-1, the reduced rate ofDnof CL-20-based aluminized explosive is slightly higher than that of RDX-based aluminized explosive.

Keywords:explosion mechanics; aluminized explosive; non-ideal detonation;Dn(κ) relation; CL-20; high speed scanning; steady wave shape

作者簡介:沈飛(1983-)，男，工程師，從事炸藥爆轟性能試驗與理論研究。

收稿日期:2014-07-06；修回日期:2014-08-17

中圖分類號：TJ55；O389

文獻(xiàn)標(biāo)志碼：A

文章編號：1007-7812(2015)01-0008-04

DOI：10.14077/j.issn.1007-7812.2015.01.002