純銅擠扭工藝數(shù)值模擬與實驗研究

王 成，李 萍，薛克敏，李 曉，章 凱

(合肥工業(yè)大學(xué) 材料科學(xué)與工程學(xué)院，合肥 230009)

純銅擠扭工藝數(shù)值模擬與實驗研究

王 成，李 萍，薛克敏，李 曉，章 凱

(合肥工業(yè)大學(xué) 材料科學(xué)與工程學(xué)院，合肥 230009)

運用有限元模擬軟件 Msc.Marc對純銅擠扭工藝變形特征進行研究，分析擠扭變形過程中應(yīng)變分布以及應(yīng)變分量的變化趨勢。采用自行設(shè)計的擠扭模具，在室溫下進行了純銅的一道次擠扭試驗。結(jié)果表明：變形可分為變形開始、完全充滿、逐步擠出3個階段；應(yīng)變呈中心低、邊緣高的分布；變形過程中正應(yīng)變分量小，切應(yīng)變分量較大，變形以剪切變形為主。擠扭過程中存在滑移和孿生兩種變形方式；擠扭后試樣的硬度得到顯著提高，在橫截面呈現(xiàn)邊緣高、中心低的分布趨勢。

純銅；擠扭變形；數(shù)值模擬；顯微組織；顯微硬度

目前，大塑性變形法(SPD)已經(jīng)成為制備超細(xì)晶材料的重要方法[1?4]。擠扭(Twist extrusion, TE)法是一種新型的大塑性變形工藝(見圖1)，該工藝將非圓形截面試樣經(jīng)螺旋通道擠出。由于一道次變形前后試樣橫截面的形狀和尺寸均不發(fā)生改變，可以重復(fù)變形以累積大的塑性變形來達(dá)到細(xì)化晶粒和提高材料性能的目的[5]。

近年來，國外學(xué)者對擠扭工藝開展了相關(guān)理論和實驗研究，但對擠扭的變形方式和細(xì)化機理研究還不夠深入[6?11]，國內(nèi)關(guān)于擠扭工藝的研究還少見報道[12?13]。本文作者運用有限元分析和實驗相結(jié)合的方法對純銅的擠扭一道次變形過程進行仿真，分析了擠扭變形過程以及變形時應(yīng)變分量變化情況；通過純銅一道次擠扭實驗進一步研究擠扭變形后材料組織及性能變化。

圖1 擠扭變形原理Fig.1 Schematic diagram of twist extrusion: (a) Sample; (b)Twist channel

1 有限元模型

采用MSC.Marc有限元軟件進行模擬仿真。沖頭與凹模采用剛性面簡化處理，螺旋通道與直通道連接處采用圓弧過渡，螺旋通道扭轉(zhuǎn)角度為 90°，具體尺寸如圖2所示。沖頭的壓下速度為1 mm/s, 模擬過程中不考慮溫度的影響。試樣與模具的接觸均采用庫倫摩擦模型，摩擦因數(shù)取0.2。

圖2 有限元模型Fig.2 FEM model (mm)

試樣為橫截面尺寸10 mm×10 mm、長度40 mm的純銅長方體坯料，材料選用如下Johnson-Cook本構(gòu)模型：

采用六面體實體單元對坯料進行劃分，獲得單元數(shù)為4 000，節(jié)點數(shù)4 961。

表1 銅的Johnson-Cook模型參數(shù)值Table 1 Parameters of pure copper to Johnson-Cook mold

2 模擬結(jié)果與分析

2.1 變形過程分析

擠扭變形過程決定了試樣的應(yīng)變分布，對晶粒細(xì)化和性能改善有重要影響。

圖3所示為一道次擠扭變形過程示意圖，可分為3個階段：第一階段為變形開始，此時端部被擠入螺旋通道，開始發(fā)生剪切變形，擠壓力迅速增加(最大為45 kN)，直至試樣完全充滿螺旋通道；第二階段為試樣完全充滿螺旋通道階段后，端部被擠出螺旋通道進入方形通道，擠壓力進一步增加(約為80 kN)；第三階段為試樣開始被完全擠出，擠壓力開始逐漸降低。擠壓力—行程曲線如圖4所示。

圖3 擠扭變形過程Fig.3 Procedure steps of TE: (a) Initial stage; (b) Steady stage;(c) Final stage

圖4 擠壓力—行程曲線Fig.4 Load—stroke curve during TE process

2.2 應(yīng)變分布

圖5所示為擠扭變形后試樣縱截面的應(yīng)變分布圖。由圖5可知，經(jīng)一道次擠扭后試樣整體發(fā)生較大變形。通過縱截面應(yīng)變分布圖分析可知：試樣邊緣處變形量較中心處大，中心塑性應(yīng)變?yōu)?0.16，邊緣為0.64。試樣兩端存在難變形區(qū)域(圓圈所示)，與沖頭接觸的左側(cè)一端由于受到摩擦作用，該區(qū)金屬難以流動，而另一端為自由面，僅邊緣受到模具的阻力，近乎剛性平移，變形較小。

圖5 縱截面應(yīng)變分布圖Fig.5 Effective strain distribution of longitudinal section

2.3 應(yīng)變分量

在MSC.Marc中可以提取各個應(yīng)變分量的變化，包含與法向和切向位移分量相關(guān)的comp11，comp22，comp33，comp12，comp23，comp31。前3個分別為垂直于XY，YZ，ZX面的正應(yīng)變分量，后3個為平行于XY，YZ，ZX面的切應(yīng)變分量。

按圖6所示選取節(jié)點(選取原則為避開端部難變形區(qū))，正應(yīng)變分量和切應(yīng)變分量的變化趨勢如圖所示。由圖6可知，在整個變形過程中3個正應(yīng)變分量值均較小(最大絕對值均在0.02左右)，而切應(yīng)變值相對較大(最大絕對值分別為0.15、0.30、0.35)，這表明擠扭過程以剪切變形為主。而Comp23、31值也說明擠扭過程中剪切變形不僅存在于橫截面上，垂直橫截面方向上也存在。這與 BEYGELZIMER等[14]指出擠扭變形具有兩個剪切面的觀點是相符合的。

切應(yīng)變分量的變化與擠扭變形過程有關(guān)，變形開始時選取點在直通道內(nèi)做剛性平移，各切應(yīng)變分量幾乎沒有變化，當(dāng)進入螺旋通道后各應(yīng)變分量開始發(fā)生變化。直至試樣完全被擠出螺旋通道后，各應(yīng)變分量又趨于不變。

圖6 節(jié)點應(yīng)變分量的變化Fig.6 Change of strain component of node: (a) Normal strain;(b) Shear strain

3 實驗結(jié)果與分析

3.1 實驗過程

為了驗證模擬分析的可靠性，以及進一步研究擠扭變形的特征及組織演化的規(guī)律，試驗條件盡量與模擬相一致。在室溫下進行純銅(Cu含量大于99.89%，質(zhì)量分?jǐn)?shù))的一道次擠扭試驗，對擠扭后試樣進行顯微組織觀察并對試樣橫截面進行顯微硬度測量。

圖7所示為擠扭變形前后試樣的宏觀照片。變形前試樣為10 mm×10 mm×40 mm的長方體坯料，經(jīng)擠扭變形后試樣產(chǎn)生了輕度的扭曲，主要是由于擠扭變形過程中試樣發(fā)生了扭轉(zhuǎn)變形，模具擠出端的矩形通道的長度不足(僅為30 mm)，導(dǎo)致試樣不能完被全矯正到初始形狀，這也從另一個角度反映了擠扭變形的劇烈程度。

圖7 擠扭變形前后試樣的宏觀照片F(xiàn)ig.7 Photos of samples before TE and after one pass

3.2 顯微組織

將試樣沿剖切面切開，取橫截面進行磨制、拋光、腐蝕觀察顯微組織，圖8(a)所示為擠扭變形前試樣組織，基體組織呈近似等軸晶狀，晶粒排列均勻規(guī)則。圖8(b)所示為擠壓后試樣中心的顯微組織，較初始態(tài)，晶粒并未發(fā)生明顯的扭轉(zhuǎn)、破碎變形。根據(jù)模擬結(jié)果，中心的應(yīng)變小，變形程度有限，試樣中心區(qū)域晶粒形狀和尺寸變化較小。銅屬于FCC晶格金屬，晶粒變形以滑移為主，但從圖8(b)可以觀察到一定數(shù)量的孿晶組織(圈中所示)，這是由于純銅屬于層錯能較低(僅70 mJ/m2)的金屬，層錯能越低, 位錯越容易分解為層錯，束集越困難，越難以產(chǎn)生交滑移，越容易發(fā)生孿生變形。因此，晶粒細(xì)化是孿生與位錯滑移共同協(xié)調(diào)的變形機制。這與文獻(xiàn)[15]的結(jié)果相似。圖8(c)所示為變形后試樣的邊緣組織，由于該區(qū)域變形量較大，晶粒沿扭轉(zhuǎn)方向被嚴(yán)重拉長，等軸晶完全消失，呈纖維狀分布的織構(gòu)組織，由于晶粒完全變形成纖維狀，因此，觀察不到孿晶。

圖8 擠扭變形前后顯微組織圖Fig.8 Microstructures of samples before and after TE:(a) Initial; (b) Center; (c) Edge

3.3 顯微硬度

圖9所示為擠扭后試樣橫截面硬度實測值的三維分布圖。純銅初始硬度為62.7。由圖9可以看出，擠扭變形后中心處硬度值最小(約為87)，邊緣硬度值大，在4個棱角處達(dá)到最大(約為130)。硬度的增加主要與變形量有關(guān)，變形導(dǎo)致晶粒細(xì)化和加工硬化，由于擠扭過程中晶粒通過滑移和孿生的變形方式，晶粒細(xì)化導(dǎo)致硬度提高。并且隨著變形的加劇，其內(nèi)部產(chǎn)生較多位錯，形成了纏結(jié)，并構(gòu)成胞狀亞結(jié)構(gòu)，位錯之間相互干擾, 使滑移困難，變形抗力增加導(dǎo)致金屬硬化。這也解釋了剪切應(yīng)變大的區(qū)域硬度值也較高的原因。

從試樣橫截面組織形貌以及硬度分布，結(jié)合模擬結(jié)果分析可知，應(yīng)變量大的區(qū)域組織變形越劇烈，硬度值也得到顯著提高，應(yīng)變量小的區(qū)域則變化不明顯。模擬分析很好地解釋了實驗結(jié)果。

圖9 擠扭后試樣橫截面硬度實測值的三維分布圖Fig.9 Measured microhardness profile of cross-section of samples after one pass TE

4 結(jié)論

1) 擠扭變形過程分為3個階段，即變形開始、完全充滿、逐步擠出階段；擠扭變形后，試樣等效塑性應(yīng)變呈中心低(0.16)、邊緣高(0.64)的分布；擠扭變形過程中正應(yīng)變分量較小，而切應(yīng)變分量相對較大。

2) 擠扭后，試樣邊緣晶粒變形程度大于中心處的，邊緣處硬度(130)高于中心處硬度(87)；擠扭變形過程中，存在滑移和孿生兩種變形方式，純銅晶粒變形為兩種變形方式的共同協(xié)調(diào)機制。

REFERENCES

[1] SEGAL V M. Severe plastic deformation: simple shear versus pure shear[J]. Materials Science and Engineering A, 2002, 338:331?344.

[2] 李 萍, 黃科帥, 薛克敏, 周明智, 韓國民. 純鋁粉末多孔燒結(jié)材料等通道轉(zhuǎn)角擠壓[J]. 中國有色金屬學(xué)報, 2009, 19(5):881?886.LI Ping, HUANG Ke-shuai, XUE Ke-min, ZHOU Ming-zhi,HAN Guo-min. Equal channel angular pressing of pure Al powder sintered material[J]. The Chinese Journal of Nonferrous Metals, 2009, 19(5): 881?886.

[3] VALIEV R Z, ISLAMGALIEV R K, ALEXANDROV I V. Bulk nanostructured materials from severe plastic deformation [J].Progress in Materials Science, 2000, 45: 103?189.

[4] 周明智, 薛克敏, 李 萍. 粉末多孔材料等徑角擠壓過程熱力耦合有限元數(shù)值分析[J]. 中國有色金屬學(xué)報, 2006, 16(9):1510?1516.ZHOU Ming-zhi, XUE Ke-min, LI Ping. Coupled thermomechanical finite element analysis of metal with porosities during equal channel angular pressing process[J]. The Chinese Journal of Nonferrous Metals, 2006, 16(9): 1510?1516.

[5] BEYGELZIMER Y, VARYUKHIN V, SYNKOV S, DRLOV D.Useful properties of twist extrusion[J]. Material Science and Engineering A, 2009, 503(1/2): 14?17.

[6] BEYGELZIMER Y, RESHETOV A, SYNKOV S. Kinematics of metal flow during twist extrusion investigated with a new experimental method [J]. Journal Materials Processing Technology, 2008, 209(7): 3650?3656.

[7] VARYUKHIN V, BEYGELZIMER Y, SYNKOV S,PROKOFEVA O, KULAQIN R. Application of twist extrusion[J]. Materials Science Forum, 2006, 503/504: 335?340.

[8] AKBARIMOUSAVI S A A. Influence of strain accumulation on microstructure of aluminum 1100 in the twist extrusion[J].International Journal of Modern Physics B, 2008, 22(18/19):2858?2865.

[9] BERTA M, ORLOV D, PRANGNELL P. Grain refinement response during twist extrusion of an Al-0.13% Mg alloy[J].International Journal of Materials Research, 2007, 98:3200-3204

[10] AKBARIMOUSAVI S A A, SHAHAB A R, MASTOORI M.Computational study of Ti-6Al-4V flow behaviors during the twist extrusion process[J]. Materials and Design, 2008, 29:1316?1329.

[11] ORLOV D, BEYGELZIMER Y, SYNKOV S, et al. Evolution of microstructure and hardness in pure Al by twist extrusion[J].Materials Transactions, 2008, 49(1): 2?6.

[12] 薛克敏, 吳戰(zhàn)立, 李 萍, 王曉溪, 賈建磊. 純鋁等徑角擠扭新工藝數(shù)值模擬及實驗研究[J]. 中國機械工程, 2010, 21(12):1485?1489.XUE Ke-min, WU Zhan-li, LI Ping, WANG Xiao-xi, JIA Jian-lei.Simulate and experiment research on a new technique named equal channel angular extrusion and twist prepared by pure aluminum[J]. China Mechanical Engineering, 2010, 21(12):1485?1489.

[13] 薛克敏, 王曉溪, 李 萍. 超細(xì)晶材料制備新工藝—擠扭[J].塑性工程學(xué)報, 2009, 16(5): 130?136.XUE Ke-min, WANG Xiao-xi, LI Ping. A new technique for preparing bulk ultrafine-grained materials though twist extrusion[J]. Journal of Plasticity Engineering, 2009, 16(5):130?136.

[14] BEYGELZIMER Y, VARYUKHIN V, SYNKOV S. Shears,vortice, and mixing during twist extrusion[J]. International Journal of Material Forming, 2008, 1(s1): 443?446.

[15] 汪建敏, 周孔亢, 陸 晉, 許曉靜, 姜銀方. 層錯能在劇烈剪切變形時對影響[J]. 機械工程學(xué)報, 2008, 44(11): 126?131.WANG Jian-min, ZHOU Kong-kang, LU Jin, XU Xiao-jing,JIANG Yin-fang. Influence of stack fault energy on grain refining during sever shear deforming[J]. Chinese Journal of Mechanical Engineering, 2008, 44(12): 126?131.

Numerical simulation and experimental research on pure copper by twist extrusion

WANG Cheng, LI Ping, XUE Ke-min, LI Xiao, ZHANG Kai
(School of Materials Science and Engineering, Hefei University of Technology, Hefei 230009, China)

The finite element simulation about pure copper deformation during twist extrusion (TE) was carried out by the simulation software MSC.Marc, the procedure steps of TE, the distribution of effective strain and the change of strain component were analyzed. The experiment about pure copper at room temperature during TE was carried out using this own die. The results show that the TE deformation can be divided into there steps. The value of the effective strain in the edge area is much larger than that of the center area. The value of the normal strain is extraordinary small, while the shear strain is relative large. There are slipping and twining during the TE process. The hardness of sample significantly increase, the hardness in the edge area is visibly larger than that in the center.

pure Cu; twist extrusion deformation; numerical simulation; microstructure; microhardness

TG376

A

1004-0609(2011)12-3071-05

國家自然科學(xué)基金資助項目(50875072)；霍英東教育基金會資助項目(121053)；安徽省優(yōu)秀青年科技基金計劃(10040606Y21)

2010-09-07；

2011-10-15

李 萍，教授，博士；電話：0551-2904758；E-mail: cisi_1314@126.com

(編輯 龍懷中)