Stability Research Considering Non-Linear Change in the Machining of Titanium Thin-Walled Parts

Aiming to solve the problem whereby the damping process effect is significant and difficult to measure during low-speed machining of titanium alloy thin-walled parts, the ploughing coefficient of the flank face is obtained based on the frequency-domain decomposition (FDD) of the measured vibration signal and the energy balance principle, and then the process-damping prediction model is obtained. Aiming to solve the problem of non-linear change of dynamic characteristics of a workpiece caused by the material removal effect in the machining of titanium alloy thin-walled parts, a prediction model of dynamic characteristics of a workpiece is established based on the structural dynamic modification method. Meanwhile, the effect of material removal on the process-damping coefficient is studied, and the internal relationship between the process-damping coefficient and the dynamic characteristics of the workpiece is revealed. The stability lobe diagram is obtained by the full discretization in the titanium alloy milling process. The correctness of the model and stability prediction is verified by experiments under different working conditions. It is found that the coupling characteristics of process-damping and workpiece dynamic characteristics control the stability of the milling process. The research results can provide theoretical support for accurate characterization and process optimization of titanium alloy thin-walled workpiece milling.


Introduction
Titanium alloy thin-walled parts have found an increasingly wide use in the aviation manufacture industry because they can meet the requirements of high performance and high stability. Due to the low stiffness and high material removal rate, chattering occurs very easily in the machining process, which reduces the surface precision of the workpiece and the service life of cutters and machine tools, limiting production efficiency. It has become an urgent problem to be solved in the aviation manufacturing industry.
Titanium alloys are mostly machined at low speed. In the process of low-speed machining, the periodic friction and extrusion between the tool flank and the machined surface will produce a ploughing effect, which improves the processing stability [1,2]. Sisson et al [3] and Peters et al [4] studied the process-damping phenomenon in a low-speed cutting process with an experimental method. It was found that the ultimate cutting depth obtained by considering process-damping has been effectively improved, especially in the low-speed range. Wu et al [5] constructed a ploughing force model which can accurately characterize the process-damping. It was assumed that ploughing force was a linear function of the total volume of the cutter-worker interference and ploughing coefficient. Subsequently, researchers have conducted in-depth research on the accurate acquisition of the total characteristics of the workpiece and the tool. The shear force is produced on the rake face due to shear effect, and the ploughing force is produced on the flank face due to friction and extrusion effect in the process of cutting material. The shear force and ploughing force can be broken into two sections-the static milling force and the ploughing milling force. Ignoring the modal coupling effect of the process system, the dynamic equation of cutting process is established, as shown in Equation (1). characteristics of the workpiece and the tool. The shear force is produced on the rake face due to shear effect, and the ploughing force is produced on the flank face due to friction and extrusion effect in the process of cutting material. The shear force and ploughing force can be broken into two sectionsthe static milling force and the ploughing milling force. Ignoring the modal coupling effect of the process system, the dynamic equation of cutting process is established, as shown in Equation 1.   Static cutting force has no influence on the regeneration effect of machining process, so it is often ignored in the process of dynamic modeling.
The dynamic equation of processing is obtained as shown in Equation 3. ..
are the mass, structural damping, and stiffness matrices of the cutting system, respectively. [m yy ], [c yy ], [k yy ] are the mass, structural damping, and stiffness matrices of the workpiece, respectively. F i (t) is the cutting force in the X and Y directions. i(t) is the response caused by the cutting force. i = x or y.
The calculation formula of the cutting force in Y direction is shown in Equation (2).
F ys,st (t) and F yp,st (t) are the static cutting forces produced by shearing and friction on rank and flank face. F ys,dy (t) and F yp,dy (t) are the dynamic cutting force produced by shearing and friction on rank and flank face.
Static cutting force has no influence on the regeneration effect of machining process, so it is often ignored in the process of dynamic modeling.
The dynamic equation of processing is obtained as shown in Equation (3).
g(φ i,j (t)) is a window function for judging the processing state. Its expression is shown in Equation (6). K Ts and K Rs are the shear force coefficients on the rake face. K T f and K R f are the ploughing force coefficients on the flank face.
l w is the boundary length between static indented and dynamic indented volumes.
The angle of cutting rotation for the cutter tooth is φ i,j (t). The expression of the tool rotation angle considering the tool helix angle, is shown in Equation (8).
where Ω is the spindle speed, N is the number of cutter teeth, β is the helix angle of the cutting tool, dz is the length of the axial unit. Similarly, the dynamic shear force F xs,dy (t) and the dynamic ploughing force F xp,dy (t) in the X direction are obtained. We can get the dynamic equations with equivalent process-damping as follows.
x . y = 1 2 a p a xx a xy a yx a yy ∆x ∆y where c px and c py are the process-damping coefficients in machining process.

Prediction of Part Dynamics
The change of dynamic characteristics caused by material removal effect can be regarded as structural modification. Combining the original frequency-response function of the workpiece with the change value of the frequency-response function caused by material removal, the inherent properties of the workpiece after material removal can be obtained. The response of process system composed of workpiece and fixture under the cutting force in frequency domain can be expressed as: where τ is the unit imaginary number. The receptance matrix, [α] is defined by: The change of inherent properties of workpiece caused by material removal effect can be expressed as [∆M], [∆C] and [∆K]. The dynamic structure-modification matrix is shown in Equation (14).
The receptance matrix is obtained by simultaneous Equation (13) and Equation (14). [ The modified receptance matrix of workpiece structure can be expressed as follows.

Solving the Process-Damping Coefficient
Side milling is a discontinuous cutting style. Because of the inconvenience of the experimental acquisition equipment and the periodic change of the cutting signal, the process-damping coefficient cannot be obtained directly in milling experiments. In this paper, four eddy current displacement sensors are used to obtain the vibration signal of the workpiece and the cutting system during the machining process, as shown in Figure 2. Then the frequency-domain decomposition method is used to analyze the vibration signal of the process system, and the total damping coefficient of the process system is obtained. Finally, the ploughing coefficients of the process system are obtained.
The workpiece deformation y(t) in the Y direction can be expressed by the mode shape U and the mode displacement q(t).
Converting the power spectrum of the time domain deformation of experimental points y1, y2 into the frequency domain, the formula is as follows. Y( jω) is the Fourier spectrum of workpiece deformation, * and T represents the conjugate of a complex number and transpose of a matrix.
Q( jω) is the Fourier spectrum of q(t). Q( jω)Q * T ( jω) is the 2*2 matrix. Due to the orthogonality of the modal, the asymmetry element is 0. In addition, because the modal of the system can be well separated (as shown in Figure 3), the Equation (19) can be simplified as follows.
The magnitude of direct frequency-response functions (FRF) at any point can be obtained by the ratio of the power spectrum of the dynamic displacement of the point to the power spectrum of the exciting force. Therefore, the frequency-response function of position 1 is calculated as shown in Equation (21).
where S f f is a fixed value [23]. The power spectrum of the deflection at point y1 can be transformed into the following form.
When the frequency range is the same, the frequency-response function can be expressed by modal parameters.
The extreme points occur at t = iπ/ω d . The absolute value of F −1 can be expressed as a linear function. Substituting t = iπ/ω d in Equation (25), we get the peaks and valleys of F −1 .
Logarithmic operation is conducted to linearize Equation (26).
The total damping ratio can be calculated by Equation (2).
It is impractical to obtain the frequency-response function in machining process, so the autocorrelation of q1 is resolved by the vibration signal measured in cutting process.
The power spectrum matrix of the workpiece vibration measured at points y1, y2 can be expressed as follows.
The similarity between eigenvectors and modal shapes is determined on each frequency line using modal assurance criteria (MAC).
When the value of MAC(ω i ) is close to 1, Q 1 Q * 1 ≈ κ 1 . The difference between the total damping coefficient and the structural damping coefficient is used to calculate the process-damping coefficient.
In the same way, the X-directional process-damping coefficient can be solved by Equation (33).
The energy of the average process-damping effect is equal to the energy consumed by the dynamic plough force in the rotational period of the spindle, thus the ploughing force coefficients can be obtained.

Experimental Setup
The experiment was carried out on a three-axis machining center manufactured by Dalian machine tool(Dalian machine tool group, Dalian, Liaoning, China). Its model is VDL-1000E. Solid carbide cutting tools with four teeth, diameter 10 and helix angle 30 • are used for processing thin-walled titanium alloy parts. Its model is GM-4E-D10.0. The coating material is TiAlN. The extended length of the cutting tool is 110 mm. The geometric dimension of thin-walled titanium alloy parts is 200 × 200 × 5 mm. The workpiece is fixed on the vise. The vice is bolted to the workbench. The extended length of the workpiece is 100 mm. The milling mode is down milling and dry cutting. The rotary dynamometer produced by Kistler Company (Winterthur, Switzerland) was used to collect cutting force data. Its model is 5236B. Three-component acceleration sensors produced by PCB Sensor Company was used to collect the acceleration data of cutting process. Its model is 368F. The three-dimensional sensitivity is 2.462 mv/g, 2.534 mv/g and 2.487 mv/g, respectively. The initial modal parameters of the cutting system and workpiece are obtained by using the unidirectional acceleration sensor produced by PCB Sensor Company and the impact hammer. The vibration signals of the workpiece and the tool are measured using four ST-2-U-05-00-20-KH07 eddy current sensors. The experimental machining site is shown in Figure 2.

Experimental Setup
The experiment was carried out on a three-axis machining center manufactured by Dalian machine tool(Dalian machine tool group, Dalian, Liaoning, China). Its model is VDL-1000E. Solid carbide cutting tools with four teeth, diameter 10 and helix angle 30° are used for processing thinwalled titanium alloy parts. Its model is GM-4E-D10.0. The coating material is TiAlN. The extended length of the cutting tool is 110 mm. The geometric dimension of thin-walled titanium alloy parts is 200 × 200 × 5 mm. The workpiece is fixed on the vise. The vice is bolted to the workbench. The extended length of the workpiece is 100 mm. The milling mode is down milling and dry cutting. The rotary dynamometer produced by Kistler Company (Winterthur, Switzerland) was used to collect cutting force data. Its model is 5236B. Three-component acceleration sensors produced by PCB Sensor Company was used to collect the acceleration data of cutting process. Its model is 368F. The threedimensional sensitivity is 2.462 mv/g, 2.534 mv/g and 2.487 mv/g, respectively. The initial modal parameters of the cutting system and workpiece are obtained by using the unidirectional acceleration sensor produced by PCB Sensor Company and the impact hammer. The vibration signals of the workpiece and the tool are measured using four ST-2-U-05-00-20-KH07 eddy current sensors. The experimental machining site is shown in Figure 2.

The Acquisition of Dynamic Characteristics
Modal parameters of cutting system are usually obtained by hammer impact experiment. Generally, the real part and imaginary part of the frequency-response function are fitted by curve fitting method to obtain the modal parameters of the process system, as shown in Figure 3. Machining

The Acquisition of Dynamic Characteristics
Modal parameters of cutting system are usually obtained by hammer impact experiment. Generally, the real part and imaginary part of the frequency-response function are fitted by curve fitting method to obtain the modal parameters of the process system, as shown in Figure 3. Machining process is mainly controlled by first-order modal parameters. The initial first-order modal parameters of the unmachined workpiece and the tool are shown in Table 1. process is mainly controlled by first-order modal parameters. The initial first-order modal parameters of the unmachined workpiece and the tool are shown in Table 1.  The thickness of thin-walled titanium alloy parts is reduced from 5 mm to 0 mm by machining along the wall thickness. The process diagram is shown in Figure 4. The comparison results between the experimental and predicted values of the modal parameters of the workpiece at different machining positions is shown in Table 2.  Table 2 showed that with the removal of the material, the natural frequency of the workpiece increased, and the stiffness and damping ratio decreased gradually. The vibration resistance of workpiece is weakened. In addition, the changes of them displayed a non-linear feature. Therefore,  The thickness of thin-walled titanium alloy parts is reduced from 5 mm to 0 mm by machining along the wall thickness. The process diagram is shown in Figure 4. The comparison results between the experimental and predicted values of the modal parameters of the workpiece at different machining positions is shown in Table 2.  The thickness of thin-walled titanium alloy parts is reduced from 5 mm to 0 mm by machining along the wall thickness. The process diagram is shown in Figure 4. The comparison results between the experimental and predicted values of the modal parameters of the workpiece at different machining positions is shown in Table 2.  Table 2 showed that with the removal of the material, the natural frequency of the workpiece increased, and the stiffness and damping ratio decreased gradually. The vibration resistance of workpiece is weakened. In addition, the changes of them displayed a non-linear feature. Therefore,   Table 2 showed that with the removal of the material, the natural frequency of the workpiece increased, and the stiffness and damping ratio decreased gradually. The vibration resistance of workpiece is weakened. In addition, the changes of them displayed a non-linear feature. Therefore, if the damping ratio is assumed to be a fixed value or an equal proportion change, the prediction of dynamic force and stability will deviate.
By comparing the predicted value and measured value of the modal parameters of the workpiece, it can be seen that the maximum error of natural frequency, stiffness and damping ratio is 10.07%, 8.85% and 7.90%, respectively. The error is within acceptable range to verify the accuracy of the workpiece dynamic model.

Ploughing Force Coefficient
With ap = 6 mm, ae = 0.5 mm, n = 900 r/min, f = 0.1 mm/z, the steps to obtain the total damping ratio in the feed direction and normal direction of position A11 are shown in Figure 5.
By comparing the predicted value and measured value of the modal parameters of the workpiece, it can be seen that the maximum error of natural frequency, stiffness and damping ratio is 10.07%, 8.85% and 7.90%, respectively. The error is within acceptable range to verify the accuracy of the workpiece dynamic model.

Ploughing Force Coefficient
With ap = 6 mm, ae = 0.5 mm, n = 900 r/min, f = 0.1mm/z，the steps to obtain the total damping ratio in the feed direction and normal direction of position A11 are shown in Figure 5.   By analyzing the cutting vibration obtained by the cutting process shown in Figure 4, the variation law of the vibration amplitude and the process-damping coefficient with the material removal effect is obtained, as shown in Figure 6.  The overall damping ratios are estimated by substituting the slope of the fitted lines in Equation (27), ξ r,x = 0.0801, ξ r,y = 0.0798.The process-damping coefficient is evaluated as c p,x = 2m x ω n,x (ξ r,x − ξ s,x ) = 53.39(Ns/m), c p,y = 2m y ω n,y (ξ r,y − ξ s,y ) = 31.09(Ns/m). The ploughing force coefficients are obtained to be K T f = 3.735 × 10 13 N/m 2 , K R f = 1.208 × 10 13 N/m 2 .
By analyzing the cutting vibration obtained by the cutting process shown in Figure 4, the variation law of the vibration amplitude and the process-damping coefficient with the material removal effect is obtained, as shown in Figure 6.  By analyzing the cutting vibration obtained by the cutting process shown in Figure 4, the variation law of the vibration amplitude and the process-damping coefficient with the material removal effect is obtained, as shown in Figure 6.   Figure 6a shows that under the same cutting parameters, the vibration amplitude of the machining system in the Y directions increases with the removal of workpiece materials. The change of vibration in the X direction is just the opposite to that in the Y direction. Meanwhile, it is found that the change rate of vibration amplitude in Y direction is greater than that in X direction.
It can be seen from Figure 6b that the process-damping coefficient in the Y directions increases with the removal of workpiece materials. The process-damping coefficient considering the dynamic characteristics of the workpiece (A) is smaller than that without considering the dynamic characteristics of the workpiece (B). The process-damping coefficient in the X directions decreases with the removal of workpiece materials.

Milling Stability
The frequency-domain method [24], discrete method and numerical method [25] are often used to solve the dynamic equation in machining process, and then then the stability of the cutting process is obtained. The frequency-domain method has the highest calculation efficiency, but the prediction accuracy is the lowest, especially it is not applicable to small radial depth of cut. The numerical method has the highest prediction accuracy, but its computational efficiency is the lowest because many equations need to be solved directly. The discrete method can be divided into the semi-discretization method [26], the full-discretization method [27] and the time FEM [28]. The semi-discretization method and the full-discretization method are most widely used. The matrix exponential function involved in the full-discretization method is only dependent on the spindle speed and independent of the axial depth of cut, so it is more efficient than the semi-discretization method. Therefore, we use the full discrete method to obtain the stability of the machining process.
The research results in Section 5.3 show that the process-damping coefficient has a non-linear relationship with the dynamic characteristics of the process system. Taking A15 machining position as an example, the influence of coupling characteristics between process-damping and dynamic characteristics of process system on milling stability is investigated.
The stability lobe diagrams with the coupled and the uncoupled process-damping, and the dynamic characteristics of the process system are shown in Figure 7a. A series of experiments (ae = 0.5 mm, f = 0.1 mm/tooth) were carried out to verify the accuracy of the stable lobe diagram. The remaining cutting parameters are detailed in Table 3. Two points A and B are selected for experimental verification to verify the correctness of the predicted results. Figure 7b,c are acceleration signals in time domain. Figure 7d,e are Fourier transform results of acceleration signals.  The maximum amplitude of point A acceleration signal is 18 m/s 2 , and its amplitude variation has better convergence, as shown in Figure 7b. Fourier transform is applied to the acceleration signal in time domain, and the result is shown in Figure 7d. It is found that the spectrum energy mainly concentrates on the cutter teeth passing frequency and its higher harmonics. The maximum amplitude of point B acceleration signal is 48 m/s 2 , and its amplitude change law is from small to large, which does not have convergence, as shown in Figure 7c. Fourier transform is applied to the acceleration signal in time domain, and the result is shown in Figure 7e. It is found that the spectrum energy is mainly concentrated near the first natural frequency of the workpiece.
Meanwhile, the white light interferometer is used to measure the micromorphology of machined surface with A and B parameters and the results are shown in Figure 8a   The maximum amplitude of point A acceleration signal is 18 m/s 2 , and its amplitude variation has better convergence, as shown in Figure 7b. Fourier transform is applied to the acceleration signal in time domain, and the result is shown in Figure 7d. It is found that the spectrum energy mainly concentrates on the cutter teeth passing frequency and its higher harmonics. The maximum amplitude of point B acceleration signal is 48 m/s 2 , and its amplitude change law is from small to large, which does not have convergence, as shown in Figure 7c. Fourier transform is applied to the acceleration signal in time domain, and the result is shown in Figure 7e. It is found that the spectrum energy is mainly concentrated near the first natural frequency of the workpiece.
Meanwhile, the white light interferometer is used to measure the micromorphology of machined surface with A and B parameters and the results are shown in Figure 8a,b. Two-dimensional Fourier transform (2DFFT) is applied to the image of micro-surface topography, and the results are shown in Figure 8c,d. The maximum amplitude of point A acceleration signal is 18 m/s 2 , and its amplitude variation has better convergence, as shown in Figure 7b. Fourier transform is applied to the acceleration signal in time domain, and the result is shown in Figure 7d. It is found that the spectrum energy mainly concentrates on the cutter teeth passing frequency and its higher harmonics. The maximum amplitude of point B acceleration signal is 48 m/s 2 , and its amplitude change law is from small to large, which does not have convergence, as shown in Figure 7c. Fourier transform is applied to the acceleration signal in time domain, and the result is shown in Figure 7e. It is found that the spectrum energy is mainly concentrated near the first natural frequency of the workpiece.
Meanwhile, the white light interferometer is used to measure the micromorphology of machined surface with A and B parameters and the results are shown in Figure 8a   The machined surface obtained at point A has lower level of surface roughness (Sa = 0.85 um, Ra = 0.47 um) and surface waviness, as shown in Figure 8a. The only significant spectral property of the machined surface at point A is related to the marking belonging to feed-rate, as shown in Figure  8c. The machined surface obtained at point A has the higher level of surface roughness (Sa = 1.63um, Ra = 1.03 um) and surface waviness, as shown in Figure 8b. The spectral characteristics of the surface topography generated by point B are inclined modes related to the chatter frequency, as shown in Figure 8d. In conclusion, A is the cutting stability point and B is the cutting chatter point. It is found that the stability obtained by considering the coupling characteristics of dynamic characteristics of workpiece and process-damping has higher prediction accuracy.
Milling stability of the coupled process-damping and dynamic characteristics of process system is lower than that of the uncoupled, as shown in Figure 7a. The reason is that the process-damping coefficient obtained from the coupling process-damping and the dynamic characteristics of the process system is small. Meanwhile, the difference between them decreases with the increase of spindle speed. The reason is that the damping coefficient descends with the increase of spindle speed.
The stability lobe diagrams with and without process-damping are shown in Figure 9. A series of experiments (ae = 0.5 mm, f = 0.1 mm/tooth) were carried out to verify the accuracy of the stable lobe diagram. The remaining cutting parameters are detailed in Table 4. It is found that the milling stability considering process-damping is higher than that without process-damping in titanium alloy milling process, especially low spindle speed. The milling stability decreases with the increase of spindle speed. The reason is that the process-damping coefficient is inversely proportional to the spindle speed.  The machined surface obtained at point A has lower level of surface roughness (Sa = 0.85 um, Ra = 0.47 um) and surface waviness, as shown in Figure 8a. The only significant spectral property of the machined surface at point A is related to the marking belonging to feed-rate, as shown in Figure 8c. The machined surface obtained at point A has the higher level of surface roughness (Sa = 1.63 um, Ra = 1.03 um) and surface waviness, as shown in Figure 8b. The spectral characteristics of the surface topography generated by point B are inclined modes related to the chatter frequency, as shown in Figure 8d. In conclusion, A is the cutting stability point and B is the cutting chatter point. It is found that the stability obtained by considering the coupling characteristics of dynamic characteristics of workpiece and process-damping has higher prediction accuracy.
Milling stability of the coupled process-damping and dynamic characteristics of process system is lower than that of the uncoupled, as shown in Figure 7a. The reason is that the process-damping coefficient obtained from the coupling process-damping and the dynamic characteristics of the process system is small. Meanwhile, the difference between them decreases with the increase of spindle speed. The reason is that the damping coefficient descends with the increase of spindle speed.
The stability lobe diagrams with and without process-damping are shown in Figure 9. A series of experiments (ae = 0.5 mm, f = 0.1 mm/tooth) were carried out to verify the accuracy of the stable lobe diagram. The remaining cutting parameters are detailed in Table 4. It is found that the milling stability considering process-damping is higher than that without process-damping in titanium alloy milling process, especially low spindle speed. The milling stability decreases with the increase of spindle speed. The reason is that the process-damping coefficient is inversely proportional to the spindle speed.  The 3D stability lobe diagram in milling titanium alloy thin-walled parts is shown in Figure 10. It is found that the ultimate axial cutting depth decreases with material removal. The reason is that with the removal of the material, the natural frequency of the workpiece increased, and the stiffness and damping ratio decreased gradually. The vibration resistance of workpiece is weakened.  The 3D stability lobe diagram in milling titanium alloy thin-walled parts is shown in Figure 10. It is found that the ultimate axial cutting depth decreases with material removal. The reason is that with the removal of the material, the natural frequency of the workpiece increased, and the stiffness and damping ratio decreased gradually. The vibration resistance of workpiece is weakened.  The 3D stability lobe diagram in milling titanium alloy thin-walled parts is shown in Figure 10. It is found that the ultimate axial cutting depth decreases with material removal. The reason is that with the removal of the material, the natural frequency of the workpiece increased, and the stiffness and damping ratio decreased gradually. The vibration resistance of workpiece is weakened. Figure 10. Three-dimensional stable lobe diagram.

Conclusions
In this paper, two prominent problems in the processing of titanium alloy thin-walled partsthe non-linear changes of dynamic characteristics of the parts caused by material removal effect and the damping process effect is difficult to measure-are investigated in depth. Based on the frequencydomain decomposition of the measured vibration signal and the principle of energy balance, a process-damping prediction model is obtained. Based on the structural dynamic modification method, a prediction model of workpiece dynamic characteristics is established. The effect of material removal on process-damping coefficient is studied. The full discrete method is used to solve the

Conclusions
In this paper, two prominent problems in the processing of titanium alloy thin-walled parts-the non-linear changes of dynamic characteristics of the parts caused by material removal effect and the damping process effect is difficult to measure-are investigated in depth. Based on the frequency-domain decomposition of the measured vibration signal and the principle of energy balance, a process-damping prediction model is obtained. Based on the structural dynamic modification method, a prediction model of workpiece dynamic characteristics is established. The effect of material removal on process-damping coefficient is studied. The full discrete method is used to solve the stability of the milling process. The correctness of the model and stability prediction is verified by experiments with different working conditions. The research results can provide theoretical support for accurate characterization and process optimization of titanium alloy thin-walled workpiece milling. The specific conclusions are as follows: (1) Under the same cutting parameters, the vibration amplitude of the workpiece increases with the material removal, which leads to the corresponding increase of the process-damping coefficient. The variation of vibration amplitude and damping coefficient of cutting system is just the opposite.
(2) The process-damping coefficient obtained by considering the dynamic characteristics of the workpiece is less than that without considering the dynamic characteristics of the workpiece, and the gap between them increases with the removal of materials.
(3) Milling stability of the coupled process-damping and dynamic characteristics of process system is lower than that of the uncoupled. The reason is that the process-damping coefficient is relatively small. Meanwhile, the difference between them decreases with the increase of spindle speed. The reason is that the damping coefficient decreases with the increase of spindle speed.
(4) The milling stability is gradually reduced with the material is removed and the spindle speed is increased. Meanwhile, it is found that the effect of material removal on milling stability in the low-speed region is less than that in high speed region. The reason is that with the removal of materials, the vibration amplitude of the workpiece increases so that the damping coefficient of the process increases in the low-speed region, which further weakens the effect of material removal.
Author Contributions: The idea of this project was conceived by H.G. H.G. consulted the related papers on the milling stability of thin-walled parts, then sorted out the existing problems in this field, and finally wrote this paper. X.L. reviewed this project and proposed constructive guidance to make the article more complete.
Funding: This project is supported by the Projects of International Cooperation and Exchanges NSFC (51720105009).

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