CN110298052B - 串列叶片优化实现方法 - Google Patents

串列叶片优化实现方法 Download PDF

Info

Publication number
CN110298052B
CN110298052B CN201810902368.5A CN201810902368A CN110298052B CN 110298052 B CN110298052 B CN 110298052B CN 201810902368 A CN201810902368 A CN 201810902368A CN 110298052 B CN110298052 B CN 110298052B
Authority
CN
China
Prior art keywords
grid
blade
blades
points
row
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
CN201810902368.5A
Other languages
English (en)
Other versions
CN110298052A (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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201810902368.5A priority Critical patent/CN110298052B/zh
Publication of CN110298052A publication Critical patent/CN110298052A/zh
Application granted granted Critical
Publication of CN110298052B publication Critical patent/CN110298052B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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

Abstract

一种串列叶片模型优化实现方法,根据选定的串列静子叶片的轴向与周向相对位置的参数变量,在原始几何模型的基础上计算出更新几何模型,在更新几何模型上对前后排叶片在不同径向叶高截面处的轮廓点的集合,逐层逐点检测设计参数是否合理,最后计算通过检测的设计参数与若干参考样本的差异,选取差异最小的样本所对应的数值模拟结构化网格点的分布方案用于生成网格;本发明设计合理,能够提高串列叶片的相对位置参数的设计效率,降低出错率,减少人工耗时,缩短设计周期并降低成本,适合二次开发和拓展应用。

Description

串列叶片优化实现方法
技术领域
本发明涉及的是一种适用于叶轮机械领域的技术,具体是一种关于串列叶片的几何重合检测与数值模拟网格自动生成方法。
背景技术
串列叶片技术能够在不增加太多重量和整体长度的前提下,提高航空发动机的工作负荷,提高气流的折转角度,进而提高航空发动机压气机的气动载荷,在一些航空发动机的型号中得到了有效的应用。
在设计串列叶片时,前后排叶片的相对位置对流场和压气机效率、压比、稳定工作裕度等性能有很大的影响。因此在串列叶片的设计过程中,找到一个合适的相对位置是十分重要的。但在传统的设计方法里,所需的考察的设计样本数量较大,逐一人工检查几何是否发生重合效率低下,容易出错。此外,较大数量的样本也增加了手动绘制网格的工作量。因此,几何重合的检测以及网格自动化生成的需求十分迫切。
发明内容
本发明针对现有技术存在的上述不足,提出一种串列叶片优化实现方法,通过对给定的串列叶片轴向和周向的位置参数进行数学几何分析,求解出这种设计参数是否发生了几何重合的状况,并通过现有的商业计算流体力学软件的造型模块,利用其内置代码指令,实现了较高质量的数值模拟网格的自动化生成。
本发明是通过以下技术方案实现的:
本发明根据选定的串列静子叶片的轴向与周向相对位置的参数变量,在原始几何模型的基础上计算出更新几何模型,在更新几何模型上对前后排叶片在不同径向叶高截面处的轮廓点的集合,逐层逐点检测设计参数是否合理,最后计算通过检测的设计参数与若干参考样本的差异,选取差异最小的样本所对应的数值模拟结构化网格点的分布方案用于生成网格。
所述的检测设计参数是否合理,即前排叶片的轮廓点是否在后排叶片的轮廓点的顺序连线所构成的封闭区域外。
所述的参数变量为:后排叶片相对前排叶片的轴向位置和周向位置,其坐标系设定为:以叶轮机械转动中心轴为轴向,设定叶轮机械出口至叶轮机械入口的轴向为正,轴向位置的单位为米。环绕正向轴向,右手旋转的方向为周向正向方向,周向位置的单位为弧度。叶片从叶根到叶尖的方向为径向位置正方向,单位为米;以后排叶片的轮毂截面上的前缘点和前排叶片的轮毂截面上的尾缘点为基准,确定前后两排叶片的相对位置。
技术效果
与现有技术相比,本发明提高了航空发动机串列叶片相对位置设计的效率,提高了自动化程度,减少了人力成本与工作时间,降低了出错的概率。基于上述方法开发的程序,具有良好的扩展性,可以对数量很大的样本进行批量处理,适合二次开发与迁移。
附图说明
图1为本发明流程示意图;
图2为串列叶片相对位置变量的示意图;
图中:1串列叶片前排、2串列叶片后排、ΔZ串列叶片轴向位置、Δθ串列叶片周向位置、
Figure BDA0001759742140000021
相对轴向位置、
Figure BDA0001759742140000022
相对周向位置、ZFrontBlade是前排叶片轴向长度、θBlade串列叶片单个通道角度;
图3为串列叶片发生几何重合时的相对位置示意图;
图4为串列叶片几何重合检测的结果示意图。
具体实施方式
如图1所示,本实施例为轴流压气机的一级串列静子叶片的几何相对位置的设计与优化过程中,对几何重合的检测以及随后的网格自动生成流程,具体包括以下过程:
步骤1、根据选定的串列静子叶片的轴向与周向相对位置的参数变量,在原始几何模型的基础上计算出更新几何模型。
如图2所示,所述的参数变量为以下任意一种组合:
①轴向位置ΔZ和周向位置Δθ,或
②相对轴向位置
Figure BDA0001759742140000023
和相对周向位置
Figure BDA0001759742140000024
所述的原始几何模型是指:叶片的数个径向高度的截面内的轮廓点的坐标的集合,其坐标系是以叶轮机械转轴为中心轴的圆柱坐标系。
步骤2、对前后排叶片在不同径向叶高截面处的轮廓点的集合,逐层逐点检测设计参数是否合理,即前排叶片的轮廓点是否在后排叶片的轮廓点的顺序连线所构成的封闭区域内。
由于本实施例中为直叶片,故叶片不同径向叶高的截面为该串列叶片造型时所需要的控制面为轮毂和机匣截面上的两个控制面。
所述的逐层逐点检测是指:以前排叶片上任意一点作为被检测的射线起点,向任意方向做一射线,计算此射线与后排叶片轮廓点顺序连线所构成的所有线段的相交次数:当相交次数为偶数时被检测点在后排叶片轮廓点顺序连线的封闭区域之外,否则被检测点在后排叶片轮廓点顺序连线的封闭区域之内;当存在任一截面内的任一前排叶片轮廓点位于后排叶片的轮廓区域内,即前、后排叶片存在几何重合时,该设计参数判定为不合理,应重新选取设计参数,当所有截面内的所有前排叶片轮廓点都位于后排叶片的轮廓区域外,即前、后排叶片不存在几何重合时,该设计参数判定为合理。
如图4所示,串列叶片相对位置的设计参数空间经过几何重合的检查,划分为发生几何重合的参数区间(白色区域),和不发生几何重合的参数区间(黑色区域),唯有在该区域内所选取的参数是合理的(白色圆圈)。
相比于人工目视检测,采用这种几何重合自动检测的方法可以迅速判断某个设计参数的可行性,从而无需在三维几何可视化软件中目视逐个判断所有样本,实现了该过程的全自动化。此外对于单个样本,人工目视一般需要数分钟的工作时间,而该方法的计算机耗时在数秒之内。
步骤3、对于通过检测的设计参数,计算其与若干参考样本的差异,选取差异最小的样本所对应的数值模拟结构化网格点的分布方案用于生成网格。
所述的参考样本为设计参数的变化空间内预先均匀给定的样本以及对应样本的数值模拟结构化网格点的分布方案。
所述的差异,即设计参数与参考样本的参数之间的差值所对应的范数。
所述的范数采用但不限于平方范数。
所述的数值模拟结构化网格点的分布方案为O4H型结构化网格,包括:
1)叶片上方H型网格进口截面周向网格点数;
2)叶片上方H型网格出口截面周向网格点数;
3)叶片下方H型网格进口截面周向网格点数;
4)叶片下方H型网格出口截面周向网格点数;
5)叶片前方H型网格进口截面周向网格点数;
6)叶片后方H型网格出口截面周向网格点数;
7)叶片前方H型网格轴向网格点数;
8)叶片后方H型网格轴向网格点数;
9)叶片贴体O型网格吸力面网格点数;
10)叶片贴体O型网格压力面网格点数;
11)叶片贴体O型网格厚度方向网格点数。
所述的生成网格是指:将设计参数对应的数值模拟结构化网格点的分布方案采用但不限于流体力学软件Numeca中的网格绘制模块Autogrid所支持的指令生成,并写入代码文件中执行后得到网格。
相比于人工绘制网格,该网格自动生成方法实现了网格生成的自动化。通常对于单个样本,人工绘制网格的参数设置需要数分钟,并且需要等待数十分钟网格生成完毕之后,进行保存网格等后续操作,全程需要人工参与。而采用新的网格自动生成方法之后,仅需启动程序,交由计算机完成所有步骤。
所述的生成网格,优先对多个设计参数进行几何重合检测,并将各个参数对应的设计方案的网格生成指令写入一个代码文件中并执行,则可计算获得所有的设计参数所对应的网格文件,从而适用于数值模拟计算的过程。
上述具体实施方式的应用环境为压气机串列静子叶片设计和优化的过程,具体操作中所涉及的参数包括串列静子前后排叶片的轴向位置和周向位置两个参数。使用该方法可以对给定的串列叶片相对位置的设计参数进行几何重合的检测,并自动化生成对应的数值模拟网格。与现有的手动生成网格,并目视观察是否发生几何重合的操作相比,本发明可以大大提高工作效率,降低出错概率,并且能够实现大量样本的批量化处理,大大提高了程序的自动化程度。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。

Claims (5)

1.一种串列叶片模型优化实现方法,其特征在于,根据选定的串列静子叶片的轴向与周向相对位置的参数变量,在原始几何模型的基础上计算出更新几何模型,在更新几何模型上对前后排叶片在不同径向叶高截面处的轮廓点的集合,逐层逐点检测设计参数是否合理,最后计算通过检测的设计参数与若干参考样本的差异,选取差异最小的样本所对应的数值模拟结构化网格点的分布方案用于生成网格;
所述的参数变量为以下任意一种组合:
①轴向位置ΔZ和周向位置Δθ,或
②相对轴向位置
Figure FDA0003924105550000011
和相对周向位置
Figure FDA0003924105550000012
所述的原始几何模型是指:叶片的数个径向高度的截面内的轮廓点的坐标的集合,其坐标系是以叶轮机械转轴为中心轴的圆柱坐标系;
所述的逐层逐点检测是指:以前排叶片上任意一点作为被检测的射线起点,向任意方向做一射线,计算此射线与后排叶片轮廓点顺序连线所构成的所有线段的相交次数:当相交次数为偶数时被检测点在后排叶片轮廓点顺序连线的封闭区域之外,否则被检测点在后排叶片轮廓点顺序连线的封闭区域之内;当存在任一截面内的任一前排叶片轮廓点位于后排叶片的轮廓区域内,即前、后排叶片存在几何重合时,该设计参数判定为不合理,应重新选取设计参数,当所有截面内的所有前排叶片轮廓点都位于后排叶片的轮廓区域外,即前、后排叶片不存在几何重合时,该设计参数判定为合理;
所述的生成网格是指:将设计参数对应的数值模拟结构化网格点的分布方案采用但不限于流体力学软件Numeca中的网格绘制模块Autogrid所支持的指令生成,并写入代码文件中执行后得到网格;
所述的生成网格,对多个设计参数进行几何重合检测,并将各个参数对应的设计方案的网格生成指令写入一个代码文件中并执行,则可计算获得所有的设计参数所对应的网格文件,从而适用于数值模拟计算的过程。
2.根据权利要求1所述的方法,其特征是,所述的检测设计参数是否合理,即前排叶片的轮廓点是否在后排叶片的轮廓点的顺序连线所构成的封闭区域外。
3.根据权利要求1所述的方法,其特征是,所述的参考样本为设计参数的变化空间内预先均匀给定的样本以及对应样本的数值模拟结构化网格点的分布方案。
4.根据权利要求1所述的方法,其特征是,所述的差异,即设计参数与参考样本的参数之间的差值所对应的范数。
5.根据权利要求1所述的方法,其特征是,所述的数值模拟结构化网格点的分布方案为O4H型结构化网格,包括:
1)叶片上方H型网格进口截面周向网格点数;
2)叶片上方H型网格出口截面周向网格点数;
3)叶片下方H型网格进口截面周向网格点数;
4)叶片下方H型网格出口截面周向网格点数;
5)叶片前方H型网格进口截面周向网格点数;
6)叶片后方H型网格出口截面周向网格点数;
7)叶片前方H型网格轴向网格点数;
8)叶片后方H型网格轴向网格点数;
9)叶片贴体O型网格吸力面网格点数;
10)叶片贴体O型网格压力面网格点数;
11)叶片贴体O型网格厚度方向网格点数。
CN201810902368.5A 2018-08-09 2018-08-09 串列叶片优化实现方法 Active CN110298052B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810902368.5A CN110298052B (zh) 2018-08-09 2018-08-09 串列叶片优化实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810902368.5A CN110298052B (zh) 2018-08-09 2018-08-09 串列叶片优化实现方法

Publications (2)

Publication Number Publication Date
CN110298052A CN110298052A (zh) 2019-10-01
CN110298052B true CN110298052B (zh) 2022-12-23

Family

ID=68026292

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810902368.5A Active CN110298052B (zh) 2018-08-09 2018-08-09 串列叶片优化实现方法

Country Status (1)

Country Link
CN (1) CN110298052B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113090582B (zh) * 2021-03-30 2022-04-01 南京工业大学 一种基于前、后缘修形产生非定常射流的串列叶片
CN113309737B (zh) * 2021-06-22 2023-08-11 西北工业大学 压气机三段式可调串列叶片

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009013399A1 (de) * 2009-03-16 2010-09-23 Mtu Aero Engines Gmbh Tandemschaufelkonstruktion
CN104317997A (zh) * 2014-10-17 2015-01-28 北京航空航天大学 一种高负荷风扇/压气机端壁造型优化设计方法
CN105257590B (zh) * 2015-10-20 2017-12-15 南京航空航天大学 半串列叶片及其设计方法
CN107679319B (zh) * 2017-09-29 2020-07-07 北京航空航天大学 一种叶轮机通流模型中周向脉动应力项的代数建模方法

Also Published As

Publication number Publication date
CN110298052A (zh) 2019-10-01

Similar Documents

Publication Publication Date Title
US8881392B2 (en) Method of repairing machined components such as turbomachine blades or blades of blisks
CN110298052B (zh) 串列叶片优化实现方法
Scharfenstein et al. Probabilistic CFD analysis of high pressure turbine blades considering real geometric effects
CN114254460A (zh) 基于图卷积神经网络的透平机械气动鲁棒性优化方法
US20040030666A1 (en) Method of designing a multi-stage compressor rotor
EP3617994A1 (en) Blade sentencing
Shahpar Building digital twins to simulate manufacturing variation
CN110276115B (zh) 基于燃机叶片型线参数的气路故障诊断方法
EP4258074A1 (en) Systems and methods for aerodynamic analysis for inspected bladed rotors
Zhou et al. Parametric modeling method for integrated design and manufacturing of radial compressor impeller
US20230315937A1 (en) Partial repair systems and methods for integrally bladed rotors
CN111353249B (zh) 涡轮封严盘非圆通气孔集成设计优化方法
Meckstroth et al. Point Cloud to Parameter: An Inverse Geometric Approach to Probabilistic Design
Kaszynski et al. Uncertainties of an automated optical 3D geometry measurement, modeling, and analysis process for mistuned ibr reverse engineering
Arkhipov et al. Creation of models of fan blades according to airfoil measurements
CN117010093B (zh) 一种汽车造型设计方法及设备
EP4257800A1 (en) Inspection data management systems and methods
EP4257803A1 (en) Data transfer systems and methods for providing repair analysis inputs
US20230313685A1 (en) Bladed rotor inspection, analysis and repair systems and methods
EP4257802A1 (en) Structural repair analysis system and method for inspected bladed rotors
US20230311252A1 (en) Blend approach based inspection and analysis systems and methods
CN117217013A (zh) 一种涡轮叶片气膜孔综合设计方法及系统
CN117454721B (zh) 基于数字仿真实验的风电场尾流叠加效应评估方法和介质
US20220048145A1 (en) Turbine airfoil design
US20230315039A1 (en) Stack based repair systems and methods

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