CN111680460B - 水力自动翻板闸门流量系数的确定方法 - Google Patents
水力自动翻板闸门流量系数的确定方法 Download PDFInfo
- Publication number
- CN111680460B CN111680460B CN202010535337.8A CN202010535337A CN111680460B CN 111680460 B CN111680460 B CN 111680460B CN 202010535337 A CN202010535337 A CN 202010535337A CN 111680460 B CN111680460 B CN 111680460B
- Authority
- CN
- China
- Prior art keywords
- gate
- flow
- formula
- flow coefficient
- weir
- 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.)
- Active
Links
Images
Classifications
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B7/00—Barrages or weirs; Layout, construction, methods of, or devices for, making same
- E02B7/20—Movable barrages; Lock or dry-dock gates
- E02B7/40—Swinging or turning gates
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Structural Engineering (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- Civil Engineering (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Mechanical Engineering (AREA)
- Architecture (AREA)
- Computational Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
水力自动翻板闸门流量系数的确定方法,属于水利设施流量系数测试和研究领域,尤其涉及一种应用于测试研究水力自动翻板闸门综合流量系数的方法。本发明方法包括流量系数计算公式的推导,并计算流量系数值,以及建立试验模型,并计算实测流量系数值,将计算值和实测值进行比较,验证公式计算精度。本发明方法解决不同倾角下的水力翻板闸门泄流能力的计算难题。并可用于其他相同结构的水力自动翻板闸门的过流流量计算,为水力自动翻板闸门的设计和运行提供科学依据。
Description
技术领域
本发明属于水利设施流量系数测试和研究领域,尤其涉及一种应用于测试研究水力自动翻板闸门综合流量系数的方法。
背景技术
水力自动翻板闸门凭借水压力、水重、自重来实施启闭。闸门开启后,门上门下分两层泄流。门顶上面的水流有自由表面,属于堰流,且水流流态随着闸门开启角度的增大,逐渐由薄壁堰流过渡到宽顶堰流。门板下面的水流呈孔流状态流向下游。水力自动翻板闸门闸顶过流后,门顶及下部孔口同时泄水,这两部分水流又相互混合,水流处于一种闸孔泄流与堰流相叠加的复杂水力现象,这样就无法直接套用堰流公式或闸孔出流公式进行闸门流量的计算。目前此种情况下的闸门过流流量计算,在理论上没有精度高的计算方法,因此,需要通过模型测试并验证推导出闸门综合流量系数,提出一种测试水力自动翻板闸门流量系数的方法。
发明内容
本发明目的是解决目前没有一种精度高的测试水力翻板闸门流量系数计算方法的问题,提出综合考虑闸孔泄流与堰流相叠加的复杂水力现象,解决不同倾角下的水力翻板闸门泄流能力的计算方法。
本发明水力自动翻板闸门流量系数的确定方法,其特征在于该方法包括流量系数计算公式的推导,并计算流量系数值,以及建立试验模型,并计算实测流量系数值,将计算值和实测值进行比较,验证公式计算精度;
所述的流量系数计算值采用的计算公式为:
式(1)中,d为闸门在对应倾角下的阻挡系数;
P1为上游堰高,单位m;
H为闸前堰顶水头,单位m;
L为闸门直立关闭挡水时的闸门门高,单位m;
所述的式(1)中d值通过以下公式计算:
式(2)中,Q1为水力翻板闸门在某一闸门倾角下的实测流量,单位m3/s;
Q2为未考虑闸门阻挡作用时的宽顶堰流量,即理论流量,单位m3/s;
所述的式(2)中Q2通过以下公式计算:
式(3)中,n为闸孔数目;
b为闸孔净宽,单位m;
H0=H+v0 2/2g为包括行进水头在内的上游总水头,单位m;v0为闸前行进流速;g 为重力加速度;
所述的(2)中Q1通过以下公式计算:
所述的试验模型由三角形量水堰、闸前水库库区、水力翻板闸门通过粘接或栓接组装而成;三角形量水堰设于库区上游,水堰岸边翻板闸门前布置水位测针钢尺,水力翻板闸门面板、支腿由有机玻璃制成。
所述的流量系数实测值采用的计算公式为:
式(5)中,b=0.40m;
n=1。
本发明与现有技术相比,其有益效果为:
(1)水力翻板闸门闸顶过流后,门顶及下部孔口同时泄水,水流处于一种闸孔泄流与堰流相叠加的复杂水力现象。目前此种情况下的闸门过流流量系数,在理论上没有现成成熟精度高的计算公式。通过实测流量系数验证推导的不同特定闸门倾角下水力自动翻板闸门流量系数提出可靠计算公式,可解决不同倾角下的水力翻板闸门泄流能力的计算难题。并可用于其他相同结构的水力自动翻板闸门的过流流量计算,为水力自动翻板闸门的设计和运行提供科学依据。
(2)模型试验实测值对推导的流量系数计算公式进行了验证,结果表明计算值与实测值的相对误差一般都在5%以内,精度较高,在允许误差范围内,能够满足实际设计工程的需要。
附图说明
图1为自动翻板闸门开启示意图;
图2为自动翻板闸门过流示意图;
图3水位测针钢尺。
图中:1-库区;2-测针钢尺;3-水力翻板闸门;
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合附图本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
实施例1:本发明首先推导水利工程上使用的水利自动翻板闸门的流量系数的一种计算公式,并根据室内闸门水工模型进行了试验以获取流量系数的实测值,对推导的流量系数计算公式进行了精度验证,最后提出一种可靠性好的水利翻板闸门流量系数计算方法。
本发明的技术方案如下:
步骤1,推导流量系数计算公式,具体过程如下:
首先将翻板闸门过流假定为堰流,用堰流公式计算流量时,将翻板闸门看作堰上障碍物,需考虑闸门对过闸水流的阻挡作用,且随着闸门倾角和堰前水头的变化,闸门对水流的阻挡程度不同,故用堰流公式计算闸门流量时,需在流量公式中加入闸门对水流的阻挡系数d。
闸门实测流量用Q1表示。考虑到翻板闸门建在宽顶堰上,闸前的理论流量用宽顶堰流量公式计算,其流量用Q2表示。堰流时的闸前水头为H,闸门实际过流时的闸前水头可看作是闸前水头H减去闸门在过流断面上的投影高度,并假定这两个流量之比与其对应的闸前水头呈正比例关系,先设定它们之间的函数关系式的形式为:
式中,Q1为水力翻板闸门在某一闸门倾角下的实测流量,单位m3/s;Q2为未考虑闸门阻挡作用时的宽顶堰流量,即理论流量,单位m3/s;H为闸前堰顶水头,单位m;L为闸门直立关闭挡水时的闸门门高,单位m;φ为闸门倾角;d为闸门在对应倾角下的阻挡系数。
理论流量Q2由宽顶堰流量堰流公式计算得出,即:
式中,P1为上游堰高,单位m;n为闸孔数目;b为闸孔净宽,单位m;H0=H+v0 2/2g 为包括行进水头在内的上游总水头,单位m;v0为闸前行进流速;g为重力加速度,一般取9.8。
由公式(1-1)、(1-2)可以得出:
则推导出的翻板闸门综合流量系数mf的计算公式可表示为:
以上推导出的闸门流量系数的计算公式,是在宽顶堰流量系数的基础上乘以一个阻挡系数和闸门相对过流面积,与闸门倾角和堰前水头有关,这个阻挡系数大小直接反映了翻板闸门的对水流的阻挡程度。如果直接用该公式,则需计算出阻挡系数d。
由公式(1-1)可推导出特定闸门倾角下对应于不同水位的阻挡系数d值:
通过试验表明:在同一闸门倾角下阻挡系数d变化不大,故可用一个平均综合阻挡系数d来表示该闸门倾角下闸门的阻挡能力。d值由模型试验得出,对于不同的闸门倾角,阻挡系数d值不同。经试验得出的翻板闸门在10°、35°、60°、80°这4个不同闸门倾角下的平均综合阻挡系数d值分别为0.81、0.99、1.10、1.03。
步骤(2),模型试验,对流量系数实测值与计算值比较,验证公式计算精度。
试验模型主要由三角形量水堰、闸前水库库区、翻板闸门组成。三角形量水堰设于库区上游,由灰塑板做成,流量大小是用布置于量水堰岸边处的测针钢尺读出水位值,再利用现有的三角堰堰流公式自动计算出上游来流量值。库区的闸门上游,水位测针钢尺布置在闸门前1.2m(原型为30.0m)处的玻璃水槽墙边,此处水流流态平稳,用于实时测读闸门上游水库水位,测量精度为0.1mm。翻板闸门模型是按重力相似(佛汝德数相似)准则设计的正态模型,与原型水利翻板闸门的比例为1:25,在40cm宽的玻璃水槽内测试单扇水利翻板闸门在不同倾角(开度)下的过流能力。闸门门高28cm,门宽 40cm。闸门面板、支腿用有机玻璃制成,模型制作误差±0.1mm,模型装配形式采用粘接和栓接。
试验研究可选取不同的闸门倾角,分别作为某一特定工况进行。在模型实际测试时,当闸门直立挡水时,闸门倾角为0°,闸门高度L为28cm,闸门全开时,闸门倾角为 80°(即翻板闸门最大倾角),闸门下部开度e为0.178m。
闸门实测流量系数m由试验中测得的水位—流量值,根据公式计算出,实测闸门流量系数计算公式为:
上式中,试验中用的b=0.40m,n=1,H0中的v0利用闸前旋浆流速仪测出。
针对10°、35°、60°、80°这4种不同闸门倾角工况,由公式(1-4)分别计算出闸门流量系数计算值mf,与公式(1-6)计算出闸门流量系数实测值m进行比较,并计算其相对误差。经对试验数据进行分析,翻板闸门流量系数计算值与实测值的最大相对误差为3.74%。
当闸门全开(倾角为80°)时的工况,对比闸门流量系数计算值与实测值,并计算其相对误差,计算结果见表1。
表1闸门全开时流量系数计算值与实测值比较
水位(米) | 366.0 | 366.2 | 366.3 | 366.5 | 366.7 | 366.8 | 367.0 | 367.2 | 367.3 | 367.5 | 367.7 | 367.8 | 368.0 |
计算流量系数m<sub>f</sub> | 0.310 | 0.311 | 0.312 | 0.313 | 0.315 | 0.315 | 0.317 | 0.318 | 0.318 | 0.032 | 0.321 | 0.321 | 0.322 |
实测流量系数m | 0.315 | 0.315 | 0.316 | 0.316 | 0.317 | 0.317 | 0.316 | 0.316 | 0.316 | 0.315 | 0.315 | 0.317 | 0.315 |
相对误差(%) | 1.47 | 1.15 | 1.09 | 0.94 | 0.76 | 0.65 | -0.17 | -0.67 | -0.92 | -1.43 | -1.94 | -1.33 | -2.20 |
三角形量水堰设于库区上游,流量大小是用布置于量水堰岸边处的测针钢尺读出水位值,再利用现成的上述三角堰堰流公式(1-1)和(1-2)自动计算出上游来流量值;库区的闸门上游水位用测针钢尺读取,闸前流速用闸前旋浆流速仪自动测出。
进行模型试验时,在某一特定闸门倾角下,测量上游不同来流量情况下,读取量水堰和上游库区水位测针的读数以及测出闸前流速,建立在该特定闸门倾角下的闸前水位和过水流量的闸门泄流曲线,利用上述公式(1-3)计算在一特定闸门倾角下水力翻板闸门实测的综合流量系数m,并与推导的公式(1-4)计算出的在这一特定倾角下的流量系数mf进行比较,进而验证推导出的翻板闸门综合流量系数计算公式的精确度和可靠性。
根据模型试验对推导的流量系数计算公式的验证结果,表明计算值与实测值的相对误差一般都在5%以内,精度较高,在允许误差范围内,能够满足实际设计工程的需要。本试验最终提出可靠的流量系数计算公式,以解决不同倾角下的水力翻板闸门泄流能力的计算,为水力自动翻板闸门的设计和运行提供科学依据。
尽管这里参照本发明的多个解释性试验实例对本发明进行了描述,但是,应该理解,本领域技术人员可以设计出很多其他的修改和实施方式,这些修改和实施方式将落在本申请公开的原则范围和精神之内。更具体地说,在本申请公开、附图和权利要求的范围内,可以对主题组合布局的组成部件和/或布局进行多种变型和改进。除了对组成部件和/或布局进行的变形和改进外,对于本领域技术人员来说,其他的用途也将是明显的。
Claims (2)
1.水力自动翻板闸门流量系数的确定方法,其特征在于该方法包括流量系数计算公式的推导,并计算流量系数值,以及建立试验模型,并计算实测流量系数值,将计算值和实测值进行比较,验证公式计算精度;
所述的流量系数计算值采用的计算公式为:
式(1)中,d为闸门在对应倾角下的阻挡系数;
P1为上游堰高,单位m;
H为闸前堰顶水头,单位m;
L为闸门直立关闭挡水时的闸门门高,单位m;
所述的式(1)中d值通过以下公式计算:
式(2)中,Q1为水力翻板闸门在某一闸门倾角下的实测流量,单位m3/s;
Q2为未考虑闸门阻挡作用时的宽顶堰流量,即理论流量,单位m3/s;
所述的式(2)中Q2通过以下公式计算:
式(3)中,n为闸孔数目;
b为闸孔净宽,单位m;
H0=H+v0 2/2g为包括行进水头在内的上游总水头,单位m;v0为闸前行进流速;g为重力加速度;
所述的(2)中Q1通过以下公式计算:
2.如权利要求1所述的水力自动翻板闸门流量系数的确定方法,其特征在于所述的试验模型由三角形量水堰、闸前水库库区、水力翻板闸门通过粘接或栓接组装而成;三角形量水堰设于库区上游,水堰岸边翻板闸门前布置水位测针钢尺,水力翻板闸门面板、支腿由有机玻璃制成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010535337.8A CN111680460B (zh) | 2020-06-12 | 2020-06-12 | 水力自动翻板闸门流量系数的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010535337.8A CN111680460B (zh) | 2020-06-12 | 2020-06-12 | 水力自动翻板闸门流量系数的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111680460A CN111680460A (zh) | 2020-09-18 |
CN111680460B true CN111680460B (zh) | 2022-12-02 |
Family
ID=72454735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010535337.8A Active CN111680460B (zh) | 2020-06-12 | 2020-06-12 | 水力自动翻板闸门流量系数的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111680460B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114896566A (zh) * | 2022-05-26 | 2022-08-12 | 长江勘测规划设计研究有限责任公司 | 一种计算闸门过闸流量的智能高效方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393228A (zh) * | 2011-09-24 | 2012-03-28 | 福建四创软件有限公司 | 基于水力学原理推求实时过闸流量应用方法 |
CN106874622A (zh) * | 2017-03-13 | 2017-06-20 | 中国水利水电科学研究院 | 一种串联多渠池闸门过闸流量系数率定方法 |
-
2020
- 2020-06-12 CN CN202010535337.8A patent/CN111680460B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393228A (zh) * | 2011-09-24 | 2012-03-28 | 福建四创软件有限公司 | 基于水力学原理推求实时过闸流量应用方法 |
CN106874622A (zh) * | 2017-03-13 | 2017-06-20 | 中国水利水电科学研究院 | 一种串联多渠池闸门过闸流量系数率定方法 |
Non-Patent Citations (1)
Title |
---|
水力翻板闸门流量系数的试验研究;张月霞等;《水力发电学报》;20101025(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111680460A (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ramamurthy et al. | Characteristics of square-edged and round-nosed broad-crested weirs | |
CN104265652B (zh) | 一种用于大流量低扬程泵站的协同优化设计方法 | |
CN106599374B (zh) | 一种适用于径流潮汐河口段的河相关系及其推导方法 | |
CN111581828B (zh) | 一种感潮河段闸下水位流量关系的计算方法 | |
CN111680460B (zh) | 水力自动翻板闸门流量系数的确定方法 | |
Josiah et al. | An experimental investigation of head loss through trash racks in conveyance systems | |
CN112227292A (zh) | 基于一维水动力模型的空腹拱桥壅水智能分析方法 | |
Xu et al. | Flow pattern and anti-silt measures of straight-edge forebay in large pump stations. | |
CN109063370B (zh) | 漫滩弯曲河道水深平均水流流向预测方法 | |
Wang et al. | The flow patterns and sediment deposition in the forebay of a forward-intake pumping station | |
Ramamurthy et al. | Momentum model of flow past weir | |
Khanahmadi et al. | Hydraulic of curved type-B piano key weirs characteristics under free flow conditions | |
CN103258085A (zh) | 一种污水深海排放扩散器上升管数量确定的方法 | |
US11781942B2 (en) | Method for determining flow velocity distribution in roughness sublayers | |
CN114841020A (zh) | 带导流屏船闸下游引航道口门区横流极值计算方法及应用 | |
CN111121854B (zh) | 测量挑流消能工消能率的装置及使用方法 | |
CN109211603A (zh) | 一种用于模拟气塞形成和积气排出特性的试验装置 | |
CN114757049A (zh) | 流域污水处理厂提标改造的必要性分析验证方法 | |
CN111598757B (zh) | 一种平原河网区城市水系流动确定方法 | |
CN110130471B (zh) | 一种应对泵站前池大扩散角的楔形导流墩及其研究方法 | |
CN208937343U (zh) | 含气塞管道进气和排气模拟试验装置 | |
Pathirana et al. | Discharge coefficient for sharp-crested side weir in supercritical flow | |
CN208167640U (zh) | 掺气水流窄缝消能工水力模型试验装置 | |
Crispino et al. | Junction chamber at vortex drop shaft: case study of Cossonay | |
Angulo et al. | Experimental validation of a low-head turbine intake designed by CFD following Fisher and Franke guidelines |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |