CN114676533B - 一种与空气系统耦合的涡轮动叶冷却敏感性分析方法 - Google Patents
一种与空气系统耦合的涡轮动叶冷却敏感性分析方法 Download PDFInfo
- Publication number
- CN114676533B CN114676533B CN202210578070.XA CN202210578070A CN114676533B CN 114676533 B CN114676533 B CN 114676533B CN 202210578070 A CN202210578070 A CN 202210578070A CN 114676533 B CN114676533 B CN 114676533B
- Authority
- CN
- China
- Prior art keywords
- blade
- throttling element
- temperature
- cooling
- turbine
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明提供了一种与空气系统耦合的涡轮动叶冷却敏感性分析方法,包括如下步骤:S1:判断涡轮动叶冷却气流量和气膜孔逆流裕度是否满足要求;S2:对气膜孔逆流裕度进行敏感性分析;S3:判断温度场是否满足要求;S4:对不满足材料长期使用温度处的温度进行敏感性分析;S5:计算涡轮动叶在S3获取的温度场中的强度;S6:对不满足强度要求处的温度进行敏感性分析。通过空气系统与叶片冷却耦合敏感性分析,获得冷却效果与空气系统节流元件及叶片冷却结构的几何尺寸之间的敏感性关系,得到影响叶片冷却的关键参数,通过调节关键参数,快速而精准的解决设计过程中的不满足因素,从而解决了涡轮动叶冷却设计。
Description
技术领域
本发明属于航空发动机技术领域,具体涉及一种与空气系统耦合的涡轮动叶冷却敏感性分析方法。
背景技术
为了实现航空发动机与日俱增的推力需求,燃烧室出口温度不断提高,导致热端部件的工作环境变得非常恶劣。由于主流道压力较高,高压涡轮动叶冷却一般采用温度较高的压气机出口引气,冷却效果很大程度上受空气系统布局及叶片冷却结构的影响。尤其是在高马赫数航空发动机中,随着冷气温度进一步提高,冷却介质的冷却品质进一步下降,使得热端部件的冷却需求和供气流量之间的矛盾日益凸显。
目前涡轮动叶冷却效果优化多集中在叶片冷却结构的改进,可调元件有限,很难达到预期的效果,由于相同的引气压力下,动叶供气系统的结构布局对进入叶片的冷气压力和流量存在很大影响,进而影响叶片的冷却效果和叶片表面气膜孔逆流裕度,因此若能将空气系统元件与叶片冷却设计元件耦合进行调整,更能接近甚至超越预期冷却效果。然而空气系统和叶片冷却节流元件数量多、流动复杂,想要精准的找出最有效的调节元件很难,同时需要花费大量的时间和人力进行反复迭代。
因此,需要提供一种能将空气系统元件和叶片冷却元件耦合分析敏感性的方法以解决上述问题。
发明内容
为了解决上述问题,本发明的目的在于提供一种与空气系统耦合的涡轮动叶冷却敏感性分析方法,所述分析方法包括如下步骤:
S1:获取空气系统节流元件几何参数和涡轮动叶冷却节流元件的几何参数,对前述参数进行耦合计算,得到涡轮动叶内流计算结果,判断涡轮动叶冷却气膜孔逆流裕度是否满足要求,若不满足,则进行S2,若满足,则进行S3;
S2:对气膜孔逆流裕度进行敏感性分析,得到气膜孔逆流裕度与相关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,以冷气流量为限制条件,将满足要求的最小逆流裕度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整空气系统和叶片冷却节流元件几何参数,调整后返回S1;
S3:根据S1获取的内流计算结果进行热分析计算,判断温度场是否满足要求,若不满足,则进行S4,若满足,则进行S5;
S4:对不满足材料长期使用温度处的温度进行敏感性分析,得到不满足材料长期使用温度处的温度与相关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,以冷气流量为限制条件,将材料长期使用温度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整耦合空气系统和叶片冷却节流元件几何参数,调整后返回S1:
S5:计算涡轮动叶在S3获取的温度场中的强度,若不满足,则进行S6,若满足,则涡轮动叶冷却结构满足要求,敏感性分析结束;
S6:对不满足强度要求处的温度进行敏感性分析,得到该温度与关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,计算满足强度的最高温度,以冷气流量为限制条件,将计算的最高温度作目标进行优化,获得满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整空气系统和叶片冷却节流元件几何参数,调整后返回S1。
本发明所提供的与空气系统耦合的涡轮动叶冷却敏感性分析方法,还具有这样的特征,关键节流元件为敏感性关系曲线斜率绝对值大于0.1的节流元件。
本发明所提供的与空气系统耦合的涡轮动叶冷却敏感性分析方法,还具有这样的特征,所述S4中的长期使用温度为根据发动机工作包线内最恶劣的状态工作时间确定的材料许用温度。
本发明所提供的与空气系统耦合的涡轮动叶冷却敏感性分析方法,还具有这样的特征,所述分析方法利用仿真模型进行,所述仿真模型为根据初步的流路布局及叶片冷却结构建立的空气系统与叶片冷却耦合仿真模型。
本发明所提供的与空气系统耦合的涡轮动叶冷却敏感性分析方法,还具有这样的特征,所述空气系统节流元件包括与涡轮动叶冷却相关的预旋喷嘴、引气孔、篦齿、缝隙和管路;所述涡轮动叶冷却节流元件包括气膜孔、冲击孔、扰流肋、尾缝、扰流柱和榫头进气孔。
本发明所提供的与空气系统耦合的涡轮动叶冷却敏感性分析方法,还具有这样的特征,所述S1中的耦合计算包括如下步骤:
A、读取空气系统计算输入文件,给定向叶片供气的节流元件初始流通面积SA0;
B、计算空气系统节流元件进出口压力、温度、流量;
C、通过空气系统计算结果获得涡轮动叶冷却的进口压力P0、进口温度T0和进口流量G0;
D、对进口压力P0和进口温度T0进行涡轮动叶内流计算,获得涡轮动叶内流计算结果,读取涡轮动叶冷却引气量G1;
E、给定流量计算残差要求值ε,进行判断,
若|G1-G0|/G0≥ε,则将修正后的面积SA1=SA0+(G1-G0)/G0*SA0设置为空气系统文件中向叶片供气的节流面积,返回步骤A;
若|G1-G0|/G0<ε ,则结束耦合计算。
本发明所提供的与空气系统耦合的涡轮动叶冷却敏感性分析方法,还具有这样的特征,所述敏感性分析包括如下步骤:
a、选取空气系统节流元件和涡轮动叶的所有节流元件的几何参数;
b、设置每个节流元件几何参数基于初始值的变化范围;
c、开展试验设计分析,获得多个由节流元件几何参数组合的样本;
d、针对c获得的每个样本,进行空气系统与叶片流动耦合计算分析;
e、绘制主因素图获得输出变量与节流元件之间的敏感性关系;
f、根据敏感性关系选择影响输出变量的参数。
本发明所提供的与空气系统耦合的涡轮动叶冷却敏感性分析方法,还具有这样的特征,所述S3包括如下步骤:
S3.1:通过求解流过涡轮动叶型面的粘性流控制微分方程,获得叶片燃气边的外换热系数;
S3.2:根据内流计算结果对S3.1获得的外换热系数计算结果进行气膜修正;
S3.3:将经过气膜修正后的燃气边换热系数分布、温度分布和压力分布作为叶片内流计算气膜出口边界条件,重新进行内流计算;
S3.4:根据气膜修正后的外换热系数、外表面燃气温度分布和内流计算结果,采用线性差值的方法,分别获得叶片内外表面的换热系数和温度分布云图;
S3.5:将得到的叶片内边界的换热系数、外边界的换热系数和温度分布加载在叶片模型上,最终得到叶片的温度场结果;
S3.6:判断温度场结果是否满足要求。
本发明所提供的与空气系统耦合的涡轮动叶冷却敏感性分析方法,还具有这样的特征,所述S1中的内流计算结果包括涡轮动叶内部各腔流量分配数、换热系数、压力分布和温度分布。
有益效果
本发明所提供的与空气系统耦合的涡轮动叶冷却结构敏感性分析方法通过与空气系统耦合设计,避免了人为迭代带来的计算误差,同时提高了设计效率;通过空气系统与叶片冷却耦合敏感性分析,获得冷却效果与空气系统节流元件及叶片冷却结构的几何尺寸之间的敏感性关系,得到影响叶片冷却的关键参数,通过调节关键参数,快速而精准的解决设计过程中的不满足因素,从而有效的解决了涡轮动叶冷却设计分析问题。
附图说明
图1为本发明实施例所提供的与空气系统耦合的涡轮动叶冷却结构敏感性分析方法的流程图;
图2为本发明实施例所提供的分析方法的耦合计算模型;
图3为本发明实施例所提供的分析方法中的逆流裕度与相关联节流元件参数的敏感性关系图;
图4为本发明实施例提供的分析方法中的热分析结果图;
图5为本发明实施例提供的分析方法中的温度与相关联节流元件参数的敏感性关系;
图6为本发明实施例提供的分析方法中的温度场结果图;
图7为本发明实施例所提供的分析方法中的温度与相关联节流元件参数的敏感性关系图。
具体实施方式
下面结合附图与实施例对本发明作进一步的详细说明,但应当说明的是,这些实施方式并非对本发明的限制,本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代,均属于本发明的保护范围之内。
在本发明实施例的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明创造和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明创造的限制。
此外,术语“第一”、“第二”、“第三”等仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”等的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明创造的描述中,除非另有说明,“多个”的含义是两个或两个以上。
术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以通过具体情况理解上述术语在本发明创造中的具体含义。
本发明实施例提供了一种与空气系统耦合的涡轮动叶冷却敏感性分析方法,如图1所示,所述分析方法包括如下步骤:
S1:获取空气系统节流元件几何参数和涡轮动叶冷却节流元件的几何参数,对前述参数进行耦合计算,得到涡轮动叶内流计算结果,判断涡轮动叶冷却气膜孔逆流裕度是否满足要求,若不满足,则进行S2,若满足,则进行S3;其中,气膜孔逆流裕度为气膜孔进出口压力之比;
S2:对气膜孔逆流裕度进行敏感性分析,得到气膜孔逆流裕度与相关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,以冷气流量为限制条件,将满足要求的最小逆流裕度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整空气系统和叶片冷却节流元件几何参数,调整后返回S1,关键节流元件为敏感性关系曲线斜率绝对值大于0.1的节流元件;
S3:根据S1获取的内流计算结果进行热分析计算,判断温度场是否满足要求,若不满足,则进行S4,若满足,则进行S5;
S4:对不满足材料长期使用温度处的温度进行敏感性分析,得到不满足材料长期使用温度处的温度与相关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,以冷气流量为限制条件,将材料长期使用温度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整空气系统和叶片冷却节流元件几何参数,调整后返回S1,关键节流元件为敏感性关系曲线斜率绝对值大于0.1的节流元件;
S5:计算涡轮动叶在S3获取的温度场中的强度,若不满足,则进行S6,若满足,则涡轮动叶冷却结构满足要求,敏感性分析结束;
S6:对不满足强度要求处的温度进行敏感性分析,得到该温度与关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,计算满足强度的最高温度,以冷气流量为限制条件,将计算的最高温度作目标进行优化,获得满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整空气系统和叶片冷却节流元件几何参数,调整后返回S1,关键节流元件为敏感性关系曲线斜率绝对值大于0.1的节流元件。其中,满足要求的关键节流元件几何参数最优组合是程序根据优化目标在给定的样本变化范围内自动匹配的最优组合。
在部分实施例中,所述S4中的长期使用温度为根据发动机工作包线内最恶劣的状态工作时间确定的材料许用温度。长期使用温度指发动机在该温度下和低于这个温度时,材料在发动机全寿命周期内,都能安全可靠的工作,超过这个温度,材料的寿命,可靠性,稳定性等参数可能不满足使用要求,该温度一般在材料手册上会有限定。
在部分实施例中,所述分析方法利用仿真模型进行,所述仿真模型为根据初步的流路布局及叶片冷却结构建立的空气系统与叶片冷却耦合仿真模型。
在部分实施例中,所述空气系统节流元件包括与涡轮动叶冷却相关的预旋喷嘴、引气孔、篦齿、缝隙和管路;所述涡轮动叶冷却节流元件包括气膜孔、冲击孔、扰流肋、尾缝、扰流柱和榫头进气孔。
在部分实施例中,所述S1中的耦合计算包括如下步骤:
A、读取空气系统计算输入文件,给定向叶片供气的节流元件初始流通面积SA0;
B、计算空气系统节流元件进出口压力、温度、流量;
C、通过空气系统计算结果获得涡轮动叶冷却的进口压力P0、进口温度T0和进口流量G0;
D、对进口压力P0和进口温度T0进行涡轮动叶内流计算,获得涡轮动叶内流计算结果,读取涡轮动叶冷却引气量G1;
E、给定流量计算残差要求值ε,进行判断,
若|G1-G0|/G0≥ε,则将修正后的面积SA1=SA0+(G1-G0)/G0*SA0设置为空气系统文件中向叶片供气的节流面积,返回步骤A;
若|G1-G0|/G0<ε ,则结束耦合计算。
在部分实施例中,所述敏感性分析包括如下步骤:
a、选取空气系统节流元件和涡轮动叶的所有节流元件的几何参数;
b、设置每个节流元件几何参数基于初始值的变化范围;
c、开展试验设计分析,获得多个由节流元件几何参数组合的样本;
d、针对c获得的每个样本,进行空气系统与叶片流动耦合计算分析;
e、绘制主因素图获得输出变量与节流元件之间的敏感性关系;
f、根据敏感性关系选择影响输出变量的参数。
在部分实施例中,所述S3包括如下步骤:
S3.1:通过求解流过涡轮动叶型面的粘性流控制微分方程,获得叶片燃气边的外换热系数;
S3.2:根据内流计算结果对S3.1获得的外换热系数计算结果进行气膜修正;
S3.3:将经过气膜修正后的燃气边换热系数分布、温度分布和压力分布作为叶片内流计算气膜出口边界条件,重新进行内流计算;
S3.4:根据气膜修正后的外换热系数、外表面燃气温度分布和内流计算结果,采用线性差值的方法,分别获得叶片内外表面的换热系数和温度分布云图;
S3.5:将得到的叶片内边界的换热系数、外边界的换热系数和温度分布加载在叶片模型上,最终得到叶片的温度场结果;
S3.6:判断温度场结果是否满足要求。
在部分实施例中,所述S1中的内流计算结果包括涡轮动叶内部各腔流量分配数、换热系数、压力分布和温度分布。
在部分实施例中,根据压气机进口流量60kg/s、压气机出口压力450kPa、压气机出口温度980K、涡轮气动数据S1和S2等输入参数进行敏感性分析,步骤如下:
步骤1:初步提出叶片供气系统结构及叶片冷却结构;
步骤2:针对步骤1提出的供气系统结构和叶片冷却结构建立空气系统与叶片冷却耦合计算模型,如图2所示;
步骤3:根据步骤2建立的模型,进行耦合计算,计算结果如表1:
表1
从上表可知,冷气量为3.462%,满足设计要求,最小逆流裕度为0.998,出现燃气倒灌,不满足设计要求;
步骤4:对叶片气膜孔最小逆流裕度进行敏感性分析,得到逆流裕度与相关联节流元件参数的敏感性关系,如图3所示,其中,图3中的APRE为空气系统预旋系统接收孔面积;Ayx为预旋喷嘴的面积;B-ANI1为第一叶片腔进口面积;B-AIN2为第二叶片腔进口面积。根据敏感性关系选取关键节流元件,以冷气流量为限制条件,将符合要求的最小逆流裕度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合;
根据最优组合参数调整空气系统及叶片冷却耦合输入模型,进行耦合计算,得到表2的数据,
表2
计算结果表明在总冷气流量不变的情况下,尾缝的冷气流量降低,其余部分的冷气流量有不同程度的增加;对比各部分最小逆流裕度的变化,可以看出逆流裕度有不同程度的增加,特别是前缘区域的逆流裕度,由原来的0.998提升到1.084,大幅度改善了涡轮叶片前缘的气膜出流情况,增强了涡轮叶片工作可靠性;
步骤5:对满足要求的叶片冷却设计结果进行热分析计算,得到如图4所示的温度场结果,可以看出叶片最高温度在前缘处,为1412K,最低温度为1114K,最高温度超过材料许用温度1393K,不满足要求;
此时,对叶片前缘温度进行敏感性分析,得到如图5所示的温度与相关联节流元件参数的敏感性关系图,其中,图5中的B-DFP1为叶片根截面气膜孔径;B-DFP2为叶片中截面气膜孔径;B-DFP3为叶片尖截面气膜孔径;B-DOUT1为第一叶片腔叶尖出气孔直径。根据敏感性关系选取关键节流元件,以冷气流量为限制条件,将材料许用温度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合;
根据最优组合参数调整空气系统及叶片冷却耦合输入模型,进行耦合计算、热分析计算,调整后的温度场结果如图6所示,表明叶片的最低温度为1089K,最高温度在前缘处,为1388K,未超过材料的许用温度,满足设计要求,可进行强度计算;
步骤6:进行强度计算,计算结果如表3所示
表3
计算结果表明,涡轮工作叶片在该工况载荷作用下根截面最小持久储备不满足1.3的强度设计要求;
对叶片根截面处温度进行敏感性分析,得到温度与相关联节流元件参数的敏感性关系图,如图7所示,其中,图7中A118-9为卸荷腔放气面积;APRE为预旋系统接受孔面积;Ayx为预旋喷嘴的面积。根据敏感性关系选取关键节流元件,计算根截面处满足要求的最高温度,以冷气流量为限制条件,将该最高温度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合;
根据最优组合参数调整空气系统及叶片冷却耦合输入模型,进行耦合计算、热分析计算、强度计算,结果如表4所示,
表4
计算结果表明,调整后根截面的最小持久储备为1.794,满足强度设计要求,则完成与空气系统耦合的涡轮动叶冷却设计与分析。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变型,这些改进和变型也应视为本发明的保护范围。
Claims (9)
1.一种与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,所述分析方法包括如下步骤:
S1:获取空气系统节流元件几何参数和涡轮动叶冷却节流元件的几何参数,对前述参数进行耦合计算,得到涡轮动叶内流计算结果,判断涡轮动叶冷却气膜孔逆流裕度是否满足要求,若不满足,则进行S2,若满足,则进行S3;
S2:对气膜孔逆流裕度进行敏感性分析,得到气膜孔逆流裕度与相关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,以冷气流量为限制条件,将满足要求的最小逆流裕度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整空气系统和叶片冷却节流元件几何参数,调整后返回S1;
S3:根据S1获取的内流计算结果进行热分析计算,判断温度场是否满足要求,若不满足,则进行S4,若满足,则进行S5;
S4:对不满足材料长期使用温度处的温度进行敏感性分析,得到不满足材料长期使用温度处的温度与相关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,以冷气流量为限制条件,将材料长期使用温度作为目标进行优化,得到满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整空气系统和叶片冷却节流元件几何参数,调整后返回S1:
S5:计算涡轮动叶在S3获取的温度场中的强度,若不满足,则进行S6,若满足,则涡轮动叶冷却结构满足要求,敏感性分析结束;
S6:对不满足强度要求处的温度进行敏感性分析,得到该温度与关联节流元件参数的敏感性关系,根据敏感性关系选取关键节流元件,计算满足强度的最高温度,以冷气流量为限制条件,将计算的最高温度作目标进行优化,获得满足要求的关键节流元件几何参数最优组合,根据最优组合参数调整空气系统和叶片冷却节流元件几何参数,调整后返回S1。
2.根据权利要求1所述的与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,关键节流元件为敏感性关系曲线斜率绝对值大于0.1的节流元件。
3.根据权利要求1所述的与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,所述S4中的长期使用温度为根据发动机工作包线内最恶劣的状态工作时间确定的材料许用温度。
4.根据权利要求1所述的与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,所述分析方法利用仿真模型进行,所述仿真模型为根据初步的流路布局及叶片冷却结构建立的空气系统与叶片冷却耦合仿真模型。
5.根据权利要求1所述的与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,所述空气系统节流元件包括与涡轮动叶冷却相关的预旋喷嘴、引气孔、篦齿、缝隙和管路;所述涡轮动叶冷却节流元件包括气膜孔、冲击孔、扰流肋、尾缝、扰流柱和榫头进气孔。
6.根据权利要求1所述的与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,所述S1中的耦合计算包括如下步骤:
A、读取空气系统计算输入文件,给定向叶片供气的节流元件初始流通面积SA0;
B、计算空气系统节流元件进出口压力、温度、流量;
C、通过空气系统计算结果获得涡轮动叶冷却的进口压力P0、进口温度T0和进口流量G0;
D、对进口压力P0和进口温度T0进行涡轮动叶内流计算,进行涡轮动叶内流计算,获得涡轮动叶内流计算结果,读取涡轮动叶冷却引气量G1;
E、给定流量计算残差要求值ε,进行判断,
若|G1-G0|/G0≥ε,则将修正后的面积SA1=SA0+(G1-G0)/G0*SA0设置为空气系统文件中向叶片供气的节流面积,返回步骤A;
若|G1-G0|/G0<ε ,则结束耦合计算。
7.根据权利要求1所述的与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,所述敏感性分析包括如下步骤:
a、选取空气系统节流元件和涡轮动叶的所有节流元件的几何参数;
b、设置每个节流元件几何参数基于初始值的变化范围;
c、开展试验设计分析,获得多个由节流元件几何参数组合的样本;
d、针对c获得的每个样本,进行空气系统与叶片流动耦合计算分析;
e、绘制主因素图获得输出变量与节流元件之间的敏感性关系;
f、根据敏感性关系选择影响输出变量的参数。
8.根据权利要求1所述的与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,所述S3包括如下步骤:
S3.1:通过求解流过涡轮动叶型面的粘性流控制微分方程,获得叶片燃气边的外换热系数;
S3.2:根据内流计算结果对S3.1获得的外换热系数计算结果进行气膜修正;
S3.3:将经过气膜修正后的燃气边换热系数分布、温度分布和压力分布作为叶片内流计算气膜出口边界条件,重新进行内流计算;
S3.4:根据气膜修正后的外换热系数、外表面燃气温度分布和内流计算结果,采用线性差值的方法,分别获得叶片内外表面的换热系数和温度分布云图;
S3.5:将得到的叶片内边界的换热系数、外边界的换热系数和温度分布加载在叶片模型上,最终得到叶片的温度场结果;
S3.6:判断温度场结果是否满足要求。
9.根据权利要求1所述的与空气系统耦合的涡轮动叶冷却敏感性分析方法,其特征在于,所述S1中的内流计算结果包括涡轮动叶内部各腔流量分配数、换热系数、压力分布和温度分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210578070.XA CN114676533B (zh) | 2022-05-26 | 2022-05-26 | 一种与空气系统耦合的涡轮动叶冷却敏感性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210578070.XA CN114676533B (zh) | 2022-05-26 | 2022-05-26 | 一种与空气系统耦合的涡轮动叶冷却敏感性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114676533A CN114676533A (zh) | 2022-06-28 |
CN114676533B true CN114676533B (zh) | 2022-08-02 |
Family
ID=82079916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210578070.XA Active CN114676533B (zh) | 2022-05-26 | 2022-05-26 | 一种与空气系统耦合的涡轮动叶冷却敏感性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114676533B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115307811B (zh) * | 2022-10-10 | 2023-03-24 | 中国航发四川燃气涡轮研究院 | 基于壁面静压的叶尖泄漏流动测试方法 |
CN115391960B (zh) * | 2022-10-27 | 2023-03-07 | 中国航发四川燃气涡轮研究院 | 考虑分散系数和多轴载荷的主轴复合疲劳寿命分析方法 |
CN115788598B (zh) * | 2023-02-10 | 2023-06-30 | 中国航发四川燃气涡轮研究院 | 一种涡轮叶片气膜孔参数化控制与设计方法 |
CN116756882B (zh) * | 2023-08-23 | 2023-10-31 | 中国航发四川燃气涡轮研究院 | 基于几何参数敏感性分析的涡轮叶片冷却结构设计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0742347A2 (en) * | 1995-05-10 | 1996-11-13 | Allison Engine Company, Inc. | Turbine blade cooling |
CN106837431A (zh) * | 2016-12-26 | 2017-06-13 | 北京航空航天大学 | 一种用于提高航空燃气涡轮发动机过渡过程涡轮叶片气膜冷却气流品质的方法 |
CN107977523A (zh) * | 2017-12-14 | 2018-05-01 | 中国航发沈阳发动机研究所 | 一种面向参数化模型的涡轮气冷叶片设计方法 |
CN114151139A (zh) * | 2021-10-20 | 2022-03-08 | 中国航发四川燃气涡轮研究院 | 一种采用渗透模型模拟涡轮叶片表面气膜孔冷气层流动的方法 |
CN114282323A (zh) * | 2021-12-27 | 2022-04-05 | 北京航空航天大学 | 一种涡轮叶片层板冷却结构的流量分配预测方法 |
-
2022
- 2022-05-26 CN CN202210578070.XA patent/CN114676533B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0742347A2 (en) * | 1995-05-10 | 1996-11-13 | Allison Engine Company, Inc. | Turbine blade cooling |
CN106837431A (zh) * | 2016-12-26 | 2017-06-13 | 北京航空航天大学 | 一种用于提高航空燃气涡轮发动机过渡过程涡轮叶片气膜冷却气流品质的方法 |
CN107977523A (zh) * | 2017-12-14 | 2018-05-01 | 中国航发沈阳发动机研究所 | 一种面向参数化模型的涡轮气冷叶片设计方法 |
CN114151139A (zh) * | 2021-10-20 | 2022-03-08 | 中国航发四川燃气涡轮研究院 | 一种采用渗透模型模拟涡轮叶片表面气膜孔冷气层流动的方法 |
CN114282323A (zh) * | 2021-12-27 | 2022-04-05 | 北京航空航天大学 | 一种涡轮叶片层板冷却结构的流量分配预测方法 |
Non-Patent Citations (4)
Title |
---|
气膜冷却孔对涡轮发动机非定常流动的影响;张玲等;《拖拉机与农用运输车》;20160215(第01期);第37-41页 * |
浅谈涡轮转子叶片冷却技术;王旭伟;《山东工业技术》;20181210(第01期);第11页 * |
燃气轮机燃烧室网络式流量分配设计;吴晶峰等;《中国民航大学学报》;20160815(第04期);第53-58页 * |
航空发动机空气系统和热分析的耦合计算与试验验证;梁津华、赵维维、徐连强;《燃气涡轮试验与研究》;20190630;第24-28页、52页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114676533A (zh) | 2022-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114676533B (zh) | 一种与空气系统耦合的涡轮动叶冷却敏感性分析方法 | |
CN111914362B (zh) | 一种研发阶段涡扇发动机模型自适应方法 | |
US9255492B2 (en) | Gas turbine engine having a multi-variable closed loop controller for regulating tip clearance | |
CN112594064B (zh) | 一种基于轴流压气机级间测量参数的s2流场诊断方法 | |
CN105512429A (zh) | 一种三轴燃气轮机的总体方案计算方法 | |
CN110647052B (zh) | 一种变循环发动机模式切换自适应身份证模型构建方法 | |
Wellborn et al. | Redesign of a 12-stage axial-flow compressor using multistage CFD | |
US11210435B2 (en) | Method of designing and producing a turbine | |
Hu et al. | Thermodynamics investigation and spike-stall identification based on energy loss of centrifugal compressor | |
CN117195760A (zh) | 基于径向掺混的轴流风扇或压气机子午面通流计算方法 | |
CN116204983A (zh) | 一种航空发动机涡轮盘腔实时气温的计算方法 | |
Chana et al. | An investigation on turbine tip and shroud heat transfer | |
CN107849925A (zh) | 具有冷却流量调节特征的涡轮机定子叶片和/或涡轮机转子叶片以及适应叶片的对应方法 | |
Elfner et al. | Comparative experimental investigation of leading edge cooling concepts of turbine rotor blades | |
CN109738193A (zh) | 燃气轮机燃烧室试验用空冷式测量段结构 | |
CN111783306B (zh) | 涡轮叶片冷气量与冷效特性分析方法 | |
CN114720145A (zh) | 一种带整流叶片的低压涡轮性能试验方法 | |
Parvizinia et al. | Numerical and Experimental Investigations into the Aerodynamic performance of a supersonic turbine blade profile | |
Yao et al. | Numerical study of axial turbine flow with variable gas property | |
Rehman et al. | Non-Axisymmetric Endwall Profiling of a Stator Row in the Presence of the Rotor in a High Pressure Turbine | |
CN116611174B (zh) | 一种发动机涡轮导向器流量计算仿真模型的构建方法 | |
Kulkarni | Development of a Methodology to Estimate Aero-Performance and Aero-Operability Limits of a Multistage Axial Flow Compressor for Use in Preliminary Design | |
CN113137282B (zh) | 一种反动式透平级组连接结构及其结构参数设计方法 | |
Maksiuta et al. | Deviation Angle in a Turbine Nozzle Cascade with Convergent Meridional Shape of Flow Path | |
Diehl | Mitigation of tip leakage induced phenomena in a low Reynolds number centrifugal compressor via blade loading distribution |
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 |