2. 徐州工程机械集团有限公司 徐工基础工程机械有限公司, 江苏 徐州 221004
2. XCMG Foundation Construction Machinery Co., Ltd., Xugong Construction Machinery Group, Xuzhou 221004, China
难变形金属薄带一般无法采用常规冷轧方法进行最终成品轧制, 但随着轧制温度的提高, 部分难变形金属的变形性能得到明显改善, 因此可以采用温轧方法实现其轧制生产[1-3].东北大学轧制技术及连轧自动化国家重点实验室在前期大量研究工作的基础上, 设计了两侧配备热卷取箱的温轧装备, 实现采用温轧生产难变形金属薄带的工艺路径.
在温轧过程中, 轧材提温和轧制温度精确控制, 是该工艺路线能否实现的关键.本文重点针对金属薄带卷取加热过程温度计算开展研究, 为整个温轧过程的温度精确控制奠定基础.
1 工艺设备概况新型温轧机组的工艺设备概况如图 1所示, 在温轧机的两侧卷取机外部各配备一台热卷取箱.热卷取箱由保温材料、耐火材料和电加热器组成, 利用电加热器通过辐射加热可以将热卷取箱内温度控制在一定范围, 从而保证卷取机中带材的温度.热卷取箱的主要工艺参数如表 1所示.
温轧过程中, 薄带在热卷取箱中的卷取机卷筒上进行卷取保温, 卷取后的带卷近似为圆柱体.带卷温度计算可以看作圆柱体温度计算, 并认为沿圆柱体的径向和轴向传热, 周向温度保持一致.由于卷取后带卷轴向左右对称, 且周向温度一致, 因此, 以轴向长度中心为对称面, 并沿周向4均分后取其中一份, 即取整个圆柱体的1/8作为温度计算的对象, 如图 2所示, 采用二维圆柱坐标系建模.带卷在热卷取箱中温度计算的基础是二维圆柱坐标系导热微分方程, 如式(1)所示:
(1) |
式中:T为待求温度; τ为时间; r, z分别为圆柱坐标系中径向和轴向坐标; ρ为材料密度; c为材料比热容; λr, λz分别为径向和轴向导热系数.其中轴向导热系数λz为薄带材料本身导热系数, 而带卷在径向为多层结构, 径向导热系数λr需要重点研究.
2.2 带卷径向传热带卷温度计算时, 一般认为带卷沿径向各层之间的结构示意图如图 3所示[4-9], 薄带表面存在氧化层, 且两层薄带之间为有限点接触, 存在间隙层, 间隙层中充满气体介质.一个径向单元层结构上下对称, 各部分厚度分别为Hs/2, Ho, Hi, 其中, 上下层薄带厚度各取其厚度的一半.则其总厚度H由式(2)确定.
(2) |
式中:Hs为单层薄带厚度; Ho为氧化层厚度; Hi为间隙层厚度.
径向导热系数计算考虑薄带层、氧化层和间隙层的各部分热阻Rs, Ro和Ri, 间隙层热阻Ri认为是带材接触传导热阻Ri, s、空气传导热阻Ri, a及辐射热阻Ri, rd三部分热阻并联, 如图 4所示.径向导热系数计算采用文献[4]中公式, 如式(3)所示.
(3) |
式中: λs, λo, λa分别为带钢、氧化层和空气的导热系数; ε为钢卷黑度; σ为玻尔兹曼常数; m, n为与表面粗糙度相关的系数; A为接触比例系数, 该参数与钢卷的材质及卷取时的层间压应力有关, 按照式(4)计算.
(4) |
其中:P0为带材的微观硬度;P为层间压应力,层间压应力P需根据初始卷取张力及加热过程卷径变化造成的压力变化确定.
2.3 带卷温度求解为求解式(1), 采用有限差分方法进行数值求解.分别建立轴向显式和隐式差分方程, 以及径向显式和隐式差分方程.并采用交替显式和隐式方向迭代法(ADI)进行求解, 将时间步长一分为二, 在前半个时间步长内, z方向(轴向)采用显式格式, r方向(径向)采用隐式格式; 在后半个时间步长内, z方向采用隐式格式, r方向采用显式格式[10].
带卷加热过程的边界条件, 即热卷取箱内部环境与钢卷表面的传热系数(h)是钢卷加热速度的关键影响因素.钢卷外径受到碳棒辐射、炉膛辐射以及炉内气体对流换热等的综合影响.将边界导热系数简化为辐射换热形式, 并引入修正系数考虑其他因素的影响, 见式(5).在后续试验中, 利用实测温度数据拟合获得修正系数.
(5) |
式中:α为修正系数; θs为钢卷温度; θE为热卷取箱内环境温度; ε为钢卷黑度; σ为玻尔兹曼常数.
3 薄带卷取加热测温实验 3.1 实验过程为了进行测温实验, 需要在常温下进行薄带卷取和焊接热电偶, 直接使用难变形材料无法完成该操作, 因此采用普碳钢作为替代材料.
将实验薄带卷取在热卷取箱卷轴上, 在各层之间焊接K型热电偶线, 并通过专用输入模块接入PLC中, 实现温度的快速测量.带钢选用Q195材质普碳钢, 厚度1.7 mm, 宽度120 mm, 长度8 770 mm.将热卷取箱加热提温, 期间卷轴带着钢卷缓慢旋转, 保证加热均匀.在热卷取箱壁上也安装有热电偶用于检测热卷取箱环境温度.实验照片如图 5所示, 实验过程中热卷取箱门关闭.
炉内环境温度实测点标记为θE, 带卷内部由里向外实测温度点标识θ1, θ2, θ3.分析温度实测值, 发现各测温点温度实测值θE, θ1, θ2, θ3与时间t近似呈线性规律变化, 对各实测温度进行一次线性拟合, 获得热卷取箱内部环境温度以及带卷各实测温度点温度随时间变化规律如式(6)~式(9)所示.各拟合直线的线性相关系数值分别为0.997 99, 0.999 25, 0.998 70和0.993 85, 可以看出, 实测温度的线性拟合精度较高, 见图 6,拟合函数可代替实测温度进行后续计算分析.
(6) |
(7) |
(8) |
(9) |
根据实测温度对式(5)中的修正系数α进行优化.修正系数α从1.0~2.5每隔0.1选取一个值, 采用ADI方法[10]分别计算各点温度随时间变化值, 并与实测温度值比较, 获得最佳的修正系数α为2.0.图 7为θ2, θ3计算温度和实测温度对比图, 温度偏差在±2 ℃以内.
1) 设计开发两侧配备热卷取箱的温轧机组, 用于难加工金属薄带成卷温轧生产.
2) 建立薄带热卷取加热过程的温度计算模型.采用二维圆柱坐标系, 考虑沿带卷的径向和轴向的传热; 将带卷沿径向分为薄带层、氧化层和间隙层, 在计算径向导热系数时考虑这三部分结构; 利用圆柱坐标系二维有限差分法, 分别建立轴向和径向的显式差分方程和隐式差分方程, 采用交替显式和隐式方向迭代法进行数值求解; 边界传热系数以辐射换热为主, 并引入修正系数以考虑其他因素的影响.
3) 进行薄带卷取加热测温实验研究.回归实测温度与加热时间的一次线性函数; 利用实测温度对传热修正系数进行优化, 带卷计算温度和实测温度偏差可以控制在±2 ℃以内.
[1] |
Ye F, Liang Y F, Wang Y L.
Fe-6.5wt%Si high silicon steel sheets produced by cold rolling[J]. Materials Science Forum, 2009, 638/639/640/641/642: 1428–1433.
|
[2] |
Fu H D, Zhang Z H, Pan H J, et al.
Warm/cold rolling processes for producing Fe-6.5wt% Si electrical steel with columnar grains[J]. International Journal of Minerals, Materiallurgy and Materials , 2013, 20(6):535-540.
|
[3] |
矫志杰, 孙涛, 李建平.
难变形材料轧制实验机开发及实验研究[J]. 东北大学学报(自然科学版), 2017, 38(2): 229–233.
( Jiao Zhi-jie, Sun Tao, Li Jian-ping. Pilot mill development for the difficult-to-deform material rolling and experimental research[J]. Journal of Northeastern University(Natural Science), 2017, 38(2): 229–233. DOI:10.3969/j.issn.1005-3026.2017.02.016 ) |
[4] |
Park S J, Hong B H, Baik S C.
Finite element analysis of hot rolled coil cooling[J]. ISIJ International, 1998, 38(11): 1262–1269.
DOI:10.2355/isijinternational.38.1262 |
[5] |
Saboonchi A, Hassanpour S.
Heat transfer analysis of hot-rolled coils in multi-stack storing[J]. Journal of Materials Processing Technology, 2007, 182(1/2/3): 101–106.
|
[6] |
Saboonchi A, Hassanpour S.
Simulation-based prediction of hot-rolled coil forced cooling[J]. Applied Thermal Engineering, 2008, 28(13): 1630–1637.
DOI:10.1016/j.applthermaleng.2007.10.029 |
[7] |
Cho H H, Cho Y G, Im Y R.
A finite element analysis for asymmetric contraction after coiling of hot-rolled steel[J]. Journal of Materials Processing Technology, 2010, 210(6/7): 907–913.
|
[8] |
Zhang X, Yu F, Wu W.
Application of radial effective thermal conductivity for heat transfer model of steel coils in HPH furnace[J]. International Journal of Thermophysics, 2003, 24(5): 1395–1405.
DOI:10.1023/A:1026115521232 |
[9] |
祁卫东.罩式退火炉传热模型与带钢氧化控制研究[D].沈阳: 东北大学, 2013.
( Qi Wei-dong.Research on heat transfer model of bell-type annealing furnace and oxidization control of strip steel[D].Shenyang: Northeastern University, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10145-1016011415.htm ) |
[10] |
杨世铭, 陶文铨.
传热学[M]. 北京: 高等教育出版社, 2006: 43-46.
( Yang Shi-ming, Tao Wen-quan. Heat transfer[M]. Beijing: Higher Education Press, 2006: 43-46. ) |