胡金海，范東方

(1.福州機(jī)場(chǎng)二期高速公路有限公司，福建福州 350002;2.同濟(jì)大學(xué)地下建筑與工程系，上海 200092)

連拱隧道圍巖變形和破壞形態(tài)的DDA分析

胡金海1，范東方2

(1.福州機(jī)場(chǎng)二期高速公路有限公司，福建福州 350002;2.同濟(jì)大學(xué)地下建筑與工程系，上海 200092)

連拱隧道圍巖一般為由節(jié)理結(jié)構(gòu)面相互切割的非連續(xù)巖體，利用非連續(xù)性分析方法研究這類圍巖的變形和破壞形態(tài)可以更好地反映工程實(shí)際情況。采用非連續(xù)變形分析方法DDA，對(duì)金雞山連拱隧道圍巖的變形和破壞過(guò)程進(jìn)行了模擬，將整個(gè)變形破壞過(guò)程劃分為了3個(gè)階段:中墻上方巖體變形、地表下沉、滑移面產(chǎn)生階段、隧道左右洞兩側(cè)滑移帶(或滑移面)的形成階段以及中墻頂部塊體失穩(wěn)、隧道上方巖體快速塌落階段。研究了金雞山隧道淺埋圍巖和深埋圍巖的變形破壞特征，其中邊墻部位、靠近中墻的內(nèi)側(cè)拱頂或拱肩部位、外側(cè)拱肩部位會(huì)首先受到變形破壞。

連拱隧道;DDA分析;圍巖變形;節(jié)理;破壞形態(tài)

0 引言

連拱隧道一般屬于淺埋中短長(zhǎng)度隧道，隧道埋深一般不超過(guò)60.0 m［1］，此種條件下的隧道巖體在各種復(fù)雜的地質(zhì)作用下，會(huì)產(chǎn)生永久的變形和構(gòu)造破壞形跡，形成如褶皺、斷層、節(jié)理、層理等結(jié)構(gòu)弱面。這些結(jié)構(gòu)面的存在，不僅削弱了巖體的力學(xué)強(qiáng)度，而且控制著巖體的變形、破壞機(jī)制和力學(xué)法則，了解和掌握連拱隧道在特定節(jié)理展布特征下圍巖潛在變形、破壞形態(tài)對(duì)隧道穩(wěn)定性分析極為重要，是隧道穩(wěn)定分析與穩(wěn)定控制的基礎(chǔ)。

對(duì)于連拱隧道，采用連續(xù)介質(zhì)力學(xué)分析方法研究隧道圍巖變形和破壞形態(tài)難以反映節(jié)理面的力學(xué)特性和其對(duì)工程穩(wěn)定性的影響，分析結(jié)果也常與實(shí)際情況不符。而采用非連續(xù)性分析方法則可以更好地反映工程實(shí)際情況。隧道工程圍巖穩(wěn)定性非連續(xù)性分析方法主要有物理相似模型試驗(yàn)和數(shù)值分析方法。Everling［2］利用相似模型試驗(yàn)對(duì)有支護(hù)結(jié)構(gòu)條件下節(jié)理巖體的變形特性進(jìn)行了探討和分析;Goodman等［3］利用相似模型試驗(yàn)對(duì)節(jié)理巖體中隧道的變形特性進(jìn)行了研究、分析;劉濤、沈明榮［4］按彈性階段相似原則進(jìn)行連拱隧道室內(nèi)模型試驗(yàn)，模擬連拱隧道的施工工況，得出施工過(guò)程中隧道圍巖位移、應(yīng)力和塑性區(qū)分布規(guī)律。

物理相似模型試驗(yàn)雖然能對(duì)特定條件下節(jié)理巖體中隧道圍巖的穩(wěn)定性進(jìn)行反映，但無(wú)法直觀的表現(xiàn)出圍巖變形和破壞的完整過(guò)程，而數(shù)值分析方法則可以模擬出節(jié)理巖體中隧道圍巖從開(kāi)始變形到破壞的整個(gè)過(guò)程。非連續(xù)巖體數(shù)值分析方法有兩種:離散單元法DEM(distinct element method)［5］和非連續(xù)變形分析方法DDA(discontinuous deformation analysis)［6］。本文采用非連續(xù)變形分析方法DDA，對(duì)節(jié)理巖體中連拱隧道圍巖的變形和破壞過(guò)程進(jìn)行了模擬，分析總結(jié)了圍巖變形和破壞形態(tài)特征以及首先受到變形破壞的特征部位。

1 DDA基本理論

DDA以塊體內(nèi)某一點(diǎn)的位移、應(yīng)變?yōu)槲粗?，每個(gè)塊體有6個(gè)自由度，包括3個(gè)剛體運(yùn)動(dòng)項(xiàng)和3個(gè)應(yīng)變項(xiàng)，即(u0，v0，r0，εx，εy，γxy)T，其中u0，v0分別為塊體沿x，y方向的剛體平動(dòng)位移;r0為繞塊體形心(x0，y0)的剛體轉(zhuǎn)動(dòng)角;εx，εy，γxy為點(diǎn)(x0，y0)處的應(yīng)變狀態(tài)。整個(gè)系統(tǒng)根據(jù)最小勢(shì)能原理進(jìn)行求解，塊體系統(tǒng)的總體方程組為:

如果該塊體系統(tǒng)由n個(gè)塊體組成，則有:

式中:Kij(i，j=1，2，…，n)——6×6的子矩陣的元素，Kij與塊體i的材料屬性有關(guān)，Kij(i≠j)與塊體i和j之間的相互作用情況有關(guān);Di——塊體i的位移矢量;Fi——作用于塊體i上的荷載矢量。

2 計(jì)算模型及結(jié)果

利用DDA程序?qū)痣u山隧道中節(jié)理力學(xué)屬性較差的不同埋深隧道圍巖變形、破壞過(guò)程進(jìn)行數(shù)值模擬。

2.1 金雞山隧道概況

金雞山隧道為雙向八車道連拱隧道，設(shè)計(jì)行車速度為100 km/h，限界寬度為17.25 m，凈高5 m。隧道場(chǎng)址區(qū)內(nèi)發(fā)育有2條斷裂構(gòu)造，隧道洞身圍巖受其影響，主要表現(xiàn)為帶內(nèi)圍巖較破碎，節(jié)理裂隙發(fā)育，易形成排水通道，使圍巖級(jí)別變差。

2.2 計(jì)算模型

選取金雞山隧道K22+402(Ⅴ級(jí)圍巖)斷面和K22+562(Ⅳ級(jí)圍巖)斷面為計(jì)算斷面。其中K22 +402斷面主要分布70°∠59°、51°∠175°兩組節(jié)理，且左洞埋深15.0 m，右洞埋深14.0 m;K22+ 562斷面分布67°∠58°、59°∠195°兩組節(jié)理，且左洞埋深31.0 m，右洞埋深28.8 m。兩個(gè)斷面計(jì)算模型尺寸(長(zhǎng)×高)分別為80.0 m×44.0 m、80.0 m× 59.0 m。圍巖參數(shù):重度為20 kN/m3;彈模為1.5 GPa;泊松比為0.28;節(jié)理面(弱化)參數(shù):內(nèi)摩擦角為22°，粘聚力100 kPa，抗拉強(qiáng)度50 kPa。計(jì)算模型如圖1所示。

圖1 金雞山隧道計(jì)算模型

2.3 計(jì)算結(jié)果

金雞山隧道圍巖變形和破壞的主要過(guò)程如下。

(1)隧道拱頂圍巖最先開(kāi)始出現(xiàn)明顯的向下變形，隨即距隧道底部最近的孤立、臨空的塊體出現(xiàn)剪切滑移。

(2)隨著隧道上部圍巖不斷變形，隧道邊墻壁上臨空塊體滑落;在左洞和右洞中的頂部圍巖出現(xiàn)拉裂破壞，拉裂破壞區(qū)域靠近隧道中墻，拉裂破壞形態(tài)呈現(xiàn)正“八”字型。其中，埋深為14.0 m的斷面K22+402比埋深較大的斷面K22+562明顯(圖2a、圖3a)，說(shuō)明隨著埋深的增加，在相同節(jié)理展布特征下，埋深較大的隧道比埋深較淺的隧道較穩(wěn)定。

(3)隧道圍巖塊體變形、破壞繼續(xù)發(fā)展，埋深為14.0 m的K22+402斷面右側(cè)邊墻塊體以沿結(jié)構(gòu)面向下滑移為主，相鄰塊體間的位移差值不明顯。但在距隧道壁一定距離后，滑移結(jié)構(gòu)面上兩相鄰塊體間位移差值增大，出現(xiàn)了裂隙，此種現(xiàn)象可以看成塊體間的碎裂現(xiàn)象;左側(cè)邊墻則以沿節(jié)理面整體向下滑移為主，滑移塊體間仍然緊密接觸(圖2b)。

隨著埋深的增加，邊墻這種變形特征出現(xiàn)了變化，其中，在金雞山隧道右側(cè)邊墻這種向下滑移、碎裂現(xiàn)象在距隧道壁較近的地方(2.0～5.0 m)處已十分嚴(yán)重，而外層圍巖此時(shí)的變形還比較穩(wěn)定，尚未破壞;左側(cè)邊墻臨空塊體出現(xiàn)滑落，但滑移層上的塊體變形不是很明顯(圖3b)。

(4)隨著隧道上方塊體的不斷脫落，埋深較淺的隧道在左右拱肩(或邊墻)處所形成的裂隙不斷發(fā)展，最終形成了兩條滑移帶。

滑移帶的下邊界由邊墻上距隧道底最近的結(jié)構(gòu)面控制，上邊界一般由拱肩部位的結(jié)構(gòu)面控制;兩條滑移帶使隧道上方成倒“八”字型下沉(圖2c)。

而埋深較深時(shí)，隧道左右兩側(cè)滑移性狀出現(xiàn)了變化和轉(zhuǎn)移，在隧道右洞，滑移帶不再是沿著滑移面延伸至地表，而是在距隧道壁一定距離的位置轉(zhuǎn)到與左側(cè)相同方向的(第2組)節(jié)理面，此后沿該節(jié)理面延至地表(圖3c)。

(5)當(dāng)中墻頂部塊體失穩(wěn)后，左右洞拱頂部位的塊體開(kāi)始脫落、下沉。

圖2 K22+402斷面圍巖變形破壞過(guò)程

圖3 K22+562斷面圍巖變形破壞過(guò)程

3 金雞山連拱隧道的變形和破壞形態(tài)特征

3.1 金雞山連拱隧道變形、破壞形態(tài)特征

通過(guò)對(duì)金雞山隧道K22+402(Ⅴ級(jí)圍巖)斷面和K22+562(Ⅳ級(jí)圍巖)斷面圍巖變形破壞過(guò)程的分析，可用圖4、圖5來(lái)說(shuō)明金雞山連拱隧道圍巖變形、破壞的形態(tài)特征。

圖4 淺埋隧道圍巖破壞區(qū)域分布圖

圖5 深埋隧道圍巖破壞區(qū)域分布圖

金雞山隧道圍巖在兩組節(jié)理控制下，隧道圍巖的變形破壞全過(guò)程可分為3個(gè)階段。

第一階段:中墻上方巖體變形、地表下沉、滑移面產(chǎn)生階段，該階段主要以中墻正上方節(jié)理巖體正“八”字型變形、下沉為主。在該階段中，中墻的寬度、強(qiáng)度對(duì)隧道上方巖體變形區(qū)域、變形形態(tài)具有決定性作用。此外，隧道邊墻臨空的塊體出現(xiàn)脫落。

第二階段:隧道左右洞兩側(cè)滑移帶(或滑移面)的形成階段。在超淺埋段(K22+402斷面)左右洞兩側(cè)滑移面分別為兩組不同節(jié)理面，從破壞形態(tài)上呈現(xiàn)為倒“八”字型(如圖4);而埋深較深段(K22+ 562斷面)左右洞滑移面從總體上看為由同一組(節(jié)理間距最小的一組)的兩個(gè)節(jié)理面，即兩個(gè)滑移面的主體相互平行(如圖5)。

第三階段:中墻頂部塊體失穩(wěn)、隧道上方巖體快速塌落階段。隨著隧道兩側(cè)滑移帶的生成，隧道上方巖體不斷下沉，在中墻頂部塊體失穩(wěn)后，隧道頂部圍巖開(kāi)始大面積塌落。

3.2 金雞山連拱隧道變形顯著特征部位

金雞山連拱隧道圍巖變形破壞的3個(gè)階段中，圍巖首先變形破壞的部位受節(jié)理展布特征影響和自身結(jié)構(gòu)特征的影響，就整個(gè)破壞過(guò)程而言，首先受到變形破壞的特征部位有:

(1)邊墻部位。受節(jié)理面切割影響，邊墻部位易形成沿節(jié)理面向臨空面滑移的塊體(節(jié)理面與隧道外徑向方向夾角＜90°)，這些塊體在變形破壞過(guò)程中容易較早出現(xiàn)滑移破壞。

(2)靠近中墻的內(nèi)側(cè)拱頂或拱肩部位。此處部位圍巖受上覆圍巖重力影響較為嚴(yán)重，而隧道中墻由于有限厚度的制約又難以提供足夠支撐力，所以此處圍巖較早的產(chǎn)生變形破壞(即為正“八”字型變形破壞區(qū)域)。

(3)外側(cè)拱肩部位。該部位圍巖處于隧道上覆圍巖松散塌落的邊界處，受到節(jié)理面切割影響和重力影響，此處圍巖也會(huì)較早發(fā)生變形、破壞。

4 結(jié)論

利用非連續(xù)性變形分析方法(DDA)對(duì)金雞山隧道圍巖的變形和破壞過(guò)程進(jìn)行了模擬，對(duì)隧道圍巖的變形和破壞形態(tài)特征進(jìn)行了研究，得出了如下結(jié)論。

(1)金雞山連拱隧道圍巖在兩組節(jié)理的作用下，變形和破壞過(guò)程可以分為3個(gè)階段。第一階段:中墻上方巖體變形、地表下沉、滑移面產(chǎn)生階段;第二階段:隧道左右洞兩側(cè)滑移帶(或滑移面)的形成階段;第三階段:中墻頂部塊體失穩(wěn)、隧道上方巖體快速塌落階段。

(2)第一階段以中墻正上方節(jié)理巖體正“八”字型變形、下沉為主;第二階段中超淺埋圍巖左右洞兩側(cè)滑移面分別為兩組不同節(jié)理面，而埋深較深段左右洞滑移面從總體上看為由同一組(節(jié)理間距最小的一組)的兩個(gè)節(jié)理面，即兩個(gè)滑移面的主體相互平行。

(3)在圍巖的整個(gè)破壞過(guò)程中，首先受到變形破壞的特征部位為:邊墻部位、靠近中墻的內(nèi)側(cè)拱頂或拱肩部位、外側(cè)拱肩部位。

［1］劉洪州，黃倫海.連拱隧道設(shè)計(jì)施工技術(shù)研究現(xiàn)狀［J］.西部探礦工程，2001，(1):54－55.

［2］Everling G.Model study of rock-joint deformation［J］.Int.J.Rock.Mech.Min.Sci＆Geomech.Abstr.，1964，(1):319－326.

［3］Goodman R.E.，Heuze H.E，Bureau G.J.On modeling techniques for the study of tunnels in jointed rock［C］.Fourteenth Symposium on Rock Mechanics:1972.441－479.

［4］劉濤，沈明榮，陶履彬.連拱隧道動(dòng)態(tài)施工模型試驗(yàn)與三維數(shù)值仿真模擬研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2006，25(9):1802－1808.

［5］Cundall，P.A.Acomputermodel for simulating progressive，large scalemovements in blocky rock system［C］.Symposium of International Society of Rock Mechanics.Nancy，F(xiàn)rance:1971.11－18.

［6］Shi，G-H.Discontinuous deformation analysis—a newmodel for the statics and dynamics of block systems［D］.Department of Civil Engineering，University of California，Berkley:1988.

Analysis on Deformation and Failure Mode Characteristics of Surrounding Rock of Multi-arch Tunnel by DDA Method

HU Jin-hai1，F(xiàn)AN Dong-fang2(1.The Second Phase Project of Fuzhou International Airport Highway Corporation，F(xiàn)uzhou Fujian 350002，China;2.Department of Geotechnical Engineering，Tongji University，Shanghai 200092，China)

The surrounding rock of multi-arch tunnel is intermittence jointed rock masses in most cases;it’s more realistic to analyze the deformation and failure mode characteristics of this kind of rock with discontinuous deformation analysis (DDA)method.The deformation and failure process of surrounding rock of Jinjishan multi-arch tunnel was simulated and the whole failure process was divided into 3 stages:rock deformation，surface subsidence and slip plane appearing above the mid-partition;slip planes(or bands)forming on both sides of tunnel;the rock mass instability above the mid-partition and rapid collapsing of rock mass above the tunnel.The paper studied the deformation and failure characters of shallow and deep surrounding rock in Jinjishan tunnel with the conclusion that the sidewall of the tunnel，the dome or the inside spandrel close to mid-partition and the outside spandrel would be first deformed.

multi-arch tunnel;DDA;surrounding rock deformation;joint;failure mode

U45

A

1672－7428(2012)04－0072－04

2011－08－27

胡金海(1966－)，男(漢族)，福建莆田人，福州機(jī)場(chǎng)二期高速公路有限公司董事長(zhǎng)，公路與城市道路專業(yè)，從事工程建設(shè)管理工作，福建省福州市，fz.hjh@163.com。