楊宏印, 陳志軍, 彭旭民
(1. 中鐵大橋 科學研究院有限公司,湖北 武漢 430034;2. 華中科技大學 土木工程與力學學院,湖北 武漢 430074)
?
基于子區(qū)間的車-線-橋振動不平順影響分析
楊宏印1, 陳志軍2, 彭旭民1
(1. 中鐵大橋 科學研究院有限公司,湖北 武漢 430034;2. 華中科技大學 土木工程與力學學院,湖北 武漢 430074)
提出了任意子波長區(qū)間的軌道不平順生成方法,討論了相應(yīng)時域樣本的特性。建立了耦合系統(tǒng)振動方程,對比分析了大波長區(qū)間及其子區(qū)間不平順的影響,研究了系統(tǒng)響應(yīng)的敏感波長。結(jié)果表明:軌道不平順樣本的功率譜和解析值一致,不平順生成方法是合理有效的;軌道不平順對橋梁位移影響很小,橋梁加速度的敏感波長區(qū)間為1~5 m; 0.05~1 m的短波不平順是軌道加速度的主要激勵源,而對車體加速度影響很小。線路養(yǎng)護和維修中應(yīng)加強對5 m以下波長不平順的控制;而為提高乘車舒適性,有必要對各種波長不平順進行控制。
軌道不平順;車-線-橋耦合振動;敏感波長;共振;動力響應(yīng)
軌道幾何不平順是引起輪軌動力作用變化進而誘發(fā)整個列車-軌道-橋梁系統(tǒng)耦合振動的主要激勵,是直接影響列車運行安全性和乘車舒適性的關(guān)鍵因素之一[1]。選用確定的不平順激擾函數(shù),文獻[2]分析了不平順波長和幅值對高速鐵路橋上列車走行性的影響。而軌道不平順是一個隨機過程[3],其功率譜密度函數(shù)能夠從波長、幅值來描述和揭示不平順的統(tǒng)計特征和規(guī)律。不同于確定性激擾函數(shù),隨機不平順采用連續(xù)波長區(qū)間來描述,已有研究[4]指出可通過功率譜密度函數(shù)得到其時域樣本,進而分析其對系統(tǒng)動力響應(yīng)的影響。盡管有關(guān)部門開展過不少研究工作,并結(jié)合實際工程給出了一些軌道譜表達式[5],但明確高速鐵路軌道不平順敏感波長仍是線路運營亟待解決的問題。
針對軌道隨機不平順對車-線-橋耦合振動的影響,本文提出了任意子波長區(qū)間的軌道隨機不平順生成方法,并討論了相應(yīng)時域樣本特性。然后結(jié)合我國高速鐵路線路實際,建立了耦合系統(tǒng)振動模型,分析了高速行車下各波長區(qū)間高低不平順對系統(tǒng)響應(yīng)的影響,進而對系統(tǒng)響應(yīng)的敏感波長進行了研究,得到了一些有益的結(jié)論。
軌道隨機不平順包含許多頻率和幅值成分,通常采用功率譜密度函數(shù)來描述[1]。考慮適用于德國時速250 km以上的高速鐵路低干擾譜[5],其高低不平順功率譜密度表達式為
(1)
式中,Sv(Ω)為功率譜密度;Ω為空間頻率(rad/m);Av為粗糙度常數(shù),其值為4.032×10-7m2·rad/m;Ωc、Ωr為截斷頻率,其值分別為0.824 6 rad/m和0.020 6 rad/m。
空間頻率Ω、時間頻率f、車速v及不平順波長λ滿足如下關(guān)系
根據(jù)2種譜密度在對應(yīng)譜帶寬度內(nèi)的均方值相等[6]和式(2),將式(1)變換為時間頻率表達式
式中,fr和fc均為時間截斷頻率。
根據(jù)離散傅里葉變換法[3],在式(3)功率譜密度函數(shù)上直接離散采樣,構(gòu)造出頻譜的幅值和相位,然后對其進行傅里葉逆變換(IFFT),從而得到時域的軌道不平順模擬樣本。對整個頻率范圍進行等帶寬離散,只需令不在截斷頻率范圍的譜密度為0,便可得到任意波長區(qū)間的時域樣本。
德國低干擾譜以中長波為主[7],而短波不平順所激發(fā)的高頻輪軌相互作用對系統(tǒng)響應(yīng)影響很大[8],且相應(yīng)的輪軌高頻噪聲是高速鐵路噪聲的重要來源[9],因而對于波長小于1 m的高低不平順采用軌道短波不平順來模擬,其功率譜密度函數(shù)為[10]
式中,S(fk)的單位為mm2/(1/m)。
根據(jù)子區(qū)間的思想,在相同的隨機相位下,分別生成波長區(qū)間為0.05~150 m及其所包含的子區(qū)間0.05~30 m、30~60 m、60~90 m、90~150 m共5個隨機不平順樣本,并對相關(guān)特性進行分析。波長區(qū)間為0.05~30 m和90~150 m的軌道不平順樣本分別見圖1和圖2。波長區(qū)間為0.05~30 m的不平順由2部分疊加而成,即0.05~1 m的短波長不平順和1~30 m的中波長不平順。對于生成的子波長區(qū)間不平順,波長越短,不平順變化越激烈,激勵頻率越大;波長越長,變化越平緩。
圖1 波長區(qū)間0.05~30 m不平順樣本
圖2 波長區(qū)間90~150 m不平順樣本
圖3 波長區(qū)間0.05~150 m不平順
圖3為大波長區(qū)間(0.05~150 m)的不平順樣本,也給出了所有子波長區(qū)間不平順疊加的結(jié)果。由圖3可見,大波長區(qū)間不平順和其子波長區(qū)間不平順疊加的結(jié)果完全一致,說明所提方法的模擬精度很高。分析發(fā)現(xiàn)生成的軌道不平順樣本的功率譜密度和解析值符合很好,表明各隨機不平順樣本均是合理有效的。
2.1 車-線-橋耦合振動方程
圖4為列車過橋時豎向耦合系統(tǒng)模型,橋上鋪設(shè)雙塊式無砟軌道,支承間距為0.65 m。列車采用CRH2動車組,由8輛四軸10自由度車輛編組而成。采用Euler-Bernoulli梁模擬軌道和橋梁,每個單元有4個自由度。將軌道和橋梁離散,通過Hertz彈簧來模擬輪軌接觸,根據(jù)文獻[11-12]方法,可得到車-線-橋耦合系統(tǒng)動力方程
式(5)中,下標bb、rr和vv分別代表橋梁、軌道和列車;Kcr、Krv、Kvr、Fcr和Fcv為與輪軌接觸有關(guān)的項,是時變的,而其它項均是時不變的。故在求解時可先令時變部分為0,組裝時不變部分,作為初始方程,然后在每個積分步疊加時變部分便可得到系統(tǒng)動力方程,再采用Newmark-β方法直接積分求解[13]。采用MATLAB語言編寫了相應(yīng)計算程序,對圖4所示8輛車作用下6跨簡支梁橋豎向耦合系統(tǒng)進行動力分析。
圖4 車-線-橋耦合系統(tǒng)模型
結(jié)合我國高速鐵路現(xiàn)狀,將二期恒載考慮為附加質(zhì)量,橋梁結(jié)構(gòu)參數(shù)如下:計算跨距L=32 m;彈性模量Eb=34.5 GPa;慣性矩Ib=11.1 m4;單位長質(zhì)量mb=43 628 kg/m;阻尼比ζ=0.02。列車車輛和軌道模型參數(shù)詳見文獻[12]。
2.2 高速下隨機不平順的影響
行車速度分別為250 km/h和300 km/h時,考慮前面生成的5種不平順的影響,表1給出了耦合系統(tǒng)最大動力響應(yīng)??梢?,對橋梁豎向位移響應(yīng)而言,各不平順響應(yīng)同光滑時結(jié)果基本一致,說明不平順影響很小。對橋梁豎向加速度響應(yīng)而言,30~150 m各子區(qū)間不平順響應(yīng)與光滑時結(jié)果相差不大;而其它2種不平順結(jié)果符合一致,且明顯大于光滑時,說明敏感波長子區(qū)間為0.05~30 m。對車體豎向加速度響應(yīng)而言,不平順影響很大,即使90~150 m不平順也可使其增加2倍,說明應(yīng)重視各種波長不平順的控制以提高乘車舒適性。
表1 系統(tǒng)最大動力響應(yīng)
動車和拖車前輪最大輪軌接觸力結(jié)果見圖5,可見0.05~30 m不平順結(jié)果和0.05~150 m不平順結(jié)果一致,明顯大于軌面光滑時結(jié)果,說明敏感波長子區(qū)間為0.05~30 m。圖6為跨中軌道豎向加速度時程結(jié)果,可見車輛依次經(jīng)過時引起了軌道有規(guī)律的強迫振動;30~60 m不平順結(jié)果和軌面光滑時結(jié)果變化趨勢一致,前者影響較小。
2.3 中短波長隨機不平順的影響
由前面分析知,波長區(qū)間0.05~30 m的隨機不平順對系統(tǒng)動力響應(yīng)影響顯著,故有必要對其進行進一步分析,以研究各響應(yīng)的敏感波長區(qū)間。采用任意子波長區(qū)間軌道不平順生成方法,分別生成0.05~1 m、1~5 m、5~12 m、12~20 m和20~30 m波長區(qū)間軌道不平順,對比分析這些不平順對車-線-橋耦合系統(tǒng)動力響應(yīng)的影響。分析發(fā)現(xiàn),軌道不平順對橋梁位移影響很小。橋梁跨中最大豎向加速度對比見圖7,可見大區(qū)間0.05~30 m和子區(qū)間1~5 m不平順影響顯著,且兩者結(jié)果基本一致,說明橋梁加速度敏感波長區(qū)間為1~5 m;0.05~1 m不平順有一定的影響,但高速時這種影響較小。
圖8為跨中軌道最大豎向加速度對比,可見0.05~30 m和0.05~1 m不平順會使響應(yīng)顯著增大,如375 km/h時可使響應(yīng)增加達7.5倍,且兩者結(jié)果相差較小,說明軌道加速度敏感波長區(qū)間為0.05~1 m;1~5 m區(qū)間不平順也會產(chǎn)生一定的影響,其它子區(qū)間不平順影響相對較小。動車最大豎向加速度對比見圖9,可見大區(qū)間0.05~30 m不平順的影響明顯大于其它子區(qū)間,但并不是各子區(qū)間響應(yīng)的疊加,說明應(yīng)采用前者進行最不利響應(yīng)分析;而子區(qū)間0.05~1 m不平順結(jié)果與光滑時符合一致,這是由于車輛懸掛系統(tǒng)較柔,很好地隔離了短波不平順激發(fā)的高頻振動。圖10為動車前輪最大接觸力對比,可見0.05~30 m和1~5 m不平順影響非常大,且兩者結(jié)果相差不大,說明相應(yīng)敏感波長區(qū)間為1~5 m;0.05~1 m不平順也會產(chǎn)生一定的影響,其它子區(qū)間不平順影響相對較小。
圖5 最大動車和拖車輪軌接觸力
圖6 跨中軌道加速度時程
圖7 橋梁最大加速度對比
圖8 軌道最大加速度對比
圖9 車體最大加速度對比
圖10 最大輪軌接觸力對比
圖11 橋梁加速度響應(yīng)時程
橋梁的一階頻率為4.65 Hz,理論共振車速[14]為431 km/h。圖11為431 km/h時,橋梁跨中加速度響應(yīng)時程,可見響應(yīng)隨著通過車輛數(shù)目的增加而急劇增大,形成了明顯的“拍”,出現(xiàn)了共振現(xiàn)象[1];同圖7結(jié)果一致,20~30 m不平順影響很小。
針對軌道高低不平順對車-線-橋耦合振動的影響,提出了任意子波長區(qū)間不平順時域樣本生成方法,對比分析了大波長區(qū)間及其子區(qū)間不平順對系統(tǒng)響應(yīng)的影響,得到如下結(jié)論:
(1) 生成的隨機不平順時域樣本的功率譜與解析值符合一致,且大波長區(qū)間不平順與其包含的所有子區(qū)間不平順的疊加一致,說明提出的軌道不平順生成方法是合理有效的。
(2) 共振車速附近出現(xiàn)了共振現(xiàn)象,可使橋梁響應(yīng)顯著增大,各波長區(qū)間不平順對橋梁位移影響均較小。橋梁豎向加速度主要由波長區(qū)間為1~5 m的不平順控制,0.05~1 m的短波長不平順也能有一些影響,其它波長不平順影響相對較小。
(3) 輪軌接觸力對波長區(qū)間為1~5 m的不平順很敏感,其它波長不平順影響有限。軌道豎向加速度對0.05~1 m的短波長不平順非常敏感,其它波長不平順影響相對較小。車體豎向加速度對0.05~1 m的短波長不平順不敏感,而對其它1~150 m各波長區(qū)間不平順均比較敏感,特別是12~30 m不平順。
(4) 高速鐵路養(yǎng)護和維修中應(yīng)重點關(guān)注5 m以下波長軌道不平順;而為提高乘車舒適性,有必要對各種波長不平順進行控制。
[1]翟婉明,夏禾.列車-軌道-橋梁動力相互作用理論與工程應(yīng)用[M].北京:科學出版社,2011.
[2] 林玉森,李小珍,強士中.軌道不平順激勵下高速鐵路橋上列車走行性研究[J].鐵道學報,2005,27(6):96-100.
[3] 楊宏印,張海龍,陳志軍,等.考慮徐變影響的車-線-橋耦合有限元分析[J].中國鐵道科學,2013,34(6):13-20.
[4] 陳春俊,李華超.頻域采樣三角級數(shù)法模擬軌道不平順信號[J].鐵道學報,2006,28(3):38-42.
[5] 翟婉明.車輛-軌道耦合動力學[M].3版.北京:科學出版社,2007.
[6] 劉寅華,李芾,黃運華.軌道不平順數(shù)值模擬方法[J].交通運輸工程,2006,6(1):29-33.
[7] 程小平,練松良,陳敏敏.軌道交通車輛車輪踏面短波不平順測試與分析[J].石家莊鐵道學院學報,2010, 23(1):6-10.
[8] Zhai W M, Wang K Y, Cai C B. Fundamentals of vehicle-track coupled dynamics[J]. Vehicle System Dynamics, 2009, 47: 1349-1376.
[9] 楊宏印,陳志軍,張海龍.彈性接觸下車輛-變截面梁橋豎向振動分析[J].公路,2013(10):97-101.
[10]王瀾. 軌道結(jié)構(gòu)隨機振動理論及其在軌道結(jié)構(gòu)減振中的應(yīng)用[D]. 北京: 鐵道科學研究院, 1988.
[11] Yang H Y, Chen Z J, Zhang H L. Vibration of train-rail-bridge interaction considering rail Irregularity with arbitrary wavelength[J]. International Journal of Engineering, 2015, 28(4): 516-522.
[12] 楊宏印. 高速鐵路車-線-橋動力分析理論及橋上軌道激勵影響研究[D]. 武漢: 華中科技大學, 2014.
[13] Yang H Y, Chen Z J, Zhang H L. Dynamic analysis of train-rail-bridge interaction considering concrete creep of a multi-span simply supported bridge[J]. Advances in Structural Engineering, 2014, 17(5): 709-720.
[14] 王慧東,馬其森.韓江大橋(48+80+80+88+48) m道岔連續(xù)梁車橋耦合動力響應(yīng)分析[J]. 石家莊鐵道大學學報: 自然科學版, 2015, 28(1): 12-16.
Influence of Track Irregularity on Train-rail-bridge Vibration Based on Sub-interval
Yang Hongyin1, Chen Zhijun1, Peng Xumin2
(1.China Railway Bridge Science Research Institute Ltd., Wuhan 430034, China;2. School of Civil Engineering & Mechanics, Huazhong University of Science & Technology, Wuhan 430074, China)
The generation method for rail irregularity with arbitrary wavelength is proposed, and the characteristics of the samples in time domain are discussed. The vibration equations of the interaction system are established, and the comparison of the influences of large wavelength and its sub-intervals irregularity is analyzed. The sensitive wavelengths of dynamic responses are investigated. The results demonstrate that the power spectra of rail irregularity samples are consistent with analytical values, and the generation method is reasonable and effective. The irregularity has little effect on bridge displacement, and the sensitive wavelength interval for bridge acceleration is 1~5 m. The irregularity with wavelength 0.05~1 m should be the main excitation source for the rail acceleration while it has little effect on the car body acceleration. The irregularity with wavelength less than 5m should be particularly controlled in track maintenance, and to improve the riding comfort, the irregularity with any wavelength should be controlled.
track irregularity; train-rail-bridge coupled vibration; sensitive wavelength; resonance; dynamic response
2015-07-20 責任編輯:車軒玉
10.13319/j.cnki.sjztddxxbzrb.2016.02.02
國家自然科學基金(51078164); 湖北省自然科學基金(2013CFA135)
楊宏印(1986-),男,博士,工程師,E-mail: yanghongyin@hust.edu.cn
TU311.3; U213.2
A
2095-0373(2016)02-0006-06
楊宏印,陳志軍,彭旭民.基于子區(qū)間的車-線-橋振動不平順影響分析[J].石家莊鐵道大學學報:自然科學版,2016,29(2):6-11.