四方折疊梁聲子晶體低頻帶隙特性研究

高南沙，李沛霖，周文林，侯 宏

（1.西北工業大學 航海學院，西安 710072；2.西北工業大學 聲學信息感知工信部重點實驗室，西安 710072）

彈性波在聲子晶體中傳播時，受其內部周期結構作用，形成特殊色散關系（能帶結構），色散關系曲線之間的頻率范圍稱為帶隙。理論上，帶隙頻率范圍的彈性波傳播被抑制，然而其它頻率范圍（稱為通帶）的彈性波將在色散關系的作用下無能量損耗地傳播。目前比較成熟的聲子晶體帶隙形成機理有兩種，即Bragg散射機理和局域共振機理。一開始被研究的聲子晶體所具有的帶隙都屬于Bragg散射機理[1-3]。若基體為流體，該聲子晶體最低聲波帶隙的中心頻率可由基體聲速c和晶格常數a確定，即c/2a。這種波長與晶格常數成兩倍的關系與Bragg所發現的晶體中X射線衍射行為類似，因此將這類帶隙的形成機理稱為Bragg散射機理[4]。該機理不利于在低頻范圍內的展開應用[5-6]。局域共振機理的出現和發展改變了這種局面[7]。與Bragg散射機理明顯不同，局域共振機理強調單個原胞的運動模式。通過適當的重構，一般是在基體和散射體中間再加入柔軟的包覆層，使原胞所構成的振子在遠低于晶格尺度波長對應的頻率范圍內發生共振，并阻斷正常的波傳播過程[8]。至今為止，在對聲子晶體和聲學超材料的理論研究上已經有了很大的突破，但是針對于實際的工程問題仍然存在著一些急需而難以解決的問題，對于中高頻振動和噪聲問題，雖然Bragg散射型聲子晶體能夠表現出優于傳統減振降噪方法的帶隙特性，但是傳統的減振降噪方法足以處理這些頻段的問題，而且在制造成本嚴格控制或空間有限的場合，Bragg散射型聲子晶體存在著體積過大等一系列缺點，所以在低頻范圍內，從實用性的角度看，Bragg散射型聲子晶體具有很大的局限性。2016年，張佳龍[9-10]提出正八邊形孔狀和雙包覆層局域共振聲子晶體模型，該結構在500 Hz以上的中高頻范圍內存在帶隙，依然沒有很好地解決中低頻的減振降噪問題，最近，Lai Y[11]和Mei J[12]等基于局域共振機理設計了相關的聲學超材料結構，解決了低頻降噪減振問題，但是同時都存在著工藝復雜，帶隙頻率范圍窄以及調節不靈活等問題。

通過分割中心諧振子從而增大自由表面來降低帶隙，基于低頻多重振動耦合機理提出了一種新型二維三組元四方折疊梁聲子晶體結構，通過改變基體和諧振子之間的幾何結構，在不增加單元結構和質量的情況下打開低頻帶隙，具有良好的降噪效果。通過調節單元結構特性參數，得出影響該結構帶隙位置、寬度的主要因素，從而總結出控制該結構低頻減振降噪效果的影響規律。

1 多重振動耦合機理

1.1 二維四方折疊梁聲子晶體模型

二維四方折疊梁聲子晶體的結構如圖1所示，利用四方折疊的彈性梁來替代傳統局域共振型聲子晶體中的包覆層，然后將四方折疊梁分別連接在正方形外框和正方形諧振子上。圖1(c)個給出了所對應Brillouin區。圖1(d)和圖1(e)分別是二維四方折疊梁局域共振型聲子晶體的兩種結構，即模型1和2。可以看出，在這兩種模型中，中間層四方折疊梁是一個互補的對應關系，而圖1(f)則給出了對比結構—模型3。

模型1、2和3中的不同顏色區域A、B和C分別代表不同的材料，四方折疊梁使用金屬鋁，正方形諧振子使用金屬鎢，而正方形外框采用硅橡膠，材料參數如表1所示。

幾何結構參數晶格常數a=34 mm，中心正方形諧振子的邊長H是20 mm，四方折疊梁寬度b=1.43 mm，四方折疊梁之間的空隙寬度e=2.14 mm，內外層開口寬度C1=C2=2.2 mm。

1.2 帶隙特性分析

圖1 二維四方折疊梁聲子晶體基本結構及其對應的Brillouin區

本文基于COMSOL Multiphysics軟件，詳述二維固-固型聲子晶體的能帶結構計算過程。在COMSOL Multiphysics軟件3.5a版本的模型導航窗口中，選擇結構力學中的模塊-特征頻率分析。在進入分析界面后，畫好幾何模型后進行常數設定，為了便于在以后的分析中進行參數掃描，設定初始常數值kx=1，ky=1，在這里kx和ky是倒格子空間的基矢。在求解域的設定中輸入材料的彈性參數。在聲子晶體的概念中，單元之間互相聯通的材料稱為基體，在單元中心質量比較大的材料稱為散射體，又可以稱為諧振子，而它們之間的部分稱為包覆層。在本算例中基體采用環氧樹脂，包覆層采用硅橡膠，而散射體采用鉛。材料參數詳見附錄。根據一般的有限元法步驟確定形函數、建立單元內的剛度矩陣和質量矩陣，獲得單元的廣義本征值方程

表1 材料參數

式中：Ks——單元剛度矩陣；Ms——質量矩陣；u(v)——整體節點位移陣列。

圖2 模型1、2和3的能帶結構

周期結構的聲子晶體，可以采用Block周期性條件將整體的計算縮減到一系列的單元中進行。所以在軟件中取X和Y方向的兩條邊作為源項，它們的表達式中分別輸入單元內的特征向量u和v，而目的端選擇的兩條邊是和源項兩條邊相互對應的。在目的端表達式中，輸入u(v)?exp(i?pi?kx(ky))，這樣就完成了周期性邊界條件的設置。在求解器參數中設置所求特征頻率的階數，然后采用平面應變三角形單元進行有限元網格劃分，完成求解。在軟件中，選擇Eigevalue求解器，即給定一個k值，即可求得一組結構本征頻率，最后進行波矢參數掃描。本節研究對象是聲子晶體具有代表性的正方晶格，所以掃描的邊界應該遍歷整個Brillouin的邊界。當波矢沿著M-Γ-X-M的邊界完成參數掃描后生成txt文檔，導入MATLAB軟件中就可以得出其能帶結構。

圖2、圖3給出了模型1、2和3的能帶結構和模型1的透射譜。在圖2中，灰色區域表示帶隙，可以很清楚地看出，對比結構模型3在1 200 Hz以下沒有帶隙，模型1和2在400 Hz以下都可以打開帶隙，不同的是模型2打開了四條完全帶隙，而模型1只打開兩條，但是模型1的兩條帶隙的寬度比模型2的大，第一帶隙從100 Hz到106.7 Hz，第二帶隙從193 Hz到297帶隙，其第二帶隙寬度達到100 Hz。因為模型1和2在XY方向上屬于不對稱的。所以其所對應的能帶結構中存在著很多方向帶隙[13]，方向帶隙在特定的單一方向上有很多應用，但本論文重點關注的是在全方向上能夠發揮減振降噪作用的完全帶隙。在400 Hz以下的頻率范圍，模型1的帶隙范圍大概為26.5%，與具有相同尺寸和材料的模型3相比，不僅可以打開帶隙，而且帶隙范圍很寬，同時其結構的面密度減小，幾何結構參數變多，能夠對帶隙進行調節。同樣屬于四方折疊梁局域共振型聲子晶體的模型2，其帶隙數量雖然多，但是寬度很窄，最寬的帶隙只能達到5.1 Hz，分布在300 Hz和350 Hz之間。但是在50 Hz以下，甚至在低于10 Hz的范圍內，圖2(b)中有很多平直帶，這說明模型2在超低頻處存在更多的共振作用，這為振動能量回收提供了可行性。在圖3中，模型1的透射譜其最低的峰值低于-250 dB，在190 Hz到300 Hz之間有比較高的隔聲性能，該區間和圖2中的帶隙相互對應，驗證了本文能帶結構計算結構的準確性。

圖3 模型1的透射譜

1.3 物理機理分析

在研究單元的模態位移時取2×2的大原胞作為研究對象，這樣也可以更加清楚地反映出各個單元之間的相互影響。圖4(a)－圖4(d)給出了模型1在第一、二帶隙的上下邊界處的位移振型圖。

首先討論模型1第一帶隙的上下邊界關鍵點A、B的模態，A、B兩點所對應的頻率為100 Hz和106.7 Hz。在A點的模態中，振動幾乎都集中在中心諧振子上面，在X或Y方向上，外邊框呈現出旋轉振動，且有一定的變形。在左上和右下小單元中，外框的變形是順時針的；而在左下和右上小單元中，外框的變形是逆時針的。B點的模態完全不同于A點的模態，這些振動的中心都在每一個小單元的諧振子中心上，且方向是順時針的，大原胞外框的振動是反對稱的彎曲振動。C點是模型1第二帶隙下邊界的點，其模態類似B點的模態，振動中心依然在諧振子上，外框也有明顯地彎曲變形，但是不同的是其中有兩個小單元的旋轉振動方向不同于B點中的旋轉方向。D點，是模型1第二帶隙上邊界的點，可以看出諧振子的振動位移不大，變形主要發生在每一個小單元的右下外邊框處，每一個小單元的上下四個外邊框都是反對稱彎曲振動。為了更加清楚地看出振型方向，圖4(e)畫出了這四個模態示意圖，在2×2的大原胞中，重點關注的是節點處的位移方向（圖4(e)中的黑色虛線框）。

圖4 模型1前兩條帶隙的上下邊界關鍵點處的模態振型圖及其關鍵點的模態示意圖

圖中表明A、C兩點在節點處的模態振動位移類似；而B、D兩點在節點處的模態振動位移類似。在B、D點的模態中，在節點處位移方向是相反的，這樣使得彈性波在每一個節點處不能前行；而在A、C點的模態中，當彈性波傳播到節點時，會向垂直于彈性波振動的方向進行分流，從而使得彈性波無法繼續向前傳播，產生帶隙。其背后的物理機理為：中心諧振子的質量相對很大，導致相鄰諧振子之間的相互作用較大，外邊框的材料非常軟，所以外邊框的振動模式不同于中心諧振子，諧振子（散射體）在第1階Mie散射峰附近呈現出剛體共振模式，這是因為諧振子在振動中充當了“質量”的作用，而外邊框基體材料充當了“彈簧”的作用。在本節的聲子晶體結構中，諧振子的材料是高密度和高剛度的金屬鎢，但是由于四方折疊梁的設計，導致填充率不高。因此，在傳播過程中，縱波將轉化成橫波向前傳播，諧振子的振動模式受到外邊框振動作用和相鄰諧振子之間的相互作用，多種振動模式在單元之間的節點處發生耦合，如圖4(e)所示，所以在這些頻率處，完全帶隙被打開。

從圖2中可以看出，四方折疊梁聲子晶體的帶隙附近的頻散曲線沒有平直帶，因此該結構不屬于典型的局域共振聲子晶體。所以，準確來說，多種振動的耦合，特別是在節點處的耦合，才是打開低頻帶隙的關鍵因素。四方折疊梁的引入，很大程度上減小了整體結構的等效剛度，而等效質量雖然也在減少，但其帶來的影響遠低于剛度的減小。這直接會導致整體結構的共振頻率向下移動，從幾何結構上知道，模型1的四方折疊梁的剛度比模型2低，相對應的帶隙頻率也低于模型2的帶隙頻率。

2 二維四方折疊梁聲子晶體帶隙影響因素

2.1 幾何尺寸對帶隙的影響

本節討論模型1中四方折疊梁的尺寸對于帶隙的影響，其中晶格常數a和中心正方形諧振子的尺寸不變。在圖5(a)中，b e的比值從0.3增加到1.6，第一和第二帶隙寬度不斷減小，當b e等于0.9時，第一帶隙消失，所以減小b的尺寸可以有效地增加低頻帶隙寬度。尺寸b是四方折疊梁的寬度，減小其尺寸可以有效地減小等效剛度，從而使得相關帶隙向低頻移動。從圖5(b)和圖5(c)中，可以看出C1和C2尺寸對于帶隙的影響具有一致性。隨著C1和C2尺寸的不斷增加，第一帶隙向低頻移動，但是第二帶隙的寬度和位置幾乎不變。從圖5(a)中看出，C1、C2尺寸對于整體結構來說非常小，其改變不會對整體結構的振動形式有很大的影響，對于等效質量和等效剛度的影響也可以忽略不計，所以C1、C2對帶隙的影響非常小。

2.2 排列方式對帶隙的影響

從結構上可以看出，二維四方折疊梁聲子晶體（模型1或者2）不是完全對稱的結構，因此不同排列組合的模型一定會對應不同的模態振型，本節取2×2的大原胞作為研究對象，不同的排列組合方式如圖5－圖6(a)－圖6(d)所示，這4種排列方式的通帶、帶隙分布范圍如圖5－圖7所示。A排列的大原胞具有兩條帶隙，其能帶結構和圖5－圖3(a)一樣；B排列依然具有兩條帶隙，但是其第一帶隙的位置更低，第二帶隙的下邊界也更低；C排列在300 Hz以下的區域內全是通帶；D排列的結構仍然可以打開兩條帶隙，但是這兩條帶隙的位置都偏高，其寬度也比較小。可以通過不同的排列組合，來調整帶隙的位置和寬度，這對于不對稱聲子晶體結構來說，是調控帶隙的有效方法。相對于以前的研究來說[14-18]，大多數學者都集中于單個單元的幾何結構尺寸和材料組分的研究。對于很多相類似的結構，通過不同的排列組合打開或者關閉帶隙是一種較簡單的調控帶隙方法。

3 結語

圖5 幾何參數變化對于帶隙的影響

(1)本文提出了一種二維四方折疊梁局域共振型聲子晶體模型，對能帶結構和帶隙關鍵點的模態位移進行分析，詳述了帶隙背后的物理機理。其結構多種振動的耦合，特別是在節點處的耦合作用，是打開低頻帶隙的關鍵。

圖6 二維四方折疊梁聲子晶體不同的四種排列組合形式

圖7 A、B、C和D四種排列方式所對應的低頻帶隙分布

(2)二維四方折疊梁局域共振型聲子晶體的幾何尺寸b，即四方折疊梁的寬度，其增大可以導致帶隙變窄，甚至消失。折疊梁、基體材料和諧振子連接處的幾何尺寸C1和C2對于帶隙的影響具有一致性，其尺寸的增加僅僅使第一帶隙向低頻移動。不同排列組合的2×2的大原胞對應著不同的帶隙變化規律。這種結構設計為制備“低頻，寬帶，強衰減”特性的聲子晶體提供了理論依據和方法指導。