吳飛 黃威 陳文淵 肖勇 郁殿龍 溫激鴻?

1) (國防科技大學(xué)智能科學(xué)學(xué)院,裝備綜合保障技術(shù)重點實驗室,長沙 410073)

2) (西南大學(xué)工程技術(shù)學(xué)院,重慶 400715)

本文設(shè)計了一種由微穿孔板與折曲通道組成的亞波長寬帶復(fù)合吸聲結(jié)構(gòu),對該復(fù)合吸聲器低頻寬帶吸聲機理進行了詳細(xì)分析,建立了該復(fù)合吸聲結(jié)構(gòu)的理論吸聲解析模型與有限元數(shù)值分析模型,完成了吸聲理論解和數(shù)值解的相互驗證.該吸聲結(jié)構(gòu)在整體厚度為60 mm時,理論上實現(xiàn)了在200—500 Hz頻段內(nèi)平均吸聲系數(shù)達0.8的低頻高效吸聲.同時在整體厚度為90 mm時,理論上實現(xiàn)了180—350 Hz頻段內(nèi)多處峰值達0.95的準(zhǔn)完美吸聲.該復(fù)合吸聲結(jié)構(gòu)在低頻噪聲控制工程中具有一定的應(yīng)用前景.

1 引 言

中低頻噪聲具有很強的穿透能力,一直是艙室噪聲污染的主要來源.目前,傳統(tǒng)的吸聲材料如多孔材料[1]已經(jīng)被證明能對中高頻噪聲(>1000 Hz)進行有效地吸收,但在低頻范圍內(nèi)，則需要與波長相當(dāng)?shù)慕Y(jié)構(gòu)厚度，這阻礙了其在實際工程中的應(yīng)用，低頻噪聲的吸收問題亟待解決[2].

近年來,深亞波長聲學(xué)超材料在低頻噪聲(<500 Hz)控制工程中受到了極大的關(guān)注.聲學(xué)超材料 (acoustic metamaterials,AMMs)概念的提出為解決低頻噪聲問題提供了新思路[3?5].目前,大量學(xué)者基于共振吸聲原理設(shè)計了許多深亞波長吸聲材料和結(jié)構(gòu),例如薄膜型共振吸聲器[6?8](decorated membrane resonator,DMRs),亥姆霍茲共振吸聲器[9,10](Helmholtz resonators,HRs)和Fabry-Pérot (FP)通道共振吸聲器[11,12].這些吸聲結(jié)構(gòu)能在低頻范圍內(nèi)展示出良好的吸聲性能，但相對吸聲帶寬較窄。 目前,已有研究人員對基于微穿孔板與折曲通道的低頻寬帶復(fù)合吸聲結(jié)構(gòu)[13?15]進行了研究,其有效吸聲帶寬為 232—533 Hz,但該吸聲結(jié)構(gòu)的整體厚度達到了100 mm.基于微穿縫板與折曲通道的低頻寬帶復(fù)合吸聲器[16,17]厚度較小,但其仍然存在吸聲帶寬較窄的問題,實際應(yīng)用場景受限.

基于此，本文設(shè)計了一種由微穿孔板與多個折曲通道組合而成的復(fù)合吸聲結(jié)構(gòu),能夠有效解決在限制吸聲結(jié)構(gòu)厚度的條件下200—500 Hz噪聲難以吸收的問題.

2 復(fù)合吸聲結(jié)構(gòu)的吸聲理論分析

減小吸聲結(jié)構(gòu)厚度的有效方法是引入空間盤繞或迷宮式結(jié)構(gòu)[18].本文設(shè)計并驗證了一種由微穿孔板(穿孔直徑小于1 mm)和折曲FP通道組合而成的復(fù)合聲學(xué)超材料.該結(jié)構(gòu)在傳統(tǒng)微穿孔板的基礎(chǔ)上,將背腔進行折曲處理,使其可以不增加總厚度而大大增強對低頻噪聲的吸收性能.同時,將共振頻率不同的吸聲單元結(jié)構(gòu)并聯(lián)設(shè)計,使其能在不添加其他吸收性材料(吸聲棉等)[19,20]的情況下獲得較寬的吸收帶寬.

如圖1所示,該復(fù)合吸聲結(jié)構(gòu)由微穿孔板和折曲FP通道組成,入射聲波沿Z軸進入垂直于微穿孔板的通道.為探究其吸聲機理,本文首先建立一個如圖1(d)所示的5通道單元的二維模型,聲波在盤繞通道內(nèi)的傳播路徑如藍(lán)色虛線箭頭所示.其中,微穿孔板的主要結(jié)構(gòu)參數(shù)有孔徑d、板厚t和穿孔率p,周期性通道的主要結(jié)構(gòu)參數(shù)包括長度L、寬度W和高度H,等效吸聲長度D,內(nèi)部隔板厚度t.在理論分析中,折曲FP通道可簡化為串聯(lián)的多個直同軸連接子通道,各子通道內(nèi)的有效傳播長度表示為Leffi(i=1,2,3,···).

具有剛性背板的吸聲系統(tǒng)的吸聲系數(shù)α可以根據(jù)阻抗理論計算公式得出[18]:

Zsn為第n個復(fù)合吸聲結(jié)構(gòu)的表面阻抗,可以按下式計算:

其中ZM和分別表示微穿孔板和單元第一通道入口 處 的 表 面阻 抗.ξ=S0/S1為面 積 修 正 系 數(shù),S0=W×L和S1=W0×L0分別指吸聲單元和單元第一通道的橫截面積.微穿孔板的表面阻抗可通過下式計算:

圖1 復(fù)合聲學(xué)超材料的結(jié)構(gòu)及理論模型 (a) 微穿孔板;(b) 折曲通道;(c) 吸聲結(jié)構(gòu)三維模型;(d) 吸聲結(jié)構(gòu)二維理論模型Fig.1.Structure and theoretical model of composite acoustic metamaterial: (a) The top micro-perforated panel (thickness t ,diameter d ,perforation rate p );(b) multiple coiled FP channels;(c) schematic of the hybrid metamaterial absorber composed of a microperforated panel (MPP) as a top face sheet and coiled-up Fabry–Perot (FP) channels with folding number n;(d) an approximate analytical two dimensional (2 D) model of one unit cell of a space-coiled metamaterial.All the widths of the channels in the YZ plane are L.The height of the channel along the Z axis is H.

γ=指穿孔直徑與黏性邊界層的厚度之比的倍;J1和J0指一階和零階 Bessel函數(shù).第n個通道頂部入口的表面阻抗和通道內(nèi)部聲阻抗可由下式得出:

其中ν=μ/ρ0及其導(dǎo)數(shù)ν′=k/ρ0Cv可由μ,k和Cv(空氣運動黏度、導(dǎo)熱系數(shù)和恒容比熱)計算得出;αk=(k+1/2)π/W和βn=(k+1/2)π/L為 常 數(shù),P0和γ分別為氣壓與比熱值.

3 低頻寬帶吸聲機理分析

為了對復(fù)合吸聲結(jié)構(gòu)的低頻寬帶吸聲機理進行研究,本文進一步使用COMSOL Multiphsics 5.3內(nèi)嵌的聲-熱黏性聲學(xué)耦合模塊建立了該復(fù)合吸聲結(jié)構(gòu)的有限元分析數(shù)值模型,假設(shè)該復(fù)合吸聲結(jié)構(gòu)所有內(nèi)壁都為硬聲學(xué)邊界.建立如圖1(c)所示的由微穿孔板和多個盤繞通道組成的復(fù)合聲學(xué)吸聲結(jié)構(gòu)的有限元分析模型并進行數(shù)值計算.為探究其低頻吸聲機理,在圖2(a)中繪制了頻率為230 Hz的聲波傳入吸聲器時通道內(nèi)部的聲粒子速度(紅色箭頭)分布情況.在微穿孔板與第一通道入口處觀察到較大的粒子速度,可以判斷出,由于微孔處粒子振動速度較大,入射波能量主要是由聲波與小孔之間的較大摩擦而耗散的.圖2(b)中繪制了230 Hz聲波傳入吸聲器時通道內(nèi)部的多切面聲壓分布圖,可以觀察到,此時黑色虛線內(nèi)的兩個通道折疊數(shù)為5的吸聲單元對聲波的耗散作用最為明顯,即折疊通道數(shù)多的單元(等效吸聲長度大)對低頻聲波起主要吸聲作用,而折疊通道數(shù)少的單元(等效吸聲長度小)對相對高頻聲波起主要吸聲作用.

圖2 230 Hz 聲波傳入時吸聲器內(nèi)聲速與聲壓 (a)聲速分布 (m/s);(b)聲壓分布 (Pa)Fig.2.Sound velocity and sound pressure in a sound absorber when a 230 Hz sound wave is introduced: (a) Sound velocity distribution (m/s);(b) sound pressure distribution (Pa).The absorbers are constructed using the coiled-up channel with geometric parameters: H=59mm ,L 1=L2=L3=L4=L5=4.85mm and W=13.6mm and d=t0=t1=1mm.

為了進一步研究所設(shè)計的復(fù)合吸聲結(jié)構(gòu)的吸聲機理,使用圖形法研究了該復(fù)合吸聲結(jié)構(gòu)在復(fù)頻率表面內(nèi)的反射系數(shù)γ的分布[18].一般地,在無損情況下,反射系數(shù)包含共軛零點和極點,如果在某一頻率下達到完美吸聲,則零點將恰好位于實軸上.

圖3(a)中繪制了該復(fù)合吸聲結(jié)構(gòu)在復(fù)頻率表面內(nèi)的 l g|r|2的分布情況,可以看出,200—250 Hz范圍內(nèi)的零點都落在實軸附近,因此,在共振頻率下可獲得近乎完美吸收.而高于250 Hz的零點偏移實軸較遠(yuǎn),并沒有達到完美吸聲.零點和極點之間的距離也可以表征吸聲帶寬,隨著共振頻率的降低,零點和極點之間的距離逐漸減小的同時吸聲帶寬逐漸變窄.

圖3(b)中繪制了復(fù)合吸聲結(jié)構(gòu)吸聲系數(shù)的理論解、數(shù)值解及各主要單元吸聲性能的理論解,圖中黑色虛線從左到右依次為等效吸聲長度逐漸減小的單元的吸聲系數(shù)理論解,根據(jù)這一規(guī)律可通過調(diào)節(jié)折曲 FP通道的等效吸聲長度來實現(xiàn)期望頻率聲波的吸收.同時發(fā)現(xiàn),通過調(diào)整等效吸聲長度可以容易地調(diào)節(jié)吸聲系數(shù)和吸聲帶寬.圖3(c)展示了相對阻抗的實部Re(ZS/Z0)和虛部Im(Zs/Z0)分布情況.當(dāng)阻抗匹配時,可實現(xiàn)完美吸聲,這意味著相對阻抗的虛部等于0,而實部等于1.圖中虛部曲線在330 Hz處過零點,表明在圖3(b)中的吸收峰存在共振狀態(tài).相對阻抗的實部不嚴(yán)格等于1,吸收系數(shù)略小于100%.

理論解與數(shù)值解變化趨勢基本相同,但由于理論模型由一維理論公式推導(dǎo)得來而數(shù)值分析模型為三維模型,故二者有一定誤差.相對帶寬達到143%(相對帶寬指吸聲系數(shù)最大值的一半處的整個帶寬與共振頻率的比值).近完美吸聲的峰值出現(xiàn)在 260,390 和 470 Hz.吸收峰隨著折疊數(shù)n的增加而明顯降低到更低的頻率,在共振頻率230 Hz處,該復(fù)合吸聲結(jié)構(gòu)厚度僅為工作波長的1/28,具有亞波長吸聲特性.

圖3 復(fù)合吸聲結(jié)構(gòu)的吸聲特性 (a) 反射系數(shù)零、極點分布;(b)吸聲性能曲線;(c)聲阻抗的實部與虛部Fig.3.Sound absorption characteristics of composite sound absorption structure: (a) Zero and pole distributions of reflection coefficients on complex frequency plane;(b) sound absorption performance curve.The dotted line is the theoretical solution and the solid line is the numerical solution,the black dotted line is the theoretical solution of each unit(Gradually increase the equivalent sound absorption length from left to right);(c) real and imaginary parts of relative impedance.

4 復(fù)合吸聲結(jié)構(gòu)的準(zhǔn)完美吸聲

為了進一步提高該復(fù)合吸聲結(jié)構(gòu)的吸聲性能,對其結(jié)構(gòu)參數(shù)進行調(diào)整,調(diào)整后吸聲結(jié)構(gòu)總厚度為90 mm,具體結(jié)構(gòu)參數(shù)見表1,進一步研究其低頻吸聲特性.同樣地,運用 MATLAB,COMSOL軟件對該吸聲結(jié)構(gòu)的吸聲系數(shù)進行了理論和數(shù)值計算分析.圖4(a)中繪制了參數(shù)調(diào)整后吸聲結(jié)構(gòu)在復(fù)頻率表面內(nèi)反射系數(shù)的零點、極點分布情況,可以觀察到零點整體更加靠近0線,吸聲系數(shù)有所增大,且工作頻率相較于圖3(a)有所下降,表明調(diào)整后該吸聲結(jié)構(gòu)的低頻吸聲性能有明顯提高.

圖4 復(fù)合吸聲結(jié)構(gòu)的準(zhǔn)完美吸聲特性 (a) 反射系數(shù)零、極點分布;(b)吸聲性能曲線;(c)聲阻抗的實部與虛部Fig.4.Quasi-perfect sound absorption characteristics of composite sound absorption structure: (a) Zero and pole distributions of reflection coefficients on complex frequency plane;(b) sound absorption performance curve,The dotted line is the theoretical solution and the solid line is the numerical solution;(c) real and imaginary parts of relative impedance.

圖4(b)中繪制了該吸聲結(jié)構(gòu)的理論解與有限元分析結(jié)果.可以看出,理論解與有限元分析結(jié)果變化趨勢基本保持一致,在170—340 Hz的頻率范圍內(nèi)理論平均吸聲系數(shù)大于0.85.相對帶寬為121%,相較于圖3(b)中的結(jié)果,有效工作頻率范圍明顯向低頻移動,且在 185,235 和 315 Hz 等頻響范圍附近的吸聲系數(shù)均大于0.95.圖4(c)為相對阻抗的實部與虛部分布情況,其中虛部曲線在250 和 280 Hz處過零點,而實部曲線在 185 和315 Hz處達到1,與圖3(b)中吸收峰值所對應(yīng)聲波頻率一致,達到準(zhǔn)完美吸聲.

表1 準(zhǔn)完美吸聲結(jié)構(gòu)各吸聲單元結(jié)構(gòu)參數(shù)Table 1.Structural parameters of each sound absorption unit of quasi-perfect sound absorption structure.

5 結(jié) 論

本文設(shè)計了一種由微穿孔板與多個折曲通道組合而成的復(fù)合吸聲結(jié)構(gòu),對該復(fù)合吸聲器低頻寬帶吸聲機理進行了詳細(xì)分析,并據(jù)此建立了理論解析模型與有限元數(shù)值分析模型,完成了理論解和數(shù)值解的相互驗證.該復(fù)合吸聲結(jié)構(gòu)的吸聲機理主要是聲波與微穿孔板開孔處內(nèi)壁的摩擦損耗.據(jù)此,本文并聯(lián)設(shè)計多個不同等效長度的復(fù)合吸聲結(jié)構(gòu),并分析了其聲波吸收的高效性及可調(diào)吸收特性.結(jié)果表明,該復(fù)合吸聲結(jié)構(gòu)與傳統(tǒng)吸聲結(jié)構(gòu)相比,設(shè)計結(jié)構(gòu)簡單,具有深亞波長特性,對低頻噪聲(小于500 Hz)的吸收性能顯著提高,且吸聲頻帶更寬,在噪聲控制工程當(dāng)中具有很大的應(yīng)用潛力.實現(xiàn)了以下工程設(shè)計目標(biāo): 在復(fù)合吸聲結(jié)構(gòu)整體厚度為 60 mm 時,實現(xiàn)了在 200—500 Hz 頻段內(nèi),理論平均吸聲系數(shù)達到0.8的低頻高效吸聲.在復(fù)合吸聲結(jié)構(gòu)整體厚度為90 mm時,實現(xiàn)了在180—350 Hz頻段內(nèi),多處理論吸聲峰值達到0.95的準(zhǔn)完美吸聲.