CN109520717B - 一种发动机主轴疲劳试验载荷确定方法 - Google Patents
一种发动机主轴疲劳试验载荷确定方法 Download PDFInfo
- Publication number
- CN109520717B CN109520717B CN201811559464.0A CN201811559464A CN109520717B CN 109520717 B CN109520717 B CN 109520717B CN 201811559464 A CN201811559464 A CN 201811559464A CN 109520717 B CN109520717 B CN 109520717B
- Authority
- CN
- China
- Prior art keywords
- load
- cyclic
- equivalent
- cycle
- fatigue
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供的一种针对发动机轴类零件疲劳试验载荷的确定方法。通过将旋转弯矩载荷等效转换为振动扭矩载荷,解决轴类疲劳试验中旋转弯矩载荷不能在普通疲劳试验机上加载的问题,使得可以使用普通疲劳试验机开展小尺寸的轴类零件疲劳试验,不需要再针对被试轴尺寸进行专门的轴类疲劳试验器设计、重新建设或大量改造,从而降低轴类疲劳试验的成本,也降低了试验载荷加载难度,提高了载荷加载精度。
Description
技术领域
本发明属于发动机领域,涉及一种发动机主轴疲劳试验载荷确定方法。
背景技术
航空发动机轴类零件,是转子上起支持作用的零件,也是转子上传递扭矩(功率)的零件。发动机主轴是非常重要的承力件,在工作中承受包括主扭矩、轴向力和轴对称载荷的低循环疲劳载荷和包括陀螺力矩、振动扭矩、不平衡力等的高循环疲劳载荷同时作用。在发动机设计中,需要通过轴类疲劳试验对其寿命进行考核。
工程上,主轴疲劳试验载荷的确定一般直接由标准循环载荷乘以散度系数作为试验载荷谱,即直接由轴向力、扭矩、振动扭矩和旋转弯矩等的标准循环载荷分别乘以对应的散度系数,得到轴向力、扭矩、振动扭矩和旋转弯矩等的标准循环载荷试验载荷谱。
由于普通的疲劳试验机一般可同时加载轴向力、扭矩、振动扭矩载荷,不能实现旋转弯矩的加载,通过这种方法确定的试验载荷仅能在专用轴类疲劳试验器上进行试验加载。因此往往需要根据被试轴尺寸进行专门的轴类疲劳试验器设计、建设或大量改造,耗时长,成本高。且试验载荷种类多,加载要求复杂,加载精度和加载频率不高。
发明内容
本发明的目的是:提供一种对发动机轴类零件疲劳试验载荷的确定方法。通过对轴类零件的疲劳试验载荷进行等效转换,使得试验载荷能够模拟零件疲劳载荷,且可以使用常用的材料疲劳试验机加载,解决主轴疲劳试验仅能在专门轴类疲劳试验器上进行的问题。本方法确定的载荷也适用于使用轴类疲劳试验器进行的主轴疲劳试验。
为解决上述技术问题,本发明的技术方案是:
一种发动机主轴疲劳试验载荷确定方法,所述方法包括:确定主轴的标准循环载荷;通过等效转换,将所述标准循环载荷转换为等效循环载荷;将所述等效循环载荷乘以对应的散度系数,得到试验载荷谱。
可选的,所述确定主轴的标准循环载荷,具体包括:根据扭矩MT、振动扭矩ΔMT、轴向力P、旋转弯矩Mb确定主轴的标准循环载荷。
可选的,所述通过等效转换,将所述标准循环载荷转换为等效循环载荷,具体包括:
按照载荷作用的循环次数要求,将扭矩MT、轴向力P划分为低循环载荷,将振动扭矩ΔMT、旋转弯矩Mb划分为高循环载荷;
分别计算主轴在所述低循环载荷和所述高循环载荷作用下的应力,得到主轴在所述低循环载荷峰值作用下的应力幅值σm和在所述高循环载荷作用下的应力幅值σa;
使用Goodman曲线评估主轴在标准循环载荷下的高循环疲劳寿命储备α;
选取适当的振动扭矩值ΔMT′,计算主轴在载荷ΔMT′作用下寿命考核部位的应力,得到应力幅值σ'a;
通过调整振动扭矩ΔMT′值,使得寿命考核部位的高循环疲劳储备α′与α相当,则可确定出等效后的高循环载荷ΔMT′。
可选的,所述将所述等效循环载荷乘以对应的散度系数,得到试验载荷谱,具体包括:
根据标准循环的低循环载荷和等效后的高循环载荷振动扭矩ΔMT′分别乘以疲劳散度系数得到试验的低循环载荷MTs、Ps和试验高循环载荷ΔMTs:
MTs=K1×K2×K3×MT…………………………………………(4)
Ps=K1×K2×K3×P……………………………………………(5)
ΔMTs=K1'×K2×K3×ΔM'T………………………………………(6)
式中:MTs为试验载荷的扭矩;Ps为试验载荷的轴向力;ΔMTs为试验载荷的振动扭矩;K1为低循环疲劳载荷的散度系数;K1'为高循环疲劳载荷的散度系数;K2为材料强度修正系数;K3为温度修正系数;MT为标准循环载荷的主扭矩;P为标准循环载荷的轴向力;ΔMT′为等效后的振动扭矩。
可选的,所述使用Goodman曲线计算主轴在标准循环载荷下的高循环疲劳寿命储备α,具体包括:
根据公式(1)、公式(2)计算主轴在标准循环载荷下的高循环疲劳寿命储备α,其中,
本发明提供的一种针对发动机轴类零件疲劳试验载荷的确定方法。通过将旋转弯矩载荷等效转换为振动扭矩载荷,解决轴类疲劳试验中旋转弯矩载荷不能在普通疲劳试验机上加载的问题,使得可以使用普通疲劳试验机开展小尺寸的轴类零件疲劳试验,不需要再针对被试轴尺寸进行专门的轴类疲劳试验器设计、重新建设或大量改造,从而降低轴类疲劳试验的成本,也降低了试验载荷加载难度,提高了载荷加载精度。
附图说明
图1本发明的方法流程图;
图2本发明的标准循环载荷谱示意图;
图3本发明的Goodman曲线示意图。
具体实施方式
本发明创造是在原寿命散度系数法确定试验载荷的基础上优化改进,在确定主轴的标准循环载荷后,通过加入对载荷的等效转换,将标准循环载荷转换为等效循环载荷,再乘以对应的散度系数,得到试验载荷谱。以轴流式发动机为例,对其技术方案和具体实施方式进行介绍。
发动机主轴疲劳试验载荷确定的具体流程如图1所示。
(1)确定标准循环载荷
对主轴疲劳寿命有影响的载荷主要为扭矩MT、振动扭矩ΔMT、轴向力P、旋转弯矩Mb,主轴寿命考核部位的这些载荷组成了主轴的标准循环载荷(图2)。
(2)区分高低循环载荷
按照载荷作用的循环次数要求,将步骤(1)中确定的标准循环载荷划分为低循环载荷和高循环载荷。其中,扭矩MT、轴向力P为低循环载荷,振动扭矩ΔMT、旋转弯矩Mb为高循环载荷。
(3)高、低循环载荷下的应力分析
分别计算主轴在步骤(2)中划分得到的低循环载荷和高循环载荷作用下的应力,得到主轴在低循环载荷(扭矩MT、轴向力P)峰值作用下的应力幅值σm和在高循环载荷(振动扭矩ΔMT、旋转弯矩Mb)作用下的应力幅值σa。
为获得较高的计算精度,建议使用有限元分析软件对主轴进行二维或三维有限元计算。计算时需要考虑主轴在发动机工作状态下的温度场。
(4)标准循环载荷下的高循环疲劳寿命储备计算
使用Goodman曲线(图3)评估主轴在标准循环载荷(MT、P、ΔMT、Mb)下的高循环疲劳寿命储备α,计算公式见式(1)、(2)。
式中,σa,max——允许的最大高循环应力;
σm——平均应力,取为低循环峰值载荷下的应力幅值;
σ-1——对应温度下的材料疲劳强度;
σb——对应温度下的材料拉伸强度。
式中,σa——高循环载荷下的应力幅值。
(5)高循环载荷等效转换
将标准循环载荷谱中的高循环载荷进行等效转换,即:使用单一的振动扭矩ΔMT′载荷等效代替标准循环载荷谱中振动扭矩ΔMT和旋转弯矩Mb振动扭矩载荷。具体方法如下:选取适当的振动扭矩值ΔMT′,计算主轴在载荷ΔMT′作用下寿命考核部位的应力,得到应力幅值σ'a,计算方法同步骤(3)。使用σ'a和步骤(4)中的σa,max按公式(3)计算出等效循环载荷(MT、P、ΔMT′)下的高循环疲劳寿命储备α'。
通过调整振动扭矩ΔMT′值,使得寿命考核部位的高循环疲劳储备α′与α相当,则可确定出等效后的高循环载荷ΔMT′。
(6)确定试验载荷
试验载荷由试验低循环载荷MTs、Ps和试验高循环载荷ΔMTs组成。由标准循环的低循环载荷(扭矩MT、轴向力和P)和等效后的高循环载荷振动扭矩ΔMT′分别乘以疲劳散度系数得到试验的低循环载荷和高循环载荷试验高循环载荷:
MTs=K1×K2×K3×MT…………………………………………(4)
Ps=K1×K2×K3×P……………………………………………(5)
ΔMTs=K1'×K2×K3×ΔM'T………………………………………(6)
式中:MTs——试验载荷的扭矩;
Ps——试验载荷的轴向力;
ΔMTs——试验载荷的振动扭矩;
K1——低循环疲劳载荷的散度系数;
K1'——高循环疲劳载荷的散度系数;
K2——材料强度修正系数;
K3——温度修正系数;
MT——标准循环载荷的主扭矩;
P——标准循环载荷的轴向力;
ΔMT′——等效后的振动扭矩;
其中,修正系数K1、K1'、K2和K3的选取与经典的疲劳散度系数法取值方法一致。
以上主轴疲劳试验载荷确定方法,以涡喷/涡扇发动机的主轴为例,在用于涡轴发动机主轴试验载荷确定时,还需考虑1×P力F1×P和1×P力矩M1×P等载荷。
Claims (3)
1.一种发动机主轴疲劳试验载荷确定方法,其特征在于:所述方法包括:
确定主轴的标准循环载荷;
通过等效转换,将所述标准循环载荷转换为等效循环载荷;
将所述等效循环载荷乘以对应的散度系数,得到试验载荷谱;
所述确定主轴的标准循环载荷,具体包括:根据扭矩MT、振动扭矩ΔMT、轴向力P、旋转弯矩Mb确定主轴的标准循环载荷;
所述通过等效转换,将所述标准循环载荷转换为等效循环载荷,具体包括:
按照载荷作用的循环次数要求,将扭矩MT、轴向力P划分为低循环载荷,将振动扭矩ΔMT、旋转弯矩Mb划分为高循环载荷;
分别计算主轴在所述低循环载荷和所述高循环载荷作用下的应力,得到主轴在所述低循环载荷峰值作用下的应力幅值σm和在所述高循环载荷作用下的应力幅值σa;
使用Goodman曲线评估主轴在标准循环载荷下的高循环疲劳寿命储备α;
选取等效振动扭矩ΔM′T,计算主轴在等效振动扭矩ΔM′T作用下寿命考核部位的应力,得到应力幅值σ′a;
通过调整等效振动扭矩ΔM′T值,使得寿命考核部位的高循环疲劳储备α′与α相当,则可确定出等效后的等效振动扭矩ΔM′T。
2.根据权利要求1所述的方法,其特征在于:所述将所述等效循环载荷乘以对应的散度系数,得到试验载荷谱,具体包括:
根据标准循环的低循环载荷和等效后的高循环载荷振动扭矩ΔM′T分别乘以疲劳散度系数得到试验的低循环载荷MTs、Ps和试验高循环载荷ΔMTs,具体为:
MTs=K1×K2×K3×MT
Ps=K1×K2×K3×P
ΔMTs=K′1×K2×K3×ΔM′T
式中:MTs为试验载荷的扭矩;Ps为试验载荷的轴向力;ΔMTs为试验载荷的振动扭矩;K1为低循环疲劳载荷的散度系数;K′1为高循环疲劳载荷的散度系数;K2为材料强度修正系数;K3为温度修正系数;MT为扭矩;P为标准循环载荷的轴向力;ΔM′T为等效振动扭矩。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811559464.0A CN109520717B (zh) | 2018-12-20 | 2018-12-20 | 一种发动机主轴疲劳试验载荷确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811559464.0A CN109520717B (zh) | 2018-12-20 | 2018-12-20 | 一种发动机主轴疲劳试验载荷确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109520717A CN109520717A (zh) | 2019-03-26 |
CN109520717B true CN109520717B (zh) | 2020-10-20 |
Family
ID=65795794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811559464.0A Active CN109520717B (zh) | 2018-12-20 | 2018-12-20 | 一种发动机主轴疲劳试验载荷确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109520717B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112199848B (zh) * | 2020-10-16 | 2022-11-22 | 中国航发四川燃气涡轮研究院 | 一种低压涡轮轴疲劳寿命评估方法 |
CN117725802B (zh) * | 2024-02-07 | 2024-04-16 | 中国航发四川燃气涡轮研究院 | 航空发动机主轴疲劳试验标准循环载荷谱构建方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706757A (zh) * | 2012-05-08 | 2012-10-03 | 上海博览达信息科技有限公司 | 多轴疲劳分析方法及其应用 |
CN103942418A (zh) * | 2014-04-03 | 2014-07-23 | 北京工业大学 | 一种多轴载荷条件下细节疲劳强度额定值的确定方法 |
CN105247339A (zh) * | 2013-05-17 | 2016-01-13 | 斯奈克玛 | 循环疲劳或循环的优化及高循环疲劳试验台 |
KR20160027529A (ko) * | 2014-09-01 | 2016-03-10 | 주식회사 대한항공 | 복합재료 날개 구조물의 피로수명 평가 방법 |
CN106777880A (zh) * | 2016-11-16 | 2017-05-31 | 北京航空航天大学 | 一种考虑载荷与材料分散性的复合材料结构可靠性优化设计方法 |
CN106872299A (zh) * | 2017-02-06 | 2017-06-20 | 太原理工大学 | 一种预测镁合金构件疲劳极限的方法 |
CN107092728A (zh) * | 2017-03-30 | 2017-08-25 | 成都航空职业技术学院 | 一种涡扇发动机拉紧轴疲劳试验方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108918106B (zh) * | 2018-05-08 | 2019-12-20 | 中国科学院工程热物理研究所 | 一种分阶段加载的风力机叶片的疲劳测试方法 |
-
2018
- 2018-12-20 CN CN201811559464.0A patent/CN109520717B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706757A (zh) * | 2012-05-08 | 2012-10-03 | 上海博览达信息科技有限公司 | 多轴疲劳分析方法及其应用 |
CN105247339A (zh) * | 2013-05-17 | 2016-01-13 | 斯奈克玛 | 循环疲劳或循环的优化及高循环疲劳试验台 |
CN103942418A (zh) * | 2014-04-03 | 2014-07-23 | 北京工业大学 | 一种多轴载荷条件下细节疲劳强度额定值的确定方法 |
KR20160027529A (ko) * | 2014-09-01 | 2016-03-10 | 주식회사 대한항공 | 복합재료 날개 구조물의 피로수명 평가 방법 |
CN106777880A (zh) * | 2016-11-16 | 2017-05-31 | 北京航空航天大学 | 一种考虑载荷与材料分散性的复合材料结构可靠性优化设计方法 |
CN106872299A (zh) * | 2017-02-06 | 2017-06-20 | 太原理工大学 | 一种预测镁合金构件疲劳极限的方法 |
CN107092728A (zh) * | 2017-03-30 | 2017-08-25 | 成都航空职业技术学院 | 一种涡扇发动机拉紧轴疲劳试验方法 |
Non-Patent Citations (1)
Title |
---|
航空发动机低压涡轮轴疲劳寿命计算分析;崔平亮;《万方数据知识服务平台》;20140917;第55-61页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109520717A (zh) | 2019-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111460583B (zh) | 一种针对复杂几何结构件的蠕变-疲劳寿命设计方法 | |
CN109520717B (zh) | 一种发动机主轴疲劳试验载荷确定方法 | |
CN205785814U (zh) | 航空发动机主轴承的可靠性试验装置 | |
CN105677973B (zh) | 实现单轴与三轴随机振动应力等效的试验谱剪裁方法 | |
CN115640666B (zh) | 一种基于损伤等效的航空发动机加速任务试车谱编制方法 | |
CN105466672A (zh) | 风轮叶片疲劳试验方法及其在全尺寸疲劳试验中的应用 | |
CN113266474B (zh) | 一种加载条件下的航空发动机起动阻力矩测量方法 | |
CN107092728A (zh) | 一种涡扇发动机拉紧轴疲劳试验方法 | |
CN109583147A (zh) | 一种离心叶轮预旋转盘心梯度模拟件设计方法 | |
CN107545110A (zh) | 一种动态应力加速寿命试验剖面编制方法 | |
JP2017053698A (ja) | ランダムな波形の荷重を負荷するシステムおよび方法、疲労試験用プログラム、および疲労試験機 | |
CN117725802B (zh) | 航空发动机主轴疲劳试验标准循环载荷谱构建方法及系统 | |
CN113221261B (zh) | 一种航空传动系统振动限制值的制定方法 | |
CN113432785B (zh) | 一种不加载条件下的航空发动机转动惯量测量方法 | |
CN113266473A (zh) | 一种不加载条件下的航空发动机起动阻力矩测量方法 | |
Williams et al. | A methodology for predicting torsional fatigue crack initiation in large turbine-generator shafts | |
CN101964565B (zh) | 跨外转子混合动平衡方法 | |
CN110646205A (zh) | 一种长寿命航空发动机的加速等效试验方法 | |
CN113268829B (zh) | 一种机械部件正弦扫频振动疲劳损伤估算方法 | |
CN113447201B (zh) | 一种加载条件下的航空发动机转动惯量测量方法 | |
CN111473859B (zh) | 整机振动限制值的制定方法 | |
CN105868451A (zh) | 一种涡轮增压器压气机叶轮叶片振动节线精确测定方法 | |
Berger et al. | Probabilistic vibration and lifetime analysis of regenerated turbomachinery blades | |
CN115982958A (zh) | 一种基于工程损伤力学的材料蠕变疲劳寿命预测方法 | |
DeMarco et al. | An Efficient Method for the Optimization of Viscoplastic Constitutive Model Constants |
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 |