适用于预测多种控制模式循环载荷变形行为的方法技术

本发明专利技术公开了一种适用于预测多种控制模式循环载荷变形行为的方法，取相同材料的四根试样，在相同的温度下进行低周疲劳试验、应变控蠕变疲劳试验以及应力应变混合控蠕变疲劳试验，获得相应的数据；建立统一粘塑性本构模型，包括统一粘塑性本构模型的应变分解、流动率、各向同向硬化以及运动硬化；利用试验数据作为基础，分步简化提出的统一粘塑性本构模型，最终确定所有的模型参数；利用确定的模型参数确定最终统一粘塑性模型，并利用此模型来预测相同材料在其他多种控制模式循环载荷下的变形行为。本发明专利技术具有操作简单，精度高，且适用性广的特点，可以采用一组参数同时预测多种循环载荷下的变形行为。循环载荷下的变形行为。循环载荷下的变形行为。

【技术实现步骤摘要】
适用于预测多种控制模式循环载荷变形行为的方法


[0001]本专利技术涉及本构模型和循环变形预测领域，尤其涉及一种适用于预测多种控制模式循环载荷变形行为的方法。

技术介绍

[0002]本构模型是用来描述材料在循环载荷下的各种非弹性变形行为一种理论方法，如循环硬化，循环软化，蠕变变形，应力松弛等。近年来，随着计算机技术的快速发展，通过大型有限元软件对本构模型进行二次开发，可以对实际受载部件的变形行为以及损伤分布进行实时监测。但这一先进技术对本构模型的精度及普适性提出了更高的要求。
[0003]尽管目前已经发展出一些较为成熟的本构模型，例如Chaboche模型，Robinson模型等，但这些模型只能描述单纯地应变控或者单纯地应力控下的变形行为。然而在实际的航空和电力领域中，众多关键高温核心部件并不仅仅简单地承受一种控制模式的载荷。通常，频繁开停车产生的温度变化使得这些部件处于应变控制的疲劳载荷中，而装备持续稳定运行过程中的恒定内压或者离心力，使高温部件承受恒应力控制的蠕变载荷作用。因此，这些关重件往往承受着应力应变混合控制的蠕变疲劳交互载荷。目前很难找到一个合理的本构模型来描述如此复杂控制模式下的变形行为。此外，更没有一种方法能够通过一组参数描述并预测多种控制模式循环载荷(低周疲劳载荷，应变控蠕变疲劳载荷，应力应变混合控蠕变疲劳载荷，应力控棘轮蠕变载荷)下的变形行为。

技术实现思路

[0004]本专利技术的目的在于提供一种适用于预测多种控制模式循环载荷变形行为的方法，可以实现材料在多种控制模式循环载荷下的变形行为预测，且具有适用性广，精度高，模型参数少等优点。
[0005]实现本专利技术目的的技术解决方案为：一种适用于预测多种控制模式循环载荷变形行为的方法，包括以下步骤：
[0006]步骤S1，取相同材料的四根试样，在相同的温度下进行低周疲劳试验、应变控蠕变疲劳试验以及应力应变混合控蠕变疲劳试验，获得相应的数据，转入步骤S2。
[0007]步骤S2，建立统一粘塑性本构模型，包括统一粘塑性本构模型的应变分解、流动率、各向同向硬化以及运动硬化，转入步骤S3。
[0008]步骤S3，利用步骤S1中的数据作为基础，分步简化步骤S2中的统一粘塑性本构模型，最终确定所有的模型参数，转入步骤S4。
[0009]步骤S4，利用步骤S3中得到的模型参数确定最终统一粘塑性模型，通过最终统一粘塑性模型预测相同材料在其他多种控制模式循环载荷下的变形行为。
[0010]本专利技术与现有技术相比，其显著优点在于：
[0011](1)本专利技术模型参数较少，共19个，目前现有技术的模型参数基本维持在25个以上。
[0012](2)本专利技术参数确定过程较为简单，易于操作。
[0013](3)本专利技术可以通过一组参数较为精确地预测低周疲劳，应变控制的蠕变疲劳，应力应变混合控制的蠕变疲劳以及应力控棘轮蠕变载荷等多种控制模式下的循环变形行为及蠕变行为。
附图说明
[0014]图1为本专利技术的适用于预测多种控制模式循环载荷变形行为的方法的流程图。
[0015]图2为本专利技术实施例中预测的低周疲劳载荷下第一周次以及半寿命周次的滞回曲线图。
[0016]图3为本专利技术实施例中预测的应变控蠕变疲劳载荷下第一周次以及半寿命周次的滞回曲线图。
[0017]图4为本专利技术实施例中预测的应变控蠕变疲劳载荷下第一周次以及半寿命周次保载期间的应力松弛曲线图。
[0018]图5为本专利技术实施例中预测的应变控蠕变疲劳载荷下峰值应力以及松弛应力随周次的演化曲线图。
[0019]图6为本专利技术实施例中预测的应力应变混合控蠕变疲劳载荷下第一周次以及半寿命周次的滞回曲线图。
[0020]图7为本专利技术实施例中预测的应力应变混合控蠕变疲劳载荷下第一周次以及半寿命周次保载期间的蠕变应变曲线图。
[0021]图8为本专利技术实施例中预测的应力应变混合控蠕变疲劳载荷下每周蠕变应变随周次的演化曲线图。
[0022]图9为本专利技术实施例中预测的应力控棘轮蠕变载荷下保载期间的蠕变应变曲线图。
具体实施方式
[0023]以下结合附图，以具体实施例对本专利技术作进一步详细说明。应该理解，以下实施例仅用于说明本专利技术而非用于限定本专利技术的范围。
[0024]本专利技术所述的适用于预测多种控制模式循环载荷变形行为的方法是基于Walker循环本构模型的。本专利技术在此基础上，考虑材料在循环载荷下的连续软化行为，修正各向同性硬化变量，并根据不同应力应变混合控制蠕变疲劳载荷下显著区别的蠕变行为，提出了自适应的运动硬化变量。同时针对循环本构模型参数较多的特点，开发出一套简单快捷的参数确定方法。基于该专利技术提出的方法，可以预测多种控制模式循环载荷下的变形行为。
[0025]本专利技术可以通过一组参数较为精确地预测低周疲劳，应变控制的蠕变疲劳，应力应变混合控制的蠕变疲劳以及应力控棘轮蠕变载荷等多种控制模式循环载荷下的循环变形行为及蠕变行为。与现有技术相比，本专利技术具有操作简单，精度高，且适用性广的特点。
[0026]结合图1，本专利技术所述的适用于预测多种控制模式循环载荷变形行为的方法，包括以下的步骤：
[0027]步骤S1，取相同材料的四根试样，在相同的温度下进行多种控制模式的循环试验。其中，一根试样进行低周疲劳试验，遵循标准【国标GB/T 26077
‑
2010】，获得第一周次的滞
回曲线以及峰值应力随周次的演化曲线。第二根试样采用与低周疲劳试验相同的应变速率以及应变幅值，进行应变控蠕变疲劳试验，遵循标准【国标GB/T 38822
‑
2020】，获得每周松弛应力随周次的演化趋势。剩余两根试样进行应力应变混合控制的蠕变疲劳试验，遵循文献【Zhang T,Wang X,Ji Y,et al.Cyclic deformation and damage mechanisms of 9％Cr steel under hybrid stress
‑
strain controlled creep fatigue interaction loadings[J].International Journal of Fatigue,2021,151:106357.】，两根试样的疲劳加载阶段相同，依然采用与低周疲劳试验相同的应变速率以及应变幅值，而两根试样的蠕变阶段不同，第三根采用小保载应力短保载时间(小保载应力：该应力的蠕变寿命在1000小时以上，短保载时间：小于300s)，第四根采用大保载应力长保载时间(大保载应力：该应力的蠕变寿命在500小时以下，长保载时间：大于300s)，从而获得两种应力应变混合控蠕变疲劳加载下第一周次保载时蠕变应变随时间的演化以及每周蠕变应变量随周次的演化。
[0028]步骤S2，建立统一粘塑性本构模型，包括统一粘塑性本构模型的应变分解、流动率、各向同向硬化以及运动硬化，其具体步骤如下：
[0029]步骤S21，统一粘塑性本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种适用于预测多种控制模式循环载荷变形行为的方法，其特征在于，包括以下步骤：步骤S1，取相同材料的四根试样，在相同的温度下进行低周疲劳试验、应变控蠕变疲劳试验以及应力应变混合控蠕变疲劳试验，获得相应的数据，转入步骤S2；步骤S2，建立统一粘塑性本构模型，包括统一粘塑性本构模型的应变分解、流动率、各向同向硬化以及运动硬化，转入步骤S3；步骤S3，利用步骤S1中的数据作为基础，分步简化步骤S2中的统一粘塑性本构模型，最终确定所有的模型参数，转入步骤S4；步骤S4，利用步骤S3中得到的模型参数确定最终统一粘塑性模型，通过最终统一粘塑性模型预测相同材料在其他多种控制模式循环载荷下的变形行为。2.根据权利要求1所述的适用于预测多种控制模式循环载荷变形行为的方法，其特征在于，所述步骤S1中，取相同材料的四根试样，在相同的温度下进行多种控制模式的循环试验：其中，一根试样进行低周疲劳试验，获得第一周次的滞回曲线以及峰值应力随周次的演化曲线；第二根试样采用与低周疲劳试验相同的应变速率以及应变幅值进行应变控蠕变疲劳试验，获得每周松弛应力随周次的演化趋势；剩余两根试样进行应力应变混合控制的蠕变疲劳试验，上述两根试样的疲劳加载阶段相同，依然采用与低周疲劳试验相同的应变速率以及应变幅值，而两根试样的蠕变阶段不同，其中第三根采用小保载应力短保载时间，第四根采用大保载应力长保载时间，从而获得两种应力应变混合控蠕变疲劳加载下第一周次保载时蠕变应变随时间的演化以及每周蠕变应变量随周次的演化。3.根据权利要求2所述的适用于预测多种控制模式循环载荷变形行为的方法，其特征在于，所述步骤S2中，统一粘塑性本构模型的应变分解如式(1)所示：ε＝ε
e
+ε
in
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)ε
e
＝σ/E
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)其中ε表示试样承受的总应变，ε
in
表示试样承受的非弹性应变，ε
e
表示试样承受的弹性应变，σ为试样承受的应力，E为试样的弹性模量。4.根据权利要求3所述的适用于预测多种控制模式循环载荷变形行为的方法，其特征在于，所述步骤S2中，统一粘塑性本构模型的流动率如式(3)所示：其中为非弹性应变ε
in
对时间的导数，A表示材料的流动系数，n表示材料的流动指数，χ表示运动硬化，D表示各向同性硬化，函数sgn(*)表示符号函数，当*>0时，sgn(*)＝1；当*＝0时，sgn(*)＝0；当*<0时，sgn(*)＝
‑
1。5.根据权利要求4所述的适用于预测多种控制模式循环载荷变形行为的方法，其特征在于，所述步骤S2中，统一粘塑性本构模型的各向同向硬化D，其演化率如式(4)所示：
当处于零时刻时，D(t＝0)＝D0，D0表示初始状态下各向同性硬化，D0为常数，表示非弹性应变速率的绝对值，p表示循环过程中的累积非弹性变形，C
d
、γ
d
、H和D0分别对应各向同性硬化系数、各向同性恢复系数、线性恢复系数以及初始各向同性硬化。6.根据权利要求5所述的适用于预测多种控制模式循环载荷变形行为的方法，其特征在于，所述步骤S2中，统一粘塑性本构模型的运动硬化χ分解成两部分，如式(5)所示：χ＝χ1+χ2ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)第一运动硬化χ1的演化率和第二运动硬化χ2的演化率分别遵循式(6)和式(7)：分别遵循式(6)和式(7)：其中C1表示第一运动硬化的强化系数，C2表示第二运动硬化的强化系数，γ1表示第一运动硬化的恢复系数，γ2表示第二运动硬化的恢复系数；S1表示第一运动硬化的恢复因子，S2表示第二运动硬化的恢复因子，其表示为式(8)：其中，k1表示第一运动硬化的恢复底数，k2表示第二运动硬化的恢复底数，m1表示第一运动硬化的恢复指数，m2表示第二运动硬化的恢复指数，令k1＝k2＝k，m1＝m2＝m，k表示统一粘塑性本构模型的恢复底数，m表示统一粘塑性本构模型的恢复指数；式(7)中的累积非弹性变形相关的强化系数ψ(p)表示为式(9)：式(7)中的累积非弹性变形相关的强化系数ψ(p)表示为式(9)：式(7)中的累积非弹性变形相关的强化系数ψ(p)表示为式(9)：其中极限强化系数ψ
∞
(χ2)和演化率ω(χ2)分别遵循式(10)和式(11)的规律...

【专利技术属性】
技术研发人员：王小威，张天宇，周德文，姜勇，张显程，巩建鸣，涂善东，
申请(专利权)人：南京工业大学，
类型：发明
国别省市：