CN107291984A - 一种综合管廊通风系统当量摩擦阻力系数确定方法 - Google Patents
一种综合管廊通风系统当量摩擦阻力系数确定方法 Download PDFInfo
- Publication number
- CN107291984A CN107291984A CN201710364516.8A CN201710364516A CN107291984A CN 107291984 A CN107291984 A CN 107291984A CN 201710364516 A CN201710364516 A CN 201710364516A CN 107291984 A CN107291984 A CN 107291984A
- Authority
- CN
- China
- Prior art keywords
- frictional resistance
- equivalent
- pipeline
- coefficient
- comprehensive 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000004364 calculation method Methods 0.000 claims abstract description 108
- 238000012937 correction Methods 0.000 claims abstract description 53
- 238000009423 ventilation Methods 0.000 claims description 59
- 238000012360 testing method Methods 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 4
- 238000009434 installation Methods 0.000 abstract 4
- 238000005273 aeration Methods 0.000 abstract 1
- 238000004891 communication Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 231100000681 Certain safety factor Toxicity 0.000 description 1
- 229910001335 Galvanized steel Inorganic materials 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000008397 galvanized steel Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
- 239000002349 well water Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D29/00—Independent underground or underwater structures; Retaining walls
- E02D29/10—Tunnels or galleries specially adapted to house conduits, e.g. oil pipe-lines, sewer pipes ; Making conduits in situ, e.g. of concrete ; Casings, i.e. manhole shafts, access or inspection chambers or coverings of boreholes or narrow wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mining & Mineral Resources (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Architecture (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Environmental & Geological Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Paleontology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种综合管廊通风系统当量摩擦阻力系数确定方法,包括:采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力;采用数值计算方法计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力;对比同等条件下理论公式与数值计算方法的计算结果,得到第一修正系数;采用数值计算方法计算考虑管线和支架敷设情况时综合管廊舱室内的摩擦阻力;对比同等条件下忽略与考虑管线和支架敷设情况时数值计算方法的计算结果,得到第二修正系数;根据第一修正系数和第二修正系数计算综合管廊通风系统的当量摩擦阻力系数。本发明具有科学、有效和方便的优点,可广泛应用于通风技术领域。
Description
技术领域
本发明涉及通风技术领域,尤其是一种综合管廊通风系统当量摩擦阻力系数确定方法。
背景技术
综合管廊又称共同沟,它是实施统一规划、设计、施工和维护,建于城市地下用于敷设市政公用管线的市政公用设施。
综合管廊通风系统通常利用管廊舱室作为风道,管廊舱室主要包括单独设置的蒸汽热力管道舱室、燃气管道舱室、电力电缆管道舱室和综合给水、中水、热水、空调冷水管及电信管的管道舱室等种类不同的舱室,其通风系统的摩擦阻力有所不同。与矿井、地铁、公路隧道等大截面通风管道相比,综合管廊最大的特点是其舱室内部管线和敷设的支架占据了绝大部分面积,其通风摩擦阻力特性也有较大的区别。而目前国内外均未对如何计算该类管道舱室的通风摩擦阻力提出有效的解决办法。在实际的工程设计中,若忽略各管线及其支架对通风摩擦阻力的影响,则可能会出现通风设备压头偏小的情况,而若考虑这部分阻力因素乘以一定的安全系数,则可能因无依据,易造成风机选型压头偏大的现象,从而造成浪费并增加了将来的运行费用。综合管廊通风系统的摩擦阻力大小与当量摩擦阻力系数有关,故当量摩擦系数又是综合管廊通风系统设备计算选型的基础。然而,目前业内也缺乏计算或确定综合管廊通风系统当量摩擦阻力系数的方法,不够科学、有效和方便,亟待进一步完善和提高。
发明内容
为解决上述技术问题,本发明的目的在于:提供一种科学、有效和方便的,综合管廊通风系统当量摩擦阻力系数确定方法。
本发明所采取的技术方案是:
一种综合管廊通风系统当量摩擦阻力系数确定方法,包括以下步骤:
采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
采用数值计算方法计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
对比同等条件下理论公式与数值计算方法的计算结果,得到第一修正系数,所述第一修正系数用于对数值计算方法得出的摩擦阻力进行修正;
采用数值计算方法计算考虑管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
对比同等条件下忽略与考虑管线和支架敷设情况时数值计算方法的计算结果,得到第二修正系数,所述第二修正系数用于对考虑管线和支架敷设情况时数值计算方法得出的摩擦阻力进行修正;
根据第一修正系数和第二修正系数计算综合管廊通风系统的当量摩擦阻力系数。
进一步,所述采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力这一步骤,其具体为:
采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力Δpm1,所述摩擦阻力Δpm1的计算公式为:或其中,λ为摩擦阻力系数,Rs为管廊水力半径,l为管廊长度,ρ为空气密度,v为管廊内空气平均流速,A为过流断面面积,qv为体积流量。
进一步,所述数值计算方法在计算综合管廊舱室内的摩擦阻力时,将执行以下操作:
S1、建立综合管廊管道舱室内的摩擦阻力数值计算模型;
S2、以综合管廊管道为原型,搭建用于模拟原型的缩尺模型试验台并建立缩尺模型的摩擦阻力数值计算模型;
S3、根据综合管廊管道舱室内的摩擦阻力数值计算模型和输入的计算参数计算原型的摩擦阻力;
S4、根据缩尺模型的摩擦阻力数值计算模型和输入的计算参数计算缩尺模型的摩擦阻力;
S5、对比同等条件下原型的摩擦阻力与缩尺模型的摩擦阻力的差值是否在设定的误差范围内,若是,则以原型的摩擦阻力作为数值计算方法的计算结果,反之,则在输入的新计算参数后返回步骤S3。
进一步,所述步骤S2包括:
选择几何相似比,搭建用于模拟原型的缩尺模型试验台,所述缩尺模型试验台根据几何相似比选择与原型相应粗糙度的材料并使得其管线敷设情况与原型相对应;
从缩尺模型试验台获取试验数据,并建立缩尺模型的初步摩擦阻力数值计算模型;
对缩尺模型的初步摩擦阻力数值计算模型进行量纲分析和相似性分析;
根据量纲分析和相似性分析的结果对缩尺模型的初步摩擦阻力数值计算模型进行修正,得出缩尺模型的最终摩擦阻力数值计算模型。
进一步,所述根据第一修正系数和第二修正系数计算综合管廊通风系统的当量摩擦阻力系数这一步骤,其包括:
以数值计算方法计算出的考虑管线和支架敷设情况时综合管廊舱室内的摩擦阻力Δpm3作为综合管廊通风系统的当量摩擦阻力Δpm3,得到当量摩擦阻力Δpm3与第一修正系数k1和第二修正系数k2的修正关系式,所述得到的修正关系式为:
其中,Δpm2为数值计算方法计算出的忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力,λ′为综合管廊通风系统的当量摩擦阻力系数,R′s为综合管廊通风系统的当量水力半径;
根据得到的修正表达式和是否考虑管线和支架敷设情况计算综合管廊通风系统的当量摩擦阻力系数λ′:若忽略管线和支架敷设情况,则R′s=Rs,λ′=k1k2λ;若考虑管线和支架敷设情况,则R′s=(断面面积-管线和支架敷设面积)÷(断面周长+管线和支架周长),
进一步,还包括根据综合管廊通风系统的当量摩擦阻力系数λ′对管线和支架敷设的情况按设定的划分标准进行划分,得出管线和支架敷设情况与当量直径的对应关系的步骤。
进一步,所述设定的划分标准包括两种划分标准:一种是按管线类型、大小尺寸、数量和层数参数进行划分;另一种是按湿周、水力半径以及当量直径的定义进行划分。
进一步,还包括在综合管廊具有非水平直线的走向时,对综合管廊通风系统的当量摩擦阻力系数进行进一步修正的步骤,其中,具有非水平直线的走向是指具有弯度或坡度。
进一步,还包括将综合管廊各类舱室各种情况下的当量摩擦阻力系数及当量摩擦阻力进行汇总,并以图表的形式进行展现的步骤,其中,各种情况包括不同过流断面、不同管线、不同支架敷设方式以及不同弯度和坡度。
本发明的有益效果是:增设了对比同等条件下理论公式与数值计算方法的计算结果以及对比同等条件下忽略与考虑管线和支架敷设情况时数值计算方法的计算结果的步骤,综合采用了理论公式计算、数值计算和二次修正的方法来确定综合管廊通风系统的当量摩擦阻力系数,克服了现有技术因忽略各管线及其支架对通风摩擦阻力的影响而导致通风设备压头偏小或因考虑各管线及其支架的对通风摩擦阻力的影响而造成风机选型压头偏大的缺陷,更加科学、有效和方便。
附图说明
图1为本发明一种综合管廊通风系统当量摩擦阻力系数确定方法的整体流程图;
图2为实施例一广州亚运城综合管廊管道舱断面图;
图3为实施例一水力半径的求解示意图;
图4是实施例一考虑弯度时管廊参数的俯视图;
图5是实施例一考虑坡度时管廊参数的正视图;
图6是实施例一摩擦阻力的拟合曲线图;
图7是实施例一修正系数的拟合曲线图。
具体实施方式
参照图1,一种综合管廊通风系统当量摩擦阻力系数确定方法,包括以下步骤:
采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
采用数值计算方法计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
对比同等条件下理论公式与数值计算方法的计算结果,得到第一修正系数,所述第一修正系数用于对数值计算方法得出的摩擦阻力进行修正;
采用数值计算方法计算考虑管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
对比同等条件下忽略与考虑管线和支架敷设情况时数值计算方法的计算结果,得到第二修正系数,所述第二修正系数用于对考虑管线和支架敷设情况时数值计算方法得出的摩擦阻力进行修正;
根据第一修正系数和第二修正系数计算综合管廊通风系统的当量摩擦阻力系数。
其中,忽略管线和支架敷设情况是指,在计算时将管线和支架情况排除在考虑的范围之外,即在计算时不考虑管线和支架情况。
进一步作为优选的实施方式,所述采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力这一步骤,其具体为:
采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力Δpm1,所述摩擦阻力Δpm1的计算公式为:或其中,λ为摩擦阻力系数,Rs为管廊水力半径,l为管廊长度,ρ为空气密度,v为管廊内空气平均流速,A为过流断面面积,qv为体积流量。
进一步作为优选的实施方式,所述数值计算方法在计算综合管廊舱室内的摩擦阻力时,将执行以下操作:
S1、建立综合管廊管道舱室内的摩擦阻力数值计算模型;
S2、以综合管廊管道为原型,搭建用于模拟原型的缩尺模型试验台并建立缩尺模型的摩擦阻力数值计算模型;
S3、根据综合管廊管道舱室内的摩擦阻力数值计算模型和输入的计算参数计算原型的摩擦阻力;
S4、根据缩尺模型的摩擦阻力数值计算模型和输入的计算参数计算缩尺模型的摩擦阻力;
S5、对比同等条件下原型的摩擦阻力与缩尺模型的摩擦阻力的差值是否在设定的误差范围内,若是,则以原型的摩擦阻力作为数值计算方法的计算结果,反之,则在输入的新计算参数后返回步骤S3。
进一步作为优选的实施方式,所述步骤S2包括:
选择几何相似比,搭建用于模拟原型的缩尺模型试验台,所述缩尺模型试验台根据几何相似比选择与原型相应粗糙度的材料并使得其管线敷设情况与原型相对应;
从缩尺模型试验台获取试验数据,并建立缩尺模型的初步摩擦阻力数值计算模型;
对缩尺模型的初步摩擦阻力数值计算模型进行量纲分析和相似性分析;
根据量纲分析和相似性分析的结果对缩尺模型的初步摩擦阻力数值计算模型进行修正,得出缩尺模型的最终摩擦阻力数值计算模型。
进一步作为优选的实施方式,所述根据第一修正系数和第二修正系数计算综合管廊通风系统的当量摩擦阻力系数这一步骤,其包括:
以数值计算方法计算出的考虑管线和支架敷设情况时综合管廊舱室内的摩擦阻力Δpm3作为综合管廊通风系统的当量摩擦阻力Δpm3,得到当量摩擦阻力Δpm3与第一修正系数k1和第二修正系数k2的修正关系式,所述得到的修正关系式为:
其中,Δpm2为数值计算方法计算出的忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力,λ′为综合管廊通风系统的当量摩擦阻力系数,R′s为综合管廊通风系统的当量水力半径;
根据得到的修正表达式和是否考虑管线和支架敷设情况计算综合管廊通风系统的当量摩擦阻力系数λ′:若忽略管线和支架敷设情况,则R′s=Rs,λ′=k1k2λ;若考虑管线和支架敷设情况,则R′s=(断面面积-管线和支架敷设面积)÷(断面周长+管线和支架周长),
进一步作为优选的实施方式,还包括根据综合管廊通风系统的当量摩擦阻力系数λ′对管线和支架敷设的情况按设定的划分标准进行划分,得出管线和支架敷设情况与当量直径的对应关系的步骤。
进一步作为优选的实施方式,所述设定的划分标准包括两种划分标准:一种是按管线类型、大小尺寸、数量和层数参数进行划分;另一种是按湿周、水力半径以及当量直径的定义进行划分。
进一步作为优选的实施方式,还包括在综合管廊具有非水平直线的走向时,对综合管廊通风系统的当量摩擦阻力系数进行进一步修正的步骤,其中,具有非水平直线的走向是指具有弯度或坡度。
进一步作为优选的实施方式,还包括将综合管廊各类舱室各种情况下的当量摩擦阻力系数及当量摩擦阻力进行汇总,并以图表的形式进行展现的步骤,其中,各种情况包括不同过流断面、不同管线、不同支架敷设方式以及不同弯度和坡度。
下面结合说明书附图和具体实施例对本发明作进一步解释和说明。
实施例一
针对现有技术缺乏计算或确定综合管廊通风系统当量摩擦阻力系数的方法的问题,本发明提出了当量摩擦阻力系数的概念,并提供了一种采用现场实地测量、缩尺模型试验研究、数值计算模拟仿真和理论经验公式修正等多种手段对综合管廊通风系统当量摩擦阻力系数进行确定的方法。该方法可以作为通风系统设备计算选型的基础,解决了不同断面、不同管线、不同支架敷设方式以及不同弧度和坡度的通风系统当量摩擦阻力系数的选取缺乏依据的难题。
以广州亚运城综合管廊管道舱为例,其管道舱断面高3.1m,宽2.8m,内部布置管线包括通信电缆1、中水管3、消防水管2,如图2所示,其当量摩擦阻力系数的确定过程具体包括:
(一)采用现有理论经验公式计算忽略管线和支架敷设情况下综合管廊舱室内的摩擦阻力Δpm1。
Δpm1的计算公式为:或
式中,λ为摩擦阻力系数,单位为m;Rs为管廊水力半径,单位为m;l为管廊长度,单位为m;ρ为空气密度,单位为kg/m3;v为管廊内空气平均流速,单位为m/s;A为过流断面面积,单位为m2;qv为体积流量。
Δpm1的计算参数又可具体分为以下2种情况:
(1)层流:
(2)紊流
紊流时,不同区域的计算参数为:
1)紊流水力光滑区:
2)紊流过渡区:
3)紊流粗糙区:
以上各式中,Re为雷诺系数,D为管廊直径,单位为m;K为管廊的粗糙度,单位为mm,薄钢板或镀锌钢板的K=0.15~0.18;砌墙体的K=3~6;混凝土的K=1~3。
而水力半径P为湿周,单位为m。
(二)在与理论公式流速或流量相同的条件下,用数值计算方法计算忽略管线和支架敷设情况下综合管廊舱室内的摩擦阻力Δpm2。
(三)对比理论经验公式与数值计算方法的计算结果,然后对数值计算方法得出的摩擦阻力进行修正,其修正系数用k1表示,则有:
(四)在与忽略管线和支架敷设情况时的流速或流量相同的条件下,用数值计算方法分别计算考虑管线和支架敷设情况下综合管廊舱室内的摩擦阻力Δpm3。
(五)对比忽略与考虑管线和支架敷设情况下数值计算方法的计算结果,对考虑管线和支架敷设情况下数值计算方法得出的摩擦阻力进行修正,其修正系数用k2表示,则有:
则当量摩擦阻力即从而分别得到相应的速度当量直径Dv和流量当量直径DL,以及得到当量摩擦阻力系数λ′与当量水力半径Rs′的修正关系式。
本发明数值计算方法的可靠性可通过以下步骤进行保证:
1)、首先对广州亚运城综合管廊管道舱的通风摩擦阻力进行现场实地测量,同时建立相关数值计算模型并对其进行仿真计算,使计算结果与实测结果的误差控制在工程允许范围内。
2)、选择适当的几何相似比,搭建模拟原型的缩尺模型试验台,根据相似比选择与原型相应粗糙度的材料,保持模型内部管线敷设情况与原型相对应,并获取试验数据;同时建立缩尺模型相关的数值计算模型并对其进行仿真计算。
在进行相似性研究时,由于原型与缩尺模型两者惯性力与浮升力所占的比例不同,如果只保证原型与缩尺模型雷诺数相等,则势必会存在一定误差,因此要先进行量纲分析和相似性分析,再在缩尺模型所获取的数据反推到实际原型时做一定的修正。
对于不可压缩流体的定常流动来说,如果其缩尺模型流动和原型流动力学相似,则它们的弗劳德数Fr、欧拉数Eu和雷诺数Re必须各自相等,此时,Fr和Re分别要各自满足和kv和kl分别为速度比例尺和长度比例尺,当kl≠1时不成立,而欧拉数准则与上述两个准则没有矛盾,因此需要先进行取舍,再对模型结果进行修正,具体做法为:在保证雷诺数相等的条件下,对当量摩擦阻力系数进行拟合修正,得出λ关于雷诺数Re和壁面相对粗糙度的关系式或者在弗劳德数和雷诺数都不相等的情况下,对当量摩擦阻力系数进行拟合修正,得出λ关于弗雷德数Fr、雷诺数Re和壁面相对粗糙度的关系式得出最佳相似方案。
3)、将仿真计算和理论公式计算的数据进行对比,使误差控制在工程允许范围内,并且分别与原型进行对比,得出原型与缩尺模型之间的相似修正系数。
4)、对每一种工况下的数值计算结果,都通过相应工况的缩尺模型所获取的数据来验证,以保证结果的可靠性。
(六)对管线和支架敷设的情况按一定的划分标准进行划分,得出管线和支架敷设情况与当量直径的一一对应关系。
本发明的划分标准分为以下两种:
第一种是按管线类型、大小尺寸、数量和层数等参数进行划分。
其中,通信线缆按层数进行划分,用n1表示通信线缆敷设的层数。
电力电缆按层数进行划分,用n2表示表示电力电缆敷设的层数。
而给水、再生水管道,排水管渠,天然气管道,热力管道和消防水管道则按大小尺寸划分,具体划分情况如下表1所示:
表1管线划分表
公称直径DN | 数量 |
200以下 | n3 |
200-400 | n4 |
400-600 | n5 |
600-800 | n6 |
800-1000 | n7 |
另一种划分标准是按湿周、水力半径以及当量直径的定义进行划分:将管线所占的截面面积从过流断面面积中舍去,并将管线周长纳入湿周,从而对水力半径进行修正。
如图3所示,宽、高分别为a、b的管廊内部布置了一条半径为r的管线,忽略管线敷设情况时,其水力半径而考虑管线敷设情况时,其水力半径代入理论经验公式进行计算,并与数值计算方法的结果进行对比,可得出其修正系数。
(七)在管廊是具有一定弯度或坡度等非水平直线的走向时,对其通风摩擦阻力系数进行进一步修正。
1)考虑弯度的管廊摩擦阻力系数确定方法:
按照一定标准(与等长度的直线管廊对比),在不同弧长和半径的情况下得出相应的修正系数kw,如图4所示,则有:
kw=kw(L,R)
其中,L为弧长,单位为m;R为半径,单位为m。
2)考虑坡度的管廊摩擦阻力系数确定方法:
按照一定标准(与等长度的直线管廊对比),在不同角度和高度的情况下得出相应的修正系数kp,如图5所示,则有:
kp=kp(φ,H)
其中,φ为角度,单位为°;H为高度,单位为m。
(八)最后通过大量的数值计算来建立数据库,并将综合管廊各类舱室不同断面、不同管线、不同支架敷设方式以及不同水平弯度和纵向坡度的摩擦阻力系数及其摩擦阻力进行汇总以图表的形式展现,如表1、表2、表3、图6、图7所示。
表2当量直径查询表
序号 | n1 | n2 | n3 | n4 | n5 | n6 | n7 | Dv | DL |
1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | Dv1 | DL1 |
2 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | Dv2 | DL2 |
3 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | Dv3 | DL3 |
4 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | Dv4 | DL4 |
5 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | Dv5 | DL5 |
… | … | … | … | … | … | … | … | … | … |
由表2可得到:
Dv=Dv(n1,n2,n3,n4,n5,n6,n7)
DL=DL(n1,n2,n3,n4,n5,n6,n7)
表3是摩擦阻力系数查询表
Dv | v | △P | λ | DL | qv | △P | λ |
Dv1 | v1 | DL1 | qv1 | ||||
Dv1 | v2 | DL1 | qv2 | ||||
Dv2 | v1 | DL2 | qv1 | ||||
Dv2 | v2 | DL2 | qv2 | ||||
… | … | … | … |
以上是对本发明的较佳实施进行了具体说明,但本发明并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。
Claims (9)
1.一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:包括以下步骤:
采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
采用数值计算方法计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
对比同等条件下理论公式与数值计算方法的计算结果,得到第一修正系数,所述第一修正系数用于对数值计算方法得出的摩擦阻力进行修正;
采用数值计算方法计算考虑管线和支架敷设情况时综合管廊舱室内的摩擦阻力;
对比同等条件下忽略与考虑管线和支架敷设情况时数值计算方法的计算结果,得到第二修正系数,所述第二修正系数用于对考虑管线和支架敷设情况时数值计算方法得出的摩擦阻力进行修正;
根据第一修正系数和第二修正系数计算综合管廊通风系统的当量摩擦阻力系数。
2.根据权利要求1所述的一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:所述采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力这一步骤,其具体为:
采用理论公式计算忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力Δpm1,所述摩擦阻力Δpm1的计算公式为:或其中,λ为摩擦阻力系数,Rs为管廊水力半径,l为管廊长度,ρ为空气密度,v为管廊内空气平均流速,A为过流断面面积,qv为体积流量。
3.根据权利要求1所述的一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:所述数值计算方法在计算综合管廊舱室内的摩擦阻力时,将执行以下操作:
S1、建立综合管廊管道舱室内的摩擦阻力数值计算模型;
S2、以综合管廊管道为原型,搭建用于模拟原型的缩尺模型试验台并建立缩尺模型的摩擦阻力数值计算模型;
S3、根据综合管廊管道舱室内的摩擦阻力数值计算模型和输入的计算参数计算原型的摩擦阻力;
S4、根据缩尺模型的摩擦阻力数值计算模型和输入的计算参数计算缩尺模型的摩擦阻力;
S5、对比同等条件下原型的摩擦阻力与缩尺模型的摩擦阻力的差值是否在设定的误差范围内,若是,则以原型的摩擦阻力作为数值计算方法的计算结果,反之,则在输入的新计算参数后返回步骤S3。
4.根据权利要求3所述的一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:所述步骤S2包括:
选择几何相似比,搭建用于模拟原型的缩尺模型试验台,所述缩尺模型试验台根据几何相似比选择与原型相应粗糙度的材料并使得其管线敷设情况与原型相对应;
从缩尺模型试验台获取试验数据,并建立缩尺模型的初步摩擦阻力数值计算模型;
对缩尺模型的初步摩擦阻力数值计算模型进行量纲分析和相似性分析;
根据量纲分析和相似性分析的结果对缩尺模型的初步摩擦阻力数值计算模型进行修正,得出缩尺模型的最终摩擦阻力数值计算模型。
5.根据权利要求2所述的一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:所述根据第一修正系数和第二修正系数计算综合管廊通风系统的当量摩擦阻力系数这一步骤,其包括:
以数值计算方法计算出的考虑管线和支架敷设情况时综合管廊舱室内的摩擦阻力Δpm3作为综合管廊通风系统的当量摩擦阻力Δpm3,得到当量摩擦阻力Δpm3与第一修正系数k1和第二修正系数k2的修正关系式,所述得到的修正关系式为:
<mrow>
<msub>
<mi>&Delta;p</mi>
<mrow>
<mi>m</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<mi>&lambda;</mi>
<mfrac>
<mi>l</mi>
<mrow>
<mn>4</mn>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msup>
<mi>&rho;v</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>=</mo>
<msup>
<mi>&lambda;</mi>
<mo>&prime;</mo>
</msup>
<mfrac>
<mi>l</mi>
<mrow>
<mn>4</mn>
<msubsup>
<mi>R</mi>
<mi>s</mi>
<mo>&prime;</mo>
</msubsup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msup>
<mi>&rho;v</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>;</mo>
</mrow>
其中,Δpm2为数值计算方法计算出的忽略管线和支架敷设情况时综合管廊舱室内的摩擦阻力,λ′为综合管廊通风系统的当量摩擦阻力系数,R′s为综合管廊通风系统的当量水力半径;
根据得到的修正表达式和是否考虑管线和支架敷设情况计算综合管廊通风系统的当量摩擦阻力系数λ′:若忽略管线和支架敷设情况,则R′s=Rs,λ′=k1k2λ;若考虑管线和支架敷设情况,则R′s=(断面面积-管线和支架敷设面积)÷(断面周长+管线和支架周长),
6.根据权利要求5所述的一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:还包括根据综合管廊通风系统的当量摩擦阻力系数λ′对管线和支架敷设的情况按设定的划分标准进行划分,得出管线和支架敷设情况与当量直径的对应关系的步骤。
7.根据权利要求6所述的一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:所述设定的划分标准包括两种划分标准:一种是按管线类型、大小尺寸、数量和层数参数进行划分;另一种是按湿周、水力半径以及当量直径的定义进行划分。
8.根据权利要求6所述的一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:还包括在综合管廊具有非水平直线的走向时,对综合管廊通风系统的当量摩擦阻力系数进行进一步修正的步骤,其中,具有非水平直线的走向是指具有弯度或坡度。
9.根据权利要求8所述的一种综合管廊通风系统当量摩擦阻力系数确定方法,其特征在于:还包括将综合管廊各类舱室各种情况下的当量摩擦阻力系数及当量摩擦阻力进行汇总,并以图表的形式进行展现的步骤,其中,各种情况包括不同过流断面、不同管线、不同支架敷设方式以及不同弯度和坡度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710364516.8A CN107291984B (zh) | 2017-05-22 | 2017-05-22 | 一种综合管廊通风系统当量摩擦阻力系数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710364516.8A CN107291984B (zh) | 2017-05-22 | 2017-05-22 | 一种综合管廊通风系统当量摩擦阻力系数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107291984A true CN107291984A (zh) | 2017-10-24 |
CN107291984B CN107291984B (zh) | 2020-11-17 |
Family
ID=60094502
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710364516.8A Active CN107291984B (zh) | 2017-05-22 | 2017-05-22 | 一种综合管廊通风系统当量摩擦阻力系数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107291984B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110108633A (zh) * | 2019-04-26 | 2019-08-09 | 辽宁工程技术大学 | 一种矿井通风摩擦阻力系数测定的方法 |
CN110569515A (zh) * | 2018-12-24 | 2019-12-13 | 中国电力科学研究院有限公司 | 用于计算综合管廊风机压损的分析方法 |
CN111046468A (zh) * | 2019-12-12 | 2020-04-21 | 悉地(苏州)勘察设计顾问有限公司 | 地下综合管廊抗震设计方法 |
CN111220501A (zh) * | 2019-12-04 | 2020-06-02 | 西南石油大学 | 一种输油管道中减阻剂减阻效果的在线评价方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100921478B1 (ko) * | 2008-05-30 | 2009-10-13 | 한국건설기술연구원 | 수치방정식의 보정계수를 적용한 터널 내의 소요환기량산출방법 |
CN101864978A (zh) * | 2010-05-21 | 2010-10-20 | 王明年 | 曲线隧道通风沿程阻力的确定方法 |
CN102636328A (zh) * | 2012-04-11 | 2012-08-15 | 中铁西南科学研究院有限公司 | 一种曲线隧道通风沿程摩擦阻力系数模拟试验方法 |
CN103353908A (zh) * | 2013-06-20 | 2013-10-16 | 江苏大学 | 一种基于数值计算的管路阻力系数精确计算方法 |
CN104331603A (zh) * | 2014-10-13 | 2015-02-04 | 南京航空航天大学 | 一种任意组合弯管阻力系数的计算方法 |
CN105740581A (zh) * | 2016-03-15 | 2016-07-06 | 上海交通大学 | 环形流道转动下的摩擦因子修正方法 |
CN106682369A (zh) * | 2017-02-27 | 2017-05-17 | 常州英集动力科技有限公司 | 供热管网水力仿真模型辨识修正方法及系统、操作方法 |
-
2017
- 2017-05-22 CN CN201710364516.8A patent/CN107291984B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100921478B1 (ko) * | 2008-05-30 | 2009-10-13 | 한국건설기술연구원 | 수치방정식의 보정계수를 적용한 터널 내의 소요환기량산출방법 |
CN101864978A (zh) * | 2010-05-21 | 2010-10-20 | 王明年 | 曲线隧道通风沿程阻力的确定方法 |
CN102636328A (zh) * | 2012-04-11 | 2012-08-15 | 中铁西南科学研究院有限公司 | 一种曲线隧道通风沿程摩擦阻力系数模拟试验方法 |
CN103353908A (zh) * | 2013-06-20 | 2013-10-16 | 江苏大学 | 一种基于数值计算的管路阻力系数精确计算方法 |
CN104331603A (zh) * | 2014-10-13 | 2015-02-04 | 南京航空航天大学 | 一种任意组合弯管阻力系数的计算方法 |
CN105740581A (zh) * | 2016-03-15 | 2016-07-06 | 上海交通大学 | 环形流道转动下的摩擦因子修正方法 |
CN106682369A (zh) * | 2017-02-27 | 2017-05-17 | 常州英集动力科技有限公司 | 供热管网水力仿真模型辨识修正方法及系统、操作方法 |
Non-Patent Citations (4)
Title |
---|
刘承东,唐宏辉,谢艺强,胡康: "城市综合管廊通风系统风亭设计及优化", 《中国市政工程》 * |
张佳佳,沈阳: "蒸汽管道评价指标及计算方法", 《发电与空调》 * |
成锋娜,常海萍等: "90。多弯管阻力系数计算方法", 《航空发动机》 * |
罗继杰,张蔚东: "《通风管道沿程阻力计算选用表》编制及设计应用", 《暖通空调》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110569515A (zh) * | 2018-12-24 | 2019-12-13 | 中国电力科学研究院有限公司 | 用于计算综合管廊风机压损的分析方法 |
CN110108633A (zh) * | 2019-04-26 | 2019-08-09 | 辽宁工程技术大学 | 一种矿井通风摩擦阻力系数测定的方法 |
CN111220501A (zh) * | 2019-12-04 | 2020-06-02 | 西南石油大学 | 一种输油管道中减阻剂减阻效果的在线评价方法 |
CN111220501B (zh) * | 2019-12-04 | 2022-04-29 | 西南石油大学 | 一种输油管道中减阻剂减阻效果的在线评价方法 |
CN111046468A (zh) * | 2019-12-12 | 2020-04-21 | 悉地(苏州)勘察设计顾问有限公司 | 地下综合管廊抗震设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107291984B (zh) | 2020-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107291984B (zh) | 一种综合管廊通风系统当量摩擦阻力系数确定方法 | |
Gao et al. | Study of the shape optimization of a tee guide vane in a ventilation and air-conditioning duct | |
Liu et al. | Numerical modeling and parametric study of a vertical earth-to-air heat exchanger system | |
US10272494B1 (en) | Fabricating conduits | |
Grishchenko et al. | The TALL-3D facility design and commissioning tests for validation of coupled STH and CFD codes | |
Li et al. | Experimental reduced-scale study on the resistance characteristics of the ventilation system of a utility tunnel under different pipeline layouts | |
CN109948182B (zh) | 一种针对中深层地热井井间距的计算方法 | |
CN101864978B (zh) | 曲线隧道通风沿程阻力的确定方法 | |
Annan et al. | Effect of minor losses during steady flow in transmission pipelines–Case study “water transmission system upgrade in Northern Saudi Arabia” | |
Josiah et al. | An experimental investigation of head loss through trash racks in conveyance systems | |
CN111539130A (zh) | 中深层地埋管地热换热器钻孔深度的设计计算方法 | |
CN109490354A (zh) | 一种基于相似理论的电力舱缩尺试验台搭建方法 | |
CN117235878B (zh) | 一种基于bim技术的综合管廊管线避免碰撞的排布方法 | |
CN108644993B (zh) | 一种均匀送风的低流动阻力管道装置 | |
Neale et al. | Flow characteristics of PVC sewer pipe | |
CN108386983B (zh) | 一种用于通风空调风管的低阻力分流三通构件 | |
CN2220574Y (zh) | 双管预制保温管 | |
CN107862147B (zh) | 正方截面外部斜切弯头局部阻力的确定方法 | |
CN115930116A (zh) | 一种基于时间与风险维度的燃气管线泄漏定位方法 | |
CN113605953B (zh) | 基于高压空气供风的隧道施工通风系统及其构建方法 | |
CN109635380A (zh) | 一种管廊电力舱防火间隔与电缆载流量定量关系确定方法 | |
CN104850708A (zh) | 一种玻璃钢管线管输系统温度场非稳态特性分析方法 | |
CN105091967B (zh) | 一种利用管道本身或管道内障碍物的流量测量方法 | |
Meky et al. | Effect of Water Depth Change on Manning Coefficient for Partially-Filled Circular Culverts. | |
CN109482689A (zh) | 建筑钢管滚压煨弯设备和建筑钢管煨弯及安装施工方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 510060 East seat, No. 348 East Ring Road, Guangzhou, Guangdong, Yuexiu District Applicant after: Guangzhou municipal engineering design and Research Institute Co., Ltd. Address before: 510060 No. 348 East Ring Road, Guangzhou, Guangdong Applicant before: Guangzhou City Engineering Design studies total institute |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |