4. JWL状态方程与炸药材料建模
在LS-DYNA爆炸载荷分析中，炸药的爆轰过程并不是通过求解化学反应动力学得到的，而是用一套“预设的爆轰产物状态方程”来近似描述。这套状态方程就是JWL（Jones-Wilkins-Lee）方程，配合*MAT_HIGH_EXPLOSIVE_BURN材料模型，构成LS-DYNA中最常用的炸药建模方式。理解JWL参数的物理意义、关键字字段的填写规则，以及起爆点与起爆时间的控制方式，是保证爆炸载荷计算正确的前提。
一、JWL状态方程的物理意义
JWL方程描述的是爆轰产物在膨胀过程中的压力—体积—能量关系，其形式为：
P = A·(1 - ω/(R1·V))·exp(-R1·V)
+ B·(1 - ω/(R2·V))·exp(-R2·V)
+ ω·E/V其中P为爆轰产物压力，V为相对体积（当前体积与初始炸药体积之比），E为单位体积内能，A、B、R1、R2、ω为五个待定参数。前两项是高压段的主导项：A·exp(-R1·V)在V接近1时贡献极大，反映爆轰波阵面附近的高压；B·exp(-R2·V)在中等膨胀阶段起主要作用；最后一项ω·E/V是低压膨胀段的理想气体修正项，ω越大，低压段压力衰减越慢。
参数A、B的量纲是压力，R1、R2和ω是无量纲数。典型TNT的JWL参数为：A=3.712×10¹¹ Pa，B=3.231×10⁹ Pa，R1=4.15，R2=0.95，ω=0.30，初始密度ρ₀=1630 kg/m³，爆速D=6930 m/s，CJ压力P_CJ=21 GPa。这些数值并非随意给定，而是通过圆筒试验（Cylinder Test）标定得到的——圆筒试验测量炸药驱动金属圆筒的膨胀速度历史，再反推JWL参数。因此，如果使用非标准炸药，必须确认参数来源，不能随意套用TNT数据。
二、*MAT_HIGH_EXPLOSIVE_BURN关键字
在LS-DYNA关键字文件中，炸药材料通过*MAT_HIGH_EXPLOSIVE_BURN定义，其典型卡片如下：
*MAT_HIGH_EXPLOSIVE_BURN
$ MID RO D PCJ BETA K G SIGY
1 1.63 0.693 0.21 0.0 0.0 0.0 0.0
*EOS_JWL
$ EOSID A B R1 R2 OMEG E0 V0
1 3.712E11 3.231E9 4.15 0.95 0.30 7.0E9 1.0各字段含义如下：MID为材料编号；RO为炸药初始密度（单位需与模型单位制一致，上述为g/cm³-μs-Mbar制下的值）；D为爆速；PCJ为CJ爆轰压力；BETA为燃烧分数，通常取0表示未起爆，程序在起爆后自动将其置为1；K、G、SIGY在标准爆轰模型中一般置零，因为爆轰产物的强度由EOS主导。
需要特别注意的是，*MAT_HIGH_EXPLOSIVE_BURN必须与*EOS_JWL配合使用，EOSID字段将两者关联。E0为单位体积初始内能，V0为初始相对体积，通常取1.0。E0的取值直接影响爆轰能量，对于TNT，E0约为7.0×10⁹ J/m³（在g/cm³-μs-Mbar制下为0.07 Mbar）。若E0设置过小，爆轰驱动能力会明显不足。
三、炸药单元定义与网格要求
炸药单元通常采用实体单元（*SECTION_SOLID）或SPH粒子（*SECTION_SPH）。对于Lagrange网格，炸药区域需要与周围结构共节点或设置接触。由于爆轰过程中炸药单元会发生极大畸变，Lagrange网格容易出现负体积，因此实际工程中常采用以下策略之一：
- 使用
*MAT_NULL配合EOS描述爆轰产物，避免强度项引起的数值问题； - 采用ALE或SPH方法，让炸药材料在欧拉网格或粒子中流动；
- 对炸药区域使用较细网格，单元尺寸一般建议不大于最小结构壁厚的1/2，以捕捉爆轰波前压力梯度。
炸药单元的PART定义需要同时引用材料、截面和EOS。例如：
*PART
$ PID SECID MID EOSID HGID GRAV ADPOPT
10 1 1 1 0 0 0其中EOSID必须指向*EOS_JWL的编号，否则炸药不会产生爆轰压力。
四、起爆点与起爆时间控制
起爆通过*INITIAL_DETONATION关键字定义，其卡片格式为：
*INITIAL_DETONATION
$ PID X Y Z LT
10 0.0 0.0 0.0 0.0PID为炸药PART编号，X、Y、Z为起爆点坐标，LT为起爆时间。当LT=0时，计算开始即起爆；若设置LT>0，则炸药在指定时刻才被点燃。起爆点位置决定了爆轰波的传播方向和到达结构表面的时间。对于多点起爆，可以定义多个*INITIAL_DETONATION卡片，每个卡片对应一个起爆点，程序会取各起爆点最早到达的时间作为该单元的起爆时刻。
起爆时间控制在实际分析中有两个典型用途：一是模拟延时起爆，例如毫秒级延期雷管；二是用于“预加载”场景，先让结构在自重或预压下稳定，再触发爆炸。需要注意的是，LS-DYNA的起爆判据是“起爆点距离/爆速”与LT之和，即单元起爆时刻t_ign = LT + d/D，其中d为单元到起爆点的距离。因此爆速D的准确性直接影响爆轰波到达时间。
*INITIAL_DETONATION中引用的PID不存在或坐标位于炸药区域之外，程序可能不报错但炸药不会起爆，表现为压力始终为零。建议在提交计算前检查PART编号与坐标范围。五、参数设置与结果判断
以一个典型示例说明：TNT炸药尺寸为100 mm×100 mm×100 mm，起爆点位于中心，采用g/cm³-μs-Mbar单位制。设置RO=1.63，D=0.693 cm/μs，PCJ=0.21 Mbar，A=3.712 Mbar，B=0.03231 Mbar，R1=4.15，R2=0.95，ω=0.30，E0=0.07 Mbar，V0=1.0。起爆点坐标取(5.0, 5.0, 5.0) cm，LT=0。
计算开始后，可以通过d3plot观察炸药单元的压力云图。正常起爆时，起爆点附近压力迅速升至CJ压力约21 GPa，并以球面波形式向外传播。若压力峰值远低于P_CJ，应检查E0和A、B是否匹配；若压力波传播速度明显偏离D，应检查密度和EOS参数的单位一致性。对于结构响应，重点关注爆轰波到达结构表面的时刻和峰值压力，这两个量直接决定结构的变形和破坏模式。
总结来看，JWL状态方程与*MAT_HIGH_EXPLOSIVE_BURN的配合使用，核心在于三点：参数来源可靠、单位制统一、起爆定义准确。参数A、B、R1、R2、ω决定了爆轰产物压力随膨胀的变化规律，E0和V0决定了初始能量状态，而*INITIAL_DETONATION则控制爆轰波何时、从何处开始传播。只有这些要素全部正确，后续的结构响应分析才有意义。