哈工大流体力学大三下(上机)报告.docVIP

工程流体力学（上级）报告 题目：二管串联管路的水力计算 学院： 班级： 姓名： 学号： 二管串联的管路水力计算分析 已知:h, l1, l2, d1, d2, Δ1, Δ2, ν, ζ, 求:qv. 用vb编程如下： 1程序代码 ’编写公共函数 Public Function jslmd(v As Double, niu As Double, d As Double, dlt As Double) As Double re1 = 2320 re2 = 80 * d / dlt re3 = 4160 * (d / (2 * dlt)) ^ 0.85 re = v * d / niu If re re1 Then lmd = 64 / re id = 1 ElseIf re re2 Then If re 100000 Then lmd = 0.3164 * re ^ (-0.25) id = 2 ElseIf re 3000000 Then lmd = 0.0032 + 0.221 * re ^ (-0.237) id = 3 End If ElseIf re re3 Then x0 = 10 Call ddf(dlt, d, re, x0, 0.000001, x) lmd = x ^ (-2) id = 4 Else lmd = (1.74 + 0.434294481 * 2 * Log(d * 0.5 / dlt)) ^ (-2) id = 5 End If jslmd = lmd End Function Public Sub ddf(ByVal dlt, ByVal d, ByVal re, ByVal x0, ByVal eps, x) Do x = -2 * 0.434294481 * Log(dlt / (3.7 * d) + 2.51 * (x0) / re) If Abs(x - x0) eps Then Exit Do x0 = x Loop End Sub Private Sub Command1_Click() ’定义变量 Dim h As Double, niu As Double, s As Double, l1 As Double, d1 As Double, dlt1 As Double, l2 As Double, d2 As Double, dlt2 As Double, v1 As Double, v2 As Double, q As Double, lmd1 As Double, lmd2 As Double, lmd11 As Double, lmd22 As Double Const pi = 3.1415927 Dim id As Integer ’输入已知参数 h = Val(Text1.Text) niu = Val(Text2.Text) s = Val(Text3.Text) l1 = Val(Text4.Text) d1 = Val(Text5.Text) dlt1 = Val(Text6.Text) l2 = Val(Text7.Text) d2 = Val(Text8.Text) dlt2 = Val(Text9.Text) lmd11 = 0.025 lmd22 = 0.015 ’执行操作 Do v1 = (2 * 9.8 * h) ^ 0.5 * (s + lmd11 * l1 / d1 + (1 - (d1 / d2) ^ 2) ^ 2 + (lmd22 * l2 / d2) * ((d1 / d2) ^ 4 + (d1 / d2) ^ 4)) ^ (-0.5) v2 = v1 * (d1 / d2) ^ 2 lmd1 = jslmd(v1, niu, d1, dlt1) lmd2 = jslmd(v2, niu, d2, dlt2) If Abs(lmd11 - lmd1) 0.000001 And Abs(lmd22 - lmd2) 0.000001 Then Exit Do lmd11 = lmd1 lmd22 = lmd2 Loop ’输出结果 q