2. 声学与结构动力学基础回顾

结构声耦合建模的本质，是把“结构在振动”和“流体在传播”这两套物理描述放在同一个界面上求解。要做到这一点，必须先统一三件事：声波方程描述流体中的声场、结构振动模态描述弹性体的固有运动形态、量纲与单位体系保证两侧变量可以正确耦合。本节围绕这三条主线展开，为后续耦合界面条件与数值建模打基础。

一、声波方程：从三个守恒关系到一个波动方程

在均匀、静止、无黏、绝热的理想流体中，声学基本方程由三部分组成：运动方程（欧拉方程）、连续性方程和状态方程。对小幅扰动做线性化后，可消去质点速度，得到以声压 p 为未知量的齐次波动方程：

∇²p − (1/c₀²)·∂²p/∂t² = 0

其中 c₀ = √(γ·p₀/ρ₀) 为流体中的声速
空气常温常压下 c₀ ≈ 343 m/s，ρ₀ ≈ 1.21 kg/m³

若存在声源，则右端出现源项，成为非齐次方程。对简谐激励，令 p(x,t)=P(x)·e^{jωt}，波动方程化为亥姆霍兹方程 ∇²P + k²P = 0，其中波数 k = ω/c₀ = 2π/λ。这一步把时域问题转为频域问题，是结构声耦合频域分析的标准入口。

二、声压与质点速度：两个必须同时掌握的场量

声压 p 是标量，表示扰动引起的压强变化；质点速度 v 是矢量，表示流体微团的运动速度。二者并非独立，而由线性化欧拉方程联系：

ρ₀·∂v/∂t = −∇p

简谐情况下：v = −(1/(jωρ₀))·∇p

这个关系说明：质点速度由声压的空间梯度决定，方向沿声压下降最快的方向。在平面波中，二者同相，比值等于特性阻抗 ρ₀c₀（空气约 415 Pa·s/m）；在近场或驻波中，二者存在相位差，这正是耦合界面附近必须用复数描述的原因。耦合面上，流体一侧的法向质点速度必须等于结构表面的法向振速，这是声学与结构之间最关键的连续性条件。

三、结构振动模态与模态叠加法

对线性弹性结构，无阻尼自由振动可写成广义特征值问题：

K·φᵢ = ωᵢ²·M·φᵢ

K：刚度矩阵   M：质量矩阵
φᵢ：第 i 阶模态振型   ωᵢ：第 i 阶固有圆频率

模态振型具有关于质量矩阵和刚度矩阵的正交性，因此可以把物理坐标下的耦合方程组投影到模态坐标上，实现解耦。模态叠加法的核心表达式为：

u(x,t) = Σᵢ φᵢ(x)·qᵢ(t)

q̈ᵢ + 2ζᵢωᵢq̇ᵢ + ωᵢ²qᵢ = Fᵢ(t)/mᵢ

mᵢ = φᵢᵀ·M·φᵢ（模态质量）
Fᵢ = φᵢᵀ·F（模态力）

实际使用中，模态截断阶数直接决定高频段精度。经验上，若关注频段上限为 f_max，至少应保留固有频率达到 1.5～2 倍 f_max 的模态；对薄板、壳体等密集模态结构，还需检查模态参与因子，避免遗漏主导模态。

四、量纲与单位体系：耦合计算中最易出错的一环

结构声耦合涉及力学量与声学量两套单位，混用会导致数量级错误。下表列出常用量的量纲与 SI 单位，便于建模时逐项核对。

需要特别注意的是：有限元建模中若长度用 mm、力用 N、密度用 t/mm³，则弹性模量应取 MPa，此时声压单位会变为 MPa 量级，必须显式换算回 Pa 才能与声学结果对比。建议在模型建立之初就固定一套一致单位制，并在后处理阶段统一转换到 dB 或 Pa 输出。

五、耦合视角下的变量对应关系

把上述内容串起来看，结构声耦合在界面上完成一次“变量交接”：结构侧输出法向位移或法向振速，流体侧把它作为声场边界条件；流体侧输出声压，结构侧把它作为分布载荷。二者通过界面法向连续性条件闭合。下面的关系图给出这一交接过程及关键变量、单位与判断依据。

图中左侧是结构侧从特征值问题到模态叠加、再到法向振速与载荷的链条，右侧是声学侧从亥姆霍兹方程到质点速度、声压与声压级的链条，中间橙色区域是耦合界面必须满足的两条连续性条件。实际建模时，只要沿这张图逐项核对变量、单位与符号方向，就能避免绝大多数因量纲不一致或界面条件写错导致的错误结果。