Reliability analysis of soil nailing-reinforced slopes considering fuzzy randomness
-
摘要: 针对土钉加固边坡未考虑边坡模糊随机性问题,提出考虑土体参数模糊随机性与边坡模糊过渡区间的土钉加固边坡可靠度计算方法。将土体样本试验力学参数随机值转换为模糊随机变量;在模糊随机变量与模糊过渡区间的基础上,建立模糊随机极限状态方程;确定土体参数不同λ-截集上下限及与之对应的边坡滑裂面;推导出仅与土体参数变量相关的土钉加固边坡可靠度计算公式,得到土钉加固边坡模糊随机可靠度。最后与传统蒙特卡洛模拟法计算出的可靠度进行对比得出:应用蒙特卡洛模拟计算出的失效概率为零时,考虑边坡模糊性后边坡失效概率为5.9%,表明考虑边坡模糊性与模糊过渡区间的加固边坡可靠度分析更能反映边坡实际状态。Abstract: With regard to the problem of no cosideration of the fuzzy randomness of slopes reinforced by soil nailing, a method for calculating the reliability of soil nailing-reinforced slopes is proposed considering the fuzzy randomness of soil parameters and the fuzzy transition interval of the slopes.Firstly, the random values of mechanical parameters for soil sample tests are converted into the fuzzy random variables.Based on the fuzzy random variables and the fuzzy transition interval, an equation for the fuzzy random limit state is established.The upper and lower limits of the λ-cut set and the corresponding slip surface of slopes are determined.Then, the formula for calculating the reliability of soil nailing-reinforced slopes only related to the soil parameter variables is derived, and the fuzzy random reliability of soil nailing-reinforced slopes is obtained.Finally, the calculated results are compared with those calculated by the traditional Monte Carlo simulation method.It is shown that the reliability analysis of reinforced slopes considering fuzzy transition interval and fuzzy randomness of soil parameters can better reflect the actual state of reinforced slopes.
-
Keywords:
- soil parameter /
- soil nailing /
- fuzzy randomness /
- fuzzy transition interval /
- reliability
-
0. 引言
边坡安全系数通常用于判断边坡是否稳定。但实际工程中,存在安全系数满足要求的边坡发生滑塌,也存在安全系数未达到要求的边坡保持稳定,其主要原因在于边坡工程涉及大量不确定性因素。影响边坡稳定的因素主要为土体物理力学参数,而土体力学参数既有取样时的随机性,也有自身存在的模糊性。另一方面,边坡由稳定到不稳定应当是渐变的,其边界是模糊的,也就是说边坡的稳定与不稳定之间存在具有模糊性的过渡区间。
针对边坡中存在的不确定性,熊文林等[1]提出基于现有的岩石样本试验参数值,在考虑参数的模糊-随机性前提下,推导出岩石力学参数样本值的模糊随机统计公式。贾厚华等[2]在边坡稳定性评价中考虑边坡存在的随机性与模糊性,并以模糊随机变量为基本变量,提出模糊随机功能函数和模糊随机极限方程。吕玺琳等[3]在一次二阶矩法的基础上,推导了仅考虑主要变量c, φ情况下,边坡模糊随机可靠分析公式。王宇等[4]将统计矩点估计法与模糊随机理论相结合,提出边坡工程模糊随机可靠性分析的模糊点估计法。El-Ramly等[5]采用蒙特卡罗模拟方法介绍边坡可靠度分析方法的实际应用。Wong[6]采有限元法结合响应面法对边坡破坏概率进行研究,并与离心机试验结果进行对比。Griffiths等[7]采用蒙特卡罗随机有限元法对简单边坡稳定性进行可靠度分析,研究土体强度参数空间变异性影响。
但以上边坡模糊随机可靠度分析方法均针对未加固边坡,很少对加固边坡提出模糊随机可靠度分析方法[8]。近年来,由于土钉具有质量轻、经济可靠、施工简便等优点,在公路、铁路等建设中,常采用土钉支护形式进行边坡加固。本文在考虑土体参数模糊随机性与边坡模糊过渡区间的前提下,提出土钉加固边坡可靠度分析方法,并推导出土钉加固边坡可靠度计算公式。
1. 模糊随机变量
在边坡工程中,影响稳定性的主要因素为土体物理力学参数,其中主要为土体抗剪强度,即c, φ值。利用式(1)~(8),将由土体样本试验得到的土体参数c, φ随机值转化为模糊随机变量。
(1) (2) 式中,
(3) (4) 式(1), (2)中djmax, djmin分别为dji(j=1, 2;i=1, 2, …, n)中的最大值、最小值。
(5) (6) 式中,
(7) (8) 式(5), (6)中tjmax, tjmin分别为tji(j=1, 2;i=1, 2, …, n)中的最大值、最小值。
式(1)~(8)进行迭代运算后(其中迭代法中的初值为算数平均数,即
)可得在考虑参数模糊随机性时,c, φ参数的均值与方差。2. 土钉加固边坡模糊随机极限状态方程
根据边坡稳定与失稳之间存在模糊过渡区间,在极限状态方程中引入边坡模糊过渡区间,建立极限状态方程。
令功能函数中的变量为模糊随机变量
,在功能函数中仅考虑土体参数c, φ为模糊随机变量(土体重度γ变异系数较小,可将其视为定值),则功能函数为(9) (10) 式中,
为考虑参数模糊随机性的抗滑力矩, 为考虑参数模糊随机性的下滑力矩。当加固边坡处于由失稳到稳定的过渡区间时,功能函数为
(11) 式中,
为模糊随机极限状态值。当 时,边坡为稳定状态;当 时,边坡处于失稳状态。令
为零点附近的有界闭模糊数,其隶属函数为关于零点对称的三角模糊数 (Z)。(12) 可转化为经典集,进行可靠度运算,则(13) 式中,
(14) 其中,bλ为
的λ-截集,d为 的最大容许值,λ为置信水平或阈值。3. 土钉加固边坡模糊随机可靠度分析及计算方法
依据模糊集分解定理:
(15) (16) (17) 则模糊随机极限状态方程在同一λ-截集下表示为
(18) 同一λ-截集下各项变量可用最大最小值表示为
(19) (20) (21) 式中,
(22) (23) 由于土钉与周围土体之间的作用主要体现在土钉锚固体与周围土体之间的侧摩阻力,其在理论上主要表现为土体抗剪强度τ,可将τ值作为土钉与周围土体之间的极限侧摩阻力[9],故土钉与周围土体的极限侧摩阻力(极限抗拉力)为
(24) 式中,γs为土钉抗拉承载力分项系数,取1.3;lni为第j层土钉在滑裂面外穿越第i层稳定土体内的长度;∑lni总和为为lnj, lnj为滑裂面以外土钉锚固体总长度;τik为第i层土体的抗剪强度标准值。
对于一般灌浆土钉在灌浆过程中未加特殊压力,其孔壁周边的法向压力σ主要取决于地层压力,故土层抗剪强度为
(25) 式中,h为稳定土层内锚固体以上的地层覆盖厚度(m), γ为锚固体土层的重度(kN/m3), K0为锚固体孔壁的土压力系数,一般取为K0=0.5~1.0。
根据式(24), (25),采用Morgenstern-Price法可得抗滑力与下滑力:
(26) (27) 式中,n为滑动体分条数,m为土钉层数,γ0为结构重要性系数,wi为第i条带的土重,bi为第i条带的宽度,cik为第i条带滑裂面处的黏聚力标准值,φik为第i条带滑裂面处的内摩擦角标准值,θi为第i条带滑裂面处的切线与水平面夹角,αj为土钉与水平面间的夹角,Li为第i条带滑裂面处的弧长,s为滑动体单元的厚度[10]。
假定式(26), (27)中的模糊随机变量仅为c, φ,且不考虑土体参数之间的相关性,土体重度γ为定值;同时,模糊随机变量分布类型为正态模糊分布,λ-截集为随机区间。
确定可靠度指标、失效概率、可靠度以及稳定性系数的上下限值,
(28) (29) (30) (31) 式中,
(32) (33) (34) (35) 最后,加权得到土钉加固边坡模糊随机可靠度。
4. 算例分析
4.1 工程概况
甘肃兰州地区某公路边坡支护,边坡高度为7 m,边坡重要性系数为1.0,边坡倾角为80°,边坡稳定性系数为1.3,采用土钉支护型式,支护剖面见图1。
4.2 土体参数对比
采用样本力学参数值的模糊-随机处理方法将土体参数的随机值转化为模糊随机变量,土体参数均值与方差对比见表1。
表 1 土体的物理力学参数Table 1. Physico-mechanical parameters of soils参数 c/kPa φ/(°) 考虑随机性 均值 16.5633 24.0411 方差 7.3752 3.8890 考虑模糊随机性 均值 16.2960 24.0278 方差 7.2643 3.7144 土体重度γ值为定值,其值为16.5 kN/m3,内摩擦角φ(最大值为33.3°,最小值为16.8°)与黏聚力c(最大值为24.9 kPa,最小值为11.4 kPa)为模糊随机变量,令
服从正态模糊分布,可求得同一λ-截集下的c, φ上下限值。4.3 应用本文方法计算模糊随机可靠度
根据正态模糊分布求得不同λ-截集的
与 值,采用Morgenstern-Price法确定未加固边坡模型的最危险滑移面,得到式(26), (27)中除变量c, φ外的其他参数值,并结合MATLAB软件计算不同λ-截集的R, S, FS, β, Pr上下限值,加权得到该加固边坡在考虑边坡模糊性后的可靠度[3],其中d=220 kN, λ取0.75, 0.80, 0.85, 0.90, 0.95, 1.0,计算结果见表2。表 2 模糊随机可靠度Table 2. Fuzzy random reliabilitiesλ c- c+ φ- φ+ b- b+ R- R+ S- S+ λ Z- Z+ μZ σZ β- β+ Prλ 0.75 14.8504 17.7416 22.9941 25.0615 -55 55 446.3929 486.7471 300.8916 327.4468 0.75 1.3633 1.6177 118.9461 185.8555 152.4008 62.3736 1.5616 3.3251 0.9408 0.9996 0.9856 0.80 15.0228 17.5692 23.1174 24.9382 -44 44 453.7152 528.0391 329.9547 337.3319 0.80 1.3450 1.6003 116.3833 198.0844 157.2338 86.478 1.3094 2.3270 0.9048 0.9900 0.9595 0.85 15.2095 17.3825 23.2508 24.8048 -33 33 425.6826 520.0764 300.8915 334.8997 0.85 1.2711 1.7285 90.7829 219.1849 154.9839 159.2538 0.7660 1.1804 0.7782 0.8811 0.8330 0.90 15.4211 17.1709 23.4022 24.6534 -22 22 463.7094 466.4952 295.2597 329.9547 0.90 1.4054 1.5799 133.7547 171.2355 152.4951 57.73513 2.2602 3.0223 0.9881 0.9987 0.9951 0.95 15.6856 16.9064 23.5913 24.4643 -11 11 470.4071 510.8604 329.9547 337.3319 0.95 1.3945 1.5483 133.0752 180.9057 156.9905 105.5953 1.3825 1.5909 0.9166 0.9442 0.9311 1.00 16.2960 16.2960 24.0278 24.0278 0 0 448.5147 448.5147 298.0881 298.0881 1.0 1.5046 1.5046 150.4266 150.4266 152.4266 94.2872 1.5954 1.5954 0.9447 0.9447 0.9447 4.4 可靠度结果对比分析
不同λ-截集可靠度加权求和得到Pr=0.9410(94.10%),利用GeoStudio软件进行未考虑边坡模糊性的蒙特卡洛模拟可靠度分析,分析类型为Morgenstern-Price, Monte Carlo Simulation为10000次,失效概率近似为0,可靠度为Pr=1.000000。(安全系数频率见图2,蒙特卡洛模拟分析结果见表3)
表 3 蒙特卡洛模拟分析结果Table 3. Analysis results of Monte Carlo simulation平均安全系数 可靠度指标 失效概率/% 最小安全系数 最大安全系数 试验次数 1.9819 5.9791 0.0000 1.4407 2.6482 10000 由于边坡自身模糊随机性的存在,采用本文方法所计算出的可靠度值低于蒙特卡洛模拟的可靠度值(蒙特卡洛模拟并未考虑边坡土体参数不确定性),两者相差0.059。蒙特卡洛模拟出该边坡处于稳定状态,而本文计算出的该加固边坡失效概率为5.9%,说明该加固边坡实际上存在失稳的可能性。
5. 结论
本文在考虑土体参数模糊随机性与边坡模糊过渡空间基础上,提出土钉加固边坡可靠度分析方法,并推导出土钉加固边坡可靠度计算公式。
(1)在土钉加固边坡工程中,应用蒙特卡洛模拟分析出加固边坡处于稳定状态,但在考虑边坡模糊性与模糊过渡区间后,加固边坡可能会出现失效概率。
(2)边坡自身存在的模糊性与模糊过渡区间对土钉加固边坡可靠度具有较大影响,在加固边坡中应考虑边坡模糊性与模糊过渡区间对边坡稳定性的影响。
(3)本文仅分析了土体侧摩阻力对土钉加固边坡的影响,然而在实际工程中,永久边坡会遇到地震作用、降雨入渗等外部作用,在后续研究应将各类外部作用考虑到土钉加固边坡可靠度分析中。
-
表 1 土体的物理力学参数
Table 1 Physico-mechanical parameters of soils
参数 c/kPa φ/(°) 考虑随机性 均值 16.5633 24.0411 方差 7.3752 3.8890 考虑模糊随机性 均值 16.2960 24.0278 方差 7.2643 3.7144 表 2 模糊随机可靠度
Table 2 Fuzzy random reliabilities
λ c- c+ φ- φ+ b- b+ R- R+ S- S+ λ Z- Z+ μZ σZ β- β+ Prλ 0.75 14.8504 17.7416 22.9941 25.0615 -55 55 446.3929 486.7471 300.8916 327.4468 0.75 1.3633 1.6177 118.9461 185.8555 152.4008 62.3736 1.5616 3.3251 0.9408 0.9996 0.9856 0.80 15.0228 17.5692 23.1174 24.9382 -44 44 453.7152 528.0391 329.9547 337.3319 0.80 1.3450 1.6003 116.3833 198.0844 157.2338 86.478 1.3094 2.3270 0.9048 0.9900 0.9595 0.85 15.2095 17.3825 23.2508 24.8048 -33 33 425.6826 520.0764 300.8915 334.8997 0.85 1.2711 1.7285 90.7829 219.1849 154.9839 159.2538 0.7660 1.1804 0.7782 0.8811 0.8330 0.90 15.4211 17.1709 23.4022 24.6534 -22 22 463.7094 466.4952 295.2597 329.9547 0.90 1.4054 1.5799 133.7547 171.2355 152.4951 57.73513 2.2602 3.0223 0.9881 0.9987 0.9951 0.95 15.6856 16.9064 23.5913 24.4643 -11 11 470.4071 510.8604 329.9547 337.3319 0.95 1.3945 1.5483 133.0752 180.9057 156.9905 105.5953 1.3825 1.5909 0.9166 0.9442 0.9311 1.00 16.2960 16.2960 24.0278 24.0278 0 0 448.5147 448.5147 298.0881 298.0881 1.0 1.5046 1.5046 150.4266 150.4266 152.4266 94.2872 1.5954 1.5954 0.9447 0.9447 0.9447 表 3 蒙特卡洛模拟分析结果
Table 3 Analysis results of Monte Carlo simulation
平均安全系数 可靠度指标 失效概率/% 最小安全系数 最大安全系数 试验次数 1.9819 5.9791 0.0000 1.4407 2.6482 10000 -
[1] 熊文林, 李胡生. 岩石样本力学参数值的随机-模糊处理方法[J]. 岩土工程学报, 1992, 14(6): 101-108. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC199206012.htm XIONG Wen-lin, LI Hu-Sheng. A random-fuzzy method for treating the experimental data of mechanical parameters of rock sample[J]. Chinese Journal of Geotechnical Engineering, 1992, 14(6): 101-108. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC199206012.htm
[2] 贾厚华, 贺怀建. 边坡稳定模糊随机可靠度分析[J]. 岩土力学, 2003, 24(4): 657-660. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200304042.htm JIA Hou-hua, HE Huaijian. Analysis of fuzzy-random reliability of slope stability[J]. Rock and Soil Mechanics, 2003, 24(4): 657-660. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200304042.htm
[3] 吕玺琳, 钱建固, 吕龙, 等. 边坡模糊随机可靠性分析[J]. 岩土力学, 2008, 29(12): 3437-3442. doi: 10.3969/j.issn.1000-7598.2008.12.048 LÜ Xi-lin, QIAN Jian-guo, LÜ Long, et al. Fuzzy stochastic reliability analysis of slopes[J]. Rock and Soil Mechanics, 2008, 29(12): 3437-3442. (in Chinese) doi: 10.3969/j.issn.1000-7598.2008.12.048
[4] 王宇, 贾志刚, 李晓, 等. 边坡模糊随机可靠性分析的模糊点估计法[J]. 岩土力学, 2012, 33(6): 1795-1800. doi: 10.3969/j.issn.1000-7598.2012.06.030 WANG Yu, JIA Zhi-gang, LI Xiao, et al. Fuzzy random reliability analysis of slope based on fuzzy point estimate method[J]. Rock and Soil Mechanics, 2012, 33(6): 1795-1800. (in Chinese) doi: 10.3969/j.issn.1000-7598.2012.06.030
[5] El-RAMLY H, MORGENSTERN N R, CRUDEN D M. Probabilistic slopestabilityanalysisforpractice[J]. Canadian Geotechnical Journal. 2002, 39(3): 665-683. doi: 10.1139/t02-034
[6] WONG F S. Slope reliability and response surface method[J]. Journal of Geotechcal Engineering, 1985, 111(1): 32-53. doi: 10.1061/(ASCE)0733-9410(1985)111:1(32)
[7] GRIFFITHS D. V, FENTON G A. Probabilistic slope stability analysisbyfiniteelements[J]. Journalof Geotechnical and Geoenvironmental Engineering, 2004, 130(5): 507-518. doi: 10.1061/(ASCE)1090-0241(2004)130:5(507)
[8] YE S. H, FANG G W, MA X R. Reliability analysis of grillage flexible slope supportingstructurewithanchors considering fuzzy transitional interval and fuzzy randomness of soil parameters[J]. Arabian Journal for Science and Engineering, 2019, 44(10): 8849-8857.
[9] 朱彦鹏, 罗晓辉, 周勇. 支挡结构设计[M]. 北京: 高等教育出版社, 2008. ZHU Yan-peng, LUO Xiao-hui, ZHOU Yong. Design of Retaining Structures[M]. Beijing: Higher Eduation Press. (in Chinese)
[10] 朱彦鹏, 李忠. 深基坑土钉支护稳定性分析方法的改进及软件开发[J]. 岩土工程学报, 2005, 27(8): 939-943. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200508018.htm ZHU Yan-peng, LI Zhong. Improvement on stability analysis of soil nailing in foundation excavations and its software development[J]. Chinese Journal of Geotechnical Engineering, 2005, 27(8): 939-943. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200508018.htm
-
期刊类型引用(4)
1. 韩非,胡超. 基于滑移线场理论的概率边坡设计. 科学技术与工程. 2023(20): 8771-8779 . 百度学术
2. 张继旭,王林峰,黄晓明,谭国金. 双向地震动的随机性对倾倒式危岩模糊可靠度的影响. 工程科学与技术. 2023(06): 161-171 . 百度学术
3. 张继旭,王林峰,夏万春. 基于模糊失效准则的危岩稳定可靠度计算. 防灾减灾工程学报. 2022(04): 695-704 . 百度学术
4. 陈昌富,李伟,张嘉睿,廖佳卉,吕晓玺. 山区公路边坡工程智能分析与设计研究进展. 湖南大学学报(自然科学版). 2022(07): 15-31 . 百度学术
其他类型引用(4)