水底明挖隧道深基坑流固耦合數(shù)值分析*

賴(lài)金星 牛方園 樊浩博 邱軍領(lǐng) 李小宏

(長(zhǎng)安大學(xué)公路學(xué)院1) 西安 710064) (中國(guó)有色金屬工業(yè)西安勘察設(shè)計(jì)研究院2) 西安 710054)

水底明挖隧道深基坑流固耦合數(shù)值分析*

賴(lài)金星1)牛方園1)樊浩博1,2)邱軍領(lǐng)1,2)李小宏2)

(長(zhǎng)安大學(xué)公路學(xué)院1)西安 710064) (中國(guó)有色金屬工業(yè)西安勘察設(shè)計(jì)研究院2)西安 710054)

以福州市穿越晉安河水底隧道明挖法深基坑工程為依托，采用Midas/GTS建立基坑流固耦合有限元模型，分析了基坑開(kāi)挖過(guò)程中周?chē)鷳?yīng)力場(chǎng)和滲流場(chǎng)的變化規(guī)律.結(jié)果表明，基坑周邊的地層沉降、坑底隆起變形和圍護(hù)結(jié)構(gòu)的水平位均移隨著開(kāi)挖深度的增加而逐漸增大；基坑支護(hù)樁全部受壓，其軸力隨深度先增大后減小，小于設(shè)計(jì)值，且其最不利位置位于基坑最終開(kāi)挖面深度處；塑性區(qū)隨著開(kāi)挖深度的增加而逐漸下移，發(fā)生塑性破壞的位置位于開(kāi)挖深度6 m處.

水底隧道； 滲流耦合； 數(shù)值模擬；受力變形；基坑

0 引 言

隨著城市基礎(chǔ)建設(shè)的推進(jìn)，城市水底隧道正在迅速發(fā)展，在這些基坑開(kāi)挖過(guò)程中，由于開(kāi)挖面與地表面存在水頭差，促使地下水向坑內(nèi)流動(dòng)，產(chǎn)生巖土變形與其空隙孔間中流體耦合作用問(wèn)題，這種流固耦合作用會(huì)引起滑坡，基坑管涌而引起的坑底隆起等工程事故.因此，在施工之前必須對(duì)滲流引起的周邊地層和支護(hù)結(jié)構(gòu)的影響做好系統(tǒng)、精確的分析，為施工合理規(guī)劃提供理論依據(jù).

國(guó)內(nèi)外學(xué)者通過(guò)數(shù)值分析、現(xiàn)場(chǎng)監(jiān)測(cè)、理論推導(dǎo)等手段對(duì)基坑滲流做了大量的研究，并且取得了豐碩的成果.文獻(xiàn)[1]對(duì)軟土深基坑排水引起的地表沉降進(jìn)行了詳細(xì)研究.姜忻良等[2]通過(guò)不同數(shù)值模擬軟件模擬了基坑開(kāi)挖過(guò)程中周?chē)鷳?yīng)力場(chǎng)及滲流場(chǎng)的變化規(guī)律，并提出了具體的安全性評(píng)價(jià)指標(biāo)，對(duì)深基坑工程的施工具有重要的指導(dǎo)意義.文獻(xiàn)[3]通過(guò)三維數(shù)值模擬，研究了深基坑在降水過(guò)程中地層沉降的規(guī)律.王樹(shù)英等[4]為了分析基坑墊層下支護(hù)的結(jié)構(gòu)應(yīng)力，對(duì)8車(chē)道明挖湖底隧道基坑進(jìn)行了現(xiàn)場(chǎng)測(cè)試.王建秀等[5]對(duì)上海地區(qū)地鐵車(chē)站深基坑降水誘發(fā)沉降機(jī)制及計(jì)算方法進(jìn)行了詳細(xì)的研究，并對(duì)土層排水固結(jié)沉降的公式進(jìn)行了修正.李玉琦等[6]基于非穩(wěn)定滲流理論研究了坑內(nèi)外水頭變化對(duì)作用在圍護(hù)結(jié)構(gòu)上的孔隙水壓力、土壓力,以及側(cè)壓力的影響.

雖然國(guó)內(nèi)外學(xué)者對(duì)基坑滲流做了眾多的研究，但由于不同基坑在開(kāi)挖過(guò)程中面臨的地質(zhì)和水文情況不盡相同，加之水底隧道基坑的滲流情況復(fù)雜，地下水豐富，因此，研究水底明挖隧道基坑的流固耦合問(wèn)題十分必要.鑒于此，以福州市湖東東路水底明挖隧道基坑為工程依托(見(jiàn)圖1)，借助有限元軟件Midas/GTS，建立二維數(shù)值模型，對(duì)基坑開(kāi)挖過(guò)程中的周?chē)翆雍椭ёo(hù)結(jié)構(gòu)的受力的變形規(guī)律進(jìn)行分析，以期為依托工程施工提供一定的理論指導(dǎo).

1 工程概況

福州市湖東東路道路工程道路全長(zhǎng)約為1 900 m，其中水底明挖隧道段約為40 m，隧道呈雙洞雙線(xiàn)布置，東西走向，起點(diǎn)K0+260位于晉安河西側(cè)(見(jiàn)圖1)，根據(jù)地質(zhì)勘測(cè)報(bào)告，本段基坑所處區(qū)域內(nèi)土層自上而下主要有：淤泥、粘土、卵石、砂性粘性土(見(jiàn)圖2)，各種土層物理力學(xué)參數(shù)見(jiàn)表1[7].

圖1 基坑位置

圖2 基坑結(jié)構(gòu)

模型單元E/MPaμγ/(kN·m-3)K/(m·s-1)c/kPaφ/(°)淤泥 90000.3315.41.2×10-71023.6粘土 100.3314.40.00075228.6砂性粘土400.3316.20.0029236.00卵石 2000.2426.52.9×10-890045.00地連墻 200.3022.01.16×10-112941.99536.00鋼支撐 2100000.1879.0立柱 22000.2025.0

2 基坑開(kāi)挖流固耦合數(shù)值模擬

2.1 計(jì)算原理

根據(jù)太沙基原理，總應(yīng)力可以分為有效應(yīng)力和孔隙水壓力.水不能承受剪切應(yīng)力，所以有效剪切應(yīng)力和總的剪切應(yīng)力相等.即總應(yīng)力的表達(dá)式如下.

(1)

(2)

2.2 本構(gòu)模型

模型中各層土體計(jì)算中均采用摩爾-庫(kù)倫彈塑本構(gòu)模型.應(yīng)力空間中摩爾庫(kù)倫本構(gòu)關(guān)系的屈服函數(shù)f(σ1,σ2,σ3)=0表達(dá)式為

(3)

式中：c為粘聚力；Nφ為參數(shù)，由下式確定.

(4)

其中：φ為材料內(nèi)摩擦角.

2.3 初始應(yīng)力

由于模擬的基坑開(kāi)挖為淺地表工程，地表平整，故地基初始應(yīng)力場(chǎng)僅考慮重力場(chǎng)，其計(jì)算式為：

(5)

式中：h為計(jì)算位置距地表的距離；γ為土體的重度，具體大小由下式確定.

(6)

式中：γd為土體的干重度；γw為水的重度；n為土體的孔隙度；Sr為飽和度，γ具體數(shù)值見(jiàn)表1.

2.4 有限元模型和材料參數(shù)

基坑模型寬16 m，開(kāi)挖深度12 m.圍護(hù)結(jié)構(gòu)采用地下連續(xù)墻，其厚度0.6 m，插入坑底以下深度為4 m.支護(hù)樁深入坑底以下2 m.地下連續(xù)墻與周?chē)林g通過(guò)設(shè)置Goodman接觸面單元進(jìn)行連接[8-10].基坑分4次開(kāi)挖，每次開(kāi)挖深度分別為：2，3，3，4 m.設(shè)置鋼支撐4層，支撐中心距地表的深度分別為1,3.5,5.5,9 m.根據(jù)基坑開(kāi)挖影響范圍(水平方向)和影響深度分別約為開(kāi)挖深度的3～4倍和2～4倍，計(jì)算域取為60 m×25 m，而在基坑長(zhǎng)度方向僅取一個(gè)單元進(jìn)行分析(即進(jìn)行二維平面應(yīng)變分析)，具體有限元模型見(jiàn)圖3，基坑各部分材料參數(shù)見(jiàn)表1.

圖3 有限元模型

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

3.1 基坑水頭與孔隙水壓力分布

圖4為水頭等值線(xiàn)云圖，圖5為孔隙水壓力云圖.由圖4、圖5可見(jiàn)，隨著基坑開(kāi)挖深度的增加，坑內(nèi)外的水頭差逐漸增大，最大水頭差達(dá)到13 m，坑內(nèi)的孔隙水壓力值(絕對(duì)值)逐漸減小.在坑內(nèi)，由于水頭差的影響，地下水從坑底下面流向坑底開(kāi)挖面，而開(kāi)挖面上水壓不及滲流帶來(lái)的壓力，所以產(chǎn)生正了負(fù)的孔隙水壓力；而在坑外，地下水從地表向下流向地下，所以產(chǎn)生正的孔隙水壓力.

圖4 水頭等值線(xiàn)

3.2 基坑周邊土層位移

圖6為基坑周邊土層隨開(kāi)挖深度的變形規(guī)律.由圖6可見(jiàn)，隨著開(kāi)挖深度的增加，基坑周邊的地層沉降隨著開(kāi)挖深度的增加也逐漸增大，最大沉降量達(dá)到6 cm且位于坑后約5 m處.基坑開(kāi)挖前要對(duì)周?chē)翆舆M(jìn)行降水施工，土層中的自由水減小引起孔隙水壓力的減小，引起土層的豎直沉降，并且距離基坑越近，降水程度越大，沉降量也越大，但圍護(hù)結(jié)構(gòu)附近仍有部分土層向上隆起，這是因?yàn)樵诘?步開(kāi)挖結(jié)束后，為了防止基坑的水平收斂，在開(kāi)挖面上方施做了鋼支撐，使得圍護(hù)結(jié)構(gòu)附近的地層發(fā)生了向上的隆起.

圖5 孔隙水壓力云圖

圖6 基坑周邊土層隨開(kāi)挖深度的變形規(guī)律

坑底隆起變形隨著開(kāi)挖深度的增加也逐漸增大，并且變形值隨著距基坑中心的距離增加而減小，最終變形值達(dá)到9 cm.由于基坑開(kāi)挖面與地面之間存在很大的水頭差，導(dǎo)致坑外的水通過(guò)坑底流向開(kāi)挖面，并且基坑兩側(cè)有圍護(hù)結(jié)構(gòu)，使得兩側(cè)的隆起量小于基坑中心，但在第3步開(kāi)挖結(jié)束后，靠近圍護(hù)結(jié)構(gòu)處土體可能收到墻體的擠壓比較嚴(yán)重，所以隆起變形較大.

隨著開(kāi)挖深度的增加，圍護(hù)結(jié)構(gòu)的水平位移逐漸增大，并且最大水平位移所在的深度也逐漸下移，最大水平位移從圍護(hù)結(jié)構(gòu)上部的2 m處(第一步開(kāi)挖)逐漸下移到大約6 m處(第4步開(kāi)挖)，并且最大水平位移為3.2 cm.圍護(hù)結(jié)構(gòu)的水平位移曲線(xiàn)呈N形變化，其原因是基坑開(kāi)挖后，兩側(cè)土層失去支撐，在水壓力的作用下向中間靠攏，但由于收到鋼支撐的約束，水平位移在達(dá)到一定程度后不再增大，隨著鋼支撐數(shù)量的增加和作用時(shí)間的持續(xù)，水平位移逐漸減小至零，此時(shí)基坑開(kāi)挖結(jié)束，而開(kāi)挖面下方的圍護(hù)結(jié)構(gòu)在滲流的作用下仍會(huì)發(fā)生水平位移，并且水平位移值隨著深度的增加而增大.

3.3 支護(hù)樁受力特征

圖7為支護(hù)樁受力特征.由圖7a)可見(jiàn)，在基坑開(kāi)挖完成后，支護(hù)樁全部受壓(負(fù)值表示受壓)，軸力隨深度先增大后減小，在基坑最終開(kāi)挖面深度處達(dá)到最大值536.16 kN，然后逐漸減小.圖7b)顯示支護(hù)樁剪力有正有負(fù)，說(shuō)明剪力變化很復(fù)雜，但最大值也發(fā)生在基坑最終開(kāi)挖面深度處，其大小為81.3 kN.根據(jù)支護(hù)樁的尺寸，混凝土的型號(hào)和鋼筋的規(guī)格，計(jì)算出該支護(hù)樁設(shè)計(jì)最大彎矩值為83 kN·m，而模擬結(jié)果顯示，支護(hù)樁的彎矩最大值為41.7 kN·m，小于設(shè)計(jì)值，說(shuō)明該基坑的支護(hù)情況良好.因?yàn)橹ёo(hù)樁受力的最不利位置都位于基坑最終開(kāi)挖面深度處，所以在施工過(guò)程中要加強(qiáng)該處的監(jiān)控量測(cè)，及時(shí)獲得該處的受力情況，為施工安全提供保證.

圖7 支護(hù)樁的受力特征

3.4 基坑塑性應(yīng)變

根據(jù)模型計(jì)算出的塑性區(qū)(見(jiàn)圖8)可知，隨著開(kāi)挖深度的增加，塑性區(qū)所在的深度也逐漸下移，并且塑性區(qū)的大小也隨之增大.土層最初的塑性區(qū)只出現(xiàn)在坑頂上部1 m左右的區(qū)域，在第2次開(kāi)挖結(jié)束后，塑性區(qū)出現(xiàn)在開(kāi)挖面與兩側(cè)墻體交界處，并且按照扇形分布；第3次開(kāi)挖結(jié)束后，塑性區(qū)位置基本沒(méi)變，面積增大，而最終開(kāi)挖結(jié)束后，塑性區(qū)面積與位置均未發(fā)生變化.由此可知，該基坑開(kāi)挖過(guò)程中，發(fā)生塑性破壞的位置位于開(kāi)挖深度6 m處，因此，在施工過(guò)程中應(yīng)該對(duì)該區(qū)域土層加固，防止其發(fā)生塑性破壞.

圖8 基坑塑性區(qū)分布

4 結(jié) 論

1) 隨著基坑開(kāi)挖深度的增加，坑內(nèi)外的水頭差逐漸增大，最大水頭差達(dá)到13 m，坑內(nèi)的孔隙水壓力(絕對(duì)值)也逐漸減小，并且開(kāi)挖面下產(chǎn)生了負(fù)的孔隙水壓力.

2) 隨著開(kāi)挖深度的增加，基坑周邊的地層沉降、坑底隆起變形和圍護(hù)結(jié)構(gòu)的水平位均移隨著開(kāi)挖深度的增加而逐漸增大.

3) 基坑支護(hù)樁全部受壓，其軸力隨深度先增大后減小，最大值為536.16 kN，支護(hù)樁受力的最不利位置位于基坑最終開(kāi)挖面深度處，在施工過(guò)程中要加強(qiáng)該處的監(jiān)控量測(cè).

4) 隨著開(kāi)挖深度的增加，塑性區(qū)所在的深度也逐漸下移，并且塑性區(qū)的大小也隨之增大，最可能發(fā)生塑性破壞的位置位于開(kāi)挖深度6 m處.

[1]ZHOU N，PIETER A V，LOU Rongxiang，et al.Numerical simulation of deep foundation pit dewatering and optimization of controlling land subsidence [J].Engineering Geology，2010，114：251-260.

[2]姜忻良，宗金輝.基坑開(kāi)挖工程中滲流場(chǎng)的三維有限元分析[J].巖土工程學(xué)報(bào)，2006，28(5):565-568.

[3]WANG Jianxiu，F(xiàn)ENG Bo，YU Haipeng，et al.Using partial penetrating wells and curtains to lower the water level of confined aquifer of gravel[J].Engineering Geology，2013，161:16-25.

[4]王樹(shù)英，陽(yáng)軍生，唐 鵬，等.八車(chē)道明挖湖底隧道結(jié)構(gòu)受力現(xiàn)場(chǎng)測(cè)試與安全性評(píng)價(jià)[J].土木工程學(xué)報(bào)，2014，47(11):120-127.

[5]王建秀，吳林高，朱雁飛，等.地鐵車(chē)站深基坑降水誘發(fā)沉降機(jī)制及計(jì)算方法[J].巖石力學(xué)與工程學(xué)報(bào)，2009，28(5):1010-1019.

[6]李玉岐，周 健，謝康和.基坑開(kāi)挖卸載誘發(fā)的滲流分析[J].巖土工程學(xué)報(bào)，2006，28(10):1259-1262.

[7]福州城鄉(xiāng)建設(shè)發(fā)展總公司.湖東東路道路工程設(shè)計(jì)圖[R].上海：華東勘測(cè)設(shè)計(jì)研究院，2010.

[8]安關(guān)峰，張洪彬，劉添俊.旋噴群樁復(fù)合地基承載特性的數(shù)值分析[J].巖土力學(xué),2012,33(3):906-912.

[9]李金湘，蘇文聰，艾桂根.深基坑開(kāi)挖工程引發(fā)的環(huán)境巖土工程問(wèn)題探討[J].西部探礦工程，2007(10)：40-43.

[10]賈 堅(jiān).土體加固技術(shù)在基坑開(kāi)挖工程中的應(yīng)用[J].地下空間與工程學(xué)報(bào)，2007(1)：8-11.

Numerical Investigation of Seepage-stress Interaction During Foundation Pit Excavation of Under-water Tunnel with the Cut and Cover Method

LAI Jinxing1)NIU Fangyuan1)FAN Haobo1,2)QIU Junling1,2)LI Xiaohong2)

(SchoolofHighway,Chang’anUniversity,Xi’an710064,China)1)(Xi’anEngineeringInvestigationandDesignResearchInstituteofChinaNationalNon-FerrousMetalIndustry，Xi’an710054，China)2)

According to the excvation situation of the foundation pit, the seepage-stress interaction is taken into account, a FEM model of the foundation pit is established by Midas/GTS, seepage and the stress field of foundation pit in different working states is simulated by finite element method. The calculate results show that: the vertical subsidence and horizontal displacement of the soil increase, the soil at the bottom of the pit swells; the axial force of the column is compressive stress which first increases with the excavation depth and then decreases, the max press is in the ending excavation surface and less than code requirement; the value and area of plastic zone increase with the excavation of the pit, the location of the plastic failure happens in the excavation depth of 6m.

under-water tunnel; seepage-stress interaction; numerical simulation; mechanical characteristics; foundation pit

2015-01-10

*國(guó)家自然科學(xué)基金項(xiàng)目資助(批準(zhǔn)號(hào)：51378071)

TU436

10.3963/j.issn.2095-3844.2015.02.016

賴(lài)金星(1973- )：男，博士，副教授 ，主要研究領(lǐng)域?yàn)樗淼琅c地下工程