Research Article  Open Access
Gerardo Beruvides, Fernando Castaño, Rodolfo E. Haber, Ramón Quiza, Alberto Villalonga, "Coping with Complexity When Predicting Surface Roughness in Milling Processes: Hybrid Incremental Model with Optimal Parametrization", Complexity, vol. 2017, Article ID 7317254, 11 pages, 2017. https://doi.org/10.1155/2017/7317254
Coping with Complexity When Predicting Surface Roughness in Milling Processes: Hybrid Incremental Model with Optimal Parametrization
Abstract
The complexity of machining processes relies on the inherent physical mechanisms governing these processes including nonlinear, emergent, and timevariant behavior. The measurement of surface roughness is a critical step done offline by expensive quality control procedures. The surface roughness prediction using an online efficient computational method is a difficult task due to the complexity of machining processes. The paradigm of hybrid incremental modeling makes it possible to address the complexity and nonlinear behavior of machining processes. Parametrization of models is, however, one bottleneck for full deployment of solutions, and the optimal setting of model parameters becomes an essential task. This paper presents a method based on simulated annealing for optimal parameters tuning of the hybrid incremental model. The hybrid incremental modeling plus simulated annealing is applied for predicting the surface roughness in milling processes. Two comparative studies to assess the accuracy and overall quality of the proposed strategy are carried out. The first comparative demonstrates that the proposed strategy is more accurate than theoretical, energybased, and Taguchi models for predicting surface roughness. The second study also corroborates that hybrid incremental model plus simulated annealing is better than a Bayesian network and a multilayer perceptron for correctly predicting the surface roughness.
1. Introduction
Despite the progress in science and technology, manufacturing processes are still considered complex systems consisting of several interacting subsystems with noisy, timevariant, and nonlinear behavior [1, 2]. New models are then required to better deal with manufacturing companies’ threats innovating in whole value chain [3]. Machining processes, and specifically milling processes, are key elements of manufacturing value chain [4]. The roots of the machining complexity rely on the inherent physical mechanisms governing the process itself. Different modeling techniques have been explored to describe and understand the geometric and physical characteristics of the cutting process [5, 6].
Quality monitoring is still the cornerstone and the bottleneck of many industrial processes [7]. For some manufacturing companies centered on aeronautical and automotive sectors, parts quality should fulfill strict requirements for dimensional [8, 9] and surface quality [10]. Inappropriate surface quality leads to removing defective parts from the production line and remanufacturing the component [11]. It implies more energy consumption of machine tools and more scraps resulting from the machining processes. Inprocess solution for predicting quality variables is in the agenda of many companies, and many solutions are being explored to decrease the environmental impact of machining processes due to poor surface quality of parts [12].
Surface roughness prediction studies in end milling operations are usually based on three main parameters composed of cutting speed, feed rate, and depth of cut [13]. Surface roughness and specific cutting energy consumption in slot milling are very related [14]. Recent studies demonstrate that net cutting specific energy is an effective process signature for surface integrity [15].
The review of traditional and emergent modeling strategies for manufacturing processes is out the scope of this paper. Nevertheless, it is important to remark that artificial neural networks (ANN) and fuzzy systems represent the most widely applied artificial intelligence techniques for modeling, control, and supervision [16, 17]. A thorough revision of computational intelligence techniques is out of the scope of this paper. The potential of artificial neural networks is often combined with the ability of fuzzy systems to represent human thought and robustness in the presence of noise and process uncertainty [18]. Fuzzy systems can also be combined with classical modeling techniques in order to better capture the most important characteristics of processes [19].
On the other hand, among modelfree methods reported in the literature, the incremental models are proposed in [20]. Incrementally, principle can be fully exploited starting from a simple and generic form of the model. Further refining procedures can be iteratively applied using some more localized techniques to represent specific regions of the input space. This work is focused on that principle for developing a hybrid incremental model inspired in [21] and then presents a case study for surface roughness prediction in milling processes. Surface roughness is influenced by tool wear, chattering, inappropriate cutting conditions, and so forth, and therefore it is mainly used to monitor the quality of the manufactured components [22].
In this work, a regression technique is applied for fitting the global characteristics whereas the fuzzy nearestneighbors algorithm [23] is selected for the subsequent refinement of the global model. The main advantages of fuzzy nearestneighbors algorithm have been already demonstrated in various studies [24, 25]. Noticeably, none of previous works has evaluated the capability of this instancebased learning method in surface roughness prediction.
However, setting parameters of the hybrid incremental model (HIM) is not an easy task because it strongly depends on the case study. Optimal tuning of parameters requires one strategy for optimization among all gradientfree techniques [26]. Simulated annealing (SA) is a probabilistic hillclimbing technique that emulates the physical annealing and cooling process of metals, based on a simple criterion that searches the problem space by piecewise perturbations of the estimates of the parameters that are being optimized [27]. SA has been used widely as optimization technique in many fields from optimal task allocation in manufacturing systems to vehicle routing systems [28].
From the best of authors’ knowledge, the main contributions of this work are twofold. Firstly, a modelbased procedure is implemented for predicting surface roughness based on a hybrid incremental modeling strategy. Secondly, a wellknown optimization method with proven convergence properties is applied to achieve the optimal parametrization for the HIM method. The performance of the proposed method is assessed by means of two comparative studies. Firstly, the prediction capability is analyzed with regard to a theoretical, energybased, and Taguchi models. Secondly, four errorbased criteria are applied to perform the comparison with a Bayesian network and a multilayer perceptron for predicting surface roughness. The suitability of the proposed approach for minimizing the number of parts to be remanufactured due to poor surface roughness is then demonstrated. The rest of the paper is organized as follows. Section 2 presents a description of hybrid incremental modeling. The experimental platform and the comparative study are presented in Section 3. Conclusions are shown in Section 4.
2. Hybrid Incremental Modeling Based on the Optimal Setting
The behavior of a locally nonlinear system [21] can be captured by hybrid incremental modeling where a global model represents the general behavior of the system and the local model shows the local behavior. Black box approach is an appropriate strategy when no prior knowledge of the system is available. For the sake of simplicity and clarity, a linear regression is selected for representing the global behavior of the system. Therefore, least squares are chosen to fit a polynomial of degree with the following output:where is the th input data and is its output value.
Therefore, the global model consists in the computed parameters of the fitting function (i.e., the degree polynomial).
The development of the local model requires a simple and easily interpretable technique from the viewpoint of industrial informatics. Data normalization [29] is carried out because different variable ranges affect negatively the performance of the algorithms. On the other hand, fuzzy Nearest Neighbors (FkNN) approach has the required characteristics justifying its selection [30]. The similarity between and the points of is given by where is the th neighbor of the query point and is the fuzzy strength parameter.
The target value of query point is now calculated as the mean of the target values of the points of the set , weighted by the similarity :
The incremental model integrates the global and local models described above. Thus, let be the function that is the output of the basic model. Then the prediction error of the basic model is
It is important to remark that FkNN only use errors resulting from the global modeling strategy. Therefore, the refinement of the global model output in regions with localized nonlinear behavior is done by FkNN.
The incremental model evaluates a sample of data input , by adding to the output of the basic model the compensation term calculated by the local model according to (4):
The tuning parameters of hybrid incremental model, namely, degree of the polynomial (), the neighborhood size (), and the fuzzy strength (), play an important role in setting the predictive model. These parameters influence the model’s performance, and they are required to be properly computed to produce the best output estimation using new input data.
Simulated annealing (SA) is selected among the arsenal of gradientfree techniques to search for an optimal set of model parameters. The main rationale for its selection relies on it easytoimplement probabilistic algorithm able to yield very good solutions for a wide variety of problems. A thorough study on the fast convergence and relatively low complexity of the algorithm is out of the scope of this paper but it has been recently reported [31]. The SA algorithm starts with initial modeling parameters and evaluates the MAPE performance index.
The current modeling parameters are perturbed to generate another and the mean absolute percentage error (MAPE) performance index is evaluated again. The Metropolis algorithm is the basement for the acceptance and rejection criterion. SA simulates the annealing process as it searches for a solution. A random perturbation is generated on the design variables , , which generates a change in the objective function (i.e., MAPE performance index). These perturbations depend on a temperature index, , and rate at which it is lowered (). More details on this temperature control parameter for simulated annealing are presented in [32].
Figure 1 represents the temperature index decreasing per iteration of the algorithm. Each set of model parameters obtained by this method is substituted into the hybrid incremental model (global plus local models) and the performance of the resulting solution is assessed by means of whole simulation of the system. Further details can be found in [33].
Finally, the MAPE performance index is evaluated by comparing the simulated responses with the desired responses and calculating the MAPE performance index. If the performance index is lower than the previous best performance index, then the new parameters replace the previous parameters. Otherwise, the new model parameters are not immediately discarded.
The probability, , of the new parameters’ cost () relative to the previous best cost () is calculated using Boltzmann’s equation:
The temperature index is then reduced by the annealing schedule using a reduction constant. The whole process is repeated until either the MAPE performance index has reached an acceptable minimum level or the temperature value has become too small to perturb the parameters. The main goal is to derive optimal modeling parameters that yield a fast and accurate model through the minimization of the MAPE performance index. Figure 2 shows the training procedure of the model using the target data. Figure 3 depicts how the resulting model estimates or predicts the corresponding output using the input data.
3. Experimental Study
Kondia HS1000 machining center equipped with a Siemens 840D openarchitecture CNC was used in experiments (see Figure 4). In all the cases, 170 × 100 × 25 aluminum AL7075T6 (UNS A97075) workpieces were used with hardness ranging from 65 to 152 Brinell. This material is commonly used in automotive and aeronautical applications.
Figure 5 illustrates the two geometry forms: (b) pockets and (d) island selected for experiments. For a better understanding, a constant spiral strategy was selected for the pocket elaboration and a Morph spiral strategy for the island form. Figure 5 also shows the simulation and cutting trajectories for the pocket (a) and the island (c), as well as, the final workpieces (b) pocket & (d) island after the milling process. In all the cases, the strategies, tool compensations, cutting regimens, and machining time were designed and simulated in ComputerAided Design and ComputerAided Manufacturing (CAD/CAM) software and, finally, exported to the CNC machine center.
(a) Pocket simulation and trajectory
(b) Pocket
(c) Island simulation and trajectory
(d) Island
Table 1 shows the six variables selected in the experimental setup. On one hand, feed rate (), tool diameter (), radial depth of cut (ae), and spindle rotation speed (rpm) were chosen for their influence on productivity, having direct relationship with key cutting variables such as cutting speed, cutting feed rate, and material removal rate, among others. On the other hand, the geometry curvature (Geom) and material hardness (HB) are much related to appropriate tool selection and the corresponding cutting strategies. For this reason, these six parameters were selected as key process indicators to model the milling process. Furthermore, the machining time for the different geometry is included in Table 1.

The most used index to characterize the surface roughness is the roughness average, Ra, that represents the arithmetic mean of the absolute ordinate values within a sampling length () as follows:
One interesting indicator related to surface roughness is the specific cutting energy consumption (SCEC) introduced by Liu et al. [14]. SCEC is analytically calculated based on the cutting parameters and toolworkpiece couple. SCEC is defined as the cutting energy consumed at the tool tip to remove 1 mm^{3} of material. Following the definition, the equation for a slotting milling process using a flatend mill can be estimated as follows:where is the average rotation power; Kte and Ktc are the tangential specific cutting and edge force coefficients; is the helix angle; is the number of flutes of the flatend mill; and are the lowest and highest position of the infinitesimal cutting edges engaged in cutting in axial direction on the th flute, respectively.
In order to match the model output with the actual industrial requirements and standards, roughness labels were also allocated according to the average value of surface roughness (μm) defined in ISO:1302 (2002) [34]. Table 2 shows variables and the respective assigned intervals for surface roughness from mirror to smooth.

Additionally, the following cutting parameters were also considered as inputs: spindle speed (rpm), feed rate (), and the resulting forces applied to all directions of cutting plane (). The cutting forces were measured using a Kistler 9257B dynamometer. The surface roughness Ra was measured with the Carl Zeiss Surfcom 130 stylus profilometer. Overall, 992 training samples were used in order to build the model, whereas 72958 test samples were used for the validation study. The average value of surface roughness, Ra, was chosen as the output variable of the model and it was the only variable offline measured.
Figure 6 represents the influence of resultant force, feed rate, and hardness on surface roughness for each tool diameter considered in the study. According to the depicted results, the richness of the training data guarantees that the four roughness levels are well taken into account to address the modeling task.
4. Experimental Results
4.1. Optimal Model Configuration
The hybrid incremental model plus simulated annealing procedure (HIM + SA) is implemented as follows. The SA algorithm parameters are reduction constant and initial temperature . Additionally, the number of proposed random perturbations at each temperature is 100 and the number of accepted solutions to proposed random perturbations is 30. The exploration space is , , where and .
In order to achieve an unbiased comparison for execution time of the algorithms, all algorithms run with a personal computer having Intel Core i53317U CPU 1.70 GHz 4 GB RAM. One hundred executions are run. The optimization process lasts 14 CPU seconds. The optimal parameters were obtained on the basis of the procedure described in Section 2 by simulated annealing optimization procedure (see Figure 3).
This means a secondorder polynomial (linear), two neighbors for the neighborhood size, and a value of 1.92 for the fuzzy strength parameter. The global model is obtained after the training. The errors (residuals) obtained in the training phase are recorded and then they are used during the evaluation by means of (4). The output of the model is obtained by applying (5). The minimum of MAPE performance index is 0.88%. The behavior of temperature () with regard to the evaluation number (ξ) to obtain the optimal solution is shown in Figure 7. It is clearly shown how after 10,000 evaluations the algorithm converges to a stable region.
Another interesting indicator is the behavior of cluster center position (Pos) with regard to the evaluation number shown in Figure 8. There are two main regions, one close to 0 value and the other close to 0.19. Furthermore, in some iterations, the cluster center is located at 0.08 or near to 0.
4.2. Discussion on Surface Roughness Prediction Problems
A comparison with three prediction models reported in the literature is performed to assess the accuracy of the proposed model. According to Wang and Chang [35], the theoretical Ra estimation for a slotted surface using a flatend mill can be calculated aswhere is the concavity angle ( by measurement).
The second model was extracted from the empirical representation surface roughness depending on specific cutting energy consumption (SCEC) index proposed by Liu et al. [14], represented bywhere and are empirical coefficients estimated in function of ap.
Finally, Taguchi methodology is a widely used approach to optimize the parameters for obtaining the minimum surface roughness, and the following representation is also taken from Liu et al. [14]: where is the predicted ratio for the predicted cutting condition.
For the sake of coherence with the abovementioned works, the performance index or figure of merit used to carry out a first comparison is the prediction relative error. Table 3 shows the prediction relative errors using models (see (10)–(12)) and HIM + SA.

This first comparative study aims at assessing the accuracy of the four models considered in the comparison. The column labelled Ra represents the nominal values measured in workpieces and the subscripts (κ, SCEC, η, HIM) represent the predicted values using the models (see (10)–(12)). The theoretical model, the SCEC model, and the Taguchi model are evaluated using the cutting conditions already reported in [14, 35] with , , . In the work herein reported the same material (AL7075) and the same tool (a 2flute solid carbide flatend mill with a diameter of 10 mm) are also used in the experiments. However, the depth of cut, the cutting speed, and the feed rate in our experiments (see Table 1) are higher than the cutting conditions reported in [14, 35]. In order to assess the precision of the four models in their respective cutting conditions, the predicted and measured values of surface roughness are compared. HIM + SA yields the best fitting between the predicted and measured values of surface roughness with best prediction errors.
Figure 9 depicts the graphical representation of measured versus predicted Ra for all models considered in this comparison. HIM + SA shows the best fit in all the cases, although it is important to highlight that SCEC model also shows very good behavior. From these results we can conclude that a clear representation of the milling process can be obtained with the proposed method, being a powerful tool to predict inprocess surface roughness based on the cutting parameters.
In the literature, there are plenty of artificial intelligence techniques applied to surface roughness detection and prediction. For the sake of clarity, two modeling strategies based on a Bayesian network (BN) and a multilayer perceptron (MLP) are also compared with the HIM + SA. In order to carry out a thorough evaluation of the model, four errorbased performance indices and two classification criteria are considered in this study. Therefore, in the second comparative study, a total of six performance indices were applied to assess the actual behavior on the basis of experimental run as follows: number of parts correctly manufactured (NPM), the number of parts to be remanufactured (NPR), the mean absolute error (MAE), the root mean squared error (RMSE), the relative absolute error (RAE), and the root relative squared error (RRSE).
Table 4 shows the main characteristics of the three approaches considered in the second comparative study. Multilayer perceptron (MLP) is one of the most widely applied neural networks at industrial level. However, the main drawback is related to the setting and tuning of network parameters such as number of hidden layers, number of nodes in the hidden layer(s), and form of activation functions. On the contrary, BNs have an easy and fast construction procedure without tuning parameters. The rationale for the good interpretability of BNs is indeed the reasoning based on realworld models.
 
BN: Bayesian network; ANN: artificial neural network; HIM + SA: hybrid incremental model with simulated annealing; TAN: tree augmented NaiveBayes; BP: backpropagation; MLP: multilayer perceptron; LSE: least square error; IBL: instance based learning. 
However, memory requirements might limit the transfers of knowledge and realtime applications. Further details on designing issues are reported in [36]. HIM + SA has some interesting features such as simple structure and easy training with few tuning parameters enabled by an optimal setting procedure (simulated annealing). HIM + SA also outperforms BN and MLP in terms of computing time. BN requires 0.08 CPU seconds, MLP requires 12.69 CPU seconds and HIM + SA requires 0.01 CPU seconds.
Table 5 shows the results of the second comparative study. The application of the BN and ANN yielded RAE of 10.41% and 13.05%, respectively. On the contrary, HIM + SA achieved an excellent accuracy, for instance, 2.68% in RAE, almost five times less error than the BN. The excellent behavior is also endorsed with the best number of parts correctly manufactured of 98.88%. Overall, HIM + SA outperforms significantly BN and ANN with regard to all figures of merits considered in this second study; specifically the number of parts to be remanufactured is the least value of all.

5. Conclusions
The development of adequate modeling strategies to deal with complexity of machining processes is a must. The prediction of surface quality is not straightforward. The need of decreasing parts to be remanufactured due to poor surface quality is motivated by the direct impact on economy and environment. This paper presents a hybrid incremental modeling strategy with an optimal setting procedure for quality detection in milling process. The procedure for designing and implementing hybrid incremental models with optimal parameters is simple and computationally efficient. The study demonstrates how surface roughness is predicted with a basic parameters configuration (secondorder polynomial, two neighbors, and the fuzzy strength parameter near to two) using a simulated annealing optimization method.
Two comparative studies with traditional techniques (i.e., theoretical model, energy model, and Taguchibased model) and artificial intelligencebased techniques (i.e., a multilayer perceptron and a Bayesian network) for predicting surface roughness have demonstrated that the proposed strategy outperforms significantly the others techniques and models considered in this study. The number of parts to be remanufactured using hybrid incremental model with optimal parametrization is less than applying either, a multilayer perceptron and a Bayesian network (e.g., about one percentage). The hybrid incremental model yields also better errorbased performance indices for predicting the surface roughness than abovementioned Bayesian and neural network models.
Hybrid incremental modeling plus simulated annealing for optimal parametrization can be extended to a wide range of manufacturing processes for estimating part quality. Further studies will be conducted to relate model outputs to other variables such as current and power consumption.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The work reported in this paper is partially supported by the European Project “IoSENSE: Flexible FE/BE Sensor Pilot Line for the Internet of Everything” funded by the Electronic Component Systems for European Leadership (ECSEL) Joint Undertaking under Grant Agreement no. 692480. The work is also partially supported by Project “AM4G: Advanced Manufacturing 4th Generation.”
References
 R. Bucki and P. Suchánek, “Modelling decisionmaking processes in the management support of the manufacturing element in the logistic supply chain,” Complexity, vol. 2017, Article ID 5286135, 2017. View at: Google Scholar  MathSciNet
 V. Modrak and Z. Soltysova, “Novel complexity indicator of manufacturing process chains and its relations to indirect complexity indicators,” Complexity, vol. 2017, Article ID 9102824, 2017. View at: Google Scholar  MathSciNet
 J. L. GarcíaAlcaraz, A. A. MaldonadoMacías, S. I. HernándezHernández, J. L. HernándezArellano, J. BlancoFernández, and J. C. Sáenz DíezMuro, “New product development and innovation in the maquiladora industry: A causal model,” Sustainability, vol. 8, no. 8, article no. 707, 2016. View at: Publisher Site  Google Scholar
 R. E. Haber, R. Haber, A. Alique, and S. Ros, “Application of knowledgebased systems for supervision and control of machining processes,” Handbook of software engineering and knowledge engineering, vol. 2, pp. 673–710, 2002. View at: Google Scholar
 G. Beruvides, R. Quiza, R. Del Toro, and R. E. Haber, “Sensoring systems and signal analysis to monitor tool wear in microdrilling operations on a sintered tungstencopper composite material,” Sensors and Actuators A: Physical, vol. 199, pp. 165–175, 2013. View at: Publisher Site  Google Scholar
 R. E. Haber and J. R. Alique, “Nonlinear internal model control using neural networks: An application for machining processes,” Neural Computing and Applications, vol. 13, no. 1, pp. 47–55, 2004. View at: Publisher Site  Google Scholar
 R.E. Precup, P. Angelov, B. S. J. Costa, and M. SayedMouchaweh, “An overview on fault diagnosis and natureinspired optimal control of industrial process applications,” Computers in Industry, vol. 74, pp. 1–16, 2015. View at: Publisher Site  Google Scholar
 R. Calvo, R. D’Amato, E. Gómez, and R. Domingo, “Integration of error compensation of coordinate measuring machines into feature measurement: Part II—experimental implementation,” Sensors, vol. 16, no. 10, article no. 1705, 2016. View at: Publisher Site  Google Scholar
 R. Calvo, R. D'Amato, E. Gómez, and R. Domingo, “Integration of error compensation of coordinate measuring machines into feature measurement: Part I—model development,” Sensors, vol. 16, no. 10, article no. 1610, 2016. View at: Publisher Site  Google Scholar
 G. Beruvides, F. Castaño, R. Quiza, and R. E. Haber, “Surface roughness modeling and optimization of tungstencopper alloys in micromilling processes,” Measurement: Journal of the International Measurement Confederation, vol. 86, pp. 246–252, 2016. View at: Publisher Site  Google Scholar
 P. B. Huang, “An intelligent neuralfuzzy model for an inprocess surface roughness monitoring system in end milling operations,” Journal of Intelligent Manufacturing, vol. 27, no. 3, pp. 689–700, 2016. View at: Publisher Site  Google Scholar
 M. Tomov, M. Kuzinovski, and P. Cichosz, “Modeling and prediction of surface roughness profile in longitudinal turning,” Journal of Manufacturing Processes, vol. 24, pp. 231–255, 2016. View at: Publisher Site  Google Scholar
 H. ElMounayri, Z. Dugla, and H. Deng, “Prediction of surface roughness in end milling using swarm intelligence,” in Proceedings of the 2003 IEEE Swarm Intelligence Symposium, SIS 2003, pp. 220–227, usa, April 2003. View at: Publisher Site  Google Scholar
 N. Liu, S. B. Wang, Y. F. Zhang, and W. F. Lu, “A novel approach to predicting surface roughness based on specific cutting energy consumption when slot milling Al7075,” International Journal of Mechanical Sciences, vol. 118, pp. 13–20, 2016. View at: Publisher Site  Google Scholar
 M. P. Sealy, Z. Y. Liu, Y. B. Guo, and Z. Q. Liu, “Energy based process signature for surface integrity in hard milling,” Journal of Materials Processing Technology, vol. 238, pp. 284–289, 2016. View at: Publisher Site  Google Scholar
 A. Alique, R. E. Haber, R. H. Haber, S. Ros, and C. Gonzalez, “Neural networkbased model for the prediction of cutting force in milling process. A progress study on a real case,” in Proceedings of the 2000 IEEE International Symposium on Intelligent Control, pp. 121–125. View at: Google Scholar
 R. E. Haber, J. R. Alique, A. Alique, J. Hernández, and R. UribeEtxebarria, “Embedded fuzzycontrol system for machining processes: Results of a case study,” Computers in Industry, vol. 50, no. 3, pp. 353–366, 2003. View at: Publisher Site  Google Scholar
 D. M. D’Addona, A. S. Ullah, and D. Matarazzo, “Toolwear prediction and patternrecognition using artificial neural network and DNAbased computing,” Journal of Intelligent Manufacturing, vol. 28, no. 6, pp. 1285–1301, 2015. View at: Publisher Site  Google Scholar
 J. M. Andújar and A. J. Barragán, “Hybridization of fuzzy systems for modeling and control,” RIAI  Revista Iberoamericana de Automatica e Informatica Industrial, vol. 11, no. 2, pp. 127–141, 2014. View at: Publisher Site  Google Scholar
 W. Pedrycz and K.C. Kwak, “The development of incremental models,” IEEE Transactions on Fuzzy Systems, vol. 15, no. 3, pp. 507–518, 2007. View at: Publisher Site  Google Scholar
 S.B. Roh, T.C. Ahn, and W. Pedrycz, “The refinement of models with the aid of the fuzzy knearest neighbors approach,” IEEE Transactions on Instrumentation and Measurement, vol. 59, no. 3, pp. 604–615, 2010. View at: Publisher Site  Google Scholar
 G. Kiswanto, D. L. Zariatin, and T. J. Ko, “The effect of spindle speed, feedrate and machining time to the surface roughness and burr formation of Aluminum Alloy 1100 in micromilling operation,” Journal of Manufacturing Processes, vol. 16, no. 4, pp. 435–450, 2014. View at: Publisher Site  Google Scholar
 C.B. Cheng and E. S. Lee, “Nonparametric fuzzy regression  kNN and kernel smoothing techniques,” Computers & Mathematics with Applications, vol. 38, no. 3, pp. 239–251, 1999. View at: Publisher Site  Google Scholar
 H.L. Chen, B. Yang, G. Wang et al., “A novel bankruptcy prediction model based on an adaptive fuzzy knearest neighbor method,” KnowledgeBased Systems, vol. 24, no. 8, pp. 1348–1359, 2011. View at: Publisher Site  Google Scholar
 Z. Ming, Z. Feng, and J. Zhengbiao, “A study on speech emotion recognition based on fuzzy K nearest neighbor,” International Journal of Multimedia and Ubiquitous Engineering, vol. 10, no. 10, pp. 57–66, 2015. View at: Publisher Site  Google Scholar
 A. Fiannaca, G. Di Fatta, R. Rizzo, A. Urso, and S. Gaglio, “Simulated annealing technique for fast learning of SOM networks,” Neural Computing and Applications, vol. 22, no. 5, pp. 889–899, 2013. View at: Publisher Site  Google Scholar
 R. A. Rutenbar, “Simulated annealing algorithms: an overview,” IEEE Circuits and Devices Magazine, vol. 5, no. 1, pp. 19–26, 1989 (English). View at: Publisher Site  Google Scholar
 F. Zhao, Y. Hong, D. Yu, Y. Yang, Q. Zhang, and H. Yi, “A hybrid algorithm based on particle swarm optimization and simulated annealing to holon task allocation for holonic manufacturing system,” The International Journal of Advanced Manufacturing Technology, vol. 32, no. 910, pp. 1021–1032, 2007. View at: Publisher Site  Google Scholar
 F. Penedo, R. E. Haber, A. Gajate, and R. M. Del Toro, “Hybrid incremental modeling based on least squares and fuzzy KNN for monitoring tool wear in turning processes,” IEEE Transactions on Industrial Informatics, vol. 8, no. 4, pp. 811–818, 2012. View at: Publisher Site  Google Scholar
 J. Derrac, F. Chiclana, S. García, and F. Herrera, “Evolutionary fuzzy knearest neighbors algorithm using intervalvalued fuzzy sets,” Information Sciences, vol. 329, pp. 144–163, 2016. View at: Publisher Site  Google Scholar
 M. Almaraashi, R. John, A. Hopgood, and S. Ahmadi, “Learning of interval and general type2 fuzzy logic systems using simulated annealing: Theory and practice,” Information Sciences, vol. 360, pp. 21–42, 2016. View at: Publisher Site  Google Scholar
 D.Y. Sun and P.M. Lin, “The solution of time optimal control problems by simulated annealing,” Journal of Chemical Engineering of Japan, vol. 39, no. 7, pp. 753–766, 2006. View at: Publisher Site  Google Scholar
 F. Castaño, R. E. Haber, R. M. Del Toro, and G. Beruvides, “Application of hybrid incremental modeling for predicting surface roughness in micromachining processes,” in Proceedings of the 2014 IEEE Symposium on Computational Intelligence for Engineering Solutions, CIES 2014, pp. 54–59, USA, December 2014. View at: Publisher Site  Google Scholar
 Geometrical Product Specifications (GPS)—Indication of surface texture in technical product documentation, 2002.
 M.Y. Wang and H.Y. Chang, “Experimental study of surface roughness in slot end milling AL2014T6,” The International Journal of Machine Tools and Manufacture, vol. 44, no. 1, pp. 51–57, 2004. View at: Publisher Site  Google Scholar
 M. Correa, C. Bielza, and J. PamiesTeixeira, “Comparison of bayesian networks and artificial neural networks for quality detection in a machining process,” Expert Systems with Applications, vol. 36, no. 3, pp. 7270–7279, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Gerardo Beruvides et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.