基于FLAC_D_强度折减法滑_省略_究_以三峡库区白果树古滑坡群为例_陈新泽.docVIP

——————————————————————————————————————————————— 基于FLAC_3D_强度折减法滑_省略_究_以三峡库区白果树古滑坡群为例_陈新泽 DOI:10.16030/j.cnki.issn.1000-3665.2008.02.017 ·24· 水文地质工程地质 2008年第2期　 基于FLAC强度折减法滑坡三维稳定性研究 ———以三峡库区白果树古滑坡群为例 陈新泽,唐辉明,杨有成,胡　斌,倪　俊 (1.中国地质大学(武汉)工程学院,武汉　430074;2.江苏省镇江市国土资源局,镇江　212000) 摘要:利用FLAC3D的内置fish语言基于强度折剪法对其进行了二次开发,编制了滑坡体安全系数计算程序。结合工程实例,将这种方法应用于已知滑动面条件下的滑坡稳定性分析中,得到滑坡体天然状况下的二维、三维安全系数,同时对内摩擦角和粘聚力参数进行了敏感性分析,对比并分析了传统极限平衡法和FLAC3D强度折减法的计算结果,得出以下结论:①利用二维FLAC3D强度折减法计算的滑坡安全系数略高于传统极限平衡法;②利用三维FLAC3D强度折减法所得滑坡安全系数略高于二维FLAC3D的计算结果;③滑坡安全系数都是随着内摩擦角和粘聚力的增大而增大,其中内摩擦角和粘聚力的变化对三维FLAC3D强度折减法比二维求解安全系数影响更大一些。然后模拟了滑体在天然条件下处于极限平衡状态的滑动方向,并得出天然状况下滑坡整体基本稳定,仅滑坡后缘可能出现局部失稳的分析结果。同时探讨了三维FLAC3D模型建立过程中滑坡外围不动体及滑带以下基岩所取范围大小、网格剖分粗细程度对计算安全系数的影响。关键词:滑坡稳定性;强度折减法;FLAC3D 中图分类号:P642.22 　文献标识码:A 　文章编号:1000-3665(2008)02-0024-06 1 1 1 1 2 3D 1　引言 目前边坡稳定性分析的主要方法可分为基于强度 [1] 折减技术的各类数值分析方法和基于极限平衡理论的各类条分法两类。极限平衡理论中,对于给定的滑动机构,需引入各种内力简化假定,消除超静定性而使问题近似地静定可解。根据极限分析概念,由此所得到的解答没有满足严格意义上静力许可内力场的基本条件。目前的各种数值分析方法,一般只能得出边坡应力、位移、塑性区及阻滑桩的位移、剪力和最大弯距。自然界很多滑坡滑动面形态复杂,甚至主滑方向都是变化的。现有的三维极限平衡方法只能考虑单一主滑方向,部分方法对滑动面形状有特殊要求 [3] [2] 衡法原理以此确定土坡的安全系数及相应的潜在破坏模式形成了一种新的研究方法。本文尝试了利用 3D FLAC的内置fish语言基于强度折剪法进行二次开发,编制了滑坡体安全系数计算程序,可以直接计算滑坡体天然状况下的二维、三维安全系数,也能够得到滑坡体应力、变形分布,模拟滑坡的实际滑移面及其主滑方向,进而综合分析滑坡三维稳定性,同时对利用FLAC强度折减法分析滑坡三维稳定性过程中的相关问题进行进一步研究。结合工程实例,研究了内摩擦 3D 角和粘聚力对利用FLAC强度折减法求解安全系数的影响。模拟滑体在天然条件下处于极限平衡状态的滑动方向,得到符合实际的分析结果。最后提出三维 3D FLAC模型建立过程中计算安全系数的影响因素。 3D [4] 。目 前对滑坡稳定性大多都集中于二维分析,再将二维分析的结果推广到三维情况,结果很难反映滑坡的真实稳定状态。 近年来,根据极限分析基本原理,或利用各类数值分析的应力计算成果,引入强度折减技术,套用极限平 收稿日期:2007-04-06;修订日期:2007-06-26 基金项目:国家自然科学基金项目;中国博士后基 金(20060400256);2007年度湖北省博士后基金项目(鄂博管办200625) 作者简介:陈新泽(1982-),男,硕士研究生,从事岩土体稳定性 评价与数值模拟研究。2　计算原理 2.1　FLAC基本原理 3D FLAC基本原理类同于离散单元法,但它却能像有限元那样适用于多种材料模式与边界条件的非规则 3D 区域的连续问题求解。在求解过程中,FLAC又采用了离散元的动态松弛法,不需求解大型联立方程组。下面以显式拉格朗日有限差分为例,介绍FLAC法的基本原理。对于四边形单元,有: u﹒i(1)(2)(a)(2)