Next Article in Journal
Quarter of a Century after: A Glimpse at the Conformation and Mechanism of Candida antarctica Lipase B
Next Article in Special Issue
Hot Deformation Behavior and Microstructure Characterization of an Al-Cu-Li-Mg-Ag Alloy
Previous Article in Journal
Investigation of Total Phenolic Content and Antioxidant Activities of Spruce Bark Extracts Isolated by Deep Eutectic Solvents
Previous Article in Special Issue
Simulation and Prediction of the Vickers Hardness of AZ91 Magnesium Alloy Using Artificial Neural Network Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Modeling and Analysis of Tribological Properties of AA6063 Aluminum Alloy Reinforced with Fly Ash by Using Response Surface Methodology

by
Alaa Mohammed Razzaq
1,2,
Dayang Laila Majid
3,*,
Mohamad Ridzwan Ishak
3 and
Uday Muwafaq Basheer
4,5
1
Oil Products Distribution Company, Iraq Ministry of Oil, Baghdad 10022, Iraq
2
Department of Materials Technology Engineering, Technical Engineering College-Baghdad, Middle Technical University, Baghdad 10074, Iraq
3
Department of Aerospace Engineering, Faculty of Engineering, Universiti Putra Malaysia, Serdang 43400, Malaysia
4
Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, Johor Bahru 81310, Malaysia
5
Centre for Advanced Composite Materials (CACM), Institute for Vehicle Systems and Engineering, Universiti Teknologi Malaysia, Skudai 81310, Malaysia
*
Author to whom correspondence should be addressed.
Crystals 2020, 10(5), 403; https://doi.org/10.3390/cryst10050403
Submission received: 12 February 2020 / Revised: 17 March 2020 / Accepted: 1 April 2020 / Published: 16 May 2020
(This article belongs to the Special Issue Microstructural and Mechanical Characterization of Alloys)

Abstract

:
Lightweight, high-strength metal matrix composites have attracted considerable interest because of their attractive physical, mechanical and tribological properties. Moreover, they may offer distinct advantages due to good strength and wear resistance. In this research, AA6063 was reinforced with FA particles using compocasting methods. The effects of fly ash content, load, sliding speed and performance tribology of AA6063 –FA composite were evaluated. Dry sliding wear tests were carried out according to experimental design using the pin-on-disc method with three different loads (24.5, 49 and 73.5 N) and three speeds (150, 200 and 250 rpm) at room temperature. Response surface methodology (RSM) was used to analyze the influence of the process parameters on the tribological behavior of the composites. The surface plot showed that the wear rate increased with increasing load, time and sliding velocity. In contrast, the friction coefficient decreased with increasing these parameters. Optimal models for wear rate and friction coefficient showed appropriate results that can be estimated, hence reducing wear testing time and cost.

1. Introduction

Light-weight metal matrix composites (MMCs) have paved the way to replace steels and reduce greenhouse effects, hence they are used in various applications. A particular example is discontinuous particulate reinforced aluminum composites (DPAMCs), employed in automotive, aerospace and defense sectors [1]. DPAMCs using high-strength refractory particle reinforcements such as SiC, TiC, B4C, Al2O3, MgO, TiO2 were found to improve wear properties; attempts were undertaken to model wear behavior effectively. The modeling involves the process of evaluating the relationship between two or more process variables, as well as with the response variable. To be able to do so effectively, requires a mathematical model that is able to predict the response outcome based on various effects of the process variables [2]—especially in the evaluation of material properties. Techniques such as multiple regressions and artificial neural networks were among those applied to build mathematical model. Mechanical and tribological properties can be also predicted by employing other statistical tools, e.g., design of experiments, analysis of variance (ANOVA) and regression analysis [3,4,5]. Typically for wear tests, the process variables are the load, sliding speed, time and the reinforcement content.
Kumar et al. [6] investigated the size range of fly ash FA particles on dry-sliding wear behavior of aluminum composites. They observed that the wear and friction coefficient of the composites were influenced by sliding speed, reinforcement content and applied load. They developed a mathematical model for prediction of wear rate and the friction coefficient of aluminum composites reinforced with fly ash while the adequacy of the model is evaluated using ANOVA. Rajmohan et al. [7] observed that wear rate increased with increased load. They employed analysis variance ANOVA to analyze the wear behavior of the composite. They found that load is a significant factor in the wear performance of the composite. Magibalan et al. [8] reported that AA8011 composite with 12% fly ash produced lower wear rates and increased hardness with a value of 63.37 VHN. They successfully developed a regression model that can predict the wear rate of cast AA8011-12% FA composites up to 95% confidence levels, within the limits of investigation. Normal load, time and sliding velocity were shown to have a strong influence on the wear rate of the composites.
Similarly, Mandaletal et al. [9] predicted the wear behavior of metal matrix composites using regression methods and significance of wear factors using ANOVA method. The increase of parameter was also done by style trained software package. Their work suggested that response surface methodology may be suitable to evaluate the wear performance of combined rolling action and sliding. Moreover, Monikandan et al. [10] presented the tribological properties of mono AA6061-10 wt% B4C and hybrid AA6061-10 wt% B4C-7.5 wt% Gr composites which were proposed as potential substitute for aluminum alloys used in automotive engines. A full factorial experimental design was conducted, and their findings suggested that the wear loss increased with sliding distance and applied load, while the coefficient of friction also increased with increase in applied load. In addition, the ANOVA analysis revealed the statistically and physically significant factors that influence the wear loss and friction coefficient. Processing-wise, AA6061-15%TiC aluminum matrix composites were produced by the stir casting method [11]. An empirical relationship incorporating stir casting parameters was developed for predicting ultimate tensile stress. A lower or higher combination of these parameters resulted in lower ultimate tensile stress. In this work, a compocasting technique that utilized a lower processing temperature compared to stir casting method was employed, where the lower temperature aided in reducing interfacial reactions.
The aim of the current work is to develop models to predict wear rate and coefficient of friction as functions of important wear variables. The influences of control factors such as fly ash content, load and sliding speed on wear rate and friction coefficient of AA6063-FA composites prepared by compocasting technique were experimentally evaluated considering design of experimental DOE- based central composite design (CCD) technique. Full quadratic RSM based predictive models for Wr and Cr were developed and subsequently validated with experimental results.

2. Experimental

2.1. Composites Fabrication

AA6063 alloy was selected as the matrix material. The AA6063-FA composites were fabricated by compocasting technique with wt% FA (0–12%) with the average particle size of FA particles 45 µm. The manufacturing process, metallographic of the composite, mechanical properties and its XRD spectra are available in an earlier works [12,13]. For each sample the experiment was conducted three times and the average of closely repeatable test values were registered as wear rate and coefficient of friction for composites, that it has explained in previous published papers [14,15]. Chemical composition of AA6063 and FA are presented in Table 1 and Table 2 respectively [16].

2.2. Statistical Analysis Using Surface Response Method SRM

The data obtained from the wear tests on the AA6063-FA composites was evaluated in Design Expert V10, Design of Experiments/ statistical analysis software using the Response Surface Method (RSM). Two models were developed, one for the wear rate and the other for the coefficient of friction. The F value is a ratio of the model and residual variance. If the ratio is near to the value of 1, then it means that the term has a remarkable influence on the response. Likewise, a p-value larger than 0.10 indicated that the term has an insignificant effect [17].
The analysis of variance ANOVA technique was employed in analyzing the wear performance of AA6063-FA composites subjected to important controlled parameters. In current study FA was observed as a remarkable contributor in the wear performance of the composites. The robustness of the developed model is typically checked using ANOVA technique [18,19,20]. In this technique, the F ratio of the developed model that was computed must be lower than the standard tabulated value of F ratio for a desired level of confidence (say 95%) in order for the model to be satisfactory within that particular confidence limit. Furthermore, the coefficient of correlation ‘r’ was utilized in finding the closeness of predicted and experimental values which is described by the following Equations described by Kumar and Balasubramanian [18]:
r 2 = explained   variation   total   variation   = Σ ( W r p W r a ) 2 Σ ( W r e W r a ) 2
r 2 = explained   variation   total   variation   = Σ ( F r p F r a ) 2 Σ ( F r e F r a ) 2
where Wrp, Frp are the predicted wear rate and coefficient of friction respectively for the given experimental wear rate, Wre and coefficient of friction Fre, for the corresponding average of experimental wear rate Wra and average coefficient of friction Fra respectively.

2.3. Experimental Design

RSM was used to model the wear rate and friction coefficient of AA6063-FA composite material. The response function representing the variables, namely, FA content, sliding speed and applied load can be expressed as [21,22] Y = f (FA, S, L). The response surface Y for the k parameter is represented in a second order polynomial regression equation, presented as the following:
Y = b 0 + i = 1 k b i x i + i = 1 i   j k b i x i 2 + i   j k b i j x i x j
Expanding the polynomial for three factors will give the equation:
Y = b 0 + b 1 F A + b 2 L + b 3 S + b 11 F A 2 + b 22 L 2 + b 33 S 2 + b 12 F A L + b 13 F A S + b 23 L S
where b0 is the average of responses and b1, b2, b3, b12, b13 and b23 are the coefficients of response dependent upon the respective main and interaction effects of parameters.
From the factorial design, the combined effects of all parameters will result in an empirical relationship between the wear rate and coefficient of friction with the designated parameters which are the fly ash content (FA), load (L) and sliding speed (S). The multiple regression technique used to develop the mathematical model to formulate the influence of the parameters on both responses were achieved via design of experiments analysis software, Design-Expert V10.0.
Feasible limits of all three factors were determined and these are presented in Table 3.
There were 63 experimental response values of various combination of factors obtained and used to formulate the response function (wear rate and coefficient of friction) as shown in Table A1 (Appendix A) and also as input into the Design-Expert V10.0 software. The highlighted results (Appendix A), which represent 20% of the total data, was not used in the model construction process. This data were utilized in verification of the proposed empirical model. The rest of the data were introduced in full factorial design using response surface method with user-defined choice. The software then automatically will evaluate all possible combinations of the factors and produce results for the significant and aliased model terms.
The design summary of the problem is tabulated in Table 4, which shows the summary of the input information and setting required to starting the analysis. In summary, the methods of the current work encompass a quantitative (wear test) and statistical analyses (randomized RSM approach), whereby the statistical method provide a means to assess the effect of all three factors, namely the fly ash content (FA), load (L) and sliding speed (S) on the wear properties with a quadratic design model.

3. Results and Discussion

3.1. Statistical Analyses

In this section, the effect of FA content, sliding speed and applied load on wear rate and coefficient of friction are discussed. RSM was applied to evaluate and develop models of wear rate and friction coefficient of AA6063-FA composites.
In an attempt to obtain a model with better performance than the one given by the quadratic model of the process, a cubic model equation was developed, but owing to the fact that the developed cubic model was found to be aliased in nature, its model expression in terms of the actual factors was not given by the Design Expert. However, the aliased cubic model, the mathematical expression of which was not available, was still analyzed for variance to see if actually increasing order of model had any effect on the response statistically.
The models summary statistics given in Table 5 and Table 6 show that the quadratic models were suggested according to the statistical evaluation of different models. The software also revealed that the cubic order model was aliased for the ranges of the data obtained. The quadratic model was used according to the value of R-squared (0.9791, 0.9818) for wear rate and coefficient of friction model respectively.

3.2. Modeling and Analysis of Variance of Wear Rate

The design of experiment software with response surface method was used to statistically study the effect of fly ash (FA), sliding speed (S) and applied load (L) on the wear rate of AA6063-FA composites and to build an empirical model for Wear rate based on the effects of these factors. The sequential F-test was conducted for testing the significance of the regression model.
The step-wise regression method was utilized for the wear rate model with the entire potential combinations of the control factors with exception of the cubic terms since these terms were aliased for the data ranges obtained. The selected terms along with the step-wise regression method led to the elimination of the less significant model terms automatically. Table 7 shows the analysis of variance generated for wear rate model. The model F-value of 213.395 confirms the significance of the model. A p-value of 0.0001 was indicative that there was a low chance of F-value which could occur as result of noise. The predicted R2 of 0.9658, which represent the measure of the amount of variation about the mean explained by the model agrees well with the adjusted R2 of 0.9745. The adjusted R2 and predicted R2 should be within approximately 0.20 of each other so as to be in agreement, or else setbacks may arise with either the data or the model. High R2 value suggest there is a statistically significant interaction between factors. Resulting wear rate within the investigated ranges of parameters can be modeled by the final equation produced, Equation (5).
W r = 0 . 037685     0 . 10651   ×   FA + 0 . 054471   ×   S + 0 . 019936   ×   L     0 . 001056   ×   FA   ×   S   0.00337143   ×   FA   ×   L + 0.0000796096   ×   S   ×   L + 0.0070058   ×   FA 2 0.0000709692   ×   S 2 + 0.000880987   ×   L 2
The order of the factors influencing the wear rate can be established through analysis of the F-value magnitude as follows: it can be observed the F-value of FA (837.432) was the most significant parameter having the highest statistical influence on the wear rate of AA6063-FA composites followed by applied load (483.037) and sliding speed (286.805) and S, which implies that the FA content has greater effect on wear rate value than the L and S factor.
Adequate precision ratio, which represent the signal-to-noise ratio, with values higher than 4 indicate adequate model discrimination [17]. It compares the range between the predicted values at the design points to the average prediction error. As shown in Table 7, the adequate ratio was 59.006, which indicates that the model has capacity of navigating the design space. The R2 value of about 0.9791 indicates the variability of the response was 97.91% about the mean which proved the model provide a good fit for the data. This fact alongside with the remarkable residual analysis led to the conclusion that the model given in Equation (5) can predict the wear rate of the compound within the investigated range of all three parameters.
Furthermore, in Table 7, both the SL (BC) and S2 (B2) terms are considered to be insignificant as their p-value was greater than 0.05. However, these terms were automatically added in order to support the model’s hierarchy. In contrast, it is evident that the combined interaction of FA×L has the strongest effect on the response. The essential effects of the parameters on the responses, i.e., wear rate and the coefficient of friction, are further illustrated using contour lines and 3D surface plots.

Graphical Results of Wear Rate Model

Figure 1 presents the variance of collected data points about the linearized line which relates the predicted values with the actual values. This observation clearly indicates the model provide a good fit for the wear rate response.
Figure 2 further presents the perturbation plot of the wear rate data at the experimented mid values of the process parameters. The comparison of the effects of all three parameters at a specific area in the design range can be enabled by the perturbation plot. The response (wear rate) was graphically plotted by changing a single parameter over its investigated limits while other parameters were maintained constant. Similar to the significance of the p-value, this tool was also effective in identifying the interactive effects of parameters on the response [17]. The X-axis of the plot showed the relative position of the chosen levels of the parameters to the coded scale. The slopes of the curves showed the rate at which the parameter influenced the response. In Figure 2, the point selected in the design range was the central point (FA = 6 wt%, S 200 rpm and L 49 N).
The plot indicates that the fly ash content FA has the strongest inverse proportion (negative slope) on the wear rate while the load L has a smaller effect than FA, which was a direct proportional one (positive slope). The quadratic curves depicting the influence of the parameters on the wear rate were attributed to the quadratic terms in the model.
In summary, it is evident from Figure 2 that the FA content has negative effect on wear rate while both the applied load and sliding speed exhibit a positive impact on the wear rate. By increasing the FA content, the wear rate can be decreased while the opposite can be obtained when the applied load and sliding speed were increased. This trend was similar to that obtained from the experimental results.
The combined effects of two process parameters on the wear rate, while the third was kept constant, were further evaluated and these are demonstrated in Figure 3, Figure 4 and Figure 5. As described earlier for the perturbation plot, the third parameter could be changed, but again, it was kept at similar value selected for the perturbation figure. The red spheres on the graphs represent the experimental data points. Moreover, it is evident, especially on the 3D plots, that these spheres have close location to the surfaces, which indicate that the model fits properly with the data points. With respect to the effects of the combined parameters, they portrayed similar findings to those presented in the perturbation plot. However, by identifying their effects relative to the numerical values of the response, they offer a visual aid in selecting the desired ranges of the parameters. Furthermore, they were important in the optimization process of the ranges of the processing parameters so as to meet certain response criteria. Owing to these facts and to avoid redundancy, only brief discussions were given for each figure as provided in the following paragraph.
Figure 3a shows the contour plot of the effect of S and FA on the wear rate Wr. Each contour curve showed combinations of the two parameters that will predict a constant wear rate value. This wear rate value was observed inside the box situated at the middle of each contour curve. From the contour plot, S has a direct proportional effect on the predicted wear rate, while FA was inversely proportionate. The same data are further illustrated in 3D plot as shown in Figure 3b. The 3D plot has the added advantage of visual display of the effect of one of the parameter changes when the value of the other parameter is modified. For instance, consider the effect of S at two different values of FA which were at 0- and 12-wt%, it can be observed that the FA effect was stronger in the first case.
Figure 4a,b present that both S and L have directly proportional effects on the Wr. However, L has stronger effect than S. This indicated that variations in L led to larger variations in the Wr than those emanating from large variations in S.
Figure 5a,b shows that FA has an inversely proportional effect on the Wr, while L has a directly proportional effect. Similarly, FA has stronger effect on the response than L exhibits.

3.3. Modeling and Analysis of Variance for Coefficient of Friction

Likewise, the previous analyses described in Section 3.2 were repeated for the case of the coefficient of friction response. Again, the step-wise regression method was employed with the entire potential combinations of all three parameters and similarly too, the cubic terms were aliased for the ranges of the data obtained. The selection of the terms alongside the step wise regression method resulted in the elimination of the negligible model terms automatically [17]. Table 8 summarizes the computed ANOVA results with the variance for the Coefficient of friction model and all terms in the model.
In Table 8, it is shown that combined interaction of FA × L also had the strongest effects on the coefficient of friction. The effect of the parameters on the response is further demonstrated via the contour and 3D surface plots. The S2 and SL terms were proven to be insignificant as evidenced from the p-value which was greater than 0.05. Nonetheless, this term again was also to support the model hierarchy. The step-wise regression method detected that this model term was insignificant however, it was necessary by other significant model terms. Moreover, in Table 8, the adequacy measures of R2, adjusted R2 and predicted R2 were near to 1; they all agree well which is indicative of sufficient model [23]. The adequacy precision was greater than 4 which indicated sufficient model discrimination. In addition, this shows that the model has the capacity of navigating the design space [17]. The R2 value of approximately 0.9818 implies that around 98.18 of the data variability can be substantiated by model. Again, similarly to the wear rate, the model was proven to be a very good fit to the data and that the coefficient of friction, within the range of parameters investigated, can be predicted. Equation (6) presents the process model estimated from the experimental results within the range of parameters investigated.
F r = 0.74823 0.040431 × FA 0.00516264 × L 0.00103544 × S + 0.000144033 × FA × L + 0.0000293612 × FA × S + 0.00000411498 × S × L + 0.00101923 × FA 2 + 0.0000160301 × L 2 + 0.00000056855 × S 2
In summary, the model described by Equation (6) can provide good predictions of the coefficient of friction through any combinations of the fly ash content, sliding speed and load, within the range investigated.

Graphical Results of Coefficient of Friction

Figure 6 presents the variance of data points about the linearized line relating the actual and predicted values, which is an indication of the good fitness value of the model.
For the process parameters, the perturbation (see Figure 7) was obtained at a similar point of location (FA = 6 wt%, S 200 rpm and L 49 N). From the figure, it is evident that FA has inversely proportional effect and has the maximum effect on the coefficient of friction. The second maximum effect was owing to L which also had an inversely proportional effect on the coefficient of friction. The sliding speed also exhibited an inversely proportional effect.
The following Figure 8, Figure 9 and Figure 10, present the combined influences of two parameters, at a time on Fr. the design points are represented by the red spheres on the graphs; their locations on the 3D plots was a confirmation that both the model and the data points have a good fit. The contour plot for the influence of FA and S on Fr is presented in Figure 8a. It is evident from the contour plot that S and FA have an inversely proportional effect on Fr. It can be observed from the 3D plot in Figure 8b that as FA and S decrease the Fr decreases. By employing the combination of low S and low FA settings, the deepest Fr values were generated, while holding L constant, in this case at its mid-range = 49 N.
Figure 9a,b indicates that both FA and L have an inversely proportional influence on the Fr. Nonetheless, FA exhibited stronger effect than L. Slight changes in FA often affect the Fr in a magnitude greater than changes caused by even large changes in L.
Figure 10a,b shows that S and L have an inversely proportional effect on the Fr. In a similar way, S, which can be evidenced from the slopes of the contour lines and 3D surface, almost exhibited equal to L effect on the Fr.

3.4. Validation of the Statistical Models

The wear rate and coefficient of friction were numerically modeled by RSM based on FA content, applied load and sliding speed. Thereafter, the statistical model was validated with a dataset, 20% from the experimental data were used in the process of models building to verify the produced models and find the reliability of the model and the measured responses were recorded. The predicted wear rate and coefficient of friction calculated by Equations (5) and (6) were plotted against experimental data in Figure 11 and Figure 12. The figures indicated that the surface response method model provides excellent coincidence with experimental value with R-squared value of 0.9665 and 0.9893 respectively. These results were expected due to the high order model used in SRM which provide more accurate response.

4. Conclusions

AA6063 alloy-based aluminum matrix composites with FA particles have been successfully fabricated with compocasting technique. The overall results revealed that the wear resistance of the composites was positively enhanced by the inclusion of the FA particles. The wear resistance enhancement was attributed to increased hardness, generation of strain fields, homogenous distribution, spherical shape of FA particles and reduction in effective contact area. Increasing FA content restricted the deformation of the matrix material. The wear rates of the composites were increased by the applied load. FA content exhibited an inversely proportional relationship with wear rate, coefficient of friction and friction force for AA6063-FA composites. Increasing the sliding speed increased the wear rate but decreased the coefficient of friction and friction force of composites. In addition, the wear rate and fractional force of composites increased with increasing applied load, but coefficient of friction varied inversely with applied load. These control parameters were then employed to build statistical models for the wear rate and coefficient of friction. FA had the highest statistical influence on the wear rate and coefficient of friction of AA6063-FA composites followed by applied load and sliding speed. Moreover, FA content with applied load exhibited the strongest interaction effects. Finally, AA6063 aluminum matrix composites reinforced with FA can be used in applications that require excellent wear resistance.

Author Contributions

The research conceptualization and methodology are contributed by A.M.R., D.L.M. and U.M.B.; the experimental and numerical works including the data curation are performed by A.M.R. and U.M.B.; the results analysis and interpretation are conducted by A.M.R., D.L.M., U.M.B. and M.R.I.; the preparation of the original draft, reviews and editing are contributed by A.M.R. and D.L.M. and lastly the project administration, supervision and funding are managed by D.L.M. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to acknowledge the financial support from the Ministry of Education Malaysia with reference code FRGS/1/2019/STG07/UPM/02/10 and Vot No. 5540209.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Nomenclature

minMinute
rpmRevolution per minute
ASTMAmerican Society for Testing and Materials
FAFly ash
AlAluminum
AA6063AA6063 aluminum alloy
AA6063-FA AA6063 aluminum alloy-fly ash composites
MMCsMetal matrix composites
AMCsAluminum matrix composites
AMMCsAluminum metal matrix composites
Al-FAAluminum–fly ash
B4CBoron carbide
GraphiteGr
XRDX-ray diffraction
XRFX-ray fluorescence
ANOVAAnalysis of variance
RSMResponse surface method
DRAMCsDiscontinuously reinforced aluminum matrix composites
DOEDesign of experimental
CCDBased central composite design
SSSemisolid-semisolid
SLSemisolid-liquid
WrWear rate (mm3/m)
FrCoefficient of friction
µThe coefficient of friction
FTTangential force (N)
FNNormal force (N)
VVolume loss (mm3)
SSliding speed (rpm)
LApplied load (N)
WrpPredicted wear rate
FrpPredicted coefficient of friction
WreExperimental wear rate
FreExperimental coefficient of friction
WraAverage of experimental wear rate
FraAverage of experimental coefficient of friction
vol%Percentage of volume fraction
wt%Percentage of weight fraction

Appendix A

Table A1. Parameters used in mathematical models for wear analysis.
Table A1. Parameters used in mathematical models for wear analysis.
NOFA wt%Load NSliding Speed RpmWear Rate mm3/mCoefficient of Friction
1024.51506.48670.4723
2024.52009.250.4471
3024.52509.98890.3451
40491507.94840.3868
504920010.72830.346
604925011.33230.3069
7073.515011.01730.3187
8073.520012.26250.269
9073.525013.34710.2486
10224.51506.1360.4414
11224.52008.04480.3541
12224.52509.09340.2985
132491507.54840.3372
142492008.8756280.2755
1524925010.26730.2344
16273.515010.18340.2766
17273.520011.35030.2468
18273.525012.98390.1873
19424.51505.75240.3533
20424.52006.97380.3086
21424.52507.60220.2313
224491506.75110.2963
234492008.02540.2592
244492509.73570.1954
25473.51509.43890.2647
26473.520010.62260.2001
27473.525012.01390.166
28624.51505.32460.3135
29624.52006.5040.2676
30624.52506.9740.2129
316491506.02670.2589
326492006.72920.2155
336492507.51660.1707
34673.51508.08760.2391
35673.52008.84260.1754
36673.525011.03180.1415
37824.51504.54720.2694
38824.52005.838510.2409
39824.52506.42980.1954
408491505.70580.2214
418492006.22070.1991
428492507.05220.1535
43873.51507.03930.209
44873.52008.11050.1518
45873.52509.68930.1296
461024.51504.15750.2518
471024.52005.37480.2216
481024.52505.92170.1835
4910491505.18290.2093
5010492006.04140.1891
5110492506.35590.1306
521073.51506.50350.1899
531073.52007.57050.1427
541073.52508.3670.1167
551224.51503.46230.2324
561224.52004.5980.2059
571224.52505.580110.1607
5812491504.28910.1912
5912492005.85190.1553
6012492506.10120.1234
611273.51505.2213760.1701
621273.52006.32480.1217
631273.52507.83190.1085

References

  1. Sardar, S.; Pradhan, S.K.; Karmakar, S.K.; Das, D. Modeling of abraded surface roughness and wear resistance of aluminum matrix composites. J. Tribol. 2019, 141, 071601. [Google Scholar] [CrossRef]
  2. Nouri, M.; Abdollah-Zadeh, A.; Malek, F. Effect of welding parameters on dilution and weld bead geometry in cladding. J. Mater. Sci. Technol. 2007, 23, 817. [Google Scholar]
  3. Palanivel, R.; Mathews, P.K.; Murugan, N. Development of mathematical model to predict the mechanical properties of friction stir welded AA6351 aluminum alloy. J. Eng. Sci. Technol. Rev. 2011, 44, 25–31. [Google Scholar] [CrossRef]
  4. Heidarzadeh, A.; Saeid, T.; Khodaverdizadeh, H.; Mahmoudi, A.; Nazari, E. Establishing a mathematical model to predict the tensile strength of friction stir welded pure copper joints. Metall. Mater. Trans. B 2013, 44, 175–183. [Google Scholar]
  5. Lakshminarayanan, A.; Balasubramanian, V. Comparison of RSM with ANN in predicting tensile strength of friction stir welded AA7039 aluminium alloy joints. Trans. Nonferrous Met. Soc. China 2009, 19, 9–18. [Google Scholar] [CrossRef]
  6. Kumar, K.R.; Mohanasundaram, K.M.; Arumaikkannu, G.; Subramanian, R. Analysis of parameters influencing wear and frictional behavior of aluminum—Fly ash composites. Tribol. Trans. 2012, 55, 723–729. [Google Scholar] [CrossRef]
  7. Rajmohan, T.; Palanikumar, K.; Ranganathan, S. Evaluation of mechanical and wear properties of hybrid aluminium matrix composites. Trans. Nonferr. Met. Soc. China 2013, 23, 2509–2517. [Google Scholar] [CrossRef]
  8. Magibalan, S.; Senthilkumar, P.; Palanivelu, R.; Senthilkumar, C.; Shivasankaran, N.; Prabu, M. Dry sliding behavior of aluminum alloy 8011 with 12% fly ash composites. Mater. Res. Express 2018, 5, 056505. [Google Scholar] [CrossRef]
  9. Mandal, N.; Roy, H.; Mondal, B.; Murmu, N.C.; Mukhopadhyay, S. Mathematical modeling of wear characteristics of 6061 Al-alloy-SiCp composite using response surface methodology. J. Mater. Eng. Perform. 2012, 21, 17–24. [Google Scholar]
  10. Monikandan, V.V.; Jacob, J.C.; Joseph, M.A.; Rajendrakumar, P. Statistical analysis of tribological properties of aluminum matrix composites using full factorial design. Trans. Indian Inst. Met. 2015, 68, 53–57. [Google Scholar] [CrossRef]
  11. Moses, J.J.; Dinaharan, I.; Jekhar, S.J. Prediction of influence of process parameters on tensile strength of AA6061/TiC aluminum matrix composites produced using stir casting. Trans. Nonferr. Met. Soc. China 2016, 26, 1498–1511. [Google Scholar] [CrossRef]
  12. Razzaq, A.M.; Majid, D.L.; Ishak, M.R.; Uday, M.B. Microstructural characterization of fly ash particulate reinforced AA6063 aluminium alloy for aerospace applications. Mater. Sci. Eng. Conf. Ser. 2017, 270, 1. [Google Scholar]
  13. Razzaq, A.M.; Majid, D.L.; Ishak, M.R.; Uday, M.B. Influence of fly ash on the microstructure and mechanical properties of AA6063 alloy using compocasting technique. Mater. Express 2019, 9, 1–14. [Google Scholar] [CrossRef]
  14. Razzaq, A.M.; Majid, D.L.; Manan, N.H.; Ishak, M.R.; Uday, M.B. Effect of fly ash content and applied load on wear behaviour of AA6063 aluminium alloy. In IOP Conference Series: Materials Science and Engineering; IOP Publishing: Bristol, UK, 2018. [Google Scholar]
  15. Razzaq, A.M.; Majid, D.L.; Ishak, M.R.; Uday, M.B. Effects of Solid Fly Ash on Wear Behaviour of AA6063 Aluminum Alloy; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  16. Razzaq, A.M.; Majid, D.L.; Ishak, M.R.; Uday, M.B. Effect of fly ash addition on the physical and mechanical properties of AA6063 alloy reinforcement. Metals 2017, 77, 477. [Google Scholar] [CrossRef] [Green Version]
  17. Design-Expert Software. User’s Guide, Technical Manual; Stat-Ease Inc.: Minneapolis, MN, USA, 2000; Volume 6. [Google Scholar]
  18. Kumar, S.; Balasubramanian, V. Effect of reinforcement size and volume fraction on the abrasive wear behaviour of AA7075 Al/SiC p P/M composites—A statistical analysis. Tribol. Int. 2010, 43, 414–422. [Google Scholar] [CrossRef]
  19. Balasubramanian, V.; Lakshminarayanan, A.K.; Varahamoorthy, R.; Babu, S. Application of response surface methodolody to prediction of dilution in plasma transferred arc hardfacing of stainless steel on carbon steel. J. Iron Steel Res. Int. 2009, 16, 44–53. [Google Scholar] [CrossRef]
  20. Elangovan, K.; Balasubramanian, V.; Babu, S. Predicting tensile strength of friction stir welded AA6061 aluminium alloy joints by a mathematical model. Mater. Des. 2009, 30, 188–193. [Google Scholar] [CrossRef]
  21. Gopalakrishnan, S.; Murugan, N. Production and wear characterisation of AA 6061 matrix titanium carbide particulate reinforced composite by enhanced stir casting method. Compos. Part B Eng. 2012, 43, 302–308. [Google Scholar] [CrossRef]
  22. Kumar, S.; Balasubramanian, V. Developing a mathematical model to evaluate wear rate of AA7075/SiC p powder metallurgy composites. Wear 2008, 264, 1026–1034. [Google Scholar] [CrossRef]
  23. Badkar, D.S.; Pandey, K.S.; Buvanashekaran, G. Development of RSM-and ANN-based models to predict and analyze the effects of process parameters of laser-hardened commercially pure titanium on heat input and tensile strength. Int. J. Adv. Manuf. Technol. 2013, 65, 1319–1338. [Google Scholar] [CrossRef]
Figure 1. Experimental versus predicted wear rate mm3/m.
Figure 1. Experimental versus predicted wear rate mm3/m.
Crystals 10 00403 g001
Figure 2. Perturbation plot of the effect of Fly ash content, sliding speed and applied load on wear rate of AA6063-FA composites.
Figure 2. Perturbation plot of the effect of Fly ash content, sliding speed and applied load on wear rate of AA6063-FA composites.
Crystals 10 00403 g002
Figure 3. View of the interactive effect of FA and S on Wr, (a) contour lines 2D, (b) 3D.
Figure 3. View of the interactive effect of FA and S on Wr, (a) contour lines 2D, (b) 3D.
Crystals 10 00403 g003
Figure 4. View of the interactive effect of S and L on the Wr (a) contour lines 2D, (b) 3D.
Figure 4. View of the interactive effect of S and L on the Wr (a) contour lines 2D, (b) 3D.
Crystals 10 00403 g004
Figure 5. FA and L contour showing the interactive effect on Wr, (a) contour lines, (b) 3D.
Figure 5. FA and L contour showing the interactive effect on Wr, (a) contour lines, (b) 3D.
Crystals 10 00403 g005
Figure 6. Actual versus predicted scattering of the data points.
Figure 6. Actual versus predicted scattering of the data points.
Crystals 10 00403 g006
Figure 7. Perturbation plot of the process control variables’ effects on the coefficient of friction.
Figure 7. Perturbation plot of the process control variables’ effects on the coefficient of friction.
Crystals 10 00403 g007
Figure 8. View of the interactive effect of FA and S on Fr, (a) contour lines 2D, (b) 3D.
Figure 8. View of the interactive effect of FA and S on Fr, (a) contour lines 2D, (b) 3D.
Crystals 10 00403 g008
Figure 9. FA and L contour showing the interactive effect on the Fr, (a) contour lines 2D, (b) 3D.
Figure 9. FA and L contour showing the interactive effect on the Fr, (a) contour lines 2D, (b) 3D.
Crystals 10 00403 g009
Figure 10. View of the interactive effect of S and L on Fr, (a) contour lines 2D, (b) 3D.
Figure 10. View of the interactive effect of S and L on Fr, (a) contour lines 2D, (b) 3D.
Crystals 10 00403 g010
Figure 11. SRM-predicted wear rate versus experimental results.
Figure 11. SRM-predicted wear rate versus experimental results.
Crystals 10 00403 g011
Figure 12. SRM-predicted coefficient of friction versus experimental results.
Figure 12. SRM-predicted coefficient of friction versus experimental results.
Crystals 10 00403 g012
Table 1. Chemical composition of tested AA6063 in wt% [16].
Table 1. Chemical composition of tested AA6063 in wt% [16].
ElementsCuFeMgMnSiTiZnCrAl
wt%< 0.010.190.510.040.860.0080.010.01Balance
Table 2. Chemical composition of fly ash (FA) in wt% [16].
Table 2. Chemical composition of fly ash (FA) in wt% [16].
ElementsSiO2Al2O3Fe2O3TiO2K2OCaOMgONa2OCuO
wt%59.9819.092.781.141.090.630.380.340.01
Table 3. Levels of factors.
Table 3. Levels of factors.
No.FactorNotationUnitLevels
1Fly ash contentFAwt%024681012
2LoadLN24.54973.5____
3Sliding speedSRpm150200250____
Table 4. Factorial design summary (Design-Expert V10.0)
Table 4. Factorial design summary (Design-Expert V10.0)
FactorNameUnitsTypeSubtypeMinMaxCoded ValuesMeanStd. Dev.
AFly ash content wt%NumericDiscrete012−1.000 = 01.000 = 126.156864.05646
BSliding speedRpmNumericDiscrete150250−1.000 = 1501.000 = 250197.05940.5114
CLoad NNumericDiscrete24.573.5−1.000 = 24.51.000 = 73.549.480419.898
ResponseNameUnitsObsAnalysisMinMaxMeanStd. Dev.RatioTransModel
R1Wear ratemm3/m51Polynomial3.462313.34717.631282.443593.85498NoneQuadratic
R2Coefficient of friction51Polynomial0.10850.47230.2367670.08307694.353NoneQuadratic
Table 5. Statistical summary of wear rate model.
Table 5. Statistical summary of wear rate model.
SourceStd. Dev.R-SquaredAdjusted R-SquaredPredicted R-SquaredPRESSRecommendation
Linear0.600.94370.94010.933819.76
2FI0.480.96560.96090.954813.50
Quadratic0.390.97910.97450.965810.20Suggested
Cubic0.360.98580.97840.963310.97Aliased
Table 6. Statistical Summary of Coefficient of Friction Model.
Table 6. Statistical Summary of Coefficient of Friction Model.
SourceStd. Dev.R-SquaredAdjusted R-SquaredPredicted R-SquaredPRESSRecommendation
Linear0.0240.92140.91640.90430.033
2FI0.0200.94830.94120.92420.026
Quadratic0.0120.98180.97780.96990.010Suggested
Cubic0.0094060.99150.98720.97690.007958Aliased
Table 7. ANOVA results for the quadratic model of wear rate.
Table 7. ANOVA results for the quadratic model of wear rate.
Sum of MeanFp-Value
SourceSquaresdfSquareValueProb > F
Model292.32932.48213.39<0.0001Significant
A-weight fraction127.461127.46837.43<0.0001
B-sliding speed43.65143.65286.80<0.0001
C-Load73.52173.52483.04<0.0001
AB1.3611.368.960.0047
AC3.4913.4922.95<0.0001
BC0.2010.201.300.2600
A20.4810.483.180.0819
B20.3510.352.330.1348
C23.1713.1720.83<0.0001
Residual6.24410.15
Cor Total298.5650
Std. Dev.0.39 R-squared0.9791
Mean7.63 Adj R-squared0.9745
C.V. %5.11 Pred R-squared0.9658
PRESS10.20 Adequate Precision59.006
Table 8. ANOVA for Response Surface Quadratic model (coefficient of friction model).
Table 8. ANOVA for Response Surface Quadratic model (coefficient of friction model).
Sum of MeanFp-Value
SourceSquaresdfSquareValueProb > F
Model0.3490.038245.70< 0.0001Significant
A-Fly ash content0.1910.191219.47< 0.0001
B-Load0.06910.069451.39< 0.0001
C-Sliding speed0.06110.061399.07< 0.0001
AB6.376 × 10−316.376 × 10−341.61< 0.0001
AC1.054 × 10−311.054 × 10−36.880.0122
BC5.305 × 10−415.305 × 10−43.460.0700
A20.01010.01066.91< 0.0001
B21.049 × 10−311.049 × 10−36.850.0124
C22.273 × 10−512.273 × 10−50.150.7021
Residual6.282 × 10−3411.532 × 10−4
Cor Total0.3550
Std. Dev.0.012 R-squared0.9818
Mean0.24 Adj R-squared0.9778
C.V.%5.23 Pred R-squared0.9699
PRESS0.010 Adeq Precision66.605

Share and Cite

MDPI and ACS Style

Mohammed Razzaq, A.; Majid, D.L.; Ishak, M.R.; Muwafaq Basheer, U. Mathematical Modeling and Analysis of Tribological Properties of AA6063 Aluminum Alloy Reinforced with Fly Ash by Using Response Surface Methodology. Crystals 2020, 10, 403. https://doi.org/10.3390/cryst10050403

AMA Style

Mohammed Razzaq A, Majid DL, Ishak MR, Muwafaq Basheer U. Mathematical Modeling and Analysis of Tribological Properties of AA6063 Aluminum Alloy Reinforced with Fly Ash by Using Response Surface Methodology. Crystals. 2020; 10(5):403. https://doi.org/10.3390/cryst10050403

Chicago/Turabian Style

Mohammed Razzaq, Alaa, Dayang Laila Majid, Mohamad Ridzwan Ishak, and Uday Muwafaq Basheer. 2020. "Mathematical Modeling and Analysis of Tribological Properties of AA6063 Aluminum Alloy Reinforced with Fly Ash by Using Response Surface Methodology" Crystals 10, no. 5: 403. https://doi.org/10.3390/cryst10050403

APA Style

Mohammed Razzaq, A., Majid, D. L., Ishak, M. R., & Muwafaq Basheer, U. (2020). Mathematical Modeling and Analysis of Tribological Properties of AA6063 Aluminum Alloy Reinforced with Fly Ash by Using Response Surface Methodology. Crystals, 10(5), 403. https://doi.org/10.3390/cryst10050403

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop