## Articulo

• Citado por SciELO
• Accesos

• Citado por Google
• Similares en SciELO
• Similares en Google

## versión On-line ISSN 0718-221X

### Maderas, Cienc. tecnol. vol.21 no.4 Concepción dic. 2019

#### http://dx.doi.org/10.4067/S0718-221X2019005000406

ARTÍCULOS

Optimization of wood machining parameters in CNC routers: Taguchi orthogonal array based simulated angling algorithm

1Faculty of Forestry, Department of Forest Industry Engineering, Istanbul University-Cerrahpasa, Istanbul, Turkey

2Faculty of Forestry, Department of Forest Industry Engineering, Istanbul University-Cerrahpasa, Istanbul, Turkey. hkoc@istanbul.edu.tr

Abstract:

In the present study, two mathematical models were developed to optimize the surface roughness for machining condition of Cedar of Lebanon pine (Cedrus libani). Taguchi approach was applied to examine the effect of CNC processing variables. Quality characteristics parameters were selected as arithmetic average roughness (R a) and average maximum height of the profile (R z ) for wood material. Analysis of variance (ANOVA) was used to determine effective machining parameters. Developed mathematical models using response surface methodology (RSM) were optimized by a combined approach of the Taguchi’s L 27 orthogonal array based simulated angling algorithm (SA). Optimum machining levels for determining the minimum surface roughness values were carried out three stages. Firstly, the desirability function was used to optimize the mathematical models. Secondly, the results obtained from the desirability function were selected as the initial point for the simulated angling algorithm. Finally, the optimum parameter values were obtained by using simulated angling algorithm. Minimum R a value was obtained spindle speed of 17377 rpm, feed rate of 2,012 m/min, tool radius of 8 mm and depth of cut of 2,009 mm by using desirability function based simulated angling algorithm. For R z these results were found as 16980 rpm, 2,004 m/min, 8,001 mm and 2,003 mm. The R-square values of the R a and R z were 95,91 % and 96,12 %, respectively. The proposed models obtained the minimum surface roughness values and provided better results than the observed values.

Keywords: Cedrus libani; response surface method; softwood; surface roughness; wood material; optimization

Introduction

Wood surface roughness is a crucial indicator of the quality of CNC processing parameters. However, wood material has very complex structure parameters such as machining procedure, physical properties of panel, processing conditions and anatomical structure (Magoss 2008, Philbin and Gordon 2006, Hiziroglu and Kosonkorn 2006, Ozdemir and Hiziroglu 2007, Ratnasingman and Scholz 2006, Hazir and Ozcan 2019). For this reason, it is difficult and complex to investigate the optimum parameter levels. A systematic approach could be applying to provide the optimum process conditions. Therefore, a novel experimental design is required to reduce the number of experiment. Traditional experimental design was required have a long time and it increased the cost. Design of experiment (DOE) and the Taguchi methodology are powerful techniques to reduce the number of experiments. Moreover, these methods have been commonly applied in different engineering applications (Yang and Tarng 1998, Taguchi et al. 2005, Sarikaya and Güllü 2016, Rao and Murthy 2018, Selaimia et al. 2017, Azhiri et al. 2014, Majumder et al. 2017, Kant and Sangwan 2014, Deepanraj et al. 2017). In the recent years, metaheuristic algorithms have been applied to various engineering problems. Especially, experimental design based metaheuristic algorithms have become popular method for optimizing the parameters (Rao and Kalyankar 2013, Mahes et al. 2015, Samanta 2009). Coelho et al. (2008) applied the Taguchi technique to evaluate the effective parameters for wood cutting conditions. Gaitonde et al. (2008) studied the Taguchi design to optimize drilling parameters such as cutting speed and feed rate to reduce delamination of medium density fiberboard (MDF). The delamination was reduced higher cutting speed and lower feed rate. Davim et al. (2009) studied the determination of cutting factors such as cutting speed and feed rate on surface roughness in MDF milling. According to the results, the wood surface roughness decreased with an increase of spindle speed and increased with the feed rate. Prakash and Palanikumar (2010) applied the Taguchi approach to find the effective ratio on R a for drilling MDF. The drilling parameters were spindle speed, feed rate and drill diameter. According to the Taguchi analysis, the feed rate was found to be the most important factor affecting surface roughness. Sofuoglu (2017) investigated that the optimization of CNC machining parameters were examined using the Taguchi method on the surface quality of massive wooden panels made of Scots pine (Pinus slyvestris L.). Optimal cutting results of R a and R z were resulted with cutter 1, at a tool clearance strategy of a raster spindle speed of 16000 rpm, feed rate of 1000 mm/min and depth of cut of 4 mm. Hazir et al. (2018) used 2 k full factorial design to determine the significant variables and optimum processing levels to determine the lower the wood surface roughness. As a result of this study, two equations were created for radial and tangential cutting directions. Moreover, each models were solved by applying desirability function. Koç et al. (2017) applied full factorial design to find the significant variables on wood surface roughness. Surface roughness values were determined by using laser integrated measurement system and stylus type equipment. According to this study, the roughness values for these measurement methods were found similar. Wilkowski et al. (2011) used the Taguchi technique to determine CNC processing factors for wood surface roughness. The parameters such as feed rate and spindle speed were investigated. Asiltürk et al. (2016) used the Taguchi orthogonal array and RSM to determine the effective factor. Moreover, the effective parameters were optimized for CNC lathe machining. Zhou et al. (2017) applied multi objective optimization technique by using grey relation analysis (GRA), radial basis function (RBF) neural network and particle swarm algorithm. In this study, workpiece material was selected as Ni-based superalloy Inconel 718. According to results, presented method was determined as successfully for multi-axis ball-end milling process. Zhou et al. (2016) used response surface methodology and genetic algorithm to optimize the surface roughness of nickel-based single crystal superalloy using microgrinding process. According to verification results, optimum parameter levels were found as reliable. Bharathi and Baskar (2011) proposed the particle swarm optimization (PSO) algorithm to analyze the optimal processing factors for minimizing the surface roughness and processing time using aluminum material. Parameters of cutting speed, feed, and depth of cut were optimized by using PSO. They indicated that presented method was efficient for machining process. Bao et al. (2018) predicted wood surface roughness conducting the neural network for wood sanding process. They reported that presented model was usefully predicting the roughness value. According to previous studies, the effective parameters on the wood surface roughness are defined in cause and effect diagram given in (Figure 1). Figure 1: Cause and effect diagram of wood surface roughness parameters.

CNC processing parameters and their effects on wood surface roughness were determined in various past studies. However, there is little as no information evaluating the optimum cutting parameters for woodworking industries. Most of these studies presented various engineering problems not inherited to wood based products. Therefore, this study aims to present an approach for determining minimum surface roughness by integrating RSM, DF and SA in CNC processing.

Materials and methods

Cedar of Lebanon pine (Cedrus libani A. Rich) species with intensive use in the furniture industry was selected for the study. The samples were prepared with the dimension of 140 mm x 50 mm x 18 mm for each procedure. Density level of Cedar of Lebanon pine was measured randomly through 27 samples. Samples were conditioned in a climate room having a temperature of 20°C and relative humidity of 65% until they reach a moisture content of 9±1%. The value of wood density was 510 kg/m³. These samples were processed with 3-axis CNC. Alpha-CAM program was used for determination of the tool path for the sample. Roughness measurement device is a stylus-based portable profilometer that is Sutronic-25 type equipment. Parameters of Ra and Rz were used to evaluate the wood surface quality (ISO 4287-1997).

Statistical design of experiment

The statistical analysis was performed using Minitab software package 17. The Taguchi L27 orthogonal statistical designs were followed to optimize the experimental procedure. The processing variables were screened by ANOVA analysis. These factors were modeled by using response surface methodology (RSM). The results were analyzed by multiple regression analysis through least square method. All the terms of model were tested and verified statistically by F-test at probability levels (p<0,05). R-square and Adj-R-square were used to perform the developed the models. After the fitting models, 3D surface plots were applied to evaluate the process parameters. Finally, desirability function based simulated angling algorithm was used to screen the optimum parameter levels. The experimental produce is shown in (Figure 2).

Response surface method (RSM)

Response surface method (RSM) is a statistical based technique. RSM method involves three stages. The first stage is the selection of independent variables of major effects on the process through screening studies and the delimitation of the experimental region. The second stage is the mathematic-statistical apply of the obtained experimental data through the fit of a polynomial function and evaluation of the model’s fitness. The third stage is the verification of the necessity and possibility of performing a displacement in direction to the optimal region and selection of optimum values of parameters. The general form second order polynomial equation is given Equation 1:

According to Equation 1 Y is the predicted reaction or reactions (R a and R z ), X i and X j are variables, β0 a constant, βi, βii and βij of linear, quadratic and the second-order terms, respectively and k is the number of independent parameters (k=4 in this study) and ϵ is the error term.

Experimental design

The traditional experimental design is time consuming and costly. L 27 orthogonal design was used to overcome this issue. Especially, this method is used by manufacturing process and engineering analysis. These arrays of different combination for variables were applied by the Taguchi approach to analyze the experimental study. While 81 experiments are required to investigate all the effects of the parameters, 27 experiments are sufficient thanks to Taguchi orthogonal array. According to the Taguchi orthogonal array, the minimum number of experiment and coded variables were computed by using Equation 2 and Equation 3.

The Taguchi experiments trial number:

The terms of L and P are displayed number of levels and number of parameters.

In terms of X i, X o and Δx are coded variable, actual variable, center palace, and altered variable. (Table 1) shows the number of parameters and their levels.

In this study, four parameters namely spindle speed, feed rate, depth of cut and tool radius were selected as input parameters while R a and R z were selected as output variables. (Table 2) shows the experimental parameters and their results.

Results and discussion

In this section, two mathematical models used for optimization of the surface roughness were explained. The Taguchi orthogonal array integrated with response surface methodology and simulated angling algorithm for solving these models were detailed. In addition to these, the suggested models were tested by using verification test. Interactions between the variables were analyzed by using 3D surface plot. Finally, the results were compared and the suggested models were discussed.

Developing the mathematical models

The results were analyzed by conducting Minitab 17 software. Correlation of the experimental roughness value and the estimated values from the regression equation were shown in (Table 3). The experimental data were analyzed by the polynomial models such as linear, 2-way interaction and full quadratic models. The adequacy of models performed indicated that linear, interaction and quadratic models had lower p-value (<0,05). (Table 3) displays the R² and Adjusted-R² values were found low in linear and interaction models. R-square and Adj-R-square values were described with Multiple Linear Regression Analysis (MLR). Quadratic models were calculated to have maximum R² and Adj- R² values. Hence the quadratic models were chosen for R a and R z in this study.

Analysis of variance (ANOVA)

Table 4 and Table 5 indicated variance analysis for R a and R z . In this case f, n, d, r, n², d², fd, fn, nd, nr, fr are significant factors. Contribution (PC %) of each factor on the total variation is computed as percentage. The most effective parameters for R a and R z were interaction between feed rate and spindle speed, their contribution to the model values were 37,16 % and 34,47 %, respectively, on the surface roughness. The coefficient of determination R² for the R a and R z of were 95,91 % and 96,12 % respectively. The Adjusted R 2 values for R a and R z were computed as 91,15 % and 95,93 % respectively.

Table 4: ANOVA results for R a . The terms of (ª) indicates significant parameter.

Table 5: ANOVA results for R z . The terms of (ª) indicates significant parameter.

Evaluation of the models

Normal Probability Plot (NPP) was applied to evaluate the data for normality. Residual is the mean difference between the observed value and the predicted or fitted value. If the residuals fall approximately along straight line, they are then normally distributed. In contrast, if the residuals do not fall fairly close to a straight line, they are then not normally distributed (Antony 2014). If the model is correct and if the assumptions are satisfied, the residuals should be structures; in particular, they should be unrelated to any other variable including the predicted response. Figure 3a, Figure 3b, Figure 4a and Figure 4b show that the residuals generally fall on a straight line implying that the errors were disturbed normally, meaning was the experimental data come from a normal population. As a result of the residuals, no unusual structures were apparent. Figure 3: (a) Fitted value for R a and (b) NPP of residuals for R a . Figure 4: (a) Fitted value for R z and (b) NPP of residuals for R z

Surface graphs and analysis

Response surface provides the relationship between the factors with each other and the response variable. (Figure 5a) shows the surface plot for R a in CNC machining of wood material with three levels of feed rate and three levels of spindle speed for mid-level hold value of depth of cut with 4 mm and tool radius of 10 mm. The surface roughness increased with decreases the spindle speed and increases the feed rate. (Figure 5b) shows the surface plot for R a in CNC machining of wood material with three levels of depth of cut and three levels of tool radius for mid-level hold value with feed rate of 5 m/min and spindle speed of 15000 rpm. The surface roughness decreased with a decrease of tool radius and depth of cut. (Figure 5c) shows the surface plot for R a in CNC machining of wood material with three levels of feed rate and three levels of tool radius for mid-level hold value with spindle speed of 15000 rpm and 4 mm depth of cut. The surface roughness increased with the increase in tool radius and feed rate. (Figure 5d) shows the surface plot for R a in CNC machining of wood material with three levels of feed rate and three levels of depth of cut for mid-level hold value with spindle speed of 15000 rpm and tool radius of 10 mm. The surface roughness increased with the increase in feed rate and depth of cut. (Figure 5e) shows the surface plot for R a in CNC machining of wood material with three levels of spindle speed and three levels of tool radius for mid-level hold value with 5 m/min of feed rate and depth of cut of 4 mm. The surface roughness increased with the increase in tool radius and decreases with spindle speed. (Figure 5f) shows the surface plot for R a in CNC machining of wood material with three levels of depth of cut and three levels of spindle speed for mid-level hold value with 5 m/min of feed rate and tool radius of 10 mm. The surface roughness increased with the increase in depth of cut and decrease with spindle speed. Figure 5: According to 3D surface plots (a, b, c, d, e, f) of different values processing parameters for R a .

(Figure 6a) shows the surface plot for R z in CNC machining of wood material with three levels of feed rate and three levels of spindle speed for mid-level hold value of depth of cut with 4 mm and tool radius of 10 mm. The surface roughness increased with decreases the spindle speed and increases the feed rate. (Figure 6b) shows the surface plot for R a in CNC machining of wood material with three levels of depth of cut and three levels of tool radius mid-level hold value with feed rate of 5 m/min and spindle speed of 15000 rpm. The surface roughness decreased with a decrease of tool radius and depth of cut. (Figure 6c) shows the surface plot for R a CNC machining of wood material with three levels of feed rate and three levels of tool radius for mid-level hold value with 15000 rpm of spindle speed and depth of cut of 4 mm. The surface roughness increased with the increase in tool radius and feed rate. (Figure 6d) shows the surface plot for R a in CNC machining of wood material with three levels of feed rate and three levels of depth of cut for mid-level hold value with 15000 rpm of spindle speed and tool radius of 10 mm. The surface roughness increased with the increase in feed rate and depth of cut. (Figure 6e) shows the surface plot for R a in CNC machining of wood material with three levels of spindle speed and three levels of tool radius for mid-level hold value with 5 m/min of feed rate and depth of cut of 4 mm. The surface roughness increased with the increase in tool radius and decreases with spindle speed. (Figure 6f) shows the surface plot for R a in CNC machining of wood material with three levels of depth of cut and three levels of spindle speed for mid-level hold value with 5 m/min of feed rate and tool radius of 10 mm. The surface roughness increased with the increase in depth of cut and decrease with spindle speed. Figure 6: According to 3D surface plots (a, b, c, d, e, f) of different values processing parameters for R z

Parameter optimization

The aim of present work was to find the optimal processing levels leading to minimum values of R a . and R z . Two second order mathematical models were created by using RSM. Moreover, backward elimination method was carried out for determining the non-significant terms in mathematical models. To formulate the optimization problem, proposed models were presented in Equation 4 and Equation 5:

These models were solved with two different methods namely, desirability function and simulated angling algorithm. Firstly, the desirability function was applied to determine the initial optimum parameter levels. Then, the results were used the initial point for the simulated angling algorithm. The functions terms were optimized within the specified range. The range of machining conditions was selected and given as follows Equation 6a, Equation 6b, Equation 6c, Equation 6d:

Within ranges of processing parameters,

The limitations of cutting variables value of spindle speed, feed rate, depth of cut and tool radius were 12000 rpm, 2 m /min, 2 mm and 8 mm as lower limits while 18000 rpm, 8 m/min, 6 mm and 12 mm were selected as upper limits, respectively.

Parameter optimization by using desirability function (DF)

This is commonly applied engineering problems for the optimization process. The function value is range 0 from 1. In the present work, selected quality characteristic was smaller-the-better for minimizing the surface roughness value. This equation was presented in Equation 7:

Where T symbolizes the target value of the i th response, y i , L symbolizes the acceptable lower limit value, U symbolizes the acceptable upper limit, for this response and W represents the weight. In (Table 6), optimal machining parameters for minimizing R a was found with spindle speed of 17450 rpm, feed rate of 2,9 m/min, tool radius of 8,0 mm and depth of cut of 2,0 mm. optimal machining parameters for minimizing R z was spindle speed of 16990 rpm, feed rate of 2,0 m/min, tool radius of 8,0 mm and depth of cut of 2,0 mm.

From the Figure 7 and Figure 8, the optimum values for R a and R z were 3,1287 µm and 11,254 µm, respectively.

Parameter optimization by simulated angling algorithm (SA)

In this study, simulated angling algorithm was applied to optimize the machining parameters. It is the one of the most useful metaheuristics method for solving the engineering problems due to easy application. It is a widely used method in solving optimization problems. SA is a probabilistic and single solving method based on the annealing process of metallurgy. The algorithm starts with an initial point using random solution, S, then a simple random change is produced to actual solution (S’). Obtained solution was compared with new solution by creating the objective function. The new solution is accepted when the objective f(S’) is smaller than other objective. This solution is recognized with a probability of exp -(f(S’)-f(S))/T . The term of T is a control variable and it is reduced from a relatively to near zero. SA was applied by solving the fitness function formulated in Equation 4 and Equation 5. In (Table 1), the range of processing variables was defined to solve the surface roughness problem. Initial points of SA are given in (Table 7).

Matlab optimization toolbox was used to solve the surface roughness problem. The values of the SA parameters such as annealing function, reannealing interval, temperature function, initial point of temperature, and probability function. Matlab software was used to determine the best optimal results using different simulated angling parameters. The best solution of these parameters was given in (Table 8).

Table 8: Selection of SA parameters to determine the optimal solution. Minimum R a value was resulted with spindle speed of 17377 rpm, feed rate of 2,012 m/min, tool radius of 8,00 mm and depth of cut of 2,009 mm using simulated angling algorithm. Minimum R z value was obtained spindle speed of 16980 rpm, feed rate of 2,004 m/min, tool radius of 8,001 mm and depth of cut of 2,003 mm using simulated angling algorithm. The results of Matlab SA toolbox were given in Figure 9a, Figure 9b, Figure 10a, Figure 10b. According to Figure 9a, Figure 10a, the best fitness values for R a and R z were 2,474 µm and 10,674 µm, respectively. In the Figure 10a, Figure 10b, the optimal solutions of R a and R z were performed at the 3106 th and 4712 th iteration, respectively. Figure 9: (a) Variation of fitness function and best individuals (b) Converged values of parameters. Figure 10: (a) Variation of fitness function and best individuals (b) Converged values of parameters.

Comparison of results

The Taguchi orthogonal design based desirability function and desirability function integrated simulated angling algorithm results were given in (Table 9). From the (Table 9), the minimum R a and Rz values were 2,474 µm and 10,674 µm. For this reason, the RSM-DF-GA approach should present an efficient methodology in order to minimize the surface roughness.

Table 9: Results of experiments and optimum values by Taguchi-RSM-DF and Taguchi-RSM-DF-SA models for R a and R z Conclusions

In this work, effect of machining of Cedar of Lebanon pine (Cedrus libani), the surface roughness characteristics were investigated. The optimization process was adopted by a combined approach of L 27 orthogonal array based simulated angling algorithm. A three-levels and four factors experimental design were applied to identify the most significant factors, desirability function and simulated angling algorithm were used to optimize the machining procedure. The results can be summarized as follows:

Different from other studies in the literature, two mathematical models were developed based on L 27 orthogonal design and RSM method that the method has not been used to determine optimum CNC parameters for wood material.

The R a and R z were influenced by interaction between feed rate and spindle speed with PCR (%) of 37,16 % and 34,47 %, respectively.

Two different mathematical models were applied by using quadratic models. Moreover, backforward elimination method was applied to determine the non-significant variables. The non-significant terms were removed in these models.

Minimum R a value was obtained spindle speed of 17377 rpm, feed rate of 2,012 m/min, tool radius of 8 mm and depth of cut of 2,009 mm by using desirability function based simulated angling algorithm. For R z these results were found as 16980 rpm 2,004 m/min, 8,001 mm and 2,003 mm.

In this study, these models can be regarded as a method experimentally and statistically for the analysis, modeling optimization of CNC machining operations and can applied wood cutting process.

Acknowledgments

The authors thank to TUBITAK (The Scientific and Technological Research Council of TURKEY) for support and contributions during the doctoral thesis.

References:

Antony, J. 2014. Design of experiments for engineers and scientists. 2th edn. Elsevier: London. [ Links ]

Asiltürk, I.; Neseli, S.; Ince, A.M. 2016. Optimization of parameters affecting surface roughness of Co28Cr6Mo medical material during CNC lathe machining by using the Taguchi and RSM methods. Measurement 78:120-128. [ Links ]

Azhiri, R.B; Teimouri, R.; Baboly, M.G.; Leseman, Z. 2014. Application of Taguchi, ANFIS and grey relation analysis for studying, modeling and optimization of wire EDM process while using gaseous media. Int J Adv Manuf Technol 71(1-4):279-295. [ Links ]

Bao, X.; Ying, J.; Cheng, F.; Zhang, J.; Luo, B.; Li, L.; Liu, H. 2017. Research on neural network model of surface roughness in belt sanding process for Pinus koraiensis. Measurement 115:11-18. [ Links ]

Bharathi, S.R; Baskar, N. 2011. Particle swarm optimization technique for determining optimal machining parameters of different work piece materials in turning operation. The International Journal of Advanced Manufacturing Technology 54 (5-8):445-463. [ Links ]

Coelho, C.L.; Carvalho, L.M.H.; Martins, J.M.; Costa, V.A.V.; Masson, D.; Meausooner, P.J. 2008. Method for evaluating the influence of wood machining conditions on the objective characterization and subjective perception of a finished surface. Wood Sci Technol 42(3):181-195. [ Links ]

Davim, J.P.; Clemente, V.C.; Silva, S. 2009. Surface roughness aspects in milling MDF (medium density fiberboard). Int J Adv Manuf Technol 40(1-2):49-55. [ Links ]

Deepanraj, B.; Sivasubramanian, V.; Jayaraj, S. 2017. Multi-response optimization of process parameters in biogas production from food waste using Taguchi-Grey relational analysis. Energy Conversion and Management 141:429-438. [ Links ]

Gaitonde, V.N.; Karnik, S.R.; Davim, J.P. 2008. Taguchi multi-performance characteristics optimization in drilling of medium density fibreboard (MDF) to minimize delamination using utility concept. J. Mater Process Tech 196(1-3):73-78. [ Links ]

Hazir, E.; Erdinler, E.S.; Koç, K.H. 2018. Optimization of CNC cutting parameters using design of experiment (DOE) and desirability function. Journal of Forestry Research 29(5):1423-1434. [ Links ]

Hazir, E.; Ozcan, T. 2019. Response surface methodology integrated with desirability function and genetic algorithm approach for the optimization of CNC machining parameters. Arabian Journal for Science and Engineering 44(3):2795-2809. [ Links ]

Hiziroglu, S.; Kosonkorn, P. 2006. Evaluation of surface roughness of Thai medium density fiberboard (MDF). Building and Environment 41(4):527-533. [ Links ]

International Organization for Standardization. ISO. 1997. Geometrical product specifications (GPS) surface texture: profile method-terms, definitions, and surface texture profile method terms, definitions and surface texture parameters. ISO 4287. 1997. International Organization for Standardization: Geneva. [ Links ]

Kant, G.; Sangwan, K.S. 2014. Prediction and optimization of machining parameters for minimizing power consumption and surface roughness in machining. Journal of Cleaner Production 83:151-164. [ Links ]

Koç, K.H.; Erdinler, E.S.; Hazir, E.; Öztürk, E. 2017. Effect of CNC application parameters on wooden surface quality. Measurement 107:12-18. [ Links ]

Magoss, E. 2008. General regularities of wood surface roughness. Acta Silv Lign Hung 4:81-93. [ Links ]

Mahes, G.; Muthu, S.; Devadasan, S.R. 2015. Prediction of surface roughness of end milling operation using genetic algorithm. Int J Adv Manuf Technol 77(1-4):369-381. [ Links ]

Majumder, H.; Paul, T.R.; Dey, V.; Dutta, P.; Saha, A. 2017. Use of PCA-grey analysis and RSM to model cutting time and surface finish of Inconel 800 during wire electro discharge cutting. Measurement 107:19-30. [ Links ]

Ozdemir, T.; Hiziroglu, S. 2007. Evaluation of some sanding factors on the surface roughness of particleboard. Silva Fennica 41(2):373-378. [ Links ]

Philbin, P.; Gordon, S. 2006. Recent research on the machining of wood-based composite materials. International Journal of Machining and Machinability of Materials 1(2):186-201. [ Links ]

Prakash, S.; Palanikumar, K. 2010. Modeling for prediction of surface roughness in drilling MDF panels using response surface methodology. J Compos Mater 45(16):1639-1646. [ Links ]

Rao, R.V.; Kalyankar, V.D. 2013. Parameter optimization of modern machining process using teaching-learning-based optimization algorithm. Engineering Application of Artificial Intelligence 26(1): 524-531. [ Links ]

Rao, V.K.; Murthy, N.S.G.B.P. 2018. Modeling and optimization of tool vibration and surface roughness in boring of steel using RSM, ANN and SVM. J Intell Manuf 29(7):1533-1543. [ Links ]

Ratnasingman, J.; Scholz, F. 2006. Optimal surface roughness for high-quality on Rubberwood. Holzals Roh-und Werkstoff 64(4): 343-345. [ Links ]

Samanta, B. 2009. Surface roughness prediction in machining using soft computing. International Journal of Computer Integrated Manufacturing 22(3):257-266. [ Links ]

Sarikaya, M.; Güllü, A. 2016. Taguchi design and response surface methodology based analysis of machining parameters in CNC turning under MQL. Journal of Cleaner Production 65:604-616. [ Links ]

Selaimia, A.A.; Yallese, M.A.; Bensouilah, H.; Meddour, I.; Khattabi, R.; Mabrouki, T. 2017. Modeling and optimization in dry face milling of X2CrNi18-9 austenitic stainless steel using RSM and desirability approach. Measurement 107:53-67. [ Links ]

Sofuoglu, D.S. 2017. Determination of optimal machining parameters of massive wooden edge glued panels which is made of Scots pine (Pinus sylvestris L.) using Taguchi design method. Eur J Wood Prod 75(1):33-42. [ Links ]

Taguchi, G.; Chowdhury, S.; Wu, Y. 2005. Taguchi’s quality engineering handbook. Wiley: Hoboken. [ Links ]

Wilkowski, J.; Czarniak, P.; Grześkiewicz, M. 2011. Machinability evaluation of thermally modified wood using the Taguchi technique. In: COST Action FP0904 Workshop Mechano-chemical transformations of wood during Thermo-Hydro-Mechanical processing: 109-111. [ Links ]

Yang, W.H.; Tarng, Y.S. 1998. Design optimization of cutting parameters for turning operations based on the Taguchi Method. J Mater Process Technol 84(1-3):122-129. [ Links ]

Zhou, Y.; Gong, Y.; Zhu, Z.; Gao, Q.; Wen, X. 2016. Modeling and optimization of surface roughness from micro grinding of nickel-based single crystal super alloy using the response surface methodology and genetic algorithm. Int Adv Manuf Technol 85(9-12):2607-2622. [ Links ]

Zhou, Y.; Ren J.;Yao, C. 2017. Multi-objective optimization of multi-axis ball-end milling Inconel 718 via grey relational analysis coupled with RBF neural network and PSO algorithm. Measurement102:271-285. [ Links ]

Received: August 11, 2018; Accepted: June 12, 2019

Corresponding author: ender.hazir@istanbul.edu.tr This is an open-access article distributed under the terms of the Creative Commons Attribution License