王式研+陳敏+陳瑩卷

摘 要：利用有限元分析軟件ABAQUS對二輥矯直的矯直過程進行數(shù)值模擬，得到棒材在矯直過程的矯直力隨時間的變化規(guī)律，并對矯后棒材的塑性應(yīng)變進行了分析。數(shù)值分析結(jié)果中矯直力與現(xiàn)場情況基本吻合，文章仿真結(jié)果可以為實際生產(chǎn)提供指導(dǎo)。

關(guān)鍵詞：有限元；二輥矯直；塑形應(yīng)變

中圖分類號：TG333 文獻標(biāo)志碼：A 文章編號：2095-2945（2018）03-0173-02

Abstract： The finite element analysis software ABAQUS was used to simulate the straightening process of two rollers， and the rule of the bar straightening force with time was obtained. The plastic strain of the straightened bar was analyzed. The results of numerical analysis show that the straightening force is basically consistent with the field situation， and the simulation results in this paper can provide guidance for practical production.

Keywords： finite element； two-roll straightening； plastic strain

引言

隨著國內(nèi)機械行業(yè)的快速發(fā)展，人們對棒材的質(zhì)量要求越來越高，因此，對于棒材的精整顯得尤為重要[1]。二輥矯直機與其他矯直機相比具有矯直精度高、結(jié)構(gòu)簡單的優(yōu)點而被廣泛采用[2][3]。棒材彎曲過程涉及到復(fù)雜的彈塑性大撓度彎曲變形問題，有限元方法對于研究復(fù)雜問題非常有效。

本文使用大型非線性有限元軟件ABAQUS對棒材二輥矯直過程進行數(shù)值模擬，分析了矯直過程中的矯直力以及矯直后的應(yīng)力、應(yīng)變和直線度等重要參數(shù)。數(shù)值模擬結(jié)果和現(xiàn)場數(shù)據(jù)基本吻合。

1 動力顯式有限元基本理論

動力顯式有限元基本方程為

Mu+Cu+p=F （1）

式中：M-質(zhì)量矩陣；C-阻尼矩陣；p-內(nèi)力矢量；■、■-節(jié)點速度、加速度矢量。

通常采用中心差分法求解上面方程。tn時刻的速度、加速度中心差分公式分別為：

（2）

（3）

把式（2）、式（3）代入式（1），就得到求解各離散時間點的遞推公式：

顯示算法方程非耦合，可以直接求解，所有非線性包含在內(nèi)力矢量中，內(nèi)力計算是最主要的部分，不需要做收斂檢查，但是要保持穩(wěn)定狀態(tài)需要較小的時間步，是有條件穩(wěn)定的，即時間步長必須小于由該問題求解方程性質(zhì)根據(jù)Courant Friedrichs Levy穩(wěn)定準(zhǔn)則所決定的[4][5]。

2 有限元模型建立

2.1 有限元幾何模型

針對棒材二輥矯直機矯直過程，采用國內(nèi)某鋼廠提供的工藝參數(shù)，利用大型非線性有限元軟件Abaqus建立有限元模型。

在工件建模的過程中，為了縮短計算時間和分析方便，模型中工件長度只取1m，工件初始撓度8mm/m，直徑40mm。

2.2 材料模型

矯直輥、推板、擋板、入口套筒和出口導(dǎo)槽，變形極小，故設(shè)定為剛體。棒材選擇16MnCr5鋼，密度？籽=7.85×10-9·mm-3，彈性模量E=210000MPa，泊松比？滋=0.3，棒材真實應(yīng)力應(yīng)變關(guān)系如表1所示。

2.3 邊界條件和網(wǎng)格劃分

模型的邊界條件設(shè)定完全參考現(xiàn)場矯直機實際矯直情況，在初始狀態(tài)時，各矯直輥沿自己的軸線空轉(zhuǎn)，沿軋制方向分速度為25m/min。棒材與各矯直輥的摩擦系數(shù)0.2；設(shè)置棒材與擋板、入口套筒和出口導(dǎo)槽為徑向無摩擦接觸。棒材單元選擇線性六面體單元C3D8R。

3 結(jié)果分析及討論

3.1 矯直力分析

仿真計算時需確保該運動過程為理想準(zhǔn)靜態(tài)過程，在整個二輥矯直過程中上下輥的支反力隨時間變化曲線如圖1所示。從圖1中可以看出隨著棒材的逐漸咬入，矯直力逐漸增大，在0.469s處達(dá)到最大值，此時棒材端部剛剛通過矯直輥的中間段，說明此時才開始對棒材進行彎曲，而此時由于棒材具有一定初速度的沖擊作用，導(dǎo)致矯直力偏大，隨著進入穩(wěn)定矯直階段，矯直力有所降低，并且在穩(wěn)定矯直階段矯直力比較穩(wěn)定。在矯直過程的末端，隨著棒材尾部逐漸走出輥縫，矯直力也隨之減小，最大矯直力為470kN。

3.2 應(yīng)變分析

二輥矯直過程中最需要關(guān)注的是棒材的塑性應(yīng)變情況，其關(guān)乎棒材是否能夠矯直，因此對矯直過程中不同時段棒材的等效塑性應(yīng)變分布進行分析。從圖2（a）可以看出，棒材剛咬入輥縫時，發(fā)現(xiàn)只有棒材頭部與矯直輥接觸區(qū)域發(fā)生塑性變形。由于棒材具有初速度，在咬入時棒材頭部容易發(fā)生應(yīng)力集中。圖2（b）為棒材通過輥縫中段時，此時棒材頭部已經(jīng)發(fā)生較大的變形。

圖2（c）和圖2（d）分別為棒材穩(wěn)定矯直階段、棒材矯直后的等效塑性應(yīng)變分布，從中我們可以看出進入穩(wěn)定矯直階段，棒材發(fā)生更大變形，同時已經(jīng)矯直段的等效塑性應(yīng)變分布均勻；出料后棒材頭部的塑性應(yīng)變較大，尾部的塑性應(yīng)變較小，中間區(qū)域塑性應(yīng)變分布均勻。頭部塑性應(yīng)變較大是由于初速度過大造成的，而尾部變形小是由于在實際矯直過程中由于輥縫參數(shù)設(shè)置大于棒材的直徑，從而造成了棒材的兩端在通過矯直輥時未發(fā)生彎曲，該區(qū)域也是二輥矯直的盲區(qū)。

4 結(jié)束語

（1）通過有限元模型分析了棒材矯直過程中的矯直力隨時間變化規(guī)律，得出了棒材咬入時矯直力最大的結(jié)果，可以為矯直機的強度設(shè)計提供參考。

（2）分析了矯直過程及矯直后的等效塑性應(yīng)變情況，得到了棒材的塑形應(yīng)變云圖。

參考文獻：

[1]崔甫.矯直原理與矯直機械（第二版）[M].冶金工業(yè)出版社，2005，04.

[2]王云.棒材二輥矯直過程及關(guān)鍵技術(shù)研究[D].燕山大學(xué)，2013.

[3]馬自勇，馬立峰.棒材高精度二輥矯直質(zhì)量控制策略研究[J].機械工程學(xué)報，2017，03.

[4]莊茁，等.基于ABAQUS的有限元分析和應(yīng)用[M].北京：清華大學(xué)出版社，2009，1.

[5]尹紀(jì)龍，李大永，等.動力有限元法輥彎成形全流程仿真技術(shù)研究[J].塑性工程學(xué)報，2006，13（1）：36-39.

[6]王式研，馬立東，陳敏.斜輥鋼管矯直過程的數(shù)值模擬與接觸狀態(tài)分析[J].科技創(chuàng)新與應(yīng)用，2016（35）：68-69.

[7]孫維陽，王孝海.輥式型鋼矯直機的機械結(jié)構(gòu)分析[J].科技創(chuàng)新與應(yīng)用，2014（31）：107.endprint