用有限元强度折减法进行边坡稳定分析_郑颖人.pdf

《用有限元强度折减法进行边坡稳定分析_郑颖人.pdf》由会员分享，可在线阅读，更多相关《用有限元强度折减法进行边坡稳定分析_郑颖人.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2002 年 10 月第 4 卷第 10 期中国工程科学Engineering ScienceOct12002Vol14 No110研 究 报 告用有限元强度折减法进行边坡稳定分析郑颖人,赵尚毅,张鲁渝(后勤工程学院军事土木工程系,重庆 400041)摘要 通过对边坡非线性有限元模型进行强度折减,使边坡达到不稳定状态时,非线性有限元静力计算将不收敛,此时的折减系数就是稳定安全系数,同时可得到边坡破坏时的滑动面。传统条分法无法获得岩质边坡的滑动面与稳定安全系数。该方法开创了求岩质边坡滑动面与稳定安全系数的先例。文章对此法的计算精度以及影响因素进行了分析。算例表明采用摩尔-库仑等面积圆屈服准则求得

2、的稳定安全系数与简化 Bishop 法的误差为 3%8%,与 Spencer 法的误差为 1%4%,证实了其实用于工程的可行性。关键词 边坡稳定分析;有限元强度折减法;屈服准则 中图分类号TU457 文献标识码A 文章编号1009-1742(2002)10-0057-05 收稿日期 2002-06-10 作者简介 郑颖人(1933-),男,浙江镇海市人,中国工程院院士,后勤工程学院教授,博士生导师1 引言 西部开发是我国实现地区平衡发展和可持续发展的重大战略举措。然而,我国西部地区山高坡陡、沟壑纵横,城市建筑依山而立,公路、铁路翻山越岭,复杂多变的地形地貌决定了我国西部开发将面临大量滑(边)坡

3、工程,滑坡与边坡事故日益增多。例如 2001 年,重庆市云阳县就发生两次大型滑坡,重庆市武隆边坡失稳坍塌造成 79人死亡。重庆市近 20 年来累计发生地质灾害 3133 万处,仅2000 年就发生 6371 处,受灾 19133 万人,倒塌房屋 8168 万间,直接经济损失 7167 亿元。重庆市已经成为地质灾害的重灾区,尤其是随着三峡库区蓄水和新兴城市的建设,有可能诱发更大的地质灾害,隐患无穷。现在,国务院已经决定拨款 40 亿元,用于三峡库区地质灾害治理,仅重庆三峡库区计划的地质灾害治理工程就有 143 个。频发的地质灾害以及大量灾害治理资金的投入,使得滑(边)坡稳定性问题成为西部开发中的

4、热点与难点问题。边坡稳定分析是经典土力学最早试图解决而至今仍未圆满解决的课题,各种稳定分析方法在国内外水平大致相当。对于均质土坡,传统方法主要有:极限平衡法、极限分析法和滑移线场法等。就目前工程应用而言,主要还是极限平衡法,但需要事先知道滑动面位置和形状。对于均质土坡,可以通过各种优化方法来搜索危险滑动面。但是,对于岩质边坡,由于实际岩体中含有大量不同构造、产状和特性等不连续结构面(比如层面、节理、裂隙、软弱夹层、岩脉和断层破碎带等),给岩质边坡的稳定分析带来了巨大的困难。传统极限平衡方法尚不能搜索出危险滑动面以及相应的稳定安全系数。而目前的各种数值分析方法,一般只是得出边坡应力、位移、塑性区

5、,也无法得到边坡危险滑动面以及相应的安全系数。随着计算机技术的发展,尤其是岩土材料的非线性弹塑性有限元计算技术的发展,有限元强度折减法近来在国内外受到关注,对于均质土坡已经得到了较好的结论,但尚未在工程中实用。笔者采用有限元强度折减法1,对均质土坡进行了系统分析,证实了用于工程的可行性,得到了节理岩质边坡坡体的危险滑动面和相应的稳定安全系数。该方法可以对贯通和非贯通的节理岩质边坡进行稳定分析,同时可以考虑地下水、施工过程对边坡稳定性的影响,可以考虑各种支挡结构与岩土材料的共同作用,为岩质边坡稳定分析开辟了新的途径。2 有限元强度折减系数法基本原理 有限元强度折减系数法的基本原理是将坡体强度参数

6、:粘聚力 c 和内摩擦角 U值同时除以一个折减系数 Ftrial,得到一组新的 cc、U c值,然后作为新的资料参数输入,再进行试算,当计算不收敛时,对应的 Ftrial被称为坡体的最小稳定安全系数,此时坡体达到极限状态,发生剪切破坏,同时可得到坡体的破坏滑动面。cc=c/Ftrial,U c=arctan(tanU/Ftrial)3 有限元强度折减系数法精度分析311 屈服准则的选用安全系数大小与程序采用的屈服准则密切相关,不同的准则得出不同的安全系数。目前流行的大型有限元软件 ANSYS,以及美国 MSC 公司的MARC、PATRAN、NASTRAN 均采用了广义米赛斯准则:F=A I1+

7、J2=k,式中 I1,J2分别为应力张量的第一不变量和应力偏张量的第二不变量。A、k 是与岩土材料内摩擦角 U和内聚力 c 有关的常数,这是一个通用表达式,通过变换 A、k 的表达式就可以在有限元中实现不同的屈服准则,各准则的参数换算关系见表1。传统的极限平衡法采用摩尔-库仑准则,但是由于摩尔-库仑准则的屈服面为不规则六角形截面的角锥体表面,存在尖顶和菱角,给数值计算带来困难。为了与传统方法进行比较,本文采用了徐干成、郑颖人(1990)提出的摩尔-库仑等面积圆屈服准则(DP4)代替传统摩尔-库仑准则2。表 1 各准则参数换算表Table 1 T he relationship of diffe

8、rent yield criterions编号准则种类AkDP1 外角点外接 D-P 圆2sinU3(3-sinU)6ccosU3(3-sinU)DP2 内角点外接 D-P 圆2sinU3(3+sinU)6ccosU3(3+sinU)DP3 内切 D-P 圆sinU33+sin2U3ccosU33+sin2UDP4 等面积 D-P 圆2 3sinU2 3P(9-sin2U)6 3ccosU2 3P(9-sin2U)算例 分析 表明(表 2):DP4 准 则 与简 化Bishop 法所得稳定安全系数最为接近。通过对误差进行统计分析可知,当选用 DP4 准则时,误差的平均值为 517%,最大误差小

9、于 8%,且离散度很小,而 DP1 的平均误差为 2915%,同时采用DP2、DP3 准则所得计算结果的离散度也非常大。因此,在数值分析中可用 DP4 准则代替摩尔-库仑准则。表 2 不同屈服准则所得最小安全系数Table 2 Safety factor by different methodh=20 m B=45b c=42 kPaU/(b)01110253545DP10152511044117692125431051DP20152501930113321153011887DP30145401848112791149911870DP40147701896113961168921182简化 B

10、ishop 法0149401846113161162321073(DP1-Bishop)/Bishop0106301234013440135501472(DP2-Bishop)/Bishop010630109901012-01080-01090(DP3-Bishop)/Bishop-0108101002-01028-01099-01098(DP4-Bishop)/Bishop-0103401059010610104101053312 有限元计算精度分析边界范围的大小在有限元法中对计算结果的影响比传统极限平衡法表现得更为敏感,当坡角到左端边界的距离为坡高的 115 倍,坡顶到右端边界的距离为坡高

11、的 215 倍,且上下边界总高不低于 2 倍坡高时,计算精度最为理想。另外如果网格划分太粗,将会造成很大的误差,计算时必须考虑适当的网格密度。4 均质土坡稳定分析 均质土坡,坡高 h=20 m,土容重 C=25 kN/m3,粘聚力 c=42 kPa,内摩擦角 U=17 b,求坡角B=30b,35b,40b,45b,50 b时边坡的稳定安全系数。如图 1,按照平面应变建立有限元模型,边界条件为左右两侧水平约束,下部固定,上部为自由边界,计算结果见表 3。从表 3 看出,采用外接圆屈服准则计算的安全系数比传统的方法大许多,而采用摩尔-库仑等面积圆屈服准则计算的结果与传统极限平衡方法(Spencer

12、 法)计算的结果十分接近,说明采用摩尔-库仑等面积圆屈服准则来代替摩尔-库仑不等角58中国工程科学第 4 卷六边形屈服准则是可行的。图 1 有限元单元网格划分Fig11 FEM model表 3 用不同方法求得的稳定安全系数T able 3 Safety factor by different method方法坡角/(b)3035404550DP111781162114811361129DP411471134112211121106简化 Bishop 法11391126111511060199Spencer法11461132112111121104另外,通过 4 组计算方案(改变内摩擦角 U、

13、粘聚力 c、坡角 B、坡高 h 的值)共计 106 个算例的比较分析表明,用摩尔-库仑等面积圆屈服准则求得的安全系数与 Bishop 法的误差为 3%8%,与 Spencer 法的误差为 1%4%,说明了有限元强度折减法完全可以实用于土坡工程。5 岩质边坡稳定分析 根据岩体中结构面的贯通情况,可以将结构面分为贯通性、半贯通性、非贯通性三种类型。根据结构面的胶结和充填情况,可以将结构面分为硬性结构面(无充填结构面)和软弱结构面。由于岩体结构的复杂性,要十分准确地反映岩体结构的特征并使之模型化是不可能的,也没有必要使问题复杂化。基于这种考虑,对于一个实际工程来说,往往根据现场地质资料,根据结构面的

14、长度、密度、贯通率,展布方向等,着重考虑 2 3 组对边坡稳定起主要控制作用的节理组或其他主要结构面。岩体是弱面体,起控制作用的是结构面强度,对于软弱结构面,可采用低强度实体单元模拟,按照连续介质处理;对于无充填的硬性结构面可以采用无厚度的接触单元来模拟,安全系数的求解与均质土坡相同。结构面单元的设置会影响计算精度,一般可先以较大间距设置结构面,求出滑动面大概位置后再在滑动面附近将结构面加密,增加结构面单元以提高计算精度。511 具有一组平行节理面的岩质边坡算例如图 2 和图3 所示,一组软弱结构面倾角 40b,间距 10 m,岩体和结构面采用平面 6 节点三角形单元模拟。岩体以及结构面材料物

15、理力学参数取值见表 4。采用不同方法的计算结果见表 5,其中极限平衡方法计算结果是在滑动面确定的情况下算出的。图 2 几何模型Fig12 Geometry model图 3 有限元模型Fig13 FEM model表 4 计算采用的物理力学参数Table 4 Material properties材料名称容重/kN#m-3弹性模量/Pa泊松比内聚力/MPa内摩擦角/(b)岩体251E+1001211038结构面171E+7013011224通过有限元强度折减计算,当有限元计算不收敛时,程序自动找出了滑动面,如图 4。在一组平行的结构面中,只出现了一条滑动面,其余结构面59第 10 期郑颖人等:

16、用有限元强度折减法进行边坡稳定分析 没有出现塑性区和滑动。表 5 计算结果T able 5 Safety factor by different method计算方法安全系数有限元法(外接圆屈服准则)1126有限元法(等面积圆屈服准则)1103极限平衡方法(解析解)1106极限平衡方法(Spencer)1106图 4 坡体达到极限状态时形成的滑动面Fig14 The failure surface at limited state512 具有两组节理面的岩质边坡算例如图 5 所示,两组方向不同的节理,贯通率100%,第一组软弱结构面倾角 30b,平均间距 10m;第二组软弱结构面倾角 75b,

17、平均间距 10 m,岩体以及结构面计算物理力学参数见表 6。图 5 几何模型Fig15 Geometry model 按照二维平面应变问题建立有限元模型,在单元划分的过程中,在两个滑动面的交汇处形成了尖角,在建模时应在曲率突然变化的区域使用更细的网格以保证计算的准确性,或者将尖角抹圆。计算步骤同上,通过有限元强度折减,求得的滑动面如图 6a 所示,它是最先贯通的塑性区。塑性区贯通并不等于破坏,当塑性区贯通后塑性发展到一定程度,岩体发生整体破坏,同时出现第二、三条贯通的塑性区,如图 6b。程序还可以动画模拟边坡失去稳定的过程,从动画演示过程可以看出边坡的破坏过程也就是塑性区逐渐发展、最后贯通形成

18、整体破坏的过程。求得的稳定安全系数见表 7。其中,极限平衡方法计算结果是根据最先贯通的那一条滑动面求得的。表 6 物理力学参数计算取值Table 6 Material properties材料名称容重/kN#m-3弹性模量/Pa泊松比内聚力/MPa内摩擦角/(b)岩体251E+1001211038第一组节理171E+7013011224第二组节理171E+7013011224图 6 极限状态后产生的滑动面和塑性区Fig16 Failure surface and plastic zone表 7 计算结果Table 7 Safety factor by different method计算方法安

19、全系数有限元法(外接圆屈服准则)1162有限元法(等面积圆屈服准则)1133极限平衡方法(Spencer)113660中国工程科学第 4 卷6 结论 1)有限元强度折减法不需要对滑动面形状和位置做假定,通过强度折减使边坡达到不稳定状态时,非线性有限元静力计算将不收敛,此时的折减系数就是稳定安全系数,同时可得到边坡破坏时的滑动面。本文对其计算精度进行了分析,算例表明采用徐干成、郑颖人(1990)提出的摩尔-库仑等面积圆屈服准则求得的稳定安全系数与简化 Bishop法的误差为3%8%,与 Spencer 法的误差为 1%4%,证实了其实用于工程的可行性。2)目前对复杂节理岩质边坡的稳定分析尚没有好

20、的办法,传统的极限平衡方法无法得到岩质边坡的滑动面及其稳定安全系数,而各种数值分析方法只能算出应力、位移、塑性区等,无法判断边坡的稳定安全系数以及相应的滑移面。本文利用非线性有限元强度折减系数法由程序自动求得边坡的危险滑动面以及相应的稳定安全系数。通过算例分析表明了此法的可行性,为岩质边坡稳定分析开辟了新的途径。3)屈服条件的选择、计算模型的建立,包括计算范围、边界条件、网格划分密度和收敛标准等,应满足有限元计算的精度要求。如果使有限元计算保持足够的计算精度,那么有限元法较传统的方法具有如下优点:a.能够对具有复杂地貌、地质的边坡进行计算;b.考虑了土体的非线性弹塑性本构关系,以及变形对应力的

21、影响;c.能够模拟土坡的失稳过程及其滑移面形状。如图 7、图 8,可见滑移面大致在水平位移突变的地方,也是在塑性区塑性发展最充分的地方,呈条带状;d.能够模拟土体与支护的共同作用,图 9 为无锚杆(锚杆单元被杀死)时边坡稳定安全系数为111,图 10 为有锚杆支护时安全系数为 115,且塑性区后移。e1 求解安全系数时,可以不需要假定滑移面的形状,也无需进行条分。图 7 水平方向位移等值云图Fig17 Continuous contours of X displacement图 8 均质土坡达到极限状态时形成的滑动面Fig18 The failure surface at limitedsta

22、te of soil slope图 9 不加锚杆 X=111时的塑性区Fig19 The plastic strain zone without support图 10 加锚杆 X=115 时的塑性区Fig110 T he plastic strain zone with support参考文献 1 Griffiths D V,lane P A.Slope stability analysis byfinite J.elements.Geotechnique,1999,49(3):387 403 2 徐干成,郑颖人.岩土工程中屈服准则应用的研究J.岩土工程学报,1990,2:93 99(下转第

23、 78 页)61第 10 期郑颖人等:用有限元强度折减法进行边坡稳定分析 Chinese Science Bulletin,1999,44(17):1538 1548 15 杨春燕,张拥军,蔡文.可拓集合及其应用研究 J.数学的实践与认识,2002,32(2):301 208 16 合德利,原秀卿.整合策划技术 M.北京:时事出版社,2001 17 Philip H M,Mitchell Lee Marks.Managing themerger:Making it work M.New York:Simon&Schuster Inc,199218 杨春燕,张拥军.可拓市场的类型与实现方式研究J

24、.工业工程,2002,5(3):46 4919 李颖生,陈海滨.市场开拓M.北京:企业管理出版社,199920 Robert Heath 著.危机管理M.王 成,等译.北京:中信出版社,200021 陈 迅,王澍文.危机决策M.兰州:甘箫文化出版社,2002Study on Extension Strategic Tactical PlanningYang Chunyan1,Zhang Yongjun2(1.Research Institute of Extension Engineering,Guangdong University of Technology,Guangzhou 51008

25、0,China;2.Springf ield Wire(China)Inc,Dongguan,Guangdong 523910,China)Abstract Extension Strategic Tactics Planning(ESTP),based upon the science of extenics,is establishedalong a perspective of solving contradictory problems.In this paper,the fundamental thoughts,theories andmethods of EST P are int

26、roduced concisely.The highlights of the paper are that the techniques,used in case ofcontradiction to tackle such problems as how to generate the ideas of strategic tactical planning,how tocompletely analyze the resources,how to exploit extension resources,how to transfer the incompatible into theco

27、mpatible,and how to transfer the opposite into the co-existent,are developed in a formalized way.Andregarding ESTP.s practice,the applications of ESTP in resource integration,project,product innovationmarket development,and crisis prevention and handling are introduced.Key words extension strategic

28、tactical planning;contradictory problem;extension method;idea generation;extension resource(cont.from p.61)Slope Stability Analysis by Strength Reduction FEMZheng Yingren,Zhao Shangyi(Logistical Engineering University,Chongqing 400041,China)Abstract An analysis method for slope safety factor through

29、 c-Ureduction algorithm by finite elements ispresented.When the system reaches instability,the numerical non-convergence occurs simultaneously.Thesafety factor is then obtained by c-Ureduction algorithm.The same time,the critical failure surface is foundautomatically.T he traditional limit equilibri

30、um method can.t get the safety factor and failure surface of jointedrock slope.Strength Reduction FEM(SRFEM)presents a powerful alternative approach for slope stabilityanalysis,especially to jointed rock slope.This paper analyzes the precision and error caused by different soil yieldcriterions、FEM i

31、tself、slope height、slope angle、cohesion and friction angle thoroughly.Through a series of casestudies,the results show average error of safety factor between Strength Reduction FEM and traditional limitequilibrium method(Bishop simplified method)is 3%8%,the error between SRFEM and Spencer.s methodis 1%4%.T he applicability of the proposed method was clearly exhibited.Key words slope stability analysis;strength reduction by FEM;yield criterions78中国工程科学第 4 卷