颗粒破碎对堆石料填充特性缩尺效应的影响研究

徐琨; 周伟; 马刚

doi:10.11779/CJGE202006004

颗粒破碎对堆石料填充特性缩尺效应的影响研究 English Version

徐琨^1,,
周伟²,
马刚²

1.
长江勘测规划设计研究院,湖北　武汉　430010
2.
武汉大学水资源与水电工程科学国家重点实验室,湖北　武汉　430072

基金项目:

国家重点研发计划项目 2017YFC0404801

湖北省自然科学基金项目 2019CFB278

水利部土石坝破坏机理与防控技术重点实验室开放基金项目 YK319007

详细信息

作者简介:
徐琨(1989—),男,博士,工程师,主要从事堆石料缩尺效应和面板堆石坝变形控制研究。E-mail: xukun_2013@whu.edu.cn

中图分类号: TU43
计量
- 文章访问数: 364
- HTML全文浏览量: 22
- PDF下载量: 214
出版历程
- 收稿日期: 2019-08-04
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-05-31

Influence of particle breakage on scale effect of filling characteristics of rockfill material

1.
Changjiang Institute of Survey, Planning, Design and Research, Wuhan 430010, China
2.
State Key Laboratory of Water Resources and Hydropower Engineering, Wuhan University, Wuhan 430072, China

摘要

摘要: 堆石料的填充特性与其力学变形特性密切相关,颗粒破碎会改变堆石料填充特性,又是影响其缩尺效应的主要因素；以填充特性为切入点,研究堆石料缩尺效应,对深入认识堆石料复杂力学响应及其产生机理有重要意义。通过对不同颗粒破碎强度的堆石料缩尺试样开展三轴数值剪切试验,研究了颗粒破碎对堆石料填充特性缩尺效应的影响,进一步揭示了缩尺效应的产生机理。研究表明：颗粒越易破碎,缩尺效应越明显；颗粒破碎对描述填充特性的宏细观参量影响明显,主要表现为相同制样标准下,试样尺寸越大,剪切前后的孔隙率均越小,颗粒越易破碎该趋势越明显,同尺寸试样,颗粒越易破碎则剪切后的有效孔隙率越小；剪切前的有效孔隙率随试样尺寸的增大而增大,与剪切后的规律相反；试验中,大尺寸试样较小尺寸试样产生了更多的力学不稳定颗粒,是造成堆石料缩尺效应的原因之一。
- 堆石料 /
- 颗粒破碎 /
- 离散元法 /
- 填充特性 /
- 缩尺效应
Abstract: The filling characteristics of rockfill materials are closely related to their mechanical properties. They can be changed by particle breakage, which is one of the main factors of the scale effect of rockfill materials. It is of great significance to further understand the complex mechanical response and mechanism of scale effect from the perspective of filling properties. A series of numerical triaxial compression tests considering different particle breakage strengths are taken to study the influence of particle breakage on the scale effect of filling characteristics of rockfill materials, and the mechanism of the scale effect is further revealed. It is found that the lower the particle breakage strength, the more obvious the scale effect. The particle breakage has obvious influence on the macro- and micro-parameters of filling characteristics. The porosity decreases with the increasing sample size before and after compression tests, which is more significant with the lower particle breakage strength. The effective porosity of samples with the same size after compression tests decreases with the lower particle breakage strength. The effective porosity of samples before compression tests increases with the increasing sample size, which is opposite to the trend after compression tests. During the process of compression tests, the larger samples produce more mechanical unstable particles than smaller ones, which is one of the reasons for the scale effect of rockfill materials.
- rockfill material /
- particle breakage /
- discrete element method /
- filling characteristic /
- scale effect

HTML全文

0. 引言

天然土体颗粒间通常具有胶结物质,比如碳酸钙、微生物、黏粒、甲烷水合物等,称为结构性土。结构性土可看成是一种特殊的胶结颗粒材料,颗粒间胶结物质的存在使结构性土的宏观力学性质不同于重塑土^[1-2]。为了分析结构性土工程问题,需要建立合理的结构性土本构模型。在临界状态土力学基础上,相关学者引入状态变量,分别建立了边界面模型^[3]、下负荷面模型^[4]、统一硬化模型^[5]等重塑土本构模型。在加荷过程中,结构性土颗粒间的胶结逐渐破坏并退出承担应力,土体向重塑土转化。因此,在重塑土本构模型基础上,通过引入表征结构性的参量及其演化规律建立结构性土本构模型是一条可行的研究途径。不少学者提出了结构性土本构模型,在模型中假定了不同的结构性表征参量。比如：①通过比较结构性土屈服面大小 ${\bar{p}}_{c}$ 和重塑土屈服面大小 $p_{c}$ 的不同,定义 $Δ p_{c} = {\bar{p}}_{c} - p_{c}$ ^[6]或 $R^{*} = p_{c} /$ ${\bar{p}}_{c}$ ^[7]为结构性表征参量；②通过引入胶结组分（相对完整状态）和摩擦组分（完全调整状态）应力分担参量（扰动参量、结构破损参量）D^[8]或B₀^[9-10]表达结构性衰退程度。由于胶结破损演化难以通过试验测得,结构性表征参量的演化通常需要假设,缺乏微观试验结果验证。

离散单元法^[11]模拟结构性土弥补了室内微观测试在获取胶结破损演化方面的不足,已广泛应用于结构性土体宏微观力学性质研究。Jiang等^[10]通过开展结构性土的二维离散元模拟,对3种结构性土本构模型结构性参量演化假设的合理性进行了验证,基于离散元结果提出了新的胶结破损演化式,并应用到本构模型研究^[12],不足之处在于胶结破损演化式中的参数受加载路径影响。通过PFC软件植入完整胶结接触模型（考虑了颗粒间抗转动和抗扭转能力及胶结尺寸对刚度和强度的影响）建立了结构性土三维离散元试样,开展了侧限压缩、等向压缩、等应力比压缩以及常规三轴和真三轴试验的离散元数值模拟,分析了结构破损参量B₀随等效塑性应变的演化规律,并提出了应力路径相关性低的结构破损参量B_σ。

1. 胶结接触模型

离散元模拟采用软硬复合胶结接触模型框架^[13],考虑了颗粒间抗转动和抗扭转能力^[14]及胶结尺寸对刚度和强度的影响^[15]。对于通过胶结间接传力/力矩的接触（“厚胶结”）直接采用胶结接触公式计算接触力/力矩；对于颗粒直接接触并有胶结传力/力矩的接触（“薄胶结”）,按并联模式通过颗粒接触力/力矩和胶结接触力/力矩之和计算接触力/力矩。

1.1 颗粒接触

假设颗粒间通过一圆形截面接触,半径为

$R_{c} = β R$ ,

(1)

式中, $β$ 为抗转动系数, $R = 2 R_{1} R_{2} / (R_{1} + R_{2})$ , $R_{1}$ , $R_{2}$ 为两颗粒半径。

接触间力/力矩的传递通过下列公式计算：

$F_{n}^{l} = {\begin{cases} k_{n} δ_{n} - F_{a} (δ_{n} \geq 0) \\ 0 (δ_{n} < 0) \end{cases}$ ,

(2)

$F_{s}^{l} = {(F_{s}^{l})}_{o} - k_{s} Δ δ_{s}$ ,

(3)

$M_{r}^{l} = {(M_{r}^{l})}_{o} - k_{r} Δ θ_{r}$ ,

(4)

$M_{t}^{l} = {(M_{t}^{l})}_{o} - k_{t} Δ θ_{t}$ 。

(5)

式中 $F_{a}$ 为粒间引力,比如毛细力、范德华力、电子引力等； $F_{n}^{l}$ , $F_{s}^{l}$ , $M_{r}^{l}$ , $M_{t}^{l}$ 分别为颗粒法向力、切向力、弯矩和扭矩； $k_{n}$ , $k_{s}$ ,k_r,k_t分别为颗粒法向、切向、抗弯、抗扭刚度, $k_{r} = 0.25 k_{n} R_{c}^{2}$ , $k_{t} = 0.5 k_{s} R_{c}^{2}$ ； $δ_{n}$ 为颗粒重叠量, $Δ δ_{s}$ , $Δ θ_{r}$ 和 $Δ θ_{t}$ 分别为剪切位移、转动角和扭转角增量。

颗粒接触的法切向刚度可按PFC软件手册^[16]按下式计算：

$k_{n} = A E^{*} / L$ ,

(6)

$k_{s} = k_{n} / κ^{*}$ ,

(7)

式中, $A = π R^{2}$ , $E^{*}$ 是颗粒接触的等效模量, $κ^{*}$ 是颗粒接触的法切向刚度比, $L$ 为两颗粒半径之和。

接触抗剪强度、抗弯强度和抗扭强度分别为^[14-15]

$‖ F_{s}^{l} ‖ \leq μ (F_{n} + F_{a})$ ,

(8)

$‖ M_{r}^{l} ‖ \leq 0.25 ζ_{c} R_{c} (F_{n} + F_{a})$ ,

(9)

$M_{t}^{l} \leq 0.65 μ R_{c} (F_{n} + F_{a})$ ,

(10)

式中, $μ$ 是摩擦系数, $ζ_{c}$ 取2.1。颗粒接触力和力矩达到以上强度后颗粒将进入滑移、转动和扭转状态。

1.2 胶结计算

假定胶结物质呈圆形截面连接两颗粒,则胶结法向力 $F_{n}^{b}$ ,切向力 $F_{s}^{b}$ ,弯矩 $M_{r}^{b}$ 和扭矩 $M_{t}^{b}$ 分别为^[16]

$F_{n}^{b} = {(F_{n}^{b})}_{o} - k_{n}^{b} Δ δ_{n}^{}$ ,

(11)

$F_{s}^{b} = {(F_{s}^{b})}_{o} - k_{s}^{b} Δ δ_{s}^{}$ ,

(12)

$M_{r}^{b} = {(M_{r}^{b})}_{o} - k_{r}^{b} Δ θ_{r}^{}$ ,

(13)

$M_{t}^{b} = {(M_{t}^{b})}_{o} - k_{t}^{b} Δ θ_{t}$ 。

(14)

式中 $k_{n}^{b}$ , $k_{s}^{b}$ , $k_{r}^{b}$ , $k_{t}^{b}$ 分别为胶结法向、切向、抗弯和抗扭刚度, $k_{t}^{b} = 0.5 k_{s}^{b} R_{b}^{2}$ , $k_{r}^{b} = 0.25 k_{n}^{b} R_{b}^{2}$ ； $Δ δ_{n}^{}$ 是法向位移增量； $R_{b} = λ_{b} R$ 为胶结半径, $λ_{b}$ 为胶结半径系数。

为了考虑胶结厚度和宽度对胶结刚度的影响,用下式计算胶结法向刚度（图1）^[15]：

图 1 胶结刚度计算示意图

Figure 1. Schematic diagram for calculation of bond stiffness

下载: 全尺寸图片幻灯片

$k_{n}^{b} = 1 / (δ_{A 1} + δ_{A 2} + δ_{B 1} + δ_{B 2} + δ_{C 1} + δ_{C 2})$ ,

(15)

式中, $δ_{A 1}$ , $δ_{B 1}$ , $δ_{C 1}$ 分别为颗粒1的A1_,B1和C1三段的线柔度,

$δ_{A 1} = \frac{R_{1} ξ_{E}}{π E_{1} R_{b}^{2}}$ ,

(16)

$δ_{B 1} = {\begin{matrix} \frac{1}{2 π (E_{1} - E_{2}) R_{1} Λ_{E}} (\ln (\frac{Λ_{E} + 1}{Λ_{E} - 1}) - \ln (\frac{Λ_{E} + ξ_{E}}{Λ_{E} - ξ_{E}})) (E_{1} > E_{2}) \\ \frac{R_{1} - R_{1} ξ_{E}}{π E_{1} R_{b}^{2}} (E_{1} = E_{2}) \end{matrix}$ ,

(17)

$δ_{C 1} = δ_{C 2} = \frac{h_{0}^{b} / 2}{π E_{2} R_{b}^{2}}$ ,

(18)

式中, $E_{1}$ , $E_{2}$ 分别为骨架材料和胶结材料的等效变形模量, $R_{1}$ 为颗粒1的半径, $h_{0}^{b}$ 为胶结厚度,同理可计算 $δ_{A 2}$ , $δ_{B 2}$ 。在模型使用中, $E_{1} = {\bar{E}}^{*}$ 为胶结等效模量, $E_{2} = η_{E} {\bar{E}}^{*}_{}$ , $η_{E}$ 为模量折减系数。 $ξ_{E}$ , $Λ_{E}$ 计算如下：

$ξ_{E} = \sqrt{1 - {(R_{b} / R_{1})}^{2}}$ ,

(19)

$Λ_{E} = \sqrt{1 + \frac{E_{2}}{E_{1} - E_{2}} {(R_{b} / R_{1})}^{2}}$ 。

(20)

胶结切向刚度根据法切向刚度比 $κ_{b}$ 计算：

$k_{s}^{b} = k_{n}^{b} / κ_{b}$ 。

(21)

采用考虑胶结厚度和宽度的胶结强度计算方法^[15,17]。胶结抗压和抗拉强度计算如下：

$R_{nc}^{b} = π R_{b}^{2} σ_{c}^{b} \frac{c_{1}}{λ_{b}^{c_{3}}} \exp (\frac{c_{2} h_{0}^{b}}{2 R_{b}})$ ,

(22)

$R_{nt}^{b} = π R_{b}^{2} σ_{t}^{b}$ ,

(23)

式中, $σ_{c}^{b}$ , $σ_{t}^{b}$ 分别为胶结材料抗压和抗拉强度； $c_{1} = 1.45$ , $c_{2} = - 0.7$ 和 $c_{3} = 1 / 6$ 是反映胶结半径和厚度对抗压强度影响的参数（定值）。

在一定的法向力下（未达到胶结抗压和抗拉强度）,胶结可能承受剪切力、弯矩和扭矩,复合荷载作用下胶结强度符合椭球包面：

${(‖ F_{s}^{b} ‖ / R_{s}^{b})}^{2} + {(‖ M_{r}^{b} ‖ / R_{r}^{b})}^{2} + {(M_{t}^{b} / R_{t}^{b})}^{2} = 1$ ,

(24)

式中, $R_{s}^{b}$ , $R_{r}^{b}$ , $R_{t}^{b}$ 分别为胶结抗剪、抗弯和抗扭强度,

$R_{s}^{b} = S_{s}^{b} (R_{nc}^{b} + R_{nt}^{b})$ ,

(25)

$R_{r}^{b} = S_{r}^{b} (R_{nc}^{b} + R_{nt}^{b}) R_{b}$ ,

(26)

$R_{t}^{b} = S_{t}^{b} (R_{nc}^{b} + R_{nt}^{b}) R_{b}$ ,

(27)

式中, $S_{s}^{b}$ , $S_{r}^{b}$ , $S_{t}^{b}$ 为胶结抗剪、抗弯和抗扭包面的形状,

$S_{i}^{b} = m_{i} f_{n}^{b} {[\ln (1 / f_{n}^{b})]}^{3 / 5}$ ,

(28)

式中,下标i是s,r和t, $m_{s} = 0.5$ , $m_{r} = 0.3$ 和 $m_{t} = 0.36$ 为参数（定值）, $f_{n}^{b} = (F_{n}^{b} + R_{nt}^{b}) / (R_{nc}^{b} + R_{nt}^{b})$ 为正则化胶结法向力。

胶结的存在对结构性土体的宏微观力学性质有很大影响。胶结达到强度包面后会发生胶结破坏,胶结破坏后胶结部分不再承担力和力矩,接触力和力矩退化到残余值（颗粒接触部分）。

2. 离散元模拟步骤

2.1 离散元试样制备

以饱和结构性黄土^[15]为基准制备离散元试样和选取接触参数,为提高研究的适用性开展了不同胶结材料强度离散元试样的数值分析,下文仅以一种胶结材料强度为例进行分析。采用分层欠压法^[18]分5层制备离散元试样,颗粒数目42180个,如图2（a）所示。结构性土离散元试样制备采用的颗粒级配曲线如图2（b）所示。试样生成后在12.5 kPa竖向压力下压缩稳定以反映原位K₀应力状态。试样稳定后若两个颗粒之间的间距小于 $g_{c} R_{b}$ 则生成胶结, $g_{c}$ 为胶结临界厚度系数。

图 2 离散元试样及其粒径级配

Figure 2. DEM sample and grain-size distribution curve

下载: 全尺寸图片幻灯片

颗粒间引力作用能使试样以更松散状态稳定赋存^[19],试样制备过程施加弱吸引力以制备松散试样^[15]：

$F_{a} = σ_{a} d_{50}^{2}$ ,

(29)

式中,d₅₀为中值粒径, $σ_{a}$ 为弱吸引力系数,根据典型黄土的最大稳定孔隙比和饱和重塑土黏聚力,通过试算取4 kPa。颗粒与墙体间仅承担法向力作用,忽略摩擦。颗粒密度=2710 kg/m³,颗粒局部阻尼系数=0.7。颗粒之间的接触参数：颗粒等效模量=800 MPa,法切向刚度比=1.5,摩擦系数=0.5,抗转动系数=0.21,胶结等效模量=200 MPa,胶结模量折减系数=0.2,胶结法切向刚度比=2,胶结半径系数=0.35,胶结材料抗压强度=25MPa,胶结材料抗拉强度=2.5MPa,胶结临界厚度系数=0.1。

2.2 不同应力路径离散元模拟步骤

为了研究不同应力路径结构性土的胶结破损演化,对上述结构性土离散元试样进行了分级加载的侧限压缩试验、等向压缩试验和等应力比压缩试验,以及应变加载的常规三轴试验和等p真三轴试验。

对分级加载压缩试验,上下墙体从12.5 kPa开始分级施加荷载,每次施加的荷载约为上级荷载的 $\sqrt{2}$ 倍；侧墙伺服应力目标通过应力比要求计算（侧限压缩试验使4个侧墙固定）。

对三轴试验,首先将离散元试样等向压缩到不同的围压,然后对侧墙进行伺服控制,使上下两墙相向运动进行剪切试验。在常规三轴试验中,侧墙的伺服目标为围压。在真三轴试验中通过平均应力和中主应力系数b不变的条件计算侧墙伺服目标应力。

3. 离散元模拟结果

为验证离散元数值分析结果的有效性,在研究结构性土的胶结破损演化规律之前,需确保离散元试样能够反映结构性土的主要压缩和剪切特性。考虑到本文针对一般结构性土体（不同胶结强度的胶结颗粒材料）进行破损规律数学建模,因此,下面主要进行定性对比,未对某一特定的结构性土进行标定。

图3给出了结构性土离散元试样在不同压缩路径下的压缩曲线。图3中 $η = q / p$ 为应力比, $CSR = σ_{3} / σ_{1}$ 为主应力比。随着平均应力的增加,离散元试样压缩线先缓慢下降而后逐渐转变为快速下降,存在结构屈服应力,离散元试样能够定性再现结构性土体的主要压缩特性^[20]。离散元试样体积变形特征与真实结构性土定量上差别较大,主要是因为颗粒形状、颗粒破碎机制等微观结构的差别。

图 3 结构性土压缩性质

Figure 3. Compression properties of structured soils

下载: 全尺寸图片幻灯片

图4给出了结构性土离散元试样常规三轴试验应力–应变关系和孔隙比–应变关系。离散元分析能够定性再现结构性土体三轴剪切性质^[2,21]。试样在剪切初始段模量很大,近似弹性段。随着围压的增加,试样抗剪强度增加,孔隙比减小,试样从倾向于软化和剪胀向倾向于硬化和剪缩发展。离散元反映结构性土剪切破坏特性的效果优于反映压缩变形特性。

图 4 结构性土常规三轴剪切性质

Figure 4. Conventional triaxial test results of structured soils

下载: 全尺寸图片幻灯片

图5给出了结构性土离散元试样真三轴试验应力-应变关系和孔隙比–应变关系（以p=300 kPa为例）。可见离散元模拟结果能够反映室内试验中中主应力系数b（ $b = (σ_{2} - σ_{3}) / (σ_{1} - σ_{3})$ ）的影响^[22-23]。同样平均应力下随着b的增加,破坏偏应力下降。b=0时试样倾向于剪缩,随着b的增加试样剪缩程度减弱。天然土试样的变形模量相对离散元试样和人工制备结构性土试样较低,这可能与试样扰动有关。

图 5 结构性土真三轴剪切性质

Figure 5. True triaxial test results of structured soils

下载: 全尺寸图片幻灯片

4. 胶结破损演化规律

4.1 胶结颗粒材料微观力学

土体可以看成一种颗粒集合体,重塑土（无黏结）可以看成是一种散粒体材料而结构性土可以看成是一种胶结颗粒材料。对于理想的散粒体材料,土中的应力通过接触力进行传递,而土中的应变由颗粒间的相对变形产生。基于此,学者通过研究土颗粒之间的接触关系和微观变形机理,建立了理想散粒体材料宏观应力和应变与微观接触力和变形之间的基本方程^[24-25]。Jiang等^[10]在此基础上建立了胶结颗粒材料的宏微观理论框架。

对于颗粒材料,荷载通过颗粒之间的接触力进行传递,因此颗粒材料的代表性单元的平均应力张量可根据下式求得

${\bar{σ}}_{i j}^{} = \frac{1}{V^{}} \sum_{k = 1}^{N^{}} F_{i}^{k} l_{j}^{k}$ ,

(30)

式中,V为代表性单元的体积,N为代表性单元内颗粒间接触点的总数, $F_{i}^{k}$ 为第k个接触点处接触力的分量, $l_{j}^{k}$ 为第k个接触点处方向向量的分量。

对于结构性土体,根据颗粒间接触点是否存在胶结,土体代表性单元平均应力张量 ${\bar{σ}}_{i j}^{}$ 可分解为胶结接触分担应力 ${\bar{σ}}_{i j}^{in}$ 与无胶结接触分担应力 ${\bar{σ}}_{i j}^{d}$ 之和,

$\begin{array}{l} {\bar{σ}}_{i j} = \frac{1}{V^{}} (\sum_{k = 1}^{N^{in}} F_{i}^{k} l_{j}^{k} + \sum_{k = 1}^{N^{d}} F_{i}^{k} l_{j}^{k}) \\ = \frac{1}{V^{}} （ V^{in} {\bar{σ}}_{i j}^{in} + V^{d} {\bar{σ}}_{i j}^{d}) = (1 - λ) {\bar{σ}}_{i j}^{in} + λ {\bar{σ}}_{i j}^{d} \end{array}$ ,

(31)

式中, $N^{in}$ , $N^{d}$ 分别为代表性单元内胶结接触点和无胶结接触点的数目,且 $N = N^{in} + N^{d}$ 。 $V^{in}$ , $V^{d}$ 分别为代表性单元内胶结颗粒所占体积与无胶结颗粒所占体积（含相应的孔隙体积）,且 $V = V^{in} + V^{d}$ ； $λ = V^{d} / V$ 为体积破损率,即代表性单元内无胶结组分所占体积与代表性单元总体积的比值； ${\bar{σ}}_{i j}^{in}$ 和 ${\bar{σ}}_{i j}^{d}$ 可由下式计算：

$\begin{array}{l} {\bar{σ}}_{i j}^{in} = \frac{1}{V^{in}} \sum_{k = 1}^{N^{in}} F_{i}^{k} l_{j}^{k}, \\ {\bar{σ}}_{i j}^{d} = \frac{1}{V^{d}} \sum_{k = 1}^{N^{d}} F_{i}^{k} l_{j}^{k} 。 \end{array}}$

(32)

为了研究结构性土体的结构性演化,需要在本构模型中引入（或蕴含）胶结破损参量。体积破损率就是一个简单直观的胶结破损参量,考虑到胶结组分和无胶结组分的体积难以划分,可用无胶结接触数量与总接触数量比值近似表示,

$λ = \frac{V^{d}}{V} \approx \frac{N^{d}}{N}$ 。

(33)

无胶结组分分担的应力与总应力的比值也可以作为胶结破损参量,

$B_{0} = \frac{λ {\bar{σ}}_{i j}^{d}}{{\bar{σ}}_{i j}} = \frac{\sum_{k = 1}^{N^{d}} F_{i}^{k} l_{j}^{k}}{\sum_{k = 1}^{N} F_{i}^{k} l_{j}^{k}}$ 。

(34)

胶结破损参量表征胶结破损程度或结构性衰退程度。当胶结破损参量为0时,土体结构性完好,土体应力全部由胶结组分承担；随着胶结破损参量的增加,土体应力由胶结组分和无胶结组分共同分担（胶结组分承担应力逐渐变小）；当胶结破损参量为1时,土体结构性完全破坏,土体应力全部由无胶结组分承担。

胶结破损参量 $λ$ , $B_{0}$ 难以通过室内试验获得,离散单元法可以用来研究 $λ$ , $B_{0}$ 的演化。尽管 $λ$ , $B_{0}$ 参量的大小有一定不同,但二者演化规律是相似的,因此本文以破损参量 $B_{0}$ 为例进行研究。

4.2 胶结破损演化规律

随着塑性应变的发展,胶结逐渐破损。由于弹性应变较小,塑性应变采用如下简单方法计算。离散元试样的弹性能 $E_{k}$ 为颗粒弹性能 $E_{k}^{p}$ 和胶结弹性能 $E_{k}^{b}$ 之和：

$E_{k}^{p} = \frac{1}{2} (\frac{{(F_{n}^{l})}^{2}}{k_{n}} + \frac{{‖ F_{s}^{l} ‖}^{2}}{k_{s}} + \frac{{‖ M_{r}^{l} ‖}^{2}}{k_{r}} + \frac{{(M_{t}^{l})}^{2}}{k_{t}})$ ,

(35)

$E_{k}^{b} = \frac{1}{2} (\frac{{(F_{n}^{b})}^{2}}{k_{n}^{b}} + \frac{{‖ F_{s}^{b} ‖}^{2}}{k_{s}^{b}} + \frac{{‖ M_{r}^{b} ‖}^{2}}{k_{r}^{b}} + \frac{{(M_{t}^{b})}^{2}}{k_{t}^{b}})$ ,

(36)

$E_{k} = E_{k}^{p} + E_{k}^{b} = \frac{1}{2} (p ε_{v}^{e} + q ε_{s}^{e}) V_{DEM}$ ,

(37)

$K = p / ε_{v}^{e}$ ,

(38)

$G = q / ε_{s}^{e}$ ,

(39)

式中, $V_{DEM}$ 为离散元试样的体积, $ε_{v}^{e}$ 为弹性体应变, $ε_{s}^{e}$ 为弹性偏应变,假设体积模量K和剪切模量G相等,联立以上各式可求得 $ε_{v}^{e}$ , $ε_{s}^{e}$ ,从而求出塑性体应变和塑性偏应变。

图6给出了破损参量随塑性体应变和塑性偏应变的变化规律。对于三轴剪切试验,塑性体应变较小且变化复杂,而对于等应力比压缩试验（应力比较小时）,塑性偏应变较小,二者均不适合单独用于描述试样破损参量演化。因此,胶结的破损既与塑性体应变有关,又与塑性偏应变有关。

图 6 不同应力路径下破损参量B₀随塑性体应变和塑性偏应变的演化规律

Figure 6. Evolution of degradation variable B₀ versus plastic volumetric and deviator strain for DEM sample under different stress paths

下载: 全尺寸图片幻灯片

为了考虑两种塑性应变对胶结破损参量的影响,采用下式简单等效塑性应变研究胶结破损演化。

$ε_{d}^{p} = \sqrt{{(ε_{v}^{p})}^{2} + {(ε_{s}^{p})}^{2}}$ ,

(40)

式中, $ε_{v}^{p}$ 为塑性体应变, $ε_{s}^{p}$ 为塑性偏应变。

图7给出结构性土离散元试样胶结破损参量 $B_{0}$ 随等效塑性应变的变化规律。从图7中可见,对不同应力路径,破损参量 $B_{0}$ 随等效塑性应变的增加先较快增加,而后缓慢增加。但不同应力路径的破损参量 $B_{0}$ 演化在数值上有较大不同。对不同应力比压缩试验,随着应力比的增加, $B_{0}$ 表征的胶结破损速率加快；对不同中主应力系数b的真三轴试验,随着b的增加, $B_{0}$ 表征的胶结破损速率减小；对不同围压常规三轴试验,围压也对 $B_{0}$ 表征的胶结破损速率影响显著。

图 7 不同应力路径下破损参量B₀随等效塑性应变的演化规律

Figure 7. Evolution of degradation variable B₀ versus equivalent plastic strain for DEM sample under different stress paths

下载: 全尺寸图片幻灯片

考虑到胶结破损跟胶结强度相关,在宏观上,胶结破损速率跟目前应力状态和强度包线之间的相对关系相关,结构参数 $B_{0}$ 破损的速率随应力比的增加而增加。因此,将应力比和强度包线引入胶结破损参量 $B_{0}$ ,通过应力比与峰值应力比的比值 $η / M_{f}$ 调整胶结破损速率,构建一个新的胶结破损参量如下式,从而蕴含应力路径对胶结破损的影响,

$B_{σ} = \frac{2 B_{0}}{(1 + η^{2} / M_{f}^{2})}$ ,

(41)

式中, $M_{f} = q_{f} / p$ 为峰值应力比,本文通过常规三轴试验结果采用莫尔–库仑准则得到峰值偏应力 $q_{f} =$ $0.9275 p + 123.8 kPa$ 。

图8给出了结构性土离散元试样胶结破损参量 $B_{σ}$ 随等效塑性应变的破损规律。对不同应力路径,胶结破损参量 $B_{σ}$ 随等效塑性应变逐渐增加,但破损速率逐渐降低。除等向压缩试验（ $η = 0$ ）外,试验应力比、围压和中主应力系数对 $B_{σ}$ 表征的胶结破损演化影响很小。由此可见,胶结破损参量 $B_{σ}$ 应力路径相关性明显降低,可作为状态变量引入结构性土体二元介质类本构模型（比如扰动状态概念模型、损伤模型等）表征胶结破损情况,使模型的结构破损规律更为合理。通过数据分析,建议采用指数函数拟合破损参量 $B_{σ}$ 与等效塑性应变的关系,

图 8 不同应力路径下破损参量B_σ随等效塑性应变的演化规律

Figure 8. Evolution of degradation variable B_σversus equivalent plastic strain for DEM sample under different stress paths

下载: 全尺寸图片幻灯片

$B_{σ} = 1 - \exp [- c_{a} {(ε_{d}^{p})}^{c_{b}}]$ ,

(42)

式中, $c_{a}$ , $c_{b}$ 为胶结破损相关的参数,可通过侧限压缩或常规三轴试验确定。

5. 结论

通过将相对完备的三维胶结接触模型植入粒间接触,反映结构性土粒间胶结的关键特征,建立了结构性土离散元模型；开展了侧限压缩、等向压缩、等应力比压缩以及常规三轴、真三轴等不同加载路径下的离散元数值分析；再现了结构性土体的主要宏观力学性质（尤其是结构性土的破坏特性）；分析了胶结破损参量 $B_{0}$ , $B_{σ}$ 随等效塑性应变的演化规律。

（1）胶结破损参量（ $λ$ , $B_{0}$ , $B_{σ}$ ）表征胶结破损程度或结构性衰退程度。当胶结破损参量为0时,土体结构性完好,土体应力全部由胶结组分承担；当胶结破损参量为1时,土体结构性完全破坏,土体应力全部由无胶结组分承担。

（2）胶结破损参量 $B_{0}$ 的演化受应力比、围压和中主应力系数等的影响,具有应力路径相关性。将应力比和强度包线引入胶结破损参量 $B_{0}$ ,构建了一个新的胶结破损参量 $B_{σ}$ ,应力路径相关性明显降低。在本构模型应用中,胶结破损参量 $B_{σ}$ 随等效塑性应变的演化规律可用指数函数表示,并通过侧限压缩或常规三轴试验确定相关参数。

不同应力路径下 $B_{σ}$ 的演化仍然具有一定的离散性,今后将考虑改进等效塑性应变的定义和强度包面公式。因为胶结破损规律 $B_{σ}$ 基于离散元结果提出,适用于理想胶结颗粒材料（理想结构性土体）,对真实结构性土体的适用性需要进一步验证。

图 1 碎片替换模式：3种尺寸共14个颗粒^[32]

Figure 1. Fragment replacement mode: 14 particles with 3 particle sizes^[32]

下载: 全尺寸图片幻灯片

图 2 原级配曲线和缩尺后的级配曲线

Figure 2. Prototype gradation curve and scaled gradation curves

下载: 全尺寸图片幻灯片

图 3 不同缩尺级配的初始试样示意图

Figure 3. Illustration of initial samples with different particle sizes

下载: 全尺寸图片幻灯片

图 4 各组缩尺试样应力–应变曲线

Figure 4. Stress-strain behaviors of numerical tests

下载: 全尺寸图片幻灯片

图 5 各试样的颗粒破碎率 $B_{m}$

Figure 5. Particle breakage factor of numerical tests

下载: 全尺寸图片幻灯片

图 6 孔隙率演化曲线

Figure 6. Evolution curves of porosity

下载: 全尺寸图片幻灯片

图 7 试验围压1.2 MPa时试样剪切前后的孔隙率

Figure 7. Porosity of samples before and after compression tests

下载: 全尺寸图片幻灯片

图 8 有效孔隙率演化曲线

Figure 8. Evolution curves of effective porosity

下载: 全尺寸图片幻灯片

图 9 试验围压1.2 MPa时试样剪切前后的有效孔隙率

Figure 9. Effective porosity of samples before and after compression tests (CP = 1.2 MPa)

下载: 全尺寸图片幻灯片

图 10 试样有效配位数演化曲线

Figure 10. Evolution curves of effective coordination number

下载: 全尺寸图片幻灯片

图 11 初始试样的配位数频率分布图

Figure 11. Frequency distribution of coordination number of initial samples

下载: 全尺寸图片幻灯片

图 12 试样剪切前后配位数频率分布图

Figure 12. Frequency distribution of coordination number of samples before and after compression tests

下载: 全尺寸图片幻灯片

图 13 各组缩尺试样配位数频率分布图

Figure 13. Frequency distribution of coordination number of samples

下载: 全尺寸图片幻灯片

表 1 离散元模拟参数

Table 1 Input parameters for DEM simulation

类别	参数	取值
颗粒参数	密度/(kg·m^-3)	2650
	摩擦系数 $μ$	0.5
	剪切模量G/GPa	28
	泊松比 $ν$	0.15
颗粒破碎准则参数	特征粒径 $d_{0}$ /mm	28
	特征破碎强度 $σ_{0}$ /MPa	18.2,14.0
	韦伯模量m	2.65

下载: 导出CSV

参考文献(43)

[1]	郦能惠. 高混凝土面板堆石坝新技术[M]. 北京: 中国水利水电出版社, 2007. LI Neng-hui. Recent Technology for High Concrete Face Rockfill Dams[M]. Beijing: China water and Power Press, 2007. (in Chinese)
[2]	ALAN BISHOP W, HENKEL D J. The Measurement of Soil Properties in the Triaxial Test[R]. London: Edward Arnold Ltd., 1948.
[3]	MARSAL R J. Large-scale testing of rockfill materials[J]. Journal of the Soil Mechanics and Foundations Division, 1967, 93(2): 27-43. doi: 10.1061/JSFEAQ.0000958
[4]	MARACHI N D, CHAN C K, SEED H B. Evaluation of properties of rockfill mechanicals[J]. Journal of Soil Mechanics and Foundation Engineering, ASCE, 1972, 98(1): 95-114.
[5]	李翀, 何昌荣, 王琛, 等. 粗粒料大型三轴试验的尺寸效应研究[J]. 岩土力学, 2008, 29(增刊1): 563-566. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2008S1113.htm LI Chong, HE Chang-rong, WANG Chen, et al. Study of scale effect of large-scale triaxial test of coarse-grained meterials[J]. Rock and Soil Mechanics, 2008, 29(S1): 563-566. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2008S1113.htm
[6]	凌华, 殷宗泽, 朱俊高, 等. 堆石料强度的缩尺效应试验研究[J]. 河海大学学报：自然科学版, 2011, 39(5): 540-544. https://www.cnki.com.cn/Article/CJFDTOTAL-HHDX201904019.htm LIN Hua, YIN Zong-ze, ZHU Jun-gao, et al. Experimental study of scale effect on strength of rockfill materials[J]. Journal of Hohai University(Natural Sciences), 2011, 39(5): 540-544. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HHDX201904019.htm
[7]	朱俊高, 刘忠, 翁厚洋, 等. 试样尺寸对粗粒土强度及变形试验影响研究[J]. 四川大学学报(工程科学版), 2012, 44(6): 92-96. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201206016.htm ZHU Jun-gao, LIU Zhong, WEN Hou-yang, et al. Study on effect of specimen size upon strength and deformation behaviour of coarse-grained soil in triaxial test[J]. Journal of Sichuan University (Engineering Science Edition), 2012, 44(6): 92-96. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201206016.htm
[8]	武利强, 朱晟, 章晓桦, 等. 粗粒料试验缩尺效应的分析研究[J]. 岩土力学, 2016, 37(8): 2187-2197. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201608009.htm WU Li-qiang, ZHU Sheng, ZHANG Xiao-hua, et al. Analysis of scale effect of coarse-grained materials[J]. Rock and Soil Mechanics, 2016, 37(8): 2187-2197. (in chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201608009.htm
[9]	VARADARAJAN A, SHARMA K G, VENKATACHALAM K, et al. Testing and modeling two rockfill materials[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2003, 129(3): 206-218. doi: 10.1061/(ASCE)1090-0241(2003)129:3(206)
[10]	郦能惠, 朱铁, 米占宽. 小浪底坝过渡料的强度与变形特性及缩尺效应[J]. 水电能源科学, 2001, 19(2): 39-42. doi: 10.3969/j.issn.1000-7709.2001.02.012 LI Neng-hui, ZHU Tie, MI Zhan-kuan. Strength and deformation properties of transition zone material of Xiaolangdi dam and scale effect[J]. Hydroelectric Energy, 2001, 19(2): 39-42. (in Chinese) doi: 10.3969/j.issn.1000-7709.2001.02.012
[11]	王继庄. 粗粒料的变形特性和缩尺效应[J]. 岩土工程学报, 1994, 16(4): 89-95. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC404.011.htm WANG Ji-zhuang. The deformation characteristic and scaling effect of coarse-grained material[J]. Chinese Journal of Geotechnical Engineering, 1994, 16(4): 89-95. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC404.011.htm
[12]	WEI K M, ZHU S, YU X H. Influence of the scale effect on the mechanical parameters of coarse-grained soils[J]. Iranian Journal of Science and Technology, Transactions of Civil Engineering, 2014, 38(C1): 75-84.
[13]	花俊杰, 周伟, 常晓林, 等. 堆石体应力变形的尺寸效应研究[J]. 岩石力学与工程学报, 2010, 29(2): 328-335. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201002016.htm HUA Jun-jie, ZHOU Wei, CHANG Xiao-lin, et al. Study of scale effect on stress and deformation of rockfill[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(2): 328-335. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201002016.htm
[14]	高莲士, 蔡昌光, 朱家启. 堆石料现场侧限压缩试验解耦K-G模型参数分析方法及在面板坝中的应用[J]. 水力发电学报, 2006, 25(6): 26-33. doi: 10.3969/j.issn.1003-1243.2006.06.005 GAO Lian-shi, CAI Chang-guang, ZHU Jia-qi. An analysis method for uncoupled K-G model parameters in site confined compression test of rock-fill materials and its application on CFRD[J]. Journal of Hydroelectric Engineering, 2006, 25(6): 26-33. (in Chinese) doi: 10.3969/j.issn.1003-1243.2006.06.005
[15]	李凤鸣, 卞富宗. 两种粗粒土的比较试验[J]. 勘察科学技术, 1991(2): 25-29. https://www.cnki.com.cn/Article/CJFDTOTAL-KCKX199102006.htm LI Feng-ming, BIAN Fu-zong. Comparative test of two kinds of coarse-grained soils[J]. Site Investigation Science and Technology, 1991(2): 25-29. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KCKX199102006.htm
[16]	董槐三, 尹承瑶. 天生桥一级水电站面板堆石坝筑坝材料性质研究[J]. 红水河, 1996, 15(4): 7-12. https://www.cnki.com.cn/Article/CJFDTOTAL-HSHZ604.001.htm DONG Huai-san, YIN Cheng-yao. Study of construction material properties for TSQ-I concrete faced rockfill dam[J]. Hong Shui He, 1996, 15(4): 7-12. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HSHZ604.001.htm
[17]	朱俊高, 翁厚洋, 吴晓铭, 等. 粗粒料级配缩尺后压实密度试验研究[J]. 岩土力学, 2010, 31(8): 2394-2398. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201008008.htm ZHU Jun-gao, WENG Hou-yang, WU Xiao-ming, et al. Experimental study of compact density of scaled coarse- grained soil[J]. Rock and Soil Mechanics, 2010, 31(8): 2394-2398. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201008008.htm
[18]	左永振, 张伟, 潘家军, 等. 粗粒料级配缩尺方法对其最大干密度的影响研究[J]. 岩土力学, 2015, 36(增刊1): 417-422. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S1073.htm ZUO Yong-zhen, ZHANG Wei, PAN Jia-jun, et al. Effects of gradation scale method on maximum dry density of coarse-grained soil[J]. Rock and Soil Mechanics, 2015, 36(S1): 417-422. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S1073.htm
[19]	朱晟, 王永明, 翁厚洋. 粗粒筑坝材料密实度的缩尺效应研究[J]. 岩石力学与工程学报, 2011, 30(2): 348-357. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201102020.htm ZHU Sheng, WANG Yong-ming, WENG Hou-yang. Study of scale effect of density of coarse-grained dam materials[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(2): 348-357. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201102020.htm
[20]	XIAO Y, MENG M, DAOUADJI A, et al. Effect of particle size on crushing and deformation behaviors of rockfill materials[J]. Geoscience Frontiers, 2019. doi: 10.1016/j.gsf.2018.10.010
[21]	孔宪京, 刘京茂, 邹德高. 堆石料尺寸效应研究面临的问题及多尺度三轴试验平台[J]. 岩土工程学报, 2016, 38(11): 1941-1947. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201611003.htm KONG Xian-jing, LIU Jing-mao, ZOU De-gao. Scale effect of rockfill and multiple-scale triaxial test platform[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(11): 1941-1947. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201611003.htm
[22]	OVALLE C, FROSSARD E, DANO C, et al. The effect of size on the strength of coarse rock aggregates and large rockfill samples through experimental data[J]. Acta Mechanica, 2014, 225(8): 2199-2216.
[23]	ALONSO E E, TAPIAS M, GILI J. Scale effects in rockfill behaviour[J]. Géotechnique Letters, 2012, 2(3): 155-160.
[24]	马刚, 周伟, 常晓林, 等. 堆石料缩尺效应的细观机制研究[J]. 岩石力学与工程学报, 2012, 31(12): 2473-2482. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201212012.htm MA Gang, ZHOU Wei, CHANG Xiao-lin, et al. Mesoscopic mechanism study of scale effects of rockfill[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(12): 2473-2482. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201212012.htm
[25]	徐琨, 周伟, 马刚, 等. 基于离散元法的颗粒破碎模拟研究进展[J]. 岩土工程学报, 2018, 40(5): 880-889. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201805016.htm XU Kun, ZHOU Wei, MA Gang, et al. Review of particle breakage simulation based on DEM[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(5): 880-889. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201805016.htm
[26]	MCDOWELL G R, DE BONO J P. On the micro mechanics of one-dimensional normal compression[J]. Géotechnique, 2013, 63(11): 895-908.
[27]	DE BONO J P, MCDOWELL G R. DEM of triaxial tests on crushable sand[J]. Granular Matter, 2014, 16(4): 551-562.
[28]	WEIBULL W. A statistical theory of the strength of materials[J]. Proceedings of the American Mathematical Society, 1939, 151(5): 1034-1034.
[29]	BAZANT Z P. Size effect in blunt fracture: concrete, rock, metal[J]. J Eng Mech, 1984, 110(4): 518-535.
[30]	MCDOWELL G R, AMON A. The application of weibull statistics to the fracture of soil particles[J]. Soils and Foundations, 2000, 40(5): 133-141.
[31]	LI H, MCDOWELL G, LOWNDES I. Discrete element modelling of a rock cone crusher[J]. Powder Technology, 2014, 263: 151-158.
[32]	CIANTIA M, ARROYO Alvarez de Toledo M, CALVETTI F, et al. An approach to enhance efficiency of DEM modelling of soils with crushable grains[J]. Géotechnique, 2015, 65(2): 91-110.
[33]	LI G, LIU Y J, DANO C, et al. Grading-dependent behavior of granular materials: from discrete to continuous modeling[J]. Journal of Engineering Mechanics, 2014, 141(6): 04014172.
[34]	Itasca Consulting Group Inc. Particle Flow Code in 3 Dimensions (PFC3D) version 5.0[CP]. Minneapolis: Itasca Consulting Inc, 2014: 2199-2216.
[35]	TYLER S W, WHEATCRAFT S W. Fractal scaling of soil particle-size distributions: analysis and limitations[J]. Soil Science Society of America Journal, 1992, 56(2): 362-369.
[36]	ZHOU W, XU K, MA G, et al. Effects of particle size ratio on the macro-and microscopic behaviors of binary mixtures at the maximum packing efficiency state[J]. Granular Matter, 2016, 18(4): 81.
[37]	王永明, 朱晟, 任金明, 等. 筑坝粗粒料力学特性的缩尺效应研究[J]. 岩土力学, 2013, 34(6): 1799-1806. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201306041.htm WANG Yong-ming, ZHU Sheng, REN Jin-ming, et al. Research on scale effect of coarse-grained materials[J]. Rock and Soil Mechanics, 2013, 34(6): 1799-1806. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201306041.htm
[38]	NG T T, ZHOU W, CHANG X L. Effect of particle shape and fine content on the behavior of binary mixture[J]. Journal of Engineering Mechanics, 2016, 143(1): C4016008.
[39]	SITHARAM T G, NIMBKAR M S. Micromechanical modelling of granular materials: effect of particle size and gradation[J]. Geotechnical and Geological Engineering, 2000, 18(2): 91-117.
[40]	THORNTON C. Numerical simulations of deviatoric shear deformation of granular media[J]. Géotechnique, 2000, 50: 43-53.
[41]	ROTHENBURG L, KRUYT N P. Critical state and evolution of coordination number in simulated granular materials[J]. International Journal of Solids and Structures, 2004, 41(21): 5763-5774.
[42]	ZAMPONI F. Mathematical physics: packings close and loose[J]. Nature, 2008, 453: 606-607.
[43]	GU X, HUANG M, QIAN J. DEM investigation on the evolution of microstructure in granular soils under shearing[J]. Granular Matter, 2014, 16(1): 91-106.

施引文献(29)

期刊类型引用(19)

1.	沈书拯，江巍，欧阳晔，谭亿虹，吴剑，陈勇. 基于三维连续-离散耦合模型的门架式双排抗滑桩受力分析. 防灾减灾工程学报. 2025(01): 84-94 . 百度学术
2.	王洪涛，李建华，殷允腾，张红军，宣兆腾，赵明珠. 承压水作用下岩质地层基坑抗突涌安全厚度研究. 建筑结构学报. 2024(06): 213-223 . 百度学术
3.	马勤. 超挖对软土基坑工程稳定性的影响分析. 现代矿业. 2024(04): 121-123 . 百度学术
4.	张瑞程. 新型斜撑-排桩基坑支护体系工程应用分析. 施工技术(中英文). 2024(12): 162-166 . 百度学术
5.	王强，宋歌，廖利，刘国伟，蔡崇德. 全地下空间枢纽深基坑支护结构变形特征控制. 铁路技术创新. 2024(03): 178-185 . 百度学术
6.	赵升峰，汪敏营，邬喜春，李明东. 不同深度基坑中桩加支撑支护特性研究. 江苏建筑. 2024(06): 112-115+129 . 百度学术
7.	赵玉成，代兆宇，夏瑞萌，娄海成. 雄安新区深基坑双排桩支护的实践与探讨. 粉煤灰综合利用. 2024(06): 79-83+95 . 百度学术
8.	蔡红岩. 供水泵站基坑排桩支护局部失效引发连续破坏机理研究. 陕西水利. 2023(01): 100-102 . 百度学术
9.	曾诚. 基坑局部设置双排桩桩长桩径对连续变形控制效果影响研究. 中国新技术新产品. 2023(07): 81-85 . 百度学术
10.	梁文辉. 深基坑双排桩式围护结构支护施工技术研究. 江西建材. 2023(09): 244-245+248 . 百度学术
11.	黄礼明. 某邻近既有建筑的异形软土基坑变形控制. 福建建筑. 2023(11): 84-88 . 百度学术
12.	徐琳，黄程翔，邵根才，杨璞，李宝山. 软土地区基坑超挖对围护结构变形的影响. 建筑安全. 2023(12): 33-36 . 百度学术
13.	魏焕卫，李传斌，种记鑫，郑晓. 局部超挖对内支撑结构的施工效应. 计算机辅助工程. 2023(04): 56-65 . 百度学术
14.	常强. 双排桩治理边坡的数值模拟研究. 黑龙江水利科技. 2022(04): 37-39+179 . 百度学术
15.	郑刚，程雪松，周海祚，张天奇，于晓旋，刁钰，王若展，衣凡，张文彬，郭伟. 岩土与地下工程结构韧性评价与控制. 土木工程学报. 2022(07): 1-38 . 百度学术
16.	唐小军，袁仕贵. 铁路桥台基坑开挖支护方案及稳定性计算. 中国水运. 2022(08): 146-148 . 百度学术
17.	唐小军，袁仕贵. 铁路桥台基坑开挖支护方案及稳定性计算. 中国水运. 2022(15): 146-148 . 百度学术
18.	孙永梅，程海涛，米春荣，孙贝. 土体超挖对深基坑变形影响分析. 山西建筑. 2021(19): 61-63 . 百度学术
19.	曾锦秀. 基于上限法的边坡双排抗滑桩变形计算方法. 福建工程学院学报. 2021(06): 538-544 . 百度学术