Strojniški vestnik - Journal of Mechanical Engineering 57(2011)10, 730-738 D01:10.5545/sv-jme.2010.075 Paper received: 09.04.2010 Paper accepted: 23.06.2011 Optimal Control of Workpiece Thermal State in Creep-Feed Grinding Using Inverse Heat Conduction Analysis Marin Gostimirovic1* - Milenko Sekulic1 - Janez Kopac2 - Pavel Kovac1 1 University of Novi Sad, Faculty of Technical Science, Serbia 2 University of Ljubljana, Faculty of Mechanical Engineering, Slovenia Due to intensive friction between grinding particles and workpiece material, a substantial quantity of thermal energy develops during grinding. Efficient determination ofreal heat loading in the surface layer of the workpiece material in grinding largely depends on the reliability of basic principles of distribution of heat sources and the character of the temperature field within the cutting zone. Therefore, this paper takes a different approach towards the identification of the thermal state of the creep-feed grinding process by using the inverse problem to approximate heat conduction. Based on a temperature measured at any point within a workpiece, this experimental and analytical method allows the determination of a complete temperature field in the workpiece surface layer as well as the unknown heat flux on the wheel/ workpiece interface. In order to solve the inverse heat conduction problem, a numerical method using finite differences in implicit form was used. When the inverse heat conduction problem is transformed into an extreme case, the optimization of heat flux leads to an allowed heat loading in the surface layer of workpiece material during grinding. Given the state function and quality criterion, the control of workpiece heat loading allows the determination of optimal creep-feed grinding conditions for particular machining conditions. ©2011 Journal of Mechanical Engineering. All rights reserved. Keywords: creep-feed grinding, head loading, inverse problem, optimal control 0 INTRODUCTION Grinding is considered as one of the most important machining methods. In addition to the conventional fine multi-pass grinding, there has been a recent introduction of high-productivity grinding methods. These high-productivity methods use higher cutting speeds and/or larger cutting depths in order to increase the relatively low productivity which has traditionally been considered the main drawback of conventional grinding. However, the increase of grinding conditions considerably changes the grinding kinematics, i.e. the conditions in the wheel/ workpiece interface [1]. High-speed grinding is characterized by lower cutting depths and shorter time of contact between the workpiece material and grinding particles, with a more intensive friction in the cutting zone. The increased contribution of friction in the generation of the total quantity of thermal energy also increases the contact temperature. The potential decrease of contact temperature during extremely high-speed machining can be attributed to a reduced cross-section of chip and faster introduction of the next grinding particle into the wheel/worpiece interface [2]. This improves cutting conditions because of reduced friction and deformation during chip forming due to interfacing between the grinding particles and the softened material layer generated by previous grinding. In creep-feed grinding, which is characterized by large cutting depths and small workpiece speeds, there is a longer wheel/ workpiece interface as well as a prolonged time of contact with workpiece material. At the same time, this results in the generation of more intensive heat sources in grinding particles and a prolonged time of their effect. In addition, a longer interface contributes to a better evacuation of heat from the cutting zone, which results in a decreased power of heat source per unit area [3] and [4]. Also, a small workpiece speed improves the distribution of total thermal energy within the cutting zone due to its prolonged effect. All this results in an increased total quantity of thermal energy per unit area but over a longer time period. *Corr. Author's Address: University of Novi Sad, Faculty of Technical Science, Trg D. Obradovica 6, 21000 Novi Sad, Serbia, maring@uns.ac.rs Based on the previous discussion, it is evident that high-productivity grinding methods cause the development of large quantities of thermal energy within the cutting zone [5]. The generated thermal energy, located within a relatively narrow area of the cutting zone, causes high cutting temperatures in creep-feed grinding. These increased temperatures instantaneously burst to a maximum, have short duration and exert a pronounced negative effect on wheel surface, workpiece quality and accuracy. Since the main task of grinding is to achieve satisfactory part quality with as large productivity as possible, special attention is focused on the effect that grinding temperatures have on the change of material properties in the workpiece surface layer. If the temperatures thus generated are high enough to cause structural and phase transformations of the workpiece material, the machined surface shall suffer from a number of disadvantages. Should, in addition, dimensional errors appear as well, the overall effect can substantially diminish exploitation features of the finished part. Efficient control of thermal phenomena in grinding, i.e. the determination of allowed heat loading on workpiece surface layer, requires knowledge of heat development and distribution in grinding, the temperature field in the cutting zone and, finally, the influence of cutting conditions on grinding temperature [6]. For that reason, identification of the thermal state in grinding based on analytical models and experimental results has been gaining popularity. As the research so far has shown, non-stationary and non-linear technical processes involving intensive heat conduction, such as creep-feed grinding, can be successfully solved using a novel approach based on inverse problem of heat transfer [7]. The inverse problem of heat transfer allows the closest possible experimental-model approximation of thermal regimes for grinding. In the case of control over the grinding thermal regime, the extreme case of the inverse problem of heat conduction [8] is practically the only way to reliably approximate the allowed heat loading on the workpiece surface layer. For a known temperature measured at a point within the workpiece surface layer, numerical methods are used to approximate the total temperature field of the surface layer as well as the unknown heat flux density on the wheel/workpiece interface [9] and [10]. For the selected quality criteria, particular machining conditions, and a predetermined loading on the workpiece surface layer, it is possible to arrive at optimal cutting conditions in creep-feed grinding by controlling the heat flux. 1 PARAMETERS OF GRINDING HEAT SOURCE The role of mathematical theory behind thermal phenomena in grinding is to adopt the most adequate model of the workpiece, grinding wheel and their inter-relationships, Fig. 1. It can be assumed that the elementary heat source on the grinding particle is the result of friction between the grinding particle, workpiece and chip in the workpiece material shear plane. Summing up all the heat sources, i.e. grinding particles in contact with the workpiece, gives the total heat source for the entire cutting zone, qT. This total heat source, whose strength varies within a narrow range, acts continuously, shifting across the workpiece surface with constant velocity [5] and [6]. Fig. 1. Model of thermal state in creep-feed grinding The heat source in grinding is characterized by its power and length of activity. Magnitudes of these two parameters depend on the machining process, wherein cutting conditions have the predominant influence. Heat source power in the cutting zone is the basic source parameter since it predominantly impacts the thermal state of the workpiece surface layer. An analytical definition of heat source power starts from the fact that in grinding the total mechanical work applied to the cutting process is transformed into thermal energy: Qt = Ft ■ vc • te (1) where Ft is the tangential grinding force, vc is the cutting speed, and te is the time of contact between the grinding particles and workpiece. The total heat quantity, Eq. (1), is transformed into heat flux density: It = Qt Ft- JJ dxdy J dt A„ (2) If, in Eq. (2), the wheel/workpiece interface is decomposed into Ac = lc x bc = (a x Ds)l/2 x bc, and the tangential grinding force is reduced to a unit grinding width bc, that is, we introduce the specific tangential grinding force F't = Ft / bc, then: It = Ft' vc (3) On the other hand, the specific tangential grinding force can be expressed by a general Eq.: F = k h sm m > (4) where hm = vw x a/vc is the mean cutting depth, and ksm is the mean specific cutting force. The final form for total heat flux density as a function of the grinding conditions can be derived by substituting Eq. (4) in Eq. (3): iT = ksm ■ D ' (5) On the other hand, the total active time interval of the heat source has a prominent influence on the heat quantity being evacuated from the cutting zone into the workpiece. The shorter active time of the heat source diminishes the heat quantity transferred via the wheel/ workpiece interface into the workpiece surface layer, while at the same time allowing a faster supply of coolant onto the freshly ground surface. This time interval is expressed as the ratio between the interface length lc and workpiece velocity vw, that is: 2 INVERSE PROBLEM OF THE GRINDING PROCESS The process of heat transfer between solid bodies or between a system and its environment, of which heat conduction in grinding is also a part, is mostly considered from the standpoint of mutual relations between the input and output process parameters. If the input parameters u(t) are known and the output parameters z(t) define the process state in time, then the output parameters are a function of the input parameters, i.e.: z = f (u, t). (7) If, for the adopted thermal model, there exist unique conditions, then the determination of the input-output relationship is the direct task of heat conduction. Conversely, the inverse heat conduction solves the problem of finding the input characteristics of the process for the known temperature field [7]. If for every unknown parameter u there is a linear, smooth operator A which allows the determination of the output parameter z, then the general case of the inverse problem is formulated by the following Eq.: Axu = z . (8) 2.1 Case of the Inverse Problem in Grinding In grinding, considering that the cutting depth is many times smaller than the length and width of the wheel/workpiece interface, the heat source can be treated as a strip of infinite length and constant heat distribution. Furthermore, if the dissipation of heat flow in the direction of heat source movement is disregarded, then the workpiece can be approximated with a series of adiabatic thin plates, Fig. 2. Substitution of the real workpiece with the semi-infinite plate is completely justified, bearing in mind that the heat source in grinding is generated within a small volume of workpiece material while the heat loading of the surface workpiece layer is considered depth-wise. v w Fig. 2. Schematic of a uni-dimensional inverse problem in grinding For such a defined thermal model of grinding, the following is the differential equation of a one-dimensional heat conduction, for x e (0,H) and t £ (0,tj: pc(ß) ™=*-{m ¿A dt d x I dx (9) where 6 = 6(x,t) is the workpiece temperature at point coordinate x at moment t, A = A(6) is the heat transfer coefficient, pc = C(6) is the specific heat capacity (p - material density, c - specific heat), H is the thickness of the surface layer of workpiece material, and tm is the largest time increment. Now the analytical form of the inverse problem of heat conduction for surface grinding can be described with Eq. (9), in conjunction with initial conditions, additional conditions and boundary conditions. The initial temperature distribution in the workpiece for the initial moment t = 0 is: O( x, * )|,=0 =H x). (l0) An additional condition is the fact that at the point x = K(0 < K <_H), there is a known temperature, measured during a time interval: e( x, t )| x=H =s(t ). (ll) Boundaries of the considered workpiece surface layer are defined by the known heat flux density: -X(0( x, t ) = q (t ), :=H and unknown heat flux density: -i{e( x, t = q(t ). (l2) (l3) £=0 The final solution of the inverse problem is the heat flux density on wheel/workpiece interface q(t) (part of total heat flux density [3]), and the temperature field 6 = d(x,t) throughout the elementary parts of the workpiece, D = {(x,t): xe[0H], t£[0,tj}. Due to their complexity, differential equations, which describe the process of heat conduction in grinding, are mostly solved using numerical methods. The first step with every numerical method is the discretization of space, i.e. the approximation of a thin, isolated plate of the workpiece by a number of elementary pieces Ax, Fig. 3. When dealing with a non-stationary -heat conduction problem, the time of temperature change and heat flux is discretized by a At increment. Fig. 3. The shape of a uni-dimensional mesh of heat conduction inverse problem To solve the partial differential equation shown in Eq. (9) an implicit form of the finite differences method was chosen [8]. A system of linear algebraic equations is used to calculate the unknown heat flux qn+1 on workpiece surface and the temperature field of the workpiece surface layer 6h«+i (h = 0, 1, ..., K-1, K+1, ..., H) as follows: [R]{0} = {B}. (l4) Solving the matrix system shown in Eq. (14) requires the initial task to be divided into two parts. First, a standalone system is calculated: [R2]{02> = {B2}. (15) The matrix Eq. (15) is solved as the direct task of heat conduction within the area D2 = {(x,t): xe[K,H], te[0,tm]}. The solution yields the unknown temperatures 6hn+1 (K+1, ..., H). Once vector 02 is determined, the problem of inverse heat conduction in the area of D1 = {(x,t): xe[0,K], te[0,tm]} can be tackled by using the system: [RjH©!} = {BJ. (16) From the system shown in Eq. (16), starting from the first unknown temperature , vector 0j is calculated. This vector represents the unknown heat flux q«*1 and the unknown temperature 0hn+1 (h = 0, ..., K-1). 2.2 Extreme Case of the Grinding Inverse Problem One of the possible ways to solve the inverse problem is to transform it into an extreme case using some method of optimization. Given the input parameters, this would allow us to determine the allowed thermal state of the process for the analytical model in hand, so as to satisfy the state and boundary functions for the given optimization criterion. For grinding, the task of optimal heat loading control is to determine the function of the process state, 0 = 0(x,t), and the control function q = q(t), so as to satisfy the analytical form of the inverse problem for the grinding process, as defined by Eqs. (9) to (13). These functions are determined under the condition that the known temperature at a particular fixed point 0(t) = d(K, t) is consistent with the temperature d (q, K, t) which is calculated based on control function q(t): _ -|2 J(q) = j|_0(q,K,t)-0(K,t)J dt. (17) In order to minimize the differences between the known and the calculated temperatures, an iterative method of optimization is used. J(q) is calculated by applying the iterative optimization algorithm for particular values of functional q^t). The procedure is repeated, decreasing the value of the functional J(qi+j) < J(qi), until its minimum is reached. The iterative gradient method procedure is considered finished for a sufficiently small functional, which means that the calculated vector components and the measured temperatures are very close. 3 EXPERIMENTATION 3.1 Experimental Setup The experimental work was carried out on a surface creep-feed grinding machine »Majevica« type CF 412 CNC. Fig. 4. Measuring system of grinding process 0 The workpiece material was high-speed steel (DIN S 2-10-1-8 or B.S. BM 42) at 66 HRc hardness, with dimensions of the parts 40*20*16 mm. The test used an aluminum oxide wheel »Winterthur« type 53 A80 F15V PMF, size 400*50*127 mm. The depth of cut was a = 0.05; 0.1; 0,25; 0.5 and 1 mm, the workpiece speed was vw = 2.5; 5; 10; 25 and 50 mm/s and the wheel speed was vs = 30 m/s. A water-based coolant (emulsion 6%) was used during the grinding test with a flow rate of 175 l/min. The temperature was measured in the workpiece surface layer using a thermocouple (type K, 00.2 mm) built into the workpiece at a specified clearance from the wheel/workpiece interface area. Measuring, analysis and control of the temperatures and forces during the process of machining with grinding was performed with the help of a computerized measuring system, Fig. 4. Metallographic identification of the state of the workpiece material surface layer after grinding was performed with an optic microscope. 3.2 Experiment Results To determine the temperature field distribution in the workpiece surface-layer, temperatures in the workpiece were measured for various distances from the measuring point to the contact surface of the workpiece and the grinding wheel. A characteristic graphic review of how creep-feed grinding temperatures change in time in the cutting zone, obtained through gradually drawing the grindstone closer to thermocouple's hot junction, is shown in Fig. 5. Metallographic examination of changes in the microstructure of the workpiece surface layer show that secondary changes appeared in each case when the measured maximal grinding contact temperature was higher than the temperature of the previous remission, which is 6c = 550 °C for the steel that was used in this testing, Fig. 6. The measured microhardness indicates that the microhardness of the secondary hardened layer is a little higher than the microhardness of the basic material. The smallest measured value of microhardness is found in the secondary remission layer. During further inspection of the workpiece material surface layer condition, not one sample was found with microcracks in it. On the other hand, burned surfaces were noticed in all samples where the measured grinding contact temperature was above the temperature of the previous remission. / 0