[发明专利]一种基于晶体塑性理论的蠕变疲劳剩余寿命评定方法

权利要求书

1.一种基于晶体塑性理论的蠕变疲劳剩余寿命评定方法，其特征在于，包括：

步骤S1，通过ABAQUS有限元软件的单元用户子程序UEL编译基于位错密度的晶体塑性本构方程；

步骤S2，拟合蠕变疲劳迟滞回线以确定所述晶体塑性本构方程的材料参数；

步骤S3，根据所述晶体塑性本构方程以及所述晶体塑性本构方程的材料参数，对待测材料代表性体积单元进行晶体塑性有限元模拟，获得待测材料在不同蠕变疲劳工况下的力学响应；

步骤S4，提取待测材料代表性体积单元在每个循环周次下的累积能量耗散，确定待测材料的局部能量耗散体积占比以及待测材料的容许蠕变损伤和容许疲劳损伤；

步骤S5，以所述容许蠕变损伤和所述容许疲劳损伤为基本坐标轴，以所述局部能量耗散体积占比为第三轴，绘制三维蠕变疲劳剩余寿命评定图，计算待测材料的蠕变疲劳剩余寿命。

2.根据权利要求1所述的基于晶体塑性理论的蠕变疲劳剩余寿命评定方法，其特征在于，所述步骤S1包括：

步骤S11，建立晶体塑性本构方程中的流动准则方程如下：

式中，为参考塑性滑移率，k为玻尔兹曼常数，θ为绝对温度，F₀为热激活自由能，τ^α第α滑移系的分解剪切应力，表示为σ为应力值，s^α表示α滑移系的滑移面法向，m^α表示α滑移系的滑移方向，表示张量的并矢积，B^α为第α滑移系的背应力，S^α为第α滑移系的滑移阻力，p和q为材料常数；

步骤S12，根据所述流动准则方程，在所述晶体塑性本构方程中引入静态回复项的背应力方程如下：

式中，h_B为背应力硬化常数，r_D为滑移阻力相关动态回复系数，r_S为背应力静态回复系数；

步骤S13，建立晶体塑性本构方程中的滑移阻力演化方程如下：

式中，表示滑移阻力变化率，μ是剪切模量，b是柏氏矢量，h^αβ表明滑移系α和β之间的硬化矩阵，和分别为滑移系β上的统计存储位错的刃型部分、螺型部分和几何必须位错部分，K_e表示晶体在α滑移系上刃型位错产生的系数，K_s表示晶体在α滑移系上螺型位错产生的系数，d_e和d_s分别表示刃型位错和螺型位错的湮灭系数，是滑移方向m^α上的纯螺型密度分量，和分别是滑移方向t^α和n^α上的纯刃型密度分量；

步骤S14，建立无位错势垒的临界剪应力方程如下：

式中，为无位错势垒时的临界剪应力，c_τ材料常数，n表示滑移系个数。

3.根据权利要求1所述的基于晶体塑性理论的蠕变疲劳剩余寿命评定方法，其特征在于，所述步骤S2中拟合蠕变疲劳迟滞回线的方法为试参法。

4.根据权利要求1所述的基于晶体塑性理论的蠕变疲劳剩余寿命评定方法，其特征在于，所述步骤S3包括：通过ABAQUS软件计算待测材料代表性体积单元的每个积分点的应力应变值，并将每个积分点的应力应变值平均化，获取待测材料在宏观尺度上的力学响应。

5.根据权利要求1所述的基于晶体塑性理论的蠕变疲劳剩余寿命评定方法，其特征在于，所述步骤S4中确定待测材料的局部能量耗散体积占比的方法包括：

步骤S41，根据所述步骤S3中待测材料代表性体积单元的晶体塑性有限元模拟结果，提取待测材料在蠕变疲劳保载时间前后的累积能量耗散；

步骤S42，根据所述步骤S41中的累积能量耗散，计算待测材料的局部能量耗散体积占比。