我国是一个地质灾害较多的国家,其中滑坡灾害占有较大比例。尤其是随着我国水利建设的大力兴建,时常有库岸滑坡发生失稳,对社会和人民造成严重威胁。水的渗透会严重影响滑坡的稳定性。实际上,只要有水存在的地方,应立场和渗硫场就会相互影响、相互作用,处于一种复杂的动态变化过程中,构成流固耦合关系。 本论文通过对耦合作用进行实质性分析,得出应力场合渗透场相互作用的机理:一方面,应立场的改变会影响土体渗透系数,从而影响渗流场;另一方面,渗流场通过改变水荷载来影响应立场。根据岩土体骨架的平衡方程和虚功原理建立应力场平衡方程。根据地下水的连续方程和变分原理建立渗流场微分方程。 应力场平衡方程和渗流场微分方程共同组成流固耦合问题的控制方程组。该方程组通常为非线性方程组,论文借助国际大型有限元软件ABAQUS进行了滑坡体内流固耦合的计算,继而对滑坡稳定性作了定量评价。 滑体本构关系具有高度的非线性,数值计算中可采用Duncan-Chang本构模型来模拟。本文以ABAQUS有限元软件为开发平台,提出了Duncan-Chang本构模型的开发思路和实现步骤,完成了模块的编程,并且依据常规三轴试验数据对开发的程序进行了一系列严格验证。结果证明开发是成功的。 茅坪滑坡是清江隔河岩水库库区一个十分典型和重要的大型古滑坡体,具有多级、多期次的滑动特征。该滑坡一旦失稳,必将产生严重后果。论文对茅坪滑坡体地质特征及监测资料进行研究,对滑坡体的稳定性作了定性分析。另外,通过对滑坡的数值分析计算,对滑坡稳定可靠度作了定量分析。结果表明:在天然水位下,滑坡中的塑性区范围较有限,尚难危及滑坡的整体稳定,但滑坡后缘附近水力梯度较大,渗流力较大,使该位置局部出现了拉应力,导致滑坡后缘稳定性较差。当水位急剧上升时,在斜坡前缘和后缘都出现了塑性区,滑坡整体稳定性降低。尤其是滑坡前缘稳定性变差,有局部失稳的可能性。 本文的工作对滑坡稳定性分析、流固耦合分析具有积极作用,而且本文提出的模拟方法、求解技术和开发思路对水工结构和岩土工程及其他工程领域的数值仿真研究也具有广泛的借鉴意义。