🙋
"流网"就是地盘里水流的网格状图吧?那为什么要画这个呢?
🎓
是的,当地下水从大坝或导板墙下渗过时,会画成网格状的图形。地盘中的渗透流遵循"拉普拉斯方程",这是一个偏微分方程。在没有计算机的年代,需要用手工求解。所以人们想出了一种方法,让流线和等势线相互垂直,形成"正方形"的网格。这就是流网法。只要画出图来,就能读出渗透流量和危险位置。
🎓
流线就是水粒子从上游流向下游的路径。等势线是全水头相等的点的连线。全水头是"水压+位置势能"。这两族线总是成直角交叉,把渗透区域分成"正方形"的曲线方格。左边的图中,蓝色线是流线,橙色线是等势线。用滑块改变 Nf 和 Nd,网格的细密度就会改变。
🙋
明白了。那渗透流量怎么算出来的?公式是 q = k·H·(Nf/Nd)…
🎓
画出流网后,只需数出两个整数:流路数 Nf 和等势降数 Nd。Nf 是被流线包围的带的本数,Nd 是被等势线分隔的段数。然后用 q = k·H·(Nf/Nd) 计算单位宽度的渗透流量。乘以构筑物的深度就得到全部的漏水量。比如默认值的情况,每降下降时水头失失 0.5m,全体漏水量是 0.6 L/s。
🙋
不只是漏水量,还有"安全系数"对吧。这个是什么的安全系数?
🎓
这是防止管涌的安全系数。看下游的出口附近,等势线紧密地集中在一起,动水坡度就变陡了。当出口的动水坡度 i_exit 接近临界动水坡度(砂土约为 1)时,向上的渗透力会超过砂粒的有效重量,砂粒就会浮起来,地盘会像"沸腾"一样破坏。这就是管涌,也叫做沸腾。实际上很多大坝决口和开挖事故都是从这里开始的。所以我们计算 FOS = i_cr/i_exit,通常要确保 FOS ≥ 3。
🎓
最直接的办法就是"加深导板"。这样水的通路就会变长,等势降数增加,出口的坡度就会变缓。或者在下游铺上导流毯来增加出口的流场长度,或者用井点降低地下水位。看左边"出口流场长度 l"的参数,增大这个值,安全系数会显著上升。流网是用一张图同时告诉我们漏水量和"最危险的地方"的非常强大的工具。
流网是绘制流线和等势线直角相交形成"正方形"网格的图。单位宽度的渗透流量用 q = k·H·(Nf/Nd) 计算。k 是透水系数,H 是上下游水位差,Nf 是流路数(被流线包围的带数),Nd 是等势降数(被等势线分隔的段数)。Nf 和 Nd 只需绘制流网后数出即可,乘以构筑物的深度就得到全部渗透流量。
在下游出口附近,等势线密集集中,动水坡度变陡。当出口动水坡度 i_exit 接近临界动水坡度 i_cr(砂土约为 1.0)时,向上的渗透力超过砂粒的有效重量,地层浮起形成"沸腾"般的破坏。这就是管涌或沸腾。许多大坝和开挖事故都源于此处。流网用一张图清楚地显示等势线最密集、最危险的位置。
单位降水头损失为 Δh = H/Nd。设出口流场长度为 l,则出口动水坡度为 i_exit = Δh/l = (H/Nd)/l。管涌安全系数为 FOS = i_cr/i_exit。通常确保 FOS ≥ 3,当安全系数不足时,采取加深导板、铺设导流毯或井点降水等措施来增加出口流场长度。
流网法是用手工求解二维拉普拉斯方程的经典方法,在没有计算机的时代就已广泛应用。如今可用渗流有限元法进行详细分析,但流网法作为概算和验证工具仍然有效。对于边界条件简单、地盘等向均质的情况,手工绘制的流网精度足以满足工程要求。如果有限元结果与流网概算相差很大,应怀疑网格划分或透水系数、边界条件的输入错误。
土坝·堆石坝设计:用土或岩石筑造的大坝中,堤体和基础地盘中通过的渗透流用流网评估。既要抑制漏水量,也要防止下游坡脚出现管涌,保持动水坡度在临界值以下。坝心(防渗区)的厚度和基础的防渗墙(截水)的深度根据流网中能获得足够的等势降数来确定。
导板截水和开挖工事:在河川和港湾工事中,用导板挡水并在内侧进行开挖时,绕过导板的渗透流可能导致开挖底面发生"沸腾"破坏。设计时需用流网检讨导板的根入深度,确认出口动水坡度的安全系数。当安全系数不足时,需加深导板或用井点强制降低地下水位。
堤·水门·取水构造物:混凝土堤或水门的基础地盘中,从上游到下游的渗透流会产生扬压力(托起构筑物的水压)。从流网的等势线可以读出基础底面的扬压力分布,用于结构物的滑动、倾覆稳定性计算。导流毯和截水墙的配置是控制渗透路径的设计手段。
CAE渗流分析的验证:在用 SEEP/W 或 PLAXIS 等渗流有限元进行详细分析之前或之后,用流网手算来确认 q 和 i_exit 的数量级。如果有限元结果与流网概算相差很大,说明可能存在透水系数单位错误或边界条件(水头边界、不透水边界)设置错误。简单的概算是复杂数值分析的可信性基础。
首先,有个很大的误解是"流网的正方形必须是真正的正方形"。流网的方格实际上是曲线正方形,只要对边的中点长度大致相等、对角线互相直交就足够了。不必为了追求完美的正方形而反复重画。即使流路数 Nf 和等势降数 Nd 有小数点(比如 3.5),也没问题。重要的是比值 Nf/Nd,即使网格粗糙,这个比值也能稳定地求出来,这正是流网法的优势。
其次,"直接把等向均质地盘的公式用于异向地盘"。q = k·H·(Nf/Nd) 的公式是基于水平和竖直透水系数相等的等向地盘。实际的堆积地盘通常水平透水性强于竖直(kh > kv),具有成层异向性。这种情况需要在"变换断面"上绘制流网,即把水平距离缩小 √(kv/kh) 倍,流量计算用等效透水系数 k' = √(kh·kv)。忽视异向性会导致漏水量和动水坡度偏离很大。
最后,"只要出口动水坡度的安全系数足够就能防止管涌"的误解。本工具计算的是出口面的渗透力与土的重量比较的"隆起(沸腾)"型安全系数。但实际的管涌还有"后退侵蚀型",即土粒从下游逆向被冲蚀,形成扩大的管状水渗通道。这种情况可能从局部弱点、粒度不连续或动物洞穴等处开始,即使平均动水坡度安全也会发生。所以必须根据粒度分布设置滤层材料,并持续监测渗水的浑浊度和流量。