Close banner

2022-11-07 16:27:59 By : Ms. Jenny Shu

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  12, Article number: 10393 (2022 ) Cite this article

This study presents optimization and prediction of tribological behaviour of filled polytetrafluoroethylene (PTFE) composites using hybrid Taguchi and support vector regression (SVR) models. To achieve the optimization, Taguchi Deng was employed considering multiple responses and process parameters relevant to the tribological behaviour. Coefficient of friction (µ) and specific wear rate (Ks) were measured using pin-on-disc tribometer. In this study, load, grit size, distance and speed were the process parameters. An L27 orthogonal array was applied for the Taguchi experimental design. A set of optimal parameters were obtained using the Deng approach for multiple responses of µ and KS. Analysis of variance was performed to study the effect of individual parameters on the multiple responses. To predict µ and Ks, SVR was coupled with novel Harris Hawks’ optimization (HHO) and swarm particle optimization (PSO) forming SVR-HHO and SVR-PSO models respectively, were employed. Four model evaluation metrics were used to appraise the prediction accuracy of the models. Validation results revealed enhancement under optimal test conditions. Hybrid SVR models indicated superior prediction accuracy to single SVR model. Furthermore, SVR-HHO outperformed SVR-PSO model. It was found that Taguchi Deng, SVR-PSO and SVR-HHO models led to optimization and prediction with low cost and superior accuracy.

Filled polymer matrix composites (PMCs) containing fillers continue to receive significant attention from academics and industries due to their modified mechanical and tribological behaviours than virgin polymers1. Polymer based composites showed improved tribological resistance2. Of the different kinds of polymers, polytetrafluoroethylene (PTFE) filled with carbon or bronze fibres are widely used due to their high mechanical and low tribological behaviour3. It has been indicated that these composites are suitable in sectors where mechanical parts including brakes and clutches tribological behaviours are significant4,5,6,7. It has been generally agreed that tribological resistance of materials can be improved by adding more filler content to a certain limit8,9 to neat polymers. Polytetrafluoroethylene (PTFE) has been one of the commonly used thermoplastic matrices for wear conditions because of its low coefficient of friction, ease of process-ability, chemical inertness, low density and low-cost10,11.

Wear is one of the most commonly encountered problems in industries causing frequent substitution of parts especially abrasion. Abrasive wear of various polymers and filled polymers have been studied experimentally. Abrasive wear rate of different matrices were studied by12 and it was found that different polymer exhibited dissimilar wear rate. Inclusion of glass and carbon fabric into vinyl/ester were analyzed. It was reinforced vinyl/ester combination indicated lower wear rate than glass and/or carbon fabric reinforced vinyl/ester composite13. As reported by14 applied load found as the most significant process parameter; reduced wear rate was observed when performance UHMWPE was reinforced with fillers. According to15 it seen that mass loss and µ increased with increase in speed and decrease in grit sizes for betelnut filled epoxy composites.

In order to study multiple responses related to tribological behaviours of composites several decision-making methods including data development, analytic hierarchy as well as grey relational analysis (GRA) have been proposed in the literature16. Of these models, GRA proposed by Deng in 1989 is the widely used methodology especially when the nature of the information is not certain and complete17. Dharmalingam, Subramanian and Kok combined grey relational analysis (GRA) with Taguchi to optimize abrasive tribological property of aluminium hybrid metal composites. Analysis of variance (ANOVA) indicated that grit size was the parameter that had the most influence on wear rate and load was found to had the greatest effect on coefficient of friction18. Sylajakumar et al.19 used Taguchi-GRA method to study the effect of load, speed and distance on coefficient of friction and wear rate of co-long composite. ANOVA showed that speed significantly affect the wear property of the co-continuous composite. Savaran and Thanigaivelan20 optimized dimple geometry and laser parameter using principal component analysis (PCA) coupled GRA. ANOVA showed that average power contributed most while depth contributed less to performance measures. An integrated Taguchi OA and GRA method has been applied to optimize injection moulding parameters of HDPE-TiO2 nanocomposites Pervez et al.21. The work established that optimum parameters were content of TiO2 at 5%, barrel temperature of 225 °C, residence time of 30 min and holding time of 20 s. Adediran et al. optimized mechanical properties of hybrid propylene reinforced bio composites using Taguchi model. It was found that collage of 4% PSS and 10% kenaf fibre produced the optimum combination for hybrid bio composites22. Besides this, Taguchi method hybridized with grey relational grade has as well been employed for multi-response optimization of wire discharge electrical discharge machining23, turning process24 and milling parameter25.

Due to nonlinearity and complex nature of tribological behaviour of materials soft computing methods are increasingly widely accepted including support vector machine (SVM), adaptive neuro-fuzzy inference system (ANFIS) and artificial neural network (ANN). The reason being the fact that these models are capable of capturing the nonlinear and complex nature of the relation between the tribological parameters and responses as compared to conventional mathematical techniques at much cheaper running costs. Various forms of wear are encountered such as abrasion, adhesive, fretting and fatigue wear. Abrasive tribology for composites, instruments, coatings, hip implant, airplane manufacturing as well as automotive components are of essential importance as it determines parts’ performance or longevity. This in general is checked experimentally, as process parameters such as materials characteristics, surface texture, sliding speed and sliding speed. In the analysis of tribology, many mathematical modeling methods have been built. Among them are atomic and molecular kinetics, finite element method, symptom modelling, continuum mechanics, dimension reduction, analysis, boundary element system, stochastic models26. Nevertheless, since tribological behaviours are complex and nonlinear, mathematical models are limited.

Lately, the use artificial intelligent (AI) models has become widely accepted in tribology. Jones et al. pioneered the use of ANN to predict life data and tribological behaviours. Accurate prediction of tribological property by ANN gives an option to the present time, cost and energy consuming testing approaches. Since then, the method has been successfully applied in the tribology discipline that includes wear of reinforced polymer composites27,28, coefficient of friction and mechanical properties, respectively29,30, compensation of magnetic levitation using ANN based on fuzzy inference31. ANFIS and ANN were compared in the prediction of Ks PTFE and its composites. It was found that ANN performed better than ANFIS32. Prediction of abrasive wear of industrial waste and glass filled polyester composites was done using ANN and linear regression model. The results found that ANN outperformed the linear model33. SVM has been employed in the prediction of tool wear appraisal34,35. Also SVM, RBFF and ANN have been contrasted to predict diameter of PCL/gelatin materials. It was reported that ANN did better than SVR and RBFF put together36. Response surface methodology, ANN-HHO as well as model was used in the prediction of abrasive wear of ultrahigh strength martensitic steel. It was reported that hybridized ANN-HHO showed better performance than the single ANN model37.

A survey of the reported literature on the database Scopus yielded the findings that there were 450 peer-reviewed papers starting from 1989 up to date adopted over the literature using the feasibility of wide interest for the abrasive tribological behaviour of PTFE based composites. Figure 1b shows 388 keywords occurrence between those studies, indicating the deep interest and implementation of this field. In addition, the popularity of this study topic was investigated in different regions throughout the world, with the bulk of the countries producing the output being China, the United States, and India. (Fig. 1a). The motivation of this study demonstrated excellent AI techniques for predictions abrasive tribological behaviour of filled PTFE composites. Generally each study has progressed to a little higher degree of accuracy for observations and efficiency at a deeper level than the previous one. To the best knowledge of the authors, no study published in a technical literature has predicted the abrasive tribological behaviour of filled PTFE composites employing this approach using small amount of data. As a result, the goal of this work is to optimize and predict multi-response variable of coefficient of friction (µ) and specific wear rate (Ks) of abrasive wear of reinforced PTFE composites using Taguchi Deng and novel hybrid support vector regression (SVR) model.

(a) Major keywords used over the literature on the abrasive tribology of PTFE based composites field (1989–2021), (b) the investigated the research region for abrasive.

The results of the tribological experiments for the various test conditions are shown in Table 1. It was seen that the tribological behaviours of the samples indicated rise and fall trend with varying the parameters. Each trial was performed twice and the average was used for the calculations. Highest SNRs of µ and Ks occurred at 20 and 12 trials, respectively. These give the minimum tribological rate of the filled PTFE composites. In spite of the fact the temperature was not computed, the temperature increased as the sliding distance increased. Data in Table 1 was used for the calibration and validation of the SVR, SVR-PSO and SVR-HHO models.

The results of the µ and Ks are shown in Fig. 2a,b, respectively. It was observed that as the load increase the µ and Ks decrease. The low µ at maximum load is because of the formation tribolayer by the fibres at interacting state as well as temperature and visco-elastic related behaviour. This layer prevented the pin samples to be in direct contact with the abrasive surface. Similar observation found in38 when wear rate and coefficient of friction of plastic reinforced glass fibre was studied against different rough and mild steel surfaces. In the work, the results showed that µand Ks decreased as the load increased from 10 to 15 N. At 6 N µ was high perhaps due to tearing of the fibre tribo layer at the contact region. This finding was contrary to results obtained by39 when abrasive wear of reinforced carbonized bone ash particulate polypropylene was investigated. It was found that Ks increased as the load increased from 5 to 15 N. Ks was low at high load because of great increase in apparent contact area at higher loads thereby leading to increase in contact area permitting a large number of particles to meet the interface and share the stress. This, in turn, reduced the wear rate.

Main effect plot for mean (a) µ and (b) Ks of filled PTFE composites.

As shown in Fig. 2b increase in abrasive sizes decrease the both µ and Ks. High µ at small grit size is related to high roughness of SiC particles that offered significant amount of resistance while the low µ is attributed to smoothness of the SiC particles that offered little resistance to the materials all due to formation of protective layer at contact surface. Decrease in wear rate due to increase in abrasive size is related to clogging of the wear track with wear debris and reduction in cutting efficiency of the abrasives due to transfer. Similar results were reported in40 where the abrasive wear performance of elastomers was investigated. In the study, different abrasive size of 82, 125, 269 and 425 µm was used as the counterface at constant speed of 0.01 m, applied load of 10 N but varying speed of 0.25, 0.5 and 1 m/min. It was observed that as the abrasive size of the SiC increased both the µ as well as Ks reduced.

Figure 2a,b shows the relationship between the parameters and the µ and Ks, respectively. As observed in the figure, increasing the sliding distance increase the µ while a decrease in Ks observed. This is explained on the basis that distance acted as a lubricant to rubbing surfaces and therefore separated pin specimens from the counter front. Ks reduction due to increase in sliding distance is attributed to pull out or fracture of abrasives as a result of presence of tough fibres. Additionally, wear debris is transferred to counterface from the PTFE causing reduced wear rate. This result is validated by41 in which two-body abrasive wear property of silicon carbide filled glass fabric epoxy composites glass was studied. Grit sizes of 600 and 1000 mesh and sliding distance of 25, 50, 75 and 100 m were used as experimental conditions. The results revealed that Ks significantly reduced as abrading distances and grit size increased. Drastic reduction of Ks was observed at 25–50 m.

As the sliding speed increase µ as well as Ks decrease and high µ and Ks are noticed at low speed due to samples’ increased contact time with the counterface (Fig. 2a,b. As the rotation rate increases and the samples are oxidized, the temperature at the contacting surfaces changes. This aids formation of a mechanically mixed and rough coating which is laid on the parts. This coating is impervious to removal and depreciates the µ and Ks significantly. When the surface of the counterface is less touched and the protection of the hard layer is heavily mixed tribological behaviours are reduced. It has been shown that reduction in wear rate depends on production of adherent, uniform and thin hard layer on the counter front when polyphenylene reinforced with CuO, SiC, TiO2 and ZnO nano particles to study the influence of these particles on the wear rate of the composites. It was found that with the addition of these particles there was a formation of adherent and uniform tribolayer between the samples and the counterface especially with 2 wt% of CuO and TiO242. More so, it was observed in43 when the function of filler deformation, filler-bonding and counterface of Ag2S, CuS, ZnF2 and SnS inorganic fillers were introduced into polyphenylene sulfide. The results indicated that minimum wear achieved by Ag2S and CuO was due to the formation of thin, adherent and uniform transfer film formed on the counterface preventing the samples from coming into direct contact with surface while ZnF2 and SnS high wear rate was attributed to thick, non-uniform tribolayer between contact surface and samples. Similar observation was made in44 when SiC and graphite particles as secondary filler was added to epoxy polymer to study the influence of speed, distance and load on the epoxy composites. It was found that as the speed increased the wear properties of the reinforced epoxy composites decreased results. Similar result was reported by44 when 3 body abrasive tribological property of glass fibre reinforced polyester composites was investigated. It was found that as the speed increased the wear rate of the composites decreased at constant applied load as well as particles of 200 and 300 6 µm.

In all the analysis, it was found that addition of bronze and carbon fillers into PTFE improved the virgin PTFE’s tribological rate. This might be attributed to stiffness and hardness of the fillers. However, BF40 composites showed a slightly lower wear resistance than CF25 composites. This is explained on this basis of the higher weight percentage of the bronze particles that induced more hardness and larger size about 6 µm.

The experimental data in Table 1 was transformed into signal to noise ratios (SNRs) using Eq. (4) and the corresponding SNRs was as shown in Table 1. Larger SNRs indicate the minimum variation difference between wanted response and measured response. The maximum value of SNRs at the main effect plot for SNRs give the desired results. Figure 3a,b shows the mean SNRs of µ and Ks, respectively. Table 2a,b presents the computed mean SNRs for the µ and Ks, respectively. As seen in Fig. 3a, the maximum mean SNR achieved for µ were load at 9 N, grit size at 1000 mesh, sliding distance at 25 m and sliding speed at 0.14 ms−1. Thus, the estimated optimum parameters for achieving a minimum µ via Taguchi optimization can be coded as L3G1SD1SS3. For the Ks (Table 2b; Fig. 3b), the highest mean SNR obtained for Ks were load at 9 N, grit size at 1000 mesh, sliding distance at 55 m and sliding speed at 0.04 ms−1. Therefore, by Taguchi method the predicted optimum parameters are styled as L3G1SD3SS1. ANOVA depicts the parametric setting that significantly influence the abrasive behaviours. Similarly, the important parametric factor that significantly affect the µ were found as grit size followed by load, distance and speed 3(a). The percentage contribution of grit size, load, distance and speed were computed as 37.24%, 33.92%, 17.62% and 11.20%, (Table 3a). Table 2b shows the percentage contribution of the parameters on Ks. As seen grit size contributed 51.06%, load contributed 24.65% distance contributed 22.57% and speed contributed 1.72% implying that grit size most significantly influence the Ks followed by load, distance and speed, respectively.

Main effect plot for SNRs of (a) µ and (b) Ks.

As seen above, Taguchi can only optimize a parameter at a time and thus it involves more cost, time and effort. Therefore, Deng widely called grey relational analysis (GRA) is principally used to optimize multiple parameters by combining all outputs into an output. Deng is used to unravel real problems made up of a bounded amount of data. It is commonly employed to approximate the properties of indefinite systems having no black and white solution. With respect to grey system black signifies without information whereas white connotes with information. This technique is largely utilized to maximize or minimize problems which have to do with several factors and responses. Data preprocessing through GRA was executed on test data of responses in Table 1 i.e. µ as well as Ks. Table 3 shows the reference sequence obtained by normalization (Eq. 5). In due course, the deviation sequence was computed following Eq. (6) (Table 3). Grey relation coefficient (GRC) and grey relational grade (GRG) of µ as well as Ks were determined using Eqs. (6) and (7), respectively. Subsequently, the mean of GRCs is calculated to establish the GRG. Calculated values of GRGs were employed to produce equivalent SNRs. A larger magnitude of SNR is useful alluding the tests lay in proximity to the actual normalized magnitude of GRG. Figure 4 depicts the plot of GRG against SNRs. It indicates that the 21st trial possesses the highest SNR. Correspondingly, the first rank was designated to 21st trial. The straggling disposition of the GRG, below the plot of SNRs in Fig. 4, also adds to the aforementioned explanation. Ever the ranks determined (Table 5), GRG response table was contrived. Individual factor of GRG at the preferred level was chosen as well as average computed to obtain the mean GRG for separate parameters. The mean response table for the GRG is presented in Table 6.

For example, variable G at level 1 in the first, fourth and seventh runs of the test. The concomitant GRG values in Table 5 were used for computation using Eq. (8). The mean of the chosen GRGs was computed through the method aforementioned put together to generate the mean response table (Table 4). The grades in the response table is used as a degree of correlation45. Hence, from Table 4, it is possible to achieve combination of optimum parameters which maximize the overall response. As observed in Table 5, the maximum GRG exists at L3, G1, SD3 and SS3. Therefore, to wrap it up, the best parameter settings for useful abrasive tribological behaviours of filled PTFE composites are load at 9 N, grit size at 1000 mesh, distance at 55 m and sliding speed at 0.14 ms−1 coded as L3G1D3S3. ANOVA for GRG shows that grit size with 68.57% ranks as the most influential followed by load with 20.57%, followed by distance having a contribution of 7.78% and finally speed with least contribution of 3.38% for minimum tribological loss. Ramesh and Suresha used Taguchi Deng model to optimize the abrasive wear mode of carbon fabric reinforced epoxy composite filled with Al2O3 and MoS2 as fillers. They reported that optimum parameters for minimum wear rate were found to be load at level 3 (15 N), grit size at level 3 (320), filler content at level 3 (10 wt%) and sliding distance at level 3 (30 m). ANOVA revealed that filler content with 52.08% was the most significant factor affecting the wear mode of the composites46.

Having determined the optimum parameters, the final phase in Taguchi-Deng is prediction as well as validation of performance enhancement of the dual responses. The predicted GRG was calculated as per Eq. (7). Validation experiments were executed to validate the results of the analysis. Validated results showed that minimum µ and Ks were 2.0 × 10–1 and 1.5353 × 10–6 mm3 N−1 m−1, respectively. More so, it can be implied from Table 6 that the findings of the validation phase are consistent with the computed values. Besides, an enhancement of 55% in GRG was achieved (Table 6). This performance enhancement in the results obtained through the experiments over the initial design parameter confirms the validity of the Taguchi-Deng method for studying the abrasive tribological behaviours of filled PTFE composites. An improvement of 8.4% in GRG was reported by47 when the same method was used to optimize wear parameters of silicon nitride reinforced AA6063 matrix composites.

One of the aims of this work is building hybrid SVR models namely SVR-PSO and SVR-HHO models and compare their efficiency in predicting tribological behaviours of filled PTFE composites. For this objective, tribological behaviours (µ and Ks) were obtained via experimental results of Table 1. Prediction of tribological behaviours by traditional methods is time and energy consuming due to nonlinearity between tribological independents and dependents of filled polymer composites leading to inaccuracy. These issues can be addressed by nonlinear models. Subsequently, this section details the results achieved in visualized and graphical forms. Before the models simulations, the data was normalized using Eq. (16). Data normalization disallows larger values overshadowing lower values, takes care of units and improves the efficiency of the models.

The simulation process was performed in MATLAB 9.3 (R2020 (a)).Optimized structure of SVR model was chosen via trial-and-error approach. An efficient model is that which allays the prerequisites of model appraisal metrics. Prediction efficiency of the models were appraised using two goodness of fit (R2, R) and two prediction error (RMSE, MAPE) metrics in training as well as testing regimes. The simulated outcomes of the individual SVR models for the prediction of µ and Ks are quantitatively presented in Table 7. From Table 7, it can be seen that the single SVR models achieved various adequacies according to the statistical evaluation metrics. More so, SVRµ shows best results in terms of goodness of fit in both testing and training stages as compared to SVRKs model. However, with respect to prediction errors SVRKs with RMSE 5 × 10–6 and MAPE 29% proved to be a relatively adequate model in predicting the tribological behaviours of filled PTFE composites than SVRµ whose accuracy is extremely poor (61%). To have a graphical map of SVR models for the tribological behaviours, a scatter plot is used. A scatter plot gives the degree of agreement between measured and calculated values for the overall goodness of fit. Figure 5a,b depicts the scatter plot of the whole data for SVRµ as well as SVRKs models, respectively. Arising from the scatter plots, it is interesting to note here that SVRKs model indicated better fitness in comparison to SVRKs when the whole data points were put together.

Scatter plot of (a) SVRµ model; (b) SVRKs model for the whole data.µ

Nevertheless, overall prediction accuracy of the single SVR models was inadequate, especially for SVRµ model. The accuracy can be enhanced using optimization approaches namely PSO and HHO. Fundamentally, it should be considered that the promising prediction accuracy occurred in the course of the training state which is originally used to measure precisely the model based on known inputs and outputs. Nevertheless, the verification stage is significant in appraising the prediction efficiency of the models since it inspects closely the models’ prediction accuracy based upon unknown magnitudes. This advantage is not enjoyed by the training phase. Consequently, a robust model should possess determinate and balanced performance in both training and testing regimes. In general, hybridized models showed a promising ability when compared to un-hybridized models. For consistency the same model evaluation metrics are used to assess the prediction accuracy of the hybridized models. Table 9 shows the results of the hybrid models in both calibration and validation regimes. In spite of the fact it is hard to rank the models as per the model evaluation criteria, the SVR-HHO model indicated higher prediction accuracy in both conditions. From Table 8, it was observed that SVR-HHOµs indicated R2 > 90% R = 95%, 99.26%, RMSE > 5%, and MAPE of 5% = Similarly, SVR-HHOKs R2 > 95%, R > 97%, RMSE < 1% as well as MAPE = 3%. This implies SVR-HHO model performed better than SVR-PSO model for prediction of the tribological behaviours of the filled PTFE composites. The predictive superiority of HHO to others is in concord with results obtained by48. Figures 6 and 7 present the scatter plot of the SVR hybrid models. Close consistency between measured and calculated points was achieved for SVR-HHO model as compared to SVR-PSO model. More so, R values of the hybrid models lie between 85 and 99%. This agrees with conclusions drawn by49,50,51 that values of R greater than 70% are regarded as acceptable. Therefore, all the optimized hybrid models are acceptable (Table 8).

Scatter plot for (a) SVR-PSOµ and (b) SVR-HHOKs models all data sets.

Scatter plot for (a) SVR-PSOµ and (b) SVR-HHOKs models for all dataset.

Zheng et al. used SVR coupled with guided PSO to predict the wear rate of aeroengine. The results indicated better prediction accuracy when compared with single SVR52. In a related development, Kahhal et al. optimized and predicted the process parameters of friction stir welding of AH12 1050 using response surface algorithm coupled with PSO model. It was reported that that the hybrid model indicated superior prediction accuracy53. The prediction efficiency of HHO coupled with generalized neural network was reported in the literature54 to predict the abrasion resistance of ultra-strength martensitic steel. The results showed that observed and the computed values of the wear property lie within the uncertainty range of 3–4%. In a dissimilar fashion, Sammen et al. hybridized artificial neural network (ANN) with PSO, HHO and genetic algorithm (GA) to predict the scour depth downstream of ski-jump spillway. It was observed that of all the hybrid models, ANN-HHO model superior prediction accuracy with mean absolute error of 0.1760 and root mean square error of 0.253854. The findings in these literatures reinforced the results of this study that hybrid models could increase the effectiveness of single models.

The SVR model and its hybrids namely SVR-PSO and SVR-HHO models are compared via 2D Taylor’s plot as shown in Figs. 8 and 9, respectively. As seen in the Taylor’s plot SVR-HHO model indicated better fitness in both cases with values of 97% and 99% for µ and Ks, respectively in the calibration regime. Therefore, it can wrapped up that SVR, SVR-PSO and SVR-HHO models can understand and follow the intricate and nonlinear correlation between tribological input parameters and response parameters of filled PTFE composites in abrasive conditions. Additional analysis can be done using a radar plot for the prediction of the µ and Ks as shown in Fig. 10. It can as well be seen that SVR-HHOµ > SVR-PSOµ > SVRµ and SVR-HHOKs > SVR-PSOKs > SVRKs. This implies that in both cases SVR-HHO model was able of capturing the best fitting trend of the tribological behaviours of filled PTFE composites.

Taylor diagram for µ models in (a) calibration stage and (b) verification stages.

Taylor diagram for Ks models in (a) calibration and (b) verification of the models.

Radar plot for (a) µ and (b) Ks in calibration and verification regimes.

In spite of the availability of ample of predictive models, there is no particular model that can guarantee consistent optimal performance in addressing various types of issues. Yet, latest research on unique population-based as well as nature inspired optimization paradigm (HHO) models indicated better adequacy in establishing optimal solutions for multi-objective problems. The statistical analysis outcomes and comparisons revealed that SVR-HHO model yields promising and often competitive outcomes than well-established models.

Multi linear regression MLR has been widely used as conventional approach for prediction tribological behaviour of materials. Ikpambese and Lawrence used MLR and artificial neural network (ANN) to predict friction and tribological rate of reinforced palm kernel epoxy for brake pad application. The prediction accuracy in terms of RMSE for the MLR and ANN models were found to be 0.0082 and 0.00450, respectively55. These values are less effective compared to RMSE of 0.000001 for SVR-PSO and SVR-HHO in this article. Similarly, Altay et al. used LR, SVM and Gaussian process regression to predict the wear rate of ferro-alloy coating56. They reported RMSE and R2 for LR, SVM and GPR as 0.86, 0.69, 0.69 and 0.93, 0.96 and 0.96, respectively which were higher than the reported RMSE and R2 of this article was higher than as depicted in Table 8. With this it can be said that the proposed models of this study are more effective than the conventional models in predicting the wear rate of materials.

Using the proposed models in the scientific community and industry for the optimization and prediction of composite materials will lead to many benefits. Among these benefits are significant aspects such as minimization of cost, reduction in human effort, prevention of time loss during experiments. The developed models are expected to be specifically applied in the design and/or modification and development of existing/new materials for tribological applications in agricultural machinery, rolling industry, and mining industry, such as crushing and milling.

In this study, multi-response optimization and prediction of objective functions of coefficient of friction (µ) and specific wear rate (Ks) of tribological behaviour of filled polytetrafluoroethylene (PTFE) composites was presented using Taguchi Deng and two hybrid SVR models (SVR-PSO and SVR-HHO). More so, effect of load, grit size, distance and speed on the two responses were studied. It was found that as load, grit size and speed increased µ and Ks decreased. However, increase in distance led to increase in µ and decrease in Ks. Based on the Taguchi Deng optimization, the optimum parameter levels were found to be load at level 3 (9 N), grit size at level 1 (1000 mesh), distance at level 3 (55 m) and speed at level 3 (0.14 ms−1) coded as L3G1SD3SS3. ANOVA for GRGs indicated that grit size with 68.57% was the most influential parameter, followed by load with 20.57%, followed by distance with 7.78% and then speed with least influence of 3.38% affecting the tribological behaviours of filled PTFE composites. Validation test revealed that there was an enhancement of 55% in GRG from 0.4335(L1G3SD3SS1) for initial design settings to 0.9589 for the optimized levels (L3G1SD3SS3). For the prediction accuracy of the models, it was found that SVRKs outperformed SVRµ model. As regards to hybrid models, there was an increase in prediction effectiveness for both SVR-PSO and SVR-HHO over SVR. Even though both SVR-HHO and SVR-PSO models were able to accurately predict the µ and Ks, SVR-HHO model exhibited the lowest prediction error of 4.06% on the average as compared to SVR-PSO model whose prediction accuracy was found to be 10.57% on the average. The integration of Taguchi with Deng approach and SVR with PSO and HHO resulted in the optimization and prediction of tribological behaviour of reinforced PTFE composites with low experimental cost and superior accuracy.

The materials utilized in this work are polytetrafluoroethylene (PTFE), carbon-filled composites (CF25) and bronze-filled composites (BF40) because of their availability and wider applications. Abrasive tribological experiment was conducted according to ASTM G99 standard using pin-on-disc tribometer (Model: Arton Paar, Made in Switzerland) shown in Fig. 11. The counterface material for the wear test is a steel of disc 140 mm in diameter and thickness of 10 mm that has been heat treated to obtain a surface hardness of 55–60 RC. This is grounded to a surface finish of nearly 0.12 µm centerline average. The samples of dimension 20 mm long, width 10 mm and depth 6 mm were cut from compression moulded rectangular plates whose dimensions are (500 × 500 × 6) mm using computer numerical control water jet machining for the pin-on-disc abrasive experiment. A specially designed fixture for holding the square samples was designed and fabricated. The samples were inserted into the fixture, bolted and then loaded against silicon carbide (SiC) abrasive papers glued to the hardened steel holder by means of liquid adhesive. Control parameters and their levels are shown in Table 9.The experimental design is as shown in Table 10. In all the experiments, mass before (m1) and mass after (m2) was measured using digital weighing balance (Model: PS 1000.RS RADWAG, Made in Poland) with 10–3 g precision accuracy. Test was performed at room temperature (29 °C and relative humidity 55%). Samples were cleaned with a brush before and after the experiment to get rid of debris and then weighed. The loss in pin weight (WL), volume (VL) and specific wear rake (Ks) was determined through the mathematical Eqs. (1), (2) and (3), respectively.

where \({\mathrm{W}}_{\mathrm{L}}\) = weight loss (g), \({\mathrm{V}}_{\mathrm{L}}\) = volume loss (mm3) \({m}_{b}\) = mass before test (g), \({m}_{a}\) = mass after test (g), \(\uprho\) = (g cm−3) of materials, L = load in N and D = sliding distance (m). Each trial was performed twice and averaged.

Arton Paar Tribometer used for the experiment.

Single parameter optimization and percentage contribution of each parameter can be performed using Taguchi L27 (34) orthogonal array (OA). Taguchi is an optimization process to establish the best process parameter. In this study, four parameters with three levels (34) are considered for the configuration of the test. Twenty seven trials have been performed based on Taguchi L27 (34) OA as shown in Table 10. The Taguchi method was chosen in this study because of its simplicity of analysis, substantial reduction in the cost of experiment and validity over a wide region spanned by the control factors and their settings. Signal to noise ratios (SNRs) are used to establish the optimum parameters. There exist three types of quality loss function for SNRs namely: nominal the better; higher is the better; lower the better. Thus, the lower the better function was used to obtain the optimum wear-rate parameters in this study. SNRs for µ and Ks was computed using Eq. (4) for all the 27 trials in accordance with Taguchi \({L}_{27}\) (34) OA.

where n = number of experiments and \({\mathrm{y}}_{i}\) = experimental value.

Taguchi optimization is capable only of optimizing a single response. However, when two or more responses of distinct features are involved Taguchi technique is limited. Thus an optimization method called Deng popularly referred to as grey relational analysis (GRA) becomes a panacea. This approach is capable of addressing the limitation imposed by Taguchi approach by integrating multiple responses into a single response whose optimal parameters combination represents the several responses thereby minimizing cost and time spent while using Taguchi optimization method. Taguchi L27 (43) OA with Deng was used to obtain the optimum levels of tribological parameters. Data normalization is categorized as smaller or larger the better. Let the actual sequence and the comparison sequences be \({X}_{i}^{*}\left(k\right)\) and \({\varphi }_{i}\left(k\right)\) , respectively. \(\mathrm{i}\) = 1, 2, 3….; m = 1, 2, 3…and n and m represent the total number of experiments and experimental values, respectively. Data preprocessing is used to transform the actual sequence into an identical sequence. Many data preprocessing techniques can be utilized in Taguchi-Deng method, depending upon the features of the actual sequence. Generally, series is normalized between 0 and 120. For this study, the target value is “the smaller the better”. Consequently, the actual sequence is pre-processed via Eq. (5).

where \({X}_{i}^{*}\left(k\right)\) = normalized for the ith experiment and \({\varphi }_{i}\left(k\right)\) = initial sequence of the average responses. After data normalization, the succeeding phase is computation of deviation sequence of the normalized data using Eq. (6).

where \({\Delta }_{oi}\left(k\right)\) = deviation, \({X}_{0}^{*}\left(k\right)\) = normalized data and \({X}_{i}^{*}\left(k\right)\) = comparability sequence. Grey relational coefficient (GRC) is thus estimated through Eq. (7).

where \({\xi }_{i}\left(k\right)\) = GRC of each response, \({\Delta }_{min}\) and \({\Delta }_{max}\) = lowest and the highest deviations of the individual target factor, respectively. Differentiating or identification coefficient is symbolized by \(\zeta\) and is demarcated within the range of \(\zeta \epsilon \left[\mathrm{0,1}\right]\) . This is usually set at ½ to assign equivalent weights to every variable. As indicated in (Eq. 8) GRG is then determined by taking mean of GRG of each output parameter:

where \({\gamma }_{i}\) = GRG obtained for ith test run, n = summation count of performance attributes. Following the determination of the optimal levels of parameter, the last phase is to predict and validate the result using Eq. (9):

where \({\gamma }_{0}\) represents the highest of mean GRG at optimum levels of variables and \({\gamma }_{m}\) defines the average GRG. \(q\) = parameter that signifies factors influencing the target values.

In 1995 Vapnik contrived and implemented support vector machine (SVM) was contrived and implemented, which is regarded an observer-based learning approach. The minimization of structural risk as well as statistical learning theory are the most important function of the SVM. Nevertheless, the properties which distinguish SVM from ANN are complexities, minimization of error as well as gain in the network’s performance capability. SVM can be categorized into linear support regression as well as nonlinear support regression (NSVR). Several engineering fields such as have witnessed the application of SVM’s kernel function. SVR model could be thought of as SVM on the basis of layers which include kernel function weighting on the inputs as well as function weighted sum of kernel targets. By and large, SVM is codified into two codes namely Support Vector Regression (SVR) and Support Vector Classifier (SVC) models. SVR model is made up of predictions whereas SVC model treats classifications. SVR model is designated as:

where w stands for weight of the vector displayed in feature space, \(\Phi\) shows the transfer function, b is bias. Therefore, in order to show the SVR function \(f\left(x\right)\) , problem of regression is presented as:

where \(\parallel w{\parallel }^{2}\) = weight norm vector, C = penalty parameter, \({\xi }_{i}\) and \({\xi }^{*}\) = slack variables. By using Lagrange functions, the solution of the nonlinear regression function can be presented based on optimization as follows:

where \(K\left(x,{x}_{i}\right)\) shows the kernel function and are binary variables (\({\alpha }_{i }and {\alpha }_{i}^{*}>0)\) . There exist several kinds of kernel functions including sigmoid, linear, polynomial but the commonly used kernel function is the radial basis function (RBF). Consequently, the RBF kernel was used in this study and it is expressed as (Eq. 13).

where \(\gamma\) = kernel parameter. SVR model performance is affected by C, \(\gamma\) and \(\varepsilon\) (size).

HHO is a unique model worked out by simulating the hawk’s hunting process. Lately, the procedure has been used with success in solving several intricate engineering as well as science issues. The hawks mostly operate alone whereas the Harris hawks pursue and hunt through operating and cooperating together. Hence, the HHO method is similar to Harris Hawks’ natural hunting characteristic and cooperative methodology. HHO model hunting methodology entails tracing, encircling, approaching and attacking. These mechanisms are achieved in three principal phases namely: exploration, a transition from exploration to exploitation as well as exploitation (Fig. 12).

In 1995, PSO was presented by Kennedy and Eberhart. It is a search method based on population which is inspired by the social behaviour and dynamics of animals. The initial intention of the SPO philosophy was to clearly mimic animals’ social behaviour bird flocking as an instance to detect trends that control capability of birds to fly with precision at the same time and to all of a sudden alter the direction with regathering in an optimum style. Arising from this first purpose, the philosophy inspired into a simple and efficient optimization approach. PSO is initiated with a group of random particles that look into an optimum value by updating the two best values in each iteration. The first one is named the personal best (pbest). This is the best value so far obtained by any particle in the population. All the particles explore the search space and the information collected by them is utilized for finding the best particle in the swarm referred to as global best (gbest). Thereafter, the particle updates its velocity and positions according to Eqs. (14, 15):

where \({V}_{i}^{k+1}\) = the velocity of individual I at iteration \(k\) +1, \({V}_{i}^{k}\) = the velocity of individual i at iteration \(k\) , \(\omega\) stands for inertia weight parameter, \({c}_{1}\) and \({c}_{2}\) show the cognitive parameters, \({r}_{1}\) and \({r}_{2}\) = random numbers between 0 and 1, \({X}_{i}^{k}\) = position of individual i at iteration k, \({pbest}_{i}^{k}\) = the best position of individual I at iteration k and \({gbest}^{k}\) indicates the best position of the group until iteration k. Figure 13 shows the flowchart of the PSO algorithm.

Improving the SVR model’s performance needs a cautious delineation of parameters involved in the SVR model. The strength of the SVR model relies upon the precise choice of C, \(\gamma\) and \(\varepsilon\) . Yet, these parameters having a wide range make the search space very large thus making it difficult to choose precise parameters. Therefore, this issue can be addressed as optimization issue that requires sorting out via optimization methods. Integration of SVR model with PSO as well as HHO models that are algorithms inspired by nature led to the following hybrid model namely: SVR-PSO and SVR-HHO for the prediction of tribological behaviours of filled PTFE composites. The nature inspired models were utilized to choose the SVR model parameters viz: C, \(\gamma\) and \(\varepsilon\) . Proposed flow-chart of the hybrid model illustrated in Fig. 14.

Flow chart for the proposed hybrid models.

One of the significant aims of any soft computing model is to ascertain that models conform to acceptable data based on the models evaluation metrics utilized to obtain a dependable and strong computed outcome of the unknown data. Nevertheless, overfitting as well as local minima problems occur in the data validation. Hence, the performance of the learning phase might be unsatisfactory. This is especially when the analysis deals with a relatively small amount of dataset, as in this study. Various validation methods can be employed including cross-validation (k-fold), hide-out and leave one out. Here, the k-fold approach was used to repeal overfitting issues. With respect to this study, the data was split into (70%) and (30%) for training and testing, respectively. The data obtained through abrasive experiments was pre-processed and normalized according to Eq. (16). Data normalization was performed prior to model training and it usually enhances the efficiency of the predictive models. The current work introduced SVR model coupled with particle swarm optimization (PSO) and Harris Hawk optimization (HHO) models to predict abrasive tribological behaviours of filled PTFE composites. Prediction of tribological behaviours is important. However, creation of a reliable model is often challenging and difficult given the nature of the data set obtained from the experiments.

where y = normalized data, \(x\) = is the experimental data while \({x}_{max} and {x}_{min}\) are the maximum and minimum experimental data, respectively.

Generally, model efficiency performance should include at least one goodness of fit and at least one prediction error metrics58. Based on this determination coefficient (R2), correlation coefficient (R), root mean square error (RMSE) and mean absolute percentage error (MAPE) are chosen as models appraisal metrics of the soft computing methods. R2, R, RMSE and MAPE are given below. These statistical tools furnish the information on efficiency of models.

where \(x,y,\widehat{x}\) and \(\widehat{y}\) are the actual, predicted, average actual and average predicted values, respectively.

Unal, H., Sen, U. & Mimaroglu, A. Abrasive wear behaviour of polymeric materials. Mater. Des. 26, 705–710 (2005).

Bijwe, J., Logani, C. M. & Tewari, U. S. Influence of fillers and fibre reinforcement on abrasive wear resistance of some polymeric composites. Wear 138, 77–92 (1990).

Khedkar, J., Negulescu, I. & Meletis, E. I. Sliding wear behavior of PTFE composites. Wear 252, 361–369 (2002).

Suresha, B. & Kumar, K. N. S. Investigations on mechanical and two-body abrasive wear behaviour of glass/carbon fabric reinforced vinyl ester composites. Mater. Des. 30, 2056–2060 (2009).

Harsha, A. P. & Tewari, U. S. Abrsasive wear resitance of glass fibre reinforced polysulfone composites. Indian J. Eng. Mater. Sci. 9, 203–208 (2002).

Kukureka, S. N., Hooke, C. J., Rao, M., Liao, P. & Chen, Y. K. The effect of fibre reinforcement on the friction and wear of polyamide 66 under dry rolling–sliding contact. Tribol. Int. 32, 107–116 (1999).

Ulas, M., Altay, O., Gurgenc, T. & Özel, C. A new approach for prediction of the wear loss of PTA surface coatings using artificial neural network and basic, kernel-based, and weighted extreme learning machine. Friction 20, 1–14 (2019).

Unal, H., Yetgin, S. H., Mimaroglu, A. & Sumer, M. The effect of test parameters on friction and wear performance of PTFE and PTFE composites. J. Reinf. Plast. Compos. 29, 1978–1986 (2010).

Şahin, Y. & Mirzayev, H. Wear characteristics of polymer-based composites. Mech. Compos. Mater. 51, 543–554 (2015).

Article  ADS  CAS  Google Scholar 

He, R., Chang, Q., Huang, X. & Bo, J. Improved mechanical properties of carbon fiber reinforced PTFE composites by growing graphene oxide on carbon fiber surface. Compos. Interfaces 25, 995–1004 (2018).

Suh, J. & Bae, D. Mechanical properties of polytetrafluoroethylene composites reinforced with graphene nanoplatelets by solid-state processing. Compos. Part B Eng. 95, 317–323 (2016).

Shipway, P. H. & Ngao, N. K. Microscale abrasive wear of polymeric materials. Wear 255, 25 (2003).

Ravi Kumar, B. N., Suresha, B. & Venkataramareddy, M. Effect of particulate fillers on mechanical and abrasive wear behaviour of polyamide 66/polypropylene nanocomposites. Mater. Des. 30, 3852–3858 (2009).

Liu, C., Ren, L. Q. & Tong, J. Abrsaive wear behavior of particle reinforced ultrahigh molecular weight polethylene composites. Wear 225–229, 199–204 (1999).

Yousif, B. F., Nirmal, U. & Wong, K. J. Three-body abrasion on wear and frictional performance of treated betelnut fibre reinforced epoxy (T-BFRE) composite. Mater. Des. 31, 4514–4521 (2010).

Pervez, M. N., Shafiq, F., Sarwar, Z., Jilani, M. M. & Cai, Y. Multi-response optimization of resin finishing by using a taguchi-based grey relational analysis. Materials (Basel) 11, 426 (2018).

Article  ADS  CAS  Google Scholar 

Julong, D. Introduction to grey system theory. J. Grey Syst. 1989, 1 (1989).

Dharmalingam, S., Subramanian, R. & Kok, M. Optimization of abrasive wear performance in aluminium hybrid metal matrix composites using Taguchi-grey relational analysis. J. Eng. Tribol. 227, 749–760 (2013).

Sylajakumari, P. A. & Ramakrishnasamy, R. Taguchi grey relational analysis for multi-response. Materials (Basel) 11, 1–17 (2018).

Saravanan, K. G. & Rajasekaran, T. Optimisation of laser parameters and dimple geometry using PCA-coupled GRG. Strojniški Vestn. J. Mech. Eng. 20, 20 (2021).

Pervez, H., Mozumder, M. & Mourad, A.-H. Optimization of injection molding parameters for HDPE/TiO2 nanocomposites fabrication with multiple performance characteristics using the Taguchi method and grey relational analysis. Materials (Basel) 9, 710 (2016).

Article  ADS  CAS  Google Scholar 

Adediran, A. A., Akinwande, A. A., Balogun, O. A. & Bello, O. S. Mechanical and optimization studies of polypropylene hybrid biocomposites. Sci. Rep. 20, 1–17. (2022).

Lal, S., Kumar, S., Khan, Z. A. & Siddiquee, A. N. Multi-response optimization of wire electrical discharge machining process parameters for Al7075/Al2O3/SiC hybrid composite using Taguchi-based grey relational analysis. Proc. IMechE Part B J. Eng. Manuf. 229, 229–237 (2015).

Raykar, S. J., Addona, D. M. D. & Mane, A. M. Multi-objective optimization of high speed turning of Al 7075 using grey relational analysis. Proced. CIRP 33, 293–298 (2015).

Kuram, E. & Ozcelik, B. Multi-objective optimization using Taguchi based grey relational analysis for micro-milling of Al 7075 material with ball nose end mill. Measurement 46, 1849–1864 (2013).

Dhande, D. Y. Comparative analysis of abrasive wear using response surface method and artificial neural network. J. Inst. Eng. Ser. D (2021).

Zhang, Z., Friedrich, K. & Velten, K. Prediction on tribological properties of short fibre composites using artificial neural networks. Wear 252, 668–675 (2002).

Velten, K., Reinicke, R. & Friedrich, K. Wear volume prediction with artificial neural networks. Tribol. Int. 33, 731–736 (2000).

Jiang, Z., Zhang, Z. & Friedrich, K. Prediction on wear properties of polymer composites with artificial neural networks. Compos. Sci. Technol. 67, 168–176 (2007).

Jiang, Z., Gyurova, L., Zhang, Z., Friedrich, K. & Schlarb, A. K. Neural network based prediction on mechanical and wear properties of short fibers reinforced polyamide composites. Mater. Des. 29, 628–637 (2008).

Tang, J., Huang, Z., Zhu, Y. & Zhu, J. Neural network compensation control of magnetic levitation ball position based on fuzzy inference. Sci. Rep. 20, 1–17. (2022).

Haghighat Mesbahi, A., Semnani, D. & NouriKhorasani, S. Performance prediction of a specific wear rate in epoxy nanocomposites with various composition content of polytetrafluoroethylene (PTFE), graphite, short carbon fibers (CF) and nano-TiO2 using adaptive neuro-fuzzy inference system (ANFIS). Compos. Part B Eng. 43, 549–558 (2012).

Ray, S. A comparative analysis of the abrasion wear characteristics of industrial wastes filled glass / polyester composites based on the design of experiment and neural network. Polym. Compos. 20, 1–15. (2020).

Dutta, S., Pal, S. K. & Sen, R. On-machine tool prediction of flank wear from machined surface images using texture analyses and support vector regression. Precis. Eng. 43, 34–42 (2016).

Bhattacharyya, P. & Sanadhya, S. K. Support vector regression based tool wear assessment in face milling. Proc. IEEE Int. Conf. Ind. Technol. 20, 2468–2473. (2006).

Kalantary, S., Jahani, A., Pourbabaki, R. & Beigzadeh, Z. Application of ANN modeling techniques in the prediction of the diameter of PCL/gelatin nanofibers in environmental and medical studies. RSC Adv. 9, 24858–24874 (2019).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Pachuri, V. S., Sunjay, S., Om Prakash, V., Bhuvnesh, B. & Tarun Kumar, S. Nikhil Prediction and optimization of abrasive wear loss of ultrahigh strength martensitic steel using response surface methodology, Harris Hawk and artificail neural network. Int. J. Syst. Assu Eng. Manage. 20, 1–16 (2021).

Chowdhury, M. A. et al. Experimental investigation of friction coefficient and wear rate of composite materials sliding against smooth and rough mild steel counterfaces. Tribol. Ind. 35, 286–292 (2013).

Asuke, F., Abdulwahab, M., Aigbodion, V. S., Fayomi, O. S. I. & Aponbiede, O. Effect of load on the wear behaviour of polypropylene/carbonized bone ash particulate composite. Egypt. J. Basic Appl. Sci. 1, 67–70 (2019).

Hakami, F., Pramanik, A., Basak, A. K. & Ridgway, N. Effect of abrasive particle size on tribological behavior of elastomers. J. Eng. 1, 1–13 (2019).

Mohan, N., Natarajan, S. & Kumareshbabu, S. P. Investigation on two-body abrasive wear behavior of silicon carbide filled glass fabric-epoxy composites. J. Miner. Mater. Charact. Eng. 9, 231–246 (2010).

Bahadur, S. & Sunkara, C. Effect of transfer film structure, composition and bonding on the tribological behavior of polyphenylene sulfide filled with nano particles of TiO2, ZnO, CuO and SiC. Wear 258, 1411–1421 (2005).

Schwartz, C. J. & Bahadur, S. The role of filler deformability, filler–polymer bonding, and counterface material on the tribological behavior of polyphenylene sulfide (PPS). Wear 251, 1532–1540 (2001).

Basavarajappa, S., Arun, K. V. & Davim, J. P. Effect of filler materials on dry sliding wear behavior of polymer matrix composites—a Taguchi Approach. J. Miner. Mater. Charact. Eng. 8, 379–391 (2009).

Wojceichowski, S., Maruda, R. W., Krolczyk, G. M. & Nieslony, P. Application of signal noise ratio and grey relational analysis to minimze forces and vibrations during precise ball end milling. Precis. Eng. 51, 582–5596 (2018).

Ramesh, B. N. & Suresha, B. Optimization of tribological parameters in abrasive wear mode of carbon-epoxy hybrid composites. Mater. Des. Des. 59, 38–49 (2014).

Stalin, B., Kumar, P. R., Ravichandran, M., Kumar, M. S. & Meignanamoorthy, M. Optimization of wear parameters using Taguchi grey relational analysis and ANN-TLBO algorithm for silicon nitride filled AA6063 matrix composites. Mater. Res. Express 6, 106590 (2019).

Article  ADS  CAS  Google Scholar 

Yu, J., Kim, C.-H. & Rhee, S.-B. The comparison of lately proposed Harris Hawks optimization and Jaya optimization in solving directional overcurrent relays coordination problem. Complexity 20, 1–22 (2020).

LeGates, D. R. & McCabe, G. J. J. Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 35, 233–241 (1999).

Moriasi, D. N. et al. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 50, 885–900 (2007).

Musa, B., Yimen, N., Abba, S. I. & Adun, H. H. Multi-state load demand forecasting using hybridized support vector regression integrated with optimal design of off-grid energy systems—a metaheuristic approach. Proceses 20, 1166 (2021).

Zheng, B., Gao, F., Ma, X. & Zhang, X. Intelligent prediction of aeroengine wear based on the SVR optimized by GMPSO. Appl. Sci. 11, 10592 (2021).

Kahhal, P., Ghasemi, M., Kashfi, M. & Menghari, H. G. A multi-objective optimization using response surface model coupled with particle swarm algorithm on FSW process parameters. Sci. Rep. 12, 1–20 (2022).

Article  ADS  CAS  Google Scholar 

Sharma, V., Sharma, S., Prakash, O. & Bhuvnesh, V. Prediction and optimization of abrasive wear loss of ultrahigh strength martensitic steel using response surface methodology, Harris Hawk and artificial neural network. Int. J. Syst. Assur. Eng. Manage. (2021).

Ikpambese, K. K. & Lawrence, E. A. Comparative analysis of multiple linear regression and artificial neural network for predicting friction and wear of automotivebrake pads produced from palm kernel shell. Tribol. Ind. 40, 25 (2018).

Altay, O., Gurgenc, T., Ulas, M. & Özel, C. Prediction of wear loss quantities of ferro-alloy coating using different machine learning algorithms. Friction 8, 107–114 (2020).

Islam, M. Z. et al. A Harris Hawks optimization based single- and multi-objective optimal power flow considering environmental emission. Sustainability 12, 5248 (2020).

Abba, S. I., et al.. Modelling of Uncertain System : A comparison study of Linear and Non-Linear Approaches. 1–6 (2020).

Mechanical of Engineering Department, Faculty of Engineering, Kano University of Science and Technology, Wudil KM 50, Kano, Gaya Road, Wudil P.M.B 3244, Kano, Kano, Nigeria

Mechanical Engineering Department, Faculty of Engineering, Near East University, via Mersin 10, 99138, Nicosia, Turkey

Musa Alhaji Ibrahim, Hüseyin Çamur & Mahmut A. Savaş

Interdisciplinary Research Center for Membrane and Water Security, King Fahd University of Petroleum and Minerals, Dhahran, 31261, Saudi Arabia

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

M.A.I. designed, conducted and analyzed the experiments and the results. S.I.A. performed the simulation. H.C. and M.A.S. reviewed the manuscript.

Correspondence to Musa Alhaji Ibrahim.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Ibrahim, M.A., Çamur, H., Savaş, M.A. et al. Optimization and prediction of tribological behaviour of filled polytetrafluoroethylene composites using Taguchi Deng and hybrid support vector regression models. Sci Rep 12, 10393 (2022).


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.