HMX 基PBX 的作功能力及其JWL狀態(tài)方程

譚凱元，韓 勇，羅 觀，殷 明，黃毅民，盧校軍

（中國工程物理研究院化工材料研究所，四川 綿陽621900）

引 言

炸藥加速金屬的能力及其爆轟產(chǎn)物狀態(tài)方程是武器戰(zhàn)斗部設(shè)計中的重要問題。圓筒試驗則是考察炸藥作功能力的一種常用試驗方法，可以確定爆轟產(chǎn)物的JWL 狀態(tài)方程參數(shù)［1］。目前圓筒試驗已發(fā)展成為標(biāo)準(zhǔn)的系列試驗，國內(nèi)外研究人員利用圓筒試驗對各類炸藥的加速金屬能力和JWL 狀態(tài)方程做了大量工作，并獲得了其JWL 狀態(tài)方程參數(shù)［2-8］。

目前，我國已有的金屬加速炸藥還是以TNT為基的高固相熔鑄炸藥為主，這類炸藥的能量較高、金屬加速能力較強(qiáng)，但存在安全性較差的問題。目前，國外常規(guī)武器戰(zhàn)斗部裝藥已由以TNT為基的熔鑄炸藥逐步換裝為綜合性能優(yōu)異的PBX炸藥。與普通熔鑄炸藥相比，PBX 炸藥的安全性能較好，具有不敏感特征。此外，HMX 是現(xiàn)有單質(zhì)炸藥中綜合性能最好的高能炸藥，具有高爆速、高爆壓和較好的耐熱性、貯存穩(wěn)定性。本研究通過圓筒試驗結(jié)合數(shù)值計算技術(shù)評價一種新的HMX 基PBX 炸藥的作功能力并確定其JWL狀態(tài)方程參數(shù)，再利用平板飛片試驗來驗證JWL 狀態(tài)方程參數(shù)的有效性，為該炸藥今后的應(yīng)用和理論計算提供基礎(chǔ)數(shù)據(jù)。

1 實 驗

1.1 圓筒試驗

圓筒試驗裝置由高壓電雷管、傳爆藥柱（JO9159），待測HMX 基PBX炸藥（90%HMX，10%黏結(jié)劑，均為質(zhì)量分?jǐn)?shù)）、紫銅圓筒管、電探針、氬氣彈、光學(xué)窗口、高速轉(zhuǎn)鏡掃描相機(jī)等組成，如圖1所示。圓筒尺寸Φ25mm×300mm。在距引爆端約200mm 處用高速相機(jī)掃描記錄管壁的徑向膨脹過程，獲得膨脹距離對時間的試驗數(shù)據(jù)。在被測試炸藥第一和最后一個藥柱之間，放置兩對電探針，與時間記錄系統(tǒng)相連接，用于測量炸藥爆速。

圖1 圓筒試驗裝置示意圖Fig.1 Schematic diagram of the cylinder test

高壓電雷管、傳爆藥柱（Φ25mm×30mm），待測Φ25mm HMX基PBX 炸藥均為中國工程物理研究院化工材料研究所研制；SJZ-15 型轉(zhuǎn)鏡式高速掃描相機(jī)為蘇州光學(xué)儀器廠制造。

圓筒膨脹高速掃描底片如圖2所示。

圖2 圓筒壁膨脹過程的底片F(xiàn)ig.2 Photograph of the cylinder expansion process

圓筒膨脹距離R-R0與膨脹時間t之間的關(guān)系用公式（1）表示：

式中：a1、a2、a3、a4均為擬合系數(shù)，擬合的相關(guān)系數(shù)及爆速測試結(jié)果如表1所示。

表1 HMX 基PBX 炸藥擬合系數(shù)及爆速測試結(jié)果Table 1 Results of detonation velocity and fitting coefficient for the HMX based PBX

1.2 平面飛片試驗

HMX 基PBX 炸藥驅(qū)動平板飛片試驗裝置如圖3所示。

圖3 炸藥驅(qū)動平板飛片試驗裝置示意圖Fig.3 Schematic diagram of the explosive driven flyer plate test

從被測藥柱一端端面中心用雷管引爆炸藥，以驅(qū)動緊貼在藥柱另一端面（用真空脂排掉間隙空氣）的銅飛片運動，采用全光纖激光位移干涉測速系統(tǒng)（DISAR）對飛片的速度歷程進(jìn)行測量。DISAR 的激光探頭布置在藥柱中心軸線上，并對準(zhǔn)飛片中心以探測由運動飛片反射回來的帶有多普勒頻移信息的激光信號，通過對頻移信息的處理就可以計算出平板飛片的自由面速度歷史［9］。

HMX 基PBX 藥柱尺寸為Φ50mm×100mm，密度為1.721g／cm3，銅飛片為Φ50mm×1.0mm，密度為8.93g／cm3。DISAR 測速系統(tǒng)為電子科技大學(xué)研制的AFDISAR8000型測速裝置。

2 結(jié)果和討論

2.1 HMX 基PBX的作功能力

標(biāo)準(zhǔn)圓筒試驗通常用R-R0=5mm 和R-R0=19mm 時的圓筒壁速u和比動能Ecyl來描述一種炸藥的作功能力。由式（1）可推得：

將表1中的擬合系數(shù)代入式（2），即可通過式（2）和式（3）算得HMX基PBX的圓筒試驗特征值：

R-R0=5mm 時，u=1.410mm／μs；Ecyl=0.994kJ／g。

R-R0=19mm 時，u=1.615mm／μs；Ecyl=1.304kJ／g。

此外，格尼模型也是一種廣泛用于表征炸藥加速金屬能力的方法，在軸對稱膨脹情況下的格尼公式為［10］：

式中：m為炸藥質(zhì)量；M為圓筒質(zhì)量；u為圓筒壁的終速度，可由式（2）算得為1.658mm／μs，為表征炸藥加速金屬能力的格尼系數(shù)（又稱格尼速度），由式（4）最終計算出HMX 基PBX 的格尼速度為2.8mm／μs。

2.2 HMX 基PBX 的JWL狀態(tài)方程

JWL狀態(tài)方程的形式為［11］：

式中：p為爆轟產(chǎn)物的壓力；V為爆轟產(chǎn)物的相對比容。A，B，R1，R2，ω和E0為6個待定參數(shù)，由圓筒試驗來確定。

采用非線性有限元動力學(xué)軟件ANSYS／LSDYNA 對圓筒試驗進(jìn)行數(shù)值計算，通過對圓筒壁膨脹過程的試驗值和計算值進(jìn)行比較，確定JWL狀態(tài)方程參數(shù)。計算時對圓筒建立軸對稱模型，如圖4所示，OABC為HMX 基PBX炸藥，采用基于JWL狀態(tài)方程的MAT＿HIGH＿EXPLOSIVE＿BURN 模型；CDEF為銅管壁，采用GRUNEISEN 狀態(tài)方程［12-13］和MAT＿JOHNSON＿COOK 力學(xué)模型［14］；由A點起爆。

圖4 圓筒試驗數(shù)值模擬軸對稱模型Fig.4 Axisymmetric model of numerical simulation for the cylinder test

在圓筒試驗數(shù)值計算中，紫銅參數(shù)由文獻(xiàn)［11］和［13］獲得，炸藥JWL狀態(tài)方程參數(shù)在開始時輸入一組已知的與其組分相近的炸藥參數(shù)，得到圓筒膨脹距離與時間的關(guān)系曲線，再與試驗測定的膨脹距離-時間曲線進(jìn)行比較，根據(jù)它們之間的差異對參數(shù)進(jìn)行調(diào)整，重新輸入程序進(jìn)行計算，再把計算結(jié)果與試驗結(jié)果比較，調(diào)整參數(shù)值后進(jìn)行下一次計算。這樣經(jīng)過多次反復(fù)計算，不斷修正狀態(tài)方程參數(shù)，直到計算的膨脹距離-時間曲線與試驗曲線間的誤差小于1%時，即認(rèn)為使用的狀態(tài)方程參數(shù)是所要標(biāo)定的參數(shù)。按上述方法，獲得該JWL狀態(tài)方程參數(shù)，見表3。計算和實驗得到圓筒壁膨脹距離與時間的關(guān)系曲線如圖5所示，圖6為圓筒壁膨脹速度與膨脹距離之間的關(guān)系曲線。

表3 HMX 基PBX 爆轟產(chǎn)物的JWL狀態(tài)方程參數(shù)Table 3 The parameters of the JWL equation state of detonation product for the HMX-based PBX

圖5 圓筒壁膨脹距離與時間之間的關(guān)系曲線Fig.5 The expansion distance-time curve of cylinder

圖6 圓筒壁膨脹距離與膨脹速度之間的關(guān)系曲線Fig.6 The expansion distance-expansion velocity curve of cylinder

從圖5和圖6都可以看出，計算值和實驗值的重合性非常好，誤差小于1%。從圖6可以看出，圓筒壁在初始膨脹階段壁速上升很快，但經(jīng)過一小段膨脹后，壁速上升明顯減緩，并迅速出現(xiàn)了壁速平臺，說明該HMX 基PBX 爆轟反應(yīng)區(qū)較窄，能量釋放過程較快，所以初始階段作功能力較強(qiáng)，而后續(xù)作功能力較弱。

2.3 JWL狀態(tài)方程參數(shù)的驗證

為了考察HMX 基PBX 沿軸向加速金屬的能力，同時驗證由圓筒試驗所確定的JWL狀態(tài)方程參數(shù)的有效性，設(shè)計了炸藥柱推動平板銅飛片的試驗，并進(jìn)行了數(shù)值計算。炸藥柱和銅飛片建立二分之一的軸對稱模型，炸藥采用MAT＿HIGH＿EXPLOSIVE＿BURN 模型，JWL 狀態(tài)方程參數(shù)使用表3中的數(shù)據(jù)，銅飛片采用GRUNEISEN 狀態(tài)方程和MAT＿ELASTIC＿PLASTIC＿HYDRO 力學(xué)模型，紫銅飛片的性能力學(xué)的參數(shù)由文獻(xiàn)［11］和［13］獲得，測量值與計算值的結(jié)果如圖7所示。

圖7 飛片速度測量值和計算值的比較Fig.7 Comparison of experimental results and calculated ones of the plate flyer test

由圖7可看出，在整個時間歷程中，飛片速度的測量值與計算值能較好地符合，驗證了圓筒試驗所標(biāo)定的JWL狀態(tài)方程參數(shù)的有效性。從圖7還可以看出，飛片速度在剛開始的2μs內(nèi)上升很快，之后迅速變緩，到6μs左右就接近于最終的速度，這與圓筒試驗的結(jié)果是一致的，說明該PBX 炸藥的能量釋放過程較快，初始階段作功能力較強(qiáng)，后續(xù)作功能力較弱。

3 結(jié) 論

（1）利 用Φ25mm標(biāo)準(zhǔn)圓筒試驗對HMX基PBX（90%HMX，10%黏結(jié)劑）的作功能力進(jìn)行了評價，得到了其作功能力的特征參量以及格尼速度（2.80mm／μs）。

（2）通過圓筒試驗結(jié)合數(shù)值計算獲得該HMX基PBX 的爆轟產(chǎn)物JWL 狀態(tài)方程參數(shù)，通過平板飛片試驗驗證了所得參數(shù)正確，為該炸藥今后的應(yīng)用和理論計算提供基礎(chǔ)數(shù)據(jù)。

［1］Kury J W.Metal acceleration by chemical explosives［C］∥4th Symp（International）on Detonation.White Oak：Office of Naval Research，1965：3-12.

［2］Lee E，F(xiàn)inge M，Collins W.JWL equation of state coefficients for high explosives.UCID-16189［R］.Livermore：Lawrence Livermore National Laboratory，1973.

［3］Souers P C，Haselman L C.Detonation equation of state at LLNL 1993.UCRL-ID116113［R］.Livermore：Lawrence Livermore National Laboratory，1994.

［4］于川，劉文翰，李良忠，等.鈍感炸藥圓筒試驗與爆轟產(chǎn)物JWL 狀態(tài)方程研究［J］.高壓物理學(xué)報，1997，11（3）：227-233.

YU Chuan，LIU Wen-han，LI Liang-zhong，et al.Studies on cylinder test and JWL equation of state of detonation product for insensitive high explosive［J］.Chinese Journal of High Pressure Physics，1997，11（3）：227-233.

［5］陳朗，馮長根，黃毅民.含鋁炸藥圓筒試驗及爆轟產(chǎn)物JWL狀態(tài)方程研究［J］.火炸藥學(xué)報，2001（3）：13-15.

CHEN Lang，F(xiàn)ENG Chang-gen，HUANG Yi-min.The cylinder test and JWL equation of state detonation product of aluminized explosives［J］.Chinese Journal of Explosives and Propellants，2001（3）：13-15.

［6］盧校軍，王蓉，黃毅民，等.兩種含鋁炸藥做功能力與JWL狀態(tài)方程研究［J］.含能材料，2005，13（3）：

144-147.LU Xiao-jun，WANG Rong，HUANG Yi-min，et al.Study on work ability and JWL equation of state of two aluminized explosives［J］.Chinese Journal of Energetic Materials，2005，13（3）：144-147.

［7］韓勇，黃輝，黃毅民，等.含鋁炸藥圓筒試驗與數(shù)值模擬［J］.火炸藥學(xué)報，2009，32（4）：14-17.

HAN Yong，HUANG Hui，HUANG Yi-min，et al.Cylinder test of aluminized explosives and its numerical simulation［J］.Chinese Journal of Explosives and Propellants，2009，32（4）：14-17.

［8］計冬奎，高修柱，肖川，等.含鋁炸藥做功能力和JWL狀態(tài)方程尺寸效應(yīng)研究［J］.兵工學(xué)報，2012，33（5）：552-555.

JI Dong-kui，GAO Xiu-zhu，XIAO Chuan，et al.Study on dimension effect of accelerating ability and JWL equation of state for aluminized explosive［J］.Acta Armamentarii，2012，33（5）：552-555.

［9］胡紹樓.激光干涉測速技術(shù)［M］.北京：國防工業(yè)出版社，2001.

［10］張寶坪，張慶明，黃風(fēng)雷.爆轟物理學(xué)［M］.北京：兵器工業(yè)出版社.2001.

［11］Lee E L.Adiabatic expansion of high explosive detonation products，UCRL250422［R］.Livermore：Lawrence Livermore National Laboratory，1968.

［12］Daniel J S.An equation of state for polymethyl-methacrylate，UCID-16982［R］.Livermore：Lawrence Livermore National Laboratory，1975.

［13］Lee E，Breithaupt D，Mcmillan C，et al.The motion of thin metal walls and the equation of state of detonation products，UCID-91490［R］.Livermore：Lawrence Livermore National Laboratory，1985.

［14］Johnson G R，Cook W H，A constitutive model and data for metals subjected to large strain，high strain rates and high temperatures［C］∥Seventh International Symposium on Ballistics.Netherlands：Hague，1983：541-547.