Measurement-Based Modal Analysis and Stability Prediction on Turn-Milling of Hollow Turbine Blade

Hollow blades with honeycomb structures are increasingly used in the turbine engines for reducing weight and saving costs. +e hollow blade is a typical thin-walled structural part with low stiffness, the machining system of which is often unstable and likely to chatter. +e most effective solution to avoid the machining chatter is to guide the hollow blade to be machined in a stable machining zone. +is paper proposes a measurement-based approach for modal analysis and stability prediction of turn-milling hollow blade. +e impact test was carried out to achieve the FRF curves on the hollow blade and the milling tool. An extremum method was employed to obtain an equivalent FRF curve, from which the modal parameters involving the natural frequency, damping ratio, and stiffness were computed. Afterwards, the semidiscretization method was used to draw a stability lobe diagram to predict the stability when turn-milling hollow blades. +e experimental results confirm the feasibility of the predicted stability lobe diagram.


Introduction
Hollow blades with honeycomb structures (HCSs) are increasingly employed in the advanced turbine engines for lightweight and cost saving [1][2][3]. Superplastic forming (SPF)/diffusion bonding (DB) combined with numerical control machining (NCM) is the major solution to manufacture the hollow blades nowadays [4]. A near-net shape of the blade is formed by the processes of SPF and DB and then machined to its final shape by the NCM process. e hollow blade with HCSs is a typical thin-walled part, which has a curve surface and a weak stiffness. e setup of machining hollow blade with HCSs is complicated and multiaxis machining is often required (cf. Figure 1). e machining system of the hollow blade with HCSs is unstable and likely to vibrate strongly due to its weak stiffness [5]. Besides, the hollow blade with HCSs is often made of difficult-to-manufacturing titanium alloy, large cutting forces of which would intensify the vibration of the machining system. Self-excited vibration is one of the vibrations generated in the machining process of the hollow blade with HCSs. Self-excited vibration is well-known as the machining chatter, which is majorly caused by the unstable machining process. e machining chatter always results in some undesirable consequences such as poor surface finish and increased levels of tool wear [6][7][8][9]. e chatter occurring in the processing of the turbine blades with HCSs can even deteriorate the performance and life of the turbine engines. e most effective solution to avoid the chatter in machining of the hollow blade with HCSs is to make the machining process stable. e key to this solution is to find a stable machining zone and guide the process within this zone to achieve a stable output [10][11][12].
Dynamic parameters of the machining system are regularly analyzed and employed to predict the machining stability. Altintas and Budak [13] reported that the parameters of natural frequency, damping ratio, and stiffness can be used to express the dynamics of the machining system of the blade. e impact test and modal analysis are important analyzing methods for the dynamic parameters. e impact test is used to construct the relationship between the response displacement and the exciting force.
is relationship is normally expressed in a formula which is called a transferring function. e Fourier transform of this transferring function is known as the frequency response function (FRF). e fitting methods are widely used to calculate the dynamic parameters based on FRF. Liu et al. [14] used a fitting method to achieve the dynamic parameters of the tooltip from FRF. Sun and Altintas [15] and Song et al. [16] build FRF of the cutting tool by impact test and identified the modal parameters using the fitting method. From [13][14][15][16], it is known that the fitting method is easily performed in computing the modal parameters. However, the fitting method always tries to obtain a mean curve of FRF which would dismiss the peaks of the original FRF curve. e response displacements of these peaks are larger than those of the fitting FRF curve.
us, the computed modal parameters from the fitting FRF curve are incorrect which would result in a poor prediction of the machining system. e stability of the machining system can be predicted based on the dynamic parameters, on which much previous work has been carried out. e semidiscretization method (SDM) is a commonly used method to predict the machining stability, which was firstly proposed by Altintas and Budak [13] and has been developed and improved furtherly by many scholars. e zero-order Fourier term of the dynamic parameters is used to approximate the variation of the cutting forces and predict the stability of the machining system by drawing a stability lobe diagram (SLD). Budak and Altintas [17] developed a stability prediction model to increase the predicting accuracy of the stability for machining the flexible workpiece. Insperger and Stepan [18] proposed an improved SDM and build a SLD for the periodic machining system. Henninger and Eberhard [19] presented a method of increasing the efficiency of SMD by accelerating the computation of a transition matrix. Since SDM has many applications in predicting the stability of the machining process, limited work on predicting the machining stability of hollow blades has been reported. Zhou et al. [20] established a 3-dimensional SLD by using a finite element method (FEM) for machining of the hollow blade. Although FEM has well-known advantages in the experimental methods, the accuracy of FEM is affected by the boundary conditions of the machining system, which is difficult to be in accordance with the actual machining conditions. For the purpose of improving the system stability and avoiding the chatter in machining hollow blade with HCSs, this paper proposes a measurement-based approach for modal analysis and stability prediction of turn-milling hollow blade with HCSs. e impact test is carried out to achieve the FRF curves on the hollow blade and the cutting tool. An extremum method was employed to obtain an equivalent FRF curve, from which the modal parameters involving the natural frequency, damping ratio, and stiffness were computed. Afterwards, SDM was used to draw a SLD to predict the stability when machining hollow blades with HCSs. e paper is organized as follows. e setup of the impact test is presented in Section 2. Section 3 details the extremum method and an equivalent FRF curve is given. e stability prediction of machining hollow blade with HCSs is presented in Section 4, which has been validated experimentally. Conclusions are given in Section 5 finally.

Setup of Impact Test
An impact test was implemented to obtain the FRF curves of the milling tool and hollow blade, the setup of which is shown in Figure 2.
e hollow blade is fixed on a turnmilling center typed Mazak 200 IV ST with the designed fixtures, which were clamped by the chucks of the machining tool. Ten regions (①∼⑩) are distributed on the tailing and leading edges of the hollow blade, on which the hammer knocked. An accelerometer typed PCB 356A32 with a sensitivity of 12.09996 mV/m•s^2 was placed on the blade to inspect the response signals from each measured point, which was collected and analyzed by a Benstone Impaq Elite dynamic signal analyzer (cf. Figure 3). e sensitivity of the hammer was set as 2.5 mV/N. e signals collected in the signal analyzer were in the time domain, which were transformed by the Fourier transfer to obtain the FRF curves in the frequency domain.
It can be seen that vibration marks are the manifestation of the displacements of tool and blade on the surface of the blade in the machining process. Considering that the larger contribution to the vibration marks of the blade surface is the normal displacement of the blade, the measuring direction of the accelerator and the exciting direction of the hammer are both set in this normal direction (cf. Figure 4). It should be noted that the axial direction of the milling tool has a certain angle β with the tangential direction of this paper.

Modal Analysis on Machining System of Hollow Blade with HCSs
e machining system of turn-milling hollow blades with HCSs consists of the hollow blade, the cutting tool, and the fixture. e vibration marks on the blade surface are mainly caused by the normal displacement of the blade. e system is simplified as a single-degree-of-freedom system and the displacement in the normal direction along the machined surface of the hollow blade is considered. e dynamic model of the machining system involving the hollow blade and the milling tool can be formulated as follows [21]: where € x(t), _ x(t), x(t) are the acceleration, velocity, and displacement, respectively. m is the mass of the system, k is the stiffness, c is the damping, and F is the exciting force of the system. is equation can also be expressed as where the damping ratio ε and the natural frequency ω n are calculated as e harmonic force can be expressed as F(t) � F 0 e jwt , and the harmonic response is x(t) � Xe j(wt+∅) in equation (2). e equations can be transformed from the time domain to the frequency domain by using a Fourier transforming. By neglecting the effects of the initial conditions, equation (2) is transferred and expressed as follows: us, the FRF curve is expressed as

Shock and Vibration
en, the magnitude of the FRF is given by When the FRF is converted to the ratio of the acceleration to the force, the magnitude of the FRF is expressed as . e machining modes dominate the dynamic response of the entire system, which contains the tool and workpiece systems. e transfer function is expressed as where Φ t (ω) and Φ w (ω) are the transfer functions of the tool and workpiece systems, respectively. e magnitude of FRF curves on regions ①-⑤ of the hollow blade (as Figure 2), which were measured by the impact tests, is presented in Figure 5. e FRF curves on regions ①-⑤, which are shown with the magnitudes, are expressed as P1-P5 and the range of the frequency is 0-1000 Hz. e magnitudes of FRF curves on regions ①-⑤ are quite different which is due to the difference of the dynamic parameters of different regions on the hollow blade. It can be seen from Figure 5 that the magnitude of P4 is slightly larger than that of other regions, and the magnitude of P1 curve is the smallest. In the same method, the FRF curves on regions ⑥-⑩, which are shown with the magnitudes, are presented as P6-P10 in Figure 6. It is shown that the FRF curves have similar trends and two peaks of the magnitude can be found in the range of 0-1000 Hz. e natural frequencies corresponding to the peaks are approximately 300 and 600 Hz, respectively.
Modal analysis is a common method used to obtain the dynamic parameters of the machining system. e dynamic parameters of different regions on the hollow blade can be calculated based on the magnitude of FRF curves. By using the previous dynamic model, the dynamic parameters involving the natural frequency, damping ratio, and stiffness were calculated based on the FRF curves in Figures 5 and 6. e calculated stiffness at ten regions of the hollow blade is given in Figure 7. It is shown that the stiffness of P1-P5 and P6-P10 decreases gradually along the Z direction (cf. Figure 1), and the stiffness of P1-P5 is larger than that of P6-P10, respectively. According to equation (8), the higher the stiffness k is, the lower the amplitude of the response function is. It is obvious that the machining vibrations are easy to occur in the regions of the hollow blade near the right chuck. Besides, the stiffness of the hollow blade with a length of 100-350 mm in Z direction is much weaker than other regions, which is largely affected by the overhang length of the blade and the fixed position.
It is apparent that the dynamic parameters of different regions on the hollow blade are different. To simplify the estimation of the machining stability of the hollow blade, an equivalent FRF curve was obtained by using an extremum method and a fitting method which are shown in Figures 8  and 9, respectively. us, the FRF curves of different regions on the hollow blade are equivalent to one curve. e equivalent FRF curve obtained by the fitting method locates between the peaks and valleys of the FRF curves. e magnitudes of the FRF curve obtained by the fitting method are considerably lower than those of weaker regions on the hollow blade, which would decrease the calculating accuracy of the dynamic parameters and the stability prediction. e extremum method prefers using the peaks of the FRF curves which is sensitive to the weak region of the hollow blade. e extremum method is to compare the FRF curves of different regions and achieve the equivalent FRF curve by choosing the maximum value in a certain frequency. e equivalent curves of P1-P5 and P6-P10 are presented in Figure 8. It is shown that the magnitude of the equivalent FRF curve obtained with the extremum method is larger than that of the fitting method.
e frequencies corresponding to the amplitude of the equivalent FRF curves are the same for both the extremum method and the fitting method. A comparison of FRF curves shows that the equivalent FRF of P6-P10 is higher than that of P1-P5. e machining of the entire blade often requires a constant cutting parameter. us, the equivalent FRF curves of the hollow blade between P1 and P10 are obtained by the extremum method and the fitting method (cf. Figure 9). e equivalent FRF curves obtained are substituted into equation (4) to compute the modal parameters of the machining system. Table 1 shows the computed modal parameters for the fitting method and the extremum method, respectively. It is shown that the modal parameters of regions P1-P5 with two methods are close. e first-order natural frequencies obtained with the fitting method and the extremum method are 300 and 370 Hz, respectively, and the stiffnesses of two methods are 2.31 × 10 5 N/m and 1.77 × 10 5 N/m N/m. e natural frequency of P6-P10 is considerably lower than that of P1-P5 at 268 and 286 Hz, respectively. e first-order natural frequency and stiffness obtained by the extremum method are lower than those obtained with the fitting method. e first-order natural frequencies of the regions P1-P10, which are calculated by the extremum method and fitting method are 260 Hz and 308 Hz. e stiffness achieved by the fitting method is 7.2 × 10 4 N/m and the stiffness of the extremum method is 2.05 × 10 4 N/m. e stiffness of the extremum method is 71.5% smaller than that of the fitting method.
As illustrated in equation (9), the magnitude of the FRF is consistent with the magnitude of the blade and milling tool. e FRF curve of the milling tool is also measured by the impact test. Given that the magnitude of FRF is influenced by the angle β, the exciting and measuring directions on the tool were changed and the curve with smaller amplitude is shown in Figure 10. e FRF curve of the ballend milling cutter is compared with the equivalent FRF of the hollow blade obtained by the extremum method. It can be known from Figure 10 that the frequency response curve of the milling tool is considerably smaller than that of the blade. erefore, the modal parameters of the tool are negligible in predicting the machining stability. From equation (9), it is known that the higher the magnitude of FRF, the lower the stiffness is. So, the stiffness of the milling tool is higher than that of the hollow blade. e stiffness of the hollow blade is the main contributor to the vibration.

Stability Prediction and Experimental Verification
According to SDM presented by Altintas [1], the singledegree stability prediction model is expressed as where [A 0 ] is the directional cutting coefficient matrix, a is the depth of cut, K t is the tangential cutting force coefficient, and [∅(jω)] is the transfer function matrix. When the system is critically stable at the chatter frequency ω c , the root of the equation is formulated as follows: e critical axial depth of cut a lim expressed as where Λ R and Λ I stand for the real part and imaginary part of Λ. e SLD curve which contains the relationship between the critical axial depth of cut and the spindle speed can be obtained. For the differences of the dynamic parameters in ten regions, the SLD curves of ten regions are combined to be 3D SLD curves. e predicted 3D SLD curves for regions P1-P5 and P6-P10 are shown in Figures 11(a) and 11(b). e length in the figures represents the length of the blade along the Z direction, where the total length of the blade is 350 mm. e stiffness value k is relatively large which allows predicting the machining stability with the high ultimate depth of cut. e maximum critical depths of cut for P1-P5 and P6-P10 are 0.94 and 2.4 mm, respectively. e middle positions of the blade such as at P6, P7, and P8 have a low rigidity which results in a small critical depth of cut. e 3D SLD can predict the stability of different regions accurately. However, it will take too much time for the calculation of 3D SLD and the chatter will occur in the region where the critical axial depth of cut is the least. To improve the calculation efficiency, 2D SLD curves of the fitting method and extremum method are illustrated in Shock and Vibration 5 Figure 12.
e influence of the extremum and fitting methods on the stability prediction is subsequently compared. e SLD curves obtained from the extremum method on P1-P5, P6-P10, and P1-P10 and the fitting method on P1-P5 and P6-P10 are analyzed (cf. Figure 12). e predicted result of the stability on regions P1-P10 with the fitting method is higher than that of the extremum method. A comparison between the two methods shows similar results for P6-P10 and P1-P5. Figure 12 shows that the stiffness of P6-P10 is lower than that for P1-P5, and the SLD curve for P1-P10 has the smallest stable region. e cutting tests are carried out to verify the predicting results of stability on a Mazak 200 IV ST turn-milling center. A ball-end milling tool with a diameter D of 10 mm and a teeth number N of 2 is employed. e radial depth of cut of 0.2 mm, the feed rate of 500 mm/min, and the angle β of 45°a re used as the machining parameters. e reasonable axial cutting depth and spindle speed are selected to verify SLD       Shock and Vibration 7 and the spindle speed is in the range of 0-8,000 rpm. In order to confirm the reliability of SLD, the machining experiments with axial cutting depths of 0.03 mm, 0.1 mm, 0.2 mm, 0.3 mm, 0.4 mm, 0.5 mm, 0.6 mm are implemented when the spindle speed is 6,000 rpm. And the axial cutting depth is also set as 0.03 mm, 0.05 mm, 0.1 mm, 0.2 mm, and 0.3 mm when the spindle speed is 2,500 rpm. It can be known from the experiments that the axial cutting depth of 0.2 mm is the marginal zone in SLD. e machining process is unstable when the spindle speed is 6,000 rpm and axial cutting depth is higher than 0.2 mm. e opposite result will be received when the axial cutting depth is lower than 0.2 mm. e similar results can be achieved that the critical depth of cut is 0.05 mm when the spindle speed is 2,500 rpm. en, the parameters of spindle speeds of 2,500 and 6,000 rpm and the axial depth of cut of 0.3 and 0.03 mm are selected to verify that the extremum method is more reliable than the fitting method. When the spindle speed is 2,500 rpm, the prediction results of the extremum method and fitting method are both stable at the axial cutting depth of 0.03 mm, the prediction results of both methods are unstable at the axial cutting depth of 0.3 mm. e prediction results are considered to be different at the axial depth of 0.3 mm when the spindle speed is 6,000 rpm. e fitting method predicts that it is in a stable state, while the extremum method considers that the parameter is above the marginal of SLD, and the result of prediction is unstable. e conclusion can be collected that the predicting result of the extremum method is more accurate than the fitting method by comparison of SLD. Figure 13(a) shows the acceleration signal in the stable zone as collected with an acceleration sensor. e time domain acceleration signal for the unstable zone is also acquired and shown in Figure 13(b). e signals are compared when the axial depth of cut is 0.03 mm and 0.3 mm and the spindle speeds are 6000 rpm. It can be observed from Figure 13(b) that the acceleration magnitudes are higher than that of Figure 13(a) and the cutting process became unstable.

Conclusions
In order to avoid the chatter in turn-milling of hollow blades with HCSs, this paper presents a measurement-based approach for modal analysis and stability prediction of the turn-milling process of the hollow blades with HCSs. e FRF curves on the blade and the cutting tool were achieved by the impact test and were approximated as an equivalent FRF curve by using an extremum method, which is compared with the result by using a fitting method. e dynamic parameters involving the natural frequency, damping ratio, and stiffness are calculated based on the equivalent FRF curve and the dynamic model. e 2D and 3D stability lobe diagrams were drawn with the semidiscretization method to predict the machining stability of the hollow blade. e cutting tests were performed to validate the feasibility of the proposed approach. e results show that the equivalent FRF curve obtained by the extremum method is larger than that of the fitting method in the range of 0-1000 Hz. e acceleration signals with spindle speeds of 2,500 and 6,000 rpm and axial depth of cuts of 0.3 and 0.03 mm are analyzed. When the spindle speed is 2500 rpm, the predicting results of the fitting and extremum method are almost the same. When the axial depth of cut is 0.3 mm and the spindle speed is 6000 rpm, the predicting result of the extremum method is more consistent with the test result than that of the fitting method, which proves that the stability prediction using the extremum method is more reliable.

Data Availability
All the data regarding this paper are available, and the author can provide them if they are required.

Conflicts of Interest
e authors declare that they have no conflicts of interest.  Shock and Vibration