毛建紅
(中鐵第五勘察設計院集團有限公司,北京102600)
隨著我國大量已建成高速鐵路的投入運營,服役過程中的軌道病害問題開始凸顯出來。在高速列車的動荷載作用下,無砟軌道下部路基發生沉降,軌面高程改變,進一步惡化高速鐵路的平順性。為了克服這一難題,對于列車作用下路基沉降變形的理論國內外專家進行了廣泛的研究。高廣運[1]推導了一種列車動荷載下基于蠕變本構方程的長期沉降計算方法,結合分層總和法計算出列車運行引起的長期沉降。陶明安[2]在動力有限元中引入等效移動列車荷載,建立了一種高速列車移動荷載作用下飽和軟黏土地基的長期沉降計算方法。J.Kennedy 等[3]提出了聚合物加固道砟后的沉降問題,發現聚合物加固道砟能顯著地降低軌道沉降,并通過在實際鐵路現場的應用證實了實驗室的研究。
針對無砟軌道路基沉降這一病害,采用高聚物注漿抬升修復技術對無砟軌道進行注漿修復后,軌道還將長期服役,修復后的路基相當于基床表層增加了一層高聚物碎石混合物層,增加這一新結構層后路基的沉降趨勢是極其值得研究的。本文先推導基于伯格斯模型的低頻循環荷載軸向應變公式,利用實驗進行公式中力學參數的識別。建立三維中心差分法動力學模型計算修復后路基內各層的內力,高聚物碎石混合物層的沉降采用本文推導公式計算,路基內其他層采用既有公式計算,最終利用分層總和法計算修復后路基的總沉降,研究修復后無砟軌道的沉降問題。
根據王軍保教授在文獻[4]中的實驗結論,推導出高聚物碎石混合物在無側限壓力下的單軸循環荷載軸向變形公式。
根據彈性力學理論:三維應力狀態下應力張量可以分解為球應力張量和偏應力張量,偏應力張量表達式:
軸向偏應力為:

式中,σ*=σ1-σ3為主應力差,對試塊采用單軸無側限加載,所以σ1=psin(ωt),σ*(t)=psin(ωt);p 為應力幅值;t 為時間;ω 為角速度。
1)低頻荷載下Maxwell 模型軸向蠕變方程:

式中,G1,η1分別為 Maxwell 體剪切模量和黏滯系數;K 為體積模量;E1為Maxwell 體彈性模量;q 為恒定應力。
2)正弦荷載下Kelvin 模型軸向蠕變方程:

式中,G2,η2分別為 Kelvin 體剪切模量和黏滯系數;E2為Maxwell 體彈性模量。
有ε(t)=εM(t)+ εK(t),故Burgers 模型在低頻循環荷載作用下的蠕變方程即軸向應變方程為:

本節引用浙江大學邊學成教授團隊的實驗結論[5]識別公式(3)中代表高聚物碎石混合物層力學性能的 E1,η1,E2,η24 項參數。本文以實驗數據和擬合結果的誤差平方和為目標函數,將其轉化為最優化問題,使其迭代收斂速度加快,快速識別參數。為解決這一最優化問題引入多種粒子群算法[6]。根據文獻[7]中的結論,選擇具有代表性的3 個值繪圖對比,發現曲線基本相同,將時間軸集中在1s 中內時發現E1=6.248 622 9×1013時曲線在最下方應變較小,而E1=1.041 034 8×1015時在最上方應變較大,與文獻[7]結論相同。
采用FLAC3D 有限差分軟件建立無砟軌道動力學模型,軌道板為CRTSⅡ型板。模型如圖1 所示。

圖1 中心差分法模型
根據模型計算繪制圖2,從圖2 可以看出隨著速度的增大,中軸線和軌底處的高聚物碎石混合物層的動偏應力顯著增大,而基床1m 以下的動偏應力增大并不顯著。在不同速度下,在高聚物碎石混合層中軌底處的動偏應力皆大于中軸線處,而其他軌底處動偏應力皆小于中軸線處。除去高聚物碎石層的動偏應力可以發現,隨著路基垂向深度的增加動偏應力不斷增大,這一結果符合文獻[8]中的結論。

圖2 動偏應力示意圖
研究不同行車速度對路基沉降的影響,速度取值分別為150km/h、250km/h、350km/h。從圖 3 中可以看出,在初始時間內,沉降的速度近似相同,之后350km/h 速度的沉降量相比于150km/h、250km/h 的結果顯著增大,在進入到第五年后不同速度下沉降的增長速度基本相同,但350km/h 的沉降量明顯大于150km/h、250km/h。可見采用較大的行車速度會加速路基的沉降。綜上所述,路基面大部分沉降量發生在運營初期,之后增長緩慢。

圖3 路基面沉降示意圖
本文研究了注漿抬升修復條件下無砟軌道路基沉降問題。通過引入Burgers 模型推導了高聚物碎石混合物動荷載作用下軸向應變公式。研究發現,該公式能合理適用于動荷載作用下高聚物碎石混合物層的應變,在進行參數識別時發現雜交粒子群算法能夠有效地解決公式中 E1、E2、η1、η24 個未知參數的識別問題。在采用分層總和法預測總沉降量時發現,隨著時間的增加軌面沉降量也不斷增加,初始1~10 年時沉降速度較快而10~25 年的沉降速度緩慢,軌面高度中軸線處的沉降量皆大于鋼軌處。在初始時間內,沉降的速度近似相同,之后采用較大的行車速度會加速路基的沉降。