胡金海，范东方

(1.福州机场二期高速公路有限公司，福建福州 350002;2.同济大学地下建筑与工程系，上海 200092)

连拱隧道围岩变形和破坏形态的DDA分析

胡金海1，范东方2

(1.福州机场二期高速公路有限公司，福建福州 350002;2.同济大学地下建筑与工程系，上海 200092)

连拱隧道围岩一般为由节理结构面相互切割的非连续岩体，利用非连续性分析方法研究这类围岩的变形和破坏形态可以更好地反映工程实际情况。采用非连续变形分析方法DDA，对金鸡山连拱隧道围岩的变形和破坏过程进行了模拟，将整个变形破坏过程划分为了3个阶段:中墙上方岩体变形、地表下沉、滑移面产生阶段、隧道左右洞两侧滑移带(或滑移面)的形成阶段以及中墙顶部块体失稳、隧道上方岩体快速塌落阶段。研究了金鸡山隧道浅埋围岩和深埋围岩的变形破坏特征，其中边墙部位、靠近中墙的内侧拱顶或拱肩部位、外侧拱肩部位会首先受到变形破坏。

连拱隧道;DDA分析;围岩变形;节理;破坏形态

0 引言

连拱隧道一般属于浅埋中短长度隧道，隧道埋深一般不超过60.0 m［1］，此种条件下的隧道岩体在各种复杂的地质作用下，会产生永久的变形和构造破坏形迹，形成如褶皱、断层、节理、层理等结构弱面。这些结构面的存在，不仅削弱了岩体的力学强度，而且控制着岩体的变形、破坏机制和力学法则，了解和掌握连拱隧道在特定节理展布特征下围岩潜在变形、破坏形态对隧道稳定性分析极为重要，是隧道稳定分析与稳定控制的基础。

对于连拱隧道，采用连续介质力学分析方法研究隧道围岩变形和破坏形态难以反映节理面的力学特性和其对工程稳定性的影响，分析结果也常与实际情况不符。而采用非连续性分析方法则可以更好地反映工程实际情况。隧道工程围岩稳定性非连续性分析方法主要有物理相似模型试验和数值分析方法。Everling［2］利用相似模型试验对有支护结构条件下节理岩体的变形特性进行了探讨和分析;Goodman等［3］利用相似模型试验对节理岩体中隧道的变形特性进行了研究、分析;刘涛、沈明荣［4］按弹性阶段相似原则进行连拱隧道室内模型试验，模拟连拱隧道的施工工况，得出施工过程中隧道围岩位移、应力和塑性区分布规律。

物理相似模型试验虽然能对特定条件下节理岩体中隧道围岩的稳定性进行反映，但无法直观的表现出围岩变形和破坏的完整过程，而数值分析方法则可以模拟出节理岩体中隧道围岩从开始变形到破坏的整个过程。非连续岩体数值分析方法有两种:离散单元法DEM(distinct element method)［5］和非连续变形分析方法DDA(discontinuous deformation analysis)［6］。本文采用非连续变形分析方法DDA，对节理岩体中连拱隧道围岩的变形和破坏过程进行了模拟，分析总结了围岩变形和破坏形态特征以及首先受到变形破坏的特征部位。

1 DDA基本理论

DDA以块体内某一点的位移、应变为未知量，每个块体有6个自由度，包括3个刚体运动项和3个应变项，即(u0，v0，r0，εx，εy，γxy)T，其中u0，v0分别为块体沿x，y方向的刚体平动位移;r0为绕块体形心(x0，y0)的刚体转动角;εx，εy，γxy为点(x0，y0)处的应变状态。整个系统根据最小势能原理进行求解，块体系统的总体方程组为:

如果该块体系统由n个块体组成，则有:

式中:Kij(i，j=1，2，…，n)——6×6的子矩阵的元素，Kij与块体i的材料属性有关，Kij(i≠j)与块体i和j之间的相互作用情况有关;Di——块体i的位移矢量;Fi——作用于块体i上的荷载矢量。

2 计算模型及结果

利用DDA程序对金鸡山隧道中节理力学属性较差的不同埋深隧道围岩变形、破坏过程进行数值模拟。

2.1 金鸡山隧道概况

金鸡山隧道为双向八车道连拱隧道，设计行车速度为100 km/h，限界宽度为17.25 m，净高5 m。隧道场址区内发育有2条断裂构造，隧道洞身围岩受其影响，主要表现为带内围岩较破碎，节理裂隙发育，易形成排水通道，使围岩级别变差。

2.2 计算模型

选取金鸡山隧道K22+402(Ⅴ级围岩)断面和K22+562(Ⅳ级围岩)断面为计算断面。其中K22 +402断面主要分布70°∠59°、51°∠175°两组节理，且左洞埋深15.0 m，右洞埋深14.0 m;K22+ 562断面分布67°∠58°、59°∠195°两组节理，且左洞埋深31.0 m，右洞埋深28.8 m。两个断面计算模型尺寸(长×高)分别为80.0 m×44.0 m、80.0 m× 59.0 m。围岩参数:重度为20 kN/m3;弹模为1.5 GPa;泊松比为0.28;节理面(弱化)参数:内摩擦角为22°，粘聚力100 kPa，抗拉强度50 kPa。计算模型如图1所示。

图1 金鸡山隧道计算模型

2.3 计算结果

金鸡山隧道围岩变形和破坏的主要过程如下。

(1)隧道拱顶围岩最先开始出现明显的向下变形，随即距隧道底部最近的孤立、临空的块体出现剪切滑移。

(2)随着隧道上部围岩不断变形，隧道边墙壁上临空块体滑落;在左洞和右洞中的顶部围岩出现拉裂破坏，拉裂破坏区域靠近隧道中墙，拉裂破坏形态呈现正“八”字型。其中，埋深为14.0 m的断面K22+402比埋深较大的断面K22+562明显(图2a、图3a)，说明随着埋深的增加，在相同节理展布特征下，埋深较大的隧道比埋深较浅的隧道较稳定。

(3)隧道围岩块体变形、破坏继续发展，埋深为14.0 m的K22+402断面右侧边墙块体以沿结构面向下滑移为主，相邻块体间的位移差值不明显。但在距隧道壁一定距离后，滑移结构面上两相邻块体间位移差值增大，出现了裂隙，此种现象可以看成块体间的碎裂现象;左侧边墙则以沿节理面整体向下滑移为主，滑移块体间仍然紧密接触(图2b)。

随着埋深的增加，边墙这种变形特征出现了变化，其中，在金鸡山隧道右侧边墙这种向下滑移、碎裂现象在距隧道壁较近的地方(2.0～5.0 m)处已十分严重，而外层围岩此时的变形还比较稳定，尚未破坏;左侧边墙临空块体出现滑落，但滑移层上的块体变形不是很明显(图3b)。

(4)随着隧道上方块体的不断脱落，埋深较浅的隧道在左右拱肩(或边墙)处所形成的裂隙不断发展，最终形成了两条滑移带。

滑移带的下边界由边墙上距隧道底最近的结构面控制，上边界一般由拱肩部位的结构面控制;两条滑移带使隧道上方成倒“八”字型下沉(图2c)。

而埋深较深时，隧道左右两侧滑移性状出现了变化和转移，在隧道右洞，滑移带不再是沿着滑移面延伸至地表，而是在距隧道壁一定距离的位置转到与左侧相同方向的(第2组)节理面，此后沿该节理面延至地表(图3c)。

(5)当中墙顶部块体失稳后，左右洞拱顶部位的块体开始脱落、下沉。

图2 K22+402断面围岩变形破坏过程

图3 K22+562断面围岩变形破坏过程

3 金鸡山连拱隧道的变形和破坏形态特征

3.1 金鸡山连拱隧道变形、破坏形态特征

通过对金鸡山隧道K22+402(Ⅴ级围岩)断面和K22+562(Ⅳ级围岩)断面围岩变形破坏过程的分析，可用图4、图5来说明金鸡山连拱隧道围岩变形、破坏的形态特征。

图4 浅埋隧道围岩破坏区域分布图

图5 深埋隧道围岩破坏区域分布图

金鸡山隧道围岩在两组节理控制下，隧道围岩的变形破坏全过程可分为3个阶段。

第一阶段:中墙上方岩体变形、地表下沉、滑移面产生阶段，该阶段主要以中墙正上方节理岩体正“八”字型变形、下沉为主。在该阶段中，中墙的宽度、强度对隧道上方岩体变形区域、变形形态具有决定性作用。此外，隧道边墙临空的块体出现脱落。

第二阶段:隧道左右洞两侧滑移带(或滑移面)的形成阶段。在超浅埋段(K22+402断面)左右洞两侧滑移面分别为两组不同节理面，从破坏形态上呈现为倒“八”字型(如图4);而埋深较深段(K22+ 562断面)左右洞滑移面从总体上看为由同一组(节理间距最小的一组)的两个节理面，即两个滑移面的主体相互平行(如图5)。

第三阶段:中墙顶部块体失稳、隧道上方岩体快速塌落阶段。随着隧道两侧滑移带的生成，隧道上方岩体不断下沉，在中墙顶部块体失稳后，隧道顶部围岩开始大面积塌落。

3.2 金鸡山连拱隧道变形显著特征部位

金鸡山连拱隧道围岩变形破坏的3个阶段中，围岩首先变形破坏的部位受节理展布特征影响和自身结构特征的影响，就整个破坏过程而言，首先受到变形破坏的特征部位有:

(1)边墙部位。受节理面切割影响，边墙部位易形成沿节理面向临空面滑移的块体(节理面与隧道外径向方向夹角＜90°)，这些块体在变形破坏过程中容易较早出现滑移破坏。

(2)靠近中墙的内侧拱顶或拱肩部位。此处部位围岩受上覆围岩重力影响较为严重，而隧道中墙由于有限厚度的制约又难以提供足够支撑力，所以此处围岩较早的产生变形破坏(即为正“八”字型变形破坏区域)。

(3)外侧拱肩部位。该部位围岩处于隧道上覆围岩松散塌落的边界处，受到节理面切割影响和重力影响，此处围岩也会较早发生变形、破坏。

4 结论

利用非连续性变形分析方法(DDA)对金鸡山隧道围岩的变形和破坏过程进行了模拟，对隧道围岩的变形和破坏形态特征进行了研究，得出了如下结论。

(1)金鸡山连拱隧道围岩在两组节理的作用下，变形和破坏过程可以分为3个阶段。第一阶段:中墙上方岩体变形、地表下沉、滑移面产生阶段;第二阶段:隧道左右洞两侧滑移带(或滑移面)的形成阶段;第三阶段:中墙顶部块体失稳、隧道上方岩体快速塌落阶段。

(2)第一阶段以中墙正上方节理岩体正“八”字型变形、下沉为主;第二阶段中超浅埋围岩左右洞两侧滑移面分别为两组不同节理面，而埋深较深段左右洞滑移面从总体上看为由同一组(节理间距最小的一组)的两个节理面，即两个滑移面的主体相互平行。

(3)在围岩的整个破坏过程中，首先受到变形破坏的特征部位为:边墙部位、靠近中墙的内侧拱顶或拱肩部位、外侧拱肩部位。

［1］刘洪州，黄伦海.连拱隧道设计施工技术研究现状［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］刘涛，沈明荣，陶履彬.连拱隧道动态施工模型试验与三维数值仿真模拟研究［J］.岩石力学与工程学报，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，France: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，FAN Dong-fang2(1.The Second Phase Project of Fuzhou International Airport Highway Corporation，Fuzhou 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－)，男(汉族)，福建莆田人，福州机场二期高速公路有限公司董事长，公路与城市道路专业，从事工程建设管理工作，福建省福州市，fz.hjh@163.com。