Predicting Continuous Chip to Segmented Chip Transition in Orthogonal Cutting of C45E Steel through Damage Modeling

: Machining process modeling has been an active endeavor for more than a century and it has been reported to be able to predict industrially relevant process outcomes. Recent advances in the fundamental understanding of material behavior and material modeling aids in improving the sustainability of industrial machining process. In this work, the ﬂow stress behavior of C45E steel is modeled by modifying the well-known Johnson-Cook model that incorporates the dynamic strain aging (DSA) inﬂuence. The modiﬁcation is based on the Voyiadjis-Abed-Rusinek (VAR) material model approach. The modiﬁed JC model provides the possibility for the ﬁrst time to include DSA inﬂuence in chip formation simulations. The transition from continuous to segmented chip for varying rake angle and feed at constant cutting velocity is predicted while using the ductile damage modeling approach with two di ﬀ erent fracture initiation strain models (Autenrieth fracture initiation strain model and Karp fracture initiation strain model). The result shows that chip segmentation intensity and frequency is sensitive to fracture initiation strain models. The Autenrieth fracture initiation strain model can predict the transition from continuous to segmented chip qualitatively. The study shows the transition from continuous chip to segmented chip for varying feed rates and rake angles for the ﬁrst time. The study highlights the need for material testing at strain, strain rate, and temperature prevalent in the machining process for the development of ﬂow stress and fracture models.


Introduction
The machining process remains one of the critical manufacturing processes in the 21st century and it has critical engineering applications [1]. Developments in the fields of plasticity and fracture mechanics are used to improve the understanding of the machining process [2]. This improved understanding has become even more necessary with the ever-increasing reliability need for engineered components [3].
Machining process modeling has been carried out while using different approaches, such as empirical, analytical and numerical methods [4]. The finite element (FE) method has been quite extensively used in the modeling of the machining process to model the chip formation process. FE modeling of the machining process provides the ability to incorporate a newer advanced understanding of material behavior. This improved understanding is usually obtained through other methods, such as material testing or more sophisticated modeling techniques. Within FE modeling of the machining process, workpiece material modeling requires the material response to large deformations (large plastic strains) at very high strain rates and very high temperatures. The fitting of the strain hardening component is traditionally carried out using standard quasistatic testing using tensile or compressive loading. In the previously reported work, A and B were modeled using flow stresses at varying temperatures and strain rates using regression modeling approach. The strain hardening exponent, n, was fitted as a first-order function of temperature. The strain hardening behavior of the JC model itself is written in the form, as shown in Equation (2) to reduce the FE implementation complexity.
ε 0 (4) The initial yield stress parameter A is modified similarly to the VAR model to incorporate the DSA influence, leading to yield stress increase with temperature increase at specific temperature ranges [8], as in Equations (3)- (6) and shown in Figure 1. Reported experimental investigations have shown that the dynamic strain aging temperature range is a function of strain rate [6,7,14]. The value of A 0 is obtained from the room temperature quasistatic compression test. The value of A ∆ corresponds to the peak initial yield stress increase in the DSA regime. The constants A 1 and t are fitted with flow stress curves that are obtained at high strain rate conditions from El Magd et al. [6] and Hokka et al. [7]. The temperatures T l and T h represent the start and end of DSA, the regime at the reference strain rate. T l and T h are modeled as a function of strain rate instead of a constant as in the VAR model. This methodology is necessary to accommodate the observation of dynamic strain aging temperature range being dependent on strain rate [6,7]. The model assumes the DSA regime to extend to higher strain rates at which experimental results are unavailable. Nevertheless, the DSA regime for machining process modeling has been incorporated by modifying the temperature component of the Power-law model while using a regression equation by Childs et al. [26], with the limitation of the temperature range to be a constant. The constants, B/A and n, are fitted for three different temperature ranges to accommodate the influence of temperature on strain hardening. The fitting of the strain hardening component is traditionally carried out using standard quasistatic testing using tensile or compressive loading. In the previously reported work, A and B were modeled using flow stresses at varying temperatures and strain rates using regression modeling approach. The strain hardening exponent, n, was fitted as a first-order function of temperature. The strain hardening behavior of the JC model itself is written in the form, as shown in Equation (2) to reduce the FE implementation complexity.
The initial yield stress parameter A is modified similarly to the VAR model to incorporate the DSA influence, leading to yield stress increase with temperature increase at specific temperature ranges [8], as in Equations (3)- (6) and shown in Figure 1. Reported experimental investigations have shown that the dynamic strain aging temperature range is a function of strain rate [6,7,14]. The value of A0 is obtained from the room temperature quasistatic compression test. The value of AΔ corresponds to the peak initial yield stress increase in the DSA regime. The constants and are fitted with flow stress curves that are obtained at high strain rate conditions from El Magd et al. [6] and Hokka et al. [7]. The temperatures and represent the start and end of DSA, the regime at the reference strain rate. and are modeled as a function of strain rate instead of a constant as in the VAR model. This methodology is necessary to accommodate the observation of dynamic strain aging temperature range being dependent on strain rate [6,7]. The model assumes the DSA regime to extend to higher strain rates at which experimental results are unavailable. Nevertheless, the DSA regime for machining process modeling has been incorporated by modifying the temperature component of the Power-law model while using a regression equation by Childs et al. [26], with the limitation of the temperature range to be a constant. The constants, ⁄ and , are fitted for three different temperature ranges to accommodate the influence of temperature on strain hardening.

Fitting of Flow Stress Curves Using the Modified JC Model
The ability of the newly developed flow stress curves from the previous section to predict DSA has been validated using the Gleeble test data from previous work [14]. The flow stress curves are obtained for three different strain rates of 1 s −1 , 5 s −1 , and 60 s −1 with programmed temperatures between 200-700 • C. The temperature range at which DSA is exhibited is a function of the strain rate [6,7]. In the temperature-strain rate range where DSA is active, the initial yield stress increases with temperature and decreases with the strain rate. The presence of DSA has also been shown to lead to serrations in the flow stresses. The temperature at which DSA peaks has been identified to be around 650 • C at high strain rates present in machining conditions [12]. Table 1 presents the constants of the modified JC model. The parameter B and n in the modified Johnson-Cook model has been modified for three different strain hardening regimes. This approach is carried out to ease the implementation of the FE method. The constants are provided, as shown in Equation (7).  The strain rate hardening and thermal softening parameters of the original JC model for C45E steel are used in this study (C = 0.0018 and m = 1) [11]. With the parameters obtained from the MJC model, the flow stress curves at a strain of 0.1 are predicted under different temperatures and three different strain rates are shown, as in Figure 2. The flow stress prediction extrapolated to higher strain rates prevalent in machining conditions is also plotted to visualize the DSA regime being a function of strain rate.

Fitting of Flow Stress Curves Using the Modified JC Model
The ability of the newly developed flow stress curves from the previous section to predict DSA has been validated using the Gleeble test data from previous work [14]. The flow stress curves are obtained for three different strain rates of 1 s −1 , 5 s −1 , and 60 s −1 with programmed temperatures between 200-700 °C. The temperature range at which DSA is exhibited is a function of the strain rate [6,7]. In the temperature-strain rate range where DSA is active, the initial yield stress increases with temperature and decreases with the strain rate. The presence of DSA has also been shown to lead to serrations in the flow stresses. The temperature at which DSA peaks has been identified to be around 650 °C at high strain rates present in machining conditions [12]. Table 1 presents the constants of the modified JC model. The parameter B and n in the modified Johnson-Cook model has been modified for three different strain hardening regimes. This approach is carried out to ease the implementation of the FE method. The constants are provided, as shown in Equation Error! Reference source not found.).
The strain rate hardening and thermal softening parameters of the original JC model for C45E steel are used in this study (C = 0.0018 and m = 1) [11]. With the parameters obtained from the MJC model, the flow stress curves at a strain of 0.1 are predicted under different temperatures and three different strain rates are shown, as in Figure 2. The flow stress prediction extrapolated to higher strain rates prevalent in machining conditions is also plotted to visualize the DSA regime being a function of strain rate.  The modification of the initial yield stress parameter of the JC model can capture the increase in initial yield stress with temperature increase. At the reference strain rate of 1 s −1 , the DSA regime is active in the temperature range of 200 • C to 500 • C. With a strain rate increase between 10 3 s −1 to Metals 2020, 10, 519 5 of 17 10 6 s −1 , the modeled DSA regime's initial yield stress moves to higher temperatures of 650 • C and it correlates with other published experimental results [6,7].
With the ability of the MJC model's initial yield stress parameter to capture the DSA regime in Figure 2, Figure 3 shows that the MJC model can capture the strain hardening behavior observed from Gleeble tests with reasonable accuracy. In the temperature ranges below 400 • C, the material strain hardens with increasing strain. At 500 • C, dynamic recovery influences strain hardening, leading to a lowering slope, and are well captured by the MJC model. Beyond 500 • C, at 600 • C and 700 • C, the flow stress curves exhibit constant flow stress with increasing strain.
Metals 2020, 10, x FOR PEER REVIEW 5 of 17 The modification of the initial yield stress parameter of the JC model can capture the increase in initial yield stress with temperature increase. At the reference strain rate of 1 s −1 , the DSA regime is active in the temperature range of 200 °C to 500 °C. With a strain rate increase between 10 3 s −1 to 10 6 s −1 , the modeled DSA regime's initial yield stress moves to higher temperatures of 650 °C and it correlates with other published experimental results [6,7].
With the ability of the MJC model's initial yield stress parameter to capture the DSA regime in Figure 2, Figure 3 shows that the MJC model can capture the strain hardening behavior observed from Gleeble tests with reasonable accuracy. In the temperature ranges below 400 °C, the material strain hardens with increasing strain. At 500 °C, dynamic recovery influences strain hardening, leading to a lowering slope, and are well captured by the MJC model. Beyond 500 °C, at 600 °C and 700 °C, the flow stress curves exhibit constant flow stress with increasing strain.  The modification of the initial yield stress parameter of the JC model can capture the increase in initial yield stress with temperature increase. At the reference strain rate of 1 s −1 , the DSA regime is active in the temperature range of 200 °C to 500 °C. With a strain rate increase between 10 3 s −1 to 10 6 s −1 , the modeled DSA regime's initial yield stress moves to higher temperatures of 650 °C and it correlates with other published experimental results [6,7].
With the ability of the MJC model's initial yield stress parameter to capture the DSA regime in Figure 2, Figure 3 shows that the MJC model can capture the strain hardening behavior observed from Gleeble tests with reasonable accuracy. In the temperature ranges below 400 °C, the material strain hardens with increasing strain. At 500 °C, dynamic recovery influences strain hardening, leading to a lowering slope, and are well captured by the MJC model. Beyond 500 °C, at 600 °C and 700 °C, the flow stress curves exhibit constant flow stress with increasing strain.

Flow Stress Modification due to Damage
JC and MJC models both assume no change in material behavior with increasing strain. However, in reality, the flow stress is altered because of fracture. Therefore, to model chip segmentation in ductile materials at a lower cutting velocity where the possibility of adiabatic shear is less incorporation of fracture is essential. Figure 5a provides the schematic modification of the flow stress curves due to fracture from the works of Childs et al. [16,26]. The flow stress before the fracture is defined by Equation (8) and is marked by 'BEFORE FRACTURE' in Figure 5a. During plastic deformation, at the microstructure level, nucleation and growth of defects, such as micro-voids & micro-cracks and their coalescence into macro-cracks takes place leading to material damage [27]. The damage due to plastic deformation in the workpiece is accumulated through the damage factor (D) and is defined, as shown in Equations (8)

Flow Stress Modification due to Damage:
JC and MJC models both assume no change in material behavior with increasing strain. However, in reality, the flow stress is altered because of fracture. Therefore, to model chip segmentation in ductile materials at a lower cutting velocity where the possibility of adiabatic shear is less incorporation of fracture is essential. Figure 5a provides the schematic modification of the flow stress curves due to fracture from the works of Childs et al. [16,26]. The flow stress before the fracture is defined by Equation (8) and is marked by 'BEFORE FRACTURE' in Figure 5a. During plastic deformation, at the microstructure level, nucleation and growth of defects, such as micro-voids & micro-cracks and their coalescence into macro-cracks takes place leading to material damage [27]. The damage due to plastic deformation in the workpiece is accumulated through the damage factor (D) and is defined, as shown in Equations (7)  The damage parameter is defined as the ratio of accumulated plastic strain to fracture initiation strain ( ). The fracture initiation strain is presented in detail in the following section and it is the focus of the study. As the accumulated plastic strain equals fracture initiation strain, the damage parameter equals one and it is shown by "FRACTURE" in Figure 5a. The damage parameter is defined as the ratio of accumulated plastic strain to fracture initiation strain (ε f i ). The fracture initiation strain is presented in detail in the following section and it is the focus of the study. As the accumulated plastic strain equals fracture initiation strain, the damage parameter equals one and it is shown by "FRACTURE" in Figure 5a.
q(η, T, D, The modeling of the flow stress behavior after the fracture point is defined by "AFTER FRACTURE" in Figure 5a. The flow stress curve is modified based on the loading conditions, η, and temperature, T, as shown in Equation (8) using the flow stress modification factor, f (η, T). In an earlier development of Childs et al. [16], f (η, T) was assumed to be a constant. With further development [26], f (η, T) is defined as a function of stress triaxiality and temperature, as shown in Equation (8) and Figure 5c and presented in detail in the following sections.

Fracture Initiation Strain
The need for the fracture initiation strain and its use in the damage factor is presented in the previous section. The fracture initiation strain is modeled as a function of stress triaxiality (η), lode angle parameter θ , temperature (T), and strain rate . ε , as shown in Equation (13).
The function, f η, θ , is used to define the loading conditions. The stress triaxiality (η) defines the varying loading conditions in two-dimensional (2D) and the lode angle parameter, θ, extends the model to loading in three-dimensional (3D) space. In the case of the 2D cutting process (orthogonal cutting process), which falls under the plane strain condition, the lode angle parameter is zero. The stress triaxiality parameter defined by, η, as σ m σ defines the varying loading conditions from pure compression to combinations of shear/compression and shear/tension to tension, as shown in Figure 6. The modeling of the flow stress behavior after the fracture point is defined by "AFTER FRACTURE" in Figure 5a. The flow stress curve is modified based on the loading conditions, , and temperature, T, as shown in Equation (8) using the flow stress modification factor, ( , ). In an earlier development of Childs et al. [16], ( , ) was assumed to be a constant. With further development [26], ( , ) is defined as a function of stress triaxiality and temperature, as shown in Equation (8) and Figure 5c and presented in detail in the following sections.

Fracture Initiation Strain
The need for the fracture initiation strain and its use in the damage factor is presented in the previous section. The fracture initiation strain is modeled as a function of stress triaxiality ( ), lode angle parameter ( ̅ ), temperature ( ), and strain rate (), as shown in Equation (13).
The function, ( , ̅ ), is used to define the loading conditions. The stress triaxiality ( ) defines the varying loading conditions in two-dimensional (2D) and the lode angle parameter, ̅ , extends the model to loading in three-dimensional (3D) space. In the case of the 2D cutting process (orthogonal cutting process), which falls under the plane strain condition, the lode angle parameter is zero. The stress triaxiality parameter defined by, , as defines the varying loading conditions from pure compression to combinations of shear/compression and shear/tension to tension, as shown in Figure  6. In order to obtain the fracture initiation strain under different conditions of (−1 to 1), specimens of varying geometrical shapes are to be used [28]. Interestingly, the fracture initiation strain for variants of C45E steel has been studied for nearly 75 years [29]. With the advances in material testing, the variants of the C45E steel's fracture strain has been studied. Autenrieth et al. [30] obtained the fracture strain as a function of stress triaxiality using torsion and torsion-tension tests. These tests provide stress triaxiality from zero to positive stress triaxiality conditions. Recently, the same fracture strain as a function of stress triaxiality for conditions that range from negative to positive has been obtained using the Shear-compression disk specimen by Karp et al. [31]. In this study, the latter two fracture strain models identified as Autenrieth fracture strain data and Karp fracture strain data are used to evaluate their influence on chip segmentation prediction. The difference in the fracture strain between the two models can be attributed to the difference in processing history and the stress triaxiality evolution before fracture. Figure 6. Fracture initiation strain predicted using Shear compression specimen under quasi-static loading, Shear specimens and compression testing as a function of stress triaxiality and fitted using an exponential function.
In order to obtain the fracture initiation strain under different conditions of η (−1 to 1), specimens of varying geometrical shapes are to be used [28]. Interestingly, the fracture initiation strain for variants of C45E steel has been studied for nearly 75 years [29]. With the advances in material testing, the variants of the C45E steel's fracture strain has been studied. Autenrieth et al. [30] obtained the fracture strain as a function of stress triaxiality using torsion and torsion-tension tests. These tests provide stress triaxiality from zero to positive stress triaxiality conditions. Recently, the same fracture strain as a function of stress triaxiality for conditions that range from negative to positive has been obtained using the Shear-compression disk specimen by Karp et al. [31]. In this study, the latter two fracture strain models identified as Autenrieth fracture strain data and Karp fracture strain data are used to evaluate their influence on chip segmentation prediction. The difference in the fracture strain between the two models can be attributed to the difference in processing history and the stress triaxiality evolution before fracture.
To implement the fracture strain in the finite element model, the exponential function used in the Bai-Wierzbicki damage model [32] is used and it is shown in Equation (14).
The Autenrieth fracture strain data are fitted using the parameters (D X1 = 0.2339 and D X2 = −1.035) and the Karp fracture strain data are fitted using the parameters (D X1 = 0.7 and D X2 = −0.9).
The function g .
ε is used to model the influence of strain rate. h(T) is used to model the influence of temperature on the fracture initiation strain. In this study, g . ε and h(T) are obtained from the JC damage fracture model, as shown in Equation (15) and Equation (16) and D 3 = 0.0018 and D 4 = 0.58 are obtained from the literature [33].

Flow Stress Modification Factor
As previously mentioned, the flow stress is modified, as shown in Figure 5a using the flow stress modification factor, f (η, T) shown in Figure 5c. The flow stress in the 'AFTER FRACTURE' region is modeled while using a steady-state flow stress modification factor defined by Equation (9) and a transient flow stress modification factor. The steady-state is defined by D > D+ and the transient state is defined by D+ > D > 1. In this study, D+ is defined to be 1.25, as suggested by Childs et al. [26].
In the steady-state, the flow stress modification is a function of stress triaxiality and it is defined by Equation (12) below a certain temperature, T S , as shown in Figure 5c. A negative stress triaxiality condition (e.g., η = −1) characterizing compression, the flow stress is modified only slightly due to incompressibility. A positive stress triaxiality condition, e.g., η = 1 characterizing tension, the flow stress is drastically reduced to zero, similar to necking in tensile testing. At temperatures above T E , the material is assumed to be healed and continue to flow plastically with the flow stress defined by σ, as shown in Equation (8). Between temperatures T S and T E , the flow stress curve is defined using the tanh function, as shown in Equation (11). The transient flow stress modification factor is defined as Equation (9) by multiplying the steady-state flow stress modification factor with a tanh function of D and f (η, T).
The values for µ i , T S , and T E used in this work are 1 • C, 600 • C, and 700 • C, respectively. The values for T S and T E based on the strain hardening transition behavior observed from the Gleeble tests.

Experimental Investigation of Chip Formation in Orthogonal Turning
The orthogonal turning process has been carried out in a normalized C45E steel tube with a diameter of 150 mm with a thickness of 3 mm, determining the uncut chip width. The material hardness is measured and an average hardness of 220 HV is recorded with an average grain size of 13.5 µm. The tubular workpiece reduces grain size variation influence. All of the experiments were carried out under dry cutting conditions. The cutting tool material is H13A grade carbide with TiCN coating. The cutting velocity is set at 150 m/min. The rake angle is varied as −5 • , 5 • , 10 • , and 20 • , and the feed is varied as 0.05 mm.rev −1 , 0.1 mm.rev −1 , 0.15 mm.rev −1 , 0.25 mm.rev −1 , 0.4 mm.rev −1 , and 0.6 mm.rev −1 . The cutting forces and feed forces are obtained through the Kistler dynamometer attached to the tool turret. The sample chips were collected for all cutting conditions and they have been used to quantify chip segmentation, as shown in Figure 7. The experimental investigation is described in detail in previous work [34]. The chips have to be plotted with varying magnification to capture the chip segmentation features at very low feed rate conditions (e.g., 0.05 mm.rev −1 ) and similarly at very high feed rate conditions (e.g., 0.60 mm.rev −1 ). This approach with varying magnification is carried out when the chips from simulations are also presented in the following sections.

FE Modeling of Chip Formation in Orthogonal Turning Process
The FE simulation of the cutting process was carried out using a commercial finite element software, Thirdwave Advantedge [35]. The software has been specifically built ground up for metal cutting simulations while using a dynamic explicit coupled Thermo-elastoplastic Lagrangian formulation [36].
The main requirement for metal cutting simulation is the ability for adaptive remeshing to account for the large plastic strains in the primary and secondary deformation zones. The simulations are carried out in 2D. The software provides the ability to customize the adaptive remeshing parameters. The adaptive remeshing parameters control the mesh transition in the tool vicinity (primary deformation zone) from coarse to fine and also the mesh transition from fine to coarse (chip mesh). The software through the mesh refinement and coarsening factors automatically control the conditions for adaptive remeshing. With predicting chip segmentation being the main aim of the study, mesh refinement factor of 6 (1-coarse → 8-fine) and mesh coarsening factor of 3 (1-fine → 8-coarse) is chosen.
A further increase in mesh refinement factors leads to computation time to increase by a factor of 3, making it practically unviable. Besides, the suggested minimum element size is set at 0.01 mm. The friction coefficient of 0.5 is used for all cutting conditions. The software provides the possibility to input the custom material model through a FORTRAN Subroutine, which and has been employed in the study to implement the modified JC model. Two different solution algorithms are supported by the software: Newton method and the secant method. At this stage, the secant method is used, as the method does not require the derivative implementation with the compromise of a slower convergence rate.
Simulations are carried out for all rake angles (−5°, 5°, 10° and 20°) and feed from 0.05 mm.rev −1 to 0.6 mm.rev −1 . The relief angle and the cutting-edge radius are kept constant with the values 7° and 30 µm, respectively, both in the experimental conditions and simulations. Simulations are run with the JC model and MJC model combined with the damage modeling approach. Within the MJC model + damage model framework, the influence of two fracture initiation strain model's (Autenrieth model and Karp model) capability to evaluate the transition from continuous chip to segmented chip is studied.

Results
The FE simulations were carried out using two different material models (JC and MJC). With the MJC+ damage model, two different parameter sets of fracture initiation strain models' ability to predict the continuous chip to segmented chip transition are evaluated. The FE simulations were run

FE Modeling of Chip Formation in Orthogonal Turning Process
The FE simulation of the cutting process was carried out using a commercial finite element software, Thirdwave Advantedge [35]. The software has been specifically built ground up for metal cutting simulations while using a dynamic explicit coupled Thermo-elastoplastic Lagrangian formulation [36].
The main requirement for metal cutting simulation is the ability for adaptive remeshing to account for the large plastic strains in the primary and secondary deformation zones. The simulations are carried out in 2D. The software provides the ability to customize the adaptive remeshing parameters. The adaptive remeshing parameters control the mesh transition in the tool vicinity (primary deformation zone) from coarse to fine and also the mesh transition from fine to coarse (chip mesh). The software through the mesh refinement and coarsening factors automatically control the conditions for adaptive remeshing. With predicting chip segmentation being the main aim of the study, mesh refinement factor of 6 (1-coarse → 8-fine) and mesh coarsening factor of 3 (1-fine → 8-coarse) is chosen.
A further increase in mesh refinement factors leads to computation time to increase by a factor of 3, making it practically unviable. Besides, the suggested minimum element size is set at 0.01 mm. The friction coefficient of 0.5 is used for all cutting conditions. The software provides the possibility to input the custom material model through a FORTRAN Subroutine, which and has been employed in the study to implement the modified JC model. Two different solution algorithms are supported by the software: Newton method and the secant method. At this stage, the secant method is used, as the method does not require the derivative implementation with the compromise of a slower convergence rate.
Simulations are carried out for all rake angles (−5 • , 5 • , 10 • and 20 • ) and feed from 0.05 mm.rev −1 to 0.6 mm.rev −1 . The relief angle and the cutting-edge radius are kept constant with the values 7 • and 30 µm, respectively, both in the experimental conditions and simulations. Simulations are run with the JC model and MJC model combined with the damage modeling approach. Within the MJC model + damage model framework, the influence of two fracture initiation strain model's (Autenrieth ε f i model and Karp ε f i model) capability to evaluate the transition from continuous chip to segmented chip is studied.

Results
The FE simulations were carried out using two different material models (JC and MJC). With the MJC+ damage model, two different parameter sets of fracture initiation strain models' ability to predict the continuous chip to segmented chip transition are evaluated. The FE simulations were run till a steady state was achieved, and the chip morphology was obtained. The chip morphology was plotted in the chip chart form. Further on, the continuous chip-segmented chip boundary is identified in the chip chart. The cutting forces predicted by the material models are also presented.

Validation of Material Model under Continuous Chip Formation Conditions
This section presents the MJC model validation in the simulation of the chip formation process. Experimental results with continuous chips are used to avoid the damage model influence. The rake angle of 10 • was used with feeds of 0.05 mm.rev −1 to 0.6 mm.rev −1 and Figure 8 presents the cutting forces.
Metals 2020, 10, x FOR PEER REVIEW 10 of 17 till a steady state was achieved, and the chip morphology was obtained. The chip morphology was plotted in the chip chart form. Further on, the continuous chip-segmented chip boundary is identified in the chip chart. The cutting forces predicted by the material models are also presented.

Validation of Material Model under Continuous Chip Formation Conditions
This section presents the MJC model validation in the simulation of the chip formation process. Experimental results with continuous chips are used to avoid the damage model influence. The rake angle of 10° was used with feeds of 0.05 mm.rev −1 to 0.6 mm.rev −1 and Figure 8 presents the cutting forces. The experimental result showing the increase in cutting force with feed is well established within the metal cutting with the influence of uncut chip thickness increase. The cutting forces predicted by the JC model and the MJC model have a similar trend. The MJC model improves the prediction for all of the cutting conditions quantitatively. This improvement in cutting force prediction as compared to the JC model is attributed to the material testing used and the ability of the MJC model to incorporate the material behavior with improved accuracy. The highest discrepancies in quantitative prediction for both JC and MJC models are at the maximum chip thickness conditions. The MJC model's cutting force prediction error can be attributed to the extrapolation that was carried out with the fitting of the MJC model at very high strain rates. Further improvement has to be carried out, among other things, by improving the friction models.

Experimental and FE Modeling of Continuous Chip-Segmented Chip Transition for Varying Rake Angle and Feed
The chips obtained from experimental investigation for varying rake angle and feed rates at a constant cutting speed of 150 m/min are plotted in the form of chip chart in to evaluate the ability of the different fracture initiation strain models in chip segmentation prediction. Figure 7 shows that, at constant cutting speed, continuous chips are produced for rake angles 10° and 20°. With a rake angle of 5°, the cutting process is in a transition zone. Continuous chips are produced for feed of 0.05 mm.rev −1 and 0.10 mm.rev −1 and segmented chips are produced for feed from 0.15 mm.rev −1 to 0.6 mm.rev −1 . For a rake angle of −5°, the produced chips are segmented for all feeds. At a constant feed of 0.40 mm.rev −1 , a change in the rake angle from 5° to 10° changes the chips from being continuous to segmented.
Similarly, at a constant rake angle of 5°, the change in feed rate from 0.10 mm.rev −1 to 0.15 mm.rev −1 leads to continuous to segmented chip. This leads to the rake angle of 5° being a transition zone. The experimental results clearly show that chip segmentation/continuous chip is stable under certain cutting conditions (−5°, 10°, and 20°) and in transition mode under certain conditions (5°). It is also seen that, in the transition zone, the feed also plays a role in chip segmentation. This states that minor changes in the cutting zone to be highly sensitive to the machining output. This would The experimental result showing the increase in cutting force with feed is well established within the metal cutting with the influence of uncut chip thickness increase. The cutting forces predicted by the JC model and the MJC model have a similar trend. The MJC model improves the prediction for all of the cutting conditions quantitatively. This improvement in cutting force prediction as compared to the JC model is attributed to the material testing used and the ability of the MJC model to incorporate the material behavior with improved accuracy. The highest discrepancies in quantitative prediction for both JC and MJC models are at the maximum chip thickness conditions. The MJC model's cutting force prediction error can be attributed to the extrapolation that was carried out with the fitting of the MJC model at very high strain rates. Further improvement has to be carried out, among other things, by improving the friction models.

Experimental and FE Modeling of Continuous Chip-Segmented Chip Transition for Varying Rake Angle and Feed
The chips obtained from experimental investigation for varying rake angle and feed rates at a constant cutting speed of 150 m/min are plotted in the form of chip chart in to evaluate the ability of the different fracture initiation strain models in chip segmentation prediction. Figure 7 shows that, at constant cutting speed, continuous chips are produced for rake angles 10 • and 20 • . With a rake angle of 5 • , the cutting process is in a transition zone. Continuous chips are produced for feed of 0.05 mm.rev −1 and 0.10 mm.rev −1 and segmented chips are produced for feed from 0.15 mm.rev −1 to 0.6 mm.rev −1 . For a rake angle of −5 • , the produced chips are segmented for all feeds. At a constant feed of 0.40 mm.rev −1 , a change in the rake angle from 5 • to 10 • changes the chips from being continuous to segmented.
Similarly, at a constant rake angle of 5 • , the change in feed rate from 0.10 mm.rev −1 to 0.15 mm.rev −1 leads to continuous to segmented chip. This leads to the rake angle of 5 • being a transition zone. The experimental results clearly show that chip segmentation/continuous chip is stable under certain cutting conditions (−5 • , 10 • , and 20 • ) and in transition mode under certain conditions (5 • ). It is also seen that, in the transition zone, the feed also plays a role in chip segmentation. This states that minor changes in the cutting zone to be highly sensitive to the machining output. This would lead to challenges in the process capability in a production environment.

Prediction of Chip Segmentation Using MJC Model and Two Different Fracture Initiation Strain Models
Initial simulations were run with the JC flow stress model [37] and the JC fracture model [38] reported for AISI 1045 steel to predict chip segmentation. The continuous chip to segmented chip transition predicted by the JC flow stress-JC fracture model, as shown in Figure 9, with the JC fracture parameters that were obtained from the literature [38]. The JC flow stress-JC fracture model is unable to predict the influence of the rake angle and, at the same time, predicts lower chip segmentation intensity when compared to experimental investigation. With the inability of the JC flow stress-JC fracture model established, the flow stress behavior is modified while using the damage model to improve chip segmentation prediction. The chip morphology predicted by the MJC model with the damage model using the two fracture initiation strain models is presented here.

Prediction of Chip Segmentation Using MJC Model and Two Different Fracture Initiation Strain Models
Initial simulations were run with the JC flow stress model [37] and the JC fracture model [38] reported for AISI 1045 steel to predict chip segmentation. The continuous chip to segmented chip transition predicted by the JC flow stress-JC fracture model, as shown in Figure 9, with the JC fracture parameters that were obtained from the literature [38]. The JC flow stress-JC fracture model is unable to predict the influence of the rake angle and, at the same time, predicts lower chip segmentation intensity when compared to experimental investigation. With the inability of the JC flow stress-JC fracture model established, the flow stress behavior is modified while using the damage model to improve chip segmentation prediction. The chip morphology predicted by the MJC model with the damage model using the two fracture initiation strain models is presented here.  Figure 10 shows the transition of the chip morphology that was predicted by the Autenrieth fracture initiation strain model and MJC material model. The continuous chip to segmented chip transition is predicted with reasonable accuracy. At a constant feed of 0.6 mm.rev −1 , the model can predict the segmented chip for rake angles: −5° and 5° and the transition to continuous chip with a rake angle of 10°. The model can predict the transition from continuous chip to segmented chip as the feed is increased from 0.05 mm.rev −1 to 0.1 mm.rev −1 . At a feed of 0.05 mm.rev −1 , a continuous chip is predicted for both negative and positive rake angle. A qualitative comparison of chip segmentation between the experimental investigation and the Autenrieth fracture initiation strain model prediction shows that the level of chip segmentation intensity is lower when compared to the experimental results. This leads to a further quantitative evaluation in terms of chip segmentation intensity become counterproductive. Chip segmentation intensity has been predicted with improved accuracy in literature but only for a constant rake angle and varying cutting speeds [39]. The models in these studies have shown cutting speed's influence and not the rake angle's influence on chip segmentation.  Figure 10 shows the transition of the chip morphology that was predicted by the Autenrieth fracture initiation strain model and MJC material model. The continuous chip to segmented chip transition is predicted with reasonable accuracy. At a constant feed of 0.6 mm.rev −1 , the model can predict the segmented chip for rake angles: −5 • and 5 • and the transition to continuous chip with a rake angle of 10 • . The model can predict the transition from continuous chip to segmented chip as the feed is increased from 0.05 mm.rev −1 to 0.1 mm.rev −1 . At a feed of 0.05 mm.rev −1 , a continuous chip is predicted for both negative and positive rake angle. A qualitative comparison of chip segmentation between the experimental investigation and the Autenrieth fracture initiation strain model prediction shows that the level of chip segmentation intensity is lower when compared to the experimental results. This leads to a further quantitative evaluation in terms of chip segmentation intensity become counterproductive. Chip segmentation intensity has been predicted with improved accuracy in literature but only for a constant rake angle and varying cutting speeds [39]. The models in these studies have shown cutting speed's influence and not the rake angle's influence on chip segmentation.  Figure 11 shows the transition of the chip morphology predicted by the fracture initiation strain that was obtained by Karp et al. [31] while using shear compression disk specimen. The model predicts segmented chips for all rake angles for feeds above and including 0.1 mm.rev −1 . For the feed of 0.05 mm.rev −1 , the segmented chip for a negative rake angle is not predicted correctly compared to experimental results. Qualitative evaluation of the chip segmentation intensity when compared to the experimental investigation is relatively low for a negative rake angle. On the other hand, the chip segmentation frequency is higher qualitatively compared to the Autenrieth fracture initiation strain model predicted chip segmentation frequency. Figure 11. Continuous chip to segmented chip transition prediction by the fracture initiation strain obtained Shear compression specimen. Note: Image magnification not to scale. Note: Image magnification not to scale.  Figure 11 shows the transition of the chip morphology predicted by the fracture initiation strain that was obtained by Karp et al. [31] while using shear compression disk specimen. The model predicts segmented chips for all rake angles for feeds above and including 0.1 mm.rev −1 . For the feed of 0.05 mm.rev −1 , the segmented chip for a negative rake angle is not predicted correctly compared to experimental results. Qualitative evaluation of the chip segmentation intensity when compared to the experimental investigation is relatively low for a negative rake angle. On the other hand, the chip segmentation frequency is higher qualitatively compared to the Autenrieth fracture initiation strain model predicted chip segmentation frequency.  Figure 11 shows the transition of the chip morphology predicted by the fracture initiation strain that was obtained by Karp et al. [31] while using shear compression disk specimen. The model predicts segmented chips for all rake angles for feeds above and including 0.1 mm.rev −1 . For the feed of 0.05 mm.rev −1 , the segmented chip for a negative rake angle is not predicted correctly compared to experimental results. Qualitative evaluation of the chip segmentation intensity when compared to the experimental investigation is relatively low for a negative rake angle. On the other hand, the chip segmentation frequency is higher qualitatively compared to the Autenrieth fracture initiation strain model predicted chip segmentation frequency. Figure 11. Continuous chip to segmented chip transition prediction by the fracture initiation strain obtained Shear compression specimen. Note: Image magnification not to scale. Note: Image magnification not to scale. Figure 11. Continuous chip to segmented chip transition prediction by the fracture initiation strain obtained Shear compression specimen. Note: Image magnification not to scale. Note: Image magnification not to scale.

Discussion
The FE simulation of the machining process's ability to predict chip segmentation is dependent on the accuracy in which the material behavior is modeled to a large extent, as shown by Fernandez-Zeliaia et al. [40]. The other well-known influencing factors are friction behavior and thermal behavior input. Within flow stress behavior, the initial yield stress, the strain hardening behavior, the fracture initiation, and the material behavior after fracture are essential inputs. For the material under study, C45E steel, strain hardening behavior, and thermal softening behavior are more influential when compared to strain rate hardening [41]. The MJC model's ability to accurately predict three different strain hardening behavior as a function of temperature, therefore, systematically improves the accuracy of FE modeling of chip formation. The material behavior at the strain, strain rate, and temperature occurring in the chip formation conditions is vital in accurately predicting the chip morphology. In this study, the material flow stress is initially validated while using the flow stress that was obtained from compression testing through Gleeble tests (Figure 3). Although the temperatures at which the material behavior is validated through Gleeble tests are comparable, the strains and strain rates are lower by orders of magnitude. The material model validation through cutting process simulation is carried out by avoiding the cutting conditions where chip segmentation is not observed. The MJC model's DSA incorporation, which is active in temperatures occurring at the primary and secondary deformation zones, is an important outcome of this work. The ability to incorporate DSA keeping all other things constant has been shown to improve the cutting forces prediction accuracy, as shown in Figure 8. The fracture initiation models with the JC model have not been able to predict the chip segmentation boundary. This shows the influence of the MJC model on the improved prediction of chip segmentation boundary.
The damage behavior modeling involves fracture initiation strain and the flow stress after the fracture initiation strain. The results of Figure 9, Figure 10, and Figure 11 show that the continuous chip to segmented chip transition and chip segmentation intensity is sensitive to the fracture initiation strain input. The Autenrieth fracture initiation strain model better predicts the transition between the continuous chip and segmented chip when compared to the Karp fracture initiation strain model. Simulations run with the JC material model combined with the fracture initiation strain models have shown poor results in terms of prediction of the chip segmentation ( Figure 9). This shows that the ability of the fracture initiation strain models to predict chip segmentation is dependent on the components of the flow stress behavior, i.e., initial yield stress and strain hardening behavior. The strain hardening behavior and the resulting temperature increase influences the fracture initiation strain. The strain rate difference and the initial yield stress variation between the two experimental conditions (for the Autenrieth fracture initiation strain model and Karp fracture initiation strain model) have been reported as being possible causes for the significant variation in the fracture initiation strain. The significant fracture initiation strain variation correspondingly leads to significant variation in the continuous chip to segmented chip transition and the chip segmentation intensity prediction in the finite element models that were developed in this study.
The experimental investigation (Figure 7) shows that the chip segmentation at a constant cutting speed is influenced by tool geometry and chip thickness. In this study, the normal rake angle is the tool geometry parameter modified, and its influence on the chip formation is studied. With the normal rake angle being modified from negative to positive rake angle under certain uncut chip thickness (feed) conditions, the chip formation process is transformed from a continuous chip to a segmented chip, as observed with the segmented chip boundary.
With the normal rake angle modification, the mechanical loading of the workpiece ahead of the tool is modified with varying levels of shear and compression. The mechanical loading is not uniform along the shear plane from the cutting edge to the free surface in the primary deformation zone. The stress that was observed by the workpiece material close to the cutting-edge is highly compressive when compared to the free end of the primary deformation zone, where the stresses are a combination of shear and compression. The workpiece material that is close to the cutting-edge rounding is mainly under the stress state defined by combined shear and compression defined by a large negative stress triaxiality. The workpiece material on the free end of the primary deformation zone has a relatively lower negative stress triaxiality; still, they are in the negative stress triaxiality conditions. Figure 6 shows that the fracture initiation strain in the negative stress triaxiality condition to be sensitive and varies based on testing conditions. This fracture initiation strain sensitivity is attributed to be the main reason for chip segmentation to be sensitive to minor changes in cutting conditions. At a constant chip thickness and cutting speed, the tool geometry influences this mechanical loading and it is observed to influence the onset of chip segmentation. Another critical factor that influences the material behavior is the material's history. From a machining perspective, the influence of the previous cut on the strain profile has been shown by Childs et al. [26]. The tool edge geometry controls the influence of the previous cut on the stress profile on the workpiece surface, tool wear, and tool vibration. The tool edge geometry could also influence the transition zone between the continuous chip to segmented chip. In this work, the influence of the previous cut is not taken into account and it could be a source of improvement in chip segmentation prediction in future work.
It is clearly understood that the uncut chip thickness, in addition to the rake angle, influences the transition from continuous to segmented chip at constant cutting speed. This is understood from the results, showing that, at a constant rake angle, the chip thickness influences chip segmentation. The chip thickness increase leads to the heat being concentrated in the primary shear zone. This concentration of heat leads to a segmented chip, as predicted in the presented FE simulations. On the other hand, chip thickness increase needs to be significant in this case at a rake angle of 5 • up to 0.6 mm.rev −1 to produce segmented chips.

Conclusions
In this study, a new modified JC model is developed based on the Voyiadjis-Abed-Rusinek [8] constitutive modeling approach. The modified JC model with the fracture modeling approach by Childs et al. [26] has been implemented for the first time in the FE framework for the chip formation simulation. The influence of the fracture initiation strain model on the continuous to segmented chip transition prediction has been successfully demonstrated. The following conclusions are made from the study.

1.
The modified JC model can incorporate the DSA influence on flow stress curves of normalized C45E steel for varying temperatures and strain rates with reasonable accuracy.

2.
The strain hardening behavior modeling of the modified JC model can predict the temperature's influence on strain hardening accurately when compared to the JC model.

3.
The strong influence of DSA in the simulation of chip formation in machining of C45E steel is established with an improvement of cutting force prediction accuracy.

4.
The transition from continuous chip to segmented chip is established to be a function of normal rake angle and feed at a constant cutting speed in orthogonal cutting.

5.
The Autenrieth fracture initiation strain model can predict the chip segmentation boundary better as compared to the Karp fracture initiation strain model.