水中爆発解析
理论与物理
水中爆炸（UNDEX）是什么
老师，水中爆炸分析是什么样的问题？
UNDEX（Underwater Explosion）是水下爆炸冲击波作用于舰船或潜艇船体的问题。这是海军舰体设计中最重要的抗冲击评估。
水中爆炸的物理
水中爆炸的特征现象：
1. 冲击波（primary shock wave） — 从爆源球面状传播。以水中音速（1500 m/s）移动
2. 气泡脉动 — 气球的膨胀→收缩→再膨胀的重复。低频压力脉动
3. 空化 — 冲击波反射在水面附近产生负压→水沸腾（空化）→崩溃时产生二次冲击
水中和空气中不一样呢。
水的密度是空气的800倍。冲击波的能量大一个数量级，相同TNT当量下比空气中造成数倍的损伤。
FEM中的分析方法
DAA法能在Nastran里用吗？
在Nastran的SOL 109/112中定义FLUIDEX（流体外部场），用DAA近似将水中爆炸的冲击波输入到结构。这是美国海军的标准方法。
总结
要点：
- 水中爆炸比空气中能量大一个数量级 — 水的密度是800倍
- 冲击波 + 气泡脉动 + 空化 — 三种现象
- DAA法（Nastran）、ALE法（LS-DYNA）、BEM-FEM — 分析方法
- 海军舰体设计中最为重要 — 抗冲击性能评估
水中爆炸的气泡溃灭破坏结构
水中爆炸（UNDEX）中，在初始冲击波之后，爆炸气体气泡会反复膨胀收缩的“气泡脉动”对结构物造成二次冲击。这个气泡脉动的周期可以用 T＝K(W^(1/3)/(D+10)^(5/6)) 秒（K是TNT常数，D是爆炸深度m）来估算，典型的爆炸（TNT 100kg，水深30m）周期约为0.5秒。第二次世界大战中许多击沉舰船的案例后来发现是由这种气泡脉动导致的龙骨断裂。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
DAA法的实现
DAA（双渐近近似）是同时近似低频附加质量效应和高频声辐射的方法。在结构的湿表面（与水接触的面）设置DAA边界条件。
在Nastran的哪个求解序列中使用？
在SOL 109/112中添加FLUIDEX（外部流体场）。爆炸冲击波用Cole公式等计算，作为入射波输入到结构面。散射波用DAA近似计算。
ALE法（LS-DYNA）
水和炸药用欧拉（ALE）网格，船体用拉格朗日网格同时求解。直接模拟冲击波的传播、反射、空化全过程。
ALE法更准确吗？
ALE法最准确，但水的网格量巨大（数百万至数千万单元）。远距离爆炸用DAA法更高效。近距离则需要ALE法。
总结
水和炸药用欧拉（ALE）网格，船体用拉格朗日网格同时求解。直接模拟冲击波的传播、反射、空化全过程。
ALE法更准确吗？
ALE法最准确，但水的网格量巨大（数百万至数千万单元）。远距离爆炸用DAA法更高效。近距离则需要ALE法。
DAA（延迟声学近似）是流固耦合的关键
UNDEX分析的流体-结构耦合广泛使用DAA（双渐近近似）来近似无限水域的散射波。1978年DeRuntz和Geers开发的DAA2（二阶精度）连续切换远场和近场的渐近解，被美国海军的舰艇UNDEX认证代码USA（水下冲击分析）采用。USA代码作为Nastran或Abaqus的外部流体求解器，至今仍由美国海军研究实验室（NRL）管理。
线性单元（一阶单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切锁定（可通过减缩积分或B-bar法缓解）。
二阶单元（带中间节点）
可以表现曲线变形。应力精度大幅提高，但自由度约增加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 迭代法的比喻
直接法是“用笔算精确解联立方程”的方法——可靠但大规模问题太耗时。迭代法是“反复猜测逼近正确答案”的方法——最初是粗略的
相关主题
なった
詳しく
報告