土石混合体—基岩界面剪切力学特性试验研究

杨忠平; 蒋源文; 李诗琪; 李进; 胡元鑫

doi:10.11779/CJGE202010021

土石混合体—基岩界面剪切力学特性试验研究 English Version

杨忠平^{1, 2, 3,},
蒋源文^{1, 2, 3},
李诗琪^{1, 2, 3},
李进^{1, 2, 3},
胡元鑫⁴

1.
重庆大学土木工程学院,重庆　400045
2.
山地城镇建设与新技术教育部重点实验室(重庆大学),重庆　400045
3.
库区环境地质灾害防治国家地方联合工程研究中心（重庆）,重庆　400045
4.
渝武机场建设指挥部,重庆　401120

基金项目:

国家重点研发计划项目 2018YFC1504802

中央高校基本科研业务经费项目 2019CDXYTM0032

国家自然科学基金项目 41772306

详细信息

作者简介:
杨忠平(1981—),男,教授,博士生导师,主要从事环境岩土与边坡稳定性方面的教学与研究工作。E-mail：yang-zhp@163.com

中图分类号: TU431
计量
- 文章访问数: 445
- HTML全文浏览量: 33
- PDF下载量: 323
出版历程
- 收稿日期: 2020-01-07
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-09-30

Experimental study on shear mechanical properties of soil-rock mixture-bedrock interface

YANG Zhong-ping^{1, 2, 3,},
JIANG Yuan-wen^{1, 2, 3},
LI Shi-qi^{1, 2, 3},
LI Jin^{1, 2, 3},
HU Yuan-xin⁴

1.
School of Civil Engineering, Chongqing University, Chongqing 400045, China
2.
Key Laboratory of New Technology for Construction of Cities in Mountain Area (Chongqing University), Ministry of Education, Chongqing 400045, China
3.
National Joint Engineering Research Center for Prevention and Control of Environmental Geological Hazards in the TGR Area, Chongqing University, Chongqing 400045, China
4.
Headquarters of Chongqing Yuwu Airport, Chongqing 401120, China

摘要

摘要: 填方体-下伏基岩接触面间的剪切强度是控制高填方体或堆积体边坡稳定性的重要因素,界面强度参数取值是高填方工程设计的重要参数之一。通过较系统地室内大型直剪试验探讨了接触面粗糙度对土石混合料-基岩接触面剪切力学特性的影响。结果表明：在低法向应力作用下,剪应力-剪切位移曲线前期呈现出应变硬化现象,后期呈现出塑性应变现象,且接触面粗糙度越大接触面发生剪切破坏时变形越小；在高法向应力作用下,曲线呈现出应变硬化现象,无明显峰值；相同法向应力水平作用下,接触面粗糙度越大,土石混合体-基岩接触面剪切刚度越大。剪切界面上块石的破碎形态可分为完全破碎、部分破碎和表面磨损3种,随着接触面粗糙度的增加,剪切界面上块石的破碎总数也增加。接触面的抗剪强度、内摩擦角和表观黏聚力随着接触面粗糙度的增加而增大,相比于内摩擦角,接触面的表观黏聚力增大较为明显。接触面粗糙度对剪切带宽度有影响作用,表现为接触面粗糙度越大,剪切带越宽。
- 土石混合体 /
- 基岩界面 /
- 粗糙度 /
- 剪切特性 /
- 块石破碎
Abstract: The shear strength of the interface between the fill and the underlying bedrock is an important factor to control the stability of high fill or accumulation slopes. The value of the interface strength parameter is one of the important parameters for the design of high backfills. The effect of the contact surface roughness on the shear mechanical properties of the soil-rock mixture-bedrock contact surface is explored through the systematic large-scale indoor direct shear tests.The test results show that under the action of low normal stress, the shear stress-shear displacement curve shows strain hardening in the early stage and plastic strain in the later stage, and the greater the roughness of the contact surface, the more the contact surface deforms when shear failure occurs. Under the action of high normal stress, the curve shows strain hardening without obvious peaks. Under the same normal stress level, the greater the contact surface roughness, the greater the shear stiffness of the soil-rock mixture-base rock interface. The crushing morphology of the rock at the shear interface includes three types: complete crushing, partial crushing, and surface abrasion. As the contact surface roughness increases, the total number of rock crushing at the shear interface also increases. The shear strength, internal friction angle and apparent cohesion of the contact surface increase with the increase of the roughness of the contact surface. Compared with the internal friction angle, the apparent cohesion of the contact surfaces increases significantly. The roughness of the contact surface has an effect on the width of the shear band, which shows that the larger the roughness of the contact surface is, the wider the shear band is.
- soil-rock mixture /
- bedrock interface /
- roughness /
- shear characteristic /
- broken stone

HTML全文

0. 引言

岩土工程、水利工程、道路桥梁工程等的附存环境或对象相当大一部分是由离散形式的颗粒材料组成,如砂土、碎石、粗粒土等,可统称为岩土颗粒材料。岩土颗粒材料的力学特性直接影响地（路）基、堤防、堆石坝等工程的设计、施工和运行。颗粒材料看似简单,但其物理力学性质复杂,存在压硬性、剪胀性、应力路径相关性等复杂力学特性^[1-2]。这些复杂的力学行为与其离散性、多尺度和能量耗散机制有关^[1]。对颗粒材料进行宏观唯象建模不足以反映其所有力学特性^[3-4],颗粒材料的微观动力学行为、接触网络、介观结构、细观组构等为揭示颗粒材料复杂的宏细观力学特性提供了新的途径^[5-7]。

细观尺度下颗粒材料的研究手段主要包括物理实验和数值模拟两大类。物理实验方面,常用的观测手段有光弹实验、扫描电镜技术、CT（computed tomography）扫描技术、DIC（digital image correlation）成像技术等^[4]。数值分析方法有离散单元法（discrete element method,DEM）^[5-6]、连续-离散耦合分析方法（combined finite and discrete method,FDEM）^[8-9]等。其中,颗粒离散单元法（DEM）发展迅速,计算简便,已广泛应用于岩土类颗粒材料的研究中,经过不断改进能够对颗粒形状^[10]、颗粒破碎^[11]、以及颗粒材料的临界状态^[12-13]、剪胀^[13]、细观接触^[14]和组构^{[6, 14]}等特性进行表征。

岩土颗粒材料是多体相互作用的耗散体系,颗粒间摩擦是耗散外力功的主要途径之一,是颗粒材料宏观剪切强度的主要来源之一,也是其变形和破坏特性的重要内部控制因素^[4]。颗粒间摩擦特性会影响颗粒间接触力的大小分布和颗粒重排列,进而引起颗粒体系宏观力学特性的变化。Dai等^[15]采用玻璃珠进行了一系列直剪试验,发现颗粒间摩擦对颗粒材料剪切行为的影响显著。Sandeep等^[16]指出颗粒间的摩擦系数与杨氏模量和表面粗糙度成反比。Senetakis等^[17]采用破碎的石灰石颗粒进行微观力学滑动实验并对颗粒间的摩擦系数进行量化。

细观数值模拟方面,Rothenburg等^[12]采用离散单元法研究了不同摩擦系数的试样孔隙比与配位数之间的关系；Antony等^[5]分析了摩擦特性对颗粒材料宏观抗剪强度的影响；Huang等^[18]评估了颗粒材料的临界状态对颗粒间摩擦的敏感性；Barreto等^[19]认为摩擦系数影响了强接触力链的稳定性。尽管关于颗粒间摩擦系数对颗粒材料宏细观力学特性的影响研究较多,但摩擦特性对颗粒材料在三维复杂加载路径下力学行为的影响,以及不同颗粒间摩擦系数下各类细观结构对宏观抗剪强度和体积变形的作用机制,还有待进一步的研究。

颗粒间接触形成的接触力网络是一个高度非均质、强非线性的的复杂网络。颗粒体系内大部分的力是通过相对稀疏的链状网络（称之为“力链”）传递^[20]。Radjai等^[21]、Thornton等^[22]将颗粒体系的接触网络划分为两个互补的子接触网络：强接触网络和弱接触网络。大于系统平均接触力的接触形成的强接触网络是传递力的主要通道；弱接触网络中的接触力小于系统平均接触力,为强接触网络的稳定提供支撑^[22]。分别对强弱接触网络的配位数、组构等细观参量进行分析,是研究颗粒材料强度和变形的有效途径^{[5, 23-24]}。

本文将基于岩土颗粒材料真三轴离散元数值试验,分析细观滑动摩擦与颗粒体系强度、变形之间的关系,研究不同摩擦系数下强、弱接触网络的配位数分布、接触力、组构张量及各向异性的特征,探讨强、弱接触网络对颗粒材料宏观应力-应变关系的贡献,揭示细观尺度摩擦特性对宏观力学行为的作用机制。

1. DEM数值试验

1.1 数值试验及参数

本文利用开源软件LIGGGHTS进行离散元模拟^[25]。进行真三轴加载路径的DEM模拟之前,需要制备各向同性且具有一定围压的立方体数值试样。制样时,在400 mm×400 mm×400 mm的立方体空间中随机生成粒径服从Gaussian分布的圆球颗粒,共计31253个。为了避免制样产生各向异性,在试样的各个方向采用位移控制等速地压缩试样直至目标大小。随后,施加三向等压应力进行各向同性固结,直至达到预定的围压值0.5 MPa,得到试样的初始孔隙比为0.582。采用Hertz-Mindlin非线性接触模型^[26]描述颗粒间的接触力,颗粒间滑移服从库仑摩擦定律,即颗粒间发生滑动的条件为｜F_t｜≥μ｜F_n｜,μ为颗粒间的滑动摩擦系数。加载过程中的μ值分别取0.05,0.1,0.2,0.3,0.4,0.5。其它DEM模拟参数如表1所示。

表 1 数值模拟细观参数

Table 1. Microscopic parameters of numerical simulation

参数	值
密度ρ/(kg·m^-3)	2600
颗粒粒径d/mm	2.25~8.25
弹性模量E/GPa	65.0
泊松比 $ν$	0.12
恢复系数e	0.95
滑动摩擦系数μ	0.05,0.1,0.2,0.3,0.4,0.5

下载: 导出CSV

| 显示表格

1.2 真三轴加载路径

采用等p（围压）等b（中主应力系数）加载路径对各向同性试样进行加载。加载时,为了消除边界效应,避免应变局部化,采用周期性边界^[27]。加载过程中,轴向加载速度为-0.01 mm/s。采用周期性边界的伺服机制,使试样的应力状态达到预定的值以满足真三轴应力路径^{[13, 27]},试样中主应力 $σ_{2}$ 和小主应力 $σ_{3}$ 的控制条件^[13]分别为

$σ_{2} = b σ_{1} + (1 - b) \frac{3 p - σ_{1} (1 + b)}{(2 - b)}$ ,

(1)

$σ_{3} = \frac{3 p - σ_{1} (1 + b)}{(2 - b)}$ 。

(2)

颗粒材料的广义剪应力q,平均静水压力p,体积应变 $ε_{v}$ ,广义剪应变 $ε_{d}$ 均可由应力应变张量主值导出：

$p = \frac{σ_{1} + σ_{2} + σ_{3}}{3}$ ,

(3)

$q = \frac{1}{\sqrt{2}} \sqrt{{(σ_{1} - σ_{2})}^{2} + {(σ_{2} - σ_{3})}^{2} + {(σ_{3} - σ_{1})}^{2}}$ ,

(4)

$ε_{v} = ε_{1} + ε_{2} + ε_{3}$ ,

(5)

$ε_{d} = \frac{\sqrt{2}}{3} \sqrt{[{(ε_{1} - ε_{2})}^{2} +(ε_{2} - ε_{3})^{2} + {(ε_{3} - ε_{1})}^{2}]}$ 。

(6)

2. 细观摩擦与宏观强度

2.1 加载过程中的应力应变

不同的摩擦系数下,剪应力q,体积应变 $ε_{v}$ 随剪应变 $ε_{d}$ 在三轴压缩路径下（b=0）的演化曲线如图1所示。图中体积应变以剪缩为正,剪胀为负。随着颗粒间摩擦系数μ的增大,峰值剪应力逐渐增大,且峰值后的应变软化程度愈加明显。体积响应表现为：随着μ值的增加,试样更快进入剪胀状态,且试样体积膨胀程度越明显。

图 1 应力应变曲线（b=0.0）

Figure 1. Curves of stress-strain relationship（b=0.0）

下载: 全尺寸图片幻灯片

将不同摩擦系数下等p等b真三轴数值试验得到的峰值应力状态点绘制在π平面,如图2所示,为大小不同的三角锥形。π平面内的应力演化反映出不同b值,μ值条件下,试样剪应力峰值的大小,可以看出,试样的剪应力峰值随着b值的增大而减小,随着μ值的增大而不断增大,且μ值较大时,试样对滑动摩擦系数的敏感度减弱,体现在剪应力的增加幅度逐渐减小。

图 2 不同摩擦系数下剪应力峰值状态

Figure 2. Peak states of deviatoric stress of samples under different friction coefficients

下载: 全尺寸图片幻灯片

2.2 三维强度准则

对于无黏结的颗粒材料,其抗剪强度 $τ_{f}$ = $σ_{n}$ tan( $φ_{mob}$ ),其中 $σ_{n}$ 是破坏面上的垂直压力, $φ_{mob}$ =sin^-1 $(σ_{1} - σ_{3})$ / $(σ_{1} + σ_{3})$ 是宏观内摩擦角,是颗粒材料的抗剪强度指标,反映了颗粒材料的宏观摩擦特性。Skinner等^[28]指出峰值状态的宏观内摩擦角 $φ_{max}$ 随着颗粒间摩擦系数μ的增大而增大。图3总结了多个学者得出的 $φ_{max}$ 与粒间摩擦系数μ的关系^{[19, 28-32]}。结合本文所得不同b值下的峰值内摩擦角,可以发现随着摩擦系数μ的增大,峰值内摩擦角 $φ_{max}$ 的增大变化幅度反而减小,与宏观抗剪强度的变化吻合（如图2所示）。当b=0.5时,峰值内摩擦角 $φ_{max}$ 较大,且 $φ_{max}$ 在三轴拉伸路径（b=1.0）下的值大于^{[23, 32]}三轴压缩路径（b=0.0）,与已有的DEM和物理试验研究结果基本一致。

图 3 峰值内摩擦角φ_max与滑动摩擦系数μ的关系

Figure 3. Relationship between peak internal friction angle φ_max and inter-particle friction coefficient μ

下载: 全尺寸图片幻灯片

岩土材料的强度理论中,最常采用的是Mohr-Coulomb^[33]和Drucker-Prager^[34]强度准则,但Mohr-Coulomb准则没有考虑中主应力对强度的影响,Drucker-Prager准则不能反映三轴拉伸和压缩应力状态下抗剪强度的不同。于是,国内外学者^[35-38]基于岩土材料的真三轴试验结果提出了更多的三维强度准则,其中常见的有Lade-Duncan准则^[35]和Mastsuoka- Nakai准则^[36]。为了探讨不同摩擦系数下各个强度准则的预测能力,以便于选取一个预测精准且形式简单的强度准则,定义以下指标反映强度准则的预测能力：

$error = \sqrt{\frac{1}{n} \sum_{i = 1}^{n} {(\frac{σ_{predicated}^{i} - σ_{measured}^{i}}{σ_{measured}^{i}})}^{2}}$ ,

(7)

式中, $σ_{predicated}^{i}$ 为第i个中主应力系数时的预测值, $σ_{measured}^{i}$ 为相应数值试验得到的应力破坏值。

将不同摩擦系数下各强度准则的预测误差列于表2,并绘制摩擦系数为0.05,0.2和0.5时各强度准则以及数值模拟所得的破坏应力在π平面的破坏迹线^[37-38],如图4所示。图中TC和TE分别表示三轴压缩（triaxial compression）和三轴拉伸（triaxial extension）应力路径。Mohr-Coulomb准则的破坏迹线为不等边三角形,比较靠近静水压力轴,是所有强度准则中最保守的。Drucker-Prager准则在π平面的破坏迹线是一个圆,由于其过高的估计了中主应力对强度的贡献,导致其预测的强度与数值模拟结果相差较大,但随着摩擦系数的减小,其预测误差逐渐降低。Mastsuoka-Nakai准则预测的破坏点在b=0和b=1时与Mohr-Coulomb准则重合,在0<b<1区间内,Mastsuoka-Nakai准则破坏迹线随着摩擦系数的减小逐渐外凸。Lade-Duncan准则在μ=0.5时与数值试验结果最为接近,误差值只有0.95%。当μ≤0.2时,摩擦特性对不同b值下的强度影响较小,各强度准则在π平面的破坏迹线逐渐向Drucker-Prager准则靠拢,此时各强度准则的预测能力相差不大。当0.2<μ≤0.4时,Mastsuoka-Nakai准则预测能力较好。当0.4<μ≤0.5时,Lade-Duncan准则的预测精度较高。由于岩土颗粒材料属于摩擦性颗粒材料,其细观摩擦系数大约在0.3～0.5之间^{[18, 27]},上述关于强度准则的对比,进一步说明了Mastsuoka- Nakai准则和Lade-Duncan准则在一般岩土材料中的适用性。

表 2 三维强度准则预测误差表

Table 2. Prediction errors of three-dimensional strength criteria

摩擦系数	误差/%
摩擦系数	Drucker-Prager	Mohr-Coulomb	Lade-Duncan	Mastsuoka-Nakai
0.05	6.32	5.18	3.33	2.57
0.10	10.57	6.64	7.86	6.81
0.20	13.60	7.00	8.57	5.46
0.30	15.16	8.05	5.94	1.88
0.40	15.75	10.03	3.91	1.23
0.50	16.33	12.05	0.95	5.07

下载: 导出CSV

| 显示表格

图 4 破坏面在π平面上的迹线

Figure 4. Trace of failure surface on π plane

下载: 全尺寸图片幻灯片

3. 摩擦特性与细观结构演化

3.1 强弱接触网络内的配位数分布

按照Radjai等^[21]提出的方法,用平均接触力将颗粒体系的总接触网络划分为强接触网络和弱接触网络。本文以Г_all,Г_strong,Г_weak分别表示总、强、弱接触网络。图5表明不同加载路径下峰值状态强、弱接触网络中平均法向和切向接触力比值与摩擦系数μ的关系基本不受b值的影响。由接触力比值可知,强接触网络对总接触网络平均接触力的贡献较大。随着μ值的增加,强弱接触网络平均法向接触力的比值随之增大,而平均切向接触力的比值随之减小,这是因为μ值的增加使得强接触网络中平均法向接触力明显增大,而弱接触网络中平均法向接触力增加较小。

图 5 强弱接触网络之间平均接触力的比值

Figure 5. Stress ratios of average contact stress between strong and weak contact networks

下载: 全尺寸图片幻灯片

峰值状态下,强、弱和总接触网络的配位数分布情况如图6所示。图中的配位数是指单个颗粒的接触数目,与之对应的颗粒数目是指接触网络中具有相同配位数的颗粒个数。Bratberg等^[41]研究了二维状态下颗粒材料内部的接触连通性,结果表明配位数越少的颗粒所受约束少,各向异性程度高；相应配位数越多的颗粒越稳定,各向异性程度低。在总接触网络中,随着μ值的增大,配位数众数由5减小为4,且配位数大于5或小于4的分布曲线均发生左移,试样整体的平均配位数减小,各向异性程度增加。弱接触网络中颗粒配位数分布与总接触网络的分布规律相似,但配位数众数由4减小为3。

图 6 峰值状态配位数的分布

Figure 6. Distribution of coordination numbers at peak stress

下载: 全尺寸图片幻灯片

值得注意的是,强接触网络的配位数分布与总网络差异较大,μ值的改变对强接触网络配位数为1,2的颗粒数目影响甚微,这些颗粒是构成“力链”末端及“力链”主体的颗粒群^{[23, 41]},此现象间接印证了强接触网络中“力链”的存在,而且表明在峰值状态下,试样内形成“力链”的颗粒数目不随μ值的改变而波动,或许意味着此时“力链”结构已经达到最优的颗粒分配状态。同时,配位数大于3的颗粒数目减少,说明摩擦系数越大,试样内传递力的强接触网络越稀疏,接触各向异性越明显。因而试样宏观抗剪强度增加的原因可以认为是强接触网络法向接触力的明显增大。

对比不同加载路径下强、弱和总接触网络配位数的分布曲线,发现其基本表现出与加载路径无关的特性,与刘嘉英等^[6]的研究结果一致。在既定的加载路径下,总接触网络配位数的分布中弱接触网络占主导位置,弱接触网络中配位数高的颗粒数目较多,结合配位数与颗粒材料体积变形的关系,可以推断弱接触网络对试样整体的体积变形贡献较大。

3.2 组构张量及其偏值在π平面的分布

岩土颗粒材料的宏观力学特性与加载过程中细观组构的演化规律密切相关^[42]。文章采用Satake^[43]和Oda^[44]提出的组构张量表达方式：

$Φ_{i j} = \frac{1}{N_{c}} \sum_{c = 1}^{N_{c}} n_{i} n_{j}$ ,

(8)

式中,n是沿着颗粒间接触面法线方向的单位矢量,N_c为颗粒体系中的接触总数。组构张量偏张量 ${Φ^{'}}_{i j}$ 的第二不变量可以表征颗粒体系接触方向的各向异性,用等效各向异性系数A计算如下：

$A = \sqrt{\frac{3}{2} {Φ^{'}}_{i j}^{} {Φ^{'}}_{i j}^{}} = \sqrt{3 J_{2}^{Φ}}$ 。

(9)

强接触网络的组构张量定义为

$Φ_{i j}^{s} = \frac{1}{N_{s}} \sum_{s = 1}^{N_{s}} n_{i}^{s} n_{j}^{s}$ ,

(10)

式中, $n_{i}^{s}$ 是单位法向量在第s个强接触处的第i个分量。N_s是颗粒体系中的强接触总数。强接触网络的等效各向异性系数A^s如下表示：

$A^{s} = \sqrt{\frac{3}{2} Φ_{i j}^{s^{'}} Φ_{i j}^{s^{'}}} = \sqrt{3 J_{2}^{Φ^{s}}}$ 。

(11)

类比于广义剪应力q和平均应力p的定义, $Φ_{d}$ 和 $Φ_{m}$ 可以表示为

$Φ_{d} = \frac{1}{\sqrt{2}} \sqrt{{(Φ_{11} - Φ_{22})}^{2} + {(Φ_{22} - Φ_{33})}^{2} + {(Φ_{33} - Φ_{11})}^{2}}$ ,

(12)

$Φ_{m} = (Φ_{11} + Φ_{22} + Φ_{33}) / 3$ 。

(13)

而强接触网络中 $Φ_{d}^{s}$ 和 $Φ_{m}^{s}$ 分别定义为

$Φ_{d}^{s} = \frac{1}{\sqrt{2}} \sqrt{{(Φ_{11}^{s} - Φ_{22}^{s})}^{2} + {(Φ_{22}^{s} - Φ_{33}^{s})}^{2} + {(Φ_{33}^{s} - Φ_{11}^{s})}^{2}}$ ,

(14)

$Φ_{m}^{s} = (Φ_{11}^{s} + Φ_{22}^{s} + Φ_{33}^{s}) / 3$ 。

(15)

图7是峰值状态下总接触网络和强接触网络组构张量各向异性不变量A和A^s在π平面的分布。总接触网络的组构张量不变量A在π 平面内呈现出倒三角锥形,与应力张量分布（如图2）的形状相反。而强接触网络组构张量不变量A^s在π平面内的分布表现出与宏观应力响应较好的一致性,说明强接触网络对宏观应力的贡献更为显著。在剪切变形中,不同b值和μ值的峰值应力的变化主要是由强接触网络的组构各向异性主导。Thornton等^{[29, 45]}指出总接触网络的组构各向异性程度小于强接触网络（对比图7（a）,7（b））。这是由颗粒体系的传力方式造成的,强接触承担颗粒体系大部分的力,会选择有利的接触方向以平衡外荷载的施加,使得强接触的传力网络呈现出强烈的方向各向异性。

图 7 π平面峰值状态组构张量偏张量不变量

Figure 7. Anisotropy invariants of fabric tensors on

$π$ plane for peak stress states

下载: 全尺寸图片幻灯片

应力比q/p与强接触网络组构比 $Φ_{d}^{s}$ / $Φ_{m}^{s}$ 的关系见图8。不同μ值下,q/p与 $Φ_{d}^{s}$ / $Φ_{m}^{s}$ 的比值始终约为1。Sazzad等^[24]指出不同加载路径下,应力比与强接触网络的组构比亦约为1。说明当只考虑强接触网络时,宏观应力比q/p与细观组构比 $Φ_{d}^{s}$ / $Φ_{m}^{s}$ 之间的关系始终呈线性唯一,该特性与粒间摩擦和加载路径均无关。同时也证明了颗粒体系的宏观力学特性与强接触网络组构各向异性的密切性。结合图2和图7可知,颗粒体系宏观抗剪强度随粒间摩擦的增大可能是由细观强接触网络组构各向异性的加强引起的。

图 8 不同μ值下q/p和

$Φ_{d}^{s}$ /

$Φ_{m}^{s}$ 的关系（b=0.0）

Figure 8. Relationship between q/p and

$Φ_{d}^{s}$ /

$Φ_{m}^{s}$ for different values of μ（b=0.0）

下载: 全尺寸图片幻灯片

3.3 强弱接触网络各向异性的特征及其演化

Yimsiri等^[46]在定量分析颗粒体系的各向异性时,定义了两类不同机理的各向异性：几何各向异性和力学各向异性。颗粒材料的几何各向异性可以用接触法向向量和枝向量的分布表示,此外组构张量的表示形式还可以如下所示^[43-44]：

$Φ_{i j} = \int_{Θ} E (Θ) n_{i} n_{j} d Θ$ ,

(16)

式中,Θ表示域内n的方向分布,E(Θ)是概率分布函数。一般采用E(Θ)的傅立叶展开式描述接触法向：

$E (Θ) = \frac{1}{4 π} (1 + a_{i j}^{c} n_{i} n_{j})$ ,

(17)

其中二阶各向异性张量 $a_{i j}^{c}$ 是一个偏张量,表征接触方向的各向异性程度。表达方式为

$a_{i j}^{c} = 15 / 2 {Φ^{'}}_{i j}$ ,

(18)

式中, ${Φ^{'}}_{i j}$ 是 $Φ_{i j}$ 的偏张量。

由于本文数值模拟采用的是理想圆球颗粒,枝向量几乎与接触法向同向,且颗粒尺寸分布相对较窄,故枝向量各向异性对总强度的贡献很小。力学各向异性包括法向接触力引起的法向接触力各向异性aⁿ和切向接触力引起的切向接触力各向异性a^t,其具体定义详见文献[47]。

试样初始状态是各向同性的,外荷载的施加使得试样内部结构急剧变化,因而,初期试样各向异性系数均随剪应变的增加而迅速增大,且摩擦系数μ越大,强接触网络各向异性系数增加速率越快,这是由于颗粒间摩擦越大,越能积极有效地调动强接触网络的颗粒参与外力的抵抗。从强、弱接触网络中a^c,aⁿ的演化曲线（如图9）可知,强接触网络中a^c,aⁿ值均大于弱接触网络,且均随着摩擦系数的增加而增大,说明强接触力沿加载方向的支撑作用明显,强接触网络的各向异性对总接触网络的各向异性贡献较大。

图 9 强弱接触网络中各向异性系数的演化

Figure 9. Evolution of anisotropic coefficient in strong-weak contact system

下载: 全尺寸图片幻灯片

4. 结论

（1）在颗粒体系真三轴加载条件下,随着摩擦系数的增大,宏观颗粒体系的剪应力峰值增大,进入剪胀状态更迅速,且后期剪胀程度更明显。当摩擦系数μ≤0.2时,Mohr-Coulomb、Drucker-Prager、Lade- Duncan和Matsuoka-Nakai准则的预测能力相差不大；当摩擦系数为0.2<μ≤0.5时,Lade-Duncan和Matsuoka-Nakai准则能较好地拟合不同加载路径下的数据点,更适用于一般岩土颗粒材料。

（2）粒间滑动摩擦对颗粒体系中强、弱接触网络的配位数分布影响差异明显。在三维加载条件下,各接触网络的配位数分布表现出与加载路径无关的特性。在剪应力峰值状态,弱接触网络的配位数分布受μ值的影响较大,但强接触网络中参与形成“力链”的颗粒数目基本保持不变。此外,强接触网络的组构张量偏不变量在π平面内的分布表现出与宏观应力响应较好的一致性,且颗粒体系宏观应力比q/p与细观强接触网络组构比 $Φ_{d}^{s}$ / $Φ_{m}^{s}$ 的关系呈线性唯一,与粒间摩擦和加载路径均无关,说明颗粒体系中的强接触网络在应力传递过程中发挥了重要且恒定的作用。

（3）滑动摩擦系数对强、弱接触网络的细观接触力及其组构各向异性影响显著。颗粒间摩擦系数的增加显著提高了强接触网络的法向接触力以及各接触网络的各向异性系数,其中强接触网络的法向接触力各向异性对整体力学各向异性贡献显著。

图 1 填方体边坡示意图

Figure 1. Schematic diagram of fill slope

下载: 全尺寸图片幻灯片

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

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

下载: 全尺寸图片幻灯片

图 3 样本筛分结果

Figure 3. Sample sieving results

下载: 全尺寸图片幻灯片

图 4 土石混合体试样级配曲线

Figure 4. Gradation curve of soil-rock mixture samples

下载: 全尺寸图片幻灯片

图 5 台阶状灰岩试件界面示意图

Figure 5. Interface diagram of step limestone test piece

下载: 全尺寸图片幻灯片

图 6 加工后的岩石试件

Figure 6. Processed rock samples

下载: 全尺寸图片幻灯片

图 7 剪切应力-剪切位移曲线

Figure 7. Shear stress-shear displacement curves

下载: 全尺寸图片幻灯片

图 8 块石骨架形成示意图

Figure 8. Schematic diagram of block stone skeleton formation

下载: 全尺寸图片幻灯片

图 9 块石破碎后重新形成骨架示意图

Figure 9. Schematic diagram of re-formed skeleton after rock is broken

下载: 全尺寸图片幻灯片

图 11 块石破碎形态

Figure 11. Broken forms of block stone

下载: 全尺寸图片幻灯片

图 10 块石破碎模式示意图

Figure 10. Schematic diagram of block crushing mode

下载: 全尺寸图片幻灯片

图 12 染色块石破碎与粗糙度关系

Figure 12. Relationship between broken stone blocks and roughness

下载: 全尺寸图片幻灯片

图 13 不同接触面粗糙度下抗剪强度与法向应力关系曲线

Figure 13. Relationship between shear strength and normal stress under different contact surface roughnesses

下载: 全尺寸图片幻灯片

图 14 接触面强度参数与粗糙度关系曲线

Figure 14. Relationship between contact surface strength parameters and roughness

下载: 全尺寸图片幻灯片

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

Figure 15. Schematic diagram of shear zone with different roughnesses

下载: 全尺寸图片幻灯片

表 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

表 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

表 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

表 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

参考文献(21)

[1]	张强, 汪小刚, 赵宇飞, 等. 土石混合体三维细观结构随机重构及其力学特性颗粒流数值模拟研究[J]. 岩土工程学报, 2017, 41(1): 60-69. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201901009.htm ZHANG Qiang, WANG Xiao-gang, ZHAO Yu-fei, et al. Random reconstruction of three- dimensional micro structure of soil rock mixture and numerical simulation of its mechanical properties particle flow[J]. Chinese Journal of geotechnical Engineering, 2017, 41(1): 60-69. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201901009.htm
[2]	徐文杰, 胡瑞林. 土石混合体概念、分类及意义[J]. 水文地质工程地质, 2009, 36(4): 50-56. doi: 10.3969/j.issn.1000-3665.2009.04.012 XU Wen-jie, HU Rui-lin. Conception, classification and significations of soil-rock mixture[J]. Hydrogeology and Engineering Geology, 2009, 36(4): 50-56. (in Chinese) doi: 10.3969/j.issn.1000-3665.2009.04.012
[3]	李天斌, 刘吉, 任洋, 等. 预加固高填方边坡的滑动机制：攀枝花机场12#滑坡[J]. 工程地质学报, 2012, 20(5): 723-731. doi: 10.3969/j.issn.1004-9665.2012.05.011 LI Tian-bin, LIU Ji, REN Yang, et al. Sliding mechanism of pre reinforced high fill slope: Panzhihua airport No. 12 landslide[J]. Journal of Engineering Geology, 2012, 20(5): 723-731. (in Chinese) doi: 10.3969/j.issn.1004-9665.2012.05.011
[4]	刘桂琴, 王子玉, 高振鲲. 贵州某高填方滑塌变形监测分析[J]. 人民长江, 2007(11): 148-149, 158. doi: 10.3969/j.issn.1001-4179.2007.11.057 LIU Gui-qin, WANG Zi-yu, GAO Zhen-kun. Monitoring and analysis of a high fill landslide in Guizhou[J]. People's Yangtze River, 2007(11): 148-149, 158. (in Chinese) doi: 10.3969/j.issn.1001-4179.2007.11.057
[5]	张嘎, 张建民. 粗粒土与结构接触面单调力学特性的试验研究[J]. 岩土工程学报, 2004, 26(1): 21-25. doi: 10.3321/j.issn:1000-4548.2004.01.003 ZHANG Ga, ZHANG Jian-min. Experimental study on monotonic mechanical properties of interface between coarse-grained soil and structure[J]. Chinese Journal of Geotechnical Engineering, 2004, 26(1): 21-25. (in Chinese) doi: 10.3321/j.issn:1000-4548.2004.01.003
[6]	ZHANG G, ZHANG J M. Large-scale apparatus for monotonic and cyclic soil-structure interface test[J]. Geotech Test J, 2006, 29(5): 401-408.
[7]	BORANA L, YIN J H, SINGH D N, et al. Interface behavior from suction controlled direct shear test on completely decomposed granitic soil and steel surfaces[J]. Int J Geomech, 2016, 16(6): 1-14.
[8]	BORANA L, YIN J H, SINGH D N, et al. Influence of matric suction and counterface roughness on shearing behavior of completely decomposed granitic soil and steel interface[J]. Indian Geotech J, 2016, 47(2): 150-160.
[9]	石熊, 张家生, 刘蓓, 等. 红黏土土与混凝土接触面剪切特性试验研究[J]. 中南大学学报(自然科学版), 2015, 46(5): 1826-1831. SHI Xiong, ZHANG Jia-sheng, LIU Bei, et al. Experimental study on shear characteristics of interface between red clay and concrete[J]. Journal of Central South University (Natural Science Edition), 2015, 46(5): 1826-1831. (in Chinese)
[10]	陈俊桦, 张家生, 李键. 接触面粗糙度对红黏土-混凝土接触面力学性质的影响[J]. 中南大学学报(自然科学版), 2016(47): 1682-1688. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201605030.htm CHEN Jun-hua, ZHANG Jia-sheng, LI Jian. Effect of interface roughness on mechanical properties of red clay concrete interface[J]. Journal of Central South University (Natural Science Edition), 2016(47): 1682-688. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201605030.htm
[11]	张吉顺, 华斌. 土与不同桩侧表面粗糙度接触面剪切试验研究[J]. 结构工程师, 2011, 27(3): 118-122. https://www.cnki.com.cn/Article/CJFDTOTAL-JGGC201103021.htm ZHANG Ji-shun, HUA Bin. Shear test research on contact surface between soil and different pile side surface roughness[J]. Structural Engineer, 2011, 27(3): 118-122. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JGGC201103021.htm
[12]	陈静, 李邵军, 孟凡震, 等. 三峡库区滑坡土石混合体与桩的接触面力学特性试验研究[J]. 岩石力学与工程学报, 2011, 30(增刊1): 2888-2895. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2011S1040.htm CHEN Jing, LI Shao-jun, MENG Fan-zhen, et al. Experimental study on the mechanical properties of the interface between the landslide soil rock mixture and the pile in the Three Gorges Reservoir Area[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(S1): 2888-2895. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2011S1040.htm
[13]	CEN D, HUANG D, REN F. Shear deformation and strength of the interphase between the soil-rock mixture and the benched bedrock slope surface[J]. Acta Geotechnica, 2016, 12(2): 391-413.
[14]	水利水电工程粗粒土试验规程：DL/T 5356—2006[S]. 2006. Code for Coarse-Grained Soil Test for Hydropower and Water Conservancy Engineering: DL/T 5356—2006[S]. 2006. (in Chinese)
[15]	刘新荣, 涂义亮, 王鹏, 等. 基于大型直剪试验的土石混合体颗粒破碎特征研究[J]. 岩土工程学报, 2017, 39(8): 1425-1434. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201708013.htm LIU Xin-rong, TU Yi-liang, WANG Peng, et al. Study on particle breakage characteristics of soil rock mixture based on large-scale direct shear test[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(8): 1425-1434. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201708013.htm
[16]	赵春风, 龚辉, 赵程, 等. 考虑法向应力历史的茹土一混凝土界面弹塑性分析[J]. 岩石力学与工程学报, 2012, 31(4): 848-855. ZHAO Chun-feng, GONG Hui, ZHAO Cheng, et al. Elastoplastic analysis of Ru Tu concrete interface considering normal stress history[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(4): 848-855. (in Chinese)
[17]	杨忠平, 李万坤, 胡元鑫, 等. 压实系数对粗粒土剪切特性的影响[J]. 地下空间与工程学报, 2017, 13(2): 348-356. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201702010.htm YANG Zhong-ping, LI Wan-kun, HU Yuan-xin, et al. Effect of compaction coefficient on shear properties of coarse-grained soil[J]. Journal of Underground Space and Engineering, 2017, 13(2): 348-356. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201702010.htm
[18]	陆勇, 周国庆, 夏红春, 等. 中、高压下粗粒土-结构接触面特性受结构面形貌尺度影响的试验研究[J]. 岩土力学, 2013, 34(12): 3491-3499. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201312023.htm LU Yong, ZHOU Cuo-qing, XIA Chun-hong, et al. Effect of shape scale on characteristics of coarse grained soil-structural interface under medium and high pressures[J]. Rock and Soil Mechanics, 2013, 34(12): 3491-3499. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201312023.htm
[19]	张俊峰, 王协群, 邹维列, 等. 土-格栅界面强度参数和剪切刚度试验研究[J]. 长江科学院院报, 2014, 31(3): 77-83. https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201403013.htm ZHANG Jun-feng, WANG Xie-qun, ZOU Wei-li, et al. Experimental study on strength parameters and shear stiffness of soil grid interface[J]. Journal of Changjiang Academy of Sciences, 2014, 31(3): 77-83. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201403013.htm
[20]	李广信. 高等土力学[M]. 北京: 清华大学出版社, 2004. LI Guang-xin. Advanced Soil Mechanics[M]. Beijing: Tsinghua University Press, 2004. (in Chinese)
[21]	马林. 钙质土的剪切特性试验研究[J]. 岩土力学, 2016(增刊1): 309-316. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2016S1041.htm MA Lin. Experimental study on shear properties of calcareous soil[J]. Geotechnical Mechanics, 2016(S1): 309-316. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2016S1041.htm

施引文献(28)

期刊类型引用(15)

1.	杨舒涵，漆天奇，刘嘉英. 堆石料宏细观力学特性离散元分析. 人民长江. 2024(02): 203-210 . 百度学术
2.	LI Shuai，GU Tianfeng，WANG Jiading，WANG Fei，LI Pu. Regulation effect of the grille spacing of a funnel-type grating water–sediment separation structure on the debris flow performance. Journal of Mountain Science. 2024(07): 2283-2304 . 必应学术
3.	王辉，钮新强，马刚，周伟. 干湿循环作用下堆石料宏细观力学特性的离散元模拟研究. 岩土力学. 2024(S1): 665-676 . 百度学术
4.	张江浩，冀鸿兰，杨震，李志军，刘晓民. 冻融循环下黄河堤岸砂质壤土宏细观破坏过程. 河海大学学报(自然科学版). 2024(06): 69-80 . 百度学术
5.	黄志刚，王轩，傅力，童立红. 加载速率和摩擦系数对颗粒材料系统剪切强度的影响研究. 力学季刊. 2024(04): 1032-1042 . 百度学术
6.	栗培龙，宿金菲，孙胜飞，王霄，马云飞. 多级矿料-沥青体系的颗粒特性、界面效应及迁移行为研究进展. 中国公路学报. 2023(01): 1-15 . 百度学术
7.	王伟，张志义，赵博. 煤矿井下风积沙箱式充填体侧向约束机理数值研究. 新疆大学学报(自然科学版)(中英文). 2023(03): 367-372 . 百度学术
8.	肖浩波，漆天奇，杨舒涵，周伟，刘嘉英. 椭球颗粒体系宏、细观特性的3维离散元分析. 工程科学与技术. 2023(06): 78-86 . 百度学术
9.	张革，曹玲，王成汤. 考虑各向异性影响的冻土修正线性黏结接触模型开发及应用. 岩土力学. 2023(S1): 645-654 . 百度学术
10.	郑虎，牛文清，毛无卫，黄雨. 颗粒材料双轴压缩试验的光弹测试. 同济大学学报(自然科学版). 2023(11): 1719-1724 . 百度学术
11.	王怡舒，刘斯宏，沈超敏，陈静涛. 接触摩擦对颗粒材料宏细观力学特征和能量演变规律的影响. 岩石力学与工程学报. 2022(02): 412-422 . 百度学术
12.	雷云，刘源，徐同桐，何子苗. 粒间摩擦和层厚比对二维分层颗粒系统底部响应的影响. 科学技术与工程. 2022(07): 2585-2591 . 百度学术
13.	刘嘉英，周伟，姬翔，魏纲，袁思莹，李欣骏. 基于细观拓扑结构演化的颗粒材料剪胀性分析. 力学学报. 2022(03): 707-718 . 百度学术
14.	蒋明杰，栗书亚，吉恩跃，张小勇，朱俊高. 粗粒土大型静止侧压力系数测定试验的颗粒流模拟. 科学技术与工程. 2021(25): 10867-10872 . 百度学术
15.	崔溦，魏杰，王超，王枭华，张社荣. 考虑颗粒级配和形态的颗粒柱坍塌特性离散元模拟. 岩土工程学报. 2021(12): 2230-2239 . 本站查看