应用有限元法求解有渗流的边坡稳定性问题-江西理工大学学报.pdf

第 卷第 期 !( % 南方冶金学院学报 .DE’!(*8D’% 年 月 ! % 45678/9 5: ;56-=78 8;--6-= 5: ?=-/9967@A 4BC ’ *! 应用有限元法求解有渗流的边坡稳定性问题 焦永斌 吴祥松 刘小生 * * +南方冶金学院环境与建筑工程学院，江西 赣州 ( %, 摘 要：介绍了有限元方法在解决有渗流等外力作用下的边坡稳定性问题中的方法 和步骤给出了渗流方程有限元表达式，并结合岳城水库土坝发生滑坡的计算实例 描 * * 述了计算过程和结果’ 关键词：有限元；边坡稳定；渗流计算；自由面 中图分类号： 文献标识码： 文章编号： （ ） -.!!’( / %0#%!!$ ! %#($#1 ! 引 言 以往求解渗流问题的解析法只用于含水层几何尺寸简单且是均质、各向同性的情况，而实 际的水文地质条件往往比较复杂 如果遇到含水层是非均质、各向异性、厚度变化且边界条件 ’ 复杂、形状不规则等情况，用解析法求解就必须首先对实际情况进行简化，建立能利用解析法 求解的数学模型但是，如果把数值法和计算机结合起来，就能解决许多复杂的渗流问题 目前 ’ ’ 2%#3 国内外有一些学者进行了这方面的研究，但是文献较少 本文将从具体实例出发说明如何应 ’ 用有限元法求解有渗流的边坡稳定性问题’ 利用有限元法进行渗流分析的基本思路是将渗流微分方程的定解问题 （包括边值和初值） 化为变分问题，即泛函求极值问题，然后把求解区域分割成许多小单元，应用里兹法的概念，通 过区域剖分和分片插值，把求泛函极值问题化为代数方程组求解，代数方程组的解即为所求渗 流问题的解’再根据渗流作用概念，将作用在滑动面上和划分土块的表面水压力转换为等价的 体积力从而一次完成在渗流和地震等各种外力作用下的稳定分析计算 * ’ 有限元法计算步骤#$%’ 利用有限元法求解渗流问题的关键是建立其单元特性方程，无论通过哪一种途径，基本上 都有相同的步骤’ 以平面恒定渗流问题为例，求解渗流问题的计算步骤如下’ 定解区域的离散化 把连续的定解区域划分为有限个小单元，渗流流态较平稳的区域 +%, ’ 单元可划分得稀一些’单元的数目由计算精度的要求和计算机容量决定’ 插值函数的选定 根据里兹法的概念把水头函数 近似的表示为含有结点坐标的 +!, ’ ! +#$, 线性插值函数对于三角形单元，常采用如下线性多项式来近似表达单元域的待求函数： ’ 收稿日期： !!#$#% 基金项目：江西省自然科学基金资助项目 （ ） %%% %% ’ 作者简介