根据浪、流、和风的变化模型以及公式计算载荷
中提到的三个量
的变化模型，
即（wave ， ，wind ）。
2. 浪、流、和风的计算公式如下：
3. 根据上述浪、流、和风的变化模型，公式，船舶的物理参数等可以得到船舶载荷计算代码：
alfa=0.7; %海风舷角rad
beta=0.7; %海流入射角rad
khai=0.7; %海浪遭遇角，艏向零度顺时针gm2增加，流力完全推动船舶gm2=0或者2*pilmd=0.2; %海浪波长g = 9.81; %重力加速度m/s^2L = 1.70; %垂线间长m
Rar = 1.29; %空气密度约为1.29Kg/m^3
Rwt = 1040; %海水密度一般在1020～1070 kg/m^3之间Cxwd=0.6; Axwd=0.13; %纵向风阻系数和水线上面积
Cywd=0.8; Aywd=0.20; %横向风阻系数和水线上面积
Cnwd=0.5; Anwd=0.35; %艏摇方向风阻系数和水线上面积Cxcr=0.3; Axcr=0.13; %纵向流阻系数和水线下面积
Cycr=0.7; Aycr=0.31; %横向流阻系数和水线下面积
Cncr=0.5; Ancr=0.66; %艏摇方向流阻系数和水线下面积Wu0=(Axcr*Cxcr*Rwt*abs(cos(beta))*cos(beta)*cr.^2)/2 + (Axwd*Cxwd*Rar*abs(cos(alfa))*cos(alfa)*wd.^2)/2 - (L*Rwt*g*cos(khai)*((51*lmd^3)/800 - (3*lmd^2)/16 + lmd/10 - (1/20))*zt.^2)/2;Wv0=(Aycr*Cycr*Rwt*abs(sin(beta))*sin(beta)*cr.^2)/2 + (Aywd*Cywd*Rar*abs(sin(alfa))*sin(alfa)*wd.^2)/2 - (L*Rwt*g*sin(khai)*(- (221*lmd^3)/200 + (389*lmd^2)/100 + (683*lmd)/200 - (23/50))*zt.^2)/2;Wr0=(Ancr*Cncr*Rwt*abs(sin(beta))*sin(beta)*cr.^2)/2 + (Anwd*Cnwd*Rar*abs(sin(alfa))*sin(alfa)*wd.^2)/2 + (L*Rwt*g*sin(khai)*((21*lmd^3)/800 - (79*lmd^2)/400 + (17*lmd)/50 - (11/100))*zt.^2)/2;
4. 在一定载荷角情况下，通过上述代码就获得了浪、流、和风对船舶的载荷大小。
subplot(131)
plot(Wu0,'x-b')
xlabel("step");ylabel("force")
title("surge load")
subplot(132)
plot(Wv0,'d-c')
xlabel("step");ylabel("force")
title("sway load")
subplot(133)
plot(Wr0,'s-g')
title("yaw load")
xlabel("step");ylabel("moment")5. 上述代码实现了浪、流和风模型对船舶载荷大小的可视化。
欢迎交流指正。
