设为首页 加入收藏 期刊导航 网站地图

Dynamical Systems and Control
Vol. 08  No. 04 ( 2019 ), Article ID: 32134 , 8 pages
10.12677/DSC.2019.84024

Temperature Field Calculation of a Certain Type of Gas Turbine Power Turbine Stator

Renchao Ge1, Fengshou Xiongcao2, Zhaochun Tang2

1The 3rd Naval Military Representative Office in Harbin, Harbin Heilongjiang

2College of Power Engineering, Naval University of Engineering, Wuhan Hubei

Received: Aug. 17th, 2019; accepted: Aug. 26th, 2019; published: Sep. 12th, 2019

ABSTRACT

In order to ensure the reliability and strength of the power turbine, the surface temperature field of the power turbine stator (including the third and fourth casings, the third and fourth stage guides, etc.) is calculated in detail. In the calculation process, the number of elements and nodes is reduced, the calculation is simplified and the computing resources are saved by setting sector section along the radial direction of blade and setting circular symmetric boundary conditions at the end of sector section. The calculation results can be used for strength calculation and analysis of power turbine components and dimension chain calculation in cold and hot state, which can be used as a reference for structural strength design and thermodynamic performance check of gas turbines.

Keywords:Gas Turbine, Power Turbine, Temperature Field

某型燃气轮机动力涡轮静子温度场计算

葛仁超1，熊曹丰收2，唐照淳2

1海军装备部驻沈阳地区军事代表局驻哈尔滨地区第三军事代表室，黑龙江 哈尔滨

2海军工程大学动力工程学院，湖北 武汉

收稿日期：2019年8月17日；录用日期：2019年8月26日；发布日期：2019年9月12日

摘 要

针对燃气轮机动力涡轮工作环境恶劣、承受负荷大等特点，为保证动力涡轮可靠性运行、满足强度要求，对动力涡轮静子(包括第三、四机匣，第三、四级导向器等零部件)在工作时的表面温度场进行了详细的数值计算。计算过程中，通过沿叶片径向设置扇形段并在扇形段的端面处设置循环对称边界条件，减少了单元与节点数目，简化了计算、节省了计算资源。计算结果可以用于动力涡轮部件的强度计算分析和冷热态尺寸链计算，对燃气轮机的结构强度设计和热力性能校核有一定的参考作用。

关键词 :燃气轮机，动力涡轮，温度场

Copyright © 2019 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

在燃气轮机设计中，温度场计算难度较大 [1] 。本文对机匣外表面进行了分段处理 [2] [3] [4] ，计算了各段的辐射热流密度和综合换热系数，分析了不同环境温度对机匣温度以及不同接触热阻对机匣温度场的影响，并选定接触热阻，按燃气轮机实际运行状态确定了计算边界条件，并对静子温度场进行了数值计算与分析。

2. 动力涡轮静子件结构

某型燃气轮机动力涡轮的第三、四级静子件结构如图1所示。

1-第三级导向器机匣；2-前外罩壳；3-保温层；4-后外罩壳；5-第三级导叶组；6-第三级护环；7-第四级导向器机匣；8-第四级导叶组；9-第四级护环；10-过度机匣；11-转接筒；12-保温层；13-内罩壳；14-法兰；15-第四级隔板

Figure 1. Structural sketch of a certain type power turbine of gas turbine’s 3rd and 4th stators

图1. 某型燃气轮机动力涡轮第三、四级静子件结构示意图

其中，温度场计算的对象主要包括了第三、四级机匣，过渡机匣，第三、四级导叶组，内、外罩壳，转接筒，法兰，第三、四级护环，第四级隔板等零件。

3. 温度场数值计算时的重要参数计算方法

动力涡轮静子件的温度场数值计算本质为热平衡的计算分析，热边界条件的选择和材料的属性对数值计算结果影响非常重要 [5] [6] 。

动力涡轮静子件既有被高温、高速燃气包围的部分，也有在箱装体中暴露于空气中的部分。因此，在数值计算前要开展必要的热流动方向分析 [7] [8] 。第三级导向器的内、外罩壳，第三、四级导叶组的上、下缘板和叶身表面部分，以及机匣的燃气侧是主要的热量来源。

3.1. 机匣外表面换热系数

机匣外表面的换热同时包含了对流换热与辐射换热。箱装体内的引射空气可以通过对流换热的方式带走机匣外表面上的热量；同时机匣外表面与箱装体内表面会发生辐射换热 [9] [10] 。根据辐射换热计算公式，换热量与表面温度关系很大。机匣外表面温度分布非常不均匀，温差最高可以达到200℃以上，因此必须对机匣外表面进行分段处理，以便计算辐射换热量。在数值计算中将机匣外表面分成4段，如图2所示，并且分别计算了各段的辐射热流密度。

Figure 2. Piecewise sketch of case heat transfer coefficient

图2. 机匣换热系数分段简图

表1给出了设计工况各段机匣外表面的辐射热流密度。燃气轮机机匣外表面的黑度为ε1 = 0.8，箱装体内表面的黑度为ε2 = 0.657；箱装体内表面温度为T2 = 270℃；环境温度为T0 = 45℃。由表1可知，由于燃气轮机机匣外表面温度变化较大，辐射换热热流密度变化也较大，因此必须采用分段计算辐射换热的方法，以提高温度场数值计算的准确度。

Table 1. Heat flux density of radiation heat transfer of external case surface power under design conditions

表1. 设计工况机匣外表面辐射换热热流密度

经过折算，得出各段机匣外表面综合换热系数如表2所示。

Table 2. Comprehensive heat transfer coefficient of external case surface

表2. 机匣外表面综合换热系数

3.2. 环境温度对机匣外表面温度场的影响

由于外界环境温度随季节变化的不确定性，本小节分析了箱装体内环境温度分别为60℃、45℃和30℃时机匣外表面测温点的温度随环境温度变化的规律，数值计算结果如图3所示。

Figure 3. Outer surface temperature field of case under different environment temperature

图3. 不同环境温度下机匣外表面温度场计算结果

由图3可知，当环境温度在45℃ ± 15℃内变化时，机匣外表面测温点的温度变化在±10℃范围内。

3.3. 接触热阻对机匣外表面温度场的影响

接触应力、接触界面状况(粗糙度、波纹度等)、接触界面之间的介质(液体、气体等)和接触件的材料性能(表面硬度、导热系数等)因素均会影响到接触界面之间的接触热阻 [11] 。此外，目前尚未有一个能准确计算接触热阻的公式，许多经验公式本身又有一定的适用范围，因此要准确的给出动力涡轮静子件接触界面之间的接触热阻是很困难的。在本节中分析了不同接触热阻对机匣外表面温度场的影响，供动力涡轮结构设计和热膨胀量计算参考。

图4为设计工况时，接触换热系数(接触热阻的倒数)分别为2000 W/(m2·K)、1500 W/(m2·K)和1000 W/(m2·K)时机匣外表面测温点温度分布。由图4可知，当接触换热系数在1500 W/(m2·K)基础上提高到500 W/(m2·K)时，机匣外表面温度约提高4℃；当接触换热系数降低到500 W/(m2·K)时，机匣外表面温度约降低8℃。因此，接触换热系数在一定范围内变化对机匣外表面的温度分布会产生影响；而在动力涡轮结构设计以及热膨胀量计算中，需要考虑换热系数变化的影响。

Figure 4. Outer surface temperature field of case under different contact thermal resistance

图4. 不同接触热阻下机匣外表面温度场计算结果

3.4. 数值计算时的边界条件

1) 第三级导向器内、外罩壳燃气侧，各级导叶组上、下缘板以及叶身表面，过度机匣燃气侧，级间表面，与箱装体内空气接触的机匣外表面，转接筒冷气侧，第四级隔板冷气侧，第三、四级护环蜂窝密封齿面，第四级隔板蜂窝密封齿面均采用第三类边界条件；

2) 各个金属表面的换热系数计算使用《航空发动机设计手册》中的相关公式，其它取自原型机计算数据以及空气系统的计算结果；

3) 机匣外表面区域考虑辐射换热影响；

4) 静子件周向端面定义为循环对称边界条件；

5) 热边界条件考虑燃气沿径向分布的不均匀性；

6) 导叶与机匣连接处以及各级机匣之间等静子件连接处接触界面给出接触热阻；

7) 空气物性参数取自《空气物理性质表》，以烟气物性参数代替燃气物性参数；

8) 燃气流道内对应的燃气温度和压力参数参照气动计算结果。

4. 温度场数值计算结果

动力涡轮在设计工况下的温度场是冷热态尺寸链计算的基本依据，对动力涡轮结构设计具有很重要的指导意义。准确计算设计工况下动力涡轮的温度场是高性能、高可靠性动力涡轮结构设计的前提。

4.1. 边界条件

本文第3节根据经验公式以及气动计算的结果，计算了设计工况下动力涡轮静子件各换热表面的边界条件。滑油腔表面的边界条件根据动力涡轮支承环滑油的回油温度确定，考虑到动力涡轮转速较高，滑油温度增加20℃，各表面的边界条件如表3所示。

4.2. 计算结果及分析

图5和图6给出了极限边界条件时，动力涡轮在设计工况下的各静子件温度场。

Table 3. Relationship of power and fuel consumption rate under different speed

表3. 设计工况下各表面边界条件

Figure 5. Temperature field of the 3rd and 4th stators under design conditions

图5. 设计工况第三、四级静子件温度场

Figure 6. Outer surface temperature field of case under design conditions

图6. 设计工况机匣外表面温度场

由图5和图6可知，在设计工况下，动力涡轮静子件的温度分布在100℃~749℃之间，主要的高温零件是第三级导叶组，前、后外罩壳和内罩壳等。

5. 结束语

本文以燃气轮机动力涡轮的气动计算结果为基本依据，计算了设计工况下高速动力涡轮静子件的温度场，同时定量分析了环境温度和接触换热系数对机匣外表面温度场的影响。计算结果可以作为动力涡轮结构强度和冷热态尺寸量计算的参考依据，为高速动力涡轮结构设计以及其它计算提供参考。

基金项目

国家自然科学基金面上项目(51179196)。

文章引用

葛仁超,熊曹丰收,唐照淳. 某型燃气轮机动力涡轮静子温度场计算
Temperature Field Calculation of a Certain Type of Gas Turbine Power Turbine Stator[J]. 动力系统与控制, 2019, 08(04): 222-229. https://doi.org/10.12677/DSC.2019.84024

参考文献

  1. 1. Granovski, A.V., Kostege, M.K., Ivanov, M.J., et al. (2001) Simulation of Temperature Field Redistribution through Multistage Cooled Turbines. In: ASME Turbo Expo 2001: Power for Land, Sea, and Air, American Society of Mechanical Engineers, New York, 1658-1662.
    https://doi.org/10.1115/2001-GT-0576

  2. 2. 王旭, 王世安, 张文平, 等. 燃气轮机涡轮后机匣温度场及应力分析[J]. 热能动力工程, 2004, 19(6): 558-561.

  3. 3. 李朝阳, 张艳春. 基于燃机涡轮盘三维瞬态温度及应力场计算分析[J]. 动力工程, 2006, 25(2): 211-214.

  4. 4. 李朝阳, 张艳春, 何光新. 燃气轮机涡轮盘温度及应力场计算分析[J]. 汽轮机技术, 2006, 48(1): 7-10.

  5. 5. Wang, Z.F., Yan, P.G., Guo, Z.Y., et al. (2008) BEMFDM Conjugate Heat Transfer Analysis of a Two-Dimensional Air-Cooled Turbine Blade Boundary Layer. Journal of Thermal Science, 17, 199-206.
    https://doi.org/10.1007/s11630-008-0199-2

  6. 6. 许文锋, 朱惠人, 张丽. 带冷却的涡轮叶片温度场耦合计算工程方法研究[J]. 汽轮机技术, 2008, 50(3): 201-203.

  7. 7. Liu, K. and Zhang, B.C. (2013) Experiment on Outlet Temperature Field of a Gas Turbine Combustor. Journal of Aerospace Power, 28, 1387-1391.

  8. 8. Liu, K. and Xu, L. (2012) Study on Combustor Outlet Tem-perature Field of Gas Turbine. In: 2nd International Conference on Applied Mechanics and Mechanical Engineering, Trans Tech Publications, Zurich, 962-966.
    https://doi.org/10.4028/www.scientific.net/AMM.138-139.962

  9. 9. 迟重然, 任静, 蒋洪德. 燃机叶片平行肋扰流内冷通道传热特性研究Part 1 [J]. 工程热物理学报, 2013, 34(4): 624-627.

  10. 10. 迟重然, 任静, 蒋洪德. 燃机叶片平行肋扰流内冷通道传热特性研究Part 2 [J]. 工程热物理学报, 2014, 35(1): 42-45.

  11. 11. 李钰洁, 刘永葆, 高洁. 间隙流动对涡轮叶片温度场分布影响的数值研究[J]. 海军工程大学学报, 2013, 25(6): 24-29.

期刊菜单