3. 结构声耦合的数学描述

结构声耦合问题的数学描述，本质上要回答三件事：结构在流体载荷下如何振动、流体在结构振动下如何响应、两者在交界面上如何交换能量。前两章我们分别建立了结构动力学与声场的基本方程，本节把它们“焊”在一起，形成耦合波动方程、边界条件与界面连续性条件，并给出弱耦合与强耦合的判别准则。理解这套数学框架，是后续有限元/边界元建模、模态耦合分析与噪声预测的起点。

一、耦合波动方程：两个场，一个界面

设结构域为 Ω_s，流体域为 Ω_f，二者共享耦合界面 Γ。结构位移场 u(x,t) 满足弹性动力学方程：

ρ_s ∂²u/∂t² − ∇·σ(u) = f_s      (在 Ω_s 内)

其中 ρ_s 为结构密度，σ 为应力张量，f_s 为外加体积力。流体采用声压 p(x,t) 描述时，满足亥姆霍兹型波动方程：

(1/c²) ∂²p/∂t² − ∇²p = 0          (在 Ω_f 内)

c 为流体声速。若考虑流体中的体积声源 q，则右端为源项。两个方程本身并不耦合，耦合发生在界面 Γ 上：结构表面法向加速度驱动流体，流体压力又反作用于结构表面。因此，耦合波动方程组的完整形式应写成“结构方程 + 流体方程 + 界面耦合项”三部分，缺一不可。

二、三类边界条件

边界条件是方程定解的关键。结构声耦合中常见的边界有三类，它们分别约束不同的物理量。

弹性边界：结构边缘连接弹簧或支承结构，边界上力与位移通过刚度建立关系。例如沿法向 n 的弹性支承可写为：

σ_nn = −k_s u_n      (弹性边界，k_s 为分布刚度)

当 k_s→∞ 时退化为固定边界 u_n=0；k_s→0 时退化为自由边界 σ_nn=0。弹性边界常用于模拟减振器、弹性安装座或柔性连接。

阻抗边界：边界上压力与法向质点速度成比例，用于描述吸声材料、局部反应型衬里或辐射阻抗：

∂p/∂n = −(jωρ_f / Z) p      (阻抗边界，Z 为声阻抗)

Z 为实数时表示纯阻性吸收，Z 为复数时同时包含惯性或弹性效应。阻抗边界是噪声控制中吸声处理的核心数学表达。

辐射边界：当流体域在外部无限延伸时，必须保证声波向外传播而不反射回来，即满足 Sommerfeld 辐射条件：

lim (r→∞) r (∂p/∂r + jk p) = 0

在数值计算中，辐射边界通常用无限元、边界元或完美匹配层（PML）近似实现。若忽略辐射边界，外场计算会出现虚假反射，导致辐射声功率预测严重失真。

三、耦合界面连续性条件

界面 Γ 上必须同时满足两个连续性条件，它们是把结构与流体“粘”在一起的核心方程。

运动连续性：结构与流体在法向不可分离、不可穿透，因此结构法向加速度等于流体法向加速度：

∂p/∂n = −ρ_f ∂²u_n/∂t²      (运动连续性)

力连续性：界面两侧法向力大小相等、方向相反，即流体压力以面载荷形式作用于结构：

σ_n = −p n      (力连续性，n 为结构外法向)

注意符号约定：本文取结构外法向 n 指向流体域，因此流体压力对结构的作用为 −p n。若采用相反法向定义，符号需整体调整，但物理结论不变。这两个条件联立后，结构方程与流体方程不再独立，形成对称或非对称的耦合系统矩阵。

四、弱耦合与强耦合的判别准则

耦合强弱并非主观判断，而由流体附加质量与结构自身惯性的比值决定。定义无量纲耦合参数：

β = ρ_f c / (ρ_s h ω)

其中 h 为结构特征厚度，ω 为角频率。β 的物理含义是“流体载荷引起的等效附加质量”与“结构单位面积质量”之比。当 β≪1（例如空气中钢板，ρ_f≈1.2 kg/m³，ρ_s≈7850 kg/m³，β 通常小于 0.01），流体对结构振动的反作用可忽略，可先算结构振动、再算声辐射，即单向弱耦合。当 β≥1（例如水中薄板或轻质复合材料结构），流体附加质量与结构质量同量级，结构模态频率和振型发生明显漂移，必须双向迭代求解，即强耦合。

除 β 外，还有两个辅助判据：一是模态密度，当结构模态与流体腔模态频率接近时，即使 β 不大也可能出现强耦合共振；二是辐射效率，当结构尺寸与声波波长可比时，辐射阻尼显著增大，耦合增强。工程中常用“频率漂移量”作为最终判据：若考虑流体后某阶模态频率变化超过 5%，则应按强耦合处理。

五、教学示例：水中钢板耦合频率估算

取一块 1 m×1 m、厚 5 mm 的钢板，四边简支，置于水中。钢材 ρ_s=7850 kg/m³，E=210 GPa，ν=0.3；水 ρ_f=1000 kg/m³，c=1500 m/s。钢板一阶弯曲模态在真空中的频率约为：

f_0 ≈ (π/2) (1/a² + 1/b²) √(D/(ρ_s h))
D = E h³ / [12(1−ν²)] ≈ 210e9 × (0.005)³ / [12×(1−0.09)] ≈ 2.40 N·m
f_0 ≈ (π/2) × 2 × √(2.40 / (7850×0.005)) ≈ 3.14 × 0.247 ≈ 0.78 Hz

水中附加质量近似为 ρ_f a / (π n) 量级，对一阶模态附加质量约 1000×1/π≈318 kg/m²，而钢板面密度仅 39.3 kg/m²，附加质量是结构质量的 8 倍以上，β 远大于 1。因此该板在水中的一阶耦合频率会显著下降，实际应通过声固耦合有限元求解，预计频率降至 0.3 Hz 量级。这个例子说明：轻质薄壁结构在水下必须按强耦合处理，否则噪声预测会严重偏高。

六、从方程到数值模型

把上述方程离散后，结构域得到 M_sü + K_su = F_s + F_c，流体域得到 M_fp̈ + K_fp = F_f，界面耦合项通过形函数积分形成耦合矩阵 C，最终组装为：

[ M_s   0  ] { ü }   [ K_s   C  ] { u }   { F_s }
[ ρ_f Cᵀ M_f ] { p̈ } + [ 0   K_f ] { p } = { F_f }

该矩阵通常是非对称的，求解时需注意：弱耦合可忽略 C 的反馈项，强耦合必须保留完整耦合矩阵并采用迭代或直接耦合求解器。理解这一矩阵结构，就能在软件中正确设置耦合面、辐射边界与阻抗边界，并对结果是否可信做出判断。