4. 空气域与炸药建模方法

在爆炸载荷分析中，空气域与炸药是两类必须同时存在的介质：炸药提供能量释放，空气负责把爆轰产物中的能量以冲击波形式传递到结构表面。LS-DYNA 中通常采用 *MAT_NULL 配合 *EOS_LINEAR_POLYNOMIAL 描述空气，采用 *MAT_HIGH_EXPLOSIVE_BURN 配合 *EOS_JWL 描述炸药，两者都通过状态方程把压力与密度、内能联系起来。本节围绕材料模型、起爆点、网格尺寸影响和初始体积填充四个关键点展开。

一、空气材料模型

空气在爆炸问题中一般按理想气体处理，使用空材料模型加线性多项式状态方程。空材料本身不提供剪切强度，只提供密度，压力完全由状态方程给出。线性多项式状态方程形式为：

p = C0 + C1*μ + C2*μ^2 + C3*μ^3
    + (C4 + C5*μ + C6*μ^2) * E
μ = ρ/ρ0 - 1

对理想气体，取 C0=C1=C2=C3=C6=0，C4=C5=γ-1，则退化为 p=(γ-1)ρE。空气常用参数为：密度 ρ0=1.225 kg/m³，初始内能 E0=2.5×10⁵ J/m³（对应约 1 个大气压、常温），比热比 γ=1.4。初始内能必须给足，否则空气无法维持 101325 Pa 的初始压力，冲击波到达前结构表面会出现异常负压。在 *INITIAL_DETONATION 或 *SECTION_SOLID 中还需注意单位制一致，推荐使用 m-kg-s-Pa 单位制。

二、炸药材料模型与 JWL 状态方程

炸药采用高能炸药燃烧模型，需要给出密度 ρ0、爆速 D、CJ 压力 PCJ 三个参数，并用 JWL 状态方程描述爆轰产物膨胀过程：

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

其中 V 为相对体积，E 为单位体积内能。以 TNT 为例，常用参数为：ρ0=1630 kg/m³，D=6930 m/s，PCJ=21 GPa，A=3.712×10¹¹ Pa，B=3.231×10⁹ Pa，R1=4.15，R2=0.95，ω=0.30，E0=7.0×10⁹ J/m³。这些参数决定爆轰波阵面压力峰值和产物膨胀的衰减速度，直接影响近场冲量。若使用 ANFO 或乳化炸药，需替换为对应爆速与 CJ 压力，不能直接套用 TNT 参数。

三、起爆点设置

起爆点通过 *INITIAL_DETONATION 指定，关键字中给出起爆单元或节点集合以及起爆时刻。起爆点位置决定爆轰波的传播方向和到达结构各部位的时间差。对于球形装药，起爆点应设在几何中心；对于柱状装药，可设在端面中心形成轴向爆轰；对于多点起爆，可设置多个起爆点实现波形整形。起爆时刻通常设为 0，若需要延迟起爆，可给不同起爆点不同时刻。起爆点必须落在炸药单元内部，若落在空气或结构上，爆轰无法正常触发，计算会出现炸药不反应、压力不上升的现象。

四、网格尺寸对冲击波传播的影响

冲击波阵面厚度很薄，网格尺寸直接决定能否分辨波阵面。经验上，网格特征尺寸 Δx 应不大于最短波长的 1/10，对于近场爆炸，通常要求 Δx ≤ 5 mm 才能把峰值超压误差控制在 5% 以内。若 Δx 取 20 mm 以上，波阵面被数值抹平，峰值压力明显偏低，冲量也偏小，结构响应会被低估。空气域网格与炸药网格应尽量共节点或采用 tied 接触，避免出现滑移面导致能量传递损失。空气域外边界需设置无反射边界条件 *BOUNDARY_NON_REFLECTING，否则冲击波在边界反射回来会污染结果。

五、初始体积填充技术

当炸药形状为球形、柱形或椭球等规则几何时，不必按形状逐单元建模，可直接在规则六面体网格中用 *INITIAL_VOLUME_FRACTION_GEOMETRY 填充。该关键字支持 TYPE=3（球）、TYPE=4（柱）、TYPE=5（椭球）等，通过给定几何中心、半径或轴向参数，把落在几何体内的单元自动标记为炸药材料，其余保持空气。这样做的好处是空气域和炸药共用一套网格，界面自然共节点，不需要额外定义接触，也避免了复杂几何建模带来的网格质量问题。使用时注意填充几何必须完全落在空气域网格范围内，且起爆点应设在填充区域内部。

综合来看，空气与炸药建模的核心是“材料参数正确、起爆位置正确、网格足够细、填充方式合理”。这四点中任何一项出错，都会在压力时程曲线上直接反映出来：峰值偏低、到达时间偏移或出现非物理振荡。建议在正式计算前，先用一维或小尺寸二维模型验证空气初始压力和爆轰波到达时间，确认无误后再扩展到三维结构响应分析。