前n项立方和公式海浪波长以波浪力计算

2020年11月9日发(作者：苏应衡)
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 γ = 1025
Private Sub Command1_Click()
Dim r As Integer
Do While True
L1 = Val(InputBox(请输入波长L1：求解设计波长：
t = Val(InputBox(请输入设计波周期T：请输入
d = Val(InputBox(请输入设计水深d：请输入
If L1 <= 0 Then
r = MsgBox(请输入一个正数！输入错误
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 设计波长是：
Print 波数：
End Sub
Private Sub Command2_Click()
End
End Sub
Private Sub Command3_Click()
H = Val(InputBox(请输入设计波高H：请输入
D1 = Val(InputBox(请输入桩柱直径D1：请输入
l = Val(InputBox(请输入桩柱间距l：请输入
If d L2 < 0.5 Then
Print
Print 相对水深dL2:
Print 采用线性波理论计算:
Else
MsgBox 重新选择计算理论
End If
Print 波陡：
Print 相对柱径：
If D1 L2 < 0.2 Then
Print 属于小直径桩柱
Else
Print 属于大直径桩柱
End If
CD = Val(InputBox(请输入拖曳力系数：请输入
CM = Val(InputBox(请输入质量系数：请输入
Print
Print 选用拖曳力系数：
Print 选用质量系数：
Dim LD As Single
LD = l D1
Print 桩柱相对间距：
Print 群桩系数Ko：
If LD > 4 Then
Ko = 1
Print Ko
ElseIf LD < 4 And LD > 3 Then
Ko = 1.25
Print Ko
ElseIf LD < 2 Then
Ko = 1.5
Print Ko
End If
End Sub
Private Sub Command4_Click()

Dim K1 As Single, K2 As Single, K3 As Single, K4 As Single
Dim e As Single, θo As Single
K1 = (2 * k * (d + H 2) + sh(2 * k * (d + H 2))) (8 * sh(2 * k * d))
Fhdmax = CD * γ * G * D1 * (H ^ 2) * K1 2
Print
Print 值：
Print 单桩柱最大水平拖曳力Fhdmax：
K2 = th(k * d)
Fhlmax = CM * γ * G * Pi * (D1 ^ 2) * H * K2 8
Print
Print 值：
Print 单桩柱最大水平惯性力Fhlmax:
K3 = (2 * (k ^ 2) * (d + H 2) ^ 2 + 2 * k * (d + H 2) * sh(2 * k * (d + H 2)) - ch(2 * k * (d
+ H 2)) + 1) (32 * sh(2 * k * d))
Mhdmax = CD * γ * G * D1 * (H ^ 2) * L2 * K3 (2 * Pi)
Print
Print 值：
Print 单桩柱最大水平拖曳力矩Mhdmax:
K4 = (k * d * sh(k * d) - ch(k * d) + 1) ch(k * d)
Mhlmax = CM * γ * G * (D1 ^ 2) * H * L2 * K4 16
Print
Print 值：
Print 单桩柱最大水平惯性力矩Mhlmax:
If Fhlmax >= 2 * Fhdmax Then
Fhmax = Fhlmax
θo = 90
ElseIf Fhlmax < 2 * Fhdmax Then
Fhmax = Fhdmax * ((1 + (Fhlmax Fhdmax) ^ 2) 4)
θo = arcsin(Fhlmax (2 * Fhdmax))
End If
Print
Print 单桩柱最大水平波力Fhmax:
If Mhlmax >= 2 * Mhdmax Then
Mhmax = Mhlmax
ElseIf Mhlmax < 2 * Mhdmax Then
Mhmax = Mhdmax * ((1 + (Mhlmax Mhdmax) ^ 2) 4)
End If
Print
Print 单桩柱最大水平波力矩Mhmax:
Print 最大水平波力和最大水平波力矩的相位θo:θo
e = Mhmax Fhmax
Print
Print 最大水平波力作用点离海底的距离e:
End Sub
Public Function sh(n) As Single
sh = (Exp(n) - Exp(-n)) 2
End Function
Public Function ch(n) As Single
ch = (Exp(n) + Exp(-n)) 2
End Function
Public Function th(n) As Single
th = (Exp(n) - Exp(-n)) (Exp(n) + Exp(-n))
End Function
Public Function arcsin(n) As Single
arcsin = Atn(n Sqr(-n * n + 1))
End Function
Public Function FH(θ) As Single
FH = Fhdmax * Cos(θ) * Abs(Cos(θ)) + Fhlmax * Sin(θ)
End Function
Public Function MH(θ) As Single
MH = Mhdmax * Cos(θ) * Abs(Cos(θ)) + Mhlmax * Sin(θ)
End Function
Private Sub Command5_Click()
Dim i As Integer
Open 桌面不同相位水平波力.txt
Print #1, Tab(8); 相位角θSpc(3); θSpc(3); θ|cosθ|Spc(3); θ
Spc(3); θ|cosθ|θ
For i = 0 To 180 Step 15
θ = i * Pi 180
Print #1, Tab(10); i; Tab(20); Format$$(Cos(θ), );
Print #1, Tab(30); Format$$(Cos(θ) * Abs(Cos(θ)),
Print #1, Tab(40); Format$$(Sin(θ),
Print #1, Tab(55); Format(Fhdmax * Cos(θ) * Abs(Cos(θ)),
Print #1, Tab(70); Format(Fhlmax * Sin(θ),
Print #1, Tab(85); Format(FH(θ),
Next i
Close #1
End Sub
Private Sub Command6_Click()
Dim i As Integer
Open 桌面不同相位水平波力矩.txt
#2
Print #2, Tab(8); 相位角θSpc(3); θSpc(3); θ|cosθ|Spc(3); θ
Spc(3); θ|cosθ|θ
For i = 0 To 180 Step 15
θ = i * Pi 180
Print #2, Tab(10); i; Tab(20); Format$$(Cos(θ),
Print #2, Tab(30); Format$$(Cos(θ) * Abs(Cos(θ)),
Print #2, Tab(40); Format$$(Sin(θ),
Print #2, Tab(55); Format(Mhdmax * Cos(θ) * Abs(Cos(θ)),
Print #2, Tab(70); Format(Mhlmax * Sin(θ),
Print #2, Tab(85); Format(MH(θ),
Next i
Close #2
End Sub
Private Sub Command7_Click()
Dim ψ As Single
Dim ∑H As Single
Dim i As Integer
Dim t As Single, m As Integer
ψ = l * 360 L2
Print
Print 前后两桩柱的波浪位相差ψ：ψ
t = 0
For i = 0 To 180 - ψ
∑H = FH(i) + FH(i + ψ)
If t < ∑H Then
t = ∑H
m = i
End If
Next i
Print 发生最大水平合波力的相位 ：
Print 前后两桩柱的最大水平合波力为：
End Sub
Private Sub Command8_Click()
Cls
End Sub
Private Sub Command9_Click()
Dim ψ As Single
Dim ∑M As Single
Dim i As Integer
Dim t As Single, m As Integer
ψ = l * 360 L2
t = 0
For i = 0 To 180 - ψ
∑M = MH(i) + MH(i + ψ)
If t < ∑M Then
t = ∑M
m = i
End If
Next i
Print
Print 发生最大水平合波力矩的相位 ：
Print 前后两桩柱的最大水平合波力矩为：
End Sub






海洋环境作业
计算结果：