-
摘要:
土石混合体边坡在我国分布范围广, 且其材料组成复杂, 受到了众多学者的关注。为了科学合理地评估土体参数的空间变异性对土石混合体边坡稳定性的影响, 基于随机场理论, 选取有效抗剪强度参数黏聚力
c 和内摩擦角φ 作为随机变量, 采用局部平均法模拟随机场, 在MATLAB中实现随机场参数生成, 在考虑土石混合体块石真实形状和块石含量基础上, 采用Python语言脚本方式, 将随机场参数映射到有限元软件中的土石混合体边坡, 应用强度折减法计算边坡稳定安全系数。计算结果显示, 土石混合体边坡的稳定安全系数符合正态分布, 随着土石混合体边坡块石含量增加, 边坡稳定安全系数的均值从1.005增长至1.095, 边坡也由浅层破坏逐渐发展为深层破坏。块石质量分数均为35%时, 块石粒径较大土石混合体边坡稳定性安全系数为1.334, 块石粒径较小土石混合体边坡稳定安全系数为1.064。相比于确定性计算结果, 考虑土体参数空间变异性的稳定安全系数更高。因此, 在进行土石混合体边坡稳定性分析时, 应充分考虑土体有效抗剪强度参数空间变异性, 避免出现设计过于保守的情况。Abstract:The extensive distribution and complex material composition of soil-rock mixture slopes in China have attracted significant attention from scholars. Objective This study aims to scientifically and rationally assess the impact of the spatial variability of soil parameters on the stability of soil-rock mixture slopes.
Methods Based on the random field theory, the effective shear strength parameters cohesion
c and internal friction angleφ are selected as random variables. The local averaging method is used to simulate the random field, with random field parameter generation conducted in MATLAB. Python scripts are employed to map the random field parameters to the soil-rock mixture slope via finite element software, accounting for the actual shape and content of block rocks in the soil-rock mixture. The strength reduction method is then applied to calculate the slope stability safety factor.Results The results reveal that the stability safety factor of soil-rock mixture slopes follows a normal distribution. As the block stone content increases, the mean value of the stability safety factor rises from 1.005 to 1.095, reflecting a transition from shallow to deep failure. For block stone content of 35%, the stability safety factor reaches 1.334 for larger block stones and 1.064 for smaller ones. Compared to deterministic calculation results, incorporating the spatial variability of the soil parameters yields higher stability safety factor.
Conclusion Therefore, in soil-rock mixture slope stability analyses, the spatial variability of effective shear strength parameters must be fully considered to prevent overly conservative designs.
-
土石混合体(soil-rock mixture, 简称SRM)是指第四纪以来形成的,由具有一定工程尺度、较高强度的块石、细粒土体及孔隙构成且具有一定含石量的极端不均匀松散岩土介质系统[1-2],天然的土石混合体赋存于残坡积物、崩坡积物和冲洪积物这些天然堆积地层中[3]。由土石混合体组成的边坡具有材料参数空间变异性大、块石分布不均匀性强的特点。自然界中存在大量土石混合体边坡[4],在我国西南山区分布尤为广泛,这些边坡的稳定是保证人民生命财产安全的前提[5-6]。
土石混合体边坡因其分布范围广、材料组成复杂等特点,受到了众多学者的关注。在土石混合体边坡稳定性分析方面,学者将土石混合体边坡简化为块石镶嵌于土体中的二相模型,通过不同的建模方法进行边坡稳定性计算。LI等[4]通过随机生成算法建立了含有多边形块石的土石混合体边坡,采用二维修正DDA-SPH耦合方法,探究了块石尺寸和块石含量对边坡稳定安全系数的影响;邵帅等[7]采用数字图像技术建立了含有四边形块石的岩土边坡,应用强度折减法计算了块石空间分布对边坡稳定性的影响;LIU等[5]提出了一种基于数字图像处理和A-star算法的土石混合体边坡稳定性分析方法;GAO等[8]研究了超大块石空间分布对边坡稳定性影响,并提出了关键岩块的概念,位于坡脚处的大粒径岩块对边坡稳定性有利;然而这些模型并未考虑块石真实的二维轮廓形状及土体力学参数的空间变异性,与真实土石混合体边坡的情况有较大出入,根据现有结论显示,不考虑土体参数空间变异性会导致稳定性计算结果偏小[9]。
天然土体由于沉积和后沉积过程、荷载历史条件、物理风化过程和化学风化过程以及自然搬运等作用,导致土体力学参数呈现出空间变异性[10-11],因此,为了科学合理地评估土体力学参数的空间变异性对边坡稳定性的影响,VANMARCKE[12]首先提出了利用随机场理论表征土体参数的空间变异性。此后,许多学者将随机场理论应用于工程可靠度分析中。薛亚东等[13]应用局部平均理论生成二维各向异性随机场,采用强度折减法计算了边坡稳定安全系数,研究结果发现了随机场在不同方向的变异性对边坡稳定影响不一致;LIU等[14]将土石混合体边坡视为二相随机介质模型,利用高斯场的非线性平移变换来模拟土体和块石,发现土石混合体边坡的安全系数对土体体积分数敏感,此方法虽然不需要建立复杂的边坡几何模型,但模型中块石的分布情况与实际相差较远;KASAMA等[15]将土体参数的空间变异性通过随机场离散技术引入数值极限分析中,提出了一种评价岩土工程稳定性问题的概率方法;何成等[16]考虑土体参数空间变异性,分析了多种工况下三峡库区赵树岭滑坡的潜在破坏模式及稳定性。以上方法充分考虑了土体参数的空间变异性,并将所提方法通过均质土坡进行验证,但并未考虑块石存在的土石混合体边坡,本研究将随机场理论应用于土石混合体边坡稳定性分析中,填补这一研究领域的空白。
本研究基于随机场理论,选取有效抗剪强度参数黏聚力(c)和内摩擦角(φ)作为随机变量,采用局部平均法模拟随机场,在MATLAB中实现随机场参数生成,在考虑土石混合体块石真实形状和块石含量基础上,采用Python语言脚本方式,将随机场参数映射到有限元软件中的土石混合体边坡,应用强度折减法计算边坡稳定安全系数,并分析边坡稳定安全系数的概率分布情况,探讨块石含量、块石粒径对土石混合体边坡稳定性影响。
1. 局部平均法模拟随机场
本研究采用二维平稳随机场[17]来描述土体力学参数的空间变异性。如图 1所示,X(t1, t2)为二维连续平稳随机场,其均值和方差分别为μ和σ2。A=T1T2为该随机场中的一个四边形单元,(t1, t2)为其中心点坐标,则随机场在单元A内的局部平均为:
XA(t1,t2)=1A∫t1+T1/2t1−T1/2∫t2+T2/2t2−T2/2X(t1,t2)dt1 dt2, (1) 式中:X(t1, t2)对应的局部平均随机场XA(t1, t2)可用各种随机场单元的均值E(XA)、方差Var(XA)和协方差Cov(XA, XA′)来描述。
二维局部均值随机场为:
E[XA(t1,t2)]=1A∫t1+T1/2t1−T1/2∫t2+T2/2t2−T2/2E[X(t1,t2)]dt1 dt2=μ。 (2) 方差为:
Var(XA)=σ2T1T2=σ2γ(T1,T2), (3) 式中: γ(T1, T2)为X(t1, t2)的方差折减函数,表示由局部平均引起的点方差σ2到局部平均方差σ2T1T2的折减程度。
γ(T1k,T2l)=1T1kT2l∫T1k−T1k∫T2l−T2l(1−|τ1|T1k)×(1−|τ2|T2l)ρ(τ1,τ2)dτ1 dτ2, (4) 式中: τ1和τ2分别为两点间t1和t2方向的坐标差; ρ(τ1, τ2) 为随机场的相关函数,其描述的是随机场的空间平均方差随平均范围增加而逐渐衰减的特性。
二维随机场常用的相关函数如表 1,式中δ1和δ2为随机场在t1和t2方向上的相关距离。
表 1 二维随机场常用的相关函数Table 1. Correlation functions for a 2D random field函数名称 相关函数 指数函数 ρ(τ1,τ2)=exp(−2√τ21δ21+τ22δ22) 高斯函数 ρ(τ1,τ2)=exp[−π(τ21δ21+τ22δ22)] 可分离的指数函数 ρ(τ1,τ2)=exp[−2(|τ1|δ1+|τ2|δ2)] 注:字母含义见正文 对于两矩形随机场单元A和A′,协方差Cov(XA, XA′)为:
Cov(XA,XA′)=1AA′σ243∑k=03∑l=0(−1)k+I.(T1kT2l)2γ(T1kT2l) (5) 式中: T1k和T2l(k=0, …, 3;l=0, …, 3)的约定如图 1所示;A和A′是2个矩形面积;在t1方向,定义T10为2个矩形内侧边界的最小间距,T11为2个矩形左侧边界间距,T12为2个矩形外侧边界最大间距,T13为2个矩形右侧边界间距;在t2方向,定义T20为2个矩形内侧边界的最小间距,T21为2个矩形下侧边界间距,T22为2个矩形外侧边界最大间距,T23为2个矩形上侧边界间距。
由式(5)求得局部平均随机场的均值、方差和协方差等数字特征之后,即可采用Gram-Schmidt特征正交化变换方法生成符合岩土参数空间变异性特征的随机场。
尽管各种相关函数在形式上有差别,但其对应的方差折减函数却相差不大,表明随机场的局部平局法具有对相关函数不敏感的特性。因此,在随机场的相关模型不明确的情况下,选用简单的相关函数形式是可取的[18],本研究选用可分离的指数函数为相关函数。
本研究采用强度折减法[14, 19-20]求解边坡稳定安全系数只对土壤参数进行折减,以塑性区贯通作为边坡整体失稳破坏的判据过于主观,因此本研究选择以边坡位移监测点处位移突变时对应的强度折减系数作为稳定安全系数[21-22]。
2. 土石混合体边坡模型
2.1 块石生成算法
目前,在二维块石生成算法方面,比较成熟的技术有以下3种:利用随机算法生成规则形状块石[23];利用数字图像处理技术生成真实形状块石[24];建立块石形状数据库,按照要求将块石导入模型[25]。本研究采用Python脚本方式,在有限元分析软件ABAQUS中建立块石模型,块石生成算法参考文献[26],具体思路如下:首先按照给定粒径范围生成相互不接触的圆形,将360°的圆心角随机进行15~20等分,即对应15~20条边,每条边对应每1个顶点,每个顶点到圆心的距离为半径乘以1个距离折减系数α,本研究α∈[0.3, 1],最后将顶点依次连线即可生成块石,生成过程如图 2所示。
2.2 土石混合体边坡模型建立
为了研究块石质量分数对土石混合体边坡稳定性影响,在有限元软件ABAQUS中,建立含有不同块石质量分数的土石混合体边坡模型,其几何尺寸如图 3所示,边坡坡度为45°,其中,O点为边坡的位移监测点。有限元计算模型如图 4所示,其中,块石尺寸应为0.1~0.4倍的边坡宽度[27-28],本模型边坡宽度为10 m,块石直径为1.05~1.07 m,图 4f为块石质量分数35%,块石直径3.15~3.21 m的土石混合体边坡。对边坡左、右侧施加水平位移约束,下侧施加固定端约束,重力加速度g取为10 m/s2。具体块石和土体参数见表 2。
表 2 土石混合体边坡模型参数Table 2. Parameters for the soil-rock mixture slope model组成成分 密度/(kg·m-3) 弹性模量/MPa 泊松比v 黏聚力c/kPa 内摩擦角φ/(°) 本构模型 块石 2 700 15 000 0.25 — — 理想弹性模型 土体 2 000 50 0.30 12.38 20 摩尔-库伦破坏准则 2.3 基于随机场的边坡稳定性计算方法
边坡的稳定性主要由土体的强度决定,因此,本研究将有效抗剪强度参数c和φ指定为考虑空间变异性的随机变量,土体的分布特征参数见表 3,黏聚力(c)和内摩擦角(φ)在t1方向的相关距离为3.00 m,在t2方向的相关距离为1.00 m。
表 3 有效抗剪强度参数的分布特征参数Table 3. Distribution characteristic parameters of effective shear strength parameters分布特征参数 黏聚力c/kPa 内摩擦角φ/(°) 均值μ 12.38 20.00 标准差σ 2.47 4.00 根据随机场局部平均法理论,采用MATLAB脚本语言和Python语言编制了岩土参数空间变异性随机场模拟程序,生成了有效抗剪强度参数c和φ的随机场,将其作为材料参数赋值给有限元软件中的计算模型,进行数值模拟,完成蒙特卡洛模拟计算。该程序分为以下步骤:①采用Python脚本语言在ABAQUS中生成土石混合体边坡模型并划分网格,获取土体网格和节点信息;②采用MATLAB语言编制的随机场模拟程序,根据网格和节点信息、相关距离、标准差和相关函数生成协方差矩阵;③由随机场的协方差矩阵和均值生成随机场向量;④根据网格单元的位置,将随机场向量作为材料参数映射到ABAQUS的有限元单元,如图 5所示为均质土坡(不含块石)和块石质量分数为35%的土石混合体边坡有效抗剪强度参数c和φ的分布示意图;⑤由MABLAB控制整个计算过程,完成蒙特卡洛模拟计算。本研究考虑兼顾计算效率和计算精度,每个块石质量分数下边坡进行1 000次蒙特卡洛模拟。
3. 结果分析
3.1 块石质量分数对边坡稳定性影响
本节对比了边坡稳定不考虑有效抗剪强度参数空间变异性计算结果(下文统称确定性计算结果)与考虑抗剪强度参数空间变异性计算结果的差异,分析了块石对边坡剪切带分布的阻碍作用与边坡稳定安全之间的关系,从而评估块石质量分数对边坡稳定性影响。
图 6为不考虑有效抗剪强度参数空间变异性的确定性计算结果图,由强度折减法原理可知,边坡位移突变点处对应的强度折减系数ω即为所求的稳定安全系数,如表 4确定性计算结果所示,随着块石质量分数的增加,土石混合体边坡的稳定安全系数呈现增大趋势,均质土坡和块石质量分数5%的土石
表 4 稳定安全系数Table 4. Stability safety factor块石质量分数/% 确定性计算结果 均值 标准差 最小值 最大值 0 0.997 1.005 0.033 0.917 1.138 5 0.997 1.008 0.058 0.910 1.127 15 1.005 1.012 0.082 0.917 1.137 25 1.061 1.079 0.062 0.915 1.151 35 1.064 1.095 0.033 1.000 1.153 #35 1.334 — — — — 注:块石质量分数0为均质土坡;块石质量分数#35%相较于其他情况,块石粒径增大;下同 混合体边坡稳定安全系数结果相同,而块石质量分数为15%时边坡稳定安全系数较均质土坡提高了0.008。这说明在块石直径不发生变化时,块石质量分数达到一定阈值,此时块石分布在边坡剪切带区域较为密集,块石有效阻碍了边坡剪切带的发展,进而提高了边坡稳定性。
块石的存在改变了剪切带传播路径,在土石混合体中,剪切带传播方式分为2种:分叉传播和绕石传播,如图 7所示。
由图 8确定性计算结果的塑性应变云图可知,随着土体强度参数降低,由于块石的阻碍效应,土石混合体边坡从坡脚到坡顶形成了1条扇形主剪切带和若干次剪切带(图中较粗的虚线段为主剪切带,较细的点线段为次剪切带);块石质量分数为5%的土石混合体边坡只存在1条主剪切带,块石并未对塑性区发展起到很好的阻碍作用,因此和均质土坡稳定性相比并无显著变化。均质土坡的剪切带存在于边坡浅层,随着块石质量分数增加,剪切带位置逐渐向边坡内部延伸,但是当块石质量分数达到35%时,剪切带位置又回到了浅层,这可能与块石的空间分布有关,有待进一步探究。
3.2 块石大小对边坡稳定性影响
相比于均质土坡,当块石质量分数达到35%时,土石混合体边坡的稳定安全系数上升了6.72%,块石对边坡稳定性提高效果并不显著,因此,在保持块石含量不变的情况下,将块石粒径增大,建立了如图 8f的土石混合体边坡,确定性计算结果显示,此边坡稳定安全系数为1.334,较均质土坡提高了33.8%。对比图 8e, f计算结果塑性应变云图可知,较大的块石对塑性区发展起到了很好的阻碍作用。由于块石与土体强度相差较大,当剪切带遇到坚硬块石,会绕过块石沿着土体与块石交界薄弱处发展。对于图 8f块石质量分数#35%边坡,较大块石使得边坡剪切带向边坡内部发展,从而形成了距离更长的剪切带,更好的发挥了土体的抗剪强度,因此边坡稳定性得到提高。
3.3 随机场对边坡稳定性影响
图 9是在考虑有效抗剪强度参数空间变异性的随机场计算结果统计图,可以看出稳定安全系数均符合正态分布,分布的均值和标准差见表 4,随着块石质量分数升高,稳定安全系数的均值呈上升趋势。可以看到,考虑有效抗剪强度参数空间变异性的边坡稳定性计算结果要比确定性计算结果高,如果按照确定性计算结果,边坡的设计会趋于保守,增加不必要的成本。
图 10为每个土石混合体边坡1 000次模拟的结果里最大和最小稳定安全系数对应的塑性应变云图,通过对比可知,块石含量、分布,荷载及边界条件均相同的条件下,土体抗剪强度参数分布的空间变异性影响着边坡剪切带分布及边坡稳定安全。图 10a中边坡的剪切带发育充分,剪切带穿过区域土体强度较高,而图 10b中边坡剪切带还未得到充分发育,边坡已经破坏并产生大变形,导致计算不收敛,说明剪切带穿过区域的土体强度较低。这进一步证明土体抗剪强度参数的空间变异性对边坡稳定有较大影响,不可以忽略不计。即在其他条件相同的情况下,剪切带穿过区域的土体强度越高,边坡更趋于稳定,反之,边坡更容易产生破坏。
4. 结论
(1) 随块石含量升高,边坡稳定性逐渐升高,但稳定安全系数增幅较小,同时,边坡由浅层破坏逐渐发展为深层破坏,破坏模式与块石的分布情况有关。
(2) 对比相同块石含量下,不同块石大小对边坡稳定性影响结果可知,较大的块石对塑性区发展起到了很好的阻碍作用,显著提高了边坡的抗滑力。
(3) 考虑有效抗剪强度参数空间变异性的边坡稳定性计算结果要比确定性计算结果高,若将确定性计算结果作为设计依据,边坡设计会趋于保守,提高边坡加固成本[29]。
(4) 由塑性应变云图可知,剪切带分布区域主要为土石混合体边坡坡脚至坡顶这一局部区域,坡脚以下的块石是否会对土石混合体边坡稳定性产生影响还需要进一步探究,此外,本研究为二维工况下的土石混合体边坡稳定性分析,为了更加贴近工程实际,有必要开发三维土石混合体模型。
所有作者声明不存在利益冲突。 -
表 1 二维随机场常用的相关函数
Table 1. Correlation functions for a 2D random field
函数名称 相关函数 指数函数 ρ(τ1,τ2)=exp(−2√τ21δ21+τ22δ22) 高斯函数 ρ(τ1,τ2)=exp[−π(τ21δ21+τ22δ22)] 可分离的指数函数 ρ(τ1,τ2)=exp[−2(|τ1|δ1+|τ2|δ2)] 注:字母含义见正文 表 2 土石混合体边坡模型参数
Table 2. Parameters for the soil-rock mixture slope model
组成成分 密度/(kg·m-3) 弹性模量/MPa 泊松比v 黏聚力c/kPa 内摩擦角φ/(°) 本构模型 块石 2 700 15 000 0.25 — — 理想弹性模型 土体 2 000 50 0.30 12.38 20 摩尔-库伦破坏准则 表 3 有效抗剪强度参数的分布特征参数
Table 3. Distribution characteristic parameters of effective shear strength parameters
分布特征参数 黏聚力c/kPa 内摩擦角φ/(°) 均值μ 12.38 20.00 标准差σ 2.47 4.00 表 4 稳定安全系数
Table 4. Stability safety factor
块石质量分数/% 确定性计算结果 均值 标准差 最小值 最大值 0 0.997 1.005 0.033 0.917 1.138 5 0.997 1.008 0.058 0.910 1.127 15 1.005 1.012 0.082 0.917 1.137 25 1.061 1.079 0.062 0.915 1.151 35 1.064 1.095 0.033 1.000 1.153 #35 1.334 — — — — 注:块石质量分数0为均质土坡;块石质量分数#35%相较于其他情况,块石粒径增大;下同 -
[1] 徐文杰, 胡瑞林. 土石混合体概念、分类及意义[J]. 水文地质工程地质, 2009, 36(4): 50-56.XU W J, HU R L. Conception, classification and significations of soil-rock mixture[J]. Hydrogeology & Engineering Geology, 2009, 36(4): 50-56. (in Chinese with English abstract) [2] WANG S, WU W, CUI D S. On mechanical behaviour of clastic soils: Numerical simulations and constitutive modelling[J]. Géotechnique, 2022, 72(8): 706-721. doi: 10.1680/jgeot.20.P.184 [3] 林成远, 唐辉明, 汪丁建, 等. 块石定向性特征对土-石混合体强度影响的数值模拟[J]. 地质科技通报, 2020, 39(5): 38-46. doi: 10.19509/j.cnki.dzkq.2020.0519LIN C Y, TANG H M, WANG D J, et al. Influence on the strength of soil-rock mixture made by the rock block orientation features based on numerical experiment[J]. Bulletin of Geological Science and Technology, 2020, 39(5): 38-46. (in Chinese with English abstract) doi: 10.19509/j.cnki.dzkq.2020.0519 [4] LI C Z, CHEN G Q, GUO L X, et al. Slope stability and post-failure analysis of soil-rock-mixture using the modified 2D DDA-SPH method[J]. International Journal of Rock Mechanics and Mining Sciences, 2022, 157: 105170. doi: 10.1016/j.ijrmms.2022.105170 [5] LIU S Q, CHENG Z C, QU X. A novel method incorporating digital image processing and A-star algorithm for soil-rock mixture (SRM)slope stability[J]. Environmental Earth Sciences, 2023, 82(7): 163. doi: 10.1007/s12665-023-10866-7 [6] 唐军峰, 唐雪梅, 肖鹏, 等. 库水位升降与降雨作用下大型滑坡体渗流稳定性分析[J]. 地质科技通报, 2021, 40(4): 153-161. doi: 10.19509/j.cnki.dzkq.2021.0409TANG J F, TANG X M, XIAO P, et al. Analysis of seepage stability of large-scale landslide under water-level fluctuation and rainfall[J]. Bulletin of Geological Science and Technology, 2021, 40(4): 153-161. (in Chinese with English abstract) doi: 10.19509/j.cnki.dzkq.2021.0409 [7] 邵帅, 季顺迎. 块石空间分布对土石混合体边坡稳定性的影响[J]. 工程力学, 2014, 31(2): 177-183.SHAO S, JI S Y. Effects of rock spatial distributions on stability of rock-soil-mixture slope[J]. Engineering Mechanics, 2014, 31(2): 177-183. (in Chinese with English abstract) [8] GAO W W, YANG H R, TENG J B, et al. Numerical study of the influence of the spatial distribution of oversized rock blocks on the stability of soil-rock mixture slopes[J]. Advances in Civil Engineering, 2023, 2023: 4575118. [9] 蒋水华, 李典庆, 周创兵, 等. 考虑自相关函数影响的边坡可靠度分析[J]. 岩土工程学报, 2014, 36(3): 508-518.JIANG S H, LI D Q, ZHOU C B, et al. Slope reliability analysis considering effect of autocorrelation functions[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(3): 508-518. (in Chinese with English abstract) [10] 李典庆, 祁小辉, 周创兵, 等. 考虑参数空间变异性的无限长边坡可靠度分析[J]. 岩土工程学报, 2013, 35(10): 1799-1806.LI D Q, QI X H, ZHOU C B, et al. Reliability analysis of infinite soil slopes considering spatial variability of soil parameters[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(10): 1799-1806. (in Chinese with English abstract) [11] WANG M X, TANG X S, LI D Q, et al. Subset simulation for efficient slope reliability analysis involving copula-based cross-correlated random fields[J]. Computers and Geotechnics, 2020, 118: 103326. doi: 10.1016/j.compgeo.2019.103326 [12] VANMARCKE E H. Probabilistic modeling of soil profiles[J]. Journal of the Geotechnical Engineering Division, 1977, 103(11): 1227-1246. doi: 10.1061/AJGEB6.0000517 [13] 薛亚东, 方超, 葛嘉诚. 各向异性随机场下的边坡可靠度分析[J]. 岩土工程学报, 2013, 35(增刊2): 77-82.XUE Y D, FANG C, GE J C. Slope reliability in anisotropic random fields[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(S2): 77-82. (in Chinese with English abstract) [14] LIU Y, XIAO H W, YAO K, et al. Rock-soil slope stability analysis by two-phase random media and finite elements[J]. Geoscience Frontiers, 2018, 9(6): 1649-1655. doi: 10.1016/j.gsf.2017.10.007 [15] KASAMA K, WHITTLE A J. Effect of spatial variability on the slope stability using Random Field Numerical Limit Analyses[J]. Georisk: Assessment and Management of Risk for Engineered Systems and Geohazards, 2016, 10(1): 42-54. doi: 10.1080/17499518.2015.1077973 [16] 何成, 唐辉明, 申培武, 等. 应变软化边坡渐进破坏模式及稳定性可靠度[J]. 地球科学, 2021, 46(2): 697-707.HE C, TANG H M, SHEN P W, et al. Progressive failure mode and stability reliability of strain-softening slope[J]. Earth Science, 2021, 46(2): 697-707. (in Chinese with English abstract) [17] 杨杰, 陈佳莹, 尹振宇. 砂土初始孔隙比空间变异性对其力学特性及破坏模式的影响分析[J]. 计算力学学报, 2021, 38(6): 754-762.YANG J, CHEN J Y, YIN Z Y. Influence of spatial variability of initial void ratio on the mechanical response and failure mode of sand[J]. Chinese Journal of Computational Mechanics, 2021, 38(6): 754-762. (in Chinese with English abstract) [18] LI K S, LUMB P. Probabilistic design of slopes[J]. Canadian Geotechnical Journal, 1987, 24(4): 520-535. [19] 郑颖人, 赵尚毅, 张鲁渝. 用有限元强度折减法进行边坡稳定分析[J]. 中国工程科学, 2002, 4(10): 57-61.ZHENG Y R, ZHAO S Y, ZHANG L Y. Slope stability analysis by strength reduction FEM[J]. Engineering Science, 2002, 4(10): 57-61. (in Chinese with English abstract) [20] 咸玉建, 陈学军, 汪志刚, 等. 基于有限元强度折减系数法岩土边坡稳定性分析与抗滑桩设计[J]. 地质科技情报, 2015, 34(4): 176-182.XIAN Y J, CHEN X J, WANG Z G, et al. Geotechnical slope stability analysis and the design of anti-slide pile based on strength reduction FEM[J]. Geological Science and Technology Information, 2015, 34(4): 176-182. (in Chinese with English abstract) [21] ZHAO T, SUN J Z, ZHANG B, et al. Analysis of slope stability with dynamic overloading from earthquake[J]. Journal of Earth Science, 2012, 23(3): 285-296. [22] 蒋先平, 张鹏, 卢艺伟, 等. 物质点强度折减法边坡失稳判据选择方法[J]. 地质科技通报, 2022, 41(2): 113-122. doi: 10.19509/j.cnki.dzkq.2021.0075JIANG X P, ZHANG P, LU Y W, et al. Slope failure criterion for the strength reduction material point method[J]. Bulletin of Geological Science and Technology, 2022, 41(2): 113-122. (in Chinese with English abstract) doi: 10.19509/j.cnki.dzkq.2021.0075 [23] LI X, LIAO Q L, HE J M. In situ tests and a stochastic structural model of rock and soil aggregate in the Three Gorges Reservoir area, China[J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(3): 702-707. [24] 徐文杰, 胡瑞林, 岳中琦, 等. 土石混合体细观结构及力学特性数值模拟研究[J]. 岩石力学与工程学报, 2007, 26(2): 300-311.XU W J, HU R L, YUE Z Q, et al. Mesostructural character and numerical simulation of mechanical properties of soil-rock mixtures[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(2): 300-311. (in Chinese with English abstract) [25] 赵鑫曜, 陈建功, 张海权, 等. 空间与颗粒参数对块石土边坡稳定性影响数值模拟研究[J]. 工程科学与技术, 2020, 52(4): 166-175.ZHAO X Y, CHEN J G, ZHANG H Q, et al. Analysis of space and particle parameter influence to stability of soil-rock mixture slope[J]. Advanced Engineering Sciences, 2020, 52(4): 166-175. (in Chinese with English abstract) [26] 李亮, 李彦军, 赵炼恒, 等. 土石混合边坡块石随机生成方法与稳定性分析[J]. 湖南大学学报(自然科学版), 2017, 44(7): 170-178.LI L, LI Y J, ZHAO L H, et al. Method for generating random soil-rock mixed slope and stability analysis[J]. Journal of Hunan University(Natural Sciences), 2017, 44(7): 170-178. (in Chinese with English abstract) [27] LIU S Q, HUANG X W, ZHOU A Z, et al. Soil-rock slope stability analysis by considering the nonuniformity of rocks[J]. Mathematical Problems in Engineering, 2018, 2018(1): 3121604. [28] NAPOLI M L, BARBERO M, RAVERA E, et al. A stochastic approach to slope stability analysis in bimrocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2018, 101: 41-49. [29] 胡新丽, 唐辉明, 陈建平. 高速公路顺层路堑边坡优化设计方法[J]. 地球科学(中国地质大学学报), 2001, 26(4): 373-376.HU X L, TANG H M, CHEN J P. Optimum design method for bedding slope of expressway[J]. Earth Science(Journal of China University of Geosciences), 2001, 26(4): 373-376. (in Chinese with English abstract) -