应力波传播分析
应力波传播的理论基础
应力波传播是指
老师，「应力波」是什么？
当结构受到冲击时，应力（变形）以波的形式在材料内传播。这就是应力波。音速 $c = \sqrt{E/\rho}$ 传播。
弹性波的种类
| 波的种类 | 速度 | 特征 |
|---|---|---|
| 纵波（P波） | $c_L = \sqrt{(\lambda+2\mu)/\rho}$ | 在压缩·拉伸方向传播 |
| 横波（S波） | $c_S = \sqrt{\mu/\rho}$ | 在剪切方向传播。$c_S < c_L$ |
| 瑞利波 | $\approx 0.9 c_S$ | 沿表面传播。地震的S波 |
| 兰姆波 | 分散型 | 在薄板中传播。取决于板厚和频率 |
钢的纵波速度约为 $c_L \approx 5900$ m/s 。
1 m远处在 $0.17$ ms 内到达。在冲击分析中，波的传播时间决定了结构的响应时间。
FEM中的应力波分析
为了准确追踪应力波的网格要求：
$n_{ppw}$ 是每个波长的单元数。线性单元的目标是20，二阶单元是10。
1波长20个单元！为了追踪高频波，需要非常细的网格啊。
所以应力波传播是显式FEM的优势领域。$\Delta t$自动根据波的传播变小，因此波的采样自然得到保证。
应用例
总结
要点：
- 应力波以音速传播 — $c = \sqrt{E/\rho}$
- P波（纵）> S波（横）> 瑞利波（表面） — 速度顺序
- 网格要求严格 — 每波长10～20个单元
- 显式FEM最优 — $\Delta t$自动适应波的追踪
- Hopkinson杆、超声波、地震波 — 主要应用
您是否知道固体中有3种弹性波
在固体中传播的弹性波有纵波（P波：压缩·膨胀，音速约5000m/s@钢）、横波（S波：剪切，约3000m/s@钢）、表面波（Rayleigh波：沿表面传播，约2800m/s@钢）三种。地震学利用P波和S波的速度差来计算震源距离，是地震学的基础知识。1883年喀拉喀托火山爆发的音波绕地球传播3.5圈，这是通过大气传播波计算出来的。在工业上，超声波探伤试验利用P波和S波来检测材料内部缺陷。
应力波传播的数值计算方法
LS-DYNA
```
*CONTROL_TIMESTEP
0.0, 0.6 $ 安全系数0.6（波传播用。通常0.9更小）
*CONTROL_TERMINATION
0.001 $ 1 ms（波的往返时间程度）
```
为什么要降低安全系数到0.6？
在应力波传播中，为了最小化数值分散（波速依赖于网格大小的现象），需要降低CFL条件的安全系数。0.9可能导致数值分散过大。
数值分散问题
使用FEM传播波时，波长越短，速度越慢（数值分散）。物理上，所有频率应以相同速度传播，但FEM的离散化导致速度依赖于波长。
对策：
谱元法
谱元法（Spectral Element Method）采用高阶GLL（Gauss-Lobatto-Legendre）点的单元方法。与常规FEM相比，数值分散极小。地震波传播仿真（SPECFEM3D等）的标准。
总结
有限元法波动分析需要高密度网格
弹性波的FEM分析为了确保精度，每波长至少需要8～10个单元（经验法则）。例如，钢材中100kHz超声波（波长50mm）的精确捕捉需要5mm以下的单元。应用到整个1m的试样，单元数会达到数百万到数千万，计算成本爆炸性增加。因此实务中通常采用关键区域细密化、远场粗化的分级网格，以及结合PML（完全匹配层）边界条件来吸收边界反射。
应力波传播的工程应用
Hopkinson杆试验的仿真
SHPB（分裂Hopkinson压杆）试验是获取材料高速应变速率（$10^2 \sim 10^4$ /s）特性的试验。用FEM仿真入射杆→试样→透射杆的应力波传播，协助试验数据的解读。
超声波NDT的仿真
用FEM仿真超声波探伤（UT）的波传播。从裂纹反射波的图案推断裂纹的大小和位置。FEM与逆问题分析的结合。
工程检查清单
「吸收边界条件」是什么？
防止模型边界处波反射回来的人工边界条件。Lysmer-Kuhlemeyer（粘性边界）和PML（完全匹配层）是代表性方法。边界不存在就可以表示无限空间。
新干线铁路焊接部的超声波检测中的FEM应用
JR东日本利用FEM波动分析来优化铁路焊接部的超声波探伤检查系统。用Abaqus/Explicit分析0.5～5MHz超声波在铁路截面中的传播时的模式转换·散射规律，优化探针的位置和角度，使得传统人工探伤难以检测到的焊接线附近的缺陷检测率在2010年代大幅提高。
应力波传播的软件比较
应力波分析工具
选择指南
存在专门从事波动分析的求解器
除了通用FEM代码（Abaqus·LS-DYNA·ANSYS）外，市场上还存在专门用于波动分析的代码。COMSOL Multiphysics的声学模块强于流体·固体耦合声学分析，用于医用超声探针设计。WAVE3D（EPFL开发）和GPU加速SDG（谱差分法）代码用于地震工程的地基波动分析。Dassault SIMULIA Wave Propagation Analysis（Abaqus内置）的波传播分析被空客用于飞机机身SHM波传播分析验证的公开案例中。
应力波传播的前沿研究
应力波的前沿研究
兰姆波改变SHM（结构健全性监测）
在薄板中传播的兰姆波（板波）有对称模式S₀和反对称模式A₀，从缺陷处的波形变化检测裂纹·腐蚀是SHM（结构健全性监测）技术的基础。波音和空客正在研究开发在飞机外蒙皮埋入压电传感器、常时收发兰姆波的系统，2020年代FEM波形传播仿真与机器学习的结合在缺陷位置和大小自动判定研究中进展迅速。
应力波传播的故障排除
应力波故障
用ABC吸收边界条件消除数值波反射的设置方法
在有限区域FEM进行波动分析时，人工边界处的波反射会回到分析区域造成误差。针对这个问题，可以使用吸收边界条件（ABC）的粘性阻尼器边界（标准方法）或PML（完全匹配层）。LS-DYNA的*BOUNDARY_NON_REFLECTING（NRB）可以吸收P波和S波，但对斜入射波和异种材料边界反射率会增加。Abaqus的*INFINITE ELEMENT精度更高，但设置复杂，建议先用试验模型确认反射率后再用于本分析。
价值
详细
错误