三轴试样绝对体积测量方法研究

蔡阳; 李林; 路毅

doi:10.11779/CJGE202112017

三轴试样绝对体积测量方法研究 English Version

南京林业大学土木工程学院,江苏　南京　210037

详细信息

作者简介:
蔡阳(1996— ),男,硕士研究生,主要从事土工测试技术研究。E-mail：caiyang2019@njfu.edu.cn

中图分类号: TU411
计量
- 文章访问数: 230
- HTML全文浏览量: 46
- PDF下载量: 144
出版历程
- 收稿日期: 2021-03-01
- 网络出版日期: 2022-11-30
- 刊出日期: 2021-11-30

Measuring absolute volume of triaxial soil specimens

College of Civil Engineering, Nanjing Forestry University, Nanjing 210037, China

摘要

摘要: 三轴试验过程中的土样体积是用于分析土体应力应变特性所需的一个基本参数。近年来,基于图像分析的方法越来越广泛地应用于土样在三轴试验过程中的体积测量。然而,基于图像分析方法测量所得结果仅仅是试样的相对体积变化,并非试样绝对体积。为了测量三轴试验过程中试样的绝对体积,综合运用摄影测量原理、光线追踪、和最小二乘优化技术,测定了覆盖试样表面点云的三维位置；通过对点云进行三角网格划分、端部截断、插值、并扣除橡皮膜体积计算得到土样在三轴试验过程中的绝对体积；同时,还针对性地开发出了用于实现上述分析步骤以及相关测量结果后处理的配套软件,GeoTri3D。为验证该绝对体积测量方法的可行性,开展了一个标准钢柱试验、一个饱和砂土和三个非饱和粉土试样的三轴试验。钢柱与饱和砂土的三轴试验结果表明该方法能够成功实现三轴试样在试验过程中任意时刻的绝对体积测量；此外,通过对粉土试样表面三维点云的进一步分析所得的土样应变云图能够真实再现三轴试验过程中的土样全局变形过程。
- GeoTri3D /
- 摄影测量 /
- 绝对体积 /
- 三轴试样 /
- 全局变形
Abstract: The soil volume in triaxial tests is an essential parameter to characterize stress and strain behaviors of soils. In recent years, more and more image-based methods have been developed for the triaxial tests on deformation measurement of soils. However, the measured results are only the relative volume change of specimens, not their absolute volume. In order to measure the absolute soil volume in triaxial tests, the photogrammetry principle, optical ray tracing and the least square optimization are used to measure the three-dimensional position of point clouds on the soil surface. The absolute volume of the soil specimens during triaxial tests is then obtained through triangulation, end cut and interpolation based on the obtained point clouds. Meantime, a software package, GeoTri3D, is developed to implement the above-mentioned steps and post processing of the measured results. In order to validate the proposed method, triaxial tests are conducted on a standard steel cylinder, a saturated sand specimen, and three soil specimens. The results obtained from the validation tests indicate that the proposed method can successfully and accurately capture the absolute volume of the soil specimens during triaxial tests at different time of interests. In addition, based on the analysis on the triangular meshes, the strain contour plots can be generated, which facilitates the visualization of the deforming process of soils during triaxial tests.
- GeoTri3D /
- photogrammetry /
- absolute volume /
- triaxial specimen /
- full-field deformation

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. Setup of coordinate system and measurement system

下载: 全尺寸图片幻灯片

图 2 钢柱试验

Figure 2. Tests on a steel cylinder

下载: 全尺寸图片幻灯片

图 3 光线追踪过程

Figure 3. Tracing process of optical ray

下载: 全尺寸图片幻灯片

图 4 三角网格和插值

Figure 4. Generation of triangular mesh and interpolation

下载: 全尺寸图片幻灯片

图 5 饱和砂土样的三轴剪切试验结果

Figure 5. Triaxial test results on saturated sand

下载: 全尺寸图片幻灯片

图 6 土样在不同时刻的三角网格

Figure 6. Evolution of triangular mesh

下载: 全尺寸图片幻灯片

图 7 不同围压下的土样试验结果 (5, 20, 60 kPa)

Figure 7. Results of specimens tested under different confining pressures (5, 20, 60 kPa)

下载: 全尺寸图片幻灯片

图 8 轴向应变

Figure 8. Axial strain

下载: 全尺寸图片幻灯片

图 9 试样剪切带与轴向应变图

Figure 9. Failure plane and axial strain

下载: 全尺寸图片幻灯片

表 1 试样参数

Table 1 Information of specimens

试样	围压/ kPa	含水率/%	直径/mm	高度/mm
钢柱	0		61.77	124.99
饱和砂	50		71.18	118.43
粉土A	5	9.46	61.82	124.95
粉土B	20	9.58	61.82	125.05
粉土C	60	9.72	61.79	124.98

下载: 导出CSV

表 2 相机校核结果

Table 2 Calibrated results of camera

参数	理想化前	理想化后
f/mm	55.7901	55.7901
M/(pixel)	5568	5568
N/(pixel)	3712	3712
F_x /mm	23.9809	24.5611
F_y /mm	16.0000	16.3726
P_x /mm	11.8583	12.2806
P_y /mm	7.9973	8.1863
K₁/(10^-5)	6.410	0
K₂/(10^-8)	0	0
P₁/(10^-6)	2.167	0
P₂/(10^-6)	-2.283	0

下载: 导出CSV

表 3 钢柱体积

Table 3 Volumes of steel cylinder

插值倍数	空气中体积/ cm³	压力室内体积/ cm³
0	368.793	368.753
2	376.372	376.331
4	376.983	376.941
6	377.151	377.110
8	377.220	377.179
16	377.297	377.255

下载: 导出CSV

表 4 饱和砂土的顶面方程参数

Table 4 Information of top end plane of saturated sand specimens

位移/mm	a	b	c
0	-0.010	-0.010	138.43
2	-0.006	-0.016	136.63
4	-0.002	-0.018	134.66
6	0.002	-0.021	132.64
8	0.006	-0.023	130.65
10	0.012	-0.025	128.62
13	0.019	-0.029	125.62
16	0.022	-0.032	122.59
20	0.028	-0.031	118.61

下载: 导出CSV

参考文献(23)

[1]	JI E, ZHU J, CHEN S, et al. Analytical solution and experimental study of membrane penetration in triaxial test[J]. Geomechanics and Engineering, 2017, 13(6): 1027-1044.
[2]	BISHOP A W. The Experimental Study of Partly Saturated Soil in the Triaxial Apparatus[C]//Proc 5th Int Conf Soil Mech Found Eng, 1961, Paris.
[3]	ROMERO E, FACIO J A, LLORET A, et al. A New Suction and Temperature Controlled Triaxial Apparatus[C]//Proceedings of the 14th International Conference on Soil Mechanics and Foundation Engineering. 1997, Hamburg.
[4]	KLOTZ E, COOP M. On the identification of critical state lines for sands[J]. Geotechnical Testing Journal, 2002, 25(3): 289-302.
[5]	LAUDAHN A, SOSNA K, BOHAC J. A simple method for air volume change measurement in triaxial tests[J]. Geotechnical Testing Journal, 2005, 28(3): 313-318.
[6]	MACARI E, PARKER J, COSTES N. Measurement of volume changes in triaxial tests using digital imaging techniques[J]. Geotechnical Testing Journal, 1997, 20(1): 103-109. doi: 10.1520/GTJ11424J
[7]	LIN H, PENUMADU D. Strain localization in combined axial-torsional testing on Kaolin clay[J]. Journal of Engineering Mechanics, 2006, 132(5): 555-564. doi: 10.1061/(ASCE)0733-9399(2006)132:5(555)
[8]	GACHET P, GEISER F, LALOUI L, et al. Automated digital image processing for volume change measurement in triaxial cells[J]. Geotechnical Testing Journal, 2007, 30(2): 98-103.
[9]	UCHAIPICHAT A, KHALILI N, ZARGARBASHI S, et al. A temperature controlled triaxial apparatus for testing unsaturated soils[J]. Geotechnical Testing Journal, 2011, 34(5): 424-432.
[10]	邵龙潭, 郭晓霞, 刘港, 等. 数字图像测量技术在土工三轴试验中的应用[J]. 岩土力学, 2015, 36(增刊1): 669-84. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S1117.htm SHAO Long-tan, GUO Xiao-xia, LIU Gang, et al. Application of digital image processing technique to measuring specimen deformation in triaxial test[J]. Rock and Soil Mechanics, 2015, 36(S1): 669-684. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S1117.htm
[11]	邵龙潭, 刘潇, 郭晓霞, 等. 土工三轴试验试样全表面变形测量的实现[J]. 岩土工程学报, 2012, 34(3): 409-415. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201203008.htm SHAO Long-tan, LIU Xiao, GUO Xiao-xia, et al. Whole surface deformation measurement of triaxial soil specimen based on digital image processing[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(3): 409-415. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201203008.htm
[12]	BHANDARI A, POWRIE W, HARKNESS R. A digital image-based deformation measurement system for triaxial tests[J]. Geotechnical Testing Journal, 2012, 35(2): 209-226.
[13]	WANG P P, GUO X X, SANG Y, et al. Measurement of local and volumetric deformation in geotechnical triaxial testing using 3D-digital image correlation and a subpixel edge detection algorithm[J]. Acta Geotechnica, 2020, 15(10): 2891-2904.
[14]	SALAZAR S E, COFFMAN R A. Consideration of internal board camera optics for triaxial testing applications[J]. Geotechnical Testing Journal, 2015, 38(1): 20140163.
[15]	SALAZAR S E, BARNES A, COFFMAN R A. Development of an internal camera-based volume determination system for triaxial testing[J]. Geotechnical Testing Journal, 2015, 38(4): 20140249.
[16]	ZHANG X, LI L, CHEN G, et al. A photogrammetry-based method to measure total and local volume changes of unsaturated soils during triaxial testing[J]. Acta Geotechnica, 2015, 10(1): 55-82.
[17]	LI L, ZHANG X, CHEN G, et al. Measuring unsaturated soil deformations during triaxial testing using a photogrammetry-based method[J]. Canadian Geotechnical Journal, 2016, 53(3): 472-489.
[18]	LI L, ZHANG X. A new triaxial testing system for unsaturated soil characterization[J]. Geotechnical Testing Journal, 2015, 38(6): 20140201.
[19]	LI L, ZHANG X. Factors influencing the accuracy of the photogrammetry-based deformation measurement method[J]. Acta Geotechnica, 2019, 14(2): 559-574.
[20]	LI L, LU Y, CAI Y, et al. A calibration technique to improve accuracy of the photogrammetry-based deformation measurement method for triaxial testing[J]. Acta Geotechnica, 2021, 16(4): 1053-1060.
[21]	FAYEK S, XIA X, LI L, et al. A photogrammetry-based approach to determine the absolute volume of soil specimen during triaxial testing[J]. Journal of the Transportation Research Board, 2020, 2674(8): 206-218.
[22]	BALDI G, NOVA R. Membrane penetration effects in triaxial testing[J]. Journal of Geotechnical Engineering, 1984, 110(3): 403-420.
[23]	详细操作步骤参见文献[EB/OL]. https://v.youku.com/v_show/id_XNDk5NjkwNDAyMA==.html?spm=a2h0c.8166622.PhoneSokuUgc_1.dtitle.

施引文献(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 . 本站查看