方矩形管连续辊弯成型延伸规律
付志强1,2, 安子军1, 杜凤山1,2, 段利英1, 于辉1, 杨帆1,2
1.燕山大学 机械工程学院,河北 秦皇岛,066004
2.燕山大学 国家冷轧板带装备及工艺工程技术研究中心,河北 秦皇岛,066004
杜凤山(1960-),男,教授,博士生导师.研究方向:轧制设备、工艺及其数字化技术.E-mail:fsdu@ysu.edu.cn

付志强(1982-),男,博士研究生.研究方向:型钢孔型设计及过程仿真技术.E-mail:fzq6670300@163.com

摘要

为了减小方矩形管辊弯成型时纵向延伸,结合某厂实际生产工艺,建立了方矩形管辊弯成型七机架的三维弹塑性非线性有限元模型。通过仿真分析得到了方矩形管辊弯成型时横向金属流动规律及纵向延伸率的分布规律,与现场测量结果一致。分析了轧辊转速、管坯壁厚、辊径位置等因素对延伸率的影响,在此基础上优化了辊弯轧制工艺参数。结果表明:优化后的辊弯轧制工艺提高了方矩形管的成型质量,为预测最终成型尺寸精度和辊弯成形工艺参数的制定提供了理论依据。

关键词: 连续轧制; 方矩形管; 辊弯成型; 延伸率; 有限元
中图分类号:TG335 文献标志码:A 文章编号:1671-5497(2015)02-0487-07
Elongation of rectangular tube in continuous roll forming
FU Zhi-qiang1,2, AN Zi-jun1, DU Feng-shan1,2, DUAN Li-ying1, YU Hui1, YANG Fan1,2
1.School of Mechanical Engineering, Yanshan University, Qinhuangdao 066004, China
2. National Engieering Research Center for Equipment and Technology of Cold Strip Rolling, Yanshan University, Qinhuangdao 066004, China
Abstract

In order to reduce the longitudinal extension of rectangular tube in the roll forming process, a three-dimensional elastic-plastic nonlinear Finite Element Model (FEM) of the roll forming is established with the actual factory production process. The metal flow law in the roll forming of rectangular tube and the distribution of the longitudinal elongation rate are obtained by analyzing the simulation results. The obtained results are consistent with the field measurements. The effects of the roller speed, the thickness of the tube wall and the location of the roll diameter on the elongation are analyzed using the proposed FEM. Then these processing parameters are optimized, which improve the continuous rolling performance of rectangular tubes. This research may help to predict the final dimensional accuracy and determine the processing parameters of the roll forming.

Keyword: continuous rolling; rectangular tube; roll forming; elongation; finite element method (FEM)

引言

近年来, 虚拟仿真FEM技术在轧制领域得到广泛应用[1, 2, 3], 在辊弯成型方面[4, 5, 6], 国内外许多学者做了相关研究, 其中Huang[7]和Leu[8]对方形管四直辊挤压成型过程中的R/t(管坯半径与厚度的比值)、加工硬化和摩擦因数对成型的影响做了比较详细的研究。张小平等[9]对四直辊辊弯成型方形管过程进行了数值模拟和实验研究。徐树成、刘才[10]对方矩形管四曲辊辊弯成型过程中前两道次的变形和应力进行了分析。杜凤山等[11]采用连续加载方式对方管辊弯成型过程中的变形和流场做了分析。大多数研究是关于四直辊的单机架和少数机架辊弯成型时的轧制力、横截面的变形、应力应变分布等, 忽略了对连续辊弯成型时金属的流动规律及延伸规律的研究。本文采用FEM技术对四曲辊连续辊弯成型方矩形管时的金属流动规律及不同工艺参数对延伸系数分布规律的影响做了详细研究, 并与现场实测数据进行了对比分析, 仿真结果与现场测得的结果相一致。

1 有限元模型的建立
1.1 辊弯生产现场

辊弯生产方矩形管是通过顺序配置的多道次具有特定轮廓型面的成型轧辊, 以圆管作坯料, 经过几架集中传动的二辊或四辊机架逐渐弯曲变形而成方矩形管。本文研究四辊辊弯生产方矩形管的辊弯过程, 图1为某厂四辊轧制机组的现场轧制照片。

图1 现场轧制照片Fig.1 Rolling scene photo

1.2 辊弯轧辊参数及几何模型

表1为四辊辊弯180 mm× 100 mm× 8 mm方矩形管轧制工艺及轧辊参数。根据现场的工艺参数, 可以得到轧辊轮廓曲线的变化过程, 如图2(a)所示, 从图中可以看出:第一道次起咬入作用, 后两道次起精成型作用。四辊轧机的轧辊几何简图如图2(b)所示, 可以看出辊缝开在角部成形处。

表1 轧制工艺参数 Table 1 Process parameters of roll forming
1.3 有限元仿真模型

根据现场轧制工艺得到的轧辊参数及方矩形管轧制的特点, 对仿真模型进行1/4对称简化, 对称边界条件由对称面上的节点速度来确定, 在对称面法线方向上的速度为零。假定母管横向和纵向壁厚均匀, 轧件视作变形体, 轧辊视作刚性体。轧辊分为平辊和立辊, 平辊为驱动辊, 立辊为从动辊, 从动辊需要建立两个辅助节点, 一个为主控辅助节点, 一个为旋转控制辅助节点。旋转控制辅助节点控制为从动, 初始辊速为0, 靠管坯带动其转动与实际情况相符。管坯后加后推板, 后推板给管坯一个力实现管坯的顺利咬入。咬入后, 让后推板与管坯分离, 实现管坯的稳定轧制, 有限元模型如图3所示。采用八节点单元建立描述管坯的有限元网格, 网格划分时沿钢管半径方向分为两层, 圆周方向分为40段, 共123个节点, 管坯长1000 mm, 分为200段。管坯材料选用Marc材料库中的C22, 它的材料特性与我国20钢材料特性一致。管坯和轧辊之间采用库伦摩擦模型, 摩擦因数取0.3。

图2 几何模型Fig.2 Geometric model

图3 有限元仿真模型Fig.3 FEM model

2 有限元模型的实验验证
2.1 分析方法

在管坯纵向上离端头50 mm位置取点A, 管坯中点取点B, 离尾端50 mm位置取点C, 稳定段A-C, 如图4(a)所示。在管坯的横截面上外层的节点位置分布如图4(b)所示。图中1号位置到20号位置为驱动辊辊弯的长边, 20号位置到30号位置为辊缝处的角部, 30号位置到40号位置为从动辊辊弯的短边。从动辊辊弯方向为 方向, 驱动辊辊弯方向为 方向, 辊弯轧制方向为 方向。

图4 节点选择分布图Fig.4 Positions of points selected

2.2 延伸系数

管坯在 向的延伸系数可以反映矩形管成型过程中金属在 向的流动情况。在本模型中设稳定段 (见图4)的延伸系数为 如下式所示:

式中: 分别为 段辊弯轧制前和辊弯轧制后的长度。

2.3 仿真结果与实测结果对比

轧制前在坯料圆管的中间稳定段绘制400 mm长的等距线, 如图5(a)所示, 管坯经过辊弯轧制后 向上有一定的延伸, 测量其轧后长度如图5(b)所示。分别对图4(b)中1、10、20、30、40号节点附近位置的延伸量进行了大量的实测, 通过测得的平均数据和式(1)可以得到实际辊弯轧制后管坯 向上的延伸系数分布情况。将实测值与仿真值进行对比, 如图6所示。从图中可以看出:仿真计算得到的延伸系数分布规律与实际测得的延伸系数分布规律基本一致。

图5 现场测量照片Fig.5 Field measurements photos

图6 仿真计算和现场测量对比Fig.6 Simulation calculation and contrast

3 延伸不均匀性的成因仿真分析
3.1 辊弯成型时X方向位移分布规律

图7 B点横截面上节点的X向位移Fig.7 X-direction displacement of the node in the cross section of B

图7为辊弯仿真计算结果稳定段 点处横截面外层各节点从第一架到第七架 向位移分布图。从图中可以看出:在驱动辊弯面上的0到10号节点位置附近位移量为负。这是由于在辊弯过程中, 长边处的金属会向角部流动, 角部外层由于处于拉伸状态, 金属会向长边处流动, 引起金属对流导致的结果。

3.2 辊弯成型时Y方向位移分布规律

图8为辊弯仿真计算结果稳定段 点处横截面外层各节点从第一架到第七架 方向位移分布图。从图中可以看出:从25号节点到40号节点位移量为负, 在整个辊弯过程中, 从动辊辊弯曲面变形时金属向角部的流动能力要小于驱动辊, 导致纵向上延伸的不均匀。

图8 B点横截面上节点的Y向位移Fig.8 Y-direction displacement of the node in the cross section of B

3.3 辊弯成型时Z方向位移分布规律

图9为辊弯仿真计算结果稳定段 点横截面外层各节点的 方向位移量对比图。从图中可以看出:在辊弯轧制方矩形管时, 纵向前段的位移比后边的位移大, 说明在整个辊弯轧制过程中纵向上有一定的延伸, 且延伸量是不一致的, 该结论与现场轧制实测情况一致。

图9 ABC点处横截面上节点的Z向位移Fig.9 Z-direction displacement of the node in the cross section of ABC

4 延伸规律的仿真计算结果分析
4.1 辊弯成型时七机架间延伸系数分布规律

图10为辊弯轧制成型时七机架间稳定段延伸系数的分布规律。从图中可以看出:随着压下量的增大, 纵向上延伸系数逐渐增大, 在驱动辊辊弯面上延伸系数最小, 从动辊辊弯面上延伸系数最大, 角部次之。

图10 延伸系数分布规律Fig.10 Distribution of longitudinal extension coefficient

4.2 辊速对轧后纵向延伸系数的影响

只改变轧辊的转速, 保持其他轧制工艺参数不变, 其延伸系数变化规律分布如图11所示。从图中可以看出:随着轧辊转速的变大, 驱动辊辊弯面上的延伸系数逐渐变大, 而从动辊辊弯面上的延伸系数在轧辊的转速高于5 r/min后变化范围较小。

图11 不同轧速下纵向延伸系数分布图Fig.11 Distribution of longitudinal extension coefficient at different rolling speed

4.3 辊径位置对轧后纵向延伸系数的影响

保持其他工艺参数不变, 改变机架前后轧辊的直径建立位置如图12所示。仿真计算后的延伸系数分布规律如图13所示。从图中可以看出:不同部位建立轧辊辊径对延伸系数有一定的影响, 辊底建立辊径工艺下的延伸系数变化比较平缓。

图12 在不同位置建立辊径的示意图Fig.12 Schematic diagram of roll diameter established at different locations

图13 在不同位置建立辊径的纵向延伸系数分布图Fig.13 Distribution of longitudinal extension coefficient at different locations to establish roll diameter

4.4 壁厚对轧后纵向延伸系数的影响

只改变管坯的壁厚, 保持其他轧制工艺参数不变, 纵向延伸系数的分布规律如图14所示。从图中可以看出:随着壁厚的增加, 辊弯轧制后管坯纵向延伸系数逐渐增大。

图14 不同壁厚下纵向延伸系数分布图Fig.14 Distribution of longitudinal extension coefficientwith different wall thickness

5 工艺优化

纵向上延伸不均匀会导致轧后纵向上残余应力分布的不均匀, 是引起开裂的原因之一。通过仿真计算得到延伸规律, 对方矩形管工艺进行优化, 优化结果如表2所示。

表2 优化后的工艺参数 Table 2 Optimized process parameters

对照表1表2可以看出:经过优化改变了从动辊辊弯的曲面, 降低了驱动辊的转速。按照辊底建立前、后轧辊的直径差, 优化前、后延伸系数分布以及B点处的纵向残余应力对比如图15所示。从图15(a)可以看出:优化后纵向上延伸系数减小, 从驱动辊辊弯曲面上的延伸量到从动辊辊弯曲面上的延伸量的差值也减小。从图15(b)可以看出:优化后消除了纵向上16号节点到22号节点间的下峰值, 降低了16号节点到22号节点和32号节点到36号节点之间的上峰值, 角部的峰值改变不明显, 仍需要进一步研究。

图15 优化前后对比图Fig.15 Comparison chart before and after optimization Simulation on residual stress

6 结论

(1)根据现场的辊弯轧制工艺参数, 建立了与实际相符的有限元仿真模型, 根据仿真结果得到了连续辊弯成型后 向上延伸率的分布规律, 与现场实测得到的分布规律比较一致, 说明建立的有限元仿真模型具有一定的实用性。

(2)通过对方矩形管辊弯成型的仿真模拟, 分析仿真过程中的 方向的位移分布, 可知辊弯成型方矩形管时两个直边向角部流动能力的不同是导致纵向上延伸量不一致的主要原因。

(3)分析了壁厚、辊径建立位置、轧制速度对方矩形管辊弯成型中延伸率的影响, 得到了壁厚对连续辊弯成型时纵向延伸率的影响最大, 确定了建立辊径的理想位置及辊弯成型时的轧制速度的选择方案, 为方矩形管辊弯成型工艺参数的制定提供了理论依据。

(4)经过工艺优化, 消除了长边到角部之间残余应力的下峰值, 降低了长边、短边和角部之间的上峰值和开裂几率, 角部峰值改变不明显, 仍需要进一步的研究。

The authors have declared that no competing interests exist.

参考文献
[1] 赵玲玲, 安子军, 杜凤山, . 热变形对42CrMo钢奥氏体晶粒尺寸影响的试验研究[J]. 塑性工程学报, 2010, 17(4): 100-103.
Zhao Ling-ling, An Zi-jun, Du Feng-shan, et al. Experimental research on the effect of hot deformation on the growth of austenite grain size of 42CrMo steel[J]. Journal of Plasticity Engineering, 2010, 17(4): 100-103. [本文引用:1] [CJCR: 0.526]
[2] 董永刚, 张文志, 宋剑锋. 钢轨万能轧制过程金属延伸规律的理论与试验研究[J]. 机械工程学报, 2010, 46(6): 87-92.
Dong Yong-gang, Zhang Wen-zhi, Song Jian-feng. Theoretical and experimental research on the elongation law of the rail in rail rolling by a universal mill[J]. Journal of Mechanical Engineering, 2010, 46(6): 87-92. [本文引用:1] [CJCR: 0.544]
[3] 汪飞雪, 杜凤山, 于辉, . 基于有限元法的PQF连轧过程宽展规律的仿真[J]. 钢铁, 2013, 48(3): 51-55.
Wang Fei-xue, Du Feng-shan, Yu Hui, et al. Spread rule during PQF seamless tube rolling process based on FEM[J]. Iron and Steel, 2013, 48(3): 51-55. [本文引用:1] [CJCR: 0.52]
[4] Rossi B, Degée H, Boman R. Numerical simulation of the roll forming of thin-walled sections and evaluation of corner strength enhancement[J]. Finite Elements in Analysis and Design, 2013, 72: 13-20. [本文引用:1] [JCR: 1.389]
[5] 韩飞, 刘继英, 艾正青, . 辊弯成型技术理论及应用研究现状[J]. 塑性工程学报, 2010, 17(5): 53-60.
Han Fei, Liu Ji-ying, Ai Zheng-qing, et al. State of the art of research on roll forming process[J]. Journal of Plasticity Engineering, 2010, 17(5): 53-60. [本文引用:1] [CJCR: 0.526]
[6] Tehrani M S, Hartley P, Naeini H M, et al. Localised edge buckling in cold roll-forming of symmetric channel section[J]. Thin-Walled Structures, 2006, 44(2): 184-196. [本文引用:1] [JCR: 1.231]
[7] Huang Y M. Elasto-plastic finite element analysis of squaring circular tube[J]. Transactions of Nonferrous Metals Society of China, 2008, 18(3): 665-673. [本文引用:1] [JCR: 0.917] [CJCR: 0.7776]
[8] Leu D K. The shaping of a circular tube into a symmetric square-tube by finite-element modeling[J]. Journal of Materials Processing Technology, 2006, 178(1): 287-296. [本文引用:1] [JCR: 1.953]
[9] 张小平, 李志永, 张少琴, . 异型管辊挤成型过程的数值模拟与实验研究[J]. 塑性工程学报, 2005, 12(4): 38-42.
Zhang Xiao-ping, Li Zhi-yong, Zhang Shao-qin, et al. FEM analysis and experiment study of non mand rel double pass roll extrusion process of shaped tube[J]. Journal of Plasticity Engineering, 2005, 12(4): 38-42. [本文引用:1] [CJCR: 0.526]
[10] 徐树成, 刘才. 方矩形管连续辊弯成形过程的变形和应力场[J]. 机械工程学报, 2000, 36(3): 93-96.
Xu Shu-cheng, Liu Cai. Deformation and stress field of square rectangle tube in continue roll forming process[J]. Journal of Mechanical Engineering, 2000, 36(3): 93-96. [本文引用:1] [CJCR: 0.544]
[11] 杜凤山, 陆卫平, 张建宇, . 方管闭式辊弯成型过程的变形及流场[J]. 东北重型机械学院学报, 1997, 21(2): 100-103.
Du Feng-shan, Lu Wei-ping, Zhang Jian-yu, et al. Deformation and flow fields of the square tube forming of 4-roll bending[J]. Journal of Northeast Heavy Machinery Institute, 1997, 21(2): 100-103. [本文引用:1]