匀质粘性土体边坡稳定性计算.docx
- 文档编号:17613022
- 上传时间:2023-07-27
- 格式:DOCX
- 页数:11
- 大小:580.66KB
匀质粘性土体边坡稳定性计算.docx
《匀质粘性土体边坡稳定性计算.docx》由会员分享,可在线阅读,更多相关《匀质粘性土体边坡稳定性计算.docx(11页珍藏版)》请在冰点文库上搜索。
匀质粘性土体边坡稳定性计算
匀质粘性土体边坡稳定性计算中
极限平衡法与强度折减法数值分析比较
1边坡的种类和滑面类型
边坡按形成的原因大致可以分为三类:
自然边坡、人工开挖边坡和人工填筑边坡;边坡按土层的种类也大致可以分为三种:
匀质黏土边坡(人工填筑边坡一般为匀质边坡)、非匀质粘性土边坡(自然的黏砂性土边坡,或者人工幵挖黏土边坡)和存在潜在软弱面或层面强度差异较大的边坡(如存在既有滑面的滑坡和上层为黏土层下为岩层的边坡)。
使边坡失稳的外因大致有:
外荷载
(地震、列车荷载、房屋和填土等)、重力和水的渗流、岩土的膨胀力等。
边坡失稳时滑动及滑面类型主要取决于边坡的外部荷载和边坡土层的类另I」,匀质粘土边坡的失稳滑面主要为圆弧型,非匀质边坡的失稳滑面主要是曲面型(复合型、对数螺旋型等),存在潜在软弱面或层面强度差异较大的边坡一般沿既有软弱面或者沿强度差异较大的两层面层间滑动(一般多为折线型,也有直线型)。
本将对匀质边坡进行有限差分强度折减法数值分析法和经典极限平衡法进行稳定性分析计算,并就安全系数的计算及粘性土边坡潜在滑动面的确定进行比较分析。
2匀质边坡的极限平衡法计算和数值模拟
2.1匀质粘土边坡的极限平衡法概述
的边坡稳定分析方法称为进行简单概述如下O
弧滑动法,属极限平衡法⑴,本文对几种常用的经典分析方法
粘性土边坡中,危险滑动面在土体的内部,常与
弧滑动面
(1)整体圆弧滑动法。
又称瑞典圆弧法,用于分析均质黏性土边坡的稳定性,即只能分析
内摩擦角u=o时的边坡稳定问题。
边坡稳定安全系数为抗滑力矩与滑动力矩之比,
\'cli
i土
n
1(Wisin:
i)
u>0时粘性土边坡的稳定性,将土坡分成若干个条块,但只考虑作用在条块上的重力、滑
弧面上的法向力和切向抗滑力,忽略条块侧面法向力和切向力的作用。
一般公式为:
Mr
K
mS
n
Z(Wi§cosctjgan+cli)i±
、(Wi:
sin:
i)
(3)简化毕肖普(Bishop)法。
简化毕肖普法仍假设滑动面为弧,将滑动土体分为若干个条块。
在进行第1个条块的受力分析时,考虑条块侧面法向力的作用,但忽略切向力的作用。
—般公式为:
,式中rrii=cos■-)
IIeiK
、「WiSin:
i・Qi・
R
由以上简化毕肖普法的安全系数计算公式中可以看出,K不是显函数,需采用试算迭代法求解。
(4)简布(Janbu)法。
简布法又称普遍条分法,其特点是对条块间正压力和切向力都给予考
力学平衡较为合理,计算精度较高,但是计算公式较复杂,需要进行迭代试算,计算速度较慢。
本文采用Bishop法进行匀质粘性土坡稳定性的极限平衡分析。
2.2粘性土边坡稳定性的Bishop法分析
Bishop法作为一种条分法,发展了最初瑞典条分法将土坡按力与力矩极限平衡确定安全系数K的思想,考虑了条间力的作用,并将土坡稳定安全系数确定为沿整个滑裂面的抗剪强度”与实际产生
的剪应力T之比,即:
K=Tf/To应用中还需要假定各土条之间的切向条间力均略去不计,即土条间的合力水平,这就是简化的Bishop法。
Bishop的公式推导及具体计算步骤请参见文献[2]。
现以新建铁路久长至永温线DK33+84〜DK33+980段高填方路堤工点为例,采用Bishop法对
铁路人工填筑路堤边坡进行稳定性分析:
(1)工程地质概况
DK33+840-DK33+980段路堤,长140m边坡填方最大高度25m站内为段内上覆坡残积(Q4W)红黏土,下伏基岩为寒武系中上统娄[」」尖群(€ol)白云岩夹泥质白云岩、角砾状白云岩。
地质建议基底岩土物理力学参数如下:
岩土物理力学参数建议值
层号及成因
岩土名称
稠度或风化程度
天然密度P(g/cm3)
凝聚力
C
(kPa
内摩擦角
U
(°)
基底摩擦系数f
基本承载力
do
(KPa)
v3乞严
红黏土
硬塑
1.9
35
15
0.3
180
<4>
€ol
白云岩夹泥
质白云岩、
角砾状白云
W3
2.3
/
50
0.5
400
W2
2.5
/
60
0.6
600
⑵计算参数
本工点选取填方边坡最高处DK33+970为代表断面,总填土高第一级填高12m按1:
1.75填筑,
arctan(171.75)=29.74度'填土重度丫-20KN/riAu=20&'
c=10kpa;填方基底为石灰岩地基(上部粘土考虑换填硬质岩W2),厚度为20m丫二25KN/m,综合u二60度。
荷载按铁路列车活载及轨道等静载换算为宽3.5m,3.2m的土柱。
⑶计算过程
边坡稳定性分析按照单线有列车荷载、双线有荷载、及双线均无列车荷载四种工况进行稳定性分析,将稳定性安全系数最小值作为设计值。
这里仅列出铁路双线均有列车行驶时的工
0.040,51.700),滑动半径R=51.7m最危险滑面见图2:
Bishop计算成果图。
从计算成果图
上看,最危险滑面下部经过坡角附近,上部与第一个换算土柱结点相交。
图1:
Bishop计算简图
图2:
Bishop计算成果图
2.2匀质边坡的数值模拟
从微观上分析和计算边
边坡稳定的数值分析,可以得到到土体本身的全范围应力■应变尖系,
坡的稳定性。
基于有限差分的强度折减法,可以通过力和位移边界,应用合适的岩土模型和相应的屈服、破坏以及流动法则,能较准确的模拟匀质边坡的变形过程和应力分布,得到塑性位移等值
线,根据塑性位移等值线图分析潜在滑移面及求得边坡安全系数。
这里利用2.1节所述算例进行数值模拟计算,填筑土体及下部岩石均采用摩尔■库仑屈服准则与非矣联流动法则的弹塑性本构模型,应用世界著名的ITASCA岩土
咨询公司研发的Flac5.0进行建模分析。
Flac(fastlagrangiananalysisforcontinus)是可以完成“拉格朗日分析”的“显式有限差分程序”,即使对准静态问题,FLAC仍然求解完整的动力学方程,用局部阻尼及混合阻尼的方法吸收动能,以模拟系统的静态反应。
强度折减法简介强度折减法是有限元计算边坡稳定性的通用方法,计算时,首先选取初始折减系数Fs(通常为
1),折减土体强度参数,将折减后的参数作为输入,进行有限元计算,若程序收敛,则土体仍处于稳定状态,然后再增加折减系数,直到不收敛为止,此时的折减系数即为边坡的稳定安全系数,此时的滑移面即为实际滑移面,这种方法称为土体强度折减系数法。
用折减系数摩擦角及粘聚力调整计算如下:
折减的摩擦角uuarctan(tanu/Fs)折减的粘聚力Cr二C/Fs。
其中u和C为土体初始力学参数。
但是目前有限元强度折减法中土体破坏的判断标准尚未统一,赵尚毅等人
认为采用有限元计算时计算数值不收敛为判断土体破坏的依据,但是梁瑶等认为“有限元数值收敛时也不一定表明边坡处于安全状态,因此将计算的收敛性作为边坡失稳的判据不具有广泛的适用性”,另外还认为采用特征部位的位移突变性和塑性区的贯通性作为失稳判据受折减系数增幅的影响较大,会错过失稳的临界点。
采用Flac的有限差分强
度折减法计算时,由于系统本身求解的是动态方程,并采用双向折减系数进行逼近,较大值向下递减,较小值向上递增,当两个计算方向都收敛于同一数值时,认为该折减系
数就是要求的安全系数。
这样就避免了上述各种失稳判据的不足和误差,计算结果是较为可信的。
根据《工程地质手册》第四版(以下简称《地质手册》)P160页表3-1—24,对照该填土层参数,取泊松比9=0.3,压缩模量为Es=15Mpa推得弹性模量为
2»'2
e()al•叩
2
20.3
-0.74315=11Mpa
则土体弹性体积模量K
E
=9.17Mp3,
3(「2・)3(1-20.3)
APB
-1—43,对照石灰岩参数,
取v~0.3,静弹性模量为E=2.1104Mpa
根据Flac的计算分析,
求得边坡安全系数为1・08,S3及图4分别为边坡塑性剪切应变等值线图
18-Mar-0817:
51
step4160
-3217E+00 FactorofSafety1.03 Max.shearstrainincrement Contourintervai=2.50E-O2 Minimum: O.OOE+OD Maximum: 2.00E-01 Boundaryplot state-owedIcompanyguiyang CI.SOCI 1•&M 2.500 a.500 5.5M (TOM) 图3边坡塑性剪切应变等值线图 LEGEND ia-Mar-0017: 21 step4160 ■3.217E+D0cxc6112E+01 FactorofSafety1.06 Max.shearstrain-rate Contourinterval-1.00E-05 Minimum: 0.00E+00 Maximum: 7.00E-05 Boundaryplat LUU—L4JXTkJ 1E1 state-cwedcompanyguiyang 0.5DO 3.5M 1.5OD 5.5DQ 率等值线图 18-Mar-0617: 21 step41SO ■3.217E+D0CXC6 112E+Q1-1,717E+O1 contours -1 OOE+O O -8.0GE-Q1 -6.00E-01• 4ODE• Contourinterval-2.0QE-Q1匚ac+crcfQafotv/1C口 state-owedcompanyguiyang 0.500 1.500 3.500 3.MM 1.50D 5.5W 图5边坡X方向位移云图计算结果分析 (1)从剪切变形等值线图及应变率等值线图上看,两等值线的脊线上应变及应变率最大,可 以看成是填土最可能的滑移面。 该近似圆弧脊线与极限平衡法计算的最危险滑面基本相同,近一步验证了匀质边坡的圆弧型滑动类型及数值模拟强度折减法分析匀质边坡稳定性的有效性和合理性。 (2)从剪切变形等值线图及应变率等值线图上可以看出受影响的区域及受影响的大小。 但将两等值线图与X方向位移云图比较可以看出,边坡表层土体X方向发生的位移最大,但其应变值却较小,即剪切变形和X方向的位移 相尖性不大。 3结论及建议 (1)极限平衡法和数值模拟强度折减法都是工程应用中可行的通用边坡稳定性定量分析方法,且极限平衡法在圆弧型滑动计算中安全系数的计算结果偏安全。 (2)极限平衡法由于本身力学分析上的限制,对某些边界力学条件进行假设或者忽略不计(如Bishop方法中的切向条间力),使得其力学方程本身不够准确,因此计算结果一般偏于安 全。 (3) 数值模拟强度折减法尖注整个系统的力学平衡,且能够根据土体模型参数求得应力及 应变结果,能够很好的反应土体应力破坏的影响区域和变形大小。 米用 Flac进行边坡计算时,由于采用的是求解动态方程和双向强度系数法逼近,可以避免有限元强度折减法分析时各种失稳判据的不足和误差,计算结果是较为可信的。 (4)数值模拟强度折减法进一步验证了匀质边坡的圆弧型滑动类型。 (5)数值模拟计算时应注意土层本身的力学参数和模型的选用。 力学参数对数值模拟的影响 很大(主要是粘聚力和内摩擦角),实际工程中需要通过现场及室内试验,求得相应工点所需的土层力学参数,并根据现场土层及外荷载对土层影响的实际情况选取适当的本构模型。 此外,数值模拟还受单元划分的大小,几何及位移边界条件等的影响。 所以,工程实践中应根据具体情况选择简布法、Bishop法及瑞典条分法等经典稳定性计算方法结合工程经验判断数值分析结果的合理性。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 粘性 土体边坡 稳定性 计算
![提示](https://static.bingdoc.com/images/bang_tan.gif)