井筒摩阻计算原理与方法

1、井筒摩阻计算第一节水头损失及其分类实际流体具有粘性，在通道内流动时，流体内部流层之间存在相对运动和流动阻力。流动阻力和水头损失的规律，因流体的流动状态和流动的边界条件而异一、水头损失分类流体在流动的过程中，在流动的方向、壁面的粗糙程度、过流断面的形状和尺寸均不变的均匀流段上产生的流动阻力称之为沿程阻力，或称为摩擦阻力。沿程阻力的影响造成流体流动过程中能量的损失或水头损失（习惯上用单位重量流体的损失表示）。沿程阻力均匀地分布在整个均匀流段上，与管段的长度成正比，一般用hf表示。另一类阻力是发生在流动边界有急变的流场中，能量的损失主要集中在该流场及附近流场，这种集中发生的能量损失或阻力称为局部阻力

2、或局部损失，由局部阻力造成的水头损失称为局部水头损失。二、水头损失分类1.沿程阻力损失对于圆管:hfhf入lv24R2glv2d2g式中：l管长；R水力半径；d管径；v断面平均流速；g一重力加速度；沿程阻力系数，也称达西系数。一般由实验确定。式中的无量纲系数不是一个常数，它与流体的性质、管道的粗糙程度以及流速和流态有关，在大多数工程问题中，hf确实与v2成正比。此外，这样做可以把阻力损失和流速水头合并在一起，便于计算。2.局部阻力损失hj2g式中:局部阻力系数，一般由实验确定。整个管道的阻力损失，应该等于各管段的沿程损失和所有局部损失的总和。第二节粘性流体流动流态一、粘性流体流动流态当流速较小

3、时，沿程损失与流速一次方成正比，当流速较大时，沿程损失几乎与流速的平方成正比，如图所示，并且在这两个区域之间有一个不稳定区域。当阀门B慢慢打开，并打开颜色水阀门D中一条线状的颜色水。它与水流不相混合，如图此时管中的水流流速较小，可以看到玻璃管63(b)所示。从这一现象可以看出，在管中流速较小时，管中水流沿管轴方向呈层状流动,各层质点互不掺混，这种流动状态称为层流。当阀门B逐渐开大，管中的水流流速也相应增大。此时会发现，在流速增加到某一数值时，颜色水原直线的运动轨迹开始波动，线条逐渐变粗，如图63(c)所示。继续增加流速，则颜色水迅速与周围的清水混合，6-3(d)所示。这表明液体质点的运动轨迹不

4、规则，各层液体相互剧烈混合，产生随机的脉动，这种流动称为紊流。水流流速从小变大。沿程阻力曲线的走线为A-B-C一D。如图62所示。若实验时流速由大变小。则上述观察到的流动现象以相反的程序重演，但有紊流转变为'层流的流速Vc（下临界流速）要小于由层流转变为紊流的流速Vc（上临界流速）。如图62所示。沿径阻力曲线的走线为D-C-A。如图6-2所示。实验进一步表明，同一实验装置的临界流速是不固定的，随着流动的起始条件和实验条件不同，外界干扰程度不同，其上临界流速差异很大，但是，其下临流流速却基本不变。在实际工程中，扰动是普遍存在的，上临界流速没有实际意义，一般指的临界流速即指下临界流速。二、

5、流态的判别准则流态不仅与断面平均流速V有关系，而且与管径d、液体粘性、密度有关。即流态既反映管道中流体的特性，同时又反映管道的特性。将上述四个参数合成一无量纲数（无具体单位），称为雷诺数，用Re表示。vdpvdp式中：Re雷t若数,V流速，m/s,P流体密度，kg/m3,流体粘度，Pa.s;运动粘度，m2/s注，粘度单位1Pa.s=1N*s/m2=10P=103cp=1Kcps=1kg/（m*s）1N=1kg*m/s2对应于临界流速的雷诺数，称为临界雷诺数，通常用表示。大量实验表明，在不同的管道、不同的液体以及不同的外界条件下临界雷诺数不同。通常情况下，临界雷诺数总在2300附近，Rec=23

6、00【例61】有一直径d25mm的室内上水管，如管中流速v1.0m/s水温t10c。（1） .试判别管中水的流态；（2） .试求管内保持层流状态的最大流速为多少？解：（1）10c时，水的运动粘性系数v1.31106m2/'s,此时，管内雷诺数Ryd100.025191002300,故管中水流为紊流。e1.31106（2）保持层流的最大流速就是临界流速，Re空2300所以vc2300 1.31 10 60.0250.12 m s第三节沿程水头损失与切应力的关系一、均匀流动方程式沿程阻力（均匀流内部流层间的切应力）是造成沿程水头损失的直接原因。建立沿程水头损失与切应力的关系式，再找出切应力

7、的变化规律，就能解决沿程水头损失的计算问题。在圆管恒定流均匀流段上设11和22断面，如图所示。作用于流段上的外力：压力、壁面切应力重力相平衡。即：p1Ap2AAlcosw10式中w 壁面切应力,湿周，a圆管倾斜角，A圆管断面截面积，1管段长度。由几何关系得:1cosz1z2,除以A整理得：R7卫2w1zlz2AA并由断面1和断面2的能量方程得:E2z1z2hf,故:hfw1RhfR1A式中：R水力半径，R；J 水力坡度，Jhfl二、圆管过流段面上切应力分布rw0即圆管均匀过流断面上切应力呈直线分布，管轴处0,管壁处切应力达最大值三、壁剪切速度卜面在均匀流动方程式的基础上，推导沿程摩阻系数和壁面

8、切应力的关系。将JLY-代入均匀流动方程式，整理得:d2g，定义v 8具有速度的量纲，称为壁剪切速度（摩擦速度）。则:上式是沿程摩阻系数和壁面切应力的关系式，该式在紊流的研究中广为引用。第四节、圆管中的层流运动层流常见于很细的管道流动，或者低速、高粘流体的管道流动，如阻尼管、润滑油管、原油输油管道内的流动。圆管中层流运动的流动特征如前述，层流各流层质点互不掺混，对于圆管来说，各层质点沿平行管轴线方向运动。与管壁接触的一层速度为零，管轴线上速度最大，整个管流如同无数薄壁圆筒一个套着一个滑动。各流层间切应力服从牛顿内摩擦定律，即满足式du drdudyyror圆管层流的断面流动分布因讨论圆管层流运

9、动，所以可用牛顿内摩擦定律来表达液层间的切应力:dududydr式中为动力粘性，u为离管轴距离r处的切应力（即离管壁距离y处）的流速。对于均匀管流而言，在半径等于r处的切应力应为：rJ于均匀管流而言，根据式（621）,在半径等于r处的切应力应为:durdr 2积分得:Jr2C4利用管壁上的边界条件，确定上式中的积分常数Co当rro时u0,得：CJ.2J,22、uLr)66所示。过流断面上上式表明，圆管中均匀层流的流速分布是一个旋转抛物面，如图流速呈抛物面分布，这是圆管层流的重要特征之一。将r0代入上式，得到管轴处最大流速为u maxJr2 r0平均流速为:AudAr2 rdror21 01_r

10、02r J(r。2r 2)2 rdr三、圆管层流的沿程阻力损失将直径d代替式（634）中的20J()232进而可得水力坡度32以J h f / l代入上式,可得沿程阻力损失为:这就从理论上证明了圆管的均匀层流中.沿程阻力损失hf与平均流速v的一次方成正比，这与雷诺实验的结果相符。上式还可以进一步改写成达西公式的形式hf32d_264lv_264lvdvvdd2gRed2gl_v2d2g该式为达西和魏斯巴哈提出的著名公式， 雷诺数的函数，与管壁粗糙情况无关。例题6 2设有一恒定有压均匀管流.已知管径v 0.12 m/s ,水温t 10c时水的运动粘度由上式可得:64Re此公式表明圆管层流中的沿程

11、阻力系数只是d20mm,管长l20mm,管中水流流速，62v1.30610m/so求沿程阻力损失。解：Revd0.120.021.30610618382300为层流64Re6418380.0352一一2hf0.026mH2Olv20(0.12)0.035d2g0.0229.8第五节紊流运动分析一、粘性底层在紊流运动中，并不是整个流场都是紊流。由于流体具有粘滞性，紧贴管壁或槽壁的流体质点将贴附在固体边界上，无相对滑移，流速为零，继而它们又影响到邻近的流体速度也随之变小，从而在紧靠近面体边界的流层里有显著的流速梯度，粘滞切应力很大，但紊动则趋于零。各层质点不产生混掺，也就是说，在取近面体边界表面有

12、厚度极薄的层流层存在，称它为粘性底层或层流底层，如图68所示。在层流底层之外，还有一层很簿的过渡层。在此之外才是紊层，称为紊流核心区。在粘性底层中，速度按线性分布，在壁面上速度为零。粘性底层虽然很薄，但它对紊流的流速分布和流速阻力却有重大的影响。紊流紊流过流断面上流速成对数曲线分布，同层流过流断面上流速成抛物线分布相比，的流速分布均匀很多。第六节沿程水头损失系数的变化规律是计算沿程损失的关键。但由于紊流的复杂性，直到目前还不能像层流那样严格地从理论上推导出适合紊流的值来，所以值的确定，现有的方法仍然只有经验和半经验方法。一、阻力系数的影响因素在紊流中，除与反映流动状态的雷诺数有关之外，还因为突

13、入紊流核心的粗突起会直接影响流动的紊动程度，因而壁面粗糙度是影响阻力系数的另一个重要因素。可以采用一个指标即检验突起高度ks来表示壁面的粗糙程度，ks称为绝对粗糙度。绝对粗糙度具有长度量纲，所以仍感到有所不便，因而引入了量纲的相对粗糙度，即ks与直径(或半径)之比ks/d(或ks/r0),它是一个能够在不同直径的管流中用来反映管壁粗糙度的量，由以上分析可知，影响紊流沿程阻力系数的因素是雷诺数和相对粗糙度ks/d,写成函数关系式为f(Re,ks.d)二、尼古拉兹实验尼古拉兹在1933在人工粗糙管中系统地进行了沿程阻力系数和断面流速v的测定。他的实验涉及的参数范围比较大，相对粗糙度范围为ks/d1

14、3011014;雷诺数范围为Re500106纵坐标为lg(100),横坐标为lgRe,再算出和Re,取对数点绘在坐标IgRe去劣击工由血 工 L 一一 H . . A立 A3 4dA_dAS4a2.E2.8队0223.41638404一24.44.6481O5.25.45.65,Bb.口纸上，就得到f(Re,ks/d)曲线，即尼古拉兹曲线图。管道的流动可分为五个区域。第一个区域是层流区，对应的雷诺数Re2300(lgRe3.36)，试验点均落在直线ab上。表明与相对粗糙ks/d无关，只是Re的函数，并符合64/Re。还可知，沿程阻力损失hf与断面平均流速v成正比，这也与雷诺试验的结果一致。第二

15、个区域为层流与紊流之间的过渡区，Re23004000(lgRe3.363.6)试验点落在bc附近，表明与相对粗糙ks/d无关，只是Re的函数。此区是层流向紊流过渡，这个区的范围很窄，实用意义不大，不予讨论。第三个区域为紊流光滑区，Re4000(lgRe3.6)不同的相对粗糙管的试验点都先后落在同一条cd线上。表明与相对粗糙ks/d无关，是Re的函数。随着Re的增大，ks/d大的管道，实验点在Re较低时便离开此线，而ks/d较小的管道，在Re较大时才离开。第四个区域是紊流过渡区，不同的相对粗糙管实验点分别落在不同的曲线上。表明既与Re有关，又与ks/d有关。第五个区域是紊流粗糙区，不同的相对粗糙

16、管实验点分别落在不同水平直线上，表明与ks/d有关，与Re无关。在这个阻力区里，对于一定的管道(ks/d一定)，是常数。沿程水头损失与流速的平方成正比，故有称为阻力平方区。三、速度分布1 .紊流光滑区5.75lg”5.5v2 .紊流粗糙区uy半经验公式：5.75lg8.48vksU内指数公式：umaxr0四、的半经验公式1.光滑区沿程摩阻系数尼古拉兹光滑管公式21gRex2.512.粗糙区沿程摩阻系数Q.储q '. 6fL 005CL Q。，3.7dks21g尼古拉兹粗糙管公式工业管道实验曲线能甚公i£l聊则则MHISH5E!五、fl. XC1S ggs19 3 7 6 5

17、4 点 & ft* o o &uo V I i L & n4 Q o &- n- 0京味mIFr举赳4 j 5 « T IC, Hl51出粤s :?:2y*R3. -KCCE t. f1 W匚1ZC10MJ2< lfic >在尼古拉兹试验中，紊流有明显的光滑区。因为人工粗糙砂粒的直径是一致的。只要粘性底层的厚度大于砂粒直径，流动就处于光滑区。而工业管道、出于工业加工的缘故，不可能制造出粗糙度完全一致的管道。壁面的粗糙部分，从微观上讲，高低不一。因此没有明显的光滑区，或者光滑区的跨越范围很窄，无法进行对比。人工粗糙区，无沦是人工管道，还是工业

18、管道，由于粗糙面完全暴露在紊流中，其水头损失的变化规律也是一致的。因此，在相同的情况下。可用人工管道的相对粗糙度来表示工业管道的相对粗糙度，即当量粗糙度。当量粗糙度是用直径相同，在紊流粗糙区相同的人工管道的粗糙度ks,来定义该工业1939管道的粗糙度，表列出了常用工业管道的当量粗糙度。由表中数据可知，工业管道的计算方法与人工管道的计算方法一样。但尼古拉兹阻力系数公式在紊流过渡区是不适用的。年，柯列勃洛克和怀特给出了工业管道紊流区中的计算公式:211gks3.7d2.51Reks工业管道的当量粗糙度。常用工业管道的当量粗糙管道材料ksmm管道材料ksmm管道材料ksmm新氯乙烯00.002钢管0

19、.046新铸铁管0.150.5铅铜玻璃0.01涂沥青铸铁管0.12旧出铁管11.5镀锌钢管0.15混凝土管0.33.0为了将式(661)图形化，1944年，美国工程师穆迪以该公式为基础，以当量粗糙度为参数，用对数坐标绘制出工业管道摩阻损失系数曲线图，即穆迪图，见图611。六、沿程摩阻系数的经验公式布拉修斯(Blasius)公式1931年德国水力学家布拉修斯在总结前人实验的基础上总结并提出了紊流光滑区经验公式(6 62)0.3164C0.25Re该式形式简单，11算方便。在Re105范围内，有极高的精度，得到广泛的应用。希弗林松公式_.,ks、0.250.11()(663)d希弗林松粗糙区公式，

20、该式形式简单，计算方便，工程界经常采用。谢才公式和谢才系数22ghfhf.将达西公式(62)变换形式vd,以d4R,J,代入上式，整理得:llv.8gJrjc.rj(664)式中：v断面平均流速；R水力半径;J水力坡度；C谢才系数。上式最初是1769年法国工程师谢才直接根据渠道和塞纳河的实测资料提出的，是水力学最古老的公式之一，称为谢才公式。(665)式（665）给出了朗i才系数C和沿程摩阻系数的关系，谢才系数含有阻力的因素。流动阻力越大，谢才系数越小，反之亦然。1895年，爱尔兰工程师曼宁提出了计算谢才系数的经验公式：(6 66)C-R16n式中：n反映壁面粗糙性质并与流动性质无关的系数，称

21、为粗糙系数。Weymouth方程：0.0094073d为Darcy-Weisban摩阻系数，无量纲；D为管道内径，mWeymouth方程取管壁粗糙度k=0.0508mm（美国取k=0.02mm,苏联取k=0.03mm）,该方程适用于管径小、输量不大、净化程度较差的矿区集气管网。适用于管径范围254-508mm,阻力平方区Zegarola方程:114991.884*log10-Re*.该方程适用于3.2*10A4<Re<3.5*10A7Colebrook-White方程:在计算摩阻系数的方程中，CW方程最常用，是一个经验方程1k / D2 log 10( 12. 51Re *适用于各

22、个雷诺数，相对粗糙度0.00001-0.03之间，水力光滑区，混合摩擦区，阻力平方区。Gerg方程:2， 一log 10K nk/ D)n3. 71( f *1.499)0.942 * n* f 为Darcy-Weisban摩阻系数，无量纲；D为管道内径，m;f为阻力因子（等同于输送效率），无量纲；n为哥指数，描述从光滑管到粗糙管转变的变化剧烈程度；Re为雷诺数，无量纲。方程的特点：（1）对于完全粗糙流态（例如非常大的雷诺数），与Colebrook-White公式是一致的；（2）f也成为气流因子，描述完全管流中产生的额外损失；（3）对于水力光滑区（粗糙度k=0）f=1,公式退化为Zagarol

23、a方程，低雷诺数粗糙管流同样可以。显式:hioJ_k/D_1.499实例:计算数据见表L取GERG方程中参数门=1,/=1进行求解管道的实际运行数据当量粗糙度A加mRe管径门加是否有内涂层防护5082J2xW0.9664i序句什她果明炉儡重、廿*时*用«11式(1)1J7»mW9»7m4ieS1E-02QJ巾OBNARAB恚相耳法：2K出M2TMSnH22I-Q2J.0745北卜苏QWM731JO70»7K319»1E-«2J.49299E-09OMI7411414)31789941l<l-0Z轧倒ME-139Mnsi式闾U07

24、0)97099«1BSL依4.1«UL-1Sa.Wl766其他公式:计目公式适用危用Weyimjuih管径范国为254508nitn,阻力平方区Pan（A）相对粗糙度等于0.0000L雷诺数2翼Hf-1,4*1伊（水力光滑区）PtiiVB）相对粗相度0.00001.相诺数“阻力平方区）相对粗检度范围为00410050.0004之间，水力AGA光滑区相对粗糙度范围0,00005-0.01之间，咀力平方区Coicbrwk相对粗糙度范W0.000017）03C%水力光滑区，混合摩擦区，咀力平方区七、当量粗糙度当量粗糙度是将工业管道绝对粗糙度（e或者ks表示）折算成人工粗糙管的绝对

25、粗糙度。具体做法：是在流态处于完全粗糙管范围内，选用与实际管道摩阻系数（用符号f表示）值相等和管径相同的人工粗糙管的绝对粗糙度当做实际管道的当量粗糙度。表I管道当鼻相建度(后:值表皙世神美加工及使用狄加片mm1变化流用平均值驶璃.蟒.用臂斯的.犬格的-整体11剂的<J,0()1-0.01«0l»留营新的.光指的.整体打制的IJL0015-0.U6(1.113无缱讯管1,请的.清洁的.整件良好的2,用过几年后加以清洗的.跺西青后的，轻富转住的,污垢不零的0.0205n.i?-o.ft010.2岸接洞管却即展概管l口衽幡朝ar仔事网忤接的蝌句11;新的.播洁的1一轮请逸后

26、情1*不显看的/管”卜仲度精蚀的III曾，)中等第建的旧管2.大口粒的管1|酣1躺和怪继拈於锦111的13卜魏握坤梅悔修摄一"纲订门iftWt甲般M键密梅两打成青两转以上削打0,03-0.10.L-AL20.2*ft7d8L50,J*L0<LK1.2-ZRft解0,I?aji0l7L2Lft度酷捌管】.慷棒面光滑汕净的着看1短扑而一般的新管工用时几牢后的唱首0.07-0.10.L-<L20.4-Ui?OLis0.5铸帙百.新管1除钟育的世首工薪渐宵的1日管0.2-fll50.I-Ik150.12-03a3aib氾凝上背及擀好濯献上首L七林版面层你推机他1：百堂良轩.揖拜平

27、滑东模桢，施匚造发一股2. W4t出面层升若拄升1有骑浆面层1 1卜表而用蝌制刷过并蟀杆捆抹世(力表面用烟中解过.但未经技罡0,Q1.0-1.«0.35-l.fS0,7-XS>4.ad7L2ft7L28腺低就管a03卜'面惮个稠皮为例flJL除尔而L例=计鼻从井u开始第一井段*)u”距离测点的油管摩阻系数3/Zu 1> Fiu (int I 5 'I . /, = 501 > FTI . f>7= (1 & d= (k 162 nu27 3 Ttl K % r3(*7K卜(i'5-敝中菁15.KMTo：/>rr1.54：心一L54：/tL74C.LE.”一gmt2)HJ20X