Comparative analysis on recent deep learning techniques for identifying chatter in milling process

Abstract

Milling seems to be the most extensively utilized production technology in modern manufacturing industries, and it plays a significant role. Chatter is a type of disturbance in the form of vibration that has a negative impact on machining operation. Chatter recognition utilizing sensor outputs is a hot topic in academia. Although some progress has indeed been documented utilizing various featurization techniques and ml techniques, conventional approaches have a number of limitations, including manual preparation and a huge dataset need. Although, these are widely being used to evaluate milling operations in terms of production efficiency & work piece surface quality,.they are not suited for real applications due to their computing duration and require large data for training process. Therefore, in this study, three well-performing deep learning approaches such as LSTM, DTW, and Bi-LSTM are used to provide an effective way for monitoring and managing chatter in the milling processes with the Duplex 2205 material. Here, some of the parameters like acceleration is measured while the milling operation is taking place, and the measured acceleration value is processed using selected three DL techniques for identifying the presence of chatter and are tested to see which one performs the best. The Bi-LSTM outperformed other approaches in detecting chatter present, according to the data.

Keywords

Bi-directional long short term memory long short term memory dynamic time warping deep learning acceleration milling chatter detection

ANN

Artificial Neural Networks

HMM

Hidden Markov Models

Fuzzy Logic

SVM

And Support Vector Machines

STFT

Short-Time Fourier Transform

DTW

Dynamic Time Wraping

SSV

Variable Speed Machining

kNN

K-Nearest Neighbor

FFT

Fast Fourier Transform

WPT

Wavelet Packet Transform

AUC

Area Under Curve

TDA

Topological Data Analysis

National Power

Genetic Algorithms

CNN

Convolutional Neural Networks

Radial Basis Function

RBF

Radial Basis Function

MLP

Multi-Layer Perceptrons

Fault Detection

EEMD

Ensemble Empirical Mode Decomposition

MRR

Material Removal Rate

RMSE

Root Mean Square Error

Deep Leraning

Machine Learning

SLD

Stability Lobe Diagram

Bi-LSTM

Bi-Directional Long Short Term Memory

SAE

Stacked Auto Encoder

LSTM

Long Short-Term Memory

PSO

Particle Swarm Optimization

SSDAE

Stacking Sparse De-noising Auto-Encoder

LSSVM

Least Squares Support Vector Machine

BPNN

Back Propagation Neural Network

PLSR

Partial Least Squares Regression

ELM

Extreme Learning Machine

DCNN

Deep Convolutional Neural Network

CWT

Continuous Wavelet Transform

LAMM

Laser-Assisted Micro-Milling

MQM

Machining Quality Monitoring

MRP

Multi-layer Recurrence Plot

1 Introduction

One of the best common cutting procedures for machining surfaces is end-milling. Milling can be utilized to create both flat surfaces and also varied profiles. Chatter has a significant influence on the pace of material removal as well as the pace of surface creation [1]. The word “chatter” refers to the cutting tool’s self-excited intense relative dynamic movement with the work piece. Chatter has a negative impact on output quality, milling precision, tool wear rate, bearings and life of the tool. The renewal of waviness on the ‘s surface causes chatter [2, 3]. This is referred to be regenerating chatter, and certainly it is one of the most common reasons of cutting procedure instabilities. To prevent this instabilities, an industrial equipment developer should first examine the spindle unit’s rigidity and damping, as well as the tool. To avoid noise, the machine tool’s user must understand how to select cutting settings. To achieve more surface smoothness, the correct cutting parameters must be set up before the process begins [4]. A smooth (i.e. no chatter) and a highly unstable (i.e. chatter) snip may be distinguished using the axial depth-of-cut as a factor of spindle speed [5]. Cutting forces, vibrations, sound and acoustic emission, and other critical aspects are retrieved from physical data in order to monitor the machining state [6, 7]. The extracted characteristics are crucial in determining the machining condition.

The lack of dynamic stiffness inside one or even more components of a system, such as the tool holder, work piece, machine tool, or cutting tool, causes self-excited vibration, forced vibration, as well as free vibration in metal cutting / machining operations [8]. Real effect is the source of free vibration. Asymmetry within spindles, bearings, gears, and many other available equipment causes enforced vibration. Once the perpetrators have indeed been identified, preventing, reducing, or eliminating these 2 kinds of vibration is relatively simple [9]. The alternative force generated by the link between the machining material and indeed the instrument causes self-excited vibration, which is an unattenuated vibration. A machine becomes more volatile, unstable and unmanageable as a result of this type of vibration. Mode coupling chatter regenerative chatter as well as frictional chatter are the 3 types of chatter dependent on the sort of self-excitation mechanism involved in the vibration [10]. Mutual friction between the work piece and then the instrument on the very same side as the cutting force and speed causes frictional chatter. The pairing of two naturally formed vibration modes causes mode coupling chatter [11, 12], which is caused by a small difference in the hardness of vibrating components in 2 different directions. Using experimental modelling and analysis, a stable state for the system may be derived [13]. Regenerative chatter is instigated by variations in the heaviness, thickness and density of the chip caused by a phase incompatibility among the vibration pattern created throughout a previous cutting operation as well as the vibration displacement of a subsequent cutting procedure [14].

Cutting state detection requires a pattern classifier in addition to feature extraction in a chatter recognition system. ANN, FL, HMM, and SVM are only a few of the clever algorithms that have been proposed thus far [15]. Multiple local minima and over-fitting are common problems with ANNs, although SVMs may overcome these flaws [16]. Furthermore, SVMs have a straightforward geometric interpretation and provide a sparse solution. As a result, SVMs are widely regarded as a novel effective way for generating efficient and powerful classification algorithms [17]. Most of the above-mentioned chatter detection algorithms can only identify chatter that is virtually fully grown. The commencement of chatter must be detected at such an early stage that no chatter markings appear on the work piece [18]. Additionally, a suitable chatter suppression approach, such as a time-varying method based on magneto-rheological fluid, can be used to suppress chatter over time [19]. However, chatter detection remains a significant difficulty at this time, requiring a quick detection technique that includes feature extraction and smart classification [20]. In several industries, including finance, aviation-aerospace, healthcare and automotive, the Deep Learning technique, which is a branch of ML, has been extensively utilized & demonstrated to be more successful than all the other conventional computing intelligence approaches. Deep learning techniques, that have led to ground-breaking new devices and systems in these fields, are now acquiring popularity in the domain of milling, where they are being used to make milling machines extremely smart and increase productivity levels.

2 Literature survey

2.1 Related works

Yun et al. [21] The image assessment of dominating bandwidths from STFT spectrograms is used in this study to present an unique feature extraction technique for chatter identification. By implementing the squared energy operators to the synthesised FFT spectrum, dominating frequency bands with good intensity are discovered in order to reduce external disturbances and emphasise chatter relevant characteristics features. Dry milling experiments of titanium alloy Ti6Al4 V are used to validate the suggested feature extraction approach, which is then contrasted to 2 other feature extraction methods. The findings demonstrate that time-frequency image characteristics from dominant bandwidths are effective for chatter identification and outperform time domain and wavelet-relied features in aspects of distinction.

Albertelli et al. [22]. In the current study, a chatter recognition algorithms was created that is especially designed for changeable speed milling. In further context, the cutting stability evaluation is done out somewhere in the spindle angular domain using real-time estimation of a normalised chatter indication based on the cyclo-stationary theory. The vibrational signals as well as information from the spindle encoder are pre-processed using order monitoring and synchronous averaging approaches prior generating the chatter indicator. Real-world machining processes involving both constant and SSV were used to validate the developed chatter monitoring methods. In almost all of the cutting settings evaluated, the designed algorithm was shown to be suitable of recognising chatter quickly and reliably.

Zhang et al. [23]. The construction of a reliable controller for chatter reduction is the topic of this study. For durable controller architecture, the machining dynamic equations involving active control force were first created and rationalized using the approximation linear time invariant model. Then, modal parameter perturbation modelling and milling process parameter perturbation modelling were created as well as studied. The control strategy was created depending on the perturbation models. The machining process numerical simulation evaluates the efficacy of the specified controller. The flank milling & end milling experiments were carried out in order to confirm the real impact. The testing findings reveal that the suggested algorithm effectively suppresses noise and leads to a finer surface of the work piece, demonstrating the device’s efficacy and reliability.

Yesilli et al. [24]. We provide an alternate method for chatter recognition depending on the kNN algorithm for categorization and DTW as a time series measure of similarity in this work. The acceleration signals output obtained from the equipment bearer in a sequence of turning trials were taken as the time series. The findings reveal that this strategy delivers detecting accuracy rate that beat conventional approaches in the vast majority of instances. The results were contrasted with WPT and EEMD focused approaches, or even the more current TDA relied method. We showed that our DTW-based technique has the greatest average categorization rate in 3 out of 4 cutting configurations, with one example achieving 99 percent precision. Our method really doesn’t need feature extraction, can reuse a classifier over several cutting configurations, & utilizes training sets of reasonable quantity. Despite our approach’s great precision is accompanied with a significant computing expense, this is unique to the DTW implementation we utilized. As a result, additional code optimization and a lower computing burden throughout the execution stage makes this technique a feasible choice for in-process chatter identification.

Torabi et al. [25] have conducted a comprehensive literature review of current milling process modelling methodologies, with a focus on recent breakthroughs and uses of AI-based modelling tools. A comparison of the available methods will be offered, as well as the applicability of each approach for the various sorts of investigations. In addition, the flaws of each technique will be discussed, as well as open research problems. As a result, a thorough examination of current advancements in the area will serve as a reference for selecting the most appropriate modelling approach for this process’s aims, circumstances, and requirements.

Zhang et al. [26] based on normal operating data obtained onsite at a NP power station in the UK, have proposed an unique coal mill modelling approach employing GA. Modelling of an E-type vertical spindle coal mill is the topic of this research. The model performance of two distinct mills is assessed under a variety of operating situations. The simulation results reveal that the model responses and observed data are in good agreement. Without resorting to extensive mill testing, the relevant data may be gathered, and the model can be developed without difficulty in computing. As a result, the work is applicable to a wide range of situations.

Serın et al. [27] CNN, a deep learning framework, were used to identify chatter in a slot machining process of the research [31]. An AE sensor were installed to the milling machinery in this study, and it gathers acoustic emissions. Chatter is a type of vibration / disturbance that has a negative impact on milling productivity. Chatter’s unwelcome vibrations can result in poor surface quality and sometimes even tool failure. Because chatter reduces the efficiency of milling processes, it is critical to identify and respond to chatter as soon as possible. In addition to identifying chatter, acquired info was utilized to test and evaluate the estimating model.

Lamraoui et al. [28]. This research provides a new approach to detecting noise in CNC milling machines. This approach is based on artificial intelligence and vibratory signal analysis. There are five important phases in the methodology: 1) data collection, 2) signal processing, 3) feature creation, 4) feature selection, and 5) classification are the steps in the process. A multiband resonance filtering approach is presented at the pre-processing stage because chatter components emerge around system resonance frequencies. After that, envelope analysis is performed. As a result, the signal-to-noise ratio and sensitivity of created characteristics can be improved. The RBF and MLP neural network techniques are compared. The created method is used to identify noise in the Huron K2X10 milling machine. This method is put to the test on a milling machine at varying depths of cut and rotational speeds. The outcomes of the discussions reflect the efficiency of the suggested approach.

Guo [29] has talked about the three steps of developing an online crown/shape control model. Each step requires special attention and offers mill operators with a variety of data. Understanding mill rolling behaviour is made easier by the theoretical derivation of the first step. The linear crown control system is developed as a result of the use of the crown/shape concept. The second step of development is to create an offline simulation model that connects theory and practise. It’s utilised to check the theoretical model, calculate control gain factors, and even test the control algorithm’s stability. When doing mill setup calculations, the real-time process model comes in last. It must deal with the uncertainties of measured devices, probable theoretical model flaws, and real-world mill operating circumstances. Its primary task is to ensure that the mill is set up in a steady, optimum, and precise manner. Statistical approaches are applied throughout these models, from tuning through process optimization. Since 1992, this best-in-class linear crown/shape control system has been effectively used in a production hot strip mill. The creation and use of this rolling model are described in this article.

Li and Ding [30]. The study looked into the best design methodologies for detecting cumulative defects in unpredictable industrial platforms. To be more specific, the design techniques of observer-based FD systems are investigated first in a closed-loop system configuration, with the objective of boosting fault detectability and system resilience versus model instabilities at the same time. After that, an optimum design strategy for observer-based FD systems is investigated for open-loop processes. To demonstrate their efficacy, the proposed approaches are used to an industrial rolling mill benchmark. The task comprises identifying attribute errors in hydraulic reduction devices in both closed-loop & open-loop variants.

In 2020, Wang et al. [34] have attempted to settle a variety of problems. Initially, we create a uniform picture of issues affecting working conditions. The SAE model adapts to the machining signal and extracts tool wear information. For tool wear prediction, the retrieved wear characteristics and associated working circumstances are combined into a working condition comes primarily. Ultimately, the LSTM model’s benefits in solving memory build-up impacts learn the regular wear trend of the perfect working order feature sequence to accomplish tool wear prediction. The efficiency of the developed strategy is demonstrated by experimentation.

In 2021, Xie et al. [35] using the PSO, SSDAE, and LSSVM, researchers have suggested a machined surface state recognition and prediction foundation model. To achieve multi-feature signal discretization having the goal of enhancing prediction accuracy, the SSDAE approach was used, which lowers the reliance on background understanding of feature selection and considerably increases modeling effectiveness. The PSO approach was useful for adaptive kernel parameter optimization, which significantly increased LSSVM model prediction accuracy and processing power. To test the usefulness of the developed approach in enhancing prediction accuracy, a dataset from a real milling operation was used. The experimental findings demonstrated that a high correlation coefficient larger than 0.95 was employed to extract feature vectors from three different directions in the frequency domain, time domain, and time-frequency domain, and that the suggested SSDAE-PSO-LSSVM model outperforms BPNN, PLSR, and ELM with respect to prediction accuracy.

In 2021, Sener et al. [36] have introduced a DCNN-based chatter detection technique. To correctly classify regeneration chatter levels, the technique employs a cardinal model-oriented chatter solution. Vibration data was obtained during milling using a non-invasive data gathering approach. CWT was utilized as a pre-processing approach to disclose essential chatter rich information, taking into account the nonlinear as well as non-stationary aspects of chatter. Following that, the images were utilized to train and evaluate the DCNN that has been constructed. When cutting parameters were incorporated as input features to the DCNN, average accuracy reaches 99.88 percent, according to the technique’s validation.

In 2021, Dun et al. [37] have suggested an unsupervised technique based on vast unlabeled measured dynamic data to evaluate chatter stability in milling. The suggested approach was durable, not subject to measurement mistakes, and does not need labeling. Numerous milling attempts were used to get the dynamic signals. The observed signals were reduced utilizing an auto-encoding-oriented approach in the suggested technique. The compressed signals are then clustered using a hybrid clustering algorithm on the basis of both distance and density metrics. On the experimentally observed dynamic signals, the suggested technique obtains a detection accuracy of 95.6033 percent.

In 2021, Kim et al. [38] have suggested a technique for LAMM of glass MQM. Machining quality was typically impacted by tool condition and machining parameters in tool-oriented mechanical processing, especially LAMM; hence, past research have concentrated on determining monitoring tool condition and optimal machining parameters to ensure machining quality. Previous research, on the other hand, had not taken into account the decline of machining quality with time. Moreover, past research has manually created features from sensory data; however, without sufficient domain expertise in LAMM of glass, these techniques were challenging to use. Since glass materials were more susceptible to break due to mechanical contact among the tool and the work piece, MQM was more crucial in LAMM than it was in metal cutting. We use a new image-oriented deep transfer learning algorithm for MQM in LAMM of glass in this study. The method was dependent on a pre-trained model that has been trained on a large-scale image dataset and was capable of extracting useful characteristics from the images. We present a Multi-layer Recurrence Plot (MRP) that transforms cutting force data into 2D images to graphically depict milling quality. In comparison to various traditional approaches, the suggested MQM method has the highest classification accuracy of machining quality, according to the experimental verification in this study. The suggested approach should be able to anticipate the machining performance of glass micro-milling having greater precision before the machining quality deteriorates.

2.2 Review

STFT [21] is effective for the discrimination within the unstable and the stable tests. But, it is not applied to the intelligent chatter detection. Chatter detection algorithm [22] is better with respect to the capability and reliability in handling the cutting conditions that evolve faster. Still, the nonlinear control algorithm is not taken for the purpose of chatter suppression. Piezoelectric stack actuators [23] greatly suppress the chatter leading to better work piece surface. But, more appropriate milling method is not adopted for the controller design. DTW [24] returns high classification accuracy. Yet, it does not speed the run time for the similarity measures for the purpose of detecting the on-machine chatter. AI [25] updates the fields by novel ideas on the basis of fresh AI approaches and distinct sensor signal features. Still, it does not choose an accurate modelling approach for distinct research process with consideration of the milling techniques. ML [26] employs the periodic online machine learning in the predictive maintenance. Yet, it does not identify the exact meaning of associated coefficients or every coefficient. CNN [27] avoids the cutting tool breakage and the surface roughness created by the chatter in the machining techniques. But, the sensitivity of the produced features is not enhanced. NN [28] provides better chatter detection from the gathered data. Still, it does not examine the SVM and cyclo-stationarity properties in the classification and the processing stages. Adaptive learning control system [29] facilitates the absolute and relative tuning processes. Yet, the models cannot be separately treated owing to the internal coupling. Optimal detection schemes [30] are based on the norm boundedness of the controller and the uncertainties. But, it does not consider the randomized algorithms-oriented threshold setting. SAE and LSTM [34] extract the tool wear features in an adaptive manner. But, the machining quality is impacted by the machining parameters. PSO, SSDAE, and LSSVM [35] enhance the accuracy of prediction. Still, the quality is affected by the tool condition. DCNN [36] labels the regenerative chatter levels in a precise manner. Yet, the machining quality degradation is not considered over time. Unsupervised method [37] diagnoses the chatter stability in the milling. But, only manual features are designed from sensory signals. Deep transfer learning [38] extracts the meaningful features from the images. Still, cracks occur from the mechanical contact among the tool and the work piece. Hence, it is necessary to develop novel machine learning methods for attaining better performance.

Table 1
Features and challenges of traditional milling methods

Author [citation] Methodology Features Challenges

Yun et al. [21] STFT •It is effective for the discrimination within the unstable and the stable tests. •It is not applied to the intelligent chatter detection.

Albertelliet al. [22] Chatter detection algorithm •It is better with respect to the capability and reliability in handling the cutting conditions that evolve faster. •The nonlinear control algorithm is not taken for the purpose of chatter suppression.

Zhang et al. [23] Piezoelectric stack actuators •It greatly suppresses the chatter leading to better work piece surface. •More appropriate milling method is not adopted for the controller design.

Yesilliet al. [24] DTW •It returns high classification accuracy. •It does not speed the run time for the similarity measures for the purpose of detecting the on-machine chatter.

Torabiet al. [25] AI •The fields are updated by novel ideas on the basis of fresh AI approaches and distinct sensor signal features. •It does not choose an accurate modelling approach for distinct research process with consideration of the milling techniques.

Zhang et al. [26] ML •The periodic online machine learning is employed in the predictive maintenance. •It does not identify the exact meaning of associated coefficients or every coefficient.

Serınet al. [27] CNN •It avoids the cutting tool breakage and the surface roughness created by the chatter in the machining techniques. •The sensitivity of the produced features is not enhanced.

Lamraouiet al. [28] NN •It provides better chatter detection from the gathered data. •It does not examine the SVM and cyclo-stationarity properties in the classification and the processing stages.

Guo [29] Adaptive learning control system •It facilitates the absolute and relative tuning processes. •The models cannot be separately treated owing to the internal coupling.

Li and Ding [30] Optimal detection schemes •It is based on the norm boundedness of the controller and the uncertainties. •It does not consider the randomized algorithms-oriented threshold setting.

Wang et al. [34] SAE and LSTM •The tool wear features are extracted in an adaptive manner. •The machining quality is impacted by the machining parameters.

Xie et al. [35] PSO, SSDAE, and LSSVM •The accuracy of prediction is enhanced. •The quality is affected by the tool condition.

Sener et al. [36] DCNN •The regenerative chatter levels are labeled in a precise manner. •The machining quality degradation is not considered over time.

Dun et al. [37] Unsupervised method •The chatter stability in the milling is diagnosed. •Only manual features are designed from sensory signals.

Kim et al. [38] Deep transfer learning •The meaningful features are extracted from the images. •Cracks occur from the mechanical contact among the tool and the work piece.

Author [citation]	Methodology	Features	Challenges
Yun et al. [21]	STFT	•It is effective for the discrimination within the unstable and the stable tests.	•It is not applied to the intelligent chatter detection.
Albertelliet al. [22]	Chatter detection algorithm	•It is better with respect to the capability and reliability in handling the cutting conditions that evolve faster.	•The nonlinear control algorithm is not taken for the purpose of chatter suppression.
Zhang et al. [23]	Piezoelectric stack actuators	•It greatly suppresses the chatter leading to better work piece surface.	•More appropriate milling method is not adopted for the controller design.
Yesilliet al. [24]	DTW	•It returns high classification accuracy.	•It does not speed the run time for the similarity measures for the purpose of detecting the on-machine chatter.
Torabiet al. [25]	AI	•The fields are updated by novel ideas on the basis of fresh AI approaches and distinct sensor signal features.	•It does not choose an accurate modelling approach for distinct research process with consideration of the milling techniques.
Zhang et al. [26]	ML	•The periodic online machine learning is employed in the predictive maintenance.	•It does not identify the exact meaning of associated coefficients or every coefficient.
Serınet al. [27]	CNN	•It avoids the cutting tool breakage and the surface roughness created by the chatter in the machining techniques.	•The sensitivity of the produced features is not enhanced.
Lamraouiet al. [28]	NN	•It provides better chatter detection from the gathered data.	•It does not examine the SVM and cyclo-stationarity properties in the classification and the processing stages.
Guo [29]	Adaptive learning control system	•It facilitates the absolute and relative tuning processes.	•The models cannot be separately treated owing to the internal coupling.
Li and Ding [30]	Optimal detection schemes	•It is based on the norm boundedness of the controller and the uncertainties.	•It does not consider the randomized algorithms-oriented threshold setting.
Wang et al. [34]	SAE and LSTM	•The tool wear features are extracted in an adaptive manner.	•The machining quality is impacted by the machining parameters.
Xie et al. [35]	PSO, SSDAE, and LSSVM	•The accuracy of prediction is enhanced.	•The quality is affected by the tool condition.
Sener et al. [36]	DCNN	•The regenerative chatter levels are labeled in a precise manner.	•The machining quality degradation is not considered over time.
Dun et al. [37]	Unsupervised method	•The chatter stability in the milling is diagnosed.	•Only manual features are designed from sensory signals.
Kim et al. [38]	Deep transfer learning	•The meaningful features are extracted from the images.	•Cracks occur from the mechanical contact among the tool and the work piece.

3 Milling model for robust controller design

3.1 Milling dynamic process

Typical chatter management algorithm designing has seldom taken into account the interaction among spindle dynamic characteristics and machining forces. Such control approaches that employ the spindle dynamic framework as the control model, regard machining / milling forces as exterior intervention. The control block diagram is shown in Fig. 1. The previous chatter suppression approach was designed to reduce the yield of the closed loop spindle dynamic model, resulting in a poorer responsiveness than an open loop framework with the same machining forces. The entire design process is simple, and the overall enhancement may well be accomplished simply raising the stability limits. Moreover, because of the disregard for the link between spindle dynamic characteristics and machining forces, the control system could not guarantee that the machining operation will always be steady for all machining conditions. This paper investigates the relationship amongst spindle dynamic properties and machining forces. The control framework is the machining process dynamic system modelling (that comprises the spindle dynamic framework as well as machining forces), and the control system’s input is the externally active control forces. The goal is to develop a management algorithm, which will enable the machining process’s dynamic model to be stabilised. Since the particular machining parameters are integrated in the control framework architecture, the control scheme could not elevate literally the entire SLD, which could flexibly adjust the SLD and stabilise the overall machining process for the particular machining parameters. Focusing on the physical framework of active chatter reduction and also the mechanical design of the milling process, active control forces might indeed be included into the milling dynamic mathematical equations.

$\begin{matrix} N \ddot{a} (u) + D \dot{a} (u) + La (u) \\ = cI (u) [a (u - τ) - a (u)] + G_{b} (u) \end{matrix}$ (1)

Here, the dynamic cutting coefficients matrix is shown by $I (u) = [\begin{matrix} i_{yy} (u) & i_{yz} (u) \\ i_{zy} (u) & i_{zz} (u) \end{matrix}]$ , the active control forces vector in the y and z directions is shown by G_b (u) = [G_b,y (u) G_b,z (u)] ^T, $N = [\begin{matrix} n_{y} & 0 \\ 0 & n_{z} \end{matrix}]$ , $D = [\begin{matrix} d_{y} & 0 \\ 0 & d_{z} \end{matrix}]$ and $L = [\begin{matrix} l_{y} & 0 \\ 0 & l_{z} \end{matrix}]$ shows the mass, damping, and stiffness matrixes respectively. $τ = \frac{60}{O Ω}$ shows the time delay owing to the regenerative effect, in which Ω and O shows the spindle speed, tooth count, and a (u) = [a_y (u) a_z (u)] ^T shows the dynamic displacement vector in the y and z directions.

Fig. 1

Architecture of chatter control in milling process.

Dynamic milling forces G_E (u) and active control forces G_b (u) wreak havoc on the milling system. Once the controller transmits controlled voltage signals v(u) to the actuators via the amplifier, the actuation unit will function. Apply the proper active control forces to the spindle mechanism in order to eliminate chatter. Whenever the actuator unit varies, the connection across control voltage signals & active control forces alters, which would be an important factor for accurate precise measurements. The suggested control method incorporates 4 piezoelectric layers as actuators. Despite the hysteresis of piezoelectric layers, the feed control device’s transfer function is generally a constant, denoted by L_b. As a consequence, active control forces can indeed be given by $G_{b} (u) = L_{b} v (u)$ (2) The substitution of Equation (2) in Equation (1) is as follows.

$\begin{matrix} N \ddot{a} (u) + D \dot{a} (u) + La (u) = cI (u) \\ [a (u - τ) - a (u)] + L_{b} v (u) \end{matrix}$ (3)

3.2 Modelling the process

In contrast to turning, in which the milling machine is continually in contact with the material, machining is a metal processing procedure wherein the cutting tool approaches and departs the work piece only on occasional basis. Throughout the machining operation, a rotating cutting tool eliminates metal particles from a work piece. The tool rotates & travels in the feed direction at a steady speed. Figure 1 depicts a straightforward dynamic design of 2-DOF end-milling utilizing a flexible tool. This dynamic approach investigates the effect of tool as well as work piece vibrations on instant tooth feed, with the work piece considered to be rigid. A machining procedure may be defined as a sequence of contemporaneous cutting processes performed employing a number of single-point cutting tools. Throughout machining, the chip thickness varies due to the instability & vibration of the machine architecture and cutting equipment. The cutter’s z teeth are supposed to possess a rake angle of 12° and then a helix angle of 30°. Whenever an instrument is flexible and adaptable, its movement between the axes affects the chip load and overall cutting forces, that might cause chatter instability underneath specific conditions. The processing variables are the feed per tooth (also called as chip load), spindle speed, the axial depth-of-cut b and the radial depth-of-cut. The thickness of the chip is not consistent, but moreover intermittent, as a function of the cutter’s rotation as well as feed rate. The purpose of this research is to combine the idea of identical orthogonal cutting situations with the conception of regeneration force models to predict cutting force in the building of a dynamic force framework for end-milling.

3.3 Curtting force description

A machining cutter’s teeth could indeed be viewed as a single-point tool with something like a rotating coordinate framework. Therefore, a forecasting force concept relying on cutting circumstances, tool geometry, work material characteristics and the related instant dynamic regeneration un-deformed uncut thickness of the chip may be used to compute the cutting force in machining. A cyclical cutting force is created whenever a tooth enters or leaves the cutting medium, increasing mostly from 0 to max and afterwards returning to 0 in one direction for each consecutive tooth. The cyclical force looks to have reached its peak. Take the example of down machining, where the machine motion is assumed to be circular. The cutter is subjected to tangential & radial cutting force elements that could very well be reflected in the y&z coordinate axes employing the rotational direction of the cutting edge.

The tangential force G_U is proportionate to the material-dependent specified cutting energy, the feed / tooth, the width of the chip c, as well as the cutting edge’s rotational position is denoted as s. $G_{U} = L_{t} cg sin φ$ (4) The tangential force G_O is supposed to be proportional to the normal or radial force. $G_{O} = l_{s} G_{U}$ (5)

Cutter geometry determines the radial cutting force variable constant, l_s (Proportionality constant). In the y & z axes, G_U and G_U are used to define the transformed to cutting forces. The result of the direction as well as summarization for many of the teeth involved will indeed be $G_{y} = \sum_{j - 1}^{a} G_{Uj} G_{Uj} cos φ_{j} + G_{Oj} sin φ_{j}$ (6) Then, the resulting force in the direction of Z. $G_{z} = \sum_{j = 1}^{a} G_{Uj} sin φ_{j} - G_{Oj} sin φ_{j}$ (7) Where, a is the cutter’s tooth count. Equations 4 & 5 are inserted in to the Equations 6 & 7, resulting in y and z cutting forces, each and every one comprising 3 terms: sine, depth of cut c and cosine of. k s = 0 is applied instead. Equations 5 and 6 may be rewritten in the following way:

$\begin{matrix} G_{y} = L_{t} cg \sum_{j - 1}^{a} sin φ_{j} cos φ_{j} + 0.3 {sin}^{2} φ_{j} \\ = \frac{1}{2} L_{t} cg \sum_{j - 1}^{a} sin 2 φ_{j} + 0.3 - 0.3 cos 2 φ_{j} \end{matrix}$ (8)

$\begin{matrix} G_{z} = L_{t} cg \sum_{j - 1}^{a} {sin}^{2} φ_{j} - 0.3 sin φ_{j} cos φ_{j} \\ = \frac{1}{2} L_{t} cg \sum_{j - 1}^{a} 1 - cos 2 φ_{j} - 0.3 sin 2 φ_{j} \end{matrix}$ (9)

4 Chatter vibration in milling

4.1 Frequency band selection

Whenever the machining operation is unsteady, the dominant vibration frequency levels or greater chatter frequencies are frequently near to the intrinsic frequencies of the machining equipment. Recognition of notable bandwidths connected to natural frequencies could aid in reducing sound in the surroundings and increasing the sensitivity of features extracted from some of those bands. The resonance values of the FRF are used to define dominating frequency bands for chatter identification. The hammer impact experiment is a common way of calculating the FRF, although it requires the usage of additional sensors (including a hammer) for study. Equation (10) describes a synthesis FFT of vibration signals from across all instances of milling situations to find predominant bandwidths in order to avoid the usage of additional sensors. It is difficult to identify dominant bandwidths immediately out of a synthetic FFT. Smoothing the processed FFT spectrum & amplifying the localized peaks employing energy regulators depending on movable frames could help identify dominant bandwidths. Candidates for energy operators include Duplex 2205, Teager non-energy operator and energy regulator, amongst many others. Rather from the other 2 options that depend on either the window function as well as the lag variable, this research utilises the duplex 2205, which only depends on the window function.

The FFT method is used to assess specimen training vibration signals {y_l (u) } (l = 1, 2, ⋯ , L) out of each instance of milling state, regardless steady or unsteady, before summing them to build a synthesis FFT spectrum. $z (g) = \sum_{l = 0}^{L - 1} FFT [y_{l} (u)]$ (10)FFT [y_l (u)] is the variable used to define the FFT of signal y_l (u), and z (g) is the synthesised FFT spectrum, g is the frequency. The Matlab function “fft” is used to implement the FFT algorithm. The duplex 2205 is used to analyse z (g) in order to discover dominating frequency bandwidths with large energy. $Z (g) = F (z^{2} (g)) = z^{2} (g) \otimes x_{1} (g)$ (11)

Here, ⊗ is utilised to denote the convolution operator as well as x₁ (g) is employed to define a smoothing window function. Since the non-dominant bandwidth has a significantly lesser energy value, the squared and smoothed energy in Equation (11) emphasises the high-energy local frequency elements.

A longer window size in the synthesized FFT spectrum may accurately distinguish predominant over non-dominant frequencies. A low window size is inadequate to smooth the synthesised FFT spectrum, but a long window length could unite 2 adjacent dominant bandwidths. As a consequence, the smoothing window’s period must be appropriately set. Rather than the short-frequency energies, a much more broadly adopted squared energy operator is employed for frequency bandwidths selection in this study. The process for determining the upper and lowest boundaries of bandwidths is discussed. As a consequence, cut-off frequencies that are widely used in filter building, are incorporated to determine the limits within each frequency bandwidth. The cut-off frequencies would be those when the squared energy Z (g) decreases to a fraction $1 / \sqrt{2}$ of its local max value. Figure 2 depicts a pictorial representation depiction of the chatter description.

Fig. 2

Description of the chatter model.

4.2 Coefficient of the cutting force

Establishing the cutting force parameter for a cutting process is critical to precisely anticipate the relationship among cutting force as well as milling reliability. 6 Cutting parameters indicate the object’s yield toughness, tool geometry and tool-to-work-material friction. Cutting force values may be computed in 2 distinct methods: (1) employing an orthogonal cutting dataset to obtain the cutting force parameter for bevelled machining, or (2) employing mechanical procedures to promptly derive the cutting force parameter for bevelled machining processes. The initial technique employs a right-angle-relied strategy for oblique snipping transformations. By obtaining orthogonal performance data, cutting force parameters for numerous cutting devices & cutting operations can be calculated. The mean cutting force approach as well as the instant cutting force approach are subgroups of the 2nd approach for calculating cutting force parameters. The averaged cutting force method transforms the averaged cutting force as something like a linear feed / tooth function. Fitting techniques or algorithms in Matlab software could perhaps be employed to calculate force parameters in different directions. The instant cutting force technique uses the minimum fitness function to match a generated cutting force and also the experimentally measured force, using the instant cutting force parameter determined by inversion.

It discussed the laws regulating cutting force coefficient values in regard to cutting parameter values, and even the regulations driving cutting force coefficients in proportion to spindle speed, as well as highlighted however that this had a propensity to become substantially stronger when combined to the tangential force. Based to their observations, the cutting force variable parameter was greater in low-speed zones, subsequently decreased and afterwards increased in high-speed parts. The standard mean force methodology and an instant estimating computation relies on a genetic algorithm were being used to obtain the machining force parameter. Through the simulations results, it is found that the instant cutting force calculation technique was proven to be more precise & productive. In the rated speed of the spindle, the cutting force parameter adjustment regulation would be the same as that suggested somewhere else in the research. Employing the instant cutting force methodology and then a nonlinear optimization method, the impacts of several cutting variables (radial cutting breadth, feed / tooth and spindle speed) on the cutting force coefficient were explored. In respect to the cutting process parameter, the cutting force parameter has been discovered to be non - linear. This indicated that the cutting force parameter was influenced by the cutting variables. Furthermore, investigators studied the effect of cutting variables (depth of the cut and spindle speed) on the cutting parameter & discovered that cutting variables seemed to have no influence on the cutting force parameter in the 500–1500 r/min spindle rated speed.

The process typically comprises evaluating a delay differential problem using a time–frequency domain strategy, followed by the calculation of an SLD, culminating in an appropriate cutting force parameter that could really increase the prediction’s reliability. A stability evaluation method that investigates the continuous changes of the cutting force parameter corresponding to spindle speed has been proposed, considering that the cutting force parameter changes with speed. Conventional techniques to stability estimation, wherein xc utilised to represent the frequency value of the chatter, are far less accurate than this technique. A dynamic machining force parameter framework has lately been developed for machining processes. A dynamic machining force coefficient is used with an upgraded semi-discrete approach to provide an SLD for a process. When contrasted to the results of forecasts dependent on a constant coefficient, the predictability is significantly improved.

4.3 Feature assessment

First as well as second-order image characteristics features could have different capacities for carrying quite so much detail to describe the milling state. ROC curve along with the AUC would be employed to assess the value of image attributes and their ability to differentiate them. Two overlapping q_ω1 (G) and q_ω2 (G) are probability density functions that describe the distribution of a single characteristic G into two categories ω₁, and ω₂, as well as G₀ is a threshold. Values to the left of the threshold are believed to be negative. Values present in the left fit in to category ω₁, whereas values present in the right fit in to class ω₂. As a result, there are two probability of making a mistake in category ω₁ and class ω₂, respectively: α (G₀) and β (G₀). $α (G_{0}) = {smallint}_{G_{0}}^{\infty} q_{ω 1} (G), β (G 0) = {smallint}_{- \infty}^{G_{0}} q_{ω 2} (G)$ (12)

With G₀ as the variable, the ROC curve shows 1 - β (G₀) vs α (G₀) parametrically. For every threshold, α (G₀) = 1 - β (G₀) if the two probability density functions completely overlap. The associated ROC curve deviates from the straight line α (G₀) = 1 - β (G₀) as the two distributions move apart.

The AUC is the shaded area seen across ROC as well as the horizontal straight line. This area fluctuates between 0 and 1/2 for complete overlap and separation. It’s a metric for a feature’s capacity to distinguish between classes [10]. A characteristic with a high AUC is good at distinguishing between two groups, whereas one with a low AUC is bad at discriminating. As a result, the time-frequency image characteristics are sorted based on their AUCs.

5 Deep learning methods

5.1 LSTM

The LSTM [31] has a candidate layer CSNN with tanh, a forget gate FS, an input gate IS and an output gate OS that are NN with a sigmoid function, a memory state MS, and a hidden state HS that is a vector. Assume that the input blocks’ memory is cs_(ts-1), the input vector is xs_ts, and the previous output of the blocks is ms_(ts-1). The current block’s output is denoted by ms_ts, and the current block’s memory is denoted by cs_ts. Three inputs are accepted by the network. The decision is made by a single unit, which considers “prior output, present input, and previous memory.”

In the underground levels, memory bricks, which are special blocks, are employed. Input and output gates are used in every memory block. During the activation of both the input and output, control operations are carried out. In the final step, there is also a forget gate. Mapping is the process of converting an input sequence xs = (xs₁, xs₂, ⋯ , xs_TS) to an output sequence ys = (ys₁, ys₂, ⋯ , ys_TS). The activations of network units are indicated in the equations below.

$\begin{matrix} {is}_{ts} = σ (W S_{isxs} x s_{ts} + W S_{isms} m s_{ts - 1} \\ + W S_{iscs} c s_{ts - 1} + b s_{is}) \end{matrix}$ (13)

$\begin{matrix} f s_{ts} = σ (W S_{fsxs} x s_{ts} + W S_{fsms} m s_{ts - 1} \\ + W S_{fscs} c s_{ts - 1} + b s_{fs}) \end{matrix}$ (14)

$\begin{matrix} c s_{ts} = f s_{ts} \otimes c s_{ts - 1} + i s_{ts} \otimes gs (W S_{csxs} x s_{ts} \\ + W S_{csms} m s_{ts - 1} + b s_{cs}) \end{matrix}$ (15)

$\begin{matrix} o s_{ts} = σ (W S_{osxs} x s_{ts} + W S_{osms} m s_{ts - 1} \\ + W S_{oscs} c s_{ts - 1} + b s_{os}) \end{matrix}$ (16) $m s_{ts} = o s_{ts} \otimes hsc s_{ts}$ (17) $y s_{ts} = φ (W S_{ysms} m s_{ts} + b s_{ys})$ (18)

The greatest weight from the input gate to the input is denoted by WS_isxs, while the weight matrices are denoted by WS, WS_iscs, WS_fscs, and WS_oscs denote the diagonal weights associated with the peep-hole connections, respectively. The sigmoid function is defined with the help of σ, bs_is is employed to denote the input gate bias vector, the forget gate is represented by fs, the input gate is designated by the variable is, the cell activation vector is denoted with the help of cs, and the output gate is defined with the variable os, as appropriate. The letters hs, gs, and φ denote the “cell output function,” “cell input function,” and “network output initiation function,” respectively. The initiation/ activation layer is responsible for deciphering the complicated arrangement or structure. The “hyperbolic tangent tanh and sigmoid σ functions” are employed in the multi-layer LSTM. The basic LSTM architecture is shown in Fig. 3.

Fig. 3

LSTM model.

5.2 DTW

DTW [32] is a technique for evaluating the closeness or correlation across the 2 time series, regardless of respective durations. Let TT₁ and TT₂ be two time series with the following components y_j and z_k lengths n and o. $T T_{1} = y_{1}, y_{2}, \dots, y_{j}, \dots, y_{n}$ (19) $T T_{2} = z_{1}, z_{2}, \dots, z_{k}, \dots, z_{o}$ (20)

The warping path ω_l = (y_j(l), z_k(l)) among the 2 time series may be expressed by mapping the time series’ respective components or elements on a n × o matrix. The warping path is made up of points ω_l that show alignment between the time series elements y_j(l) and z_k(l). The length M of the warping route satisfies the criteria n ⩽ M ⩽ o, where o ⩾ n is assumed. The diagonal line via n × o matrix when n = o, i.e. ω_l = (y_l, z_l) is always the warping path. However, warping pathways are not always unique, and for the same two time series, several warping paths can be formed. The DTW method selects the warping route that produces the least gap across element/ component pairs of differing lengths when certain constraints are met. Although there are a variety of approaches for computing the distance among 2 time series components (y_j, z_k), the Euclidean distance d (y_j, z_k) = y_j - z_k ₂ is employed in this version. The DTW algorithm’s minimization of the distance between TT₁ and TT₂ may thus be stated as Equation (21). $DTW (T T_{1}, T T_{2}) = min (\sum_{l = 1}^{M} d (w_{l}))$ (21)

Different warping pathways are subject to a number of constraints. Slope constraint, monotonicity, boundary conditions, window condition continuity and adjustment are all examples of these. Because there are so many alternatives for warping paths without even any limits, these parameters or constrains are fixed to the alignment window to reduce the range of possible warping routes.

Monotonicity: The indices j and k should always rise or remain constant, resulting in j (l) ⩾ j (l - 1) and k (l) ⩾ k (l - 1).

Continuity: The indices j and k can only rise by one, making j (l) - j (l - 1) ⩽ 1 and k (l) - k (l - 1) ⩽ 1 respectively.

Boundary condition: The warping routes should begin with j and k equal to 1 and conclude with j and k equal to o and n.

Adjustment window condition: To prevent significant temporal differences between the two pathways, the warping path with the shortest distance is sought over a restricted portion of the alignment window. j - s ⩽ k ⩽ j + s represents the limited area.

Slope constraint: This condition prevents significant one-way movement. It cannot travel in the same direction after taking b horizontal or vertical step without first taking c diagonal steps. The slope constraint’s effective intensity can be denoted as Q = c/ b. Q = 1 was chosen since it was shown to be the best value in a voice recognition experiment.

5.3 Bi-LSTM

A Bi-LSTM [33] layer is used to understand bidirectional long-term relationships across time scales of time sequence or data series. Such connections may be useful whenever it is necessary to understand from the complete time sequence through each time unit by the network. Bidirectional LSTMs are little more than a LSTM, which could really help improve efficiency of the system when dealing with sequence categorization problems. Bidirectional LSTM networks educate 2 LSTMs rather than just one for the input sequences. The first one is focused on the actual original input sequences, whereas the other is relied on a reversed/ inverted replica. It can also provide additional perspective for the network, allowing it to comprehend the issue very promptly and comprehensively.

The characters j designate input gate in each and every BI-LSTM neurons, g define forget gate in each and every BI-LSTM neurons and p designate output gate in each and every BI-LSTM neurons. The inner anatomy of Bi-LSTM neurons are depicted in Fig. 4. The input gate j controls the quantity of prevailing connection input data y_u, which would be ended up saving to the storage unit D_u, helps to evaluate novel data to be protected utilising the sigmoid function, actually creates a fresh applicant vector ${\tilde{D}}_{u}$ the tanh function, & transfer the message to be stored to the storage / memory unit D_u. The unit will be the one who completes the enhancement. The forget gate g supervises the self-connecting unit and also it screens the data in storage device D_u-1 at the prior moment to determine the quantity of genuine info, which should be conserved in the current storage device D_u simultaneously discarding the worthless data. The output gate p describes the amount of data produced by the storage device D_u at time step u as well as regulates the storage device D_u’s influence on the actual output variable i_u.

Fig. 4

Bi-LSTM model.

The formula is defined as:

j_{u} = σ (X_{yj} y_{u} + X_{ij} i_{u - 1} + c_{j})

(22)

{\tilde{D}}_{u} = tan h (X_{yd} y_{u} + X_{id} i_{u - 1} + c_{d})

(23)

g_{u} = σ (X_{yg} y_{u} + X_{ig} i_{u - 1} + c_{g})

(24)

D_{u} = g_{u} Θ d_{u - 1} + j_{u} Θ {\tilde{D}}_{u}

(25)

p_{u} = σ (X_{yp} y_{u} + X_{ip} i_{u - 1} + c_{p})

(26)

i_{u} = p_{u} Θ tan h (D_{u})

(27) where D denotes the storage device, typically referred as the state of the cell, D_u represents the storage cell state at time step u,

{\tilde{D}}_{u}

represents the storage cell’s applicant vector at time step u, y_u represents time step u ‘s input vector at time step u, i u represents the time step u ‘s output vector, X represents the network’s weight vector, as well as c represents the offset vector. Θ is used to define a multiplication function of vector elements. The high-dimensional characteristic of the input timed signal produces the forward LSTM network vector

{\vec{i}}_{u}

, whilst the inverted LSTM network vector

i_{u}^{\leftarrow}

produces the low-dimensional feature, BiLSTM network eigenvalue output vector at time step u seems to be Q_u, as well as σ (·) is utilized to represent the sigmoid function. The following is the information.

{\vec{i}}_{u} = \vec{LSTM} (i_{u - 1}, y_{u}, D_{u - 1})

(28)

i_{u}^{\leftarrow} = \leftarrow LSTM (i_{u + 1}, y_{u}, D_{u + 1})

(29)

Q_{u} = [{\vec{i}}_{u}, {\overset{\leftarrow}{i}}_{u}]

(30)

The flow of the process taking place in the Bi-LSTM is given in the pseudo code given below,

Pseudo code of the Bi-LSTM process
Start
Initialize the parameters
Train Bi-LSTM with the input dataset
Detect the chatter in the given data
Calculate the RSME and loss function
End

6 Result and discussion

6.1 Experimental setup

The Milling process is done with the Duplex 2205 stainless steal material and the parameters namely speed, depth of the cut and feed is given as input for the milling machine. Some of the parameters like MRR, actual cutting force and acceleration are calculated from the milling process. The obtained acceleration data is converted as vibration data for better learning and easy detection of chatter. The DL methods adapted for the vibration analysis and to identify the chatter. The Matlab software utilized for the implementation purpose. The output obtained from the simulation is analysed for identifying the well performing DL among them. Here, the considered DL techniques are DTW, LSTM and Bi-LSTM.

6.2 Vibration analysis and its purpose

Vibration reduction can undoubtedly contribute to highly efficient as well as productive milling processes. That ensures machinists could manufacture many high-quality components, lower scrap levels, and increase profits. Vibration analysis plays a part in this scenario. Vibration assessment is the process of obtaining knowledge regarding vibration in order to determine the best, most consistent speeds, feeds and other process parameters for milling objects. From a news piece on Manufacturing.net, “Vibration analysis may indeed be the ears for reading the vernacular of machinery and can assist us identify wherein to spend time and money.”

It’s similar to checking the fitness of the assembly of the milling equipment, that is comprised of the milling cutter, machine spindle and toolholder. More and more vibration & noise your milling assembly experiences, the more urgent the requirement for speed and feed adjustments to optimize your components and prevent the arrangement from failure. It’s similar to checking the condition of the machining equipment, which is comprised of the machine spindle, toolholder, and milling cutter. Therefore more vibration and noise your milling assembly experiences, the more urgent the requirement for speed & feed adjustments to optimize the components and prevent the system from failures. Vibration (detected through accelerometers and sensors) is transformed to an electric current and monitored as a signal. Vibration analysis examines velocity, vibration frequency, displacement and acceleration measure and delivers insights. The signal data obtained from your machining setup may then be represented and compared using four dimensions. Time domain, joint domain, frequency domain and modal analysis are the 4 main features mentioned above. According to ScienceDirect, “each domain gives precise data on the operating circumstances and characteristics of the vibrating portion.”

6.3 Experimental results of DTW

For chatter identification, DTW using DL combines the similarity measurements of time series for the data stream. The DTW distance is used to compare the similarities of various time series. For DTW, the RMSE is 0.0035. For surface structuring purposes, creating or maintaining intermediary chatter may be preferable. The fact that the power spectrums for chatter as well as intermediate chatter are so similar makes frequency domain featurization incredibly difficult. In the temporal domain, there is a clear distinction among the 2 chatter regimens. As a result, retrieving unique features and attributes in the time domain for a 3-class categorization (no chatter, moderate chatter, and chatter) appears to be somewhat more advantageous. The accelerometer vibration data obtained by the DTW is given in the Fig. 5 (a). The Fig. 5(b) shows the results of the chatter error detection using the accelation value obtained. Figure 5(c) and 5(d) depicts the training process of the DTW and the confusion matrix, respectively.

Fig. 5

DTW results, (a) Accelerometer vibration data, (b) chatter error detection, (c) training progress of DTW and (d) Confusion matrix of DTW.

6.4 Experimental results of LSTM

To obtain higher-level characteristics w ith improved test categorization, the LSTM framework captures long-term dependencies and maintains strong classification precision for the machining operation. The value of RMSE for LSTM is 0.0030. The Fig. 6 (a) . gives the accelerometer vibration data obtained by the LSTM. From the accelation value obtained the chatter error detection is obtained and shown in Fig. 6(b). Figure 6(c) and 6(d) depicts the training process of the LSTM and the confusion matrix, respectively.

Fig. 6

LSTMresults, (a) Accelerometer vibration data, (b) chatter error detection, (c) trainning progress of LSTM and (d) Confusion matrix of LSTM.

6.5 Experimental results of Bi-LSTM

The accelerometer vibration data obtained by the Bi-LSTM is given in the Fig. 7 (a). The Fig. 7(b) shows the results of the chatter error detection using the accelation value obtained. Figure 7(c) and 7(d) depicts the training process of the Bi-LSTM and the confusion matrix, respectively.The value of RMSE for Bi-LSTM is 0.0011. The details of the different layers in the Bi-LSTM and specific training parameters of the model is tabulated and presented in the Tables 2 and 3, respectively.

Fig. 7

Bi-LSTMresults, (a) Accelerometer vibration data, (b) chatter error detection, (c) trainning progress of Bi-LSTM and (d) Confusion matrix of Bi-LSTM.

Table 2

Details of the different layers

Layer Name	Feature Size	No of layer
of the Output
Input layer	12*125	2
Bi-LSTM Layer	2	1
FullyConnectedLayer	125	1
Regression layer	3	1

Table 3

Specific training metrics of the system

Training parameters
Rate of learning	0.001
Epoch	100
Optimizer	Adam
Dropout	0.5
Batch Size	16

6.6 Comparitive analysis

The RMSE value obtained for the considered three DL techniques are presented in the Table 4. The root mean square error (RMSE) is a typical means of calculating a systems’s error in forecasting probabilistic results. The prediction error also named as residual, whose standard deviation is denoted as RMSE. Residuals are indeed a metric about how far the measured values seem to be from the line of regression. RMSE is a measurement of how scattered out these residuals are. In simple terms, it indicates how densely the values is clustered all along line of best fit. In meteorology, predicting, and simple regression, the RMSE is extensively utilised check results of the experiment. The RMSE is a popular metric for comparing forecasted and actual data by a system or predictor.The mathematically derived formula for calculating the RMSE is gven in the Equation (32).

Table 4
Comparision of RMSE of the DL methods

Dl Methods	RMSE Value
DTW	0.0035
LSTM	0.0030
Bi-LSTM	0.0011

rmse = \sqrt{\bar{{(f_{p} - o_{a})}^{2}}}

(31) where, the predicted or forecasted value obtained is denoted by the variable f_p and the observed, actual or standard value is defined as o_a. The bar given above the square of the difference between the predicted data and standard data represents the mean of it.

In the chatter detection process using the various DL process it is eviedent from the Table 5 that the Bi-prediction LSTM’s accuracy was much greater than DTW and LSTM’s. When contrasted to the existing algorithm approaches, this technique performed significantly better. The Table 5 shows the comparition ofloss function along with verification accuracy set of the three methods considered.

Table 5

Comparision of loss function & verification accuracy set of the DL methods

DL methods	Loss	Rate of Verification Precision/ Accuracy	Single Test (Time, m/Sec)	Rate of Test Accuracy%
DTW	0.2688	88.34%	2	87.57
LSTM	0.2857	87.13	20	86.36
Bi-LSTM	0.030	99.590	5	99.5

7 Conclustion

In milling processes, chatter is a self-excited vibration of components. It may be found in a variety of cutting operations and has an influence on both efficacy and effectiveness in manufacturing. The evolution of technology that can anticipate and identify chatter has been the focus of a lot of study. In this paper, three well-performing deep learning approaches such as LSTM, DTW, and Bi-LSTM were used to mill Duplex 2205 material and provide an effective way for monitoring and managing milling processes along with identifying the presence of chatter. A basic milling architecture was evaluated, and the selected three DL techniques were put to the test to see which one performs the best. The Bi-LSTM outperformed other approaches in detecting chatter present by attaining 99.59% verification accuracy rate on comparing with other two methods.

References

Jin

and Xie

, Experimental study on surface generation in vibration-assisted micro-milling of glass, The International Journal of Advanced Manufacturing Technology 81(1-4) (2015), 507–512.

Geng

, Zhang

, Xu

, He

, Liu

and Duan

, Rotary ultrasonic elliptical machining for side milling of cfrp: Tool performance and surface integrity, , Ultrasonics 59 (2015), 128–137.

Zhu

and Liu

, Recent progress of chatter prediction, detection and suppression in milling,. , Mech. Syst. Signal Process 143 (2020).

, Cao

and Chen

, Fuzzy control of milling chatter with piezoelectric actuators embedded to the tool holder,. , Mech. Syst. Signal Process 148 (2021).

Zhou

and Xue

, A multisensor fusion method for tool condition monitoring in milling, Sensors 18(11) (2018), 3866.

Tsai

N.-C.

, Chen

D.-C.

and Lee

R.-M.

, Chatter prevention for milling process by acoustic signal feedback, Int. J. Adv. Manuf. Technol. 47(9–12) (2010), 1013–1021.

Lamraoui

M.M.E.F.

Chatter detection in CNC milling processes based onWiener-SVM approach and using only motor current signals, in Vibration Engineering and Technology of Machinery. Cham, Switzerland: Springer (2014), 567–578.

Zhang

, Anders

, Löser

, Ihlenfeldt

, Czarske

and Kuschmierz

, Non-contact, bi-directional tool tip vibration measurement in CNC milling machines with a single optical sensor,. , Mech. Syst. Signal Process 139 (2020).

, Song

, Liu

, Wang

and Wan

, Size-dependent responses of micro-end mill based on strain gradient elasticity theory, Int. J. Adv. Manuf. Technol. 100(5–8) (2019), 1839–1854.

10.

Huang

, Li

, Sun

and Ge

, Milling force vibration analysis in high-speed-milling titanium alloy using variable pitch angle mill, Int. J. Adv. Manuf. Technol. 58(1-4) (2012), 153–160.

11.

, Jia

, Wang

, Liu

, Feng

, Liang

S.Y.

, Measurement and prediction of vibration displacement in micro-milling of nickelbased superalloy, , Measurement 145 (2019), 254–263.

12.

Yue

, Gao

, Liu

, Liang

S.Y.

and Wang

, A review of chatter vibration research in milling, Chin. J. Aeronaut 32(2) (2019), 215–242.

13.

Liu

, Zhang

and Wu

, An efficient full-discretization method for prediction of milling stability,. , Int. J. Mach. Tools Manuf 63 (2012), 44–48.

14.

Madoliat

, Hayati

and Ghasemi

, Ghalebahman, Investigation of chatter suppression in slender endmill via a frictional damper, Scientia Iranica 18(5) (2011), 1069–1077.

15.

Zhang

, Li

, Meng

, Tu

and Cheng

, Chatter detection in milling process based on the energy entropy of VMD and WPD,. , Int. J. Mach. Tools Manuf 108 (2016), 106–112.

16.

Tran

M.-Q.

, Liu

M.-K.

and Tran

Q.-V.

, Milling chatter detection using scalogram and deep convolutional neural network,. , Int. J. Adv. Manuf. Technol 107(3-4) (2020), 1505–1516.

17.

Paul

and Morales-Menendez

, Active control of chatter in milling process using intelligent PD/PID control, , IEEE Access 6 (2018), 72698–72713.

18.

Wang

, Song

and Liu

, Vibration suppression of thin-walled workpiece milling using a time-space varying PD control method via piezoelectric actuator,. , Int. J. Adv. Manuf. Technol 105(7-8) (2019), 2843–2856.

19.

Ding

, Zhu

, Zhang

and Ding

, A full-discretization method for prediction of milling stability,. , Int. J. Mach. Tools Manuf 50(5) (2010), 502–509.

20.

Altintas

, Stepan

, Merdol

and Dombovari

, Chatter stability of milling in frequency and discrete time domain,. , CIRP J. Manuf. Sci. Technol 1(1) (2008), 35–44.

21.

Yun Chen , Huaizhong Li , Liang Hou , Xiangjian Bu , Feature extraction using dominant frequency bands and time-frequency image analysis for chatter detection in milling, , Precision Engineering 56 (2019), 235–245.

22.

Paolo Albertelli , Luca Braghieri , Mattia Torta , Michele Monno , Development of a generalized chatter detection methodology for variable speed machining, , Mechanical Systems and Signal Processing 123 (2019), 26–42.

23.

Xingwu Zhang , Chenxi Wang , Jinxin Liu , Ruqiang Yan , Hongrui Cao , Xuefeng Chen , Robust active control based milling chatter suppression with perturbation model via piezoelectric stack actuators, , Mechanical Systems and Signal Processing 120 (2019), 808–835.

24.

Melih Yesilli

, Firas Khasawneh

, Andreas Otto , Chatter Detection in Turning Using Machine Learning and Similarity Measures of Time Series via Dynamic Time Warping, arXiv, August 2019.20.

25.

Torabi

A.J.

et al., A Survey on Artificial Intelligence-Based Modeling Techniques for High Speed Milling Processes, in IEEE Systems Journal 9(3) (2015), 1069–1080.

26.

Zhang

Y.G.

, Wu

Q.H.

, Wang

, Oluwande

, Matts

and Zhou

X.X.

, Coal mill modeling by machine learning based on onsite measurements, , in, IEEE Transactions on Energy Conversion 17(4) (2002), 549–555.

27.

Gokberk Serın , Ugur Gudelek

, Batıhan Sener , Murat Ozbayoglu

, Hakki Ozgur Unver , CHATTER DETECTION IN MILLING WITH ACOUSTIC EMISSIONSANDDEEP LEARNING, 10th INTERNATIONAL CONGRESS ON MACHINING, November 2019.

28.

Mourad Lamraoui , Barakat

, Marc , Mohamed El badaoui , Chatter detection in milling machines by neural network classification and feature selection, Journal of Vibration and Control (2013).

29.

Guo

, Development of an optimal crown/shape level-2 control model for rolling mills with multiple control devices, in, IEEE Transactions on Control Systems Technology 6(2) (1998), 172–179.

30.

and Ding

S.X.

, Optimal Detection Schemes for Multiplicative Faults in Uncertain Systems With Application to Rolling Mill Processes, in, IEEE Transactions on Control Systems Technology 28(6) (2020), 2432–2444.

31.

Muhammad Abbasi

, Anum Rashad , Anas Basalamah , Muhammad Tariq , Detection of Epilepsy Seizures in Neo-Natal EEG Using LSTM Architecture, , IEEE Access 7 (2019).

32.

Yesilli

M.C.

, Khasawneh

F.A.

and Otto

, Chatter Detection in turning using machine learning and similarity measures of time series via dynamic time warping, Journal of Manifacturing Process 77 (2022), 190–206.

33.

Qipeng Chen , Qingsheng Xie , Qingni Yuan , Haisong Huang , Yiting Li , Research on a Real-Time Monitoring Method for the Wear State of a Tool Based on a Convolutional Bidirectional LSTM Model, , Symmetry 11 (2019).

34.

Wang

, Zhou

, Gao

, Li

and Li

, Milling Tool Wear Prediction Method Based on Deep Learning Under Variable Working Conditions, IEEE Access 8 (2020), 140726–140735.

35.

Xie

, Zhang

and Liu

, Tool Wear Status Recognition and Prediction Model of Milling Cutter Based on Deep Learning , IEEE Access 9 (2021), 1616–1625.

36.

Sener

, Gudelek

M.U.

, Ozbayoglu

A.M.

and Unver

H.O.

, A novel chatter detection method for milling using deep convolution neural networks, Measurement 182 (2021).

37.

Dun

, Zhu

, Yan

and Wang

, A chatter detection method in milling of thin-walled TC4 alloy workpiece based on auto-encoding and hybrid clustering, Mechanical Systems and Signal Processing 158 (2021).

38.

Kim

, Kim

, Youn

B.D.

and Ahn

S.H.

, Machining quality monitoring (MQM) in laser-assisted micro-milling of glass using cutting force signals: an image-based deep transfer learning, Journal of Intelligent Manufacturing (2021).

Comparative analysis on recent deep learning techniques for identifying chatter in milling process

Abstract

Keywords

1 Introduction

2 Literature survey

2.1 Related works

2.2 Review

3.1 Milling dynamic process

3.3 Curtting force description

4.1 Frequency band selection

4.3 Feature assessment

5.1 LSTM

6.1 Experimental setup

6.2 Vibration analysis and its purpose

6.3 Experimental results of DTW

Table 4 Comparision of RMSE of the DL methods

References

Table 4
Comparision of RMSE of the DL methods