一种轴流压气机亚音速叶栅造型关键角度参数确定方法与流程

1.本发明涉及的是一种燃气轮机设计方法，具体地说是压气机设计方法。


背景技术：

2.压气机作为燃气轮机的关键核心部件之一，其性能优劣对燃气轮机的技术指标实现起着决定性作用。纵观整个压气机的发展历史，所有的压气机设计都是建立在对通流能力、压比、效率以及喘振裕度的追求之上的。时至今日，随着燃气轮机性能的不断提升，对压气机的气动性能需求也自然地演变成了更高的压比、更高的效率和更大的喘振裕度，如今的压气机气动设计体系也正是基于以上需求而形成的。
3.随着对内部流动的认识不断深入，压气机气动设计体系也在飞速发展，越来越多的技术手段已经被应用到压气机的气动设计当中。所谓压气机气动设计，其本质就是将压气机气动性能需求转换成几何造型的过程。而叶片的几何造型设计则是这个过程中最重要的关键环节之一。目前以s1、s2两类流面理论为基础的准三维气动设计作为现阶段压气机气动设计体系的核心部分，如何将准三维通流设计结果转换为叶片的几何造型，其中最重要的环节就在于叶片造型关键角度参数的确定。可以说，合理、准确的叶片造型角度参数直接关系着压气机内部的级间流动匹配，对压气机气动性能指标的实现起着至关重要的作用。


技术实现要素：

4.本发明的目的在于提供能解决轴流压气机气动设计中亚音速叶栅造型的关键角度参数确定问题的一种轴流压气机亚音速叶栅造型关键角度参数确定方法。
5.本发明的目的是这样实现的：
6.本发明一种轴流压气机亚音速叶栅造型关键角度参数确定方法，其特征是：
7.(1)确定进出口几何角，叶栅进口、出口几何角β
1k
、β
2k
分别为叶型中弧线在前、尾缘点切线与圆周方向的夹角，根据轴流压气机通流反问题设计结果，获得叶栅进口气流角β1与出口气流角β2，计算叶栅进口攻角i和出口落后角δ，从而确定β
1k
＝β1+i，β
2k
＝β2+δ；
8.(2)计算叶型弯角，在确定了进出口几何角之后，计算得到叶型弯角ε＝β
2k-β
1k
；
9.(3)确定进出口楔角，叶栅进口、出口楔角χ1、χ2分别为叶型中弧线在前、尾缘点切线与弦线的夹角，根据叶型中弧线形式，确定进口、出口楔角χ1、χ2值；
10.(4)计算叶型安装角，根据进口几何角β
1k
和进口楔角χ1，计算叶型安装角γ＝β
1k
+χ1。
11.本发明还可以包括：
12.1、步骤(1)中所述的计算叶栅进口攻角i和出口落后角δ，采用如下方式：
13.进口攻角i：
[0014][0015]
式中，b/t为叶型对应位置的叶栅稠度；为轮毂比；为叶型对应位置的相对半径；为相对平均半径，m
b/t＝1.0
为叶型稠度为1时攻角随弯角变化率，是b/t的函数，m
b/t＝1.0
＝-0.1686(b/t)2+0.9393(b/t)+0.2302，若计算结果小于0.7，则取m
b/t＝1.0
＝0.7；δβ为气流转折角，δβ＝β
2-β1；为叶型稠度为1时气流转折角名义值，是β2的函数，δi为攻角修正量，根据叶栅平均半径参数来计算，其中(b/t)m为叶栅平均半径处稠度，β
1m
为叶栅平均半径处进口气流角，ωm为叶栅平均半径处反动度；
[0016]
出口落后角δ：
[0017][0018]
式中，m为叶型中弧线系数，其中为叶型相对最大挠度位置；为落后角修正系数，其取值与叶栅轮毂比有关：
[0019][0020]
上式中c
1a
为叶型出口轴向速度，c
2a
为叶型进口轴向速度。
[0021]
2、步骤(3)中所述的根据叶型中弧线形式，确定进口、出口楔角χ1、χ2值，采用如下方式：
[0022]
当采用圆弧中弧线时，χ1＝χ2＝0.5ε；
[0023]
当采用抛物线中弧线时，χ1＝0.6ε，χ2＝0.4ε。
[0024]
本发明的优势在于：
[0025]
1、本发明能够快速、准确地计算出轴流压气机大部分常用亚音速叶栅在设计状态
下的造型关键角度参数，有效提高了设计精度，提高了压气机的气动性能。
[0026]
2、本发明有效减少了通过三维计算来逐级修正多级轴流压气机各级亚音速叶片角度造型参数的次数，实现了压气机叶片角度造型参数的流程化计算，可节省大量的设计迭代时间，缩短设计周期。
[0027]
3、本发明不仅局限于燃气轮机轴流压气机，同样适用于各种工业用轴流压缩机、航空发动机轴流压气机/风扇的气动设计过程。
附图说明
[0028]
图1为本发明的流程图；
[0029]
图2为轴流压气机亚音速叶栅造型关键角度参数定义图。
具体实施方式
[0030]
下面结合附图举例对本发明做更详细地描述：
[0031]
结合图1-2，本发明一种轴流压气机亚音速叶栅造型关键角度参数确定方法，通过以下步骤实现：
[0032]
步骤一：确定进出口几何角。叶栅进口、出口几何角β
1k
、β
2k
分别为叶型中弧线在前、尾缘点切线与圆周方向的夹角。根据轴流压气机通流反问题设计结果，获得叶栅进口气流角β1与出口气流角β2，计算叶栅进口攻角i和出口落后角δ，从而确定β
1k
＝β1+i，β
2k
＝β2+δ；
[0033]
步骤二：计算叶型弯角。在确定了进出口几何角之后，即可计算得到叶型弯角ε＝β
2k-β
1k
；
[0034]
步骤三：确定进出口楔角。叶栅进口、出口楔角χ1、χ2分别为叶型中弧线在前、尾缘点切线与弦线的夹角。根据叶型中弧线形式，确定进口、出口楔角χ1、χ2值；
[0035]
步骤四：计算叶型安装角。根据进口几何角β
1k
和进口楔角χ1，可以计算叶型安装角γ＝β
1k
+χ1。
[0036]
进一步的：步骤一中所述的“计算叶栅进口攻角i和出口落后角δ”，采用如下方式：
[0037]
进口攻角i：
[0038][0039]
式中，b/t为叶型对应位置的叶栅稠度；为轮毂比；为叶型对应位置的相对半径；为相对平均半径，m
b/t＝1.0
为叶型稠度为1时攻角随弯角变化率，是b/t的函数，m
b/t＝1.0
＝-0.1686(b/t)2+0.9393(b/t)+0.2302，若计算结果小于0.7，则取m
b/t＝1.0
＝0.7；δβ为气流转折角，δβ＝β
2-β1；为叶型稠度为1时气流转折角名义值，是β2的
函数，δi为攻角修正量，根据叶栅平均半径参数来计算，其中(b/t)m为叶栅平均半径处稠度，β
1m
为叶栅平均半径处进口气流角，ωm为叶栅平均半径处反动度。
[0040]
出口落后角δ：
[0041][0042]
式中，m为叶型中弧线系数，其中为叶型相对最大挠度位置；为落后角修正系数，其取值与叶栅轮毂比有关：
[0043][0044]
上式中c
1a
为叶型出口轴向速度，c
2a
为叶型进口轴向速度。
[0045]
进一步的：步骤三中所述的“根据叶型中弧线形式，确定进口、出口楔角χ1、χ2值”，采用如下方式：
[0046]
当采用圆弧中弧线时，χ1＝χ2＝0.5ε；
[0047]
当采用抛物线中弧线时，χ1＝0.6ε，χ2＝0.4ε。