Published: 30 June 2017

# Reliability of chatter stability in CNC turning process by Monte Carlo simulation method

Haitao Zhu2
Hassan A. Wahab3
1College of Mechanical and Electrical Engineering, Harbin Engineering University, No. 145 Nantong Street, Harbin 150001, China
2College of Ship Building Engineering, Harbin Engineering University, No. 145 Nantong Street, Harbin 150001, China
3, 1Faculty of Engineering and Technical Studies, University of Elimam Elmahdi, Kosti, Sudan
Corresponding Author:
Views 120

#### Abstract

In turning process, regenerative chatter stability is regarded as an outcome favor in achieving reliable cutting performance. Since nowadays there is great demand in producing very high quality parts, researchers devoted great efforts in developing theoretical and analytical means for understanding, analyzing and solving the stability of a given machining process. In general, the stability of turning process could be influenced by random variables such as machine stiffness, damping and cutting force. Therefore, in this paper, a new method for predicting the reliability of regenerative chatter stability in turning process is proposed. This involves application of Monte Carlo simulation (MCS) tool to study the statistical distribution of regenerative chatter stability during turning process. The employed dynamic model of the regenerative chatter was defined, and the relationship between cutting depth and spindle speed was obtained by Laplace transformation. The sampling distribution information of kinetic parameters of the cutting system was evaluated for the regenerative chatter model. Thus, the limit cutting width of the samples was obtained. The statistical characteristics probability showed that the actual cutting width is less than the limit cutting width. Therefore, the regenerative chatter stability prediction by Monte Carlo simulation was regarded as a criterion.

## 1. Introduction

In numerical control turning operations, cutting vibration seriously contributes in affecting the workpiece machineability, and hence, lowers the finished part quality. Previous research demonstrated that the cutting vibration is often caused by regenerative feedback inside the system which is termed as self-excited vibration or cutting chatter [1]. In general, chatter in machining process can be categorized as regenerative cutting chatter, vibration mode coupled cutting chatter and negative friction cutting chatter. It is well known that regenerative chatter is the most important form of vibration which deleteriously affecting the machining process [2, 3]. Usually for solving the machining chatter instability issue, dynamic model analysis of the system, which includes the influence of axial cutting force and gyroscopic effect on chatter stability, must be introduced. Then dynamic cutting experiments should be performed in order to verify the correctness of the theoretical model [4-6]. Similar procedure was employed in the chatter prediction of milling operation using continuous wavelet features and multi class ball support vector machines (SVM) where stiffness and damping were optimized [7]. In other research works Y. Altintas et al. [8] investigated the influence of the vibration coupling on the machining stability of CNC machine tool subsystem. P. Albertelli et al. [9] established chatter stability analysis model for heavy duty milling machine to optimize the process planning and tool selection process. However, the stability analysis model of traditional cutting chatter ignores the influence of uncertain factors, which will inevitably lead to certain analysis errors. In machining practice, usually there are few deterministic factors, most of which are stochastic factors inherent in the process of metal cutting and cannot be ignored. Accordingly, the practical machining process is very likely to have performance deviation, and it is often difficult to achieve the exact desired effect. In this study, we intended to take into consideration the effect of random factors in the stability of the CNC turning operation. Therefore, a generalized method for predicting the reliability of chatter stability for this process is intended to be achieved by using Monte Carlo numerical simulation.

## 2. Analysis of the regenerative chatter

The regenerative chatter analysis of the turning process in question is based on the difference between cutting thickness caused by the phase difference and the vibration displacement [10]. Fig. 1 represents the dynamic model of regenerative chatter system for this process. Accordingly, the differential equation of vibration [1, 3] is given as follows:

1
$m\stackrel{¨}{q}\left(t\right)+c\stackrel{˙}{q}\left(t\right)+kq\left(t\right)=∆{F}_{d}\left(t\right)\mathrm{cos}\left(\beta -\alpha \right),$

where $m$ is the equivalent mass of lathe tool system, $c$ is the equivalent resistance, $k$ is the equivalent stiffness, $\beta$ is the angle between the dynamic cutting force, $\mathrm{\Delta }{F}_{d}\left(t\right)$ is the tool vibration direction, $\alpha$ is the angle between the main vibration direction and $q\left(t\right)$ is the direction of tool vibration.

Fig. 1Dynamic and inertia model of the vibrated system

The dynamic cutting force $\mathrm{\Delta }{F}_{d}\left(t\right)$ is generally caused by the variation in the cutting thickness, and can be expressed by the follows relationship:

2
${∆F}_{d}\left(t\right)={k}_{c}{b}_{a}\left(t\right)={k}_{c}b\left[\mu y\left(t-T\right)-y\left(t\right)\right],$

where $b$ is the cutting width, $kc$ is the stiffness coefficient per unit width of cut, $y\left(t\right)=q\left(t\right)\mathrm{c}\mathrm{o}\mathrm{s}\alpha$ is the cutter vibration displacement, $0<\mu <1$ is the overlap coefficient of twice cutting.

By substituting Eq. (2) into Eq. (1) the differential equation of regenerative chatter can be rearranged as follows:

3
$\stackrel{¨}{y}\left(t\right)+2{\omega }_{n}\zeta \stackrel{˙}{y}\left(t\right)+{\omega }_{n}^{2}y\left(t\right)=K\left[y\left(t\right)-\mu y\left(t-T\right)\right],$

where ${\omega }_{n}$ is the natural angular frequency of the vibrated system, ${\omega }_{n}^{2}=k/m$, $\zeta$ is the damping ratio, $\zeta =c/\left(2m{\omega }_{n}\right)$ and $u=\mathrm{c}\mathrm{o}\mathrm{s}\left(\beta –\alpha \right)\mathrm{c}\mathrm{o}\mathrm{s}\alpha$ is the directional coefficient.

The dynamic response of turning vibration is given by the following expression:

4
$y\left(t\right)=A\mathrm{sin}\left(\omega t+\phi \right).$

Substitution Eq. (4) into Eq. (3), we can obtain expression for Eq. (5) as follows:

5
$\stackrel{¨}{y}\left(t\right)+\left(2{\omega }_{n}\zeta -K\mu \mathrm{sin}\left(\omega T\right)/\omega \right)\stackrel{˙}{y}\left(t\right)+\left[{\omega }_{n}^{2}-K\left(1-\mu \mathrm{cos}\left(\omega T\right)\right)\right]y\left(t\right)=0.$

By employing Laplace transformation to Eq. (5), we can obtain the characteristic equation of the chatter system as follows:

6
${s}^{2}+\left(2{\omega }_{n}\zeta -K\mu \mathrm{sin}\left(\omega T\right)/\omega \right)s+{\omega }_{n}^{2}-K\left(1-\mu \mathrm{cos}\left(\omega T\right)\right)=0.$

Hence, the stability state of the vibrated system can be judged by the root of Eq. (6) where $s=\sigma +i\omega$. Then if the real part of $s$ is $\sigma >0$, the turning process system is in unstable state. However, when $\sigma <0$ the turning process is considered to be in stable state. Otherwise, the turning process system is in limiting state when $\sigma =$ 0. Therefore, in Eq. (6) the value of $s$ was replaced by $i\omega$ and the equation became as follows:

7
${\omega }_{n}^{2}-{\omega }^{2}+2{\omega }_{n}\omega \zeta i=K\left(1-\mu \mathrm{cos}\left(\omega T\right)\right)+K\mu \mathrm{sin}\left(\omega T\right)i.$

If the real parts in the left and right sides of the Eq. (7) are feasible, then should be respectively equal. Thus, the following expressions were concluded:

8
$\left\{\begin{array}{l}{\omega }_{n}^{2}-{\omega }^{2}=K\left(1-\mu \mathrm{cos}\left(\omega T\right)\right),\\ 2{\omega }_{n}\omega \zeta =K\mu \mathrm{sin}\left(\omega T\right).\end{array}\right\$

Assuming that that $\lambda =\omega /{\omega }_{n}\text{,}$ then Eq. (8) can be further converted to the following expression:

9
$\left\{\begin{array}{l}{1-\lambda }^{2}=-{k}_{c}{b}_{u}\left(1-\mu \mathrm{cos}\left(\omega T\right)/k,\\ 2\zeta \lambda =-Kb\mu u\mathrm{sin}\left(\omega T\right)/k.\end{array}\right\$

By Solving Eq. (9), the limit cutting width of machine tool is obtained:

10
${b}_{lim}=\frac{1-\mu \mathrm{cos}\left(2j\pi +\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{s}\mathrm{i}\mathrm{n}\frac{2\zeta \lambda }{\mu \sqrt{{\left(2\mu \lambda \right)}^{2}+{\left(1-{\lambda }^{2}\right)}^{2}}}-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{a}\mathrm{n}\left(\frac{2\zeta \lambda }{1-{\lambda }^{2}}\right)\right)}{1+{\mu }^{2}-2\mu \mathrm{cos}\left(2j\pi +\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{s}\mathrm{i}\mathrm{n}\frac{2\zeta \lambda }{\mu \sqrt{{\left(2\mu \lambda \right)}^{2}+{\left(1-{\lambda }^{2}\right)}^{2}}}-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{a}\mathrm{n}\left(\frac{2\zeta \lambda }{1-{\lambda }^{2}}\right)\right)}$
$\bullet \frac{k\left[{\left(1-{\lambda }^{2}\right)}^{2}+{\left(2\zeta \lambda \right)}^{2}\right]}{{k}_{c}u\left({\lambda }^{2}-1\right)},\left(j=0,1,\dots \right).$

And the corresponding spindle speed $\left(N\right)$ is given by the following relationship:

11
$N=\frac{60\lambda {\omega }_{n}}{2j\pi +\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{s}\mathrm{i}\mathrm{n}\frac{2\mathit{\zeta \lambda }}{\mu \sqrt{{\left(2\mathit{\mu \lambda }\right)}^{2}+{\left(1-{\lambda }^{2}\right)}^{2}}}-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{a}\mathrm{n}\left(\frac{2\zeta \lambda }{1-{\lambda }^{2}}\right)},\left(j=0,1\dots \right).$

By substituting the dynamic parameters $k$, $c$, $m$, $\alpha$, $\beta$, $kc$, and $\mu$ in Eqs. (10) and (11), the limit cutting width for the stable machining ${b}_{lim}$ and the corresponding spindle speed $N$ at $j=$ 0, 1,…, $n$, can be obtained.

By substituting $\mu =$ 1 into Eq. (9), with a larger overlap coefficient of cutting, the chatter can be easily generated, the minimum limiting width of cut for stable machining is expressed as follows:

12
${\left({b}_{lim}\right)}_{\mathrm{m}\mathrm{i}\mathrm{n}}=\frac{2k\zeta \left(1+\zeta \right)}{{k}_{c}u}.$

## 3. Reliability expressions for Monte Carlo simulation

In general, mechanical reliability refers to the ability of mechanical products to complete the specified function in the exact required service period, and it can be a measurable value. The two-common method for measuring the reliability are the moments method and Monte Carlo simulation. Comparing Monte Carlo simulation with moment method, the former one has the advantage of high universality and high precision [11]. However, Monte Carlo method is based on a large number of simulation and calculations, which is not suitable for predicting the reliability of complex numerical models as finite element method. As the case in this study does not involve complicated numerical model, Monte Carlo simulation is an appropriate tool to be used here, since numerical control cutting width is an explicit expression. Therefore, the chatter stability of the CNC turning operation was predicted and compared in the term of cutting width. Thus, the reliability of the turning system stability can be described as a multidimensional integral, given by:

13
${p}_{s}=\int ...{\int }_{{}_{b<{b}_{\mathrm{l}\mathrm{i}\mathrm{m}}\left(x\right)}}^{}fx\left({x}_{1},{x}_{2},...,{x}_{n}\right)d{x}_{1}d{x}_{2}...d{x}_{n},$

where ${x}_{1}$, ${x}_{2}$,..., ${x}_{n}$ are random design variables ($k$, $\zeta$, ${\omega }_{n}$, $\mu$, $u$, ${k}_{c}$), $n$ is the number of random design variables, $fx$(${x}_{1}$, ${x}_{2}$,…, ${x}_{n}$) the joint probability density function of random design variables.

If the design variables are independent of each other, Eq. (13) can be rewritten as follows:

14
${p}_{s}=\int ...{\int }_{b<{b}_{\mathrm{l}\mathrm{i}\mathrm{m}}\left(x\right)}^{}f{x}_{1}\left({x}_{1}\right),f{x}_{2}\left({x}_{2}\right)...f{x}_{n}\left({x}_{n}\right)d{x}_{1}d{x}_{2}...d{x}_{n},$

where $f{x}_{i}\left({x}_{i}\right)\left(i=1,2,...,n\right)$is probability density function of ${x}_{i}$ of random variables.

The basic idea is that the reliability of turning chatter stability indeed can be solved by using MCS method where the $N$ groups of random samples ${x}_{j}$ ($j=$ 1, 2,…, $N$) are generated from the probability density function (PDF) ${f}_{xi}\left({x}_{i}\right)$ of random variable ${X}_{i}$. The $N$ random samples are substituted into the function expression of limiting width of cut for stable machining as represented by Eq. (10). Then, it is counted for the number of sample points $Nr$ that ${b}_{lim}$ is greater than the actual cutting width $b$, and the reliability $Ps$ is approximately replaced by using a frequency that is corresponding to $Nr/N$ as in the following expression:

15
${p}_{s}=\int ...{\int }_{{R}^{n}}^{}{I}_{r}\left(x\right)fx\left({x}_{1},{x}_{2},...,{x}_{n}\right)d{x}_{1}d{x}_{2}...d{x}_{n}=E\left({I}_{r}\left(x\right)\right)\approx {N}_{r}/N,$

where:

${I}_{r}\left(x\right)=\left\{\begin{array}{l}1,\left(x\in R\right),\\ 0,\left(x\notin R\right).\end{array}\right\$

Is the indicator function in the security domain ${R}^{n}$; $n$ is dimensional variable space; and $E\left(\right)$ represents of the mathematical.

## 4. Calculation and estimation of chatter stability

Fig. 2 reveals the followed steps to obtain the estimated numerical values for the reliability of chatter stability as function of cutting width and the machine spindle speed. In the first step, we need to determine the number of simulation $N$. Secondly the $N$ groups of samples are extracted according to the information distribution of the design parameters. Thirdly, the machining stability is changed by the sample generation analysis of the model to calculate the ultimate cutting width (${b}_{lim}$). In the final step, the reliability of the machine stability is calculated. This requires taking into consideration the characteristic parameters of the CNC lathe machine during cutting process Table 1. Therefore, the characteristics parameters values of the lathe machine were substituted in Eqs. (9), (10) and (11) respectively in order to obtain the output data for cutting width with the spindle speed. After then the results were plotted in the lobe graph of the CNC lathe as shown in Fig. 3. It can be seen that the limiting cutting width is varied with the spindle speeds, and the minimum limiting cutting width (${b}_{lim}$) min is equal to 0.36 mm. Then Monte Carlo simulation method was used to estimate the reliability of the process chatter stability, which is mainly based on statistical analysis of the data. Generally, the variance of random parameters can be obtained by the reliability testing or the statistical analysis of experimental data. However, in the current case there is no experimental data, and thus, the variance of the random parameters can be determined from the variation coefficient $\left(C\right)$ which is set to be equal to 0.05. In this case, the random parameters are said to obey the statistical normal distribution, if a large number of independent factors [12, 13] influence the random variables.

By following the steps in Fig. 2 using 105 Monte Carlo sampling simulations, the limit cutting widths and the equivalent reliability of chatter stability were estimated. Fig. 4 shows the result obtained from the simulation. It can be seen that the very good stability of the machine chatter was obtained at cutting width equal to 0.28 mm, while the poorest machine stability occurs at cutting width of (${b}_{lim}$) where the equivalent reliability is only 52 % (Table 2).

Table 1Characteristic parameters of CNC lathe machine

 $K$ (N/mm2) $C$ (Ns/mm) $m$ (Ns2/mm) $\alpha$ (degree) $\beta$ (degree) $Kc$ (N/mm2) $\mu$ 3100 0.6 0.01 60 45 2018 1

Table 2Calculated reliability values by the advanced MCS methods

 Cutting width (mm) Reliability (%) Chatter stability indication 0.28 98 Good stability 0.32 84 Low stability ${\left({b}_{lim}\right)}_{\mathrm{m}\mathrm{i}\mathrm{n}}$ 52 Poor stability

Fig. 2Flowchart for calculating chatter stability reliability

## 5. Conclusions

Based on Monte Carlo simulation, the reliability of regenerative chatter stability was calculated. In Comparison with the traditional cutting force prediction method, Monte Carlo simulation take into consideration the influence of the random parameters such as stiffness, damping, mass and cutting stiffness coefficient on the cutting process stability. It can well predict with very close approximation to the actual situation of the process. Furthermore, different spindle speed for turning process can be predicted where good machine stability is attained. This study provided information on how to improve the machining accuracy, ensure long tool life and increasing production rate. The output of this study is expected to have practical application in manufacturing engineering.

Fig. 3Stable lobe graph of CNC turning operation

Fig. 4Reliability of chatter stability of CNC turning operation by MCS method

#### References

• Siddhpura M., Paurobally R. A review of chatter vibration research in turning. International Journal of Machine Tools and Manufacture, Vol. 61, 2012, p. 27-47.
• Schmitz T. L., Smith K. S. Machining Dynamics: Frequency Response to Improved Productivity. Springer, New York, 2008.
• Fu Lianyu, Yu Junyi, Bao Ming Study on the phase characteristics of cutting chatter. Journal of Vibration Engineering, Vol. 13, Issue 4, 2000, p. 510-515, (in Chinese).
• Shi W., Da Q., Xianli L., Fei X. Influence of axial milling force and gyroscopic effect on milling chatter stability. Journal of Vibration, Measurement and Diagnosis, Vol. 6, 2013.
• Huang Q., Zhang G. B., Zhang X. Y., Cao D. Experimental analysis on regenerative chatter model. Journal of Vibration Engineering, Vol. 21, Issue 6, 2008, p. 547-552.
• Shi W., Xianli L., Yanxin W. Chatter prediction based on continuous wavelet features and multi-class spherical support vector machine. Journal of Vibration, Measurement and Diagnosis, Vol. 1, 2012.
• Yang Y.-Q., Liu Q., Wang M. Optimization of multiple tuned mass dampers for chatter suppression in turning. Journal of Vibration Engineering., Vol. 23, Issue 4, 2010, p. 468-474.
• Altintas Y., Stepan G., Merdol D., Dombovari Z. Chatter stability of milling in frequency and discrete time domain. CIRP Journal of Manufacturing Science and Technology, Vol. 1, Issue 1, 2008, p. 35-44.
• Albertelli P., Cau N., Bianchi G., Monno M. The effects of dynamic interaction between machine tool subsystems on cutting process stability. International Journal of Advanced Manufacturing Technology, Vol. 58, Issue 9, 2012, p. 923-932.
• Iglesias A., Munoa J., Ciurana J. Optimisation of face milling operations with structural chatter using a stability model based process planning methodology. International Journal of Advanced Manufacturing Technology, Vol. 70, Issues 1-4, 2014, p. 559-571.
• Choi S.-K., Grandhi R. V., Canfield R. A. Reliability-Based Structural Design. Springer, London, 2006.
• Zhang Y., He X., Liu Q., Wen B. Robust reliability design of banjo flange with arbitrary distribution parameters. Journal of Pressure Vessel Technology, Vol. 127, Issue 4, 2005, p. 408-413.
• O’connor P. D. T., O’Connor P., Kleyner A. Practical Reliability Engineering. John Wiley and Sons, 2012.

26 May 2017
Accepted
27 May 2017
Published
30 June 2017
SUBJECTS
Mechanical vibrations and applications
Keywords
regenerative chatter
CNC lathe
reliability
stability
Acknowledgements

This paper is funded by the international Exchange Program of Harbin Engineering University for Innovation-oriented Talented Cultivation.