楊 悅 陳孝國(guó) 劉紀(jì)峰 付曉強(qiáng)

(1.三明學(xué)院建筑工程學(xué)院，福建 三明 365004； 2.黑龍江科技大學(xué)建筑工程學(xué)院，黑龍江 哈爾濱 150022； 3.工程材料與結(jié)構(gòu)加固福建省高等學(xué)校重點(diǎn)實(shí)驗(yàn)室(三明學(xué)院)，福建 三明 365004； 4.三明學(xué)院信息工程學(xué)院，福建 三明 365004)

0 引言

TBM(Tunnel Boring Machine)全斷面掘進(jìn)機(jī)適用于軟巖和硬巖的隧道掘進(jìn)方法，掘進(jìn)、支護(hù)、出渣等施工工序并行連續(xù)作業(yè)。全斷面掘進(jìn)機(jī)是機(jī)、電、液、光、氣等系統(tǒng)集成的工廠(chǎng)化流水線(xiàn)隧道施工裝備，具有掘進(jìn)速度快、利于環(huán)保、綜合效益高等優(yōu)點(diǎn)，可實(shí)現(xiàn)傳統(tǒng)鉆爆法難以實(shí)現(xiàn)的復(fù)雜地理地貌深埋長(zhǎng)隧洞的施工，在中國(guó)鐵道、水電、交通、礦山、市政等隧洞工程中應(yīng)用正在迅猛增長(zhǎng)[1]。

在破碎巖體力學(xué)性質(zhì)方面已經(jīng)有一定的研究成果，金豐年等[2,3]利用數(shù)值模擬方法針對(duì)南京九華山隧道的破碎度研究了圍巖穩(wěn)定性。鄧國(guó)華等[4]對(duì)黃土隧道圍巖的結(jié)構(gòu)變化特征進(jìn)行了詳細(xì)分析，得到了一些有意義的結(jié)論。孫曉勇[5]對(duì)斷層破碎帶中的隧道支護(hù)結(jié)構(gòu)背后注漿層對(duì)管片受力的影響進(jìn)行了研究。楊?lèi)偟萚6-9]利用數(shù)值模擬和理論解析的方法研究了深部隧道圍巖變形和支護(hù)結(jié)構(gòu)力學(xué)性質(zhì)隨著盾構(gòu)機(jī)開(kāi)挖進(jìn)行的演化規(guī)律，得出一些重要結(jié)論，但是沒(méi)有涉及到巖體破碎性的影響。張子新等[10]針對(duì)黏性土地層中盾構(gòu)隧道開(kāi)挖面支護(hù)壓力計(jì)算方法進(jìn)行了探討。關(guān)于TBM(盾構(gòu))施工隧道的力學(xué)性質(zhì)研究成果諸多[11-13]，但是關(guān)于巖體破碎帶的影響的研究?jī)?nèi)容相對(duì)較少，本文擬對(duì)此進(jìn)行研究。

1 巖體破碎度指標(biāo)

1.1 裂隙度K

巖石的裂隙度K是取樣線(xiàn)方向上單位長(zhǎng)度上的節(jié)理數(shù)。對(duì)于一組節(jié)理，設(shè)有一條長(zhǎng)度為l的取樣線(xiàn)，其上有n條節(jié)理出現(xiàn)，則裂隙度為：

(1)

沿著該取樣線(xiàn)，節(jié)理的平均間距d為:

(2)

若巖體上有n組不同方向的節(jié)理，其裂隙度分別為Ka1,Ka2,Kb1,Kb2,…,Kn1,Kn2，取樣線(xiàn)上的節(jié)理平均間距為max,mbx,…,mnx為:

(3)

其中，ξa,ξb,…,ξn分別為各組節(jié)理與取樣線(xiàn)之間的夾角。

該取樣線(xiàn)上的裂隙度K為各組節(jié)理的裂隙度之和，即K=Ka+Kb+…+Kn。

按裂隙度的大小可以將節(jié)理分組：壓縮或糜棱化帶(K=100 m-1～1 000 m-1)；非常密集節(jié)理(K=10 m-1～100 m-1)；密節(jié)理(K=1 m-1～10 m-1)；疏節(jié)理(K=0 m-1～1 m-1)。

1.2 切割度

有些節(jié)理將巖體完全切割開(kāi)，有些節(jié)理切割一部分，節(jié)理在巖體中分離的程度用切割度Xe的大小來(lái)衡量。假設(shè)平直斷面的面積為A，節(jié)理面的面積為a，那么，切割度為:

(4)

若Xe=1，巖體被完全切割開(kāi)；Xe=0，巖體完整無(wú)裂隙；一般情況下，0

(5)

根據(jù)巖體的裂隙度和切割度的關(guān)系，可以將巖體的破碎程度劃分為不同破碎程度的巖體：似均質(zhì)性、弱節(jié)理化、節(jié)理化、強(qiáng)節(jié)理化及完全節(jié)理化。

2 巖體強(qiáng)度理論

巖石的強(qiáng)度理論是巖石在某一應(yīng)力或組合應(yīng)力的作用下，巖石產(chǎn)生的破壞準(zhǔn)則，是在大量巖石力學(xué)試驗(yàn)的基礎(chǔ)上，加以歸納、分析、描述建立起來(lái)的。由于巖石的成分與成因不同，使巖石的破壞特征會(huì)有很多差別。不同的受力狀態(tài)也會(huì)影響巖石的強(qiáng)度。目前巖石強(qiáng)度理論主要有四種：最大剪應(yīng)力理論、八面體剪應(yīng)力強(qiáng)度理論、摩爾—庫(kù)侖強(qiáng)度理論和格里菲斯強(qiáng)度理論。本文研究巖體破碎程度的影響，因此，選用格里菲斯強(qiáng)度理論。

格里菲斯(Griffith)準(zhǔn)則脆性破壞理論認(rèn)為：實(shí)際的固體在結(jié)構(gòu)構(gòu)造上既不是絕對(duì)均勻的，也不是絕對(duì)連續(xù)的，其內(nèi)部包有大量的微裂紋和微孔洞。這種固體在外力作用下，即使作用的平均應(yīng)力不大，但由于微裂紋微孔洞邊緣上的應(yīng)力集中，很可能在邊緣局部產(chǎn)生很大的拉應(yīng)力。當(dāng)這種拉應(yīng)力達(dá)到或超過(guò)抗拉強(qiáng)度時(shí)，微裂紋便開(kāi)始擴(kuò)展，當(dāng)許多這樣的微裂紋擴(kuò)展、聯(lián)合時(shí)，最后使固體沿某一個(gè)或若干個(gè)平面或曲面形成宏觀(guān)破裂。巖石就是這樣一種包含于大量微裂紋和微孔洞的固體材料，因此，格里菲斯理論為巖石破壞判據(jù)提供重要的理論基礎(chǔ)。

3 工程實(shí)例分析

3.1 工程地質(zhì)資料

某TBM施工隧道長(zhǎng)2 453 m，最大埋深為280 m，建設(shè)過(guò)程主要采用TBM施工方法。利用護(hù)盾式掘進(jìn)機(jī)，開(kāi)挖洞徑7.6 m，開(kāi)挖進(jìn)尺為1.2 m。巖體風(fēng)化程度較大，節(jié)理較發(fā)育。地下水資源豐富，開(kāi)挖難度較大。結(jié)合實(shí)際工程地質(zhì)資料和相關(guān)設(shè)計(jì)規(guī)范，取各層巖體的力學(xué)性質(zhì)參數(shù)如表1所示。

表1 巖石力學(xué)性質(zhì)參數(shù)表

3.2 數(shù)值模擬方法

采用FLAC3D有限差分軟件進(jìn)行隧道開(kāi)挖的數(shù)值模擬計(jì)算。由于洞口呈圓形，隧道所穿越的巖層不對(duì)稱(chēng)，因此取整個(gè)洞口及一定范圍內(nèi)的圍巖為研究對(duì)象。開(kāi)挖洞徑7.6 m，根據(jù)數(shù)值模擬中研究范圍的一般規(guī)定，橫向取洞周?chē)鷰r的邊界是距洞中心點(diǎn)5倍左右洞徑，取40 m，縱向沿軸線(xiàn)長(zhǎng)度也取40 m，即取40 m×40 m×40 m的圍巖，中間是直徑為7.6 m的洞口。對(duì)于研究模型上初始荷載取其上覆圍巖的原巖應(yīng)力，應(yīng)力邊界條件是上邊界施加實(shí)際上覆圍巖自重應(yīng)力，利用水土合算的辦法計(jì)算其自重應(yīng)力。隧道有一定的坡度，但是此處不考慮傾斜角度，建模時(shí)認(rèn)為隧道是水平方向。按照常規(guī)建立三維坐標(biāo)系，O點(diǎn)為隧洞圓心，水平向右為OX軸正向，隧洞徑向?yàn)镺Y軸正向，豎直向上為OZ軸正向。

圍巖是有一定破碎程度的非連續(xù)巖體，具有彈塑性變形性質(zhì)，采用FLAC3D有限差分軟件中的Solid三維六面體單元模擬均質(zhì)圍巖，支護(hù)體系采用Solid單元模擬，盾構(gòu)機(jī)采用Shell殼單元模擬。開(kāi)挖過(guò)程采用分步開(kāi)挖，每次開(kāi)挖進(jìn)尺為1.2 m。由于重點(diǎn)分析巖體裂隙度的影響，采用格里菲斯(Griffith)強(qiáng)度理論作為巖體的本構(gòu)模型。先模擬不考慮巖體破碎的情況，再模擬有巖體破碎的實(shí)際情況，將二者計(jì)算結(jié)果進(jìn)行對(duì)比分析。破碎巖體按照均勻化理論處理，將破碎巖體視為均勻連續(xù)材料，只是和非破碎巖石的力學(xué)參數(shù)不同。

根據(jù)實(shí)際的施工參數(shù)進(jìn)行開(kāi)挖模擬，從第一個(gè)開(kāi)挖步開(kāi)始，被開(kāi)挖的巖體賦予空單元來(lái)模擬開(kāi)挖掉。計(jì)算此時(shí)圍巖的徑向位移、縱向位移和最大主應(yīng)力。整個(gè)研究范圍需要32個(gè)開(kāi)挖步才能完成，那么將同樣的開(kāi)挖和計(jì)算過(guò)程循環(huán)31次即可，每次循環(huán)的初始應(yīng)力是上一個(gè)開(kāi)挖步的應(yīng)力值。第32個(gè)開(kāi)挖步完成以后的圍巖的徑向位移量是需要提取的數(shù)據(jù)。

3.3 結(jié)果分析

完成開(kāi)挖以后洞周?chē)鷰r的最終徑向位移云圖如圖1所示。由圖1可見(jiàn)，圍巖最終徑向位移量的分布形態(tài)受巖體破碎的影響較為明顯。相對(duì)于完整巖體，當(dāng)考慮破碎影響時(shí)，圍巖中每一點(diǎn)處的最大徑向位移量均有不同程度的增大。

將第1，10,20,30,32等五個(gè)開(kāi)挖步對(duì)應(yīng)的圍巖位移量列入表2，并繪制于圖2。由表2，圖2可見(jiàn)，各種工況下，隨著開(kāi)挖步的增加圍巖徑向位移量都呈增加趨勢(shì)。裂隙度越大，圍巖徑向位移量越大。

表2 開(kāi)挖支護(hù)過(guò)程中洞周?chē)鷰r的徑向位移 mm

在TBM施工過(guò)程中，隨著開(kāi)挖步的開(kāi)展，圍巖的徑向位移不斷發(fā)展；隨著巖體裂隙度K值的增加，巖體的徑向位移顯著增加；隨著開(kāi)挖步的增加，破碎巖體的位移量受裂隙度的影響越來(lái)越顯著。

4 結(jié)語(yǔ)

首先給出了巖體破碎度的兩個(gè)指標(biāo)，即裂隙度和張開(kāi)度的計(jì)算方法。根據(jù)巖體四大強(qiáng)度理論的區(qū)別，選擇格里菲斯準(zhǔn)則作為巖體破壞的判據(jù)。借助實(shí)際工程案例的地質(zhì)資料，利用FLAC3D有限差分軟件對(duì)TBM施工全過(guò)程進(jìn)行了仿真模擬，得到了不同開(kāi)挖步對(duì)應(yīng)的圍巖徑向位移分布云圖，將不同破碎程度的圍巖位移進(jìn)行了詳細(xì)對(duì)比分析。明確了圍巖的裂隙度對(duì)其位移量的影響規(guī)律，研究成果將為隧洞設(shè)計(jì)和施工技術(shù)提供理論依據(jù)。