汽轮机计算MATLAB代码.doc

《汽轮机计算MATLAB代码.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、附：汽轮机课程设计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.0

2、3;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.026D_0=(3.6*P_el*1000/(deltahmac*eta_g*eta_m)*m)/0.97%汽轮机的总进气量deltaD_1=0.03*D_0D_fw=D_0-deltaD_1deltaD_e1=-D_fw*(-h_ed+h_cw)/(h_e1-h_cw)deltaD_cw=D_fw*(h_e1-h_ed)/(h_e1-h

3、_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_1P_i1=D_1*(h_l-h_e1)/3.6%第二级组D_2=D_1-deltaD_e1P_i2=D_2*(h_e1-h_e2)/3.6%第三级组D_3=D_2-deltaD_e2P_i3=D_3*(h_e2-h_z)/3.6%整机功率P_i=P_i0+P_i1+P_i2+P_i3%机械损失deltaP_m=P_i*(

4、1-eta_m)%汽轮机轴端功率P_a=P_i-deltaP_m%发电机功率P_e=P_a*eta_g%误差计算erro1=P_el*1000-P_eerro=(P_el-P_e/1000)/P_el*100%------------------------------------------------------------------------%汽耗率d=D_0*1000/P_ed_1=D_0*1000/(D_0*(h0-h_z)/3.6-deltaP_m)/eta_gq=d*(h0-h_ed)%汽轮机装置的绝

5、对电效率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.2

6、4;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_b1deltah_n=deltah_t*(1-oum_m)h_1t=h_0-deltah_np_1=input('喷嘴出口压力')v_1t=input('喷嘴出口比体积')efthel

7、en=p_1/p_0c_1t=sqrt(2000*deltah_n)c_1=0.97*c_1tk=1.3sinaerfa1_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_1u=x_a*c_1%调节级平均直径d_m=60*x_a*c_1/(pi*n)%m%喷管出口面积A_n=D

8、_0/3.6/(0.648*sqrt(p_0*10^6/v_0))*10000%cm2e=0.2%---------------l_n=A_n./(e*pi*d_m*100*sin(aerfa_1/180*pi))deltah_n_kthe=(1-0.97^2)*deltah_n%-------------------------