Experimental Investigations and Pareto Optimization of Fiber Laser Cutting Process of Ti6Al4V

In the current study, laser cutting of Ti6Al4V was accomplished using Taguchi’s L9 orthogonal array (OA). Laser power, cutting speed, and gas pressure were selected as input process parameters, whereas surface roughness (SR), kerf width, dross height, and material removal rate (MRR) were considered as output variables. The effects of input variables were analyzed through the analysis of variance (ANOVA), main effect plots, residual plots, and contour plots. A heat transfer search algorithm was used to optimize the parameters for the single objective function including higher MRR, minimum SR, minimum dross, and minimum kerf. A multi-objective heat transfer search algorithm was used to create non-dominant optimal Pareto points, giving unique optimal solutions with the corresponding input parameters. For better understanding and ease of selection of input parameters in industry and by scientists, a Pareto graph (2D and 3D graph) is generated from the Pareto points.


Introduction
Titanium and its alloys have a good strength-to-weight ratio that is retained even at elevated temperatures [1][2][3]. They have outstanding corrosion resistance and erosion resistance, owing to which there is a wide variety of applications for these alloys in the domain of aerospace engineering, the chemical and petrochemical industries, the pharmaceutical and food processing sector, and offshore and marine engineering [2]. Ti6Al4V is one of the most widely used alloys of titanium, especially in the aerospace industry. It has a better corrosion resistance compared to materials such as stainless steel and cobalt alloys. They are highly resistant to fatigue and have a satisfactory tensile strength at higher temperatures [3].
The ever-increasing manufacturing competitiveness in the present scenario makes the adoption of newer machining technologies more desirable. Challenges such as achieving higher accuracy, better surface finish, increased production rate, and enhanced product life with minimum impacts on the environment, machine operators, and society are of great importance. Additionally, machining hard materials such as Ti6Al4V also requires intelligent machining strategies. Despite their advantageous properties, machining Ti6Al4V alloy can be challenging due to its low thermal conductivity, high reflectivity, and tendency to form an oxide layer. In addition, the high viscosity of Ti6Al4V molten material, and the formation of a brittle heat-affected zone present challenges [1]. Conventional machining of the performance of laser cutting of Ti6Al4V with laser power, gas pressure, and cutting speed as the input process parameters, while HAZ and the kerf quality of the cut edge were the output variables. HAZ was observed to be increased with increased laser power, was it decreased with increasing values of gas pressure and cutting speed. Kerf width increased with increasing values of power and gas pressure, while it decreased with an increase in speed. An optimal combination of process parameters was observed at a cutting speed of 2400 mm/min, gas pressure of 2 bars, and laser power of 1 kW. Attia Boudjemline et al. [24] performed experiments using a CO 2 laser source for Ti6Al4V with a thickness of 5 mm, with cutting speed and gas pressure as the machining variables. The cutting speed showed a significant effect for controlling the surface roughness of cut specimens while gas pressure did not show any significant effect. Multiple linear regression analysis and ANOVA were performed to validate the model and there was good agreement between the experimental results and the proposed model. L.D. Scintilla et al. [9] analyzed the effect of cutting speed and gas pressure for fiber laser cutting of Ti6Al4V. They concluded that the selection of machining parameters played a major role in achieving better machining quality. Scintilla and Tricarico [25] investigated the fiber laser cutting process using Nd-Yag and CO 2 laser sources to better understand the effects of cutting speed and assist gas pressure on the cutting quality. Higher cutting speed was achieved in fiber laser cutting, along with a lower SR value. Leonardo Daniele Scintilla and Luigi Tricarico [26] performed experiments using a continuous-mode fiber laser cutting process on aluminum, titanium, and magnesium materials. All experiments were performed using 2 kW laser power while varying the cutting speed from 10 m/min to 30 m/min and using a range of assist gas pressures from 0.3 to 0.6 MPa. They aimed to identify the effect of cutting speed and gas pressure on the cut quality. It was observed that dross and SR increased with increasing speed at 0.3 MPa pressure for titanium. Manjoth et al. [27] used Al7075-TiB 2 as the work material for fiber laser cutting. Taguchi's L 9 orthogonal array (OA) was used with standoff distance, cutting speed, and gas pressure as the input variables. SR, MRR, and dimensional accuracy were examined. For SR, cutting speed was observed to be the most significant parameter. For MRR, gas pressure and dimensional error, standoff distance was found to be the most influential factor.
At elevated temperatures, titanium alloys tend to react with oxygen and nitrogen. Proper selection of laser cutting process parameters is essential for obtaining good cut quality and for the prevention of hard and brittle martensite formation, as well as nitride and oxide phases when using nitrogen gas as an assist gas [28]. Previous studies have shown that the proper control of machining variables can eliminate or minimize the oxidation effect of laser cutting of Ti6Al4V [29,30]. An increase in the cutting speed significantly decreases the HAZ. An increase in cutting speed rapidly rejects the molten metal in front of the cutting kerf. This gives a very short time for heat conduction. Thus, HAZ decreases with increases in speed due to the shorter exposure time [31]. In a study conducted by Anderson et al. [32], a smooth cut surface was obtained along with negligible HAZ during laser cutting of Ti6Al4V using nitrogen as an assist gas. Poshyananda et al. [28] studied the fiber laser cutting of Ti6Al4V sheet with nitrogen as an assist gas. Their results revealed that the higher cutting speed resulting in smaller HAZ and higher productivity, with enhanced quality of laser cutting. HAZ can be decreased by increasing the cutting speed. However, an increase in laser power has no significant effect on HAZ. Their results also revealed that a far smaller HAZ effect was found during fiber laser cutting, and that this could be easily eliminated through polishing. B. Yilbas et al. [33] examined the effect of laser cutting of titanium alloy. The use of nitrogen as an assist gas resulted in shallow HAZ because of the low thermal diffusivity of titanium alloy. They also concluded that the use of nitrogen resulted in the formation of TiN x in the cutting section, which in turn increased the brittleness of the recast layer. André Reck et al. [34] examined the effect of laser cutting of Ti6Al4V and the effect of mechanical post-processing. They revealed that optimized parametric settings of laser cutting could significantly eliminate the HAZ. They also analyzed HAZ and found that it did not play a significant role in fatigue failure. SR was found to be more dominant with respect to fatigue failure. Nevertheless, HAZ was entirely removed simply by means of mechanical polishing. Very thin HAZ was obtained during fiber laser cutting as compared to CO 2 and Nd-YAG laser cutting [9]. This is because fiber laser cutting has a high-brightness beam and moves faster, meaning that the material receives very little heat. These past studies show that the proper control of laser cutting process parameters, the use of nitrogen as an assist gas, fiber laser cutting, and mechanical post-processing can significantly eliminate HAZ [28][29][30][31][32][33][34][35][36].
On the basis of the previous literature, it can be seen that fiber laser cutting has not been properly explored with respect to the multi-objective optimization of the machining variables of the fiber laser cutting process. Furthermore, to the best of our knowledge, multi-objective optimization using the HTS algorithm for the fiber laser cutting process parameters of Ti6Al4V has not yet been reported. On the basis of the previous literature, laser power, cutting speed, and gas pressure were selected as the input process parameters, while surface roughness, kerf width, dross height, and material removal rate were considered as output variables in the current investigation. Taguchi's experimental design, which consists of three levels (L 9 OA), was used to statistically evaluate the effect of input parameters on the analyzed output parameters. Analyses of output responses and the effects of input variables on output responses were carried out with the help of ANOVA, main effect plots, residual plots, and contour plots. Simultaneous optimization of output responses was carried out using the HTS algorithm and mathematical models generated on the basis of the regression analysis. A multi-objective heat transfer search algorithm was used to create non-dominant optimal Pareto points, giving unique optimal solutions with the corresponding input parameters. Finally, a validation study was conducted to verify the results obtained from the algorithm.

Materials and Methods
A prime fiber laser cutting machine made by Sahajanand technology was used to perform the experiments on Ti6Al4V. Table 1 shows the chemical composition of Ti6Al4V. The yield strength, ultimate tensile strength, and hardness values of the selected Ti6Al4V alloy were 955 MPa, 990 MPa, and 320 HV, respectively. Machining specimens of dimensions 30 mm × 30 mm × 3 mm were cut using the fiber laser process. A 3 kW laser source was used to conduct the experiments. In order to supply the assist gas (nitrogen as assist gas), a converging nozzle with an outlet diameter of 2 mm was used. Laser power, cutting speed, and gas pressure were selected as the input process parameters, while surface roughness, kerf width, dross height, and material removal rate were considered as the output variables for this investigation. The input parameters, along with their ranges, were selected on the basis of previous literature and preliminary experimental trials. Table 2 shows Taguchi's three-level L 9 orthogonal arrays for the three factors that were used to perform the experiments. Table 2 also shows the constant parameters of fiber laser cutting used in the current study. Taguchi's method has many advantages, such as a reduction in the number of trails, which in turn reduces the cost, time, and resources required during the operations. Taguchi's method enables many process parameters to be considered, and their effects on selected responses of the process to be analyzed. Figure 1 shows the measurement techniques for the selected response variables of the current study.
Surface roughness is a key parameter that defines the quality of the cut edge. The surface roughness was determined using the Surftest SJ-410 model (Mitutoyo, Ahmedabad, India). An enlarged view of the SR measurement is shown in Figure 1a. The arithmetic average roughness (R a ) value was measured from cut specimens in the current study. R a is the arithmetic mean of the peaks and valleys of the surface from the centerline or mean line. The cutoff length λ c was selected as 0.8 mm and the evaluation length was selected as 15 mm during the measurement. Kerf is the amount (width) of material removed from the material during the cutting process. Kerf width measurement was performed using the BX53 microscope (Olympus, Gandhinagar, India). Figure 1b shows the schematic of the kerf measurement. Dross refers to the molten material that is not excluded from the machining zone by gas flow that re-solidifies at the bottom part of the specimen. Dross measurement was performed with help of trinocular microscope MSZ-TR (Magnus, Gandhinagar, India), as shown in Figure 1c. The material removal rate was calculated as the ratio of the weight of the removed specimen to the machining time. In order to achieve higher productivity, MRR must be higher. MRR was measured in g/s, as per Equation (1): weight of the removed material machining time (1)   Table 3 shows the measured values of response variables such as SR, kerf width, dross measurement, and MRR with respect to the input machining variables. SR value   Table 3 shows the measured values of response variables such as SR, kerf width, dross measurement, and MRR with respect to the input machining variables. SR value was observed in the range of 6.83 µm, as a minimum, to 9.92 µm as a maximum on the basis of the nine experimental trials. Similarly, kerf width in the range of 330.81 µm, as a minimum, to 364.48 µm as a maximum, dross value in the range of 0.210 mm, as a minimum, to 0.297 mm as a maximum, and MRR value in the range of 2.1248 g/s, as a minimum, to 2.8877 g/s as a maximum were measured. For the analysis of these measured responses, the statistical software Minitab (Version V17, Gandhinagar, India) was used. Mathematical relations were obtained for each response variable with respect to the machining parameters by using regression analysis. The regression equations for SR, kerf width, dross and MRR are as shown in Equations (2)-(5), respectively:

Regression Equations
Kerf width = 321.27 + 0.0087 ·x 1 − 5.75·x 2 + 2.069·x 3 where x 1 is power, x 2 is speed, and x 3 is gas pressure. Analysis of variance was implemented to understand the significance of machining variables on the output responses. The contribution of the machining variables was evaluated for each response variable with the help of the ANOVA table. Finally, the R-sq. value from the ANOVA table was used to show the significance of the data fitness.

Analysis of Surface Roughness
In order to achieve optimal quality, the lowest possible value of SR is desirable. The ANOVA technique was used to investigate the significance of the selected input variable. Table 4 shows the ANOVA table for SR. The significance of the input variables was determined on the basis of the F-value or the p-value. ANOVA analysis was conducted at a 95% confidence level. At a 95% confidence level, the p-value should be less than 0.05 to indicate that it has a significant effect on the selected response [38]. Table 4 shows that all of the machining parameters have a significant impact on SR, as the p-values for laser power, speed, and gas pressure were found to be 0.045, 0.009, and 0.004, respectively. The difference of less than 20% between the values of R-sq. and Adjusted (Adj.) R-sq.
indicates that the existing model is suitable [39,40]. The R-sq. and Adj. R-sq. values for SR suggest that the model is able to provide the best fit for the existing data. The ANOVA results show that the contribution of gas pressure was the highest, at 63.44%, followed by cutting speed, with 30.3% contribution, and laser power, with 5.97% contribution. A negligible error contribution of 0.29% was observed, indicating that the present study could be effectively used as a basis for future predictions. Figure 2 shows a pie chart of the percentage contribution of machining variables for SR. The standard deviation for SR was observed to be 0.1226. This indicates that the maximum deviation from the mean value of SR is only 0.1226.  The impact of machining variables on SR on the basis of the main effect plot is shown in Figure 3. Main effect plots indicate the mean response of each level of the factors connected by the line. The Y-axis of the plot represents the mean of the means, i.e., the range of SR values obtained for different levels of the input machining parameters. It can be observed that increasing power results in a decrease in SR values, while a further increase in power results in increasing SR values. This is due to the heating supply in the machining zone. With a lower level of heat supply, the melted material has a higher viscosity and is not easily excluded from the machining zone [41]. However, at higher levels of heat supply, the higher temperature gradient may result in more molten material [26] and, as a result, the removal of the molten material from the cutting zone, thereby affecting the surface quality. Increasing cutting speed resulted in lower values of SR, which can mainly be attributed to the shorter exposure periods of the laser beam on the specimen [10,42]. An increase in SR was observed with increasing values of gas pressure. Increasing gas pressure generated more peaks and valleys in the cutting region, which in turn led to increased values of SR [42]. The impact of machining variables on SR on the basis of the main effect plot is shown in Figure 3. Main effect plots indicate the mean response of each level of the factors connected by the line. The Y-axis of the plot represents the mean of the means, i.e., the range of SR values obtained for different levels of the input machining parameters. It can be observed that increasing power results in a decrease in SR values, while a further increase in power results in increasing SR values. This is due to the heating supply in the machining zone. With a lower level of heat supply, the melted material has a higher viscosity and is not easily excluded from the machining zone [41]. However, at higher levels of heat supply, the higher temperature gradient may result in more molten material [26] and, as a result, the removal of the molten material from the cutting zone, thereby affecting the surface quality. Increasing cutting speed resulted in lower values of SR, which can mainly be attributed to the shorter exposure periods of the laser beam on the specimen [10,42]. An increase in SR was observed with increasing values of gas pressure. Increasing gas pressure generated more peaks and valleys in the cutting region, which in turn led to increased values of SR [42]. Figure 4 shows the residual plots for SR. Residual plots are important for verifying the ANOVA results [43]. The given four-in-one residual plots include a normal probability plot, a residual vs. fitted value plot, a histogram, and a residual vs. observation order plot. In the normality test, a straight line is formed by all of the residuals, as shown in Figure 4. It can be concluded that there is no residual clustering, and that these are all distributed normally. In the residual versus fitted plot, the allocation of all the residuals is random, which is a good method of statistical analysis for ANOVA. In the histogram plot, a parabolic structure can be observed, which supports the normality of the data. In the residual versus fitted plot, the absence of any pattern indicates the fitness of the ANOVA results. This verifies that all four tests confirm the satisfactory nature of future outcomes of the proposed model. the removal of the molten material from the cutting zone, thereby affecting the surface quality. Increasing cutting speed resulted in lower values of SR, which can mainly be attributed to the shorter exposure periods of the laser beam on the specimen [10,42]. An increase in SR was observed with increasing values of gas pressure. Increasing gas pressure generated more peaks and valleys in the cutting region, which in turn led to increased values of SR [42].  Figure 4 shows the residual plots for SR. Residual plots are important for verifying the ANOVA results [43]. The given four-in-one residual plots include a normal probability plot, a residual vs. fitted value plot, a histogram, and a residual vs. observation order plot. In the normality test, a straight line is formed by all of the residuals, as shown in Figure 4. It can be concluded that there is no residual clustering, and that these are all distributed normally. In the residual versus fitted plot, the allocation of all the residuals is random, which is a good method of statistical analysis for ANOVA. In the histogram plot, a parabolic structure can be observed, which supports the normality of the data. In the residual versus fitted plot, the absence of any pattern indicates the fitness of the ANOVA results. This verifies that all four tests confirm the satisfactory nature of future outcomes of the proposed model.
. The influences of two simultaneous input variables (determined by keeping the third input variable constant) on SR were determined on the basis of the contour plots. Figure  5 shows the contour plots for SR. In the contour plot for speed and power, the horizontal axis represents power, and the vertical axis represents speed. This contour plot was generated by maintaining the value of gas pressure constant at the middle level (12 bar). The range of SR, which varied from 7.5 to 9.5 µm, is depicted using different colors in Figure  5. The lowest value of SR (<7.5 µm) was obtained by keeping the cutting speed in the range of 4.5 to 5 m/min, and with any value of laser power (between 2000 to 3000 W). The reason for this is that the higher cutting speed results in insufficient exposure of the specimen to the laser beam. The contour plot of gas pressure and power was generated by holding the cutting speed value at the middle level (4 m/min). The minimum value of SR (<7.5 µm) can be obtained when the gas pressure is around 8 bar, and for any value of laser power (between 2000 to 3000 W). In the contour plot of gas pressure and speed, the laser power was kept at a constant value of 2500 W. The lowest SR can be achieved at the lowest value of gas pressure and the maximum value of cutting speed used in the current study. The influences of two simultaneous input variables (determined by keeping the third input variable constant) on SR were determined on the basis of the contour plots. Figure 5 shows the contour plots for SR. In the contour plot for speed and power, the horizontal axis represents power, and the vertical axis represents speed. This contour plot was generated by maintaining the value of gas pressure constant at the middle level (12 bar). The range of SR, which varied from 7.5 to 9.5 µm, is depicted using different colors in Figure 5. The lowest value of SR (<7.5 µm) was obtained by keeping the cutting speed in the range of 4.5 to 5 m/min, and with any value of laser power (between 2000 to 3000 W). The reason for this is that the higher cutting speed results in insufficient exposure of the specimen to the laser beam. The contour plot of gas pressure and power was generated by holding the cutting speed value at the middle level (4 m/min). The minimum value of SR (<7.5 µm) can be obtained when the gas pressure is around 8 bar, and for any value of laser power (between 2000 to 3000 W). In the contour plot of gas pressure and speed, the laser power was kept at a constant value of 2500 W. The lowest SR can be achieved at the lowest value of gas pressure and the maximum value of cutting speed used in the current study.  Table 5 shows the ANOVA results for kerf width. ANOVA analysis was conducted at a 95% confidence level to identify the most significant effects of the input variables on the kerf width. Table 5 shows that the cutting speed and gas pressure have a significant impact on the kerf width, as the p-values for the speed and gas pressure can be observed to be less than 0.05. The difference between the values of R-sq. and Adj. R-sq. was less than 20%, which shows that the existing model is suitable for predicting future outcomes with minimal deviations. The ANOVA results show that the gas pressure has a significant impact on kerf width, with a 53.17% contribution, followed by cutting speed, with a 29.37% contribution, and laser power, with a 15.09% contribution. Figure 6 shows a pie chart of the percentage contribution of machining variables on kerf width. The standard deviation for SR was observed to be 3.036. This indicates that the maximum deviation from the mean value of SR is only 3.036.   Table 5 shows the ANOVA results for kerf width. ANOVA analysis was conducted at a 95% confidence level to identify the most significant effects of the input variables on the kerf width. Table 5 shows that the cutting speed and gas pressure have a significant impact on the kerf width, as the p-values for the speed and gas pressure can be observed to be less than 0.05. The difference between the values of R-sq. and Adj. R-sq. was less than 20%, which shows that the existing model is suitable for predicting future outcomes with minimal deviations. The ANOVA results show that the gas pressure has a significant impact on kerf width, with a 53.17% contribution, followed by cutting speed, with a 29.37% contribution, and laser power, with a 15.09% contribution. Figure 6 shows a pie chart of the percentage contribution of machining variables on kerf width. The standard deviation for SR was observed to be 3.036. This indicates that the maximum deviation from the mean value of SR is only 3.036.  Table 5 shows the ANOVA results for kerf width. ANOVA analysis was conducted at a 95% confidence level to identify the most significant effects of the input variables on the kerf width. Table 5 shows that the cutting speed and gas pressure have a significant impact on the kerf width, as the p-values for the speed and gas pressure can be observed to be less than 0.05. The difference between the values of R-sq. and Adj. R-sq. was less than 20%, which shows that the existing model is suitable for predicting future outcomes with minimal deviations. The ANOVA results show that the gas pressure has a significant impact on kerf width, with a 53.17% contribution, followed by cutting speed, with a 29.37% contribution, and laser power, with a 15.09% contribution. Figure 6 shows a pie chart of the percentage contribution of machining variables on kerf width. The standard deviation for SR was observed to be 3.036. This indicates that the maximum deviation from the mean value of SR is only 3.036.   Figure 7 shows the main effect plot for kerf width. The main effect plot of power vs. kerf shows an increase in the value of kerf width with increasing power. With higher heat supply to the machining zone, radial transfer of heat takes place in the material. Because of this higher energy, more material is vaporized from the cut edge, resulting in a wider kerf width [23]. A continuous decrement in kerf width was observed with increasing cutting speed. This was due to fact that the energy input at the cutting zone decreases with higher travel speed of the laser source [44]. An increase in the value of kerf width can be seen in Figure 7 with increasing gas pressure. Higher gas pressure creates greater force on the molten material, which increases the material removal rate [24].  Figure 7 shows the main effect plot for kerf width. The main effect plot of power vs kerf shows an increase in the value of kerf width with increasing power. With higher hea supply to the machining zone, radial transfer of heat takes place in the material. Because of this higher energy, more material is vaporized from the cut edge, resulting in a wider kerf width [23]. A continuous decrement in kerf width was observed with increasing cut ting speed. This was due to fact that the energy input at the cutting zone decreases with higher travel speed of the laser source [44]. An increase in the value of kerf width can be seen in Figure 7 with increasing gas pressure. Higher gas pressure creates greater force on the molten material, which increases the material removal rate [24].   Figure 8 shows the residual plots for kerf width. All four tests of the residual plo show good agreement, satisfying ANOVA results. This also shows the fitness of the pro posed model for kerf with respect to improving future outcomes.   kerf shows an increase in the value of kerf width with increasing power. With higher supply to the machining zone, radial transfer of heat takes place in the material. Bec of this higher energy, more material is vaporized from the cut edge, resulting in a w kerf width [23]. A continuous decrement in kerf width was observed with increasing ting speed. This was due to fact that the energy input at the cutting zone decreases higher travel speed of the laser source [44]. An increase in the value of kerf width ca seen in Figure 7 with increasing gas pressure. Higher gas pressure creates greater forc the molten material, which increases the material removal rate [24].   Figure 8 shows the residual plots for kerf width. All four tests of the residual show good agreement, satisfying ANOVA results. This also shows the fitness of the posed model for kerf with respect to improving future outcomes.   Figure 9 shows the contour plot of kerf width. The effects of two simultaneous machining variables on kerf width were studied by maintaining the third input variable at a constant value. In the contour plot of speed and power, a constant gas pressure value of 12 bar was used. The range of kerf width obtained for the given levels of the machining variables ranged from 335 to 350 µm. The minimum value of kerf (<335 µm) was obtained with a higher cutting speed (between 4.2 to 5 m/min) along with a lower laser power value (in the range of 2000 to 2250 W). For the contour plot of gas pressure and laser power, the lowest kerf width was obtained with a minimum gas pressure of 8 bar, while maintaining a minimum laser power in the range of 2000 to 2250 W. For the simultaneous effect of gas pressure and speed with a constant laser power of 2500 W, the minimum kerf width was achieved by maintaining a higher cutting speed, between 4.5 to 5 m/min, along with a lower value of gas pressure, 8 bar.

Analysis of Kerf Width
Metals 2021, 11, x FOR PEER REVIEW 12 of 27 Figure 9 shows the contour plot of kerf width. The effects of two simultaneous machining variables on kerf width were studied by maintaining the third input variable at a constant value. In the contour plot of speed and power, a constant gas pressure value of 12 bar was used. The range of kerf width obtained for the given levels of the machining variables ranged from 335 to 350 µm. The minimum value of kerf (<335 µm) was obtained with a higher cutting speed (between 4.2 to 5 m/min) along with a lower laser power value (in the range of 2000 to 2250 W). For the contour plot of gas pressure and laser power, the lowest kerf width was obtained with a minimum gas pressure of 8 bar, while maintaining a minimum laser power in the range of 2000 to 2250 W. For the simultaneous effect of gas pressure and speed with a constant laser power of 2500 W, the minimum kerf width was achieved by maintaining a higher cutting speed, between 4.5 to 5 m/min, along with a lower value of gas pressure, 8 bar.

Analysis of Dross
The ANOVA technique was used for dross, as shown in Table 6, to investigate the significance of the input variables with respect to the developed model. The ANOVA results for dross height indicate that the cutting speed and gas pressure have a significant impact on dross, as the p-values for the cutting speed and gas pressure were 0.039 and 0.023, respectively, which are both below the α value of 0.05. The R-sq. and R-sq. (adj.) values suggest that the model fits the existing data. The difference between both values is less than 20%. The ANOVA results show that gas pressure (56.68%) has the greatest effect on dross, followed by cutting speed (32.22%) and laser power (9.82%). A negligible error contribution of 1.29% was observed, which indicates that the present study can be effectively used as a basis for future predictions. Figure 10 presents a pie chart of the percentage contributions of machining variables with respect to dross. The standard deviation for dross was observed to be 0.005811, indicating that the maximum deviation from the mean value of dross is much lower.  The impact of machining variables on dross is depicted in Figure 11. It can be observed that with increasing laser power, dross height decreases. Further increases in laser power showed the reverse effect, increasing the value of dross height. It can be observed that increasing cutting speed results in decreased dross height. This is due to the heat energy supplied at the cut edge, and is dependent on the laser power and cutting speed [9]. The temperature of the cut edge changes the force magnitude (i.e., the assist gas force for removing the molten material) and direction. At lower values of cut edge temperature, the magnitude of the gas force is lower, and the direction is also away from the centerline. Meanwhile, in cases with higher cut edge temperatures, the magnitude of the gas force is higher, but the direction is far away from the centerline. Therefore, the optimal selection of power and speed will effectively exclude molten material while resulting in minimal dross [9]. An increase in dross was observed with increasing gas pressure. Vibration and unstable conditions were observed at higher gas pressure, resulting in poor machining quality. Additionally, the cooling effect was comparatively greater at higher gas pressure, The impact of machining variables on dross is depicted in Figure 11. It can be observed that with increasing laser power, dross height decreases. Further increases in laser power showed the reverse effect, increasing the value of dross height. It can be observed that increasing cutting speed results in decreased dross height. This is due to the heat energy supplied at the cut edge, and is dependent on the laser power and cutting speed [9]. The temperature of the cut edge changes the force magnitude (i.e., the assist gas force for removing the molten material) and direction. At lower values of cut edge temperature, the magnitude of the gas force is lower, and the direction is also away from the centerline. Meanwhile, in cases with higher cut edge temperatures, the magnitude of the gas force is higher, but the direction is far away from the centerline. Therefore, the optimal selection of power and speed will effectively exclude molten material while resulting in minimal dross [9]. An increase in dross was observed with increasing gas pressure. Vibration and unstable conditions were observed at higher gas pressure, resulting in poor machining quality. Additionally, the cooling effect was comparatively greater at higher gas pressure, which in turn helps in the re-solidification of molten material and results in more dross [45]. which in turn helps in the re-solidification of molten material and results in more dross [45].
. Figure 11. Impact of machining variables on dross. Figure 12 shows the residual plots for dross. All four tests of the residual plot for dross show good agreement similar to SR and kerf width which satisfies the ANOVA results. This also shows a good fit of the proposed model for dross. Figure 11. Impact of machining variables on dross. Figure 12 shows the residual plots for dross. All four tests of the residual plot for dross show good agreement similar to SR and kerf width which satisfies the ANOVA results. This also shows a good fit of the proposed model for dross. Figure 11. Impact of machining variables on dross. Figure 12 shows the residual plots for dross. All four tests of the residual plo dross show good agreement similar to SR and kerf width which satisfies the ANO results. This also shows a good fit of the proposed model for dross.  The influences of two simultaneous input variables (by keeping the third input variable constant) on dross height were determined with the help of contour plots. Figure 13 shows the contour plot for dross height. The contour plot of speed and power was generated by holding the gas pressure value constant at 12 bar. The range of dross height was varied from 0.24 to 0.30 mm. The lowest value of dross (<0.24 mm) was obtained with a maximum cutting speed of 5 m/min along with a laser power in the range of 2000 to 2500 W. The contour plot of gas pressure and power was generated by holding the cutting speed value at the middle level (4 m/min). A minimum value of dross height (<0.24 mm) was obtained with a minimum gas pressure of around 8 bar and any value of laser power (between 2000 and 3000 W). In the contour plot of gas pressure and speed, the laser power was kept at a constant value of 2500 W. The lowest dross height was obtained at the maximum value of cutting speed and the lowest value of gas pressure.

Analysis of Material Removal Rate
A higher value of MRR is always desirable, as it increases productivity, which in turn results in savings with respect to cost, time, and other expenses. Table 7 shows the ANOVA results for MRR. ANOVA analysis was conducted at a 95% confidence level to identify the most significant effects of input variables on MRR. Table 6 shows that cutting

Analysis of Material Removal Rate
A higher value of MRR is always desirable, as it increases productivity, which in turn results in savings with respect to cost, time, and other expenses. Table 7 shows the ANOVA results for MRR. ANOVA analysis was conducted at a 95% confidence level to identify the most significant effects of input variables on MRR. Table 6 shows that cutting speed has a significant effect on MRR, as the p-value for gas pressure was observed to be 0.013, which is less than 0.05. The difference of less than 20% between the values of R-sq. and R-sq. (Adj.) indicates that the present model is suitable for predicting future outcomes with few deviations. The ANOVA result shows that the significance of cutting speed was the highest, with a contribution of 88.81%, followed by laser power, with a contribution of 6.26%, and gas pressure, with a contribution of 3.74%. An error contribution of 1.19% was observed, which shows that the present study can be used as the basis for future predictions. Figure 14 presents a pie chart of the percentage contributions of machining variables towards MRR. The standard deviation for MRR was observed to be 0.05862. This shows that the maximum deviation from the mean value of MRR was only 0.05862.   Figure 15 shows the main effect plot for MRR. The initial increase in laser p sulted in an increase in MRR, while a further increment in power resulted in low Heat supply to the machining zone increases with increasing power, which in tu in material melting in a shorter time, which in turn increases MRR [45]. At th laser power, more material is melted in the machining zone. This melted mate moved from the work specimen, resulting in lower values of MRR at higher value power [45]. A continuous increase in MRR was observed with increasing cuttin MRR is largely dependent on cutting speed. The main effect plot for gas pressure no significant effect on MRR.  Figure 15 shows the main effect plot for MRR. The initial increase in laser power resulted in an increase in MRR, while a further increment in power resulted in lower MRR. Heat supply to the machining zone increases with increasing power, which in turn results in material melting in a shorter time, which in turn increases MRR [45]. At the highest laser power, more material is melted in the machining zone. This melted material is removed from the work specimen, resulting in lower values of MRR at higher values of laser power [45]. A continuous increase in MRR was observed with increasing cutting speed. MRR is largely dependent on cutting speed. The main effect plot for gas pressure showed no significant effect on MRR. laser power, more material is melted in the machining zone. This melted material is removed from the work specimen, resulting in lower values of MRR at higher values of laser power [45]. A continuous increase in MRR was observed with increasing cutting speed MRR is largely dependent on cutting speed. The main effect plot for gas pressure showed no significant effect on MRR.     Figure 17 shows the contour plot of MRR. The effects of two simultaneous machining variables were studied by holding the third input variable constant with respect to MRR In the contour plot of speed and power, a constant gas pressure value of 12 bar was used The range of MRR obtained for the given levels of the machining variables ranged from 2.2 to 2.8 g/s. The maximum MRR (>2.8 g/s) was attained at a higher cutting speed along with any value of laser power in the range of 2000 to 3000 W. Regarding the contour plo of gas pressure and power, a maximum MRR of about 2.6 g/s was obtained at highe power and lower gas pressure. Regarding the simultaneous effect of gas pressure and speed with a constant laser power of 2500 W, a maximum MRR (>2.8 g/s) could b achieved when maintaining a higher cutting speed of between 4.5 and 5 m/min along with any value of gas pressure between 8 and 16 bar.  Figure 17 shows the contour plot of MRR. The effects of two simultaneous machining variables were studied by holding the third input variable constant with respect to MRR. In the contour plot of speed and power, a constant gas pressure value of 12 bar was used. The range of MRR obtained for the given levels of the machining variables ranged from 2.2 to 2.8 g/s. The maximum MRR (>2.8 g/s) was attained at a higher cutting speed along with any value of laser power in the range of 2000 to 3000 W. Regarding the contour plot of gas pressure and power, a maximum MRR of about 2.6 g/s was obtained at higher power and lower gas pressure. Regarding the simultaneous effect of gas pressure and speed with a constant laser power of 2500 W, a maximum MRR (>2.8 g/s) could be achieved when maintaining a higher cutting speed of between 4.5 and 5 m/min along with any value of gas pressure between 8 and 16 bar.
with any value of laser power in the range of 2000 to 3000 W. Regarding the contour plot of gas pressure and power, a maximum MRR of about 2.6 g/s was obtained at higher power and lower gas pressure. Regarding the simultaneous effect of gas pressure and speed with a constant laser power of 2500 W, a maximum MRR (>2.8 g/s) could be achieved when maintaining a higher cutting speed of between 4.5 and 5 m/min along with any value of gas pressure between 8 and 16 bar.

Optimization Using Heat Transfer Search Algorithm
The HTS algorithm, proposed by Patel and Savsani [18], works on heat transfer due to the interaction between the system molecules and with the surroundings to reach thermal equilibrium. Thermal equilibrium can be achieved as a result of three heat transfer phenomena: conduction, convection, and radiation. To obtain an equilibrium state, the HTS algorithm makes use of the three phases of heat transfer. The conduction, convection and radiation phases have equal opportunity to achieve heat transfer. For each generation, one of the modes of heat transfer is selected randomly. The HTS algorithm starts by selecting a randomly generated population. In a randomly generated population, the system has a number of molecules 'n' and a temperature level. The population is then updated in each generation using one of the randomly selected heat transfer modes. Moreover, the updated solution in the HTS algorithm is accepted if it has a better functional value. Subsequently, the worst solutions of the population are replaced by better solutions [40]. Figure 18 shows the flowchart of the HTS algorithm.

Conduction Phase
The solutions were updated in the conduction phase as per Equations (6) and (7), below: where X j,i is the updated solution; j = 1,2, . . . , n; k is a randomly selected solution; j = k; k ∈ (1,2, . . . , n); i is a randomly selected design variable; i ∈ (1,2, . . . , m); g max is the maximum number of generation specified; CDF is the conduction factor; R is the probability variable; R ∈ {0, 0.3333}; and r i ∈ {0, 1} is a uniformly distributed random number.

Convection Phase
The solutions were updated in the convection phase as per Equations (8) and (9), below: TCF = abs(R − r i ), i f g ≤ g max/COF round(1 + r i ), i f g > g max/COF (9) where X j,i is the updated solution; j = 1,2, . . . , n; i = 1,2, . . . , m. COF is the convection factor; R is the probability variable; R ∈ {0.6666, 1}; r i ∈ {0, 1} is a uniformly distributed random number; X s is the temperature of the surroundings; X ms is the mean temperature of the system; and TCF is the temperature change factor.

Radiation Phase
The solutions were updated in the radiation phase as per Equations (10) and (11), below: where X j,i is the updated solution; j = 1, 2, . . . , n; i = 1, 2, . . . , m; j = k; k ∈ (1, 2, . . . , n) and k is a randomly selected molecule; RDF is the radiation factor; R is the probability variable; R ∈ {0.3333, 0.6666}; r i ∈ {0, 1} is a uniformly distributed random number. During the implementation of the HTS algorithm, all of the objectives (SR, dross, kerf, MRR) were considered to be positive integers. The extreme machining limits of the used set-up were considered during the implementation of the algorithms.

•
Power: 1000 W ≤ Power ≥ 3000 W • Cutting speed: 1 m/min ≤ Speed ≥ 8 m/min • Gas pressure: 6 bar ≤ Gas pressure ≥ 20 bar By considering the extreme machining limits, single-objective optimization was carried out for SR, kerf, dross, and MRR, as seen in Table 8. Minimum surface roughness is always desirable for any machined part in order to best fulfill its function during application. The narrowness of the kerf width is a major advantage of the laser process, as it reduces the amount of scrap metal and enhances the cut quality. Minimum kerf width is essential for any cutting, slitting, or trimming process. The minimum value of dross is desirable for achieving better machining quality. A higher value of MRR is always desirable, as it increases productivity, which in turn saves cost, time, and other expenses. On this basis, single-objective optimization was carried out for minimum SR, minimum kerf, minimum dross, and maximum MRR.
Simultaneous optimization was also carried out for all of the selected output responses-SR, kerf, dross, and MRR-by giving an equal weight of 0.25 to all output responses. Equation (12) presents the objective function for simultaneous optimization.   Simultaneous optimization was also carried out for all of the selected output responses-SR, kerf, dross, and MRR-by giving an equal weight of 0.25 to all output responses. Equation (12) presents the objective function for simultaneous optimization.
The corresponding machining variables of a laser power of 1742 W, a cutting speed of 8 m/min, and a gas pressure of 6 bar were observed to provide the optimum values for The corresponding machining variables of a laser power of 1742 W, a cutting speed of 8 m/min, and a gas pressure of 6 bar were observed to provide the optimum values for the present objective function. The results obtained with these parameters were an SR value of 5.25 µm, a kerf of 302.83 µm, a dross of 0.1684 mm, and an MRR of 3.65. However, different users have different requirements. To tackle such situations, multiobjective optimization along with non-dominated optimal solutions were generated. For multi-objective optimization, the MOHTS technique was implemented. MOHTS is an acronym for "multi-objective heat transfer search", and is the multi-objective version of the HTS algorithm. It is able to handle two-or-more-output responses, and creates nondominant solutions for input variables [46]. The MOHTS algorithm was used to generate 49 non-dominant Pareto optimal points for all of the output responses, as shown in Table 9. All of these points give an optimum solution for SR, kerf, dross, and MRR for individual input variables. Table 9 shows 49 feasible optimal Pareto points. Each optimal solution shown in Table 9 is independent and unique. Therefore, depending on the requirements of the required output response values, users can select different levels of the corresponding input variables. Furthermore, whenever any decision variable creates a conflicting effect with respect to the objective function (i.e., output response), the discrete variation of that variable can be observed to be between the upper and lower limit. If the decision variable does not create a conflicting effect with respect to the output response, a constant value of that variable can be observed at all Pareto points. In the present work, the cutting speed and gas pressure create similar effects (i.e., there is no conflict effect) on the output response, and hence a constant value for both parameters (cutting speed of 8 m/min and gas pressure of 6 bar) can be observed at all Pareto points. Validation studies were conducted by randomly selecting five Pareto optimal points to verify the adequacy of the HTS algorithm. Table 10 presents the confirmatory trials conducted for the five Pareto optimal points, along with the experimentally measured values for all of the response variables. A negligible difference can be observed between the experimentally measured values and the values predicted by the HTS algorithm. This confirms that the existing HTS algorithm, along with the proposed model, is highly capable and suitable for fiber laser cutting of Ti6Al4V. 3D and 2D plots were prepared using non-dominant optimal Pareto points, as shown in Table 9. As the minimum values of SR, kerf, and dross are desirable for the laser machining of Ti6Al4V, a 3D plot for these three output response variables was generated, as shown in Figure 19. In the 3D Pareto plot, the X-axis represents the kerf values, the Y-axis represents the dross values, and the Z-axis represents the SR values. However, the values represented in the 3D plots include the effect of the remaining output response variables, i.e., MRR. Each point represented on the 3D Pareto curve gives a unique solution, with its own corresponding input process parameters. The operator can select any Pareto point on the curve, depending on the required values of output responses, to obtain the corresponding values of input process parameters that can be used for machining.
By considering the need for industrial users/operators, 2D plots were also generated. Figure 20a-c show three different 2D plots, corresponding to MRR vs. kerf, MRR vs. SR, and MRR vs. dross. However, the values presented in these 2D plots include the effect of the other two response variables. 2D Pareto curves are more useful when the other two response variables are not desirable for the users. Figure 20a represents the 2D Pareto curve of MRR vs. kerf, in which the X-axis represents the kerf values and the Y-axis represents the MRR values. The machining parameter corresponding to the red point mentioned on the curve indicates that the maximum MRR can be achieved under this machining condition. However, at the same time, the kerf value will be increased, which is not desirable. A similar situation can be observed for the machining parameters, represented by a yellow point, which can provide the desired minimum kerf value, while resulting in decreased values of MRR, which is not desirable. Therefore, the operator can select any optimal point of the curve on the basis of the required values of MRR and kerf. Figure 20b shows the 2D Pareto curve of MRR vs. SR, in which the X-axis represents the SR values and the Y-axis represents the MRR values. The green point on the curve describes the machining conditions under which the desired higher MRR and lower SR can be achieved simultaneously. Figure 20c shows the 2D Pareto curve of MRR vs. dross. The selection of the Pareto point for maximum MRR will result in a higher value of dross (represented by the red point). On the other hand, the Pareto point for lower dross will also result in a decreased value of MRR (represented by the yellow point). Therefore, again depending on the required values of MRR and dross, the optimal point on the curve can be selected, which provides the corresponding machining variables.  By considering the need for industrial users/operators, 2D plots were also generated. Figure 20a-c show three different 2D plots, corresponding to MRR vs. kerf, MRR vs. SR, and MRR vs. dross. However, the values presented in these 2D plots include the effect of the other two response variables. 2D Pareto curves are more useful when the other two response variables are not desirable for the users. Figure 20a represents the 2D Pareto curve of MRR vs. kerf, in which the X-axis represents the kerf values and the Y-axis represents the MRR values. The machining parameter corresponding to the red point mentioned on the curve indicates that the maximum MRR can be achieved under this machining condition. However, at the same time, the kerf value will be increased, which is not desirable. A similar situation can be observed for the machining parameters, represented by a yellow point, which can provide the desired minimum kerf value, while resulting in decreased values of MRR, which is not desirable. Therefore, the operator can select any optimal point of the curve on the basis of the required values of MRR and kerf. Figure 20b shows the 2D Pareto curve of MRR vs. SR, in which the X-axis represents the SR values and the Y-axis represents the MRR values. The green point on the curve describes the machining conditions under which the desired higher MRR and lower SR can be achieved simultaneously. Figure 20c shows the 2D Pareto curve of MRR vs. dross. The selection of the Pareto point for maximum MRR will result in a higher value of dross (represented by the red point). On the other hand, the Pareto point for lower dross will also result in a decreased value of MRR (represented by the yellow point). Therefore, again depending on the required values of MRR and dross, the optimal point on the curve can be selected, which provides the corresponding machining variables.

Conclusions
In the current study, the influence of different process parameters of the fiber laser cutting process was analyzed with respect to the output response variables of Ti6Al4V. Laser power, cutting speed, and gas pressure were selected as the input process parameters, while SR, kerf width, dross height, and MRR were considered as the output variables.

Conclusions
In the current study, the influence of different process parameters of the fiber laser cutting process was analyzed with respect to the output response variables of Ti6Al4V. Laser power, cutting speed, and gas pressure were selected as the input process parameters, while SR, kerf width, dross height, and MRR were considered as the output variables. The following conclusions can be drawn from the current study:

•
The ANOVA results of SR show that all three machining variables were found to be significant, with the highest contribution of 63.44% being for gas pressure, followed by the cutting speed, with a 30.3% contribution, and laser power, with a 5.97% contribution. The ANOVA results of kerf demonstrated the significant effect of cutting speed (53.17%) and gas pressure (29.37%). For dross height, the cutting speed and gas pressure were observed to have a significant effect, while for MRR, the cutting speed had a significant effect. • Regression equations were generated for each response variable that was found to be robust from ANOVA. This model is helpful for the prediction of output responses for the manufacturing of components in industry without conducting the actual experiments.

•
The residual plots for all response variables verified the fitness of the proposed model for achieving a better future outcome. The significance of the machining parameters for each response was analyzed using contour plots. • A difference of less than 20% between the values of R-sq. and Adj. R-sq. was observed for all responses, suggesting that the proposed models achieve the best fit for the existing data. • Single-objective optimization results revealed the optimal process parameter settingsmaximized MRR (3.738 g/s) and minimized SR (5.13 µm)-to be obtained at a laser power of 3000 W, a cutting speed of 8 m/min and a gas pressure of 6 bar.

•
Single-objective optimization results demonstrated that the optimized process parameter settings achieving the minimum values of kerf width (296.384 µm) and dross height (0.1661 mm) were a laser power of 1000 W, a cutting speed of 8 m/min, and a gas pressure of 6 bar. • Simultaneous-objective optimization results showed that a maximum MRR of 3.654 g/s, a minimum SR of 5.25 µm, a minimum kerf width of 302.839 µm, and a minimum dross height of 0.1684 mm could be obtained at a laser power of 1742 W, a cutting speed of 8 m/min, and a gas pressure of 6 bar.

•
The MOHTS algorithm was used to generate 49 non-dominant Pareto optimal points for all output responses. All these points provide optimum solutions for SR, kerf, dross, and MRR for individual input variables. Additionally, 3D and 2D Pareto curves were generated to allow better understanding by users. • A negligible difference was observed between experimentally measured values and the predicted values of the HTS algorithm. This confirms that the existing HTS algorithm in addition to the proposed model is highly capable of and suitable for fiber laser cutting of Ti6Al4V.
Funding: This research received no external funding.
Data Availability Statement: Data presented in this study is available in this article.

Conflicts of Interest:
The authors declare no conflict of interest.