高玄濤

(中鐵第一勘察設(shè)計院集團(tuán)有限公司，西安 710043)

近年來，我國軌道交通得到大力發(fā)展。隨著我國軌道交通的大規(guī)模建設(shè)，近接工程數(shù)量越來越多，施工難度也越來越大。特別是隧道下穿高速鐵路施工時，不可避免地對周圍土體產(chǎn)生擾動，引起周圍地層損失及路基沉降，造成高速鐵路軌道的縱、橫向不平順，從而影響列車的運(yùn)營安全［1］。以某隧道下穿高速鐵路項(xiàng)目為研究背景，主要采用數(shù)值模擬的方法對下穿隧道施工至不同位置時列車動載作用下地層動力響應(yīng)進(jìn)行研究分析，為下穿隧道施工提供理論指導(dǎo)。

1 工程概況

該隧道下穿高速鐵路A(即將投入運(yùn)營)、客運(yùn)專線B(已運(yùn)營)、客運(yùn)專線C(在建)段，需與在建的工程同步完成施工。隧道全長300 m，隧道內(nèi)為單面坡，隧道開挖拱頂?shù)缴喜胯F路軌頂面最大埋深約27.3 m。隧道下穿段與地表高速鐵路A、客運(yùn)專線B及客運(yùn)專線C在平面位置上呈斜交，斜交角度約30°。

高速鐵路A本段地基采用CFG樁加固，樁徑0.5 m，樁長 10.5 m，樁底高程約 13.275 m，樁端一般進(jìn)入強(qiáng)風(fēng)化層不少于0.5 m，樁間距1.80 m，按正方形布置，樁頂設(shè)1.0 m擴(kuò)大樁頭，樁頂設(shè)0.6 m厚墊層(0.2 m厚中粗砂+0.4 m厚碎石)，墊層內(nèi)鋪設(shè)二層雙向經(jīng)編高強(qiáng)土工格柵;土工格柵極限抗拉強(qiáng)度不小于100 kN/m。地基表層挖除換填0.5 m厚C組填料后再打樁。

客運(yùn)專線B本段地基采用CFG樁加固，樁徑0.5 m，樁長9.5～6.5 m，CFG 樁樁端進(jìn)入硬層不少于1.0 m;樁間距 1.8 m，按正方形布置，樁頂設(shè)直徑為1.0 m擴(kuò)大樁頭;樁頂設(shè)0.6 m厚墊層(0.2 m厚中粗砂+0.4 m厚碎石)，墊層內(nèi)鋪設(shè)二層雙向經(jīng)編高強(qiáng)土工格柵;土工格柵極限抗拉強(qiáng)度不小于100 kN/m。地基表層挖除換填0.5 m厚C組填料后再打樁。

客運(yùn)專線C本段地基采用CFG樁加固，樁徑0.5 m，樁間距1.6 m，樁長 3.0 ～9.0 m(樁長嵌入全風(fēng)化層內(nèi)不小于1.0 m)，按正三角形布置。CFG樁樁頂采用1 mx1 m擴(kuò)大樁頭，樁頂設(shè)0.5 m厚碎石墊層，墊層內(nèi)鋪設(shè)1層雙向經(jīng)編土工格柵，土工格柵抗拉強(qiáng)度不小于250 kN/m。

下穿隧道與地表線路路基相對位置及縱斷面分別見圖1、圖 2。

圖1 工程平面位置(單位:m)

圖2 工程縱斷面(單位:m)

2 地質(zhì)概況

隧道處地貌為秦淮河及其階地、崗地、坳谷以及剝蝕丘陵，隧道洞身位于第四系地層中。根據(jù)地質(zhì)資料揭露，上部為2～5 m厚人工填土，下部為第四系上更新統(tǒng)粉質(zhì)黏土層，硬塑，厚度0～5 m，具有弱膨脹性，下伏基巖為白堊系砂巖或侏羅系凝灰?guī)r、砂巖，節(jié)理裂隙發(fā)育，全風(fēng)化～弱風(fēng)化，地下水主要賦存于土層與風(fēng)化層、強(qiáng)風(fēng)化與弱風(fēng)化層間及巖層裂隙內(nèi)，水量較豐富。隧道拱頂以上0～3 m范圍為弱風(fēng)化巖，隧道所處地質(zhì)根據(jù)圍巖情況判定為Ⅲ～Ⅳ級，根據(jù)隧道穿越地段地表建筑物的重要性及地下水情況，將隧道按為Ⅴ級型式加強(qiáng)支護(hù)。如圖3所示。

圖3 地質(zhì)柱狀圖(單位:m)

3 動力分析本構(gòu)模型

土的動力本構(gòu)關(guān)系是極其復(fù)雜的，它在不同的荷載條件、土性條件下會表現(xiàn)出極不相同的動力本構(gòu)特性。目前具體建立的動力本構(gòu)模型已達(dá)數(shù)十種，大致可以分為線彈性模型、黏彈性模型、彈塑性模型、邊界面模型、內(nèi)時模型等［2］。

根據(jù)彈塑性理論，土體變形可以分為彈性變形和塑性變形兩部分，其中彈性變形可以應(yīng)用廣義胡克定理計算，塑性變形可以應(yīng)用塑性增量理論計算。塑性增量理論包含3個基本要素:屈服準(zhǔn)則、流動法則和硬化定律，分別規(guī)定了塑性應(yīng)變增量的產(chǎn)生條件、方向和大小。目前基于不同的基本要素已經(jīng)提出了多種彈塑性本構(gòu)模型，如Drucker-Prager模型、多重屈服面模型(Nested Surfaces Model)、邊界面模型(Bounding Surfaces Model)等［2］。

事實(shí)上，目前還沒有一種本構(gòu)模型能夠模擬各種情況下土的動力非線性特征，并具有實(shí)際應(yīng)用中必要的簡單性。因此針對具體的問題應(yīng)該選擇合理而簡單的本構(gòu)模型［2］。

ABAQUS是國際上最先進(jìn)的大型非線性有限元計算分析軟件之一，具有強(qiáng)健的非線性計算功能。其提供的非線性模型有Mohr-Coulomb模型、擴(kuò)展的Drucker-Prager模型、修正Drucker-Prager帽蓋模型、(修正)劍橋模型。其中Mohr-Coulomb模型在巖土工程中應(yīng)用最廣泛，模擬結(jié)果與實(shí)際也較為吻合。因此計算中采用 Mohr-Coulomb 模型［3］。

4 有限元計算模型的建立

采用三維有限元模型進(jìn)行計算。根據(jù)工程實(shí)際及圣維南原理，參考已有研究成果進(jìn)行建模［4-7］，隧道頂部以上覆土為23.5 m，隧道底部以下為29.7 m，左右邊界分別離隧道中線50 m，縱向延伸50 m。通過ABAQUS程序建立的計算分析模型如圖4所示。

圖4 下穿隧道計算分析模型

4.1 本構(gòu)模型

土體的非線性本構(gòu)模型:采用了彈塑性物理模型。襯砌單元采用強(qiáng)度等級為C35的鋼筋混凝土，計算時采用Mohr-Coulomb本構(gòu)模型并將混凝土與鋼筋合為一體考慮。為了建模及劃分網(wǎng)格方便，將計算區(qū)域不同地層簡化為同一地層，土體參數(shù)進(jìn)行相應(yīng)轉(zhuǎn)換。根據(jù)下穿隧道段地質(zhì)勘察資料，計算過程中土體的材料參數(shù)取值如表1所示。

表1 材料物理力學(xué)參數(shù)

4.2 邊界條件［8-10］

本工程通過對固定邊界、黏性邊界條件及無限元邊界條件進(jìn)行計算比較分析，計算結(jié)果相差不大，一方面是由于列車振動荷載不是很大，另一方面模型選取的比較大，列車振動荷載傳至邊界處時已衰減較多。而從計算所需時間上看，固定邊界所需時間要少，故選用固定邊界進(jìn)行計算。

固定邊界條件:土層底部完全固定;左右兩側(cè)面限制水平方向的位移，豎向自由;前后兩面限制軸線方向的位移，豎向自由;地面完全自由。

5 動力計算及結(jié)果分析

本文針對無砟軌道的結(jié)構(gòu)特點(diǎn)，采用彈性地基梁板模型分析了CRTSⅠ型板式無砟軌道結(jié)構(gòu)動力響應(yīng)，給出了CRTSⅠ型板式無砟軌道結(jié)構(gòu)在CRH2型列車以時速350 km經(jīng)過時軌下壓力時程曲線，確定了列車振動過程中的激勵荷載［11］。如圖5所示。

圖5 列車輪軌激振力時程曲線

本文主要計算了隧道開挖到不同位置時列車以時速350 km通過，研究地層動力響應(yīng)。計算過程中分別在地層表面、地層不同深度處及隧道拱頂設(shè)置監(jiān)測點(diǎn)，監(jiān)測該處位移變化情況。

5.1 地層表面監(jiān)測點(diǎn)位移分析

隧道分別開挖至15、25 m和35 m時，施加列車振動荷載，地層表面各監(jiān)測點(diǎn)在列車動載作用下的豎向位移時程曲線如圖6所示。

圖6 隧道開挖至不同位置時施加動載地表監(jiān)測點(diǎn)位移時程曲線

隧道開挖至不同位置施加動載地表監(jiān)測點(diǎn)最大豎向位移如表2所示。

表2 地表監(jiān)測點(diǎn)最大豎向位移 mm

通過對隧道開挖至不同部位施加動載地表監(jiān)測點(diǎn)位移時程曲線及最大位移值分析可知。

(1)在單次列車動載作用下地表監(jiān)測點(diǎn)產(chǎn)生的豎向位移值與隧道開挖引起的沉降量相比很小，可以忽略不計(隧道開挖引起的最大沉降量為2.7 mm)。但是根據(jù)國內(nèi)外工程實(shí)際在長期高速列車振動荷載的作用下，地層產(chǎn)生的沉降量是不能忽視的;

(2)隨著一組組車輪的滾過，監(jiān)測點(diǎn)位移時程曲線具有明顯的車輪滾過效應(yīng);當(dāng)車輪壓在監(jiān)測點(diǎn)時，該監(jiān)測點(diǎn)出現(xiàn)位移峰值，車輪離開監(jiān)測點(diǎn)后，位移逐漸減小;當(dāng)列車駛離后，位移逐漸回彈恢復(fù);

(3)當(dāng)?shù)谝还?jié)車廂經(jīng)過時，z=15 m處監(jiān)測點(diǎn)產(chǎn)生的位移最大;隨后，z=25 m處即下穿隧道與高速鐵路路基交叉點(diǎn)處監(jiān)測點(diǎn)位移最大，因此，在設(shè)計與施工時要注意對下穿隧道與高速鐵路路基交叉部分加強(qiáng)支護(hù);

(4)隧道中點(diǎn)前監(jiān)測點(diǎn)隨著隧道開挖的不斷進(jìn)行，施加動載產(chǎn)生的位移值在逐漸減小;隧道中點(diǎn)及中點(diǎn)后監(jiān)測點(diǎn)則隨著隧道開挖的不斷進(jìn)行，施加動載產(chǎn)生的位移值在不斷增大。

5.2 隧道拱頂監(jiān)測點(diǎn)位移分析

隧道分別開挖至15、25 m和35 m時，施加列車振動荷載，隧道拱頂各監(jiān)測點(diǎn)在列車動載作用下的豎向位移時程曲線如圖7所示。

圖7 隧道開挖至不同位置時施加動載隧道拱頂監(jiān)測點(diǎn)位移時程曲線

隧道開挖至不同位置施加動載隧道拱頂監(jiān)測點(diǎn)最大豎向位移值如表3所示。

表3 拱頂監(jiān)測點(diǎn)最大位移值 mm

通過對隧道開挖至不同部位施加動載隧道拱頂監(jiān)測點(diǎn)位移時程曲線及最大位移值分析可知:

(1)隧道拱頂位移時程曲線表現(xiàn)出和地表位移時程曲線同樣的特性;

(2)隧道開挖至某一處施加動載，該處拱頂產(chǎn)生的位移值最大，因此，隧道開挖后，特別是上臺階開挖后，要及時施加支護(hù);

(3)隨著開挖的不斷推進(jìn)，開挖面即掌子面處拱頂位移值在不斷增大。

5.3 地層不同深度監(jiān)測點(diǎn)位移分析

為研究在列車動載作用下，地層不同深度處動力反應(yīng)，取隧道開挖面即掌子面處地層不同深度監(jiān)測點(diǎn)來進(jìn)行研究，監(jiān)測點(diǎn)布置如圖8所示。隧道開挖至不同位置時施加動載，在動載作用下各監(jiān)測點(diǎn)的位移時程曲線如圖9所示。

圖8 不同深度監(jiān)測點(diǎn)分布示意

圖9 隧道開挖至不同位置時施加動載地層不同深度監(jiān)測點(diǎn)位移時程曲線

隧道開挖至不同位置施加動載地層不同深度監(jiān)測點(diǎn)最大位移值如表4所示。

通過對隧道開挖至不同部位施加動載地層不同深度監(jiān)測點(diǎn)位移時程曲線及最大位移值分析可知:

(1)隨著地層深度的增加，在動載作用下產(chǎn)生的位移不斷減小，地表處監(jiān)測點(diǎn)位移最大，隧道拱頂處位移最小。這說明振動荷載在地層中傳播時，由于地層阻尼等因素使振動波能量不斷減小。

(2)地層動力響應(yīng)從地表到隧道拱頂衰減速度逐漸減小。

(3)隧道開挖至25 m處施加動載，地層不同深度各監(jiān)測點(diǎn)的位移值較大。

5.4 地層應(yīng)力響應(yīng)分析

為分析在列車通過時地層的豎向應(yīng)力變化情況，通過對地層監(jiān)測點(diǎn)(圖8)進(jìn)行監(jiān)測分析，得到不同深度處地層在列車振動荷載作用下的豎向應(yīng)力時程曲線，如圖10所示。

表4 地層不同深度監(jiān)測點(diǎn)最大位移值 mm

隧道開挖至不同位置施加動載地層不同深度監(jiān)測點(diǎn)最大應(yīng)力值如表5所示。

表5 隧道開挖至不同位置施加動載地層不同深度監(jiān)測點(diǎn)應(yīng)力最大值 Pa

圖10 隧道開挖至不同位置時施加動載地層不同深度監(jiān)測點(diǎn)應(yīng)力時程曲線

通過對隧道開挖至不同部位施加動載地層不同深度監(jiān)測點(diǎn)應(yīng)力時程曲線及應(yīng)力最大值分析可知:

(1)隨著地層深度的增加，在動載作用下產(chǎn)生的豎向應(yīng)力不斷減小，地表處監(jiān)測點(diǎn)沉降最大，隧道拱頂處沉降最小。這說明振動荷載在地層中傳播時，由于地層阻尼等因素使振動波能量不斷減小。

(2)地層動力響應(yīng)從地表到隧道拱頂衰減速度逐漸減小，但是隧道拱頂附近衰減較快，這主要可能是由于隧道拱頂?shù)膲毫靶?yīng)造成的。

(3)隨著隧道開挖的不斷推進(jìn)，開挖面即掌子面處不同深度地層在列車動載作用下的豎向應(yīng)力在不斷減小，這主要是因?yàn)樗淼篱_挖深度越深，高速鐵路路基與下穿隧道交叉處的“臨空面”面積不斷增加，地層在列車動載作用下向隧道凈空發(fā)生位移，使應(yīng)力在一定程度上釋放。

6 結(jié)語

通過對隧道施工至不同位置時列車以時速350 km通過，研究地表及地層不同位置處位移值的變化，結(jié)論如下。

(1)在單次列車動載作用下產(chǎn)生的位移值與隧道開挖引起的沉降量相比很小，可以忽略不計。但是根據(jù)國內(nèi)外工程實(shí)際，在長期高速列車振動荷載的作用下地層產(chǎn)生的沉降量不能忽視。

(2)下穿隧道與高速鐵路路基交叉點(diǎn)處為該工程薄弱環(huán)節(jié)，因此，在設(shè)計與施工時要注意對該處進(jìn)行加強(qiáng)。

(3)隧道開挖后，特別是上臺階開挖后，應(yīng)立即施加支護(hù)。

(4)埋深對地層沉降影響比較明顯，隨著地層深度的增加，在動載作用下產(chǎn)生的沉降不斷減小，并且地層動力響應(yīng)從地表到隧道拱頂衰減速度逐漸減小。與日本E rich iTan igueh i的研究結(jié)果一致［12］。

［1］賀斯進(jìn)，藺云宏.地鐵盾構(gòu)區(qū)間穿越鐵路站場設(shè)計與施工分析［J］.現(xiàn)代隧道技術(shù)，2012，49(3):166.

［2］索然緒.雅瀘高速公路徐店子隧道洞口段地震動力響應(yīng)研究［D］.成都:西南交通大學(xué)，2007.

［3］陳衛(wèi)忠，等.ABAQUS在隧道及地下工程中的應(yīng)用［M］.北京:中國水利水電出版社，2010.

［4］傅鑫彬.淺埋軟弱圍巖大跨隧道的施工技術(shù)研究［D］.成都:西南交通大學(xué)，2006.

［5］徐林生，等.洋碰隧道CRD工法施工過程的動態(tài)仿真數(shù)值模擬研究［J］.地質(zhì)災(zāi)害與環(huán)境保護(hù)，2001，12(1):58-62.

［6］朱正國，等.鐵路隧道下穿公路引起的路面沉降規(guī)律和控制基準(zhǔn)研究［J］.巖土力學(xué)，2012，33(2):559-563.

［7］房明，等.新建隧道盾構(gòu)下穿施工對既有隧道影響的三維數(shù)值模擬［J］.鐵道科學(xué)與工程學(xué)報，2011，8(1):68-71.

［8］陳燈紅，等.基于ABAQUS的黏彈性邊界單元及在重力壩抗震分析中的應(yīng)用［J］.世界地震工程，2010，26(3):128-129.

［9］劉晶波，等.三維黏彈性靜-動力統(tǒng)一人工邊界［J］.中國科學(xué)，2005，35(9):967-979.

［10］孫海峰，等.ABAQUS中動力問題邊界條件的選取［J］.地震工程與工程振動，2011，31(3):72-76.

［11］翟婉明.列車-軌道-橋梁動力相互作用理論與工程應(yīng)用［M］.北京:科學(xué)出版社，2011.

［12］李曉霖.地鐵誘發(fā)振動對地面以及地上結(jié)構(gòu)的影響規(guī)律研究［D］.北京:北京工業(yè)大學(xué)，2003.