07881 - Artigo IPDO 2013

07881 - Artigo IPDO 2013

(Parte 1 de 2)

4th Inverse Problems, Design and Optimization Symposium (IPDO-2013)

Rogério Fernandes Britoa , Carlos Adriano Corrêa Ribeirob

, Solidônio Rodrigues de

Carvalhoc , Sandro Metrevelle Marcondes de Lima e Silva a Universidade Federal de Itajubá – UNIFEI, Rua Um, s/n, Distrito Industrial I, CEP 35903-081,

Itabira, MG, Brasil, rogbrito@unifei.edu.br b Universidade Federal de Itajubá - UNIFEI, Instituto de Engenharia Mecânica – IEM, Laboratório de

Transferência de Calor – LabTC, Av. BPS, 1303, Bairro Pinheirinho, CEP 37500-903, Caixa Postal 50,

Itajubá, MG, BrasilI, cadrianocr@unifei.edu.br metrevel@unifei.edu.br c Universidade Federal de Uberlândia – UFU, Faculdade de Engenharia Mecânica – FEMEC, Av. João

Naves de Ávila, 2160, Bairro Santa Mônica, CEP 38400-902, Uberlândia, MG, Brasil, srcarvalho@mecanica.ufu.br ;


This work proposes the use of inverse problem techniques in connection with COMSOL, with the aim to estimate the heat flux and the temperature field on a turning cutting tool in transient regime. Specification function, which is an inverse problem technique, was implemented in a program to estimate the heat flux applied on the tool, from the experimental temperature records. Once the heat flux is known, COMSOL is again utilized to obtain the temperature field on the cutting tool. The validation of the methodology is carried out by comparing the numerical results of the temperature with the experimental ones.

Nomenclature F objective function h heat transfer coefficient, W/m2 K

Mgeneral time index

k thermal conductivity, W/mK unknown heat flux, W/m 2 r number of future time steps T numerical temperature, °C t time, s

T∞ the medium temperature, °C

T0 initial temperature, °C Y experimental temperatures, °C x Cartesian coordinate, m y Cartesian coordinate, m z Cartesian coordinate, m

Greek α thermal diffusivity, m/s 2 η the outward drawn normal to the surface

Subscripts p index of points

Superscript np number of points ns number of temperature sensors

4th Inverse Problems, Design and Optimization Symposium (IPDO-2013)


Many engineering processes have their efficiency affected by high temperature values. An example is machining process in which the lifespan of the cutting tool is highly affected by high temperatures. The right knowledge about these temperature values and heat flux applied in this process enables the development of more efficient cooling techniques. However, the direct measurement of the temperature in a machining process is difficult to accomplish due to the movement of the piece as well as the presence of chips. Thus, the use of inverse heat conduction techniques conveys a good alternative to obtain these temperatures, since these techniques allow the use of experimental data obtained from accessible regions. Inverse problems consist of obtaining the value of a variable through the measurement of another variable measured directly. These techniques often use optimization algorithm in order to minimize the error between the calculated and real value of the variable in question. Nowadays, several researchers proposed the combination of inverse techniques and analytical or numerical heat transfer solutions to analyze the thermal fields during machining processes. One way to study the heat transfer problem is to know the heat flux to determine the temperature fields from the solution of the heat diffusion equation. In literature, this methodology is called direct problem in heat transfer. However, during machining the experimental heat flux is unknown and inverse techniques have to be used to predict this parameter. This proposal estimates the transient heat flux from a numerical model based on heat diffusion equation by using inverse problems and experimental temperatures measured at accessible regions of the sample.

Several engineering processes have their performance and quality affected by high temperature values. A typical example is the machining process in which cutting tool temperatures may be higher than 900ºC [1]. High temperatures cause changes in the microstructure of the tool during machining, change their physical chemical properties reducing their capacity to mechanical stress that appear during their use [2]. The direct consequence of these alterations is the reduction of their lifespan and performance. This leads to high operation costs and reduction of the quality of the end product. The right knowledge of the temperature values and applied heat flux in this kind of process, results in advantages like the development of more efficient cooling techniques as well as better specifications of the cutting parameters in machining processes. According to [3], [4] were the first to study the applied heat flux in cutting tools, as those used in turning. Since then, several other researchers have studied the subject. These studies involve experimental techniques like the use of thermocouples inserted in micro holes in the tool, [5], and experimental infrared chambers, [6]. Other ways to measure the temperature in a cutting tool are mentioned by [4], like the optical pyrometer, tool–workpiece thermocouple, among others. Along the years, the instrumentation techniques have also evolved. [7] have also used thermocouples in their studies on the temperature of the cutting tool; however these thermocouples were fixed by capacitive discharge. New optical pyrometer and infrared chambers presenting better precision and shorter response time have been used nowadays, as in [8]. However, in some situations, analytical methods have been used to predict temperature fields on a cutting tool, as in [9]. Although the means of instrumentation have also evolved, the direct measuring of the temperature in a machining process is difficult due to the relative movement of the workpiece and the presence of chips. Hence, the right instrumentation must be used to measure the temperature fields on the tool. Even though there is a variety of instruments that may be used to measure temperature fields on a cutting tool, there must be criteria to select them, for each instrument convey advantages and disadvantages when applied in a determined situation. The thermocouples, for example, may be reliable, however the presence of chips in the cutting operation may damage them, and the use of micro holes to insert them change the geometry of the tool, leading to mistaken results. The chambers and pyrometers are not subject to damages by the presence of chips and changes in the geometry of the tool, but can not capture the temperature on the exact region in contact with the cutting tool and the workpiece.

Conveying a greater availability of computational resources, the use of numerical methods gained terrain, and it did not take long before they started being used, along with experimental methods in the studies of temperature fields on cutting tools. The most common are, respectively, the finite difference method and the finite element method. Later, the finite element method also emerged. Each of these methods has its own characteristics. The finite difference method had its spread in fluid dynamics problems due to its ability to deal with non-linear terms while the finite element method is the typical method to solve purely diffusive problems, [10]. The finite element method presents a characteristic of

4th Inverse Problems, Design and Optimization Symposium (IPDO-2013) being conservative in terms of elements, that is, it balances energy, momentum and mass in all the elements of the mesh, which does not occur with the aforementioned methods. More than one numerical method may be used in a certain application, aiming to explore the strengths of each numerical method, as in [1].

Another consequence of the evolution of computers is the development and increasing use of inverse techniques. Inverse techniques consist of obtaining the value of a variable of interest through the measuring of another property which may be measured directly from a region of easy access. Therefore these techniques seem to be a good alternative for thermal analysis in machining processes once they allow, for example, the obtaining of the temperature fields and the applied heat flux on a cutting tool from the temperature signal obtained at a single point of easy access.

Inverse techniques have already been used to study temperature fields on a cutting tool, as in [6] and [12]. In such situations, the temperatures on the cutting tool are obtained by using thermocouples inserted close to the rake face as in [7], [13], [14], [15] and [16]. An alternative to the use of thermocouples is the infrared chamber, as in [15]. Due to the nature of the process, which precludes the direct measurement of the heat flux on the rake face, inverse technique such as Golden Section, Specification Function among others have been used, in several of these studies, to estimate heat flux applied on the tool.

This work proposes the use of inverse problem techniques with the commercial software COMSOL® 4.3, to estimate the heat flux and the temperature field, in a transient regime, in a turning cutting tool. A Fortran program, with implemented inverse technique, Specification Function, was developed to estimate the heat flux applied on the cutting tool, from experimental temperature records, in a determined point.

Once the heat flux is estimated, the software COMSOL® 4.3 is once again used to obtain the temperature field in the cutting tool. The validation of the methodology is done by comparing the numerical temperature results with the experimental data.


2.1 Thermal Model Set

The problem dealt with in this work is represented by Figs. 1a and 1b. The Figure represents a set consisting of a cutting tool, a hard metal, a wedge positioned under the cutting tool between the tool and the tool holder. There is also a staple and a bolt to fix the set. A perspective is shown in Figure 1a, whereas a blown up view of the set is shown in Fig. 1b. The heat diffusion equation ruling this problem may be given as:

T tzyx

T tzyx

T tzyx

Subject to the following boundary conditions

on the contact interface with the workpiece (Fig. 1b) (2)


k in the remaining regions of the set (3) and having the following as the initial condition

As the linear Function Specification method was used in this work, the thermal properties were consideed constant.

4th Inverse Problems, Design and Optimization Symposium (IPDO-2013)

a) b)

Fig. 1. a) Tool set, Shim and Tool holder, and b) Detail of the contact interface between the tool and the workpiece.

2.2 Direct Problem Solution.

The direct problem consists in solving the heat diffusion equation according to the boundary conditions mentioned above. The COMSOL® 4.3 program, which solves thermal problems by using the finite element method, is used for this purpose.

2.3 The Inverse Problem.

The inverse technique adopted in this work is the Specification Function [17]. This technique demands, first, the sensitivity coefficient calculation, which is done numerically from Duhamel Theorem [18]. The sensitivity coefficient is then obtained with the use of a numerical probe which follows the temperature change, in the point equivalent to that on which the thermocouple was placed in the experiments. Once the sensitivity coefficient is at hand, the head flux is estimated with the use of a Fortran language program. The methodology is represented in Fig. 2.

Another important parameter is the value of future time steps r. In the Specification Function technique, a determined value of future time steps r is used to estimate the heat flux at present instant. In the resolution of the inverse problem, the Specification Function searches for a heat flux value that minimizes the objective function given in Eq. 5, for each time step rp ns j pMjpMjTYF (5)

Contact area x z x z

4th Inverse Problems, Design and Optimization Symposium (IPDO-2013)

Fig. 2. Methodology for the use of the Specification Function for the tridimensional thermal model, in transient regime with the COMSOL® .


The validation of the methodology was carried out according to [18]. The experiments were accomplished under controlled condition, in which the heat flux and the temperatures were measured from a homogeneous AISI 304 stainless steel sample of 60.0 x 100.0 x 9.5 m. The experimental apparatus is represented in Fig. 3.

Fig. 3. Experimental Apparatus used to validate the proposed methodology.

Fig. 4. Position of the thermocouples on the tridimensional model.

During the experiment, the homogeneous AISI 304 stainless steel samples were heated uniformly on their surface by a resistive kapton heater, with a 15 resistence and a heating area of 50.0 x 50.0 m x y

Thermocouple 1 (5,0; 6,0; 0,0)

4th Inverse Problems, Design and Optimization Symposium (IPDO-2013) and 0.20 m thick. As the heater has a smaller area than the surface area of the sample, a tridimensional heat propagation effect is obtained. This tridimensional problem has boundary condition of heat flux imposed in the heating region and isolated in the remaining regions [17]. Figure 4 shows the position of the heater and the thermocouple on the metal sample.

A digital source Instrutemp ST-305D-I was used to generate the heat flux, through Joule effect. The control of the applied heat flux intensity is done through the variations of the voltage and current parameters. To minimize the errors in the measurement of the heat flux, a symmetrical assembly was used. Moreover, the applied voltage and current values were measured by the multimeters Instrutherm MD-380 and Minipa ET-2042C, previouly calibrated. The convection effects are disregarded due to the use of polystyrene involving the assembly, so as to keep the width of 50 m, in any direction. Furthermore, the whole set was placed inside a Marconi MA030 oven. The interstices between the samples are reduced with the use of Artic silver paste, thus improving the contact between the samples.

The solution of the tridimensional heat diffusion equation is obtained with the use of the finite element method, through the commercial software COMSOL® 4.3. For this, a computational thermal model was used to faithfully represent the experimental model of the sample. This model was discretized in a computational mesh of hexahedrical elements. In solving the discretized equations that rule the thermal model, a temperature field is obtained in any point of the sample. The mesh displays a refinement so as to have 60 elements in direction x, 100 elements in direction y and 20 elements in direction z, totalinrobablug 120.0 volumes.

(Parte 1 de 2)