一种预测带钢热连轧过程温度演变的方法

摘要

本发明公开了一种预测带钢热连轧过程温度演变的方法，包括以下步骤：建立带钢热传导偏微分方程；建立有限元模型；输入带钢热连轧工艺参数；计算带钢在除鳞区的温度；计算带钢在粗轧区的温度；计算带钢在精轧区的温度；计算带钢在层流冷却区的温度；输出带钢热连轧温度的计算结果。本发明通过建立有限元模型、求解偏微分方程来解决平面二维热传导问题，在此过程中，结合了设备及带材的各种工艺参数（设备的水量、带材速度、带材厚度、带材宽度），考虑了热连轧全程各个环节的温度演变，从而对包括带钢从出加热炉到高压水除鳞、粗轧、精轧、层流冷却全程的温降进行预测。有了这样详细的温度预测数据结果，能够为热连轧工艺制度的优化提供参考。

主权项

1.一种预测带钢热连轧过程温度演变的方法，其特征在于：包括以下步骤：A、建立带钢热传导偏微分方程假设带钢轧件沿长度方向的温度分布相对均匀，只考虑带钢轧件厚度和宽度方向的二维热传导；对于含内热源的二维问题，固体热传导微分方程的形式如下式所示：$\frac{\partial T}{\partial t}=\frac{k}{{\mathrm{\rho c}}_p}(\frac{{\partial }^{2}T}{{\partial x}^2}+\frac{{\partial }^{2}T}{{\partial y}^2}+\frac{\stackrel{·}{q}}{k})---\left(1\right)$式中,T为带钢的瞬态温度；t为热传导过程进行的时间；k为材料的导热系数；ρ为材料的密度；c_p为材料的定压比热；表示内热源强度，为每单位体积在单位时间的发热量；x和y分别为带钢轧件宽度和厚度方向的直角坐标；B、建立有限元模型假设沿带钢长度方向上温度分布均匀，按二维传热问题处理；利用带钢轧件的对称性，假设冷却方式对称，取四分之一横断面进行网格单元划分，采用四节点等参单元，在所选取的四分之一横断面上划分单元网格50×9个；求解平面二维热传导问题是偏微分方程（1）在给定的边界条件和初始条件之下求解的问题，它等效地表达为以下泛函形式：$J\left[T(x,y)\right]=\int \underset{D}{\int }\left\{\frac{k}{2}\right[{\left(\frac{\partial T}{\partial x}\right)}^2+{\left(\frac{\partial T}{\partial y}\right)}^2]-{\mathrm{\rho c}}_p\frac{\partial T}{\partial t}T\}\mathrm{dxdy}+\frac{1}{2}{\int }_{\tau }h{(T-T_{\infty })}^2\mathrm{ds}---\left(2\right)$能使泛函式（2）的一阶变分为零的温度函数，必为此问题的正确解；C、输入带钢热连轧工艺参数输入的工艺参数包括钢种、化学成分、热物性参数、坯料尺寸、出炉温度、设备参数、粗轧和精轧轧制规程和层流冷却制度，所述的层流冷却制度包括前段冷却规程、后段冷却规程和分段冷却规程；D、计算带钢在除鳞区的温度除鳞区指从加热炉出口至轧前除鳞结束时带钢轧件经历的区域，计算带钢轧件在经历这一区域过程中的温度场；首先进行除鳞区设备参数的输入，设备参数包括除鳞前空冷距离和除鳞设备长度，然后进行出炉后的空冷温度场模拟，再进行除鳞区的水冷温度场模拟；带钢轧件在辊道上运行过程中，通过自由表面和外界进行热交换，主要换热方式有辐射和对流，带钢轧制时的对流传热为自然对流，在带钢轧件温度>800℃时它引起的温降仅为辐射温降的1%，因此，在高温区段空冷过程中，忽略热对流的影响；利用Stefan-Boltzmann定律，辐射换热的热交换系数如下：$h_r=\mathrm{ϵS}(T_{\mathrm{su}}^2+T_{\infty }^2)(T_{\mathrm{su}}+T_{\infty })---\left(3\right)$式中，S为Stefan-Boltzmann常数；ε为辐射率或黑度系数；T_su为带钢表面温度；T_∞为环境温度；在对称面上无热交换发生，为绝热边界条件，其热交换量为Q＝0；对于高压水除鳞过程，带钢表面温度T_su、水压以及水量密度ω对其热交换系数h_w（W/（m²·K））的影响较大，尤其是高压水除鳞的水压高达20MPa，对热交换系数的影响很大；随着表面温度不同，各温度区间的传热现象也不同；对于这种水冷形式的计算公式:$h_w=r\times 107.2{\omega }^{0.663}\times {10}^{{-0.00147T}_{\mathrm{SU}}\times 1.163}---\left(4\right)$式中：r—水压影响系数，高压水除鳞时r>1；E、计算带钢在粗轧区的温度粗轧区包括带钢轧件从除鳞区出口到粗轧末道次出口中间经历的所有区域，该区域中带钢轧件经历了空冷、轧制前除鳞和轧制变形过程；其中空冷与轧前除鳞计算模型与公式（3）和（4）相同；轧制过程变形区中主要存在以下几种热现象：带钢轧件塑性变形转化成的变形热、带钢轧件与轧辊相对滑动产生的摩擦热、带钢轧件与轧辊接触产生的接触热传导；E1、变形热：把塑性功转变的热能看作是一种内热源，其变形热的表达式为：$W_e=\eta \int \int {\int }_V\bar{\sigma }\stackrel{\stackrel{·}{‾}}{ϵ}\mathrm{dV}---\left(5\right)$式中，－平均变形速率；－等效应力；η－塑性功转变为热的比例，取0.7~0.9；E2、摩擦热：在轧制过程中因轧辊与带钢轧件间的摩擦而产生的单位表面积、单位时间的摩擦功由下式表示：${\stackrel{·}{W}}_f=\mathrm{\mu p}\left(\phi \right)v_r---\left(6\right)$式中，μ、p(φ)和v_r分别为摩擦系数、轧制压力的法向分量和轧辊与带钢轧件的相对速度的绝对值；E3、接触热传导：在带钢轧件和轧辊接触面上，热量由带钢轧件传向轧辊，此为热流连续、温度不连续的热阻问题，热交换系数h按下式处理：$h_r=k/\sqrt{\mathrm{\pi at}}---\left(7\right)$$t=\sqrt{\mathrm{\Delta H}·R}/v---\left(8\right)$a=k/(ρ·c)   （9）式中，t为带钢轧件与轧辊接触时间，v为轧制速度，ΔH为压下量；R为轧辊半径；F、计算带钢在精轧区的温度精轧区指带钢轧件从粗轧末道次轧辊出口到精轧末道次机架轧辊出口经历的区域，中间经历空冷、保温、精轧前除鳞和机架间喷水冷却与轧制变形过程；其中机架间喷水冷却的热交换系数计算公式：$h_w=107.2{\omega }^{0.663}\times {10}^{{-0.00147T}_{\mathrm{SU}}\times 1.163}---\left(10\right)$G、计算带钢在层流冷却区的温度对于层流冷却过程，是指冷却水在无压力的情况下流向带钢，形成冲击水区和稳定水膜沸腾区，该过程是相对比较复杂的传热过程，对于此过程的热交换系数主要受设备条件、冷却水量和带钢表面温度的影响；对于此过程交换系数采用以下计算公式:$h_w=\frac{9.72\times {10}^5{\omega }^{0.355}}{(T_{\mathrm{SU}}-T_W)}{\left\{\frac{(2.5-1.5\log T_W)D}{P_L{P}_C}\right\}}^{0.645}\times 1.163---\left(11\right)$式中：ω—水量密度，m³/（min·m²）;D—喷嘴直径，m;T_SU，T_W—带钢表面温度和水温，℃；P_L，P_C—轧线方向和轧线垂直方向的喷嘴间距，m；H、输出带钢热连轧温度的计算结果输出带钢从出加热炉到高压水除鳞、粗轧、精轧、层流冷却过程中温度预测结果，输出带钢由出炉后直到层冷结束全程的温降曲线，所述的降温曲线包括心部温度、表面温度、平均温度；输出带钢头部、中部及尾部由出炉后直到层冷结束全程的典型横断面温度分布，以及断面上代表点的温度值；所述的典型横断面为1/4横断面。