，

(沈陽建筑大學(xué)交通與機械工程學(xué)院，遼寧 沈陽 110168)

壓下量對鋼筋矯直系統(tǒng)力學(xué)性能的影響分析

鄭夕健，朱禹霏

(沈陽建筑大學(xué)交通與機械工程學(xué)院，遼寧 沈陽 110168)

0 引言

鋼筋在建筑行業(yè)的工程應(yīng)用非常普遍。鋼筋在軋制、冷卻、運輸過程中由于散熱不均、殘余應(yīng)力和外力沖擊等因素影響，會出現(xiàn)不同程度的彎曲，所以鋼筋在加工使用之前必須對其進行矯直[1]。在矯直鋼筋過程中，壓下量直接決定矯直質(zhì)量的好壞，又影響矯直效率，因此，壓下量對矯直質(zhì)量的影響一直是人們重點研究的課題[2]。

鋼筋矯直過程中的彎曲變形具有壓下量小、彈塑性變形顯著和載荷分布多樣化等特點，壓下量對矯直變形的影響一般除矯直力可以被測量外，變形區(qū)的應(yīng)力、應(yīng)變大小則很難通過實驗或計算得到，造成了認識矯直機理、正確掌握矯直變形規(guī)律的困難。

1 矯直系統(tǒng)基本參數(shù)的設(shè)定

1.1 壓下量的設(shè)定

隨著矯直工作的進行，鋼筋將與矯直輥逐個接觸，鋼筋會受到矯直輥作用，改變鋼筋原始彎曲曲率，直至殘余曲率符合平直度的要求為止。上排各矯直輥壓下量將直接影響矯直過程中與矯直后鋼筋內(nèi)部的應(yīng)力、應(yīng)變的變化規(guī)律及分布情況，最終影響矯直后鋼筋的質(zhì)量。GT6-12型號矯直機的計算壓下量值和實測壓下量值如表1所示。其中，l=60時為理論值，l=85時為實測值[3]。

表1 上輥壓下量 mm

鋼筋直徑Φ8Φ10Φ12上輥壓下量δ1#l=60000l=850002#l=601．941．581．34l=853．93．22．63#l=601．71．781．15l=853．42．72．34#l=601．581．291．07l=853．22．52．15#l=601．471．21l=8532．41．86#l=601．2810．86l=852．621．5

1.2 矯直系統(tǒng)基本參數(shù)及鋼筋材料屬性

矯直輥按照上、下平行交錯方式進行配置，在矯直鋼筋時，通過上輥下壓對鋼筋進行重復(fù)塑性彎曲以達到矯直目的。選用直徑為8 mm、10mm、12mm、長度為1m、牌號為HRB335，即屈服強度為335 MPa的鋼筋進行仿真。矯直系統(tǒng)基本參數(shù)及鋼筋材料屬性如表2所示。

表2 基本數(shù)據(jù)及材料屬性

名稱參數(shù)單位輥徑100mm輥數(shù)(上6下5)11個輥距60，85mm矯直速度2．88mm/s彈性模量2．03e11Pa鋼筋屈服極限3．35e8Pa鋼筋直徑8～12mm泊松比0．3

2 有限元模型的建立

矯直輥與鋼筋均采用8節(jié)點solid164單元，該單元是一種8節(jié)點實體單元，不需定義實常數(shù)，只用在動力顯示分析中，支持所有許可的非線性特性[4]。solid164單元有2種算法：一是單點積分法，該方法對大變形問題非常有效，但缺點是有沙漏問題；二是完全積分法，該方法無沙漏問題，但求解慢，若泊松比比較大時要謹慎使用。

因此，選擇單點積分法，以便節(jié)省時間，只是需要控制沙漏問題[5]。由于矯直輥工作時幾乎不變形，所以在有限元建模時把矯直輥設(shè)為剛性體，這樣可以大大減少顯示分析的計算時間，因為無論定義了多少個節(jié)點，剛性體都只有6個自由度，剛性體內(nèi)所有節(jié)點的自由度都會自動耦合到剛性體的質(zhì)心上[6]。鋼筋為變形體，選用雙向性隨動強化材料模型，該模型可以使用2種斜率(彈性和塑性)來表示材料的應(yīng)力應(yīng)變行為。

由于矯直輥輥型比較復(fù)雜，因此，采用自由網(wǎng)格劃分的方式。因為研究對象是鋼筋，所以鋼筋的網(wǎng)格劃分盡量密一些，矯直輥的可以稀疏一些。劃分完成后，每個矯直輥單元數(shù)為16 000個，Φ8 mm、Φ10mm、Φ12mm鋼筋單元數(shù)分別為19 456,21952和25 434。

模型中，矯直輥繞中心軸轉(zhuǎn)動，所以除了繞中心軸轉(zhuǎn)動的約束釋放外，其余自由度全部約束。由于鋼筋矯直時與矯直輥輥槽有接觸，接觸時矯直力會使鋼筋表面有微小變形，使鋼筋表面與輥槽面面接觸，所以接觸類型選擇面面接觸，選擇此類型中的自動接觸類型，需定義Contact面和Target面[7]。矯直過程中鋼筋會發(fā)生彈塑性變形，將其定義為Target面，將矯直輥定義為Contact面，輥子與鋼筋之間的靜摩擦系數(shù)取0.3，動摩擦系數(shù)取0.2[8]。

創(chuàng)建原始曲率為零的鋼筋進行模擬仿真，計算結(jié)束后，在ANSYS后處理器中直接查看鋼筋在矯直過程中和矯后的應(yīng)力、應(yīng)變特征。

3 有限元仿真分析

3.1 壓下量對鋼筋應(yīng)力分布的影響

根據(jù)表2設(shè)置的矯直參數(shù)進行仿真，得到矯直過程中和矯直后的應(yīng)力等值線圖。其中Φ8鋼筋的應(yīng)力等值線如圖1,圖2所示。

圖1 l=60時Φ8鋼筋的應(yīng)力等值線

圖2 l=85時Φ8鋼筋的應(yīng)力等值線

由圖1，圖2可以看出，鋼筋在l=60的壓下量下矯直時，應(yīng)力主要分布在173～693MPa之間，矯直后應(yīng)力主要分布在150～595 MPa之間；鋼筋在l=85的壓下量下矯直時，應(yīng)力主要分布在214～856 MPa之間，矯直后應(yīng)力主要分布在232～919 MPa之間。在l=60的壓下量下矯直后，鋼筋的應(yīng)力值要小于矯直時的應(yīng)力值，而在l=85的情況下恰好相反。同時，還發(fā)現(xiàn)無論在矯直時還是在矯直后，l=60時的應(yīng)力值都小于l=85時的應(yīng)力值，應(yīng)力分布也比較均勻。

為了研究不同直徑鋼筋的應(yīng)力變化情況，列出了Φ10,Φ12鋼筋的應(yīng)力等值線，如圖3～圖6所示。

圖3 l=60時Φ10鋼筋的應(yīng)力等值線

圖4 l=85時Φ10鋼筋的應(yīng)力等值線

圖5 l=60時Φ12鋼筋的應(yīng)力等值線

圖6 l=85時Φ12鋼筋的應(yīng)力等值線

從圖3～圖6可以明顯的看到，粗鋼筋矯直后應(yīng)力值明顯低于細鋼筋的。同一直徑的鋼筋在l=60的壓下量下矯直時，應(yīng)力值比l=85時的應(yīng)力小，且在l=60時矯直后的應(yīng)力值比矯直過程中的應(yīng)力值有所減小，這一點對鋼筋的倉儲、運輸及后期的使用有重要意義，可以有效避免鋼筋放置后由于內(nèi)部應(yīng)力不平衡而再次變形。通過觀察矯后等值線圖還發(fā)現(xiàn)，同一直徑鋼筋在l=60時，應(yīng)力值明顯低于l=85的應(yīng)力值，且分布比較均勻，說明理論壓下量的矯直方案更加合理。

3.2 壓下量與彈塑性應(yīng)變的關(guān)系

為研究壓下量與彈塑性應(yīng)變的關(guān)系，以l=60時的壓下量，對3種直徑鋼筋兩端面各選取靠近中心的一節(jié)點和應(yīng)力最大處的節(jié)點進行彈塑性變化分析。圖7,圖8為Φ8鋼筋應(yīng)變隨時間變化的曲線。

圖7 l=60時Φ8鋼筋所選節(jié)點的彈性應(yīng)變隨時間變化曲線

由圖7～圖8可以看出，鋼筋在矯直過程中，所選端面節(jié)點的彈塑性應(yīng)變很小，可以忽略不計，即中性軸處既不發(fā)生彈性應(yīng)變也不發(fā)生塑性應(yīng)變，說明鋼筋矯直時中心軸位置不發(fā)生變形，與中性層不變形理論具有一致性，說明應(yīng)用ANSYS/LS-DYNA軟件對鋼筋矯直進行模擬是可行的。

圖8 l=60時Φ8鋼筋所選節(jié)點的塑性應(yīng)變隨時間變化曲線

由應(yīng)變線圖還可以看出，塑性應(yīng)變隨時間變化呈現(xiàn)梯形遞增趨勢，這是因為矯直時，鋼筋發(fā)生塑性變化后，還會在矯直輥的作用下繼續(xù)被加載及卸載后又重新加載，屈服極限發(fā)生了變化，出現(xiàn)冷作硬化現(xiàn)象，以新的強度進入下一矯直區(qū)，進而發(fā)生新的塑性變形，直至鋼筋走出矯直區(qū)。鋼筋矯直時最大應(yīng)變位置靠近表層金屬，過大的塑性變形會被釋放，所以在l=60時的壓下量下矯直時，能夠有效防止鋼筋表面開裂現(xiàn)象發(fā)生。

4 結(jié)束語

通過對矯直系統(tǒng)有限元模型的仿真分析，可以得到如下結(jié)論：

a.應(yīng)用ANSYS/LS-DYNA軟件對鋼筋矯直進行有限元動力學(xué)分析，該方法可以直觀地對比分析在矯直時和矯直后鋼筋內(nèi)部應(yīng)力應(yīng)變的分布情況及變化規(guī)律。通過分析發(fā)現(xiàn)，矯直時應(yīng)力主要分布在鋼筋表層位置；彈性應(yīng)變無變化規(guī)律，但塑性應(yīng)變隨著鋼筋前進而呈現(xiàn)梯形遞增趨勢，且彈性應(yīng)變遠小于塑性應(yīng)變；矯直時，鋼筋越粗，其應(yīng)力值越小，說明粗鋼筋矯直時所用矯直力越小，更容易被矯直。

b.Φ8鋼筋在理論壓下量下矯直時，應(yīng)力主要分布分別在173～693MPa之間，矯直后應(yīng)力主要分布在150～595 MPa之間；在實測的壓下量下矯直時，應(yīng)力主要分布分別在214～856 MPa之間，矯直后應(yīng)力主要分布在232～919 MPa之間，在理論壓下量下矯直后，鋼筋的應(yīng)力值要小于矯直時的應(yīng)力值，而實測壓下量的情況下恰好相反，觀察Φ10、Φ12鋼筋的應(yīng)力等值線圖也符合此現(xiàn)象。無論在矯直時還是在矯直后，理論壓下量下矯直的鋼筋應(yīng)力值比較小，分布也更加均勻，說明理論壓下量下的矯直方案更加合理。

c.在中性軸附近，彈塑性應(yīng)變在0～4×10-3之間變化，應(yīng)變值很小，可以認為中性軸應(yīng)變?yōu)榱?，即中性軸不發(fā)生變形；在應(yīng)力最大結(jié)點處，其塑性應(yīng)變遠大于彈性應(yīng)變，觀察不同直徑鋼筋塑性應(yīng)變曲線發(fā)現(xiàn)，無論鋼筋粗細，矯直時最大塑性應(yīng)變出現(xiàn)在靠近表層位置，過大的塑性變形會被釋放，在理論壓下量下矯直時，可以有效防止鋼筋表面發(fā)生開裂現(xiàn)象，能夠滿足對鋼筋表面質(zhì)量的要求。

[1] 崔 甫.矯直原理與矯直機械[M].北京：冶金工業(yè)出版社，2005.

[2] Zhao Binli，Yang Baolai，Zhang Banhua.Reduction of the I-beam straightening the quality of finite element analysis [J].Machinery Design and Manufacture, 2010(6):6-7.

[3] Zheng Xijian,Ma Lei,Zhu Yufei.Analysis of influence of the roller number and rolling schedule to straighten quality of the project steel bar [C]//Applied Mechanics and Materials，2013.

[4] 劉 念.輥式矯直機矯直過程有限元仿真[D].武漢：武漢科技大學(xué)，2008.

[5] 楊寶來.型鋼矯直機理的數(shù)值模擬及實驗研究[D]. 秦皇島:燕山大學(xué)，2010.

[6] 孟欣佳,張立香. 鋼筋彎曲過程的有限元模擬[J].河北工程大學(xué)學(xué)報(自然科學(xué)版)，2011，28(2)：30-32,37.

[7] 周存龍,王國棟,劉相華,等.壓彎量對軋件矯直變形影響的有限元分析[J]．塑性工程學(xué)報，2006(1)：78-81.

[8] 尚曉江,蘇建寧.ANSYS/LS-DYNA動力分析方法與工程實例[M]．北京：中國水利水電出版社，2006.

Analysis of the Influence of Rolling Schedule on Mechanical Characteristicsof Straightening System of Steel Bars

ZHENGXijian,ZHUYufei

(School of Transportation and Mechanical Engineering, Shenyang Jianzhu University, Shenyang 110168，China)

以十一輥平行輥矯直機為原型，運用有限元分析軟件ANSYS/LS-DYNA建立了鋼筋矯直系統(tǒng)有限元模型；在理論壓下量和實測壓下量下，分別對3種不同直徑的鋼筋進行了動力學(xué)仿真分析，得到了鋼筋在矯直過程中與矯直后的應(yīng)力大小及分布情況；分析并比較了在理論壓下量和實測壓下量下，端面中心節(jié)點處與應(yīng)力最大點處的彈塑性應(yīng)變變化規(guī)律，找到了相對合理的壓下量方案。該研究結(jié)果對掌握鋼筋矯直過程的力學(xué)特性變化和提高矯直質(zhì)量具有一定的參考價值和實際意義。

鋼筋矯直；壓下量；應(yīng)力；應(yīng)變

Taking eleven-roll parallel roll straightening machine as a prototype, the authorsestablished a FEM of straightening system of steel bars by ANSYS/LS-DYNA； under theoretical and measured rolling reduction, the authors simulated the process for straightening steel bars of three different kinds of diameter and did kinetic analysis, as a result, the authors obtained stress and its distribution in the process of straightening and after straightening； through analyzing and comparing the change rules of elastic and plastic strain at the center node of each face and at the biggest stress point under theoretical and measured rolling reduction, the authors found a comparatively reasonable rolling schedule. The study provides certain references for mastering mechanical characteristics of straightening process and improving the quality of straightening and it has practical significance.

straightening steel bars； rolling reduction； stress； strain

2014-05-07

“十二五”國家科技支撐計劃資助項目(2011BAJ02B03-3)

TH12

A

1001-2257(2014)09-0015-04

鄭夕健(1963-)，男，山東容城人，教授，研究方向為建筑機械設(shè)計理論、鋼結(jié)構(gòu)穩(wěn)定理論、建筑結(jié)構(gòu)動力學(xué)分析、有限元結(jié)構(gòu)分析和建筑機械產(chǎn)品開發(fā)等；朱禹霏(1990-)，女，遼寧沈陽人，碩士研究生，研究方向為工程機械現(xiàn)代設(shè)計方法。