海浪波長(zhǎng)以及波浪力計(jì)算

1、Op ti on Exp licitDim L1 As Sin gle, L2 As Sin gle, t As Sin gle, d!, k!, kd!, thkd!, H!, D1!Dim CD As Sin gle, CM As Sin gle, l As Sin gle, Ko As Sin gleDim Fhdmax As Sin gle, Fhlmax As Si ngle, Mhdmax As Sin gle, Mhlmax!, Fhmax!, Mhmax!Dim 0 As SgieConst Pi = 3.141592653Co nst G = 9.8Co nst Y = 10

2、25Private Sub Comma nd1_Click()Dim r As In tegerDo While True", "100")T :","請(qǐng)輸入", "6")d:","請(qǐng)輸入", "20")L1 = Val(InputBox("請(qǐng)輸入波長(zhǎng)L1 :","求解設(shè)計(jì)波長(zhǎng): t = Val(InputBox("請(qǐng)輸入設(shè)計(jì)波周期 d = Val(InputBox("請(qǐng)輸入設(shè)計(jì)水深I(lǐng)f L1 <= 0 The

3、 n",5,"輸入錯(cuò)誤”)r = MsgBox("請(qǐng)輸入一個(gè)正數(shù)！If r = 2 ThenEndEnd IfElseExit DoEnd IfLoopk = 2 * Pi / L1kd = k * dthkd = (Ex p(kd) - Exp(-kd) / (Ex p(kd) + Ex p(-kd)L2 = G * (t A 2) * thkd / (2 * Pi)Do Un til Abs(L2 - L1) < 0.001L1 = L2k = 2 * Pi / L1kd = k * dthkd = (Exp (kd) - Exp (-kd) / (Ex

4、 p(kd) + Exp (-kd)L2 = G * (t A 2) * thkd / (2 * Pi)LoopPrint"設(shè)計(jì)波長(zhǎng)是："L2Print "波數(shù)："Format$(k, "0.0000")End SubPrivate Sub Comma nd2_Click()EndEnd SubPrivate Sub Comma nd3_Click()H = Val(lnputBox("請(qǐng)輸入設(shè)計(jì)波高 H :","請(qǐng)輸入", "3")D1 = Val(InputBox(&qu

5、ot;請(qǐng)輸入樁柱直徑 D1 :","請(qǐng)輸入", "2") l = Val(InputBox("請(qǐng)輸入樁柱間距 I:","請(qǐng)輸入", "15")If d / L2 < 0.5 ThenPrintPrint "相對(duì)水深 d/L2:" d / L2Print "采用線性波理論計(jì)算:"ElseMsgBox "重新選擇計(jì)算理論”End IfPrint"波陡："H / L2Print"相對(duì)柱徑："D1

6、/ L2If D1 / L2 < 0.2 The nPrint "屬于小直徑樁柱”ElsePrint "屬于大直徑樁柱”End IfCD = Val(InputBox("請(qǐng)輸入拖曳力系數(shù)：","請(qǐng)輸入", "1.0")CM = Val(InputBox("請(qǐng)輸入質(zhì)量系數(shù)：","請(qǐng)輸入", "2.0")Print”;CD”；CMPrint"選用拖曳力系數(shù)：Print"選用質(zhì)量系數(shù)：Dim LD As Si ngleLD = l /

7、D1"LDH.Print"樁柱相對(duì)間距：Print"群樁系數(shù)Ko :'If LD > 4 The nKo = 1Print KoElseIf LD < 4 And LD > 3 The nKo = 1.25Print KoElseIf LD < 2 The nKo = 1.5Print KoEnd IfEnd SubPrivate Sub Comma nd4_Click()Dim K1 As Sin gle, K2 As Si ngle, K3 As Si ngle, K4 As Si ngleDim e As Single,0 o

8、 As SingleK1 = (2 * k * (d + H / 2) + sh(2 * k * (d + H / 2) / (8 * sh(2 * k * d)Fhdmax = CD * 丫 * G * D1 * (H 2) * K1 / 2PrintPrin t "K1 值："Format$(K1, "0.0000")Print "單樁柱最大水平拖曳力Fhdmax: " FhdmaxK2 = th(k * d)Fhlmax = CM * 丫 * G * Pi * (D1 人 2) * H * K2 / 8PrintPrint &q

9、uot;K2 值："Format$(K2, "0.0000")Print "單樁柱最大水平慣性力Fhlmax:" FhlmaxK3 = (2 * (k A 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)PrintPrint "K3 值：&q

10、uot;Format$(K3, "0.0000")Print "單樁柱最大水平拖曳力矩Mhdmax:" MhdmaxK4 = (k * d * sh(k * d) - ch(k * d) + 1) / ch(k * d)Mhlmax = CM * 丫 * G * (D1 人 2) * H * L2 * K4 / 16PrintPrint "K4 值："Format$(K4, "0.0000")Print "單樁柱最大水平慣性力矩Mhlmax:" MhlmaxIf Fhlmax >= 2 *

11、 Fhdmax The nFhmax = Fhlmax0 o = 90Elself Fhlmax < 2 * Fhdmax The nFhmax = Fhdmax * (1 + (Fhlmax / Fhdmax)人 2) / 4)0 o = arcsin(Fhlmax / (2 * Fhdmax)End IfPrintPrint "單樁柱最大水平波力Fhmax:" FhmaxIf Mhlmax >= 2 * Mhdmax The nMhmax = MhlmaxElself Mhlmax < 2 * Mhdmax The nMhmax = Mhdmax *

12、(1 + (Mhlmax / Mhdmax)人 2) / 4)End IfPrintPrint "單樁柱最大水平波力矩Mhmax:" MhmaxPrint"最大水平波力和最大水平波力矩的相位0o:" 0 oe = Mhmax / FhmaxPrintPrint "最大水平波力作用點(diǎn)離海底的距離e:" eEnd SubPublic Fun cti on sh( n) As Sin glesh = (Ex p(n) - Exp (-n) / 2End FunctionPublic Fun cti on ch( n) As Sin glec

13、h = (Exp(n) + Exp (-n) / 2End FunctionPublic Fu nction th( n) As Si ngleth = (Exp(n) - Exp (-n) / (Ex p(n) + Exp (-n)0)End FunctionPublic Fun cti on arcsin(n) As Sin gle arcsin = Atn(n / Sqr(-n * n + 1)End FunctionPublic Fu nction FH( FH = Fhdmax * Cos(End FunctionPublic Fu nction MH(MH = Mhdmax * C

14、os(End FunctionPrivate Sub Comma nd5_Click()Dim i As In teger0 ) As Single0 ) * Abs(Cos(0 ) As Single0 ) * Abs(Cos(0 ) + Fhlmax * Sin(0 ) + Mhlmax * Sin(0)Print #1, Tab(8);"相位角 0 ”； Spc(3);桌面不同相位水平波力.txt" For Out put As #1 "cos 0 " Spc(3); "cos 0 |cos0 I" Spc(3); "

15、sin 0 "0 " Spc(6); "FH"Spc(3); "Fhdmaxcos 0 |cos0 |" Spc(3); "FhlmaxsinFor i = 0 To 180 Ste p 150 = i * Pi / 1800 ), "0.0000");)* Ab s'(C0000"); 0 ),"0.0000");0 ) * Abs(Cos( 0 ), "0.00"); 0 ), "0.00");),"0.00&qu

16、ot;)Print #1, Tab(10); i; Tab(20); Format$(Cos(Print #1, Tab(30); Format$(Cos(Print #1, Tab(40); Format$(S in(Print #1, Tab(55); Format(Fhdmax * Cos(Print #1, Tab(70); Format(Fhlmax * Sin(Print #1, Tab(85); Format(FH(Next iClose #1End SubPrivate Sub Comma nd6_Click()Dim i As In teger桌面 不同相位水平波力矩.txt

17、" For Out put As#2Print #2, Tab(8);"相位角 0 " Spc(3); "cos 0 " Spc(3); "cos 0 |cos0 |" Spc(3); "sin 0 " Spc(3); "Fhdmaxcos 0 |cos0 |" Spc(3); "Mhlmaxsin 0 " Spc(6); "MH"For i = 0 To 180 Ste p 150 = i * Pi / 180Print #2, Tab(10)

18、; i; Tab(20); Format$(Cos(0 ), "0.0000");Print #2, Tab(30); Format$(Cos(0 ) * Abs(Cos( 0 ), "0.0000");Print #2, Tab(40); Format$ (Sin( 0 ), "0.0000");Print #2, Tab(55); Format(Mhdmax * Cos(0 ) * Abs(Cos(0 ), "0.00");Print #2, Tab(70); Format(Mhlmax * Sin(0 ),

19、"0.00");Print #2, Tab(85); Format(MH(0 ), "0.00")Next iClose #2End SubPrivate Sub Comma nd7_Click()Dim e As SingleDim 刀 H As SingleDim i As In tegerDim t As Sin gle, m As In tegere = l * 360 / L2PrintPrint "前后兩樁柱的波浪位相差e：t = 0For i = 0 To 180 - e刀 H = FH(i) + FH(i + e )If t < 刀 H Thent =刀Hm = iEnd IfNext i"m"；tPrint "發(fā)生最大水平合波力的相位:Print "前后兩樁柱的最大水平合波力為:End SubPrivate Sub Comma nd8_C