闸孔出流计算程序
闸孔出流计算程序
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
If ht
σ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
σ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