• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

正断型断裂模拟及其对山岭隧道影响研究

汪振, 钟紫蓝, 赵密, 杜修力, 黄景琦

汪振, 钟紫蓝, 赵密, 杜修力, 黄景琦. 正断型断裂模拟及其对山岭隧道影响研究[J]. 岩土工程学报, 2020, 42(10): 1876-1884. DOI: 10.11779/CJGE202010013
引用本文: 汪振, 钟紫蓝, 赵密, 杜修力, 黄景琦. 正断型断裂模拟及其对山岭隧道影响研究[J]. 岩土工程学报, 2020, 42(10): 1876-1884. DOI: 10.11779/CJGE202010013
WANG Zhen, ZHONG Zi-lan, ZHAO Mi, DU Xiu-li, HUANG Jing-qi. Simulation of normal fault rupture and its impact on mountain tunnels[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(10): 1876-1884. DOI: 10.11779/CJGE202010013
Citation: WANG Zhen, ZHONG Zi-lan, ZHAO Mi, DU Xiu-li, HUANG Jing-qi. Simulation of normal fault rupture and its impact on mountain tunnels[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(10): 1876-1884. DOI: 10.11779/CJGE202010013

正断型断裂模拟及其对山岭隧道影响研究  English Version

基金项目: 

国家自然科学基金项目 51678015

广东省地震工程与应用技术重点实验室开放基金项目 2017B030314068

详细信息
    作者简介:

    汪振(1993—),男,博士研究生,主要从事隧道等地下结构性能研究。E-mail:wangzhenSKG@126.com

    通讯作者:

    赵密, E-mail:zhaomi@bjut.edu.cn

  • 中图分类号: U45

Simulation of normal fault rupture and its impact on mountain tunnels

  • 摘要: 活动断层的运动和与其相交山岭隧道的震害密切相关,是隧道抗震设计所面临的严峻挑战之一。为此,以滇中引水工程香炉山隧洞为研究背景,运用断裂力学中的黏性界面模型结合有限元方法模拟正断层破裂过程。通过震害调查和试验结果对三维自由场的响应进行验证。进一步建立跨断层隧道三维数值分析模型探讨了不同断层错动量和倾角对隧道响应的影响规律,并引入损伤指数建立隧道安全评价的分类标准对结构的震害进行初步评估。结果表明:正断层错动所引起的地表破裂机制主要有弯曲陡坎和断裂陡坎;隧道衬砌的轴向拉应变和环向剪应变在其与断层滑动面相交位置处达到最大;断层错动量和倾角的变化对隧道不同震害状态沿纵向分布有明显影响;随断层倾角的增大,隧道衬砌处于严重损伤和完全损毁状态的长度要明显减小,断层倾角在50°~70°时对结构安全更为不利。
    Abstract: The seismic damage of mountain tunnels is closely associated with the movement of active faults. Seismic design of tunnels crossing active faults is one of the great challenges nowadays. Based on the engineering prototype of the Xianglu mountain tunnel, the water diversion project in central Yunnan Province, a numerical method to simulate the propagation of normal fault rupture is proposed using the finite element method incorporated with the cohesive interface model in fracture mechanics. The proposed method is verified against the post-earthquake reconnaissance and experimental results using the three-dimensional free-field model. It is used to simulate a tunnel crossing a normal fault, and the effects of fault displacement and dip angle on the response of the tunnel linings are discussed. Besides, the damage indices and safety assessment criteria are introduced to preliminarily evaluate the damage of the tunnel linings subjected to fault movement. The results show that the mechanisms of surface rupture exhibit the forms of folding or fault scarps under normal faulting. The axial tensile strain and hoop shear strain of the tunnel linings reach the maximum at the position where they intersect the fault slip surface. The seismic damage state of tunnel along the longitudinal direction is significantly affected by the fault displacement and dip angle. The length of the tunnel linings in a severely damaged and completely damaged state is significantly reduced with the increase of the dip angle. Dip angles of 50° to 70° are more detrimental to structural safety.
  • 土石混合体(soil-rock mixture,S-RM)是一种极不均匀的松散材料,如崩积层、残积层和洪积层,主要组成成分是碎块石和土[1]。在中国西南山区的基础设施建设过程中,土石混合体被广泛地用作土工材料[2]填充到低洼场地,形成了众多的高填方体边坡。高填方体边坡滑坡破坏制约着中国西南山区基础设施建设的发展,如攀枝花机场填方体滑坡[3]、贵州某高填方体边坡滑坡[4],其原因之一为土石混合料填筑体与下伏基岩界面间的抗滑阻力不足[3]。鉴于此,为了提供更大抗滑阻力,基岩坡面常采用台阶式开挖(如图1所示),形成了长宽比不同的台阶形状基岩界面。尽管如此,对大多数高填方体边坡来说,其填方体与下伏基岩接触面仍然是边坡的第一不连续面,受地震、堆积体重力、降雨作用等因素的影响,此接触面常常是高填方体边坡失稳不容忽视的潜在滑移面之一。

    图  1  填方体边坡示意图
    Figure  1.  Schematic diagram of fill slope

    受施工场地地形条件约束,不同地形下基岩放坡开挖所采用的坡率亦有所差异,从宏观上可认为不同坡率开挖所形成的台阶状基岩界面其粗糙度不同。大量研究表明接触面粗糙度对土-基岩界面的剪切强度及变形特征有重要影响。Zhang等[5-6]对颗粒土-钢板接触面进行的大型接触面静力和循环加载试验结果表明,接触面强度符合莫尔-库仑强度准则;随着钢板表面粗糙度的增加,接触面摩擦角呈现增大的趋势。Borana等[7-8]对全风化花岗岩土壤-钢板接触面进行的直剪试验也表明,钢板表面粗糙度越大,接触面抗剪强度越高。并指出法向应力变化历史对接触面变形特征和剪切破坏强度参数影响较大。石熊等[9]研究揭示了红黏土与混凝土接触面粗糙度对接触面抗剪及残余强度有较大的影响。陈俊桦等[10]指出红黏土与混凝土接触面剪切破坏、变形与接触面粗糙度有较大关系,接触面内摩擦角受其影响较小,接触面黏聚力随接触粗糙度增大而增大并趋近红黏土黏聚力。张吉顺等[11]重点研究了粗糙度对黏性土与混凝土结构接触面的影响规律,试验制作了表面为规则的锯齿凹凸混凝土试块,指出粗糙度明显影响其力学特性。以上学者在研究接触面粗糙度对颗粒材料和结构物接触面的剪切特性影响机制中,常常采用砂土、黏土和粗粒土,几乎很少考虑土石混合体(填方体)这种极不均匀、松散的材料。虽然陈静等[12]通过土石混体与混凝土桩接触面室内大型直剪研究了含水率和含石率对其接触面力学特性的影响;Cen等[13]对土石混合体-台阶状界面接触面开展了数值模拟剪切试验,研究了含石率和土石混体最大块石粒径对其接触界面的力学特性的影响。但接触面粗糙度对土石混合体-基岩接触面剪切力学特性的影响尚不多见,仍需深入研究。鉴于此,本文通过土石混合体-基岩接触面室内大型直剪试验,较系统地研究了接触面粗糙度对填方体-下伏基岩界面的剪切强度及变形特征的影响。可为高填方体或堆积体边坡稳定性分析提供参考。

    室内大型直剪试验采用四川大学华西岩土研究所研制的大型粗粒土压缩直剪仪ZJ50-2G(见图2),仪器主要由刚性框架、上下剪切盒、水平加载装置、垂直加载装置和数据采集装置等组成,直剪试验过程中,垂直方向加载到法向压力设计值,固定上剪切盒并水平推动下剪切盒,使试样均匀受剪。仪器采用数字控制系统,可以自动化采集数据。试验所用上、下剪切盒长×宽×高为300 mm×300 mm×200 mm。

    图  2  ZJ50-2G大型粗粒土压缩直剪仪示意图
    Figure  2.  ZJ50-2G large coarse soil compression direct shear apparatus

    土石混合料试样取自重庆武隆某机场填方区具有代表性的强夯区域,为粉质黏土、黏土以及灰岩碎块石组成的混合物,试样基本物理参数指标如表1所示。

    表  1  土石混合体及灰岩基本物理参数指标
    Table  1.  Basic physical parameter indexes of soil-rock aggregate
    土体类型物理参数指标
    干密度/(g·m-3)孔隙比天然含水率/%天然密度/(kg·m-3)c/kPaφ/(°)弹性模量/GPa单轴抗压强度/MPa
    土石混合体17880.249.32211023.910.54
    灰岩2730143335.6729.1468.09
    下载: 导出CSV 
    | 显示表格

    由于室内试验剪切盒尺寸限制,将现场原土石混合体级配按相同比例缩小,根据《水电水利工程粗粒土试验规程》[14]用相似级配法可将原样粒径按照几何相似条件等比例缩小,各粒组的相对比例保持不变。

    将场地原土石级配采用相似级配法等比例缩小后,试验级配中最大块石粒径为25 mm(如图3)。试验土石混合体级配选择5 mm作为土石阈值[15],试验级配曲线见图4所示。

    图  3  样本筛分结果
    Figure  3.  Sample sieving results
    图  4  土石混合体试样级配曲线
    Figure  4.  Gradation curve of soil-rock mixture samples

    试验所用基岩取自重庆武隆仙女山机场现场,地层为二叠系(P)上统吴家坪组地层,岩性主要是泥质灰岩。采用WDAJ-600型微机控制电液伺服多功能试验机测得岩石试件的单轴抗压强度、抗剪强度指标和弹性模量如表1所示。

    台阶状灰岩试件设计坡率分别为1∶2,1∶1.75和1∶1.5。台阶高度h不变,台阶长度L发生变化,使其满足设计的坡率,制作了表面粗糙度不同的3组试件,每组试件长×宽×高为300 mm×300 mm×200 mm。参考张吉顺等[11]、赵春风等[16]对不同混凝土试件粗糙度的定义,本试验所用灰岩试件粗糙度定义为(如图5):

    Y=R/S, (1)
    图  5  台阶状灰岩试件界面示意图
    Figure  5.  Interface diagram of step limestone test piece

    式中,R为台阶峰谷距离,S为两峰的距离。

    为更好地描述土石混合体-基岩接触面的剪切机理,定义台阶J斜面角度α

    α=arctan(h/L), (2)

    式中,h为台阶高度,L为台阶长度(如图5)。

    根据粗糙度与J斜面的定义,3组试件的粗糙度及台阶斜面J的角度如表2所示。

    表  2  台阶基岩界面力学参数指标
    Table  2.  Mechanical parameter indexes of step bedrock interface
    试件编号坡率台阶高/cm台阶宽/cm粗糙度Y/mmJ斜面α/(°)
    11∶225500.39926.56
    21∶1.752543.750.43129.74
    31∶1.52537.440.46233.69
    下载: 导出CSV 
    | 显示表格

    将重庆武隆仙女山机场现场不规则的灰岩试块,按照表2设计方案加工后的试件如图6所示。

    图  6  加工后的岩石试件
    Figure  6.  Processed rock samples

    将台阶状灰岩试件放入下剪切盒中,在上剪切盒中放入土石混合体。室内试验所用土石混合体试验级配如图4所示,将筛分好的不同粒径块石与土颗粒先拌和至粗细颗粒分布均匀,按天然含水率(9.32%)分3次加入水,搅拌至试样干湿均匀;称取适量拌和均匀的土样分3次逐层均匀装填入上剪切盒内,每层都采用静压法压实[17],按压实系数0.9控制,并对各层接触面凿毛处理,防止层面过于明显,影响剪切效果。考虑到在剪切过程中由于台阶与块石之间的咬合作用,部分块石很难翻越相邻块石或台阶,随着剪切位移的增加,这些块石可能会被剪碎。为了探究剪切面上块石的破碎现象,本文在土样装填完成后,用红色油漆将剪切面块石染色,在试验结束后统计染色块石破碎数量,以此来判断剪切面上块石的破碎情况。

    夏红春等[18]建议直剪试验剪切速率采用0.02~1.2 mm/min,本文采用0.8 mm/min的剪切速率控制。如表3 所示,每组试件在法向应力为 200,400,600和800 kPa下进行土石混合体-基岩接触面的剪切试验。根据《水电水利工程粗粒土试验规程》[14]当剪切位移为60 mm(试样长度的20%)时终止试验。

    表  3  室内大型直剪试验方案
    Table  3.  Indoor large-scale direct shear test schemes
    试件编号粗糙度C法向压力/kPa
    10.399200,400,600,800
    20.431200,400,600,800
    30.462200,400,600,800
    下载: 导出CSV 
    | 显示表格

    不同法向应力水平作用下,土石混合体与不同粗糙度试块接触面的剪应力-剪切位移曲线如图7所示。

    图  7  剪切应力-剪切位移曲线
    Figure  7.  Shear stress-shear displacement curves

    图7可知,剪应力-剪切位移曲线在低法向应力(如图7(a),(b))作用下曲线前期呈现出应变硬化现象,后期呈现出塑性应变现象。曲线有明显的峰值,并且接触面粗糙度越大,峰值剪应力对应的位移越小。而在高法向应力(如图7(c),(d))作用下,曲线均呈现出应变硬化现象,当剪切位移在35~45 mm之前时,剪应力随剪切位移增长较快,之后剪应力仍持续增长但增长速度开始变缓,整个剪切过程无明显峰值;相同法向应力水平作用下,接触面粗糙度越大,剪应力-剪切位移曲线的斜率越大,土石混合体-基岩接触面剪切刚度[19]越大。

    在剪切过程中,随着剪切位移不断增加,台阶凹槽内块石与接触面附近一定范围内块石不断错动和旋转,凹槽内块石与接触面上部一定范围内土石混合体中的块石相互咬合、闭锁形成了块石骨架(如图8),块石骨架又会与台阶产生咬合力,随着剪切位移的增加,咬合作用力逐渐上升,接触面抗剪强度也在逐渐提高。当剪切继续进行,凹槽内块石、上部土石混合体中块石与台阶三者之间的咬合作用逐渐增强到某一值时,块石骨架中部分块石被剪碎,接触面上剪应力发生跌落。在低法向应力作用下,剪应力跌落后剪切面发生塑性流动破坏(如图7(a),(b));而在高法向应力作用下,块石被剪碎,剪切面上的块石数量增加,并能重新发生错动、咬合,形成新的更多的咬合块石骨架(如图9),因此剪切应力-剪切位移曲线表现出剪应力短暂跌落后又上升的现象(如图7(c),(d))。

    图  8  块石骨架形成示意图
    Figure  8.  Schematic diagram of block stone skeleton formation
    图  9  块石破碎后重新形成骨架示意图
    Figure  9.  Schematic diagram of re-formed skeleton after rock is broken

    由于台阶状灰岩的弹性模量远大于土石混合体的弹性模量几个数量级,可将灰岩试件近似看作刚性体[10],因此接触面剪切区的切向变形主要由灰岩试件剪切面上部土石混体和灰岩试件台阶凹槽内土石混合体产生。凹槽内的土石混体沿着台阶斜面发生爬坡现象,将凹槽内的土石混体沿斜面的变形分解为竖向变形分量和切向变形分量,粗糙度大时,剪切斜面J的倾角α越大,在产生相同切向变形下,将产生更大的竖向变形,需要更大的克服法向应力做功,因此剪切应力增长更快,在低法向应力作用下,能更快地达到峰值。

    接触面附近的块石破碎量反映了土石混合填方体中块石与基岩接触面相互咬合的紧密程度,咬合作用越强烈,则在剪切时块石越容易被剪碎。如图11所示,块石破碎形态可分为完全破碎、部分破碎和表面磨损3种。块石发生完全破碎或部分破碎意味着块石所受的咬合作用较大,不易翻越相邻块石或台阶,因此块石容易被剪碎(如图10(a))。块石发生表面磨损说明此块石虽然与其它块石或台阶产生了接触,但是咬合作用很小,在剪切过程中主要以滑移为主,所以只在块石表面留下了摩擦的痕迹(如图10(b))。

    图  11  块石破碎形态
    Figure  11.  Broken forms of block stone
    图  10  块石破碎模式示意图
    Figure  10.  Schematic diagram of block crushing mode

    染色的块石是随机布设在剪切面上的,因此在一定程度上染色块石的破碎情况可反映整个剪切面破碎情况。试验结束后挑出事先随机分布于剪切面上的染色块石进行观察,按照图11对块石破碎形态的分类方法,将不同粗糙度下破碎的染色块石进行分类统计,结果如图12所示。由图12可知,相同法向应力作用下,接触面粗糙度越大,完全破碎块石数量越多,完全破碎和部分破碎块石总数也越多。

    图  12  染色块石破碎与粗糙度关系
    Figure  12.  Relationship between broken stone blocks and roughness

    由于试件表面粗糙度较大时,台阶斜面J也越陡,台阶凹槽内土石混体在沿着台阶斜面J发生爬坡运动时,斜面越陡其法向位移分量越大,剪切破坏区厚度增加,剪切破坏区内有更多的块石相互接触,形成强度较高的块石骨架,随着剪切位移不断增加,紧密咬合的块石更容易被台阶剪断,所以试件粗糙度增加时,块石破碎的现象会加剧。完全破碎或者部分的块石数量可以从侧面体现出剪切面上块石与块石、块石与台阶咬合作用和嵌固作用的强弱程度,在剪切过程中咬合作用和嵌固作用力越强,块石越不易翻越相邻块石或台阶,则块石容易被剪碎。

    根据《水电水利工程粗粒土试验规程》[14],抗剪强度取剪应力与水平位移关系曲线上峰值或稳定值,如无明显峰值,则取水平位移达到试样直径或长度10%处的剪应力作为抗剪强度。统计不同粗糙度接触面在不同法向压力下的抗剪强度如表4所示,从表中可见,在相同的法向应力作用下,随着接触面粗糙度的增加,接触面抗剪强度显著提高。

    表  4  界面抗剪强度
    Table  4.  Interface shear strengths
    项目试件编号
    123
    抗剪强度σn=200 kPa179.0500224.2300264.8700
    σn=400 kPa356.9700401.5600460.8900
    σn=600 kPa490.3700536.6700575.6700
    σn=800 kPa605.5600694.8200794.5600
    粗糙度0.39900.43100.4620
    相关系数R0.98990.99750.9874
    表观黏聚力/kPa54.750077.600098.0400
    内摩擦角/(°)35.240037.720040.4300
    下载: 导出CSV 
    | 显示表格

    在剪切过程中,接触面上的抗剪强度主要来自于两部分:一是克服块石与块石之间、块石与台阶之间的咬合作用力,台阶粗糙度越大,块石与块石之间、块石与台阶之间的咬合作用力越强,在整体上则表现为接触面的抗剪强度越大;二是克服块石与块石、块石与基岩结构面之间的滑移阻力,台阶凹槽内的土石混合体会沿着台阶斜面J发生爬坡现象,爬坡现象增加了剪切破坏区内土石混合体在法向方向上的挤压作用,粗糙度越高,凹槽内土石混合体法向位移分量越大,法向挤压作用越明显,滑移阻力越大,因此需要更大的克服法向应力做功,从而接触面抗剪强度得以提高。同时块石的破碎以及破碎后块石的重新排列需要外力额外做功[20],这些额外做的功也将会使接触面抗剪强度增大。

    不同接触面粗糙度下抗剪强度-法向应力关系曲线如图13所示。进行拟合后其线性相关系数R均大于0.98,这表明在本文试验的法向应力范围内,接触面的抗剪强度与法向应力关系具有较好的线性相关。参考已有研究[10-13]中研究接触面抗剪强度参数的计算方法,本文采用莫尔-库仑准则对不同粗糙度接触面的抗剪强度与法向应力关系进行描述。

    图  13  不同接触面粗糙度下抗剪强度与法向应力关系曲线
    Figure  13.  Relationship between shear strength and normal stress under different contact surface roughnesses

    采用莫尔-库仑强度准则对试验数据进行拟合,得到其表观黏聚力[20-21]和内摩擦角,表观黏聚力和内摩擦角与粗糙度关系曲线如图14所示。由图14可知,随着粗糙度的增大,内摩擦角和表观黏聚力呈线性增长趋势。粗糙度由0.399增加到0.462时,接触面的内摩擦角和表观黏聚力分别提高了14.72%和79.07%,表明粗糙度对表观黏聚力的影响更为显著。

    图  14  接触面强度参数与粗糙度关系曲线
    Figure  14.  Relationship between contact surface strength parameters and roughness

    在剪切过程中沿接触面发生剪切破坏主要克服摩擦阻力和表观黏聚力。在本文试验中摩擦阻力可分为两部分,一是台阶凹槽内土石混合体与剪切面上部土石混合体之间的滑动摩擦阻力,二是台阶凹槽内土石混合体与台阶凹槽表面的滑动摩擦阻力。由于试验所用土石混合体级配、含水率均保持一致,不同试件的台阶凹槽表面几乎无任何差异,因此不同试件接触面上的摩擦阻力也应该差别不大,但是随着剪切位移的增加,剪切面上的块石发生错动和破碎使这种理想的状态发了改变。李广信[20]认为颗粒的破碎和重排列均需要外力额外做功,这将会增加剪切面上的内摩擦角。随着粗糙度的增加剪切面上块石破碎越严重,正是这种现象导致了接触面的内摩擦角随着粗糙度的增加有轻微的上升。

    沿接触面发生剪切变形需要克服的另一种因素是表观黏聚力。马林[21]认为表观黏聚力并非来源于颗粒间的胶结和各种化学键,表观黏聚力表现为剪切面上不规则土颗粒之间的咬合力。在本试验中,接触面附近的块石与块石、块石与台阶相互咬合。当剪切不断进行时,剪切面上的块石要么翻越相邻块石或者台阶,要么被剪断,发生剪切破坏。参考前人研究成果[21],本文认为表观黏聚力产生的原因是块石与块石之间、块石与台阶之间的咬合作用力。随着粗糙度的增加,块石与块石之间、块石与台阶之间的咬合作用力随之增强,导致接触面的表观黏聚力显著提高。需要补充的一点是,接触面上不仅有表观黏聚力,还有台阶凹槽内土石混合体与剪切面上部土石混合体中土颗粒间的黏聚力,但由于本试验土石混合体含石率较大,作为胶结材料的土颗粒之间的接触较少,所以这部分黏聚力很小,接触面上主要由表观黏聚力占主导作用。

    Cen等[13]认为剪切带是接触面附近变形较大的区域,而变形较小的区域被称为非剪切带。在本试验中,凹槽内的土石混体沿着台阶斜面发生爬坡现象,将凹槽内的土石混体沿斜面的变形分解为竖向变形分量和切向变形分量,粗糙度大时,剪切斜面J的倾角α越大,在产生相同切向变形下,将产生更大的竖向变形;同时接触面粗糙度较大时,凹槽内土石混合体法向位移分量越大,法向挤压作用越明显,这使得凹槽内块石与接触面上部一定范围内块石相互咬合、闭锁形成的块石骨架范围更广,当块石骨架在剪切过程中被破坏时,变形较大的区域随之增加。因此接触面粗糙度较大时,剪切过程中形成的剪切带较宽,如图15所示。

    图  15  不同粗糙度下剪切带示意图
    Figure  15.  Schematic diagram of shear zone with different roughnesses

    通过对土石混体—基岩接触面室内大型直剪试验,分析了基岩在不同粗糙度下接触面的剪切变形特征。

    (1)在低法向应力作用下,剪应力-剪切位移曲线前期呈现出应变硬化现象,后期呈现出塑性应变现象,且接触面粗糙度越大接触面发生剪切破坏时变形越小;在高法向应力作用下,剪应力-剪切位移曲线呈现出应变硬化现象,无明显峰值。相同法向应力水平作用下,接触面粗糙度越大,剪应力-剪切位移曲线的斜率越大,土石混合体-基岩接触面剪切刚度越大。

    (2)剪切界面上块石的破碎形态可分为完全破碎、部分破碎和表面磨损3种。相同法向应力作用下,接触面粗糙度越大,完全破碎块石数量越多,完全破碎和部分破碎块石总数也越多。说明在剪切过程中剪切面上块石与块石、块石与台阶咬合作用和嵌固作用力越强,块石越不易翻越相邻块石或台阶,则块石容易被剪碎。

    (3)采用莫尔-库仑强度准则对试验数据进行了线性拟合,获得了其表观黏聚力和内摩擦角。结果表明:在相同法向应力作用下,接触面抗剪强度随着接触面粗糙度增加而增大;土石混体-基岩接触面的内摩擦角和表观黏聚力随着接触面粗糙度的增加也都有所增大,其中接触面的表观黏聚力增大较为明显。

    (4)接触面粗糙度对剪切带宽度有影响作用,表现为接触面粗糙度越大,剪切带越宽。

  • 图  1   香炉山隧洞引水工程地质剖面图

    Figure  1.   Geological profile of Xianglu mountain tunnel

    图  2   钻爆法施工隧洞支护结构图

    Figure  2.   Schematic of tunnel supporting system constructed by drilling and blasting method

    图  3   黏性界面单元COH3D8空间示意图

    Figure  3.   Spatial representation of a three-dimensional cohesive interface element COH3D8

    图  4   黏性界面单元双线性本构模型

    Figure  4.   Bilinear constitutive model of cohesive interface element

    图  5   三维自由场数值分析模型

    Figure  5.   Three-dimensional numerical model of free-field condition

    图  6   正断层错动下自由场工况的等效塑性应变云图

    Figure  6.   Equivalent plastic strain contours of free-field condition under normal fault rupture

    图  7   汶川地震中“正断层”性质的地表破裂图

    Figure  7.   Surface rupture with features of “normal fault” in Wenchuan earthquake

    图  8   文献[5]的有限元分析与试验结果对比

    Figure  8.   Comparison between finite element analysis and test results of Reference [5]

    图  9   断层—隧道三维数值分析模型

    Figure  9.   Three-dimensional numerical model for a tunnel crossing a fault

    图  10   不同错动量下隧道直径变形率

    Figure  10.   Deformation ratios of tunnel diameter at different fault displacements

    图  11   不同错动量下隧道轴向应变云图

    Figure  11.   Axial strain contours of tunnel at different fault displacements

    图  12   不同错动量下隧道环向剪应变云图

    Figure  12.   Hoop shear strain contours of tunnel at different fault displacements

    图  13   不同错动量下OLDT变化规律

    Figure  13.   Distribution of OLDT at different fault displacements

    图  14   不同倾角下隧道直径变形率

    Figure  14.   Deformation ratios of tunnel diameter at different dip angles

    图  15   不同倾角下隧道轴向应变云图

    Figure  15.   Axial strain contours of tunnel at different dip angles

    图  16   不同倾角下隧道环向剪应变云图

    Figure  16.   Hoop shear strain contours of tunnel at different dip angles

    图  17   不同倾角下OLDT变化规律

    Figure  17.   Distribution of OLDT at different dip angles

    表  1   鹤庆—洱源断裂岩石的力学参数

    Table  1   Geomechanical parameters of Heqing-Eryuan fault

    类型密度ρ/(103kg·m-3)弹性模量E/GPa剪切模量G/GPa泊松比内摩擦角φ/(°)黏聚力c/MPa抗拉强度ft/MPa抗剪强度τ/MPa断裂能G/(N·m-1)
    完整岩石2.97.50.28451.1
    断层破碎带2.11.50.560.33290.150.10.320
    下载: 导出CSV

    表  2   混凝土损伤塑性模型力学参数

    Table  2   Mechanical parameters of concrete damaged plasticity model

    参数数值
    密度ρ/(kg·m-3)2.5×103
    弹性模量E/GPa30
    泊松比0.2
    剪胀角ψ/(°)36.31
    压缩屈服应力fc/MPa20.1
    拉伸屈服应力ft/MPa2.01
    下载: 导出CSV
  • [1] 洪开荣. 中国隧道及地下工程近两年的发展与展望[J]. 隧道建设, 2017, 37(2): 123-134. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201702002.htm

    HONG Kai-rong. Development and prospects of tunnels and underground works in china in recent two years[J]. Tunnel Construction, 2017, 37(2): 123-134. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201702002.htm

    [2] 钮新强, 张传健. 复杂地质条件下跨流域调水超长深埋隧洞建设需研究的关键技术问题[J]. 隧道建设, 2019, 39(4): 523-536. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201904001.htm

    NIU Xin-qiang, ZHANG Chuan-jian. Some key technical issues on construction of ultra-long deep-buried water conveyance tunnel under complex geological conditions[J]. Tunnel Construction, 2019, 39(4): 523-536. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201904001.htm

    [3] 何川, 李林, 张景, 等. 隧道穿越断层破碎带震害机理研究[J]. 岩土工程学报, 2014, 36(3): 427-434. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201403006.htm

    HE Chuan, LI Lin, ZHANG Jing, et al. Seismic damage mechanism of tunnels through fault zones[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(3): 427-434. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201403006.htm

    [4] 崔光耀, 王明年, 于丽, 等. 汶川地震断层破碎带段隧道结构震害分析及震害机理研究[J]. 土木工程学报, 2013, 46(11): 122-127. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201311019.htm

    CUI Guang-yao, WANG Ming-nian, YU Li, et al. Study on the characteristics and mechanism of seismic damage for tunnel structures on fault rupture zone in Wenchuan seismic disastrous area[J]. China Civil Engineering Journal, 2013, 46(11): 122-127. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201311019.htm

    [5]

    KIANI M, AKHLAGHI T, GHALANDARZADEH A. Experimental modeling of segmental shallow tunnels in alluvial affected by normal faults[J]. Tunnelling and Underground Space Technology, 2016, 51: 108-119. doi: 10.1016/j.tust.2015.10.005

    [6]

    CAI Q P, PENG J M, NG C W W, et al. Centrifuge and numerical modelling of tunnel intersected by normal fault rupture in sand[J]. Computers and Geotechnics, 2019, 111: 137-146. doi: 10.1016/j.compgeo.2019.03.010

    [7] 王道远, 崔光耀, 袁金秀, 等. 断裂黏滑错动下隧道减错措施作用效果模型试验研究[J]. 岩土工程学报, 2018, 40(8): 1515-1521. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201808021.htm

    WANG Dao-yuan, CUI Guang-yao, YUAN Jin-xiu, et al. Model tests on effect of dislocation reducing measures of stick-slip fault of tunnels[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(8): 1515-1521. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201808021.htm

    [8] 孙飞, 张志强, 易志伟. 正断层黏滑错动对地铁隧道结构影响的模型试验研究[J]. 岩土力学, 2019, 40(9): 3037-3044, 3053. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201908020.htm

    SUN Fei, ZHANG Zhi-qiang, YI Zhi-wei. Model experimental study on the influence of normal fault with stick-slip dislocation on the structure of subway tunnel[J]. Rock and Soil Mechanics, 2019, 40(9): 3037-3044, 3053. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201908020.htm

    [9]

    SABAGH M, GHALANDARZADEH A. Centrifugal modeling of continuous shallow tunnels at active normal faults intersection[J]. Transportation Geotechnics, 2020, 22: 100325.

    [10]

    YAN G M, SHEN Y S, GAO B, et al. Damage evolution of tunnel lining with steel reinforced rubber joints under normal faulting: an experimental and numerical investigation[J]. Tunnelling and Underground Space Technology, 2020, 97: 103223.

    [11] 刘学增, 林亮伦, 桑运龙. 逆断层粘滑错动对公路隧道的影响[J]. 同济大学学报(自然科学版), 2012, 40(7): 1008-1014. https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ201207009.htm

    LIU Xue-zeng, LIN Liang-lun, SANG Yun-long. Effect of thrust fault stick-slip rupture on road tunnel[J]. Journal of Tongji University(Natural Science), 2012, 40(7): 1008-1014. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ201207009.htm

    [12]

    CHANG Y Y, LEE C J, HUANG W C, et al. Evolution of the surface deformation profile and subsurface distortion zone during reverse faulting through overburden sand[J]. Engineering Geology, 2015, 184: 52-70.

    [13]

    BAZIAR M H, NABIZADEH A, MEHRABI R, et al. Evaluation of underground tunnel response to reverse fault rupture using numerical approach[J]. Soil Dynamics and Earthquake Engineering, 2016, 83: 1-17.

    [14] 焦鹏飞, 来弘鹏. 不同倾角逆断层错动对隧道结构影响理论分析[J]. 土木工程学报, 2019, 52(2): 106-117. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201902011.htm

    JIAO Peng-fei, LAI Hong-peng. Theoretical analysis on the influence of different dip angle reverse faults' dislocation on tunnel structure[J]. China Civil Engineering Journal, 2019, 52(2): 106-117. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201902011.htm

    [15]

    ZHAO K, CHEN W Z, YANG D S, et al. Mechanical tests and engineering applicability of fibre plastic concrete used in tunnel design in active fault zones[J]. Tunnelling and Underground Space Technology, 2019, 88: 200-208.

    [16]

    ZAHERI M, RANJBARNIA M, DIAS D, et al. Performance of segmental and shotcrete linings in shallow tunnels crossing a transverse strike-slip faulting[J]. Transportation Geotechnics, 2020, 23: 100333.

    [17]

    FOSSEN H. Structural Geology[M]. New York: Cambridge University Press, 2010.

    [18]

    ZHANG Z Q, CHEN F F, LI N, et al. Influence of fault on the surrounding rock stability of a tunnel: location and thickness[J]. Tunnelling and Underground Space Technology, 2017, 61: 1-11.

    [19] 宋佳佳, 孙建孟, 王敏, 等. 断层内部结构研究进展[J]. 地球物理学进展, 2018, 33(5): 1956-1966. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ201805024.htm

    SONG Jia-jia, SUN Jian-meng, WANG Min, et al. Research progress in the internal structure of the fault[J]. Progress in Geophysics, 2018, 33(5): 1956-1966. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ201805024.htm

    [20]

    HOEK E, BROWN E T. Practical estimates of rock mass strength[J]. International Journal of Rock Mechanics and Mining Sciences, 1997, 34(8): 1165-1186.

    [21] 中国地震局地壳应力研究所. 2008年汶川8.0级地震科学考察图集[M]. 北京: 地震出版社, 2009.

    The Institute of Crustal Dynamics, CEA. Scientific Investigation Album of the Wenchuan MS 8.0 Earthquake of 2008[M]. Beijing: Seismological Press, 2009. (in Chinese)

    [22]

    Code of Practice for Temporary Works Procedures and the Permissible Stress Design of Falsework: BS5975—2008[S]. 2008.

    [23] 公路隧道设计细则:JTG/T D70—2010[S]. 2010.

    Guidelines for Design of Highway Tunnel: JTG/T D70—2010[S]. 2010. (in Chinese)

    [24] 地下结构抗震设计标准:GB/T 51336—2018[S]. 2019.

    Standard for Seismic Design of Underground Structures: GB/T 51336—2018[S]. 2019. (in Chinese)

    [25]

    CHEN Z Y, WEI J S. Correlation between ground motion parameters and lining damage indices for mountain tunnels[J]. Natural Hazards, 2013, 65: 1683-1702.

    [26]

    WANG Z Z, ZHANG Z. Seismic damage classification and risk assessment of mountain tunnels with a validation for the 2008 Wenchuan earthquake[J]. Soil Dynamics and Earthquake Engineering, 2013, 45: 45-55.

    [27]

    Federal Emergency Management Agency. Multi-Hazard Loss Estimation Methodology: Earthquake Model HAZUS-MH 2.1 User Manual[R]. Washington D C: FEMA Mitigation Division, 2010.

    [28] 刘国庆, 肖明, 陈俊涛. 基于增量动力分析的隧洞结构抗震性能评估[J]. 工程科学与技术, 2019, 51(3): 92-100. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201903012.htm

    LIU Guo-qing, XIAO Ming, CHEN Jun-tao. Seismic performance assessment of tunnel structure based on incremental dynamic analysis[J]. Advanced Engineering Sciences, 2019, 51(3): 92-100. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201903012.htm

  • 期刊类型引用(19)

    1. 杜岩,张洪达,谢谟文,蒋宇静,张明,贾北凝. 矿山采空区边坡动态稳定性评价方法. 工程科学学报. 2025(02): 215-223 . 百度学术
    2. 潘网生,马宗源,傅良同. 贵州三叠纪垄头组软岩优势结构面研究意义和进展. 绿色科技. 2024(06): 236-240 . 百度学术
    3. 陶通铭. 采动作用下岩质斜坡变形破坏机制研究——以贵州龙场崩塌为例. 中国水运(下半月). 2024(11): 113-114+145 . 百度学术
    4. 陶通铭. 采动作用下岩质斜坡变形破坏机制研究——以贵州龙场崩塌为例. 中国水运. 2024(22): 113-114+145 . 百度学术
    5. 梁博,杨更社,冯伟,潘振兴,孙杰龙,刘慧,陈奇. 冻融诱发平面滑移型岩质边坡失稳模型试验研究. 西安科技大学学报. 2024(06): 1118-1126 . 百度学术
    6. 于群群,孙朝燚. 顺层岩质边坡滑剪破坏规律研究. 河南城建学院学报. 2023(01): 15-21 . 百度学术
    7. 杨忠平,向宫固,赵茜,刘新荣,赵亚龙. 水动力-溶蚀作用下灰岩结构面剪切力学特性. 岩土工程学报. 2023(08): 1555-1563 . 本站查看
    8. 李华. 基于无人机航摄三维模型的地质信息提取和数值模拟应用. 中国水运(下半月). 2023(09): 31-33 . 百度学术
    9. 牛犇,冯春,丛俊余,孙子正,张一鸣. 基于CDEM颗粒流的三维高速远程滑坡成灾范围分析. 岩石力学与工程学报. 2023(S2): 4018-4027 . 百度学术
    10. 杨小龙,王刚. 滑坡堆积体反粒序现象的离散元数值分析. 工程地质学报. 2023(06): 1941-1950 . 百度学术
    11. 李华. 基于无人机航摄三维模型的地质信息提取和数值模拟应用. 中国水运. 2023(18): 31-33 . 百度学术
    12. 穆成林,裴向军,王睿,王超. 基于物理模型试验的含多层软弱夹层顺层开挖高边坡变形破坏特征分析. 中国地质灾害与防治学报. 2022(03): 61-67 . 百度学术
    13. 朱彦鹏,杜一博,杨校辉,张卫雄,朱鋆川. 甘肃舟曲河那滑坡变形特征及孕灾机理. 科学技术与工程. 2022(25): 10884-10895 . 百度学术
    14. 朱赛楠,殷跃平,王猛,朱茂,王晨辉,王文沛,李俊峰,赵慧. 金沙江结合带高位远程滑坡失稳机理及减灾对策研究——以金沙江色拉滑坡为例. 岩土工程学报. 2021(04): 688-697 . 本站查看
    15. 易连兴,李瑜. 岩溶及水动力对鸡尾山滑坡影响作用研究. 工程地质学报. 2021(03): 583-592 . 百度学术
    16. 刘宁波,钟立力,龙森. 山区公路单斜坡回头曲线路段滑坡成因及治理方案研究. 路基工程. 2020(03): 190-195 . 百度学术
    17. 李巧学,周洪福,冉涛,铁永波. 中倾侧向坡中一种岩体破坏过程的模拟研究. 路基工程. 2020(04): 62-65+72 . 百度学术
    18. 崔芳鹏,李滨,杨忠平,吴乐乐,李宁,彭健全. 贵州纳雍普洒滑坡动力触发机制离散元模拟分析. 中国岩溶. 2020(04): 524-534 . 百度学术
    19. 何忠明,杨煜,曾新发,刘森峙. 土洞演化过程中路基变形的响应分析. 中南大学学报(自然科学版). 2018(12): 3068-3076 . 百度学术

    其他类型引用(28)

图(17)  /  表(2)
计量
  • 文章访问数:  321
  • HTML全文浏览量:  25
  • PDF下载量:  183
  • 被引次数: 47
出版历程
  • 收稿日期:  2020-02-02
  • 网络出版日期:  2022-12-07
  • 刊出日期:  2020-09-30

目录

/

返回文章
返回