1. Introduction
Fibre metal laminates, often known as FMLs, are a pioneering kind of advanced composite material that has spurred a revolution in the use of composite materials in aerospace and structural engineering applications. FMLs have emerged as a viable alternative to fulfil the ever-growing need for lightweight, high-strength materials [
1]. By combining the great mechanical capabilities of both metal alloys and high-performance fibers, FMLs have been able to combine the exceptional mechanical properties of both types of materials. The mechanical properties of FMLs are considerably improved with the addition of nano fillers into the polymer matrix materials. Recently, several researchers have studied the impact of incorporating nano fillers such as nano clay, carbon nano tubes (CNTs), titanium oxide (TiO
2), graphene oxides (GO and r-GO), and silicon oxide (SiO
2), and their concentrations, on the mechanical properties of FMLs. Among the various nano fillers, reduced graphene oxide (r-GO), a novel two-dimensional derived member of nanocarbon, has gained significant attention from researchers due to its outstanding properties, such as higher surface area, exceptional elastic modulus, and improved electrical conductivity [
2,
3,
4]. Because of the one-of-a-kind synergy that exists between metallic layers and fibre-reinforced composites, FMLs provide a harmonic balance of stiffness, strength, and damage tolerance. As a result, FMLs are very helpful in a variety of important sectors [
5]. The inherent complexity of FMLs results in a challenging set of problems that must be overcome before FMLs can be considered successful for end-use applications. One of the most significant obstacles is the heterogeneity of FMLs, which occurs when alternating layers of metal and fibre-reinforced composites display significantly different mechanical characteristics. This heterogeneity frequently leads to concerns such as differential wear on cutting tools and uneven material removal rates, both of which may impact the dimensional accuracy of machined components [
6]. Furthermore, the abrasive nature of the reinforcing fibers inside FMLs may result in quick tool wear and delamination, which poses substantial challenges to obtaining the necessary surface finish and limiting damage to the material. As FMLs continue to gain popularity in aerospace and structural applications, a more in-depth knowledge of these machining issues becomes important. This drives the need for creative machining processes and tactics to unleash the full potential of these sophisticated materials [
7].
Abrasive waterjet (AWJ) machining has emerged as a promising unconventional machining process for the precise shaping and trimming of FMLs. FMLs consist of alternating layers of metal and fibre-reinforced composites, requiring machining processes that can effectively minimize delamination and fibre damage while preserving the structural integrity of the metallic layers. AWJ machining, known for its non-thermal and non-contact characteristics, has demonstrated significant advantages, such as precision and versatility for intricate designs and shapes, which reduce the need for secondary finishing processes. Additionally, it is environmentally friendly, as it produces minimal waste and does not emit harmful fumes or gases, setting it apart as a sustainable choice for modern manufacturing when compared with other non-traditional machining processes. By utilizing a high-velocity flow of water mixed with abrasive particles, AWJ machining enables precise cutting and trimming of FMLs while minimizing the formation of heat-affected zones and reducing the risk of mechanical damage [
8]. The versatility of this technique allows for the machining of intricate geometries with exceptional accuracy, making it highly suitable for the complex structures and architectures commonly encountered in FMLs. On the other hand, achieving effective cut quality characteristics in AWJ cutting of novel engineering materials, especially in stacked composite laminates, is still difficult due to the abundance of process related variables and the dominance of quality deficiencies such as delamination, fibre pull-out, a higher kerf ratio, and inferior surface quality [
9]. Therefore, in order to achieve defect-free components that fulfil functional requirements, particularly for FMLs, it is necessary to conduct in-depth research into the effects of AWJ cutting parameters on these response properties.
In the recent past, several researchers have extensively investigated the impact of AWJ cutting on composite laminates. The AWJ cutting characteristics of a composite consisting of lanthanum phosphate and yttria were investigated by K. Balamurugan et al. [
10] using a response surface approach. The impacts of kerf taper, surface quality, and material removal rate, along with cutting speed, stand-off distance, and water pressure, have been taken into consideration as input factors. The results of their study indicate that the cutting speed has a detrimental effect on surface quality, while water pressure has a beneficial effect on both kerf taper and material removal rate. The effects of significant parameters of AWJ cutting on the kerf taper and surface roughness of a hybrid composite made of glass and coconut fibre were investigated by Kalirasu et al. [
11] using the Taguchi methodology. Based on the findings of the investigation, it can be concluded that the chosen response characteristics are primarily influenced by the size of the abrasive particles. Alberdi et al. [
12] conducted a drilling operation through AWJ machining on Ti/CFRP stacked composite laminates to investigate their qualitative traits. The research findings indicate that the arrangement of stacks, water pressure, and traverse speed of the focusing tube play a crucial role in determining the taper of drilling and surface quality. Hutyrová et al. [
13] conducted an experimental investigation to examine the surface topography of wood plastic composites machined using AWJ. Based on their findings, the utilization of AWJ cutting was determined to be an effective method for machining plastic composites while avoiding the melting of the matrix materials. Kalirasu et al. [
14] employed an analytical technique known as response surface methodology (RSM) and multi-objective optimization based on ratio analysis (MOORA) to evaluate the performance of jute/polyester composites in terms of AWJ cutting assessment. An investigation was conducted to examine the impact of stand-off distance, cutting speed, and jet pressure on the kerf taper and surface roughness of machined surfaces. The results indicate that the proposed method is capable of effectively handling fibre-reinforced composites with a maximum thickness of 3 mm. Pahuja et al. [
15] conducted experimental and statistical analyses to investigate the surface and kerf characteristics of stacked titanium (Ti) and carbon fibre reinforced polymer (CFRP) fibre metal laminates (FMLs) using AWJ cutting. It has been asserted that the uppermost layer of the skin is the site where micro buckling and fractures initiate the erosion mechanism of FMLs. Furthermore, the influence of the metal skin and polymer composite designs on the failures of the FML interface is significant. Ramulu [
16] conducted an analytical investigation of the influence of AWJ cutting factors on the kerf quality of machined Ti/CFRP composite stacks. An empirical model was developed to predict the depth of penetration and material removal mechanism resulting from water jet pressure. It was found that this model can effectively guide the adjustment of parameters in AWJ cutting to achieve defect-free cutting zones. However, only a few research studies are available on the multi-response optimization of AWJ cutting on FMLs, especially when adopting metaheuristic algorithms.
This study focuses on the analysis of the AWJ cutting capabilities of a newly developed fibre intermetallic laminate (FIL). The FIL is composed of reduced graphene oxide (r-GO)-filled epoxy prepregs, which are reinforced with a combination of carbon and aramid fibres. Additionally, Nitinol shape memory alloy sheets are incorporated as the skin materials in the FIL structure. The independent variables considered in this research include traverse speed, waterjet pressure, and nozzle height, along with different wt% values of r-GO nano fillers. The influence of these independent variables on three response characteristics, namely material removal rate, delamination factor, and kerf deviation, were investigated through statistical analysis and response surface plots. A metaheuristic-based moth–flame optimization (MFO) algorithm was used to optimize the AWJ cutting parameters to improve the cut quality features, with its results compared to those of other, more well-known metaheuristic algorithms, such as the genetic algorithm (GA), particle swarm algorithm (PSO), grey-wolf algorithm (GWO), and dragonfly algorithm (DFO).
4. Results and Discussion
In order to assess the effectiveness of the experimental approach, it was necessary to undertake statistical analysis on the conducted experiments and the resulting response variables. The study used a multi-parametric analysis of variance (ANOVA) to examine several statistical measures, including the coefficient of determination, sum of squares, lack of fit, individual, interaction, quadratic effects, and F-statistics. This analysis was conducted on chosen response features, namely MRR, FD, and KD, of the test samples processed using the AWJ cutting process. The investigation was carried out in three distinct phases. Initially, ANOVA was used to assess the statistical significance and the impact of selected independent parameters on the performance and quality characteristics of cut specimens. The subsequent phase included the formulation of second-order polynomial equations for each response feature, with the aim of establishing the correlation between the selected dependent variables and the features of the response. During the third step, the metaheuristic algorithms were used to identify the most optimal parametric combinations.
4.1. Statistical Analysis of Developed Polynomial Models
Multi-parametric ANOVA is used to conduct comprehensive investigations on the statistical implications of response characteristics acquired from experimental data.
Table 3,
Table 4 and
Table 5 provide the findings of the ANOVA for selected output responses, including
MRR,
FD, and
KD. With a 95% level of confidence, the AWJ cutting experiments were carried out, and Design Expert™ 13 software was used to perform the statistical analysis. The development of quadratic polynomial models for subsequent studies followed an analysis of statistical indicators such individual and interaction effects, coefficient of determination, lack of fit, and sum of squares. The stated AWJ cutting variables have a considerable influence on the response features, according to the statistical results from ANOVA. The lack-of-fit values were also found to be statistically significant within the chosen range of processing parameters, indicating the significance of the experiments that were carried out. In addition, the backward elimination technique was used to exclude the insignificant parametric terms with probabilities greater than 0.05 from the ANOVA in order to develop polynomial models. Multiple-regression coefficients (R
2) of 0.9955, 0.9954, and 0.9992 for
MRR,
FD, and
KD, respectively, indicate the closeness of the derived models to the actual data. Equations (11)–(13) represent the obtained polynomial regression models following the elimination of insignificant parameter combinations. Residual plots were used in addition to ANOVA to statistically confirm the results of the accomplished AWJ cutting experimental trials.
Figure 3a–c provide a summary of the statistical analysis conducted for each response characteristic. Based on the summary reports, it has been observed that the distribution of data points along the center line provides evidence of the statistical significance of the developed models, as demonstrated by the normal probability plots. Additionally, the distributional frequency of data, as depicted by the histogram, shows that the measured responses are well constructed at a 95% confidence level. These findings indicate that the measured experimental data can be effectively utilized for further investigations.
4.2. Influence of AWJ Cutting Parameters on Response Characteristics
Increased MRR with reduced FD and KD are desirable during the mechanical-based AWJ cutting process to enhance cutting quality. An enhanced MRR is often sought by manufacturing industries to improve the production rate while also reducing manufacturing costs, and FD and KD are quality indices of processed end-use components that should be decreased. The refining of the delamination factor and the reduction in kerf deviation have considerably improved the precision of the cutting process, reducing inaccuracy during assembly and decreasing substrate material depletion. From the statistical investigation results, the considered AWJ cutting parameters and the addition of r-GO nano fillers to the fabrication of FILs has a significance influence on the cut quality characteristics.
Figure 4a–f provide valuable insights into the intricate relationship between the AWJ cutting parameters and the cut quality characteristics, where two variables are varied while the other two variables are kept unchanged. From the statistical investigation, it is found that all the selected process parameters have a significant influence on
MRR (
Table 3). The influence of the addition of r-GO and the nozzle height on the
MRR is indicated in
Figure 4a. From the response surface plot, the maximized
MRR was attained at an augmented nozzle height (4 mm) and r-GO addition (2 wt%). An increase in nozzle height causes depth craters by increasing the momentum of hard abrasive particles that contact the substrate. In turn, this increases substrate surface erosion, allowing for the recovery of increased
MRR [
19]. The impact of traverse speed and the waterjet pressure on the
MRR is indicated in
Figure 4b. The surface plot indicates that the
MRR was increased linearly by increasing the traverse speed as well as the waterjet pressure. However, the maximized
MRR was attained at higher traverse speed (600 mm/min) or higher waterjet pressure (200 MPa). This may be attributed to the hard erosion of abrasive particles at higher jet pressure, along with maximized cutting speed, which leads to an improved removal of substrate [
14].
Figure 4c,d shows the influence of AWJ parameters on the delamination factor of the machined FILs at different cutting conditions. From the plots, it is perceived that the
FD was increased with an increase in all the selected parameters, whereas the quantified
FD was seen at lower waterjet pressure, nozzle height, and unreinforced FILs. This may be attributed to the complete fracture and debonding of plies in FILs due to higher energy levels associated with the pressurized waterjet, which leads to an augmented delamination [
20]. Furthermore, it is noted that the waterjet bends at greater traverse speeds, which produces a curved cutting front that results in the highest normal force. As a result, the metal-composite interface sustains significant damage, which increases the factor of delamination [
21].
The effects of AWJ cutting parameters on the kerf deviation are shown in response surface plots (
Figure 4e,f). The kerf deviation is found linear increase with an increase in traverse speed and addition of r-GO fillers up to a certain level, and then decreasing, whereas it is found to linearly increase with an increase in nozzle height and waterjet pressure. A faster traverse speed during the AWJ cutting process permits abrasive particles to travel swiftly across the substrate, resulting in an incorrect cutting kerf and decreased quality characteristics. Therefore, keeping a consistent traverse speed is critical to achieve a lower kerf variation [
22]. Moreover, the peak velocity of the waterjet expands (the jet diverges) when it leaves the mixing chamber at a greater nozzle height. When the jet swerves, kinetic energy is lost during penetration of the substrate surface, and it may not have enough force on its visible surface to cut effectively. This effect causes kerf deviation to rise as nozzle height increases [
23].
4.3. Optimization of AWJ Cutting Process
The objective of this study is to address an optimization problem that involves the simultaneous optimization of multiple objectives; specifically, to maximize the
MRR while minimizing the
FD and
KD. To achieve this, the decision variables associated with the AWJ cutting process, including traverse speed, waterjet pressure, and nozzle height, along with the impact of incorporating reduced graphene oxide (r-GO) in the composite laminates, were considered for optimization. In general, a process parameter set for one objective function cannot be used for any other objective functions. Since the objective functions chosen for this study conflict with each other, the optimization process is difficult. Usually, there are two ways to solve such complicated optimization problems. The first is to turn multiple objectives into a single objective by giving each objective a weight or a utility function. The second is to find non-dominated Pareto optimum settings for the decision variables [
24]. This study utilizes non-dominated solutions for each response characteristic to obtain the optimal processing parameters.
The objective functions for the present optimization work are formulated using the second-order polynomial equations obtained through the statistical analysis (Equations (11)–(13)) and the boundary conditions of the selected processing parameters (
Table 1). For the multi-objective optimization, a highly versatile metaheuristic-based moth–flame optimization (MFO) algorithm was considered to simultaneously improve the
MRR and minimize the
FD and
KD of AWJ-processed FILs. In order to assess the performance characteristics of the proposed MFO algorithm, several established metaheuristic algorithms were taken into consideration. These algorithms include the genetic algorithm (GA), particle swarm optimization (PSO) algorithm, dragonfly algorithm (DFO), and grey-wolf optimization algorithm (GWO). By comparing various performance features, such as convergence plot, diversity, spacing, and hypervolume values, a comprehensive evaluation of the MFO algorithm can be achieved.
The proposed optimization algorithms were executed in a MATLAB 2022b™ environment with a population size and number of iterations of 30 and 100, respectively, for each algorithm. The parameters used for executing different optimization algorithms to obtain non-dominated Pareto optimal solutions are shown in
Table 6. Each optimization algorithm was executed thirty-six times to obtain the number of Pareto optimal solutions. For each execution, an optimal combination of process parameters and their corresponding response values were determined and considered as the best solution for that execution. Similarly, each algorithm yielded thirty-six sets of optimal parameters for each response characteristic, which were considered for further investigations.
Figure 5a–c depict a sample convergence plot derived for each response characteristic using various algorithms. Among the obtained thirty-six non-dominated Pareto optimal solutions, Deng’s similarity-based analytical hierarchy approach [
25] was adopted for each algorithm to obtain the optimal solution based on the overall performance index. The optimal AWJ cutting parameters and their corresponding response characteristics, such as
MRR,
FD, and
KD, achieved through Deng’s approach, are listed in
Table 7.
4.4. Comparison of MFOA Performance with State-of-the-Art Metaheuristic Algorithms
In this section, the effectiveness of MFO is compared to that of other metaheuristic strategies such as GA, PSO, GWO, and DFO. The results of this comparison analysis show how MFO is more efficient than other algorithms in AWJ cutting applications. From the convergence plots shown in
Figure 5a–c, MFO is found to outperform the similar algorithms in terms of the convergence with minimal iterations. MFO converges quickly (i.e., the number of iterations for MRR is 5, that for FD is 12, and that for KD is 3) to obtain the optimal solutions; hence, the execution time of MFO can be significantly minimized for such complex optimization problems relative to the compared metaheuristics.
4.4.1. Diversity Values
Diversity values are used to quantify the difference or similarity of the solutions produced by each algorithm throughout the optimization or decision-making process. Diversity values not only help to determine which strategy is best for a particular optimization problem, but also provide insight into the flexibility and resilience of certain algorithms in situations when complicated decisions must be made [
26]. The diversity values between the two successive optimal solutions were calculated based on the following mathematical relation:
The diversity values for selected algorithms are depicted in
Table 8. The diversity values for all the algorithms exhibited a high degree of similarity. The statistical analysis, including the Anderson–Darling normality test, normal probability analysis, and analysis of variance, was conducted on the diversity values obtained in order to assess the performance of the algorithms. Based on the data presented in
Table 9, it can be observed that the statistical indicators of MFO algorithms demonstrate superior performance compared to other optimization algorithms. Traditionally, diversity values have been evaluated using a higher-the-better approach. The
p-value for the MFO algorithm is determined to be lower (0.006) in comparison to that of other metaheuristic algorithms. Therefore, the MFO algorithm demonstrates superior performance in optimizing the AWJ cutting parameters. The significance of the proposed algorithms is confirmed by the normal probability plots depicted in
Figure 6a–e at a 95% confidence interval. Furthermore, the statistical significance of the proposed algorithms was confirmed by Friedman’s ANOVA, as depicted in
Figure 7, with a probability of less than 0.05.
4.4.2. Spacing Values
The evaluation of the algorithms’ performance via the adjustment of spacing values is an essential process in enhancing their efficiency and efficacy. This research assessed five algorithms, and the findings provided insights into the individual strengths and limitations of each algorithm. The spacing values, which denote the distance between data points or items inside the algorithmic process, were significant in influencing the resulting performances. The spacing values can be calculated by the following relation [
27]:
where
denotes the minimum value of the sum of the absolute difference between the
ith run and all other runs,
indicates the mean value of
, and
indicates the number of runs.
and
can be calculated as follows:
Table 10 illustrates the spacing values for selected algorithms. It is observed that the spacing values for all the algorithms were close to each other. Hence, the statistical analysis, such as the Anderson–Darling normality test, normal probability analysis, and analysis of variance, was performed for the obtained spacing values to identify the performance of the algorithms. From
Table 11, it is perceived that the statistical indicators of the MFO algorithms are found to be better than those of the other optimization algorithms. In general, lower spacing values are considered to be better. The
p-value of MFO is found to be less (0.077) that that of the compared similar metaheuristics. Hence, the MFO algorithms outperform when optimizing the AWJ cutting parameters. Moreover, the normal probability plots, as shown in
Figure 8a–e, at a 95% confidence interval, confirm the significance of the proposed algorithms. In addition to these indicators, Friedman’s ANOVA, as shown in
Figure 9, proved the statistical significance of the proposed algorithms with a probability of less than 0.05.
4.4.3. Hypervolume Values
The hypervolume metric has significant importance within the field of optimization techniques, especially in the context of multi-objective optimization. It measures how well an algorithm can explore and cover the Pareto front, which is a collection of non-dominated solutions that reflect the best trade-offs between many competing goals. It is a key indicator of the quality of a solution. One may evaluate the coverage, variety, and effectiveness of an algorithm’s solutions by computing the hypervolume. A more thorough examination of the goal space is indicated by higher hypervolume values, which also indicate better solution quality performance [
28]. The mathematical expression for calculating the hypervolume can be defined as follows:
where
PF is defined as the Pareto front, the point of approximation of
PF can be expressed as
S, the generalization of a volume is defined by
, and the domination relation may be defined by
.
Table 12 presents the hypervolume values corresponding to the chosen algorithms. The hypervolume values for all the algorithms exhibited a high degree of similarity. The statistical analysis included the Anderson–Darling normality test, normal probability analysis, and analysis of variance. These tests were conducted on the obtained hypervolume values in order to assess the performance of the algorithms. Based on the data presented in
Table 13, it can be observed that the statistical indicators of GA and DFO algorithms demonstrate superior performance compared to other optimization algorithms. Traditionally, the higher-the-better approach has been favored when evaluating hypervolume values. The
p-values for GA and DFO were determined to be lower (0.022 and 0.024, respectively) when compared to those of other metaheuristics of a similar nature. Therefore, when evaluating the performance metric as the hypervolume, the GA and DFO algorithms were surpassed. The significance of the proposed algorithms was confirmed by the normal probability plots, as depicted in
Figure 10a–e, at a 95% confidence interval. Furthermore, the statistical significance of the proposed algorithms was confirmed by Friedman’s ANOVA, as depicted in
Figure 11, with a probability of less than 0.05.
The mean rank values of the optimization algorithms were assessed based on the performance indices, such as diversity, spacing, and hypervolume values, by conducting Friedman’s statistical analysis [
29]. To summarize the performance of distinct metaheuristics, the mean rank values were determined. The table displays the mean rank values for selected algorithms, and the MFO algorithm was found to provide higher performance index values, such as 4.9444 for diversity and 2.5833 for spacing, whereas GWO exhibits an improved hypervolume value of 3.5. When comparing the overall performance index, it can be inferred that the MFO algorithm exhibits superior performance in comparison to GA, PSO, GWO, and DFO while optimizing the AWJ cutting process based on the selected matrices, except for the hypervolume. Among the optimized parameters and their corresponding response characteristics obtained through different optimization algorithms, as shown in
Table 14, the MFO algorithm’s optimal values were considered as the best set of parameters to improve the cut quality and performance features of the AWJ cutting process. The corresponding optimal parameters were unreinforced laminate (i.e., 0 wt% r-GO), traverse speed of 600 mm/min, jet pressure of 253.36 MPa, and nozzle height of 2 mm, for improved response features, namely the corresponding
MRR,
FD, and
KD values are 3.096 g/s, 1.928, and 1.833 mm, respectively. The optimal parameters show that the inclusion of nano fillers does not provide any significance regarding the selected cutting responses, whereas the lower nozzle height with the combination of higher traverse speed and medium jet pressure may enhance the selected cutting characteristics within the range of selected parameters.
4.5. Confirmation Experiments
In order to evaluate the rationality of the proposed optimization approach, a series of confirmation experiments were conducted. These experiments aimed to validate the effectiveness and efficiency of the proposed approach. The confirmation experiments were conducted thrice using the optimal AWJ cutting parameters obtained through the MFO algorithm. The average values of these experiments are shown in
Table 15. The presented table demonstrates a relatively strong correlation between the predicted and experimentally measured response values. The average error for
MRR is found to be 1.94%, while the average errors for
FD and
KD are 2.75% and 3.33%, respectively. These results indicate a satisfactory level of accuracy in the prediction of the response values.
5. Conclusions
In the present study, novel fiber intermetallic laminates were cut using AWJ by modifying the wt% of reduced graphene oxide in the laminates, traverse speed, jet pressure, and nozzle height. The influence of AWJ parameters on the cut quality characteristics was investigated, and multi-response optimization using five distinct metaheuristic algorithms, i.e., GA, PSO, MFO, GWO, and DFO, was accomplished. The performance features of the algorithms were also investigated. The findings of this study are summarized as follows.
The statistical analysis reveals that JP and r-GO addition have a significant influence on MRR, followed by TS, whereas r-GO addition and NH have a substantial impact on FD. In addition, it was observed that KD is significantly influenced by WP and NH.
MRR is found to increase with an increase in NH and JP due to the increased momentum of hard abrasive particles that contact the substrate, which in turn increases substrate surface erosion, whereas it leads to augmented FD and KD, which deteriorate the cut quality. Therefore, it is necessary to maintain a consistent traverse speed, nozzle height, and jet pressure to achieve improved cut quality features.
From the convergence plots, the MFO algorithm is found to converge more rapidly than similar metaheuristic algorithms due to its simplicity and the presence of minimal control parameters. Therefore, the space and time complexity can be significantly reduced by the MFO algorithm for optimizing such complex engineering problems.
The optimal AWJ cutting parameters for improved quality features of FILs are unreinforced composite laminates (0 wt% of r-GO), 600 mm/min of TS, 253.36 MPa of JP, and 2 mm of NH. These parameter values were achieved using the moth–flame optimization algorithm.
The suggested MFO algorithm also showed its efficacy in improving the quality and performance features of AWJ cutting by predicting appropriate parameter settings because the error variation between predicted and experimental measures was determined to be less than 3.5% for all the response characteristics of the cut specimens.
The effectiveness of MFO was evaluated against GA, PSO, GWO, and DFO using the most widely used key metric indicators, including diversity, spacing, and hypervolume values. The comparative analysis results showed that the MFO algorithm produced lower spacing values and higher diversity and hypervolume values, which shows the efficiency of the algorithm.
The present work encompassed multi-response optimization of AWJ cutting parameters via different metaheuristic optimization algorithms with their standard control parameters. Future works may concentrate on tuning the algorithm parameters to improve their efficiency for better exploration and exploitation of similar optimization problems. Moreover, the effectiveness of the proposed algorithms may be improved by considering a range of workpiece materials, AWJ cutting parameters, and performance characteristics.