2. 声学波动方程与模态理论基础
声腔模态仿真的全部数学根基，可以压缩成一条演化链：从三维波动方程出发，经时谐假设化为 Helmholtz 方程，再叠加边界条件构成定解问题，最后通过分离变量得到一系列特征频率与特征函数，即模态。理解这条链，才能明白仿真软件里"求解模态"到底在解什么，也才能判断结果是否可信。
一、三维波动方程：声压随时间和空间如何演化
在均匀、静止、无黏、绝热的理想流体中，对连续性方程、动量方程和状态方程做线性化，消去质点振速，得到关于声压 p 的三维波动方程：
∂²p/∂t² = c² ∇²p
其中 c 为介质声速，∇² 为拉普拉斯算子
直角坐标下 ∇²p = ∂²p/∂x² + ∂²p/∂y² + ∂²p/∂z²它的物理含义很直接：某点声压对时间的二阶变化率，正比于该点声压分布的空间曲率。空间曲率越大（声压"凹凸"越剧烈），时间上振荡越快。声速 c 是唯一介质参数，空气在 20 ℃、1 atm 下约 343 m/s，水中约 1480 m/s。c 越大，同样的空间起伏对应更高的振荡频率。
波动方程是线性的，因此满足叠加原理：若 p₁、p₂ 都是解，则 a₁p₁ + a₂p₂ 也是解。这一性质是后面"模态叠加"能够成立的根本原因，也是仿真中可以把复杂激励分解成若干简单声源分别求解再相加的依据。
二、Helmholtz 方程：从时域走向频域
实际声场多由单频或窄带声源驱动，此时对时间做傅里叶变换，令 p(x,y,z,t) = P(x,y,z)·e^{jωt}，代入波动方程，时间导数 ∂²/∂t² 变为 −ω²，得到 Helmholtz 方程：
∇²P + k²P = 0
k = ω/c = 2πf/c 为波数（rad/m）
ω = 2πf 为角频率（rad/s）这一步把"时空耦合"的偏微分方程降为"纯空间"的方程，代价是引入频率参数 k。Helmholtz 方程的解 P 是复数，其模为声压幅值，辐角为相位。仿真中所谓"求解某频率下的声场"，就是在这个频率对应的 k 下解 Helmholtz 方程。
需要强调：Helmholtz 方程本身只描述介质内部的传播规律，它有无穷多解。要得到唯一确定的声场，必须补充边界条件——这正是下一节的核心。
三、边界条件分类：决定声场"长什么样"的关键
声腔的几何形状决定模态的"骨架"，边界条件决定模态的"血肉"。四类边界条件在数学形式与物理含义上各不相同：
| 类型 | 数学表达 | 物理含义 | 典型对象 |
|---|---|---|---|
| 刚性壁面 | ∂P/∂n = 0 | 法向振速为零，声波全反射 | 厚混凝土墙、刚性封闭腔 |
| 柔性/压力释放 | P = 0 | 声压恒为零，界面"软" | 自由液面、开口通大气 |
| 阻抗边界 | ∂P/∂n = −jωρP/Z | 部分吸收、部分反射 | 吸声衬里、多孔材料 |
| 辐射边界 | ∂P/∂n = −jkP（远场） | 声波单向向外传播不反射 | 开口辐射、无限域外边界 |
其中 n 为边界外法向，ρ 为介质密度，Z 为边界法向声阻抗。刚性边界对应 Z→∞，柔性边界对应 Z→0，阻抗边界是两者的连续过渡。辐射边界（Sommerfeld 条件）用于截断无限计算域，避免人工边界产生虚假反射。
边界条件的选择直接改变模态频率。以边长 1 m 的立方体空气腔为例：六面刚性时，第一阶模态频率约为 c/(2L) = 343/(2×1) ≈ 171.5 Hz；若一面改为压力释放，最低模态频率会降到约 c/(4L) ≈ 85.8 Hz。同一几何、不同边界，结果相差一倍，说明边界条件不是"细节"，而是定解问题的组成部分。
四、模态叠加原理与声腔模态的物理含义
对封闭声腔，在齐次边界条件下求解 Helmholtz 方程，得到特征值问题：
∇²Φₙ + kₙ²Φₙ = 0 （腔内）
边界条件（齐次）
kₙ = ωₙ/c 为第 n 阶特征波数
Φₙ 为第 n 阶模态振型（特征函数）特征值 kₙ 对应模态频率 fₙ = kₙc/(2π)，特征函数 Φₙ 描述该频率下腔内声压的空间分布形态。矩形刚性腔有解析解：
fₙ = (c/2)·√[(nx/Lx)² + (ny/Ly)² + (nz/Lz)²]
nx, ny, nz = 0,1,2,…（不全为零）
Lx, Ly, Lz 为腔体三边尺寸模态振型具有正交性：不同模态在腔体上的空间积分（加权）为零。正是正交性使得任意声场可以唯一展开为模态级数：
P(x,y,z,ω) = Σₙ aₙ(ω)·Φₙ(x,y,z)
aₙ(ω) = [∫∫∫ P₀·Φₙ dV] / (ωₙ² − ω² + jηωₙω)式中 P₀ 为声源分布，η 为阻尼损耗因子。分母揭示了共振本质：当激励频率 ω 接近某阶 ωₙ 时，该阶系数 aₙ 急剧增大，该模态在总声场中占主导。这就是"模态叠加"的工程含义——不必直接解复杂声场，只需算出各阶模态系数再求和。
声腔模态的物理含义可归纳为三点：其一，模态频率是腔体的固有属性，只由几何、边界和介质决定，与外部声源无关；其二，模态振型给出声压的"驻波图案"，节线处声压恒为零，腹点处声压最大；其三，模态密度随频率升高而增大，高频段模态重叠，单模态主导的假设失效，需考虑统计能量或宽带叠加。
五、从理论到仿真：参数与判断依据
把上述理论落到仿真操作上，关键输入与判断点如下。以某汽车乘员舱声腔为例（示例参数，非实测）：舱内空气 c = 343 m/s，ρ = 1.21 kg/m³，等效腔体尺寸约 2.4 m × 1.5 m × 1.2 m，壁面按刚性处理，座椅与内饰以阻抗边界近似（Z ≈ 800 Pa·s/m）。
按矩形腔公式估算前三阶：f₁₀₀ = 343/(2×2.4) ≈ 71.5 Hz，f₀₁₀ = 343/(2×1.5) ≈ 114.3 Hz，f₀₀₁ = 343/(2×1.2) ≈ 142.9 Hz。仿真求解后应重点核对：模态频率是否与解析估算同量级；振型中节线位置是否与腔体对称性一致；刚性壁面处声压梯度是否接近零。若某阶频率明显偏低，常见原因是边界被误设为柔性或阻抗过小；若出现大量非物理高频模态，则需检查网格是否过粗——通常要求每波长至少 6 个单元，即网格尺寸 h ≤ c/(6f_max)。
最后要区分两个容易混淆的概念：模态频率是腔体的固有属性，与声源无关；而某频率下的实际声压级分布，是声源通过模态系数 aₙ(ω) 加权后的结果。仿真中"扫频"看到的是后者，"模态分析"得到的是前者。只有先算准模态，后续的频响、声压级预测才有意义。