降雨入渗条件下边坡的稳定性分析.pdfVIP

水 利 学 报 2008年9月 SHUIUXUEBAO 第39卷第9期 文章编号：0559—9350(2008)09—1066．08 降雨入渗条件下边坡的稳定性分析 周家文1，徐卫亚1，邓俊哗1，刘兴宁2 (1．河海大学土木工程学院，江苏南京210098；2．中水顾问集团昆明勘测设计研究院，云南昆明650051) 摘要：把降雨入渗条件下的边坡稳定性问题转化为图论中的最短路问题，基于饱和一非饱和渗流有限元计算得到 的应力场和渗流场来计算边坡的最危险滑动面及其安全系数，在降雨入渗过程中，考虑基质吸力对滑动面抗剪强 度的影响，采用修改的Mohr-Coulomb破坏准则来计算非饱和边坡的抗剪强度。降雨入渗过程及其对边坡力学行为 的影响采用一个二维的有限元渗流一应力耦合程序进行分析，研究降雨入渗过程对边坡的孔隙水压力、位移以及 安全系数的影响，计算结果表明，在降雨入渗过程中，边坡的安全系数在不断减小，降雨停止之后，边坡的安全系数 还会持续减小一段时间，随着排水的不断进行，边坡的安全系数又会逐渐增大。 关键词：边坡稳定性；降雨入渗；饱和一非饱和渗流；极限平衡有限元方法；图论；排水 中图分类号：TU457 文献标识码：A 1研究背景 边坡稳定性问题一直是岩土工程的一个重要课题。统计资料表明，边坡失稳破坏的最主要诱发因 素是降雨入渗。在日本，平均每年由于降雨导致失稳破坏的边坡超过l万以上甚至多达10万个…。在 美国，非饱和特殊土问题所造成的损失超过风灾、洪灾、地震、泥石流等自然灾害损失的总和，年损失近 100亿美元性1。在我国，由于降雨入渗引起的边坡失稳造成的损失也非常巨大。降雨入渗条件下的边 坡稳定性问题受到了许多学者的关注。非饱和边坡降雨人渗过程中，降雨人渗导致边坡地下水位线上 升和边坡浅层含水率上升，由此导致孔隙水压力的增大或者说基质吸力的减小，边坡潜在危险滑动面的 抗剪强度在不断减小，直到不能满足抗剪的要求为止，边坡从而产生破坏。降雨入渗条件下边坡的稳定 性不仅仅受到初始边界条件、降雨强度和渗透系数的影响，还取决于降雨时间、材料参数和初始含水率 等。由于基质吸力(It,。一fw)的存在和大小对边坡的稳定性致关重要，所以基质吸力的现场量测对于研 究降雨型边坡非常重要，文献[3]通过对现场得到的边坡基质吸力进行分析，得到了几条有益的结论。 为了研究降雨入渗条件下边坡不同高程位移的变化情况，边坡位移的现场量测对于研究降雨型边坡也 非常重要，文献[4]依据试验数据对非饱和土中基质吸力变化产生的应变进行了研究，对降雨引起的边 坡位移机制和规律进行了探讨。文献[5]采用强度折减有限元法对渗流条件下边坡稳定性进行了计算 入渗条件下边坡的稳定性问题，采用比较多的是基于饱和一非饱和渗流计算得到的应力结果来计算边 坡的最危险滑动面及其安全系数’8’，并且对降雨入渗条件下边坡稳定性的影响参数进行敏感性分 析‘引。 本文基于饱和一非饱和渗流计算过程中得到的节点应力数据和孔隙水压力数据，通过引入图论方 法，把降雨人渗条件下的边坡稳定性问题转化为图论中的最短路问题，采用Fredlund¨21提出的降雨入渗 收稿13期：200／-05．17 基金项目：国家自然科学基金重点项目 ·_——1066--—- 万方数据 条件下的抗剪强度准则计算滑动面的剪应力r和抗剪强度o，从而可以确定降雨入渗条件下边坡的最 危险滑动面及其安全系数。 2降雨人渗条件下的边坡稳定性分析 2．1饱和、非饱和渗流非饱和渗流与饱和渗流一样，也符合达西定律： ”以㈩蠡。器 (1) 式中：吼为水的达西流速；k。为饱和渗透张量；H为总水头。 水流运动在饱和区和非饱和区的主要区别在于渗透系数的不同，在饱和区，一般认为渗透系数为一 常数，而在非饱和区，渗透系数为饱和度、基质吸力的一个关系式。在非饱和区，孔隙水压力为一负值， 而在饱和区，孔隙水压力为一正值。即使有这些区别，但是局部流动微分方程还是一样的。根据质量守 恒原理，饱和一非饱和渗流控制微分方程为 毒【k