黃 阜, 何 瀚

(1. 長(zhǎng)沙理工大學(xué) 土木與建筑學(xué)院, 長(zhǎng)沙 410114; 2. 長(zhǎng)沙理工大學(xué) 橋梁工程安全控制省部共建教育部重點(diǎn)實(shí)驗(yàn)室, 長(zhǎng)沙 410114)

大斷面高速鐵路隧道施工方法的三維數(shù)值研究

黃 阜1,2, 何 瀚1

(1. 長(zhǎng)沙理工大學(xué) 土木與建筑學(xué)院, 長(zhǎng)沙 410114; 2. 長(zhǎng)沙理工大學(xué) 橋梁工程安全控制省部共建教育部重點(diǎn)實(shí)驗(yàn)室, 長(zhǎng)沙 410114)

針對(duì)大斷面高速鐵路隧道施工過(guò)程中出現(xiàn)的復(fù)雜施工力學(xué)問(wèn)題, 以花甲山高速鐵路隧道為依托工程, 采用三維數(shù)值軟件對(duì)花甲山隧道的施工動(dòng)態(tài)過(guò)程進(jìn)行數(shù)值模擬, 得到了隧道圍巖的應(yīng)力、位移分布特征和地表變形規(guī)律, 通過(guò)與現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比, 總結(jié)出大斷面高速鐵路隧道的動(dòng)態(tài)施工力學(xué)特征, 為大斷面高速鐵路隧道施工方案優(yōu)化提供依據(jù).

大斷面鐵路隧道; 數(shù)值模擬; 監(jiān)控量測(cè); 力學(xué)特征; 施工方法

引言

隨著我國(guó)高速鐵路建設(shè)進(jìn)入高峰期, 越來(lái)越多的高速鐵路隧道出現(xiàn)在工程實(shí)踐中, 尤其是在我國(guó)西南部的山地丘陵地帶, 隧道工程更為常見(jiàn). 高速鐵路隧道由于要克服高速列車在通過(guò)隧道時(shí)產(chǎn)生的空氣動(dòng)力學(xué)問(wèn)題, 一般都采用單洞雙線的斷面型式. 這樣的單洞雙線鐵路隧道具有跨度大、斷面扁平的特點(diǎn),而且在爆破過(guò)程中對(duì)圍巖進(jìn)行多次擾動(dòng)以及施工期間諸多工序的相互影響, 導(dǎo)致施工過(guò)程中的動(dòng)態(tài)施工力學(xué)極為復(fù)雜, 易發(fā)生失穩(wěn)或者坍塌事故. 因此, 有必要對(duì)大斷面高速鐵路隧道的施工力學(xué)行為開(kāi)展研究,為其施工方案優(yōu)化提供理論依據(jù).

目前, 廣泛應(yīng)用于大斷面隧道中的施工方法主要有: 全斷面法、CRD 法、雙側(cè)壁導(dǎo)坑法、上下短臺(tái)階法等. 國(guó)內(nèi)已經(jīng)有一些學(xué)者對(duì)大斷面隧道的施工力學(xué)進(jìn)行過(guò)研究, 并取得了一些研究成果[1～3]. 基于數(shù)值模擬技術(shù), 高峰和尹小寶[4]分析了關(guān)子嶺隧道在不同施工方法下圍巖的應(yīng)力和位移分布特征, 為該隧道最優(yōu)施工方案的確定奠定了基礎(chǔ). 劉光唯和丁志亮[5]針對(duì)哈大客運(yùn)專線筆架山隧道的施工特點(diǎn), 優(yōu)化了該隧道施工中采用的交叉中隔墻法, 保障了該隧道施工的順利進(jìn)行. 孫兆遠(yuǎn)等[6]采用數(shù)值模擬方法對(duì)客運(yùn)專線鐵路大斷面隧道的力學(xué)特征進(jìn)行了對(duì)比分析, 總結(jié)了不同施工方法的技術(shù)特點(diǎn), 為大斷面隧道的安全施工提供了參考. 黃嘉億[7]通過(guò)建立三維數(shù)值計(jì)算模型, 綜合分析多種圍巖條件下臺(tái)階法不同參數(shù)對(duì)高速鐵路隧道穩(wěn)定性的影響, 為采用臺(tái)階法施工的高速鐵路隧道參數(shù)選擇提供了理論依據(jù). 模型試驗(yàn)是研究隧道穩(wěn)定性的有效方法之一, 吳夢(mèng)軍和黃倫海[8]利用這一方法研究了四車道公路隧道的施工力學(xué)特征, 并提出了最有利于該隧道圍巖穩(wěn)定的施工方案.

花甲山隧道屬于高速雙線鐵路隧道, 具有跨度大、斷面大的特點(diǎn), 且因地形、地質(zhì)條件復(fù)雜多變, 同一條隧道內(nèi)軟弱圍巖和不良地質(zhì)種類繁多, 造成施工難度大、安全風(fēng)險(xiǎn)高. 為了能使隧道施工在保證安全的前提下滿足工期要求, 利用數(shù)值模擬技術(shù), 對(duì)該隧道的動(dòng)態(tài)施工過(guò)程進(jìn)行了模擬, 得到了隧道圍巖的變形規(guī)律和應(yīng)力分布特征. 通過(guò)對(duì)數(shù)值模擬結(jié)果進(jìn)行分析, 提出有利于花甲山隧道圍巖穩(wěn)定的施工方案, 為工程的順利實(shí)施提供了有效保障.

1 工程概況

花甲山隧道淺埋段圍巖風(fēng)化較為嚴(yán)重, V級(jí)圍巖段原設(shè)計(jì)采用大拱腳臺(tái)階法和CRD法施工. 結(jié)合花甲山地質(zhì)情況, 為了在保證施工安全的前提下加快施工進(jìn)度, 大拱腳臺(tái)階法變更為三臺(tái)階七步施工法. 在圍巖松散破碎或土質(zhì)地段, 為了減少工序轉(zhuǎn)換次數(shù), 同時(shí)保證支護(hù)強(qiáng)度, CRD法在實(shí)際施工中變更為三臺(tái)階預(yù)留核心土法. 三臺(tái)階七步法開(kāi)挖法的具體施工工序如圖1所示.

圖1 V級(jí)圍巖地段三臺(tái)階七步施工工序橫斷面圖

2 數(shù)值模擬計(jì)算過(guò)程

2.1 數(shù)值模型的建立過(guò)程

基于巖體和支護(hù)結(jié)構(gòu)的力學(xué)特征, 隧道圍巖及注漿加固圈都采用Mohr-Coulomb本構(gòu)模型進(jìn)行模擬;砂漿錨桿采用錨桿單元模擬; 鋼拱架和噴射混凝土采用殼單元模擬; 二次襯砌采用實(shí)體單元. 各材料參數(shù)根據(jù)花甲山地質(zhì)勘探報(bào)告以及相應(yīng)的設(shè)計(jì)規(guī)范進(jìn)行取值, 具體數(shù)值如表1～4所示. 隧道模型左右邊界為4倍洞徑, 上下邊界為3.5倍洞徑, 由于實(shí)際工程中采用了三臺(tái)階七步法, 模擬過(guò)程中開(kāi)挖面也分為7個(gè)部分進(jìn)行開(kāi)挖, 每循環(huán)開(kāi)挖進(jìn)尺定為0.6m, 各臺(tái)階之間開(kāi)挖不同步. 圖2顯示了具體計(jì)算模型和開(kāi)挖步驟.

表1 實(shí)體單元力學(xué)參數(shù)

表2 拱墻砂漿錨桿cable結(jié)構(gòu)單元力學(xué)參數(shù)

表3 鋼拱架+網(wǎng)噴初支shell結(jié)構(gòu)單元力學(xué)參數(shù)

表4 臨時(shí)鋼支撐beam結(jié)構(gòu)單元力學(xué)參數(shù)

圖2 隧道模型及開(kāi)挖過(guò)程

2.2 三臺(tái)階法開(kāi)挖數(shù)值模擬結(jié)果

1) 圍巖塑性區(qū)分析

第一輪開(kāi)挖完成后的隧道縱向塑性區(qū)分布, 如圖3所示. 從圖中可以看出, 洞室臺(tái)階面與圍巖交界區(qū)域內(nèi)出現(xiàn)小范圍的剪切破壞區(qū)域; 由于采用小斷面開(kāi)挖, 上部環(huán)形掌子面處于較為穩(wěn)定的狀態(tài), 拱頂位置出現(xiàn)小范圍的受拉破壞; 臺(tái)階核心土受到開(kāi)挖的影響, 出現(xiàn)大范圍的塑性區(qū)域. 完整開(kāi)挖完成后, 隧道塑性區(qū)分布如圖4所示. 從圖中可以看出, 洞門位置出現(xiàn)較大區(qū)域的塑性破壞, 在拱腳位置出現(xiàn)小范圍的剪切塑性區(qū).

從以上分析可知, 隨著隧道臺(tái)階的開(kāi)挖, 隧道拱頂和仰拱處的塑性區(qū)得到了控制. 雖然噴錨支護(hù)阻止了圍巖發(fā)生大面積的屈服破壞, 但是由于噴錨支護(hù)屬于柔性支護(hù), 不能完全限制圍巖的變形, 部分圍巖出現(xiàn)了少量的塑性變形. 因此, 圍巖在噴錨支護(hù)的保護(hù)下, 變形達(dá)到穩(wěn)定后應(yīng)立即施做二次襯砌, 確保圍巖不發(fā)生過(guò)大的變形.

圖3 第一輪開(kāi)挖完成后隧道縱向塑性區(qū)分布

圖4 完全開(kāi)挖后隧道縱向塑性區(qū)分布

2) 圍巖應(yīng)力特征分析

圖5～8為第一輪和最后一輪開(kāi)挖完成后圍巖最大和最小的主應(yīng)力圖. 從這些圖中可以看出, 在隧道開(kāi)挖后, 圍巖應(yīng)力呈對(duì)稱分布, 最大和最小主應(yīng)力在隧道周邊區(qū)域的變化較大, 而在遠(yuǎn)離隧道的區(qū)域變化較小.

圖5 第一輪開(kāi)挖完成后的最大主應(yīng)力圖

圖6 第一輪開(kāi)挖完成后的最小主應(yīng)力圖

圖7 開(kāi)挖完成后的最大主應(yīng)力圖

圖8 開(kāi)挖完成后的最小主應(yīng)力圖

從最大主應(yīng)力等值線圖可以看出, 隧道周邊區(qū)域基本以壓應(yīng)力為主. 此外, 圍巖最大主應(yīng)力在隧道斷面形狀發(fā)生變化的位置變化比較顯著. 在隧道的拱頂和拱底的局部位置出現(xiàn)了拉應(yīng)力. 最小主應(yīng)力等值線的分布特征與最大主應(yīng)力等值線分布特征相似, 在隧道拱腳處出現(xiàn)了應(yīng)力集中現(xiàn)象.

3) 圍巖位移特征分析

由于隧道的開(kāi)挖采用了鉆爆施工, 隧道周邊圍巖被擾動(dòng), 從而導(dǎo)致巖體的初始應(yīng)力場(chǎng)和位移場(chǎng)發(fā)生變化, 為了分析隧道開(kāi)挖對(duì)圍巖初始位移場(chǎng)的影響, 繪制了不同開(kāi)挖狀態(tài)下圍巖位移場(chǎng)的分布圖, 如圖9～12所示.

圖9 第一輪開(kāi)挖完成后圍巖豎向位移分布云圖

圖10 第一輪開(kāi)挖后圍巖水平位移分布云圖

圖11 開(kāi)挖完成后圍巖豎向位移分布圖

圖12 開(kāi)挖完成后圍巖水平位移分布

圖9和圖10為第一輪開(kāi)挖完成后, 隧道周邊圍巖豎向和水平位移分布圖. 從圖中可以看出, 隧道豎向位移從上臺(tái)階掌子面向后發(fā)展, 掌子面處拱頂沉降值為6.88mm, 而拱底隆起值為22.8mm. 由于隧道在施工過(guò)程中采用了超前小導(dǎo)管注漿對(duì)圍巖進(jìn)行加固, 圍巖的最大水平位移為28mm, 出現(xiàn)在拱腰處. 開(kāi)挖完成后, 隧道拱頂?shù)淖畲蟪两禐?2.1mm, 拱底隆起值為25mm. 從以上數(shù)據(jù)可知, 采用三臺(tái)階七步開(kāi)挖方法,能較好地控制隧道周邊圍巖位移的發(fā)展.

4) 隧道初期支護(hù)內(nèi)力分析

隧道的初期支護(hù)由錨桿和噴射混凝土組成, 錨桿和噴射混凝土施做到開(kāi)挖后的圍巖中, 與圍巖形成了共同受力體系. 在錨噴支護(hù)的限制下圍巖能夠產(chǎn)生適量的變形, 從而釋放掉原始巖體中儲(chǔ)存的巨大原始應(yīng)力, 對(duì)維持隧道圍巖的穩(wěn)定性起到重要作用. 隧道施工的進(jìn)行對(duì)初期支護(hù)結(jié)構(gòu)受力影響很大, 通過(guò)分析施工過(guò)程中初期支護(hù)的受力情況, 可以為初期支護(hù)的優(yōu)化設(shè)計(jì)提供有效參考.

圖13～15為隧道完全開(kāi)挖之后, 噴射混凝土及錨桿的受力情況. 從噴射混凝土受力圖中可以看出, 噴射混凝土所受到的剪切應(yīng)力值在隧道縱向分布一致, 而軸力值受臺(tái)階法開(kāi)挖方法的影響在隧道縱向有所變化. 噴射混凝土內(nèi)力的最大值出現(xiàn)在拱腰位置, 因此, 拱腰位置的附近圍巖受到更為明顯的擠壓作用.從圖15可以看出, 錨桿受壓區(qū)和受拉區(qū)基本一致, 受壓區(qū)基本分布在隧道起拱線附近, 其他部分位置的錨桿均受到拉應(yīng)力.

圖13 噴射混凝土軸力圖

圖14 噴射混凝土剪力圖

圖15 錨桿軸力分布情況

3 基于數(shù)值模擬結(jié)果的花甲山隧道施工方案優(yōu)化

現(xiàn)場(chǎng)監(jiān)控量測(cè)數(shù)據(jù), 真實(shí)地反映了實(shí)際施工過(guò)程中圍巖變形的發(fā)展情況, 對(duì)隧道施工安全的監(jiān)控起到了至關(guān)重要的作用. 根據(jù)監(jiān)控量測(cè)采集數(shù)據(jù)結(jié)果并整理出斷面監(jiān)測(cè)點(diǎn)的位移變化情況, 與數(shù)值模擬結(jié)果進(jìn)行對(duì)比, 可以驗(yàn)證數(shù)值模擬結(jié)果的有效性. 圖16為現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)與數(shù)值模擬得到的拱頂沉降曲線對(duì)比圖.

圖16 現(xiàn)場(chǎng)實(shí)測(cè)與數(shù)值模擬拱頂沉降曲線對(duì)比圖

從圖中可以看出, 隨著隧道開(kāi)挖的進(jìn)行, 基于數(shù)值模擬得到的拱頂位移變形曲線與根據(jù)現(xiàn)場(chǎng)監(jiān)控?cái)?shù)據(jù)繪制的拱頂位移變形曲線變化趨于一致. 這說(shuō)明, 基于數(shù)值模擬得到的隧道拱頂位移變化特征與現(xiàn)場(chǎng)實(shí)測(cè)結(jié)果是基本吻合的. 且數(shù)值模擬得到的拱頂位移最終值為32.2mm, 而實(shí)測(cè)得到的拱頂位移最終值為36.6mm, 數(shù)值計(jì)算結(jié)果略小于監(jiān)測(cè)值, 證明了本文的數(shù)值模擬是可靠和有效的.

花甲山隧道監(jiān)測(cè)的最大拱頂沉降量在鐵路隧道施工規(guī)范(TB 10204-2002)的許可范圍之內(nèi). 這說(shuō)明三臺(tái)階七步開(kāi)挖法, 有效限制了圍巖變形的發(fā)展, 并且最大化地發(fā)揮了圍巖的自身承載能力, 為花甲山隧道的安全施工提供了有效保障.

4 結(jié)論

基于數(shù)值模擬技術(shù), 模擬了大斷面高速鐵路隧道的動(dòng)態(tài)施工過(guò)程, 通過(guò)分析圍巖應(yīng)力、位移分布和支護(hù)結(jié)構(gòu)受力等施工力學(xué)特征, 得到了以下結(jié)論:

(1) 三臺(tái)階七步開(kāi)挖法將大斷面隧道分成了7個(gè)部分分步開(kāi)挖, 充分發(fā)揮了圍巖的自身承載能力, 限制了圍巖變形的發(fā)展, 為隧道施工的安全進(jìn)行提供了有效保障;

(2) 為了驗(yàn)證數(shù)值模擬的正確性, 將數(shù)值模擬結(jié)果與現(xiàn)場(chǎng)監(jiān)控量測(cè)數(shù)據(jù)進(jìn)行了對(duì)比. 對(duì)比分析結(jié)果表明, 基于數(shù)值模擬的拱頂位移變形曲線與根據(jù)現(xiàn)場(chǎng)監(jiān)控?cái)?shù)據(jù)繪制的拱頂位移變形曲線基本一致, 證明了本文的數(shù)值模擬是可靠和有效的.

[1] 喬建崗, 關(guān)法志, 高 超, 等. 武廣高鐵獅子石隧道CRD法施工技術(shù)研究[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2010, 7(2): 98~101

[2] 王樂(lè)明, 馬志富, 許占良. 高速鐵路隧道穿越富水全風(fēng)化花崗巖地層施工方法研究[J]. 石家莊鐵道大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 26(zk): 105~114

[3] 孫雪蓮, 金美海. 高速鐵路隧道穿越富水大斷層施工監(jiān)控分析研究[J]. 鐵道建筑, 2008, (12): 43~47

[4] 高 峰, 尹小寶. 蘭渝鐵路關(guān)子嶺隧道施工方法的數(shù)值分析[J]. 現(xiàn)代隧道技術(shù), 2011, 48(2): 37~43

[5] 劉光唯, 丁志亮. 高速鐵路超大斷面隧道CRD施工技術(shù)[J]. 鐵道建筑, 2010, (5): 43~47

[6] 孫兆遠(yuǎn), 羅 瓊, 耿 偉, 等. 鐵路客運(yùn)專線大斷面隧道開(kāi)挖方法選擇[J]. 隧道建設(shè), 2007, 27(zk): 9~15

[7] 黃嘉億. 高速鐵路隧道臺(tái)階法施工技術(shù)參數(shù)分析與應(yīng)用[J]. 現(xiàn)代隧道技術(shù), 2012, 49(3): 77~82

[8] 吳夢(mèng)軍, 黃倫海. 四車道公路隧道動(dòng)態(tài)施工力學(xué)研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2006, 25(zk): 3057~3062

Study of Construction Method for Large Cross Section High Speed Railway Tunnel Using Three-dimensional Numerical Simulation

HUANG Fu1,2, HE Han1
(1.School of Civil Engineering and Architecture, Changsha University of Science and Technology, Changsha 410114, China; 2. Key Laboratory for Safety Control of Bridge Engineering, Ministry of Education and Hunan Province, Changsha University of Science & Technology, Changsha 410114, China)

As the construction mechanics characteristics of large cross section high speed railway tunnel is complex, the dynamic construction process of Hua jia shan tunnel is simulated by using three-dimensional numerical software. The stress and displacement characteristics of rock mass and surface deformation law are obtained. By comparing with field measured data, the dynamic construction mechanical characteristics of large cross section high speed railway tunnel are concluded which provides support for construction scheme optimization of this project.

large cross section railway tunnel; numerical simulation; monitoring measurement; mechanical characteristics, construction method

U45

A

1672-5298(2014)03-0070-06

2014-07-20

國(guó)家自然科學(xué)基金青年科學(xué)基金(51308072); 長(zhǎng)沙理工大學(xué)橋梁工程安全控制省部共建教育部重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(13KC04)

黃 阜(1983 ? ), 男, 湖南岳陽(yáng)人, 博士, 長(zhǎng)沙理工大學(xué)土木與建筑學(xué)院講師. 主要研究方向: 隧道與地下工程