边坡支挡结构岩土荷载的分项系数计算方法

2022-06-22 08:43王桂林陈瑶杨洋
土木建筑与环境工程 2022年5期
关键词:可靠度

王桂林 陈瑶 杨洋

摘 要:巖土参数的变异特性给边坡支挡结构岩土荷载的确定带来困难,分项系数法考虑了可靠度,且计算简单,为此提供了解决方案。基于极限状态方程,将抗剪强度参数视为正态分布,采用一般分离法,通过分项系数的求解,推导了边坡支挡结构岩土荷载的分项系数以及荷载值计算公式;对单一滑面和折线型滑面两类典型滑面边坡案例分析表明:利用该方法能较容易地获得不同边坡目标可靠度指标值对应的支挡结构岩土荷载值;根据分项系数法计算得到的荷载值所设计的边坡,其由蒙特卡罗法获得的可靠度指标与目标可靠度指标存在一定偏差;其中,单一滑面边坡主要是负偏差,偏差最大约6%,对于折线型滑面边坡,当目标可靠度指标为1.65时,是负偏差,偏差约7%,但当目标可靠度指标大于1.65时出现正偏差。由分项系数法计算得到的支挡结构荷载值总体上安全可控。

关键词:支挡结构;岩土荷载;分项系数法;可靠度;蒙特卡罗法

中图分类号:TU432 文献标志码:A 文章编号:2096-6717(2022)05-0001-07

收稿日期:2020-08-27

基金项目:国家重点研发计划(2018YFC1505504)

作者简介:王桂林(1970- ),男,博士,教授,主要从事岩土工程研究,E-mail:glw@cqu.edu.cn。

Received:2020-08-27

Foundation item:National Key R&D Program of China (No. 2018YFC1505504)

Author brief:WANG Guilin (1970- ),PhD, professor, main research interest: geotechnical engineering, E-mail: glw@cqu.edu.cn.

Partial coefficient determination method for loadings of slope retaining structure

WANG Guilin, CHEN Yao, YANG Yang

(a. School of Civil Engineering; b. National Joint Engineering Research Center of Geohazards Prevention in the Reservoir Areas; c. Key Laboratory of New Technology for Construction of Cities in Mountain Area, Ministry of Education, Chongqing University, Chongqing 400045, P.R. China)

Abstract:The geotechnical parameters of slope have great variability which makes it difficult to determine the load of the support structure accurately. The partial coefficient method,which considers design uncertainty and is simple to calculate, is capable of solving this problem effectively. Based on the limit state equation, this study takes the shear strength parameter as normal distribution,derived the calculation formulas of the partial coefficients and the loading of the slope support structure via the general separation method.Then it analyzed two typical slope cases: single sliding surface and polyline sliding surface, by using this proposed method. The results show that the loading values corresponding to the target reliability index of different slopes can be easily obtained;and there is a deviation between the reliability index obtained by the Monte Carlo method and the target reliability value by the presented method; for single slide slope, the deviation is mainly negative about 6%, while for polyline sliding type, when the target reliability index is 1.65, the negative deviation is about 7%, but when the target reliability index is greater than 1.65, the deviation became positive. Therefore, the load value of the retaining structure calculated by the partial coefficient method is generally unconservative.

Keywords:retaining structure; rock and soil load; partial coefficient method; reliability; Monte Carlo method

第5期   王桂林,等:边坡支挡结构岩土荷载的分项系数计算方法

目前,边坡(滑坡)支挡结构岩土荷载常采用单一安全系数法来确定,但安全系数受工程经验影响较大,不能充分反映支挡结构设计的岩土参数变异性与安全裕度。因此,学者们对边坡支挡结构岩土荷载计算进行了深入研究。

方玉树分析了支挡结构岩土荷载取值方法,认为应根据稳定性进行计算,并指出支挡结构岩土荷载既不是滑坡推力也不是主动土压力。付文光等指出相关标准中稳定安全系数计算公式存在着5类抗力与荷载错位、不符合安全系数定义的现象,导致有时工程安全储备不足或较大浪费。Wu等首次提出利用可靠度理论评价边坡稳定性,此后,利用可靠度理论求解边坡支挡结构岩土荷载的方法逐渐发展。祝玉学对可靠度用于边坡工程做了专门的解释;唐小松等基于Copula理论研究了有限信息情况下的边坡可靠度研究方法;周泽华系统研究了抗剪强度参数的联合分布与边坡可靠度的关系。但可靠度方法计算繁琐复杂,未能在实际工程中得到推广使用。分项系数极限状态设计方法是一种与规定目标可靠度相联系的设计方法,通过各变量的分项系数反映其变异水平,使结构功能函数通过显示方式表达,提高了准确度,广泛应用于桥梁、水利、灾害防护等多个结构工程领域的安全性研究。在边坡工程方面,陈祖煜根据安全系数与可靠指标提出了相对安全率的概念,并研究其作为安全判据在重力坝边坡等工程中的适用性。黎康平通过对比基于可靠度的安全系数标准和分项系数标准,研究了重力墙及土石坝的稳定性。李昂等研究了几何和强度参数对分项系数的影响规律,提出了对应小、中、大3种变异水平的分项系数推荐值。但现有边坡工程设计对分项系数法的应用与规定还较为简单,有待进一步研究。

为了进一步探索分项系数法在边坡支挡结构设计方面的适用性,笔者采用一般分离法推导出边坡的分项系数以及荷载值的计算公式,将抗剪强度参数黏聚力c和摩擦系数f视为随机变量,通过算例比较分项系数法与蒙特卡罗法在单一型滑面和折线型滑面边坡可靠度的偏差,并分析其原因。

1 分项系数法

分项系数γ通常为变量的设计值与标准值的比值。一般情况下,影响边坡稳定性的因素分为边坡的抗滑力R与边坡的下滑力S,描述边坡稳定性的函数定义为其功能函数Z,表示为Z=g(R,S)=R-S。现将支挡结构岩土荷载R加入到功能函数,得到边坡的功能函数在极限状态下的表达式,见式(1),示意图见图1。

Z=R+R-S=0(1)

式中:R为支挡力,kN,其方向与所在条块滑面倾向相同;S为边坡下滑力;R为边坡抗滑力,计算公式见式(2)、式(3)。

S=∑ni=1Wsin α(2)

R=∑ni=1lc+∑ni=1(Wcos α)f(3)

式中:∑ni=1l为折线型滑面总长度,m;c为黏聚力,kPa;W为各土条的重力,kN;α为各土条滑面法线与铅垂线的夹角,(°);f为摩擦系数,f=tan φ,φ为内摩擦角,(°)。

图1中,P为第i滑块的剩余下滑力,为该滑块的下滑力与抗滑力差值;R为支挡结构对最底部滑块的总支挡力;R和T分别为沿着底部滑块倾向和垂直底部滑块倾向的支挡分量。

根据一般分离法思想,可对任意变量X、X进行式(4)所示的公式变换。

X+X=X+XX+X=ΦX+ΦX(4)

式中:Φ、Φ为分離函数。对于变量X(i=1,2…n),其分离函数为Φ=X(∑ni=1X)。

对抗剪参数c、φ的分项系数γ和γ进行求解。假定下滑力S和抗滑力R均服从正态分布,根据一般分离法及可靠度指标的表达式,可得式(5)。

μ-μ+R=β(σ+σ)=

βΦσ+βΦσ(5)

式中:μ为抗滑力均值;μ为下滑力均值;β为可靠度;σ为抗滑力标准差;σ为下滑力标准差;Φ、Φ为分离函数,表达式分别见式(6)、式(7)。

Φ=σσ+σ(6)

Φ=σσ+σ(7)

将式(2)和式(3)代入式(5),可得式(8)。

∑ni=1lμ+∑ni=1Wcos αμ+R-∑ni=1Wsin α=

β((∑ni=1l)σ+(∑ni=1Wcos α)σ)=

βΦσ∑ni=1l+βΦσ(∑ni=1Wcos α)(8)

根据变异系数的定义,抗剪参数均值与标准差存在关系

σ=δμ,σ=δμ

式中:δ、δ分别为黏聚力和摩擦系数的变异系数,将其代入式(8),可得式(9)。

∑ni=1lμ+(∑ni=1Wcos α)μ-∑ni=1Wsin α+R=

βΦδμ∑ni=1l+βΦδμ(∑ni=1Wcos α) (9)

移项整理后得到边坡支挡结构岩土荷载的计算式,见式(10)。

R=∑ni=1Wsin α-∑ni=1lγμ-

(∑ni=1Wcos α)γμ(10)

式中:γ为黏聚力的分项系数,γ=1-βΦδ,当γ<0时,取γ=0;γ为摩擦系数的分项系数,γ=1-βΦδ,当γ<0时,取γ=0;Φ、Φ分别为黏聚力和摩擦系数的分离函数,表达式见式(11)、式(12)。

Φ=∑ni=1lσ(∑ni=1l)σ+(∑ni=1Wcos α)σ(11)

Φ=(∑ni=1Wcos α)σ(∑ni=1l)σ+(∑ni=1Wcos α)σ(12)

2 与目标可靠度指标印证及讨论

由于分项系数可以反映支擋设计后边坡需要达到的可靠度指标,因此,按分项系数法计算所得到的荷载进行边坡支护,支护后边坡的可靠度指标理论上应等于计算时设定的目标可靠度指标,即分项系数设定的目标可靠度与实际可靠度之间存在相互印证,印证的流程如图2所示。

印证时,边坡实际可靠度采用蒙特卡罗法计算,通过输入黏聚力与摩擦系数的分布区间以及分项系数法得到的支挡荷载,计算判断边坡是否失稳破坏,取重复次数N=5 000,得到失稳破坏次数与N的比值即为失效概率,进而得到可靠度β。

选取单一滑面和折线型滑面两类典型滑面的边坡案例进行计算印证。根据概率统计知识,当功能函数服从正态分布时,失效概率P与可靠度指标β的关系见表1。

2.1 算例一

取单一滑面的边坡模型,上部为土层,下部为岩体,边坡几何参数有上覆土层厚度H,滑面倾角α,如图3所示。模型参数取值分别为H=5 m、α=45°,土体重度γ=18 kN/m。抗剪强度参数黏聚力c和摩擦系数f为随机变量,其统计特征见表2。

2.1.1 边坡支挡结构岩土荷载计算

分离函数只与抗剪强度参数标准差有关,与目标可靠度指标无关,根据式(11)和式(12)计算分离函数。

黏聚力的分离函数

Φ=∑ni=1l×σ(∑ni=1l)σ+(∑ni=1Wcos θ)σ=0.646 2

摩擦系数的分离函数

Φ=(∑ni=1Wcos θ)×σ(∑ni=1l)σ+(∑ni=1Wcos θ)σ=0.763 2

按照式(10)进行支挡结构岩土荷载的计算。不同目标可靠度指标情况下,分项系数及支挡结构岩土荷载见表3。

2.1.2 支挡后可靠度与目标可靠度对比

根据表3的边坡支挡结构岩土荷载对边坡进行支挡设计,应用蒙特卡罗法求得支挡后边坡的可靠度指标β,并与分项系数法目标可靠度β进行比较,结果见表4。

表4给出了目标可靠度指标β和对应的支挡后边坡可靠度指标β以及二者之间的相对误差。支挡后边坡可靠度指标β与目标可靠度指标β相对误差绝对值在0.9%~5.8%之间,其原因有:1)可靠度指标的计算方法不同,分项系数法计算过程是线性运算,而蒙特卡罗法的可靠度是由概率计算获得;2)抗剪强度参数变异性处理方式不同,分项系数法只考虑了抗剪强度参数概率上的变异性,而蒙特卡罗法采用随机场理论,考虑了抗剪强度参数的空间变异性。

2.2 算例二

某实际折线型滑面边坡上部主要为次生红黏土,少量碎石土和素填土,下面为岩体,滑动面为土岩分界面。边坡剖面高约420 m,长约120 m,土体重度γ=18 kN/m,黏聚力c和内摩擦角φ视为随机变量,统计特征见表5,边坡的计算剖面见图4,本算例不考虑地下水作用。

根据式(11)和式(12),得到黏聚力的分离函数Φ=0.306 2;摩擦系数的分离函数Φ=0.952 0。

2.2.1 边坡支挡结构岩土荷载计算

分别取目标可靠度指标进行支挡结构岩土荷载的计算,在不同目标可靠度指标情况下,其分项系数及支挡结构岩土荷载见表6。

由表6可以看出,随着目标可靠度指标β取值的增加,支挡结构岩土荷载大幅度增加,需要支护。目标可靠度指标β的选取对支挡结构的设计有重要影响。

2.2.2 支挡后可靠度与目标可靠度对比及讨论

根据表6得到的支挡结构岩土荷载对边坡进行支挡设计,应用蒙特卡罗法求得支挡后边坡的可靠度指标β,并与分项系数法目标可靠度β进行比较,结果见表7。

表7给出了目标可靠度指标β和对应的支挡后边坡可靠度指标β,以及二者之间的相对误差。可以看出,相比算例1,算例2中目标可靠度β和支挡后边坡可靠度指标β相差较大,当目标可靠度β为1.65时,存在负偏差-7.39%,当目标可靠度β大于1.65时则为正偏差,可见总體上是安全可控的。

β和β产生误差的原因:除可靠度指标计算方法和抗剪强度参数变异性处理方式不同外,与计算过程结构功能函数是显式表达式Z=R+R-S且忽略了条块间力传递的假定有关;除此之外,假定抗剪强度参数服从正态分布,未考虑各参数之间的关联性,然而,实际工程中岩土体参数分布类型的拟合并不都是正态分布,且各参数之间存在一定的关联性。

后续研究可从以上几个方面改进分项系数计算方法,并通过大量实例对相互印证过程进行对比分析,通过对比结果确定其相互印证程度,并采取有效措施改进分项系数计算方法和分项系数极限状态设计方法。

3 结论

采用一般分离法推导了边坡支挡结构岩土荷载分项系数表达式,并经算例进行了目标可靠度与实际可靠度的对比。结果表明:

1)通过分项系数法能方便地获得不同边坡目标可靠度指标值对应的支挡结构岩土荷载值。

2)在指定目标可靠度指标值时,单一滑面边坡的可靠度主要出现负偏差,偏差绝对值最大约6%;对于折线型滑面边坡,当目标可靠度指标为1.65时,是负偏差,偏差约7%,但目标可靠度指标大于1.65时,出现正偏差,总体上是安全可控的。

3)导致偏差的主要原因有:分项系数法与可靠度法可靠度指标的具体计算方法不同;抗剪强度参数变异性处理方式不同;结构功能函数忽略了条块间力;岩土体参数分布类型仅假定为正态分布且参数之间无关联。

4)单一滑面和折线型滑面边坡可靠度分析中,目标可靠度与支挡后边坡可靠度相对大小并不是固定的,即便如此,由于分项系数法可以考虑工程的安全裕度,计算简单,符合工程设计习惯,在工程中具有应用价值。

参考文献:

[1] 建筑边坡工程技术规范: GB 50330—2013 [S]. 北京: 中国建筑工业出版社, 2013.

Technical code for building slope engineering: GB 50330-2013 [S]. Beijing: China Architecture & Building Press, 2013. (in Chinese)

[2] 滑坡防治设计规范: GB/T 38509—2020 [S]. 北京: 中国标准出版社, 2020.

Code for the design of landslide stabilization: GB/T 38509—2020 [S]. Beijing: Standards Press of China, 2020. (in Chinese)

[3] 赵尚毅, 郑颖人, 王建华, 等. 基于强度折减安全系数的边坡岩土侧压力计算方法探讨[J]. 岩石力学与工程学报, 2010, 29(9): 1760-1766.

ZHAO S Y, ZHENG Y R, WANG J H, et al. Discussion on calculation method of slope rock-soil lateral pressure based on factor of safety with strength reduction [J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(9): 1760-1766. (in Chinese)

[4] 赵尚毅, 郑颖人, 何平, 等. 边坡侧向荷载计算方法讨论[J]. 地下空间与工程学报, 2017, 13(2): 434-441.

ZHAO S Y, ZHENG Y R, HE P, et al. Discussion on calculation methods of slope lateral load [J]. Chinese Journal of Underground Space and Engineering, 2017, 13(2): 434-441. (in Chinese)

[5] 李广信. 关于《建筑边坡工程技术规范GB 50330—2013》的讨论[J]. 岩土工程学报, 2016, 38(12): 2322-2326.

LI G X. Discussions on “Technical code for building slope engineering (GB 50330-2013)” [J]. Chinese Journal of Geotechnical Engineering, 2016, 38(12): 2322-2326. (in Chinese)

[6] 方玉树. 边坡与滑坡抗滑支护结构岩土荷载既不是主动土压力也不是剩余下滑力[J]. 重庆建筑, 2016, 15(9): 27-33.

FANG Y S. Rock and soil load of anti-sliding support structure of side slope and landslide is neither active earth pressure nor residual slide force [J]. Chongqing Architecture, 2016, 15(9): 27-33. (in Chinese)

[7] 付文光, 陈双, 王振威. 稳定安全系数计算公式中荷载与抗力错位影响探讨[J]. 岩土工程学报, 2021, 43(3): 556-563.

FU W G, CHEN S, WANG Z W,Discussion of dislocation phenomena of resistance and load in formula for stability safety factor [J]. Chinese Journal of Geotechnical Engineering, 2021, 43(3): 556-563. (in Chinese)

[8] WU T H, KRAFT L M Jr. Safety analysis of slopes [J]. Journal of the Soil Mechanics and Foundations Division, 1970, 96(2): 609-630.

[9] 祝玉學. 边坡可靠性分析[M]. 北京: 冶金工业出版社, 1993.

ZHU Y X.Slope reliability analysis [M]. Beijing: Metallurgical Industry Press, 1993. (in Chinese)

[10] 唐小松, 李典庆, 周创兵, 等. 基于Copula函数的抗剪强度参数间相关性模拟及边坡可靠度分析[J]. 岩土工程学报, 2012, 34(12): 2284-2291.

TANG X S, LI D Q, ZHOU C B, et al. Modeling dependence between shear strength parameters using Copulas and its effect on slope reliability [J]. Chinese Journal of Geotechnical Engineering, 2012, 34(12): 2284-2291. (in Chinese)

[11] 周泽华. 岩土参数的联合分布识别与边坡稳定性动力可靠度分析[D]. 西安: 西安理工大学, 2017.

ZHOU Z H.Identification of the joint distribution of rock and soil parameters and analysis for slope dynamic reliability. [D]. Xi'an: Xi'an University of Technology, 2017. (in Chinese)

[12] 李夢, 高晓军, 王启鸿. 基于分项系数法对拱坝抗滑稳定的分析[J]. 水利与建筑工程学报, 2010, 8(3): 136-137, 141.

LI M, GAO X J, WANG Q H. Analysis on stability against sliding of arch dam based on partial coefficient method [J]. Journal of Water Resources and Architectural Engineering, 2010, 8(3): 136-137, 141. (in Chinese)

[13] DE BIAGI V, MARCHELLI M, PEILA D. Reliability analysis and partial safety factors approach for rockfall protection structures [J]. Engineering Structures, 2020, 213: 110553.

[14] BLOMFORS M, LARSSON IVANOV O, HONF D, et al. Partial safety factors for the anchorage capacity of corroded reinforcement bars in concrete [J]. Engineering Structures, 2019, 181: 579-588.

[15] 陈祖煜. 建立在相对安全率准则基础上的岩土工程可靠度分析与安全判据[J]. 岩石力学与工程学报, 2018, 37(3): 521-544.

CHEN Z Y.Reliability analysis and safety criterion in geotechnical engineering based on the index of safety margin [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(3): 521-544. (in Chinese)

[16] 黎康平. 基于可靠度理论的土工构筑物抗滑稳定安全标准研究[D]. 西安: 西安理工大学, 2018.

LI K P.Research on safety standards in earth works stability against sliding based on reliability theory [D]. Xi'an: Xi'an University of Technology, 2018. (in Chinese)

[17] 李昂, 罗强, 张文生, 等. 土坡稳定分析中考虑参数变异水平的分项系数取值[J]. 铁道科学与工程学报, 2019, 16(2): 341-350.

LI A, LUO Q, ZHANG W S, et al. Discussion on the value of partial factors considering variation level of geotechnical parameter in slope stability analysis [J]. Journal of Railway Science and Engineering, 2019, 16(2): 341-350. (in Chinese)

(编辑 胡玲)

猜你喜欢
可靠度
具有旋转对称结构零部件的可靠度计算方法研究
基于BIM技术的建筑日照与风环境研究
非线性Mohr—Coulomb破坏准则下边坡可靠度上限
覆冰荷载作用下高压输电塔可靠性研究
边坡稳定可靠度分析的新型四阶矩法
公路水泥混凝土路面设计问题解决策略
Weibull分布下复杂系统可靠度的Bayes估计
考虑失效准则和沉降控制条件的基桩可靠度分析
水泥混凝土疲劳寿命可靠度研究