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

爆炸是能量在极短时间内以冲击波、高速破片、热辐射等形式向外释放的物理过程。在工程分析中，我们关心的不是爆炸本身的化学细节，而是它施加到结构上的载荷，以及结构在这种载荷下如何变形、损伤乃至失效。LS-DYNA 之所以成为爆炸仿真领域的主流工具，正是因为它把"载荷生成"与"结构响应"放在同一个显式动力学框架内求解。本章先建立整体认知：爆炸现象怎么分类、爆炸载荷有哪些特征、结构响应分析要回答什么问题，以及 LS-DYNA 在其中处于什么位置、典型工作流程如何展开。

一、爆炸现象的分类

从载荷来源看，工程中常见的爆炸现象大致分为四类。第一类是凝聚炸药爆炸，如 TNT、RDX、HMX 及各类军用/工业炸药，特点是能量密度高、爆速可达 6000～9000 m/s，爆轰波后压力可达数十 GPa，是爆炸仿真中最典型的载荷源。第二类是气体爆炸，如甲烷、氢气、粉尘云在受限空间内的爆燃，超压通常只有几十 kPa 到几百 kPa，但持续时间长、空间约束强，常见于石化与煤矿安全评估。第三类是物理爆炸，如压力容器爆破、蒸汽爆炸，本质是高压介质突然释放，载荷由容器内能决定。第四类是高速撞击/侵彻引发的类爆炸载荷，如动能弹撞击、碎片云冲击，虽然不涉及化学反应，但产生的应力波与爆炸载荷在数学形式上高度相似。

分类的意义在于：不同爆炸现象对应的载荷模型、状态方程和网格尺度完全不同。用 TNT 的 JWL 状态方程去模拟甲烷爆燃，或用理想气体模型去模拟凝聚炸药，都会得到错误的量级。

二、爆炸载荷的特征

爆炸载荷区别于地震、风载、静载的核心特征有三点。

第一，强非线性与瞬时性。以自由场 TNT 爆炸为例，冲击波到达某点时压力在微秒量级内从零跃升到峰值，随后按指数规律衰减。典型的 Friedlander 波形可写为 p(t)=p_so(1−t/t₊)e^−bt/t₊，其中 p_so 为入射超压峰值，t₊ 为正压持续时间。峰值可达 MPa 级，而作用时间只有毫秒级，这种"高幅值、短持时"决定了必须用显式时间积分求解。

第二，传播与反射的耦合。冲击波遇到结构表面会发生反射，反射超压可达入射超压的 2～8 倍（正入射时理论上可达 8 倍，实际受马赫反射影响）。因此结构迎爆面的实际载荷往往远大于自由场入射压力，必须考虑入射角、地面反射和结构遮挡。

第三，载荷与结构变形的相互作用。当结构在载荷作用下发生大变形或破坏时，会反过来改变流场，形成流固耦合。近场爆炸、密闭空间内爆炸尤其明显，此时不能简单地把爆炸载荷当作预设压力曲线施加到结构上。

三、结构响应分析的目标

爆炸载荷下的结构响应分析，通常要回答四类问题：一是变形与位移，结构最大挠度是否超过允许值；二是应力与内力，关键截面是否进入塑性、是否发生剪切破坏；三是损伤与失效，是否出现撕裂、层裂、连接失效；四是防护效果评估，如防爆墙、夹芯板、吸能层的削波与耗能能力。

与静力分析不同，爆炸响应必须考虑应变率效应。钢材在 10³ s⁻¹ 量级的应变率下屈服强度可提高 30%～80%，LS-DYNA 中常用 Cowper-Symonds 模型描述：σ_d/σ_s=1+(ε̇/D)^1/q，其中 D、q 为材料常数。忽略应变率会显著高估结构变形。

四、LS-DYNA 的定位与典型工作流程

LS-DYNA 是显式动力学有限元程序，采用中心差分时间积分，天然适合处理毫秒级、强非线性、大变形和接触问题。在爆炸仿真中，它同时提供三类能力：用 JWL 状态方程描述炸药爆轰产物、用 ALE 或 S-ALE 方法描述爆炸产物与空气的流动、用 Lagrange 网格描述结构变形，并通过流固耦合算法把两者连接起来。

典型工作流程可以概括为：建立几何与网格 → 定义材料与状态方程 → 设置爆炸源与流场 → 定义流固耦合与接触 → 施加边界与初始条件 → 设置求解控制与输出 → 提交计算 → 后处理判读。下面这张图把这一流程与各环节的关键输入、输出对应起来，便于建立整体框架。

需要强调的是，流程中的第 8 步不是终点。爆炸仿真中网格尺寸、人工体积黏性、沙漏控制都会显著影响结果。工程上通常要求沙漏能不超过总内能的 5%，流场网格在爆源附近应足够细以分辨冲击波前沿，结构网格在预期塑性区应至少保证 3～5 层单元。若这些指标不满足，应回到第 1 步重新划分网格，而不是直接采信结果。

五、一个可参照的入门算例框架

为便于后续章节展开，这里给出一个钢板在近场爆炸下的分析框架（参数为教学示例，非实测值）：