秦佳良 謝 毅 周冠南 周 濤 王建立 劉林芽 李秋義
(1.華東交通大學(xué), 南昌 330013;2.中鐵二院工程集團(tuán)有限責(zé)任公司, 成都 610031;3.中國(guó)鐵建大橋工程局集團(tuán)有限公司, 天津 300300;4.中鐵五局集團(tuán)有限公司, 長(zhǎng)沙 410007;5.隔而固(青島)振動(dòng)控制有限公司, 山東 青島 266108;6.中鐵第四勘察設(shè)計(jì)院集團(tuán)有限公司, 武漢 430063)
近年來(lái),隨著我國(guó)高速鐵路的快速發(fā)展,我國(guó)自主研發(fā)的CRTSⅢ型板式無(wú)砟軌道由于高平順、高穩(wěn)定性和少維修等諸多優(yōu)點(diǎn),在鐵路上廣泛應(yīng)用。
對(duì)于車(chē)輛-無(wú)砟軌道結(jié)構(gòu)的振動(dòng)問(wèn)題,國(guó)內(nèi)外學(xué)者進(jìn)行了大量的研究。Guigou-Carter等[1]人開(kāi)發(fā)了一種二維的板式軌道系統(tǒng)預(yù)測(cè)模型,軌道系統(tǒng)由鋼軌、軌道板和混凝土底座板組成,分析了軌枕墊的動(dòng)態(tài)剛度對(duì)軌道系統(tǒng)隔振效果的影響;Khajehdazfuly[2]建立了高速列車(chē)-CRTSⅡ板式無(wú)砟軌道垂向數(shù)值模型,研究了諧波不平順下軌道墊板剛度對(duì)輪軌力的影響;翟婉明等[3]運(yùn)用車(chē)輛-軌道耦合動(dòng)力學(xué)理論,建立了高速列車(chē)與板式軌道相互作用的垂向動(dòng)力學(xué)模型,針對(duì)日本新干線(xiàn)板式軌道,分析了CA砂漿的剛度和阻尼變化對(duì)軌道結(jié)構(gòu)動(dòng)力特性的影響;蔡成標(biāo)[4]建立了路基上無(wú)砟軌道空間耦合動(dòng)力學(xué)模型,考慮了垂、橫向運(yùn)動(dòng),研究了遂渝線(xiàn)無(wú)砟軌道綜合試驗(yàn)段路基上板式軌道的動(dòng)力學(xué)特性;向俊等[5-6]基于彈性系統(tǒng)動(dòng)力學(xué)總勢(shì)能不變值原理,提出了一種橫向有限條與無(wú)砟軌道板段單元分析模型,針對(duì)CRTSⅠ型板式軌道研究了車(chē)速和軌道高低不平順對(duì)系統(tǒng)豎向振動(dòng)響應(yīng)的影響;雷曉燕等[7-8]提出了一種新型板式軌道單元,建立了車(chē)輛-板式軌道-路基耦合系統(tǒng)的動(dòng)力學(xué)模型,針對(duì)CRTS Ⅱ型板式無(wú)砟軌道研究了板式軌道結(jié)構(gòu)及過(guò)渡段的動(dòng)力學(xué)行為,評(píng)估了列車(chē)速度和軌道剛度對(duì)軌道振動(dòng)的影響。對(duì)于CRTSⅢ型板式無(wú)砟軌道,相關(guān)研究主要集中在結(jié)構(gòu)的設(shè)計(jì)以及施工方面,而對(duì)于CRTSⅢ型板式無(wú)砟軌道動(dòng)力學(xué)的研究還較少。
基于以上研究成果,本文采用有限元法,建立了列車(chē)-減振型CRTSⅢ板式軌道耦合系統(tǒng)豎向動(dòng)力響應(yīng),以中國(guó)高速鐵路軌道不平順譜作為激勵(lì)源,分析列車(chē)荷載作用下減振型CRTSⅢ板式軌道的動(dòng)力特性,并對(duì)減振層剛度變化對(duì)軌道結(jié)構(gòu)振動(dòng)響應(yīng)的影響進(jìn)行了研究。
采用有限元法建立車(chē)輛-CRTS Ⅲ型板軌道耦合系統(tǒng)豎向振動(dòng)模型,應(yīng)滿(mǎn)足以下假設(shè):
(1)車(chē)輛系統(tǒng)簡(jiǎn)化為附有二系彈簧阻尼的整車(chē)模型,軌道系統(tǒng)簡(jiǎn)化為三層板式軌道單元模型,輪軌之間采用非線(xiàn)性彈性接觸。
(2)鋼軌視為離散點(diǎn)支承的二維Euler梁,扣件墊板的彈性系數(shù)和阻尼系數(shù)分別用Kp和Cp表示。
(3)軌道板和自密實(shí)混凝土視為連續(xù)粘彈性支承的二維梁?jiǎn)卧?,減振層的彈性系數(shù)和阻尼系數(shù)分別用Km和Cm。
(4)底座板也視為連續(xù)粘彈性支承的二維梁?jiǎn)卧?,底座板下部基礎(chǔ)支承的彈性系數(shù)和阻尼系數(shù)分別用Ki和Ci。
路基上CRTS Ⅲ型板式無(wú)砟軌道系統(tǒng)結(jié)構(gòu)如圖1所示,基于上述的假設(shè)可以得到建立的三層軌道單元模型,如圖2所示。
圖1 路基上CRTS Ⅲ型板式軌道圖
圖2 三層梁有限元模型圖
圖2中,v1和v4表示鋼軌的豎向位移;θ1和θ4表示鋼軌的轉(zhuǎn)角;v2和v5表示軌道板+自密實(shí)混凝土的豎向位移;θ2和θ5表示軌道板+自密實(shí)混凝土的轉(zhuǎn)角;v3和v6表示底座板的豎向位移;θ3和θ6表示底座板的轉(zhuǎn)角。
定義CRTS Ⅲ型板軌道單元節(jié)點(diǎn)位移向量為:
(1)
板式無(wú)砟軌道單元的剛度矩陣可表示為:
(2)
板式無(wú)砟軌道單元的阻尼矩陣可表示為:
(3)
板式無(wú)砟軌道單元的質(zhì)量矩陣可表示為:
(4)
無(wú)砟軌道單元等效節(jié)點(diǎn)荷載向量為:
(5)
其中:Fi(i=1~4)——輪軌接觸力,可由輪軌接觸關(guān)系求得。
根據(jù)上述軌道單元模型,利用有限元“對(duì)號(hào)入座法”組集軌道單元?jiǎng)偠?、阻尼和質(zhì)量矩陣以及荷載向量,即可得到軌道系統(tǒng)振動(dòng)方程為:
(6)
車(chē)輛系統(tǒng)簡(jiǎn)化為一個(gè)附有二系懸掛彈簧阻尼的整車(chē)模型,由1個(gè)車(chē)體、2個(gè)轉(zhuǎn)向架和4個(gè)輪對(duì)組成??紤]了車(chē)體的沉浮和點(diǎn)頭運(yùn)動(dòng),前后構(gòu)架的沉浮和點(diǎn)頭運(yùn)動(dòng),以及4個(gè)輪對(duì)的垂向振動(dòng)共10個(gè)自由度。定義車(chē)輛單元節(jié)點(diǎn)位移向量為:
(7)
其中:vc——車(chē)體沉浮振動(dòng)的豎向位移;
θc——車(chē)體點(diǎn)頭振動(dòng)的角位移;
vti和θti(i=1,2)——分別為前、后轉(zhuǎn)向架沉浮振動(dòng)的豎向位移和點(diǎn)頭振動(dòng)的角位移;
vwi(i=1,2,3,4)——第i個(gè)車(chē)輪的豎向位移。
由Hamilton原理,可得到車(chē)輛系統(tǒng)振動(dòng)方程為:
(8)
Mv——車(chē)輛系統(tǒng)的質(zhì)量;
Kv——車(chē)輛系統(tǒng)的剛度;
Cv——車(chē)輛系統(tǒng)的阻尼矩陣。這些矩陣的顯式表達(dá)式參見(jiàn)文獻(xiàn)[9]。
Qv為車(chē)輛系統(tǒng)的廣義力向量:
(9)
式中:g——重力加速度;
Fi(i=1~4)——輪軌接觸力,可用輪軌非線(xiàn)性Hertz接觸公式求得。
(10)
(11)
(12)
仿真分析中,車(chē)輛選用高速動(dòng)車(chē)CRH3,具體參數(shù)如表1所示,軌道選擇減振型CRTSⅢ板式無(wú)砟軌道,其參數(shù)如表2所示。
表1 CRH3高速客車(chē)的動(dòng)力學(xué)參數(shù)表
表2 減振型CRTSⅢ板式軌道參數(shù)表
減振型CRTSⅢ板式無(wú)砟軌道的特點(diǎn)是在自密實(shí)混凝土和底座板之間鋪設(shè)1層橡膠減振墊層。為研究不同減振層剛度和阻尼變化對(duì)軌道結(jié)構(gòu)變形的影響,本文分別取減振層剛度km為100 kN/mm、400 kN/mm、700 kN/mm、1 000 kN/mm。軌道不平順類(lèi)型選取我國(guó)的高速鐵路無(wú)砟軌道不平順譜,計(jì)算減振層剛度變化對(duì)車(chē)輛、軌道系統(tǒng)振動(dòng)的影響。
減振層剛度變化對(duì)車(chē)體和輪軌力影響的對(duì)比如圖3所示。從圖3可以看出,隨著減振層剛度的增大,車(chē)體加速度有所減小但變化不明顯,當(dāng)減振層剛度從100 kN/mm增加到 1 000 kN/mm,車(chē)體加速度減小 0.001 04 m/s2,影響較??;輪軌力也有所減小,當(dāng)減振層剛度從100 kN/mm增加到 1 000 kN/mm,車(chē)體加速度減小 0.001 04 m/s2,影響較小。
圖3 減振層剛度對(duì)車(chē)體和輪軌力的影響圖
減振層剛度變化對(duì)鋼軌動(dòng)力特性影響的對(duì)比如圖4所示。從圖4可以明顯看出,隨著減振層剛度的增大,鋼軌位移整體有所減小,當(dāng)減振層剛度為100 kN/mm時(shí),鋼軌最大位移為1.21 mm;當(dāng)減振層剛度變?yōu)?00 kN/mm、700 kN/mm和 1 000 kN/mm時(shí),鋼軌最大位移為0.98 mm、0.94 mm和0.93 mm,鋼軌最大位移降幅為22.93%、3.64%、1.48%,降低的速率隨著減振層剛度的增大而減小。隨著減振層剛度的增大,鋼軌加速度也有所減小,當(dāng)減振層剛度為100 kN/mm時(shí),鋼軌最大加速度為8.50 m/s2;當(dāng)減振層剛度變?yōu)?00 kN/mm、700 kN/mm和 1 000 kN/mm時(shí),鋼軌最大加速度為8.44 m/s2、8.35 m/s2和8.28 m/s2,鋼軌加速度降幅為0.72%、1.11%、0.79%。
圖4 減振層剛度對(duì)鋼軌動(dòng)力特性的影響圖
減振層剛度變化對(duì)軌道板動(dòng)力特性影響的對(duì)比如圖5所示。從圖5可以明顯看出,隨著減振層剛度的增大,軌道板位移整體減小,當(dāng)減振層剛度為100 kN/mm時(shí),軌道板最大位移為0.789 mm;當(dāng)減振層剛度變?yōu)?00 kN/mm、700 kN/mm和 1 000 kN/mm時(shí),軌道板最大位移為0.588 mm、0.556 mm和0.5 44 mm,軌道板最大位移降幅為34.36%、5.58%、2.20%,降低的速率有所減緩。隨著減振層剛度的增大,軌道板加速度也有所減小,當(dāng)減振層剛度為100 kN/mm時(shí),軌道板最大加速度為3.01 m/s2;當(dāng)減振層剛度變?yōu)?00 kN/mm、700 kN/mm和 1 000 kN/mm時(shí),鋼軌最大加速度為2.46 m/s2、2.23 m/s2和2.10 m/s2,軌道板加速度降幅為21.95%、10.22%、6.35%。
圖5 減振層剛度對(duì)軌道板動(dòng)力特性的影響圖
減振層剛度變化底座板動(dòng)力特性影響的對(duì)比如圖6所示。從圖6可以看出,隨著減振層剛度的增大,底座板位移不斷增大但變化較小,減振層剛度從100 kN/mm增加到 1 000 kN/mm,底座板位移增大 0.001 09 mm,影響較小。底座板加速度隨著減振層剛度的增大也不斷增大,底座板的加速度增幅為15.08%、4.03%、0.69%。
圖6 減振層剛度對(duì)底座板動(dòng)力特性的影響圖
本文基于有限元理論,建立了列車(chē)-CRTS Ⅲ型減振板式軌道耦合動(dòng)力學(xué)模型,分析了不同減振層剛度對(duì)軌道結(jié)構(gòu)動(dòng)力特性的影響[10],主要結(jié)論如下:
(1)減振層剛度的變化對(duì)車(chē)體的振動(dòng)影響較小,對(duì)鋼軌、軌道板和底座板的振動(dòng)影響較大。
(2)減振層剛度的增大能有效的降低減振層以上結(jié)構(gòu)的振動(dòng),引起下部結(jié)構(gòu)的劇烈振動(dòng)。
(3)從設(shè)計(jì)的角度來(lái)看,減振層剛度不宜設(shè)置的太大,否則將削弱板式軌道的彈性,不利于軌道結(jié)構(gòu)減振;同時(shí)考慮減振層上部和下部結(jié)構(gòu)的動(dòng)力效應(yīng),建議減振層剛度應(yīng)在 400~700 kN/mm。