ABAQUS在饱和非饱和渗流分析中的应用.pdfVIP

ABAQUS 在饱和－非饱和渗流分析中的应用 徐海奔 河海大学水工结构工程专业，南京 (210024) E-mail：hohaixhb@163.com 摘 要：本文首先对大型通用有限元软件ABAQUS 在土石坝渗流分析中的应用进行分析， 着重从多孔介质的饱和渗流，非饱和渗流及二者的混合问题（渗流自由面的计算）等方面论 述。结合一个土石坝库水位下降时二维渗流计算实例，考虑流体重力作用下，采用非线性定 律求解总孔隙压力及库水位下降过程渗流自由面变化过程。 关键词：非饱和；渗流；ABAQUS；土石坝；自由面 1．引言 ABAQUS 大型通用有限元软件,在我国土木工程结构分析方面应用日益广泛。本文对它 在土石坝渗流计算分析中的应用进行评述。 近年来，在国内外随着孔隙介质非饱和渗流和土体饱和渗流理论的发展，人们逐渐意识 到堤坝稳定性与非饱和区渗流作用密切相关。在研究堤坝非饱和渗流问题时，主要采用数值 模拟的方法。长期蓄水的土坝，当库水位以太快的速度下降时，坝体内孔隙水压力常常不能 很快消散，因而坝体的浸润线高于上游库水水位。在这种情况下，渗流的动水压力或渗透力 的作用对上游坝坡造成浮起及下滑的趋势，甚至酿成滑坡事故。因此在实际工程中必须防止 因库水位下降速度太快而导致这类事故发生。为进行上游坝坡的稳定分析，需要确定库水位 下降过程中各时段坝体浸润线的位置，也就是通常所说的进行土坝不稳定渗流计算。 坝体浸润线下降的速度，一般决定于库水位下降的速度V 、土坝坝体渗透系数k 以及土 体的给水度u等因素[1]，与坝体的结构形式特别是坝体及地基上游面的排水条件也有很大关 系。 2 ．ABAQUS在均质土坝饱和－非饱和渗流计算原理 在饱和土壤中，引起水分转移的力是重力和水的压力。在非饱和土中，支配着土壤水在 液态下整体转移的是重力和水的表面张力。Richards等曾在1931年就证明非饱和土中的渗流 [2] 与饱和土一样符合达西定律和连续方程 。若将达西定律代入连续方程（忽略渗透过程中总 应力的改变和土颗粒骨架的变形）并以总水头h作为未知量，当渗透的主方向与坐标轴一致 时，非饱和土渗流的二维微分方程就可表示为： ∂ ⎛ ∂h ⎞ ∂ ⎛⎜ ∂h ⎞ ∂θ⎟ w ⎜k ⎟+ k （1） x ⎜ y ⎟ ∂x ⎝ ∂x ⎠ ∂x ⎝ ∂y ⎠ ∂t 式中，k ，k 分别为x ，y方向的渗透系数；θ 为体积含水量；h为总水头；t为时间。 x y w u 令 y 为位置水头，则： h w +y ，若m 为土水特征曲线斜率，则： w γw ∂θ m ∂u m γ ∂(h −y )。式（1）就可以写为： w w w w w ( ) ∂ ⎛ ∂h ⎞ ∂ ⎛ ∂h ⎞ mwγw ∂h −y ⎜ ⎟ （2 ） ⎜k ⎟+ k x ⎜ y ⎟ ∂x ⎝ ∂x ⎠ ∂x ⎝ ∂y ⎠ ∂t 因为y 为常数，式（2 ）可简化为： -1-