Document Type: Regular Article
Authors
^{1} Department of Mining and Metallurgical Engineering, Urmia University of Technology, Urmia, Iran
^{2} Young Researchers and Elite Club, Rasht Branch, Islamic Azad University, Rasht, Iran
^{3} Department of Mining Engineering, Hacettepe University, Ankara, Turkey
^{4} Graduate of Civil Engineering, Department of Civil Engineering, Islamic Azad University, Astara Branch, Astara, Iran
Abstract
Highlights
Keywords
Main Subjects
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 mineralogicalpetrographical 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 Fabrasivity 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 inﬂuence 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 GAANN 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.
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 GAANN and MLP
Step 4: Evaluation of results
A flowchart followed in this study is illustrated in Figure 1.
Start 
Quarries studies 
Laboratory studies 
Cutting of dimension stone with diamond wire saw 
Determination of wear rate 
Sampling of stone blocks 
Preparation of cylindrical specimens from stone samples 
Determination of physical and mechanical properties 
Investigation of relationship between wear rate and characteristics of rock 
Genetic AlgorithmArtificial Neural Network 
Multilayer Perceptron 
Evaluation of results 
Fig. 1. Flowchart of study
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 Fabrasivity (SFa) 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
Cutting No. 
Sample 
Cutting performance 
Rock characteristics 

WR (mm/m_{2}) 
UCS (MPa) 
SH

YM (GPa) 
SFa (N/mm) 

1 
Andesite 
Andesite1 
0.00150 
28.05 
33.00 
6.80 
0.104 
2 
Andesite2 
0.00230 
47.39 
39.93 
7.60 
0.134 

3 
Andesite3 
0.00360 
77.25 
65.00 
20.80 
0.331 

4 
Andesite4 
0.00340 
84.02 
63.62 
23.50 
0.290 

5 
Andesite5 
0.00210 
26.55 
49.93 
7.90 
0.043 

6 
Andesite6 
0.00560 
67.41 
67.00 
15.50 
0.531 

7 
Andesite7 
0.00160 
57.75 
43.70 
8.00 
0.171 

8 
Andesite8 
0.01570 
87.53 
61.38 
30.56 
0.313 

9 
Andesite9 
0.01680 
75.75 
63.70 
24.60 
0.254 

10 
Andesite10 
0.01590 
81.35 
62.48 
25.40 
0.292 

11 
Andesite11 
0.00850 
78.75 
60.30 
26.40 
0.272 

12 
Andesite12 
0.00800 
82.50 
61.20 
28.30 
0.320 

13 
Andesite13 
0.00360 
27.23 
42.50 
7.30 
0.133 

14 
Andesite14 
0.00220 
51.92 
41.20 
7.76 
0.203 

15 
Andesite15 
0.00360 
56.25 
43.50 
8.01 
0.221 

16 
Limestone 
Usak White 
0.00370 
69.00 
47.00 
12.20 
0.021 
17 
Kozagac White 
0.00280 
42.00 
40.00 
12.10 
0.034 

18 
Milas Lilac 
0.00260 
55.00 
46.00 
11.00 
0.154 

19 
Afyon Cream 
0.00300 
64.00 
46.00 
11.80 
0.014 

20 
Kutahya Lilac 1 
0.00820 
52.26 
42.90 
17.02 
0.004 

21 
Kutahya Lilac 2 
0.00980 
79.00 
43.05 
17.50 
0.004 

22 
Kutahya Violet 
0.00820 
63.49 
43.25 
21.14 
0.004 

23 
Afyon Violet 1 
0.00650 
74.19 
45.23 
21.43 
0.003 

24 
Afyon Violet 2 
0.00490 
51.84 
41.60 
15.96 
0.003 

25 
Afyon Gray 1 
0.00440 
49.02 
41.55 
13.07 
0.004 

26 
Afyon Gray 2 
0.00440 
45.57 
39.85 
15.72 
0.003 

27 
Mugla Nacre 
0.00290 
28.68 
50.45 
12.74 
0.003 

28 
Mugla White 
0.00150 
30.00 
32.90 
9.86 
0.002 

29 
Real marble 
Yesilova Beige 
0.00310 
70.50 
56.00 
9.90 
0.014 
30 
Sivrihisar beige1 
0.00370 
72.00 
60.00 
12.50 
0.014 

31 
Sivrihisar beige2 
0.00380 
70.00 
62.00 
13.20 
0.013 

32 
Sivrihisar beige3 
0.00390 
68.00 
58.00 
12.20 
0.012 

33 
Antalya Beige 1 
0.00610 
55.30 
58.70 
22.83 
0.061 

34 
Antalya Beige 2 
0.00850 
65.80 
58.15 
25.24 
0.104 

35 
Antalya Beige 3 
0.00710 
59.77 
58.00 
17.43 
0.122 

36 
Antalya Beige 4 
0.00820 
56.75 
59.48 
20.35 
0.370 

37 
Toros Black 1 
0.01990 
105.48 
65.30 
20.25 
0.122 

38 
Toros Black 2 
0.01573 
110.77 
64.75 
17.94 
0.121 
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^{2}) 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^{2} and STD should be close to 0, 1 and 0, respectively.
The soft computing acts as huge incentive to solve complex problems [1517]. 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 [1819]. 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 multilayer perceptron neural network and has the maximum consistency with this network. Therefore, for learning weights of a multilayer 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 [2021]. 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.
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 Fabrasivity (SFa), 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_{0} =1 are considered as number of input neuron and number of output neuron, respectively. Furthermore hidden layer size is used as a range of 110 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]
Researchers 
Heuristic 
HechtNielsen [23] 

Hush [24] 

Kaastra and Boyd [25] Kannellopoulas and Wilkinson [26] 

Ripley [27] 

Paola [28] 

Wang [29] 

Masters [30] 

N_{i}: Number of input neuron, N_{0}: Number of output neuron 
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^{2} and STD. are shown in Tables (3) and ranking of each model are listed in Table (4) based on a simple ranking method [3233].
Table 3. Effects of hidden layer size on statistical functions performance in MLP network
Model No. 
Hidden Layer Size (HLS) 
The Results of Network for R^{2} 
The Results of Network for RMSE 
The Results of Network for STD. 

Training 
Validation 
Testing 
Training 
Validation 
Testing 
Training 
Validation 
Testing 

1 
1 
0.54 
0.99 
0.7 
0.0026 
0.0021 
0.0048 
0.0027 
0.0024 
0.0042 
2 
2 
0.73 
0.81 
0.56 
0.0027 
0.0028 
0.0018 
0.0027 
0.0031 
0.0019 
3 
3 
0.7 
0.21 
0.14 
0.0028 
0.0042 
0.0023 
0.0029 
0.0044 
0.0018 
4 
4 
0.8 
0.84 
0.8 
0.0021 
0.0036 
0.003 
0.0021 
0.0039 
0.0028 
5 
5 
0.83 
0.78 
0.73 
0.002 
0.0037 
0.0023 
0.002 
0.0031 
0.0023 
6 
6 
0.96 
0.97 
0.54 
0.0011 
0.0028 
0.0049 
0.0011 
0.0022 
0.0042 
7 
7 
0.76 
0.84 
0.77 
0.0024 
0.0022 
0.0027 
0.0024 
0.0025 
0.0025 
8 
8 
0.95 
0.55 
0.92 
0.0011 
0.0035 
0.0021 
0.0011 
0.0038 
0.0017 
9 
9 
0.7 
0.92 
0.88 
0.0028 
0.0011 
0.0022 
0.0028 
0.0011 
0.0021 
10 
10 
0.67 
0.99 
0.91 
0.0027 
0.0018 
0.002 
0.0026 
0.002 
0.0021 
Table 4. Ranking of each model using MLP network
Model No. 
Hidden Layer Size (HLS) 
The Ranking of Network for R^{2} 
The Ranking of Network for RMSE 
The Ranking of Network for STD. 
Total rank 

Training 
Validation 
Testing 
Training 
Validation 
Testing 
Training 
Validation 
Testing 

1 
1 
2 
10 
4 
6 
8 
3 
5 
7 
3 
48 
2 
2 
5 
6 
3 
5 
6 
10 
5 
5 
8 
53 
3 
3 
4 
3 
1 
4 
2 
6 
3 
2 
9 
34 
4 
4 
7 
7 
6 
8 
4 
4 
8 
3 
5 
52 
5 
5 
8 
5 
5 
9 
3 
6 
9 
5 
7 
57 
6 
6 
10 
9 
2 
10 
6 
2 
10 
8 
3 
60 
7 
7 
6 
7 
7 
7 
7 
5 
7 
6 
6 
58 
8 
8 
9 
4 
10 
10 
5 
8 
10 
4 
10 
70 
9 
9 
4 
8 
8 
4 
10 
7 
4 
10 
4 
59 
10 
10 
3 
10 
9 
5 
9 
9 
6 
9 
4 
64 
Furthermore, Fig (2) shows correlation of determination between measured and predicted wear rate that there is a reasonable R^{2} 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^{2} 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
In recent years, Metaheuristic 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 [3436]. Genetic algorithm is a populationbased algorithm like particle swarm optimization, firstly proposed by John Holland in 1975 at the University of Michigan [37]. The GA is a Metaheuristic 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 Metaheuristic 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 geneticalgorithmsbased 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 GAANN structure.
Fig. 5. The structure of hybrid GAANN flowchart [40]
In recent years, different studies were conducted in different scientific areas using a hybrid GAANN 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 GAANN and a model with high efficiency was developed for anticipation with the minimum possible error and the maximum correlation coefficient [43].
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 GAANN 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 GAbased 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 pseudocode of the hybrid GAANN 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 crossover 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 50500 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 GAANN algorithm, 10 hybrids GAANN models are constructed for the optimum G_{Max} = 400, the results of root mean square error (RMSE), coefficient of determination (R^{2}) 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 GNANN algorithm
Model No. 
Population Size 
The Results of Network for R^{2} 
The Results of Network for RMSE 
The Results of Network for STD. 

Training 
Testing 
Training 
Testing 
Training 
Testing 

1 
50 
0.89 
0.25 
0.168 
0.489 
0.459 
0.435 
2 
75 
0.69 
0.64 
0.261 
0.375 
0.376 
0.384 
3 
100 
0.67 
0.81 
0.2 
0.55 
0.264 
0.337 
4 
150 
0.64 
0.62 
0.299 
0.357 
0.389 
0.375 
5 
200 
0.68 
0.68 
0.31 
0.264 
0.418 
0.442 
6 
250 
0.89 
0.24 
0.174 
0.495 
0.443 
0.389 
7 
300 
0.75 
0.56 
0.264 
0.41 
0.467 
0.595 
8 
350 
0.7 
0.37 
0.297 
0.417 
0.444 
0.447 
9 
400 
0.65 
0.7 
0.276 
0.348 
0.374 
0.407 
10 
500 
0.49 
0.69 
0.319 
0.318 
0.455 
0.467 
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^{2} 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 GNANN algorithm
Model No. 
Population Size 
The Ranking of Network for R^{2} 
The Ranking of Network for RMSE 
The Ranking of Network for STD. 
Total rank 

Training 
Testing 
Training 
Testing 
Training 
Testing 

1 
50 
10 
2 
10 
3 
2 
5 
32 
2 
75 
7 
6 
7 
6 
8 
8 
42 
3 
100 
5 
10 
8 
1 
10 
10 
44 
4 
150 
3 
5 
3 
7 
7 
9 
34 
5 
200 
6 
7 
2 
10 
6 
4 
35 
6 
250 
10 
1 
9 
2 
5 
7 
34 
7 
300 
9 
4 
6 
5 
1 
1 
26 
8 
350 
8 
3 
4 
4 
4 
3 
26 
9 
400 
4 
9 
5 
8 
9 
6 
41 
10 
500 
2 
8 
1 
9 
3 
2 
25 
Fig. 7. R^{2} between measured and predicted wear rate using hybrid GAANN model for Test data
Fig. 8. The graphical comparison between measured and predicted wear rate using hybrid GAANN model for test data
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 nonlinear 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 GAANN 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^{2}, RMSE and STD. of ANN and GAANN 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 GNANN
Optimization Techniques 
Model No. 
Value of R^{2} 
Score of R^{2} 
Value of RMSE 
Score of RMSE 
Value of STD. 
Score of STD. 
Total Score 
MLP 
Train 1. 
0.8 
3 
0.0014 
3 
0.0016 
3 
9 
Train 2. 
0.88 
4 
0.001 
4 
0.001 
4 
12 

Train 3. 
0.62 
1 
0.0034 
2 
0.0034 
2 
5 

Train 4. 
0.65 
2 
0.0049 
1 
0.0048 
1 
4 

Test 1. 
0.91 
4 
0.0013 
4 
0.0015 
4 
12 

Test 2. 
0.77 
2 
0.0044 
2 
0.0017 
3 
7 

Test 3. 
0.81 
3 
0.006 
1 
0.007 
1 
5 

Test 4. 
0.5 
1 
0.004 
3 
0.0033 
2 
6 

GAANN 
Train 1. 
0.968 
3 
0.11 
4 
0.88 
1 
8 
Train 2. 
0.95 
4 
0.15 
3 
0.54 
4 
11 

Train 3. 
0.95 
4 
0.18 
2 
0.73 
3 
9 

Train 4. 
0.889 
2 
0.31 
2 
0.74 
2 
6 

Test 1. 
0.652 
2 
0.31 
4 
0.55 
2 
8 

Test 2. 
0.991 
4 
0.61 
1 
0.53 
3 
8 

Test 3. 
0.601 
1 
0.43 
3 
0.11 
4 
8 

Test 4. 
0.881 
3 
0.56 
2 
0.6 
1 
6 
Table 8. Final ranking of each model using two optimization techniques
Optimization Techniques 
Simulation No. 
Total Score 
MLP 
1 
21 
2 
19 

3 
10 

4 
10 

Hybrid GAANN 
1 
16 
2 
19 

3 
17 

4 
12 
Based on the results of Tables (7) and (8), the first model in MLP simulation and the second model in GAANN 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 GAAAN simulation obtained the most score in other GAANN simulations with rank=19. In comparison between the best MLP and hybrid GAANN simulations, the high and reasonable R^{2} values between the model predictions and the measured data for training = 0.95 and testing =0.991 using the hybrid GAANN method describes its high capability in the prediction of diamond wire saw performance.
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 GAANN algorithm. This work compared the application and performance of hybrid GAANN 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 Fabrasivity, 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 GAANN algorithm is superior to MLP in terms of some model performance indices such as RMSE, R^{2} and STD. The comparison was made between the three simulations based upon the performance indices, hybrid GAANN algorithm with coefficient of determination (R^{2}) 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 GAANN 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 ICAANN, PSOANN, Hybrid Harmony Search (HSBFGS) and Grey Wolf Optimizer (GWO) and other performance indicators such as the mean absolute percentage error (MAPE) and value account for (VAF).
We would like to express our deepest thanks to Professor Mahdi Ghaem for his excellent advice.
[1] Ozcelik, Y., Polat, E., Bayram, F., & Ay, A. M. (2004). Investigation of the effects of textural properties on marble cutting with diamond wire. International Journal of Rock Mechanics and Mining Sciences, 41, 228234.
[2] Luo, S. Y., & Liao, Y. S. (1993). Effects of diamond grain characteristics on sawblade wear. International Journal of Machine Tools and Manufacture, 33(2), 257266.
[3] Luo, S. Y., & Liao, Y. S. (1995). Study of the behaviour of diamond sawblades in stone processing. Journal of materials processing technology, 51(1), 296308.
[4] Luo, S. Y. (1996). Characteristics of diamond sawblade wear in sawing. International Journal of Machine Tools and Manufacture, 36(6), 661672.
[5] Luo, S. Y. (1997). Investigation of the worn surfaces of diamond sawblades in sawing granite. Journal of materials processing technology, 70(1), 18.
[6] Xu, X. (2001). Study on the thermal wear of diamond segmented tools in circular sawing of granites. Tribology Letters, 10(4), 245250.
[7] Wei, X., Wang, C. Y., & Zhou, Z. H. (2003). Study on the fuzzy ranking of granite sawability. Journal of materials processing technology, 139(1), 277280.
[8] Ersoy, A., Buyuksagic, S., & Atici, U. (2005). Wear characteristics of circular diamond saws in the cutting of different hard abrasive rocks. Wear, 258(9), 14221436.
[9] Özçelik, Y., Kulaksız, S., (2000). Investigation of the relationship between cutting angles and wear on beads in diamond wire cutting method, 9th Mine Planning and Equipment Selection Symposium, Panagiotou & Michalakopoulos (Ed.), pp. 661666.
[10] Özçelik, Y., Kulaksız, S., & Çetin, M. C. (2002). Assessment of the wear of diamond beads in the cutting of different rock types by the ridge regression. Journal of Materials Processing Technology, 127(3), 392400.
[11] Özçelik, Y. (2003). Multivariate statistical analysis of the wear on diamond beads in the cutting of andesitic rocks. In Key Engineering Materials (Vol. 250, pp. 118130). Trans Tech Publications.
[12] Mikaeil, R., Ozcelik, Y., Ataei, M., & Shaffiee Haghshenas, S. (2016). Application of harmony search algorithm to evaluate performance of diamond wire saw. Journal of Mining and Environment.
[13] Almasi, S. N., Bagherpour, R., Mikaeil, R., & Ozcelik, Y. (2017). Analysis of bead wear in diamond wire sawing considering the rock properties and production rate. Bulletin of Engineering Geology and the Environment, 115.
[14] Mikaeil, R., Haghshenas, S.S., Ataei, M., Haghshenas, S.S., & Shafiee Haghshenas, A., (2017). Ranking and Assessment of Tunneling Projects Risks Using Fuzzy MCDM (Case Study: Toyserkan Doolayi Tunnel). 25th International Mining Congress and Exhibition of Turkey, P. 122128.
[15] Haghshenas, S. S., Neshaei, M. A. L., Pourkazem, P., & Haghshenas, S. S. (2016). The Risk Assessment of Dam Construction Projects Using Fuzzy TOPSIS (Case Study: Alavian Earth Dam). Civil Engineering Journal, 2(4), 158167.
[16] Haghshenas, S. S., Haghshenas, S. S., Barmal, M., & Farzan, N. (2016). Utilization of Soft Computing for Risk Assessment of a Tunneling Project Using Geological Units. Civil Engineering Journal, 2(7), 358364.
[17] Rad, M. Y., Haghshenas, S. S., & Haghshenas, S. S. (2014). Mechanostratigraphy of cretaceous rocks by fuzzy logic in East Arak, Iran. In The 4th International Workshop on Computer Science and EngineeringSummer, WCSE.
[18] Onyari, E. K., & Ilunga, F. M. (2013). Application of MLP neural network and M5P model tree in predicting streamflow: A case study of Luvuvhu catchment, South Africa. International Journal of Innovation, Management and Technology, 4(1), 11.
[19] Hornik, K., Stinchcombe, M., & White, H. (1989). Multilayer feedforward networks are universal approximators. Neural networks, 2(5), 359366.
[20] Daniels, H., & Kamp, B. (1999). Application of MLP networks to bond rating and house pricing. Neural Computing & Applications, 8(3), 226234.
[21] Khotanzad, A., & Chung, C. (1998). Application of multilayer perceptron neural networks to vision problems. Neural Computing & Applications, 7(3), 249259.
[22] Sonmez, H., Gokceoglu, C., Nefeslioglu, H. A., & Kayabasi, A. (2006). Estimation of rock modulus: for intact rocks with an artificial neural network and for rock masses with a new empirical equation. International Journal of Rock Mechanics and Mining Sciences, 43(2), 224235.
[23] HechtNielsen R (1987) Kolmogorov’s mapping neural network existence theorem. In: Proceedings of the first IEEE international conference on neural networks, San Diego, CA, pp 11–14.
[24] Hush DR (1989) Classification with neural networks: a performance analysis. In: Proceedings of the IEEE international conference on systems engineering. Dayton, OH, pp 277–280
[25] Kaastra, I., & Boyd, M. (1996). Designing a neural network for forecasting financial and economic time series. Neurocomputing, 10(3), 215236.
[26] Kanellopoulos, I., & Wilkinson, G. G. (1997). Strategies and best practice for neural network image classification. International Journal of Remote Sensing, 18(4), 711725.
[27] Ripley, B. D. (1993). Statistical aspects of neural networks. Networks and chaos—statistical and probabilistic aspects, 50, 40123.
[28] Paola, J. D. (1994). Neural network classification of multispectral imagery. Master Tezi, The University of Arizona, USA.
[29] Wang C (1994) A theory of generalization in learning machines with neural application. PhD thesis, The University of Pennsylvania.
[30] Masters T (1994) Practical neural network recipes in C++. Academic Press, Boston MA.
[31] Hagan, M. T., & Menhaj, M. B. (1994). Training feedforward networks with the Marquardt algorithm. IEEE transactions on Neural Networks, 5(6), 989993.
[32] Armaghani, D. J., Hasanipanah, M., & Mohamad, E. T. (2016). A combination of the ICAANN model to predict airoverpressure resulting from blasting. Engineering with Computers, 32(1), 155171.
[33] Zorlu, K., Gokceoglu, C., Ocakoglu, F., Nefeslioglu, H. A., & Acikalin, S. (2008). Prediction of uniaxial compressive strength of sandstones using petrographybased models. Engineering Geology, 96(3), 141158.
[34] Mikaeil, R., Haghshenas, S. S., Shirvand, Y., Hasanluy, M. V., & Roshanaei, V. (2016). Risk Assessment of Geological Hazards in a Tunneling Project Using Harmony Search Algorithm (Case Study: ArdabilMianeh Railway Tunnel). Civil Engineering Journal, 2(10), 546554.
[35] Mikaeil, R., Haghshenas, S. S., Haghshenas, S. S., & Ataei, M. (2016). Performance prediction of circular saw machine using imperialist competitive algorithm and fuzzy clustering technique. Neural Computing and Applications, 110.
[36] Haghshenas,S.S., Haghshenas, S.S., Mikaeil, R., Sirati Moghadam, P., & HaghshenasA.S., “A New Model for Evaluating the Geological Risk Based on Geomechanical Properties —Case Study: The Second Part of Emamzade Hashem Tunnel” Electronic Journal of Geotechnical Engineering, 2017. (22.01), pp 309320. Available at ejge.com.
[37] Holland, J. H. (1975). Adaptation in natural and artificial systems: an introductory analysis with applications to biology, control, and artificial intelligence. Ann Arbor, MI: University of Michigan Press, 441393.
[38] Terranova, O. G., Gariano, S. L., Iaquinta, P., & Iovine, G. G. (2015). GA SAKe: forecasting landslide activations by a geneticalgorithmsbased hydrological model. Geoscientific Model Development, 8(7), 19551978.
[39] Woodward, M., Kapelan, Z., & Gouldby, B. (2014). Adaptive flood risk management under climate change uncertainty using real options and optimization. Risk Analysis, 34(1), 7592.
[40] Saemi, M., Ahmadi, M., & Varjani, A. Y. (2007). Design of neural networks using genetic algorithm for the permeability estimation of the reservoir. Journal of Petroleum Science and Engineering, 59(1), 97105.
[41] Khandelwal, M., & Armaghani, D. J. (2016). Prediction of drillability of rocks with strength properties using a hybrid GAANN technique. Geotechnical and Geological Engineering, 34(2), 605620.
[42] Momeni, E., Nazir, R., Armaghani, D. J., & Maizir, H. (2014). Prediction of pile bearing capacity using a hybrid genetic algorithmbased ANN. Measurement, 57, 122131.
[43] Monjezi, M., Khoshalan, H. A., & Varjani, A. Y. (2012). Prediction of flyrock and backbreak in open pit blasting operation: a neurogenetic approach. Arabian Journal of Geosciences, 5(3), 441448.
[44] Aghajanloo, M. B., Sabziparvar, A. A., & Talaee, P. H. (2013). Artificial neural network genetic algorithm for estimation of crop evapotranspiration in a semiarid region of Iran. Neural Computing and Applications, 23(5), 13871393.