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

(19)中华人民共和国国家知识产权局(12)发明专利(10)授权公告号(10)授权公告号CNCN103028615103028615B(45)授权公告日2014.12.10(21)申请号201210501286.2(22)申请日2012.11.29(73)专利权人一重集团大连设计研究院有限公司地址116600辽宁省大连市经济技术开发区东北大街96号专利权人中国第一重型机械股份公司(72)发明人任德方赵华国刘付强(74)专利代理机构大连东方专利代理有限责任公司21212代理人李洪福(51)Int.Cl.B21B37/74(2006.01)审查员陈成权权利要求书3页利要求书3页说明书6页说明书6页附图5页附图5页(54)发明名称一种预测带钢热连轧过程温度演变的方法(57)摘要本发明公开了一种预测带钢热连轧过程温度演变的方法，包括以下步骤：建立带钢热传导偏微分方程；建立有限元模型；输入带钢热连轧工艺参数；计算带钢在除鳞区的温度；计算带钢在粗轧区的温度；计算带钢在精轧区的温度；计算带钢在层流冷却区的温度；输出带钢热连轧温度的计算结果。本发明通过建立有限元模型、求解偏微分方程来解决平面二维热传导问题，在此过程中，结合了设备及带材的各种工艺参数（设备的水量、带材速度、带材厚度、带材宽度），考虑了热连轧全程各个环节的温度演变，从而对包括带钢从出加热炉到高压水除鳞、粗轧、精轧、层流冷却全程的温降进行预测。有了这样详细的温度预测数据结果，能够为热连轧工艺制度的优化提供参考。CN103028615BCN1032865BCN103028615B权利要求书1/3页1.一种预测带钢热连轧过程温度演变的方法，其特征在于：包括以下步骤：A、建立带钢热传导偏微分方程假设带钢轧件沿长度方向的温度分布相对均匀，只考虑带钢轧件厚度和宽度方向的二维热传导；对于含内热源的二维问题，固体热传导微分方程的形式如下式所示：式中,T为带钢的瞬态温度；t为热传导过程进行的时间；k为材料的导热系数；ρ为材料的密度；cp为材料的定压比热；表示内热源强度，为每单位体积在单位时间的发热量；x和y分别为带钢轧件宽度和厚度方向的直角坐标；B、建立有限元模型假设沿带钢长度方向上温度分布均匀，按二维传热问题处理；利用带钢轧件的对称性，假设冷却方式对称，取四分之一横断面进行网格单元划分，采用四节点等参单元，在所选取的四分之一横断面上划分单元网格50×9个；求解平面二维热传导问题是偏微分方程(1)在给定的边界条件和初始条件之下求解的问题，它等效地表达为以下泛函形式：式中，T∞为环境温度(K)；能使泛函式(2)的一阶变分为零的温度函数，必为此问题的正确解；C、输入带钢热连轧工艺参数输入的工艺参数包括钢种、化学成分、热物性参数、坯料尺寸、出炉温度、设备参数、粗轧和精轧轧制规程和层流冷却制度，所述的层流冷却制度包括前段冷却规程、后段冷却规程和分段冷却规程；D、计算带钢在除鳞区的温度除鳞区指从加热炉出口至轧前除鳞结束时带钢轧件经历的区域，计算带钢轧件在经历这一区域过程中的温度场；首先进行除鳞区设备参数的输入，设备参数包括除鳞前空冷距离和除鳞设备长度，然后进行出炉后的空冷温度场模拟，再进行除鳞区的水冷温度场模拟；带钢轧件在辊道上运行过程中，通过自由表面和外界进行热交换，主要换热方式有辐射和对流，带钢轧制时的对流传热为自然对流，在带钢轧件温度>800℃时它引起的温降仅为辐射温降的1％，因此，在高温区段空冷过程中，忽略热对流的影响；利用Stefan-Boltzmann定律，辐射换热的热交换系数如下：式中，S为Stefan-Boltzmann常数；ε为辐射率或黑度系数；Tsu为带钢表面温度；T∞为环境温度；在对称面上无热交换发生，为绝热边界条件，其热交换量为Q＝0；对于高压水除鳞过程，带钢表面温度Tsu、水压以及水量密度ω对其热交换系数hw(W/2CN103028615B权利要求书2/3页(m2·K))的影响较大，随着表面温度不同，各温度区间的传热现象也不同；对于这种水冷形式的计算公式:式中：r—水压影响系数，高压水除鳞时r>1；E、计算带钢在粗轧区的温度粗轧区包括带钢轧件从除鳞区出口到粗轧末道次出口中间经历的所有区域，该区域中带钢轧件经历了空冷、轧制前除鳞和轧制变形过程；其中空冷与轧前除鳞计算模型与公式(3)和(4)相同；轧制过程变形区中主要存在以下几种热现象：带钢轧件塑性变形转化成的变形热、带钢轧件与轧辊相对滑动产生的摩擦热、带钢轧件与轧辊接触产生的接触热传导；E1、变形热：把塑性功转变的热能看作是一种内热源，其变形热的表达式为：式中，－平均变形速率；－等效应力；η－塑性功转变为热的比例，取0.7～0.9；E2、摩擦热：在轧制过程中因轧辊与带钢轧件间的摩擦而产生的单位表面积、单位时间的摩擦功，由下式表示：式中