莫里森方程求柱体的动力响应

1、Monrison理论公式F=Fcos0cos0HHDmax+FHImaxsin0(1)其中FHDmaxCrDH2K'=D12(2)FHImaxCr兀D2HK'M28(3)HH2k(d+H)+sh2k(d+H)(4)K二K'二21 i8sh2kdshkdK=K'=(5)2 2chkdr=pg(6)液k=w2/g二1(7)0gTr式中FH是总水平波力，D是柱体的直径，H是波咼，d是水深，T是波的周期，p液是液体的密度，CD是拖拽系数，CM是质量系数Matlab程序functionFH=monrison(DD,H,d,T,p,Cd,Cm)%用莫里森公式计算柱体动力响应

2、%DD是柱体的直径，H是波高，d是水深，T波的周期，p是液体的密度%Cd是拖拽系数，CM是质量系数%使用语法%FH=monrison(DD,H,d,T,p,Cd,Cm)pp=pi;k=1/9.8/T/T;r=p*9.8;kk1=(2*k*(d+H/2)+sinh(2*k*(d+H/2)/(8*sinh(2*k*d);kk2=sinh(k*d)/cosh(k*d);FHD=Cd*r*DD*H*H*kk1/2;FHI=Cm*r*pp*DD*DD*H*kk2/8;q=0:pp/20:2*pp;FH=zeros(1,size(q,1);fori=1:41FH(1,i)=FHD*cos(q(1,i)*a

3、bs(cos(q(1,i)+FHI*sin(q(1,i);endplot(q,FH)gridonend调用函数计算>>FH=monrison(0.3,3,20,8,1,1,2)q=Columns1through1000.15710.31420.47120.62830.78540.94251.09961.25661.4137Columns11through201.57081.72791.88502.04202.19912.35622.51332.67042.82742.9845Columns21through303.14163.29873.45583.61283.76993.9270

4、4.08414.24124.39824.5553Columns31through404.71244.86955.02655.18365.34075.49785.65495.81195.96906.1261Column416.2832FH=Columns1through103.55453.47793.23562.85202.36541.82411.28170.79160.40240.1524Columns11through200.0662-0.0216-0.2764-0.6736-1.1745-1.7304-2.2875-2.7919-3.1946-3.4572Columns21through30-3.5545-3.4779-3.2356-2.8520-2.3654-1.8241-1.2817-0.7916-0.4024-0.1524Columns31through40-0.0662