CN106709179A - 一种适用于有限元自由网格的叶片气动力分布数据处理方法 - Google Patents
一种适用于有限元自由网格的叶片气动力分布数据处理方法 Download PDFInfo
- Publication number
- CN106709179A CN106709179A CN201611193176.9A CN201611193176A CN106709179A CN 106709179 A CN106709179 A CN 106709179A CN 201611193176 A CN201611193176 A CN 201611193176A CN 106709179 A CN106709179 A CN 106709179A
- Authority
- CN
- China
- Prior art keywords
- module
- aerodynamic
- simulation lattice
- unit
- blade
- 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.)
- Granted
Links
Classifications
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及数据处理方法,具体涉及一种适用于有限元自由网格的航空发动机叶片气动力分布数据处理方法,在进行航空发动机叶片强度设计时,实现了方便快捷的叶片有限元自由网格气动力数据获取及处理,结果可靠,所用时间少。
Description
技术领域
本发明涉及数据处理方法,具体涉及一种适用于有限元自由网格的航空发动机叶片气动力分布数据处理方法。
背景技术
叶片是航空发动机的主要功能零件,用于实现整流、增压、膨胀做功等功能,由于气体作用力、离心力及温度载荷的综合作用,叶片的强度、振动及疲劳寿命问题比较突出,其中气动力的表达对强度和疲劳寿命分析结果具有了决定性的影响,因此需要尽可能表述准确。
叶片所受气动力依托流体力学方法,建立流道控制体,根据动量守恒定律及牛顿第三定律推导出数学表达式。
由于叶片结构形式的多样性及复杂性,在进行强度和寿命分析时通常采用有限元的方法。网格的多样性使叶片气动力的有限元表达增加了复杂性。
一般在处理四边形有限元网格的气动力分布计算时,可采用将控制体的气动力按节点总数平均分配到各节点上的方法进行简化处理。
但在工程设计中,叶片的有限元网格形式自由,无法形成规则网格,导致无法采用按节点总数平均分配的简化处理方法加载气动力。
发明内容
本发明的目的是针对叶片气动力平均分配数据处理方法不适用于自由网格的数据处理问题,提供一种适用于有限元自由网格的叶片气动力分布数据处理方法。
为了实现所述目的,本发明采取了如下技术方案:
一种适用于有限元自由网格的叶片气动力分布数据处理方法,包括以下步骤:
步骤1原始参数载入模块获取气动数据截面半径,获取速度、温度、密度、压力,同时获取叶片表面单元信息、叶片表面节点坐标、进排气边节点坐标,并将获取的单元信息发送至单元类型适应性转化模块、边界曲线拟合模块和单元形函数计算模块、模拟网格模块;
步骤2单元类型适应性转化模块,将步骤1中获取的单元信息转化为以三角形单元为基础的统一信息存储结构和计算标准,并将转化后的信息发送至单元形函数计算模块;
步骤3边界曲线拟合模块,将步骤1中获取的进、排气边节点坐标拟合出边界曲线;
步骤4单元形函数计算模块,根据步骤1中获取的叶片表面节点坐标和步骤2转化后的单元信息,计算自由网格的单元形函数及最小单元长度;
步骤5模拟网格构建,模拟网格模块根据步骤4确定的自由网格最小单元长度,结合步骤3拟合的边界曲线和步骤1获取的气动数据截面半径等参数,构建模拟网格;
步骤6模拟网格气动力计算,模拟网格模块以步骤5构建的模拟网格为基础,由步骤1中获取的气动数据计算出模拟网格的气动力;
步骤7自由网格气动力分配模块将步骤6计算得到的模拟网格气动力,按照步骤4计算到的单元形函数进行分配;
步骤8气动力计算结果格式化输出模块将步骤7完成的自由网格气动力分配数据格式化输出至存储模块和显示模块。
本发明的有益效果:
在进行航空发动机叶片强度设计时,实现了方便快捷的叶片有限元自由网格气动力数据获取及处理,结果可靠,所用时间少。
附图说明
图1是本发明的流程示意图
具体实施方式
一种适用于有限元自由网格的叶片气动力分布数据处理方法,包括原始参数载入模块、单元类型适应性转化模块、边界曲线拟合模块、单元形函数计算模块、模拟网格模块、自由网格气动力分配模块、气动力计算结果格式化输出模块。
所述的原始参数载入模块,其作用在于:将叶片气动力计算所涉及的截面半径、气流环量或周向速度、气流轴向速度、温度、气体密度、静压、转速等S2流场数据自动输入并暂存,同时将叶片表面单元信息、叶片表面节点坐标、进排气边节点坐标等自由网格信息自动输入并形成数据链表暂存。
所述的单元类型适应性转化模块,其作用在于:将自由网格可能涉及的8节点四边形面网格、4节点四边形面网格和6节点三角形面网格统一转化为3节点三角形面网格,形成统一的信息存储结构和计算标准。
所述的边界曲线拟合模块,其作用在于:将叶片进、排气边的节点坐标分别进行多项式曲线拟合,根据节点数量自动判定多项式阶数,最高不超过20阶。
所述的单元形函数计算模块,其作用在于:根据单元数据链表中的节点坐标信息,将自由网格单元中各节点对应的形函数参数分别按照3节点三角形单元或6节点等参数三角形单元计算获取并在数据链表中暂存,计算自由网格的最小单元长度。
所述的模拟网格模块,其作用在于:取自由网格最小单元长度的一半作为模拟网格的单元长度,同时设置防止单元长度过小的限制措施,根据气动截面半径、自由网格几何边界及模拟网格的单元长度确定规则化的模拟网格数量并构建模拟网格;根据动量守恒定律及牛顿第三定律计算得到模拟网格对应各流道控制体的气动力,并通过平均分配的原则将其分配至对应的模拟网格节点上。
所述的自由网格气动力分配模块,其作用在于:搜索模拟网格中的节点坐标所在的自由网格单元,并通过单元的形函数将模拟网格节点力转换到自由网格节点上。
所述的气动力计算结果格式化输出模块,其作用在于:将自由网格节点力按照有限元软件ANSYS的APDL程序格式输出至存储模块和显示模块。
Claims (1)
1.一种适用于有限元自由网格的叶片气动力分布数据处理方法,其特征在于,包括以下步骤:
步骤1原始参数载入模块获取气动数据截面半径,获取速度、温度、密度、压力,同时获取叶片表面单元信息、叶片表面节点坐标、进排气边节点坐标,并将获取的单元信息发送至单元类型适应性转化模块、边界曲线拟合模块和单元形函数计算模块、模拟网格模块;
步骤2单元类型适应性转化模块,将步骤1中获取的单元信息转化为以三角形单元为基础的统一信息存储结构和计算标准,并将转化后的信息发送至单元形函数计算模块;
步骤3边界曲线拟合模块,将步骤1中获取的进、排气边节点坐标拟合出边界曲线;
步骤4单元形函数计算模块,根据步骤1中获取的叶片表面节点坐标和步骤2转化后的单元信息,计算自由网格的单元形函数及最小单元长度;
步骤5模拟网格构建,模拟网格模块根据步骤4确定的自由网格最小单元长度,结合步骤3拟合的边界曲线和步骤1获取的气动数据截面半径等参数,构建模拟网格;
步骤6模拟网格气动力计算,模拟网格模块以步骤5构建的模拟网格为基础,由步骤1中获取的气动数据计算出模拟网格的气动力;
步骤7自由网格气动力分配模块将步骤6计算得到的模拟网格气动力,按照步骤4计算到的单元形函数进行分配;
步骤8气动力计算结果格式化输出模块将步骤7完成的自由网格气动力分配数据格式化输出至存储模块和显示模块。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611193176.9A CN106709179B (zh) | 2016-12-21 | 2016-12-21 | 一种叶片气动力分布数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611193176.9A CN106709179B (zh) | 2016-12-21 | 2016-12-21 | 一种叶片气动力分布数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106709179A true CN106709179A (zh) | 2017-05-24 |
CN106709179B CN106709179B (zh) | 2020-04-07 |
Family
ID=58938674
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611193176.9A Active CN106709179B (zh) | 2016-12-21 | 2016-12-21 | 一种叶片气动力分布数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106709179B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108090292A (zh) * | 2017-12-26 | 2018-05-29 | 中国航发四川燃气涡轮研究院 | 一种宽弦风扇叶片二维有限元建模方法 |
CN110162823A (zh) * | 2019-03-19 | 2019-08-23 | 北京机电工程研究所 | 考虑气动面曲面效应和法向运动的非定常气动力计算方法 |
CN111539082A (zh) * | 2020-04-26 | 2020-08-14 | 重庆水轮机厂有限责任公司 | 一种水轮机叶片三维模型快速重建方法 |
CN115688314A (zh) * | 2022-11-03 | 2023-02-03 | 北京全四维动力科技有限公司 | 基于有限元分析的汽轮机叶片预扭仿真设计方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102479269A (zh) * | 2010-11-29 | 2012-05-30 | 北京临近空间飞行器系统工程研究所 | 一种流-固载荷的转换方法 |
CN105787199A (zh) * | 2015-07-06 | 2016-07-20 | 中航空天发动机研究院有限公司 | 一种叶轮机转子叶片从热态到冷态的叶型转换方法 |
CN106021827A (zh) * | 2016-07-14 | 2016-10-12 | 湖北工业大学 | 一种考虑气动载荷的风力机叶片内部结构拓扑设计方法 |
-
2016
- 2016-12-21 CN CN201611193176.9A patent/CN106709179B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102479269A (zh) * | 2010-11-29 | 2012-05-30 | 北京临近空间飞行器系统工程研究所 | 一种流-固载荷的转换方法 |
CN105787199A (zh) * | 2015-07-06 | 2016-07-20 | 中航空天发动机研究院有限公司 | 一种叶轮机转子叶片从热态到冷态的叶型转换方法 |
CN106021827A (zh) * | 2016-07-14 | 2016-10-12 | 湖北工业大学 | 一种考虑气动载荷的风力机叶片内部结构拓扑设计方法 |
Non-Patent Citations (1)
Title |
---|
刘良玉: "叶片三维有限元数据处理方法研究", 《中国优秀博硕士学位论文全文数据库 (硕士) 工程科技Ⅱ辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108090292A (zh) * | 2017-12-26 | 2018-05-29 | 中国航发四川燃气涡轮研究院 | 一种宽弦风扇叶片二维有限元建模方法 |
CN108090292B (zh) * | 2017-12-26 | 2021-08-03 | 中国航发四川燃气涡轮研究院 | 一种宽弦风扇叶片二维有限元建模方法 |
CN110162823A (zh) * | 2019-03-19 | 2019-08-23 | 北京机电工程研究所 | 考虑气动面曲面效应和法向运动的非定常气动力计算方法 |
CN111539082A (zh) * | 2020-04-26 | 2020-08-14 | 重庆水轮机厂有限责任公司 | 一种水轮机叶片三维模型快速重建方法 |
CN115688314A (zh) * | 2022-11-03 | 2023-02-03 | 北京全四维动力科技有限公司 | 基于有限元分析的汽轮机叶片预扭仿真设计方法及装置 |
CN115688314B (zh) * | 2022-11-03 | 2024-01-23 | 北京全四维动力科技有限公司 | 基于有限元分析的汽轮机叶片预扭仿真设计方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106709179B (zh) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106709179A (zh) | 一种适用于有限元自由网格的叶片气动力分布数据处理方法 | |
Ferrer et al. | : A high-order discontinuous Galerkin solver for flow simulations and multi-physics applications | |
KR101612506B1 (ko) | 전산 유체 역학을 이용한 항공기 공력해석 시스템 및 방법 | |
US7941300B1 (en) | Process for the design of an airfoil | |
Restemeier et al. | Numerical and experimental analysis of the effect of variable blade row spacing in a subsonic axial turbine | |
Pinelli et al. | On the numerical evaluation of tone noise emissions generated by a turbine stage: an in-depth comparison among different computational methods | |
Bontempo et al. | A ring-vortex free-wake model for uniformly loaded propellers. Part I-Model description. | |
Lengyel et al. | Design of a counter rotating fan-an aircraft engine technology to reduce noise and CO2-emissions | |
Barton et al. | Numerical optimization of a vaned shroud design for increased operability margin in modern centrifugal compressors | |
JP6684330B2 (ja) | インテークダクト設計方法、インテークダクト設計プログラム及びインテークダクト設計装置 | |
Akaydin et al. | Implementation of a body force model in OVERFLOW for propulsor simulations | |
Xu et al. | A study of Nash-EGO algorithm for aerodynamic shape design optimizations | |
Patrao et al. | An optimization platform for high speed propellers | |
Yanhui et al. | An introduction of aerodynamic shape optimization platform for compressor blade | |
Park et al. | Optimization of a 3-Stage Booster: Part 1—The Axisymmetric Multi-Disciplinary Optimization Approach to Compressor Design | |
Anisimov et al. | Aerodynamic optimization of airplane propulsion system within the framework of AGILE project | |
Verstraete | Toward gradient-based optimization of full gas turbines | |
Anisimov et al. | Propulsion system-airframe integration and optimization of civil aircraft-AGILE EU project | |
Dawes | Multi-blade row navier-stokes simulations of fan-bypass configurations | |
CN113987694A (zh) | 一种基于空间推进算法的旋转爆震发动机流场参数分布预测方法 | |
CN107784188A (zh) | 一种基于matlab的增压器压气机叶轮设计优化方法 | |
CN113935114A (zh) | 一种用于水上飞机着水数值仿真的规则波造波方法 | |
Lotz | Aerodynamic Optimization Process for Turbocharger Compressor Impellers | |
Xu et al. | Implicit and conservative mixing-plane method for multistage turbomachinery aerodynamic analysis | |
Wu et al. | Design Optimization of Turbo-Machinery Components With Independent FEA and CFD Tools in an Optimization Software Environment: A Mid-Span Shroud Ring Study Case |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 610500 Xindu Xuefu Road, Xindu District, Chengdu, Sichuan Patentee after: AECC SICHUAN GAS TURBINE Research Institute Address before: 621703 P.O. Box 305, Mianyang City, Sichuan Province Patentee before: CHINA GAS TURBINE EST |