3. 数值模拟方法总览

振动与声波问题在数学上最终都归结为对波动方程（或弹性动力学方程）的求解。解析解只对极少数规则几何与简单边界条件存在，工程中的实际对象——梁、板、壳体、管道、消声器、房间、土壤——几乎都必须借助数值方法。本节不推导公式，而是从“离散思想—网格与未知量—边界处理—适用场景”四个维度，把有限差分法（FDM）、有限元法（FEM）、边界元法（BEM）和谱方法（Spectral Method）放在同一张“方法地图”上对比，帮助你在具体问题面前做出合理选择。

四种方法的基本思想

有限差分法从微分方程出发，用差商代替导数。以一维波动方程 ∂²u/∂t² = c²∂²u/∂x² 为例，把空间离散为步长 Δx 的节点，二阶空间导数用中心差分 (u_i+1−2u_i+u_i−1)/Δx² 近似，时间同样离散后得到显式递推格式。它的核心优势是实现直接、计算量小、易于并行，对均匀介质中的波传播非常自然；缺点是规则网格难以贴合复杂曲面边界，阶梯状近似会引入虚假散射，且显式格式受 CFL 稳定性条件 cΔt/Δx ≤ 1 约束。

有限元法从等效积分弱形式出发，把求解域剖分为单元，在每个单元内用形函数插值未知场，组装成整体刚度、质量矩阵后求解 K u + M ü = f。它的最大优势是几何适应性强，三角形、四面体、壳单元可以贴合任意复杂结构，边界条件（自由、固定、阻抗）处理规范，因此成为结构振动与声固耦合分析的主流工具。代价是需要对整个体积域划分网格，高频问题网格量急剧膨胀，且存在数值色散——网格过粗时波速被高估。

边界元法基于边界积分方程，利用格林函数把域内问题转化为边界上的积分方程，只需对边界离散。对无限域外声场问题（如管道辐射、水下声辐射）它天然满足 Sommerfeld 辐射条件，无需人工截断边界，这是它相对 FEM 的独特优势。但 BEM 生成的是满阵，规模增大时存储与求解成本按 O(N²)～O(N³) 增长，且对多域、非均匀介质处理不便，特征频率处还存在非唯一性问题。

谱方法用全局光滑基函数（傅里叶级数、Chebyshev 多项式）展开解，而非局部的分段插值。对光滑解，其误差随基函数个数呈指数收敛，即“谱精度”，在同样精度下所需自由度远少于低阶 FEM/FDM。它特别适合规则几何、均匀介质、光滑解的波传播与稳定性问题，例如均匀管道声场、简支梁高阶模态。缺点是几何适应性差，解中出现激波或间断时会出现 Gibbs 振荡，边界处理也不如 FEM 灵活。

离散对象与未知量：理解差异的关键

四种方法最本质的区别在于“离散什么”。FDM 离散的是微分算子，未知量是网格节点上的场值；FEM 离散的是弱形式积分，未知量是单元节点自由度，物理意义更接近“广义位移”；BEM 离散的是边界积分方程，未知量是边界上的声压与法向振速；谱方法离散的是解的全局展开系数，未知量是各阶模态或傅里叶分量的幅值。这一差异直接决定了网格规模、矩阵结构和边界处理方式。

以矩形房间内声场为例：若用 FDM，取 0.05 m 网格、房间 5 m×4 m×3 m，约需 60×48×36 ≈ 10 万节点；用 FEM 四面体网格量级相近但可贴合吸声壁面；用 BEM 只需离散约 100 m² 的壁面，单元数可能仅数千；若房间是规则矩形且壁面均匀，谱方法用几十个模态即可获得很高精度。这个数量级对比说明：方法选择首先取决于几何规则性与域的类型，而不是精度口号。

适用场景对比与选型建议

实操层面的几个判断要点

第一，先看域是有限还是无限。外场辐射、散射问题优先考虑 BEM；若必须用 FEM，则要在截断边界上施加吸收边界条件（如 PML 或一阶 Sommerfeld 条件），并验证截断距离是否足够远。第二，看几何能否用规则网格描述。若结构含大量曲面、倒角、加强筋，FDM 的阶梯近似会带来明显虚假反射，此时 FEM 更稳妥。第三，看分析频段。以 FEM 为例，经验上每个波长至少需要 6～10 个单元才能控制色散误差；若最高分析频率为 2 kHz、空气中波长约 0.17 m，单元尺寸应取 0.02 m 量级，网格量会迅速上升，这时可考虑谱方法或 BEM。第四，看是否需要时域瞬态。FDM 显式格式在时域波传播中实现简单、并行效率高；FEM 时域分析需配合 Newmark 等时间积分，BEM 时域实现则相对复杂。

最后需要强调，任何数值方法的结果都必须经过验证：网格收敛性检查（加密网格看结果是否稳定）、与解析解或实验数据对比、检查能量守恒或互易性。方法只是工具，对物理问题的判断力才是数值模拟真正的核心。