田文昊,吳錦武,李 威,韓 偉
(1.南昌航空大學(xué) 飛行器工程學(xué)院,南昌 330063;2.江蘇航空職業(yè)技術(shù)學(xué)院 航空工程學(xué)院,江蘇 鎮(zhèn)江 212134)
我國(guó)著名聲學(xué)專家馬大猷院士在20世紀(jì)60年代提出了微穿孔板吸聲結(jié)構(gòu)理論,該結(jié)構(gòu)只需要在原來(lái)面板上進(jìn)行微穿孔就能產(chǎn)生與大氣聲阻相匹配的聲阻,同時(shí)聲質(zhì)量也較小,不需要添加任何聲阻尼材料就可以對(duì)聲能進(jìn)行消耗,具有清潔、無(wú)污染和耐高溫等優(yōu)點(diǎn),既可以用于航空發(fā)動(dòng)機(jī)燃燒室,又可以用于飛機(jī)客艙進(jìn)行減振降噪。
單層微穿孔板結(jié)構(gòu)缺點(diǎn)在于其吸聲頻帶較窄,因此設(shè)計(jì)雙層串聯(lián)微穿孔板結(jié)構(gòu),可在拓寬頻帶同時(shí)獲得較高的吸聲系數(shù),本文以此為出發(fā)點(diǎn),通過(guò)理論和仿真設(shè)計(jì)合適的參數(shù),從而得到較優(yōu)的吸聲效果。
根據(jù)馬大猷院士提出的微穿孔板吸聲體的準(zhǔn)確理論,單層微穿孔板吸聲體的聲阻抗率除以空氣的特性阻抗可得到微穿孔板吸聲體的相對(duì)聲阻抗為
式中:r為微穿孔板的相對(duì)聲阻,ω為角頻率,f為入射聲音的頻率,m為相對(duì)聲質(zhì)量,D為微穿孔板后空腔厚度,ρ為空氣密度,c為聲速,η為黏滯系數(shù),t、d和p分別為微穿孔板的厚度、穿孔直徑和穿孔率,k為穿孔常數(shù)。其中t、d單位為mm,p為百分?jǐn)?shù)。式中:D單位為m,其他均為標(biāo)準(zhǔn)單位。則聲波垂直入射時(shí)微穿孔板結(jié)構(gòu)的吸聲系數(shù)公式為
圖1是雙層串聯(lián)微穿孔板結(jié)構(gòu)示意圖,其中t、d和代表參數(shù)和前文相同。
圖1 雙層串聯(lián)微穿孔板結(jié)構(gòu)示意圖
用傳遞矩陣法計(jì)算結(jié)構(gòu)的聲阻抗,其中微穿孔板和空腔的傳遞矩陣分別為
其中:Zs為微穿孔板的聲阻抗率,除以ρc即為微穿孔板相對(duì)聲阻抗z=r+jωm,計(jì)算公式參考文獻(xiàn)[3] 。對(duì)于雙層串聯(lián)微穿孔板結(jié)構(gòu),將微穿孔板的傳遞矩陣[P] 和空腔的傳遞矩陣[S] 按順序連乘,得到總的傳遞矩陣
根據(jù)文獻(xiàn)[4],建立微穿孔板左側(cè)與剛性壁面的聲壓與振動(dòng)速度的矩陣關(guān)系為
其中:p1、u1分別為微穿孔板左側(cè)的聲壓與振動(dòng)速度,p3、u3分別為剛性壁面處的聲壓與振動(dòng)速度。由于振動(dòng)速度在剛性壁面為零,即u3=0,則總聲阻抗率為
展開(kāi)T總得到T11及T21,代入上式可得Zs,再除以空氣的特性阻抗ρc,即為雙層串聯(lián)微穿孔板結(jié)構(gòu)的相對(duì)聲阻抗
則聲波垂直入射雙層串聯(lián)微穿孔板結(jié)構(gòu)的吸聲
系數(shù)
以上為吸聲系數(shù)理論計(jì)算公式,用MATLAB編程繪圖便可得到理論吸聲系數(shù)曲線圖,與COMSOL的仿真結(jié)果對(duì)比驗(yàn)證其可靠性,再深入分析其參數(shù)對(duì)整體吸聲系數(shù)的影響。
首先對(duì)單層微穿孔板結(jié)構(gòu)和雙層串聯(lián)微穿孔板結(jié)構(gòu)進(jìn)行研究對(duì)比,其吸聲系數(shù)曲線圖見(jiàn)圖2。
圖2以及后面圖中各穿孔板的具體參數(shù)詳見(jiàn)表1。
表1 各微穿孔板參數(shù)表
由圖2可以得知理論計(jì)算曲線和仿真曲線基本重合,證明理論計(jì)算的準(zhǔn)確性,其次雙層串聯(lián)微穿孔板結(jié)構(gòu)與單層微穿孔板結(jié)構(gòu)相比,其整體的吸聲頻帶得到拓寬,吸聲系數(shù)也有明顯增大。
圖2 單層穿孔板和雙層串聯(lián)微穿孔板吸聲系數(shù)曲線圖
進(jìn)而改變穿孔參數(shù)使2層微穿孔板的穿孔情況不同,其吸聲系數(shù)曲線見(jiàn)圖3。
圖3 不同穿孔直徑和穿孔率對(duì)吸聲系數(shù)的影響曲線圖
根據(jù)圖3可知,當(dāng)串聯(lián)的雙層微穿孔板穿孔參數(shù)不同時(shí),整體的吸聲系數(shù)有明顯增大,其中穿孔直徑的影響比穿孔率的影響大。
將雙層串聯(lián)微穿孔板的2層穿孔板進(jìn)行位置互換,其吸聲系數(shù)曲線如下。
圖4 對(duì)調(diào)兩層穿孔板對(duì)吸聲系數(shù)的影響曲線圖
對(duì)比圖3和圖4可以得出,當(dāng)?shù)?層穿孔板的穿孔直徑和穿孔率都大于第2層時(shí),才會(huì)有同時(shí)拓寬吸聲頻帶和提高吸聲系數(shù)的效果。
然后成比例改變穿孔直徑和穿孔率,其吸聲系數(shù)曲線見(jiàn)圖5。
由圖5可知,當(dāng)?shù)?層穿孔板和第2層穿孔板的穿孔參數(shù)成一定比例時(shí),減小第2層的穿孔直徑和穿孔率可以引起第2個(gè)吸聲峰值凸起并增大吸聲系數(shù)。
進(jìn)一步在板后空腔總體厚度不變的情況下改變2層穿孔板的相對(duì)位置,其吸聲系數(shù)曲線見(jiàn)圖6。
圖5 按比例減小穿孔直徑和穿孔率對(duì)吸聲系數(shù)的影響曲線圖
圖6結(jié)果表明,在穿孔板腔體總體厚度及穿孔參數(shù)不變的情況下,第1層穿孔板和第2層穿孔板后空腔厚度比為1:1~1:3時(shí),結(jié)果較優(yōu),兩者相比,前者總體吸聲系數(shù)較大但吸聲頻帶寬度不如后者。
圖6 穿孔板相對(duì)位置對(duì)吸聲系數(shù)的影響曲線圖
在板后空腔總體厚度不同條件下,其吸聲系數(shù)曲線見(jiàn)圖7。
圖7 板后空腔總體厚度對(duì)吸聲系數(shù)的影響曲線圖
結(jié)果表明,第1層、第2層穿孔板板后空腔之比為1:1時(shí),增加板后空腔厚度,吸聲系數(shù)大小變化不大,但整體吸聲系數(shù)頻帶變窄并向低頻移動(dòng)。
使用傳遞矩陣法得到了雙層串聯(lián)微穿孔板機(jī)構(gòu)的相對(duì)聲阻抗從而得出其吸聲系數(shù)曲線,通過(guò)軟件仿真驗(yàn)證了其準(zhǔn)確性,再改變穿孔參數(shù)得出當(dāng)?shù)?層穿孔板的穿孔直徑和穿孔率比第2層大且第2層穿孔板的穿孔直徑和穿孔率都較小時(shí),能夠一定程度上拓寬吸聲頻帶,提高整體的吸聲系數(shù);當(dāng)板后空腔厚度不變,第1層穿孔板的穿孔直徑和穿孔率接近第2層2倍且第2層穿孔板的穿孔直徑達(dá)到0.2 mm,穿孔率接近2%,2層板后空腔厚度之比為1:1到1:3時(shí)(1:1時(shí)吸聲頻帶內(nèi)吸聲系數(shù)最大),2個(gè)吸聲波峰之間幾乎無(wú)波谷出現(xiàn),吸聲頻帶內(nèi)的整體吸聲系數(shù)較高;當(dāng)板后空腔總體厚度增加且兩層板后空腔厚度比例不變時(shí),吸聲系數(shù)峰值不變,但吸聲頻帶向低頻移動(dòng),吸聲帶寬變窄。