www.pudn.com > ship.zip > luoxuanjiang.asv, change:2016-08-31,size:918b


function F=luoxuanjiang(X) 
global L B Cp Cb rou V n Dp Cpp wp; 
 
 
%% 需要输入的变量量:u,v,r,δ 
u=X(1); 
v=X(2); 
r=X(3)*pi/180; 
rn=r*L/V;%无因次r 
xpn=-0.4;%无因次xp 
beta=-atan(v/u);%漂角β 
delta=X(4)*pi/180;%转化为弧度 
 
%% --推力系数的计算-- 
kp(1)=(6*Cb+0.9*L/B-5.7)*10^(-3); 
kp(2)=(14*Cb-0.5*L/B-0.3)*10^(-3); 
 
wp0=0.55*Cb-0.2; 
for i=1:2 
  wp(i)=wp0-kp(i)*(beta-xpn*rn); 
  ua(i)=(1-wp(i))*u; 
  JP(i)=ua(i)/(n*Dp); 
  kT(i)=-0.4315*JP(i)+0.6179; 
  T(i)=rou*n^2*Dp^4*kT(i); 
end 
 
%% --推力减额系数的计算-- 
tp0=0.5*Cp-0.18; 
BR=beta+0.95*L*(r/V); 
xB=-0.068; 
lcb=xB*100/L; 
A=(B/L)/(1.3*(1-Cb)-3.1*lcb); 
Kt=0.00023*(A*L/Dp)-0.028; 
tp=tp0-Kt*BR; 
%% --桨力和力矩计算-- 
if delta>=0 
    a=1; 
else 
    a=-1; 
end 
XP=(1-tp)*(T(1)+T(2)); 
YP=0; 
NP=a*(1-tp)*Cpp*(T(1)-T(2)); 
F=[XP YP NP]; 
%XPn=XP/(0.5*rou*V^2*L^2); 
%YPn=YP/(0.5*rou*V^2*L^2); 
%NPn=NP/(0.5*rou*V^2*L^3); 
%F=[XPn YPn NPn];