波动基本性质 — 波长・频率・干涉与声学有限元分析
1. 波动与CAE的联系
超声波探伤是怎么工作的？听说可以检测金属内部的裂缝，但看不见摸不着的裂缝怎么检测到的？
利用声波在不同界面处会反射的原理。探头发出超声波（MHz频率），在材料内部传播，遇到裂缝（实际上是一个气-固界面）时声阻抗突变，一部分波被反射回来。通过测量回波的到达时间和幅值，就能推算裂缝的深度和大小。就像蝙蝠用超声波"看"障碍物一样。
波动现象在CAE中的应用范围极广：弹性波传播（冲击响应）、声学分析（NVH）、超声探伤、电磁波（天线设计）、地震波（建筑抗震）。掌握波动基础是理解这些高级分析的前提。
2. 波的分类：横波与纵波
按振动方向与传播方向的关系：
| 波型 | 振动方向 | 传播介质 | 工程应用 |
|---|---|---|---|
| 纵波（P波） | 平行于传播方向 | 固、液、气 | 超声探伤（直探头） |
| 横波（S波） | 垂直于传播方向 | 固体（不能在流体中传播） | 超声探伤（斜探头）、地震波 |
| 表面波（Rayleigh波） | 椭圆运动（P+S耦合） | 固体表面 | 地震破坏、SAW传感器 |
| 导波（Lamb波） | 板中复杂振动模态 | 薄板 | 板材大面积无损检测 |
在固体中，P波（纵波）速度比S波（横波）快，地震学中"P先到，S后到"用于定位震源。
3. 波动基本参数
描述波动的基本参数：
| 参数 | 符号 | 单位 | 物理含义 |
|---|---|---|---|
| 波速 | $c$ | m/s | 波形传播的速度 |
| 频率 | $f$ | Hz | 单位时间内完整振动次数 |
| 波长 | $\lambda$ | m | 一个完整波形的空间长度 |
| 周期 | $T$ | s | 完成一次完整振动的时间 |
| 圆频率 | $\omega = 2\pi f$ | rad/s | 角速度形式的频率 |
| 波数 | $k = 2\pi/\lambda$ | rad/m | 单位长度内的相位变化量 |
平面行波的数学表达：
4. 波动方程
一维波动方程（弹性体/声波）：
三维波动方程：
在弹性体中，应力波的传播遵循波动方程，对于各向同性材料：
其中 $\lambda, \mu$ 是Lamé常数，$\rho$ 是密度。
5. 不同介质中的波速
| 介质 | 纵波速度 (m/s) | 横波速度 (m/s) | 声阻抗 Z = ρc (MRayl) |
|---|---|---|---|
| 空气（20°C） | 343 | — | 0.00041 |
| 水（25°C） | 1496 | — | 1.49 |
| 铝 | 6320 | 3130 | 17.1 |
| 钢（结构） | 5920 | 3250 | 46.5 |
| 钛合金 | 6070 | 3140 | 27.0 |
| 混凝土 | 3000-4500 | — | 7.2-10.8 |
声阻抗差异是波在界面处反射的根本原因。反射系数 $R$ 和透射系数 $T$：
钢-空气界面的反射系数：$R = (0.00041 - 46.5)/(0.00041 + 46.5) \approx -0.9999998$——几乎完全反射！这就是超声探伤能检测内部裂缝（气体夹杂）的物理基础。
6. 反射、折射与透射
波在界面处的行为遵循Snell定律：
其中 $\theta_i$ 是入射角，$\theta_r$ 是反射角，$\theta_t$ 是折射角。
当 $c_2 > c_1$ 时（波速从慢介质进入快介质），存在临界角 $\theta_c = \arcsin(c_1/c_2)$，超过临界角发生全反射。这一原理在超声探伤的斜探头设计中被广泛利用。
超声探伤用斜探头是为了什么？正对着零件打进去不好吗？
竖向裂缝就需要斜探头。直探头的波束垂直射入，只能反射回水平方向的缺陷。但焊缝里常见的是竖向或斜向裂缝，波束要打到它才能有反射信号。斜探头把波束以一定角度射入，经过折射后在材料内部斜向传播，正好能"扫到"这些方向的裂缝。探伤工程师要根据缺陷取向来选择探头角度——常用45°、60°、70°。
7. 超声波探伤原理
超声波无损检测（UT, Ultrasonic Testing）的工作原理：
- 探头发出超声脉冲（频率1~20 MHz）
- 声波在工件内部传播
- 遇到缺陷（裂缝、气孔、夹杂）时产生回波
- 探头接收回波，测量传播时间 $t$ 和幅值 $A$
- 计算缺陷深度：$d = c \cdot t / 2$（往返距离）
缺陷深度计算：
相控阵超声（PAUT）通过控制多个晶片的时序，可以实现波束转向和聚焦，相当于超声CT扫描。目前已广泛用于压力容器、核电管道、航空复合材料的检测。
| 探伤方法 | 频率范围 | 可检最小缺陷 | 适用材料 |
|---|---|---|---|
| 普通脉冲反射 | 1-10 MHz | ~1 mm² | 金属、陶瓷 |
| 高频超声 | 10-50 MHz | ~0.1 mm² | 薄件、涂层 |
| 相控阵（PAUT） | 1-20 MHz | ~0.5 mm² | 焊缝、复合材料 |
| 全矩阵捕获（FMC/TFM） | 1-10 MHz | ~0.3 mm² | 复杂几何件 |
8. 声学有限元：Helmholtz方程
声压场满足Helmholtz方程（频域波动方程）：
声学FEM将声压 $p$ 作为节点未知量，单元刚度矩阵和质量矩阵：
频域方程（$\omega^2[M_a] = k^2[K_a]$）：
声学FEM关键的建模要点：
- 每个波长至少需要6个网格节点（6PPW，Points Per Wavelength）
- 高频分析（>1000 Hz）网格量急剧增大，常改用统计能量法（SEA）
- 吸声材料用Delany-Bazley等多孔介质模型
- 无反射边界条件（吸收边界）防止计算域边界的人为反射
9. 汽车NVH仿真
NVH（Noise, Vibration, Harshness）是汽车舒适性的核心指标。CAE在NVH开发中的流程：
| 频段 | 分析方法 | 关注问题 | 软件 |
|---|---|---|---|
| 0-100 Hz（低频） | 结构FEM + 声学FEM | 车身共振、低频轰鸣 | Nastran, Abaqus |
| 100-500 Hz（中频） | 混合FEM/SEA | 路噪、风噪 | VA One（ESI） |
| 500+ Hz（高频） | 统计能量法（SEA） | 高频噪声传递 | VA One, AutoSEA |
典型汽车NVH设计目标（驾驶室内）：
- 怠速（700 rpm）：< 40 dB(A)
- 100 km/h巡航：< 65 dB(A)
- 全速加速：< 75 dB(A)
- 变速器啸叫（Whine）：主观评分 > 7/10
总结
波动现象（弹性波、声波）在CAE中的应用从超声无损检测到汽车NVH，从建筑声学到航空结构颤振，无处不在。理解波长、波速、反射和干涉等基本概念，是进入声学FEM、弹性波仿真等高级分析领域的必要基础。