CN106202913B - 时间相关的蠕变疲劳损伤评定方法 - Google Patents
时间相关的蠕变疲劳损伤评定方法 Download PDFInfo
- Publication number
- CN106202913B CN106202913B CN201610531611.8A CN201610531611A CN106202913B CN 106202913 B CN106202913 B CN 106202913B CN 201610531611 A CN201610531611 A CN 201610531611A CN 106202913 B CN106202913 B CN 106202913B
- Authority
- CN
- China
- Prior art keywords
- mrow
- damage
- creep
- fatigue
- msub
- 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
- 230000006735 deficit Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000006378 damage Effects 0.000 claims abstract description 167
- 238000009825 accumulation Methods 0.000 claims abstract description 77
- 230000003993 interaction Effects 0.000 claims abstract description 37
- 230000001186 cumulative effect Effects 0.000 claims abstract description 33
- 239000000463 material Substances 0.000 claims abstract description 24
- 238000002474 experimental method Methods 0.000 claims abstract description 23
- 239000013077 target material Substances 0.000 claims abstract description 19
- 238000012417 linear regression Methods 0.000 claims description 29
- 238000012360 testing method Methods 0.000 claims description 22
- 230000035939 shock Effects 0.000 claims description 9
- 238000009661 fatigue test Methods 0.000 claims description 5
- 230000014759 maintenance of location Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 230000000750 progressive effect Effects 0.000 abstract description 3
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 8
- 238000013461 design Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 4
- 229910052759 nickel Inorganic materials 0.000 description 4
- 229910000601 superalloy Inorganic materials 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 239000007769 metal material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000009975 flexible effect Effects 0.000 description 1
- 239000013521 mastic Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000009183 running Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种时间相关的蠕变疲劳损伤评定方法，包括步骤：对一目标材料进行多次实验获取所需的一参数组；根据参数组计算获得各周次的蠕变损伤dc和各周次的疲劳损伤df；计算第n周次时的累积蠕变损伤Dc(n)和第n周次时的累积疲劳损伤Df(n)；根据所述第n周次时的累积蠕变损伤Dc(n)和第n周次时的累积疲劳损伤Df(n)绘制一蠕变疲劳损伤交互图；结合所述蠕变疲劳损伤交互图和至少一损伤累积法则评估所述目标材料的实时累积损伤情况。本发明的一种时间相关的蠕变疲劳损伤评定方法，能够更好地实现材料在蠕变疲劳交互作用下的实时损伤累积评估，具有精确度高、灵活性强、可实施评估和直观的优点。
Description
技术领域
本发明涉及损伤评定领域，尤其涉及一种时间相关的蠕变疲劳损伤评定方法
背景技术
在能源动力、石油化工和航空航天等领域中，许多结构部件长期运行在高温条件下受应变波形控制，循环失效寿命往往不足105次，即，在高温低周疲劳载荷下，其寿命往往受到蠕变、疲劳和蠕变疲劳交互作用等多种机制的制约。在蠕变疲劳交互作用的条件下，材料的损伤评估非常重要，而利用线性损伤累积准则进行较为准确的实时损伤累积评估是亟需突破的重点和难点。
自20世纪70年代以来，国内外学者推演出了许多蠕变疲劳交互作用下的损伤评估方法。但是在这些损伤评估方法中，基于线性损伤累积准则，即可以实现损伤的每周次实时评估。经典的利用线性损伤累积准则并写入美国ASME设计准则(American Society ofMechanical Engineers.Boiler and Pressure Vessel Code,Section III,subsection-NH 2005)和英国R5设计准则(British Energy.Assessment procedure for the hightemperature response of structures.R5Issue 3,2003)的模型有时间分数法和延性耗竭法。在时间分数法中，主要以应力水平σ和蠕变断裂时间tR为基本参量推演蠕变损伤公式；在延性耗竭法中，主要以蠕变延性εf和非弹性应变率为基本参量推演蠕变损伤公式。在两种方法中，疲劳损伤的定义则为固定值1/N0，其中N0为特定工况纯疲劳的寿命。
但是这两种经典的方法都有各自的缺陷。在时间分数法中，计算损伤的结果通常偏于非保守；而在延性耗竭法中，计算损伤的结果又过于保守。并且，在这两种模型中，疲劳损伤仅为一个固定的数值，对疲劳损伤的表述过于简单粗略。
发明内容
针对上述现有技术中的不足，本发明提供一种时间相关的蠕变疲劳损伤评定方法，能够更好地实现材料在蠕变疲劳交互作用下的实时损伤累积评估，具有精确度高、灵活性强、可实施评估和直观的优点。
为了实现上述目的，本发明一种时间相关的蠕变疲劳损伤评定方法，包括步骤：
S1：对一目标材料进行多次实验获取所需的一参数组；
S2：根据所述参数组计算获得各周次的蠕变损伤dc和各周次的疲劳损伤df；
S3：根据所述各周次的蠕变损伤dc计算第n周次时的累积蠕变损伤Dc(n)和第n周次时的累积疲劳损伤Df(n)，n为大于0的自然数；
S4：以所述第n周次时的累积蠕变损伤Dc(n)为纵坐标并以和第n周次时的累积疲劳损伤Df(n)为横坐标绘制一蠕变疲劳损伤交互图；
S5：结合所述蠕变疲劳损伤交互图和至少一损伤累积法则评估所述目标材料的实时累积损伤情况。
进一步地，所述S2步骤包括利用公式(1)和公式(2)计算获得各周次的蠕变损伤dc和各周次的疲劳损伤df：
其中，dc为蠕变损伤，t为时间，D为第一线性回归参数，n1为第一线性回归幂指数，wf0为平台值，为非弹性应变能密度耗散率，σ0为最大应力，σm为平均应力，Δεp为塑性应变范围参数，A为第一松弛参数，B为第二松弛参数；E为弹性模量；th为蠕变疲劳保持时间；
其中，df为疲劳损伤，α为第二线性回归参数，β为第二线性回归幂指数。
进一步地，在所述S3步骤中，根据公式(3)计算所述第n周次时的累积蠕变损伤Dc(n)和所述第n周次时的累积疲劳损伤Df(n)：
其中，dc(i)为第i周次的蠕变损伤，df(i)为第i周次的疲劳损伤，i为大于零的自然数。
进一步地，所述损伤累积法则包括一线性损伤累积法则，所述线性损伤累积法则的数学表达式为公式(4)：
Dc+Df＝1 (4)；
其中，Dc为累积蠕变损伤，Df为累积疲劳损伤，当公式(4)成立时判定所述目标材料失效。
进一步地，所述损伤累积法则还包括一双线性损伤累积法则，所述双线性损伤累积法则的数学表达式为公式(5)：
其中，所述双线性损伤累积法则的折点坐标为当公式(5)成立时判定所述目标材料失效。
进一步地，所述损伤累积法则还包括一简化的非线性累积损伤法则，所述简化的非线性累积损伤法则的数学表达式为公式(6)：
Dc m+Df m＝1 (6)；
其中，m为简化的非线性损伤累积法则系数，当公式(6)成立时判定所述目标材料失效。
进一步地，所述S1步骤进一步包括步骤：
S11：通过一不同应力水平下的高温蠕变试验，获得所述目标材料的每个应力水平下的一蠕变延性数据εf、一蠕变断裂时间tR、所述第一线性回归参数D、所述第一线性回归幂指数n1和所述平台值wf0；
S12：通过一相同温度不同应变水平的应力松弛试验，获得每个应变水平下应力与松弛时间的变化关系，并根据所述变化关系获得所述第一松弛参数A和所述第二松弛参数B；
S13：通过一相同温度不同应变幅的应变控制下的纯疲劳试验，获得各应变幅下的所述塑性应变范围参数Δεp和疲劳寿命参数N0，根据纯疲劳实验条件下的所述塑性应变范围参数Δεp和疲劳寿命参数N0获得所述第二线性回归参数α和第二线性回归幂指数β；
S14：通过一相同温度的应变控制下的蠕变疲劳试验，记录各周次的应力应变关系，获得各周次的所述最大应力σ0、各周次的所述平均应力σm以及各周次的蠕变疲劳实验条件下的所述塑性应变范围参数Δεp。
进一步地，所述步骤S12中，根据所述变化关系和公式(7)计算获得所述第一松弛参数A和所述第二松弛参数B：
σ(t)＝σ0-(A·logΔεp+B)·log(1+t) (7)；
其中，t为时间，σ(t)为t时的应力水平。
进一步地，所述步骤S13中，根据所述塑性应变范围参数Δεp和疲劳寿命参数N0并通过公式(8)计算获得所述第二线性回归参数α和第二线性回归幂指数β：
N0＝α·(Δεp)β (8)；
其中，N0为疲劳寿命参数，Δεp为塑性应变范围参数。
本发明由于采用了以上技术方案，使其具有以下有益效果：
(1)本发明利用应变能耗散的概念，具有清晰的物理意义，并有较为精确的预测结果；
(2)本发明利用损伤交互图，能够实时评估材料的累积损伤，也可以直观观察材料剩余的蠕变疲劳持久性能；
(3)本发明具有很强的灵活性，在不同的部件、不同的设计要求下，可以利用不同的损伤累积法则，进行实时累积损伤的评估。
附图说明
图1为本发明实施例的时间相关的蠕变疲劳损伤评定方法的流程图；
图2为本发明实施例的实验中的应力应变关系图；
图3为本发明实施例的在自然坐标损伤交互图中的线性损伤累积法则、双线性损伤累积法则以及简化的累积损伤法则示意图；
图4为本发明实施例的高温镍基GH4169超合金650℃下不同应力水平下高温蠕变试验参数拟合结果图；
图5为本发明实施例的高温镍基GH4169超合金650℃下不同应变水平的应力松弛试验参数拟合结果图；
图6为本发明实施例的高温镍基GH4169超合金650℃下不同应变幅下的纯疲劳试验参数拟合结果图；
图7为本发明实施例一的对数坐标下损伤交互图中的蠕变疲劳实时损伤累积评估结果图；
图8为本发明实施例二的对数坐标下损伤交互图中的蠕变疲劳实时损伤累积评估结果图；
图9为本发明实施例三的对数坐标下损伤交互图中的蠕变疲劳实时损伤累积评估结果图。
具体实施方式
下面根据附图1-9，给出本发明的较佳实施例，并予以详细描述，使能更好地理解本发明的功能、特点。
请参阅图1，本发明的一种时间相关的蠕变疲劳损伤评定方法：
首先，进行S1步骤：使用高温蠕变试验机，对一目标材料进行多次实验获取所需的一参数组。
其中，S1步骤进一步包括步骤：
S11：使用高温蠕变试验机，通过一不同应力水平下的高温蠕变试验，获得目标材料的每个应力水平下的一蠕变延性数据εf、一蠕变断裂时间tR、第一线性回归参数D、第一线性回归幂指数n1和平台值wf0；如以应变百分比数值为横轴，应力为纵轴作一应力应变关系图，第一线性回归参数D为应力应变关系图中实线的截距、n1线性回归幂指数为蠕变能量耗散图中实线的斜率；
S12：使用高温蠕变试验机，通过一相同温度不同应变水平的应力松弛试验，获得每个应变水平下应力与松弛时间的变化关系，根据变化关系和公式(7)即应力松弛公式，计算获得第一松弛参数A和第二松弛参数B：
σ(t)＝σ0-(A·logΔεp+B)·log(1+t) (7)；
其中，t为时间，σ(t)为t时的应力水平。
S13：使用高温蠕变试验机，通过一相同温度不同应变幅的应变控制下的纯疲劳试验，获得各应变幅下的塑性应变范围参数Δεp和疲劳寿命参数N0，根据纯疲劳条件下塑性应变范围参数Δεp和疲劳寿命参数N0并通过公式(8)计算获得第二线性回归参数α和第二线性回归幂指数β：
N0＝α·(Δεp)β (8)；
其中，N0为疲劳寿命参数，Δεp为塑性应变范围参数。
S14：通过一相同温度的应变控制下的蠕变疲劳试验，记录各周次的应力应变关系，获得各周次的最大应力σ0、各周次的平均应力σm以及各周次的蠕变疲劳条件下塑性应变范围参数Δεp。
接着，进行S2步骤：在蠕变疲劳工况下，根据参数组并利用公式(1)和公式(2)计算获得各周次的蠕变损伤dc和各周次的疲劳损伤df：
其中，dc为蠕变损伤，t为时间，D为第一线性回归参数，n1为第一线性回归幂指数，wf0为平台值，为非弹性应变能密度耗散率，σ0为最大应力，σm为平均应力，Δεp为塑性应变范围参数，A为第一松弛参数，B为第二松弛
参数；E为弹性模量；th为蠕变疲劳保持时间；
df为疲劳损伤，α为第二线性回归参数，β为第二线性回归幂指数，以应变百分比数值为横轴，应力为纵轴作一应力应变关系图，应力应变关系图可如图2所示，显示了某一个周次下的蠕变疲劳试验中的应力应变关系。
请继续参加图1，其中，失效应变能密度wf与非弹性应变能密度耗散率之间的函数关系表示为公式(9)：
D为第一线性回归参数，为非弹性应变能密度耗散率，n1为第一线性回归幂指数，其中，
wf＝σ·εf (1O)；
其中，σ表示蠕变试验中施加的应力水平，εf表示蠕变试验中获取的蠕变延性，tR表示蠕变试验中获取的蠕变断裂时间。
然后，进行S3步骤：根据各周次的蠕变损伤dc计算第n周次时的累积蠕变损伤Dc(n)和第n周次时的累积疲劳损伤Df(n)，n为大于0的自然数；
具体地，根据公式(3)计算第n周次时的累积蠕变损伤Dc(n)和第n周次时的累积疲劳损伤Df(n)：
其中，dc(i)为第i周次的蠕变损伤，df(i)为第i周次的疲劳损伤，i为大于零的自然数。
之后，进行S4步骤：根据第n周次时的累积蠕变损伤Dc(n)和第n周次时的累积疲劳损伤Df(n)绘制一蠕变疲劳损伤交互图。
最后，进行S5步骤：结合蠕变疲劳损伤交互图和至少一损伤累积法则评估目标材料的实时累积损伤情况。
其中，损伤累积法则包括一线性损伤累积法则，线性损伤累积法则的数学表达式为公式(4)：
Dc+Df＝1 (4)；
其中，Dc为累积蠕变损伤，Df为累积疲劳损伤，当公式(4)成立时判定目标材料失效。
损伤累积法则还包括一双线性损伤累积法则，双线性损伤累积法则的数学表达式为公式(5)：
其中，双线性损伤累积法则的折点坐标为当公式(5)成立时判定目标材料失效。
损伤累积法则还包括一简化的非线性累积损伤法则，简化的非线性累积损伤法则的数学表达式为公式(6)：
Dc m+Df m＝1 (6)；
其中，m为简化的非线性损伤累积法则系数，当公式(6)成立时判定目标材料失效。
针对不同的材料、部件以及设计要求，选取不同的非线性损伤累积法则，根据美国ASME以及英国R5准则，越是重要的部件以及苛刻的设计要求，所选取的损伤累积法则越是保守，当进行累积损伤的评估时，设计准则可以表示在在蠕变疲劳损伤交互图中。
请参阅图3，图示了一种在自然坐标损伤交互图中的线性损伤累积法则、双线性损伤累积法则以及简化的累积损伤法则示意图。
本发明的一种材料的蠕变疲劳实时损伤累积评估方法基于线性损伤累积准则，结合应变能密度耗散的概念，并引用蠕变疲劳损伤交互设计图。不同于时间分数法和延性耗竭法，该方法主要以失效应变能密度wf和非弹性应变能耗散率为基本参量，并且修正了每周次疲劳损伤df的表达式。通过蠕变疲劳损伤交互图，该方法能够实时评估材料的累积损伤，也可以直观观察材料剩余的蠕变疲劳持久性能。
以下实施例1～3中，采用本发明的材料蠕变疲劳实时损伤累积评估方法对650℃下的高温镍基GH4169超合金进行验证。通过10个应力水平，650℃下的不同应力水平下高温蠕变试验结果，并利用公式(9)～(11)求出的纯蠕变试验中的三个参数：σ蠕变试验中施加的应力水平，εf蠕变试验中获取的蠕变延性，tR蠕变试验中获取的蠕变断裂时间。如图4所示；通过2个应变水平，650℃下的不同应变水平的应力松弛试验结果，求出公式(7)中应力松弛公式中的两个参数A和B，如图5所示；通过4个应变水平，650℃下的不同应变幅下的纯疲劳试验，求出式(2)中每周次疲劳损伤df的两个参数α和β。如图6所示。另外，该材料在650℃下的弹性模量为161.64MPa。
实施例1：
选取数据为总应变范围为2.0％，在应变最大处的保载时间为120s，温度为650℃下的轴向等幅蠕变疲劳交互试验，参照标准为GB/T 15248-2008《金属材料轴向等幅低循环疲劳试验方法》。其蠕变疲劳实验寿命为255周次，实时损伤累积评估所用的参数包括：失效应变能密度wf、非弹性应变能密度耗散率蠕变试验中施加的应力水平σ、时间t、疲劳寿命参数N0、塑性应变范围参数Δεp，各参数的关系请参阅图4、图5和图6。
实验中可得塑性应变范围参数Δεp、最大应力σ0以及平均应力σm随着循环周次n的变化，根据公式(1)依次计算每个周次的蠕变损伤dc，根据公式(2)依次计算每个周次的疲劳损伤df，然后根据公式(3)计算循环至某个周次n的累积损伤，其中，第n周次时的累积蠕变损伤Dc(n)表示为在损伤交互图中某个周次n的纵坐标，第n周次时的累积疲劳损伤Df(n)表示为损伤交互图中的横坐标。每个循环周次坐标点在蠕变疲劳损伤交互图中所连成轨迹线则为该材料在该工况下实时的累积损伤，当轨迹线与损伤累积法则线相交时，则判定为在该工况、该设计准则下材料失效，即不能进行实际工况的服役。在科学研究中，一般采用经典的线性损伤累积法则进行寿命评估和预测。但是在实际工程运用中，双线性损伤准则考虑了蠕变疲劳的交互作用，这使得该准则更为保守也更加适用于工程。在核电以及压力管道领域中，双线性损伤准则得到广泛使用。21世纪以来，英国能源科技评估研究小组提出了简化的损伤累积法则，它因更为保守、更为简明的表达而被广泛认可。依据美国ASME和英国R5准则，这里给出的公式(5)中双线性累积损伤法则的折点推荐坐标值为(0.3,0.3)，公式(6)中简化的非线性损伤累积法则的系数m＝0.576。
图7表示了本实例中在蠕变疲劳损伤交互途中实时损伤累积评估的结果，其中，黑点表示选取的某个循环周次n的累积损伤(Df(n),Dc(n))的位置。已知实验失效时的寿命为255周次，在放大图中可以看出，公式(4)中线性损伤累积法则计算失效时的寿命为236周次，公式(5)中双线性损伤累积法则计算失效时的寿命为226周次，公式(6)中简化的非线性损伤累积法则计算失效时的寿命为197周次。结果显示，预估的寿命较为准确并且均略微保守。
实施例2：
选取数据为总应变范围为1.6％，在应变最大处的保载时间为1800s，温度为650℃下的轴向等幅蠕变疲劳交互试验，参照标准为GB/T 15248-2008《金属材料轴向等幅低循环疲劳试验方法》。其蠕变疲劳实验寿命为353周次，实时损伤累积评估所用的参数包括：失效应变能密度wf、非弹性应变能密度耗散率蠕变试验中施加的应力水平σ、时间t、疲劳寿命参数N0、塑性应变范围参数Δεp，各参数的数值请参阅图4、图5和图6。
实验中可得塑性应变范围参数Δεp、最大应力σ0以及平均应力σm随着循环周次n的变化，根据公式(1)依次计算每个周次的蠕变损伤dc，根据公式(2)依次计算每个周次的疲劳损伤df，然后根据公式(3)计算循环至某个周次n的累积损伤，其中，第n周次时的累积蠕变损伤Dc(n)表示为在损伤交互图中某个周次n的纵坐标，第n周次时的累积疲劳损伤Df(n)表示为损伤交互图中的横坐标。每个循环周次坐标点在蠕变疲劳损伤交互图中所连成轨迹线则为该材料在该工况下实时的累积损伤，当轨迹与损伤累积法则线相交时，则判定为在该工况下材料发生失效，即不能进行实际工况的服役。依据美国ASME和英国R5准则，这里给出的公式(5)中双线性累积损伤法则的折点坐标为(0.3,0.3)，公式(6)中简化的非线性损伤累积法则的系数n＝0.576。
图8表示了本实例中在蠕变疲劳损伤交互途中实时损伤累积评估的结果，其中，黑点表示选取的某个循环周次n的累积损伤(Df(n),Dc(n))的位置。已知实验失效时的寿命为353周次，在放大图中可以看出，在实际材料失效时并未达到公式(4)中的线性损伤累积法则的实现，这表示线性损伤累积法则并不保守，但若是按照累积损伤轨迹线的趋势来看，失效时的寿命约为411周次，公式(5)中双线性损伤累积法则计算失效时的寿命为296周次，公式(6)中简化的非线性损伤累积法则与公式(5)的计算结果非常接近，为295周次。结果显示，预估的寿命较为准确，但是运用线性累积损伤法则会使得结果略微保守，而双线性累积损伤法则和简化的非线性损伤累积法则依然使得结果略微保守。
实施例3：
选取数据为总应变范围为1.0％，在应变最大处的保载时间为300s，温度为650℃下的轴向等幅蠕变疲劳交互试验。其蠕变疲劳实验寿命为2862周次，实时损伤累积评估所用的参数包括：失效应变能密度wf、非弹性应变能密度耗散率蠕变试验中施加的应力水平σ、时间t、疲劳寿命参数N0、塑性应变范围参数Δεp，各参数的数值请参阅图4、图5和图6。
实验中可得塑性应变范围Δεp、最大应力σ0以及平均应力σm随着循环周次n的变化，根据公式(1)依次计算每个周次的蠕变损伤dc，根据公式(2)依次计算每个周次的疲劳损伤df，然后根据公式(3)计算循环至某个周次n的累积损伤，其中，第n周次时的累积蠕变损伤Dc(n)表示为在损伤交互图中某个周次n的纵坐标，第n周次时的累积疲劳损伤Df(n)表示为损伤交互图中的横坐标。每个循环周次坐标点在蠕变疲劳损伤交互图中所连成轨迹线则为该材料在该工况下实时的累积损伤，当轨迹线与损伤累积法则线相交时，则判定为在该工况下材料发生失效，即不能进行实际工况的服役。依据美国ASME和英国R5准则，这里给出的公式(5)中双线性累积损伤法则的折点推荐坐标值为(0.3,0.3)，公式(6)中简化的非线性损伤累积法则的系数n＝0.576。
图9表示了本实例中在蠕变疲劳损伤交互途中实时损伤累积评估的结果，其中，黑点表示选取的某个循环周次n的累积损伤(Df(n),Dc(n))的位置。由于在小应变范围的工况，每个循环贡献的损伤数值较小，因此在30周及以下的累积损伤尚未达到10-4，超出了图9的横纵坐标范围，故这里标出的累积损伤从第40周开始。已知实验失效时的寿命为2862周次，在放大图中可以看出，公式(4)中线性损伤累积法则计算失效时的寿命为2661周次，公式(5)中双线性损伤累积法则计算失效时的寿命为2037周次，公式(6)中简化的非线性损伤累积法则计算失效时的寿命为1815周次。结果显示，预估的寿命较为准确并且均略微保守。
以上记载的，仅为本发明的较佳实施例，并非用以限定本发明的范围，本发明的上述实施例还可以做出各种变化。即凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰，皆落入本发明专利的权利要求保护范围。
Claims (5)
1.一种材料时间相关的蠕变疲劳损伤评定方法，包括步骤：
S1：对一目标材料进行多次实验获取所需的一参数组；
S2：根据所述参数组计算获得各周次的蠕变损伤dc和各周次的疲劳损伤df；其中，所述S2步骤包括利用公式(1)和公式(2)计算获得各周次的蠕变损伤dc和各周次的疲劳损伤df：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<msub>
<mi>t</mi>
<mi>h</mi>
</msub>
</msubsup>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mfrac>
<mi>M</mi>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mi>M</mi>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>log</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>w</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
</msubsup>
<mo>,</mo>
<msub>
<mi>w</mi>
<mrow>
<mi>f</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mfrac>
<mi>M</mi>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mi>M</mi>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>log</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>w</mi>
<mrow>
<mi>f</mi>
<mn>0</mn>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>log&Delta;&epsiv;</mi>
<mi>p</mi>
</msub>
<mo>+</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>E</mi>
<mo>&CenterDot;</mo>
<mi>ln</mi>
<mn>10</mn>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>N</mi>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>log&Delta;&epsiv;</mi>
<mi>p</mi>
</msub>
<mo>+</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mi>E</mi>
<mo>&CenterDot;</mo>
<mi>ln</mi>
<mn>10</mn>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，dc为蠕变损伤，t为时间，D为第一线性回归参数，n1为第一线性回归幂指数，wf0为平台值，为非弹性应变能密度耗散率，σ0为最大应力，σm为平均应力，Δεp为塑性应变范围参数，A为第一松弛参数，B为第二松弛参数；E为弹性模量；th为蠕变疲劳保持时间；
<mrow>
<msub>
<mi>d</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;&epsiv;</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>&beta;</mi>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，df为疲劳损伤，α为第二线性回归参数，β为第二线性回归幂指数；
S3：根据所述各周次的蠕变损伤dc计算第n周次时的累积蠕变损伤Dc(n)和第n周次时的累积疲劳损伤Df(n)，n为大于0的自然数；其中，在所述S3步骤中，根据公式(3)计算所述第n周次时的累积蠕变损伤Dc(n)和所述第n周次时的累积疲劳损伤Df(n)：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>d</mi>
<mrow>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>d</mi>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，dc(i)为第i周次的蠕变损伤，df(i)为第i周次的疲劳损伤，i为大于零的自然数；
S4：以所述第n周次时的累积蠕变损伤Dc(n)为纵坐标并以第n周次时的累积疲劳损伤Df(n)为横坐标绘制一蠕变疲劳损伤交互图；
S5：结合所述蠕变疲劳损伤交互图和至少一损伤累积法则评估所述目标材料的实时累积损伤情况；其中，所述损伤累积法则包括一线性损伤累积法则，所述线性损伤累积法则的数学表达式为公式(4)：
Dc+Df＝1 (4)；
其中，Dc为累积蠕变损伤，Df为累积疲劳损伤，当公式(4)成立时判定所述目标材料失效；所述损伤累积法则还包括一双线性损伤累积法则，所述双线性损伤累积法则的数学表达式为公式(5)：
其中，所述双线性损伤累积法则的折点坐标为当公式(5)成立时判定所述目标材料失效。
2.根据权利要求1所述的时间相关的蠕变疲劳损伤评定方法，其特征在于，所述损伤累积法则还包括一简化的非线性累积损伤法则，所述简化的非线性累积损伤法则的数学表达式为公式(6)：
Dc m+Df m＝1 (6)；
其中，m为简化的非线性损伤累积法则系数，当公式(6)成立时判定所述目标材料失效。
3.根据权利要求1-2任一项所述的时间相关的蠕变疲劳损伤评定方法，其特征在于，所述S1步骤进一步包括步骤：
S11：通过一不同应力水平下的高温蠕变试验，获得所述目标材料的每个应力水平下的一蠕变延性数据εf、一蠕变断裂时间tR、所述第一线性回归参数D、所述第一线性回归幂指数n1和所述平台值wf0；
S12：通过一相同温度不同应力水平的应力松弛试验，获得每个应力水平下应力与松弛时间的变化关系，并根据所述变化关系获得所述第一松弛参数A和所述第二松弛参数B；
S13：通过一相同温度不同应变幅的应变控制下的纯疲劳试验，获得各应变幅下的所述塑性应变范围参数Δεp和疲劳寿命参数N0，根据纯疲劳实验条件下的所述塑性应变范围参数Δεp和疲劳寿命参数N0获得所述第二线性回归参数α和第二线性回归幂指数β；
S14：通过一相同温度的应变控制下的蠕变疲劳试验，记录各周次的应力应变关系，获得各周次的所述最大应力σ0、各周次的所述平均应力σm以及各周次的蠕变疲劳实验条件下的所述塑性应变范围参数Δεp。
4.根据权利要求3所述的时间相关的蠕变疲劳损伤评定方法，其特征在于，所述步骤S12中，根据所述变化关系和公式(7)计算获得所述第一松弛参数A和所述第二松弛参数B：
σ(t)＝σ0-(A·logΔεp+B)·log(1+t) (7)；
其中，t为时间，σ(t)为t时的应力水平。
5.根据权利要求3所述的时间相关的蠕变疲劳损伤评定方法，其特征在于，所述步骤S13中，根据所述塑性应变范围参数Δεp和疲劳寿命参数N0并通过公式(8)计算获得所述第二线性回归参数α和第二线性回归幂指数β：
N0＝α·(Δεp)β (8)；
其中，N0为疲劳寿命参数，Δεp为塑性应变范围参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610531611.8A CN106202913B (zh) | 2016-07-07 | 2016-07-07 | 时间相关的蠕变疲劳损伤评定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610531611.8A CN106202913B (zh) | 2016-07-07 | 2016-07-07 | 时间相关的蠕变疲劳损伤评定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106202913A CN106202913A (zh) | 2016-12-07 |
| CN106202913B true CN106202913B (zh) | 2018-03-06 |
Family
ID=57472511
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610531611.8A Active CN106202913B (zh) | 2016-07-07 | 2016-07-07 | 时间相关的蠕变疲劳损伤评定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106202913B (zh) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106840877B (zh) * | 2017-01-22 | 2019-02-12 | 北京工业大学 | 一种基于应力的多轴小裂纹全寿命预测方法 |
| CN106896133B (zh) * | 2017-02-14 | 2019-04-05 | 北京工业大学 | 一种基于等温疲劳和蠕变疲劳的多轴热机械疲劳寿命预测方法 |
| CN106934168B (zh) * | 2017-03-21 | 2018-08-28 | 中国石油大学(华东) | 一种材料多轴蠕变失效应变预测方法 |
| CN107103140A (zh) * | 2017-04-28 | 2017-08-29 | 电子科技大学 | 一种基于双线性累积损伤的时变疲劳可靠性分析方法 |
| CN107391903A (zh) * | 2017-06-07 | 2017-11-24 | 上海电力学院 | 用于马氏体钢的蠕变‑疲劳寿命预测及其可靠性分析方法 |
| CN108052698B (zh) * | 2017-11-21 | 2020-12-25 | 中国石油大学(华东) | 一种滞回能密度预测材料蠕变疲劳交互寿命的方法 |
| CN108170905B (zh) * | 2017-12-08 | 2020-02-04 | 南昌航空大学 | 一种用于镍基高温合金叶片热机械疲劳载荷下的寿命预测方法 |
| CN108182327B (zh) * | 2017-12-30 | 2021-06-11 | 北京工业大学 | 一种基于线性损伤累积的多轴热机械疲劳寿命预测方法 |
| CN109885874B (zh) * | 2019-01-11 | 2022-12-23 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
| CN110261469A (zh) * | 2019-04-25 | 2019-09-20 | 北京工业大学 | 一种用于铁磁性材料疲劳损伤检测的电磁混频检测方法 |
| CN110208086B (zh) * | 2019-06-05 | 2021-07-09 | 东北石油大学 | 脉动压裂岩石累积损伤求解方法 |
| CN113495032B (zh) * | 2020-03-19 | 2022-11-01 | 上海电气电站设备有限公司 | 材料蠕变疲劳寿命评估方法及评估系统 |
| CN111879636A (zh) * | 2020-07-31 | 2020-11-03 | 华东理工大学 | 一种材料的蠕变-疲劳-氧化实时损伤累积评估方法 |
| CN112651164A (zh) * | 2021-01-20 | 2021-04-13 | 华东理工大学 | 一种基于机器学习的蠕变疲劳寿命预测方法 |
| CN114065576B (zh) * | 2021-11-05 | 2024-03-08 | 华东理工大学 | 一种蠕变疲劳载荷下缺口部件的缺口效应评价方法 |
| CN116052813B (zh) * | 2022-12-30 | 2023-10-13 | 南京工业大学 | 一种全数值驱动的实时蠕变疲劳损伤监测方法 |
| CN115931532B (zh) * | 2022-12-30 | 2023-10-13 | 南京工业大学 | 基于修正时间分数法的电力高温部件损伤状态判定方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102466597A (zh) * | 2010-11-05 | 2012-05-23 | 华东理工大学 | 一种金属构件/材料剩余寿命的无损检测和评价方法 |
| CN103761365A (zh) * | 2013-12-28 | 2014-04-30 | 合肥通用机械研究院 | 一种基于寿命的高温压力容器蠕变疲劳强度设计方法 |
| CN105158084A (zh) * | 2015-09-15 | 2015-12-16 | 华东理工大学 | 一种材料的蠕变-疲劳寿命预测方法 |
-
2016
- 2016-07-07 CN CN201610531611.8A patent/CN106202913B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102466597A (zh) * | 2010-11-05 | 2012-05-23 | 华东理工大学 | 一种金属构件/材料剩余寿命的无损检测和评价方法 |
| CN103761365A (zh) * | 2013-12-28 | 2014-04-30 | 合肥通用机械研究院 | 一种基于寿命的高温压力容器蠕变疲劳强度设计方法 |
| CN105158084A (zh) * | 2015-09-15 | 2015-12-16 | 华东理工大学 | 一种材料的蠕变-疲劳寿命预测方法 |
Non-Patent Citations (1)
| Title |
|---|
| P91钢蠕变-疲劳交互作用研究现状;纪东梅等;《压力容器》;20111231;第28卷(第6期);第39页右栏第3段至42页左栏第2段，图3 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106202913A (zh) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106202913B (zh) | 时间相关的蠕变疲劳损伤评定方法 | |
| Yu et al. | Visco-plastic constitutive modeling on Ohno–Wang kinematic hardening rule for uniaxial ratcheting behavior of Z2CND18. 12N steel | |
| CN105158084A (zh) | 一种材料的蠕变-疲劳寿命预测方法 | |
| Xing et al. | Cyclic deformation of 316L stainless steel and constitutive modeling under non-proportional variable loading path | |
| CN103712865A (zh) | 高温部件的蠕变损伤评价方法以及损伤评价系统 | |
| Nath et al. | Evaluation of ratcheting behaviour in cyclically stable steels through use of a combined kinematic-isotropic hardening rule and a genetic algorithm optimization technique | |
| Xu et al. | Characterizing crack growth behavior and damage evolution in P92 steel under creep-fatigue conditions | |
| Zhang et al. | Fatigue life prediction of FV520B with internal inclusions | |
| Zhao et al. | Modelling creep-fatigue behaviours using a modified combined kinematic and isotropic hardening model considering the damage accumulation | |
| Kim et al. | Probabilistic assessment of creep crack growth rate for Gr. 91 steel | |
| CN107391903A (zh) | 用于马氏体钢的蠕变‑疲劳寿命预测及其可靠性分析方法 | |
| Beesley et al. | A novel simulation for the design of a low cycle fatigue experimental testing programme | |
| Khutia et al. | Development of new cyclic plasticity model for 304LN stainless steel through simulation and experimental investigation | |
| Wang et al. | Creep-fatigue life prediction in nickel-based superalloy GH4169 based on microstructural damage quantification with the help of electron backscatter diffraction | |
| Hossain et al. | Reliability prediction of 304 stainless steel using sine-hyperbolic creep-damage model with Monte Carlo simulation method | |
| Chen et al. | Integrity assessment of a 3D tubeplate using the linear matching method. Part 1. Shakedown, reverse plasticity and ratchetting | |
| Karvan et al. | Time-dependent ratcheting assessment of steel samples undergoing multi-step loading cycles through use of hardening rules | |
| Shimada et al. | Creep-fatigue life prediction of 316H stainless steel by utilizing non-unified constitutive model | |
| JP6582753B2 (ja) | 耐熱鋼材の寿命予測方法 | |
| Jin et al. | Mixed-mode I&II fatigue crack growth behaviors of 16MND5 steel: The role of crack driving forces and crack closure | |
| Hares et al. | The effect of creep strain rate on damage accumulation in type 316H austenitic stainless steel | |
| Shahraini et al. | Ratcheting Analysis of Steel Plate under Cycling Loading using Dynamic Relaxation Method Experimentally Validated | |
| Abe et al. | Statistical model of micro crack growth for the evaluation of accumulated fatigue in NPPs | |
| Belodedenko et al. | Fatigue lifetime model under a complex loading with application of the amalgamating safety indices rule | |
| Hou et al. | Notch Effect on Creep-Fatigue Behavior of Alloy 617 at Elevated Temperature |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address |
Address after: 200237 No. 130, Meilong Road, Shanghai, Xuhui District Patentee after: EAST CHINA University OF SCIENCE AND TECHNOLOGY Country or region after: China Patentee after: AECC COMMERCIAL AIRCRAFT ENGINE Co.,Ltd. Address before: 200237 No. 130, Meilong Road, Shanghai, Xuhui District Patentee before: EAST CHINA University OF SCIENCE AND TECHNOLOGY Country or region before: China Patentee before: AVIC COMMERCIAL AIRCRAFT ENGINE Co.,Ltd. |
|
| CP03 | Change of name, title or address |