1.引言
我们在前一篇文章中介绍了水闸闸区饱和渗流计算的一般方法，从本篇开始我们将开始介绍另一种岩土渗流问题：非饱和渗流。
如前篇所述，稳定渗流问题的实质是基于边界条件和初始条件的拉普拉斯方程的定解问题。不同于完全饱和渗流，非饱和渗流的边界条件（顶部自由水面边界条件）是未知的，它与流经岩土体多孔介质的水流有关。以堤防和土石坝为代表的非饱和渗流问题如图1所示。
图1 非饱和渗流示意图
确定自由水面边界条件可用多种方法通过迭代得到，作为最一般方法，可通过变形网格，更新网格节点坐标使饱和区域的顶部边界和自由水面一致，这个自由水面就是浸润线，渗流只发生于浸润线以下的饱和区域。
2.工程案例
下面分别以土石坝中的心墙坝和斜墙坝为例进行说明。
2.1.心墙坝
本大坝最大坝高22m，顶宽8.5m，梯形心墙底宽17m，顶宽8.5m，高18m。假定建基面以下为相对不透水地基，大坝的典型断面结构图及网格划分图如图2，图3所示。
图2 心墙坝典型断面结构图
图3 心墙坝典型断面网格划分图
在上游水位18m，下游水位1.8m，按照给定的心墙和坝壳渗透系数的情况下，计算得到的稳定渗流场浸润线、渗压水头分布、单宽渗流流量如下图所示。
图4 心墙坝稳定渗流场单宽流量、浸润线及渗压水头云图
总水头分布、渗流速度分布、边界流量、材料编号、节点编号、单元编号分别如图5-10所示。
图5 心墙坝稳定渗流场总水头云图
图6 心墙坝稳定渗流场渗流速度云图
图7 心墙坝稳定渗流场边界流量
图8 心墙坝材料编号（心墙1，坝壳2）
图9 心墙坝节点编号
图10 心墙坝单元编号
用户可点击点采样工具，查询任意位置的属性值，如图11所示。
图11 点采样
用户可根据坝体的设计材料物理力学参数或者根据反演得到的参数进行任意上下游水位条件的渗流计算，如下图12所示，当点击提交计算，将与服务器之间建立websocket长连接，实时的反馈计算进度，如图13所示。
图12 计算输入控制
图13 实时反馈计算进度
2.2斜墙坝
本大坝为粘土斜墙坝，最大坝高41m。上游坝坡设有一级马道和一处变坡，坡比自上而下分别为 1:2.5、1:2.5 和 1:3.0。下游坝坡设有两级马道，坡比自上而下分别为 1:2.25、1:2.25 和 1:2.5。上下游坝坡分别采用干砌石和碎石进行护坡。大坝采用粘土斜墙防渗，斜墙顶部厚 2m，底部厚 18m；坝基采用粘土铺盖加截水槽防渗，铺盖长 43.5m，厚 2m；截水槽中心线位于坝轴线上游 40m，槽底与基岩相连，底宽 10m，深 10m，并与粘土铺盖和粘土斜墙连为一体。大坝典型断面的结构图及网格划分图分别见图14，图15所示。
图14 斜墙坝典型断面结构图
图15 斜墙坝典型断面网格划分图
在上游水位62m，下游水位38m，按照给定的斜墙和坝壳渗透系数的情况下，计算得到的稳定渗流场浸润线、渗压水头分布、单宽渗流流量如下图所示。
图16 斜墙坝稳定渗流场单宽流量、浸润线及渗压水头云图
斜墙坝其他变量的云图展示，查询，以及实时计算过程与心墙坝类似，限于文档篇幅，不再赘述，读者可移步至文末视频了解。
3.结语
本文以心墙坝与斜墙坝这两大典型土石坝结构为研究载体，系统介绍了非饱和渗流问题的求解过程，以及计算结果的前端展现方法。该方法作为一种物理分析手段，可与测管水位、集水堰渗漏量等监测数据相互佐证，为堤坝渗流安全提供有力保障。
非饱和渗流稳定渗流场中的渗流作用力，是一种空间分布不均匀的体力。它作用于土颗粒之上，与土体自重类似，同样会对堤坝的抗滑稳定性产生显著影响。
下一步，我们将介绍非饱和渗流与坝体变形、稳定性分析的耦合计算，进而构建更为完整的工程安全评估体系。
