CN106934168B - 一种材料多轴蠕变失效应变预测方法 - Google Patents

Abstract

一种材料多轴蠕变失效应变预测方法，属于应变预测技术领域。其特征在于：包括如下步骤：步骤（1），获得材料蠕变速率和应变速率之间的关系；步骤（2），由幂律蠕变控制孔洞长大理论，获得单轴及多轴应力状态参数；步骤（3），表示单轴及多轴应力作用下的蠕变失效应变，并获得多轴蠕变延性因子；步骤（4），获得不同应力状态下的多轴蠕变失效应变拟合参数，从而得到多轴蠕变延性因子预测方程；步骤（5），利用有限元软件预测材料在多轴应力状态下的蠕变失效应变及其寿命。本材料多轴蠕变失效应变预测方法能够通过材料在多轴应力状态下的蠕变失效应变计算方法，更准确地预测材料在高温状态下的多轴蠕变的失效应变。

Description

一种材料多轴蠕变失效应变预测方法

技术领域

一种材料多轴蠕变失效应变预测方法，属于应变预测技术领域。

背景技术

在核电、石油化工和航空航天等领域中，许多结构部件如换热器等长期在高温高压下工作，整个结构处于复杂的多轴应力状态，蠕变及其引起的损伤是结构的主要失效方式之一。多轴应力状态下的蠕变-损伤失效研究是结构完整性评定中最为重要的环节之一，因此对多轴应力状态下材料的蠕变-损伤进行研究，对高温高压工作条件下部件的寿命预测具有积极的意义。

针对模型参数繁冗的问题，基于应变的连续损伤力学模型日渐得到人们的重视。基于应变损伤的模型又称延性耗竭模型，认为局部蠕变应变累积达到蠕变延性（蠕变失效应变）值时，材料将完全损伤，出现裂纹直至失效。多轴蠕变失效应变的测试较为困难，成本也较高，其值通常利用单轴蠕变失效应变转换而来。经典的单轴与多轴蠕变失效应变的转化关系是Cocks-Ashby提出的转换关系。根据转换关系获得的经验多轴蠕变延性因子（MCDF）被广泛的应用于蠕变损伤和蠕变寿命的预测。虽然多轴蠕变延性因子可以很容易地被用来模拟蠕变裂纹扩展并给出可接受的预测结果，但它缺乏应有的物理意义，并且在某些情况下，Cocks-Ashby MCDF对多轴蠕变失效应变的预测并不合理，因而亟需改进。

发明内容

本发明要解决的技术问题是：克服现有技术的不足，提供一种能更准确的预测材料在高温状态下的多轴蠕变的失效应变的材料多轴蠕变失效应变预测方法。

本发明解决其技术问题所采用的技术方案是：该材料多轴蠕变失效应变预测方法，其特征在于：包括如下步骤：

步骤（1），基于应变损伤准则，根据微观孔洞长大机制，获得材料蠕变速率和应变速率之间的关系；

步骤（2），由幂律蠕变控制孔洞长大理论，获得单轴及多轴应力状态参数；

步骤（3），对微观孔洞进行积分获得材料恒载下的失效时间，进而表示单轴及多轴应力作用下的蠕变失效应变，并获得多轴蠕变延性因子；

步骤（4），对多轴蠕变延性因子进行参数拟合，获得不同应力状态下的多轴蠕变失效应变拟合参数，从而得到多轴蠕变延性因子预测方程；

步骤（5），利用有限元软件预测材料在多轴应力状态下的蠕变失效应变及其寿命。

优选的，步骤（1）中所述的微观孔洞长大机制在裂尖区域为粘塑性孔洞长的机制。

优选的，步骤（1）中所述的材料蠕变速率和应变速率之间的关系为：

，

其中，为包含一个孔洞的圆柱体的轴向应变速率，为包含一个孔洞的圆柱体的径向应变速率，为圆柱体不包含孔洞时的稳态蠕变速率。

优选的，步骤（2）中所述的单轴及多轴应力状态参数利用能量原理，并通过New-Raphson方法获得。

优选的，步骤（2）中所述的单轴及多轴应力状态参数的计算公式为：

，

，

其中，为多轴应力状态参数，为单轴应力状态参数，与为与材料相关的常数，为确定的材料常数，为应力三轴度。

优选的，步骤（3）中所述的多轴蠕变延性因子的计算公式如下：

，

其中，为单轴应力作用下的蠕变失效应变，为多轴应力作用下的蠕变失效应变，和均为与材料相关的常数，为确定的材料常数，为应力三轴度。

优选的，利用缺口蠕变试验获得不同应力状态下的蠕变失效应变拟合参数和，进而得到步骤（4）中所述的多轴蠕变延性因子预测方程。

优选的，步骤（4）中所述的多轴蠕变延性因子预测方程结合蠕变-损伤本构方程，并通过Fortran语言嵌入到有限元软件中。

与现有技术相比，本发明所具有的有益效果是：

1、本材料多轴蠕变失效应变预测方法能够通过材料在多轴应力状态下的蠕变失效应变计算方法，定义多轴蠕变延性因子，获得多轴蠕变延性因子计算方程，进而可获得材料的多轴蠕变失效应变，根据建立的多轴蠕变损伤本构方程，利用有限元软件预测材料在多轴应力状态下的蠕变失效应变，能更准确地预测材料在高温状态下的多轴蠕变的失效应变。

2、本发明根据幂率蠕变控制孔洞长大理论，利用能量原理，定义多轴蠕变延性因子，为多轴蠕变延性因子赋予相应的物理意义，获得新的多轴蠕变延性因子预测公式，并利用Fortran语言，根据建立的多轴蠕变损伤本构方程编制子程序并嵌入到有限元软件中，从而能更准确地预测材料在高温状态下的多轴蠕变的失效应变。

附图说明

图1为材料多轴蠕变失效应变预测方法的流程图。

图2为受拉伸载荷作用的含裂纹试样模型图。

图3为含裂纹试样裂尖区域模型图。

图4为多轴载荷作用下的理想晶粒模型图。

图5为包含一个孔洞的圆柱体的放大示意图。

图6为600℃下，9Cr-1Mo合金不同转化模型获得的多轴蠕变失效应变与试验值比较图。

图7为圆柱缺口拉伸试样的几何尺寸图。

具体实施方式

图1~7是本发明的最佳实施例，下面结合附图1~7对本发明做进一步说明。

如图1所示，一种材料多轴蠕变失效应变预测方法，其特征在于：包括如下步骤：

步骤（1），基于应变损伤准则，根据微观孔洞长大机制，获得材料蠕变速率和应变速率之间的关系；

微观孔洞长大机制认为孔洞主要在晶界面上（尤其是与拉应力垂直的晶界面上）形核和长大，充分长大的孔洞将聚合进而形成晶粒尺寸大小的微裂纹（孔化晶界面）。最后，微裂纹合并导致宏观蠕变裂纹的扩展。

尽管空位凝聚，晶界滑移及位错堆积通常被认为是空洞形核的可能驱动力，但孔洞形核背后的真正机制还不清晰。所以本实施例中将主要考虑晶界上孔洞长大与合并的过程。

孔洞长大的机制也有许多种。其中，粘塑性孔洞长大，扩散控制孔洞长大和受约束扩散孔洞长大是三种被广泛接受的孔洞长大模型。而具体哪一种机制起关键作用则取决于材料性能、温度和应力水平等。例如，在高应变速率和高应力下，孔洞的长大更倾向于受周边晶粒的蠕变或塑性变形控制；而在低应力水平下，晶界上的空位扩散则可能是主要原因。如图2~3所示，在裂尖附近的区域，由于宏观裂纹的存在，局部应力和应变总体上保持在一个很高的水平。因此，粘塑性孔洞长大机制将在裂尖区域起到主导作用。

如图4所示：在多轴载荷作用下的理想晶粒模型的晶界上一系列孔洞。若以下假设成立，孔洞的长大可以通过包含孔洞的板层的体积变化来度量：

(i) 蠕变变形中材料不可压缩，且其总体积保持不变；

(ii) 球形孔洞只改变体积而不改变形状；

(iii) 包含孔洞的板层的宽度要远大于其厚度；

(iv) 晶界滑移使得两边的晶粒刚性位移和板层体积变化相适应；

(v) 静水压力对于无孔洞时的蠕变变形没有影响；

(vi) 晶界孔洞长大仅受幂律蠕变控制。

图5所示：包含一个孔洞的圆柱体的为晶粒尺寸，为孔洞的半径，为孔洞间的距离，是进入边界计算的距离，是轴向应力，是附加的三轴应力。那么，晶界上孔洞的面积分数可表示为：

（1）

单轴应力作用下，圆柱体的体积变化率与的关联如下：

（2）

其中，代表圆柱体的体积，是不含孔洞时的稳态蠕变速率。另一方面，圆柱体的体积变化还可以由蠕变速率来定义：

（3）

其中，为包含一个孔洞的圆柱体的轴向应变速率，为包含一个孔洞的圆柱体的径向应变速率。

步骤（2），利用能量原理，由幂律蠕变控制孔洞长大理论，通过New-Raphson方法，获得单轴及多轴应力状态参数；

利用能量原理，Cocks和Ashby给出了多轴应力状态下的上边界具有如下形式：

（4）

式中G的定义为：

（5）

其中，，，且为应力三轴度，为确定的材料常数。

合并式（2）、式（3）和式（4），可以得到描述幂律蠕变控制孔洞长大速率的复杂数学表达式：

（6）

应当注意的是，从式（6）中无法直接获得孔洞长大速率和应力三轴度之间的关系，为了使结果更方便实用，Cocks和Ashby给出了一个拟合的半经验公式，如下所示：

（7）

其中，多轴应力状态参数定义为：

（8）

然而，公式(8) 在材料处于压缩状态，即应力三轴度时，多轴应力状态参数为负值或出现数值奇异现象，这显然与实际情况不符。此外，随应力三轴度的增加而逐渐减小，当时，近乎为0，这与实际观察到的空洞变化理论不相符。为了解决此问题，本发明中提出了另外一种近似公式，如下式所示：

（9）

其中，与为与材料相关的常数，在简单拉伸的情况下，即时，，，式（9）退化为下式：

（10）

其中，为单轴应力状态参数。

步骤（3），对微观孔洞进行积分获得材料恒载下的失效时间，进而表示单轴及多轴应力作用下的蠕变失效应变，并获得多轴蠕变延性因子；

为了获得多轴加载情况下的蠕变失效应变，将式（7）在如下所示的上下限进行积分：

（11）

其中，为初始孔洞面积分数，为孔洞合并发生时的面积分数，为多轴应力状态下从初始到孔洞合并发生所需的时间。恒载荷下积分所得结果便是失效时间：

（12）

多轴应力作用下的蠕变失效应变可通过下式得到：

（13）

用同样的方法可以得到单轴应力状态下的蠕变失效应变：

（14）

其中，为单轴应力状态下从初始到孔洞合并发生所需的时间。

可定义新的多轴蠕变延性因子MCDF：

（15）

对于具体的材料，在某一应力状态下，为确定的材料常数，因而公式（15）可以简化为下式：

（16）

其中，为与材料相关的常数。

步骤（4），对多轴蠕变延性因子进行参数拟合，获得不同应力状态下的多轴蠕变失效应变拟合参数，从而得到多轴蠕变延性因子预测方程；

对公式（16）中的参数进行拟合：

公式（16）中，和均为与材料有关的常数，可通过缺口蠕变试验获得不同应力三轴度下的蠕变失效应变，然后利用公式（16）进行拟合，获得参数和，从而得到多轴蠕变延性因子预测方程。

步骤（5），利用有限元软件预测材料在多轴应力状态下的蠕变失效应变及其寿命；

根据步骤（4）中建立的多轴蠕变延性因子预测方程，结合蠕变-损伤本构方程，利用Fortran语言，编制子程序并嵌入到有限元软件ABAQUS中，从而实现材料在多轴应力状态下的蠕变失效应变及寿命的预测。

表1 不同转化模型获得的多轴蠕变失效应变与试验值误差分析

表2 圆柱缺口拉伸试样缺口敏度比的数据变化表

表3 不同蠕变损伤本构模型及拉伸载荷下，各缺口敏度比试样的蠕变寿命预测值及试验值的对比数据表

预测9Cr-1Mo合金在600°C下的多轴蠕变失效应变。图6所示为在不同的应力三轴度下不同模型预测的多轴蠕变失效应变与试验值的比较曲线。本发明中所提的模型中的参数通过试验数据拟合获得，拟合过程中蠕变指数，单轴蠕变失效应变，获得的材料常数，。从图3中可以看出，本发明中提出的转化模型有效的避免了因材料中应力三轴度时出现的比值过小的问题。表1所示为利用Cocks and Ashby (C-A)模型、Wen and Tu (W-T)模型以及本发明中提出的转化模型获得的多轴蠕变失效应变与试验值得误差分析。从表1中可以看出，本发明中提出的模型获得多轴蠕变失效应变的误差最小。从图6及表1中的数据表明了本发明中提出的计算多轴蠕变失效应变的合理性和可靠性。改变圆柱缺口拉伸试样的缺口敏度比，利用有限元软件ABAQUS模拟预测在拉伸载荷为130MPa、150MPa、170MPa、210MPa四种情况下，圆柱缺口拉伸试样的蠕变寿命。缺口蠕变试样的几何形状如图7所示，敏度比的数据变化如表2所示。计算过程中分别采用K-R蠕变损伤本构模型、W-T蠕变损伤本构模型，对不同拉伸载荷下，不同缺口敏度比试样的蠕变寿命进行模拟预测。试验值及模拟值的对比数据如表3所示。从表中可以看出，与试验蠕变寿命对比，本发明提出的拟合修正模型估计的误差是最小的，然后是K-R蠕变损伤模型，Wen-Tu蠕变损伤模型计算结果误差最大。除了个别数据外，所提出的模型预测的误差控制在50%的范围内。从表1和表3的结果反映，本发明所提出的蠕变本构模型可以合理地预测高温材料的蠕变和失效行为，因此，本发明提出的蠕变失效应变预测方法可以准确计算材料的多轴蠕变失效应变及结构的寿命。

改变圆柱缺口拉伸试样的缺口敏度比，利用有限元软件ABAQUS模拟预测在拉伸载荷为130MPa、150MPa、170MPa、210MPa四种情况下，圆柱缺口拉伸试样的蠕变寿命。缺口蠕变试样的几何形状如图7所示，敏度比的数据变化如表2所示。计算过程中分别采用K-R蠕变损伤本构模型、W-T蠕变损伤本构模型，对不同拉伸载荷下，不同缺口敏度比试样的蠕变寿命进行模拟预测。试验值及模拟值的对比数据如表3所示。从表中可以看出，与试验蠕变寿命对比，本发明提出的拟合修正模型估计的误差是最小的，然后是K-R蠕变损伤模型，Wen-Tu蠕变损伤模型计算结果误差最大。除了个别数据外，所提出的模型预测的误差控制在50%的范围内。从表1和表3的结果反映，本发明所提出的蠕变本构模型可以合理地预测高温材料的蠕变和失效行为，因此，本发明提出的蠕变失效应变预测方法可以准确计算材料的多轴蠕变失效应变及结构的寿命。

以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

Claims (6)

1.一种材料多轴蠕变失效应变预测方法，其特征在于：包括如下步骤：

步骤（1），基于应变损伤准则，根据微观孔洞长大机制，获得材料蠕变速率和应变速率之间的关系；

步骤（2），由幂律蠕变控制孔洞长大理论，获得单轴及多轴应力状态参数；

步骤（3），对微观孔洞进行积分获得材料恒载下的失效时间，进而表示单轴及多轴应力作用下的蠕变失效应变，并获得多轴蠕变延性因子；

步骤（4），对多轴蠕变延性因子进行参数拟合，获得不同应力状态下的多轴蠕变失效应变拟合参数，从而得到多轴蠕变延性因子预测方程；

步骤（5），利用有限元软件预测材料在多轴应力状态下的蠕变失效应变及其寿命；

步骤（2）中所述的单轴及多轴应力状态参数利用能量原理，并通过New-Raphson方法获得；

步骤（2）中所述的单轴及多轴应力状态参数的计算公式为：

，

，

其中，为多轴应力状态参数，为单轴应力状态参数，与为与材料相关的常数，为确定的材料常数，为应力三轴度。

2.根据权利要求1所述的材料多轴蠕变失效应变预测方法，其特征在于：步骤（1）中所述的微观孔洞长大机制在裂尖区域为粘塑性孔洞长的机制。

3.根据权利要求1所述的材料多轴蠕变失效应变预测方法，其特征在于：步骤（1）中所述的材料蠕变速率和应变速率之间的关系为：

，

其中，为包含一个孔洞的圆柱体的轴向应变速率，为包含一个孔洞的圆柱体的径向应变速率，为圆柱体不包含孔洞时的稳态蠕变速率。

4.根据权利要求1所述的材料多轴蠕变失效应变预测方法，其特征在于：步骤（3）中所述的多轴蠕变延性因子的计算公式如下：

，

其中，为单轴应力作用下的蠕变失效应变，为多轴应力作用下的蠕变失效应变，和均为与材料相关的常数，为确定的材料常数，为应力三轴度。

5.根据权利要求4所述的材料多轴蠕变失效应变预测方法，其特征在于：利用缺口蠕变试验获得不同应力状态下的蠕变失效应变拟合参数和，进而得到步骤（4）中所述的多轴蠕变延性因子预测方程。

6.根据权利要求1所述的材料多轴蠕变失效应变预测方法，其特征在于：步骤（4）中所述的多轴蠕变延性因子预测方程接合蠕变-损伤本构方程，并通过Fortran语言嵌入到有限元软件中。