1. 声学仿真工程概览
声学仿真并不是“把声压算出来”这么简单。工程中真正要回答的问题通常是：某台设备在特定转速下，辐射噪声在哪个频段超标？一间会议室的语言清晰度是否足够？水下声呐罩在特定频率下是否形成透声窗口？这些问题分别落在振动噪声（NVH）、室内声学和水下声学三类典型场景中，而有限元法（FEM）与边界元法（BEM）正是支撑这些场景的两大数值工具。理解它们的角色分工，比记住公式更重要。
一、两类方法的核心差异：域离散 vs 边界离散
有限元法把整个声学域（空气腔、水腔、结构本体）剖分为体单元，在单元内构造形函数，最终求解稀疏、带状的大型线性方程组。它的优势是能同时处理结构—声耦合：结构振动作为声源的边界条件，声压反作用于结构形成附加质量与阻尼。因此，FEM 特别适合封闭腔体内部声场、结构模态与声腔模态耦合、以及中低频段（通常 20 Hz～1 kHz）的振动噪声问题。
边界元法只离散辐射面或散射面，把三维域问题降为二维边界问题，天然满足无穷远 Sommerfeld 辐射条件，不需要人为截断外场计算域。因此 BEM 更适合外声场辐射与散射：发动机辐射噪声、水下目标声散射、开口或半开口结构的远场指向性。代价是系数矩阵为满阵，频率升高时求解规模增长快，通常用于中频段（几百 Hz～数 kHz）的外场问题。
二、三类场景中的角色分工
振动噪声（NVH）是最典型的“结构—声耦合”问题。以汽车车内噪声为例：发动机悬置点输入激励，经车身板件振动向车内辐射。工程上通常先用 FEM 计算白车身模态与声腔模态，检查是否有模态耦合风险；再用 FEM 或 FEM-BEM 混合方法计算 20～500 Hz 的驾驶员耳旁声压级。若关注车外通过噪声，则把车身表面振速作为 BEM 的边界条件，计算远场声压。
室内声学在低频段（房间尺寸与波长可比时）必须用 FEM 或 BEM 求解波动方程，才能捕捉驻波、简并模态和角落声压聚集；高频段则转向几何声学（射线追踪、声线法）或统计能量分析（SEA）。一个音乐厅的座椅吸声、墙面扩散体、舞台耦合空间，在 125 Hz 以下用 FEM 建模是合理选择。
水下声学的介质是水，声速约 1500 m/s，同频率下波长比空气长 4 倍以上，因此相同网格尺寸可覆盖更低频段。潜艇声呐罩透声、鱼雷自噪声、海底管线泄漏声辐射，常用 BEM 或 FEM-BEM 混合方法。水的重质负载使结构—流体耦合更强，FEM 中的附加质量效应不可忽略。
三、典型工作流：从几何到结论
无论哪类场景，声学仿真都遵循一条可复现的主线：几何清理 → 网格划分 → 材料与边界条件 → 求解 → 后处理判读。关键差异在于边界条件的物理含义。
以车内噪声 FEM 分析为例，一个可执行的教学示例参数如下：声腔空气密度 1.225 kg/m³，声速 343 m/s；结构板件钢材料弹性模量 210 GPa，密度 7850 kg/m³，阻尼损耗因子 0.02；声腔体网格尺寸按最高分析频率 500 Hz 控制，每波长至少 6 个单元，即单元尺寸 ≤ 343/(500×6) ≈ 0.114 m，实际取 0.08 m 留裕量。激励为悬置点单位力 1 N，频率扫描 20～500 Hz，步长 2 Hz。
求解后重点观察三条曲线：驾驶员右耳声压级频谱、板件表面法向振速云图、声腔模态频率表。若某阶声腔模态与板件模态频率差小于 5%，且该频率处声压级出现峰值，则判定存在模态耦合，需通过加强筋、阻尼层或改变声腔容积来错频。
四、软件生态与工具链
工程中很少只用一款软件。常见组合是：前处理用 ANSA、HyperMesh 或 Simcenter 3D 完成几何清理与网格划分；求解用 Actran、Virtual.Lab、COMSOL、LMS Samtech 或自编 MATLAB/Python 脚本；后处理用 ParaView、Tecplot 或软件自带模块。BEM 专用工具包括 LMS Virtual.Lab Acoustics、Actran、WAON 等。
选择工具时优先看三件事：是否支持结构—声耦合、是否支持所需频率范围的求解器、是否能与现有 CAD/CAE 流程对接。不要因为某软件“功能全”就忽略网格规模与计算资源——一个 500 Hz 的车内声腔模型，体网格量级常在百万级，直接求解内存需求可能超过 64 GB。
五、常见误区与排障方向
初学者最容易犯的错误是“网格越细越好”。声学网格只需满足每波长 6～10 个单元，过密会显著增加求解时间而不提升精度。第二个误区是忽略边界条件的物理一致性：BEM 中若把刚性壁面误设为声软边界，远场指向性会完全错误。第三个误区是频率范围与单元尺寸不匹配，导致高频结果出现数值色散。
排障时按顺序检查：几何是否有微小缝隙或重叠面；网格是否有自由边、负体积；材料参数单位是否统一（特别注意 Hz 与 rad/s、Pa 与 MPa）；边界条件是否覆盖所有外表面；求解频率步长是否足够捕捉峰值。若结果异常，先做 1 个单元的单频验证，再逐步放大模型。