pavailable3.m

用MATLAB编写关于飞机设计的matlab程序

function  [Pky,T]=Pavailable3(Vin,vf)
%   该函数求出给定电压 Vin 和前飞速度 vf 时的螺旋桨稳定转速nps,
%   利用nps计算此时的可用功率和拉力
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     RHO=1.225;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     PROpara;   %调用螺旋桨的特性参数文件，含有进矩比Jar，功率系数Cp，效率Eta,直径Dp
     MOTORpara; %调用马达的特性参数文件，含有Kv,Kt,Ra,I0
     GR=3.2;    %减速比
     Etagear=0.95;%减速器效率
     nstep=1;
     np=99;
     %螺旋桨转速步长
     for  i=100:1:20000 %求稳定转速的循环
            np=np+nstep;
            nm=np*GR;%马达转速
            %马达在nm转速下可提供的功率计算语句
            E=nm/Kv;
            Ia=(Vin-E)/Ra;
            Tout=Kt*(Ia-I0);
            Poffer=Tout*(2*pi*nm/60)*Etagear;
            %马达在nm转速下可提供的功率计算语句结束
            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            %螺旋桨在np转速下所需的驱动功率计算语句
            Jp=60*vf/(np*Dp);
            kk=1;
           
               while Jp>Jar(kk+1)%内部小循环得到Jp的所在区间
                     kk=kk+1;
               end
           
            Cpdrive=CP(kk)+(Jp-Jar(kk))*(CP(kk+1)-CP(kk))/(Jar(kk+1)-Jar(kk));
            Pdrive=Cpdrive*RHO*((np/60)^3)*(Dp^5);
            %螺旋桨在np转速下所需的驱动功率计算语句结束
            DetaP=Poffer-Pdrive;
            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            %根据DetaP的大小和正负确定npstep的值
            if (DetaP<1)
                nstep=0.01;
            end
            
            
            if (DetaP<0.001)
               
                   break
               
            end
     end
     Etap=ETA(kk)+(Jp-Jar(kk))*(ETA(kk+1)-ETA(kk))/(Jar(kk+1)-Jar(kk));
     Pky=Poffer*Etap;
     T=Pky/vf;