DragScreen.rar_OCEAN engineering_fish net_matlab_海洋水动力_网箱

clear all;clc;close all; U=1.2;%来流速度 rho=1025;%来流密度 a1=15.75;%中间网长16.555 b1=13.65;%中间网高13.741 a2=17;%上/下网边长 A1=a1*b1/2;%中间每片面积 A2=1.5*sqrt(3)*a2^2;%上/下面积 Sn=[0.111,0.0396,0.2755];%密实度：内网0.111，外网0.0396，减速网0.2755 theta=[0,pi/6,pi/3,pi/2];%来流方向与x的夹角 alpha=[pi/6,pi/6,pi/2,pi/2,pi/6,pi/6,pi/6,pi/6,pi/2,pi/2,pi/6,pi/6;...%0度来流与平面夹角 pi/3,pi/3,pi/3,pi/3, 0, 0,pi/3,pi/3,pi/3,pi/3, 0, 0;...%30度来流与平面夹角 pi/2,pi/2,pi/6,pi/6,pi/6,pi/6,pi/2,pi/2,pi/6,pi/6,pi/6,pi/6;...%60度来流与平面夹角 pi/3,pi/3, 0, 0,pi/3,pi/3,pi/3,pi/3, 0, 0,pi/3,pi/3]; %90度来流与平面夹角 % alpha=[-5*pi/6,-5*pi/6,-pi/2,-pi/2,-pi/6,-pi/6,pi/6,pi/6,pi/2,pi/2,5*pi/6,5*pi/6;...%0度来流与平面夹角 % -2*pi/3,-2*pi/3,-pi/3,-pi/3, 0, 0,pi/3,pi/3,2*pi/3,2*pi/3, pi, pi;...%30度来流与平面夹角 % -pi/2,-pi/2,-pi/6,-pi/6,pi/6,pi/6,pi/2,pi/2,5*pi/6,5*pi/6,-5*pi/6,-5*pi/6;...%60度来流与平面夹角 % -pi/3,-pi/3, 0, 0,pi/3,pi/3,2*pi/3,2*pi/3, pi, pi,-2*pi/3,-2*pi/3]; %90度来流与平面夹角 for i=1:3 for j=1:4 for k=1:12 Cd1 (i,j,k)=0.04+(-0.04+Sn(i)-1.24*Sn(i)^2+1.37*Sn(i)^3)*cos(alpha(j,k)); Cl1(i,j,k)=(0.57*Sn(i)-354*Sn(i)^2+10.1*Sn(i)^3)*sin(2*alpha(j,k)); end end end % 1 2 3 4 5 6 7 8 9 10 11 12 beta(:,1,:)=[0.8 1 1 1 1 0.8 0.8 1 1 1 1 0.8;... 0.8 1 1 1 1 0.8 0.8 1 1 1 1 0.8;... 1 0 0 0 0 0.8 0.8 0 0 0 0 1];%0度速度衰减系数 beta(:,2,:)=[0.8 1 1 0.8 0.8 1 1 1 1 1 1 0.8;... 0.8 1 1 0.8 0.8 1 1 1 1 1 1 0.8;... 1 0 0 0 0 1 1 0 0 0 0 1];%30度速度衰减系数 beta(:,3,:)=[1 1 0.8 1 1 1 1 1 1 1 1 0.8;... 1 1 0.8 1 1 1 1 1 1 1 1 0.8;... 1 0 0 0 0 1 1 0 0 0 0 1];%60度速度衰减系数 beta(:,4,:)=[0.8 1 1 1 1 0.8 0.8 1 1 1 1 0.8;... 0.8 1 1 1 1 0.8 0.8 1 1 1 1 0.8;... 0.8 0 0 0 0 0.8 1 0 0 0 0 1];%90度速度衰减系数 Ue=U*beta;%有效速度 Fdnet1=0.5*rho*A1*Ue.*Ue.*Cd1;%中间网拖曳力 Flnet1=0.5*rho*A1*Ue.*Ue.*Cl1;%中间网升力 for j=1:4 Lift(3*(j-1)+1:3*(j-1)+3,:)=sum(sum(Flnet1(:,j,:),2),3); Drag(3*(j-1)+1:3*(j-1)+3,:)=sum(sum(Fdnet1(:,j,:),2),3); end Cd2=0.04; Fdnet2=0.5*rho*A2*U*U*Cd2;%上/下网拖曳力 for j=1:4 fdnet1(j,:)=sum(Fdnet1(:,j,:),1); Fdnet(j)=sum(fdnet1(j,:))+4*Fdnet2;%网衣拖曳力 end % Fdstr=[181982.4,177110.8,181982.4,177110.8];%钢结构拖曳力 Fdstr=[175213,180591,175213,180591];%钢结构拖曳力from AQWA Fd=Fdnet+Fdstr; %竖直网等效系数------------------------------------------------------------- w1=0.035;w2=0.15;w3=0.02; m1=15.75/w1;m2=13.75/w2;m3=17/2/w3; n1=13.65/w1;n2=13.65/w2;n3=13.65/w3; S1h=15.75*0.002*n1;S2h=13.75*0.003*n2;S3h=17/2*0.003*n3; S1v=13.65*0.002*m1;S2v=13.65*0.003*m2;S3v=13.65*0.003*m3; Sh=[S1h*0.5,S2h*0.5,S3h]; Sv=[S1v*0.5,S2v*0.5,S3v]; for j=1:4 pha(j,:)=[pi/3,pi/3,0,0,2*pi/3,2*pi/3,pi/3,pi/3,0,0,2*pi/3,2*pi/3]-theta(j)*ones(1,12); end for i=1:3 for j=1:4 for k=1:12 Se(i,j,k)=Sh(i).*abs((sin(pha(j,k)))^3)+Sv(i); fo(i,j,k)=0.5*rho*U*U*Se(i,j,k); Cde(i,j,k)=Fdnet1(i,j,k)/fo(i,j,k); end end end for i=1:3 for k=1:12 Cd00(i,k)=Cde(i,1,k);%0度等效系数（内网、外网、减速） Cd30(i,k)=Cde(i,2,k);%30度等效系数 Cd60(i,k)=Cde(i,3,k);%60度等效系数 Cd90(i,k)=Cde(i,4,k);%90度等效系数 end end %上下网等效系数----------------------------------------------------------- garma1=0:pi/m1/3:pi;%内网径向 garma2=0:pi/m2/3:pi;%外网径向 garma3=[pi/3,0,2*pi/3];%内外网周向 l=17/10:17/10:17;%周向边长 for j=1:4 phai(j,:)=garma1-theta(1)*ones(1,length(garma1)); phao(j,:)=garma2-theta(1)*ones(1,length(garma2)); phac(j,:)=garma3-theta(j)*ones(1,length(garma3)); end Sr1o=17*0.002.*abs(sin(phai).*sin(phai).*sin(phai)); Sr2o=17*0.003.*abs(sin(phao).*sin(phao).*sin(phao)); Sr1=sum(Sr1o,2);%内网径向等效面积 Sr2=sum(Sr2o,2);%外网径向等效面积 for j=1:4 for m=1:3 for n=1:length(l) Sc1o(j,m,n)=abs((sin(phac(j,m)))^3)*l(n)*0.002; Sc2o(j,m,n)=abs((sin(phac(j,m)))^3)*l(n)*0.003; end end end Sc1=sum(sum(Sc1o,3),2);%内网周向等效面积 Sc2=sum(sum(Sc2o,3),2);%外网周向等效面积 Se1=2*(Sr1+Sc1);%内网等效面积 Se2=2*(Sr2+Sc2);%外网等效面积 Cde1=Fdnet2*ones(4,1)./(0.5*rho*U*U*Se1);%上下内网等效系数（0 30 60 90） Cde2=Fdnet2*ones(4,1)./(0.5*rho*U*U*Se2);%上下外网等效系数 % F90=squeeze(Fdnet1(:,4,:)); % fo90=squeeze(fo(:,4,:)); % i=3;j=4;k=1; % Fdnet1(i,j,k)/fo(i,j,k)