爆轰产物组成和爆轰参数计算方法的理论研究

杜明燃; 汪旭光; 郭子如; 颜事龙

doi:10.11883/1001-1455(2015)04-0449-05

爆轰产物组成和爆轰参数计算方法的理论研究

doi: 10.11883/1001-1455(2015)04-0449-05

1.
安徽理工大学土木建筑学院，安徽淮南 232001
2.
北京矿冶研究总院，北京 100044
3.
安徽理工大学化学工程学院，安徽淮南 232001

基金项目: 国家自然科学基金项目(51134012)

详细信息

作者简介:
杜明燃(1987-), 男, 博士研究生, dumingranaust@163.com

中图分类号: O381
计量
- 文章访问数: 4767
- HTML全文浏览量: 668
- PDF下载量: 847
- 被引次数: 0
出版历程
- 收稿日期: 2013-04-12
- 修回日期: 2014-02-27
- 刊出日期: 2015-08-10

Theoretical studies for calculating the detonation products and properties of explosives

1.
School of Civil Engineering and Architecture, Anhui University of Science and Technology, Huainan 232001, Anhui, China
2.
Beijing General Research Institute of Mining and Metallurgy, Beijing 100044, China
3.
School of Chemical Engineering, Anhui University of Science and Technology, Huainan 232001, Anhui, China

摘要

摘要: 为实现爆轰产物组成和爆轰参数的计算，采用拉格朗日乘数法和牛顿迭代的方法预测爆轰产物组成，利用BKW状态方程预测爆轰参数，在0~600 GPa和300~15 000 K压力温度范围内选取金刚石作为碳的生成相；对爆轰产物系统采用最小自由能原理，结合牛顿迭代法求解爆轰产物的化学平衡方程组；对BKW状态方程参数提出修订，取α=0.5，β=0.298，θ=6 620，κ=9.50；采用自编程序实现计算过程。使用此方法和Hugoniot关系计算密度为1.77 g/cm³的PETN爆轰CJ点爆轰参数验证计算精度，结果显示计算与实验结果的误差均小于1%。利用此方法结合Hugoniot关系预测出爆轰CJ点的产物密度为2.43 g/cm³。
- 爆炸力学 /
- 爆轰参数 /
- 热力学性质 /
- 状态方程 /
- 吉布斯自由能
Abstract: In order to calculate the detonation products and parameters, Lagrange multiplier and Newton iterative method were used to predict detonation products. The state equation of BKW was used to predict detonation parameters. In a range of pressure from 0 to 600 GPa and temperature from 300 to 15 000 K, diamond was intended as the elemental carbon product. Based on the principle of minimum free energy, the equilibrium compositions of detonation products were calculated by using Newton iterative method, which need not calculate the free energy of each composition. The parameters of the state equation of BKW were modified.α=0.5; β=0.298; θ=6 620; κ=9.50. Using self-made program, the detonation properties at CJ point of PETN, whose density is 1.77 g/cm³, were calculated with the theory in this paper and the equation of Hugoniot. The results show satisfactory agreement with the experimental data, with the error less than 1%. The density of detonation products is also predicted easily. When the density of PETN is 1.77 g/cm³, the density of detonation products is 2.43 g/cm³.
- mechanics of explosion /
- detonation parameters /
- thermodynamic property /
- equation of state /
- Gibbs free energy

HTML全文

表 1 BKW状态方程参数取值

Table 1. The parameter values of BKW equition of state

炸药	α	β	θ	κ
RDX	0.5	0.16	400	10.91
工业炸药	0.5	0.16	2 890	10.91
本文中	0.5	0.298	6 620	9.50

下载: 导出CSV

表 2 每摩尔PETN炸药CJ点爆轰产物组成

Table 2. Compositions of detonation products per mole of PETN explosive at CJ point

方法	n(CO₂)/mol	n(CO)/mol	n(N₂)/mol	n(O₂)/mol	n(H₂)/mol
BKW^{[1, 2]}	3.950 9	0.096 2	1.999 2	1.359 1×10^-5	0.568 9×10^-4
LJD^{[1, 2]}	3.779 0	0.446 9	1.969 5	0	0.480 4×10^-1
WCA^[13]	3.934 9	0.154 6	1.991 7	4.494 3×10^-5	3.989 5×10^-4
本文	3.999 5	0.001 1	1.999 9	1.190 4×10^-5	2.077 7×10^-4
方法	n(H₂O)/mol	n(NH₃)/mol	n(NO)/mol	n(CH₄)/mol	n(C)/mol
BKW^{[1, 2]}	3.998 6	1.327 5×10^-4	2.244 1×10^-4	0.537 3×10^-5	0.951 6
LJD^{[1, 2]}	3.929 1	0	0.619 5×10^-1	1.197 9×10^-2	0.759 8
WCA^[13]	3.975 2	1.629 0×10^-2	3.141 7×10^-4	0	0.910 5
本文	3.999 7	7.354 8×10^-8	4.530 7×10^-6	5.482 0×10^-6	0.999 1

下载: 导出CSV

表 3 PETN(1.77 g/cm³)炸药爆轰CJ点的爆轰参数

Table 3. Detonation parameter of CJ for PETN explosive(1.77 g/cm³)

方法	p/GPa	D/(m·s^-1)	T/K	方法	p/GPa	D/(m·s^-1)	T/K
WCA^[2]	31.5	8 416	4 349	LJD^[6]	28.5	8 087	4 378
BKW^[5]	31.8	8 421	2 833	WCA^[7]	33.49	8 540	3 237
实验^[2]	33.5	8 300	3 400	本文	33.53	8 342	3 473

下载: 导出CSV

参考文献(13)

[1]	Chirat R, Pittion-Rossillon G. A new equation of state for detonation products[J]. Journal of Chemical Physics, 1981, 74(8): 4634-4645. doi: 10.1063/1.441653
[2]	Mader C L. Numerical modeling of explosives and propellants[M]. 3rd Edtion. Boca Raton, FL, USA: CRC Press, 2008: 31-63.
[3]	Thiel M V, Ree F H. Nonequilibrium effects of slow diffusion controlled reactions on the properties of explosives[C]∥Proceedings of 9th International Symposium on Detonation. Portland, USA: Oregon, 1991: 743-750.
[4]	Wu X. BKW equation of state for detonation products[C]∥Proceedings of 8th International Symposium on Detonation. Portland, USA: Oregon, 1985: 438.
[5]	李德华, 程新路, 杨向东, 等. PETN、RDX和HMX炸药爆轰参数的数值模拟[J].爆炸与冲击, 2005, 25(4): 325-328. Li De-hua, Cheng Xin-lu, Yang Xiang-dong, et al. Numerical simulation of detonation parameters for PETN, RDX and HMX explosives[J]. Explosion and Shock Waves, 2005, 25(4): 325-328.
[6]	赵艳红, 刘海风, 张弓木. PETN炸药爆轰产物状态方程的理论研究[J].高压物理学报, 2009, 23(2): 143-149. Zhao Yan-hong, Liu Hai-feng, Zhang Gong-mu. Equation of state of detonation products for PETN explosive[J]. Journal of High Pressure Physics, 2009, 23(2): 143-149.
[7]	Fried L E, Howard W M. Explicit Gibbs free energy equation of state applied to the carbon phase diagram[J]. Physical Review, 2000, B61(13): 8734-8743.
[8]	赵艳红, 刘海风, 张弓木.基于统计物理的爆轰产物物态方程研究[J].物理学报, 2007, 56(8): 4791-4797. Zhao Yan-hong, Liu Hai-feng, Zhang Gong-mu. Equation of state of detonation products based on statistical mechanical theory[J]. Acta Physica Sinica, 2007, 56(8): 4791-4797.
[9]	杨祖一.国外民爆器材法规和技术资料汇编[M].北京: 中国爆破器材行业协会, 2010: 290-298.
[10]	杨向东, 谢文, 武保剑.液氮的冲击压缩理论计算[J].高压物理学报, 1998, 12(1): 1-5. Yang Xiang-dong, Xie Wen, Wu Bao-jian. Theoretical calculation for the Hugoniot curves of liquid nitrogen[J]. Journal of High Pressure Physics, 1998, 12(1): 1-5.
[11]	刘福生, 陈先猛, 陈攀森, 等.液态CO₂高温高密度状态方程研究[J].高压物理学报, 1998, 12(1): 28-33. Liu Fu-sheng, Chen Xian-meng, Chen Pan-sen, et al. Equation of state of liquid CO₂ at high temperatures and high densities[J]. Journal of High Pressure Physics, 1998, 12(1): 28-33.
[12]	Ihsan B, Gregor P. Thermochemical Data of Pure Substances: Volumn 2[M]. 3rd Edtion. New York: VCH, 1995.
[13]	李德华, 程新路, 杨向东, 等. PETN炸药爆轰参数的数值模拟[J].四川师范大学学报, 2005, 28(4): 448-451. Li De-hua, Cheng Xin-lu, Yang Xiang-dong, et al. Numerical simulation of detonation parameters for PETN explosive[J]. Journal of Sichuan Normal University, 2005, 28(4): 448-451.