鄧 躒，施 洲，勾紅葉

(1.中鐵二院工程集團(tuán)有限責(zé)任公司 環(huán)境工程研究院，成都 610031;2.西南交通大學(xué) 土木工程學(xué)院，成都 610031)

隨著高速鐵路的飛速發(fā)展，在帶來方便快捷的運(yùn)輸條件時(shí)，也給鐵路沿線的環(huán)境造成一定影響，其中噪聲問題將日益突出。鑒于目前高速鐵路最高運(yùn)營速度為350 km/h，因此350 km/h以上高速鐵路相關(guān)列車脈動(dòng)風(fēng)荷載的相關(guān)研究資料相對(duì)較少，更無試驗(yàn)實(shí)測資料。時(shí)速350 km及以下高速鐵路列車脈動(dòng)風(fēng)荷載的既有研究資料已有不少，其中最為系統(tǒng)的屬于德國的相關(guān)資料。相應(yīng)的資料為研究380 km/h高速鐵路列車脈動(dòng)風(fēng)荷載的特性提供良好的基礎(chǔ)。因此，在既有350 km/h及以下速度相關(guān)資料的基礎(chǔ)上，通過理論分析，CFD仿真分析的方法，系統(tǒng)研究380 km/h高速鐵路列車脈動(dòng)風(fēng)荷載的特性。

1 既有高速鐵路列車脈動(dòng)風(fēng)荷載研究資料

1.1 ICE3列車脈動(dòng)風(fēng)荷載研究成果

高速列車的脈動(dòng)風(fēng)荷載取決于以下參數(shù):列車的類型(特別是空氣動(dòng)力特性和車廂長度)、列車的速度、軌道中心線和聲屏障之間的距離、聲屏障的形狀和高度。德國的研究資料表明，由16車編組ICE3列車經(jīng)過產(chǎn)生的脈動(dòng)風(fēng)荷載見圖1所示。最大壓力和吸力出現(xiàn)在列車前面。車廂的每個(gè)接合處引起較小的壓力和吸力載荷。16車編組ICE3的掛鉤引起較大的風(fēng)力循環(huán)載荷。在列車車尾經(jīng)過時(shí)，出現(xiàn)下一次高吸力和壓力載荷。8車編組列車經(jīng)過時(shí)，會(huì)出現(xiàn)類似的脈動(dòng)風(fēng)載荷。德國的既有研究是基于滿載列車的試驗(yàn)測試結(jié)果，測試并記錄了在德國高速鐵路線上高速列車ICE3運(yùn)行速度300 km/h下的風(fēng)力。屏障至軌道中心線的距離為3.8 m，軌面上屏障高度為3.0 m。測量數(shù)據(jù)表明，最大脈動(dòng)風(fēng)壓力為450 Pa，最小值為－343 Pa。

圖1 16車編組ICE3列車經(jīng)過產(chǎn)生的脈動(dòng)風(fēng)荷載示意

在德國(圖1)，已在現(xiàn)有高速鐵路聲屏障脈動(dòng)風(fēng)壓測量結(jié)果的基礎(chǔ)上研發(fā)了脈動(dòng)風(fēng)荷載的公式。對(duì)于ICE3列車脈動(dòng)風(fēng)壓力和吸力，脈動(dòng)風(fēng)荷載的定義如下

式中，wk為脈沖風(fēng)力荷載(Pa);vTrain為列車的速度(m/s);ρAir為空氣密度(kg/m3);cz為軌面以上的高度系數(shù)，在聲屏障的上端為0.6，并在上部2/3高度處線性變化至1.0，下部1/3高度為1.0;cp為脈動(dòng)風(fēng)壓力系數(shù)，取決于列車類型和到軌道中心線的距離。

式中，ag為聲屏障和軌道中心線的距離;cp(3.8)即為聲屏障和軌道中心線之間3.8 m距離的脈動(dòng)風(fēng)壓力系數(shù)。該系數(shù)取決于時(shí)間并基于測量結(jié)果確定。對(duì)于其它的速度vk，時(shí)間比值將按以下公式轉(zhuǎn)換

1.2 中國高速鐵路列車脈動(dòng)風(fēng)荷載

脈動(dòng)風(fēng)荷載主要取決于列車的類型、速度、空氣密度、屏障的高度和形狀、屏障和軌道中心線之間的距離。因此ICE3的脈動(dòng)風(fēng)荷載定義原則上也可適用于其它類型的列車。其它類型列車的列車系數(shù)只有系數(shù)cp(ag)可以通過測量決定。中國的高速列車CRH系列列車是基于不同的國際高速列車系統(tǒng)發(fā)展而成。參考德國的 DB cp曲線，可以推導(dǎo)得出高速鐵路列車CRH2和CRH3的插值 cp曲線。與標(biāo)準(zhǔn) ICE3列車編制的德國DB cp曲線相比，考慮到CRH列車的不同長度和寬度的實(shí)際情況。德國研究資料對(duì)CHR3列車cp曲線的建議:CHR3列車與ICE3列車在一個(gè)截面上十分相似。在縱向方向，幾何形狀相同。在橫向方向，其寬度更寬。相對(duì)于ICE3列車，CRH3聲屏障與線路中心的距離ag相對(duì)變短，因而系數(shù)cp(ag)大約大7.3%，由于CHR3的尺寸較大，排換的空氣量較大。與一列ICE3列車的cp曲線相比，兩種因素引起cp(3.8)系數(shù)增加了10%。德國研究資料對(duì)CHR2列車cp曲線的建議:CHR2列車的寬度同樣寬于ICE3列車，因此，系數(shù)cp(ag)大約大10.2%，由于CHR2的尺寸較大，排換的空氣量比ICE3較大。CHR2列車的空氣動(dòng)力形狀比CHR3列車和ICE3列車稍差。CHR2列車(201.4 m)的長度僅僅比CHR3列車(206.8 m)長0.4%。這個(gè)因素可以忽略不計(jì)。CHR2列車cp(3.8)與一列ICE3列車的cp曲線相比，所有因素導(dǎo)致cp(3.8)系數(shù)增加了16%。

盡管德國研究資料給出了良好的建議，在實(shí)際應(yīng)用中，應(yīng)在試驗(yàn)段采用實(shí)際列車運(yùn)行來進(jìn)行現(xiàn)場測量以便檢驗(yàn)該插值曲線。對(duì)于最終的聲屏障動(dòng)力設(shè)計(jì)，只有進(jìn)行真正的實(shí)際測量才能保證脈動(dòng)風(fēng)荷載的準(zhǔn)確性。

2 高速列車作用聲屏障風(fēng)荷載的CFD理論與模型

高速列車通行聲屏障區(qū)域時(shí)，在列車與聲屏障之間產(chǎn)生氣動(dòng)力作用，可通過計(jì)算流體力學(xué)的方法來模擬高速列車運(yùn)行導(dǎo)致的脈動(dòng)風(fēng)荷載作用，即 CFD計(jì)算。在CFD數(shù)值計(jì)算中，高速列車從接近并駛?cè)肼暺琳虾笏鸬穆暺琳蟽?nèi)空氣流動(dòng)是復(fù)雜的三維、非定常、可壓縮的紊流流動(dòng)?？刹捎么笮土黧w數(shù)值模擬軟件FLUENT，可針對(duì)列車行駛通過聲屏障前后時(shí)，準(zhǔn)確確定聲屏障上受到的風(fēng)壓。

2.1 數(shù)值計(jì)算方程

為了保證計(jì)算精度，在高速列車及聲屏障的脈動(dòng)力計(jì)算中，采用黏性流體計(jì)算，并將空氣視為理想氣體，求解納維—斯托克斯(Navier-Stokes)方程(以下簡稱為N-S方程)，得到流場的數(shù)值解。

N-S方程包括質(zhì)量守恒方程、動(dòng)量守恒方程和能量守恒方程。

2.2 CFD計(jì)算理論與計(jì)算模型

FLUENT提供了分離式和耦合式兩類求解器，分離式是順序的、逐一地求解各方程(關(guān)于 u，v，w，p，T的方程);耦合式是同時(shí)求解連續(xù)方程、動(dòng)量方程、能量方程的耦合方程，然后再逐一地求解湍流等標(biāo)量方程。耦合式求解器，可采用隱式和顯式兩種方案進(jìn)行控制方程的線性化。顯式對(duì)于給定變量，每一個(gè)單元內(nèi)的未知量用只包含已知值的關(guān)系式來計(jì)算;隱式對(duì)于給定變量，單元內(nèi)的未知量用鄰近單元的已知和未知值來計(jì)算。

在高速列車及聲屏障的脈動(dòng)力計(jì)算中空氣介質(zhì)取黏性、可壓縮理想氣體，紊流方程采用 k-ε模型，采用耦合式求解器隱式方案對(duì)三維Navier-Stocks方程求解。采用動(dòng)網(wǎng)格法，模擬列車進(jìn)入聲屏障的過程。

在高速列車作用于聲屏障的CFD分析中，線路考慮為雙線，線路間距為5.0 m，聲屏障距離線路中心3.8 m。分析中考慮聲屏障距離為400 m，列車長為200 m。列車通過聲屏障，測點(diǎn)布置圖形和計(jì)算長度圖形如圖2所示。CFD仿真分析模型見圖3。

圖2 計(jì)算模型布置示意(單位:m)

圖3 列車進(jìn)入聲屏障區(qū)域示意

3 高速列車作用于聲屏障風(fēng)荷載的CFD分析結(jié)果

3.1 豎直聲屏障CFD分析結(jié)果

當(dāng)列車以350 km/h，380 km/h的速度通行聲屏障區(qū)域，通過CFD分析聲屏障承受動(dòng)態(tài)脈動(dòng)風(fēng)壓力作用。計(jì)算結(jié)果中，主要提取兩側(cè)聲屏障自列車入口處50，100，200，300，400 m 處的計(jì)算結(jié)果?？拷熊囘\(yùn)行一側(cè)的聲屏障，即近側(cè)聲屏障區(qū)間內(nèi)壓力波傳播情況分別如圖4、圖5所示，不同車速下近側(cè)聲屏障上不同測點(diǎn)壓力峰值如表1。另外一側(cè)，即遠(yuǎn)側(cè)聲屏障承受的脈動(dòng)壓力相對(duì)較小，近側(cè)、遠(yuǎn)側(cè)聲屏障的風(fēng)壓力結(jié)果見表2。

圖4 列車350 km/h通行50 m處測點(diǎn)的壓力時(shí)程

圖5 列車380 km/h通行50 m處測點(diǎn)的壓力時(shí)程

表1 聲屏障上不同測點(diǎn)壓力峰值計(jì)算結(jié)果Pa

表2 380 km/h下近遠(yuǎn)側(cè)聲屏障上不同測點(diǎn)壓力峰值計(jì)算結(jié)果Pa

從高速列車以350 km/h，380 km/h的速度通行聲屏障區(qū)域的CFD計(jì)算分析結(jié)果可見，作用于聲屏障的最大風(fēng)壓力出現(xiàn)在聲屏障的下部，且最大正壓力稍大于最大負(fù)壓力。350 km/h速度下最大風(fēng)壓力為1 474 Pa，380 km/h速度下最大風(fēng)壓力為1 707 Pa。在單一列車通行時(shí)，線路兩側(cè)聲屏障承受的列車脈動(dòng)風(fēng)壓相差顯著，遠(yuǎn)側(cè)的聲屏障承受風(fēng)壓力遠(yuǎn)小于近側(cè)聲屏障風(fēng)荷載。

3.2 聲屏障CFD分析與討論

高速列車通行聲屏障區(qū)域的CFD結(jié)果表明，在聲屏障的高度方向上，風(fēng)荷載分布有一定的差異，350 km/h列車通行時(shí)，沿聲屏障高度方向的最大風(fēng)壓力分布如圖6所示，380 km/h列車通行時(shí)，沿聲屏障高度方向的最大風(fēng)壓力分布見圖7。

圖6 350 km/h下聲屏障測點(diǎn)正壓力豎向分布

圖7 380 km/h下聲屏障測點(diǎn)正壓力豎向分布

從圖中可見，聲屏障底部承受的風(fēng)荷載最大，并沿高度向上至聲屏障一半高度 H/2處緩慢變小，并在H/2至頂部附近較快減小至底部風(fēng)荷載的1/2左右，并在頂部大約30 cm高度范圍內(nèi)保持不變。

在列車通行聲屏障區(qū)域時(shí)，沿聲屏障縱向的最不利風(fēng)荷載也有一定的差異。列車以380 km/h通行過程中，聲屏障承受最大正壓力沿縱向分布見表3、圖8。分析結(jié)果可見，在聲屏障的頂部，縱向最不利風(fēng)壓結(jié)果相對(duì)較小且沿縱向基本不變，而聲屏障的底部及中部沿縱向分布呈現(xiàn)在列車入口處最小，沿著列車前進(jìn)方向50 m處迅速增大，后稍減小并在100～400 m處即聲屏障末尾保持較平穩(wěn)值。

3.3 改進(jìn)聲屏障的CFD分析結(jié)果

為進(jìn)一步分析不同聲屏障結(jié)構(gòu)形式對(duì)列車脈動(dòng)風(fēng)壓力結(jié)果的影響，分析了兩種改進(jìn)的聲屏障結(jié)構(gòu)在380 km/h下的CFD結(jié)果。一種是在聲屏障的開始端與末尾端沿縱向在聲屏障的高度采用三角形過渡，另外一種是在頂端設(shè)置向線路內(nèi)側(cè)傾斜的弧形頂部裝置的聲屏障。改進(jìn)后的聲屏障最不利脈動(dòng)風(fēng)壓力結(jié)果見表4。其中弧形頂部的風(fēng)壓脈動(dòng)時(shí)程見圖9。

表3 聲屏障上測點(diǎn)壓力縱向比較(工況Ⅰ，V=380 km/h) Pa

圖8 聲屏障上測點(diǎn)壓力縱向比較(工況Ⅰ，V=380 km/h)

表4 近端遠(yuǎn)端聲屏障上不同測點(diǎn)壓力峰值計(jì)算結(jié)果Pa

從改進(jìn)的聲屏障CFD分析結(jié)果可見，改進(jìn)的三角形聲屏障可一定程度降低列車通行時(shí)的脈動(dòng)風(fēng)壓力值。增加弧形頂部后，聲屏障的脈動(dòng)風(fēng)壓力變得更為復(fù)雜一些，在縱向及聲屏障的高度上呈現(xiàn)一些并不規(guī)律的變化，但其量值并未顯著增加。

4 結(jié)論

詳細(xì)分析高速鐵路列車脈動(dòng)風(fēng)荷載的特性，進(jìn)行高速列車以350 km/h，380 km/h的速度通行聲屏障區(qū)域的CFD計(jì)算分析。結(jié)果表明，作用于聲屏障的最大列車脈動(dòng)風(fēng)壓力出現(xiàn)在聲屏障的下部，且最大正壓力稍大于最大負(fù)壓力。350 km/h速度下最大風(fēng)壓力為1 474 Pa，380 km/h速度下最大風(fēng)壓力為1 707 Pa。在單一列車通行時(shí)，線路遠(yuǎn)側(cè)的聲屏障承受風(fēng)壓力遠(yuǎn)小于近側(cè)聲屏障風(fēng)荷載。

聲屏障底部承受的風(fēng)荷載最大，并沿高度向上至聲屏障一半高度H/2處緩慢變小，并在H/2至頂部附近較快減小至底部風(fēng)荷載的1/2左右，并在頂部大約30 cm高度范圍內(nèi)保持不變。

沿聲屏障縱向的最不利風(fēng)荷載也有一定的差異，在聲屏障的頂部，縱向最不利風(fēng)壓結(jié)果相對(duì)較小且沿縱向基本不變，而聲屏障的底部及中部沿縱向分布呈現(xiàn)在列車入口處最小，沿著列車前進(jìn)方向50 m處迅速增大，后稍減小并在100～400 m處即聲屏障尾保持較平穩(wěn)值。

端部三角形聲屏障可一定程度降低列車通行時(shí)的脈動(dòng)風(fēng)壓力值。增加弧形頂部后，聲屏障的脈動(dòng)風(fēng)壓力變得更為復(fù)雜一些，在縱向及聲屏障的高度上呈現(xiàn)一些并不規(guī)律的變化，但其量值并未顯著增加。

［1］德國PEC+S工程設(shè)計(jì)咨詢服務(wù)有限公司.客運(yùn)專線聲屏障咨詢報(bào)告［R］.北京:德國PEC+S工程設(shè)計(jì)咨詢服務(wù)有限公司，2007.

［2］馬筠.對(duì)高速鐵路聲屏障降噪效果影響因素的探討［J］.鐵道勞動(dòng)安全衛(wèi)生與環(huán)保，2008，35(1):5-8.

［3］焦長洲，高波，王廣地.聲屏障結(jié)構(gòu)的列車脈動(dòng)風(fēng)致振動(dòng)分析［J］.西南交通大學(xué)學(xué)報(bào)，2007，429(5):531-536.

［4］王德利.利用聲屏障降低鐵路噪聲對(duì)居民環(huán)境的影響［J］.沈陽師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，2007，25(3):381-383.

［5］李艷.聲屏障設(shè)計(jì)中的結(jié)構(gòu)驗(yàn)算［J］.噪聲與振動(dòng)控制，2007(3):97-99.

［6］趙麗濱，龍麗平，蔡慶云.列車風(fēng)致脈動(dòng)力下聲屏障的動(dòng)力學(xué)性能［J］.北京航空航天大學(xué)學(xué)報(bào)，2009，35(4):505-508.

［7］呂堅(jiān)品，張繼文，廖建州，等.既有鐵路橋梁聲屏障的高速列車脈動(dòng)風(fēng)致響應(yīng)［J］.西南交通大學(xué)學(xué)報(bào)，2009，44(4):547-581.

［8］董宇，成志強(qiáng)，朱正清.鐵路橋梁聲屏障氣動(dòng)力初步分析［J］.鐵道建筑，2010(7):47-49.