CN113513373B - 一种变几何涡轮一维气动设计方法 - Google Patents

一种变几何涡轮一维气动设计方法 Download PDF

Info

Publication number
CN113513373B
CN113513373B CN202110772097.8A CN202110772097A CN113513373B CN 113513373 B CN113513373 B CN 113513373B CN 202110772097 A CN202110772097 A CN 202110772097A CN 113513373 B CN113513373 B CN 113513373B
Authority
CN
China
Prior art keywords
guide vane
outlet
vane
variable geometry
flow
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110772097.8A
Other languages
English (en)
Other versions
CN113513373A (zh
Inventor
高杰
霍东晨
李彦静
赵天笑
闫睿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202110772097.8A priority Critical patent/CN113513373B/zh
Publication of CN113513373A publication Critical patent/CN113513373A/zh
Application granted granted Critical
Publication of CN113513373B publication Critical patent/CN113513373B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D9/00Stators
    • F01D9/02Nozzles; Nozzle boxes; Stator blades; Guide conduits, e.g. individual nozzles
    • F01D9/04Nozzles; Nozzle boxes; Stator blades; Guide conduits, e.g. individual nozzles forming ring or sector
    • F01D9/041Nozzles; Nozzle boxes; Stator blades; Guide conduits, e.g. individual nozzles forming ring or sector using blades
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D17/00Regulating or controlling by varying flow
    • F01D17/10Final actuators
    • F01D17/12Final actuators arranged in stator parts
    • F01D17/14Final actuators arranged in stator parts varying effective cross-sectional area of nozzles or guide conduits
    • F01D17/16Final actuators arranged in stator parts varying effective cross-sectional area of nozzles or guide conduits by means of nozzle vanes
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D21/00Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D25/00Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Mechanical Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Control Of Turbines (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)

Abstract

本发明的目的在于提供一种变几何涡轮一维气动设计方法,利用常见的平均中径法,完成常规涡轮的气动设计;进而利用流量与可转导叶转角关联式估算出可转导叶转角范围,通过一维总静流函数法并结合可调导叶间隙泄漏损失模型,计算出变几何涡轮在各个转角下效率,最后将各转角下效率平均得出变几何涡轮平均效率;若变几何涡轮平均效率满足设计要求,则变几何涡轮一维设计完成;若不符合调整变几何涡轮一维气动部分设计参数,直至变几何涡轮性能满足设计要求。本发明可综合判断变几何涡轮各工况效率,从而得到一种各工况下均具有较高效率的变几何涡轮气动设计方案,从而提高变几何涡轮变工况性能;此外,本发明也可明显提高变几何涡轮一维设计速度。

Description

一种变几何涡轮一维气动设计方法
技术领域
本发明涉及的是一种燃气轮机设计方法,具体地说是涡轮设计方法。
背景技术
燃气轮机涡轮的气动设计工作是一个从低维到高维逐步设计及优化的过程,低维设计结果作为高维设计的初值及基础,为了完成变几何涡轮气动设计工作,首先要解决的问题在于如何选取适合于变几何涡轮的气动设计参数,这就需要在一维空间中进行设计及参数选取工作,而涡轮一维损失模型作为一维计算的核心及基础,其准确性直接决定了一维设计的成功与否,通过修正可调导叶间隙泄漏损失模型,从而得到适用于并且准确预估变几何涡轮损失模型。
常规涡轮的一维参数选取规律研究通常采用平均中径法,将一维计算设置在叶片排进出口平均中径上,根据已知的几何参数和边界条件,通过基本气动关系式对叶片排逐排求解,并结合损失模型建立涡轮效率及涡轮气动参数如载荷系数,流量系数,反动度,轴向速比及几何参数,如导叶及动叶展弦比,叶栅稠度等联系,从而分析涡轮效率与上述参数相关性,并绘制出Smith等关系图直观的寻找出涡轮设计参数中的最优选择。变几何涡轮与常规涡轮有着极大的区别,设计工况涡轮效率并不足以描述变几何涡轮性能的好坏,变几何涡轮除常规设计工况外,还必须保证其非设计工况依旧保持较高的涡轮效率,因此与常规涡轮一维参数选取规律研究不同的是,变几何涡轮一维参数选取规律研究必须综合考虑不同转角下涡轮效率,即通过一维参数完成对设计工况下涡轮气动设计后,再利用考虑可调导叶间隙泄漏损失模型的一维性能预测方法,评估不同转角下涡轮效率,求取各个转角工况涡轮效率平均值,平均效率最高的区域即为最佳一维参数区域。
为了探究变几何涡轮一维参数选取规律,首先需要解决的问题在于变几何涡轮导叶旋转之后涡轮气动设计参数如反动度等的变化规律,从而不同工况下一维参数与转角建立联系,最终获得不同转角下涡轮效率值,通过分析设计参数与不同转角下的平均效率进而确定变几何涡轮一维参数的选取范围,进而探讨和发展变几何涡轮一维气动设计方法。
发明内容
本发明的目的在于提供能使变几何涡轮所有工况平均效率达到较高水平的一种变几何涡轮一维气动设计方法。
本发明的目的是这样实现的:
本发明一种变几何涡轮一维气动设计方法,其特征是:
(1)给定初始几何参数与气动参数,包括动静叶展弦比、最大厚度比弦长、尾缘厚度比喉部宽度、轮缘功、进口总温总压、导叶进口气流角、载荷系数、流量系数、反动度、轴向速比、进出口中径比等参数。
(2)根据涡轮一维几何参数求解动静叶轮毂比、叶高等几何参数,并通过一维总静流函数法分析导叶转角转动后,涡轮级反动度、载荷系数、流量系数、轴向速比、气流角的变化。
(3)将导叶转动后,气动参数与几何参数输入到可调导叶间隙泄漏损失模型中,计算各种转角变化后动静叶速度损失系数,即可求出变几何涡轮各转角下涡轮效率。
本发明还可以包括:
1、(a)利用可转导叶转角与变几何涡轮流量关联式
Figure GDA0003986670380000021
预估所需的可转导叶旋转角,其中α1导叶未旋转的出口气流角,seita导叶旋转角,m1为变几何涡轮导叶转角变化前的流量,m1z为变几何涡轮导叶转角变化后的流量;
(b)将可转导叶旋转角分为n份,利用一维总静流函数法并综合可调导叶间隙泄漏损失模型分别计算各转角工况效率;
(c)基于步骤(b)计算变几何涡轮在不同转角下的平均效率;
(d)如果步骤(c)得到的性能指标不符合预定设计目标,则修改常规涡轮一维设计循环中的一维设计参数,并循环(a)-(c),直至涡轮气动性能达到预定设计目标。
2、对于变几何涡轮,一维总静流函数法的推导过程为:
一维总静流函数法分别包含了导叶流量m1和动叶流量m2随导叶总静膨胀比
Figure GDA0003986670380000022
的变化规律,将涡轮级通道看作一维喷管,导叶入口总压P0 *,导叶入口总温T0 *,动叶出口静压P1。导叶流量定义为:
m1=ρ1c1A1sinα1
其中ρ1为导叶通道内气体密度,c1为导叶出口绝对速度,A1为导叶出口截面积,α1为导叶绝对出口气流角;
根据理想气体状态方程P0 *=ρ0 *RgT0 *及多变过程
Figure GDA0003986670380000023
及T0 *=T1 *得公式:
Figure GDA0003986670380000031
Figure GDA0003986670380000032
其中ns,s为导叶气体多变指数,T1 *为导叶出口总温,ρ0 *为导叶进口密度,T1为导叶出口静温,Cp为气体定压比热容,h0 *为导叶进口滞止焓,h1为导叶出口焓值,k为绝热指数。
综合以上公式,可得出
Figure GDA0003986670380000033
即可得到导叶流量随导叶总静膨胀比变化规律;
通过下式建立导叶多变过程指数与速度损失系数
Figure GDA0003986670380000034
的关系:/>
Figure GDA0003986670380000035
根据连续方程及能量守恒方程,可得到λc1和动叶进口静温T1
Figure GDA0003986670380000036
T* 1=T* 0
T1=T* 1τ(λc1)
其中Rg为气体常数,λc1为导叶出口绝对速度因数,ε(λc1)为基于导叶出口绝对速度因数的密度气动函数,τ(λc1)为基于导叶出口绝对速度因数的温度气动函数。
根据动叶进口连续方程,可由λc1得到λω1,并根据λω1即可得出动叶进口相对总温及相对总压:
Figure GDA0003986670380000037
P* ω1=P1/π(λω1)
T* ω1=T1/τ(λω1)
其中λω1为基于动叶进口相对速度因数,P* ω1为动叶进口相对总压,T* ω1为动叶进口相对总温,π(λω1)为基于动叶进口相对速度因数的压力气动函数,τ(λω1)为基于动叶进口相对速度因数的温度气动函数。
并且根据上述公式及气动函数可求出导叶出口速度因数
Figure GDA0003986670380000041
静温T1,对于动叶而言,可以在相对坐标系下同样通过下述公式推导出动叶流量随导叶总静膨胀比,速度损失系数,出口相对气流角的变化规律:
m2=ρ2ω2A2sinβ2
其中ρ2为动叶出口气体密度,ω2为动叶出口相对速度,A2为动叶出口截面积。
根据理想气体状态方程及P1=ρ1RgT1及多变过程
Figure GDA0003986670380000042
可得
Figure GDA0003986670380000043
Figure GDA0003986670380000044
Figure GDA0003986670380000045
其中ns为气体多变指数,P2为动叶出口静压,其中ns,r为动叶气体多变指数,ω2为动叶出口相对速度,ns,r为动叶气体多变指数,hω1 *为动叶进口相对滞止焓,h2动叶出口焓值。
联立上述公式即可得到动叶流量随涡轮导叶总静膨胀比的变化规律:
Figure GDA0003986670380000046
其中β2为动叶几何出气角,其他参数均可由已知参数求得,因此可得到动叶流量随导叶出口总静膨胀比变化规律。动叶流量与导叶流量交点即为涡轮级工作点。
3、对于变几何涡轮,变几何涡轮导叶转角变化与流量关联式的推导过程为:
导叶几何调节后,直接改变导叶出口几何出气角,通过几何关系得到旋转后导叶几何出气角α1z,α1z为α1与旋转角seita之和,其中seita为正值时表示导叶开大,seita为负值则相反,在忽略落后角的条件下旋转后导叶绝对出口气流角与动叶相对出口气流角如下式所示:
α1z=α1+seita
β2z=γ2
对于变几何涡轮,变几何涡轮可调导叶间隙泄漏损失模型的修正公式为:
YTI=Kz(Ygap+Ymix)
Figure GDA0003986670380000047
Figure GDA0003986670380000051
其中YTI为间隙泄漏损失,Kz为旋转轴修正系数,Ygap为间隙内损失,Ymix为泄漏流与主流掺混损失,C为叶片弦长,S为叶片栅距,H为叶片高度,μd为叶尖间隙流量系数,CL升力系数,αm为通过叶片的气流角平均值,KE为湍动能损失系数,KG为叶尖间隙损失系数;
旋转轴修正系数Kz为:
Figure GDA0003986670380000052
其中变量x为无量纲旋转轴位置,变量y为无量纲轴径,p00、p10、p01、p20、p11、p02、p30、p21、p12、p03、p31、p22、p13、p04均为拟合多项式的系数,其中字母p后的第一位数字表示变量x的次幂,字母p后的第二位数字表示变量y的次幂。
4、变几何涡轮导叶出口气流角的范围在15°-25°度之间,变几何涡轮最佳反动度区域基本在0.25-0.6之间,流量系数在0.3-0.5之间,轴向速比选择范围为1-1.4之间。
本发明的优势在于:本发明根据变几何涡轮设计特点,通过保持导叶旋转前后导叶入口总压、总温和动叶出口静压不变,在一维总静流函数法中仅改变中的导叶出口气流角,通过二分法即可求出不同转角下导叶与动叶工作点的流量、气流角、及各个气动参数。此设计方法可以方便、准确评估变几何涡轮在不同转角下的气动性能,给出变几何涡轮速度三角形关键参数的选取方法或一维参数选取规律,进而探讨和发展变几何涡轮一维气动设计方法。
附图说明
图1为本发明的流程图;
图2为叶片排一维管流模型示意图;
图3为可转导叶与动叶流量匹配示意图;
图4为可转导叶旋转示意图;
图5为拟合后与拟合前不同轴位置及轴径下Kz;
图6为不同转角工况平均效率与设计效率随反动度变化曲线图。
具体实施方式
下面结合附图对本发明做更详细地描述:
结合图1-6,本发明提供了一种变几何涡轮的一维气动设计方法,该方法中,综合考虑可调导叶间隙泄漏损失模型的一维性能预测方法、一维总静流函数法,以及几何参数求解效率计算公式编写变几何涡轮一维参数选取规律程序,具体步骤如下:
1、给定初始几何参数与气动参数,包括动静叶展弦比、最大厚度比弦长、尾缘厚度比喉部宽度、轮缘功、进口总温总压、导叶进口气流角、载荷系数、流量系数、反动度、轴向速比、进出口中径比等参数。
2、根据涡轮一维几何参数求解动静叶轮毂比、叶高等几何参数,并通过一维总静流函数法分析导叶转角转动后,涡轮级反动度、载荷系数、流量系数、轴向速比、气流角的变化。
3、将导叶转动后,气动参数与几何参数输入到可调导叶间隙泄漏损失模型中,计算各种转角变化后动静叶速度损失系数,即可求出变几何涡轮各转角下涡轮效率。
在所述常规涡轮一维设计循环之后还包括如下步骤:
(1)根据实际流量变化需求,并利用公式(16),求得变几何涡轮可转导叶转角变化范围。
(2)将可转导叶转角分为若干转角工况,利用一维总静流函数法,求取变几何涡轮
可转导叶转角变化后变几何涡轮各转角工况效率。
(3)将步骤(2)求得的各转角工况效率加权平均,求得变几何涡轮平均效率。
(4)如果步骤(3)中得到的性能指标不符合预定设计目标,则重复步骤(1)~(3),
直至步骤(3)得到的性能指标与预定设计目标相符合。
一维总静流函数法分别包含了导叶流量m1和动叶流量m2随导叶总静膨胀比
Figure GDA0003986670380000061
的变化规律,将涡轮级通道看作一维喷管,导叶入口总压P0 *,导叶入口总温T0 *,动叶出口静压P1。导叶流量定义为:
m1=ρ1c1A1sinα1 (1)
其中ρ1为导叶通道内气体密度,c1为导叶出口绝对速度,A1为导叶出口截面积,α1为导叶绝对出口气流角;
根据理想气体状态方程P0 *=ρ0 *RgT0 *及多变过程
Figure GDA0003986670380000062
及T0 *=T1 *得公式:/>
Figure GDA0003986670380000063
Figure GDA0003986670380000071
其中ns,s为导叶气体多变指数,T1 *为导叶出口总温,ρ0 *为导叶进口密度,T1为导叶出口静温,Cp为气体定压比热容,h0 *为导叶进口滞止焓,h1为导叶出口焓值,k为绝热指数。
综合以上公式,可得出
Figure GDA0003986670380000072
即可得到导叶流量随导叶总静膨胀比变化规律;
通过下式建立导叶多变过程指数与速度损失系数
Figure GDA0003986670380000073
的关系:
Figure GDA0003986670380000074
根据连续方程及能量守恒方程,可得到λc1和动叶进口静温T1
Figure GDA0003986670380000075
Figure GDA0003986670380000076
其中Rg为气体常数,λc1为导叶出口绝对速度因数,ε(λc1)为基于导叶出口绝对速度因数的密度气动函数,τ(λc1)为基于导叶出口绝对速度因数的温度气动函数。
根据动叶进口连续方程,可由λc1得到λω1,并根据λω1即可得出动叶进口相对总温及相对总压:
Figure GDA0003986670380000077
Figure GDA0003986670380000078
其中λω1为基于动叶进口相对速度因数,P* ω1为动叶进口相对总压,T* ω1为动叶进口相对总温,π(λω1)为基于动叶进口相对速度因数的压力气动函数,τ(λω1)为基于动叶进口相对速度因数的温度气动函数。
并且根据上述公式及气动函数可求出导叶出口速度因数
Figure GDA0003986670380000079
静温T1,对于动叶而言,可以在相对坐标系下同样通过下述公式推导出动叶流量随导叶总静膨胀比,速度损失系数,出口相对气流角的变化规律:
m2=ρ2ω2A2sinβ2 (10)
其中ρ2为动叶出口气体密度,ω2为动叶出口相对速度,A2为动叶出口截面积。
根据理想气体状态方程及P1=ρ1RgT1及多变过程
Figure GDA0003986670380000081
可得
Figure GDA0003986670380000082
Figure GDA0003986670380000083
Figure GDA0003986670380000084
其中ns为气体多变指数,P2为动叶出口静压,其中ns,r为动叶气体多变指数,ω2为动叶出口相对速度,ns,r为动叶气体多变指数,hω1 *为动叶进口相对滞止焓,h2动叶出口焓值。
联立上述公式即可得到动叶流量随涡轮导叶总静膨胀比的变化规律:
Figure GDA0003986670380000085
其中β2为动叶几何出气角,其他参数均可由已知参数求得,因此可得到动叶流量随导叶出口总静膨胀比变化规律。动叶流量与导叶流量交点即为涡轮级工作点。
如图3,动叶流量与导叶流量交点即为涡轮级工作点。本发明应用上述一维总静流函数法可迅速确定涡轮级工作点,此方法可迅速确定导叶旋转后各个转角工况下流量、反动度、等启动参数的变化规律。
图4为导叶旋转示意图,导叶几何调节后,直接改变导叶出口几何出气角。通过简单的几何关系易得到旋转后导叶几何出气角α1z为α1与旋转角seita之和,其中seita为正值时表示导叶开大,seita为负值则相反,在忽略落后角的条件下旋转后导叶绝对出口气流角与动叶相对出口气流角如下式所示:
Figure GDA0003986670380000086
导叶旋转大幅度改变了导叶的流量特性曲线,并且由于导叶旋转后出口静压等气动参数的改变,动叶流量特性曲线也会有小幅度改变。流量的改变主要由两个因素决定,一是导叶出口气流角、二是导叶总静膨胀比,并且导叶出口气流角影响远大于导叶总静膨胀比,因此可以用公式(16)表示变几何涡轮导叶转角变化对流量的影响。
Figure GDA0003986670380000091
根据实际流量变化需求,并利用公式(16),求得变几何涡轮可转导叶转角变化范围。
对于变几何涡轮,变几何涡轮可调导叶间隙泄漏损失模型的修正公式为:
YTI=Kz(Ygap+Ymix) (17)
Figure GDA0003986670380000092
Figure GDA0003986670380000093
其中YTI为间隙泄漏损失,Kz为旋转轴修正系数,Ygap为间隙内损失,Ymix为泄漏流与主流掺混损失,C为叶片弦长,S为叶片栅距,H为叶片高度,μd为叶尖间隙流量系数,CL升力系数,αm为通过叶片的气流角平均值,KE为湍动能损失系数,KG为叶尖间隙损失系数;
旋转轴修正系数Kz为:
Figure GDA0003986670380000094
其中p00=0.8145、p10=-2.135、p01=16.55、p20=4.16、p11=-0.2775、p02=-81.45、p30=-2.36、p21=-30.58、p12=37.67、p03=142.9、p31=28.77、p22=-17.22、p13=-20.86、p04=-93.6。变量x为无量纲旋转轴位置,变量y为无量纲轴径,p00、p10、p01、p20、p11、p02、p30、p21、p12、p03、p31、p22、p13、p04均为拟合多项式的系数,其中字母p后的第一位数字表示变量x的次幂,字母p后的第二位数字表示变量y的次幂。
以某型变几何涡轮为例,给定折合转速
Figure GDA0003986670380000095
折合流量
Figure GDA0003986670380000096
轮缘功给定263.792(kW/(kg/s)),动静叶展弦比给定1.5,静叶叶栅稠度给定1.1、动叶稠度给定1.7,静叶最大厚度比弦长选定0.15,动叶最大厚度比弦长选定0.2,静叶尾缘厚度比喉部宽度选定0.12,动叶尾缘厚度比弦长给定0.2,静叶间隙高度比叶高选定1.8%,动叶间隙高度比叶高选定1%,旋转轴选择60%轴向弦长位置,轴径比弦长选择34%。以此为条件,研究不同转角工况下,反动度、轴向速比、载荷系数、流量系数的选取规律,且不同转角下边界条件保持相同。
从图6为流量系数为0.4,不同载荷系数(HT)下,不同转角工况平均效率与设计效率随反动度变化曲线图,总体上来说,反动度较小时,各工况下平均效率小于设计工况下效率,反之各工况下平均效率大于设计工况下效率,效率随反动度先增加后降低,存在最佳反动度,最佳反动度随载荷系数的增加而逐渐增大。总体来说,载荷系数为1时平均效率对应的最佳反动度比设计工况下的最佳反动度约小10%,随着载荷系数的增大,平均效率对应的最佳反动度比设计工况下的最佳反动度逐渐接近,变几何涡轮最佳反动度区域基本在0.25-0.6之间。
采用本发明设计的单级涡轮,在涡轮流量保持基本不变的前提下,其气动效率可提高0.7个百分点。
本发明也可用于多级变几何涡轮的高效气动设计。

Claims (5)

1.一种变几何涡轮一维气动设计方法,其特征是:
(1)给定初始几何参数与气动参数,包括动静叶展弦比、最大厚度比弦长、尾缘厚度比喉部宽度、轮缘功、进口总温总压、导叶进口气流角、载荷系数、流量系数、反动度、轴向速比、进出口中径比;
(2)根据涡轮一维几何参数求解动静叶轮毂比、叶高,并通过一维总静流函数法分析导叶转角转动后,涡轮级反动度、载荷系数、流量系数、轴向速比、气流角的变化;
(3)将导叶转动后,气动参数与几何参数输入到损失模型中,计算各种转角变化后动静叶速度损失系数,即可求出变几何涡轮各转角下涡轮效率。
2.根据权利要求1所述的一种变几何涡轮一维气动设计方法,其特征是:
(a)利用可转导叶转角与变几何涡轮流量关联式
Figure FDA0004125555510000011
预估所需的可转导叶旋转角,其中α1导叶未旋转的出口气流角,seita导叶旋转角,m1为变几何涡轮导叶转角变化前的流量,m1z为变几何涡轮导叶转角变化后的流量;
(b)将可转导叶旋转角分为n份,利用一维总静流函数法并综合可调导叶间隙泄漏损失模型分别计算各转角工况效率;
(c)基于步骤(b)计算变几何涡轮在不同转角下的平均效率;
(d)如果步骤(c)得到的性能指标不符合预定设计目标,则修改常规涡轮一维设计循环中的一维设计参数,并循环(a)-(c),直至涡轮气动性能达到预定设计目标。
3.根据权利要求1所述的一种变几何涡轮一维气动设计方法,其特征是:
对于变几何涡轮,一维总静流函数法的推导过程为:
一维总静流函数法分别包含了导叶流量m1和动叶流量m2随导叶总静膨胀比
Figure FDA0004125555510000012
的变化规律,将涡轮级通道看作一维喷管,导叶入口总压P0 *,导叶入口总温T0 *,动叶出口静压P1;导叶流量定义为:
m1=ρ1c1A1sinα1
其中ρ1为导叶通道内气体密度,c1为导叶出口绝对速度,A1为导叶出口截面积,α1为导叶绝对出口气流角;
根据理想气体状态方程P0 *=ρ0 *RgT0 *及多变过程
Figure FDA0004125555510000013
及T0 *=T1 *得公式:
Figure FDA0004125555510000021
Figure FDA0004125555510000022
其中ns,s为导叶气体多变指数,T1 *为导叶出口总温,ρ0 *为导叶进口密度,T1为导叶出口静温,Cp为气体定压比热容,h0 *为导叶进口滞止焓,h1为导叶出口焓值,k为绝热指数;
综合以上公式,可得出
Figure FDA0004125555510000023
即可得到导叶流量随导叶总静膨胀比变化规律;
通过下式建立导叶多变过程指数与速度损失系数
Figure FDA0004125555510000024
的关系:
Figure FDA0004125555510000025
根据连续方程及能量守恒方程,可得到λc1和动叶进口静温T1
Figure FDA0004125555510000026
T* 1=T* 0
T1=T* 1τ(λc1)
其中Rg为气体常数,λc1为导叶出口绝对速度因数,ε(λc1)为基于导叶出口绝对速度因数的密度气动函数,τ(λc1)为基于导叶出口绝对速度因数的温度气动函数;
根据动叶进口连续方程,可由λc1得到λω1,并根据λω1即可得出动叶进口相对总温及相对总压:
Figure FDA0004125555510000027
P* ω1=P1/π(λω1)
T* ω1=T1/τ(λω1)
其中λω1为基于动叶进口相对速度因数,P* ω1为动叶进口相对总压,T* ω1为动叶进口相对总温,Δπ(λω1)为基于动叶进口相对速度因数的压力气动函数,τ(λω1)为基于动叶进口相对速度因数的温度气动函数;
并且根据上述公式及气动函数可求出导叶出口速度因数
Figure FDA0004125555510000036
动叶进口静温T1,对于动叶而言,可以在相对坐标系下同样通过下述公式推导出动叶流量随导叶总静膨胀比,速度损失系数,出口相对气流角的变化规律:
m2=ρ2ω2A2sinβ2
其中ρ2为动叶出口气体密度,ω2为动叶出口相对速度,A2为动叶出口截面积;
根据理想气体状态方程及P1=ρ1RgT1及多变过程
Figure FDA0004125555510000031
可得
Figure FDA0004125555510000032
Figure FDA0004125555510000033
/>
Figure FDA0004125555510000034
其中ns为气体多变指数,P2为动叶出口静压,ns,r为动叶气体多变指数,ω2为动叶出口相对速度,hω1 *为动叶进口相对滞止焓,h2动叶出口焓值;
联立上述公式即可得到动叶流量随涡轮导叶总静膨胀比的变化规律:
Figure FDA0004125555510000035
其中β2为动叶几何出气角,其他参数均可由已知参数求得,因此可得到动叶流量随导叶出口总静膨胀比变化规律,动叶流量与导叶流量交点即为涡轮级工作点。
4.根据权利要求1所述的一种变几何涡轮一维气动设计方法,其特征是:对于变几何涡轮,变几何涡轮导叶转角变化与流量关联式的推导过程为:
导叶几何调节后,直接改变导叶出口几何出气角,通过几何关系得到旋转后导叶几何出气角α1z,α1z为α1与旋转角seita之和,其中seita为正值时表示导叶开大,seita为负值则相反,在忽略落后角的条件下旋转后导叶绝对出口气流角与动叶相对出口气流角如下式所示:
α1z=α1+seita
β2z=β2
对于变几何涡轮,变几何涡轮可调导叶间隙泄漏损失模型的修正公式为:
YTI=Kz(Ygap+Ymix)
Figure FDA0004125555510000041
Figure FDA0004125555510000042
其中YTI为间隙泄漏损失,Kz为旋转轴修正系数,Ygap为间隙内损失,Ymix为泄漏流与主流掺混损失,C为叶片弦长,S为叶片栅距,H为叶片高度,μd为叶尖间隙流量系数,CL升力系数,αm为通过叶片的气流角平均值,KE为湍动能损失系数,KG为叶尖间隙损失系数;
旋转轴修正系数Kz为:
Figure FDA0004125555510000043
其中变量x为无量纲旋转轴位置,变量y为无量纲轴径,p00、p10、p01、p20、p11、p02、p30、p21、p12、p03、p31、p22、p13、p04均为拟合多项式的系数,其中字母p后的第一位数字表示变量x的次幂,字母p后的第二位数字表示变量y的次幂。
5.根据权利要求1所述的一种变几何涡轮一维气动设计方法,其特征是:变几何涡轮导叶出口气流角的范围在15°-25°度之间,变几何涡轮最佳反动度区域在0.25-0.6之间,流量系数在0.3-0.5之间,轴向速比选择范围为1-1.4之间。
CN202110772097.8A 2021-07-08 2021-07-08 一种变几何涡轮一维气动设计方法 Active CN113513373B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110772097.8A CN113513373B (zh) 2021-07-08 2021-07-08 一种变几何涡轮一维气动设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110772097.8A CN113513373B (zh) 2021-07-08 2021-07-08 一种变几何涡轮一维气动设计方法

Publications (2)

Publication Number Publication Date
CN113513373A CN113513373A (zh) 2021-10-19
CN113513373B true CN113513373B (zh) 2023-05-30

Family

ID=78067104

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110772097.8A Active CN113513373B (zh) 2021-07-08 2021-07-08 一种变几何涡轮一维气动设计方法

Country Status (1)

Country Link
CN (1) CN113513373B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114017201B (zh) * 2021-11-16 2023-04-07 中国航发沈阳发动机研究所 一种航空燃气涡扇发动机喷口控制方法及装置
CN114491417B (zh) * 2022-04-07 2022-07-26 中国航发四川燃气涡轮研究院 一种基于cdfs模态变化性能的一维输入修正方法
CN115186443B (zh) * 2022-06-15 2023-06-13 中国船舶重工集团公司第七0三研究所 一种倒置浴盆载荷船用发电型燃气轮机多级动力涡轮气动设计方法
CN115238370B (zh) * 2022-06-15 2023-06-27 中国船舶重工集团公司第七0三研究所 一种船用燃气轮机单级涡轮通流尺度快速评估方法
CN115186441B (zh) * 2022-06-15 2023-07-07 中国船舶重工集团公司第七0三研究所 一种递增载荷船用发电型燃气轮机多级动力涡轮气动设计方法
CN115186442B (zh) * 2022-06-15 2023-08-08 中国船舶重工集团公司第七0三研究所 一种递减载荷船用发电型燃气轮机多级动力涡轮气动设计方法
CN115186440B (zh) * 2022-06-15 2023-06-27 中国船舶重工集团公司第七0三研究所 一种船用发电型燃气轮机双级高速动力涡轮气动设计方法
CN116702511B (zh) * 2023-08-01 2023-10-31 中国航发四川燃气涡轮研究院 一种可调导叶落后角的计算方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002256810A (ja) * 2001-03-05 2002-09-11 Toshiba Corp 軸流タービン
EP1900908A2 (en) * 2006-09-12 2008-03-19 Iveco Motorenforschung AG Variable geometry turbine
WO2013080795A1 (ja) * 2011-11-30 2013-06-06 三菱重工業株式会社 ラジアルタービン
CN105626158A (zh) * 2016-03-03 2016-06-01 哈尔滨工程大学 一种带有动叶片前部消涡孔结构的变几何涡轮
CN108729958A (zh) * 2018-04-24 2018-11-02 哈尔滨工程大学 一种具有低稠度无升力叶型的可调静叶的反转变几何涡轮
CN112528474A (zh) * 2020-11-30 2021-03-19 北京动力机械研究所 一种无导叶对转涡轮一维气动设计方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2184882A1 (en) * 1995-09-08 1997-03-09 Hideomi Harada Turbomachinery with variable-angle flow guiding vanes
US5794432A (en) * 1996-08-27 1998-08-18 Diversitech, Inc. Variable pressure and variable air flow turbofan engines
EP1881155A1 (en) * 2006-07-21 2008-01-23 Ansaldo Energia S.P.A. Apparatus for designing turbines and turbine blades
JP6866019B2 (ja) * 2014-06-24 2021-04-28 コンセプツ エヌアールイーシー,エルエルシー ターボ機械の流動制御構造及びその設計方法
US9890700B2 (en) * 2014-11-21 2018-02-13 Ford Global Technologies, Llc Systems and methods for a variable geometry turbine nozzle
CN105781626B (zh) * 2016-03-03 2017-05-24 哈尔滨工程大学 一种大子午扩张涡轮的三维正交非定常设计方法
CN205422837U (zh) * 2016-03-03 2016-08-03 哈尔滨工程大学 一种带有动叶片前部消涡孔结构的变几何涡轮
CN109779972B (zh) * 2019-01-25 2021-01-12 杭州汽轮动力集团有限公司 一种兼顾低雷诺数工况性能的压气机二维叶型优化方法
CN111680372B (zh) * 2020-06-10 2023-06-20 大连海事大学 考虑自然预旋时离心风机叶轮做功能力的一维计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002256810A (ja) * 2001-03-05 2002-09-11 Toshiba Corp 軸流タービン
EP1900908A2 (en) * 2006-09-12 2008-03-19 Iveco Motorenforschung AG Variable geometry turbine
WO2013080795A1 (ja) * 2011-11-30 2013-06-06 三菱重工業株式会社 ラジアルタービン
CN105626158A (zh) * 2016-03-03 2016-06-01 哈尔滨工程大学 一种带有动叶片前部消涡孔结构的变几何涡轮
CN108729958A (zh) * 2018-04-24 2018-11-02 哈尔滨工程大学 一种具有低稠度无升力叶型的可调静叶的反转变几何涡轮
CN112528474A (zh) * 2020-11-30 2021-03-19 北京动力机械研究所 一种无导叶对转涡轮一维气动设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Jie Gao等.Aerothermal characteristics of a transonic tip flow in a turbine cascade with tip clearance variations.Applied Thermal Engineering.2016,(第107期),271-283. *
吴中野 ; 方祥军 ; .变几何低压涡轮气动性能研究.推进技术.2018,(第02期),全文. *

Also Published As

Publication number Publication date
CN113513373A (zh) 2021-10-19

Similar Documents

Publication Publication Date Title
CN113513373B (zh) 一种变几何涡轮一维气动设计方法
CN109815624B (zh) 一种考虑进气总压畸变影响的压气机稳定边界判断方法
CN110701086B (zh) 压气机全工况性能预测方法
CN112685968B (zh) 基于空间负荷定制思想的轴流压气机气动设计方法
CN114444196B (zh) 一种自适应缝隙调节可变弯度导叶设计方法及其装置
CN114444331B (zh) 一种多级轴流压气机的级特性匹配方法
CN112464392B (zh) 一种压气机防喘扩稳的多列可转导/静叶联合调控转角规律设计方法
CN109635512B (zh) 一种基于修正控制方程的离心叶轮进口设计方法
CN111435399A (zh) 风扇组件的造型方法
CN106939902B (zh) 节能型直壁前、后盘变曲率曲线元素三元叶轮及采用该叶轮的离心式风机
CN111832189A (zh) 一种离心式蒸汽压缩机选型方法
CN109505790B (zh) 高负荷高通流能力的轴流风机
CN117195760A (zh) 基于径向掺混的轴流风扇或压气机子午面通流计算方法
CN114444415B (zh) 径流涡轮定子喷嘴导流叶片设计优化方法和系统
CN115221650A (zh) 一种基于径向速比寻优的轴流透平三维叶片设计方法
CN112052633B (zh) 一种小容积流量下汽轮机末级叶片脱流高度的确定方法
CN109711047B (zh) 一种氦气压气机s2流面气动的设计方法
CN114117665B (zh) 一种s2流面框架下轴流压气机经验模型标定方法
CN115841085A (zh) 基于无量纲数的多级轴流跨声速涡轮特性一维预估方法
CN113137282B (zh) 一种反动式透平级组连接结构及其结构参数设计方法
CN111911423B (zh) 一种基于piv的离心泵叶轮极限外径测量及优化方法
CN118194472A (zh) 一种航空发动机多涵道风扇多点协调设计方法
CN115270474A (zh) 一种基于变比热计算的轴流涡轮一维性能正问题评估方法
CN115270341A (zh) 基于快速流量匹配计算的轴流多级涡轮一维特性计算方法
CN114707264A (zh) 一种变循环发动机压缩部件特性解耦方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant