技术领域
[0001] 本
发明涉及数值模拟领域，尤其涉及一种多轴蠕变疲劳预测方法。
背景技术
[0002] 随着高温旋转部件大型化、高性能、长寿命的需求日益提升，包含几何不连续的结构在复杂苛刻环境下的结构完整性评估成为亟需解决的关键技术
瓶颈之一。由几何不连续和复杂加载历史导致的多轴应
力应变状态不可避免地成为限制了此类部件的使用寿命。近些年，有限元
软件的发展可以很好地满足人们对于复杂
应力应变行为的理解以及提供了该状态下精确寿命预测的可行性。
[0003] 现今针对复杂结构的蠕变疲劳分析及寿命预测的主要可以分为三类。第一类是近年基于微观分析方法的晶体塑性有限元理论，其通过累积塑性滑移带和存储的
能量耗散等疲劳指示因子表征蠕变疲劳演化过程。第二类是能够描述包含裂纹萌生和扩展阶段的连续损伤力学理论，其通过损伤变量引入统一的蠕变疲劳本构的方式描述材料在循环
载荷下损伤累积直至断裂的过程。第三类流行于当前的设计准则，其通过非统一本构模型描述蠕变疲劳行为，并通过分别计算的蠕变损伤和疲劳损伤方程以及唯象的包络线来预测蠕变寿命。但是这三类方法有各自的缺点：第一类方法仅适用于描述微观层面的应力应变行为，对于宏观层面的大型高温部件并不适用；第二类方法侧重于描述裂纹扩展阶段的蠕变疲劳行为，其编程复杂、收敛性差、计算成本高等特点决定了这类方法不具有很强的普适性；第三类方法虽然具有操作性强的特点，但是其多用于分析稳态情况下单轴应力应变行为，对于复杂应力应变状态和复杂加载历史情况下的蠕变疲劳分析和预测并不准确。
[0004] 基于此，期望获得一种新的蠕变疲劳预测方法，以更好地实现几何不连续结构在多轴应力应变状态下的蠕变疲劳分析，获得更为直观、精确度高的结果，且该蠕变疲劳预测方法实用性强。
发明内容
[0005] 本发明的目的在于提供一种基于ABAQUS的多轴蠕变疲劳预测方法，该多轴蠕变疲劳预测方法可以更好地实现几何不连续结构在多轴应力应变状态下的蠕变疲劳分析，获得更为直观、精确度高的结果，且该蠕变疲劳预测方法实用性强。
[0006] 根据上述发明目的，本发明提出了一种基于ABAQUS的多轴蠕变疲劳预测方法，其包括步骤：
[0007] S1：建立ABAQUS有限元模型，通过用户子程序UMAT定义待测试样材料在循环载荷过程中的粘塑性本构方程；
[0008] S2：确定所述粘塑性本构方程所需的模型参数；
[0009] S3：建立待测试材料的多轴应力应变状态的疲劳损伤计算模型和蠕变损伤计算模型；
[0010] S4：建立多轴应力应变状态下的ABAQUS有限元模型，基于步骤S1中用户子程序UMAT所定义的粘塑性本构方程和步骤S2中的模型参数，计算得到每个循环周次的应力应变张量；
[0011] S5：通过用户子程序USDFLD计算等效应力和等效塑性应变，基于步骤S3中的疲劳损伤计算模型和蠕变损伤计算模型并结合步骤S4中得到的应力应变张量，通过线性累计损伤准则
叠加每个循环周次的疲劳损伤和蠕变损伤，得到待测试材料的裂纹萌生寿命。
[0012] 在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，步骤S4中的用户子程序UMAT可以计算每个
节点的应力应变张量。步骤S5中的用户子程序USDFLD设置的作用为：(1)提取步骤S4的应力应变张量，并进行标量化运算，计算应力应变等效数值；(2)将步骤S3的疲劳损伤计算模型和蠕变损伤计算模型输入用户子程序USDFLD，以获得每个循环周次的疲劳损伤和蠕变损伤，上述损伤的计算可以基于等效应力应变的；(3)对每个循环周次的损伤进行线性累计损伤的叠加，以得到裂纹萌生寿命。
[0013] 进一步地，在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，在步骤S2中，在给定
温度下对待测试材料进行单轴拉伸试验以及该给定温度下不同应变幅和保载时间的单轴蠕变疲劳试验，以得到高温拉伸曲线、循环
软化曲线、应力松弛曲线和
磁滞回线，用以确定所述粘塑性本构方程所需的模型参数。
[0014] 进一步地，在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，在步骤S2中，通过试参法模拟ABAQUS有限元模型的高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线，使之与实验结果吻合，也就是说，通过试参法所获得的曲线与实验所获得的曲线具有较好的拟合程度，以获取粘塑性本构方程所需的模型参数。
[0015] 进一步地，在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，在步骤S1中，粘塑性本构方程包括：粘塑性本构的主控方程、粘塑性本构的粘塑性方程、粘塑性本构的背应力张量非弹性随动强化方程和粘塑性本构的各项同性强化方程。
[0016] 更进一步地，在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，S1步骤进一步包括步骤：
[0017] S11：采用下述公式(1)和公式(2)描述粘塑性本构的主控方程：
[0018] εt＝εe+εin (1)；
[0019]
[0020] 其中，εt为总应变张量，εe为弹性应变张量，εin为非弹性应变张量，E为
弹性模量，v为泊松比，σ为应力张量，trσ为应力张量的迹，I为二阶单位张量；
[0021] S12：采用下述公式(3)、公式(4)和公式(5)描述所述粘塑性本构的粘塑性方程：
[0022]
[0023]
[0024]
[0025] 其中， 为非弹性应变率张量，为累积非弹性应变率，s为应力张量的偏量，a为背应力张量的偏量，J(σ-α)为Von-Mises应力空间距离，α为背应力张量，K和n为率相关的材料参数，R为各向同性
变形抗性，κ为弹性区域初始尺寸，“：”表示张量的内积运算；
[0026] S13：采用下述公式(6)和公式(7)描述所述粘塑性本构的背应力张量非弹性随动强化方程：
[0027]
[0028] m(q)＝φ1e-q/ω+φ2 (7)；
[0029] 其中，αi表示若干个背应力张量部分中的每一个部分，ζi和ri为背应力张量每个部分的材料参数，γ为描述静力回复项的材料参数，m(q)为描述静力回复项的指数方程，q为塑性应变幅值，φ1，φ2和ω为指数方程中三个材料参数，J(αi)表示背应力第二不变量，e表示以自然常数为底的指数函数， 表示背应力张量每个部分的应力变化率；
[0030] S14：采用下述公式(8)描述所述粘塑性本构的各项同性强化方程：
[0031]
[0032] 其中，Q为在第一阶段快速软化的各项同性抗力渐近值，b为接近渐近值的速度参数，H为第二阶段线性软化的斜率相关的参数，p为累积非弹性应变， 表示各项同性强化率。
[0033] 进一步地，在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，其中背应力张量分为8个部分，即
[0034] 进一步地，在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，在所述S3步骤中，多轴应力应变状态的所述纯疲劳损伤计算模型为：
[0035]
[0036] 其中，“()max”表示临界平面上的最大疲劳损伤因子，τmax为临界平面上的最大
剪切应力，τ′f为剪切疲劳强度常数，Δγ/2为临界平面上的
剪切应变幅，σn,max为临界平面上的最大正应力，σ′f为疲劳强度常数，Δεn/2为临界平面上的正应变幅，G为
剪切模量，df为一个循环周次的疲劳损伤，b0为疲劳强度指数，γ′f为剪切疲劳延性常数，c0为疲劳延性指数；
[0037] 多轴应力应变状态的所述蠕变损伤计算模型为：
[0038]
[0039] 其中，dc为一个循环周次的蠕变损伤，th为一个循环周次的保载时间，Z为弹性跟随因子，t表示一个循环周次内从保载开始的时间， 为蠕变损伤第一线性回归参数，MDF为多轴延性因子，n1为蠕变损伤第二线性回归参数，wf,trans为失效应变能
密度平台值， 为一个循环周次保载前最大等效应力，A为松弛第一参数，B为松弛第二参数， 为等效弹性模量， 为一个循环周次由疲劳产生的等效塑性应变范围， 为一个循环周次的等效平均应力，n2为稳态蠕变指数，σH为静
水应力， 表示等效应力。
[0040] 进一步地，在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，在步骤S4中，建立多轴应力应变状态下的ABAQUS有限元模型，施加边界条件和外部载荷，划分模型网格，以得到各个积分点每个循环周次的应力应变张量。
[0041] 进一步地，在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中，步骤S5进一步包括步骤：
[0042] S51：通过用户子程序USDFLD提取ABAQUS模型中每个节点的应力张量和应变张量；
[0043] S52：通过用户子程序USDFLD，并结合所述步骤S51，通过标量化运算获得每一个时刻的等效应力、等效塑性应变和每一个循环周次保载时间内的弹性跟随因子，最终得到每一个循环周次的疲劳损伤和蠕变损伤；
[0044] S53：通过用户子程序USDFLD采用下述公式(11)计算多轴蠕变疲劳条件下的总损伤：
[0045]
[0046] 其中，D(n)为前n个循环周次的累积总损伤， 为第i个循环周次产生的疲劳损伤， 为第i个循环周次产生的蠕变损伤；
[0047] 其中，当某个节点的总损伤叠加率先到失效数值1时，即可定义为最危险节点并确定裂纹萌生寿命ni。
[0048] 需要说明的是，在本发明所述的技术方案中，裂纹萌生寿命采用循环周次ni进行表征。
[0049] 本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法具有如下所述的优点和有益效果：
[0050] (1)本发明所述的多轴蠕变疲劳预测方法利用基于用户子程序UMAT定义待测试材料的粘塑性本构方程，从而得到多轴应力应变状态下的蠕变疲劳行为；
[0051] (2)本发明所述的多轴蠕变疲劳预测方法利用基于用户子程序USDFLD计算等效应力和等效塑性应变，从而获得每个循环周次中各个积分点的蠕变损伤、疲劳损伤和总损伤数值；
[0052] (3)本发明所述的多轴蠕变疲劳预测方法具有很强的直观性，可以直观地获得几何不连续结构的裂纹萌生
位置和该位置的裂纹萌生寿命。
附图说明
[0053] 图1为本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的
流程图。
[0054] 图2示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴拉伸试验和模拟曲线的拟合结果图。
[0055] 图3示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验循环软化数据和模拟曲线的拟合结果图。
[0056] 图4示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验应力松弛数据和模拟曲线的拟合结果图。
[0057] 图5示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验静力回复项的指数方程拟合结果图。
[0058] 图6示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的两个潜在危险点每周次疲劳损伤和每周次蠕变损伤随循环周次的轨迹图。
[0059] 图7示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的某一个潜在危险点前100个循环周次的磁滞回线。
[0060] 图8示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的另一个潜在危险点前100个循环周次的磁滞回线。
[0061] 图9示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的缺口根部次表面最危险积分点的蠕变疲劳损伤轨迹图和裂纹萌生寿命预测。
[0062] 图10示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的两个潜在危险点每周次疲劳损伤和每周次蠕变损伤随循环周次的轨迹图。
[0063] 图11示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的某一个潜在危险点前100个循环周次的磁滞回线。
[0064] 图12示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的另一个潜在危险点前100个循环周次的磁滞回线。
[0065] 图13示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的缺口根部表面的蠕变疲劳损伤轨迹图和裂纹萌生寿命预测。
具体实施方式
[0066] 下面将结合
说明书附图和具体的
实施例对本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法做进一步的详细说明。
[0067] 图1为本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的流程图。
[0068] 如图1所示，在本实施方式中，基于ABAQUS的多轴蠕变疲劳预测方法包括步骤：
[0069] S1：建立ABAQUS有限元模型，通过用户子程序UMAT定义待测试材料在循环载荷过程中的粘塑性本构方程；
[0070] S2：确定所述粘塑性本构方程所需的模型参数；
[0071] S3：建立待测试材料的多轴应力应变状态的疲劳损伤计算模型和蠕变损伤计算模型；
[0072] S4：建立多轴应力应变状态下的ABAQUS有限元模型，基于步骤S1中用户子程序UMAT所定义的粘塑性本构方程和步骤S2中的模型参数，计算得到每个循环周次的应力应变张量；
[0073] S5：通过用户子程序USDFLD计算等效应力和等效塑性应变，基于步骤S3中的疲劳损伤计算模型和蠕变损伤计算模型并结合步骤S4中得到的应力应变张量，通过线性累计损伤准则叠加每个循环周次的疲劳损伤和蠕变损伤，得到待测试材料的裂纹萌生寿命。
[0074] 其中，在步骤S2中，在给定温度下对待测试材料进行单轴拉伸试验以及该给定温度下不同应变幅和保载时间的单轴蠕变疲劳试验，以得到高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线，用以确定所述粘塑性本构方程所需的模型参数，并且可以通过试参法模拟ABAQUS有限元模型的高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线，以获取粘塑性本构方程所需的模型参数。
[0075] 其中，在步骤S1中，粘塑性本构方程包括：粘塑性本构的主控方程、粘塑性本构的粘塑性方程、粘塑性本构的背应力张量非弹性随动强化方程和粘塑性本构的各项同性强化方程。其可以通过下述步骤获得：
[0076] S11：采用下述公式(1)和公式(2)描述所述粘塑性本构的主控方程：
[0077] εt＝εe+εin (1)；
[0078]
[0079] 其中，εt为总应变张量，εe为弹性应变张量，εin为非弹性应变张量，E为弹性模量，v为泊松比，σ为应力张量，trσ为应力张量的迹，I为二阶单位张量；
[0080] S12：采用下述公式(3)、公式(4)和公式(5)描述所述粘塑性本构的粘塑性方程：
[0081]
[0082]
[0083]
[0084] 其中， 为非弹性应变率张量，为累积非弹性应变率，s为应力张量的偏量，a为背应力张量的偏量，J(σ-α)为Von-Mises应力空间距离，α为背应力张量，K和n为率相关的材料参数，R为各向同性变形抗性，κ为弹性区域初始尺寸，“：”表示张量的内积运算；
[0085] S13：采用下述公式(6)和公式(7)描述所述粘塑性本构的背应力张量非弹性随动强化方程：
[0086]
[0087] m(q)＝φ1e-q/ω+φ2 (7)；
[0088] 其中，αi表示若干个背应力张量部分中的每一个部分，ζi和ri为背应力张量每个部分的材料参数，γ为描述静力回复项的材料参数，m(q)为描述静力回复项的指数方程，q为塑性应变幅值，φ1，φ2和ω为指数方程中三个材料参数，J(αi)表示背应力第二不变量，e表示以自然常数为底的指数函数， 表示背应力张量每个部分的应力变化率；
[0089] S14：采用下述公式(8)描述所述粘塑性本构的各项同性强化方程：
[0090]
[0091] 其中，Q为在第一阶段快速软化的各项同性抗力渐近值，b为接近渐近值的速度参数，H为第二阶段线性软化的斜率相关的参数，p为累积非弹性应变， 表示各项同性强化率。
[0092] 需要说明的是，上述步骤中，背应力张量分为8个部分，即
[0093] 而在S3步骤中，多轴应力应变状态的所述纯疲劳损伤计算模型为：
[0094]
[0095] 其中，“( )max”表示临界平面上的最大疲劳损伤因子，τmax为临界平面上的最大剪切应力，τ′f为剪切疲劳强度常数，Δγ/2为临界平面上的剪切应变幅，σn,max为临界平面上的最大正应力，σ′f为疲劳强度常数，Δεn/2为临界平面上的正应变幅，G为剪切模量，df为一个循环周次的疲劳损伤，b0为疲劳强度指数，γ′f为剪切疲劳延性常数，c0为疲劳延性指数；
[0096] 多轴应力应变状态的所述蠕变损伤计算模型为：
[0097]
[0098] 其中，dc为一个循环周次的蠕变损伤，th为一个循环周次的保载时间，Z为弹性跟随因子，t表示一个循环周次内从保载开始的时间， 为蠕变损伤第一线性回归参数，MDF为多轴延性因子，n1为蠕变损伤第二线性回归参数，wf,trans为失效应变能密度平台值， 为一个循环周次保载前最大等效应力，A为松弛第一参数，B为松弛第二参数， 为等效弹性模量， 为一个循环周次由疲劳产生的等效塑性应变范围， 为一个循环周次的等效平均应力，n2为稳态蠕变指数，σH为静水应力， 表示等效应力。
[0099] 而在步骤S4中，建立多轴应力应变状态下的ABAQUS有限元模型，可以施加边界条件和外部载荷，划分模型网格，以得到各个积分点每个循环周次的应力应变张量[0100] 在本实施方式中，步骤S5可以进一步包括步骤
[0101] S51：通过用户子程序USDFLD提取ABAQUS模型中每个节点的应力张量和应变张量；
[0102] S52：通过用户子程序USDFLD，并结合所述步骤S51，通过标量化运算获得每一个时刻的等效应力、等效塑性应变和每一个循环周次保载时间内的弹性跟随因子，最终得到每一个循环周次的疲劳损伤和蠕变损伤；
[0103] S53：通过用户子程序USDFLD采用下述公式(11)计算多轴蠕变疲劳条件下的总损伤：
[0104]
[0105] 其中，D(n)为前n个循环周次的累积总损伤， 为第i个循环周次产生的疲劳损伤， 为第i个循环周次产生的蠕变损伤；
[0106] 其中，当某个节点的总损伤叠加率先到失效数值1时，即可定义为最危险节点并确定裂纹萌生寿命ni。
[0107] 为了更好地阐述本案的基于AVAQUS的多轴蠕变疲劳预测方法的预测效果，将采用缺口半径为8毫米的单边缺口试样进行验证。
[0108] 验证所采用的单边缺口试样的材料为高温镍基GH4169超
合金，蠕变疲劳试验在650℃空气环境下进行。试验时，试样两端施加的外载和为整体的应变控制，由于单边缺口试样几何不连续性的影响，其缺口根本最薄弱的部位处于多轴应力应变状态。在此之前，需要进行650℃空气环境下的单轴拉伸试验以及该环境下不同应变幅和保载时间的单轴蠕变疲劳试验，获得的试验结果用于确定公式(1)～(8)粘塑性本构方程所需的材料参数。
[0109] 通过试参法调整单轴拉伸试验的模拟结果，使之与单轴拉伸试验能够较好地吻合。图2示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴拉伸试验和模拟曲线的拟合结果图。
[0110] 如图2所示，I表示单轴拉伸试验数据，II表示采用试参法模拟得到的单轴拉伸试验曲线，通过试参法调整单轴拉伸试验曲线II的模拟结果，使之与单轴蠕变疲劳试验所得的单轴拉伸试验数据能够较好地吻合，其中，弹性模量E＝177GPa,泊松比v＝0.33,弹性区域初始尺寸κ＝815MPa,各个背应力张量每个部分的随动强化材料参数：ζ1＝6130,ζ2＝1807,ζ3＝892,ζ4＝352,ζ5＝150,ζ6＝88.2,ζ7＝75.0,ζ8＝28.4,r1＝23.4,r2＝68.0,r3＝
75.9,r4＝48.0,r5＝43.4,r6＝25.4,r7＝54.5,r8＝28.0,率相关的粘塑性材料参数：K＝
400,n＝2.0。
[0111] 图3示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验循环软化数据和模拟曲线的拟合结果图。
[0112] 如图3所示，III表示循环软化试验数据，IV表示采用试参法模拟得到的循环软化曲线，通过试参法调整循环软化曲线IV的模拟结果，使之与单轴蠕变疲劳试验所得的循环软化数据能够较好地吻合，其中，各项同性强化的材料参数：第二阶段线性软化的斜率相关的参数H＝-8.5，接近渐近值的速度参数b＝4.1，第一阶段快速软化的各项同性抗力渐近值Q＝618MPa。
[0113] 图4示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验应力松弛数据和模拟曲线的拟合结果图。
[0114] 如图4所示，V表示单轴蠕变疲劳试验应力松弛数据，VI表示采用试参法模拟得到的应力松弛曲线，通过试参法调整第一个循环周次的应力松弛曲线的模拟结果，使之与单轴蠕变疲劳试验所得的应力松弛数据能够较好地吻合，其中，静力回复项中独立材料参数：γ＝4.0×10-7。
[0115] 与此同时，拟合公式(7)中的静力回复项的指数方程，拟合结果图如图5所示。图5示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验静力回复项的指数方程拟合结果图。
[0116] 如图5所示，VII表示单轴蠕变疲劳试验数据，VIII表示模拟结果曲线，其中，静力回复项中指数方程的材料参数：φ1＝0.37，φ2＝2.82，ω＝6.6×10-4。
[0117] 在本案的一个实施方式中，选取数据为总应变范围为1.0％，每个循环周次应变最大处的保载时间为3600s，650℃空气环境下的单边缺口蠕变疲劳试验。其蠕变疲劳裂纹萌生于缺口根部次表面的位置，裂纹萌生寿命为76个循环周次。
[0118] 图6给出了该实施方式下两个典型位置每周次疲劳损伤和蠕变损伤的轨迹图。
[0119] 如图6所示，曲线IX表示缺口根部表面每周次疲劳损伤曲线，曲线X表示缺口根部次表面每周次疲劳损伤曲线，曲线XI表示缺口根部表面每周次蠕变损伤曲线，曲线XII表示缺口根部次表面每周次蠕变损伤曲线，其中，选取缺口根部表面的积分点，因其在循环载荷的过程中受到应力集中的影响，疲劳裂纹通常萌生于表面，所以缺口根部表面的积分点通常被认为是疲劳过程中的潜在危险点；而选取缺口根部次表面的积分点，是由于次表面通常具有相比于表面更高的应力三轴度，蠕变裂纹通常萌生于内部，所以缺口根部次表面的积分点通常被认为是蠕变过程中的潜在危险点。而从图6中可以看出：此实施方式下由于保载时间较长，循环载荷的过程中蠕变损伤占主导地位，所计算最大总损伤的位置在缺口根部次表面，与实验所观测的现象完全一致。
[0120] 图7和图8给出了图6所示的缺口根部表面和次表面两个位置前100个循环周次的磁滞回线。其中，图7表示缺口根部表面积分点处的磁滞回线，图8表示缺口根部次表面积分点处的磁滞回线，图7和图8中的A1表示第一个循环周次，A2表示第一百个循环周次。
[0121] 结合图7和图8中可以看出：虽然缺口根部表面积分点具有更大的非弹性应变范围，但是次表面积分点具有更加明显的非弹性应变的累积。这是由于受到多轴应力应变状态的影响，几何不连续处在保载阶段处于应力应变混合控制的模式。结合图7和图8还可以看出，次表面积分点更加接近应力控
制模式，导致其蠕变损伤更加明显，进而解释了蠕变主导的循环载荷其裂纹萌生位置出现在缺口根部次表面。
[0122] 图9示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的缺口根部次表面最危险积分点的蠕变疲劳损伤轨迹图和裂纹萌生寿命预测。
[0123] 借助蠕变疲劳损伤交互图，追踪缺口根部附近最危险的次表面积分点，如图9所示，可以获得该实施方式下的裂纹萌生寿命为105个循环周次，这与实验测得的76个循环周次较为接近，在1.5倍误差带范围之内证明此数值模拟方法具有较高的可靠性。
[0124] 需要说明的是，图9中的曲线XIII表示
[0125] 在本案的另一个实施方式中，选取数据为总应变范围为1.0％，每个循环周次应变最大处的保载时间为60s，650℃空气环境下的单边缺口蠕变疲劳试验。其蠕变疲劳裂纹萌生于缺口根部表面的位置，裂纹萌生寿命为480个循环周次。
[0126] 图10给出了该实施方式下两个典型位置每周次疲劳损伤和蠕变损伤的轨迹图。
[0127] 如图10所示，曲线XIV表示缺口根部表面每周次疲劳损伤曲线，曲线XV表示缺口根部次表面每周次疲劳损伤曲线，曲线XVI表示缺口根部表面每周次蠕变损伤曲线，曲线XVII表示缺口根部次表面每周次蠕变损伤曲线，其中，选取缺口根部表面的积分点，是因其在循环载荷的过程中受到应力集中的影响，疲劳裂纹通常萌生于表面，所以缺口根部表面的积分点通常被认为是疲劳过程中的潜在危险点；选取缺口根部次表面的积分点，是由于次表面通常具有相比于表面更高的应力三轴度，蠕变裂纹通常萌生于内部，所以缺口根部次表面的积分点通常被认为是蠕变过程中的潜在危险点。而从图10中可以看出：在此实施方式下由于保载时间较短，循环载荷的过程中疲劳损伤占主导地位，所计算最大总损伤的位置在缺口根部表面，与实验所观测的现象完全一致。
[0128] 图11和图12给出了图10所示的缺口根部表面和次表面两个位置前100个循环周次的磁滞回线。其中，图11表示缺口根部表面积分点处的磁滞回线，图12表示缺口根部次表面积分点处的磁滞回线，图11和图12中的A1表示第一个循环周次，A2表示第一百个循环周次。
[0129] 结合图11和图12可以看出，缺口根部表面积分点具有更大的非弹性应变范围，而且缺口根本表面和次表面两个位置的积分点具有几乎没有非弹性应变的累积。结合图11和图12可以看出，该实施方式下保载时间较短导致的蠕变/松弛的现象几乎可以忽略不计，在这种情况下，由非弹性应变范围导致的疲劳损伤使得裂纹萌生位置出现在缺口根部表面。
[0130] 图13示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的缺口根部表面的蠕变疲劳损伤轨迹图和裂纹萌生寿命预测。
[0131] 借助蠕变疲劳损伤交互图，追踪缺口根部表面最危险的积分点，如图13所示，可以获得该实施方式下的裂纹萌生寿命为346个循环周次，这与实验测得的480个循环周次较为接近，在1.5倍误差带范围之内证明此数值模拟方法具有较高的可靠性。
[0132] 需要说明的是，图13中的曲线XVIII表示
[0133] 综上所述可以看出，本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法利用基于用户子程序UMAT定义待测试材料的粘塑性本构方程，从而得到多轴应力应变状态下的蠕变疲劳行为。
[0134] 此外，本发明所述的多轴蠕变疲劳预测方法利用基于用户子程序USDFLD计算等效应力和等效塑性应变，从而获得每个循环周次中各个积分点的蠕变损伤、疲劳损伤和总损伤数值；
[0135] 另外，本发明所述的多轴蠕变疲劳预测方法具有很强的直观性，可以直观地获得几何不连续结构的裂纹萌生位置和该位置的裂纹萌生寿命。
[0136] 需要说明的是，本发明的保护范围中
现有技术部分并不局限于本
申请文件所给出的实施例，所有不与本发明的方案相矛盾的现有技术，包括但不局限于在先
专利文献、在先公开出版物，在先公开使用等等，都可纳入本发明的保护范围。
[0137] 此外，本案中各技术特征的组合方式并不限本案
权利要求中所记载的组合方式或是具体实施例所记载的组合方式，本案记载的所有技术特征可以以任何方式进行自由组合或结合，除非相互之间产生矛盾。
[0138] 还需要注意的是，以上所列举的实施例仅为本发明的具体实施例。显然本发明不局限于以上实施例，随之做出的类似变化或变形是本领域技术人员能从本发明公开的内容直接得出或者很容易便联想到的，均应属于本发明的保护范围。
