一个化工厂的管道系统,从泵出口到末端用户,总压降可能有几十米液柱。如果设计时算小了——泵选小了,流量不够,产线开不起来。算大了——泵选大了,白花电费,阀门常年憋压。一个年产 10 万吨的化工厂,泵的电费一年 300-500 万。泵的效率差 5 个百分点,一年就多花 15-25 万。
管道压降计算是工艺工程师的基本功。但说实话,很多工程师只会用软件点按钮,不明白公式里每个参数在物理上对应什么。这导致参数输错了也不知道,算出来的结果跟现场差 30% 也解释不通。
这篇文章把 Darcy-Weisbach 公式拆开讲透,配上实际算例,让你现场调泵的时候心里有数。
核心公式:Darcy-Weisbach
$$\Delta P = f \cdot \frac{L}{D} \cdot \frac{\rho v^2}{2}$$
或者用压头(液柱高度)表示:
$$h_f = f \cdot \frac{L}{D} \cdot \frac{v^2}{2g}$$
每个符号的含义:
| 符号 | 名称 | 单位 | 工程上的意义 |
|——|——|——|————|
| ΔP | 摩擦压降 | Pa | 这段管子吃掉的压力 |
| h_f | 摩擦压头损失 | m | 同一段管子,用液柱高度表示 |
| f | Darcy 摩擦因子 | 无量纲 | 管子的”粗糙程度”对流动的影响 |
| L | 管长 | m | 直管段长度(弯头阀门另算当量长度) |
| D | 管内径 | m | 注意是内径不是公称直径! |
| ρ | 流体密度 | kg/m³ | 水 1000,汽油 ~750,浓硫酸 ~1840 |
| v | 平均流速 | m/s | 流量除以截面积 |
| g | 重力加速度 | m/s² | 9.81 |
第一步:算摩擦因子 f —— Colebrook 方程
Darcy 摩擦因子 f 不是常数,它取决于两个东西:雷诺数 Re(流态)和相对粗糙度 ε/D(管壁有多糙)。
雷诺数:
$$Re = \frac{\rho \cdot v \cdot D}{\mu}$$
– Re < 2300:层流,f = 64/Re(简单公式)
– Re > 4000:湍流,f 用 Colebrook 方程(需要迭代算)
– 2300 < Re < 4000:过渡区,按最不利工况选大的
Colebrook 方程(湍流):
$$\frac{1}{\sqrt{f}} = -2.0 \log_{10}\left(\frac{\varepsilon/D}{3.7} + \frac{2.51}{Re\sqrt{f}}\right)$$
这个方程右边也有 f,所以需要迭代求解。工程上有两种做法:
1. Moody Chart(查图)——传统做法,精度一般但够用
2. Swamee-Jain 显式公式——不需要迭代,精度在 ±1%以内
Swamee-Jain 公式(推荐直接用这个):
$$f = \frac{0.25}{\left[\log_{10}\left(\frac{\varepsilon/D}{3.7} + \frac{5.74}{Re^{0.9}}\right)\right]^2}$$
适用范围:10⁻⁶ ≤ ε/D ≤ 10⁻²,5000 ≤ Re ≤ 10⁸。绝大多数工业管道都落在这个范围里。
常用管材的绝对粗糙度 ε:
| 管材 | ε (mm) | 相对粗糙度(DN100) |
|——|——–|——————-|
| 无缝钢管(新) | 0.02-0.05 | 0.0002-0.0005 |
| 焊接钢管(新) | 0.04-0.10 | 0.0004-0.0010 |
| 镀锌钢管 | 0.10-0.15 | 0.0010-0.0015 |
| 不锈钢管(新) | 0.015-0.03 | 0.00015-0.0003 |
| 铸铁管(新) | 0.25-0.50 | 0.0025-0.0050 |
| 铸铁管(旧/锈蚀) | 1.0-3.0 | 0.01-0.03 |
| PVC/PE 塑料管 | 0.001-0.005 | 0.00001-0.00005 |
| 玻璃钢管 | 0.005-0.01 | — |
| 混凝土管 | 0.3-3.0 | — |
注意:旧管道的粗糙度可以比新管大 5-20 倍。如果你在改造一个运行了 15 年的化工厂,按新管的 ε 算压降,现场实测会比你算的大 20-40%。这就是为什么老厂改造时经常发现泵”越开越不够”——不是泵老了,是管壁锈了。
第二步:弯头、阀门、三通——当量长度法
管道系统不止有直管,还有管件。管件的压降用”当量长度”法——每个管件相当于多长直管的摩擦损失。
当量长度 L_eq = K × (D/f),或者直接用经验值:
| 管件 | 当量长度(Le/D) | DN100 当量长度 (m) |
|——|—————–|——————-|
| 90° 长半径弯头 | 16-20 | 1.6-2.0 |
| 90° 短半径弯头 | 25-30 | 2.5-3.0 |
| 45° 弯头 | 10-12 | 1.0-1.2 |
| 三通(直通) | 18-22 | 1.8-2.2 |
| 三通(分支) | 50-65 | 5.0-6.5 |
| 全开闸阀 | 7-8 | 0.7-0.8 |
| 全开截止阀 | 300-340 | 30-34 |
| 全开球阀 | 3-5 | 0.3-0.5 |
| 全开蝶阀 | 20-30 | 2.0-3.0 |
| 止回阀(旋启式) | 50-135 | 5-13.5 |
| 全开调节阀 | 130-200 | 13-20 |
关键认知:一个截止阀的压降相当于 30 米直管。设计泵的出口管路时,泵出口的第一个阀——如果空间有限只能装截止阀而不是闸阀或球阀——等于在泵口上直接加了 30 米DN100的管子。所以泵出口优先用闸阀或球阀,不要用截止阀。
第三步:系统总压降 = 直管 + 管件 + 高程 + 速度头
一条完整管路的伯努利方程:
$$h_{total} = h_f(直管) + h_f(管件) + \Delta z + \frac{\Delta v^2}{2g}$$
– Δz:出口和进口的高程差(m)。往上打 10 米就是 10 米压头
– Δv²/2g:速度头变化。进出口管径一样则这一项等于 0,出口管径变小则速度头增加
工程上常被忽略的三项附加损失:
1. 入口损失:流体从储罐进入管道,突然收缩造成的损失。K_entrance ≈ 0.5(齐平入口)到 0.78(伸入入口)。
2. 出口损失:流体从管道排入储罐,突然扩大造成的损失。K_exit ≈ 1.0(所有的速度头都损失掉了)。
3. 过滤器/滤网的压降:泵入口的 Y 型过滤器或 T 型过滤器,清洁状态下压降约 1-3 kPa(0.1-0.3 m),堵塞后可以到 50-100 kPa。过滤器堵了导致泵汽蚀,是工业现场最常见的泵故障原因之一。
实战算例:一个循环水系统的压降计算
场景:循环水泵从水池吸水,送到 15 米高的冷却塔布水器,流量 100 m³/h,管径 DN150,管道总长 80 米(含来回)。
已知条件:
– 流体:水,30°C,ρ = 995 kg/m³,μ = 0.80 cP = 8.0×10⁻⁴ Pa·s
– 管道:DN150 无缝钢管(新),内径 154.1 mm(查管径表),ε = 0.045 mm
– 流量 Q = 100 m³/h = 0.0278 m³/s
– 管长 L = 80 m
– 高程差 Δz = 15 m(从水池液位到布水器出口)
管件清单:
– 90° 长半径弯头 × 6
– 闸阀 × 3(全开)
– 止回阀 × 1(旋启式)
– 三通(直通)× 2
– 入口(从水池)× 1
– 出口(到布水器)× 1
计算步骤
1. 流速:
v = Q / A = 0.0278 / (π × 0.1541² / 4) = 0.0278 / 0.01864 = 1.49 m/s
(在推荐流速范围内,DN150 循环水建议 1.2-2.0 m/s ✓)
2. 雷诺数:
Re = ρvD/μ = 995 × 1.49 × 0.1541 / (8.0×10⁻⁴) = 2.86×10⁵
湍流(>4000)
3. 相对粗糙度:
ε/D = 0.045 / 154.1 = 0.000292
4. 摩擦因子 f(Swamee-Jain):
f = 0.25 / [log₁₀(0.000292/3.7 + 5.74/Re^0.9)]²
= 0.25 / [log₁₀(0.0000789 + 5.74/1.26×10⁵)]²
= 0.25 / [log₁₀(0.0000789 + 0.0000456)]²
≈ 0.0162
5. 直管摩擦压头:
h_f,straight = f × (L/D) × (v²/2g)
= 0.0162 × (80/0.1541) × (1.49²/(2×9.81))
= 0.0162 × 519.2 × 0.1132
= 0.95 m
6. 管件当量长度:
| 管件 | Le/D | 数量 | 当量长度 (m) |
|——|——|——|————-|
| 90°弯头 | 18 | 6 | 18×6×0.1541 = 16.6 |
| 闸阀(全开) | 8 | 3 | 8×3×0.1541 = 3.7 |
| 止回阀 | 80 | 1 | 80×1×0.1541 = 12.3 |
| 三通(直通) | 20 | 2 | 20×2×0.1541 = 6.2 |
| 管件合计 | | | 38.8 m |
7. 管件摩擦压头:
h_f,fittings = f × (L_eq/D) × (v²/2g)
= 0.0162 × (38.8/0.1541) × 0.1132
= 0.46 m
8. 入口 + 出口损失:
h_entrance = 0.5 × v²/2g = 0.5 × 0.1132 = 0.06 m
h_exit = 1.0 × v²/2g = 1.0 × 0.1132 = 0.11 m
9. 总压头:
h_total = 0.95 + 0.46 + 0.06 + 0.11 + 15
= 16.58 m
≈ 16.6 m
10. 泵选型:
Q = 100 m³/h,H ≈ 18 m(加 10% 裕量)
选泵时 H_design = 1.1 × 16.6 ≈ 18.3 m,取整选 20 m 的泵。
现场调泵速查
如果泵打不上量,按这个顺序排查:
1. 过滤器堵了 → 拆开清洗,压降立刻下降(最常见,占 60%)
2. 泵入口有气 → 检查吸入管道的”U”形弯(高点集气),装排气阀
3. 泵汽蚀 → 听泵的声音——噼里啪啦像石子撞击的声音就是汽蚀。提高吸入液位或降低泵安装高度
4. 阀门没全开 → 尤其是泵入口阀(截止阀?换成闸阀!)
5. 管壁结垢/锈蚀 → 长期运行后压降逐年增加,考虑管道清洗或更换
6. 实际流量远超设计流量 → 泵的工作点偏移,出口压力不够。关小出口阀让泵回到设计点
推荐工具
– 手机上:Pipe Flow Wizard(iOS/Android),现场快速估算
– 电脑上:AFT Fathom(流体分析)、Pipe-Flo(全管网建模)、EPANET(给水管网,免费)
– Excel:自己建个压降计算表,把 Colebrook 和当量长度固化进去。一个下午的事,用十年
总结
Darcy-Weisbach 是管道压降计算的灵魂。公式本身不难,难在三个细节:
1. 用内径不是公称直径——DN150 的内径不是 150 mm
2. 管件压降不是零头——一条复杂管路的管件损失可以占直管损失的 30-50%
3. 旧管道的粗糙度是变量——运行十年后和新管完全不是一个东西
搞清这三点,你算的压降就能和现场对得上。对得上,泵就选得对。泵选得对,厂里的电费就对。