CN113255030A - 一种矩形断面明渠水力最优断面设计方法 - Google Patents
一种矩形断面明渠水力最优断面设计方法 Download PDFInfo
- Publication number
- CN113255030A CN113255030A CN202110459238.0A CN202110459238A CN113255030A CN 113255030 A CN113255030 A CN 113255030A CN 202110459238 A CN202110459238 A CN 202110459238A CN 113255030 A CN113255030 A CN 113255030A
- Authority
- CN
- China
- Prior art keywords
- section
- channel
- sand
- water
- hydraulic
- 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
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Geometry (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Architecture (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Underground Or Underwater Handling Of Building Materials (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种矩形断面明渠水力最优断面设计方法,以最大流量模数原理为理论基础,在计算过程中以水流连续性方程、水流阻力方程以及推移质输沙率为控制方程,通过迭代计算的方法,综合考虑了矩形断面明渠内泥沙输运与床面变形对渠道过水断面大小与形状的影响,通过引入水力学中的物理量—流量模数K,全面考虑实际渠道内水流阻力大小对其水流通过能力的影响,是一种基于最大流量模数原理的矩形断面人工明渠水力最优断面设计方法,解决了以往设计方法中床面粗糙程度取值随意以及无法考虑水流‑泥沙‑河床相互作用的问题,弥补了人工明渠水力最优断面设计方法的理论缺陷。
Description
技术领域
本发明涉及一种明渠断面设计方法,特别是一种矩形断面明渠水力最优断面设计方法。
背景技术
明渠是人工修建或自然形成的具有自由表面的渠道。对于自然形成的明渠,其横断面往往分为主槽和滩地,呈不规则形状。人工修建的渠道,其断面往往根据地质条件来确定。土质地基在水利工程实践中最为常见,且考虑到渠道岸坡的稳定性与施工工艺限制,土质地基渠道的断面形状一般不会做成圆形或半圆形,多为矩形或梯形,后续又逐渐发展出U型断面以及抛物型断面等。
在进行明渠水力计算时,当流量Q、比降S以及床面曼宁系数n已知的情况下,进行合理设计可以得到最小的过水断面面积A,或者当比降S以及床面曼宁系数n已知的情况下,对于确定的过水断面面积A,其通过的流量Q能够达到最大,达到减小工程施工量的目的,此时的过水断面就称为水力最优断面或水力最经济断面。迄今为止,国内外学者对于明渠水力最优断面已经做过许多研究,但成果较为零散,应用的理论基础也不尽相同,尚未具有普遍性的指导意义。尤其是在以往的研究中,均未考虑渠道内输沙强度与床面形态对水流通过能力的影响。实际上渠道内的水流、泥沙与床面三者之间相互作用相互影响,水流受到的阻力分布十分复杂,假定渠道过水断面中各处的床面曼宁系数n处处相等是不符合实际情况的。
发明内容
本发明为解决公知技术中存在的技术问题而提供一种矩形断面明渠水力最优断面设计方法,该方法全面考虑实际渠道内水流阻力大小对其通过能力的影响。
本发明为解决公知技术中存在的技术问题所采取的技术方案是:一种矩形断面明渠水力最优断面设计方法,利用迭代法采用计算机程序进行计算,采用以下步骤:
S1:根据研究渠道的来水来沙条件以及推移质级配分布,确定渠道内实际的水流流量Q,断面输沙率Qs和床沙中值粒径d50;
S2:确定假定的渠道平均宽度B;
S3:确定假定的渠道平均水深H;
S4:根据渠道断面的几何形态关系求解水力半径R,计算公式如下:
R=BH/(B+2H) (1)
S5:根据水流连续方程求解平均流速U,计算公式如下:
U=Q/(BH) (2)
S6:在公式(3)中输入假定的与沙粒阻力对应的水力半径R′,求解水面比降S:
S7:求解与沙粒阻力对应的无量纲剪切应力Θ′,并判断其值是否满足1>Θ′>0.06,若不满足,返回步骤S6;若满足,则进入步骤S8,Θ′的计算公式如下:
其中,γs与γ分别为泥沙颗粒与水的容重,取值为2650kg/m3和1000kg/m3;
S8:根据步骤S7求解得到的Θ′,判断其值与0.55的相对大小,然后根据公式(5),求解床面上总的无量纲剪切应力Θ,Θ的计算公式如下:
S9:根据步骤S8中求解得到的Θ计算水力半径R″,判断其是否等于R,若不相等,返回步骤S6;若相等,则进入步骤S10;R″的计算公式如下:
S10:根据Meyer-Petter和Müller输沙率公式计算推移质输沙率Q′s,并判断其是否等于Qs,若不相等,返回步骤S3;若相等,则进入步骤11,Meyer-Petter和Müller输沙率公式如下:
S11:输出B、K、H、λ;其中,断面宽深比λ、流量模数K的计算公式分别如下:
λ=B/H (8)
S12:重复步骤S2~S11,得到特定来水来沙条件以及推移质级配分布条件下的一系列B、K、H、λ;
S13:找出步骤S12中输出的流量模数K的最大值Kmax,以及与Kmax对应的B0、H0、λ0;B0、H0、λ0即为在特定来水来沙条件以及推移质级配分布条件下水力最优断面的水面平均宽度、平均水深和断面宽深比。
进一步地,所述步骤S2和S3中的B和H取值范围为10-1m~103m,迭代计算过程中二者均以10-2m为增量从最小值10-1m开始取值计算。
进一步地,所述步骤S6中的R′取值范围为10-3m~103m,迭代计算过程中以10-4m为增量从最小值10-3m开始取值计算。
本发明具有的优点和积极效果是:以最大流量模数原理为理论基础,在计算过程中以水流连续性方程、水流阻力方程以及推移质输沙率为控制方程,通过迭代计算的方法,综合考虑了矩形断面明渠内泥沙输运与床面变形对渠道过水断面大小与形状的影响,解决了以往设计方法中床面粗糙程度取值随意以及无法考虑水流-泥沙-河床相互作用的问题,弥补了人工明渠水力最优断面设计方法的理论缺陷。
附图说明
图1为本发明的流程图;
图2为流量Q=100m3/s,输沙率Qs=0.002m3/s、0.004m3/s、0.005m3/s和0.006m3/s,d50=0.3mm四种工况条件下水面平均宽度B~平均水深H分布曲线;
图3为流量Q=100m3/s,输沙率Qs=0.002m3/s、0.004m3/s、0.005m3/s和0.006m3/s,d50=0.3mm四种工况条件下水面平均宽度B~流量模数K分布曲线。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下:
请参阅图1,一种矩形断面明渠水力最优断面设计方法,利用迭代法采用计算机程序进行计算,包括以下步骤:
S1:根据研究渠道的来水来沙条件以及推移质级配分布,确定渠道内实际的水流流量Q,断面输沙率Qs和床沙中值粒径d50。
S2:确定假定的渠道平均宽度B。
S3:确定假定的渠道平均水深H。
S4:根据渠道断面的几何形态关系求解水力半径R。对于矩形断面明渠,水力半径R的计算公式如下:
R=BH/(B+2H) (1)
S5:根据水流连续方程求解平均流速U。对于矩形断面明渠,平均流速U的计算公式如下:
U=Q/(BH) (2)
S6:在公式(3)中输入假定的与沙粒阻力对应的水力半径R′,求解水面比降S。
S7:求解与沙粒阻力对应的无量纲剪切应力Θ′,并判断其值是否满足1>Θ′>0.06,若不满足,返回步骤S6;若满足,则进入步骤S8。Θ′的计算公式如下:
其中,γs与γ分别为泥沙颗粒与水的容重,一般取值为2650kg/m3和1000kg/m3。
S8:根据上一步骤S7中求解得到的Θ′,判断其值与0.55的相对大小,然后根据公式(5),求解床面上的总的无量纲剪切应力Θ,Θ的计算公式如下:
S9:根据上一步骤S8中求解得到的Θ计算水力半径R″,判断其是否等于R,若不相等,返回步骤S6;若相等,则进入步骤S10。R″的计算公式如下:
S10:根据Meyer-Petter和Müller输沙率公式计算推移质输沙率Q′s,并判断其是否等于Qs,若不相等,返回步骤S3;若相等,则进入步骤11。Meyer-Petter和Müller输沙率公式如下:
S11:输出B、K、H、λ;其中,断面宽深比λ、流量模数K的计算公式分别如下:
λ=B/H (8)
S12:重复步骤S2~S11,得到特定来水来沙条件以及推移质级配分布条件下的一系列B、K、H、λ。
S13:找出上一步骤S12中输出的流量模数K的最大值Kmax,以及与Kmax对应的B0、H0、λ0;B0、H0、λ0即为在特定来水来沙条件以及推移质级配分布条件下水力最优断面的水面平均宽度、平均水深和断面宽深比。
在本实施例中,所述步骤S2和S3中的B和H取值范围为10-1m~103m,迭代计算过程中二者均以10-2m为增量从最小值10-1m开始取值计算。所述步骤S6中的R′取值范围为10-3m~103m,迭代计算过程中以10-4m为增量从最小值10-3m开始取值计算。
上述设计方法考虑了渠道内泥沙输运和床面形态变化对渠道水流通过能力的影响,通过引入水力学中的物理量—流量模数K,全面考虑实际渠道内水流阻力大小对其水流通过能力的影响,是一种基于最大流量模数原理的矩形断面人工明渠水力最优断面设计方法。
请参见图2和图3,图2和图3分别为在流量Q=100m3/s,输沙率Qs=0.002m3/s、0.004m3/s、0.005m3/s和0.006m3/s,d50=0.3mm四种工况条件下,采用本发明所述设计方法,计算得到B~H、B~K分布曲线。由图2可以发现,在四种来水来沙条件以及推移质级配分布条件下,满足条件的B和H的组合均有多组,且B和H均是一一对应呈反比的关系。由图3可以发现,尽管满足特定来水来沙条件和推移质级配分布条件的B与H有多组,但有且仅有一组B与H的组合能够使得K=Kmax,即流量模数K达到最大值。按照最大流量模数原理,与该组B与H对应的过水断面即为研究渠道的水力最优断面。
尽管上面结合附图对本发明的优选实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,并不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围的情况下,还可以做出很多形式,这些均属于本发明的保护范围之内。
Claims (3)
1.一种矩形断面明渠水力最优断面设计方法,其特征在于,利用迭代法采用计算机程序进行计算,采用以下步骤:
S1:根据研究渠道的来水来沙条件以及推移质级配分布,确定渠道内实际的水流流量Q,断面输沙率Qs和床沙中值粒径d50;
S2:确定假定的渠道平均宽度B;
S3:确定假定的渠道平均水深H;
S4:根据渠道断面的几何形态关系求解水力半径R,计算公式如下:
R=BH/(B+2H) (1)
S5:根据水流连续方程求解平均流速U,计算公式如下:
U=Q/(BH) (2)
S6:在公式(3)中输入假定的与沙粒阻力对应的水力半径R′,求解水面比降S:
S7:求解与沙粒阻力对应的无量纲剪切应力Θ′,并判断其值是否满足1>Θ′>0.06,若不满足,返回步骤S6;若满足,则进入步骤S8,Θ′的计算公式如下:
其中,γs与γ分别为泥沙颗粒与水的容重,取值为2650kg/m3和1000kg/m3;
S8:根据步骤S7求解得到的Θ′,判断其值与0.55的相对大小,然后根据公式(5),求解床面上总的无量纲剪切应力Θ,Θ的计算公式如下:
S9:根据步骤S8中求解得到的Θ计算水力半径R″,判断其是否等于R,若不相等,返回步骤S6;若相等,则进入步骤S10;R″的计算公式如下:
S10:根据Meyer-Petter和Müller输沙率公式计算推移质输沙率Qs′,并判断其是否等于Qs,若不相等,返回步骤S3;若相等,则进入步骤11,Meyer-Petter和Müller输沙率公式如下:
S11:输出B、K、H、λ;;其中,断面宽深比λ、流量模数K的计算公式分别如下:
λ=B/H (8)
S12:重复步骤S2~S11,得到特定来水来沙条件以及推移质级配分布条件下的一系列B、K、H、λ;
S13:找出步骤S12中输出的流量模数K的最大值Kmax,以及与Kmax对应的B0、H0、λ0;B0、H0、λ0即为在特定来水来沙条件以及推移质级配分布条件下水力最优断面的水面平均宽度、平均水深和断面宽深比。
2.根据权利要求1所述的矩形断面明渠水力最优断面设计方法,其特征在于,所述步骤S2和S3中的B和H取值范围为10-1m~103m,迭代计算过程中二者均以10-2m为增量从最小值10-1m开始取值计算。
3.根据权利要求1所述的矩形断面明渠水力最优断面设计方法,其特征在于,所述步骤S6中的R′取值范围为10-3m~103m,迭代计算过程中以10-4m为增量从最小值10-3m开始取值计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110459238.0A CN113255030B (zh) | 2021-04-27 | 2021-04-27 | 一种矩形断面明渠水力最优断面设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110459238.0A CN113255030B (zh) | 2021-04-27 | 2021-04-27 | 一种矩形断面明渠水力最优断面设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113255030A true CN113255030A (zh) | 2021-08-13 |
CN113255030B CN113255030B (zh) | 2022-05-06 |
Family
ID=77221846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110459238.0A Active CN113255030B (zh) | 2021-04-27 | 2021-04-27 | 一种矩形断面明渠水力最优断面设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113255030B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7039565B1 (en) * | 2003-01-03 | 2006-05-02 | Haestad Methods, Inc. | Method and system for developing a numerical dynamic sanitary sewer and storm water drainage simulation model |
CN101418564A (zh) * | 2008-12-10 | 2009-04-29 | 天津大学 | 一种稳定输水渠道 |
CN104535123A (zh) * | 2014-12-29 | 2015-04-22 | 太原理工大学 | 在线式渠道流量自动测量方法 |
CN104765977A (zh) * | 2015-04-27 | 2015-07-08 | 河海大学 | 考虑局部水头损失的河段综合糙率的计算方法 |
JP6047738B1 (ja) * | 2016-08-16 | 2016-12-21 | 有限会社北沢技術事務所 | 開水路流量計測装置 |
CN106503344A (zh) * | 2016-10-31 | 2017-03-15 | 济南大学 | 一种冰盖下明渠输水渠道及其水力最优断面的求解方法 |
CN106595777A (zh) * | 2016-12-01 | 2017-04-26 | 广西师范大学 | 一种非接触式探测河流断面流量的计算方法 |
CN107480384A (zh) * | 2017-08-24 | 2017-12-15 | 北方民族大学 | 河道水流泥沙平面二维数值模拟方法和系统 |
CN108593023A (zh) * | 2018-02-21 | 2018-09-28 | 李银才 | 利用雷达测速仪自动测量明渠流量的方法 |
KR20180115631A (ko) * | 2017-04-13 | 2018-10-23 | 부산대학교 산학협력단 | 유체시스템 최적설계방법 |
CN109165402A (zh) * | 2018-07-11 | 2019-01-08 | 济南大学 | 一种求解通用幂函数形明渠水力最优断面的方法 |
CN109255097A (zh) * | 2018-07-11 | 2019-01-22 | 济南大学 | 一种3.3471次方幂函数形明渠输水断面及最优断面求解方法 |
CN110008513A (zh) * | 2019-03-06 | 2019-07-12 | 济南大学 | 一种新型平底悬链线形明渠及其水力最优断面的求解方法 |
CN111400974A (zh) * | 2020-04-27 | 2020-07-10 | 中国水利水电科学研究院 | 一种估算矩形河渠壁面和床面切应力的方法 |
-
2021
- 2021-04-27 CN CN202110459238.0A patent/CN113255030B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7039565B1 (en) * | 2003-01-03 | 2006-05-02 | Haestad Methods, Inc. | Method and system for developing a numerical dynamic sanitary sewer and storm water drainage simulation model |
CN101418564A (zh) * | 2008-12-10 | 2009-04-29 | 天津大学 | 一种稳定输水渠道 |
CN104535123A (zh) * | 2014-12-29 | 2015-04-22 | 太原理工大学 | 在线式渠道流量自动测量方法 |
CN104765977A (zh) * | 2015-04-27 | 2015-07-08 | 河海大学 | 考虑局部水头损失的河段综合糙率的计算方法 |
JP6047738B1 (ja) * | 2016-08-16 | 2016-12-21 | 有限会社北沢技術事務所 | 開水路流量計測装置 |
CN106503344A (zh) * | 2016-10-31 | 2017-03-15 | 济南大学 | 一种冰盖下明渠输水渠道及其水力最优断面的求解方法 |
CN106595777A (zh) * | 2016-12-01 | 2017-04-26 | 广西师范大学 | 一种非接触式探测河流断面流量的计算方法 |
KR20180115631A (ko) * | 2017-04-13 | 2018-10-23 | 부산대학교 산학협력단 | 유체시스템 최적설계방법 |
CN107480384A (zh) * | 2017-08-24 | 2017-12-15 | 北方民族大学 | 河道水流泥沙平面二维数值模拟方法和系统 |
CN108593023A (zh) * | 2018-02-21 | 2018-09-28 | 李银才 | 利用雷达测速仪自动测量明渠流量的方法 |
CN109165402A (zh) * | 2018-07-11 | 2019-01-08 | 济南大学 | 一种求解通用幂函数形明渠水力最优断面的方法 |
CN109255097A (zh) * | 2018-07-11 | 2019-01-22 | 济南大学 | 一种3.3471次方幂函数形明渠输水断面及最优断面求解方法 |
CN110008513A (zh) * | 2019-03-06 | 2019-07-12 | 济南大学 | 一种新型平底悬链线形明渠及其水力最优断面的求解方法 |
CN111400974A (zh) * | 2020-04-27 | 2020-07-10 | 中国水利水电科学研究院 | 一种估算矩形河渠壁面和床面切应力的方法 |
Non-Patent Citations (2)
Title |
---|
张玮等: ""基于最小能耗原理的明渠水利最佳断面数值研究"", 《科学技术与工程》 * |
闵志华: ""山区小型灌溉渠道横断面的宽深比分析"", 《浙江水利水电学院学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113255030B (zh) | 2022-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Annandale | Reservoir sedimentation | |
CN102002925B (zh) | 多沙河流河工动床模型人工转折方法及人工转折导流槽 | |
CN108629055B (zh) | 一种基于饱和输沙原理的沙质内河航道回淤量预报方法 | |
CN112989565B (zh) | 一种基于水流挟沙能力的造床流量计算方法 | |
Wormleaton et al. | Flow structures in a two-stage channel with a mobile bed | |
CN106908349B (zh) | 确定筑坝原级配砂砾料最大干密度的方法 | |
Kalathil et al. | Control of sediment inflow into a trapezoidal intake canal using submerged vanes | |
CN113255030B (zh) | 一种矩形断面明渠水力最优断面设计方法 | |
CN106320255B (zh) | 一种漫滩水流滩槽过流量的计算方法 | |
CN104195979A (zh) | 河道交汇水流停滞区楔锥体及其构建方法与应用 | |
CN111241467A (zh) | 一种多沙河流水库排沙比的计算方法 | |
Moghadam et al. | Numerical study of energy dissipation and block barriers in stepped spillways | |
CN116084337B (zh) | 用于产卵场修复中宽级配河床糙率确定的水槽试验方法 | |
CN103392630B (zh) | 用于形成有利于中华鲟产卵地形的结构及底质优化方法 | |
Eshev et al. | Calculation of parameters of subsurface ridges in a steady flow of groundwater channels | |
Gupta et al. | Design of minimum cost earthen channels having side slopes riveted with different types of riprap stones and unlined bed by using particle swarm optimization | |
Singh et al. | Optimal design of hydraulic structures with hybrid differential evolution multiple particle swarm optimization | |
CN111666620B (zh) | 一种颗粒材料侧向压力分布特征定量描述方法 | |
Sobir et al. | Calculation of parameters of subsurface ridges in a steady flow of groundwater channels | |
CN117874897B (zh) | 一种堰塞体高度预测方法 | |
CN118070686B (zh) | 一种基于机器学习构造河流床沙级配分布的方法 | |
Elsayed et al. | Comparative study of different scenarios for the morphological evolution in a river stream | |
Burele et al. | Numerical Simulation-Optimization for Channelization of River Kosi | |
CN108755785B (zh) | 一种阶梯形变径空心桩内力计算方法 | |
Prasetyo et al. | Analysis on The Effect of Groyne Type Impermeable Placement on Sediment Distribution in Lariang River Bend |
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 |