汽轮机转子的传热分析和热应力计算_唐军

( 1. Automation Department, School of Electr onics, Informati on and Electrical Engineering, Shanghai Jiaotong University, Shanghai 200030, China;
2. Shanghai Electrical Power Generation Equipment Co., Ltd. Shanghai Turbi ne Plant, Shanghai 200240, China)
计算的方程个数 n。子表面的轴向 和切向热应力 S 的大小正比于转子表面温度 T S
与体积平均温度T 的差值, 即:
S=
-
E 1-
v(T S-
T)
( 8)
其中, E 为弹性模 量, 为热膨胀 系数, v 为泊松
比。
在上节中推导了各体积单元的温度计算式, 由此不难得到体积平均温度T 。
=2
rS -
r 2
l, l 为选取的某段转子长度, 代入
( 5) 式有:
t
c ( r rS - r 2/ 4) 2 hrS + k( r S / r- 0. 5)
( 6)
用( 6) 可以计算出满足要求的 r 和 t 。
n= ( r S - rB ) / r
( 7)
由( 7) 计算出所转子分层的个数即计算机要
第 39 卷 第 3 期 2010 年 9 月
热 力 透平
THERMAL TURBIN E
V ol. 39 N o. 3 Sept. 2010
汽轮机转子的传热分析和热应力计算
唐 军1, 2 , 季 佳2 , 刘 慧2
( 1. 上海交 通大学电子信息与电气学院自动化系, 上海 200030; 2. 上海电气电站设备有限公司上海汽轮机厂, 上海 200240)
的传热系数。研究表明, 正是转子表面温度与其体 积平均温度的差异导致了表面热应力。
从目前的技术条件看, 直接测量高速转子内 外表面的温差还是有困难的, 只能按照传热学的 原理, 建立数学模型, 以蒸汽的温度、压力、放热系 数及转子材料 的物理性能参数 等为计算原始 数
据, 利用计算机进行转子温度分布的实时在线跟 踪计算, 并以此为依据优化机组的运行[ 2] 。本文 详细介绍和推导了某 300M W 机型的汽轮机的温 度分布计算模型及算法实现。
; SC3 =
E 1-
v[
1+
k] = SC1 + SC2 ; T U : 转子最靠近沟槽底部的单元
点温度;
k
=
KT
r r
2 S 2 G
-
1; K T : 轴向应力集中因子。
设:
=
S S
C1 C3
=
( 1+ k ) -
1+
k 2
( 1+ k )
rS- rG Ur
,
并且
K TEE =
KT
r r
2 S
2 G
摘 要: 详细给出了某 300M W 机型的汽轮机转子温度 分布计算 方法, 推导了 应力计 算模型 并定义 了有效 温
差的概念, 由此可将汽轮机应力控制问题转化成 有效温差的控制问题, 对实际应用具有重要意义。
关键词: 转子温度; 应力控制; 汽轮机自动控制
中图分类号: T K 262
文献标识码: A
类型的叶片固定沟槽几何结构的不同, 它们相应 的 K TEE 和 值通常是不一样的。实际上, 现有的 侧进式沟槽转子具有整体轮盘, 其转子中心线到 叶片沟槽底部的距离 r G 大约等于转子表面的半 径 r S 。如果设定 rG = rS , 则式( 9) 中的 = 1, K TEE = K T , 侧进式沟槽底部的切向热应力由下式给 定:
,
则
(
8)
式
变为:
185
第3期
汽轮机转子的传热分析和热应 力计算
Z=
E 1-
vK
TEE [ T -
T S + ( - 1) T U ]
( 10)
2. 2 侧进式( Side- Ent ry Groov e) 沟槽转子
对侧进式沟槽转子, 叶轮固定槽是轴向的, 因 而关注的应力在圆周方向或切线方向。因为两种
Abstract: In this paper, the det ailed calculation of tur bine r oto r temperatur e distr ibution is g iv en, the
models for ther mal str ess are der ived, and the concept of effect ive temper ature differ ence is defined. T hus, the pr oblem of contro lling t he therma l stress of t ur bine r oto r can be conver ted t o a pro blem of contr olling the effective temper ature difference, which is o f gr eat significant fo r actual applicatio n. Key words: roto r temperatur e; stress contro l; auto matic turbine co nt rol
式转子底部的裂纹都与轴向垂直, 因此这个位置 上的轴向 应力是最主要的。此时 轴向热应力 Z 的计算方程式为:
Z = - S C1 T S - S C2 T U + S C3 T
( 9)
式中: S C1 =
E 1-
v
( 1+ k ) -
1+
k 2
r S- rG Ur
;
SC2=
E 1-
v
1+
k 2
rS - r G Ur
1 转子温度计算模型
转子的示意图见图 1, 按 径向将转 子分成 n 个部分( 其半径差为 r ) , 转子表面半径为 r S, 轴 孔半径为 rB , SS 为表面面积, T S 为表面温度, SB
图 1 转子径向分割示意图
收稿日期: 2009- 11- 06 修改日期: 2010- 02- 10 作者简介: 唐军( 1975- ) , 男 , 2000 年毕业于燕山大学机械学院, 工程师, 长期从事大型汽轮机调节保安系统设计开发工作。目前就读上
3 结论
本文详细给出了某 300M W 机型的汽轮机转 子温度分布计算模型, 针对两种类型的热应力计 算方法, 即 T 型根槽式( T- Root Groov e) 转子和 侧进式( Side- Ent ry Groo ve) 沟槽转子, 分别推导 了计算其相应应力的计算公式。从这些计算公式 可知, 转子表面应力的大小是与相应的有效温差 成正比的, 在一定的转子材料情况下, 对应力的控 制问题可以转变为对有效温差的控制。
考虑到实际汽轮机转子的温度分布计算十分 复杂, 影响因素很多, 如果采用多维计算模型, 势 必导致计算量巨大[ 3] , 难以满足实时在线的要求, 因此合理地简化模型是必要的。
目前对实际运转的汽轮机转子, 采用直接的 应力检测还有比较大的困难[ 4] , 因此通过某些温 度点的检测以及模型的计算, 控制其有效温差, 同 样可以达到对应力控制的要求, 从而实现对汽轮 机从启动、变负荷、停机等全过程的自动控制, 并 可在此基础上进行优化。
1-
hS VS
S
c
+
S1 k VS c
r
t0
( 4)
适当整理( 4) 式, 可得:
t
1/
hS VS
S
c
+
S1 k VS c
r
( 5)
其中 SS =
2r s
l, V S =
(
r
2 S
-
r
2 1
)
l, VS=
(
r
2 S
-
r21 )
l, r1= rs-
r 2
,
V
S
=
( rS
r-
r2 / 4) l , S 1 = 2r1 l
CF =
( 3+ v) 4
g w2
r
2 B
+
( 1( 3+
v) v)
r 2S
+
(
r
2 S
+
r2B )
rs
r
2 S
-
r 2B
( 14)
将离心力导致的应力和热应力相加, 即得到
转子沟槽底部总的切向应力 T 为:
T= KS- E
E 1-
v(T-
T S )+
CF
NR NR
2
( 15)
让=
CF /
E 1-
vN
T=
E 1-
vK S- E(T -
T S)
( 11)
式中, K S - E = K T = 应力集中因子。
对空心旋转圆柱转子, 由离心力导致的切向
应力 T1 可表示为:
T1=
(3+ v) 8
g w~ 2
r
2 S
+
r
2 B
+
r2S r r2
2 B
-
( 1+ ( 3+
3v ) v)
r
2
( 12) 式中, 为转子材料的比重, g 为重力加速度, w~ 为
转子旋转速度。