炸藥格尼系數(shù)的一種簡易估算法

沈 飛，王 輝，袁建飛，田清政

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

引 言

格尼系數(shù)（又稱格尼速度）可以反映炸藥對金屬的加速作功能力，是炸藥的一個(gè)重要參數(shù)。工程上應(yīng)用格尼系數(shù)可以求出炸藥對殼體破片加速的最大速度，同時(shí)也作為戰(zhàn)斗部裝藥選評的一個(gè)重要參考［1－3］。目前，炸藥的格尼系數(shù)主要通過圓筒試驗(yàn)進(jìn)行確定，此外，也有一些簡易的估算方法，較為典型的是Har desty和Kennedy提出的H＆K算法及Kamlet和Finger提出的 K＆F算法［1，4］。這兩種估算方法對于CHNO型炸藥的計(jì)算誤差都較小，但前提是必須知道爆轟產(chǎn)物的成分及摩爾數(shù)，而工程上常用的很多高聚物黏結(jié)炸藥、混合炸藥等，其配方成分復(fù)雜，很難準(zhǔn)確判斷出爆轟產(chǎn)物的組成情況，因此無法通過H＆K算法及K＆F算法較為準(zhǔn)確地估算出炸藥的格尼系數(shù)。

在工程實(shí)際中，炸藥的初始密度及爆速比較容易準(zhǔn)確獲得，本研究基于圓筒裝藥結(jié)構(gòu)的格尼模型及炸藥爆轟產(chǎn)物的γ律狀態(tài)方程得出了一種根據(jù)炸藥密度及爆速估算其格尼系數(shù)的方法，并采用該方法對多種炸藥的格尼系數(shù)進(jìn)行估算，以檢驗(yàn)其準(zhǔn)確性。

1 圓筒裝藥結(jié)構(gòu)的格尼模型

在圓筒裝藥結(jié)構(gòu)格尼公式的推導(dǎo)過程中，首先假定圓筒為無限長，忽略其作用過程中的端部效應(yīng)，且圓筒膨脹過程中，整個(gè)圓筒壁具有相同的速度［5］。那么，單位長度的圓筒壁所具有的動(dòng)能為：

式中：Re和Ri分別為圓筒壁外表面和內(nèi)表面的半徑；ρm為圓筒材料的密度；um為圓筒壁的徑向膨脹速度；M為單位長度圓筒壁的質(zhì)量。

對于爆轟產(chǎn)物，假定其具有相同的密度，且徑向速度沿半徑呈線性分布［5］，即：

則單位長度的炸藥產(chǎn)生的爆轟產(chǎn)物所具有的動(dòng)能為：

式中：ρp為爆轟產(chǎn)物的密度；C為單位長度炸藥的質(zhì)量。

根據(jù)公式（1）和（3），可得出單位質(zhì)量炸藥釋放能量轉(zhuǎn)化成的動(dòng)能為：

當(dāng)圓筒壁的膨脹速度幾乎不再增加時(shí)［6－7］，可將此時(shí)的筒壁速度記為uL，且此時(shí)的Ed便等于炸藥的格尼能Eg，進(jìn)而可求出該炸藥的格尼系數(shù)：

2 簡易估算法

2.1 估算法的理論模型

在圓筒膨脹過程中，根據(jù)能量守恒關(guān)系可認(rèn)為，隨著圓筒壁的膨脹，爆轟產(chǎn)物的內(nèi)能不斷轉(zhuǎn)換成圓筒壁和爆轟產(chǎn)物的動(dòng)能，即：

式中：E0為圓筒膨脹初始時(shí)刻爆轟產(chǎn)物的內(nèi)能；Es為膨脹過程中爆轟產(chǎn)物的內(nèi)能。對于凝聚炸藥的爆轟產(chǎn)物，可采用 Grüneisen物態(tài)方程［8－9］描述，即：

式中：冷壓Pk（ρ）取為AV－γ，熱壓f（V）T＝BT／V，代入式（7）得到：

式中：T為溫度；P和V分別為爆轟產(chǎn)物的壓力和比容；γ為爆轟產(chǎn)物的多方指數(shù)；A和B均為與炸藥性質(zhì)有關(guān)的常數(shù)。對于實(shí)際中常用的凝聚炸藥，其密度一般大于1.0 g／c m3，爆轟產(chǎn)物中熱壓對壓力的作用相對于冷壓要小得多，特別是在爆轟產(chǎn)物膨脹的初始階段，熱壓的作用可以忽略［9］，則式（8）可簡化為：

式（9）即為爆轟產(chǎn)物的γ律狀態(tài)方程。

此外，C－J面上的參數(shù)有如下關(guān)系［8］：

式中：D為炸藥的爆速；ρ0和V0分別為炸藥的初始密度及比容，ρ0＝1／V0。

根據(jù)公式（9）、（10）和（11），可得出：

將公式（1）、（3）、（13）和（14）代入式（6），則可得出：

當(dāng)爆轟產(chǎn)物膨脹到一定程度時(shí)，圓筒壁的速度接近最大值，而此時(shí)爆轟產(chǎn)物的相對比容V一般都在7.0以上［6－7］，則≈1。那么，由式（5）和（15）可得出：

式（16）中多方指數(shù)γ是爆轟產(chǎn)物的體積和溫度的函數(shù)，與炸藥的成分和密度相關(guān)。有關(guān)多方指數(shù)γ的計(jì)算方法很多，但大部分需要已知爆轟產(chǎn)物的成分及摩爾分?jǐn)?shù)，這對于工程上經(jīng)常用到的一些混合炸藥而言，其計(jì)算不方便，可以采用僅需要知道炸藥密度便可求出γ的經(jīng)驗(yàn)公式。本研究對于γ值的計(jì)算主要根據(jù)Johansson和Persson提出的公式（17），其形式簡單，對于密度大于1.0 g／c m3的凝聚炸藥具有較好的精度［10］。

2.2 算 例

選擇部分已知密度和爆速的炸藥，采用公式（16）和（17）計(jì)算其格尼系數(shù)，并與實(shí)驗(yàn)值進(jìn)行對比，結(jié)果如表1和表2所示。其中，表1中所列的格尼系數(shù)實(shí)驗(yàn)值均來自文獻(xiàn)［11］，表2中所列格尼系數(shù)實(shí)驗(yàn)值是將文獻(xiàn)［12］中Φ25 mm圓筒試驗(yàn)值代入式（5）得到的。計(jì)算過程中，取圓筒壁膨脹距離為19 mm時(shí)（此時(shí)爆轟產(chǎn)物相對比容約為7）圓筒壁的膨脹速度作為uL。從表1和表2可以看出，格尼系數(shù)的計(jì)算值與實(shí)驗(yàn)值相差較小。

表1 炸藥格尼系數(shù)的估算值與文獻(xiàn)值的對比Table 1 Co mparison of calculated Gur ney coefficient with literature one

表2 格尼系數(shù)的估算值與Φ25 mm圓筒試驗(yàn)值的對比Table 2 Co mparison of calculated Gur ney coefficient with values measured fr o mΦ25 mm cylinder test

3 結(jié) 論

（1）基于圓筒裝藥結(jié)構(gòu)的格尼模型及凝聚炸藥爆轟產(chǎn)物的γ律狀態(tài)方程推導(dǎo)出了炸藥的格尼系數(shù)關(guān)于

（2）該方法對于密度大于1.0 g／c m3的凝聚炸藥具有較好的精度，能夠滿足工程應(yīng)用的需要。

［1］ 孫業(yè)斌，惠君明，曹欣茂.軍用混合炸藥［M］.北京：兵器工業(yè)出版社，1995.

［2］ 韓秀鳳，嚴(yán)楠，蔡瑞嬌.對炸藥驅(qū)動(dòng)飛片速度的理論計(jì)算方法的分析與評價(jià)［J］.火炸藥學(xué)報(bào)，2005，28（1）：63－66.HAN Xiu－feng，YAN Nan，CAI Rui－jiao.Analysis and assess ment on theoretical calculation methods for the velocity of explosive－driven flyer plate［J］.Chinese Journal of Explosives and Propellants，2005，28（1）：63－66.

［3］ Wang M F，Lu F Y，Li X Y，et al.A for mula for calculating the velocities of frag ments from velocity enhanced war head［J］.Propellants，Explosives，Pyr otechnics，2013，38：232－237.

［4］ Keshavarz M H.An alter native method for esti mation Gurney velocity based on assumed detonation products［J］.Chinese Jour nal of Ener getic Materials，2006，14（6）：499－452.

［5］ Walters W P，Zukas J A.Fundamentals of Shaped Charges［M］.Canada：John Wiley ＆Sons，1989.

［6］ Hornberg H，Volk F.The cylinder test in the context of physical detonation measurement methods［J］.Propellants，Explosives，Pyrotechnics，1989，14：199－211.

［7］ Souers P C，Lauderbach L，Garza R，et al.Upgraded analytical model of the cylinder test［J］.Pr opellants，Explosives，Pyrotechnics，2013，38：419－424.

［8］ 奧爾連科 ЛП.爆炸物理學(xué)（上冊）［M］.孫承緯，譯.北京：科學(xué)出版社，2011.

［9］ 孫承緯，衛(wèi)玉章，周之奎.應(yīng)用爆轟物理［M］.北京：國防工業(yè)出版社，2000.

［10］Johansson C H，Persson P.Density and pressure in the Chap man Jouguet Plane as f unctions of initial density of explosive［J］.Nature，1966，212：1230－1231.

［11］Dobratz B M，Crawford P C.LLNL Explosives Handbook，Properties of chemical explosives and explosives si mulants，Report UCRL－52997［R］.Lawrence：Liver more National Laboratory，1985.

［12］董海山，周芬芬.高能炸藥及相關(guān)物性能［M］.北京：科學(xué)出版社，1989.