闸孔出流计算程序

For personal use only in study and research; not for commercial use

闸孔出流计算程序

Function lq(h0, h2, n, e) 'h0上游水位，h2下游水位，n开闸孔数，e闸孔开度

zmk = 14 '闸门的宽度

dbgc = 1.5 '底板高程

eopen = e * 0.2 '闸门开档数×每档0.2米

h = h0 - 1.5 '上游水深

ht = h2 - 1.5 '下游水深

eh = eopen / h 'e/h

μ = (0.6 - 0.176 * eh) '流量系数

ε2 = 0.3568 * eh ^ 3 - 0.2043 * eh ^ 2 + 0.0822 * eh + 0.6088 '垂直收缩系数

hc = ε2 * eopen '收缩水深

vc = μ / ε2 * (2 * 9.81 * h) ^ 0.5 '收缩断面流速

hcc = hc / 2 * ((1 + 8 * vc ^ 2 / (9.81 * hc)) ^ 0.5 - 1) '收缩断面的共轭水深

qlb = (ht - hcc) / (h - hcc) '潜流比

If eh <= 0.65 Then '闸孔出流

If ht <= hcc Then '闸孔自由出流

σs = 1

Else ' 闸孔淹没出流

σs = -6.0158 * qlb ^ 5 + 14.425 * qlb ^ 4 - 13.59 * qlb ^ 3 + 6.0832 * qlb ^ 2 - 1.8909 * qlb + 1.0002

'σs闸孔淹没系数

End If

lq = σs * μ * zmk * n * eopen * (2 * 9.81 * h) ^ 0.5

Else '堰流状态

If ht / h >= 0.9 Then '堰流高淹没状态

μ0 = 0.877 + (ht / h - 0.65) ^ 2

lq = μ0 * zmk * n * ht * (2 * 9.81 * (h - ht)) ^ 0.5

Else '堰流低淹没状态

If ht / h <= 0.8 Then

σs = 1 '堰流淹没系数

Else

σs = -96.007 * (ht / h) ^ 3 + 235.32 * (ht / h) ^ 2 - 193.23 * (ht / h) + 54.137

End If

lq = σs * 0.95 * 0.385 * 14 * n * (2 * 9.81 * h ^ 3) ^ 0.5 '0.95：侧收缩系数，0.385：堰流流量系数

End If

End If

End Function

仅供个人用于学习、研究；不得用于商业用途。

For personal use only in study and research; not for commercial use.

Nur für den persönlichen für Studien, Forschung, zu kommerziellen Zwecken verwendet werden.

Pour l 'étude et la recherche uniquement à des fins personnelles; pas à des fins commerciales.

только для людей, которые используются для обучения, исследований и не должны использоваться в коммерческих целях.

以下无正文