Influence of Constitutive Models and the Choice of the Parameters on FE Simulation of Ti6Al4V Orthogonal Cutting Process for Different Uncut Chip Thicknesses "2279

The constitutive model and its pertinent set of parameters are important input data in finite element modeling to define the behavior of Ti6Al4V during machining process. The present work focusses on comparing different constitutive models and the parameters sets available in literatures and investigating the quality of the predictions when varying uncut chip thickness (40 μm, 60 μm, 100 μm and 280 μm). In addition, temperature-dependent strain hardening factor along with strain softening phenomenon based reconstructed material model is proposed. The results from the numerical simulations are compared with experimental results available in literature. The comparison shows that the force values are highly influenced by constitutive models and the choice of parameters sets, whereas the chip morphologies are mainly influenced by the uncut chip thickness and constitutive models. This work justifies the need for an appropriate set of parameters and constitutive model that replicate the machining behavior of Ti6Al4V alloy for different cutting conditions.


Introduction
The titanium alloy, especially Ti6Al4V, is one of the most important and commonly used alloys for its attractive properties in aerospace, biomedical, submarine applications, etc. as it exhibits elevated strength, low density, good corrosion resistance, excellent high temperature properties [1,2]. Most of mechanical components with titanium alloys, are produced with machining operation [1,3]. The conventional machining process of Ti6Al4V has proved to be particularly complex due to its inherent properties such as, low thermal conductivity and high reactivity along with the diversity of physical phenomena involved, including large elasto-plastic deformation, complex contact/friction conditions, thermo-mechanical coupling and chip separation mechanisms [4][5][6]. Highspeed machining of Ti6Al4V alloy is still the interest of many researchers even after many decades of research [4]. The direct measurements of the stress, the strain, the strain rate and the temperature in the small deformation zone during chip formation are extremely difficult to assess experimentally [5]. So numerical techniques are often used. Finite element modeling is the most significant and extensively used numerical technique that is employed to investigate thermomechanical phenomena and nonlinearities involved in the simulation of machining process [6][7][8].
Since the early 1980s the numerical models are usually based on Lagrangian and Eulerian formulations [9,10]. In Eulerian formulation, the mesh is fixed in space and the material flows through the element that avoids critical element distortion. The main limitation of Eulerian model is that the chip geometry must be known in advance. The Lagrangian formulation has been extensively used by the researchers to simulate the chip geometry from the incipient stage to a steady state despite its drawback on mesh distortion problems [4,9,10]. These Lagrangian models require a separation criterion to remove the chip from the work piece. The separation criterions are based on either the geometrical considerations or the physical considerations. The geometrical method criterion is based on the chip separation along a predefined cutting plane and the physical method criterion is based on the chip separation using element deletion [11,12]. The Lagrangian models are also used with mesh adaptivity and automatic remeshing to eliminate the need for the chip separation criterion [13]. Consequently, this leads to an increase in the computation time [9,10].
More recent advanced numerical formulations such as the Arbitrary Lagrangian-Eulerian (ALE) formulation and Coupled Eulerian-Lagrangian (CEL) that take the best features of Lagrangian and Eulerian formulations have been implemented by Ducobu et al. [14,15] for modelling the orthogonal cutting process of Ti6Al4V alloy. In [15], the authors implemented Lagrangian, ALE and CEL formulations for different uncut chip thicknesses that produce continuous chips and showed that the results from ALE and CEL models are in good agreement with the experimental results with the reduction in computation time. The progression in the computational environments led the researchers to move from simple 2D to complex 3D finite element model of cutting process. Ducobu et al. [16] extensively discussed on 3D FE models that are implemented in the study of cutting process and developed a 3D finite element Coupled Eulerian-Lagrangian model that replicates the experimental cutting operation of Ti6Al4V alloy. For this study, an explicit Lagrangian finite element model with 2D assumption [7,17] is considered for the simulation of cutting process.
The mechanics of cutting process examines the deformation of a workpiece subjected to machining, such as chip formation, machining forces, residual stresses, temperature effects, strain rate effects and surface integrity [4]. A reliable constitutive material model is essential to relate to the large plastic strains at the very high strain rates and the very high temperatures observed during machining process [4,18,19]. So far, several empirical and physical based models have been developed and are employed by the researchers when machining of Ti6Al4V alloy, in which empirical models are mostly recommended for their flexibility, robustness and large availability of data's [4,20]. The Johnson-Cook (JC) [21] material model coupled with JC damage criterion [22] are highly employed for machining simulation. Mostly, the parameters identified for this law is by the experiments for testing the dynamic stress-strain response of materials that are fitted to stress strain curve [23]. These experiments can reach maximum strain of 0.5 and strain rate near to 10 3 s −1 which is well below the strain of 3 or even higher and the strain rate above 10 5 that are actually evident during cutting process. The resulting extrapolation leads to loss of some important information [24]. To solve this issue, many modified JC empirical models are proposed [25][26][27].
The parameters values that are implemented in constitutive models is the another most significant input that directly influence the simulated results. The availability of different sets of parameters in literatures raises the concern [28,29] which initiates the discussion on the influence of parameters sets and their impacts on the results. In [29], the authors justify the influence of the parameters set for a single cutting condition with three different constitutive model. In this present article the work is further extended with three other cutting condition along with the influence of four different constitutive models.
This paper investigates the finite element simulated results from 1.

2.
The Modified Johnson-Cook model from Calamaz et al. [25] that considers the strains softening phenomenon.

3.
The Modified Johnson-Cook model from Hou et al. [30] that takes into account temperature dependent hardening aspect and its coupled effects between strain and temperature.

4.
A reconstructed material model proposed in this study to figure out the influence of temperature hardening effect on the model proposed by Calamaz et al. [25]. The model combines strain softening phenomenon from Calamaz et al. [25] and temperature dependent hardening effect from Hou et al. [30].
These constitutive models are employed with three different sets of parameters acquired from Seo et al. [31], Calamaz et al. [25] and Hou et al. [30] Each constitutive model is analyzed with three different sets of parameters making a total of twelve models. These twelve models are investigated with four different uncut chip thicknesses (40 µm, 60 µm, 100 µm and 280 µm) while the rest of the cutting conditions are unaltered throughout the study. This work aims to highlight the effect on the results of the choice of the constitutive models and the parameters sets for uncut chip thicknesses ranging from small uncut chip thicknesses (h = 40 µm, 60 µm and 100 µm) that produce continuous chips to large uncut chip thickness (h = 280 µm) that produces serrated chip for the cutting speed of 30 m/min. This is necessary as it is currently still a challenge to predict the change in the chip morphology and to model the cutting forces without changing the parameters of the finite element model. The cutting force, feed force, chip morphologies and chip thickness measured from numerical simulations are compared and are validated with the experimental results that are acquired from Ducobu et al. [32].

Material Models
The empirical or phenomenological constitutive model is considered by many researchers for its mathematical simplicity and flexibility, which reflects the behavior of the material observed experimentally. These models consider the variables of deformations such as the plastic strain ε, plastic strain rate . ε and temperature T under macroscopic scale. The general form is:

Johnson-Cook Constitutive Model (JC)
The Johnson-Cook model [21] is a well-known empirical formulation which is heavily exploited in modeling and simulation of orthogonal machining. The JC model is widely employed by many researchers to determine the behavior of the material during machining process for its simple form and robustness. The Johnson-Cook flow stress equation links the plastic term, viscous term and softening term. Its flow stress equation is represented by the following equation: where A is the yield stress at reference temperature, B is the modules of strain hardening, n is the strain-hardening exponent, C the strain rate sensitivity and m is the thermal softening exponent. T is the current temperature, T melt and T room are the melting temperature and the room temperature, respectively, while . ε 0 is the reference strain rate. A, B, C, n, m are the constants that depend on the material and are determined by flow stress data obtained from material tests.
Even though the JC model is globally accepted, there are some noticeable drawbacks when considered for machining process simulations. The JC model describes the strain hardening phenomenon as an increasing function and it does not account for the softening behavior observed at the strains and temperatures in the primary shear zone regardless of the material. This leads to modified or updated versions of Johnson-Cook model to capture the unique behavior exhibited by certain alloys such as Ti6Al4V.

Modified Johnson-Cook Model by CALAMAZ (JC-Calamaz)
A modified version has been suggested by Calamaz et al. [25] taking into account the strain-softening phenomenon related to the recrystallization and dynamic recovery mechanism. The model includes novel multiplicative strain and temperature dependency by adding a hyperbolic tangent (TANH) term into the Johnson-Cook flow stress model. The material flow stress for TANH model is expressed by the following equation: where, D = 1 − (T/T melt ) d and S = [1 − T/T melt ) b , parameters A, B, C, n, m have the same meaning as for JC model and the new constants a, b, c and d are introduced by TANH law. The authors concluded that the TANH law significantly improves the prediction of chip morphology and suggested for the optimization of parameters for accurate prediction of the cutting force feed forces and the shear band frequencies.

Modified Johnson-Cook Model by HOU (JC-Hou)
A modified JC model is proposed by Hou et al. [30] by considering a temperature dependent strain hardening factor. JC model is altered by coupling temperature with the strain hardening factor in the strain function of JC model. This implies that when the temperature increases, the strain hardening rate decreases. The authors mentioned that the strain hardening rate Q of Ti6Al4V alloy has no noticeable strain rate sensitivity but has apparent temperature sensitivity. Strain or work hardening rate Q is a function of temperature: The flow stress equation is expressed by the following equation: where B is the hardening modulus, m 1 is the thermal sensitivity coefficient with the increasing strain and the rest of the parameters have the same meaning as for JC model. The author concluded that the strain-hardening rate is dependent on the temperature.

Modified TANH Material Model (JC-MTANH)
A reconstructed JC-Calamaz model is introduced in this paper. It combines the effect of temperature dependent hardening phenomenon proposed by Hou et al. [30] along with hyperbolic function introduced by Calamaz et al. [25] observed during high-speed machining process. The equation is constructed by introducing a temperature function into the work hardening term (proposed by Hou et al. [30]) along with flow softening at higher strain (proposed by Calamaz et al. [25]). The flow stress equation is expressed by the following equation: where, m 1 is the thermal sensitivity coefficient with the same meaning as for JC-Hou and The rest of all other parameters meaning is the same as for JC model.

Constitutive Models Parameters for Modeling of Machining Process of Ti6Al4V
The behavior of the material during machining process is determined by the constitutive model and its parameters. The values of the parameters are significantly important in predicting the results. The parameters for the models are usually identified by curve fitting of experimental stress-strain plot, this leads to some discrepancies in the parameters identification. In fact, there are plenty of parameters sets available in the literatures for the same material. This leads to the discussion on identification of right parameters set to simulate the results.
In this context, Ducobu et al. [28] extensively analyzed 20 different sets of Johnson-Cook parameters available in the literature for Ti6Al4V by comparing the outputs of a Coupled Eulerian-Lagrangian orthogonal cutting model with the experimental results. They concluded that the parameters set from Seo et al. [31] is well capable of producing near accurate results for the forces and chip thickness. From the above discussion, the set of parameters for JC model is adopted from Seo et al. [31], set of parameters for JC-Hou comes from Hou et al. [30] and the set of parameters for JC-Calamaz is adopted from the work of Calamaz et al. [25]. The sets of parameters are mentioned as 1, 2 and 3, respectively, throughout the article and their values are reported in Table 1. Table 1. Sets of parameters considered for this study-set 1 [31], set 2 [30] and set 3 [25].
The values of m 1 , a, b, c, d identified by Hou et al. [30] and Calamaz et al. [25] are kept constant for all models if the constitutive equation uses these parameters. Only the common parameters are such as A, B, C, n, m T melt and T room are interchanged and are represented correspondingly as 1, 2 and 3. In order to highlight the influence of parameters sets on numerical modeling of cutting process simulation the three different sets of parameters are considered with the four constitutive models.
The evolution of stress-strain curve is plotted at fixed temperatures of T = 573 K and T = 973 K for strain rate of From the evolution of stress-strain curves in Figures 1 and 2. It is evident that at high strain and high temperature the JC-Calamaz law takes the softening behavior into account and the JC-Hou law that considers the temperature dependent hardening effect shows a significant difference in the evolution of flow stress at higher temperature and strain rate. The difference in the stress levels might influence the result. Indeed, it is essential to note that the initial stress value of JC-Calamaz_3 and JC-MTANH_3 is very low when compared with initial stress value of the other two constitutive models considered in this study and this could significantly influence the numerical results. From the stress strain curves, it is seen that the influence of introducing temperature dependent strain hardening term in JC-Calamaz law is less significant, but this is further discussed with numerical simulations.  The stress-strain evolutions are observed at two different temperatures. When temperature increases from 573 K to 973 K the stress level decreases as expected, confirming the presence of thermal softening phenomenon. Figure 2 shows the flow stress evolution curves for the four constitutive models with three different sets of parameters, respectively. When compared, the stress-strain evolution with set 3 parameters shows a low initial stress value regardless of the constitutive laws, while the other sets of parameters show approximately the same initial stress value. However, JC-Hou_1 shows higher stress value with respect to increasing strain. This confirms the influence of temperature dependency in strain hardening.
The stress strain evolution curves confirm the influence of the choice of the sets of parameters for each constitutive model and their capability in influencing the results.

Experimental Reference
The Experimental reference is considered from the work performed by Ducobu et al. [32] on orthogonal cutting of Ti6Al4V on a high-speed milling machine shown in Figure 3.  Table 2. The chips were observed with an optical microscope. Globally, continuous chips were observed for the uncut chip thicknesses of 40 µm, 60 µm and 100 µm even if very small and irregular teeth along its entire length are noticed. For uncut chip thickness of 280 µm saw toothed or serrated chips are observed with similar teeth. The chip morphology for 280 µm with a magnification factor 50 and for 40 µm, 60 µm and 100 µm with a magnification factor of 200 is shown in Figure 4. The cutting forces are measured in the three directions with a Kistler 9257B dynamometer. The RMS cutting force (Fc), feed force (Ff) and cut chip thickness (h ) measured from the experimental results are given in Table 3. The chip morphology of h = 40 µm, 60 µm and 100 µm is considered as continuous chips. Whereas the chip morphology of h = 280 µm is clearly a saw-toothed chip. For h = 280 µm the undeformed tooth length, L, the tooth height, H and the valley, C, are calculated and are given in Table 4.

Finite Element Model
In this finite element model developed with Abaqus, the work piece is fixed in space and the tool moves over the workpiece with a prescribed velocity. A two-dimensional thermo-mechanical model with plain strain assumption is adopted for this work. An explicit algorithm with Lagrangian finite element formulation is considered for the simulation of orthogonal cutting process. The workpiece is modeled as an elastoplastic material, with the dimension of a rectangular block of 1.5 mm by 1 mm. The regions of interest in the work piece are modeled with an element size of 5 µm by 5 µm (that agrees with the size of grains) with square linear quadrilateral elements of type CPE4RT. The material properties adopted for the workpiece are given in Table 5. The interaction between tool and workpiece is modeled as Coulomb friction law with a friction coefficient of 0.2 [33]. Tungsten carbide is selected as tool material and the cutting tool is defined as an elastic material with the properties given in Table 5. The cutting conditions and tool geometry are the same as in the experiments and are given in Table 2. The tool and the workpiece basic geometry along with the mesh are kept constant for all the finite element models. Only the uncut chip thickness and the associated fracture energy regions dimensions are altered accordingly.
The boundary conditions and the basic geometry for the uncut chip thickness h = 100 µm are illustrated in Figure 5. The predefined temperature field for both tool and workpiece is defined according to the parameters set chosen for simulation. For this study, an energy-based failure criterion is adopted for FE modeling of orthogonal cutting simulation using Lagrangian formulation [37]. For chip separation criteria, a two-stage process is considered: damage or fracture initiation and fracture evolution. For damage initiation, Johnson-Cook damage model [23] is considered.
where D is an internal state variable which links the ratio between the incremental plastic strain (∆ε) and the effective plastic strain to fracture ε f . D is a monotonically increasing function and the damage initiation occurs when D reaches the values of 1. The term ε f , takes the form.
The material parameters d 1 , d 2 , d 3 indicates the fracture strain dependency to the triaxial stress, d 4 the strain-rate dependency and d 5 accounts for temperature softening these values are obtained through mechanical properties experiments. The JC failure model parameters adopted for Ti6Al4V are reported in Table 6. An energy-based fracture criterion is applied to the FE model for failure evolution. The failure evolution law is specified in terms of fracture energy dissipation G f [38]. G f is defined as the energy required to open a unit area of crack. After damage initiation instead of a stress-strain response, a stress-displacement response is followed. The fracture energy is defined as: where L is the characteristic length of the element (it depends on the element geometry) and u is the equivalent plastic displacement at failure. Before damage initiation, u = 0 and after damage initiation u = Lε [39]. From the concept of fracture mechanics, it is considered that during orthogonal cutting process two different values of fracture energies are considered as an input in Abaqus explicit code. In orthogonal cutting, the crack propagation is a combination of Mode I (opening mode) and Mode II (shearing mode). The Mode II will occur along shear plane accompanied by considerable shear deformation that are evident in region 1 and are given by fracture energy (G f ) II . The (G f ) I denotes the fracture energy of Mode I which is a tensile mode (opening mode normal to the plane of the fracture). The layer of chip is defined by region 1 and the tool passage layer is defined by region 2 [35,37,40], as shown in Figure 5. The thickness of the regions is given by the equation: The fracture energy is given by Equation (9) [27]: where K c is fracture toughness and the value of K Ic = 55 MPa √ m [41], E is the Young's modulus and ν is Poisson's ratio. The evolution of damage D f in chosen to be an expo-nential. When D = 1 and D f = 0, damage is initiated, and the material fails when D f reaches 1.
In this paper a total of 48 finite element simulations has been performed by considering three sets of parameters for four different uncut chip thickness with each of the four constitutive models.

Numerical Results and Discussions
Initially, JC_1, JC-Hou_2, JC-Calalamz_3 and JC-MTANH_3 (i.e., the numerical simulation of JC-MTANH_3 is performed with set of parameters identified by Calamaz et al. [25] with m 1 from Hou et al. [30]) for uncut chip thickness of 60 µm are analyzed and compared with the experimental reference. Further they are investigated with different parameters sets of parameters with the three other uncut chip thickness (40 µm, 100 µm and 280 µm) considered in this study.

Comparison of Constitutive Models for Uncut Chip Thickness of 60 µm
The temperature distribution of JC_1, JC-Hou_2, JC-Calamaz_3 and JC-MTANH_3 chips are shown in Figure 6a-d. Heat transfer in the machining process takes place primarily in the shear zones, where the plastic work is converted into heat and the chip-tool interface where frictional heat is generated [36]. The maximum temperature observed in the secondary deformation zones is expected. Indeed, differences in temperature are seen between the simulated chips from these models. The temperature in the secondary deformation zone of chips produced by JC_1 and JC-HOU_2 is quite close and higher when compared with JC-Calamaz_3 and JC-MTANH_3. That can be explained by the level of stresses shown in Figures 1 and 2. The temperature in the numerically simulated chip from JC-MTANH_3 model is notably high when compared with JC-Calamaz_3 model. It is also important to note that in JC-Calamaz model the temperature distribution is somewhat uniform throughout the chip. In addition, highly deformed and elongated elements are observed in JC-Calamaz_3 and JC-MTANH_3. The contour plots of equivalent plastic strain for JC_1, JC-Hou_2, JC-Calalamz_3 and JC-MTANH_3 is shown in Figure 7a-d. In JC_1 and JC-Hou_2 models, high plastic strain is located in the chip side in contact with the tool rake face.  Table 7.
To conclude on the results of the comparison on the four phenomenological constitutive models for the above cutting condition, the predictive finite element model with JC_1 and JC-Hou_2 models can accurately predict forces, chip morphology and chip thickness. The results from JC-MTANH model are nearly the same as from the JC-Calamaz model, even though JC-MTANH model includes temperature dependent strain hardening factor. Number*-the difference is in the uncertainty range of the measurement.

Influence of Parameters Sets for Different Uncut Chip Thickness
This part will focus on the combined influence of parameters sets and constitutive models on the results, which is further extended to the different uncut chip thickness considered in this study. It is important to note that for Ti6Al4V there is large availability of parameters sets in the literature. The constitutive models are fitted even on the physical parameters which can be seen as problematic. Determining the most suitable parameters sets is always a concerning point in numerical simulation of machining process. In this context three sets of parameters are considered as given in Table 1. To validate the influence of the choice of the parameters and effect of uncut chip thickness on constitutive models, only parameters sets are interchanged in the constitutive models for 4 uncut chip thickness (i.e., 4 models are employed with 3 sets of parameters for 4 different uncut chip thickness making a total of 48 (4 × 3 × 4) simulations).

Forces
Cutting and feed forces are generated using these models for the uncut chip thickness of 40 µm. The RMS value of cutting force calculated from JC and JC-Hou with set 1 and set 2 parameters are quite close around 86 N/mm, but with set 3 parameters significant differences are observed. This can be explained by the fact that the set 3 parameters are identified by Calamaz et al. [17] for a particular cutting condition. Comparing with the two other constitutive models, JC-Calamaz and JC-MTANH, with any set of parameters predicts cutting force value well below the experimentally calculated values and includes oscillation in the force evolution. This can be explained by the presence of hyperbolic function in JC-Calamaz and JC-MTANH models.
The RMS value of feed force calculated from the JC and JC-Hou models with set 1 and set 2 parameters are in variation of +9 N/mm to +25 N/mm when compared with experimental result. The feed force calculated from JC-Calamaz_3 and JC-MTANH_3 is higher when compared with JC and JC-Hou. The RMS values cutting force and feed force are given in Table 8. Number*-the difference is in the uncertainty range of the measurement.
When the uncut chip thickness increases, the RMS value of cutting force increases for all the constitutive models with any sets of parameters in the same pattern as for the uncut chip thickness of 40 µm. Indeed, the RMS values of cutting force for JC and JC-Hou with set 1 and set 2 parameters are quite close, but with set 3 parameters differences are quite high when compared with experiments. The JC-Calamaz and JC-MTANH with any set of parameters give lower cutting force values when compared with the two other models. The RMS values of forces for h = 100 µm and 280 µm are given in Tables 9 and 10.  Number*-the difference is in the uncertainty range of the measurement. * these chips are fully developed serrated chips and the corresponding measurements are given in Table 11. Table 11. Numerical saw-toothed chip measurement of tooth length, L, the tooth height, H, the valley, C and the standard deviation of the length x for h = 280 µm, σ x . Experiment  206  17  288  14  157  21  JC-Calamaz_1  171  12  316  16  108  26  JC-Calamaz_2  188  16  302  26  82  17  JC-Calamaz_3  168  9  288  5  104  18  JC-MTANH_1  178  11  312  12  112  21  JC-MTANH_2  180  19  306  28  101  8  JC-MTANH_3  166  7  287  7  96  16 The RMS values of feed force for uncut chip thickness of 40 µm and 60 µm is overestimated by all the four models with any set of parameters. However, for uncut chip thickness of 100 µm and 280 µm the RMS feed force value from JC and JC-Hou with any set of parameters is underestimated. On the contrary, JC-Calamaz and JC-MTANH with set 3 parameters overestimate the value of feed force for uncut chip thickness of 100 µm and 280 µm.
From this investigation, the force values generated by set 2 parameters with the four constitutive models show improvement towards the prediction of cutting force and feed force near to the experimental value.

Predicted Chip Morphology
The chip morphologies for uncut chip thickness of 40 µm, 60 µm, 100 µm and 280 µm are analyzed in this section. Continuous chips are produced by JC and JC-Hou for any uncut chip thickness and with any sets of parameters. Whereas the JC-Calamaz and JC-MTANH model with any set of parameters for uncut chip thickness of 40 µm, 60 µm and 100 µm produces slightly serrated chip. This serration becomes more pronounced with increase of uncut chip thickness to 280 µm.
Initially the temperature contours of uncut chip thickness of 280 µm are analyzed and are plotted in Figure 8. The temperature rise in the secondary shear zone is evident for all models with any set of parameters. When the uncut chip thickness increases the temperature in the secondary deformation zone also increases. It is also observed that the large increase of temperature in the primary shear zone leading to the formation of the adiabatic shear band is clearly noticed for 280 µm with JC-Calamaz and JC-MTANH. In addition to above discussion the evolution of cutting force clearly depicts the oscillation in the curve which confirms the presence of serrated chip. The JC and JC-Hou models produce continuous chip with maximum temperature value localized along the tool-chip interface with steady state evolution of forces. It is also observed that temperature is higher and quite close for JC and JC-Hou models with parameters sets 1 and 2. But with set 3 parameters the temperature is comparatively lower. The JC-Calamaz and JC-MTANH models follow the same trend: with set 1 and set 2 parameters, the temperature is high, and with set 3 parameters, the temperature is low. Nevertheless, the temperature is notably higher for JC and JC-Hou models regardless of the parameters sets.
For the three other cutting conditions (h = 40 µm, 60 µm and 100 µm), the temperature is higher and quite close for JC and JC-Hou models with set 1 and set 2 parameters, whereas with set 3 parameters, the temperature is less. This follows the same trend as for h = 280 µm. Continuous chips are observed in all the cases, and it is obvious that temperature rises mainly in the deformation zone of chip region, with maximum values localized along the tool-chip interface.
The JC-Calamaz and JC-MTANH model for h = 40 µm and 60 µm the temperature distribution is relatively homogenous in the chip for any set of parameters and the maximum temperature is not clearly observed in the secondary deformation zone. When the uncut chip thickness increases to 100 µm, the temperature variations are observed with maximum temperature in the secondary deformation zone. This might be due to the incorporation of hyperbolic tangent function in the behavior law that alters the stress level and influences the temperature.
The adiabatic shear bands are not clearly noticed for this cutting condition. Nevertheless, the chips produced by these two constitutive models are slightly serrated and the serration is further confirmed by investigating the equivalent plastic strain contour. It is also worth noting that JC and JC-Hou models produce a continuous chip with less distorted or elongated elements in the chip region whereas with the JC-Calamaz and JC-MTANH models, highly elongated and distorted elements are clearly visible. This could possibly influence the results to some extent.
Overall, when the uncut chip thickness increases, the temperature in the secondary deformation zone also increases. This is because when h, increases the tool experiences more cutting resistance resulting in an increased temperature. Some improvements in the forces and temperature values are seen from the results produced by all the models with set 2 parameters. The temperature contours of h = 40 µm, 60 µm and 100 µm with set 2 parameters for the four constitutive models are given in Figure 9. Equivalent plastic strain contours are analyzed to clearly understand the chip morphologies. For continuous chips, the equivalent plastic strain (PEEQ) lies in the range of 1 to 3 and for serrated chip the equivalent plastic strain is approximately in the range of 1-5. However, in some cases PEEQ reaches the value of 10 or even more where elements are highly distorted or elongated.
For all the continuous chips produced by the JC and JC-Hou models, high plastic strain is observed at the chip side in contact with the tool rake face. For slightly serrated chips from JC-Calamaz and JC-MTANH two unequal segments with high strain localized areas and areas with low deformation is observed. This further confirms the presence of serration in the chips.
The saw tooth is clearly observed for h = 280 µm and are presented in Figure 10. The chip thickness values of the continuous chips from JC and JC-Hou for h = 40 µm, 60 µm and 100 µm with any set of parameters are within the deviation range of 1% to 10%. In some cases, a thinning of chips is observed with the increase of feed force.   Table 11.
The parameters used to characterize chip serration are given in Figure 10. The length of the tooth is represented by L. The maximum tooth height or peak and the valley from the base of the chip are represented by H and C, respectively.

Discussion
For uncut chip thickness h = 40 µm, 60 µm and 100 µm, the RMS values of cutting force from JC and JC-Hou models with set 1 and set 2 parameters are very close to the experimental values within the range of ±5% and feed forces are in the range of ±20%. With set 3 parameters, the cutting forces values are notable less with a difference of 15% when compared with experimental reference, whereas the feed force are in the same range. The JC-Calamaz and JC-MTANH models with any set of parameters underestimate the cutting force value (with a difference of 25%) and overestimate feed force value (in the deviation range of 25%).
For the case of h = 280 µm, the trend is somewhat different. The JC and JC-Hou models with sets 1 and 2 parameters cutting force are in the same range at ±5%, but underestimate feed force (with a difference of 50%). The cutting and feed force values are underestimated (with a difference of 35%) by JC-Calamaz and JC-MTANH models with set 1 and set 2 parameters when compared with experimental results. Based on the above discussions, the cutting force values predicted by JC-Calamaz and JC-MTANH models with any sets of parameters are very low when compared with experimental force values. Indeed, the metallurgical state and operating conditions of Ti6Al4V considered by Calamaz et al. [25] is unknown and may be different from the Ti6Al4V employed in experiments. This might influence the force calculation.
The temperature in the secondary deformation zone increases while uncut chip thickness increases, and a direct link was observed between temperature and cutting force value for all the cases. The comparison of simulated chips with the experimental chips shows that JC and JC-Hou constitutive models with any sets of the parameters for h = 40 µm, 60 µm and 100 µm produce continuous chip as the experimental reference.
Considering uncut chip thickness h = 280 µm, the chip morphology observed from the experiments is a saw-toothed chip, but the JC and JC-Hou models with any set of parameters give continuous chip. This is because in JC and JC-Hou law flow stress increases with increasing strain regardless of temperature. In addition, Guo et al. [42] demonstrated that these models lack in predicting the adiabatic phenomenon responsible for serrated chip formation observed during machining of these alloys.
On the contrary, the JC-Calamaz and JC-MTANH models produce slightly serrated chips for all cutting conditions and with any set of parameters considered in this study. This serration of chip is more prominent with increase of uncut chip thickness.
From the above discussions, JC-Hou with set 2 parameters predicts cutting force and feed force values very close to the experiments for all uncut chip thicknesses considered in this study. This model can predict the chip morphologies for small uncut chip thickness i.e., h = 40 µm, 60 µm and 100 µm but, fails to predict the chip morphology for h = 280 µm. This is because this model does not include an important behavior of Ti6Al4V when machined: the adiabatic and strain softening phenomenon observed. Even though the JC-Calamaz and JC-MTANH models include strain softening phenomenon and temperature dependent strain hardening factor, the chips produced by these models are all serrated for any uncut chip thickness, although this is not the case with the experimental chips.
The experimental results from Ducobu et al. [31] highlight that the chip morphology at a constant cutting speed of 30 m/min is influenced by the uncut chip thickness. The increase in chip thickness leads to the accumulation of heat in the primary shear zone which promotes the serrated chip. The FE simulation and constitutive models fails to predict the transition from continuous chip to serrated chip for any set of parameters and the cutting conditions considered in this study.
The JC-MTANH model proposed in this work which incorporates temperature dependent strain hardening and hyperbolic tangent law together does not show a notable improvement in predicting the forces and chip thickness values. Meanwhile, JC law with temperature dependent hardening factor can significantly improves the prediction of forces and chip thickness.
Overall, the main objective of this work is to employ the exact same FE model for different cutting condition to predict forces and chip morphologies. These models work well for some cutting conditions and fails for other cutting condition. Consequently, for the cutting conditions considered in this work, all the models examined in this study fail to capture the transition of continuous chip to serrated chip. The transition from continuous chip to serrated chip and the influence of strain softening phenomenon and their temperature dependency is still a valid question. In addition, identification of the right parameters sets is of prime importance in predicting the forces and chip thickness through numerical simulation of cutting process.

Conclusions
In this research work, comparison of constitutive models and the influence of parameters sets on different uncut chip thickness was performed through finite element analysis of orthogonal cutting of Ti6Al4V. In the meantime, a reconstructed constitutive model is proposed. The forces and chip morphologies from FE simulations are compared with experimental references to highlight the differences. The following conclusions are drawn from this study: • A direct link is observed between temperature and cutting force values. Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.