CN109520717A - 一种发动机主轴疲劳试验载荷确定方法 - Google Patents
一种发动机主轴疲劳试验载荷确定方法 Download PDFInfo
- Publication number
- CN109520717A CN109520717A CN201811559464.0A CN201811559464A CN109520717A CN 109520717 A CN109520717 A CN 109520717A CN 201811559464 A CN201811559464 A CN 201811559464A CN 109520717 A CN109520717 A CN 109520717A
- Authority
- CN
- China
- Prior art keywords
- load
- fatigue
- torque
- main shaft
- cycle
- 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
- 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
Abstract
本发明提供的一种针对发动机轴类零件疲劳试验载荷的确定方法。通过将旋转弯矩载荷等效转换为振动扭矩载荷,解决轴类疲劳试验中旋转弯矩载荷不能在普通疲劳试验机上加载的问题,使得可以使用普通疲劳试验机开展小尺寸的轴类零件疲劳试验,不需要再针对被试轴尺寸进行专门的轴类疲劳试验器设计、重新建设或大量改造,从而降低轴类疲劳试验的成本,也降低了试验载荷加载难度,提高了载荷加载精度。
Description
技术领域
本发明属于发动机领域,涉及一种发动机主轴疲劳试验载荷确定方法。
背景技术
航空发动机轴类零件,是转子上起支持作用的零件,也是转子上传递扭矩(功率)的零件。发动机主轴是非常重要的承力件,在工作中承受包括主扭矩、轴向力和轴对称载荷的低循环疲劳载荷和包括陀螺力矩、振动扭矩、不平衡力等的高循环疲劳载荷同时作用。在发动机设计中,需要通过轴类疲劳试验对其寿命进行考核。
工程上,主轴疲劳试验载荷的确定一般直接由标准循环载荷乘以散度系数作为试验载荷谱,即直接由轴向力、扭矩、振动扭矩和旋转弯矩等的标准循环载荷分别乘以对应的散度系数,得到轴向力、扭矩、振动扭矩和旋转弯矩等的标准循环载荷试验载荷谱。
由于普通的疲劳试验机一般可同时加载轴向力、扭矩、振动扭矩载荷,不能实现旋转弯矩的加载,通过这种方法确定的试验载荷仅能在专用轴类疲劳试验器上进行试验加载。因此往往需要根据被试轴尺寸进行专门的轴类疲劳试验器设计、建设或大量改造,耗时长,成本高。且试验载荷种类多,加载要求复杂,加载精度和加载频率不高。
发明内容
本发明的目的是:提供一种对发动机轴类零件疲劳试验载荷的确定方法。通过对轴类零件的疲劳试验载荷进行等效转换,使得试验载荷能够模拟零件疲劳载荷,且可以使用常用的材料疲劳试验机加载,解决主轴疲劳试验仅能在专门轴类疲劳试验器上进行的问题。本方法确定的载荷也适用于使用轴类疲劳试验器进行的主轴疲劳试验。
为解决上述技术问题,本发明的技术方案是:
一种发动机主轴疲劳试验载荷确定方法,所述方法包括:确定主轴的标准循环载荷;通过等效转换,将所述标准循环载荷转换为等效循环载荷;将所述等效循环载荷乘以对应的散度系数,得到试验载荷谱。
可选的,所述确定主轴的标准循环载荷,具体包括:根据扭矩MT、振动扭矩ΔMT、轴向力P、旋转弯矩Mb确定主轴的标准循环载荷。
可选的,所述通过等效转换,将所述标准循环载荷转换为等效循环载荷,具体包括:
按照载荷作用的循环次数要求,将扭矩MT、轴向力P划分为低循环载荷,将振动扭矩ΔMT、旋转弯矩Mb划分为高循环载荷;
分别计算主轴在所述低循环载荷和所述高循环载荷作用下的应力,得到主轴在所述低循环载荷峰值作用下的应力幅值σm和在所述高循环载荷作用下的应力幅值σa;
使用Goodman曲线评估主轴在标准循环载荷下的高循环疲劳寿命储备α;
选取适当的振动扭矩值ΔMT′,计算主轴在载荷ΔMT′作用下寿命考核部位的应力,得到应力幅值σ'a;
使用所述σ'a和σa,max按公式(3)计算出等效循环载荷(MT、P、ΔMT′)下的高循环疲劳寿命储备α';
通过调整振动扭矩Δ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)式中,σa,max为允许的最大高循环应力;σm为平均应力,取为低循环峰值载荷下的应力幅值;σ-1为对应温度下的材料疲劳强度;σb为对应温度下的材料拉伸强度;
公式(2)式中,σa为高循环载荷下的应力幅值。
本发明提供的一种针对发动机轴类零件疲劳试验载荷的确定方法。通过将旋转弯矩载荷等效转换为振动扭矩载荷,解决轴类疲劳试验中旋转弯矩载荷不能在普通疲劳试验机上加载的问题,使得可以使用普通疲劳试验机开展小尺寸的轴类零件疲劳试验,不需要再针对被试轴尺寸进行专门的轴类疲劳试验器设计、重新建设或大量改造,从而降低轴类疲劳试验的成本,也降低了试验载荷加载难度,提高了载荷加载精度。
附图说明
图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 (5)
1.一种发动机主轴疲劳试验载荷确定方法,其特征在于:所述方法包括:
确定主轴的标准循环载荷;
通过等效转换,将所述标准循环载荷转换为等效循环载荷;
将所述等效循环载荷乘以对应的散度系数,得到试验载荷谱。
2.根据权利要求1所述的方法,其特征在于:所述确定主轴的标准循环载荷,具体包括:
根据扭矩MT、振动扭矩ΔMT、轴向力P、旋转弯矩Mb确定主轴的标准循环载荷。
3.根据权利要求2所述的方法,其特征在于:所述通过等效转换,将所述标准循环载荷转换为等效循环载荷,具体包括:
按照载荷作用的循环次数要求,将扭矩MT、轴向力P划分为低循环载荷,将振动扭矩ΔMT、旋转弯矩Mb划分为高循环载荷;
分别计算主轴在所述低循环载荷和所述高循环载荷作用下的应力,得到主轴在所述低循环载荷峰值作用下的应力幅值σm和在所述高循环载荷作用下的应力幅值σa;
使用Goodman曲线评估主轴在标准循环载荷下的高循环疲劳寿命储备α;
选取适当的振动扭矩值ΔM′T,计算主轴在载荷ΔM′T作用下寿命考核部位的应力,得到应力幅值σ'a;
使用所述σ'a和σa,max按公式(3)计算出等效循环载荷(MT、P、ΔM′T)下的高循环疲劳寿命储备α';
通过调整振动扭矩ΔM′T值,使得寿命考核部位的高循环疲劳储备α′与α相当,则可确定出等效后的高循环载荷ΔM′T。
4.根据权利要求3所述的方法,其特征在于:所述将所述等效循环载荷乘以对应的散度系数,得到试验载荷谱,具体包括:
根据标准循环的低循环载荷和等效后的高循环载荷振动扭矩ΔM′T分别乘以疲劳散度系数得到试验的低循环载荷MTs、Ps和试验高循环载荷ΔMTs,具体为:
MTs=K1×K2×K3×MT…………………………………………(4)
Ps=K1×K2×K3×P……………………………………………(5)
ΔMTs=K′1×K2×K3×ΔM'T………………………………………(6)
式中:MTs为试验载荷的扭矩;Ps为试验载荷的轴向力;ΔMTs为试验载荷的振动扭矩;K1为低循环疲劳载荷的散度系数;K′1为高循环疲劳载荷的散度系数;K2为材料强度修正系数;K3为温度修正系数;MT为标准循环载荷的主扭矩;P为标准循环载荷的轴向力;ΔM′T为等效后的振动扭矩。
5.根据权利要求3所述的方法,其特征在于:所述使用Goodman曲线计算主轴在标准循环载荷下的高循环疲劳寿命储备α,具体包括:
根据公式(1)、公式(2)计算主轴在标准循环载荷下的高循环疲劳寿命储备α,其中,
公式(1)式中,σa,max为允许的最大高循环应力;σm为平均应力,取为低循环峰值载荷下的应力幅值;σ-1为对应温度下的材料疲劳强度;σb为对应温度下的材料拉伸强度;
公式(2)式中,σa为高循环载荷下的应力幅值。
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 true CN109520717A (zh) | 2019-03-26 |
CN109520717B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112199848A (zh) * | 2020-10-16 | 2021-01-08 | 中国航发四川燃气涡轮研究院 | 一种低压涡轮轴疲劳寿命评估方法 |
CN117725802A (zh) * | 2024-02-07 | 2024-03-19 | 中国航发四川燃气涡轮研究院 | 航空发动机主轴疲劳试验标准循环载荷谱构建方法及系统 |
Citations (8)
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 | 成都航空职业技术学院 | 一种涡扇发动机拉紧轴疲劳试验方法 |
CN108918106A (zh) * | 2018-05-08 | 2018-11-30 | 中国科学院工程热物理研究所 | 一种分阶段加载的风力机叶片的疲劳测试方法 |
-
2018
- 2018-12-20 CN CN201811559464.0A patent/CN109520717B/zh active Active
Patent Citations (8)
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 | 成都航空职业技术学院 | 一种涡扇发动机拉紧轴疲劳试验方法 |
CN108918106A (zh) * | 2018-05-08 | 2018-11-30 | 中国科学院工程热物理研究所 | 一种分阶段加载的风力机叶片的疲劳测试方法 |
Non-Patent Citations (2)
Title |
---|
崔平亮: "航空发动机低压涡轮轴疲劳寿命计算分析", 《万方数据知识服务平台》 * |
邱宣怀等: "《机械设计》", 31 July 1997, 高等教育出版社 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112199848A (zh) * | 2020-10-16 | 2021-01-08 | 中国航发四川燃气涡轮研究院 | 一种低压涡轮轴疲劳寿命评估方法 |
CN112199848B (zh) * | 2020-10-16 | 2022-11-22 | 中国航发四川燃气涡轮研究院 | 一种低压涡轮轴疲劳寿命评估方法 |
CN117725802A (zh) * | 2024-02-07 | 2024-03-19 | 中国航发四川燃气涡轮研究院 | 航空发动机主轴疲劳试验标准循环载荷谱构建方法及系统 |
CN117725802B (zh) * | 2024-02-07 | 2024-04-16 | 中国航发四川燃气涡轮研究院 | 航空发动机主轴疲劳试验标准循环载荷谱构建方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109520717B (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mansour et al. | Mechanical and dynamic behavior of fused filament fabrication 3D printed polyethylene terephthalate glycol reinforced with carbon fibers | |
Song et al. | Imperfection sensitivity of thin elastic cylindrical shells subject to partial axial compression | |
EP2983099B1 (en) | Method for determining reduction factor of axial load bearing capacity of a cylindrical shell structure in a rocket | |
CN109520717A (zh) | 一种发动机主轴疲劳试验载荷确定方法 | |
Bisagni | Composite cylindrical shells under static and dynamic axial loading: An experimental campaign | |
CN109614715B (zh) | 一种多轴载荷作用下考虑缺口效应的能量场强法及其应用 | |
Li et al. | Experiments on the geometrically nonlinear vibration of a thin-walled cylindrical shell with points supported boundary condition | |
Martini et al. | Elastodynamic behavior of balanced closed-loop mechanisms: numerical analysis of a four-bar linkage | |
CN111339700B (zh) | 核电汽轮机叶片的疲劳损伤评估方法、装置和存储介质 | |
Hilburger et al. | Buckling behavior of compression-loaded quasi-isotropic curved panels with a circular cutout | |
Becker et al. | Vibration isolation for parameter-varying rotor systems using piezoelectric actuators and gain-scheduled control | |
Wertz et al. | An energy-based torsional-shear fatigue lifing method | |
Wu et al. | Crankshaft High Cycle Bending Fatigue Research Based on a 2D Simplified Model and Different Mean Stress Models | |
JP2001050863A (ja) | 軸受試験装置及び軸受試験方法 | |
JP5721226B2 (ja) | 多軸疲労寿命評価方法 | |
CN107291988A (zh) | 一种动量轮安装界面等效激励力获取方法 | |
Livieri et al. | An approximation in closed form for the integral of Oore–Burns for cracks similar to a star domain | |
Sadati et al. | Identification of a nonlinear joint in an elastic structure using optimum equivalent linear frequency response function | |
Peridas et al. | Crack initiation: the choice of tests available to calibrate Dang Van's criterion | |
Wu et al. | Numerical simulation of fatigue damage and shape instability behavior of steel 40Cr by the damage-coupled crystal plastic model | |
CN112824867A (zh) | 一种测试金属材料疲劳极限的方法以及装置 | |
CN112035781B (zh) | 橡胶隔振器耐久寿命预测方法、装置、设备及存储介质 | |
CN104677383B (zh) | 一种直接输出频率谐振式陀螺研究系统 | |
Demenko et al. | Evaluation of spacecraft composite structure load-carrying capability by shock test | |
Zhang et al. | FEA based Dissipation Energy and Temperature Distribution of Rubber Bushing |
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 |