純銅棒材三輥螺旋軋制過(guò)程數(shù)值模擬

張宏軒，王瑞雪，謝君賢，程 明,3

（1.中國(guó)科學(xué)院 金屬研究所，遼寧 沈陽(yáng) 110016；2.廣東冠邦科技有限公司，廣東 佛山 528311；3.廣東精藝金屬股份有限公司，廣東 佛山 528311）

純銅棒材三輥螺旋軋制過(guò)程數(shù)值模擬

張宏軒1，王瑞雪1，謝君賢2，程 明1,3

（1.中國(guó)科學(xué)院 金屬研究所，遼寧 沈陽(yáng) 110016；2.廣東冠邦科技有限公司，廣東 佛山 528311；3.廣東精藝金屬股份有限公司，廣東 佛山 528311）

三輥螺旋軋機(jī)適合于多品種、小批量、難變形材料的小尺寸管棒材開坯。通過(guò)對(duì)三輥螺旋軋制輥形設(shè)計(jì)和軋制過(guò)程進(jìn)行研究，根據(jù)軋輥輥形計(jì)算方程，獲得三輥螺旋軋機(jī)軋輥各段輥形參數(shù)。采用Simufact. Forming軟件模擬了純銅棒材三輥螺旋軋制過(guò)程。結(jié)果顯示，純銅棒材在三輥螺旋軋制過(guò)程中的變形規(guī)律為圓形—三角形—圓形的連續(xù)過(guò)程，軋件在變形段溫度可達(dá)360℃左右，在軋輥與軋件接觸區(qū)域內(nèi)的最大應(yīng)力為590MPa，最大等效塑性應(yīng)變?yōu)?.8。軋輥在Z軸方向上所受最大載荷為43kN，而在X、Y軸方向上的最大載荷為350kN。Z方向的軸向力有助于銅棒的軸向咬入和拋出，X、Y方向的徑向力主要起到減徑和定徑的作用。

三輥螺旋軋制；輥形；純銅棒材；軋制力；數(shù)值模擬

1 引言

近二十年來(lái)，水平連鑄-行星軋制的生產(chǎn)方法在制冷空調(diào)用管的生產(chǎn)上取得了可喜的成果，生產(chǎn)效率、綜合成品率等指標(biāo)都有了較大幅度的提高，而制造成本和綜合能耗有所下降，而且水平連鑄更換合金牌號(hào)困難，適用于具有較大產(chǎn)量的規(guī)模化生產(chǎn)，而對(duì)于多牌號(hào)、小規(guī)模的產(chǎn)品可能會(huì)得不償失。此外，由于該生產(chǎn)方式生產(chǎn)設(shè)備和工藝技術(shù)的限制，所能生產(chǎn)的產(chǎn)品規(guī)格受到一定的局限。

針對(duì)于三輥行星軋機(jī)的復(fù)雜結(jié)構(gòu)與運(yùn)動(dòng)特征，三輥螺旋軋機(jī)更適合于多品種、小批量、難變形材料的小尺寸管棒材開坯。與傳統(tǒng)的擠壓法相比，加工效率和生產(chǎn)成本具有顯著優(yōu)勢(shì)。目前市場(chǎng)上對(duì)一些多規(guī)格、小批量棒材的需求，為這種軋機(jī)的推廣提供了市場(chǎng)機(jī)會(huì)。但國(guó)內(nèi)外關(guān)于三輥螺旋軋制這項(xiàng)技術(shù)的研究還不多［1］。因此，本文通過(guò)數(shù)值模擬的方法來(lái)對(duì)這一軋制過(guò)程進(jìn)行研究，探討其中的工藝規(guī)律。

2 三輥螺旋軋輥輥形設(shè)計(jì)

［2］中介紹的斜軋軋輥輥形計(jì)算方程和文獻(xiàn)［3］提出的輥形設(shè)計(jì)流程，用Matlab設(shè)計(jì)程序求解軋輥的各段輥形曲線，如圖1所示。結(jié)合文獻(xiàn)［4-5］給出的以直線相連的三輥螺旋軋制輥形的基本尺寸（螺旋角16.5°、后錐角3°、軋輥長(zhǎng)度105mm、軋輥?zhàn)畲笾睆?14mm、均整段長(zhǎng)度66mm）。通過(guò)模擬試軋，獲得以弧線相連的三輥螺旋軋制軋輥輥形尺寸，如圖2所示。這與文獻(xiàn)［6］設(shè)計(jì)的三輥斜軋機(jī)輥形類似，可用于進(jìn)一步的數(shù)值模擬計(jì)算。

圖1 軋輥各段輥形曲線

圖2 弧線相連的三輥螺旋軋輥輥形

3 三輥螺旋軋制的有限元模型建立及參數(shù)優(yōu)化

3.1 軋輥及坯料有限元模型的建立

三輥螺旋軋機(jī)模型主要由三個(gè)120°平均分布的軋輥、棒坯、喂料小車與套管構(gòu)成，如圖3所示。模型中考慮到了軋輥?zhàn)陨磔S線與棒坯軸線調(diào)整形成旋轉(zhuǎn)角α和傾斜角β，如圖4所示。

圖3 三輥螺旋軋制空間模型

圖4 旋轉(zhuǎn)角α和傾斜角β示意圖

坯料頭部采用錐形方便咬入。在Simufact中，除坯料外，模具默認(rèn)H13模具鋼，剛性體。套管在實(shí)際中不存在，模擬仿真時(shí)棒坯長(zhǎng)度遠(yuǎn)小于實(shí)際軋制長(zhǎng)度，故需對(duì)棒坯進(jìn)行徑向支撐與約束。截取棒坯長(zhǎng)度120mm，直徑50mm。軋輥的剪切摩擦摩擦系數(shù)0.95，套管摩擦系數(shù)0.01。軋輥轉(zhuǎn)速50r/min，喂料小車速度40mm/s。初始溫度20℃，工件和軋輥間熱傳導(dǎo)系數(shù)401W/（m2·K），工件和環(huán)境間熱輻射與熱對(duì)流系數(shù)0.25。軋制材料模型選取材料庫(kù)自帶的銅材料AISI為C11000。

3.2 三輥螺旋軋輥輥形參數(shù)優(yōu)化

三輥軋機(jī)軋輥的輥形和傾斜角α、螺旋角β、棒材初始直徑和成品直徑相關(guān)。

圖5為初選螺旋角β=18°［2］、棒料初始直徑50mm、不同傾斜角和斷面縮減率的模擬軋制過(guò)程，由圖可見，傾斜角α=14°斷面縮減率50%時(shí)，軋件不能軋出，其他兩種參數(shù)都可以順利軋制，考慮到效率情況，選取傾斜角α=10.5°斷面縮減率40%。

圖5 不同傾斜角和斷面縮減率的模擬結(jié)果

在傾斜角α=10.5°、斷面縮減率40%和不同螺旋角的條件下模擬軋制過(guò)程，圖6為在不同螺旋角下，分析同一橫截面上等效應(yīng)力、應(yīng)變沿徑向分布的情況，可以看出，等效應(yīng)力和等效應(yīng)變隨著螺旋角的加大而增大。雖然增大螺旋角會(huì)使變形抗力增加，但由于三輥軋制時(shí)，坯料容易出現(xiàn)心部疏松，所以適當(dāng)增大螺旋角，有助于改善軋件心部組織。表1為不同螺旋角所對(duì)應(yīng)的平均扭矩值，可以看出，軋輥的扭矩隨著螺旋角的增加而增大，三輥螺旋軋制過(guò)程中，軋件的一次軋制壓下量較大，軋機(jī)工作時(shí)承受的載荷較大，為降低軋輥扭矩和軋機(jī)功率，所以要選取相對(duì)較小的螺旋角。綜上所述，三輥螺旋軋機(jī)在設(shè)計(jì)中應(yīng)該選擇合理的螺旋角，過(guò)大會(huì)降低軋輥使用壽命，過(guò)小又不利于改善心部組織。本文根據(jù)模擬結(jié)果分析選擇螺旋角β=18°。

圖6 螺旋角對(duì)應(yīng)力、應(yīng)變的影響

表1 不同螺旋角所對(duì)應(yīng)的平均扭矩值

4 純銅棒材三輥螺旋軋制模擬結(jié)果分析

模擬的具體工況參數(shù)與上一節(jié)相同，圖7和圖8三輥螺旋軋制的等效應(yīng)力、等效應(yīng)變分布圖，可以看出，等效應(yīng)力在接觸區(qū)由外向內(nèi)逐漸減小，變形慢慢滲透到金屬的芯部，在接觸中心應(yīng)力達(dá)到最大，為590MPa（圖7）；變形區(qū)會(huì)隨著軋輥的轉(zhuǎn)動(dòng)而在圓周方向移動(dòng)，應(yīng)變?cè)茍D分布形狀與三角形相似，最大應(yīng)變?yōu)?.8（圖8）。

圖7 等效應(yīng)力分布圖

圖8 等效應(yīng)變分布圖

圖9為軋輥的軋制力變化情況，可以看出，軋輥在Z軸方向上所受最大軋制力為43kN，見圖9（a），在X、Y軸方向上的最大軋制力為350kN，見圖9（b），（c）。Z方向的軸向力有助于銅棒的軸向咬入和拋出，X、Y方向的徑向力主要起到減徑和定徑的作用，對(duì)設(shè)備剛度和軋輥材料提出一定的要求。

圖9 軋輥各個(gè)方向的軋制力變化情況

圖10為模擬軋制過(guò)程中溫度場(chǎng)的變化情況。隨著棒坯被咬入減徑后，由于變形熱以及摩擦的作用，棒坯溫度逐漸上升；當(dāng)通過(guò)集中變形段時(shí)，在軋輥輥縫的集中變形段溫度可達(dá)360℃左右，在集中變形段溫度的升高為充分的動(dòng)態(tài)再結(jié)晶的發(fā)生創(chuàng)造了重要的條件。

圖10 純銅棒材溫度場(chǎng)分布圖

圖11為不同位置截面形狀的變化情況。可以觀察到，純銅銅棒經(jīng)歷了從棒坯開始咬入、減徑、平整直到軋出的過(guò)程，棒坯的形狀從圓形—三角形－圓形的順序不斷循環(huán)。在棒坯被軋輥咬入之后，受到徑向的軋制力作用壓縮減徑，初步的形成一個(gè)接近三角形的斷面（圖11B），當(dāng)棒坯達(dá)到輥縫最小的軋輥?zhàn)冃螀^(qū)后，棒材被強(qiáng)制性減徑，此時(shí)變形最大，出現(xiàn)三角形壓扁（圖11C）。隨后經(jīng)過(guò)均整段，軋輥母線近似的平行于棒坯的軸線，此時(shí)棒坯受力均勻，橫向變形加劇，三角形壓扁逐漸消失（圖11D）。

圖11 不同位置截面形狀變化

5 結(jié)論

本文利用Matlab計(jì)算三輥行星軋機(jī)軋輥的各段輥形，通過(guò)數(shù)值模擬獲得優(yōu)化的輥形參數(shù)。通過(guò)對(duì)純銅棒材三輥螺旋軋制變形過(guò)程的模擬分析，得到以下結(jié)果：

（1）確定采用的三輥螺旋軋制軋輥傾斜角α=10.5°、螺旋角β=18°、軋輥長(zhǎng)度120mm、后錐角5°、最大直徑119mm。

（2）純銅棒材在三輥螺旋軋制過(guò)程中的變形規(guī)律為圓形－三角形－圓形的連續(xù)過(guò)程，在軋輥輥縫的集中變形段溫度可達(dá)360℃左右，超過(guò)材料的再結(jié)晶溫度，能夠有效細(xì)化晶粒，提高軋件性能。

（3）在模擬軋制過(guò)程中，在軋輥與軋件接觸的區(qū)域內(nèi)，軋件的最大應(yīng)力為590MPa，最大等效塑性應(yīng)變?yōu)?.8。軋輥在Z軸方向上所受最大載荷為43kN，而在X、Y軸方向上的最大軋制力為350kN。Z方向的軸向力有助于銅棒的軸向咬入和拋出，X、Y方向的徑向力主要起到減徑和定徑的作用，對(duì)設(shè)備剛度和軋輥材料提出一定的要求。

參考文獻(xiàn)：

［1］李臻熙, 郝夢(mèng)一. 國(guó)內(nèi)首臺(tái)中型高精度三輥螺旋軋機(jī)試制成功并通過(guò)驗(yàn)收[J]. 材料工程, 2012(1):98.

［2］朱之超, 王鋼, 趙文. 行星斜軋管機(jī)的輥形設(shè)計(jì)及實(shí)驗(yàn)研究[J]. 鋼管技術(shù), 1983(4):1-19.

［3］郭佛印, 郭琳. 基于VB語(yǔ)言的三輥行星軋機(jī)軋輥精整段曲線的設(shè)計(jì)[J]. 有色金屬加工, 2012, 41(6):39-41.

［4］張新旺. 三輥螺旋軋制成形數(shù)值模擬與實(shí)驗(yàn)研究[D]. 機(jī)械科學(xué)研究總院, 2011:15-22.

［5］陸文林, 張新旺, 吳曉煒, 等. 三輥螺旋軋機(jī)軋輥參數(shù)[J]. 塑性工程學(xué)報(bào), 2013, 20(2):40-43.

［6］李振學(xué), 姚東紅. 三輥斜軋過(guò)程數(shù)值模擬與實(shí)驗(yàn)對(duì)比研究[J]. 精密成形工程, 2011, 4(2):36-39.

Numerical Simulation on Three-roll Helical Rolling Process of Pure Copper Bar

ZHANG Hong-xuan1,WANG Rui-xue1, XIE Jun-xian2, CHENG Ming1,3
（1. Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, Liaoning, China; 2. Guangdong Guanbang Co. Ltd., Foshan 528311, Guangdong, China; 3. Guangdong Jingyi Metal Co. Ltd., Foshan 528311, Guangdong, China）

The three-roll helical rolling mill is suitable for breaking down pipe and rod that is for the material of varieties, small batch and difficult deformation. The roller profile design and rolling process of the three-roll helical rolling mill is investigated in the paper. The calculation equation of roller profile design is established and programmed by Matlab software. After the optimized parameters of the roller profile is obtained, the three-roll helical rolling process is investigated by using Simufact.Forming finite element simulation software. During the three-roll helical rolling process, the cross sections of the pure copper bar start from circle to triangle and then back to circle. And the highest temperature of the deformed billet increases to around 360℃. The maximum stress and strain in the contact area of the roller and billet is about 590MPa and 3.8 respectively. The maximum load along the Z axis of rollers is 43kN. The maximum load along the X and Y axis of rollers is around 350kN. The force in Z axis is contributed to nipping and throwing of the copper bar. The force in X and Y axisis mainly helpful to the diameter reducing and sizing of the bar.

three-roll helical rolling；roller profile；pure copper bar；rolling force；numerical simulation

TG339

A

1009-3842（2016）06-0015-04

2016-07-01

張宏軒（1993-），男，遼寧丹東人，碩士研究生，研究方向?yàn)楣馨舨能堉萍夹g(shù)。E-mail: 747839116@qq.com