Comparative Analysis of 3 Constitutive Models for 93 Tungsten Alloy
-
摘要: 在温度范围296~1 273 K、应变率范围0.000 5~6 000 s-1内, 对93钨合金材料(93%W-4.9%Ni-2.1%Fe)的力学行为进行了系统研究。实验结果表明:钨合金的塑性流动应力对温度和应变率非常敏感,它随温度的降低和应变率的增加而提高;真实塑性应变达到0.6时未出现剪切破坏。针对实验得到的应力-应变关系,分别进行了两种唯象本构模型(JC本构模型和KHL本构模型)和一种基于物理概念的本构模型(PB本构模型)的流动应力拟合。通过误差分析和应变率阶跃实验,对这3种本构模型的预测效果的精度进行了分析和评价。Abstract: In this research, we carried out a systematic study of the mechanical behaviors of the 93 tungsten alloy (93%W-4.9%Ni-2.1%Fe) with a temperature range of 296-1 273 K and a strain rate range of 0.000 5-6 000 s-1.The results demonstrate that the flow stress is highly dependent on the temperature and the strain rate:the flow stress increases with the decrease of temperature and the increase of strain rate, and the shear failure does not occur when the true plastic strain is up to 0.6.Then we introduced and established 3 constitutive models, including two phenomenological constitutive models, JC model and KHL model and one physically-based constitutive model, to fit the stress-strain relationship through a procedure of regression analysis and constrained optimization.Finally, we evaluated the fitted results of the 3 constitutive models through a fitting error analysis and a specific strain-rate jump experiment.
-
Key words:
- tungsten alloy /
- plastic flow stress /
- constitutive model
-
1. 引言
钨合金材料因其高密度、高强度特性被广泛应用于国防军事领域,主要是制造侵彻动能弹和代替铀制造穿甲弹。在动能侵彻过程中,材料的塑性变形通常包含应变、应变率和温度耦合效应,从而使预测材料的塑性流动应力变得极具挑战性。范景莲等[1]综述了钨合金的力学性能研究,得出钨合金对应变率敏感而对绝热剪切不敏感的结论。魏志刚等[2]通过实验研究了在动态载荷作用下预扭钨合金的微观结构对绝热剪切带位置的影响,对不同形状的钨合金试样绝热剪切带形成机理进行了讨论。Kim等[3]研究了钨含量为89%的钨合金试样在两端装有限制环时,材料的绝热剪切带产生的情况。Yadav等[4]研究了两种不同的钨合金在103~105 s-1应变率范围内的压缩和剪切力学特性,并给出了其中一种钨合金的唯象本构模型。Coates等[5]测量了90%W-7%Ni-3%Fe钨合金从准静态10-4 s-1直到动态7 000 s-1范围内的塑性流动行为,发现在高应变率加载下几乎没有应变硬化效应。Lee等[6-7]较为全面地测试了几种钨合金在不同温度和应变率下的力学性能,发现钨合金颗粒变形与应变率水平有关。目前针对钨合金材料的研究主要集中在其绝热剪切性能上,而对这种典型弹体材料在高温高应变率耦合作用下的力学性能和本构模型研究还不够全面,不能满足弹丸优化设计、侵彻过程数值模拟的需要。所以,全面研究钨合金材料在不同状态(包括极端温度环境和应变率条件)下的力学行为,给出其相适应的本构模型,对于在国防军事领域上更好地应用钨合金材料具有实际指导意义。
针对本研究的钨合金材料,首先进行不同温度和应变率下的压缩实验,系统地研究钨合金的力学特性和变形机制,实验的温度变化范围为296~1 273 K,应变率范围为0.000 5~6 000 s-1;然后,针对实验得到的应力-应变关系,分别采用两种唯象本构模型(JC本构模型[8]和KHL本构模型[9-10])和一种基于物理概念的本构模型(PB本构模型[11-12])进行拟合,并将拟合结果与实验结果进行对比;最后通过误差分析,评价这3种本构模型的预测效果。
2. 实验过程与结果
2.1 材料与实验
实验材料为93%W-4.9%Ni-2.1%Fe钨合金。在准静态压缩实验中,试样的尺寸为∅5 mm×8 mm;在动态压缩实验中,试样的尺寸为∅4 mm×4 mm。为了减小端面摩擦的影响,每次实验前试样两端均涂上润滑剂。
准静态压缩实验采用CSS4410型电子万能材料试验机及自带的加热炉,开展两个应变率0.000 5和0.05 s-1下的加载。动态压缩实验(800 s-1和6 000 s-1)采用分离式Hopkinson压杆,高温实验装置为自行设计的加热炉以及温度反馈控制仪。实验温度为296、673、1 073和1 273 K。当温度高于573 K时,动态压缩实验必须采用同步组装装置,即在实验加热过程中,入射杆和透射杆置于加热炉外面,以避免加热试样过程中入射杆和透射杆杆端温升过高,对加载波的传播造成干扰。该同步组装装置的结构及工作原理参见文献[13-14]。另外,由于动态压缩实验应变率很高,使试样在变形过程中与入射杆和透射杆接触的时间很短,以应变率800 s-1、最大应变0.6计,整个变形过程不到1 ms,在此时间内,试样上的热量损失可不予考虑。
2.2 实验结果
图 1所示为实验获得的钨合金在不同应变率、不同温度下的流动应力-塑性应变曲线。从图 1中可以看出,在不同的条件下,流动应力差别很大。在相同应变条件下,流动应力随着温度升高而急剧减小,随着应变率增加而增大,表现出明显的温度软化效应和应变率硬化效应,这也是体心立方(bcc)材料典型的性能特征。
从图 1中还可以看出,钨合金材料在低应变率加载下,温度越低应变硬化越明显:温度为296 K时,钨合金的应变硬化特别明显,流动应力从屈服时的600 MPa增加到接近1 600 MPa(应变达到0.6),增幅达到160%左右;随着温度的升高,应变硬化程度逐渐减弱,温度为1 273 K时,钨合金应变硬化几乎消失,如图 1(a)所示,这是由于温度软化对材料造成的影响,即温度越高,软化效果越明显。在高应变率加载下,钨合金材料却表现出不一样的行为,如图 1(d)所示,当温度为296 K、塑性应变超过0.2之后出现应变软化现象,在高温下却始终没有出现应变软化阶段。这是因为在同一高应变率下,温度较低时材料的流动应力较大,因此同一应变下其单位体积的塑性功较大,且在高应变率加载下,塑性功转化产生的热量来不及扩散,瞬间使材料产生绝热温升,从而导致材料的热软化效应更加明显,宏观上表现为应变软化效应更明显。从图 1(d)可以看出,常温下钨合金的塑性流动应力将近为高温下的两倍,其从塑性功转化的热量导致的温升要高得多,这可能就是导致常温下出现应变软化现象的原因。
3. 本构模型
金属材料的塑性流动应力强烈依赖于材料所经受的应变率(˙ε)、温度(T)和应变。因此,要想准确地描述材料的塑性流动应力,本构模型中必须考虑到应变率、温度和应变这3个最主要的影响因素。在描述金属材料的塑性流动应力的本构模型中,JC本构模型[8]、KHL本构模型[9-10]和PB本构模型[11-12]应用较广泛, 原因有二:一是这3个模型均考虑了应变率、温度和应变对流动应力的影响;二是相比于其他模型,这3个模型形式较简单,参数较少,且容易通过实验数据获得。本研究分别采用这3个模型来描述钨合金的塑性流动应力。在动态实验中,塑性变形产生的热无法很快扩散,变形的绝热过程引起温升,温升ΔT可以根据塑性功转变式得出[15]
ΔT=ηρcV∫σdε (1) 式中:ρ是材料密度(17.04 g/cm3); cV是材料的定容比热容(134 J/(kg·K)); η是塑性功转化为热的比例系数,本研究假设所有的塑性功全部转换成热,即η取1;ε为等效塑性应变;σ为塑性流动应力。
3.1 JC本构模型
在JC本构模型[8]中,塑性流动应力σ的表达式如下
σ=(A+Bεn)(1+Cln˙ε∗)(1−T∗m) (2) 式中:˙ε∗=˙ε/˙ε0,˙ε0=0.000 5 s-1是参考应变率; T*=(T-Tr)/(Tm-Tr),Tm是材料的熔点(钨合金材料的熔点为1 730 K),Tr=296 K是参考温度; A、B、n、C、m是5个待定参数,取值列于表 1。图 2为实验数据与JC模型拟合曲线的对比。
表 1 JC模型参数Table 1. Determined values of JC model parametersA/(MPa) B/(MPa) n C m 600.8 1 200 0.494 4 0.059 0 0.820 3 3.2 KHL本构模型
σ=[A+B(1−ln˙εlnD0)n1εn0](˙ε˙ε0)C(Tm−TTm−Tr)m (3) 式中:D0=106 s-1, 是一个常量[9-10]; A、B、n1、n0、C、m是6个模型参数,取值列于表 2; 其他参数含义与JC模型一致。图 3为实验数据与KHL模型拟合曲线的对比。
表 2 KHL模型参数Table 2. Determined values of KHL model parametersA/(MPa) B/(MPa) n1 n0 C m 264.5 523.4 2.296 0.299 6 0.113 4 1.063 3.3 PB本构模型
σ=σa+σ∗ (4) 式中:σa代表塑性流动应力的非热部分,σ*代表塑性流动应力的热激活部分。σ*的表达式如下
σ∗={σ0[1−(−kTG0ln˙ε˙ε0)1/q]1/pT≤Tc0T>Tc (5) Tc=−G0kln(˙ε/˙ε0) (6) 式中:Tc为临界温度;σ0代表 0 K时的临界应力;˙ε0可视为参考应变率;参数k/G0表征材料的温度敏感性;0<p≤1和1≤q≤2定义了短程障碍对位错的阻力构型,对于绝大部分bcc材料,p=0.5和q=1.5是适用的[11-12],本研究也选择这两个值。塑性流动应力的非热部分σa通常以指数形式描述
σa=aεn (7) 式中:a和n是材料参数,需要通过实验数据拟合得到。因此,PB模型的塑性流动应力表达式整理如下
σ=aεn+σ0[1−(−kTG0ln˙ε˙ε0)1/q]1/p (8) 钨合金材料的PB本构模型参数值列于表 3。图 4为实验数据与PB模型拟合曲线的对比。
表 3 PB模型参数值Table 3. Determined values of PB model parametersa/(MPa) n σ0/(MPa) (k/G0)/(K-1) ˙ε0/(s-1) p q 1 054.4 0.263 5 3 061.4 3.89×10-5 2×1010 0.5 1.5 4. 本构模型的比较与验证
4.1 本构模型的比较
为了比较上述3个本构模型的拟合效果,采用误差计算公式进行评价,其形式如下
er=1NN∑i|σiexp−σimodelσiexp|×100% (9) 式中:σexp和σmodel分别是实验测得的流动应力和本构模型预测的流动应力,er为预测值与实验值之间的误差,N为实验数据点总数。根据(9)式,3个本构模型的预测效果如图 5所示。可见,3个模型的预测效果都比较接近,预测误差值均在10%左右;但相对来说,PB模型的预测效果最好,误差值最低。
图 5还列出了Xu等[16]的研究结果。可见,PB模型的预测效果最好,误差值最低,与本研究结论一致,证实了本研究结论的普适性。而Xu等[16]的每个本构模型的预测误差值均比本研究小,是因为他们研究的实验应变率范围(0.001~3 000 s-1)和温度范围(173~873 K)均更窄。
4.2 本构模型的验证
对于得到的本构模型,在其适用范围内,期望它能描述任意工况下的材料塑性流变行为。为此,采取应变率跳跃试验验证3个本构模型的预测能力。
首先,在实验温度873 K和应变率1 000 s-1条件下,使用分离式Hopkinson压杆进行动态加载,得到试样的应力-应变曲线;然后,在实验温度296 K和应变率4 000 s-1条件下,采用同样方法得到试样的应力-应变曲线。前者取塑性应变小于0.15的数据点,后者取塑性应变大于0.15的数据点。最终,将实验得到的数据与本构模型预测的结果进行比较。
图 6显示了实验结果与本构模型预测结果的对比。在第1个工况下,PB模型预测值能非常准确地贴合实验测量值,而JC模型和KHL模型的预测值均高于实验测量值;在第2个工况下,KHL模型预测值略低于实验测量值,PB模型预测值刚开始偏高,随后慢慢贴合实验测量值,而JC模型预测值一直明显高于实验测量值。为了定量地评价3个模型预测效果的好坏,通过(9)式计算得到, JC模型、KHL模型和PB模型的预测误差值依次是18.59%、16.20%和2.29%,如图 7所示。可以看出,PB模型的预测效果明显好于JC模型和KHL模型的预测效果。
5. 结论
系统全面地研究了钨合金在不同应变率(0.000 5~6 000 s-1)和不同温度(296~1 273 K)下的塑性流动行为和本构关系,发现钨合金材料表现出明显的温度软化效应和应变率硬化效应,其应变硬化率随着应变率的增加或温度的升高而急剧变小。
基于实验数据,给出了两种唯象本构模型(JC本构模型和KHL本构模型)和一种基于物理概念模型(PB本构模型)的拟合结果。可以得出,3个模型对实验数据的预测效果都比较接近,平均误差值均在10%左右;但相对来说,PB模型的预测效果最好,平均误差值最低。该结论与Xu等[16]的结论一致。在对于应变率跳跃试验的预测中,JC模型和KHL模型的预测误差均在15%以上,而PB模型的预测误差仅为2.29%。因此,可以得出:相比于JC模型和KHL模型,PB模型能更准确地预测钨合金的塑性流动行为。
造成3种本构模型预测误差值不同的原因在于:PB模型考虑了材料变形的微观机制,能更贴近材料的真实变形行为;而JC模型和KHL模型是基于对材料塑性变形规律的归纳总结而得到的,未考虑材料变形的微观机制,因此预测效果不如PB模型好。但对于工程应用来说,JC模型和KHL模型形式简单,参数获取较PB模型容易,使用更广泛。
-
表 1 JC模型参数
Table 1. Determined values of JC model parameters
A/(MPa) B/(MPa) n C m 600.8 1 200 0.494 4 0.059 0 0.820 3 表 2 KHL模型参数
Table 2. Determined values of KHL model parameters
A/(MPa) B/(MPa) n1 n0 C m 264.5 523.4 2.296 0.299 6 0.113 4 1.063 表 3 PB模型参数值
Table 3. Determined values of PB model parameters
a/(MPa) n σ0/(MPa) (k/G0)/(K-1) ˙ε0/(s-1) p q 1 054.4 0.263 5 3 061.4 3.89×10-5 2×1010 0.5 1.5 -
[1] 范景莲, 刘涛, 成会朝, 等.钨合金在高速加载条件下的动态力学性能和失效机理[J].稀有金属材料与工程, 2006, 35(6):841-844. doi: 10.3321/j.issn:1002-185X.2006.06.001FAN J L, LIU T, CHENG H C, et al.Dynamic properties and failure mechanism of tungsten heavy alloy under high-speed loading[J].Rare Metal Materials and Engineering, 2006, 35(6):841-844. doi: 10.3321/j.issn:1002-185X.2006.06.001 [2] 魏志刚, 胡时胜, 李永池, 等.微细观结构对预扭转钨合金材料绝热剪切的影响[J].弹道学报, 1999, 11(1):15-19. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199900042055WEI Z G, HU S S, LI Y C, et al.Structural influence on the adiabatic shearing failure of pre-torqued tungsten heavy alloy[J].Journal of Ballistics, 1999, 11(1):15-19. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199900042055 [3] KIM D S, NEMAT-NASSER S, ISAACS J B, et al.Adiabatic shear band in WHA in high-strain-rate compression[J].Mech Mater, 1998, 28(1):227-236. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=78c0a5aa25b27f88859bed95e1323e56 [4] YADAV S, RAMESH K T.The mechanical properties of tungsten-based composites at very high strain rates[J].Mater Sci Eng A, 1995, 203(1/2):140-153. http://www.sciencedirect.com/science/article/pii/0921509395098658 [5] COATES R S, RAMESH K T.The rate-dependent deformation of a tungsten heavy alloy[J].Mater Sci Eng A, 1991, 145(2):159-166. doi: 10.1016/0921-5093(91)90245-I [6] LEE W S, CHIOU S T.The influence of loading rate on shear deformation behaviour of tungsten composite[J].Compos Part B Eng, 1996, 27(2):193-200. doi: 10.1016/1359-8368(95)00051-8 [7] LEE W S, XIEA G L, LIN C F.The strain rate and temperature dependence of the dynamic impact response of tungsten composite[J].Mater Sci Eng A, 1998, 257(2):256-267. doi: 10.1016/S0921-5093(98)00852-1 [8] JOHNSON G R, COOK W H.A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures[C]//Proceedings of the 7th International Symposium on Ballistics.AIP, 1983. [9] KHAN A S, LIANG R.Behaviors of three BCC metal over a wide range of strain rates and temperatures:experiments and modeling[J].Int J Plasticity, 1999, 15(10):1089-1109. doi: 10.1016/S0749-6419(99)00030-3 [10] KHAN A S, SUH Y S, KAZMI R.Quasi-static and dynamic loading responses and constitutive modeling of titanium alloys[J].Int J Plasticity, 2004, 20(12):2233-2248. doi: 10.1016/j.ijplas.2003.06.005 [11] NEMAT-NASSER S, GUO W.Flow stress of commercially pure niobium over a broad range of temperatures and strain rates[J].Mater Sci Eng A, 2000, 284(1):202-210. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1f51f05b708b795911767017d4b0e6e9 [12] NEMAT-NASSER S, GUO W G.Thermomechanical response of DH-36 structural steel over a wide range of strain rates and temperatures[J].Mech Mater, 2003, 35(11):1023-1047. doi: 10.1016/S0167-6636(02)00323-X [13] 李玉龙, 索涛, 郭伟国, 等.确定材料在高温高应变率下动态性能的Hopkinson杆系统[J].爆炸与冲击, 2005, 25(6):487-492. doi: 10.3321/j.issn:1001-1455.2005.06.002LI Y L, SUO T, GUO W G, et al.Determination of dynamic behavior of materials at elevated temperatures and high strain rates using Hopkinson bar[J].Explosion and Shock Waves, 2005, 25(6):487-492. doi: 10.3321/j.issn:1001-1455.2005.06.002 [14] LI Y, GUO Y, HU H, et al.A critical assessment of high-temperature dynamic mechanical testing of metals[J].Int J Impact Eng, 2009, 36(2):177-184. doi: 10.1016/j.ijimpeng.2008.05.004 [15] KAPOOR R, NEMAT-NASSER S.Determination of temperature rise during high strain rate deformation[J].Mech Mater, 1998, 27(1):1-12. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=7c17e09746f2eec0386fd35b3d3e02c8 [16] XU Z, HUANG F.Thermomechanical behavior and constitutive modeling of tungsten-based composite over wide temperature and strain rate ranges[J].Int J Plasticity, 2013, 40(1):163-184. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=4a31669ec7b7a87b9c1b72c76c7f39a2 -