4. 声固耦合有限元建模原理
声固耦合有限元建模的核心任务，是把结构域中的弹性振动与声腔域中的声压波动放在同一个离散系统中求解。结构域通常用位移向量 u 描述，声腔域可以用声压 p 描述，也可以用声学位移 u_a 描述。由此形成两条主流路线：位移-压力格式与位移-位移格式。两者在界面耦合条件、矩阵结构和求解规模上差异明显，选择哪一种格式，直接决定后续声学单元与结构单元如何匹配、耦合矩阵如何组装以及是否需要进行对称化处理。
一、位移-压力格式：最常用的声固耦合列式
位移-压力格式以结构位移 u 和声压 p 为基本未知量。结构域满足弹性动力学方程，声腔域满足亥姆霍兹方程。在流固交界面上，需要同时满足两个条件：一是法向位移连续，即结构表面法向位移等于流体法向位移；二是法向力平衡，即结构表面受到流体声压的作用。将这两个条件代入虚功原理，可得到耦合系统的分块方程：
[ M_s 0 ] { ü } [ K_s -S ] { u } { F_s }
[ ρS^T M_f] { p̈ } + [ 0 K_f] { p } = { 0 }其中 M_s、K_s 为结构质量阵与刚度阵，M_f、K_f 为声腔质量阵与刚度阵，S 为耦合矩阵，ρ 为流体密度。耦合矩阵 S 的物理含义是：结构法向位移在声腔边界上产生等效体积速度，从而激励声压；反过来，声压又对结构表面施加法向载荷。由于 S 与 ρS^T 互为转置关系，系统矩阵整体是非对称的。非对称矩阵会带来两个后果：一是特征值求解需要采用非对称算法，计算量增大；二是模态结果可能出现复模态，物理上对应能量在结构与声腔之间的交换与辐射阻尼。
二、位移-位移格式：对称化的一条路径
位移-位移格式把声腔也当作一种“弹性介质”，用声学位移 u_a 作为未知量。此时声压与声学位移的关系为 p = -ρ c² ∇·u_a，声腔的等效“刚度”来自流体压缩性，等效“质量”来自流体惯性。界面条件直接写成位移连续，即 u·n = u_a·n，不再出现压力与位移的交叉项。系统方程变为：
[ M_s 0 ] { ü } [ K_s K_c ] { u } { F_s }
[ 0 M_a ] { ü_a} + [ K_c^T K_a ] { u_a} = { 0 }其中 K_c 是界面耦合刚度阵，M_a、K_a 是声学质量阵与刚度阵。由于 K_c 与 K_c^T 成对出现，整体矩阵是对称的，可以直接使用标准对称特征值求解器。代价是声学位移场比声压场多出三个自由度，对于大体积声腔，规模会显著增加。因此位移-位移格式更适合声腔较小、结构自由度占主导、且希望利用对称求解器的场景。
三、声学单元与结构单元的匹配
无论采用哪种格式，界面上的声学单元与结构单元必须满足几何匹配和形函数协调。实际建模中，结构常用壳单元或实体单元，声腔常用六面体或四面体声学单元。界面处要求：结构单元的面与声学单元的面共享同一组节点，或者至少保证法向位移插值一致。若结构网格与声学网格在界面上不共节点，就需要引入界面映射或罚函数耦合，把结构节点位移映射到声学积分点上，再组装耦合矩阵。映射会引入额外误差，因此优先推荐共节点匹配。
匹配时还要注意法向方向的一致性。耦合矩阵 S 或 K_c 中会出现界面法向 n，如果相邻声学单元法向定义相反，耦合项会相互抵消，导致声腔无法被激励。工程上常用的检查方法是：在界面施加一个已知法向位移，观察声腔是否产生预期的声压响应；若响应接近零，应优先排查法向方向。
四、耦合矩阵的组装与对称化处理
耦合矩阵的组装通常分三步完成。第一步，遍历界面上的所有声学单元面，计算每个面的形函数矩阵 N_s 和 N_f（或 N_a）；第二步，在面内做高斯积分，得到局部耦合矩阵 S_e = ∫ N_s^T n N_f dΓ；第三步，按照节点自由度编号，把 S_e 累加到全局耦合矩阵 S 的对应位置。由于界面节点同时拥有结构自由度和声学自由度，组装时需要建立“结构节点编号—声学节点编号”的映射表，避免自由度错位。
对于位移-压力格式，组装完成后系统矩阵是非对称的。若希望使用对称求解器，可进行对称化处理。常见做法是引入变量替换，令 p = ρ c² q̇，把声压方程改写为以 q 为未知量的形式，使耦合项变为 S^T 与 S 的对称组合。另一种做法是直接对非对称矩阵做相似变换，但会改变特征向量的物理含义，后处理时需要还原。实际工程中，若声腔规模不大，直接使用非对称求解器往往更简单可靠；若声腔自由度占比很高，对称化带来的求解效率提升才值得付出实现代价。
综合来看，位移-压力格式与位移-位移格式并不是互相替代的关系，而是面向不同规模与求解需求的两种列式。理解耦合矩阵 S 与 K_c 的积分来源、界面匹配条件和对称化代价，才能在具体模型中做出合理选择，并为后续的模态分析与声压响应预测打下可靠基础。