3. 多孔吸声材料建模基础

多孔吸声材料的声学行为，本质上由“骨架 + 孔隙 + 空气”三者共同决定。声波进入材料后，空气在孔隙中往复运动，与孔壁摩擦生热，同时材料骨架在声压作用下发生微小变形，二者耦合形成宏观吸声。要在仿真中复现这一过程，必须把材料抽象成一组可计算的参数，并选择合适的等效流体模型。本节围绕 Delany-Bazley 模型、Johnson-Champoux-Allard（JCA）模型，以及流阻率、孔隙率、曲折度等关键参数的获取方法展开。

一、从现象到模型：两类等效流体思路

多孔材料内部结构极其复杂，直接建立真实几何几乎不可行。工程上普遍采用“等效流体”思路：把材料视为一种具有等效复密度 ρ(ω) 和等效复体积模量 K(ω) 的均匀介质，声波在其中传播时按这些等效参数衰减。区别在于，不同模型对 ρ(ω)、K(ω) 的表达式不同，所需输入参数也不同。

Delany-Bazley 模型是最早被广泛使用的经验模型。它只依赖一个参数——流阻率 σ，通过大量纤维材料实测数据拟合出幂律关系：

Zc = ρ0·c0·[1 + 0.0571·(ρ0·f/σ)^(-0.754) - j·0.087·(ρ0·f/σ)^(-0.732)]
k  = (ω/c0)·[1 + 0.0978·(ρ0·f/σ)^(-0.700) - j·0.189·(ρ0·f/σ)^(-0.595)]

其中 ρ0 为空气密度，c0 为声速，f 为频率，σ 为流阻率。该模型形式简单、参数少，适合纤维类材料（玻璃棉、岩棉、聚酯纤维）在中高频段的快速估算。但它有两个明显局限：一是低频段误差偏大，二是无法体现材料厚度方向上的结构差异，也无法描述孔隙率、曲折度对声速和衰减的独立影响。

Johnson-Champoux-Allard 模型则从物理机制出发，把黏性效应和热效应分开描述。它需要五个参数：流阻率 σ、孔隙率 φ、曲折度 α∞、黏性特征长度 Λ、热特征长度 Λ′。其复密度和复体积模量分别由 Johnson 等和 Champoux-Allard 给出的表达式计算，能够较准确地覆盖低频到高频的宽频范围，尤其适合对精度要求较高的仿真场景。

二、五个参数的物理含义与获取方法

流阻率 σ定义为单位厚度材料两侧压差与气流速度之比，单位 Pa·s/m²。它是 Delany-Bazley 模型唯一输入，也是 JCA 模型中最敏感的参数之一。获取方式主要有三类：一是直流气流法，用恒定低速气流穿过样品，测压差和流量，按 σ = Δp/(v·d) 计算，d 为厚度；二是声学反演法，用阻抗管测得的吸声系数或表面阻抗反推 σ；三是经验估算，对玻璃棉等纤维材料，σ 与容重、纤维直径存在近似关系，可作初值。教学示例中，若取 50 mm 厚玻璃棉、容重 24 kg/m³，σ 常在 8000～15000 Pa·s/m² 量级，具体值应以实测或厂家数据为准。

孔隙率 φ是材料内部连通空气体积与总体积之比，通常 0.85～0.99。测量可用压汞法、气体膨胀法或密度法：由骨架密度 ρs、材料表观密度 ρb 估算 φ = 1 − ρb/ρs。注意只有“连通孔隙”才对声波有效，闭孔部分应扣除。

曲折度 α∞描述孔隙通道的弯曲程度，直通道为 1，实际材料通常 1.2～3.0。它影响高频声波的惯性效应。获取较难，常用超声法测高频折射率反演，或用 JCA 模型对阻抗管数据做多参数拟合。工程初值可按纤维材料取 1.2～1.6。

黏性特征长度 Λ 与热特征长度 Λ′分别表征黏性耗散和热耗散发生的特征尺度，量级与孔径相当，通常几十到几百微米。二者可通过逆方法从阻抗管数据拟合得到，也可由孔径分布估算。Λ′ 一般略大于 Λ。

三、参数获取的实操路径与判断

实际建模时，推荐按“先易后难、先粗后精”的顺序推进。第一步，用密度法测 φ，用直流气流法测 σ，这两个参数成本低、重复性好。第二步，若只需中高频估算，直接代入 Delany-Bazley 模型，观察 500 Hz～4 kHz 的吸声系数曲线是否与阻抗管实测趋势一致。第三步，若低频偏差明显或需要宽频精度，切换到 JCA 模型，用阻抗管实测数据对 α∞、Λ、Λ′ 做最小二乘拟合。

下表给出一组教学示例参数，用于演示两类模型的输入差异，数值为合理假设，不代表实测结果：

把参数代入模型后，最终要落到可验证的观测量上：表面阻抗 Zs、吸声系数 α(f)、以及背后空腔深度 d 共同决定的共振峰位置。若仿真得到的 α(f) 峰值频率与阻抗管实测相差超过约 10%，应优先回查 σ 和厚度是否准确，再检查 α∞ 是否偏离合理范围。只有参数、模型、边界条件三者一致，多孔吸声材料的仿真结果才具备工程参考价值。