各向异性饱和土体的渗流耦合分析和数值模拟-中国科技论文在线.PDF

__________________________________________________________________________ 第21 卷 第10 期 岩石力学与工程学报 21(10) ：1447～1451 2002 年10 月 Chinese Journal of Rock Mechanics and Engineering Oct.，2002 各向异性饱和土体的渗流耦合分析和数值模拟* 杨林德 杨志锡 (同济大学地下建筑与工程系 上海 200092) (上海交通大学建筑工程学院 上海 200030) 摘要 将饱和土体视为均质、连续的各向异性弹塑性多孔介质，根据虚位移原理推导出饱和土体内各向异性渗流 直接耦合的有限元法计算公式。针对直接耦合法所生成的病态方程采用 MATLAB 语言编写出平面条件下的计算 程序，对各向异性弹性多孔介质中 Mandel 效应进行数值模拟分析。计算结果验证了本文直接耦合有限元法的正 确性和适用性。 关键词 耦合，渗流，各向异性，有限元法 分类号 O 357.3，TU 44，O 241.82 文献标识码 A 文章编号 1000-6915(2002)10-1447-05 1 引 言 2 各向异性耦合场的基本方程 对于饱和软土的渗流耦合分析，Sandhu 和 由于饱和土体内位移场和渗流场是两个具有不 Wilson(1969)假定孔隙水及土体颗粒不可压缩，根 同运动规律的物理力学环境，所以要描述其耦合响 [1] 据变分原理推出Biot 方程的有限元法计算公式 ； 应的数学模型也应包含位移场和渗流场的控制微分 Christian 和 Boehmer(1970)采用有限元法和有限差 方程及其对应的边界条件和初始条件。下面，分别 [2] 分法相结合的方法求解 Biot 方程 。国内沈珠江 介绍本文所用到的应力场的基本方程、渗流场的基 (1977)首先将Biot 理论的有限元法应用于土体固结 本方程、边界条件和初始条件。 [3] 分析 ；殷宗泽等(1978)根据流量平衡的概念，结合 2.1 应力场的基本方程 [4] 虚位移原理得到了类似的结果 ；龚晓南(1981)假 在直角坐标系中，以张量形式表示的饱和土体 设等价接点流量等于等价接点压缩量，推出 Biot 的平衡方程为 [5] 方程 ；张有天(1989)等对岩土渗流力学的研究也作 [6] s ij ，j Fi 0 (1) 出了重要贡献 。上述学者的工作均以Biot 理论为 基础，假定岩土介质为均质、连续、各向同性体。 同时，各向异性土体在外部条件作用下的力学 考虑到土体的地质沉积作用，本文将饱和土体 响应还受到土体本构关系、小变形条件下几何方程 视为均质、连续、各向异性的弹塑性多孔介质，将 的约束，并遵循与渗流场相互作用的广义 Terzaghi 有效应力系数视为二阶张量，根据虚位移原理推导 有效应力原理，即 出各向异性饱和土体直接耦合的有限元法计算公 s D ep e a p (2)