水力学水面曲线计算VB代码.doc

程序共有2个文件，Hydraulics.bas和Form1.frm。Hydraulics.bas文件里定义了水力学计算所需的相关函数，后面的计算全部基于这些函数。相关函数有：（参数已省略）author:吴昱驹e-mail:328522073platform: VB6.0fxA()计算过水面积fxv()计算断面平均流速fxEs()计算断面单位能量fxv_average()计算平均速度fxR_average()计算平均水力半径fxC_average()计算平均谢才系数fxX()计算湿周fxR()计算水利半径fxC()计算谢才系数fxJ()计算水力坡度fxdeltaS()计算SfxHc()计算临界水深fxH0()计算正常水深Hydraulics文件内容如下Public Const g As Single = 9.81定义重力加速度常量gFunction fxA(b As Single, m As Single, h As Single) As Single 计算过水面积 fxA = (b + m * h) * hEnd FunctionFunction fxv(q As Single, a As Single) As Single 计算断面平均流速 fxv = q / aEnd FunctionFunction fxEs(h As Single, alpha As Single, v As Single) As Single 计算断面单位能量 fxEs = h + alpha * v * v / (2 * g)End FunctionFunction fxv_average(v1 As Single, v2 As Single) As Single 计算平均速度 fxv_average = 0.5 * (v1 + v2)End FunctionFunction fxR_average(r1 As Single, r2 As Single) As Single 计算平均水力半径 fxR_average = 0.5 * (r1 + r2)End FunctionFunction fxC_average(c1 As Single, c2 As Single) As Single 计算平均谢才系数 fxC_average = 0.5 * (c1 + c2)End FunctionFunction fxX(b As Single, h As Single, m As Single) As Single 计算湿周 fxX = b + 2 * h * Sqr(1 + m * m)End FunctionFunction fxR(a As Single, x As Single) As Single 计算水力半径 fxR = a / xEnd FunctionFunction fxC(n As Single, r As Single) As Single 计算谢才系数 fxC = (r (1 / 6) / nEnd FunctionFunction fxJ(v As Single, c As Single, r As Single) As Single 计算水力坡度 fxJ = v * v / (c * c * r)End FunctionFunction fxdeltaS(es1 As Single, es2 As Single, i As Single, J As Single) As Single 计算S fxdeltaS = (es2 - es1) / (i - J)End FunctionFunction fxHc(alpha As Single, q As Single, htest As Single, b As Single, m As Single) As Single 试算法计算临界水深 Dim temp As Single, ht As Single, temph As Single temp = alpha * q * q / g For ht = htest To 0 Step -0.0001 temph = (b + m * ht) * ht) 3 / (b + 2 * m * ht) If Abs(temph - temp) 0 Then fxHc = ht Else fxHc = -1 End IfEnd FunctionFunction fxH0(q As Single, b As Single, m As Single, n As Single, i As Single) As Single 计算正常水深 Dim qt As Single, ht As Single, c As Single, a As Single, r As Single, x As Single For ht = 0 To 10 Step 0.0001 x = fxX(b, ht, m) a = fxA(b, m, ht) r = fxR(a, x) c = fxC(n, r) qt = c * a * Sqr(r * i) If Abs(qt - q) 0.1 Then Exit For End If Next ht If ht H0 Then h = h + s End If Next ii For ii = 2 To 7 v_ave = fxv_average(v(ii), v(ii - 1) C_ave = fxC_average(c(ii), c(ii - 1) R_ave = fxR_average(r(ii), r(ii - 1) J(ii) = fxJ(v_ave, C_ave, R_ave) deltaS(ii) = fxdeltaS(Es(ii), Es(ii - 1), i, J(ii) Next ii h = h1 For ii = 1 To 7 If ii = 1 Then Picture1.Print h X R C v Es S End If If deltaS(ii) 0 Then Picture1.Print Format(h, 0.00 ); Format(x(ii), 0.00 ); Format(r(ii), 0.00 ); Format(c(ii), 0.00 ); Format(v(ii), 0.000 ); Format(Es(ii), 0.000 ); Format(deltaS(ii), 0) Else Picture1.Print Format(h, 0.00 ); Format(x(ii), 0.00 ); Format(r(ii), 0.00 ); Format(c(ii), 0.00 ); Format(v(ii), 0.000 ); Format(Es(ii), 0.000 ) & - End If If h H0 Then h = h + s End If Next ii Dim sumS sumS = 0 For ii = 1 To 7 If deltaS(ii) 0 Then sumS = sumS + deltaS(ii) End If Next ii Picture1.Print Picture1.Print S= & CStr(sumS)End SubPrivate Sub Command2_Click() Picture1.Cls Picture2.ClsEnd SubPrivate Sub Form_Load() b = CSng(Text1.Text) m = CSng(Text2.Text) n = CSng(Text3.Text) q = CSng(Text4.Text) h1 = CSng(Text5.Text) alpha = CSng(Text6.Text) i = CSng(Text7.Text) L = CSng(Text8.Text) s = CSng(Text9.Text)End SubPrivate Sub printHcH0() If Hc = -1 Then Picture2.Print 临界水深5，超出计算范围 Else Picture2.Print 临界水深: Picture2.Print Hc= & CStr(Hc) End If If H0 = -1 Then Picture2.Print