袁守利，辛超，劉志恩，，肖生浩

（1.武漢理工大學(xué)汽車工程學(xué)院，武漢430070；2.現(xiàn)代汽車零部件技術(shù)湖北省重點實驗室，武漢430070）

汽車排氣消聲器三維聲場分析

袁守利1,2，辛超1,2，劉志恩，1,2，肖生浩1,2

（1.武漢理工大學(xué)汽車工程學(xué)院，武漢430070；2.現(xiàn)代汽車零部件技術(shù)湖北省重點實驗室，武漢430070）

采用三維有限元法對某車型排氣消聲器進行優(yōu)化設(shè)計，根據(jù)傳遞導(dǎo)納理論對消聲器穿孔管和穿孔板進行處理，建立數(shù)值模型并進行三維聲場仿真分析，獲得主副消聲器總成的傳遞損失；運用雙負載四傳聲器法測試消聲器傳遞損失，測試結(jié)果表明三維有限元法預(yù)測消聲器聲學(xué)性能有較高的精度，根據(jù)仿真結(jié)果和消聲器設(shè)計原理，對主消聲器進行優(yōu)化，可提高排氣系統(tǒng)聲學(xué)性能，滿足汽車噪聲排放法規(guī)的要求。

聲學(xué)；排氣消聲器；三維有限元；傳遞損失

為有效降低汽車噪聲給城市環(huán)境帶來的影響，機動車噪聲控制法規(guī)日益嚴格，汽車降噪技術(shù)與測量分析技術(shù)進入深層次的研究。目前運用GTPower軟件對汽車排氣噪聲的研究已是一種常用的方法，GT-Power軟件基于一維理論，在計算復(fù)雜消聲器或分析中高頻率噪聲時并不準確，會有高次模式波的出現(xiàn)，且一維模型僅局限于頻域范圍內(nèi)的橫波。國內(nèi)外已有人采用三維有限元法預(yù)測了消聲器的聲學(xué)性能[1―3]，經(jīng)過試驗驗證三維有限元法可以更為準確的計算消聲器內(nèi)部聲場，并且可以通過聲場分布計算得到消聲器腔體內(nèi)部各頻率段的聲壓分布云圖，更直觀表現(xiàn)出消聲器結(jié)構(gòu)對聲學(xué)性能的影響，為汽車排氣消聲器的設(shè)計及優(yōu)化提供理論指導(dǎo)。因此，本文采用三維聲學(xué)有限元法對某1.8 L排量車型排氣消聲器進行聲學(xué)性能評估，據(jù)此優(yōu)化消聲器結(jié)構(gòu)，并采用試驗證明了分析結(jié)果的有效性。

1 穿孔聲阻抗和吸聲材料的處理

傳遞損失是評價消聲器消聲性能的重要參量。采用有限元法計算消聲器傳遞損失時，有兩個難點需要處理：一是消聲器穿孔管和穿孔板上小孔的簡化；二是阻性消聲結(jié)構(gòu)中，吸聲材料的處理。

在進行消聲器設(shè)計時，為了減少流體的阻力和增大隔聲效果，會在消聲器內(nèi)部的插入管管壁或隔板上設(shè)計許多穿孔，聲波在通過這些小孔時，會有能量的衰減，由于這些小孔的存在，使得劃分有限元的網(wǎng)格變得非常困難；在建立有限于模型時必須對其進行簡化。為降低中高頻的噪聲，消聲器中需要填充吸聲材料，因此，如何在有限元模型中正確定義其屬性是一難點。

本文在建立消聲器聲學(xué)有限元網(wǎng)格時，忽略小孔，在穿孔板兩邊的網(wǎng)格之間通過定義一種傳遞導(dǎo)納關(guān)系，建立起穿孔板兩側(cè)的振動速度和聲壓的線性關(guān)系，間接模擬這些小孔。傳遞導(dǎo)納關(guān)系可以用公式（1）表達[4]

式中vn1和vn2為穿孔板兩側(cè)的法向振動速度；p1和p2為穿孔板兩側(cè)的聲壓；α1，α2，α4和α5為傳遞導(dǎo)納系數(shù)，α3和α6為由聲源引起來的系數(shù)，在消聲器計算中都是零，α1，α2，α4和α5通常需要試驗測定，聲音通過穿孔板的阻抗為Zp，采用季振林[5]穿孔聲阻抗計算公式，在無吸聲材料填充時，表示如下

式中t為穿孔板壁厚，dh為穿孔直徑，Φ為穿孔率，i為虛數(shù)，k0為波數(shù)，α為穿孔板的聲學(xué)厚度修正系數(shù)。

這樣可以得到

如果消聲器中添加吸聲材料，在進行三維有限元數(shù)值模擬分析時，吸聲材料用等效流體來替代，設(shè)吸聲材料的復(fù)阻抗和復(fù)波數(shù)分別為Zb、Kb，由Utsuno和Tanaka[6]計算公式得到吸聲材料的特性阻抗和波數(shù)的表達式為（空氣密度ρ0=1.156 kg/m3，聲速

式中σ為材料的流阻率，Z0為空氣的特性阻抗，K0為空氣中的波數(shù)。由于本文副消聲器中的吸聲材料為長纖維玻璃絲棉，其密度為100 kg/m3，參考吸聲材料流阻率測量方法[7]，測得該材料的流阻率為4 896 Rayls/m。

吸聲材料的復(fù)聲速和復(fù)密度可由下式得到

根據(jù)Kirby和Cumm ings[8]的研究，有吸聲材料填充時穿孔的特性聲阻抗表達式被修改為

2 聲學(xué)有限元模型的建立

對某汽車排氣系統(tǒng)中的主副消聲器總成分別進行分析，副消聲器為三腔結(jié)構(gòu)，兩隔板距副消聲器進氣口處截面分別為255 mm，445 mm，前隔板上打有200個直徑為4 mm的孔，后隔板上打有300個直徑為3 mm的孔，三個腔體中都有穿孔，其中第一腔中打有500個直徑為3 mm的孔，第二腔中打有324個直徑為3 mm的孔，第三腔中打有142個直徑為3.5 mm的孔，且第三腔中加入吸聲棉；主消聲器為三腔抗性結(jié)構(gòu)，兩隔板距消聲器進氣口端面分別為150 mm，450 mm，第二腔中有穿孔，其中進氣管上打有210個直徑為3.5 mm的孔，出氣管上打有60個直徑為3.5 mm的孔。副消聲器總成和主消聲器總成三維實體模型如圖1和2所示。

圖1 副消聲器總成

圖2 主消聲器總成

導(dǎo)入有限元軟件中建立其離散模型。有限元法計算聲場時常使用四面體和六面體單元，由于后者的計算精度要高，也為了節(jié)省單元數(shù)量和計算時間，選用六面體網(wǎng)格，單元大小為5 mm，對于內(nèi)部復(fù)雜的主消聲器，特別是帶有穿孔板的腔室，采用分塊建立，保證管壁內(nèi)外表面上的節(jié)點對應(yīng)，以便傳遞導(dǎo)納可以正確使用，隔板上的很多小孔也采用傳遞導(dǎo)納的關(guān)系來處理，網(wǎng)格劃分后的副消聲器和主消聲器總成如圖3和4所示：

圖3 副消聲器總成網(wǎng)格模型

圖4 主消聲器總成網(wǎng)格模型

計算消聲器傳遞損失時，在消聲器的入口端定義-1m/s速度振動邊界條件來給予激勵；消聲器的出口直接接著大氣，聲音通過出口直接傳播到大氣中，由于聲學(xué)有限元網(wǎng)格是實體網(wǎng)格，其壁面如果不做任何處理，聲音傳播到出口的壁面上后就會全部反射回來，為了模擬聲音在出口處沒有反射的效果，在出口端定義一個全吸聲的邊界條件，在消聲器的出口處定義聲阻抗Z0=ρ0c0=399.7 kg/(m2s)。

3 計算結(jié)果分析

通過設(shè)定聲學(xué)網(wǎng)格、流體材料、流體屬性、前處理操作、穿孔板、傳遞導(dǎo)納屬性、入口端和出口端、邊界條件、吸聲屬性等參數(shù)，設(shè)置計算頻率為10～3 000 Hz，步長為20 Hz，后進行三維聲學(xué)數(shù)值模擬。計算獲得主副消聲器總成的聲壓分布如圖5—圖10所示，分別為低、中、高三個頻率下消聲器的聲壓云圖分析。

圖5 150 Hz副消聲器總成聲壓分布圖

圖6 1 050 Hz副消聲器總成聲壓分布圖

圖7 2 450 Hz副消聲器總成聲壓分布圖

由圖5—圖7可以看出，在150 Hz時，第一腔和第二腔內(nèi)聲壓值較小為80～170 N/m2，而第三腔及副消聲器出氣管中，副消聲器聲壓值要稍大些為170～210 N/m2；在1 050 Hz時，副消聲器第二腔中的聲壓值要大，最大達到246 N/m2，第三腔中的聲壓值很小，最小為7.58 N/m2；在2 450 Hz時，第一腔和第二腔中的聲壓值更大，且不均勻，局部達到了952N/m2，第三腔中的聲壓值也很小，最小為2.57 N/m2。這是因為第一腔和第二腔為共振腔，主要是為了消除中低頻的噪聲，由于該車型的排氣副消聲器采用了三腔結(jié)構(gòu)，所以在150 Hz時，排氣噪聲的控制還較為理想，而第三腔中加有吸聲材料，主要是為了消除中高頻的噪聲，氣流經(jīng)過第三腔后，中高頻的噪聲得到了有效的控制。

圖8 150 Hz主消聲器總成聲壓分布圖

圖9 1 050 Hz主消聲器總成聲壓分布圖

圖10 2 450 Hz主消聲器總成聲壓分布圖

由圖8—圖10可以看出，在150 Hz時，主消聲器三腔內(nèi)聲壓值都較小，為0.2 N/m2～80.7 N/m2；在1 050 Hz，主消聲器第一腔和第二腔聲壓值較大，大部分區(qū)域是在40 N/m2～306 N/m2，而在第三腔中聲壓值要小些為20 N/m2～131 N/m2；在2 450 Hz時，主消聲器三個腔內(nèi)聲壓值都較大，且局部區(qū)域出現(xiàn)較大值，最大達到1340 N/m2，這是由于主消聲器第一腔為共振和擴張結(jié)構(gòu)，第二腔為共振結(jié)構(gòu)，第三腔為擴張結(jié)構(gòu)，都是為了消除中低頻的噪聲，所以中低頻的噪聲得到了很好的控制，而高頻噪聲沒有得到有效的控制。

為進一步了解副消聲器和主消聲器在各個頻段的消聲效果，在聲學(xué)有限元軟件中計算得到副消聲器總成和主消聲器總成的傳遞損失曲線。計算結(jié)果如圖11所示。

圖11 副消聲器和主消聲器總成傳遞損失

由圖11可以看出，副消聲器的消聲效果較理想，考慮到副消聲器體積有限，結(jié)構(gòu)已相對復(fù)雜，并且副消聲器的結(jié)構(gòu)變化對消聲器最后的消聲性能影響較小，因此對副消聲器的結(jié)構(gòu)不作優(yōu)化；主消聲器在10～3 000 Hz所有頻率下，消聲量都不高，最大只有36.38 dB，且在150 Hz以下，消聲量較小，為10～18.8 dB。通過改變主消聲器的結(jié)構(gòu)來改善排氣系統(tǒng)的消聲性能。

4 確定主消聲器優(yōu)化方案

根據(jù)聲學(xué)理論知識和消聲器設(shè)計原理，對于原主消聲器，進出氣管上在第二腔中均布置有穿孔高速氣流會通過出氣管中的小孔直接排出而沒有得到處理，且原主消聲器在中高頻降噪效果不好，所以需要在消聲器中增加吸聲材料來降低中高頻的噪聲，通常加在末腔中，以便平順消聲器所產(chǎn)生的酷利噪聲。優(yōu)化后的主消聲器采用三腔結(jié)構(gòu)，其中前隔板上打有220個直徑為3.5 mm的孔，在第一腔中，進氣管上打有90個直徑為3.5 mm的孔，第三腔中出氣管上打有190個直徑為3.5 mm的孔，且加入吸聲棉，結(jié)構(gòu)如圖12所示，以同樣的方法建立有限元模型，并分析計算得到其聲壓分布如圖13—圖15所示。

圖12 主消聲器優(yōu)化方案

圖13 150 Hz主消聲器優(yōu)化方案聲壓分布圖

圖14 1 050 Hz主消聲器優(yōu)化方案聲壓分布圖

圖15 2 450 Hz主消聲器優(yōu)化方案聲壓分布圖

由圖13—圖15可以看出，優(yōu)化后的主消聲器總成在150 Hz時，內(nèi)部聲壓值有明顯減小，為1.05 N/m2～50.4 N/m2，在1 050 Hz和2 450 Hz時，優(yōu)化后的主消聲器在第一腔和第二腔中，中高頻的噪聲并沒有得到提高，但是由于第三腔中添加了吸聲材料，所以在第三腔中的聲壓值較小，最大為132 N/m2，由此可見，經(jīng)過優(yōu)化后的主消聲器降噪量在各頻段都有所提高，為進一步體現(xiàn)優(yōu)化后的主消聲器的消聲效果，與原主消聲器對比，傳遞損失計算結(jié)果如圖16所示：

圖16 原主消聲器與主消聲器優(yōu)化方案傳遞損失對比

由圖16可以看出，優(yōu)化后的主消聲器的消聲量在10 Hz～3 000 Hz幾乎所有頻率段都有提高，特別是710 Hz以下和1 410 Hz～3 000 Hz有較大提高，710 Hz以下最大提高了26.9 dB，1 410 Hz～3 000 Hz最大提高了25.5 dB，主消聲器經(jīng)過優(yōu)化后，低頻和高頻噪聲都得到了有效提高，主消聲器的消聲性能相應(yīng)提高。

5 試驗驗證

由于兩負載法用于消聲器傳遞損失的測量操作方便，而且能測量到較精確的傳遞損失值，所以本文將采用兩負載法進行傳遞損失的測量[9]。傳聲損失測量原理如圖17所示。

圖17 測量系統(tǒng)原理圖

由圖17可以看出，布置在上游管一端的揚聲器發(fā)出穩(wěn)態(tài)寬帶隨機白噪聲給整個測量系統(tǒng)，被測試樣件布置在管中間部位，在上、下游管相對于被測試樣件的合適位置分別布置兩個傳聲器安裝位置，用于安裝壓力場型傳聲器。通過進行兩種不同典型邊界條件（下游管末端開口和使用海綿堵塞）進行測量，從而獲得消聲器傳聲損失。

5.1 消聲器傳遞損失試驗設(shè)備及儀器

本試驗主要采用B&K公司4206 T阻抗管實驗測量系統(tǒng)，在其基礎(chǔ)上來測量無流情況下消聲器傳遞損失。此套系統(tǒng)為固定件，所以聲源位置固定不變，采用兩負載法來進行消聲器傳遞損失的測量。試驗儀器設(shè)備參數(shù)見表1。

5.2 消聲器傳遞損失的測量與驗證

對原排氣主消聲器和副消聲器及優(yōu)化后的主消聲器均按此試驗方法進行試驗：由頻率發(fā)生器作為噪聲源，采用基于傳遞函數(shù)的4傳聲器位置、雙末端邊界條件的傳聲損失測量方法[10]。實驗溫度為20°C，副消聲器傳遞損失測試如圖18所示。

表1 儀器設(shè)備參數(shù)

圖18 副消聲器總成傳遞損失測試

圖18為副消聲器和主消聲器傳遞損失測量裝置圖。兩對傳聲器分別被安裝在被測聲學(xué)單元的上游和下游端。傳聲器被平嵌在管道上，所以傳聲器的存在不會破壞管道內(nèi)的原有聲場。聲源信號被四個傳聲器拾取，經(jīng)過功率放大器信號放大后輸入多通道分析儀進行分析運算。實驗要求在測量過程中信噪比一定要盡可能地高，至少大于10 dB。為了減少因信號的不確定性而產(chǎn)生的實驗錯誤，實驗中對信號采用了500次線性平均[9]。

測試時，消聲器前后連接管按實際排氣系統(tǒng)管路進行測試：副消聲器和主消聲器傳聲損失測試結(jié)果見圖19和圖20。

圖19 副消聲器總成傳遞損失模擬與試驗對比

由圖19和圖20可以看出，副消聲器和主消聲器總成的預(yù)測結(jié)果和實驗測量結(jié)果在整個頻域內(nèi)吻合很好，但在高頻出現(xiàn)稍微偏差，這是因為穿孔管和隔板采用的是聲阻抗計算公式，以及吸聲材料采用等效流體，這些經(jīng)驗公式在高頻還不夠準確，為了使三維有限元法可以更加準確的預(yù)測消聲器的聲學(xué)性能，還需要對穿孔聲阻抗計算公式講行進一步研究，尋找更加合適的修正參數(shù)。

圖20 主消聲器總成模擬與試驗對比

6 結(jié)語

通過對消聲器隔板或管壁上的小孔采用傳遞導(dǎo)納理論來模擬，消聲器中的吸聲材料采用等效流體的方法處理，可以快速準確實現(xiàn)三維聲學(xué)有限元方法預(yù)測消聲器的傳遞損失，由分析結(jié)果可以看出，副消聲器和主消聲器總成的預(yù)測結(jié)果和實驗測量結(jié)果在整個頻域內(nèi)吻合很好，對分析和改進消聲器的聲學(xué)性能提供了一定的依據(jù)。

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

[2]季振林.直通穿孔管阻性消聲器消聲性能計算及分析[J]振動工程學(xué)報，2005，18（4）：453-457.

[3]劉飛，陸森林，劉紅光，等.擴張室式消聲器有聲學(xué)特性的有限元的分析[J].噪聲振動與控制，2010，6（3）：144-147.

[4]李增剛，詹福良.聲學(xué)仿真計算高級應(yīng)用實例[M].北京：國防工業(yè)出版，2010，62-72.

[5]康鐘緒，季振林.穿孔板的聲學(xué)厚度修正[J].聲學(xué)學(xué)報，2008，33(4)：327-333.

[6]Utsuno H,Tanaka T.Transfer function method for measuring characteristic impedance and propagation constant of porous materials[J].Journal of the Acoustical Society of America,1989,86(2):637-643.

[7]ASTM:C522-03.Standard test method for air flow resistance of acoustical materials[S].American Society for Testing and Materials,2003.

[8]Kirby R,Cumm ings A.The impedance of perforated plates subjected to grazing gas flow and backed by porous media[J].Journal of Sound and Vibration，1998,217(4): 619-636.

[9]Munjal M L.Acoustic of ducts and mufflers[M].New York:AWiley-Interscience Publication,1987.

[10]康鐘緒.消聲器及穿孔元件聲學(xué)特性研究[D].哈爾濱：哈爾濱工程大學(xué)，2009.

Analysis of Three-dimensional Sound Field of s Exhaust Mufflers Automobile’

YUAN Shou-li1,2,XIN Chao1,2,LIU Zhien1,2,XIAO Sheng-hao1,2

(1.School ofAutomotive Engineering,Wuhan University of Technology,Wuhan 430070,China; 2.Hubei Key Laboratory ofAdvanced Technology of Automotive Parts,Wuhan 430070,China)

Three-dimensional finite element method is used to optimize the design of a vehicle’s exhaust mufflers.Based on the transfer admittance theory,behaviors of the perforated pipe and the perforated plate of a muffler are analyzed.The numerical model is established and the three-dimensional acoustic fields of the primary and secondary mufflers are simulated and analyzed,and their transmission losses are obtained.Then,the transmission losses of the primary and secondary mufflers are measured with the double-load and four-microphone method.The experiment results show that the three-dimensional finite element method has high accuracy to predict the muffler’s acoustic performance.According to the simulation results and the principle of muffler’s design,the primary muffler is optim ized,and the acoustic performance of its exhaust system is improved.The demands of laws and regulations for the vehicle’s noise emission limitation are satisfied.

acoustics;exhaust muffler;three-dimensional FEM;transm ission loss

TB5

A

10.3969/j.issn.1006-1335.2014.01.026

1006-1355(2014)01-0113-05

2013-04-08

國家科技支撐計劃—生物燃氣凈化提純技術(shù)裝備及城市公交系統(tǒng)應(yīng)用研究（基金編號：2012BAC18B02）

袁守利（1966-），男，武漢人，副教授，目前從事車輛新能源及排放控制研究。

E-mail:YSL0202@163.com