有限元强度折减法在高速公路滑坡稳定分析中的运用.docVIP

有限元强度折减法在高速公路滑坡稳定分析中的运用 　　摘要：随着有限元理论的广泛应用和计算机技术的发展，有限元极限分析方法在滑坡稳定性分析中日益受到重视。本文采用有限元强度折减法与传统极限分析法结合，分析某高速公路K18+200～300路段的滑坡稳定性。结果表明，有限元强度折减法与传统方法相比，不需进行滑面假设，就能获得边坡的滑面及安全系数，能够对具有复杂地貌、地质的高速公路滑坡进行失稳计算，可为滑坡治理方案的确定提供依据。 　　关键词：滑动面安全系数 滑坡稳定性有限元强度折减法 极限状态 　　[中图分类号]TU457文献标识号A 　　 　　一、前言 　　滑坡一般是指土方边坡在一定范围内整体沿某一滑动面向下或向外移动而丧失其稳定性。边坡的稳定，主要由土体的抗滑能力来保持。当土体下滑力超过抗滑力，边坡就会失去稳定而发生滑动。滑坡往往是在外界不利因素影响下触发和加剧的。这些外界因素往往导致土体剪应力的增加或抗剪强度的降低，使土体中剪应力大于土的抗剪强度而造成土体滑坡。 　　由于影响因素较多，精确地计算边坡稳定尚有困难，因此，在施工中一般工程目前多是综合考虑影响边坡稳定的各种因素，根据经验确定土方边坡，保证边坡大小，使坡顶荷载符合规范要求，或设置必要的支护，以防滑坡。由于滑坡计算的复杂性，包括要查清坡体的地质状况及其强度参数，同时又要有科学合理的分析方法，所以这个问题一直都没有妥善的解决。传统方法主要有：极限平衡法，极限分析法等，就目前工程应用而言，主要还是极限平衡法，但需要事先知道滑动面位置和形状〔1〕。 　　随着有限元理论的广泛应用和计算机的发展，有限元强度折减法近来在国内外受到关注，对于均质土坡的滑坡问题已经得到了较好的结论。该方法可以考虑地下水、施工过程对边坡稳定性的影响，可以考虑各种支挡结构与岩土材料的共同作用，为滑坡分析开辟了新的途径〔2-5〕。本文以某高速公路滑坡为例，采用有限元极限分析方法与极限分析方法结合，研究复杂边坡稳定性，研究方法可为复杂滑坡稳定分析提供参考。 　　二、计算原理 　　1.有限元强度折减法原理[3] 　　所谓强度折减法就是将土体的抗剪强度指标C和Ф,用一个折减系数Fs,如式(1)和(2)所示,然后用折减后的抗剪强度指标CF和ФF取代原来的抗剪强度指标C和Ф。直到其达到失稳临界状态为止,同时得到安全系数Fs。 　　 (1) 　　 (2) 　　式中,C是土体的粘聚力,Ф是土体的内摩擦角。 　　2. 屈服准则的选取 　　有限元强度折减法中岩土材料本构模型采用理想弹塑性模型。本构模型常选用摩尔-库仑准则(M-C)、Drucker-Prager准则以及摩尔-库仑等面积圆准则（5-6）。 　　Drucker-Prager 准则可表述为如下形式: 　　 (3) 　　式中I1,J2分别为应力张量的第一不变量和应力偏张量的第二不变量。α、k是与岩土材料内摩擦角φ和内聚力c 有关的常数,不同的α、k 在π平面上代表不同的圆。 　　M-C准则较为可靠，它的缺点在于三维应力空间中的屈服面存在尖顶和棱角的不连续点,导致数值计算不收敛。而D-P准则在偏平面上是一个圆,不存在尖顶处的数值计算问题,更适合数值计算。徐干成、郑颖人(1990)提出的摩尔-库仑等面积圆屈服准则代替传统摩尔-库仑准则实际上是将M-C准则转化成近似等效的D-P准则形式（6-7）。 　　3.有限元收敛准则 　　有限元强度折减法分析边坡稳定性的关键问题是如何根据有限元计算结果来判别边坡是否处于失稳状态。目前的失稳判据主要有两类: 　　1)在有限元计算过程中采用力和位移的不收敛作为边坡失稳的标志;2)以广义塑性应变或者等效塑性应变从坡脚到坡顶贯通作为边坡失稳的标志。 　　三、算例分析 　　1.工程概况 　　某高速公路K18+200～+300路段边坡施工正在进行时，下了一场大暴雨，造成边坡失稳，使该地段的斜坡坡面产生不同程度的变形，部分土体出现滑动，使修建的方格状水泥混凝土护坡设施和水沟受到破坏，如图1。 　　 　　图1 高速公路边坡失稳现状图 　　结合现场实际情况，计算采用的力学参数见表1。 　　表1 计算采用物理力学参数 　　材料 　　名称 饱和重度 弹性模量 泊松比 内聚力 内摩擦角 　　 kN/m3 Mpa kPa ° 　　残积粘土层 　　（滑体） 20 80 0.3 17.8 11 　　全风化砂岩 　　（滑床） 25 1500 0.2 100 28 　　C25砼 24 29×103 0.2 按线弹性材料处理 　　2.模型建立 　　计算采用的软件为美国ANSYS公司的大型有限元计算程序ANSYS，按照平面应变问题处理，岩土材料用6节点三角形平面单元PLANE2模拟。边界条件：底