汽轮机计算MATLAB代码
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
delta_r=0.5;
delta_t=1.5;
delta=delta_r+delta_t
l_b11=l_n*10+delta_r+delta_t
miu_b=input('查图1-11，0.93')
A_b=D_0/3.6/(miu_b*w_2t/v_1)*10000
l_b1=A_b./(e*pi*d_m*sin(beta_2/180*pi))/10
l_b21=l_gb+delta
oum_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_2t1
P_u1_2=deltah_t-deltah_u
delta_eta=(P_u1_2-P_u1_1)/P_u1_2*100
%--------------------------------------------------------------------------
%轮周效率
eta_u=deltah_u/deltah_t
miu_b1=input('查图11-11')
A_b1=(D_0/3.6-deltaG_bt1)/(miu_b1*w_2t1/v_2t1)*10000
l_b2=A_b1/(e*pi*d_m*sin(beta_21/180*pi))/10
l_b2-l_b21
%--------------------------------------------------------------------------
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
%-------------------------------------------------------------------------
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
%轮周功校核
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))/1000
deltah_u=deltah_n_kthe+deltah_bkthe+deltah_gb1+deltah_b11+deltah_c21
w_2t1=sqrt(2000*deltah_b1+w_11^2)
fai_b1=input('查图1-18')
w_21=fai_b1*w_2t1
delta_beta2=input('7~18')
beta_21=beta_11-delta_beta2
c_21=sqrt(w_21^2+u^2-2*w_21*u*cos(beta_21/180*pi))
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;
deltah_n_kthe=(1-0.97^2)*deltah_n
%--------------------------------------------------------------------------
%第一列动叶热力计算
w_1=sqrt(c_1^2+u^2-2*c_1*u*cos((aerfa_1+s1)/180*pi))
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)
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
%--------------------------------------------------------------------------
c_11=fai_gb*c_1t1
delta_aerfa=input('输入deltaaerfa 5~10')
aerfa_11=aerfa_2-delta_aerfa
l_gb1=l_b1+delta
deltaG_gbt=e*0.6*pi*(d_m+l_gb1/1000)*0.001*sqrt(2000*oum_gb*deltah_t)/v_11
aerfa_21=asin((w_21*sin(beta_21/180*pi))/c_21)*180/pi
deltah_b11=w_2t1^2/2000*(1-fai_b1^2)%动叶损失
deltah_c21=c_21^2/2000
h_21_shiji=h_2t1+deltah_c21
v_21_shiji=input('动叶出口比体积')
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
deltah_b1=oum_b1*deltah_t
h_2t1=h_1_shiji-deltah_b1
p_21=input('动叶出口压力')
v_2t1=input('动叶出口比体积')
w_11=sqrt(c_11^2+u^2-2*c_11*u*cos((aerfa_11)/180*pi))
beta_11=asin(c_11*sin((aerfa_11)/180*pi)/w_11)/pi*180
二、汽轮机调节级热力计算的代码：
%汽轮机调节级的计算
%--------------------------------------------------------------------------
n=5600;
D_0=11.36;
p_0=2.256;
h_0=3220.17;
v_0=0.1259;
c_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*180
h_1=h_1t+deltah_n_kthe %动叶进口比焓
v_1=input('动叶进口比体积，1/5.81')
l=(l_n+l_gb1+l_gb+l_b11+l_b1+l_b2+l_b21)/7
deltah_e=2/l*P_u1_2
E_0=deltah_t
x_a1=u/sqrt(2000*deltah_t)
kthe_w=0.55*1/e*(1-e-e_c/2)*x_a1^3
%导叶热力计算
deltah_gb=oum_gb*deltah_t
h_1t1=h_2-deltah_gb
p_11=input('导叶出口压力mpa')
v_11=input('导叶出口比体积')
c_1t1=sqrt(2000*deltah_gb+c_2^2)
fai_gb=input('查图1-18,fai_gb')
beta_1=asin(c_1*sin((aerfa_1+s1)/180*pi)/w_1)/pi*180
w_2t=w_1
fai_b=-0.07/480*w_2t+0.92+0.7/24
w_2=fai_b.*w_2t
delta_beta_1=input('动叶出口气流角变化，3~5')
beta_2=beta_1-delta_beta_1
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
%--------------------------------------------------------------------------
%级内损失计算
deltaP_f=1.2*(u/100)^3*d_m^2*0.5*(1/v_1+1/v_21_shiji)
deltah_f=3600*deltaP_f/(D_0*1000)
l_b1-l_b11
pause
deltah_bkthe=w_2t^2/2000*(1-fai_b^2)
h_2=h_1+dFra bibliotekltah_bkthe
v_2=input('动叶出口比体积v_2')
%--------------------------------------------------------------------------
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=A_n./(e*pi*d_m*100*sin(aerfa_1/180*pi))
附：汽轮机课程设计MATLAB计算程序代码：
一、汽轮机流量计算的代码：
%-------------------------------------------------------------------------
%汽轮机流量的计算
P_el=2.25;%汽轮机的的设计功率
deltah_n=200; %汽轮机调节级的比焓降
%第二级组
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)
%汽轮机轴端功率
G_gb=D_0/3.6-deltaG_gbt
miu_gb=input('查图11-11')
A_gb=G_gb/(miu_gb*c_1t1/v_11)*10000
l_gb=A_gb/(e*pi*d_m*10*sin(aerfa_11/180*pi))
deltah_gb1=c_1t1^2/2000*(1-fai_gb^2)
epidmsinbeta2180pi10lb1lb11pausedeltahbkthew2t220001faib2h2h1deltahbkthev2input动叶出口比体积v2导叶热力计算deltahgboumgbdeltahth1t1h2deltahgbp11input导叶出口压力mpav11input导叶出口比体积c1t1sqrt2000deltahgbc22faigbinput查图118faigbc11faigbc1t1deltaaerfainput输入deltaaerfa510aerfa11aerfa2deltaaerfalgb1lb1deltadeltaggbte06pidmlgb110000001sqrt2000oumgbdeltahtv11ggbd036deltaggbtmiugbinput查图1111agbggbmiugbc1t1v1110000lgbagbepidm10sinaerfa11180pideltahgb1c1t1220001faigb2导叶出口实际参数h1shijih1t1deltahgb1v1shijiinput由实际出口焓查焓熵图实际出口比体积第二列动叶热力计算deltahb1oumb1deltahth2t1h1shijideltahb1p21input动叶出口压力v2t1input动叶出口比体积w11sqrtc112u22c11ucosaerfa11180pibeta11asinc11sinaerfa11180piw11pi180w2t1sqrt2000deltahb1w112faib1input查图118w21faib1w2t1deltabeta2input718beta21beta11deltabeta2c21sqrtw212u22w21ucosbeta21180piaerfa21asinw21sinbeta21180pic21180pideltahb11w2t1220001faib12动叶损失deltahc21c2122000h21shijih2t1deltahc21v21shijiinput动叶出口比体积lb21lgbdeltaoumbr111oumb1dmdmlb211000oumbt111oumbr11lb2110001lb211000deltagbt1e06pidmlb2110000001sqrt2000oumbt1deltahtv2t1miub1input查图1
%导叶出口实际参数
h_1_shiji=h_1t1+deltah_gb1
v_1_shiji=input('由实际出口焓查焓熵图—实际出口比体积')
delta_t=1.5;
delta=delta_r+delta_t
l_b11=l_n*10+delta_r+delta_t
miu_b=input('查图1-11，0.93')
A_b=D_0/3.6/(miu_b*w_2t/v_1)*10000
l_b1=A_b./(e*pi*d_m*sin(beta_2/180*pi))/10
l_b21=l_gb+delta
oum_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_2t1
P_u1_2=deltah_t-deltah_u
delta_eta=(P_u1_2-P_u1_1)/P_u1_2*100
%--------------------------------------------------------------------------
%轮周效率
eta_u=deltah_u/deltah_t
miu_b1=input('查图11-11')
A_b1=(D_0/3.6-deltaG_bt1)/(miu_b1*w_2t1/v_2t1)*10000
l_b2=A_b1/(e*pi*d_m*sin(beta_21/180*pi))/10
l_b2-l_b21
%--------------------------------------------------------------------------
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
%-------------------------------------------------------------------------
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
%轮周功校核
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))/1000
deltah_u=deltah_n_kthe+deltah_bkthe+deltah_gb1+deltah_b11+deltah_c21
w_2t1=sqrt(2000*deltah_b1+w_11^2)
fai_b1=input('查图1-18')
w_21=fai_b1*w_2t1
delta_beta2=input('7~18')
beta_21=beta_11-delta_beta2
c_21=sqrt(w_21^2+u^2-2*w_21*u*cos(beta_21/180*pi))
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;
deltah_n_kthe=(1-0.97^2)*deltah_n
%--------------------------------------------------------------------------
%第一列动叶热力计算
w_1=sqrt(c_1^2+u^2-2*c_1*u*cos((aerfa_1+s1)/180*pi))
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)
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
%--------------------------------------------------------------------------
c_11=fai_gb*c_1t1
delta_aerfa=input('输入deltaaerfa 5~10')
aerfa_11=aerfa_2-delta_aerfa
l_gb1=l_b1+delta
deltaG_gbt=e*0.6*pi*(d_m+l_gb1/1000)*0.001*sqrt(2000*oum_gb*deltah_t)/v_11
aerfa_21=asin((w_21*sin(beta_21/180*pi))/c_21)*180/pi
deltah_b11=w_2t1^2/2000*(1-fai_b1^2)%动叶损失
deltah_c21=c_21^2/2000
h_21_shiji=h_2t1+deltah_c21
v_21_shiji=input('动叶出口比体积')
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
deltah_b1=oum_b1*deltah_t
h_2t1=h_1_shiji-deltah_b1
p_21=input('动叶出口压力')
v_2t1=input('动叶出口比体积')
w_11=sqrt(c_11^2+u^2-2*c_11*u*cos((aerfa_11)/180*pi))
beta_11=asin(c_11*sin((aerfa_11)/180*pi)/w_11)/pi*180
二、汽轮机调节级热力计算的代码：
%汽轮机调节级的计算
%--------------------------------------------------------------------------
n=5600;
D_0=11.36;
p_0=2.256;
h_0=3220.17;
v_0=0.1259;
c_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*180
h_1=h_1t+deltah_n_kthe %动叶进口比焓
v_1=input('动叶进口比体积，1/5.81')
l=(l_n+l_gb1+l_gb+l_b11+l_b1+l_b2+l_b21)/7
deltah_e=2/l*P_u1_2
E_0=deltah_t
x_a1=u/sqrt(2000*deltah_t)
kthe_w=0.55*1/e*(1-e-e_c/2)*x_a1^3
%导叶热力计算
deltah_gb=oum_gb*deltah_t
h_1t1=h_2-deltah_gb
p_11=input('导叶出口压力mpa')
v_11=input('导叶出口比体积')
c_1t1=sqrt(2000*deltah_gb+c_2^2)
fai_gb=input('查图1-18,fai_gb')
beta_1=asin(c_1*sin((aerfa_1+s1)/180*pi)/w_1)/pi*180
w_2t=w_1
fai_b=-0.07/480*w_2t+0.92+0.7/24
w_2=fai_b.*w_2t
delta_beta_1=input('动叶出口气流角变化，3~5')
beta_2=beta_1-delta_beta_1
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
%--------------------------------------------------------------------------
%级内损失计算
deltaP_f=1.2*(u/100)^3*d_m^2*0.5*(1/v_1+1/v_21_shiji)
deltah_f=3600*deltaP_f/(D_0*1000)
l_b1-l_b11
pause
deltah_bkthe=w_2t^2/2000*(1-fai_b^2)
h_2=h_1+dFra bibliotekltah_bkthe
v_2=input('动叶出口比体积v_2')
%--------------------------------------------------------------------------
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=A_n./(e*pi*d_m*100*sin(aerfa_1/180*pi))
附：汽轮机课程设计MATLAB计算程序代码：
一、汽轮机流量计算的代码：
%-------------------------------------------------------------------------
%汽轮机流量的计算
P_el=2.25;%汽轮机的的设计功率
deltah_n=200; %汽轮机调节级的比焓降
%第二级组
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)
%汽轮机轴端功率
G_gb=D_0/3.6-deltaG_gbt
miu_gb=input('查图11-11')
A_gb=G_gb/(miu_gb*c_1t1/v_11)*10000
l_gb=A_gb/(e*pi*d_m*10*sin(aerfa_11/180*pi))
deltah_gb1=c_1t1^2/2000*(1-fai_gb^2)
epidmsinbeta2180pi10lb1lb11pausedeltahbkthew2t220001faib2h2h1deltahbkthev2input动叶出口比体积v2导叶热力计算deltahgboumgbdeltahth1t1h2deltahgbp11input导叶出口压力mpav11input导叶出口比体积c1t1sqrt2000deltahgbc22faigbinput查图118faigbc11faigbc1t1deltaaerfainput输入deltaaerfa510aerfa11aerfa2deltaaerfalgb1lb1deltadeltaggbte06pidmlgb110000001sqrt2000oumgbdeltahtv11ggbd036deltaggbtmiugbinput查图1111agbggbmiugbc1t1v1110000lgbagbepidm10sinaerfa11180pideltahgb1c1t1220001faigb2导叶出口实际参数h1shijih1t1deltahgb1v1shijiinput由实际出口焓查焓熵图实际出口比体积第二列动叶热力计算deltahb1oumb1deltahth2t1h1shijideltahb1p21input动叶出口压力v2t1input动叶出口比体积w11sqrtc112u22c11ucosaerfa11180pibeta11asinc11sinaerfa11180piw11pi180w2t1sqrt2000deltahb1w112faib1input查图118w21faib1w2t1deltabeta2input718beta21beta11deltabeta2c21sqrtw212u22w21ucosbeta21180piaerfa21asinw21sinbeta21180pic21180pideltahb11w2t1220001faib12动叶损失deltahc21c2122000h21shijih2t1deltahc21v21shijiinput动叶出口比体积lb21lgbdeltaoumbr111oumb1dmdmlb211000oumbt111oumbr11lb2110001lb211000deltagbt1e06pidmlb2110000001sqrt2000oumbt1deltahtv2t1miub1input查图1
%导叶出口实际参数
h_1_shiji=h_1t1+deltah_gb1
v_1_shiji=input('由实际出口焓查焓熵图—实际出口比体积')