JB9014炸藥超壓爆轟產物的狀態方程*

陳 軍，曾代朋，孫承緯，張震宇，譚多望

（1.中國工程物理研究院流體物理研究所，四川 綿陽 621900；

2.國防科學技術大學理學院技術物理研究所，湖南 長沙 410073）

JB9014炸藥超壓爆轟產物的狀態方程*

陳 軍1，曾代朋1，孫承緯1，張震宇2，譚多望1

（1.中國工程物理研究院流體物理研究所，四川 綿陽 621900；

2.國防科學技術大學理學院技術物理研究所，湖南 長沙 410073）

根據P.K.Tang等提出的對炸藥爆轟產物超壓狀態方程建模時只對JWL狀態方程CJ等熵線中高壓指數項做修正的研究思路，首先給定超壓狀態下內能等熵線的修正項，再根據熱力學定律對內能等熵線求微分而得到沿壓力等熵線的修正項。對JB-9014炸藥超壓爆轟沖擊Hugoniot實驗數據和聲速實驗數據同時進行擬合，得到了3個JWL狀態方程在超壓爆轟狀態下的修正項，并進行了分析與比較。得到超壓修正項的方法簡單，3組超壓修正項與P.K.Tang的修正項一樣，都能很好地擬合超壓Hugoniot數據。在實驗數據范圍外，對超壓狀態下的聲速－壓力實驗數據的擬合精度有所差別。

爆炸力學；狀態方程；超壓爆轟；炸藥

1 引 言

在高速飛片撞擊炸藥、球形會聚爆轟等情況下，炸藥中會產生高于CJ爆轟的超壓爆轟狀態。由于標準的JWL方程不能很好地描述超壓爆轟狀態下的產物行為，很多學者都在尋找能夠描述爆轟產物在超壓狀態下特征的狀態方程［1－3］。P.K.Tang［3］建立了PBX-9501和 PBX-9502炸藥超壓爆轟產物狀態方程模型。他們假設爆轟產物的標準JWL狀態方程在低壓部分是可以接受的，而對超壓爆轟區中的狀態進行描述時有所欠缺。為此他們建議增加Hugoniot曲線在CJ點以上的斜率，以更好地匹配超壓爆轟Hugoniot實驗數據，同時也希望增加超壓爆轟產物的聲速。建模的思路類似于可壓縮理想氣體在高壓區的概念，只對JWL狀態方程中的高壓指數項做修正，首先給定超壓狀態下的壓力等熵線的修正項，而沿內能等熵線的修正項則根據熱力學定律對壓力等熵線進行積分求得。這樣做雖然使沿壓力等熵線和沿內能等熵線的修正項具有相同的修正系數，但是不能取較復雜的沿壓力等熵線的修正項，因為很難根據熱力學定律對壓力等熵線積分求得沿內能等熵線的修正項。

本文中，根據P.K.Tang對炸藥超壓爆轟產物狀態方程研究的建模思路，首先給定3個超壓狀態下內能等熵線高壓項的修正項，而沿壓力等熵線的修正項可以根據熱力學定律對內能等熵線微分求得。通過對JB-9014炸藥超壓爆轟沖擊Hugoniot實驗數據和聲速實驗數據同時進行擬合，得到3個可以描述JB-9014炸藥超壓爆轟狀態的JWL狀態方程，并對這3個JWL狀態方程在超壓爆轟狀態下的修正項進行比較。

2 超壓狀態下對標準JWL爆轟產物狀態方程的修正方法

根據Grüneisen狀態方程，及其在Hugoniot曲線上的應用，有

將式（3）代入式（2），解得炸藥爆轟產物的Hugoniot曲線

式中：帶下標S的量為等熵線上的變量，帶下標H的量為Hugoniot曲線上的變量，帶下標CJ的量為CJ狀態參量；V是相對比容。標準JWL狀態方程的CJ等熵線方程為

假設標準的JWL爆轟產物狀態方程可以很好地描述爆轟產物在CJ狀態以下的等熵膨脹過程，而對超壓爆轟區中的狀態進行描述時有所欠缺。因此，需要對標準JWL狀態方程的高壓指數項做修正。首先，對標準JWL狀態方程的CJ內能等熵線方程高壓指數項做如下修正

式中：Fe（V）為沿內能等熵線的修正項。

其中沿壓力等熵線的修正項Fp（V）可由內能等熵線的修正項Fe（V）求得

沿超壓等熵線可得修正后的聲速表達式為

Hugoniot聲速的表達式為

方程（9）、（10）和（12）對CJ等熵線方程和聲速的修正只在CJ狀態點以上成立，在CJ點和CJ點以下仍用標準的CJ等熵線方程和聲速方程。因此，為保證在CJ點上的連續性，沿內能等熵線修正項Fe（V）表達式的基本選取原則是，在CJ狀態點上，Fe（V）及其1、2階導數為0。

3 JB-9014炸藥超壓爆轟產物狀態方程的擬合

為了對比，給出如下3種沿內能等熵線修正項Fe（V）的解析式

分別對以上3種Fe（V）進行求導，由式（10）得出與之相匹配的Fp（V）的表達式分別為

在標準JWL狀態方程中的所有修正項，只當V＜VCJ時，即在超壓狀態下才起作用。

用于擬合其中修正項系數A0、B0和C0的JB-9014炸藥超壓Hugoniot實驗數據和超壓聲速實驗數據分別取自文獻［4-5］。由于文獻中所做的超壓實驗是利用化爆加載，壓力可調范圍受到限制，實驗數據不多。考慮到JB-9014炸藥與PBX-9502炸藥的主要成分、比例和典型壓制密度基本相同，為提高擬合精度、增加實驗數據點，擬合同時使用了文獻［3］中關于PBX-9502炸藥相應的實驗數據（通過讀取像素點獲得）。圖1給出了用于擬合超壓修正項的JB-9014炸藥超壓爆轟Hugoniot曲線和超壓Hugoniot聲速-Hugoniot壓力關系曲線的實驗數據點。

圖1 JB-9014炸藥超壓爆轟Hugoniot實驗數據Fig.1The Hugoniot data for the overdriven-detonation products of JB-9014explosive

圖2 JB-9014炸藥超壓爆轟產物的模型結果比較Fig.2 Comparison of the model results for the overdriven-detonation products of JB-9014explosive

因為JB-9014炸藥的參數不完全，在以下的計算中，JB-9014炸藥爆轟產物標準JWL狀態方程參數采用PBX-9502炸藥的相應參數［6］，其中E0＝6.9GPa，wg＝0.5，Ag＝1 361.77GPa，Bg＝71.99GPa，r1＝6.2，r2＝2.2，Cvg＝1.0MPa／K。

對修正的等熵線壓力和內能表達式中的修正項系數A0、B0和C0，利用非線性優化方法，同時擬合CJ點之上的超壓Hugoniot實驗數據和超壓聲速實驗數據得到。目標函數取圖1（a）所示實驗點與式（5）在相同相對比容處計算結果之差的最小平方和，加上圖1（b）所示實驗點與式（12）在相同Hugoniot壓力處計算結果之差的最小平方和。擬合所得3組模型最優化修正項系數的數值分別為

式中：給出了目標函數f的最小值，它代表了3組修正模型對實驗數據擬合的精度。

3組模型的擬合超壓Hugoniot曲線和聲速－壓力曲線如圖2所示，圖中實線為標準JWL爆轟產物狀態方程的計算結果。

4 結 論

通過對JB-9014炸藥標準JWL狀態方程中的高壓指數項做修正，建立了3組可以描述JB-9014炸藥超壓爆轟狀態和正常爆轟產物等熵膨脹狀態的JWL狀態方程，得出以下結論：

（1）首先給出沿內能等熵線超壓修正項Fe（V）的表達形式，然后根據熱力學第1定律得到沿壓力等熵線的超壓修正項Fp（V）的具體表達式。這種方法不僅使沿壓力等熵線和沿內能等熵線的超壓修正項有相同的修正系數，而且從Fe（V）的表達形式推導Fp（V）表達式更容易。

（2）對超壓爆轟Hugoniot曲線與超壓爆轟狀態下聲速－壓力曲線同時進行非線性優化擬合，確定修正項系數。此方法不僅精度高，而且方便、快捷。

（3）給出的3組超壓修正模型與P.K.Tang的修正模型一樣，都能很好地擬合超壓Hugoniot數據，而對超壓狀態下的聲速－壓力實驗數據的擬合，在實驗數據范圍之外精度有所差別。模型1的超壓聲速－壓力擬合曲線不通過CJ點，這是因為本文中給出的所有內能等熵線修正項Fe（V）都是（VCJ－V）的2次以上冪函數，修正模型只當相對比容V小于CJ點相對比容VCJ時，即在超壓狀態下才起作用。按照公式（10），模型1對壓力等熵線的修正項Fp（V）的最低冪次是1次，由公式（8）～（9），標準JWL壓力等熵線和內能等熵線與模型1的超壓壓力等熵線和內能等熵線在CJ點上的值相等，但由公式（11）確定的標準JWL聲速與超壓聲速在CJ點的值不相等。

［1］van Thiel M，Lee E L，Cochran S.Effects of overdriven shock states on the equation of state of 9404explosive［J］.Journal of Applied Physics，1983，54（11）：6760-6763.

［2］鄧全農.過加載下固體炸藥的實驗研究和爆轟產物的高壓狀態方程［D］.綿陽：中國工程物理研究院研究生部，1989.

［3］Tang P K.A study of overdriven behaviors of PBX-9501and PBX-9502［C］∥Proceedings of the 11th International Detonation Symposium.Colorado，USA：Office of the Naval Research，1998：1058-1064.

［4］陳軍.超壓爆轟產物狀態及相關現象研究［D］.綿陽：中國工程物理研究院研究生部，2008.

［5］曾代朋，陳軍，譚多望，等.JB-9014超壓爆轟雨貢紐實驗測量［C］∥第5屆全國計算爆炸力學會議.青海西寧，2008.

［6］Tarver C M，McGuire E M.Reactive flow modeling of the interaction of TATB detonation waves with inert materials［C］∥Proceedings of the 12th International Detonation Symposium.San Diego，California，2002.

Equations of state for overdriven-detonation products of JB-9014explosive＊

CHEN Jun1，ZENG Dai-peng1，SUN Cheng-wei1，ZHANG Zhen-yu2，TAN Duo-wang1
（1.Institute of Fluid Physics，China Academy of Engineering Physics，Mianyang 621900，Sichuan，China；
2.Institute of Technical Physics，College of Science，National University of Defense Technology，Changsha 410073，Hunan，China）

Based on P K Tang’s opinion that only the high-pressure term of Jones-Wilkins-Lee equation of state（JWL EOS）needed to be modified when used for overdriven detonation，the modified term of the internal energy isentropic curve was presented and the modified term of the pressure isentropic curve was obtained from the differential form of the internal energy isentropic curve by applying the first law of thermodynamics.And three modified JWL equations of state for overdriven detonation products of JB-9014explosive were gained by fitting simultaneously the shock Hugoniot pressure and sound velocity data.Results show that，the modified terms gained in this paper can model Hugoniot data for overdriven detonation reasonably as well as P K Tang’s，and the precision is a little worse out of the range of experimental data.

mechanics of explosion；equation of state；overdriven detonation；explosive

19November 2009；Revised 29March 2010

CHEN Jun，chengjun＠caep.ac.cn

（責任編輯 張凌云）

O381 國標學科代碼：130·3510

A

1001-1455（2010）06-0583-05

2009-11-19；

2010-03-29

國防科工局基礎科研項目（B1520110002）

陳 軍（1973— ），男，博士，副研究員。

Supported by the Bureau of Science，Technology and Industry for National Defence of China（B1520110002）