第47卷第5期
四辊冷轧机2021年5月水力发电
流动及传热特性研究
张杨1，刘国峰1，张维欢1，王玲花2
(1.中国电建集团西北勘测设计研究院有限公司，陕西西安710065；
2.华北水利水电大学电力学院，河南郑州450000)
摘要：为解决推力轴承瓦温过高的问题，采用数值模拟方法分析了某水轮发电机组推力轴承内循环冷却系统油槽内润滑油的流动及传热特性#结果表明，靠近镜板处油流速度大，远离镜板位置油流速度小；镜板内缘压力最低! 沿着推力轴瓦径向向外压力逐渐增大；初始瓦温、初始油温一定时，随着水温降低，整体油温下降幅度并不大#根据研究结果，对推力轴承内循环冷却系统提出适当增大冷却器U型管直线段的长度、增大冷却器管径、提高冷却水的压力和流量、加大油冷却器铜管间距的改进建议# 关键词：水轮发电机组；推力轴承；油槽；冷却器；传热特性
Study on Flow and Heaf Transfer Characteristics in Thrust Bearing Oil Tank of Hydroelectric Generating Und
ZHANG Yang1,LII Guofeng1,ZHANG Weihuan1,WANG Linghua2
(1.PowerChina Xibei Engineering Corporation Limited,Xiin710065,Shaanxi,China;
2.SchooooOEoecteicPowee,Noeth ChinaUnieeesityoOWateeResouecesand EoecteicPowee,Zhengzhou450000,Henan,China)
Abstraci:In order O solve the problem of high bearing temperature of thrust bearings,the numerical simulation method is used toanaoyzetheOoow and heat t eansOe echa eacte eistics oO oub eicating oi oin oiotank oOtheustbeaeinginteenaociecuoation coooing system of a hydro-yenerator unit.The results show that,(a)the Oow rate of oil near the thrust runner ccllar is W/c,and the Oow rate of oil away from the thrust runner ccllar is smal l;(b)the pressure on the inner edge of the thrust runner ccllar is the oowest,and thepeessueeinceeaseseapidoyoutwaed aoongthetheustshaOtbushing;and(c)when theinitiaotempeeatueeoOshaOt bushingand the initiaooiot
empeeatuee aee constant,the oeeea o oiotempeeatuee wi o notdeceease geeatoy with the watee temperature decreasing.Acccrding O the research results,some suggestions are put fo/ard for the internal circulation ccoling system oOtheustbeaeing,such asappeopeiateoyinceeasingtheoength oOthesteaightoineoOtheU-typecoooeepipe,inceeasingthe diamete eoOthe coo oe epipe,inc eeasing the p eessu ee and Ooow oOcoooingwatee,and inceeasingthedistancebetween coooee capper tubes.
Key Word t;hydLoe oect ic geneLating unit;thLust beaLing;oi otank;coo oeL;heat t ansOeLchaLacteListics
中图分类号：TH133.3文献标识码：A文章编号：0559-9342%2021)05-0095-06
0引言
水轮发电机组的推力轴承主要承担机组转动部分的重量和轴向水推力，承担的载荷往往比较大，其工作性能的好坏，直接影响到水轮发电机能否长期、安全、可靠运行，而推力轴承的油循环冷却系统又是保证推力轴承安全可靠运行的关键。水轮发电机组在运行过程中对推力轴承温度的测量只是在局部位置布置若干测点，对轴瓦采取温度监视，但这种温度监视并不能够避免轴瓦烧损，更不能反映出推力轴承油槽中整体润滑油的流动和传热规律，机组运行过程中推力轴承瓦温过高甚至烧瓦现象时有发
生。葛洲坝二江电厂1号机组在并网送电时曾发生推力瓦温过高的情况，造成1号机开关自动断开，机组与电网解列，停机后对推力轴瓦进行检查，
""收稿日期：2020-12-30
基金项目：河南省高等学校重点科研项目资助(16B57UUU3)
作者简介：张杨%1992—)，男，河南南阳人，工程师，硕士，研究方向为水力机组特性及其控制技术•
Warr%Powr%Vo.+47No+5
水力发电2021年5月
发现瓦的径向中部造成了严重烧损，有14块瓦的温度越限口+。鉴于水轮发电机组推力轴承工作的重要性，许多专家学者对推力轴承的性能进行了专门研究。王青华等*2+采用计算流体力学的方法，对推力瓦的润滑特性、受力和变形情况进行流固耦合分析，揭示了最小油膜厚度和瓦面倾角对油膜承载力、瓦块等效应力及变形的影响；赖喜德等⑶针对水电机组镜板泵流体动力设计过程中孔的布置形式以及集油槽的断面形状、面积和出口管路直径等几何参数优化问题，通过不同几何参数组合，分析镜板泵的主要过流部件的关键尺寸、形状等对其性能的影响，并在此基础上提出一种基于整个循环系统
协同的镜板泵水力优化设计方法。然而，这些研究主要是针对推力轴承本身的承载力以及推力轴承本体的结构设计方面，缺少对推力轴承整个油循环冷却系统油槽内部润滑油的流动特性及传热特性的分析。
水轮发电机组推力轴承油槽内的油循环是一个比较复杂的问题，目前尚没有较完整的试验数据，难以准确定量地计算，通常根据经验估算确定*48+。按照经验与规范，在常规设计中轴承系统安全系数裕度较大，往往引起不必要的浪费。本文以某电站水轮发电机组推力轴承内循环冷却系统为例，运用Fluent软件对推力轴承油槽内的速度分布、压力分布、温度分布进行计算，总结了油槽内润滑油的流动及传热规律，并对推力轴承油循环冷却系统的结构优化提出了改进建议。
1模型及计算方法
1.1模型的建立及网格划分
推力轴承内循环冷却系统的特点是推力轴承和油冷却器浸于同一个油槽，油的循环主要依靠轴承的转动部件旋转使油在轴承与冷却器之间流动，进行热交换，由冷却水将轴承损耗在油中的热量带走，以保证推力轴承在热平衡状态下，油温控制在规定的温度下运行⑹。本文以某水电站水轮发电机组推力轴承内循环冷却系统为研究对象，该电站水轮发电机组推力瓦的支承方式采用弹性油箱支承，弹性油箱沿圆周方向均匀分布，共18个；推力瓦为扇形双层轴瓦，固定在弹性油箱上。油冷却器采用抽屉式
油冷却器，沿圆周均匀分布，机组额定转速392.5 ed/min#
三维实体、流场和温度场耦合分析难点在于计算很难收敛，而且需要占用大量的计算机资源，对计算机的配置要求很高，大量的时间都花费在模型的简化以及网格的划分上，工作繁重，技术难度大*7+。由于推力轴承油槽内的油流场具有周期性变化的特点，结合电站具体参数，对模型做了进一步的修改和简化，只建立模型的1/18进行分析。冷却器水管为铜管，冷却器管内通入冷却水，在分析过程中将水看作低速不可压缩流体，不考虑结构变形对流体产生的影响，只传递热量，简化后的模型见图1#采用有限体积法对计算域进行网格划分，对流动复杂的部位进行网格加密处理，最终网格划分效果见图2#
图11/18模型结构示意
图2网格划分效果示意
本文模拟的是定常流动，在计算中默认油的密度是一个常数，不随温度的变化而变化。推力轴承油循环冷却系统中油的传热方式属于对流传热，在假定油流没有相变的基础上，同时对连续方程、动量方程、能量方程进行求解。
1.2计算工况及边界条件
该电站坝前取水为推力轴承冷却用水的主水源，蜗壳取水为推力轴承冷却用水的备用水源。电站所在河流常年平均水温为18.2｛，电站推力轴承运行平均油温44.1｛，平均瓦温55.2｛。结合电站运行数据及运行规程规定，取不同工况进行数值计算，计算工况如表1所示。
本文研究的是周期性流动与换热问题，选取推力轴承油循环冷却系统的1/18进行计算分析，上、
Watrr Poorr Vol47No.
5
第47卷第5期张杨，等：水轮发电机组推力轴承油槽内流动及传热特性研究
下游边界设定为周期性边界条件。镜板转动角速度为392.5rad/min，镜板不与推力轴瓦接触的外圆柱面属于拖动边界。油面边界条件设置成对称边界SYMMETRY，其余边界设置成WALL，轴瓦四周与瓦底及油沟面为强制对流面。
工况推力轴承瓦温油温水温
表1计算工况°C
155.244.124
255.244.118.2
355.244.110
推力轴承油槽内油的流动属于低速不可压缩流动，采用压力基求解器进行求解，湍流模型选用标准k_模型，近壁区的流动计算采用标准壁面函数Standard Wai Functi理，其余参数采取默认设置#选
用coupled算法，同时求解流体运动方程和能量方程，选用混合初始化Hyb/d Initialization方式对流场进行初始化处理#在计算过程中不考虑镜板与轴瓦之间的油膜，认为镜板与推力轴瓦相接。
2计算结果及分析
2.1速度和压力的计算结果及分析
不同计算工况下机组转速均为392.5rad/min，油槽内的速度分布与压力分布基本一致，以工况1为例对1/18油槽内速度分布进行分析，三维速度流线图见图3。由图3可知，油流的速度和方向比较复杂，在推力轴瓦外径处，由镜板甩出的油流向外扩散，碰到冷却器后，经冷却器冷却、热交换，冷油下沉，热油上浮，速度矢量方向发生变化，向内侧流动，补充因镜板向外甩油造成的推力瓦内圆的缺油；向内侧流动的油碰到上层薄瓦和厚托瓦后速度方向又发生了改变，进而形成了进油侧漩涡。在瓦的出油侧，轴瓦外径处速度向下，在内径，进油侧和出油侧的速度均向上。
速度/Hrs」
H1.23990e+001
1.10213e+001
9.64367e+000
8,26600e+000
6.88833e+000
烘手机5.51067e+000
4.13300e+000
2.75533e+000
1.37767e+000
0.00000e+000
图3三维速度流线
1/18模型中间轴面的速度分布和压力分布如图4所示。由图4a可知，推力轴瓦上方区域整体速度大于推力轴瓦下方区域整体速度，由于润滑油的粘滞力作用，靠近镜板处油流速较大，沿着径向向外，油流速逐渐减小，油流经冷却器，受到冷却器的阻力，冷却器U型管外缘到圆柱形弹性油箱之间的区域速度最小。由图4b可知，中间轴面镜板内缘压力最低，沿着径向向外压力均匀增大，这是由镜板的泵效应引起的。
取［=0.9m（［为三维坐标系纵轴）与［=0.55 m处水平截面进行研究，这两个水平截面速度和压力分布如图5、6所示。［=0.9m水平截面距离镜板比较近，从整体来看，［=0.9m水平截面上的速度比［=0.55m水平截面上的速度要大。由图5a可
速度/nrs-i
а.6B464e+00Q
7.816iae+000
б.94771e+000
6.07925e+000
5.21073e+000
4.34232e+000
3.47386e+000
2.60539e+000
1.73693e+000
8.6B464e-001
0.00000e+000
a速度分布
3.35593e+003
压力/Pa
2.19988e+003
1.04384e+003
-1.12206e+002
-126825e+003
-2.42429e+003
-
3.58034e+O03
-4.73638e+O03
-5.89242e+003
-7.04847e+003
-8.20451e+003
b压力分布
图4中间轴面速度分布和压力分布
Watrr Powrr Vol.47No.
5
水力发电2021年5月
知，在［ = 0.9 m 水平截面中冷却油沟内油流速度较
推力轴瓦四周的油流速度低，速度变化比较均匀, 轴瓦进油侧和出油侧外缘比其他位置的油流速大,
瓦两侧油流的速度场近似对称分布。由图5 b 可知,
在［ = 0. 9 m 水平截面中压力沿着径向向外非线性增 大, 压力场在 推 力 轴 瓦 两 侧 近 似 对 称 分 布 , 故 轴 瓦
高苏宁的受力相对而言比较均匀。
由图6可知,［ = 0.55 m 处的水平截面中，油 的速度场和 压 力 场 近 似 对 称 分 布 , 压 力 沿 着 径 向 向
外基本呈现出逐渐增大的趋势。在圆柱形弹性油箱
的两侧沿着径向指向机组中心的方向，速度成条状
分布，且该处速度较大，弹性油箱周围其他位置速 度较小。
2.2温度场的计算结果及分析
2.2. 1不同工况下周期性平面温度场的计算结果 及分析
各工况下周期性平面温度场计算结果如图7所
示。由图7可知，初始瓦温、初始油温一定时，随 着水温按照24、18.2、10 ｛降低，周期性平面内最
低油温逐渐下降，但最低油温相差并不大。整体来
看，不同工况下周期性平面内温度变化趋势基本一 致 ， 从推力轴瓦外缘沿着径向向外 ， 油温逐渐降 低 ， 靠近冷却器处的油温最低，推力轴瓦下方区域油温
也在均匀减小，但整个周期性平面内最高油温和最
低油温相差不大，说明经过推力轴瓦和冷却器之后 的润滑油油温分布比较均匀。
速度/Hrs 」
6.O7078e+OOO 5.47090e+000 4.86303e+000 4.25515e+000 3.64727e+000 3.03939e+000 2.431516+000 1.823636+000 1.21576e+000 6.07878e-001 0.00000e+000
a 速度分布
压力/Pa
&66807e+0037.23698e+0035.80589e+0034.37479e+0032.94370e+0031.51260e+0038.15103e+001-1.34958e+003-2.78063e+003-4.21177e+003-S.64287e+003
b 压力分布
图5
% = 0. 9 N 处水平截面速度分布和压力分布
速度/nrs-i
2.29290e+000 2.06361e+000 1.83432e+000 1.60503e+000 1.37574e+000 1.14645e+000 9.17161e-001 6.87B71e-001 4.58580&-001 2.29290e-001 0.00000e+000
a 速度分布
压力/Pa
.2.97207e+003 fl 2.60756e+0032.24306e+0031.87B55e+0031.51405e+003I 1J4955e+003 7.85042e+0024.20538e+002j 5.60334e+001I -3.084710+002■ -6.72975B+002
b 压力分布
图6
% = 0. 55 N 处水平截面速度分布和压力分布
Watrr Poorr Vol 47 No.哺乳服装
5
第47卷第5期张 杨，等：水轮发电机组推力轴承油槽内流动及传热特性研究
c 工况3
a 工况1
b 工况2
图7 工况1 ~3周期性平面温度场
温度/K
3.2330+002 3.216e+002 3.198e+002 3.181e+002 3.163e+002 3.146e+002 3.128e+002 3.111e+002 3.094e+002 3.076M002 3.059e+002 3.041 e+002 3.024e+002 3.006e+002 2.989e+002 2.971 e+002a 工况1温度/K
3.233e+002 3.212e+002 3.190e+002 3.169e+002 3.148e+002 3.127e+002 3.105e+002 3.084e+002 3.063e+002 3.041 e +002 3.020e+002 2.999e+002 2.977e+002 2.956e+002 2.935e+002 2.913e+002b 工况2
温度/K ■ 3.233e+002■ 3.206e+002■ 3.180e+002 3.153e+0023.126e+002 3.099e+0023.072e+002 U 3.046e+002 H 3.019e+0022.992e+0022.965e+002 2.939e+002■ 2.912e+002■ 2.885e+002■ 2.858e+002■ 2.831 e+002
c 工况3条码检测
图8 工况1~3条件下％ = 0. 9 N 处水平截面温度分布
温度/K
3.216e+0023.200e+0023.183e+0023.167e+0023.151e+0023.134e+0023.118e+0023.102e+0023.086e+0023.069e+0023.053e+0023.037e+0023.020e+0023.004e+0022.988e+0022.971e+002温度/K
3.2288+002 3.207e+002 3.186e+002 3.165e+002 3.1449+002 3.123e+002 3.102e+002 3.081e +002 3.060e+002 3.039M002 3.018e+002 2.997e+002 2.976e+002 2.955e+002 2.934e+002 2.913e*002b 工况2
温度/K
3.246e+002 3.218e+002 3.191e+002 3.163e+002 3.135e+002 3.108e+002 3.080e+002 3.053e+002 3.025e+002 29970+002 2.970e+002 2.942e+002 2.914e+002 2.887e+002 2.859e+002 2.831e +002
C 工况3
a 工况1图9 工况1~3条件下％ = 0. 55 N 处水平截面温度分布
2.2.2 不同工况下水平截面温度场的计算结果及 分析
为进一步研究1/18油槽内的温度分布情况，取 ［ = 0.9 m 与［ = 0. 55 m 处水平截面的计算结果进行
分析，各工况下［ = 0.9 m 处水平截面的温度分布如 图8所示，各工况下［ = 0.55m 处水平截面的温度 分布如图9所示。
3种工况初始瓦温均高于初始油温，由图8可
知，在［ = 0.9 m 处水平截面中，每种工况下冷却油
沟中的油温均高于该水平截面内其他区域的油温,
在冷却器附近，油温下降梯度较大，其余区域油温 比较均匀；由图9可知，在［二0.55 m 处水平截面
中，每种工况下该水平截面整体油温分布都比较均
匀，冷却器附近油温下降梯度较大。
对3种工况中每种工况下［二0. 9 m 与［二0.55
m 处水平截面的计算结果进行对比可知，［ = 0.9 m
Warr Powr Vo .+47
家庭视频电话
No+5
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议