Nonlinear dynamics of a regenerative cutting process

We examine the regenerative cutting process by using a single degree of freedom nonsmooth model with a friction component and a time delay term. Instead of the standard Lyapunov exponent calculations, we propose a statistical 0-1 test analysis for chaos detection. This approach reveals the nature of the cutting process signaling regular or chaotic dynamics. For the investigated deterministic model, we are able to show a transition from chaotic to regular motion with increasing cutting speed. For two values of time delay showing the different response, the results have been conﬁrmed by the means of the spectral density and the multiscaled entropy.


Introduction
A cutting process is a basic machining technology to obtain the surface of the assumed parameters. In certain working conditions, it can be disturbed by chatter appearing as unexpected waves on the machined surface of a workpiece. The appearance of chatter was noticed and described by Taylor in the beginning of the twentieth century [1]. But the first approaches toward explanations of this phenomenon came about 50 years later through the analysis of self-sustained vibrations [2], regenerative effects [3], structural dynamics [4,5], and finally the dry friction phenomenon [6,7]. Consequently, elimination and stabilization of the associated oscillations have become of high interest in science and technology [8][9][10]. The plausible adaptive control concept, based on relatively short time series [11], has been studied to gain deeper understanding.
Recently, apart from the widely developed chatter vibrations chaotic oscillations caused by various system nonlinearities were predicted and detected [12][13][14][15][16][17][18][19][20]. The recent technological demand is to improve the final surface properties of the workpiece and to minimize the production time with higher cutting speeds [21]. Thus, a better understanding of the physical phenomena associated with a cutting process becomes necessary [22]. In this paper, we will continue the work on chaotic instabilities in cutting processes proposing the 0-1 test [23][24][25][26][27] as a tool identifying a possible chaotic solution [28]. This paper is organized as follows. After the present introduction (Sect. 1), we describe the model in Sect. 2. In Sect. 3, we provide the results of the simulations and corresponding power spectral densities (PSD) while in Sect. 4 the 0-1 test is applied, and subsequently the findings are confirmed by means of the multiscale entropy (Sect. 5). The paper ends with conclusions (Sect. 6).

The model
A regenerative cutting process may exhibit a wide range of complex behavior due to frictional effects [9,12], structural nonlinearities [29], and delay dynamics [17,18,[30][31][32]. Moreover, it may also involve loss of contact between the tool and the workpiece. The following equations model the regenerative cutting process and the mentioned properties.
After the first pass of the tool, the cutting depth can be expressed as where y(t − τ ) corresponds to the position of the workpiece during the previous pass, and τ is the time delay scaled by the period of revolution of the workpiece 2π/Ω 0 (Fig. 1). The motion of the workpiece can be determined from the model proposed by Stépán [30] y + 2γẏ + ω 2 where ω 0 = √ k/m is the frequency of free vibration, v 0 is the feed velocity, and 2γ = c/m is the damping coefficient. F y (h) is the thrust force, which is the Fig. 1 Physical model of a regenerative cutting process [17] horizontal component of the cutting force, and m is the effective mass of the workpiece. The thrust force F y is based on dry friction between the tool and the chip. It is assumed to have a power law dependence on the actual cutting depth h and to be proportional to the chip width w and a friction coefficient c 1 . Θ(·) denotes the Heaviside step function. The restitution parameter β = 0.75 is associated with the impact after contact loss, while t − and t + denotes the time instants before and after the impact. Substituting (1) into (2), we derive a delay differential equation (DDE) for the workpiece motion y(t). Plugging its solution into (1) results in the history of cutting depth h(t).

Simulation results
The nonsmooth model equations are solved by a simple Euler integration scheme. The used parameters [17,32] are presented in Table 1. Furthermore, the feed velocity v 0 has been assumed to be fairly large so that v 0 >ẏ. Note that, in this case, the system nonlinearities are limited to the exponential dependence of the cutting force on the chip thickness and to the contact loss between the tool and workpiece.
The corresponding time series for two choices of the time delay parameter τ = 1.8 and 2.1 ms are presented in Fig. 2. These series have been plotted with points. On the first sight, one can notice that both solutions are complex but Fig. 2a shows points grouped in selected lines while the distribution of time history points of Fig. 2b looks more random. In Fig. 2b, h reaches negative values that signal that the contact between the tool and the workpiece is lost. The power spectral densities (PSD) of cutting depth S(ω) = 2π/T |F {h(t)}| 2 for the two chosen delay times (τ = 1.8 ms and τ = 2.1 ms) 1 indicate a   A broad band spectral density indicates chaotic/stochastic dynamics whereas sharp peaks imply regular motion transition from regular to chaotic motion. The sharp peaks in Fig. 3 belong to a high-periodic orbit (regular motion) whereas the broad spectrum indicates chaotic dynamics.
Both power spectra are dominated by a main peak. In case of regular motion, its position belongs to the delay time τ = 1.8 ms while in case of chaotic dynamics the time scale belonging to the peak (t p ≈ 2.0 ms) is smaller than the delay time τ = 2.1 ms. This smaller value could be a consequence of a tool-workpiece contact loss. Based on that, we take a closer look on other measures to characterize the model's dynamics and use a 0-1 test for chaos to display a possible transition from regular to chaotic motion with increasing delay time τ .

Application of 0-1 test
Based on the time series {h j } which is a discretization of the solution h(t) of the DDE normalized by its standard deviation, we define dimensionless displacements in the (p, q)-plane in the following way [23,24,28]: where c 0 is a constant. In this way, regular dynamics is related to a bounded motion while any chaotic dynamics leads to an unbounded motion in the (p, q)-plane [23]; see Fig. 4a.
To obtain a quantitative description of the examined system, we perform calculations of the asymptotic properties defined by the total mean square displacement (MSD) M(n), Fig. 4b, and finally we obtain the growth rate K in the limit of large times For almost all values of the constant c 0 , the parameter K is approaching asymptotically 0 or 1 for regular or chaotic motion, respectively. Note, practically, one has to truncate the sums in (4). Thus, we derived K ≈ 0.21 for τ = 1.8 ms and K ≈ 1.09 for τ = 2.1 ms, which supports the first impression gained from the time series themselves, Fig. 2a and b. Note further that for delay time τ = 1.8 ms, K decays with increasing n on much smaller values, Fig. 4c, which corroborates the result pointing toward regular motion.
Note that the parameter c 0 acts like a frequency in a spectral calculation, cp. (3). If it is badly chosen, c 0 / t resonates with one frequency of the process dynamicsh(t). Such a frequency belongs to a peak in the PSD, Fig. 3. In the 0-1 test, regular motion would yield a ballistic behavior in the (p, q)-plane and the corresponding quadratic growth of MSD results in an asymptotic growth rate. The disadvantage of the test, its strong dependence on the chosen parameter c 0 , could be overcome by a proposed modification. Gottwald and Melbourne [24,26] suggest to take several randomly chosen values of c 0 and compute the median of the belonging K-values. Particularly, in [26], the problems of averaging over c 0 as well as sampling the data points are discussed extensively. We followed this approach [26,33], which improves the convergence of the test (Fig. 4c) without the consideration of longer time series, to find the time delay τ leading to chaos (see Fig. 5). We defined a modified square displacement D(n) which exhibits the same asymptotic growth where the oscillatory term V osc (n, c 0 ) can be expressed by and E[h] denotes the average of examined time series h i where N max is the number ofh i elements. Consequently, the oscillatory behavior is subtracted from the MSD M(n, c 0 ) and the regression analysis of the linear growth of D(n, c 0 ) (see (6)) with increasing n is performed using the linear correlation coefficient which determines the value of K c 0 .
where vectors X = [1, Finally, the median is taken of K c 0 -values (see (9)) corresponding to 100 different values of c 0 ∈ (0, π). The results of K for different delay times τ , Fig. 5, in the window between 1.75 ms and 2.3 ms indicates a transition from regular to chaotic dynamics with increasing delay time in the region of 1.9 ms.
As a consequence, we conclude that in the investigated window increasing cutting speed leads to a transition from chaotic chatter dynamics to regular motion with improved surface quality.

Multiscale entropy
To characterize the solutions of the DDE, (1) and (2), with regard to information production rate and complexity, we aim to calculate multiscale entropy (MSE) [34]. This method was successfully applied to analyze the complexity of biological signals [34,35]. It is suitable for short and noisy time series. As a consequence, the chosen procedure would be applicable to experimental data as well. We use an algorithm provided by PhysioNet [36]. First, we compute coarsegrained time series {x (N ) } using nonoverlapping intervals containing N equidistant data points h i , In the next step, we calculate sample entropy S (N ) E [37] for these coarse-grained time series. Sample entropy is the negative of the logarithm of the conditional probability that sequences of m consecutive data points x i+m−1 ) and x (N ) j close to each other will also be close to each other when one more point is added to them. Hence, it is estimated as follows: where U (N ) m (r) represents the relative frequency that a vector x (N ) i is close to a vector x (N ) j (i = j ). Close to each other in the sense that their infinity norm distance is less than ε = rσ . By σ , we denote the standard deviation of the data. In the limit of m → ∞ and r → 0, sample entropy is equivalent to order-2 Rényi entropy K 2 and is suitable to characterize the system's dynamics [38]. For independent variables {ξ }, the entropy follows from S To clear up the characteristics of the cutting process, we look at MSE depending on box size r for the two chosen delay times, Fig. 6. For regular motion, we expect the entropy to approach zero with decreasing r. This is observed for the time series with delay time τ = 1.8 ms. For chaotic dynamics, the entropy should stay finite, observed for τ = 2.1 ms. For the sake of completeness, it should be mentioned that in the case of stochastic dynamics the entropy would diverge with decreasing spatial resolution r, cp. (13). In Figs. 6 and 7, we further analyze the scale factor dependence of MSE. The entropic measure is always larger for the chaotic time series since it is the more complex one. MSE for small scale factor, Fig. 7a, indicates that there is no characteristic time scale, comparable to 1/f -noise [34]. But for larger scale factors, Fig. 6 Multiscale entropy S E depending on scale factor N and box size r. For regular motion, we expect the entropy to approach zero with decreasing r. For chaotic dynamics, the entropy should stay finite. Since S E is not decreasing with scale factor significantly, it seems there is no characteristic time scale present in the data MSE is decaying comparable to Gaussian white noise, Fig. 7b. Thus, even in the chaotic case, there exists a characteristic time scale which is close to the delay time. The frequencies dominating S(ω), Fig. 3, are also present in S (N ) E (m, r). They belong to minima in Fig. 7. We learn coarse-graining of the data over multiple of time scale t p belonging to structures in the cutting process dynamics leads to less complex time series and contains less information.

Conclusions and last remarks
Concluding, the 0-1 test differentiates between the two types of motion. Depending on the chosen delay time for the investigated DDE, (1) and (2), regular or chaotic motion is observed and a transition from chaotic to regular motion is detected with increasing cutting speed. The nature of solutions has been also confirmed by the corresponding power spectral densities and multiscale entropies. The latter reveals more insights into the process dynamics but is of much higher computational cost than the 0-1 test and the spectral calculations.
The 0-1 test appeared to be relatively simple, and consequently, useful for systems with delay and discontinuities. A huge advantage of the test is its low computational effort and the possibility to compute it "on the fly" while the data is still growing. One of the useful aspects of the 0-1 test is that the result can be plotted against the parameter τ .
The presented method gives a quantitative criterion for chaos similar to the maximum Lyapunov exponent.
As demonstrated by Falconer et al. [25] and Krese and Govekar [33], the method can be used on experimental data as well. Unfortunately, in case of the cutting process, experimental data are often characterized by a relatively high level of noise [20]. In the examined system, we waived the possibility of additive noise. It was shown that the 0-1 test could be applied on dynamical systems with additive noise and a good signal to noise ratio [24].