INCONEL ® Alloy Machining and Tool Wear Finite Element Analysis Assessment: An Extended Review

: Machining INCONEL ® presents significant challenges in predicting its behaviour, and a comprehensive experimental assessment of its machinability is costly and unsustainable. Design of Experiments (DOE) can be conducted non-destructively through Finite Element Analysis (FEA). However, it is crucial to ascertain whether numerical and constitutive models can accurately predict INCONEL ® machining. Therefore, a comprehensive review of FEA machining strategies is presented to systematically summarise and analyse the advancements in INCONEL ® milling, turning, and drilling simulations through FEA from 2013 to 2023. Additionally, non-conventional manufacturing simulations are addressed. This review highlights the most recent modelling digital solutions, prospects, and limitations that researchers have proposed when tackling INCONEL ® FEA machining. The genesis of this paper is owed to articles and books from diverse sources. Conducting simulations of INCONEL ® machining through FEA can significantly enhance experimental analyses with the proper choice of damage and failure criteria. This approach not only enables a more precise calibration of parameters but also improves temperature ( T ) prediction during the machining process, accurate Tool Wear (TW) quantity and typology forecasts, and accurate surface quality assessment by evaluating Surface Roughness (SR) and the surface stress state. Additionally, it aids in making informed choices regarding the potential use of tool coatings.


Introduction 1.Material Numerical Modelling
INCONEL ® alloys are hard-to-cut metals that pose machinability challenges yet to be overcome [1].An application upsurge for these Ni-superalloys has been observed across several industries, including aeronautics, automotive, and energy power plants [2].Renowned for their ability to withstand high temperatures (T) without succumbing to creep, these superalloys proved to be highly attractive and well suited for producing items like jet engines and steam turbines.Following the work of Pedroso et al. [2] and Thornton et al. [3], it becomes evident that understanding the chip formation mechanism is one of the crucial factors that will lead to a comprehensive understanding of the entire cutting process in conventional manufacturing (CM) [4,5].When a material exceeds its yield point (σ y , yield stress), it experiences a strain-hardening phenomenon, leading to increased load and energy [6] required for material deformation [7].Two hardening models address elastoplastic scenarios: (1) isotropic hardening, which refers to a material's ability to uniformly increase its yield strength and stiffness under repeated loading without changing the shape or orientation of its yield surface [8], and (2) kinematic hardening, which, in turn, involves changes in the material's yield surface shape and orientation, but not its size, in response to plastic deformation [9].Understanding the depiction of material behaviour is paramount, particularly concerning how strain (ε), the strain rate ( .ε), and T impact the material's flow stress (σ).Table 1 displays several integrated material constitutive models in use.For a better grasp of the additional variables outlined in Table 1, it is recommended to refer to the research by Iturbe et al. [10], Lewis et al. [11], Lin et al. [12], and Rudnytskyj et al. [13].
In FEA simulations [19][20][21][22], INCONEL ® alloys are typically represented as materials exhibiting elastoplastic behaviour and isotropic hardening, with σ influenced by ε, .ε, and T [23].The prevalent analytical model used to reproduce the elastoplastic behaviour of INCONEL ® alloys is the JC constitutive model [14][15][16][17][18].The JC model has some derivations: (1) the Baumann-Chiesa-Johnson (BCJ) model [24], which incorporates .ε and T sensitivity, as well as damage, through a yield surface [25], and (2) the Steinberg-Cochran-Guinan (SCG) model [26], which describes the shear modulus (G) and σ y at high .ε. Due to the inherent flaws of the JC model, as addressed in this section, some improvements have been developed, such as the derivations in [25][26][27].The JC model is deemed a rate-independent constitutive model [28], along with the other models detailed in Table 1, and the expression for the stress in the JC constitutive model (σ JC ) is provided in Equation (1): where ε p is the equivalent strain; . ε p is the equivalent strain rate; . ε 0 p is the reference strain rate; T 0 is the ambient temperature; and T m is the melting temperature.The yield stress (A), strain-hardening constant (B), modulus of strain rate hardening (C), strain-hardening coefficient (n), and thermal softening coefficient (m) [29,30] are material constants and exponents derived by fitting the data obtained from the tensile test, Split Hopkinson Pressure Bar (SHPB) test [30,31] (not following any agreed-upon standards regarding the operation or design of the apparatus), and dilatometer test under various T and .ε. Equation (2) provides insight into the linear relationship between ln(σ JC − A) and ln(ε p ).In quasistatic conditions, viscous-plastic deformation and thermal softening hold little significance, resulting in this trivial expression: Figure 1 illustrates the JC constitutive model equation segregated into three phenomena, each associated with their specific experimental tests and conditions.In order to establish the A, B, C, m, and n constants of the JC model, an initial tensile test is conducted following ASTM A 370-17 [32] guidelines.As shown in Equation (3), the thermal softening effect is deemed negligible for determining C, leading to a rearranged expression.The value of C is derived through a linear regression model using different .ε values obtained from the SHPB test [33].
The JC model can be reorganised to determine the m parameter, overlooking .ε strengthening effects, as demonstrated in Equation (4) [33].Estimating m involves performing the push-rod dilatometer test under ASTM E 228-17 [34], analysing the linear thermal expansion, and evaluating m in Equation (4) using the linear regression model graph [35].
The JC material model remains extensively employed in metal-cutting simulations [36].Equation (1) describes strain-hardening (B × ε p n ) as an ascending function [37] that tends towards infinity and does not account for strain-softening effects [10].Softening phenomena are crucial for instigating and amplifying the chip plastic state by forming adiabatic shear bands during machining operation simulations.In recent times, parameter determination has seen advancements, notably by Škrlec and Klemenc [39], who devised optimisation algorithms to streamline the determination of material constants for E185 structural steel.Murugesan and Jung [15] formulated surrogate models based on constitutive relations, and the model constants were estimated using experimental data.Souza, et al. [33] successfully predicted the behaviour of INCONEL ® 625 by integrating the JC and Avrami models [40], as demonstrated in Equation ( 5), for the resulting stress (s JC + A ): where σ JC saturation stress (σ s ), the fraction of microstructure (X, 0 ≤ X ≤ 1) that underwent dynamic recrystallisation (DRX) [41], and the steady-state stress (σ ss ) are acknowledged.
The FEA-predicted behaviour surpassed the existing data in the material database.Lewis et al. [11] aimed to establish suitable constitutive models for simulating the mechanical behaviour of INCONEL ® 625 components fabricated using Laser Powder Bed Fusion (LPBF) for aerospace applications.SHPB testing results were used to calculate the coefficients for five constitutive models: JC, modified JC, Hensel-Spittel, modified Hensel-Spittel, and modified Zerilli-Armstrong models.It was determined that the last model referred to has the lowest absolute relative error values of the models used: 2.88% for as-printed and 2.71% for annealed specimens.Hokka et al. [42] studied the mechanical behaviour of Ti-6246 and an INCONEL ® 625-like alloy while obtaining JC material model parameters from the experimental data.The model was used to describe the plastic behaviour of the studied alloys in simulations of the orthogonal cutting of the material.The JC model for the INCONEL ® 625-like alloy was improved by introducing an additional strain-softening term (∆T in Equations ( 6) and ( 7)) that allowed a decrease in .ε hardening at large deformations.
For better comprehension of the Equation ( 7) variables, it is suggested to consult Hokka et al. [42].Table 2 provides the A, B, C, m, and n values for INCONEL ® alloys.Notably, information is abundant regarding INCONEL ® 718, while there is a dearth of information about the INCONEL ® 625 alloy.Heat treatment is a variable that significantly influences JC model parameters [10]; hence, substantial variations can be observed between research groups, particularly for INCONEL ® 718.INCONEL ® machinability modelling in FEA must rely on a proper damage evolution model.Since INCONEL ® 718 and 625 present elastoplastic behaviour, as shown in Figure 2, from T 0 through different T values until reaching T ≈ 1200 • C, a ductile material damage model is the most suitable [52][53][54][55][56][57][58] to predict the damage caused by machining.
In ductile materials, the ultimate failure unfolds through three distinct phases.While true-stress curves are attractive to study materials in some ways, FEA software works essentially with engineering stress curves, and in Figure 3, a typical uniaxial engineering σ-ε response for a ductile material is illustrated.Phase one is the linear behaviour of the elastic regime deformation that ends in B; in phase two, the material transitions into plastic deformation (curve BCDE), where C, or equivalent plastic strain at damage initiation (ε p 0 ), is the point of damage (D) initiation in the plastic regime (D = 0, and cumulative damage parameter ω = 1).D is an intermediate point of damage evolution, and point E, a boundary between the plastic regime (phase two) and complete fracture (phase three), indicates the point of ε p f , and the straight line d' depicts ε hardening [60].The principle of damage classification in a plasticised material is the characterisation of a degraded Young's modulus (E * in Equation ( 8)), which is a function of the pristine Young's modulus (E), and degraded material flow stress, or mechanical strength (σ * in Equation ( 9)), which is calculated using the intact material flow stress (σ) [55].
The JC shear fracture damage model [61,62] can forecast the D sustained by tracking the increment in equivalent plastic strain (∆ε p ) and equivalent plastic strain at fracture (ε p f ).Equation (10) demonstrates the ω calculus, which is updated in every FEA solving step.
The material element removal criterion embedded within FEA programs, such as ABAQUS ™ or ANSYS ® , is the trivial relation between = ε p f and equivalent plastic strain (ε p ) [62], as described in Equation (11), (11) and in Equation (12), which shows the calculus of ε p f [16,63,64]: where σ m is the average stress, and D 1 and D 3 depict how stress triaxiality influences the initial failure ε and the impact of the equivalent plastic strain rate ( .ε p ). D 2 is an exponential factor of JC, and the thermal effects are denoted by D 4 and D 5 [65].Table 3   The employment of the JC shear fracture damage model in isolation corresponds to an instance of abrupt failure (Figure 3, curve described by ABCε p 0 ), not capturing D evolution.Instead, when combined with the D evolution approach, the JC shear fracture damage model is exclusively employed as a criterion to represent the onset of D. In contrast, an evolution criterion controls the behaviour of the material through curve CE (Figure 3).Table 4 compiles some of the fracture energy-based criteria that characterise the propagation of damage in a crack [60].Cohesive Zone Model (CZM) and R-curve approach are also mentionable fracture-energy-based criteria.
These fracture-energy-based criteria have been widely used for machining simulations, particularly Hillerborg's fracture energy criterion [54][55][56].The model's second form of the equation is solved in order to obtain u p , a change of variable that is very useful so that u p f can be calculated according to Equation (13): where σ y0 stands for yield strength at the damage initiation state (point C in Figure 3).The D evolution depends on the rate of stiffness degradation ( .D, from point C to point E in Figure 3) and can be expressed either linearly (Equation ( 14)) or exponentially (Equation (15)): .
Table 4. Fracture-energy-based criteria for damage progression.

Model Equation
Griffith's critical stress intensity factor [57] -displacement gradient vector; a-crack dimension; A f -fracture area; G c -critical energy release rate; G f -fracture dissipation energy; J-energy release rate; K IC -critical stress intensity factor; L-characteristic length; n-normal vector to the boundary; t-coordinate directions; u p -equivalent plastic displacement; u As an example of the applicability and culmination of material damage models, Liu et al. [60] investigated the residual stress evolution in INCONEL ® 718 affected by multiple cutting operations, making use of the JC shear fracture damage model alongside the Hillerborg's fracture energy criterion to predict the serrated chip morphology, as depicted in Figure 4.The resulting FEA simulations showed that the residual stress level might be controlled by optimising the previous cuts to obtain the desired surface integrity.In addition to material modelling, damage initiation, and propagation characterisation, it is essential to model tool wear (TW) to predict tool life (TL).While some authors may treat the tool as a rigid body in FEA when their primary aim is to study INCONEL ® 's chip morphology, this review also aims to shed light on the existing mathematical models for TW and findings on improving TL.

Tool Numerical Modelling
TW significantly impacts both quality and productivity, particularly in the case of INCONEL ® alloys like 718 and 625, whether considering coated tools [66], ceramics, or even PcBN.Sousa et al. [67] scrutinised recent advancements in TiAlN-based coatings, including nanolayered, nanocomposite, and Ru-, Mo-, and Ta-doped coatings, evaluating their mechanical properties and comparing their cutting behaviours during turning and milling processes.Later on, Sousa et al. [68] experimentally investigated the wear characteristics of multilayered TiN/TiAlN-coated end mills through Physical Vapour Deposition (PVD) [69] and High-Power Impulse Magnetron Sputtering (HiPIMS).These tools were employed in finishing operations on INCONEL ® 718 to advance the comprehension of the wear patterns exhibited by coated tools during the machining of these alloys.Silva et al. [70] evaluated machined surfaces' integrity and TW resistance using cutting tools coated through PVD HiPIMS with TiAlYN during the end milling of INCONEL ® 718.These efforts demonstrate that experimental work is time-consuming and simulation can save resources and time for researchers.Three-dimensional (3D) modelling is essential to analyse and predict flank wear (VB), and FEA enables predictive insights into the sequence of TW mechanisms [2].The ISO 8688-2:1989(E) [71] standard defines a milling tool as worn when a uniform VB = 300 µm is achieved or the local maximum flank wear (VB max ) is 500 µm.The ISO 3685:1993(E) [72] standard provides relevant turning tool guidelines, considering the same VB limits.Table 5, summarised by Pedroso, et al. [2], presents typical mathematical models characterising various TW mechanisms.For a deeper understanding of the additional variables outlined, it is advisable to consult Wang et al. [73].

TW Model Comments
Taylor Taylor's empirical TL model.
Surface wear model.
Halila [84,85] TW model is dependent on the material removal rate.
Pálmai [86] dW TW model, considering the effects of wear-induced cutting, force, and T rise on TW.
The JC model, among those presented in Table 1, has enhanced the ability to predict the behaviour of INCONEL ® alloys during machining, and presently, the academic community can predict elastoplastic phenomena [89]; nonetheless, TW behaviour predictability while machining is an ongoing challenge in the academic community.Computationally, the authors must pre-process TW via Computer-Aided Design (CAD) according to experimentally obtained geometries.A simulation of the entire TL is unreasonable due to the long computing time (t) taken [90].Figure 5 presents a schematic of the actual procedure used by many authors to simulate machining when considering TW.From the point of view of studying the T distribution, considering TW is paramount because cutting becomes less effective with each machining pass, and therefore, T rises as the tool deteriorates.Although it is not an explicit simulation where material machining and TW coexist, speaking in computational terms, a quasi-transient situation can still be examined.Thus, the continuous progress of TW is discretised into finite steps for simulation [90].
For instance, Liang et al. [92] studied the experimental orthogonal cutting process of Ti-6Al-4V, and VB was considered the only tool failure standard despite both the crater wear and VB appearing simultaneously.Afterwards, a numerical model was built with three distinct tool CAD models: (1) new, VB = 0 mm; (2) semi-worn, VB = 150 µm; and (3) worn, VB = 300 µm; these were used to determine initial, regular, and failure stages.Both simulated and experimental findings revealed that the evolution of local plastic behaviour was primarily driven by the thermo-mechanical loads resulting from TW effects.Shi et al. [37] carried out a study employing Usui's model [46,76,77], analysing output variables like the tool-chip contact area (w), pressure, sliding, velocity, and initial T. Hosseinkhani and Ng [90] developed a new empirical model (Equation ( 16)) to predict TW in the initial or break-in period as a function of the Von Mises stress (σ VM ) field while experimentally turning AISI 1045 annealed steel, comparing it afterwards with the experimental results and Usui's model [46,76,77].Experimental validation demonstrates that the recently devised model has significantly enhanced the initial TW rate in both its pattern and magnitude.
Zhang et al. [93] also used Usui's model [46,76,77] to study adhesion and diffuse wear mechanisms, which are found to be the main wear mechanisms during the machining of INCONEL ® 718 with TiAlN-coated carbide tools.The predictive inaccuracy is below 15% compared to the experimental data, proving an enhancement in the precision of TW crater predictions.Binder et al. [94] proposed a better understanding of Usui's model [46,76,77] by considering normal stress (σ n 0 ) in the equation instead of the maximum shear flow stress (τ max ) that can be exerted on the tool.The authors experimentally assessed the orthogonal cutting of AISI 1045 steel with Physical-Vapour-Deposited (PVD) [95] TiAlN-coated and uncoated carbide tools.The predicted crater depth aligns well with experimental data, and the projected position of the crater on the rake face appears to be near the cutting edge; nevertheless, the outcomes of this study indicate that the challenges in predicting the crater width stem from an underestimated contact length in the chip formation simulation.Liu et al. [96] validated numerical simulations performed for six different VB values using Usui's model [46,76,77], culminating in a study on an expression for VB dependent on the spindle (s), feed per tooth (f z ), axial depth of cut (a p or ADOC), and radial depth of cut (a e or EDOC) while end milling TC4 alloy.The authors stated that FEA can directly determine TW and morphology.However, due to the limitations of t and complex boundary conditions, FEA is not commonly used in TL prediction.The TW model's efficacy was confirmed by comparing the experimental findings with the simulation results.Additionally, to expedite tool life prediction, an empirical formula for forecasting TW was formulated, allowing for the characterisation of the evolution of TW over time.Tool compatibility with the workpiece and other parameters unaccounted for in the model have spurred recent research efforts to further enhance machining and TW behaviour predictability.The necessity of having numerous models, as depicted in Table 5, demonstrates how thermodynamically, tribologically, and material-dependent milling, turning, and drilling are, not to mention the complexity of each inherent subject.

Assembly Numerical Modelling
The numerical formulation is only complete if an assembly between the workpiece and a cutting tool exists [97] in order to assess, for instance, the cutting forces (F cut ) [98] inherent to the machining process.This model must encompass the constitutive material, TW models, adequate friction models, mathematically simulated elastoplastic deformations [99], and boundary conditions.CAD software conducts numerical analyses, simulating the tool and workpiece with specified boundary conditions and meshing generation [36].Table 6 offers insights into the available meshing approaches when conducting FEA machining and their advantages and disadvantages.
In traditional practice, the LAG formulation typically modelled metal-cutting processes [100].The computational mesh grid deforms with the material [107]; nevertheless, the concentrated plastic deformation close to the cutting edge significantly distorts the finite elements (FEs), resulting in numerical and computational inaccuracies [108].The Eulerian approach [109] is commonly employed to analyse processes involving substantial deformation [100] in plastic and visco-plastic materials, with minimal consideration for elastic deformation [108] and the mesh grid in which it is fixed in space [107].While the disparity in errors between the Eulerian and LAG approaches varies based on the specific problem, errors within the Eulerian framework generally tend to be more significant than those within the LAG framework [110].The ALE method [111,112] finds extensive use among researchers, as it aids in overcoming certain limitations inherent in the LAG method [108].Integrating the Eulerian and LAG approaches offers a viable solution to address the distortion of elements and effectively model the extensive plastic deformation characteristic of machining [100].Employing an adaptive re-meshing technique aids the process.Nonetheless, it is essential to note that the Eulerian formulation models the local deformation zone, resulting in typically continuous chip formation.Consequently, accurately predicting processes that yield segmented chips proves challenging with this approach.The CEL formulation [54,101] describes the workpiece using the Eulerian formulation and the tool using the LAG approach [113].The CEL meshing approach was not applied to metal cutting until 2016 in a study by Ducobu et al. [113], typically being only used to study fluid-structure interactions until that time.It constitutes a specific variant of the ALE framework.Figure 6 depicts Xu et al. [114]'s CEL [54,101] application to the properties of each domain, namely, the tool and the workpiece, within the numerical model in ABAQUS ™ .The workpiece, including the void and material regions, is meshed using EC3D8RT elements (eight-node thermally coupled linear brick, multimaterial, reduced integration with shear locking effect control [115], also called hourglass effect control), and the tool is meshed as C3D8RT elements (eight-node thermally coupled brick, trilinear displacement and T, reduced integration, shear locking effect control [115]).The boundary conditions and mesh size are also described.While not explicitly stated, the tool is considered a LAG domain.Figure 7 demonstrates the chip morphological difference between the CEL and LAG meshing approaches.Due to the different mathematical formulations in each approach, the two approaches were expected to reveal differences when focusing on the ε p and T distributions.When analysing machining via FEA, it is vital to understand the T distribution, which refers to the varying heat levels during machining, impacting material properties and TW, and ε p , which measures u p , aiding in identifying potential failure areas.Accurate FEA simulations consider cutting parameters, material models, and heat generation to optimise machining processes.During heat generation, a friction model [102] should accurately simulate the interactions between the tool and workpiece during machining.These thermodynamic considerations are critical assessments for carrying out the machining simulation.Table 7 shows some of the most used friction models in machining simulations.
Table 7.Some of the existing friction models for tool-workpiece interaction simulation.

Model Equation
Amonton-Coloumb [115,116] Usui and Shirakashi [115]  For µ calculus, Equation ( 17) presents the dependency of this coefficient on the tangential and radial cutting force components (F t and F r , respectively) and on γ [119].
While the friction models allow for the simulation of the heat and thus T distribution, it is also interesting to study how the generated heat is shared between the tool and workpiece.The heat partition coefficient (HPC) [120] represents the proportion of heat generated during machining that is transferred to the workpiece and the cutting tool.As demonstrated in Equation ( 18), the HPC is a ratio between heat entering the workpiece (Q workpiece ) and the total heat generated during machining (Q total ).
For instance, if HPC = 0.8, 80% of the total heat generated during machining enters the workpiece, while the remaining 20% goes into the tool.Higher HPC values might lead to thermal distortion in the workpiece, whereas lower HPC values lead to increased TW and potentially shorter TL.Once the HPC is calculated for a specific machining process, it can be applied to optimise and improve the machining operation, such as (1) tool material and coating optimisation, (2) machining parameter optimisation, (3) TL prediction, (4) material selection and machining strategy, (5) process improvement and innovation, and (6) heat management and cooling strategies.Table 8 demonstrates no unified model for the HPC since it depends on several material parameters and the affinity between the tool and the workpiece.To better comprehend the proposed models, it is suggested to consult the work of Zhao et al. [121] and Hao and Liu [122].[121,122].

A Predictive Model for HPC Equation Establishment Basis
Loewen Dry-cutting process of AISI 1113 steel with K 2 S cemented carbide tool, v c = 30 − 182 m/min.Shaw [124] HPC SH = Surface grinding process of stainless steel/carbon steel with Al-oxide wheel.

Gecim-Winer [127]
Based on the thermal behaviour of the two-dimensional transient T distribution near a small, stationary, circular heat source equation of the average T of the moving and stationary heat sources between frictional contacts.
Based on the Green function to analyse the chip deformation and friction work along the tool rake face.Berliner-Krajnov [129] HPC Considers Peclet numbers for the tool and workpiece materials in sliding tribological contact.
A SF -area shape factor; l c -tool-chip contact length; q F -frictional heat flux generated in the secondary deformation zone (SDZ); R chip -heat partition coefficient representing heat entering the moving chip from the SDZ; v ch -chip moving speed; α T -tool thermal diffusivity; α W -workpiece thermal diffusivity; ∆θ p max -maximum tool-chip interface temperature rise due to heat generation in PDZ; θ 0 -environmental temperature; λ T -tool thermal conductivity; λ W -workpiece thermal conductivity.
In the extensive experimental work of Zhao et al. [121], it was found that in INCONEL ® 718 machining, the HPC could be calculated by some of the models from Table 8 with a reasonably low sum of relative errors (Σ|δ|) for the different tested v c with the predictive models compared to the measured T. Considering three different v c = 40, 80, and 120 m/min for a PVD AlTiN-coated carbide tool, the lowest Σ|δ| = 52.61%value was from the R R model, and for tungsten carbide (WC), the lowest Σ|δ| = 31.83%was obtained from the R G-W model.The higher Σ|δ| value for the PVD AlTiN-coated tool is explained by the change in the thermal characteristics at the tool-chip contact interface compared with the WC tool.The numerical calculation of conventional processes is complete once the simulation has been constituted with all these considerations.

Hybrid Manufacturing Process Modelling: Thermally Assisted Machining (TAM)
Hybrid non-conventional processes with CM as a base factor have been investigated numerically to efficiently improve Laser-Assisted Machining (LAM) and Induction-Assisted Machining (IAM).Before experiments, thermal FEA of LAM is usually conducted to determine the preheating T and effective a p [130].Taking INCONEL ® 718, for instance, the preferable range of preheating T is between 700 and 900 • C, where σ u sharply decreases [131].On the other hand, the INCONEL ® 625 preheating T was chosen by Parida and Maity [132] to be 600 • C. Thermal FEA for laser heating simulation is governed by Equation (19) [133]: where α and .
Q are thermal diffusivity (Equation ( 20)) and power generation per unit volume, and the constants ρ, c p , and k represent the volumetric mass density, specific heat, and thermal conductivity, respectively.Equation (21) gives the initial condition at t = 0 s: Equation ( 22) [133] gives the boundary condition: where q, h, ε e , σ SB , T s , and T sur represent the heat flux, heat transfer coefficient, emissivity, Stefan-Boltzmann constant, surface temperature, and radiated surface temperature, respectively.The Gaussian Equation (Equation ( 23) [134]) is employed to elucidate the distribution of laser power on the focal plane of the laser beam: where P Laser is the laser's total power input, and r is the laser spot radius.This numerical approach to LAM is applied to hybridised milling and turning operations.Figure 8 illustrates a 3D CAD model to plan the laser heat source path.Electromagnetic FEA is a suitable type of numerical analysis for IAM.The primary aspects for the investigation of IAM are the dwell t required to attain the intended preheating T and the EDOC.The governing equation of the electromagnetic analysis can be expressed using Maxwell's Equation ( 24) [131,135,136]: where µ 0 is the vacuum magnetic permeability, µ r is the relative magnetic permeability, A mag is the magnetic vector potential, σ elc is the electrical conductivity, and J s is the source current density.As described in the LAM section, the preferable range of preheating T for INCONEL ® 718 and INCONEL ® 625 is the same.Following the presentation of the theoretical framework in Section 1, Section 2 delineates the methodology employed in this study, which is based on the Systematic Literature Review (SLR) approach [137] aimed at summarising how the research was conducted.In Section 3, the state of the art is demonstrated for different types of machining, depicting evolutionary trends and remarks from the researchers' work.Section 4 discusses findings derived from content analysis, and Section 5 succinctly summarises the findings and offers a brief outlook on INCONEL ® alloy machining simulation.

Materials and Methods
The research and information-compiling phases were carried out through the SLR approach since it is based on a systematic, method-driven, and replicable approach [138,139].The platform used for SLR was Dimensions.ai, which is connected to all data in Scopus through quality criteria by consulting each journal's influence within the academic community on the basis of the impact factor (IF). Thanks to Elsevier, Springer, AIP Publishing, ASME, MDPI, and IMeche articles and books from Woodhead Publishing, Butterworth-Heinemann, Elsevier, Academic Press, CRC Press, and Publindustria, this article collects information from 125 reports, 25 book chapters, and 5 standards from 266 articles and 20 books researched.The research procedure is described next: 1.
Information was searched with the "INCONEL ® 718", "INCONEL ® 625", and "Johnson-Cook criteria" keywords to gather a broad range of information about material modelling.

2.
The keywords "FEA" and "Traditional machining processes" were added to enable a search for information about the numerical modelling of traditional manufacturing processes.To refine the data even further, "Numerical models", "ABAQUS ™ ", "ANSYS ® ", and "DEFORM ® " keywords were included to gather the desired information.

3.
After collecting the articles, the journal's influence was evaluated with its Web of Science score from 2022 (ignoring quartiles).All journals with an IF value less than three were excluded, although rounding to the unit was allowed.

4.
The abstracts and conclusions from the collected articles were analysed.

5.
Knowledge from 2013 to 2023 about the modelling of traditional INCONEL ® 718 and 625 processes was compiled.
The review article regarding INCONEL ® alloy's numerical analysis on traditional machining predictability was written down.

CM Processes
This section addresses FEA applied to more traditional manufacturing processes, such as milling, turning, and drilling.

Milling
Jia et al. [140] developed a simulation model using ABAQUS ™ to replicate a milling experiment involving the micro-milling of thin-walled parts made of INCONEL ® 718.The aim was to predict the cutting force (F c ) and anticipate wall deflection in these parts using an element birth/death technique.The mesh type of the tool model was set as C3D4 (four-node tetrahedral element for general purpose), and the mesh type of the thin-walled part model was set as C3D8R (eight-node brick element with reduced integration and shear locking effect control [141]).Comparing the simulation output for F c with experimentally measured values, the model showcased maximum and average relative errors of 7.8% and 2.2%, respectively, validating its accuracy.The results underscored the accuracy of the proposed simulation model for micro-milling thin-walled parts, as depicted in Figure 9, and highlighted the precision of the F c prediction.
Okafor and Sultan [142] constructed a mechanistic model aimed at predicting F c during the high-speed end milling of INCONEL ® 718, employing a wavy-edge, bull-nose, helical end mill (WEBNHE) as illustrated in Figure 10.The authors proceeded to conduct end-milling experiments to validate the output of the mechanistic model.The MATLAB ® mathematical model of the end-mill geometry was simulated.Upon comparing the predicted F c values with the experimentally measured ones, the authors concluded that the values in Table 9 accurately supported the mechanistic model.Ducroux et al. [143] (Figure 11) presented numerical and experimental milling comparisons between wrought and additively manufactured (AM) INCONEL ® 718 specimens.The machinability was investigated through microstructure observation and F c analysis, followed by TW characterisations for both specimens.Novel formulations of the F c model in milling were developed and modelled with a fully parameterised mechanistic approach.Additionally, the tool geometry and the model concerning local forces consider the VB effect.The additively manufactured INCONEL ® 718 showed better machinability.TL was double compared with the wrought-stock.Considering that TW significantly and directly affects F c , a local model formulation, as proposed with the tool geometry evolution, has to be considered to improve the model precision; however, it is challenging to consider VB max > 0.25 mm in a predictive model, as it creates too many random TW profiles.

Turning
Díaz-Álvarez et al. [45] evaluated the progression of TW and the T distribution during the turning of INCONEL ® 718 in a TiN-coated insert.This assessment encompassed both experimental and FEA approaches in DEFORM ® 3D software, as depicted in Figure 13, varying the tool cutting edge angle (κ r ).The numerical model was validated through turning experiments, and the principal wear mechanisms, such as chipping, notching, and built-up edge (BUE), were compared with variables forecasted by the numerical model.Zhang et al. [93] devised an innovative numerical forecasting approach for the tool VB rate of PVD TiAlN-coated carbide tools.This method is grounded on the positive feedback interconnection between the tool geometry and TW rate during the dry orthogonal cutting of INCONEL ® 718, featuring various scenarios, as demonstrated in Figure 14.The numerical model made in AdvantEdge ® v7.0 from Third Wave System © has a margin of error below 15% when juxtaposed with experimental findings, highlighting the potential of this method to predict TW accurately.Yadav et al. [46] experimentally and numerically assessed the material removal rate (MRR) and VB during the turning of INCONEL ® 718.A Chemical Vapour Deposition (CVD)-coated WC tool was used for experimental studies.In parallel, a 3D model was developed using the LAG approach in the DEFORM ® 3D software.The comparative error between the experimental and simulated results for the MRR and VB remained under 7%.The simulations consistently projected higher values for MRR and VB than those observed in the experimental data (Figure 15).The authors concluded that the simulations could significantly economise valuable t and material resources owing to the demonstrated accuracy.Xu et al. [47] devised an experimental and numerical investigation of the turning process for INCONEL ® 718 utilising worn tools.The study encompassed an assessment of chip formation, as depicted in Figure 16, force variation, TW, and thermal distribution in ABAQUS ™ .The element type used for the workpiece was EC3D8RT, and the tool was C3D8T (an eight-node trilinear thermally coupled brick element).The numerical meshing approach was CEL.Empirical findings revealed that, at higher v c , more distinct serrated chips were evident, with a narrower shear band compared to lower v c .The study demonstrated a strong agreement between the experimentally observed results and the simulated outputs, particularly concerning F c , chip morphology, and thermal distribution, which substantiates the capability of the numerical model to replicate the machining process, even with worn tools, accurately.Jafarian et al. [145] proposed a new and efficient method based on an evolutionary optimisation algorithm to identify JC material constants for INCONEL ® 718.The subsequent orthogonal cutting process of INCONEL ® 718 was numerically simulated in DEFORM ® 2D using the new material model and two tool edges (chamfered and honed), and the outcomes were compared with experimental results.The numerical F c , maximum temperature (T max ), and chip geometry, the last one shown in Figure 17, were validated.The overall error of the simulations in the first and second steps was reported at 12.8% and 11.3%, respectively, enhancing the knowledge of modelling surface integrity in machining processes.López-Gálvez and Soldani [48] focused on determining optimal numerical parameters, such as mesh size and element deletion criterion, to model the chip-removing process, as Figure 18 shows.The numerical model was built in ABAQUS ™ with a LAG approach.The thermo-mechanical analysis was realised using C3D8R elements.The comparison between numerical and experimental measurements exhibited excellent accuracy, manifested by an average relative error of 2% between numerical and experimental F c , considering v c = 300 m/min.Gong et al. [146] applied DOE to a developed numerical model in AdvantEdge ® to assess a turning simulation of INCONEL ® 718 in dry and cryogenic N 2 (l) cooling/lubrication environments and then compared the results with T cut and F c obtained in experimental tests carried out with the same input parameters.As illustrated in Figure 19, numerically, it was found that the cryogenic strategy reduced T cut by 11.68%, 9.42%, and 11.98% for v c = 70 m/min, 100 m/min, and 130 m/min, respectively.On the other hand, the peak tool temperature (T t ) was reduced in cryogenic conditions by 9.36%, 8.81%, and 9.23% for v c = 70 m/min, 100 m/min, and 130 m/min, respectively.These values enhance the predictability of numerical modelling for improving machining performance.
Jafarian et al. [147] adopted new strategies to improve the accuracy of FEA modelling of the INCONEL ® 718 cutting process in the software DEFORM ® 2D v10.0.A novel hybrid strategy was established to simultaneously calibrate controllable simulation parameters, implemented based on DOE, intelligent systems, and the FEA of the cutting process.The numerical data were validated based on experimental results, verifying significant improvements in predicting the chip geometry, F c , and T max (the last two portrayed in Figure 20).Experimental grain size and hardness variation results were compared with the corresponding numerical counterparts obtained using the JC and new material models.It was shown that implementing a grain-size-based method increased the simulation accuracy.Qiu et al. [148] built a 3D FEA model in AdvantEdge ® v7.1, and the simulation's reliability was confirmed through a comparative analysis of the outcomes with a hybrid orthogonal experiment, as shown in Figure 21.A sensitivity analysis of JC parameters and µ in the simulation results under low-and high-v c conditions was carried out.The authors found that the simulation accuracy of INCONEL ® 718 is susceptible to ε hardening and thermal softening parameter evaluation with the JC constitutive model, which significantly influences the residual stress, h ch , F c , and T. The µ parameter only significantly affects F a and F r in the high-v c condition.
Liu et al. [54] experimentally and numerically analysed the effect of the tool geometry on the thermal-mechanical load and residual stresses in the orthogonal machining of INCONEL ® 718.In ABAQUS ™ , the CEL method was selected to simulate the effect of the tool geometry on T, F c , PEEQ, and residual stresses.The element type used in the Eulerian-part domain was EC3D8RT, and the tool was defined as rigid by using an eightnode thermally coupled brick element.The experimental results validated the numerical model, and the CEL model was able to cope with reality regarding the T distribution, F c, and residual stress profiles.Figure 22     Tian et al. [18] conducted a reverse identification of JC parameters for INCONEL ® 718, considering specimens in both the solution-annealed and precipitation-hardened states.Experimental orthogonal cutting was employed to obtain the average F c and h ch , validating the FEA model produced in DEFORM ® 2D v12.0.In the experimental analysis, the Waldorf model [149] was integrated to describe the deformation characteristics (ε, T, .ε, and σ) of the primary shear zone, incorporating the influence of the cutting-edge radius on F c during the cutting process.Figure 23 demonstrates an evaluation of how TW affects the material separation point during turning.The study determined JC constitutive constants for specimens in both solution-annealed and precipitation-hardened states, revealing a notable divergence in the mechanical response as evidenced in the constitutive parameters, F c , and chip morphology.The error of F c was found to be within 5%, indicating that the proposed method significantly improves the accuracy of the simulation results.Furthermore, it was established that the type of TW influenced the distribution pattern of the thermalmechanical load during cutting by altering the length of the tool-workpiece contact and the position of the separation point.However, the authors stated that the proposed method only applies to low v c in challenging-to-machine materials.Peng et al. [150] focused on proposing a novel empirical equation for predicting residual stress on the turned surface of INCONEL ® 718, considering tool parameters such as κ r , γ, and inclination angle (λ s ).Initially, the accuracy of the parameters set for the 3D FEA model through AdvantEdge ® v7.4015 was validated by comparing the residual stresses and chips with the experimental results (Figure 24).According to the obtained results, the predictive equation of surface residual stress accurately predicted turned surface residual stress for INCONEL ® 718 materials.However, some flaws arose: γ and λ s cannot be zero at the same time, and |γ| cannot equal |λ s | when κ r = 45 • or 135 • .The absolute mean discrepancy and the Pearson correlation coefficient (R) were 13.40% and 0.9624, respectively, indicating that this model can be used in the real-time monitoring of turning INCONEL ® 718 and can also predict the residual stress level to guide process planning.Tu et al. [151] systematically investigated the TW characteristic of cubic Boron Nitride (cBN) at different v c values in the dry turning of INCONEL ® 718 in experiments and FEA simulations (Figure 25).The FEA model developed in AdvantEdge ® with initial VB was established and validated.The parameter µ as a function of VB between the tool and workpiece in the turning process was determined to be in the range of 0.3-0.7 by F c component measurements.Differences were attributed to the varied TW and F c .FEA simulations visually presented the tool flank T distribution dependent on VB.T max was located at the tool cutting edge and tooltip and decreased in the tool substrate.The dominant TW mechanisms were described as adhesive wear and cutting-edge chipping at the flank face and rake face.The tool failed due to the catastrophic fracture of the cutting edge.Kortabarria et al. [152] estimated machining-induced residual stresses in INCONEL ® 718 orthogonal cutting.An FEA model was constructed using DEFORM ® 2D v10.2 software, and a sensitivity analysis was performed to evaluate how the input data within the model affected the predicted residual stresses, as depicted in Figure 26.After analyses, it was concluded that only the material constitutive law was significant in predicting residual stresses.Moreover, the material behaviour at a high heating rate was crucial for realising accurate predictions.
Pervaiz et al. [153] conducted a machinability analysis for INCONEL ® 718 using a numerical approach undertaken in AdvantEdge ® , substantiated by experimental validations.Dry and N 2 (l) methods were compared based on relevant parameters, including the chip compression ratio, shear angle, contact length, F c , and energy consumption (P in ) for the primary deformation zone (PDZ) and SDZ.In addition, parameters related to chip morphology were also investigated.Observations from experimentation and numerical analysis (Figure 27) demonstrated a lower chip compression ratio when using N 2 (l) cooling, accompanied by a larger shear plane angle.Consequently, this configuration reduced the tool-to-chip contact length and improved lubrication.4) upon contact, the stress direction at point P1 shifted towards tensile, followed by a rapid cooling process and the relaxation of compressive stresses; (5) subsequently, due to a gradual cooling process, the stress at point P1 transitioned towards tensile [152].Razanica et al. [53] suggested an innovative model, the Progressive Ductile Damage (PDD) model, to depict and simulate the chip formation process of INCONEL ® 718, incorporating a novel approach based on rigid visco-plastic flow and continuous damage evolution.FEA was performed using DEFORM ® 2D software.The resultant material response with this approach is depicted in Figure 28.The JC model was employed at the onset of damage initiation, while a modified CL failure criterion controlled the damage evolution.The accuracy of numerical outcomes was validated through experimental machining tests by observing, analysing, and comparing the formation of continuous and serrated chips (depending on v c ), F c , chip shapes, and tool-chip contact lengths.Bücker et al. [49] presented a new approach that combines high-pressure fluid supply and high cooling performance using a Deep Temperature Emulsion (DTE) to improve cooling effectiveness.DEFORM ® 3D FEA was carried out on the high-speed turning of INCONEL ® 718, and the model (Figure 29) was validated through experimental measurements of F c , TW, chip formation, and surface integrity.Machining and heat transfer FEA at the cutting zone indicated the importance of cutting fluid penetrability into the cutting zone due to the low thermal conductivity of the workpiece, chips, and cutting tool, which augments conductive heat transfer.

Drilling
In the context of traditional drilling processes, Chenegrin et al. [103] carried out computational calculations in a numerical model focused on INCONEL ® 718 drilling.The aim was to analyse heat transfer as a function of drilling parameters and analytically depict chip geometries (as shown in Figure 30) using process parameters and the drill tool's geometry.In the 2D Merchant's software, the R-ALE method modelling approach was applied to simulate the drilling of a 17 mm deep hole in a cylindrical workpiece.Notably, the numerical model successfully predicted the T distribution at the hole's mid-depth and the drilling operation's beginning and end.Attanasio et al. [154] aimed to model (Figure 31a,b) and simulate TW (Figure 31c,d) in the drilling of INCONEL ® 718 using an innovative numerical procedure to update the worn tool's geometry in DEFORM ® 3D software.Experimental tests were performed to measure TW in drilling and to validate the numerical models.Conventional Metal Working Fluids (MWFs) [155] and N 2 (l) cryogenic cooling environments were used in both approaches.The comparison between simulated and measured results demonstrated the suitability of the developed drilling model to predict TW under MWF and N 2 (l) cooling conditions.Thus, the formulated model is adept at effectively assessing the impact of the cutting process and cooling conditions on TW, minimising the number of expensive tests.

Non-Conventional Manufacturing Processes
This section addresses FEA applied to non-conventional manufacturing processes.

Electrochemical Machining (ECM)
In the context of non-conventional manufacturing, such as ECM, Niu et al. [156] concluded that enhanced tool designs for the initial machining phase could be achieved by modifying the configuration and arrangement of tool-sidewall outlet holes in the electrochemical mill grinding (ECMG) of INCONEL ® 718 using the FLUENT 15 computational fluid dynamics software.Four distinct tools, each featuring varying numbers of rows of tool-sidewall outlet holes, were conceived and subjected to numerical analysis.The FEA outcomes demonstrated that an abrasive tool with four rows of tool-sidewall outlet holes yielded a superior MRR, f parameter, and electrolyte pressure due to an optimal flow velocity distribution, as depicted in Figure 32.Experimental assessments of machining a slot using this tool indicated an amplified MRR with a higher applied voltage (U), electrolyte pressure, and f.In the rough machining phase, the original tool yielded an average sidewall flatness and an arithmetic average of profile height deviation (R a ) of 549.6 µm and 2.509 µm, respectively.The results of the average sidewall flatness and sidewall R a attained with the newer tool were markedly reduced to 340.5 µm and 1.65 µm, respectively.The new tool was also employed for finishing machining, resulting in an average sidewall flatness and sidewall R a of 69.5 µm and 0.648 µm, respectively.

Thermally Assisted Machining (TAM)
Hybrid manufacturing processes like TAM can feature conventional milling and lathing.Kim and Lee [130] investigated the deployment of LAM for the machining of an INCONEL ® 718 workpiece with a 3D curved shape by employing Non-Uniform Rational B-Spline (NURBS) techniques and employing different tool paths, such as ramping and contouring, while varying machining conditions.The authors utilised ANSYS ® v.18 software for the numerical evaluation of the varied tool paths, employing different meshing approaches, as depicted in Figure 33.Furthermore, the authors examined the T resulting from laser heating on the INCONEL ® 718 workpiece.Experimentally, they scrutinised a p obtained in thermal FEA, F c , Specific Cutting Energy (SCE), and R a .These parameters were meticulously measured, analysed, and set against CM and LAM.Remarkably, LAM demonstrated a significant enhancement in these machining characteristics compared to conventional machining.The numerical results obtained in the reported study employed experimental LAM endeavours involving a diversity of 3D shapes.
Kim and Lee [131] assessed machining efficiency in the context of INCONEL ® 718 machining, utilising P in as a metric.The study encompassed three machining approaches: CM, LAM, and IAM, as depicted in Figure 34a,b.Utilising ANSYS ® software, the researchers conducted numerical thermal and thermal-electromagnetic analyses to validate the EDOC and dwell time resulting from preheating, presented in Figure 34c-f, revealing that the optimal dwell t = 20 s, and the effective a p = 0.5 mm.Parameters including P in , F c , and R a were thoroughly experimentally analysed under LAM and IAM approaches' machining conditions.P in increased with an increase in f, primarily due to the necessity to sustain the preheating T. Notably, IAM exhibited the highest P in as it necessitated a lengthier preheating duration due to lower thermal energy concentration.Compared with CM, a 32% increase in P in was observed in LAM, accompanied by 41% and 51% reductions in F c and R a , respectively.Similarly, when compared with CM, a 66% increase in P in was noted in LAM, alongside 45% and 32% decreases in F c and R a , respectively.It was observed that both LAM and IAM consumed more power than CM; nevertheless, they exhibited enhanced machinability.Ultimately, this study identified LAM as the most suitable hybrid machining process.
Jeong and Lee [133] made significant strides in enhancing the machinability of INCONEL ® 718 by applying LAM and a heat shield, as illustrated in Figure 35a,b.The authors conducted experimental analyses, specifically TL experiments, with a uniform MRR while varying the usage of the heat shield.A numerical thermal analysis was also conducted in ANSYS ® software to ascertain the EDOC.The study provided a comprehensive account of the alterations in TL and machining efficiency resulting from implementing the heat shield, elucidating its effectiveness.The FEA thermal analysis determined that an optimal a p = 0.3 mm could be achieved with a preheating T = 900 • C. In terms of T t , it was observed that, without the heat shield, T t = 643.73• C (Figure 35c), whereas with the heat shield, T t = 535.56• C (Figure 35d), effectively mitigating the dissipation of laser-generated heat to the tool.Concerning TL, LAM increased it by approximately 53% compared to CM.Conversely, using LAM in conjunction with the heat shield resulted in a remarkable 78.3% increase in TL compared to CM, attributed to the protective function of the heat shield, which reduced the detrimental effects of heat transferred to the tool, thereby minimising abrasions, cracks, and fractures in the tool.Zhang et al. [51] employed LAM to enhance the machinability of INCONEL ® 718 and focused on investigating the Heat-Affected Zone (HAZ).The paper details a proposed and developed theoretical model to describe the gradient of σ within the HAZ during the LAM process on INCONEL ® 718 based on the workpiece's strengthening mechanisms, microstructure, and T distribution.Comparative analysis with the JC equation revealed that the model successfully predicts the HAZ range, quantitatively characterises the EDOC within the HAZ, and also predicts σ of the HAZ and determines its range during LAM, as depicted in Figures 36a and 36b, respectively.The T amplitude at different depths decreased with increasing distance from the surface.The error between the σ predicted by the proposed model at T 0 and the experimental data remained within a 5% tolerance, confirming the accuracy of the predicted results, as they aligned with the experimental outcomes.Additionally, it was observed that the gradient of σ decreased with higher P Laser and reduced laser scanning speed.During the heating process, when T < 700 • C, the precipitation strengthening component (σ p ) contributed to over 50% of σ, surpassing the contributions of the solid solution strengthening component (σ SS ), grain boundary strengthening component (σ D ), and Ni-matrix strength component (σ Ni ).As the T amplitude increased, σ p significantly decreased, σ SS increased, and σ D and σ Ni experienced a gradual reduction.The slow rise in σ during the cooling stage was attributed to σ D and σ Ni .
Zhang et al. [157] assessed the Laser-Assisted Micro-Milling (LAMM) process of INCONEL ® 718 using a combination of experiments and FEA simulations in ABAQUS ™ (Figure 37).First, a 3D thermal-mechanical coupled FEA model of LAM was established, followed by experiments and FEA simulations to investigate the impact of LAM on F c , chip morphology, TW, and surface topography.The milling tool was discretised by R3D4 (rigid tetrahedral elements).C3D8RT elements were considered to discretise the workpiece, and the Coloumb [115,116] friction model described the friction behaviour between the milling tool and the workpiece with a µ = 0.2.The results indicated that LAM enables σ u and G reduction, which improved the cutting performance, reducing F c by 40.5% while enhancing surface integrity and chip continuity and prolonging TL, compared with CM.

Laser-Beam Machining and Laser Drilling Machining (LBM and LDM)
In the context of LBM, Pan et al. [134] conducted an experimental investigation to characterise the shape of the melting zone for INCONEL ® 718.The authors explored the impact of the laser scanning speed, utilising a trochoidal path and varying P Laser , on the absorption ratio.Additionally, a 3D FEA was proposed to predict the T distribution in parallel with experiments by utilising ANSYS ® , as portrayed in Figure 38.This FEA successfully predicted the T distribution field of the coaxial laser preheating system for INCONEL ® 718 LBM, closely approximating the experimental conditions.The study observed that the melting zone area, depth, and width (MZA, MZD, and MZW, respectively) consistently decreased with increasing laser scanning speed for each P Laser value utilised (400, 600, and 1000 W).Furthermore, higher P Laser values were found to positively influence MZA's increase.

Discussion
Considering all the reviewed content, some previous conclusions can be drawn and discussed.Starting with the work of Chenegrin et al. [103], the authors demonstrated that chip formation is crucial for drilling and milling; nonetheless, pre-modelling may be necessary to start the material cutting.Additionally, it was noted that most of the reviewed works were based on orthogonal cutting.This situation is tied to the fact that the measuring equipment used in turning is more accessible for the practical validation of numerical values than the equipment used in milling or drilling, posing a challenge for researchers attempting milling operations in association with FEA.Regarding numerical strategies, Outeiro [6] and Mir et al. [100] highlighted a significant limitation of both ALE and CEL formulations for metal-cutting simulations-the inability to accurately simulate chip geometry, which is particularly problematic for the serrated chip formation observed when machining challenging-to-cut materials like INCONEL ® alloys and is a subject that needs future improvements in order to fine-tune the accuracy of the outcomes.Iturbe et al. [10] emphasised the importance of considering thermal softening phenomena in initiating and amplifying chipping morphology, leading to the formation of adiabatic shear bands in machining simulations, and it was found that Hokka et al. [42] improved the numerical results by applying an additional ε softening term to the INCONEL ® 625-like alloy JC model, allowing a decrease in .ε hardening at large deformations.On the other hand, Wang and Liu [65] demonstrated that JC parameters have a notable impact on the chip shape for Ti6Al4V, underscoring the need to adjust JC constants through experimental results, which, in this case, can be very much extrapolated to the INCONEL ® alloy study.

Conclusions
A concise overview of FEA strategies is presented to systematically summarise and analyse the recent advancements in predicting INCONEL ® machining using FEA from 2013 to 2023.This review also highlights the most recent numerical solutions, prospects, and limitations that researchers have faced.The numerical predictive models developed by these researchers were validated through experimental testing.Notably, a prevalent trend among researchers was observed:

•
Seventeen reviews for orthogonal cutting or turning; • Seven reviews for milling; • Two reviews for drilling; • One review for LBM/LDM.
Hybrid manufacturing processes were considered within each category, depending on CM.Regarding the preference for FEA turning in the context of INCONEL ® machining, it can be explained as follows: • FEA models of the turning process are more straightforward to simulate than those for milling or drilling processes.
Regarding the model design, two primary approaches were identified: • The tool is fixed and rigid while the material is moved against it; • The material is secured in the chuck (as in real-life scenarios) while the tool moves along the workpiece.
Nowadays, there is still a necessity to have numerous models for TW behaviour predictability since machining is a very highly thermodynamically, tribologically, and material-dependent process.In this regard, the following is noticeable:

•
Most of the authors prefer Usui's diffusive wear model since Taylor's empirical TL model requires extensive calibration work, and material adhesion is difficult to measure (considering all adhesion models addressed); • The JC constitutive model does not adequately address thermal softening phenomena, which may lead to divergent outcomes when compared to experimental results; • Although it has been determined that the major setback in terms of numerical approaches is the material constitutive model chosen, FEA applied to the machining of INCONEL ® alloys has shown promise in providing accurate results, effectively addressing certain drawbacks inherent in JC application in numerical material modelling, a preferred approach for many researchers; • There is high confidence that future studies will increasingly utilise FEA to predict and support experimental tests by improving INCONEL ® alloy machining and material modelling and considering all thermodynamical and tribological phenomena in the tool-workpiece interface.
INCONEL ® machining simulation through FEA effectively enhances experimental analyses.It provides a more precise calibration of parameters, the improved prediction of process T, a more accurate forecast of TW and surface quality assessment, and a more informed selection of potential coatings for tools.This knowledge is poised to facilitate broader accessibility to INCONEL ® machining for researchers and practitioners alike.Funding: This work was developed under the "DRIVOLUTION-Transition to the factory of the future", with the reference DRIVOLUTION C644913740-00000022 research project, supported by European Structural and Investments Funds with the "Portugal2020" programme scope.

Figure 1 .
Figure 1.Subdivision of the JC constitutive model equation according to distinct deformation phenomena and the corresponding experimental trials and conditions [38].
p f -equivalent plastic displacement at failure; Y-geometrical constant; w e -specific essential work of fracture; W e -essential work of fracture; W f -work of fracture; W n1 -strain energy density; w p -specific work of propagation; W p -work of propagation; β-shape parameter; Γ-chosen boundary; [σ]-Cauchy stress tensor; σ u -tensile strength.

Figure 4 .
Figure 4.The development of chip morphology is impacted by prior machining operations, including (a) the influence of the initial chip thickness (h ch ), (b) the impact of the rake angle (γ) of the tool, (c) the ramifications of the tool's edge radius (r β ), and (d) the effects of intervening cuts [60].

Figure 6 .
Figure 6.Orthogonal cutting model with a minimal width of cut (merely 4 µm) employing the CEL approach: (a) mesh structure and (b) specifications of boundary conditions [114].

Figure 7 .
Figure 7. Difference in numerical behaviour between the CEL and LAG approaches when studying ε p (a,b) and T (c,d).Cutting speed, v c = 250 m/min [114].

Figure 8 .
Figure 8.The 3D model and conceptualisation of a mobile laser heat source [133].

Figure 10 .
Figure 10.(a) The WEBNHE partially immersed radially during up-milling; (b) the axial segmentation of the end mill and the various F t , F r , and axial cutting force (F a ) components experienced at the cutting edge of a representative disk; (c) the distinct F t , F r , and F a [142].

Figure 11 .
Figure 11.Graphical summary of Ducroux et al. [143]'s work.Zheng et al. [144] numerically (Figure 12) and experimentally assessed the machinability of INCONEL ® 718 in milling.The FEA model, produced in AdvantEdge ® v7.1 from Third Wave System © , was used to analyse the stress field, T distribution, and F c .Milling experiments on INCONEL ® 718 were designed, the results obtained were compared with the FEA model to validate it, and MATLAB ® was used to further optimise the experimental processing parameters.FEA revealed notable T increases in the machining zone with rising s, a e , a p , and f z .Validation against experiments and the calculated results from the milling force model and FEA exhibited minor deviations.The optimal experimental parameters were s = 3199.2rpm, a feed rate (f ) of 80 mm/min, and a p = 0.25 mm.

Figure 14 .
Figure 14.Distribution of σ n 0 on a worn tool: (a) v c = 50 m/min and cutting time (t cut ) 100 s; (b) v c = 70 m/min and t cut = 35 s [93].

Figure 15 .
Figure 15.(a) Comparative analysis of simulated and experimental outcomes for VB; (b) comparative assessment of simulated and experimental results for MRR (adapted from [46]).

Figure 17 .
Figure 17.FEA simulation illustrating the creation of serrated chips during the orthogonal cutting procedure of INCONEL ® 718 alloy [145].

Figure 19 .
Figure 19.Fluctuation of T t in different T distributions concerning v c and different lubrication/cooling conditions, maintaining a constant f of 0.2 mm/rev and a p = 0.2 mm constant: (a) dry v c = 70 m/min, (b) cryogenic v c = 70 m/min, (c) dry v c = 100 m/min, (d) cryogenic v c = 100 m/min, (e) dry v c = 130 m/min, (f) cryogenic v c = 130 m/min [146].
demonstrates the T distribution around the tool edge and along the workpiece depth, influenced by different tool geometries when turning INCONEL ® 718.

Figure 21 .
Figure 21.h ch compared experimentally and numerically in low-and high-v c conditions [148].

Figure 23 .
Figure 23.Effect of different wear types on material separation point with (a) new tool, (b) rounded cutting edge, (c) worn rake face and (d) worn flank face [18].

Figure 24 .
Figure 24.The comparison of h ch between experiments and simulations [150].

Figure 25 .
Figure 25.Stress profile in INCONEL ® 718 emanating from the machined surface using a cBN tool at different VB and (a) v c = 200 m/min and (b) v c = 350 m/min [151].

Figure 26 .
Figure 26.Progression of the near-surface point P1 during machining at T 0 : (1) the material adjacent to point P1 experienced compressive stress; (2) point P1 attained the peak compressive stress; (3) point P1 came into contact with the tool, reaching T max = 600 • C; (4) upon contact, the stress direction at point P1 shifted towards tensile, followed by a rapid cooling process and the relaxation of compressive stresses; (5) subsequently, due to a gradual cooling process, the stress at point P1 transitioned towards tensile[152].

Figure 27 .
Figure 27.Analysis encompassing the comparison of chip morphologies and related assessments between experimentally obtained chips and their simulated counterparts [153].

Figure 28 .
Figure 28.Comparison of T, CL, and damage fields in the PDD modelling approach at varying v c together with the experimental chip shapes [53].

Figure 30 .
Figure 30.Configuration of the chip based on the drill tool geometry and the process parameters [103].

Figure 31 .
Figure 31.(a) A 3D FEA model representing the drilling process, (b) elaborate mesh details concerning the cutting edge of the drill bit.Simulations were conducted to analyse (c) T and (d) σ eq distributions within the drill under cryogenic cooling after t cut = 6 min duration [154].

Figure 32 .
Figure 32.Numerical analysis on flow velocity distributions of the different initial arrangements of ECMG tools with distances between the hole centres of any two adjacent rows of (a) 6.8 mm, (b) 3.4 mm, (c) 2.3 mm, and (d) 1.7 mm [156].

Figure 35 .
Figure 35.Conceptual diagram of (a) LAM and (b) LAM with heat shield, (c) thermal FEA of LAM, (d) thermal FEA of LAM with heat shield (adapted from [133]).

Figure 36 .
Figure 36.Distribution of T and σ along the depth direction during laser scan, (a) T amplitude at different depths decreases with the increase of the distance from the surface, (b) as time changes, σ y at different depths decreases first and then increases slightly with the T change [51].

Figure 38 .
Figure 38.(a) The top-section and (b) cross-section perspectives of the T distribution generated by a laser rotating at 3500 rpm.(c) The top-section and (d) cross-section perspectives at a rotational speed of 7000 rpm, featuring r = 0.2 mm and a moving speed of 1000 mm/min (following a trochoidal path, with T denoted in units of Kelvin).[134].

Table 2 .
Values of A, B, C, m, and n for INCONEL ® 718 and 625 were obtained from different authors.
provides values of the INCONEL ® 718 damage constants for the JC shear fracture damage model.For INCONEL ® 625, a literature gap exists since no values were found.

Table 3 .
The constant parameters of INCONEL ® alloys for JC shear fracture damage model.

Table 8 .
Predictive models based on physics and experiments for HPC

Table 9 .
[142]um percentile prediction error for the cutting forces according to the orthogonal axis (F x , F y , and F z ) components at 62 and 93 rpm[142].