1. 爆炸载荷与结构响应分析概述

爆炸载荷与结构响应分析，研究的是“爆炸能量如何转化为作用于结构的压力脉冲，以及结构在该脉冲下如何变形、破坏或保持完整”这一完整链条。在LS-DYNA中开展此类分析，核心难点并不在于软件操作本身，而在于对爆炸现象、冲击波传播规律和结构响应类型的准确判断——只有先明确物理过程，才能正确选择材料模型、状态方程、网格尺度与求解策略。

一、爆炸现象的分类

从力学分析角度，爆炸可按能量释放方式和约束条件分为几类，它们对载荷建模的要求差异明显。

• 凝聚炸药爆炸（化学爆炸）：TNT、RDX、乳化炸药等，能量在微秒量级释放，产生高温高压爆轰产物。分析中通常用JWL状态方程描述爆轰产物膨胀，用炸药单元或*LOAD_BLAST关键字施加压力。
• 气体爆炸（可燃气体云爆）：甲烷、氢气泄漏后遇火源爆燃或爆轰，压力峰值较低但持续时间长，常需考虑火焰加速与障碍物扰动。
• 粉尘爆炸：可燃粉尘在受限空间内被点燃，压力上升速率受粉尘浓度、粒径和湍流度影响，多用于工业安全评估。
• 物理爆炸：压力容器破裂、蒸汽爆炸等，本质是压缩能或相变能释放，载荷由容器失效压力和介质膨胀决定。
• 接触爆炸与非接触爆炸：按爆源与结构表面是否接触区分。接触爆炸产生局部冲切和层裂，非接触爆炸以冲击波整体加载为主，二者在网格与算法选择上完全不同。

分类的意义在于：它直接决定载荷是“由炸药单元自行计算”还是“由经验公式或映射方式施加”。例如近场接触爆炸通常必须建立炸药与空气网格，而远场非接触爆炸可用CONWEP经验模型简化。

二、冲击波传播特征

空气中爆炸后，爆轰产物迅速膨胀压缩周围空气，形成向外传播的冲击波。冲击波阵面处压力、密度、温度发生突跃，随后按指数规律衰减。描述其时间历程的关键参数包括：峰值超压Δp_max、正压作用时间t₊、比冲量i（压力-时间曲线下面积）以及到达时间t_a。

冲击波传播具有几个必须注意的特征：其一，超压随距离衰减，自由场中近似服从比例距离Z=R/W^1/3的相似律，W为TNT当量质量；其二，反射增强，冲击波遇到刚性壁面时正反射可使超压升高约2倍（强激波极限约8倍），斜反射则更复杂；其三，绕射与遮挡，结构后方和凹角处会出现压力降低或局部聚焦；其四，马赫反射，当入射角超过临界角时形成马赫杆，近地面爆炸中尤为常见。

在LS-DYNA中，若采用ALE或S-ALE方法，冲击波由空气与炸药材料共同求解，能自然捕捉反射、绕射和马赫反射；若采用*LOAD_BLAST或*LOAD_BLAST_ENHANCED，则通过经验曲线直接施加到结构表面，计算量小但难以反映复杂遮挡效应。

三、结构响应类型

结构在爆炸载荷下的响应，可按变形程度和失效模式划分。

• 弹性响应：载荷较小，结构应力低于屈服强度，卸载后恢复原状，通常用模态叠加或线性瞬态分析即可。
• 塑性大变形：薄板、壳体在冲击波下产生永久挠曲，需考虑应变率效应（如Cowper-Symonds模型）和几何非线性。
• 局部破坏：接触爆炸或近场爆炸造成冲孔、撕裂、层裂，需高分辨率网格和失效准则（如等效塑性应变、最大主应力）。
• 整体失稳与倒塌：框架、网壳结构在冲击波下发生构件屈曲、节点断裂乃至连续倒塌，需考虑连接非线性与接触。
• 动态响应中的惯性效应：爆炸载荷持续时间短，结构响应往往由惯性主导，准静态方法不再适用，必须进行显式动力学求解。

判断响应类型的关键，是比较载荷特征时间与结构自振周期。当t₊远小于结构周期时，冲量起控制作用；当二者接近时，需同时考虑峰值与冲量。

四、LS-DYNA在爆炸分析中的定位与典型工作流程

LS-DYNA是显式动力学有限元程序，其优势在于：支持炸药爆轰、空气冲击波、结构大变形、接触、失效与流固耦合的统一求解；提供ALE、S-ALE、SPH、DEM、CONWEP等多种爆炸载荷建模途径；具备丰富的材料本构与状态方程库。

典型工作流程可概括为：问题定义→几何与网格→材料与状态方程→载荷与边界→求解控制→结果判读。其中，炸药与空气区域常用欧拉网格，结构用拉格朗日网格，二者通过*CONSTRAINED_LAGRANGE_IN_SOLID耦合；若采用CONWEP，则只需结构网格和*LOAD_BLAST关键字。

求解控制方面，爆炸分析通常需要设置较小的时步、控制沙漏能、监控能量平衡。结果判读关注：超压时程、结构位移、塑性应变、失效单元比例以及能量曲线是否合理。

上图把本节的核心逻辑串成一条链：爆源类型决定状态方程与载荷建模方式，冲击波传播特征决定是否需要捕捉反射与绕射，结构响应类型决定网格尺度与失效准则，最终在LS-DYNA中通过ALE、CONWEP或流固耦合实现求解，并以超压、位移、塑性应变和能量曲线作为判读依据。理解这条链，是后续各章开展具体建模与调试的基础。