第四节管路系统的计算-欢迎访问河南城建学院化学与化学工程

1.4 流体流动阻力
本节重点：直管阻力与局部阻力的计算，摩擦系数的影响因素。

难点：用因次分析法解决工程实际问题。

流动阻力的大小与流体本身的物理性质、流动状况及壁面的形状等因素有关。

化工管路系统主要由两部分组成，一部分是直管，另一部分是管件、阀门等。

相应流体流动阻力也分为两种：
直管阻力：流体流经一定直径的直管时由于内摩擦而产生的阻力；
局部阻力：流体流经管件、阀门等局部地方由于流速大小及方向的改变而引起的阻力。

1.4.1 流体在直管中的流动阻力
1. 阻力的表现形式
如图1-24所示，流体在水平等径直管中作定态流动。

在1-1′和2-2′截面间列柏努利方程，
f W p u
g z p u g z +++=++ρ
ρ222212112121 因是直径相同的水平管，21u u = 21z z =
ρ2
1p p W f -=∴ （1-34）
若管道为倾斜管，则
)()(22
11
g z p g z p W f +-+=ρρ （1-34a ）
由此可见，无论是水平安装，还是倾斜安装，流体的流动阻力均表现为静压能的减少，仅当水平安装时，流动阻力恰好等于两截面的静压能之差。

2. 直管阻力的通式
在图1-24中，对1-1′和2-2′截面间流体进行受力分析：
由压力差而产生的推动力为 ()4221d p p π- 与流体流动方向相同
流体的摩擦力为 dl A F τπτ== 与流体流动方向相反。

流体在管内作定态流动，在流动方向上所受合力必定为零。

dl d p p τππ=-4)(2
21
整理得 τd
l p p 421=- （1-35） 将式(1-35)代入式（1-34）中，得
τρ
d l W f 4= （1-36） 将式（1-36）变形，把能量损失f W 表示为动能2
2
u 的某一倍数。

2
82
2u d l u W f ρτ= 令 28u
ρτλ= 则 2
2
u d l W f λ= （1-37） 式（1-37）为流体在直管内流动阻力的通式，称为范宁（Fanning ）公式。

式中λ为无因次系数，称为摩擦系数或摩擦因数，与流体流动的Re 及管壁状况有关。

根据柏努利方程的其它形式，也可写出相应的范宁公式表示式：
压头损失 g
u d l h f 22
λ= （1-37a ） 压力损失 2
2
u d l p f ρλ=∆ （1-37b ） 值得注意的是，压力损失f p ∆是流体流动能量损失的一种表示形式，与两截面间的压力差)(21p p p -=∆意义不同，只有当管路为水平时，二者才相等。

应当指出，范宁公式对层流与湍流均适用，只是两种情况下摩擦系数λ不同。

以下对层流与湍流时摩擦系数λ分别讨论。

3. 层流时的摩擦系数
流体在直管中作层流流动时，管中心最大速度如式（1-35）所示。

将平均速度max 21u u =及2
d R =代入上式中，可得 22132)(d
lu p p μ=- 232d lu p f μ=
∆ （1-38） 式（1-38）称为哈根-泊谡叶（Hagen-Poiseuille ）方程，是流体在直管内作层流流动时压力损失的计算式。

结合式（1-34），流体在直管内层流流动时能量损失或阻力的计算式为
2
32d lu W f ρμ= （1-39） 表明层流时阻力与速度的一次方成正比。

式（1-39）也可改写为
2Re 64264322
22u d l u d l u d d
lu W f ⋅⋅=⋅⋅==ρμρμ （1-39a ） 将式（1-39a ）与式（1-37）比较，可得层流时摩擦系数的计算式
Re
64=λ （1-40） 即层流时摩擦系数λ是雷诺数Re 的函数。

4．湍流时的摩擦系数
（1）因次分析法
层流时阻力的计算式是根据理论推导所得，湍流时由于情况要复杂得多，目前尚不能得到理论计算式，但通过实验研究，可获得经验关系式，这种实验研究方法是化工中常用的方法。

在实验时，每次只能改变一个变量，而将其它变量固定，如过程涉及的变量很多，工作量必然很大，而且将实验结果关联成形式简单便于应用的公式也很困难。

若采用化工中常用的工程研究方法——因次分析法，可将几个变量组合成一个无因次数群（如雷诺数Re 即是由d 、ρ、u 、μ四个变量组成的无因次数群），用无因次数群代替个别的变量进行实验，由于数群的数目总是比变量的数目少，就可以大大减少实验的次数，关联数据的工作也会有所简化，而且可将在实验室规模的小设备中用某种物料实验所得的结果应用到其它物料及实际的化工
设备中去。

因次分析法的基础是因次一致性原则，即每一个物理方程式的两边不仅数值相等，而且每一项都应具有相同的因次。

因次分析法的基本定理是白金汉（Buckinghan ）的π定理：设影响某一物理现象的独立变量数为n 个，这些变量的基本因次数为m 个，则该物理现象可用N ＝（n －m ）个独立的无因次数群表示。

根据对摩擦阻力性质的理解和实验研究的综合分析，认为流体在湍流流动时，由于内摩擦力而产生的压力损失f p ∆与流体的密度ρ、粘度μ、平均速度u 、管径d 、管长l 及管壁的粗糙度ε有关，即
()εμρ,,,,,l d u f p f =∆ （1-41）
7个变量的因次分别为：
[p]=M θ-2L -1 [ρ]=ML -3 [u]=M θ-1
[d]=L [l]=L [ε]=L [μ]=M θ-1L -1
基本因次有3个。

根据π定理，无因次数群的数目
N=n-m=7-3=4个
将式（1-41）写成幂函数的形式：
f e d c b a f u l kd p εμρ=∆
因次关系式：
f e d c b a L ML ML L L L L M )1()()(13112-=-----θθθ
根据因次一致性原则：
对于M ： e d ＋=1
对于L ： f e d c b a +--+=-31＋
对于θ： e c --=-2
设b ，e ，f 已知，解得：
e
d e c f
e b a -=-=---=12
f e e e b f e b f u l kd p εμρ-----=∆12
f e b f d u d d l k u p ⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛⎪⎭
⎫ ⎝⎛=∆-εμρρ2 即 ⎪⎪⎭
⎫ ⎝⎛=∆d d l u d u p f εμρφρ,,2 （1-42） 式中 μρu
d ——雷诺数R
e ，
2u p f
ρ∆——欧拉（Euler ）准数，也是无因次数群。

d l 、d
ε均为简单的无因次比值，前者反映了管子的几何尺寸对流动阻力的影响，后者称为相对粗糙度，反映了管壁粗糙度对流动阻力的影响。

式（1-42）具体的函数关系通常由实验确定。

根据实验可知，流体流动阻力与管长l 成正比，该式可改写为：
⎪⎭⎫ ⎝⎛=∆d d l u p f
εφρRe,2
（1-43） 或 2Re,u d d l p W f f ⎪⎭⎫ ⎝⎛=∆=εψρ
（1-43a ） 与范宁公式（1-37）相对照，可得
)(Re,d
εψλ= （1-44） 即湍流时摩擦系数λ是Re 和相对粗糙度
d
ε的函数，如图1-25所示，称为莫狄（Moody ）摩擦系数图。

根据Re 不同，图1-25可分为四个区域；
（1）层流区 （Re ≤2000），λ与d ε无关，与Re 为直线关系，即Re
64=
λ，此时u W f ∝，即f W 与u 的一次方成正比。

（2）过渡区（2000<Re<4000），在此区域内层流或湍流的λ～Re 曲线均可应用，对于阻力计算，宁可估计大一些，一般将湍流时的曲线延伸，以查取λ值。

（3）湍流区（Re ≥4000以及虚线以下的区域），此时λ与Re 、d ε都有关，当 d ε一定时，λ随Re 的增大而减小，Re 增大至某一数值后，λ下降缓慢；当R e 一定时，λ随d ε的增加而增大。

（4）完全湍流区 （虚线以上的区域），此区域内各曲线都趋近于水平线，即λ与Re 无
关，只与d ε有关。

对于特定管路d ε一定，λ为常数，根据直管阻力通式可知，2u W f ∝，所以此区域又称为阻力平方区。

从图中也可以看出，相对粗糙度d ε愈大，达到阻力平方区的Re 值愈低。

对于湍流时的摩擦系数λ，除了用Moody 图查取外，还可以利用一些经验公式计算。

这里介绍适用于光滑管的柏拉修斯（Blasius ）式：
25
.0Re 3164.0=λ （1-45） 其适用范围为Re ＝5×103～105。

此时能量损失f W 约与速度u 的1.75次方成正比。

图1-25 摩擦系数λ与雷诺数Re 及相对粗糙度d ε的关系
考莱布鲁克（Colebrook ）式
此式适用于湍流区的光滑管与粗糙管直至完全湍流区。

（2）管壁粗糙度对摩擦系数的影响
光滑管：玻璃管、铜管、铅管及塑料管等称为光滑管；
粗糙管：钢管、铸铁管等。

管道壁面凸出部分的平均高度，称为绝对粗糙度，以ε表示。

绝对粗糙度与管径的比值即d
ε，称为相对粗糙度。

工业管道的绝对粗糙度数值见教材。

管壁粗糙度对流动阻力或摩擦系数的影响，主要是由于流体在管道中流动时，流体质点与管壁凸出部分相碰撞而增加了流体的能量损失，其影响程度与管径的大小有关，因此在摩擦系数图中用相对粗糙度d ε，而不是绝对粗糙度ε。

流体作层流流动时，流体层平行于管轴流动，层流层掩盖了管壁的粗糙面，同时流体的流动速度也比较缓慢，对管壁凸出部分没有什么碰撞作用，所以层流时的流动阻力或摩擦系数与管壁粗糙度无关，只与Re 有关。

流体作湍流流动时，靠近壁面处总是存在着层流内层。

如果层流内层的厚度δL 大于管壁的绝对粗糙度ε，即δL >ε时，如图1-26（a ）所示，此时管壁粗糙度对流动阻力的影响与层流时相近，此为水力光滑管。

随Re 的增加，层流内层的厚度逐渐减薄，当δL <ε时，如图1-26（b ）所示，壁面凸出部分伸入湍流主体区，与流体质点发生碰撞，使流动阻力增加。

当Re 大到一定程度时，层流内层可薄得足以使壁面凸出部分都伸到湍流主体中，质点碰撞加剧，致使粘性力不再起作用，而包括粘度μ在内的Re 不再影响摩擦系数的大小，流动进入了完全湍流区，此为完全湍流粗糙管。

图1-26 流体流过管壁面的情况
⎪⎪⎭⎫ ⎝
⎛+-=λελRe 7.182log 274.11d
例 分别计算下列情况下，流体流过φ76×3mm 、长10m 的水平钢管的能量损失、压头损失及压力损失。

（1） 密度为910kg/、粘度为72cP 的油品，流速为1.1m/s ；
（2） 20℃的水，流速为2.2 m/s 。

解：（1）油品：
200097310
721.191007.0Re 3<=⨯⨯⨯==-μρu
d 流动为层流。

摩擦系数可从图1-25上查取，也可用式（1-40）计算：
0658.0973
64Re 64===λ 所以能量损失 J/kg 69.52
1.107.0100658.022
2===u d l W f λ 压头损失 m 58.081
.969.5===g W h f
f 压力损失 Pa 517869.5910=⨯==∆f f W p ρ
（2）20℃水的物性：3kg/m 2.998=ρ，310005.1-⨯=μ Pa ·s
531053.110
005.12.22.99807.0⨯=⨯⨯⨯==-μρu
d R
e 流动为湍流。

求摩擦系数尚需知道相对粗糙度d ε，查表1-2，取钢管的绝对粗糙度ε为0.2mm ，则 00286.070
2.0==d ε
根据Re =1.53×105及d ε＝0.00286查图1-25，得λ＝0.027
所以能量损失 J/kg 33.92
2.207.010027.022
2===u d l W f λ 压头损失 m 95.081
.933.9===g W h f
f 压力损失 Pa 931333.92.998=⨯==∆f f W p ρ
5. 非圆形管道的流动阻力
对于非圆形管内的湍流流动，仍可用在圆形管内流动阻力的计算式，但需用非圆形管道的当量直径代替圆管直径。

当量直径定义为
∏
⨯⨯=A d e 44＝润湿周边流通截面积 (1-46) 对于套管环隙，当内管的外径为d 1，外管的内径为d 2时，其当量直径为
()
1212212
244d d d d d d d e -=+-=πππ
对于边长分别为a 、b 的矩形管，其当量直径为
b
a a
b b a ab d e +=+=2)(24 在层流情况下，当采用当量直径计算阻力时，还应对式（1-40）进行修正，改写为 Re
C =
λ （1-47） 式中C 为无因次常数。

一些非圆形管的C 值见教材。

注意，当量直径只用于非圆形管道流动阻力的计算，而不能用于流通面积及流速的计算。

1.4.2 局部阻力
局部阻力有两种计算方法：阻力系数法和当量长度法。

1. 阻力系数法
克服局部阻力所消耗的机械能，可以表示为动能的某一倍数，即 22
'
u W f
ζ= （1-48） 或 g u h f
22
'
ζ= （1-48a ） 式中ζ称为局部阻力系数，一般由实验测定。

常用管件及阀门的局部阻力系数见教材。

注意表中当管截面突然扩大和突然缩小时，式（1-48）及（1-48a ）中的速度u 均以小管中的速度计。

当流体自容器进入管内，5.0=进口ζ，称为进口阻力系数；当流体自管子进入容器或从
管子排放到管外空间，1=出口ζ，称为出口阻力系数。

当流体从管子直接排放到管外空间时，管出口内侧截面上的压强可取为与管外空间相同，但出口截面上的动能及出口阻力应与截面选取相匹配。

若截面取管出口内侧，则表示流体并未离开管路，此时截面上仍有动能，系统的总能量损失不包含出口阻力；若截面取管出口外侧，则表示流体已经离开管路，此时截面上动能为零，而系统的总能量损失中应包含出口阻力。

由于出口阻力系数1=出口ζ，两种选取截面方法计算结果相同。

2．当量长度法
将流体流过管件或阀门的局部阻力，折合成直径相同、长度为e l 的直管所产生的阻力即
22
'
u d l W e f
λ= （1-49） 或 g u d l h e f
22
'
λ= （1-49a ） 式中e l 称为管件或阀门的当量长度。

同样，管件与阀门的当量长度也是由实验测定，有时也以管道直径的倍数d l e 表示。

见教材。

1.4.3 流体在管路中的总阻力
前已说明，化工管路系统是由直管和管件、阀门等构成，因此流体流经管路的总阻力应是直管阻力和所有局部阻力之和。

计算局部阻力时，可用局部阻力系数法，亦可用当量长度法。

对同一管件，可用任一种计算，但不能用两种方法重复计算。

当管路直径相同时，总阻力：
22'
u d l W W W f f f ⎪⎭⎫ ⎝⎛∑+=+=∑ζλ （1-55） 或 22
'
u d l l W W W e f f f ∑+=+=∑λ （1-55a ） 式中ζ∑、e l ∑分别为管路中所有局部阻力系数和当量长度之和。

若管路由若干直径不同的管段组成时，各段应分别计算，再加和。

例 如附图所示，料液由敞口高位槽流入精馏塔中。

塔内进料处的压力为30kPa （表压），输送管路为φ45×2.5mm 的无缝钢管，直管长为10m 。

管路中装有180º回弯头一个，90º标准弯头一个，标准截止阀（全开）一个。

若维持进料量为5m 3
/h ，问高位槽中的液面至少高出进料口多少米？
操作条件下料液的物性：3kg/m 890=ρ，3102.1-⨯=μ Pa ·s
解：如图取高位槽中液面为1-1′面，管出口内侧为2-2′截面，且以过2-2′截面中心线的水平面为基准面。

在1-1′与2-2′截面间列柏努利方程： f W p u g z p u g z ∑+++=++ρ
ρ222212112121 其中： z 1=h ； u 1≈0； p 1=0（表压）；
z 2=0； p 2=20kPa （表压）；
1.104
.0785.0360054222=⨯==d V u s
π
m/s 管路总阻力 22'
u d l W W W f f f ⎪⎭⎫ ⎝⎛∑+=+=∑ζλ 431001.310
3.11.18900
4.0Re ⨯=⨯⨯⨯==-μρu
d 取管壁绝对粗糙度3.0=εmm ，则 0075.040
3.0==d ε
从图1-25中查得摩擦系数036.0=λ
由表1-4查得各管件的局部阻力系数：
进口突然缩小 5.0=ζ
180º回弯头 5.1=ζ
90º标准弯头 75.0=ζ
标准截止阀（全开） 4.6=ζ
15.94.675.05.15.0=+++=∑∴ζ
98.102
1.115.904.010036.0222=⎪⎭⎫ ⎝⎛+⨯=⎪⎭⎫ ⎝⎛∑+=∑u d l W f ζλJ/kg 所求位差
m 62.481.9/)98.1021.18901030(/)2(2
3222
=++⨯=∑++=g W u p h f ρ 本题也可将截面2-2′取在管出口外侧，此时流体流入塔内，2-2′截面速度为零，无动能项，但应计入出口突然扩大阻力，又1=出口ζ，所以两种方法的结果相同。