CN117386664A - 分段积叠的轴流压气机的转子叶片及其设计方法 - Google Patents
分段积叠的轴流压气机的转子叶片及其设计方法 Download PDFInfo
- Publication number
- CN117386664A CN117386664A CN202210787403.XA CN202210787403A CN117386664A CN 117386664 A CN117386664 A CN 117386664A CN 202210787403 A CN202210787403 A CN 202210787403A CN 117386664 A CN117386664 A CN 117386664A
- Authority
- CN
- China
- Prior art keywords
- blade
- point
- rotor blade
- section
- stacking line
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000005484 gravity Effects 0.000 claims abstract description 74
- 238000003475 lamination Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 abstract description 17
- 230000000694 effects Effects 0.000 abstract description 11
- 238000005457 optimization Methods 0.000 description 14
- 230000007704 transition Effects 0.000 description 9
- 239000012530 fluid Substances 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 4
- 230000002349 favourable effect Effects 0.000 description 4
- 238000012886 linear function Methods 0.000 description 4
- 230000035939 shock Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000306 component Substances 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/26—Rotors specially for elastic fluids
- F04D29/32—Rotors specially for elastic fluids for axial flow pumps
- F04D29/38—Blades
- F04D29/388—Blades characterised by construction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Mechanical Engineering (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明公开了一种分段积叠的轴流压气机的转子叶片及其设计方法,该设计方法包括:确定积叠线的轴向位置;确定积叠线的周向位置:积叠线在周向上为S型,积叠线沿高度方向由叶根至叶尖依次包括A、B、C和D四个控制点;B点至C点的第一区间中的各基元叶型截面的重心位置点的坐标,符合特定公式;根据确定后的积叠线的轴向位置和周向位置,确定积叠线的线型。本发明通过设置转子叶片的积叠线在周向上为S型,以在径向方向产生逆压力梯度,使得两端的气体向叶中迁移,起到改善端区流动的目的。从而减弱了叶尖区域的二次流影响;并通过设置特定的函数,以通过构造简便的方程来控制积叠线的周向偏移量,提高优化叶型的计算效率。
Description
技术领域
本发明涉及一种分段积叠的轴流压气机的转子叶片及其设计方法。
背景技术
随着科技的发展,现代航空发动机对经济性和可靠性的要求越来越高,航空发动机各部件的性能遇到了前所未有的挑战。轴流压气机是航空发动机的核心部件之一,其性能对整个发动机有着举足轻重的影响。然而,随着航空发动机性能指标的不断提高,发动机内部的温度和压比逐渐升高,压气机内流系统的损失问题日趋严重。多级高压比轴流压气机内流损失严重影响其气动性能,进而影响整个发动机性能。
压气机损失主要来源为叶型损失、端壁损失和泄漏损失。高负荷压气机端壁区域流动十分复杂,端壁附面层损失和叶尖泄漏损失是压气机转子叶片损失的主要组成部分。
压气机是航空发动机内部承受逆压梯度最大的旋转部件,实现高效的气流增压是压气机设计的一大难题。在转子的轮毂和叶片叶根吸力面组成的角区区域,由于轮毂附面层和气体潜流的影响,易发生分离,产生气动损失,造成压气机的效率下降。在叶尖区域,转子和外机匣之间会存在间隙,转子叶片压力面的气流会通过该间隙泄漏到叶片的吸力面,产生泄漏涡等复杂的二次流特征,加之转子叶片的切线速度较高,在压气机的前面级转子叶尖会产生激波,激波和泄漏涡及外机匣附面层相互干扰,使得叶尖间隙区域的流动及其复杂,存在着多种复杂特征的涡系结构。大量研究表明,压气机的失速通常发生在叶尖和轮毂区域。
因此,如何有效对叶尖区域激波和流场进行控制,从而减弱叶尖区域二次流,是本领域有待解决的一个难题。
另一方面,压气机的设计通常需要多轮优化迭代才能获得理想的压气机气动叶型。为了节省优化时间和避免过渡依赖设计人员的设计经验,压气机气动优化方法近年来获得了设计者的青睐。压气机设计优化方法是将计算流体力学与数值优化方法相结合应用的一种跨学科设计方法,其基本思想是将气动性能作为目标函数值,借助CFD三维数值仿真工具,利用优化算法不断迭代寻找最优的压气机气动叶型。叶型设计的参数化选取方法是压气机气动优化的最为关键技术之一。如何提取压气机周向和径向积叠形式的关键参数并用合适的数学方程式对其进行描述,使其兼顾压气机叶型设计特征和计算方法可行性是压气机叶型积叠形式优化设计的一个难题。
发明内容
本发明要解决的技术问题是为了克服现有技术中分段积叠的轴流压气机的转子叶片存在气流损失的缺陷,同时提高叶型优化设计效率,提供一种分段积叠的轴流压气机的转子叶片及其设计方法。
本发明是通过下述技术方案来解决上述技术问题:
一种分段积叠的轴流压气机的转子叶片的设计方法,其特点在于,所述转子叶片包括叶根和叶尖,所述转子叶片沿高度方向从所述叶根到所述叶尖具有若干个基元叶型截面,将各所述基元叶型截面的重心作为积叠点沿所述转子叶片的高度方向进行叠加得到积叠线,所述设计方法包括以下步骤:
S1:确定所述积叠线的轴向位置;
S2:确定所述积叠线的周向位置:所述积叠线在周向上为“S”型,所述积叠线沿高度方向由所述叶根至所述叶尖依次包括A、B、C和D四个控制点,A点为所述叶根的基元叶型截面的重心位置,且A点为原点;B点位于所述积叠线在X轴正向上的最大周向偏移量位置点,且B点为预设点;C点位于所述积叠线在X轴负向上的最大周向偏移量位置点,且C点为预设点;D点为所述叶尖的基元叶型截面的重心位置,且D点在X轴上位于B点和C点之间;
B点至C点的第一区间中的各所述基元叶型截面的重心位置点的坐标(x1,y1),符合以下公式:
其中,N=(xb–xc)/2;α=Pi/(yc–yb),Pi=3.1415926;参数和M根据B点的坐标(xb,yb)和C点的坐标(xc,yc)代入公式得到;
S3:根据确定后的所述积叠线上的每一所述基元叶型截面的重心的轴向位置和周向位置,确定所述积叠线的线型;将基元叶型沿所述积叠线的线型进行延伸构造曲线组,得到三维的转子叶片。
在本技术方案中,通过设置转子叶片的积叠线在周向上为“S”型,可以在径向方向会产生逆压力梯度,有利于叶尖区域的低能流体往叶中迁移,进而减弱了叶尖区域二次流影响,改善了叶尖流动状况,提升了压气机转子叶片的稳定工作范围。进一步,依据损失最小或效率最佳原则确定最优的最大周向偏移量xb和xc及相应高度和延伸长度(定义为周向往叶盆方向偏移为正,往叶背方向为负),并通过设置第一区间中的各基元叶型截面的重心位置点的坐标符合特定的函数,以通过构造简便的方程来控制位于的第一区间中的各基元叶型截面的重心位置的周向偏移量,提高优化叶型的计算效率。也就是说,这相对于现有技术中,通过试错的方法来构造积叠线的方式,大大地提高了计算效率。需要说明的是,转子叶片的积叠线在周向上为“S”型,也包括转子叶片的积叠线在周向上为近似“S”型的情况。
较佳地,在步骤S2中,包括:
将所述第一区间在B点沿高度方向向下延伸第一延伸高度Δyb,得到第一中间点E;
A点至E点的第二区间中的各所述基元叶型截面重心的周向位置点的坐标(x2,y2),符合以下公式:
其中,参数h、i、j通过E点的坐标(xe,ye)和E点的斜率ke,以及A点的坐标(xa,ya)代入公式得到。
在本技术方案中,通过将所述第一区间在B点沿高度方向向下延伸第一延伸高度Δyb,能够使积叠线在转向位置保证曲率连续,从而达到平滑地过渡,有利于保持较高的压气机效率情况下,转子叶片保持较佳强度水平;并且,通过设置第二区间中的各基元叶型截面的重心位置点的坐标符合特定的函数,以通过构造简便的方程来控制位于的第二区间中的各基元叶型截面的重心位置的周向偏移量,提高优化叶型的计算效率。
较佳地,所述第一延伸高度Δyb为所述第一区间的高度的5%~20%。
在本技术方案中,通过设置第一延伸高度Δyb范围,以有效保证曲率的延续以及产生有效的逆压力梯度,即避免该高度过小不能达到保证曲率连续,而过大影响在径向方向产生有效的逆压力梯度。
较佳地,在步骤S2中,包括:
将所述第一区间在C点沿高度方向向上延伸第二延伸高度Δyc,得到中间点F;
F点至D点的第三区间中的各所述基元叶型截面的重心位置点的坐标(x3,y3),符合以下公式:
x3=kfy3+Q
其中,kf为F点的斜率,参数Q通过F点的坐标(xf,yf)和F点的斜率kf代入公式得到。
在本技术方案中,通过将所述第一区间在C点沿高度方向向上延伸第一延伸高度Δyc,能够使积叠线在转向位置保证曲率连续,从而达到平滑地过渡,有利于保持较高的压气机效率情况下,转子叶片保持较佳强度水平;并且,通过设置第三区间中的各基元叶型截面的重心位置点的坐标符合特定的函数,以通过构造简便的方程来控制位于的第二区间中的各基元叶型截面的重心位置的周向偏移量,提高优化叶型的计算效率。此外,通过C点向上延伸高度搭配线性函数的方法,可以使得叶型在叶尖部位形成类似“兜风”墙的效果,避免叶尖叶盆处的高压气体潜移至叶背,减弱叶尖二次流效果。
较佳地,所述第二延伸高度Δyc为所述第一区间的高度的5%~20%。
在本技术方案中,通过设置第二延伸高度Δyc范围,以有效保证曲率的延续以及产生有效的逆压力梯度,即避免该高度过小不能达到保证曲率连续,而过大影响在径向方向产生有效的逆压力梯度。
较佳地,在步骤S2中,在B点的坐标(xb,yb)中,yb为叶高的25%~35%;和/或,
在步骤S2中,在C点的坐标(xc,yc)中,yc为叶高的70%~85%。
在本技术方案中,通过设置B点的坐标(xb,yb)中yb的取值范围,有利于保持较高的压气机效率情况下,兼顾转子叶片保持较佳强度水平。
通过设置C点的坐标(xc,yc)中yc的取值范围,有利于保持较高的压气机效率情况下,兼顾转子叶片保持较佳强度水平。
较佳地,在步骤S1中,所述确定所述积叠线的轴向位置,包括:
在轴向上,将各所述基元叶型截面的重心位置以所述叶根的所述基元叶型截面的重心位置为基准点沿高度方向移至一条直线上;设定叶尖的所述基元叶型截面的重心位置的轴向偏移量为最大轴向偏移量,以确定叶尖的所述基元叶型截面的重心位置;其余所述基元叶型截面的重心位置的轴向偏移量采用线性插值的方式得到。
在本技术方案中,通过上述方式得到各基元叶型截面的重心位置的轴向偏移量,提高优化叶型的计算效率。
较佳地,在步骤S1中,在轴向上所述叶尖的所述基元叶型截面的重心位置比叶根的所述基元叶型截面的重心位置靠前。
在本技术方案中,叶尖的基元叶型截面的重心位置在轴向上比叶根的基元叶型截面的重心位置靠前,这样会形成前掠的转子叶型,有助于提高转子叶片的裕度。
较佳地,在步骤S3之后,还包括:
S4:重复步骤S1至S3多次,获得多个不同的所述转子叶片,获得多个不同的所述转子叶片所对应的压气机效率,选择最高的所述压气机效率所对应的所述转子叶片作为最终的转子叶片。
在本技术方案中,在步骤S4中,将多个不同的最大轴向积叠偏移量代入步骤S1,以及将多个不同的最大周向积叠偏移量代入步骤S2,重复步骤S1至S3多次,得到多个不同的所述转子叶片,即通过迭代的方法,得到多个不同的所述转子叶片,该些转子叶片对应不同的最大轴向积叠偏移量和不同的最大周向积叠偏移量。
对多个不同的所述转子叶片进行压气机效率计算,选择最高的所述压气机效率所对应的所述转子叶片作为最终的转子叶片,以得到压气机效率最优的转子叶片,即得到最优的周向积叠和轴向积叠的转子叶片。
一种分段积叠的轴流压气机的转子叶片,其特点在于,所述转子叶片通过如上所述的分段积叠的轴流压气机的转子叶片的设计方法得到。
本发明的积极进步效果在于:
本发明通过设置转子叶片的积叠线在周向上为S型,可以在径向方向会产生逆压力梯度,有利于叶尖区域的低能流体往叶中迁移,进而减弱了叶尖区域二次流影响,改善了叶尖流动状况,提升了压气机转子叶片的稳定工作范围。进一步,依据损失最小或效率最佳原则确定最优的最大周向偏移量xb和xc及相应高度和延伸长度,并通过设置第一区间中的各基元叶型截面的重心位置点的坐标符合特定的函数,以通过构造简便的方程来控制位于的第一区间中的各基元叶型截面的重心位置的周向偏移量,提高优化叶型的计算效率。也就是说,这相对于现有技术中,通过试错的方法来构造积叠线的方式,大大地提高了计算效率。
另一方面,现有技术中用于描述近似“S”型曲线的方程多为单方程形式;而在压气机的积叠叶型优化中,通常涉及多个变量的优化,如最大周向偏移量及其位置,和根尖位置的倾斜角度等,此时单一方程形式无法对压气机周向积叠的关键特征参数进行优化。本发明通过提取压气机周向积叠最为关键的设计参数,最大正、负向周向偏移量及其高度位置,和根尖位置的倾斜角度,并将其作为输入量即自变量,采用分段函数形式对其余高度位置的偏移量进行封闭求解,并保证了分段函数连接点处的光滑过渡,以及至少一阶导数连续。在保证压气机叶型周向积叠关键参数能够被完全描述的基础上,通过构造形式简单的分段函数,使得方程易于求解,提升了叶型优化效率。也就是说,本发明给出了一种将压气机周向积叠最大正、负向周向偏移量及其高度位置,和根尖位置的倾斜角度作为控制参数的参数化方法,采用分段函数形式对其余高度位置的偏移量进行封闭求解,保证了分段函数连接点处的光滑过渡,以及至少一阶导数连续。
附图说明
图1为本发明一较佳实施例分段积叠的轴流压气机的转子叶片的设计方法的流程图。
图2为本发明一较佳实施例分段积叠的轴流压气机的转子叶片的叶根、叶尖、以及中间位置的基元叶型截面及其上的重心位置的结构示意图。
图3为本发明一较佳实施例分段积叠的轴流压气机的转子叶片的积叠线轴向偏移量示意图。
图4为本发明一较佳实施例分段积叠的轴流压气机的转子叶片的积叠线周向偏移量示意图。
图5为本发明一较佳实施例分段积叠的轴流压气机的转子叶片的设计方法的流程图。
图6为本发明一较佳实施例分段积叠的轴流压气机的转子叶片的最终叶型的结构示意图。
图7为本发明一较佳实施例分段积叠的轴流压气机的转子叶片的原始叶型的结构示意图。
具体实施方式
下面通过实施例的方式进一步说明本发明,但并不因此将本发明限制在所述的实施例范围之中。
本实施例提供一种分段积叠的轴流压气机的转子叶片的设计方法。
所述转子叶片包括叶根和叶尖,所述转子叶片沿高度方向从所述叶根到所述叶尖具有若干个基元叶型截面,将各所述基元叶型截面的重心作为积叠点沿所述转子叶片的高度方向进行叠加得到积叠线。
如图1所示,本实施例的分段积叠的轴流压气机的转子叶片的设计方法包括以下步骤:
S1:确定所述积叠线的轴向位置;
S2:确定所述积叠线的周向位置:所述积叠线在周向上为“S”型,所述积叠线沿高度方向由所述叶根至所述叶尖依次包括A、B、C和D四个控制点,A点为所述叶根的基元叶型截面的重心位置,且A点为原点;B点位于所述积叠线在X轴正向上的最大周向偏移量位置点,且B点为预设点;C点位于所述积叠线在X轴负向上的最大周向偏移量位置点,且C点为预设点;D点为所述叶尖的基元叶型截面的重心位置,且D点在X轴上位于B点和C点之间;
B点至C点的第一区间中的各所述基元叶型截面的重心位置点的坐标(x1,y1),符合以下公式:
其中,N=(xb–xc)/2;α=Pi/(yc–yb),Pi=3.1415926;参数和M根据B点的坐标(xb,yb)和C点的坐标(xc,yc)代入公式得到;
S3:根据确定后的所述积叠线上的每一所述基元叶型截面的重心的轴向位置和周向位置,确定所述积叠线的线型;将基元叶型沿所述积叠线的线型进行延伸构造曲线组,得到三维的转子叶片。
这样,通过设置转子叶片的积叠线在周向上为S型,可以在径向方向会产生逆压力梯度,有利于叶尖区域的低能流体往叶中迁移,进而减弱了叶尖区域二次流影响,改善了叶尖流动状况,提升了压气机转子叶片的稳定工作范围。进一步,依据损失最小或效率最佳原则确定最优的最大周向偏移量xb和xc及相应高度和延伸长度(定义为周向往叶盆方向偏移为正,往叶背方向为负),并通过设置第一区间中的各基元叶型截面的重心位置点的坐标符合特定的函数,以通过构造简便的方程来控制位于的第一区间中的各基元叶型截面的重心位置的周向偏移量,提高优化叶型的计算效率。也就是说,这相对于现有技术中,通过试错的方法来构造积叠线的方式,大大地提高了计算效率。需要说明的是,转子叶片的积叠线在周向上为“S”型,也包括转子叶片的积叠线在周向上为近似“S”型的情况。
具体地,在步骤S1中,所述确定所述积叠线的轴向位置,包括:
将各所述基元叶型截面的重心位置以所述叶根的所述基元叶型截面的重心位置为基准点沿高度方向移至一条直线上;设定叶尖的所述基元叶型截面的重心位置的轴向偏移量为最大轴向偏移量,以确定叶尖的所述基元叶型截面的重心位置;其余所述基元叶型截面的重心位置的轴向偏移量采用线性插值的方式得到。
具体地,如图2所示,图中展示了叶根的基元叶型截面111及其上的重心位置1111,叶尖的基元叶型截面121及其上的重心位置1211,以及位于叶根和叶尖之间的中间位置上的叶中的基元叶型截面131及其上的重心位置1311。
如图3所示,在步骤S1中,确定所述积叠线的轴向位置包括:在轴向O上,将各所述基元叶型截面的重心位置以所述叶根的所述基元叶型截面的重心位置1111为基准点沿高度方向移至一条直线上,即将各所述基元叶型截面的重心位置均放置在初始位置上(图3中用虚线表示);设定叶尖的所述基元叶型截面的重心位置1211的轴向偏移量为最大轴向偏移量Ax1,以确定叶尖的所述基元叶型截面的重心位置1311,即,得到叶尖的所述基元叶型截面的重心位置1311的最终位置(图3中用实线表示)。其中,叶尖的所述基元叶型截面的重心位置1211的轴向偏移量Ax1为预设的值,即,根据需要预先设定叶尖的所述基元叶型截面的重心位置1211的轴向偏移量Ax1。
其余所述基元叶型截面的重心位置的轴向偏移量采用线性插值的方式(线性函数)得到,即通过线性插值的方式得到其余所述基元叶型截面的重心位置的最终位置,也就是说通过求解周向的线性函数,确定整个叶高的轴向积叠偏移量。
这样,通过上述方式得到各基元叶型截面的重心位置的轴向偏移量,提高优化叶型的计算效率。
其中,请复参阅图3,在步骤S1中,在轴向上所述叶尖的所述基元叶型截面的重心位置1211比叶根的所述基元叶型截面的重心位置1111靠前,即,所述叶尖的所述基元叶型截面的重心位置1211向转子叶片的吸力面偏移,以形成前掠的转子叶型,有助于提高转子叶片的裕度;并且,能够起到类似“兜风墙”的效果,进而有效抑制叶尖泄漏。
具体地,如图4所示,在步骤S2中,所述积叠线20在周向上为S型,所述积叠线20沿高度方向H由所述叶根至所述叶尖依次包括A、B、C和D四个控制点,A点位于所述叶根的基元叶型截面的重心位置1111,且A点为原点,即,在二维坐标系中,A点为原点,A点的坐标(xa,ya)为(0,0)。
B点位于所述积叠线在X轴正向上的最大周向偏移量位置点,且B点为预设点,C点位于所述积叠线在X轴负向上的最大周向偏移量位置点,且C点为预设点。也就是说,该参数化方法配合CFD计算,依据损失最小或效率最佳原则确定最优的最大周向偏移量及该位置的相应高度,即,B点的坐标(xb,yb)和C点的坐标(xc,yc)依据损失最小或效率最佳原则进行确定。积叠线在B点位置和C点位置均转向。
D点为位于所述叶尖的基元叶型截面的重心位置1211,且D点在X轴上位于B点和C点之间。
需要说明的是,所述积叠线在周向上为S型,是指积叠线在周向上为近似S型,其具体特征为,沿Y轴正向(即叶片的高度方向)依次设置在积叠线上的A点、B点、C点和D点;同时,B点位于所述积叠线在X轴正向上的最大周向偏移量位置点,C点位于所述积叠线在X轴负向上的最大周向偏移量位置点,A点和D点在X轴上的位置均位于B点和C点在X轴上的位置之间。
B点至C点的第一区间21中的各所述基元叶型截面的重心位置点的坐标(x1,y1),符合以下公式:
其中,N=(xb–xc)/2;α=Pi/(yc–yb),Pi=3.1415926;参数和M根据B点的坐标(xb,yb)和C点的坐标(xc,yc)代入公式得到。
较佳地,在步骤S2中,包括:
将所述第一区间在B点沿高度方向向下延伸第一延伸高度Δyb,得到第一中间点E;
A点至E点的第二区间22中的各所述基元叶型截面重心的周向位置点的坐标(x2,y2),符合以下公式:
其中,参数h、i、j通过E点的坐标(xe,ye)和E点的斜率ke,以及A点的坐标(xa,ya)代入公式得到。E点的斜率ke通过以下公式求得:
其中,参数N、α、从上文的公式中得到。
这样,通过将所述第一区间在B点沿高度方向向下延伸第一延伸高度Δyb,也就是说B点到E点之间的各所述基元叶型截面的重心的周向位置点的坐标能够使积叠线在转向位置保证曲率连续,从而达到平滑地过渡,有利于保持较高的压气机效率情况下,转子叶片保持较佳强度水平;并且,通过设置第二区间中的各基元叶型截面的重心位置点的坐标符合特定的函数,以通过构造简便的方程来控制位于的第二区间中的各基元叶型截面的重心位置的周向偏移量,提高优化叶型的计算效率。其中,将所述第一区间在B点沿高度方向向下延伸第一延伸高度Δyb,得到第一中间点E,也就是说B点到E点之间的各所述基元叶型截面的重心的周向位置点的坐标亦符合B点至C点的第一区间21中的各所述基元叶型截面的重心位置点的坐标所符合的公式。
其中,所述第一延伸高度Δyb为所述第一区间的高度的5%~20%。这样,通过设置第一延伸高度Δyb范围,以有效保证曲率的延续以及产生有效的逆压力梯度,即避免该高度过小不能达到保证曲率连续,而过大影响在径向方向产生有效的逆压力梯度。
较佳地,在步骤S2中,包括:
将所述第一区间在C点沿高度方向向上延伸第二延伸高度Δyc,得到中间点F;
F点至D点的第三区间23中的各所述基元叶型截面的重心位置点的坐标(x3,y3),符合以下公式:
x3=kfy3+Q
其中,kf为F点的斜率,参数Q通过F点的坐标(xf,yf)和F点的斜率kf代入公式得到。F点的斜率kf通过以下公式求得:
其中,参数N、α、从上文的公式中得到。
这样,通过将所述第一区间在C点沿高度方向向上延伸第一延伸高度Δyc,能够使积叠线在转向位置保证曲率连续,从而达到平滑地过渡,有利于保持较高的压气机效率情况下,转子叶片保持较佳强度水平;并且,通过设置第三区间中的各基元叶型截面的重心位置点的坐标符合特定的函数,以通过构造简便的方程来控制位于的第二区间中的各基元叶型截面的重心位置的周向偏移量,提高优化叶型的计算效率。此外,通过C点向上延伸高度搭配线性函数的方法,可以使得叶型在叶尖部位形成类似“兜风”墙的效果,避免叶尖叶盆处的高压气体潜移至叶背,减弱叶尖二次流效果。
其中,将所述第一区间在C点沿高度方向向上延伸第二延伸高度Δyc,得到中间点F,也就是说C点到F点之间的各所述基元叶型截面的重心的周向位置点的坐标亦符合B点至C点的第一区间21中的各所述基元叶型截面的重心位置点的坐标所符合的公式。
其中,所述第二延伸高度Δyc为所述第一区间的高度的5%~20%。这样,通过设置第二延伸高度Δyc范围,以有效保证曲率的延续以及产生有效的逆压力梯度,即避免该高度过小不能达到保证曲率连续,而过大影响在径向方向产生有效的逆压力梯度。
较佳地,在步骤S2中,在B点的坐标(xb,yb)中,yb为叶高的25%~35%。
这样,通过设置B点的坐标(xb,yb)中yb的取值范围,有利于保持较高的压气机效率情况下,兼顾转子叶片保持较佳强度水平。
较佳地,在步骤S2中,在C点的坐标(xc,yc)中,yc为叶高的70%~85%。
这样,通过设置C点的坐标(xc,yc)中yc的取值范围,有利于保持较高的压气机效率情况下,兼顾转子叶片保持较佳强度水平。
也就是说,A点到E点、E点到F点、以及F点到D点的各所述基元叶型截面的重心的周向位置点的坐标所符合的公式各不相同,需要分段求解对应的公式,以确定整个叶高的周向积叠偏移量。
接续,在步骤S3之后,还包括:
S4:重复步骤S1至S3多次,获得多个不同的所述转子叶片,获得多个不同的所述转子叶片所对应的压气机效率,选择最高的所述压气机效率所对应的所述转子叶片作为最终的转子叶片。
在步骤S4中,将多个不同的最大轴向积叠偏移量代入步骤S1,以及将多个不同的最大周向积叠偏移量代入步骤S2,重复步骤S1至S3多次,得到多个不同的所述转子叶片,即通过迭代的方法,得到多个不同的所述转子叶片,该些转子叶片对应不同的最大轴向积叠偏移量和不同的最大周向积叠偏移量。
对多个不同的所述转子叶片进行压气机效率计算,选择最高的所述压气机效率所对应的所述转子叶片作为最终的转子叶片,以得到压气机效率最优的转子叶片,即得到最优的周向积叠和轴向积叠的转子叶片。
即,将气动性能作为目标函数值,借助网格划分、CFD三维数值仿真工具计算气压机效率,利用优化算法不断迭代寻找最优的压气机气动叶型。
需要说明的是,在其他实施方式中,也可以是先完成步骤S2的确定所述积叠线的周向位置,再完成步骤S1的确定所述积叠线的轴向位置,或者步骤S1和步骤S2同时完成。
具体如图5所示,图5为包括上述所有步骤的方法的流程图。
如图6所示,本实施例还提供一种分段积叠的轴流压气机的转子叶片10,所述转子叶片10通过上述分段积叠的轴流压气机的转子叶片的设计方法得到。
需要说明的是,图6中的最终叶型的转子叶片10,是在图7中的原始叶型的转子叶片10’的基础上通过上述分段积叠的轴流压气机的转子叶片的设计方法得到的。可以从图7中看到,原始叶型的转子叶片10’的积叠线20’为一直线,而本实施例的最后得到的转子叶片10的积叠线20在周向上为S型。
本实施例通过设置转子叶片的积叠线在周向上为S型,可以在径向方向会产生逆压力梯度,有利于叶尖区域的低能流体往叶中迁移,进而减弱了叶尖区域二次流影响,改善了叶尖流动状况,提升了压气机转子叶片的稳定工作范围。进一步,依据损失最小或效率最佳原则确定最优的最大周向偏移量xb和xc及相应高度和延伸长度,并通过设置第一区间中的各基元叶型截面的重心位置点的坐标符合特定的函数,以通过构造简便的方程来控制位于的第一区间中的各基元叶型截面的重心位置的周向偏移量,提高优化叶型的计算效率。也就是说,这相对于现有技术中,通过试错的方法来构造积叠线的方式,大大地提高了计算效率。另一方面,现有技术中用于描述近似“S”型曲线的方程多为单方程形式;而在压气机的积叠叶型优化中,通常涉及多个变量的优化,如最大周向偏移量及其位置,和根尖位置的倾斜角度等,此时单一方程形式无法对压气机周向积叠的关键特征参数进行优化。本方法通过提取压气机周向积叠最为关键的设计参数,最大正、负向周向偏移量及其高度位置,和根尖位置的倾斜角度,并将其作为输入量即自变量,采用分段函数形式对其余高度位置的偏移量进行封闭求解,并保证了分段函数连接点处的光滑过渡,至少一阶导数连续。在保证压气机叶型周向积叠关键参数能够被完全描述的基础上,通过构造形式简单的分段函数,使得方程易于求解,提升了叶型优化效率。也就是说,本发明给出了一种将压气机周向积叠最大正、负向周向偏移量及其高度位置,和根尖位置的倾斜角度作为控制参数的参数化方法,采用分段函数形式对其余高度位置的偏移量进行封闭求解,保证了分段函数连接点处的光滑过渡,以及至少一阶导数连续。
虽然以上描述了本发明的具体实施方式,但是本领域的技术人员应当理解,这仅是举例说明,本发明的保护范围是由所附权利要求书限定的。本领域的技术人员在不背离本发明的原理和实质的前提下,可以对这些实施方式做出多种变更或修改,但这些变更和修改均落入本发明的保护范围。
Claims (10)
1.一种分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,所述转子叶片包括叶根和叶尖,所述转子叶片沿高度方向从所述叶根到所述叶尖具有若干个基元叶型截面,将各所述基元叶型截面的重心作为积叠点沿所述转子叶片的高度方向进行叠加得到积叠线,所述设计方法包括以下步骤:
S1:确定所述积叠线的轴向位置;
S2:确定所述积叠线的周向位置:所述积叠线在周向上为“S”型,所述积叠线沿高度方向由所述叶根至所述叶尖依次包括A、B、C和D四个控制点,A点为所述叶根的基元叶型截面的重心位置,且A点为原点;B点位于所述积叠线在X轴正向上的最大周向偏移量位置点,且B点为预设点;C点位于所述积叠线在X轴负向上的最大周向偏移量位置点,且C点为预设点;D点为所述叶尖的基元叶型截面的重心位置,且D点在X轴上位于B点和C点之间;
B点至C点的第一区间中的各所述基元叶型截面的重心位置点的坐标(x1,y1),符合以下公式:
其中,N=(xb–xc)/2;α=Pi/(yc–yb),Pi=3.1415926;参数和M根据B点的坐标(xb,yb)和C点的坐标(xc,yc)代入公式得到;
S3:根据确定后的所述积叠线上的每一所述基元叶型截面的重心的轴向位置和周向位置,确定所述积叠线的线型;将基元叶型沿所述积叠线的线型进行延伸构造曲线组,得到三维的转子叶片。
2.如权利要求1所述的分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,在步骤S2中,包括:
将所述第一区间在B点沿高度方向向下延伸第一延伸高度Δyb,得到第一中间点E;
A点至E点的第二区间中的各所述基元叶型截面重心的周向位置点的坐标(x2,y2),符合以下公式:
其中,参数h、i、j通过E点的坐标(xe,ye)和E点的斜率ke,以及A点的坐标(xa,ya)代入公式得到。
3.如权利要求2所述的分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,所述第一延伸高度Δyb为所述第一区间的高度的5%~20%。
4.如权利要求1所述的分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,在步骤S2中,包括:
将所述第一区间在C点沿高度方向向上延伸第二延伸高度Δyc,得到中间点F;
F点至D点的第三区间中的各所述基元叶型截面的重心位置点的坐标(x3,y3),符合以下公式:
x3=kfy3+Q
其中,kf为F点的斜率,参数Q通过F点的坐标(xf,yf)和F点的斜率kf代入公式得到。
5.如权利要求4所述的分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,所述第二延伸高度Δyc为所述第一区间的高度的5%~20%。
6.如权利要求1所述的分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,
在步骤S2中,在B点的坐标(xb,yb)中,yb为叶高的25%~35%;和/或,
在步骤S2中,在C点的坐标(xc,yc)中,yc为叶高的70%~85%。
7.如权利要求1所述的分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,
在步骤S1中,所述确定所述积叠线的轴向位置,包括:
在轴向上,将各所述基元叶型截面的重心位置以所述叶根的所述基元叶型截面的重心位置为基准点沿高度方向移至一条直线上;设定叶尖的所述基元叶型截面的重心位置的轴向偏移量为最大轴向偏移量,以确定叶尖的所述基元叶型截面的重心位置;其余所述基元叶型截面的重心位置的轴向偏移量采用线性插值的方式得到。
8.如权利要求7所述的分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,
在步骤S1中,在轴向上所述叶尖的所述基元叶型截面的重心位置比叶根的所述基元叶型截面的重心位置靠前。
9.如权利要求1-8中任意一项所述的分段积叠的轴流压气机的转子叶片的设计方法,其特征在于,在步骤S3之后,还包括:
S4:重复步骤S1至S3多次,获得多个不同的所述转子叶片,获得多个不同的所述转子叶片所对应的压气机效率,选择最高的所述压气机效率所对应的所述转子叶片作为最终的转子叶片。
10.一种分段积叠的轴流压气机的转子叶片,其特征在于,所述转子叶片通过权利要求1-9中任意一项所述的分段积叠的轴流压气机的转子叶片的设计方法得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210787403.XA CN117386664A (zh) | 2022-07-04 | 2022-07-04 | 分段积叠的轴流压气机的转子叶片及其设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210787403.XA CN117386664A (zh) | 2022-07-04 | 2022-07-04 | 分段积叠的轴流压气机的转子叶片及其设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117386664A true CN117386664A (zh) | 2024-01-12 |
Family
ID=89467281
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210787403.XA Pending CN117386664A (zh) | 2022-07-04 | 2022-07-04 | 分段积叠的轴流压气机的转子叶片及其设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117386664A (zh) |
-
2022
- 2022-07-04 CN CN202210787403.XA patent/CN117386664A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112528575B (zh) | 一种船舶燃气轮机多级轴流压气机多维度协同设计方法 | |
US7547186B2 (en) | Nonlinearly stacked low noise turbofan stator | |
US10519980B2 (en) | Turbomachine component or collection of components and associated turbomachine | |
CN114444196B (zh) | 一种自适应缝隙调节可变弯度导叶设计方法及其装置 | |
US10352330B2 (en) | Turbomachine part with a non-axisymmetric surface | |
CN112685852B (zh) | 一种轴流压气机负荷定制化气动优化方法 | |
CN108829970A (zh) | 基于涡动力学的轴流风扇叶片优化设计方法及优化系统 | |
CN105332948A (zh) | 改进型压气机叶片及其实现方法 | |
US10626728B2 (en) | Blade for a turbomachine impeller and method for modelling the same | |
CN114186513A (zh) | 一种具有反s型前缘的轴流压气机叶片造型设计方法 | |
CN109815590B (zh) | 一种基于端区附面层的多级轴流压气机三维叶片造型方法及叶片 | |
CN116090129A (zh) | 耦合非轴对称端壁和叶片弯曲的涡轮造型方法 | |
CN110030038B (zh) | 考虑bli进气畸变效应的叶尖跨音风扇非对称静子设计方法 | |
CN214247795U (zh) | 一种带有鼓包和凹缝结构的跨声速压气机转子叶片 | |
CN101158991A (zh) | 大小叶片压气机的气动设计方法 | |
CN109779971B (zh) | 基于曲率控制的高负荷压气机叶型径向积叠造型优化方法 | |
CN117386664A (zh) | 分段积叠的轴流压气机的转子叶片及其设计方法 | |
CN115270318B (zh) | 一种船用燃气轮机轴流压气机跨音级动叶片造型方法 | |
CN116579088A (zh) | 一种基于贝塞尔曲线的压气机静子叶片仿生前缘设计方法 | |
CN113569498A (zh) | 一种轴流压气机端部弯曲静叶片设计方法 | |
CN105298912A (zh) | 鼓包前缘进口导向器叶片 | |
CN112446107A (zh) | 一种压气机端区流动控制建构的建立方法 | |
CN117386663A (zh) | 轴流压气机的转子叶片及其积叠线的获取方法 | |
CN112648235B (zh) | 一种带有鼓包和凹缝结构的跨声速压气机转子叶片 | |
CN114046269B (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 |