• 大小: 11KB
    文件类型: .zip
    金币: 2
    下载: 0 次
    发布日期: 2024-02-13
  • 语言: 其他
  • 标签:

资源简介

船舶建立数学模型可以分为两种方法[ ]:第一种方法是从基本运动方程出发,以日本学派为代表的船舶运动分离型数学模型和以欧美学派为代表的整体型船舶运动数学模型。

资源截图

代码片段和文件信息

% %RungeKuttta4.m

% function dy=MMG1(ty)
% u=y(1);v=y(2);r=y(3);xg=y(4);yg=y(5);course=y(6);xB=0;% xB 是漂心
% % Cr=(0.0067*u^4-0.166*u^3+1.553*u^2-6.3699*u+10.02)*10^(-3);


% %假设
% % u=1;v=0.01;r=0.0001;
% xC=0;% xf 是浮心

% u0=13.5*1852/3600;%=0.6945
% p=1.025;M=14278.12*p;Lw=124;B=20.8;dm=8;Cb=0.681;Cp=0.692;panbi=0.67;%盘面比
% Vpai=14278.12;%排水体积
% Dp=4.6;P=3.66;L=126;n=154/60;
% zhanbi=1.72;%展铉比
% Ar=18.8;Hr=6.1;%舵高
% mper=M/(0.5*p*L^2*dm);%mxper=0.0133;myper=0.2022;Jzzper=0.046;Izzper=0.0113;
% %   =0.2248
% V=(u^2+v^2)^0.5;
% Izzper=(1+Cb^4.5)*p*M/24*(L^2+B^2)/(0.5*p*L^4*dm);  % p145 杨论文
% %     =0.0116
% % mxper=mper/100*(0.398+11.97*Cb*(1+3.73*dm/B)-2.89*Cb*L/B*(1+1.113*dm/B)+0.175*Cb*(L/B)^2*(1+0.541*dm/B)-1.107*L/B*dm/B); 
% %    =0.0133
% mxper=0.0112;

% myper=1.129*mper*(0.882-0.54*Cb*(1-1.6*dm/B)-0.156*(1-0.673*Cb)*L/B+0.826*dm/B*L/B*(1-0.678*dm/B)-0.638*Cb*dm/B*L/B*(1-0.669*dm/B));
% %    =0.2283
% Jzzper=0.01*(33-76.85*Cb*(1-0.784*Cb)+3.43*L/B*(1-0.63*Cb))*mper/2.6;%0.0177
% % Jzzper=mper*(Jzz/L)^2;%     =0.0094

% %计算Xuu
%   %计算Ct
%     %计算水湿面积
% %       k=5.5;
%      k=2.77;
%       
%     S=k*(Vpai*Lw)^0.5;%桑地公式p149杨
%     %计算Cf p151杨
%       vv=1.18831*10^(-6);
%       Rn=u*L/vv;%=7*10^8
%     Cf1=0.066/(log10(Rn)-2.03)^2;% Rn雷诺数
%     Cf2=0.455/(log10(Rn))^2.58;
%     Cf3=0.075/(log10(Rn)-2.03)^2;%=0.0016
%     Cf4=0.455/(log10(Rn))^2.58-1700/Rn;
%     Car=0.000362;
%     Cr=(0.0067*V^4-0.166*V^3+1.553*V^2-6.3699*V+10.02)*10^(-3)*0.85/0.95;%=6*10^(-4)
%  %    Cr=(0.0067*V^4-0.166*u^3+1.553*u^2-6.3699*u+10.02)*10^(-3)
%   Ct4=Cf3+Cr+Car;%=0.0032
%   Ct=Ct4*1.001;
% Xuu=-S/(L*dm)*Ct;%=-0.0089
% Xvv=0.4*B/L-0.006*L/dm;%=-0.0285
% Xrr=0.0005*L/dm;%=0.0079
% Xvr=(1.11*Cb-0.07)*myper;%=0.1387
% XH=Xuu*(u/u0)^2+Xvv*(v/u0)^2+Xvr*(v/u0)*(r*L/u0)+Xrr*(r*L/u0)^2;% 0--0.0123


% %计算Xp=(1-tp)*T
%    %计算 tp 旋转中的推力减额系数
%        tp0=0.5*Cp-0.12;
%        tp0per=0.04+(tp0-0.04)*u/u0;%=0.08--0.1
%        lR=-0.95;
%        pjiao=atan(-v/u);
%        pjiaoR=pjiao-lR*r*L/u0;%pjiaoP=atan(-(v+lp*r)/u) p109 杨
% % pjiaoR=atan(-(v+lR*r)/u);

%  % 求 tp  p111 杨       
%        lcb=xC/L*100;
%        rA=(B/dm)*(1.3*(1-Cb)-3.1*lcb);
%        kt=0.00023*(rA*L/Dp)-0.028;%=-0.0212
%        f=kt*pjiaoR;
%     tp=tp0per-f; %=-0.0173
%    %计算T=p*n^2*Dp^4*Kt
%      % Kt Kq的计算
% %       pjiaoP=pjiao+0.5*L*r/V;
%      lp=-0.5*L;
% %       pjiaoP=atan(-(v+lp*r)/u);%p109 杨
%         pjiaoP=pjiao+0.5*L/u0*r;%p104 人民交通出版社 操纵与耐波性  吴秀恒
%       wp0=0.5*Cb-0.05 ; % =0.2905 p105 杨
%       wp=wp0*exp(-4*(pjiaoP)^2);%p104 人民交通出版社 操纵与耐波性  吴秀恒 =0.295
%       J=(1-wp)*u/(n*Dp);%= 0.3818p148 杨 论文
%      Kt=0.0821*J^3-0.2308*J^2-0.2768*J+0.3499;%=0.2918
%      Kq=0.0085*J^3-0.027*J^2-0.0252*J+0.0416;%=0.0361
% XP=(1-tp)*p*n^2*Dp^4*Kt/(p*L*dm*u0^2);%=0.0376--0.016
% %计算 Xr=(1-tR)*Fn*sin(rudder)
%         %计算 Fn=-0.5*p*Ar*fa*Vr^2*sin(aR);
%           zx=-22.2*(Cb*B/L)^2+0.02*(Cb*B/

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       11152  2017-10-17 16:00  MMG模型\MMG1.m
     文件        2533  2012-10-09 20:31  MMG模型\new_MMG.m
     文件       26888  2012-10-09 09:16  MMG模型\z_type_rh.mdl
     目录           0  2018-01-23 19:33  MMG模型\

评论

共有 条评论

相关资源