CN113514350B - 适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法 - Google Patents

适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法 Download PDF

Info

Publication number
CN113514350B
CN113514350B CN202110325998.2A CN202110325998A CN113514350B CN 113514350 B CN113514350 B CN 113514350B CN 202110325998 A CN202110325998 A CN 202110325998A CN 113514350 B CN113514350 B CN 113514350B
Authority
CN
China
Prior art keywords
creep
strain
fatigue
stress
damage
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
CN202110325998.2A
Other languages
English (en)
Other versions
CN113514350A (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.)
Nanjing Tech University
Original Assignee
Nanjing Tech University
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 Nanjing Tech University filed Critical Nanjing Tech University
Priority to CN202110325998.2A priority Critical patent/CN113514350B/zh
Publication of CN113514350A publication Critical patent/CN113514350A/zh
Application granted granted Critical
Publication of CN113514350B publication Critical patent/CN113514350B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/0071—Creep
    • 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)

Abstract

本发明公开了一种适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，首先在同一温度下开展多组不同的应变控纯疲劳试验以及蠕变试验；然后确定应变幅值与寿命的关系以及非弹性应变能密度与非弹性应变能密度率的关系；接着开展多组应变控以及应力应变混合控蠕变疲劳试验；进而根据应变幅值确定每周的疲劳损伤，根据保载时的应力演化或应变演化，确定每周蠕变损伤；基于损伤预测所述材料在蠕变疲劳载荷下的寿命；最后，根据试样的累积疲劳和蠕变损伤以及美国ASME III‑NH标准规定的双线性损伤准则，判定试样在任意时刻的损伤状态。本发明可以很好地预测材料在应变控以及混合控蠕变疲劳载荷下的寿命，还可以判定试样在任意时刻的损伤状态。

Description

适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法
技术领域
本发明涉及寿命预测和损伤评定领域，尤其涉及适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法。
背景技术
在航空和电力领域中，众多关键高温核心部件往往承受着复杂的蠕变、疲劳以及蠕变疲劳交互的损伤作用。通常，频繁开停车产生的温度变化使得这些部件处于应变控制的疲劳载荷中，而装备持续稳定运行过程中的恒定内压或者离心力，使高温部件承受恒应力控制的蠕变载荷作用。因此，目前实验室采用的传统应变控制蠕变疲劳试验很难反映上述高温部件真实复杂的运行工况。另外，由于传统应变控蠕变疲劳载荷是通过应变保载期间的应力松弛引入蠕变损伤，导致蠕变损伤容易饱和，难以形成蠕变损伤主导的蠕变疲劳失效模式。所以有必要引入一种能够产生宽范围蠕变疲劳损伤的应力应变混合控制蠕变疲劳载荷加载波形。
蠕变疲劳载荷下的损伤识别是装备设计和制造最为关心的问题之一。目前常用的是线性损伤累积方法。美国ASME III-NH(锅炉及压力容器设计规范)以及法国RCC-MR(核岛机械设施设计和建造规范)均采用这种方法来评估循环加载过程中的蠕变和疲劳损伤，其中蠕变损伤用时间分数表示，疲劳损伤用寿命分数表示。然而，不同于传统的应变控蠕变疲劳载荷，在应力应变混合控制蠕变疲劳载荷下，应力保载期间的蠕变变形显著且应变幅值随循环周次不断变化，这些势必会对现有的寿命预测及损伤评定方法提出挑战，因此需要发明一种适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，以满足不同载荷作用下蠕变疲劳损伤状态的识别。
发明内容
本发明的目的在于提供一种适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，可以实现材料在应变控制以及应力应变混合控制蠕变疲劳载荷下的寿命预测及损伤评估，且具有适用性广，精度高，操作简单等优点。
实现本发明目的的技术解决方案为：一种适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，包括以下步骤：
步骤S1，取相同材料的若干试样，在同一温度下对部分试样进行同应变速率不同应变幅值的应变控制疲劳试验，获取疲劳试验数据，转入步骤S2；在与应变控制疲劳试验相同的温度下对部分试样开展不同应力的蠕变试验，获取蠕变试验数据，转入步骤S3；
步骤S2，据疲劳试验数据，确定任意应变幅值下疲劳试验的循环寿命，转入步骤S4；
步骤S3，根据蠕变试验数据，确定蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的定量关系，转入步骤S4；
步骤S4，在同一温度及同一应变速率下，开展多组应变控制蠕变疲劳试验以及应力应变混合控制蠕变疲劳试验，获得蠕变疲劳试验的半寿命周次，这些数据将用于验证本发明的准确性，同时转入步骤S5和步骤S6；
步骤S5，以蠕变疲劳试验的半寿命周次为特征周次，结合步骤S2中得到的任意应变幅值下疲劳试验的循环寿命，根据蠕变疲劳试验半寿命周次的应变幅值确定每一周的疲劳损伤，转入步骤S7；
步骤S6，以蠕变疲劳试验的半寿命周次为特征周次，结合步骤S3中得到的蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的定量关系，根据蠕变疲劳试验半寿命周次应变保载期间应力演化或应力保载期间应变演化，确定非弹性应变能密度及非弹性应变能密度率，从而得到每一周的蠕变损伤，转入步骤S7；
步骤S7，利用线性损伤累积法则，根据前述每一周的疲劳损伤以及每一周的蠕变损伤预测该材料在蠕变疲劳载荷下的循环寿命，转入步骤S8；
步骤S8，根据每一周的疲劳和蠕变损伤，确定循环加载过程中累积蠕变损伤和累积疲劳损伤，基于确定的累积疲劳损伤和累积蠕变损伤以及美国ASME III-NH标准规定的双线性蠕变疲劳包络线，判定材料在蠕变疲劳任意时刻下的损伤状态：若损伤状态点在双线性包络线以外，则判定为材料失效，反之则为未失效。
本发明与现有技术相比，其显著优点在于：
(1)本发明基于线性损伤累积法则，耦合能量的观点，并且考虑平均应力的作用，有明确的物理意义；
(2)本发明操作简单，计算方便，适用性广，能够准确预测不同材料在应变控以及应力应变混合控蠕变疲劳载荷下的循环寿命；
(3)基于本发明提出的蠕变和疲劳损伤，不同材料所有的蠕变疲劳失效数据点几乎都落在标准规定的包络线以外。
(4)本发明可以判定材料在应变控以及应力应变混合控蠕变疲劳载荷的任意时刻下的损伤状态。
附图说明
图1为本发明的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法的流程图。
图2为本发明一个实施例应力应变混合控制蠕变疲劳载荷下的应变-时间以及应力-时间关系图。
图3为本发明实施例1的预测寿命结果图。
图4为实施例的损伤交互图中的线性损伤累积法则、双线性损伤累积法则示意图。
图5为基于现有准则的蠕变疲劳损伤交互图。
图6为本发明实施例2的蠕变疲劳损伤交互图。
具体实施方式
以下结合附图，以具体实施例对本发明作进一步详细说明。应该理解，以下实施例仅用于说明本发明而非用于限定本发明的范围。
本发明所述的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，基于线性损伤累积，认为总损伤等于累积疲劳损伤与累积蠕变损伤之和。本发明在此基础上，基于应变能耗散的观点，考虑应力应变混合控制蠕变疲劳载荷下非线性蠕变变形以及变幅值疲劳的特点，提出了同时适用于应变控制以及应力应变混合控制蠕变疲劳载荷的蠕变及疲劳损伤的表达式。基于该发明提出的蠕变及疲劳损伤，不仅能够预测不同蠕变疲劳载荷下的循环寿命还可以判定材料在应变控以及应力应变混合控蠕变疲劳载荷下任意时刻的损伤状态。
本发明可以很好地预测不同材料在应变控以及应力应变混合控蠕变疲劳载荷下的循环寿命，且基于本发明的累积蠕变及累积疲劳损伤，可以判定蠕变疲劳试样在任意时刻的损伤状态，结果表明不同控制模式下的蠕变疲劳失效试样，其蠕变疲劳损伤预测值均落在标准规定的包络线以外，证明了本方法相比传统方法的优越性。
结合图1，本发明所述的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，包括以下的步骤：
步骤S1，取相同材料的若干试样，在同一温度下对部分试样进行同应变速率不同应变幅值的应变控制疲劳试验，获取疲劳试验数据【国标GB/T26077-2010】，转入步骤S2；在与应变控制疲劳试验相同的温度下对部分试样开展不同应力的蠕变试验【国标GB/T 2039-2012】，获取蠕变试验数据，转入步骤S3。
步骤S2，根据疲劳试验数据，确定任意应变幅值下疲劳试验的循环寿命，具体如下：
建立应变幅值与应变控制疲劳寿命的关系，如式(1)所示：
Figure BDA0002994688490000041
Figure BDA0002994688490000042
Figure BDA0002994688490000043
式(1)中，等式右边第一项表示弹性应变幅值Δεe与寿命Nf的关系，其中E表示材料在试验温度下的弹性模量，由试验得到的弹性应变幅值Δεe与寿命Nf，拟合得疲劳强度系数σ′和疲劳强度指数n1。式(1)等式右边第二项表示非弹性应变幅值Δεin与寿命Nf的关系，如式(3)，同样，根据试验的非弹性应变幅值Δεin与寿命Nf，通过拟合得疲劳延性系数ε′和疲劳延性指数m1。
基于上述应变幅值与应变控制疲劳寿命的关系，即可求出任意应变幅值下疲劳试验的循环寿命。
转入步骤S4。
步骤S3，根据蠕变试验数据，确定蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的定量关系，具体如下：
根据蠕变试验数据，建立蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的关系，如式(4)所示：
Figure BDA0002994688490000051
win＝σ·εf (5)
Figure BDA0002994688490000052
其中蠕变失效应变能系数C和蠕变失效应变能指数b由所有蠕变试验失效时的非弹性应变能密度与非弹性应变能密度率计算得到，式(4)中蠕变失效时的非弹性应变能密度win表示为式(5)，即蠕变应力σ与蠕变断裂应变εf的乘积；蠕变失效时的非弹性应变能密度率
Figure BDA0002994688490000053
表示为式(6)，即失效时的非弹性应变能密度win除以蠕变断裂时间tf。
转入步骤S4。
步骤S4，在同一温度及同一应变速率下，开展多组应变控制蠕变疲劳试验以及应力应变混合控制蠕变疲劳试验，获得蠕变疲劳试验的半寿命周次，这些数据将用于验证本发明的准确性，同时转入步骤S5和步骤S6。
其中，应变控制蠕变疲劳试验采用应变控制的梯形波进行加载【美国标准ASTME2714-13】，保载时间td为输入参量，保载应力σd随保载时间演化，应力应变混合控制蠕变疲劳试验的疲劳加载部分采用应变控制，应力应变混合控制蠕变疲劳试验的蠕变加载部分采用应力控制，图2给出了应力应变混合控制蠕变疲劳载荷下的应变-时间以及应力-时间关系图。
首先以试验输入的应变速率
Figure BDA0002994688490000054
加载到试验输入的目标拉伸应变幅值+εa；接着以同样的试验输入应变速率
Figure BDA0002994688490000055
进行卸载；当应力达到试验输入保载应力σd时停止卸载，与此同时试验控制模式切换至应力控制并进行应力为σd的恒应力保载，当保载时间达到试验输入的保载时间td时，保载结束。试验控制模式切换至应变控制，以同样的试验输入应变速率
Figure BDA0002994688490000056
进行压缩，压缩至试验输入的目标压缩应变幅值-εa。
步骤S5，以蠕变疲劳试验的半寿命周次为特征周次，结合步骤S2中得到的任意应变幅值下疲劳试验的循环寿命，根据蠕变疲劳试验半寿命周次的应变幅值确定每一周的疲劳损伤，具体如下：
以蠕变疲劳试验的半寿命周次为特征周次，则蠕变疲劳试验每一周的疲劳损伤
Figure BDA0002994688490000061
表示为式(7)，
Figure BDA0002994688490000062
其中N为应变控制纯疲劳试验的循环寿命，该应变控制纯疲劳试验的应变幅值等于蠕变疲劳试验特征周次的应变幅值，即蠕变疲劳试验半寿命周次的应变幅值，根据式(1)建立的应变幅值与应变控制疲劳寿命关系可求得N。
转入步骤S7。
步骤S6，以蠕变疲劳试验的半寿命周次为特征周次，结合步骤S3中得到的蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的定量关系，根据蠕变疲劳试验半寿命周次应变保载期间应力演化或应力保载期间应变演化，确定非弹性应变能密度及非弹性应变能密度率，从而得到每一周的蠕变损伤，具体如下：
步骤S61，应变控制蠕变疲劳试验应变保载期间的保载应力与时间的演化σd(t)以及混合控制蠕变疲劳试验应力保载期间的保载应变与时间的演化εd(t)均满足式(8)的非线性演化规律，
Figure BDA0002994688490000063
Figure BDA0002994688490000064
其中，应力松弛系数p1和应力松弛指数q1通过拟合应变保载期间的应力松弛曲线得到，蠕变变形系数p2和蠕变变形指数q2通过拟合应力保载期间的蠕变曲线得到；本发明在保证精度的前提下，提出新的应力松弛曲线以及蠕变曲线的表达式，方便本发明的推广及工程应用。
步骤S62，考虑平均应力作用，保载阶段t时刻的非弹性应变能密度
Figure BDA0002994688490000065
如式(9-1)和式(9-2)所示，
Figure BDA0002994688490000066
其中E表示材料在试验温度下的弹性模量，σ0表示蠕变疲劳试验半寿命周次保载开始时的应力，σm表示蠕变疲劳试验半寿命周次的平均应力；σ表示保载阶段的蠕变应力。
步骤S63，通过应变控制蠕变疲劳试验应变保载阶段的保载应力与时间的演化σd(t)、混合控制蠕变疲劳试验应力保载阶段的保载应变与时间的演化εd(t)，以及保载阶段t时刻的非弹性应变能密度
Figure BDA0002994688490000071
得到非弹性应变能密度率
Figure BDA0002994688490000072
即对式(9-1)和式(9-2)进行求导，并与式(8-1)和式(8-2)联立，得到式(10-1)和式(10-2)：
Figure BDA0002994688490000073
σd表示混合控蠕变疲劳载荷的保载应力。
初始保载阶段较高非弹性应变速率会造成对蠕变损伤的高估。为解决这一问题，现有技术提出设置非弹性应变速率阈值来限制蠕变损伤，但阈值的确定是建立在大量试验数据的基础上的，这大大提高了寿命预测的成本。
步骤S64，本发明将最小的非弹性应变能密度率，即保载结束时的非弹性应变能密度率作为特征非弹性应变能密度率，在考虑蠕变损伤高估的基础上，极大地简化了预测方法，利于工程应用及推广。因此本发明中蠕变疲劳试验每一周的蠕变损伤
Figure BDA0002994688490000074
可以表示为式(11)，其中
Figure BDA0002994688490000075
表示蠕变疲劳失效时材料的非弹性应变能密度。
Figure BDA0002994688490000076
其中
Figure BDA0002994688490000077
表示蠕变疲劳失效时材料的非弹性应变能密度。
步骤S65，根据蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的关系，以及非弹性应变能密度率
Figure BDA0002994688490000078
得到最终的蠕变疲劳试验每一周的蠕变损伤
Figure BDA0002994688490000079
如式(12)所示：
Figure BDA00029946884900000710
Figure BDA00029946884900000711
转入步骤S7。
步骤S7，利用线性损伤累积法则，根据前述每一周的疲劳损伤以及每一周的蠕变损伤预测该材料在蠕变疲劳载荷下的循环寿命Nc-f可以表示为式(13)
Figure BDA0002994688490000081
转入步骤S8。
步骤S8，根据每一周的疲劳损伤和蠕变损伤，确定蠕变疲劳试样加载过程中的累积疲劳损伤Df(n)以及累积蠕变损伤Dc(n)：
Figure BDA0002994688490000082
Figure BDA0002994688490000083
其中n表示循环周次数，i表示当前循环周次数。
将确定的累积疲劳损伤和累积蠕变损伤绘制于蠕变疲劳损伤评定图中，其中每一个损伤状态点横坐标为累积疲劳损伤，纵坐标为累积蠕变损伤。发现失效试样的损伤状态点均位于美国ASME III-NH标准(锅炉及压力容器设计规范，2019版)规定的双线性损伤包络线之外，从而验证了本发明损伤方法的准确性和合理性。接着将其他承受蠕变疲劳载荷试样对应的损伤状态点绘制于损伤评定图上，判定材料在此蠕变疲劳载荷任意时刻下的损伤状态：若损伤状态点在双线性包络线以外，则判定为材料失效，反之则为未失效。
实施例1采用本发明的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法预测了不同温度以及不同材料在应变控以及应力应变混合控蠕变疲劳载荷下的循环寿命。
选取数据为15组P92钢在625℃下的应力应变混合控制蠕变疲劳试验，保载时间为300s,600s,1800s，保载应力为140MPa-180MPa，其试验方法如说明书中所述。4组P92钢在650℃下的应变控制蠕变疲劳试验，以及5组P92钢在600℃下的应变控制蠕变疲劳试验，试验方法如说明书中所述。不同温度以及应变幅值下P92钢的疲劳试验采用发表论文的数据【Wang X,Zhang W,Zhang T,et al.A New Empirical Life Prediction Model for 9–12％Cr Steels under Low Cycle Fatigue and Creep Fatigue Interaction Loadings[J].Metals,2019,9(2).和Zhang T,Wang X,Ji Y,et al.P92 steel creep-fatigueinteraction responses under hybrid stress-strain controlled loading and alife prediction model[J].International Journal of Fatigue,2020:105837.】。不同温度以及不同应力下P92钢的蠕变试验采用论文数据【Kimura K,Takahashi Y.Evaluationof Long-Term Creep Strength of ASME Grades 91,92,and 122Type Steels[C]//ASME2012Pressure Vessels and Piping Conference.2012.】。9组P91钢在550℃下的应变控制蠕变疲劳试验以及纯蠕变和纯疲劳试验采用论文数据【Asayama T.Update and ImproveSubsection NH-Alternative Simplified Creep-Fatigue Design Methods[R].ASMEStandards Technology,LLC.2009.】。
首先按照本发明的步骤S2，根据所述的疲劳试验数据，计算出任意应变幅值下疲劳试验的循环寿命，然后根据步骤S3，确定蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的定量关系，接着按照步骤S5，以蠕变疲劳试验的半寿命周次为特征周次，根据半寿命周次的应变幅值确定每一周的疲劳损伤；按照步骤S6，计算非弹性应变能密度及非弹性应变能密度率，最终得到每一周的蠕变损伤。最后利用线性累积损伤法则，预测出不同材料在不同试验载荷下的循环寿命。
由图3的结果可见，无论是应变控制蠕变疲劳还是应力应变混合控制的蠕变疲劳试验，利用本发明预测的寿命均在2倍误差带以内。因此本发明的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法对蠕变疲劳载荷，试验材料，试验温度都具有普适性。
实施例2采用本发明的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，判定材料在应变控以及应力应变混合控蠕变疲劳载荷的任意时刻下的损伤状态。
线性损伤累积没有考虑蠕变疲劳交互的作用，常常会低估蠕变疲劳损伤的情况。因此在现行蠕变疲劳设计中，引入了双线性设计准则。对于9Cr-1Mo材料，ASME III-NH(锅炉及压力容器设计规范)定义累积蠕变损伤以及累积疲劳损伤均达到0.3为双线性准则中的交互转折点，如图4所示。本实施例选取数据与实施例1相同，基于本发明提出的蠕变及疲劳损伤表达式，本实施例判定材料在应变控以及应力应变混合控蠕变疲劳载荷的任意时刻下的损伤状态。
首先按照本发明的步骤S2，根据所述的疲劳试验数据，计算出任意应变幅值下疲劳试验的循环寿命，然后根据步骤S3，确定蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的定量关系，接着按照步骤S5，以蠕变疲劳试验的半寿命周次为特征周次，根据半寿命周次的应变幅值确定每一周的疲劳损伤；按照步骤S6，计算非弹性应变能密度及非弹性应变能密度率，最终得到每一周的蠕变损伤。按照步骤S8，计算蠕变疲劳加载过程中的累积蠕变损伤以及累积疲劳损伤，验证本发明损伤方法的准确性和合理性。最后，取多个试验损伤状态点，判定材料在该蠕变疲劳以及该时刻下的损伤状态。
由图5和图6的结果可见，基于本发明提出的蠕变及疲劳损伤，大部分的蠕变疲劳失效点几乎都落在标准规定的包络线以外，而基于ASME III-NH标准中蠕变及疲劳损伤，相当多的失效点落在标准规定的包络线以内，尤其对于应力应变混合控蠕变疲劳载荷，大部分失效点都落在标准规定的包络线以内。因此验证本发明的蠕变疲劳损伤是合理的且准确性高于美国ASME III-NH标准提出的损伤方法。
除此以外，由图6可见，多个试验损伤状态点位于双线性蠕变疲劳包络线以内，表明这些试验点对应的试样未发生失效，与实际损伤状态相符，因此本发明可以判定材料在该蠕变疲劳以及该时刻下的损伤状态。
从实施例1和实施例2的结果看出：采用本发明的方法，可以很好地预测不同材料在应变控以及应力应变混合控蠕变疲劳载荷下的循环寿命，且基于本发明的蠕变及疲劳损伤，几乎所有的蠕变疲劳失效点都落在标准规定的包络线以外，同时本发明还可以判定材料在蠕变疲劳载荷下任意时刻下的损伤状态。

Claims (7)

1.一种适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，其特征在于，包括以下步骤：
步骤S1，取相同材料的若干试样，在同一温度下对部分试样进行同应变速率不同应变幅值的应变控制疲劳试验，获取疲劳试验数据，转入步骤S2；在与应变控制疲劳试验相同的温度下对部分试样开展不同应力的蠕变试验，获取蠕变试验数据，转入步骤S3；
步骤S2，根据疲劳试验数据，确定任意应变幅值下疲劳试验的循环寿命，转入步骤S4；
步骤S3，根据蠕变试验数据，确定蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的定量关系，转入步骤S4；
步骤S4，在同一温度及同一应变速率下，开展多组应变控制蠕变疲劳试验以及应力应变混合控制蠕变疲劳试验，获得蠕变疲劳试验的半寿命周次，同时转入步骤S5和步骤S6；
步骤S5，以蠕变疲劳试验的半寿命周次为特征周次，结合步骤S2中得到的任意应变幅值下疲劳试验的循环寿命，根据蠕变疲劳试验半寿命周次的应变幅值确定其每一周的疲劳损伤，转入步骤S7；
步骤S6，以蠕变疲劳试验的半寿命周次为特征周次，结合步骤S3中得到的蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的定量关系，根据蠕变疲劳试验半寿命周次应变保载期间应力演化或应力保载期间应变演化，确定非弹性应变能密度及非弹性应变能密度率，从而得到每一周的蠕变损伤，具体如下：
步骤S61，应变控制蠕变疲劳试验应变保载阶段的保载应力与时间的演化σd(t)以及混合控制蠕变疲劳试验应力保载阶段的保载应变与时间的演化εd(t)分别满足下式的非线性演化规律，
Figure FDA0003539232370000011
Figure FDA0003539232370000012
其中，应力松弛系数p1和应力松弛指数q1通过拟合应变保载期间的应力松弛曲线得到，蠕变变形系数p2和蠕变变形指数q2通过拟合应力保载期间的蠕变曲线得到；
步骤S62，考虑平均应力作用，保载阶段t时刻的非弹性应变能密度
Figure FDA0003539232370000021
如式(9-1)和式(9-2)所示，
Figure FDA0003539232370000022
其中E表示材料在试验温度下的弹性模量，σ0表示蠕变疲劳试验半寿命周次保载开始时的应力，σm表示蠕变疲劳试验半寿命周次的平均应力；σ表示保载阶段的蠕变应力；
步骤S63，通过应变控制蠕变疲劳试验应变保载阶段的保载应力与时间的演化σd(t)、混合控制蠕变疲劳试验应力保载阶段的保载应变与时间的演化εd(t)，以及保载阶段t时刻的非弹性应变能密度
Figure FDA0003539232370000023
得到非弹性应变能密度率
Figure FDA0003539232370000024
Figure FDA0003539232370000025
σd表示混合控蠕变疲劳载荷的保载应力；
步骤S64，将最小的非弹性应变能密度率，即保载结束时的非弹性应变能密度率作为特征非弹性应变能密度率，因此蠕变疲劳试验每一周的蠕变损伤
Figure FDA0003539232370000026
表示为式(11)，
Figure FDA0003539232370000027
其中
Figure FDA0003539232370000028
表示蠕变疲劳失效时材料的非弹性应变能密度，td为保载时间；
步骤S65，根据蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的关系，以及非弹性应变能密度率
Figure FDA0003539232370000029
得到最终的蠕变疲劳试验每一周的蠕变损伤
Figure FDA00035392323700000210
如式(12)所示：
Figure FDA00035392323700000211
Figure FDA00035392323700000212
转入步骤S7；
步骤S7，利用线性损伤累积法则，根据前述每一周的疲劳损伤以及每一周的蠕变损伤预测该材料在蠕变疲劳载荷下的循环寿命Nc-f表示为式(13)：
Figure FDA0003539232370000031
转入步骤S8；
步骤S8，根据每一周的疲劳和蠕变损伤，确定循环加载过程中累积蠕变损伤和累积疲劳损伤，基于确定的累积疲劳损伤和累积蠕变损伤以及美国ASME III-NH标准规定的双线性蠕变疲劳包络线，判定材料在蠕变疲劳任意时刻下的损伤状态：
若损伤状态点在双线性包络线以外，则判定为材料失效，反之则为未失效。
2.根据权利要求1所述的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，其特征在于：所述步骤S2中，根据疲劳试验数据，建立应变幅值与应变控制疲劳寿命的关系，如公式(1)，
Figure FDA0003539232370000032
Figure FDA0003539232370000033
Figure FDA0003539232370000034
式(1)中，等式右边第一项表示弹性应变幅值Δεe与寿命Nf的关系，其中E表示材料在试验温度下的弹性模量，由试验得到的弹性应变幅值Δεe与寿命Nf，拟合得疲劳强度系数σ′和疲劳强度指数n1；式(1)等式右边第二项表示非弹性应变幅值Δεin与寿命Nf的关系，如式(3)，同样，根据试验的非弹性应变幅值Δεin与寿命Nf，通过拟合得疲劳延性系数ε′和疲劳延性指数m1；
基于上述应变幅值与应变控制疲劳寿命的关系，求出任意应变幅值下疲劳试验的循环寿命。
3.根据权利要求1所述的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，其特征在于：所述步骤S3中，根据蠕变试验数据，建立蠕变失效时的非弹性应变能密度与非弹性应变能密度率之间的关系，如式(4)所示：
Figure FDA0003539232370000041
win＝σ·εf (5)
Figure FDA0003539232370000042
其中，蠕变失效应变能系数C和蠕变失效应变能指数b由所有蠕变试验失效时的非弹性应变能密度与非弹性应变能密度率计算得到，式(4)中蠕变失效时的非弹性应变能密度win表示为式(5)，即蠕变应力σ与蠕变断裂应变εf的乘积；蠕变失效时的非弹性应变能密度率
Figure FDA0003539232370000043
表示为式(6)，即失效时的非弹性应变能密度win除以蠕变断裂时间tf。
4.根据权利要求1所述的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，其特征在于：所述步骤S4中，应变控制蠕变疲劳试验采用应变控制的梯形波进行加载，保载时间td为输入参量，保载应力σd随保载时间演化；应力应变混合控制蠕变疲劳试验的疲劳加载部分采用应变控制，蠕变加载部分采用应力控制。
5.根据权利要求4所述的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，其特征在于，应力应变混合控制蠕变疲劳试验的疲劳加载部分采用应变控制，保载蠕变阶段采用应力控制，具体如下：
首先以试验输入的应变速率
Figure FDA0003539232370000044
加载到试验输入的目标拉伸应变幅值+εa；接着以同样的试验输入应变速率
Figure FDA0003539232370000045
进行卸载；当应力达到试验输入保载应力σd时停止卸载，与此同时试验控制模式切换至应力控制并进行应力为σd的恒应力保载，当保载时间达到试验输入的保载时间td时，保载结束；试验控制模式切换至应变控制，以同样的试验输入应变速率
Figure FDA0003539232370000046
进行压缩，压缩至试验输入的目标压缩应变幅值-εa。
6.根据权利要求1所述的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，其特征在于：所述步骤S5中，以蠕变疲劳试验的半寿命周次为特征周次，确定蠕变疲劳试验每一周的疲劳损伤
Figure FDA0003539232370000047
Figure FDA0003539232370000048
其中N为应变控制纯疲劳试验的循环寿命，该应变控制纯疲劳试验的应变幅值等于蠕变疲劳试验特征周次的应变幅值，即蠕变疲劳试验半寿命周次的应变幅值，根据式(1)建立的应变幅值与应变控制疲劳寿命关系求得N。
7.根据权利要求1所述的适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法，其特征在于，所述步骤S8中，根据每一周的疲劳损伤和蠕变损伤，确定蠕变疲劳试样加载过程中的累积疲劳损伤Df(n)以及累积蠕变损伤Dc(n)：
Figure FDA0003539232370000051
Figure FDA0003539232370000052
其中n表示循环周次数；
基于确定的累积疲劳损伤Df(n)、累积蠕变损伤Dc(n)以及美国ASME III-NH标准规定的双线性蠕变疲劳包络线，判定材料在蠕变疲劳任意时刻下的损伤状态：若损伤状态点在双线性包络线以外，则判定为材料失效，反之则为未失效。
CN202110325998.2A 2021-03-26 2021-03-26 适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法 Active CN113514350B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110325998.2A CN113514350B (zh) 2021-03-26 2021-03-26 适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110325998.2A CN113514350B (zh) 2021-03-26 2021-03-26 适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法

Publications (2)

Publication Number Publication Date
CN113514350A CN113514350A (zh) 2021-10-19
CN113514350B true CN113514350B (zh) 2022-05-13

Family

ID=78062096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110325998.2A Active CN113514350B (zh) 2021-03-26 2021-03-26 适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法

Country Status (1)

Country Link
CN (1) CN113514350B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114295491B (zh) * 2021-12-14 2024-06-11 南京航空航天大学 一种蠕变损伤及变形随时间演化行为的预测方法
CN114239175B (zh) * 2021-12-17 2022-11-15 天津大学 基于非线性蠕变-疲劳交互作用的高温结构件寿命预测方法
CN115017745B (zh) * 2022-08-09 2022-10-21 天津大学 一种用于火电机组结构件的蠕变疲劳寿命预测方法及系统
CN115982958B (zh) * 2022-12-07 2023-10-13 南京工业大学 一种基于工程损伤力学的材料蠕变疲劳寿命预测方法
CN115931532B (zh) * 2022-12-30 2023-10-13 南京工业大学 基于修正时间分数法的电力高温部件损伤状态判定方法
CN116052813B (zh) * 2022-12-30 2023-10-13 南京工业大学 一种全数值驱动的实时蠕变疲劳损伤监测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1457994A (en) * 1974-12-02 1976-12-08 Allied Chem Method of determining dynamic strains in composite structures
CN103344515A (zh) * 2013-07-05 2013-10-09 北京航空航天大学 一种基于局部应力应变法的低周疲劳和高强度冲击耦合的损伤计算方法
CN105158084A (zh) * 2015-09-15 2015-12-16 华东理工大学 一种材料的蠕变-疲劳寿命预测方法
CN111428361A (zh) * 2020-03-24 2020-07-17 南京工业大学 一种适用于多种循环载荷的寿命预测方法
CN111879636A (zh) * 2020-07-31 2020-11-03 华东理工大学 一种材料的蠕变-疲劳-氧化实时损伤累积评估方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1457994A (en) * 1974-12-02 1976-12-08 Allied Chem Method of determining dynamic strains in composite structures
CN103344515A (zh) * 2013-07-05 2013-10-09 北京航空航天大学 一种基于局部应力应变法的低周疲劳和高强度冲击耦合的损伤计算方法
CN105158084A (zh) * 2015-09-15 2015-12-16 华东理工大学 一种材料的蠕变-疲劳寿命预测方法
CN111428361A (zh) * 2020-03-24 2020-07-17 南京工业大学 一种适用于多种循环载荷的寿命预测方法
CN111879636A (zh) * 2020-07-31 2020-11-03 华东理工大学 一种材料的蠕变-疲劳-氧化实时损伤累积评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
保载条件下FGH95材料的疲劳特性及寿命建模;魏大盛等;《航空动力学报》;20070328(第03期);第425-430页 *
基于能量密度耗散准则的蠕变—疲劳寿命预测模型及应用;王润梓;《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》;20190815;第C031-4页 *

Also Published As

Publication number Publication date
CN113514350A (zh) 2021-10-19

Similar Documents

Publication Publication Date Title
CN113514350B (zh) 适用于判定应力应变混合控蠕变疲劳损伤状态的通用方法
CN107843510B (zh) 基于室温布氏硬度预测超临界机组t/p91耐热钢剩余持久寿命评估方法
Pejkowski et al. Stress-strain response and fatigue life of four metallic materials under asynchronous loadings: Experimental observations
CN110274826A (zh) 一种基于单轴疲劳s-n曲线的硬质金属材料多轴高周疲劳失效预测方法
CN111428361A (zh) 一种适用于多种循环载荷的寿命预测方法
CN107843509B (zh) 基于室温布氏硬度预测超临界机组t/p92耐热钢剩余持久寿命评估方法
CN112730061B (zh) 一种多级变温变载蠕变寿命评价方法
CN112284921A (zh) 基于高温液压鼓胀试样测试的材料单轴应力-应变关系确定方法
Hyde et al. Some considerations on specimen types for small sample creep tests
CN112417580A (zh) 一种考虑保载效应的高温低周疲劳寿命预测方法
JP4676323B2 (ja) 金属材料の破壊靭性値測定方法
CN111859619B (zh) 一种利用硬度对热老化材料低周疲劳寿命的无损预测方法
Stelzer et al. Comparison of quasi-static and cyclic fatigue delamination resistance of carbon fiber reinforced polymer-matrix laminates under different mode loading
McMurtrey Report on the FY18 creep rupture and creep-fatigue tests on the first commercial heat of Alloy 709
CN110220784B (zh) 相变诱导塑性钢马氏体相变强度增量表征方法及系统
Cheng et al. Geometric discontinuity effect on creep-fatigue behaviors in a nickel-based superalloy hole structure considering ratcheting deformation
Socie Multiaxial fatigue damage assessment
CN115931532B (zh) 基于修正时间分数法的电力高温部件损伤状态判定方法
CN114065576B (zh) 一种蠕变疲劳载荷下缺口部件的缺口效应评价方法
Dzioba et al. Temperature dependency of fracture toughness of high-strength ferritic steel Hardox-400
WO2003054521A2 (en) Method for tasting structural materials fatigue
Pohja et al. Recommendation for Creep and Creep-fatigue assessment for P91 Components
Song et al. Numerical evaluation of fatigue crack growth rate of hot-rolled 16Mn steel based on compliance method
Schopf et al. Investigations on Multi-Stage Tests and Transient Endurance Limit Behavior Under Low-, High-and Very High Cycle Fatigue Loads
CN210923422U (zh) 温度应力耦合条件下合金应力腐蚀开裂敏感性评定装置

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