Application of model order reduction to a finite element model of cryogenic turning

Abstract

Metastable austenitic steels offer the opportunity of a surface layer hardening integrated in the machining process. The hardening effect is achieved by a deformation induced austenite-martensite phase transformation, for which high mechanical loads and low temperatures are necessary, typically below room temperature. These conditions can be accomplished during cryogenic turning, allowing a phase transformation in the surface layer of the workpiece. To study the austenite-martensite transformation behavior, information about the temperatures in the contact zone between tool and workpiece during machining is necessary, which can hardly be measured. Therefore, FE simulations of the process are utilized to calculate the transient process temperatures. In this paper the idea of projection-based model order reduction and its implementation into the finite element program FEAP are described, allowing a solution speed-up while maintaining accuracy of the full-order model. Motivation and ModelIn [1, 2] a model has been presented to calculate the temperatures arising during dry and cryogenic turning of metastable austenitic steel AISI 347. Pre-calculations were necessary to determine the magnitude of boundary heat fluxes first, which are mainly induced by the tool and the cryogenic cooling. This happens in an iterative way in terms of multiple forward solutions in order to correlate measured and simulated temperatures. However, the repeated calculation of such a finite element problem may take a considerable amount of time due to the number of equations to be solved. To counteract this drawback, model order reduction techniques may be considered, which reduce the number of equations and thus solution time, while maintaining accuracy of the solution. In the present case projection-based model order reduction, which is based on the method of snapshots [5], has been implemented into the finite element code of FEAP [6].The underlying equation is the heat equation, which in the present case considers the convective time derivative as well, resulting in an Eulerian representation, c.f. [2]. The finite element system of equations can then be formulated by K T + MṪ = F, which may be referred to as the full-order model (FOM) of size K, M ∈ R m×m and T, F ∈ R m .Model order reduction states the existence of an approximation T ≈ U rT , where the projection matrix U r is calculated from snapshots of FOM solutions. In more detail, n snapshots are gathered in a snapshot matrix ∈ R m×n , which is analyzed using singular value decomposition. The sought projection matrix stems from the left singular vectors U r = U[1 : d], d ∈ [1, n]. Additionally, further projection of the FOM is needed to overcome an over-determined system of equations. In case of Galerkin projection, the projection matrix U r is used for this second projection, like it was utilized in the present approach. In case of Petrov-Galerkin projection, a matrix different from U r performs the second projection. Here the system of equations after Galerkin projection with U r re...