A MODELING STUDY TO EVALUATE THE QUALITY OF WOOD SURFACE

The goal of this study was to develop a model to predict sanding conditions of different type of materials such as Lebnon cedar (Cedrus libani) and European Black pine (Pinus nigra). Specimens were prepared using different values of grit size, cutting speed, feed rate, and sanding direction. Surface quality values of specimens were measured employing a laserbased robotic measurement system and stylus type measurement equipment. Full factorial design based Analysis of Variance was applied to determine the effective factors. These factors were used to develop the Artificial Neural Networks models for two different measurement systems. The MATLAB Neural Network Toolbox was used to predict the Artificial Neural Networks models. According to the results, the Artificial Neural Networks models were performed using Mean Absolute Percentage Error and R-square values. Mean Absolute Percentage Error values for laser and stylus equipment were found as 2,405 % and 3,766 %, respectively. R-square values were determined as 96,2% and 92,7 % for laser and stylus measurement equipment, respectively. These results showed that the proposed models can be successfully used to predict the surface roughness values.


INTRODUCTION
Sanding process is applied in different wood manufacturing applications.In furniture industry, this process presents some advantages such as low surface roughness, high productivity, appearance of wood products and high wood coating performance (Richter et. al 1995, Cool and Hernandez 2011, Scrinzi et al. 2011, Landry and Blanchet 2012, Hiziroglu et al. 2014, Gurau et al. 2015, Ugulino and Hernandez 2016, Sogutlu et al. 2016).Therefore, the surface roughness is the major indicator of wood surface quality.It is mainly a result of various controllable or uncontrollable sanding parameters.Anatomical structure, hardness, density, annual ring variation, cell structure early-late wood ratio are uncontrollable variables, while the sanding parameters such as feed rate, cutting speed, grit size, pressure, depth of cut, sandpaper type and cutting direction are controllable variables (Tan et al. 2012, Gurau et al. 2013, Magoss 2015, Ramananantoandro et al. 2017).
There are different roughness measuring techniques such as pneumatic, laser and light scattering techniques to determine the surface quality of wood and wood products (Hiziroglu and Suziki 2007, Hazir 2013, Zhong et al. 2013).Stylus type profilometer is commonly applied due to its usefulness and advantage in obtaining accurate numerical results (Sandak andTanaka 2003, Hiziroglu et al. 2014).Laser measurement system provides crucial advantages such as measuring the complex surface structures with non-contact equipment, decreasing the time consumption, gathering more data from the surfaces in a short time and obtaining online measurements in the real production process (Sandak andTanaka 2003, Koc et al. 2017).The full factorial design is a powerful technique to determine the significant factors.It involves all possible combinations between input and output variables.Therefore, this design has been widely used in various engineering applications.
In the recent years, artificial intelligence algorithms such as Artificial Neural Network (ANN), Genetic Algorithm (GA), Fuzzy Logic (FL) and Particle Swarm Optimization (PSO) have applied to different engineering problems (Ozsahin and Aydin 2014, Mahes et al. 2015, Jain and Raj 2017).Carrano et al. (2002) investigated sanding process of hard maple, white oak and eastern white pine as function of spindle speed, feed rate, depth of cut, grit size, tooling resilience and grain orientation.The results showed the grit size, tooling resilience and grain orientation were significant for all species.The feed rate was found as a significant factor for white oak and eastern white pine.Zhong et al.
(2013) evaluated surface quality using stylus type profilometer and 3D image analyzer different wood materials such as particle board, medium density fiberboard, plywood and ten different solid wood.According to the results, these methods can be successfully applied to determine the surface quality.Tiryaki et al. (2014) modeled planing and sanding process of Spruce and Beech wood as function of spindle speed, cutting depth, feed rate, number of cutter, wood zone and grain size of abrasive.The results indicated the surface roughness decreased with increasing the grit number and number of the cutter.It was shown that ANN method can also be used successfully for modeling of surface roughness.A study was carried out by Laina et al. (2017) investigated sanding process of beech, oak and pine the function of grain direction, wood hardness and machining conditions such as planing and sanding process.The results showed the surface roughness was decreased from 60 to 180 grit size.Hardness was found as a significant factor for wood surface roughness.Hazir et al. (2017) developed a mathematical model to evaluate optimum sanding conditions of European black pine (Pinus nigra).Samples were sanded using different grit sizes, feed rates, cutting speeds and depths of cut.Response surface methodology (RSM) was used to determine the optimum parameters values.
The objective of this study was to develop a model to predict sanding conditions using ANN of two different wood species such as Lebnon cedar (Cedrus libani A.Rich) and European Black pine (Pinus nigra Arnold) for two different measurement equipment, namely laser and stylus.

MATERIALS AND METHODS
Lebnon cedar (Cedrus libani A.Rich) and European Black pine (Pinus nigra Arnold) species are extensively used in the furniture industry.The samples were prepared with the dimension of 200 mm x 100 mm x 30 mm for each test.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 10±1%.The density of Lebnon cedar and Black pine was found as 570 kg/m³ and 680 kg/m³, respectively.The samples were processed with wide-belt sanding machine equipped with open coat aluminum oxide abrasive paper.

Evaluation of wood surface quality
In this study, wood surface quality was determined by using two different methods.(1) Laser-based system was employed to evaluate the surface quality.Cartesian robot integrated with laser sensor was used to gather data for evaluating the surface quality (Figure 1a).In order to determine the surface quality with the laser sensor, the robot was gathered with 500 data every 7 mm movement and the average values were calculated.Laser sensor was used to measure the wood surface for determining the various sanding parameters.The robot was controlled in X and Y axis for evaluating the wood surface quality.The laser sensor gave 500 measurements and the results were transferred into the MATLAB program.
(2) Another measurement equipment for determining the surface roughness values of machined wood material was the Taylor-Hobsan Suftronic type equipment (Figure 1b).This device is a stylus-A modeling study to evaluate..: Hazir and Koc based portable profilometer, equipped with a diamond stylus with a 5 µm radius 12,5 mm length of measurement 2,5 mm length of sampling, 15 mm travel of stylus and 90° contact angle running at a speed of 0,5 mm/s.were taken from sample surfaces.With reference to ISO 4287:1997, average roughness (R a ) and mean peak to valley height (R z ) are admitted as roughness parameters.In this study, R a parameter was selected to evaluate the surface roughness of the samples.

Statistical design of experiment
The statistical analysis was performed using Minitab software package 17.The full factorial experimental design was used to obtain the results.This design is one of the most important methods to investigate two or more parameters (Montgomery 1997).ANOVA was applied to experimental data in order to determine effective factors for both laser and stylus type equipment.Each independent variable had two and three levels which were coded as (-1), ( 0) and (+1).The low (-1), medium (0) and high (+1) levels are given in (Table 1).

Artificial neural networks (ANN)
Artificial neural networks (ANN) are developed with inspiration from information processing model of the biological neural system of the human brain.It is used non-linear and linear model for prediction and optimization of the data.ANNs are applied for engineering applications such as pattern recognition, forecasting and data processing (Karazi et al. 2009).This model consists of inputs, which are multiplied by weights.These weights are computed using mathematical function determining the activation of the neuron.This model learns the correlation between the input and output factors by using recorded data.An ANN system depends on neurons connected with the number of weighted links.Every piece of information is transferred to other neurons.Artificial neural network structure is given in (Figure 2).2: Where the terms of net j , w ij , j, x i , θ j and y i are sums of information, weight factors, neuron, layer information, bias of the layer and output values, respectively.In this study, a feedforward and backpropagation multilayer ANN was carried out predicting the wood surface roughness for two different measurement systems.Moreover, the hyperbolic sigmoid function (tansig) and the linear transfer function were selected as transfer function.Levenberg-Marquardt algorithm (trainlm) was applied as training algorithm and gradient descent with a momentum back-propagation algorithm (traingdm) was selected as learning rule.To evaluate the ANN model, Mean Squared Error (MSE), Mean Absolute Percentage Error (MAPE) and R-square (R 2 ) values were used to test the accuracy of the results.MAPE, MSE and R 2 were computed with Equation 3, Equation 4 and Equation 5.The terms of A t , F t and t F indicate the actual, predicted and the average of predicted values, respectively.
(3) (4) (5) Maderas.Ciencia y tecnología 20(4): 691 -702, 2018 A modeling study to evaluate..: Hazir and Koc In order to determine an equal contribution of each variable, these parameters were normalized using Equation 6.This process ensured the best generalization performance of ANN model.Training, testing and validation data were normalized by carrying out their minimum and maximum values within the range of [-1, 1].The model parameters were normalized by computing Equation 6. (6) The terms of X norm , X max, and X min are normalized value of a variable X, maximum and minimum values of X, respectively.

RESULTS AND ANALYSIS
The experiment consisted of five factors and one replicate, making a total of 108 runs ( ). (Table 2) shows the experimental parameters and their recoded laser and stylus measurement roughness values.In both measurement methods, the average roughness values were calculated by taking measurements from three different points on the wooden surface.
ANOVA based F-test was applied to evaluate the significance factor on the surface roughness.This analysis investigates the following for each parameter given in Equation 7and Equation 8: ...

:
for at least one pair , The F value is calculated by: (8) The terms of (α-1) and (N-α) are the degrees of freedom and the error degrees of freedom for the parameter A, respectively.MS A and MS E are indicated the sum squares of means and errors for the variable A, respectively.The null hypothesis is rejected when the F 0 is higher than the critical value of F a,a-1,N-a , where α is the level of the significance (Antony 2014). (Table 3 and Table 4) displayed the P -value is less than 0,05 showing the model is significant at 95% confidence level.Value of "prob>F" are lower than 0,05 indicating that the model terms are significant.In this case, the terms of feed rate, cutting speed, material type, sanding direction and grit size resulted in significant factors for stylus and laser measurement equipment.DF: degrees of freedom, SS: Sum of squares, F: F-test value and P:error variance ª At a given response, parameters belonging to the filled cells are effective within 95 % reliability interval.

Evaluation of the models
The normal probability plot of the residuals and residuals versus the predicted for R a are shown in (Figure 3).Evaluating on the normal probability plots (NPP) are shown in Figure 3a, Figure 3b) depicts that the residuals generally fall on a straight line implying that the errors are distributed normally.In addition to these Figure 3c, Figure 3d) show that the residuals versus the fitted values for the surface roughness data.As result of the residuals no unusual structure is apparent.This implies that the models proposed are adequate and there is no reason to suspect any violation of the independence or constant variance assumption (Montgomery 1997).

Parameter prediction by using ANN
According to the results obtained from the analysis of variance, the sanding parameters of grit size, feed rate, cutting speed, material type and sanding directions were found as effective factors on surface quality for laser and stylus type equipment.In this reason, these variables were selected as input parameters while the surface roughness was selected as output parameter for ANN.The gathered data making of 108 was used as 76 samples for the training, 16 samples for the validation and 16 samples for the testing.The accuracy of models was performed by using correlation coefficient (R 2 ) and MSE values.These results were given in (Table 5).; 1,38 and 1,97 for stylus-based measurement equipment, respectively.These results were used to evaluate the performance for predicting models.Because of the high values of R 2 and the low level of errors, these models were satisfactory.A value of R 2 obtained from laser-based measurement was better than stylus type equipment.(Figure 5 and Figure 6) showed the relationship between the measured and predicted values for laser and stylus type equipment.Moreover, the MAPE values for all data are given in (Figure 5 and Figure 6).MAPE values for laser and stylus measurement equipment were computed as 2,405 % and 3,766 %, respectively.MAPE values for laser measurement equipment lower than stylus type equipment.
(a) Laser-based measurement system (b) stylus-based measurement equipment.

Figure 2 .
Figure 2. Artificial neural network structure.An ANN is performed by using transfer function type, training algorithm, training and testing data size and values of weights and biases.This algorithm is formulated in Equation 1 and Equation 2:

Figure 3 .
Figure 3. (a-b) NPP of residuals and (c-d) plot of residuals fitted for R a -Stylus and Laser results.

Figure 4
Figure4displays the relationship between the measured and predicted values for training, validation and testing data.(Figure4a) displays R 2 values for training, validation and testing data sets in predicting stylus measure of 0,9450; 0,92462 and 0,87836 respectively.(Figure4b) displays R 2 values for training, validation and testing data sets in predicting laser measuring of 0,97528; 0,94573 and 0,90692 respectively.

Figure 4 .
Training, validation and testing results for (a) stylus and (b) laser type equipment.According to Figure 4, R 2 values are close to 1 for training, validation and testing data.It means that there is a good relationship between the measured values and the predicted values.Moreover, MSE values of training, validation and testing parameters were found as 1,45; 7,44 and 6,25 for laser-based measurement system, respectively, whereas MSE values of training, validation and testing parameters were computed as 8,59

Figure 5 .
Figure 5.The relationship between the measured and predicted values for laser measurement.

Figure 6 .
Figure 6.The relationship between the measured and predicted values for stylus measurement.

Table 1 .
Sanding procedure parameters and levels.

Table 2 .
Experimental parameters and the recorded roughness values of laser and stylus measurement.

Table 3 .
ANOVA for R a -Stylus measurement results.At a given response, parameters belonging to the filled cells are effective within 95 % reliability interval.
DF: degrees of freedom, SS: Sum of squares, F: F-test value and P:error variance ª

Table 4 .
ANOVA for -R a Laser measurement results.

Table 5 .
ANN performance results.