有限元强度折减法和极限平衡法适用性研究

《有限元强度折减法和极限平衡法适用性研究》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、有限元强度折减法和极限平衡法适用性研究曾红丽戚明军陕西省建筑设计研究院有限责任公司710003西安摘要：简述了有限元强度折减法原理和两种极限平衡理论进行边坡稳定性分析的方法，分别利用ANSYS和SLOPE/W程序对•一均质土坡进行了稳定性分析。建立了同样尺寸的两种il•©模型，得岀了分别采用有限元强度折减法、Bishop法和Janbu法在未考虑地下水的情况下边坡安全系数。将计算结果进行比较分析，表明基于极限平衡理论的分析方法分析得到的结果偏小，而有限元强度折减法巾于考虑了土体内部应力应变关系所得结果更切合工程实际。关键词：有限元强度折减法；极限平衡法；适用性；稳定

2、性1引言边坡工程是公路、铁路、水利水电以及矿山工程一个不可或缺的组成部分，因此边坡稳定性问题的研究就成为岩土工程界研究的热点问题。国内外学者，已经取得了大量研究成果。0前，研究边坡稳定性的传统方法主要有：极限平衡法，极限分析法，滑移线场法等。随着计算机技术的发展，有限元强度折减法在边坡稳定性分析中己经得到了较好的应用"~51。本文通过算例分析，比较分析有限元强度折减法和基于极限平衡理论的方法在边坡工程屮适用性的差别。2有限元强度折减法原理［6j所谓强度折减，就是在理想弹塑性有限元计算中将边坡岩土体抗剪强度参数逐渐降低直到其达到破坏状态力止，程序可以自动根据弹塑性计

3、算结果得到破坏滑动面，同时得到边坡的强度储备安全系数K，也称为强度折减系数。通常，边坡的稳定性安全系数定义为沿滑动而的抗剪强度与滑动而的实际剪力的比位，公式表示为K=f(ctan(p)dA/ftcIA(1)将式(1)两边同时除以K，得(2)式中，c*=c/K,(p'=arctan(tan(pjK),为滑动面上的抗剪强度。通过逐步调整系数K，得到不同的c1,#为粘聚力，p为内摩擦角，将#代入有限元程序,cr为滑动面上的法向应力，r反复分析边坡，直到坡体达到临界状态，坡体达到临界状态时的K值即作为边坡稳定性安全系数。3极限平衡法利用极限平衡法对边坡的稳定性进行了分析，

4、该方法的基本特点是，只考虑静力平衡条件和土的Mohr-Coulomb破坏准则，也就是说，通过分析土体的破坏那一刻的平衡来求得问题的解［7］。极限平衡理论的主要思想是将滑动土体进行条分，由极限状态下土条所受力和力矩的平衡来分析边坡稳定性。它是目前应用最多的一种分析方法。文中所用极限平衡理论的方法有以下儿种(1)Bishop法：该方法考虑了土条间的作用力，这是对传统的瑞典条分法的重要改进。该方法忽略各土条之间的切向条间力，认为条间力的合力是水平的，同时假设破坏面是圆弧面，且定义边坡安全系数为沿整个滑动面上的抗剪强度与实际产生的剪应力的比值。该方法的计算结果比较接近实际

5、，常用于土质、软岩质及碎岩边坡的稳定性分析。(2)Janbu法：Janbu提出了同时满足力和力矩平衡的“通用条分法”。这一方法区别于其他方法的一个重要方面，就是通过假定土条侧向力的作用点而不是作用方向来求解安全系数的。采用非圆弧曲*，按条块滑动平衡确定条间力，按推力线确定法向力的作川点，简化计算条间切向力为零，然后再对稳定系数进行修正。4算例分析4.1算例M某均质土边坡工程，坡角接近45%坡高20m。土体物理力学参数见表1。表1土体物理力学参数表粘聚力c/kPa内摩擦角小/°剪胀角弹性模SE/MPa重度P/kN•m-3泊松比P17.030.0050.017.00.

6、354.2有限元强度折减法求解［1()］计算边界范围的人小在有限元法中对计算结果的影响比在传统极限平衡法中表现得更为敏感。如图12示，建立模型吋坡角到左端边界的距离为坡高的1.5倍(30m),坡顶到右端边界的距离力坡高的2.5倍(50m),上下边界总高为坡高的2倍(40m),这时计算精度较为理想有限元网格剖分采用四边形八节点等参单元，将边坡中鞞易发生滑动的区域单元进一步细化。边界条件：左右两侧水平约束，下部X、Y方向均约束，上部边界为自由边界。图2有限单元网格采用理想弹塑性模型和Drucker-Prager屈服准则］1。Drucker-Prager屈服准则是-•

7、种经过修正的Mises屈服准则，他考虑丫平均应力^或b对屈服强度的影响，对土体材料有较好的适用性。该准则形式如下:cxl}+=k(3)在平而应变状态下，式屮.•sm(p5/373+sin(pa/3ccos^9a/3+sin(p人为应力张量第一不变量，A为应力偏量第二不变量。考虑到土体的材料非线性，有限元分析屮采用非关联流动法则，应用Newton—Kaphson迭代法对边坡进行大位移大变形分析。通过调整折减系数后反复计算，根据传统的计算不收敛即为破坏的判断标准，当K=1.54吋边坡处于极限平衡状态。得到当折减系数为1.54时，迭代过程收敛，但是当折减系数为1.55时

8、，迭代过程