材料工程  2016, Vol. 44 Issue (5): 15-21   PDF    
http://dx.doi.org/10.11868/j.issn.1001-4381.2016.05.003
0

文章信息

张施琦, 冯定, 张跃, 洪继要
ZHANG Shi-qi, FENG Ding, ZHANG Yue, HONG Ji-yao
新型超高强度热冲压用钢的热变形行为及本构关系
Hot Deformation Behavior and Constitutive Model of Advanced Ultra-high Strength Hot Stamping Steel
材料工程, 2016, 44(5): 15-21
Journal of Materials Engineering, 2016, 44(5): 15-21.
http://dx.doi.org/10.11868/j.issn.1001-4381.2016.05.003

文章历史

收稿日期: 2015-10-26
修订日期: 2016-03-24
新型超高强度热冲压用钢的热变形行为及本构关系
张施琦1,2, 冯定1,2 , 张跃3, 洪继要4    
1. 非常规油气湖北省协同创新中心, 武汉 430100 ;
2. 长江大学 机械工程学院, 湖北 荆州 434023 ;
3. 北京科技大学 材料科学与工程学院, 北京 100083 ;
4. 汽车用钢开发与应用技术国家重点实验室(宝钢), 上海 201900
摘要: 利用 Gleeble-1500D热模拟机对新型超高强度热冲压用钢22MnB5Nb进行等温单向拉伸实验,研究了其在变形温度为650~950℃,应变速率为0.1,1.0,10s-1下的热变形行为,并采用3种本构分析方法,即基于传统拟合回归方法的Arrhenius 型、考虑材料常数应变补偿的Arrhenius 型和本工作新提出的基于Quasi-Newton BFGS算法的Arrhenius 型本构方程来描述22MnB5Nb钢的热变形行为。结果表明:22MnB5Nb钢表现出典型的加工硬化和动态回复软化行为,变形温度与应变速率均对其流变应力有较大影响;3种方程均可以准确预测实验钢的峰值流变应力,其中,Quasi-Newton BFGS算法具有可一次性求解所有材料参数、求解步骤简单和预测精度最高(R=0.99578,Re=11.03MPa,E=2.48%)的特点,考虑材料常数应变补偿的Arrhenius 型本构方程预测精度相对较低,但能直接预测不同变形条件下的流变应力曲线且可以较好地预测变形过程中的加工硬化效应、动态回复软化效应和应变速率强化效应。
关键词: 热冲压用钢    热变形行为    本构模型    Nb   
Hot Deformation Behavior and Constitutive Model of Advanced Ultra-high Strength Hot Stamping Steel
ZHANG Shi-qi1,2, FENG Ding1,2 , ZHANG Yue3, HONG Ji-yao4    
1. Hubei Cooperative Innovation Center of Unconventional Oil and Gas, Wuhan 430100, China ;
2. School of Mechanical Engineering, Yangtze University, Jingzhou 434023, Hubei, China ;
3. School of Materials Science and Engineering, University of Science and Technology Beijing, Beijing 100083, China ;
4. State Key Laboratory of Development and Application Technology of Automotive Steels(Baosteel), Shanghai 201900, China
Abstract: The hot deformation behavior of the advanced ultra-high strength hot stamping steel 22MnB5Nb was studied through the isothermal uniaxial tensile tests at 650-950℃ and strain rates of 0.1, 1.0s-1 and 10s-1 by Gleeble 1500D system. The conventional Arrhenius-type hyperbolic sine equation, the Arrhenius-type model considering the material constant strain compensation and the new Arrhenius-type model based on Quasi-Newton BFGS algorithm were established to describe the high-temperature deformation behavior of 22MnB5Nb. The results indicate that 22MnB5Nb steel shows typical work hardening and dynamic recovery softening behavior during hot tensile. And the strain rate and deformation temperature have significant effects on the flow stress. The peak flow stress values predicted by these models are highly consistent with the experimental values, and the Quasi-Newton BFGS algorithm can solve all the material parameters in one time and it is simpler in calculation process and has the highest accurate(R=0.99578, Re=11.03MPa, E=2.48%), while, the Arrhenius-type model considering material constant strain compensation with lower accuracy, but can directly predicts not only the flow stress curve under different deformation conditions, but also the work hardening behavior, the dynamic recovery behavior and the strain rate strengthening effect of the experimental steel during the deformation process.
Key words: hot stamping steel    hot deformation behavior    constitutive model    Nb   

先进高强钢在汽车结构件中的应用是实现汽车轻量化和提高车身碰撞安全性的有效途径之一。高强度钢板热冲压成形技术克服了传统冷冲压过程中高强度钢板易开裂、回弹严重且成形困难的缺点,所制成的零件具有较高的成形强度和成形精度,在汽车高强度零部件的制造中已得到广泛应用[1, 2, 3]。在热冲压过程中,钢板被预先加热到奥氏体化温度以上,然后在模具中迅速冲压成形同时进行淬火,其间钢板经历复杂的温度、应力、组织变化,这些变化都会影响最终成形性能;因而,研究热冲压钢板的热变形行为并建立本构方程对热冲压工艺优化及组织性能控制具有十分重要的意义。

22MnB5钢是目前热冲压工艺中最常用的钢材,对于其热变形行为和本构模型的研究也已大量开展,相继建立了Johnson-Cook[1, 2]、Voce-Kocks[2]、Arrhenius[3]、Norton-Hoff型的本构模型[4, 5, 6]。Naderi等[2]基于等温单轴压缩实验数据建立了Voce-Kocks模型和Molinari-Ravichandran模型来描述22MnB5钢的热变形行为。Li等[3]较为系统地研究了热冲压硼钢在各种微观组织状态下的热变形行为,并基于Arrhenius和Johnson-Cook方程建立其本构模型。最近,Zhang等[7]在此基础上开发出了1800MPa级新型热冲压用钢22MnB5Nb,由于其相比传统热冲压钢具备更好的强韧性、更优异的耐延迟开裂性能等特性,使其不仅在汽车制造中广泛应用,而且在页岩气开采装备、航空航天装备等的承力结构件上也具有广阔的应用前景。而关于这种新型热冲压用钢的热变形行为及本构关系的研究鲜有报道,这在一定程度上制约了新型热冲压钢的研发和应用,因此,本工作以新型超高强度热冲压用钢22MnB5Nb为研究对象,采用Gleeble 1500D 热模拟试验机研究其热变形行为,并建立反映该材料流动特性的本构模型。

1 实验材料与方法

实验钢采用50kg真空感应炉冶炼,其化学成分见表 1。热拉伸试样的尺寸如图 1所示,标距为60mm。等温单向热拉伸实验在Gleeble-1500D 热模拟试验机上进行,具体流程为:将试样以5℃/s 的速率加热到1000℃,保温300s使材料充分奥氏体化,然后以40℃/s的速率冷却至变形温度后保温10s,使试件的温度稳定,再进行等温拉伸。变形温度为600,700,800,900,950℃,变形速率为0.1,1.0,10s-1,变形结束后进行空冷。

表 1 新型热冲压用钢22MnB5Nb化学成分(质量分数/%) Table 1 Chemical compositions of 22MnB5Nb steel (mass fraction/%)
CSiMnSPAlTiCrBNbFe
0.23000.30001.2000≤0.0100≤0.01500.03000.02500.15000.00250.0500Bal
图1 热拉伸试样尺寸 Fig.1 Dimension of specimen for hot tensile test
2 结果与分析 2.1 真应力-真应变曲线

图 2为22MnB5Nb钢在不同变形温度与应变速率下的真应力-真应变曲线,可见,实验钢的流变应力主要与应变速率、变形温度及应变量有关。在同一应变速率和应变量一定的条件下,其流变应力随变形温度升高而降低(图 2(a)),在同一变形温度和应变量一定的条件下,其流变应力随应变速率增加而增大(图 2(b))。同时,实验钢的流变应力随着应变量的增加而逐渐增大到某一峰值,此后便趋于稳定,未表现出明显下降,表现为典型的动态回复型变形行为。这是由于在变形初始阶段,实验钢因发生加工硬化而使得流变应力不断提高;随着变形的进行,位错密度不断增加,储存能不断增大,回复的驱动力也变大,由此引起的动态回复软化作用逐渐抵消部分加工硬化效果,从而在真应力-真应变曲线上表现出真应力上升变缓的趋势,最终在达到峰值应力后,实验钢中的动态回复软化作用与加工硬化作用处于动态平衡而表现为稳态流变特征[3, 6]

图2 22MnB5Nb钢的真应力-真应变曲线 (a)温度;(b)应变速率 Fig.2 True stress-strain curves of steel 22MnB5Nb at different temperatures (a) and different strain rates (b)
2.2 基于传统拟合回归方法的热变形本构方程

Arrhenius 本构方程是目前使用最广泛的材料高温热变形本构方程之一,常用于描述材料在高温下的流变应力与应变速率、变形温度之间的关系[8, 9, 10]。它包括指数方程、幂指数方程和双曲正弦方程3种形式[11, 12, 13]

(1)
(2)
(3)
(4)

式中:σ为流变应力,MPa;为应变速率,s-1R为气体常数,8.3145J·mol-1·K-1T为温度,K;Q为热变形激活能,J·mol-1A (A1,A2),αn (n1),β为材料常数,α=β/n1。在传统的本构方程中,Zener-Hollomon参数Z被用于表示应变速率和变形温度对材料热变形行为的影响(式(1)),幂指数形式的本构方程适用于流变应力较低时(ασ<0.8)(式(2)),指数形式的本构方程适用于流变应力较高时(ασ>1.2)(式(3)),而双曲正弦本构方程适用于较宽应力范围,且能更好地描述参数Z与应力的关系[14, 15]。上述3个方程式(式(2)~(4))的描述是不完整的,这是由于式中并未指明流变应力对应的应变量,但在目前研究中广泛选用峰值应力来建立本构方程,本工作也是如此。

对式(2),(3)两边分别取自然对数,可得:

(5)
(6)

当变形温度为恒定值时,利用偏微分法可求得

(7)
(8)

对ln-lnσP和lnP进行线性回归,如图 3所示,由斜率可以得到n1β值。对不同温度下n1β值取平均值,可得n1=9.11911,β=0.03155,α=β/n1=0.00346。

图3 峰值应力与应变速率的关系曲线 (a)ln-lnσP;(b)lnP Fig.3 Relationships between peak stress and strain rate (a)ln-lnσP;(b)lnP

对式(4)两边取对数可得

(9)

当分别假设变形温度和应变速率恒定时,利用偏微分可得式(10),(11):

(10)
(11)

根据ln[sinh(ασ)]-ln和ln[sinh(ασ)]-1/T关系,如图 4所示,可以求得n的平均值为6.7156,Q值为182.683kJ/mol,A值为2.609×108。将求解值代入式(1),(4)中即得到22MnB5Nb钢的双曲正弦本构方程为:

(12)
图4 峰值应力与应变速率以及温度的关系曲线 (a)ln[sinh(ασ)]-ln;(b)ln[sinh(ασ)]-1/T Fig.4 Relationships between peak stress,strain rate and deformation temperature (a)ln[sinh(ασ)]-ln;(b)ln[sinh(ασ)]-1/T
2.3 基于Quasi-Newton BFGS算法的新型热变形本构方程

目前报道的双曲正弦本构方程均采用上文所述传统方法进行求解[11, 12, 13, 16],即通过线性回归求得局部范围内的参数n1β值,然后对α进行求解,此时获得的α并非全局解,且参数nQ也是在各变形参数下的近似平均值,并非最优解;同时,利用传统方法求解不能一次求得所有参数,在求解过程中须进行多次拟合,步骤较复杂。

Quasi-Newton BFGS算法是求解非线性优化问题最有效的方法之一,在数学问题求解优化中得到了广泛应用[17, 18],但鲜有研究者将其应用于本构方程求解。为此,本工作提出一种基于Quasi-Newton BFGS算法的新型求解方法,即利用Matlab软件,采用最小二乘法中的Quasi-Newton BFGS算法,对双曲正弦本构方程进行一次性拟合求得各参数的全局优化解。此方法求解简单且能更好地考虑各参数与全局数据之间的关系。

鉴于Arrhenius双曲正弦本构方程的高度非线性特征,首先对式(4)变形以提高其拟合精度:

(13)

设:y=σPx1=lnx2=1/Tk1=1/nk2=Q/nRk3=lnA/nk4=1/α,则有:

(14)

在Matlab软件中导入实验数据,利用Quasi-Newton BFGS算法对式(14)进行拟合求解。在求解过程中需要设定合适初值以便于计算收敛,已有研究表明各类钢铁材料的Arrhenius本构方程的αn等系数均处于同一数量级[9, 10, 11],故任意选定某一钢种的系数来赋初值均能够得到很好的收敛且求解结果几乎无差别。最终求解得到:k1=0.1307,k2=2.9481,k3=2.627,k4=293.2551,α=0.00341,n=7.653,Q=187.582,lnA=20.1081。将求解值代入式(1),(4)中,可得:

(15)
2.4 考虑材料常数应变补偿的Arrhenius本构方程

传统计算方法只考虑峰值应变对应的峰值应力,而并未将其他应变量考虑在内,这在一定程度上限制了其应用。为此,有学者在传统Arrhenius型本构方程基础上开发出了考虑材料常数应变补偿的Arrhenius型本构方程,并得到广泛采用[11, 13, 16]。在本工作中,利用上述求解峰值应力材料常数的方法,分别对不同应变量(0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50)的材料常数值进行求解,并利用式(16)~(19)对求解数值和应变量进行5次多项式拟合分析,得到拟合系数如表 2所示,拟合结果如图 5所示。

(16)
(17)
(18)
(19)
表 2 多项式拟合系数值 Table 2 Coefficients of the polynomial for α,lnA,n and Q
αlnAnQ
α1=0.00608A1=27.60187n1=8.82911Q1=264691.7
α2=-0.02916A2=-132.84704n2=-37.69934Q2=-1343360.0
α3=0.15040A3=1433.02609n3=422.57911Q3=14347100.0
α4=-0.44154A4=-7080.36721n4=-2113.80460Q4=-70293100.0
α5=0.69953A5=15037.02974n5=4524.67240Q5=148624000.0
α6=-0.45202A6=-11339.80718n6=-3437.14550Q6=-111832000.0
图5 多项式拟合实验钢的材料常数与真应变的关系(a)α;(b)n;(c)lnA;(d)Q Fig.5 Relationship between α (a),n(b),lnA(c),Q(d) and ε
2.5 本构方程的预测精度分析

本工作建立的3种本构方程均可用于材料的热变形本构特征的研究。利用基于传统拟合回归方法的本构方程和基于Quasi-Newton BFGS算法的本构方程可以求得某一应变量所对应的应力(如峰值应力、稳态应力等,在本工作中为峰值应力)和变形参数间的关系。而利用考虑材料常数应变补偿的Arrhenius本构方程则能预测不同变形条件下的流变应力曲线,得到任意应变量对应的应力和变形参数间的关系[13]。故在本工作中选用3种本构方程都能计算求得的不同变形参数下的理论峰值应力来对本构方程的精度进行分析。

为定量评价所建立方程的精度,引入相关系数(R)、均方根误差(Re)和平均相对误差(E)参数进行评价[10],如图 6所示。由图 6可知,3种本构方程所预测的流变应力与实验值均存在较好的相关性,其中,基于Quasi-Newton BFGS算法的本构方程具有最高预测精度(R=0.99578,Re =11.03MPa,E=2.48%),传统双曲正弦形式的Arrhenius 型本构方程其次(R=0.99543,Re =11.26MPa,E=2.60%),考虑材料常数应变补偿的Arrhenius本构方程精度最低(R=0.99517,Re =11.52MPa,E=2.66%)。但应当注意的是,考虑材料常数应变补偿的Arrhenius本构方程由于考虑了应变量的影响且能直接得到应力-应变曲线而具有独特优势,在实际应用中(如热加工工艺制定、有限元模拟等)需要获取多个应变量下应力与参数关系时,则优先选用本方程。图 7为采用考虑材料常数应变补偿的Arrhenius本构方程预测得到的曲线与实验测得的真应力-真应变曲线对比,可见在不同变形条件下两者均有较好吻合,且该方程可以较好地预测变形过程中的加工硬化效应、动态回复软化效应和应变速率强化效应。

图6 实验钢不同本构方程预测的峰值应力值与实验值的相关性 (a)基于传统拟合回归方法的本构方程;(b)基于Quasi-Newton BFGS算法的本构方程;(c)考虑应变补偿的本构方程 Fig.6 Correlation between measured peak stress and predicted data by the three models (a)the conventional hyperbolic sine equation;(b)the model based on Quasi-Newton BFGS algorithm;(c)the Arrhenius-type model with the strain compensation
图7 考虑应变补偿的本构方程的预测应力-应变曲线与实验值的比较 (a)0.1s-1;(b)1.0s-1;(c)10s-1 Fig.7 Comparison between predicted and experimental stress-strain curves by the Arrhenius-type model with the strain compensation (a)0.1s-1;(b)1.0s-1;(c)10s-1
3 结论

(1)22MnB5Nb钢在变形温度为650~950℃,应变速率为0.1,1.0,10s-1的条件下的热变形行为表现为典型的加工硬化和动态回复型,变形温度与应变速率均对其流变应力有较大影响。

(2)建立了实验钢的传统双曲正弦形式的Arrhenius 型和考虑材料常数应变补偿的Arrhenius 型本构方程,并提出了基于Quasi-Newton BFGS算法的Arrhenius 型本构方程求解方法,所提出的新算法具有可一次性求解所有材料参数和求解步骤简单的特点。

(3)3种本构方程均可以对实验钢的峰值流变应力给出准确的预测。其中,基于Quasi-Newton BFGS算法的本构方程具有最高预测精度。考虑材料常数应变补偿的Arrhenius本构方程精度相对较低,但其能直接预测不同变形条件下的流变应力曲线。

参考文献(References)
[1] KARBASIAN H, TEKKAYA A E. A review on hot stamping[J]. Journal of Materials Processing Technology,2010, 210 (15) : 2103 –2118.
[2] NADERI M, DURRENBERGER L, MOLINARI A, et al. Constitutive relationships for 22MnB5 boron steel deformed isothermally at high temperatures[J]. Materials Science & Engineering: A,2008, 478 (1-2) : 130 –139.
[3] LI H, HE L, ZHAO G, et al. Constitutive relationships of hot stamping boron steel B1500HS based on the modified Arrhenius and Johnson-Cook model[J]. Materials Science & Engineering: A,2013, 580 : 330 –348.
[4] NIKRAVESH M, NADERI M, AKBARI G H. Influence of hot plastic deformation and cooling rate on martensite and bainite start temperatures in 22MnB5 steel[J]. Materials Science & Engineering: A,2012, 540 : 24 –29.
[5] NADERI M, KETABCHI M, ABBASI M, et al. Analysis of microstructure and mechanical properties of different high strength carbon steels after hot stamping[J]. Journal of Materials Processing Technology,2011, 211 (6) : 1117 –1125.
[6] ZHOU J, WANG B, HUANG M. Two constitutive descriptions of boron steel 22MnB5 at high temperature[J]. Materials & Design,2014, 63 (2) : 738 –748.
[7] ZHANG S, HUANG Y, SUN B, et al. Effect of Nb on hydrogen-induced delayed fracture in high strength hot stamping steels[J]. Materials Science & Engineering: A,2015, 626 : 136 –143.
[8] SAMANTARAY D, MANDAL S, BHADURI A K. Constitutive analysis to predict high-temperature flow stress in modified 9Cr-1Mo (P91) steel[J]. Materials & Design,2010, 31 (2) : 981 –984.
[9] XIAO X, LIU G Q, HU B F, et al. A comparative study on Arrhenius-type constitutive equations and artificial neural network model to predict high-temperature deformation behavior in 12Cr3WV steel[J]. Computational Materials Science,2012, 62 : 227 –234.
[10] LI H Y, HU J D, WEI D D, et al. Artificial neural network and constitutive equations to predict the hot deformation behavior of modified 2.25Cr-1Mo steel[J]. Materials & Design,2012, 42 : 192 –197.
[11] 王进, 褚忠, 张琦. 38MnVS6非调质钢两种高温本构模型的对比[J]. 材料工程,2014 (2) : 81 –86. WANG Jin, CHU Zhong, ZHANG Qi. A comparative study of two high-temperature constitutive models of 38MnVS6 microalloyed forging steel[J]. Journal of Materials Engineering,2014 (2) : 81 –86.
[12] ZENG Z, ZHANG Y, JONSSON S. Microstructure and texture evolution of commercial pure titanium deformed at elevated temperatures[J]. Materials Science & Engineering: A,2009 : 83 –90.
[13] 魏海莲, 刘国权, 肖翔, 等. 表观的和基于物理的35Mn2钢奥氏体热变形本构分析[J]. 金属学报,2013, 49 (6) : 731 –738. WEI Hai-lian, LIU Guo-quan, XIAO Xiang, et al. Apparent and physically based constitutive analyses for hot deformation of austenite in 35Mn2 steel[J]. Acta Metallurgica Sinica,2013, 49 (6) : 731 –738.
[14] MIRZADEH H, CABRERA J M, PRADO J M, et al. Hot deformation behavior of a medium carbon microalloyed steel[J]. Materials Science & Engineering: A,2011, 528 : 3876 –3882.
[15] MEYSAMI M, MOUSAVI S. Study on the behavior of medium carbon vanadium microalloyed steel by hot compression test[J]. Materials Science & Engineering: A,2011, 528 : 3049 –3055.
[16] 郑漫庆, 王高潮, 喻淼真, 等. 应变速率循环法构建TC4-DT钛合金本构方程[J]. 材料工程,2014 (8) : 32 –35. ZHENG Man-qing, WANG Gao-chao, YU Miao-zhen, et al. Superplastic constitutive relationship of TC4-DT titanium alloy with strain rate circulation method[J]. Journal of Materials Engineering,2014 (8) : 32 –35.
[17] ZHANG Y, MU B, ZHENG H. Link between and comparison and combination of Zhang neural network and Quasi-Newton BFGS method for time-varying quadratic minimization[J]. IEEE Trans Syst Man Cybern B Cybern,2012, 43 (2) : 490 –503.
[18] LEITHEAD W E, ZHANG Y N. O(N2)-operation approximation of covariance matrix inverse in Gaussian process regression based on Quasi-Newton BFGS method[J]. Communications in Statistics Simulation & Computation,2007, 36 (2) : 367 –380.