鐘庭生，李亞光，圣小珍，肖新標(biāo)，杜幾平，劉謀凱

（西南交通大學(xué) 牽引動力國家重點實驗室，成都 610031）

運用模態(tài)疊加法研究擴張式消聲器的傳遞損失

鐘庭生，李亞光，圣小珍，肖新標(biāo)，杜幾平，劉謀凱

（西南交通大學(xué) 牽引動力國家重點實驗室，成都 610031）

采用聲學(xué)模態(tài)疊加法建立單腔擴張式消聲器傳遞損失計算模型，然后通過Matlab編程實現(xiàn)單腔擴張式消聲器傳遞損失的數(shù)值計算。在此基礎(chǔ)上，比較聲學(xué)模態(tài)疊加法、有限元法和基于平面波假定的經(jīng)典公式法在計算單腔擴張式消聲器傳遞損失上的差別，研究單腔擴張式消聲器膨脹段尺寸對傳遞損失的影響。結(jié)果表明，對于平面入射波，聲學(xué)模態(tài)疊加法可用于單腔擴張式消聲器各頻段傳遞損失的計算；增大膨脹段的半徑能有效提高低頻段的傳遞損失，但對高頻段的影響較?。浑S著膨脹段寬度的增大，傳遞損失的峰值向低頻移動，傳遞損失最大的頻段向高頻移動。

聲學(xué)；單腔擴張式消聲器；聲學(xué)模態(tài)疊加法；傳遞損失；Matlab

飛機和汽車發(fā)動機的進排氣系統(tǒng)、建筑物的通風(fēng)和空調(diào)系統(tǒng)由許多管路組成，管路在換氣、輸送冷熱氣體的時候，也會傳播噪聲。消聲器主要用于消除管道中的噪聲，一般可以分為抗性消聲器和阻性消聲器兩大類[1]。抗性消聲器主要以聲抗性原理使原在管道中傳播的噪聲遇到加裝的消聲器后，被反射回去，而不是將傳播中的噪聲能量消耗掉[1]。因為它具有明顯的頻率選擇特性，所以特別適用于在管道噪聲的頻譜中具有明顯的峰值特征的情況，一般適用于低中頻段[1]。擴張式消聲器是工程上常用的抗性消聲器，其原理是通過膨脹段來改變聲波傳播路徑上的阻抗，從而使入射波發(fā)生反射，反射得越多，則傳到下游段的透射波就越少，就越能達到降噪的目的。

消聲器的傳遞損失定義為消聲器出口無反射時，消聲器進口處的入射聲功率級與出口處的透射聲功率級之差[2]。消聲器的傳遞損失反映了其本身的聲波傳遞特性，不受聲源、管道系統(tǒng)的影響?；趯θ肷洳?、反射波和透射波均作平面波假定的單腔擴張式消聲器的傳遞損失經(jīng)典公式為

式中R1為管道半徑，R2為膨脹段半徑，L為膨脹段的寬度，k=ω/c為聲波的波數(shù)，ω為聲波的圓頻率，c為聲速。但是，如果入射波不是平面波，且其頻率高于消聲器的截止頻率，則基于平面波假定的公式就不能用，這時需要用聲學(xué)模態(tài)疊加法、聲學(xué)有限元法、聲學(xué)邊界元法等方法來計算它的傳遞損失。

國內(nèi)外學(xué)者對消聲器的傳遞損失進行了大量研究，Davis等將一維解析法引入到擴張式消聲器的計算當(dāng)中，并且給出了消聲器傳遞損失的計算公式[3]。Munjal等利用傳遞矩陣法計算了消聲器的傳遞損失[4]。Miles采用多維解析法計算了消聲器的傳遞損失[5]。Ji等應(yīng)用邊界元方法計算了考慮運流效應(yīng)的消聲器傳遞損失，對于尺寸較大的結(jié)構(gòu)，提出了基于阻抗矩陣的子結(jié)構(gòu)邊界元法[6]。Selamet等應(yīng)用模態(tài)匹配方法計算了擴張式消聲器的傳遞損失[7]。Wu等通過建立擴張式消聲器傳遞損失的理論模型，探討了不同進口彎角對消聲器聲學(xué)性能的影響[8]。宮建國等利用聲傳遞矩陣法，計算了消聲器的傳遞損失，分析了汽車消聲器的聲學(xué)特性，討論了不同參數(shù)對消聲器聲學(xué)性能的影響[9]。王鍵等利用有限元法分析了擴張比、內(nèi)插管長度以及擴張腔的個數(shù)等因素對抗性消聲器聲學(xué)性能的影響，表明改進消聲器的結(jié)構(gòu)參數(shù)可以有效提高消聲器的傳遞損失[10]。

本文采用聲學(xué)模態(tài)疊加法建立單腔擴張式消聲器傳遞損失的計算模型，然后利用有限元法和基于平面波假定的經(jīng)典公式法來驗證該模型，最后利用所建立的模型討論膨脹段的尺寸參數(shù)對消聲器傳遞損失的影響。

1 運用聲學(xué)模態(tài)疊加法計算傳遞損失

1.1 基本假設(shè)[1]

為使問題簡化，對媒質(zhì)和聲波過程作如下假設(shè)和簡化：

（1）媒質(zhì)為理想流體；

（2）沒有聲擾動時，媒質(zhì)在宏觀上是靜止的，同時媒質(zhì)是均勻的；

（3）聲傳播時，媒質(zhì)中的稠密和稀疏過程是絕熱的；

（4）媒質(zhì)中傳播的是小振幅聲波，各聲學(xué)量是一級微量。

1.2 基本參數(shù)

如圖1所示，假設(shè)圓形截面單腔擴張式消聲器上游段和下游段的半徑為R1，膨脹段的半徑為R2，寬度為L，各段的軸線共線且各段內(nèi)壁均為光滑剛性壁面。管內(nèi)聲場用柱坐標(biāo)(r,θ,z)來描述。

圖1 單腔擴張式消聲器示意圖

對于管內(nèi)頻率為ω的簡諧波，其聲壓和空氣質(zhì)點速度幅值均是θ的周期函數(shù)，周期是2π，可以利用傅里葉級數(shù)分解為所有環(huán)向諧波分量的疊加，即

其中m是環(huán)向諧波指數(shù)。

所以本文只考慮入射波是環(huán)向諧波指數(shù)為m、徑向模態(tài)階次為s的模態(tài)聲波[1–2]。

式中κ1,s為上游段的第s(s=0,1,2,…,)階徑向波數(shù)，可通過式(4)求得，即κ1,s為式(4)中κ的第(s+1)個非負實根[1–2]，為上游段的第s階軸向波數(shù)為聲波的波數(shù)，vi1z為上游段入射波的軸向速度，Jm為m階柱貝塞爾函數(shù)。

1.3 傳遞損失計算

一般來說，膨脹段對上游段是一個聲負載，下游段對膨脹段是一個聲負載，因而會引起部分聲波的反射和透射。設(shè)在上游段有一入射波pi1和一反射波pre1，膨脹段有一入射波pi2和一反射波pre2，而下游段無限延伸，僅有透射波pt3。入射波是環(huán)向諧波指數(shù)為m、徑向模態(tài)階次為s的模態(tài)聲波，所以反射到上游段、在膨脹段和下游段中的波關(guān)于環(huán)向角度θ的變化規(guī)律不變，但關(guān)于徑向r和軸向z的變化規(guī)律要改變，反射到上游段、在膨脹段和下游段中的波可以看成是相應(yīng)的聲學(xué)模態(tài)波的疊加。如圖1所示，假定坐標(biāo)原點取在上游段和下游段的接口處，可以寫出以下五種波的聲壓和質(zhì)點速度表達式。

(1)上游段各處聲壓p1和質(zhì)點軸向速度v1z

式中N1為上游段反射波所取的模態(tài)波的階數(shù)，An為系數(shù)，κ1,n為上游段的第n階徑向波數(shù)，為上游段的第n階軸向波數(shù)，ρ0為空氣密度，其它同上。

(2)膨脹段各處聲壓p2和質(zhì)點軸向速度v2z

式中N2為膨脹段入射波和反射波所取的模態(tài)波的階數(shù)，Bn、Cn為系數(shù)，κ2,n為膨脹段的第n階徑向波數(shù)，可通過式(8)求得，即κ2,n為式(8)中κ的第(n+1)個非負實根，為膨脹段第n階軸向波數(shù)，其它同上。

(3)下游段各處聲壓p3和質(zhì)點軸向速度v3z

式中Dn為系數(shù)，其它同上。

(4)邊界條件

在上游段和膨脹段接口處(z=0)由聲壓連續(xù)條件得

將式(5)和式(7)代入式(12)得

由軸向質(zhì)點速度連續(xù)條件得

將式(6)和式(9)代入式(14)得

在膨脹段和下游段接口處(z=L)，由聲壓連續(xù)條件得

將式(7)和式(10)代入式(16)得

由軸向質(zhì)點速度連續(xù)條件得

將式(9)和式(11)代入式(18)得

根據(jù)上下游段半徑R1和需要計算的最大頻率fmax從r=0到r=R1沿徑向等間距（間距為dr）取N1個點，再從r=R1+dr到r=R2以同樣的間距dr取(N2-N1)個點，再將這N2個點的徑向坐標(biāo)r(i)代入式(13)、式(15)、式(17)和式(19)得到 2(N1+N2)個含有2(N1+N2)個未知量的線性代數(shù)方程組，通過數(shù)值計算可求得這些未知量，再代入式(5)、式(6)、式(7)、式(9)、式(10)、式(11)即可得到不同位置的聲壓和質(zhì)點軸向速度。

(5)計算傳遞損失

入射波和透射波的平均聲功率分別為

單腔擴張式消聲器的傳遞損失為

將式(3)代入式(20)得入射波的平均聲功率，將式(10)、式(11)代入式(21)得透射波的平均聲功率，再代入式(22)得單腔擴張式消聲器的傳遞損失

2 數(shù)值計算與結(jié)果討論

本文計算的單腔擴張式消聲器各段內(nèi)壁均為光滑剛性壁面，上下游段半徑R1=0.05 m，空氣密度ρ0=1.21 kg/m3，聲速c=343 m/s，距離單腔擴張式消聲器右端接口的距離Lz=0.01 m，只考慮入射波為平面波(m=0,s=0)的情況。

2.1 三種方法計算所得的單腔擴張式消聲器傳遞損失的比較

圖2給出了分別用聲學(xué)模態(tài)疊加法、有限元法（LMS Virtural.Lab仿真計算）和基于平面波假定的經(jīng)典公式法計算所得的單腔擴張式消聲器傳遞損失的對比情況。其中膨脹段的截止頻率上下游段的截止頻率4 184Hz，單腔擴張式消聲器的截止頻率

圖2 三種方法計算所得的單腔擴張式消聲器傳遞損失的比較

由圖2可知，對于平面入射波：

（1）聲學(xué)模態(tài)疊加法與有限元方法計算所得的單腔擴張式消聲器的傳遞損失基本一致，所以聲學(xué)模態(tài)疊加法可用于計算單腔擴張式消聲器在各頻段的傳遞損失；

（2）低于單腔擴張式消聲器截止頻率時，起初聲學(xué)模態(tài)疊加法與基于平面波假定的經(jīng)典公式法計算所得的傳遞損失差別非常較小，隨著頻率接近fcut-off，差距增大；

（4）在頻率小于單腔擴張式消聲器的截止頻率fcut-off之前，除了平面波外主要是非傳播的高次波在起作用，且隨著頻率接近fcut-off，峰值向右偏移量越大，峰值的高度也有所增加，在頻率大于fcut-off時，可傳播的行波模式的高次波開始起作用，平面波下的傳遞損失的周期規(guī)律不復(fù)存在。

2.2 膨脹段半徑對傳遞損失的影響

圖3給出了膨脹段半徑R2對單腔擴張式消聲器(L=0.3 m)傳遞損失的影響云圖。其中膨脹段截止頻率最小的單腔擴張式消聲器(R2=0.2 m)的截止頻率fcut-off-R2=3.832c∕2 πR2=1 046Hz ，上下游段的截止頻率fcut-off-R1=3.832c∕2 πR1=1 046Hz ，截止頻率最小的單腔擴張式消聲器的截止頻率fcut-off-R=min(fcut-off-R1,fcut-off-R2)=1 046Hz。

圖3 膨脹段半徑對單腔擴張式消聲器傳遞損失的影響

由圖3可知，單腔擴張式消聲器的其他參數(shù)不變的情況下，對于平面入射波：

（1）低于截止頻率最小的單腔擴張式消聲器的截止頻率前一個零點頻率可用([2Lfcut-off-R/c]c/2L)來估算時，消聲器的傳遞損失隨著膨脹段半徑的增大而增大；

（2）在消聲器上下游段截止頻率內(nèi)，膨脹段半徑對消聲器的傳遞損失的影響較大，但隨著頻率的增大膨脹段半徑對消聲器傳遞損失的影響越來越小，高出一定頻率后，膨脹段半徑對消聲器的傳遞損失基本沒影響；

（3）低于消聲器的截止頻率的峰值頻率不隨膨脹段半徑變化而變化。

2.3 膨脹段寬度對傳遞損失的影響

圖4給出了膨脹段寬度L對單腔擴張式消聲器(R2=0.08 m )的傳遞損失的影響云圖。

圖4 膨脹段寬度對單腔擴張式消聲器傳遞損失的影響

由圖4可知，單腔擴張式消聲器其他參數(shù)不變的情況下，對于平面入射波：

（1）隨著膨脹段寬度的增大，消聲器傳遞損失的峰值向低頻移動；

（2）隨著膨脹段寬度的增大，消聲器傳遞損失最大的頻段向高頻移動。

2.4 膨脹段寬度與膨脹段半徑的比值對傳遞損失的影響

圖5給出了膨脹段寬徑比L/R2對單腔擴張式消聲器(R2=0.08 m，Lz=0.01 m )的傳遞損失的影響曲線。

圖5 膨脹段寬徑比對單腔擴張式消聲器傳遞損失的影響

由圖5可知，當(dāng)L/R2＜0.82時沒有完整的近似周期存在。根據(jù)基于平面波假定的經(jīng)典公式可知，出現(xiàn)峰值的頻率為

所以只有當(dāng)L∕R2＞0.82時才可能出現(xiàn)完整的周期，模態(tài)疊加法的結(jié)果與此吻合。

3 結(jié)語

本文利用聲學(xué)模態(tài)疊加法計算了擴張式消聲器的傳遞損失，并與有限元法和基于平面波假定的經(jīng)典公式法計算結(jié)果進行比較，研究了消聲器膨脹段的尺寸對其傳遞損失的影響。對于平面入射波，研究表明：

（1）基于平面波假定的經(jīng)典公式法只能用于計算單腔擴張式消聲器低于其截止頻率時的傳遞損失，頻率低于時，不但能簡化計算還能得到較好的結(jié)果，聲學(xué)模態(tài)疊加法（入射波可以是模態(tài)波）可用于計算單腔擴張式消聲器各頻段的傳遞損失，且可通過數(shù)值計算來輔助單腔擴張式消聲器降噪的設(shè)計；

（3）隨著單腔擴張式消聲器膨脹段寬度的增大，消聲器傳遞損失峰值向低頻移動，消聲器傳遞損失最大的頻段向高頻移動，可以通過改變膨脹段寬度來使其傳遞損失較大的頻段包含噪聲源的主要頻段，從而達到降噪的目的。

[1]何琳，朱海潮，邱小軍，等.聲學(xué)理論與工程應(yīng)用[M].北京：科學(xué)出版社，2006：145-153.

[2]季振林.消聲器聲學(xué)理論與設(shè)計[M].北京：科學(xué)出版社，2015：14-37.

[3]DAVIS D D,STOKES G M,MOORE D,et al.Theoretical and experimental investigations of mufflers with comments on engine exhaust muffler design[R].NACA,1192.Washington,DC,United States:US Government Printing Office,1954.

[4]MUNJAl M L.Veocity ratio-cum-transfer matrix method for the evaluation of a muffler with mean flow[J].Journal of Sound and Vibration,1975,39(1):105-119.

[5]MILES J.The reflection of sound due to a change in cross section ofacirculartube[J].TheJournalofthe Acoustical Society ofAmerica,1994,16:14-19.

[6]JI Z L,MA Q,ZHANG Z.Application of the boundary element method to prediction acoustic performance of expansion chamber mufflers with mean flow[J].Journal of Sound and Vibration,1994,173(1):57-71.

[7]SELAMET A,JI Z L.Acoustic attenuation performance of circular expansion chambers with offset inlet/outlet:I.Analytical approach[J].Journal of Sound and Vibration,1998,213(4):601-617.

[8]WU C H,WANG C N.Attenuation for the simple expansion chamber muffler with a right angle inlet[J].Journal of Mechanics,2011,27(3):287-292.

[9]宮建國，馬宇山，崔巍升，等.汽車消聲器聲學(xué)特性的傳遞矩陣分析[J].振動工程學(xué)報，2010，23(6)：636-641.

[10]王鍵，陳凌珊.結(jié)構(gòu)參數(shù)對抗性消聲器消聲性能的影響分析[J].小型內(nèi)燃機與摩托車，2010，39(3)：86-88.

[11]詹福良.Virtual.Lab Acoustics聲學(xué)仿真計算從入門到精通[M].西安：西北工業(yè)大學(xué)出版社，2013.

Study on the Transmission Loss of an Expansion Muffler Using Modal Superposition Method

ZHONG Ting-sheng,LI Ya-guang,SHENG Xiao-zhen,XIAO Xin-biao,DU Ji-ping,LIU Mou-kai
(State Key Laboratory of Traction Power,Southwest Jiaotong University,Chengdu 610031,China）

The acoustic modal superposition method is used to establish the trasmission loss calculation model for the single chamber expansion muffler.Its transmission loss is computed through Matlab programming.On this basis,the difference of the acoustic modal superposition method,the finite element method and the classical formula method in calculating the transmission loss of the single chamber expansion muffler is found.The influence of the size of the single chamber expansion muffler on the transmission loss is studied.The research indicates that,for the incident plane wave,the acoustic mode superposition method can be used to calculate the transmission loss of the single chamber expansion muffler in whole frequency range.Increasing the radius of the expanding section can effectively raise the transmission loss in low frequency range,but has small influence on the transmission loss in high frequency range.With increasing of the width of the expansion section,the peak values of the transmission loss shift to low frequency side and the frequency band of the maximum transmission loss moves to the high frequency side.

acoustics;single chamber expansion muffler;acoustic modal superposition;transmission loss;Matlab

TB535+.1

ADIO編碼10.3969/j.issn.1006-1355.2017.06.041

1006-1355(2017)06-0206-05+215

2017-07-01

鐘庭生(1990-)，男，江西省贛州市人，碩士研究生，主要研究方向為軌道交通振動與噪聲。

圣小珍，男，博士生導(dǎo)師。E-mail:shengxiaozhen@hotmail.com