孫會(huì)朝,周俊陵,許榮昌,王文生,康永林

(1.萊蕪鋼鐵集團(tuán)有限公司技術(shù)研發(fā)中心,山東萊蕪,271104;2.北京科技大學(xué)材料科學(xué)與工程學(xué)院,北京,100083)

隨著我國H型鋼需求量的不斷增長,國內(nèi)H型鋼產(chǎn)量也在不斷增加,H型鋼的品種規(guī)格得到進(jìn)一步完善。隨著H型鋼軋制理論研究不斷深化,國內(nèi)眾多研究者采用有限元方法對H型鋼軋制過程進(jìn)行研究。朱國明等[1-4]采用彈塑性熱力耦合有限元方法,完成大型H型鋼全軋程三維熱力耦合仿真分析,并在全軋程仿真結(jié)果的基礎(chǔ)上,對大型H型鋼軋后殘余應(yīng)力進(jìn)行了仿真分析,結(jié)合仿真結(jié)果提出了控制大型H型鋼殘余應(yīng)力的方法和手段;奚鐵等[5-6]借助有限元分析軟件SuperForm對H型鋼開坯軋制和萬能軋制過程進(jìn)行了模擬;徐旭東等[7-8]利用有限元分析方法對H型鋼軋制和軋后冷卻過程進(jìn)行二維溫度場的模擬,同時(shí)采用顯式動(dòng)力學(xué)有限元分析的方法,模擬了不同變形參數(shù)下H型鋼的萬能軋制過程。

H型鋼的斷面比較復(fù)雜,軋制過程中軋件內(nèi)的力場對產(chǎn)品的性能起決定性作用,力場直接決定成型過程,同時(shí)在力的作用下產(chǎn)生的變形功又會(huì)迅速引起溫度的變化。為此,本文利用ANSYS有限元軟件,以大型H型鋼生產(chǎn)線為基礎(chǔ),采用熱力耦合彈塑性有限元方法建立萬能軋制有限元仿真分析模型,并對不同摩擦系數(shù)下的多種工況進(jìn)行了仿真分析,以期為H型鋼萬能軋制提供理論依據(jù)。

1 模擬計(jì)算條件

選擇HN800 mm×300 mm連軋過程的第4連軋道次的UF道次作為研究道次。圖1為UF軋機(jī)孔型圖和軋件尺寸示意圖。表1為連軋道次的主要軋制條件。該道次軋制仿真分析的初始溫度場通過全軋程熱力耦合的方法得到,圖2為軋件斷面溫度場的分布云圖。

圖1 UF軋機(jī)孔型和軋件尺寸(mm)示意圖Fig.1 Sketch of the pass system and stock dimension of UF rollingm ill

表1 連軋道次的主要計(jì)算邊界條件Table 1 Main calculation border conditions for continuous rolling process

圖2 軋件斷面溫度場分布云圖Fig.2 Temperature distribution of stock section

2 仿真分析

2.1 工況

在H型鋼的仿真分析過程中采用庫侖摩擦,在軋制仿真過程中材料的摩擦系數(shù)對軋制力、金屬流動(dòng)等均會(huì)產(chǎn)生影響。在熱軋過程中材料的摩擦系數(shù)一般定為0.35左右。為了針對摩擦系數(shù)對H型鋼軋制過程的影響進(jìn)行研究,特人為改變摩擦系數(shù),確定的不同工況如表2所示。

表2 不同工況下的摩擦系數(shù)Table 2 Friction coefficient of each model

2.2 X方向金屬的流動(dòng)

在建立有限元模型過程中,X為腹板寬度方向,Y為翼緣寬展方向,Z為軋制方向。圖3為調(diào)整摩擦系數(shù)后各工況X方向位移的計(jì)算結(jié)果。由圖3可看出,摩擦系數(shù)的改變對H型鋼穩(wěn)定軋制階段X方向的位移影響不大,其斷面整體位移分布狀態(tài)基本不變。

2.3 Y方向翼緣寬展仿真

圖4為各工況Y方向穩(wěn)定軋制階段斷面節(jié)點(diǎn)位移云圖。由圖4可看出,摩擦系數(shù)的改變對H型鋼斷面Y方向的位移分布趨勢影響不大,但其寬展量發(fā)生明顯變化。按照翼緣寬展量和摩擦系數(shù)的關(guān)系作曲線,不難看出,隨著摩擦系數(shù)的增加,翼緣寬展量降低,且呈線性關(guān)系。

圖3 X方向位移云圖Fig.3 Contours of X-displacement

圖4 Y方向位移云圖Fig.4 Contours of Y-displacement

摩擦系數(shù)的增加,導(dǎo)致翼緣寬展的增大,這一點(diǎn)可以通過翼緣內(nèi)側(cè)金屬與水平輥端部側(cè)面相接觸進(jìn)行分析。在變形區(qū)內(nèi)的金屬,其翼緣內(nèi)側(cè)與水平輥端部側(cè)面存在接觸,水平輥的旋轉(zhuǎn)會(huì)造成翼緣部位金屬向R角方向流動(dòng),并導(dǎo)致腿部寬展量減少。而摩擦系數(shù)越大,金屬向R角方向流動(dòng)越強(qiáng),其寬展量也就越小。同時(shí),在R角部位的翼緣內(nèi)側(cè),隨著摩擦系數(shù)的增加,向R角部位的流動(dòng)越大。圖5為翼緣寬展量和摩擦系數(shù)之間的關(guān)系曲線。

圖5 摩擦系數(shù)與翼緣寬展量的關(guān)系Fig.5 Relationship between friction coefficient and flange spreading

2.4 前后滑

圖6為不同摩擦系數(shù)下穩(wěn)定軋制階段變形區(qū)內(nèi)部金屬在軋制方向節(jié)點(diǎn)速度的等值面云圖。由圖6可看出,隨著摩擦系數(shù)的逐漸增加,前滑區(qū)域逐漸增大。同時(shí),隨著摩擦系數(shù)的減小,前滑區(qū)首先從腹板中心部位開始減小,并向腿腰連接部位逐漸收縮。因此,H型鋼腹板和翼緣是一個(gè)整體,不能單獨(dú)考慮。

圖6 變形區(qū)內(nèi)軋制方向節(jié)點(diǎn)速度等值面云圖Fig.6 Isosurface contours of node velocity along rolling direction

2.5 軋制力仿真

圖7 不同摩擦系數(shù)下水平輥軋制力的變化曲線

Fig.7 Rolling force curve of horizontal roller

圖8 不同摩擦系數(shù)下立輥軋制力的變化曲線Fig.8 Rolling force curve of vertical roller

圖9 摩擦系數(shù)對水平輥軋制力的影響Fig.9 Effect of friction coefficient on rolling force of horizontal roller

提取不同摩擦系數(shù)下的軋制力進(jìn)行仿真,仿真結(jié)果如圖7、圖8所示。由圖7、圖8可看出,隨著摩擦系數(shù)的增加,無論水平輥軋制力還是立輥軋制力,均明顯增強(qiáng),但軋制力曲線形狀仍舊保持原來的狀態(tài)。圖9、圖10為不同摩擦系數(shù)下穩(wěn)定軋制階段水平輥和立輥軋制力的變化曲線。由圖9、圖10可看出,穩(wěn)定軋制階段水平輥和立輥軋制力與摩擦系數(shù)呈正比,即摩擦系數(shù)越大,穩(wěn)定軋制階段水平輥和立輥軋制力越大。

圖10 摩擦系數(shù)對立輥軋制力的影響Fig.10 Effect of friction coefficien t on rolling force of vertical roller

3 結(jié)論

(1)摩擦系數(shù)的改變對H型鋼斷面Y方向的位移分布趨勢影響不大,但其寬展量發(fā)生明顯變化。

(2)隨著摩擦系數(shù)的逐漸增加,前滑區(qū)域逐漸增大。

(3)隨著摩擦系數(shù)的增加,無論水平輥軋制力還是立輥軋制力,均明顯增強(qiáng)。

[1] Zhu Guoming,Kang Yonglin,Chen Wei.3D thermalmechanical coup led elasto-p lastic finite element analysis in the w hole rolling p rocessof H-beam[J].Materials Science Fo rum,2008,575:532-538.

[2] 朱國明,康永林,陳偉.H型鋼多道次可逆開坯軋制過程的三維熱力耦合仿真分析[J].中國機(jī)械工程,2007,18(14):1 747-1 751.

[3] 朱國明,康永林,陳偉,等.H型鋼空冷過程中殘余熱應(yīng)力的有限元分析[J].機(jī)械工程材料,2008,32(4):77-80.

[4] 朱國明,康永林,陳偉,等.降低H型鋼殘余熱應(yīng)力的三維有限元仿真分析[J].軋鋼,2007,24(4):22-26.

[5] 奚鐵,錢奕峰,章靜.H型鋼開坯軋制變形分析[J].軋鋼,2004,21(6):47-49.

[6] 曹杰,奚鐵,章靜,等.H型鋼萬能軋制變形分析[J].重型機(jī)械,2005(1):23-25.

[7] 徐旭東,王秉新,劉相華,等.H型鋼控制冷卻的有限元模擬[J].鋼鐵研究學(xué)報(bào),2005,17(2):30-33.

[8] 徐旭東,王秉新,劉相華,等.H型鋼軋制力的數(shù)值模擬分析[J].鋼鐵,2005,40(7):56-59.