Issue 
Wuhan Univ. J. Nat. Sci.
Volume 27, Number 5, October 2022



Page(s)  439  452  
DOI  https://doi.org/10.1051/wujns/2022275439  
Published online  11 November 2022 
Information Technology
CLC number: TG 506.9
Study on the Predictive of Dynamic Milling Force of Milling Process Based on Data Mining
School of Mechanical and Electrical Engineering, Lanzhou University of Technology, Lanzhou 730050, Gansu, China
^{†} To whom correspondence should be addressed. Email: 18846431603@163.com
Received:
11
May
2022
To improve surface accuracy of the workpiece and obtain potentially valuable information, a dynamic milling force prediction model was proposed based on data mining. In view of the current dynamic milling force obtained through finite element simulation and analytical calculation, in the finite element modeling, the model built is inevitably different from the actual working conditions, and the analytical calculation is slightly cumbersome and complex, and a dynamic milling force prediction model based on data mining is proposed. The model was established using a combination of regression analysis and Radial Basis Function (RBF) neural network. Using data mining as a means, the internal relationship between milling force, cutting parameters, temperature, vibration and surface quality is deeply analyzed, and the influence of dynamic milling force changes on different situations is extracted and summarized by the methods of cluster analysis and correlation analysis. The results show that the proposed dynamic milling force model has a good prediction effect, ensures the production quality, reduces the occurrence of flutter, improves the surface accuracy of the workpiece, and provides a more accurate basis for the selection of process parameters.
Key words: data mining / data modeling / dynamic milling force / prediction of milling force
Biography: JIN Lan, female, Associate professor, research direction: system dynamics. Email: lan_jane@lut.edu.cn
Fundation item: Supported by Gansu Science and Technology Program(21YF5GA080)
© Wuhan University 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
0 Introduction
In recent years, the deep integration of information technology and manufacturing industry has become one of the ultimate goals of the manufacturing industry to promote "smart manufacturing"^{[1]}. Datadriven predictive and proactive manufacturing is beginning to permeate all aspects of the machining process^{[2]}.
Judging and predicting workpiece surface quality by monitoring signals related to vibration, cutting forces, temperature and etc. is an important aspect of intelligent cutting^{[3]}. The quality, efficiency and competitiveness of milling can be greatly improved by converting the large amount of relevant data contained in the milling process, such as vibration, cutting forces and temperature, into useful information, exploring its intrinsic relevance to actual machining problems, and applying methods such as cluster analysis, correlation analysis, regression analysis, and artificial neural networks to generate new and unusual information on the basis of the available data sets to guide production^{[4]}.
As early as 1985, Wang et al^{[5]} used fuzzy cluster analysis to identify the cutting state, and this method has been widely used in data mining. Cai et al^{[6]} used cluster analysis to classify and identify the energy efficiency status by applying image processing, and invented a method to locate and trace the source of low energy efficiency problems occurring in cutting and machining systems, providing a theoretical basis and a potential technical solution for online monitoring of the energy efficiency status of cutting and machining systems. Bracke et al^{[7]} outlined an integrated approach to detect grinding parameter settings which ensures repeatability on the one hand and economy on the other hand. Cluster analysis was used to identify groups of grinding parameters and grinding results with similar characteristics.
Based on the connection function, Pei et al^{[8]} analyzed the correlation of force signal, vibration signal and roughness signal in milling process, and got a better result. Rodrigues et al^{[9]} applied correlation analysis to explore the degree of correlation between roughness and cutting speed and feed rate when turning 45gauge steel in order to find the optimal configuration of the parameters. Thasthakeer et al^{[10]} investigated the correlation between the variation of different parameters and the rate of change of cutting forces and surface roughness after machining. They found that this correlation was higher under wet machining conditions. Based on the multiple regression analysis, Lin et al^{[11]} established the prediction model of the cutting parameters and the machining vibration of the machined parts in face milling by the surface roughness, and got a good prediction result. In Ref. [12], to explore the relationship between machining parameters and cutting temperature, the analysis of variance was used to test the significance of the parameters and a regression equation was developed to demonstrate the accuracy of the predictive model by comparing predicted and experimental values of cutting temperature. Rajesh et al^{[13]} used multivariate regression analysis, which included a variance analysis, to establish a predictive model for the relative effects of process parameters on incision width and surface roughness. Yeganefar et al^{[14]} used support vector machine, the neural network and regression analysis to predict and optimize the surface roughness and cutting forces during milling, and compared them with each other in order to obtain a selfuse and optimal selection scheme. Maher et al^{[15]} used an adaptive neurofuzzy inference system model to study the degree of correlation between cutting parameters, cutting forces and surface roughness. Basak et al^{[16]} applied Radial Basis Function (RBF) neural network to milling, and fitted RBF neural network model by 27 kinds of experimental data under different cutting conditions, which showed that this method can be applied to milling.
The advantage of data mining is the introduction of cluster analysis, correlation analysis and regression analysis which are efficient, convenient and accurate in solving practical machining problems. This method could start from the actual milling process and take cluster analysis, correlation analysis and regression analysis and RBF neural network as the main means to build a dynamic milling force prediction model. Results demonstrate that the proposed dynamic milling force model has a better prediction effect, which ensures production quality and reduces the occurrence of chattering to provide a more accurate process parameters selection basis.
1 Milling Data Gain with Multiple Sensors
In order to obtain the necessary experimental data, in this paper, the complex milling experiments including surface quality detection before milling, milling process, milling temperature detection, milling vibration detection, dynamic milling force detection and milling process are designed. The experimental scheme is shown in Fig. 1.
Fig. 1 Experimental scheme design idea
Experimental scheme design idea 
1.1 Design of the Milling Experimental Programme
In this paper, 12CrNi2 was selected as the experimental material (as shown in Fig. 2). This material belongs to alloy steel, which contains nickel and has certain sticky characteristics. Therefore, in order to make the cut surface have a high surface finish, a fourflute tungsten steel monolithic alloy end mill with the size of D12 mm×75 mm×4F is selected, as shown in Fig. 3.
Fig. 2 Workpiece after rough machining
Workpiece after rough machining 
Fig. 3 Integral alloy end mill
Integral alloy end mill 
Considering the influence of chatter on the actual cutting process before setting up the process parameters, the cutting stability region is found through the milling stability analysis in advance (Fig. 4). The range of process parameters is determined by combining the process manual.
Fig. 4 Stability lobe diagram
Stability lobe diagram 
The technological parameters are as follows: cutting width a_{e}=0.56 mm, cutting depth a_{p}=0.12 mm, cutting speed v_{c}=72140 m/min, feed per tooth f_{z} =0.0250.063 mm, rotational speed n=1 0004 000 r/min, and f is the feed speed. Then the response surface method is used to design the concrete parameter scheme. The specific cutting process parameters are shown in Table 1.
Process parameter scheme
1.2 Experimental Data Acquisition
Before milling, one end of the three acceleration sensors is fixed at the spindle handle and fixture, respectively. The other end is connected to an LMS tester connected to a laptop computer, so as to collect the vibration of the spindle tool holder and the vibration of the workpiece itself during the milling process, and prepare for subsequent data processing. In order to measure the change of cutting temperature with time during the milling process, we focus the infrared thermal image on the tool tip, set a moderately sized temperature collection frame, select the realtime highest temperature in the display collection frame, and use this temperature as the tool tip for milling temperature to realize realtime monitoring of milling temperature. The milling method is down milling. To avoid the influence of cutting fluids on cutting temperature variations, a dry cutting machining method is chosen. During cutting, the pretreated workpiece is tightly fixed to the fixture of the pressuretype threeway force test sensor, the surface of the workpiece is gently tapped with a copper rod to keep its level, and the data on the realtime variation of the cutting force is recorded by this sensor. The experimental conditions are shown in Fig. 5.
Fig. 5 Experimental equipment installation
Experimental equipment installation 
When the experiments were completed, roughly 3.5×10^{8} experimental data were gained. The resulting data were collated to obtain the experimental results. Taking the first set of experiments as an example, the data images are shown in Figs. 69.
Fig. 6 Roughness curve graph
Roughness curve graph 
Fig. 7 Dynamic cutting force
Dynamic cutting force 
Fig. 8 Infrared thermal imaging (a) and temperature change curve(b)
Infrared thermal imaging (a) and temperature change curve(b) 
Fig. 9 Original vibration time domain graph
Original vibration time domain graph "Amplitude" is acceleration amplitude measured by accelerometer 
Figure 6 is the first set of roughness curves for the first piece of workpiece. According to the distribution of cutting traces, the surface of each workpiece is divided into 15 areas, and the roughness is measured separately,and the data records are saved. The SJ400 measuring in strument under the Surftest brand is used to measure surface roughness as shown in Fig. 10. The dynamic cutting force curve of the first set of experiments is shown in Fig. 7. In the first set of experimental cutting processes, the maximum temperature is 103.7 ℃, the minimum temperature is 50.6 ℃, and the average temperature is 87.8 ℃. Figure 8 shows the infrared thermal imaging and the temperature change curve. Figure 9 is a timedomain plot of the original vibration in the X, Y, Z directions measured by the accelerometer attached to the spindle.
Fig. 10 Surface roughness measuring instrument
Surface roughness measuring instrument 
1.3 Experimental Data Preprocessing
Before data mining can be performed on a sample of data, the data are usually needed to preprocess, i.e. clean, integrate, transform and normalize.
As the samples used for data mining in this paper are derived from experiments where the data acquisition sources are the respective test equipment, the acquired roughness, temperature, cutting force and vibration data have been preclassified and placed at the time of acquisition. These data are imported into EXCLE tables separately and integrated to complete the integration work; as the data samples are collected under specific conditions and the collection process is supervised by personal, there is no missing data, misstorage and duplication, but the data signal itself has trend terms and noise, so the data with numerical anomalies are needed to eliminate from the samples through the data screening function of EXCLE before analysis to obtain a more reasonable database, and then eliminate the trend items and reduce the noise. In order to complete the subsequent data mining work efficiently, workstations were chosen to be used for data processing where the data was within the computer's reach and, therefore, no data imputation was required.
As the data samples have various characteristics and their extends are different from each other, the data needs to be standardized. There are many types of experimental data in this paper, such as the temperature data, vibration data, the surface roughness data and cutting force data. The data itself has a large numerical span, so the Zscore standardization is used. After the Zscore normalization process, the mean value of the data could become 0, while the standard deviation becomes 1. The calculation is as follows:
where x is the original data of the sample, is the mean of the original data, is the standard deviation of the original data, and x^{*} is the new sample data after the standardization process. Taking the cutting force data of the first group of experiments as an example, its cutting force data after standardization is shown in Fig. 11.
Fig. 11 Cutting force after standardized treatment
Cutting force after standardized treatment 
Considering that data preprocessing will make some features of the original data impossible to observe, the original data is qualitatively and quantitatively analyzed.
1) It is known from the temperature data that the influence of cutting width on cutting temperature in the cutting process parameters is more significant, and the cutting instability, severe roughness and machining surface hardening will make the overall temperature of the cutting process rise significantly.
2) By analyzing vibration data, it was found that higher amplitudes often mean higher temperatures. When the spindle speed reaches a certain value, the increase in cutting depth at this time will lead to increased vibration.
3) From the surface roughness data, it can be seen that if the spindle speed and feed speed can be reasonably increased, and the cutting depth can be reduced, it will effectively improve the quality of the processing surface and improve the processing efficiency, otherwise it will be counterproductive.
2 Cluster Analysis and Correlation Analysis of Milling Data
2.1 Cluster Analysis
The cutting force data has three characteristics: X, Y and Z direction components. In this paper, we propose to use the preprocessed dynamic milling force data as a classification, and use the kmeans algorithm as the core to filter and cluster the dynamic milling force and its corresponding cutting process parameters and time according to the principle of "the highest degree of similarity".
Taking the first group of experiments as an example, after the sample data is preprocessed and normalized, its distribution is shown in Fig. 12. When performing the cluster analysis, the best number of classifications for the group of samples was first determined by drawing a line graph, and then the iteration value was set. The value of k is usually picked up by a large number of experiments, d is the sum of the distances under different k values from a large number of experiments. The appropriate k value can be determined by making a line chart with d as the ordinate and k as the abscissa. The corresponding k value at the inflection point is the best choice. Here, the number of clustering centers k=4 and the number of iterations is 100. After 100 iterations, the clustering result suiting to the smallest selection of d is shown in Fig. 13, where different colors represent different categories.
Fig. 12 Distribution of experimental samples in Group 1
Distribution of experimental samples in Group 1 
Fig. 13 Clustering of Group 1 samples
Clustering of Group 1 samples 
The four data center coordinates are (0.692, 0.809, 0.210), (0.491, 0.142, 0.981), (0.866, 0.253, 1.340) and (0.730, 0.912, 0.115), and the sums of the distances from each centre to the data in the corresponding class are 666.203 5, 513.901 4, 548.703 8 and 635.161 0, respectively. As can be seen from Fig.14, the boundaries of the four data groups are obvious and the classification effect is good.
Fig. 14 Scatter plots of the correlation of the first set of experimental
Scatter plots of the correlation of the first set of experimental 
The remaining 45 groups of milling force data are clustered using the same method. Eventually it was found the majority (21 groups in total) of the experimental groups had milling force data that could be classified into 3 categories; 9 experimental groups had data that could be classified into 4 categories; 6 experimental groups had data that could be classified into 2 categories; few (2 groups in total) had data that could be classified into 5 categories; and the rest of the data was unclassifiable data that each belonged to a category. Cluster analysis of data, the sample with a large amount of data, is divided into different subcategories according to certain rules, so as to avoid problems such as large errors and insufficient accuracy when different types of data are analyzed together in the later stage.
2.2 Correlation Analysis
The clustering method described above classified the available 46 sets of milling force data into 126 categories, which shows that the main influencing factors for the variation of the milling forces are still different, even when they are generated under the same cutting process parameters. However, the clustering is of little significance for dynamic stiffness calculations at a later stage. In order to clarify the differences and characteristics of each type of data, further categorize them according to their main influencing factors, and improve the efficiency of data processing, the correlation analysis method was chosen to process each type of data and find their respective "labels".
From the above analysis it was found that temperature, roughness and vibration were the main influencing factors on dynamic milling forces, from which a preliminary judgement can be made that each group of experimental values of milling forces can be subdivided into multiple categories and their division may be related to these influencing factors. However, we know that dynamic milling forces are influenced by more factors than this. Therefore, in order to find a corresponding "label" for each group of data, a correlation analysis was chosen to correlate each category of data with temperature, roughness and vibration data, and to carry out a twotailed test for significance, while data that were not significantly correlated with these three factors were grouped into a new category "Other".
Before conducting the correlation analysis, a scatter plot was drawn to make a preliminary judgement on the linearity or nonlinearity of the relationship between the data, with a view to determining the type of correlation coefficient chosen. Taking the first six sets of data as an example, Fig. 14 shows the scatter plot of correlations corresponding to the first six sets of experiments.
As shown in Fig. 14, scatter plots of the correlation between the different classes of the first set of experimental threeway milling forces and the vibration data in the Z1 direction, the vibration data in the Y2 direction, the roughness data and the temperature data were plotted, respectively, and it was found that all relationships were not linear and the data itself did not obey a normal distribution, so the Spearman coefficient was chosen for further analysis.
Table 2 is an excerpt from the results of the correlation analysis performed on the different classes of the 46 groups of experimental milling forces (since the original data were divided into 126 classes of more categories, they are not listed here).
It is important to note that the correlations are grouped in the same category, but in different directions, e.g. two categories are significantly correlated with vibration data, but the first category is positively correlated and the second is negatively correlated.
The correlation analysis revealed that 40 classes were significantly correlated with vibration data, 72 classes were significantly correlated with temperature data, 67 classes were significantly correlated with roughness data, and only 16 classes were classified as correlated with other factors, some of which were significantly correlated with two or three kinds of vibration, temperature and roughness data at the same time. Therefore, the milling force data obtained from all experimental groups can be classified into eight major classes as a whole: 1) vibration as the main factor; 2) temperature as the main factor; 3) roughness as the main factor; 4) vibration and temperature as the main factors; 5) vibration and roughness as the main factors; 6) temperature and roughness as the main factors; 7) vibration, temperature and roughness all as the main factors; 8) vibration, temperature and roughness are not major factors. The data modeling should be performed later to match with them and be handled separately.
Results of correlation analysis
3 Data Modeling
As mentioned above, the dynamics of milling forces are divided into eight categories. The relationship between the dynamic milling force and corresponding independent variables was also different because of the different main categories of each type of data. According to the characteristics of different types of data, the dynamic milling force prediction model with high fitting accuracy is obtained by modeling the classified data with multielement regression analysis and RBF neural network, respectively.
3.1 Regression Analysis of Dynamic Cutting Force
Regression analysis is an important analysis tool for data mining, which has been widely used in modeling of various fields for its high efficiency, intuitionistic and accurate characteristics. In practical processing, the change of a certain index is caused by the coupling of several parameter variables, so multiple regression analysis is the most commonly used regression method in engineering application.
The current experimental formula of milling force is usually expressed in exponential form, which consists of correction factor, cutting process parameter and exponential factor^{[17]}. The available process parameters are cutting width, cutting depth, cutting speed, feed per tooth, feed, spindle speed, milling cutter diameter and number of teeth. However, there is an absolute correlation between the feed speed, the number of teeth per tooth and the spindle speed, and between the cutting speed, spindle speed and the diameter of the milling cutter.
With consideration of degree of close correlation between dynamic cutting forces and time, the time parameter is introduced into the model at the same time. The regression model for the dynamic cutting forces was initially determined to contain the following six parameters: depth of cut a_{p}, feed per tooth f_{z}, cutting width a_{e}, cutting speed v_{c}, number of teeth of the milling tool z and time t. The model expression is:
Compared with the general model of multilinear regression, the original model is in the form of polynomial summation, and the model is in the form of exponential class product, so the model is needed to deform, which are obtained from equations as follows:
Let y = lgF, b = lgAz, x_{1} = lga_{p}, x_{2} = lgf_{z}, x_{3} = lga_{e}, x_{4} = lgt, x_{5} = lgv_{c}, then the regression model is transformed into
The independent variable and dependent variable data were processed using the least squares. To find the constant term B of a polynomial, add a column of 1 to the data of the independent variable, and the first coefficient is the constant term.
Through correlation analysis, take the data of "vibration as the main factor" as an example, carry on logarithmic analysis to the data sample, then carry on the regression analysis. The results show that the judgment parameters of the regression model are not ideal and from the residual analysis graph it is proved that the confidence intervals of many residuals do not contain zeros.
(i) When the first regression analysis of the data is not satisfactory, the reasonable expectation of goodness of fit R2 is given first, and then the iteration is repeated. In the process of iteration, the original data without zero in all residual confidence intervals are eliminated, and the remaining data are retained to form new independent and dependent data samples. And then regression analysis is carried out again, until the goodness of fit is met.
(ii) At this point, according to the latest residual analysis image, adjust the fitting degree properly, and repeat the iterative process in step (1) until the evaluation index can judge the model more accurately.
Figure 15 shows the residual analysis of the first major category of data after the residuals are eliminated. The red lines in the two yellow circles in the graph indicate abnormal data, while the remaining green lines indicate the zero confidence interval from normal data to the residuals. The small circle in the centre of each line indicates the value of the residuals, while the whole line indicates the confidence interval of the residuals.
Fig. 15 Residual analysis graph
Residual analysis graph 
As shown in Fig. 15, after the exception data are culled, there are only 2 exception data left, which are negligible compared with the nearly 6 000 data samples. The analysis of the other scoring criteria showed that the R2 of fit was 0.776 6, the Fvalue was 3 989.293, and the probability of error P was 0. Therefore, it is clear that the constructed regression model is successful. The final regression analysis results for the eight categories of data are shown in Table 3.
After judging and analyzing the final results of the above eight types of multiple linear regression models, it is found that the first five dynamic cutting force multiple linear regression models, "vibration was the main factor", "vibration and roughness were the main factors", "vibration and temperature were the main factors", "vibration, roughness and temperature were the main factors", and "vibration, temperature and roughness were not the main factors", fit better, all at 0.6 F values are large, the probability of error is zero, and the number of confidence intervals for the residuals without zeros is less than 10, which can be directly applied. Their linear fit equations are:
1) Vibration is the main factor affecting the dynamic milling force variation:
2) Vibration and roughness are the main factors affecting the dynamic milling force variation:
3) Vibration and temperature are the main factors affecting the dynamic milling force variation:
4) Vibration, roughness and temperature are the main factors affecting the dynamic milling force variation:
5) Vibration, temperature and roughness are not the main factors influencing the dynamic milling force variation (i. e. "other"):
However, multiple linear regression models were constructed for the main factors of roughness, temperature, and temperature and roughness. Although the Fvalues, Pvalues and residual analysis were good, the fit was low at 0.266 5, 0.184 1 and 0.272 1, respectively, indicating that the three models were poorly fitted and not very accurate. Therefore, there was no significant linear or fitted linear relationship between the logarithms of the independent and dependent variables for the latter three types of data. By combining their scatter plots with the number of independent variables, it was found that the scatter plots did not fit the common nonlinear regression model. This result supports the conjecture that there may be a nonlinear relationship between the independent and dependent variables of these three types of data that cannot be expressed directly, or that there is no direct correlation between the three types of data. Therefore, the latter three types of data could be further analyzed by using RBF neural network modeling.
Results of regression analysis of milling force
3.2 RBF Neural Network Modeling of Dynamic Cutting Force
RBF neural network has a more intuitive structure, including the input layer, the implicit layer and the output layer. The input layer and the output layer respectively place the data samples of the independent variable and the dependent variable, and the number of input neurons is the number of independent variables. The implicit layer can be seen as a "bridge". When building a network, there are two transformations: the nonlinear transformation from the input layer to the implicit layer, and the linear transformation from the implicit layer to the output layer. It can be seen that the implicit layer plays a "bearing" role in the construction of the RBF neural network. In fact, when the vector is at a low latitude, there will be a linear inseparability, and the mapping of the vector from low latitude to high latitude can be realized through the implicit layer, so as to achieve the purpose of linear separability.
Assuming that there is a nonlinear relationship between the independent and dependent variables of the three types of data mentioned above and that the form of the relationship is complex and the deterministic relational equation can be obtained, we choose to use RBF neural networks for prediction and error analysis modelling. To build the model, a total of six neurons were set up in the input layer and lga_{p}, lgf_{z}, lga_{e}, lgv_{c}, lgz and lgt were entered into each neuron. The output layer has only one neuron and stores the milling force data in it. The three categories of data were divided into two separate parts. The last 200 data sets of each category were used as test samples and the rest of the data were used as training samples for the construction of the neural network.
To improve the accuracy of the neural network, relatively reasonable expansion coefficients were set for each of the three types of data. The schematic diagram of the neural network model is shown in Fig. 16. Finally, RBF neural network models were constructed for the three main factors of roughness (expansion coefficient of 11), temperature (expansion coefficient of 11) and temperature and roughness (expansion coefficient of 11) respectively. The fitted curves and error curves for the experimental values of the test samples are shown in Figs. 17 and 18 for the network model with roughness and temperature as the main factors. The red curve in Fig. 17 indicates the fitted value of the neural network, and the blue line curve indicates the experimental value, both of which have similar trends and relatively small differences. Combined with the error plot, it can be tentatively judged that the constructed model fits well and the error is relatively small and within the acceptable
Fig. 16 Schematic diagram of RBF neural network model
Schematic diagram of RBF neural network model 
Fig. 17 Fitting curve of RBF neural network
Fitting curve of RBF neural network 
Fig. 18 RBF neural network error profile
RBF neural network error profile 
range. To further judge the reasonableness of the model, mean squared error (known as MSE) is introduced, which was calculated as:
The MSE of the neural network model is 0.045 9 with temperature and roughness as influencing factors, 0.040 3 with temperature as influencing factor, and 0.040 4 with roughness as influencing factor. The MSE of the three models is very small, and the RNF neural network model has good prediction accuracy and can be applied in practice.
Comparing the modeling of the above regression analysis, it can be found that not all the relationships between dynamic milling force data and cutting process parameters can be obtained by multiple linear regression analysis, especially when considering the dynamic milling force variation at different times.
For the dynamic milling force data, where the vibration factor is not one of the main influencing factors, it is more meaningful to construct a model using RBF neural network. At present, the constructed model can directly predict the dynamic milling force by cutting process parameters and time, and then combine with the dynamic displacement change to obtain the dynamic stiffness of the milling tool system, but the specific expression of the dynamic stiffness of the milling tool system cannot be obtained.
4 Conclusion
A combination of data mining and cluster analysis with correlation analysis studies is conducted in this paper. The results show that the dynamic cutting force data is divided into eight main categories based on correlation as follows. "Vibration is the main factor", "vibration and temperature are the main factors", "vibration and roughness are the main factors", "vibration, temperature and roughness are the main factors", "vibration, temperature and roughness are not the main factors", "temperature is the main factor", "roughness is the main factor", "both temperature and roughness are the main factors.
Multilinear regression analysis was applied to model the data of "vibration as the main factor", "vibration and temperature as the main factor", "vibration and roughness as the main factor", "vibration, temperature and roughness as the main factor", "vibration, temperature and roughness are not the main factors", and fit the expressions to the five types of data. The three types of data, "temperature is the main factor", "roughness is the main factor", and "temperature and roughness are the main factors", were modeled by RBF neural network, and achieved high prediction accuracy.
Finally, with two prediction models, better predictions were obtained, providing a more accurate basis for selecting process parameters, ensuring production quality and reducing the occurrence of chatter.
References
 LvL S, DengZ H, LiuT, et al. Intelligent technology in grinding process driven by data: A review [J]. Journal of Manufacturing Processes, 2020, 58: 10391051. [CrossRef] [Google Scholar]
 KangH S, LeeJ Y, ChoiS, et al. Smart manufacturing: Past research, present findings, and future directions [J]. International Journal of Precision Engineering and ManufacturingGreen Technology, 2016, 3(1): 111128. [CrossRef] [Google Scholar]
 LiuX L, LiuQ, YueC X. Intelligent machining technology in cutting process [J]. Journal of Mechanical Engineering, 2018, 54(16): 4561(Ch). [Google Scholar]
 KantardzicM. Data Mining: Concepts, Models, Methods, and Algorithms [M]. Hoboke: John Wiley & Sons, Inc., 2011. [CrossRef] [Google Scholar]
 WangM, ZhuJ Y, ZhangY Z. Fuzzy pattern recognition of metal cutting state (II)—Fuzzy ISODATA clustering analysis method [J]. Journal of Aeronautics, 1985, 6(2): 173179(Ch). [Google Scholar]
 CaiY, ShaoH, YuanJ. [P]. CN:110428000A, 20191108. [Google Scholar]
 BrackeS, RadetzkyM, RosebrockC, et al. Efficiency and effectivity of high precision grinding manufacturing processes: An approach based on combined DEA and cluster analyses [J]. Procedia CIRP, 2019, 79: 292297. [CrossRef] [Google Scholar]
 PeiH J, ChenY Y, LiG A, et al. Correlation between milling force, vibration and surface roughness based on copula function [J]. Aeronautical Manufacturing Technology, 2019, 62(9): 5967(Ch). [Google Scholar]
 RodriguesJ R P, Gon?alvesR S, AraújoA SJr, et al. Correlation analysis of the roughness as function of cutting speed and feed rate (the roughness as function of machining parameters) [J]. Acta Scientiarum Technology, 2019, 41(1): 37228. [CrossRef] [Google Scholar]
 ThasthakeerA T, FaridA A, SengC T, et al. Analysis of the correlation between fractal structure of cutting force signal and surface roughness of machined workpiece in end milling operation [J]. Fractals, 2019, 27(2): 1950013. [CrossRef] [Google Scholar]
 LinY C, WuK D, ShihW C, et al. Prediction of surface roughness based on cutting parameters and machining vibration in end milling using regression method and artificial neural network[J]. Applied Sciences, 2020, 10(11): 3941. [Google Scholar]
 ShihabS K, KhanZ A, MohammadA, et al. RSM based study of cutting temperature during hard turning with multilayer coated carbide insert [J]. Procedia Materials Science, 2014, 6: 12331242. [CrossRef] [Google Scholar]
 RajeshK, Murali Krishnam RajuV V, RajeshS, et al. Effect of process parameters on machinability characteristics of CO_{2} laser process used for cutting SS304 stainless steels [J]. Materials Today: Proceedings, 2019, 18: 20652072. [CrossRef] [Google Scholar]
 YeganefarA, NiknamS A, AsadiR. The use of support vector machine, neural network, and regression analysis to predict and optimize surface roughness and cutting forces in milling [J]. The International Journal of Advanced Manufacturing Technology, 2019, 105(1): 951965. [CrossRef] [Google Scholar]
 MaherI, EltaibM E H, SarhanA A D, et al. Cutting forcebased adaptive neurofuzzy approach for accurate surface roughness prediction in end milling operation for intelligent machining [J]. The International Journal of Advanced Manufacturing Technology, 2015, 76(5): 14591467. [CrossRef] [Google Scholar]
 BasakS, DixitU S, DavimJ P. Application of radial basis function neural networks in optimization of hard turning of AISI D2 coldworked tool steel with a ceramic tool [J]. Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2007, 221(6): 987998. [CrossRef] [Google Scholar]
 WangP B. Radial Basis Function (RBF) neural network in gearbox mechanical fault diagnosis [C] // Abstracts of the 4th Advanced Materials, Mechanics and Structural Engineering Conference (AMMSE 2017). Beijing: IOP Publishing, 2017: 321325. [Google Scholar]
All Tables
All Figures
Fig. 1 Experimental scheme design idea
Experimental scheme design idea 

In the text 
Fig. 2 Workpiece after rough machining
Workpiece after rough machining 

In the text 
Fig. 3 Integral alloy end mill
Integral alloy end mill 

In the text 
Fig. 4 Stability lobe diagram
Stability lobe diagram 

In the text 
Fig. 5 Experimental equipment installation
Experimental equipment installation 

In the text 
Fig. 6 Roughness curve graph
Roughness curve graph 

In the text 
Fig. 7 Dynamic cutting force
Dynamic cutting force 

In the text 
Fig. 8 Infrared thermal imaging (a) and temperature change curve(b)
Infrared thermal imaging (a) and temperature change curve(b) 

In the text 
Fig. 9 Original vibration time domain graph
Original vibration time domain graph "Amplitude" is acceleration amplitude measured by accelerometer 

In the text 
Fig. 10 Surface roughness measuring instrument
Surface roughness measuring instrument 

In the text 
Fig. 11 Cutting force after standardized treatment
Cutting force after standardized treatment 

In the text 
Fig. 12 Distribution of experimental samples in Group 1
Distribution of experimental samples in Group 1 

In the text 
Fig. 13 Clustering of Group 1 samples
Clustering of Group 1 samples 

In the text 
Fig. 14 Scatter plots of the correlation of the first set of experimental
Scatter plots of the correlation of the first set of experimental 

In the text 
Fig. 15 Residual analysis graph
Residual analysis graph 

In the text 
Fig. 16 Schematic diagram of RBF neural network model
Schematic diagram of RBF neural network model 

In the text 
Fig. 17 Fitting curve of RBF neural network
Fitting curve of RBF neural network 

In the text 
Fig. 18 RBF neural network error profile
RBF neural network error profile 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.