25MW汽轮机课程设计计算书.doc

汽轮机课程设计
汽轮机参数：
容量：25MW
蒸汽初参数：压力:3.43Mpa 温度:435℃
排汽参数：冷却水温20℃背压：0.005～0.006Mpa (取0.005 Mpa)
前轴封漏汽与轴封加热器耗汽量为0.007D○,轴封加热器焓升 21KJ/Kg
加热器效率ηjr=0.98
设计功率:Pr=25MW
最大功率P=25*(0.2～0.3)

1. 近拟热力过程图
在焓熵图上选取进口参数P0=3.43MPa,t0=435℃,可得h0=3304kJ/Kg.设进汽机构的节流损失△P0=0.04P0，可得调节级压力=3.3 MPa,并确定调节级前蒸汽状态点1（3.3 MPa, 435℃）过1点作等比熵线向下交PZ线于2点，查得h2t=2128KJ/Kg,整机理想比焓降（△htmac）’=h0-h2t=3304-2128=1176 KJ/Kg.选取汽轮机的内效率η=0.85，有效比焓降△himac=（△htmac）’*ηri=999.6 KJ/Kg,排气比焓和hz=2304kj/kg.在焓熵图上得排汽点Z，用直线连接1，Z，去两点的中点沿等压线下移21-25Kj/Kg，用光滑曲线连接1，3两点，得热力过程曲线的近似曲线见图1，
T P0 P0’ t0
H0 1 2

△Hi
△Ht

Pc0’
Pc

S
图1
选取给水温度T=160℃ 回热级数:5
内效率η=0.85
主汽门和调节阀中节流损失△P0=(0.03～0.05)PO
排汽管中压力损失 △PC=(0.02～0.06)PC
回热抽汽管中的压力损失 △PE=(0.04～0.08)PE
2.汽轮机进汽量D○
ηm=0.99 ηg=0.97 m=1.15 △D=0.03DO
D0=/ himacηmηg*m+△D=3.6*20000*1.15/(93*0.99*0.97)
+0.03△D
=107.19 t/h
2. 抽汽压力确定
采用大气式除氧器 压力为0.118 MPA 饱和温度为104.3℃

CY

轴封加热器
1# 2# 3# 4#

加热
器号
抽汽
压力Pe
抽汽比焓he
抽汽管压损

加热器工作压力Pe’
饱和水温度te’
饱和水比焓he’
出口端差
给水出口温度tw2
给水出口比焓hw2
H1
0.9296
3024
8
0.8552
172.5
730.17
7
165.5
697.4
H2
0．48
2888
8
0．4416
147
619.27
7
140
592.4
Hd
0.1917
2764
8
0.1764
115.8
437.77
0
115.8
486.7
H3
0.09553
2644
8
0.08789
95.9
402.2
5
92.9
381.5
H4
0．03643
2492
8
0.03354
71.9
300.9
5
66.7
276.2

3. 回热抽汽流量的计算
(1) H1高加
给水量 △De=0.5 △DL1=0.77 △DC=1
Dfw=D0-△DC+△DL1+△Dej
=107.19-1+0.77+0.5=107.46 t/h
抽汽量
△De1(he1-he1’) ηjr= Dfw(hW2-hw1)
（t/h）
（2）H2高加
抽汽量
（t/h）
H1疏水流入H2放热
(t/h)
考虑前轴封漏汽
(t/h)
(t/h)
(3) Hd除氧器

(4)H3低加
(5)H4低加
回热系统的校验

4. 流经各级组蒸汽量及其内功率
调节级

第一级组

第二级组

第三级组

第四级组

第五级组

第六级组

整机内功率

5. 计算汽机装置的热经济性
机械损失:
汽机轴端功率:
发电机功率:
内功率大于25000KW,合格
汽耗率:
不抽汽估计汽耗率:

汽轮机装置的热耗率

绝对电效率
25MW凝汽式汽轮机热平衡计算数据
基本数据
汽机初压
P0
MPa
3.43
汽机初温
T0
℃
435
汽机初比焓
H0
Kj/kg
3304
工作转速
n
r/min
3000
冷却水温度
Tc1
℃
20
射汽抽气器耗汽量
△ej
t/h
0.5
射汽抽气器比焓降
△hej
Kj/kg
21
汽轮机总进汽量
D0
t/h
107.19
前轴封漏汽
△DT
t/h
0.75
流入凝汽器汽量
DC
t/h
70.272
汽机背压
PC’/PC
MPa
0.005/0.0052
凝汽器出口水温
tc
℃
36.1
给水泵压头
Pfp
MPa
6.4
抽汽冷却器出口水温
tej
℃
40
凝结泵压头
pcp
MPa
1.2

6. 双列速度级的热力计算
（1） 速度级的选择
选择双列速度级（195-250KJ/Kg）选择焓降为250kj/kg.故速度级的参数为：
℃

1. 喷嘴热力计算
(1) 喷嘴理想焓降

(2) 喷嘴进口状态参数

(3) 喷嘴出口状态参数
由△hn可以从H-S图上查得:

(4) 喷嘴形状的确定
前后压比:
选用渐缩型喷嘴.
(5) 喷嘴出口速度
理想速度:
速度系数
实际速度:
喷嘴出口汽流偏转角
喷嘴出口汽流方向角

(6) 轮周速度u

(7) 速度级的平均直径dm

(8) 喷嘴出口面积An

(9) 喷嘴出口高度ln

选取部分进汽度e=0.6则叶高ln=16mm>15mm
(10) 喷嘴损失

2. 第一列动叶热力计算
(1) 动叶进口汽流的相对速度
(2) 根据C1,U1作速度三角形，由余弦定理可得：

(3) 动叶出口汽流相对速度
因为 则
查图,

复速级动叶出口汽流角
取
(4) 动叶绝对速度

(5) 动叶进口状态参数
喷嘴出口实际状态点参数
动叶比焓
由H-S图查得动叶进口密度
(5)动叶进口高度 (△r △t由表1-1查得)

(6)动叶出口面积

( 由图1-11查得)
(7)动叶出口高度

(8)动叶损失

(9)动叶出口汽流状态参数
动叶出口比焓
查H-S图得:出口密度
因为则
3. 导叶热力计算
(1) 导叶中汽流的理想比焓降

(2)导叶出口汽流理想状态参数
由导叶进口状态( 第一列动叶出口状态)参数和△hgb从H-S图查得
导叶出口压力
导叶出口比焓
导叶出口密度
(3)导叶出口汽流理想速度

导叶出口实际速度

(由图1-18查取)
导叶出口汽流角

(4)导叶进口高度

(6) 导叶顶部漏汽量

(7) 导叶出口面积

(8) 导叶出口高度

(9) 导叶损失

(10) 导叶出口汽流实际状态参数
导叶出口焓
由H-S图查得导叶出口密度
4. 第二列动叶热力计算
(1) 动叶中汽流的理想比焓降

(2) 动叶出口汽流理想状态参数

由H-S图查得动叶出口压力
动叶出口密度
(3) 动叶进口相对速度

(4) 动叶出口汽流相对速度
相对理想速度:

相对实际速度:

(由图1-18查得)
动叶出口汽流相对速度角

(5) 动叶出口汽流绝对速度

(6) 动叶损失

(7) 余速损失

(8) 动叶出口汽流实际状态参数
动叶出口实际比焓
(9) 动叶进口高度

(10) 动叶顶部漏汽量

由于,
根部反动度

顶部反动度

(11) 动叶出口面积

(由图1-11查得)
(12) 动叶出口高度

5. 轮周功校核
1KG蒸汽所做的轮周功

计算符合要求
6. 轮周效率

7. 级内损失的计算
(1) 叶轮摩擦损失

(2) 叶高损失

(3) 部分进汽损失
鼓风损失

斥汽损失

(4) 导叶及动叶顶部漏汽损失

8. 级的内功率

9. 级的内效率

7. 压力级的确定及焓降的分配
1. 第一压力级的平均直径
=1.11m
2. 凝汽式汽轮机末级直径的估算

3. 平均理想焓降的计算
各级组的直径及反动度
L1
55.5mm
d1
1100mm
X1
0.421
L2
57.5mm
d2
1150mm
X2
0.431
L3
62.5mm
d3
1250mm
X3
0.441
L4
67.5mm
d4
1350mm
X4
0.4446
L5
73.5mm
d5
1470mm
X5
0.456
L6
83mm
d6
1660mm
X6
0.5

各级的理想焓降估算

级的平均理想焓降

级数目的确定

比焓降分配辅助表格
级号
平均直径
速比
计算比焓
取定比焓
级后压力
dm
Xa
△ht
△ht’
Pi’
1
1110
0.421
85.67
68
1.3
2
1130
0.43
85.10
72
0．8
3
1160
0.44
86.56
76
0.58
4
1210
0.45
89.10
80
0.4
5
1270
0.455
96.01
80
0.26
6
1340
0.46
104.58
84
0.14
7
1380
0.47
106.21
84
0.08
8
1440
0.48
110.96
104
0.03
9
1540
0.49
121.73
104
0.015
10
1600
0.5
135.84
12
0.0051
8. 回热系统抽汽压力的重新确定
加热
器号
抽汽
压力Pe
抽汽比焓he
抽汽管压损

加热器工作压力Pe’
饱和水温度te’
饱和水比焓he’
出口端差
给水出口温度tw2
给水出口比焓hw2
H1
1
3074
8
0.8552
177
749.9
7
170
723
H2
0.58
2904
8
0.5336
154
649.6
7
147
622.83
Hd
0．26
2748
17
0.2392
125.6
527.6
0
125.6
531.6
H3
0.08
2608
8
0.0736
94
393.78
5
89
372
H4
0.035
2500
8
0.032
66.13
276.75
5
61
276.75

(1) H1高加
给水量
Dfw=D0-△DC+△DL1+△Dej
=107.19-0.75+0.58+0.5
=107.52 t/h
抽汽量
△De1(he1-he1’) ηjr= Dfw(hW2-hw1)
（t/h）
(2) H2高加

(t/h)
(t/h)

(3) Hd除氧器

(4) H3低加

(5) H4低加

回热系统的校验

流经各级组流量及其内功率
调节级

第一级组

第二级组

第三级组

第四级组

第五级组

第六级组

整机内功率

装置热经济性
机械损失
汽机轴端损失
发电机功率
汽耗率

不抽汽估计汽耗率

汽机装置热耗率

绝对电效率

9.压力级第九级第十级的详细热力计算演示
1.级内的比焓降分配
(1)焓降 104kj/kg
初焓 2500 初压 0.037MP
初速 反动度
等熵滞止焓降
(2) 蒸汽在动叶的理想比焓降：

2.喷管的热力计算
⑴ 喷管前后的蒸汽参数
根据,以由h-s图得
喷管滞止压力=0.037 滞止比焓=2540.3
滞止密度=0.22 喷管前比焓2500
喷管后压力=0.017MP 理想密度 =0.125
理想比焓 =2418
⑵ 喷管截面积形状的确定
等熵指数 k=1.035+0.1=1.129
临界压比 ==0.566
喷管前后压力比 =0.016/0.035=0.457
因为0.457，所以汽流在喷管出口为超声速流动但是>0.3~0.4
故喷管应该是渐缩型超音速斜切部分达到超音速。
⑶ 临界参数的计算
临界压力 =*=0.037*0.566=0.021MP
临界焓 =2440
临界密度 =0.147
临界速度 ==358.6m/s
⑷ 喷管出口汽流速度
喷管出口汽流理想速度 ==416.25m/s
喷管出口的实际汽流速度
=0.97*416.25=403.77m/s
喷管出口汽流出口角
因为喷管出口压力，斜切部分中汽流产生膨胀，
发生偏转，所以喷管汽流出口角应为喷管出气角加上
而

⑸ 隔板漏气量计算

⑹ 流经喷管的流量
、
⑺ 喷管叶栅出口面积

⑻ 喷管出口高度

⑼ 喷管损失

⑽ 喷管出口蒸汽参数
喷管出口实际比焓值

喷管出口的实际密度值

3 动叶栅热力计算
（1） 动叶栅进口速度
圆周速度
动叶进口相对速度
的方向角

（2） 动叶出口速度
动叶进口汽流焓

动叶出口汽流理想相对速度

动叶出口汽流的实际相对速度

动叶出口相对速度汽流角
一般选取时应使动叶出口高度不低于进口高度
现取
根据ω2，u，可作出动叶出口速度三角形
动叶出口汽流绝对速度C2

=114.2m/s
动叶出口汽流的速度方向角α2：

=83.840
(3)动叶损失Δhbε=ω2t2(1-φ2)/2=4.56kj/kg
(4)动叶出口蒸汽参数
动叶出口理想比焓降h2t=h1-Δhb=2423-21.66=2401.34kj/kg
动叶出口理想密度ρ2t=0.111kg/m3
动叶出口压力P2=0.016MPa

动叶出口实际比焓
动叶出口实际密度ρ2=0.115 kg/m3
动叶出口干度x2=0.917
（5）动叶出口高度
（6）动叶定不漏气量积计算

（7）动叶出口截面积

（8）动叶出口高度
lb>lb’ 故β2选择符合要求。
（9）余速动能Δhc2

4. 轮周功及功率的计算
通过轮周的有效焓降计算轮周功

通过速度三角形计算的轮周功

Pu与Pu’基本相同，符合要求。
5.级内的各项损失的计算
（1）叶高损失

（2）叶轮摩擦损失
摩擦耗功
用能量表示摩擦损失
（3） 隔板漏汽损失

（4）顶部漏气

（5）湿气损失

（6）级的内功率

10.第十级的详细热力计算
1.级内的比焓降分配
(1)焓降 102kj/kg
初焓 2404 初压 0.016MP
初速 反动度
等熵滞止焓降
(3) 蒸汽在动叶的理想比焓降：

2.喷管的热力计算
⑴ 喷管前后的蒸汽参数
根据,以由h-s图得
喷管滞止压力=0.037MP a 滞止比焓=2540.3 kj/kg
滞止密度=0.117 喷管前比焓2404
喷管后压力=0.0055 理想密度 =0.05
理想比焓 =2328 kj/kg
⑵ 喷管截面积形状的确定
等熵指数 k=1.035+0.1=1.1235
临界压比 ==0.579
喷管前后压力比 =0.016/0.035=0.35
因为0.457，所以汽流在喷管出口为超声速流动但是>0.3~0.4
故喷管应该是渐缩型超音速斜切部分达到超音速。
⑶ 临界参数的计算
临界压力 =*=0.037*0.566=0.009246MPa
临界焓 =2352kj/kg
临界密度 =0.0769
临界速度 ==340.57m/s
⑷ 喷管出口汽流速度
喷管出口汽流理想速度 ==450.72m/s
喷管出口的实际汽流速度
=0.97*416.25=436.6m/s
喷管出口汽流出口角
因为喷管出口压力，斜切部分中汽流产生膨胀，
发生偏转，所以喷管汽流出口角应为喷管出气角加上
而

⑸ 隔板漏气量计算

⑹ 流经喷管的流量
、
⑺ 喷管叶栅出口面积

⑻ 喷管出口高度

⑼ 喷管损失

⑽ 喷管出口蒸汽参数
喷管出口实际比焓值

喷管出口的实际密度值

3 动叶栅热力计算
（4） 动叶栅进口速度
圆周速度
动叶进口相对速度
的方向角

（5） 动叶出口速度
动叶进口汽流焓

动叶出口汽流理想相对速度

动叶出口汽流的实际相对速度

动叶出口相对速度汽流角
一般选取时应使动叶出口高度不低于进口高度
现取
根据ω2，u，可作出动叶出口速度三角形
动叶出口汽流绝对速度C2
=142.3
动叶出口汽流的速度方向角α2：

=82.30
(3)动叶损Δhbε=ω2t2(1-φ2)/2=5.4kg
(4)动叶出口蒸汽参数
动叶出口理想比焓降h2t=h1-Δhb=2324-25=2299kj/kg
动叶出口理想密度ρ2t=0.045kg/m3
动叶出口压力P2=0.0045 MPa

动叶出口实际比焓
动叶出口实际密度ρ2=0.005 kg/m3
动叶出口干度x2=0883
（5）动叶出口高度
（6）动叶定不漏气量积计算

（7）动叶出口截面积

（8）动叶出口高度
lb>lb’ 故β2选择符合要求。
（9）余速动能Δhc2

5. 轮周功及功率的计算
通过轮周的有效焓降计算轮周功

通过速度三角形计算的轮周功

Pu与Pu’基本相同，符合要求。
5.级内的各项损失的计算
（1）叶高损失

（2）叶轮摩擦损失
摩擦耗功
用能量表示摩擦损失
（6） 隔板漏汽损失

（4）顶部漏气

（5）湿气损失

（6）级的内功率

级数

1

2

3

4
焓降△ht
68
72
76
80
初焓h0
6146
3074
2990
2904
初压p0
1.7
1.3
0.8
0.58
初速C0
198.99
91.89
88
92.577
反动度
0.08
0.08
0.09
0.09
等熵滞止焓降
87.79
76.22
79.872
84.258
喷嘴等熵滞止焓降
80.77
70.124
72.683
76.699
喷嘴理想出口速度
401.93
374.497
381.27
391.66
等熵出口焓值
3065.52
3003
2917.31
2827
V1t
0.19
0.31
0.27
0.6
P1
1.3
0.78
0.75
0.4
出口面积An
144.47
195
207
433
等熵滞止焓
3165.789
3003.87
2993
2908
查得P0(0)
1.7
1.4
1
0.95
压比
0.764
0.7222
0.8
0.421
隔板漏气量
0.00365
0.329
0.31
0.143
流量
107.19
102.46
102.46
98.85
有效流量
106.725
102.13
102.46
98.76
级速比
0.42
0.43
0.419
0.45
⑻ 级理想速度
419.04
390.44
399.67
410.57
⑼ 圆周速度
175.99
167.88
167.46
184.75
⑽ 平均直径动叶
1065
1069
1067
1177
⑽ 平均直径喷嘴
1063
1068
1066
1176
⑾ 喷嘴出口角a1
13
13
13
13
叶片数
46
46
46
46
估算玄长
110
110
110
110
喷嘴节距
76.45
72.9
72.76
80.27
相对节距
0.695
0.66
0.66
0.729
⑿ 喷嘴高度
18
26
28
52
⒀ 动叶高度
20
28
30
54
⒁ 喷嘴实际出口速度
389.87
363.26
369.82
379.91
⒂ 喷嘴损失
4.77
4.14
4.29
4.53
⒃ 喷嘴实际出口焓
3069.999
3008.02
2921.61
2831
⒄ 动叶进口角B1
0.4301
0.439
0.4313
0.4609
(反正切）角度B1
23.275
23.7
23.33
24.74
⒅ 动叶进口相对速度
251.61
203.28
210.06
204.16
⒆ 理想相对速度
221.94
203.21
210.06
238.44
h2t
3064.5
3002
2914
2824
查P2
1.3
1.1
0.82
0.4
V2t
0.19
0.25
0.29
0.6
[21]动叶出口面积Ab
237.65
325.53
361.15
732.41
[22] 动叶出口角B2
0.3375
0.3463
0.3593
0.3669
（反正玄）角度B2
19.72
20.265
21.05
21.53
需动叶数目
212
203
202
223
Tb
0.647
0.647
0.647
0.647
Bb
25.6
25.6
25.6
25.6
[23] 动叶出口相对速度
235.759
216.71
226.63
223.42
[24] 动叶损失
3.86
3.26
3.569
3.46
[25]动叶绝对出口角a2
1.733
2.11
1.84
3.55
(反正切）角度
60.015
64.74
61.59
74.284
[26] 动叶绝对出口速度
91.88
82.99
92.57
85.17
[27] 余速损失
4.22
3.44
4.28
3.627
[28] 叶高损失
4.99
3.01
2.9
1.677
[29] 轮周焓降
74.94
65.369
67.72
72.65
[30] 轮周效率
85.355
85.76
84.78
86.2
[31]扇形损失
0.019
0.03664
0.0441
0.1241
[32] 摩擦耗功
1.216
0.7624
0.649
0.53109

[34] 隔伴汽封漏汽损失
1.519
0.9807
0.9575
0.49
[35] 叶顶漏汽损失
1.8989
1.08
1.046
0.5628
顶部反动度
0.0632
0.055
0.06
0.04624
当量间隙
0.64
0.64
0.64
0.64
[37] 有效焓降
65.29
59.494
62.12
69.269
[38] 焓降H2
3080.7
3014
2927
2834
[39] 效率
65.29
81.16
91.49
85.3
利用系数
0.77
0.81
0.85
0.85
[40] 做功
1944
1693
1768
1902

级数
5
6
7
8
焓降△ht
80
84
84
104
初焓h0
2811
2748
2652
2608
初压p0
0.4
0.3
0.19
0.202
初速C0
85.1786
85.0055
80.0622
92.45
反动度
0.1
0.1
0.114
0.13
等熵滞止焓降
83.6276
87.6129
87.2049
108.2736
喷嘴等熵滞止焓降
75.2649
78.8516
77.2636
94.198
喷嘴理想出口速度
387.9817
197.1188
393.0995
434.0462
等熵出口焓值
2735.735
2669.1483
2574.7363
2513.8019
V1t
1
1
1
1.6
P1
0.18
0.2
0.18
0.1
出口面积An
728.9814
693.5174
700.6148
961.3623
等熵滞止焓
2814.6277
2751.6129
2655.205
2612.2736
查得P0(0)
0.32
0.32
0.195
0.21
压比
0.5625
0.625
0.92312
0.48
隔板漏气量
0.08523
0.08724
0.08636
0.0596
流量
98.85
96.26
96.26
91.13
有效流量
98.7647
96.1727
96.1736
91.07
级速比
0.455
0.46
0.47
0.48
⑻ 级理想速度
408.9687
418.5999
417.6241
465.3463
⑼ 圆周速度
186.0807
192.5559
196.2833
223.3662
⑽ 平均直径动叶
1185
1226
1250
1423
⑽ 平均直径喷嘴
1184
1225
1249
1422
⑾ 喷嘴出口角a1
13
13
13
13
叶片数
54
54
56
58
估算玄长
110
112
114
116
喷嘴节距
68.8474
71.2314
70.0332
76.9841
相对节距
0.6258
0.6359
0.6143
0.6637
⑿ 喷嘴高度
87
80
79
96
⒀ 动叶高度
89
82
180
200
⒁ 喷嘴实际出口速度
376.3423
385.2052
381.3065
421.0248
⒂ 喷嘴损失
4.4481
4.6601
4.5662
5.5671
⒃ 喷嘴实际出口焓
2740.1832
2673.8084
2579.3026
2519.369
⒄ 动叶进口角B1
0.468721
0.47489
0.489444
0.506829
(反正切）角度B1
25.1135
25.3651
26.0792
26.88
⒅ 动叶进口相对速度
199.4722
202.2767
195.1155
209.4982
⒆ 理想相对速度
237.7282
241.7404
240.7338
268.4038
h2t
2732.1832
2665.4084
2569.7266
2505.8491
查P2
0.18
0.2
0.18
0.1
V2t
1
0.92
1.1
1.8
[21]动叶出口面积Ab
1225.0893
1079.2874
1295.8601
1800.9744
[22] 动叶出口角B2
0.3699
0.3419
0.4075
0.41129
（反正玄）角度B2
21.7118
19.9929
24.0541
24.2858
需动叶数目
225
232
237
270
Tb
0.647
0.647
0.647
0.647
Bb
25.6
25.6
25.6
25.6
[23] 动叶出口相对速度
222.7513
226.5108
225.5676
254.1784
[24] 动叶损失
4.8064
3.5656
3.536
3.7169
[25]动叶绝对出口角a2
3.9488
3.8142
9.4821
12.567
(反正切）角度
75.7893
75.309
83.9797
85.45
[26] 动叶绝对出口速度
85.0055
80.0621
92.451
104.8712
[27] 余速损失
3.6129
3.2049
4.2736
6.5095
[28] 叶高损失
0.9947
1.1427
1.1366
1.1686
[29] 轮周焓降
72.1182
76.1822
74.829
93.4905
[30] 轮周效率
86.2373
86.9531
85.8082
86.2465
[31]扇形损失
0.3302
0.2743
0.2563
0.3594
[32] 摩擦耗功
0.3297
0.4366
0.402
0.4956

[34] 隔伴汽封漏汽损失
0.2897
0.3218
0.3128
0.2849
[35] 叶顶漏汽损失
0.2633
0.3746
0.3817
0.4342
顶部反动度
0.0269
0.03548
0.05261
0.0656
当量间隙
0.64
0.64
0.64
0.64
[37] 有效焓降
69.9104
73.6321
72.3395
90.7478
[38] 焓降H2
2741.0895
2674.3679
2579.6605
2517.2522
[39] 效率
86.7841
86.7395
86.559
87.59485
利用系数
0.85
0.85
0.85
0.85
[40] 做功
1919.6245
1968.8398
1934.2768
2297.1797

参考文献

1.《工程热力学》……………………中国建筑工