CN109883709B - 一种基于相对等效应变的随机多轴热机计数方法 - Google Patents

一种基于相对等效应变的随机多轴热机计数方法 Download PDF

Info

Publication number
CN109883709B
CN109883709B CN201910174125.9A CN201910174125A CN109883709B CN 109883709 B CN109883709 B CN 109883709B CN 201910174125 A CN201910174125 A CN 201910174125A CN 109883709 B CN109883709 B CN 109883709B
Authority
CN
China
Prior art keywords
point
equivalent strain
relative
strain
heat engine
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
Application number
CN201910174125.9A
Other languages
English (en)
Other versions
CN109883709A (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201910174125.9A priority Critical patent/CN109883709B/zh
Publication of CN109883709A publication Critical patent/CN109883709A/zh
Application granted granted Critical
Publication of CN109883709B publication Critical patent/CN109883709B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种基于相对等效应变的随机多轴热机计数方法,该方法通过相对等效应变把多轴热机转化为单轴热机,基于单轴情况下的雨流计数原理,只计数相对等效应变从零到达最高点的半循环(反复),然后在计算蠕变损伤时把大反复上的小反复剔除出去,使得每个时间历程对应的蠕变损伤只被计算一次,避免了反复计算一个时间历程的蠕变损伤。该方法使用方便,便于编程,能够广泛应用于变幅多轴热机计数。通过验证,采用该方法进行的变幅多轴热机械疲劳寿命估算取得较好的预测效果。

Description

一种基于相对等效应变的随机多轴热机计数方法
技术领域
本发明涉及疲劳强度领域,特指一种基于相对等效应变的随机多轴热机计数方法。
背景技术
很多工程构件,例如航空发动机涡轮盘、燃气轮机等在服役时不仅受到多轴载荷,而且往往是在变幅、变温环境下,从而产生不可逆的热机械疲劳损伤。而这些热疲劳损伤往往使关键部位零件失效,大大降低设备的可靠性和使用寿命,因此迫切需要解决随机载荷下的疲劳寿命预测问题。对于变幅情况下的寿命预测,首先要解决的问题就是多轴计数。因此,解决随机热机情况下载荷谱的多轴计数问题,就可以真正的解决其疲劳寿命预测问题,从而解决实际工程所遇到的问题。
目前有很多关于变幅多轴计数的研究,但是多数都是在常温情况下,并不涉及热机械载荷。由于常温下没有蠕变等与时间相关的损伤,所以在计数时并不需要考虑时间的影响,但热机情况下的多轴计数却需要把时间因素考虑进去。因此,在热机计数时综合考虑蠕变等与时间相关的因素考量进去,从而能够很合理的计算蠕变损伤,最终使总的热机疲劳损伤的计算更加合理。因此,研究一种基于相对等效应变的随机多轴热机计数方法具有重要意义。
发明内容
本发明目的在于针对变幅多轴热机械疲劳的发展要求,提出了一种基于相对等效应变的随机多轴热机计数方法。
本发明所采用的技术方案为一种基于相对等效应变的随机多轴热机计数方法,该方法的实现步骤如下:
步骤1):确定热机循环历程中的最大或峰值等效应变点以及其等效应变值;
步骤2):以此点作为第一个最大值参考点,计算相对于此点的等效应变点的各点的相对值;
步骤3):新的相对载荷时间历程从此点开始,此处原为最大等效应变点,现为相对等效应变为零的起始点。
步骤4):计算各处的相对于起始点的等效应变趋势直到开始下降的点,一次类推,直到找到相对于起始点的最大等效应变点;
步骤5):计数出每个反复的时间历程,当一个反复里包含其他的反复,称这个反复为大反复,被包含的反复称为小反复。计算大反复的蠕变损伤时,剔除小反复所造成的蠕变损伤;
与现有技术相比,本发明具有如下有益效果:
本发明提出一种基于相对等效应变的随机多轴热机计数方法。该方法通过相对等效应变把多轴热机转化为单轴热机,基于单轴情况下的雨流计数原理,只计数相对等效应变从零到达最高点的半循环(反复),然后在计算蠕变损伤时把大反复上的小反复剔除出去,使得每个时间历程对应的蠕变损伤只被计算一次,避免了反复计算一个时间历程的蠕变损伤。该方法使用方便,便于编程,能够广泛应用于变幅多轴热机计数。通过验证,采用该方法进行的变幅多轴热机械疲劳寿命估算取得较好的预测效果。
附图说明
图1为变幅多轴热机计数方法流程图。
图2为该方法的计数结果。
图3为利用本发明所提出的多轴热机计数方法进行寿命预测,预测结果与试验结果比较。
具体实施方式
一种基于相对等效应变的随机多轴热机计数方法,具体实施方式如下:
步骤1):确定热机循环历程中的最大或峰值等效应变点M以及其等效应变值
Figure BDA0001988980420000031
步骤2):对于任意一点A,以等效应变点M作为第一个最大值参考点,计算相对于点M的等效应变点的各点的相对值,计算公式如下:
Figure BDA0001988980420000032
Figure BDA0001988980420000033
其中,v为依赖于弹塑性应变比的有效泊松比,
Figure BDA0001988980420000034
为A点的相对等效应变。εxy分别为x、y方向上的应变;εij表示A点的应变,εij max表示参考点M的应变值;i、j分别表示三阶张量矩阵的各个行列;γ表示剪切应变。
步骤3):确定新的载荷时间历程。新的相对载荷时间历程从点M开始,此处原为最大等效应变点,现为相对等效应变为零的起始点。
步骤4):对于新的载荷时间历程,计算各处的相对于起始点的等效应变趋势直到开始下降的点,一次类推,直到找到相对于起始点的最大等效应变点;
步骤5):计数出每个反复的时间历程,当一个反复里包含其他的反复,称这个反复为大反复,被包含的反复称为小反复。计算大反复的蠕变损伤时,剔除小反复所造成的蠕变损伤。大反复的蠕变损伤计算公式如下:
Figure BDA0001988980420000041
其中,
Figure BDA0001988980420000042
是大反复上剔除小反复之后所造成的蠕变损伤,
Figure BDA0001988980420000043
是大反复上整个时间历程所造成的蠕变损伤,
Figure BDA0001988980420000044
是大反复上的小反复所造成的蠕变损伤。

Claims (2)

1.一种基于相对等效应变的随机多轴热机计数方法,其特征在于:该方法的实现步骤如下,
步骤1):确定热机循环历程中的最大或峰值等效应变点M以及其等效应变值
Figure FDA0003072220420000011
步骤2):以此点作为第一个最大值参考点,计算相对于此点的等效应变点的各点的相对值,计算公式如下:
Figure FDA0003072220420000012
Figure FDA0003072220420000013
εxy分别为x、y方向上的应变;εij表示任意时刻的应变,εij max表示任意时刻最大应变值;i、j分别表示二阶张量矩阵的各个行列;γ表示剪切应变;v为依赖于弹塑性应变比的有效泊松比;
步骤3):新的相对载荷时间历程从此点开始,此处原为最大等效应变点,现为相对等效应变为零的起始点;
步骤4):计算各处的相对于起始点的等效应变趋势直到开始下降的点,依次类推,直到找到相对于起始点的最大等效应变点;步骤5):计数出每个反复的时间历程,当一个反复里包含其他的反复,称这个反复为大反复,被包含的反复称为小反复;计算大反复的蠕变损伤时,剔除小反复所造成的蠕变损伤。
2.根据权利要求1所述的一种基于相对等效应变的随机多轴热机计数方法,其特征在于:所述步骤5)中在计算蠕变损伤时,虽然大反复包含小反复上的数据,但小反复的蠕变损伤不计入大反复的蠕变损伤。
CN201910174125.9A 2019-03-08 2019-03-08 一种基于相对等效应变的随机多轴热机计数方法 Active CN109883709B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910174125.9A CN109883709B (zh) 2019-03-08 2019-03-08 一种基于相对等效应变的随机多轴热机计数方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910174125.9A CN109883709B (zh) 2019-03-08 2019-03-08 一种基于相对等效应变的随机多轴热机计数方法

Publications (2)

Publication Number Publication Date
CN109883709A CN109883709A (zh) 2019-06-14
CN109883709B true CN109883709B (zh) 2021-06-25

Family

ID=66931268

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910174125.9A Active CN109883709B (zh) 2019-03-08 2019-03-08 一种基于相对等效应变的随机多轴热机计数方法

Country Status (1)

Country Link
CN (1) CN109883709B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111159871B (zh) * 2019-12-23 2024-03-26 北京工业大学 一种基于路径曲线积分的随机多轴循环计数方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6294623B2 (ja) * 2013-10-01 2018-03-14 住友ゴム工業株式会社 高分子材料の寿命予測方法
CN103926152B (zh) * 2014-04-09 2016-08-24 北京工业大学 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法
CN107423540A (zh) * 2017-04-09 2017-12-01 北京工业大学 一种基于权平均最大剪切应变幅平面的多轴疲劳寿命预测方法
CN109426685B (zh) * 2017-08-28 2023-09-22 长城汽车股份有限公司 零部件疲劳分析方法

Also Published As

Publication number Publication date
CN109883709A (zh) 2019-06-14

Similar Documents

Publication Publication Date Title
CN111624116B (zh) 基于权平均最大剪切应力平面的疲劳寿命预测方法及装置
CN111090957B (zh) 一种高温结构危险点应力-应变计算方法
CN105628511B (zh) 一种高温合金蠕变寿命预测方法
CN109142083A (zh) 一种变载荷历程下的蠕变损伤计算方法及模型
CN103926152A (zh) 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法
CN102663201B (zh) 一种考虑电气互连结构可靠性的电子产品可靠性预计方法
CN108984977B (zh) 一种确定多阶段任务焊点振动损伤累积规则的有限元方法
CN115640666B (zh) 一种基于损伤等效的航空发动机加速任务试车谱编制方法
CN115329544B (zh) 一种多轴变幅载荷下的疲劳寿命预测方法
CN109883709B (zh) 一种基于相对等效应变的随机多轴热机计数方法
CN108254250B (zh) 一种考虑动态应变时效的热机多轴应力应变关系确定方法
CN111881603A (zh) 一种考虑失效相关性的机械结构疲劳可靠性评估方法
CN110096841B (zh) 多轴热机械加载下缺口根部应力-应变状态评估方法
CN107748817B (zh) 一种考虑非比例附加强化的高温多轴本构关系确定方法
CN114139276A (zh) 一种盘轴一体式整体叶盘结构疲劳寿命分析方法
CN113109190B (zh) 一种多轴热机械载荷下基于短裂纹的寿命预测方法
Dileep et al. Effect of fatigue damage parameter on the cumulative life of a turbine rotor under multiaxial loading
Nakane et al. Thermal ratcheting of solder-bonded layered plates: cyclic recovery and growth of deflection
O’Nora et al. Constitutive modeling of TMF and creep-fatigue of a Ni-base alloy
Shao et al. The Principle and application of physics-of-failure based reliability technology
CN109829222A (zh) 一种宽应变率宽温域的复合推进剂压缩力学本构方法
CN114021357A (zh) 基于应力概率密度法的结构随机疲劳寿命估算方法
CN115952628B (zh) 一种计算航空发动机涡轮盘盘心等效应力的方法
SHI et al. Constitutive Modeling of Hardening and Creep Response of a Nickel-Based Superalloy Udimet 720Li
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