陳長征,白曉天,于慎波, 張 磊

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

?

2015200

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

陳長征1,2,白曉天1,于慎波1, 張 磊1

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

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

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

前言

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

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

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

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

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

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

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

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

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

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

(1)

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

kx2+ky2+kz2=k2

(2)

式中k為空氣中聲波波數(shù)。入射聲波首先通過穿孔管，為計算穿孔管內(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)

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

(7)

(8)

(9)

K(ω)= 101320×

(10)

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

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

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

(11)

(12)

(13)

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

(14)

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

(15)

2 阻抗復(fù)合式消聲器消聲性能實驗設(shè)計

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

3 實驗與計算結(jié)果分析

3.1 理論計算值與實測值對比

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

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

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

3.2 穿孔率對消聲器消聲量的影響

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

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

3.3 孔徑對消聲器消聲量的影響

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

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

3.4 吸聲材料流阻率對消聲器消聲量的影響

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

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

4 結(jié)論

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

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

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

[2] 劉晨,季振林,徐航手. 穿孔管消聲器聲學(xué)性能三維時域計算及分析[J]. 機械工程學(xué)報, 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] 宮建國,馬宇山,崔巍升,等. 汽車消聲器聲學(xué)特性的聲傳遞矩陣分析[J]. 振動工程學(xué)報, 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]. 汽車工程，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] 王偉,李國祥,倪計民,等. 不同結(jié)構(gòu)消聲器內(nèi)部流場的比較[J]. 內(nèi)燃機學(xué)報,2007(3):277-280.

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

[11] 畢嶸,劉正士,劉煥進，等. 阻抗復(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].長春：吉林大學(xué),2009.

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

[16] 方智,季振林. 直通穿孔管消聲器聲學(xué)特性預(yù)測的數(shù)值模態(tài)匹配法[J]. 聲學(xué)學(xué)報,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

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

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