考虑渗流与蠕变耦合作用的水库滑坡变形数值分析

水库滑坡土体的变形主要是流固耦合作用的结果,但是变形并不是立刻发生的,大都经历了一个时间发展的过程,具有典型的流变特性。本文基于某滑坡土体的蠕变试验结果,推导了扩展的Burgers蠕变模型M–2K,并将其推广至3维形式,应用蠕变试验数据对推导的M–2K蠕变模型进行拟合,求得了蠕变模型参数。在传统的流固耦合分析中引入了时间变量关系,构建了渗流与蠕变耦合作用的数学模型,对其进行有限元求解,并应用MATLAB软件编写了有限元程序。以三峡库区白家包滑坡为例,应用该程序分别对滑坡体进行了考虑蠕变及不考虑蠕变的流固耦合计算。研究结果表明:3维形式的M–2K蠕变模型拟合程度较高,能反映土体的蠕变性质;有限元计算结果说明考虑渗流与蠕变耦合作用下孔压减少值大于不考虑耦合作用下时孔压减小值;选取关键点的耦合及非耦合条件位移计算结果与实际的GPS位移进行比较,两种计算值与实际监测值变形趋势基本一致,考虑渗流与蠕变耦合条件下两个监测点的位移值比不考虑耦合时大,而且与实际值比较接近。研究结果表明:所建立的渗流与蠕变耦合数学模型是合理的,有限元数值求解程序是正确的,该成果为水库滑坡的预测与预报提供了新途径。