CN104088671B - 基于多参数影响下的气膜冷却效果正交预测方法 - Google Patents
基于多参数影响下的气膜冷却效果正交预测方法 Download PDFInfo
- Publication number
- CN104088671B CN104088671B CN201410326373.8A CN201410326373A CN104088671B CN 104088671 B CN104088671 B CN 104088671B CN 201410326373 A CN201410326373 A CN 201410326373A CN 104088671 B CN104088671 B CN 104088671B
- Authority
- CN
- China
- Prior art keywords
- factor
- test
- mistake
- cause
- sum
- 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.)
- Expired - Fee Related
Links
Landscapes
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明是一种基于多参数影响下的气膜冷却正交预测方法,其特点是,包括的步骤有:选定涡轮叶片冷却效果影响因素;根据正交表安排试验,制定试验计划后进行试验;进行试验结果直观分析,得出较优的生产条件;进行试验结果的方差分析,区分因素水平变化引起的试验结果间的差异与误差波动引起的试验结果间的差异,用最少的试验次数得到最优的方案,简单易行,同时降低了研究时可能出现较大误差的几率,使预测结果更准确。
Description
技术领域
本发明涉及涡轮动力装置中的气膜冷却技术领域,尤其涉及一种基于多参数影响下的气膜冷却效果下的正交预测方法,适用于暴露在高温环境中且需要被冷却的涡轮叶片。
背景技术
燃气轮机的热效率随着燃气初温的提高而增加。现代燃气轮机设计的进口运行温度很高,远超出当前材料的温度极限,受到部件材料耐高温性能的限制,需要对材料表面进行有效冷却。研究和改善燃气轮机冷却技术,以降低高温部件工作温度和延长金属材料的使用寿命,是当前改进燃气轮机安全性的重要手段。
燃气轮机的冷却一般是通过抽取部分压缩空气并旁通燃烧室直接进入透平来实现的。气膜冷却技术作为保护高温部件非常有效的冷却手段,在现代高负荷航空发动机中被广泛采用。气膜冷却是让冷却气体从内部通过射流孔直接喷射到燃气通道表面,以达到保护入射区域以及射流下游区域表面的目的。涡轮叶片气膜冷却主要依赖于冷却工质和高温燃气主流的压比、温比,主流湍流度和气膜冷却叶片上的冷却孔的位置、形状与分布。
本领域针对叶片冷却效果的研究已经做了大量工作,例如:专利公开(公告)号CN103244196A,提出了一种离散气膜冷却孔型,比较了不同孔型、不同吹风比对冷却效果的影响;专利公开(公告)号CN101832154A,提出了将固态工质NH4HCO3作为冷却工质来源,比较NH4HCO3与传统冷却气体的冷却效果。目前大多数的研究都集中在一、两种参数对冷却效果的影响,并未对三种及以上参数共同作用下的冷却效果进行研究,影响因素影响大小的顺序也未曾进行研究。为了设计出更好的叶片气膜冷却模式,本发明将正交法应用于涡轮叶片气膜冷却的研究,用最少的试验得到最准确的结果以及影响因素影响大小的顺序,同时对结果进行进一步的直观分析与统计分析,使得计算结果更具有代表性。
发明内容
本发明的目的在于设计出更好的叶片气膜冷却模式,提出一种通用性强且简单可靠的基于多参数影响下的气膜冷却正交预测方法。
正交设计方法是处理多因素试验的一种科学的试验方法,它利用正交表安排试验,只做较少次数的试验便可判断出较优的条件,若再对结果进行简单的统计分析,还可以更全面、更系统地掌握试验结果,做出正确判断。对涡轮叶片冷却效果进行正交预测可以明确以下三个问题:
1.对指标的影响,哪个因素重要,哪个因素不重要;
2.每个因素以哪个水平为好;
3.各因素以什么样的水平搭配起来,指标最好。
本发明的目的是由以下技术方案来实现的:一种基于多参数影响下的气膜冷却正交预测方法,其特征在于,它包括以下步骤:
(1)涡轮叶片冷却效果受多参数共同作用,首先选定影响因素,给出各因素的水平变化范围,各因素的水平需要根据实际情况选择,影响涡轮叶片冷却效果的因素有:①几何参数,孔的几何参数喷射角度、孔间距、孔径比、孔长和孔出口形状,以及叶片几何参数的叶片前缘形状、曲率和表面粗糙度;②气动参数,吹风比、不稳定尾流、自由流湍流、主流密度比和压力梯度;③其他因素,间隙泄露;
(2)根据正交表安排试验,安排试验时应遵循“均匀分散性,整齐可比”的原则,要求每一列中,不同的数字出现的次数相等;任意两列中数字的排列方式齐全而且均衡,制定试验计划后进行试验;
(3)进行试验结果直观分析,得出较优的生产条件:
Step1:根据结果选则冷却效果最好的工况为最优工况,
Step2:计算每一因素在不同水平条件下的冷却效果平均值,
Step3:计算各因素不同水平下的极差值,因素极差越大,对冷却效果的影响越大,极差越小,对冷却效果影响就越小,根据极差大小判断因素的重要程度,
Step4:综合分析;
(4)进行试验结果的方差分析,区分因素水平变化引起的试验结果间的差异与误差波动引起的试验结果间的差异,
Step1:试验的误差分析
S误=各因素(数据-平均值)2的和(1)
S误称为误差的偏差平方和,为消除数据个数的影响,采用误差的平均偏差平方和V误表示:
f误称为误差的自由度,f误=各条件下(数据个数-1)之和(3)
Step2:试验的因素水平变化分析
S因=因素各水平[重复数×(平均值-总平均值)2]之和(4)
S因称为因素的变动平方和,为消除数据个数的影响,可采用平均变动平方和V因表示:
f因称为因素的自由度,f因=因素的水平数-1(6)
Step3:总的偏差平方和及总自由度
S总=各(数据-总平均)2之和=S误+S因(7)
f总=f误+S因=试验数据总个数-1(8)
为了方便起见,利用下式进行计算:
S总=各(数据)2之和-CT(10)
CT称为修正项:
Step4:显著性检验
根据计算结果判断因素对指标影响的显著性,所谓显著,是指因素水平改变时,确使试验结果的真值有改变,比较V因与V误,计算两者比值F
查找F分布表,临界值记为Fa,当F>Fa,有(1-a)的把握说明该因素对指标有显著性影响,a值叫做显著性水平。
所述的步骤(1)中,燃气轮机叶片中压比为1.02-1.10,相应的喷气比为0.5-2.0,温比为0.5-0.85,相应的密度比为2.0-1.5。
本发明的基于多参数影响下气膜冷却正交预测方法的优点在于:
1.采用正交方法安排试验,对多参数共同作用下的冷却效果进行预测,不仅能够得到单一因素对冷却效果的影响,还能得到多因素共同作用下的最优方案;
2.对涡轮叶片冷却效果进行直观分析与方差分析,可以得到不同参数影响的重要性大小,同时还可以区分试验结果间的差异变化原因;
3.采用正交方法预测冷却效果,可以用最少的试验次数得到最优的方案,简单易行,同时降低了研究时可能出现较大误差的几率,使预测结果更准确。
附图说明
图1为基于多参数影响下的气膜冷却正交预测方法流程图;
图2为应用于平板的开槽姊妹冷却孔三维示意图;
图3为应用于平板的姊妹冷却孔三维示意图;
图4为应用于平板的开槽圆柱冷却孔三维示意图。
具体实施方式
下面结合附图与实施例对本发明进行详细说明。
如图1所示,本发明的一种基于多参数影响下的气膜冷却正交预测方法,包括以下步骤:
(1)涡轮叶片冷却效果受多参数共同作用,首先选定影响因素,给出各因素的水平变化范围,各因素的水平需要根据实际情况选择,影响涡轮叶片冷却效果的因素有:①几何参数,孔的几何参数喷射角度、孔间距、孔径比、孔长和孔出口形状,以及叶片几何参数的叶片前缘形状、曲率和表面粗糙度;②气动参数,吹风比、不稳定尾流、自由流湍流、主流密度比和压力梯度;③其他因素,间隙泄露;
(2)根据正交表安排试验,安排试验时应遵循“均匀分散性,整齐可比”的原则,要求每一列中,不同的数字出现的次数相等;任意两列中数字的排列方式齐全而且均衡,制定试验计划后进行试验;
(3)进行试验结果直观分析,得出较优的生产条件:
Step1:根据结果选则冷却效果最好的工况为最优工况,
Step2:计算每一因素在不同水平条件下的冷却效果平均值,
Step3:计算各因素不同水平下的极差值,因素极差越大,对冷却效果的影响越大,极差越小,对冷却效果影响就越小,根据极差大小判断因素的重要程度,
Step4:综合分析;
(4)进行试验结果的方差分析,区分因素水平变化引起的试验结果间的差异与误差波动引起的试验结果间的差异,
Step1:试验的误差分析
S误=各因素(数据-平均值)2的和(1)
S误称为误差的偏差平方和,为消除数据个数的影响,采用误差的平均偏差平方和V误表示:
f误称为误差的自由度,f误=各条件下(数据个数-1)之和(3)
Step2:试验的因素水平变化分析
S因=因素各水平[重复数×(平均值-总平均值)2]之和(4)
S因称为因素的变动平方和,为消除数据个数的影响,可采用平均变动平方和V因表示:
f因称为因素的自由度,f因=因素的水平数-1(6)
Step3:总的偏差平方和及总自由度
S总=各(数据-总平均)2之和=S误+S因(7)
f总=f误+S因=试验数据总个数-1(8)
为了方便起见,利用下式进行计算:
S总=各(数据)2之和-CT(10)
CT称为修正项:
Step4:显著性检验
根据计算结果判断因素对指标影响的显著性,所谓显著,是指因素水平改变时,确使试验结果的真值有改变,比较V因与V误,计算两者比值F
查找F分布表,临界值记为Fa,当F>Fa,有(1-a)的把握说明该因素对指标有显著性影响,a值叫做显著性水平。
所述的步骤(1)中,燃气轮机叶片中压比为1.02-1.10,相应的喷气比为0.5-2.0,温比为0.5-0.85,相应的密度比为2.0-1.5。
下面是本发明的一个具体实施例:
本发明的基于多参数影响下的气膜冷却正交预测方法以三个因素,每个因素取三个水平为例对涡轮首级静叶冷却效果进行研究。平板气膜冷却是研究复杂流动问题的基础,为了简化计算,以平板代替叶片进行冷却效率的研究。在对三因素三水平的冷却效果进行研究时,传统方法需要做27次试验,采用正交方法只需9次试验,就能得到最优叶片冷却模式,能够大大节省工作时间。
(1)选择对首级静叶冷却效果影响较大的吹风比、孔型与湍流度这三个因素进行研究,吹风比变化范围为1.0、1.5和2.0,湍流度变化范围为10%、20%和30%,孔型选用开槽姊妹孔、姊妹孔与开槽圆柱孔,如图2、3、4所示,分别为应用于平板的开槽姊妹孔、姊妹孔与开槽圆柱孔的三维示意图,三种孔型开孔率相同,孔轴线与平面的夹角都为45°,姊妹孔次孔中心到主孔纵向中心线距离为0.75D,主孔中心到次孔横向中心线的距离在0.75D,横槽深度为D并且与流动方向夹角为45°。
(2)根据试验要求选用L9(34)正交表,表头设计及试验方案如表1所示,第四列没有安排因素直接省略,选用fluent软件进行模拟,并按照正交方案设计的工况进行计算。
表1
(3)进行试验直观分析
Step1:按照9次试验结果可知,9号试验区域平均冷却效率最高,为0.278。试验条件为A3B3C2;
Step2:计算每一因素在不同水平下的平均冷却效率,
Step3:计算各因素不同水平的极差值,
R1=|0.273-0.121|=0.152
R2=|0.222-0.199|=0.023
R3=|0.214-0.201|=0.013;
Step4:综合分析。
根据极差大小,可知孔型对区域平均冷却效率影响最大,取三水平最好。吹风比影响较小,应取一水平最好。而湍流度影响最小,说明改变湍流度对结果影响不大,从试验结果看取二水平最好。经综合分析最好的叶片气膜冷却模式为A3B1C2,这是9次试验中未做过的试验,因此需要对A3B1C2条件进行试验。试验结果表明A3B1C2条件下的区域平均冷却效率为0.267,小于A3B3C2的冷却效率0.278,因此最终确定A3B3C2为最优气膜冷却模式。
(4)进行试验结果的方差分析
Step1:试验的误差分析
计算误差的偏差平方和可用正交表中未安排因素的空白列计算。本例中用第四列的Ⅰ4,Ⅱ4,Ⅲ4的平均值分别与之差平方和的三倍计算S空,因为没有安排因素,所以计算S空的偏差平方和中,也就没有因素水平间的差异造成的偏差,该列仅仅反应了试验误差的大小,因此
f误=空列水平数-1=3-1=2
Step2:试验的因素水平变化分析
因素A的偏差平方和
同样可求出SB=4.6×10-4SC=0.1×10-4
自由度fA=fB=fC=3-1=2
Step3:总的偏差平方和及总自由度
f总=9-1=8;
Step4:显著性检验
采用单因素试验的方差分析方法进行显著性检验,计算V因与V误的比值,而后与F分布表的临界值进行比较,根据以上计算结果可分别计算VA,VB,VC以及V误,
由计算结果可知,VB,VC与V误相差不大,说明他们的偏差中,由因素水平变化的影响部分很小,他们的偏差实际上主要是由误差干扰造成的,将SB,SC与S误合并在一起,用以估计误差影响的大小,而且误差的自由度越大,进行显著性检验时越灵敏,
S误′=SB+SC+S误=2.67×10-4
f误′=fB+fC+f误=6
检验余下因素A的显著性:
因素A的自由度fA=2,S误′的自由度f误′=6,查a=0.01显著性水平的F分布表,F0.01(2,6)=10.9:
FA=422.2>F0.01(2,6)=10.9
所以,因素A是高度显著的。
a)通过本例的方差分析可知:孔型对叶片冷却效率有高度显著的影响,取三水平最好;
b)吹风比与湍流度的改变对冷却效率影响不大,综合考虑孔型、吹风比与湍流度三因素的影响,最优气膜冷却模式为A3B3C2,即选用开槽圆柱孔,吹风比为2.0,湍流度为15%是最优冷却方案。
Claims (2)
1.一种基于多参数影响下的气膜冷却正交预测方法,其特征在于,它包括以下步骤:
(1)涡轮叶片冷却效果受多参数共同作用,首先选定影响因素,给出各因素的水平变化范围,各因素的水平需要根据实际情况选择,影响涡轮叶片冷却效果的因素有:①几何参数,孔的几何参数喷射角度、孔间距、孔径比、孔长和孔出口形状,以及叶片几何参数的叶片前缘形状、曲率和表面粗糙度;②气动参数,吹风比、不稳定尾流、自由流湍流、主流密度比和压力梯度;③其他因素,间隙泄露;
(2)根据正交表安排试验,安排试验时应遵循“均匀分散性,整齐可比”的原则,要求每一列中,不同的数字出现的次数相等;任意两列中数字的排列方式齐全而且均衡,制定试验计划后进行试验;
(3)进行试验结果直观分析,得出较优的生产条件:
Step1:根据结果选则冷却效果最好的工况为最优工况,
Step2:计算每一因素在不同水平条件下的冷却效果平均值,
Step3:计算各因素不同水平下的极差值,因素极差越大,对冷却效果的影响越大,极差越小,对冷却效果影响就越小,根据极差大小判断因素的重要程度,
Step4:综合分析;
(4)进行试验结果的方差分析,区分因素水平变化引起的试验结果间的差异与误差波动引起的试验结果间的差异,
Step1:试验的误差分析
S误=各因素(数据-平均值)2的和(1)
S误称为误差的偏差平方和,为消除数据个数的影响,采用误差的平均偏差平方和V误表示:
f误称为误差的自由度,f误=各条件下(数据个数-1)之和(3)
Step2:试验的因素水平变化分析
S因=因素各水平[重复数×(平均值-总平均值)2]之和(4)
S因称为因素的变动平方和,为消除数据个数的影响,可采用平均变动平方和V因表示:
f因称为因素的自由度,f因=因素的水平数-1(6)
Step3:总的偏差平方和及总自由度
S总=各(数据-总平均)2之和=S误+S因(7)
f总=f误+S因=试验数据总个数-1(8)
为了方便起见,利用下式进行计算:
S总=各(数据)2之和-CT(10)
CT称为修正项:
Step4:显著性检验
根据计算结果判断因素对指标影响的显著性,所谓显著,是指因素水平改变时,确使试验结果的真值有改变,比较V因与V误,计算两者比值F
查找F分布表,临界值记为Fa,当F>Fa,有(1-a)的把握说明该因素对指标有显著性影响,a值叫做显著性水平。
2.根据权利要求1所述的一种基于多参数影响下的气膜冷却正交预测方法,其特征在于,所述的步骤(1)中,燃气轮机叶片中压比为1.02-1.10,相应的喷气比为0.5-2.0,温比为0.5-0.85,相应的密度比为2.0-1.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410326373.8A CN104088671B (zh) | 2014-07-10 | 2014-07-10 | 基于多参数影响下的气膜冷却效果正交预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410326373.8A CN104088671B (zh) | 2014-07-10 | 2014-07-10 | 基于多参数影响下的气膜冷却效果正交预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104088671A CN104088671A (zh) | 2014-10-08 |
CN104088671B true CN104088671B (zh) | 2016-01-20 |
Family
ID=51636416
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410326373.8A Expired - Fee Related CN104088671B (zh) | 2014-07-10 | 2014-07-10 | 基于多参数影响下的气膜冷却效果正交预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104088671B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109614728A (zh) * | 2018-12-18 | 2019-04-12 | 济南大学 | 一种齿轮摩擦因数计算公式的分析方法 |
CN112084597B (zh) * | 2020-09-08 | 2021-06-15 | 北京航空航天大学 | 基于钟形曲线的单排气膜冷却效率二维分布ai预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2025870A2 (en) * | 2007-07-31 | 2009-02-18 | Ansaldo Energia S.P.A. | Method for modifying a gas turbine blade |
CN101655416A (zh) * | 2009-09-11 | 2010-02-24 | 北京航空航天大学 | 超音速气膜冷却实验系统 |
CN103076362A (zh) * | 2013-01-08 | 2013-05-01 | 中国科学院工程热物理研究所 | 一种测量含化学热沉气膜冷却效率的方法和装置 |
CN103398835A (zh) * | 2013-08-21 | 2013-11-20 | 中国人民解放军国防科学技术大学 | 基于高超声速炮风洞的气膜冷却瞬态热流测试系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9188504B2 (en) * | 2012-10-31 | 2015-11-17 | General Electric Company | Methods for testing turbine blades |
-
2014
- 2014-07-10 CN CN201410326373.8A patent/CN104088671B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2025870A2 (en) * | 2007-07-31 | 2009-02-18 | Ansaldo Energia S.P.A. | Method for modifying a gas turbine blade |
CN101655416A (zh) * | 2009-09-11 | 2010-02-24 | 北京航空航天大学 | 超音速气膜冷却实验系统 |
CN103076362A (zh) * | 2013-01-08 | 2013-05-01 | 中国科学院工程热物理研究所 | 一种测量含化学热沉气膜冷却效率的方法和装置 |
CN103398835A (zh) * | 2013-08-21 | 2013-11-20 | 中国人民解放军国防科学技术大学 | 基于高超声速炮风洞的气膜冷却瞬态热流测试系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104088671A (zh) | 2014-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Enhanced film cooling effectiveness with new shaped holes | |
US9483605B2 (en) | Probabilistic high cycle fatigue (HCF) design optimization process | |
Bacci et al. | Flow field and hot streak migration through a high pressure cooled vanes with representative lean burn combustor outflow | |
US9982540B2 (en) | Interwoven channels for internal cooling of airfoil | |
Lee et al. | Shape optimization of a laidback fan-shaped film-cooling hole to enhance cooling performance | |
CN104088671B (zh) | 基于多参数影响下的气膜冷却效果正交预测方法 | |
Li et al. | Experimental investigation of endwall heat transfer with film and impingement cooling | |
Lee et al. | Numerical analysis of film-cooling performance and optimization for a novel shaped film-cooling hole | |
Khalil et al. | Numerical Predictions of Three-Dimensional Unsteady Turbulent Film-Cooling for Trailing Edge of Gas-Turbine Blade Using Large Eddy Simulation | |
Pichler et al. | LES and RANS Analysis of the End-Wall Flow in a Linear LPT Cascade: Part I—Flow and Secondary Vorticity Fields Under Varying Inlet Condition | |
Sperling et al. | Film Cooling Hole Shape Effects on Turbine Blade Heat Transfer–Part I: Computational Comparison to Experiment | |
Leo´ n De Paz et al. | A numerical study of an impingement array inside a three dimensional turbine vane | |
Zhang et al. | Discharge coefficients and aerodynamic losses for cylindrical and cratered film-cooling holes with various coolant crossflow orientations | |
Verma et al. | Modeling of combustor and turbine vane interaction | |
Wang et al. | Experimental and numerical study of heat transfer performance for an engine representative two-pass rotating internal cooling channel | |
Saha et al. | Suction and Pressure Side Film Cooling Influence on Vane Aero Performance in a Transonic Annular Cascade | |
Yang et al. | A new multi-objective reliability-based robust design optimization method | |
Song et al. | Knowledge-based aero-thermal multi-disciplinary design optimization of a high temperature blade | |
KR20120013501A (ko) | 축류압축기의 케이싱 그루브 최적설계방법 | |
Andreini et al. | Investigation of Circular and Shaped Effusion Cooling Arrays for Combustor Liner Application—Part 2: Numerical Analysis | |
Ghodake et al. | Modification of noise sources by virtue of blade sweep in low speed fan | |
Pawar et al. | Evaluating the effects of multiple hole profiles on gas turbine blade cooling rate: a computational study | |
Zhang et al. | Co-Optimization of Turbine Blade Aero and Thermal Designs Based on Computational Fluid Dynamics (CFD) Models | |
Yang | Numerical study of effusion cooling flow and heat transfer | |
김윤기 | Variable Fidelity Optimization of Film Cooling Hole Arrangements on High Pressure Turbine Considering Conjugate Heat Transfer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 Termination date: 20200710 |
|
CF01 | Termination of patent right due to non-payment of annual fee |