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

1、Cib+Wj tan Q匀质粘性土体边坡稳定性计算中极限平衡法与强度折减法数值分析比较1边坡的种类和滑面类型边坡按形成的原因大致可以分为三类：自然边坡、人工开挖边坡和人工填筑边坡；边 坡按土层的种类也大致可以分为三种：匀质黏土边坡（人工填筑边坡一般为匀质边 坡）、 非匀质粘性土边坡（自然的黏砂性土边坡，或者人工幵挖黏土边坡）和存在潜在软弱面或 层面强度差异较大的边坡（如存在既有滑面的滑坡和上层为黏土层下为岩层的 边坡）。使 边坡失稳的外因大致有：外荷载（地震、列车荷载、房屋和填土等）、重力和 水的渗流、 岩土的膨胀力等。边坡失稳时滑动及滑面类型主要取决于边坡的外部荷载和边坡土层的类 别，匀质粘

2、土边坡的失稳滑面主要为圆弧型，非匀质边坡的失稳滑面主要 是曲面型（复合 型、对数螺旋型等），存在潜在软弱面或层面强度差异较大的边坡一般 沿既有软弱面或者 沿强度差异较大的两层面层间滑动（一般多为折线型，也有直线型）。本将对匀质边坡进 行有限差分强度折减法数值分析法和经典极限平衡法进行稳定性分 析计算，并就安全系数 的计算及粘性土边坡潜在滑动面的确定进行比较分析。2匀质边坡的极限平衡法计算和数值模拟弧滑动面2.1匀质粘土边坡的极限平衡法概述粘性土边坡中，危险滑动面在土体的内部，常与圆弧面相似。经典的基于的边坡稳定分析方法称为弧滑动法，属极限平衡法，本文对几种常用的经典分析方法进行简单概述如下O（

3、1）整体圆弧滑动法。又称瑞典圆弧法，用于分析均质黏性土边坡的稳定性，即只能分析内摩擦角U=0时的边坡稳定问题。边坡稳定安全系数为抗滑力矩与滑动力矩之比, 'clii 土n1 (Wisin：i)i 土（2）简单圆弧条分法。又称瑞典条分法,仍假定滑动面为一个弧面。为分析摩擦角u>0时粘性土边坡的稳定性，将土坡分成若干个条块，但只考虑作用在条块上的重力、滑弧面上的法向力和切向抗滑力，忽略条块侧面法向力和切向力的作用。一般公式为：Z (Wi §cos ctj gan + cl i)i ±mS、(Wi：sin : i)i=1弧，将滑动土体分为若（3）简化毕肖普（Bish

4、op）法。简化毕肖普法仍假设滑动面为干个条块。在进行第1个条块的受力分析时，考虑条块侧面法向力的作用，但忽略切向力 的作用。一般公式为：tan 3in imi，式中 rrii = cos - jII eiK、rWiSin : i- Qi-R由以上简化毕肖普法的安全系数计算公式中可以看出,K不是显函数，需采用试算迭 代法求解。(4) 简布(Janbu)法。简布法又称普遍条分法，其特点是对条块间正压力和切向力都 给予考虑，由于滑动面不必是圆弧而可以是任何滑动面，其稳定安全系数也需要采用迭代方 式求解(公式略)。上面四种方法中，瑞典圆弧法和瑞典条分法由于不考虑条间力的影响，公式简单，计算 速度快，但

5、安全系数计算偏低10%20%；后两种方法由于考虑条间力的作用，力学平衡较为合理，计算精度较高，但是计算公式较复杂，需要进行迭代试算，计算速度较 慢。本文采用Bishop法进行匀质粘性土坡稳定性的极限平衡分析。2.2粘性土边坡稳定性的Bishop法分析Bishop法作为一种条分法，发展了最初瑞典条分法将土坡按力与力矩极限平衡确定安全 系数K的思想，考虑了条间力的作用，并将土坡稳定安全系数确定为沿整个滑裂面的抗剪强 度"与实际产生的剪应力T之比，即：K=Tf/T。应用中还需要假定各土条之间的切向条间力 均略去不计，即土条间的合力水平，这就是简化的Bishop法。Bishop的公式推导及具

6、体计 算步骤请参见文献2。现以新建铁路久长至永温线DK33+84DK33+980段高填方路堤工点为例，采用 Bishop法对铁路人工填筑路堤边坡进行稳定性分析：(1) 工程地质概况DK33+840-DK33+980段路堤，长140m边坡填方最大高度25m站内为段内上覆坡 残积 (Q4W)红黏土，下伏基岩为寒武系中上统娄尖群( ol)白云岩夹泥质白云 岩、角砾状白 云岩。地质建议基底岩土物理力学参数如下：岩土物理力学参数建议值层号及成因岩土名称稠度或 风化程度天然密度P( g/cm3 )凝聚力C(kPa内摩擦角U(°)基底摩 擦系数f基本承载力do(KPa)<3乞严红黏土硬塑1.

7、935150.3180<4> ol白云岩夹泥质白云岩、角砾状白云W32.3/500.5400W22.5/600.6600(2)计算参数本工点选取填方边坡最高处DK33+970为代表断面，总填土高第一级填高12m 按1:1.75填筑,中间设2m平台，第二级填高8m按1:1.5填筑至坡顶，路基面宽度12.9m。 坡脚换算角度 a= arctan(171.75)=29.74 度填土重度 丫 二 20KN/riA u= 20 度时的圆心步长0.5m,半径步长0.5m。求得安全系数为1.01，最危险滑面滑动0=c=10kpa；填方基底为石灰岩地基（上部粘土考虑换填硬质岩W2）,厚度为20m

8、丫二25KN/m, 综合u二60度。荷载按铁路列车活载及轨道等静载换算为宽3.5m,3.2m的土柱。计算过程边坡稳定性分析按照单线有列车荷载、双线有荷载、及双线均无列车荷载四种工况进行稳定性分析，将稳定性安全系数最小值作为设计值。这里仅列出铁路双线均有列车行驶时的工搜索况，计算简图如图1。每条分宽度im采用自动搜索最危险滑面的办法，0.040,51.700）,滑动半径R=51.7m最危险滑面见图2： Bishop计算成果图。从计算成果图上看，最危险滑面下部经过坡角附近，上部与第一个换算土柱结点相交。匕4图1： Bishop计算简图图2： Bishop计算成果图2.2匀质边坡的数值模拟边坡稳定的

9、数值分析，可以得到到土体本身的全范围应力应变尖系， 从微观上分 析和计算边坡的稳定性。基于有限差分的强度折减法，可以通过力和位移边界，应用合适的 岩土模型和相应的屈服、破坏以及流动法则，能较准确的模拟匀质边坡的变形过程和应力分 布，得到塑性位移等值线，根据塑性位移等值线图分析潜在滑移面及求得边坡安全系数。这 里利用2.1节所述算例进行数值模拟计算，填筑土体及下部岩石均采用摩尔库仑屈服准则 与非矣联流动法则的弹塑性本构模型，应用世界著名的ITASCA岩土 咨询公司研发的Flac5.0进行建模分析。Flac （fast lagrangian analysis for continus）是可以 完成

10、“拉格朗日分析”的“显式有限差分程序”，即使对准静态问题，FLAC仍然求解完整 的动力学方程，用局部阻尼及混合阻尼的方法吸收动能，以模拟系统的静态反应。强度折减法简介强度折减法是有限元计算边坡稳定性的通用方法，计算时，首先选取初始折减系数 Fs（通常为1），折减土体强度参数，将折减后的参数作为输入，进行有限元计算，若程序收 敛，则土体仍处于稳定状态，然后再增加折减系数，直到不收敛为止，此时的折减系数即为 边坡的稳定安全系数，此时的滑移面即为实际滑移面，这种方法称为土体强度折减系数法。 用折减系数摩擦角及粘聚力调整计算如下：折减的摩擦角u r=arctan（tan u/Fs）折减的粘聚力Cr二C

11、/Fs。其中u和C为土体初 始力 学参数。但是目前有限元强度折减法中土体破坏的判断标准尚未统一，赵尚毅等人认为采用有限元计算时计算数值不收敛为判断土体破坏的依据，但是梁瑶等认为“有限 元数值收敛时也不一定表明边坡处于安全状态，因此将计算的收敛性作为边坡失稳的判据 不具有广泛的适用性”，另外还认为采用特征部位的位移突变性和塑性区的贯通性作为失稳 判据受折减系数增幅的影响较大，会错过失稳的临界点。采用Flac的有限差分强度折减法计算时，由于系统本身求解的是动态方程，并采用双向折减系数进行逼近，较大值 向下递减，较小值向上递增，当两个计算方向都收敛于同一数值时，认为该折减系数就是要求的安全系数。这样

12、就避免了上述各种失稳判据的不足和误差，计算结果是较为可 信的。根据工程地质手册第四版（以下简称地质手册）P160页表31 24,对照该 填土层参数，取泊松比9=0.3，压缩模量为Es=15Mpa推得弹性模量为2 »'2E0AL'叩220.3)151- 0.3二0.743 15=11 Mpa则土体弹性体积模量KE3(1-2.)=9.17 Mp3,3 (1-20.3)土体弹性剪切模量 G=4.23Mpa；根据地质手 P170页表32（1 ）APB-1 43,对照石灰岩参数，取v二0.3，静弹性模量为E=2.1 104 Mpa根据Flac的计算分析，求得边坡安全系数为1.0

13、8,S3及图4分别为边坡塑性剪切应变 等值线图及塑性剪切应变率等值线图，图5为X方向位移云图。01E 1state-owedl company guiyangCI.SOCI1&M2.500(TOM)a.5005.5MLEGEND图3边坡塑性剪切应变等值线图ia-Mar-00 17:21step 4160 3.217E+D0 cxc6 112E+01.717E+01 4.717E+01Factor of Safety 1.06 Max. shear strain-rateContour interval-1.00E-05Minimum: 0.00E+00Maximum: 7.00E-05

14、Boundary platLUUL4J XT IJ1E 1图4边坡塑 性剪切应变state-cwed company guiyang0.5DO3.5M1.5OD5.5DQ率等值线图18-Mar-08 17:51step 4160-3 217E+00 <JC< 6-1iaE+01 -1 717E+01、严 4.717E+D1Factor of Safety 1.03Max. shear strain incrementContour intervai= 2.50E-O2Minimum: O.OOE+ODMaximum: 2.00 E-01Boundary plot18-Mar-06

15、17:21step 41 SO 3.217E+D0 CXC 6112E+Q1 -1,717E+O1contours-1OOE+OO-8.0GE-Q1-6.00E-01 4 ODEContour interval- 2.0QE-Q1Factor of Safety 1 .OBstate-owed company guiyang0.5001.5003.5003.MM1.50D5.5W图5边坡X方向位移云图计算结果分析(1 )从剪切变形等值线图及应变率等值线图上看，两等值线的脊线上应变及应变率 最大，可以看成是填土最可能的滑移面。该近似圆弧脊线与极限平衡法计算的 最危险滑面基本相同，近一步验证了匀质

16、边坡的圆弧型滑动类型及数值模拟强 度折减法分析匀质边坡稳定性的有效性和合理性。(2) 从剪切变形等值线图及应变率等值线图上可以看出受影响的区域及受影响的大 小。但将两等值线图与X方向位移云图比较可以看出，边坡表层土体 X方向发生的位移最大，但其应变值却较小，即剪切变形和X方向的位移相尖性不大。3结论及建议(1) 极限平衡法和数值模拟强度折减法都是工程应用中可行的通用边坡稳定性 定量 分析方法，且极限平衡法在圆弧型滑动计算中安全系数的计算结果偏安全。(2) 极限平衡法由于本身力学分析上的限制，对某些边界力学条件进行假设或者 忽 略不计(如Bishop方法中的切向条间力)，使得其力学方程本身不够准确，因此计算 结果一般偏于安全。(3) 数值模拟强度折减法尖注整个系统的力学平衡，且能够根据土体模型参数求 得应力及应变结果，能够很好的反应土体应力破坏的影响区域和变形大小。米用Flac进行边坡计算时，由于采用的是求解动态方程和双向强度系数法逼近， 可以避 免有限元强度折减法分析时各种失稳判据的不足和误差，计算结果