基于极限平衡法的非均质边坡稳定性分析.docx
《基于极限平衡法的非均质边坡稳定性分析.docx》由会员分享，可在线阅读，更多相关《基于极限平衡法的非均质边坡稳定性分析.docx（10页珍藏版）》请在冰豆网上搜索。
基于极限平衡法的非均质边坡稳定性分析
基于极限平衡法的非均质边坡稳定性分析
摘要：
概述了边坡稳定性分析的4种极限平衡法及适用范围，建立计算模型，采用GEO-SLOPE/W程序对某非均质、人工高填土边坡进行稳定性分析。
结果表明：
非均质、人工高填土边坡稳定性分析不适宜采用瑞典条分法；当滑面为圆弧型时，简布法、摩根斯坦-普莱斯法、简化毕肖普法安全系数计算值偏差极小；当滑面为折线形时，简化毕肖普法安全系数计算值偏差较大，简布法、摩根斯坦-普莱斯法偏差较小。
滑面为圆弧型或折线形，通用简布法、摩根斯坦-普莱斯法各自偏差均较小。
上部土层为孔隙率低、渗透率小的人工填土边坡，分析其非常运行工况时，宜按岩体从上而下逐渐饱和至一定深度范围、且该范围内的土体按饱和状态进行强度折减，不宜按稳定渗流场考虑，否则偏差很大。
关键词：
极限平衡法、非均质、稳定性分析、人工高填土边坡、滑面、偏差
中图分类号：
TV16文献标志码：
A
Stabilityanalysisofananisotropicslopebasedonlimitequilibrium
HUYanpeng，LIUBin
（PowerChinaZhongnanEngineeringCorporationLimited,Changsha,Hunan,410014）
Abstract:
Summarized4limitequilibriumandtheapplicationrangeinstabilityanalysis,establishedacalculationmodelbasedonengineeringpractice,applyedGEO-SLOPE/Wsoftware,andanalysedthestabilityofaanisotropicandhighfilledslope.TheresultsindicatedthatOrdinaryisunfitforstabilityanalysisoftheanisotropicandhighfilledslope.Whentheslipsurfaceiscircular,securitycoefficientdeviationisslightestbyJanbu,Mogenstern-priceandBishop.Whentheslipsurfaceisfoldline,securitycoefficientdeviationbyBishopislargerbyJanbu,andsmallerbyMogenstern-price.Whentheslipsurfaceiscircularorfoldline,securitycoefficientdeviationissmallerbyJanbugeneralizedandMogenstern-price.Tohighfilledslopewhoseuppersoillayerwithlowerporosityandpermeability,inordertoanalysizetheunconventionalworkingcondition,itisadvisabletoconsiderthattherock-soilbodyisgraduallysaturatetoacertaindepthanditsstrengthisreductiveaccordingtothesaturationstate,anditisunadvisabletoconsiderasteadyseepagefield.Ifnot,thedeviationismuchlarger.
Keywords:
limitequilibrium；anisotropic；stabilityanalysis；highfilledslope；slipsurface；deviation
0概述
随着国家基础建设的不断发展，越来越多的人工高填土边坡在城郊、城市出现，该类边坡的稳定关系着人民群众生命和财产安全，对其稳定性分析显得尤为重要。
边坡种类及其稳定性分析方法众多，国内学者做了诸多研究：
胡海［1］等通过GeoStudio为该排土场在生产上提供理论上的技术基础，降低了安全事故发生的概率；刘建伟
［2］等
针对弃渣场边坡易发生细沟侵蚀问题，采集弃渣样品、测定物理指标，对弃渣场边坡稳定性特征进行研究，得到了弃渣场坡面上部土体较中、下部土体稳定性弱和弃渣细粒物质组成对边坡稳定性有重要影响的结论；陈勋辉［3］等采用三种极限平衡法对比分析了某土石坝下游坝坡滑动面各条块底部法向力、切向力分布，得出了毕肖普法、M-P法计算条间力、安全系数分布规律基本一致，瑞典条分法因忽略条间力导致安全系数偏小的结论；韩松江［4］应用强度折减法探讨了不同浸润线高度条件下边坡的安全系数、塑性区范围的变化趋势；李亮［5］基于尾矿库概化模型，假定三种平均坡度、四种浸润线深度，采用FLAC2D软件对尾矿坝的抗滑稳定性程度进行了模拟分析，得到了不同条件下尾矿坝的剪应变云图和相对应的抗滑稳定安全系数；郭方琴［6］等对非饱和土体进行了二维渗流和边坡稳定计算，分析了不同降雨强度、不同降雨持时对边坡稳定性的影响；贺炜［7］等对修建于蓄洪区范围内的公路边坡引入非饱和渗流理论，进行了浸润线及安全系数变化规律的研究，得到了浸润线稳定状态与蓄洪时段的关系，土体渗透系数、蓄洪时段与安全系数升降速率的关系；徐晗［8］等建立了非饱和土流固耦合有限元模型，对降雨条件下非饱和土边坡渗流场和应力场耦合进行了数值模拟，得到了非饱和土边坡变形与应力的若干规律。
前述研究对边坡稳定性的分析主要分为三类：
各类边坡正常运行工况下的稳定性；求解浸润线和稳定渗流场状态下坝体边坡的稳定性；运用非饱和渗流理论、求解暂态渗流场状态下山体、路堤等边坡的稳定性。
但是，目前针对水电工程中人工填土形成的弃渣场非常运用工况下边坡的稳定性研究较少。
本文利用GEO-SLPOE/W程序，得到了弃渣场正常运行工况下4种极限平衡法求解的安全系数，并分析了各计算方法的合理性。
同时，对弃渣场非常运行工况（降雨工况）下边坡稳定性进行了分析，在考虑土体孔隙水压力时进行了两种假定：
1.具有稳定的渗流场和浸润线；2.从上而下一定深度范围内视为饱和区、岩体物理力学参数按饱和土体进行折减；分析、对比计算结果，得到了非均质、人工填土边坡非常运行时稳定性分析宜按第二种假定考虑。
1计算方法介绍
1.1瑞典条分法（Ordinary）
瑞典条分法［9］是最早、最简单的粘性土坡分析方法，假定滑动面为圆弧形，其忽略条块间切向力和法向力，满足整体力矩平衡条件，不满足条间静力平衡条件，其计算结果一般偏于安全，是工程上常用的方法之一。
瑞典条分法的安全系数计算公式为：
（1）
式中：
ci为土条粘聚力；li为土条滑弧长度；Wi为土条自重；θi为土条底面倾斜角；φi为土条内摩擦角。
1.2简化毕肖普法（Bishop）
简化毕肖普法［10,11］仅考虑条块间法向力、不考虑条块间切向力，满足整体力矩平衡条件，满足各条块力的多边形闭合条件、不满足条块的力矩平衡条件，其计算结果较瑞典条分法精度高，适用于均质粘性土的圆弧形或近似圆弧形滑动滑坡，是工程上使用较多的一种方法。
其安全系数计算公式为：
（2）
（3）
式中：
ci为土条粘聚力；bi为土条宽度；Wi为土条自重；θi为土条底面倾斜角；φi为土条内摩擦角。
1.3通用简布法（Janbu）
通用简布法［12］同时考虑条块间法向力和切向力，且假定条块间法向力的作用位置，满足静力平衡条件、极限平衡条件、整体力矩平衡条件，计算准确但较复杂，主要适用于复合破坏面的边坡，既可用于圆弧滑动，也可用于非圆弧滑动。
其安全系数计算公式为：
（4）
（5）
式中：
ci为土条粘聚力；bi为土条宽度；Wi为土条自重；θi为土条底面倾斜角；φi为土条内摩擦角；Hi+1、Hi为相邻土条切向力，△H为相邻土条切向力差值，当△H=0时，则称简化简布法［13］。
1.4摩根斯坦-普莱斯法（Mogenstern-price，简称M-P法）
摩根斯坦-普莱斯法［14］考虑条块间法向力和切向力，其假定条块间切向力与法向力之间存在一个函数关系，该函数关系可用水平坐标x表达。
该方法满足静力平衡条件、整体力矩平衡条件，是极限平衡法中常见方法。
其安全系数计算公式为：
（6）
（7）
式中：
ci为土条粘聚力；bi为土条宽度；Wi为土条自重；θi为土条底面倾斜角；φ
i为土条内摩擦角；u为孔隙水压力；Ni为条块底部法向力；D为集中荷载；d为集中荷载力矩；β为折算系数；R为滑弧半径；f为摩擦系数；F（x）为切向力与法向力的函数；λ为函数所用的百分数。
2算例分析
2.1计算模型、工况及方法
（1）模型建立
西南地区某电站临河型弃渣场地，边坡总高度约80m、人工填土高度约45m，坡脚设置混凝土挡渣墙。
岩土层从上至下依次为填渣层、覆盖层、强风化层、弱风化层。
填渣主要为强风化砂岩、粉砂岩或泥质粉砂岩，遇水泥化；地表覆盖层主要为残坡积粉质粘土，局部夹碎块石，坡面零星分布有大小不等的崩石；下伏基岩地层为侏罗系蓬莱镇组灰绿色细砂岩、紫红色泥岩、粉砂岩，概化后计算模型如图1所示。
图1边坡计算模型
Fig.1Calculationmodelofslope
（2）计算工况及参数选取
本文分析3种工况。
工况1：
正常运用工况，各岩土层处于天然状态，取天然容重及强度指标。
工况2：
非常运用工况，假定边坡内部有稳定的渗流场，即假定浸润线高度，浸润线以上岩土层处于天然状态、取天然状态参数；浸润线以下取浮容重γsat和饱和强度指标。
工况3：
非常运用工况，假定浸润线高度，浸润线以上取饱和容重及强度指标，浸润线以下取天然容重及强度指标。
填渣体及覆盖层参数见表1所示。
需要说明的是，本文浸润线高度的假定基于边坡非常运行工况时（降雨工况），雨水下渗致使部分岩层饱和。
假定2个浸润线高度，分别位于：
a、填渣层-覆盖层接触面；b、覆盖层-强风化层接触面，见图2。
对于该2个浸润线高度，采取以下两种处理方式：
①按照水头边界条件给定边坡左、右侧水头，即视为边坡上下游浸水、按稳定渗流场考虑，接触面就假定为该渗流场的已知浸润线，这种方式即为工况2；②考虑降雨有一定强度、持续时间以及雨水渗流规律，将接触面上、下岩体分别作为饱和土体、天然状态土体考虑，这种方式即为工况3。
表1岩土体物理力学参数
Tab.1Physicalandmechanicalparametersofrock-soil
岩土体
容重γ
kN/m³
内摩擦角φ
°
凝聚力c
kPa
天然状态
饱和状态
天然状态
饱和状态
天然状态
饱和状态
覆盖层
17
19
20
18
25
12
填渣
17
19.5
23
19
15
10
图2假定浸润线a和b
Fig.2Hypothetical seepagelineaandb
（3）计算程序及方法
本文采用GEO-SLOPE/W程序，该程序包括了前述4种极限平衡法，便于计算结果的对比、分析。
2.2计算结果分析
2.2.1工况1稳定性分析
工况1土条划分及滑面位置见图3。
不
升级会员 