CN107991200A - 一种钛合金叶轮的疲劳寿命预测方法 - Google Patents
一种钛合金叶轮的疲劳寿命预测方法 Download PDFInfo
- Publication number
- CN107991200A CN107991200A CN201711189300.9A CN201711189300A CN107991200A CN 107991200 A CN107991200 A CN 107991200A CN 201711189300 A CN201711189300 A CN 201711189300A CN 107991200 A CN107991200 A CN 107991200A
- Authority
- CN
- China
- Prior art keywords
- impeller
- mrow
- titanium alloy
- fatigue life
- stress
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/32—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0005—Repeated or cyclic
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0073—Fatigue
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明提供了一种钛合金叶轮的疲劳寿命预测方法,包括:步骤(1)叶轮仿真中的主要参数选择;步骤(2)流固耦合下的叶轮静态结构强度分析;步骤(3)流固耦合下的叶轮动态结构强度分析;步骤(4)叶轮仿真载荷谱的处理;步骤(5)钛合金叶轮疲劳寿命预测,本发明能够准确预测钛合金叶轮的疲劳寿命。本发明针对钛合金叶轮载荷工况,提出的疲劳寿命预测模型更为精确,能够适应超高寿命的疲劳预测。
Description
技术领域
本发明涉及疲劳寿命预测的技术领域,具体涉及一种钛合金叶轮的疲劳寿命预测方法,适用于钛合金材料的回转类工作零件的疲劳寿命预测。
背景技术
近年来,随着航空技术的发展,各种高性能涡轮轴发动机的应用越来越广泛,而随之而来的发动机叶片共振、震颤、疲劳断裂等问题日渐凸显,占到发动机故障的60%。目前,国内航空发动机传动生产企业引进国外先进生产加工设备,已经可以加工出满足设计图纸的航空发动机叶轮,初步具备航空发动机叶轮的批量生产能力。但是,由于叶轮其结构及形面复杂、加工精度要求高、材料难去除以及机械加工空间有限等原因,而且国内航空发动机传动生产企业缺乏一定的加工理论技术作为指导,国产航空发动机叶轮的抗振抗疲劳性能指标远远落后与国外水平。因此,高性能的发动机叶轮在设计及制造过程中,急切的需要一种准确的疲劳寿命预测方法来指导叶轮的设计和制造。
发明内容
本发明主要解决的技术问题为:现今,航空发动机叶轮设计制造中采用的寿命预测方法仍为有限寿命方法和无限寿命方法,远远不能满足航空发动机的设计制造需求,本发明针对钛合金材料,结合发动机叶轮工况特点和材料特性,提出了一种适用于钛合金材料的叶轮疲劳寿命预测方法,能够为现在的发动机叶轮设计制造提供依据。
本发明的目的是针对现有技术的不足,提供一种钛合金叶轮的疲劳寿命预测方法,其方法流程如下:
步骤(1)、叶轮仿真中的主要参数选择;
步骤(2)、流固耦合下的叶轮静态结构强度分析;
步骤(3)、流固耦合下的叶轮动态结构强度分析;
步骤(4)、叶轮仿真载荷谱的处理;
步骤(5)、钛合金叶轮疲劳寿命预测。
进一步的、所述步骤(1)中需要确定的仿真参数为加叶轮工作温度t(℃)、工作转速v(RPM)、进口流量m(kg/s)、出口压力a(Mpa)以及钛合金的弹性模量E、泊松比υ。
进一步的、所述步骤(2)需要进行流固耦合下的叶轮静态结构强度分析,选取叶轮中应力最大的点作为疲劳载荷监测点。
进一步的、所述步骤(3)需要进行流固耦合下的叶轮动态结构强度分析,提取疲劳载荷检测点的应力变化曲线作为叶轮的疲劳载荷谱。
进一步的、所述步骤(4)需要把叶轮疲劳载荷谱通过雨流法整理出整个历程中所有大小的完整应力循环,并进行分类。
进一步的、所述步骤(5)中各应力循环对应的疲劳寿命计算公式为:
其中Ni表示第i个应力所对应的循环次数,δa表示各应力循环对应的应力幅值,R表示各应力循环对应的应力比,β=0.2579,λ=0.005,μ=E/2(1+υ)。
进一步的、所述步骤(5)中叶轮寿命的计算公式为:
其中N表示钛合金叶轮的疲劳寿命,Ni表示第i个应力所对应的循环次数,n表示雨流法整理出的完整应力循环的个数。
本发明的原理:通过流固耦合方法确定叶轮疲劳载荷,使用建立的钛合金寿命预测模型进行叶轮疲劳寿命预测。
本发明与现有技术相比的有益效果是:首先,提供了一套完整、可靠的寿命预测流程;其次,本发明针对钛合金叶轮载荷工况,提出的疲劳寿命预测模型更为精确,能够适应超高寿命的疲劳预测。本发明可操作性强、精确度高,只需经过仿真分析和数据计算,即可得到接近真实寿命的钛合金叶轮寿命,为钛合金材料的叶轮的设计和制造提供可靠的数据支撑。
附图说明
图1为本发明的一种钛合金叶轮的疲劳寿命预测方法流程图。
具体实施方式
下面结合附图以及具体实施例进一步说明本发明。
本发明提供一种钛合金叶轮的疲劳寿命预测方法,其方法流程如下:
步骤(1)、叶轮仿真中的主要参数选择;
所述步骤(1)中需要确定的仿真参数为加叶轮工作温度t(℃)、工作转速v(RPM)、进口流量m(kg/s)、出口压力a(Mpa)以及钛合金的弹性模量E(Gpa)、泊松比υ。
步骤(2)、流固耦合下的叶轮静态结构强度分析;
所述步骤(2)需要进行流固耦合下的叶轮静态结构强度分析,选取叶轮中应力最大的点作为疲劳载荷监测点。
步骤(3)、流固耦合下的叶轮动态结构强度分析;
所述步骤(3)需要进行流固耦合下的叶轮动态结构强度分析,提取疲劳载荷检测点的应力变化曲线作为叶轮的疲劳载荷谱。
步骤(4)、叶轮仿真载荷谱的处理;
所述步骤(4)需要把叶轮疲劳载荷谱通过雨流法整理出整个历程中所有大小的完整应力循环,并进行分类。
步骤(5)、钛合金叶轮疲劳寿命预测。
所述步骤(5)中各应力循环对应的疲劳寿命计算公式为:
其中Ni表示第i个应力所对应的循环次数,δa表示各应力循环对应的应力幅值,R表示各应力循环对应的应力比,β=0.2579,λ=0.005,μ=E/2(1+υ)。
所述步骤(5)中叶轮寿命的计算公式为:
其中N表示钛合金叶轮的疲劳寿命,Ni表示第i个应力所对应的循环次数,n表示雨流法整理出的完整应力循环的个数。
实施例1
下面以某发动机叶轮为例,具体说明本发明方法,但本发明的保护范围不限于下述实例:
步骤(1),确定的仿真参数为加叶轮工作温度t(℃)、工作转速v(RPM)、进口流量m(kg/s)、出口压力a(Mpa)以及钛合金的弹性模量E、泊松比υ,具体数据如表1所示:
表1仿真主要参数
步骤(2),进行流固耦合下的叶轮静态结构强度分析,找出叶轮应力最大值对应的节点编号为19149。
步骤(3),进行流固耦合下的叶轮动态结构强度分析,提取19149号节点所对用的疲劳载荷谱,具体数据如表2所示:
表2疲劳载荷谱
步骤(4)把叶轮疲劳载荷谱通过雨流法整理出整个历程中所有大小的完整应力循环,并进行分类,具体数据如表3所示:
表3处理后的疲劳载荷
步骤(5)根据疲劳寿命计算公式计算各循环应力所对应的疲劳寿命,数据如下:
表4各载荷对应的疲劳寿命
最终疲劳寿命预测结果:
Claims (7)
1.一种钛合金叶轮的疲劳寿命预测方法,其特征在于,该方法步骤如下:
步骤(1)、叶轮仿真中的主要参数选择;
步骤(2)、流固耦合下的叶轮静态结构强度分析;
步骤(3)、流固耦合下的叶轮动态结构强度分析;
步骤(4)、叶轮仿真载荷谱的处理;
步骤(5)、钛合金叶轮疲劳寿命预测。
2.根据权利要求1所述钛合金叶轮的疲劳寿命预测方法,其特征在于:所述步骤(1)中需要确定的仿真参数为加叶轮工作温度t(℃)、工作转速v(RPM)、进口流量m(kg/s)、出口压力a(Mpa)以及钛合金的弹性模量E、泊松比υ。
3.根据权利要求1所述钛合金叶轮的疲劳寿命预测方法,其特征在于:所述步骤(2)需要进行流固耦合下的叶轮静态结构强度分析,选取叶轮中应力最大的点作为疲劳载荷监测点。
4.根据权利要求1所述钛合金叶轮的疲劳寿命预测方法,其特征在于:所述步骤(3)需要进行流固耦合下的叶轮动态结构强度分析,提取疲劳载荷检测点的应力变化曲线作为叶轮的疲劳载荷谱。
5.根据权利要求1所述钛合金叶轮的疲劳寿命预测方法,其特征在于:所述步骤(4)需要把叶轮疲劳载荷谱通过雨流法整理出整个历程中所有大小的完整应力循环,并进行分类。
6.根据权利要求1所述钛合金叶轮的疲劳寿命预测方法,其特征在于:所述步骤(5)中各应力循环对应的疲劳寿命计算公式为:
<mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>{</mo>
<mn>0.0051</mn>
<mo>*</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>R</mi>
</mrow>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
<msup>
<mo>)</mo>
<mi>&beta;</mi>
</msup>
</mrow>
<msub>
<mi>&delta;</mi>
<mi>a</mi>
</msub>
</mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mn>32</mn>
<msup>
<mi>&mu;</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>}</mo>
</mrow>
<mn>23.3</mn>
</msup>
</mrow>
其中Ni表示第i个应力所对应的循环次数,δa表示各应力循环对应的应力幅值,R表示各应力循环对应的应力比,β=0.2579,λ=0.005,μ=E/2(1+υ)。
7.根据权利要求1所述钛合金叶轮的疲劳寿命预测方法,其特征在于:所述步骤(5)中叶轮寿命的计算公式为:
<mrow>
<mi>N</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<mi>n</mi>
</munderover>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</mfrac>
</mrow>
</mfrac>
</mrow>
其中N表示钛合金叶轮的疲劳寿命,Ni表示第i个应力所对应的循环次数,n表示雨流法整理出的完整应力循环的个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711189300.9A CN107991200B (zh) | 2017-11-24 | 2017-11-24 | 一种钛合金叶轮的疲劳寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711189300.9A CN107991200B (zh) | 2017-11-24 | 2017-11-24 | 一种钛合金叶轮的疲劳寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107991200A true CN107991200A (zh) | 2018-05-04 |
CN107991200B CN107991200B (zh) | 2020-02-14 |
Family
ID=62032783
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711189300.9A Expired - Fee Related CN107991200B (zh) | 2017-11-24 | 2017-11-24 | 一种钛合金叶轮的疲劳寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107991200B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110717229A (zh) * | 2019-10-16 | 2020-01-21 | 北京航空航天大学 | 一种基于试验数据的虚拟载荷谱自更新构造方法 |
CN112069616A (zh) * | 2020-08-20 | 2020-12-11 | 西安空天能源动力智能制造研究院有限公司 | 一种退役航空发动机再利用的智能寿命延长控制方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102436515A (zh) * | 2010-08-27 | 2012-05-02 | 通用电气公司 | 用于估算涡轮机翼型的剩余寿命的方法和系统 |
US20120288373A1 (en) * | 2011-05-13 | 2012-11-15 | Hamilton Sundstrand Corporation | Rotor with asymmetric blade spacing |
CN104236915A (zh) * | 2014-09-19 | 2014-12-24 | 华北电力大学 | 发电机组轴系-叶片扭振疲劳寿命损耗评估方法及系统 |
CN104268335A (zh) * | 2014-09-23 | 2015-01-07 | 工业和信息化部电子第五研究所 | 微组装组件振动疲劳寿命预测方法和系统 |
KR20160067402A (ko) * | 2014-12-04 | 2016-06-14 | 한국생산기술연구원 | 부등피치 재생 블로워 및 이의 최적화 설계 방법 |
CN106644464A (zh) * | 2016-11-18 | 2017-05-10 | 南京工业大学 | 一种基于载荷谱分析的轧机传动系统关键零部件的疲劳寿命预警方法 |
CN106844958A (zh) * | 2017-01-19 | 2017-06-13 | 沈阳航空航天大学 | 基于改进雨流计数法的薄壁结构热声疲劳寿命预估方法 |
CN107228761A (zh) * | 2017-06-02 | 2017-10-03 | 中国科学院金属研究所 | 一种复杂载荷和环境下的构件寿命预测方法 |
-
2017
- 2017-11-24 CN CN201711189300.9A patent/CN107991200B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102436515A (zh) * | 2010-08-27 | 2012-05-02 | 通用电气公司 | 用于估算涡轮机翼型的剩余寿命的方法和系统 |
US20120288373A1 (en) * | 2011-05-13 | 2012-11-15 | Hamilton Sundstrand Corporation | Rotor with asymmetric blade spacing |
CN104236915A (zh) * | 2014-09-19 | 2014-12-24 | 华北电力大学 | 发电机组轴系-叶片扭振疲劳寿命损耗评估方法及系统 |
CN104268335A (zh) * | 2014-09-23 | 2015-01-07 | 工业和信息化部电子第五研究所 | 微组装组件振动疲劳寿命预测方法和系统 |
KR20160067402A (ko) * | 2014-12-04 | 2016-06-14 | 한국생산기술연구원 | 부등피치 재생 블로워 및 이의 최적화 설계 방법 |
CN106644464A (zh) * | 2016-11-18 | 2017-05-10 | 南京工业大学 | 一种基于载荷谱分析的轧机传动系统关键零部件的疲劳寿命预警方法 |
CN106844958A (zh) * | 2017-01-19 | 2017-06-13 | 沈阳航空航天大学 | 基于改进雨流计数法的薄壁结构热声疲劳寿命预估方法 |
CN107228761A (zh) * | 2017-06-02 | 2017-10-03 | 中国科学院金属研究所 | 一种复杂载荷和环境下的构件寿命预测方法 |
Non-Patent Citations (3)
Title |
---|
SHU LIU 等: "Fatigue life assessment of centrifugal compressor impeller based on FEA", 《ENGINEERING FAILURE ANALYSIS》 * |
丛建华 等: ""基于ABAQUS的离心叶轮力学性能研究"", 《工程与试验》 * |
段远飞: "某航空发动机离心叶轮疲劳使用寿命与裂纹扩展分析", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110717229A (zh) * | 2019-10-16 | 2020-01-21 | 北京航空航天大学 | 一种基于试验数据的虚拟载荷谱自更新构造方法 |
CN112069616A (zh) * | 2020-08-20 | 2020-12-11 | 西安空天能源动力智能制造研究院有限公司 | 一种退役航空发动机再利用的智能寿命延长控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107991200B (zh) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Torresi et al. | Accurate numerical simulation of a high solidity Wells turbine | |
CN104481601B (zh) | 汽轮机转子低周疲劳与高周疲劳的裂纹扩展寿命监控方法 | |
CN111382544B (zh) | 一种基于等效疲劳寿命的压气机叶片加速寿命试验方法 | |
CN103268388B (zh) | 一种航空发动机涡轮叶盘系统可靠性分析方法 | |
Krishnakanth et al. | Structural and thermal analysis of gas turbine blade by using FEM | |
CN109033488B (zh) | 一种基于水机电耦合模型的抽水蓄能机组寿命分析方法 | |
CN107991200A (zh) | 一种钛合金叶轮的疲劳寿命预测方法 | |
CN104929944B (zh) | 基于间隙流的多级泵转子动力密封的设计方法 | |
Mischo et al. | Numerical and experimental FSI-study to determine mechanical stresses induced by rotating stall in unshrouded centrifugal compressor impellers | |
Sladojevic et al. | Investigation of the influence of aerodynamic coupling on response levels of mistuned bladed discs with weak structural coupling | |
CN105697420B (zh) | 部分处理机匣性能预估模型 | |
Anding et al. | Development of a Novel Axial Compressor Generation for Industrial Applications: Part 2—Blade Mechanics | |
Huang et al. | Fully clocking effect in a two-stage compressor | |
Umamaheswararao et al. | Design and Analysis of turbine blade by using fem | |
Wang et al. | Mistuning identification for integrally bladed disks based on the SNM technique | |
Agnalt | Rotor-Stator Interaction in Low-Specific Speed Francis Turbines | |
Dragan et al. | Global analysis of tip clearance influence on centrifugal compressor performance | |
Xie et al. | Strength design and numerical analysis of radial inflow turbine impeller for a 100kw microturbine | |
Patel et al. | Fe analysis of runner blade for wells turbine | |
Binama et al. | Runner Blade Number Influencing the RPT Runner Upstream Flow Characteristics: A CFD Numerical Simulation | |
Repetckii et al. | Mathematical modeling and the computer analysis of strength characteristics and sensitivity of high-loaded parts of power turbines | |
Sheard et al. | The aerodynamic and mechanical performance of a high-pressure turbine stage in a transient wind tunnel | |
CN114151146B (zh) | 多联带冠涡轮转子叶片气流激振力参数的获取方法 | |
Wang et al. | Inlet compatibility and fan aeromechanics of HBP turbofan engine | |
Zhao et al. | Influence of blade vibration on part-span rotating stall |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200214 Termination date: 20201124 |
|
CF01 | Termination of patent right due to non-payment of annual fee |