1. 课程导入与声学结构耦合概述
声学结构耦合仿真，研究的不是“声音在空气里怎么传播”这一单独问题，而是流体介质中的声压场与固体结构振动之间的双向作用。当结构振动时，它会推动相邻流体产生声压；而流体中的声压又会反过来在结构表面施加压力，改变结构的振动状态。这种“你推我、我推你”的相互作用，就是声固耦合（Acoustic-Structure Interaction）的物理本质。在COMSOL Multiphysics中，这一本质被表达为声学域与固体力学域在耦合界面上的法向位移连续与法向应力平衡两个条件。
理解这一点非常关键：如果只做单向计算，比如先算结构振动再把它当作声源，那只是“声辐射”的近似；真正的耦合仿真必须让声压和位移在同一个求解步中相互反馈。COMSOL的“声-结构相互作用”多物理场接口正是为此设计的，它自动在声学域与固体域的交界面上建立耦合关系，无需手动逐点定义。
耦合界面的两个核心条件：
- 法向位移连续：结构表面的法向位移等于流体质点的法向位移，保证介质不分离、不穿透。
- 法向应力平衡：结构表面受到的法向应力等于声压（方向相反），即 σ·n = −p·n。
典型工程场景中，声固耦合几乎无处不在。水下声呐换能器工作时，压电陶瓷的振动通过水介质向外辐射声波，同时水的反作用力会改变换能器的谐振频率；汽车NVH分析中，车身板件振动向车内空腔辐射噪声，而车内声压又会让板件产生附加振动，形成“轰鸣”问题；扬声器振膜推动空气发声，空气负载会降低振膜的共振频率；超声换能器中，压电片、匹配层与负载介质之间的声阻抗匹配直接决定能量传输效率。这些场景的共同点是：流体负载不可忽略，且结构与流体之间存在明显的能量交换。
在COMSOL中建立声固耦合模型，通常遵循一条清晰的流程。首先是几何与材料：分别建立流体域和固体域，并指定各自的材料属性——流体需要密度和声速，固体需要弹性模量、泊松比和密度。对于压电材料，还需输入压电耦合矩阵和介电常数。其次是物理场选择：流体域添加“压力声学，频域”或“压力声学，瞬态”，固体域添加“固体力学”，然后在多物理场节点下添加“声-结构相互作用”耦合。这一步是自动完成的，COMSOL会识别所有声学域与固体域的公共边界，并在其上施加耦合条件。
接下来是耦合界面设置。虽然耦合是自动建立的，但用户仍需检查界面是否被正确识别。例如，如果流体域和固体域之间存在薄壳或刚性连接，可能需要手动指定耦合边界。对于扬声器这类带有振膜、音圈和磁路的复杂结构，通常还需要考虑预紧力或非线性大变形的影响，此时耦合条件需要与几何非线性选项配合使用。
网格与求解是决定计算成败的关键。声学域网格密度通常要求每个波长至少包含6～10个单元，而固体域在弯曲波传播时也需要类似的分辨率。对于频域分析，COMSOL使用直接求解器或迭代求解器处理耦合矩阵；对于瞬态分析，则采用隐式时间步进。一个常见的教学示例是：一块铝板（E=70 GPa，ν=0.33，ρ=2700 kg/m³）浸没在水中（ρ=1000 kg/m³，c=1500 m/s），板厚5 mm，尺寸300 mm×300 mm，在中心施加1 N的简谐点力，频率范围100～2000 Hz。通过耦合仿真可以观察到，水的附加质量效应使板的基频从空气中的约180 Hz下降到水中的约95 Hz，同时辐射声压峰值出现在板的一阶弯曲模态附近。这些数值属于合理假设的示例，用于说明耦合效应，并非实测结果。
学习建议：本章后续将逐步展开从几何建模到结果判读的完整操作。建议先理解耦合界面的两个条件，再动手搭建第一个模型。不要急于追求复杂几何，先用一块平板和一个水槽验证耦合是否生效——观察水中基频是否明显低于空气中基频，这是最直观的验证方法。
后处理阶段，除了查看声压级分布、结构位移云图外，还应关注耦合界面的能量传递。COMSOL允许计算声强、辐射功率和结构表面法向应力，这些量能帮助判断耦合是否充分、是否存在数值反射或虚假模态。对于NVH问题，通常还需要提取场点声压级曲线，与实验或经验数据对比；对于换能器问题，则关注电-声转换效率和指向性图。