劉林芽，呂 銳，許群峰，邱仕輝

(1．華東交通大學(xué)鐵路環(huán)境振動與噪聲教育部工程研究中心，江西南昌330013;2．廣州市地下鐵道總公司，廣東廣州510310)

隨著中國鐵路和城市軌道交通的飛速發(fā)展，列車運行的速度不斷提高，而由此帶來的環(huán)境與安全問題也日益突出。其中高速鐵路輪軌噪聲就嚴(yán)重干擾了人們的正常生活、工作和身體健康［1］。每一種噪聲均是由相對應(yīng)結(jié)構(gòu)的振動所產(chǎn)生，而輪軌噪聲產(chǎn)生的原因則是車輪和鋼軌之間相互作用產(chǎn)生的振動，隨之向外輻射出聲波而形成的。輪軌噪聲按產(chǎn)生機理可分為三種噪聲:滾動噪聲、沖擊噪聲、摩擦噪聲。輪軌滾動噪聲的產(chǎn)生是［2］因為軌道系統(tǒng)中鋼軌表面的短波不平順激發(fā)，當(dāng)車輪通過此不平順的鋼軌時，車輪將產(chǎn)生振動，振動通過空氣傳播而產(chǎn)生噪聲，產(chǎn)生機理如圖1所示。沖擊噪聲是車輪經(jīng)過鋼軌接縫處或表面呈波紋狀的鋼軌、車輪轍叉、車輪扁疤所產(chǎn)生的噪聲;摩擦噪聲也稱尖嘯聲，是當(dāng)列車通過小半徑曲線鋼軌時產(chǎn)生的高音調(diào)噪聲。

圖1 輪軌滾動噪聲產(chǎn)生機理及其傳播途徑Fig．1 The rolling wheel/rail noise generation mechanism and the transmission way

在高速鐵路中，輪軌滾動噪聲是高速鐵路的主要噪聲，而滾動噪聲實際上就是輪軌之間的高頻振動造成的［3］。在軌道系統(tǒng)的聲輻射問題上，國外一些學(xué)者如Remington［4］從輪軌相互作用角度出發(fā)，在輪軌系統(tǒng)動力特性分析的基礎(chǔ)上，建立了輪軌滾動噪聲預(yù)測模型;緊接著Thompson在Remington的基礎(chǔ)上建立了軌道聲輻射二維模型，他利用有限元法和邊界元法，研究了車輪與軌道的聲輻射特性［5－6］;Nielsen［7］利用邊界元法預(yù)測軌枕聲輻射。國內(nèi)學(xué)者中，徐志勝［8］運用車輛 －軌道耦合動力學(xué)與聲輻射理論，建立了基于時域的軌道結(jié)構(gòu)噪聲預(yù)測模型;吳天行［9］利用有限元法和邊界元法研究了鋼軌的振動模態(tài)及其聲輻射特性;方銳［10］結(jié)合有限元和邊界元方法對在不同軌道結(jié)構(gòu)參數(shù)下，軌道結(jié)構(gòu)的聲輻射特性;另外，國內(nèi)還有其他學(xué)者［11－13］也對輪軌高頻振動與噪聲的理論計算進(jìn)行了大量研究。

文章以高速鐵路無砟軌道為研究對象，建立了軌道系統(tǒng)垂向高頻振動聲輻射計算模型，運用有限元和邊界元相結(jié)合的方法，將文獻(xiàn)［12］計算出軌道系統(tǒng)的垂向高頻振動響應(yīng)作為聲輻射計算邊界條件，研究線路見圖2，最后可以得到鋼軌與軌下結(jié)構(gòu)的聲輻射特性。

1 輪軌噪聲預(yù)測直接邊界元法

邊界元法是在有限元法之后發(fā)展起來的一種較精確有效的工程數(shù)值分析方法，它與有限元法相比，由于降低了問題的維數(shù)，而顯著降低了自由度數(shù)，所以邊界元法被公認(rèn)為在分析結(jié)構(gòu)聲學(xué)數(shù)值問題中有巨大的優(yōu)勢［14－16］。

對于聲學(xué)波動方程，由傅里葉變換可知，隨時間變化的任意振動都可以看作是多個簡諧振動的疊加或積分，可以推導(dǎo)出鋼軌聲輻射的Helmholtz方程，其微分方程式為:

式中:▽2為拉普拉斯算子;為聲波波數(shù);c為空氣介質(zhì)中的聲速;ω為角頻率;p為聲壓。

針對鋼軌輻射外聲場，采用直接邊界元法，將鋼軌結(jié)構(gòu)表面S離散成M個邊界單元，其節(jié)點數(shù)為K，1個單元上的節(jié)點數(shù)為 L，p′和 v′n表示節(jié)點L的聲壓和法向振速。則鋼軌表面單元上任意一點的聲壓和法向振速可用單元節(jié)點的聲壓和法向振速表示為:

圖2 研究線路圖Fig．2 The circuit diagram

式中:Nl(ξ，η)為插值函數(shù)，依次將每個節(jié)點作為源點，即可得直接邊界元法的求解方程，其表達(dá)式為:

式中:［A］、［B］為K×K階系數(shù)矩陣;{}p 為結(jié)構(gòu)(鋼軌)表面節(jié)點聲壓向量;{vn}為結(jié)構(gòu)(鋼軌)表面節(jié)點法向速度向量。通過求解方程(3)就可以求得鋼軌結(jié)構(gòu)表面節(jié)點聲壓向量{p}。

對于空間聲場中任意一點r(r?S)的聲壓，可由鋼軌表面上的聲壓 {p}和法向振動速度{vn}積分得到:

式中:{a}和為插值系數(shù)向量。

2 軌道系統(tǒng)邊界元計算模型

由于模型結(jié)構(gòu)和邊界條件的對稱性，在軌道中心設(shè)置一個對稱面來研究軌道系統(tǒng)的聲輻射情況。利用邊界元來計算鋼軌的噪聲輻射量時，邊界元網(wǎng)格對最后結(jié)果的求解有很大影響，模型選取越長，網(wǎng)格劃分的越細(xì)，數(shù)目越多，計算時對計算機的要求就越高［17］。

采用直接邊界元法對軌道系統(tǒng)進(jìn)行聲輻射特性分析，由于混凝土底座對軌道系統(tǒng)聲輻射貢獻(xiàn)很少，所以，軌下結(jié)構(gòu)只考慮軌道板，重新劃分網(wǎng)格并計算滿足精度要求后，邊界元模型如圖3所示。

圖3 軌道系統(tǒng)邊限元計算模型Fig．3 Transient boundary element model of the rail system

通過軌道系統(tǒng)輻射聲場云圖可以非常直觀地形象地看出軌道系統(tǒng)周圍空間聲輻射特性的輻射情況與變化情況。為了了解軌道系統(tǒng)的振動聲輻射特性，在軌道系統(tǒng)附近位置建立了1個面場點和4個點場點。其中面場點位于軌道系統(tǒng)的左方，與列車運行方向平行，距離鋼軌中心線1．0 m，其平行軌道方向長度為5．0 m，高度為2．0 m，如圖4(a)和4(b)所示;另4個點場點為高出軌面1．5 m布置，距離鋼軌的距離依次為7．5，15，30和60 m，如圖4(c)所示。

應(yīng)用降脂湯結(jié)合針灸治療高脂血癥127例患者，其療效為顯效92例，有效28例，無效及惡化7例?？傆行?4.5%。

圖4 場點布置圖Fig．4 Layout drawings of field point

3 計算結(jié)果與討論

將文章［12］中計算得到的振動響應(yīng)結(jié)果作為位移邊界條件，利用直接邊界元計算軌道系統(tǒng)的輻射噪聲，通過聲輻射效率、輻射聲功率及輻射聲場云圖對此進(jìn)行綜合分析。計算中取空氣密度ρ=1．21 kg/m3，空氣中聲速 c=340 m/s。

經(jīng)過計算，分別得到了鋼軌與軌道板的聲輻射效率，如圖5和圖6所示。

由圖5可知:隨著頻率的增大，鋼軌的聲輻射效率也越來越大，這是因為聲波的頻率越高，波長也越來越小，輻射的能量相對越大。在800－3000 Hz內(nèi)聲輻射效率有震蕩趨勢，在2 000 Hz時出現(xiàn)1個明顯的峰值，說明此時鋼軌振動引起的聲輻射效率較高，高頻時，聲輻射效率近似為1。由圖6可知:當(dāng)?shù)陀?00 Hz時，軌道板的聲輻射效率幾乎成線性增長，在500 Hz時，軌道板的聲輻射效率已經(jīng)達(dá)到1，說明在這個頻率軌道板的聲輻射效率很高。

圖5 鋼軌聲輻射效率Fig．5 Radiation efficiency of rail

圖6 軌道板聲輻射效率Fig．6 Radiation efficiency of track slab

計算整個時間內(nèi)軌道系統(tǒng)的輻射聲功率平均值，如圖7所示。并提取4個點場點整個時間內(nèi)的最大聲壓，如圖8所示。

圖7 軌道系統(tǒng)輻射聲功率Fig．7 Acoustic radiation power of rail system

圖8 點場點聲壓值Fig．8 The point field point sound pressure value

由圖7可知:分析鋼軌和軌道板對噪聲的貢獻(xiàn)量，發(fā)現(xiàn)軌道板的貢獻(xiàn)主要在500 Hz以下，鋼軌主要貢獻(xiàn)是在800 Hz以上，這與上面得到的結(jié)論是一樣的，即在800～3 000 Hz時主要是以鋼軌的聲輻射為主，在0～500 Hz時主要以軌道板的聲輻射為主。

由圖8可知，該圖反映了軌道系統(tǒng)聲輻射隨距離變化的情況，隨著距離的增長，軌道系統(tǒng)的聲輻射呈線性遞減趨勢，其中最大聲壓為102 dB。

圖9 研究線路圖Fig．9 Acoustic pressure nephogram at the start time

圖10 結(jié)束時刻處聲壓云圖Fig．10 Acoustic pressure nephogram at the end time

從圖9和10可知:軌道系統(tǒng)的聲輻射隨著時間的變化，有明顯的傳遞規(guī)律。當(dāng)力剛開始作用于軌道系統(tǒng)時，即鋼軌剛接受隨機力作用時，如圖9所示，聲輻射主要貢獻(xiàn)區(qū)就是軌道系統(tǒng)接受隨機力作用的區(qū)域;當(dāng)力作用于軌道系統(tǒng)的最后時刻，如圖10所示，情況與開始時刻的一樣。

另外，從軌道系統(tǒng)表面聲壓云圖來看，所計算軌道系統(tǒng)，其鋼軌部位聲輻射要比軌道板部位的顯著，平均大15 dB左右。并且觀察鋼軌的表面聲壓云圖，發(fā)現(xiàn)軌腰的聲輻射量要比其他地方的顯著。所以，若能有效降低鋼軌軌腰部位的聲輻射，即可有效降低鋼軌的聲輻射，同時也能降低軌道系統(tǒng)的總噪聲量［18］。這也是現(xiàn)有許多低噪聲鋼軌設(shè)計從控制鋼軌軌腰振動、聲輻射著手的理論依據(jù)。

4 結(jié)論

(1)運用直接邊界元法計算軌道系統(tǒng)的聲輻射特性，得到了鋼軌與軌道板的聲輻射效率、輻射聲功率和軌道系統(tǒng)的聲場輻射云圖。分析鋼軌和軌道板對噪聲的貢獻(xiàn)量，發(fā)現(xiàn)軌道板的貢獻(xiàn)主要在500 Hz以下，鋼軌主要貢獻(xiàn)是在800 Hz以上，即在800～3 000 Hz時主要是以鋼軌的聲輻射為主，在0～500 Hz時主要以軌道板的聲輻射為主;隨著距離的增長，軌道系統(tǒng)的聲輻射呈線性遞減趨勢;鋼軌部位聲輻射要比軌道板部位的顯著，平均大15 dB左右，軌腰的聲輻射量要比鋼軌其他部位的顯著;

(2)首先采用有限元法求解非線性的多剛體耦合動力學(xué)，接著運用直接邊界元理論計算輪軌的聲輻射情況，這樣計算出的輪軌滾動噪聲更加接近現(xiàn)實情況;

(3)計算出軌道系統(tǒng)輻射噪聲，不管與既有的且己被專家公認(rèn)比較合理的模型計算結(jié)果相比，還是與實測的結(jié)果相比，本文軌道系統(tǒng)噪聲預(yù)測的結(jié)果與這些結(jié)果都有很好的一致性，說明了本文模型與做法是合理可行的。

［1］雷曉燕，圣小珍．鐵路交通噪聲與振動［M］．北京:科學(xué)出版社，2004．LEI Xiao－yan，SHENG Xiao－zhen．Rail noise and vibration［M］．Beijing:Science Press，2004．

［2］Kitagawa T．The influence of wheel and track parameters on rolling noise［J］．QR of RTRI，2009，50:32 －38．

［3］Thompson D J，Jones C J C．A review of the modelling of wheel/rail noise generation ［J］．Journal of Sound and Vibration，2000，231(3):519 －536．

［4］Remington P J．Wheel/rail noise，Partl:ehara eterization of the wheel/rail dynamie system［J］．Journal of Sound and Vibration，1976(46):359 －379．

［5］Thompson D J，Jones C JC，Turner N．Investigation into the validity of two－ dimensional models for sound radiation from waves in rails［J］．Acoustical Society of America，2003，113(41):1965 －1974．

［6］Thompson D J，Jones C J C．Sound radiation from avibrating railway wheel［J］．Journal of Sound and Vibration，2002，253(2):401 －419．

［7］Nielsen J C O．Acoustic optimization of railway sleepers［J］．Journal of Sound and Vibration，2000，231(3):753－764．

［8］徐志勝，翟婉明．軌道交通輪軌噪聲預(yù)測模型［J］．交通運輸工程學(xué)報，2005，5(3):14 －19．XU Zhi－sheng，ZHAI Wan－ming．Prediction model of wheel－ rail noise for rail transit［J］．Journal of Traffic and Transportation Engineering，2005，5(3):14 －19．

［9］萬淑敏，吳天行．鐵路鋼軌垂向振動的聲輻射分析［J］．噪聲與振動控制，2009(3):86 －89．WAN Shu－min，WU Tian－xing．Analysis of rail acoustic radiation due to vertical vibration［J］．Noise and Vibration Control，2009(3):86 －89．

［10］方 銳，肖新標(biāo)．軌道結(jié)構(gòu)參數(shù)對鋼軌和軌枕振動特性的影響［J］．鐵道學(xué)報，2011，33(3):71 －76．FANG Rui，XIAO Xin－biao．Effect of structure parameters of railway track on dynamic behavior of rail and sleeper［J］．Journal of the China Railway Society，2011，33(3):71－76．

［11］魏 偉，翟婉明．輪軌系統(tǒng)高頻振動響應(yīng)［J］．鐵道學(xué)報，1999，20(2):33 －36．WEIWei，ZHAIWan－ming． Dynamicresponseof wheel/rail system to high frequency excitation［J］．Journal of the China Railway Society，1999，20(2):33 －36．

［12］劉林芽，呂 銳，劉海龍．無砟軌道垂向高頻振動響應(yīng)分析［J］．鐵道科學(xué)與工程學(xué)報，2011，8(6):1 －6．LIU Lin－ya，LV Rui，LIU Hai－long．Vertical high frequency vibration response analysis of ballastless track［J］．Journal of Railway Science and Engineering，2011，8(6):1－6．

［13］劉林芽，劉海龍，雷曉燕．輪軌系統(tǒng)車輪的高頻振動特性［J］．交通運輸工程學(xué)報，2011，11(6):44 －49．LIU Lin－ya，LIU Hai－long，LEI Xiao－yan．High － frequency vibration analysis of wheels with wheel/rail system［J］．Journal of Traffic and Transportation Engineering，2011，11(6):44 －49．

［14］楊德慶，王德禹，劉洪林．艦艇振動聲學(xué)特性數(shù)值分析［J］．上海交通大學(xué)學(xué)報，2002，36(11):87 －93．YANG De－qing，WANG De－Yu，LIU Hong－lin．Numerical analysis of vibro － acoustic characters of ship［J］．Journal of Shanghai Jiaotong University，2002，36(11):87－93．

［15］黎 勝，趙德有．半空間內(nèi)結(jié)構(gòu)聲輻射研究［J］．船舶力學(xué)，2004，8(2):45 －51．LI Sheng，ZHAO De－you．Research on acoustic radiation in a three － dimensional half space［J］．Journal of Ship Mechanics，2004，8(2):45 －51．

［16］Moyne S L，Tebec J L．Ribs effects in acoustic radiation of a gear box－their modelling in a boundary element method［J］．Applied Acoustics，2002，63(5):223 －233．

［17］房建英，肖新標(biāo)，金學(xué)松．行車速度對高速列車車輪振動聲輻射特性的影響［J］．機械工程學(xué)報，2010，46(22):97－104．FANG Jian－ying，XIAO Xin－biao，JIN Xue－song．Effect of train speed on acoustic radiation characteristics of high－ speed train wheel vibration［J］．Journal of Mechanical Engineering，2010，46(22):97 －104．

［18］李再幃，雷曉燕，張鵬飛．減振型阻尼鋼軌有限元分析［J］．噪聲與振動控制，2009，29(4):64 －66．LI Zai－wei，LEI Xiao－yan，ZHANG Peng－fei．FEM analysis of damped rails for vibration reduction［J］．Noise and Vibration Control，2009，29(4):64 －66．