变参数分空腔单层微穿孔板低频降噪性能研究*

杨菲，沈新民，王强，张蕉蕉

（陆军工程大学野战工程学院 南京，210007）

引言

随着声学侦测技术的高速发展，军用装备的战场生存能力受到威胁［1］。低频噪声因其具有较强的绕射性和透射性，容易穿越障碍物，在传播过程中随着距离的增加不会明显衰减，增加了被敌方武器侦测和袭击的危险［2］。因此，弱化军用装备的低频噪声，可以降低被敌方发现和打击的可能性。由于军事装备种类和型号的多样性，存在的低频噪声频带带宽也不尽相同，为节约制作成本，选取常见的100～500 Hz，100～1 000 Hz，100～1 500 Hz和100～2 000 Hz进行研究，针对不同的带宽设计相匹配的吸声器，为军用装备低频降噪控制技术的实际应用提供理论指导。

微穿孔板（microperforated panel，简称MPP）吸声结构是一种防水防潮、表面强度高及无二次污染的共振吸声结构，可以胜任极端的工况环境，能够很好地适应复杂多变的战场环境［3］。微穿孔板吸声结构的吸声性能主要取决于自身的4个结构参数（孔径d、板厚t、孔间距b、空腔深度D）。在前人的研究中，大多数是理论上的分析而未将现实约束条件和目标频率范围考虑在内。笔者所研究的变参数分空腔单层微穿孔板结构是对4个结构参数进行设计与优化，以获得最佳的吸声效果。

通过优化算法搜索得到最优的结构参数组合可以使变参数分空腔单层微穿孔板结构的吸声性能最大化。布谷鸟搜索算法具有利用参数少、操作简单、易于实现、优良的随机搜索路径以及出色的优化能力等优点，已被用于优化声学结构参数［4-5］。Duan等［4］通过布谷鸟搜索算法优化了MPP和多孔金属复合结构的参数，在2 000～6 000 Hz的频率范围内获得的平均声音吸收系数为0.976 5。Yang等［5］通过布谷鸟搜索算法优化了具有有限尺寸的标准化多层MPP的结构参数，层数为1～4的标准化多层MPP在频率范围100～6 000 Hz内平均吸声系数分别为0.574 5，0.708 5，0.719 9和0.722 8。所以，笔者选取布谷鸟搜索算法作为工具，对变参数分空腔单层微穿孔板的结构参数进行优化。

首先，基于微穿孔板理论［6］构建了变参数分空腔单层微穿孔板结构的理论吸声模型，分析结构参数与吸声系数的关系；其次，确定优化目标和约束条件，通过布谷鸟搜索算法对变参数分空腔单层微穿孔板的结构参数进行优化，为后续样品的制备和吸声性能的评估提供了有效指导；然后，利用有限元仿真初步评估了最佳变参数分空腔单层微穿孔板结构的吸声性能，并通过驻波管测试验证有限元仿真的准确性和有效性；最后，对于每个最佳的变参数分空腔单层微穿孔板结构吸声性能进行了比较，为不同频率范围找到最佳的吸声器，这为装备吸声降噪提供了理论指导，具有重要的军事意义和应用价值。

1 变参数分空腔单层微穿孔板结构

微穿孔板的结构参数包括板的厚度、微孔的直径、相邻孔间距离以及空腔长度，这些参数影响其吸声性能。笔者研究的变参数分空腔单层微穿孔板结构是由多组结构参数不相同的微穿孔板并联拼接而成，其基本结构如图1所示（以4组结构参数微穿孔板为例）。微孔为圆形、呈线性等距排列，并且不同微孔的板后空腔需要用隔板隔开，使各不同孔径微孔独立的吸声特性有效地结合起来［7］。焦风雷等［8］已经证明，不同孔径组成的单层变孔径微穿孔板结构可以拓宽吸声频带。

图1 变参数分空腔单层微穿孔板结构示意图Fig.1 Schematic diagram of a single-layer microperforated panel structure with variable parameters and separated cavities

笔者所研究的变参数分空腔单层微穿孔板结构中不同结构参数的组数为1～4，优化目标是在100～500 Hz，100～1 000 Hz，100～1 500 Hz和100～2 000 Hz的频率范围内获得最佳的吸声性能。变参数分空腔单层微穿孔板结构的理论模型可以通过声电类比法构建。在声学研究领域，很多研究已经通过简化的电路图对吸声结构的吸声特性进行分析，此方法被称为声电类比法。因此，图1的结构可以理解为微穿孔板上不同结构参数的部分微穿孔板的声阻抗与各自的板后空腔串联再相互并联，由此可以得出等效电路如图2所示。

图2 变参数分空腔单层微穿孔板结构等效电路图Fig.2 The equivalent circuit diagram of a single-layer microperforated panel structure with variable parameters and separated cavities

对于笔者研究的变参数分空腔单层微穿孔板结构，其第i组结构参数的声阻抗率Zsi可以根据微穿孔板理论［6］通过式（1）计算得到

其中：Ri为微穿孔板的声阻；ω为声波角频率；Mi为微穿孔板的声质量。

Ri和Mi可以分别由式（2）和（3）计算得到

其中：ρ为空气密度；μ为空气的黏度系数；υ为金属板的导热系数；εi为穿孔率；ti和di分别为微穿孔板的厚度和微孔的直径；ki为穿孔板常数。

其中：di和bi分别为微孔的直径和相邻微孔之间的距离；f为声波频率。

根据微穿孔板理论［6］，可以推导出变参数分空腔单层微穿孔板结构的总声阻抗率ztotal为

其中：Di为板后空腔长度。

由此可得变参数分空腔单层微穿孔板结构的吸声系数为

2 布谷鸟搜索优化算法

根据第1节中理论吸声模型，选取布谷鸟搜索优化算法，优化获得变参数分空腔单层微穿孔板结构的最佳结构参数。

2.1 优化流程

对于变参数分空腔单层微穿孔板结构，其吸声性能由各组微穿孔板的结构参数共同决定。基于布谷鸟搜索算法的变参数分空腔单层微穿孔板结构的优化设计步骤如下。

1）建立优化模型。以变参数分空腔单层微穿孔板结构在给定的频率范围[fmin，fmax]平均吸声系数最大为目标，建立目标函数，即

2）确定决策变量和约束条件。考虑实际加工工艺条件限制，对各结构参数设定如下约束条件（数值单位为m）

3）确定布谷鸟搜索算法的运算参数。宿主鸟窝种群数量n=20，最大发现概率pa=0.25，步长控制因子的最大值βmax≤1，最小值βmin≥0，步长控制因子通过式（11）计算得到。设置种群进化最大迭代次数为10 000，生成初始的随机种群，并令迭代次数初始值为1。随机生成符合条件的自变量，计算每一个宿主鸟窝的目标函数值并记录，不断迭代直至满足结束条件，输出最优解。

其中：z为迭代的次数；a为β的衰减率的初始值。

2.2 优化结果

布谷鸟搜索优化得到的最佳结构参数以及相应的平均吸声系数理论值分别如表1～4所示。优化结果表明，对于每个目标频率范围，变参数分空腔单层微穿孔板结构的平均吸声系数随着结构参数组数从1～4的增加而增长。从变参数分空腔单层微穿孔板结构的最佳平均吸声系数可以看出，通过优化结构参数可以改善其吸声性能。

表1 目标频率范围为100～500 Hz时的最佳结构参数Tab.1 The optimal structural parameters for the target frequency range of 100～500 Hz

尽管通过布谷鸟搜索算法获得了目标频率范围的变参数分空腔单层微穿孔板结构的理论最佳结构参数，但仍需通过有限元仿真和驻波管测试，检验理论模型与布谷鸟搜索算法的有效性和可靠性。

表2 目标频率范围为100～1 000 Hz时的最佳结构参数Tab.2 The optimal structural parameters for the target frequency range of 100～1 000 Hz

表3 目标频率范围为100～1 500 Hz时的最佳结构参数Tab.3 The optimal structural parameters for the target frequency range of 100～1 500 Hz

表4 目标频率范围为100～2 000 Hz时的最佳结构参数Tab.4 The optimal structural parameters for the target frequency range of 100～2 000 Hz

3 有限元仿真

为了提高验证效率并降低实验成本，采用有限元仿真对目标频率范围可变的变参数分空腔单层微穿孔板结构的吸声性能进行了初步验证。在虚拟声学实验室［9］中建立如图3所示的有限元仿真模型，驻波管的尺寸设置为60 mm×60 mm×300 mm，将其划分为5 mm的三角形网格。传声器1到微穿孔板的距离为90 mm，传声器2到微穿孔板的距离为40 mm。将驻波管模型的左侧面定义为声波入射面，微穿孔板为一层与表1～4中微穿孔板板厚相同的空气层。由于微穿孔板小孔的存在，使得划分有限元的网格变得非常困难，而且网格质量很差，网格的数量和计算量非常巨大。为了解决小孔的问题，在建立声学有限元网格时，通常忽略小孔，在声学计算的时候，在微穿孔两边的网格之间通过定义一种传递导纳关系，间接模拟这些小孔。所以，本次仿真根据表1～4中的最佳结构参数计算出微穿孔板两侧的传递导纳关系来模拟微穿孔板吸声机理［10-12］。

图3 单层变孔径微穿孔板结构的有限元仿真模型Fig.3 Finite element simulation model of the single-layer microperforated panel structure with variable parameters and separated cavities

将变参数分空腔单层微穿孔板结构平均吸声系数理论值和仿真值进行比较，可以评估最佳结构参数的准确性和布谷鸟搜索优化算法的有效性。

4 驻波管测试

为检验有限元仿真的有效性和准确性，根据表1～4中的最佳结构参数，利用精密激光束加工技术制备出不同频率范围内优化后的变参数分空腔单层微穿孔板结构，通过驻波管测量样品的垂直入射吸声系数［13-16］。实验中使用的是AWA6128A型驻波管测试仪，结构如图4所示。将测试样品安装在被测材料管中，调整刚性活塞获得板后空腔长度，驻波管的另一端连接扬声器，将正弦波信号通过扬声器放大垂直射向待测样品，通过变参数分空腔单层微穿孔板结构反射的回波则由安装在滑车上的测试探头接收，通过调节滑车在滑道上的位置，可以获得一定频率的峰值声级和相应的谷值声级，根据驻波比可以测出样品的吸声系数。

图4 AWA6128A型驻波管测试仪结构示意图Fig.4 Structure diagram of AWA6128A standing wave tube tester

变参数分空腔单层微穿孔板结构的实际吸声系数可以通过计算机中的控制软件计算得到［13-16］。通过这种方式，获得具有可变目标频率范围的变参数分空腔单层微穿孔板结构的实际吸声系数，可以为构建的理论吸声模型、选择的布谷鸟搜索优化算法以及利用的有限元仿真方法提供实验验证。

5 结果与讨论

图5中展示了这些最佳变参数分空腔单层微穿孔板结构在可变目标频率范围内的吸声系数理论值、仿真值和实验值。为避免多组数据之间出现重叠和混淆，图5根据结构参数组数来划分。优化变参数分空腔单层微穿孔板结构的平均吸声系数理论值、仿真值和实验值对比如表5所示。

表5 优化变参数分空腔单层微穿孔板结构平均吸声系数对比Tab.5 Comparison of average sound absorption coefficient of single-layer microperforated panel structure with optimized variable parameters and separated cavities

图5 最佳变参数分空腔单层微穿孔板结构的吸声系数Fig.5 The sound absorption coefficient of the optimal single-layer microperforated panel structure with variable parameters and separated cavities

由图5可以看出，对于不同的组数和频率范围，优化后的吸声系数理论值、仿真值和实验值的一致性高，证明了本研究所建立的理论吸声模型、选择的布谷鸟搜索优化算法和使用的有限元仿真方法是准确的。由表5可以看出，4种结构参数微穿孔板吸声结构的平均吸声系数更优于2种以及3种结构参数微穿孔板吸声结构，这说明变参数分空腔单层微穿孔板结构的吸声性能与不同结构参数的组数多少密切相关。

6 结论

1）理论数据、仿真数据和实验数据之间的一致性高，证明了理论吸声模型、布谷鸟搜索优化算法和有限元仿真的准确性。

2）在总厚度为20 mm的限制条件下，对于目标频率范围100～500 Hz，100～1 000 Hz，100～1 500 Hz和100～2 000 Hz，相应的变参数分空腔单层微穿孔板结构的最佳实际平均吸声系数分别为0.662 6，0.745 6，0.787 5和0.826 7，表现出了出色的低频吸声性能。

3）随着结构参数种类的增加，变参数分空腔单层微穿孔板结构的吸声频带变宽，3组结构参数分空腔单层微穿孔板结构和4组结构参数分空腔单层微穿孔板结构在400～2 000 Hz的频率范围内的吸声系数超过0.8。

4）与传统微穿孔板吸声结构相比，优化变参数分空腔单层微穿孔板结构的低频吸收性能显著提高，吸声频带宽，因而在装备吸声降噪领域具有很大的实际应用价值。