三轴排水剪切下钙质砂的颗粒破碎特性

张季如; 华晨; 罗明星; 张弼文

doi:10.11779/CJGE202009003

三轴排水剪切下钙质砂的颗粒破碎特性 English Version

武汉理工大学土木工程与建筑学院,湖北　武汉　430070

基金项目:

国家自然科学基金面上项目 41272334

详细信息

作者简介:
张季如(1964—),男,教授,博士生导师,主要从事岩土工程方面的教学和科研工作。E-mail：zhangjr@whut.edu.cn

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

Behavior of particle breakage in calcareous sand during drained triaxial shearing

School of Civil Engineering and Architecture, Wuhan University of Technology, Wuhan 430070, China

摘要

摘要: 颗粒破碎是影响粒状土的变形和强度机理的重要因素。为了研究钙质砂在剪切过程中的颗粒破碎特性及其对变形和强度性质的影响,对3种不同初始分布的钙质砂进行了不同围压下的三轴排水剪切试验。结果显示：初始分形的粒径分布在三轴剪切过程中始终保持着较为严格的分形特性,该现象与各粒组中的破碎颗粒主要向相邻的下一级粒组中迁移的机制有关。钙质砂的应力–应变特性与围压大小和初始粒径分布有关,围压越低,初始粒径分布越不均匀,钙质砂的剪胀效应越显著。随着围压的增大,钙质砂的剪胀倾向减少,并逐渐过渡到剪缩状态。钙质砂的破碎率随剪切过程中的应力和应变的增长而增大,其峰值内摩擦角随着破碎率的增大而降低,最后趋于定值。用非线性的指数函数来描述峰值内摩擦角与破碎率的相关关系,揭示了颗粒破碎对钙质砂抗剪强度的影响规律。
- 钙质砂 /
- 三轴排水剪试验 /
- 颗粒破碎 /
- 粒径分布 /
- 抗剪强度
Abstract: Particle breakage is an important factor affecting the deformation and strength mechanism of granular soil. In order to study the particle breakage characteristics of calcareous sand during shearing and their influences on the deformation and shear strength, several drained triaxial shear tests are performed on three kinds of calcareous sand with different initial particle-size distributions under varied confining pressures. The results show that the initial fractal particle-size distribution maintains a relatively strict fractal characteristic in the triaxial shear process. This phenomenon can be well explained by the behavior that the broken particles in each particle fraction are mainly added to the adjacent next-level particle fraction. The stress-strain characteristics of calcareous sand are related to confining pressure and distribution of initial particle size. The lower the confining pressure, the more uneven the distribution of initial particle size, the more significant the dilatancy effect of calcareous sand. As the confining pressure increases, the dilatancy tendency of calcareous sand decreases and gradually transits to the shrinkage state. The breakage index of calcareous sand increases with the increase of stress and strain in the shearing process. The peak angle of internal friction decreases with the increase of breakage index, and finally tends toward a fixed value. A non-linear exponential function is used to describe the correlation between the peak angle of internal friction and the breakage index, and it represents the effect of particle crushing on the shear strength of calcareous sand.
- calcareous sand /
- drained triaxial shear test /
- particle breakage /
- particle-size distribution /
- shear strength

HTML全文

0. 引言

柱状节理常见于玄武岩岩体,节理面独特的空间结构分布使得岩体强度表现出强各向异性。已有学者针对柱状节理岩体开展了现场原位试验和大量的相似材料室内试验^[1-5],验证了柱状节理岩体的强度各向异性并总结了材料强度参数的各向异性分布规律。但目前对于柱状节理岩体各向异性本构研究的报道还不多,而已有的一些各向同性模型则不能很好的描述柱状节理岩体的强度各向异性特性。

柱状节理岩体为岩块和结构面的综合体,其表征单元体较室内试验可研究的尺度要大许多,很难通过室内试验来研究柱状节理岩体的变形和强度特征。而原位试验虽然可以在一定程度上反映节理岩体的工程力学特性,但受限于场地、设备、成本等因素,一些物理力学参数难以获取。由于柱状节理岩体具有结构周期性,可以应用数值均匀化方法开展多尺度分析计算,如多重有限元（FE²）或跨尺度有限元离散元耦合（FEM/DEM）。该类方法避免了复杂本构方程的建立,适用于满足均匀化条件的复杂细观结构材料。但多尺度方法计算量巨大,难以直接应用于大型水电工程的计算。因此建立可以准确描述柱状节理岩体力学特性的本构模型具有重要意义。柱状节理岩体显著的各向异性是本构的建模的难点。微结构张量作为材料微观结构的一种度量,可用来描述岩体介质中孔隙或裂纹的空间分布。本文在传统力学模型的基础上引入微结构张量,建立柱状节理岩体各向异性本构模型。

Kanatani^[6-7]最早系统研究了材料空间分布与微结构张量的关系,初步建立了空间分布函数的表达式。Pietruszczak等^[8]通过组合强度参数的空间分布,把经典各向同性准则拓展为各向异性,建立了基于微结构张量的各向异性屈服准则。Gao等^[9]将微结构张量参数引入屈服函数,提出了一种具有交叉各向异性的岩土材料的广义破坏准则,该准则在描述大范围材料的强度随加载方向的变化时具有通用性和鲁棒性。钟世英等^[10]采用多组节理面局部坐标系与整体坐标系的投影关系定义了材料空间的微结构张量表达,将其引入莫尔库仑强度准则,得到了多组节理岩体的各向异性强度准则。

然而现阶段基于微结构张量的各向异性本构模型尚存在不足,或以页岩等层状岩石为研究对象^[8-9],模型不能较好反映柱状节理岩体物理力学试验特性,或模型较为复杂,参数不易确定^[10],实际应用较为困难。

为此,本文提出了一种既保留了数学的严谨性又具有一定实用性的方法。首先,结合微结构张量理论推导了柱状节理岩体各向异性弹塑性本构积分公式,然后基于FLAC^3D对自定义本构进行了二次开发,并结合柱状节理相似结构材料三轴岩石力学试验对本文提出的自定义本构模型进行了验证。应用所开发各向异性弹塑性本构模型对白鹤滩水电站坝基开挖工程进行了数值模拟分析,数值模拟结果与坝基岩体开挖工程监测结果基本相符。研究结果对白鹤滩水电工程建设具有一定应用价值。

1. 各向异性弹塑性本构

1.1 弹性刚度矩阵及空间旋转变换

理想弹性体的本构关系可表示为

$ε = C σ$ 。

(1)

柱状节理岩体（如图1）可简化为横观各向同性材料,由广义胡克定律,在材料空间的局部坐标系下其柔度矩阵可表示为

图 1 全局坐标系中的柱状节理材料

Figure 1. Cylindrical joint materials in global coordinate system

下载: 全尺寸图片幻灯片

$\begin{matrix} C = [\begin{matrix} \frac{1}{E_{1}} & - \frac{ν_{12}}{E_{1}} & - \frac{ν_{13}}{E_{3}} & 0 & 0 & 0 \\ - \frac{ν_{12}}{E_{1}} & \frac{1}{E_{1}} & - \frac{ν_{13}}{E_{3}} & 0 & 0 & 0 \\ - \frac{ν_{13}}{E_{3}} & - \frac{ν_{13}}{E_{3}} & \frac{1}{E_{3}} & 0 & 0 & 0 \\ 0 & 0 & 0 & \frac{1}{G_{13}} & 0 & 0 \\ 0 & 0 & 0 & 0 & \frac{1}{G_{13}} & 0 \\ 0 & 0 & 0 & 0 & 0 & \frac{2 (1 + ν_{12})}{E_{1}} \end{matrix}] \\ = {[D^{e}]}^{- 1} \end{matrix}$ ,

(2)

式中, $E_{1}$ 为平行于各向同性面的弹性模量, $E_{3}$ 为垂直于各向同性面的弹性模量, $G_{13}$ 为垂直平面的剪切模量, $ν_{12}$ 为平行面内的泊松比, $ν_{13}$ 为垂直方向的泊松比, $D^{e}$ 为弹性刚度矩阵。

全局坐标系下的柔度矩阵 $C^{'}$ 可由局部坐标系下的柔度矩阵 $C$ 旋转得到

$C^{'} = Q^{T} C Q$ ,

(3)

式中, $Q$ 为一个包含局部坐标轴与全局坐标轴夹角 $α$ 的二阶张量,可表示为

$Q = [\begin{matrix} \cos^{2} α & \sin^{2} α & 0 & \frac{\sin 2 α}{- 2} & 0 & 0 \\ \sin^{2} α & \cos^{2} α & 0 & \frac{\sin 2 α}{2} & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 & 0 \\ \sin 2 α & - \sin 2 α & 0 & \cos 2 α & 0 & 0 \\ 0 & 0 & 0 & 0 & \cos α & - \sin α \\ 0 & 0 & 0 & 0 & \sin α & \cos α \end{matrix}]$ 。

(4)

对于横观各向同性材料,需要给出两个夹角,分别为材料的倾向 $α$ 和倾角 $β$ 。后续推导若不加说明,默认在材料空间中的局部坐标系展开。

1.2 各向异性屈服函数及塑性势函数

定义微结构张量 $a_{i j}$ 主轴与材料空间的局部坐标系同轴,三维空间中微结构张量 $a_{i j}$ 可由单位向量 $e_{i}$ 表示为

$\begin{matrix} a_{i j} = a_{1} e_{i}^{(1)} e_{j}^{(1)} + a_{2} e_{i}^{(2)} e_{j}^{(2)} + e_{i}^{(3)} e_{j}^{(3)} \\ = a_{1} m_{i j}^{(1)} + a_{2} m_{i j}^{(2)} + a_{3} m_{i j}^{(3)} \end{matrix}$ ,.

(5)

式中, $a_{i}$ 为微结构张量的主值, $m_{i j}$ 为结构取向张量。

定义广义加载单位向量 $l_{i}$ 为

$l_{i} = \frac{L_{i}}{\sqrt{L_{k} L_{k}}} = \frac{L_{i}}{\sqrt{σ_{i j} σ_{i j}}} = \frac{L_{i}}{\sqrt{tr σ^{2}}}$ ,

(6)

$L_{i} = \sum_{x = 1}^{3} L_{x} e_{i}^{(x)}$ ,

(7)

$L_{x} = \sqrt{σ_{1 x}^{2} + σ_{2 x}^{2} + σ_{3 x}^{2}}$ 。

(8)

考虑微结构张量 $a_{i j}$ 在 $l_{i}$ 方向上的投影,即

$\begin{matrix} η = a_{i j} l_{i} l_{j} = a_{1} \frac{tr(m^{(1)} σ^{2})}{tr σ^{2}} + a_{2} \frac{tr(m^{(2)} σ^{2})}{tr σ^{2}} + \\ a_{3} \frac{tr(m^{(3)} σ^{2})}{tr σ^{2}} \end{matrix}$ 。

(9)

式（9）所定义参数 $η$ 是一个无量纲的标量,综合考虑了材料结构面取向对加载空间的影响,反映了材料强度的方向依赖性,因此可用来描述材料强度的各向异性空间分布。式（9）也可表示为偏量的形式：

$η = \hat{η} (1 + A_{i j} l_{i} l_{j})$ ,

(10)

$\hat{η} = a_{k k} / 3, A_{i j} = (a_{i j} - \hat{η} δ_{i j}) / \hat{η}$ 。

(11)

式（10）最早由Kanatani^[6-7]提出。更一般的,该式可表示为

$\begin{matrix} η = \hat{η} (1 + A_{i j} l_{i} l_{j} + A_{i j k l} l_{i} l_{j} l_{k} l_{l} + \dots) \\ = \hat{η} (1 + A_{i j} l_{i} l_{j} + b_{1} {(A_{i j} l_{i} l_{j})}^{2} + \dots) \end{matrix}$ 。

(12)

本文采用一种修正的莫尔库仑准则^[11]来描述岩石柱体的塑性行为,并考虑各向异性和塑性硬化：

$\begin{matrix} f = \sqrt{{(q K (θ))}^{2} + {(e_{f} c \cos ϕ)}^{2}} + \\ η_{f} h (ξ) (p \sin ϕ - c \cos ϕ) \end{matrix}$ ,

(13)

$\begin{matrix} p = \frac{σ_{k k}}{3},_{}^{} {_{}^{}}_{}^{}_{}^{} q = \sqrt{\frac{3}{2} J_{2}}, \\ J_{2} = s_{i j} s_{i j},_{}^{} {_{}^{}}_{}^{}_{}^{} s_{i j} = σ_{i j} - \frac{σ_{k k}}{3} δ_{i j} \end{matrix}$ 。

(14)

式中 $c$ 为材料黏聚力, $ϕ$ 为材料内摩擦角,可通过试验确定； $e_{f}$ 为子午面上拟合屈服线的双曲偏心率,可通过调整 $e_{f}$ 来调整抗拉强度大小； $η_{f}$ 为微结构张量参数, $h (ξ)$ 为塑性硬化函数, $ξ$ 为塑性内变量； $p$ 为平均应力, $q$ 为偏应力, $J_{2}$ 为应力偏量第二不变量, $θ$ 为Lode角, $K (θ)$ 为关于Lode角的函数。

Lode角函数采用分段形式对屈服面角点进行平滑修正：

$K (θ) = {\begin{matrix} \cos θ - \frac{\sin θ \cos ϕ}{\sqrt{3}} & | θ | \leq θ_{t} \\ K_{1} - K_{2} \sin 3 θ & | θ | > θ_{t} \end{matrix}$ ,

(15)

$\begin{matrix} K_{1} = \frac{1}{3} \cos θ_{t} (3 + \tan θ_{t} \tan 3 θ_{t} + \\ \frac{1}{\sqrt{3}} (\tan 3 θ_{t} - 3 \tan θ_{t}) sign(θ) \sin ϕ) \end{matrix}$ ,

(16)

$K_{2} = \frac{1}{3 \cos 3 θ_{t}} (sign(θ) \sin θ_{t} + \frac{1}{\sqrt{3}} \sin θ_{t} \sin ϕ)$ ,

(17)

式中, $θ_{t}$ 为容许Lode角,当 $| θ | \leq θ_{t}$ 时,与经典Mohr–Coulomb准则一致,当 $| θ | > θ_{t}$ 时,对屈服函数迹线进行平滑处理。

取等效塑性剪应变 ${\bar{γ}}^{p}$ 作为塑性内变量,定义塑性硬化函数为^[12]

$h (ξ) = h_{0} - (h_{1} - h_{0}) \frac{ξ}{B + ξ}$ ,

(18)

$ξ = {\bar{γ}}^{p} = \int d {\bar{γ}}^{p},_{}^{} {_{}^{}}_{}^{} d {\bar{γ}}^{p} = \sqrt{\frac{2}{3} d e_{i j}^{p} d e_{i j}^{p}}$ ,

(19)

式中, $h_{0}$ 为初始屈服阈值, $h_{1}$ 为峰值屈服阈值, $B$ 为描述塑性硬化率大小的参数。

塑性势函数与屈服函数形式一致,即

$g = \sqrt{q^{2} K^{2} (θ) + {(e_{f} c \cos ψ)}^{2}} + η_{f} p \sin ψ$ ,

(20)

式中, $ψ$ 为膨胀角,当 $ψ \neq ϕ$ 时即为非关联流动。

1.3 屈服函数及塑性势函数求导

由于本构数值积分中将涉及到屈服函数及塑性势函数的偏导数,故先行给出偏导数的求解。首先计算屈服函数对应力张量的偏导数,根据式（13）有

$\frac{\partial f}{\partial σ_{i j}} = (\frac{\partial f}{\partial I_{1}} \frac{\partial I_{1}}{\partial σ_{i j}} + \frac{\partial f}{\partial J_{2}} \frac{\partial J_{2}}{\partial σ_{i j}} + \frac{\partial f}{\partial θ} \frac{\partial θ}{\partial σ_{i j}}) + (\frac{\partial f}{\partial η_{f}} \frac{\partial η_{f}}{\partial σ_{i j}})$ ,

(21)

式中,

$\begin{matrix} \frac{\partial I_{1}}{\partial σ_{i j}} = δ_{i j},_{}^{} {_{}^{}}_{}^{} \frac{\partial J_{2}}{\partial σ_{i j}} = s_{i j}, \\ \frac{\partial θ}{\partial σ_{i j}} = \frac{\sqrt{3}}{2 \sqrt{J_{2}^{3}} \cos 3 θ} (\frac{3 J_{3}}{2 J_{2}} s_{i j} - s_{i k} s_{k j} + \frac{2}{3} J_{2} s_{i j}) \end{matrix}$ ,

(22)

$\frac{\partial η_{f}}{\partial σ_{i j}} = \frac{2 η_{0} (A_{k i} σ_{k j} σ_{p q} σ_{p q} - A_{p k} σ_{p q} σ_{k q} σ_{i j})}{{(σ_{m n} σ_{m n})}^{2}}$ ,

(23)

$\frac{\partial f}{\partial I_{1}} = \frac{1}{3} η_{f} h \sin ϕ, \frac{\partial f}{\partial J_{2}} = \frac{1}{2} α_{f} K, \frac{\partial f}{\partial θ} = α_{f} J_{2} \frac{d K}{d θ}$ ,

(24)

$\frac{\partial f}{\partial η_{f}} = h (ξ) (p \sin ϕ - c \cos ϕ)$ ,

(25)

$α_{f} = 3 K / \sqrt{3 J_{2} K^{2} + {(e_{f} c \cos ϕ)}^{2}}$ ,

(26)

$\frac{d K}{d θ} = {\begin{matrix} - \sin θ - \frac{\cos θ \sin ϕ}{\sqrt{3}} & | θ | \leq θ_{t} \\ - 3 K_{2} \cos 3 θ & | θ | > θ_{t} \end{matrix}$ 。

(27)

然后计算塑性势函数对应力张量的偏导数。由于式（20）采取了与屈服函数相一致的形式,因此塑性势函数对应力张量的一阶偏导数只需将内摩擦角 $ϕ$ 替换为膨胀角 $ψ$ ,即

$\frac{\partial g}{\partial σ_{i j}} = (\frac{\partial g}{\partial I_{1}} \frac{\partial I_{1}}{\partial σ_{i j}} + \frac{\partial g}{\partial J_{2}} \frac{\partial J_{2}}{\partial σ_{i j}} + \frac{\partial g}{\partial θ} \frac{\partial θ}{\partial σ_{i j}}) + (\frac{\partial g}{\partial η_{f}} \frac{\partial η_{f}}{\partial σ_{i j}})$ 。

(28)

式中,

$\frac{\partial g}{\partial I_{1}} = \frac{1}{3} η_{f} h \sin ψ, \frac{\partial g}{\partial J_{2}} = \frac{1}{2} α_{g} K, \frac{\partial g}{\partial θ} = α_{g} J_{2} \frac{d K}{d θ}$ ,

(29)

$\frac{\partial g}{\partial η_{f}} = h (ξ) p \sin ψ$ ,

(30)

$α_{g} = 3 K / \sqrt{3 J_{2} K^{2} + {(e_{f} c \cos ψ)}^{2}}$ 。

(31)

1.4 弹塑性刚度矩阵

根据塑性一致性条件,对于式（13）中的屈服函数,其增量形式满足：

$d f = \frac{\partial f}{\partial σ_{i j}} d σ_{i j} + \frac{\partial f}{\partial h} \frac{\partial h}{\partial ξ} d ξ = 0$ 。

(32)

此外,考虑到应力增量和应变增量满足

$d σ_{i j} = D_{i j k l}^{e} d ε_{k l}^{e} = D_{i j k l}^{e} (d ε_{k l} - d ε_{k l}^{p})$ 。

(33)

根据塑性力学,塑性应变增量可由塑性因子 $λ$ 以流动法则表示为

$d ε_{k l}^{p} = d λ \frac{\partial g}{\partial σ_{i j}}$ 。

(34)

等效塑性剪应变 ${\bar{γ}}^{p}$ 也可以流动法则表示为

$d ξ = d {\bar{γ}}^{p} = d λ \sqrt{\frac{2}{3} dev (\frac{\partial g}{\partial σ_{i j}}) dev (\frac{\partial g}{\partial σ_{i j}})}$ 。

(35)

将式（32）～（35）联立,可解得

$d λ = H^{- 1} \frac{\partial f}{\partial σ_{i j}} D_{i j k l}^{e} d ε_{k l}^{p}$ ,

(36)

$\begin{matrix} H = \frac{\partial f}{\partial σ_{i j}} D_{i j k l}^{e} \frac{\partial g}{\partial σ_{k l}} - \\ \frac{\partial f}{\partial h} \frac{\partial h}{\partial ξ} \sqrt{\frac{2}{3} dev (\frac{\partial g}{\partial σ_{i j}}) dev (\frac{\partial g}{\partial σ_{i j}})} \end{matrix}$ ,

(37)

式中,H为塑性硬化模量。

将式（36）回代入式（33）,可得应变增量与应力增量的一致性表达：

$d σ_{i j} = D_{i j k l}^{e} (d ε_{k l} - d λ \frac{\partial g}{\partial σ_{i j}}) = D_{i j k l}^{ep} d ε_{k l}$ ,

(38)

$D_{i j k l}^{ep} = D_{i j k l}^{e} - H^{- 1} D_{i j p q}^{e} \frac{\partial f}{\partial σ_{p q}} \frac{\partial g}{\partial σ_{m n}} D_{m n k l}^{e}$ ,

(39)

式中, $D_{i j k l}^{ep}$ 即为一致性弹塑性刚度矩阵。

2. 数值算例

基于建立的柱状节理岩体的各向异性弹塑性本构模型,结合FLAC^3D软件对该模型进行了二次开发,并对柱状节理相似结构材料三轴岩石力学试验成果进行了模拟。根据国际岩石力学学会（ISRM）推荐标准,圆柱样直径为50 mm,柱体高度为100 mm,网格数量10000个,底面固定,顶面控制加载速度为5×10^-7 m/s,柱体四周施加固定围压4 MPa。

可由不同空间方向的柱状节理岩体的峰值强度进行拟合确定微结构张量系数。图2给出了强度参数关于方向角的微结构张量参数拟合曲线。

图 2 不同拟合阶数下单轴抗压强度随柱体倾角变化规律

Figure 2. Variation of uniaxial compressive strength under different fitting orders with column inclination

下载: 全尺寸图片幻灯片

为了解释固有的各向异性,假设单轴抗压强度 $f_{c}$ 受样本方向的影响,并结合与式（12）相似的分布函数来描述其变化。为了能反映材料各向异性屈服强度,并防止欠拟合和过拟合,本文采用二阶微结构张量来描述柱状节理岩体强度参数的空间分布屈服函数：

$f_{c} = \hat{f} {1 + A (1 - 3 \cos^{2} β) + b {[A (1 - 3 \cos^{2} β)]}^{2}}$ ,

(40)

式中,

$\hat{f} = 12.76 MPa, A = 1.894, b = 0.3526$ 。

(41)

其余计算参数可根据室内试验曲线^[1]数据,应用单纯形最优化方法在MATLAB中构造寻优函数编程计算得到。模型计算参数如表1所示。

表 1 柱状节理岩体基本力学参数表

Table 1. Mechanical parameters of columnar joints

编号	^#1	^#2	^#3	^#4	^#5	^#6	^#7	^#8	^#9	^#10	^#11	^#12	^#13	^#14	^#15	^#16
参数	E₁/MPa	E₂/MPa	G₁₃/MPa	$v_{12}$	$v_{13}$	c/MPa	$φ$ /(°)	$e_{f}$	$\overset{⌢}{η}$	A	b	$h_{0}$	$h_{1}$	$e_{f}$	B	$θ_{t}$
数值	9.81	22.1	4.69	0.22	0.06	0.72	35.2	0.75	0.23	1.89	0.35	0.77	1.21	0.53	0.0018	25°

下载: 导出CSV

| 显示表格

应用表1给出的力学参数,在FLAC^3D中可以调用自定义本构模型来模拟柱状节理岩体的弹塑性变形行为。数值模拟结果如图3和图4所示。

图 3 三轴数值模拟位移云图（倾角75°,围压4 MPa）

Figure 3. Displacement nephogram by triaxial numerical simulation (angle of 75° and confining pressure of 4 MPa)

下载: 全尺寸图片幻灯片

图 4 三轴试验与模拟结果对比（倾角75°,围压4 MPa）^[13]

Figure 4. Comparison between triaxial test and simulated results (angle of 75°and confining pressure of 4 MPa)

下载: 全尺寸图片幻灯片

本文将柱状节理岩体概化为横观各向同性连续体,建立描述岩体力学特性的宏观本构模型。该模型采用了一种考虑抗拉强度和屈服面角点平滑的修正Mohr–Coulomb准则,将其扩展到了空间各向异性,包含一个各向异性标量参数,该参数用应力和结构取向张量的混合不变量表示。

结果表明,柱体倾角为75°时节理岩体的变形具有明显的方向性,模拟所得应力–应变曲线与试验数据^[13]基本吻合。

应该强调的是,在对问题进行分析时,式（18）中硬化特性仅通过调用一个标量参数B来描述,由于式（18）仅限于硬化状态,不包括应变软化分支。当然也可以采用一种更复杂的方法,并假设这些特征受到围压的影响,然而这样将不可避免地生成一组额外的参数及公式。此外,该方法目的在于研究工程尺度柱状节理岩体的物理力学特性,应用了一种折中的处理策略,当研究尺度较小时,则需要进一步考虑对岩体内部的结构面开展精细模拟。

3. 白鹤滩高坝坝基开挖工程数值分析

3.1 白鹤滩坝区工程地质特征及网格建模

白鹤滩水电站坐落于金沙江下游云南巧家县镜内,是金沙江下游干流河段梯级开发的第二个梯级电站。主坝为混凝土双曲拱坝,最大坝高289 m,计划2022年建设完工。白鹤滩水电站坝区属中山峡谷地貌,坝址为单斜地层,主要为二叠系上统峨眉山组玄武岩,如图5,可分为11个岩流层,岩层走向N30～50°E,倾向SE,倾角15°～25°。岩层岩性主要分为斜斑玄武岩、隐晶玄武岩、柱状节理玄武岩、杏仁状玄武岩、角砾熔岩和凝灰岩等。右岸坡顶顶部出露三叠系下统飞仙关组泥质粉砂岩、粉砂质泥岩及少量砂岩^[2]。其中出露 $P_{2} β_{3}$ 层为柱状节理玄武岩,柱体长度为2～3 m,柱体倾角70°～85°,柱体截面主要为不规则五边形和四边形,截面边长为10～20 cm。柱体内部隐节理发育,开挖前节理闭合紧密^[14]。

图 5 拱坝轴线地质剖面示意图

Figure 5. Illustration of engineering geological section of arch axis

下载: 全尺寸图片幻灯片

图6（a）为开挖模拟所用网格模型,坐标原点位于拱坝轴线,建模范围为1300 m×1500 m×1296 m。在模型局部区域,对单元进行了细化剖分,网格平均边长约为1.5 m。采用Hypermesh剖分三维实体,并形成网格文件,然后将数据转换为相应的软件格式,建立了数值计算网格模型,单元总数145万个,节点26万个。坝基开挖边坡及柱状节理玄武岩体为计算分析重点关注的部位。除了层间错动带,该网格模型还考虑了断层F16、F17,以及柱状节理玄武岩层 $P_{2} β_{3}$ 中的层内错动带LS3319、LS3319-1、RS3315和RS336,如图6（b）所示。

图 6 白鹤滩坝基开挖计算网格模型

Figure 6. Calculation grid model for excavation of foundation of …Baihetan Dam

下载: 全尺寸图片幻灯片

3.2 模型参数及工况设置

对模型前后左右4个侧面施加法向约束,对底面施加x,y,z三向约束作为边界条件。模型自下而上共划分有8个地层,模型计算参数根据工程地质条件和岩石力学参数试验等确定,主要岩体力学参数如表2所示。

表 2 主要岩体基本物理力学参数

Table 2. Basic mechanical parameters of main rock mass

岩层分类	$ρ$ /(kg·m^-3)	E/GPa	$ν$	c/MPa	$φ$ /(°)
C₄层面以上	2700	11.0	0.23	1.10	47.7
C₄~C₃层面	2700	13.2	0.23	1.37	51.8
C_3-1~柱状节理顶面	2600	10.2	0.24	1.07	48.0
柱状节理层	2700	10.0	0.25	0.90	47.7
角砾熔岩层	2600	9.0	0.24	1.08	47.2
微新岩岩层	2850	22.0	0.22	2.10	57.2
C₄	2130	0.4	0.35	0.10	28.0
C₃	2130	0.4	0.34	0.10	28.0
C_3-1	2050	0.7	0.35	0.04	28.0
F₁₆	2130	1.5	0.35	0.20	29.2
F₁₇风化上段	2000	0.7	0.35	0.04	15.6
F₁₇风化下段	2130	1.2	0.35	0.05	20.8
F₁₇微新岩体	2130	2.0	0.35	0.20	29.2
LS₃₃₁	2150	1.5	0.35	0.30	30.0
LS₃₃₁₉	2110	0.3	0.35	0.10	26.0
LS_3319-1	2000	0.3	0.35	0.15	27.0
RS₃₃₆	2000	0.3	0.35	0.15	27.0
RS₃₃₁₅	2000	0.3	0.35	0.15	26.0

下载: 导出CSV

| 显示表格

考虑左右岸坝基岩石工程现场实际施工工序,依照表3对边坡开挖过程设置了10个开挖步。

表 3 模型开挖工况设置

Table 3. Setting of model excavation conditions

开挖步编号	开挖台阶高程	现场开挖完成时间
^#1	左岸800 m	2013.10.24
^#2	左岸700 m	2014.08.19
^#3	右岸800 m	2014.11.10
^#4	左岸650 m	2014.12.03
^#5	右岸700 m	2015.03.30
^#6	右岸650 m	2015.06.24
^#7	右岸600 m	2016.01.27
^#8	左岸600 m	2016.06.03
^#9	600 m高程以下	2016.07.26
^#10	保护层	2016.08.14

下载: 导出CSV

| 显示表格

3.3 计算结果分析

图7为白鹤滩开挖模型左岸坝基不同高程处测点位移变化曲线图。左岸坝基建基面岩体位移随开挖进程逐渐增大；开挖完成时,岩体位移随岩体所在高程减小而先增大后减小,坝基底部柱状节理岩体位移大于上部非柱状节理岩体。

图 7 左岸坝基不同高程处测点位移变化曲线图

Figure 7. Curves of displacement change at measurement points at different elevations of dam foundation at left bank

下载: 全尺寸图片幻灯片

图8为考虑柱状节理岩体各向异性时白鹤滩坝基边坡开挖过程位移变化云图。

图 8 考虑柱状节理岩体各向异性模型开挖过程位移变化云图

Figure 8. Cloud charts of displacement change in excavation process considering anisotropic model of columnar jointed rock mass

下载: 全尺寸图片幻灯片

开挖完成时,建基面最大位移超过50 mm,平均位移为25～40 mm；坝基开挖全过程中,最大位移出现在各开挖步开挖梯段的边坡坡脚及开挖台阶面上；岩体变形范围和变形程度随开挖进程逐渐增大；表层岩体的变形较大,变形程度随岩体深度增加而逐渐减小；岩体变形围绕断层和错动带发展,断层和错动带周围的岩体变形显著大于不含断层和错动带的岩体（实际现场开挖中对坝基出露断层采取了置换处理）。

图9为左岸坝基658 m高程处多点变位计Mzjc-2监测所得位移时程曲线,以及分别采用各向异性模型和各向同性模型计算所得的变形曲线。从图中可以看到,相比于只采用各向同性模型,采用各向异性模型的位移计算结果与现场监测结果更为接近。

图 9 左岸坝基多点变位计Mzjc-2监测结果与计算结果对比

Figure 9. Comparison between monitoring and calculated results of Mzjc-2 multi-point displacement meter

下载: 全尺寸图片幻灯片

表4为坝基开挖完成时刻左岸坝基柱状节理试验区及左岸坝肩下游边坡处多点变位计变形监测结果与数值模型计算结果的对比。

表 4 坝基开挖结束时监测变形与计算位移对比

Table 4. Comparison of monitored deformations and calculated displacements at completion end of excavation of dam foundation

测点编号	监测结果/mm	各向异性模型		各向同性模型
测点编号	监测结果/mm	计算值/mm	误差/%	计算值/mm	误差/%
Mzjc-2	25.4	23.7	6.7	19.2	24.4
Mzjc-3	24.9	21.8	12.4	18.0	27.7
M_ZBJ-3	6.6	5.6	15.2	4.9	25.8

下载: 导出CSV

| 显示表格

结果表明：各向异性模型位移计算误差明显小于各向同性模型,这在柱状节理区域表现更为显著。

为生成质量较高的三维网格模型,本文所用开挖模型中开挖梯段高差（50～100 m）与实际工程施工台阶（5～10 m）存在较大差距,这在模拟计算结果中主要体现在变形的时间节点与实际监测结果存在一定出入。另一方面,本文计算所采用的各向异性弹塑性本构模型不能完整体现节理岩体卸荷松弛的时效性。在后续的工作中,考虑引入各向异性时效损伤以本文的相关工作进一步深化完善。

4. 结论

（1）结合微结构张量理论扩展了经典莫尔库仑模型,完整推导了柱状节理岩体的各向异性弹塑性本构模型。考虑各向异性的本构模型可以直观反映岩体强度的各向异性空间分布,计算所需参数均可由试验结果计算得到。

（2）结合FLAC^3D二次开发接口对自定义本构模型进行了开发。对所开发模型进行了数值验证,数值计算结果与柱状节理相似材料室内三轴压缩试验结果基本一致,表明该模型可以体现柱状节理材料的空间各向异性特性。

（3）应用所开发本构模型对白鹤滩水电站坝基边坡开挖过程进行了模拟。计算结果表明,与未考虑各向异性相比,考虑柱状节理岩体强度各向异性时坝基边坡位移更接近于现场多点位移计监测结果,本文研究工作对白鹤滩水电工程建设有一定参考价值。

图 1 初始粒径分布曲线

Figure 1. Initial particle-size distribution curves

下载: 全尺寸图片幻灯片

图 2 粒径分布曲线

Figure 2. Particle-size distribution curves

下载: 全尺寸图片幻灯片

图 3 粒径分布曲线的演化

Figure 3. Evolution of particle-size distribution curves

下载: 全尺寸图片幻灯片

图 4 不同粒组含量的演化

Figure 4. Evolution of different particle fraction contents

下载: 全尺寸图片幻灯片

图 5 土样FCS的粒组相对含量的演化

Figure 5. Evolution of relative fraction content of particles of FCS samples

下载: 全尺寸图片幻灯片

图 6 偏应力与轴向应变的关系

Figure 6. Relationship between deviator stress and axial strain

下载: 全尺寸图片幻灯片

图 7 体应变与轴向应变的关系

Figure 7. Relationship between volumetric strain and axial strain

下载: 全尺寸图片幻灯片

图 8 主应力比与轴向应变的关系

Figure 8. Relationship between principal stress ratio and axial strain

下载: 全尺寸图片幻灯片

图 9 破碎率对分形维数的影响

Figure 9. Influences of breakage index on fractal dimension

下载: 全尺寸图片幻灯片

图 10 轴向应变对破碎率的影响

Figure 10. Influences of axial strain on breakage index

下载: 全尺寸图片幻灯片

图 11 体应变对破碎率的影响

Figure 11. Influences of volumetric strain on breakage index

下载: 全尺寸图片幻灯片

图 12 偏应力对破碎率的影响

Figure 12. Influences of deviator stress on breakage index

下载: 全尺寸图片幻灯片

图 13 主应力比对破碎率的影响

Figure 13. Influences of principal stress ratio on breakage index

下载: 全尺寸图片幻灯片

图 14 峰值主应力比与破碎率的关系

Figure 14. Relationship between peak principal stress ratio and breakage index

下载: 全尺寸图片幻灯片

图 15 峰值内摩擦角与破碎率的关系

Figure 15. Relationship between peak angle of internal friction and breakage index

下载: 全尺寸图片幻灯片

表 1 试样的物理参数

Table 1 Physical parameters of samples

试样编号	初始干密度ρ_d0/(g·cm^-3)	初始孔隙比e₀	最大孔隙比e_max	最小孔隙比e_min	相对密实度D_r
SCS	1.24	1.21	1.41	0.94	0.43
DCS	1.35	1.03	1.28	0.79	0.51
FCS	1.45	0.90	1.20	0.62	0.52

下载: 导出CSV

表 2 分形维数和拟合相关系数

Table 2 Values of fractal dimension and correlation coefficient

围压 $σ_{3}$ /MPa	轴向应变ε_a/%	分形维数D	相关系数R²
0	0	2.150	1.000
0.4	5	2.158	0.999
	10	2.168	0.999
	15	2.177	0.998
	20	2.185	0.997
0.8	5	2.163	0.996
	10	2.171	0.997
	15	2.206	0.998
	20	2.223	0.998
1.6	5	2.181	0.994
	10	2.217	0.992
	15	2.247	0.996
	20	2.275	0.998
3.2	5	2.196	0.998
	10	2.257	0.993
	15	2.307	0.993
	20	2.347	0.992

下载: 导出CSV

参考文献(28)

[1]	沈建华, 汪稔. 钙质砂的工程性质研究进展与展望[J]. 工程地质学报, 2010, 18(增刊1): 26-32. https://cpfd.cnki.com.cn/Article/CPFDTOTAL-GCDZ201011001008.htm SHEN Jian-hua, WANG Ren. Study on engineering properties of calcareous sand[J]. Journal of Engineering Geology, 2010, 18(S1): 26-32. (in Chinese) https://cpfd.cnki.com.cn/Article/CPFDTOTAL-GCDZ201011001008.htm
[2]	尹振宇, 许强, 胡伟. 考虑颗粒破碎效应的粒状材料本构研究：进展及发展[J]. 岩土工程学报, 2012, 34(12): 2170-2180. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201212006.htm YIN Zhen-yu, XU Qiang, HU Wei. Constitutive relations for granular materials considering particle crushing: review and development[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(12): 2170-2180. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201212006.htm
[3]	张家铭, 汪稔, 石祥锋, 等. 侧限条件下钙质砂压缩和破碎特性试验研究[J]. 岩石力学与工程学报, 2005, 24(8): 3327-3331. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200518021.htm ZHANG Jia-ming, WANG Ren, SHI Xiang-feng, et al. Compression and crushing behavior of calcareous sand under confined compression[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(8): 3327-3331. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200518021.htm
[4]	秦月, 姚婷, 汪稔, 等. 基于颗粒破碎的钙质沉积物高压固结变形分析[J]. 岩土力学, 2014, 35(11): 3123-3128. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201411014.htm QIN Yue, YAO Ting, WANG Ren, et al. Particle breakage-based analysis of deformation law of calcareous sediments under high-pressure consolidation[J]. Rock and Soil Mechanics, 2014, 35(11): 3123-3128. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201411014.htm
[5]	张季如, 张弼文, 胡泳, 等. 粒状岩土材料颗粒破碎演化规律的模型预测研究[J]. 岩石力学与工程学报, 2016, 35(9): 1898-1905. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201609019.htm ZHANG Ji-ru, ZHANG Bi-wen, HU Yong, et al. Predicting the particle breakage of granular geomaterials[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(9): 1898-1905. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201609019.htm
[6]	吕亚茹, 李治中, 李浪. 高应力状态下钙质砂的一维压缩特性及试验影响因素分析[J]. 岩石力学与工程学报, 2019, 38(增刊1): 3142-3150. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2019S1055.htm LÜ Ya-ru, LI Zhi-zhong, LI Lang. One-dimensional compression behavior of calcareous sand and its experimental technology under high stress conditions[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(S1): 3142-3150. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2019S1055.htm
[7]	COOP M R, SORENSEN K K, BODAS FREITAS T, et al. Particle breakage during shearing of a carbonate sand[J]. Géotechnique, 2004, 54(3): 157-163. doi: 10.1680/geot.2004.54.3.157
[8]	MIAO G, AIREY D. Breakage and ultimate states for a carbonate sand[J]. Géotechnique, 2013, 63(14): 1221-1229. doi: 10.1680/geot.12.P.111
[9]	何建乔, 魏厚振, 孟庆山, 等. 大位移剪切下钙质砂破碎演化特性[J]. 岩土力学, 2018, 39(1): 165-172. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801021.htm HE Jian-qiao, WEI Hou-zhen, MENG Qing-shan, et al. Evolution of particle breakage of calcareous sand under large displacement shearing[J]. Rock and Soil Mechanics, 2018, 39(1): 165-172. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801021.htm
[10]	黄宏翔, 陈育民, 王建平, 等. 钙质砂抗剪强度特性的环剪试验[J]. 岩土力学, 2018, 39(6): 2082-2088. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201806021.htm HUANG Hong-xiang, CHEN Yu-min, WANG Jian-ping, et al. Ring shear tests on shear strength of calcareous sand[J]. Rock and Soil Mechanics, 2018, 39(6): 2082-2088. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201806021.htm
[11]	纪文栋, 张宇亭, 王洋, 等. 循环单剪下珊瑚钙质砂和普通硅质砂剪切特性对比研究[J]. 岩土力学, 2018, 39(增刊1): 282-288. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2018S1035.htm JI Wen-dong, ZHANG Yu-ting, WANG Yang, et al. Comparative study of shear performance between coral sand and siliceous sand in cycles simple shear test[J]. Rock and Soil Mechanics, 2018, 39(S1): 282-288. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2018S1035.htm
[12]	蔡正银, 侯贺营, 张晋勋, 等. 考虑颗粒破碎影响的珊瑚砂临界状态与本构模型研究[J]. 岩土工程学报, 2019, 41(6): 989-995. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201906002.htm CAI Zheng-yin, HOU He-ying, ZHANG Jin-xun, et al. Critical state and constitutive model for coral sand considering particle breakage[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(6): 989-995. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201906002.htm
[13]	蔡正银, 侯贺营, 张晋勋, 等. 密度与应力水平对珊瑚砂颗粒破碎影响试验研究[J]. 水利学报, 2019, 50(2): 184-191. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201902004.htm CAI Zheng-yin, HOU He-ying, ZHANG Jin-xun, et al. Experimental study on the influence of density and stress level on particle breakage of coral sand[J]. Journal of Hydraulic Engineering, 2019, 50(2): 184-191. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201902004.htm
[14]	王刚, 叶沁果, 查京京. 珊瑚礁砂砾料力学行为与颗粒破碎的试验研究[J]. 岩土工程学报, 2018, 40(5): 802-810. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201805006.htm WANG Gang, YE Qin-guo, ZHA Jing-jing. Experimental study on mechanical behavior and particle crushing of coral sand-gravel fill[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(5): 802-810. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201805006.htm
[15]	陈火东, 魏厚振, 孟庆山, 等. 颗粒破碎对钙质砂的应力–应变及强度影响研究[J]. 工程地质学报, 2018, 26(6): 1490-1498. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201806011.htm CHEN Huo-dong, WEI Hou-zhen, MENG Qing-shan, et al. The study on stress-strain-strength behavior of calcareous sand with particle breakage[J]. Journal of Engineering Geology, 2018, 26(6): 1490-1498. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201806011.htm
[16]	刘汉龙, 肖鹏, 肖杨, 等. MICP胶结钙质砂动力特性试验研究[J]. 岩土工程学报, 2018, 40(1): 38-45. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201801003.htm LIU Han-long, XIAO Peng, XIAO Yang, et al. Dynamic behaviors of MICP-treated calcareous sand in cyclic tests[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(1): 38-45. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201801003.htm
[17]	王刚, 查京京, 魏星. 循环三轴应力路径下钙质砂颗粒破碎演化规律[J]. 岩土工程学报, 2018, 41(4): 755-760. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201904025.htm WANG Gang, ZHA Jing-jing, WEI Xing. Evolution of particle crushing of carbonate sand under cyclic triaxial stress path[J]. Chinese Journal of Geotechnical Engineering, 2018, 41(4): 755-760. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201904025.htm
[18]	张家铭, 张凌, 蒋国盛, 等. 剪切作用下钙质砂颗粒破碎试验研究[J]. 岩土力学, 2008, 29(10): 2789-2793. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200810039.htm ZHANG Jia-ming, ZHANG Ling, JIANG Guo-sheng, et al. Research on particle crushing of calcareous sands under triaxial shear[J]. Rock and Soil Mechanics, 2008, 29(10): 2789-2793. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200810039.htm
[19]	张晨阳, 谌民, 胡明鉴, 等. 细颗粒组分含量对钙质砂抗剪强度的影响[J]. 岩土力学, 2019, 40(增刊1): 195-202. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2019S1029.htm ZHANG Chen-yang, CHEN Min, HU Ming-jian, et al. Effect of fine particles content on shear strength of calcareous sand[J]. Rock and Soil Mechanics, 2019, 40(S1): 195-202. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2019S1029.htm
[20]	沈扬, 沈雪, 俞演名, 等. 粒组含量对钙质砂压缩变形特性影响的宏细观研究[J]. 岩土力学, 2019, 40(10): 3733-3740. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910006.htm SHEN Yang, SHEN Xue, YU Yan-ming, et al. Macro-micro study of compressive deformation properties of calcareous sand with different particle fraction contents[J]. Rock and Soil Mechanics, 2019, 40(10): 3733-3740. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910006.htm
[21]	TURCOTTE D L. Fractals in geology and geophysics[J]. Pure Application Geophysics, 1989, 131(1/2): 171-196.
[22]	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.
[23]	张季如, 胡泳, 张弼文, 等. 石英砂砾破碎过程中粒径分布的分形行为研究[J]. 岩土工程学报, 2015, 37(5): 784-791. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201505004.htm ZHANG Ji-ru, HU Yong, ZHANG Bi-wen, et al. Fractal behavior of particle-size distribution during particle crushing of quartz sand and gravel[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(5): 784-791. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201505004.htm
[24]	HARDIN B O. Crushing of soil particles[J]. Journal of Geotechnical Engineering, 1985, 111(10): 1177-1192.
[25]	TURCOTTE D L. Fractals and fragmentation[J]. Journal of Geophysical Research, 1986, 91(B2): 1921-1926.
[26]	纪文栋, 张宇亭, 裴文斌, 等. 加载方式和应力水平对珊瑚砂颗粒破碎影响的试验研究[J]. 岩石力学与工程学报, 2018, 37(8): 1953-1961. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201808018.htm JI Wen-dong, ZHANG Yu-ting, PEI Wen-bin, et al. Influence of loading method and stress level on the particle crushing of coral calcareous sand[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(8): 1953-1961. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201808018.htm
[27]	张季如, 祝杰, 黄文竞. 侧限压缩下石英砂砾的颗粒破碎特性及其分形描述[J]. 岩土工程学报, 2008, 30(6): 783-789. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200806001.htm ZHANG Ji-ru, ZHU Jie, HUANG Wen-jing. Crushing and fractal behaviors of quartz sand-gravel particles under confined compression[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(6): 783-789. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200806001.htm
[28]	VESIC A S, CLOUGH G W. Behavior of granular materials under high stresses[J]. Journal of the Soil Mechanics and Foundations Division, ASCE, 1968, 94(SM3): 661-688.

施引文献(35)

期刊类型引用(24)

1.	夏晶晶，贺姣姣. 基于直剪试验土石混合体路基的稳定性问题研究. 山西建筑. 2025(01): 125-128+172 . 百度学术
2.	郜力君，郝捷. 含石率及坡度对土石混合路堑边坡稳定性的影响分析. 交通世界. 2025(Z2): 154-156 . 百度学术
3.	邹弈，朱碧堂，吴颖彪，周宇航. 断层破碎带黏性夹泥岩体与基岩接触界面大型直剪试验研究. 华东交通大学学报. 2025(01): 45-51 . 百度学术
4.	李利萍，余泓浩，李秋雨，潘一山. 砂岩不同含水特性对超低摩擦效应影响试验研究. 力学学报. 2025(03): 687-700 . 百度学术
5.	程虎，李重情，穆朝民. 冻结温度对不同粒径冻土石混合体劈裂特性的影响. 煤矿安全. 2024(01): 160-166 . 百度学术
6.	代兴先，宋杨，张志彬，齐子怡，刘棋瑞. 土石混合料力学特性和颗粒破碎研究. 河北水利电力学院学报. 2024(01): 1-8 . 百度学术
7.	李慎刚，石云方，刘晋宁，蒋琛. 碎石土路基填料压实及渗透特性. 工程科学学报. 2024(05): 918-926 . 百度学术
8.	SHI Yunfang，LI Shengang，JIANG Chen，LIU Jinning. Gravel hardness effect on compaction characteristics of gravelly soil. Journal of Mountain Science. 2024(04): 1432-1443 . 必应学术
9.	谢周州，赵炼恒，李亮，黄栋梁，张子健，周靖. 基于振动台试验的不同含石率土-石混合体边坡地震动响应差异性研究. 岩土力学. 2024(08): 2324-2337 . 百度学术
10.	王楠楠，高霞，张吉哲，张保勇，吴强. 基于平行黏结模型的含瓦斯水合物煤体宏细观力学性质研究. 煤炭学报. 2024(S1): 314-326 . 百度学术
11.	魏东旭，彭雄志，张佳，范帅，冯瑨. 石灰岩碎石土大型直剪试验研究. 勘察科学技术. 2024(04): 1-5+43 . 百度学术
12.	李刚，尹小涛. 基于数字化施工的山区公路土石混合弃渣工程尺度参数测定方法研究. 水利与建筑工程学报. 2024(06): 76-82 . 百度学术
13.	陈小翔. 碎石土路基填料压实及渗透特性研究. 江西建材. 2024(10): 52-54 . 百度学术
14.	杨忠平，李进，刘浩宇，张益铭，刘新荣. 土石混合体-基岩界面剪切力学特性块石尺寸效应. 岩土力学. 2023(04): 965-974 . 百度学术
15.	刘建平，周花玉，何天奎，余镜南，张坤，潘玉丛，刘泉声. 含根量对根–土石复合体的抗剪强度影响试验初探——以垂丝海棠为例. 岩石力学与工程学报. 2023(S1): 3618-3628 . 百度学术
16.	刘飞禹，孔剑捷，姚嘉敏. 含石量和压实度对格栅-土石混合体界面剪切特性的影响. 岩土工程学报. 2023(05): 903-911 . 本站查看
17.	韩志洋，曹志翔，黄开放. 基于离散元模拟的土石混合体剪切与变形特性研究. 中国农村水利水电. 2023(05): 238-244 . 百度学术
18.	唐丽云，黄涛，汪卫兵，金龙，孙强，李国玉，罗滔. 冻融循环下土石混合体-混凝土界面剪切特性及孔隙结构演化特征试验研究. 中南大学学报(自然科学版). 2023(05): 1954-1969 . 百度学术
19.	张俊云，张乐，高福洲，唐永吉，何卓岭，王鹰. 干湿循环下红层土石混合料强度及变形特性的试验研究. 西南交通大学学报. 2023(06): 1394-1404 . 百度学术
20.	程晓颖，乔婷，秦建敏，季顺迎. 复杂形态碎石颗粒的三维离散元模拟及试验验证. 计算机辅助工程. 2023(04): 40-47 . 百度学术
21.	YANG Zhong-ping，LI Shi-qi，TIAN Xin，HU Yuan-xin，LI Wan-kun. Cumulative damage effect on debris slopes under frequent microseisms. Journal of Mountain Science. 2022(03): 781-797 . 必应学术
22.	崔熙灿，张凌凯，王建祥. 高堆石坝砂砾石料的细观参数反演及三轴试验模拟. 农业工程学报. 2022(04): 113-122 . 百度学术
23.	宋颖能. 朱家涧水库大坝土石方及渗控工程实施探讨. 内蒙古煤炭经济. 2022(08): 160-162 . 百度学术
24.	覃国强. 不同含石率及坡度条件下土石混合体路堑高边坡失稳机制分析. 福建交通科技. 2022(06): 21-25 . 百度学术