多频带声学超构材料的透射性能

郝丽梅 杨宵 闫小乐 程红 李育佳 王雅喆 解忧

摘 要：为了研究多频带声学超构材料的透射性能，采用声学有限元方法（FEM）模拟不同长度的OE-CE（一端封闭一端连通）空心方管结构单元的声透射谱，分析讨论管长对透射谱的影响，设计出一种占空比较小的由不同长度OE-CE空心方管并联组合的多频带声学超构材料板。结果表明：嵌有9个空心方管声学超构材料的声透射曲线出现声透射阻带，随着空心方管的长度增加，透射阻带峰值和负等效质量密度频带向低频方向移动；透射阻带的个数与不同长度的空心方管数目有关，即2种不同长度的空心方管并联组合的声学超构材料板可以实现2个透射阻带，3种不同长度的则可以实现3个透射阻带；宽频带声学超构材料板则需要在较低频率区域增加对应长度的空心方管数目，在较高频率区域减小相邻空心方管的长度梯度。所提出的采用OE-CE结构单元构建的声学超构材料具有结构简单、响应频带低的优点，并在多个响应频带具有优异的声学透射阻带特性。

关键词：并联型声学超构材料；透射曲线；多频带；OE-CE空心方管

中图分类号：TB 132文献标志码：A

文章编号：1672-9315（2023）03-0631-06

DOI：10.13800/j.cnki.xakjdxxb.2023.0321开放科学（资源服务）标识码（OSID）：

Transmission property of multi-band acoustic metamaterialsHAO Limei，YANG Xiao，YAN Xiaole，CHENG Hong，LI Yujia，WANG Yazhe，XIE You

（College of Science，Xian University of Science and Technology，Xian 710054，China）

Abstract：In order to study the transmission properties of parallel acoustic metamaterial with multi band，the structure unit of the OE-CE（one end closed and the other connected）hollow square tube was simulated by the finite element method（FEM），and the influence of tube length on the transmission spectrum was discussed.Multi-band acoustic metamaterials plate with low duty cycle was designed by structural units of OE-CE hollow square tubes of different lengths in parallel.The results show that the acoustic transmission spectrum of acoustic metamaterials embedded with nine hollow square tubes has a stop-band.As the length of the hollow square tube increases，the peak of the stop band and frequency region of negative effective mass density would move to low frequency.It is found that the number of the stopbands is related to the number of hollow square tubes with different lengths，that is，two stopbands can be achieved by the acoustic metamaterials with parallel hollow square tubes of two different lengths，and three stopbands done by three different lengths.The broadband acoustic metamaterial needs to increase the number of the corresponding hollow square tubes in the lower frequency region，and reduce the length gradient of adjacent hollow square tubes in the higher frequency region.The OE-CE unit has the advantages of simple structure and low response frequency band，which exhibits an excellent transmission stopband property in multiple frequency bands.

Key words：parallel acoustic metamaterial；transmission；multiband；OE-CE hollow square tube

0 引 言

聲学超构材料是一种人工合成的周期性结构材料，它的结构尺寸远远小于材料的工作波长，具有很多奇异特性，如负等效模量^[1]、负等效质量^[2]和负折射率^[3]等。重要的是它还可以控制声波的传播方向，近年来在声隐身^[4]、亚波长成像^[5]、超透镜^[6]、反常多普勒效应^[7]、超吸收^[8]、聚焦^[9]、振动隔声^[10]等方面有显著的应用潜力。近20 a，声学超构材料的研究受到人们的广泛关注。

声学超构材料的核心是设计人工微结构实现对声波的有效调控。2000年，LIU等将铅球包裹在硅橡胶中，该局域共振型结构可以等效为质量-弹簧模型，在400 Hz和1 350 Hz这2个频率附近实现了声透射阻带和负等效质量密度^[2]。YAO等基于质量-弹簧模型，设计丁苯橡胶填充钢格栅弹性超材料实现不同频带的声透射阻带^[11]。这些材料均采用局域共振型结构原理设计微结构单元实现负等效质量密度。2010年，YANG等将质量块固定在张紧的薄膜上制得薄膜型声学超构材料，在200～300 Hz实现了声透射阻带^[12]。MEI等利用带有多块小板的薄膜，实现较宽频的声透射阻带^[13]。FAN等提出一种开孔薄膜的管状一维声学超构材料，在400～600 Hz附近出现透射阻带^[14]。LIN等设计一种橡胶薄膜型超材料，该薄膜可以看作是同时具有阻尼和弹簧双重作用的负质量密度薄膜。薄膜型超材料表现出了其结构单元偶极共振的物理学特征^[15]。

声学超构材料的声波传输调控依赖于结构单元的形态及参数变化，通常采用偶极共振模式^[16]。目前，研究者们旨在优化超材料异构结构单元数目来实现多频或宽频的声波透射调控及负质量密度。例如，LEE等利用多规格薄膜阵列表现出宽频的声透射阻带^[17]。YANG等通过将多个模板堆叠实现宽频带的声透射阻带^[18]。CHEN等以同心环结构来构建二维单振子、双振子和三振子带有耦合作用的声学超构材料，在多个频段出现声透射阻带^[19]。REN等提出一种轻质亚波长超薄多缝的结构，在降低系统整体质量的基础上能够在较宽的频率范围内吸收声音^[20]。FAN等提出了一种基于膜涂层穿孔板的声学超构材料，以实现在低频范围内的隔音^[21]。XU等提出一种双层多孔薄板声学超构材料，通过调节孔径大小来调控声透射阻带的位置^[22]。ZHAO等将很多不同长度的空心管组成双层结构声学超构材料，其中在介质分别为水和空气的条件下，均实现了宽频带的声透射阻带^[23-24]。通过增加额外质量结构单元的形式拓宽了声学超构材料的响应频带，但存在制备复杂、应用不便等问题。

文中基于空心管结构^[23-26]，将两端均与介质连通的空心管优化为一端封闭一端连通的（OE-CE）空心方管，在薄的声学板上内嵌OE-CE空心方管结构构建多频带声学超构材料。采用声学有限元方法（FEM）模拟多种嵌有OE-CE空心方管的声学超构材料板，仿真研究不同空心方管數目对透射性能的影响规律。最后通过组合不同长度的空心方管设计出并联型多频带声学超构材料板。

1 模型构建

1.1 结构单元设计与模拟

将OE-CE空心方管并列组成声学超构材料板单元，板内部嵌有多个空心方管，一侧与空气连通，另一侧是封闭的薄壁，如图1所示。OE-CE空心方管的初始结构参数设定为：边长a₁=a₂=8 mm，薄壁厚度t=1 mm，相邻OE-CE空心方管的间距u=2 mm。声学薄板的几何参数：长L₁=100 mm，宽L₂=90 mm，厚度d=10 mm，这里的厚度仅有大约最短工作波长的八分之一，因此称之为声学薄板。

在外界声波的激励下，OE-CE空心方管中的空气在谐振频率附近发生谐振，从而使得该结构单元呈现出负等效质量密度的声学特性。

利用COMSOL Multiphysics的声-结构模块在频域内对结构单元进行模拟计算。在波导管中，设置幅值为1 Pa的平面波从波导管一端入射，另一端设为匹配边界，可实现波导管内声波传递的完全吸收，避免声波来回振荡影响计算结果。4个边界为周期性边界条件，这时每个板可以被单独看作一个在空间内无限衍生的单元结构。当声波经过样品时，通过4个不同的探测面采集对应处声波的幅值和相位，利用传递函数法解耦得到透射系数T和反射系数R，再对其数据进行后处理得到等效质量密度等声学参数。

在COMSOL仿真中，设置固体域材料为Steel AISI 4340，其余域为流体域，设置流体域材料为空气，声速为343 m/s，密度为1.29 kg/m³。

1.2 长度对结构单元透射性能的影响

为了探究空心方管长度L₁对结构单元透射性能的影响规律，该部分以嵌有9个空心方管的声学薄板为研究对象，选用能反映出谐振频率随空心方管长度L₁变化规律的长度，即L₁分别为30，40，50，60，70 mm，并利用Comsol仿真软件对其进行仿真研究，最终得出不同长度声学超构材料板的频率-透射曲线（图2（a））。

图2（a）可以看出，不同长度声学超构材料板的透射谱中均出现一个明显的透射阻带。随着空心方管长度的增加，透射阻带出现红移，且阻带峰值逐渐升高。例如当L₁分别为30，40，50，60，70 mm时，对应的声学薄板谐振频率分别为2 740，2 070，1 665，1 395，1 200 Hz，透射阻带最低峰值分别为0.23，0.28，0.33，0.40，0.46。其中出现红移现象的原因可能是：空心方管的长度增加，使得等效声感（L_eff∝ρL₁/S）增加；同时管腔体积的增加引起等效声容（C_eff∝Vρc²）的增加。因此对应的谐振频率（f∝1L_effC_eff）降低^[25]。

图2（b）为不同长度声学超构材料板的等效质量密度曲线。不同长度的情况下，空心方管的等效质量密度在各自的谐振频率附近（即透射阻带）均为负值。其原因可能是空心方管内存储了大量声能，其中的空气柱以谐振频率振动，而外界声场驱动力与空气柱加速度方向相反^[23]。

2 结果与讨论

当空心方管一端与空气连通，另一端用薄壁封闭时，空气柱只能在连通端进入空心方管，在管腔内来回振动。结合声电类比方法分析，对比两端都连通的空心管，OE-CE空心方管具有更高的电感，谐振频率更低，选用OE-CE空心方管设计在低频范围内多频带声学超构材料有较大的优势。因此，文中主要研究不同长度OE-CE空心方管的数目对声学超构材料板的频带个数的影响规律。

2.1 双频带并联型声学超构材料的透射性能

对于由2种不同长度的空心方管并联组合形成双频带声学超构材料，其薄壁厚度t为1 mm，薄壁将连通的腔体分割为2个独立的OE-CE空心方管，板长L₁为2种空心方管的长度之和（图3）。

例如将40 mm和70 mm长度的空心方管进行组合，则声学超构材料的板长L₁为110 mm。其余参量和上一节相同。当板长L₁的长度和薄壁位置发生变化时，空心方管的组合也随之改变。

图4为2种不同长度组合下的声学薄板的透射曲线。由图可知，5种不同组合的情况下，其透射曲线均出现2个透射阻带，且相比于单独同长度的空心方管，2种长度组合的板型声学超构材料的谐振频率点均出现了少许红移。以40 mm和50 mm这2种长度组合的声学薄板为例，其透射曲线出现2个透射阻带，阻带峰值分别为2 025 Hz和1 640 Hz，与单独的40 mm（2 070 Hz）和50 mm（1 665 Hz）声学薄板相比，谐振频率出现少许红移。结合声电类比分析其原因可能为相邻的空心方管出现相互耦合作用，使得等效声感增加，从而导致谐振频率的降低。

类似的，由40 mm和60 mm、40 mm和70 mm、50 mm和60 mm、50 mm和70 mm组合的声学薄板对应的谐振频率分别为2 020 Hz和1 375 Hz，2 005 Hz和1 185 Hz，1 630 Hz和1 370 Hz，1 630 Hz和1 185 Hz。另外，空心方管的长度越长，透射阻带峰值的红移越明显，这主要是由于长度增加使得耦合作用增强所致。

因此，2种不同长度的空心方管组合的声学薄板能够实现双频带的声透射阻带。

2.2 三频带并联型声学超构材料的透射性能

基于双频带声学薄板，设计了3种不同长度，即40，50 mm和60 mm的空心方管組合构建的声学薄板（图5）。其中，板长L₁为100 mm，薄壁将每个腔体分割为2个独立的OE-CE空心方管（即谐振腔），声学薄板中共有18个谐振腔，中间6个谐振腔的长度为50 mm，左上及右下的6个谐振腔长度为40 mm，右上及左下的6个谐振腔长度为60 mm。

图6为3种不同长度组合的声学薄板的透射曲线，由图可知，每一种声学薄板的透射曲线均出现3个透射阻带，阻带峰值分别为1 380，1 620，2 020 Hz。从图6还可以看出，在每个透射阻带附近都伴随相位突变，这进一步说明空气柱在阻带峰值附近发生谐振。另外，与单独的40 mm（2 070 Hz），50 mm（1 665 Hz）和60 mm（1 395 Hz）声学薄板相比，其峰位出现少许红移，这是由于相邻的空心方管出现相互耦合作用所致。

因此，3种不同长度的空心方管组合的声学薄板能够实现三频带的声透射阻带。

2.3 多频带并联型声学超构材料的透射性能

为研究多频带声学薄板超材料，设计10种不同长度，即32，34，36，38，40，60，62，64，66 mm和68 mm的空心方管组合构建的声学薄板，如图7所示。

图8为多长度组合下的声学薄板的透射曲线。由图可知，低频处出现1个宽的透射阻带，高频处出现5个窄的透射阻带，这6个透射阻带与单独不同长度的空心方管对应的谐振频率大致相符。从图8还可以看出，在每个透射阻带附近都伴随相位突变，说明在每一个透射阻带附近均发生了谐振。

进一步分析多长度组合下的声学薄板透射曲线。在低频处，透射阻带的最低峰值较高；在高频处，透射阻带的最低峰值较低。其原因可能是：每种尺寸的谐振腔只有2个，不同尺寸的空心方管对声波的吸收能力并不相同，且谐振腔在低频处的吸声能力弱，在高频处的吸声能力强，使得相同数目的空心方管在低频处的透射阻带幅值较小，导致透射阻带的最低峰值较高。因此，在设计组合声学薄板时需要增加较低谐振频率处对应尺寸的空心方管数目。另外，从图8还可以看出，间隔为2 mm的空心方管在高频处出现独立的透射阻带。因此，在高频处设计一定带宽的透射阻带时则需要缩小相邻空心方管的长度梯度。

综上，在设计多频带或宽频带的声学超构材料板时，在较低频率区域需要增加对应尺寸的空心方管数目，在较高频率区域需要减小相邻空心方管的长度梯度。

3 结 论

1）对于嵌有不同长度的空心方管的声学超构材料板：不同管长的空心方管的声学薄板均能够表现出透射阻带和负等效质量密度特性。随着方管长度的增加，透射阻带出现红移，且阻带峰位升高。

2）透射阻带的频带个数与不同长度的空心方管种类数目有关。由2种不同长度的空心方管组合的声学薄板，在其对应的2个谐振频率附近均表现出明显的声透射阻带。其中，空心方管的长度越长，透射阻带的红移现象越明显。而由3种不同长度的空心方管组合的声学薄板则能够在对应的3个频带附近实现声透射阻带。

3）对于实现宽频带声学薄板：在较低频率区域需要增加对应尺寸的空心方管数目，而在较高频率区域则需要减小相邻空心方管的长度梯度。

參考文献（References）：

[1] FANG N，XI D J，XU J Y，et al.Ultrasonic metamaterials with negative modulus[J].Nature Materials，2006，5（6）：452-456.

[2]LIU Z Y，ZHANG X X，MAO Y W，et al.Locally resonant sonic materials[J].Science，2000，289（5485）：1734-1736.

[3]SONG X P，CHEN T N，HUANG W K，et al.Frequency-selective modulation of reflected wave fronts using a four-mode coding acoustic metasurface[J].Physics Letters A，2021，394：127145.

[4]ZHOU H T，FU W X，WANG Y F，et al.Ultra-broadband passive acoustic metasurface for wide-angle carpet cloaking[J].Materials and Design，2021，199：109414.

[5]ZHU J，CHRISTENSEN J，JUNG J，et al.A holey-structured metamaterial for acoustic deep-subwavelength imaging[J].Nature Physics，2011，7（1）：52-55.

[6]KAINA N，LEMOULT F，FINK M，et al.Negative refractive index and acoustic superlens from multiple scattering in single negative metamaterials[J].Nature，2015，525（7567）：77-81.

[7]ZHAI S L，ZHAO X P，LIU S，et al.Inverse doppler effects in broadband acoustic metamaterials[J].Scientific Reports，2016，6（1）：32388.

[8]GUO J W，ZHANG X，FANG Y，et al.An extremely-thin acoustic metasurface for low-frequency sound attenuation with a tunable absorption bandwidth[J].International Journal of Mechanical Sciences，2022，213：106872.

[9]WU Y，LAI Y，ZHANG Z Q.Elastic metamaterials with simultaneously negative effective shear modulus and mass density[J].Physical Review Letters，2011，107（10）：105506.

[10]CUI S C，HARNE R L.Soft materials with broadband and near-total absorption of sound[J].Physical Review Applied，2019，12（6）：064059.

[11]YAO S S，ZHOU X M，HU G K.Investigation of the negative-mass behaviors occurring below a cut-off frequency[J].New Journal of Physics，2010，12（10）：103025.

[12]YANG Z，MEI Jun，YANG Min，et al.Membrane-type acoustic metamaterial with negative dynamic mass[J].Physical Review Letters，2010，101（20）：204301.

[13]MEI J，MA G C，YANG M，et al.Dark acoustic metamaterials as super absorbers for low-frequency sound[J].Nature Communications，2012，3（1）：1-7.

[14]FAN L，GE H，ZHANG S Y，et al.Research on pass band with negative phase velocity in tubular acoustic metamaterial[J].Journal of Applied Physics，2012，112（5）：1734-1736.

[15]LIN G C，CHEN S Q，LI Y L，et al.Transmission loss of membrane-type acoustic metamaterial with negative effective mass[J].Materials Science Forum，2017，898：1749-1756.

[16]LIU Z Y，CHAN C T，SHENG P.Analytic model of phononic crystals with local resonances[J].Physical Review B，2005，71（1）：014103.

[17]LEE S H，PARK C M，SEO Y M，et al.Acoustic metamaterial with negative density[J].Physics Letters A，2009，373（48）：4464-4469.

[18]YANG Z，DAI H M，CHAN N H，et al.Acoustic metamaterial panels for sound attenuation in the 50～1 000 Hz regime[J].Applied Physics Letters，2010，96（4）：041906.

[19]CHEN M，MENG D，JIANG H，et al.Investigation on the band gap and negative properties of concentric ring acoustic metamaterial[J].Shock and Vibration，2018，2018：1369858.

[20]REN S W，MENG H，XIN F X，et al.Ultrathin multi-slit metamaterial as excellent sound absorber：Influence of micro-structure[J].Journal of Applied Physics，2016，119（1）：014901.

[21]FAN L，CHEN Z，ZHANG S Y，et al.An acoustic metamaterial composed of multi-layer membrane-coated perforated plates for low-frequency sound insulation[J].Applied Physics Letters，2015，106（15）：151908.

[22]XU Y C，WU J H，CAI Y Q，et al.Study on dynamic effective parameters of bilayer perforated thin-plate acoustic metamaterials[J].Modern Physics Letters B，2021，35（2）：2150048.

[23]CHEN H J，ZHAI S L，DING C L，et al.Meta-atom cluster acoustic metamaterial with broadband negative effective mass density[J].Journal of Applied Physics，2014，115（5）：054905.

[24]CHEN H J，ZHAI S L，DING C L，et al.Acoustic metamaterial with negative effective mass density in water[J].Journal of Applied Physics，2015，118（9）：094901.

[25]闫小乐，郝丽梅，叶临峰，等.模型参数对开口空心球超材料透射行为的影响[J].西安科技大学学报，2020，40（3）：512-517.

YAN Xiaole，HAO Limei，YE Linfeng，et al.Effects of model parameters on the transmission property of acoustic metamaterial with split hollow sphere[J].Journal of Xian University of Science and Technology，2020，40（3）：512-517.

[26]CHEN H J，ZENG H C，DING C L，et al.Double-negative acoustic metamaterial based on hollow steel tube meta-atom[J].Journal of Applied Physics，2013，113：104902.

（責任编辑：高佳）