发明内容

本发明的发明目的在于：针对目前材料多应变率本构模型硬化曲线建立过程中所存在的问题，提出获得多应变率下材料高精度硬化模型参数的方法，主要解决曲线建立过程中不能同时兼顾多应变率下材料性能和最终应用效果较差的问题，采用多目标多参数的优化方法对多应变率下的材料试样进行材料硬化模型曲线建立，以克服现有技术所存在的不足。

本发明采用的技术方案如下：一种获得多应变率下材料高精度硬化模型参数的方法，包括以下步骤：

S1、对材料进行高速拉伸试验，获得材料在多应变率单向拉伸下的工程应力-工程应变曲线；

S2、计算材料试样多应变率下的真应力塑性应变曲线；

S3、对多应变率下真应力塑性应变曲线进行拟合外延，得到各应变率所对应的外延应力应变曲线；

S4、将各应变率所对应的外延应力应变曲线组合成外延应力应变表曲线表，曲线所对应值为曲线所对应应变率；

S5、通过赋予各应变率下外延应力应变曲线加权系数α不同的值，调整外延应力应变表曲线表的形状，α的取值范围为0-1；

S6、对各应变率的材料试样建立数值模型，统一调用S5所得的外延应力应变表曲线表，在有限元软件中进行仿真计算，对比试验及仿真结果中的力-变形曲线；

S7、返回S5，优化加权系数α的值，直到S6中所有材料试样数值模型对标结果满足误差5％以内的要求，最终得到多应变率高精度硬化模型参数。

在本发明的方法中，力除以试样平行段截面积获得材料试样的工程应力，变形除以标距获得材料试样的工程应变。进一步，在S1中，根据国标GB/T228.1-2010金属材料拉伸试验第1部分的室温试验方法，找出试验所得工程应力-应变中弹性段部分，获得材料试样的屈服强度、抗拉强度、弹性模量数据。

在本发明的方法中，计算材料的工程应力和工程应变，剔除屈服强度之前、颈缩点之后的数据，并通过公式转换得到材料试样的真应力塑性应变曲线。材料的真应力和真应变及塑性应变的计算采用如下公式：

真应力计算：

σ_T＝σ*(1+ε) 公式(1)

真应变计算：

ε_T＝ln(1+ε) 公式(2)

塑性应变计算：

ε_pl＝ln(1+ε_T-σ/E) 公式(3)

公式(1)和公式(2)中，σ和ε分别为工程应力应变。公式(3)中，ε_pl为塑性应变，ε_T为真应变，E为弹性模量。

在本发明中，由于试验所获得并处理的真应力塑性应变曲线仅为缩颈前数据。对于缩颈后数据，因实际截面缩小导致所测得应力失真，因此需采用硬化模型对处理后真应力塑性应变曲线数据进行拟合外推。常用的硬化模型分为饱和硬化模型和非饱和硬化模型，具体公式如下所示。为保证曲线有较大的调整范围，选用饱和硬化模型和非饱和硬化模型各一种采用加权系数进行混合得到的混合硬化模型，从而获得较大的调整空间。

非饱和型硬化模型:

Hollmon本构方程：

简化J-C本构方程：

Swift本构方程：

Ghosh本构方程：

Voce++本构方程：

饱和型硬化模型:

Hockett-Sherby本构方程：

Voce本构方程：

上述公式中，a、b、c、d为未知参数，需拟合得到。

在本发明中，选用Voce++和Hockett-Sherby本构方程(当然，也可以选择其他方程组合，不同方程之间可以自由组合，此处仅以这两个方程为例)拟合得到外延后的应力应变曲线，Voce++和Hockett-Sherby本构方程结合后的方程如公式(11)所示：

Voce++-Hockett-Sherby本构方程：

式中：σ为真应力，α为加权系数，取值为(0-1)，a₅、a₆、b₅、b₆、c₅、c₅、d₅、d₆为未知参数，需拟合得到，ε_pl为塑性应变。

进一步，所述多应变率用试样平行段的应变率进行描述，其中，多应变率的范围为0.1/s-1000/s，所述多应变率对应的试样为高速拉伸试样。

作为优选，所述有限元软件为LS-DYNA仿真分析软件。

综上所述，由于采用了上述技术方案，本发明的有益效果是：

1、本发明先通过材料高速拉伸试验并通过数据处理获得多应变率下材料在颈缩点前的真应力塑性应变曲线，基于该曲线，采用硬化模型公式进行拟合，得到初始硬化模型参数表，然后基于不同应变率下的材料高速拉伸试验，建立相应加载条件的有限元模型，同时对标各应变率试验中的力-变形曲线，并优化硬化模型参数，优化过程中不断调整硬化模型参数，直到所有条件下试样的对标结果满足要求，从而得到高精度的多应变率材料硬化模型参数；

2、本发明将不同应变率下的材料试样进行同时对标，优化迭代得到高精度的多应变率硬化模型参数，解决了现有硬化模型参数建立过程中不能同时兼顾不同应变率下材料性能表征的问题，在实际试样对标中的应用效果良好，误差在5％以内，克服现有方法所存在的不足。