📖 人物：亨利·达西（Henry Darcy）与朱利叶斯·魏斯巴赫 达西是法国工程师，在 19 世纪研究城市供水管道时，通过大量实验总结出沿程损失的规律；德国学者魏斯巴赫把它整理成了现在这个形式。达西公式至今仍是全世界管道工程计算的标准工具——一百多年过去了，还在用。
🎯 本章学习目标
- 理解沿程损失与局部损失的区别，以及它们各从哪来。
- 掌握达西-魏斯巴赫公式，理解每一项的含义。
- 知道沿程阻力系数 λ 在层流（λ=64/Re）和湍流中如何确定，会查莫迪图。
- 掌握局部损失的计算（hj = ζ·v²/2g）。
- 会用「总水头损失 = 沿程损失之和 + 局部损失之和」做管路计算。
- 会用当量长度法折算局部损失，能处理串联、并联管路的流量分配，理解管路特性曲线与工况点。
7.1水头损失：沿程损失与局部损失
第 6 章讲过：真实流体有粘性，机械能沿流动方向不断耗散成热。这个耗散在「水头」语言里就是水头损失 hw。它分两类，来源完全不同。
| 类型 | 符号 | 来源 | 特点 |
|---|---|---|---|
| 沿程损失 | hf | 流体与管壁、流体层间的粘性摩擦 | 沿整个管长均匀产生，管越长越大 |
| 局部损失 | hj | 流经阀门、弯头、突变等「局部障碍」时的涡流、冲击 | 集中在局部位置，与管长无关 |
可以这样理解：沿程损失像「走长路脚底磨损」，越走越多；局部损失像「过路口被绊一下」，集中发生在特定位置。整段管路的总水头损失就是两者相加：
hw = Σhf + Σhj
于是，有粘性的伯努利方程（能量方程）要「带上损失」：
实际流体的能量方程 z₁ + p₁/(ρg) + v₁²/(2g) = z₂ + p₂/(ρg) + v₂²/(2g) + hw。即：上游总水头 = 下游总水头 + 损失（水头沿流「缩水」了）。
损失的物理本质：机械能的「不可逆贬值」
要真正理解 hw，得回到热力学。能量方程的三项 z、p/(ρg)、v²/(2g) 分别是单位重量流体的位置势能、压强势能、动能，合称机械能。粘性耗散把机械能变成热——热能虽然还留在流体里，却再也无法无偿变回「能推动叶轮」的机械能，热力学第二定律保证了这一过程单向、不可逆。所以能量方程里损失项只能出现在下游一侧、恒为正值：任何声称让「越流能量越多」的装置，都等价于永动机。
再看 hw 的单位。它是「米」，准确含义是「每牛顿重的流体损失了多少焦耳」（J/N）。用「水头」而不是「焦耳」记账的好处是直观：把总水头 H = z + p/(ρg) + v²/(2g) 沿管长画成曲线，就是总水头线（能量线）；把 z + p/(ρg) 单独画出来，是测压管水头线。理想流体里总水头线是一条水平线；实际流体里它沿程下降——沿程损失让它「缓降」，局部损失让它「陡降」。一条曲线画完，整段管路哪里费劲、费多少，一目了然。
「缓降 + 陡降」正是 hw = Σhf + Σhj 在几何上的样子。而「把损失一路加起来」之所以合法，也是热力学给的许可：耗散恒为正，同一条确定流线上各段的损失就是沿途耗散的简单叠加，不存在「前面的损失改变后面的损失」的耦合。下面这道例题，先用能量方程把总损失「称」出来——这是一切管路计算的第一步。
📐 例题 1 用能量方程「称」出总水头损失
一高位水池，自由液面比管道出口中心高 12 m，水经长管从出口流入大气。已测得出口流速 v₂ = 3 m/s；水池液面很大，可认为 v₁ ≈ 0；两端均为大气压（表压为零）。求整段管路的总水头损失 hw，并换算成压降。
思路 在「水池液面（断面 1）→ 管道出口（断面 2）」之间列带损失的能量方程，未知量恰好只有 hw 一个。
解 z₁ + p₁/(ρg) + v₁²/(2g) = z₂ + p₂/(ρg) + v₂²/(2g) + hw。取出口中心为基准：z₂ = 0、z₁ = 12 m；p₁ = p₂ = 0（表压）；v₁ ≈ 0。于是 hw = z₁ − v₂²/(2g) = 12 − 3²/(2×9.81) = 12 − 0.459 = 11.54 m。换算成压降：Δpw = ρg·hw = 1000×9.81×11.54 ≈ 1.13×10⁵ Pa ≈ 0.113 MPa。
答：总水头损失约 11.54 m 水柱，合压降约 0.113 MPa。 ⚠ 易错点：一是别忘了扣掉出口的速度水头 0.459 m——损失是「总水头之差」，不是「高度之差」；二是两个断面必须同取表压或同取绝对压，混用会凭空多出一个大气压的误差。
工程视角：管路就是一张「能量预算表」
城市供水、供热管网、消防系统，以及飞机上的燃油与液压管路，本质都是同一张预算表：泵或高位水池注入能量，管路沿途消耗能量，末端剩下的能量必须够用。设计管路的三大任务——定流量、定管径、选泵——全部围绕这张表展开。以飞机液压系统为例（任何机型都遵循同样的物理）：泵的出口压力有限，从泵到作动器之间每一段管路、每一个接头、每一个弯头都在「吃」压力；油温降低时粘度升高、雷诺数下降，沿程损失进一步增大，末端压强可能跌到作动器推不动负载的程度。所以管路设计规范要求逐段核算压降，并给低温工况留足余量。
辨析：水头（m）还是压降（Pa）？ 同一份损失有两种记法：hw 用「米流体柱」，Δpw = ρg·hw 用「帕」。换算时别忘密度：同样 hw = 10 m，水对应 Δp ≈ 98.1 kPa，密度 850 kg/m³ 的液压油对应 83.4 kPa。液体管路习惯用水头；气体管路密度小、折算成水头是天文数字，习惯直接用压降（Pa）。读手册先看单位体系，「米」和「帕」绝不能直接比大小。
7.2达西公式：计算沿程损失
沿程损失怎么定量算？答案是达西-魏斯巴赫公式（简称达西公式）：
达西公式 hf = λ · (L/d) · (v²/2g)
逐项看它的含义：
- L/d：管长与管径之比，反映「管道有多长（相对粗细）」。管越长，损失越大。
- v²/2g：速度水头（第 4 章见过）。流速越大，损失越大。
- λ：沿程阻力系数（摩擦因子），是一个无量纲系数，由流动状态（雷诺数）和管壁粗糙度决定。
达西公式为什么长这样？
直觉上它很合理：沿程损失应当和「管长」成正比、和「动压（速度水头）」成正比。达西公式的骨架其实可以从不严谨但直观的推理看出来——损失正比于壁面摩擦力的功，而摩擦力正比于动压 ρv²/2 乘以接触面积（∝ L·d）。真正难的是那个系数 λ，它把层流、湍流、粗糙度的影响全「藏」进去了。所以沿程损失计算的核心，就是确定 λ。
🧠 一句话记忆 沿程损失 hf = λ·(L/d)·(v²/2g)：管越长、流越快、λ 越大，损失越大。λ 是灵魂。
从壁面切应力看达西公式：λ 的真面目
达西公式还能「倒过来」理解——直接从壁面切应力出发把它推出来。取一段长 L、直径 d 的水平等径直管，把管内整段流体当作一个受力对象：上游压强 p₁、下游压强 p₂（p₁ > p₂），压差 Δp = p₁ − p₂ 作用在过流断面 A = πd²/4 上，是推动力；管壁通过切应力 τw「拽」住流体，作用在整个侧面积 πdL 上，是阻力。定常流动没有加速度，推力与阻力必须平衡：
Δp·(πd²/4) = τw·πdL ⟹ Δp = 4τwL/d
把它与达西公式的压降形式 Δpf = λ·(L/d)·(ρv²/2) 放在一起比对，立刻得到：
λ = 8τw/(ρv²)
这个式子揭示了 λ 的真面目：λ 就是把壁面切应力用动压归一化的结果，度量「壁面拽流体的力」相对「流体自身惯性」有多大。λ 大，说明摩擦相对惯性大、单位动能被摩擦吃掉得多。于是「查 λ」这件事，物理上等价于「问壁面切应力有多大」。整个第 7 章的核心问题——λ 怎么定——其实就是第 6 章遗留的问题「壁面切应力由什么决定」的继续。
辨析：达西 λ 与范宁 f，差 4 倍的坑 文献里还有另一个摩擦系数——范宁系数 f = τw/(ρv²/2)，它把同一个壁面切应力除以动压（ρv²/2）而不是 ρv²/8。两者关系是 λ = 4f：层流时 λ = 64/Re，相应 f = 16/Re。查国外手册、对照文献时务必先看清用的是哪一个系数，否则结果会差整整 4 倍——这是管路计算最经典的「翻车点」之一。
「二次方」的真相：损失到底跟流速几次方成正比？
达西公式表面写着 v²，很多人据此记住「损失与流速平方成正比」。更严谨的说法是：hf ∝ v² 的前提是 λ 为常数，而 λ 本身随 v（通过 Re）变化。把三种流态逐一展开：
- 层流（λ = 64/Re ∝ 1/v）：hf = (64ν/vd)·(L/d)·(v²/2g) = 32νLv/(gd²)，损失与流速一次方成正比——流速翻倍，损失只翻一倍。
- 湍流光滑区（Blasius 公式 λ = 0.3164/Re0.25）：hf ∝ v1.75，接近平方但略低。
- 湍流完全粗糙区（λ 只取决于 ε/d，与 Re 无关）：hf 严格 ∝ v²，「平方阻力区」的名字就是这么来的。
所以「损失与流速的 1~2 次方成正比」不是含糊其辞，而是有明确物理来源的区间：Re 越大、壁面越粗糙，越接近平方律。这也解释了日常经验：水龙头从开一半拧到全开，管网里的损失不是按线性增长，而是按接近平方的速度飙升。
📐 例题 2 光滑管湍流：Blasius 公式 + 达西公式
20 ℃ 的水（ν = 1.0×10⁻⁶ m²/s）在内壁光滑的直管内流动，d = 100 mm，L = 50 m，v = 1.0 m/s。先验证流态，再用 Blasius 公式取 λ，求沿程损失 hf 与压降 Δpf。
思路 Re 判流态与公式适用范围 → Blasius 取 λ → 代达西公式 → 换算压降。
解 ① Re = vd/ν = 1.0×0.1/1.0×10⁻⁶ = 1×10⁵ > 4000，为湍流，且恰在 Blasius 公式适用上限（约 4000 ≤ Re ≤ 10⁵）附近，工程上可直接使用（Re 更高时换用科尔布鲁克公式，见 7.3 节）。 ② λ = 0.3164/Re0.25 = 0.3164/(1×10⁵)0.25 = 0.3164/17.78 = 0.0178。 ③ hf = λ(L/d)(v²/2g) = 0.0178×(50/0.1)×(1²/19.62) = 0.0178×500×0.0510 = 0.454 m。 ④ Δpf = ρg·hf = 1000×9.81×0.454 ≈ 4.45×10³ Pa ≈ 4.45 kPa。
答：hf ≈ 0.454 m 水柱，压降约 4.45 kPa。 ⚠ 易错点：(L/d) = 500 别算错数量级；Blasius 公式只对「水力光滑管」成立——若这根管子是旧的镀锌管（ε/d ≈ 0.0015），λ 可能要按莫迪图查到 0.024 左右，损失会大出六成。
工程视角：为什么「加粗一档」这么值钱
把 v = 4Q/(πd²) 代入达西公式，损失可以改写成对流量的形式：hf = λ·(L/d)·(8Q²/(π²gd⁴)) = 8λLQ²/(π²gd⁵)。注意分母上的 d 的五次方：管径只增大 50%（d → 1.5d），同样流量下损失就降到原来的 1/1.5⁵ ≈ 13%（λ 变化暂不计）。这就是为什么市政与工业设计宁可初投资买粗管、再省下几十年泵的电费——管路「全寿命成本」里，电费往往远超管材本身。反过来看飞机：液压管径每加粗一档，阻力省下的泵功率有限，但管子重量、占用空间、安装成本立刻上升；于是航空液压更倾向「高压小流量」路线——系统压力越取越高，本质是用高「压强势能」替代大「动能」，让同样的功率流过更细的管子。同样的 d⁵ 规律，在地面上和在天上给出相反的设计取向：公式不变，约束变了，最优解就变了。
衔接 到这里，沿程损失的计算流程已经完整：判流态（第 6 章的 Re）→ 定 λ（下一节）→ 套达西公式。7.3 节专门解决中间那一步——λ 究竟怎么定、为什么层流能算而湍流只能查。
7.3沿程阻力系数 λ 与莫迪图
λ 怎么定？分两种情况。
层流：λ 有精确解
层流时，速度剖面是抛物线，可以严格推导出（用牛顿内摩擦定律 + 受力平衡）：
λ = 64 / Re
这个结论漂亮而精确：层流的沿程阻力系数只和雷诺数有关，与管壁粗糙度无关。原因是层流里流体「分层滑动」，贴壁那层把壁面粗糙度「抹平」了，粗糙度影响不到主流。
这个 64/Re 并不是实验凑出来的经验数，而是一步步推出来的解析解。推导值得完整走一遍——它是「牛顿内摩擦定律 + 受力平衡」联手解题的典范，也是本章最重要的一段数学。
层流速度剖面：完整推导（哈根–泊肃叶流动）
第一步：切应力沿半径的分布。在半径 R 的水平直管内，取一个半径 r（r < R）、长 L 的同心圆柱形流体做受力分析：两端压差 Δp 作用在圆断面 πr² 上是推力，侧面切应力 τ(r) 作用在 2πrL 的圆柱面上是阻力。定常层流无加速度，二者平衡：
Δp·πr² = τ(r)·2πrL ⟹ τ(r) = Δp·r/(2L)
切应力沿半径线性分布：轴线处为零，壁面处最大，τw = ΔpR/(2L) = Δp·d/(4L)。注意这个结论对层流、湍流都成立（只要定常且等径），它只用了受力平衡。
第二步：用牛顿内摩擦定律积分出速度。层流中 τ = −μ·du/dr（负号因为速度 u 从轴线向壁面递减）。代入上式并积分：
−μ·(du/dr) = Δp·r/(2L) ⟹ u(r) = Δp·(R² − r²)/(4μL)
积分中用了壁面不滑移条件 u(R) = 0。结果是一条抛物线——层流管流的速度剖面是旋转抛物面，最大速度在轴线上：umax = ΔpR²/(4μL)。
第三步：流量与平均速度。对剖面沿断面积分得流量 Q = πΔpR⁴/(8μL)，这就是著名的哈根–泊肃叶公式（Hagen–Poiseuille，1839/1840 年由两人各自独立发现）。断面平均速度为：
v = Q/A = ΔpR²/(8μL) = umax/2
第四步：反解压降，对出 λ。由上式 Δp = 8μLv/R² = 32μLv/d²。代入达西公式的压降形式 Δp = λ(L/d)(ρv²/2)，两边约去公因子：
32μLv/d² = λ(L/d)(ρv²/2) ⟹ λ = 64μ/(ρvd) = 64/Re
推导同时送出两个「免费」的工程判据：① 层流剖面 umax/v = 2——实测管中心流速与平均流速之比接近 2，即可判断是层流（湍流剖面饱满得多，这个比值只有约 1.2~1.3）；② 层流压降 Δp ∝ v 一次方——与 7.2 节「名义二次律」的辨析完全对上。
📐 例题 3 液压管路层流损失：λ = 64/Re 全链条
某液压系统供油管内径 d = 25 mm、长 L = 20 m；油液密度 ρ = 850 kg/m³、运动粘度 ν = 1.0×10⁻⁴ m²/s；平均流速 v = 3.0 m/s。判断流态，求沿程损失 hf 与压降 Δpf，并用哈根–泊肃叶公式复核。
思路 Re 判流态 → 层流用 λ = 64/Re → 达西公式算 hf → Δpf = ρg·hf，最后用 Δp = 32μLv/d² 独立验证一遍。
解 ① Re = vd/ν = 3.0×0.025/1.0×10⁻⁴ = 750 < 2000，为层流。 ② λ = 64/750 = 0.0853。 ③ hf = λ(L/d)(v²/2g) = 0.0853×(20/0.025)×(3²/19.62) = 0.0853×800×0.459 = 31.3 m（油柱）。 ④ Δpf = ρg·hf = 850×9.81×31.3 ≈ 2.61×10⁵ Pa ≈ 0.26 MPa。 ⑤ 验证：μ = ρν = 850×10⁻⁴ = 0.085 Pa·s；Δp = 32μLv/d² = 32×0.085×20×3/0.025² = 2.61×10⁵ Pa ✓ 两条路完全一致。
答：层流；hf ≈ 31.3 m 油柱，Δpf ≈ 0.26 MPa。 ⚠ 易错点：一是 31.3 m 是油柱不是水柱，换算压降必须用油的密度 850 而不是 1000；二是注意机理方向——若油温下降（ν 变大），Re 变小、λ 变大，损失反而更大，低温液压系统「压降暴涨」正是这个原因，与「粘稠=润滑好」的直觉相反。
湍流：λ 依赖 Re 和粗糙度
湍流的情况复杂得多，λ 同时依赖雷诺数 Re 和相对粗糙度（管壁粗糙凸起高度 ε 与管径 d 之比 ε/d）。湍流时流体剧烈掺混，壁面粗糙凸起会「刺」进流核，显著增加阻力。工程上把这些关系画成一张著名的图——莫迪图（Moody chart）：横轴 Re，纵轴 λ，一簇曲线对应不同相对粗糙度 ε/d。查图流程：
- 算出 Re 和相对粗糙度 ε/d。
- 在莫迪图上找到对应曲线，读出 λ。
关键对比 层流 λ=64/Re（只与 Re 有关，可精确算）；湍流 λ 依赖 Re 和 ε/d（查莫迪图或用经验公式）。记住：粗糙度只在湍流里才重要。
湍流的 λ：三个阻力区与一个看不见的「底层」
莫迪图上那「一簇曲线」看似杂乱，其实有清晰的分区。理解分区的钥匙是一个看不见的结构：粘性底层（层流底层）。即使在剧烈的湍流中，紧贴壁面也有一层极薄的流体几乎保持层流状态（它的来历第 9 章专门讲），厚度通常只有几十微米到零点几毫米。拿它与壁面粗糙凸起高度 ε 比一比，湍流就分出三个区：
- 水力光滑区：粘性底层完全盖住粗糙凸起，主流根本「感觉不到」粗糙度，λ 只与 Re 有关，Blasius 公式适用。
- 过渡区：Re 增大使粘性底层变薄，部分凸起「露出头来」刺进湍流核心，λ 同时依赖 Re 和 ε/d——莫迪图上的曲线簇正是从这里开始分开的。
- 完全粗糙区（平方阻力区）：凸起几乎全部裸露、层层激起小涡，λ 只取决于 ε/d，与 Re 无关，损失严格与 v² 成正比。
所以「光滑管 / 粗糙管」不是看管子摸起来光不光，而是看粘性底层厚度与凸起高度的相对大小——同一根管子，Re 小时可能处于光滑区，Re 大时进入过渡区甚至粗糙区。在莫迪图上沿同一根 ε/d 曲线从左下往右上走，就是这一过程本身。
| 管材 | ε（mm） | 备注 |
|---|---|---|
| 玻璃、塑料拉制管 | ≈ 0.0015 或更小 | 常直接按水力光滑处理 |
| 普通商用钢管 | 0.045 | 莫迪图基准数据，最常用 |
| 镀锌铁管 | 0.15 | 通风、消防水管常见 |
| 铸铁管 | 0.26（旧管可达 1 以上） | 城市老旧管网损失逐年增大 |
| 混凝土管 | 0.3 – 3 | 施工质量离散度大 |
📐 例题 4 商用钢管湍流：查莫迪图并用科尔布鲁克公式核对
水（ν = 1.0×10⁻⁶ m²/s）以 Q = 20 L/s 流过 d = 200 mm 的普通商用钢管（ε = 0.045 mm），管长 L = 300 m。求沿程损失。
思路 Q → v → Re → ε/d → 查莫迪图（或用科尔布鲁克公式迭代）得 λ → 达西公式。
解 ① v = 4Q/(πd²) = 4×0.02/(π×0.2²) = 0.637 m/s。 ② Re = vd/ν = 0.637×0.2/1.0×10⁻⁶ = 1.27×10⁵，湍流，处于过渡区。 ③ 相对粗糙度 ε/d = 0.045/200 = 2.25×10⁻⁴。 ④ 用科尔布鲁克公式 1/√λ = −2lg[ε/(3.7d) + 2.51/(Re√λ)] 迭代（或查莫迪图）：取初值 λ = 0.02 迭代两轮收敛，得 λ ≈ 0.0184，莫迪图上读数约 0.018~0.019，一致。 ⑤ hf = λ(L/d)(v²/2g) = 0.0184×(300/0.2)×(0.637²/19.62) = 0.0184×1500×0.0207 = 0.571 m。
答：沿程损失约 0.57 m 水柱。 ⚠ 易错点：① Re 的特征长度用管径 d，别把流量当速度直接代入；② 科尔布鲁克公式对 λ 是隐式——λ 出现在方程两边，必须迭代求解，这也是为什么工程上更愿意直接查莫迪图：一条曲线顶一个隐式方程；③ 若忘了单位换算（mm vs m），ε/d 会差 1000 倍，查图结果完全不可用。
为什么工程公式会「查无解析解」？ 层流能用 64/Re，是因为牛顿内摩擦定律给出封闭方程组；湍流的掺混应力（雷诺应力）至今没有普适的解析表达式，λ 只能靠实验 + 半经验公式（如科尔布鲁克公式）或直接查图。1944 年莫迪把当时积累的实验数据整理成一张图，使管路计算「一张纸搞定」，这张图被称为工程水力学最著名的图之一。近年也常用显式近似公式（如 Haaland 公式）替代隐式迭代，误差在工程精度内。
衔接第 9 章 「粘性底层」是提前出场的主角：它其实是边界层贴壁最内层的结构。为什么湍流反而有一个层流内层、它随 Re 如何变薄、它和边界层分离有什么关系，第 9 章专门展开。这里只需记住一句：莫迪图上 λ 的一切复杂行为，几乎都是「粘性底层厚度 vs 粗糙凸起」这一个问题在主导。
7.4局部损失与管路计算
流体流经弯头、阀门、管径突变等「局部障碍」时，流线被迫急剧改变，产生涡流和冲击，消耗额外的能量，这就是局部损失 hj。
局部损失公式
局部损失通常写成速度水头乘以一个系数：
局部损失公式 hj = ζ · v²/(2g)，其中 ζ 是局部阻力系数，由障碍物的几何形状决定，通常查表得到。
举例：一个 90° 弯头 ζ ≈ 0.9；一个全开的闸阀 ζ ≈ 0.2；管径突然缩小/扩大也有各自的 ζ 值。ζ 越大，这个障碍越「费劲」。这些系数都是工程师通过实验测出来、列成手册的。
局部损失的「重灾区」：突扩与突缩
所有局部障碍里，突扩最值得单独研究——它是局部损失中少有的解析解，推导本身就是一次漂亮的「动量方程 + 能量方程」联合作战。设管径从 d₁ 突扩到 d₂，流速从 v₁ 降到 v₂，扩口死角里形成环形旋涡区。取「突扩断面 1 → 扩后流动恢复平直的断面 2」为控制体，做两个关键近似：① 旋涡区很短、壁面摩擦可忽略；② 突扩断面上（含旋涡区环形面）压强仍近似为 p₁——死角里的「死水」把上游压强传递到扩口端面，实验证实误差不大。
动量方程给出压差与流速变化的关系：(p₁ − p₂)/(ρg) = v₂(v₂ − v₁)/g。能量方程（带损失）给出 hj = (p₁ − p₂)/(ρg) + (v₁² − v₂²)/(2g)。把前者代入后者：
hj = [2v₂(v₂ − v₁) + v₁² − v₂²]/(2g) = (v₁ − v₂)²/(2g) （Borda–Carnot 公式）
损失恰好是「速度突降」对应的动能——流体撞进死角减速时，动能没有回收成压能，而是全部耗散成涡流与热。把 Borda–Carnot 公式分别用 v₁、v₂ 改写成 hj = ζ·v²/2g 的形式，就得到两种系数记法：
ζ₁ = (1 − A₁/A₂)²（配大流速 v₁） ζ₂ = (A₂/A₁ − 1)²（配小流速 v₂）
突缩则相反：流体先收缩成一个比小管断面还小的「缩脉」（vena contracta，锐缘时面积约为小管的 0.6 倍），然后再扩散充满小管。所以突缩损失的大头不在「缩」而在缩脉之后的「内突扩」段——本质上还是突扩损失。工程上突缩通常按下游小管流速 v₂ 记：锐缘突缩/进口 ζ ≈ 0.5；进口做成圆角（渐缩）可降到 0.04~0.2——这就是为什么好的设计总要给进口「倒圆角」。
📐 例题 5 突扩损失：同一个损失、两套系数
管径从 d₁ = 100 mm 突扩到 d₂ = 200 mm，小管流速 v₁ = 2.0 m/s。求突扩局部损失，并分别用「小管系数 ζ₁」「大管系数 ζ₂」各算一遍，验证两者一致。
思路 面积比 A₁/A₂ = (d₁/d₂)² = 0.25；先由 Borda–Carnot 公式算出损失真值，再检验两种系数记法。
解 ① 由连续性 v₂ = v₁·A₁/A₂ = 2.0×0.25 = 0.5 m/s。 ② Borda–Carnot：hj = (v₁ − v₂)²/(2g) = (2.0 − 0.5)²/19.62 = 2.25/19.62 = 0.115 m。 ③ 小管系数：ζ₁ = (1 − A₁/A₂)² = (1 − 0.25)² = 0.5625；hj = ζ₁·v₁²/2g = 0.5625×(4/19.62) = 0.115 m ✓ ④ 大管系数：ζ₂ = (A₂/A₁ − 1)² = (4 − 1)² = 9；hj = ζ₂·v₂²/2g = 9×(0.25/19.62) = 0.115 m ✓
答：突扩损失约 0.115 m 水柱；ζ₁ = 0.5625（配大流速 v₁）与 ζ₂ = 9（配小流速 v₂）给出同一物理结果。 ⚠ 易错点：局部阻力系数必须与「配套流速」成对使用——手册中突扩的 ζ 既可以按上游也可以按下游速度定义，查表时务必看清；ζ₂ = 9 看着吓人，只是因为配的是小流速。另注意方向：突缩的 ζ 习惯上按下游（小管）流速计。
管路计算：把所有损失加起来
一段真实的管路，往往由若干直管段 + 若干弯头阀门组成。总水头损失就是：
hw = Σ λ·(L/d)·(v²/2g) + Σ ζ·(v²/2g)
有了总损失，就能用「带损失的能量方程」解出需要的泵扬程、允许流速、或管道尺寸。这是给水、供热、油气输送等工程的日常工作。
当量长度法：把「绊脚的路口」折算成「多长的直路」
沿程损失与局部损失的公式里都有 v²/2g，只是系数不同——λ(L/d) 与 ζ。这提示一种很实用的记账法：令某个局部障碍的 ζ 等值于「某段长度直管的 λ(L/d)」，即 ζ = λ·(Le/d)，解出：
Le = ζ·d/λ
把每个弯头、阀门都折算成 Le 米「等效直管」，加到实际管长上，整段管路就只需算一次沿程损失——这就是当量长度法。它特别适合初步设计时的快速估算，手册常直接给出「管件当量长度表」（以管径倍数计，湍流时 90° 弯头约 30~45 d、全开闸阀约 10 d）。两点注意：① 折算依赖 λ，当量长度不是管件的固有属性——同一弯头在层流管路（λ 大）和湍流管路（λ 小）里折出的 Le 相差可达十几倍；② 当量长度法与「ζ 相加法」只是两种记账方式，对同一个障碍只能用其中一种，重复计入会多算损失。
📐 例题 6 当量长度法：弯头折算成直管
d = 100 mm 的水管，λ = 0.025。管路中有 2 个 90° 弯头（ζ = 0.9）和 1 个全开闸阀（ζ = 0.2），管内流速 v = 1.5 m/s。用当量长度法把三个管件折算成等效直管长度，并求它们总共造成的水头损失。
思路 每个管件折算 Le = ζd/λ，求和得 ΣLe；再并入达西公式还原损失；最后用 ζ 相加法独立验证。
解 ① 每个弯头：Le = ζd/λ = 0.9×0.1/0.025 = 3.6 m；闸阀：Le = 0.2×0.1/0.025 = 0.8 m。 ② ΣLe = 2×3.6 + 0.8 = 8.0 m——三个管件相当于凭空多出 8 m 直管。 ③ 损失：hj = λ·(ΣLe/d)·(v²/2g) = 0.025×(8.0/0.1)×(1.5²/19.62) = 0.025×80×0.1147 = 0.229 m。 ④ 验证（ζ 相加法）：Σζ = 2×0.9 + 0.2 = 2.0；hj = 2.0×v²/2g = 2×0.1147 = 0.229 m ✓ 两条路必然一致。
答：折算总当量长度 8.0 m；这些管件共造成约 0.23 m 水头损失。 ⚠ 易错点：折算长度必须乘回同一个 λ 才能还原损失——管路各段 λ 不同（如变径管）时应逐段折算；当量长度法的价值在于把「杂项」并入直管一次算完，而不是算得更准。
管路的串联与并联：像电路一样算管网
真实管网是管段的组合，最基本的两种接法是串联与并联。它们与电路的串并联惊人地相似：把「水头」当「电压」、「流量」当「电流」，直觉就能直接搬过来用。
串联管路：各管段流量相同（连续性方程的要求），总损失 = 各段损失之和。相当于电阻串联：Q 相同、水头逐段跌落。由于损失正比于 1/d⁵，串联系统中「细而长」的管段贡献绝大部分损失——改粗哪一段最划算，一目了然。
并联管路：两节点之间并联几条支路，节点间的压差（即水头损失）只有唯一定值，所以各支路损失必然相等；总流量 = 各支路流量之和。流量分配自动调节到「损失相等」为止：阻抗大的支路分到的流量小。把每条支路的损失写成 h = S·Q²，其中 S 称管路阻抗，由达西公式可推出：
S = 8λL/(π²gd⁵)
则串联时 S总 = ΣSi（同一 Q 下损失相加）；并联时 1/√S总 = Σ(1/√Si)，且各支路流量比 Q₁ : Q₂ = √S₂ : √S₁——按阻抗平方根反比分配。类比也有边界：欧姆定律是线性的（U = RI），而湍流管路是平方律（h = SQ²）；只有层流管路 h ∝ Q 才真正线性、「像欧姆定律」。把类比用对场合，室内供水、暖通水系统、消防环网的分析就有了统一框架。
📐 例题 7 串联管路全链条：从吸水池到高位水箱的泵扬程
水泵把水从吸水池送进高位水箱，两液面高差 Δz = 25 m，流量 Q = 10 L/s。管路两段串联：L₁ = 120 m、d₁ = 150 mm、λ₁ = 0.02；L₂ = 80 m、d₂ = 100 mm、λ₂ = 0.022。局部损失：吸水口（带滤网，ζ = 0.5，在大管段）、90° 弯头 2 个（ζ = 0.9，在大管段）、d₁→d₂ 突缩（ζ = 0.28，按小管流速）、淹没出流出口（ζ = 1.0，按小管流速）。泵效率 η = 0.75。求泵扬程与轴功率。
思路 串联 → Q 相同 → 各段算流速 → 沿程 + 局部逐项求和得 hw → 两液面间列能量方程得 H = Δz + hw → 轴功率 P = ρgQH/η。
解 ① 流速：A₁ = π×0.15²/4 = 0.01767 m²，v₁ = 0.01/0.01767 = 0.566 m/s，v₁²/2g = 0.0163 m；A₂ = π×0.1²/4 = 0.007854 m²，v₂ = 0.01/0.007854 = 1.273 m/s，v₂²/2g = 0.0826 m。 ② 沿程：hf1 = λ₁(L₁/d₁)(v₁²/2g) = 0.02×800×0.0163 = 0.261 m；hf2 = 0.022×(80/0.1)×0.0826 = 0.022×800×0.0826 = 1.454 m。 ③ 局部（大管段三项，配 v₁）：hj1 = (0.5 + 2×0.9)×0.0163 = 2.3×0.0163 = 0.038 m；局部（小管段两项，配 v₂）：hj2 = (0.28 + 1.0)×0.0826 = 1.28×0.0826 = 0.106 m。 ④ 总损失：hw = 0.261 + 1.454 + 0.038 + 0.106 ≈ 1.86 m。 ⑤ 扬程（两液面均为大气压、液面速度≈0）：H = Δz + hw = 25 + 1.86 = 26.86 m。 ⑥ 轴功率：P = ρgQH/η = 1000×9.81×0.01×26.86/0.75 ≈ 3511 W ≈ 3.5 kW。
答：泵扬程约 26.9 m，轴功率约 3.5 kW。 ⚠ 易错点：① 串联各段流量相同但流速不同，每个 ζ 必须配它所在位置的流速；② 突缩 ζ 按下游小管流速计；③ 本例中 25 m 用于「抬升」、只有 1.86 m 是损失——高扬程系统里损失占比小，但长距离、低高差输水时损失反而占大头，两类场景都要心里有数。
📐 例题 8 并联管路：流量按阻抗分配
两节点间并联两条支路，d 同为 100 mm、λ 同为 0.02：支路 1 长 L₁ = 100 m，支路 2 长 L₂ = 25 m（局部损失忽略）。总流量 Q = 6 L/s。求两支路各自的流量与共同的水头损失。
思路 并联 → 损失相等：S₁Q₁² = S₂Q₂²；加总流量 Q₁ + Q₂ = Q；阻抗 S = 8λL/(π²gd⁵)。
解 ① 阻抗比：同 d、同 λ 时 S ∝ L，故 S₁/S₂ = L₁/L₂ = 4。 ② 损失相等 → Q₁/Q₂ = √(S₂/S₁) = √(1/4) = 1/2——长管只分到短管一半的流量。 ③ 分配：Q₁ + Q₂ = 6 L/s → Q₁ = 2 L/s，Q₂ = 4 L/s。 ④ 共同损失核算：S₁ = 8λL₁/(π²gd⁵) = 8×0.02×100/(9.870×9.81×10⁻⁵) ≈ 1.65×10⁴ s²/m⁵；h = S₁Q₁² = 1.65×10⁴×(0.002)² = 0.066 m。支路 2：S₂ ≈ 0.413×10⁴ s²/m⁵，S₂Q₂² = 0.413×10⁴×(0.004)² = 0.066 m ✓ 两支路损失严格相等。
答：Q₁ = 2 L/s、Q₂ = 4 L/s，共同水头损失约 0.066 m。 ⚠ 易错点：并联分配的是「损失相等」，不是「流量相等」——阻抗不同就自动按 √(S₂/S₁) 分流量；另外注意 d 的五次方主导：若支路 2 管径改成 75 mm（只细 25%），S₂ 会增大 (0.1/0.075)⁵ ≈ 4.2 倍，长而粗的支路 1 反而成为「主路」，流量分配立刻逆转。
工程视角：管路特性曲线与「工况点」
把例题 7 的思路一般化：对固定管路，所需扬程随流量按 H = Δz + S总Q² 变化（把全部沿程与局部损失折进一个系数 S总）——这条抛物线叫管路特性曲线。泵出厂时也带一条 H–Q 曲线（流量越大、能给的扬程越低）。两条曲线的交点就是工况点：管路「要多少」与泵「给得起多少」在此一拍即合。把阀门关小（S 增大），管路曲线上翘、交点左移、流量下降——这就是用阀门「调流量」的原理，代价是能量白白耗散在阀门里；变频调速则让泵曲线整体平移下移、交点沿管路曲线滑动，比节流阀节能得多。设计验收时若实测流量偏离设计值，第一件事就是把两条曲线画出来，看是管路阻抗变了还是泵磨损了。
🧠 一句话记忆 沿程损失 λ(L/d)(v²/2g) 管长摊；局部损失 ζ(v²/2g) 一处一处加。总损失 = 全部加起来。
物理根源 局部损失的本质也是粘性耗散：障碍物使流动产生涡流，涡流里流体剧烈旋转、相互摩擦，机械能被大量转成热。所以局部损失的「锅」最终还是要粘性来背——没有粘性，就不会有涡流耗散。
本章小结
- 水头损失分沿程损失 hf（管壁摩擦，沿管长均匀）和局部损失 hj（弯头阀门等局部障碍），总损失 hw = Σhf + Σhj。
- 达西公式 hf = λ·(L/d)·(v²/2g)，核心是确定阻力系数 λ。
- 层流 λ = 64/Re（精确，与粗糙度无关）；湍流 λ 依赖 Re 和相对粗糙度 ε/d，查莫迪图。
- 局部损失 hj = ζ·v²/2g，ζ 查表；其本质仍是粘性导致的涡流耗散。
- 管路计算 = 用「带损失的能量方程」把沿程和局部损失全部计入，求泵扬程、流速或管径。
- 当量长度法 Le = ζd/λ：局部损失可折算成等效直管长度——沿程与局部损失本质同源。
- 串联：流量相同、损失相加（S总=ΣSi）；并联：损失相等、流量按 √(S₂/S₁) 分配；阻抗 S = 8λL/(π²gd⁵)，管径占五次方。
读完本章，你应该能……
- 分清沿程损失和局部损失，并说出各自的来源。
- 默写达西公式，解释每一项的物理含义。
- 写出层流的 λ=64/Re，并说明为何层流与粗糙度无关。
- 说明湍流里 λ 由什么决定、如何用莫迪图。
- 用总损失公式做一个简单管路的损失估算。
- 会用当量长度法 Le = ζd/λ 把局部损失折算成直管长度。
- 会算串联管路（流量相同、损失相加）与并联管路（损失相等、流量按 √(S₂/S₁) 分配）。
- 能用「管路特性曲线 + 泵的 H–Q 曲线」解释工况点与阀门调流量的代价。
思考题
1. 达西公式 hf = λ(L/d)(v²/2g) 里，为什么损失正比于「管长」而不是「管径」？
提示：沿程损失来自壁面摩擦力沿整段管的累积，管越长、摩擦路径越长，损失自然越大；(L/d) 表示「相对管径有多长」。
2. 为什么层流的沿程阻力系数与管壁粗糙度无关，而湍流却强烈依赖粗糙度？
提示：层流分层滑动，贴壁层把粗糙凸起「抹平」、隔绝了主流；湍流剧烈掺混，粗糙凸起刺入流核、直接扰动主流，所以粗糙度影响显著。
3. 局部损失的本质是什么？它和沿程损失有什么共同点？
提示：本质都是粘性耗散——障碍物诱发涡流，涡流内剧烈摩擦把机械能转成热。共同点：都源于粘性，都表现为机械能（水头）不可逆损失。
4. 一段管道有两个 90° 弯头和一个阀门，总水头损失怎么写出来？
提示：hw = λ(L/d)(v²/2g)（直管沿程）+ ζ弯1(v²/2g) + ζ弯2(v²/2g) + ζ阀(v²/2g)，把沿程和各个局部损失逐项相加。
5. 当量长度法的思想是什么？为什么说它体现了「沿程损失与局部损失同质」？
提示：令 ζ(v²/2g) = λ(Le/d)(v²/2g)，得 Le = ζd/λ——把局部损失折算成等效直管长度。既然二者可以互相折算，说明本质上都是粘性耗散，只是空间分布不同（沿管均匀分布 vs 集中在局部）。
6. 两条支路并联，流量一定按各占一半分配吗？
提示：不对。并联支路损失必然相等，流量按 Q₁/Q₂ = √(S₂/S₁) 分配，只有阻抗相等时才各占一半；管径的影响是五次方的，远比管长的影响大。
7. 为什么突扩损失能用理论公式算出来，而弯头、阀门的 ζ 只能靠实验？
提示：突扩后的流动分离区结构简单，可近似「突扩断面压强均匀、壁面摩擦可忽略」，动量方程与能量方程联立即得 Borda–Carnot 解析式；弯头、阀门内部是复杂三维分离流，理论无从下手，只能实验测定后列成手册。