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

为实现爆轰产物组成和爆轰参数的计算,采用拉格朗日乘数法和牛顿迭代的方法预测爆轰产物组成,利用BKW状态方程预测爆轰参数,在0~600 GPa和300~15 000 K压力温度范围内选取金刚石作为碳的生成相;对爆轰产物系统采用最小自由能原理,结合牛顿迭代法求解爆轰产物的化学平衡方程组;对BKW状态方程参数提出修订,取α=0.5,β=0.298,θ=6 620,κ=9.50;采用自编程序实现计算过程。使用此方法和Hugoniot关系计算密度为1.77 g/cm3的PETN爆轰CJ点爆轰参数验证计算精度,结果显示计算与实验结果的误差均小于1%。利用此方法结合Hugoniot关系预测出爆轰CJ点的产物密度为2.43 g/cm3。