Next Article in Journal
Porous Carbon–Carbon Composite Materials Obtained by Alkaline Dehydrochlorination of Polyvinyl Chloride
Previous Article in Journal
Comparative Study of Carbon Nanotube Composites as Capacitive and Piezoresistive Pressure Sensors under Varying Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Minimization of Surface Roughness and Temperature during Turning of Aluminum 6061 Using Response Surface Methodology and Desirability Function Analysis

by
Endalkachew Mosisa Gutema
1,
Mahesh Gopal
1 and
Hirpa G. Lemu
2,*
1
Department of Mechanical Engineering, College of Engineering and Technology, Wollega University, Nekemte P.O. Box 395, Ethiopia
2
Faculty of Science and Technology, University of Stavanger, N-4036 Stavanger, Norway
*
Author to whom correspondence should be addressed.
Submission received: 20 September 2022 / Revised: 27 October 2022 / Accepted: 27 October 2022 / Published: 30 October 2022

Abstract

:
Aluminum alloy is the second most abundant metal on Earth, known for its wide range of utilization in commercial goods due to its heat capacity and tensile strength. This study examines the effect of nose radius on the turning process. Further, it explores the implications of cutting parameters such as the cutting speed, the rate of feed, the cutting depth, and the nose radius of the tool. The trials were carried out with an Al 6061 workpiece and an Al2O3-coated carbide tool as the cutter, utilizing the response surface methodology. A mathematical model was developed to investigate the performance characteristics of the turning operation using the analysis of variance method. The multi-response desirability function analysis combines individual desirability values to create a composite desirability value. The ideal parameter levels were determined using the composite desirability value, and the significant influence of parameters was assessed. The obtained optimum surface roughness and temperature parameters are at a cutting speed of 116.37 m/min, a rate of feed of 0.408 mm/rev, a cutting depth of 0.538 mm, and a tool nose radius of 0.20 mm. The related ideal surface roughness and temperature values are 0.374 µm and 27.439 °C. The optimal overall desirability value is 0.829, close to the target response.

1. Introduction

The ability to produce lightweight materials has created the path for the automotive and aerospace industries to reduce, for instance, the overall weight of vehicles [1]. To achieve this weight reduction, aluminum materials are used in various engineering fields to replace iron and steel, due to their superior mechanical properties including a high stiffness-to-weight ratio, physical and thermal qualities, corrosion resistance, and recyclability [2]. Aluminum and its alloys, on the other hand, have poor machinability. As a result, extensive research into studying performance attributes and efficient machining technologies are underway to lower the total costs.
To assess flank wear of the cutting tool in machining [3], experiments were conducted to investigate the cutting performance of aluminum alloy LM 25, which was reinforced with green-bonded silicon carbide. The study employed the response surface methodology (RSM) and desirability function analysis (DFA) to improve the surface roughness. A drilling process in aluminum-silicon carbide was performed [4] using the analysis of variance (ANOVA) and Taguchi methodologies to increase material removal rate (MRR) and to decrease surface roughness, where the RSM was employed to assess the tool’s surface waviness, thrust force, burr height, and wear. Rajmohan and Palanikumar [5] used Al 356 aluminum bonded with silicon carbide composition in the drilling experiment. In the study reported in [6], the composite desirability (CD) value of numerous performance attributes, such as surface quality and cutting force, was transformed into a single performance feature. Taguchi and DFA were used to optimize the parameters of aluminum-silicon carbide. The research work in [7] reported that the RSM and DFA enhanced the machining settings to maximize tool life and minimize energy usage during the Al 7075 alloy machining. The study used ANOVA to investigate the aluminum alloy, while the RSM was utilized to develop the mathematical model. The cutting speed was the most critical parameter that influenced the experiment compared to other parameters; further, sensitivity analysis was carried out to study the cutting force [8]. Mugendiran et al. [9] improved an aluminum alloy’s surface quality and wall thickness using the RSM and ANOVA. The RSM and the genetic algorithm were employed for the experimentation to optimize the aluminum sample. The input parameters considered were rotation speed, the rate of feed, the axial cutting and radial cutting depth, and tool rake angle [10]. Furthermore, the author in [11] optimized surface quality using the RSM and DFA of aluminum alloys. The RSM and DFA were also used to maximize aluminum alloy energy utilization, surface waviness, and material removal rate (MRR) by taking the cutting depth, the rate of feed, and the cutting speed into account [1,12,13].
The RSM is frequently used in conjunction with desirability functions to forecast more optimum results [14,15,16,17]. The research reported in [18] constructed an empirical model to anticipate overcut in machining using the Box Behnken design and compared the empirical model with the RSM model. The model corresponds well with the experimental data, and the model was compared with the RSM. DFA is a prominent method used in industry and academia to optimize several responses [19]. The author in [20] introduced the RSM and DFA as an optimizing tool to select the best option in turning replies into desirability indices ranging from 0 to 1. The response of predetermined maximum or minimum values was identified with process parameters inside a given range. The response was then transformed into the desired value via one-sided transformation in the desirability-based method. A numerical optimization method and desirability analysis were used to optimize the wear rate [21]. The author [22] employed the Taguchi method to investigate and determine the influencing variables using ANOVA. DFA was used as the performance index of the output by computing the composite DFA values. The research reported in [23] carried out three optimizations: (1) quality optimization, (2) economic optimization and (3) combined optimization, using ANOVA, the RSM, and DFA to minimize surface unevenness and wear in the tool. The desirability function combined with the Taguchi approach allows optimization of the multi-response concerns such as surface irregularity and MRR [24], optimization of the coating breadth, thickness, and interface temperature [25] and optimizing machining pressure and surface roughness [26]. An experimental study was conducted by [27] to optimize the microhardness and surface quality of IN625 material using Taguchi L9 optimization and Super Ranking techniques. The input parameters were laser power, scan speed, and hatch distance. The authors in [28] used the ANOVA, the RSM, fuzzy MCDM, fuzzy AHP, and fuzzy TOPSIS to optimize titanium (Ti6Al4V) alloy using a electrical wire discharge machining process to optimize the cutting speed, MRR, and surface roughness. The result emphasizes that the RSM is an effective tool for the design of experiments. The author [29] experimented by considering layer thickness, building orientation, raster angle, raster width, and air gap as input parameters and the output parameter as impact strength, flexural strength and tensile strength during processing by using fused-filament fabrication methodology in AM technology. The naked mole-rat algorithm (NMRA) was applied to solve the optimization process. The ball end milling operation was performed by [30] using hardened 55NiCrMoV6 steel material considering the cutting speed and surface inclination angle as input parameter to optimize the cutting force by using the RSM. The results indicated that the surface inclination angle has a significant influencing parameter compared to others. Aluminum 6061 is a precipitation-hardened alloy with excellent properties such as high ultimate tensile strength and yield strength that are desirable for many applications. It also contains significant elements of magnesium and silicon. Therefore, careful consideration should be taken during the machining process. The literature indicates that the desirability function technique is a proper method to enhance the multi-response optimization processes, such as those experienced in machining aluminum and its alloys [31,32].
Utilizing ANOVA and RSM, a lot of research have been reported on aluminum 6061 as the workpiece material. However, the application of ANOVA, the RSM, and DFA to optimize the input and output parameters during the turning of aluminum 6061 is not well documented in the literature. DFA application, in particular, has significant benefits to optimize and reduce several criteria simultaneously and to bring out the process changes that lead to better product quality and higher productivity.
In this study, experiments were conducted using the design of experiment (DoE) strategy and optimized using the RSM. The DoE is necessary to study the relationship between the multiple input and output variables and gain knowledge to estimate the best-operating conditions. Furthermore, the DoE was used to determine the individual and interactive effects of variables influencing the output result during measurements. In addition, DFA was used to convert the multiple response characteristics into single response characteristics and calculate the ideal machining environment of the process parameters to reduce surface roughness and temperature.

2. Methods, Materials and Process Parameters

2.1. Materials and Experimentation Setup

In this experiment, untreated aluminum 6061 with a diameter of 50 mm and a length of 100 mm was used as the working material. Table 1 shows the chemical properties of aluminum 6161 material.
The trials were performed in dry conditions on an XLTURN-CNC lathe (MTAB, Tamilnadu, India), Figure 1a, with an Al2O3-coated carbide tool by using a Sandvik Coromant T-Max P Turning Tool Holder (Sandviken, Sweden). The machining parameters considered were the cutting speed, the rate of feed, the cutting depth, and the nose radius of the tool. A 1 mm hole was bored into the aluminum material, and the sample was positioned 10 mm beneath the cutting surface. As illustrated in Figure 1b, a K-type thermocouple (Shanghai MKYD Instrument, Shanghai, China) was used to detect the temperature, and a SURFTEST SJ-201 surface roughness tester (Mitutoyo America Corporation, Aurora, IL, USA) was used to assess the roughness of the surface.

2.2. Process Parameter Identification for Experiment

For the experimental investigation, the following process factors have an impact on surface roughness and temperature were identified as per the recommendations of machining conditions and input parameters in [33]. The output responses were formulated as follows.
-
Minimization of the surface roughness (Ra in μm); the cutting factor equation is:
Minimize: Ra = c (Vck1 Fz k2 Dc k3 Rn k4)
-
Minimization of the temperature (Temp, °C) using the cutting parameters
Minimize: Temp = c (Vck1 Fz k2 Dc k3 Rn k4)
where k1, k2, k3, and k4 are the model parameters (estimated from experimental data) and c is the response error.

2.3. The Creation of a Design Matrix and the Selection of Parameter Levels

The DoE was used in the creation of a design matrix, selection of parameter levels and the experiments, while the operating limits of all evaluation criteria were used to determine the levels of parameters. As the DoE is a systematic approach used to solve engineering problems, its application reduces the number of experiments and gives a high level of control. The operating ranges of all parameters were determined through trial runs, with one parameter modified while the others remained fixed. Trial runs before the stated parameters were used to identify the upper-bound (+2) and lower-bound (2) levels of all five independent variables, resulting in the precise predicted values shown in Table 2 and Table 3. The intermediate levels of all variables, 1, 0, and +1, were calculated via interpolation [34]. The design matrix chosen to conduct the experiments using the DoE method was a four-factor central composite rotatable design (CCD) consisting of 30 sets of coded conditions and a full replication, where 24 are non-center points and 6 are center points.

2.4. Surface Roughness Prediction Using a Response Surface Model

The experimental results were examined methodically using the Design-Expert software V11 (StatEase, Minneapolis, MN, USA). A second-order quadratic model is developed to predict surface roughness. As reported in [35], ANOVA has been used to determine the model’s suitability. The ANOVA table for Ra prediction is shown in Table 4.
The F-value of 19.96 suggests that the model is significant, where this large F-value might occur due to noise of 0.01% of the time. p-values less than 0.0500 indicate that the model terms, i.e., the model terms Vc, Dc, Vc, Fz, and Vc2 are significant. On the other hand, values over 0.1000 suggest that the model terms are not significant. Model reduction may help for a model with many insignificant terms (except those necessary to enable hierarchy).
The F-value of 0.2783 for the lack of fit suggests that it is minimal compared to the pure error. A substantial lack of fit F-value may occur 95.96% of the time due to noise, and a non-significant lack of fit is preferred. The following Design-Expert software’s (StatEase, Minneapolis, MN, USA) regression equation of the objective factors was obtained upon running the regression.
Ra = +0.786250 − 0.006694 × Vc + 0.476852 × Fz + 0.089583 × Dc − 0.110417 × Rn
− 0.007407 × Vc × Fz − 3.69886 × 10−17× Vc × Dc − 0.000417 × Vc × Rn v
0.208333 × Fz Dc + 0.277778 × Fz × Rn − 1.50104 × 10−15 × Dc Rn + 0.000033Vc2 +
0.462963 × Fz2 − 7.93479 × 10−16 Dc2 + 0.062500 × Rn2

2.5. Temperature Prediction Using a Response Surface Model

A second-order quadratic model was developed to predict temperature. ANOVA analysis was used to determine the model’s fitness. The ANOVA table for the temperature prediction is shown in Table 5.
The F-value of 786.09 for the model indicates that it is significant. An F-value this large may be attributable to noise merely 0.01 percent of the time. Significant model terms have p-values less than 0.0500. The following model terms are significant: Vc, Fz, Dc, Rn, Vc Dc, Vc Rn, Fz Dc, Fz Rn, Dc Rn, Vc2, Fz2, Dc2, and Rn2. Values greater than 0.1000 indicate that the model terms are insignificant.
As can be observed from the table, the F-value for lack of fit is 0.2783, which indicates that the lack of fit is negligible compared to the pure error. A significant lack of fit F-value due to noise has a 96.06% chance of occurring, but a non-significant lack of fit was preferred. The regression equation for the Design-Expert software (StatEase, Minneapolis, MN, USA) in terms of actual variables is provided below.
Temp = +123.32917 − 1.86083 × Vc − 19.67593 × Fz + 29.35417 × Dc − 43.10417
× Rn − 0.023148 × Vc × Fz − 0.260417 × Vc × Dc + 0.431250 × Vc × Rn − 23.26389
× Fz × Dc + 23.95833 × Fz × Rn − 28.59375 × Dc × Rn + 0.009097 × Vc2 + 44.36728
× Fz2 + 14.60937 × Dc2 + 7.42187 × Rn2

3. Discussion of Results

3.1. Interaction Effect of Surface Roughness

Figure 2a depicts the interaction plot between the cutting speed (Vc) and the rate of feed (Fz) in terms of surface roughness (µm). The cutting speed increases consistently, the surface waviness decreases, and the rate of feed increases; the surface roughness increases correspondingly. The cutting speed is between 90 and 105 m/min, the surface waviness value is low, and the rate of feed is low, at 0.09 mm/rev. The result shows that the surface waviness is negligible at a medium cutting speed and a low rate of feed. On the other hand, changes in the feed rate significantly impact both low and high cutting speeds. The cutting speed and the rate of feed determine surface finish and MRR. The ANOVA was also used to validate the result (Table 4).
The interaction impact of the cutting speed and the cutting depth on surface roughness is shown in Figure 2b. The cutting depth is low at a higher cutting speed substantially impacts surface roughness. The cutting depth is a critical cutting parameter that influences the whole process, including stability, cutting forces, vibrations, and spindle load. The surface roughness is lower at a higher cutting speed and a lower depth of cut. The effect of the cutting speed and tool nose radius on surface roughness is shown in Figure 2c. Lower cutting speeds result in higher surface roughness, while lower tool nose radius results in significantly lower surface roughness. The surface roughness significantly increases between the tool nose radius of 0.2 mm and 0.4 mm. At all cutting speeds, the same pattern emerges. The results revealed that the industry deserves a suitable surface roughness; the nose radius should be 0.8 mm and 1 mm.
Figure 2d depicts a graph of observed values vs. predicted values, and it assists in identifying the observations. The 45° line evenly divides the data points.

3.2. Interaction Effect of Temperature

The temperature interaction plot between the cutting speed and the rate of feed is depicted in Figure 3a using the interaction effect. When the cutting speed is between 90 and 105 m/min, the temperature is relatively low, and it is also inadequate when the rate of feed is 0.27 to 0.36 mm/rev. At a medium cutting speed and rate of feed, the temperature is low. The ANOVA may also be used to validate the result (Table 5).
Figure 3b shows the effect of temperature on the cutting speed and the cutting depth. A smaller cutting depth significantly influences the temperature reduction. The cutting force increases when there is an increase in the cutting depth and the rate of feed, reducing with an increasing cutting speed.
Figure 3c demonstrates the effect of temperature interaction on the cutting speed and tool nose radius. The temperature rises dramatically with increasing the cutting speed and falls significantly with increasing tool nose radius, with the temperature decreasing noticeably between 0.6 and 1 mm.
The result indicates that to maintain a low temperature in the work piece, the nose radius should be between 0.8 and 1 mm. Figure 3d depicts a graph of observed vs. predicted values, which aids in detecting observations that the model fails to predict. The 45° line should split the data points evenly.

3.3. Desirability Function Approach

Researchers have employed various optimization approaches to increase product quality and productivity [7,31,32]. The studies indicate that a desirability function analysis improves the process parameters. The following procedures were used to optimize the input parameters using DFA and the RSM.
The work reported in [20] suggested the individual desirability index (di) for the comparable responses, where two desirability functions were considered based on the response characteristics, namely (a) the smaller, the better and (b) the larger, the better. The experiment was carried out to determine the features of electro-discharge machining (EDM) using an optimization-based desirability technique [36]. The experiment analysis was performed on hardened 55NiCrMoV6 steel to minimize cutting force and vibration. The signal-to-noise ratio (S/N) and grey relational analysis optimization techniques were performed, considering surface inclination angle and tool overhang as input parameters. According to the optimization results, changes in tool overhang increase the cutting force and vibration [37].
(a) The smaller the better: This method aims to reduce the output variable to the smallest possible value; hence, the smaller, the better characteristic is utilized to establish distinct desirability levels. To normalize their values and evaluate individual desirability indexes, the smaller, the better desirability function is utilized within the range of [0, 1]. When the response must be lowered, according to Harrington [38], Equation (5) defines the smaller, the better sort of quality characteristic.
The desirability functions di can be expressed as.
di = 1 0   y m   y max y tgt   y max r  
where ym is minimized and
y m   <     y tar ,     y tar   <   y m   <   y max ,   r     0 ,   and   y m   <     y max
In this case, ytar is the response’s lowest value ym, i.e., ytar = ymin.
Equation (5) can be modified as follows (Equation (6)).
di = 1 0   y m   y max y min   y max r
where y m < y min , y min < y m < y max , r ≥ 0 and y m < y max
The condition is most attractive when ym is smaller than ymin and the individual desire score is 1. When ym deviates from ymin, the value of di declines until it approaches “0” when yi exceeds ymax. In all other cases, the di values produced are in the range [0, 1]. The exponent w represents the weights allocated to each response depending on their relevance. Equation (7) computes the composite desirability as the geometric mean of each experimental circumstance’s desires.
(b) Composite Desirability (CD):
The composite desirability is expressed as:
CD = D 1 W 1 ×   D 2 W 2 ×   D 3   W 3   ×     1 N
where N is the number of responses, D1, D2, and D3 are individual desirability indexes, and W1, W2, and W3 and weight assigned responses, where w = i N W i   = 1.
If any of the replies is entirely unacceptable, the value of CD becomes zero, i.e., Di = 0.
(c) Choosing the optimum parameter combination and level: A developed composite desirability value suggests that the product is of quality. As a result, the composite desirability is used to evaluate the parameter result and ideal level for each manageable parameter.

3.3.1. Desirability Function Approach for Surface Roughness

The desirability functions presented in Equation (6) are used to calculate the individual desirability of surface unevenness. The minimum surface unevenness value is 0.32 µm, and the maximum value is 0.49 µm. Surface roughness increase is given equal the weightage (w1 = w2 = 1).
d temp = 1 0   y m   0.49 0.32 0.49 1
where y m < 0.32, 0.32 < y m   < 0.49, r ≥ 0 and y m < 0.49.
Equation (9) was employed to compute the surface roughness based on the responses under consideration.
CD = D µ m i   1 2 where   i   =   1   to   30

3.3.2. Desirability Function Approach for Temperature

Individual desirability of temperature is calculated using the desirability functions as given in Equation (6). Surface roughness increases by giving equal weightage (w1 = w2 = 1), the minimum temperature value is 24.2°, and the maximum temperature value is 38.8°.
d temp = 1 0   y m 38.8 24.2 38.8 1
where y m < 24.2, 24.2 < y m   < 38.8, r ≥ 0 and y m < 38.8.
Equation (9) computes the surface roughness based on the responses under consideration.
CD = D temp 1 1 2 where   i   =   1   to   30

3.3.3. Process Parameter Multi-Response Optimization for Lowest Surface Imperfection and Temperature

Świercz et al. [39] used a multi-response optimization (MRO) technique for DFA. The optimization process was used, and the use of a statistical alternative to specifying the test to be performed earlier. Nowadays, the desirability function is extensively used to reduce multi-response situations into a single response [40]. The author [41] used the RSM and DFA to optimize the results. The highest overall desirability factor settings are found under the optimal parameter situations. The simultaneous goal function is a geometric mean of all updated replies. To minimize surface roughness and temperature, a multi-response optimization analysis of process parameters (Table 6) was performed using Design Expert V11 (StatEase, Minneapolis, MN, USA). Models were used to optimize surface roughness and temperature. The measure of the solution to fulfill the specified objectives for all responses are reviewed in MRO.
The objectives of the two optimization approaches, as shown in Table 6 and Figure 4a,b, include criteria for reducing surface roughness and temperature, respectively. For instantaneous optimization, each response must have a low and high value given to each objective. The optimization is performed by a set of objectives, which have been applied to the variables and replies. The goal of the replies was to “minimize”. A weight can be used for a purpose to affect the shape of its desirability function. The results remain organized in order of preference, with the most desirable appearing first.

3.4. Desirability Function Analysis

To achieve low surface roughness and temperature, desirability function analysis optimization was used. The two responses have the same goal, i.e., ‘minimize’ surface roughness and temperature. The depreciation seeks to get high-quality finished machined products with lower temperatures for more excellent performance and cost savings. These investigations targeted surface roughness and temperature characteristics are ‘smaller-is-better’. The desirability-based approach finds many best solutions, and the option with the most significant desirability is preferred: the best 63 solutions were acquired in the optimization.
Table 7 shows that only five higher desirability alternatives were chosen for the investigation, with 1 having tremendous desirability (D = 0.829). The optimal turning conditions are a cutting speed of 116.37 m/min, a rate of feed of 0.408 mm/rev, a cutting depth of 0.538 mm, and a tool nose radius of 0.20 mm. The ideal surface roughness and temperature are 0.374 µm and 27.439°. The smaller, the better principle was selected.
The estimated surface over desirability for Vc vs. Fz, Vc vs. Dc and Vc vs. Rn are shown in Figure 5a–c, while Figure 6a–c shown the contour plots on the desirability, respectively. The surface plots on the desirability, shown in Figure 5, are two-dimensional representations of a three-dimensional connection, where variables are given on the x- and y-axes and a smooth surface defining the dependent variables is shown on the z-axis. The contour plots on desirability are also shown in Figure 6a–c. A contour plot diagram is employed to investigate the relationship between the three variables. Two independent variables are displayed on the x-axis and y-axis, while one dependent variable is displayed on the z-axis in these plots. Contour plots aid in the identification of combinations that provide satisfactory results. These statistics also show that the best response was obtained using a higher cutting speed, a lower rate of feed, and a low cutting depth.
Figure 7 depicts the ramps of the optimized results. The red dots represent the precise values of the components, while the blue dots represent how well the objectives are satisfied. It is preferable if the ramp is higher up. The relevant desirability bar graph for the provided cutting circumstances, replies, and the combined desirability of 0.829 is shown in Figure 8. The blue bottom bar represents the overall desirability of all the parameters and responses. However, it is essential to note that the objective of an optimization process is to identify an appropriate combination of conditions that meet the goals, not to get a D value of 1. As a result, attractiveness is best accomplished at a low rate of feed and quick speeds while keeping a consistent low cutting depth. Additionally, the ideal solutions and contour map give first-hand data for investigators and industrialists to select optimum machining settings based on design or customer demands. Table 8 shows the optimum values during the turning of aluminum.

4. Model Validation

Desirability function analysis is used to estimate the optimal conditions and then validated by physical measurements. The model’s validity is demonstrated using the error percentage of less than ±2%. Surface roughness and temperature testing findings indicate the model is in good agreement with the optimal cutting settings (as predicted by DSA), as shown in Table 9.

5. Conclusions

The design of experiment was developed utilizing the central composite design of the RSM to predict surface roughness and temperature rise. Further, desirability function analysis optimization was performed to find the optimum value. The following conclusions have been drawn from this study:
  • The cutting speed is the most important influencing parameter compared to the other parameters.
  • The surface roughness is minimum at a cutting speed range of 90 m/min to 105 m/min, a rate of feed of 0.09 mm/rev and 0.6 mm and a tool nose radius of 1 mm.
  • The temperature is low at a 90 105 m/min cutting speed, and it is much lower when the rate of feed is 0.27 to 0.36 mm/rev. It is observed to be even lower when the nose radius is between 0.6 and 1 mm.
  • The desired function optimization strategy is suggested to obtain the optimal tuning parameters. The low surface roughness and temperature value are achievable at a cutting speed of 116.37 m/min, a rate of feed of 0.408 mm/rev, a cutting depth of 0.538 mm, and a tool nose radius of 0.200 mm The related ideal surface roughness and temperature values are 0.374 µm and 27.439°, respectively, with a desirability of 0.829 for the value of 1.

Author Contributions

Conceptualization, E.M.G. and M.G.; methodology, E.M.G. and M.G.; software, M.G.; validation, E.M.G. and H.G.L.; formal analysis, M.G. and E.M.G.; investigation, E.M.G. and M.G.; resources, E.M.G. and H.G.L.; data curation, E.M.G.; writing—original draft preparation, M.G.; writing—review and editing, E.M.G. and H.G.L.; visualization, E.M.G. and H.G.L.; supervision, E.M.G. and H.G.L.; project administration, H.G.L.; funding acquisition, E.M.G. and H.G.L. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by INDMET project granted by NORHED II Program (Grant Nr. 62862).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available upon request from the authors.

Acknowledgments

This research performed through the INDMET Project granted by NORHED II Program (Grant Nr. 62862). The authors will gratefully acknowledge the financial support provided for the Publication.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nataraj, M.; Balasubramanian, K. Parametric optimization of CNC turning process for hybrid metal matrix composite. Int. J. Adv. Manuf. Technol. 2017, 93, 215–224. [Google Scholar] [CrossRef]
  2. Miller, W.S.; Zhuang, L.; Bottema, J.; Wittebrood, A.J.; De Smet, P.; Haszler, A.; Vieregge, A. Recent development in aluminium alloys for the automotive industry. Mater. Sci. Eng. A 2000, 280, 37–49. [Google Scholar] [CrossRef]
  3. Seeman, M.; Ganesan, G.; Karthikeyan, R.; Velayudham, A. Study on tool wear and surface roughness in machining of particulate aluminum metal matrix composite-response surface methodology approach. Int. J. Adv. Manuf. Technol. 2010, 48, 613–624. [Google Scholar] [CrossRef]
  4. Dhavamani, C.; Alwarsamy, T. Optimization of Machining Parameters for Aluminum and Silicon Carbide Composite Using Genetic Algorithm. Procedia Eng. 2012, 38, 1994–2004. [Google Scholar] [CrossRef] [Green Version]
  5. Rajmohan, T.; Palanikumar, K. Application of the central composite design in optimization of machining parameters in drilling hybrid metal matrix composites. Measurement 2013, 46, 1470–1481. [Google Scholar] [CrossRef]
  6. Muthukrishnan, N.; Babu, T.S.M.; Ramanujam, R. Fabrication and turning of Al/SiC/B4C hybrid metal matrix composites optimization using desirability analysis. J. Chin. Inst. Ind. Eng. 2012, 29, 515–525. [Google Scholar] [CrossRef]
  7. Bhushan, R.K. Optimization of cutting parameters for minimizing power consumption and maximizing tool life during machining of Al alloy SiC particle composites. J. Clean. Prod. 2013, 39, 242–254. [Google Scholar] [CrossRef]
  8. Joardar, H.; Das, N.S.; Sutradhar, G.; Singh, S. Application of response surface methodology for determining cutting force model in turning of LM6/SiCP metal matrix composite. Measurement 2014, 47, 452–464. [Google Scholar] [CrossRef]
  9. Mugendiran, V.; Gnanavelbabu, A.; Ramadoss, R. Parameter Optimization for Surface Roughness and Wall Thickness on AA5052 Aluminium Alloy by Incremental Forming Using Response Surface Methodology. Procedia Eng. 2014, 97, 1991–2000. [Google Scholar] [CrossRef] [Green Version]
  10. Mahesh, G.; Muthu, S.; Devadasan, S.R. Prediction of surface roughness of end milling operation using genetic algorithm. Int. J. Adv. Manuf. Technol. 2015, 77, 369–381. [Google Scholar] [CrossRef]
  11. Horváth, R.; Drégelyi-Kiss, Á. Analysis of surface roughness of aluminum alloys fine turned: United phenomenological models and multi-performance optimization. Measurement 2015, 65, 181–192. [Google Scholar] [CrossRef]
  12. Camposeco-Negrete, C. Optimization of cutting parameters using Response Surface Method for minimizing energy consumption and maximizing cutting quality in turning of AISI 6061 T6 aluminum. J. Clean. Prod. 2015, 91, 109–117. [Google Scholar] [CrossRef]
  13. Kumar, V.; Kharub, M.; Sinha, A. Modeling and Optimization of Turning Parameters during Machining of AA6061 composite using RSM Box-Behnken Design. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1057, 012058. [Google Scholar] [CrossRef]
  14. Senthil, S.M.; Parameshwaran, R.; Ragu Nathan, S.; Bhuvanesh Kumar, M.; Deepandurai, K. A multi-objective optimization of the friction stir welding process using RSM-based-desirability function approach for joining aluminum alloy 6063-T6 pipes. Struct. Multidiscip. Optim. 2020, 62, 1117–1133. [Google Scholar] [CrossRef]
  15. Srivastava, M.; Maheshwari, S.; Kundra, T.; Rathee, S. Multi-Response Optimization of Fused Deposition Modelling Process Parameters of ABS Using Response Surface Methodology (RSM)-Based Desirability Analysis. Mater. Today Proc. 2017, 4, 1972–1977. [Google Scholar] [CrossRef]
  16. Verma, S.; Gupta, M.; Misra, J.P. Optimization of process parameters in friction stir welding of armor-marine grade aluminium alloy using desirability approach. Mater. Res. Express 2018, 6, 026505. [Google Scholar] [CrossRef]
  17. Pandiyarajan, R.; Maran, P.; Murugan, N.; Marimuthu, S.; Sornakumar, T. Friction stir welding of hybrid AA 6061-ZrO2-C composites FSW process optimization using desirability approach. Mater. Res. Express 2019, 6, 066553. [Google Scholar] [CrossRef]
  18. Rouniyar, A.K.; Shandilya, P. Semi-empirical modeling and optimization of process parameters on overcut during MFAPM-EDM of Al6061 alloy. Proc. Inst. Mech. Eng. Part E J. Process. Mech. Eng. 2021, 235, 1784–1796. [Google Scholar] [CrossRef]
  19. Devarajaiah, D.; Muthumari, C. Evaluation of power consumption and MRR in WEDM of Ti–6Al–4V alloy and its simultaneous optimization for sustainable production. J. Braz. Soc. Mech. Sci. Eng. 2018, 40, 400. [Google Scholar] [CrossRef]
  20. Derringer, G.; Suich, R. Simultaneous Optimization of Several Response Variables. J. Qual. Technol. 1980, 12, 214–219. [Google Scholar] [CrossRef]
  21. Canute, X.; Majumder, M.C. Investigation of tribological and mechanical properties of aluminium boron carbide composites using response surface methodology and desirability analysis. Ind. Lubr. Tribol. 2018, 70, 301–315. [Google Scholar] [CrossRef]
  22. Pradhan, S.; Singh, G.; Bhagi, L.K. Study on surface roughness in machining of Al/SiCp metal matrix composite using desirability function analysis approach. Mater. Today Proc. 2018, 5, 28108–28116. [Google Scholar] [CrossRef]
  23. Jumare, A.I.; Abou-El-Hossein, K.; Abdulkadir, L.N.; Liman, M.M. Predictive modeling and multiobjective optimization of diamond turning process of single-crystal silicon using RSM and desirability function approach. Int. J. Adv. Manuf. Technol. 2019, 103, 4205–4220. [Google Scholar] [CrossRef]
  24. Kumar, P.; Misra, J.P. Optimization of duplex stainless steel dry turning parameters using desirability function. Mater. Today Proc. 2020, 26, 1580–1584. [Google Scholar] [CrossRef]
  25. Sahoo, D.K.; Mohanty, B.S.; John, D.F.; Pradeep, A.M.V. Multi response optimization and desirability function analysis on friction surfaced deposition of AISI 316 stainless steel over EN8 medium carbon steel. Mater. Today Proc. 2020, 40, S1–S9. [Google Scholar] [CrossRef]
  26. Sarath, V.V.N.; Tamiloli, N. Optimization of multiple quality characteristics for end milling under dry cutting environment using desirability function. J. Phys. Conf. Ser. 2021, 2070, 012218. [Google Scholar] [CrossRef]
  27. Sheshadri, R.; Nagaraj, M.; Lakshmikanthan, A.; Chandrashekarappa, M.P.G.; Pimenov, D.Y.; Giasin, K.; Prasad, R.V.S.; Wojciechowski, S. Experimental investigation of selective laser melting parameters for higher surface quality and microhardness properties: Taguchi and super ranking concept approaches. J. Mater. Res. Technol. 2021, 14, 2586–2600. [Google Scholar] [CrossRef]
  28. Fuse, K.; Dalsaniya, A.; Modi, D.; Vora, J.; Pimenov, D.Y.; Giasin, K.; Prajapati, P.; Chaudhari, R.; Wojciechowski, S. Integration of Fuzzy AHP and Fuzzy TOPSIS Methods for Wire Electric Discharge Machining of Titanium (Ti6Al4V) Alloy Using RSM. Materials 2021, 14, 7408. [Google Scholar] [CrossRef]
  29. Chohan, J.S.; Mittal, N.; Kumar, R.; Singh, S.; Sharma, S.; Dwivedi, S.; Saxena, A.; Chattopadhyaya, S.; Ilyas, R.; Le, C.; et al. Optimization of FFF Process Parameters by Naked Mole-Rat Algorithms with Enhanced Exploration and Exploitation Capabilities. Polymers 2021, 13, 1702. [Google Scholar] [CrossRef]
  30. Wojciechowski, S.; Maruda, R.; Barrans, S.; Nieslony, P.; Krolczyk, G.M. Optimisation of machining parameters during ball end milling of hardened steel with various surface inclinations. Measurement 2017, 111, 18–28. [Google Scholar] [CrossRef]
  31. Azlan, A.; Mohd, A.L.; Nur, K.Y.; Syaiful, N.A.R. Statistical Optimization by the Response Surface Methodology of Direct Recycled Aluminum-Alumina Metal Matrix Composite (MMC-AlR) Employing the Metal Forming Process. Processes 2020, 8, 805. [Google Scholar] [CrossRef]
  32. Qazi, M.; Abas, M.I.; Khan, R.; Saleem, W.; Pruncu, C.I.; Omair, M. Experimental Investigation and Multi-Response Optimization of Machinability of AA5005H34 Using Composite Desirability Coupled with PCA. Metals 2021, 11, 235. [Google Scholar] [CrossRef]
  33. Hindustan Machine Tools (HMT). Production Technology; Tata McGraw-Hill Education: Noida, India, 2001. [Google Scholar]
  34. Box, G.E.P.; Hunter, J.S.; Hunter, W.G. Statistics for Experimenters: Design, Innovation, and Discovery; Wiley & Sons: New York, NY, USA, 2005. [Google Scholar]
  35. Myers, R.H.; Montgomery, D.C.; Anderson-Cook, C.M. Response Surface Methodology: Process and Product Optimization Using Designed Experiments; John Wiley & Sons: New York, NY, USA, 2016. [Google Scholar]
  36. Gopalakannan, S.; Senthilvelan, T. Optimization of machining parameters for EDM operations based on central composite design and desirability approach. J. Mech. Sci. Technol. 2014, 28, 1045–1053. [Google Scholar] [CrossRef]
  37. Wojciechowski, S.; Maruda, R.W.; Krolczyk, G.M.; Niesłony, P. Application of signal to noise ratio and grey relational analysis to minimize forces and vibrations during precise ball end milling. Precis. Eng. 2018, 51, 582–596. [Google Scholar] [CrossRef]
  38. Harrington, E.C. The Desirability Function. Ind. Qual. Control 1965, 21, 494–498. [Google Scholar]
  39. Świercz, R.; Oniszczuk-Świercz, D.; Chmielewski, T. Multi-Response Optimization of Electrical Discharge Machining Using the Desirability Function. Micromachines 2019, 10, 72. [Google Scholar] [CrossRef] [Green Version]
  40. Padilla-Atondo, J.M.; Limon-Romero, J.; Perez-Sanchez, A.; Tlapa, D.; Baez-Lopez, Y.; Puente, C.; Ontiveros, S. The Impact of Hydrogen on a Stationary Gasoline-Based Engine through Multi-Response Optimization: A Desirability Function Approach. Sustainability 2021, 13, 1385. [Google Scholar] [CrossRef]
  41. Chabbi, A.; Yallese, M.A.; Meddour, I.; Nouioua, M.; Mabrouki, T.; Girardin, F. Predictive modeling and multi-response optimization of technological parameters in turning of Polyoxymethylene polymer (POM C) using RSM and desirability function. Measurement 2017, 95, 99–115. [Google Scholar] [CrossRef]
Figure 1. Experimental setup.
Figure 1. Experimental setup.
Materials 15 07638 g001
Figure 2. Interaction effect in terms of surface roughness of (a) Vc vs. Fz, (b) Vc vs. Dc, (c) Vc vs. Rn and (d) plot of actual vs. predicted values. Plots produced using Design-Expert software V11 (StatEase, Minneapolis, MN, USA).
Figure 2. Interaction effect in terms of surface roughness of (a) Vc vs. Fz, (b) Vc vs. Dc, (c) Vc vs. Rn and (d) plot of actual vs. predicted values. Plots produced using Design-Expert software V11 (StatEase, Minneapolis, MN, USA).
Materials 15 07638 g002
Figure 3. Interaction effect of temperature (a) Vc vs. Fz, (b) Vc vs. Dc, (c) Vc vs. Rn and (d) actual vs. predicted plot. Plots produced using Design-Expert software V11 (StatEase, Minneapolis, MN, USA).
Figure 3. Interaction effect of temperature (a) Vc vs. Fz, (b) Vc vs. Dc, (c) Vc vs. Rn and (d) actual vs. predicted plot. Plots produced using Design-Expert software V11 (StatEase, Minneapolis, MN, USA).
Materials 15 07638 g003
Figure 4. Criteria to (a) minimize µm and (b) minimize temp (T°).
Figure 4. Criteria to (a) minimize µm and (b) minimize temp (T°).
Materials 15 07638 g004
Figure 5. Estimated surface plots on desirability (a) Vc vs. Fz, (b) Vc vs. Dc and (c) Vc vs. Rn. Plots produced using Design-Expert software V11 (StatEase, Minneapolis, MN, USA).
Figure 5. Estimated surface plots on desirability (a) Vc vs. Fz, (b) Vc vs. Dc and (c) Vc vs. Rn. Plots produced using Design-Expert software V11 (StatEase, Minneapolis, MN, USA).
Materials 15 07638 g005
Figure 6. Contour plot on desirability for (a) Vc vs. Fz, (b) Vc vs. Dc and (c) Vc vs. Rn. Plots produced using Design-Expert software V11 (StatEase, Minneapolis, MN, USA).
Figure 6. Contour plot on desirability for (a) Vc vs. Fz, (b) Vc vs. Dc and (c) Vc vs. Rn. Plots produced using Design-Expert software V11 (StatEase, Minneapolis, MN, USA).
Materials 15 07638 g006
Figure 7. Optimal parameter ramps function graphs and combined optimization.
Figure 7. Optimal parameter ramps function graphs and combined optimization.
Materials 15 07638 g007
Figure 8. Desirability bar graph for combined optimization.
Figure 8. Desirability bar graph for combined optimization.
Materials 15 07638 g008
Table 1. Chemical properties of aluminum 6061.
Table 1. Chemical properties of aluminum 6061.
Al 6061MgSiFeMnCuCrZnTiAl
Weight
(%)
0.8–
1.2
0.40–
0.80
0.0–0.700.150.15–
0.40
0.04–
0.35
0.0–0.250.0–0.15Bal
Table 2. The factorial levels of output process parameters.
Table 2. The factorial levels of output process parameters.
Sl. NoParametersFactorial Levels
−2−1012
1Cutting speed, Vc (m/min)75.090.0105.0120.0135.0
2Rate of feed, Fz (mm/rev)0.090.180.270.360.45
3Cutting depth, Dc (mm)0.200.400.600.801.00
4Tool nose radius, Rn (mm)0.200.400.600.801.00
Table 3. Responses to experimental values.
Table 3. Responses to experimental values.
Sl.
No
Cut.
Speed
(Vc)
Rate of Feed
(Fz)
Cut. Depth
(Dc)
Tool Nose Radius (Rn)Rough. (µmob)
(Practical)
Rough. (µm RSM)
(Pred Using RSM)
Temp. (°Cob)
Tob
(Practical)
Temp. (°CRSM)
(Pred Using RSM)
11200.360.80.40.380.372128.828.77
2900.180.80.40.440.447129.829.84
31050.271.00.60.410.411726.626.55
41050.270.60.60.390.398325.825.77
5900.360.80.80.460.458824.924.76
61050.270.60.20.410.410026.926.80
71200.360.80.80.370.377929.930.10
81050.270.60.60.390.398326.225.77
9750.270.60.60.490.486729.529.58
10900.180.40.40.430.426326.826.71
11900.360.80.40.450.447928.628.61
121050.270.60.60.420.398325.625.77
131050.270.20.60.390.385029.629.66
141200.180.80.80.400.397129.829.73
151200.180.40.40.390.390430.130.12
161050.270.60.60.390.398325.625.77
17900.360.40.80.450.452927.827.88
18900.360.40.40.440.442127.227.16
191050.270.61.00.410.406727.127.11
201200.180.40.80.370.376334.234.30
21900.180.40.80.410.417125.825.71
221050.450.60.60.410.411727.627.61
231050.090.60.60.420.415026.826.80
241200.360.40.40.360.366330.430.45
251050.270.60.60.390.398325.825.77
261200.180.80.40.410.411230.130.13
271350.270.60.60.370.370038.438.33
28900.180.80.80.440.437924.224.26
291200.360.40.80.380.372136.536.34
301050.270.60.60.410.398325.625.77
Table 4. Surface roughness—ANOVA table.
Table 4. Surface roughness—ANOVA table.
SourceSum of Squares ValuedfMean Square ValueF-Valuep-Value
Model0.0256140.001819.96<0.0001Significant
Vc0.021910.0219239.28<0.0001
Fz0.000010.00000.17450.6820
Dc0.001010.001011.170.0045
Rn9.000 × 10−619.000 × 10−60.09820.7583
Vc Fz0.001610.001617.450.0008
Vc Dc0.000010.00000.00001.0000
Vc Rn0.000010.00000.27270.6091
Fz Dc0.000210.00022.450.1380
Fz Rn0.000410.00044.360.0542
Dc Rn0.000010.00000.00001.0000
Vc20.001510.001516.830.0009
Fz20.000410.00044.210.0581
Dc20.000010.00000.00001.0000
Rn20.000210.00021.870.1916
Residual0.0014150.0001
Lack of Fit0.0005100.00000.27830.9596Not Significant
Pure Error0.000950.0002
Cor Total0.027029
Table 5. ANOVA table of temperature.
Table 5. ANOVA table of temperature.
SourceSum of Squares ValuedfMean Square ValueF-Valuep-Value
Model311.211422.23786.09<0.0001Significant
Vc49.56149.561752.53<0.0001
Fz1.0111.0135.72<0.0001
Dc9.7019.70343.14<0.0001
Rn0.429010.429015.170.0014
Vc Fz0.015610.01560.55260.4688
Vc Dc9.7719.77345.35<0.0001
Vc Rn26.78126.78947.06<0.0001
Fz Dc2.8112.8199.22<0.0001
Fz Rn2.9812.98105.23<0.0001
Dc Rn20.93120.93740.18<0.0001
Vc2114.921114.924063.88<0.0001
Fz23.5413.54125.27<0.0001
Dc29.3719.37331.24<0.0001
Rn22.4212.4285.49<0.0001
Residual0.4242150.0283
Lack of Fit0.1508100.01510.27590.9606not significant
Pure Error0.273350.0547
Cor Total311.6329
Table 6. Range of parameters and responses for desirability.
Table 6. Range of parameters and responses for desirability.
Sl.NoInput ParameterGoalLower Limit Upper Limit
1VcIn range75130
2FzIn range0.090.45
3DcIn range0.21
4RnIn range0.21
5Surface roughness observed (µm)Minimize0.320.49
6Temperature rise observed (Tob) DegreeMinimize24.238.8
Table 7. Best global solution for an optimization.
Table 7. Best global solution for an optimization.
NumberVcFzDcRnSurface RoughnessTemperatureDesirability
1116.370.410.540.200.3827.440.829Selected
2116.350.410.540.200.3827.430.829
3116.370.410.540.200.3727.430.829
4116.340.410.530.200.3827.430.829
5116.280.410.540.200.3827.430.828
Table 8. Optimum values during turning of aluminum.
Table 8. Optimum values during turning of aluminum.
Sl.NoInput ParameterGoalOptimum Value
1Vc (m/mm)In range116.371
2Fz (mm/rev)In range0.408
3Dc (mm)In range0.538
4Rn (mm)In range0.200
5Surface roughness (µm)Minimize0.374
6Temperature rise Observed (Tob) (°C)Minimize27.439
7Overall Desirability 0.829
Table 9. Optimized process parameter—validation model.
Table 9. Optimized process parameter—validation model.
Sl. NoVcFzDcRnObservationsConfirmatory Test% Error
Opt. Value by DFAExper. Value
1116.370.410.540.20Surface roughness (µm)0.3740.37001.06
Temperature (°C)27.43927.43100.02
2114.570.380.490.200Surface roughness (µm)0.3790.37301.58
Temperature (°C)27.40627.39800.03
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gutema, E.M.; Gopal, M.; Lemu, H.G. Minimization of Surface Roughness and Temperature during Turning of Aluminum 6061 Using Response Surface Methodology and Desirability Function Analysis. Materials 2022, 15, 7638. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15217638

AMA Style

Gutema EM, Gopal M, Lemu HG. Minimization of Surface Roughness and Temperature during Turning of Aluminum 6061 Using Response Surface Methodology and Desirability Function Analysis. Materials. 2022; 15(21):7638. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15217638

Chicago/Turabian Style

Gutema, Endalkachew Mosisa, Mahesh Gopal, and Hirpa G. Lemu. 2022. "Minimization of Surface Roughness and Temperature during Turning of Aluminum 6061 Using Response Surface Methodology and Desirability Function Analysis" Materials 15, no. 21: 7638. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15217638

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