基于强度折减的离散单元法在岩质边坡稳定性分析中的应用.docxVIP

第25卷第5期长?江?科?学?院?院?报Vol.25No.5

2008年10月JournalofYangtzeRiverScientificResearchInstituteOct.2008

??文章编号:1001?5485(2008)05?0107?04

基于强度折减的离散单元法在岩质边坡

稳定性分析中的应用

周先齐1,2,徐卫亚1,钮新强2,郑文棠1

(1.河海大学土木工程学院,南京?210098;2.长江勘测规划设计研究院,武汉?430010)

摘要:离散元方法是一种用于非连续介质的数值计算方法,适用于岩体断层、节理、裂隙比较发育的岩质高边坡稳

定性分析。强度折减法是一种在边坡稳定性分析中应用比较广泛的数值方法,其所认定的强度储备系数跟极限平

衡所提出的稳定安全系数是等效的。将离散元方法与强度折减准则相结合,采用两套离散网格:密网格代表节理

裂隙发育的工况;稀网格近似等效为岩体经过加固后完整性得到提高的工况,对构皮滩水垫塘高边坡进行了计算,

分析了自然边坡及开挖边坡的稳定性,并对支护效果进行了评价,为工程的设计施工提供了参考。

关?键?词:离散元法;强度折减;非连续介质;安全系数;位移矢量差

中图分类号:TU457???文献标识码:A

??现代水利工程规模越来越大,产生的边坡越来

对边坡稳定性进行分析评价是一种较好的分析方

越高,其稳定性日益得到重视。特别是富含节理、裂法。可以借助离散元较有效地模拟节理岩体大变形

隙、断层的高边坡,其稳定性应该有一个更深入的评特性,同时由强度折减理论得到边坡的安全系数,以

价。对于结构面岩质边坡,传统的用于土质边坡的有效地指导工程施工。

滑弧面搜索方法已很难直接应用[1],极限平衡方法本文以强度折减和离散元为手段,用于构皮滩

则应用广泛,已被很多工程设计人员接受,但是对于水电站节理岩体边坡稳定性分析中去,对边坡安全

复杂的节理岩质边坡,极限平衡法也不能很好的评稳定性进行了评价。

估。借助数值分析方法评估边坡稳定性方法,在边

坡稳定性分析中应用很广。数值模拟中的连续介质

方法具有较好的计算精度,但多是基于小变形,这与

1?相关理论

实际相差较大。近年来迅速发展的非连续介质方法

1.1?离散元计算原理

是一种比较适合节理岩体的数值方法,其中离散元离散单元法[5]是由Cundall于l971年提出来的

就是一种应用比较广泛的非连续介质方法。刚体离一种岩石力学计算分析方法,其理论基础是牛顿第

散单元法可以考虑块体之间接触面的滑移、分离和

二运动定理,结合不同本构关系,采用动态或者静态

倾倒等大变形,可变形离散元法还可以计算岩体内松弛方法进行求解,其基本运动方程为

部的变形和应力分布。因此,对于节理岩体,离散元

具有很好的优势。强度折减法在边坡稳定性分析中

m

d2u

2u

2(t)+c

dt

du

dt

+ku(t)=f(t)。(1)

应用也非常广泛。强度折减法的思想就是将结构面式中:m是单元的质量;u是位移;t是时间;c是

粘性阻尼;k是刚度系数;f是单元所受的外荷载。

的强度参数按某一规则进行折减,直到边坡出现破

坏,这个时候的折减系数即为边坡的安全系数。用

该方法是用来解决不连续介质问题的一种数值

这个方法求出的安全系数基本上跟极限平衡法得出模拟方法,比较适用于节理发育的岩体。其把节理

的结果一致[2]。Duncan[3]将这种思想用于土坡的

岩体看作是离散的岩块和岩块间的节理面这样一个

稳定性分析之中。赵尚毅等[4]用强度折减非线性二元介质,各岩块之间的接触可以看成角-角接触、

有限元模型来模拟分析节理岩质边坡的稳定性并得角-边接触、边-边接触[6]。由于其离散性,块体

到了边坡的安全系数,为节理岩质边坡稳定性分析之间可以产生平移、转动和变形,节理面可被压缩、

开辟了新的途径。分离或者滑动。

由此可见,将强度折减法与其它方法进行结合,求解思路当块体所受的作用力系和边界约束条

收稿日期:2008?06?18

作者简介:周先齐(1981?),男,湖北麻城人,博士研究生,主要从事岩石力学与工程的研究,(电话)(电子信箱)xianqizhou@

163.com。

108长江科学院院报2008年

件发生变化时,