Modelling of Grain Size Evolution with Different Approaches via FEM When Hard Machining of AISI 4140

: Thermo-mechanical loads during hard turning lead to the formation of so-called White Layers on the machined surface. Characterized by a very ﬁne microstructure and high hardness, White Layers have a negative effect on the fatigue life of a component. The fundamental mechanism for the White Layer formation is the dynamic recrystallization (DRX). Therefore, in the current work, two different DRX models, Helmholtz free energy and Zener-Hollomon, are implemented into Abaqus/Explicit to predict the thickness of the White Layer when hard turning quenched/tempered AISI 4140 and the results are compared with each other. For the simulation of the machining process a Finite Element Method (FEM) model based on the Coupled-Eulerian-Lagrangian (CEL) method is built up. Although both DRX models achieved a very good match between predicted and measured White Layer thickness and grain size evolution on the workpiece rim zone, the Zener-Hollomon model produced more closer agreement.


Introduction
Cr/Mo-alloyed steels are mainly used in the automotive and aviation constructions due to their high strength combined with high toughness. For example, AISI 4140 steel (42CrMo4 in Europe) in quenched/tempered condition has a hardness of approximately 48 HRC in addition to their high strength and toughness, which make this type of material challenging to machine. At this point, hard turning can be considered as a possible material removal process [1] and provides many advantages, among them are direct machining of the workpiece in hardened state, higher flexibility and cutting performance and requiring no grinding [2,3]. High temperatures, strains and strain rates that occur during hard machining, on the other hand, can be regarded as disadvantages which make this process highly interesting for research on Surface Integrity (i.e., general term for microstructural alternations, residual stress states, hardness, chemical decomposition and corrosion resistance on the workpiece rim zone) [4,5].
According to the detailed literature review by Jawahir et al. [5] and the recent research of Borchers et al. [6] about Surface Integrity in machining, the functionality and fatigue life of a machined part is strongly related to its topographical, metallurgical and mechanical properties in the rim zone. Hard turning is known for the induced thermo-mechanical loads that lead to altered microstructures, hardness, residual stresses, dislocation densities and grain sizes at the surface layers.
Typical thermo-mechanical related microstructural modifications are the White and Dark Layers [7,8]. Specifically, White Layer is a layer which has a few microns of thickness. It is characterized by a hard and brittle structure [8][9][10]. In contrast, Dark Layer is identified as a soft and ductile structure with a higher thickness than White Layer (20-40 microns) [10]. The combination of White and Dark Layers can significantly reduce the fatigue life, due to the inhomogeneous stress distributions at cyclic loads [2,3,11]. White Layers pose a high probability of formation and propagation of crack [12,13], while dark layers influence the magnitude and location of the residual stresses [10]. Griffiths summarized the three main theories for the White Layer formation; dynamic recrystallization (DRX) (severe plastic deformation), phase transformation (rapid heating and quenching) and surface reaction with the environment [7]. Dark Layer, however, is formed due to the microstructural changes triggered in the heat affected zone during the hard turning [14].
Most recent study by Fang-yuan et al. investigated the microstructure, phase transformation and carbon distribution of White and Dark Layers when machining AISI 52100 steel to understand their formation mechanisms [15]. They found that White Layer formation is related to rapid austenite transformation and quenching process, while Dark Layer formation is caused by a tempering process. They also showed that the carbon and chromium concentrations are different in White (the highest) and Dark Layers (the lowest).
In another recent study, Brown et al. focused on the characterization of super chrome molybdenum vanadium steel in terms of Surface Integrity (White Layer formation) [16]. They showed that in their case White Layer formation is triggered by severe plastic deformation, not by phase transformation, as demonstrated by the missing retained austenite peaks in the X-ray diffraction patterns and compressive residual stress profiles.
For many years, several researchers have been using different approaches to investigate machining-induced surface modifications, microstructural changes and grain size development. The most common ones are based on the following (DRX) models: Ramesh and Melkote presented a Finite Element approach in Abaqus/Explicit for the modelling of White Layer formation based on the martensite phase transformation for orthogonal cutting of AISI 5200 in hardened state [8]. They showed that White Layer leads to a compressive residual stress on the rim zone.
In another study, Akcan et al. examined the microstructure and mechanical properties of the White Layer on machining of AISI 52100 steel cylinders which were hardened and tempered [9]. Analyzes of the machined surface showed a White Layer thickness of 8.2 µm at 50 m/min cutting speed, 0.2 mm depth of cut and 0.1 mm/rev of feeding rate. An increase in the thickness of the White Layer was observed with the increasing cutting speed and flank wear.
For the hard machining of AISI 52100, Umbrello et al. [22] proposed a hardness-based flow rule, which is extended in [10] by advanced empirical models for phase transformation to simulate microstructural changes on Surface Integrity, such as White and Dark Layers with related changes in hardness. By considering the phase transformation, it was tacitly assumed that the austenitizing temperature would be exceeded during the process, which makes it difficult to implement this model for other materials. Caruso et al. applied the mentioned hardness-based flow stress and coupled it with the Zener-Hollomon model to predict the grain size [23]. They showed that the increase in cutting speed causes a thicker White Layer, while the thickness of the Dark Layer decreases with increasing cutting speed. In addition, a higher cutting speed leads to a Dark Layer dominated by grain refining. Caruso et al. also pointed out that the average grain size in the White Layer increases with increasing cutting speed [23]. The FEM software SFTC Deform ® -2D with Updated Lagrangian method and re-meshing were utilized in their simulations.
In a separate study, Ambrosy et al. implemented the Zener-Hollomon approach in Abaqus/Implicit with user-defined continuous re-meshing algorithm to model nanocrystalline grains in orthogonal cutting of AISI 4140 steel [24]. They revealed that the grain refinement and thickness of DRX-layer increase with the increasing relative roundness r β /h, which is the ratio of cutting edge radius r β to undeformed chip thickness h.
To simulate the effect of tool geometry and cutting conditions on microstructural alterations in turning of Ti-6Al-4V alloys, Arısoy and Özel used the JMAK model in SFTC Deform ® -3D [25]. They showed that regarding to change of microhardness on the rim zone tool coating does not play a decisive role at low cutting speeds (v c = 55 m/min), while the TiAlN coated tool leads for the investigated process parameter to larger grains in contrast to uncoated tool. Furthermore, in their study, the grain size decreases with smaller cutting radius at low cutting speeds. Similarly, the average grain size decreases with increasing feed rate.
The study by Estrin et al. introduced a dislocation density-based model to describe strain hardening behavior [19], which is then employed in [26] to model grain refinement in equal channel angular pressing for copper. Ding et al. applied this dislocation density-based model for the simulation of dislocation density and grain size evolution of titanium in [27] and of Al6061-T6 alloys in [28] after orthogonal cutting, respectively. This approach was later utilized by Atmani et al. [29] to simulate the grain refinement in machining of oxygen-free high conductivity copper. This DRX model has also been adapted for (hard) machining of steels latest by Li et al. for H13 steel [30] and earlier by Ding and Shin for AISI 52100 steel [31]. In these studies [26][27][28][29][30], the FE software Abaqus/Explicit with CEL along with ALE methods and in [31] AdvantEdge FEM were used.
Recently, Buchkremer and Klocke have presented a thermodynamically motivated DRX model based on the variation of Helmholtz free energy to predict the dynamically recrystallized grain size during the orthogonal cutting of AISI 4140 with SFTC Deform ® -2D [21]. Their experimental and simulative results showed that the austenitizing temperature is not reached or exceeded for cutting speeds v c = 50-150 m/min and undeformed chip thicknesses h = 0.05-0.2 mm. Thus, a phase transformation is not expected for the White Layer formation, which is also proven by their X-ray diffraction. Moreover, the results of the finite element model showed a White Layer thickness of 2 µm for v c = 50 m/min and h = 0.05 mm, 2.5 µm for v c = 150 m/min and h = 0.05 mm, 1.8 µm for v c = 50 m/min and h = 0.20 mm, 2.2 µm for v c = 150 m/min and h = 0.20 mm, which corresponds well to the presented experimental results. It can be deduced from these findings that the White Layer thickness increases with increasing cutting speed, while it decreases with increasing undeformed chip thickness.
Zhang and Zhuang analyzed the effect of cutting edge geometry on White Layer formation when cutting AISI 52100 steel [32]. Under the assumption that White Layer formation is related to rapid heating during the process, they implemented a phase transition temperature depending on stress, strain and current temperature to estimate the White Layer thickness with a hybrid approach. The 2D chip formation model was built up in Abaqus/Explicit using ALE method. They showed that White Layer thickness and process temperature increase with increasing cutting edge radius at constant cutting speed and undeformed chip thickness. This was also observed in experiments.
It can be concluded that machining-induced surface layer modifications are functions of thermo-mechanical process parameters and can be controlled or manipulated by skillful combinations. This paper focuses on FEM based modelling of White Layers in the hard machining of AISI 4140 steel, which can be used to determine the process parameters where the formation of White Layers is reduced or does not occur at all. The modelling results could be highly interesting from an industrial perspective as well. Since the FEM analyses are relatively time-consuming and cannot be used for monitoring a real process with respect to Surface Integrity, the accurate simulations can be utilized to calibrate a so-called Soft-Sensor, which is able to work with fast analytical empirical relations between the process forces and temperatures [33][34][35], in order to predict or avoid the possible formation of White Layers during the machining process. The calibrated Soft-Sensor could then be used to control the process and adjust related process parameters to increase the quality of machined parts [35]. This concept has not been well researched yet, but offers great potential.
Based on the research conducted by Buchkremer and Klocke, dynamic recrystallization is considered as a fundamental mechanism of the White Layer formation [21]. The above-mentioned empirical DRX models, Yada, Zener-Hollomon and JMAK, were very similar in terms of strain-based formulation of the DRX onset criteria (i.e., current plastic strain ≥ critical plastic strain) and recrystallized grain size calculation. Still, there are some differences. The main difference is that Yada model did not consider the recrystallized volume fraction as a function of time, while Zener-Hollomon and JMAK models took the DRX kinetics based on Arrhenius and Avrami type equations into account, respectively. In this work, the authors chose the Zener-Hollomon model among others because of its relatively simple implementation and calibration. It is also powerful enough to model White Layer formation due to DRX. In the case of the dislocation density-based model there are no specific DRX criteria and grain size refinement is calculated based on the total dislocation density which is given by a rule of mixtures. This approach will be analyzed in the upcoming work. Therefore, two different DRX models, Zener-Hollomon (empirical) and Helmholtz free energy (physics-based), were investigated and implemented in a Finite Element model and the results were compared and validated with different set of process parameters.

Material
In this work, the investigated material was a low alloyed quenched/tempered AISI 4140 steel, which is used for the manufacturing of steering knuckles, connecting rods, crankshafts, gear shafts or gear wheels in automotive and aircraft industries. Its chemical composition, mechanical and physical properties are given in Table 1. Table 1. Chemical composition, mechanical and physical properties of quenched/tempered AISI 4140 [21].

Property Notation Value Unit
Chemical composition

Finite Element Method-Based Modelling of Hard Machining
Hard machining is a complicated process to model due to the high temperatures, strain rates and extensive plastic deformations. The main problem hereby is the distorted elements during chip formation. There are several techniques to handle this problem in FE simulations, some of them are listed below considering the FEM software Abaqus.
In modelling of DRX on the rim zone temperature is one of the relevant state variables. Since element deletion does not involve temperature transfer from the deleted element to the element under it, this method is not suitable for the DRX modelling. Not to mention that due to element deletion, the machined surface layer is lost. In the CEL and UL methods there is a continuous contact between the elements themselves and the cutting tool, which makes them suitable for the calculation of dynamically recrystallized grain size on the workpiece surface. Due to the fact that the implemented re-meshing algorithm in Abaqus is not good enough to model the machining process, this paper focuses on the CEL approach.

Modelling with CEL Method in Abaqus/Explicit
A 3D model with a thickness of one element (2 µm) based on the CEL approach for orthogonal cutting was developed in Abaqus/Explicit 2018, as shown in Figure 1. In this approach, the Euler space has to be specified in terms of material (workpiece) and void to ensure that the removal chip flows out into the void. For a better contact behavior, the tool (rigid) should be placed in the Euler space. Besides, a plate (rigid), made of the same material as the tool, was placed over the workpiece in Euler space, to prevent the rolling chip from touching the surface of the workpiece. This method is called Coupled-Euler-Lagrangian, because tool and workpiece were meshed according to Lagrangian and Euler approaches, respectively. The Euler mesh is fixed in space and the material moves through it, while the nodes of an element in Lagrangian mesh move with the material and a mesh deformation takes place. For meshing of the workpiece and void the element type EC3D8RT was used, which stands for thermally coupled linear brick elements with eight nodes and one Gauß point. The tool and plate were meshed by utilizing the element type C3D8T, which describes an eight-node trilinear displacement and temperature element with eight Gauß points. Tool geometry and cutting parameters used in the simulations and experiments are listed in Table 2.

Euler space
To represent the DRX on the surface layer and the resulting chip geometry in Euler space successfully, the mesh size must be fine enough. For this reason, the mesh size of the workpiece in this model was set to 0.9 µm, while the mesh size of the tool on the cutting radius was 4 µm. For the modelling of DRX in orthogonal cutting, a full coupled thermo-mechanical analysis step was employed.
The boundary conditions in Euler space were formulated in terms of velocity, where the material flow in z-direction was prohibited in the whole model, while the material flow in y-direction was allowed for chip formation. The tool and plate were fixed in terms of velocity and displacement boundary conditions, respectively. Only workpiece moves with a prescribed cutting velocity v c against the tool.
Friction has a great influence on the simulation results in terms of temperature. In this work, Coulomb friction coefficient was formulated as a function of cutting speed v c . Figure 2 shows the measured Coulomb friction coefficient at different cutting speeds [21] and their linear interpolation, which is implemented into Abaqus. Since the CEL simulations were time consuming, it was necessary to implement a mass scaling to increase the explicit time step so that the computing time could be reduced. The modelling related material parameters of workpiece and tool are depicted in Tables 1 and 3, respectively. Note that the utilized mass scaling factor of 10 4 was not taken into account when specifying the properties, such as density and heat capacity.

Modelling of Dynamic Recrystallization
In this section, the material model and two different approaches to modelling dynamic recrystallization are introduced briefly.

Material Model
The visco-plastic material behavior of AISI 4140 was described based on the phenomenological constitutive material model of Johnson-Cook (JC) [37]. The yield potentialσ Y is given as a function of temperature, strain rate and plastic strain in Equation (1). The influence of stress state (stress triaxiality & Lode angle) and damage are not considered in this formulation. Furthermore, the strengthening mechanism through dislocation damping at high strain rates (>500 1/s) is also not taken into account.
where A is the yield stress at reference temperature T 0 and strain rateε p 0 , B stands for the strain hardening coefficient, n denotes the exponent of strain hardening, C and m are the material constants that represent the coefficient of strain rate hardening and the thermal softening exponent, respectively. While T m refers the melting temperature, T represents the current temperature. The calibrated model parameters for the material AISI 4140 were taken from Buchkremer and Klocke [21] and listed in Table 4. VUHARD subroutine was employed for coding Johnson-Cook model in Abaqus/Explicit, because this subroutine allows to call the values of the state variables from previous t − 1 and current t time steps, which are needed later when implementing the DRX models. VUHARD subroutine requires the given partial derivatives of the yield potential in Equations (2)-(4) as output to calculate the hardening behavior. Note that, the temperature was calculated by Abaqus during the explicit dynamic temperature-displacement coupled analysis.

Zener-Hollomon Model
Zener-Hollomon model is an empirical model for the prediction of dynamically recrystallized grain size. A detailed review of this model is given in literature [18,24]. Briefly, the Zener-Hollomon parameter Z H gives the relation between the working temperature and the equivalent plastic strain rate, which can be calibrated by performing compression tests at different temperatures and strain rates.
whereε p denotes the equivalent plastic strain rate, Q in J stands for the activation energy, R in J/mol/K is the universal gas constant and T in K represents the current temperature. The dynamically recrystallized grain size d DRX in µm can be calculated according to [18,24] as follows, where κ i represent the material related constants. According to this model, a critical plastic strainε p crit depending on the Zener-Hollomon Z H parameter and grain size from previous time step d t−1 has to be determined, as shown in Equation (7) [18].
In this model, the driving force for the DRX initiation is the equivalent plastic strainε p , so dynamic recrystallization takes place when the current plastic strain reaches or exceeds the critical one. Note that, the same criterion has been adopted for the initiation of White Layer formation. The DRX kinetics model in incremental form represented in Equation (8) and based on the Avrami function, is used to determine the recrystallized volume fraction triggered by the plastic strain [18]. To update the volume fraction, Forward Euler (explicit) integration scheme can be used, which is given in Equation (9).
In the equation, X DRX and dX DRX describe the volume fraction and its increment, respectively. While t represents the current time step, t − 1 is the previous time step. J is defined as the rate of dynamic recrystallization and depends on the strainε 50% at 50% DRX.
Average grain size d t in µm is composed of recrystallized d DRX and non-recrystallized d t−1 grains, as shown in Equation (11).
Based on the Hall-Petch relation the change in hardness can be calculated using current grain size according to the following equation, where κ HP is a material parameter and HV 0 in HV is the initial hardness [38].
The material parameters were iteratively fitted using a Python script developed by Mr. Alexander Schowtjak at IUL/TU-Dortmund. The parameters from Ambrosy et al. [24] were used as initial guess, because these parameters showed a very good agreement for the investigated material AISI 4140 with d 0 = 10 µm in their work. Yet, it was not possible to minimize the objective function with those values as initial guess. Therefore, this model was additionally implemented into MATLAB to optimize the initial guess based on an element on the workpiece rim zone.
The calibrated parameters at the cutting condition v c = 150 m/min and h = 0.05 mm are summarized in Table 5. The Zener-Hollomon model was implemented into the commercial software Abaqus/Explicit utilizing the user-defined subroutine VUHARD which was programmed in FORTRAN. The corresponding algorithm is given in Algorithm 1.

Helmholtz Free Energy Model
In this section, the Helmholtz free energy model is introduced. A full derivation of this approach can be found in the research of Buchkremer and Klocke [21]. This framework for DRX is thermodynamically driven and based on the variation of Helmholtz free energy dψ.
where u and s represent the specific internal energy and specific entropy, respectively. T stands for the temperature. The variation in specific internal energy du can be calculated as sum of the variation in mechanical de M , thermal de T and interface energy de I .
The variation of mechanical energy de M is given by the integration of the yield potential (here, of JC) with respect to equivalent plastic strainε p , as presented in Equation (15).
The driving force here is the stored dislocation energy, which is equal to the undissipated part of the specific mechanical energy (= 0.1 · de M ) [21]. The onset of DRX is characterized by the following three conditions: 1. Critical temperature must be exceeded: Critical specific mechanical energy must be exceeded: 0.1 · de M ≥ e M crit 3. Variation of the Helmholtz free energy should be negative: dψ < 0 All these three criteria must be fulfilled for grain size evolution due to dynamic recrystallization according to current model. The physical meaning of the 3rd criterion is that grain refinement takes place spontaneously, when Helmholtz free energy decreases [21]. Since the DRX is the fundamental mechanism for the formation of White Layers, same three criteria were interpreted from above as White Layer initiation conditions.
The new grain size d t in mm after DRX can be calculated as follows: where e gb (T) denotes the grain boundary energy and λ(T) represents the part of energy, which is transformed from stored dislocation energy into interface energy [21]. Grain refinement due to DRX has also an influence on the current plastic strainε p t which can be considered according to Equation (17) based on the Johnson-Cook material model (see Equation (1)) [21].
The calibrated model parameters are given in Table 6. Table 6. Helmholtz free energy model parameters of AISI 4140 based on [21].

Property Notation Value Unit
Grain boundary energy e gb (T) − 4.72 × 10 −7 · T + 4.098 × 10 −3 mJ/mm 2 Fraction of interface energy Critical specific mechanical energy e M crit 582 MPa Critical temperature T crit 887.55 K Initial grain size d 0 680 × 10 −6 mm The Helmholtz free energy model was implemented using FORTRAN in user-defined subroutine VUHARD in Abaqus/Explicit according to the Algorithm 2.

Results and Discussion
In this section, first, the validation of the introduced Finite Element orthogonal cutting model with respect to experimentally measured process forces at different cutting conditions is presented. Next, the predictions of the grain size refinement on the workpiece rim zone with Zener-Hollomon and Helmholtz free energy models are shown and finally, these results are compared and validated through experiments.

Experimental Validation of the Finite Element Orthogonal Cutting Model
To validate the FEM model, orthogonal cutting experiments at different cutting speeds v c = 50, 100 and 150 m/min and undeformed chip thicknesses h = 0.05, 0.1 and 0.2 mm with a constant cutting tool geometry (α = 3 • , γ = 6 • and r β = 66 µm) were performed. Figure 3a summarizes all predicted and experimentally measured and averaged cutting F c and feed forces F f , which were taken from a steady state condition and normalized by the thickness δ of the workpiece (δ sample = 3.3 mm and δ CEL = 2 µm). Figure 3b,c exemplary show the evolution of the force components from experiment and simulation with a cutting speed of v c = 50 m/min and an undeformed chip thickness of h = 0.2 mm. It can be seen from Figure 3a, on the one hand, the feed forces F f were generally underestimated with the CEL approach besides the cutting condition with v c = 150 m/min and h = 0.2 mm. On the other hand, the cutting forces F c which were depending on the cutting condition either slightly overestimated or underestimated showed a better correlation with the experiments (except v c = 150 m/min and h = 0.2 mm). The maximum deviations in cutting and feed forces were observed at a cutting depth of h = 0.2 mm and were (+)17 % and (−)41 % respectively. The general difference between CEL method and experiment in case of the feed forces has been also noted in literature [40,41]. As a possible explanation for this finding, the rapid deformation of the cutting edge in experiments based on the research [42,43] was given. Another reason could be the modeled length of the workpiece, which has an influence on the stress state and thus on the process forces in the CEL formulation.

Prediction of the White Layer Thickness Using Helmholtz Free Energy and Zener-Hollomon Models
For the validation of the Zener-Hollomon and Helmholtz free energy models the case studies from Table 7 were performed. Table 7. Hard turning experimental conditions for DRX-model verification [39].  It can be seen, on the one hand, that at lower cutting speed v c = 50 m/min there was a distinct difference between the two DRX models regarding to the calculated recrystallized grain volumes at the chip surface (Figure 4a,b) but the prediction of DRX on the chip surface is not an interest of this study. On the other hand, both models predicted at lower cutting speed independent of the undeformed chip thickness a stagnation zone in front of the cutting edge radius where no DRX occurred (Figures 4a,b and 5a,b). The stagnation zone is characterized by a lower relative velocity and also observed by Buchkremer and Klocke in their simulations [39,44]. Moreover, Figure 4 shows that at lower cutting depth h = 0.05 mm the Zener-Hollomon model could not predict the grain size evolution as accurate as Helmholtz free energy model and underestimates the grain size which can be recognized by the displayed color transitions blue(fine grains)-green-red(initial grain size) on the machined surface. In summary, as expected, both DRX models achieved similar results.  As can be observed from Figures 4 and 5 different cutting conditions cause different chip geometries. Therefore, an additional experimental validation of the used CEL model regarding to chip thickness ratio (=chip thickness/undeformed chip thickness) is provided in Figure 6. From the diagram in Figure 6 it can be concluded that the presented Finite Element model was able to represent the chip thickness at higher cutting speed better than at lower cutting speed. In literature Ducobu et al., who compared the CEL and ALE methods when machining Ti6Al4V, also reported that the CEL method overestimates the chip thickness [45]. Figure 7 summaries the measured and with two different DRX models simulated thickness of the White Layer at various cutting conditions (see Table 7). Both Zener-Hollomon (ZH) and Helmholtz free energy (HE) models while slightly overestimated the White Layer thickness in Case 1 and 2, a perfect agreement was achieved in Case 3. In Case 4, however, a more pronounced deviation was observed between simulations and experiment. A reason for this result is the high cutting forces that were obtained in the simulation under the condition v c = 150 m/min and h = 0.2 mm (see Figure 3a).  High cutting forces lead to higher temperatures and plastic strains on the machined surface and sublayer. This causes thicker White Layer due to the implemented DRX criteriaε p ≥ε p crit (ZH model) and T ≥ T crit (HE model). Figure 7 indicated that the Zener-Hollomon model predicted the White Layer thickness slightly accurate than the Helmholtz free energy model. This can be explained by the many material constants which has to be fitted in case of ZH model in a time-consuming calibration procedure. Note that the whole depth of the dynamically recrystallized layer from the simulations where the grain size was smaller than the initial one was interpreted as a White Layer. Figure 8 demonstrates the comparison of the predicted grain size evolution via Zener-Hollomon and Helmholtz free energy models on the workpiece rim zone after cutting in steady state with the experimental data. The path plots were created 0.4 mm behind the cutting tool after a step time of 5 ms. It can be seen that both DRX models performed very similar results which were in good agreement with the experiments. Furthermore, Figure 8 illustrates that the Helmholtz free energy model delivered better correlations at lower undeformed chip thickness (h = 0.05 mm), while Zener-Hollomon model showed better agreement at high cutting depth (h = 0.2 mm). Note that grain growth was not considered in the simulations due to rapid cooling rates of the machined surfaces.   Figure 8. Comparison of the predicted and experimentally measured grain size evolution with Zener-Hollomon and Helmholtz energy models at different cutting speeds and undeformed chip thicknesses.

Conclusions
In this study, Helmholtz free energy and Zener-Hollomon models were implemented into Abaqus/Explicit 2018 to predict the workpiece rim zone modifications based on the Surface Integrity in form of White Layer formation due to dynamic recrystallization (DRX) during hard turning of quenched/tempered AISI 4140. The aim was to compare the results from both models to determine, which one is more suitable for the investigated material. In addition, a 3D model for orthogonal cutting was developed using the Coupled-Eulerian-Lagrangian approach. Both DRX models and Finite Element model were validated experimentally.
To conclude, although Zener-Hollomon model tended to show a slightingly closer match, both Helmholtz free energy and Zener-Hollomon models achieved a very good agreement regarding to the predicted and experimentally measured White Layer thickness as well as the grain size evolution on the machined surface.
In further investigations, an advanced material model with damage coupling, more DRX cases and an additional DRX model based on dislocation density will be analyzed.