📚 课程目录第3章爆炸载荷的物理描述与数学建模

3. 爆炸载荷的物理描述与数学建模

在LS-DYNA中做爆炸分析，第一步不是画网格，而是把“炸药变成压力”这件事说清楚。爆炸载荷的数学建模本质上要回答三个问题：炸药内部爆轰波如何传播、爆轰产物如何膨胀做功、以及能量如何传入空气形成冲击波并作用到结构上。本节围绕JWL状态方程、CJ爆轰理论、产物膨胀和空气冲击波超压/冲量计算展开，给出可直接用于关键字输入和结果判读的建模思路。

一、爆轰波与CJ理论：给爆速和压力定基准

炸药被起爆后，爆轰波以稳定速度D向前传播。CJ（Chapman–Jouguet）理论给出一个关键结论：稳定爆轰时，波后产物相对于波阵面的速度恰好等于当地声速，即马赫数等于1。由此可确定CJ面上的爆轰压力P_CJ、密度ρ_CJ和粒子速度u_CJ。工程上常用经验关系估算：

P_CJ ≈ ρ0 * D^2 / (γ + 1)
u_CJ ≈ D / (γ + 1)
ρ_CJ ≈ ρ0 * (γ + 1) / γ

其中ρ₀为炸药初始密度，γ为产物绝热指数（TNT常取约3）。以TNT为例，ρ₀=1630 kg/m³、D=6930 m/s，可估算P_CJ≈19.6 GPa，u_CJ≈1.73 km/s。这些量是JWL参数标定的物理锚点：无论用*EOS_JWL还是*MAT_HIGH_EXPLOSIVE_BURN，爆速D和CJ压力都必须与实验值一致，否则后续冲击波超压会系统性偏离。

二、JWL状态方程：描述产物膨胀的压力-体积-能量关系

JWL（Jones–Wilkins–Lee）状态方程把爆轰产物的压力写成相对体积V和单位体积内能E的函数：

P = A*(1 - ω/(R1*V))*exp(-R1*V)
  + B*(1 - ω/(R2*V))*exp(-R2*V)
  + ω*E/V

式中A、B、R₁、R₂、ω为拟合参数，V=V₀/V_s为相对体积。第一项主导高压段（爆轰初期），第二项主导中压段（产物膨胀），第三项是理想气体型的热贡献。LS-DYNA中典型TNT参数为：A=3.712×10¹¹ Pa、B=3.231×10⁹ Pa、R₁=4.15、R₂=0.95、ω=0.30，初始内能E₀=7.0×10⁹ J/m³。这些参数不是随便填的，它们共同决定了产物膨胀到不同体积时的压力衰减曲线。

JWL参数必须与爆速、CJ压力自洽。若只改D而不调A、B，会出现“爆速对但超压错”的典型问题。建议先用一维楔形模型验证CJ面压力，再用于三维结构。

三、爆炸产物膨胀与空气冲击波的形成

爆轰完成后，产物以极高压力向外膨胀，推动周围空气。空气被强烈压缩形成冲击波，波阵面压力突跃到峰值，随后按指数规律衰减。冲击波在自由场中传播时，超压峰值ΔP和正压持续时间t₊随距离R变化。工程上常用比例距离Z=R/W^1/3（R为距离，W为TNT当量质量）来统一描述。例如，在标准大气压下，自由场球形装药的超压峰值可用经验公式估算：

ΔP_max ≈ 0.084/Z + 0.27/Z^2 + 0.7/Z^3  (MPa, Z单位m/kg^(1/3))

冲量I则是超压-时间曲线下的面积，I=∫ΔP(t)dt。对结构响应而言，冲量往往比峰值更关键：当结构自振周期远小于t₊时，峰值主导；当周期远大于t₊时，冲量主导。LS-DYNA中可通过*LOAD_BLAST或*LOAD_BLAST_ENHANCED直接施加经验超压曲线，也可用*ALE+*CONSTRAINED_LAGRANGE_IN_SOLID做流固耦合，让冲击波自然生成。

四、超压与冲量的计算与LS-DYNA实现

在LS-DYNA中，空气冲击波载荷有三种常用施加方式。第一种是*LOAD_BLAST，直接输入TNT当量W、起爆点坐标和比例距离，程序按经验公式生成超压时程，适合快速评估。第二种是*LOAD_BLAST_ENHANCED，支持更复杂的反射和绕射修正。第三种是*ALE流固耦合，用*SECTION_POINT_SOURCE或*INITIAL_VOLUME_FRACTION_GEOMETRY定义炸药，空气用*MAT_NULL+*EOS_LINEAR_POLYNOMIAL，炸药区域用*MAT_HIGH_EXPLOSIVE_BURN+*EOS_JWL，通过*CONSTRAINED_LAGRANGE_IN_SOLID把压力传给结构。

以一块3 m×3 m、厚10 mm的Q345钢板为例，TNT当量10 kg，爆距5 m，比例距离Z=5/10^1/3≈2.32 m/kg^1/3。按经验公式估算ΔP_max约0.05 MPa，t₊约10 ms。若用*LOAD_BLAST，只需在关键字中写入：

*LOAD_BLAST
$#  W    R    X    Y    Z    T    ...
   10.0  5.0  0.0  0.0  0.0  0.0
$#  BLAST  ...
     1

若用ALE方法，则需注意网格尺寸：空气域网格在爆源附近应足够细，一般建议最小单元尺寸不超过装药半径的1/5，否则冲击波峰值会被数值耗散抹平。同时，*CONTROL_TIMESTEP中时间步长应满足CFL条件，避免爆轰波传播失真。

ALE方法中，炸药和空气的网格必须共节点或采用*ALE_MULTI-MATERIAL_GROUP，否则会出现材料界面穿透。流固耦合的*CONSTRAINED_LAGRANGE_IN_SOLID中，NQUAD和CTYPE的选择直接影响压力传递精度，建议先用简单平板验证。

五、结果判读：从压力曲线到结构损伤

计算完成后，应提取结构迎爆面的压力时程和位移时程。判断是否达到设计目标，通常看三个指标：峰值超压是否超过材料屈服对应的临界值、冲量是否引起足够大的塑性变形、以及是否出现撕裂或穿透。对于钢板，若等效塑性应变超过0.05～0.10，通常认为进入显著塑性阶段；若超过0.20，需警惕撕裂风险。这些阈值应结合材料本构和应变率效应（*MAT_PLASTIC_KINEMATIC或*MAT_JOHNSON_COOK）综合判断。

最后强调一点：JWL参数、CJ理论、超压经验公式三者必须相互校核。建议在正式三维计算前，先用一维或二维轴对称模型验证爆轰波到达时间和峰值压力，确认无误后再扩展到全尺寸模型。这样既能节省计算成本，也能避免因参数错误导致整个分析失效。