Next Article in Journal
Green Hydrogen Generation from Eco-Friendly and Cost-Effective Red Sea Water Using a Highly Photocatalytic Nanocomposite Film, As2O3/Poly-3-methylaniline
Previous Article in Journal
Preparation and Properties of Flexible CuI/Polyvinylpyrrolidone Nanocomposite Thermoelectric Film
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Response Optimization of Abrasive Waterjet Cutting on r-GO-Reinforced Fibre Intermetallic Laminates through Moth–Flame Optimization Algorithm

by
Devaraj Rajamani
1,*,
Mahalingam Siva Kumar
1 and
Esakki Balasubramanian
2
1
Centre for Advanced Materials Processing, Department of Mechanical Engineering, Vel Tech Rangarajan Dr. Sagunthala R&D Institute of Science and Technology, Chennai 600062, India
2
Department of Mechanical Engineering, National Institute of Technical Teachers Training and Research, Chennai 600113, India
*
Author to whom correspondence should be addressed.
J. Compos. Sci. 2023, 7(11), 462; https://doi.org/10.3390/jcs7110462
Submission received: 10 October 2023 / Revised: 24 October 2023 / Accepted: 31 October 2023 / Published: 3 November 2023
(This article belongs to the Section Composites Manufacturing and Processing)

Abstract

:
Laminated metal-composite structures, also known as fibre metal laminates (FMLs), have emerged as prominent engineering materials in various industries, particularly in the domains of aircraft and automobile manufacturing. These materials are sought after due to their enhanced impact and fatigue resistance capabilities. The machining of FMLs plays a crucial role in achieving near-net shapes for the purpose of joining and assembling components. Delamination is a prevalent issue encountered during the process of conventional machining, thus rendering FMLs are challenging materials to machine. This study aims to investigate the cutting process of novel fibre intermetallic laminates (FILs) using the abrasive water jet (AWJ) cutting technique. The FILs consists of carbon and aramid fibers that are adhesively bonded with a resin matrix filled with reduced graphene oxide (r-GO) nano fillers. Moreover, these laminates contain embedded Nitinol shape memory alloy sheets as the skin materials. Specifically, the study aims to investigate the impact of different factors, such as the addition of reduced graphene oxide (r-GO) in the laminates (ranging from 0 to 2 wt%), traverse speed (ranging from 400 to 600 mm/min), waterjet pressure (ranging from 200 to 300 MPa), and nozzle height (ranging from 2 to 4 mm), on the material removal rate (MRR), delamination factor (FD), and kerf deviation (KD). ANOVA was used in the statistical analysis to determine the most influential parameters and their effects on the selected responses. The optimal AWJC parameters are determined using a metaheuristic-based moth–flame optimization (MFO) algorithm in order to enhance cut quality. The efficacy of MFO is subsequently compared with similar well-established metaheuristics such as the genetic algorithm, particle swarm algorithm, dragonfly algorithm, and grey-wolf algorithm. MFO was found to outperform in terms of several performance indices, including rapid divergence, diversity, spacing, and hypervolume values, among the algorithms compared.

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 (TiO2), graphene oxides (GO and r-GO), and silicon oxide (SiO2), 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).

2. Methodologies

2.1. Response Surface Methodology

This study investigates the effects of selected AWJ cutting variables and the incorporation of r-GO in FILs on selected quality and performance characteristics including MRR, FD, and KD. Since there are numerous process-related parameters that are dynamically engaged in the material removal mechanism, the AWJ cutting process has become exceedingly complicated. Thus, the experimentation and prediction of optimal process variables requires a viable and systematic experimental plan in order to minimize the number of experimental trials, and thereby reduce processing costs and time consumption. Response surface methodology (RSM) has been found to be an effective technique for designing experiments with a minimum number of experimental trials in order to determine the interaction between dependent and predictor variables of complex processes such as machining, additive manufacturing, and thermo-mechanical systems [17]. In order to construct higher-order mathematical models for optimization, the current study employs Box–Behnken design (BBD), a subset of RSM, to design and conduct AWJ cutting experiments. The formulation of quadratic polynomial models by means of RSM is as follows:
R = β 0 i = 1 k β i x i + i = 1 k β i i x i 2 + i j β i j x i x j + ε
where R is the response and x i is the value of the ith AWJ cutting parameter; β0 is the model constant; β i , β i j , β i i represent the linear, interaction, and quadratic coefficients, respectively; and ε indicates the statistical distribution error. The relevant data are typically acquired through the design of experiments, followed by the statistical multiple regression approach, to develop the empirical models. To explain the relevance of the empirical model, a more popular statistical method such as analysis of variance (ANOVA) is used.

2.2. Moth–Flame Optimization Algorithm

The moth–flame optimization (MFO) algorithm is a novel stochastic population-based approach which is primarily inspired by moths’ natural transverse orientation navigational strategy [18]. Moths can only migrate in a straight line using this strategy when the light source is far away. The moth’s objective is to get to the flame. The flame’s objective is to achieve its ideal location and improve flame density. According to the flowchart shown in Figure 1, MFO begins by creating moths at random in the solution space, calculates each moth’s fitness values (i.e., location), and tags the ideal position with a flame. The next step is to update the moths’ positions using a spiral movement function to obtain better positions that are tagged by a flame, update the new best individual positions, and then repeat the previous steps (updating the moths’ positions and generating new positions) until the termination criteria are satisfied.
The overall structure of the MFO algorithm consists of three stages, namely initialization, iteration, and the stopping criterion. The step-by-step implementation of MFO algorithms is as follows:

2.2.1. Stage 1: Initializing the Population of Moths

The location of the moth in space must be solved as a variable in the MFO algorithm, which assumes that the moth represents a potential solution to the considered problem. Moths can fly in one, two, three, and hyperdimensional space by modifying their position vectors. Given that MFO is basically a swarm intelligence optimization technique, the moth population can be represented as follows:
M = m 11 m 12 m 21 m 22 m n 1 m n 2 m 1 d m 2 d m n d
where the number of moths can be denoted as n and the number of process parameters (dimensions of the problem) can be represented as d . The corresponding fitness values for the selected number of moths can be represented and memorized in an array as follows:
O M = O m 1 O m 2 O m n
where the fitness function value of the corresponding moth is denoted as O m i , which is determined according to the actual situation.
The remaining elements in MFO are flames, which can be denoted in D-dimensional space as a matrix [ F ] followed by its fitness function vector [ O F ].
F = f 11 f 12 f 21 f 22 f n 1 f n 2 f 1 d f 2 d f n d
O F = O f 1 O f 2 O f n
It must be noted that the MFO algorithm provides solutions for both moths and flames. The way they are handled and updated in each iteration is what makes them different from one another. The moths are the real searchers; they travel around the search area, and their current best position is among the flames. In other words, flames may be thought of as flags or pins that moths drop while scouring the search area. In order to update a flag (flame) if a better solution is discovered, each moth examines the area surrounding it. A moth never loses its ideal solution with this process.

2.2.2. Stage 2: Updating the Position of Moths

To obtain the global optimal solutions of the optimization problems, MFO uses three distinct functions. The following is a definition of these functions:
M F O = I , P , T
where I indicates the first random locations of the moths, P indicates the movement of the moths in the search area, T indicates the end of the search process. The populations of MFO are initialized as follows:
M i , j = u i l j × r a n d ( ) + l i
where the upper and lower bounds of the algorithm parameters are defined as u & l , respectively. The moths fly in a transverse direction in the search area, as was previously indicated. When using a logarithmic spiral, the following three requirements must be met:
  • The initial point of the spiral should start from the moth;
  • The end point of the spiral should be at the position of the flame;
  • The movement of the spiral should not exceed the search space.
To fulfil the above conditions, the logarithmic spiral for the MFO algorithm can be defines as follows:
S M i , F j = D i . e b t . cos 2 π t + F j
where the distance between the moth and flame F j M i is defined by D i , the shape of the logarithmic spiral is denoted by b , and a random number between −1 to 1 is defined by t . The spiraling motion of the moth towards the flame in the search space ensures that exploitation and exploration are balanced in MFO. Furthermore, utilizing the OF and OM matrices, the moths fly around the flames (i.e., each moth flies encircling the closest flame) in order to avoid the probability of obtaining a local optimum.

2.2.3. Stage 3: Updating the Number of flames

The number of flames should be less during the iteration for obtaining the global optimal solutions. Hence, an adaptive mechanism has been selected for enhancing the exploitation of the MFO algorithm by simultaneously reducing the number of flames. The mathematical expression for updating the number of flames can be represented as follows:
f N o = r o u n d N l × N l T
where the maximum number of flames can be represented as N , and the present iteration number and the maximum number of iterations are denoted as l   and   T , respectively.

3. Experimentation and Measurements

3.1. Fabrication of Fibre Intermetallic Laminates

The AWJ cutting experiments were conducted using FILs. These laminates were composed of Nitinol foils as the skin materials, epoxy resin as the matrix, and carbon/Kevlar as the prepreg materials. To fabricate these laminates, a vacuum-assisted resin infusion process was employed. Three different weight percentages of r-GO nano fillers (0%, 1%, and 2%) were incorporated into the laminates. The r-GO fillers with varying weight percentages were evenly distributed into the acetone-diluted epoxy resin using mechanical stirring. The stirring process lasted for approximately 15 min at a probe speed of 120 rpm. To achieve a homogeneous mixture, a sonication process was conducted using an ultrasonic bath (Model: Sonoplus HD 2000, manufactured by Bandelin Electronic GmbH & Co. KG, Germany). The process lasted for a duration of 45 min, with a power output of 70 W. The temperature of resin and r-GO mixtures was reduced during the sonication process by placing the vessel in an ice water bath. Following the sonication process, the acetone present in the resin mixture was removed by subjecting the mixture to heat in a vacuum oven at a temperature of 70 °C for a duration of 60 min. The hardener was combined with the resin in a ratio of 10:1 using a magnetic stirrer set at a stirring speed of 180 rpm for a duration of 15 min. This process ensures improved dispersion of r-GO particles and the removal of any trapped air voids within the mixture. The r-GO-reinforced FILs were produced using a vacuum-assisted resin infusion process. This involves stacking woven carbon and aramid fibers in an alternating pattern, with two layers of outer-skinned Nitinol foils. The stacking sequence is as follows: NiTi/[carbon/epoxy/aramid]8/NiTi. The fabricated laminates underwent a curing process at an ambient temperature for a duration of 24 h, while being subjected to a pressure of 6 bar. The fabricated FIL had cross-sectional dimensions of 300 × 300 × 3.5 ± 0.3 mm3.

3.2. Abrasive Waterjet Cutting of FILs

Using a high-precision cantilever type computer-numerical-controlled AWJ cutting machine (Model: S3015, Waterjet Germany, India), the cutting experiments were carried out on three different kinds of FILs made with 0, 1, and 2 wt% of r-GO nano fillers. To facilitate effective machining operations, the AWJ machine was outfitted with a high-pressure intensifier pump with an operating pressure of up to 450 MPa and a 0.76 mm diameter carbide nozzle. Throughout the studies, the machining parameters, including the abrasive flow rate of 200 g/min, the abrasive particle size of 177 µm, and the jet impact angle of 90°, were held constant. For the investigation of material removal rate, delamination factor, and kerf deviation, a total of 29 straight parallel slots with a cutting length of 50 mm were made in the various FILs in accordance with the RSM-BBD experimental design strategy for four dependent variables with three levels (Table 1). Figure 2 illustrates the AWJC experimental setup and the processed FIL samples. The AWJ cutting system has many process-related dependent variables, but the most important ones, such TS, WP, and NH, were chosen for the present study, and were obtained from exhaustive pilot studies based on the capability of the system and the available existing studies. Following the execution of the AWJ cutting experiments, the response variables MRR, FD, and KD were measured three times for each cut, and their mean values were considered for further studies to eliminate the statistical prediction errors. MRR represents the mass loss of substrate material during each experimental run and is calculated as follows:
M R R g / s = M a t e r i a l   r e m o v a l   f o r   e a c h   c u t × c u t t i n g   s p e e d L e n g t h   o f   c u t
A video measuring system (made by Accurate Gauge, India) outfitted with a high-resolution scanning system and M2-2D software was used to measure the delamination factor (FD), which is defined as the peel-up and push-out delamination along the slot of through cuts. The ratio between the maximum width of the delaminated area and the actual width of the cut along its length is used to calculate the delamination factor. The top kerf width of the cut surface was measured using a table top co-ordinate measuring machine. The measured response characteristics for the proposed experimental design are presented in Table 2.

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 (R2) 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.
M R R g / s = 6.433 0.028 × A 0.0014 × B 0.0302 × C + 0.721 × D 0.0006 × A B + 0.0005 × A C + 0.0712 × A D 0.00002 × B C 0.0008 × B D + 0.0178 × A 2 + 7.97 × 10 6 × B 2 + 0.00007 × C 2 0.067 × D 2
F D = 2.877 0.434 × A 0.004 × B + 0.006 × C 0.161 × D + 0.001 × A B + 0.068 × A D 0.000009 × B C + 0.0005 × B D + 0.0007 × C D 0.125 × A 2 + 3.529 × 10 6 × B 2 5.583 × 10 6 × C 2 0.0427 × D 2
K D   ( m m ) = 3.345 + 1.285 × A + 0.015 × B + 0.00041 × C + 0.8086 × D 0.0001 × A B 0.002 × A C 0.149 × A D + 1.35 × 10 6 × B C 0.0002 × B D 0.001 × C D 0.096 × A 2 1.548 × 10 5 × B 2 + 7.91 × 10 6 × C 2 0.044 × D 2

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:
D = j = 1 k f j max f j min 2
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]:
S p a c i n g = 1 n r 1 i = 1 n r S ¯ S i 2
where S i denotes the minimum value of the sum of the absolute difference between the ith run and all other runs, S ¯ indicates the mean value of S i , and n r indicates the number of runs. S i and S ¯ can be calculated as follows:
S i = min i = 1 , 2 n r j = 1 n o a b s O i j O k j k = i + 1 , n r
S ¯ = i = 1 n r S i n r
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:
H V P F = Λ s P F S : S S S n a d i r
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.

Author Contributions

Conceptualization, D.R. and M.S.K.; Methodology, D.R., M.S.K. and E.B.; Software, D.R. and M.S.K.; Validation, D.R., M.S.K. and E.B.; Formal analysis, M.S.K. and D.R.; Investigation, D.R., E.B. and M.S.K.; Data curation, D.R. and E.B.; Writing—original draft preparation, D.R. and M.S.K.; Writing—review and editing, M.S.K. and E.B.; Supervision, M.S.K.; Project administration, D.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research does not receive any funding from external sources.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

They are available in the same article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Costa, R.D.F.S.; Sales-Contini, R.C.M.; Silva, F.J.G.; Sebbe, N.; Jesus, A.M.P. A Critical Review on Fiber Metal Laminates (FML): From Manufacturing to Sustainable Processing. Metals 2023, 13, 638. [Google Scholar] [CrossRef]
  2. Jamali, N.; Khosravi, H.; Rezvani, A.; Tohidlou, E. Mechanical properties of multiscale graphene oxide/basalt fiber/epoxy composites. Fiber. Polym. 2019, 20, 138–146. [Google Scholar] [CrossRef]
  3. Saeimi Sadigh, M.A.; Marami, G. Investigating the effects of reduced graphene oxide additive on the tensile strength of adhesively bonded joints at different extension rates. Mater. Des. 2016, 92, 36–43. [Google Scholar] [CrossRef]
  4. Emran, M.Y.; Miran, W.; Gomaa, H.; Ibrahim, I.; Belessiotis, G.V.; Abdelwahab, A.A.; Othman, M.B. Biowaste Materials for Advanced Biodegradable Packaging Technology. In Handbook of Biodegradable Materials; Springer International Publishing: Cham, Switzerland, 2022; pp. 1–37. [Google Scholar]
  5. Smolnicki, M.; Lesiuk, G.; Duda, S.; de Jesus, A.M.P. A Review on Finite-Element Simulation of Fibre Metal Laminates. Arch. Comput. Methods Eng. 2023, 30, 749–763. [Google Scholar] [CrossRef]
  6. Marques, F.; Silva, F.G.A.; Silva, T.E.F.; Rosa, P.A.R.; Marques, A.T.; de Jesus, A.M.P. Delamination of Fibre Metal Laminates Due to Drilling: Experimental Study and Fracture Mechanics-Based Modelling. Metals 2022, 12, 1262. [Google Scholar] [CrossRef]
  7. Pai, A.; Kini, C.R.; Shenoy, B.S. Scope of Non-Conventional Machining Techniques for Fibre Metal Laminates: A Review. Mater. Today 2022, 52, 787–795. [Google Scholar] [CrossRef]
  8. Siva Kumar, M.; Rajamani, D.; El-Sherbeeny, A.M.; Balasubramanian, E.; Karthik, K.; Hussein, H.M.A.; Astarita, A. Intelligent Modeling and Multi-Response Optimization of AWJC on Fiber Intermetallic Laminates through a Hybrid ANFIS-Salp Swarm Algorithm. Materials 2022, 15, 7216. [Google Scholar] [CrossRef]
  9. Doğankaya, E.; Kahya, M.; Özgür Ünver, H. Abrasive Water Jet Machining of UHMWPE and Trade-off Optimization. Mater. Manuf. Process. 2020, 35, 1339–1351. [Google Scholar] [CrossRef]
  10. Balamurugan, K.; Uthayakumar, M.; Sankar, S.; Hareesh, U.S.; Warrier, K.G.K. Predicting Correlations in Abrasive Waterjet Cutting Parameters of Lanthanum Phosphate/Yttria Composite by Response Surface Methodology. Measurement 2019, 131, 309–318. [Google Scholar] [CrossRef]
  11. Kalirasu, S.; Rajini, N.; Winowlin Jappes, J.T.; Uthayakumar, M.; Rajesh, S. Mechanical and Machining Performance of Glass and Coconut Sheath Fibre Polyester Composites Using AWJM. J. Reinf. Plast. Compos. 2015, 34, 564–580. [Google Scholar] [CrossRef]
  12. Alberdi, A.; Artaza, T.; Suárez, A.; Rivero, A.; Girot, F. An Experimental Study on Abrasive Waterjet Cutting of CFRP/Ti6Al4V Stacks for Drilling Operations. Int. J. Adv. Manuf. Technol. 2016, 86, 691–704. [Google Scholar] [CrossRef]
  13. Hutyrová, Z.; Ščučka, J.; Hloch, S.; Hlaváček, P.; Zeleňák, M. Turning of Wood Plastic Composites by Water Jet and Abrasive Water Jet. Int. J. Adv. Manuf. Technol. 2015, 84, 1615–1623. [Google Scholar] [CrossRef]
  14. Kalirasu, S.; Rajini, N.; Rajesh, S.; Jappes, J.T.W.; Karuppasamy, K. AWJM Performance of Jute/Polyester Composite Using MOORA and Analytical Models. Mater. Manuf. Process. 2017, 32, 1730–1739. [Google Scholar] [CrossRef]
  15. Pahuja, R.; Ramulu, M.; Hashish, M. Surface Quality and Kerf Width Prediction in Abrasive Water Jet Machining of Metal-Composite Stacks. Compos. B Eng. 2019, 175, 107134. [Google Scholar] [CrossRef]
  16. Pahuja, R. Ramulu Abrasive Water Jet Machining of Titanium (Ti6Al4V)–CFRP Stacks–A Semi-Analytical Modeling Approach in the Prediction of Kerf Geometry. J. Manuf. Process. 2019, 39, 327–337. [Google Scholar] [CrossRef]
  17. Hazir, E.; Ozcan, T. Response Surface Methodology Integrated with Desirability Function and Genetic Algorithm Approach for the Optimization of CNC Machining Parameters. Arab. J. Sci. Eng. 2019, 44, 2795–2809. [Google Scholar] [CrossRef]
  18. Mirjalili, S. Moth-Flame Optimization Algorithm: A Novel Nature-Inspired Heuristic Paradigm. Knowl. Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  19. Palleda, M. A Study of Taper Angles and Material Removal Rates of Drilled Holes in the Abrasive Water Jet Machining Process. J. Mater. Process. Technol. 2007, 189, 292–295. [Google Scholar] [CrossRef]
  20. Mm, I.W.; Azmi, A.I.; Lee, C.C.; Mansor, A.F. Kerf Taper and Delamination Damage Minimization of FRP Hybrid Composites under Abrasive Water-Jet Machining. Int. J. Adv. Manuf. Technol. 2018, 94, 1727–1744. [Google Scholar] [CrossRef]
  21. Kumar, D.; Gururaja, S. Abrasive Waterjet Machining of Ti/CFRP/Ti Laminate and Multi-Objective Optimization of the Process Parameters Using Response Surface Methodology. J. Compos. Mater. 2020, 54, 1741–1759. [Google Scholar] [CrossRef]
  22. Rajamani, D.; Balasubramanian, E.; Murugan, A.; Tamilarasan, A. AWJC of NiTi Interleaved R-GO Embedded Carbon/Aramid Fibre Intermetallic Laminates: Experimental Investigations and Optimization through BMOA. Mater. Manuf. Process. 2023, 38, 1144–1158. [Google Scholar] [CrossRef]
  23. Li, M.; Huang, M.; Yang, X.; Li, S.; Wei, K. Experimental Study on Hole Quality and Its Impact on Tensile Behavior Following Pure and Abrasive Waterjet Cutting of Plain Woven CFRP Laminates. Int. J. Adv. Manuf. Technol. 2018, 99, 2481–2490. [Google Scholar] [CrossRef]
  24. Ciurana, J.; Arias, G.; Ozel, T. Neural Network Modeling and Particle Swarm Optimization (PSO) of Process Parameters in Pulsed Laser Micromachining of Hardened AISI H13 Steel. Mater. Manuf. Process. 2009, 24, 358–368. [Google Scholar] [CrossRef]
  25. Ibrahim, M.A.; Çamur, H.; Savaş, M.A.; Abba, S.I. Optimization and Prediction of Tribological Behaviour of Filled Polytetrafluoroethylene Composites Using Taguchi Deng and Hybrid Support Vector Regression Models. Sci. Rep. 2022, 12, 10393. [Google Scholar] [CrossRef]
  26. Khalilpourazari, S.; Pasandideh, S.H.R. Multi-Objective Optimization of Multi-Item EOQ Model with Partial Backordering and Defective Batches and Stochastic Constraints Using MOWCA and MOGWO. Oper. Res. 2020, 20, 1729–1761. [Google Scholar] [CrossRef]
  27. Khalilpourazari, S.; Naderi, B.; Khalilpourazary, S. Multi-Objective Stochastic Fractal Search: A Powerful Algorithm for Solving Complex Multi-Objective Optimization Problems. Soft Comput. 2020, 24, 3037–3066. [Google Scholar] [CrossRef]
  28. Myszkowski, P.B.; Laszczyk, M. Diversity Based Selection for Many-Objective Evolutionary Optimisation Problems with Constraints. Inf. Sci. 2021, 546, 665–700. [Google Scholar] [CrossRef]
  29. Devaraj, R.; Mahalingam, S.K.; Esakki, B.; Astarita, A.; Mirjalili, S. A Hybrid GA-ANFIS and F-Race Tuned Harmony Search Algorithm for Multi-Response Optimization of Non-Traditional Machining Process. Expert Syst. Appl. 2022, 199, 116965. [Google Scholar] [CrossRef]
Figure 1. Flowchart for the execution of the MFO algorithm.
Figure 1. Flowchart for the execution of the MFO algorithm.
Jcs 07 00462 g001
Figure 2. Abrasive waterjet cutting on fabricated fibre intermetallic laminates.
Figure 2. Abrasive waterjet cutting on fabricated fibre intermetallic laminates.
Jcs 07 00462 g002
Figure 3. Statistical validation of proposed models: (a) MRR, (b) FD, and (c) KD.
Figure 3. Statistical validation of proposed models: (a) MRR, (b) FD, and (c) KD.
Jcs 07 00462 g003aJcs 07 00462 g003b
Figure 4. Effect of AWJ cutting parameters and wt% of r-GO on the response characteristics: (a,b) MRR, (c,d) FD, and (e,f) KD.
Figure 4. Effect of AWJ cutting parameters and wt% of r-GO on the response characteristics: (a,b) MRR, (c,d) FD, and (e,f) KD.
Jcs 07 00462 g004
Figure 5. Convergence plots for: (a) MRR, (b) FD, and (c) KD obtained through different metaheuristic algorithms.
Figure 5. Convergence plots for: (a) MRR, (b) FD, and (c) KD obtained through different metaheuristic algorithms.
Jcs 07 00462 g005
Figure 6. Normal probability plots for diversity values of: (a) GA, (b) PSO, (c) MFO, (d) GWO, and (e) DFO.
Figure 6. Normal probability plots for diversity values of: (a) GA, (b) PSO, (c) MFO, (d) GWO, and (e) DFO.
Jcs 07 00462 g006
Figure 7. Friedman’s statistical results for diversity values obtained through different algorithms.
Figure 7. Friedman’s statistical results for diversity values obtained through different algorithms.
Jcs 07 00462 g007
Figure 8. Normal probability plots for spacing values of: (a) GA, (b) PSO, (c) MFO, (d) GWO, and (e) DFO.
Figure 8. Normal probability plots for spacing values of: (a) GA, (b) PSO, (c) MFO, (d) GWO, and (e) DFO.
Jcs 07 00462 g008
Figure 9. Friedman’s statistical results for spacing values obtained through different algorithms.
Figure 9. Friedman’s statistical results for spacing values obtained through different algorithms.
Jcs 07 00462 g009
Figure 10. Normal probability plots for hypervolume values of: (a) GA, (b) PSO, (c) MFO, (d) GWO, and (e) DFO.
Figure 10. Normal probability plots for hypervolume values of: (a) GA, (b) PSO, (c) MFO, (d) GWO, and (e) DFO.
Jcs 07 00462 g010
Figure 11. Friedman’s statistical results for hypervolume values obtained through different algorithms.
Figure 11. Friedman’s statistical results for hypervolume values obtained through different algorithms.
Jcs 07 00462 g011
Table 1. Process parameters for AWJ cutting experimentation.
Table 1. Process parameters for AWJ cutting experimentation.
S. No.Dependent VariablesCodeLevels
LowMediumHigh
1Weight percentage of r-GO (wt%)r-GO012
2Traverse speed (mm/min)TS400500600
3Waterjet Pressure (MPa)WP200250300
4Nozzle height (mm)NH234
Table 2. Experiment details based on RSM-BBD and their measured responses.
Table 2. Experiment details based on RSM-BBD and their measured responses.
Exp. No.Dependent VariablesOutput Responses
r-GO (wt%)Traverse Speed (mm/min)Waterjet Pressure (MPa)Nozzle Height (mm)Material Removal Rate (g/s)Delamination FactorKerf Deviation (mm)
1040025033.0892.1411.857
2240025033.3382.1321.937
3060025033.1051.8221.825
4260025033.1152.3171.851
5150020023.2332.092.148
6150030022.9932.0482.136
7150020043.3232.1572.196
8150030043.1362.2481.897
9050025022.9941.9021.857
10250025022.9452.0062.201
11050025042.9371.9242.057
12250025043.1732.3011.802
13140020033.3892.1962.094
14160020033.4572.2172.019
15140030033.3292.3221.937
16160030033.0812.1531.889
17050020033.3291.91.977
18250020033.3962.1752.245
19050030033.0581.9422.057
20250030033.2332.1961.884
21140025023.0092.1931.977
22160025023.0732.0061.958
23140025043.2452.2591.912
24160025042.9772.281.826
25150025033.0452.182.126
26150025033.0572.192.117
27150025033.0592.2132.116
28150025033.0732.2042.112
29150025033.0792.1882.122
Table 3. Statistical analysis of MRR.
Table 3. Statistical analysis of MRR.
SourcesSSDFMSFProb > F
Model0.6130140.04379224.815<0.0001Significant
A-r-GO0.0394410.03944202.510<0.0001
B-TS0.0291010.02910149.432<0.0001
C-WP0.1401810.14018719.696<0.0001
D-NH0.0246610.02466126.609<0.0001
AB0.0142810.0142873.313<0.0001
AC0.0029110.0029114.9700.0017
AD0.0203010.02030104.251<0.0001
BC0.0249610.02496128.163<0.0001
BD0.0275510.02755141.470<0.0001
A20.0020610.0020610.5800.0058
B20.0412010.04120211.532<0.0001
C20.1951410.195141001.864<0.0001
D20.0286210.02862146.934<0.0001
Residual0.00272140.00019
Lack of Fit0.00199100.000191.0910.5085Insignificant
Pure Error0.0007340.00018
R299.55% Adj. R299.11%
Table 4. Statistical analysis of FD.
Table 4. Statistical analysis of FD.
SourcesSSDFMSFProb > F
Model0.51333140.03666223.756<0.0001Significant
A-r-GO0.1865010.186501138.112<0.0001
B-TS0.0167210.01672102.065<0.0001
C-WP0.0025210.0025215.3960.0015
D-NH0.0711410.07114434.175<0.0001
AB0.0635010.06350387.528<0.0001
AD0.0186310.01863113.702<0.0001
BC0.0090210.0090255.074<0.0001
BD0.0108110.0108166.003<0.0001
CD0.0044210.0044226.9860.0001
A20.1008710.10087615.606<0.0001
B20.0080710.0080749.301<0.0001
C20.0012610.001267.7120.0148
D20.0118310.0118372.200<0.0001
Residual0.00229140.00016
Lack of Fit0.00159100.000150.90350.5945Insignificant
Pure Error0.0007040.000176
R299.54% Adj. R299.12%
Table 5. Statistical analysis of KD.
Table 5. Statistical analysis of KD.
SourcesSSDFMSFProb > F
Model0.48251140.034461343.006<0.0001Significant
A-r-GO0.0070010.00700273.089<0.0001
B-TS0.0099710.00997388.742<0.0001
C-WP0.0643810.064382508.923<0.0001
D-NH0.0287110.028711118.886<0.0001
AB0.0007210.0007228.4060.0001
AC0.0486210.048621894.559<0.0001
AD0.0897010.089703495.301<0.0001
BC0.0001810.000187.1010.0185
BD0.0011210.0011243.730<0.0001
CD0.0205910.02059802.407<0.0001
A20.0602010.060202346.004<0.0001
B20.1555110.155516060.049<0.0001
C20.0025310.0025398.923<0.0001
D20.0126810.01268494.165<0.0001
Residual0.00035142.5 × 10−5
Lack of Fit0.00024102.4 × 10−50.8050.6457Insignificant
Pure Error0.0001140.00002
R299.92% Adj. R299.85%
Table 6. Parameters and their values used in GA, PSO, MFO, DFO, and GWO.
Table 6. Parameters and their values used in GA, PSO, MFO, DFO, and GWO.
AlgorithmParametersValue/Range of Parameters
GAMethod of Selection for ReproductionRoulette Wheel Selection Method
Cross over probability 0.5
Cross over operatorSingle Point Cross Over Technique
Mutation probability0.04
Mutation operatorRight-side swapping
Replacement strategy 100% replacement Strategy
PSOLearning factors (C1 and C2)2 and 2
Inertia weight (ω)0.9
MFOPosition of moth close to the flame (t)−1 to −2
Update mechanism Logarithmic spiral
DFOInertia weight (IW) 0.2
Separation weight s w = 0.1 0.1     i t r n i t r
Alignment weight a w = 0.1 0.1     i t r n i t r
Cohesion weight a w = 0.1 0.1     i t r n i t r
Food factor f f = 2 r a n d
Enemy factor e f = 0.1 0.1     i t r n i t r
GWOScale Factor (SF)2
Table 7. Optimal AWJ cutting parameter achieved through different algorithms.
Table 7. Optimal AWJ cutting parameter achieved through different algorithms.
Algorithmsr-GOTSJPNHMRRFDKD
GA0.004561.74277.892.022.9631.9451.972
PSO0.038592.68261.922.073.0451.9511.893
MFO0.000600.00253.362.003.0961.9281.833
GWO0.026595.67255.662.043.0741.9471.865
DFO0.000600.00246.802.003.1271.9351.819
Table 8. Diversity values for different algorithms.
Table 8. Diversity values for different algorithms.
R. No.GAPSOMFOGWODFO
12.1686492.3048932.4996542.0438352.375572
22.2134752.2186712.4793632.2280792.171757
31.9690532.2176242.4955592.0092522.135089
42.1819772.2282352.4634192.0145812.228036
51.9551332.2894232.4954922.1869531.928915
62.014372.2627222.4500181.9984522.443005
72.2275462.2893072.496162.1971052.282027
82.3582512.089062.4991232.2377592.282308
92.2018082.3251282.4716082.186952.29433
102.1705612.2238382.4931662.2626422.189975
112.3686792.3618172.5020112.2144972.388704
122.1496182.2179162.4653042.1260081.981838
132.2675972.3165092.483312.1259322.184926
142.1000452.2805082.5045732.1552232.269794
152.2317282.3032612.4927612.3197212.492793
162.1500872.067872.4905072.2609062.452572
172.1549842.28252.4916262.282242.437811
182.12292.2576342.490072.0723232.236772
192.2485692.1762582.4870192.214592.217735
201.9590282.1823782.4965552.2496772.430681
212.1760412.2197052.5139952.0991852.336751
222.1822762.246782.505392.042432.468178
232.1521852.2997862.4948642.2422182.505339
242.2221582.1745972.5054022.2216462.289126
252.0580782.2227792.5156582.3182322.426425
262.142182.1956572.4997042.0795012.412103
272.0183712.3081232.4998312.1134462.262402
282.2752562.1660782.4902782.179571.997895
292.1299462.1109842.5096372.1149112.247057
302.1829252.2733632.498271.9826842.411676
312.0454382.174862.50082.1840172.274243
322.0498432.2205352.5110272.0967992.042195
332.2650882.2243452.4913872.201412.337978
342.1152092.1103372.4981012.1449632.450462
352.4492852.1772142.4864782.1073562.183684
361.9869942.1969142.5054732.3455332.080497
Table 9. Statistical analysis results for diversity values (the higher the better).
Table 9. Statistical analysis results for diversity values (the higher the better).
MeanMedianStDevp-ValueA-Squared
GA2.15742.16180.11520.4430.35
PSO2.22832.22330.06930.3020.42
MFO2.49372.49590.0140.0061.11
GWO2.16282.18180.09580.7720.24
DFO2.2822.28220.15170.1650.53
Table 10. Spacing values for different algorithms.
Table 10. Spacing values for different algorithms.
R. No.GAPSOMFOGWODFO
10.1718610.126020.0953470.1402960.22027
20.0775130.1488090.139610.1483310.175924
30.1799120.1290010.1685050.1005610.153444
40.1373940.118050.1401850.1417230.199079
50.0593740.1532990.0951160.1293410.129182
60.0899240.1337240.1385160.1558510.187791
70.0318280.1237910.1199120.1150530.160384
80.1218460.1078090.1151560.0883890.2012
90.1687610.1305160.0876490.1168460.19233
100.1034020.1078630.1213930.1158570.134854
110.1274150.1397830.1383810.0775220.081749
120.0759330.1096260.1265150.1332770.211131
130.1460910.0924080.1154360.1469420.171441
140.1492020.1236440.1253370.1619280.161768
150.1770370.0986770.1366480.1592450.168211
160.1176930.1353840.1088840.1062430.181798
170.1312980.0978190.1145630.1200590.194169
180.0735940.1465260.1280280.1451660.228594
190.1568170.151690.1476220.1217310.16827
200.1326610.1400570.1149890.1326360.17157
210.0878260.1415180.1003010.1295330.153153
220.097490.0984010.082950.1120570.176316
230.0333610.1554460.1370770.0896750.172437
240.0714030.1305270.1159460.0844840.1369
250.0941190.1362010.1425860.0893930.191336
260.0430070.1354370.1063770.126310.11356
270.1163290.1247050.1167740.143020.116459
280.1352790.1173420.1185170.1554350.19377
290.10780.0960620.1118020.1541880.130144
300.1871680.1097870.1220270.1136830.171007
310.1464620.1210020.151110.1267610.136029
320.1883340.1076050.1900770.1050730.169524
330.200720.1589270.1176470.1302790.199426
340.0642190.1033740.1173060.1255850.13606
350.1864870.0955410.1083070.1524310.15704
360.1415570.1013330.1136520.1050750.169872
Table 11. Statistical analysis results for spacing values (the lower the better).
Table 11. Statistical analysis results for spacing values (the lower the better).
MeanMedianStDevp-ValueA-Squared
GA0.120310.124630.046840.7030.26
PSO0.123550.124250.019410.2820.44
MFO0.123060.118050.021430.0770.66
GWO0.1250.126540.023030.5580.3
DFO0.167120.170440.031760.4470.35
Table 12. Hypervolume values for different algorithms.
Table 12. Hypervolume values for different algorithms.
R. No.GAPSOMFOGWODFO
10.0020.0070.0030.0120.004
20.0010.0060.0030.0040.003
30.0010.0040.0030.0050.002
40.0040.0080.0040.0020.006
50.0030.0030.0050.0020.006
60.0020.0010.0020.0040.004
70.0030.0050.0070.0060.005
80.0020.0030.0040.0050.006
90.0010.0010.0060.0070.001
100.0020.0010.0040.0050.001
110.0030.0020.0030.0030.008
120.0040.0010.00700.002
130.0040.0020.0060.0030.002
140.0030.0050.0070.0020
150.0030.0030.0020.0030.005
160.0020.0040.0040.0050.003
1700.0020.0030.0050.009
180.0040.0060.0020.0030.003
190.0040.0060.0040.0060.002
200.0030.0060.0050.0030.003
210.0020.0030.0050.0070.005
220.0020.0050.0090.0030.007
230.0010.0040.0040.0040.002
240.0010.0030.0070.0070.003
250.0020.0070.0040.0040.004
260.0060.0040.0010.0050.003
270.0060.0040.0050.0040.002
2800.0020.0030.0050.003
290.0040.0040.0070.0070.002
300.0020.0010.0010.0050.006
310.0030.0030.0070.0040.007
320.001000.0040.006
330.0020.0040.010.0030.01
3400.0010.0050.0010.004
350.0020.0060.0040.0060.001
3600.0020.0090.0020.002
Table 13. Statistical analysis results for hypervolume values (the higher the better).
Table 13. Statistical analysis results for hypervolume values (the higher the better).
MeanMedianStDevp-ValueA-Squared
GA0.0023610.0020.0015150.0220.88
PSO0.0035830.00350.0020480.1080.6
MFO0.0045830.0040.0023470.1250.58
GWO0.0043330.0040.0021650.0460.75
DFO0.0039440.0030.002390.0240.86
Table 14. Mean rank values of different algorithms in Friedman’s statistical analysis.
Table 14. Mean rank values of different algorithms in Friedman’s statistical analysis.
Performance IndicatorsGAPSOMFOGWODFO
Diversity2.08332.75004.94442.02783.1944
Spacing2.66672.58332.58332.77784.3889
Hypervolume2.11112.91673.43063.50003.0417
Table 15. Confirmation experimental results of AWJ cutting responses at optimal parameters obtained through the MFO algorithm.
Table 15. Confirmation experimental results of AWJ cutting responses at optimal parameters obtained through the MFO algorithm.
ResponsesPredictedExperimentalError %
MRR (g/s)3.0963.1561.94
FD1.9281.9812.75
KD (mm)1.8331.7723.33
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rajamani, D.; Kumar, M.S.; Balasubramanian, E. Multi-Response Optimization of Abrasive Waterjet Cutting on r-GO-Reinforced Fibre Intermetallic Laminates through Moth–Flame Optimization Algorithm. J. Compos. Sci. 2023, 7, 462. https://doi.org/10.3390/jcs7110462

AMA Style

Rajamani D, Kumar MS, Balasubramanian E. Multi-Response Optimization of Abrasive Waterjet Cutting on r-GO-Reinforced Fibre Intermetallic Laminates through Moth–Flame Optimization Algorithm. Journal of Composites Science. 2023; 7(11):462. https://doi.org/10.3390/jcs7110462

Chicago/Turabian Style

Rajamani, Devaraj, Mahalingam Siva Kumar, and Esakki Balasubramanian. 2023. "Multi-Response Optimization of Abrasive Waterjet Cutting on r-GO-Reinforced Fibre Intermetallic Laminates through Moth–Flame Optimization Algorithm" Journal of Composites Science 7, no. 11: 462. https://doi.org/10.3390/jcs7110462

APA Style

Rajamani, D., Kumar, M. S., & Balasubramanian, E. (2023). Multi-Response Optimization of Abrasive Waterjet Cutting on r-GO-Reinforced Fibre Intermetallic Laminates through Moth–Flame Optimization Algorithm. Journal of Composites Science, 7(11), 462. https://doi.org/10.3390/jcs7110462

Article Metrics

Back to TopTop