CN104331603A - 一种任意组合弯管阻力系数的计算方法 - Google Patents
一种任意组合弯管阻力系数的计算方法 Download PDFInfo
- Publication number
- CN104331603A CN104331603A CN201410539203.8A CN201410539203A CN104331603A CN 104331603 A CN104331603 A CN 104331603A CN 201410539203 A CN201410539203 A CN 201410539203A CN 104331603 A CN104331603 A CN 104331603A
- Authority
- CN
- China
- Prior art keywords
- curved
- coefficient
- zeta
- resistance
- bend pipe
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Rigid Pipes And Flexible Pipes (AREA)
Abstract
本发明公开了一种组合弯管阻力系数计算方法,其特征是通过参数修正的方法获得每一个弯头及连接段直管的阻力系数,然后对各个部分的阻力系数进行叠加,从而计算出整个组合弯管的阻力系数。本发明能对不同相对曲率半径和连接段长度的组合弯管阻力系数值进行计算,结果更加精确可靠。
Description
技术领域
本发明涉及弯曲管道阻力系数的计算方法,具体地说是一种同一个平面内90°任意组合弯管阻力系数的计算方法。
背景技术
流体在弯管内的流动作为一种常见现象,广泛应用于如发动机系统、暖气装置、通风设备、煤气和天然气输送、热电厂及核电冷却水循环系统等工业设备中。阻力系数是管道设计中极为重要的水力参数,正确掌握其规律有助于提高设计的合理性。
由量纲分析原理可知,对圆截面管中的不可压缩流而言,形变件造成的流体阻力系数为:
ζ=△P/(1/2ρ·V2)=ζ(δ,R/D,Re) (1)
式中,△P为流体的总压压力降[Pa];ρ为流体密度[kg/m3];V为流体速度[m/s];δ为转角[°];R/D为相对曲率半径;R为弯管轴线的转弯半径[m];D为管径[m];Re为雷诺数,Re=V·D·ρ/μ,μ为流体动力粘滞系数[Pa·s]。
已发表的文献主要针对单个弯管、S型弯管、U型弯管和一些几何参数确定的组合弯管阻力系数的研究,其中最为经典的研究结果来自前苏联学者Idelichik,《Handbook ofhydraulic resistance(3rd Edition)》,2008,其认为弯管的总阻力系数由转弯区域内的摩擦阻力系数ζTP和局部阻力系数ζM之和来确定:
ζ=ζTP+ζM (2)
式中,ζTP与直管段计算公式相同,为:ζTP=λl/D,其中λ与Re和相对粗糙度有关,针对2×104≤Re≤105且管壁光滑的圆截面管子,λ值按下列公式计算:
紊流状态(4000<Re<105)
见图1中的a图。
紊流状态(Re>105)
此外,λ值也可按照下表一选取
表一:稳定流动;管壁光滑的圆截面管子摩擦阻力系数表
l为弯管的轴线长度,l/D的值按下式计算:
因此:
针对不同弯管类型,局部阻力系数ζM的值不同,使得阻力系数ζ也不同:
1、单弯管阻力系数(图2),光滑壁(△=0)并且Re=V·D·ρ/μ≥2×105
式中
A1——考虑弯管转角δ影响的系数;
B1——考虑弯管相对曲率半径R/D影响的系数;
C1——考虑弯管相对高宽比影响的系数。
A1按照Nekrasov,Hydraulics,1954,290p的数据选取:
当δ=90°时,A1=1;
当δ<70°时,A1=0.9sinδ;
当δ>100°时,
B1按照近似公式选取:
当R/D<1时,B1=0.21/(R/D)0.25;
当R/D≥1时,B1=0.21/(R/D)0.5,
C1可按照图表选取,本例中只考虑圆截面管子,因此C1=1。
此外,考虑到Re数的影响:
kRe=f(Re),见图3或者下表二:
表二:弯管(0<δ≤180°;l0/D≥10和任意R/D)kRe值
2、S型弯管(图4),流体在同一平面内流动,光顺的(R/D>1.0):
式中ζM′为式(6)中单个弯管的ζM值;
见图5或者下表三(当Re≥2×104时适用)。
式中:lk为连接段直管长度。
表三:R/D>1.0的S型弯管的局部阻力系数修正系数A值
3、U型弯管(图6),在同一平面内;光顺的(R/D≥1.0);0≤δ≤180°)
式中ζM′为式(6)中单个弯管的ζM值;
见图7或者下表四(当Re≥2×104时适用)。
表四:R/D>1.0的U型弯管的局部阻力系数修正系数A值
4、特定组合弯管(图8),在同一平面内;光顺的(R/D≥1.0);0≤δ≤180°
式中ζ1为式(7)的ζ;
ζTP′=λlk′/D;λ见图1或表一;当λ=0.02时,ζTP=0.02lk′/D;
lk′见图8。
由上可知Idelichik.Handbook of hydraulic resistance(3rd Edition)[B].New York:JaicoPublishing House,2008,虽然给出一些几何参数确定的特定组合弯管阻力系数的计算方法,但是对于任意组合弯管的阻力系数计算并未提出具体实用的计算方法。
此外,为了验证本发明的正确性,还采用CFD技术对弯管阻力系数值进行了数值分析。图9为用CFD数值计算得到的R/D=1的S型弯管阻力系数随连接段长度的变化,和Idelichik,《Handbook of hydraulic resistance(3rd Edition)》,2008中得到的阻力系数对比,可以看出两者之间吻合很好,特别是lk/D≥4以后两者的差距基本都在5%之内,证明了数值计算的准确性。
发明内容
本发明是为了克服现有技术的不足之处,提供一种具有两个弯以上的任意组合弯管阻力系数的计算方法,通过参数修正的方法确定各弯头的局部阻力系数,进一步使总阻力系数的计算结果更加精确可靠。
本发明采用的技术方案是:一种任意组合弯管阻力系数的计算方法,基于同一平面的90°单弯管、两个90°单弯管加之间的连接直管构成的S型弯管以及U型弯管的阻力系数计算公式,其特征在于,对于具有两个以上同一个平面内的90°弯及其连接直管构成的任意组合弯管,lkj/D≥4,lkj为相邻两个弯之间的连接段长度,D为管径;相对曲率半径1≤R/D≤3,R为弯管轴线的转弯半径;雷诺数2×104≤Re≤105且管壁光滑的圆截面管,从流体入口端到出口端通过参数修正的方法获得各弯头及连接段直管的阻力系数,然后对各个部分的阻力系数进行叠加,从而计算出整个组合弯管的总阻力系数:
为第i弯的单弯管局部阻力系数,采用单弯管局部阻力系数值式(6)计算,即:
为第i弯管弯曲部分的摩擦阻力系数;
为连接段的总摩擦阻力系数,lkj为第j个连接段;
Ai为第i弯的局部阻力系数的修正系数;
计算中需要确定的值只有局部阻力系数的修正系数Ai值,流体入口第二个弯以后的局部阻力系数Ai只与前一个弯和连接段长度有关,取值原则如下:
1)当连接段长度lkj/D>25时,该连接段之后那个弯管的局部阻力采用流体入口后第一个弯的计算公式计算;
2)第一个弯的局部阻力系数修正系数A1与相对曲率半径有关:
R/D=1时,A1=0.6AS/U;
R/D=2时,A1=0.7AS/U;
R/D≥2.5时,A1=0.5AS/U。
其中AS/U为由第一个弯和第二个弯组成的S/U型弯管的局部阻力系数修正系数A值,在Idelichik,《Handbook of hydraulic resistance(3rd Edition)》,2008中查得,即图6、图8或者表三、表四,其值是lkj/D的函数,若两弯的相对曲率半径不同,按照第一个弯的相对曲率半径计算;
3)第二个弯及其后面的弯管局部阻力系数的修正系数Ai值如下:
R/D=1时,Ai=0.4AS/U;
R/D=2时,Ai=0.3AS/U;
R/D≥2.5时,Ai=0.5AS/U;
其中AS/U为第i弯及其前一个弯组成的S/U型弯管的局部阻力系数修正系数A值,在Idelichik,《Handbook of hydraulic resistance(3rd Edition)》,2008中查得,即图6、图8或者表三、表四。若两弯的相对曲率半径不同,按照第i弯的相对曲率半径计算。
与已有技术相比,本发明的优点及有益效果体现在:
1、计算任意组合弯管阻力系数时考虑了连接段长度的影响以及邻弯之间曲率半径不同,适用范围更广。
2、本发明认为弯头局部阻力系数值的值只与前弯及其相邻的连接段长度相关,计算方法更加简便、准确。
3、对任意组合弯管每个弯头的局部阻力系数值都进行了合理的修正,从而使对弯头和直管的阻力系数进行叠加后的总阻力系数更加精确。
附图说明
图1为稳定流动;管壁光滑的圆截面管子摩擦阻力系数;
图2为单弯管示意图;
图3为Re数对弯管阻力系数的影响规律;
图4为S型弯管示意图;
图5为R/D>1.0的S型弯管的局部阻力系数修正系数A值;
图6为U型弯管示意图;
图7为R/D>1.0的U型弯管的局部阻力系数修正系数A值;
图8为特定组合弯管示意图;
图9为用CFD数值计算得到的R/D=1的S型弯管阻力系数随连接段长度的变化,与Idelichik,《Handbook of hydraulic resistance(3rd Edition)》,2008中得到的阻力系数对比,可以看出两者之间吻合很好,特别是lk/D≥4以后两者的差距基本都在5%之内,证明了数值计算的准确性;
图10为数值计算的三弯管阻力系数和采用本文提出的公式值之间的对比,可以看出数值计算值和公式值非常接近,均在3%之内;
图11某一平面90°任意组合多弯管,其中转弯的方向、相对曲率半径以及连接段长度都在计算方法所适用的参数范围内任意选取;
具体实施方式
如图11所示,是一个具有9个弯和8个连接段长度的任意组合弯管,为了直观仅对其轴线进行显示。则从左端入口到出口的分别为第1,2…i…9弯,弯头编号与相对曲率半径的对应关系如下表所示:
第一个弯和第二个弯之间直管段定义为第一个连接段长度,则从左端入口到出口的分别为第1,2…j…8个连接段长度,连接段编号与相对连接段长度对应关系如下表所示:
假设弯管内的流体流动已经不再受Re数影响,即Re≥2×105,则其阻力系数计算方法按如下过程进行:
其中ζ为总阻力系数;
为第i弯的单弯管局部阻力系数,可采用Idelichik,《Handbook of hydraulicresistance(3rd Edition)》,2008中的单弯管局部阻力系数值(式6);
为第i弯管弯曲部分的摩擦阻力系数;
为连接段的总摩擦阻力系数,lkj为第j个连接段;
Ai为第i弯的局部阻力系数的修正系数。
按照本发明的方法对图11的组合弯管的阻力系数进行计算如下(未出现相对连接段长度为达到lkj/D>25情况,故无需进行截断):
第一个弯和第二个弯以及其中间的连接段直管构成了一个S型弯管,且第一个弯的相对曲率半径为R/D=1,两个弯之间的相对连接段长度lk1/D=4,则第一个弯的局部阻力系数计算如下:
A1=0.6*AS=0.6*1.38=0.828,
摩擦阻力系数计算如下:
第二个弯和第一个弯以及其中间的连接段直管构成了一个S型弯管,且第二个弯的相对曲率半径为R/D=2,两个弯之间的相对连接段长度lk1/D=4,则第二个弯的局部阻力系数计算如下:
A2=0.3*AS=0.3*1.38=0.414,
摩擦阻力系数计算如下:
第三个弯和第二个弯以及其中间的连接段直管构成了一个S型弯管,且第三个弯的相对曲率半径为R/D=1,两个弯之间的相对连接段长度lk2/D=6,则第三个弯的局部阻力系数计算如下:
A3=0.3*AS=0.4*1.47=0.588,
摩擦阻力系数计算如下:
第四个弯和第三个弯以及其中间的连接段直管构成了一个S型弯管,且第四个弯的相对曲率半径为R/D=3,两个弯之间的相对连接段长度lk3/D=8,则第四个弯的局部阻力系数计算如下:
A4=0.5*AS=0.5*1.55=0.775,
摩擦阻力系数计算如下:
第五个弯和第四个弯以及其中间的连接段直管构成了一个U型弯管,且第五个弯的相对曲率半径为R/D=2,两个弯之间的相对连接段长度lk4/D=20,则第五个弯的局部阻力系数计算如下:
A5=0.3*AS=0.3*1.7=0.51,
摩擦阻力系数计算如下:
第六个弯和第五个弯以及其中间的连接段直管构成了一个U型弯管,且第六个弯的相对曲率半径为R/D=1,两个弯之间的相对连接段长度lk5/D=18,则第六个弯的局部阻力系数计算如下:
A6=0.4*AU=0.4*1.67=0.667,
摩擦阻力系数计算如下:
第七个弯和第六个弯以及其中间的连接段直管构成了一个U型弯管,且第七个弯的相对曲率半径为R/D=2,两个弯之间的相对连接段长度lk6/D=6,则第七个弯的局部阻力系数计算如下:
A7=0.3*AU=0.3*1.45=0.435,
摩擦阻力系数计算如下:
第八个弯和第七个弯以及其中间的连接段直管构成了一个S型弯管,且第八个弯的相对曲率半径为R/D=3,两个弯之间的相对连接段长度lk7/D=4,则第八个弯的局部阻力系数计算如下:
A8=0.5*AS=0.5*1.38=0.69,
摩擦阻力系数计算如下:
第九个弯和第八个弯以及其中间的连接段直管构成了一个U型弯管,且第九个弯的相对曲率半径为R/D=1,两个弯之间的相对连接段长度lk8/D=4,则第九个弯的局部阻力系数计算如下:
A9=0.4*AU=0.4*1.38=0.54,
摩擦阻力系数计算如下:
对以上结果进行叠加:
采用数值计算的方法得到多弯管阻力系数值为2.50,两者之间的相差在3%之内,证明本发明提出的方法精确可靠。
Claims (1)
1.一种任意组合弯管阻力系数的计算方法,基于同一平面的90°单弯管、两个90°单弯管加之间的连接直管构成的S型弯管以及U型弯管的阻力系数计算公式,其特征是,对于具有两个以上同一个平面内的90°弯及其连接直管构成的任意组合弯管,lkj/D≥4,lkj为相邻两个弯之间的连接段长度,D为管径;相对曲率半径1≤R/D≤3,R为弯管轴线的转弯半径;雷诺数2×104≤Re≤105且管壁光滑的圆截面管,从流体入口端到出口端通过参数修正的方法获得各弯头及连接段直管的阻力系数,然后对各个部分的阻力系数进行叠加,从而计算出整个组合弯管的总阻力系数:
为第i弯的单弯管局部阻力系数,采用单弯管局部阻力系数值计算,即:
为第i弯管弯曲部分的摩擦阻力系数;
为连接段的总摩擦阻力系数,lkj为第j个连接段;
Ai为第i弯的局部阻力系数的修正系数;
计算中需要确定的值只有局部阻力系数的修正系数Ai值,流体入口第二个弯以后的局部阻力系数Ai只与前一个弯和连接段长度有关,取值原则如下:
1)当连接段长度lkj/D>25时,该连接段之后那个弯管的局部阻力采用流体入口后第一个弯的计算公式计算;
2)第一个弯的局部阻力系数修正系数A1与相对曲率半径有关:
R/D=1时,A1=0.6AS/U;
R/D=2时,A1=0.7AS/U;
R/D≥2.5时,A1=0.5AS/U。
其中AS/U为由第一个弯和第二个弯组成的S/U型弯管的局部阻力系数修正系数A值,在Idelichik,《Handbook of hydraulic resistance(3rd Edition)》,2008中查得,其值是lkj/D的函数,若两弯的相对曲率半径不同,按照第一个弯的相对曲率半径计算;
3)第二个弯及其后面的弯管局部阻力系数的修正系数Ai值如下:
R/D=1时,Ai=0.4AS/U;
R/D=2时,Ai=0.3AS/U;
R/D≥2.5时,Ai=0.5AS/U;
其中AS/U为第i弯及其前一个弯组成的S/U型弯管的局部阻力系数修正系数A值,在Idelichik,《Handbook of hydraulic resistance(3rd Edition)》,2008中查得,若两弯的相对曲率半径不同,按照第i弯的相对曲率半径计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410539203.8A CN104331603B (zh) | 2014-10-13 | 2014-10-13 | 一种任意组合弯管阻力系数的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410539203.8A CN104331603B (zh) | 2014-10-13 | 2014-10-13 | 一种任意组合弯管阻力系数的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104331603A true CN104331603A (zh) | 2015-02-04 |
CN104331603B CN104331603B (zh) | 2018-01-02 |
Family
ID=52406326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410539203.8A Expired - Fee Related CN104331603B (zh) | 2014-10-13 | 2014-10-13 | 一种任意组合弯管阻力系数的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104331603B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107291984A (zh) * | 2017-05-22 | 2017-10-24 | 广州市市政工程设计研究总院 | 一种综合管廊通风系统当量摩擦阻力系数确定方法 |
CN107862147A (zh) * | 2017-11-16 | 2018-03-30 | 华北电力科学研究院有限责任公司 | 正方截面外部斜切弯头局部阻力的确定方法 |
CN107917089A (zh) * | 2017-11-17 | 2018-04-17 | 上海交通大学 | 循环水泵泵叶开度优化调节方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000002371A (ja) * | 1998-06-15 | 2000-01-07 | Hitachi Ltd | エルボ配管,曲げ管及びその製造方法並びに曲げ管の製造装置 |
CN103353908A (zh) * | 2013-06-20 | 2013-10-16 | 江苏大学 | 一种基于数值计算的管路阻力系数精确计算方法 |
CN103983313A (zh) * | 2014-05-30 | 2014-08-13 | 重庆长厦安基建筑设计有限公司 | 一种管网阻力系数确定方法及管网风量测量方法 |
-
2014
- 2014-10-13 CN CN201410539203.8A patent/CN104331603B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000002371A (ja) * | 1998-06-15 | 2000-01-07 | Hitachi Ltd | エルボ配管,曲げ管及びその製造方法並びに曲げ管の製造装置 |
CN103353908A (zh) * | 2013-06-20 | 2013-10-16 | 江苏大学 | 一种基于数值计算的管路阻力系数精确计算方法 |
CN103983313A (zh) * | 2014-05-30 | 2014-08-13 | 重庆长厦安基建筑设计有限公司 | 一种管网阻力系数确定方法及管网风量测量方法 |
Non-Patent Citations (1)
Title |
---|
杨任: "方形弯管内部流动与阻力系数的数值计算", 《管道技术与设备》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107291984A (zh) * | 2017-05-22 | 2017-10-24 | 广州市市政工程设计研究总院 | 一种综合管廊通风系统当量摩擦阻力系数确定方法 |
CN107291984B (zh) * | 2017-05-22 | 2020-11-17 | 广州市市政工程设计研究总院有限公司 | 一种综合管廊通风系统当量摩擦阻力系数确定方法 |
CN107862147A (zh) * | 2017-11-16 | 2018-03-30 | 华北电力科学研究院有限责任公司 | 正方截面外部斜切弯头局部阻力的确定方法 |
CN107862147B (zh) * | 2017-11-16 | 2021-03-02 | 华北电力科学研究院有限责任公司 | 正方截面外部斜切弯头局部阻力的确定方法 |
CN107917089A (zh) * | 2017-11-17 | 2018-04-17 | 上海交通大学 | 循环水泵泵叶开度优化调节方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104331603B (zh) | 2018-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103353908B (zh) | 一种基于数值计算的管路阻力系数精确计算方法 | |
Ramos et al. | Surge damping analysis in pipe systems: Modelling and experiments | |
Oka et al. | Energy losses at tees with large area ratios | |
CN101881640A (zh) | 涡街质量流量计 | |
CN104331603A (zh) | 一种任意组合弯管阻力系数的计算方法 | |
CN109855691B (zh) | 一种差分式层流流量测量方法及装置 | |
CN104933269A (zh) | 一种油气管道穿越地震断层的设计方法 | |
CN106126847B (zh) | 一种热塑性增强复合管道设计系数的确定方法 | |
CN103074873A (zh) | 一种非理想流场下的多声道超声流量计声道布置方法 | |
Michael et al. | Comparison of plastic limit and collapse loads in pipe bends with shape imperfections under in-plane bending and an internal pressure | |
CN101493861A (zh) | 一种自然循环热水锅炉水动力数值计算方法 | |
CN102630426A (zh) | 文丘里施肥器的结构优化方法 | |
Kabaciński et al. | Numerical and experimental research on new cross-sections of averaging Pitot tubes | |
CN113221354A (zh) | 一种管道弯曲变形拟合算法 | |
Rennels et al. | Pipe flow | |
CN201707087U (zh) | 涡街质量流量计 | |
CN106289415A (zh) | 一种管路流量计算方法、装置及管路系统 | |
CN103837259B (zh) | 裸露油气管道最低管壁温度的测量方法及装置 | |
CN105160062B (zh) | 一种同程管网水力校核方法 | |
CN205037930U (zh) | 一种双差压气液两相流弯管流量计 | |
CN104421558A (zh) | 一种三通管 | |
Yengale et al. | Analysis of Water Distribution Network for Karanja Village-A Case Study | |
Michalos | A review of available minor loss coefficients for combining and dividing flow tees | |
Zheng et al. | Research on flow diagnosis of multipath ultrasonic flowmeter | |
CN204286514U (zh) | 一种差压式流量计 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180102 Termination date: 20191013 |
|
CF01 | Termination of patent right due to non-payment of annual fee |