謝仁青
福建巖土工程勘察研究院有限公司,福建 龍巖 361000
隧道開挖面的穩(wěn)定性關(guān)系到作業(yè)安全,稍有不慎就可能引起塌方。膨脹性巖土主要為親水性的礦物,容易發(fā)生脹縮變化,成為影響開挖面穩(wěn)定性的重要因素,其脹縮機理較為復(fù)雜,當(dāng)前缺乏可靠的數(shù)學(xué)或物理模型。有限元分析法可用于建立隧道及膨脹土模型,進(jìn)而確定支護(hù)力的設(shè)置要求。國內(nèi)針對膨脹性巖土的隧道開挖面穩(wěn)定性研究較為匱乏,故結(jié)合某工程案例對其開展探究。
某城市地鐵1號線設(shè)計為雙線圓形隧道,采用土壓平衡式盾構(gòu)機進(jìn)行施工,隧道設(shè)計外徑為6.0 m,襯砌管片的厚度為0.3 m,因而內(nèi)徑為5.4 m。該項目部分盾構(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)示意圖
在盾構(gòu)隧道開挖面穩(wěn)定性分析中,通常采用理論計算、現(xiàn)場實測以及計算機數(shù)值模擬3種方法。理論計算以極限平衡法和經(jīng)驗公式法為主,數(shù)值模擬主要采用有限元分析法。考慮膨脹巖土的開挖面穩(wěn)定性數(shù)值模擬,本次研究利用ABAQUS軟件對開挖面進(jìn)行有限元分析,評價其穩(wěn)定性。
2.1.1 模型尺寸
隧道建立在x、y、z坐標(biāo)系內(nèi),x方向代表隧道寬度(取50 m),y方向代表隧道的掘進(jìn)方向(取90 m),z方向代表土層厚度(取50 m),模型計算范圍是隧道直徑的5倍,分別在x、y、z這3個方向設(shè)置位移約束。模型中的土體材料設(shè)計為三維六面體單元,襯砌設(shè)計為殼體結(jié)構(gòu)(厚度為0.3 m)[1]。
2.1.2 設(shè)置模型初始條件
初始條件是指初始地應(yīng)力的平衡條件,可在模型四周設(shè)置邊界條件,并且這些條件應(yīng)與工程實際情況保持一致,包括土體載荷、土體產(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ù)是開挖過程的必要條件,而支護(hù)力是開挖面穩(wěn)定性分析的重要影響因素。在開挖的初始階段,靜止土體所產(chǎn)生的側(cè)壓力σ0為:
σ0=K0γH
(3)
式中:H表示地面到開挖面中心點的垂直距離,可理解為土層厚度。
支護(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)較大幅度的水平位移,研究過程以位移突變?yōu)橹饕袛鄺l件,同時利用地面沉降、塑性發(fā)展區(qū)云圖等因素進(jìn)行輔助判斷[3]。
2.2.1 工程特性指標(biāo)
膨脹性巖土是指遇水后體積膨脹、失水后體積收縮的地質(zhì)巖土,巖土脹縮過程非常復(fù)雜,當(dāng)前主要采取4個指標(biāo)對其進(jìn)行工程描述,分別為膨脹力、自由膨脹率、無載荷膨脹率以及有載荷膨脹率[4]。
2.2.2 濕度應(yīng)力場
膨脹土受到雨水、地下水的影響,體積容易出現(xiàn)脹縮,成為影響隧道開挖面穩(wěn)定性的重要因素。濕度應(yīng)力場理論可用于描述膨脹性巖土的應(yīng)力變化,為開挖面穩(wěn)定性分析提供理論依據(jù)[5]。國內(nèi)科研人員在室內(nèi)進(jìn)行膨脹性巖土試驗,求出濕度線性膨脹系數(shù)α、膨脹土自由膨脹率δef、側(cè)限無載荷膨脹率δH之間的關(guān)系,具體為:
(5)
式中:Δw為膨脹過程中的含水率變化;μ是與含水率相關(guān)的參數(shù)。
2.2.3 膨脹性巖土的數(shù)值模型
膨脹性巖土的膨脹機理非常復(fù)雜,為了便于進(jìn)行數(shù)值建模,需要設(shè)置簡化的條件,降低建模難度,基本假設(shè)如下。①將含水率的變化區(qū)域限定在開挖面前方的特定距離內(nèi)。②土體線性膨脹系數(shù)穩(wěn)定不變,膨脹率與含水率呈正相關(guān)[6]。③土體泊松比和彈性模量穩(wěn)定不變。
隧道模型的幾何參數(shù)保持不變,在原本的巖土模型中設(shè)置膨脹土的范圍,增加膨脹系數(shù),將膨脹土對開挖面的影響進(jìn)深設(shè)置為2.0 m。
2.2.4 確定膨脹系數(shù)
巖土的膨脹程度與載荷的大小高度相關(guān),在全風(fēng)化膨脹巖土上方存在人工填土層和黏土層,載荷限制了膨脹土的膨脹率,建模時認(rèn)為膨脹率與載荷大小呈線性相關(guān)。該地鐵工程膨脹性巖土膨脹力和自由膨脹率統(tǒng)計結(jié)果如表1所示。根據(jù)表述數(shù)據(jù),膨脹土的自由膨脹率取均值38.7%,膨脹力取均值52.6 kPa。將不同載荷下的側(cè)限膨脹率記為δ,則為:
表1 地鐵項目沿線膨脹性巖土工程特性統(tǒng)計結(jié)果
(6)
式中:P為膨脹土所受載荷(與支護(hù)力相等)。在式(5)中,將巖土參數(shù)Δw和μ分別設(shè)定為42、0.3,將相關(guān)參數(shù)代入,可計算出濕度線性膨脹系數(shù)α為0.306。
2.3.1 支護(hù)應(yīng)力比與開挖面中心點水平位移模擬結(jié)果
模擬結(jié)果如圖2所示,將開挖面中心點的水平位移變化分為3個階段。在第一階段,支護(hù)應(yīng)力比介于1.0~0.3,開挖面中心點的水平位移在0~200 mm,位移量較小,并且位移增幅較為緩慢。在第二階段,支護(hù)應(yīng)力比介于0.3~0.2,中心點水平位移量在200~310 mm,水平位移增幅加快,臨近失穩(wěn)狀態(tài)。在第三階段,支護(hù)應(yīng)力比介于0~0.2,開挖面中心點水平位移量最大可達(dá)到760 mm,位移增幅快速升高,進(jìn)入失穩(wěn)狀態(tài)。
圖2 開挖面中心點水平位移量與支護(hù)應(yīng)力比的關(guān)系
另外,圓點線條部分的含水率變化值為0%,其含義為不考慮膨脹因素,當(dāng)含水率變化值為10%時,意味著膨脹土吸收水分,發(fā)生膨脹作用。從模擬結(jié)果來看,膨脹因素加劇了開挖面失穩(wěn),中心點水平位移量明顯增加。根據(jù)以上分析,支護(hù)應(yīng)力比不可低于0.2。
2.3.2 隧道開挖面支護(hù)力與地面沉降的關(guān)系
隧道開挖面失穩(wěn)常伴隨地面沉降,在模型地表設(shè)置一系列觀測點,分別沿橫向、縱向進(jìn)行布置,模擬并觀察支護(hù)力變化與地面沉降量之間的關(guān)系。以坐標(biāo)系原點O為參考點,沿y方向布置縱向觀測點,沿x方向布置橫向觀測點。將支護(hù)應(yīng)力比分別設(shè)置為0、0.1、0.2、0.3、0.5、0.8,觀察橫縱向沉降量,以縱向觀測點為例,其沉降量模擬結(jié)果如表2。沉降量中的負(fù)值表示地面下沉,正值表示地面隆起。該結(jié)果表明,地表沉降幅度與支護(hù)應(yīng)力比呈負(fù)相關(guān),提高支護(hù)應(yīng)力比有助于提高開挖面的穩(wěn)定性。
表2 縱向觀測點沉降量模擬結(jié)果示例
2.3.3 含水率變化對開挖面穩(wěn)定性的影響
將膨脹性巖土的含水率分別設(shè)定為0%、5%、10%、15%,支護(hù)應(yīng)力比從1.0梯度下降至0(單次下降0.1),觀察開挖面中心點的水平位移量。模擬結(jié)果顯示,當(dāng)支護(hù)應(yīng)力比在1.0~0.3時,4種含水率對應(yīng)的中心點水平位移曲線基本重合。當(dāng)支護(hù)應(yīng)力比介于0.3~0.2時,水平位移量出現(xiàn)分叉,與含水率的增加呈正相關(guān)。當(dāng)支護(hù)應(yīng)力比介于0.2~0時,水平位移量達(dá)到最大值。0%含水率對應(yīng)的最大位移量為650 mm,5%含水率對應(yīng)的最大位移量為700 mm,10%含水率對應(yīng)的最大位移量為760 mm,15%含水率對應(yīng)的最大位移量為830 mm,不同含水率下開挖面的極限支護(hù)力如表3所示。
表3 不同含水率下開挖面極限支護(hù)力
將地鐵盾構(gòu)隧道的實際工況參數(shù)設(shè)置為模型的邊界條件(見表4),同時考慮該項目的全風(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時,中心點位移量為0~100 mm,當(dāng)支護(hù)應(yīng)力比介于0.5~0時,中心點位移量為100~920 mm。判斷失穩(wěn)的條件為開挖面中心點位移超過100 mm,故該項目支護(hù)應(yīng)力比不應(yīng)低于0.5,否則將產(chǎn)生開挖面失穩(wěn)的風(fēng)險。
表4 地鐵項目材料參數(shù)
膨脹性巖土對隧道開挖面的穩(wěn)定性具有較大的危害,由于此類地質(zhì)條件容易出現(xiàn)脹縮,有可能引起較大的應(yīng)力變化,因此必須設(shè)計足夠的開挖面支護(hù)力。本文建立了隧道開挖面和膨脹土模型,在不同支護(hù)應(yīng)力比條件下模擬開挖面水平位移和地表沉降。經(jīng)過對比,發(fā)現(xiàn)含水率越高時,膨脹土引起的位移量就越大。由此可見,作業(yè)時應(yīng)該設(shè)置足夠的支護(hù)應(yīng)力比。