Milling Mechanism and Chattering Stability of Nickel-Based Superalloy Inconel 718

Nickel-based superalloy Inconel 718 is widely used in the aerospace industry for its excellent high-temperature strength and thermal stability. However, milling Inconel 718 presents challenges because of the significantly increased cutting force and vibration, since Inconel 718 is a typical difficult-to-machine material. This paper takes the milling process of Inconel 718 as the research object, initially, and a milling force model of Inconel 718 is established. Subsequently, the finite element analysis method is used to analyze the stress field, temperature field, and milling force in the milling process of Inconel 718. Building upon this, a dynamic equation of the milling of Inconel 718 is established, and based on the modal experiment, stability lobe diagrams are drawn. Furthermore, milling experiments on Inconel 718 are designed, and the results calculated using the milling force model and finite element analysis are verified through comparison to the experimental results; then, the fmincon optimization algorithm is used to optimize the processing parameters of Inconel 718. Eventually, the results of the multi-objective optimization illustrate that the best processing parameters are a spindle speed of 3199.2 rpm and a feed speed of 80 mm/min with an axial depth of cut of 0.25 mm. Based on this, the best machining parameters are determined, which point towards an improvement of the machining efficiency and quality of Inconel 718.


Introduction
Nickel-based superalloy Inconel 718 is widely used in the aerospace industry because of its extraordinary high-temperature strength, thermal stability, heat fatigue resistance, corrosion resistance, creep resistance, and oxidation resistance, especially in the manufacturing of key components for the aerospace industry, such as turbine blades, combustion chamber components, and gas turbines [1][2][3][4][5][6].Common processing methods for Inconel 718 are mainly casting [7,8], forging [9], additive manufacturing [10,11], and electrical discharge machining [12].However, each of these machining methods suffers from flaws, like low dimensional accuracy; inability to process complex geometry; high machining costs; and containing heat-affected areas, slow machining speed, and poor surface quality, respectively.Although utilizing the milling method can effectively avoid the abovementioned problems, in the cutting process Inconel 718 is prone to serious tool wear [13], high thermal stress, high cutting forces [14], vibration, and other phenomena [15].To address these challenges and improve the machinability of Inconel 718, the milling mechanisms and chattering stability of Inconel 718 need to be investigated.Based on this, the optimized processing parameters can be determined by analyzing the influence of different processing parameters, such as spindle speed, feed per tooth, radial depth of cut, and axial depth of cut on the cutting force, temperature, stress, and deformation, including machining quality during manufacturing.
Accurate modeling of the milling force can help adjust the cutting parameters to ensure high-quality machining results; therefore, it is especially necessary to establish an accurate milling force model.At present, the milling force models proposed by scholars mainly include empirical models, physical analytical models, and finite element models.
The empirical model of milling force is a model derived by fitting milling test data and is mainly used to predict the magnitude of milling force.Ding [16] and Lei [17] used orthogonal tests to obtain the milling force coefficients and established an empirical model of milling force.Bergmann [18] developed a nonlinear empirical model of milling force and improved the prediction of stability limits achieved by parameterizing a milling force model for stability analysis.Zhao [19] established an empirical model of milling force for ultrasonic vibration milling of titanium alloy and explored the effects of different machining parameters on milling force.Liang [20] proposed an empirical model of cutting forces in micro-end-milling operations that takes into account the parameters of feed per tooth per tool radius, the true trochoidal nature of the tool edge path, and the chip thickness.In fact, the empirical model of milling force relies on a large amount of experimental data and has obvious limitations, which cannot reveal the dynamic characteristics and mechanism of cutting and can only be used to calculate the average milling force, rather than accurately predict the instantaneous cutting force during the milling process.
A physical analytical model of milling force is a model that describes the cutting forces generated during the milling process with the help of mathematical analytical methods.The physical analytical model of milling force is usually based on the basic theories of machining principles and material mechanics.Jiang [21] proposed that the projected area of the milling shear surface in three directions determines the milling force in three directions, and established an analytical model of milling force based on the shear surface.Madajewski [22] presented a computational method based on a combination of finite element analysis and classical analytical methods for predicting the cutting force components in face milling processes.Wang [23] discretized the end milling cutter along the circumferential direction and established an analytical cutting force model for the end milling cutter and simplified analytical indexes for predicting and estimating cutting force fluctuations.Kao [24] developed a mathematical model of the cutting forces of a ball-end milling cutter that includes tangential, radial, and axial cutting force components and characterized the relationship between average cutting force and feed per groove as a linear function.Kazuki [25] and Li [26] modeled the milling force of a ball-end milling cutter considering changes in the microelement friction angle, shear angle, and shear stress constraints using microelement diagonal shear theory.Falta [27] and Dong [28] used an analytical modeling method for predicting the static milling force of a ball-end milling cutter by establishing a tool-workpiece engagement region.The removal mechanism of materials can be deeply analyzed through the physical analytical model; however, the physical analytical model is closely related to material constitutive parameters, which are difficult to obtain.Therefore, the physical analytical model is not convenient for practical application, and its modeling process is complex and cumbersome.
A finite element model of milling force is a model that predicts milling forces and contact forces during the cutting process by discretizing the material regions being milled into small finite element regions.Soo [29,30] introduced a sophisticated finite element model developed for simulating the intricate process of high-speed ball nose end milling, specifically applied to Inconel 718, while successfully correlating predictions with experimental data.Filho [31] delved into the challenges posed by the low machinability of Inconel 718, and distinct plastic deformation levels between up and down milling, a factor attributed to the work-hardening behavior of Inconel 718, was revealed with finite element analysis and validated using experimental data.Hu [32] modeled the thermal-mechanical coupling property of ultrasonic torsional vibration-assisted micromilling (UTVAM) using ABAQUS software.Pratap [33] used the Johnson-Cook constitutive equation to develop a finite element model for microchannel micro-end milling of titanium alloys using ABAQUS that takes into account the strain, strain rate, and temperature effects on material properties.Molaiekiya [34] employed various methods, including tool life analysis, wear pattern investigation, and finite element analysis, to understand the behavior of SiAlON ceramic tools in high-speed face milling of Inconel 718.Deng [35] derived the milling mechanism of superalloy GH4169 under different processing parameters using the finite element method, and the results of the finite element analysis were successfully verified through milling experiment.Azeem [36] conducted molecular dynamics simulations to the investigate oxide-dislocation interaction mechanisms of oxide-dispersed strengthened (ODS) steel, and the interaction mechanisms between dislocations and oxides in ODS steel models were investigated.Azeem [37] utilized molecular dynamics (MD) simulations to gain insight into the irradiation effects on oxide-dispersed strengthened steel (ODSS) alloys, shedding light on the intricate microstructural changes and radiation resistance of these materials.Through using ABAQUS finite element software, Zhang [38] elucidated the impact of Al content on milling temperature, cutting forces, residual stress, and surface roughness during the milling process of CoCrFeNiAl X high-entropy alloys.Bouzakis [39] introduced an innovative methodology for assessing the effectiveness of coated tools in milling challenging aerospace alloys; using analytical methods, impact tests, and finite element simulations, the study established a strong correlation between coating impact resistance and cutting performance under varying conditions.Kim [40] delved into the realm of the induction-assisted milling (IAM) of Inconel 718 and through the utilization of the Taguchi method and finite element analysis, the research identified optimal machining conditions for both coated and uncoated carbide tools.The finite element model of milling force can simulate the real milling process, consider the contact and interaction between the tool and the workpiece, and is able to predict the mechanical properties and deformation behavior of material.However, the use of different finite element models will lead to different accuracy results.Therefore, it is recommended that a higher precision, more targeted finite element model be chosen to simulate the milling process accurately and intuitively.
Chattering is an unstable state of motion in metal cutting processing, mainly caused by cutting force, cutting heat, and material deformation.Inconel 718 is difficult-to-machine material, and its high hardness and heat resistance make it more susceptible to chattering during the milling process.At present, the theoretical analysis methods for milling stability mainly include frequency domain analysis and time domain analysis.
Frequency domain analysis is a method used to evaluate and solve the chattering problem in the milling process by converting the dynamic response of the milling system to the frequency domain.The frequency domain analysis method can be further divided into the zero-order frequency domain method and the multifrequency method.Altintaş and Budak [41] proposed a zero-order frequency domain method of milling stability in which the oriented dynamic milling force coefficient matrix is transformed into the frequency domain by Fourier transform to solve for the stability of the system.Based on this, the critical axial depth of cut without causing chattering is determined.Merdol [42] considered the higher-order expansion terms of the Fourier series of the oriented dynamic milling force coefficient matrix and developed a multifrequency method that can be used in cases of different radial depths of cut.Time domain analysis is another method used to evaluate and solve the chattering problem in the milling process by observing and analyzing the time response.The time domain analysis method can be further divided into the fully discrete and semi-discrete methods.Ding [43,44] verified the accuracy of the method using the single-degree-offreedom and two-degree-of-freedom milling models and then proposed the second-order full-discrete method.Insperger [45][46][47] proposed a semi-discrete method to solve a system of time-delay differential equations containing a matrix of periodic coefficients.Subsequently, the method is improved so that it cannot only solve the stability limit of a single-degree-of-freedom vibration system but also a two-degree-of-freedom vibration system, and a first-order semi-discrete method with better convergence is proposed.Based on the semi-discrete method, Sun [48] established a stability lobe diagram of the rotary ultrasonic milling of titanium alloy by defining the ultrasonic function angle, and the results show that the rotary ultrasonic milling can suppress the machining vibration and improve the stability of the milling of titanium alloy.Gajera [49] proposed a study focusing on optimizing selective laser melting parameters for Invar-36 material, considering hardness and surface roughness as key performance indicators.
In summary, scholars have conducted research and a series of studies in the field of metal cutting and have achieved results in the modeling of milling forces and milling stability analysis.However, studies that concentrate on the milling mechanism and chattering stability of nickel-based superalloy Inconel 718 are relatively rare, and there exists no thorough and systematic research on the optimization of the processing parameters of nickel-based superalloy Inconel 718.Therefore, it is imperative to carry out in-depth research on the influence of different processing parameters on the milling of nickel-based superalloy Inconel 718 and the chattering stability during the milling process, as well as the optimization of the process parameters.In this paper, a milling force model of a ball-end milling cutter is established, and the milling force coefficients are deduced.According to the material failure criteria, a finite element model of both the tool and workpiece is created based on a finite element analysis, and the influence of different processing parameters on the milling of Inconel 718 is analyzed.Through applying the zero-order frequency domain method, a dynamic equation of the milling of Inconel 718 is established, and based on the modal experiment, stability lobe diagrams are drawn, which can be used as an essential reference for the selection of stable cutting parameters.Eventually, single factor experiments and orthogonal experiments are carried out to verify the correctness of the milling force model and results of the finite element analysis.After the reliability of the milling force model and results of the finite element analysis is demonstrated, the fmincon optimization algorithm is used to optimize the processing parameters of Inconel 718, and the best processing parameters under different processing requirements are determined.
This article presents a novel exploration of the machining of nickel-based superalloy Inconel 718, addressing a significant gap in the existing literature.While studies in metal cutting and machining are prevalent, there is a distinct lack of comprehensive research concerning the milling mechanism and chattering stability of Inconel 718.A distinct milling force model is formulated for ball-end milling cutters applied specifically to Inconel 718.This model enables the precise comprehension of the forces involved during machining, enhancing insights into the effects of different parameters.Incorporating finite element analysis, the intricate relationship among various processing parameters and the material behavior of Inconel 718 during machining are explored.A novel aspect of this work is the meticulous examination of the chattering stability during the machining process.This exploration is paramount to achieving stable machining outcomes for challenging materials like Inconel 718 so that the appropriate processing parameters can be selected to prevent chattering.Lastly, this study introduces a multi-objective optimization strategy, where conflicting objectives, such as surface roughness and material removal rate, are concurrently considered.Based on the multi-objective optimization algorithm, the best machining parameters are determined, which point to an improvement in the machining efficiency and quality of Inconel 718, which has not been addressed by current studies.

Modeling of Milling Forces of Ball-End Milling Cutter
The milling state of the helix cutting edge of a ball-end milling cutter varies with different cutting positions.Thus, it is essential to develop a geometric model of the ball-end milling cutter to investigate the change in the milling force when the cutting edge interacts with the workpiece.The schematic of the proposed model is illustrated in Figure 1.
According to the literature [24], neglecting the tip radius of the ball-end milling cutter, the cutting force on any microelement of the cutting edge can be described using Equation (1).   dF t,j φ j (z) = K te dS φ j (z) + K tc h j φ j (z) db dF r,j φ j (z) = K re dS φ j (z) + K rc h j φ j (z) db dF a,j φ j (z) = K ae dS φ j (z) + K ac h j φ j (z) db According to the literature [24], neglecting the tip radius of the ball-end m ter, the cutting force on any microelement of the cutting edge can be described u tion (1).
, t j j te j tc j j r j j re j rc j j a j j ae j ac j j The instantaneous chip thickness at the immersion angle ( ) ( ) = sin sin The instantaneous edge length of the cutting segment The instantaneous chip width db can be expressed as: By transforming the cutting force microelement from coordinate system The instantaneous chip thickness at the immersion angle h j φ j (z) can be expressed as: The instantaneous edge length of the cutting segment dS φ j (z) can be expressed as: The instantaneous chip width db can be expressed as: By transforming the cutting force microelement from coordinate system O − tra to O − f nz and integrating it, the cutting force of the milling cutter along the feed, normal, and axial directions can be obtained as: Subsequently, the average milling force per rotation of the ball-end milling cutter can be calculated using the Equation (6): Equation ( 7) is finally obtained by simplifying the equation of the microelement of the milling force and substituting it into Equation ( 6): The parameters C 1 − C 16 in Equation (7) were provided in Reference [24].

Recognition Experiment of Milling Force Coefficients
In order to accurately determine the milling force generated during milling processing, it is necessary to determine the milling force coefficients K te , K re , K ae , K tc , K rc , and K ac in Equation (7); therefore, recognition experiments of the milling force coefficients are required.
For this experiment, the 1B240-0800-XA 1630 carbide ball-end milling cutter (developed by Sandvik in Germany) was chosen, the diameter of the cutter was 8 mm, the overall length of the cutter was 100 m, the depth of cut was set to 0.5 mm, and the rotating speed of spindle was set to 1000 rpm.To facilitate the calculation of the milling force coefficients, slot milling was used as the milling method, and the selected machining parameters and the measured average milling force data are recorded in Table 1.Based on the data in Table 1, the milling force coefficients were obtained using the coefficient recognition equation of the milling force (Equation ( 8)).
The coefficients of the milling force calculating through Equation ( 8) are represented in Table 2.

Finite Element Analysis of Milling of Inconel 718
In this study, a finite element model is established, as shown in Figure 2, including the workpiece and the tool.In the following, the temperature, stress, and milling force in the milling process are calculated.Finally, the influence of different machining parameters on the machining condition of Inconel 718 is evaluated through simulation results.In this study, a finite element model is established, as shown in Figu the workpiece and the tool.In the following, the temperature, stress, and m the milling process are calculated.Finally, the influence of different machin on the machining condition of Inconel 718 is evaluated through simulation

Methodology of Finite Element Analysis
Nickel-based superalloy Inconel 718 was chosen as the research obje element analysis, and AdvantEdge (version 7.1, developed by Third W America) was utilized as the analysis software.
To study the effect of different processing parameters, such as spind per tooth, and radial and axial depths of cut, on the temperature, stress, an

Methodology of Finite Element Analysis
Nickel-based superalloy Inconel 718 was chosen as the research object of the finite element analysis, and AdvantEdge (version 7.1, developed by Third Wave System in America) was utilized as the analysis software.
To study the effect of different processing parameters, such as spindle speed, feed per tooth, and radial and axial depths of cut, on the temperature, stress, and cutting force generated during the milling process of Inconel 718, the parameters shown in Table 3 were used for the finite element analysis.The constitutive model is a critical prerequisite of the accuracy and reliability of the finite element analysis.The behavior of materials and results of an engineering problem can be better provided by the right chosen constitutive model [50].The constitutive models applicable to metal cutting are mainly the power law and Johnson-Cook models [51].The power law constitutive model is usually applicable to the analysis of the deformation behavior of materials at high temperatures and high strain rates, while the Johnson-Cook constitutive model is more suitable for the analysis of the deformation behavior of materials at low to medium temperatures and low to medium strain rates.Therefore, the deformation behavior of Inconel 718 at high temperatures and high strain rates can be described using the power law constitutive model.The mathematical equation of the power law model can be expressed as: In Equation ( 9), the strain enhancement function g(ε p ) can be further expressed as: The strain enhancement parameters of Inconel 718 are shown in Table 4.In Equation ( 9), the thermal softening function Θ(T) can be expressed as: The thermal softening parameters of Inconel 718 are shown in Table 5.In Equation ( 9), the strain rate effect function Γ .
ε can be expressed as: The strain rate parameters of Inconel 718 are shown in Table 6.The damage model for Inconel 718 can be determined using the following equation: In Equation ( 13), the parameter ∆ε p i can be expressed as: In Equation ( 14), the parameters d 0 − d 5 are the strain values deduced using the tensile test.

Temperature Field Analysis of Milling Process
Inconel 718 is a difficult-to-machine material, which has low thermal conductivity.The phenomenon of thermal mechanical coupling could occur during the machining process, which results in residual stress and reduces the surface quality.Therefore, to effectively control tool wear and improve the surface integrity of the workpiece, the milling temperatures must be systematically studied and analyzed.
The fitted curves of the maximum temperature that occurred at the tool-workpiece contact region at different spindle speeds, feed per tooth, and radial and axial depths of cut are plotted according to the results of the finite element analysis, and these curves are shown in Figure 3a-d.
As shown in Figure 3, the temperature increases as the spindle speed rises, and this is mainly due to the fact that the friction phenomenon between the front (or back) of the tool and workpiece becomes more intense as the spindle speed increases, and this increasing trend is especially conspicuous near the cutting edge.Moreover, with the increase in the radial depth of cut, axial depth of cut, and feed per tooth, the temperature also increases, because increases in all of these parameters cause the tool to cut more material per unit of time, so the heat continues to accumulate and more heat diffuses into the tool.

Stress Field Analysis of Milling Process
Because of the high hardness and strength properties of Inconel 718, the milling process of this material is susceptible to stress concentration and crack formation.Stress field simulation enables the visualization of the stress distribution on both the tool and workpiece surfaces during machining, facilitating a deeper understanding of the physical phenomena involved.Consequently, the impact of various machining parameters on stress can be predicted to mitigate stress concentration and prevent crack formation.
To investigate the impact of processing parameters on stresses during machining, the stress values in the top 10% of elements with the highest stresses on the tool are averaged and plotted as a function of time.Fitted curves representing stress variations at different spindle speeds, feed per tooth, and radial and axial depths of cut are presented in Figure 4a-d.
From Figure 4, it is evident that an increase in the spindle speed does not significantly affect the stress, as the material volume cut by the tool remains unchanged despite the enlargement of the spindle speed.Conversely, with an increment in feed per tooth, radial depth of cut, and axial depth of cut, there is a simultaneous increase in stress.This observation can be attributed to the larger volume of material being removed because of changes in these parameters.Therefore, more force is required to remove the material from the surface of the workpiece, and the stress acting on the tool tends to increase.As shown in Figure 3, the temperature increases as the spindle speed rises, and this is mainly due to the fact that the friction phenomenon between the front (or back) of the tool and workpiece becomes more intense as the spindle speed increases, and this increasing trend is especially conspicuous near the cutting edge.Moreover, with the increase in the radial depth of cut, axial depth of cut, and feed per tooth, the temperature also increases, because increases in all of these parameters cause the tool to cut more material per unit of time, so the heat continues to accumulate and more heat diffuses into the tool.

Force Analysis of Milling Process
Given that the magnitude of the milling force directly influences the wear and fracture condition of the tool, a comprehensive analysis of the milling force generated during the Inconel 718 milling process becomes crucial.
piece surfaces during machining, facilitating a deeper understanding of the physical phe-nomena involved.Consequently, the impact of various machining parameters on stress can be predicted to mitigate stress concentration and prevent crack formation.
To investigate the impact of processing parameters on stresses during machining, the stress values in the top 10% of elements with the highest stresses on the tool are averaged and plotted as a function of time.Fitted curves representing stress variations at different spindle speeds, feed per tooth, and radial and axial depths of cut are presented in Figure 4a-d  In order to better study the influence trend of different machining parameters on milling force, the spindle speed was set to 1000 rpm, 1300 rpm, and 1600 rpm, and the other parameters were kept unchanged.The milling force data of each axis are shown in Figure 5.
Given that the magnitude of the milling force directly influences the wear and fracture condition of the tool, a comprehensive analysis of the milling force generated during the Inconel 718 milling process becomes crucial.
In order to better study the influence trend of different machining parameters on milling force, the spindle speed was set to 1000 rpm, 1300 rpm, and 1600 rpm, and the other parameters were kept unchanged.The milling force data of each axis are shown in Figure 5. Figure 5 illustrates that the effect of different spindle speeds on the magnitude of the milling force for each axis is not significant.This can be attributed to the nearly constant volume of material being removed, regardless of the spindle speed variation.Actually, Inconel 718 is known for its high toughness and work-hardening characteristics.As the spindle speed increases, the cutting edge may experience slightly different interactions with the material due to the variations in chip thickness and heat generation.However, the material's inherent properties, including its resistance to deformation and heat, might minimize the impact of spindle speed changes on the cutting forces.
The feed per tooth was varied at values of 0.2 mm, 0.3 mm, 0.4 mm, and 0.5 mm while keeping the remaining parameters unchanged.The corresponding milling force data for each axis are presented in Figure 6.
Figure 6 demonstrates a clear and consistent increasing trend in the milling forces for each axis with the increase in the feed per tooth.This observation can be attributed to the larger volume of workpiece material being removed by the cutter as the feed per tooth increases, because Inconel 718 is known for its high strength, toughness, and resistance to deformation.As the feed per tooth increases, more material is being removed with each tooth engagement.This requires the cutting tool to exert greater force to cut through the tough material.When machining Inconel 718, heat generated at the cutting interface can influence the material's behavior.Higher feeds can lead to more heat being generated due to increased friction, and this can affect the material's strength and work-hardening tendencies, potentially increasing the forces needed to cut through the material.The chip formation also accounts for this phenomenon of increasing milling force, since the increase in feed per tooth can lead to the generation of larger chips or greater chip thickness.Inconel 718 tends to form short, tightly curled chips that can exert backpressure on the cutting tool, causing additional cutting forces.Figure 5 illustrates that the effect of different spindle speeds on the magnitude of the milling force for each axis is not significant.This can be attributed to the nearly constant volume of material being removed, regardless of the spindle speed variation.Actually, Inconel 718 is known for its high toughness and work-hardening characteristics.As the spindle speed increases, the cutting edge may experience slightly different interactions with the material due to the variations in chip thickness and heat generation.However, the material's inherent properties, including its resistance to deformation and heat, might minimize the impact of spindle speed changes on the cutting forces.
The feed per tooth was varied at values of 0.2 mm, 0.3 mm, 0.4 mm, and 0.5 mm while keeping the remaining parameters unchanged.The corresponding milling force data for each axis are presented in Figure 6.The radial depth of cut was varied at values of 2 mm, 3 mm, 4 mm, and 5 mm while keeping the remaining parameters unchanged.The corresponding milling force data for each axis are presented in Figure 7.
Figure 7 reveals that the milling forces in the x-axis, y-axis, and z-axis exhibit a notable increase with the rise in the radial depth of cut.A larger radial depth of cut means that a larger portion of the cutting edge is engaged with the material.This results in a larger contact area between the tool and the workpiece.Therefore, the tool has to remove material from a wider area, requiring more force to achieve this.Inconel 718 has high strength and hardness, which means that cutting through it requires more force.In addition, as the cutting tool engages with the material, it generates heat due to the friction.This heat can cause work hardening in the material, making it even tougher to cut.With a higher radial depth of cut, more heat is generated, leading to increased work hardening and subsequently higher cutting forces.Apart from this, Inconel 718 tends to produce shorter, more tightly curled chips because of its toughness and heat resistance.These chips can put additional stress on the cutting tool and increase the cutting forces.With a higher radial depth of cut, there is more material being removed and, consequently, more chips are generated, contributing to higher forces.In addition, cutting forces can also increase because of accelerated tool wear.The high temperatures generated during the machining process, especially with a higher radial depth of cut, can lead to faster tool wear.As the tool wears, it becomes less effective, requiring higher forces to achieve the same level of material removal.increases, because Inconel 718 is known for its high strength, toughness, and resistance to deformation.As the feed per tooth increases, more material is being removed with each tooth engagement.This requires the cutting tool to exert greater force to cut through the tough material.When machining Inconel 718, heat generated at the cutting interface can influence the material's behavior.Higher feeds can lead to more heat being generated due to increased friction, and this can affect the material's strength and work-hardening tendencies, potentially increasing the forces needed to cut through the material.The chip formation also accounts for this phenomenon of increasing milling force, since the increase in feed per tooth can lead to the generation of larger chips or greater chip thickness.Inconel 718 tends to form short, tightly curled chips that can exert backpressure on the cutting tool, causing additional cutting forces.The radial depth of cut was varied at values of 2 mm, 3 mm, 4 mm, and 5 mm while keeping the remaining parameters unchanged.The corresponding milling force data for each axis are presented in Figure 7.  Figure 7 reveals that the milling forces in the x-axis, y-axis, and z-axis exhibit a notable increase with the rise in the radial depth of cut.A larger radial depth of cut means that a larger portion of the cutting edge is engaged with the material.This results in a larger contact area between the tool and the workpiece.Therefore, the tool has to remove Similarly, the value of the axial depth of cut was set to 2 mm, 3 mm, 4 mm, and 5 mm, and the rest of the parameters were kept unchanged; the milling force data of each axis are shown in Figure 8.The data shown in Figure 8a unveil that the milling forces of the x-axis increase with an augmentation of the value of the axial depth of cut.This phenomenon can be explained by the abovementioned reasons.Figure 8b,c disclose that the milling forces in the y-and z-axes do not show a tangible rising tendency with the increase in the radial depth of cut.This observation can be attributed to the fact that a higher radial depth of cut primarily affects the x-axis (radial) direction forces, while the y and z forces are influenced by the tool's interaction with the material in their respective directions.If the tool's engagement angle remains relatively consistent, it might lead to stable forces.Furthermore, in the y-and z-directions, the material might flow more smoothly as the tool progresses through the cut.This controlled material flow can lead to more predictable forces, even with variations in the radial depth of cut.
tool wears, it becomes less effective, requiring higher forces to achieve the same level of material removal.
Similarly, the value of the axial depth of cut was set to 2 mm, 3 mm, 4 mm, and 5 mm, and the rest of the parameters were kept unchanged; the milling force data of each axis are shown in Figure 8.The data shown in Figure 8a unveil that the milling forces of the xaxis increase with an augmentation of the value of the axial depth of cut.This phenomenon can be explained by the abovementioned reasons.Figure 8b,c disclose that the milling forces in the y-and z-axes do not show a tangible rising tendency with the increase in the radial depth of cut.This observation can be attributed to the fact that a higher radial depth of cut primarily affects the x-axis (radial) direction forces, while the y and z forces are influenced by the tool's interaction with the material in their respective directions.If the tool's engagement angle remains relatively consistent, it might lead to stable forces.Furthermore, in the y-and z-directions, the material might flow more smoothly as the tool progresses through the cut.This controlled material flow can lead to more predictable forces, even with variations in the radial depth of cut.

Establishment of Dynamic Equation for Milling Inconel 718
Inconel 718 belongs to difficult-to-machine materials and is susceptible to the phenomenon of chattering during the milling process.Chattering in the milling process will lead to a positive feedback effect between the cutting force and vibration, which will result in degradation of the milling quality, cutting edge wear, tool fracture, and other problems.Through a stability analysis of the chattering, the stability condition of the milling process can be predicted, and whether chattering will occur during the milling process can also be determined.Therefore, it is necessary to carry out a chattering stability analysis of Inconel 718 to improve the milling efficiency and machining quality, as well as to reduce the production cost.
As shown in Figure 9, the machine-tool system is simplified as a vibration system with two degrees of freedom, where the vibrations in the x-and y-directions can be described using the kinetic differential equations:

Analysis of Milling Stability 4.1. Establishment of Dynamic Equation for Milling Inconel 718
Inconel 718 belongs to difficult-to-machine materials and is susceptible to the phenomenon of chattering during the milling process.Chattering in the milling process will lead to a positive feedback effect between the cutting force and vibration, which will result in degradation of the milling quality, cutting edge wear, tool fracture, and other problems.Through a stability analysis of the chattering, the stability condition of the milling process can be predicted, and whether chattering will occur during the milling process can also be determined.Therefore, it is necessary to carry out a chattering stability analysis of Inconel 718 to improve the milling efficiency and machining quality, as well as to reduce the production cost.
As shown in Figure 9, the machine-tool system is simplified as a vibration system with two degrees of freedom, where the vibrations in the x-and y-directions can be described using the kinetic differential equations: Inconel 718 belongs to difficult-to-machine materials and is susceptible to the phenomenon of chattering during the milling process.Chattering in the milling process will lead to a positive feedback effect between the cutting force and vibration, which will result in degradation of the milling quality, cutting edge wear, tool fracture, and other problems.Through a stability analysis of the chattering, the stability condition of the milling process can be predicted, and whether chattering will occur during the milling process can also be determined.Therefore, it is necessary to carry out a chattering stability analysis of Inconel 718 to improve the milling efficiency and machining quality, as well as to reduce the production cost.
As shown in Figure 9, the machine-tool system is simplified as a vibration system with two degrees of freedom, where the vibrations in the x-and y-directions can be described using the kinetic differential equations: In the milling process, as shown in Figure 10, the current cutting process will be affected by the previous cutting process, since the previous machining process leaves vibrations on the machined surface, and when the tool cuts the surface again, this can lead to In the milling process, as shown in Figure 10, the current cutting process will be affected by the previous cutting process, since the previous machining process leaves vibrations on the machined surface, and when the tool cuts the surface again, this can lead to changes in the cutting thickness.Such changes in the cutting thickness can cause fluctuations in the cutting force, leading to further vibrations in the system.The total cutting thickness generated by the superposition of the static cutting thickness and the dynamic cutting thickness caused by the tool vibration can be expressed as: In Equation ( 16), the parameter ( ) j g φ is a unit step function, which is expresse below: In Equation ( 16), the parameter g φ j is a unit step function, which is expressed below: Therefore, the total cutting force on the cutter can be expressed as: In Equation ( 18), the parameters α xx , α xy , α yx , and α yy are the average directional force coefficient: Thus, the equation of the dynamic milling force can be expressed as: Equation ( 20) can be further simplified as: In Equation ( 21), the parameter [A 0 ] can be expressed as: The stability of the system is solved using the zero-order frequency domain method, such that the frequency response transfer function matrix of the tool-workpiece contact zone is:

Methodology of the Modal Experiment
To solve the dynamic equation, it is necessary to obtain the modal parameters of the tool, including the natural frequency, damping ratio, and vibration mode.The approach of the force hammer excitation was used in this experiment to obtain accurate modal parameters.
The equipment utilized for conducting the modal experiment included a 1B240-0800-XA carbide ball-end milling cutter (developed by Sandvik in Germany), a three-way acceleration sensor (developed by Kistler in Switzerland), a TST5928 Dynamic Strain Gauge, and a computer equipped with the TST5928 Distributed Dynamic Signal Test and Analysis System.A schematic of the hardware installation of the system for the modal experiment is shown in Figure 11, and photos of the connection of the experimental equipment is shown in Figure 12.
XA carbide ball-end milling cutter (developed by Sandvik in Germany), a three-way acceleration sensor (developed by Kistler in Switzerland), a TST5928 Dynamic Strain Gauge, and a computer equipped with the TST5928 Distributed Dynamic Signal Test and Analysis System.A schematic of the hardware installation of the system for the modal experiment is shown in Figure 11, and photos of the connection of the experimental equipment is shown in Figure 12.  Figure 13 shows the distribution of the measurement points of the tool, where a force hammer was used to excite the tool at different measurement points, and the response of the tool was measured.Since the stiffness of the tool in the z-direction was much greater than in the x-and y-directions, only modal experiments in the x-and y-directions were performed.The equipment utilized for conducting the modal experiment included a 1B240-0800 XA carbide ball-end milling cutter (developed by Sandvik in Germany), a three-way acce eration sensor (developed by Kistler in Switzerland), a TST5928 Dynamic Strain Gauge, an a computer equipped with the TST5928 Distributed Dynamic Signal Test and Analysis Sys tem.A schematic of the hardware installation of the system for the modal experiment i shown in Figure 11, and photos of the connection of the experimental equipment is show in Figure 12.  Figure 13 shows the distribution of the measurement points of the tool, where a forc hammer was used to excite the tool at different measurement points, and the response o the tool was measured.Since the stiffness of the tool in the z-direction was much greate than in the x-and y-directions, only modal experiments in the x-and y-directions wer performed.Figure 13 shows the distribution of the measurement points of the tool, where a force hammer was used to excite the tool at different measurement points, and the response of the tool was measured.Since the stiffness of the tool in the z-direction was much greater than in the x-and y-directions, only modal experiments in the x-and y-directions were performed.

The Result of Modal Experiment
The results of modal experiments are recorded in Table 7.

Modal Order
Natural Frequency (Hz) Damping Ratio (%) First order of x-direction 935.06 0.049

The Result of Modal Experiment
The results of modal experiments are recorded in Table 7.

The Result of Modal Experiment
The results of modal experiments are recorded in Table 7.   Combining the results of Figures 14 and 15, as well as Table 7, it is possible to derive that the modal shapes corresponding to the first-order natural frequency in the x-and ydirections of the tool are both bent and deformed in the middle part of the tool.The modal shapes corresponding to the second-order natural frequency in the x-and y-directions of the tool are both bent and deformed in the upper part of the tool, and the degree of deformation in the y-direction is slightly larger than that in the x-direction.The bending defor- Combining the results of Figures 14 and 15, as well as Table 7, it is possible to derive that the modal shapes corresponding to the first-order natural frequency in the x-and y-directions of the tool are both bent and deformed in the middle part of the tool.The modal shapes corresponding to the second-order natural frequency in the x-and y-directions of the tool are both bent and deformed in the upper part of the tool, and the degree of deformation in the y-direction is slightly larger than that in the x-direction.The bending deformation of the modal shapes corresponding to third-order natural frequency in the x-direction occurs in the middle and upper parts of the tool, and the deformation of the bottom of the tool is small, while the bending deformation of the modal shapes corresponding to the third-order natural frequency in the y-direction occurs in the lower part of the tool.
As it is necessary to select suitable processing parameters to avoid chattering during the milling process, therefore, there is a need to investigate the effect of the tool radius, radial depth of cut, natural frequency ω n , and damping ratio ζ on the cutting stability.The stability lobes are plotted for different tool radii, radial depths of cut, natural frequencies ω n , and damping ratios ζ, as shown in Figures 16-19.Combining the results of Figures 14 and 15, as well as Table 7, it is possible to derive that the modal shapes corresponding to the first-order natural frequency in the x-and ydirections of the tool are both bent and deformed in the middle part of the tool.The modal shapes corresponding to the second-order natural frequency in the x-and y-directions of the tool are both bent and deformed in the upper part of the tool, and the degree of deformation in the y-direction is slightly larger than that in the x-direction.The bending deformation of the modal shapes corresponding to third-order natural frequency in the x-direction occurs in the middle and upper parts of the tool, and the deformation of the bottom of the tool is small, while the bending deformation of the modal shapes corresponding to the third-order natural frequency in the y-direction occurs in the lower part of the tool.
As it is necessary to select suitable processing parameters to avoid chattering during the milling process, therefore, there is a need to investigate the effect of the tool radius,       Figure 16 illustrates the influence of the tool radius on stability.As the tool radius increases, the system stability rises.This is because a larger tool radius provides greater stability by reducing the tendency for the tool to vibrate or chatter during machining.Inconel 718′s tough nature can amplify vibration problems, so using a larger tool radius helps dampen vibrations and leads to smoother cutting.Increasing the tool radius effectively decreases the depth of cut for a given engagement.This leads to a reduction in cutting forces as the material removal becomes more gradual, and the vibration can also be alleviated during the milling process.
Figure 17 demonstrates the influence of the radial depth of cut on stability, and it can be seen that the stability of the system decreases as the radial depth of cut increases, because a higher radial depth of cut leads to a larger chip thickness and engagement of the cutting tool with the workpiece.This generates higher cutting forces that can lead to tool deflection, chatter, and vibrations, reducing the stability of the machining process.With larger radial depths of cut, the likelihood of chattering and vibrations increases, and it can be especially problematic in milling hard-to-machine materials like Inconel 718.Additionally, larger radial depths of cut can lead to the formation of larger chips.Proper chip evacuation becomes challenging, and chip clogging can cause sudden force variations and affect system stability.Figure 16 illustrates the influence of the tool radius on stability.As the tool radius increases, the system stability rises.This is because a larger tool radius provides greater stability by reducing the tendency for the tool to vibrate or chatter during machining.Inconel 718 s tough nature can amplify vibration problems, so using a larger tool radius helps dampen vibrations and leads to smoother cutting.Increasing the tool radius effectively decreases the depth of cut for a given engagement.This leads to a reduction in cutting forces as the material removal becomes more gradual, and the vibration can also be alleviated during the milling process.
Figure 17 demonstrates the influence of the radial depth of cut on stability, and it can be seen that the stability of the system decreases as the radial depth of cut increases, because a higher radial depth of cut leads to a larger chip thickness and engagement of the cutting tool with the workpiece.This generates higher cutting forces that can lead to tool deflection, chatter, and vibrations, reducing the stability of the machining process.With larger radial depths of cut, the likelihood of chattering and vibrations increases, and it can be especially problematic in milling hard-to-machine materials like Inconel 718.Additionally, larger radial depths of cut can lead to the formation of larger chips.Proper chip evacuation becomes challenging, and chip clogging can cause sudden force variations and affect system stability.
Figure 18 illustrates the influence of the natural frequency of the machine-tool system on stability, and it can be deduced that as the natural frequency increases, there is no significant change in the stability diagram, but the curve in the stability diagram is shifted to the right as the natural frequency increases.When the stability curve shifts to the right as the natural frequency increases, this suggests that the system's inherent ability to resist vibrations or chatter improves.In other words, higher natural frequencies indicate that the system is less prone to self-excited vibrations that lead to chatter.This is because a higher natural frequency often implies that the system has better inherent damping characteristics.With increased damping, the system is more effective at absorbing and dissipating the energy generated during machining, leading to reduced vibrations and improved stability.
Figure 19 demonstrates the influence of the damping ratio on stability, and it can be concluded that an increasing damping ratio of the machine-tool system can improve the system stability.With increased damping, the system is more effective at absorbing and dissipating the energy generated during machining, leading to reduced vibrations and improved stability, because a higher damping ratio effectively increases the effective rigidity of the machine-tool system.This added rigidity counters the deflection caused by cutting forces and enhances the overall stability of the machining process.

Methodology of Milling Experiment
To verify the correctness of the established milling force model and finite element analysis, a series of experiments including single factor experiments and orthogonal experiments are designed.After, the milling force measured from the experiments with that calculated from the theoretical model, as well as the milling forces obtained from the finite element analysis, are compared.
A DMC635V vertical machining center (developed by DMG in Germany), 1B240-0800-XA 1630 carbide ball-end milling cutter (developed by Sandvik in Germany), 9527B force gauge (developed by Kistler in Switzerland), and 5073A charge amplifier and acquisition card, as well as Dynoware computer analysis software (version 3.2.5.0), were selected to build the experimental platform.A schematic diagram of the connection of experimental equipment is shown in Figure 20, and a photo of the device's connections is shown in Figure 21.
Materials 2023, 16, x FOR PEER REVIEW 23 of 39 improved stability, because a higher damping ratio effectively increases the effective rigidity of the machine-tool system.This added rigidity counters the deflection caused by cutting forces and enhances the overall stability of the machining process.

Methodology of Milling Experiment
To verify the correctness of the established milling force model and finite element analysis, a series of experiments including single factor experiments and orthogonal experiments are designed.After, the milling force measured from the experiments with that calculated from the theoretical model, as well as the milling forces obtained from the finite element analysis, are compared.
A DMC635V vertical machining center (developed by DMG in Germany), 1B240-0800-XA 1630 carbide ball-end milling cutter (developed by Sandvik in Germany), 9527B force gauge (developed by Kistler in Switzerland), and 5073A charge amplifier and acquisition card, as well as Dynoware computer analysis software (version 3.2.5.0), were selected to build the experimental platform.A schematic diagram of the connection of experimental equipment is shown in Figure 20, and a photo of the device's connections is shown in Figure 21.The single factor experiments and orthogonal experiments were designed as shown in Tables 8 and 9.The single factor experiments and orthogonal experiments were designed as shown in Tables 8 and 9.

Results of the Analysis of the Single Factor Experiment
Based on the above-proposed milling force model of the ball-end milling cutter, the MATLAB program, which is illustrated in Appendix A, was used to calculate the theoretical milling force.The results of the theoretical analysis were compared to the experimentally measured results to verify the correctness of the proposed milling force model.The same manufacturing parameters were set for both the theoretical and experimental analyses, namely, the spindle speed was 1000 rpm, the axial depth of cut was 0.4 mm, and the feed per tooth was 0.04 mm/z.The results of the comparative analysis of the theoretical and experimental values of milling forces of the x-, y-, and z-axes are shown in Figure 22.Observing the milling force values for each axis, it is evident that both the theoretical and experimental results exhibit a similar change trend, validating the accuracy of the proposed theoretical model for calculating the ball-end milling cutter's milling force.Although a slight degree of error is present in the figures, they remain within an acceptable range.The milling forces measured through milling experiment are presented in Figure 23a, and the milling forces calculated using the theoretical milling force model are presented in Figure 23b.The milling forces of each axis obtained from the theoretical model align closely with the experimental results, validating the accuracy and correctness of the proposed model.To investigate the influence of the spindle speed on the milling of Inconel 718, the parameter of spindle speed was varied at 1000 rpm, 1200 rpm, 1400 rpm, 1600 rpm, and 1800 rpm while keeping the other machining parameters unchanged.The axial depth of cut was set to 0.4 mm, and the feed per tooth was set to 0.02 mm/z.The experimental results of the milling force obtained under these conditions are presented in Figure 24.In Figure 24, the average milling force refers to the force calculated by selecting a stable range of cutting area presented in Dynoware software and making average the milling force in The milling forces measured through milling experiment are presented in Figure 23a, and the milling forces calculated using the theoretical milling force model are presented in Figure 23b.The milling forces of each axis obtained from the theoretical model align closely with the experimental results, validating the accuracy and correctness of the proposed model.To investigate the influence of the spindle speed on the milling of Inconel 718, the parameter of spindle speed was varied at 1000 rpm, 1200 rpm, 1400 rpm, 1600 rpm, and 1800 rpm while keeping the other machining parameters unchanged.The axial depth of cut was set to 0.4 mm, and the feed per tooth was set to 0.02 mm/z.The experimental results of the milling force obtained under these conditions are presented in Figure 24.In Figure 24, the average milling force refers to the force calculated by selecting a stable range To investigate the influence of the spindle speed on the milling of Inconel 718, the parameter of spindle speed was varied at 1000 rpm, 1200 rpm, 1400 rpm, 1600 rpm, and 1800 rpm while keeping the other machining parameters unchanged.The axial depth of cut was set to 0.4 mm, and the feed per tooth was set to 0.02 mm/z.The experimental results of the milling force obtained under these conditions are presented in Figure 24.In Figure 24, the average milling force refers to the force calculated by selecting a stable range of cutting area presented in Dynoware software and making average the milling force in that area.Figure 24 illustrates that the average milling forces for the x-, y-, and z-axes show minimal changes with an increasing spindle speed, which aligns with the conclusions drawn from the finite element analysis.
To unveil the influence of the axial depth of cut on the milling of Inconel 718, the parameter of the axial depth of cut was varied at 0.2 mm, 0.3 mm, 0.4 mm, 0.5 mm, and 0.6 mm while keeping the other machining parameters constant.The spindle speed was set to 1000 rpm, and the feed per tooth was set to 0.02 mm/z.The experimental results of the milling force obtained under these conditions are presented in Figure 25.   Figure 24 illustrates that the average milling forces for the x-, y-, and z-axes show minimal changes with an increasing spindle speed, which aligns with the conclusions drawn from the finite element analysis.
To unveil the influence of the axial depth of cut on the milling of Inconel 718, the parameter of the axial depth of cut was varied at 0.2 mm, 0.3 mm, 0.4 mm, 0.5 mm, and 0.6 mm while keeping the other machining parameters constant.The spindle speed was set to 1000 rpm, and the feed per tooth was set to 0.02 mm/z.The experimental results of the milling force obtained under these conditions are presented in Figure 25. Figure 24 illustrates that the average milling forces for the x-, y-, and z-axes show minimal changes with an increasing spindle speed, which aligns with the conclusions drawn from the finite element analysis.
To unveil the influence of the axial depth of cut on the milling of Inconel 718, the parameter of the axial depth of cut was varied at 0.2 mm, 0.3 mm, 0.4 mm, 0.5 mm, and 0.6 mm while keeping the other machining parameters constant.The spindle speed was set to 1000 rpm, and the feed per tooth was set to 0.02 mm/z.The experimental results of the milling force obtained under these conditions are presented in Figure 25.   Figure 25 reveals that the average milling force for the x-, y-, and z-axes all demonstrated an increasing trend with the augmentation of the axial depth of cut.The z-axis exhibited the most significant increase, while the average milling forces in the x-and y-axes showed relatively smaller increments.These observations are consistent with the conclusions drawn from the finite element analysis.
To investigate the influence of the feed per tooth on the milling of Inconel 718, the parameter of feed per tooth was varied at 0.02 mm/z, 0.04 mm/z, 0.06 mm/z, 0.08 mm/z, and 0.10 mm/z while keeping the other machining parameters unchanged.The spindle speed was set to 1000 rpm, and the axial depth of cut was set to 0.4 mm.The experimental results of the milling force obtained under these conditions are presented in Figure 26.exhibited the most significant increase, while the average milling forces in the x-and yaxes showed relatively smaller increments.These observations are consistent with the conclusions drawn from the finite element analysis.
To investigate the influence of the feed per tooth on the milling of Inconel 718, the parameter of feed per tooth was varied at 0.02 mm/z, 0.04 mm/z, 0.06 mm/z, 0.08 mm/z, and 0.10 mm/z while keeping the other machining parameters unchanged.The spindle speed was set to 1000 rpm, and the axial depth of cut was set to 0.4 mm.The experimental results of the milling force obtained under these conditions are presented in Figure 26.It can be seen from Figure 26 that the average milling force for the x-, y-, and z-axes all reveal an increasing trend with the increase in feed per tooth, where the increase for the z-axis was the largest, and the increase for the x-and y-axes was relatively small, which is compatible with the conclusions drawn from the finite element analysis.

Results Analysis of the Orthogonal Experiments
The orthogonal experiments were analyzed based on the parameters shown in Table 9, and the orthogonal experimental results obtained are shown in Table 10.

Experiment
No.  It can be seen from Figure 26 that the average milling force for the x-, y-, and z-axes all reveal an increasing trend with the increase in feed per tooth, where the increase for the z-axis was the largest, and the increase for the x-and y-axes was relatively small, which is compatible with the conclusions drawn from the finite element analysis.

Results Analysis of the Orthogonal Experiments
The orthogonal experiments were analyzed based on the parameters shown in Table 9, and the orthogonal experimental results obtained are shown in Table 10.
An extreme variance analysis of the results shown in Table 10 lead to the following conclusions: 1.
Axial depth of cut was the most influential factor on the average milling force for the x-axis, followed by feed per tooth and spindle speed.

2.
Feed per tooth was the most influential factor on the average milling force for the y-axis, followed by axial depth of cut and spindle speed.

3.
Axial depth of cut was the most influential factor on the average milling force for the z-axis, followed by the feed per tooth and spindle speed.
Considering that the effect of the feed per tooth and axial depth of cut on the average milling force in the previous analysis was large, while the effect of the spindle speed was relatively small, different parameters of the axial depth of cut and feed per tooth were selected for the verification experiments.Groups 1-4, as shown in Table 9, were selected, the theoretical and finite element analysis value of the average milling force of these groups were calculated in order to compare with the value of the experiments shown in Table 10, and the results are recorded in Tables 11 and 12. Figure 27 illustrates the comparison of the theoretical, finite element analysis simulation, and experimental results of the average milling force for the x-, y-, and z-axes in four different groups.Tables 11 and 12, combined with Figure 27, reveal that the minimum relative error between the milling force calculated by the theoretical model and the experimental data is 6.48%, while the minimum relative error between the milling force obtained from finite element analysis and the experimental data is 4.24%.These results further demonstrate the reliability and validity of both the theoretical model and the finite element analysis.Figure 27 illustrates the comparison of the theoretical, finite element analysis simulation, and experimental results of the average milling force for the x-, y-, and z-axes in four different groups.Tables 11 and 12, combined with Figure 27, reveal that the minimum relative error between the milling force calculated by the theoretical model and the experimental data is 6.48%, while the minimum relative error between the milling force obtained from finite element analysis and the experimental data is 4.24%.These results further demonstrate the reliability and validity of both the theoretical model and the finite element analysis.

Establishment of Optimization Objective Function
To enhance the machining efficiency and tool life expectancy of Inconel 718, optimizing the processing parameters becomes essential.In this study, surface roughness, material removal rate, and their combination were considered the optimization objectives.The fmincon algorithm, a nonlinear minimization method employing the interior point approach, was utilized to determine the optimal processing parameters.Known for its high accuracy in solving nonlinear optimization problems, the fmincon algorithm is well suited for this task.The optimization procedure employing the fmincon algorithm is depicted in Figure 28.
Surface roughness is an important evaluation index in the machining process of Inconel 718; therefore, a model with surface roughness as the optimization objective was developed, and the optimization objective function can be expressed as:

Establishment of Optimization Objective Function
To enhance the machining efficiency and tool life expectancy of Inconel 718, optimizing the processing parameters becomes essential.In this study, surface roughness, material removal rate, and their combination were considered the optimization objectives.The fmincon algorithm, a nonlinear minimization method employing the interior point approach, was utilized to determine the optimal processing parameters.Known for its high accuracy in solving nonlinear optimization problems, the fmincon algorithm is well suited for this task.The optimization procedure employing the fmincon algorithm is depicted in Figure 28.
The parameter n in Equation ( 31) can be expressed as: To ensure the machining stability, the processing parameters of Inconel 718 need to be selected within the stability region of the stability lobes, as shown in Figures 16-19.Surface roughness is an important evaluation index in the machining process of Inconel 718; therefore, a model with surface roughness as the optimization objective was developed, and the optimization objective function can be expressed as: In Equation ( 24), the parameter C is influenced by both the geometry of the milling cutter and material of the workpiece.Parameters c 1 , c 2 , and c 3 are the influencing coefficients of the milling parameters.
Taking logarithms on both sides of Equation ( 24) yields: Taking y = ln R a , c 0 = ln C, x 1 = ln n, x 2 = ln f z , and x 3 = ln a p and substituting these parameters into Equation (25) yields: Based on this, multiple linear regression equations are built sequentially: Equation ( 27) is changed to a matrix representation: The parameters in Equation ( 28) can be expressed as: x 11 x 12 x 13 x 14 1 x 21 x 22 x 23 x 24 . . .
Using the experimental results in Reference [52], the optimization objective Function (24) can be changed to: In metal milling, improving productivity is critical, and the material removal rate is usually used to measure productivity.Therefore, a model with material removal rate as the optimization objective was developed: The parameter n in Equation ( 31) can be expressed as: To ensure the machining stability, the processing parameters of Inconel 718 need to be selected within the stability region of the stability lobes, as shown in Figures 16-19.

Optimization Results Analysis
When the minimum surface roughness is used as the optimization objective, it can be observed from Figure 29 that the solution of the optimization objective function gradually converges as the number of iterations increases.According to the optimal milling parameters recorded in Table 13, it can be concluded that the surface roughness reached the optimal value at the 19th generation, and the optimal surface roughness is 0.43 µm.

Optimization Results Analysis
When the minimum surface roughness is used as the optimization objective, it can be observed from Figure 29 that the solution of the optimization objective function gradually converges as the number of iterations increases.According to the optimal milling parameters recorded in Table 13, it can be concluded that the surface roughness reached the optimal value at the 19th generation, and the optimal surface roughness is 0.43 μm.When the maximum material removal rate was used as the optimization objective, it can be seen from Figure 30 that the solution of the optimization objective function gradually converged as the number of iterations increased.According to the optimal milling parameters recorded in Table 14, it can be concluded that the material removal rate reached the optimal value at the 11th generation, and the optimal material removal rate is 58,788.32(mm 3 /min).When the maximum material removal rate was used as the optimization objective, it can be seen from Figure 30 that the solution of the optimization objective function gradually converged as the number of iterations increased.According to the optimal milling parameters recorded in Table 14, it can be concluded that the material removal rate reached the optimal value at the 11th generation, and the optimal material removal rate is 58,788.32(mm 3 /min).If the maximum material removal rate and minimum surface roughness are taken as the research objective of the multi-objective optimization, the mathematical model of the multi-objective optimization is as follow: min f (x) = min(ω 1 f 1 (x) + ω 2 f 2 (x)) a pmin ≤ a p ≤ a pmax (33) In Equation (33), the parameters ω 1 and ω 2 are weighting values of the machining efficiency and surface roughness, respectively.
It can be seen from Figure 31 that the solution of the optimization objective function gradually converges as the number of iterations increases.According to the optimal milling parameters recorded in Table 15, it can be concluded that the optimal material removal rate reached 4199.2 mm 3 /min, while the optimal surface roughness reached 3.5 µm.
In Equation (33), the parameters 1 ω and 2 ω are weighting values of the machining efficiency and surface roughness, respectively.It can be seen from Figure 31 that the solution of the optimization objective function gradually converges as the number of iterations increases.According to the optimal milling parameters recorded in Table 15, it can be concluded that the optimal material removal rate reached 4199.2 mm 3 /min, while the optimal surface roughness reached 3.5 μm.

Conclusions
This study addresses a critical gap in the realm of metal cutting, specifically focusing on the machining of nickel-based superalloy Inconel 718.While previous research has made strides in modeling milling forces and stability analysis, a distinct gap persists in the understanding of milling mechanisms and chattering stability of Inconel 718.This work introduced a tailored milling force model for ball-end milling cutters, offering insights into machining forces for Inconel 718.By coupling finite element analysis with material failure criteria, the intricate interplay between processing parameters and material behavior during Inconel 718 machining is elucidated.A noteworthy contribution of this study lies in its meticulous examination of chattering stability during the machining process.The establishment of stability lobe diagrams through dynamic equations enhances the selection of stable cutting parameters, bolstering machining outcomes for challenging materials like Inconel 718.Furthermore, this study pioneered a multi-objective optimization approach, concurrently considering conflicting objectives such as surface roughness and material removal rate.By employing this algorithm, optimal machining parameters were derived, signifying a path towards improved machining efficiency and quality for Inconel 718.
The main conclusions are as follow: 1.
The finite element analysis revealed notable temperature increases in the machining zone with rising spindle speed, radial and axial depths of cut, and feed per tooth.Tool stress, influenced mainly by feed per tooth and radial and axial depths of cut, exhibited significant changes.While spindle speed moderately affected milling force for each axis, the feed per tooth and radial depth of cut substantially increased forces in the x-, y-, and x-axes.2.
Chattering stability was investigated, with tool radius enhancement aiding stability, while radial depth of cut undermined it.An improved damping ratio enhances stability, while the system's inherent ability to resist vibrations or chatter improves with the increase in natural frequency.

3.
Validation against experiments and calculated results from the milling force model and finite element analysis exhibited minor deviations.

Figure 1 .
Figure 1.Milling force model of ball-end milling cutter.

Figure 2 .
Figure 2. Finite element model of tool and workpiece.

Figure 2 .
Figure 2. Finite element model of tool and workpiece.

Figure 3 .
Figure 3. Influence of processing parameters on temperature.(a) Influence of spindle speed on temperature; (b) Influence of feed per tooth on temperature; (c) Influence of radial depth of cut on temperature; (d) Influence of axial depth of cut on temperature.

Figure 3 .
Figure 3. Influence of processing parameters on temperature.(a) Influence of spindle speed on temperature; (b) Influence of feed per tooth on temperature; (c) Influence of radial depth of cut on temperature; (d) Influence of axial depth of cut on temperature.

Figure 4 .
Figure 4. Influence of machining parameters on stress.(a) Influence of spindle speed on stress; (b) Influence of feed per tooth on stress; (c) Influence of radial depth of cut on stress; (d) Influence of axial depth of cut on stress.

Figure 5 .
Figure 5.Comparison of milling forces of each axis at different spindle speeds.(a) Influence of spindle speed on milling force of x-axis; (b) Influence of spindle speed on milling force of y-axis; (c) Influence of spindle speed on milling force of z-axis.

Figure 5 .
Figure 5.Comparison of milling forces of each axis at different spindle speeds.(a) Influence of spindle speed on milling force of x-axis; (b) Influence of spindle speed on milling force of y-axis; (c) Influence of spindle speed on milling force of z-axis.

Figure 6 .
Figure 6.Comparison of milling forces of each axis at different feeds per tooth.(a) Influence of feed per tooth on milling force of x-axis; (b) Influence of feed per tooth on milling force of y-axis; (c) Influence of feed per tooth on milling force of z-axis.

Figure 6
Figure6demonstrates a clear and consistent increasing trend in the milling forces for each axis with the increase in the feed per tooth.This observation can be attributed to the larger volume of workpiece material being removed by the cutter as the feed per tooth

Figure 6 .
Figure 6.Comparison of milling forces of each axis at different feeds per tooth.(a) Influence of feed per tooth on milling force of x-axis; (b) Influence of feed per tooth on milling force of y-axis; (c) Influence of feed per tooth on milling force of z-axis.

Figure 7 .
Figure 7.Comparison of milling forces of various axes under different radial depths of cut.(a) Influence of radial depth of cut on milling force of x-axis; (b) Influence of radial depth of cut on milling force of y-axis; (c) Influence of radial depth of cut on milling force of z-axis.

Figure 7 .
Figure 7.Comparison of milling forces of various axes under different radial depths of cut.(a) Influence of radial depth of cut on milling force of x-axis; (b) Influence of radial depth of cut on milling force of y-axis; (c) Influence of radial depth of cut on milling force of z-axis.

Figure 8 .
Figure 8.Comparison of the milling forces of different axes under different axial depths of cut.(a) Influence of axial depth of cut on milling force of x-axis; (b) Influence of axial depth of cut on milling force of y-axis; (c) Influence of axial depth of cut on milling force of z-axis.

Figure 8 .
Figure 8.Comparison of the milling forces of different axes under different axial depths of cut.(a) Influence of axial depth of cut on milling force of x-axis; (b) Influence of axial depth of cut on milling force of y-axis; (c) Influence of axial depth of cut on milling force of z-axis.

Figure 9 .
Figure 9. Dynamic model of the machine-tool system.

Figure 9 .
Figure 9. Dynamic model of the machine-tool system.

Figure 10 .
Figure 10.Dynamic change in the cutting thickness.

Figure 11 .
Figure 11.Schematic diagram of the equipment connection of the modal experiment.

Figure 12 .
Figure 12.Equipment connection for the modal experiment.

Figure 11 .
Figure 11.Schematic diagram of the equipment connection of the modal experiment.

Figure 11 .
Figure 11.Schematic diagram of the equipment connection of the modal experiment.

Figure 12 .
Figure 12.Equipment connection for the modal experiment.

Figure 12 .
Figure 12.Equipment connection for the modal experiment.

Figure 13 .
Figure 13.Distribution of measuring points on tool.

Figures 14 and 15
Figures 14 and 15 show the tool modal of the x-and y-directions, respectively.

Figure 13 .
Figure 13.Distribution of measuring points on tool.

Figure 14 .
Figure 14.Tool modal of x-direction.(a) First order of x-direction; (b) Second order of x-direction; (c) Third order of x-direction.

Figure 14 . 39 Figure 15 .
Figure 14.Tool modal of x-direction.(a) First order of x-direction; (b) Second order of x-direction; (c) Third order of x-direction.Materials 2023, 16, x FOR PEER REVIEW 20 of 39

Figure 15 .
Figure 15.Tool modal of the y-direction.(a) First order of y-direction; (b) Second order of y-direction; (c) Third order of y-direction.

Figure 15 .
Figure 15.Tool modal of the y-direction.(a) First order of y-direction; (b) Second order of y-direction; (c) Third order of y-direction.
radial depth of cut, natural frequency n ω , and damping ratio ζ on the cutting stability.The stability lobes are plotted for different tool radii, radial depths of cut, natural frequencies n ω , and damping ratios ζ , as shown in Figures 16-19.

Figure 17 .
Figure 17.Influence of radial depth of cut on stability.Figure 17.Influence of radial depth of cut on stability.

Figure 17 .
Figure 17.Influence of radial depth of cut on stability.Figure 17.Influence of radial depth of cut on stability.

Figure 17 .
Figure 17.Influence of radial depth of cut on stability.

Figure 18 .
Figure 18.Influence of natural frequency on stability.Figure 18. Influence of natural frequency on stability.

Figure 18 . 39 Figure 19 .
Figure 18.Influence of natural frequency on stability.Figure 18. Influence of natural frequency on stability.Materials 2023, 16, x FOR PEER REVIEW 22 of 39

Figure 19 .
Figure 19.Influence of damping ratio on stability.

Figure 20 .
Figure 20.Schematic diagram of the connection of the experimental equipment.Figure 20.Schematic diagram of the connection of the experimental equipment.

Figure 20 .
Figure 20.Schematic diagram of the connection of the experimental equipment.Figure 20.Schematic diagram of the connection of the experimental equipment.

Figure 20 .
Figure 20.Schematic diagram of the connection of the experimental equipment.

Materials 2023 , 39 Figure 22 .
Figure 22.Comparison of the theoretical and experimental milling forces.

Figure 23 .
Figure 23.Comparison of the experimental and theoretical milling forces: (a) experimental milling force; (b) theoretical milling force.

Figure 22 .
Figure 22.Comparison of the theoretical and experimental milling forces.

Materials 2023 , 39 Figure 22 .Figure 23 .
Figure 22.Comparison of the theoretical and experimental milling forces.The milling forces measured through milling experiment are presented in Figure23a, and the milling forces calculated using the theoretical milling force model are presented in Figure23b.The milling forces of each axis obtained from the theoretical model align closely with the experimental results, validating the accuracy and correctness of the proposed model.

Figure 23 .
Figure 23.Comparison of the experimental and theoretical milling forces: (a) experimental milling force; (b) theoretical milling force.

Figure 24 .
Figure 24.Influence of the spindle speed on the average milling force of each axis.

Figure 25 .
Figure 25.Influence of the axial depth of cut on the average milling force of each axis.

Figure 25
Figure 25 reveals that the average milling force for the x-, y-, and z-axes all demonstrated an increasing trend with the augmentation of the axial depth of cut.The z-axis

Figure 24 .
Figure 24.Influence of the spindle speed on the average milling force of each axis.

Materials 2023 , 39 Figure 24 .
Figure 24.Influence of the spindle speed on the average milling force of each axis.

Figure 25 .
Figure 25.Influence of the axial depth of cut on the average milling force of each axis.

Figure 25
Figure 25 reveals that the average milling force for the x-, y-, and z-axes all demonstrated an increasing trend with the augmentation of the axial depth of cut.The z-axis

Figure 25 .
Figure 25.Influence of the axial depth of cut on the average milling force of each axis.

Figure 26 .
Figure 26.Influence of the feed per tooth on the average milling force of each axis.

Figure 26 .
Figure 26.Influence of the feed per tooth on the average milling force of each axis.

Figure 28 .
Figure 28.Flowchart of the procedure of optimization using the fmincon algorithm.Figure 28.Flowchart of the procedure of optimization using the fmincon algorithm.

Figure 28 .
Figure 28.Flowchart of the procedure of optimization using the fmincon algorithm.Figure 28.Flowchart of the procedure of optimization using the fmincon algorithm.

Figure 29 .
Figure 29.Iteration curve of the surface roughness.

Figure 30 .
Figure 30.Iterative curve of material removal rate.

Table 1 .
Data from the recognition experiments of the milling force coefficients.

Table 2 .
Coefficients of the milling force.

Shearing Force Coefficient (N/mm 2 ) Value Edge Force Coefficient (N/mm 2 )
ae K

Table 3 .
Processing parameters used in the finite element analysis.

Table 5 .
Thermal softening parameters of Inconel 718. c

Strain Rate Sensitivity Factor m 1 High Strain Rate Sensitivity Factor m 2 Reference Plastic Strain Rate . ε 0 Strain Rate Critical Value .
1 × 10 7 s −1

Table 7 .
Results of the modal experiments.

Table 7 .
Results of the modal experiments.

Table 7 .
Results of the modal experiments.
Figures 14 and 15 show the tool modal of the x-and y-directions, respectively.

Table 8 .
Parameters of the single factor experiments.

Table 9 .
Parameters of the orthogonal experiments.

Table 10 .
Result of the orthogonal experiments.

Table 10 .
Result of the orthogonal experiments.

Table 11 .
Comparison between the theoretical and experimental results of the milling force.

Table 12 .
Comparison between the finite element analysis and experimental results of the milling force.

Table 13 .
Results of the optimization of the surface roughness.

Table 14 .
Results of the optimization of the material removal rate.

Table 15 .
Results of the multi-objective optimization.

Table 15 .
Results of the multi-objective optimization.