CN116052821A - 一种基于温度影响的轮盘材料利用系数确定方法 - Google Patents

一种基于温度影响的轮盘材料利用系数确定方法 Download PDF

Info

Publication number
CN116052821A
CN116052821A CN202310340757.4A CN202310340757A CN116052821A CN 116052821 A CN116052821 A CN 116052821A CN 202310340757 A CN202310340757 A CN 202310340757A CN 116052821 A CN116052821 A CN 116052821A
Authority
CN
China
Prior art keywords
wheel disc
wheel
utilization coefficient
material utilization
different temperatures
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
Application number
CN202310340757.4A
Other languages
English (en)
Other versions
CN116052821B (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.)
AECC Sichuan Gas Turbine Research Institute
Original Assignee
AECC Sichuan Gas Turbine Research Institute
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 AECC Sichuan Gas Turbine Research Institute filed Critical AECC Sichuan Gas Turbine Research Institute
Priority to CN202310340757.4A priority Critical patent/CN116052821B/zh
Publication of CN116052821A publication Critical patent/CN116052821A/zh
Application granted granted Critical
Publication of CN116052821B publication Critical patent/CN116052821B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • 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
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及航空发动机技术领域,公开了一种基于温度影响的轮盘材料利用系数确定方法,以不同温度下的平均周向应力、轮盘材料拉伸强度以及对应的实际破裂转速值进行不同温度下轮盘的初始材料利用系数值计算,采用公式
Figure ZY_1
进行曲线拟合得到系数
Figure ZY_2
Figure ZY_3
的值,最终得到轮盘的材料利用系数
Figure ZY_4
的表达式。本发明轮盘材料利用系数的确定方法综合考虑了轮盘材料性能的差异的影响,基于轮盘材料性能数据和轮盘破裂试验结果建立了材料利用系数计算的关系式,利用该轮盘材料利用系数进行轮盘设计时,可提高轮盘强度设计精度,保障轮盘设计的可靠性及安全性。

Description

一种基于温度影响的轮盘材料利用系数确定方法
技术领域
本发明涉及航空发动机技术领域,公开了一种基于温度影响的轮盘材料利用系数确定方法。
背景技术
破裂是轮盘的一种常见故障模式,一旦出现总是灾难性的,因此为保证轮盘在规定超转转速下的变形在允许范围内、破裂转速下不破坏,必须进行破裂转速预估和破裂试验验证。国军标《航空涡轮喷气和涡轮风扇发动机通用规范》中要求:当轮盘材料承受最大温度梯度和最高工作温度时,轮盘的设计破裂转速不低于最高允许稳态转速的122%。
目前工程应用中周向破裂转速预测广泛使用平均应力法,即当轮盘截面的平均周向应力与材料拉伸强度(单轴拉伸应力)相等时轮盘就会破裂,但该方法的基本假设是轮盘为理想塑性材料,且应力沿盘的厚度不变。为考虑轮盘实际材料与理想塑性材料间应力应变关系的差异(如图2),及几何形状的影响,引入了材料利用系数这一修正参数,通常轮盘强度设计时取经验系数0.9,导致破裂转速预测误差大,不能满足工程设计。因此可靠的轮盘材料拉伸强度和材料利用系数是准确预测破裂转速的关键因素,有必要研究和建立基于大量材料实测数据与轮盘破裂试验数据的材料利用系数确定方法,以提高设计精度,保障轮盘设计的可靠性及安全性。
发明内容
本发明的目的在于提供一种基于温度影响的轮盘材料利用系数确定方法,可提高轮盘强度设计精度,保障轮盘设计的可靠性及安全性。
为了实现上述技术效果,本发明采用的技术方案是:
一种基于温度影响的轮盘材料利用系数确定方法,包括:
分析获取同一材质不同设计尺寸的轮盘在不同温度下应力分布,计算对应温度下各轮盘平均周向应力;
开展或仿真轮盘旋转破裂试验,获得各轮盘不同温度下的实际破裂转速值;
根据不同温度下的平均周向应力、轮盘材料拉伸强度以及对应的实际破裂转速值,计算不同温度下轮盘的初始材料利用系数值;
根据不同温度下轮盘材料拉伸强度、屈服强度性能数据,分析得到不同温度的轮盘材料拉伸强度置信度为50%的-3,以及屈服强度置信度为50%的-3
采用公式拟合得到轮盘的材料利用系数表达式,其中均为常数;
根据拟合得到轮盘的材料利用系数表达式以及轮盘设计工作温度,计算设计工作温度下的轮盘材料利用系数。
进一步地,计算对应温度下各轮盘平均周向应力,包括:
建立二维轮盘有限元计算模型,并对轮盘有限元计算模型施加转速、温度和轴向位移约束,分析得到轮盘应力分布,并根据轮盘应力分布求解得到平均周向应力,其中为有限元单元的平均周向应力,为轮盘子午面上第个有限元单元面积,为轮盘子午面上有限元单元总数。
进一步地,根据公式计算得到不同温度下轮盘的初始材料利用系数值;其中为材料利用系数,为轮盘实际破裂转速,为轮盘试验基准转速,为轮盘平均周向应力,为轮盘材料在轮盘按面积加权平均温度下的拉伸强度极限值。
进一步地,不同温度的轮盘材料拉伸强度置信度为50%的-3是根据公式计算获得;
其中为多个温度下的轮盘材料在轮盘按面积加权平均温度下的拉伸强度极限值的样本平均值,为样本标准差,为单侧容限系数,为存活率为99.87的标准正态偏量,取值为-3;为置信度为50%的t分布百分位值,为子样数。
与现有技术相比,本发明所具备的有益效果是:
本发明通过综合考虑轮盘材料性能的差异的影响,基于轮盘材料性能数据和轮盘破裂试验结果建立了材料利用系数计算的关系式,利用该轮盘材料利用系数进行轮盘设计时,可提高轮盘强度设计精度,保障轮盘设计的可靠性及安全性。
附图说明
图1为实施例中基于温度影响的轮盘材料利用系数确定方法流程图;
图2为理想塑性材料与塑性材料的应力—应变关系差异示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
实施例1
参见图1-图2,一种基于温度影响的轮盘材料利用系数确定方法,包括:
分析获取同一材质不同设计尺寸的轮盘在不同温度下应力分布,计算对应温度下各轮盘平均周向应力;
开展或仿真轮盘旋转破裂试验,获得各轮盘不同温度下的实际破裂转速值;
根据不同温度下的平均周向应力、轮盘材料拉伸强度以及对应的实际破裂转速值,计算不同温度下轮盘的初始材料利用系数值;
根据不同温度下轮盘材料拉伸强度、屈服强度性能数据,分析得到不同温度的轮盘材料拉伸强度置信度为50%的-3,以及屈服强度置信度为50%的-3
采用公式拟合得到轮盘的材料利用系数表达式,其中均为常数;
根据拟合得到轮盘的材料利用系数表达式以及轮盘设计工作温度,计算设计工作温度下的轮盘材料利用系。
在本实施例中,同一材质不同设计尺寸的轮盘在不同温度下应力分布可以通过仿真软件仿真获得,也可以采用其他力学分析的方式获得;根据不同温度下轮盘的应力分布计算对应温度下各轮盘平均周向应力,通过仿真或者试验的方式获得各轮盘不同温度下的实际破裂转速值后,以不同温度下的平均周向应力、轮盘材料拉伸强度以及对应的实际破裂转速值进行不同温度下轮盘的初始材料利用系数值计算。采用公式进行曲线拟合,得到的值,得到轮盘的材料利用系数的表达式;根据拟合得到轮盘的材料利用系数表达式以及轮盘设计工作温度,计算设计工作温度下的轮盘材料利用系数。通常轮盘材料表现为较好的塑性性能,轮盘破裂试验过程中子午面整体首先达到屈服强度,然后再达到拉伸强度破坏,该过程中塑性材料会不断强化,往往屈服强度低于拉伸强度。因此本实施例中轮盘材料利用系数的确定方法综合考虑了轮盘材料性能差异的影响,基于轮盘理想塑性材料的基本假设,采取拉伸拉伸强度与屈服强度的比值来反映轮盘的塑性行为,并通过轮盘实测材料性能数据(拉伸拉伸强度和屈服强度)和轮盘破裂试验结果建立了材料利用系数计算的关系式,利用该轮盘材料利用系数进行轮盘设计时,可满足不同温度下同种材料的轮盘破裂转速计算,提高轮盘强度设计精度,保障轮盘设计的可靠性及安全性。
实施例2
本实施例以粉末高温合金轮盘材料利用系数确定流程为例,对本发明的方法做更进一步的解释,主要包括以下步骤:
1)计算对应温度下各轮盘平均周向应力,包括:
建立二维轮盘有限元计算模型,并对轮盘有限元计算模型施加转速、温度和轴向位移约束,分析得到轮盘应力分布,并根据轮盘应力分布求解得到平均周向应力,其中为有限元单元的平均周向应力,为轮盘子午面上第个有限元单元面积,为轮盘子午面上有限元单元总数。本实施例中的轮盘材料利用系数确定方法中,在考虑轮盘材料性能的差异的基础上,通过有二维轮盘限元分析的方式也考虑了不同尺寸的轮盘几何形状的影响,可进一步提高轮盘强度设计精度并保障轮盘设计的可靠性及安全性。
2)开展或仿真轮盘旋转破裂试验,获得各轮盘不同温度下的实际破裂转速值;
3)根据不同温度下的平均周向应力、轮盘材料拉伸强度以及对应的实际破裂转速值,计算不同温度下轮盘的初始材料利用系数值;本实施例中,根据公式计算得到不同温度下轮盘的初始材料利用系数值;其中为材料利用系数,为轮盘实际破裂转速,为轮盘试验基准转速,为轮盘平均周向应力,为轮盘材料在轮盘按面积加权平均温度下的拉伸强度极限值。
4)收集整理随炉入厂复验材料性能数据,获得轮盘材料不同温度下的拉伸强度、屈服强度性能数据;根据不同温度下轮盘材料拉伸强度、屈服强度性能数据,分析得到不同温度的轮盘材料拉伸强度置信度为50%的-3,以及屈服强度置信度为50%的-3;其中,不同温度的轮盘材料拉伸强度置信度为50%的-3是根据公式计算获得;
其中为多个温度下的轮盘材料在轮盘按面积加权平均温度下的拉伸强度极限值的样本平均值,为样本标准差,为单侧容限系数,为存活率为99.87的标准正态偏量,取值为-3;为置信度为50%的t分布百分位值,为子样数。
5)利用最小二乘法拟合,采用公式拟合得到粉末高温合金轮盘材料利用系数表达式;
6)根据拟合得到轮盘的材料利用系数表达式以及轮盘设计工作温度,计算设计工作温度下的轮盘材料利用系数。
以上仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种基于温度影响的轮盘材料利用系数确定方法,其特征在于,包括:
分析获取同一材质不同设计尺寸的轮盘在不同温度下应力分布,计算对应温度下各轮盘平均周向应力;
开展轮盘旋转破裂试验,获得各轮盘不同温度下的实际破裂转速值;
根据不同温度下的平均周向应力、轮盘材料拉伸强度以及对应的实际破裂转速值,计算不同温度下轮盘的初始材料利用系数值;
根据不同温度下轮盘材料拉伸强度、屈服强度性能数据,分析得到不同温度的轮盘材料拉伸强度置信度为50%的-3,以及屈服强度置信度为50%的-3
采用公式拟合得到轮盘的材料利用系数表达式,其中均为常数;
根据拟合得到轮盘的材料利用系数表达式以及轮盘设计工作温度,计算设计工作温度下的轮盘材料利用系数。
2.根据权利要求1所述基于温度影响的轮盘材料利用系数确定方法,其特征在于,计算对应温度下各轮盘平均周向应力,包括:
建立二维轮盘有限元计算模型,并对轮盘有限元计算模型施加转速、温度和轴向位移约束,分析得到轮盘应力分布,并根据轮盘应力分布求解得到平均周向应力,其中为有限元单元的平均周向应力,为轮盘子午面上第个有限元单元面积,为轮盘子午面上有限元单元总数。
3.根据权利要求2所述的基于温度影响的轮盘材料利用系数确定方法,其特征在于,根据公式计算得到不同温度下轮盘的初始材料利用系数值;其中为材料利用系数,为轮盘实际破裂转速,为轮盘试验基准转速,为轮盘平均周向应力,为轮盘材料在轮盘按面积加权平均温度下的拉伸强度极限值。
4.根据权利要求2所述的基于温度影响的轮盘材料利用系数确定方法,其特征在于,不同温度的轮盘材料拉伸强度置信度为50%的-3是根据公式计算获得;
其中为多个温度下的轮盘材料在轮盘按面积加权平均温度下的拉伸强度极限值的样本平均值,为样本标准差,为单侧容限系数,为存活率为99.87的标准正态偏量,取值为-3;为置信度为50%的t分布百分位值,为子样数。
CN202310340757.4A 2023-04-03 2023-04-03 一种基于温度影响的轮盘材料利用系数确定方法 Active CN116052821B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310340757.4A CN116052821B (zh) 2023-04-03 2023-04-03 一种基于温度影响的轮盘材料利用系数确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310340757.4A CN116052821B (zh) 2023-04-03 2023-04-03 一种基于温度影响的轮盘材料利用系数确定方法

Publications (2)

Publication Number Publication Date
CN116052821A true CN116052821A (zh) 2023-05-02
CN116052821B CN116052821B (zh) 2023-06-30

Family

ID=86131700

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310340757.4A Active CN116052821B (zh) 2023-04-03 2023-04-03 一种基于温度影响的轮盘材料利用系数确定方法

Country Status (1)

Country Link
CN (1) CN116052821B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105488252A (zh) * 2015-11-22 2016-04-13 沈阳黎明航空发动机(集团)有限责任公司 一种无螺栓结构涡轮盘挡板的破裂转速仿真分析方法
CN109543213A (zh) * 2018-10-10 2019-03-29 南京航空航天大学 一种非等温轮盘破裂转速的简化预测方法
CN111238804A (zh) * 2018-11-28 2020-06-05 中国航发商用航空发动机有限责任公司 航空发动机转子的超转破裂试验的试验件构型方法
EP3960994A1 (en) * 2020-08-31 2022-03-02 General Electric Company Hybrid electric aircraft engine
CN114293964A (zh) * 2021-12-30 2022-04-08 中国石油大学(北京) 一种室内评价空间应力场主动利用效果的系统及方法
CN114706920A (zh) * 2022-03-31 2022-07-05 中国航发沈阳发动机研究所 一种基于多参数耦合的轮盘破裂转速预测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105488252A (zh) * 2015-11-22 2016-04-13 沈阳黎明航空发动机(集团)有限责任公司 一种无螺栓结构涡轮盘挡板的破裂转速仿真分析方法
CN109543213A (zh) * 2018-10-10 2019-03-29 南京航空航天大学 一种非等温轮盘破裂转速的简化预测方法
CN111238804A (zh) * 2018-11-28 2020-06-05 中国航发商用航空发动机有限责任公司 航空发动机转子的超转破裂试验的试验件构型方法
EP3960994A1 (en) * 2020-08-31 2022-03-02 General Electric Company Hybrid electric aircraft engine
CN114293964A (zh) * 2021-12-30 2022-04-08 中国石油大学(北京) 一种室内评价空间应力场主动利用效果的系统及方法
CN114706920A (zh) * 2022-03-31 2022-07-05 中国航发沈阳发动机研究所 一种基于多参数耦合的轮盘破裂转速预测方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
A. SHINOTAKE 等: "高炉寿命与利用系数之间的关系", 《鞍钢技术》, pages 51 - 54 *
GUO, YUE 等: "Design Optimization and Burst Speed Prediction of a Ti2AlNb Blisk", 《INTERNATIONAL JOURNAL OF AEROSPACE ENGINEERING》, pages 1 - 12 *
MARUTHI B H 等: "Finite Element Formulation for Prediction of Over-speed and burst-margin limits in Aero-engine disc", 《INTERNATIONAL JOURNAL OF SOFT COMPUTING AND ENGINEERING》, vol. 2, no. 3, pages 172 - 176 *
吴长波 等: "某高压涡轮整体叶盘破裂转速计算方法及试验验证", 《燃气涡轮试验与研究》, vol. 19, no. 3, pages 33 - 36 *
李根 等: "FGH96合金高压涡轮转子超转破裂试验验证及分析", 《中国新技术新产品》, pages 61 - 64 *
秦仕勇 等: "粉末冶金涡轮盘破裂转速分析与验证", 《燃气涡轮试验与研究》, vol. 30, no. 5, pages 36 - 41 *

Also Published As

Publication number Publication date
CN116052821B (zh) 2023-06-30

Similar Documents

Publication Publication Date Title
Niu et al. Probabilistic modeling of uncertainties in fatigue reliability analysis of turbine bladed disks
WO2021098169A1 (zh) 一种几何不连续结构的疲劳寿命预测方法
CN105628511B (zh) 一种高温合金蠕变寿命预测方法
CN109583147B (zh) 一种离心叶轮预旋转盘心梯度模拟件设计方法
CN112906281B (zh) 一种基于拟蒙特卡洛抽样的涡轮盘裂纹扩展可靠性分析方法
US11181897B2 (en) Damage probability calculation device and damage probability calculation method
CN114330055A (zh) 基于转子稳态变形假设的盘轴一体转子破裂转速预测方法
CN114492107A (zh) 一种考核涡轮盘轮缘凸块的疲劳试验模拟件设计方法
CN111090957A (zh) 一种高温结构危险点应力-应变计算方法
Fu et al. Study on the fatigue life and damage accumulation of a compressor blade based on a modified nonlinear damage model
Polanka et al. Turbine Tip and Shroud Heat Transfer and Loading: Part B—Comparisons Between Prediction and Experiment Including Unsteady Effects
CN116052821B (zh) 一种基于温度影响的轮盘材料利用系数确定方法
CN110135006B (zh) 镍基单晶涡轮冷却叶片气膜孔的损伤判定方法
Bhachu et al. Application of 3D fracture mechanics for improved crack growth predictions of gas turbine components
SNEDDEN The Application of Non-Axisymmetric Endwall Contouring in a 1½ Stage, Rotating Turbine
CN115062378A (zh) 一种基于晶体塑性的激光冲击强化粉末高温合金低周疲劳寿命预测方法
Vasilyev et al. NUMERICAL METHOD OF SINGLE-CRYSTAL TURBINE BLADE STATIC STRENGTH ESTIMATION TAKING INTO ACCOUNT PLASTICITY AND CREEP EFFECTS.
Arkhipov et al. Modeling of cyclic life for compressor rotor of gas turbine engine taking into account production deviations
CN114706920A (zh) 一种基于多参数耦合的轮盘破裂转速预测方法
EP3161692A1 (en) Method for analysis of 3d features using a 2d probabilistic analysis
CN114139276A (zh) 一种盘轴一体式整体叶盘结构疲劳寿命分析方法
Mazur et al. Evaluation of Creep Damage in a Gas Turbine First Stage Blade
Nozhnitsky et al. Numerical simulation of spin testing for turbo machine disks using energy-based fracture criteria
Vasilyev Stress–Strain State Prediction of High-Temperature Turbine Single Crystal Blades Using Developed Plasticity and Creep Models
Botto et al. A method for on-line temperature calculation of aircraft engine turbine discs

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