Research Article
A feature selection method based on an improved fruit fly optimization algorithm in the process of numerical control milling Min Yuan
Advances in Mechanical Engineering 2018, Vol. 10(5) 1–10 Ó The Author(s) 2018 DOI: 10.1177/1687814018778227 journals.sagepub.com/home/ade
and Mei Wang
Abstract Automatic control is the key to improved production quality and efficiency of numerical control milling operations. Because the milling cutter is the most important tool in milling operations, the automatic monitoring of the tool wear state is of great significance. This work establishes a set of time domain and time-frequency domain features based on measurements of the cutting force for a computer numerical control milling machine and develops a method incorporating the Fisher criterion in an improved fruit fly optimization algorithm for selecting features most indicative of the tool wear state. A back propagation neural network was employed to test the effectiveness of the proposed feature selection method. Experimental comparisons with three other feature selection methods demonstrate that the proposed improved fruit fly optimization algorithm offers the advantages of the selection of a small number of significant features, easy implementation, precise optimization, rapid training, and good back propagation network performance. The proposed method has great potential for facilitating the practical monitoring of the milling tool wear state. Keywords Diagnostics, milling, numerical control machine tools, pattern recognition, tool wear
Date received: 28 March 2017; accepted: 18 April 2018 Handling Editor: Yangmin Li
Introduction In numerical control (NC) milling operations, the manufacturing quality of workpieces, such as surface roughness and dimensional accuracy, is largely influenced by the tool wear state. Moreover, most milling machine tool failures are generated from the tool system. Therefore, the development of online tool condition monitoring (TCM) is of great significance. TCM seeks to provide the essential information required for the predictive maintenance of milling cutters via signal processing methods and the extraction of the most significant signal features is indicative of tool condition (i.e. feature selection) from the time domain, frequency domain, and time-frequency domain. This represents
an extensive research focus worldwide. Feature selection, in particular, has raised considerable interest among researchers. Yum et al.1 adopted a new two-step combination feature selection method, which improved the performance of all classifiers. The maximum classification rate was 98.3%, which was an improvement of 4.2% compared with the best single-step feature School of Manufacturing Science and Engineering, Sichuan University, Chengdu, China Corresponding author: Mei Wang, School of Manufacturing Science and Engineering, Sichuan University, Chengdu 610065, Sichuan, China. Email:
[email protected]
Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License (http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).
2 selection method. Yu2 exploited the adaptive Gaussian mixture model (AGMM) for the effective assessment of tool wear. Good feature extraction results were achieved using Daubechies wavelet of order 5. Wang and Sun3 proposed a feature selection method based on ant colony optimization (ACO). This method, which modeled the feature selection process according to the behavior of ants searching for food, achieved good results. Lin et al.4 developed an efficient feature selection method by regression analysis. However, the experimental results represented a non-precise assessment because regression analysis cannot precisely describe non-linear relationships and decision variants for tool wear state monitoring. Zhao et al.5 used a clustering evaluation model to select significant features, although the model largely relied on professional knowledge, which is not suitable for a dynamic milling process. Goldberg6 employed genetic algorithms (GAs) in searching, optimization, and machine learning, which resulted in low efficiency owing to the absence of a common rule for parameter selection. To predict flank wear in drilling, Garg et al.7 employed particle swarm optimization (PSO) with a trained artificial neural network (ANN). Their experiments provided good prediction results in conjunction with rapid computation. In general, correlation analysis (i.e. the application of a correlation coefficient and fuzzy classification) is a good feature selection method with acceptable results.8 Moreover, intelligent optimization algorithms have attracted the attention of numerous scholars in recent years. These algorithms, such as the fruit fly optimization algorithm (FOA), ACO, and PSO, which are among the most effective, have become the most widely employed methods in optimization problems. In TCM, the signal acquisition methods employed to predict tool wear can be classified as direct and indirect methods. In the direct method, the actual tool wear is measured directly, which requires a temporary halt in manufacturing. Thus, most researchers have investigated indirect methods employing monitoring signals that can be obtained online, including the milling force, vibration, sound, acoustic emission (AE), temperature, spindle power, and surface roughness. In actual practice, a worn tool requires more force than a sharp tool to remove an equivalent amount of material, and the milling force is therefore considered one of the most effective parameters for monitoring tool wear based on previous experimental work.9 It has been observed that indirect measurements are subject to several experimental limits. For example, the computational effort involved in correlating process parameters with flank wear is high. Thus, significant effort has been devoted to improving computational models. Ren et al.10 unitized milling force measurements in a Takagi–Sugeno–Kang (TSK) fuzzy approach for TCM. However, it was observed that
Advances in Mechanical Engineering such models had difficulty estimating approximation errors and therefore required development to capture uncertainties during the turning process. Al-Habaibeh et al.11 and Fang et al.12 used multiscale methods, both of which were based on a milling force signal. Their research demonstrated the usefulness of the milling force, and the designed systems were able to predict tool wear successfully. However, the milling force is also sensitive to other parameters and can vary with cutting speed, depth of cut, and workpiece hardness, making correlation with wear more complicated. In this article, the time domain features of milling force signals are extracted by sensors. The features in the time-frequency domain are then extracted by wavelet analysis. Feature selection is conducted by an improved fruit fly optimization algorithm (IFOA), and the selected features are then input to a back propagation (BP) neural network to monitor the tool wear state. To verify the advantages of the proposed IFOA, experiments comparing four feature selection methods (the proposed IFOA, ACO, a correlation analysis method, and PSO) were conducted. The results verify that the proposed IFOA exhibits good adaptation and good optimization effectiveness. Therefore, this method is suitable for feature extraction in TCM.
Experimental milling cutter feature extraction scheme Experimental signal acquisition of milling force Milling force signals directly reflect the state of tool wear. Moreover, milling force signals respond quickly to changes in the state of tool wear, and are easily extracted, resulting in easily achieved TCM. Therefore, the milling force signal was selected as the feature extraction object and formed the basis for the judgment of the tool wear state. The overall work flow of the proposed TCM system is illustrated in Figure 1. The experiments were conducted on a Makino computer numerical control (CNC) machine equipped with an EGD 4440R milling cutter, A30N cutting material, and an ASSAB 718 HH machining workpiece with dimensions 206 mm 3 43 mm 3 106 mm, the ways of milling is face down milling. It has been shown that the milling force signal has no simple linear relationship to the tool wear state.13,14 Therefore, the milling force signal extracted by a single sensor cannot reflect the tool wear state accurately. A Kistler 9257B three-phase dynamometer and a Kistler 5019 multi-channel charge amplifier were employed to evaluate the milling force, and a NI-DAQ PCI1200 data acquisition board, Olympus microscope, and a Panasonic digital camera are employed to record the cutting process. Feature selection experiments were conducted in MATLAB R2013, equipped with a Windows 10 system, 64-bit
Yuan and Wang
3
Figure 1. Flow chart illustrating the tasks performed by a TCM system.
Figure 2. Schematic diagram of the experimental setup.
Figure 3. Images of the four different states representative of progressive tool wear.
operating system, processor frequency of 2.7 GHz, and 8 GB of memory. The cutting conditions tested involved spindle speeds of 600, 800, 1000, and 1200 r/min; feed rates of 100, 150, 200, and 300 mm/min; and a depth of cut of 1 mm. Figure 2 illustrates the experimental setup.
Figure 3 presents images of cutting tool wear at four stages: initial wear, normal wear, severe wear, and tool failure. Because tool wear mainly consists of various stages of flank wear (VB), the ISO 08688-1 standard defines 0.5 mm of VB as the tool failure limit.
4
Advances in Mechanical Engineering 2.
Mean value (N = 20) Fa (j) =
3.
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi u N u1 X F(j, i)2 Frms (j) = t N j=1
4.
ð3Þ
Fa(j) is the average milling force over Tj, which has a close relationship with tool wear, and has therefore been widely applied in TCM. F(j, i) represents the jth mechanical sampling in the ith circle. Root-mean-square value
Figure 4. Experimental tool wear process curves for a variety of tools.
Figure 4 presents the experimental tool wear curves of seven tools. Initial tool wear represents VB values of 0.00–0.10 mm, and the extent of VB increases rapidly in this stage. Normal tool wear represents VB values of 0.10–0.40 mm. In this stage, the wear resistance of the cutter increases, and the wear rate is not as rapid as that of the initial tool wear stage. Severe tool wear represents VB values of 0.40–0.50 mm. During this phase, the vibration resistance of the cutter is reduced, and the extent of VB increases rapidly. Finally, tool failure occurs when the extent of VB exceeds 0.50 mm.
N 1X F(j, i) N i=1
ð4Þ
Frms(j) represents the average energy of milling force signals. The total energy of the milling force increases with increasing tool wear, so Frms(j) can be used as a monitoring index reflecting the cutting tool state. Standard variance vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi u N u 1 X F(i, j) Fa (i) s(i) = t N 1 j=1
ð5Þ
The value of s (j) describes the quantitative fluctuation of signals over Tj and indicates the dynamic milling force. 5. Peak value
Extraction of time domain milling force features An end mill is employed as the research subject in the experiment, and six time-domain milling force features related to the Z direction were extracted through the dynamometer, including the maximum value (X1), peak amplitude (X2), given as the difference between the maximum and minimum values, mean value (X3), rootmean-square value (X4), standard variance (X5), and peak value (X6). Collection of the milling force signal (denoted as F(j, t)) was divided into 12 time periods Tj, for j = 1, 2, ..., 12, varying from the early wear to severe wear stages, and 20 data points i = 1, 2, ..., 20 were collected for each Tj. The time domain analysis and calculations are described as follows: 1.
Fpeak (j) =
N 1X f (i, j) N j=1
ð6Þ
Here, f (i, j) represents the ith peak value found in the sequence F(i, j), indicating the signal amplitude.
Extraction of time-frequency domain milling force features
fm (j)(X1 ) = max jF(j, i)j
ð1Þ
Wavelet analysis can decompose signals into independent frequency regions orthogonally without gaps or overlaps. These signals in the frequency domain are useful information for TCM. This study employed Daubechies wavelet of order 5 to decompose the milling force into four layers. Each node obtains a time-frequency feature Enj after wavelet analysis, and 16 features can be obtained after application of wavelet analysis to all four layers.
fa (j)(X2 ) = max½F(j, i) min½F(j, i)
ð2Þ
Theoretical background
Maximum value and peak amplitude of the milling force t2Tj
t2Tj
t2Tj
The two features, respectively, represent the steady state and transient state of the milling force.
Basic FOA The FOA is a global optimization searching algorithm proposed by Pan15 that models the food searching
Yuan and Wang
5
behavior of fruit flies. A schematic of the FOA is given in Figure 5, representative of the extraction of features X and Y. A fruit fly can discover far away food sources by gathering scent information over long distances due to its sensitive osphresis. Over shorter distances, a fruit fly can locate food sources and the clustering of other fruit flies using its acute vision. This search process is denoted as an interactive optimization evolutionary process. FOA has the advantages of rapid convergence rate and ease of implementation. Moreover, FOA performs well for solving mathematical extremum problems. Currently, the FOA has been applied to adjust financial warning models, locate mathematical extrema, and optimize the parameters of general regression neural networks and vector machines.16 Wu and Li17 compared the performance of FOA with five other evolutionary algorithms (GA, ACO, PSO, fish school algorithm, and immune algorithm) using the Schaffer formula. FOA was found to be superior to the other algorithms in terms of its reduced calculational burden. Moreover, FOA can optimize non-negative parameters easily. However, the disadvantage of FOA is equally obvious: low optimization precision owing to a tendency to converge to local optima. The basic steps involved in FOA include the following: 1. 2.
3.
4.
Randomly establish initial position of the fruit fly cluster: InitX ; InitY . Randomly establish food searching position (Xi, Yi) according to the osphresis of individual i: Xi = X + RaV ; Yi = Y + RaV , where RaV represents a random number sampled from a uniform distribution. Because the location of the food source is unknown, the distance to the origin (Disti) of the ith individual is first calculated as an estimate of the distance between it and the food source, which is then employed as the smell concentration judgment value (Si). Substitute Si into the smell concentration judgment function to determine the smell concentration of the ith individual (Smelli). The formula of Disti, Function Si, and Smelli are listed below Disti = = Xi2 + Yi2 Function Si = Disti1 Smelli = Function(Si )
5.
Determine the maximal Smelli ½bestSmell, bestindex = maxðSmellÞ
where bestSmell records the maximal Smelli and bestIndex records the value of the index i.
Figure 5. Illustration of the fruit fly optimization algorithm (FOA).
6.
Record the value of bestSmell as the present maximum Smellbest, record the coordinates corresponding to bestIndex, that is, (X1, Y1)bestIndex, and move the fruit fly cluster to this position by vision Smellbest = bestSmell X = X (bestindex); Y = Y (bestindex)
7.
Conduct an iterative optimization by repeating steps (2)–(5). If the new value of bestSmell is superior to Smellbest, implement step (6); otherwise, stop.
Preliminary testing using FOA demonstrated an unacceptable level of instability because tool wear is a random process, and premature convergence owing to its low optimization precision, making it unsuitable for the milling feature selection process of TCM. As such, the optimization ability of FOA requires improvement for application to the feature selection process. Numerous successful efforts have been made to improve the optimization performance of FOA. Marko et al.18 proposed chaotic FOA (CFOA) based on an investigation of FOA and another 10 different chaotic systems, and the method demonstrated a superior global optimal reliability and a high success rate. Wu et al.19 improved the FOA global searching ability by adjusting the entropy parameter to amplify the searching radius in a cloud model–based FOA. Wang20 optimized a wavelet neural network based on an IFOA to predict the melt index of industrial polypropylene. In the experiment, inertial weight parameters were employed to balance global and local searching abilities, resulting in an improved global searching ability. Wang21 added randomized mutation and group cooperation in FOA to optimize complicated functions and
6
Advances in Mechanical Engineering
Figure 6. Experimental tool wear process curves for a variety of tools.
solve joint replenishment problems (JRPs), resulting in an efficient FOA with an improved global searching ability. Experiments demonstrated the good performance of the IFOA.
Basic theory of Fisher screening Fisher linear discrimination analysis is one of the most effective methods of feature extraction. The Fisher discriminant achieves high discriminant efficiency because it can maximize and minimize different sample diversity. The Fisher criterion is given by the following expression22
dimensionality of features. Assuming the initial coefficient of each feature is 1, after FOA optimization, features with small coefficients are screened, and nearby fruit flies are clustered at the selected features. Then, the Fisher discrimination criterion is employed as a second optimization standard. Here, if the Fisher discriminant value of optimized features satisfies the output standard, the process is terminated; otherwise, we return to step (2) until the standard condition is fulfilled. A flow chart of the proposed IFOA is shown in Figure 6.
Experimental data analysis
T
X n SnX n Jn (X n ) = nT bn n , (X 6¼ 0) X Sw X
ð7Þ
Here, Sbn is the distribution matrix between classifications, Swn is the distribution matrix within a classification, and X n is the pattern vector in Rn , where Rn is the original feature space with n dimensions. Jn (X n ) represents the classification ability of an optimum classification identification vector X n . Rn + 1 is formed after a and related feature is added to Rn , n n+1 ). Thus, the classification ability is Jn (X ) = Jn + 1 (X not affected whether relevant features are added or deleted.
Improved fruit fly algorithm Because tool wear is a random process, a self-adapting FOA was adopted, where Fisher screening is added as a second decision after step (7), resulting in a reduced
Feature selection based on IFOA A total of 22 features were selected in the experiments, consisting of the 6 time domain features and 16 timefrequency domain features. The fixed factors are spindle speed of 800 r/min, feed of 150 mm/min, cutting depth of 1 mm, and sampling frequency of 2 kHZ. Authors set the initial coefficients of all 22 features as 1, which is then optimized by the proposed IFOA. The population of the fruit fly cluster is 22, and the maximum number of interactive optimization is 100. According to FOA theory, individuals with the largest Smell values fly in the same direction, which optimizes the vector coefficient. Features with small Fisher criterion values are removed, which completes the feature selection process for those features. In this article, training samples are divided into m and n classes, and the evaluation index is established by the Fisher method23 as follows. We assume that all feature sets
Yuan and Wang
7
Table 1. The results of feature selection by the proposed IFOA. Results of feature selection
Fisher discriminant values
Training time (s)
MSE
X2, X4, X5, X6, X11, X12, X17, X22 X1, X2, X6, X18, X19, X20 X1, X2, X4, X20, X22 X2, X4, X5, X10, X11, X12, X18, X19, X20, X22 All 22 features selected X2, X4, X5, X6, X12, X13, X14, X17, X18, X19, X22 X7, X8 X3, X4, X5, X6, X10, X11, X12, X13, X17, X18, X19, X20, X22 X2, X4, X6, X10, X11, X12, X13, X17, X18, X19, X22 X3, X10, X14, X17
1.7, 14.3, 5.5, 5.8, 2.2, 8.5, 1.6, 1.1 5.5, 1.7, 5.8, 4.3, 1.3, 3.1 5.5, 1.8, 14.3, 3.1, 1.1 1.7, 14.3, 5.5, 1.4, 2.2, 8.5, 4.3, 1.3, 3.2, 1.1
6.22 6.19 6.10 6.18
0.0831 0.0831 0.0831 0.0820
27.5 6.34
0.0611 0.0470
6.18 19.75
0.0382 0.0166
7.31
0.0156
6.29
0.0044
1.7, 14.3, 5.5, 5.8, 8.5, 13.6, 1.6, 1.6, 4.3, 1.3, 1.1 5.5, 1.3 13.6, 14.3, 5.5, 5.8, 1.4, 2.2, 8.5, 13.6, 1.6, 4.3, 1.3, 3.1, 1.1 1.7, 14.3, 5.8, 1.4, 2.2, 8.5, 13.6, 1.6, 4.3, 1.3, 1.1 13.6, 1.4, 1.6, 1.6
IFOA: improved fruit fly optimization algorithm; MSE: mean square error.
are ffj jj = 1, 2, :::g, and all features are ranked according to J as follows J ðfz1 Þ.J ðfz2 Þ. .J ðfzM Þ The selected features should clearly distinguish between the initial wear state and the severe wear state. Therefore, the evaluation index is defined by the Fisher criterion as follows J (fj ) =
ujm ujn 2 d2jm + d2jn
neural network, where the tool wear condition was obtained as the output of neurons. According to Kolmogorov’s theorem, the number of hidden neurons is 2N + 1. The hyperbola tangent sigmoid (sigmoid) function is suitable for application as the training function of a BP neural network. The logsine function was employed as the transfer function of output neurons, where (0,1) is the mapping range of neurons, which sufficiently satisfies the output requirements of the tool wear state.
ð8Þ
Here, ujm and ujn are the mean values of fj , and d2jm and d2jn are the variances of fj , where fj represents the corresponding feature in classes m and n. Defining classes m and n as the initial tool wear state and the tool failure state, respectively, J (fj ) reflects the distinguish ability of the two states of feature fj . As such, the greater the value of J (fj ), the stronger the ability to monitor the tool wear state. Therefore, the IFOA process involves calculating the Fisher criterion values of all features according to equation (8), ranking them from large to small, and then selecting the features with larger J (fj ) values to diagnose the tool wear state. Through trial and error, it can be concluded that the result of step (7) of the IFOA should be output only if J (fj ).1; otherwise, return to step (2) to choose again.
Diagnose tool wear state by BP neural network ANNs have demonstrated a strong capacity for fault identification.24 This study established a BP neural network with three layers to evaluate the effectiveness of feature selection. First, the data were normalized, and the selected cutting feature set was input into the BP
Experimental results Because IFOA is a random optimization searching algorithm for feature selection, the number of selected features varies after each optimization, but the good performance of the BP neural network is retained each time and is therefore suitable for monitoring the tool wear state after training. Table 1 reflects the relationship among the various feature subsets, the training time of the BP neural network, and the simulation errors (mean square error (MSE)), where the feature subsets have been listed according to decreasing MSE values. Clearly, the feature selection method based on IFOA demonstrated good performance. Compared to the results with all 22 features selected, the proposed method reduced the training time significantly and demonstrated low prediction error and fast training speed when the dimension was less than 10. Figure 7 illustrates the structure of BP neural network. N represents the number of input nodes (X3, X10, X14, X17), so the number of hidden layer is 9 and the final BP neural network is 4 3 9 3 1. Figure 8 shows the BP neural network performance of selected feature set (X3, X10, X14, X17), which provided the lowest MSE.
8
Advances in Mechanical Engineering foraging behavior that obtain high optimization efficiency and optimal positioning on the basis of evolutionary particle position changes. As for ACO and the correlation analysis method, the distinguish ability between features is not as definite as in the case of IFOA and PSO, and a greater number of features are therefore selected, possibly resulting in redundancy.
Figure 7. The structure of BP neural network.
Comparison experiments Comparison experiments were conducted to compare the results of feature selection methods employing the proposed IFOA, ACO, a correlation analysis method, and PSO on equivalent sets of the experimental milling force data. Feature selection results. The number of selected features directly affects the progress of TCM. As shown in Figure 9, the number of selected features obtained from equivalent data varied for each method. Obviously, IFOA achieved the best feature selection result, under the condition of equivalent selection effectiveness. While PSO selected a similar number of features as IFOA (i.e. 6), this is not unexpected because they are both evolutionary algorithms deduced from natural
Fisher discriminant values. The Fisher discriminant values of each selected feature were arranged from low to high for each of the selection methods considered, as is shown in Figure 10. In this experiment, the Fisher discriminant value represents the relationship between a feature and the tool wear condition, where the higher the value, the more closely the feature reflects the tool wear state. The feature sets selected by ACO and the correlation method both include X21, which has a Fisher discriminant value less than 1. While PSO and IFOA selected features with similar Fisher discriminant values, the average Fisher discriminant value obtained by PSO is 3.98, which is lower than that of IFOA with 4.55. Training time. Figure 11 presents the comparative results of computing time. It can be observed from the figure that ACO requires the longest computing time, which is greater than 140 s, while the other three feature selection methods have similar times of around 6 s and do not differ by more than 0.2 s. MSE. MSE is the most important index for the diagnostic effectiveness of the BP neural network. As shown in Figure 12. The following conclusions can be drawn from the comparison analysis. ACO selected the largest number
Figure 8. BP neural network performance of the selected feature set (X3, X10, X14, X17).
Yuan and Wang
9
Figure 11. Comparison of the training times required for various feature selection methods. Figure 9. Comparison of the number of selected features obtained from various feature selection methods.
Figure 12. Comparison of MSE values obtained for various feature selection methods.
Figure 10. Fisher discriminant values of the selected feature sets obtained from various feature selection methods.
of features, required the longest training time, and provided the highest MSE, which corresponds with the lowest optimization efficiency of all methods considered. The correlation method provided both a low training time and MSE, but it selected a relatively large number of potentially redundant features. While PSO performed well, the number of features selected, training time, and MSE were all a little greater than those obtained for IFOA. Therefore, we can conclude that IFOA provided the best optimization performance compared with the other three algorithms.
Conclusion An IFOA employing the Fisher criterion was developed and applied for feature selection in the process of
tool wear state monitoring for a CNC milling machine based on milling force data. The optimized set of features provided was demonstrated to provide for effective monitoring of the tool wear state. The following conclusions can be drawn from the experimental results: 1.
2.
The proposed IFOA realizes easy implementation, precise optimization, and rapid training by selecting a small number of significant features, resulting in good BP neural network performance. The proposed IFOA demonstrates good effectiveness and is suitable for use in TCM.
Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
10
Advances in Mechanical Engineering
Funding The author(s) received no financial support for the research, authorship, and/or publication of this article.
12.
ORCID iDs Mei Wang Min Yuan
https://orcid.org/0000-0002-4484-9714 https://orcid.org/0000-0003-4001-7651
References 1. Yum J, Kim TH, Kannatey-Asibu E, et al. A two-step feature selection method for monitoring tool wear and its application to the coroning process. Int J Adv Manuf Technol 2013; 64: 1355–1364. 2. Yu J. Machine tool condition monitoring based on an adaptive Gaussian mixture model. ASME J Manuf Sci Eng 2012; 134: 1–13. 3. Wang M and Sun QH. A feature selection method based on ant colony algorithm. J Vib Meas Diagn 2014; 34: 372–378. 4. Lin JT, Bhattacharyya D and Kecman V. Multiple regression and neural networks analyses in composites machining. Compos Sci Technol 2003; 63: 539–548. 5. Zhao DB, Song LL and Yan JH. Feature recognition method based on fuzzy clustering analysis and its application. Comput Integr Manuf 2009; 15: 2417–2423. 6. Goldberg DE. Genetic algorithms in search, optimization, and machine learning. Reading, MA: Addison Wesley, vol. 3, 1988, pp.95–99. 7. Garg S, Patra K and Pal SK. Particle swarm optimization of a neural network model in a machining process. Indian Acad Sci 2014; 39: 533–548. 8. Mei W. A study of tool wear condition identification and tool remain useful life prediction. Chengdu: Sichuan University, 2011. 9. Siddhpura A and Paurobally R. A review of flank wear prediction methods for tool condition monitoring in a turning process. Int J Adv Manuf Technol 2013; 65: 371–393. 10. Ren Q, Balazinski M, Baron L, et al. TSK fuzzy modeling for tool wear condition in turning processes: an experimental study. Eng Appl Artif Intel 2011; 24: 260–265. 11. Al-Habaibeh A, Al-Azmi A, Radwan N, et al. The application of force and acoustic emission sensors for detecting
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23. 24.
tool damage in turning processes. Key Eng Mater 2010; 419–420: 381–384. Fang N, Pai P and Mosquea S. Effect of tool edge wear on the cutting forces and vibrations in high-speed finish machining of Inconel 718: an experimental study and wavelet transform analysis. Int J Adv Manuf Technol 2011; 52: 65–77. Huang H and Li AP. Wear detection on cutting tools based on wavelet neural work. Trans Chin Soc Agric Mach 2008; 39: 173–177. Huang CJ. Using multi-stage data mining technique to build forecast model for Taiwan stocks. Neural Comput Appl 2012; 21: 2057–2063. Pan WT. A new fruit fly optimization algorithm: taking the financial distress model as an example. Knowl: Based Syst 2012; 26: 69–74. Han JY and Liu CZ. Fruit fly optimization algorithm with adaptive mutation. Appl Res Comput 2013; 30: 2641–2644. Wu XW and Li Q. Research of optimizing performance of fruit fly optimization algorithm and five kinds of intelligent algorithm. Fire Control Command Control 2013; 38: 17–25. Marko M, Vukovic´b N, Petrovic´a M, et al. Chaotic fruit fly optimization algorithm. Knowl: Based Syst 2015; 89: 446–458. Wu LH, Zuo CL and Zhang HQ. A cloud model based fruit fly optimization algorithm. Knowl: Based Syst 2015; 89: 603–617. Wang WC. Improved fruit fly optimization algorithm optimized wavelet neural network for statistical data modeling for industrial polypropylene melt index prediction. J Chemometr 2015; 29: 506–513. Wang L. An improved fruit fly optimization algorithm and its application to joint replenishment problems. Expert Syst Appl 2015; 42: 4310–4323. Lu L, Wang TY and Jiang YX. Tool wear state recognition based on hyper-sphere support vector machine. Trans Chin Soc Agric Mach 2011; 42: 218–222. Fukunaga K. Introduction to statistical pattern recognition. San Diego, CA: Academic Press Inc., 1992. Zhang J, Lu SL and He WX. Vibrating diagnosis of rolling bearings based on wavelet packet energy feature. Trans Chin Soc Agric Mach 2007; 38: 178–181.