平面叶栅吹风实验

最新推荐文章于 2023-07-11 14:35:21 发布
最新推荐文章于 2023-07-11 14:35:21 发布
阅读量497 收藏 2
点赞数
分类专栏： 杂货铺 文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/dasdkhaksd/article/details/115676293
版权

目录

clc;clear;
k = 1.4; %空气绝热指数
R = 287.13; %气体常数(J/(kg*K))
matrix_P_in = zeros(41,13); %记录进口总压(kPa)
matrix_P_out = zeros(41,13); %记录出口总压(kPa)
matrix_Angle_in = zeros(41,13); %记录进气角度(°)
matrix_P0 = zeros(41,13); %记录大气压(kPa)
matrix_T_in = zeros(41,13); %记录进气温度(℃)
matrix_P2 = zeros(41,13); %记录测点2表压(kPa)
matrix_P3 = zeros(41,13); %记录测点3表压(kPa)
matrix_P4 = zeros(41,13); %记录测点4表压(kPa)
matrix_P5 = zeros(41,13); %记录测点5表压(kPa)
matrix_P6 = zeros(41,13); %记录测点6表压(kPa)
matrix_P7 = zeros(41,13); %记录测点7表压(kPa)
matrix_P9 = zeros(41,13); %记录测点9表压(kPa)
matrix_P10 = zeros(41,13); %记录测点10表压(kPa)
matrix_P11 = zeros(41,13); %记录测点11表压(kPa)
matrix_P13 = zeros(41,13); %记录测点13表压(kPa)
matrix_P14 = zeros(41,13); %记录测点14表压(kPa)
matrix_P15 = zeros(41,13); %记录测点15表压(kPa)
matrix_P16 = zeros(41,13); %记录测点16表压(kPa)
matrix_P17 = zeros(41,13); %记录测点17表压(kPa)
matrix_P19 = zeros(41,13); %记录测点19表压(kPa)
vecter_arcLength = [7.65 11.58 15.68 19.66 24.295 26.915 32.655 35.395 37.96 41.61 ...
                    45.96 49.27 51.906 55.565 63.978]; %记录15个测点对应的弧长(mm)
vecter_P = zeros(1,15); %记录15个测点的表压值(kPa)
vecter_sigma = zeros(1,15); %记录15个测点的标准误差值
vecter_coefficient_P = zeros(1,15); %记录15个测点的压力系数值
vecter_relative_arcLength = zeros(1,15); %记录15个测点的相对弧长值
vecter_Angle_probe = [38.3 35.1 31.2 30.4 30.4 30.7 32.2 31.2 31.0 31.0 32.1 35.5 36.5]; %记录探针角度(°)

%读取数据
for i = 0:12
    filename = ['t',num2str(i),'.xlsx'];
    matrix_P_in(:,i+1) = xlsread(filename,'B5:B45');
    matrix_P_out(:,i+1) = xlsread(filename,'U5:U45');
    matrix_Angle_in(:,i+1) = xlsread(filename,'W5:W45');
    matrix_P0(:,i+1) = xlsread(filename,'V5:V45');
    matrix_T_in(:,i+1) = xlsread(filename,'X5:X45');
    matrix_P2(:,i+1) = xlsread(filename,'C5:C45');
    matrix_P3(:,i+1) = xlsread(filename,'D5:D45');
    matrix_P4(:,i+1) = xlsread(filename,'E5:E45');
    matrix_P5(:,i+1) = xlsread(filename,'F5:F45');
    matrix_P6(:,i+1) = xlsread(filename,'G5:G45');
    matrix_P7(:,i+1) = xlsread(filename,'H5:H45');
    matrix_P9(:,i+1) = xlsread(filename,'J5:J45');
    matrix_P10(:,i+1) = xlsread(filename,'K5:K45');
    matrix_P11(:,i+1) = xlsread(filename,'L5:L45');
    matrix_P13(:,i+1) = xlsread(filename,'N5:N45');
    matrix_P14(:,i+1) = xlsread(filename,'O5:O45');
    matrix_P15(:,i+1) = xlsread(filename,'P5:P45');
    matrix_P16(:,i+1) = xlsread(filename,'Q5:Q45');
    matrix_P17(:,i+1) = xlsread(filename,'R5:R45');
    matrix_P19(:,i+1) = xlsread(filename,'S5:S45');
end

%计算平均值
for i = 1:41
    for j = 1:13
        if(matrix_P_in(i,j) == 0)
            matrix_P_in(i,j) = NaN;
        end
        if(matrix_P_out(i,j) == 0)
            matrix_P_out(i,j) = NaN;
        end
        if(matrix_Angle_in(i,j) == 0)
            matrix_Angle_in(i,j) = NaN;
        end
        if(matrix_P0(i,j) == 0)
            matrix_P0(i,j) = NaN;
        end
        if(matrix_T_in(i,j) == 0)
            matrix_T_in(i,j) = NaN;
        end
        if(matrix_P2(i,j) == 0)
            matrix_P2(i,j) = NaN;
        end
        if(matrix_P3(i,j) == 0)
            matrix_P3(i,j) = NaN;
        end
        if(matrix_P4(i,j) == 0)
            matrix_P4(i,j) = NaN;
        end
        if(matrix_P5(i,j) == 0)
            matrix_P5(i,j) = NaN;
        end
        if(matrix_P6(i,j) == 0)
            matrix_P6(i,j) = NaN;
        end
        if(matrix_P7(i,j) == 0)
            matrix_P7(i,j) = NaN;
        end
        if(matrix_P9(i,j) == 0)
            matrix_P9(i,j) = NaN;
        end
        if(matrix_P10(i,j) == 0)
            matrix_P10(i,j) = NaN;
        end
        if(matrix_P11(i,j) == 0)
            matrix_P11(i,j) = NaN;
        end
        if(matrix_P13(i,j) == 0)
            matrix_P13(i,j) = NaN;
        end
        if(matrix_P14(i,j) == 0)
            matrix_P14(i,j) = NaN;
        end
        if(matrix_P15(i,j) == 0)
            matrix_P15(i,j) = NaN;
        end
        if(matrix_P16(i,j) == 0)
            matrix_P16(i,j) = NaN;
        end
        if(matrix_P17(i,j) == 0)
            matrix_P17(i,j) = NaN;
        end
        if(matrix_P19(i,j) == 0)
            matrix_P19(i,j) = NaN;
        end
    end
end

vecter_P_in = nanmean(matrix_P_in);
sigma_P_in = sqrt( mean( (vecter_P_in - mean(vecter_P_in,2)).^2,2 ) * length(vecter_P_in) ...
             / (length(vecter_P_in) - 1) );
P_in = mean(vecter_P_in,2);
disp('(1) 13次测量进口总压(kPa)：');
disp(vecter_P_in);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\n',P_in,sigma_P_in);

vecter_P_out = nanmean(matrix_P_out);
sigma_P_out = sqrt( mean( (vecter_P_out - mean(vecter_P_out,2)).^2,2 ) * length(vecter_P_out) ...
              / (length(vecter_P_out) - 1) );
P_out = mean(vecter_P_out,2);
disp('(2) 13次测量出口总压(kPa)：');
disp(vecter_P_out);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\n',P_out,sigma_P_out);

vecter_Angle_in = nanmean(matrix_Angle_in);
sigma_Angle_in = sqrt( mean( (vecter_Angle_in - mean(vecter_Angle_in,2)).^2,2 ) * length(vecter_Angle_in) ...
                 / (length(vecter_Angle_in) - 1) );
Angle_in = mean(vecter_Angle_in,2);
disp('(3) 13次测量进气角度(°)：');
disp(vecter_Angle_in);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\n',Angle_in,sigma_Angle_in);

vecter_P0 = nanmean(matrix_P0);
sigma_P0 = sqrt( mean( (vecter_P0 - mean(vecter_P0,2)).^2,2 ) * length(vecter_P0) ...
           / (length(vecter_P0) - 1) );
P0 = mean(vecter_P0,2);
disp('(4) 13次测量大气压(kPa)：');
disp(vecter_P0);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\n',P0,sigma_P0);

vecter_T_in = nanmean(matrix_T_in);
sigma_T_in = sqrt( mean( (vecter_T_in - mean(vecter_T_in,2)).^2,2 ) * length(vecter_T_in) ...
             / (length(vecter_T_in) - 1) );
T_in = mean(vecter_T_in,2);
disp('(5) 13次测量进气温度(℃)：');
disp(vecter_T_in);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\n',T_in,sigma_T_in);

vecter_P2 = nanmean(matrix_P2);
vecter_P(1) = mean(vecter_P2,2);
vecter_sigma(1) = sqrt( mean( (vecter_P2 - mean(vecter_P2,2)).^2,2 ) * length(vecter_P2) ...
                  / (length(vecter_P2) - 1) );
vecter_coefficient_P(1) = vecter_P(1) / ( P_in - P0 );
vecter_relative_arcLength(1) = vecter_arcLength(1) / 72.2;
disp('(6) 13次测量测点2表压(kPa)：');
disp(vecter_P2);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(1),vecter_sigma(1),vecter_coefficient_P(1),vecter_arcLength(1),vecter_relative_arcLength(1));

vecter_P3 = nanmean(matrix_P3);
vecter_P(2) = mean(vecter_P3,2);
vecter_sigma(2) = sqrt( mean( (vecter_P3 - mean(vecter_P3,2)).^2,2 ) * length(vecter_P3) ...
                  / (length(vecter_P3) - 1) );
vecter_coefficient_P(2) = vecter_P(2) / ( P_in - P0 );
vecter_relative_arcLength(2) = vecter_arcLength(2) / 72.2;
disp('(7) 13次测量测点3表压(kPa)：');
disp(vecter_P3);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(2),vecter_sigma(2),vecter_coefficient_P(2),vecter_arcLength(2),vecter_relative_arcLength(2));

vecter_P4 = nanmean(matrix_P4);
vecter_P(3) = mean(vecter_P4,2);
vecter_sigma(3) = sqrt( mean( (vecter_P4 - mean(vecter_P4,2)).^2,2 ) * length(vecter_P4) ...
                  / (length(vecter_P4) - 1) );
vecter_coefficient_P(3) = vecter_P(3) / ( P_in - P0 );
vecter_relative_arcLength(3) = vecter_arcLength(3) / 72.2;
disp('(8) 13次测量测点4表压(kPa)：');
disp(vecter_P4);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(3),vecter_sigma(3),vecter_coefficient_P(3),vecter_arcLength(3),vecter_relative_arcLength(3));
    
vecter_P5 = nanmean(matrix_P5);
vecter_P(4) = mean(vecter_P5,2);
vecter_sigma(4) = sqrt( mean( (vecter_P5 - mean(vecter_P5,2)).^2,2 ) * length(vecter_P5) ...
                  / (length(vecter_P5) - 1) );
vecter_coefficient_P(4) = vecter_P(4) / ( P_in - P0 );
vecter_relative_arcLength(4) = vecter_arcLength(4) / 72.2;
disp('(9) 13次测量测点5表压(kPa)：');
disp(vecter_P5);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(4),vecter_sigma(4),vecter_coefficient_P(4),vecter_arcLength(4),vecter_relative_arcLength(4));
    
vecter_P6 = nanmean(matrix_P6);
vecter_P(5) = mean(vecter_P6,2);
vecter_sigma(5) = sqrt( mean( (vecter_P6 - mean(vecter_P6,2)).^2,2 ) * length(vecter_P6) ...
                  / (length(vecter_P6) - 1) );
vecter_coefficient_P(5) = vecter_P(5) / ( P_in - P0 );
vecter_relative_arcLength(5) = vecter_arcLength(5) / 72.2;
disp('(10) 13次测量测点6表压(kPa)：');
disp(vecter_P6);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(5),vecter_sigma(5),vecter_coefficient_P(5),vecter_arcLength(5),vecter_relative_arcLength(5));
    
vecter_P7 = nanmean(matrix_P7);
vecter_P(6) = mean(vecter_P7,2);
vecter_sigma(6) = sqrt( mean( (vecter_P7 - mean(vecter_P7,2)).^2,2 ) * length(vecter_P7) ...
                  / (length(vecter_P7) - 1) );
vecter_coefficient_P(6) = vecter_P(6) / ( P_in - P0 );
vecter_relative_arcLength(6) = vecter_arcLength(6) / 72.2;
disp('(11) 13次测量测点7表压(kPa)：');
disp(vecter_P7);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(6),vecter_sigma(6),vecter_coefficient_P(6),vecter_arcLength(6),vecter_relative_arcLength(6));
    
vecter_P9 = nanmean(matrix_P9);
vecter_P(7) = mean(vecter_P9,2);
vecter_sigma(7) = sqrt( mean( (vecter_P9 - mean(vecter_P9,2)).^2,2 ) * length(vecter_P9) ...
                  / (length(vecter_P9) - 1) );
vecter_coefficient_P(7) = vecter_P(7) / ( P_in - P0 );
vecter_relative_arcLength(7) = vecter_arcLength(7) / 72.2;
disp('(12) 13次测量测点9表压(kPa)：');
disp(vecter_P9);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(7),vecter_sigma(7),vecter_coefficient_P(7),vecter_arcLength(7),vecter_relative_arcLength(7));
   
vecter_P10 = nanmean(matrix_P10);
vecter_P(8) = mean(vecter_P10,2);
vecter_sigma(8) = sqrt( mean( (vecter_P10 - mean(vecter_P10,2)).^2,2 ) * length(vecter_P10) ...
                  / (length(vecter_P10) - 1) );
vecter_coefficient_P(8) = vecter_P(8) / ( P_in - P0 );
vecter_relative_arcLength(8) = vecter_arcLength(8) / 72.2;
disp('(13) 13次测量测点10表压(kPa)：');
disp(vecter_P10);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(8),vecter_sigma(8),vecter_coefficient_P(8),vecter_arcLength(8),vecter_relative_arcLength(8));
    
vecter_P11 = nanmean(matrix_P11);
vecter_P(9) = mean(vecter_P11,2);
vecter_sigma(9) = sqrt( mean( (vecter_P11 - mean(vecter_P11,2)).^2,2 ) * length(vecter_P11) ...
                  / (length(vecter_P11) - 1) );
vecter_coefficient_P(9) = vecter_P(9) / ( P_in - P0 );
vecter_relative_arcLength(9) = vecter_arcLength(9) / 72.2;
disp('(14) 13次测量测点11表压(kPa)：');
disp(vecter_P11);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(9),vecter_sigma(9),vecter_coefficient_P(9),vecter_arcLength(9),vecter_relative_arcLength(9));
    
vecter_P13 = nanmean(matrix_P13);
vecter_P(10) = mean(vecter_P13,2);
vecter_sigma(10) = sqrt( mean( (vecter_P13 - mean(vecter_P13,2)).^2,2 ) * length(vecter_P13) ...
                  / (length(vecter_P13) - 1) );
vecter_coefficient_P(10) = vecter_P(10) / ( P_in - P0 );
vecter_relative_arcLength(10) = vecter_arcLength(10) / 72.2;
disp('(15) 13次测量测点13表压(kPa)：');
disp(vecter_P13);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(10),vecter_sigma(10),vecter_coefficient_P(10),vecter_arcLength(10),vecter_relative_arcLength(10));
    
vecter_P14 = nanmean(matrix_P14);
vecter_P(11) = mean(vecter_P14,2);
vecter_sigma(11) = sqrt( mean( (vecter_P14 - mean(vecter_P14,2)).^2,2 ) * length(vecter_P14) ...
                  / (length(vecter_P14) - 1) );
vecter_coefficient_P(11) = vecter_P(11) / ( P_in - P0 );
vecter_relative_arcLength(11) = vecter_arcLength(11) / 72.2;
disp('(16) 13次测量测点14表压(kPa)：');
disp(vecter_P14);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(11),vecter_sigma(11),vecter_coefficient_P(11),vecter_arcLength(11),vecter_relative_arcLength(11));
    
vecter_P15 = nanmean(matrix_P15);
vecter_P(12) = mean(vecter_P15,2);
vecter_sigma(12) = sqrt( mean( (vecter_P15 - mean(vecter_P15,2)).^2,2 ) * length(vecter_P15) ...
                  / (length(vecter_P15) - 1) );
vecter_coefficient_P(12) = vecter_P(12) / ( P_in - P0 );
vecter_relative_arcLength(12) = vecter_arcLength(12) / 72.2;
disp('(17) 13次测量测点15表压(kPa)：');
disp(vecter_P15);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(12),vecter_sigma(12),vecter_coefficient_P(12),vecter_arcLength(12),vecter_relative_arcLength(12));
    
vecter_P16 = nanmean(matrix_P16);
vecter_P(13) = mean(vecter_P16,2);
vecter_sigma(13) = sqrt( mean( (vecter_P16 - mean(vecter_P16,2)).^2,2 ) * length(vecter_P16) ...
                  / (length(vecter_P16) - 1) );
vecter_coefficient_P(13) = vecter_P(13) / ( P_in - P0 );
vecter_relative_arcLength(13) = vecter_arcLength(13) / 72.2;
disp('(18) 13次测量测点16表压(kPa)：');
disp(vecter_P15);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(13),vecter_sigma(13),vecter_coefficient_P(13),vecter_arcLength(13),vecter_relative_arcLength(13));
    
vecter_P17 = nanmean(matrix_P17);
vecter_P(14) = mean(vecter_P17,2);
vecter_sigma(14) = sqrt( mean( (vecter_P17 - mean(vecter_P17,2)).^2,2 ) * length(vecter_P17) ...
                  / (length(vecter_P17) - 1) );
vecter_coefficient_P(14) = vecter_P(14) / ( P_in - P0 );
vecter_relative_arcLength(14) = vecter_arcLength(14) / 72.2;
disp('(19) 13次测量测点17表压(kPa)：');
disp(vecter_P17);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(14),vecter_sigma(14),vecter_coefficient_P(14),vecter_arcLength(14),vecter_relative_arcLength(14));
    
vecter_P19 = nanmean(matrix_P19);
vecter_P(15) = mean(vecter_P19,2);
vecter_sigma(15) = sqrt( mean( (vecter_P19 - mean(vecter_P19,2)).^2,2 ) * length(vecter_P19) ...
                  / (length(vecter_P19) - 1) );
vecter_coefficient_P(15) = vecter_P(15) / ( P_in - P0 );
vecter_relative_arcLength(15) = vecter_arcLength(15) / 72.2;
disp('(20) 13次测量测点19表压(kPa)：');
disp(vecter_P19);
fprintf('\t\t1.平均值为：%8.4fkPa\n\t\t2.标准误差为：%8.4f\n\t\t3.压力系数为：%8.4f\n\t\t4.弧长为：%8.4fmm\n\t\t5.相对弧长为：%8.4f\n\n',...
        vecter_P(15),vecter_sigma(15),vecter_coefficient_P(15),vecter_arcLength(15),vecter_relative_arcLength(15));

%叶形表面压力分布曲线
figure;
subplot(221);
hold on;
axis tight;
grid on;
plot(vecter_relative_arcLength,vecter_coefficient_P,'pr','markersize',10);
x = min(vecter_relative_arcLength):0.0005:max(vecter_relative_arcLength);
y = interp1(vecter_relative_arcLength,vecter_coefficient_P,x,'spline');
plot(x,y,'b','linewidth',2);
title('叶形表面压力分布曲线');
xlabel('相对弧长');
ylabel('压力系数');
legend('数据点','拟合曲线');

%叶栅出口气流角沿栅距的变化曲线
vecter_Angle_out = 90 - Angle_in + vecter_Angle_probe; %出口气流角
disp('(21) 出口气流角(°)：');
disp(vecter_Angle_out);
vecter_relative_gridSpacing = [0 2 4 6 8 10 12 14 16 18 20 22 24] / 22.65; %相对栅距
disp('(22) 相对栅距(°)：');
disp(vecter_relative_gridSpacing);
subplot(222);
hold on;
axis tight;
grid on;
plot(vecter_relative_gridSpacing,vecter_Angle_out,'pr','markersize',10);
x = min(vecter_relative_gridSpacing):0.0005:max(vecter_relative_gridSpacing);
y = interp1(vecter_relative_gridSpacing,vecter_Angle_out,x,'spline');
plot(x,y,'b','linewidth',2);
title('叶栅出口气流角沿栅距的变化曲线');
xlabel('相对栅距');
ylabel('出口气流角(°)');
legend('数据点','拟合曲线');

%损失系数沿相对栅距的变化曲线
loss_coefficient = ( (P0./(vecter_P_out+P0)).^((k-1)/k) - (P0./(vecter_P_in+P0)).^((k-1)/k) ) ./ ( 1 - (P0./(vecter_P_in+P0)).^((k-1)/k) ); %损失系数
for i = 1:13
    if(loss_coefficient(i) < 0)
        loss_coefficient(i) = 0;
    end
end
disp('(23) 损失系数：');
disp(loss_coefficient);
subplot(223);
hold on;
axis tight;
grid on;
plot(vecter_relative_gridSpacing,loss_coefficient,'pr','markersize',10);
x = min(vecter_relative_gridSpacing):0.0005:max(vecter_relative_gridSpacing);
y = interp1(vecter_relative_gridSpacing,loss_coefficient,x,'spline');
plot(x,y,'b','linewidth',2);
title('损失系数沿相对栅距的变化曲线');
xlabel('相对栅距');
ylabel('损失系数');
legend('数据点','拟合曲线');

%叶栅出口气流速度沿相对栅距的变化曲线
vecter_V_out = sqrt(1-loss_coefficient).*sqrt( 2*k*R*T_in/(k-1).*(1-(P0./(P0+vecter_P_in))).^((k-1)/k) ); %出口气流速度
disp('(24) 出口气流速度(m/s)：');
disp(vecter_V_out);
subplot(224);
hold on;
axis tight;
grid on;
plot(vecter_relative_gridSpacing,vecter_V_out,'pr','markersize',10);
x = min(vecter_relative_gridSpacing):0.0005:max(vecter_relative_gridSpacing);
y = interp1(vecter_relative_gridSpacing,vecter_V_out,x,'spline');
plot(x,y,'b','linewidth',2);
title('叶栅出口气流速度沿相对栅距的变化曲线');
xlabel('相对栅距');
ylabel('叶栅出口气流速度(m/s)');
legend('数据点','拟合曲线');
wibibaboo
关注
专栏目录
fluent二维叶型仿真_CFX案例 | 混合器流场仿真——热和流动
weixin_39862484的博客
10-22 2932
点击上方蓝字，关注并设为星标\ 知识分享 · 精品课程 · 工程仿真 \本文由iCFD授权转载CFX和Fluent都是ANSYS旗下专门用于流体力学仿真的两个软件。能够同时被ANSYS保留下来，他们在流体仿真方面是有其各自优点的。由于Fluent的普及度和市场占有率非常大，是大哥大，这里就不介绍了。下面说说CFX的一些亮点：CFX采用基于有限元的有限体积法，推出全隐式多网格耦合算...
CFX计算环形叶栅时，圆柱面展开为平面的Matlab 程序
Kevin的小屋
05-05 4272
在CFD后处理过程中，有时需要将一个圆柱面上的参数展开为一个平面并导入到Tecplot中进行出云图或流线图，这就需要坐标变换，下面的程序就是采用matlab编写的柱面展开程序。如果有建议，欢迎大家共同交流。 %% 压力数据展开 % 数据读取 clc clear all close all cyl_data = csvread('cyl_data_Pre_50H.csv'); % 读取CFX中E...
发动机叶栅的实验研究和数值预测 (2005年)
05-14
在暂冲式跨音速平面透平叶栅风洞上对某叶栅进行了吹风试验。通过数值预估实验叶栅流场,提供了流场大体结构,设计了实验方案。为选择压力传感器量程和布置测点位置提供了依据。实验测量得到了设计安装角及变工况下叶片表面压力、叶片表面等熵马赫数及波系结构随叶栅进口气流角的变化,并对叶栅内部及出口流场进行了纹影照相。数值预测与试验结果相比较,在大部分情况下,吻合得很好,说明了计算程序的可靠性。反过来,试验结果又可用于检验和改进N.S方程计算程序。
ICEM平面叶栅网格划分
最新发布
qq_41817886的博客
07-11 365
定义各边网格点数目及边上网格大小，保证各块之间网格过渡连续，检查网格质量并继续调整拓扑顶点位置，必要时可继续分割块，直至网格质量达到要求，最后输出网格。定义叶栅几何，文件开头写上【数据点个数，1】，表明是由N个数据点构成的1条曲线。分割边，并移动分割点，使得叶片Block和叶片紧贴。依次点击blocking顶点和对应的几何顶点。导入几何数据，并定义组成计算域的点、线。关联Blocking和几何，关联点。定义Blocking周期性。分割块，生成O型block。分割并关联分割后的顶点。画O型网格，继续分割。
航空发动机-网格技术篇1之叶栅流域正交画法
mp615739004的博客
06-03 3639
航空发动机-网格技术篇1之叶栅流域正交画法1 本篇目的2几何和边界条件简介2.1叶栅几何特征2.2边界条件简介3结构化画法简介4网格划分结果展示5总结与展望 1 本篇目的 本篇主要介绍一种ICEM软件划分叶片式(叶孔式也适用)叶栅流体域的结构化划分方法。划分要求如下： 尽可能保证流道的网格正交性； 在叶栅周围生成适体边界层。 本篇使用的软件为CAXA CAD电子图板2016，ICEM 2016和...
高性能亚音速平面叶栅风洞设计 (2010年)
05-29
论述了用于测量压气机和涡轮叶片在亚声速状态下的二维气动特性的平面叶栅风洞设 计问题。风洞采取暂冲式直流型式,由压缩机和储气罐提供高压气流,稳定段前采用大开角扩散 段,并在扩散段进口段设置整流锥,锥后加装两层...
平面叶栅轴流涡轮设计冲角的确定 (2009年)
06-18
通过分析最高效率点的实际液流平均进口方向,提出了流道中线决定实际设计的液流平均流动方向的观点,分析了实际机械性能与设计存在偏差的原因,较好解决了涡轮叶片在传统设计过程中冲角选择的问题,从而提高了涡轮设计的...
NUMECA FINE/Turbo™6.x 平面叶栅培训教程
06-03
非常实用的numeca教学案例，numeca越来越受欢迎了，相信这篇教程肯定会受到大家的喜欢，尤其对于初学者，希望给大家带来方便
采用平面叶栅模拟压气机动叶叶尖间隙流.pdf
07-06
本研究采用了平面叶栅作为实验模型来模拟压气机动叶叶尖间隙流。通过对比三种不同工况下的实验结果，探讨了不同条件下叶尖间隙流的特性及其对整体性能的影响。具体来说： - **工况1**：平面叶栅端壁面静止，即没有...
墙角导热问题数值解
Aiden Lee's blog
11-30 3575
墙角导热问题数值解(1) 物理问题(2) 用热平衡法建立节点代数方程(2.1) 内节点(2.2) 外部平直边界节点(2.3) 外部角点(2.4) 内部角点(3) 高斯-赛德尔(Gauss-Seidel)迭代法(4) 显化节点方程(4.1) 等温边界条件(4.2) 对流边界条件(5) 导热量换热面积选取(6) FORTRAN迭代(6.1) 等温边界条件(6.2) 对流边界条件(7) MATLAB数据可视化 (1) 物理问题 已知一个建筑物墙角如图所示，其中墙体为二维、稳态、无内热源、常物性导热。现在要求用数
平板边界层内的流速分布实验
Aiden Lee's blog
03-19 2462
目录(1) 代码(2) 运行结果 (1) 代码 clc;clear; Rg = 287; %气体常数（J/(kg·K)）。 g = 9.8; %重力加速度（m/s^2）。 waterDensity = 1000; %水的密度（kg/m^3）。 %对已知实验参数的输入。 atmosphericPressure = input('请输入大气压：(Pa)'); temperature = input('请输入温度：(℃)'); KVC = input('请输入流体运动粘性系数：(10^-7m^2/s)');
对流换热实验
Aiden Lee's blog
03-19 1563
目录(1) 空气横掠单管强制对流换热实验(1.1) 代码(1.2) 运行结果(2) 水平管外自然对流换热实验(2.1) 代码(2.2) 运行结果 (1) 空气横掠单管强制对流换热实验 (1.1) 代码 clc;clear; d = 3.73e-3; l = 99.6e-3; tw = [40.57 42.41 44.27 46.48 49.72]; tf = [21.12 21.59 21.77 21.77 21.99]; deltaP = [213.5 180.3 144.0 108.5 74.4]; U
压力表校验及压力变送器标定实验
Aiden Lee's blog
04-02 1520
目录(1) 代码(2) 运行结果 (1) 代码 %压力表数据处理 clc;clear; %大气压(MPa) p0 = 0.096; %活塞式压力机标准压力值(MPa) pStd_Upstroke = [0 3 6 9 12]; pStd_Downstroke = [9 6 3 0]; %压力表第一组读数(MPa) p1_Upstroke = [0.0 3.2 6.2 9.2 12.2]; p1_Downstroke = [9.2 6.2 3.2 0.0]; %压力表第二组读数(MPa) p2_Upstrok
Fortran集成VS2019
Aiden Lee's blog
05-20 1170
只下载VS2019是不能编译Fortran的，还需要下载另外两个东西。建议安装community版本，适合个人使用，并且免费。
零维系统热自燃问题求解
Aiden Lee's blog
01-15 1014
目录(1) 物理问题(2) 求解目标(3) Euler法求解微分方程(4) Fortran求解微分方程(5) Matlab数据可视化 (1) 物理问题 热平衡方程： VρcvdTdτ=Q1−Q2V{\rho}c_{v}\frac{dT}{d\tau}=Q_{1}-Q_{2}Vρcv​dτdT​=Q1​−Q2​ 产热Q1Q_{1}Q1​为： Q1=k0e−ERT⋅Cn⋅V⋅QQ_{1}=k_{0}e^{-\frac{E}{RT}}{\cdot}C^{n}{\cdot}V{\cdot}QQ1​=k0​e−R
换热器综合实验
Aiden Lee's blog
03-19 888
目录(1) 代码(2) 运行结果 (1) 代码 clc;clear; heatExchangeArea = 0.14; %换热面积（m^2）。 flowArea = (pi / 4) * ((16 - 1.5 * 2) * 0.001)^2; %内管流体横截面面积（m^2）。 %顺流。 parallelFlowTi1 = [60.2 59.4 60.0 60.1 60.0]; %热流体进口温度（℃）。 parallelFlowTo1 = [54.0 52.7 52.4 51.6 50.7]; %热流体出
CSDN博客排版技巧
Aiden Lee's blog
03-14 819
目录(1) 基本排版操作(1.1) 添加目录(1.2) 添加标题(1.3) 调整图片大小(1.4) 插入gif图片(2) LateX语法(2.1) 希腊字母(2.2) 数学运算符(2.3) 上下标(2.4) 特殊符号(2.5) 函数符号 (1) 基本排版操作 (1.1) 添加目录 @[TOC](目录名称) (1.2) 添加标题 n个# + 空格 + 标题名称 (1.3) 调整图片大小 原图：直接复制过来，不做改变。 按原图比例缩放：在图片链接后面加上“空格=尺寸x”。 按固定大小缩放：在图片链接后
高性能亚音速平面叶栅风洞设计及其关键技术
"高性能亚音速平面叶栅风洞设计 (2010年) - 论文 - 工程技术" 高性能亚音速平面叶栅风洞是专门用于研究压气机和涡轮叶片在亚声速状态下二维气动特性的关键设施。这种风洞设计采用了暂冲式直流型式，其工作原理是...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值