海浪波长以及波浪力计算整理版

2021-06-14 14:37:30本页面

海浪波长以及波浪力计算整理版


【正文】

. OptionExplicit DimL1AsSingle,L2AsSingle,tAsSingle,d!,k!,kd!,thkd!,H!,D1! DimCDAsSingle,CMAsSingle,lAsSingle,KoAsSingle DimFhdmaxAsSingle,FhlmaxAsSingle,MhdmaxAsSingle,Mhlmax!,Fhmax!,Mhmax!DimθAsSingle ConstPi=3.141592653 ConstG=9.8 Constγ=1025 PrivateSubCommand1Click() DimrAsInteger DoWhileTrue L1=Val(InputBox("。

请输入波长L1:","求解设计波长:","100")) t=Val(InputBox("请输入设计波周期T:","请输入","6")) d=Val(InputBox("请输入设计水深d:","请输入","20")) IfL1<=0Then r=MsgBox("请输入一个正数!",5,"输入错误") Ifr=2Then End EndIf Else ExitDo EndIf Loop k=2*Pi/L1 kd=k*d thkd=(Exp(kd)Exp(kd))/(Exp(kd)+Exp(kd)) L2=G*(t^2)*thkd/(2*Pi) DoUntilAbs(L2L1)<。

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") EndSub PrivateSubCommand2Click() End EndSub PrivateSubCommand3Click() H=Val(InputBox("请输入设计波高H:","

请输入","3")) D1=Val(InputBox("请输入桩柱直径D1:","请输入","2")) l=Val(InputBox("请输入桩柱间距l:","请输入","15")) . . Ifd/L2<0.5Then Print Print"相对水深d/L2:"d/L2 Print"采用线性波理论计算:" Else MsgBox"重新选择计算理论"。

EndIf Print"波陡:"H/L2 Print"相对柱径:"D1/L2 IfD1/L2<0.2Then Print"属于小直径桩柱" Else Print"属于大直径桩柱" EndIf CD=Val(InputBox("请输入拖曳力系数:","请输入","1.0"))CM=Val(InputBox("请输入质量系数:","请输入","2.0")) Print Print"。

选用拖曳力系数:"CD Print"选用质量系数:"CM DimLDAsSingle LD=l/D1 Print"桩柱相对间距:"LD Print"群桩系数Ko:" IfLD>4Then Ko=1 PrintKo ElseIfLD<4AndLD>3Then Ko=1.25 PrintKo ElseIfLD<2Then Ko=1.5 PrintKo EndIf EndSub PrivateSubCommand4Click() DimK1AsSingle,K2AsSingle。

K3AsSingle,K4AsSingleDimeAsSingle,θoAsSingle 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"K1值:"Format$(K1,"0.0000") Print"单桩柱最大水平拖曳力Fhdmax:"Fhdmax K2=th(k*d) . . Fhlmax=CM*γ*G*Pi*(D1^2)*H*K2/8 Print Print"K2值:"Format$(K2。

"0.0000") Print"单桩柱最大水平惯性力Fhlmax:"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"K3值:"Format$(K3,"0.0000") Print"单桩柱最大水平拖曳力矩Mhdmax:"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"。

K4值:"Format$(K4,"0.0000") Print"单桩柱最大水平惯性力矩Mhlmax:"Mhlmax IfFhlmax>=2*FhdmaxThen Fhmax=Fhlmax θo=90 ElseIfFhlmax<2*FhdmaxThen Fhmax=Fhdmax*((1+(Fhlmax/Fhdmax)^2)/4) θo=arcsin(Fhlmax/(2*Fhdmax)) EndIf Print Print"单桩柱最大水平波力Fhmax:"Fhmax IfMhlmax>=2*MhdmaxThen Mhmax=Mhlmax ElseIfMhlmax<。

2*MhdmaxThen Mhmax=Mhdmax*((1+(Mhlmax/Mhdmax)^2)/4) EndIf Print Print"单桩柱最大水平波力矩Mhmax:"Mhmax Print"最大水平波力和最大水平波力矩的相位θo:"θo e=Mhmax/Fhmax Print Print"最大水平波力作用点离海底的距离e:"e EndSub PublicFunctionsh(n)AsSingle sh=(Exp(n)Exp(n))/2 EndFunction PublicFunctionch(n)AsSingle ch=(Exp(n)+Exp(n))/2 EndFunction PublicFunctionth(n)AsSingle th=(Exp(n)Exp(n))/(Exp(n)+Exp(n)) 。

. EndFunction PublicFunctionarcsin(n)AsSingle arcsin=Atn(n/Sqr(n*n+1)) EndFunction PublicFunctionFH(θ)AsSingle FH=Fhdmax*Cos(θ)*Abs(Cos(θ))+Fhlmax*Sin(θ) EndFunction PublicFunctionMH(θ)AsSingle MH=Mhdmax*Cos(θ)*Abs(Cos(θ))+Mhlmax*Sin(θ) EndFunction PrivateSubCommand5Click() DimiAsInteger Open"。

C:\DocumentsandSettings\AllUsers\桌面\不同相位水平波力.txt"ForOutputAs#1Print#1,Tab(8);"相位角θ"Spc(3);"cosθ"Spc(3);"cosθ|cosθ|"Spc(3);"sinθ" Spc(3);"Fhdmaxcosθ|cosθ|"Spc(3);"Fhlmaxsinθ"Spc(6);"FH" Fori=0To180Step15 θ=i*Pi/180 Print#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(θ),"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") Nexti Close#1 EndSub PrivateSubCommand6Click() DimiAsInteger Open"C:\DocumentsandSettings\AllUsers\桌面\不同相位水平波力矩.txt"ForOutputAs#2 Print#2,Tab(8);"相位角θ"Spc(3);"cosθ"Spc(3);"cosθ|cosθ|"Spc(3);"sinθ"Spc(3);"Fhdmaxcosθ|cosθ|"。

Spc(3);"Mhlmaxsinθ"Spc(6);"MH" Fori=0To180Step15 θ=i*Pi/180 Print#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.0000"); Print#2,Tab(55);Format(Mhdmax*Cos(θ)*Abs(Cos(θ))。

家长会发言稿相关推荐  
三九文库 www.999doc.com
备案图标苏ICP备2020069977号