汽轮机计算MATLAB代码

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.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.026D_0=(3.6*P

2、_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_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

3、_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*(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_e

4、d)%汽轮机装置的绝对电效率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_b1deltah_n=deltah_t*(1-oum_m)h_1t=h_0-deltah_np_1=input(喷嘴出口压力)v_1t=input(

5、喷嘴出口比体积)efthelen=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_0/3.6/(0.648*sqrt(p_0*106/v_0)*1

6、0000 %cm2e=0.2% -l_n=A_n./(e*pi*d_m*100*sin(aerfa_1/180*pi)deltah_n_kthe=(1-0.972)*deltah_n%-%第一列动叶热力计算w_1=sqrt(c_12+u2-2*c_1*u*cos(aerfa_1+s1)/180*pi)beta_1=asin(c_1*sin(aerfa_1+s1)/180*pi)/w_1)/pi*180w_2t=w_1fai_b=-0.07/480*w_2t+0.92+0.7/24w_2=fai_b.*w_2tdelta_beta_1=input(动叶出口气流角变化，35)beta_2=beta

7、_1-delta_beta_1c_2=sqrt(w_2.2+u.2-2.*w_2.*u.*cos(beta_2/180*pi)aerfa_2=asin(w_2*sin(beta_2/180*pi)/c_2)/pi*180h_1=h_1t+deltah_n_kthe %动叶进口比焓v_1=input(动叶进口比体积，1/5.81)delta_r=0.5;delta_t=1.5;delta=delta_r+delta_tl_b11=l_n*10+delta_r+delta_tmiu_b=input(查图1-11，0.93)A_b=D_0/3.6/(miu_b*w_2t/v_1)*10000l_b1

8、=A_b./(e*pi*d_m*sin(beta_2/180*pi)/10l_b1-l_b11pausedeltah_bkthe=w_2t2/2000*(1-fai_b2)h_2=h_1+deltah_bkthev_2=input(动叶出口比体积v_2)%-%导叶热力计算deltah_gb=oum_gb*deltah_th_1t1=h_2-deltah_gbp_11=input(导叶出口压力mpa)v_11=input(导叶出口比体积)c_1t1=sqrt(2000*deltah_gb+c_22)fai_gb=input(查图1-18,fai_gb)c_11=fai_gb*c_1t1delta

9、_aerfa=input(输入deltaaerfa 510)aerfa_11=aerfa_2-delta_aerfal_gb1=l_b1+deltadeltaG_gbt=e*0.6*pi*(d_m+l_gb1/1000)*0.001*sqrt(2000*oum_gb*deltah_t)/v_11G_gb=D_0/3.6-deltaG_gbtmiu_gb=input(查图11-11)A_gb=G_gb/(miu_gb*c_1t1/v_11)*10000l_gb=A_gb/(e*pi*d_m*10*sin(aerfa_11/180*pi)deltah_gb1=c_1t12/2000*(1-fai_

10、gb2)%导叶出口实际参数h_1_shiji=h_1t1+deltah_gb1v_1_shiji=input(由实际出口焓查焓熵图实际出口比体积)%-%第二列动叶热力计算deltah_b1=oum_b1*deltah_th_2t1=h_1_shiji-deltah_b1p_21=input(动叶出口压力)v_2t1=input(动叶出口比体积)w_11=sqrt(c_112+u2-2*c_11*u*cos(aerfa_11)/180*pi)beta_11=asin(c_11*sin(aerfa_11)/180*pi)/w_11)/pi*180w_2t1=sqrt(2000*deltah_b1+

11、w_112)fai_b1=input(查图1-18)w_21=fai_b1*w_2t1delta_beta2=input(718)beta_21=beta_11-delta_beta2c_21=sqrt(w_212+u2-2*w_21*u*cos(beta_21/180*pi)aerfa_21=asin(w_21*sin(beta_21/180*pi)/c_21)*180/pideltah_b11=w_2t12/2000*(1-fai_b12)%动叶损失deltah_c21=c_212/2000h_21_shiji=h_2t1+deltah_c21v_21_shiji=input(动叶出口比体

12、积)l_b21=l_gb+deltaoum_br1=1-(1-oum_b1)*d_m/(d_m-l_b21/1000)oum_bt1=1-(1-oum_br1)*(1-l_b21/1000)/(1+l_b21/1000)deltaG_bt1=e*0.6*pi*(d_m+l_b21/1000)*0.001*sqrt(2000*oum_bt1*deltah_t)/v_2t1miu_b1=input(查图11-11)A_b1=(D_0/3.6-deltaG_bt1)/(miu_b1*w_2t1/v_2t1)*10000l_b2=A_b1/(e*pi*d_m*sin(beta_21/180*pi)/1

13、0l_b2-l_b21%-%轮周功校核P_u1_1=u*(c_1*cos(aerfa_1+s1)/180*pi)+c_11*cos(aerfa_11/180*pi)+c_2*cos(aerfa_2/180*pi)+c_21*cos(aerfa_21/180*pi)/1000deltah_u=deltah_n_kthe+deltah_bkthe+deltah_gb1+deltah_b11+deltah_c21P_u1_2=deltah_t-deltah_udelta_eta=(P_u1_2-P_u1_1)/P_u1_2*100%-%轮周效率eta_u=deltah_u/deltah_t%-%级内

14、损失计算deltaP_f=1.2*(u/100)3*d_m2*0.5*(1/v_1+1/v_21_shiji)deltah_f=3600*deltaP_f/(D_0*1000)l=(l_n+l_gb1+l_gb+l_b11+l_b1+l_b2+l_b21)/7deltah_e=2/l*P_u1_2E_0=deltah_tx_a1=u/sqrt(2000*deltah_t)kthe_w=0.55*1/e*(1-e-e_c/2)*x_a13deltah_w=kthe_w*E_0kthe_s=0.016*1/e*2/d_m*x_a1deltah_s=kthe_s*E_0deltah_e1=delta

15、h_w+deltah_sdeltahh=deltah_t-deltah_n_kthe-deltah_bkthe-deltah_gb1-deltah_b11-deltah_e-deltah_c21deltah_t1=(deltaG_gbt+deltaG_bt1)/(D_0/3.6)*deltahh%-%级的内功率deltah_i=deltahh-deltah_f-deltah_t1-deltah_e1P_i=D_0/3.6*deltah_i%-%级的内效率eta_i=deltah_i/E_0%-三、 汽轮机压力级热力计算的代码：d_m=input(级的平均直径);c_0=input(流入该级的蒸

16、汽初速);D=input(流量);p_0=input(级前蒸汽压力);x_0=input(干度);deltah_t=input(级的理想比焓);oum_m=input(级的平均反动度);oum_t=input(顶部反动度);aerfa_1=input(喷管出气角);miu_1=input(余速利用系数);%喷管的理想滞止比焓降deltah_c0=c_02/2000deltah_n_xing=(1-oum_m)*(deltah_t+deltah_c0)deltah_t_xing=deltah_t+deltah_c0deltah_b=oum_m*deltah_t_xing%喷管热力计算p_0_xi

17、ng=input(喷管前的滞止压力);h_0_xing=input(滞止比焓);v_0_xing=input(滞止比体积);h_0=input(喷管前比焓);p_1=input(喷管后压力);v_1t=input(理想比体积);h_1t=h_0-deltah_tk=1.035+0.1*x_0efthlen_cr=(2/(k+1)(k/(k-1)efthlen_n=p_1/p_0_xingpausep_cr=p_0_xing*efthlen_crh_cr=input(临界焓);v_cr=input(临界比体积);c_cr=sqrt(2000*(h_0_xing-h_cr)c_1t=sqrt(20

18、00*deltah_n_xing)fai=0.97;c_1=fai*c_1taerfa_jia_delta=asin(sin(aerfa_1/180*pi)*c_cr*v_1t/(c_1t*v_cr)/pi*180miu_p=input(0.80.7);d_p=input(汽封齿的平均直径);delta_p=input(汽封间隙);z_p=input(齿数);deltaD_p=miu_p*pi*d_p*delta_p*sqrt(2*deltah_n_xing*1000)/(v_1t*sqrt(z_p+1/z_p)D_n=D/3.6-deltaD_pA_n=D_n/(0.648*sqrt(p_0

19、_xing*106/v_0_xing) l_n=A_n/(1*pi*d_m*sin(aerfa_1/180*pi)*1000deltah_nkthe=(1-fai2)*deltah_n_xingh_1=h_1t+deltah_nkthev_1=input(喷管出口实际比体积);%-%动叶栅热力计算n=3000u=pi*d_m*n/60w_1=sqrt(c_12+u2-2*c_1*u*cos(aerfa_jia_delta/180*pi)beta_1=atan(c_1*sin(aerfa_jia_delta/180*pi)/(c_1*cos(aerfa_jia_delta/180*pi)-u)/

20、pi*180deltah_w1=w_12/2000w_2t=sqrt(2000*(deltah_b+deltah_w1)faifai=input(速度系数);w_2=faifai*w_2tbeta_2=input(beta2);c_2=sqrt(w_22+u2-2*w_2*u*cos(beta_2/180*pi)aerfa_2=asin(w_2*sin(beta_2*pi/180)/c_2)/pi*180deltah_bkthe=w_22/2000*(1-faifai2)h_2t=h_1-deltah_bv_2t=input(动叶出口理想比体积);p_2=input(动叶出口压力);efthl

21、en_b=p_2/p_1pauseh_2=h_2t+deltah_bkthev_2=input(动叶出口实际比体积);x_2=input(动叶出口实际干度);l_b1=l_n+1+2deltaD_t=1*0.6*pi*(1.442+0.075)*0.002/v_2t*sqrt(oum_t*deltah_t_xing*1000)miu_b=0.935;A_b=(D_n-deltaD_t)/(miu_b*w_2t/v_2t)d_b=input(输入直径db)l_b=A_b/(pi*d_b*sin(beta_2*pi/180)*1000deltah_c2=c_22/2000%-%轮周效率P_u=deltah_t_xing-