張哲，王小鵬，陳天寧，奚延輝

（1.西安交通大學(xué)機(jī)械工程學(xué)院，710049，西安；

2.西安交通大學(xué)機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室，710049，西安）

金屬多孔材料在新型汽車、高速列車、航空航天器等高能耗裝備的消聲降噪技術(shù)中均具有重要的應(yīng)用價(jià)值［1－3］。目前，多孔金屬材料的吸聲性能研究大多在常溫常壓條件下進(jìn)行，而對高溫或材料前、后端面間有溫度梯度條件下的吸聲性能研究則相對不足［4］，國內(nèi)外已發(fā)表的相關(guān)研究成果寥寥無幾，其中一個(gè)主要原因是缺少相應(yīng)的測試技術(shù)。由于目前市場上沒有可以工作在高溫環(huán)境（高于150℃）的揚(yáng)聲器，因此在用阻抗管測量吸聲材料高溫吸聲系數(shù)時(shí)，總是揚(yáng)聲器處于低溫，而吸聲材料處于高溫，從而在阻抗管內(nèi)的空氣中形成了溫度梯度。傳統(tǒng)的傳遞函數(shù)要求阻抗管內(nèi)空氣溫度均勻［5］，因此當(dāng)阻抗管內(nèi)存在溫度梯度時(shí)，就必須考慮能否繼續(xù)應(yīng)用傳遞函數(shù)法測量吸聲系數(shù)。馮凱、薛昌意等分別設(shè)計(jì)了高溫、溫度梯度阻抗管［6－7］，但是由于管徑較小、管長較大且一端封閉，使溫度測量較為困難，導(dǎo)致測量結(jié)果誤差較大，誤認(rèn)為管內(nèi)溫度符合傳遞函數(shù)法的均溫要求，因而沒有對阻抗管內(nèi)部能否產(chǎn)生平面波進(jìn)行驗(yàn)證。本文利用自制的溫度傳感器支架得到阻抗管內(nèi)的準(zhǔn)確溫度，對自制阻抗管內(nèi)能否產(chǎn)生平面波進(jìn)行校核，同時(shí)考慮軸向溫度梯度的影響對傳統(tǒng)傳遞函數(shù)法進(jìn)行修正，以期為高溫環(huán)境下傳遞函數(shù)法的使用提供理論基礎(chǔ)。

1 高溫阻抗管中平面波的產(chǎn)生條件

只有當(dāng)阻抗管內(nèi)傳播平面波時(shí)，才可以用傳聲器測點(diǎn)處的聲壓值來代表整個(gè)波的能量，否則每點(diǎn)聲壓值不同，將使得測量失去意義。因此，傳遞函數(shù)法應(yīng)用于高溫阻抗管的要求依然是管中要能產(chǎn)生平面波。現(xiàn)已知阻抗管中的空氣有溫度梯度，那么該溫度梯度能否產(chǎn)生平面波？這是本文首先要解決的問題。首先對自制的高溫阻抗管進(jìn)行溫度場仿真，并與實(shí)驗(yàn)測量值進(jìn)行比較，得到一個(gè)準(zhǔn)確的溫度場，然后將該溫度場導(dǎo)入comsol軟件中進(jìn)行仿真，以驗(yàn)證在該溫度場下阻抗管中能否產(chǎn)生平面波，并且探索產(chǎn)生平面波的條件。

1.1 高溫阻抗管溫度場分析

自制高溫阻抗管的原理圖如圖1［7］所示。該裝置是在傳統(tǒng)阻抗管上集成了加熱裝置的結(jié)果。從圖1可以看出，加熱方案是在材料后端面加熱（右熱源），在材料前端面周向進(jìn)行水冷，以便在材料兩端形成溫度梯度，并且在材料前端一段阻抗管上進(jìn)行加熱保溫（左熱源）。本文僅考慮材料在高溫環(huán)境下材料前端空氣的溫度分布，并未考慮材料兩端溫度梯度的實(shí)現(xiàn)，因此在仿真和實(shí)驗(yàn)中均不考慮水冷。

圖1 自制高溫阻抗管原理圖

阻抗管的內(nèi)徑為50mm，壁厚為3mm，多孔材料厚度為50mm，較近的傳聲器距多孔材料前端面50mm，較遠(yuǎn)的傳聲器距材料前端面80mm，其他結(jié)構(gòu)參數(shù)可參考文獻(xiàn)［7］。分析可知阻抗管內(nèi)部的空氣產(chǎn)生了自然對流換熱，因此，采用Fluent軟件對阻抗管內(nèi)的溫度場進(jìn)行仿真分析。首先對模型進(jìn)行簡化，由于對稱性取阻抗管的一半進(jìn)行分析，用ICEM CFD軟件對其進(jìn)行結(jié)構(gòu)網(wǎng)格劃分，網(wǎng)格數(shù)為65 529，并給定各邊界條件，最終的計(jì)算模型如圖2所示，各部件的物性參數(shù)見表1，阻抗管中的空氣為不可壓縮理想氣體。

圖2 Fluent計(jì)算模型

表1 高溫阻抗管各部件的物性參數(shù)

重力取9.8m/s2，選擇k－ε湍流模型和分離隱式求解器并進(jìn)行計(jì)算。定義阻抗管中空氣柱的中軸線，以及距離中軸線上、下、左、右均為16mm的4條直線的位置來對比仿真與實(shí)驗(yàn)的結(jié)果。仿真與實(shí)驗(yàn)（300℃）所得中軸線上和中軸線下的溫度分布結(jié)果如圖3所示。

從圖3中的實(shí)驗(yàn)數(shù)據(jù)可以看出：管內(nèi)空氣在加熱后產(chǎn)生了自然對流，溫度高的氣體密度小，在浮力的作用下處在較高的位置，溫度較低的氣體則處在較低的位置，在阻抗管中的上、下位置間存在較大的溫度差，即徑向溫度梯度，其中材料前端面到傳聲器1這一段主要是軸向溫度差（約為60℃），徑向溫度差較小。由圖3可見，仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)吻合較好，這為下一步仿真聲場做好了準(zhǔn)備。

圖3 溫度場的仿真與實(shí)驗(yàn)結(jié)果對比

1.2 高溫下的阻抗管聲場分析

將Fluent所得到的空氣溫度場導(dǎo)入comsol軟件中，計(jì)算區(qū)域?yàn)樽杩构軆?nèi)的圓柱空氣柱。通過行波來研究聲波在空氣軸向溫度梯度下是否滿足平面波的傳播特性，將2個(gè)圓柱端面設(shè)為平面波輻射邊界條件，即平面波無反射，采用單頻信號，入射波幅值為1Pa，管壁為硬邊界條件，頻率分別設(shè)為800Hz和4 800Hz，管壁為硬邊界條件。仿真結(jié)果見圖4、圖5，圖中同時(shí)畫出了中軸線上、下的溫度差，用來代表徑向溫度梯度的大小。

圖4 入射聲波頻率為800Hz時(shí)的管內(nèi)聲場分布

圖5 入射聲波頻率為4 800Hz時(shí)的管內(nèi)聲場分布

從圖4和圖5可以看出：低頻時(shí)溫度梯度對聲波的影響只體現(xiàn)在聲壓幅值的變化上，管內(nèi)仍為平面波，體現(xiàn)了低頻聲波的越障性；當(dāng)入射波為高頻時(shí)，在徑向溫度梯度大的地方，如在加熱區(qū)域和揚(yáng)聲器之間，同一截面不同位置的聲壓及其幅值產(chǎn)生了較大的變化，不再是平面波，而在徑向溫度梯度小的地方，如傳聲器附近，則產(chǎn)生了近似平面波。對于阻抗管中遠(yuǎn)離傳聲器的非平面波，可以認(rèn)為是產(chǎn)生平面波的聲源部分，故只需要考慮2個(gè)傳聲器處的平面波。因此，可將阻抗管內(nèi)空氣的均溫要求擴(kuò)寬到僅有軸向溫度梯度。

由于實(shí)驗(yàn)測得材料前端空氣幾乎只有軸向溫度梯度，因此可認(rèn)為自制裝置滿足產(chǎn)生平面波的要求，能夠使用傳遞函數(shù)法計(jì)算吸聲系數(shù)。為了考慮軸向溫度梯度對傳遞函數(shù)法的影響，需要在該梯度下對傳遞函數(shù)法進(jìn)行修正。

2 空氣有軸向溫度梯度時(shí)對傳遞函數(shù)法的修正

2.1 修正公式推導(dǎo)

文獻(xiàn)［8－10］研究了在介質(zhì)的密度和聲速僅沿一個(gè)方向變化時(shí)平面波的解析解。但是，文獻(xiàn)［8］中推導(dǎo)解析公式時(shí)采用了均溫波導(dǎo)管理論［11］中的處理辦法，將公式（8）中的f（z）當(dāng)作常項(xiàng)，因而得到的結(jié)果是錯(cuò)誤的。與文獻(xiàn)［9］相比，文獻(xiàn)［10］給出的方法更為簡單易懂——只需要通過溫度場計(jì)算出密度場和聲速場，然后用選定的密度和聲速擬合函數(shù)進(jìn)行擬合，就可以得到一個(gè)簡單的解析解。密度擬合函數(shù)和聲速擬合函數(shù)分別為

式中：Q、a、α、c1、δ均為待擬合參數(shù)。本文采用文獻(xiàn)［10］中的方法，最終獲得的空氣軸向溫度梯度下的平面波解析解為

式中：A、B為艾里函數(shù)；η為位置x的函數(shù)；C、D為與溫度場有關(guān)的常數(shù)。將正向波、負(fù)向波分別表示為

式中：C1、D1、C－1、D－1均為與溫度場有關(guān)的常數(shù)。由此，可寫出空氣中駐波的形式為

式中：r為反射系數(shù)。從傳聲器2處到傳聲器1處的傳遞函數(shù)為

從而可求出反射系數(shù)的修正公式為

2.2 仿真驗(yàn)證

在comsol中分別用聲強(qiáng)法、傳遞函數(shù)法以及修正傳遞函數(shù)法計(jì)算某多孔材料的吸聲系數(shù)。聲強(qiáng)法是通過進(jìn)入材料的聲強(qiáng)與總聲強(qiáng)之比來計(jì)算吸聲系數(shù)。利用comsol計(jì)算全場節(jié)點(diǎn)的聲壓及振速，進(jìn)而計(jì)算出實(shí)際的聲強(qiáng)，將其與入射平面波的總聲強(qiáng)（設(shè)定值）相比可計(jì)算出吸聲系數(shù)。仿真模型如圖6所示。多孔材料的參數(shù)值如下：孔隙率為0.8；流阻率為10.5Pa·s/m2；黏性特征長度參數(shù)為0.49；熱特征長度為470μm；黏性特征長度為240μm；曲折度為1.005 9。采用Johnson－Champoux－Allard等效流體模型計(jì)算其等效密度和等效聲速。

圖6 聲場仿真模型

空氣的密度及聲速用實(shí)驗(yàn)測得的材料前端到較遠(yuǎn)傳聲器之間的溫度進(jìn)行計(jì)算［12］，用式（1）、式（2）進(jìn)行擬合可得到相應(yīng)的曲線，如圖7所示。

圖7 密度和聲速擬合結(jié)果

設(shè)多孔材料后端面及圓柱面為硬邊界條件，空氣柱另一端面為平面波入射，采用單頻信號，幅值為1Pa。仿真結(jié)果如圖8所示。

圖8 吸聲系數(shù)計(jì)算結(jié)果對比

從圖8可知，聲強(qiáng)法和修正傳遞函數(shù)法所得結(jié)果完全重合，證明了修正傳遞函數(shù)法的正確性。同時(shí)，修正與未修正傳遞函數(shù)法所得曲線的差異很小，說明材料前端空氣有較小軸向溫度梯度（本例為60℃/0.140m≈428.6℃/m）對吸聲系數(shù)的影響不大。此結(jié)論從另一方面也說明，在設(shè)計(jì)高溫阻抗管時(shí)，阻抗管前端空氣有軸向溫度梯度的工況是允許的，可由上述修正后的公式計(jì)算吸聲系數(shù)。

3 結(jié) 論

本文對使用傳遞函數(shù)法時(shí)的均溫要求做了推廣，允許在傳聲器測量段空氣僅有軸向溫度梯度，并給出了修正公式。研究結(jié)果表明，當(dāng)材料前端空氣軸向溫度梯度較小時(shí)，其對吸聲系數(shù)計(jì)算的影響可以忽略不計(jì)。這一結(jié)論擴(kuò)寬了高溫阻抗管的設(shè)計(jì)空間，具有一定的潛在應(yīng)用價(jià)值。

［1］盧天健，何德坪，陳常青，等.超輕多孔金屬材料的多功能特性及應(yīng)用 ［J］.力學(xué)進(jìn)展，2006，36（4）：517－535.LU Tianjian，HE Deping，CHEN Changqing，et al.The multi－functionality of ultra－light porous metals and their applications［J］.Advances in Mechanics，2006，36（4）：517－535.

［2］彭鋒，王曉林，孫艷，等.高聲壓級時(shí)多孔金屬板的吸聲特性研究 ［J］.聲學(xué)學(xué)報(bào)，2009，24（3）：266－274.PENG Feng，WANG Xiaolin，SUN Yan，et al.Investigation on the sound absorption characteristics of porous metal plates at high sound pressure levels［J］.Acta Acustica，2009，24（3）：266－274.

［3］張波，陳天寧，馮凱，等.燒結(jié)金屬纖維多孔材料的高溫吸聲性能 ［J］.西安交通大學(xué)學(xué)報(bào)，2008，42（11）：1327－1331.ZHANG Bo，CHEN Tianning，F(xiàn)ENG Kai，et al.Sound absorption properties of sintered fibrous metals under high temperature conditions ［J］.Journal of Xi’an Jiaotong University，2008，42（11）：1327－1331.

［4］張波.多孔金屬材料吸聲性能研究 ［D］.西安：西安交通大學(xué)，2008：3－8.

［5］全國聲學(xué)標(biāo)準(zhǔn)化技術(shù)委員會(huì).GB/T 18696.2－2002聲學(xué)阻抗管中吸聲系數(shù)和聲阻抗的測量：第2部分傳遞函數(shù)法［S］.北京：中國標(biāo)準(zhǔn)出版社，2002.

［6］馮凱.多孔金屬材料高溫吸聲性能實(shí)驗(yàn)裝置及吸聲性能研究 ［D］.西安：西安交通大學(xué)，2008.

［7］薛昌意.變參數(shù)溫度場下多孔金屬材料吸聲性能研究［D］.西安：西安交通大學(xué)，2010.

［8］田曉培，金韜.聲波在非均勻介質(zhì)波導(dǎo)中的傳播 ［J］.復(fù)旦學(xué)報(bào)：自然科學(xué)版，2011，50（4）：527－532.TIAN Xiaopei，JIN Tao.Acoustic wave propagates in nonuniform medium waveguides［J］.Journal of Fudan University：Natural Science，2011，50（4）：527－532.

［9］SUJITH R I，WALDHERR G A，ZINN B T，An exact solution for one－dimensional acoustic fields in ducts with an axial temperature gradient ［J］.Journal of Sound and Vibration，1995，184（3）：389－402

［10］Robins A J.Reflection of a plane wave from a fluid layer with continuously varying density and sound speed［J］.Journal of the Acoustical Society of America，1991，89（4）：1686－1695.

［11］杜功煥，朱哲民，龔秀芬.聲學(xué)基礎(chǔ) ［M］.3版.南京：南京大學(xué)出版社，2012：190－191.

［12］S A′NCHEZ－DEHESA J，ANGELOV M I，CERVERA F，et al.Sound control by temperature gradients［J］.Applied Physics Letters，2009，95：204102（1－3）.

［本刊相關(guān)文獻(xiàn)鏈接］

張健，趙桂平，盧天健.泡沫金屬在沖擊載荷下的能量吸收特性.2013，47（11）：105－112.［doi：10.7652/xjtuxb201311019］谷偉，張虎，李增耀，等.混合氣體在典型多孔介質(zhì)內(nèi)擴(kuò)散過程的數(shù)值模擬.2012，46（3）：107－112.［doi：10.7652/xjtuxb 201203019］

張虎，谷偉，李增耀，等.氧化硅納米多孔材料表面的水蒸氣吸附動(dòng) 力 學(xué) 研 究.2012，46（3）：1－8.［doi：10.7652/xjtuxb 201203001］

于雪梅.氣體靜壓軸承用多孔SiC陶瓷的制備及靜態(tài)性能.2011，45（3）：117－120.［doi：10.7652/xjtuxb201103022］

張健，趙桂平，盧天健.閉孔泡沫鋁應(yīng)變率效應(yīng)的試驗(yàn)和有限元分 析.2010，44（5）：97－101.［doi：10.7652/xjtuxb201005 020］

張波，陳天寧.高聲壓激勵(lì)下多孔金屬材料的吸聲性能數(shù)值計(jì)算.2010，44（3）：58－62.［doi：10.7652/xjtuxb201003012］

屈治國，徐治國，陶文銓，等.通孔金屬泡沫中的空氣自然對流傳熱 實(shí) 驗(yàn) 研 究.2009，43（1）：1－4.［doi：10.7652/xjtuxb 200901001］

張波，陳天寧，馮凱，等.燒結(jié)金屬纖維多孔材料的高溫吸聲性能.2008，42（11）：1327－1331.［doi：10.7652/xjtuxb2008 11003］

張波，陳天寧.燒結(jié)金屬纖維材料的吸聲模型研究.2008，42（3）：328－332.［doi：10.7652/xjtuxb200803016］

DUPèRE Iain D J，盧天健，DOWLING Ann P.聲學(xué)多孔材料的孔結(jié)構(gòu)優(yōu)化.2007，41（11）：1251－1256.［doi：10.7652/xjtuxb200711001］

DOWLING Ann P，盧天健，DUPèRE Iain D J.多孔泡沫材料的聲 吸收 特性.2007，41（9）：1003－1011.［doi：10.7652/xjtuxb200709001］