用户名: 密码: 验证码:
中厚板轧后冷却过程温度场解析解研究与应用
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
中厚板轧后冷却过程自动控制是控冷技术得到有效、稳定、精确实施的重要保障,其中冷却过程温度模型是提高温度控制精度的核心和关键。随着计算机技术的快速发展,数值解法得到了广泛的应用,目前温度场有限差分计算模型已成功实现了在线应用,温度控制精度有了明显提高。温度场解析解法由于只能用于求解比较简单的问题,对于实际冷却过程的复杂情况无能为力,因而在工程实际中极少应用。但是解析解具有不可替代的理论价值,在认识问题本质,分析不同变量之间的物理关系,以及作为标准解用于检验和发展数值解法等方面都起到至关重要的作用。为此,本文以中厚板轧后冷却过程为研究对象,从简化条件下的温度场解析解研究入手,对解析解描述的温度分布规律进行了深入研究,结合实际重点探讨表面温度、心层温度与平均温度之间的对应关系,并基于此,对在线温度模型及离线模拟优化模型提出了改进方案。主要研究进展如下:
     (1)针对中厚板轧后冷却过程的特点,详细研究冷却过程温度场各种解法,包括:一维、二维、三维解析解法和一维显式、隐式、C-N格式,二维显式、ADI格式以及三维显式格式的有限差分解法。调试开发这些解法的计算程序,在统一计算条件下,解析与差分程序的计算结果具有很好的一致性。在此基础上,以解析解为基准,研究不同差分格式的稳定性和计算误差,并进一步通过程序运行时间的分析,探讨适于在线应用的差分格式及步长。研究结果对于在线温度模型差分方法的选定以及模型精度与实时性的分析有重要的参考价值。
     (2)研究非稳态热传导过程正规阶段的温度分布规律,提出正规解概念。正规解的温度分布、心表温差及平均温度由表面温度和Bi数综合确定,在正规阶段,可通过正规解法来求解温度分布。给出正规阶段的判定方法:当正规解的计算偏差减小到工程误差允许范围以内时,标志导热过程进入正规阶段,并定义相应的无量纲时间Fo_m作为非正规阶段和正规阶段的分界点。分析表明,Fo_m并不是某一定值,具体根据偏差要求、正规解计算偏差初值和Bi值共同确定。
     (3)针对边界条件和物性条件变化的实际情况,提出近似正规阶段概念。同样通过分析正规解计算偏差确定Fo_m值来判定近似正规阶段,在近似正规阶段,可应用正规解近似描述钢板的温度分布。经分析证实,钢板在水冷阶段前后的空冷过程中存在近似正规阶段,并且该阶段里的正规解温度分布是每一时刻实际温度分布的可能下限。特别得出,返红结束时刻可代替Fo_m值作为判定到达近似正规阶段的标志。
     (4)结合解析解描述的温度分布规律,研究钢板水冷后返红过程发生机理,模拟计算出不同钢种规格钢板在不同冷却条件下的返红结束位置,为终冷温度控制及测温仪位置确定提供参考依据;同时得出较厚钢板需在矫直之前停留时间的参考值。进一步针对实际的返红过程,研究钢板内部温度值等于平均温度的位置变化情况,得出返红发生的明显区域为温度值等于平均温度的稳定位置至表面,即无量纲位置X=0.57~1部分,这一结果为分析返红过程温度回复对钢板组织性能的影响提供了参考数据。此外,研究钢板在水冷及之后返红过程心表温差和均表温差的变化规律,为合理制定冷却策略以及返红结束后心表温差和均表温差的估算提供理论指导。
     (5)提出改进中厚板轧后冷却过程控制系统的思路和方法。确定返红结束后就近测温仪的实测表面温度作为终冷温度控制依据;将正规解计算嵌入到在线模型中,为模型设定计算和自学习计算提供较可信的温度分布和平均温度,可提高模型预报精度;离线模型扩展三个新的功能:返红过程分析、心表温差和均表温差分析以及正规解偏差分析,为在线模型控制精度提高起到促进作用。
     本文基于温度场解析解所做的研究工作,对于中厚板轧后冷却过程温度模型精度的提高有重要的理论指导意义。其中定义的正规解为实时深入了解物体内部的温度分布提供了可靠的计算方法,具有实用价值。在中厚板生产过程中,中间坯在粗轧、精轧之间待温过程,钢板在冷床上的空冷过程等,都具备条件应用正规解法计算温度分布。钢板越厚,采用正规解法计算温度分布的意义越大。
Automation control system for accelerated cooling process of plate is essential for controlled cooling technology to put into practice effectively, stably and accurately. Specifically, accuracy of temperature control for cooling process is subject to mathematical models. With the development of computer technology, numerical solutions have been applied widely. At present, the model based on finite difference method has been used on line to calculate temperature evolution, which increased accuracy of temperature control significantly. Due to the complexity of actual accelerated cooling process, analytical solution can be exactly solved only by simplifying conditions of the problem, which has limited its application in industrial practice. However, analytical solutions have irreplaceable theoretical value. They play very important roles in understanding essence of the problems, analyzing physical relation between different variables and acting as standard solutions to verify and develop numerical solutions. Therefore, concentrating on accelerated cooling process of plate, and commencing on the study of analytical solution for temperature field under the simplified conditions, the regularities of temperature distribution is researched in detail. In combination with the practical cooling process, the relation of surface, center and average temperatures is discussed. Based on the researches in this dissertation, improved projects for on-line temperature calculation model and off-line simulation model are put forward. The main progresses are as follows:
     (1) To obtain the temperature field for accelerated cooling process of plate, the solving methods consisting of one, two and three dimensional analytical solutions, as well as finite difference solutions including one-dimensional explicit, implicit and C-N (Crank-Nicolson) scheme, two-dimensional explicit and ADI (Alternating Direction Implicit) scheme, and three-dimensional explicit scheme are researched in detail. Computer programs of all these solutions are developed, and the results of them are good agreement under identical computational conditions. Then, compared with analytical solution, the computational stability and truncation error of different finite difference schemes with various step-sizes are investigated. The scheme of finite difference method adequate for on-line calculation with appropriate spacing and time step is discussed. The conclusion provides significant reference for analyzing the accuracy and real-time performance of on-line temperature model based on finite difference method.
     (2) Through analyzing the characteristic of temperature field in regular regime of unsteady heat conduction process, the concept of regular solution is defined. Regular solution depends on surface temperature and Bi number, and it is analytical solution in regular regime. Regular regime will be indicated by the calculation deviation of regular solution within the accepted tolerances. And Fo_m number, namely, the dimensionless time is defined as the dividing point between initial regime and regular regime. It is obtained that Fo_m is not a fixed value corresponding to same calculation deviation but depends on Bi number and the initial value of calculation deviation of regular solution.
     (3) In consideration of the actual variations of boundary conditions and thermophysical properties, the concept of approximate regular regime is defined. Approximate regular regime will also be indicated by the calculation deviation of regular solution within the accepted tolerances or by Fo_m. In approximate regular regime, regular solution is approximate to actual temperature field. It is proved that approximate regular regime exists in two air cooling processes which before and after water cooling process. And the regular solution of air-cooling in this regime is lower limit of the actual temperature field. In particular, it is drawn that the time of surface temperature recovery finishing after water-cooling can be considered as the sign of arriving in approximate regular regime instead of Fo_m .
     (4) Combined with the regularities of temperature distribution described by analytical solution, the mechanism of temperature recovery after water-cooling is introduced. The position of surface temperature recovery finishing is obtained by simulated calculation given different thickness, steel grade and cooling conditions, which provides reference for the control of finish cooling temperature and the location of pyrometers. And the reference values of stay time for thicker plates before leveling are presented. Furthermore, the position where the value of temperature equals to average temperature is analyzed according to temperature recovery process after water-cooling. It shows that the region of X = 0.57~1 is the distinct area of temperature recovery which will provide reference data to discuss the influence of temperature recovery on the structure and performance of plate. In addition, the difference between center and surface temperatures and that between average and surface temperatures in the water-cooling process and temperature recovery process after that are researched. The conclusion will provide theoretical direction to establish cooling strategies and estimate the difference between center and surface temperatures and that between average and surface temperatures after temperature recovery finishing.
     (5) The improved projects of process control system for accelerated cooling of plate are presented. Surface temperature measured by the pyrometer which close to the position of temperature recovery finishing will be the basis of finishing cooling temperature control. The regular solution is applied into on-line model, and provides reliable temperature distribution and average temperature to the preset model and self-learning model to increase the preset accuracy essentially. Three new functions can be extended into off-line model to promote the control accuracy of on-line model, which involve analysis of temperature recovery process after water-cooling, estimate of the difference between center and surface temperatures and that between average and surface temperatures, as well as research on the calculation deviation of regular solution.
     The researches on analytical solution for temperature field are of great theoretical significance to improve accuracy of temperature control in accelerated cooling of plate. The regular solution defined in this dissertation creates a reliable calculation method for real-time acquiring temperature distribution of things, which exhibiting significant utility value. In plate mills, holding process of inter-slab between rough rolling and finish rolling and cooling process of plate on cooling bed, etc., both meet conditions to apply regular solution to calculate temperature distribution. Application of the regular solution is more meaningful for thicker plate.
引文
[1]孙本荣,王有铭,陈瑛.中厚钢板生产[M].北京:冶金工业出版社,1993.
    [2]王占学.控轧控冷与控制冷却[M].北京:冶金工业出版社,1988.
    [3]王有铭,李曼云,韦光.钢材的控制轧制和控制冷却[M].北京:冶金工业出版社,1995.
    [4]李曼云,孙本荣.钢的控制轧制和控制冷却技术手册[M].北京:冶金工业出版社,1990.
    [5]小指军夫著,李伏桃等译.控制轧制控制冷却-改善材质的轧制技术的发展[M].北京:冶金工业出版社,2002.
    [6]王国栋,刘相华.日本中厚板生产技术的发展和现状(一)-随中国金属学会代表团访问日本观感之二[J].轧钢,2007,24(2):1-5.
    [7]王笑波,王仲初.中厚板轧后控制冷却技术的发展及现状[J].轧钢,2000,17(3):44-47.
    [8]张松亭,孙跃洲.厚板在线快速冷却技术和设备[J].宽厚板,1996,2(4):1-9.
    [9]陈松年.快速冷却和直接淬火技术[J].宽厚板,1996,2(5):1-8.
    [10]蔡晓辉.中厚板控制冷却数学模型的研究及应用[D].东北大学,2004.
    [11]龚彩军.中厚板控制冷却过程温度均匀性的研究[D].东北大学,2005.
    [12]Bernhard Ludwig.Systems for the accelerated cooling of plates[J].MPT,1988(4):10-17.
    [13]S.J.Chen,A.A.Tseng.Spray and jet cooling in steel rolling[J].International Journal of Heat and Fluid Flow.1992,13(4):358-369.
    [14]Fred C.Kohring.Waterwall water-cooling systems[J],Iron and Steel Engineer,1985,62(6):30-36.
    [15]T.Hope.Water curtain cooling systems in hot strip mills[J].Steel Time International,1984(4):14-22.
    [16]常崇民,陈继虞.气雾式在线冷却装置在厚板轧机上的应用[J].轧钢,1997(6):45-47.
    [17]万雪明.ADCO系统的冷却控制策略[J].轧钢,1999(6):16-18.
    [18]王崇尚.中厚板轧后在线快速冷却系统ADCO的自动控制[J].冶金自动化,2002(4):64-66.
    [19]J.Chen,T.Nijhuis.热轧机的新型冷却系统[J].钢铁,2005,40(11):84-86.
    [20]N.Champion.西门子奥钢联的中厚板轧机技术与能力[J].钢铁,2006,41(4):88-90.
    [21]刘相华,余广夫,焦景民等.超快速冷却装置及其在新品种开发中的应用[J].钢铁,2004,39(8):71-74.
    [22]王昭东,袁国,王国栋等.热带钢超快速冷却条件下的对流换热系数研究[J].钢铁,2006,41(7):54-64.
    [23]Y.V.Leeuwen,M.Onink,J.Sietsma.The γ-α transformation kinetics of low carbon steels under ultra-fast cooling conditions[J].ISIJ International,2001,41(9):1037-1046.
    [24]G.Buzzichelli,E.Anelli.Present status and Perspectives of European research in the field of advanced structural steels[J].ISIJ International,2002,42(12):1354-1363.
    [25]Aurelie Lucas,Pierre Simon,Guillaume Bourdon,et al.Metallurgical aspects of ultra fast cooling in front of the down-coiler[J].Steel research international,2004,75(2):139-146.
    [26]袁伟刚.日本JFE公司高性能中厚板生产技术介绍[J].冶金管理,2006(10):49-51.
    [27]Hiroshi Kagechika.Recent progress and future trends in the research and development of steel[J].NKK Technical Review,2003,22:6-9.
    [28]王国栋.均匀化冷却技术与板带材板形控制[J].上海金属,2007,29(6):1-5.
    [29]王国栋,刘相华,王君.我国中厚板生产设备、工艺技术的发展[J].中国冶金,2004(9):1-8.
    [30]孔金满.钢铁厂自动化系统的开发[J].鞍钢技术,1997(5):32-37,41.
    [31]孙一康.带钢热连轧的模型与控制[M].北京:冶金工业出版社,2002.
    [32]杨节.轧制过程数学模型(修订版)[M].北京:冶金工业出版社,1993.
    [33]丁修堃.轧制过程自动化[M].北京:冶金工业出版社,1986.
    [34]金以慧.过程控制[M].北京:清华大学出版社,1993.
    [35]何纯玉,吴迪,王君等.中厚板轧制过程计算机控制系统结构的研制[J].东北大学学报,2006,27(2):173-176.
    [36]F.Buschmann,A.Geisler,T.Heimake,et al.Framework-based software architectures for process automation systems[J].Annual Reviews in Control,2000,24:163-175.
    [37]D.Auzinger.Recent development in process optimization for laminar cooling in hot strip mills[J].Ironmaking and Steelmaking,1996,23(1):84-87.
    [38]K.Yahiro,J.Yamasaki,M.Furukawa,et al.Development of coiling temperature control system on hot strip mill[J].Kawasaki Steel Technical Report,1991,24(4):32-40.
    [39]T.Hope.Control of run-out table cooling[J].Steel Time International,1994(1):21-23.
    [40]D.Auzinger,F.Parzer and G.Posch.Process optimization for laminar cooling[J].MPT International,1996(5):68-75.
    [41]L.Henrich,R.Holz,G.Kneppe.Physically based cooling line model to meet growing demands in temperature and flexibility[J].Ironmaking and Steelmaking,1996,23(1):79-81.
    [42]韩曾晋.自适应控制[M].北京:清华大学出版社,1995.
    [43]胡贤磊.中厚板轧制过程控制模型的研究[D].东北大学,2003.
    [44]B.Fazan,D.Boubel,P.Ratte,et al.Optimum computer control of a plate mill[J].Iron and Steel Engineer,1980,57(11):58-64.
    [45]A.F.MacAlister.Modeling and adaptive techniques for rolling mill automation[J].Iron and Steel Engineer,1989,66(12):38-46.
    [46]P.A.Atack,I.S.Robinson.Adaptation of hot mill process models[J].Journal of Materials Processing Technology,1996,60(6):535-542.
    [47]J.F.Evans,L.D.Roebuck,H.R.Watkins.Numerical modeling of hot strip mill runout table cooling[J].Iron and Steel Engineer,1993,70(1):50-55.
    [48]W.J.Lawrence,A.F.MacAlister,P.J.Reeve.On line modeling and control of strip cooling [J].Ironmaking and Steelmaking,1996,23(1):74-78.
    [49]姜启源,谢金星,叶俊.数学模型(第三版)[M].北京:高等教育出版社,2003.
    [50]王洪.中板在线加速冷却装置及工艺[J].轧钢,1992(4):31-34.
    [51]杨增堂,董希满,蔡庆伍等.中厚板水幕冷却温度预报和冷却参数设定模型的研究[J].钢铁,1996,31(9):34-38.
    [52]冯光宏,蔡庆伍,宋佩莼.中厚板轧后水幕冷却控制模型的研究[J].钢铁,1997,32(7):27-30.
    [53]侯安全,蔡庆伍,米振莉等.中厚板轧后高密度管层流快冷控制模型[J].北京科技大学学报,2002,24(5):544-546.
    [54]蔡晓辉,张殿华,谢丰广等.中厚板控制冷却数学模型[J].东北大学学报,2003,24(10):923-926.
    [55]刘华强,唐荻,杨荃等.中厚板轧后加速冷却过程的温度场模型[J].上海金属,2007,29(1):48-51.
    [56]C.Devadas,I.V.Samarasekera.Heat transfer during hot rolling of steel strip[J].Ironmaking and Steelmaking,1986,13(6):311-321.
    [57]J.Filipovic,Raymond Viskanta,E P.Incropera,et al.Thermal behaviour of a moving steel strip cooled by an array of planar water jets[J].Steel Research,1992,63(10):438-446.
    [58]N.Hatta,Y.Yanaka,H.Takuda,et al.A numerical study on cooling process of hot steel plates by a water curtain[J].ISIJ International,1989,29(8):673-679.
    [59]N.Hatta,H.Osakabe.Numerical modeling for cooling process of a moving hot plate by a laminar water curtain[J].ISIJ International,1989,29(11):919-925,
    [60]M.M.Prieto,L.S.Ruiz,and J.A.Menendez.Thermal performance of numerical model of hot strip mill runout table[J].Ironmaking and Steelmaking,2001,28(6):474-480.
    [61]K.Laber,H.Dyja,D.Rydz.Analytical and numerical methods of determining the distribution of temperature of air-cooled strip[J].Metalurgija.2005,44(1):31-35.
    [62]张鹏云,徐维峰.ACC的应用及功能开发[J].宽厚板,2005,11(3):6-9.
    [63]詹胜利,张华等.中厚板轧后控冷系统及其应用[J1.轧钢,2005,22(4):33-35.
    [64]杨世铭,陶文铨.传热学(第三版)[M].北京:高等教育出版社,1998.
    [65]赵镇南.传热学[M].北京:高等教育出版社,2002.
    [66]贾力,方肇洪,钱兴华.高等传热学[M].北京:高等教育出版社,2003.
    [67]陶文铨.计算传热学的近代进展[M].北京:科学出版社,2000.
    [68]陶文铨.数值传热学(第二版)[M].西安:西安交通大学出版社,2001.
    [69]曹玉璋,邱绪光.实验传热学[M].北京:国防工业出版社,1998.
    [70]辛放斌.材料成形计算机模拟[M].北京:冶金工业出版社,2006.
    [71]陆全康,赵蕙芬.数学物理方法[M].北京:高等教育出版社,2003.
    [72]梁昆淼.数学物理方法(第二版)[M].北京:人民教育出版社,1978.
    [73]俞昌铭.热传导及其数值分析[M].北京:清华大学出版社,1981.
    [74]陆金甫,关治.偏微分方程数值解法(第二版)[M].北京:清华大学出版社,2004.
    [75]李立康,於崇华,朱政华.微分方程数值解法[M],上海:复旦大学出版社,1999.
    [76]刘庄,吴肇基,吴景之等.热处理过程的数值模拟[M].北京:科学出版社,1996.
    [77]孔祥谦.有限单元法在传热学中的应用(第三版)[M].北京:科学出版社,1998.
    [78]D.R.克罗夫特著,张风禄译.传热的有限差分方程计算[M].北京:冶金工业出版社,1982.
    [79]陆金甫,顾丽珍,陈景良.偏微分方程差分方法[M].北京:高等教育出版社,1988.
    [80]胡祖炽,雷功炎.偏微分方程初值问题差分方法[M].北京:北京大学出版社,1988.
    [81]A.B.雷柯夫著,裘烈钧,丁履德译.热传导理论[M].北京:高等教育出版社,1955.
    [82]B.П.伊萨琴科等著,王丰,冀守礼,周筠清等译.传热学[M].北京:高等教育出版社,1987.
    [83]E.R.G.埃克尔特,R.M.德雷克著,徐明德译.传热与传质[M].北京:科学出版社,1963.
    [84]钱壬章,俞昌铭,林文贵.传热分析与计算[M].北京:高等教育出版社,1987.
    [85]辛荣昌,陶文铨.非稳态导热充分发展阶段的分析解[J].工程热物理学报,1993,14(1):80-83.
    [86]胡贤磊,李建民,王昭东等.中厚板轧制过程温度计算的二次曲线建模法[J].轧钢,2002,19(6):12-14.
    [87]韩斌,周存龙,王国栋等.带钢控冷过程温度计算二次曲线方法分析与应用[J].东北大学学报,2005,26(10):975-978.
    [88]王国栋,吴国良等译.板带轧制理论与实践[M].北京:中国铁道出版社,1990.
    [89]Y.S.Touloukian.Specific heat(Metallic elements and alloys)-Thermophysical properties of Matter(Volume 4)[M].New York:IFI/Plenum Press,1970.
    [90]Goldsmith,Alexander.Handbook of thermophysical properties of solid materials(Volume 1and 2)[M].London:Pergamon,1961.
    [91]中国金属学会,中国有色金属学会编.金属材料物理性能手册(第一册)--金属物理性能及测试方法[M].北京:冶金工业出版社,1987.
    [92]谭真,郭广文.工程合金热物性[M].北京:冶金工业出版社,1994.
    [93]谭延昌.金属材料物理性能测量及研究方法[M].北京:冶金工业出版社,1989.
    [94]宋学孟.金属物理性能分析[M].北京:机械工业出版社,1981.
    [95]徐京娟,邓志煜,张同俊.金属物理性能分析[M].上海:上海科学技术出版社,1988.
    [96]王润.金属材料物理性能[M].北京:冶金工业出版社,1993.
    [97]田莳.材料物理性能[M].北京:北京航空航天大学出版社,2004.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700