根据浪、流、和风的变化模型以及公式计算载荷

1. 浪、流、和风的变化模型就是https://blog.csdn.net/ZLH_ZLH/article/details/134387003?spm=1001.2014.3001.5502https://blog.csdn.net/ZLH_ZLH/article/details/134387003?spm=1001.2014.3001.5502中提到的三个量的变化模型，

即（wave height，current velocity，wind velocity）。

2. 浪、流、和风的计算公式如下：

3. 根据上述浪、流、和风的变化模型，公式，船舶的物理参数等可以得到船舶载荷计算代码：

alfa=0.7;  %海风舷角rad
beta=0.7;  %海流入射角rad
khai=0.7;  %海浪遭遇角，艏向零度顺时针gm2增加，流力完全推动船舶gm2=0或者2*pi
 lmd=0.2;  %海浪波长

  g = 9.81;       %重力加速度m/s^2
  L = 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. 上述代码实现了浪、流和风模型对船舶载荷大小的可视化。

欢迎交流指正。