4. 声学材料属性与边界条件设置

在声腔模态与声场分布仿真中，几何模型只决定了“声波能在哪里传播”，而真正决定声场形态、模态频率和衰减特性的，是介质属性与边界条件。同一腔体，空气温度从 20 ℃ 变为 200 ℃，或把刚性壁换成 30% 穿孔率的吸声衬里，模态频率与声压分布都会明显改变。本节围绕三件事展开：流体介质参数如何取值、阻抗边界如何表达、吸声材料与多层结构如何等效建模。

一、空气与流体介质参数

对理想气体中的小振幅声波，声速与介质密度由温度决定：

c = sqrt(γ·R·T)
ρ = p_abs / (R·T)

γ = 1.4（空气比热比）
R = 287.05 J/(kg·K)（空气气体常数）
T 为热力学温度，单位 K

以 20 ℃（293.15 K）为例，c ≈ 343.2 m/s，ρ ≈ 1.204 kg/m³。若腔内温度升到 200 ℃（473.15 K），c ≈ 436 m/s，密度降到约 0.746 kg/m³。声速变化会直接平移所有模态频率，因此仿真前必须明确工作温度，而不是默认 343 m/s。对于非空气介质，如烟气、氮气或水，应分别给出 γ、R 或直接给出实测声速与密度；水中声速约 1480 m/s，密度约 998 kg/m³，与空气相差三个数量级，阻抗边界的效果也完全不同。

除声速和密度外，还需关注体粘滞系数与热传导损耗。在窄缝、小孔或高频段，边界层损耗不可忽略，此时可在材料或流体属性中引入复声速或复波数，把损耗折算为等效衰减。

二、阻抗边界条件

声学边界通常用比声阻抗 Z_s 描述，定义为界面声压与法向质点速度之比。工程上更常用归一化阻抗：

z = Zs / (ρ·c) = r + j·x

r：归一化阻性分量（耗散）
x：归一化抗性分量（惯性/弹性储能）
刚性壁：|z| → ∞，法向速度为零
压力释放边界：z = 0，声压为零
理想匹配层：z = 1 + j0，无反射

在有限元或边界元中，阻抗边界一般以法向速度或声压梯度形式施加。需要注意符号与法向方向约定：若法向指向流体外部，则阻抗关系写作 p = −Z_s·v_n；若法向指向流体内部，符号相反。全文与模型必须统一，否则会出现“吸声边界反而放大反射”的错误结果。实际设置时，刚性壁可直接用零法向速度，吸声衬里则输入 r、x 随频率变化的曲线或表格。

三、吸声材料建模：Delany-Bazley 与 Miki 模型

多孔吸声材料（玻璃棉、岩棉、聚酯纤维）的声学行为主要由流阻率 σ（Pa·s/m²）和层厚决定。Delany-Bazley 模型基于大量纤维材料实测数据拟合，用无量纲量 X = ρf/σ 给出特性阻抗与传播常数：

Delany-Bazley（经验拟合）：
Zc/(ρc) = 1 + 0.0571·X^(-0.754) − j·0.087·X^(-0.732)
k/k0   = 1 + 0.0978·X^(-0.700) − j·0.189·X^(-0.595)
X = ρ·f / σ

该模型在 X 较大（中高频、流阻率适中）时较准，但在低频或 X 较小时可能给出负实部，导致非物理结果。Miki 模型对 Delany-Bazley 做了修正，保证实部为正，更适合宽频与低频段：

Miki 模型：
Zc/(ρc) = 1 + 0.0699·X^(-0.632) − j·0.107·X^(-0.632)
k/k0   = 1 + 0.109·X^(-0.618) − j·0.160·X^(-0.618)

例如取 σ = 20000 Pa·s/m²、层厚 50 mm、频率 500 Hz，则 X = 1.204×500/20000 ≈ 0.0301。代入 Miki 公式可算得归一化特性阻抗与传播常数，再结合刚性背衬的传递矩阵，即可得到该层的表面阻抗与吸声系数。实际建模时，可把多孔层等效为具有复声速和复密度的流体域，也可直接折算成阻抗边界施加在壁面上，后者计算量更小，但只适用于层内无显著横向波传播的情形。

四、穿孔板与多层结构的等效参数

穿孔板由小孔和板后空腔组成，常用 Maa 模型等效。单孔声阻抗包含短管惯性项与孔口黏性、热损耗项：

Z_hole = j·ω·ρ·t_eff + 损耗项
t_eff = t + 0.85·d（d 为孔径，t 为板厚）
穿孔率 p = 孔面积 / 板面积
板整体归一化阻抗 z = Z_hole / (p·ρ·c)

孔径越小、板越厚、穿孔率越低，阻性损耗越强，吸声峰越宽但峰值频率向低频移动。多层结构（穿孔板 + 多孔层 + 空腔）可用传递矩阵法逐层串联：每层用其特性阻抗和传播常数构造 2×2 矩阵，最终由总矩阵求出表面阻抗，再换算吸声系数。这样既保留了各层物理参数，又避免了对每根纤维或每个孔做几何建模。

五、建模设置要点与结果判读

把上述参数落到仿真中，建议按以下顺序检查。第一，确认温度与介质，避免用默认 343 m/s 覆盖实际工况；第二，明确每个边界的类型与法向方向，刚性壁、阻抗面、对称面、压力释放面分别对应不同的数学条件；第三，吸声层优先用 Miki 模型，若已有实测吸声系数曲线，可直接以阻抗边界形式输入，减少经验模型误差；第四，穿孔板与多层结构用传递矩阵先算出表面阻抗，再施加到壁面，比逐孔建模高效得多。

结果判读时，重点看三处：模态频率是否随温度、声速变化而整体平移；阻抗边界的实部 r 是否为正、是否在目标频段内产生合理吸收；吸声层厚度与流阻率变化时，吸声峰是否按预期移动。若出现低频吸声系数异常偏高或模态频率不随介质变化，通常说明边界条件施加位置或法向约定有误，应回到模型检查边界定义与材料参数的单位一致性。