蒸汽发生器课程设计matlab文本.pdfVIP

formatlong

%热力计算

D=126;%蒸汽产量

x=0.99;%蒸汽干度

nr=0.99;%蒸汽发生器热效率

p1=15.0;%一回路侧额定工作压力

ps1=1.25*p1*100/9.8;%一回路侧设计压力

t11=310;%一回路侧冷却剂入口温度

t111=290;%一回路侧冷却剂出口温度

ta=(t11+t111)/2;%一回路侧冷却剂平均温度

tf=220;%二回路侧给水温度

ps=5;%二回路侧额定工作压力

ps2=1.25*ps*100/9.8;%二回路侧设计压力

lw=17.4;%传热管壁导热系数

o1=18;%传热管壁许用应力

o2=18;%下筒体许用应力

o3=18;%上筒体许用应力

o4=14.5;%球形下封头许用应力

o5=1800;%管板许用应力

d0=22;%传热管外径

t=1.4*d0;%传热管最小节距,取1.4

h11=1394.21;%一回路进口焓值

h111=1284.45;%一回路出口焓值

hf=944.38;%二回路给水焓

ts=263.94;%二回路饱和水温

hs=1154.50;%二回路饱和水焓

hs11=2794.23;%二回路饱和蒸汽焓值

r=hs11-hs;%二回路汽化潜热

Cs=0.01;%排放系数

Dd=Cs*D;%排污量

Q=D*r*x+(D+Dd)*(hs-hf);%换热量

G1=Q/(nr*(h11-h111));%一回路水流量

%管径的选取以及传热管数目的确定

Rmin=2*t;%最小节圆半径

O=1.102;%负公差修正系数

OR=1+d0/(4*Rmin);%弯曲减薄系数

s111=ps1*d0/(200*o1+0.8*ps1);%直管计算壁厚

s11=s111*O*OR;%直管壁厚

s1=1.4;%取大于s11的值

di=d0-2*s1;%传热管内径

a=0.25*pi*di^2;%单管流通面积

u1=5.5;%回路侧水流速度

[~,v1,~,~]=PT(p1,ta);%一回路水平均比容

A=G1*v1/u1;%回路侧流通面积

n=A/(a*1e-6);%U型管数目

%换热面积的计算

l1=0.561436;%一回路侧水导热系数

Prf=0.861567509;%一回路侧水普朗特数

N1=8.833336e-005;%一回路侧水动力粘度

Ref=u1*di/(N1*v1*1000);%一回路侧雷诺数

a1=0.023*l1*1000/di*Ref^0.8*Prf^0.4;%一回路侧换热系数

Rw=d0/(1000*2*lw)*log(d0/di);%传热管壁热阻

Rf=2.6e-5;%污垢热阻

dtmax=t11-ts;

dtmin=t111-ts;

dtln=(dtmax-dtmin)/log(dtmax/dtmin);%传热温差

k0=500;%假设传热系数k为500

fori=0:10000%对k进行迭代

q=k0*dtln;%热负荷

a2=0.557*(ps*1000000)^0.15*q^0.7;%二回路侧换热系数

kt=1/(d0/(di*a1)+Rw+1/a2+Rf);%传热系数

ifabs(kt-k0)1%若计算值与假设值相差小于1

k=kt;%输出传热系数

break

else

k0=kt;%否则另计算值为假设值迭代

end

end

F=Q*1000/q;%计算传热面积

C=1.08;%传热裕度系数

Fs=C*F;%设计传热面积

%管束结构的计算

Lz=Fs*1000/(pi*d0);%传热管总长

Dj=4*t/1000;%最小U型管节圆直径

nx=5;%最大节圆半径上所能排列传热管根数假设值

fori=0:100

nt=0;%1/4最大节圆内的传热管根数

forj=2:nx%j表示该管束的纵坐标

form=1:(nx-1)%m表示该管束的横坐标

ifj^2+m^2=(nx)^2%该管束在圆内

nt=nt+1;

end

end