匀质粘性土体边坡稳定性计算

合集下载

边坡稳定性计算书 5#

边坡稳定性计算书 5#

边坡稳定性计算书计算依据:1、《建筑基坑支护技术规程》JGJ120-20122、《建筑边坡工程技术规范》GB50330-20023、《建筑施工计算手册》江正荣编著一、基本参数边坡稳定计算方式平面滑动法边坡工程安全等级三级边坡边坡土体类型粘性土土的重度γ(KN/m3) 20土的内摩擦角φ(°)15 土的粘聚力c(kPa) 12边坡高度H(m) 9 边坡斜面倾角α(°)40坡顶均布荷载q(kPa) 0.2二、边坡稳定性计算计算简图滑动体自重和顶部所受荷载:W= (1/2γH+q)×H×(ctgω-ctgα)=1/2(γH+2q)×H×sin(α-ω)/sinω/sinα边坡稳定性系数为:K s=(W×cosω×tanφ+H/sinω×c)/(W×sinω)=cotω×tanφ+2c/(γH+2q)×sinα/(sin(α-ω)×sinω)滑动面位置不同,Ks值亦随之而变,边坡稳定与否根据稳定性系数的最小值Ksmin 判断,相应的最危险滑动面的倾角为ω0。

求K smin值,根据dKs/dω=0,得最危险滑动面的倾角ω0的值:ctgω=ctgα+(a/(tanφ+a))0.5×cscα式中:a=2c/(γH+2q)= 2×12/(20×9+2×0.2)= 0.133ctgω=ctgα+(a/(tanφ+a))0.5×cscα= ctg(40°)+(0.133/(tan(15°)+0.133))0.5×csc(40°) = 2.088 则边坡稳定性最不利滑动面倾角为:ω0= 25.591°K smin=(2a+tanφ)×ctgα+2×(a(tanφ+a))0.5×cscα=(2×0.133+tan(15°))×ctg(40°)+2×(0.133×(tan(15°)+0.133) )0.5×csc(40°)=1.355≥1.25满足要求!。

(整理)边坡稳定性计算方法

(整理)边坡稳定性计算方法

一、边坡稳定性计算方法在边坡稳定计算方法中,通常采用整体的极限平衡方法来进行分析。

根据边坡不同破裂面形状而有不同的分析模式。

边坡失稳的破裂面形状按土质和成因不同而不同,粗粒土或砂性土的破裂面多呈直线形;细粒土或粘性土的破裂面多为圆弧形;滑坡的滑动面为不规则的折线或圆弧状。

这里将主要介绍边坡稳定性分析的基本原理以及在某些边界条件下边坡稳定的计算理论和方法。

(一)直线破裂面法化计算这类边坡稳定性分析采用直线破裂面法。

能形成直线破裂面的土类包括:均质砂性土坡;透水的砂、砾、碎石土;主要由内摩擦角控制强度的填土。

图 9 - 1 为一砂性边坡示意图,坡高 H ,坡角β,土的容重为γ,抗剪度指标为c、φ。

如果倾角α的平面AC面为土坡破坏时的滑动面,则可分析该滑动体的稳定性。

沿边坡长度方向截取一个单位长度作为平面问题分析。

图9-1 砂性边坡受力示意图已知滑体ABC重 W,滑面的倾角为α,显然,滑面 AC上由滑体的重量W= γ(ΔABC)产生的下滑力T和由土的抗剪强度产生的抗滑力Tˊ分别为:T=W · sina和则此时边坡的稳定程度或安全系数可用抗滑力与下滑力来表示,即为了保证土坡的稳定性,安全系数F s 值一般不小于 1.25 ,特殊情况下可允许减小到 1.15 。

对于C=0 的砂性土坡或是指边坡,其安全系数表达式则变为从上式可以看出,当α =β时,F s 值最小,说明边坡表面一层土最容易滑动,这时当 F s =1时,β=φ,表明边坡处于极限平衡状态。

此时β角称为休止角,也称安息角。

此外,山区顺层滑坡或坡积层沿着基岩面滑动现象一般也属于平面滑动类型。

这类滑坡滑动面的深度与长度之比往往很小。

当深长比小于 0.1时,可以把它当作一个无限边坡进行分析。

图 9-2表示一无限边坡示意图,滑动面位置在坡面下H深度处。

取一单位长度的滑动土条进行分析,作用在滑动面上的剪应力为,在极限平衡状态时,破坏面上的剪应力等于土的抗剪强度,即得式中N s =c/ γ H 称为稳定系数。

匀质粘性土体边坡稳定性计算

匀质粘性土体边坡稳定性计算

匀质粘性土体边坡稳定性计算边坡稳定性计算是土力学领域的一项重要内容,对于工程设计和施工具有重要的指导意义。

本文将对匀质粘性土体边坡稳定性计算方法进行详细介绍。

1.边坡稳定性计算的基本原理边坡稳定性计算是根据土体力学性质和力学平衡原理,对边坡的稳定性进行评价和判断。

边坡稳定性计算的基本原理包括力学平衡原理、摩擦角承载力理论、土体抗剪强度理论等。

2.边坡稳定性计算的步骤(1)确定边坡几何参数:包括边坡坡度、高度、底宽等几何参数,根据实际工程情况确定。

(2)确定土体力学参数:包括土体重度、内摩擦角、黏聚力等力学参数,通过室内试验或现场取样测试获得。

(3)计算边坡最大抗滑力:根据土体力学参数和边坡几何参数,使用土体力学理论计算边坡最大抗滑力。

(4)计算边坡滑动力:根据边坡上部土体的自重和外荷载,使用力学平衡原理计算边坡滑动力。

(5)判断边坡稳定性:比较边坡最大抗滑力与边坡滑动力的大小,判断边坡的稳定性。

3.土体力学参数的确定确定土体力学参数是边坡稳定性计算的重要基础,一般可以通过以下方法获得:(1)室内试验:如承压试验、剪切试验等,可以通过试验数据获得土体的内摩擦角、黏聚力等力学参数。

(2)现场取样测试:通过实际工程现场取样,并进行室内试验,获得土体力学参数。

(3)经验法:根据相似工程经验数据,利用统计分析方法获得土体力学参数的估计值。

4.边坡稳定性计算中的应力分析边坡稳定性计算中的应力分析是一个重要的环节,主要包括荷载分析和力的平衡分析。

荷载分析是指对边坡上的各种荷载进行分析,包括自重、活动荷载等。

力的平衡分析是指根据力学平衡原理,对边坡上部和下部的各种力进行分析,包括抗滑力、摩阻力、重力等。

5.边坡稳定性计算中的安全系数安全系数是衡量边坡稳定性的一个重要指标,一般通过比较抗滑力和滑动力的大小来确定。

安全系数的计算是边坡稳定性计算的最终结果,一般要求安全系数大于1,即抗滑力大于滑动力,以确保边坡的稳定。

土质边坡破坏模式与稳定性计算公式

土质边坡破坏模式与稳定性计算公式
土质边坡破坏模式与稳定性计算公式
添加副标题
主讲人: ppt制作与资料汇总: 资料收集:
CONTENTS
目 录
1
WORKREVIEW
添加标题
2
添加标题
UNDERWORK
一、土质边坡的结构类型
土质边坡:泛指具有倾斜面的土体,根据土体结构可分为三类:
类型
结构特征
稳定性影响因素
均质土边坡
整个坡体由均质土构成,基本不含节理、裂隙、没有贯通性的结构面
人为因素:人为因素的影响主要考虑施工步骤对边坡稳定性产生的影响,主要是坡形
四、稳定性影响因素分析
直线滑动面法:
01
松散的砂类土边坡,渗水性强,粘性差,边坡稳定主要靠其内摩擦力。失稳土体的滑动面近似直线状态,故直线滑动面法适用于砂类土:如下图所示,验算时,先通过坡脚或变坡点假设一直线滑动面,将边坡斜上方分割出下滑土楔体ABD,沿假设的滑动面AD滑动,其稳定系数K按下式计算(按边坡纵向单位长度计):
1、内在因素 临空条件:
四、稳定性影响因素分析
①、水对土体有软化作用,会使土体的抗剪强度降低,而且还会对结构面起到润滑作用;
2




②、在土体结构面中会形成一定的地下静水压力,随着时间的推移会促使滑面变形而使土体失稳;
③、会产生冻胀现象。所以水文条件也是影响边坡稳定性的重要因素。
④、增加了土体的重度,也就增加了土体的下滑力;
主要是水的作用,因水湿陷,或对边坡浸泡,水下渗使下垫隔水粘土层泥化等
1.崩塌; 2.张裂: 3.湿陷; 4.高或超高边坡可能出现高速滑坡
三、均质土边坡各种破坏模式
均质土边坡各种破坏模式 根据上表,可以看出土质边坡影响稳定性的因素主要是土体强度和水的作用,而产生的破坏形式以滑坡为多,崩塌和坍塌是开挖边坡过程中常见的(该处应该加上坡高、坡角、坡形的影响)

土质边坡破坏模式与稳定性计算公式

土质边坡破坏模式与稳定性计算公式
土体的下滑力; ⑤、当水发生渗流时会产 素 生动水压力;
四、稳定性影响因素分析
2、外在因素
地 震:地震作用对边坡的稳定性不利,主要是由于地震力的水平 分力,在横波的作用下,使边坡岩土体产生向临空面的拉力,易造 成边坡失稳破坏。另外,地震可能引起砂土液化、地下水位骤变、 坡体应力急剧改变、岩土体松动等多种不利因素,所以地震造成边 坡失稳的情况非常复杂,要做到合理的防护,还有待开展更多的地 震工程的研究。地震的发震时间、地点及地震震级都是很难准确预 测 人为因素:人为因素的影响主要考虑施工步骤对边坡稳定性产生的土体中的结构面对边坡 稳定性影响较大。边坡稳定系数K随结构面粘聚力 c和内摩擦系数tanφ的增大均呈对数关系增大,变 化较明显。结构面的强度与结构面的贯通度及面 内填充物相关。(不同土层面间的粘聚力和内摩 擦角怎么确定)
结构面倾角:结构面的存在,降低了土体的整体强度, 增大了坡体的变形性能,加强了土体的流变力学特性, 加深了岩土体的不均匀性、各向异性和非连续性等性 质。边坡的稳定系数随接触面倾角的增大而减小(稳 定系数与接触面倾角的tan值呈线性关系)
② 、坡度:坡度对边坡的稳定 性影响最大,边坡的稳定系数 随坡度增大呈幂函数减小,坡 度小于50 ° 时变化较大,大于 50 °变化趋缓。随着坡度的增 大,坡体位移与剪应变急剧增 加,发生突变,由此可将其视 为坡体即将发生破坏的判据;
③、坡形:边坡开挖时阶梯型 边坡比一坡到顶的边坡稳定性 好,并且台阶越宽、台阶数越 多稳定性越好,设计时结合造 价和旅工难易程度进行取舍。 坡表形态中,微凹型、直线型、 微凸型边坡稳定性依次减小。
坡面侵蚀
三、均质土边坡各种破坏模式
1、各种破坏形式形成条件及岩土结构
形成条件(即出现下列情况应采取措施)

土质边坡破坏模式与稳定性计算公式

土质边坡破坏模式与稳定性计算公式

② 、坡度:坡度对边坡的稳定 性影响最大,边坡的稳定系数 随坡度增大呈幂函数减小,坡 度小于50 ° 时变化较大,大于 50 °变化趋缓。随着坡度的增 大,坡体位移与剪应变急剧增 加,发生突变,由此可将其视 为坡体即将发生破坏的判据;
③、坡形:边坡开挖时阶梯型 边坡比一坡到顶的边坡稳定性 好,并且台阶越宽、台阶数越 多稳定性越好,设计时结合造 价和旅工难易程度进行取舍。 坡表形态中,微凹型、直线型、 微凸型边坡稳定性依次减小。
程度;
2.含水情况; 3.振动; 4.地表水及地下水
振动力作用下,易产生液化滑 坡; 2.管涌、流土; 3.坍塌和剥落
作用
以粉粒为主、质地均一。一 主要是水的作用,因 1.崩塌;
黄 土 边
般含钙量高,无层理,但柱 状节理发育,天然水含水量 低,干时坚固,部分黄土遇 水湿陷,有时呈固结状,有
水湿陷,或对边坡浸 2.张裂:
气候下施工,如雨季、寒冬季节。
1、地形条件:坡度一般要大于岩屑的休止角,要大于33°; 坍塌产生于易风化的土质边坡
坡地的相对高度大于50米时,可发生大型崩塌.2、地质条 和类土质边坡,尤其在膨胀土边坡
坍 塌
件:软弱面与坡面的倾向和倾角的关系不同,斜坡发生崩 或处于冻胀作用强烈区的边坡,一 塌的可能性也不一样.3、气候条件:温差较大,降水较多的 般发生在坡度大于20°时,随坡度 地区易发生崩塌.4、地震,强烈的融冰化雪.5、人工开挖边 增大发生坍塌的几率也越大,在暴
坡.。
雨季节,边坡表层岩土强度迅速降
低,也会促使坍塌破坏发生
我们实际中考虑滑坡的破坏性
三、均质土边坡各种破坏模式
• 2、土质边坡破坏机理
边坡的失稳破坏主要是由于边坡内所受的应力超过岩土体或 结构面的强度,从而导致边坡结构破坏。边坡变形表现为卸荷 回弹和蠕变两种主要方式。(具体破坏机理见下表。)

边坡稳定性计算方法

边坡稳定性计算方法

一、边坡稳定性计算方法在边坡稳定计算方法中,通常采用整体的极限平衡方法来进行分析。

根据边坡不同破裂面形状而有不同的分析模式。

边坡失稳的破裂面形状按土质和成因不同而不同,粗粒土或砂性土的破裂面多呈直线形;细粒土或粘性土的破裂面多为圆弧形;滑坡的滑动面为不规则的折线或圆弧状。

这里将主要介绍边坡稳定性分析的基本原理以及在某些边界条件下边坡稳定的计算理论和方法。

(一)直线破裂面法所谓直线破裂面是指边坡破坏时其破裂面近似平面,在断面近似直线。

为了简化计算这类边坡稳定性分析采用直线破裂面法。

能形成直线破裂面的土类包括:均质砂性土坡;透水的砂、砾、碎石土;主要由内摩擦角控制强度的填土。

图 9 - 1 为一砂性边坡示意图,坡高 H ,坡角β,土的容重为γ,抗剪度指标为c、φ。

如果倾角α的平面AC面为土坡破坏时的滑动面,则可分析该滑动体的稳定性。

沿边坡长度方向截取一个单位长度作为平面问题分析。

已知滑体ABC重 W,滑面的倾角为α,显然,滑面 AC上由滑体的重量W= γ(Δ ABC)产生的下滑力T和由土的抗剪强度产生的抗滑力Tˊ分别为:T=W · sina和则此时边坡的稳定程度或安全系数可用抗滑力与下滑力来表示,即为了保证土坡的稳定性,安全系数F s 值一般不小于 1.25 ,特殊情况下可允许减小到 1.15 。

对于C=0 的砂性土坡或是指边坡,其安全系数表达式则变为从上式可以看出,当α =β时,F s 值最小,说明边坡表面一层土最容易滑动,这时图9-1 砂性边坡受力示意图当 F s =1时,β=φ,表明边坡处于极限平衡状态。

此时β角称为休止角,也称安息角。

此外,山区顺层滑坡或坡积层沿着基岩面滑动现象一般也属于平面滑动类型。

这类滑坡滑动面的深度与长度之比往往很小。

当深长比小于 0.1时,可以把它当作一个无限边坡进行分析。

图 9-2表示一无限边坡示意图,滑动面位置在坡面下H深度处。

取一单位长度的滑动土条进行分析,作用在滑动面上的剪应力为,在极限平衡状态时,破坏面上的剪应力等于土的抗剪强度,即得式中N s =c/ γ H 称为稳定系数。

土质边坡破坏模式与稳定性计算公式

土质边坡破坏模式与稳定性计算公式

1、对于粘聚力为零的砂性均质 土边坡,发生滑坡破坏时,表现为 平面滑动,破坏面在截面上为一条 通过坡脚的直线;对于存在粘聚力 的粘性均质土边坡,则表现为圆弧 面滑动,破坏面在截面上为通过坡 脚的圆弧
2、黄土边坡坡体破坏形式主 要为滑坡和崩塌。滑坡或崩塌的形 成与边坡原始坡度有关:对于小于 50°的不稳定边坡,其破坏模式主 要是滑坡;51°~70°的不稳定边 坡破坏模式以滑坡为主,并伴有崩 塌;当边坡大于70°时,基本不发 生滑坡,主要破坏模式为崩塌。
沿圆弧形滑动面 滑移
岩土体沿莫一弱 面或朝向坡外的结构面 整体向下滑移
剪切-滑移。人工开挖增大坡角, 或地表水入渗使内摩擦角和内聚力降低, 达到临界值沿圆弧形滑动面滑移
拉裂、剪切-滑移。层面或贯通性 结构面形成滑动面,结构面临空,坡脚 岩层被切断或坡脚岩层挤压剪切
坍塌
坡体松弛带内的岩土由于震动、或侧 向卸荷、与坡面加载以及四季中时干时湿等 使松弛带内岩土的结合密实度在不断变化而 塌坡,塌至与其相适应的斜率为止
主要特征
影响稳定的主要 因素
可能的主要变形模式
以粘粒为主,一般干时坚硬, 1.矿物成分,特别 1.裂隙性粘土常沿光滑裂隙
粘 遇水膨胀崩解。某些粘土具 是亲水、膨胀、溶滤 面形成滑面,含膨胀性亲水矿
性 大孔隙性(如山西南部的粘土),性矿物含量;
物粘土易产生滑坡,巨厚层半
土 某些粘土甚坚固(如南方网纹 2.节理裂隙的发育 成岩粘土高边坡因坡脚蠕变可
② 、坡度:坡度对边坡的稳定 性影响最大,边坡的稳定系数 随坡度增大呈幂函数减小,坡 度小于50 ° 时变化较大,大于 50 °变化趋缓。随着坡度的增 大,坡体位移与剪应变急剧增 加,发生突变,由此可将其视 为坡体即将发生破坏的判据;
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

匀质粘性土体边坡稳定性计算中 极限平衡法与强度折减法数值分析比较1边坡的种类和滑面类型边坡按形成的原因大致可以分为三类:自然边坡、人工开挖边坡和人工填筑边坡;边坡按土层的种类也大致可以分为三种:匀质黏土边坡(人工填筑边坡一般为匀质边坡)、非匀质粘性土边坡(自然的黏砂性土边坡,或者人工开挖黏土边坡)和存在潜在软弱面或层面强度差异较大的边坡(如存在既有滑面的滑坡和上层为黏土层下为岩层的边坡)。

使边坡失稳的外因大致有:外荷载(地震、列车荷载、房屋和填土等)、重力和水的渗流、岩土的膨胀力等。

边坡失稳时滑动及滑面类型主要取决于边坡的外部荷载和边坡土层的类别,匀质粘土边坡的失稳滑面主要为圆弧型,非匀质边坡的失稳滑面主要是曲面型(复合型、对数螺旋型等),存在潜在软弱面或层面强度差异较大的边坡一般沿既有软弱面或者沿强度差异较大的两层面层间滑动(一般多为折线型,也有直线型)。

本将对匀质边坡进行有限差分强度折减法数值分析法和经典极限平衡法进行稳定性分析计算,并就安全系数的计算及粘性土边坡潜在滑动面的确定进行比较分析。

2匀质边坡的极限平衡法计算和数值模拟2.1匀质粘土边坡的极限平衡法概述粘性土边坡中,危险滑动面在土体的内部,常与圆弧面相似。

经典的基于圆弧滑动面的边坡稳定分析方法称为圆弧滑动法,属极限平衡法[1],本文对几种常用的经典分析方法进行简单概述如下。

(1)整体圆弧滑动法。

又称瑞典圆弧法,用于分析均质黏性土边坡的稳定性,即只能分析内摩擦角υ=0时的边坡稳定问题。

边坡稳定安全系数为抗滑力矩与滑动力矩之比,一般公式为:11(sin )nii r nsiii clM K M W α====∑∑(2)简单圆弧条分法。

又称瑞典条分法,仍假定滑动面为一个圆弧面。

为分析摩擦角υ>0时粘性土边坡的稳定性,将土坡分成若干个条块,但只考虑作用在条块上的重力、滑弧面上的法向力和切向抗滑力,忽略条块侧面法向力和切向力的作用。

一般公式为:11(cos tan )(sin )niii i r nsiii W cl M K M W αφα==+==∑∑(3)简化毕肖普(Bishop)法。

简化毕肖普法仍假设滑动面为圆弧,将滑动土体分为若干个条块。

在进行第i个条块的受力分析时,考虑条块侧面法向力的作用,但忽略切向力的作用。

一般公式为:1[tan ]sin i i i i iiii ic b W m K e W QRφα+=+∑∑∑,式中tan sin cos i ii i m Kφαα=+由以上简化毕肖普法的安全系数计算公式中可以看出,K 不是显函数,需采用试算迭代法求解。

(4)简布(Janbu )法。

简布法又称普遍条分法,其特点是对条块间正压力和切向力都给予考虑,由于滑动面不必是圆弧而可以是任何滑动面,其稳定安全系数也需要采用迭代方式求解(公式略)。

上面四种方法中,瑞典圆弧法和瑞典条分法由于不考虑条间力的影响,公式简单,计算速度快,但安全系数计算偏低10%~20%;后两种方法由于考虑条间力的作用,力学平衡较为合理,计算精度较高,但是计算公式较复杂,需要进行迭代试算,计算速度较慢。

本文采用Bishop 法进行匀质粘性土坡稳定性的极限平衡分析。

2.2粘性土边坡稳定性的Bishop 法分析Bishop 法作为一种条分法,发展了最初瑞典条分法将土坡按力与力矩极限平衡确定安全系数K 的思想,考虑了条间力的作用,并将土坡稳定安全系数确定为沿整个滑裂面的抗剪强度τf 与实际产生的剪应力τ之比,即:K=τ f /τ。

应用中还需要假定各土条之间的切向条间力均略去不计,即土条间的合力水平,这就是简化的Bishop 法。

Bishop 的公式推导及具体计算步骤请参见文献[2]。

现以新建铁路久长至永温线DK33+840~DK33+980段高填方路堤工点为例,采用 Bishop 法对铁路人工填筑路堤边坡进行稳定性分析: (1)工程地质概况DK33+840~DK33+980段路堤,长140m ,边坡填方最大高度25m ,站内为段内上覆坡残积(Q 4dl+el )红黏土,下伏基岩为寒武系中上统娄山关群(∈ol )白云岩夹泥质白云岩、角砾状白云岩。

地质建议基底岩土物理力学参数如下:岩土物理力学参数建议值(2)计算参数本工点选取填方边坡最高处DK33+970为代表断面,总填土高H=20m ,第一级填高12m 按1:1.75填筑,中间设2m 平台,第二级填高8m 按1:1.5填筑至坡顶,路基面宽度12.9m 。

坡脚换算角度 α=arctan(1/1.75)=29.74度,填土重度γ=20KN/m 3, υ=20度,c=10kpa ;填方基底为石灰岩地基(上部粘土考虑换填硬质岩W2),厚度为20m ,γ=25KN/m 3, 综合υ=60度。

荷载按铁路列车活载及轨道等静载换算为宽3.5m ,高3.2m 的土柱。

(3)计算过程边坡稳定性分析按照单线有列车荷载、双线有荷载、及双线均无列车荷载四种工况进行稳定性分析,将稳定性安全系数最小值作为设计值。

这里仅列出铁路双线均有列车行驶时的工况,计算简图如图1。

每条分宽度1m ,采用自动搜索最危险滑面的办法,搜索时的圆心步长0.5m ,半径步长0.5m 。

求得安全系数为1.01,最危险滑面滑动圆心O =(-0.040,51.700), 滑动半径R=51.7m ,最危险滑面见图2:Bishop 计算成果图。

从计算成果图上看,最危险滑面下部经过坡角附近,上部与第一个换算土柱结点相交。

X图2:Bishop 计算成果图2.2匀质边坡的数值模拟 边坡稳定的数值分析,可以得到到土体本身的全范围应力-应变关系,从微观上分析和计算边坡的稳定性。

基于有限差分的强度折减法,可以通过力和位移边界,应用合适的岩土模型和相应的屈服、破坏以及流动法则,能较准确的模拟匀质边坡的变形过程和应力分布,得到塑性位移等值线,根据塑性位移等值线图分析潜在滑移面及求得边坡安全系数。

这里利用2.1节所述算例进行数值模拟计算,填筑土体及下部岩石均采用摩尔-库仑屈服准则与非关联流动法则的弹塑性本构模型,应用世界著名的ITASCA 岩土咨询公司研发的Flac5.0进行建模分析。

Flac (fast lagrangian analysis for continus)是可以完成“拉格朗日分析”的“显式有限差分程序”,即使对准静态问题,FLAC 仍然求解完整的动力学方程,用局部阻尼及混合阻尼的方法吸收动能,以模拟系统的静态反应。

强度折减法简介强度折减法是有限元计算边坡稳定性的通用方法,计算时,首先选取初始折减系数Fs(通常为1),折减土体强度参数,将折减后的参数作为输入,进行有限元计算,若程序收敛,则土体仍处于稳定状态,然后再增加折减系数,直到不收敛为止,此时的折减系数即为边坡的稳定安全系数,此时的滑移面即为实际滑移面,这种方法称为土体强度折减系数法。

用折减系数摩擦角及粘聚力调整计算如下:折减的摩擦角υr=arctan(tan υ/Fs);折减的粘聚力Cr =C/Fs 。

其中υ和C 为土体初始力学参数。

但是目前有限元强度折减法中土体破坏的判断标准尚未统一,赵尚毅等人认为采用有限元计算时计算数值不收敛为判断土体破坏的依据,但是梁瑶等认为“有限元数值收敛时也不一定表明边坡处于安全状态,因此将计算的收敛性作为边坡失稳的判据不具有广泛的适用性”,另外还认为采用特征部位的位移突变性和塑性区的贯通性作为失稳判据受折减系数增幅的影响较大,会错过失稳的临界点。

采用Flac 的有限差分强度折减法计算时,由于系统本身求解的是动态方程,并采用双向折减系数进行逼近,较大值向下递减,较小值向上递增,当两个计算方向都收敛于同一数值时,认为该折减系数就是要求的安全系数。

这样就避免了上述各种失稳判据的不足和误差,计算结果是较为可信的。

根据《工程地质手册》第四版(以下简称《地质手册》)P160页表3-1-24,对照该填土层参数,取泊松比ν=0.3,压缩模量为Es=15Mpa ,推得弹性模量为022220.3(1)15110.3s E E νν⨯⨯=-⨯⨯)=(1---=0.743⨯15=11Mpa ,则土体弹性体积模量119.173(12)3(120.3)EK ν===⨯-⨯-⨯Mpa,土体弹性剪切模量G=112(1)2(10.3)Eν=⨯+⨯+=4.23Mpa ;根据《地质手册》P170页表3-1-43,对照石灰岩参数,取ν=0.3,静弹性模量为E=42.110⨯Mpa 。

根据Flac 的计算分析,求得边坡安全系数为1.08,图3及图4分别为边坡塑性剪切应变等值线图及塑性剪切应变率等值线图,图5为X 方向位移云图。

图3 边坡塑性剪切应变等值线图图4 边坡塑性剪切应变率等值线图图5 边坡X方向位移云图计算结果分析(1)从剪切变形等值线图及应变率等值线图上看,两等值线的脊线上应变及应变率最大,可以看成是填土最可能的滑移面。

该近似圆弧脊线与极限平衡法计算的最危险滑面基本相同,近一步验证了匀质边坡的圆弧型滑动类型及数值模拟强度折减法分析匀质边坡稳定性的有效性和合理性。

(2)从剪切变形等值线图及应变率等值线图上可以看出受影响的区域及受影响的大小。

但将两等值线图与X方向位移云图比较可以看出,边坡表层土体X方向发生的位移最大,但其应变值却较小,即剪切变形和X方向的位移相关性不大。

3结论及建议(1)极限平衡法和数值模拟强度折减法都是工程应用中可行的通用边坡稳定性定量分析方法,且极限平衡法在圆弧型滑动计算中安全系数的计算结果偏安全。

(2)极限平衡法由于本身力学分析上的限制,对某些边界力学条件进行假设或者忽略不计(如Bishop方法中的切向条间力),使得其力学方程本身不够准确,因此计算结果一般偏于安全。

(3)数值模拟强度折减法关注整个系统的力学平衡,且能够根据土体模型参数求得应力及应变结果,能够很好的反应土体应力破坏的影响区域和变形大小。

采用Flac进行边坡计算时,由于采用的是求解动态方程和双向强度系数法逼近,可以避免有限元强度折减法分析时各种失稳判据的不足和误差,计算结果是较为可信的。

(4)数值模拟强度折减法进一步验证了匀质边坡的圆弧型滑动类型。

(5)数值模拟计算时应注意土层本身的力学参数和模型的选用。

力学参数对数值模拟的影响很大(主要是粘聚力和内摩擦角),实际工程中需要通过现场及室内试验,求得相应工点所需的土层力学参数,并根据现场土层及外荷载对土层影响的实际情况选取适当的本构模型。

此外,数值模拟还受单元划分的大小,几何及位移边界条件等的影响。

所以,工程实践中应根据具体情况选择简布法、Bishop法及瑞典条分法等经典稳定性计算方法结合工程经验判断数值分析结果的合理性。

相关文档
最新文档