CN112364535B - 一种基于晶体塑性的蠕变疲劳寿命预测方法 - Google Patents
一种基于晶体塑性的蠕变疲劳寿命预测方法 Download PDFInfo
- Publication number
- CN112364535B CN112364535B CN202011227885.0A CN202011227885A CN112364535B CN 112364535 B CN112364535 B CN 112364535B CN 202011227885 A CN202011227885 A CN 202011227885A CN 112364535 B CN112364535 B CN 112364535B
- Authority
- CN
- China
- Prior art keywords
- creep
- fatigue
- slip
- creep fatigue
- equation
- 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
- 239000013078 crystal Substances 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000000463 material Substances 0.000 claims abstract description 56
- 229920003023 plastic Polymers 0.000 claims abstract description 53
- 230000000977 initiatory effect Effects 0.000 claims abstract description 39
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims abstract description 37
- 102100021503 ATP-binding cassette sub-family B member 6 Human genes 0.000 claims abstract description 11
- 101100000375 Homo sapiens ABCB6 gene Proteins 0.000 claims abstract description 11
- 230000008859 change Effects 0.000 claims abstract description 7
- 230000010354 integration Effects 0.000 claims abstract description 5
- 238000012935 Averaging Methods 0.000 claims abstract description 4
- 238000010998 test method Methods 0.000 claims abstract description 3
- 230000003068 static effect Effects 0.000 claims description 17
- 230000001186 cumulative effect Effects 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 12
- 238000011084 recovery Methods 0.000 claims description 12
- 230000005489 elastic deformation Effects 0.000 claims description 9
- 238000009661 fatigue test Methods 0.000 claims description 8
- 125000004122 cyclic group Chemical group 0.000 claims description 7
- 230000000737 periodic effect Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012805 post-processing Methods 0.000 claims description 5
- 108010014173 Factor X Proteins 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 230000003213 activating effect Effects 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000014759 maintenance of location Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 229920001887 crystalline plastic Polymers 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002003 electron diffraction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 229910000601 superalloy Inorganic materials 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000007725 thermal activation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种基于晶体塑性理论的蠕变疲劳寿命预测方法,包括:基于电子背散射衍射技术建立ABAQUS的代表性单元模型;将背应力模型修正并写入子程序UMAT,以得到蠕变疲劳迟滞回线;通过试参法拟合蠕变疲劳迟滞回线,来获取材料参数;计算每个积分点的应力应变值并将其平均化,获取蠕变疲劳迟滞回线和后处理云图;从中提取最大的塑性滑移和能量耗散,分析它们随循环周次的变化规律,提出蠕变和疲劳指示因子;根据指示因子预测蠕变疲劳裂纹萌生寿命。本发明的蠕变疲劳寿命预测方法利用塑性滑移和能量耗散作为疲劳和蠕变指示因子,能更好地反映蠕变疲劳损伤演化规律,准确预测裂纹萌生位置,具有直观、适用性强、精确度高的优点。
Description
技术领域
本发明涉及蠕变疲劳的寿命预测领域,尤其涉及一种基于晶体塑性理论的蠕变疲劳寿命预测方法。
背景技术
现代航空发动机、工业燃气机、超临界发电机等装置的工作环境日趋复杂,其中的关键限寿部件除了承受装置稳态运行的恒定荷载之外,通常还承受装置起停和温度波动等引起的交变载荷的作用,其服役过程伴随着严重的蠕变疲劳载荷交互作用,降低了这些关键限寿部件的服役寿命。因此,为了保证关键限寿部件的安全性和可靠性,构建准确的蠕变疲劳寿命预测模型是非常重要的。
为了准确地评估高温结构在蠕变疲劳载荷下的耐久性问题,自20世纪50年代以来,国内外学者已先后提出了几十种蠕变疲劳寿命预测模型。近些年,有限元软件的发展为了更好地理解复杂工况下的应力应变行为以及实现该工况下的寿命预测。ABAQUS软件允许用户通过二次开发接口补充ABAQUS蠕变前后处理模块中无法实现的功能。其中,基于FORTRAN语言的材料用户自定义的子程序UMAT扩展了ABAQUS在本构方程方面上的应用。用户可以基于UMAT根据自己的材料以及需求,来编译所需的本构模型并通过ABAQUS软件实现模拟。
现今针对蠕变疲劳分析及寿命预测的主要能够描述包含裂纹萌生和扩展阶段的连续损伤力学理论,其通过损伤变量引入统一粘塑性理论的方式描述材料在循环载荷下损伤累积直至断裂的过程。这类方法侧重于宏观层面的类构件或者构件的损伤分析,无法描述微观层面的应力应变响应以及损伤机理演化。由于蠕变疲劳裂纹萌生通常是在微观层面上,使用这类方法很好精确地判定裂纹萌生以及揭示损伤机理。因此,如何有效地评估裂纹萌生寿命以及分析蠕变疲劳损伤机理,对关键限寿部件是非常重要的。
发明内容
针对上述现有技术中的不足,本发明提供一种基于晶体塑性的蠕变疲劳寿命预测方法,能够更好地实现蠕变疲劳损伤分析以及蠕变疲劳剩余寿命预测,具有直观、适用性强、精确度高的优点。
为了实现上述目的,本发明提供一种基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,包括以下步骤:
S1:建立基于材料的电子背散射衍射信息的ABAQUS的代表性单元模型;
S2:将背应力模型修正为含有静态回复项的背应力模型并写入晶体塑性本构方程,该晶体塑性本构方程通过用户自定义的子程序UMAT定义,将该晶体塑性本构方程用于对所述S1中的ABAQUS代表性单元模型进行有限元计算,以此模拟得到所述材料在给定载荷条件下的蠕变疲劳迟滞回线;
S3:通过试参法拟合蠕变疲劳迟滞回线,来获取所述晶体塑性本构方程的材料参数;
S4:通过ABAQUS软件计算每个积分点的应力应变值,并将应力应变值平均化,获取材料的蠕变疲劳迟滞回线和后处理云图;
S5:从后处理云图中提取最大的塑性滑移和能量耗散,将这两个参数视为材料在循环加载过程中损伤驱动力,分析它们随着循环周次的变化规律,从而提出一种蠕变指示因子和疲劳指示因子;
S6:结合所述S5的蠕变指示因子和疲劳指示因子和线性累积损伤准则其中,/>和/>分别表示第j周次的蠕变损伤和疲劳损伤,预测不同保载时间下的蠕变疲劳裂纹萌生寿命。
进一步地,所述S1包括:
S11:通过matlab语言将材料的电子背散射衍射信息转化为ABAQUS可读取的input文件中的节点和单元以及相应的晶粒取向文件;
S12:通过matlab语言读取所述S11中的晶粒取向文件,从而编译每个晶粒的多个材料常数;
S13:通过编写ABAQUS的脚本程序获得代表性体积单元的周期性边界条件和蠕变疲劳加载条件;
S14:将S11中的节点和单元、S12中的材料属性、S13中的周期性边界条件和蠕变疲劳加载条件手动组成ABAQUS可读取的input文件,从而获得ABAQUS代表性体积单元模型。
优选地,在所述S11中,所述晶粒的取向参数采用三个欧拉角φ、/>表示。
进一步地,在所述S2中,晶体塑性本构方程包括主控方程,流动准则方程,滑移阻力演化方程以及引入静态回复项的背应力方程。。
进一步地,所述S2包括:
S21:建立描述蠕变疲劳行为的晶体塑性本构方程中变形速率梯度Lp的主控方程,所述的变形速率梯度Lp的主控方程为:
其中,Lp为非弹性变形速率梯度,Fe为弹性变形梯度,Fp为非弹性变形梯度,为第α滑移系的塑性滑移速率,sα为第α滑移系滑移方向的单位矢量,mα为第α滑移系法向方向的单位矢量;
S22:建立晶体塑性本构方程中的流动准则方程,流动准则方程为:
其中,为参考塑性滑移速率,k为玻尔兹曼常数,T为绝对温度,F0为热激活自由能,τα第α滑移系的分解剪切应力,表示为/>σ为应力值,Bα为第α滑移系的背应力,Sα为第α滑移系的滑移阻力,τ0、p、q为材料常数;
S23:建立晶体塑性本构方程中的滑移阻力演化方程,滑移阻力演化方程为:
hαβ=hs[w+(1-w)δαβ],
其中,hαβ表明滑移系α和β之间的硬化矩阵,Ssat为饱和滑移阻力,S0为初始滑移阻力,hs为静态硬化模量,w代表晶格硬化率,δαβ代表克罗内克函数;
S24:建立晶体塑性本构方程中引入静态回复项的背应力方程,引入静态回复项的背应力方程为:
其中,hB为背应力硬化常数,rD为滑移阻力相关动态回复系数,rS为背应力静态回复系数。
进一步地,在所述S3中,通过试参法拟合蠕变疲劳迟滞回线是通过将S2得到的材料在给定载荷条件下的蠕变疲劳迟滞回线与该材料的蠕变疲劳试验的数据对比来实现的。
进一步地,在在所述S4中,结合所述S1中建立的ABAQUS的代表性单元模型、S2中构建的晶体塑性本构方程、S3中确定的晶体塑性本构方程的材料参数,并通过ABAQUS软件来计算得到每个积分点的应力应变值。
进一步地,所述步骤S5包括:
S51:从ABAQUS后处理云图中提取蠕变疲劳保载时间前后的最大塑性滑移和能量耗散,分析保载时间前后的塑性滑移和能量耗散的变化规律;
S52:将所述步骤S51中的塑性滑移和/或能量耗散视为蠕变指示因子和疲劳指示因子。
进一步地,所述步骤S6中,预测不同保载时间下的蠕变疲劳裂纹萌生寿命包括:构建蠕变疲劳裂纹萌生寿命预测模型和确定蠕变疲劳裂纹萌生寿命。
进一步地,所述S6包括:
S61:构建蠕变疲劳裂纹萌生寿命预测模型,所述蠕变疲劳裂纹萌生寿命预测模型为:
其中,为第j循环的蠕变损伤,/>为第j循环的疲劳损伤,Ni为蠕变疲劳裂纹萌生寿命,/>是第j循环的蠕变指示因子,/>是第j循环的疲劳指示因子,Xc,crit是蠕变指示因子的临界值,Xf,crit是疲劳指示因子的临界值,X表示塑性滑移p和/或能量耗散W;
塑性滑移p和能量耗散W的公式为:
其中,为塑性滑移速率,Lp为非弹性变形速率梯度,τα为第α滑移系临界剪切应力,/>为第α滑移系的塑性滑移速率,n为滑移系个数;
S62:确定蠕变疲劳裂纹萌生寿命,所述蠕变疲劳裂纹萌生寿命为:
其中,Xc,cycle是任意循环周次下的蠕变指示因子,Xf,cycle是任意循环周次下的疲劳指示因子。
进一步地,所述基于晶体塑性的蠕变疲劳寿命预测方法还包括步骤S7:基于累积塑性滑移占比构建蠕变疲劳剩余寿命评定图。
本发明的基于晶体塑性的蠕变疲劳寿命预测方法利用塑性滑移和能量耗散作为疲劳指示因子和蠕变指示因子,可以将微观尺度的参数用做损伤参量来进行寿命的评估,这种方法的好处在于微观尺度上的预测,可以有效地获得蠕变疲劳裂纹萌生位置和精确地预测蠕变疲劳裂纹萌生寿命。本发明不仅可以预测寿命,还能够获得每个循环周次中的塑性滑移损伤积累图,由此预测裂纹萌生的位置。
本发明的基于晶体塑性的蠕变疲劳寿命预测方法能够得到材料在蠕变疲劳载荷作用下的应力应变行为以及损伤演化规律,其中晶体塑性理论耦合了修正的背应力模型并写入材料用户自定义的子程序UMAT。
附图说明
图1为根据本发明的一个实施例的基于晶体塑性的蠕变疲劳寿命预测方法的流程图;
图2为根据本发明的一个实施例的基于晶体塑性的蠕变疲劳寿命预测方法的ABAQUS代表性体积单元模型图;
图3为根据本发明的一个实施例的基于晶体塑性的蠕变疲劳寿命预测方法的有限元模型的周期性边界条件和加载条件;
图4为根据本发明的一个实施例的基于晶体塑性的蠕变疲劳寿命预测方法的蠕变疲劳模拟曲线的拟合结果图;
图5为根据本发明的一个实施例的基于晶体塑性的蠕变疲劳寿命预测方法的累积非弹性滑移和能量耗散随着循环周次线性增加的轨迹图;
图6给出了累积非弹性滑移和能量耗散随着蠕变疲劳循环周次变化的轨迹图,其中,图6(a)示出了累积非弹性滑移随着蠕变疲劳循环周次变化的轨迹,图6(b)示出了能量耗散随着蠕变疲劳循环周次变化的轨迹;
图7为根据本发明的一个实施例的基于晶体塑性的蠕变疲劳寿命预测方法的试验寿命和预测寿命的对比图,其中,图7(a)示出了基于非弹性滑移的预测结果;图7(b)示出了基于能量耗散的预测结果。
具体实施方式
下面根据附图,给出本发明的较佳实施例,并予以详细描述,使能更好地理解本发明的功能和特点。
请参阅图1为本发明公开的一种基于晶体塑性的蠕变疲劳寿命预测方法,其包括以下步骤:
S1:通过背散射电子衍射技术获得原始材料的EBSD(电子背散射衍射)信息,EBSD信息主要包括晶体取向和晶体结构信息,基于材料的EBSD信息建立ABAQUS代表性单元模型;
其中,S1具体包括:
S11:通过matlab语言将材料的EBSD信息转化为ABAQUS可读取的input文件中的节点和单元以及相应的晶粒取向文件,其中,晶粒的取向参数采用三个欧拉角φ、/>表示,input文件中包含了有限元节点和单元,这里的节点和单元就是通过网格划分生成的节点和单元。
S12:通过matlab语言读取所述步骤S11中的晶粒取向文件,从而编译每个晶粒的多个材料常数;在本实施例中,每个晶粒的材料常数的数量为15个,这些材料常数用于在ABAQUS计算过程中供UMAT调用,用于控制UMAT中背应力、滑移阻力和滑移系等的开启,并用于反映不同晶粒的不同取向;
S13:通过编写ABAQUS的脚本程序获得代表性体积单元的周期性边界条件和蠕变疲劳加载条件;
S14:将S11中的节点和单元、S12中的材料常数、S13中的周期性边界条件和蠕变疲劳加载条件手动组成ABAQUS可读取的input文件,从而获得ABAQUS代表性体积单元模型,用于蠕变疲劳有限元模拟。
S2:将背应力模型修正为含有静态回复项的背应力模型并写入晶体塑性本构方程,晶体塑性本构方程通过ABAQUS有限元软件提供的材料用户自定义子程序UMAT编译,由此,修正后的晶体塑性本构方程包括主控方程,流动准则方程,滑移阻力演化方程以及引入静态回复项的背应力方程,该晶体塑性本构方程用来描述材料的蠕变疲劳迟滞曲线;最后,通过ABAQUS软件调用UMAT子程序,将该晶体塑性本构方程用于对步骤S1中的ABAQUS代表性单元模型进行有限元计算,一方面得到材料在给定载荷条件下的蠕变疲劳迟滞回线,另一方面获得材料的塑性滑移和能量耗散;
其中,步骤S2进一步包括步骤:
S21:建立描述蠕变疲劳行为的晶体塑性本构方程中变形速率梯度Lp的主控方程:
其中,Lp为非弹性变形速率梯度,Fe为弹性变形梯度,Fp为非弹性变形梯度,为第α滑移系的塑性滑移速率(单位为s-1),sα为第α滑移系滑移方向的单位矢量,mα为第α滑移系法向方向的单位矢量;
S22:建立晶体塑性本构方程中的流动准则方程:
其中,为参考塑性滑移速率(单位为s-1),k为玻尔兹曼常数,T为绝对温度(单位为K),F0为热激活自由能(单位为kJ·mol-1),τα第α滑移系的分解剪切应力(单位为MPa),表示为/>σ为应力值(单位为MPa),Bα为第α滑移系的背应力(单位为MPa),Sα为第α滑移系的滑移阻力(单位为MPa),τ0、p、q为材料常数,运算符号<>表示:当x>0时,<x>=x;当x≤0时,<x>0。
S23:建立晶体塑性本构方程中的滑移阻力演化方程:
hαβ=hs[w+(1-w)δαβ] (4);
其中,hαβ表明滑移系α和β之间的硬化矩阵,Ssat为饱和滑移阻力(单位为MPa),S0为初始滑移阻力(单位为MPa),hs为静态硬化模量(单位为MPa),w代表晶格硬化率,δαβ代表克罗内克函数;
S24:由于考虑蠕变回复效应,建立引入静态回复项的背应力方程:
其中,hB为背应力硬化常数(单位为MPa),rD为滑移阻力相关动态回复系数(单位为MPa),rS为背应力静态回复系数(单位为s-1)。
S3:通过将S2得到的材料在给定载荷条件下的蠕变疲劳迟滞回线与该材料的蠕变疲劳试验的数据对比以通过试参法拟合蠕变疲劳迟滞回线,从而通过试参法拟合蠕变疲劳迟滞回线,来获取晶体塑性本构方程所需的材料参数k,F0,τ0,p,q,Ssat,S0,hs,hB,rD,rS。由此,通过与蠕变疲劳试验获得的迟滞回线对比来确定S3中拟合得到的材料参数的合理性。
S4:通过ABAQUS软件计算每个积分点的应力应变值,并将应力应变值平均化,获取材料的蠕变疲劳迟滞回线和后处理云图。由此,得到了材料的宏观的力学响应。
在所述S4中,结合所述S1中建立的ABAQUS的代表性单元模型、S2中构建的晶体塑性本构方程、S3中确定的晶体塑性本构方程的材料参数,并通过ABAQUS软件来计算得到每个积分点的应力应变值。
S5:从后处理云图中提取最大的塑性滑移和能量耗散,将这两个参数视为材料在循环加载过程中损伤驱动力,分析它们随着循环周次的变化规律,从而提出一种蠕变指示因子和疲劳指示因子。
其中,所述塑性滑移和能量耗散作为ABAQUS的状态变量通过后处理输出。
所述步骤S5包括:
S51:从ABAQUS后处理云图中提取蠕变疲劳保载时间前后的最大的塑性滑移和能量耗散,分析保载时间前后的塑性滑移和能量耗散的变化规律;
S52:将所述步骤S51中的塑性滑移和/或能量耗散视为蠕变指示因子和疲劳指示因子(即蠕变疲劳损伤失效的指示因子)。
S6:基于所述S5的蠕变指示因子和疲劳指示因子和线性累积损伤准则,预测不同保载时间下的蠕变疲劳裂纹萌生寿命;
其中,线性累积损伤准则就是把累积疲劳损伤和蠕变损伤进行线性叠加,即其中,/>和/>分别表示第j周次的蠕变损伤和疲劳损伤。预测不同保载时间下的蠕变疲劳裂纹萌生寿命包括构建蠕变疲劳裂纹萌生寿命预测模型和确定蠕变疲劳裂纹萌生寿命。
其中,步骤S6进一步包括:
S61:构建蠕变疲劳裂纹萌生寿命预测模型,所述蠕变疲劳裂纹萌生寿命预测模型为:
其中,为第j循环的蠕变损伤,/>为第j循环的疲劳损伤,Ni为蠕变疲劳裂纹萌生寿命,/>是第j循环的蠕变指示因子,/>是第j循环的疲劳指示因子,Xc,crit是蠕变指示因子的临界值,Xf,crit是疲劳指示因子的临界值,X表示塑性滑移p和/或能量耗散W;
塑性滑移p和能量耗散W的公式为:
为塑性滑移速率,Lp为非弹性变形速率梯度,τα为第α滑移系临界剪切应力,/>为第α滑移系的塑性滑移速率,n为滑移系个数。
塑性滑移p和能量耗散W的计算公式编写在用户子程序UMAT中,通过ABAQUS有限元软件计算,提取蠕变疲劳每个周次下以及保载时间前后的塑性滑移p和能量耗散W,用于获得蠕变指示因子的临界值Xc,crit和疲劳指示因子的临界值Xf,crit。
S62:结合S61所述的蠕变疲劳裂纹萌生寿命预测模型,确定蠕变疲劳裂纹萌生寿命,蠕变疲劳裂纹萌生寿命为:
其中,Xc,cycle是任意循环周次下的蠕变指示因子,Xf,cycle是任意循环周次下的疲劳指示因子。
此外,还可以包括S7:基于累积塑性滑移占比构建蠕变疲劳剩余寿命评定图。
实验结果
下面采用本发明提供的基于晶体塑性的蠕变疲劳寿命预测方法针对镍基高温合金GH4169,并采用实验的裂纹萌生寿命信息来验证本发明的有效性。其中,蠕变疲劳试验在高温650℃下进行,棒状试样的加载方式为应变控制,应变比为-1,加载的应变范围为1%、1.4%、1.6%、2%,保载时间施加在最大拉应变处,保载时间分别为120s、300s、1800s。
本发明的基于晶体塑性的蠕变疲劳寿命预测方法需要对同样材质的试样进行高温下的不同应变幅的单轴蠕变疲劳试验,获得的试验结果用于在本发明的蠕变疲劳寿命预测方法的步骤S3中确定晶体塑性本构方程所需的材料参数。首先基于EBSD信息构建如图2所示的ABAQUS代表性体积单元模型,代表性体积单元模型的约束条件和加载方式如图3所示。通过试参法调整蠕变疲劳试验的模拟结果,使之与蠕变疲劳试验的滞回曲线和应力松弛曲线相吻合,拟合结果如图4和图5所示。其中,拟合结果为:F0=295kJ·mol-1,p=0.96,q=1.12,τ0=702MPa,hs=360MPa,Ssat=317MPa,S0=340MPa,hB=850MPa,rD=8MPa,rS=0.001/s。
通过采用本发明的基于晶体塑性的蠕变疲劳寿命预测方法,可以根据蠕变疲劳损伤的演化规律,确定蠕变指示因子的临界值Xc,crit和疲劳指示因子的临界值Xf,crit,然后提取蠕变疲劳一周次的塑性滑移和能量耗散,根据公式(12)可以确定蠕变疲劳裂纹萌生寿命。
图6给出了累积非弹性滑移(即累积塑性滑移)和能量耗散随着蠕变疲劳循环周次变化的轨迹图,其中,其中图6(a)示出了累积非弹性滑移随着蠕变疲劳循环周次变化的轨迹,图6(b)示出了能量耗散随着蠕变疲劳循环周次变化的轨迹。可以看出,随着保载时间的增加,累积非弹性滑移和能量耗散也相应增加。此外,累积非弹性滑移和能量耗散均与循环周次呈线性关系,因此可以用步骤S6所述公式(12)计算不同工况下的蠕变疲劳萌生寿命。
图7给出了不同蠕变疲劳工况下的裂纹萌生寿命预测图,其中,图7(a)示出了基于非弹性滑移的预测结果;图7(b)示出了基于能量耗散的预测结果。可以看出,本发明的蠕变疲劳寿命预测方法通过提出的蠕变指示因子和疲劳指示因子来预测萌生寿命的结果,与实验的结果较为接近,在2倍误差带范围内,此外,基于能量耗散的寿命预测结果比基于非弹性滑移的预测结果更加精确。因此,证明了此数值模拟方法的具有较高的可靠性。
以上记载的,仅为本发明的较佳实施例,并非用以限定本发明的范围,本发明的上述实施例还可以做出各种变化。即凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰,皆落入本发明专利的权利要求保护范围。
Claims (8)
1.一种基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,包括以下步骤:
S1:建立基于材料的电子背散射衍射信息的ABAQUS的代表性单元模型;
S2:将背应力模型修正为含有静态回复项的背应力模型并写入晶体塑性本构方程,该晶体塑性本构方程通过用户自定义的子程序UMAT定义,将该晶体塑性本构方程用于对所述S1中的ABAQUS代表性单元模型进行有限元计算,以此模拟得到所述材料在给定载荷条件下的蠕变疲劳迟滞回线;
S3:通过试参法拟合蠕变疲劳迟滞回线,来获取所述晶体塑性本构方程的材料参数;
S4:通过ABAQUS软件计算每个积分点的应力应变值,并将应力应变值平均化,获取材料的蠕变疲劳迟滞回线和后处理云图;
S5:从后处理云图中提取最大的塑性滑移和能量耗散,将这两个参数视为材料在循环加载过程中损伤驱动力,分析它们随着循环周次的变化规律,从而提出一种蠕变指示因子和疲劳指示因子;
S6:结合所述S5的蠕变指示因子和疲劳指示因子和线性累积损伤准则其中,/>和/>分别表示第j周次的蠕变损伤和疲劳损伤,预测不同保载时间下的蠕变疲劳裂纹萌生寿命;
所述S2包括:
S21:建立描述蠕变疲劳行为的晶体塑性本构方程中变形速率梯度Lp的主控方程,所述的变形速率梯度Lp的主控方程为:
其中,Lp为非弹性变形速率梯度,Fe为弹性变形梯度,Fp为非弹性变形梯度,为第α滑移系的塑性滑移速率,sα为第α滑移系滑移方向的单位矢量,mα为第α滑移系法向方向的单位矢量;
S22:建立晶体塑性本构方程中的流动准则方程,流动准则方程为:
其中,为参考塑性滑移速率,k为玻尔兹曼常数,T为绝对温度,F0为热激活自由能,τα第α滑移系的分解剪切应力,表示为/>σ为应力值,Bα为第α滑移系的背应力,Sα为第α滑移系的滑移阻力,τ0、p、q为材料常数;
S23:建立晶体塑性本构方程中的滑移阻力演化方程,滑移阻力演化方程为:
hαβ=hs[w+(1-w)δαβ],
其中,hαβ表明滑移系α和β之间的硬化矩阵,Ssat为饱和滑移阻力,S0为初始滑移阻力,hs为静态硬化模量,w代表晶格硬化率,δαβ代表克罗内克函数;
S24:建立晶体塑性本构方程中引入静态回复项的背应力方程,引入静态回复项的背应力方程为:
其中,hB为背应力硬化常数,rD为滑移阻力相关动态回复系数,rS为背应力静态回复系数;
所述S6包括:
S61:构建蠕变疲劳裂纹萌生寿命预测模型,所述蠕变疲劳裂纹萌生寿命预测模型为:
其中,为第j循环的蠕变损伤,/>为第j循环的疲劳损伤,Ni为蠕变疲劳裂纹萌生寿命,/>是第j循环的蠕变指示因子,/>是第j循环的疲劳指示因子,Xc,crit是蠕变指示因子的临界值,Xf,crit是疲劳指示因子的临界值,X表示塑性滑移p和/或能量耗散;
塑性滑移p和能量耗散W的公式为:
其中,为塑性滑移速率,Lp为非弹性变形速率梯度,τα为第α滑移系临界剪切应力,/>为第α滑移系的塑性滑移速率,n为滑移系个数;
S62:确定蠕变疲劳裂纹萌生寿命,所述蠕变疲劳裂纹萌生寿命为:
其中,Xc,cycle是任意循环周次下的蠕变指示因子,Xf,cycle是任意循环周次下的疲劳指示因子。
2.根据权利要求1所述的基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,所述S1包括:
S11:通过matlab语言将材料的电子背散射衍射信息转化为ABAQUS可读取的input文件中的节点和单元以及相应的晶粒取向文件;
S12:通过matlab语言读取所述S11中的晶粒取向文件,从而编译每个晶粒的多个材料常数;
S13:通过编写ABAQUS的脚本程序获得代表性体积单元的周期性边界条件和蠕变疲劳加载条件;
S14:将S11中的节点和单元、S12中的材料属性、S13中的周期性边界条件和蠕变疲劳加载条件手动组成ABAQUS可读取的input文件,从而获得ABAQUS代表性体积单元模型。
3.根据权利要求2所述的基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,在所述S11中,所述晶粒的取向参数采用三个欧拉角φ、/>表示。
4.根据权利要求1所述的基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,在所述S2中,晶体塑性本构方程包括主控方程,流动准则方程,滑移阻力演化方程以及引入静态回复项的背应力方程。
5.根据权利要求1所述的基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,在所述S3中,通过试参法拟合蠕变疲劳迟滞回线是通过将S2得到的材料在给定载荷条件下的蠕变疲劳迟滞回线与该材料的蠕变疲劳试验的数据对比来实现的。
6.根据权利要求1所述的基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,在所述S4中,结合所述S1中建立的ABAQUS的代表性单元模型、S2中构建的晶体塑性本构方程、S3中确定的晶体塑性本构方程的材料参数,并通过ABAQUS软件来计算得到每个积分点的应力应变值。
7.根据权利要求1所述的基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,所述步骤S5包括:
S51:从ABAQUS后处理云图中提取蠕变疲劳保载时间前后的最大塑性滑移和能量耗散,分析保载时间前后的塑性滑移和能量耗散的变化规律;
S52:将所述步骤S51中的塑性滑移和/或能量耗散视为蠕变指示因子和疲劳指示因子。
8.根据权利要求1所述的基于晶体塑性的蠕变疲劳寿命预测方法,其特征在于,在所述S6中,预测不同保载时间下的蠕变疲劳裂纹萌生寿命包括:构建蠕变疲劳裂纹萌生寿命预测模型和确定蠕变疲劳裂纹萌生寿命。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011227885.0A CN112364535B (zh) | 2020-11-06 | 2020-11-06 | 一种基于晶体塑性的蠕变疲劳寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011227885.0A CN112364535B (zh) | 2020-11-06 | 2020-11-06 | 一种基于晶体塑性的蠕变疲劳寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112364535A CN112364535A (zh) | 2021-02-12 |
CN112364535B true CN112364535B (zh) | 2024-03-29 |
Family
ID=74509709
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011227885.0A Active CN112364535B (zh) | 2020-11-06 | 2020-11-06 | 一种基于晶体塑性的蠕变疲劳寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112364535B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113358678B (zh) * | 2021-05-11 | 2022-08-09 | 哈尔滨工业大学(深圳) | α钛变形过程介观应力和织构的半定量预测及可视化方法 |
CN113252479A (zh) * | 2021-05-14 | 2021-08-13 | 北京理工大学 | 一种考虑加工表面完整性的能量法预测疲劳寿命方法 |
CN113591340B (zh) * | 2021-06-24 | 2023-09-01 | 内蒙古工业大学 | 随机载荷下齿轮弯曲疲劳时变可靠性分析方法 |
CN113611377B (zh) * | 2021-07-05 | 2023-08-22 | 南京工业大学 | 一种利用晶体塑性模型模拟混合控制蠕变疲劳变形的方法 |
CN113792446B (zh) * | 2021-07-28 | 2024-03-29 | 华东理工大学 | 一种基于晶体塑性理论的蠕变疲劳剩余寿命评定方法 |
CN114047210B (zh) * | 2021-10-28 | 2023-12-12 | 北京理工大学 | 一种考虑表面完整性的疲劳裂纹萌生预测方法 |
CN113962135B (zh) * | 2021-11-29 | 2024-10-25 | 上海索辰信息科技股份有限公司 | 一种用于航空发动机部件蠕变寿命预测分析方法 |
CN114091315B (zh) * | 2022-01-24 | 2022-04-15 | 西北工业大学 | 高温合金应力及损伤演化方法、装置、存储介质及电子设备 |
CN114662356B (zh) * | 2022-02-09 | 2024-09-20 | 南京航空航天大学 | 一种基于ebsd表征和晶体塑性的疲劳微裂纹扩展预测方法 |
CN116052813B (zh) * | 2022-12-30 | 2023-10-13 | 南京工业大学 | 一种全数值驱动的实时蠕变疲劳损伤监测方法 |
CN116296820B (zh) * | 2023-03-31 | 2024-04-05 | 四川大学 | 双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品 |
CN117973126B (zh) * | 2024-01-29 | 2024-08-27 | 北京科技大学 | 一种基于位错不可逆运动耗散能准则的疲劳寿命预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885874A (zh) * | 2019-01-11 | 2019-06-14 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
CN110826285A (zh) * | 2019-11-22 | 2020-02-21 | 华东理工大学 | 一种几何不连续结构的疲劳寿命预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8515688B2 (en) * | 2009-03-12 | 2013-08-20 | GM Global Technology Operations LLC | Systems and methods to predict fatigue lives of aluminum alloys under multiaxial loading |
-
2020
- 2020-11-06 CN CN202011227885.0A patent/CN112364535B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885874A (zh) * | 2019-01-11 | 2019-06-14 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
CN110826285A (zh) * | 2019-11-22 | 2020-02-21 | 华东理工大学 | 一种几何不连续结构的疲劳寿命预测方法 |
Non-Patent Citations (2)
Title |
---|
卿华 ; 江和甫 ; 温卫东 ; 吴长波 ; 胡仁高 ; 覃志贤 ; 孙景国 ; .DD6单晶合金循环蠕变性能研究.燃气涡轮试验与研究.2007,(03),全文. * |
李聪成 ; 荆洪阳 ; 徐连勇 ; 韩永典 ; 张文 ; .蠕变疲劳交互作用下裂纹萌生的有限元模拟.焊接学报.2016,(08),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN112364535A (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112364535B (zh) | 一种基于晶体塑性的蠕变疲劳寿命预测方法 | |
CN110826285B (zh) | 一种几何不连续结构的疲劳寿命预测方法 | |
Seifert et al. | Mechanism-based thermomechanical fatigue life prediction of cast iron. Part I: Models | |
CN112966347B (zh) | 一种非连续结构双尺度蠕变疲劳寿命预测方法 | |
Ahmed et al. | Constitutive modeling for thermo-mechanical low-cycle fatigue-creep stress–strain responses of Haynes 230 | |
Beesley et al. | A novel simulation for the design of a low cycle fatigue experimental testing programme | |
CN113792446B (zh) | 一种基于晶体塑性理论的蠕变疲劳剩余寿命评定方法 | |
Seifert et al. | Constitutive Modeling and Thermomechanical Fatigue Life Predictions of A356-T6 Aluminum Cylinder Heads Considering Ageing Effects | |
Saltsman et al. | An update of the total-strain version of strainrange partitioning | |
Khoei et al. | Numerical simulation of cyclic behavior of ductile metals with a coupled damage–viscoplasticity model | |
Cheng et al. | Optimal sequential ALT plans for systems with mixture of one-shot units | |
Messner et al. | Development of the Technical Basis of a Unified Viscoplastic Model of 316H Stainless Steel for Incorporation into ASME Division 5 | |
Messner et al. | FY17 Status Report on the Initial EPP Finite Element Analysis of Grade 91 Steel | |
Lindström et al. | Accounting for anisotropic, anisothermal, and inelastic effects in crack initiation lifing of additively manufactured components | |
Evans | A new statistical framework for the determination of safe creep life using the theta projection technique | |
CN117454735A (zh) | 模型生成方法、可靠性分析方法、系统、设备和介质 | |
Messner et al. | Development of Grade 91 inelastic model for incorporation in ASME Division 5 | |
JP2010008396A (ja) | 係数算出装置、係数算出方法、及び係数算出プログラム | |
Cathcart et al. | Probabilistic lifing methods for digital assets | |
O’Nora et al. | Constitutive modeling of TMF and creep-fatigue of a Ni-base alloy | |
Dileep et al. | Effect of fatigue damage parameter on the cumulative life of a turbine rotor under multiaxial loading | |
Nagode et al. | The Damage Operator Approach: Fatigue, Creep and Viscoplasticity Modeling in Thermo-Mechanical Fatigue | |
Messner et al. | A Computer Design Tool for Ceramic Receivers | |
Lytwyn et al. | Ratchet analysis of structures under a generalised cyclic load history | |
Green et al. | A Mean Stress Correction to a Ductility Exhaustion Approach for Gas Turbine Component Life Assessment |
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 |