含孔压间断面的饱与一非饱与渗流数值模拟.pdfVIP

含孔压间断面的饱和一非饱和渗 流数值模拟 吴梦喜1 丁留谦2 (1中国科学院力学研究所北京 100190； 2中国水利水电科学研究院 北京 100038) 【摘 要】 非饱和区渗流介质中存在吸力，其大小取决于介质本身的孔隙特征与含水量。 非饱和介质存在一个最小水连通含水量如，含水量低于靠时，吸力消失，介质中的水流动 不符合一般的渗流规律或无渗流。非饱和区两种介质接触面孔隙水压力一般是不连续的。心 墙土石坝工程中的渗流是典型的情况，土质心墙与下游坝壳的接触面上孔压间断，下游坝壳 砂卵石或堆石等粗颗粒介质中的％基本上等于饱和含水量，这些介质中非饱和渗流不符合 一般的渗流规律。本文提出了含甜压问断面渗流场的有限元模拟方法，并给出了二维心墙土 石坝稳定渗流的算佣，模拟结果与土石坝中的渗流现象相符。 【关键词】 非饱和渗流 孔压间断 有限单元法 心墙土石坝 吸力 1 引言 在渗流非饱和区，孑L隙水压力在两种介质的接触面上一般是不连续的，这就是所谓的 考虑了接触界面两侧孔隙水压力的不连续问题。为避免同一节点存在两个孑L压值，在界面 上不设置节点，而是将差分节点对称布置于界面两侧，界面两侧节点之间的渗透系数，综 合两种介质渗透系数取值，而界面上的孔压根据上下两侧节点外推而获得。这是最早考虑 分层介质非饱和渗流界面孔压不连续问题的文献。目前，此方法的发展仍然主要局限于分 层介质的一维非饱和渗流差分计算∽_3]。就笔者所知，对于非饱和渗流计算的有限元方 法，目前很少有孔压间断面模拟研究的报导，更未见有商业软件提供孔压间断面模拟 功能。 本文的主要目的是要介绍一种模拟饱和一非饱和渗流场孔隙水压力间断面的有限元模 拟方法。该方法已经应用于实际土石坝工程三维稳定渗流问题的有限元计算。 2饱和非饱和渗流问题的描述 渗流的基本方程一般采用基于达西(Darcy)定律和质量平衡原理推导出的理查兹 319 (Richards)方程‘4] [愚4(驴)驴。，]，t+忌，3．t(驴)=口 (1) 式中：f=1，2，3为空间直角坐标下标，量纲为L，j为高程；惫。为土体渗透系数张量， 量纲为LT～，曲为孔隙水压力水头，量纲为L；口为体积含水量，无量纲，变量上边的一 点代表对时间t的偏导数，￡的量纲为T。 理查兹方程中孔隙水压力妒和土体含水量口均为未知变量，非饱和区中渗透系数也随 着含水量的变化而变化。确定了孔隙水压力与含水量之间的关系和渗透系数与含水量之间 的关系后式(1)才能求解。 土体含水量与饱和含水量的比值称为饱和度。孔隙水压力与饱和度关系常采用Van Genuchten公式(VG公式)[5]表示，其公式为 S。(≯)：骅：j‘1+I口”驴㈩一_’Ko (2) ’ 伊，一以 【 1， 驴≥0 式中：S。为有效饱和度；以为饱和体积含水量；只为残余体积含水量；口。为和平均粒径 大小有关的参数；咒。为和粒径均匀性有关的参数，m。=1—1／n。。 非饱和土透水性与饱和度的关系常使用Mualem模型[63描述，其公式为 K(S。)一K；S：胆[1一(1一S：厢t)””]2 (3) 式中：K为渗透系数；K，为饱和渗透系数；其余同前。 基于式(1)～式(3)，可以构建有限差分、有限体积和有限单元等数值算法。对于 饱和一非饱和渗流问题，一般采用≯作为未知变量，式(1)中口也转化为妒的函数表示。 一般数值方法包括有限单元法的求解区域也是在整个介质区域进行的，隐含整个介质区域 均有渗流存在，而且渗流满足式(1)关系的假定。这种方法的优点是不需要在迭代过程 中确定求解区域。对于复杂的水电站地下厂房和坝区稳定渗流分析问题，虽然我们不太关 心非饱和区的渗流情况，但我们仍然喜欢采用饱和一非饱和渗流的分析方法。 时，S，一O，惫一O。渗流介质中非饱和区的孔隙水压力驴为负值，其绝对值称为吸力，由 基