海浪波长以及波浪力计算

1、.option explicitdim 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 singledim fhdmax as single, fhlmax as single, mhdmax as single, mhlmax!, fhmax!, mhmax! dim as singleconst pi = 3.141592653const g = 9.8const = 1025private sub c

2、ommand1_click()dim r as integerdo while truel1 = val(inputbox(请输入波长 l1：, 求解设计波长：, 100)t = val(inputbox(请输入设计波周期 t：, 请输入, 6)d = val(inputbox(请输入设计水深 d：, 请输入, 20)if l1 = 0 thenr = msgbox(请输入一个正数！, 5, 输入错误)if r = 2 thenendend ifelseexit doend ifloopk = 2 * pi / l1kd = k * dthkd = (exp(kd) - exp(-kd) /

3、(exp(kd) + exp(-kd)l2 = g * (t 2) * thkd / (2 * pi)do until abs(l2 - l1) 0.001l1 = l2k = 2 * pi / l1kd = k * dthkd = (exp(kd) - exp(-kd) / (exp(kd) + exp(-kd)l2 = g * (t 2) * thkd / (2 * pi)loopprint 设计波长是：; l2print 波数：; format$(k, 0.0000)end subprivate sub command2_click()endend subprivate sub comm

4、and3_click()h = val(inputbox(请输入设计波高 h：, 请输入, 3)d1 = val(inputbox(请输入桩柱直径 d1：, 请输入, 2)l = val(inputbox(请输入桩柱间距 l：, 请输入, 15).if d / l2 0.5 thenprintprint 相对水深 d/l2:; d / l2print 采用线性波理论计算:elsemsgbox 重新选择计算理论end ifprint 波陡：; h / l2print 相对柱径：; d1 / l2if d1 / l2 4 thenko = 1print koelseif ld 3 thenko =

5、 1.25print koelseif ld = 2 * fhdmax thenfhmax = fhlmaxo = 90e lseif fhlmax = 2 * mhdmax thenmhmax = mhlmaxelseif mhlmax 2 * mhdmax thenmhmax = mhdmax * (1 + (mhlmax / mhdmax) 2) / 4)end ifprintprint 单桩柱最大水平波力矩 mhmax:; mhmaxprint 最大水平波力和最大水平波力矩的相位o:; oe = mhmax / fhmaxprintprint 最大水平波力作用点离海底的距离 e:; e

6、end subpublic function sh(n) as singlesh = (exp(n) - exp(-n) / 2end functionpublic function ch(n) as singlech = (exp(n) + exp(-n) / 2end functionpublic function th(n) as singleth = (exp(n) - exp(-n) / (exp(n) + exp(-n).end functionpublic function arcsin(n) as singlearcsin = atn(n / sqr(-n * n + 1)en

7、d functionpublic function fh() as singlefh = fhdmax * cos() * abs(cos() + fhlmax * sin()end functionpublic function mh() as singlemh = mhdmax * cos() * abs(cos() + mhlmax * sin()end functionprivate sub command5_click()dim i as integeropen c:documents and settingsall users桌面不同相位水平波力.txt for output as

8、 #1 print #1, tab(8); 相位角; spc(3); cos; spc(3); cos|cos|; spc(3); sin;spc(3); fhdmaxcos|cos|; spc(3); fhlmaxsin; spc(6); fhfor i = 0 to 180 step 15 = i * pi / 180print #1, tab(10); i; tab(20); format$(cos(), 0.0000);print #1, tab(30); format$(cos() * abs(cos(), 0.0000);print #1, tab(40); format$(sin

9、(), 0.0000);print #1, tab(55); format(fhdmax * cos() * abs(cos(), 0.00);print #1, tab(70); format(fhlmax * sin(), 0.00);print #1, tab(85); format(fh(), 0.00)next iclose #1end subprivate sub command6_click()dim i as integero pen c:documents and settingsall users 桌面不同相位水平波力矩.txt for output as #2p rint

10、 #2, tab(8); 相位角; spc(3); cos; spc(3); cos|cos|; spc(3); sin; spc(3); fhdmaxcos|cos|; spc(3); mhlmaxsin; spc(6); mhfor i = 0 to 180 step 15 = i * pi / 180print #2, tab(10); i; tab(20); format$(cos(), 0.0000);print #2, tab(30); format$(cos() * abs(cos(), 0.0000);print #2, tab(40); format$(sin(), 0.00

11、00);print #2, tab(55); format(mhdmax * cos() * abs(cos(), 0.00);print #2, tab(70); format(mhlmax * sin(), 0.00);print #2, tab(85); format(mh(), 0.00)next iclose #2end subprivate sub command7_click().dim as singledim h as singledim i as integerdim t as single, m as integer = l * 360 / l2printprint 前后两桩柱的波浪位相差：; t = 0for i = 0 to 180 - h = fh(i) + fh(i + )if t h thent = hm = iend ifnext iprint 发生最大水平合波力的相位 ：; m print 前后两桩柱的最大水平合波力为：; tend subprivate sub command8_click()clsend subprivate sub command9_click()dim as singledim m as singledim i as integerdim t as