CN109883709A - 一种基于相对等效应变的随机多轴热机计数方法 - Google Patents
一种基于相对等效应变的随机多轴热机计数方法 Download PDFInfo
- Publication number
- CN109883709A CN109883709A CN201910174125.9A CN201910174125A CN109883709A CN 109883709 A CN109883709 A CN 109883709A CN 201910174125 A CN201910174125 A CN 201910174125A CN 109883709 A CN109883709 A CN 109883709A
- Authority
- CN
- China
- Prior art keywords
- point
- equivalent strain
- relative
- heat engine
- strain
- 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
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 2
- 230000006735 deficit Effects 0.000 abstract 3
- 238000005050 thermomechanical fatigue Methods 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000930 thermomechanical effect Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于相对等效应变的随机多轴热机计数方法,该方法通过相对等效应变把多轴热机转化为单轴热机,基于单轴情况下的雨流计数原理,只计数相对等效应变从零到达最高点的半循环(反复),然后在计算蠕变损伤时把大反复上的小反复剔除出去,使得每个时间历程对应的蠕变损伤只被计算一次,避免了反复计算一个时间历程的蠕变损伤。该方法使用方便,便于编程,能够广泛应用于变幅多轴热机计数。通过验证,采用该方法进行的变幅多轴热机械疲劳寿命估算取得较好的预测效果。
Description
技术领域
本发明涉及疲劳强度领域,特指一种基于相对等效应变的随机多轴热机计数方法。
背景技术
很多工程构件,例如航空发动机涡轮盘、燃气轮机等在服役时不仅受到多轴载荷,而且往往是在变幅、变温环境下,从而产生不可逆的热机械疲劳损伤。而这些热疲劳损伤往往使关键部位零件失效,大大降低设备的可靠性和使用寿命,因此迫切需要解决随机载荷下的疲劳寿命预测问题。对于变幅情况下的寿命预测,首先要解决的问题就是多轴计数。因此,解决随机热机情况下载荷谱的多轴计数问题,就可以真正的解决其疲劳寿命预测问题,从而解决实际工程所遇到的问题。
目前有很多关于变幅多轴计数的研究,但是多数都是在常温情况下,并不涉及热机械载荷。由于常温下没有蠕变等与时间相关的损伤,所以在计数时并不需要考虑时间的影响,但热机情况下的多轴计数却需要把时间因素考虑进去。因此,在热机计数时综合考虑蠕变等与时间相关的因素考量进去,从而能够很合理的计算蠕变损伤,最终使总的热机疲劳损伤的计算更加合理。因此,研究一种基于相对等效应变的随机多轴热机计数方法具有重要意义。
发明内容
本发明目的在于针对变幅多轴热机械疲劳的发展要求,提出了一种基于相对等效应变的随机多轴热机计数方法。
本发明所采用的技术方案为一种基于相对等效应变的随机多轴热机计数方法,该方法的实现步骤如下:
步骤1):确定热机循环历程中的最大或峰值等效应变点以及其等效应变值;
步骤2):以此点作为第一个最大值参考点,计算相对于此点的等效应变点的各点的相对值;
步骤3):新的相对载荷时间历程从此点开始,此处原为最大等效应变点,现为相对等效应变为零的起始点。
步骤4):计算各处的相对于起始点的等效应变趋势直到开始下降的点,一次类推,直到找到相对于起始点的最大等效应变点;
步骤5):计数出每个反复的时间历程,当一个反复里包含其他的反复,称这个反复为大反复,被包含的反复称为小反复。计算大反复的蠕变损伤时,剔除小反复所造成的蠕变损伤;
与现有技术相比,本发明具有如下有益效果:
本发明提出一种基于相对等效应变的随机多轴热机计数方法。该方法通过相对等效应变把多轴热机转化为单轴热机,基于单轴情况下的雨流计数原理,只计数相对等效应变从零到达最高点的半循环(反复),然后在计算蠕变损伤时把大反复上的小反复剔除出去,使得每个时间历程对应的蠕变损伤只被计算一次,避免了反复计算一个时间历程的蠕变损伤。该方法使用方便,便于编程,能够广泛应用于变幅多轴热机计数。通过验证,采用该方法进行的变幅多轴热机械疲劳寿命估算取得较好的预测效果。
附图说明
图1为变幅多轴热机计数方法流程图。
图2为该方法的计数结果。
图3为利用本发明所提出的多轴热机计数方法进行寿命预测,预测结果与试验结果比较。
具体实施方式
一种基于相对等效应变的随机多轴热机计数方法,具体实施方式如下:
步骤1):确定热机循环历程中的最大或峰值等效应变点M以及其等效应变值
步骤2):对于任意一点A,以等效应变点M作为第一个最大值参考点,计算相对于点M的等效应变点的各点的相对值,计算公式如下:
其中,v为依赖于弹塑性应变比的有效泊松比,为A点的相对等效应变。εx,εy分别为x、y方向上的应变;εij表示A点的应变,εij max表示参考点M的应变值;i、j分别表示三阶张量矩阵的各个行列;γ表示剪切应变。
步骤3):确定新的载荷时间历程。新的相对载荷时间历程从点M开始,此处原为最大等效应变点,现为相对等效应变为零的起始点。
步骤4):对于新的载荷时间历程,计算各处的相对于起始点的等效应变趋势直到开始下降的点,一次类推,直到找到相对于起始点的最大等效应变点;
步骤5):计数出每个反复的时间历程,当一个反复里包含其他的反复,称这个反复为大反复,被包含的反复称为小反复。计算大反复的蠕变损伤时,剔除小反复所造成的蠕变损伤。大反复的蠕变损伤计算公式如下:
其中,是大反复上剔除小反复之后所造成的蠕变损伤,是大反复上整个时间历程所造成的蠕变损伤,是大反复上的小反复所造成的蠕变损伤。
Claims (2)
1.一种基于相对等效应变的随机多轴热机计数方法,其特征在于:该方法的实现步骤如下,
步骤1):确定热机循环历程中的最大或峰值等效应变点M以及其等效应变值
步骤2):以此点作为第一个最大值参考点,计算相对于此点的等效应变点的各点的相对值,计算公式如下:
εx,εy分别为x、y方向上的应变;εij表示任意时刻的应变,εij max表示任意时刻最大应变值;i、j分别表示三阶张量矩阵的各个行列;γ表示剪切应变;
步骤3):新的相对载荷时间历程从此点开始,此处原为最大等效应变点,现为相对等效应变为零的起始点;
步骤4):计算各处的相对于起始点的等效应变趋势直到开始下降的点,一次类推,直到找到相对于起始点的最大等效应变点;
步骤5):计数出每个反复的时间历程,当一个反复里包含其他的反复,称这个反复为大反复,被包含的反复称为小反复;计算大反复的蠕变损伤时,剔除小反复所造成的蠕变损伤。
2.根据权利要求1所述的一种基于相对等效应变的随机多轴热机计数方法,其特征在于:所述步骤5)中在计算蠕变损伤时,虽然大反复包含小反复上的数据,但小反复的蠕变损伤不计入大反复的蠕变损伤。
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 true CN109883709A (zh) | 2019-06-14 |
CN109883709B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111159871A (zh) * | 2019-12-23 | 2020-05-15 | 北京工业大学 | 一种基于路径曲线积分的随机多轴循环计数方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926152A (zh) * | 2014-04-09 | 2014-07-16 | 北京工业大学 | 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法 |
JP2015072141A (ja) * | 2013-10-01 | 2015-04-16 | 住友ゴム工業株式会社 | 高分子材料の寿命予測方法 |
CN107423540A (zh) * | 2017-04-09 | 2017-12-01 | 北京工业大学 | 一种基于权平均最大剪切应变幅平面的多轴疲劳寿命预测方法 |
CN109426685A (zh) * | 2017-08-28 | 2019-03-05 | 长城汽车股份有限公司 | 零部件疲劳分析方法 |
-
2019
- 2019-03-08 CN CN201910174125.9A patent/CN109883709B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015072141A (ja) * | 2013-10-01 | 2015-04-16 | 住友ゴム工業株式会社 | 高分子材料の寿命予測方法 |
CN103926152A (zh) * | 2014-04-09 | 2014-07-16 | 北京工业大学 | 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法 |
CN107423540A (zh) * | 2017-04-09 | 2017-12-01 | 北京工业大学 | 一种基于权平均最大剪切应变幅平面的多轴疲劳寿命预测方法 |
CN109426685A (zh) * | 2017-08-28 | 2019-03-05 | 长城汽车股份有限公司 | 零部件疲劳分析方法 |
Non-Patent Citations (2)
Title |
---|
孙斌斌: "多轴随机载荷循环计数及疲劳寿命预测方法研究", 《中国优秀硕士学位论文全文数据库》 * |
斯海维: "《结构与材料的疲劳》", 31 January 2014, 航空工业出版社 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111159871A (zh) * | 2019-12-23 | 2020-05-15 | 北京工业大学 | 一种基于路径曲线积分的随机多轴循环计数方法 |
CN111159871B (zh) * | 2019-12-23 | 2024-03-26 | 北京工业大学 | 一种基于路径曲线积分的随机多轴循环计数方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109883709B (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111624116B (zh) | 基于权平均最大剪切应力平面的疲劳寿命预测方法及装置 | |
CN112364535B (zh) | 一种基于晶体塑性的蠕变疲劳寿命预测方法 | |
CN109885874A (zh) | 一种基于abaqus的多轴蠕变疲劳预测方法 | |
CN111090957B (zh) | 一种高温结构危险点应力-应变计算方法 | |
CN109142083A (zh) | 一种变载荷历程下的蠕变损伤计算方法及模型 | |
CN103926152A (zh) | 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法 | |
CN110987675B (zh) | 一种基于临界面损伤的多轴变幅热机械疲劳寿命预测方法 | |
CN110220805B (zh) | 一种基于蠕变疲劳损伤累积的变幅多轴热机疲劳寿命预测方法 | |
CN111475972B (zh) | 一种高温疲劳-蠕变交互损伤评估方法、系统及存储介质 | |
CN114547943A (zh) | 火箭发动机阀门的寿命计算方法、装置及电子设备 | |
CN110705019A (zh) | 一种高温蠕变损伤等效加速方法 | |
CN108254250B (zh) | 一种考虑动态应变时效的热机多轴应力应变关系确定方法 | |
CN114186444A (zh) | 航空发动机高温部件的载荷谱编制方法、介质、终端、应用 | |
CN109883709B (zh) | 一种基于相对等效应变的随机多轴热机计数方法 | |
CN106653099B (zh) | 一种聚变堆偏滤器的疲劳寿命分析评判方法 | |
CN118364666A (zh) | 基于热弹塑性损伤本构模型的微动疲劳寿命预测方法 | |
CN110096841B (zh) | 多轴热机械加载下缺口根部应力-应变状态评估方法 | |
CN113109190B (zh) | 一种多轴热机械载荷下基于短裂纹的寿命预测方法 | |
CN109357957B (zh) | 一种基于极值窗口的疲劳监测计数方法 | |
Roy et al. | A stochastic model for evolution of creep damage in engineering materials | |
Dileep et al. | Effect of fatigue damage parameter on the cumulative life of a turbine rotor under multiaxial loading | |
CN114139276A (zh) | 一种盘轴一体式整体叶盘结构疲劳寿命分析方法 | |
O’Nora et al. | Constitutive modeling of TMF and creep-fatigue of a Ni-base alloy | |
CN114152527B (zh) | 一种基于监测数据的三维雨流疲劳分析方法 | |
SHI et al. | Constitutive Modeling of Hardening and Creep Response of a Nickel-Based Superalloy Udimet 720Li |
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 |