降雨条件下非饱和土坡稳定性数值模拟

1、降雨条件下非饱和土坡稳定性数值模拟论文摘要：为了研究非饱和土坡在降雨入渗条件下的稳定性，运用饱和-非饱和渗流有限元法模拟降雨条件下饱和-非饱和土坡暂态渗流场的变化情况，分析了降雨强度、降雨持时以及土壤饱和渗透系数等参数对非饱和土坡稳定性的影响。分析结果说明：土坡的平安系数随降雨强度的增大而增大，随饱和渗透系数的增大而减小。由于降雨的进行，雨水入渗量逐渐增加，基质吸力逐渐丧失，孔隙水压力逐渐增大，因而土坡的平安系数随降雨持时的增加而减小。论文关键词：降雨,土坡稳定性,渗流,平安系数引言诱发边坡失稳的因素有很多，主要是降雨、地震、人类工程活动等。边坡中的渗透水压力和渗流场分布是影响边坡稳定和变形的

2、重要因素之一，而降雨，特别是历时长、强度大的暴雨更是导致边坡失稳破坏的常遇诱导因素。以往的边坡渗流场研究通常采用稳定渗流模型，对降雨入渗补给的作用，仅考虑多年平均降雨量对应的入渗条件，且入渗边界假设在地下水面上，忽略降雨在非饱和区的运动过程。这种假设对于那些地下水位较高，非饱和带较浅的地区影响并不大；但对于地下水位较低，具有深厚非饱和区的边坡，那么显得不尽合理。传统的边坡稳定分析方法也是建立在饱和土假设根底之上的，理论上忽略基质吸力对边坡稳定的奉献。这样的考虑是偏于平安，同时，也是由于非饱和土的研究还不够深入造成的。近年来，随着非饱和土力学的开展，为更加合理的进行边坡稳定分析提供了新的理论根底

3、。因此，本文考虑建立降雨入渗条件下的饱和-非饱和土坡计算模型，对雨水入渗引起的土坡平安系数下降进行定量分析，对于评估边坡的稳定性、预报滑坡灾害和边坡治理均具有重要的指导意义。1饱和-非饱和渗流的控制方程1.1饱和-非饱和渗流的控制方程及定解条件饱和-非饱和渗流问题所满足的微分方程和边界条件为：(1)式中：S水头分布规律的边界，H为边界水头，称为第一类边界条件；S流量情况己知的边界，q为单位时间边界法向流量，为边界面外法线方向余弦，称为第二类边界条件；S渗出面边界，z为渗出面节点坐标，可以归为第一类边界条件；H初始渗流场水头，称为初始条件。1.2Galerkin法有限元列式将渗流场离散为有限个单

4、元。假设单元内水头与结点水头的关系为：式中，N为单元形函数；n为单元节点数。将式将上式以矩阵形式表示并积分得：式中，为总渗透矩阵，由单元渗透矩阵叠加而成，=dA；、为形函数对坐标的微分向量，；为单元渗透矩阵；h为节点总水头列向量，h=h,h，,h；为总贮水矩阵，由各单元贮水矩阵叠加而成，=dA，；F为反映边界条件的流量矢量，F=NqdS。对于时间的离散可以采用有限差分法。采用中心差分为：采用向后差分为：一般来说，中心差分法较向后差分法得到的结果精确；而向后差分法那么在减小数值波动方面更为有效。1.3有限元计算中几个问题的处理在非稳定渗流计算中，时间步的取值不能过大，否那么将不收敛或收敛太慢；时

5、间步也不能太小，否那么将产生过大的累积误差。在实际计算中，应先进行一定量的试算，根据计算的收敛情况选定一初始时间步进行计算。如果收敛，那么取此值为该时段的步长，否那么，时间步减半，直到获取收敛的时间步为止。边界条件确实定是解决问题的关键之一。本文中采用的典型土坡剖面如图1所示，非饱和区边界AH、DE的状态很难确定。事实上，它是由上游和下游区域渗流情况确定的，具有不确定性，因此需要一定的假设，本文假设该区域不因降雨而和相邻区域发生水交换，即为不透水边界。根据雨水入渗的实际特点，入渗边界可以考虑为两种模型：1积水模型。当降雨强度小于土壤的入渗能力时，雨水将全部入渗，按第二类边界条件处理；2降雨模型

6、。当降雨强度超过土壤的入渗能力时，地表产生积水或径流，按第一类边界条件处理。本文根据初始条件和入渗过程的不同，将入渗边界AB、BC、CD按积水模型;和降雨模型;相互转化考虑，即按式当q式中，r 降雨持时 t/h 饱和渗透系数 k / No.1 60 8 10 No.2 60 8 10 No.3 60 8 10 No.4 30 16 10 No.5 30 16 10 No.6 30 16 10 No.7 15 32 10 No.8 15 32 10 No.9 15 32 10 3土坡稳定结果分析3.1降雨强度与饱和渗透系数的影响降雨强度与土壤饱和渗透系数是相互影响、相互制约的，需要将它们结合起来

7、考虑。在进行土坡稳定分析时，同样需要将它们结合起来考虑。图2土坡平安系数与降雨强度和土壤饱和渗透系数的关系Fig.2Therelationofslopesafetyfactorandrainfallintensityandthepermeabilitycoefficientofunsaturatedsoil土坡平安系数与降雨强度和土壤饱和渗透系数的关系如图2所示。总体来说，土坡的平安系数随降雨强度的增大而增大，随饱和渗透系数的增大而减小。土壤饱和渗透系数为k=10m/s时，当降雨强度逐渐增大时，雨后土坡孔隙水压力的变化量逐渐减小，相应的抗剪强度逐渐增大，因此，土坡的平安系数逐渐增大；土壤饱和渗

8、透系数k=10m/s时，三种降雨强度下的孔隙水压力变化相同，相应的抗剪强度变化相同，那么土坡平安系数相同。降雨强度为r=60mm/h时，当土壤饱和渗透系数逐渐增大时，雨后土坡孔隙水压力的变化量逐渐增大，相应的抗剪强度逐渐减小，因此，土坡的平安系数逐渐减小。当土壤饱和渗透系数为10m/s时，三种降雨强度下，雨后土坡始终是不稳定的，说明较大的土壤饱和渗透系数对非饱和土坡稳定的影响很大。3.2降雨持时的影响土坡平安系数与降雨持时的关系如图3。总体来说，土坡的平安系数随降雨持时的增加而减小。这是由于随着降雨的进行，雨水入渗量逐渐增加，基质吸力逐渐丧失，孔隙水压力逐渐增大造成的。在降雨初期，平安系数下降

9、较快，特别是土壤饱和渗透系数较大的土坡，这主要是由于基质吸力对土壤含水率比拟敏感，这从土-水特性曲线见表1可以看出。如图3a，在降雨强度为60mm/h，土壤饱和渗透系数为k=10m/s时，4小时以后，土坡将不在稳定；如图3c，在降雨强度为15mm/h，土壤饱和渗透系数为k=10m/s时，2小时以后，土坡已经不稳定。随着降雨的进行，土坡平安系数的下降减缓，并逐渐趋于一稳定值，这是由于此时孔隙水压力变化减慢造成的。图3a降雨强度为60mm/h图3b降雨强度为30mm/h图3c降雨强度为15mm/h图3土坡平安系数与降雨持时的关系Fig.3Therelationofslopesafetyfactor

10、andtherelationshipwithrain4结论降雨入渗是边坡失稳破坏的常遇诱导因素。究其原因，主要是因为这些土坡的近地表浅层土大多为非饱和土，在大气降雨时，雨水人渗使土体饱和度增加，基质吸力锐减并引起土壤抗剪强度大幅度下降，所以当持续降雨的历时和强度超过某种限度时，雨水入渗到达一定深度，那么可导致土坡失稳。总之，土坡的平安系数随降雨强度的增大而增大，随饱和渗透系数的增大而减小；土坡的平安系数随降雨持时的增加而减小。参考文献3 Fredlund D G,Rahardio H.Hillside slope stability assessment in unsaturated residual soils.IKRANM Swminar on the4 SpinerenburgSEJ,van Esch JM.Slope stability during infiltration.Proc-eedings of Inte