爆风荷载响应分析
爆风荷载响应的理论基础
爆风荷载是什么
老师，爆风荷载是什么？
爆炸（化学炸药、气体爆炸等）产生的冲击波（爆炸波）作用在结构上的压力荷载。具有峰值压力→指数衰减→负压的典型时间波形。
Friedlander波形
爆风压力的标准近似式（改进Friedlander式）：
- $P_s$ — 峰值过压
- $t_d$ — 正压持续时间
- $b$ — 衰减参数
- $p_0$ — 大气压
它以指数方式衰减，同时也转向负压啊。
正压阶段（推力）之后是负压阶段（拉力）。建筑物的墙体受到正压（外向）→负压（内向）的反复作用。
爆风参数的估算
使用ConWep（Conventional Weapons Effects Program）的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方法（Arbitrary Lagrangian-Eulerian） — 用欧拉网格模拟爆风传播，用拉格朗日网格模拟结构变形，同时求解。包含反射和绕射
ALE方法更精确吗？
ALE方法自动计算爆风的反射、绕射和叠加，更精确。但需要三维的空气网格，计算成本很大。对于简单形状，直接输入压力就足够了。
总结
要点：
- Friedlander波形 — 峰值过压→指数衰减→负压
- 用ConWep（Kingery-Bulmash式）估算参数 — $Z = R/W^{1/3}$
- 直接输入压力或ALE方法 — 根据复杂程度选择
- LS-DYNA的*LOAD_BLAST_ENHANCED — ConWep自动计算
爆炸压力随距离的三次方衰减
TNT当量换算缩放规律（Hopkinson-Cranz规律）规定，爆炸最大过压由距离R与药量W的立方根之比Z＝R/W^(1/3)整理，Z＝1时最大过压约为0.3MPa。1944年Brode式后来精化为Kingery-Bulmash数据库（美军TM5-855-1），目前仍被ANSYS Autodyn的ATBLAST功能采用。
爆风荷载响应的数值计算方法
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（Arbitrary Lagrangian-Eulerian）方法是主流，空气、炸药用欧拉网格处理，结构用拉格朗日单元处理。LS-DYNA的ALE多物质代码在1990年代由LLNL（劳伦斯·利弗莫尔国家实验室）开发，1 kg TNT的自由场爆炸在100万节点、5ms内的模拟目前用HPC约需30分钟。
爆风荷载响应的实际应用
爆风分析的实践
在防卫、石油化工厂、反恐结构设计中使用。
应用案例
| 应用 | 目的 |
|---|---|
| 防爆墙设计 | 隔挡爆风。将变形量控制在允许值内 |
| 建筑耐爆设计 | 玻璃破坏、结构响应 |
| 车辆耐爆性能 | 对IED（简易爆炸装置）的耐性 |
| 厂房安全距离 | 爆炸事故时对结构的影响 |
实践检查清单
UFC 3-340-02是什么？
美国国防部对爆炸荷载的结构设计手册。规定了峰值压力、持续时间、结构允许变形。防爆设计的国际标准。
防护设计的基准是UFC 4-010-01
美国防部UFC（统一设施规范）4-010-01是政府设施爆炸防护设计基准，规定了防护距离和结构响应限值。2001年9·11事件后，2012年修订版要求将RC板的韧性比（μ）控制在10以内，并明确要求使用动态分析（SDOF或FEM）进行验证。
爆风荷载响应的软件对比
爆风分析的工具
选型指南
爆炸分析中LS-DYNA和Autodyn并驾齐驱
爆炸冲击分析市场中LS-DYNA（Ansys）和ANSYS Autodyn是两强。Autodyn由Century Dynamics公司于1986年开发，2005年被Ansys收购。SPH、ALE、拉格朗日混合分析在该领域特别评价高。Abaqus/Explicit在军需防护领域采用例较少，但在汽车OEM的爆炸阀分析中广泛使用，其CEL耦合功能有效。
爆风荷载响应的先端研究
SPH方法的爆风模拟
近距离爆炸时空气网格变形很大。用SPH方法将空气作为粒子处理，可避免网格扭曲问题，相关研究活跃。
多物质ALE
炸药→爆炸生成气体→空气→结构的多物质ALE方法可模拟爆炸全过程。用JWL状态方程描述炸药膨胀。
TBI（外伤性脑损伤）的模拟
用人体有限元模型（HBM）模拟爆风对人体（特别是脑）的损伤。爆风→颅骨变形→脑部压力传播→TBI链的预测。
总结
CONWEP函数快速计算爆风荷载
LS-DYNA的*LOAD_BLAST_ENHANCED卡片参考CONWEP（Conventional Weapons Effects）数据库，无需ALE流体解析直接将爆风压力映射到结构面。计算成本仅为ALE方法的1/10以下，但在近距离（Z＜0.5 m/kg^(1/3)）或结构-爆炸相互作用较大的情况下误差增加，因此用途决定选择很重要。
爆风荷载响应的故障处理
ConWep压力不正确
ALE方法中空气网格不正确
总结
忽略负压阶段的设计很危险
爆炸压力在正压阶段之后必然产生负压阶段（吸引压）。玻璃窗向外破裂事故多数由此负压（－0.01〜－0.05 MPa）引起，但简易SDOF设计中常被忽略。Autodyn和LS-DYNA的CONWEP实现默认截断负压，需明确开启*LOAD_BLAST_ENHANCED的NEGATIVE标志。
错误