考慮微穿孔管消聲器結(jié)構(gòu)參數(shù)的共振頻率預(yù)估模型

第一作者左曙光男，博士，教授，博士生導(dǎo)師，1968年生

考慮微穿孔管消聲器結(jié)構(gòu)參數(shù)的共振頻率預(yù)估模型

左曙光，龍國，吳旭東，相龍洋，張孟浩，胡佳杰(同濟大學(xué) 新能源汽車工程中心，上海201804)

摘要：引入微穿孔板聲阻抗理論模型，通過微穿孔管消聲器傳聲損失測試驗證基于該模型的傳聲損失數(shù)值計算方法的正確性。利用數(shù)值方法計算微穿孔管消聲器傳聲損失獲得共振頻率。單因素法分析穿孔段長度、穿孔率、穿孔直徑和空氣腔厚度等結(jié)構(gòu)參數(shù)變化對共振頻率影響發(fā)現(xiàn)，穿孔段長度及空氣腔厚度對共振頻率影響顯著。利用均勻設(shè)計結(jié)合回歸分析法所得共振頻率預(yù)估模型，能直接反映微穿孔管消聲器共振頻率與結(jié)構(gòu)參數(shù)之關(guān)系，對微穿孔管消聲器優(yōu)化設(shè)計具有指導(dǎo)意義。

關(guān)鍵詞：微穿孔管消聲器；聲阻抗模型；共振頻率；單因素分析；預(yù)估模型

基金項目：國家重大儀器設(shè)備專項(2012YQ150256)

收稿日期：2014-02-26修改稿收到日期：2014-04-24

中圖分類號:TB535文獻標(biāo)志碼：A

Resonant frequency predicting model for a micro-perforated tube muffler considering its structural parameters

ZUOShu-guang,LONGGuo,WUXu-dong,XIANGLong-yang,ZHANGMeng-hao,HUJia-jie(Clean Energy Automotive Engineering Center, Tongji University, Shanghai 201804, China)

Abstract:The resonant frequency of a micro-perforated tube muffler (MPTM) is closely related to its parameters, such as, perforation segment length, perforation rate, perforation diameter and air-cavity thickness etc. Here, a theoretical acoustic impedance model of a mirco-perrorated plate (MPP) was introduced firstly and a measurement for transmission loss (TL) of a micro-perforated tube muffler was conducted to validate the correctness of the TL numerical calculation method based on the model introduced above. Then the MPTM’s resonant frequency was obtained from the TL results calculated with the numerical calculation method. The effects of the MPTM’s structural parameters on its resonant frequency were analyzed with the single factor analysis method. The results showed that perforation segment length and air-cavity thickness have significant effects on a MPTM’s resonant frequency. Moreover, a resonant frequency predicting model for a MPTM was derived with the uniform design method and the regression analysis, it could reflect directly the relationship between a MPTM’s resonant frequency and its structural parameters. The predicting model provided a guide for the optimal design of MPTMs.

Key words:micro-perforated tube muffler (MPTM); acoustic impedance model; resonant frequency; single factor analysis; predicting model

微穿孔管消聲器為基于微穿孔板吸聲理論[1]用于管道消聲的共振式消聲器。由穿孔管及其外部空氣腔組成的腔體類似多個亥姆霍茲共振器并聯(lián)，具有良好的寬頻消聲效果[2-3]。聲波在微穿孔管消聲器內(nèi)激起共振時聲能損耗最大，消聲效果最顯著。由于微穿孔管消聲器共振頻率附近消聲帶較窄，共振峰偏移易導(dǎo)致消聲量急劇下降，準(zhǔn)確確定共振頻率至關(guān)重要。馬大猷最先分析獲得微穿孔板吸聲結(jié)構(gòu)的共振頻率，但該隱式復(fù)雜表達式不能直接求值。而微穿孔板吸聲結(jié)構(gòu)共振頻率亦不能反映集阻抗復(fù)合消聲特性于一體的微穿孔管消聲器共振頻率。Selamet等[4]利用二維分析方法建立穿孔管消聲器傳聲損失的理論模型獲得共振頻率，并定性分析消聲器參數(shù)對共振頻率影響。高林[5]利用有限元數(shù)值仿真計算亥姆霍茲共振式消聲器共振頻率，發(fā)現(xiàn)與實驗存在一定誤差并分析誤差產(chǎn)生原因。畢嶸[6]將集中參數(shù)法、理論及數(shù)值仿真法所得亥姆霍茲共振式消聲器共振頻率與實驗結(jié)果對比發(fā)現(xiàn)，數(shù)值計算結(jié)果更接近實驗值，相對誤差約3%。孟新等[7]通過求解亥姆霍茲共振消聲器的氣流再生噪聲獲得共振頻率，但缺少實驗證明。羅虹等[8-9]利用有限元方法計算穿孔管消聲器的傳聲損失并分析歸納穿孔管結(jié)構(gòu)參數(shù)對共振頻率影響。由此可見，共振式消聲器共振頻率求解大多基于頸部直徑較大的單腔亥姆霍茲共振器模型，利用集中參數(shù)法或理論方法，所求結(jié)果與實驗結(jié)果偏差較大，數(shù)值方法則更繁瑣。而針對微穿孔管消聲器共振頻率研究極少，仍未獲得較準(zhǔn)確反映共振頻率與結(jié)構(gòu)參數(shù)關(guān)系計算式。

本文用數(shù)值計算方法求解微穿孔管消聲器傳聲損失獲得共振頻率，分析主要結(jié)構(gòu)參數(shù)對共振頻率影響。引入微穿孔板聲阻抗理論模型，通過實驗驗證基于該模型的微穿孔管消聲器傳聲損失數(shù)值計算方法。通過單因素法分析結(jié)構(gòu)參數(shù)對共振頻率影響。基于均勻設(shè)計結(jié)合回歸分析方法獲得微穿孔管消聲器共振頻率預(yù)估模型。

1微穿孔板聲阻抗理論模型

圖1為考慮中低聲強下基于短管(管長遠小于聲波長)的聲傳播模型，其短管聲阻抗[10]為

(1)

(2)

圖1　考慮管壁粘滯作用的 短管內(nèi)媒質(zhì)運動模型 Fig.1 Model of medium motion in short tube considering viscous effect

式中：j為虛數(shù)單位；ω為聲波波動圓頻率；ρ為流體介質(zhì)密度；t為短管長度；d為短管直徑；J0，J1分別為0、1階第一類貝塞爾函數(shù)；η為流體切變粘滯系數(shù)。

當(dāng)管長t不比管徑d大的多時，需加末端修正值[11]，則修正后的短管聲阻抗為

(3)

式中：Δl0為末端修正等效長度，即

Δl0=0.85d

(4)

考慮微穿孔板小孔間距大于孔徑，忽略穿孔間相互影響。單個小孔(視為短管)聲阻抗除以穿孔率及空氣特性阻抗得微穿孔板相對聲阻抗率為

(5)

式中：p為微穿孔板穿孔率；c為空氣聲速；ρc為空氣特性阻抗。

式(3)～式(5)構(gòu)成微穿孔板聲阻抗理論模型。需指出的是，在計算機技術(shù)尚不能滿足計算要求年代，文獻[1]為方便制作參數(shù)設(shè)計表格，對模型中貝塞爾函數(shù)近似簡化。現(xiàn)可借助科學(xué)軟件直接計算貝塞爾函數(shù)，無需簡化。

2數(shù)值計算方法的實驗驗證

本文研究的微穿孔管消聲器各結(jié)構(gòu)參數(shù)表示見圖2(穿孔段長度l、空氣腔厚度D、前后內(nèi)插管長度l1及l(fā)2、內(nèi)管管徑D0分別對應(yīng)圖2標(biāo)示；其它參數(shù)包括穿孔孔徑d、穿孔層厚度t、穿孔率p)。此處簡要介紹其傳聲損失的數(shù)值計算方法[12]。

圖2　微穿孔管消聲器示意圖 Fig.2 Schematic diagram of a micro-perforated tube muffler

利用有限元軟件對微穿孔管消聲器傳聲損失進行數(shù)值求解。有限元仿真模型見圖3。實體建模時忽略消聲器內(nèi)部聲傳播過程的流固耦合作用，僅建立消聲器內(nèi)部空氣腔模型。建模后在有限元軟件中分別定義入口激勵、介質(zhì)屬性、出口無反射邊界條件、小孔傳遞導(dǎo)納關(guān)系[13]、求解器，仿真并進行后處理獲得傳聲損失。仿真過程不計溫度、氣流對聲傳播影響；小孔傳遞導(dǎo)納矩陣須結(jié)合微穿孔板聲阻抗理論模型編程求得。

圖3　微穿孔管消聲器有限元模型示意圖 Fig.3 Schematic diagram of a micro-perforated tube muffler FEM model

為驗證基于微穿孔板聲阻抗理論模型的消聲器傳聲損失有限元數(shù)值計算方法，據(jù)雙負(fù)載法[14]利用阻抗管測量消聲器傳聲損失。實驗布置見圖4。所用阻抗管型號為SW466-60，傳聲器間距45 mm，頻率測試范圍400～3 150 Hz。測量樣件結(jié)構(gòu)參數(shù)為：D0=60 mm、D=19 mm、l=220 mm、d=0.8 mm、t=1 mm、p=2%、l1=l2=10 mm。給定聲激勵為白噪聲信號，為忽略聲阻抗非線性效應(yīng)[15],實驗中保證各傳聲器測得聲壓值低于100 dB。所有測量均在相同條件下重復(fù)三次，取結(jié)果平均值。將微穿孔管消聲器傳聲損失測量結(jié)果與傳聲損失有限元數(shù)值計算結(jié)果對比，見圖5。由圖5看出，微穿孔管消聲器傳聲損失數(shù)值計算結(jié)果與實驗數(shù)據(jù)吻合較好，在消聲器共振頻率處基本重合(相差20 Hz，相對誤差0.8%)，能較準(zhǔn)確反映實際消聲器的消聲性能。故可利用該數(shù)值方法仿真計算微穿孔管消聲器傳聲損失以確定其共振頻率。

圖4　微穿孔管消聲器傳聲損失測量實驗布置 Fig.4 Experimental arrangement of TL measurement

圖5　微穿孔管消聲器傳聲損失實驗值與仿真值對比 Fig.5 Comparison between experiment and FEM

3微穿孔管消聲器結(jié)構(gòu)參數(shù)對共振頻率影響分析

據(jù)文獻[9]，穿孔管消聲器共振頻率與穿孔段長度、穿孔率、穿孔直徑、空氣腔厚度、穿孔層厚度等參數(shù)聯(lián)系密切，需利用單因素法分析各結(jié)構(gòu)參數(shù)對微穿孔管消聲器共振頻率影響。微穿孔管消聲器共振頻率即傳聲損失達到峰值時頻率值記為fmax，見圖6。

圖6　某微穿孔管消聲器傳聲損失 Fig.6 TL curve of a micro-perforated tube muffler

為簡化數(shù)值建模及分析，本文不考慮內(nèi)插管長度、穿孔層厚度及內(nèi)管徑對共振頻率影響，即令仿真模型前后內(nèi)插管長度l1=l2=0、穿孔層厚度t=0.5 mm、內(nèi)管管徑D0=50 mm保持不變，下同。仿真計算時頻率步長取20 Hz。

3.1穿孔段長度對共振頻率影響

令穿孔率p=2.2%、穿孔直徑d=0.5 mm、空氣腔厚度D=20 mm保持不變，穿孔段長度l(mm)分別取90、96、102、108、114、120，分析其變化對共振頻率影響。共振頻率隨穿孔段長度變化趨勢見圖7。由圖7看出，穿孔段長度由90 mm增加到120 mm時,共振頻率近似線性降低且變化明顯。故可通過調(diào)整微穿孔管消聲器的穿孔段長度有效改變共振頻率。

圖7　不同穿孔段長度對應(yīng)的共振頻率 Fig.7 Resonant frequencies of micro-perforated tube mufflers with different perforation segment lengths

3.2穿孔率對共振頻率影響

令穿孔段長度l=102 mm、穿孔直徑d=0.5 mm、空氣腔厚度D=20 mm保持不變，穿孔率p分別取1.2、1.5、1.8、2.1、2.4、2.7，分析其變化對共振頻率的影響。共振頻率隨穿孔率變化趨勢見圖8。由圖8看出，穿孔率由1.2%增加到2.7%時,共振頻率變化較復(fù)雜，但整體變化范圍較小(在40 Hz幅度變化)，可認(rèn)為穿孔率在1.2%～2.7%范圍變化對共振頻率影響有限。考慮傳統(tǒng)微穿孔板的穿孔率一般在3%以下，故不宜通過改變微穿孔管消聲器穿孔率改變共振頻率。

圖8　不同穿孔率對應(yīng)的共振頻率 Fig.8 Resonant frequencies of micro-perforated tube mufflers with different perforations

3.3穿孔直徑對共振頻率影響

令穿孔段長度l=102 mm、穿孔率p=2.2%、空氣腔厚度D=20 mm保持不變，穿孔直徑d(mm)分別取0.2、0.3、0.4、0.5、0.6、0.7，分析其變化對共振頻率影響。共振頻率隨穿孔直徑變化趨勢見圖9。由圖9看出，穿孔直徑由0.2 mm增大到0.7 mm時,共振頻率呈現(xiàn)一定波動，但幅度較小(約40 Hz)。考慮傳統(tǒng)微穿孔板穿孔直徑一般不超過1 mm，故不宜通過改變微穿孔管消聲器的穿孔直徑改變共振頻率。

圖9　不同穿孔直徑對應(yīng)的共振頻率 Fig.9 Resonant frequencies of micro-perforated tube mufflers with different perforation diameters

3.4空氣腔厚度對共振頻率影響

令穿孔段長度l=102 mm、穿孔率p=2.2%、穿孔直徑d=0.5 mm保持不變，空氣腔厚度D(mm)分別取10、15、20、25、30、35，分析其變化對共振頻率影響。共振頻率隨空氣腔厚度變化趨勢見圖10。由圖10看出，空氣腔厚度由10 mm增加到35 mm時,共振頻率逐漸降低且降幅較大，且空氣腔厚度越大共振頻率變化幅度越小。故當(dāng)微穿孔管消聲器的空氣腔厚度較小時，調(diào)整其值可顯著、有效改變共振頻率。

圖10　不同空氣腔厚度對應(yīng)的共振頻率 Fig.10 Resonant frequencies of micro-perforated tube mufflers with different air-cavity thickness

4微穿孔管消聲器共振頻率預(yù)估模型

利用均勻設(shè)計試驗方法設(shè)計微穿孔管消聲器的有限元數(shù)值試驗?zāi)Ｐ筒⒎抡媲蠼鈧髀晸p失，通過回歸分析方法獲得考慮結(jié)構(gòu)參數(shù)的共振頻率預(yù)估模型。

4.1均勻設(shè)計

均勻設(shè)計法[16]作為科學(xué)試驗方法，將分析問題化為有同樣解答的概率問題，用統(tǒng)計模擬方法處理。均勻設(shè)計法特點為：①試驗點均勻分散，具有較好代表性；②試驗次數(shù)較正交設(shè)計等試驗設(shè)計方法明顯減少，尤其適合多因素多水平及系統(tǒng)模型完全未知情況；③通過回歸分析、最優(yōu)化及關(guān)聯(lián)度分析等方法處理試驗結(jié)果能估計出模型因素的主效應(yīng)與交互效應(yīng)。均勻設(shè)計作為優(yōu)于正交設(shè)計的新方法，已廣泛應(yīng)用于化工、醫(yī)藥、軍事等行業(yè)。

表1　均勻設(shè)計試驗方案

影響微穿孔管消聲器共振頻率的因素較多且無法確定各因素間有無交互作用。本文利用均勻設(shè)計方法，選U15(158)均勻表并結(jié)合其使用表設(shè)計了15組數(shù)值試驗分別求解傳聲損失獲得各組共振頻率，見表1。微穿孔管消聲器的內(nèi)管管徑均為50 mm，穿孔層厚度均為0.5 mm。

4.2多因素回歸分析

對試驗所得共振頻率，見表2。用回歸分析方法[17]進行處理，獲得共振頻率與消聲器主要結(jié)構(gòu)參數(shù)的回歸模型。利用數(shù)學(xué)統(tǒng)計分析軟件對表2數(shù)據(jù)進行二次多項式回歸分析，所得回歸方程為

fmax=379+84.6X1+5544X2-633X4-

80.5X1X2-15.8X1X3+4.61X1X4+2052X2X3+

70.9X2X4-52.1X3X4-1567X32

式中：fmax，X1，X2，X3，X4分別對應(yīng)共振頻率(Hz)、穿孔段長度(mm)、穿孔率(%)、穿孔直徑(mm)、空氣腔厚度(mm)取值。

表2　試驗所得共振頻率f max

圖11為共振頻率仿真試驗值與回歸預(yù)估值對比。由圖11可知，共振頻率試驗值與回歸擬合值吻合良好，平均誤差3.76%。所得回歸模型的復(fù)相關(guān)系數(shù)R=98.7%，修正決定系數(shù)Ra=95.5%，方差分析及回歸系數(shù)顯著性檢驗分別見表3、表4。由量表看出，回歸方程在α=0.01上能反映出較好的顯著度及預(yù)估精度。故認(rèn)為微穿孔管消聲器共振頻率參數(shù)回歸模型滿足預(yù)估要求。

圖11　共振頻率的仿真試驗值與回歸預(yù)估值 Fig.11 Resonant frequencies calculated by FEM and prediction model

來源自由度方差SS均方差MS方差與均方差比值F概率P值回歸 1012070162120701630.610.002殘差誤差415773139433合計 1412180343

表4　回歸系數(shù)的顯著性檢驗

為進一步檢驗回歸方程的精確度，本文另取均勻設(shè)計試驗?zāi)Ｐ鸵酝獾奈⒋┛坠芟暺饔邢拊Ｐ筒⒂嬎闫鋫髀晸p失獲得共振頻率值。計算結(jié)果及回歸方程預(yù)估結(jié)果見表5。由表5知，誤差不超過1%，考慮仿真計算的頻率步長20 Hz會導(dǎo)致一定誤差，故認(rèn)為回歸方程預(yù)估結(jié)果可以接受。

表5　共振頻率回歸預(yù)估值檢驗

在已知內(nèi)管管徑50 mm(受消聲管路接口限制)及穿孔層厚度0.5 mm(材料限制)前提下，優(yōu)化設(shè)計無內(nèi)插管單層直通微穿孔管消聲器，可據(jù)經(jīng)驗先確定微穿孔管消聲器的穿孔率及穿孔直徑(二者對共振頻率影響相對較小)，基于共振頻率預(yù)估模型，據(jù)目標(biāo)消聲頻率及降噪量并綜合考慮實際情況(消聲器空間布置及加工難易等因素)，調(diào)整穿孔段長度及空氣腔厚度參數(shù)即可。

5結(jié)論

本文通過分析微穿孔管消聲器主要結(jié)構(gòu)參數(shù)對共振頻率影響，結(jié)論如下：

(1)引入較準(zhǔn)確的微穿孔板聲阻抗理論模型，并通過實驗驗證基于該聲阻抗理論模型的消聲器傳聲損失有限元數(shù)值計算方法。

(2)在一定范圍內(nèi)改變穿孔率及穿孔直徑對共振頻率影響不大，而改變穿孔段長度及空氣腔厚度則可顯著改變共振頻率。增大穿孔段長度使共振頻率近似線性降低；增大空氣腔厚度使共振頻率減小。

(3)所得考慮微穿孔管消聲器結(jié)構(gòu)參數(shù)的共振頻率預(yù)估模型，一定程度上可指導(dǎo)微穿孔管消聲器的優(yōu)化設(shè)計。

參考文獻

[1]馬大猷. 微穿孔板吸聲結(jié)構(gòu)的理論和設(shè)計[J]. 中國科學(xué), 1975, 18(1): 38-50.

MA Da-you. Theory and design of micro-perforated panel absorber[J]. Science China, 1975, 18(1): 38-50.

[2]阮登芳. 共振式進氣消聲器設(shè)計理論及其應(yīng)用研究[D]. 重慶: 重慶大學(xué), 2005.

[3]左言言,周晉花,劉海波,等. 穿孔板吸聲結(jié)構(gòu)的吸聲性能及其應(yīng)用[J]. 中國機械工程, 2007, 18(7): 778-780.

ZUO Yan-yan, ZHOU Jin-hua, LIU Hai-bo, et al. Absorption performance and application of perforated panel absorber[J]. China Mechanical Engineering,2007,18(7): 778-780.

[4]Selamet A, Xu M B, Lee I J. Analytical approach for sound attenuation in perforated dissipative silencers with inlet/outlet extension[J].Journal of Acoustical Society of America, 2005, 117(4): 2078-2089．

[5]高林. 多線譜可調(diào)頻水消聲器仿真設(shè)計研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2011.

[6]畢嶸. 復(fù)合式消聲器聲學(xué)特性的分析方法和實驗研究[D]. 合肥:合肥工業(yè)大學(xué), 2012.

[7]孟新,陸森林,劉紅光,等. 共振消聲器共振頻率的數(shù)值分析[J]. 重慶交通大學(xué)學(xué)報 (自然科學(xué)版),2012,31(5):1081-1085.

MENG Xin, LU Sen-lin, LIU Hong-guang, et al. Numerical analysis of resonant frequency of resonant muffler[J]. Journal of Chongqing Jiaotong University (Natural Science), 2012, 31(5): 1081-1085.

[8]羅虹,鄧海濤,董紅亮,等. 消聲器共振腔及穿孔隔板消聲特性數(shù)值分析[J]. 機械設(shè)計與制造, 2008 (3): 124-126.

LUO Hong, DENG Hai-tao, DONG Hong-liang, et al.Numerical analysis of mufflering characteristic for resonator and perforated clapboard in muffler[J]. Machinery Design and Manufacture, 2008 (3): 124-126.

[9]徐磊,劉正士,畢嶸. 結(jié)構(gòu)參數(shù)對直通穿孔管消聲器消聲性能影響的數(shù)值分析[J]. 合肥工業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2009, 32(11): 1637-1641.

XU Lei, LIU Zheng-shi, BI Rong. Numerical analysis of the effect of structural parameters on the acoustic attenuation performance of straight-through perforated tube silencers[J]. Journal of Hefei University of Technology(Natural Science), 2009, 32(11): 1637-1641.

[10]馬大猷. 現(xiàn)代聲學(xué)理論基礎(chǔ)[M]. 北京: 科學(xué)出版社, 2004: 210-220.

[11]Ingard K. On the theory and design of acoustic resonators[J]. Journal of Acoust. Soc. Am, 1953, 25: 1037-1061.

[12]徐貝貝,季振林. 穿孔管阻性消聲器聲學(xué)特性的有限元分析[J]. 振動與沖擊, 2010, 29(3): 58-62.

XU Bei-bei, JI Zhen-lin.Analysis on acoustical performance of perforated tube muffler by FEM[J]. Journal of Vibration and Shock, 2010, 29(3): 58-62.

[13]賀巖松,李景,盧會超. 組合式穿孔管消聲器聲學(xué)仿真[J]. 噪聲與振動控制, 2012, 32(2): 151-154.

HE Yan-song, LI Jing, LU Hui-chao.Acoustic simulation of mufflers with multiple perforated tubes[J]. Noise and Vibration Control, 2012, 32(2): 151-154.

[14]Munjal M L. Acoustics of ducts and mufflers[M]. John Wiley & Sons, 2014.

[15]劉克,焦風(fēng)雷,喬五之. 微穿孔板吸聲體非線性聲學(xué)特性初探[J]. 聲學(xué)技術(shù), 2001, 20(3): 135-136.

LIU Ke, JIAO Feng-lei, QIAO Wu-zhi. Rudimentary researches of nonlinear properties of microperforated-panel [J]. Acoustic Technology, 2001, 20(3): 135-136.

[16]何為,薛衛(wèi)東,唐斌. 優(yōu)化試驗設(shè)計方法及數(shù)據(jù)分析[M]. 北京: 化學(xué)工業(yè)出版社, 2011: 191-200.

[17]謝宇. 回歸分析[M].北京: 社會科學(xué)文獻出版社, 2013.