为了研究边坡安全系数的影响因素,考虑边坡土体粘聚力、内摩擦角、体积模量、密度、抗拉强度及边坡倾角,采用FLAC3D有限差分软件计算边坡安全系数,在分析模型网格疏密程度对计算结果影响规律的基础上,运用正交试验法对计算结果分析,研究影响边坡安全系数的各因素敏感性。研究结果表明,内聚力、内摩擦角、密度及边坡倾角对边坡安全系数的影响非常显著,敏感性影响程度从大到小依次为:内聚力、内摩擦角、密度、边坡倾角,且内聚力、内摩擦角与安全系数成正相关,边坡倾角、密度与安全系数成负相关;体积模量与抗拉强度对边坡安全系数的影响不显著。经过全长锚固支护后的边坡安全系数从0.96提高到1.22,锚杆支护在提高边坡稳定性方面有较好的效果。 In order to study the influencing factors of slope safety factor, considering the cohesive force of soil mass, internal friction angle, volume modulus, density, tensile strength and slope inclination, FLAC3Dfinite difference software was used to calculate the slope safety factor. On the basis of analyzing the influence law of model mesh density on the calculation results, orthogonal test method was used to analyze the calculation results. The sensitivity of factors affecting the safety factor of slope is studied. The results show that cohesion, internal friction angle, density and slope inclination have significant effects on slope safety factor, and the degree of sensitivity is as follows from large to small: cohesion, internal friction angle, density and slope inclination, and cohesion and internal friction angle are positively correlated with safety factor, while slope inclination and density are negatively correlated with safety factor. The influence of volume modulus and tensile strength on slope safety factor is not significant. The safety factor of the slope increased from 0.96 to 1.22 after the full length anchorage support, and the bolt support has a good effect on improving the stability of the slope.
为了研究边坡安全系数的影响因素,考虑边坡土体粘聚力、内摩擦角、体积模量、密度、抗拉强度及边坡倾角,采用FLAC3D有限差分软件计算边坡安全系数,在分析模型网格疏密程度对计算结果影响规律的基础上,运用正交试验法对计算结果分析,研究影响边坡安全系数的各因素敏感性。研究结果表明,内聚力、内摩擦角、密度及边坡倾角对边坡安全系数的影响非常显著,敏感性影响程度从大到小依次为:内聚力、内摩擦角、密度、边坡倾角,且内聚力、内摩擦角与安全系数成正相关,边坡倾角、密度与安全系数成负相关;体积模量与抗拉强度对边坡安全系数的影响不显著。经过全长锚固支护后的边坡安全系数从0.96提高到1.22,锚杆支护在提高边坡稳定性方面有较好的效果。
边坡稳定,安全系数,正交试验,敏感性分析,锚杆支护
Jiale Hu1, Yang Hu1,2, Wei Jiang1,2, Jian Xie3, Hubin Xia1
1School of Resources and Civil Engineering, Suzhou University, Suzhou Anhui
2National Coal Mine Water Hazard Prevention Engineering Technology Research Center (Suzhou University), Suzhou Anhui
3Wugou Coal Mine, Anhui Hengyuan Coal Power Group Co., Ltd., Huaibei Anhui
Received: Apr. 1st, 2024; accepted: Apr. 21st, 2024; published: Apr. 30th, 2024
In order to study the influencing factors of slope safety factor, considering the cohesive force of soil mass, internal friction angle, volume modulus, density, tensile strength and slope inclination, FLAC3Dfinite difference software was used to calculate the slope safety factor. On the basis of analyzing the influence law of model mesh density on the calculation results, orthogonal test method was used to analyze the calculation results. The sensitivity of factors affecting the safety factor of slope is studied. The results show that cohesion, internal friction angle, density and slope inclination have significant effects on slope safety factor, and the degree of sensitivity is as follows from large to small: cohesion, internal friction angle, density and slope inclination, and cohesion and internal friction angle are positively correlated with safety factor, while slope inclination and density are negatively correlated with safety factor. The influence of volume modulus and tensile strength on slope safety factor is not significant. The safety factor of the slope increased from 0.96 to 1.22 after the full length anchorage support, and the bolt support has a good effect on improving the stability of the slope.
Keywords:Slope Stability, Safety Factor, Orthogonal Test, Sensitivity Analysis, Bolt Support
Copyright © 2024 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
边坡稳定性是指边坡在重力或其他外力作用下保持稳定的能力。边坡的稳定性除与边坡的几何构造、岩土体物理力学性质(内聚力、内摩擦角、密度、弹性模量等)等内部因素直接相关 [
强度折减法中边坡稳定的安全系数定义为:使边坡刚好达到临界破坏状态时,对岩土体的抗剪强度进行折减的程度,即定义安全系数为岩土体的实际抗剪强度与临界破坏时的折减后剪切强度的比值。强度折减法即采用式(1)、(2)来调整岩土体的强度指标c和φ
C F = C / F trial (1)
φ F = tan − 1 [ ( tan φ ) / F trial ] (2)
式中:c——折减前的粘结力,(Pa);
φ——折减前的摩擦角,(˚);
cF——折减后的粘结力,(Pa);
φF——折减后的摩擦角,(˚);
Ftrial——折减系数。
然后对边坡稳定性进行数值分析,不断地增加折减系数,反复计算,直至边坡发生临界破坏状态,此时得到的折减系数Ftrial即为安全系数FS。
假设边坡模型为各向同性的均质弹塑性材料,在求解安全系数中以最大不平衡力比作为计算终止标准,即最大不平衡力比R ≤ 10−5即达到数值计算平衡状态。模型采用摩尔–库伦本构模型,上表面为自由边界,底部为固定位移边界,(对x、y、z三个方向进行位移约束),模型的四周为法向位移约束(按照其所处位置分别对x、y两个方向进行位移约束)。
数值模型中网格划分采用六面块体网格(图1所示),笔者在对边坡安全系数的数值计算过程中发现,网格的疏密程度对计算结果具有一定的影响,故对边坡安全系数影响因素的敏感性分析前,首先进行了网格疏密程度对安全系数计算结果的影响规律分析。模型的网格数量按照基础网格数量划分(238个)、加密一倍网格数量(476个)、稀疏一倍网格数量(119个)三种进行选取,数值模型大小保持一致,所选取的边坡体物理力学参数如表1所示。
图1. 数值模型网格划分
黏聚力/Pa | 体积模量/m3 | 抗拉强度/Pa | 边坡倾角/(˚) | 内摩擦角/(˚) | 密度/(kg/m3) | 抗剪强度/Pa |
---|---|---|---|---|---|---|
12000.0 | 1.0e8 | 1.0e8 | 45 | 20 | 2000 | 3.0e7 |
表1. 边坡体物理力学参数
不同网格数量条件下所得剪切应变增量及安全系数云图如图2所示。由图2可知,在数值模型大小及其余条件保持不变的情况下,当模型网格数量为119、238、476时,对应安全系数结果分别为:1.07、1.03、1.01,即随着网格数量的增大,所得安全系数结果呈现减小的趋势,误差率分别为:3.7%、1.9%,误差率随着网格数量的增大也呈现不断减小的趋势,因此从该点上看,网格数量越大其所得计算精度越高,但网格数量的增加会大大增加数值计算所用的时间,本实例中当模型网格数量为119、238、476时,计算所用时间分别为1 min、3 min、23 min,因此合理选择网格疏密程度对提高计算效率与减小计算误差意义重大,本算例中按照基础网格数量划分为238个与加密一倍网格数量476个,所得结果误差仅为1.9%,而时间多耗费近10倍,故在下述计算过程中,数值模型的网格按照基础网格数量模型进行划分。
(a) 网格加密一倍
图2. 不同网格数量模型剪切应变增量与安全系数
正交试验设计常用于分析涉及多因素的材料、结构等,该试验方法已被证明十分有效 [
选取黏聚力、内摩擦角、体积模量、边坡倾角、密度以及抗拉强度这6个参数作为敏感性因素。参数具体取值如表2所示。
水平 | 黏聚力/kPa | 体积模量/m3 | 抗拉强度/Pa | 边坡倾角/(˚) | 内摩擦角/(˚) | 密度/(kg/m3) |
---|---|---|---|---|---|---|
1 | 6.00 | 1.0e7 | 0.5e8 | 38 | 16 | 1500 |
2 | 12.00 | 1.0e8 | 1.0e8 | 45 | 20 | 2000 |
3 | 24.00 | 1.0e9 | 1.5e8 | 51 | 24 | 2500 |
表2. 边坡因素水平
选择正交表为L18(37),坡体稳定性正交试验方案见表3,正交试验计算结果见表4。其中K1、K2、K3分别为任一列上各因素水平为1、2、3时各自对应的试验结果之和, G = ∑ k = 1 18 y k ;C为校正数, C = G 2 / 18 ;Si为各因素离差平方和, S i = ( k 1 2 + k 2 2 + k 3 2 ) / 6 − C ;ST为总离差平方和, S T = ∑ k = 1 18 y k 2 − C ;Se为试验误差的离差平方和, S e = S T - S i ,i = 1, 2, …, 6。由于极差的值看不出显著性,故引入方差来判断各影响因素的显著性。利用正交设计助手软件,取显著性水平α = 0.05,计算结果见表4。表中各主控因素的平均偏差平方和与误差的平均偏差平方和之比为F,采用F检验法将各主控因素的F值与F表进行对比分析,进而确定各主控因素对底板破坏深度影响程度大小,其计算公式为
F = S i / f i S e / f e (3)
式中:Si——主控因素偏方和;
fi——自由度;
Se——试验误差偏差平方和;
fe——为试验误差自由度。
将表4中相应数据代入式中,可得主控因素的F值,如表5所示。
由表4所得结果可知;黏聚力、边坡倾角、内摩擦角、密度对边坡安全系数的影响为非常显著,体积模量和抗拉强度为不显著。根据敏感性从大到小依次为黏聚力、内摩擦角、密度、边坡倾角、体积模量、抗拉强度。因此在实际施工或治理边坡稳定性方面,要优先着重分析边坡岩土体的黏聚力大小,影响岩土体的黏聚力因素有颗粒组成(土类别)、颗粒大小、颗粒形状、含水量等,所以对于实际应用中的边坡,优先选黏聚力大的土类,最优含水率的选取;有时仅仅考虑提高岩土体黏聚力,比较有限,特别对于自然边坡,要做好排水措施,辅助增设人工防护措施,如增设锚杆、锚索等,来提高岩土体的黏聚力与内摩擦角等因素,减少、防范边坡的坍塌。
试验编号 | 黏聚力/kPa | 体积模量/m3 | 抗拉强度/Pa | 边坡倾角/(˚) | 内摩擦角/(˚) | 密度/(kg/m3) | 空列 | 安全系数 |
---|---|---|---|---|---|---|---|---|
1 | 1 (6) | 1 (1.0e7) | 1 (0.5e8) | 1 (38) | 1 (16) | 1 (1500) | 1 | 0.87 |
2 | 1 | 2 (1.0e8) | 2 (1.0e8) | 2 (45) | 2 (20) | 2 (2000) | 2 | 0.78 |
3 | 1 | 3 (1.0e9) | 3 (1.5e8) | 3 (51) | 3 (24) | 3 (2500) | 3 | 0.74 |
4 | 2 (12) | 1 | 1 | 2 | 2 | 3 | 3 | 0.94 |
5 | 2 | 2 | 2 | 3 | 3 | 1 | 1 | 1.19 |
6 | 2 | 3 | 3 | 1 | 1 | 2 | 2 | 1.04 |
7 | 3 (24) | 1 | 2 | 1 | 3 | 2 | 3 | 1.81 |
8 | 3 | 2 | 3 | 2 | 1 | 3 | 1 | 1.17 |
9 | 3 | 3 | 1 | 3 | 2 | 1 | 2 | 1.61 |
10 | 1 | 1 | 3 | 3 | 2 | 2 | 1 | 0.70 |
11 | 1 | 2 | 1 | 1 | 3 | 3 | 2 | 0.99 |
12 | 1 | 3 | 2 | 2 | 1 | 1 | 3 | 0.76 |
13 | 2 | 1 | 2 | 3 | 1 | 3 | 2 | 0.74 |
14 | 2 | 2 | 3 | 1 | 2 | 1 | 3 | 1.35 |
15 | 2 | 3 | 1 | 2 | 3 | 2 | 1 | 1.16 |
16 | 3 | 1 | 3 | 2 | 3 | 1 | 2 | 1.89 |
17 | 3 | 2 | 1 | 3 | 1 | 2 | 3 | 1.23 |
18 | 3 | 3 | 2 | 1 | 2 | 3 | 1 | 1.47 |
表3. 正交试验方案
试验编号 | 黏聚力/kPa | 体积模量/m3 | 抗拉强度/Pa | 边坡倾角/(˚) | 内摩擦角/(˚) | 密度/(kg/m3) | 空列 | 安全系数 |
---|---|---|---|---|---|---|---|---|
K1 | 4.84 | 6.95 | 6.80 | 7.53 | 5.81 | 7.67 | 6.56 | G = 20.44 G2= 417.79 C = 23.21 |
K2 | 6.42 | 6.71 | 6.75 | 6.70 | 6.85 | 6.72 | 7.05 | |
K3 | 9.18 | 6.78 | 6.89 | 6.21 | 7.78 | 6.05 | 6.83 | |
K 1 2 | 23.43 | 48.30 | 46.24 | 56.70 | 33.76 | 58.83 | 43.03 | |
K 2 2 | 41.22 | 45.02 | 45.56 | 44.89 | 46.92 | 45.16 | 49.70 | |
K 3 2 | 84.27 | 45.97 | 47.47 | 38.56 | 60.53 | 36.60 | 46.65 | |
Si | 1.61 | 0.005 | 0.002 | 0.148 | 0.324 | 0.222 | 0.020 | |
ST | 2.33 | |||||||
Se | 0.02 |
表4. 正交试验结果
来源 | 偏差平方和 | 自由度 | 均偏差平方和 | F | F临界值 | 显著性 |
---|---|---|---|---|---|---|
黏聚力 | 1.61 | 2 | 0.805 | 178.89 | F0.05(2, 6) | ** |
体积模量 | 0.005 | 2 | 0.0025 | 0.556 | F0.05(2, 6) | |
抗拉强度 | 0.002 | 2 | 0.001 | 0.222 | F0.05(2, 6) | |
边坡倾角 | 0.148 | 2 | 0.074 | 16.44 | F0.05(2, 6) | ** |
内摩擦角 | 0.324 | 2 | 0.162 | 36.00 | F0.05(2, 6) | ** |
密度 | 0.222 | 2 | 0.111 | 24.67 | F0.05(2, 6) | ** |
误差e | 0.020 | 2 | 0.01 | |||
误差 | 0.027 | 6 | 0.0045 |
表5. 方差分析表
注:① F0.01(2, 6) = 10.93;F0.05(2, 6) = 5.14;② 当Si< 2Se时,可将因素或交互作用的偏差平方和、自由度并入误差的偏差平方和、自由度,记为误差eΔ,使得误差的偏差平方和、自由度增大,提高了F检验的灵敏度;③ F > F0.01(fi, fe),该因素高度显著,以**表示;④ F0.05(fi, fe) < F < F0.01(fi, fe),该因素显著,以*表示;⑤ F < F0.05(fi, fe),该因素不显著。
由上述试验结果可得,黏聚力、边坡倾角、内摩擦角、密度这四种影响因素对边坡安全系数的影响非常显著,文章在此进一步讨论这四个影响因素对边坡安全系数变化的影响规律。边坡安全系数随岩土体密度、边坡倾角变化规律分别如图3、图4所示。
图3. 边坡安全系数与岩土体密度关系
由图3可知,在其余条件保持不变的前提下,当岩土体的密度为1000 kg/m3时,所得边坡安全系数为2.11;当岩土体的密度取值为1500 kg/m3,求得边坡安全系数为1.60;当密度取值为2000 kg/m3时,求得边坡安全系数为1.33,随着密度的不断增大,边坡安全系数在不断减小,且减小率不断减小,以上结果表明岩土体密度与边坡安全系数成负相关关系,且当岩土体密度增大到一定程度后,这种负相关关系逐渐降低。
图4. 边坡安全系数与边坡倾角关系
由图4可知,在其余条件保持不变的情况下,当边坡倾角为38˚时,求得安全系数为1.04;当边坡倾角为40˚时,安全系数为0.99;当边坡倾角为42˚时,安全系数为0.96。即随着边坡倾角的增大,安全系数不断减小,边坡倾角与安全系数近似呈现线性负相关关系,边坡倾角越大,安全系数越小,即边坡越缓,边坡稳定性越强。
内聚力、内摩擦角作为土体抗剪强度中两个重要参数,其对边坡安全系数的影响规律如图5所示。
图5. 边坡安全系数与内聚力、内摩擦角关系
根据图5可知,当内聚力为6 Pa、内摩擦角取值为16˚时,边坡安全系数为0.67;当黏聚力取值为8 Pa,内摩擦角取值为17˚时,边坡安全系数为0.79;当黏聚力取值为10 Pa,内摩擦角取值为18˚时,边坡安全系数为0.89,即随着黏聚力与内摩擦角的不断增大,边坡安全系数也不断增大,即边坡岩土体的内聚力、内摩擦角与边坡安全系数成正相关关系。在其余条件保持不变的情况下,岩土体的内聚力、内摩擦角越大,边坡安全系数越大,边坡越稳定。
工程实践表明 [
采用全锚法锚杆支护加固边坡,边坡倾角为45˚,锚杆杆体采用直径为8 mm的钢筋,锚杆长度为取4 m,锚杆间距为1 m,锚固剂单位长度上的刚度为1.0e8,锚固剂单位长度上的黏聚力为10e6Pa。锚杆支护前后边坡剪切应变增量与安全系数如图6所示。由图6可知,天然边坡在未支护状态时安全系数为0.96,经过全锚支护后的安全系数提高至1.22,通过全锚支护边坡得到了加固。
图6. 锚杆支护前后边坡剪切应变增量及安全系数
(1) 网格疏密程度对边坡安全系数数值计算结果有着较大影响,基于FLAC计算边坡安全系数需要根据实际情况合理选择网格疏密程度。
(2) 边坡安全系数影响因素敏感性从大到小依次为黏聚力、内摩擦角、密度、边坡倾角、体积模量、抗拉强度。
(3) 黏聚力、内摩擦角、密度、边坡倾角对边坡安全系数的影响为非常显著,体积模量与抗拉强度为不显著。黏聚力与内摩擦角和边坡安全系数成正相关,密度与边坡倾角和边坡安全系数成负相关。
(4) 天然状态下的边坡经过锚杆支护后安全系数从0.96增加到1.22,锚杆支护在提高边坡稳定性方面有较好的效果。
宿州学院博士启动金项目(2021BSK014);宿州学院博士后启动金项目(2022BSH003);宿州学院校级重点项目(2021yzd04);安徽省大学生创新创业训练计划项目(S202110379198)。
扈家乐,胡 洋,姜 玮,解 建,夏胡斌. 基于正交试验的边坡安全系数影响因素敏感性分析Sensitivity Analysis of Influencing Factors of Slope Safety Factor Based on Orthogonal Test[J]. 土木工程, 2024, 13(04): 490-499. https://doi.org/10.12677/hjce.2024.134054
http://kns.cnki.net/kcms/detail/34.1055.TD.20211118.2010.002.html, 2022-09-03.
https://doi.org/10.27014/d.cnki.gdnau.2022.000425
https://doi.org/10.16031/j.cnki.issn.1003-8035.2020.03.05
https://doi.org/10.13301/j.cnki.ct.2021.01.017
https://doi.org/10.1016/j.ijprt.2017.03.008
https://doi.org/10.13646/j.cnki.42-1395/u.2023.05.032