针对燃气轮机动力涡轮工作环境恶劣、承受负荷大等特点,为保证动力涡轮可靠性运行、满足强度要求,对动力涡轮静子(包括第三、四机匣,第三、四级导向器等零部件)在工作时的表面温度场进行了详细的数值计算。计算过程中,通过沿叶片径向设置扇形段并在扇形段的端面处设置循环对称边界条件,减少了单元与节点数目,简化了计算、节省了计算资源。计算结果可以用于动力涡轮部件的强度计算分析和冷热态尺寸链计算,对燃气轮机的结构强度设计和热力性能校核有一定的参考作用。 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.
葛仁超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. 机匣换热系数分段简图
表1给出了设计工况各段机匣外表面的辐射热流密度。燃气轮机机匣外表面的黑度为ε1 = 0.8,箱装体内表面的黑度为ε2 = 0.657;箱装体内表面温度为T2 = 270℃;环境温度为T0 = 45℃。由表1可知,由于燃气轮机机匣外表面温度变化较大,辐射换热热流密度变化也较大,因此必须采用分段计算辐射换热的方法,以提高温度场数值计算的准确度。
段号 | I | II | III | IV |
---|---|---|---|---|
T1 | 350℃ | 550℃ | 420℃ | 420℃ |
A1/A2 | 0.34 | 0.42 | 0.45 | 0.47 |
q | 2530 W/m2 | 14,349 W/m2 | 5487 W/m2 | 5449 W/m2 |
表1. 设计工况机匣外表面辐射换热热流密度
经过折算,得出各段机匣外表面综合换热系数如表2所示。
段号 | I | II | III | IV |
---|---|---|---|---|
h | 18.3 W/(m2·K) | 38.4 W/(m2·K) | 24.6 W/(m2·K) | 24.5 W/(m2·K) |
表2. 机匣外表面综合换热系数
由于外界环境温度随季节变化的不确定性,本小节分析了箱装体内环境温度分别为60℃、45℃和30℃时机匣外表面测温点的温度随环境温度变化的规律,数值计算结果如图3所示。
图3. 不同环境温度下机匣外表面温度场计算结果
由图3可知,当环境温度在45℃ ± 15℃内变化时,机匣外表面测温点的温度变化在±10℃范围内。
接触应力、接触界面状况(粗糙度、波纹度等)、接触界面之间的介质(液体、气体等)和接触件的材料性能(表面硬度、导热系数等)因素均会影响到接触界面之间的接触热阻 [
图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℃。因此,接触换热系数在一定范围内变化对机匣外表面的温度分布会产生影响;而在动力涡轮结构设计以及热膨胀量计算中,需要考虑换热系数变化的影响。
图4. 不同接触热阻下机匣外表面温度场计算结果
1) 第三级导向器内、外罩壳燃气侧,各级导叶组上、下缘板以及叶身表面,过度机匣燃气侧,级间表面,与箱装体内空气接触的机匣外表面,转接筒冷气侧,第四级隔板冷气侧,第三、四级护环蜂窝密封齿面,第四级隔板蜂窝密封齿面均采用第三类边界条件;
2) 各个金属表面的换热系数计算使用《航空发动机设计手册》中的相关公式,其它取自原型机计算数据以及空气系统的计算结果;
3) 机匣外表面区域考虑辐射换热影响;
4) 静子件周向端面定义为循环对称边界条件;
5) 热边界条件考虑燃气沿径向分布的不均匀性;
6) 导叶与机匣连接处以及各级机匣之间等静子件连接处接触界面给出接触热阻;
7) 空气物性参数取自《空气物理性质表》,以烟气物性参数代替燃气物性参数;
8) 燃气流道内对应的燃气温度和压力参数参照气动计算结果。
动力涡轮在设计工况下的温度场是冷热态尺寸链计算的基本依据,对动力涡轮结构设计具有很重要的指导意义。准确计算设计工况下动力涡轮的温度场是高性能、高可靠性动力涡轮结构设计的前提。
本文第3节根据经验公式以及气动计算的结果,计算了设计工况下动力涡轮静子件各换热表面的边界条件。滑油腔表面的边界条件根据动力涡轮支承环滑油的回油温度确定,考虑到动力涡轮转速较高,滑油温度增加20℃,各表面的边界条件如表3所示。
图5和图6给出了极限边界条件时,动力涡轮在设计工况下的各静子件温度场。
单位 | 表面编号 | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 2 | 3前 | 3后 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | ||
环境温度 | ℃ | 380 | 750 | 45 | 45 | 750 | 750 | 750 | 665 | 750 | 440 | 750 | 750 |
换热系数 | W/(m2·K) | — | 438 | 18.3 | 38.4 | 438 | 1500 | 386 | 1200 | 438 | 100 | 1640 | 2310 |
单位 | 表面编号 | ||||||||||||
12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | ||
环境温度 | ℃ | 750 | 750 | 750 | 45 | 45 | 631 | 630 | 630 | 630 | 630 | 631 | 522 |
换热系数 | W/(m2·K) | 796 | 796 | 755 | 24.6 | 24.5 | 934 | 1920 | 548 | 548 | 522 | 254 | 861 |
单位 | 表面编号 | ||||||||||||
24 | 25 | 26 | 27 | 28 | 29 | 30 | 31 | 32 | 33 | 34 | 35 | ||
环境温度 | ℃ | 487.2 | 630 | 414 | 340 | 350 | 45 | 45 | 45 | 45 | 45 | 45 | 487.2 |
换热系数 | W/(m2·K) | 164 | 1090 | 350 | — | 200 | 10 | 25.4 | 24.3 | 25.4 | 10 | 24.3 | 164 |
单位 | 表面编号 | ||||||||||||
36 | 37 | 38 | 39 | 40 | 41 | 42 | 43 | 45 | |||||
环境温度 | ℃ | 486 | 414 | 351 | 351 | 351 | 351 | 351 | 250 | 487.3 | |||
换热系数 | W/(m2·K) | 152 | 30 | 35 | 35 | 35 | 35 | 30 | — | 164 |
表3. 设计工况下各表面边界条件
图5. 设计工况第三、四级静子件温度场
图6. 设计工况机匣外表面温度场
由图5和图6可知,在设计工况下,动力涡轮静子件的温度分布在100℃~749℃之间,主要的高温零件是第三级导叶组,前、后外罩壳和内罩壳等。
本文以燃气轮机动力涡轮的气动计算结果为基本依据,计算了设计工况下高速动力涡轮静子件的温度场,同时定量分析了环境温度和接触换热系数对机匣外表面温度场的影响。计算结果可以作为动力涡轮结构强度和冷热态尺寸量计算的参考依据,为高速动力涡轮结构设计以及其它计算提供参考。
国家自然科学基金面上项目(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