Considering multiple process observables to determine material model parameters for FE-cutting simulations

Analyzing the chip formation process by means of the finite element method (FEM) is an established procedure to understand the cutting process. For a realistic simulation, different input models are required, among which the material model is crucial. To determine the underlying material model parameters, inverse methods have found an increasing acceptance within the last decade. The calculated model parameters exhibit good validity within the domain of investigation, but suffer from their non-uniqueness. To overcome the drawback of the non-uniqueness, the literature suggests either to enlarge the domain of experimental investigations or to use more process observables as validation parameters. This paper presents a novel approach merging both suggestions: a fully automatized procedure in conjunction with the use of multiple process observables is utilized to investigate the non-uniqueness of material model parameters for the domain of cutting simulations. The underlying approach is two-fold: Firstly, the accuracy of the evaluated process observables from FE simulations is enhanced by establishing an automatized routine. Secondly, the number of process observables that are considered in the inverse approach is increased. For this purpose, the cutting force, cutting normal force, chip temperature, chip thickness, and chip radius are taken into account. It was shown that multiple parameter sets of the material model can result in almost identical simulation results in terms of the simulated process observables and the local material loads.


Introduction
The use of simulation techniques, such as the finite element method (FEM), to model the machining process enables the prediction of difficult-to-measure process quantities, such as stresses, strains, and temperature fields [1,2]. The knowledge about these quantities plays a pivotal role to improve the process understanding and the identification of the acting mechanisms. Machining simulations offer the possibility to enhance the process design by overcoming classical, empirical "trial-and-error" approaches [3,4]. In metal cutting, the backbone of such a simulation is the chip formation process. To model the chip formation process by means of FE simulations, different parameters and input models are required. Among them, the material and friction models are essential for the reliability and accuracy of simulated results [5][6][7].
In the state of the art, the determination of material model parameters is conducted by two different approaches: & By direct determination of the model parameters from conventional tests (e.g., tensile tests) and nonconventional material tests (e.g., split Hopkinson pressure bar (SHPB) tests), or & By inverse techniques.
For direct tests, the thermo-mechanical loads are limited, wherefore the occurring loads are far away from those of the cutting process [4,8] and extrapolation into the regime of metal cutting becomes necessary [9][10][11][12]. This extrapolation into the regime of metal cutting conditions can cause inaccuracies of the simulated results [13]. To overcome this issue, Venuvinod and Jin suggested to use the cutting process itself as a material test [14], since a process can be represented best by itself in terms of the occurring loads. For this inverse parameter determination, the model parameters are iteratively modified until the simulations of the process are in agreement with the experimental measurements. Therefore, different process observables are assessed to evaluate the deviation between simulations and experiments, such as the process forces or the chip form.
The inverse parameter determination has found increased application within the last decade. Özel and Altan presented one of the first approaches for the inverse parameter determination from the cutting process [15]. To calibrate the material model parameters for AISI P20 mold steel, flow stress data from strain and strain rate tests were used as starting point to model the chip formation process, followed by an iterative fitting of the parameters. To evaluate the match between experimental and simulative data, only the cutting force was considered. However, Bäker reported that it is not sufficient to consider just the cutting force, since the cutting force is less sensitive to an accurate material model than other integral process observables [16]. Further studies on the inverse parameter determination for metal cutting simulations were presented by Shrot and Bäker [17][18][19][20]. The authors developed a procedure that bases on the Levenberg-Marquardt algorithm, which was utilized for the inverse determination of the Johnson-Cook (JC) material model parameters A and B. Their results showed that it was possible to re-identify the material model parameters within a relatively small number of iterations. However, it remains uncertain if the procedure is applicable for the inverse determination of all JC material model parameters. Klocke et al. presented a method for the inverse parameter determination of the JC parameters C and m [21,22]. Therefore, lower and upper values of the parameters were guessed and their results were interpolated in order to find the best fit with the experimental results. Later, Klocke et al. presented a pure iterative procedure for the inverse determination of material and friction parameters of Inconel 718 Direct Aged [23]. Albeit, their procedure is dependent on the experience of the user, wherefore the robustness remains questionable. Zhang et al. developed a novel procedure for the inverse determination by taking cutting forces and deformation fields into account, which were measured by means of the digital image correlation (DIC) [24]. Due to technical limitation of the utilized DIC system, only low cutting speeds of up to v c ≤ 0.1 mm/min were considered. Further, the effect of strain rate hardening and thermal softening was neglected, which are two important effects that influence the material behavior in the cutting process. Besides the outlined drawbacks of the inverse parameter determination, Arrazola et al. [3] and Shrot and Bäker [25] reported that an absolute solution of the parameter set in not necessarily given, because many combinations of parameters can result for the investigated domain in the same simulation results, which means that the material model parameters are not unique [2].
To reduce the computational effort and to increase the robustness of the solution, optimization algorithms were applied to the inverse problem. In the field of sheet metal forming, Chaparro et al. utilized a genetic and a gradient-based algorithm for the inverse determination of material model parameters [26]. The authors showed that it was possible to fit the numerical with experimental data by using optimization algorithms. For machining simulations, Özel and Karpat used the Particle Swarm Optimization to determine JC parameters [27]. As experimental data basis, the authors utilized SHPB tests, which required extrapolation into the domain of metal cutting.
To improve the inverse parameter determination, the authors of this paper developed a novel approach by using the Downhill simplex algorithm (DSA). In previous papers, this optimization algorithm was examined regarding its potentials for the inverse problem [28], its applicability for experimental results [29], and its influencing factors [30]. However, the results of these studies resulted in non-unique parameter sets for the domain of investigation. Further, the evaluation of the simulations required a manual comparison between the simulations and the experimentally derived chip form, which was characterized by high efforts and limited accuracies [30]. Different approaches for this inverse re-identification were used, all suggesting that the parameters sets are not unique. However, the simulated process observables deviated just slightly. For the developed procedure, two drawbacks were eminent: the effort to evaluate the simulations and the accuracy of the manually determined simulation results. For the iterative procedure, the effort to evaluate the chip formation simulations can easily sum up to multiple days. When analyzing the results from FE simulations, the determined results and their accuracy depend on simulation-inherent parameters, such as the mesh size as well as on the user, who analyzes the simulations. Some relevant simulation results that can be used for the comparison between simulations and experiment can only be extracted with high efforts, such as the chip radius.
The novelty of this paper is two-fold: in the first part, an approach for an automatized inverse determination of material model parameters from coupled Eulerian-Lagrangian (CEL) FE chip formation simulations is developed. Further, the number of process observables to assess the deviations between experiments and simulations is increased. To extract information on the chip form from CEL simulations, a method to reconstruct the chip contour was developed and validated.
In the second part, the automatized evaluation of chip formation simulations is employed in an optimization-based inverse parameter determination. The non-uniqueness is evaluated in terms of the considered process observables as well as the local thermo-mechanical material loads.
The paper is organized as follows: in the following chapter, the CEL model of orthogonal cutting is presented. Thereafter, the theoretical background of the two-fluid interface problem is outlined and applied to orthogonal cutting simulations. Further, a methodology to extract process observables from the FE simulations is developed. The automatized analysis is integrated in an inverse parameter determination procedure. Inverse re-identification approaches are conducted and the calculated parameter sets are evaluated. Thereafter, the results are discussed. At the end of the paper, a summary will be given and an outlook will be drawn.

Coupled Eulerian-Lagrangian model of orthogonal cutting
To model the chip formation process, a CEL model of orthogonal cutting was developed within the commercial simulation software Abaqus 6.14-6. Orthogonal cutting was chosen as cutting process since it can be considered as the most elementary process and allows to consider the process as two-dimensional. The CEL formulation was firstly applied by Klocke et al. [31] and Ducobu et al. [32] to the field of machining. In the CEL model, the tool was discretized by using the Lagrangian formulation. For the Lagrangian tool, 8-node thermally coupled elements of type C3D8RT were utilized. The workpiece material on the other side was modeled by the Eulerian formulation so that the material can freely flow through the mesh. Therefore, 8-node thermally coupled linear Eulerian elements with reduced integration and hourglass control (C3D8RT) were used. Additionally, an area for the formation of the chip and an area of the machined material were discretized by the Eulerian formulation. Thereby, the chip was not restricted in its form, but only in its length, which in turn was limited by the size of the modeled domain. Within the Eulerian and the Lagrangian domain, the minimal mesh size was 5 μm. For the simulations, the explicit formulation was utilized, which is favorable for non-linear dynamic problems such as the metal cutting process [3]. The setup of the CEL model that was used is illustrated in Fig. 1. It resembles a model used in a previous study by the authors [28].
An initial inflow of the material was used to model the cutting speed v c . The inflowing material exited the Euler domain either in form of the chip or as machined material. For all conducted simulations, a cutting length of l c = 3.33 mm was utilized. This cutting length showed to be sufficient to reach the steady state in the simulations. To reduce the computational time of the simulations, a mass scaling factor of 1000 was applied. The simulations were conducted on a virtual machine with 16 CPU of type Intel® Xeon® CPU E5-2698v4 with 2.2 GHz. For each simulation, 4 CPU were used. Therefore, a simulated cutting time of t c = 0.002 s required approximately 2 h and 20 min of computational time. Further details on the used material parameters, material properties, and boundary conditions are given in [28].
The JC material model was utilized to model the material behavior of AISI 1045 [33]. The JC model is one of the most widely used material models to describe the constitutive workpiece material behavior in metal cutting simulations, Eq. (1). AISI 1045 was chosen as material of investigation, since this material tends to form a continuous chip over a wide range of cutting conditions. Therefore, a damage model to simulate chip segmentation was not considered within this study.
In the JC model, the effects of strain hardening, strain rate hardening, and thermal softening are considered by the three separate terms in the brackets. The variables A, B, C, n, and m are the model parameters, which are assumed to be constant,

Tool (Lagrangeformulation)
Initial workpiece T m is the material's melting temperature, and the symbolsε 0 and T 0 denote reference values. In the present paper it is aimed to evaluate the developed methodology for the inverse parameter determination. Therefore, the routine is applied to the inverse re-identification of material model parameters. The material model parameters to be inversely re-identified were taken from the literature [22] since they are unknown a priori when using experimental data (Table 1). This allows an evaluation of the determined parameter sets in comparison to the target parameter set and, therefore, an assessment regarding the uniqueness of material model parameters. Since numerical errors, such as rounding errors, can occur in FE simulations, non-uniqueness is assumed in this paper if the evaluated results differ by less than 1%.

Two-fluid interface problem
The determination of the chip form from CEL simulations is a two-phase problem, where the interface of the two phases needs to be determined. In order to calculate the chip form, established methods from the field of fluid mechanics are transferred to the given task. In CEL cutting simulations, the two-phase flow is given by material that flows in the void of the Euler domain. The method aims to decrease the evaluation time and to increase the accuracy of the results. Simplified, the problem of two immiscible fluids in a mesh can be sketched in a fixed coordinate system, where some cells are fully filled with one fluid (Fig. 2a, blue region) and others with a different fluid (Fig. 2a, gray region). The cells filled with both fluids define the two-fluid interface. The volume-of-fluid (VOF) method is a widespread technique to analyze interfaces in two-phase flows [34]. Within this study, the concept of the VOF is applied to determine the chip form.
For the reconstruction of the two-fluid interface, volumetracking methods were applied, which differ by the features of the interface reconstruction algorithm and by the time integration of the volume fraction advection equation. Rider and Kothe [35] proposed a piecewise linear interface calculation (PLIC) method and a multi-dimensional unsplit time integration scheme, in which material volume fluxes are computed by using a set of simple geometrical tasks.
Within the CEL simulation in Abaqus, the concept of the Eulerian volume fraction (EVF) is used to describe the fraction volume of a cell occupied by a fluid. In this concept, the EVF is equal to one in the case that the cell is filled with the fluid and zero if the cell does not contain the fluid or a different fluid. Cells that contain a portion of the fluid in its volume have an EVF proportional to the share of the fluid and are called mixed cells. For these cells, the border between the two fluids can be approximated by polygons. This can be done with the piecewise linear interface calculation (PLIC) [35]. The algorithm bases on a second-order geometric solution of a volume evolution equation. The method uses local discrete material volume and velocity data to track the surface of a topology. In a single cell of a two-dimensional regime, the geometric function for the interface is assumed to be a straight line. For the application to CEL-chip formation simulations, the assumption of a straight line instead of a function of higher order appeared to be reasonable, as long as the mesh size is small. The assumption of the straight line will be evaluated in section 4.1.
Afterwards, the resulting polygon can be defined. The vertices of the polygon of a mixed cell (x v , y v ) are those lying inside the fluid (points 1 and 2) and the intersection points of the cell and the straight line (points 3 and 4) (Fig. 2b). The nsided polygon in a Cartesian x/y-geometry surrounds the fluid in each cell, whereby the matrix of the enclosed area A is calculated according to: Based on the given equations, the linear interface has to be reconstructed for every mixed cell of the mesh. The line constant ρ can be calculated by volume conservation, and the normal n is calculated from volume fraction gradients. The value of ρ is limited in such a way that the volume of the polygon is equal to the volume of the fluid within the cell. This relationship can be described in the 2D case by: where A c is the total cell area and A ρ the area of the polygon. An established method to calculate the normal vector for general meshes is the extension of Youngs' second method   [36] that was developed by Kothe et al. [37]. In the algorithm, the finite-difference approximations to determine the volume fraction gradient ∇f are utilized. This procedure is derived from the work of Barth [38], who developed the linear and quadratic reconstruction of discrete data on unstructured meshes based on least square algorithms. For this approach, volume fraction Taylor series expansions f TS i are formed from each reference cell volume fraction f i at x i to each neighbor cell The result is the volume fraction gradient ∇f i , Eq.
The normal vector is computed as: where 4 Reconstruction of the chip geometry from CEL-cutting simulations The EVF values and the position of the elements are extracted from the simulations to apply the theory of the two-fluid interface problem and to reconstruct the chip form. The centroid position of the elements was calculated in the region of the CEL model, where the chip forms and only those elements with 0 < EVF < 1 were considered. To determine the normal vector n of each mixed-cell element, the PLIC method in conjunction with Youngs' least square method was utilized. Thereby, the results from each cell were linked with the results from its neighbor cells.
The information on the EVF values, total element area, and normal vectors were used to iteratively solve Eq. (3). For the iterative procedure, the Wijgaarden-Dekker-Brent method [39], also known as Brent's method, was used. Brent's method combines parts of the bisection method, secant method, and inverse quadratic interpolation. The advantage of using Brent's method is that the algorithm converges for functions computable within a fixed interval [40]. Here, the interval is defined by the lower and upper boundaries ρ min , ρ max . These boundaries are estimated based on the position of the cell vertices and their normal direction (Fig. 3). As convergence criterion, a remaining error in Eq. (3) of 0.01 % was defined.
After the algorithm runs through all of the mixed cells, the geometry of the chip can be reconstructed (Fig. 4).

Determination of chip form parameters
To determine the chip form parameters, the simulated length must be limited in such a way that the chip does not hit the surface of the workpiece and curls up. This would affect the chip radius. The two chip form parameters that are evaluated within this study are the chip thickness h and the chip radius r c . The reconstruction provides the circumference of the chip with an outer radius r out at the tool side and an inner radius r in at the free chip surface. The calculation neglects a certain distance in the x-direction within the Euler domain to avoid the influence of the chip curl. Based on the calculation of the centroid coordinates, the value of the chip radius can be determined using the equation of a circle. The results of the calculated radii are averaged in the inner and outer chip radius r c . A least squares approach is employed to calculate the centroid of the circumferences. The generalized matrix to compute the coordinates of the chip centroid (O x , O y ) is given in Eq. (11), where P i are the determined points on the outer circumference of the chip.
Using the least squares approach leads to a numerical residue, which is also given to evaluate the results. The numerical Finally, the chip thickness h' is calculated as the difference between the outer and the inner chip radius, Eq. (12).

Validation and evaluation of the determined chip form parameters
To evaluate the automatized determination of the chip form parameters from FE chip formation simulations, the results from the algorithm were compared with the manual measured ones from the simulated chip form and to each other by the resulting residue. The chip radius was not considered, since it is impossible to measure it accurately by hand. When using the automated approach, the residue was in every case lower than R max < 5 × 10 −5 mm. This is smaller than the smallest mesh size within the Euler domain (l mesh,min = 0.005 mm). This suggests a reasonable good agreement with the underlying data points. It also confirms that linear elements are sufficient to reconstruct the borderline.
The comparison of the manually and the automatically determined chip thicknesses is shown in Fig. 5 for different cutting speeds and undeformed chip thicknesses. The results emphasize that the values of the automatic procedure and of their manually determined counterparts are in a close agreement. This can be attributed to the small mesh size enabling a high resolution of the chip contour. However, the standard deviations of the manual measurements were higher.

Extraction of process observables from FE-chip formation simulations
Besides the chip form, other process observables, such as the cutting force, cutting normal force, or thermal measurements, were used to validate cutting simulations. When evaluating the cutting force components from cutting simulations, the forces are usually averaged within the region of steady state.
To have single values representing the mechanical reaction forces in each coordinate direction and to establish a consistent procedure to evaluate the force signals, a data filter was applied to the automatized evaluation of FE chip formation simulations. The filter function bases on the temporal averages and its differences [41]. Therefore, the simulated signals of the cutting force and cutting normal force are firstly passed through a Wiener filter to reduce undesired noise and eliminate outliers that would cumber the definition of equilibrium [42]. Afterwards, a moving average with a step size of 5% of the total signal length is applied. The resulting list of averages is tested for the minimal discrepancy to its neighbors. If the resulting discrepancy is lower than the tolerated one, which was set to 5%, the function is considered to be in the steady state. Additionally, the end of the signal is cut off, since it can include undesired residues. An example of the force evolution is shown in Fig. 6.
As a further process observable, the chip temperature was taken into account in the inverse re-determination, since it was possible to measure the chip temperature experimentally [29]. For the proper evaluation of the chip temperature, the measuring position as well as the diameter of the measuring spot was taken into account (Fig. 7). The tool-sided chip temperature was averaged within a spot of 25 μm diameter. The developed procedure to analyze the results from FE chip formation simulations was utilized in an inverse determination of material model parameters. Therefore, an initial parameter set from the literature was used [22] (Table 1). With this parameter set, the chip formation process was simulated and the process observables were extracted from the simulations as target values for the inverse parameter determination.
To evaluate the deviation between the iterative simulation results x sim and the target values x target of the considered process observables, the observables are weighted in a single objective function: where ω i are weighting factors that weight the individual process observables. Their values were set based on established empirical values to ω F c ¼ 0:30, ω F cN ¼ 0:15, ω h 0 ¼ 0:20, ω r c ¼ 0:20, and ω T = 0.15. The optimization was performed with the Downhill simplex algorithm (DSA) [28][29][30]. The DSA, also known as Nelder-Mead simplex algorithm [43], is an unconstrained direct search method for multi-dimensional optimization [44]. For the optimization in an n-dimensional search space, the algorithm compares n + 1 vertices and replaces the worst vertex in terms of the evaluation function ξ. Five different operators are applied within the algorithm that modify the simplex spanned by the vertices. The operators are reflection, expansion, internal and external contraction, and shrinkage. The parameters of the operators that influence the step size of the algorithm were set to ρ DSA = 1 (reflection), ϵ DSA = 0.5 (expansion), κ DSA = 0.5 (internal and external contraction), and σ DSA = 0.5 (shrinkage) [28,30].
The initial simplex for the conduction of the algorithm was randomly generated and the boundaries of the to-bedetermined JC parameters were limited to physical meaningful ranges ( Table 2). The same initial simplexes were used as in [30] to study if multiple process observables lead to a reproducible (i.e., "unique") identification of material model parameters for the domain of investigation.
Three conditions were defined as termination criteria for the optimization: (i) Error criterion: the error value, Eq. (13), undercuts a value of ξ ≤ 1% for all cutting conditions. (ii) Convergence criterion I: the individual parameters of two consecutive sets deviate by less than 1%. (iii) Convergence criterion II: the improvement in the error function value is less than Δξ ≤ 0.1% for two consecutive iterations.
The initial JC parameter set ("target parameter set") was used to simulate the orthogonal cutting process for different cutting conditions (Table 3). Multiple cutting conditions were used for the inverse re-identification since it was expected to enhance the reliability of the model parameters due to the increased domain of investigation and the increased spectrum of occurring thermo-mechanical loads [16,20,25].   Figure 8 shows the results of the DSA starting with simplex 1, where the resulting error values ξ for the three cutting conditions are plotted over the iterations N. The diagram shows an exponentially decreasing error value. The optimization terminated after 82 iterations due to convergence criterion I (less than 1 % change in the parameters for consecutive iterations). For the conducted approach, the lowest error values were obtained at iteration 43, wherefore the average error value for the three cutting conditions was ξ aver = 1.1%. The error values of the individual cutting conditions were ξ 50 = 1.66%, ξ 100 = 0.98%, and ξ 150 = 0.80%. Besides the good fit, the obtained parameters deviated distinctively from the target parameters. For the domain of investigation, this does not support the hypothesis of a unique material model parameter set nor that using multiple cutting conditions foster unique parameter sets.
Analyzing the evolution of the error values reveals that the algorithm gets stuck in a local minimum. The global minimum is given by the target parameter set. The determination of a local minimum can mainly be attributed to two influencing factors: to the DSA, which is a local optimization algorithm, and to the existence of multiple local minima.
For the initial simplex 2, the results are shown in Fig. 9. It shows lower error values from the beginning. It seems that this set of random parameters better describes the three different cutting conditions. The approach was terminated after 92 iterations, due to convergence criterion II. The lowest error values were calculated for iteration 80, with an average error value of ξ aver = 1.6%. The individual error values of the three investigated cutting condition were ξ 50 = 1.43%, ξ 100 = 1.23%, and ξ 150 = 2.12%.
Again, the target parameter set was not obtained. By using two different initial simplexes, two different parameter sets were determined, which both produced low error values. Figure 10 provides the evolution of the DSA for the initial simplex 3. The inverse identification was terminated after 92 iterations, due to the second termination criteria (convergence criterion I). In this case, the lowest parameter set was determined by iteration 91. For the calculated parameter set, the average error value was ξ aver = 0.99%. The error values of the underlying cutting conditions were ξ 50 = 1.11%, ξ 100 = 0.66%, and ξ 150 = 1.18%.
The parameter set deviated from the target parameter set as well as from the determined parameter sets of the other two

Discussion
In the iterative procedure of this study, the deviation of a process observable from its target value was assessed by means of the merged error value ξ, Eq. (13). To evaluate the accuracy of the developed method, the relative errors of the individual process observables in comparison to their target values are given in Table 4 for the cutting condition v c = 100 m/min, and h = 0.05 mm. The highest deviations occurred in parameter set 3, wherefore the cutting normal force deviated by 1.2% from its target value. Numerical inaccuracies can only be held accustomed for minor differences in the order of a tenth of a percent. The chip form parameter showed the largest deviations. The chip radius deviated at most by 3.8% from its target value. The approaches revealed that the chip radius r c is very sensitive to the material model parameters.
Even slight changes in the material model parameters have a significant effect on the chip radius. The investigations suggested that the model parameter sets were not unique-at least in the considered domain of cutting conditions and for the considered process observables. In all cases, the error value was smaller than 2% suggesting the existence of multiple suitable local minima. To assess the accuracy of the determined parameter sets, the calculated material model parameters and the target parameter set are summarized in Table 5. The largest deviations occurred in the material model parameters A, C, and m. For these parameters, the differences can be larger than 100%. This results in a drastically different material behavior-respectively flow stress-with regard to the thermo-mechanical loads. There are multiple parameter sets, which may serve equally good to model the material behavior for the specific loads for the investigated domain. Figure 11 shows the resulting flow stresses separately listing the three influences on the flow stress that are considered by the JC model (strain hardening, strain rate hardening, and thermal softening). In the first term of the JC material model, the flow stress is calculated for different strains, whereas the   Parameter set 2 Parameter set 3 Fig. 11 Comparison of the material behavior for the three different terms of the Johnson-Cook material model representing the strain hardening, strain rate hardening, and thermal softening for the target parameter set and the determined parameter sets effect of strain rate hardening and thermal softening is considered by a hardening factor Kε and a softening factor K T . The curves show that for the identified parameter sets, the effect of strain hardening, strain rate hardening, and thermal softening are modeled differently. Among the three determined parameter sets, parameter set 2 resulted in the lowest flow stress for increasing strains. This is mainly attributed to the low value of parameter A, which defines the initial flow stress. On the contrary, parameter set 2 resulted in the highest strain rate hardening factor, since parameter C was higher than for the other parameter sets. Parameter set 3 on the other side resulted in the highest flow stress for increasing strains and in the lowest thermal softening factor K T . The differences in the magnitudes emphasize that the three effects can compensate each other within the investigated regime of machining conditions. Thereby, increased strain or strain rate hardening can be compensated by a higher thermal softening effect and vice versa. This observation can be considered as a drawback of the JC model for machining simulations, due to its uncoupled nature [45]. Material models with a coupled term, which takes the mutual influence of the effects into account, could solve this drawback. However, the reciprocal compensation of the three effects is expected to be reduced by considering experimental results from different domains, such as from quasi-static, dynamic material tests (e.g., SHPB tests), and cutting tests. From conventional material tests, different parameters can be determined directly by means of a curve fitting routine. This procedure would further reduce the number of parameters to be inversely calculated.
The inverse determination based on the minimization of the error value function ξ, Eq. (13). To compare the simulated results of the determined parameter sets regarding the occurring loads, Fig. 12, Fig. 13, and Fig. 14 show the results of the chip formation simulations for the same cutting conditions. Figure 12 shows the results regarding the local temperature field, Fig. 13 the plastic strains, and Fig. 14 the equivalent stresses. The comparison of simulated local temperatures, Fig.  12, revealed small deviations in terms of the absolute temperature and the expansion of the temperature fields for the different material model parameter sets. Thereby, the deviations between parameter set 1 and the target parameter set as well as between parameter set 2 and parameter set 3 were almost not identifiable. In comparison with the temperature field of the target parameter set, parameter set 2 and parameter set 3 resulted in a spatially more pronounced temperature profile.
The simulated chip form tended to deviate between the approaches, especially at the beginning of the chip formation. This is due to the individual differences regarding the three considered hardening and softening effects within the nonsteady state of the machining process. The differences of the chip form became smaller for longer cutting times and the more the system gets into the steady state.
As further material load, the plastic strain ϵ affects the flow stress. The comparison of the volume-averaged plastic strain, Fig. 13, revealed just slight differences. For the target parameter set and parameter set 1, the deviations regarding the plastic strain are not identifiable. Both, the magnitude and the local position of the volume-averaged plastic strain match each other. On the other side, small variations can be seen for the results of parameter set 2 and parameter set 3, which again are very similar. Based on the absolute differences regarding the material model parameters, the deviations regarding the local plastic strains ϵ are only small. This contributes to the assumption  Figure 14 shows the equivalent stress σ eq in and around the shear zone. In contrast to the comparisons of the temperature and the strain fields, the images from the simulations in Fig. 14 zoom into the cutting zone in order to enable a comparison of the local, highly resolved stress fields for the different parameter sets. As before, the equivalent stresses did not differ noticeable between the target parameter set and parameter set 1. On the other hand, small deviations can be identified for parameter set 2 and parameter set 3, which are expected to cause the differences in the considered process observables. A better agreement regarding the local stress fields is presumed for smaller error values.
The comparison of the simulation results showed that the different material parameter sets can all accurately model the locally resolved material loads. This emphasizes that uniqueness of material model parameters is no requirement to model the chip formation process within the domain of investigation of this study.

Summary, conclusion, and outlook
Within this paper, a novel approach to automatize the evaluation of chip formation simulations was developed. The methodology covers the determination of the chip form parameters chip thickness and chip radius from CEL chip formation simulations. The calculation of the chip form parameters proved a high accuracy. Additional procedures to average the force signal within the steady state and to compute local temperature information were established and integrated into an inverse procedure to determine material model parameters.
In the second part of the paper, the uniqueness of material model parameters for machining simulations was investigated. For this, an inverse re-identification by means of an optimization algorithm (Downhill simplex) was used on the automatized analysis of simulation results. In combination, they enable an automatic, inverse calibration of material model parameters from FE-cutting simulations. It was shown that different parameter sets can result in almost identical temperature, strain, and stress profiles. Based on the numerical results of this paper, the hypothesis of a unique set of JC material model parameters for the domain of cutting conditions cannot be supported. This was reinforced by the evaluated material loads.
The relevant findings are: & The identification is sensitive to the initial parameter set (simplex) but manages to converge to an accurate solution in all considered cases. & In spite of considering five process observables, the identification remains non-unique. This was attributed to the search space, which exhibits multiple local minima. & The differences in the sets affect all effects modeled by the JC material model: strain hardening, strain rate hardening, and thermal softening and thus may compensate each other. & The uniqueness of material model parameters is no requirement to model the chip formation process within the domain of investigation of this study. & The non-uniqueness of the JC material model may be considered as an inherent drawback.
Based on the results of this study, the need for further research activities was identified. Firstly, the proposed methodology needs to be applied to other material models than the JC model in order to investigate the generic nature of the methodology. In this study, the JC material model was chosen due its wide applicability to model the material behavior under metal cutting conditions. However, in the state of the art, limitations for the JC model were identified, which led to the development of new material models [46]. Secondly, the automatized evaluation of the chip form was limited to the continuous chip formation under the assumption of an isotropic material. However, materials that tend to the formation of continuous chips represent only a fraction of the machined materials in industry. To extend the methodology of the inverse parameter determination to materials that tend to segmented chip formation, the automatized evaluation needs to be enhanced and the inverse routine needs to be expanded to take the parameters of a damage model into account. For anisotropic materials, a separate calibration of the material model parameters for each phase will be necessary [47]. In future studies, the methodology will be applied to experimental data in order to determine the material models and to evaluate the accuracy of a material model to describe the material behavior