用有限元强折减法进行边坡稳定分析.pptx
1 1、前言 边坡稳定分析是经典土力学最早试图解决而至今仍未圆满解决的课题，各种稳定分析方法在国内外水平大致相当。对于均质土坡，传统方法主要有：极限平衡法，极限分析法，滑移线场法等，就目前工程应用而言，主要还是极限平衡法，但需要事先知道滑动面位置和形状。对于均质土坡，可以通过各种优化方法来搜索危险滑动面，但是对于岩质边坡，由于实际岩体中含有大量不同构造、产状和特性的不连续结构面（比如层面、节理、裂隙、软弱夹层、岩脉和断层破碎带等），这就给岩质边坡的稳定分析带来了巨大的困难，传统极限平衡方法尚不能搜索出危险滑动面以及相应的稳定安全系数，而目前的各种数值分析方法，一般只是得出边坡应力、位移、塑性区，而无法得到边坡危险滑动面以及相应的安全系数。第1页/共53页 随着计算机技术的发展，尤其是岩土材料的非线性弹塑性有限元计算技术的发展，有限元强度折减法近来在国内外受到关注，对于均质土坡已经得到了较好的结论，但尚未在工程中实用，本文采用有限元强度折减法，对均质土坡进行了系统分析，证实了其实用于工程的可行性，对节理岩质边坡得到了坡体的危险滑动面和相应的稳定安全系数。该方法可以对贯通和非贯通的节理岩质边坡进行稳定分析，同时可以考虑地下水、施工过程对边坡稳定性的影响，可以考虑各种支挡结构与岩土材料的共同作用，为节理岩质边坡稳定分析开辟了新的途径。第2页/共53页2.有限元强度折减系数法的基本原理 进行强度折减非线性有限元分析要有一个过硬的非线性有限元程序和收敛性能良好的本够模型。因为收敛失败可能表明边坡已经处于不稳定状态，也可能仅仅是有限元模型中某些数值问题造成计算不收敛。第3页/共53页3.有限元强度折减系数法精度分析 3.1 岩土本构关系的影响 3.1.1 3.1.1 屈服准则的影响 3.1.2 3.1.2 流动法则的影响3.2 有限元法引入的误差3.3 边坡参数的影响第4页/共53页运筹学解决问题的过程 1）提出问题：认清问题。2）寻求可行方案：建模、求解。3）确定评估目标及方案的标准或方法、途径。4）评估各个方案：解的检验、灵敏性分析等。第5页/共53页运筹学解决问题的过程运筹学解决问题的过程 5）选择最优方案：决策。6）方案实施：回到实践中。7）后评估：考察问题是否得到完满解决。1）2）3）形成问题；4）5）分析问题：定性分析与定量分析相结合，构成决策。第6页/共53页3.1.1 3.1.1 屈服准则的影响屈服准则的影响用折减系数法求解实际边坡稳定问题时，通常将土体假设成理想弹塑性体，其中屈服准则常选用摩尔-库仑准则（M-C）、德鲁克普拉格准则（D-P）以及摩尔-库仑等面积圆准则。第7页/共53页M-C准则较为可靠，它的缺点在于三维应力空间中的屈服面存在尖顶和棱角的不连续点，导致数值计算不收敛，所以有时也采用抹圆了的M-C修正准则，它是用光滑连续曲线来逼进摩尔-库仑准则，此法虽然方便了数值计算，但往往公式复杂不利用实际应用。而D-P准则在偏平面上是一个圆，更适合数值计算。通常取M-C准则的外角点外接圆、内角点外接圆或其内切圆作为屈服准则，以利数值计算由徐干成、郑颖人（1990）提出的摩尔-库仑等效面积圆准则实际上是将M-C准则转化成近似等效的D-P准则形式。该准则要求平面上的摩尔-库仑不等边六角形与相同静水压力条件下D-P圆面积相等。计算表明它与摩尔-库仑准则十分接近。第8页/共53页 表1 各准则参数换算表 编号编号准则种类准则种类DP1DP1外角点外接外角点外接D-PD-P圆圆DP2DP2内角点外接内角点外接D-PD-P圆圆DP3DP3内切内切D-PD-P圆圆DP4DP4等面积等面积D-PD-P圆圆第9页/共53页算例分析表明：摩尔库仑等面积圆准（M-C EAC）则与简化Bishop法所得稳定安全系数最为接近（图1）。对有效算例（0）的误差进行统计分析可知，当选用 M-C EAC 准则时，误差的平均值为5.7%，且离散度很小（图2）。而外角点外接D-P圆准则的平均误差为29.5，同时采用内角点外接D-P圆准则、内切D-P圆准则准则所得计算结果的离散度非常大，均不可用。因此在数值分析中可用 M-C EAC 准则代替摩尔库仑准则。第10页/共53页图图1 1 折减系数曲线折减系数曲线 Fig.1 Fig.1 Reduction coefficient curveReduction coefficient curve第11页/共53页图图2 2 摩尔库仑等面积圆准准则的计算误差分布摩尔库仑等面积圆准准则的计算误差分布Fig.2Fig.2 The error distribution of computation by M-C EACThe error distribution of computation by M-C EAC 第12页/共53页3.1.2 3.1.2 不同流动法则的影响不同流动法则的影响有限元计算中，采用关联还是非关联流动法则，取决于值（剪胀角）：=，为关联流动法则,=0，为非关联流动法则.表表2 2 不同流动法则的影响不同流动法则的影响=45=45，C=40KPaC=40KPa，H H20m20m，DP4DP4准则准则材料参材料参数数=10=10 =17 =17=25=25非关联非关联0.8710.8711.105 1.105 1.3631.363关联关联0.8870.8871.1371.1371.4251.425相对误相对误差差0.0180.0180.0290.0290.0450.045第13页/共53页3.2 3.2 有限元法引入的误差有限元法引入的误差 3.2.1 网格的疏密 有限元单元网格划分第14页/共53页表3 网格疏密对计算结果的影响H H=20m=20m =45=45=17=17c c=10000Pa=10000Pa节点数节点数5775771111111122502250DP4DP40.6610.6610.6180.6180.5930.593简化简化BishopBishop法法0.5830.5830.5830.5830.5830.583(DP4-(DP4-Bishop)/BishopBishop)/Bishop0.1340.1340.0600.0600.0170.017第15页/共53页3.2.2 3.2.2 模型边界范围模型边界范围 表4 边界条件对折减系数的影相对相对边距比边距比0 00.50.51.01.01.51.52.02.02.52.53.03.0L/HL/H1.1291.1291.1241.1241.1241.1241.1201.1201.1221.1221.1211.1211.1291.129R/HR/H1.0971.0971.0781.0781.1211.1211.1221.1221.1221.1221.1201.1201.1231.123B/HB/H1.1061.1061.1171.1171.1201.1201.1311.1311.1241.1241.1321.1321.1311.131L坡脚到左端边界的距离（左边距），R坡顶到右端边界的距离（右边距），B坡脚到底端边界的距离（底边距），H坡高 第16页/共53页3.3 3.3 边坡参数的影响边坡参数的影响3.3.1 坡高H 3.3.2 坡角3.3.3 粘聚力C 3.3.4 摩擦角 第17页/共53页3.3.1 3.3.1 坡高坡高H H 的影响的影响表5 H为变量时的最小安全系数（节点数1190个）图5 H折减系数曲线 =45=45 c c=42KPa=42KPa=17=17H H(m)(m)10102020303040405050DP4DP41.7331.7331.1281.1280.9230.9230.8200.8200.7350.735简化简化BishopBishop法法1.6121.6121.0641.0640.8670.8670.7640.7640.6980.698(DP4-(DP4-Bishop)/BBishop)/Bishopishop0.0750.0750.0600.0600.0650.0650.0730.0730.0530.053第18页/共53页3.3.2 3.3.2 坡角坡角的影响的影响表6 为变量时的最小安全系数（节点数1210）图6 折减系数曲线 H H=20m=20mc c=42KPa=42KPa=17=17坡角坡角（）30303535404045455050DP4DP41.4551.4551.3231.3231.2141.2141.1281.1281.0441.044简化简化BishopBishop法法1.3981.3981.2691.2691.1561.1561.0641.0640.9870.987(DP4-(DP4-Bishop)/BBishop)/Bishopishop0.0410.0410.0430.0430.0500.0500.0600.0600.0580.058第19页/共53页3.3.3 3.3.3 粘聚力粘聚力C C的影响的影响表 7 c为变量时的最小安全系数（节点数1111个）图 7 C折减系数曲线 H H=20m=20m =45=45=17=17C C(KPa)(KPa)2020404060609090DP4DP40.7930.7931.1011.1011.3791.3791.7811.781简化简化BishopBishop法法0.7520.7521.0361.0361.3021.3021.6851.685(DP4-Bishop)/Bishop(DP4-Bishop)/Bishop0.0550.0550.0630.0630.0590.0590.0570.057第20页/共53页3.3.4 3.3.4 摩擦角摩擦角的影响的影响H=20m=45 C=42KPa()0.110253545DP10.5251.0441.7692.2543.051DP20.5250.9301.3321.5301.887DP30.4540.8481.2791.4991.870DP40.4770.8961.3961.6892.182简化Bishop法0.4940.8461.3161.6232.073(DP1-Bishop)/Bishop0.0630.2340.3440.3550.472(DP2-Bishop)/Bishop0.0630.0990.012-0.080-0.090(DP3-Bishop)/Bishop-0.0810.002-0.028-0.099-0.098(DP4-Bishop)/Bishop-0.0340.0590.0610.0410.053第21页/共53页4.4.均质土坡稳定分析 均质土坡，坡高 ，土容重 ，粘聚力 ，内摩擦角 ，求坡角 时边坡的稳定安全系数。第22页/共53页 坡角（度）3035404550有限元法（外接圆屈服准则）1.781.621.481.361.29有限元法（莫尔-库仑等面积圆屈服准则）1.471.341.221.121.06简化Bishop法1.3941.2591.1531.0620.99Spencer法1.4631.3181.2121.1151.04计算结果计算结果第23页/共53页从表中计算结果可以看出，采用外接圆屈服准则计算的安全系数比传统的方法大许多，而采用莫尔-库仑等面积圆屈服准则计算的结果与传统极限平衡方法（Spencer法)计算的结果十分接近，说明采用莫尔-库仑等面积圆屈服准则来代替莫尔-库仑不等角六边形屈服准则是可行的。通过4组计算方案（改变内摩擦角、内聚力、坡角、坡高H的值）共计106个算例的比较分析表明，用莫尔-库仑等面积圆屈服准则求得的安全系数与Bishop法的误差为4%-8%，与Spencer法的误差为0.5%-4%，说明了有限元强度折减法完全可以实用于土坡工程。第24页/共53页5.5.岩质边坡稳定分析 岩体中的结构面，根据结构面的贯通情况，可以将结构面分为贯通性、半贯通性、非贯通性三种类型。根据结构面的胶结和充填情况，可以将结构面分为硬性结构面（无充填结构面）和软弱结构面。由于岩体结构的复杂性，要十分准确地反映岩体结构的特征并使之模型化是不可能的，也没有必要使问题复杂化，基于这种考虑，对于一个实际工程来说，往往根据现场地质资料，根据结构面的长度、密度、贯通率，展布方向等着重考虑2-3组对边坡稳定起主要控制作用的节理组或其它主要结构面。第25页/共53页5.1 5.1 有限元模型极其安全系数的求解有限元模型极其安全系数的求解 （1）软弱结构面 岩体是弱面体，起控制作用的是结构面强度，对于软弱结构面，可采用低强度实体单元模拟，按照连续介质处理，材料本构关系采用理想弹塑性模型，屈服准则为广义米赛斯准则。第26页/共53页有限元模型以及变形后产生的塑性区第27页/共53页 （2）硬性结构面。无充填的硬性结构面，不能按照传统连续介质原理进行处理，本文采用美国ANSYS程序提供的无厚度接触单元来模拟硬性结构面的不连续性。第28页/共53页接触面的接触摩擦行为服从库仑定律:n安全系数:第29页/共53页5.2 5.2 折线型滑动面边坡稳定分析算折线型滑动面边坡稳定分析算第30页/共53页有限元强度折减法Spencer法C=160kPa,1.000.99C=160kPa,2.112.11C=320kPa,2.332.33C=160kPa,2.091.98C=03.082.94不同方法求得的稳定安全系数第31页/共53页5.3 5.3 具有一组平行节理面的岩质边坡算例 如图所示，一组软弱结构面倾角40度，间距10m(在求出初始滑动面后，可在滑动面附近将结构面间距加密）,岩体和结构面采用平面6节点三角形单元模拟，岩体以及结构面材料物理力学参数取值见表5.3.1。采用不同方法的计算结果见表5.3.2，其中极限平衡方法计算结果是在滑动面确定的情况下算出的。通过有限元强度折减计算,当有限元计算不收敛时,程序自动找出了滑动面.在一组平行的结构面中,只出现了一条滑动面,其余结构面没有出现塑性区和滑动。第32页/共53页 表5.3.1 计算采用物理力学参数材料名称重度弹性模量泊松比内聚力内摩擦角kN/m3MPaMPa度岩体25100000.21.038结构面17100.30.1224 表5.3.2 计算结果计算方法安全系数有限元法（外接圆屈服准则）1.26有限元法（等面积圆屈服准则）1.03极限平衡方法(解析解)1.06极限平衡方法(Spencer)1.06第33页/共53页5.4 5.4 具有两组平行节理面的岩质边坡算例 如图所示，两组方向不同的节理，贯通率100%，第一组软弱结构面倾角30度，平均间距10m，第二组软弱结构面倾角75度，平均间距10m，岩体以及结构面计算物理力学参数见表5.4.1。按照2维平面应变问题建立有限元模型.计算步骤同上，通过有限元强度折减，求得的滑动面如图5.4(a)所示,它是最先贯通的塑性区。塑性区贯通并不等于破坏，当塑性区贯通后塑性发展到一定程度，岩体发生整体破坏，同时出现第二、三条贯通的塑性区，如图5.4(b)，程序还可以动画模拟边坡失去稳定的过程，从动画演示过程可以看出边坡的破坏过程也整体破坏的过程。第34页/共53页材料名称重度弹性模量泊松比内聚力内摩擦角kN/m3MPaMPa度岩体25100000.21.038第一组节理17100.30.1224第二组节理17100.30.1224 表5.4.1 计算采用物理力学参数 表5.4.2 计算结果 计算方法安全系数有限元法（外接圆屈服准则）1.62有限元法（等面积圆屈服准则）1.33极限平衡方法(Spencer)1.36第35页/共53页通过有限元强度折减，求得的滑动面如图6-(a)所示,它是最先贯通的塑性区。塑性区贯通并不等于破坏，当塑性区贯通后塑性发展到一定程度，岩体发生整体破坏，同时出现第二、三条贯通的塑性区，如图6-(b)，程序还可以动画模拟边坡失去稳定的过程，从动画演示过程可以看出边坡的破坏过程也整体破坏的过程。图5.4(a)首先贯通的滑动面 图5.4(b)滑动面继续发展第36页/共53页第37页/共53页5.5 5.5 带软弱夹层的土坡稳定分析算例带软弱夹层的土坡稳定分析算例这个算例最早是由Frelund和Krahn(1977)提出的，随后被广泛引证。该土坡在坡底1m深处含有一个0.5m厚的软弱夹层。第38页/共53页计算方法计算方法安全系数安全系数澳大利亚莫纳什大学土坡稳定澳大利亚莫纳什大学土坡稳定分析软件分析软件GWEDGEMGWEDGEM1.2881.288陈祖煜和邵长明陈祖煜和邵长明1.2421.242DonaldDonald和和GiamGiam1.271.27加拿大边坡计算程序加拿大边坡计算程序slope/wslope/w（Spencer)Spencer)1.2521.252本文有限元强度折减法本文有限元强度折减法1.2561.256不同方法求得的稳定安全系数第39页/共53页用有限元强度折减法求得的滑动面第40页/共53页5.6 福宁高速公路福宁高速公路A15-2标段二埔塘标段二埔塘2号深路堑高边号深路堑高边坡算例坡算例 碎石土强风化晶屑凝灰岩堆土弱风化晶屑凝灰岩未开挖前的原始地貌（断面一）第41页/共53页福鼎至宁德高速公路A15-2标段K102+720-K102+900右侧路堑高边坡位于霞浦县盐田乡二铺村西侧，边坡设计最高为50m，长180m，工点附近属丘陵剥蚀地貌，地形陡峭，路堑穿过的山体平均坡度大于25度。该路堑的施工是先清除表层植被，设置坡顶截水沟后从上向下逐级开挖，当挖到第二级边坡后（标高约103m，路基设计标高85m），2000年10月12日坡顶开始出现裂缝；10月15日，k102+740-k102+770段第五级边坡跨塌，此后坡顶裂缝发展较快，裂缝达到10-50cm。自出现裂缝后便停止施工，进行了详细的地质勘察，然后修改设计，清除已经发生滑坡破碎体，放缓边坡坡度，同时进行防滑加固，主要措施为预应力锚索格子梁加固，修筑截排水天沟。第42页/共53页原设计开挖断面 第43页/共53页修改设计后开挖断面及加固措施 第44页/共53页不同工况下的稳定安全系数稳定安全系数稳定安全系数未开挖前未开挖前1.311.31按照原设计开挖按照原设计开挖0.920.92按照修改设计开挖（未加锚索）按照修改设计开挖（未加锚索）1.061.06按照修改设计开挖（加锚索）按照修改设计开挖（加锚索）2.452.45第45页/共53页按照原设计开挖后的滑动面 滑动面第46页/共53页按照修改设计进行二次开挖没有支护情况下的滑动面第47页/共53页6.6.结论 目前对复杂节理岩质边坡的稳定分析尚没有好的办法,传统的极限平衡方法无法得到岩质边坡的滑动面及其稳定安全系数，而各种数值分析方法只能算出应力、位移、塑性区等,无法判断边坡的稳定安全系数以及相应的滑移面.利用非线性有限元强度折减系数法可以由程序自动求得边坡的危险滑动面以及相应的稳定安全系数，通过算例分析表明了此法的可行性，为岩质边坡稳定分析开辟了新的途径。所求安全系数的大小与所采用的屈服准则有关，莫尔-库仑屈服准则是目前边坡分析中被广泛采用的屈服准则，但是莫尔-库仑准则的屈服面存在尖顶和棱角，给数值计算带来困难。算例表明采用徐干成、郑颖人（1990）提出的莫尔-库仑等面积圆屈服准则求得的稳定安全系数与简化Bishop法的误差为4-8%，与Spencer法的误差为0.5-4%，证实了其实用于土坡工程的可行性.第48页/共53页如果使有限元计算保持足够的计算精度，那么有限元法较传统的方法具有如下优点：(a)能够对具有复杂地貌、地质的边坡进行计算;(b)考虑了土体的非线性弹塑性本构关系，以及变形对应力的影响；(c)能够模拟土坡的失稳过程及其滑移面形状。由图可见滑移面大致在水平位移突变的地方，也是在塑性区塑性发展最充分的地方，呈条带状；第49页/共53页(d)能够模拟土体与支护的共同作用，图7为无锚杆（锚杆单元被杀死）时边坡稳定安全系数为1.1，图8为有锚杆支护时安全系数为1.5，且塑性区后移。(e)求解安全系数时，可以不需要假定滑移面的形状，也无需进行条分。不加锚杆时的塑性区 加锚杆时的塑性区第50页/共53页 均质土坡破坏过程均质土坡破坏过程第51页/共53页谢谢大家！谢谢大家！第52页/共53页感谢您的观看！第53页/共53页
