張鵬飛，雷曉燕，高亮，劉慶杰

(1.北京交通大學(xué) 土木建筑工程學(xué)院，北京 100044;2.華東交通大學(xué)鐵路環(huán)境振動與噪聲教育工程研究中心，江西 南昌 330013)

近年來，鐵路和城市軌道交通在我國得到了快速發(fā)展，軌道交通引起的環(huán)境振動問題引起了政府和民眾的普遍關(guān)注。鐵路環(huán)境振動會對居民的日常生活和工作、古建筑的結(jié)構(gòu)安全以及工廠、學(xué)校內(nèi)精密儀器設(shè)備的生產(chǎn)和使用產(chǎn)生很大的影響［1］。為了解決軌道交通引起的環(huán)境振動問題，國內(nèi)外學(xué)者進(jìn)行了大量的研究工作。Sheng等［2－6］采用解析的波數(shù)－頻率域法建立了軌道結(jié)構(gòu)單層或多層梁模型，分析了高速列車引起的軌道和大地振動。國外學(xué)者［7－8］對精密儀器的環(huán)境振動標(biāo)準(zhǔn)進(jìn)行了詳細(xì)討論，Xia等［9－11］對地鐵列車運行引起的振動對精密儀器的影響作了大量的現(xiàn)場測試和理論分析工作，這些研究取得了很多成果，得出了大量有意義的結(jié)論。本文以新建龍煙鐵路沿線附近廠房內(nèi)的精密儀器為研究對象，利用有限元法，建立軌道—橋梁及橋墩－基礎(chǔ)－大地的計算模型，利用這2個模型計算貨車在不同運行速度下引起的大地振動響應(yīng)，得出振動加速度、速度及位移沿距離的衰減規(guī)律，分析鐵路環(huán)境振動對沿線廠房內(nèi)精密儀器的影響。

1 背景資料

1.1 龍煙鐵路概況

龍煙鐵路西起大萊龍鐵路的龍口西站(不含)，經(jīng)龍口市、蓬萊市、經(jīng)濟(jì)開發(fā)區(qū)、福山區(qū)、芝罘區(qū)至珠璣站，新建正線全長112.735 km;福山聯(lián)絡(luò)線13.166 km。與規(guī)劃德大鐵路、既有大萊龍鐵路共同形成太原煤炭基地至環(huán)渤海港口煙臺(龍口)港的便捷通道，是環(huán)渤海鐵路、沿海鐵路的重要組成部分。線路為Ⅰ級單線(預(yù)留雙線條件)鐵路，旅客列車采用城際動車組，最高運行速度為160 km/h，貨物列車采用SS4和HXD1型電力機車牽引，最高運行速度為120 km/h，全線采用有砟軌道。

1.2 精密儀器

在煙臺開發(fā)區(qū)，龍煙鐵路以高架形式近距離經(jīng)過煙臺勝利汽車零部件制造有限公司，該公司內(nèi)設(shè)有兩臺精密儀器－光譜儀和制芯機，根據(jù)線路設(shè)計平面圖和現(xiàn)場測量，光譜儀到線路中心線的垂直距離為96.93 m，振動參數(shù)限值要求如下:加速度＜1.98 m/s2;速度 ＜8 ×10－3m/s;位移 ＜8 ×10－5m;制芯機到線路中心線的垂直距離為124.75 m，振動參數(shù)限值要求如下:x和y方向的加速度 ＜1.15 m/s2，z方向的加速度 ＜0.8 m/s2。

2 有限元模型

為模擬列車運行產(chǎn)生的振動，建立軌道－橋梁耦合有限元模型，將車輛(機車)簡化為一系列以速度v移動的荷載，計算作用于橋梁支座上的反力，建立橋墩－基礎(chǔ)－大地三維有限元模型，將第1個模型計算的力加載到第2個模型上，計算與線路中心線不同距離的地面振動加速度、速度和位移。

2.1 有限元振動計算理論

利用有限元方法進(jìn)行結(jié)構(gòu)振動計算，可以得到在移動荷載作用下結(jié)構(gòu)隨時間變化的節(jié)點加速度、速度和位移的響應(yīng)，基本方程為:

式中:［Me］，［Ce］和［Ke］分別為結(jié)構(gòu)質(zhì)量矩陣、阻尼矩陣和剛度矩陣;{}，{}和{u}分別為e結(jié)構(gòu)振動加速度、速度和位移;{Fe}為荷載向量。對于上述微分方程的求解可采用Newmark時間積分方法在離散的時間點上求解動力學(xué)方程，系統(tǒng)采用Raleigh阻尼。

2.2 軌道－橋梁有限元模型

該模型從上到下依次為:60 kg/m鋼軌、扣件、軌枕、道砟、橋梁。

鋼軌采用梁單元Beam3來模擬，軌枕用質(zhì)量點單元Mass21等效，墊板、扣件以及道床可以采用彈簧阻尼單元combin14等效。軌道－橋梁有限元力學(xué)模型如圖1所示。

圖1 軌道－橋梁有限元力學(xué)模型Fig.1 Track－bridge finite element mechanics model

根據(jù)一般軌道的實際支撐情況取鋼軌梁單元的長度為0.6m，鋼軌的每個節(jié)點處用彈簧阻尼單元連接一個軌枕質(zhì)量點單元，再由彈簧阻尼單元連接到橋梁節(jié)點上。模型的參數(shù)取值見表1。

2.3 橋墩—基礎(chǔ)—大地有限元模型

橋墩—基礎(chǔ)—大地模型考慮5跨簡支梁橋，橋墩尺寸由線路設(shè)計單位提供。

橋墩及承臺采用solid45單元，最小單元尺寸0.5 m;樁基礎(chǔ)采用beam4模擬，樁頂單元的節(jié)點與承臺節(jié)點耦合，為方便計算，樁長均設(shè)為42 m，大地邊界采用粘彈性吸收邊界，用以吸收邊界上的反射波［12］。橋墩－樁基礎(chǔ)模型如圖2所示，橋墩－基礎(chǔ)－大地有限元模型如圖3所示，土層參數(shù)由線路勘探設(shè)計單位提供，見表2。

表1 軌道－橋梁模型參數(shù)表Table 1 Track－bridge model parameters

圖2 橋墩－樁基礎(chǔ)有限元模型Fig.2 Pier－ pile foundation finite element model

圖3 橋墩－基礎(chǔ)－大地有限元模型Fig.3 Pier－foundation－ground finite element model

表2 土層參數(shù)表Table 2 Soil parameters

3 計算工況

為了模擬振動最不利的情況，計算時考慮兩列貨車在煙臺勝利汽車零部件制造有限公司處交會，貨車速度分別取120，100及80 km/h。各計算工況見表3。

表3 計算工況Table 3 Calculation modes

4 計算結(jié)果分析

4.1 橋墩支座反力計算

在高架軌道誘發(fā)環(huán)境振動的研究中，車輛(機車)的振動并不是主要考慮因素，因此將其簡化為一系列以速度v移動的荷載。在考慮軌道不平順對輪軌力的影響時，可以將不平順導(dǎo)致的慣性力疊加在輪對靜荷載上，軌道不平順采用美國6級不平順譜樣本，如圖4所示。整個分析計算過程采用全瞬態(tài)的分析方法，時間積分步長取為0.005 s，計算的總時間7.8 s。為取得有效的計算精度，整列車的車輛數(shù)取為8節(jié)。

圖4 不平順樣本Fig.4 Irregularities samples

利用軌道—橋梁模型，帶入表1中的參數(shù)，可以得到3種工況下的橋墩支座反力，如圖5～7所示。

圖5 工況1橋墩支座反力Fig.5 Piers supporting force of Mode 1

圖6 工況2橋墩支座反力Fig.6 Piers supporting force of Mode 2

4.2 振動傳播規(guī)律分析

圖7 工況3橋墩支座反力Fig.7 Piers supporting force of Mode 3

圖8 加速度沿距離的變化曲線Fig.8 Acceleration curve along the distance

圖9 速度沿距離的變化曲線Fig.9 Speed curve along the distance

將橋墩支座反力計算結(jié)果加載到橋墩－基礎(chǔ)－大地有限元模型中，采用瞬態(tài)分析方法對距離線路中心線不同位置的地面振動進(jìn)行分析，圖8～10所示為各工況振動加速度、速度及位移沿距離的變化規(guī)律曲線。從圖8～10可以看出:貨物列車速度對振動響應(yīng)的影響較大，速度越高，地面上相同位置上的振動響應(yīng)越大;各工況振動指標(biāo)沿距離衰減特性基本一致，呈現(xiàn)指數(shù)曲線的衰減規(guī)律，在距離線路中心線80 m左右的范圍內(nèi)，振動迅速衰減，但在20－40 m范圍內(nèi)有振動“反彈增大”現(xiàn)象出現(xiàn)，超過80 m之后，振動表現(xiàn)為緩慢衰減的趨勢。

4.3 振動對精密儀器的影響分析

光譜儀與線路中心線之間的距離為96.93 m，分別提取光譜儀位置處各工況下的振動加速度、速度和位移，見表4。制芯機與線路中心線之間的距離為124.75 m，分別提取制芯機位置處各工況下的振動加速度、速度和位移，見表5。

表4 光譜儀的振動參數(shù)值Table 4 Spectrometer vibration parameters

表5 制芯機的振動參數(shù)值Table 5 Core－making machine vibration parameters

從表4和表5可以看出:當(dāng)2貨車以120，100和80 km/h的速度在精密儀器附近交會時，光譜儀的加速度、速度及位移的預(yù)測值滿足儀器振動限值要求;制芯機的加速度預(yù)測值滿足儀器的限值要求。

5 結(jié)論

(1)貨物列車速度對振動響應(yīng)的影響較大，速度越高，地面上相同位置上的振動響應(yīng)越大。

(2)振動加速度、速度和位移沿距離的衰減特性基本一致，呈現(xiàn)指數(shù)曲線的衰減規(guī)律，在距離線路中心線80 m左右的范圍內(nèi)，振動迅速衰減，但在20～40 m范圍內(nèi)有振動“反彈增大”現(xiàn)象出現(xiàn)，超過80 m之后，振動表現(xiàn)為緩慢衰減的趨勢。

(3)當(dāng)兩貨車以120，100和80 km/h的速度在精密儀器附近交會時，光譜儀的加速度、速度及位移的預(yù)測值滿足儀器振動限值要求;制芯機的加速度預(yù)測值滿足儀器的限值要求。

［1］Crispino M，Dapuzzo M.Measurement and prediction of traffic－ induced vibrations in a heritage building［J］.Journal of Sound and Vibration，2001，246(2):319 －335.

［2］Sheng X，Jones C J C，Petyt M.Ground vibration generated by a load moving along a railway track［J］.Journal of Sound and Vibration，1999，228(1):129 －156.

［3］Picoux B，Houedec D L.Diagnosis and prediction of vibration from railway trains［J］.Soil Dynamics and Earthquake Engineer－ing，2005(25):905 －921.

［4］謝偉平，胡建武，徐 勁.高速移動荷載作用下的軌道－地基系統(tǒng)的動力響應(yīng)［J］.巖石力學(xué)與工程學(xué)報，2002，21(7):1075 －1078.XIE Weiping，HU Jianwu，XU Jin.Dynamic response of track－ground systems under high speed moving lord［J］.Chinese Journal of Rock Mechanics and Engineering，2002，21(7):1075 －1078.

［5］邊學(xué)成，陳云敏.列車荷載作用下軌道和地基的動響應(yīng)分析［J］.力學(xué)學(xué)報，2005，37(4):477 －484.BIAN Xuecheng，CHEN Yunmin.Dynamic analyses of track and ground coupled system with high－speed train loads［J］.Chinese Journal of Theoretical and Applied Mechanics，2005，37(4):477 －484.

［6］雷曉燕.高速列車誘發(fā)地面波與軌道強振動研究［J］.鐵道學(xué)報，2006，28(3):78 －82.LEI Xiaoyan.Study on ground waves and track vibration boom induced by high speed trains［J］.Journal of the China Railway Society，2006，28(3):78 －82.

［7］Amick H.On generic vibration criteria for advanced technology facilities［J］.Journal of the Institute of Environmental Sciences，1997，40(5):35 －44.

［8］Gordon C G.Generic vibration criteria for vibration sensitive equipment［C］.Proceedings of the 1999 Optomechanical Engineering and Vibration Control.Denver，Colorado，USA:Society of Photo－Optical Instrumentation Engineers，1999:22 －33.

［9］Xia He.Vibration effects of light－rail train－viaduct system on surrounding environment［J］.Int J of Structural Stability and Dynamics.2002，2(2):227－ 240.

［10］馬 蒙，劉維寧，丁德云，等.地鐵列車振動對精密儀器影響的預(yù)測研究［J］.振動與沖擊，2011，30(3):185－190.MA Meng，LIU Weining，DING Deyun，et al.Prediction of influence of metro trains induced vibrations on sensitive instruments［J］.Journal of Vibration and Shock，2011，30(3):185－190.

［11］栗潤德，張鴻儒，劉維寧.地鐵引起的地面振動及其對精密儀器的影響［J］.巖石力學(xué)與工程學(xué)報，2008，27(1):206－213.LI Runde，ZHANG Hongru，LIU Weining.Metro－induced ground vibrations and their impacts on precision instrument［J］.Chinese Journal of Rock Mechanics and Engineering，2008，27(1):206－213.

［12］雷曉燕.城市軌道交通環(huán)境振動預(yù)測與評價［J］.城市軌道交通研究，2008(11):12－17.LEI Xiaoyan.Prediction and assessment of environment vibration induced by urban mass transit［J］.Urban Mass Transit，2008(11):12 －17.