結(jié)構(gòu)參數(shù)對(duì)阻抗復(fù)合消聲器消聲量的影響分析*

陳長(zhǎng)征,白曉天,于慎波, 張 磊

(1.沈陽(yáng)工業(yè)大學(xué)機(jī)械工程學(xué)院，沈陽(yáng) 110870； 2.遼寧省振動(dòng)噪聲控制工程技術(shù)研究中心，沈陽(yáng) 110870)

?

2015200

結(jié)構(gòu)參數(shù)對(duì)阻抗復(fù)合消聲器消聲量的影響分析*

陳長(zhǎng)征1,2,白曉天1,于慎波1, 張 磊1

(1.沈陽(yáng)工業(yè)大學(xué)機(jī)械工程學(xué)院，沈陽(yáng) 110870； 2.遼寧省振動(dòng)噪聲控制工程技術(shù)研究中心，沈陽(yáng) 110870)

本文中應(yīng)用橢圓柱坐標(biāo)系下的三維解析算法對(duì)橢圓形阻抗復(fù)合式消聲器的聲學(xué)特性進(jìn)行了研究，分析了不同參數(shù)對(duì)消聲器消聲量的影響規(guī)律。結(jié)果表明，穿孔管穿孔率、孔徑和吸聲材料流阻率均對(duì)消聲器的消聲特性有較大影響，不同參數(shù)下消聲器在高、低頻段的消聲量有較大差異。實(shí)際應(yīng)用中建議根據(jù)不同工況下排氣噪聲頻率特性來調(diào)整消聲器結(jié)構(gòu)參數(shù)，以獲得滿意的效果。

阻抗復(fù)合式消聲器；橢圓形斷面；消聲性能；結(jié)構(gòu)參數(shù)

前言

隨著汽車工業(yè)的發(fā)展，汽車的振動(dòng)噪聲粗糙度(noise vibration harshness, NVH)性能逐漸成為了衡量汽車品質(zhì)的重要指標(biāo)之一。排氣噪聲是汽車噪聲總成中重要的一部分，如何針對(duì)汽車排氣噪聲設(shè)計(jì)全頻段消聲器,成為當(dāng)今汽車行業(yè)中的重要課題[1-3]。鑒于汽車排氣噪聲的寬頻帶特性，采用結(jié)構(gòu)單一的單腔消聲器不能達(dá)到好的效果，因此根據(jù)不同工況下汽車排氣噪聲頻帶特性研制多級(jí)阻抗復(fù)合式消聲器已經(jīng)成為研究熱點(diǎn)[4-6]。

國(guó)內(nèi)外學(xué)者針對(duì)阻抗復(fù)合式消聲器進(jìn)行了大量研究，對(duì)于不同結(jié)構(gòu)的消聲器消聲性能的計(jì)算方法主要有一維平面波法、傳遞矩陣法、有限元與邊界元法和解析法。一維平面波法和傳遞矩陣法基于一維平面波理論，可對(duì)簡(jiǎn)單結(jié)構(gòu)消聲器進(jìn)行定性分析，但由于其理論基礎(chǔ)忽略了高階次聲波的影響，在對(duì)中高頻率和復(fù)雜結(jié)構(gòu)的消聲器計(jì)算中結(jié)果誤差較大。有限元和邊界元法計(jì)算結(jié)果準(zhǔn)確，可用于結(jié)構(gòu)復(fù)雜的消聲器聲學(xué)性能分析，但其計(jì)算量大，計(jì)算時(shí)間長(zhǎng)，模型改變后需要重新計(jì)算，不適用于消聲器優(yōu)化設(shè)計(jì)。解析法在滿足計(jì)算精度的同時(shí)，計(jì)算量較小，計(jì)算時(shí)間較短，可以較方便地進(jìn)行變參計(jì)算，適用于消聲器的性能分析和優(yōu)化設(shè)計(jì)[7-9]。在相關(guān)問題的研究中，文獻(xiàn)[10]中基于平面波傳播假設(shè)并使用解耦方法，在三維子結(jié)構(gòu)邊界元法的基礎(chǔ)上導(dǎo)出阻抗矩陣，通過阻抗矩陣預(yù)測(cè)和分析穿孔管消聲器的消聲性能。平面波理論計(jì)算方便，然而用于計(jì)算阻抗復(fù)合式消聲器時(shí)不能滿足計(jì)算精度要求；文獻(xiàn)[11]中采用二維分析法推導(dǎo)了阻抗復(fù)合式消聲器的聲學(xué)模型，并通過計(jì)算得出阻抗復(fù)合式消聲器吸聲材料的填充位置對(duì)消聲器在不同頻段上的消聲性能會(huì)產(chǎn)生一定的影響，提出可考慮在消聲器內(nèi)填充不同吸聲材料以拓寬消聲頻帶。然而二維分析法具有一定局限性，無(wú)法用于計(jì)算橢圓形聲腔的消聲器聲學(xué)性能；文獻(xiàn)[12]中基于空間場(chǎng)格林函數(shù)推導(dǎo)得出可用于計(jì)算橢圓消聲器消聲性能的傳遞矩陣，針對(duì)消聲器進(jìn)出口位置對(duì)傳遞損失的影響做了細(xì)致分析。然而其消聲器模型為抗性，未能考慮阻抗復(fù)合情況；文獻(xiàn)[13]中提出了使用穿孔管的多級(jí)消聲器的傳遞矩陣解析方法，并將解析解與有限元法結(jié)果進(jìn)行對(duì)比，驗(yàn)證了解析方法的有效性，之后在此基礎(chǔ)上分別對(duì)不同結(jié)構(gòu)的圓柱形消聲器進(jìn)行了傳遞損失計(jì)算與分析，得出了優(yōu)化參數(shù)，但其研究模型具有一定局限性，未考慮添加阻性材料后聲波波數(shù)的變化。

汽車排氣系統(tǒng)噪聲在不同工況下表現(xiàn)出的頻率特性相差很大，因此研究結(jié)構(gòu)參數(shù)對(duì)消聲器消聲性能影響，并根據(jù)不同工況排氣噪聲聲源特性調(diào)整消聲器內(nèi)部結(jié)構(gòu)參數(shù)，以提升全工況排氣系統(tǒng)NVH性能十分必要。本文中采用基于橢圓柱坐標(biāo)系的三維解析法計(jì)算消聲器的聲學(xué)特性，著重研究了穿孔管孔徑、穿孔率等參數(shù)對(duì)消聲器消聲特性的影響趨勢(shì)，以提升全工況下消聲器的全頻段消聲性能，拓寬消聲頻帶。

1 阻抗復(fù)合式消聲器的建模與理論分析

1.1 橢圓形阻抗復(fù)合式消聲器建模

一種橢圓形阻抗復(fù)合式消聲器三維示意圖如圖1所示。消聲器由3個(gè)消聲腔組成，其外殼及內(nèi)部消聲腔均為橢圓形，消聲器兩端為抗性腔，中間為阻性腔，阻性腔中吸聲材料填充在穿孔管外殼體內(nèi)的空間內(nèi)。各腔之間有隔板，以保證結(jié)構(gòu)上的相對(duì)獨(dú)立。這種結(jié)構(gòu)的優(yōu)點(diǎn)在于保證消聲量的同時(shí)流阻較小，且故障率低和便于維修與更換。

1.2 腔中聲波傳遞規(guī)律計(jì)算

消聲器的結(jié)構(gòu)簡(jiǎn)圖如圖2所示。橢圓形消聲器左端為入口，右端為出口，半長(zhǎng)軸與半短軸長(zhǎng)度分別為ra與rb，消聲器各腔長(zhǎng)度分別為L(zhǎng)1，L2和L3，3段穿孔管直徑均為d。設(shè)消聲器軸向方向?yàn)閤向，半長(zhǎng)軸與半短軸方向分別為y向與z向。

假設(shè)消聲器內(nèi)部空氣為無(wú)吸收衰減、無(wú)黏性、無(wú)熱傳導(dǎo)的理想流體，根據(jù)赫姆霍茲方程，可得消聲器入口處聲壓表達(dá)式[14]為

pI=Ae-j(kxx+kyy+kzz)

(1)

式中：A為實(shí)數(shù)；kx，ky和kz分別為沿x，y和z軸方向的聲波波數(shù)，且

kx2+ky2+kz2=k2

(2)

式中k為空氣中聲波波數(shù)。入射聲波首先通過穿孔管，為計(jì)算穿孔管內(nèi)聲壓，引入圓柱坐標(biāo)系ρφx，其拉普拉斯算子為

(3)

將式(3)代入赫姆霍茲方程，可解得穿孔管內(nèi)部聲壓解為

p(ρ,x)= [AmJm(kρρ)cos(mφ)+

BmJm(kρρ)sin(mφ)]ej(ωt-kxx)

(4)

(5)

將式(5)帶入赫姆霍茲方程，解得抗性腔中聲壓解為

(6)

聲波通過抗性腔后，進(jìn)入阻性腔。阻性腔中聲波傳遞規(guī)律可分兩部分來討論，即穿孔管內(nèi)介質(zhì)為空氣的圓柱形腔與穿孔管外介質(zhì)為吸聲材料的橢圓柱形腔。其中，穿孔管內(nèi)部聲波傳遞規(guī)律與式(4)形式相同，而穿孔管外吸聲材料中的聲壓解為

(7)

(8)

(9)

K(ω)= 101320×

(10)

式中：ω為聲波頻率;ρ0為空氣密度;σ為吸聲材料流阻率。

聲波在通過穿孔管時(shí)存在一定聲壓降，根據(jù)聲壓連續(xù)性條件，有

p+=p--ρ0c0ξpV-

(11)

(12)

(13)

由式(4)及式(6)～式(13)可知，消聲器內(nèi)聲壓分布與穿孔率、孔徑和吸聲材料流阻率均有關(guān)。由于消聲器由3個(gè)消聲腔組成，消聲器的總消聲量為3個(gè)消聲腔的總和。每個(gè)消聲腔的消聲量為

(14)

式中：N為消聲腔的序號(hào)，本文中N=1,2,3，并設(shè)L0=0。則消聲器的總消聲量為

(15)

2 阻抗復(fù)合式消聲器消聲性能實(shí)驗(yàn)設(shè)計(jì)

實(shí)驗(yàn)原理簡(jiǎn)圖如圖3所示。消聲器結(jié)構(gòu)與圖1相同，ra=115mm，rb=80mm，各腔長(zhǎng)度L1=135mm，L2=190mm，L3=130mm。穿孔管直徑d=50mm，壁厚tw=0.7mm。通過采集消聲器入口和出口處的測(cè)量數(shù)據(jù)來計(jì)算消聲器的消聲量。

3 實(shí)驗(yàn)與計(jì)算結(jié)果分析

3.1 理論計(jì)算值與實(shí)測(cè)值對(duì)比

在采用相同參數(shù)(各腔穿孔管穿孔率均為Φ=4%，孔徑均為dh=4.2mm，吸聲材料流阻率為σ=10000Rayls/m)的條件下，阻抗復(fù)合式消聲器消聲量的三維解析法和有限元法的計(jì)算值與實(shí)際測(cè)量值的對(duì)比如圖4所示。

由圖4可見，采用三維解析法的計(jì)算值與實(shí)測(cè)值數(shù)據(jù)在全頻域內(nèi)吻合較好，總的趨勢(shì)一致。在低頻段有一定誤差，這與三維解析法的計(jì)算方法有關(guān)。三維解析法將所有隔板處設(shè)置為剛性反射壁面進(jìn)行計(jì)算，而實(shí)際消聲器中隔板厚度僅為1mm，在反射的同時(shí)還存在聲波透射現(xiàn)象，因此實(shí)際消聲量要小于理論計(jì)算值。此外，在中頻段出現(xiàn)了消聲頻率峰值處消聲量的差異，這可能是由于算法未考慮聲振耦合效應(yīng)對(duì)消聲器消聲量的影響，聲音傳播導(dǎo)致消聲器結(jié)構(gòu)產(chǎn)生振動(dòng)，而振動(dòng)模態(tài)與聲腔模態(tài)產(chǎn)生了耦合現(xiàn)象，從而影響了消聲器的消聲量。通過比較三維解析法與有限元法的計(jì)算結(jié)果可知，三維解析法與有限元法計(jì)算結(jié)果相近，與實(shí)驗(yàn)測(cè)量數(shù)據(jù)相比精度大致相同。

總體來說，采用三維解析法計(jì)算值與實(shí)驗(yàn)數(shù)據(jù)吻合較好，計(jì)算精度較高，可以用來進(jìn)行橢圓形阻抗復(fù)合式消聲器的消聲量模擬計(jì)算。下面用三維解析法研究不同參數(shù)對(duì)消聲器消聲量的影響規(guī)律，以得到進(jìn)一步改善消聲器全頻帶消聲性能的優(yōu)化方向。

3.2 穿孔率對(duì)消聲器消聲量的影響

圖5為采用不同穿孔率的橢圓形阻抗復(fù)合式消聲器消聲量的理論計(jì)算值，消聲器參數(shù)選取孔徑dh=4.2mm，吸聲材料流阻率σ=10000Rayls/m，穿孔管穿孔率Φ分別為4%，6%，8%。

由圖5可見，在低頻段，使用穿孔率較小的穿孔管的消聲器消聲量較大，在高頻段，使用穿孔率較大的穿孔管的消聲器消聲量較大。而且隨著穿孔率的增大，消聲器在全頻段的消聲量最大值與最小值的差異逐漸增大。這是由于在阻抗復(fù)合式消聲器中，阻性腔主要消聲原理為吸聲材料吸聲，其高頻消聲效果較好，低頻消聲效果較差，而抗性腔主要消聲原理為共振消聲，其低頻消聲效果較好，高頻消聲效果較差。穿孔率決定了聲波與吸聲材料接觸的概率，在其他參數(shù)不變的情況下，隨著穿孔率的增大，穿孔管內(nèi)傳播的聲波與吸聲材料的接觸變得頻繁，消聲腔也由較封閉向較開放逐漸轉(zhuǎn)變，阻性腔的消聲量比重增大，抗性腔的消聲量比重減小，因此消聲器的消聲量在低頻段較小，在高頻段較大。

3.3 孔徑對(duì)消聲器消聲量的影響

圖6為采用不同孔徑的橢圓形阻抗復(fù)合式消聲器消聲量的理論計(jì)算值，消聲器參數(shù)選取穿孔管穿孔率Φ=4%，吸聲材料流阻率σ=10000Rayls/m，孔徑分別為3，4.2和6mm。

由圖6可見，在低頻段，孔徑對(duì)消聲器消聲量的影響不是太有規(guī)律，在高頻段(400Hz以上)，隨著孔徑的加大，消聲器的消聲量增大，孔徑對(duì)消聲量影響趨勢(shì)與穿孔率對(duì)消聲量影響趨勢(shì)基本相似。因?yàn)榭讖經(jīng)Q定了穿孔管內(nèi)聲波與吸聲材料接觸的面積，在其他參數(shù)不變的情況下，隨著孔徑的增大，穿孔管內(nèi)聲波能夠更多地由孔傳入到吸聲材料中，使吸聲材料的作用發(fā)揮得更充分，從而使排氣噪聲中的高頻部分得到更好地衰減，同時(shí)消聲腔也由較封閉向較開放逐漸轉(zhuǎn)變，阻性腔的消聲量比重增大。

3.4 吸聲材料流阻率對(duì)消聲器消聲量的影響

圖7為采用不同流阻率的吸聲材料的橢圓形阻抗復(fù)合式消聲器的理論計(jì)算值，消聲器參數(shù)選取穿孔管穿孔率Φ=4%，孔徑dh=4.2mm，穿孔管外空腔內(nèi)吸聲材料流阻率分別為σ=5000Rayls/m，σ=10000Rayls/m，σ=15000Rayls/m。

由圖7可見，吸聲材料流阻率對(duì)消聲量的影響規(guī)律和孔徑類似。在高頻段(200Hz以上)，隨著選用吸聲材料流阻率的增大，消聲器的消聲量加大，且峰值稍向低頻方向移動(dòng)，總體而言，選用流阻率大的吸聲材料的消聲器的消聲量在全頻段消聲量較大。由于吸聲材料流阻率與穿孔管外消聲腔內(nèi)阻抗成正比，且吸聲材料緊貼穿孔管放置，對(duì)聲波透過穿孔管的壓力損失也有所影響。在其他條件不變的情況下，隨著吸聲材料流阻率的增大，消聲腔內(nèi)聲波傳遞阻抗增大，消聲腔逐漸由較開放向較封閉轉(zhuǎn)變，抗性腔消聲量比重增大，所以消聲量峰值會(huì)隨著吸聲材料流阻率增大而向低頻方向偏移。

4 結(jié)論

(1)本文中應(yīng)用橢圓柱坐標(biāo)系下的三維解析法對(duì)阻抗復(fù)合式消聲器的聲學(xué)性能進(jìn)行了理論推導(dǎo)與實(shí)驗(yàn)研究，將消聲器消聲量理論值與實(shí)測(cè)值進(jìn)行對(duì)比，證明了三維解析法的全頻段適用性，并將解析法計(jì)算值與有限元法計(jì)算值進(jìn)行了對(duì)比，證明解析法精度可以滿足要求。

(2)采用三維解析法分析了橢圓形阻抗復(fù)合式消聲器的消聲量隨穿孔管穿孔率、孔徑和消聲腔內(nèi)吸聲材料流阻率對(duì)消聲量的影響規(guī)律。結(jié)果表明，增大穿孔率、增大孔徑、減小吸聲材料流阻率都會(huì)增加高頻段的消聲量，但可能會(huì)相應(yīng)減少較低頻段的消聲量，在實(shí)際應(yīng)用中可根據(jù)不同工況下排氣噪聲聲源特性選取不同的內(nèi)部結(jié)構(gòu)參數(shù)。本研究為現(xiàn)階段消聲器內(nèi)部結(jié)構(gòu)參數(shù)的選取提供了依據(jù)，并為未來消聲器的設(shè)計(jì)提供了參考。

[1] 畢嶸，劉正士，陸益民,等.多入口多出口抗性消聲器的聲學(xué)性能研究[J]. 汽車工程，2014,36(2):243-253.

[2] 劉晨,季振林,徐航手. 穿孔管消聲器聲學(xué)性能三維時(shí)域計(jì)算及分析[J]. 機(jī)械工程學(xué)報(bào), 2012,48(10):7-13.

[3] Munjal M L. Analysis and Design of Pod Silencers[J]. Journal of Sound and Vibration, 2003,262:497-507.

[4] Carl Q Howard, Richard A Craig. Noise Reduction Using a Quarter Wave Tube with Different Orifice Geometries[J]. Applied Acoustics,2014,76:180-186.

[5] 宮建國(guó),馬宇山,崔巍升,等. 汽車消聲器聲學(xué)特性的聲傳遞矩陣分析[J]. 振動(dòng)工程學(xué)報(bào), 2010, 23(6):636-641.

[6] Torregrosa A J, Fajardo P, Gil A, et al. Development of Non-reflecting Boundary Condition for Application in 3D Computational Fluid Dynamics Codes[J]. Engineering Applications of Computational Fluid Mechanics, 2012, 6(3):447-460.

[7] 侯鎖軍,杜艷霞. 基于整車噪聲控制的消聲器分析與改進(jìn)[J]. 汽車工程，2013,35(10):944-954.

[8] Min-Chie Chiu. Shape Optimization of Multi-Chamber Mufflers with Plug-inlet Tube on a Venting Process by Genetic Algorithms[J]. Applied Acoustics, 2010, 71:494-505.

[9] 王偉,李國(guó)祥,倪計(jì)民,等. 不同結(jié)構(gòu)消聲器內(nèi)部流場(chǎng)的比較[J]. 內(nèi)燃機(jī)學(xué)報(bào),2007(3):277-280.

[10] 季振林.直通穿孔管消聲器聲學(xué)性能計(jì)算及分析[J].哈爾濱工程大學(xué)學(xué)報(bào), 2005, 26(3):302-306.

[11] 畢嶸,劉正士,劉煥進(jìn)，等. 阻抗復(fù)合式消聲器聲學(xué)特性的研究[J].汽車工程, 2012, 34(7):633-637.

[12] Mimani A,Munjal M L. 3-D Acoustic Analysis of Elliptical Chamber Mufflers Having an End-inlet and a Side-outlet: an Impedance Matrix Approach[J]. Wave Motion, 2012, 49:271-295.

[13] Vijayasree N K,Munjal M L. On an Integrated Transfer Matrix Method for Multiply Connected Mufflers[J]. Journal of Sound and Vibration, 2012, 331:1926-1938.

[14] 馬家義. 車用發(fā)動(dòng)機(jī)排氣噪聲有源主動(dòng)控制管道模擬研究[D].長(zhǎng)春：吉林大學(xué),2009.

[15] 劉剛田,吉曉民. 汽車駕駛室吸聲材料的吸聲系數(shù)計(jì)算方法[J]. 機(jī)械工程學(xué)報(bào),2014(12):104-109.

[16] 方智,季振林. 直通穿孔管消聲器聲學(xué)特性預(yù)測(cè)的數(shù)值模態(tài)匹配法[J]. 聲學(xué)學(xué)報(bào),2013(5):607-614.

An Analysis on the Influence of Structural Parameters on the Sound Attenuation of Impedance Compound Muffler

Chen Changzheng1,2, Bai Xiaotian1, Yu Shenbo1& Zhang Lei1

1.SchoolofMechanicalEngineering,ShenyangUniversityofTechnology,Shenyang110870; 2.EngineeringResearchCenterforNoiseandVibrationControlofLiaoningProvince,Shenyang110870

In this paper, the acoustic performance of impedance compound muffler with an elliptical cross section is studied by using three-dimensional analytic algorithm derived in an elliptical cylindrical coordinates to analyze the effects of different parameters on the sound attenuation of muffler. The results show that the perforation rate, the aperture and the flow resistivity of sound absorption materials all have relatively significant effects on the acoustic attenuation characteristics of muffler. The magnitudes of sound attenuation at high and low frequencies have rather great difference for mufflers with different parameters. So it is strongly suggested in practical application to adjust the structural parameters of muffler according to the frequency characteristics of exhaust noise in different conditions for obtaining satisfactory results.

impedance compound muffler; elliptical cross section; acoustic attenuation performance; structural parameters

*國(guó)家自然科學(xué)基金(51175350)資助。

原稿收到日期為2015年1月8日，修改稿收到日期為2015年4月8日。