波动与声学基础
波动方程、声阻抗、NVH与结构声学耦合仿真
声学与波动力学涵盖了从喷气发动机轰鸣到豪华轿车车厢内几乎察觉不到的嗡嗡声等一切现象。在现代CAE实践中，噪声和振动分析——统称NVH——是重要的设计驱动因素，受到的工程关注丝毫不亚于结构强度分析。从物理第一原理出发理解这些问题，使工程师能够准确诊断噪声问题、设计更安静的系统，并对仿真结果充满信心。
1. 声波方程
教授，声波方程是从哪里来的？它跟Navier-Stokes方程是完全独立的东西，还是从那里推导出来的？
是从Navier-Stokes方程推导出来的——它是完整N-S方程的线性化、无粘、无旋简化版本。推导步骤是：从连续性方程和动量方程出发，假设扰动量相对于静止平均态很小（无平均流），忽略粘性和非线性项，最终得到声压的经典波动方程。这个线性化是关键：声学研究的是小振幅波的物理。一旦振幅很大——比如在喷气发动机喷管附近——就需要非线性声学，简单波动方程就不再适用了。
1.1 从流体力学推导
将压力和密度分解为均值与扰动：$p = p_0 + p'$，$\rho = \rho_0 + \rho'$，速度 $\mathbf{v} = \mathbf{v}'$（零均值流）。对连续性方程和欧拉动量方程线性化：
利用等熵关系 $p' = c^2 \rho'$，对动量方程取散度并结合连续性方程对时间的导数：
这就是经典线性声波方程——二阶双曲型偏微分方程。对位移势 $\phi$（其中 $\mathbf{v}' = \nabla\phi$）具有相同的形式。
1.2 亥姆霍兹方程（频域）
对于时谐场 $p'(\mathbf{x},t) = \hat{p}(\mathbf{x})e^{i\omega t}$，波动方程退化为亥姆霍兹方程：
其中 $k$ 为波数。亥姆霍兹方程在空间域是椭圆型偏微分方程，对于窄频带问题比时域波动方程更易数值求解。大多数商业声学FEM软件（ACTRAN、Comsol Acoustics、Nastran SOL 108）均在频域工作。
1.3 含源波动方程
当存在体积声源时（如振动边界作为单极子声源）：
其中 $q(\mathbf{x},t)$ 为单位体积质量源率。自由空间Green函数解为：
2. 声速与声压级
为什么声音在水中比在空气中传播更快？水明显更重，按直觉应该更慢才对……
直觉不错，但密度只是故事的一半。声速取决于恢复力（体积弹性模量 $K$）与惯性（密度 $\rho$）的比值：$c = \sqrt{K/\rho}$。水的体积弹性模量约2.2 GPa——极其坚硬，几乎不可压缩。空气的体积弹性模量只有142 kPa。虽然水的密度是空气的800倍，但其刚度约为空气的15000倍，所以 $K/\rho$ 在水中约大20倍。钢更夸张——声音在钢中约5000 m/s，因为钢的刚度极高。
2.1 声速公式
对于理想气体，等熵（绝热）体积弹性模量 $K_s = \gamma P$。20℃空气中：$c = \sqrt{1.4 \times 287 \times 293} = 343$ m/s。
| 介质 | $c$ [m/s] | $\rho_0$ [kg/m³] | $Z = \rho_0 c$ [Pa·s/m] |
|---|---|---|---|
| 空气（20℃，1 atm） | 343 | 1.204 | 413 |
| 水（20℃） | 1482 | 998 | 1.48 × 10⁶ |
| 钢 | 5960 | 7850 | 4.68 × 10⁷ |
| 铝 | 6420 | 2700 | 1.73 × 10⁷ |
| 混凝土 | 3100 | 2300 | 7.1 × 10⁶ |
| 软组织 | 1540 | 1060 | 1.63 × 10⁶ |
2.2 声压级（SPL）
声压级使用以人耳听觉阈值 $p_{\text{ref}} = 20\,\mu$Pa 为基准的对数刻度：
参考值：0 dB = 听觉阈值；65 dB = 正常对话；85 dB = 长期暴露可致听力损伤；130 dB = 痛觉阈值；194 dB = 空气中无畸变声音的理论上限。
3. 声阻抗与反射
比声阻抗为 $Z = \rho_0 c$ [Pa·s/m = 瑞利]。在介质1与介质2的交界面处（法向入射），声压反射系数为：
声强透射系数为：
做B超的时候为什么要在探头和皮肤之间涂耦合剂？感觉很奇怪……
这正是声阻抗失配的问题。空气的 $Z \approx 413$ Pa·s/m，软组织的 $Z \approx 1.63 \times 10^6$ Pa·s/m。代入透射系数公式：$4 \times 413 \times 1.63\times10^6 / (413 + 1.63\times10^6)^2 \approx 0.001$——不到0.1%的声能能透过去！耦合凝胶的声阻抗与组织接近，消除了空气间隙，实现近乎完美的透射。这个阻抗匹配原理同样用于声学换能器设计、建筑隔声设计和声纳系统。
3.1 斜入射与斯涅尔定律
对于以角度 $\theta_i$ 入射界面法线的波，斯涅尔定律成立：
当 $c_2 > c_1$ 时，存在全反射临界角 $\theta_c = \arcsin(c_1/c_2)$，超过此角时无能量透射——类似光的全内反射。这一现象用于声学波导设计。
4. 驻波与共振
我们的车厢在50 Hz附近有明显的轰鸣声，声学仿真显示有驻波。怎么判断这是空腔声学模态还是结构模态引起的？
先做个快速估算：矩形腔第一纵向声学模态约为 $f \approx c/(2L)$，$L$ 取最长尺寸。典型车厢约2米长，$f \approx 343/(2\times2) = 86$ Hz。你的50 Hz轰鸣低于此值，说明很可能是结构-声学耦合模态——车身白板件参与了振动。做振动声学联合分析：先看50 Hz附近有哪些结构模态，再查看这些模态振型是否能高效地把能量泵入空腔。如果发现地板在50 Hz附近与空腔声压同相振动，那就是罪魁祸首。对该板件进行加强或增加阻尼是通常的解决办法。
4.1 一维管道驻波
对于两端为刚性壁（速度节点）、长度为 $L$ 的管道，共振频率为：
压力振型：$\hat{p}_n(x) = A\cos(n\pi x/L)$。壁面处为压力波腹，对 $n=1$ 中心处为压力波节。
对于一端开口的管道（四分之一波长共振器），仅存在奇次谐波：
4.2 三维声学空腔模态
对于 $L_x \times L_y \times L_z$ 的矩形腔，固有频率为：
其中 $m, n, l = 0, 1, 2, \ldots$（不全为零）。$(0,0,0)$ 模态不存在——那将意味着无波动的均匀压力变化。
5. 结构声学与振动声学
结构振动时，有些模态是"高效声辐射体"，有些却不是。是什么决定了一个结构模态的声辐射效率？
关键概念是辐射效率，更准确地说是结构波数与声学波数是否匹配。这样想：一个结构模态有特定的位移空间分布。如果板的相邻区域反向运动（异相），它们辐射的声音会相互抵消——这个模态"短路"了声辐射。但如果结构波长远大于该频率下的声波长（即结构是"超声速"的），辐射无法抵消，模态就成为高效辐射体。这就是吻合频率（临界频率）的概念——低于此频率，辐射效率低；高于此频率，板可以成为非常高效的辐射体。
5.1 辐射效率
对于障板中的板，辐射效率 $\sigma$ 表征结构模态与等面积活塞相比的辐射能力：
其中 $W_{\text{辐射}}$ 为辐射声功率，$\langle v^2 \rangle$ 为均方表面速度。大多数结构模态在临界频率以下时 $\sigma \ll 1$。
5.2 吻合频率（临界频率）
对于厚度为 $h$ 的薄板，吻合频率为：
对于2 mm钢板：$f_c \approx 2500$ Hz——高于此频率，钢板辐射效率很高。对于50 mm混凝土墙：$f_c \approx 120$ Hz——低于此频率，隔声性能差。这解释了为什么混凝土墙对低频低音的隔绝效果很差。
5.3 结构-声学耦合方程
结构位移 $\mathbf{u}$ 与声压 $p$ 的耦合方程为：
耦合矩阵 $\mathbf{C}$ 代表共享界面处的流固相互作用：声压对结构加载（声力 = $p \cdot \mathbf{n} \cdot A$），结构速度驱动流体（界面处结构法向速度 = 声学质点速度）。
6. NVH：噪声、振动与声振粗糙度
NVH（Noise, Vibration, Harshness，噪声-振动-声振粗糙度）是汽车和航空航天工程中研究声学和振动响应主观质量的学科。"噪声"指空气传声，"振动"指结构传振感觉，"声振粗糙度"指令人不适的高频结构粗糙感，而非声音本身。
6.1 源-路径-接收者框架
每个NVH问题都可分解为三要素：
- 声源： 发动机燃烧、路面不平度、风噪、电机谐波
- 传递路径： 结构传声（通过悬置、底盘）或空气传声（通过缝隙、车身板件）
- 接收者： 驾乘人员耳部，或法规测量麦克风
NVH的CAE工作系统地针对每个环节：通过封闭力测量进行声源表征、传递路径分析（TPA）量化各路径贡献、接收端灵敏度分析。
6.2 传递路径分析（TPA）
第 $k$ 个声源经第 $j$ 条路径传递至接收端的声压为：
其中 $H_j(\omega)$ 为路径 $j$ 的声学传递函数（Pa/N），$F_j$ 为该路径的封闭力。TPA使工程师能够按贡献量排列各路径的优先级，仅针对主导路径采取治理措施。
6.3 声品质评价指标
| 指标 | 单位 | 说明 |
|---|---|---|
| A计权声压级 | dB(A) | 按人耳听觉灵敏度加权的总声级 |
| 响度 | sone | 心理声学响度（Zwicker法） |
| 尖锐度 | acum | 高频能量集中程度 |
| 粗糙度 | asper | 15~300 Hz幅值调制程度 |
| 音调度 | tu | 离散纯音相对于背景噪声的突出程度 |
7. 声学FEM与BEM仿真
同事说对于外场声学问题，BEM比FEM更好。为什么外场问题BEM更有优势？
关键在于索末菲辐射条件。在外场声学中，声场延伸至无穷远——没有可以用作边界的远处壁面。用FEM时，你必须对有限体积的空气建立网格，然后在截断边界上施加吸收边界条件（如完美匹配层PML）来防止反射——这增加了计算量，PML本身也会引入近似误差。BEM直接利用满足整个无界域波动方程的Green函数进行面积分——只需对结构表面划网格，无需对周围空气建模。计算汽车车身振动辐射声场这类问题，BEM效率高得多。对于乘客舱等内腔问题，FEM更好，因为BEM在内场问题上有困难。
7.1 声学FEM弱形式
亥姆霍兹方程的弱形式，含结构耦合边界条件 $\partial p/\partial n = -\rho_0\omega^2 u_n$：
声学网格的关键规则：
- 每个声波波长 $\lambda = c/f$ 至少需要6个线性单元（或3个二次单元）
- 1000 Hz空气中：$\lambda = 343/1000 = 0.343$ m → 最小单元约57 mm（线性）
- 10 kHz时需要6 mm单元——高频声学FEM计算量极大！
7.2 BEM公式
外场问题的BEM积分方程（Kirchhoff-Helmholtz积分方程）：
其中 $G(\mathbf{x},\mathbf{y}) = e^{-ik|\mathbf{x}-\mathbf{y}|}/(4\pi|\mathbf{x}-\mathbf{y}|)$ 为自由空间Green函数，光滑边界点处 $c(\mathbf{x}) = 1/2$。BEM组装的是全填充（稠密）矩阵，大型问题需使用迭代求解器和快速多极子方法（FMM）。
7.3 高频段统计能量分析（SEA）
对于全车模型，FEM和BEM在约2~4 kHz以上时因网格尺寸要求变得不切实际。统计能量分析将系统建模为相互交换振动和声学能量的子系统网络：
其中 $E_i$ 为子系统能量，$n_i$ 为模态密度，$\eta_i$ 为损耗因子，$\eta_{ij}$ 为耦合损耗因子。SEA给出频带平均能量水平而非空间细节，但对高频问题计算成本极低。常用工具：VA One（ESI）、AutoSEA、ANSYS Mechanical SEA扩展模块。