We study the dynamics of targeted energy transfers in suppressing chatter instability in a single-degree-of-freedom (SDOF) machine tool system. The nonlinear regenerative (time-delayed) cutting force is a main source of machine tool vibrations (chatter). We introduce an ungrounded nonlinear energy sink (NES) coupled to the tool, by which energy transfers from the tool to the NES and efficient dissipation can be realized during chatter. Studying variations of a transition curve with respect to the NES parameters, we analytically show that the location of the Hopf bifurcation point is influenced only by the NES mass and damping coefficient. We demonstrate that application of a well-designed NES renders the subcritical limit cycle oscillations (LCOs) into supercritical ones, followed by Neimark–Sacker and saddle-node bifurcations, which help to increase the stability margin in machining. Numerical and asymptotic bifurcation analyses are performed and three suppression mechanisms are identified. The asymptotic stability analysis is performed to study the domains of attraction for these suppression mechanisms which exhibit good agreement with the bifurcations sets obtained from the numerical continuation methods. The results will help to design nonlinear energy sinks for passive control of regenerative instabilities in machining.

## Introduction

In machining processes, the undesired vibration of the tool relative to the workpiece gives rise to a low-quality product. One of the most important causes of the instability in the cutting process is the so-called regenerative effect, which arises from the fact that the cutting force exerted on a tool is influenced not only by the current tool position but also by that in the previous revolution. Hence, the equation of motion for the tool is a delay differential equation, which renders even a single-degree-of-freedom (SDOF) dynamical system to be infinite-dimensional. In practical machining processes, regenerative limit cycle oscillations (LCOs) create adverse effects on machining quality, and many papers deal with the stability and bifurcations of machining processes, as well as with various passive and active means to improve their stability (see, for example, Refs. [1–8]).

Direct use or variations of linear/nonlinear-tuned mass dampers (TMDs, [7,8]) are probably the most popular approach to passive chatter suppression. However, even if a TMD is initially designed (tuned) to eliminate resonant response near the eigenfrequency of a primary system, the mitigating performance may become less effective over time due to aging of the system, temperature, or humidity variations, thus, requiring additional adjustment or tuning of parameters. It is only recently that passively controlled spatial/dynamic transfers of vibrational energy in coupled oscillators to a targeted point where the energy eventually localizes were studied by utilizing a nonlinear energy sink (NES, see Ref. [9] for the summary of developments); and this phenomenon is simply called targeted energy transfer (TET). The NES is basically a device that interacts with a primary structure over broad frequency bands; indeed, since the NES possesses essential stiffness nonlinearity (no linear stiffness term), it may engage in (transient) resonance capture [10] with any mode of the primary system. This is why an NES, more effective than a TMD [11,12], can be designed to extract broadband vibration energy from a primary system, engaging in transient resonance with a set of “most energetic” modes [13]. In Particular, Lee et al. [14–16] applied an ungrounded NES to an aero-elastic system, and numerically and experimentally demonstrated that a well-designed NES can even completely eliminate aero-elastic instability. Three suppression mechanisms were identified; that is, recurrent burstouts and suppressions, intermediate and complete elimination of self-excited instability in the aero-elastic system. Such mechanisms were investigated by means of bifurcation analysis and complexification-averaging (CxA) technique [17].

Kalmár-Nagy et al. [3] analytically proved, by means of the center manifold theorem [18], the existence of subcritical Hopf bifurcation in an SDOF machine tool model with the regenerative cutting force. Furthermore, practical stability limit in turning process was investigated by considering contact loss issues in the regenerative cutting force [4], which can predict stable, steady-state periodic tool vibrations (or limit cycle oscillations—LCOs). Nankali et al. [19,20], proposed application of the NES system to a time-delayed machine tool system. They utilized numerical techniques to obtain basin of attraction for different TET mechanisms. Recently, Gourc et al. [21] investigated the different response regimes of a cutting tool on a lathe strongly coupled to a nonlinear energy sink. They derived the equation of the slow invariant manifold (SIM) and explained the behavior of the system by studying the location of the fixed points of the slow flow on this manifold. Moreover, Gourc et al. [22] have studied the passive control of chatter instability in turning processes using a vibro-impact nonlinear energy sink.

In this work, we present a comprehensive study of TET mechanisms in suppressing regenerative chatter instability in a turning process. For this purpose, we first review nonlinear dynamics of a SDOF machine tool model; then, perform a linear stability analysis to explore the effects of NES parameters on the occurrence of Hopf bifurcation (i.e., stability boundary on the plane of cutting depth and rotational speed of a workpiece). Furthermore, the numerical [23] and asymptotic [24] techniques are applied to properly understand the suppression mechanisms that appear by the NES. Moreover, applications of an NES to a practical machine tool model are introduced.

## SDOF Machine Tool Dynamics: Summary

Kalmár-Nagy et al. [3] studied a nonlinear dynamics of a single-degree-of-freedom machine tool model depicted in Fig. 1. The equation of motion can be expressed by a delay-differential equation with nonlinear regenerative terms retained up to cubic order

where $x\tau \u225cx(t\u2212\tau )$ is the delayed variable with a time delay $\tau =2\pi /\Omega $; Ω is the rotating speed of the workpiece; $\omega n=k/m$ is the (linearized) natural frequency; $\omega 0=k1/m$, where *k*_{1} is the cutting force coefficient; *ζ* is the damping factor of the machine tool; and *f*_{0} is the chip thickness at the steady-state cutting.

where the differentiation is now with respect to the new time variable. Note that the parameter $p=(\omega 0/\omega n)2=(k1/k)2$ indicate the effect of stiffness-hardening due to machining conditions.

Kalmár-Nagy et al. [3] calculated the transition curves (cf. Fig. 2(a) for $\zeta =0.1$) as a function parameterized by the eigenfrequency *ω*, $(\Omega (\omega ),p(\omega ))$, through which the steady-state cutting loses stability through a Hopf bifurcation. Since the periodic solution that is born from the equilibrium point is unstable, this bifurcation is a subcritical one. For the damping factor $\zeta =0.1$ (which will be assumed throughout this work), the stiffness ratio *p* has the minimum value, $pmin=2\zeta (1+\zeta )=0.22$, and the minimum eigenfrequency of the limit cycle oscillation (LCO) is $\omega min=1.0$ such that $p(\omega )>0$ for $\omega >\omega min$. Moreover, Fig. 2(b) depicts the minimum and maximum eigenfrequencies (denoted by $\omega min(n)$ and $\omega max(n)$, respectively) that the $nth$ lobe yields. Whereas, the first two or three lobes exhibit relatively broad distribution of eigenfrequencies, the rest possess harmonic components concentrated near $\omega =1.1$. Note that these eigenfrequencies act as “seed” harmonic bases for the new-born periodic motion, and the triggering of a machine tool chatter appears as a result of competition between the eigenfrequency and the rotational speed of the workpiece.

## Nonlinear Energy Sink

and $f=(0,0,\u2212C(x1\u2212x2)3+(3p/10)[(x1\u2212x1\tau )2\u2212(x1\u2212x1\tau )3],\u2212C(x2\u2212x1)3/\u03f5)T$.

**I**is an identity matrix. For a nontrivial eigenvector

**X**, we derive the characteristic equation as $det(\lambda I\u2212A\u2212Re\u2212\lambda \tau )=0$. That is, we write

where *n* is the order of the lobe in the stability chart, and $K(\omega )$ and $\varphi $ are defined in Eqs. (11) and (12).

where $b=1\u22124\zeta 12/\u03f52\u22124\zeta 12/\u03f5$ and $c=\u22124\zeta 12/\u03f52$. It can be analytically shown that $\omega min=1$ in the limit of $\u03f5\u21920$ or $\zeta 1\u21920$.

Figure 3 depicts the minimum eigenfrequency $\omega min$ with respect to the mass ratio *ϵ* and the damping factor *ζ*_{1} of the NES. As *ϵ* and *ζ*_{1} increase, $\omega min$ tends to decrease; however, it remains near unity but decreasing for small mass ratios (e.g., $\u03f5<0.1$) and damping factors (e.g., $\zeta 1<0.1$), which is usually the case in practice. For fixed mass ratios, $\omega min$ decreases as *ζ*_{1} increases; and after a while the eigenfrequency seem to converge a constant (top of Fig. 3(b)). On the other hand, for fixed damping factors, the variation of $\omega min$ with respect to a mass ratio differs from the order of the damping factor (bottom of Fig. 3(b)). That is, for small *ζ*_{1}, there exist a minimum of $\omega min$ and then it increases; but for larger *ζ*_{1} it decreases monotonically. This observation suggests not only that the frequency of the subcritical LCO when no NES is applied should be higher than unity, but also that the LCO with a frequency altered by an NES can be smaller than unity. Note that the linearized natural frequency in the rescaled nondimensional equation (4) is unity.

Now, we compute the transition curves based on Eqs. (10) and (13). Figure 4 depicts the changes of the stability boundary in (Ω, *p*)-plane by varying the mass ratio *ϵ* and fixing the other two NES parameters. The stability enhancement due to the application of an NES can be measured by directly calculating the point-wise shift amount as $\Delta p=(p\u2032\u2212p)/p\xd7100$ (%), where *p* and $p\u2032$ denote the values at the stability boundary with respect to each Ω without and with an NES, respectively. Upward shift of the stability boundary occurs more significantly near the valley than near the cusp points of the lobes, which will be useful in practical applications of chatter suppression. The shifting amount of the transition curve does not appear to be significant with a small NES mass (about 5% improvement near valley of the lobes); however, the upward shift becomes increasing monotonically as the mass ratio increases. The ranges of the eigenfrequencies at the transition curves tend to become lower as the mass ratio increases; and above certain mass ratio the eigenfrequency intervals are shifted upward (cf. Fig. 4(b) when $\u03f5=0.6$). Note that the bifurcation occurring on the transition curve can be referred to as a degenerate Hopf bifurcation in delay-differential equations, because one of the eigenvalues is always zero.

We remark that, although we can delay the occurrence of Hopf bifurcations by adding NES, this is not all one can achieve with an NES as a passive broadband vibration controller. The application of an NES can also alter the topology of local bifurcations such that it produces Neimark–Sacker and saddle-node bifurcations as well as Hopf bifurcation. The former two bifurcations are essential in discussing targeted energy transfer (TET) mechanisms in suppressing any types of instabilities introduced in a dynamical system. We deal with this in more detail in Sec. 4 by utilizing a numerical continuation technique for a system of delay-differential equations.

## Bifurcation Analysis and Robustness

### TET Mechanisms.

As in the previous aero-elastic applications [14], three distinct TET mechanisms are identified in suppressing regenerative chatter instability; that is, recurrent burstouts and suppression, intermediate and complete elimination of regenerative instability (cf. see Fig. 5 for typical time history for each suppression mechanism).

*The first suppression mechanism* is characterized by a recurrent series of suppressed burstouts of the tool response, followed by eventual complete suppression of the regenerative instabilities. The beating-like (quasiperiodic) modal interactions observed during the recurrent burstouts turn out to be associated with Neimark–Sacker bifurcations of a periodic solution (cf. Fig. 6) and is critical for determining domains of robust suppression [16]. To investigate this mechanism in more detail, Fig. 7 depicts the displacements of both the tool and NES and their wavelet transforms. Energy exchanges between the two modes are evidenced in Fig. 7, through which a series of 1:1 transient resonance captures and escapes from resonance occurs (see Fig. 9(b).

*The second suppression mechanism* is characterized by intermediate suppression of LCOs, and is commonly observed when partial LCO suppression occurs. The initial action of the NES is the same as in the first suppression mechanism. Targeted energy transfer to the NES then follows under conditions of 1:1 transient resonance capture, followed by conditions of 1:1 permanent resonance capture where the tool mode attains constant (but nonzero) steady-state amplitudes. We note that, in contrast to the first suppression mechanism, the action of the NES is nonrecurring in this case, as it acts at the early phase of the motion stabilizing the tool and suppressing the LCO.

Finally, in the *third suppression mechanism*, energy transfers from the tool to the NES are caused by nonlinear modal interactions during 1:1 RCs. Both tool mode and the NES exhibit exponentially decaying responses resulting in complete elimination of LCOs.

### A Practical Tool Model With Contact Loss.

We note that the numerical and analytical studies for TET mechanisms above are valid only for vibrations with small amplitudes; in particular, the permanent contact model with truncated nonlinear terms cannot predict any stable steady-state periodic vibrations of high amplitudes. That is, the truncated nonlinearity in the regenerative cutting force will not predict the existence of a saddle-node bifurcation point right after contact loss occurs. The details of machine tool dynamics can be found in Kalmár-Nagy [4], where stable periodic motions are predicted. Performing numerical continuation analysis, we obtain the LCO surfaces for without and with NES being applied (Fig. 8). Still three distinct TET mechanisms for the model with contact loss are observed.

### Complexification-Averaging (CxA) Technique.

where $\varphi \u2261\beta 1\u2212\beta 2$.

*β*

_{2}to yield

Figure 9(a) directly compares the approximate (Eq. (19)) and exact (Eq. (4)) solutions for the tool displacement, which demonstrates a good agreement; furthermore, the nontime-like patterns (i.e., spirals) of the phase difference $\varphi $ depicts that the underlying TET mechanism for the first suppression mechanism involves a series of 1:1 transient resonance captures and escapes from resonance. Also, Fig. 9(b) depicts 1:1 resonance capture in the slow-flow phase plane ($\beta ,\beta \u0307$) in which $\beta =\beta 1\u2212\beta 2$.

## Analytical Study: Basins of Attraction for Suppression Mechanisms

We note that through this approximation only external nonlinearities due to the NES are left in the system. Further, the NES stiffness *C* is dropped out from the equations which indicates indecency of the suppression mechanisms from this parameter.

*v*and

*w*are the physical quantities for the center of mass (with a factor of $1+\epsilon $) and the relative displacement, respectively. Then, the equations of motion (22) become

*e*, we obtain the slow-flow equation

^{iwt}*R*and

*F*represent real amplitude modulations while $\delta =\delta =\delta 1\u2212\delta 2$ is the phase difference. Equation (29) represents the slow dynamics of the original system (Eq. (21)). The slow dynamics is analyzed using the method of multiple scales. Based on this method, we consider different scales of time by defining $t=\epsilon n\tau n$, where $n=0,1,2,\u2026$ and consequently, $(\u2202/\u2202t)=(\u2202/\u2202\tau 0)+\epsilon (\u2202/\u2202\tau 1)+\cdots $. Moreover, we introduce the perturbation-series of the variables as

*ε*, we derive subproblems governing the solution of slow dynamics. In this study, we calculate the responses for the “slow” time scale (

*τ*

_{0}) and “super-slow” time scale (

*τ*

_{1}). The first-order approximation is computed by considering $\epsilon 0$ subproblem

*R*

_{0}is fixed with respect to slow time variable

*τ*

_{0}, but not with respect to slower time variable

*τ*

_{1}. No limit cycle oscillations are possible for system (31) and the only steady-state solutions are in the form of equilibrium points. Equilibrium points of the slow dynamics with respect to slow time (

*τ*

_{0}) can be calculated as

*ζ*and

*ω*, there are either one or three branches, when plotting $R\u030202$ versus $F\u030202$. For the case of three branches (Figs. 10–12), performing linearized stability analysis on the slow dynamics with respect to slow time, reveals that the middle branch is unstable while the other two are stable. Equilibrium points of stable and unstable branches are in the form of node (stable) and saddle node, respectively. Therefore, at the leading-order approximation the dynamics will be attracted to either of stable nodes.

*ε*should be considered to compute possible LCO's for the slow flow. To this end, we plug Eq. (30) into Eq. (31) and match coefficients of $\epsilon 1$ for the first equation to get

*τ*

_{1}). This is called super-slow flow (SSF). The intersection points of the SIM and SSF can be calculated by plugging $(\u2202R\u03020/\u2202\tau 1)$ from Eq. (32) into Eq. (34) and solving for $(\u2202F\u030220/\u2202\tau 1)\u2261f(F\u030202)$

*First suppression mechanism:* This mechanism corresponds to the existence of stable LCO on slow dynamics. It refers to the condition in which stable equilibrium point of slow dynamics $(F\u03020e22)$ lies on the unstable (middle) branch of SIM, while $F\u03020e12$ is unstable. In fact, the attractor is in the form of LCO rather than equilibrium point for the slow dynamics. Cycle of relaxation oscillation of slow dynamics creates a quasi-periodic oscillation in full-order system which are illustrated in Fig. 10. We remark that none of the equilibrium points of slow dynamics are complex numbers, and the trivial equilibrium point is unstable for this mechanism.

*Second suppression mechanism:* This mechanism corresponds to stable LCO for full dynamics. According to the slow, fast dynamics separation (Eq. (25)), second suppression mechanism refers to existence of stable nontrivial equilibrium point of slow dynamics. As derived in Eq. (35), there are two nontrivial equilibrium points $(F\u03020e22,F\u03020e32)$ for the slow dynamics. It can be shown that $F\u03020e22$ is stable when it lies on the stable branches of SIM which indicates second suppression mechanism. Note that none of the slow dynamics equilibrium points are complex, and only the second equilibrium point is stable for this mechanism. The SIM and super-slow flow corresponding to this mechanism are shown in Fig. 11, along with the time simulation of the full dynamics.

*Third suppression mechanism:* This mechanism corresponds to the stable trivial equilibrium point of slow dynamics. It can be shown that if $f\u2032(F\u030220e1)<0$, then origin is the only attractor of the slow dynamics and third suppression mechanism is predicted (Fig. 12).

*No suppression:* In addition to the three suppression mechanisms explained above, there is possibility of no suppression. This happens when trivial equilibrium point is the only real equilibrium point of the slow dynamics; e.g., other two equilibrium points are complex. In this case, if trivial equilibrium point is unstable $(f\u2032(F\u0302e02)>0)$, no attractor exists for the slow dynamics which predicts instability of the original system. Super-slow flow and SIM for this case are shown in Fig. 13.

Table 1 summarizes mathematical conditions for the SIM and SSF configuration, corresponding to different cases discussed above. Also, basin of attraction for all suppression mechanisms, in $\zeta 1\u2212p$ plane for given set of parameters, is depicted in Fig. 14.

Sup. mechanism type | Mathematical condition | Figures |
---|---|---|

First | $F20n<F\u030220e2<F20p\u2009and\u2009F20u<F\u030220e3$ | 10 |

Second | $0\u2264F\u030220e2\u2264F20n\u2009or\u2009F20p\u2264F\u030220e2$ | 11 |

Third | $f\u2032(0)<0\u2009and\u2009Im(F\u030220e2,3)\u22600$ | 12 |

No suppression | $f\u2032(0)>0\u2009and\u2009Im(F\u030220e2,3)\u22600$ | 13 |

Sup. mechanism type | Mathematical condition | Figures |
---|---|---|

First | $F20n<F\u030220e2<F20p\u2009and\u2009F20u<F\u030220e3$ | 10 |

Second | $0\u2264F\u030220e2\u2264F20n\u2009or\u2009F20p\u2264F\u030220e2$ | 11 |

Third | $f\u2032(0)<0\u2009and\u2009Im(F\u030220e2,3)\u22600$ | 12 |

No suppression | $f\u2032(0)>0\u2009and\u2009Im(F\u030220e2,3)\u22600$ | 13 |

Figure 15 compares the numerical and analytical calculated basin of attractions for TET mechanisms of SDOF machine tool. Dashed lines depict boundary of suppression mechanism computed numerically utilizing DDEBIFTOOL, while solid lines correspond to asymptotic analysis results. Numerical boundaries introduce Neimark–Sacker, saddle-node, and Hopf points for the first-, second-, and third-suppression mechanism, respectively. It is obvious that the numerical and analytical results are well matched. We note at this point that the accuracy of asymptotic stability analysis dependents on the value of mass ratio (*ε*). The smaller the mass ratio is the more accurate basin of attraction is computed. That is because we consider *ε* as the perturbation parameter in our study.

As it is plotted in Fig. 15, there is a type of bifurcation we could not detect through numerical study. It is not observed from analytical study for large values of *ε* either. This bifurcation which is called Shilnikov homoclinic bifurcation is predicted when strongly modulated response (SMR) disappears due to the existence of an unstable equilibrium point in an SMR cycle. In other words, it occurs when an unstable equilibrium point $(F\u030220e3)$ meets a point of an SMR cycle $(F\u030220u)$ while the second equilibrium point $(F\u030220e2)$ is unstable. Figure 16 illustrates this situation.

Domains of attraction obtained from asymptotic analysis and associated points on the bifurcation surface are plotted in Fig. 17. A similar technique can be used to compute basins of attraction for suppression mechanisms in the $\epsilon \u2212p$ plane. Figure 18 illustrates the bifurcation surface computed numerically (DDEBIFTOOL), and the corresponding analytical basin of attraction in the $\epsilon \u2212p$ plane. As seen, for small values of *ε* there is a good agreement between the numerical and analytical results. However, larger values of *ε* decrease the accuracy of the asymptotic analysis. That is because *ε* is used as perturbation parameter and has to be small ($\epsilon \u226a1$).

## Concluding Remarks

In this paper, we studied targeted energy transfer phenomena in suppressing chatter instability in a single-degree-of-freedom machine tool system, to which an ungrounded nonlinear energy sink is connected. Two models were considered for the tool dynamics: permanent contact model and contact loss model. The limit cycle oscillation due to the regenerative instability in a tool model which appeared as being subcritical for permanent contact model were (locally) eliminated or attenuated at a fixed rotational speed of a workpiece (i.e., a delay period) by TETs to the NES. It was shown that there should be an optimal value of damping for a fixed mass ratio to shift the stability boundary for stably cutting more material off by increasing chip thickness. Also, magnitude of NES nonlinear stiffness does not have any effect on stability boundary while increasing mass ratio improves stability. Three suppression mechanisms have been identified and each mechanism was investigated numerically by time histories of displacements, and wavelet transforms and instantaneous modal energy exchanges. Furthermore, we extend the CxA analysis to perform asymptotic analysis by introducing a reduced-order model and partitioning slow-fast dynamics. The resulting singular perturbation analysis yields parameter conditions and regions for the three suppression mechanisms, which exhibit good agreement with the bifurcations sets obtained from numerical continuation methods.

## Acknowledgment

This work was supported in part by National Science Foundation of United States, Grant Numbers CMMI-0928062 (YL) and CMMI-0846783 (TKN).