#### Abstract

Wire cut EDM is a quite regularly used machining process in mechanical and electronic industries. This research has attempted to machine aluminum alloy for which experimental design was prepared using Box-Behnken design. Different combinational options of pulse-on time (P_{1}), pulse-off time (P_{2}), servo wire feed (WF), and current (I) were investigated and surface roughness after machining was observed. Collected 27 datasets were further used in Adaptive Neuro Fuzzy Inference System (ANFIS) to produce about 500 datasets. These 500 datasets are approximated data derived from experimental datasets, known as synthetic data. Data model was further developed and used in Gorilla Troops Optimizer (GTO) to locate the optimum machining parameters. With the excellent three search operators: move towards other gorillas, migrate towards unknown places, and migrate towards known places, GTO has produced the lowest surface roughness value of 0.500953 *μ*m when the machining parameters of pulse-on time, pulse-off time, wire feed, and current values were set as 121 *μ*s, 52 *μ*s, 3 m/min, and 166A, respectively. To ensure the accuracy of the synthetic data-based model and optimality, verification and validation were conducted. Wilcoxon signed rank test was conducted for the pairwise comparison of GTO with each of its competing algorithms at the significance level of *σ* = 0.05. Friedman test was conducted to calculate the average ranking of each algorithm and to detect the global differences between all compared algorithms. Outperforming performance by GTO algorithm in machining of the selected material is found.

#### 1. Introduction

Among unconventional machining processes available, wire cut electric discharge machining (WEDM) is much effective process to prepare intricate shapes and delicate objects. This thermal erosion method uses the localized thermal energy to remove the material from the substrate metal. The supply of current and voltage to electrode produces plasma between the substrate metal and tool. A small amount of heat is transferred to the substrate metal to melt and evaporate. Besides, the electrolyte used in the process removes the machined metals away from the machine. Swift toggling of pulse-on time (P_{1}) and pulse-off time (P_{2}) causes fast heating and cooling effects over the surface. This advantage of WEDM has attracted this process to use for supper alloys, hard-cutting metals, and conductive materials. Khan et al. [1] investigated machining of ASTM A572-grade 50 steel, where grey rational analysis (GRA) was used to locate the optimal machining condition. Abbasi et al. [2] also attempted the same material and analyzed it to report that pulse-on time is the predominant parameter for surface finish. Takayama et al. [3] developed a high precision rotary table with many sensors to detect wire breakage, water leakage, and machining quality. The erosion of metal destabilizes at some depth of cutting due to accumulation of debris and it affects the process time. Ayesta et al. [4] used camera vision system to monitor the debris removal. They reported that debris removal is affected by the process parameters and hence optimal process parameters must be used not only to have high surface finish, but also to have high debris removal rate. Since electrolyte is important in this process, kerosene oil, EDM oil, and distilled water were also evaluated by hazard and operability analysis [5]. They revealed that aerosol concentration and dielectric consumption affect the environment and hence new dielectric fluid must be invented. Liu et al. [6] developed a pulse power source to avoid back-off error, which occurs in WEDM of semiconductors as discharge of silicon affects the servo controller. Newly developed power source has stabilized the discharge probability. Significant research contributions [7–25] in the recent years are reviewed and listed in Table 1.

Alam et al. [26] recently reviewed optimization of WEDM parameters, in which they detailed the past published research using GRA and RSM methods. Abidi et al. [27] used MOGA-II algorithm to optimize the parameters to machine nickel-titanium-based shape memory alloy. At this end, it is acknowledged that only minimal research has been conducted in finding the best machining condition using optimization algorithm.

Metaheuristic algorithms are more efficient in locating the optimal solutions for single or multiobjective problems [28, 29] and hence Gorilla Troops Optimizer (GTO) is attempted in this research to search for optimal solution in machining of Al alloy. In a brief of method of research, experiments were conducted using L_{27} orthogonal array. Surface roughness of machined samples and their respective cutting parameters were used to develop an Adaptive Neuro Fuzzy Inference System (ANFIS) from which about 500 synthetic data were further generated. Response surface model (RSM) was developed with these data and applied to GTO algorithm to locate the optimal process parameters. At this end, the paper is organized as below; Section 2 presents the material, experimentation, ANFIS model, and GTO algorithm, while Section 3 presents and discusses the results, while conclusion of current research is presented in Section 4.

#### 2. Materials and Methods

##### 2.1. Materials

Aluminum 7075 alloy was the material investigated in this research, as this material is the most suitable material for structural application including transport, aerospace, marine, and defence. Aluminum alloys are famous for their mechanical properties and weight-strength ratio. Machining this alloy using a conventional machining process does not result in the dimensional accuracy every time. Besides, rejection of parts is also there due to no align with near net shape. WEDM is one of the unconventional machining processes that perfectly suits this scenario and gives better surface finish. Mechanical properties and chemical properties of Al7075 alloy are listed in Tables 2 and 3, respectively.

##### 2.2. Experimentation

Al7075 alloy plate in 200 mm × 150 mm × 10 mm and WEDM SPRINTCUT model with CNC control was used for experimentation. Brass wire electrode of 0.25 mm diameter was used as cathode and deionised water was used as the dielectric fluid. Initially, the plate was machined using arbitrary chosen parameter setting to identify the sparking condition and wire breakage. After conducting a few preliminary experiments, the parameter setting range was fixed for pulse-on time (P_{1}), pulse-off time (P_{2}), servo wire feed (WF), and current (I) as shown in Table 4. Box-Behnken design (BBD) is the best method for threelevel design [25] and hence design of experiments (L_{27} DoE) was prepared using BBD. The machining was carried out as per DoE and surface roughness was measured in three different locations using surface roughness tester (Mitutoyo make Surftest 211). Figure 1 shows the machining and machined samples, while Table 5 reports the experimental results.

ANOVA, a statistical tool, was used to find the goodness of fit of the experimental data. The two-way ANOVA with 95% confidence level was carried out as the value of probability <5% is considered significant. Table 6 shows the results from ANOVA. Each parameter and cross product of parameters are found significant. The model is significant with *R*^{2} = 99.41% and standard deviation = 0.0201349.

##### 2.3. Adaptive Neurofuzzy Inference System (ANFIS)

ANFIS is a soft computing technique that utilizes artificial neural network (ANN) and fuzzy logic theory to predict the data samples. It is a simple and accurate adaptive technique to build a prediction model. The reason for using ANFIS in this research is that we wanted to generate as many datasets as possible from a small size of experimental data (27 datasets). Two-tier approach of using ANFIS and RSM was tested and reported as a better method [30]. Optimal searching was more accurate when data population was large. Hence, the above tabulated twenty-seven datasets were used in ANFIS for training the Neurofuzzy Interference System, which could result in plentiful datasets.

ANFIS has different types of membership function (MF) such as triangular shaped membership function (trimf), trapezoidal-shaped membership function (trapmf), generalized bell-shaped membership function (gbellmf), Gaussian curve membership function (gaussmf), Gaussian combination membership function (gauss2mf), P-shaped membership function (pimf), difference between two sigmoidal membership functions (dsigmf), and product of two sigmoid membership functions (psigmf). We attempted all these MFs to build predictive models and finally selected generalized bell-shaped membership function (gbellmf) based on minimum root mean square error (RMSE). This MF is Sugeno type IF-THEN rules to predict the outcomes. RMSE value observed in each model is listed as Trimf model = 1.1394, Trapmf model = 1.1696, Gbelmf model = 0.9605, Gausmf model = 0.9701, Gaus2mf model = 1.1677, Sigmf model = 0.832, Dsigmf model = 1.1672, Psigmf model = 1.1666, Pimf model = 1.7280, Smf model = 1.1728, and Zmf model = 1.1728.

Training of model was done with 27 experimental data. Six additional experiments were conducted with arbitrarily selected input parameters and these data as shown in Table 7 were used for testing the ANFIS model. Once training and testing of ANFIS model were completed, about 500 datasets were generated from the ANFIS model. These 500 datasets were used to build a response surface methodology (RSM) model shown below.

##### 2.4. Gorilla Troops Optimizer (GTO)

Gorilla Troops Optimizer (GTO) is one of the emerging metaheuristic search algorithms proposed by Abdollahzadeh et al. [31] in 2021 to solve global optimization problems, where its search mechanisms are essentially inspired by the collective behavior and social intelligence of gorilla life in nature. Similar to other species of apes, gorillas live in the troops consisting of an adult male gorilla known as silverback, multiple female and male gorillas, and offspring. Silverback refers to the strongest gorilla with unique silver-colored hair on its back and it is considered as troop leader that is responsible for decision making and conflict resolution, as well as spearheading other troop activities such as food searching and strategic retreat for troop safety. On the contrary, other male gorillas that do not have silver-colored hair on their back are relatively weaker and hence affiliated to silverback as the backup defenders of troops. Five typical behaviors are generally observed from gorilla troops to enable the formulation of search operators with different exploration and exploitation strengths. These behaviors include the moving to other gorillas, migration to unknown areas, migration to the surrounding of known areas, staying in original troop to follow silverback, and competition to mate with female gorillas upon the decease of silverback. Similar to other metaheuristic search algorithms, the optimization process using GTO can be expressed in three main stages, i.e., population initialization, exploration phase, and exploitation phase as explained in next subsections.

###### 2.4.1. Population Initialization

Suppose that a group of gorillas with population size of *N* is used to solve a given optimization problem with dimensional size of *D*. The position of each *i*-th gorilla in solution space at any *t*-th iteration is considered as a potential solution and expressed as , where and represent the population index of gorilla and dimensional index of problem, respectively. Let and be the lower and upper boundary limits of solution, respectively. At , the position value of every *i*-th gorilla, i.e., can be initialized aswhere is a real-valued number randomly generated from uniform distribution.

###### 2.4.2. Exploration Phase

From the perspective of metaheuristic search algorithms, an exploration process is essential in early stage of optimization to enable the population for discovering the unvisited regions in solution space. For the exploration phase of GTO, three search operators are incorporated by emulating different behaviors of gorilla troops, i.e., move towards other gorillas, migrate towards unknown places, and migrate towards known places.

Define as a constant parameter used to indicate the likelihood of *i*-th gorilla to migrate towards unknown places and as a real-valued number randomly generated from uniform distribution. If , the *i*-th gorilla chooses migration towards an unknown region to further enhance the exploration strength of GTO. For the scenario of , the *i*-th gorilla tends to move towards a randomly selected gorilla denoted as to attain good balancing of exploration and exploitation searches. Finally, the *i*-th gorilla chooses migration towards a known region to maintain good population diversity when . Let be candidate position of a randomly selected gorilla created in the *t*-th iteration, whereas refers to the candidate position of each *i*-th gorilla created in the next -th iteration. Define as two real-valued numbers randomly generated from uniform distribution; then

The parameters of , , and presented in equation (3) can be calculated as follows:where is a cosine operator; is the maximum iteration numbers set for GTO; and are the real-valued numbers randomly generated from uniform distribution. Equation (4) is used to simulate the tendency of GTO to perform searching with larger interval of changes in the earlier stage of optimization process and these search ranges are gradually reduced in later stage. Meanwhile, equation (5) is designed to emulate the variation of silverback’s leadership throughout the optimization process.

Upon the completion of exploration phase, the candidate position of each *i*-th gorilla is obtained as and its fitness value is evaluated as . The current fitness value of each *i*-th gorilla, i.e.,, is then compared with . The new candidate position is used to update the current position if the former solution has more superior fitness. Otherwise, the solution is discarded if it has more inferior fitness. Similar mechanism is used to update the position of silverback by comparing its fitness value of with that of .

###### 2.4.3. Exploitation Phase

On the contrary to exploration, exploitation has crucial role in the later stage of optimization process because it can further refine the solutions found in promising regions of search space. For the exploitation phase of GTO, two search operators are designed by simulating the behaviors of gorilla troops to follow the silverback for food searching and to compete for mating with the adult female gorillas.

Suppose that *W* is a parameter used to determine the behavior of each *i*-th gorilla in exploitation phase and *C* is a parameter obtained from equation (4). When , the *i*-th gorilla chooses to obey the instructions of silverback in exploitation phase and its candidate position at the next -th iteration; i.e., , is updated as follows:where *L* is a parameter obtained from equation (5); refers to the position of silverback at the *t*-th iteration. Meanwhile, the parameter *M* in equation (7) is calculated aswhere *N* refers to population size of GTO; represents the candidate position of each *i*-th gorilla at the current *t*-th iteration.

On the contrary, the *i*-th gorilla tends to compete with other male gorillas for mating with the adult female gorillas during exploitation phase if . The mathematical models used to simulate this behavior are formulated aswhere *Q* represents the impact force; *A* is a vector used to indicate the degree of violence encountered in competition. Both of *Q* and *A* can be computed aswhere are two real-valued numbers randomly generated from uniform distribution; is a constant; is used to simulate the violence effect on each dimension of solution based on two normal distributions of and If, is assigned as an array with size of, where is used to randomly generate different values for each dimension. On the other hand, only a single random number is generated by and assigned to when .

Similar with exploration phase, the fitness value of each *i*-th gorilla’s candidate position is evaluated as and then compared with the current fitness . If has more superior fitness than , the latter solution will be replaced by the former one. Otherwise, will be discarded due to its inferior fitness. Similar mechanism is also used to update the position of silverback by comparing its fitness value of with that of .

###### 2.4.4. Overall Search Mechanisms of GTO

The overall search mechanisms of GTO are presented in Figure 2. Initially, the population of GTO is randomly generated by using uniform distribution. Both of the exploration and exploitation phases are switched alternately during the optimization process of GTO. For exploration phase, three different search operators can be used to determine the candidate position of every *i*-th gorilla. Meanwhile, two different search operators are devised in exploitation phase to update the candidate position of every *i*-th gorilla. It is notable that greedy selection is used to update the current position of each *i*-th gorilla and position of silverback for both exploration and exploitation phases. The iterative search processes of GTO in exploration and exploitation phases are repeated until the predefined termination condition, i.e., current iteration exceeding maximum iteration (), is satisfied. At the end of optimization process, the silverback’s position, i.e., , is returned as the best solution to solve given problem. Due to the promising optimization performance of GTO, it has been applied to solve some real-world optimization problems such as feature selection [32], image segmentation [33], and parameter estimation of fuel cell [34] and solar panel [35]. To the best of authors’ knowledge, the feasibility of GTO to solve machining optimization problem is one of the areas that remains unexplored.

#### 3. Performance Evaluation of Proposed WEDM Machining Optimization Problem

##### 3.1. Simulation Settings of All Compared Algorithms

Extensive simulation studies are conducted to evaluate the performance of GTO in solving the proposed WEDM machining optimization problem. The simulation results produced by GTO are then compared with another six metaheuristic search algorithms known as particle swarm optimization (PSO) [36], differential evolution (DE) [37], teaching learning-based optimization (TLBO) [28], grey wolf optimizer (GWO) [38], sine cosine algorithm (SCA) [39], and arithmetic optimization algorithm (AOA) [40]. The optimal parameter settings of all compared algorithms are determined based on the recommendations in their respective literature as shown in Table 8 and original source codes are obtained from their authors to ensure fair comparison. The same population size of *N* = 10 and maximum iteration numbers of are also set to ensure all seven metaheuristic algorithms are compared in fair manner. When solving real-world optimization problem such as WEDM machining considered in current study, it is essential to achieve proper tradeoff between the solution accuracy and computational overhead of algorithm. Although it is more likely for metaheuristic search algorithm such as GTO to obtain the solutions with better accuracy when larger population size is set, the computational complexity of algorithm is estimated as *O* (*N* × (1 + + *D*) × 2) [31] using Big-O notation that tends to increase with population size, where *N*, , and *D* refer to the population size, maximum iteration numbers, and total dimensional size of problems, respectively. The increasing computational complexity of GTO is undesirable and not computationally feasible for real-world optimization problems. Furthermore, our simulation studies indicate that GTO can solve the WEDM machining optimization problems with better accuracy than other competing algorithms when the population size of *N* = 10 is set. This is a crucial finding to reveal the competitive optimization performances of GTO when solving real-world optimization problems. In order to alleviate the random discrepancy issue, all compared algorithms are simulated independently for 20 runs under the same simulation environment of Matlab (R2021a) installed in a personal computer with Intel® Core™ i5-7400 CPU processor with 24.0 GB RAM when solving the proposed WEDM machining optimization problems.

##### 3.2. Simulation Results and Discussion

The optimal machining parameters (i.e., pulse-on time, pulse-off time, wire feed, and current values) obtained by GTO and their corresponding response variable of surface roughness (*R*_{a}) when solving the proposed WEDM machining optimization problems for 20 consecutive times are presented in Table 9. The number of iteration numbers consumed by GTO to achieve the best machining result in every simulation run is also recorded in Table 9. Accordingly, GTO is able to consistently search for the optimal values of pulse-on time, pulse-off time, wire feed, and current that lead to the low values of surface roughness. It is also noteworthy that GTO has consumed an average iteration number of 24 to find the optimal machining parameters in these 20 simulations runs, implying the promising search efficiency of this algorithm. From Table 9, it is observed that GTO has produced the lowest surface roughness value of 0.500953 when the machining parameters of pulse-on time, pulse-off time, wire feed, and current values are set as 121 , 52, 3 m/min, and 166 A, respectively, at the 18th simulation run with 8 iterations and 19th simulation run with 15 iterations. These optimal machining parameters are further validated by measuring the experimental value of surface roughness as 0.500, i.e., only 0.2% of deviation from the predicted value of surface roughness. Given the small performance difference, it can be concluded that there is good consistency between the simulated and actual results.

Table 10 presents the best, worst, mean, and standard deviation (SD) values of surface roughness obtained by all compared seven metaheuristic search algorithms when solving the proposed WEDM machining optimization problems for 20 consecutive times. The best (i.e., lowest) surface roughness values produced by each algorithm are highlighted in bold. Although there are four algorithms (i.e., PSO, DE, TLBO, and GTO) that are able to produce the lowest values, GTO is the only algorithm that can solve the proposed WEDM machining optimization problems with the lowest values of worst and mean surface roughness. The promising performance of GTO to consistently solve the WEDM machining problem with low surface roughness value is also reflected from the lowest SD value. On the contrary, the optimization performances of algorithms such as PSO, DE, TLBO, and AOA are observed to be inconsistent as revealed by the inferior values of worst surface roughness and SD. These findings imply the high tendencies of PSO, DE, TLBO, and AOA to be trapped into local optima and suffer with premature convergence issues in certain trials of simulations. Although GWO and SCA can deliver more consistent performances when dealing with the proposed WEDM machining optimization problem, the overall surface roughness values obtained by these two competing algorithms are generally more inferior than those of GTO.

Figure 3 depicts the convergence curves obtained by all metaheuristic search algorithms when solving the proposed WEDM machining problem to further analyze their optimization performances in real-world application. Accordingly, each compared algorithm has exhibited different accuracy and efficiency in searching for the best combinations of machining parameters (i.e., pulse-on time, pulse-off time, wire feed, and current) that aim to minimize the surface roughness value. Among all compared metaheuristic search algorithms, AOA is identified to have worst performance because it has very slow convergence speed during the early stage of search process and its population only starts to converge towards the promising solution regions at the middle stage of optimization. The poor convergence characteristic of AOA is considered as a main contributing factor that prohibits its population to locate the global optimum in earlier stage, hence resulting in its inferior performance to solve the proposed WEDM machining problems with undesirably large surface roughness value. Although both PSO and DE have started the optimization processes with relatively lower surface roughness values, these two algorithms are observed to suffer with premature convergence issues as demonstrated by the plateau regions of their convergence curves in early stage of optimization, i.e., before 20th iteration. The high tendency of PSO and DE populations to be trapped into the local optima regions can be justified by their relatively simple search mechanisms that are unable to achieve good balancing of exploration and exploitation when dealing with fitness landscapes of real-world optimization problems. The inferior performances of PSO and DE to avoid misleading information of local optima regions tend to prevent these two algorithms consistently solving the proposed WEDM machining problems with lower surface roughness values. The remaining algorithms such as TLBO, GWO, SCA, and GTO have better robustness to handle premature convergence issues as demonstrated by their convergence curves with more promising convergence characteristics and lower surface roughness values. Among these four algorithms, GTO is proven as the best performing algorithm to solve the proposed WEDM machining optimization problem because it shows fastest convergence rate and is able to locate the solution that is nearest to the global optimum. The competitive optimization performances demonstrated by GTO against TLBO, GWO, and SCA can be justified by its inherent mechanisms used to achieve proper balancing of exploration and exploitation searches of algorithm. As compared to the latter three algorithms equipped with lesser numbers of search operators, GTO is designed to have two optimization phases (i.e., exploration and exploitation) with five search operators with different exploration and exploitation strengths. Depending on the optimization stage and location of each individual solution in search space, one of the embedded search operators can be triggered and play dominant role to guide the individual solution searching towards the promising regions of search space. The excellent capability of GTO in balancing its exploration and exploitation searches enables its population to locate the global optimum of WEDM machining optimization problem rapidly without requiring large numbers of population size for effective searching.

Different nonparametric statistical procedures [41, 42] are further used for performance evaluation of all competing metaheuristic search algorithms by referring to their surface roughness values obtained. Wilcoxon signed rank test [41, 42] is employed for the pairwise comparison of GTO with each of its competing algorithms at the significance level of *σ* = 0.05. The results of Wilcoxon signed rank test are reported in terms of the sum of rank (i.e., and) and the associated *p*-values. and shown in Table 11 indicate the sum of rank for GTO to perform better and perform worse than a given competing algorithm, respectively. Meanwhile, *p*-value is the minimum level of significance for detecting performance differences between algorithms. If , it implies the existence of strong evidence to reject null hypothesis and the best results achieved by better performing algorithms are statistically significant. The pairwise comparison results reported in Table 12 show that GTO can perform significantly better than PSO, GWO, SCA, and AOA because their corresponding *p*-values are all smaller than the threshold significance value of *σ* = 0.05. In other words, the competitive performances of GTO against PSO, GWO, SCA, and AOA are evident and not achieved by any random chances. On the other hand, there are no significance performance differences detected between GTO and DE as well as GTO and TLBO as indicated by their relatively large *p*-values. This implies that the poor performances of DE and TLBO might only happen in certain simulation runs. GTO is still considered as a better performing algorithm due to excellent consistency to solve the proposed WEDM machining optimization problem with low surface roughness values.

Multiple comparison analyses [41, 42] are also conducted for more thorough evaluations of GTO and its competing algorithms. Friedman test was first employed to calculate the average ranking of each algorithm and detect the global differences between all compared algorithms. As reported in Table 12, GTO with lowest average rank value emerges as the best performing algorithm to solve the proposed WEDM problem with lowest surface roughness values. The remaining compared algorithms are ranked as follows: TLBO, DE, SCA, GWO, PSO, and AOA. The *p*-value presented in Table 12 also detects the significant global differences between all compared algorithms because , where *σ* = 0.05 is defined as the threshold significance level. Referring to the findings from Friedman test, another three post hoc analyses known as Bonferroni-Dunn, Holm, and Hochberg methods are employed to further compare the performance differences between GTO and each peer algorithm [41, 42]. Table 13 presents the results of each post hoc analysis in terms of *z* values, unadjusted *p*-values, and adjusted *p*-values. It is observed that all post hoc analyses have validated that GTO can outperform AOA and PSO significantly in terms of surface roughness value because all adjusted *p*-values obtained are smaller than *σ* = 0.05. If the threshold significance level is adjusted to *σ* = 0.10, GTO is verified by all post hoc analyses to perform significantly better than GWO. Finally, both of Holm and Hochberg analyses can verify the significant performance of GTO against SCA when *σ* = 0.10.

#### 4. Conclusion

This research was aimed at using an emerging metaheuristic algorithm known as Gorilla Troops Optimizer (GTO) to solve WEDM machining problem by searching for the optimal machining parameters that can lead to the best (i.e., lowest) surface finish.

The machining of aluminum alloy was done with brass wire as a cathode tool and deionised water as a dielectric medium. ANFIS was used to generate about 500 datasets from 27 experimental datasets. On the contrary to many existing metaheuristic search algorithms, GTO has two different optimization phases with five search operators. It was able to search for the optimal machining parameters in solution space effectively and efficiently without having large population number. At the end of optimization processes, the optimal machining parameters as stored in silverback solution of GTO were pulse-on time = 121 *μ*s, pulse-off time = 52 *μ*s, wire feed = 3 m/min, and current value = 166 A. These parameters are the best process parameters to achieve the minimum surface roughness of 0.500953 microns.

The current metaheuristic algorithm was able to solve WEDM machining problem in 20 independent simulation runs with small standard deviation value of *SD* = 0.0060295, implying the high consistency of the algorithm to search for the optimal machining parameters. Various statistical analyses such as Wilcoxon signed rank test, Friedman test, and post hoc analyses were also performed to verify the significant performance gain of GTO against other metaheuristic search algorithms when solving the current machining problem. To ensure the practicability of optimization results produced by GTO, validation experiments were also performed and found deviation of ∆ = 0.2% between the predicted and measured surface roughness values. The negligible error value implies the feasibility of GTO to solve real-world WEDM machining optimization problems.

#### Nomenclature

ANFIS: | Adaptive Neurofuzzy Inference System |

WEDM: | Wire cut electric discharge machining |

GRA: | Grey rational analysis |

HSS: | High speed steel |

MQL: | Minimum quantity lubrication |

MRR: | Material removal rate |

EWR: | Electrode wear rate |

OC: | Overcut |

RSM: | Response surface model |

ARAS: | Additive ratio assessment |

AHP: | Analytic hierarchy process |

TOPSIS: | Technique for order preference by similarity to ideal situation |

BBD: | Box-Behnken design |

MOPSO: | Multiobjective particle swarm optimization |

GTO: | Gorilla Troops Optimizer |

DoE: | Design of experiment |

ANN: | Artificial neural network |

MF: | Membership function |

trimf: | Triangular shaped membership function |

trapmf: | Trapezoidal-shaped membership function |

gbellmf: | Generalized bell-shaped membership function |

gaussmf: | Gaussian curve membership function |

gauss2mf: | Gaussian combination membership function |

pimf: | P-shaped membership function |

dsigmf: | Difference between two sigmoidal membership functions |

psigmf: | Product of two sigmoid membership functions |

gbellmf: | Generalized bell-shaped membership function |

RMSE: | Root mean square error |

PSO: | Particle swarm optimization |

DE: | Differential evolution |

TLBO: | Teaching learning-based optimization |

GWO: | Grey wolf optimizer |

SCA: | Sine cosine algorithm |

AOA: | Arithmetic optimization algorithm |

#### Indices

i: | Population index of gorilla in GTO |

d: | Dimension index |

t: | Iteration index |

: | Cosine operator |

: | Operator used to calculate the fitness value of a given position value |

O : | Complexity of algorithm estimated using Big-O notation |

#### Parameter and Variables

P_{1}: | Pulse-on time |

P_{2}: | Pulse-off time |

WF: | Servo wire feed |

I: | Current |

R _{a}: | Surface roughness |

D: | Total dimension size of optimization problem |

N: | Population size |

: | Position value of ith gorilla at the tth iteration |

: | Position value of ith gorilla in dth dimension at the tth iteration |

, , , ,, , : | Real-valued number within [0, 1] randomly generated from uniform distribution |

: | Real-valued number within [-1, 1] randomly generated from uniform distribution. |

: | Lower boundary limits |

: | Upper boundary limits |

p: | Constant parameter within [0, 1] used to indicate the likelihood of gorilla to migrate towards unknown places |

: | Position value of a randomly selected gorilla at the tth iteration |

: | Candidate position value of a randomly selected gorilla at the tth iteration |

: | Candidate position value of ith gorilla at the tth iteration |

: | Search range of gorilla |

L: | Variation of silverback’s leadership |

: | Maximum iteration numbers |

W: | Parameter used to determine the behavior of gorilla in exploitation phase |

M: | Mean position of gorilla population |

Q: | Impact force |

A: | Degree of violence encountered by gorilla in competition |

: | Constant parameter used to determine violence effect of gorilla |

: | Parameter used to simulate the violence effect on each dimension of solution |

, : | Normal distributions used to model the violence effect |

: | Inertia weight of particle swarm optimization |

: | Acceleration coefficient of particle swarm optimization |

: | Scaling factor of differential evolution |

: | Crossover rate of differential evolution |

: | Teaching factor of teaching learning based optimization |

: | Convergence coefficient of grey wolf optimizer |

: | Constant used to balance exploration and exploitation of sine cosine algorithm |

: | Control parameter used to adjust the search process of arithmetic optimization algorithm |

: | Control parameter used to define exploitation accuracy of arithmetic optimization algorithm over iterations |

: | Iteration number consumed by an algorithm to achieve best result |

SD: | Standard deviation |

σ: | Significant level of nonparametric statistical procedures |

: | Sum of rank for GTO to perform better than a given competing algorithm |

: | Sum of rank for GTO to perform worse than a given competing algorithm. |

#### Data Availability

The data associated with this research can be obtained from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest.