周成濤，陳俊濤

（1.重慶市勘測(cè)院，重慶400020；2.重慶交通大學(xué) 土木建筑學(xué)院，重慶400074）

渝州隧道位于重慶市渝北區(qū)龍王洞背斜東翼，巖層傾角稍陡，線路走向?yàn)?4°～48°，與地質(zhì)構(gòu)造線小角度斜交，屬構(gòu)造剝蝕淺丘地貌，地形坡度一般為10°～20°，局部由于人工開挖形成陡坎，坡度達(dá)60°～70°，上覆土層為殘坡積粉質(zhì)粘土，一般厚為0m～1m，下伏基巖為侏羅系中統(tǒng)沙溪廟組砂巖、砂質(zhì)泥巖。

經(jīng)調(diào)查，隧道現(xiàn)狀分左、右雙洞，均為單向兩車道隧道，平面成喇叭型布置，進(jìn)口段兩隧道結(jié)構(gòu)間凈距約10m，出口段為15m。隧道凈寬10m，凈高6.6m，圓拱直墻，拱厚75cm，邊墻厚110cm。采用礦山法施工，復(fù)合式襯砌，初期支護(hù)以噴射混凝土，錨桿和鋼筋網(wǎng)為主要支護(hù)手段；二次襯砌采用鋼筋混凝土結(jié)構(gòu)，現(xiàn)洞壁穩(wěn)定，無變形跡象，局部有水痕。

根據(jù)設(shè)計(jì)方案，現(xiàn)擬對(duì)原有的渝州隧道進(jìn)行擴(kuò)挖，但不改變?cè)淼雷呦蚣伴L(zhǎng)度。擴(kuò)挖后，隧道洞凈跨17.5m，洞高約8m，雙向四車道，擬采用復(fù)合式襯砌，礦山法施工。

上述工程為典型的小凈距隧道施工問題。在采用鉆爆法進(jìn)行隧道拓寬施工過程，將對(duì)圍巖進(jìn)行二次擾動(dòng)，這是一個(gè)周邊松動(dòng)圍巖發(fā)生進(jìn)一步變形及圍巖應(yīng)力場(chǎng)、應(yīng)變場(chǎng)重新建立的過程［1－5］。一般小凈距隧道施工對(duì)隧道圍巖受力分布及圍巖變形有重要的影響［6－11］。針對(duì)渝州隧道工程特點(diǎn)，本文基于有限元理論方法對(duì)小凈距隧道拓寬前后圍巖應(yīng)力應(yīng)變分析，為類似工程安全施工提供理論支持。

1 隧道有限元模型

模型取隧道受偏壓最嚴(yán)重的最不利斷面進(jìn)行分析；為了減小邊界效應(yīng)，設(shè)置隧道邊緣至模型邊緣為隧道寬度的3倍，隧道低部至模型邊緣設(shè)置為隧道高度的3倍，模型的尺寸為72m×152m［12］。

計(jì)算模型左、右邊界為X方向約束，底部邊界為Y方向約束，頂部邊界為自由面。計(jì)算模型詳見圖1［13］。

另外，在模擬巖土體開挖時(shí)，采用單元的生死狀態(tài)模擬計(jì)算。根據(jù)“殺死”單元部分的應(yīng)力值，計(jì)算得到開挖邊界面各節(jié)點(diǎn)處，由 “死”單元作用的節(jié)點(diǎn)力，通過把與上述節(jié)點(diǎn)力等值反向的力施加到開挖邊界對(duì)應(yīng)的節(jié)點(diǎn)上。“殺死”單元時(shí)，程序用一個(gè)很小的因數(shù)乘以單元的剛度系數(shù)，并從總的質(zhì)量矩陣中消去單元的質(zhì)量來實(shí)現(xiàn)“殺死”單元［14］。

圖1 有限元分析模型

2 計(jì)算結(jié)果分析

隧道圍巖開挖后，改變了之前原始的應(yīng)力狀態(tài)，選取圍巖作為分析對(duì)象，可以更加直觀地反應(yīng)出隧道拓寬前后圍巖的穩(wěn)定情況［15］。

2.1 圍巖位移對(duì)比分析

如圖2、圖3所示分別為原有隧道及隧道拓寬后圍巖所發(fā)生的水平位移云圖。從圖2中我們可以看出，由于受偏壓影響，圍巖的最大水平位移主要發(fā)生在地表上部、隧道左洞右上角：向左方向的水平位移最大值為0.5mm，向右的水平位移最大值約0.4mm；由圖3我們可看到隧道拓寬后圍巖水平位移明顯增大，向左方向最大值達(dá)到1mm。同時(shí)，我們可以發(fā)現(xiàn)左洞隧道拱頂右側(cè)圍巖松動(dòng)圈的范圍明顯擴(kuò)大。

圖2 原有隧道施工完成時(shí)水平位移云圖（m）

圖3 隧道拓寬后水平位移云圖（m）

如圖4、圖5所示分別為原有隧道及隧道拓寬后圍巖所發(fā)生的豎向位移云圖。從圖4可以看出，由于受偏壓影響，隧道右洞的豎向位移大于左洞，并且右洞拱頂圍巖的松動(dòng)范圍也比左洞更大。其中，最大拱頂沉降位移為2.7mm。由圖5可看到隧道拓寬后圍巖豎向位移明顯增大，最大拱頂沉降位移為5.0mm。

2.2 圍巖拉應(yīng)力對(duì)比分析

如圖6、圖7所示分別為原有隧道及隧道拓寬后圍巖拉應(yīng)力云圖。從圖6中可以看出最大拉應(yīng)力主要發(fā)生在右洞拱頂右側(cè)附近，其中最大值為0.33MPa。從圖7中可以看出最大拉應(yīng)力主要發(fā)生在右洞右邊墻附近，其中最大值為0.49MPa。

圖4 原有隧道施工完成時(shí)豎向位移云圖（m）

圖5 隧道拓寬后豎向位移云圖（m）

圖6 原有隧道圍巖拉應(yīng)力分布云圖（Pa）

圖7 隧道拓寬后圍巖拉應(yīng)力分布云圖（Pa）

2.3 圍巖壓應(yīng)力對(duì)比分析

如圖8、圖9所示分別為原有隧道施工完成時(shí)、隧道拓寬后圍巖壓應(yīng)力云圖。從圖5中可以看出最大壓應(yīng)力主要發(fā)生在兩隧道的左右邊墻附近，其中右洞壓應(yīng)力明顯大于左洞，最大值為－2.86MPa。從圖5可以看出圍巖壓應(yīng)力分布情況基本同隧道拓寬以前一樣，但是最大值增大為－3.22MPa。

2.4 圍巖塑性區(qū)對(duì)比分析

如圖10所示為原隧道圍巖塑性區(qū)主要發(fā)生在兩隧道的右邊墻附近，但是塑性變形并不大，為0.21×10－3。從圖11中，我們可以看到塑性應(yīng)變的分布狀況基本不變，但是塑性變形增大到0.39×10－3。

圖8 原有隧道圍巖壓應(yīng)力分布云圖（Pa）

圖9 隧道拓寬后圍巖壓應(yīng)力分布云圖（Pa）

圖10 原有隧道塑性變形云圖

圖11 隧道拓寬后塑性變形云圖

3 結(jié) 論

通過以上對(duì)隧道拓寬前后圍巖位移、應(yīng)力應(yīng)變及塑性區(qū)的分析，可以得出以下結(jié)論：

1）隧道拓寬后，圍巖豎向位移明顯增大；左右隧道拱頂松動(dòng)圈擴(kuò)展也較迅速，甚至可能貫通至地表，建議在拓寬施工過程中及時(shí)對(duì)拱頂進(jìn)行支護(hù)。

2）隧道拓寬后，圍巖的拉應(yīng)力及壓應(yīng)力明顯增大；受上部偏壓影響，其最大值主要是分布在中央巖柱及右洞右邊墻附近。

3）通過對(duì)塑性區(qū)的分析，可以發(fā)現(xiàn)右洞右邊墻塑性應(yīng)變較大，在隧道拓寬施工過程中需要加強(qiáng)支護(hù)力度。

［1］關(guān)寶樹.隧道工程維修管理要點(diǎn)集［M］.北京：人民交通出版社，2004.

［2］汪敏，廖少明，侯學(xué)淵.近距離重疊地鐵隧道盾構(gòu)法施工地面變形預(yù)測(cè)［J］.同濟(jì)大學(xué)學(xué)報(bào)：自然科學(xué)版，2004，32（9）：1141－1146.Wang M，Liao S M，Hou X Y.Ground deformation prediction of shield tunnelling in closely overlapped subway tunnel construction［J］.Journal of Tongji University：Natural Science，2004，32（9）：1141－1146.

［3］畢繼紅，江志峰，常斌.近距離地鐵施工的有限元數(shù)值模擬［J］.巖土力學(xué)，2005，26（2）：277－281.Bi J H，Jiang Z F，Chang B.Numerical simulation for construction of memo tunnels with short distance between them［J］.Rock and Soil Mechanics，2005，26（2）：277－281.

［4］龔建伍.扁平大斷面小凈距公路隧道施工力學(xué)研究［D］.上海：同濟(jì)大學(xué)，2008.

［5］李曉紅，靳曉光，王宏圖.采動(dòng)下臥煤層對(duì)深埋隧道結(jié)構(gòu)穩(wěn)定性的影響［J］.巖土力學(xué)，2005，26（9）：1448－1455.Li X H，Jin X G，Wang H T.Study on structural stability of deep buried highway tunnel with underlying coal bed exploiting［J］.Rock and Soil Mechanics，2005，26（9）：1448－1455.

［6］靳曉光，李曉紅.某深埋長(zhǎng)隧地應(yīng)力演化及圍巖應(yīng)力位移模擬研究［J］.水文地質(zhì)工程地質(zhì)，2004（1）：40－43.Jin X G，Li X H.Numerical modeling of in－situ stress and surrounding rock displacement of deeplong tunnel ［J］.Hydrogeology and Engineering Geology，2004（1）：40－43.

［7］曾小清，曹志遠(yuǎn).半解析數(shù)值法在地鐵工程雙線隧道分析中的應(yīng)用［J］.工程力學(xué)，1998，15（1）：46－52.Zeng X Q，Cao Z Y.Application of semi analytical numerical method to double－tube tunnel structures ［J］.Engineering Mechanics，1998，15（1）：46－52.

［8］張冬梅，黃宏偉，王箭明.盾構(gòu)推進(jìn)引起地面沉降的黏彈性分析［J］.巖土力學(xué)，2001，22（3）：311－314.Zhang D M，Huang H W，Wang J M.The viseo·elastic analysis of the settlement of the soil above the tunnel caused by shield tunneling［J］.Rock and Soil Mechanics，2001，22（3）：3 11－3 14.

［9］于寧，朱合華.盾構(gòu)隧道施工地表變形分析與三維有限元模擬［J］.巖土力學(xué)，2004，25（8）：1331－1334.Yu N，Zhu H H.Analysis of earth deformation caused by shield tunnel construction and 3D－FEM simulation［J］.Rock and Soil Mechanics，2004，25（8）：133l－1334.

［10］劉艷青，鐘世航，盧汝綏，等.小凈距并行隧道力學(xué)狀態(tài)的試驗(yàn)研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2000，19（5）：590－594.Liu Y Q，Zhong S H，Lu R S，et al.Experimental study of mechanical characteristics of twin tunnels with small spacing［J］.Chinese Journal of Rock Mechanics and Engineering，2000，19（5）：590－594.

［11］劉波，陶龍光，丁城剛，等.地鐵雙隧道施工誘發(fā)地表沉降預(yù)測(cè)研究與應(yīng)用［J］.中國礦業(yè)大學(xué)學(xué)報(bào)，2006，35（3）：356－361.Liu B，Tao L G，Ding C G，et al.Prediction for ground subsidence induce by subway tube tunneling ［J］.Journal of China University of Mining＆Technology，2006，35（3）：356－361.

［12］朱漢華.公路隧道設(shè)計(jì)與施工新法［M］.北京：人民交通出版社.2003.

［13］王文正.公路雙連拱隧道開挖方法及施工過程數(shù)值模擬研究［D］.西安：長(zhǎng)安大學(xué)，2003.

［14］龔曙光.有限元分析與ANSYS APDL編程及高級(jí)應(yīng)用［M］.北京：機(jī)械工業(yè)出版社，2009.

［15］徐干成，白洪才，鄭穎人，等.地下工程支護(hù)結(jié)構(gòu)［M］.北京：中國水利水電出版社，2003.