爆風荷重応答解析
理论与物理
什么是爆炸载荷
老师，什么是爆炸载荷？
是爆炸（化学炸药、气体爆炸等）产生的冲击波（爆炸波）作用在结构上的压力载荷。具有峰值压力→指数衰减→负压这种特征性的时间波形。
Friedlander波形
爆炸压力的标准近似式（修正Friedlander式）：
- $P_s$ — 峰值超压
- $t_d$ — 正压持续时间
- $b$ — 衰减参数
- $p_0$ — 大气压
是指数衰减并转为负压呢。
正压阶段（推力）之后有负压阶段（拉力）。对于建筑物的墙壁，是正压（向外）→负压（向内）的重复作用。
爆炸参数的估算
使用ConWep（常规武器效应程序）的Kingery-Bulmash式估算爆炸参数：
$Z$ 是比例距离（m/kg^{1/3}），$R$ 是距爆源的距离，$W$ 是TNT当量。由 $Z$ 估算 $P_s, t_d, I$。
仅凭炸药量和距离就能决定参数吗？
这是“比例定律”（Hopkinson-Cranz定律）。相同的$Z$就有相同的压力波形。TNT 1 kg在10 m处与TNT 1000 kg在100 m处具有相同的$Z$和相同的超压。
FEM中的设置
两种方法：
1. 直接输入压力时程 — 将Friedlander式的压力作为时间函数施加在结构面上。最简单
2. ALE法（任意拉格朗日-欧拉法） — 用欧拉网格求解爆炸波的传播，同时用拉格朗日网格求解结构的变形。包含反射、绕射
ALE法更精确吗？
ALE法能自动计算爆炸波的反射、绕射和叠加，所以更精确。但需要空气的3D网格，计算成本高。对于简单形状，直接输入压力就足够了。
总结
要点：
- Friedlander波形 — 峰值超压→指数衰减→负压
- 使用ConWep（Kingery-Bulmash式）估算参数 — $Z = R/W^{1/3}$
- 压力直接输入 或 ALE法 — 根据复杂程度选择
- LS-DYNA的*LOAD_BLAST_ENHANCED — 自动计算ConWep
爆炸压力随距离的三次方衰减
根据TNT当量比例定律（Hopkinson-Cranz定律），爆炸的最大超压由距离R与装药质量W的1/3次方的比值Z＝R/W^(1/3)来整理，当Z＝1时，最大超压可达约0.3MPa。1944年的Brode式后来被精炼为Kingery-Bulmash数据库（美军TM5-855-1），至今仍被ANSYS Autodyn的ATBLAST功能采用。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，那是“因为缓慢施加力所以加速度可以忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长性”就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高＝强度高”是不对的。刚度是“不易变形的程度”，强度是“不易破坏的程度”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内容物上的力”（体积力），轮胎压路面的力是“只作用在表面的力”（表面力）。风压、水压、螺栓紧固力…全都是外力。这里容易犯的错误：弄错载荷的方向。本想“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时，实际上确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力或弦的内部摩擦变成了热。汽车的减震器也是同样的原理——故意吸收振动能量来提高乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会那样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷、弹性模量等也要统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢为= 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
LS-DYNA中的爆炸分析
LS-DYNA有专用的爆炸功能：
```
*LOAD_BLAST_ENHANCED
1, 100., 10., 0., 0., 0., 1, 1.0 $ W=100kg TNT, R=10m
```
自动计算ConWep的爆炸压力并作用在结构的各个面上。反射角度也自动考虑。
内置了ConWep真方便啊。
*LOAD_BLAST_ENHANCED根据炸药量、爆源位置、结构面的朝向自动计算反射压力。比ALE法轻量得多，对许多问题精度足够。
ALE法的设置
1. 在结构周围布置空气的欧拉网格
2. 用*INITIAL_DETONATION等定义爆源
3. 用拉格朗日网格定义结构
4. 通过ALE耦合（*CONSTRAINED_LAGRANGE_IN_SOLID）连接流体与结构
总结
ALE法求解爆炸-结构耦合
在爆炸分析中，ALE（任意拉格朗日-欧拉）法是主流，用欧拉网格处理空气、炸药，用拉格朗日单元处理结构物。LS-DYNA的ALE多物质代码于1990年代在LLNL（劳伦斯利弗莫尔国家实验室）开发，用当前的HPC模拟TNT 1kg的自由场爆炸（100万节点，5ms时长）大约需要30分钟。
线性单元（1次单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切锁定（可通过减缩积分或B-bar法缓解）。
2次单元（带中间节点）
可以表现曲线变形。应力精度大幅提高，但自由度约增加2~3倍。推荐：应力评价很重要时使用。
完全积分 vs 减缩积分
完全积分：有过约束（锁定）风险。减缩积分：有沙漏模式（零能量模式）风险。根据情况选择。
自适应网格
基于误差指标（ZZ估计量等）的自动细化。有效提高应力集中部位的精度。有h法（单元细分）和p法（增加阶次）。
牛顿-拉弗森法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿-拉弗森法
切线刚度矩阵使用初始值或每隔几次迭代更新。每次迭代成本低，但收敛速度为线性。
收敛判定标准
力残差范数: $||R|| / ||F_{ext}|| < \epsilon$（通常 $\epsilon = 10^{-3}$〜$10^{-6}$）。位移增量范数: $||\Delta u|| / ||u|| < \epsilon$。能量范数: $\Delta u \cdot R < \epsilon$
载荷增量法
不一次性施加全部载荷，而是分小步增加。弧长法（Riks法）可以超越载荷-位移关系的极值点进行追踪。
直接法 vs 迭代法的比喻
直接法是“用笔算精确解联立方程”的方法——可靠但大规模问题太耗时。迭代法是“反复猜测逼近正确答案”的方法——最初答案粗糙，但每次迭代精度提高。就像查字典时，从第一页开始顺序查找（直接法）不如先估计位置翻开，再前后调整（迭代法）来得高效，原理相同。
网格阶次与精度的关系
1次单元是“用直尺近似曲线”——用直线折线表现，精度有限。2次单元是“柔性曲线”——可以表现曲线变化，即使网格密度相同，精度也显著提高。但是，每个单元的计算成本增加，需要根据总体的成本效益来判断。
实践指南
爆炸分析实务
用于国防、石油化工厂、反恐结构设计。
应用实例
| 应用 | 目的 |
|---|---|
| 防爆墙设计 | 屏蔽爆炸冲击波。将变形量控制在允许值内 |
| 建筑物抗爆设计 | 玻璃破坏、结构响应 |
| 车辆抗爆性能 | 对IED（简易爆炸装置）的耐受性 |
| 工厂安全距离 | 爆炸事故对结构的影响 |
实务检查清单
UFC 3-340-02是什么？
是美国国防部的针对爆炸载荷的结构设计手册。规定了峰值压力、持续时间、结构的允许变形。是防爆设计的世界标准。
防护设计的标准是UFC 4-010-01
美国国防部UFC（统一设施标准）4-010-01是政府设施的爆炸防护设计标准，规定了安全距离和结构响应限值。2001年9.11事件后修订的2012年版要求将RC板的延性比（μ）控制在10以下，并明确规定验证时使用动态分析（SDOF或FEM）。
分析流程的比喻
分析流程其实和烹饪非常相似。首先采购食材（准备CAD模型），进行预处理（相关主题
なった
詳しく
報告