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

《ABAQUS在饱和非饱和渗流分析中的应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、http://www.paper.edu.cnABAQUS在饱和－非饱和渗流分析中的应用徐海奔河海大学水工结构工程专业，南京(210024)E-mail：hohaixhb@163.com摘要：本文首先对大型通用有限元软件ABAQUS在土石坝渗流分析中的应用进行分析，着重从多孔介质的饱和渗流，非饱和渗流及二者的混合问题（渗流自由面的计算）等方面论述。结合一个土石坝库水位下降时二维渗流计算实例，考虑流体重力作用下，采用非线性定律求解总孔隙压力及库水位下降过程渗流自由面变化过程。关键词：非饱和；渗流；ABAQUS；土石坝；自由面1．引言ABAQUS大型通用有限元软件,在

2、我国土木工程结构分析方面应用日益广泛。本文对它在土石坝渗流计算分析中的应用进行评述。近年来，在国内外随着孔隙介质非饱和渗流和土体饱和渗流理论的发展，人们逐渐意识到堤坝稳定性与非饱和区渗流作用密切相关。在研究堤坝非饱和渗流问题时，主要采用数值模拟的方法。长期蓄水的土坝，当库水位以太快的速度下降时，坝体内孔隙水压力常常不能很快消散，因而坝体的浸润线高于上游库水水位。在这种情况下，渗流的动水压力或渗透力的作用对上游坝坡造成浮起及下滑的趋势，甚至酿成滑坡事故。因此在实际工程中必须防止因库水位下降速度太快而导致这类事故发生。为进行上游坝坡的稳定分析，需要确定库水位下降过程中

3、各时段坝体浸润线的位置，也就是通常所说的进行土坝不稳定渗流计算。坝体浸润线下降的速度，一般决定于库水位下降的速度V、土坝坝体渗透系数k以及土[1]体的给水度u等因素，与坝体的结构形式特别是坝体及地基上游面的排水条件也有很大关系。2．ABAQUS在均质土坝饱和－非饱和渗流计算原理在饱和土壤中，引起水分转移的力是重力和水的压力。在非饱和土中，支配着土壤水在液态下整体转移的是重力和水的表面张力。Richards等曾在1931年就证明非饱和土中的渗流[2]与饱和土一样符合达西定律和连续方程。若将达西定律代入连续方程（忽略渗透过程中总应力的改变和土颗粒骨架的变形）并以总水头

4、h作为未知量，当渗透的主方向与坐标轴一致时，非饱和土渗流的二维微分方程就可表示为：∂⎛∂h⎞∂⎛∂h⎞∂θw⎜kx⎟+⎜⎜ky⎟⎟=（1）∂x⎝∂x⎠∂x⎝∂y⎠∂t式中，k，k分别为x，y方向的渗透系数；θ为体积含水量；h为总水头；t为时间。xywuw令y为位置水头，则：h=+y，若m为土水特征曲线斜率，则：wγw∂θ=m∂u=mγ∂()h−y。式（1）就可以写为：wwwww∂⎛∂h⎞∂⎛∂h⎞mwγw∂(h−y)⎜kx⎟+⎜⎜ky⎟⎟=（2）∂x⎝∂x⎠∂x⎝∂y⎠∂t因为y为常数，式（2）可简化为：-1-http://www.paper.edu.cn∂⎛∂h

5、⎞∂⎛∂h⎞mwγw∂h⎜kx⎟+⎜⎜ky⎟⎟=（3）∂x⎝∂x⎠∂x⎝∂y⎠∂t[3]上式即为二维饱和－非饱和渗流方程。从式（3）也可以看出，非饱和土中的渗透系数不再是常数，而是含水量的一个函数，此函数被称为非饱和土的渗透系数函数。因此，要进行非稳定渗流分析，需具备以下基本条件（1）材料的渗透系数函数，包括负孔压区的渗透系数函数；（2）库水位下降速度；（3）确定边界条件。当上述条件已知时，就可以对式（3）进行求解，从而得到饱和－非饱和土的非稳定渗流场。3．均质土包饱和－非饱和渗流分析某水库大坝为均质坝，坝高12m，坝顶宽4m，坝底宽52m，坝体边坡坡度为1:2，

6、在坝下游设有排水棱体（底宽4m），底部其他为不透水层，下游没有水。水库上游水位随时间的变化过程如图1，假定历时6个月时间将水库蓄满水，此时水位高11m，并保持该水位5年，5年后又花6个月时间将水库蓄水排空。由于该土坝足够长，所以视为平面应变问题，采用孔压/位移耦合的CPE8RP平面应变单元，模型有限元网格见图2。图1水位随时间变化过程图2有限元网格3水的重度为10KN/m，土坝土体在饱和时的渗透系数为0.06912m/d，土体非饱和时，土的渗透系数满足以下基本假定：渗透系数随饱和度的减小而降低，并和饱和度之间呈^[4]3次方函数关系。其渗透系数k=kk,其中k是基

7、质吸力或饱和度的函数，如图3所示。ss非饱和土对应的孔隙水压力与饱和度关系可由土水特征线（吸湿/干燥曲线）来描述，如图4所示，假定吸湿曲线和干燥曲线重合。初始状态假定其饱和度为0.167，初始孔压为负值，[5]并随高度线性变化，堤顶孔压为-120kPa，初始孔隙比为1.0。-2-http://www.paper.edu.cn图3渗透系数随饱和度的变化曲线图4土水特征线图5为堤坝土体初始孔压分布图，图6－图8为蓄水不同时期坝体孔隙压力分布图，图9为渗流到达稳定时坝体孔隙压力分布图，图10－图13为排水不同时期坝体孔隙水压分布图：（其中零孔压线即为渗流自由面）图5：初

8、始孔压分布