Shock and Vibration
Volume 2018 (2018), Article ID 3831825, 10 pages
https://doi.org/10.1155/2018/3831825

Research Article

Location-Dependent Prediction of Dynamic Stability Limit for Peripheral Milling of Surfaces with Variable Curvatures

Baosheng Wang,¹ Hongyan Hao,¹ Mulan Wang,² and Junming Hou²

¹Research Department of Intelligent Manufacturing Equipment, Nanjing Institute of Technology, Nanjing 211167, China
²Jiangsu Key Laboratory of Advanced Numerical Control Technology, Nanjing 211167, China

Correspondence should be addressed to Baosheng Wang; wbaosh@163.com

Received 13 September 2017; Revised 15 January 2018; Accepted 21 February 2018; Published 3 May 2018

Academic Editor: Mickaël Lallart

Copyright © 2018 Baosheng Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The stability limit may change with the cutter’s location due to effect of curvature during the milling of a complex surface. The method for calculating the actual radial cutting depth is presented by accounting for the effects of curvature on the actual cutting parameters. The computed radial cutting depth is in turn used to determine the entrance/exit angles. Moreover, a milling system dynamic model is established based on the instantaneous milling force coefficients, and the stability limit is determined by means of the time-domain semidiscretization method. In addition, a location-dependent method for predicting the stability associated with the peripheral milling of a complex surface is put forward and simulation is carried out to generate a stability limit diagram. The effectiveness of the proposed method is verified through milling tests.

1. Introduction

Cutting vibration may influence the quality and efficiency of the machining process or even lead to damage of cutters or milling machines. For this reason, it is desired to suppress the cutting vibration to the extent possible during the machining process. It is a traditional practice to select conservative processing parameters for preventing the cutting vibration, whose disadvantages include reduced production efficiency and machine utilization. Therefore, characterizing the milling system’s cutting vibration pattern and determining the range of operating conditions for ensuring the stable milling operation have become an important research area concerning the high speed milling.

Currently, scholars at home and abroad put forward many methods for predicting the chatter stability of a milling system. Altintaş and Budak first propose a single-frequency method [1, 2], which utilizes the zero-order term of the Fourier series associated with the coefficient matrix in the milling force direction to linearize the mathematical model so as to allow analytical calculation of the stability boundary. Merdol put forward a multifrequency method, which collectively takes into consideration a multitude of harmonics associated with the Fourier series of the coefficient matrix in the milling force direction [3]. Insperger presents a semidiscretization method in time domain, which, by discretizing the delay term, converts the delay differential equations governing the milling dynamics into a set of ordinary differential equations; this allows the construction of a system state transition matrix, based on which the stability of the system can be assessed using Floquet’s theory [4]; moreover, Zatarain employs this method to analyze the impact of milling cutter’s helix angle on the stability boundary by means of this method [5], and Long examines the effect of speed variation on the chatter stability [6]. Guo studies the milling stability limit associated with the multiregenerative effects based on the 3rd-order semidiscretization method [7]. Ding proposes a full-discretization method, which, by fully discretizing the delay and state terms, allows the application of Floquet’s theory to determine the stability limit of the system [8]. Ozoegwu carries out analysis concerning the impact of the discretization order on the prediction accuracy for the full-discretization method [9]. Zatarain proposes an implicit subspace iteration approach to enhance the efficiency of calculating the stability limit diagram, which is compatible with both semidiscretization and full-discretization methods in time [10]. Constant milling force coefficients are utilized for calculating the dynamic milling force in [1–10]; hence, the effect of feed per tooth on the chatter stability is neglected. Landers establishes a dynamic model based on the instantaneous milling force coefficients in exponential form so as to analyze the effect of feed per tooth on the chatter stability [11]; moreover, Ahmadi studies the linear milling stability with varying radial depth of cut based on the instantaneous milling force coefficients in exponential form [12].

It should be noted that the previous studies cited above [1–12] are all focused on linear milling. Kardes studies the stability limit of the circular milling process using single-frequency method [13], and Peng analyzes the chatter stability during a circular corner milling process based on the time-domain simulation method [14]. Unlike the linear and circular milling, during the peripheral milling of a surface with variable curvature, tool feed direction and the actual radial cutting depth are changing continuously along tool path. Then, milling dynamics and stability characteristics are also changing along tool path. Yang presented a method to predict chatter for the peripheral milling of thin-walled workpieces with curved surfaces [15]. In the paper, the milling force coefficients are simplified as constants. So, the feed per tooth cannot be included in dynamic model, and limit of feed per tooth cannot be predicted.

In the present study, effective radial cut depth is proposed, and a new milling force model is established. With the milling force model, a dynamic model of cutting tools is constructed. Subsequently, the time-domain semidiscretization method is employed to calculate the stability limit. Considering that the actual radial cutting depth and tool feed direction are changing along tool path, a new method for location-dependent prediction of dynamic stability limit for peripheral milling of surfaces with variable curvatures is presented. Moreover, the proposed algorithm is experimentally validated. The advantage of the proposed approach lies in that stability during milling process can be predicted. Thus, stable milling will be ensured.

2. Dynamic Model for Peripheral Milling of Complex Surface

2.1. Geometric Description of Peripheral Milling of Complex Surface

The dynamic model of a milling system for peripheral milling of a complex surface is schematically shown in Figure 1, where , , and represent the parameterized tool path, contour prior to processing, and target contour, respectively. and denote the fixed and follow-up coordinate systems, respectively. represents the tool radius. represents the nominal radial depth of cut. represents the instantaneous feed angle.

Figure 1: Schematic diagram of the proposed dynamic model describing the peripheral milling system for complex surface machining.

Regarding the parametric geometry of the target contour for reference, the tool path and contour prior to processing are expressed as

where and represent the derivatives of the parametric curve, respectively.

During a machining process, the congruent relationship between the machining time and the parameter is calculated by

where represents the spindle speed. represents the nominal feed per tooth. represents the number of cutter teeth. represents the parameter spacing. represents the spacing between two adjacent parameter points on a tool path.

2.2. Dynamic Model

Our model, as schematically shown in Figure 2, consists of a rigid workpiece and a tool vibration system with two structural degrees of freedom. By accounting for only a single vibration mode in lieu of the two coupled vibration modes associated with the motions in two directions, the tool vibration can be described by the following differential equation:

where ; ; ; ; and . , , and represent the modal mass, damping coefficient, and stiffness in the direction, respectively. Similarly, , , and represent the modal mass, damping coefficient, and stiffness in the direction, respectively. and represent the instantaneous position at time () in the fixed coordinate system (). and represent the instantaneous milling forces in the and directions, respectively.

Figure 2: Schematic diagram of the proposed dynamic model describing a peripheral milling system.

By taking the tool’s helix angle into account and dividing the milling cutter equally into an array of infinitesimal elements along the axial direction, the milling force associated with the element at time can be expressed as

where and represent the tangential and radial milling forces, respectively. represents the instantaneous undeformed chip thickness. and represent the tangential and radial milling force coefficients, respectively. represents the axial length, and represents the number of elements per tooth.

The tool vibration typically leads to complicated chatter marks on the surface of the workpiece, and therefore the undeformed chip thickness for the next tooth may be influenced. For this reason, the instantaneous undeformed chip thickness needs to include two components, namely, the static undeformed chip thickness and the regenerative thickness due to cutting tool vibration, which can be calculated by the following equations:

where and represent the static undeformed chip thickness and the regenerative undeformed chip thickness, respectively; represents the Heaviside step function; represents the angular position at time in the follow-up coordinate system ; represents the tool’s helix angle; and represents the angular position at time in the follow-up coordinate system .

It is found in our previous studies that the effect of the curvature on the actual feed per tooth is negligibly small during the peripheral milling process on a surface with variable curvature, and therefore the time-dependent feed per tooth can be approximated using a constant nominal feed per tooth [16]. Thus, the static undeformed chip thickness may be expressed as

Figure 2 indicates that the regenerative thickness corresponds to the difference between the tool’s dynamic displacements at time (Point ) and time (Point ), which may be described as

where represents the regenerative delay, and one can get by neglecting the effect of the surface curvature on the delay. represents the angular position at time in the fixed coordinate system (), and .

During the peripheral milling process on a surface with variable curvature, the variation of the actual radial depth of cut causes the entrance/exit angles of the cutter teeth to change accordingly. Whether an element is involved in milling can be determined by the following equations:

where and represent the instantaneous entrance and exit angles, respectively. represents the instantaneous actual radial depth of cut, which corresponds to DH in Figure 3. and represent the coordinates of Point .

Figure 3: Actual radial cutting depth.

Based on [16], the instantaneous milling force coefficients can be expressed as

where , , , and are constants.

Substituting (5) and (13) into (4), milling force of the element can be expressed as

Then, linearizing the milling force by applying Taylor expansion to the infinitesimal elements while neglecting the high order terms and the static part yields

where and .

By expressing in the fixed coordinate system () the cutting force associated with the infinitesimal elements and conducting summation over all the elements, one can obtain the dynamic milling force:

where represents the direction coefficient matrix. , , , and represent the direction coefficients which are related to the instantaneous tool position and are time-dependent while milling a surface with varying curvature.

3. Analysis of Stability Based on the Time-Domain Semidiscretization Method

The time-domain semidiscretization method is based on discretizing the delay terms and converting the delay differential equation into a first-order differential equation in the state space, which is in turn solved [4]. Fully accounting for the cutting and noncutting stages of the entire milling process, this method can be widely applied to those cases with various radial depths of cut. Also, it can greatly facilitate the milling stability analysis with varying parameters.

Substituting the dynamic milling force derived from (15) into (3), the delay differential equation of the milling system can be described as

By assuming , the state equation is expressed as

where represents the system state matrix. represents the system input matrix. ; ; and .

The direction coefficient matrix is time-dependent for the variable-curvature peripheral milling process; thus, the system is a linear time-dependent system. The delay period is discretized into small time elements (). As for the sampling time (), , , and represent the corresponding system state, as well as the state and input matrices, respectively; and represents . If is small enough, , and (18) is expressed as

The solution to (19) is as follows:

where is constant and depends on the initial condition of vibration. For , the following equation is obtained:

Equation (21) can be rewritten as

For , the following equation is obtained:

Equation (23) can be expressed in a matrix form as

where

In accordance with the recursive principle, the transfer matrix () of the system within a single period () can be constructed from ():

According to Floquet’s theory, the system’s stability state depends on the modulus of the maximum eigenvalue of ; that is,

4. Location-Dependent Prediction of Stability Limit and Experimental Validation

Milling system parameters are as follows:

 Milling machine: DMC 70V Hi-Dyn Toolholder: HSKA63-ER25-100 Cutter: Sandvick 1P240-1000-XA1630

 Diameter: 10 mm Number of teeth: 4 Helix angle: 35°

 Cutter’s exposed length postinstallation: 45 cm

With hammering method as shown in Figure 4, the experimental results with DH5920 are shown in Figure 5. The modal parameters of the present tool system can be obtained as follows:

Figure 4: Test of modal parameters of tool system.

(a) Frequency response functions in direction

(b) Time response curve in direction

(c) Frequency response functions in direction

(d) Time response curve in direction

Figure 5: Results of modal experiments with hammering method.

Workpiece material is Al6061-T6.

Experimental identification of the milling force coefficient is linear milling.

Testing parameters are as follows:

 Spindle speed: 3000 r/min Axial depth of cut: 5 mm Nominal radial depth of cut: 3 mm Feed per tooth: 0.04 mm/r

Milling mode is dry cutting and down milling.

With the method mentioned in [17], milling force coefficients can thereby be calculated as follows:

4.1. Stability Limit Diagrams

During a variable-curvature peripheral milling process, the feed direction and the actual radial depth of cut are subject to a constant change due to the variation of surface curvature, which in turn causes the associated stability limit to vary accordingly. For a given set of nominal milling parameters, it is possible to encounter a situation where the system is stable during certain periods of milling while experiencing vibration during some other periods of milling. Suppose that the workpiece’s contour can be expressed using the parametric equation below:

The milling parameters are set as  mm/r and  mm. Based on point-wise scanning method [12], one can compute the stability limits associated with the cutter locations of (, 0.3, 0.7, and 0.9); the corresponding feed direction angles and the actual radial depths of cut are listed below:

The results are shown in Figure 6. A comparison of results between and indicates that, for the same feed direction angle, the surface curvature affects the effective radial depth of cut, which in turn affects the system’s stability limit. Moreover, a comparison of results between and shows that their stability limit diagrams differ dramatically despite the fact that their effective radial depths of cut are quite close, which indicates that the feed direction also affects the stability to some extent.

Figure 6: Stability limit diagram when the tool cutting tool position.

For a given set of spindle speed, nominal radial depth of cut, and feed per tooth, one can place the cutting tool at different locations and use (27) to calculate the corresponding maximum radial depths of cut based on instantaneous effective radial depth of cut and feed direction, with the results shown in Figure 7 ( mm/r,  mm, and ). It can be seen in the figure that, for the given milling parameters, the axial depths of cut must be less than 2.8 mm to ensure the stability during the entire milling process.

Figure 7: Maximum axial cutting depth versus tool position.

Similarly, the maximum nominal radial depths of cut corresponding to various cutter locations are calculated for a given set of spindle speed, axial depth of cut, and feed per tooth, with the results shown in Figure 8 ( mm/r,  mm, and  r/min). It can be seen in the figure that the maximum nominal radial depth of cut must be less than 3.6 mm to ensure the stability during the entire milling process.

Figure 8: Maximum nominal radial cutting depth versus tool position.

For a given combination of spindle speed, axial depth of cut, and nominal radial depth of cut, the maximum feeds per tooth corresponding to various cutter positions are calculated, with the results shown in Figure 9 ( r/min,  mm, and  mm). It can be seen in the figure that the feed per tooth must be more than 0.028 mm/r to ensure the stability during the entire milling process.

Figure 9: Maximum feed per tooth versus tool position.

4.2. Experimental Validation

For verifying the correctness of location-dependent prediction of the variable-curvature milling stability, a test is designed in accordance with the predicted results (Figure 7). The experimental scene is shown in Figure 10. The test parameters are  mm/r,  mm, and  r/min, and the axial depths of cut are  mm and  mm, respectively. The maximum axial cutting depth reaches its lowest level for a cutter position of and the corresponding machining duration is 1.95 s. Thus, analysis is carried out concerning the signals of -direction milling force near  s and the results are shown in Figure 11. As we know, when milling is stable, frequency arising in frequency spectrum of cutting force signal is cutting frequency or it is harmonic frequency. According to the milling parameters, cutting frequency is 320 Hz. It can be seen in Figure 11 that the milling stability is achieved at those cutting frequencies marked in the frequency spectrum of milling force (i.e., 320 Hz and its multiplications) when the axial depth of cut is 2.7 mm. When the axial depth of cut is 2.9 mm, it is obvious that the cutting vibration frequencies are 4147 Hz, 4707 Hz, and 5027 Hz. Their spacing is the integer multiples of the cutter’s rotating frequency. This indicates that an unstable stage is encountered during the milling process under the given machining conditions. The measured results at two axial depths of cut indicate that our method for conducting location-dependent prediction concerning the stability of peripheral milling over surfaces with variable curvatures is accurate and feasible.

Figure 10: Milling tests.

(a) Milling force in -direction (ap = 2.7 mm)

(b) Milling force in -direction (ap = 2.9 mm)

(c) Frequency spectrum of -direction milling force (ap = 2.7 mm)

(d) Frequency spectrum of -direction milling force (ap = 2.9 mm)

Figure 11: Frequency spectrum of -direction milling force.

It should be noted that the cutter runout has a certain effect on milling stability, because runout makes the uncut chip thickness deviate from the nominal value. Wan has analyzed the effect of runout on milling stability in [18]. In this paper, the main purpose is to provide a method for location-dependent prediction of dynamic stability. To simplify the dynamic model solution, the runout is neglected.

5. Conclusions

By considering the impact of surface curvature on the effective radial depth of cut, the present study proposes a method for computing the effective radial depth of cut, which can be employed to rapidly determine the entrance/exit angles based on the effective radial depth of cut. As a result, the computational overhead for predicting the milling force is greatly reduced, and the computational efficiency is substantially enhanced; moreover, a dynamic milling model is established based on the instantaneous milling force coefficients, which is solved using the time-domain semidiscretization method to determine the milling stability; in addition, a method for conducting location-dependent prediction of stability related to the peripheral milling of complex surfaces is put forward by fully taking into consideration the complex surface milling stability’s dependency on cutter location. Furthermore, simulation concerning the peripheral milling of complex surfaces is carried out, based on which the stability limit diagrams pertaining to axial and radial depths of cut as well as the feed per tooth are plotted. Also, the effectiveness of the proposed method is validated with the milling experiments.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work was financially supported by the Joint Research Project of Jiangsu Province, China (BY2016008-04), the University Science Research Project of Jiangsu Province, China (16KJD460003), and the Scientific Foundation of Nanjing Institute of Technology (ZKJ201201).

References

Y. Altintaş and E. Budak, “Analytical prediction of stability lobes in milling,” CIRP Annals - Manufacturing Technology, vol. 44, no. 1, pp. 357–362, 1995.

E. Budak and Y. Altıntaş, “Analytical prediction of chatter stability in milling—part I: general formulation,” Journal of Dynamic Systems, Measurement and Control—Transactions of the ASME, vol. 120, no. 1, pp. 22–30, 1998.

S. D. Merdol and Y. Altintas, “Multi frequency solution of chatter stability for low immersion milling,” Journal of Manufacturing Science and Engineering, vol. 126, no. 3, pp. 459–466, 2004.

T. Insperger and G. Stepan, “Semi-discretization method for delayed systems,” International Journal for Numerical Methods in Engineering, vol. 55, no. 5, pp. 503–518, 2002.

M. Zatarain, J. Muñoa, G. Peigne, and T. Insperger, “Analysis of the influence of mill helix angle on chatter stability,” CIRP Annals - Manufacturing Technology, vol. 55, no. 1, pp. 365–368, 2006.

X. Long and B. Balachandran, “Stability of up-milling and down-milling operations with variable spindle speed,” Journal of Vibration and Control, vol. 16, no. 7-8, pp. 1151–1168, 2010.

Q. Guo, Y. Sun, Y. Jiang, and D. Guo, “Prediction of stability limit for multi-regenerative chatter in high performance milling,” International Journal of Dynamics and Control, vol. 2, no. 1, pp. 35–45, 2014.

Y. Ding, L. Zhu, X. Zhang, and H. Ding, “A full-discretization method for prediction of milling stability,” The International Journal of Machine Tools and Manufacture, vol. 50, no. 5, pp. 502–509, 2010.

C. G. Ozoegwu, S. N. Omenyi, and S. M. Ofochebe, “Hyper-third order full-discretization methods in milling stability prediction,” The International Journal of Machine Tools and Manufacture, vol. 92, pp. 1–9, 2015.

M. Zatarain, J. Alvarez, I. Bediaga, J. Munoa, and Z. Dombovari, “Implicit subspace iteration as an efficient method to compute milling stability lobe diagrams,” The International Journal of Advanced Manufacturing Technology, vol. 77, no. 1-4, pp. 597–607, 2015.

R. G. Landers and A. G. Ulsoy, “Nonlinear feed effect in machining chatter analysis,” Journal of Manufacturing Science and Engineering, vol. 130, pp. 1707-1708, 2008.

K. Ahmadi and F. Ismail, “Modeling chatter in peripheral milling using the Semi Discretization Method,” CIRP Journal of Manufacturing Science and Technology, vol. 5, no. 2, pp. 77–86, 2012.

N. Kardes and Y. Altintas, “Mechanics and dynamics of the circular milling process,” Journal of Manufacturing Science and Engineering, vol. 129, no. 1, pp. 21–31, 2007.

C. Peng, L. Wang, Z. Li, and Y. Yang, “Time-domain simulation and experimental verification of dynamic cutting forces and chatter stability for circular corner milling,” Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, vol. 229, no. 6, pp. 932–939, 2015.

Y. Yang, W.-H. Zhang, Y.-C. Ma, and M. Wan, “Chatter prediction for the peripheral milling of thin-walled workpieces with curved surfaces,” The International Journal of Machine Tools and Manufacture, vol. 109, pp. 36–48, 2016.

B. Wang, H. Hao, M. Wang, J. Hou, and Y. Feng, “Identification of instantaneous cutting force coefficients using surface error,” The International Journal of Advanced Manufacturing Technology, vol. 68, no. 1–4, pp. 701–709, 2013.

H. Hao, B. Wang, and W. Tang, “Prediction of instantaneous milling force taking runout into account in peripheral milling of curved surface,” The International Journal of Advanced Manufacturing Technology, vol. 79, no. 1-4, pp. 49–56, 2015.

M. Wan, W. H. Zhang, J. W. Dang, and Y. Yang, “A unified stability prediction method for milling process with multiple delays,” in Proceedings of the International Journal of Machine Tools Manufacture 50, vol. 50, pp. 29–41, 2010.