冷軋鋁帶卷取內(nèi)部應(yīng)力有限差分仿真研究

張波，王偉,2，李世苗

（1. 福州大學(xué) 機(jī)械工程及自動化學(xué)院，福州 350100；2. 福建省高端裝備制造協(xié)同創(chuàng)新中心，福州 350100）

冷軋鋁帶卷取內(nèi)部應(yīng)力有限差分仿真研究

張波1，王偉1,2，李世苗1

（1. 福州大學(xué) 機(jī)械工程及自動化學(xué)院，福州 350100；2. 福建省高端裝備制造協(xié)同創(chuàng)新中心，福州 350100）

目的 研究冷軋鋁帶卷取內(nèi)部應(yīng)力的分布情況以及影響因素，為鋁帶卷筒的設(shè)計(jì)提供參考。方法 通過冷軋鋁板壓縮實(shí)驗(yàn)，測定鋁帶徑向壓縮系數(shù)與應(yīng)力的關(guān)系曲線；采用非線性有限差分法對卷取過程和內(nèi)部應(yīng)力進(jìn)行仿真；分析線性、一次非線性和二次非線性材料模型以及卷筒剛度對帶卷內(nèi)部應(yīng)力分布的影響。結(jié)果 徑向壓縮系數(shù)的理論值與實(shí)驗(yàn)值有較好的一致性，卷筒剛度對鋁帶卷內(nèi)層應(yīng)力影響較大，而對外層影響很小，卷筒剛度較小時，在卷芯處容易產(chǎn)生塌卷缺陷，卷筒剛度較大時，容易出現(xiàn)扁卷現(xiàn)象。結(jié)論 計(jì)算帶卷內(nèi)部應(yīng)力分布時，需要考慮鋁帶卷的各向異性以及層與層之間的緊密程度，合理選擇卷筒的剛度對提高卷取質(zhì)量有很大的實(shí)際意義。

冷軋鋁帶；有限差分法；徑向壓縮系數(shù)

鋁帶卷取是帶逐層卷到卷筒上的過程，涉及帶與卷筒之間以及帶層之間的接觸[1—2]。卷取過程不僅會在卷筒上施加很大的壓力，影響卷筒的使用；而且隨著卷取層數(shù)的增加，帶卷內(nèi)部會發(fā)生復(fù)雜的應(yīng)力應(yīng)變變化。帶卷內(nèi)部應(yīng)力分布不合理會造成卷芯塌陷、扁卷、層間滑移等缺陷[3—6]，因此，鋁帶卷取過程中內(nèi)部應(yīng)力的研究不但能揭示卷取過程中的力學(xué)現(xiàn)象，而且將為鋁帶卷筒的選擇提供參考。

鋁帶的卷取過程是一個高度非線性的過程，涉及材料非線性、幾何非線性以及接觸非線性[7]。對工程上大多數(shù)非線性問題，很少可直接利用解析法獲得其解，少數(shù)情況可利用簡化假設(shè)來求解，但是過多的假設(shè)可能導(dǎo)致結(jié)果不正確甚至錯誤的解，因此數(shù)值解是一個好的替代方法，有限差分和有限元方法是兩種常用的數(shù)值計(jì)算方法[8—11]。燕山大學(xué)的李海亮利用有限元軟件ANSYS/LS-DYNA，建立了冷軋卷取過程實(shí)體模型，對卷取過程進(jìn)行動態(tài)模擬仿真分析，得出在卷取過程中張力分布規(guī)律[12]。常鐵柱等人在現(xiàn)場實(shí)測緊密系數(shù)實(shí)驗(yàn)的基礎(chǔ)上，通過卷取過程應(yīng)力場的差分求解以及卷取結(jié)束后卸卷過程中機(jī)械應(yīng)力和熱應(yīng)力耦合應(yīng)力場的差分法求解，獲得了鋼卷內(nèi)部應(yīng)力場的分布規(guī)律[13]。

為研究鋁帶卷取過程中帶卷內(nèi)部應(yīng)力分布情況以及影響因素。文中利用有限差分法對鋁帶卷取過程中帶卷內(nèi)部應(yīng)力進(jìn)行了數(shù)值仿真，分析了線性、一次非線性和二次非線性材料模型以及卷筒剛度對帶卷內(nèi)部應(yīng)力分布的影響。

1 冷軋鋁帶徑向壓縮系數(shù)

鋁帶卷曲的徑向變形除了一般的彈性形變外，還應(yīng)考慮層間粗糙面接觸引起的附加壓縮變形，所以引入了鋁帶徑向壓縮系數(shù)概念。緊密系數(shù)m是材料彈性模量E與徑向壓縮系數(shù)Er的比值，其值取決于帶材表面的粗糙度、壓力以及厚度，所以確定m就需要知道徑向壓縮系數(shù)Er。下面通過實(shí)驗(yàn)和理論對Er的值進(jìn)行研究。

鋁帶徑向壓縮系數(shù)的測定采用電子萬能試驗(yàn)機(jī)CMT5305，最大的試驗(yàn)力為300 kN，精度為1級。壓板采用淬火軸承鋼Φ150 mm，試驗(yàn)對象取某廠0.38 mm規(guī)格的1050鋁合金，尺寸為130 mm×130 mm，取 160層。通過試驗(yàn)得到的數(shù)據(jù)為壓力-位移的關(guān)系曲線，轉(zhuǎn)換為應(yīng)力-應(yīng)變曲線，其斜率為徑向壓縮系數(shù)Er，由此可得到層疊板的徑向壓縮系數(shù)與應(yīng)力之間的關(guān)系。

根據(jù)文獻(xiàn)[14]所述，徑向壓縮系數(shù)的計(jì)算公式為：

式中：h為鋁帶的平均厚度；k0為修正系數(shù)，一般取0.3～0.6，這里取為0.45；Δ為鋁帶材接觸表面不平度的最大值，這里取為5.87 μm；p為加載壓力；E為鋁帶的彈性模量。

由式(1)可知，徑向壓縮系數(shù)Er的值與壓力和厚度成正比關(guān)系，與表面不平度成反比。將板疊壓縮試驗(yàn)中的壓力實(shí)測值代入式(1)，可得徑向壓縮系數(shù)與應(yīng)力的關(guān)系，徑向壓縮系數(shù)的理論值與實(shí)驗(yàn)實(shí)測的對比見圖1，可以看出，兩者有較好的一致性，但存在一定誤差。

圖1 徑向壓縮系數(shù)的理論值與實(shí)驗(yàn)實(shí)測的對比Fig.1 Comparison between theoretical value of radial compression coefficient and experimental measurement

2 鋁帶卷取平面軸對稱有限差分模型

2.1 鋁帶卷內(nèi)部應(yīng)力解析模型

根據(jù)卷取機(jī)的實(shí)際工作情況來看，鋁帶是逐層纏繞在卷筒上的，分析模型見圖2。當(dāng)卷取完成時，將每層帶材的徑向壓力增量在徑向方向進(jìn)行疊加，即可得到鋁帶卷層間壓力隨徑向的變化，同樣，求得每層帶材切向應(yīng)力增量后，也可求得切向應(yīng)力沿徑向分布。在卷取半徑為r的單層鋁帶時產(chǎn)生的壓力稱為壓力增量，用δp來表示。

圖2 卷取解析模型Fig.2 Coiling analysis model

由文獻(xiàn)[15]知，用壓力增量表示的微分方程為：

壓力增量滿足兩個邊界條件：① 最外層壓力增量與卷取張力滿足式(3)；② 鋁帶卷內(nèi)表面徑向位移滿足連續(xù)性條件，即帶卷內(nèi)部徑向位移等于卷筒徑向位移，即式(4)。

式中：rmax為鋁帶卷外徑；Ec為卷筒徑向剛度；r0為鋁帶卷內(nèi)表面半徑，即為卷筒外表面半徑；δT為帶卷內(nèi)部切向應(yīng)力增量；ν為鋁帶泊松比；Tw為張力；Eθ為鋁帶切向彈性模量。

根據(jù)不考慮剪切時用極坐標(biāo)表示的平衡方程，可得到每層鋁帶切向應(yīng)力增量方程式(5)。

式(2), (3)和(4)組成了卷取單層鋁帶時求解徑向壓力增量的邊值求解模型，求解出鋁帶卷每層的徑向壓力增量后，再根據(jù)式(5)可得到切向應(yīng)力增量。

2.2 有限差分法列式

假設(shè)鋁帶卷共有N層，則卷取第i層時的內(nèi)表面半徑ri為在半徑ri處，由卷取N+1層引起的壓力增量用δpi表示。將徑向壓力的一階與二階基本差分方程代入式(2)中，得徑向壓力增量的差分方程(6)。

邊界條件的差分方程分別為：

求解方程組可得每層鋁帶的徑向壓力增量，利用以下疊加公式得到沿半徑方向分布的鋁帶卷內(nèi)部徑向壓力與切向應(yīng)力。

式中：δpij為第i層由于第j層卷取產(chǎn)生的徑向壓力增量；δTij為第i層由于第j層產(chǎn)生的張力增量。

2.3 鋁帶卷材料模型和有限差分求解

鋁帶卷為各向異性材料，徑向壓縮系數(shù)由于卷取的緊密程度不同而與切向彈性模量不相等。根據(jù)實(shí)驗(yàn)所得的徑向壓縮系數(shù)與壓力的關(guān)系曲線，可得到一次和二次擬合方程式。

某四輥鋁帶冷軋卷取機(jī)Ec=61 MPa/mm，卷筒外徑為r0=332 mm。軋件設(shè)定張力Tw=35 MPa，鋁帶板厚h=0.72 mm，切向彈性模量Eθ=77 GPa，泊松比ν=0.33，共卷取1000層。有限差分計(jì)算流程見圖3，利用MATLAB[16—17]編制計(jì)算程序。

圖3 差分法計(jì)算鋁帶卷內(nèi)部應(yīng)力流程Fig.3 Flow chart for calculating internal stress of aluminum coil with difference method

3 鋁帶卷內(nèi)部徑向壓力和切向應(yīng)力結(jié)果分析

3.1 鋁帶卷內(nèi)部應(yīng)力分布規(guī)律

由以上工況計(jì)算鋁帶卷內(nèi)部應(yīng)力分布情況，圖4為考慮徑向壓縮系數(shù)Er為一次擬合時的帶卷內(nèi)部徑向壓力與切向應(yīng)力的分布，可得到如下規(guī)律。

帶卷內(nèi)層徑向壓力最大，層間徑向壓力隨卷徑增大而減小，中部區(qū)域變化較為平緩，最外層為自由表面，其徑向壓力為0，在鋁帶卷內(nèi)部區(qū)域與外部區(qū)域徑向壓力沿半徑方向下降較快，而中部區(qū)域變化較為平緩。帶卷內(nèi)部切向應(yīng)力較大，隨后減小，在外層一定區(qū)域開始再次增大，最外層切向應(yīng)力與卷取張應(yīng)力相等，在帶卷中部區(qū)域切向應(yīng)力出現(xiàn)負(fù)值，表明帶卷出現(xiàn)張力松弛，這與實(shí)際經(jīng)驗(yàn)相符。從圖4還可看出帶卷內(nèi)層徑向壓力最大，對卷筒強(qiáng)度影響較大；外部區(qū)域徑向壓力較小，開卷時有滑移的趨勢。

3.2 材料模型對內(nèi)部應(yīng)力的影響分析

圖5為線性材料、一次非線性和二次非線性3種情況下的鋁帶卷徑向壓力與切向應(yīng)力分布情況，從這些分布情況可得出，線性材料模型計(jì)算出的徑向壓力值以及切向應(yīng)力值均比非線性模型結(jié)果要大得多，這是由于線性模型雖然考慮了鋁帶卷的各向異性特性，但沒有考慮鋁帶卷層與層之間的緊密程度對內(nèi)部應(yīng)力分布的影響。兩種非線性模型的計(jì)算結(jié)果基本接近，說明對徑向壓縮系數(shù)采用一次非線性擬合就足夠了。

圖4 鋁帶卷應(yīng)力分布Fig.4 Stress distribution of aluminum coil

圖5 材料模型對鋁帶卷壓力的影響Fig.5 Influences of material model on stress of aluminum strip

3.3 卷筒剛度對內(nèi)部應(yīng)力的影響分析

卷筒作為鋁帶卷取的支撐機(jī)構(gòu)，其剛度對鋁帶卷應(yīng)力分布影響很大。現(xiàn)研究卷筒剛度Ec分別取5, 28,61 MPa/mm時，鋁帶卷應(yīng)力分布的不同。從圖6可以看出，卷筒剛度對鋁帶卷內(nèi)層應(yīng)力影響很大，而對外層影響很小，且外層幾乎按相同曲線變化。卷筒剛度較小時，雖然可以降低卷芯處的徑向壓力，但切向應(yīng)力卻變成壓縮應(yīng)力，此時卷芯處容易產(chǎn)生塌卷缺陷，卷筒剛度較大時，內(nèi)部徑向壓力與切向應(yīng)力都比較大，鋁帶材容易出現(xiàn)扁卷現(xiàn)象，所以需要合理選擇卷筒的剛度。

圖6 卷筒剛度對鋁帶卷應(yīng)力的影響Fig.6 Influences of drum stiffness on stress of aluminum strip

4 結(jié)論

1) 冷軋鋁帶卷取內(nèi)層徑向壓力最大，并且隨卷徑的增大而減小，其中內(nèi)部區(qū)域與外部區(qū)域徑向壓力沿半徑方向下降較快，而中部區(qū)域變化較為平緩；帶卷內(nèi)部切向應(yīng)力較大，隨后減小，在外層一定區(qū)域開始再增大，最外層切向應(yīng)力與卷取張應(yīng)力相等。

2) 徑向壓縮系數(shù)的理論值與實(shí)驗(yàn)值相符；計(jì)算內(nèi)部應(yīng)力分布時，需要考慮鋁帶卷的各向異性以及層與層之間的緊密程度。

3) 卷筒剛度對鋁帶卷內(nèi)層應(yīng)力影響較大，而對外層影響很小，卷筒剛度較小時，在卷芯處容易產(chǎn)生塌卷缺陷，卷筒剛度較大時，容易出現(xiàn)扁卷現(xiàn)象，合理選擇卷筒的剛度，有利于鋁帶材的卷取。

[1] 連家創(chuàng), 李春滿. 軋鋼設(shè)備及工藝(第一分冊鋼板軋機(jī))[M]. 哈爾濱: 東北重型機(jī)械學(xué)院出版社, 1985:134—138.LIAN Jia-chuang, LI Chun-man. Rolling Equipment and Technology (The First Volume Plate Mill)[M]. Harbin:Northeast Heavy Machinery Institute Press, 1985: 134—138.

[2] 周國盈. 帶鋼精整設(shè)備[M]. 北京: 機(jī)械工業(yè)出版社,1982.ZHOU Guo-ying. Strip Finishing Equipment[M]. Beijing:Machinery Industry Press, 1982.

[3] PARK W W, KIM D K, IM Y T, et al. Effects of Processing Parameters on Elastic Deformation of the Coil During the Thin-strip Coiling Process[J]. Metals & Materials International, 2014, 20(4): 719—726.

[4] 白振華, 司紅鑫, 周慶田, 等. 熱軋帶鋼卷取過程中鋼卷內(nèi)部應(yīng)力模型的研究[J]. 機(jī)械工程學(xué)報(bào), 2014, 50(2):110—115.BAI Zhen-hua, SI Hong-xin, ZHOU Qing-tian, et al. Research on Steel Roll Internal Stress Model in Hot Rolled Strip Coiling Process[J]. Chinese Journal of Mechanical Engineering, 2014, 50(2): 110—115.

[5] PARK K T, YONG H P, PARK H C, et al. Stress Analysis of Cold Rolled Strip Coiling Process[J]. Transactions of the Korean Society of Mechanical Engineers A, 2017,41(5): 409—414.

[6] SUN W Q, ZHANG X Y, HE A R, et al. Research on Ridge Buckle of High Accuracy Cold Rolled Thin Strip at Local High Points During the Coiling Process[J]. La Metallurgia Italiana, 2016, 108(1): 45—52.

[7] 白振華, 連家創(chuàng), 吳彬, 等. 冷軋鋼卷卷取過程內(nèi)部應(yīng)力三維分布的研究[J]. 鋼鐵, 2001, 36(9): 42—45.BAI Zhen-hua, LIAN Jia-chuang, WU Bin, et al. Research of Three-dimensional Distribution of Internal Stress in Strip Coiling Process[J]. Steel, 2001, 36(9):42—45.

[8] WANG Yong-qin, LI Li, YAN Xing-chun, et al. Modeling of Stress Distribution During Strip Coiling Process[J].Journal of Iron and Steel Research (International), 2012,19(8): 6—11.

[9] HOOG F R D, COZIJNSEN M, YUEN W Y D, et al.Predicting Winding Stresses for Wound Coils of Linear Orthotropic Material[J]. ARCHIVE Proceedings of the Institution of Mechanical Engineers Part C Journal of Mechanical Engineering Science, 2004, 218(1): 27—38.

[10] LIU M L. A Nonlinear Model of Center-wound Rolls Incorporating Refined Boundary Conditions[J]. Computers & Structures, 2009, 87(9/10): 552—563.

[11] 胡小虎. 帶鋼熱連軋地下卷取機(jī)卷筒壓力的有限元解法[D]. 秦皇島: 燕山大學(xué), 2001.HU Xiao-hu. Finite Element Solution of Scroll Pressure of Down-coiler for Hot Strip Mill[D]. Qinhuangdao:Yanshan University, 2001.

[12] 李海亮. 冷軋卷取機(jī)振動及張力分布的研究[D]. 秦皇島: 燕山大學(xué), 2008.LI Hai-liang. The Research on Coiling Vibration and Tension Distribution[D]. Qinhuangdao: Yanshan University, 2008.

[13] 常鐵柱, 張清東, 丁軍. 冷軋帶鋼卷取過程中受力與變形研究[J]. 鋼鐵, 2008, 43(4): 51—56.CHANG Tie-zhu, ZHANG Qing-dong, DING Jun. Study on Deformation of Cold-rolled Steel Strip During Coiling[J]. Steel, 2008, 43(4): 51—56.

[14] 連家創(chuàng). 卷取機(jī)卷筒單位壓力的計(jì)算[J]. 燕山大學(xué)學(xué)報(bào), 1979, 16(Z1): 4—12.LIAN Jia-chuang. Calculation of Unit Pressure in Drum of Coiler[J]. Journal of Yanshan University, 1979, 16(Z1):4—12.

[15] HAKIEL Z. Nonlinear Model for Wound Roll Stress[J].Tappi Journal, 1987, 70.

[16] Stormy Attaway. MATLAB編程與工程應(yīng)用[M]. 北京:電子工業(yè)出版社, 2013.Stormy Attaway. MATLAB Programming and Engineering Applications[M]. Beijing: Electronics Industry Press,2013.

[17] 丁修堃. 軋制過程自動化[M]. 北京: 冶金工業(yè)出版社,2009.DING Xiu-kun. Rolling Process Automation[M]. Beijing:Metallurgical Industry Press, 2009.

Finite Difference Simulation of Internal Stress for Cold Rolled Aluminum Stripcoil

ZHANG Bo1,WANG Wei1,2,LI Shi-miao1
(1. College of Mechanical Engineering and Automation, Fuzhou University, Fuzhou 350100, China;2. Collaborative Innovation Center of High-End Equipment Manufacturing in Fujian, Fuzhou 350100, China)

The paper aims to study distribution of internal stress in coiling of cold rolled aluminum strip and influencing factors to provide reference for design of aluminum strip drum. Through the cold rolling of aluminum compression test, relation curves for radial compression coefficient and stress of aluminum strip were determined; nonlinear finite difference method was used to simulate coiling process and internal stress; effects of linear, first-order nonlinear and second-order nonlinear material models and stiffness of drum on stress distribution inside the roll were analyzed. The calculated radial compressive coefficient was in good agreement with the experimentally values; drum stiffness had great influence on inner stress of aluminum coil, but had little effect on the outer layer. When stiffness of drum was small, it was easy to produce collapse defect at core, when drum stiffness was large, it was easy to flat rolling phenomenon. When calculating internal stress distribution of a coil, it is necessary to consider the anisotropy of the aluminum coil and the tightness between the layers, reasonable choice of rigidity of drum has great practical significance on improving the coil quality.

cold rolled aluminum strip; finite difference method; radial compression coefficient

2017-12-08

國家自然科學(xué)基金（51375094）；福建省科技計(jì)劃（2018H0101）

張波（1993—），男，碩士研究生，主要研究方向?yàn)榻饘俨牧暇艹尚喂に嚰把b備設(shè)計(jì)。

王偉（1970—），男，博士，教授，主要研究方向?yàn)楝F(xiàn)代機(jī)械設(shè)計(jì)理論及方法、冶金設(shè)備力學(xué)及控制。

10.3969/j.issn.1674-6457.2018.01.016

TG335

A

1674-6457(2018)01-0127-05