两种不同表达式的水闸淹没宽顶堰流量计算公式计算结相一致的讨论

1、两种不同表达式的水闸淹没宽顶堰流量计算公式计算结相一致的讨论水闸淹没宽顶堰流量计算公式有两种表达形式：一种是 Q=pssBh （0-1）另一种是 Q=O BH03/2 （0-2）式中：Q 过闸流量； 一侧收缩系数；B 一闸室净过水宽度；h一闸室水深，比势能；g 一重力加速度；H。一闸上游翼墙前河道末收缩断面（后简称断面1-1）单位水体总能量；m一流量系数； s-淹没流速系数；o s-淹没流量系数。从水力学知，式（01）是由闸室过水断面（后简称断面2-2）与断面1一1建立能量 关系方程H0=h+ 羹 +& 源 （0-3）_h_推导而来；而式（0 2）又是引进参数K=丑口 ,m*KJl-尺由式（0

2、1）演变而得，二 式同根同源。然而水闸设计工程师都知道，此二公式在相同条件下计算结果是不等的。对某 闸过闸流量核算淹没度hs/h0=0.965,用【参1】按式（01）计算得Q=470m3/s;用【参2】 按式（0-2）计算得Q=394 m3/s，相差近20%。式（01）、（0 2）本同根同源，它们计算结果却不一致，这是不合理的，也不是必 然的。对一个具体的水闸来说，其闸室qh关系曲线只有一条，即在某一水深只能通过一个 流量。一、室矩形过水断面的水力特性特性方程：E=h+ 2g =h+ 2*（1-1）s式中：E一闸室收缩断面22单位水体总能量；v2一平均流速；q 单宽流量；a2动能改正系数；类一

3、断面22比动能;h一同前。式（1 1）即式（03）等号右边的前两项。式（1 1） Es = f（h,q），令q=常量，使其变为平面问题，（如q=5，10，15, 25）可作Eh关系曲线，见图1sI1HL17图1中相应于每一个流量q的曲线就是一条Esh关系曲线。该曲线以横坐标和与横坐 标成45的线oa为渐近线，并且有一断面单位能量最小的点k,该点将曲线分为上、下图1；建形姓水断面$-11曲线翱伽W99Esnh=2.0EE3Tnh=5,1b厘*两支。在下支为急流，Es随h增大而减小；在上支属缓流，Es随h增大而增大。对应 K点的水深即为由急流向缓流转变的水深称其为临界水深hk。故对hk有定义：矩形

4、 过水断面，流量一定，对应于最小断面单位水体能量的水深称为临界水深。令=常 s量（如 Es=2.06，3.26，4.26，5.16，5.99m），可作 qh 关系曲线如图 2。h儿0246010 12 11 16 18 2D 22 2426图？E$定值矩形过水断面q-h曲线km /s)由图2可知，在下支，h增大q随之增大；当q=qmax后，在上支，h增大q随之减小。在q 处只有一个水深与之对应，此水深也定义为临界水深。定义：矩形过水断面， max断面单位水体能量为定值时，对应于最大流量的水深称临界水深。根据上面两定义都可以通过微分方程证明hk=&Es。在图1上也可以看到，连接各流量的K点，可以

5、得到通过坐标原点的直线Ob，直线Ob将45线Oa上任一点的纵坐标分割为1： 2。可见给定Es时，hk= 3 Es。hk的定义与定义是一致的，只是从不同的限定条件即不同的平面角度来表达hk。从图1、图2可以看出式（1 1）的图形是一个以图2曲线为横断面，以图1横坐标（Es）轴为一边，以Oa线为另一边的锥形曲线。断面2 2的淹没流态就是断面22实际进入的缓流区流态。二、式（1-1）可得：(岛-幻V2=Y 的(21)f (Es-h)Q=J% (22)理论上讲，计算淹没宽顶堰的根本公式应是式（2 2），因为水流通过断面2 2 时，必须遵守断面22的水力特性。由于Es及h在设计情况下是未知的，且是互为

6、影响的变量，式（2 2）在设计情况下无法直接使用，需要将其变为只含一个自变量 的公式，建立断面1 1、2 2间的能量关系式（0 3），令 = j彼（1+G整理得:州制 T）V2(23)q=h引类偶T）(24)对比式（2-1）、（23）及式（22）、（24）可见，用H0代替Es后，的参与不 可忽视。式（2 3）也清楚的说明，用用H0代替Es后，根号内计算的速度，必须乘 以修正系数p才能使之等于断面2 2的真实速度V2，所以不能简单地令p = 1.0o水体从断面1 1携带能量H0流到断面2 2,能量要有损失和转换。损失包括沿程损 失（次要的），墩头阻力损失（主要的）水体到达断面2 2后，所带有的能

7、量已经 不是H0而是小于H0的Es。断面1 1、22间的水位差一部分是损失，另一部分转 化为动能。为方便讨论，将断面1 1、22能量关系协成H0=Es+Ahi 2 （2-5）h2 一水流从断面1 1到2 2间的能量损失或称能量差。即（03）式中的。0，当hi-2 = 0，Es=H0,这时断面 22 的 hk=hka=H0。将式（2 1）中的 Es 用 H0 代 之，以h为自变量可作曲线A,见图3。曲线A是断面2 2淹没流的上限，也是一条 只能无限接近的理想曲线。坪MXMA 41）-; 两券*W岳mt .am耳心关方翻d畴I日忍帷& 爰中旅正Eg煞战CfMM 未妁舞f正启拘卦jH ）sa 她觑同

8、岫眺桐心爵幡-琅项虐MWf-449-412用初酷顽|如果在H0 一定的情况下，能通过试验测得断面2 2实际的hk,譬如别列金斯基 试验hk=0.61H0=hkb （任何实测的hk绝对不可能大于或等于H0）。这时断面2 2的Es= 5hkb=0.915H0,h1-2=.085H0。用 Es0 代入式（21），以 h 为自变量可作曲线 B, 见图3。曲线B是断面2 2淹没流淹没开始时刻的h-V关系曲线。三、分析曲线A与曲线B围成的区域可知：h在hkbH0间变化；Es在EsH0之间 变化；V2在V2b0一H0之间变化。在这个淹没流区域里，应该存在一条断面2一 2真实的唯一的V2h关系曲线C。分析曲线

9、C的下端点b0, b0点的纵坐标是实测的hkb，因而它的横坐标V2b0也 应属实测的确定值，所以b0点就是曲线C的起始点。再分析式（21），当Es=H0, hH0时，V20,因此曲线C的上端点趋近于曲线A与纵轴交点a10。曲线C中间各点走向这样确定：将曲线的纵坐标hkb-Es0段（数值等于1/3Es0）分成若 干等分（如10等分），落在曲线B上各点b0、blb10将曲线A纵坐标hka-H0 段（数值等于1/3H0）分成10等分，落在曲线A上各点a0、ala10。连接a0b0,a 1b1,a9b9。将 Es0-H0段（数值等于 0.085H0）纵坐标分成 10 等分，即有Es1,Es2.。 用E

10、s1,Es2再作九条V2h关系曲线。这样就把曲线A与曲线B围成的区域打成网格。从b0点开始，连接Es1与albl的交点，连接Es2与a2b2的交点.直到a10, 曲线C即作成。利用V+AV=| 修 + 敬)- (h + Ah)列表计算一样可以作出曲线C。曲线C体现随h增加，Es增加一个量级，速度V2减小一个量级，每次减小不是 一个常数。分析图1可见，在缓流区，随h增加Es增加，但动能减小，即V2减小。 曲线C符合这一规律。由于曲线C首尾是确定的，中间各点由网格法确定，所以整个 曲线是唯一的。根据q=hV，同样可以作图4曲线C。A：阻了41券并B： Ah 1-2=0.085甘时关彖曲线 笑的营储

11、系由皇& W（M）计费中时关系由姬”国 4S2-A/心Ma 话 *网 玷 的 sb4所昵-己时摆W成曲线C在图1、图2上的表示即曲线MN。四、式（21）是曲线C，式（2 3）根号部分是曲线A。要使式（23）根号 部分的计算结果符合式（21），就要对其加以修正。在同一水深曲线A上的流速V 2a需要乘一个小于1.0的p值，使之等于断面22的实际流速值即曲线C上的速度全V2c，所以咔矿弛。在图3上，利用各淹没度时的水深的V2c、V2a即可得到与淹没冬岛F 竺上七度对应的p值。另据V2=知=，令g=。，K=可得p=V知“,取不同K值、g值即可得p值。用此式与前图解求得p值结果是一致的。同样利用图4曲线

12、C上各淹没度对应水深的qi与b0点的qb0的比值。=，可求出与 各淹没度对应的。值。根据hkb = 0.61H0,结合实历本文求得p值、。值如表1。表1中ps、os为式（01）、（0 2）用值。表1hs0.800.820.840.860.880.900.920.940.960.980.990.880.880.880.880.880.880.890.890.890.890.89456789012341.001.011.000.990.970.940.900.850.790.720.6801831004555如1.000.960.950.940.930.920.920.910.900.900.90

13、074468159311.000.990.970.950.900.840.780.700.590.400.28表1中当=0.82、0.84时，o值大于1.0，这与前述hk的定义（2）没有矛盾。Hk 的定义（2）是在Es为定值，只有h 一个变量的情况。曲线C是q=f（Es,h）,Es、h均 为变量的情况。五、公式（01）、（02）存在的问题式（01）、（0 2）存在的问题总的来说就是qh关系曲线的起始点都在曲线A上。这意味着在临界水深hkb时,h1-2=0, p=1.0,速度没有修正，倒是h增加，V2减小反而给予修正。实例：H0=4.92m, hkb=0.61H0=3.00m,Es=4.50。临

14、界水深hkb时的真实流量应是b0点 的流量 qb=3.0 x4.43x=16.28m3/s。式（0 1）计算水深 hkb 时的流量 qka=1.0 x 3.0 x4.43xJ492-3.=18.415 m3/s;式 （0 1）计算水深 hka 时的流量 qka= 1.0 x3.28x4.43xJ4.92-3.28 =18.61 m3/s。在曲线A上论临界水深是hka,这时流量最大； 当hhka,随 h减小q减小。hka hka,所以qka qka，但qka仍然属于曲线A向下延伸线上的点。比较图4曲线D与曲线C可见，曲线D计算结果一直偏大，但在高淹没度二者非常 接近。式（0 2）的计算思想是：以

15、临界水深的流量为基准，其上各淹没度对应水深时的流 量是用临界水深时的流量乘一个比值系数os而得。M*k,k=。以平底闸m=0.385看， 显然是令p=1.0，h=hka而得。就是说m底 H03/2是曲线A的下端点a0点的流量。hs实例 qka=osm H03/2= 1.0 x0.385x4.43x10.913= 18.61 m3/s。当丑 =0.80，os=1.0，用式（02）计算的qh关系E的下端点是a0点的流量垂直下移至hkb的高度。 这显然也是不经修正的流量。令m =k亦耳，取不同k值可得m =f（k）曲线，该曲 线形状同式（2 2）曲线形状完全一样。当k=时有m max=0.385,因

16、此可以说m的 物理意义是H0=1.0时的流量。p值前已论述，在淹没流开始时，p不可能是1.0，所 以m也不应是0.385。hs见图4,对比曲线E与曲线C，在丑0=0.92时，二曲线相交，其他各处相差甚大。 在交点以下式（0 2）计算结果偏大；交点以上式（0 2）计算结果偏小。六、修正公式的表达式以曲线C与曲线A的关系为基础，建立五个公式表达式如下: TOC o 1-5 h z 1、Q= Bhp 独（丑口 -幻（6-1）式中：一流速改正系数，见表1。此式表达与式（0-1 ）同，但p与桥有明显不同。2、 Q=o E Bh/jgH （6-2）式中：a 一淹没流量系数见表1。hkb=0.61H0, %

17、=0.9565, pk意义是将图3中a0 点的速度修正为b0点的速度。3、 Q=a E BhkJ （6-3）式中：#=0.915,氏 H0=Es0。4、 Q=a E Bm *（声 H0）3/2 （6-4）5、Q=a E Bm(6-5)式中：m为不同坎高（丑）堰顶临界水深时的最大流量系数。打=0.8748,的意义是将图 4中a0点的流量修正为b0点的流量。本文认为E在淹没流过程是个变量，隐含在损失中。篇幅所限，此处按传统写法仍写于此。现以实例取式（61）、（65）与式（01）、（02）计算结果进行比较，见表2。从 表2可见修正公式计算结果是一致的。表2Hs/H00.800.820.840.860

18、.880.900.920.940.960.980.99h3.003.253.463.643.803.954.094.224.354.474.53公式(0 1)700.0683.7671.4654.4633.6607.8577.7543.8502.5455.7429.0公式(0 2)707.1700.0685.9671.8636.4594.0551.5495.0417.2282.8198.0公式(6 1)618.6625.7623.6614.9601.1582.2558.3529.5493.2450.8425.8公式(6 5)618.6625.5623.6614.3600.7581.5556.8528.3491.8448.5423.8主要参考书1、水闸设计谈松曦水利电力出版社1986年第一版2、水闸设计规范（试行）SD133-863、水力学华东水利学院教研室黄文煌主编人民教育出