日常学习练手—1045螺旋桨的气动定常分析

clear all
close all
clc
%%  导入
v_first=importdata('v_135.txt');
%V_5=importdata('E:\maker_BBS\2016_6_25\cfx\8_20\v_5.txt');
%r_80100=importdata('E:\maker_BBS\2016_6_25\cfx\8_20\r_80100.txt');
lamda=importdata('lamda.txt');
%%  整理
force_v_1=v_first(1,:);force_v_3=v_first(2,:);force_v_5=v_first(3,:);
turque_v_1=-v_first(4,:);turque_v_3=-v_first(5,:);turque_v_5=-v_first(6,:);
%F_v_5=V_5(1,:);T_v_5=V_5(2,:);
%force_r_80=r_80100([1:6],1);force_r_100=r_80100([1:10],2);
%turque_r_80=r_80100([11:16],1);turque_r_100=r_80100([11:20],2);
%%  常数
%v=(1:2:19);     %%风速m/s   
%r=(40:20:460);  %%转速r/s
D=0.254;        %%直径m
ro=1.205;
%%  计算第一个文件夹
    %%  前进比
    r_135=[40:20:120];
    la_v_1=lamda(1,[1:5]);la_v_3=lamda(2,[1:5]);la_v_5=lamda(3,[1:5]);
    %% 拉力系数
    f_v_1=force_v_1./(ro.*(r_135.^2).*(D^4));
    f_v_3=force_v_3./(ro.*(r_135.^2).*(D^4));
    f_v_5=force_v_5./(ro.*(r_135.^2).*(D^4));
    %% 阻力矩系数
    t_v_1=turque_v_1./(ro.*(r_135.^2).*(D^5));
    t_v_3=turque_v_3./(ro.*(r_135.^2).*(D^5));
    t_v_5=turque_v_5./(ro.*(r_135.^2).*(D^5));
    %% 总功率
    P_v_1=turque_v_1.*2*pi.*r_135;
    P_v_3=turque_v_3.*2*pi.*r_135;
    P_v_5=turque_v_5.*2*pi.*r_135;
    %% 效率
    g_v_1=f_v_1.*1./(P_v_1);
    g_v_3=f_v_3.*3./(P_v_3);
    g_v_5=f_v_5.*5./(P_v_5);
%%  绘图
figure(1)
plot(r_135,force_v_1,'-o',r_135,force_v_3,'-*',r_135,force_v_5,'-^')
title('转速与拉力')
xlabel('转速')
ylabel('拉力')
legend('v=1m/s','v=3m/s','v=5m/s');
  
figure(2)
plot(r_135,turque_v_1,'-o',r_135,turque_v_3,'-*',r_135,turque_v_5,'-^')
title('转速与阻力矩')
xlabel('转速')
ylabel('阻力矩')
legend('v=1m/s','v=3m/s','v=5m/s');
  
figure(3)
plot(r_135,P_v_1,'-o',r_135,P_v_3,'-*',r_135,P_v_5,'-^')
title('转速与总功率')
xlabel('转速')
ylabel('总功率')
legend('v=1m/s','v=3m/s','v=5m/s');
  
figure(4)
plot(r_135,g_v_1,'-o',r_135,g_v_3,'-*',r_135,g_v_5,'-^')
title('转速与效率')
xlabel('转速')
ylabel('效率')
legend('v=1m/s','v=3m/s','v=5m/s');
  
figure(5)
plot(la_v_1,g_v_1,'-o',la_v_3,g_v_3,'-*',la_v_5,g_v_5,'-^')
title('前进比与效率')
xlabel('前进比')
ylabel('效率')
legend('v=1m/s','v=3m/s','v=5m/s');