王榮杰 陳虹微

(龍巖學(xué)院機電系，龍巖 364012)

引言

消聲器是控制噪聲的重要舉措，阻抗復(fù)合消聲器又是應(yīng)用廣泛的消聲器之一.人們越來越關(guān)注消聲器在實際工作的使用，國內(nèi)外關(guān)于消聲器研究很多，并已取得不少成果，使消聲器可靠性、經(jīng)濟性不斷提高.聲學(xué)性能和空氣動力性性能是衡量消聲器的兩個重要指標(biāo)，而壓力損失和阻力系數(shù)又是評價消聲器空氣動力性能的兩個重要內(nèi)容.當(dāng)前對消聲器的聲學(xué)性能分析較多，但對空氣動力性性能研究較少［1-3］;對消聲器的傳遞損失研究較多，但對壓力損失分析較少［4-6］;文獻［7］研究了消聲器的聲學(xué)特性，文獻［8］研究了消聲器的傳遞損失，文獻［9］研究了消聲器的壓力損失.文獻［10］研究了消聲器的空氣動力性能.但是對消聲器空氣動力性中的阻力系數(shù)的研究幾乎沒有.本文以阻抗復(fù)合型消聲器為研究對象，分析其結(jié)構(gòu)和特點，建立分析模型，以計算流體力學(xué)(CFD)方法［11-13］，對消聲器的阻力系數(shù)進行研究，重點分析流速、結(jié)構(gòu)等對消聲器空氣動力性能的影響，為消聲器的優(yōu)化設(shè)計提供理論基礎(chǔ).

1 阻抗復(fù)合消聲器結(jié)構(gòu)和特點

阻性消聲器適用于中、高頻噪聲的消除，而抗性消聲器對低、中頻噪聲的消除效果較好.為此，常將二者結(jié)合，可在寬頻帶范圍有良好的消聲，在噪聲控制工程中，阻抗復(fù)合式消聲器應(yīng)用廣泛.

1.1 阻抗復(fù)合消聲器的類型

阻抗復(fù)合消聲器一般由阻性和抗性兩部分組成.根據(jù)工作原理又分為:阻性-擴張室復(fù)合消聲器、阻性-共振腔復(fù)合消聲器、阻性-共振腔-擴張室復(fù)合消聲器、微穿孔板消聲器等.

1.2 阻抗復(fù)合式消聲器的結(jié)構(gòu)

以微穿孔板消聲器為例，是襯裝微穿孔板結(jié)構(gòu)的消聲器.一般是用厚度小于1mm的薄板制作，在薄板上用孔徑小于1mm的鉆頭穿孔，穿孔率為1%一3%.選擇不同的穿孔率和板厚不同的腔深，就可以控制消聲器的頻譜性能，使其在需要的頻率范圍內(nèi)獲得良好的消聲效果.

1.3 阻抗復(fù)合式消聲器的特點

阻力小，能在較寬的頻帶范圍內(nèi)消除氣流噪聲，而且具耐高溫、耐油污、耐腐蝕的性能，即使在氣流中帶有大量水分，也不影響工作.由于穿孔直徑小、板面光滑，因此該類消聲器壓力損失比一般消聲器要?。?4］.

2 消聲器的阻力系數(shù)建模及計算

聲學(xué)性能和空氣動力性性能是評價消聲器的重要指標(biāo).聲學(xué)性能通常用消聲量大小和消聲頻譜特性來表示，主要是傳聲損失和插入損失.空氣動力性能是衡量消聲器好壞的重要指標(biāo)，反映消聲器對氣流阻力的大小，主要是壓力損失和阻力系數(shù).消聲器的壓力損失一般由兩部分組成:一是局部損失;二是沿程損失.局部損失發(fā)生在消聲器內(nèi)收縮、擴張等截面突變處，由于氣流速度因截面變化而變化，導(dǎo)致流體相互碰撞，進一步加劇流體間的相互摩擦，造成損失.與局部損失相對應(yīng)的是局部阻力系數(shù)，其大小取決于消聲器結(jié)構(gòu)形式、管道直徑和氣流速度.沿程損失發(fā)生在氣流通道處，由于消聲器通道壁面摩擦而造成的損失.與沿程損失相對應(yīng)的是沿程阻力系數(shù)，其大小取決于消聲器壁面粗糙度及氣流速度.本文所研究的阻力系數(shù)并不是上述局部阻力系數(shù)和沿程阻力系數(shù)，而是與壓力損失和動壓有關(guān)的一個比值.阻力系數(shù)定義為:通過消聲器前后的壓力損失與氣流動壓之比值.它體現(xiàn)了消聲器對氣流阻力的大小，能全面反映消聲器的空氣動力性能.下面以阻抗復(fù)合消聲器中常見的微穿孔板消聲器為例，利用CFD來計算阻力系數(shù)，并分析空氣流速和消聲器結(jié)構(gòu)對阻力系數(shù)的影響.

2.1 模型基礎(chǔ)

為便于計算和分析，對消聲器的工作條件假設(shè)如下:一是氣體流動為定常流動中的流;二是消聲器物理性能參數(shù)不發(fā)生變化;三是入口速度為勻速;四是不計重力影響.消聲器的物理模型如圖1.

圖1 消聲器模型圖Fig.1 Muffler model

2.2 數(shù)學(xué)模型

根據(jù)定義，阻力系數(shù):

式中a為修正因子;pm為某點處的動壓值.

式中N為測量點的個數(shù);v1為每一個點的速度值，m/s;vm為截面中心處氣流速度值，m/s.

式中p為空氣密度;v為平均空氣速度;g為全力加速度.

［12-13］，控制方程為:

式中φ為通用變量;Γ和S分別為與φ相對應(yīng)的廣義擴散系數(shù)及廣義源項，三者對應(yīng)關(guān)系見表1;V為速度矢量;ρ為空氣密度.

表1 φ、Γ和S對應(yīng)關(guān)系Table 1 The corresponding relations ofφ、Γand S

u、v、w分別為x、y、z三個方向的速度，T為絕對溫度，k為湍流粘性系數(shù)，G為湍流應(yīng)力，Cu、C1、C2、σε、σT和σk分別為量綱一常數(shù)，Pr為普朗特數(shù)，Su、Sv、和Sw分別為x、y、z三個方向的源項.

各常數(shù)取值如下:Cu=0.09，σT=0.95，σk=1.0，C1=1.44，C2=1.92，σε=1.3，g=9.8m/s2，消聲器的外面空氣流速u=v=w=0.

特定截面的全壓為:

式中:ps為截面的氣流靜壓.

消聲器的壓力損失為:

2.3 簡化模型

阻抗復(fù)合消聲器類型較多，簡化模型如圖1.

2.4 基本條件

模型消聲器入口長度L=60mm，入口直徑D=50mm，出口長度L1=60mm，出口直徑D=50mm，擴張腔長度L=200mm，直徑D=100mm，空氣密度為1.23kg/m.邊界條件與實踐工作條件基本相同:輸入流體溫度803K，輸出流體溫度300K，消聲器材料為Q235A，厚度1.5mm，溫度700K.

2.5 計算結(jié)果

利用CFD，將模型的數(shù)值編制成計算機程序上機計算，具體計算過程不闡述，計算結(jié)果見表2.

表2 阻力系數(shù)求解結(jié)果Table 2 Numerical solution of drag coefficient

2.6 實驗驗證

由于消聲器阻力系數(shù)測量比較繁瑣，而阻力系數(shù)又與壓力損失有關(guān)，因此利用壓力損失的計算值和測量值對比來驗證，也能說明問題.為驗證CFD方法的準(zhǔn)確性，參照文獻［15］，分別在距離消聲器出入口5mm的截面上，均勻的選取5個點，以此5個點的全壓值作為出入口的全壓值，兩者之差即為消聲器的壓力損失，入口空氣流速為35m/s.并與計算結(jié)果進行對比，實測值與計算值見表3.從表中看出，二者有一定的誤差，這主要因為實際測量中環(huán)境的影響、儀器本身的精度、人為讀數(shù)的偏差以及CFD計算的理想化與實物的差別等均會造成實測值與計算值的不太吻合，計算值整體小于測量值.但這誤差率在10%以下，在工程中是可以接受的.這也說明利用CFD方法在消聲器動力學(xué)分析中是有效的.

表3 壓力損失計算值與測量值對比表Table 3 Comparison of pressure loss between measurement value and calculation

3 阻力系數(shù)影響因素分析

消聲器必須要有良好的空氣動力性能，特別是安裝消聲器后輸氣是否通暢，對風(fēng)量有無影響，風(fēng)壓有無變化等.而阻力系數(shù)給全面反映消聲器的空氣動力性能.而通過計算發(fā)現(xiàn)，阻力系數(shù)與下列因素有關(guān).

3.1 速度

從表2可以看出阻力系數(shù)隨速度增大而增大.這是因為消聲器的壓力損失隨流體速度的增大而增大.流體在流動過程中，由于沿程損失增大，壓力損失隨之增大.圖2為不同速度下的阻力系數(shù)變化圖.

圖2 不同速度下的阻力系數(shù)變化曲線圖Fig.2 The curve changes of resistance coefficient under different speed

3.2 結(jié)構(gòu)

不同結(jié)構(gòu)的消聲器其阻力系數(shù)大小也不一樣.通過不同模型消聲器阻力系數(shù)的計算，發(fā)現(xiàn)消聲器的結(jié)構(gòu)對阻力系數(shù)有不同影響.

3.2.1 插入管

消聲器可分為無內(nèi)插管和帶內(nèi)插入管單腔消聲器，從計算結(jié)果看，二者的阻力系數(shù)都隨膨脹腔長度的增大而增大.但是，有插入管的消聲器的阻力系數(shù)小于無插入管.這是因為有插入管使出口速度比同位置無插入管的速度高，動壓增加，阻力系數(shù)變小.而無插入管的消聲器其壓力損失大于有插入管的，阻力系數(shù)也大.如圖3.

圖3 有、無內(nèi)插入管消聲器阻力系數(shù)變化曲線圖Fig.3 The curve changes ofresistance coefficient with and without inserting duct in muffler

3.2.2 膨脹腔長度

消聲器的膨脹腔長度不同，其阻力系數(shù)也不一樣.通過計算發(fā)現(xiàn)，膨脹腔長度增大，阻力系數(shù)增大.這是因為膨脹腔長度增大，壓力損失相應(yīng)增大.如圖4.

圖4 阻力系數(shù)隨膨脹腔長度變化曲線圖Fig.4 The curve changes ofresistance coefficient with the length of the inflation cavity

3.2.3 穿孔率

對于微穿孔消聲器，其穿孔率對阻力系數(shù)也有影響.在穿孔孔徑不變的情況下，對橫流微穿孔消聲器，隨著穿孔率的增大，阻力系數(shù)變小.這是因為，穿孔率逐漸增大，穿孔處流通面積增大，氣流速度低，穿孔造成的壓力損失變小，阻力系數(shù)隨之減小，如圖5.對直流微穿孔消聲器，隨著穿孔率的增大，阻力系數(shù)變大.這是因為穿孔率變大，進入消聲器膨脹腔的氣體加多，氣流的流程增加，沿程壓力損失變大，總的壓力損失變大，從而阻力系數(shù)變大，如圖6.

圖5 橫流穿孔消聲器阻力系數(shù)隨穿孔率變化曲線圖Fig.5 The change curve of resistance coefficient of cross-flow perforated muffler with perforation rate

圖6 直流穿孔消聲器阻力系數(shù)隨穿孔率變化曲線圖Fig.6 The change curve of resistance coefficient of direct-flow perforated muffler with perforation rate

3.2.4 穿孔直徑

穿孔直徑對阻力系數(shù)也有影響.在穿孔率不變的情況下，對橫流微穿孔直徑變大，流速變小，壓力損失變小，如圖7.對直流穿孔消聲器增加，阻力系數(shù)變大.這是因為穿孔直徑變大，進入消聲器膨脹腔的氣體增加，沿程壓力損失增大，總的壓力損失增加，阻力系數(shù)變大，如圖8.但這種變化范圍較小，也就是相對橫流微穿孔消聲器，變化不大明顯.

圖7 橫流穿孔消聲器阻力系數(shù)隨穿孔直徑變化曲線圖Fig.7 The change curve of resistance coefficient of cross-flow perforated muffler with perforation diameter

圖8 直流穿孔消聲器阻力系數(shù)隨穿孔直徑變化曲線圖Fig.8 The change curve of resistance coefficient of direct-flow perforated muffler with perforation diameter

3.3 內(nèi)部流場

消聲器的內(nèi)部流場主要表現(xiàn)在內(nèi)部氣流的速度場，壓力場的分布特點以及相互關(guān)系.

借鑒文獻［15］，消聲器的速度場與壓力場具有相似的分布規(guī)律，消聲器煙筒內(nèi)的氣流速度大于入口流速，而入口流速又高于出口流速.入口處的氣流速度對消聲器的阻力系數(shù)影響較大.

4 應(yīng)用實例

以龍巖機械廠生產(chǎn)的30D叉車消聲器為例.利用CFD計算阻力系數(shù)，分析消聲器的性能，找出不足和原因，提出相應(yīng)的改進方法，優(yōu)化消聲器的空氣動力性能，再次通過CFD計算加以驗證，把理論研究轉(zhuǎn)化為實際應(yīng)用.

4.1 物理模型

LG30D消聲器是以微穿孔為主的復(fù)合消聲器，其三個腔長度分別為:前腔335mm、中腔300mm、后腔275mm.膨脹長度200mm，直徑420mm，內(nèi)插入管長管為160mm，短管為40mm，穿孔管內(nèi)徑100mm，穿孔個數(shù)20個，穿孔直徑10mm，穿孔率9%.消聲器結(jié)構(gòu)如圖9.

圖9 消聲器結(jié)構(gòu)圖Fig.9 Muffler structure

4.2 阻力系數(shù)計算

該叉車發(fā)動機排量3L，轉(zhuǎn)速2600r/min，消聲器入口速度30m/s，利用CFD進行計算，阻力系數(shù)為1.69，壓力損失2932Pa.

4.3 動力性能分析

當(dāng)氣體從進氣管流入前腔，經(jīng)插入管到中腔、后腔，最后經(jīng)出氣管排出，途中氣體通流面積產(chǎn)生變化，造成渦流現(xiàn)象，產(chǎn)生能量消耗，導(dǎo)致壓力損失.從該消聲器相關(guān)圖表中可以看出，插入管和穿孔區(qū)域是整個消聲器壓力損失最為嚴(yán)重的部位，也是能量消耗最為集中的部分.消聲器每個腔的壓力有差別.從進氣開始到排氣，順著氣流方向，各腔的壓力逐漸變小，如前腔氣體壓力范圍為2500～5000Pa，中腔壓力范圍為1500～1900Pa，后腔壓力范圍為700～1000Pa.壓力損失在前腔較大，為3342Pa，主要集中在有穿孔結(jié)構(gòu)的區(qū)域;從內(nèi)部流場選取幾個穿孔位置分析，發(fā)現(xiàn)穿孔的入口氣流速度對阻力系數(shù)影響較大.也就是入口速度高，阻力系數(shù)大，壓力損失大，反之，變小.因此，降低消聲器入口處的速度，可降低壓力損失和阻力系數(shù)，也可降低氣流對消聲器氣體和穿孔管的沖擊，提高消聲器的使用壽命和消聲性能.另外，穿孔直徑和穿孔的配置對阻力系數(shù)和壓力損失影響較大，這也驗證3.2.3和3.2.4所述.

該消聲器阻力系數(shù)較小，壓力損失不大，這在實際使用中已得到驗證.但仔細分析，還有一些不足.如消聲器中長插入管偏長，已進入到后腔，氣流緩沖、回轉(zhuǎn)路徑偏短，不利于氣體流動.前腔與中腔長度較接近，消聲頻率基本一致，不利于消聲.穿孔直徑和配置還可以進一步優(yōu)化.因此，在保證消聲的前提下，對消聲器進行改進:一是將三個的長度分別調(diào)整為285mm、315mm、280mm.二是縮短長插入管的長度，使其插入到后腔的出口長度為零.三是將原先的10mm，20排穿孔改為8mm，10排的結(jié)構(gòu).再對改進后的消聲器進行CFD仿真計算，其阻力系數(shù)為1.13，下降了33%;壓力損失為2150Pa，下降了27%.三個腔的速度也比較均勻，消聲器的消聲性能和空氣動力性能更好.

5 結(jié)論

(1)利用CFD方法對消聲器進行計算和分析是可行的，可以為消聲器設(shè)計和優(yōu)化提供理論基礎(chǔ).

(2)阻力系數(shù)能全面反映消聲器的空氣動力性能.消聲器的阻力系數(shù)小，壓力損失小，空氣動力性能好.

(3)相同條件下，空氣流速增加，阻力系數(shù)增大，壓力損失增加.

(4)消聲器的結(jié)構(gòu)對阻力系數(shù)的影響較大.有插入管的消聲器的阻力系數(shù)小于無插入管.膨脹腔長度增大，阻力系數(shù)增大.

(5)對于微穿孔消聲器其穿孔直徑和穿孔率對阻力系數(shù)有很大影響.穿孔的優(yōu)化配置可以使消聲器有比較好的空氣動力性能.

(6)消聲器的速度場與壓力場具有相似的分布規(guī)律，消聲器煙筒內(nèi)的氣流速度大于入口流速，而入口流速又高于出口流速.入口處的氣流速度對消聲器的阻力系數(shù)影響較大.

參考文獻

1 Dokumaci E.Matrizant approach to acoustic analysis of perforated multiple pipe mufflers carrying mean flow.Journal of Sound and Vibration，1996，191(4):505～518

2 Wu T W，Zhang P.Boundary element analysis of mufflers with an improved method for deriving the four-pole parameters.Journal of Sound and Vibration，1998，217(4):767～77

3 Tsuji T，Tsuchiya T，Kagawa Y.Finite element and bounary element modeling for the acoustic wave transmission in mean flow medium.Journal of Sound and Vibration，2002，255(5):849～866

4 Munjal M L，Behera B K，Thawani P T.Transfer matrix model for the reverse-flow，three-duct，open end perforated element muffler.Applied Acoustics，1998，54(3):229～238

5 Yeh L J，Chang Y C，Chiu M C，et al.GA Optimiza-tion on multi-segments muffler under space con-straints.Applied Acoustics，2004，65(5):521～543

6 Selamet A，Denia F D，Besa A J.Acoustic behavior of circular dual-chamber mufflers.Journal of Sound and Vibration，2003，265(5):967～985

7 季振林.直通穿孔消聲器聲學(xué)性能計算及分析.哈爾濱工程大學(xué)學(xué)報，2005，26(3):301～306(Ji Z L.Calculation and analysis on acoustic performance of the straightthrough perforated muffler.Journal of Harbin Engineering University，2005，26(3):301～306(in Chinese)).

8 王劍，王伯良，周建立等.穿孔插入式消聲器消聲頻率特性的四端網(wǎng)絡(luò)法計算.洛陽工學(xué)院學(xué)報，1999，20(2):70～74(Wang J，Wang B L，Zhou L L，et al.Computing with four-terminal network on perforated-insert muffler of noise elimination frequency characteristic.Journal of Luoyang Institute of Technology，1999，20(2):70～74(in Chinese))

9 方建華，周以齊，胡效東等.內(nèi)插管抗性消聲器的CFD仿真及壓力損失研究.武漢理工大學(xué)學(xué)報，2009，33(4):795～798(Fang J H，Zhou Y Q，Hu X D，et al.CFD computatin of pressure loss of insert-duct resistance mufflers.Journal of Wuan University of Science and Technology，2009(4):795～798(in Chinese))

10 Hu X D，Zhou Y Q，F(xiàn)anf J H，et al.Computational fluid dynamics research on pressure loss of cross-flow perforated muffler.Chinese Journal of Mechanical Engineering，2007，20(2):88～93

11 Tamer WM，Ahmed K N.Vinsualization of CFD results in immersive virtual environment.Advances in Engineering Software，2001，32(9):717～730.

12 EI-Sayed M，Sun T，Berry J.Shape optimization with computational fluid dynamics.Advances in Engineering Software，2005，36(9):607～613.

13 ShaoL，Riffat S B.Accuracy of CFD for predicting pressure loss in HVACduct fittings.Applied Energy，1995，51(3):233～248

14 顧強.噪聲控制工程.北京:煤炭工業(yè)出版社，2002(Gu Q.Noise control engineering.Beijing:Coal Industry Publishing House，2002(in Chinese))

15 胡效東，周以齊，方建華.單雙腔抗性消聲器壓力損失CFD研究.中國機械工程，2006，17(24):2567～2572(Hu X D，Zhou Y Q，F(xiàn)ang J H.CFD computation of pressure loss of single and dual-chamber resistance mufflers.Chinese Journal of Mechanical Engineering，2006，17(24):2567～2572(in Chinese))