甘春雷，郑开宏，王海艳，戚文军，周 楠，蔡 畅

广东省工业技术研究院(广州有色金属研究院)，广东 广州 510650

随着人们环保意识的增强以及对自身健康的关注，无铅易切削黄铜合金的开发和应用研究逐渐得到重视.相关研究者对无铅易切削黄铜合金开展了广泛的研究[1-4]，并取得了一些研究成果.然而，由于无铅易切削黄铜合金在热变形后易出现组织和性能缺陷，极大地限制了其实用性，因此，有必要对无铅易切削黄铜合金的热变形行为进行深入的研究.迄今为止，对无铅黄铜合金的研究主要集中在切削性能、拉伸性能和腐蚀性能等方面[5-7]，而对其热变形行为的研究报道较少.

材料在热变形过程中的流变应力通常采用本构方程来进行表征，通过本构方程将流变应力、变形温度和应变速率联系起来，分析材料在热变形过程中的本构特征.目前已有多种流变应力本构模型被用来模拟材料在热变形过程中的流变应力变化情况[8-12].在这些模型中，Sellars C M等人[13]提出的现象本构模型被广泛用来预测材料的热变形行为.近年来，该模型被进一步修正，发展出了考虑应变补缩的本构模型来模拟不同合金的热变形行为[14-17]，对热变形行为的预测更加精确.

本研究以无铅易切削黄铜合金为对象，通过在Gleeble-1500热模拟机上进行圆柱体压缩试验，研究在热压缩变形时，变形温度和应变速率对合金的流变应力的影响，获得无铅易切削黄铜合金在不同变形工艺参数条件下的应力值，建立该合金在高温下变形时的流变应力预测本构模型(考虑应变补缩)，为开展无铅易切削黄铜合金加工过程的数值模拟、加工工艺的制订以及成形过程组织性能精确控制提供实验数据和理论基础.

1 实验材料及方法

试验材料为采用普通浇铸方法制备的无铅易切削黄铜合金，其成分为(质量分数，%)：Cu-36.3Zn-1.76Mg-0.5Bi-0.5Sb，铸态微观组织如图1所示.采用线切割机，将样品制成直径10 mm，长度15 mm的圆柱体.压缩试验在Gleeble-1500热模拟机上进行.进行压缩试验时，在压头与试样两端接触处夹一层1 mm厚的石墨薄片进行润滑，以减少摩擦对应力及变形状态的影响.设定压缩试验变形温度范围为823～973 K，应变速率范围为0.01～1 s-1，总压缩变形量为60%.压缩试验完成后，立即水淬.采用LEICA-DMI3000M型金相显微镜对试样中部进行金相组织观察，所用浸蚀剂配比为：40 mL硝酸+40 mL冰醋酸+20 mL水.

图1 无铅易切削黄铜合金铸态微观组织

Fig.1Microstructure of the as-cast lead-free machinable brass

2 实验结果与讨论

2.1 流变应力变化

通过计算机采集数据，获得无铅易切削黄铜合金压缩变形时的变形温度、应变速率、真应变和流变应力等数据，绘制出在不同条件下，无铅易切削黄铜合金的压缩变形真应力-真应变曲线(图2).由图2可见：变形温度和应变速率对合金的流变应力具有很大的影响，随着应变速率的增加和变形温度的降低，流变应力增加.进一步观察可以发现，随着应变的增加，流变应力快速增加到峰值，然后缓慢减少至稳态值，在某些条件下随着应变的增加，流变应力会减少至更低值.这些特征表明，在本文的试验条件下，无铅易切削黄铜合金在热压缩变形过程中出现了明显的再结晶现象.

图2 无铅易切削黄铜合金在温度823～973 K和不同应变速率下的真应力-真应变曲线

Fig.2True stress-true strain curves of the lead-free machinable brass alloy in a temperature range of 823-973 K with strain rate of (a) 0.01 s-1, (b) 0.1 s-1, (c) 1 s-1

2.2 微观组织观察

无铅易切削黄铜合金在热变形条件下，其流变应力出现峰值，这与合金中发生了动态再结晶密切相关.在变形初始阶段，由于加工硬化占主导地位，合金中的硬化作用大大超过软化作用；随变形程度的继续增加，位错密度不断增高，加快了动态回复和动态再结晶，软化作用逐渐增强，使得流变应力达到峰值后开始下降.这种动态再结晶现象可以通过试样热变形后的微观组织得到证明(图3).由图3可见，在变形温度为823 K，应变速率为1 s-1的条件下，在箭头所指区域内出现了新的尺寸较小的动态再结晶晶粒.

图3823 K，1 s-1条件下压缩试样组织

Fig.3Microstructure of the compression specimen: 823 K，1 s-1

另外，无铅易切削黄铜合金在热变形条件下出现动态软化与发生了β→α相变也有关系.图4是在应变速率0.01～1 s-1，温度823～973 K的条件下，合金的典型金相组织照片.由图4可见，在温度低于923 K时，由于试样未完全发生动态再结晶，因而存在再结晶组织和变形带组织共存的现象(见图4(e)).当温度超过923 K时，由于发生了相变，合金的微观组织出现了显著的变化.例如，在变形温度为973K时，出现了细小的针状α相(见图4(d)和(f)).由于α相(面心立方结构)比β相(体心立方结构)具有更好的塑性变形能力，因此，β相增多有利于提高无铅易切削黄铜合金在热变形过程中的变形协调能力，从而降低变形阻力[18].

图4 典型的压缩试样微观组织

2.3 流变应力本构模型的建立

从无铅易切削黄铜合金的热压缩流变曲线可知，流变应力的大小取决于变形温度和应变速率，这种力学变形行为可采用Arrhenius方程进行描述：

(1)

其中：

(2)

将式(2)中的F(σ)代入式(1)并求对数得到：

(3)

(4)

(5)

热变形条件下通常采用温度补偿的变形速率因子Zener-Hollomon(Z)参数，有

(6)

得到适用于工艺分析的流变应力模型：

(7)

变形程度对变形激活能和材料常数具有重要的影响，为了提高流变应力模型的预测精度，将材料常数表示成真应变ε的多项式函数.在真应变0.05～0.85范围内进行拟合，从1至9次进行多项式拟合，发现6次多项式具有较好的拟合结果(图5).获得了如式(8)所示的多项式方程组，式(8)中相应的系数列于表1.

图5 无铅易切削黄铜合金6次多项式拟合获得的(a) α′, (b) n, (c) Q, (d) lnA与真应变的关系

Fig.5Relationships between (a)α′, (b)n, (c)Q, (d) lnAand true strain by sixth polynomial fit of the lead-free machinable brass alloy

(8)

表1 α′，n，Q和lnA的多项式拟合系数Table 1 Coefficients of the polynomial for α′,n, Q and lnA

为了验证所建立的流变应力模型的可靠性，采用相关系数R和平均绝对误差AARE进行表征：

(9)

(10)

图6 流变应力预测值与实验值比较

Fig.6Correlation between predicted and experimental stresses

3 结 论

(1)无铅易切削黄铜合金热压缩流变应力随应变速率的减小和变形温度的增加而减小，其相互关系可用Z参数进行表征.

(2)由金相组织观察可知，在应变速率为0.01～1 s-1的条件下，当变形温度低于923 K时发生动态再结晶，在变形温度超过923 K时发生β→α相变，这两种作用使流变应力出现动态软化.

(3)所建立的无铅易切削黄铜合金流变应力模型的相关系数为0.99815，平均绝对误差为6.105%，表明该模型具有较好的预测精度.

参考文献：

[1] MOHAMED A T, NAHED A E, RAWIA M H, et al. Machinability characteristics of lead free-silicon brass alloys as correlated with microstructure and mechanical properties[J]. Ain Shams Engineering Journal,2012,3(4):383-392.

[2] LI S F, KONDOH K, IMAI H, et al. Fabrication and properties of lead-free machinable brass with Ti additive by powder metallurgy[J]. Powder Technology, 2011,205(1-3):242-249.

[3] TOSHIKAZU M, TAKAYUKI O. Cutting of lead-free copper alloy “Eco Brass”[J]. Journal of Japanese Research Institute Advanced Copper-Base Material Technology,2006, 45:250-255.

[4] EMELINA N B, ALABIN A N, BELOV N A. Influence of bismuth and lead on the formation of the structure of experimental alloys with a Cu-30%Zn based composition during crystallization, deformation, and thermal treatment[J]. Russian Journal of Non-ferrous Metals, 2010,51(6):476-482.

[5] XIAO L R, SHU X P, YI D Q, et al. Microstructure and properties of unleaded free-cutting brass containing stibium[J]. Trans Nonferrous Met Soc China, 2007, 17: 1055-1059.

[6] ATSUMI H, IMAI H, LI S F, et al. High-strength, lead-free machinable-duplex phase brass Cu-40Zn-Cr-Fe-Sn-Bi alloys[J]. Materials Science and Engineering A, 2011, 529:275-281.

[7] YOU S J, CHOI Y S, KIM J G, et al. Stress corrosion cracking properties of environmentally friendly unleaded brasses containing bismuth in Mattsson’s solution[J]. Materials Science and Engineering A, 2003, 345:207-214.

[8] MOMENI A, DEHGHANI K. Characterization of hot deformation behavior of 410 martensitic stainless steel using constitutive equations and processing maps[J]. Materials Science and Engineering A, 2010, 527(21-22):5467-5473.

[9] XIAO M L, LI F G, ZHAO W, et al. Constitutive equation for elevated temperature flow behavior of TiNiNb alloy based on orthogonal analysis[J]. Materials and Design, 2012, 35:184-193.

[10] SHI L, YANG H, GUO L G, et al. Constitutive modeling of deformation in high temperature of a forging 6005A aluminum alloy[J]. Materials and Design, 2014, 54: 576-581.

[11] GAO L, LUO A A. Hot deformation behavior of as-cast Mg-Zn-Mn-Ce alloy in compression[J]. Materials Science and Engineering A, 2013, 560: 492-499.

[12] LIN Y C, DING Y, CHEN M S, et al. A new phenomenological constitutive model for hot tensile deformation behaviors of a typical Al-Cu-Mg alloy[J]. Materials and Design, 2013, 52:118-127.

[13] SELLARS C M, MCTEGART W J. On the mechanism of hot deformation[J]. Acta Metallurgical, 1996, 14(9):1136-1138.

[14] MIRZAEI A, ZAREI-HANZAKI A, HAGHDADI N, et al. Constitutive description of high temperature flow behavior of Sanicro-28 super-austenitic stainless steel[J]. Materials Science and Engineering A, 2014, 589: 76-82.

[15] ZHONG J, XIAO Y H, GUO C. Constitutive modelling for high temperature behaviour of 12CrNiMoWV martensitic stainless steel[J]. Materials Science Technology, 2012, 28(6):719-726.

[16] YIN F, HUA L, MAO H J, et al. Constitutive modeling for flow behavior of GCr15 steel under hot compression experiments[J]. Materials and Design, 2013, 43:393-401.

[17] LUAN J, SUN C, LI X, et al. Constitutive model for AZ31 magnesium alloy based on isothermal compression test[J]. Materials Science Technology, 2014, 30(2):211-219.

[18] XIAO Y H, GUO C, GUO X Y. Constitutive modeling of hot deformation behavior of H62 brass[J]. Materials Science and Engineering A, 2011, 528(21): 6150-6158.