一種基于局域共振的低頻超寬帶隙瓣狀聲學超材料

陳琳，吳衛國，周榕

?

一種基于局域共振的低頻超寬帶隙瓣狀聲學超材料

陳琳，吳衛國，周榕

(江蘇大學土木工程與力學學院，江蘇鎮江212013)

針對低頻聲波的衰減問題，設計了一種基于局域共振(Local Resonance，LR)機理的帶瓣型結構的聲學超材料，并運用COMSOL 有限元軟件MULTIPHYSICS 4.3計算分析了所設計的結構元胞帶隙特性及其振動特性，研究表明，所設計的聲學超材料比無瓣型聲學超材料帶隙更寬，并且僅用單層結構同種元胞組合就可以在一定頻率下形成較寬的多個完全帶隙，而不需要通過不同種元胞組合或是多層結構的復合而達到一定帶隙特性。在此基礎上，進一步對所設計的聲學超材料的幾何尺寸進行了優化。

聲學超材料；局域共振；聲波帶隙

1 引言

低頻噪聲因其超強的穿透能力，一直是噪聲控制領域的難點之一，傳統的被動隔聲設備對500 Hz以上的中高頻噪聲有較好的隔聲效果，對于低頻噪聲由于受質量定律的限制，其隔聲效果不太理想，如對于300 Hz以下的低頻噪聲，根據質量定律，隔聲部件的厚度達到1 m才能滿足隔聲部件尺寸與聲波波長處于同一數量級，起到隔聲效果。這樣的隔聲部件太笨重，難以滿足航空、軍事、高速列車中隔聲設備小巧高效的要求。1993年Kushwaha M S和Halevi P[1]等人提出了聲子晶體，文獻[2]表明聲子晶體能帶結構有帶隙存在，在帶隙范圍內可以產生很好的隔聲效果，但這些研究工作仍然沒有打破質量定律，其低頻減振降噪效果仍不夠理想，一般在1~4 kHz范圍內效果較好。2000年，Liu[3]等人提出了局域共振機理，其設計的聲學結構可以控制波長大于晶格尺寸兩個數量級的聲波，實現了毫米級結構對大波長低頻聲波的有效控制。

近年來研究者們設計了多種基于局域共振機理的聲學超材料結構，包括一維膜型[4]，二維二組元膜型[5]或是板型[6-7]，二維三組元膜型[8]，三維三組元體型[3]，二維Helmhotz腔型[9]，一維超材料梁型[10]，以及二維膜型和腔型的結合[11]等。這些結構在一定頻率范圍內都可以產生禁帶，但是這些聲學超材料結構如果僅憑借一種元胞而組成隔聲裝置其帶隙并不是很寬，相對帶隙(?f/f)也僅為15%左右，很多都是通過多層結構堆疊或是改變元胞組成元素從而形成復合元胞達到較寬帶隙的[8,12]。本文提出了一種基于局域共振機理的瓣型聲學超材料微結構，并研究其元胞的帶隙特性以及帶隙優化規律。

2 模型的建立

本文設計的模型如圖1所示，模型設計借鑒了向內彎曲的花瓣受到外力時產生較大的彈性曲率，能量從而減小，花朵振動進而保護花蕊的機理。圖1(a)是本文設計的一個超材料元胞的模型，由半球、“瓣”、薄膜、柵格支架四部分組成。元胞為隔聲裝置中阻隔低頻噪聲的最小聲學超材料單元，圖1(b)是元胞背面和柵格；圖1(c)是元胞陣列后的超材料結構；圖1(d)是陣列后的背面結構圖。其中，柵格支架由正方形格子沿方向和方向周期性延拓而成，起著固定支撐的作用，使得其上的質量有能力振動。中間的彈性薄膜相當于“彈簧-質量”系統中的彈簧，半球和“瓣”相當于質量塊。

(a)???????(b)

(c) ???????(d)

圖1 “瓣”狀結構的聲學超材料模型

Fig.1 Schematic of the petal-like acoustic metamaterial model

3 數值計算及分析

由于該聲學超材料是周期性的，本文首先計算了一個元胞的特征頻率及帶隙特性。所研究元胞結構尺寸為：半球半徑的范圍為2.6~3.4 mm，“瓣”的厚度統一為0.5 mm，薄膜厚度為0.4~0.7 mm，剛性支架長度也就是晶格常數為9~14 mm，支架高度和厚度分別為1 mm和0.75 mm。粘貼于彈性薄膜的半球是金屬鎢，彈性薄膜是硅膠，剛性支架是硬質塑料。材料參數：硅膠密度=1300 kg/m3；彈性模量=0.1175 MPa；泊松比0.469；塑料密度=1190 kg/m3；彈性模量=2.2 GPa；泊松比0.375；鎢密度=17800 kg/m3；彈性模量=360 GPa；泊松比0.27。

3.1 能帶結構計算

本文采用的元胞尺寸為=3.4 mm，=10 mm，=0.6 mm，元胞四周的邊界均設置為Bloch周期性邊界條件，其他邊界為自由邊界，本文利用有限元軟件COMSOL MULTIPHYSICS 4.3對元胞的整個不可約布里淵(Brillouin)區邊界---進行掃描得到元胞的能帶圖，如圖2所示。從圖2中可看出，在相同尺寸、相同材料下，有“瓣”和無“瓣”[13]元胞的帶隙特性相差很大，對比圖2(a)和圖2(b)可知，有“瓣”模型產生的完全帶隙的個數和帶寬都比無“瓣”模型要優越，有“瓣”模型可以產生3條完全帶隙，而無“瓣”模型只產生很窄的一條完全帶隙。圖2(a)中黃色區域是包括彎曲波[13]帶隙在內的第一條帶隙，用(1)表示，藍色區域表示以“平帶”[7]為上邊界的第二條完全帶隙，用(2)來表示，綠色區域為最寬的以“平帶”作為下邊界的第三條完全帶隙，用(3)來表示。由圖2中的能帶圖可計算出瓣型結構產生的帶隙寬度是無“瓣” 模型的16倍([(305-140)+(92-76)]/(57-46)=16.4)，并且在400 Hz以下有三條完全帶隙，帶隙寬度達254 Hz，占作用總頻率的65%，相對帶隙(?f/f)高達到60%。

(a)有“瓣”元胞

(b)無“瓣”元胞

圖2 有“瓣”和無“瓣”元胞的能帶圖

Fig.2 Sonic band structures of the metamaterial cell: a. with "petals" and b. without "petals"

3.2 結構模態分析

由圖2可見，瓣型結構比較寬的帶隙，尤其是第三完全帶隙(圖中用(3)表示)的形成與其結構的振動模態特性密不可分。圖3所示4張圖分別為圖2(a)所示能帶圖中曲線A、B、F、G對應固有頻率在點的振動模態圖。圖3(a)對應圖2(a)中曲線A在點的振動模態，對應于元胞結構的第一階模態，固有頻率為62 Hz，從圖中可以看出半球和瓣都沿方向振動，而四周框架幾乎保持靜止。這表明，硬質塑料框架可以看做剛性基礎，起到隔離每個元胞的作用，使每個元胞的振動都局域化，該階模態與文獻[6]中的“elongation”模態相對應。圖3(b)對應圖2(a)中曲線B在點的振動模態，對應于元胞結構的第二階模態，固有頻率為 76 Hz，半球和花瓣在水平方向或方向振動。由于曲線C在點的的振動模態與曲線B在點的相似，都是半球和花瓣在水平或方向振動，故只用圖3(b)代表元胞的第二階振動模態。由于結構在方向和方向的對稱性，結構在方向和方向的振動也相似，只是方向不同，所以曲線B、C在遠離點處幾乎是重合的，曲線D、E在遠離Γ點處也幾乎是重合的，該模態與文獻中的“shear”模態相對應；圖3(c)對應圖2(a)中曲線F在點的振動模態，對應于元胞結構的第三階模態，固有頻率為146 Hz，從圖中可以看出半球和瓣繞著幾何對稱軸做扭動。曲線F代表“平帶”，其值不隨波矢的改變而變化，該模態與其“breath”模態相對應。最后，圖3(d)所示模態與曲線G對應，對應于胞元結構的第四階模態，固有頻率為332 Hz，只有瓣的水平相向振動，半球保持不動。

(a)

(b)

(c)

(d)

圖3 圖2(a)中點A、B、F、G對應的振動模態

Fig.3 The vibration modes corresponding to the points of A, B, F, G shown in Fig.2

從結構模態分析可以看出，前三種振型都是框架保持不動，只有半球、瓣和薄膜振動。這樣，當整個結構受到振動干擾時，塑料框架不動，薄膜與其上的結構振動耗散能量，從而達到減振的效果；當結構受到來自空氣的垂直入射聲波激勵時，如果激勵頻率與結構的固有振動頻率接近時，則聲波與結構發生強烈的耦合作用，從而達到降噪的效果。圖3(d)所示模態與曲線G對應，只有瓣的水平相向振動，半球保持不動，該振型特征是超寬帶隙形成的主要原因之一。

3.3 結構參數優化研究

由圖2(a) 和2(b) 的對比得知，瓣型模型比無“瓣”模型的帶隙特性優越。但是，有“瓣”不是產生優越帶隙特性的充分條件，元胞中各結構幾何尺寸在一定范圍內滿足相互匹配才能達到最佳效果。

(a) 晶格常數對帶隙特性的影響

首先分析晶格常數對帶隙特性的影響，計算中薄膜厚度=0.6 mm和半球半徑=3.4 mm都保持不變，晶格常數分別取9、10、11、12、13、14 mm時，計算的能帶結構如圖4所示。

從圖4(a)~4(f)可以看出，隨著晶格常數的逐漸變大，能帶結構中(3)區域寬度越來越小，從232 Hz逐漸減小到8 Hz，當晶格常數達到12 mm時，(3)區域消失，而(2)區域從無到有再到無，(1)區域寬度基本不變，帶隙中心頻率逐漸降低，具體變化數量及趨勢如圖5所示，圖中主要標注了第三完全帶隙的寬度，第一、二完全帶隙寬度變化不大，都在25Hz以下，圖上沒有標注。

(a)=9 mm

(b)=10 mm

(c)=11 mm

(d)=12 mm

(e)=13 mm

(f)=14 mm

圖4 超材料晶格常數變化時結構的能帶圖

Fig.4 The sonic band structures of the metamaterial cells with different lattice constants

(b) 薄膜厚度和半球半徑對帶隙特性的影響

半球在元胞結構中相當于質量塊，半球半徑的大小等價于質量塊質量的大小。彈性薄膜在結構中相當于緩沖振動的彈簧，薄膜的厚度影響著薄膜的彈性，所以也會對結構帶隙特性產生一定影響。在“彈簧-質量”系統中，彈簧和質量任何一方發生改變，其振動特性就會受到影響，固有頻率就會有所變化，對應元胞的能帶結構就會發生變化。本文在保持晶格常數=10 mm不變的前提下，分別計算了薄膜厚度和半球半徑變化時元胞的帶隙特性。計算結果如圖6所示，從圖6(a)中可看出，隨著薄膜厚度的增加，(3)區域寬度和中心頻率都逐漸增大，(1)、(2)區域寬度基本保持不變，其中心頻率均隨薄膜厚度增加而逐漸增加。

從圖6(b)中可看出，(3)區域寬度隨半球半徑R增大變化較大，從9 Hz迅速變到141 Hz；(2)區域寬度基本不變，其中心頻率反而逐漸降低；(1)區域寬度基本不變，中心頻率逐漸升高。

由以上的分析知，要得到最寬的帶隙，就要優先考慮(3)區域(第三完全帶隙)的寬度，因為(3)區域占總完全帶隙的比例(86%)最大，從(3)區域寬度變化趨勢來看，就要盡量減小晶格常數，增大薄膜厚度，增大半球半徑。但是，縮小晶格常數和增大半球半徑是兩個矛盾的過程。因為半球半徑增大，附屬在其外的瓣范圍也將擴大，有可能超出晶格常數的范圍，所以必須找到最合適半球半徑和晶格常數的搭配，才能得到最優帶隙特性。本研究中結構尺寸為=9 mm，=3.4 mm，=0.6 mm時(3)區域達到最寬(232 Hz)，如圖4(a)所示。此時晶格常數較小，瓣形所占空間比相對增大，瓣對帶隙的影響隨之增大，故禁帶增寬也在預料之中。

(a)帶隙隨薄膜厚度的變化

(b)帶隙隨半球半徑的變化

圖6 薄膜厚度和半球半徑對帶隙特性的影響

Fig.6 The relationships of bandwidth with (a) film thickness and (b) hemisphere radius

4 結論

本文設計了一種基于局域共振(LR)機理的帶瓣型結構的聲學超材料，模型設計借鑒了向內彎曲的花瓣受到外力時產生較大的彈性曲率能量從而減小花朵振動進而保護花蕊的機理。通過COMSOL計算研究表明本文所設計的聲學超材料比無瓣型聲學超材料帶隙更寬，并且僅用單層結構同種元胞組合就可以在一定頻率下形成較寬的多個完全帶隙，而不需要通過不同種元胞組合或是多層結構的復合而達到一定帶隙特性，在此基礎上進一步對所設計的聲學超材料幾何尺寸進行了優化，當結構尺寸為(=9 mm，=3.4 mm，=0.6 mm)時(3)區域達到最寬(232 Hz)，此時晶格常數較小，瓣形所占空間比相對增大，瓣對帶隙的影響隨之增大。

[1] Kushwaha M S, Halevi P, Dobrzynski L, et al. Acoustic band structure of periodic elastic composites [J]. Phys. Rev. Lett, 1993, 71(13): 2022-2025.

[2] Economou E N, Sigalas M M.Classical wave propagation in periodic structure[J]. Phys. Rev. B, 1993, 48(18): 13434-13438.

[3] LIU Zhengyou, ZHANG Xixiang, MAO Yiwei, et al. Locally resonant sonic materials[J]. Science, 2000, 289(5485): 1734-1734-6.

[4] 高東寶. 基于聲學超材料的新型隔聲技術研究[D]. 長沙: 國防科技大學, 2013, 74-87.

GAO Dongbao. New acoustical technology of sound insulation based on acoustic metamaterials[D]. Changsha: National University of Defense Technology, 2013, 74-87.

[5] MEI Jun, MA Guancong, YANG Min, et al. Dark acoustic metamaterials as super absorbers for low-frequency sound[J]. Nature Communication, 2012, 10(41): 425-433.

[6] Badreddine Assouar M, Mourad Oudich. Enlargement of a locally resonant sonic band gap by using double-sides stubbed phononic plates[J]. Appl. Phys. Lett, 2012, 100(12): 1235061-123506-4.

[7] Mourad Oudich1, Li Yong, Badreddine M Assouar, et al.A sonic band gap based on the locally resonant phononic plates with stubs[J].New Journal of Physics, 2010, 12(08): 0830490-0830499.

[8] 吳健, 王剛, 白曉春, 等. 一種輕質低頻寬帶薄膜超材料隔聲裝置: 中國, 103594080[P]. 2014-02-19.

WU Jian, WANG Gang, BAI Xiaochun. A kind of lightweight membrane-type sound insulation device of metamaterials with broadband gap: China, 103594080[P]. 2014-02-19.

[9] WANG Yanfeng, WANG Yuesheng, SU Xiaoxing. Large bandgaps of two-dimensional phononic crystals with cross-like holes[J]. Appl. Phys, 2011, 110(11): 1135201-11352013.

[10] 劉亞峰, 周蕭明, 胡更開. 彎曲波在薄板超材料中的傳播特性研究[C]//北京力學會第18屆學術年會論文集. 北京: 北京力學會, 2012: 47-48.

LIU Yafeng, ZHOU Xiaoming, HU Gengkai. The research of propagation characteristics of bending wave in thin plate of metamaterials [C]∥18th GAMM in Beijing, 2012: 47-48

[11] Yang Min, Ma Guancong, Yang Zhiyu, et al. Coupled Membranes with Doubly Negative Mass Density and Bulk Modulus[J]. Phys. Rev. Lett, 2013, 110(13): 1342301-1342309.

[12] Yang Z, Dai H M, Chan N H, et al. Acoustic metamaterial panels for sound attenuation in the 50-1000Hz regime[J]. Appl. Phys. Lett, 2012, 96(04): 041906-041909

[13] MA Fuyin, WU Jiuhui, MENG Huang, et al. A purely flexible lightweight membrane-type acoustic metamaterial[J]. Phys. D: Appl. Phys, 2015, 48(17): 175105-175108.

A petal-like acoustic metamaterial structure based on local resonance with ultra-wide sonic band gap in low frequency range

CHEN Lin, WU Wei-guo, ZHOU Rong

(,212013,,)

This paper proposes a petal-like acoustic metamaterial structure based on local resonance to solve the low-frequency acoustic attenuation problem. The band structures and vibration characteristics of the designed material cell are analyzed with COMSOL MULTIPHYSICS 4.3. It is shown that the proposed petal-like metamaterial structure can generate band gap much wider than that of the structure without petals in low frequency range. Furthermore, the structure can produce several complete band gaps by only one kind of cell of a single layer rather than different cells and several layers. The paper also studies the problem on size optimization of the acoustic metamaterial structure.

acoustic metamaterial; local resonance; sonic band gap

O422.7

A

1000-3630(2016)-03-0222-06

10.16300/j.cnki.1000-3630.2016.03.007

2015-11-13;

2016-02-13

鎮江市科技支撐項目(GY2013052)

陳琳(1988－), 女, 甘肅定西人,碩士, 研究方向為聲學超材料、低頻噪聲控制。

吳衛國, E-mail: wuwg@ujs.edu.cn