Multi-Objective Optimization of Drilling GFRP Composites Using ANN Enhanced by Particle Swarm Algorithm

: This paper aims to optimize the quality characteristics of the drilling process in glass ﬁber-reinforced polymer (GFRP) composites. It focuses on optimizing the drilling parameters with drill point angles concerning delamination damage and energy consumption, simultaneously. The effects of drilling process parameters on machinability were analyzed by evaluating the machinability characteristics. The cutting power was modeled through drilling parameters (speed and feed), drill point angle, and laminate thickness. The response surface analysis and artiﬁcial neural networks enhanced by the particle swarm optimization algorithm were applied for modeling and evaluating the effect of process parameters on the machinability of the drilling process. The most inﬂuential parameters on machinability properties and delamination were determined by analysis of variance (ANOVA). A multi-response optimization was performed to optimize drilling process parameters for sustainable drilling quality characteristics. The obtained models were applied to predict drilling process characteristics, and exhibited an excellent harmony with the experiment results. The optimal drilling process factors were the highest spindle speed and the lowest feed, with a drill point angle of 118 ◦ for the laminate of 4.75 mm thickness.


Introduction
Fiber-reinforced polymer composite materials attract great interest in the industry because of their properties that make them superior to other traditional materials in the manufacture of structures for aircraft, cars, trains, ships, and others [1,2]. These structures are usually assembled by screw connections or by any other means, through the holes drilled in them. The quality of these holes is critical, as non-matching of quality characteristics such as roundness, roughness, straightness, parallelism, and perpendicularity can create excessive stress on the screws, leading to their failure.
Many researchers have reported their studies on the impact of drilling process parameters such as cutting speed, feed rate, and geometry of drilling tools on machining process characteristics such as critical thrust force, drilling torque, induced temperature, and damage in laminates of fiber-reinforced polymer composites like Glass-FRP and Carbon-FRP. The delamination damage is more intricate, due to the need for systematic repair increasing manufacturing costs, since repairing a defective hole caused by such delamination takes 5-6 h per hole [3,4]. Moreover, Davim and Reis [5] investigated the effect of the machining parameters on the cutting power in drilling carbon fiber-reinforced plastic composite experimentally.
Manufacturing industries require large inputs of energy, while sometimes the manufacturing industry's participation in the energy consumed in the industry is greater than that for its economic return from the industry [6]. Regarding sustainable machining, Ali Bastas [7] examined 46 research articles targeting sustainable manufacturing, where it was noted that turning and milling are the most common manufacturing processes in techniques based on artificial intelligence and machine learning can be used as an effective alternative to overcome these obstacles [41].
The optimization process aims to find the optimal value or the best solution. Optimization problems involve finding the maximum or minimum value using a single or multiple objective. Multi-objective optimization refers to finding the optimal solution values for more than one desired objective. Multi-objective optimization became crucial as a problem-solving tool in industrial processes [42]. The controlled elitist multi-objective genetic algorithm was applied for optimizing cell size and density to produce foam with the best mechanical properties [43], in addition to using the finite element method in analyzing and improving the characteristics of operating processes [44].
Machining processes are responsible for a big part of the manufacturing environmental footprint, due to their power consumption. The cutting power could be minimized by optimizing the cutting condition in machining processes [45]. Thus, manufacturers can save money and gain safer sustainability by minimizing their power demands to a minimal level. Moreover, the quality of a product achieved by minimizing the delamination by controlling the drilling parameters will be improved, and reflects economic sustainability. To the author's knowledge and in line with the literature review, an investigation of optimizing the drilling process of woven GFRP composites with a high-accuracy CNC machine by minimizing cutting power consumption with minimal delamination damage has not been addressed by controlling drilling parameters (feed and speed) with various drill point angles and for different thickness of laminates.
Statistical methods were used in analyzing the results of the full factorial experiments and extracting the regression models for process outcomes. The proposed machine learning model PSO-NN was trained and tested using the experimental data with an optimizing training process using the PSO algorithm. Different measures were used to evaluate the prediction accuracy of the models obtained.
In the end, the drilling operations were optimized in order to minimize the level of delamination and the power consumption in cutting with the lowest thrust force and maximum drilling torque by controlling the drilling operation parameters. The optimal operating parameters of the feed, spindle speed, angular drilling point, and lamination thickness of the woven GFRE material are thus determined.

Material and Methods
The samples used in the experiments were 3 laminates of woven GFRP fabricated by the hand lay-up technique, produced in different thicknesses of 2.6 mm, 5.3 mm, and 7.7 mm. The fiber volume fraction was 40% in the matrix fabricated from Araldite LY5138-2 epoxy polymer with the use of HY5138 Hardener. The specimens were prepared by cutting to standard dimensions using an abrasive water-jet CNC machine [46]. The tensile properties of the samples such as tensile strength, modulus of elasticity, and Poisson's ratio were 203.86 MPa, 16.05 GPa, and 0.295, respectively.
Full factorial design was used to plan experiments for four control factors (spindle speed, feed, drill point angle, and laminate thickness) within the given levels shown in Table 1. The experimental results extracted from [46] were used in the investigation of the effect of the mentioned machining parameters on the machinability proprieties, drilled hole quality, and cutting-power consumption. The delamination factors could be assessed by applying several methods [47]. In the current work, the peel-up and push-out surface delamination were assessed using the AutoCAD technique, which is appropriate for such materials under study, as they are quasi-transparent composite materials. This technique is reported in detail by [48,49].
The delamination factor F d was evaluated by where D o is the nominal hole diameter and D max is the maximum diameter of delamination that can be drawn from the center point of the D o.
The consumed cutting power shown in Table 2 was computed by applying the maximum drilling torque, measured by the Kistler dynamometer, in Equation (2).
where T is the torque (N.m), and s is the spindle speed (r/min).

PSO-NN Model
ANNs are known as an arrangement of mathematical algorithms mimicking the human brain in learning by applying simple computations. An ANN is a vigorous tool used in modeling non-linear relationships between the inputs and output(s) of the process [50], and it is one of the means of carrying out machine learning. As a most common type of feed-forward network, a multi-layer perceptron is used. It consists of three types of layers, the input, hidden, and output layers. Every layer consists of neurons containing weights, and biases. Activation functions such as linear, sigmoid, and logistic are used for converting inputs into outputs, as well as the learning algorithm [51].
To develop and apply ANN, three sequential processes are required, which are training, validation, and testing. Training the networks aims to find the optimum values of weights and biases through the network's layers. The objective function used to find the optimum training result of ANN is minimizing mean square error (MSE), using Equation (3), where I is the measured value as taken in the training data set and O is the ANN output, and n is the number of examples used in the training process. The best trained network is selected for the minimal MSE between the measured and predicted responses. For evaluating the performance of the best-obtained network, the correlation coefficient (R) is applied between the actual and the predicted values.
PSO as a metaheuristic improvement technique simulates the behavior of birds while searching for food [52,53]; the PSO algorithm was simplified, and it was observed that it leads to improvement, by applying the intelligent swarm method.
Each particle stores in its own memory its best position, along with the global best position. Given that information and the particle's velocity v i , the new particle's position becomes x k+1 i . Inside the range of the search domain with upper (U) and lower (L) boundaries, the population set of moving particles is generated using Equation (4).
where w is the weight of inertia, c 1 and c 2 are the acceleration coefficients, r 1 and r 2 are random numbers, pBest and gBest are best particle position and best global solution, respectively, and the iteration number is k. Then, determine the new particle's position by applying Equation (6).
And then compute the fitness of each particle. After that, determine the pBest and the gBest. Repeat these steps until the preset parameter value is reached.
NN training is a miniaturization problem in which the network weights and biases are determined, to achieve minimal prediction error. Enhancing neural networks with methods such as metaheuristic optimization is often required to avoid the hurdles that classical training algorithms, face such as the possibility of long training time or deception at local minimum points, and thus improve the prediction performance of these integrated networks. The neural network is enhanced by the PSO algorithm. The PSO-NN is here to utilize the PSO algorithm in the optimization of the neural network performance by minimizing the prediction errors. That is performed by identifying the optimal both weight and bias values of the net. The minimal MSE expressed by Equation (2) is the objective function of the optimum performance of the trained neural network applying the PSO algorithm, which updates the neural network even until reaching the stop limit.

Results and Discussion
In this section, the focus is on analyzing the effect of processing parameters on the cutting power consumed, and comparing it with the torque generated by drilling, which has not been addressed in previous research. In addition, we present the performance results of predictive models using neural networks and comparing them with multiple regression models. Then the cutting power consumed is included as a new dimension in the multi-objective optimization process, the focus of this work.

Effect of Drilling Parameters on the Torque
It is evident from the torque main effects plot shown in Figure 1 that the torque is significantly increased by increasing the feed, due to the increase in the cross-sectional area of the uncut chip. It is decreased by 15% by increasing the spindle speed from 400 r/min to 1600 r/min. NN training is a miniaturization problem in which the network weights and biases are determined, to achieve minimal prediction error. Enhancing neural networks with methods such as metaheuristic optimization is often required to avoid the hurdles that classical training algorithms, face such as the possibility of long training time or deception at local minimum points, and thus improve the prediction performance of these integrated networks. The neural network is enhanced by the PSO algorithm. The PSO-NN is here to utilize the PSO algorithm in the optimization of the neural network performance by minimizing the prediction errors. That is performed by identifying the optimal both weight and bias values of the net. The minimal MSE expressed by Equation (2) is the objective function of the optimum performance of the trained neural network applying the PSO algorithm, which updates the neural network even until reaching the stop limit.

Results and Discussion
In this section, the focus is on analyzing the effect of processing parameters on the cutting power consumed, and comparing it with the torque generated by drilling, which has not been addressed in previous research. In addition, we present the performance results of predictive models using neural networks and comparing them with multiple regression models. Then the cutting power consumed is included as a new dimension in the multi-objective optimization process, the focus of this work.

Effect of Drilling Parameters on the Torque
It is evident from the torque main effects plot shown in Figure 1 that the torque is significantly increased by increasing the feed, due to the increase in the cross-sectional area of the uncut chip. It is decreased by 15% by increasing the spindle speed from 400 r/min to 1600 r/min.
Meanwhile, the drill point angle influences the generated torque differently. The torque produced with the use of drill point angles 100° and 140° is higher than at drill point angle 118°. It means that the highest average torque is in the cases at the lower point angle (100°) and the lowest average torque is at the angle 118°. It is evident from Figure 1 that the effect of thickness on torque contrasts with that of the drill point angle. As shown, the thickness increase from 2.6 to 5.3 mm was accompanied by an increase in torque of ~45%, but as the thickness increased from 5.3 to 7.7 mm, the torque slightly decreased.

Effect of Drilling Parameters on Power Consumption
Ali Bastas, 2021 [7] examined 46 research articles targeting sustainable manufacturing, where it was noted that turning and milling are the most common manufacturing processes in sustainable manufacturing technology research, with 67% of the total examined articles. A total of 17% of research articles studied cutting condition optimization, in which the optimal parameters were investigated, including speed, feed, depth of cut, and material removal rates for improving sustainability, including optimizing power Meanwhile, the drill point angle influences the generated torque differently. The torque produced with the use of drill point angles 100 • and 140 • is higher than at drill point angle 118 • . It means that the highest average torque is in the cases at the lower point angle (100 • ) and the lowest average torque is at the angle 118 • . It is evident from Figure 1 that the effect of thickness on torque contrasts with that of the drill point angle. As shown, the thickness increase from 2.6 to 5.3 mm was accompanied by an increase in torque of 45%, but as the thickness increased from 5.3 to 7.7 mm, the torque slightly decreased.

Effect of Drilling Parameters on Power Consumption
Ali Bastas, 2021 [7] examined 46 research articles targeting sustainable manufacturing, where it was noted that turning and milling are the most common manufacturing processes in sustainable manufacturing technology research, with 67% of the total examined articles. A total of 17% of research articles studied cutting condition optimization, in which the optimal parameters were investigated, including speed, feed, depth of cut, and material removal rates for improving sustainability, including optimizing power consumption. In general, energy consumption is reduced by reducing the thrust force and cutting torque [54]. Figure 2 shows a graph of the main effects of the power consumed, and it is observed from the curves represented in the figure that the behavior of the power function and torque function in Figure 1 are similar concerning the feed, the angle of the drilling point and the thickness of the laminate. This is due to the mathematical Equation (2) used to calculate the cutting power (P) for the drilling process, where the cutting power is proportional to the generated torque. The cutting power is significantly increased by increasing the feed. The computed cutting power at the use-of-drill point angles 100 • and 140 • is higher than at the use-of-drill point angle 118 • . It means that the highest cutting power is at the case of the point angle of 100 • and the lowest power is at the case of the point angle of 118 • , while it is evident that the computed cutting power is linearly proportional to cutting speed, based on the formula used in computing the power. It can be observed that the lowest value of cutting power is generated at the lowest values of speed and feed. Similar results were reported in [55]. consumption. In general, energy consumption is reduced by reducing the thrust force and cutting torque [54]. Figure 2 shows a graph of the main effects of the power consumed, and it is observed from the curves represented in the figure that the behavior of the power function and torque function in Figure 1 are similar concerning the feed, the angle of the drilling point and the thickness of the laminate. This is due to the mathematical Equation (2) used to calculate the cutting power (P) for the drilling process, where the cutting power is proportional to the generated torque. The cutting power is significantly increased by increasing the feed. The computed cutting power at the use-of-drill point angles 100° and 140° is higher than at the use-of-drill point angle 118°. It means that the highest cutting power is at the case of the point angle of 100° and the lowest power is at the case of the point angle of 118°, while it is evident that the computed cutting power is linearly proportional to cutting speed, based on the formula used in computing the power. It can be observed that the lowest value of cutting power is generated at the lowest values of speed and feed. Similar results were reported in [55].
Similar to torque, power consumption increases with more feed, while the power consumption differs from the torque with the increase in the spindle speed. In other words, as the spindle speed increases, the power consumption increases and the torque tends to decrease slightly.
Regression analysis and factorial interaction plots are used to examine the interaction between the torque and spindle speed, which is illustrated in Figure 3.  Delamination of FRP composites affects the safety and reliability of any built structure. Therefore, the delamination induced during the drilling process has a significant economic impact, especially keeping in mind the different stages through which a Similar to torque, power consumption increases with more feed, while the power consumption differs from the torque with the increase in the spindle speed. In other words, as the spindle speed increases, the power consumption increases and the torque tends to decrease slightly.
Regression analysis and factorial interaction plots are used to examine the interaction between the torque and spindle speed, which is illustrated in Figure 3.
Delamination of FRP composites affects the safety and reliability of any built structure. Therefore, the delamination induced during the drilling process has a significant economic impact, especially keeping in mind the different stages through which a component goes when it reaches the final assembly line. These defects that occurred during manufacturing could mean the loss of the sustainable advantage of FRP composites. This delamination exhibited a convoluted failure mode that was produced by a combination of thermal and mechanical damage. The delamination may be the entry (peel-up) or exit (push-down) of the composite laminate. Entry delamination occurs due to fracture modes II and III. Fractures of mode II occurs when the fibers of the upper layers are not sufficiently cut. Fractures of mode III occur when the twisted drill edges face the laminates and a peeling force is generated by the slope of the drill flute that separates the upper layers [4]. The bottom-surface exit delamination occurs due to both mode I and mode II of fracture, and is attributed to the exposure of the material to axial forces and bending forces. The previous literature showed that exit delamination is larger and more critical than entry delamination, due to the absence of backup support, which can neutralize thrust force during the permeation of the drill [44,45].
words, as the spindle speed increases, the power consumption increases and the torque tends to decrease slightly.
Regression analysis and factorial interaction plots are used to examine the interaction between the torque and spindle speed, which is illustrated in Figure 3.  Delamination of FRP composites affects the safety and reliability of any built structure. Therefore, the delamination induced during the drilling process has a significant economic impact, especially keeping in mind the different stages through which a  This delamination exhibited a convoluted failure mode that was produced by a combination of thermal and mechanical damage. The delamination may be the entry (peel-up) or exit (push-down) of the composite laminate. Entry delamination occurs due to fracture modes II and III. Fractures of mode II occurs when the fibers of the upper layers are not sufficiently cut. Fractures of mode III occur when the twisted drill edges face the laminates and a peeling force is generated by the slope of the drill flute that separates the upper layers [4]. The bottom-surface exit delamination occurs due to both mode I and mode II of fracture, and is attributed to the exposure of the material to axial forces and bending forces. The previous literature showed that exit delamination is larger and more critical than entry delamination, due to the absence of backup support, which can neutralize thrust force during the permeation of the drill [44,45]. Figure 4 shows the main effect plot for Fd-exit. The Fd-exit factor is increased significantly by increasing feed, by reason of an increasing thrust force. Also, Fd-exit is increased by increasing the drill point angle. Despite the thrust force associated with the 5.3 mm thickness being higher than that of the 2.6 mm laminate, the Fd-exit of the 2.6 mm laminate is higher than those of the 5.3 mm one. This result is attributed to the poor stiffness of the thinnest laminate, leading to a bending in the last layer larger than those in the thickness of 5.3 mm. In the reverse direction, the Fd-exit of the thicker specimen of 7.7 mm is larger than those of 5.3 mm. This result is due to the accumulation of thermal and mechanical damages together in the 7.7 mm specimen, and with the highest drilling temperature, as detailed in [43]. The effect of the feed rate may be due to the increase in chip area accompanying the increase in the feed rate, which increases the material's resistance to cutting and thus the effort expended to remove the chip [56].

ANOVA
ANOVA tables are summarized in Table 3, as well as the R 2 (coefficient of determination) value of each model. From ANOVA results, the influencing control factors and all terms of RSA models on the process outputs are evaluated by their contribution percentage in variation. The p-value infers the significance value of the model term's effect on the responses if its value is less than 0.05, as in most experimental work. Consequently, the process parameters and other control factors studied influence the drilling process outcomes, as illustrated in Table 3.
ANOVA analysis showed that feed rate has the greatest influence and participation in building the mathematical models obtained for the drilling process in glass fiber-reinforced polymer composites. The contribution of each main effect factor which resulted Despite the thrust force associated with the 5.3 mm thickness being higher than that of the 2.6 mm laminate, the Fd-exit of the 2.6 mm laminate is higher than those of the 5.3 mm one. This result is attributed to the poor stiffness of the thinnest laminate, leading to a bending in the last layer larger than those in the thickness of 5.3 mm. In the reverse direction, the Fd-exit of the thicker specimen of 7.7 mm is larger than those of 5.3 mm. This result is due to the accumulation of thermal and mechanical damages together in the 7.7 mm specimen, and with the highest drilling temperature, as detailed in [43]. The effect of the feed rate may be due to the increase in chip area accompanying the increase in the feed rate, which increases the material's resistance to cutting and thus the effort expended to remove the chip [56].

ANOVA
ANOVA tables are summarized in Table 3, as well as the R 2 (coefficient of determination) value of each model. From ANOVA results, the influencing control factors and all terms of RSA models on the process outputs are evaluated by their contribution percentage in variation. The p-value infers the significance value of the model term's effect on the responses if its value is less than 0.05, as in most experimental work. Consequently, the process parameters and other control factors studied influence the drilling process outcomes, as illustrated in Table 3. ANOVA analysis showed that feed rate has the greatest influence and participation in building the mathematical models obtained for the drilling process in glass fiber-reinforced polymer composites. The contribution of each main effect factor which resulted from the ANOVA of the general linear model of each response is shown in Figure 5. The feed participations were 84.6%, 43.5%, 15%, and 44.4% for thrust force, torque, power consumption, and delamination factor, respectively. The impact of the feed on power consumption is due to its proven effect on the drilling torque, which in turn is a key factor in computing the power consumed.
Processes 2023, 11, x FOR PEER REVIEW 9 of 18 consumption, and delamination factor, respectively. The impact of the feed on power consumption is due to its proven effect on the drilling torque, which in turn is a key factor in computing the power consumed. From all of the above, and in order to achieve a sustainable drilling process that takes into account the quality of the holes produced, represented by the delamination factor, and that it keeps the cost of the energy consumed low, controlling the feed rate and the cutting speed simultaneously is necessary to achieve this goal.  Figure 5. Contribution of main factors to the drilling process attributes.
From all of the above, and in order to achieve a sustainable drilling process that takes into account the quality of the holes produced, represented by the delamination factor, and that it keeps the cost of the energy consumed low, controlling the feed rate and the cutting speed simultaneously is necessary to achieve this goal.
The correlation between drilling process responses and the different process control factors was evaluated utilizing RSA, which gave the mathematical equations of thrust force, torque, and delamination factor mentioned below. The significance index showed that the models obtained were statistically significant in analyzing each response, with a confidence level of 95% with acceptable goodness. As the parameters were in multiple levels, the used quadratic model was expressed as follows: where Y refers to response, f, s, a and t are the control factors of experimental design, and the β , s are the variation resource coefficients of the predictive model. The quadratic equations of thrust force, torque, and delamination factor are expressed by Equations (7)- (9).
The surface response for the drilling process properties can be plotted, and can predict the response values for any combination of the drilling parameters within the space of the designed experiment. In a later section, these predictive models will be compared to other models obtained using neural networks.

Predictive PSO-NN Models
The experimental data are used in training the proposed PSO-NN model. The PSO-NN model learns to estimate the relationship between drilling parameters (feed, cutting speed, drill point angle, and laminate thickness) and various responses (F t , T, P, and F d ).
The data used in training are divided into the training dataset (70%) and the test dataset (30%). This distribution is applied for the thrust force, torque, cutting power, and exit delamination factor, and then the training is performed, utilizing the PSO algorithm to obtain optimal parameters for the neural network, which minimizes prediction errors. Increasing the number of neurons in the hidden layer reduces the error while taking care not to overfit [50].
The proposed PSO-NN structures are four inputs and eight neurons for one hidden layer, C 1 and C 2 are 1.4, and 2.5, respectively, and the search boundaries LB and UB are −1.273 and 1.273, respectively. Table 4 elucidates the performance of trained PSO-NN models, which are measured by MSE, along with the correlation coefficient (R) between predicted and measured values.
It is evident that the predicted data from trained PSO-NN models coincide with the experimental results in relation to all drilling properties under study. This is also illustrated by the strong correlation between measured and predicted responses for the four drilling properties represented by the scatter diagram shown in Figure 6. It is evident that the predicted data from trained PSO-NN models coincide with the experimental results in relation to all drilling properties under study. This is also illustrated by the strong correlation between measured and predicted responses for the four drilling properties represented by the scatter diagram shown in Figure 6. Figure 6a,b,d of the Thrust force, Torque, and Power PSO-NN models show high accuracy of the models, where the predicted values appear closely scattered around the straight line, while Figure 6c shows the accuracy of the Fd PSO-NN model is lower than the other models mentioned above. Nevertheless, the high values of R 2 of 0.986, 0.982, 0.977, and 0.963 for Ft, T, P, and Fd, respectively, are higher than that obtained using the traditional training technique (the Levenberg-Marquardt algorithm), which was examined too.
The response surface can be plotted by applying the obtained predictive models and used in optimizing the drilling process by determining the optimal machining parameters. The ability of developed models to relearn for enhancing their predictive capability with increasing data enables them to facilitate the exploitation of artificial intelligence in the FRP manufacturing.    Figure 6c shows the accuracy of the F d PSO-NN model is lower than the other models mentioned above. Nevertheless, the high values of R 2 of 0.986, 0.982, 0.977, and 0.963 for F t , T, P, and F d , respectively, are higher than that obtained using the traditional training technique (the Levenberg-Marquardt algorithm), which was examined too.
The response surface can be plotted by applying the obtained predictive models and used in optimizing the drilling process by determining the optimal machining parameters. The ability of developed models to relearn for enhancing their predictive capability with increasing data enables them to facilitate the exploitation of artificial intelligence in the FRP manufacturing. Figures 7-9 show the response surface plots for the four drilling process characteristics plotted using the models obtained, using both RSA regression and PSO-NN as functions in the f, s, a, and t variables. These models can be used to estimate response values for any appropriate set of drilling parameters. On other hand, Figure 10 shows the response surface plot of cutting power versus feed, spindle speed, and drill point angle plotting by Equation (6).
Processes 2023, 11, x FOR PEER REVIEW 12 of 18  show the response surface plots for the four drilling process characteristics plotted using the models obtained, using both RSA regression and PSO-NN as functions in the f, s, a, and t variables. These models can be used to estimate response values for any appropriate set of drilling parameters. On other hand, Figure 10 shows the response surface plot of cutting power versus feed, spindle speed, and drill point angle plotting by Equation (6).     show the response surface plots for the four drilling process characteristics plotted using the models obtained, using both RSA regression and PSO-NN as functions in the f, s, a, and t variables. These models can be used to estimate response values for any appropriate set of drilling parameters. On other hand, Figure 10 shows the response surface plot of cutting power versus feed, spindle speed, and drill point angle plotting by Equation (6).

Comparison of the Obtained Predictive Models
Performance statistical measures: R 2 , MSE, and mean absolute percentage error (MAPE) of predictive models are used to evaluate and compare the two modeling techniques used. MAPE is calculated using Equation (10).
where and O are the measured and predicted results of i run, and n is the total number of runs. These measures illustrated in Table 5 evidence the capability of obtained predictive models to describe the experimental results. The superiority of the PSO-NN models is evident.

Comparison of the Obtained Predictive Models
Performance statistical measures: R 2 , MSE, and mean absolute percentage error (MAPE) of predictive models are used to evaluate and compare the two modeling techniques used. MAPE is calculated using Equation (10).
where I and O are the measured and predicted results of i run, and n is the total number of runs. These measures illustrated in Table 5 evidence the capability of obtained predictive models to describe the experimental results. The superiority of the PSO-NN models is evident. Within the range of low error rates mentioned in Table 5, it can be said that the models successfully describe the drilling process and predict process responses with high accuracy, and can be used for optimization.

Optimizing Process Outputs
According to Table 6 and Figure 11, the optimal inputs of the drilling process of WGFRP (f = 0.025 mm/r, s = 655 rpm, a = 117.7 • , and t = 3.88 mm) resulted in an optimal thrust force, torque, power consumption, and delamination factor of 32.3 N, 10.66 N.cm, 3.88 W, and 1.387, respectively, with the desirability of 0.91. Paying more attention to the quality of the hole produced and the power consumed, we find that the second solution is optimal for the application, resulting in thrust force, torque, power consumption, and delamination factors of 35.2 N, 15 N.cm, 2.3 W, and 1.35, respectively, at combination (f = 0.025 mm/r, s = 401 rpm, a = 112 • , and t = 4.75 mm), with a desirability of 0.896. It is worth taking advantage of such models to decide which is the best thickness during the product design phase, as well as the best process plan for drilling holes in WGFRP composites.

Conclusions
The impact of drilling parameter feed rate and cutting speed using different drill point angles through different laminate thicknesses on the power consumption during cutting, as well as delamination, were investigated. A machine learning model was developed based on an ANN enhanced by the PSO algorithm to predict the critical thrust force, drilling torque, delamination factor, and cutting power consumption. Predictive models developed were used in multi-responses optimization of the drilling process. The results of this study were as follows: • Multivariate regression explains that the effect of feed rate on the power consumed is directly proportional.

•
Feed rate is prevalent in influencing the delamination factor among the studied control factors, followed by drilling point angle. Also, the speed had a statistically significant effect on the delamination, but with a slight contribution.

•
Obviously, a combination of a lower feed rate and a higher spindle speed reduces the

Conclusions
The impact of drilling parameter feed rate and cutting speed using different drill point angles through different laminate thicknesses on the power consumption during cutting, as well as delamination, were investigated. A machine learning model was developed based on an ANN enhanced by the PSO algorithm to predict the critical thrust force, drilling torque, delamination factor, and cutting power consumption. Predictive models developed were used in multi-responses optimization of the drilling process. The results of this study were as follows: • Multivariate regression explains that the effect of feed rate on the power consumed is directly proportional.
• Feed rate is prevalent in influencing the delamination factor among the studied control factors, followed by drilling point angle. Also, the speed had a statistically significant effect on the delamination, but with a slight contribution. • Obviously, a combination of a lower feed rate and a higher spindle speed reduces the delamination, as with the torque generated. On the contrary, however, the higher speed level increases the power consumption.

•
The effect of feed rate came out on top for thrust force (84.6%), torque (43.5%), and delamination factor (44.4%). Meanwhile, it was found that the rotational speed provided the largest contribution to the power consumption, followed by the feed rate, due to the dependence of the calculated power on the rotational speed and the torque generated during cutting. • It has been observed that the drill point angle has an effect on the critical thrust force, and hence this effect on the push-out delamination is amplified, especially with the highest feed rates and speeds.

•
In order to achieve a sustainable drilling process that takes into account the quality of the holes produced, represented by the delamination factor, and keeps the cost of the energy consumed low, controlling the feed rate and cutting speed simultaneously is necessary. • Models of RSA regression and PSO-NN developed were able to predict drilling process characteristics showing a very high consistency with the measured data. However, the PSO-NN models were more accurate than the others. • An optimal combination of factors (f = 0.025 mm/r, s = 401 rpm, a = 112 • , and t = 4.75 mm) was found to obtain the optimum produced hole quality as well as a low cutting power consumed, with an overall desirability factor of 91%.