32RM92.中厚板轧制过程高精度温度数学模型研究-RAL刘刚

2007年全国塑性加工理论与新技术学术研讨会2007年5月沈阳183中厚板轧制过程高精度温度数学模型研究*刘刚,李长生，胡贤磊，刘相华（东北大学轧制技术及连轧自动化国家重点实验室，辽宁沈阳110004）摘要：中厚板轧制过程中温度数学模型对于过程设定起到重要作用，采用有限元方法计算温度具有较高的计算精度。本文采用有限元方法开发了轧制过程温度数学模型，通过对中厚板厂实际生产过程中温度场的计算验证了该方法具有较高的精度。为了实现该技术的应用，提出了有限元温度数学模型的在线方案。关键词：温度场；数学模型；有限元；中厚板；在线应用High-PrecisionTemperatureMathematicsModelforPlateRollingProcessLIUGang,LIChangSheng,HUXianLei,LIUXianghua(StateKeyLaboratoryofRollingandAutomation,NortheasternUniversity,Shenyang110004,China)Abstract:Temperaturemathematicsmodelhasimportantinfluenceonprocesscontrolduringtheplaterollingprocess.HighpreciseresultscanbeobtainedwhiletemperatureiscomputedbyFEM.AtemperaturemathematicsmodelwasdevelopedontheplaterollingprocessbyFEMinthepaper.Whilethetemperaturefieldofapracticalplaterollingprocesswasanalysedthroughthemathematicsmodel,thecalculatedresultsshowhighprecisioncomparedwiththemeasuredtemperature.AnonlineapplicationschemeofFEMmathematicsmodelwasproposedfinally.Keywords:temperaturefield;mathematicsmodel;finiteelementmethod;plate;onlineapplication1前言在中厚板生产过程中，温度是最重要的参数之一。由于温度直接影响到轧制力，因此精确预报生产过程中各阶段的温度是保证厚度、宽度和板形等数学模型命中率的关键；另外，温度也为产品组织性能的预测和控制提供了必要的数据支持。过去中厚板生产过程中使用的温度计算模型或受其求解机理限制或过于简化，难以适应生产环境的变化，计算偏差较大。因此，开发高精度的温度数学模型具有十分重要的意义。作为一种数值计算方法，有限元法能够准确地计算轧制过程中轧件的温度场。过去由于受有限元法计算时间的限制和有限元计算瞬态温度场时震荡现象等原因的影响，在线应用受到限制。随着计算机技术的迅猛发展，有限元分析瞬态温度场基本理论的完善和一些针对有限元方法快速算法的提出[1～4]，用有限元法在线计算热轧过程中的温度成为可能。本文以传热学基本理论为基础,采用有限元方法开发了中厚板轧制过程中二维温度高精度数学模型；该模型可以计算出任意时刻钢板横断面上各点的温度值。以某中厚板厂实际生产过程为例计算了中厚板温度场的演变，并与实测温度进行了比较。为了实现高精度模型技术的应用，讨论和制定了在线方案。2数学模型2.1有限元方程的建立考虑到中厚板生产的实际情况，在计算时采用如下假设：*基金项目：国家自然科学基金重点项目,项目编号：505340202007年全国塑性加工理论与新技术学术研讨会2007年5月沈阳1841）几何对称面两侧的温度分布对称，在对称面上没有热量交换；2）钢板的黑度考虑成钢板厚度的函数[5]：ε=f1(H)（1）式中，ε为钢板的黑度，H为钢板的厚度；3）钢板的密度和比热按钢种考虑成钢板温度的函数：ρ=f2(T)c=f3(T)（2）式中，ρ钢板的密度，c钢板的比热，T为钢板的温度；4）忽略因摩擦力引起的钢板表面生热；5）钢板与工作辊接触时，工作辊中心温度保持不变。将变形热视为内热源，并将自由表面和接触面都归为第三类边界条件，二维非稳态导热问题的数学描述为：tTkqyTxT∂∂=+∂∂+∂∂⋅α12222（3）0)()(=−+∂∂+∂∂∞TThlyTlxTkyx（4）ckρα/=（5）式中，ρ为板坯的密度；c为板坯的比热容；k为导热系数；α为导温系数；⋅q为内热源强度；xl,yl为边界外法向的方向余弦；h为板坯与周围流体间的换热系数；T为温度；∞T为周围流体的温度；t为时间。利用欧拉方程上述导热问题可以等效为以下泛函求极值问题[6]：∫∫⎥⎥⎦⎤⎢⎢⎣⎡⎟⎟⎠⎞⎜⎜⎝⎛∂∂−−⎥⎥⎦⎤⎢⎢⎣⎡⎟⎟⎠⎞⎜⎜⎝⎛∂∂+⎟⎠⎞⎜⎝⎛∂∂=dsTtTcqyTxTkISρ&22122()∫∞−+dlTThL321（6）令所研究区域所划分的E个单元的泛函一阶变分分别为零，并将所有单元组装成整体表达式可得：[]{}[]{}ptTKTKT=⎭⎬⎫⎩⎨⎧∂∂+3（7）式中[]TK为温度刚度矩阵，[]3K为热容矩阵，{}p为载荷矩阵。2.2数值震荡问题对于（7）式，往往通过在时间域内使用有限差分法进行求解。其求解过程为从初始温度场开始,每隔一个时间步长，求解下一时刻的温度场，一步一步向前推进,称为步进积分。该方法的优点是节省计算机内存，但是容易产生振荡现象。由于震荡常常使计算结果偏离解析解，为了提高计算精度，应加以解决。关于震荡产生的原因，国内外的学者进行了广泛而深入的研究，并给出了相应的降低或避免震荡的对策。经过对各种对策进行编程比较发现采用集中热容矩阵的向后差分格式可以有效的抑制震荡现象，而且计算效率较高。采用向后差分法时（7）式可以转化为下面的有限元方程：[][]{}[]{}{}pTKtTKtKtttT+Δ=⎟⎟⎠⎞⎜⎜⎝⎛⎥⎦⎤⎢⎣⎡Δ+Δ−3311（8）2007年全国塑性加工理论与新技术学术研讨会2007年5月沈阳185其中，tΔ为迭代时间步长。2.3边界条件(1)辊道上的待温过程在待温过程中，轧件自由表面与外界的热交换方式主要有热辐射和热对流两种，在此期间的综合热交换系数可表示为[7]：[]))(()(0,202,3/10,TTTTTThjijijia+++−=εσα（9）其中,α为修正系数，jiT,为节点处的温度；ε为钢材的辐射率；σ为斯蒂芬-玻尔兹曼常数；0T为室内温度。(2)除鳞过程在高压水除鳞过程中，带钢表面温度、水压以及水量密度对热交换系数的影响较大。文献[8]根据大量数据回归出喷水冷却在各温度区间的热交换系数:163.1102.10700147.0663.0×××=−sTwrhω（10）其中，sT为钢板表面温度，r水压影响系数，ω为水量密度。(3)轧制过程在轧制过程中，水平表面与轧辊表面接触，垂直表面通过对流和辐射而冷却。在钢材表面与轧辊发生接触时，总的热交换系数可通过式(11)来计算[7]:)/(2απβtkhr=（11）其中，β为修正系数，t为轧件与轧辊接触时间。2.4塑性变形升温单位时间、单位体积内因塑性应变产生的热量,取决于钢材的瞬时屈服应力和应变。单位时间、单位体积的塑性变形功转化的热量可以表示为：∫∫•−−=dSWeSeεση)((12)其中，η为修正系数，−σ为单元的平均变形抗力，•−ε为单元的平均变形速率。2.5网格划分为了节约计算时间，考虑到钢板的对称性，取其横断面四分之一进行计算。采用四边形等参单元，可以根据精度要求输入钢板宽度方向的单元数m1和厚度方向的单元数m2，程序将自动划分单元,并对单元进行编号。网格划分得越密，计算精度越高，同时计算时间也越长。根据生产的实际需要，对于一次轧制过程，模型的计算时间不可超过0.5秒。因此，网格的划分也要考虑到生产实际的需要。3实验及结果分析为了验证模型的准确性和实用性，选取某厂生产的中厚板为研究对象加以说明。其工艺参数如下：钢种为Q235B，初始板坯厚度为220mm，轧制后钢板厚度为20.5mm，出加热炉温度为1200℃，钢坯从加热炉到除鳞机的运行时间为22s，除鳞时间为3s，工作辊半径为525mm，计算中用到的其2007年全国塑性加工理论与新技术学术研讨会2007年5月沈阳186它轧制工艺制度如表1所示。表1轧制工艺制度Table.1rollingpassschedule模型的开发语言为Fortran77，模型的运行环境为：FortranPowerStation4.0软件平台，CPU主频3.06MHZ，内存512M。长度方向取20个单元，宽度方向上取10个单元。空冷过程每次迭代时间取1s，除鳞过程每次迭代时间取0.1s，轧制过程每次迭代时间取轧件与轧辊接触时间的十分之一。道次平均运行时间为0.265秒，满足工程计算的需要。图1板坯温度分布的演变Fig.1Evolutionoftemperaturedistributionofslab102010301040105010601070108010901100020030040050060070080090计算温度/℃图2各道次轧制温度计算值与实测值的比较Fig.2Thecomparisonofrollingtemperatureeachpassbetweencalculatedandmeasured道次123456789101112间隔时间/s3988915779912711实测出口厚度/mm20118216242114.588.567.250.838.228.622.620.5轧辊转速/（r·min-1）4030393550476060606062620204060801001201401601802008008509009500000500050200时间/s心部温度平均温度表面温度2007年全国塑性加工理论与新技术学术研讨会2007年5月沈阳1874模型的在线应用方案在模型应用之前，通过某中厚板厂大量实测温度数据对黑度和换热系数等影响计算精度的因子进行优化，保证一定的计算精度。在模型应用时，结合该厂的实际生产情况，在不影响生产正常进行的情况下制定了有限元温度数学模型的在线应用方案。由于该厂的计算机控制系统是二级计算机系统，可以由二级过程控制系统调入PDI数据和工艺参数，结合模型参数进行计算，并将计算结果用于轧制过程中的参数设定与控制。另外，为保证模型的计算精度，还应将有限元温度数学模型的计算结果和实测的温度值在线作进一步比较，对模型的相关参数进一步优化，确保模型具有较高的精度。该方案的数据流如图3所示。图3温度数学模型的数据流图Fig.3Dataflowoftemperaturemathematicsmodel5结语(1)采用有限元法建立了中厚板轧制过程中温度高精度数学模型。(2)对某中厚板厂实际生产中的温度场进行了计算，通过与实测温度进行比较，可以看出该模型具有较高的计算精度。(3)结合某中厚板厂的实际生产情况，提出了有限元温度数学模型的在线应用方案。参考文献[1]欧阳华江,等.一类抛物型方程有限元算法的计算准则应用数学和力学[J].应用数学和力学,1989,10(12):1115-1120.[2]牛山廷,等.淬冷过程三维温度场有限元模拟数值振荡问题的研究[J].金属热处理,2006,31(6):71-74.[3]马向平,等.瞬态温度场有限元法求解的研究[J].装甲兵工程学院报.2002,16(2):22-26.[4]张永杰,孙秦.有限元线性方程组改进的全稀疏存贮方法[J].航空计算技术.2006,36(2):117-119.[5]胡贤磊,等.中厚板精轧过程的高精度温度预测模型[J].东北大学学报(自然科学版).2003,24(1):71-74.[6]刘相华.刚塑性有限元及其在轧制中的应用[M].冶金工业出版社,1994.[7]周晓光,等.中厚板热轧