1. Introduction

There are several methods in dimension stone block. Nowadays, diamond wire saw with bead diamond is the most widely used process for dimension stone quarries.

In diamond wire cutting operations, the cutting action primarily includes the pulling of continuous loops of spinning wire mounted with diamond beads through the dimension stone. In this cutting operation, firstly the horizontal cutting is done for avoiding the gravity effect of dimension stone block. Then, vertical cuttings are done. The initial step for a vertical cutting is to drill two holes that intersect at a 90° angle. Then, the diamond wire threaded through these holes, and over the drive wheel, clamps at the two ends to form a continuous loop. The diamond wire cutting machine is mounted on temporary track, along which it reverses to maintain tension in the wire as it cuts through the stone. Water is applied with the spin direction of the wire as a coolant and as a means of removing the particles of stone [1].

Diamond wire saw wear in rock cutting is one of the major criteria in determining the diamond wire saw life, energy consumption, production cost, and determine the cutting method selected for a given rock type. There are some important factors, which need to be considered to evaluate the wear rate of diamond wire saw. These factors can be divided into three key categories: (1) the characteristics of diamond wire saw, (2) the operating parameters and (3) the characteristics of the cut rock. Many researchers have attempted to investigate the effect of these parameters on wear up to now. Some researchers have studied the wear of circular diamond saw blade and diamond wire saw in rock cutting process [1–12]. In the field of diamond wire saw wear, Özçelik et al. [1] studied the effects of textural properties on marble cutting with diamond wire. They evaluated the relationships between textural characteristics and wear rate. The results showed that decreasing grain size increases the wear rate. In addition, there is a significant relationship between the texture coefficients and wear on diamond beads. This study indicated that textural characteristics can be considered in selection and design of diamond beads in marble industry. Özçelik and Kulaksız studied the relationship between cutting angles and wear on diamond beads in diamond wire cutting process [9]. Özçelik et al. investigated the wear rate of diamond beads in the cutting of different rock types. They used ridge regression method to evaluate the wear of beads in cutting of different rock with different physical, chemical, mechanical and mineralogical-petrographical properties. They concluded that the developed statistical models can be used to determine diamond wire life and cutting efficiency [10]. Similarly, Özçelik applied the multivariate statistical analysis of the wear on diamond beads in the cutting of andesitic rocks according to physical and mechanical properties of rock [11]. Mikaeil et al. predicted the wear of diamond wire saw with respect to the uniaxial compressive strength, Schimazek F-abrasivity factor, Shore hardness, and Young's modulus using the harmony search algorithm. The results showed that the applied algorithm can be used to evaluate the wear of diamond bead [12]. Almasi et al., carried out an investigation for the 11 types of igneous rocks based on the rock properties and production rate. For this purpose, they used linear and nonlinear regressions for analysis. The results indicated that the developed model can be suitable system to predict the wear rate of diamond beads [13]. Mikaeil et al., investigated different carbonate rocks in some famous quarries located in Iran, according to some important mechanical and physical properties of stone such as elasticity modulus, equivalent quartz content and uniaxial compressive strength. They used the application of multivariate regression analysis to evaluate the performance of diamond wire saw [14].

All of these studies were simply studied the diamond bead wear with statistical analysis and metaheuristic algorithm. No study has been found relating to the influence of rock characteristics on the diamond bead wear rate in diamond wire sawing with soft computing such as artificial and intelligence algorithm. In this research, it is aimed to develop an accurate, practical and applicable model for predicting the wear rate of diamond bead based on rock properties using intelligent systems. The remainder of this paper is organized as follows. In Section 2, methodology is briefly summarized. Section 3 presents the rock properties and laboratory testing of case study. In Section 4, the development of the MLP and a combination of GA-ANN models for wear rate prediction are explained. Section 5 discusses and assesses the results and performances of moldings. Finally, Section 6 gives conclusions and recommendations for future work.

2. Methodology

The methodology of this study is organized as following steps.

Step 1: Quarries studies (Cutting of dimension stone with diamond wire saw and determination of wear rate and sampling of stone blocks)

Step 2: Laboratory studies (Preparation of cylindrical specimens from stone samples and determination of physical and mechanical properties)

Step 3: Investigation of relationship between wear rate and characteristics of rock with GA-ANN and MLP

Step 4: Evaluation of results 

A flowchart followed in this study is illustrated in Figure 1.

Fig. 1. Flowchart of study

3. Quarries and laboratory studies

This study has been performed at 38 quarries including different type of rocks including andesites, limestones and real marbles in eleven areas in Turkey (Table 1). The extractions of andesite, limestones and real marble blocks have been achieved by diamond wire sawing. Wear rates of diamond beads for any rock types from 38 different locality have been recorded.

After quarries studies, experimental studies were done on rock block samples. To determine the main physical and mechanical properties of studied rocks, laboratory studies were done. Uniaxial compressive strength (UCS), Shore hardness (SH), Young modulus (YM) and Schmiazek F-abrasivity (SF-a) were selected as major rock characteristics. Results of laboratory tests for studied rock and quarries studies are given in Table (1).

Table 1. Diamond bead wear rate and characteristics of studied rock for evaluated deferent rocks

4. Prediction of wear rate through statistical and intelligent techniques

In the present study, two intelligent systems, namely GA- ANN and MLP are proposed to create a precise equation for the prediction of diamond wire saw performance and then a comparison of their performances are conducted and discussed. Some statistical functions indices, i.e., root mean square error (RMSE), coefficient of determination (R²) and standard deviation (STD) were computed to check for assessment and evaluating the accuracy and performance capacity of models as shown in Eqs.1 to 4, respectively.

Where n explains the number of data sets. The y_i and x_i are the forecasted and measured wear rate values, respectively. Note that, in a modeling with high and acceptable accuracy, the values of RMSE, R² and STD should be close to 0, 1 and 0, respectively.

4.1. Multilayer perceptron (MLP)

The soft computing acts as huge incentive to solve complex problems [15-17]. Artificial neural networks have a special place among soft computing methods considering their high ability in complex and imprecise data analysis and processing. The performance of human brain and neural systems considering million years of evolution can be used as the most complete and efficient pattern for the recognition of the surrounding events. In recent decades, neural networks have had a great impact on the development and modeling of industrial problems, as well as the control and optimization of the production process. One of the most practical and appropriate types of neural networks is the multilayer perceptron network used with a special type of learning algorithm in optimization problems. In a multilayer perceptron network, the linear relation between input and output vectors is shown through connections between neurons in each node and previous and next layers. The weight of network is determined through the minimum error between input and output data and or through the end of number of teachings to a predetermined value [18-19]. Different methods are used for teaching artificial neural networks, among which the back propagation algorithm is one of the most efficient and appropriate methods for teaching the multi-layer perceptron neural network and has the maximum consistency with this network. Therefore, for learning weights of a multi-layer perceptron network, the back propagation rule is used. This method was proposed by Wiliams Rumelhalt in 1986. In this method, using the gradient descent, it is attempted to minimize the square error between network outputs and objective function [20-21]. In fact, the error produced by the comparison between output data and estimated data must be smaller than the mean square error (MSE) or the root mean square error (RMSE); otherwise, the network must be propagated back in order to correct weights and reduce errors. The computation of output sensitivity to weights is started from the end of network and finally weights are updated at once. The network’s output error is computed in Eq (5) based on BP:

                                                        (5)

where  is the total output error, d is training samples, k is output data set,  and  are kth value of objective function (corresponding to the kth output unit) for dth training sample and kth value of output function (corresponding to the kth output unit) for dth training sample, respectively.

4.1.1. MLP modeling

In this work, Multilayer Perceptron (MLP) is utilized to build a prediction model for the assessment of diamond wire saw performance from the samples of stones as data set using the MATLAB software. The same datasets is used in three simulations in this study. Wear rate was considered as the dependent variable (output) and the uniaxial compressive strength (UCS), Schmiazek F-abrasivity (SF-a), Shore hardness (SH), and Young's modulus (YM) were considered as the independent variables (input). The dataset of 38 different varieties of dimension stones from Turkey quarry mines is considered in the current study, while 26 data points (70%) are utilized for constructing the model as train data, 8 data points (20%) are used as test data and the rest data points (4 data points) are considered as validation data for evaluation of the degree of accuracy and robustness.

The number of hidden layers and the number of neurons in each hidden layers are two important factors in MLP structure. Hence; in this modeling by contribution of experimental equations and after several simulations conducted, N_i =4 and N₀ =1 are considered as number of input neuron and number of output neuron, respectively. Furthermore hidden layer size is used as a range of 1-10 with one hidden layer for more accurate computing. Some of these equations are shown in Table (2).

Table 2. The equations for determining the number of neuron in hidden layer [22]

Different studies are conducted about conventional gradient descent techniques. Levenberg–Marquardt (LM) is one of the most effective and accurate algorithm based on suggestion of Hagan and Menhaj [31]. Therefore, Levenberg–Marquardt (LM) learning algorithm is considered in constructing MLP models for training net. In this study, the tansig and purelin are considered as transfer functions of the hidden layers and output, respectively. The effects of hidden layer size on the results of RMSE, R² and STD. are shown in Tables (3) and ranking of each model are listed in Table (4) based on a simple ranking method [32-33].

Table 3. Effects of hidden layer size on statistical functions performance in MLP network

Table 4. Ranking of each model using MLP network

Furthermore, Fig (2) shows correlation of determination between measured and predicted wear rate that there is a reasonable R² with a coefficient higher than 0.88. Figs(3) and (4) present RMSE values  and the histograms of errors for training, validation, and testing steps for all data set in the eighth model using MLP, respectively. According to the statistical functions and their total rank in Table (4), the model number 8 indicates higher performance capacities compared to other models with hidden layer size of 8 and a total rank of 70.

Fig. 2. R² of predicted and measured WR values for all data set using MLP model

Fig. 3. RMSE values for training, validation and testing steps

Fig. 4. Histograms of errors for all data set using MLP model

4.2. Hybrid GA-ANN algorithm

In recent years, Meta-heuristic algorithms have attracted the interest of researcher in many engineering fields and industry. These algorithms are the precise scientific tools instead of statistical methods to deal with the uncertain systems [34-36]. Genetic algorithm is a population-based algorithm like particle swarm optimization, firstly proposed by John Holland in 1975 at the University of Michigan [37]. The GA is a Meta-heuristic algorithm that is suitable for dealing with complex problems, especially when the goal is to find an optimization result. The genetic algorithm can model the qualitative and quantitative aspects of uncertain systems in industry. The genetic algorithm (GA) was inspired by Darwin's principle of natural evolution. The optimization and searching method in this algorithm is based on the principle of natural biological evolution and inheritance rules.

In the genetic algorithm, numbers are expressed in terms of binary strings and converged toward the range of solutions during the stepwise algorithm’s implementation using the probability distribution function. The population to population searching is a technique for obtaining an optimal solution. In addition, in problems with a complex hypothesis space with the unknown effect of components on the general hypothesis, GA can be used for searching and finding an approximate solution for an optimal answer. GA has a significant flexible nature compared to other Meta-heuristic algorithms and is formed based on the natural selection mechanism and stochastic techniques. Furthermore, in this algorithm, differentiation is not required and only the objective function and basic information fitting method is used. In GA, each set of chromosomes and each replication of algorithm are called population and generation, respectively. GA first fits the existing population in each replication by determining the initial population and using the fitness function and then starts optimization. In fact, the compatibility of initial population is computed and assessed through the objective function. Next, the new generation (population) is produced based on GA operators, i.e. reproduction, crossover and mutation. The fitness steps for answers and production of new generations continue until an optimal answer is reached. GA has a wide range in the optimization and solution of complex and uncertain problems. Landslide was evaluated and studied by Terranova et al. The results and the subsequent validation showed that the genetic-algorithms-based hydrological model was a reliable approach for the purpose of their research. [38]. The flood risk management was done by Woodward et al using a multi objective genetic algorithm. The results indicated that the simulations were very suitable [39].

In addition, one of the most important applications of GA is teaching neural networks. Since GA is able to run away from trapping in local optimums, does not depend on any special structure of the network and is applied for any defined structure, it can be considered as a proper and efficient tool for being combined with neural networks and teaching neural networks. Therefore, in this research, the wear rate is anticipated using a GA and ANN combination for optimizing the performance of diamond wire saw. Figure (5) shows a combination of GA-ANN structure.

Fig. 5. The structure of hybrid GA-ANN flowchart [40]

In recent years, different studies were conducted in different scientific areas using a hybrid GA-ANN algorithm. A model was proposed by Armaghani and Khandelwal for the anticipation of the drilling rate index using GA and neural network and used on the rock strength characteristics [41]. A good anticipation of the pile bearing capacity was done by Momeni et al by developing a hybrid algorithm model. In their study, they obtained answers with a very good accuracy by combining GA and ANN [42]. A model for the anticipation of flyrock and backbreak in open pit mines was proposed by Monjezi et al based on GA-ANN and a model with high efficiency was developed for anticipation with the minimum possible error and the maximum correlation coefficient [43].

4.2.1. Hybrid GA-ANN algorithm modeling

As mentioned, the predicting diamond wire saw performance in this study is based on four important measures of rocks which the same datasets performed in the assessment of MLP simulations were applied. A hybrid GA-ANN algorithm is considered as a flexible predicting method. This technique is based on artificial intelligence for solving complex issues and uncertain systems, which is one of the most efficient soft computing methods. Hagan and Menhaj (1994) introduced more details for the hybrid GA-based ANN model [31].

In order to obtain a high level of precision in data analysis and predicting process, it is necessary to determine the appropriate control parameters. Hence; firstly the pseudo-code of the hybrid GA-ANN algorithm is written in MATLAB software. Some parameters can define based on visual observations and suggestion of previous studies [31][44]. The recombination percent (RP) was determined at 15%, the mutation percent (MP) and cross-over percent (CP) were fixed at 35% and 50%. The maximum number of generation (G_Max) and population size are two of the effective factors in during algorithm implementation process. In the next step, in order to obtain the optimum G_Max value, the efficacy of the number of generation on the network performance for RMSE as statistical functions is carried out for deferent population size as a range of 50-500 with G_Max = 500.The result of analysis is illustrated in Fig (6). Based on the results, it’s obvious that the optimum G_Max was set to be 400, because the network performance is unchanged after this value of G_Max. Furthermore, the ANN structure is determined based on Table (2).

Fig. 6. The efficacy of the number of generation on the network performance based on RMSE

In the final step, in order to determine the optimum population size in GA-ANN algorithm, 10 hybrids GA-ANN models are constructed for the optimum G_Max = 400, the results of root mean square error (RMSE), coefficient of determination (R²) and standard deviation (STD) are listed in Table (5). The models are ranked according to suggestion of Zorlu et al. (2008) as a simple ranking approach [33].

Table 5. Effects of population size on statistical functions performance in Hybrid GN-ANN algorithm

The results of ranking shown in Table (6). As shown in Table (6), the third model has highest rank among other models with rank of 44. Therefore, Figs (7) and (8) illustrate the value of R² and the graphical comparison between measured and predicted wear rate using hybrid model for test data set of the third simulation, respectively.

Table 6. Ranking of each model using Hybrid GN-ANN algorithm

Fig. 7. R² between measured and predicted wear rate using hybrid GA-ANN model for Test data

Fig. 8. The graphical comparison between measured and predicted wear rate using hybrid GA-ANN model for test data

5. Evaluation of the results and discussion

Assessment of diamond wire saw performance has been considered as one of the most notable topics to study in mining engineering and rock mechanics. Investigation of effective parameters on diamond wire saw performance is frequently encountered with complex and non-linear problems; hence, soft computing techniques are the approaches that are suitable for dealing with complex and uncertain processes. In this research, in order to investigate the applications and performances of each optimization technique for predicting diamond wire saw performance, after implementation simulations, the most appropriate structures of MLP and hybrid GA-ANN model were determined. Based on the optimum models, all the datasets were randomly selected to 4 various datasets. Also, 80% of samples were used randomly for training dataset and the remaining 20 % was considered as testing data in each simulation. In Tables (7), the results of R², RMSE and STD. of ANN and GA-ANN methods and final ranking are shown. Finally, Table (8) listed shows the final ranking of all simulations.

Table 7. Ranking of each model for 4 datasets randomly selected using MLP and hybrid GN-ANN

Table 8. Final ranking of each model using two optimization techniques

Based on the results of Tables (7) and (8), the first model in MLP simulation and the second model in GA-ANN simulation obtained the highest score according to statistical functions. The first MLP simulation has highest total rank with 21 among other MLP simulations. Also, the second GA-AAN simulation obtained the most score in other GA-ANN simulations with rank=19. In comparison between the best MLP and hybrid GA-ANN simulations, the high and reasonable R² values between the model predictions and the measured data for training = 0.95 and testing =0.991 using the hybrid GA-ANN method describes its high capability in the prediction of diamond wire saw performance.

6. Conclusion

In this research, the aim is to develop prediction models for assessment of diamond wire saw performance using two optimization approaches including MLP and hybrid GA-ANN algorithm. This work compared the application and performance of hybrid GA-ANN algorithm with MnLR on the database of 38 different varieties of dimension stones from Turkey quarry mines. Four important physical and mechanical rock characteristics on cutting process, namely wear rate from diamond wire cutting machine was set as output data and uniaxial compressive strength, Schmiazek F-abrasivity, Shore hardness, and Young's modulus were considered as input data set. From the results found in this study, it can be concluded that the performance of hybrid GA-ANN algorithm is superior to MLP in terms of some model performance indices such as RMSE, R² and STD. The comparison was made between the three simulations based upon the performance indices, hybrid GA-ANN algorithm with coefficient of determination (R²) of training = 0.95 and testing = 0.991 was selected as the best predictive model. In addition, it definitely outperforms MLP based on robustness and solution quality for simulation some problems involved in rock mechanics engineering. From what has been discussed above, it can be concluded that hybrid GA-ANN algorithm is a reliable system simulation technique for predicting the performance of diamond bead with highly acceptable level of accuracy and it can be applied as an appropriate alternative which has a wide application in management and planning for costs and designs of quarries. In future research, prediction of diamond wire saw performance can also be investigated and improved the using ICA-ANN, PSO-ANN, Hybrid Harmony Search (HS-BFGS) and Grey Wolf Optimizer (GWO) and other performance indicators such as the mean absolute percentage error (MAPE) and value account for (VAF).

Acknowledgment

We would like to express our deepest thanks to Professor Mahdi Ghaem for his excellent advice.