传递矩阵法：多孔材料特征阻抗与波数分析

这是一篇针对 Bryan H. Song 与 J. Stuart Bolton 发表于《美国声学学会会刊》（J. A. S. A.）的经典文献《基于传递矩阵法的柔软与刚性多孔材料特征阻抗及波数估算方法》（A transfer-matrix approach for estimating the characteristic impedance and wave numbers of limp and rigid porous materials）的深度学术解读。本文系统梳理了该文献的研究背景、核心理论框架、实验系统设计以及交叉验证结论，旨在为从事噪声控制、声学超材料设计以及声学有限元/边界元（FEM/BEM）仿真的科研人员提供严谨、详实的文献参考。

一、 研究背景与声学测试的科学瓶颈

在现代噪声控制工程中，玻璃纤维、聚合物纤维及各类泡沫等多孔材料被广泛应用于吸声与隔声设计。为了在设计阶段充分利用有限元、边界元或统计能量分析（SEA）等现代数值仿真工具对这些材料的声学性能进行预测，科研人员必须精确获取材料内部的宏观声学参数：特征阻抗（Characteristic impedance）和波数（Wave number）。

多孔材料的声学建模本质上是多相流体力学与固体力学的耦合问题。多数噪声控制材料属于多孔弹性体，理论上可支持两种纵波和一种横波的传播。然而，在特定的频率范围或材料属性下，许多均匀且各向同性的多孔材料可以被简化为“等效流体模型”：

  1. 刚性骨架（Rigid）：材料的固相密度或刚度极高，流体与固体的耦合极弱，固体骨架运动可忽略不计。
  2. 柔软骨架（Limp）：例如低密度、高流阻的非增强型玻璃纤维，其本体固体相的真空刚度相对于饱和流体（空气）可忽略不计。
    在这两种极限假设下，多孔材料内部仅支持单一的膨胀波（dilational wave）传播，因而可以完全由复数特征阻抗和复数波数这两个物理量进行全面表征。

早期测量方法的局限性：
历史上，声学界提出了多种提取这两种参数的测量手段，但均存在显著的工程局限。Scott提出的探针法会对多孔材料内部声场造成破坏性干扰；Ferrero 和 Sacerdote 提出的“双厚度法（two-thickness method）”以及 Yaniv提出的“双腔法（two-cavity method）”均依赖于对表面法向阻抗的测量。当测试材料具有高流阻和高耗散特性时，声波在材料内部衰减剧烈，导致表面阻抗对背板边界条件或材料厚度的变化极不敏感，从而引入巨大的计算误差。为了克服这一问题，学界开始引入穿透样品的传递函数测量（如 Champoux 和 Stinson 的方法），但这些方法通常需要预设特定的背腔深度或理想的刚性背板条件，测试流程依然繁琐。

二、 核心理论：传递矩阵法与对称互易性约束

本文的最大理论创新在于：巧妙地将流体管道元件（如消声器）测试中常用的传递矩阵法（Transfer matrix method）引入到多孔材料参数测试中，并结合均匀各向同性材料的对称性（Symmetry）与互易性（Reciprocity），成功将测量工作量减半，同时摆脱了对特定边界条件的依赖。

1. 四端网络与传递矩阵构建

研究人员将放置在阻抗管中的多孔材料层视为一个声学四端网络（two-port system），其上游（x=0）和下游（x=d）的声压（P）与法向质点速度（V）可以通过一个2×2的复数传递矩阵相联系：

其中，T11 表示下游处于零速度状态（刚性背板）时的上下游声压比，是一个无量纲量；T12 具有阻抗的量纲；T21 具有导纳的量纲。

使用MathJax渲染数学公式
\[ \left[ \begin{array}{l} P \ V \end{array} \right] _ {x = 0} = \left[ \begin{array}{l l} T _ {1 1} & T _ {1 2} \ T _ {2 1} & T _ {2 2} \end{array} \right] \left[ \begin{array}{l} P \ V \end{array} \right] _ {x = d} \]
2. 对称性与互易性的数学降维

在常规的“双负载法（Two-load method）”中，求解上述4个未知的矩阵元素（ T11 , T12 , T21 , T22 ）需要构建 4 个独立的方程，这意味着必须在阻抗管末端更换两种不同的物理边界（例如消声末端和刚性末端）进行两次完整的测试。
然而，对于均匀且各向同性的多孔材料层，本文指出其满足两大物理约束：

  • 对称性约束：材料层从两个方向看声学特性相同，因此 T11 = T22 。
  • 互易性约束：对于线性无源互易声学网络，传递矩阵的行列式恒为1，即 T11 T22 – T12 T21 = 1。

通过引入这两个约束方程，4个未知数被缩减至仅需两个额外的方程即可求解。这意味着，测试者只需在一种任意的、甚至未知的下游末端边界条件下，进行一次单次实验，即可完全解析出多孔材料的传递矩阵。

3. 材料宏观参数的解析反演

一旦利用上下游的声压及速度（由阻抗管中的四位置传声器通过平面波分离技术求得）计算出矩阵元素，多孔材料的核心属性便可直接以闭式解析解的形式提取：

使用MathJax渲染数学公式
\[ 复数波数（k_p）：由方程 k_p = \frac{1}{d} \cos^{-1}(T_{11}) 或 k_p = \frac{1}{d} \sin^{-1}(\sqrt{-T_{12}T_{21}}) 求得。 \]
使用MathJax渲染数学公式
\[ 特征阻抗（Z_c）：由方程 Z_c = \sqrt{T_{12}/T_{21}} 提取。 \\ 进而，这些参数可以进一步转化为有限元分析所需的复数声速（c_p）和复数等效密度（\rho_p）。 \]

三、 实验系统设计与系统误差控制

为了验证上述理论，作者基于商业化的 Brüel & Kjær Type 4206 双传声器阻抗管进行了系统性硬件改造。

  1. 四位置麦克风阵列：在传统阻抗管的样品管下游，增加了一段与上游等径的测量管，布置了第3和第4个麦克风测量点。上游测量点用于提取入射波与反射波（系数 A 和 B ），下游测量点用于提取透射波（系数 C 和 D ）。
  2. 末端配置：在系统末端采用3M Thinsulate吸声材料构造了一个近似的消声末端（Anechoic termination）。虽然理论上本方法不要求已知末端阻抗，但使用消声末端可使下游声场接近纯行波状态，大幅降低传递函数估算中的驻波畸变误差。
  3. 克服麦克风相位失配：作者在实验中采用了单探头麦克风扫描（Roving microphone）的策略。由于所有的传递函数均参考同一个扬声器输入电信号，使用同一只麦克风在四个位置依次测量，从物理根源上消除了多麦克风系统极其敏感的相位标定误差（Phase-calibration error）。

实验选用的材料为航空级玻璃纤维（A组密度6.7 kg/m³，B组密度9.6 kg/m³），并将材料切成2.5 cm厚的圆柱块组合测试。

Schematic of the measurement setup.

四、 实验结果与多维交叉验证

本文最核心的学术价值在于其通过极其严密的对照实验，证明了该测试方法的高鲁棒性。

1. 与半经验公式的高度吻合

实验测得的波数实部/虚部（表征相位与衰减）、特征阻抗、复数相速度（归一化）以及复数密度，均与多孔纤维材料界经典的 Delany-Bazley 半经验理论公式预测值高度吻合（特别是在 1 kHz以上的中高频段）。

2. 独立于材料厚度的验证

作为一个物理内秉属性，特征阻抗和波数不应随材料的测试厚度而改变。作者分别测试了 2.5 cm、5.0 cm 和 7.5 cm 厚度的玻璃纤维组合。结果表明，5.0 cm与 7.5 cm样本解析出的波数和特征阻抗曲线几乎完全重合。仅在 2.5 cm 极薄样本中出现轻微偏差，作者严谨地指出这可能是由于轻薄样本在管内边缘密封不严导致的声学泄漏（edge leakage）或层间不均匀性引起的。

3. 独立于末端边界条件的验证

这是证明该方法优越性的关键实验。作者分别将下游阻抗管末端设置为：(1) 消声末端；(2) 刚性刚性背板（全反射）；(3) 自由开口边界。在这三种声场反射特性截然不同的边界下，所提取出的特征阻抗、波数以及材料自身的传递损失（Transmission Loss, TL）曲线几乎完全一致。这确凿地证明了：该传递矩阵解耦算法彻底剥离了测试系统的边界串扰。

4. 效率对比：与“双负载法”的比较

作者将本方法提取的结果与传统的“双负载法（Two-load method）”进行了背靠背对比。结果显示，两者的曲线高度重合。但由于传统双负载法不利用互易性和对称性，必须在消声末端和刚性末端下分别测试两次，而本方法不仅精度相当（甚至在高频段由于下游声场无极端驻波而表现出更少的数值振荡），且测试时间与工作量直接减半。

5. 方法的物理极限与误差溯源

文章也客观地探讨了该方法在特定频段的测量畸变及其物理根源：

  • 低频偏差（< 1000 Hz）：在低频区，特征阻抗及传递损失的测量值偏离理想无限大平板理论。这是经典的“边缘约束效应（Sample-edge constraint）”。低频下固相与流相间存在强烈的粘性耦合，样品在阻抗管壁的摩擦约束阻碍了流体膨胀，导致有效声学参数发生改变。
  • 高频突变（~4.6 – 4.9 kHz）：在高频段部分数据出现奇异峰值。作者进行了结构模态分析，发现夹具套筒在该频段存在管壁的“椭圆化模态（Ovalling mode）”。这种管壁的机械共振强烈耦合了内部声场的高阶周向模态，导致了一维平面波假设的局部失效。

五、 研究亮点与主要结论

综上所述，Song 与 Bolton 的这项研究在声学材料测试领域具有里程碑式的改进意义。其主要研究亮点和结论可归纳如下：

  1. 理论降维与效率翻倍：打破了必须进行两次独立阻抗管测试的传统认知。通过引入均匀各向同性材料的“对称性（ T11=T22 ）”和“互易性（ det[T]=1 ）”宏观约束，将求解传递矩阵的自由度减半，实现了一次测试提取全部复数声学参数的创举。
  2. 极高的系统鲁棒性：该方法不对下游系统的边界条件做任何要求，不仅极大地简化了实验装置的设计（无需精密的刚性背板或机械调腔），更有效避免了强耗散材料表面阻抗不敏感带来的计算发散问题。
  3. 为声学仿真提供完备输入：该算法不但能计算传统的吸声系数与反射系数，还能直接输出复数等效密度（ ρp ）和复数声速（ cp ）。这两项指标恰恰是现代商业声学软件（如 SYSNOISE、Actran 等）建模多孔介质等效流体域时最为依赖的核心数据。
  4. 适用范围限定提示：作者明确指出，本方法完美适用于能够等效为流体的刚性骨架或极软骨架材料（如玻璃棉）。然而，若应用于支持显著横波模态的某些聚氨酯泡沫塑料（复杂多孔弹性体）时，由于传递矩阵的结构远比二维流体矩阵复杂，需谨慎评估其等效适用性。

整体而言，这篇文章兼具深刻的物理洞察力与极高的工程实用价值。其确立的“四麦克风传递矩阵法（Four-microphone transfer-matrix method）”现已成为评估隔声垫、吸音棉及汽车座舱声学包（Acoustic package）等被动噪声控制材料的行业标准测试范式之一。在声学材料设计与仿真分析中，特征阻抗与波数的准确获取直接决定了模型预测的可靠性。传递矩阵法为多孔材料声学参数的高效估算提供了重要理论基础，但在实际应用中，对测试系统精度与数据稳定性的要求同样严格。针对这一需求，苏州东原电子推出的声学材料测试系统与阻抗管解决方案，基于成熟的测量方法与稳定的数据采集能力，可实现多孔材料关键声学参数的精准获取，广泛服务于声学仿真、材料研发及工程应用场景，为声学设计提供可靠的数据支撑。

B. H. Song and J. S. Bolton, “A transfer-matrix approach for estimating the characteristic impedance and wave numbers of limp and rigid porous materials,” J. Acoust. Soc. Am., 2000

admin

类似文章

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注