Abstract

The turning, which consists of the removal of metal from the outer diameter of a rotating cylindrical workpiece, is one of the most common techniques for cutting, especially when finishing the product. The values of the cutting parameters, such as feed rate, cutting speed, and depth of cut, in a turning operation must be selected carefully to improve the profit of operations by enhancing productivity and reducing the total manufacturing cost for each component. A high vibration leads to poor surface finish and reduced productivity and shortens the tool life; therefore, this parameter should be controlled. In this study, an experiment is conducted to investigate the effects of these cutting parameters in the turning process of a workpiece, composed of AISI 1040 steel, using the response surface method. Statistical tools were used to design the experiments. These parameters are optimized by using analysis of variance, regression, and optimization techniques to achieve the condition of minimum vibration and chip frequency, therefore improving the surface roughness after the turning process.

1. Introduction

Machining processes, such as turning, milling, and drilling, are important in manufacturing industries. Turning is a process in which the metal from the outer diameter of the rotating cylindrical workpiece is removed, reducing the diameter of the workpiece, usually to a specific dimension, with a smooth finish on the metal. In manufacturing studies, vibration is one of the main issues that is bounded to machining processes. Vibration is an unavoidable phenomenon in machining that affects the quality of the product, machining accuracy, cutting tool, and machine tool life, increasing operation cost.

The goal of the industry has always been the manufacturing of low-cost and high-quality products in a short time. Modern industries try to enhance these characteristics. To achieve a high cutting performance, the industry needs to operate close to the optimal cutting parameters [1]. Consequently, cutting conditions such as feed rate, speed, and depth of cut in the turning process must be selected appropriately to improve productivity, reduce the total manufacturing cost for each component, or achieve a predefined appropriate criterion. Because of the high cost of digitally controlled machine tools, when compared with their traditional counterparts, their operation is required to be efficiently conducted to manufacture the required product [2].

In the turning process, vibrations occur because of friction between the workpiece and the tool. Consequently, the performance of the machine depends heavily on vibration cutting processes, in which a vibration-free scenario is ideally desired. The rate of deterioration and inaccuracy that increase with the usage of the machine tool can be determined by monitoring vibration [2]. Theories of coupling situation are theories of vibration that have become prominent [3]. Wiercigroch and Budak [4] stated that the coupling situation, resulting from the vibration in the direction of momentum, generates the vibration in the direction of the cutting force and vice versa, leading to a multidirectional path (in the directions of the cutting thrust and cutting). Hamdan and Bayoumi [5] suggested that the increase in the rake angle and flank clearance has increased the stability of the dynamic cutting process because of variation effects. Moreover, a mathematical model, a linear differential equation suitable for cutting speeds, low and high dynamic processes of cutting (considering the frictional force on the flank), and the tool of the tool rake have also been proposed. Moreover, these parameters have affected the strength of the initial friction.

Dimla [6] proposed that the tool wear directly affects the bad edge caused by the vibration of the tool, which is, thereafter, increased. Because of the plastic deformation, the high value of tool vibration leads to a decrease due to the linear increase of the contact area. Meanwhile, Dimla [7] investigated the effect of the cutting conditions on the cutting forces and vibrations resulting from the wear of the cutting tool. The results show that variations in cutting conditions affect the cutting forces and vibration signals. Furthermore, Asiltürk and Akkuş [8] suggested that the feed rate and depth affect the stability of the machines, increasing the vibration and resulting in the bad finishing of the surface. Moreover, the increase in speed led to a reduction in the vibration, improving the stability of the manufacturing process.

Accelerometers are commonly used as a sensor for recording vibration signals. Dan and Mathew [9] proposed vibration control techniques [1014] that are more practical and cost-effective and discussed their limitations and advantages. Sick [15] also reviewed various vibration control techniques to predict the flank wear [1618] at the preprocessing level.

Ambhore et al. [19] developed a mathematical prediction model to determine vibration acceleration and surface roughness with specific machining parameters such as cutting speed, feed rate, and cutting depth. The results show that the conditions under which the vibration signal is cut have a significant mixed effect. For the estimation of vibration and surface roughness, the regression model was found to be appropriate. Sofuoğlu et al. [20] used an advanced machining method to research the effects of cuts, cutting velocity and tool overhang lengths, surface ruggedness, stable cuts depth, and tool temperature in turning operations of Hastelloy-X alloys. Çakır et al. [21] studied the machining characteristics of Hastelloy-X and developed a numerical model for the turning of Hastelloy-X. The operation varied from two input parameters (cutting speed and feed rate) to results based on process outputs such as cutting forces, cutting temperature, effective stresses, and chip morphology. The results were evaluated. The results showed that the model was appropriate for turning Hastelloy-X at low and medium speeds of processing. Sofuoğlu et al. [22] carried out 2D finite element (FE) analysis to investigate the effects of these machining methods on titanium and Hastelloy-X alloys in terms of cutting forces, cutting tool temperatures, and effective stresses. DEFORM-2D software was used during analyses. They confirmed that the hot ultrasonic-assisted turning technique reduces cutting forces and effective stress significantly but cutting temperature increases compared to conventional and ultrasonic-assisted turning. A new method of multicriteria decision-making (MCDM) has been developed by Gürgen et al. [23], to investigate the various turning procedures (conventional, ultrasonic-assisted, and hot ultrasonic-assisted) of Ti6Al4V alloy. According to MCDM results, higher cutting speeds, lower tool overhang lengths, and ultrasonic-assisted turning are the appropriate levels of the cutting parameters. The findings have contributed to understanding the MCDM experience of modern nontraditional methods of machining. Parida et al. [24] studied the effect of microdimple-textured tool in turning of Ti6Al4V alloy which is assessed by considering cutting force, chip morphology, and chip flow angle using the finite element method (FEM). The proposed 3D FE model helps in reducing the number of experimental trials that require optimizing the machining process.

The turning process of Inconel 718 was investigated in a study by Guergen et al. [25] in terms of workpiece surface quality and the wear behavior of cutting tools. The study showed the correlation between input and output parameters such as feed rate, insert nose radii, and the insert of coating methods. The analysis gave short results for the proper selection of instruments in Inconel 718’s turning process. In the study of Şahinoğlu and Rafighi [26], turning tests were conducted on three different hardened AISI 1040 steels at three different cutting depths, feeding rates, and noncoolant cutting speeds. The study aimed mainly at minimizing both surface ruggedness and power consumption by estimating optimized machining parameters. The effects of cutting parameters and workpiece hardness on surface roughness, sound level, and power consumption were examined. The findings showed that the feed rate was the most significant machining parameter affecting surface roughness and power consumption. The depth of cut and cutting speed also have significant effects following the feed rate.

The study of Haque et al. [27] found that feed rate, depth of cut, and cutting condition with water are the most significant factors affecting the surface roughness for turning of AISI 1040 steel. Prakash and Krishnaiah [28], Pavani et al. [29], and Nair et al. [30] studied the effect of cutting speed, feed rate, and depth of cut on surface roughness for turning of AISI 1040 steel. However, Nair et al. [30] suggested the same three factors to study the effect of tool vibration on flank wear and surface roughness during high-speed machining of 1040 steel. Kumar et al. [31] carried out Taguchi’s experimental design to obtain the effect of cutting force, cutting tool temperature, and tool flank wear on surface roughness in turning of AISI 1040 steel.

In this work, the response surface method (RSM) was used to obtain the optimal value of response variables (vibration and chip frequency). RSM is a useful set of mathematical and statistical techniques for modeling and analyzing the effects of any variable on a specific problem, to improve the response, as conducted by Montgomery [32]. RSM is widely used in process improvement experiments regarding the turning process. Gaitonde et al. [33] analyzed the effects generated by the depth of the cut and machining time with this method in the turning of AISI D2 cold work tool steel. Moreover, by using the RSM, Asiltürk and Neseli [34] obtained a model for the surface roughness, as a function of the cutting parameters, and analyzed the impact of these parameters on the surface quality. Camposeco-Negrete [35] optimized variables related to the turning of AA 6061 T6 to obtain the conditions for the minmum energy consumption. Bhushan [36] used the RSM analysis to investigate the effects of the cutting parameters in the turning of SiC 7075 Al alloy to reduce the energy consumption of the machine and increase the life of the tool.

According to the literature review, the RSM has been proved to be practical and effective; therefore, it was used in this research to determine the effects of parameters of the turning process on the surface roughness. Furthermore, previous research concluded that the selection of the cutting parameters is essential for machining, especially in the turning process; however, the combination of multiple effects from different parameters on vibration and chip frequency has not been considered yet for the turning process. An experiment was, therefore, conducted to investigate the effects of depth of cut, cutting speed, and feed rate in the turning process of the AISI 1040 steel using RSM.

It can be seen from the literature that few papers have been reported on the optimization of the factors affecting vibration characteristics of turning process of AISI 1040 steel by using RSM in terms of experimental and mathematical analysis. Furthermore, none of the papers considers the combination of the effect of the real conditions such as cutting speed, feed rate, and depth of cut to study the vibration and chip frequency simultaneously by applying the statistical analysis and heuristic methods. Therefore, this study has contributed to studying the effect of cutting speed, feed rate, and depth of cut on vibration and chip frequency simultaneously in turning of AISI 1040 steel using RSM, and it found the optimized parameters using two heuristic algorithms.

The methodology used to investigate the parameter optimization of the turning process of the AISI 1040 steel using the RSM is shown in Figure 1. Furthermore, the proposed methodology steps are defined as follows:Step 1: literature review on the turning process, vibration on machining, RSM, and required terms.Step 2: designing of the experiment. In this step, the RSM is designed to evaluate the effects of parameters on the surface roughness after the turning process.Step 3: experimental processes consist in selecting the appropriate material, which is AISI 1040 steel, and executing the turning process using a Harrison M390 machine.Step 4: data analysis using the RSM to optimize the parameters of the turning process.Step 5: results and discussion to verify if the implementation of the experiment was satisfactory.

2. Experimental Setup for Machining

2.1. Material

In this research, a commercial and commonly used AISI 1040 steel [26, 30, 37], whose chemical composition and mechanical properties are shown in Table 1, was used.

2.2. Mounting Workpiece on Lathe Machine

The experiment has been implemented using a lathe machine Harrison M390, as shown in Figure 2, which has 5.5 kW of nominal power and 2000 rpm of maximum spindle speed. The material of the workpiece was a medium AISI 1040 steel, as described in the previous section. The length and diameter of the workpiece were 500 and 100 mm, respectively.

The workpiece was placed on the tailstock, which needed more support, as the spindle bearings could be harmful to the experiment because of the heaviness of the workpiece. After completing the setup of the lathe machine, the surfacing of the workpiece was conducted to eliminate the outer layer of the workpiece, which was rusty, as shown in Figure 3.

2.3. Tool Holder and Insert

Figure 4 shows the selected insert, which is in the steel cutting grade. Moreover, tungsten carbide and titanium carbide of the insert bound together because of cobalt, CO. The insert was made with face-centered cubic structures to enhance the cutting speed and precision of the process, therefore improving the smoothness of the workpiece. Because of their unique properties, such as low friction and high resistance against crack, diffusion, and wear, coated tools can be used at high cutting speeds, reducing the time and cost required for machining operations.

After the insert was clamped to the tool holder, the clamp was tightened with a clamp screw, as shown in Figure 5. Two contcgqzurations of tool overhang (70 and 120 mm) were used. Afterward, the tool holder was tightened on the turret and the cutting tool was accurately centered, considering the center of the tailstock tip, to avoid imprecision problems in the process, which could hinder the smoothness of the final workpiece.

3. Design of the Experiment Using RSM

The RSM was used to estimate the significance of the effects of the process parameters considering the surface roughness. Speed, feed rate, and depth of cut, which are independent variables, were selected as these parameters based on the RSM methodology and according to the previous literature. Figure 6 shows the dependence relationship between the independent variables and dependent variables (vibration and chip frequency variables), which are essential to analyze the surface roughness. The expertise of the user is classified as a moderate variable which has a strong contingent effect on the relationship of the variables. By contrast, the accuracy is used as a mediating variable that influences the results.

In this research, central composite design (CCD) was the methodology chosen to design the experiment. This methodology can be coupled with RSM to build a quadratic model for the variable response.

The CCD is the most popular among the categories of RSM design because of the following reasons:(i)CCD can be divided into two group analyses that represent different effects and, therefore, can be individually performed. The first one estimates the linear effect, and the second one estimates the curvature effect. Consequently, when the analysis data points of the first subset have no significant curvature effect, for example, the second subset is not required.(ii)CCD is effective, as it provides information on experiment variables and overall experimental error with minimum required interactions.(iii)CCD can be used in various experimental areas. CCD contains a design factorial or fractional factorial of the central points that are supported with a set of optimum points, allowing an accurate estimation of the curvature.

In CCD, two types of experimental designs can be chosen (full or small). For this experiment, the small CCD was chosen, as shown in Table 2.

3.1. Vibration Detector Setup

The easy-to-use DasyLab was applied to easily and quickly analyze vibration data by working with a display directly on the screen. An accelerometer was used to sense and obtain vibration data. The accelerometer will be placed on the part to be tested. Accelerometers are installed using special wax because they were a good mounting tool for lightweight sensors in temporary installations. The wire connected to the accelerometers must be properly connected to support the wire weight, as shown in Figure 7(a). Unit icons were placed on the screen and connected to wires in a diagram representing the flow of data across the system. In this case, the vibration sensor was connected to the tool holder and connected to the computer. The captured and analyzed signal was further adapted using the DasyLab 5.6 software in the Fast Fourier transformation (FFT) field, as shown in Figure 7(b). From the FFT plots of the vibration data (the sample is shown in Figure 7)(b), the vibration amplitude was easily determined from the program. The dynamics of the turning process is as follows: the tangential or cutting force runs down at the tip of the tool allowing the workpiece to deviate upward, the feed rate of axial force works in the longitudinal direction, the thrust or radial force works in the radial direction and tends to push the tool away from the workpiece [38].

4. Results and Discussion

This section provides the RSM results based on the measurements of vibration and chip frequency. Furthermore, their relation with the cutting parameters is quantified with interaction, main effect, and surface plots. Moreover, mathematical models and optimization are presented. Table 3 presents the configuration of the experiments and the vibration and chip frequency results.

4.1. Development of Vibration Model
4.1.1. Analysis of Variance

The analysis of variance, shown in Table 3, indicates that the model is adequate as the value of the model is lower than 0.05; therefore, the terms selected to compose the model, such as cutting speed (mm/s), depth of cut (mm), feed rate (ft/rot), A2, B2, C2, A B, A C, and B C, could represent the behavior of the vibration. Consequently, the model could accurately predict the vibration behavior because of the high value of R-sq. (adj) (99.98%). R-squared is a statistical measure that shows the proportion of the variance for a dependent variable which is explained by an independent variable or variables in a regression model.

The estimated regression coefficients for vibration, which are presented in Table 4, indicate that their values were significantly different. The main effects can be seen from the value, so if the value is less than 0.05, it can be said that the result is significant. Furthermore, the cutting speed, depth of cut, and feed rate are significant to the response model at 0.05.

4.1.2. Residual Plots for Vibration

Figure 8 shows the residual plots against the percentage and fitted values, the frequency in the dataset, and the observation index. The results were distributed randomly based on the fitted values and observation index. Moreover, the results show that the value of the residuals (0.223) was higher than 0.05 for the vibration. A linear relation is observed when the residuals are plotted against the percent values.

Table 5 shows the contributions of each cutting parameter and their combinations that have been calculated from equation (1). The combination of cutting speed and depth of cut had the highest contribution for vibration, followed by the depth of cut and feed rate, cutting speed, depth of cut, cutting speed and feed rate, and feed rate:

4.1.3. Interaction Plot for Vibration

Figure 9 shows the influences of the depth of cut and the feed rate on the vibration. Consequently, the increase of depth of cut (B) and/or feed rate (C) contributes to the increase of vibration. The interaction between feed rate and depth of cut, as shown in Figure 9, also influences the vibration.

4.1.4. Contour Plot for Interaction Factor

Figure 10 displays the vibration contours at all different combinations of cutting speed parameters. Consequently, the increase in depth of cut led to an increase in vibration increased.

4.1.5. Surface Plot for Vibration

Figure 11 displays a surface plot of the vibration based on the feed rate and depth of cut. The target vibration (red area) is achieved, therefore, at a low value of depth of cut and high value of feed rate.

4.1.6. Mathematical Model of Vibration

The vibration (V) model obtained from the procedure in the previous section is represented in equation (2).

Vibration:

4.2. Development of Chip Frequency Model
4.2.1. Analysis of Variance

The analysis of variance shown in Table 6 shows that each term of the model (A, B, A B, A C, B C, A2, and C2) was important to describe the chip frequency due to its low value. Furthermore, the results of the model indicate a high accuracy as the R-sq. value was high (98.73%).

Table 7 shows the estimated regression coefficients for the chip frequency. Consequently, the values of the coefficients were significantly different. As the values were lower than 0.05 for the linear term, this set of parameters could represent the behavior of the chip frequency. The cutting speed and depth of cut are significant to the response model at .

4.2.2. Residual Plots for Chip Frequency

Figure 12 displays the residual plots for chip frequency against the percentage and fitted values, the frequency of the values in the dataset, and the observation index. The results were distributed randomly based on the fitted values and observation index. Moreover, the results show that the value of the residuals (0.08) was higher than 0.05 for the chip frequency. A linear relation is observed when the residuals are plotted against the percent values.

Table 8 indicates the contribution of each significant parameter and their combination. The cutting speed had the highest contribution for chip frequency, followed by the combination of cutting speed and depth of cut, feed rate, cutting speed and feed rate, and feed rate.

4.2.3. Interaction Plot for Chip Frequency

Figure 13 shows the influences of the combinations of the cutting speed and feed rate, and depth of cut and feed rate on the chip frequency. The interaction plot shows that the increase in cutting speed led to a decrease in chip frequency. By contrast, the chip frequency increased with the depth of cut.

4.2.4. Contour Plot for Chip Frequency

Figure 14 displays the chip frequency contours for all different combinations of cutting parameters. Consequently, the chip frequency increases with the feed rate and decreases with the cutting speed.

Figure 15 shows the chip frequency contours for all different combinations of cutting parameters. Consequently, increasing the depth of cut and decreasing feed rate leads to increase in the chip frequency

4.2.5. Surface Plot for Chip Frequency

Figures 16 and 17 show surface plots for chip frequency based on the combination of cutting speed and feed rate, and depth of cut and feed rate, respectively. In the first case, the target chip frequency can be achieved by decreasing the cutting speed and increasing the feed rate. Meanwhile, in the second case, the target chip frequency can be achieved by decreasing the depth of cut and increasing the feed rate.

4.2.6. Mathematical Model for Chip Frequency

The chip frequency model obtained from the procedure in the previous section is represented in the following equation:

4.3. Optimization of Vibration and Chip Frequency

The goal of optimization is to obtain minimum vibration and chip frequency with the assist of numerical optimization technique by modifying the turning conditions. Thus, this study will solve the optimization problem by combining the RSM vibration and chip frequency models with the tabu search (TS) and simulated annealing (SA) algorithms, which has since been very popular in the solution of optimization problems [39]. TS uses memory (called tabu list) for saving search information, which helps avoid previous solutions that have been visited. The tabu list is updated at each iteration of the TS. Due to restrictions on the tabu list, the risk of rejecting solutions which are not yet generated may occur.

Tabu solutions are therefore inspected for certain condition, known as the criteria for aspiration. If the tabu solution is better than the best objective value found so far, the tabu solution will be accepted and its move will get removed from the tabu list [40]. Like any local search method, TS has various elements to be tuned before the search process starts. Initial solution, tabu list, neighboring structure, aspiration criteria, and stopping criteria are the main elements of TS [41].

In case of a minimization problem, the general TS algorithm can be presented as shown in Algorithm 1 [41].

; %Initial solution
 Initialize the tabu list, medium-term and long-term memories;
Repeat
  Generate a set “” of solutions;
  Find best neighbor of ; %nontabu or aspiration criterion holds
  ;
  Update tabu list and aspiration conditions;
If;
Until stopping criteria satisfied
Output: is the best-found solution.

Simulated annealing is an approach of solving combinatorial optimization problems, and this approach points to a straight analogy in which liquids freeze and crystallize or that metals cool and anneal. Simulated annealing approach is based on the work of Metropolis et al. [42] and Kaid et al. [43] for solving combinatorial optimization problems. Basically, the change in energy of a system when it is converges to a steady “frozen” state by applying a cooling process on that system has been simulated by Metropolis’s algorithm. SA explores the set of all possible solutions, minimizing the chance of being cohesive to local optima by accepting moves that may worsen the value of the objective function to escape from the local optima and move toward a new area in the solution space. A better move is always accepted. Assuming a minimization problem with solution space S, objective function f, and neighborhood structure N, the general SA algorithm can be stated as shown in Algorithm 2 [43, 44].

; % Initial solution
  
  
  
  
  
Output: is the best-found solution.

The optimization of vibration and chip frequency will be formulated in the standard mathematical formats.Vibration”(i)Find: A (cutting speed), B (depth of cut), C (feed rate)(ii)Minimize V (A, B, C)(iii)Subject to(iv)Amin ≤ A ≤ Amax(v)Bmin ≤ B ≤ Bmax(vi)Cmin ≤ C ≤ CmaxChip frequency:(i)Find: A (cutting speed), B (depth of cut), C (feed rate)(ii)Minimize CF (A, B, C)(iii)Subject to(iv)Amin ≤ A ≤ Amax(v)Bmin ≤ B ≤ Bmax(vi)Cmin ≤ C ≤ Cmax

The proposed TS algorithm was implemented on MATLAB R2015a. A PC with Intel (R) Core (TM) i7-4702MQ CPU @ 2.20 GHz, 16 GB RAM, and running on Windows 10, 64-bit operating system. The parameter tuning is an important problem in the design of the TS and SA algorithms. The values of the TS and SA algorithm parameters influence its performance. For this purpose, test runs have been performed for a variety of TS and SA algorithm parameter values, and the best possible setting is selected and presented in Table 9.

By solving the optimization problem, Minitab optimizer, TS, and SA predicted the optimum vibration and chip frequency for the machining of AISI 1040 steel in the selected turning condition range. Table 10 shows the optimal conditions and for minimum vibration and chip frequency. In this condition, the cutting speed was 50 mm/s, depth of the cut was 0.79 mm, and feed rate was 0.80 ft/rot. The optimum vibration and chip frequency were 0.076 Volt and 10193 Hz, respectively. In addition, Figures 18 and 19 show the final graph of MATLAB operation results for the parameter convergence curves in the vibration and chip frequency optimization.

With physical measurements, the optimum conditions predicted by the Minitab optimizer, TS, and SA were further validated. The results of the experiments with optimum turning parameters for the vibration and chip frequency (as predicted by Minitab optimizer, TS, and SA) are good.

5. Conclusion

In this study, the RSM based on TS and SA was applied to predict the influences of key parameters, such as cutting speed, depth of cut, feed rate, vibration, and frequency, on the surface roughness. The workpiece of the experiment was made of a commercial and commonly used AISI 1040 steel, and the experiment was performed using a lathe machine Harrison M390, which has 5.5 kW of nominal power and 2000 rpm of maximum spindle speed. The parameters required to achieve the optimum cutting, regarding surface roughness in a turning process, were obtained through an optimization procedure. The findings and conclusions are presented as follows:(i)The RSM methodology combined with the design of the experiments is a useful technique for analyzing the surface roughness, which is influenced by the vibration and chip frequency of the process.(ii)The results obtained using the RSM provided an effective, systematic, and practical method for improving the turning process.(iii)The cutting speed, depth of cut, and feed rate were statistically associated with vibration and chip frequency. The vibration, therefore, increased when the cutting speed (A) increased from 50 to 200 mm/s or when the depth of cut (B) increased from 1 to 2 mm. Meanwhile, the chip frequency decreased when the cutting speed (A) increased from 50 to 200 mm/s and increased when the depth of cut increased from 1 to 2 mm.(iv)Three-dimensional surfaces generated by RSM can help the user to visualize the effects of the parameters on the surface roughness, considering the specified range.(v)The optimization is generated from the RSM and TS and SA algorithms. The optimal conditions to minimize the vibration and chip frequency were a cutting speed of 50 mm/s, depth of the cut of 0.79 mm, and feed rate of 0.80 ft/rot.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Adel Al-Shayea, Fawaz M. Abdullah, and Mohammed A. Noman were responsible for the conceptualization of the study. Husam Kaid developed the software. Adel Al-Shayea, Fawaz M. Abdullah, Mohammed A. Noman, and Husam Kaid were responsible for the methodology of the study. Fawaz M. Abdullah, Mohammed A. Noman, and Husam Kaid performed the formal analysis and validated the study. Adel Al-Shayea, Fawaz M. Abdullah, and Mohammed A. Noman wrote the original draft. Adel Al-Shayea, Fawaz M. Abdullah, Mohammed A. Noman, Husam Kaid, and Emad Abouel Nasr edited and reviewed the draft. Fawaz M. Abdullah, Mohammed A. Noman, and Husam Kaid were responsible for visualization of the data and for data curation. Adel Al-Shayea supervised the study. Emad Abouel Nasr was responsible for acquisition of the fund. All authors have read and agreed to the publication of the manuscript.

Acknowledgments

This study was supported by Researchers Supporting Project number RSP-2020/164, King Saud University, Riyadh, Saudi Arabia.