深部节理岩体开挖卸荷动力破坏机制的离散元研究

薄英鋆; 王华宁; 蒋明镜

doi:10.11779/CJGE2020S2035

深部节理岩体开挖卸荷动力破坏机制的离散元研究 English Version

薄英鋆^1,,
王华宁^2, ,,
蒋明镜^{1, 3, 4}

1.
同济大学地下建筑与工程系,上海　200092
2.
同济大学航空航天与力学学院,上海　200092
3.
天津大学建筑工程学院土木工程系,天津　300350
4.
天津大学水利工程仿真与安全国家重点实验室,天津　300350

基金项目:

国家自然科学基金项目 11872281

国家自然科学基金项目 51639008

同济大学土木工程防灾国家重点实验室自主课题 SLDRCE19-A-06

详细信息

作者简介:
薄英鋆（1995— ）,男,上海人,硕士研究生,从事岩石开挖的动力效应研究。E-mail：1832349@tongji.edu.cn

通讯作者:
王华宁, E-mail：wanghn@tongji.edu.cn

中图分类号: TU43
计量
- 文章访问数: 196
- HTML全文浏览量: 51
- PDF下载量: 247
出版历程
- 收稿日期: 2020-08-06
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-10-31

Cracking mechanism of deep jointed rock mass induced by transient excavation via DEM

1.
Department of Geotechnical Engineering, College of Civil Engineering, Tongji University, Shanghai 200092, China
2.
School of Aerospace Engineering and Applied Mechanics, Tongji University, Shanghai 200092, China
3.
Department of Civil Engineering, School of Civil Engineering, Tianjin University, Tianjin 300350, China
4.
State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300350, China

摘要

摘要: 对于高地应力下的岩体开挖,岩体自身的缺陷以及开挖诱发的动力效应对围岩的破坏有十分重要的影响。对节理的分布进行简化与量化,采用二维离散元方法模拟了断续节理岩体中隧道开挖过程,研究不同连通率下岩体的宏观与微观动力破坏模式,以及节理面附近点处的动力响应。结果表明考虑动力效应后,岩体破坏范围以及程度显著增大；高连通率下动力计算模式更容易产生翼型裂纹的聚合；节理面的存在放大了岩桥中心的动力效应,并对节理面后的岩体起到缓冲作用,但是节理面之前的岩体将会受到反射应力波的影响。
- 离散元 /
- 隧道 /
- 动力效应 /
- 断续节理
Abstract: For the excavation of rock mass under high stress, the defects of the rock mass and the dynamic effects induced by the excavation have a great impact on the damage of the surrounding rock. The distribution of joints is simplified and quantified. Based on the two-dimensional discrete element method, the process of excavation of rock masses with non-persistent joints is simulated. The macro-and micro-dynamic failure modes of the rock mass under different connectivity rates and the dynamic response near the joint are studied. The results show that after considering the dynamic effect, the range and degree of damage increase significantly. The dynamic model is more likely to produce wing crack coalescence under high connectivity. The joint amplifies the dynamic effect at the rock bridge center and acts as a cushion to absorb the shock behind the joint, but the rock mass before the joint will be affected by the reflected stress waves.
- discrete element method /
- tunnel /
- dynamic effect /
- non-persistent joint

HTML全文

0. 引言

土石混合体（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]对土石混合体-台阶状界面接触面开展了数值模拟剪切试验,研究了含石率和土石混体最大块石粒径对其接触界面的力学特性的影响。但接触面粗糙度对土石混合体-基岩接触面剪切力学特性的影响尚不多见,仍需深入研究。鉴于此,本文通过土石混合体-基岩接触面室内大型直剪试验,较系统地研究了接触面粗糙度对填方体-下伏基岩界面的剪切强度及变形特征的影响。可为高填方体或堆积体边坡稳定性分析提供参考。

1. 室内大型直剪试验

1.1 试验设备

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

图 2 ZJ50-2G大型粗粒土压缩直剪仪示意图

Figure 2. ZJ50-2G large coarse soil compression direct shear apparatus

下载: 全尺寸图片幻灯片

1.2 试样制备

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

表 1 土石混合体及灰岩基本物理参数指标

Table 1. Basic physical parameter indexes of soil-rock aggregate

土体类型	物理参数指标
土体类型	干密度/(g·m^-3)	孔隙比	天然含水率/%	天然密度/(kg·m^-3)	c/kPa	φ/(°)	弹性模量/GPa	单轴抗压强度/MPa
土石混合体	1788	0.24	9.32	2110	23.9	10.54	—	—
灰岩	—	—	—	2730	1433	35.67	29.14	68.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/mm	J斜面α/(°)
1	1∶2	25	50	0.399	26.56
2	1∶1.75	25	43.75	0.431	29.74
3	1∶1.5	25	37.44	0.462	33.69

下载: 导出CSV

| 显示表格

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

图 6 加工后的岩石试件

Figure 6. Processed rock samples

下载: 全尺寸图片幻灯片

1.3 试验方案

将台阶状灰岩试件放入下剪切盒中,在上剪切盒中放入土石混合体。室内试验所用土石混合体试验级配如图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
1	0.399	200,400,600,800
2	0.431	200,400,600,800
3	0.462	200,400,600,800

下载: 导出CSV

| 显示表格

2. 试验结果与分析

2.1 剪应力-剪切位移曲线分析

不同法向应力水平作用下,土石混合体与不同粗糙度试块接触面的剪应力-剪切位移曲线如图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的倾角 $α$ 越大,在产生相同切向变形下,将产生更大的竖向变形,需要更大的克服法向应力做功,因此剪切应力增长更快,在低法向应力作用下,能更快地达到峰值。

2.2 粗糙度对接触面附近块石破碎的影响

接触面附近的块石破碎量反映了土石混合填方体中块石与基岩接触面相互咬合的紧密程度,咬合作用越强烈,则在剪切时块石越容易被剪碎。如图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发生爬坡运动时,斜面越陡其法向位移分量越大,剪切破坏区厚度增加,剪切破坏区内有更多的块石相互接触,形成强度较高的块石骨架,随着剪切位移不断增加,紧密咬合的块石更容易被台阶剪断,所以试件粗糙度增加时,块石破碎的现象会加剧。完全破碎或者部分的块石数量可以从侧面体现出剪切面上块石与块石、块石与台阶咬合作用和嵌固作用的强弱程度,在剪切过程中咬合作用和嵌固作用力越强,块石越不易翻越相邻块石或台阶,则块石容易被剪碎。

2.3 粗糙度对接触面抗剪强度的影响

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

表 4 界面抗剪强度

Table 4. Interface shear strengths

项目		试件编号
项目		1	2	3
抗剪强度	$σ_{n}$ =200 kPa	179.0500	224.2300	264.8700
	$σ_{n}$ =400 kPa	356.9700	401.5600	460.8900
	$σ_{n}$ =600 kPa	490.3700	536.6700	575.6700
	$σ_{n}$ =800 kPa	605.5600	694.8200	794.5600
粗糙度		0.3990	0.4310	0.4620
相关系数R		0.9899	0.9975	0.9874
表观黏聚力/kPa		54.7500	77.6000	98.0400
内摩擦角/(°)		35.2400	37.7200	40.4300

下载: 导出CSV

| 显示表格

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

2.4 粗糙度对接触面剪切强度参数的影响

不同接触面粗糙度下抗剪强度-法向应力关系曲线如图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],本文认为表观黏聚力产生的原因是块石与块石之间、块石与台阶之间的咬合作用力。随着粗糙度的增加,块石与块石之间、块石与台阶之间的咬合作用力随之增强,导致接触面的表观黏聚力显著提高。需要补充的一点是,接触面上不仅有表观黏聚力,还有台阶凹槽内土石混合体与剪切面上部土石混合体中土颗粒间的黏聚力,但由于本试验土石混合体含石率较大,作为胶结材料的土颗粒之间的接触较少,所以这部分黏聚力很小,接触面上主要由表观黏聚力占主导作用。

2.5 粗糙度对剪切带的影响

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

图 15 不同粗糙度下剪切带示意图

Figure 15. Schematic diagram of shear zone with different roughnesses

下载: 全尺寸图片幻灯片

3. 结论

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

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

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

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

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

图 1 DEM模型示意图

Figure 1. Schematic diagram of DEM specimen

下载: 全尺寸图片幻灯片

图 2 不同连通率胶结破坏分布图

Figure 2. Distribution of cracks under different joint persistencies

下载: 全尺寸图片幻灯片

图 3 不同连通率下胶结破坏数量

Figure 3. Number of cracks under different joint persistencies

下载: 全尺寸图片幻灯片

图 4 破坏模式示意图

Figure 4. Schematic diagram of failure modes

下载: 全尺寸图片幻灯片

图 5 靠近洞口的破坏（考虑动力效应）

Figure 5. Failures around tunnel (dynamic solution)

下载: 全尺寸图片幻灯片

图 6 远离洞口的破坏（考虑动力效应）

Figure 6. Failures away from tunnel (dynamic solution)

下载: 全尺寸图片幻灯片

图 7 靠近洞口的破坏（不考虑动力效应）

Figure 7. Failures around tunnel (static solution)

下载: 全尺寸图片幻灯片

图 8 远离洞口的破坏（不考虑动力效应）

Figure 8. Failures away from tunnel (static solution)

下载: 全尺寸图片幻灯片

图 9 近孔口的节理布置示意图

Figure 9. Schematic diagram of joints around tunnel

下载: 全尺寸图片幻灯片

图 10 第一组观测点颗粒速度时程变化

Figure 10. Change of velocity of first group of particles

下载: 全尺寸图片幻灯片

图 11 第二组观测点颗粒速度时程变化

Figure 11. Change of velocity of second group of pariticles

下载: 全尺寸图片幻灯片

图 12 第三组观测点颗粒速度时程变化

Figure 12. Change of velocity of third group of particles

下载: 全尺寸图片幻灯片

表 1 DEM岩石试样模型参数

Table 1 Material parameters of DEM rock specimens

类型	微观参数	数值
颗粒部分	颗粒密度 $ρ$ /(kg·m^-3)	2700
	颗粒法向刚度 $k_{n}$ /(N·m^-1)	1.5×10¹⁰
	颗粒切向刚度 $k_{s}$ /(N·m^-1)	1×10¹⁰
	颗粒摩擦系数 $μ$	0.7
	颗粒抗转动系数 $β$	1.2
胶结部分	最大胶结厚度 $h_{\max}$ /m	1×10^-5
	胶结弹性模量 $E_{b}$ /N	2.3×10⁸
	胶结抗压强度 $σ_{c}$ /Pa	2.7×10⁹
	胶结抗拉强度 $σ_{t}$ /Pa	2.7×10⁸
	胶结延伸率 $s_{p}$	0.15

下载: 导出CSV

表 2 无节理岩石宏观参数

Table 2 Macro-parameters of intact rock

对比	弹性模量^[13]/ GPa	泊松比^[13]	单轴抗压强度^[14]/MPa	巴西劈裂强度^[15]/MPa	内摩擦角^[13]/(°)
试验	10	0.152	38	2.55	40.1
DEM	10.2	0.15	36	2.552	39.7
误差/%	2	1.32	5	0.07	1.00

下载: 导出CSV

参考文献(18)

[1]	钱七虎. 岩爆、冲击地压的定义、机制、分类及其定量预测模型[J]. 岩土力学, 2014, 35(1): 1-6. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201401001.htm QIAN Qi-hu. Definition, mechanism, classification and quantitative forecast model for rockburst and pressure bump Rock and Soil Mechanics[J]. Rock and Soil Mechanics, 2014, 35(1): 1-6. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201401001.htm
[2]	SHEMYAKIN E I, FISENKO G L, KURLENYA M V, et al. Zonal disintegration of rocks around underground workings, Part 1: data of in situ observations[J]. Soviet Mining, 1986, 22(3): 157-168. doi: 10.1007/BF02500863
[3]	MIKLOWITZ J. Plane-stress unloading waves emanating from a suddenly punched hole in a stretched elastic plate[J]. Journal of Applied Mechanics, 1960, 27(1): 165-171. doi: 10.1115/1.3643892
[4]	CAI M. Influence of stress path on tunnel excavation response: Numerical tool selection and modeling strategy[J]. Tunnelling and Underground Space Technology, 2008, 23(6): 618-628. doi: 10.1016/j.tust.2007.11.005
[5]	LI X, CAO W, ZHOU Z, et al. Influence of stress path on excavation unloading response[J]. Tunnelling and Underground Space Technology, 2014, 42: 237-246. doi: 10.1016/j.tust.2014.03.002
[6]	FAN Y, LU W B, YAN P, et al. Transient characters of energy changes induced by blasting excavation of deep-buried tunnels[J]. Tunnelling and Underground Space Technology, 2015, 49: 9-17. doi: 10.1016/j.tust.2015.04.003
[7]	CHEN M, LU W B, YAN P, et al. Blasting excavation induced damage of surrounding rock masses in deep-buried tunnels[J]. KSCE Journal of Civil Engineering, 2016, 20(2): 933-942. doi: 10.1007/s12205-015-0480-3
[8]	YAN P, LU W, CHEN M, et al. Contributions of in-situ stress transient redistribution to blasting excavation damage zone of deep tunnels[J]. Rock Mechanics and Rock Engineering, 2015, 48(2): 715-726. doi: 10.1007/s00603-014-0571-3
[9]	谢和平, 高峰, 鞠杨, 等. 深部开采的定量界定与分析[J]. 煤炭学报, 2015, 40(1): 1-10. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201501001.htm XIE He-ping, GAO Feng, JU Yang, et al. Quantitative definition and investigation of deep mining[J]. Journal of China Coal Society, 2015, 40(1): 1-10. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201501001.htm
[10]	蒋明镜, 金树楼, 张宁. 不同胶结尺寸的粒间胶结强度统一表达式[J]. 岩土力学, 2015, 36(9): 2451-2457, 2466. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201509004.htm JIANG Ming-jing, JIN Shu-lou, ZHANG Ning. Unified expression for bonding strength of cemented granules with different bond sizes[J]. Rock and Soil Mechanics, 2015, 36(9): 2451-2457, 2466. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201509004.htm
[11]	IVARSD MPOTYONDYDPIERCEMThe smooth-joint contact model8th World Congress on Computational Mechanics 5th European Congress on Computational Methods in Applied Sciences and Engineering200823 IVARS D M, POTYONDY D, PIERCE M, et al. The smooth-joint contact model[J]. 8th World Congress on Computational Mechanics 5th European Congress on Computational Methods in Applied Sciences and Engineering, 2008: 2-3.
[12]	苏鹏. 砂岩、花岗岩阻尼特性及机制的试验研究[D]. 西安: 西安理工大学, 2017. SU Peng. Experimental Study on Damping Characteristic and Mechanism of Sandstone and Granite[D]. Xi'an: Xi'an University of Technology, 2017. (in Chinese)
[13]	邹航, 刘建锋, 边宇, 等. 不同粒度砂岩力学和渗透特性试验研究[J]. 岩土工程学报, 2015, 37(8): 1462-1468. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201508020.htm ZOU Hang, LIU Jian-feng, BIAN Yu, et al. Experimental study on mechanical and permeability properties of sandstone with different granularities[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(8): 1462-1468. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201508020.htm
[14]	TAN Y L, NING J G, LI H T. In situ explorations on zonal disintegration of roof strata in deep coalmines[J]. International Journal of Rock Mechanics and Mining Sciences, 2012, 49: 113-124. doi: 10.1016/j.ijrmms.2011.11.015
[15]	WANG J J, LI Y Q, JIAN F X, et al. Rate dependence of splitting tensile behaviors of sandstone and mudstone[J]. Geotechnical and Geological Engineering, 2019, 37(4): 3469-3475. doi: 10.1007/s10706-018-00797-7
[16]	JIANG M J, KONRAD J M, LEROUEIL S. An efficient technique for generating homogeneous specimens for DEM studies[J]. Computers and Geotechnics, 2003, 30(7): 579-597. doi: 10.1016/S0266-352X(03)00064-8
[17]	BAUD P, MEREDITH P, TOWNEND E. Permeability evolution during triaxial compaction of an anisotropic porous sandstone: Permeability Evolution in Sandstone[J]. Journal of Geophysical Research: Solid Earth, 2012, 117(B5): 5023.
[18]	ZHOU X, SHENG Q, CUI Z. Dynamic boundary setting for discrete element method considering the seismic problems of rock masses[J]. Granular Matter, 2019, 21(3): 66.

施引文献(16)

期刊类型引用(7)

1.	袁炀，彭丰，袁希强，徐洪伟. 基于连续-非连续数值模拟方法的层状围岩大变形分析. 甘肃水利水电技术. 2025(01): 38-42+48 . 百度学术
2.	侯公羽，邵耀华，刘春雷，王一哲，陈钦煌，胡志宇，赵庆儒. 节理巷道围岩卸荷变形试验研究. 岩土工程学报. 2024(03): 616-623 . 本站查看
3.	李壮，徐奴文，孙志强，刘军，李彪，孙悦鹏，朱建林. 基于微震监测与数值模拟的高应力软岩隧道围岩大变形特征分析. 岩石力学与工程学报. 2024(11): 2725-2737 . 百度学术
4.	邓华锋，陈天楠，李建林，肖瑶，张恒宾，李涛，陈兴周. 峰前卸荷损伤灰岩的再加载力学性能研究. 岩石力学与工程学报. 2023(06): 1301-1311 . 百度学术
5.	王梦，范立峰. 卸载效应对节理裂隙岩体内应力波能量耗散影响研究. 应用力学学报. 2022(05): 869-878 . 百度学术
6.	董平，刘婷婷，李新平，丁鹿阳，王刚. 冲击荷载作用下非贯通节理岩体细观破坏模式研究. 爆破. 2021(04): 17-25 . 百度学术
7.	陈锐，陈高生，杨静，黄锋，刘星辰. 节理岩体小净距地铁隧道暗挖施工响应数值分析. 现代隧道技术. 2021(S1): 232-239 . 百度学术

其他类型引用(9)

资源附件(0)

图(12) / 表(2)

计量

文章访问数: 196
HTML全文浏览量: 51
PDF下载量: 247
被引次数: 16

0. 引言
1. 室内大型直剪试验
1.1 试验设备
1.2 试样制备
1.3 试验方案
2. 试验结果与分析
2.1 剪应力-剪切位移曲线分析
2.2 粗糙度对接触面附近块石破碎的影响
2.3 粗糙度对接触面抗剪强度的影响
2.4 粗糙度对接触面剪切强度参数的影响
2.5 粗糙度对剪切带的影响
3. 结论

0. 引言
1. 室内大型直剪试验
1.1 试验设备
1.2 试样制备
1.3 试验方案
2. 试验结果与分析
2.1 剪应力-剪切位移曲线分析
2.2 粗糙度对接触面附近块石破碎的影响
2.3 粗糙度对接触面抗剪强度的影响
2.4 粗糙度对接触面剪切强度参数的影响
2.5 粗糙度对剪切带的影响
3. 结论

参考文献(18)

施引文献

资源附件(0)

深部节理岩体开挖卸荷动力破坏机制的离散元研究 English Version

作者简介: 薄英鋆（1995— ）,男,上海人,硕士研究生,从事岩石开挖的动力效应研究。E-mail：1832349@tongji.edu.cn

通讯作者: 王华宁, E-mail：wanghn@tongji.edu.cn

计量

出版历程

Cracking mechanism of deep jointed rock mass induced by transient excavation via DEM

0. 引言

1. 室内大型直剪试验

1.1 试验设备

1.2 试样制备

1.3 试验方案

2. 试验结果与分析

2.1 剪应力-剪切位移曲线分析

2.2 粗糙度对接触面附近块石破碎的影响

2.3 粗糙度对接触面抗剪强度的影响

2.4 粗糙度对接触面剪切强度参数的影响

2.5 粗糙度对剪切带的影响

3. 结论

期刊类型引用(7)

其他类型引用(9)

计量

出版历程

目录

0. 引言

1. 室内大型直剪试验

1.1 试验设备

1.2 试样制备

1.3 试验方案

2. 试验结果与分析

2.1 剪应力-剪切位移曲线分析

2.2 粗糙度对接触面附近块石破碎的影响

2.3 粗糙度对接触面抗剪强度的影响

2.4 粗糙度对接触面剪切强度参数的影响

2.5 粗糙度对剪切带的影响

3. 结论

作者简介:
薄英鋆（1995— ）,男,上海人,硕士研究生,从事岩石开挖的动力效应研究。E-mail：1832349@tongji.edu.cn

通讯作者:
王华宁, E-mail：wanghn@tongji.edu.cn