Abstract
The present study focuses on the experimental investigation and optimization of compressive residual stress on ball burnished Mg Ze41A alloy. The optimum compressive residual stress value of −182 MPa is obtained at burnishing force of 50 N, speed of 1200 r/min, feed of 150 mm/min and 4 number of passes. Burnishing force and speed are found to be the most sensitive parameters in affecting the compressive residual stress, followed by number of passes and burnishing feed. The demarcation lines in residual stress-burnishing maps of Mg Ze41A alloy are non-parallel and non-uniform, indicating significant interactions among burnishing parameters. Ball burnishing process imparted desired amount of work hardening, which resulted in homogenized (σx1 = σx2) and isotropic stress (σ1 = σ2) state. Additionally, the fuzzy logic model has been deployed to estimate the compressive residual stress values at different burnishing conditions. Fuzzy logic model estimated the results with an error of 6.11%. The results demonstrate that ball burnishing process is an effective method in inducing compressive residual stresses on Mg alloy. The outcome of this study can be of great help in enhancing the fatigue performance of small size products.
Introduction
The behaviour of metallic components against fatigue performance in automotive, aerospace and power generation industries is determined by fatigue strength. Hence, these industrial sectors always aim at improving the fatigue limit to avoid crack nucleation and fatigue failure of critical components. 1 The most popular methods to improve the service life of the components against cyclic loading are mechanical surface treatments. 2 They are being used widely to enhance the physical, mechanical and metallurgical properties of cyclic loaded components. The plastic deformation during mechanical treatments results in compressive residual stresses (CRSs), work hardening effect and microstructure changes; thereby improving the fatigue strength significantly. Eventually, these mechanical methods hinder the crack initiation and propagation. On the other side, they also improve corrosion and wear resistance. Among the mechanical surface treatments, burnishing, shot peening, laser shock peening and laser cladding are the popular methods. 3 Ball burnishing is a simple, rapid, flexible and economical mechanical surface treatment, which can be performed by using a computer numerical control (CNC) or vertical machining centre (VMC) machine. The process involves the plastic flow of material from “peaks or ridges” to “valleys or depressions” by virtue of small plastic deformation. Ball burnishing process has multiple benefits, as depicted in Figure 1. Ball burnishing process improves the surface finish, surface hardness and induces high CRSs on the surface. As a result, significant improvement in fatigue life, corrosion and wear behaviour can be observed. Also, it can be realized that ball burnishing process can replace shot peening, laser shock peening and laser cladding in inducing CRS. 4

Benefits of ball burnishing process.
Ball burnishing process found to induce CRS up to 1 mm, as depicted in Figure 2. The CRS is not maximum on the surface but at a certain depth. In order to balance the stresses in the material, tensile residual stresses are formed after the CRS. 5 Numerous researchers have explored the potentiality of burnishing process in inducing the CRS on a material. For instance, the effect of burnishing parameters in inducing CRS has been investigated on AISI 8620 steel. 6 It has been observed that the burnishing force and feed have strong influence on CRS while the effect of speed is nominal. The CRS is highest at highest burnishing force and lowest feed. This might be due to the predominant plastic deformation at higher burnishing force and desirable surface roughness at lower burnishing feed. 6 Similarly, ball burnishing process has been deployed on 17-4 PH stainless steel to impart CRS on the surface. Burnishing force is found to be the highly influential factor affecting CRS. 7 A maximum CRS value of −994.6 MPa is observed at a depth of 0.1 mm and also, it has been realized that the residual stresses are more compressive in circumferential direction than axial direction for a cylindrical workpiece. 7 Similar conclusions are obtained during slide diamond burnishing of 41Cr4 steel. Burnishing force has the strongest effects on CRS and subsequently on fatigue strength and fatigue life. 8 The state of residual stress during roller burnishing of hot-rolled Mg–12Gd–3Y alloy has been studied to determine its significance against fatigue properties. 9 It is observed that the maximum CRS (−410 MPa) induced on the hot-rolled surface is greater than that of tensile strength (362 MPa).

Compressive residual stress (CRS) of ball burnishing process.
The number of passes during burnishing process has an effect on CRS as well. For instance, ball burnishing of EN series steels have revealed that the CRS increases with increase in passes up to certain extent. A maximum CRS of 280 MPa has been observed after two passes in EN series steels. 10 Apart from passes, burnishing speed also has its prominence on CRS. This has been observed during burnishing of TA2 alloy, where the CRS value increases significantly with decrease in burnishing speed. The maximum CRS of −501 MPa is obtained at burnishing speed of 150 mm/min. 11 Not only burnishing parameters but also the direction of motion of the burnishing tool on the surface has an effect on CRS. The amount of CRS induced on the surface or sub-surface during burnishing process depends on nature of the material. The magnitude of CRS will be higher for harder materials and lower for soft materials. But, the depth of CRS will be higher for softer materials than that of harder materials. 12 The assistance of ultrasonic vibration during burnishing will have a significant change on the intensity of CRS than the conventional burnishing. For example, ultrasonic assisted ball burnishing of AA6061-T6 alloy induced a maximum CRS of −217 MPa, while the conventional ball burnishing process induced −163 MPa. 13 In another example, ultrasonic assisted ball burnishing has induced a maximum CRS of 110 MPa on AA6061-T6 alloy at a penetration depth of 0.7 mm, static load of 120 N, vibration amplitude of 14 μm and tungsten carbide ball diameter of 3 mm. 14 Ball burnishing process can be a useful process in transforming the TRS induced during machining process such as turning, milling and grinding etc., into CRS. For instance, the TRS of 1223 MPa formed during turning of AISI 4130 steel has been transformed into CRS of −223.7 MPa by deploying ball burnishing process. 15 Also, it has been observed that the ball burnished surface delayed the fatigue crack propagation in stress concentration zones, thereby improving the fatigue strength drastically. 16
In recent times, machine learning techniques have become quite popular in modelling complex engineering problems precisely and are being regarded as a substitution for classical modelling methods. 17 The popular machine learning techniques can be artificial neural network, adaptive neuro-fuzzy interface system, Fuzzy Logic, Cellular Automata, Hybrid Systems, etc. The fuzzy logic model is an inexact interpretation of an engineering problem by using human linguistic conditions. The fuzzy logic model predicts the response characteristics by a set of fuzzy rules. The fuzzy logic rules are developed by establishing the inter-connection between input variables and response characteristics. Therefore, the fuzzy rules are formulated by the degree of truth rather than Boolean logic. 18 With the ability of estimating the response characteristics quite effectively, the fuzzy logic model has been instrumental for wide range of industrial applications. For instance, the material removal rate of high speed steel during electrochemical machining has been estimated by fuzzy logic model effectively. The model predicted the material removal rate with an accuracy of 97.25%. 19 The choice of linguistic terms and membership function (MF) decide the precision of fuzzy logic model. Fuzzy logic modelling of multi-responses of AISI 1095 steel during electric-discharge machining has revealed that the best predicted results are obtained by considering triangular MF for input parameters and trapezoidal MF for output responses. 20 The UTS and elongation of friction stir welded copper joints has been estimated by fuzzy modelling. The maximum UTS of 276.1 MPa and elongation of 44.6% is obtained by fuzzy model. 21 Similarly, mechanical properties of abaca fibre reinforced polymer composites is predicted by fuzzy logic with an accuracy of 87%. The results obtained are slightly inaccurate as the input data is not sufficient enough to model the problem. Therefore, sufficient data sets are to be created before adopting fuzzy logic modelling. 22
From the state of the art review, it can be understood that the experimental investigation on CRS during ball burnishing process is very limited. In fact, there are no studies related to the optimization of ball burnishing process for maximum CRS. Also, the effect of major burnishing parameters on CRS has been studied to an extent, but deep insights into the relationship are yet to be reported. The rare earth base Mg alloy, Ze41A alloy was originally developed for gear box applications in the aerospace industry. However, with the nature of biodegradability and biocompatibility, it is now being explored for biomedical applications such as orthopaedic implants and cardiovascular stents. In the current study, attention is paid to examine the capability of ball burnishing process in inducing CRS on Mg Ze41A alloy at different burnishing conditions. Also, the effect of each burnishing parameter on CRS has been explored. An attempt has been made to optimize the ball burnishing process for maximum CRS using Taguchi method. Finally, the popular soft computing technique, fuzzy logic is deployed to estimate the CRS value of ball-burnished magnesium alloy in different test cases. The predictive capability of the technique has been determined by a validation study.
Methodology
CNC burnishing process
The ball burnishing process has been carried out in a vertical head CNC machine. The process is deployed on rare earth base magnesium (Ze41A) alloy, which has a chemical composition of Zirconium 0.54%, Cerium 0.61%, Zinc 5.4%, Manganese 0.02%, Iron 0.003%, and rest Magnesium. The workpiece of dimensions 70 × 70 × 25 mm3 has been parted off from the casted block of Mg Ze41A alloy with a plan of executing 25 burnishing experiments, as shown in Figure 3(a). Each burnishing experiment is performed on 8 × 8 mm2 and parted by 5 mm on the workpiece consecutively. After the ball burnishing operation, each burnishing region is separated by wire-EDM cut, as shown in Figure 3(b). Now, each sample has been evaluated for CRS on the burnished surface (Figure 3(c)). It is important to mention that the workpiece is milled before performing burnishing operation, such that the surface roughness of Ra = 0.941 μm has been obtained initially. Figure 4 represents the innovative ball burnishing tool used for experimentation. The shank with circular cross section has been designed for mounting the tool into the tool holder and the mandrel is designed in such a way that it holds the rotating tungsten carbide ball (ϕ8 mm) firmly. The rotating ball will be assisted during the contact with workpiece by three tungsten balls (ϕ4 mm).

(a) Design of workpiece, (b) burnished samples, (c) burnished sample geometry.

Ball burnishing tool for burnishing operation.
The ball burnishing experiments are conducted on ‘CNC Fanuc series oi Mate-MC’, which is having an accuracy of 1 μm. Ball burnishing process demands for such high precession CNC machines as burnishing force is highly sensitive parameter. The other specifications include spindle power 11 kW and maximum spindle speed 4000 r/min. The preliminary experiments have provided great assistance in selecting the range of each burnishing parameter. The burnishing parameters and their corresponding levels are provided in Table 1. Each burnishing parameter has been varied in five levels to apprehend the non-linear interactions in ball burnishing process appropriately. 4 An illustration of experimental set-up has been showcased in Figure 5(a)–(c). The workpiece is fixed firmly on the CNC machine to ensure no fluctuations during operation. The initial milled and burnished surfaces are cleaned with alcohol to remove foreign contaminants and hard particles from the surface. The initial milling experiments are also performed on the same CNC machine with a speed of 1000 r/min, feed of 100 mm/min and depth of cut 1 mm. Surface integrity measurements after milling operation has revealed the microhardness of 75.2 HV and surface roughness (Ra) of 0.941 μm.

An illustration of (a), (b) experimental set-up in VMC (c), (d) workpiece and burnished samples after wire-EDM cut.
The levels of the burnishing parameters.
CRS measurement
Compressive residual stress during ball burnishing process has been evaluated by PANalytical X-pert diffractometer. X-rays of wavelength (λ = 0.1542 nm) are being generated by applying current of 35 mA and voltage of 40 kV with copper (Cu) as the anode material. During burnishing operation, the surface will be subjected to stress, which results in crystal lattice structure changes due to the plastic flow of the material. This results in change of crystallographic interplanar spacings between the planes. These crystallographic spacings behave as internal strain gauges and change the strain consequently based on Bragg's law of diffraction. Bragg's angle of 2θ = 104.37° and the corresponding crystallographic planes {1 1 2} have been used to determine residual strains. The X-rays are diffracted on the sample at different tilt angles (ψ), and the corresponding crystallographic interplanar distance (dψ) is recorded and then sin2 ψ method is deployed to evaluate residual stresses. The residual stress from sin2 ψ method has been evaluated by equation (1).
Results and discussion
The Taguchi method is a powerful and robust tool for process optimization in an industry. It designs the experiments using an orthogonal array and subsequently, the analysis is carried out by using signal-to-noise ratio (S/N). The quantities ‘signal’ and ‘noise’ indicate the desirable and undesirable value of the performance characteristics respectively. The optimum process condition derived from the S/N ratios of Taguchi method is found to be stable, robust, and results in minimum process variation. In the current investigation, Taguchi method with L25 orthogonal array has been used for ball burnishing process optimization. The orthogonal array is designed by considering the levels of parameters, as depicted in Table 1. The main focus of the study is to maximize the CRS induced during ball burnishing process. Therefore, ‘larger the better’ type characteristic is chosen for obtaining the best CRS value. Positive value of CRS has been considered for computing S/N ratio of CRS. The experimental data and S/N ratio of each test have been presented in Table 2. The S/N ratio of CRS has been evaluated from equation (2).
Taguchi results of ball burnishing process.
Role of burnishing force in inducing CRS
The influence of burnishing force on CRS has been depicted in Figure 6. It can be observed that the CRS value increases with an increase in burnishing force up to 50 N and decreases further. When the burnishing force is low, the penetration of the rotating ball into the work surface will be very little, thereby cause incomplete deformation of the asperities. As a result, the strain hardening induced on the burnished surface is low and hence, the CRS value is low. As the force increases steadily, the protrusion of the rotating ball will be high and the deformation zone widens, resulting in high-strain hardening. This phenomenon has been observed up to a burnishing force of 50 N. Therefore, the maximum value of CRS has been obtained at 50 N burnishing force. However, the CRS value decreases as the rotating ball advances further, due to severe plastic deformation and excessive strain hardening of the burnished layers of Mg Ze41A surface. This leads to distortion and flaking at extremely high burnishing forces and apparently the value of CRS is low. Therefore, the CRS value decreases beyond 50 N burnishing force. The earlier reports on the burnishing of non-ferrous materials have revealed that the CRS value has been reduced up to a certain burnishing force. Experimental exploration on ball burnishing of Mg Ze41A alloy also compliments these findings. The interactions between the burnishing force and other burnishing parameters are presented in Figure 7. The crossed lines in the interaction plot represent that the burnishing force interacts with the other burnishing parameters significantly. Similar findings are also observed during analysis of variance (ANOVA) analysis, as described in the ‘ANOM and ANOVA of CRS’ section. Hence, the burnishing force has to be rightly chosen for achieving the best CRS value. Also, a comprehensive analysis related to interaction of burnishing force and burnishing speed has been provided in section ‘RsB maps’ as residual stress-burnishing (RsB) map.

Main effects plot of compressive residual stress (CRS).

Interaction effects plot of compressive residual stress (CRS).
Role of burnishing speed in inducing CRS
The effect of burnishing speed on CRS value can be found in Figure 6. The CRS value at the initial burnishing speed is observed to be lower. Although the penetration time of the burnishing ball is high at low speeds the plastic deformation induced on the surface is quite less. The burnishing effect has been found to impart elastic deformation at low burnishing speeds. As a result, the elastic recovery is found to be higher and subsequently low work hardening and thereby low residual stress. However, with the increase in burnishing speed, the CRS value increases. This is perhaps due to the increase in plastic deformation and the lateral material flow of the strain hardened surface. Also, the interfacial temperature at the junction of the rotating ball and the workpiece rises. All these factors contribute for the high residual stresses on the burnished surface. The results reveal that the burnishing speed of 1200 r/min has yielded the best CRS value. But, the CRS value decreases with increase in burnishing speed further. This might be due to excessive strain hardening and lateral material flow. At higher burnishing speeds, the lateral material flow will be higher due to heavy centrifugal action. As a result, the material gets dispersed laterally to long distances, which results in excessive strain hardening. Consequently, the surface gets damaged and results in low CRS value. Therefore, it is advisable to adopt medium burnishing speed for achieving the best CRS value. The reduction in CRS value at higher burnishing speeds has been observed during burnishing of AISI 8620 Steel 6 and Mg-0.8Ca alloy, which closely concur with the results of Mg Ze41A alloy. The interaction plot (Figure 8) of Mg Ze41A alloy showcases very severe interactions between burnishing speed and other process parameters as the lines are crossed at multiple locations. Therefore, the burnishing speed has to be chosen accurately for achieving the best CRS value. In order to elucidate the interaction further, a detailed discussion has been presented in the ‘RsB maps’ section as RsB map.

RsB maps of ball burnished Mg Ze41A alloy (the unit of CRS is MPa). RsB: residual stress-burnishing; CRS: compressive residual stress.
Role of burnishing feed in inducing CRS
Along with burnishing force and speed, burnishing feed also has an influence on surface roughness, microhardness, and CRS during ball burnishing operation. The role of burnishing feed in inducing CRS has been explicitly represented in Figure 6. The lowest CRS value at low burnishing feed is pertaining due to excessive cold work and large plastic deformation. This might be due to the burnishing ball deformations at a particular place are higher at lower burnishing feed. As a result, the longitudinal material flow is higher, which resulted in severe strain hardening. However, with increase in burnishing feed, the CRS value increases and attains maximum value at burnishing feed of 150 mm/min. The longitudinal material flow decreases with increase in burnishing feed, which results in desirable plastic deformation and cold work. Therefore, the ball burnishing process yields the best CRS value at medium burnishing feed. On the other hand, the CRS value is lower at higher burnishing feeds. This could be due to insufficient burnishing at higher burnishing speeds. This results in less material flow in longitudinal direction and inadequate cold work. Also, at high burnishing feeds, incoherent deforming action of the rotating ball has been observed, which resulted in instability of burnishing tool. Burnishing feed has been found to interact with the most influential burnishing parameters such as burnishing force and speed profoundly, as depicted in Figure 7. Therefore, the burnishing feed has to be chosen in synchronize to burnishing force and speed. Similar research work related to ball burnishing of engineering materials support the findings of present study on Mg Ze41A alloy.
Role of passes in inducing CRS
The elastic recovery during plastic deformation of ball burnished surface can be monitored by number of passes. Hence, it is always recommended to adopt multiple passes during burnishing of metals. However, it is essential to know the influence of multiple passes on CRS. Figure 6 shows the state of residual stress with increase in passes. The elastic recovery in plastic zone can not be prevented by single pass. Therefore, the intensity of strain hardening will be lowered, resulting in lower cold work. Apparently, the lower CRS value has been attained at single pass. As the passes increase, the elastic recovery increases and thereby the intensity of strain hardening increases, subsequently the cold work will be higher. Therefore, increase in the CRS value has been observed for multiple passes and this phenomenon continues till four passes. This shows that the elastic recovery in plastic zone has been prevented completely by four passes. On the other hand, the CRS value decreases after four passes. This could be due to over-strain hardening effect and flaking of the surface. The interaction plot (Figure 7) reveals that the interactions of passes with other burnishing parameters is reasonable and has to considered during ball burnishing operation. The experimental investigation on ball burnished Mg Ze41A alloy establishes that the intensity of CRS increases with increase in burnishing force, speed, feed and passes, reaches the highest value and decreases further. Also, it has been found that the basic characteristics of CRS for burnished MgCa0.8 alloy in the literature can be transferable to Mg Ze41A alloy.
Analysis of means and ANOVA of CRS
The analysis of means (ANOM) has been carried out by using S/N ratios of CRS to decide the optimal level of parameters in the ball burnishing process. The optimum level of the burnishing parameter is chosen based on the highest S/N ratio. The ANOM results based on the S/N ratios of CRS have been presented in Table 3. It can be revealed that the best CRS value has been obtained at 50 N burnishing force, 1200 r/min burnishing speed, 150 mm/min burnishing feed and four passes. Further, the ANOVA has been carried out based on the S/N ratio of CRS value to elucidate the prime importance of each burnishing parameter and the results are tabulated in Table 4. The analysis is performed at confidence level of 95%, such that the level of significance is maintained as 0.05. In the present investigation, the significance level is found to be less than 0.05, which indicates that the burnishing parameters will have a direct influence on CRS value. The results emphasize that burnishing force (26.78%) and burnishing speed (26.08%) have the high relative importance in maximizing the CRS value; while burnishing feed and number of passes have low-relative importance.
ANOM of CRS.
ANOM: analysis of mean; CRS: compressive residual stress.
ANOVA of CRS.
ANOVA: analysis of variance; CRS: compressive residual stress.
Evaluation of optimum results
The optimum level of burnishing parameters has been determined from ANOM results of CRS. However, the optimum burnishing condition obtained for the best CRS value has not been in the list of L25 orthogonal array. Therefore, the CRS value at optimum burnishing condition has been estimated using S/N ratios. The S/N ratio (ηopt) of CRS at optimum CRS condition is estimated by equation (3).
Optimum burnishing parameters for compressive residual stress (CRS).
RsB maps
The CRS of Mg Ze41A alloy is profoundly affected by burnishing force and speed, as depicted in ANOVA table. Further, the literature lacks in analytical modelling of residual stresses during burnishing process. Hence, an effort is made to correlate the residual stresses with regards to highly influencing burnishing parameters. Therefore, the RsB maps are designed based on burnishing force and speed. This map is designed by considering burnishing force on abscissa and burnishing speed on ordinate. When a hard rotating ball moves over a surface repeatedly, the plastic deformation occurs, which in turn induces strain hardening effect and results in residual stress. However, the intensity of residual stress depends on mechanical load, material deformation and material flow of work piece. This residual stress will be isotropic in nature and it is subjected on all the grains because of the grain constraints in each direction. The RsB map of ball burnished Mg Ze41A alloy is depicted in Figure 8. This map basically represents the interactions encountered during ball burnishing operation. The demarcation lines are highly non-parallel and non-uniform, indicating the dominant interactions during burnishing of Mg Ze41A alloy.
The value of CRS is varied from one burnishing condition to another. Therefore, it is essential to identify the range of parameters for best CRS value. An attempt has been made to propose the range of parameters for best CRS value by designing burnishing zones, as represented in Figure 9. The burnishing zones are categorized based on the range of CRS values marked as Unlikely zone, Likely zone and Most Likely zone. The unlikely zone is observed to be predominant at extreme burnishing force and speed conditions. At low extreme conditions, the intensity of strain hardening is quite low because of the poor burnishing effect. In contrary, the burnished surface gets distorted due to excessive strain hardening at higher extreme conditions. As a result, the burnished Mg Ze41A surface exhibited low CRS value at extreme burnishing conditions. However, a clear transition has been observed from unlikely to likely zone, as moved away from unlikely zone. This could be possible because of improvement in the strain hardening, when shifted from extreme conditions. However, the intensity of strain hardening is not at the optimum level in likely zone, which resulted in moderate values of CRS. Meanwhile, the shift from likely to most likely zone is inevitable, when drifted towards the centre. The amount of strain hardening induced in likely zone is of desirable value, which resulted in best CRS value. It can be concluded that the best CRS value can be achieved, when the ball burnishing operation is performed in most-likely zone. Further, these burnishing zones can assist in process parameters selection for any equivalent strength material. The most-likely zone for an equivalent strength material (σy = 140 MPa) can be found in between 45 and 55 N of burnishing force and 1150–1250 r/min for burnishing speed.

Burnishing zones based on residual stress.
Fuzzy logic modelling
The information in FIS is processed by linguistic reasoning. A linguistic term describes the state of a phenomenon by a word rather than numerical values. However, the words are classified further into numerical values by the help of fuzzy sets. In the present study, the CRS of magnesium alloy is modelled by considering four major ball burnishing parameters. The linguistic values of ball burnishing parameters are represented in Table 6. The input parameters are partitioned according to the levels of Taguchi method. Similarly, the linguistic terms of CRS are presented in Table 7. CRS has been fuzzified into seven linguistic terms to ensure great predictive response. The fuzzy sets have been represented graphically by MFs. Also, the MFs can be used to present the degree of fuzziness in a fuzzy set by mapping the value between 0 and 1 by a curve in the input space. The selection of MF and its value is exclusively based on the process knowledge and the experimental data. Triangular MFs are generally preferred as it has been gradually decreasing on either sides with only one definite value. In the current study, each burnishing parameter and CRS has been modelled by triangular MF. A triangular MF, μA(x), is defined in equation (4). The triangular MF is found to be suitable as it suffices to evaluate the slope and the central value. As a result, burnishing parameters and CRS values are constrained to two fuzzy sets and hence, their summation always remains unity. In addition, the triangular MF outperformed the other MF's in estimating the results of multi-performance characteristic problems. Fuzzy logic model in Matlab-R2016a has been developed by mamdani type inference system.
Linguistic values of burnishing parameters.
Linguistic values of compressive residual stress (CRS).
The linguistic terms of burnishing parameters and CRS are interrelated by heuristic rules as If ‘condition’ then ‘consequence’. A rule base comprising 25 fuzzy rules has been formulated by considering the experimental data. For example, Rule 15: If burnishing force is medium, burnishing speed is very large, burnishing feed is small and passes are large then the surface roughness is good.
Prior to defuzzification process, the outcome of the fuzzy rule base has been integrated by aggregation process. The aggregation method used for unification of fuzzy rule bank is max (maximum) method for its easiness and adoptability. Additionally, the MF of the ‘consequence’ part of fuzzy rule base is unified by fuzzy-union operator. Further, the cumulative response of fuzzy rules obtained through unification process undergoes defuzzification. The term ‘defuzzification’ represents the conversion of fuzzified value into real value. The cumulative response will be transformed into crisp value. The crisp value represents the CRS of ball burnished MgZe41A alloy. The commonly used defuzzification methods are mean of maxima method, bisector method and centroid method. In the present investigation, centroid method is used, for its adoptability, execution and accuracy. The governing equation for centroid method has been provided in equation (5). The experimental data and aforesaid linguistic terms with triangular MF have been used to simulate the fuzzy logic model in the Matlab software on the basis of mamdani interface system. The generated model has been validated by predicting the values at three unknown test cases and the outputs generated are represented in Figure 10.

Prediction results of fuzzy logic model.
where xi is centroid value, x*is crisp value, μ(xi) is firing level of the ith rule among ‘n’ rules fired.
Validation of FIS model
After the successful generation of fuzzy logic in Matlab software, the models are tested for accuracy. The CRS values at three unknown burnishing conditions have been predicted by the models, as represented in Table 8. The predicted CRS values are compared with experimental values to validate the models. The percentage error for the results has been calculated to compare accuracy between the models. The average absolute percentage error for the developed fuzzy logic model is found to be 6.11%, as depicted in Table 8. The least errors prove that the fuzzy logic results are nearer to the actual experimental results. The close assent of predicted and measured values of CRS emphasizes that the fuzzy logic model can be effective in predicting the residual stresses under different burnishing parameters satisfactorily. Thus, the developed fuzzy logic model can be a promising alternative for predicting the CRS value in the specified range of burnishing parameters. The number of arithmetic operations is substantially low during implementation of fuzzy logic model. Similar estimation capabilities of fuzzy logic models were observed during machining of ferrous materials.
Validation of FIS model.
CRS: compressive residual stress; FIS: fuzzy interface system.
Conclusion
Ball burnishing process was found to impart CRSs on the surface successfully. Based on the results, the following conclusions are derived:
Burnishing force and speed are found to be the most sensitive parameters in affecting the CRS with a contribution of 26.78% and 26.08%, respectively, followed by burnishing passes and feed with a contribution of 21.97% and 19.97%, respectively. The optimum CRS value of −182 MPa has been obtained at burnishing force of 50 N, speed of 1200 r/min, feed of 150 mm/min and 4 passes. The demarcation lines in RsB maps are found to be non-parallel and non-uniform, indicating significant interactions among burnishing parameters during ball burnishing of Mg Ze41A alloy. Based on RsB maps, it has been proposed that the ball burnishing process for an equivalent material (σy = 140 MPa) can be performed at burnishing force of 45–55 N and burnishing speed of 1150–1250 r/min for obtaining the best CRS value. The fuzzy logic model with triangular MF and centroid method of defuzzification has yielded the results with the lowest error of 6.11%. The FIS model can be of great help in generating the global database for process improvements during CRS estimation of ball burnished Mg Ze41A alloy. The ball burnishing process is effective and economical in inducing the CRS and hence, can be recommended for industrial applications, where fatigue strength and corrosion resistance are of high importance. However, further research is required to explore the applicability and feasibility of the burnishing process in the case of complicated parts consisting of free-form surfaces.
Footnotes
Acknowledgements
The authors would like to acknowledge Meda Hemalatha for her support during fuzzy logic modelling
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
