汽轮机计算MATLAB代码.docVIP

附：汽轮机课程设计MATLAB计算程序代码： 汽轮机流量计算的代码： %------------------------------------------------------------------------- %汽轮机流量的计算 P_el=2.25;%汽轮机的的设计功率 deltah_n=200; %汽轮机调节级的比焓降 h0=3224.5; h_l=h0-deltah_n; %调节级后的比焓 h_e1=2711.4; h_e2=2637.6; h_z=2429.74; h_ed=437; h_cw=316.6; h_w2=304.03; h_w1=171.73; deltahmac=h0-h_z ; %通流部分的理想比焓降 eta_ri=0.775; eta_m=0.97; eta_g=0.97; eta_h=0.98; m=1.026 D_0=(3.6*P_el*1000/(deltahmac*eta_g*eta_m)*m)/0.97 %汽轮机的总进气量 deltaD_1=0.03*D_0 D_fw=D_0-deltaD_1 deltaD_e1=-D_fw*(-h_ed+h_cw)/(h_e1-h_cw) deltaD_cw=D_fw*(h_e1-h_ed)/(h_e1-h_cw) deltaD_e2=deltaD_cw*(h_w2-h_w1)/((h_e2-h_cw)*eta_h) %调节级 P_i0=D_0*deltah_n/3.6 %第一级组 D_1=D_0-deltaD_1 P_i1=D_1*(h_l-h_e1)/3.6 %第二级组 D_2=D_1-deltaD_e1 P_i2=D_2*(h_e1-h_e2)/3.6 %第三级组 D_3=D_2-deltaD_e2 P_i3=D_3*(h_e2-h_z)/3.6 %整机功率 P_i=P_i0+P_i1+P_i2+P_i3 %机械损失 deltaP_m=P_i*(1-eta_m) %汽轮机轴端功率 P_a=P_i-deltaP_m %发电机功率 P_e=P_a*eta_g %误差计算 erro1=P_el*1000-P_e erro=(P_el-P_e/1000)/P_el*100 %------------------------------------------------------------------------ %汽耗率 d=D_0*1000/P_e d_1=D_0*1000/(D_0*(h0-h_z)/3.6-deltaP_m)/eta_g q=d*(h0-h_ed) %汽轮机装置的绝对电效率 eta_el=3600/q %------------------------------------------------------------------------- 汽轮机调节级热力计算的代码： %汽轮机调节级的计算 %-------------------------------------------------------------------------- n=5600; D_0=11.36; p_0=2.256; h_0=3220.17; v_0=0.1259; deltah_t=200; x_a=0.24; aerfa_1=14; oum_b=0; oum_gb=0.03; oum_b1=0.06; delta_t=1; miu_t=0.6; e_c=0.4 %-------------------------------------------------------------------------- oum_m=oum_b+oum_gb+oum_b1 deltah_n=deltah_t*(1-oum_m) h_1t=h_0-deltah_n p_1=input(喷嘴出口压力) v_1t=input(喷嘴出口比体积) efthelen=p_1/p_0 c_1t=sqrt(2000*deltah_n) c_1=0.97*c_1t k=1.3 sinaerfa1_jia_s1=(2/(k+1))^(1/(k-1))*sqrt((k-1)/(k+1))/(efthelen^(1/k)*sqrt(1-efthelen^((k-1)/k)))*sin(aerfa_1/180*pi) s1=asin(sinaerfa1_jia_s1)*180/pi-aerfa_1 u=x_a*c_1 %调节级平均直径 d_m=60*x_a*c_1/(pi*n) %m %喷管出口面积 A_n=D_0/3.6/(0.648*sqrt(p_0*10^6/v_0))*10000 %cm2 e=0.2% --------------- l_n=