蒸汽发生器课程设计说明书
学院:核科学与技术学院
09级
目录
第一章绪论...................................................... 第二章蒸汽发生器的设计与计算......................................
2.1 根据热平衡确定换热量........................................
2.2 管径的选取以及传热管数目的确定 ..............................
2.3 换热面积的计算..............................................
2.4 管束结构的计算..............................................
2.5 强度计算....................................................
2.6 主要管道内径的计算..........................................
2.7 一回路水阻力计算............................................
2.8 二回路水循环阻力计算........................................
2.9 运动压头计算................................................
2.10 循环倍率的确定.............................................. 第三章结论与评价..................................................
附录1 蒸汽发生器热力计算表
附录2 蒸汽发生器水力计算表
附录3 蒸汽发生器强度计算表
第一章绪论
蒸汽发生器是产生汽轮机所需蒸汽的换热设备。在核反应堆中,核裂变产生的能量由冷却剂带出,通过蒸汽发生器将热量传递给二回路工质,使其产生具有一定温度一定压力和一定干度的蒸汽。此蒸汽再进入汽轮机中做功,转换为电能或机械能。在这个能量转换过程中,蒸汽发生器既是一回路的设备,又是二回路的设备,所以被称为一、二回路的枢纽。蒸汽发生器作为一回路主设备,主要功能有:
1、将一回路冷却剂的热量通过传热管传递给二回路给水,加热给水至沸腾,经过汽水分离后产生驱动汽轮机的干饱和蒸汽;
2、作为一回路压力边界,承受一回路压力,并与一回路其他压力边界共同构成防止放射性裂变产物溢出的第三道安全屏障;
3、在预期运行事件、设计基准事故工况以及过度工况下保证反应堆装置的可靠运行。
实际运行经验表明,蒸汽发生器能否安全、可靠的运行,对整个核动力装置的经济性和安全可靠性有着十分重要的影响。据压水堆核电厂事故统计显示,蒸汽发生器在核电厂事故中居重要地位。一些蒸汽发生器的可靠性是比较低的,它对核电厂的安全性、可靠性和经济效益有重大影响。因此,各国都把研究与改进蒸汽发生器当做完善压水堆核电厂技术的重要环节,并制定了庞大的科研计划,主要包括蒸汽发生器热工水力分析;腐蚀理论与传热管材料的研制;无损探伤技术;振动、磨损、疲劳研究;改进结构设计,减少腐蚀化学物的浓缩;改进水质控制等。
第二章蒸汽发生器的设计与计算
蒸汽发生器的设计计算包括热力计算、水动力计算、强度计算、结构设计等。热力学计算主要通过传热学知识计算传热管传热面积,然后与管束结构设计相结合可以得到传热管长度、管束直径等结构参数。强度计算通过压力校验等用于选取蒸汽发生器结构材料、确定结构尺寸等。蒸汽发生器中要通过强度计算得出的参数有传热管、上下筒体、球形下封头管板等的壁厚。在强度计算得出参数后要留取一定余量,以满足变工况下出现超压情况的需求。水动力计算是最后一步,因为只有结构尺寸确定,运动状态已知的情况下水力情况才得以确定。水力计算包括一回路水阻力计算、二回路水循环阻力计算、运动压头计算等。其中一回路水阻力计算相对简单主要包括单相水U型管管内摩擦阻力和局部阻力两项;二回路水循环阻力计算及运动压头计算比较复杂。因为蒸汽发生器内二回路侧流体的水力特性取决于流体工质的性质和状态、流道的结构和几何形态,以及工质的流动形式。而且计算过程误差较大,往往需要在试验中进一步修正。二回路水循环阻力包括下降空间阻力、上升空间阻力、汽水分离器阻力等。而上升空间阻力又包括摩擦阻力、局部阻力、弯管区阻力、加速阻力、流量分配孔阻力五项。设计中常用图解法来确定循环倍率,即先假设几个不同的循环倍率分别计算其运动压头和总阻力,在直角坐标系作出相应曲线,两根曲线交点即为稳定工况的循环倍率值。循环倍率值一般取2~5为宜,其值过小会导致传热恶化,腐蚀加剧等;而过大则会增大汽水分离器负荷,使蒸汽干度降低,危机汽轮机安全。
计算过程中水力计算是在结构选型和热力计算之后进行,但是结构设计和热力计算又需要水力计算数据,因此三者往往要反复交替进行,以使设计逐步完善。
2.1根据热平衡确定换热量
一回路进口焓值:'1i =1394.21 kj/kg (15.0MPa ,310℃); 一回路出口焓值:''1i =1284.45 kj/kg (15.0MPa ,290℃); 二回路给水焓:
f
i =944.38 kj/kg (5.0MPa ,220℃);
二回路饱和水温:s t =264 ℃ (5.0MPa ); 二回路饱和水焓:s i =1154.50kj/kg (5.0MPa ); 二回路饱和蒸汽焓值:''s i =2793.6 kj/kg (5.0MPa );
二回路汽化潜热: r=''s i -s i =2793.6-1154.2=1639.73kj/kg (5.0MPa ); 排污量:0.01126 1.26/d s D C D kg s ==?=; 干度:x=0.99:
换热量:()()231278.98=++-=d s f Q Drx D D i i ; 一回路水流量: 1'''112128.40()
η=
=-Q
kg G s i i ; η=0.99为蒸汽发生器的热效率。
2.2 管径的选取以及传热管数目的确定
选取传热管的外径为:d 0=22mm ;
节距选取为:t=1.4*d 0=30.8mm ; 则最小节圆半径为:min 261.6R t mm == ; 负公差修正系数:?=1.102; 弯曲减薄系数:min
1 1.0824o
R d R ?=+
=;
一回路侧设计压:221,14
1.2515*1000000/9.8*10/191.3265/设-===P P kg cm kg cm
传热管壁厚:,111,1
1.34200[]0.8设设?φσ=
=+o
R P d S mm P ,
此处取1 1.5=S mm ;
传热管内径:1219=-=i o d d S mm ;
单管流通面积:22211
1928444
ππ=
=?=i a d mm 选取一回路侧水流速度为:1 5.0/u m s = 一回路水平均比容:310.001378/v m kg =(15.0MPa ,300℃,饱和水)
则一回路侧流通面积为:211
1
0.5866m ==G v A u U 型管数目为:111/2128.40*0.001378/52069284
====A G v u n a a 根据排管最后确定U 型管数目为: n=2086
2.3 换热面积的计算
一回路侧水导热系数:10.5614/(℃)λ=?W m (15.0MPa ,300℃,饱和水) 一回路侧水普朗特数:Pr 0.8616=f (15.0MPa ,300℃,饱和水)
一回路侧水动力粘度:-5
18.83310kg/m s η=?g
(15.0MPa ,300℃,饱和水) 一回路侧雷诺数:5111
Re 7.810η=
=?i
f u d v 一回路侧换热系数:0.80.41
10.023
Re Pr 33128.72λα==f f i
d
传热管导热系数:17.4/(℃)w W m λ=? (给定)
传热管壁热阻:52ln 9.2710/2℃λ-=
=??o o w w i
d d
R m W d 污垢热阻:5
2
2.610/℃-=??f R m W (按I-600选定)
传热温差: 'max 1''
min 1max min ln max min
31026446℃
29026426℃
4626
35℃
46
ln ln
26s s t t t t t t t t t t t ?=-=-=?=-=-=?-?-?=
==??
1
21
211q 5258假设一个值,进行迭代可以求得试取αα-??
=+++=? ???o w f i d k R R W m K
d 52
ln 0.150.722 1.84610/0.55727357.7/()
则此时的热负荷为:℃α=??=?==?s q k t W m p q W m
2
ln
1252.80所需传热面积:=
=??Q
F m k t 传热裕度系数:C=1.1
21378.08设设计传热面积:=?=F C F m
2.4 管束结构的计算
19939设
总传热管总长:π=
=o
F L m d 传热管排列方式:按正方形排列
U 40.1232节
最小型管节圆直径:==D
t m
实际布管数:n=2086根
D 2.34tb 管束直径:=m max 1.17弯管段高为:=D mm
平均直径:0.5(0.1232 1.17)0.6466m =?+=D 弯管总长:0.54036.87弯π=?=L D n m
直管总长:199394036.8715902.13直总弯=-=-=L L L m
管束直段高: 3.82直
直=
=L H m n
管束弯段高:max 0.5 1.17弯==H D m
管束总高: 4.98tb 直弯=+=H H H m
2.5 强度计算
12 2.38管束()衬筒内径:δ=+=wi t D D m
其中δt 是装配间隙,约10~20mm ,取20mm 。
2 2.404衬筒外径:δ=+=wo wi D D m ,其中δ是衬筒壁厚,给定为12mm 。 ,22 2.58下
()下筒体内径:=+=wo i D
D B m
其中B 为下降流道宽度,取为88mm 。
22
1.2563.7755/设,二回路侧设计压力:==s P P kg cm 下筒体许用应力:[σ]=18kg/mm2
下筒体计算壁厚:,2'
2
46.72200[] 1.2下设,设,下筒体计算壁厚:σ?=
=-i P D S mm P
47设计壁厚为:=S mm
,,22675下下下筒体外径为:=+=o i D D S mm
(3)上筒体许用应力:[σ]=18kg/mm2
,上上筒体内径:3200i D mm = (给定) 2,'
2
57.92200[] 1.2设,上设,上筒体计算壁厚:σ?=
=-i P D S mm P
设计壁厚取为:58S mm =
,上,上则上筒体外径为:2 3.316o i D D S m =+=
(4)球形下封头许用应力:2
/5.14][mm kg =σ
2675球形封头外径:==o s D D mm
'11
83.82400[] 1.6设,设,计算壁厚:σ?=
=+o
P D S mm P
设计壁厚取为:84=S mm
(5)管板许用应力:2
/1800][mm kg =σ
2581承压部分直径:==i D D mm
'1437.572
计算壁厚:=
=S mm
设计壁厚选取为:440=S mm 堆焊层厚度:10堆=S mm
6210.44总
总()传热管实际平均长度:=
+=l l S m n
max max 2215.86m 直最长管子长:π=++=l R H S max min 228.89直最短管子长:π=++=l R H S m
2.6 主要管道内径的计算
主管道计算流速:'
1010/u m s = (选取,
8~12m/s)
'10.61主管道计算内=
=i d m 10.62主管道设计内径:=i d m
11
10214 9.72/主管道设计流速:π=
=i
G v u m s d '''32(1)0.039064/新蒸汽比容:=+-=v v x v x m kg
蒸汽管计算流速:'
235/=u m s (选取,30~40m/s )
'
20.4231蒸汽管计算内径:=
=i d m
蒸汽管设计内径:20.43=i d m
22
222433.89/蒸汽管设计流速:π=
=i
G v u m s d 二回路给水比容:330.001187/v m s =
给水管计算流速:'
34/=u m s (选取,2~5m/s )
'30.219给水管计算内径:=
=i d m 给水管设计内径:30.22=i d m
33
32
34 3.97/给水管设计流速:π=
=i
G v u m s d 2.7 一回路水阻力计算
(1)U 型管内摩擦阻力计算
考虑堵管后流速'
111.05 1.05 5.00 5.25/u u m s ==?=
'111
Re 819327.34一回路侧水雷诺数:η==i u d v
摩擦阻力系数:0.25
0.3164Re
0.010826λ-==
平均壁温:0.5()0.5(300264)282℃a s t t t =+=?+=
在此壁温下的动力粘度:'5
19.54410η-=??Pa s
温度修正系数:0.14
1'
1(
)0.989温度修正系数:ηφη== 2
'
11
54139.772摩擦阻力:λφ??==?f i l u P Pa d v
(2)局部阻力计算
下封头内径: 2.507=l D m
222.478
水室截面积:π
=
=c l F D m
进口管内径:10.62=i d m (与主管道相同) 进口管截面积:222110.6200.3024
4
π
π
=
=
?=i A d m
突扩阻力系数:2
11(1)0.7703ξ=-
=c
A F 一回路侧水入口处比容:310.001421m /kg i v = (15.0MPa ,310℃)
1112149.72/入口管内流速:π=
=i
i i
G v u m s d 2111
125585.252从入口管至水室阻力:ξ?==i
i
u P Pa v 水室转弯45°阻力系数:20.9ξ= (查表,线性插值)
212214529892.472水室转弯阻力:ξ??==i
i
u P Pa v
考虑堵管后截面积:'
20.5591.05
=
=A
A m 系数:'
0.23=c
A F
查图得传热管入口阻力系数:30.4ξ=
传热管入口阻力:2
'
133
13878.632ξ?==i
u P Pa v U 型管转180°阻力系数:4
0.5ξ
=
U 型管转180°阻力:2'1
44
1
4999.522ξ?==u
P Pa v 传热管出口阻力系数:50.53ξ=
出口处水比容(15.0MPa ,290℃):320.001341/
v m kg =
传热管出口阻力:2
'155
2
5444.302ξ?==u P Pa v 出口管内流速:12
21
9.46/=
=G v u m s A 水室转弯45°阻力系数:60.9ξ= (查表,线性插值)
水室转弯阻力:22
66
2
30004.972ξ?==u P Pa v 出口管突缩阻力系数:70.46ξ=
出口管突缩阻力:22
77
2
15335.882ξ?==u P Pa v 总阻力:7
1
169280.78=?=?+
?=∑f i
i P P P Pa
设计阻力: 1.1186208.86设?=?=P P Pa
2.8 二回路水循环阻力计算
(1)下降空间阻力
下降空间水比容:30.001286/d v m kg =
入口阻力系数:in 1ξ= 出口阻力系数:out 1ξ=
定位装置阻力系数:f 1ξ=
下降空间高度:
0 5.78=H m
衬筒外径:0 5.78=H m 下筒体内径:, 2.58下
=i D
m
下降空间当量直径:,0.176下=-=wo i De D D m 绝对粗糙度:0.15mm ?= 摩擦系数:22
176(1.742lg )(1.742lg )0.012220.15
λ--=+=+?=??d De 下降空间截面积:222,()0.6894
下π
=
-=d wo i F D D m
下降空间水流速:1260.001286
[3,4,5][0.71,0.94,1.18]/0.68
?=
=?=R d d d C Dv u m s F 下降空间阻力:20()[700.56,1038.51,1622.68]2
ρλξ?=+∑=d d
d H u P Pa D
e (2)上升空间阻力
a) 摩擦阻力
饱和水比容:'
30.001286/v
m kg =
饱和蒸汽比容:''
30.039439/v m kg = 饱和蒸汽密度:''
3''
1
25.36/kg m v
ρ
==
衬筒内径: 2.3808=wi D m 传热管外径:0.022=o d m 支撑板定位拉杆数量:'
10=n
上升空间流通面积:2'22[(2)] 2.864
π
=
-+=u wi o F D n n d m
上升空间当量直径: 040.0386[(2')]
π=
=++u
wi F de D n n d
循环速度:'01260.001286
[3,4,5][0.170,0.227,0.283]/2.81?==?=R u C Dv u m s F
出口水相折算速度:'
'02
(1)[0.1133,0.1699,0.2265]/-==R u
C Dv u m s F
液相平均折算速度:'
'0021()[0.1416,0.1982,0.2548/2
=
+=o u u u m s 出口汽相折算速度:''
''02
1.737/==u
Dv u m s F
汽相平均折算速度:''
''
00210.87/2
=
=u u m s 液相动力粘度:2
0.000100003/=l v m s
汽相动力粘度:52
1.80310/-=?g v m s
液相雷诺数:'
Re [42491.64,59488.29,76484.95]==o lo l u de
v 汽相雷诺数:''
Re 47130.34==o go g
u de
v 判别流型:都为紊流。 管束直段高: 3.81直=H m
液相摩阻系数:0.2
0.184Re 0.0220λ-==lo lo
汽相摩阻系数:0.25
0.3164Re 0.0214λ-==go go
按折算速度计算的液相摩擦阻力:'2
'
()[16.95,30.54,47.41]2直ρλ?==o f lo lo H u P Pa de
按折算速度计算的汽相摩擦阻力:''2
'1
() 6.7532直ρλ?==o
f go go H u P Pa de
参量X
:[1.58,2.13,2.65]=
=X
2l Φ:22
2011[14.02,10.62,8.69]Φ=++=l
X X
2g
Φ:22
120[35.20,48.06,61.02]Φ=++=g X X
液相摩擦阻力:2
()()[237.68,324.53,412.04]?=Φ?=f l l f lo P P Pa 汽相摩擦阻力:2
()()[237.68,324.53,412.04]?=Φ?=f g g f go P P Pa 摩擦阻力:1
[()()][237.68,324.53,412.04]2
?=?+?=f f l f g P P P Pa b )局部阻力:
支撑板数目:N=6
上升流道单元面积:2
22568.54
π
=-=u o A t d mm
支撑板单元开孔面积:2222
(1310)113.13
π=-=u a mm 面积比:0.2
局部阻力系数:50ξ=l
按折算速度计算的液相局部阻力: '2
'
()[2337.36,4581.22,7573.04]2ρξ?==o l lo l
u P N Pa 按折算速度计算的汽相局部阻力:''2
''
1()955.6332
ρξ?=
=o l go l u P N Pa 参量X
:[1.56,2.19,2.82]=
=X
二回路工作压力: 5.0s P MPa = 临界压力:22.12c P MPa =
参数1
(0.190.92
) 2.513s R c
P Z P -=+= 参数K :1
2.911R R
K Z Z =
+
=
2221
1[3.27,2.54,2.16]:ΦΦ=+
+=l l K X X
222
1[8.00,12.17,17.12]:ΦΦ=++=g g KX X
液相局部阻力:2()()[7643.30,11627.28,16359.22]?=Φ?=l l l l lo P P Pa
汽相局部阻力:2
()()[7643.30,11627.28,16359.22]?=Φ?=g l g l go P P Pa 局部阻力:[7643.30,11627.28,16359.22]?=l P Pa
c )弯管区阻力
管束弯头最大节圆直径: 2.34=b d mm
弯管区重心至圆心距离:0.21220.4967==s b y d 计算冲刷排数:'116=
-?s
y t
N
系数:12 1.4o
t
x x d ==
= 系数:1
1.13
0.43 1.24n x =+
=
液相雷诺数:'02Re [33993.31,50989.96,67986.62]==lo l
u de
v
汽相雷诺数:''2Re 94260.68==o go g
u de
v
液相摩擦阻力:0.15
210.084[0.044]Re [0.328,0.308,0.295](1)
λ-=+
=-lo lo n
x x 汽相摩擦阻力:0.15
1.5715Re 0.281λ-==go go
液相阻力:'2'
'
2()[26.15,55.37,94.28]2
ρλ?==o b lo lo
u P N Pa 汽相阻力: ''2'''
2()57.352
ρλ?==o b go go
u P N Pa 参量X
:[0.67,0.98,1.28]=
=X
22o o 2
201:1[7.5,5.00,3.88]ΦΦ=+
+=l l X X 222
:120[3.42,4.83
6.38],ΦΦ=++=g g X X 液相阻力:2()()[196.219,276.733,365.645]?=Φ?=b l l b lo P P Pa
汽相阻力:2
()()[196.219,276.733,365.645]?=Φ?=b g g b go P P Pa 弯管区阻力: 1
[()()][196.219,276.733,365.645]2
?=?+?=b b l b g P P P Pa d )加速阻力
管束出口质量含汽率:21
[0.333,0.25,0.20]R
x C =
=
管束出口体积含汽率:
2
''
22
2''
'
[0.939,0.911,0.885](1)
x x x ρβρρ=
=-+
系数:60.8330.05ln(
)0.91310
=+=s
p C 管束出口截面含汽率:22[0.857,0.832,0.808]φβ==C 质量流速:'2
0[132.07,176.095,220.119]/ρ?
==G u kg m s
加速阻力:2
22
22''''
22(1)1
[][136.72,185.61,240.11](1)ρφρφρ
?-?=+-=-a x x P G Pa e )流量分配孔板阻力:
阻力系数:
8ξ=h
孔板局部阻力:22
[89.75,159.56,249.32]ρξ?==o u h h
P Pa
上升空间阻力:[8303.67,12573.71,17626.34]?=?+?+?+?+?=r f l b a h P P P P P P Pa f )汽水分离器阻力:[12600,14900,17090]?=s P Pa
g )循环总阻力:[21604.23,28512.23,36339.02]总?=?+?+?=d r s P P
P P Pa 2.9 运动压头计算
(1)预热段高度计算:
二回路给水焓:944380/=f i J kg 二回路饱和水焓:1154200/=s i J kg 液面高度: 5.78=o H m
下降空间下端压力: 5.04ρ=+=low s d o p p g H MPa
low
p 压力下的饱和水焓:1157255/=sl i J kg
0.062/a -?==?-g sl s low s
i i i
J kg P p p p 预热段高度:()
[0.552,0.558,0.564]2ρπρ-?+
+??=
=?+?s f
d o d R p o d i i i
gH P C p
H m d nq i
g G p
(2)运动压头计算
蒸发段高:1 4.43=-=r tb p H H H m 管束上方区段高度:20.8=-=r o tb H H H m
蒸发段平均质量含汽率:12
1
[0.167,0.125,0.100]2x x m =
= 蒸发段平均体积含汽率:1
''
11
1''
'
[0.860,0.814,0.773](1)
x x x ρβρρ=
=-+
蒸发段平均截面含汽率:11[0.78,0.74,0.71]φβ==C 管束上方区段平均截面含汽率:22[0.857,0.832,0.808]φφ==
蒸发段运动压头: '''
111()[25640.40,24278.86,23054.63]ρρφ=-=m r P g H Pa
管束上方区段压头:'''222()[5055.84,4905.66,4764.14]ρρφ=-=m r P g H Pa
运动压头:12[30696.24,29184.52,27818.77]=+=m m m P P P Pa
2.10 循环倍率的确定
将循环倍率、循环运动压头和循环运动压头对应数值列出表并做曲线,得以下结果:
循环倍率 3 4 5
总阻力21604.23358 28512.22731 36339.02244
运动压头30696.24145 29184.52 27818.77
由上图可近似确定循环倍率为3.8,满足2到5之间的要求。
第三章结论与评价
本设计说明书主要对蒸汽发生器的设计计算进行说明解释。基本完成了课程设计要求的对蒸汽发生器换热、水力和强度等项目的计算任务。
通览本次设计,本人对蒸汽发生器的结构、原理有了进一步理解。同时也对蒸汽发生器“设计”的大体步骤与过程初步了解。应该说,此次设计是十分不成熟的,特别是在管束直径和弯管计算以及传热管排布等环节做了较大简化处理,另外,由于对机械加工、材料等方面知识了解不够,结构设计也存在诸多不不严谨设计,这一方面是因为设计者的经验不足,二是时间所限,计算结果没有得到进一步仿真验证,故设计结果实际应用价值不大。