謝仁青

福建巖土工程勘察研究院有限公司,福建 龍巖 361000

0 引言

隧道開挖面的穩(wěn)定性關(guān)系到作業(yè)安全,稍有不慎就可能引起塌方。膨脹性巖土主要為親水性的礦物,容易發(fā)生脹縮變化,成為影響開挖面穩(wěn)定性的重要因素,其脹縮機(jī)理較為復(fù)雜,當(dāng)前缺乏可靠的數(shù)學(xué)或物理模型。有限元分析法可用于建立隧道及膨脹土模型,進(jìn)而確定支護(hù)力的設(shè)置要求。國(guó)內(nèi)針對(duì)膨脹性巖土的隧道開挖面穩(wěn)定性研究較為匱乏,故結(jié)合某工程案例對(duì)其開展探究。

1 工程概況

某城市地鐵1號(hào)線設(shè)計(jì)為雙線圓形隧道,采用土壓平衡式盾構(gòu)機(jī)進(jìn)行施工,隧道設(shè)計(jì)外徑為6.0 m,襯砌管片的厚度為0.3 m,因而內(nèi)徑為5.4 m。該項(xiàng)目部分盾構(gòu)施工區(qū)間的地質(zhì)結(jié)構(gòu)較為復(fù)雜,按照自上而下的順序,其土層分別為人工填土(厚度5～13.7 m)、淤泥質(zhì)軟土(厚度小于等于5.4 m),黏土和粉質(zhì)黏土(厚度0.4～3.5 m)、全風(fēng)化泥質(zhì)粉砂巖和泥巖(0.6～20 m)。盾構(gòu)隧道埋深在-12～-30 m,位于全風(fēng)化膨脹巖土層內(nèi),其地質(zhì)結(jié)構(gòu)與隧道的關(guān)系如圖1所示。

圖1 某地鐵盾構(gòu)隧道地質(zhì)結(jié)構(gòu)示意圖

2 盾構(gòu)隧道開挖面穩(wěn)定性分析

在盾構(gòu)隧道開挖面穩(wěn)定性分析中,通常采用理論計(jì)算、現(xiàn)場(chǎng)實(shí)測(cè)以及計(jì)算機(jī)數(shù)值模擬3種方法。理論計(jì)算以極限平衡法和經(jīng)驗(yàn)公式法為主,數(shù)值模擬主要采用有限元分析法?？紤]膨脹巖土的開挖面穩(wěn)定性數(shù)值模擬,本次研究利用ABAQUS軟件對(duì)開挖面進(jìn)行有限元分析,評(píng)價(jià)其穩(wěn)定性。

2.1 盾構(gòu)隧道建模

2.1.1 模型尺寸

隧道建立在x、y、z坐標(biāo)系內(nèi),x方向代表隧道寬度(取50 m),y方向代表隧道的掘進(jìn)方向(取90 m),z方向代表土層厚度(取50 m),模型計(jì)算范圍是隧道直徑的5倍,分別在x、y、z這3個(gè)方向設(shè)置位移約束。模型中的土體材料設(shè)計(jì)為三維六面體單元,襯砌設(shè)計(jì)為殼體結(jié)構(gòu)(厚度為0.3 m)[1]。

2.1.2 設(shè)置模型初始條件

初始條件是指初始地應(yīng)力的平衡條件,可在模型四周設(shè)置邊界條件,并且這些條件應(yīng)與工程實(shí)際情況保持一致,包括土體載荷、土體產(chǎn)生的側(cè)壓力等,初始條件設(shè)置為:

σx=σy=K0σz

(1)

σz=γz′

(2)

式中:σx、σy、σz分別為x、y、z方向的初始地應(yīng)力,z′表示應(yīng)力在豎直方向的深度,將土體的靜止壓力系數(shù)記為K0,土體重度記為γ[2]。

2.1.3 開挖面支護(hù)力表達(dá)

隧道支護(hù)是開挖過(guò)程的必要條件,而支護(hù)力是開挖面穩(wěn)定性分析的重要影響因素。在開挖的初始階段,靜止土體所產(chǎn)生的側(cè)壓力σ0為:

σ0=K0γH

(3)

式中:H表示地面到開挖面中心點(diǎn)的垂直距離,可理解為土層厚度。

支護(hù)應(yīng)力σs可表示為:

σs=λK0γH

(4)

式中:λ為支護(hù)應(yīng)力比。

2.1.4 隧道開挖面失穩(wěn)判斷條件

隧道開挖面失穩(wěn)具有多種判斷方法,常用的判斷條件為位移突變,即支護(hù)應(yīng)力比(σs/σ0)基本穩(wěn)定,但開挖面突然出現(xiàn)較大幅度的水平位移,研究過(guò)程以位移突變?yōu)橹饕袛鄺l件,同時(shí)利用地面沉降、塑性發(fā)展區(qū)云圖等因素進(jìn)行輔助判斷[3]。

2.2 膨脹性巖土理論分析

2.2.1 工程特性指標(biāo)

膨脹性巖土是指遇水后體積膨脹、失水后體積收縮的地質(zhì)巖土,巖土脹縮過(guò)程非常復(fù)雜,當(dāng)前主要采取4個(gè)指標(biāo)對(duì)其進(jìn)行工程描述,分別為膨脹力、自由膨脹率、無(wú)載荷膨脹率以及有載荷膨脹率[4]。

2.2.2 濕度應(yīng)力場(chǎng)

膨脹土受到雨水、地下水的影響,體積容易出現(xiàn)脹縮,成為影響隧道開挖面穩(wěn)定性的重要因素。濕度應(yīng)力場(chǎng)理論可用于描述膨脹性巖土的應(yīng)力變化,為開挖面穩(wěn)定性分析提供理論依據(jù)[5]。國(guó)內(nèi)科研人員在室內(nèi)進(jìn)行膨脹性巖土試驗(yàn),求出濕度線性膨脹系數(shù)α、膨脹土自由膨脹率δef、側(cè)限無(wú)載荷膨脹率δH之間的關(guān)系,具體為:

(5)

式中:Δw為膨脹過(guò)程中的含水率變化;μ是與含水率相關(guān)的參數(shù)。

2.2.3 膨脹性巖土的數(shù)值模型

膨脹性巖土的膨脹機(jī)理非常復(fù)雜,為了便于進(jìn)行數(shù)值建模,需要設(shè)置簡(jiǎn)化的條件,降低建模難度,基本假設(shè)如下。①將含水率的變化區(qū)域限定在開挖面前方的特定距離內(nèi)。②土體線性膨脹系數(shù)穩(wěn)定不變,膨脹率與含水率呈正相關(guān)[6]。③土體泊松比和彈性模量穩(wěn)定不變。

隧道模型的幾何參數(shù)保持不變,在原本的巖土模型中設(shè)置膨脹土的范圍,增加膨脹系數(shù),將膨脹土對(duì)開挖面的影響進(jìn)深設(shè)置為2.0 m。

2.2.4 確定膨脹系數(shù)

巖土的膨脹程度與載荷的大小高度相關(guān),在全風(fēng)化膨脹巖土上方存在人工填土層和黏土層,載荷限制了膨脹土的膨脹率,建模時(shí)認(rèn)為膨脹率與載荷大小呈線性相關(guān)。該地鐵工程膨脹性巖土膨脹力和自由膨脹率統(tǒng)計(jì)結(jié)果如表1所示。根據(jù)表述數(shù)據(jù),膨脹土的自由膨脹率取均值38.7%,膨脹力取均值52.6 kPa。將不同載荷下的側(cè)限膨脹率記為δ,則為:

表1 地鐵項(xiàng)目沿線膨脹性巖土工程特性統(tǒng)計(jì)結(jié)果

(6)

式中:P為膨脹土所受載荷(與支護(hù)力相等)。在式(5)中,將巖土參數(shù)Δw和μ分別設(shè)定為42、0.3,將相關(guān)參數(shù)代入,可計(jì)算出濕度線性膨脹系數(shù)α為0.306。

2.3 數(shù)值模擬過(guò)程及結(jié)果分析

2.3.1 支護(hù)應(yīng)力比與開挖面中心點(diǎn)水平位移模擬結(jié)果

模擬結(jié)果如圖2所示,將開挖面中心點(diǎn)的水平位移變化分為3個(gè)階段。在第一階段,支護(hù)應(yīng)力比介于1.0～0.3,開挖面中心點(diǎn)的水平位移在0～200 mm,位移量較小,并且位移增幅較為緩慢。在第二階段,支護(hù)應(yīng)力比介于0.3～0.2,中心點(diǎn)水平位移量在200～310 mm,水平位移增幅加快,臨近失穩(wěn)狀態(tài)。在第三階段,支護(hù)應(yīng)力比介于0～0.2,開挖面中心點(diǎn)水平位移量最大可達(dá)到760 mm,位移增幅快速升高,進(jìn)入失穩(wěn)狀態(tài)。

圖2 開挖面中心點(diǎn)水平位移量與支護(hù)應(yīng)力比的關(guān)系

另外,圓點(diǎn)線條部分的含水率變化值為0%,其含義為不考慮膨脹因素,當(dāng)含水率變化值為10%時(shí),意味著膨脹土吸收水分,發(fā)生膨脹作用。從模擬結(jié)果來(lái)看,膨脹因素加劇了開挖面失穩(wěn),中心點(diǎn)水平位移量明顯增加。根據(jù)以上分析,支護(hù)應(yīng)力比不可低于0.2。

2.3.2 隧道開挖面支護(hù)力與地面沉降的關(guān)系

隧道開挖面失穩(wěn)常伴隨地面沉降,在模型地表設(shè)置一系列觀測(cè)點(diǎn),分別沿橫向、縱向進(jìn)行布置,模擬并觀察支護(hù)力變化與地面沉降量之間的關(guān)系。以坐標(biāo)系原點(diǎn)O為參考點(diǎn),沿y方向布置縱向觀測(cè)點(diǎn),沿x方向布置橫向觀測(cè)點(diǎn)。將支護(hù)應(yīng)力比分別設(shè)置為0、0.1、0.2、0.3、0.5、0.8,觀察橫縱向沉降量,以縱向觀測(cè)點(diǎn)為例,其沉降量模擬結(jié)果如表2。沉降量中的負(fù)值表示地面下沉,正值表示地面隆起。該結(jié)果表明,地表沉降幅度與支護(hù)應(yīng)力比呈負(fù)相關(guān),提高支護(hù)應(yīng)力比有助于提高開挖面的穩(wěn)定性。

表2 縱向觀測(cè)點(diǎn)沉降量模擬結(jié)果示例

2.3.3 含水率變化對(duì)開挖面穩(wěn)定性的影響

將膨脹性巖土的含水率分別設(shè)定為0%、5%、10%、15%,支護(hù)應(yīng)力比從1.0梯度下降至0(單次下降0.1),觀察開挖面中心點(diǎn)的水平位移量。模擬結(jié)果顯示,當(dāng)支護(hù)應(yīng)力比在1.0～0.3時(shí),4種含水率對(duì)應(yīng)的中心點(diǎn)水平位移曲線基本重合。當(dāng)支護(hù)應(yīng)力比介于0.3～0.2時(shí),水平位移量出現(xiàn)分叉,與含水率的增加呈正相關(guān)。當(dāng)支護(hù)應(yīng)力比介于0.2～0時(shí),水平位移量達(dá)到最大值。0%含水率對(duì)應(yīng)的最大位移量為650 mm,5%含水率對(duì)應(yīng)的最大位移量為700 mm,10%含水率對(duì)應(yīng)的最大位移量為760 mm,15%含水率對(duì)應(yīng)的最大位移量為830 mm,不同含水率下開挖面的極限支護(hù)力如表3所示。

表3 不同含水率下開挖面極限支護(hù)力

3 工程應(yīng)用效果

將地鐵盾構(gòu)隧道的實(shí)際工況參數(shù)設(shè)置為模型的邊界條件(見(jiàn)表4),同時(shí)考慮該項(xiàng)目的全風(fēng)化膨脹巖土,自由膨脹率取38.7%,膨脹力取101 kPa。按照上文所建立的盾構(gòu)隧道模型開展開挖面穩(wěn)定性分析,將支護(hù)應(yīng)力比設(shè)置在0～1.0。模擬結(jié)果顯示,當(dāng)支護(hù)應(yīng)力比在1.0～0.5時(shí),中心點(diǎn)位移量為0～100 mm,當(dāng)支護(hù)應(yīng)力比介于0.5～0時(shí),中心點(diǎn)位移量為100～920 mm。判斷失穩(wěn)的條件為開挖面中心點(diǎn)位移超過(guò)100 mm,故該項(xiàng)目支護(hù)應(yīng)力比不應(yīng)低于0.5,否則將產(chǎn)生開挖面失穩(wěn)的風(fēng)險(xiǎn)。

表4 地鐵項(xiàng)目材料參數(shù)

4 結(jié)束語(yǔ)

膨脹性巖土對(duì)隧道開挖面的穩(wěn)定性具有較大的危害,由于此類地質(zhì)條件容易出現(xiàn)脹縮,有可能引起較大的應(yīng)力變化,因此必須設(shè)計(jì)足夠的開挖面支護(hù)力。本文建立了隧道開挖面和膨脹土模型,在不同支護(hù)應(yīng)力比條件下模擬開挖面水平位移和地表沉降。經(jīng)過(guò)對(duì)比,發(fā)現(xiàn)含水率越高時(shí),膨脹土引起的位移量就越大。由此可見(jiàn),作業(yè)時(shí)應(yīng)該設(shè)置足夠的支護(hù)應(yīng)力比。