炸药JWL状态方程的参数拟合.pdf

节选段落一：
炸药 JWL 状态方程的参数拟合
cucral
研究爆炸力学方向的版友经常用到 LS-DYNA 中的 JWL 状态方程。但是由于其参数不
好查找，给大家的模拟带来了一定的困难。本文用比较简单的“仅考虑分子间排斥作用的凝
聚体炸药状态方程”（以下简称 k 方程）来拟合得出 JWL 方程的参数。
1. K 方程
k
Ap
v
=
1( 1) k
Ae
k v −=
−
2. 将 k
Ap
v
= 中的 v 用 JWL 中的相对体积 V 来表示。 0v v V=
则 0
0( )
k
k k
AA Mp
v V V V
ρ
= = = k 其中 0
kM Aρ= ——常数
3.


节选段落二：
根据 k
Mp
V
= V∈
1
k
k +
～3.0 的数据，利用 1stOpt（或 Matlab）拟合出
1 2
1 2
1 1R V R V Ep A e B e
RV R V V
ω ω ω− −⎛ ⎞ ⎛ ⎞
= − + − +⎜ ⎟ ⎜ ⎟
⎝ ⎠ ⎝ ⎠
中的 ，A B ， 1R ， 2R ，ω值。
（其中的 E 是与爆热成正比的常量。可由密度 1.64 的 TNT 的 E=0.07 近似求得）。
5. 例：密度 1.64，爆速 0.693，k=3 的 TNT，拟合其 JWL 参数。
1） 计算出 Gpa 比较符合实验值。


节选段落三：
p4*x)+0.07*p5/x;
data;
//x y
0.75 0.1965
0.85 0.135
0.95 0.0967
…………
…………
2.71 0.00417375
2.73 0.00408269
4） 采用差分进化法拟合结果
p1 0.0100000000000042 A
p2 B 2.93869361309839
p3 1R 7.9999999999995
p4 2R 3.69601343652259
p5 ω 0.246325358893465
5） 结果分析
Hp 的相对误差为