【正文】 e Sub Cal_D_Click() = ( ) / 2 = ( ) / 2 = 下游坡面計算 = 下游水位H2(m) = 下游壩面坡度m2End SubPrivate Sub Doc_Click() = ( ) / 2 = ( ) / 2End SubPrivate Sub Exit_Click() EndEnd SubPrivate Sub Form_Load() = ( ) / 2 = ( ) / 2ed = 1End SubPrivate Sub Option_Click() = ( ) / 2 = ( ) / 2End SubOption ExplicitPrivate Sub Command1_Click() If = Then MsgBox 請正確輸入計算所需數(shù)據(jù)！, 48, 錯誤提示 Else ed = Val() Unload Me End IfEnd SubPrivate Sub Command2_Click() Unload MeEnd SubPrivate Sub Form_Load() = Str(ed)End SubOption ExplicitPrivate Sub Command1_Click() Dim H As Single, WH As Single, m As Single, R As Single, RF As Single, D As Single, am As Single, ah As Single, Wr As Boolean Dim xa As Single, ya As Single, xb As Single, xc As Single, xc_max As Single, xd As Single, yd As Single, ye As Single, xf As Single, xg As Single, a1 As Single, a2 As Single Dim h1 As Single, h2 As Single, Alf1 As Single, Alf2 As Single, W1 As Single, W2 As Single, F1 As Single, F2 As Single Dim a As Single, b As Single, c As Single, Dlt As Double, k As Single, kc As Single, xd_f As Single, yd_f As Single kc = 5 Wr = FalseIf = True Then If = 7 Then am = 3 ah = ElseIf = 8 Then am = ah = ElseIf = 9 Then am = 2 ah = Else Wr = True End If End If If = Or = Or = Or = Or = Or = Or Wr = True Then MsgBox 請正確輸入計算所需數(shù)據(jù)！, 48, 錯誤提示 Else H = Val() WH = Val() m = Val() R = Val() RF = Val() D = Val() xb = H * m For xa = 0 To xb ed Step ed ya = xa / m For xd = xa + ed To xb Step ed For yd = ya + ed To xd / m ed Step ed a2 = Atn((yd ya) / (xd xa)) ye = xd / m xc_max = xb + D If xa + (H ya) / Tan(a2) xc_max Then xc_max = xa + (H ya) / Tan(a2) End If For xc = xb To xc_max Step ed If xd xc Then a1 = Atn((H yd) / (xc xd)) If WH ya Or WH = ya Then W1 = ((ye yd) * (xb xd) + (xc xb) * (H yd)) * R / 2 W2 = (ye yd) * (xd xa) * R / 2 ElseIf WH ya And WH yd Or WH = yd Then xf = WH * m xg = xa + (WH ya) / Tan(a2) W1 = ((ye yd) * (xb xd) + (xc xb) * (H yd)) * R / 2 W2 = ((xg xf) * (WH ya) * RF + ((xg xf) * (ye WH) + (ye yd) * (xd xg)) * R) / 2 ElseIf WH yd And WH ye Or WH = ye Then xf = WH * m xg = xd + (WH yd) / Tan(a1)W1 = (((ye WH) * (xb xd) + (xc xb) * (H WH) + (xg xd) * (H WH)) * R +(xg xd) * (WH yd) * RF) / 2W2 = ((xd xf) * (ye WH) * R + ((xd xf) * (WH ya) + (WH yd) * (xd xa)) * RF) / 2 ElseIf WH ye Then xf = WH * m xg = xd + (WH yd) / Tan(a1) W1 = (((xg xf) + (xc xb)) * (H WH) * R + ((xg xf) * (WH ye) + (ye yd) * (xg xd)) * RF) / 2 W2 = (ye yd) * (xd xa) * RF / 2 End If If = True Then a = W1 * Sin(a1) * Cos(a1 a2) + W2 * Sin(a2) b = W1 * Cos(a1) * Cos(a1 a2) * W1 * Sin(a1) * Sin(a1 a2) * W2 * Cos(a2) * c = W1 * Cos(a1) * Sin(a1 a2) * * ElseIf = True Then h1 = (ya + yd + ye) / 3 If xc = xb Then h2 = (ye + yd + H) / 3 Else h2 = (ye + yd + H + H) / 4 End If If H 40 Or H = 40 Then Alf1 = 1 + h1 / H * (am 1) Alf2 = 1 + h2 / H * (am 1) Else If h1 * H Or h1 = * H Then Alf1 = 1 + (am 1) / 3 * h1 / ( * H) Else Alf1 = 1 + (am 1) / 3 + (h1 * H) / ( * H) * (am 1 (am 1) / 3) End If If h2 * H Or h2 = * H Then Alf2 = 1 + (am 1) / 3 * h2 / ( * H) ElseAlf2 = 1 + (am 1) / 3 + (h2 * H) / ( * H) * (am 1 (am 1) / 3) End If End If F1 = ah * * W1 * Alf1 F2 = ah * * W2 * Alf2 a = (W1 * Sin(a1) + F1 * Cos(a1)) * Cos(a1 a2) + W2 * Sin(a2) + F2 * Cos(a2) b = W1 * Cos(a1) * Cos(a1 a2) * (W1 * Sin(a1) + F1 * Cos(a1)) * Sin(a1 a2) * W2 * Cos(a2) * c = W1 * Cos(a1) * Sin(a1 a2) * * End If Dlt = b ^ 2 4 * a * c k = (b + Sqr(Dlt)) / (2 * a) If k kc Then kc = k xd_f = xd yd_f = yd End If End If Next xc Next yd Next xd Next xa = amp。 水位 amp。 Str(WH) amp。 m， amp。 坡度1: amp。 Str(m) amp。 時，滑坡面折點坐標(biāo)為 amp。 ( amp。 Str(xd_f) amp。 , amp。 Str(yd_f) amp。 ) amp。 ， amp。 安全系數(shù)為 amp。 Str(kc) amp。 。 amp。 Chr(13) amp。 Chr(10) End IfEnd SubPrivate Sub Command2_Click() = 8 = 62 = = = = = 7 = End Sub Private Sub Form_Load() = 8 = 62 = = = = = 7 = = TrueEnd SubPrivate Sub Form_Unload(Cancel As Integer) End SubPrivate Sub Form_Load() = 、參數(shù)等均遵照《水工建筑物荷載設(shè)計規(guī)范(DL+50771997)》、《水工建筑物抗震設(shè)計規(guī)范(DL50732000)》、《碾壓式土石壩設(shè)計規(guī)范(SL2742001)》的相關(guān)解釋說明； amp。 Chr(13) amp。 Chr(10) amp。 Chr(13) amp。 Chr(10) amp。 ，模型有所簡化； amp。 Chr(13) amp。 Chr(10) amp。 Chr(13) amp。 Chr(10) amp。 ，過高的計算精度會嚴(yán)重影響您的計算機的響應(yīng)速度，由此我們推薦您的計算機配置有雙(多)核心的中央處理器； amp。 Chr(13) amp。 Chr(10) amp。 Chr(13) amp。 Chr(10) amp。 ，本程序不負(fù)任何責(zé)任。 amp。 Chr(13) amp。 Chr(10) amp。 Chr(13) amp。 Chr(10) amp。 2015/5/7End Sub附錄2 GeoStudio計算結(jié)果1. 河床中間22斷面滲流圖21 正常蓄水位滲流計算圖注：縱坐標(biāo)壩高，橫坐標(biāo)壩長圖22 設(shè)計洪水位位滲流計算圖注：縱坐標(biāo)壩高，橫坐標(biāo)壩長圖23 校核洪水位位滲流計算圖注：縱坐標(biāo)壩高，橫坐標(biāo)壩長2. 最危險滑裂面穩(wěn)定計算