高等土力学课程论文

昆明理工大学

高等土力学课程论文

降雨入渗过程中的土质边坡稳定分析

姓名：***

年纪：2012

学号：**********

专业：建筑与土木工程

20130428

1摘要

边坡的稳定性由内在因素和外在因素共同决定。内在因素是边坡本身所固有的，外界因素主要包括降雨、地震、开挖等。近年来，降雨导致边坡失稳的情况越来越多。因此，研究分析降雨对边坡稳定性的影响以及准确预报事故发生的可能性，是当前亟待解决的一个复杂的工程问题。常规的边坡稳定性分析方法主要基于饱和土理论和稳定渗流场情况，难以反映降雨作用对边坡稳定的影响。

本文针对受降雨影响的边坡，运用饱和一非饱和土壤水分运动的理论和二维非稳态渗流有限元模型，模拟雨水入渗引起的暂态渗流场，将计算所得到的暂态孔隙水压力和渗流力分布用于考虑基质吸力影响的边坡稳定安全系数的计算中，建立了考虑非饱和土边坡从雨水入渗到出现滑坡危险全过程的计算模型；通过实例分析，研究了降雨对非饱和边坡渗流场分布、发展和对边坡稳定性的影响、状态预测等进行了较为深入的研究，得出了一些有意义的结论。

结合不同的降雨条件，计算不同降雨持时情况对土质边坡稳定的影响，得出随着降雨持时的增加，边坡稳定性逐渐降低，降低幅度随降雨持时增加而增加的结论。

本文的研究工作在改进和完善边坡稳定性数值分析方法方面进行了有益尝试。通过

对数值算例较为深入细致的分析，对降雨引发非饱和土边坡失稳的机理和规律有了进一步的认识，可为非饱和土边坡稳定分析、滑坡灾害预报和边坡治理等方面提供有用的参考数据，具有一定的理论意义和实用价值。

关键词：降雨入渗；边坡稳定；非饱和土；饱和一非饱和渗流；有限元分析

2论文研究背景

降雨和滑坡都是人类生活中常见的自然现象，在自然规律作用下，两者有着紧密的关系。滑坡是土木、水利、交通、矿山等基本建设工程常见的事故和灾害，触发滑坡的因素多种多样，水是诱发滑坡的主要因素，而降雨，特别是历时长、强度大的暴雨更是导致边坡失稳破坏的常遇诱导因素。在我国特别是南方地区，尤其以降雨入渗的影响最为频繁和明显。

降雨型滑坡产生的机理主要表现在两个方面：一是降雨动能外部作用，二是降雨入渗后坡体内部产生的相关作用。对于由降雨因素导致的边坡失稳机制，目前，较为普遍的认识是：降雨入渗使得边坡体内的地下水潜水面升高，滑面处土体软化，从而降低边坡的稳定性，导致滑坡的发生。基于这一认识的评价、分析方法是，运用饱和渗流理论，模拟降雨引起的潜水面升高；运用刚体极限平衡方法，计算潜水面升高后的边坡稳定性，并以此作为设计依据。这是一种基于饱和土力学理论的边坡稳定性分析方法。

但是这种评价方法往往很难被实践经验丰富的工程师所采信。主要是因为：①降雨对边坡稳定性的影响与具体地区有关，不同地区诱发滑坡的临界降雨量有一定的差别。例如，四川盆地滑坡主要发生在日降雨强度超过200mm／d的分布区，香港滑坡发生的日降雨强度大都在100mm／d以上；②即使对于同一地区的边坡，不同专家推荐的确

定临界降雨量的标准也有所不同。例如，有建议以前期总降雨量值来确定临界降雨量的；有建议以日降雨量值来确定临界降雨量的；也有建议小时最大降雨量值来确定临界降雨量，等等。事实上，边坡失稳与总降雨量的大小、日降雨强度以及降雨持续时间的长短等均有着直接关系。

故已有的工程设计均采用对孔隙水压力进行假定的方法。例如美国一些工程采

用水面达地表的静水压力分布。这一假定基于以下认识：历时长的降雨使边坡裂隙完全饱水，地下水位达到地表。对于高边坡工程，采用这一水压力分布使加固设计过于保守。

我国一些边坡工程常将静水压力乘以折减系数来作边坡设计，例如漫湾折减系数0．4；

五强溪水电站则采用类似Hock建议图形，但取值较小；三峡工程采用折减系数O．3。目前采用强排水体系而不考虑暂态水压力。图1．1是不同方法采用的边坡暂态水压力分布，彼此差别极大，且均缺少理论分析及实测数据的支持。

降雨过程中边坡的稳定分析是岩土工程实践中最受关注的一个课题之一。对于一个

大型工程边坡的设计，若能科学合理地确定边坡中的孔隙水压力，首先，将节省大量的

工程投资；其次，工程的安全稳定运营将带来巨大的社会效益以及经济效益。同时，这

一方面的研究，将使我们对于由降雨因素导致边坡失稳的作用机制有一个更加深入的认识。

3降雨入渗对非饱和边坡渗流场变化的影响及算例分析

3.1均质各向同性边坡渗流场算例分析

某边坡坡顶高程为15．00m，地面高程为0．00m，边坡坡率为砚=3：4，边坡上游水

位为10．00m，下游水位为1．00m，取边坡计算长度80．00m，假设模型的下边界为不透水层。其剖面如图1所示。

假定边坡是均质各向同性的，饱和渗透系数k。=3．7x10～cm／s。

利用大型有限元软件ABAQUS进行模拟，根据己建立的模型，采用CPE4P四边形

单元，进行结构化网格划分，整个模型共划分2436个节点，2523个单元，如图2所示。

图1 边坡的断面图图2边坡有限元计算网络计算所得边坡断面的孔隙水压力等值线和饱和度等值线分布情况如图3、图4 所示。

图3边坡孔隙水压力等值线图图4边坡饱和度等值线图

利用ANSYS热分析功能模拟二维渗流渗流进行复核。计算条件与ab边，总水头H=30m，efgh边，总水头H=21m；bcd边和ha边为第二类边界，流量q=0；de 边假定为溢出边，总水头H=Y。计算所得断面的孔隙水压力等值线和总水头等值线分

布情况如图5、图6所示。

图5边坡孔隙水压力等值线图图6边坡总水头等值线图

ABAQUS的计算结果与ANSYS热分析模拟渗流计算所得的比较吻合，边坡孔隙水压力的大小以及分布都非常的相近，饱和区和非饱和区的大小、位置以及逸出点的位置也相同。

可以判断出，对于均质各向同性边坡渗流的情况，ABAQUS计算所得的孔隙水压力、饱和度分布较为合理，计算结果能真实的反应实际自然状态的土质边坡渗流场的状态。

3．3降雨入渗非饱和边坡的渗流分析

以一均质各向同性的边坡为研究对象，图7给出了均质土体断面示意图，设地表的高程为Om，坡顶高程15m，边坡坡比3：4，地下水位高程一4m，取计算底面高程一20m。

图8给出了有限元网格示意图，整个计算过程中地下水位不变。

图7计算边坡示意图图8有限元网格图

土体为各向同性的均质细砂土，水的体积模量为2e9。土体介质的饱和渗透系数t=3.7e-3cm／s，渗透系数随饱和度变化曲线如图9所示，土一水特征曲线选取如图10所示。

图 9渗透系数随饱和度变化曲线图10土一水特征曲线