Mathematical Problems in Engineering

Volume 2017 (2017), Article ID 4091342, 6 pages

https://doi.org/10.1155/2017/4091342

## Optimization of Roller Velocity for Quenching Machine Based on Heat Transfer Mathematical Model

School of Environment and Municipal Engineering, Qingdao University of Technology, Qingdao 266033, China

Correspondence should be addressed to Jingkui Zhang

Received 28 May 2017; Accepted 22 June 2017; Published 6 August 2017

Academic Editor: Marcello Vasta

Copyright © 2017 Yunfeng He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

During quenching process of steel plate, control parameters are important to product quality. In this work, heat transfer mathematical model has been developed for roller-type quenching machine to predict the temperature field of plate at first, and then an optimization schedule considering quenching technology and equipment limitations is developed firstly based on the heat transfer mathematical model with considering the shortest quenching time. A numerical simulation is performed during optimization process to investigate the effects of roller velocity on the temperature of representative plate. Based on the optimization method, study is also performed for different thickness of plate to obtain the corresponding roller velocity. The results show that the optimized roller velocity can be achieved for the roller-type continuous quenching machine based on the heat transfer mathematical model. With the increasing of plate’s thickness, the optimized roller velocity decreases exponentially.

#### 1. Introduction

Quenching is a common manufacturing process to obtain unique service performance properties of metal by controlling the amount and the morphology of microstructural constituents [1]. Roller-type continuous quenching machine is an important device which is commonly used in quenching heat treatment of plate in the iron and steel industry. In recent years, a number of works have been conducted on quenching machine referring to the heat transfer coefficient, critical cooling rate, plate temperature field, flow field, stress distribution, and microstructure distribution. Tang et al. [2], Şimşir and Gür [3], and Huiping et al. [4] simulated the temperature, stress, and microstructure field considering nonlinear material properties using finite element method during the quenching process. Kang and Im [5] established a 3D thermo-elastic-plastic model for plain-carbon steel along with phase transformation to improve the accuracy of numerical simulation. Zhu et al. [6] developed a three-dimensional finite element flow field model of spraying pipe of roller-type quenching machine. Chantasiriwan [7] treated the problem of estimating the transient heat transfer coefficient at the surface of steel bars subjected to quenching using the sequential function specification method. Wang et al. [8] conducted a study on the heat transfer during quenching for plate roller quenching machine. Fu et al. [9] studied the critical cooling rate of plate during quenching. In their study, three types of quenching critical cooling rate models of medium carbon steel, low-carbon steel, and low-alloy construction steel were established with the modified Maynier model, Ei-dis model, and isothermal curve model, respectively.

In recent years, few researches on control parameters’ optimization of continuous quenching machine were performed. Wang et al. [10] proposed a calculation method of optimal speed profile for online control system. Lü et al. [11] carried out a calculation model based on coupled metallo-thermo-mechanics theory. With the simulated results, they proposed an optimum quenching scheme. Huiping et al. [12, 13] performed a step function model with respect to optimized technologic parameters during gas quenching process. The heat transfer plays an important role on the quenching heat treatment. However, studies on the optimization of control parameters based on the heat transfer model are few for the continuous quenching machine up to now. Therefore, in this work, the heat transfer mathematical model is developed during quenching process firstly. And then, an optimization method of roller velocity is proposed based on the model considering the quenching technology and water consumption. The optimized roller velocity is also obtained for different plate thickness.

#### 2. Heat Transfer Mathematical Model

Heat transfer during quenching occurs via all possible heat transfer mechanisms, that is, conduction, convection, and radiation. Basically, heat is removed from the surface of the plate by convection to the water on the boundary during quenching process, which results in thermal gradients driving the conduction inside the plate. Because the surface temperature of the plate is higher than 800°C before quenching and higher than 100°C when the plate leaves the high-pressure water injection zone, the form of heat transfer between plate surface and cooling water is big tank boiling convection at the high-pressure water injection zone. The cooling water sprays on the surface of hot plate at constant velocity and angle. The result is that water near the hot plate’s surface vaporizes quickly, taking the heat away in the form of latent heat of vaporization. The water away from the steel surface takes away the sensible heat through conduction and convection. When the plate passes the low-pressure water injection zone, the main heat transfer form is forced convection.

##### 2.1. Governing Equations

The transient heat transfer within the plate during quenching can be described mathematically by an appropriate form of Fourier’s heat conduction equation (1). The equation can be expressed in its most general form aswhere is the thermal conductivity. is the density of material. is the constant pressure specific heat, and is quenching time.

##### 2.2. Initial and Boundary Conditions

The initial conditions at time can be described aswhere is the initial temperature of plate.

The boundary conditions are related to the heat transfer between quenching plate and cooling water. The relevant boundary conditions are given aswhere is the heat flux of plate’s surface. is the total heat transfer coefficient including convection and radiation. and are the thickness and width of the plate. is the total quenching time. is the water temperature.

The total heat transfer coefficient is related to the flow flux and velocity of the cooling water, the surface temperature, roughness and emissivity of the plate, and so on. The empirical equation is employed to express the total heat transfer coefficient which is available in [14, 15], and it can be written aswhere , , and are regression coefficients. is the flow flux of the cooling water on the plate surface. is the surface temperature of the plate.

##### 2.3. Numerical Method and Codes Validation

The heat transfer mathematical model developed in this work can be used to predict the evolution of temperature during quenching process for the roller-type quenching machine. The governing equation (energy equation) with the boundary conditions is solved by finite volume method (FVM). The FVM is a common method to calculate the heat conduction problem. The codes validation on solving of energy equation has been done by the FVM in our former works which are available in [16, 17].

#### 3. Optimization Method

##### 3.1. Objective Function

The cooling rate of plate during quenching process is determined by the roller velocity under constant water injection condition. Therefore, the roller velocity is the direct factor for quenching. The optimization of roller velocity plays an important role for the quenching machine. In this work, the optimization method of roller velocity is developed on the basis of the heat transfer mathematical model. The shortest quenching time is set as the optimization objective. Because the water injection condition is constant when the roller velocity is optimized, the shortest quenching time corresponds to the least water consumption. The objective function is given aswhere and are the length of high-pressure water injection zone and low-pressure water injection zone of the quenching machine. is the roller velocity.

##### 3.2. Constraints

Constraints mainly include the quenching technology requirements and quenching machine equipment limitations. They can be described aswhere and are the actual cooling rate of plate from 800°C to 500°C and critical cooling rate, respectively. and are the maximal velocity and the minimal velocity of roller. is the temperature at the centre of plate.

As the flow diagram shown in Figure 1, the roller velocity is optimized with the objective function in meeting the constraint conditions. Quenching under the optimized roller velocity corresponding to the shortest quenching time is also a process consuming the least water in the case of constant water injection condition.