海浪波长以波浪力计算

Option Explicit Dim L1 As Single, L2 As Single, t As Single, d!, k!, kd!, thkd!, H!, D1! Dim CD As Single, CM As Single, l As Single, Ko As Single Dim Fhdmax As Single, Fhlmax As Single, Mhdmax As Single, Mhlmax!, Fhmax!, Mhmax! Dim θ As Single Const Pi = 3.141592653 Const G = 9.8 Const γ = 1025Private Sub Command1_Click() Dim r As Integer Do While True L1 = Val(InputBox(" 请输入波长L1：", " 求解设计波长：", "100")) t = Val(InputBox(" 请输入设计波周期T：", " 请输入 ", "6")) d = Val(InputBox(" 请输入设计水深d：", " 请输入 ", "20")) If L1 <= 0 Then r = MsgBox(" 请输入一个正数！", 5, " 输入错误 ") If r = 2 Then End End If Else Exit Do End If Loop k = 2 * Pi / L1 kd = k * d thkd = (Exp(kd) - Exp(-kd)) / (Exp(kd) + Exp(-kd)) L2 = G * (t ^ 2) * thkd / (2 * Pi) Do Until Abs(L2 - L1) < 0.001 L1 = L2 k = 2 * Pi / L1 kd = k * d thkd = (Exp(kd) - Exp(-kd)) / (Exp(kd) + Exp(-kd)) L2 = G * (t ^ 2) * thkd / (2 * Pi) Loop Print " 设计波长是： "; L2 Print " 波数： "; Format$(k, "0.0000") End Sub Private Sub Command2_Click() End End Sub Private Sub Command3_Click() H = Val(InputBox(" 请输入设计波高H：", "请输入 ", "3")) D1 = Val(InputBox(" 请输入桩柱直径D1：", " 请输入 ", "2")) l = Val(InputBox(" 请输入桩柱间距l：", "请输入 ", "15")) If d / L2 < 0.5 Then Print Print " 相对水深 d/L2:"; d / L2 Print " 采用线性波理论计算:" Else MsgBox " 重新选择计算理论" End If Print " 波陡： "; H / L2 Print " 相对柱径： "; D1 / L2 If D1 / L2 < 0.2 Then Print " 属于小直径桩柱" Else Print " 属于大直径桩柱" End If CD = Val(InputBox(" 请输入拖曳力系数：", "请输入 ", "1.0")) CM = Val(InputBox(" 请输入质量系数：", "请输入 ", "2.0")) Print Print " 选用拖曳力系数："; CD Print " 选用质量系数："...