Next Article in Journal
Influence of Hole Localization on Local and Global Dynamic Response of Thin-Walled Laminated Cantilever Beam
Next Article in Special Issue
Analysis and Optimization of Dimensional Accuracy and Porosity of High Impact Polystyrene Material Printed by FDM Process: PSO, JAYA, Rao, and Bald Eagle Search Algorithms
Previous Article in Journal
Evaluation and Multi-Objective Optimization of Lightweight Mortars Parameters at Elevated Temperature via Box–Behnken Optimization Approach
Previous Article in Special Issue
Study of a Multicriterion Decision-Making Approach to the MQL Turning of AISI 304 Steel Using Hybrid Nanocutting Fluid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integration of Fuzzy AHP and Fuzzy TOPSIS Methods for Wire Electric Discharge Machining of Titanium (Ti6Al4V) Alloy Using RSM

1
Department of Mechanical Engineering, School of Technology, Pandit Deendayal Energy University, Raysan, Gandhinagar 382007, India
2
Department of Automated Mechanical Engineering, South Ural State University, 454080 Chelyabinsk, Russia
3
School of Mechanical and Design Engineering, University of Portsmouth, Portsmouth PO1 3DJ, UK
4
Faculty of Mechanical Engineering and Management, Poznan University of Technology, 60-965 Poznan, Poland
*
Authors to whom correspondence should be addressed.
Submission received: 3 November 2021 / Revised: 16 November 2021 / Accepted: 16 November 2021 / Published: 3 December 2021

Abstract

:
Titanium and its alloys exhibit numerous uses in aerospace, automobile, biomedical and marine industries because of their enhanced mechanical properties. However, the machinability of titanium alloys can be cumbersome due to their lower density, high hardness, low thermal conductivity, and low elastic modulus. The wire electrical discharge machining (WEDM) process is an effective choice for machining titanium and its alloys due to its unique machining characteristics. The present work proposes multi-objective optimization of WEDM on Ti6Al4V alloy using a fuzzy integrated multi-criteria decision-making (MCDM) approach. The use of MCDM has become an active area of research due to its proven ability to solve complex problems. The novelty of the present work is to use integrated fuzzy analytic hierarchy process (AHP) and fuzzy technique for order preference by similarity to ideal situation (TOPSIS) to optimize the WEDM process. The experiments were systematically conducted adapting the face-centered central composite design approach of response surface methodology. Three independent factors—pulse-on time (Ton), pulse-off time (Toff), and current—were chosen, each having three levels to monitor the process response in terms of cutting speed (VC), material removal rate (MRR), and surface roughness (SR). To assess the relevance and significance of the models, an analysis of variance was carried out. The optimal process parameters after integrating fuzzy AHP coupled with fuzzy TOPSIS approach found were Ton = 40 µs, Toff = 15 µs, and current = 2A.

1. Introduction

Titanium and its alloys possess a high strength-to-weight ratio which can be retained at high temperatures [1]. Owing to very high corrosion and erosion resistance, these metals are versatile in nature and they find applications in the pharmaceutical, aerospace, marine, chemical engineering, and food industries [2,3]. They have excellent bio-compatibility and as a consequence, they have been broadly utilized in biomedical applications and surgical implants [4]. Ti6Al4V alloy is used as an armor material for military vehicles, which require excellent ballistic performance. The titanium implants in a patient’s body permit examination with MRIs and NMRIs due to their non-ferromagnetic property [5,6]. However, titanium has very poor thermal conductivity, which leads to the localization of heat at the point of contact of the tool with the chip, resulting in high thermal gradients within the machining zone [7,8]. This in turn leads to increased tool wear rate and eventually tool failure. Moreover, titanium is chemically reactive at elevated temperatures, which may cause the tool to weld with the metal, leading to its premature failure [9]. Titanium has a low elastic modulus responsible for the deflection of job, chatter, and vibrations while machining. Given these limitations of conventional machining of titanium, there is a crucial need to explore the machinability of titanium with non-traditional processes, and wire electrical discharge machining (WEDM) is one of the process to investigate for effective machining.
The WEDM process is used to machine materials that are electrically conductive despite their mechanical properties [10,11,12]. The material removal in WEDM occurs with the input of the pulsating DC power supply, which discharges the electric current along the narrow inter-electrode gap between the workpiece and the tool [13,14,15]. The fine, clean, and precise cuts on the workpiece are possible because the wire, which acts as a cutting tool, possesses a small diameter and significant mechanical properties [16,17].
For sustainable manufacturing, the selection of optimum process parameters is indeed important. Several researchers worked on the optimization of the WEDM parameters for the machining of the Ti6Al4V alloy. The systematic design of experiments is essential to extract information from the experiments with minimal loss of resources in terms of time, money, raw materials, etc. Response surface methodology (RSM), Taguchi techniques using an orthogonal array, and fractional factorial are some of the methods for experimental design. The recorded output from the various experiments using one of the above methods is further used to optimize input process parameters. Multi-response optimization is needed in contemporary manufacturing. Various methods such as gray relational analysis (GRA) [11,15,18], heat transfer search (HTS) algorithm [2,19], teacher learning-based algorithm [20], particle swarm algorithm [21], genetic algorithm [22], artificial neural networks [23], etc., have been attempted to check their feasibility in finding the tread-off solution in terms of optimized process parameters.
Developing mathematical models of the WEDM process using various techniques is discussed in several studies. Mathematical models are helpful in predicting responses without performing extensive experiments. Sharma et al. [24] investigated WEDM parameters’ effect on dimensional deviation and cutting speed using a high-strength, low-alloy steel workpiece. They developed a mathematical relation of response surface using second-order CCD of RSM. Based on the mathematical model of cutting speed, they revealed that Ton, Toff, spark gap voltage, peak current, Ton × Toff, Ton × IP, Ton × wire tension, and spark gap voltage × wire tension have significant effects on cutting speed. Furthermore, a mathematical model of dimensional deviation suggests that the main effects of Ton, Toff, SV, IP, and WT, and interaction effects between Ton and SV, Ton and IP are statistically significant for the analysis. Kavimani et al. [25] proposed a mathematical model for predicting the response data of MRR and surface roughness by using regression analysis with orthogonal array. The mathematical model proposed by Bose and Nandi [26] is based on two factors and four levels design of experiments on output response, such as surface roughness. They used desirable gray relational analysis algorithm for optimization of the WEDM process and model building. Daniel et al. [27] analyzed the effects of several WEDM process parameters on MRR and developed a mathematical model using regression analysis for predicting the response.
Arikatla et al. [28] optimized the parameters during the WEDM of titanium (Ti6Al4V) alloy using “RSM” method. The effects of pulse-on time (Ton), wire tension, servo voltage, and pulse-off time (Toff) were studied on the kerf width, surface roughness (SR), and MRR of the material. An analysis of variance (ANOVA) was conducted to check the impact of the process variables on the desired responses. Chaudhari et al. [29] attempted machining of Ni55.8Ti super-elastic shape memory alloy with the WEDM process. To exhibit the viability for industrial applications, a scientific approach consisting of RSM and HTS algorithm was planned and prepared for optimization. Saedon et al. [30] studied optimization of kerf and MRR during cutting of titanium (Ti6Al4V) alloy using WEDM. The experiments were performed with varying peak current, Toff, wire tension, and wire feed. The design of the process parameters for the experiment was determined using Taguchi’s L9 design. The GRA approach was utilized for emphasizing the study for multiple performance characteristics. Payal et al. [31] employed the Taguchi fuzzy integration for parametric optimization of EDM process with multiple response measures. A fuzzy model was formed through which the optimal blend of parameters was obtained on the basis of multi-performance fuzzy index values.
Nowadays, multi-criteria decision making (MCDM) mechanisms are emerging in the manufacturing field. These are powerful tools for solving problems with multiple criteria. Some MCDM techniques help in assigning systematic weights to the criteria and some assign hierarchy of the alternatives [32]. Ananthakumar et al. [33] used the MCDM method, namely Technique for Order Preference by Similarity to Ideal Solution (TOPSIS), to determine the desired cutting conditions of plasma arc cutting. They defined the arc current of 45A, standoff distance of 2 mm, gas pressure of 3 bar, and cutting speed of 2400 mm/min as the best conditions for superior quality. Prabhu et al. [34] revealed that multiple quality characteristics of friction stir welded aluminum matrix composite can be improved using the TOPSIS approach. Tamjidy et al. [35] selected the best optimal solution in friction stir welding (FSW) of AA7075-T6 and AA6061-T6 using TOPSIS and Shannon entropy method. Sudhagar et al. [36] compared the ranking performance of TOPSIS and GRA for FSW of Al2024 aluminum alloy. The analytical hierarchy process (AHP) is an MCDM technique used to assign systematic weights to the criterion. Gaidhani and Kalamani [37] selected the most influential processing parameters in abrasive water jet machining based on weights obtained with the help of the AHP technique of MCDM. Babu and Venkataramaiah [38] optimized input parameters in WEDM of Al6061/SiCP composite using the AHP-TOPSIS approach. They performed experiments using Taguchi’s L18 design considering sensitivity, Ton, wire feed rate (WF), Toff, and wire type as input parameters. A similar study was performed using the AHP–TOPSIS hybrid approach during WEDM by Nayak and Mahapatra [39]. The current research on optimization is focused on developing a hybrid approach of optimization by integrating different MCDM techniques with a fuzzy approach. The fuzzy approach is a method of solving problems that are associated with uncertainty and vagueness. Different types of uncertainty can be observed in different optimization and decision-making problems. Chou et al. [40] used the fuzzy AHP and fuzzy TOPSIS to assess the performance of human resources in science and technology (HRST) in Southeast Asian countries. A fuzzy AHP approach was implemented to decide the preference weights for various performance measures while fuzzy TOPSIS was used to identify the best tread-off alternatives to accomplish the ideal HRST levels. They concluded that countries such as Singapore, South Korea, and Taiwan have better HRST performances compared to other Southeast Asian countries. Sirisawat and Kiatcharoenpol [41] employed the fuzzy AHP and fuzzy TOPSIS techniques to classify the reverse logistics barriers and also to prioritize and rank the solutions to implementing reverse logistics in the electronics industry. A hybrid decision-making method with a fuzzy approach was also found effective in deciding the best alternative in the manufacturing field. Roy and Dutta [42] studied the working of the EDM process considering duty cycle, current, Ton, and gap voltage as controllable parameters. They used Taguchi’s L9 technique for experimental design and used integrated fuzzy AHP and fuzzy TOPSIS methods to determine the optimal set of controllable parameters. Furthermore, ANOVA analysis of closeness coefficient index (CCi) revealed that the current was the highest contributing factor among the selected parameters.
Previous studies suggest that many researchers attempted different MCDM techniques such as TOPSIS, AHP, fuzzy TOPSIS, and fuzzy AHP separately or in combination in different fields, including logistics, electronic industries, additive manufacturing, abrasive water jet machining, plasma arc cutting, etc. However, a study about investigating the use of response surface methodology using fuzzy AHP and fuzzy TOPSIS hybrid approach for WEDM machining of Ti6A4V has not yet been explored.
The present work intended to optimize the process parameters of Ton, Toff, and current during the WEDM of Ti6Al4V alloy. Cutting speed (VC), MRR, and SR were considered for the analysis as output response variables. The experiments were systematically designed using central composite design (CCD) of RSM and mathematical models were developed between responses and input parameters. The appropriateness of developed models was checked using ANOVA by analyzing the R-Seq values, lack-of-fit, and residual plots. Furthermore, the weights were assigned to considered criterion using fuzzy AHP, and optimal process parameter conditions were predicted using fuzzy TOPSIS. The confirmatory study was performed to verify the results of optimization.

2. Materials and Methods

2.1. Experimental Methodology

In the present study, a WEDM machine (Concord WEDM, DK7732, Concord Limited, Bangalore, India) was used to perform the experiments using Ti6Al4V as the work material (with a dimension of 90 × 50 × 6 mm3) and molybdenum wire (180 µm diameter) as the tool electrode. The dielectric fluid utilized was deionized water. Table 1 shows the chemical composition of the selected work material of Ti6Al4V. Ton, Toff, and current were selected as input process parameters based on recent literature to investigate their effects on WEDM machining in terms of cutting speed, MRR, and SR. Full factorial CCD of RSM was selected to design the experimental plan. CCD is very popular amongst multiple variants of RSM approaches because it offers great flexibility and allows sequential operations and effectiveness by providing the optimum solution in a minimum number of iterations. CCD design of RSM for three factors at three levels was implemented to prepare the experimental matrix. The levels of factors were selected based on preliminary trials and literature study. Table 2 shows the input parameters of the WEDM process at the three levels. Table 3 shows the full factorial CCD composed of the six axial and central points, and eight factorial points. During each experimental run, the wire was fed along a width direction of the work material. Each experiment was repeated three times; the average values of the three trials were reported and considered for the analysis. Minitab 17 software was utilized for the RSM design and analyzing the experimental data.
In total, 20 experiments were run randomly to avoid an experimental error. The output characteristic cutting speed refers to the speed in mm/min at which the wire cuts the workpiece. It is calculated using Equation (1):
V c = L t
where L is the distance of 50 mm width of the workpiece which was cut in each pass and t is the machining time in minutes.
The material removal rate is considered as the volumetric material removed per unit time, which is calculated using Equation (2):
MRR = (wt.initial − wt.final)/(ρ ∗ t)
where wt.initial is the weight of the workpiece measured before the cut and wt.final is the weight of workpiece measured after the cut. ρ is the density of the titanium alloy, 4.42 gm/cm3. t is the cutting time measured for a particular cut in one minute.
Surface roughness is an important parameter that defines the quality of the machined surfaces. The SR of the WEDM part was determined with the use of the Surftest SJ-410 model. The arithmetic average roughness (Ra) value was determined in µm from cut specimens in the current study. The cutoff length of 0.8 mm and the evaluation length of 20 mm were used for the measurement of SR.

2.2. Fuzzy Analytical Hierarchy Process

To use fuzzy TOPSIS for ranking of alternatives, each response needs to provide priority values or weights. These weights may vary from person to person. To overcome this limitation, T. L. Saaty [43] developed a technique called the analytic hierarchy process (AHP). This technique decomposes the decision-making situation into a systematic hierarchy of objectives, attributes, and alternatives. AHP is a compelling tool for complicated decision-making situations and helps decision makers define objectives and reach the best possible choice. By simplifying the complex choices to a progression of pairwise observations, and then integrating the outcomes, the AHP catches both abstract and target parts of a decision [44]. The limitation of conventional AHP is that it works with crisp information derived from linguistic responses. The scale of converting linguistic responses into crisp data in AHP is very unbalanced. The judgment of decision makers largely affects the result of conventional AHP due to aleatory uncertainty present in the natural language. This attracted researchers to integrate fuzzy theory with AHP and hence, fuzzy AHP was developed. The steps of fuzzy AHP in this research are as follows [42]:
Step 1:
Construct the various leveled structure of objective, criterion, and alternatives of the problem.
Step 2:
Construct a pairwise comparison matrix from the criteria/options available. Furthermore, assign linguistic terms using Figure 1 to the pairwise comparisons collected from decision makers. Convert linguistic terms into fuzzy numbers using Table 4. The generalized pairwise comparison matrix will be of the form shown in Equation (3).
A ˜ = ( 1 y ˜ 1 n 1 y ˜ 1 n 1 )
where y ˜ i j = 1; i = j.
Step 3:
Fuzzification is used to convert the linguistic term into a membership term. The fuzzification of the linguistic term can be possible using various functions such as triangular, bell-shaped, and trapezoidal functions. For this study, we used the triangular membership function, as shown in Figure 2. The assumed fuzzy numbers are shown in Equations (4) and (5).
X = ( r , s , t )
Y = ( l , m , n )
where r, s, and t denote the lower, middle, and upper bounds of fuzzy number X, respectively, and l, m, and n denote the lower, middle, and the upper bounds of fuzzy number Y, respectively.
Fuzzy weights can be found using fuzzy addition and multiplication [45]. The generalized fuzzy addition and fuzzy multiplication formulas are expressed by Equations (6) and (7).
Fuzzy addition:
X ˜ Y ˜ = ( r , s , t ) ( l , m , n ) = ( r + l , s + m , t + n )
Fuzzy multiplication:
X ˜ Y ˜ = ( r , s , t ) ( l , m , n ) = ( r × l , s × m , t × n )
Step 4:
Determine the fuzzy mean geometric value (FMGV) of each criteria using the geometric mean method. Equation (8) can be used for calculating FMGV. The fuzzy weights can be determined by using Equation (9).
r ˜ j = ( y ˜ j 1 y ˜ j 2 y ˜ j n ) 1 n
w ˜ j = r ˜ j ( r ˜ 1 r ˜ 2 r ˜ n ) 1
where y ˜ i j is the comparison of fuzzy value from criterion i to j; r ˜ j is the geometric mean value for comparison of the fuzzy value of criterion j to every other criterion; w ˜ j is the fuzzy weight of each criterion.

2.3. Fuzzy TOPSIS

Hwang and Yoon developed a multi-criteria decision analysis technique in 1981 called the Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS). It was further modified by Yoon in 1987, and in 1993, it was further improved by Hwang Lai and Liu [46]. TOPSIS helps in choosing the best alternative with the least proximity from the positive ideal solution and the furthest from the negative ideal solution [47]. It is the strategy that compares the set of responses by identifying the contribution (weights) of each criterion. As the criteria are generally of random measurements, this might create problems in assessments. To avoid this situation, the need for fuzzy numbers is fundamental. Utilizing the fuzzy numbers in TOPSIS for multi-criteria decision making makes it easy for assessment [48]. The proposed steps for fuzzy TOPSIS implementation in this research are given below:
Step 1:
Normalization of response: The normalization is important for converting measured outputs into the fuzzy number. The process of normalization was carried out considering the output based on the benefit criteria or the cost criteria. The VC and MRR were normalized using the benefit criteria using Equation (10), whereas SR was normalized using the cost criteria using Equation (11).
For benefit criteria:
r i j ( x ) = x i j x m i n x m a x x m i n
For cost criteria:
r i j ( x ) = x m i n x i j x m i n x m a x
where r i j ( x ) is the normalized value of output, x m a x is the maximum of x i j   and x m i n is the minimum of x i j .
After normalization of the ratings, the responses of experiments were converted to fuzzy linguistic variables using Table 5. The five-level fuzzy linguistic variables are represented using a triangular fuzzy number [49].
Step 2:
Fuzzification of normalized decision matrix: The decision matrix normalized in Step 1 can be converted to a fuzzified normalized decision matrix by assigning a sub-criteria grade to each alternative using Table 5 of the K membership function scale. Additionally, assign the weights to each sub-criteria grade.
The weight of criteria:
w ˜ j k = ( w ˜ j 1 k , w ˜ j 2 k , w ˜ j 3 k )
Step 3:
Calculate the weighted normalized fuzzy decision matrix: The weights obtained from fuzzy AHP are required to construct this matrix. The weighted normalized values can be calculated as:
V ˜ = ( v ˜ i j )   Where   v ˜ i j = r ˜ i j × w j
Step 4:
Identify the positive ideal (V+) and negative ideal (V) solutions: The fuzzy positive ideal solutions (FPIS, V+) and the fuzzy negative ideal solutions (FNIS, V) must be calculated using Equations (14) and (15).
V + = ( v ˜ 1 + , v ˜ 2 + , v ˜ 3 + ) , where:
v ˜ j + = max i { v i j 3 }
V = ( v ˜ 1 , v ˜ 2 , v ˜ 3 ) , where:
v ˜ j = min i { v i j 1 }
Consideration of the maximum and minimum of Vij does not necessarily result in triangular fuzzy numbers, but we can obtain the ideal solutions as the fuzzy numbers using Equation (16).
d ( p ˜ ,   q ˜ ) = 1 3 { ( p 1 q 1 ) 2 + ( p 2 q 2 ) 2 + ( p 3 q 3 ) 2 }
where p ˜ = ( p 1 , p 2 ,   p 3 ) and q ˜ = ( q 1 ,   q 2 ,   q 3 ).
Step 5:
Calculate separation measures: The separation measure d i + is the summation of the distance of each response to the FPIS and d i is the summation of the distance of each response to the FNIS. The distance can be calculated by using the following equations.
d i + = j = 1 n d ( v ˜ i j , v ˜ j + ) ;   i = 1 ,   2 ,   ,   m
d i = j = 1 n d ( v ˜ i j , v ˜ j ) ;   i = 1 ,   2 ,   ,   m
Step 6:
Calculate the similarities to the ideal solution: To solve the similarities, compute the closeness coefficient CCi for each alternative [48].
C C i = d i d i + + d i ;   C C i [ 0 , 1 ]   i = 1 ,   2 ,   ,   n

3. Results and Discussions

3.1. Regression Equations

Evaluated values of the selected output response variables such as cutting speed, MRR, and SR are mentioned in Table 6. These measured response variables are then normalized using Equations (10) and (11) for further analysis. Table 6 shows the uncoded actual values of input process parameters as per “RSM” design. Minitab 17 software was used to find regression coefficients, and subsequently, ANOVA was assessed. The significance of the coefficients was tested at a 95% confidence level which is essential to recognize the most influencing model terms [50]. With the help of regression analysis, a mathematical relationship was obtained for VC, MRR, and SR in terms of input process variables. Equations (20)–(22) show the regression equations for VC, MRR, and SR, respectively.
The developed regression models consist of linear, quadratic, and interaction terms. It becomes is essential to study the significance of these terms on the output parameters. This can be executed by performing ANOVA. Thus, the ANOVA investigation was further studied to predict the significant terms of the proposed study. ANOVA tests the hypothesis based on the equality of means when several factors are considered. It is a statistical inference technique used to determine the influence of the hypothesis made for the model [51]. The study of the influence of input parameters at various levels is important for single-objective optimization. The main effect plots of the response highlight the optimum factor-level combinations for a given response. A detailed discussion is provided on mathematical model development, model adequacy checking, and the main effect plot in the following section for each response.
VC = 1.429 + 0.00567(Ton) − 0.166(Toff) + 1.613(Current) − 0.00005(Ton × Ton) + 0.00233(Toff × Toff) − 0.1699(Current × Current) + 0.000265(Ton × Toff) + 0.001450(Ton × Current) − 0.01280(Toff × Current)                                           
MRR = 1.37 + 0.0153(Ton) − 0.118(Toff) + 1.552(Current) − 0.000111(Ton × Ton) + 0.00062(Toff × Toff) − 0.1345(Current × Current) + 0.000267(Ton × Toff) + 0.00183(Ton × Current) − 0.01450(Toff × Current)                                            
SR = −1.48 + 0.1273(Ton) − 0.341(Toff) + 2.66(Current) − 0.000664(Ton × Ton) + 0.0133(Toff × Toff) − 0.318(Current × Current) − 0.00167(Ton × Toff) + 0.01204(Ton × Current) − 0.0443(Toff × Current)                                          

3.2. Analysis of Cutting Speed

ANOVA was carried out for cutting speed assuming a 95% confidence level. The ANOVA results are summarized in Table 7. F-value of the model (112.30) is much larger than the F-stat value (3.02). This suggests that the model is largely significant. There is just a 0.01% probability that a model F-value this large could be due to noise. The lack of fit F-value of 3.38 implies there is a 10.4% chance that a lack of fit F-value this large could be due to noise. For a confidence level of 95%, the p-value for any input parameter should be less than 0.05 to consider that parameter as significant [52,53]. It can be observed from Table 7 that all the input process variables have a p-value less than 0.05, which shows that all input process parameters are significant for VC. Toff shows the highest significance for obtaining higher value of VC followed by current and Ton. The p-value of the model is also less than 0.05, highlighting that the model is significant and best fitted for the selected range of process parameters. Insignificant lack-of-fit reveals the adequacy and fitness of the model [10]. The value of R2 indicates that 99.02% of the variation of cutting speed is contributed by the control factors and only 0.98% of total variation cannot be described by the quadratic model. The adjacent R-squared’ and predicted R-squared values of 98.14% and 94.04%, respectively, are in reasonable agreement. Figure 3 shows the normal probability plot of residuals for VC. ANOVA results are considered to be valid depending on the analysis of these plots. We observed that the developed regression model fit well with the observed values.
Figure 4 shows the influence of the three parameters at various levels on cutting speed. The increase in current elevates the discharge energy, which causes the rise in VC. The increase in the value of Ton signifies a rise in the duration of a spark, which causes the discharge energy to increase [2]. An increase in Ton and current significantly increases the spark intensity, which in turn escalates the melting and vaporization of the material from workpiece [7]. Hence, an increase in both Ton and current increases VC. However, a continuous decrease in the value of VC has been observed with an increase in the value of Toff due to the absence of the spark during the machining [10].

3.3. Analysis of MRR

The ANOVA results for the quadratic model developed between considered inputs and MRR at a 95% confidence level are shown in Table 8. The results indicate that the F-value of the model is 63.96 with a corresponding p-value of 0.000, which is less than 0.05, implying that the quadratic model is significant at the considered confidence level. A lack of fit of 4.46 with a corresponding p-value of 0.063 implies that it is not significant and suggests the fitness of the model. All the input process parameters (Ton, Toff and current) show p-values of less than 0.05, suggesting that all are having a significant effect on MRR. The regression equation and ANOVA results of MRR revealed that MRR is significantly dependent on Toff. Toff is most contributing factor (41.73%) on MRR followed by Ton and current. Saedon et al. [30] also reported Toff as the most significant factor (58%) while investigating WEDM machining on Ti-6Al-4V alloy. The R2 value of 98.29% indicates that 98.29% of the variation of MRR can be explained by the empirical model and only 1.71% of total variation cannot be described by the developed model. When predicted R2 and adjusted R2 values are in reasonable agreement, this confirms a strong correlation between observed and predicted values. Here, adjusted R2 is 96.76%, and the predicted value of R2 is 88.69%. The closeness of the values depicts a strong correlation between them. Figure 5 shows the residual plots for MRR in terms of normal probability plot versus residual, residual versus fitted values, histogram, and residual versus observation order. All the residual plots indicate that the regression model fits well with the observed values.
From Figure 6, it is inferred that with the rise in Ton and current, the influence of clearance form on MRR exhibits an increasing tendency. An increase in Ton and current significantly increases the spark intensity which in turn escalates the melting and vaporization of the material from the workpiece [10]. This further increases the MRR by a large amount. Saedon et al. [30] attributed increasing MRR with increasing Ton and current to the reduced dynamic shear strength of Ti alloy due to higher thermal influence in the machining region. The outcome of Toff on MRR shows a decreasing trend with a rise in Toff because of reduced spark ejection time and less MRR [9,54]. Thus, with an increase in Toff, MRR is decreasing. Moreover, the slope indicates that it has a great effect on MRR; a slight increase in Toff leads to a decrease in MRR [9]. The positive dependency/correlation/relationship of MRR with Ton and negative dependency with Toff is also reported by Arikatla et al. [28].

3.4. Analysis of SR

The ANOVA results for the quadratic model developed between considered inputs and machining time at a 95% confidence level are shown in Table 9. The results indicate that the F-value of the model is 11.67 with a corresponding p-value of 0.000, which is less than 0.05, implying that the quadratic model is significant at the considered confidence level. Ton and current were found to be the significant input process parameters for SR, with a higher contribution of Ton followed by current. The ANOVA results of SR revealed Ton as the main contributing factor (50.09%). Similarly, the largest contributing effect of Ton on SR was also reported by Arikatla et al. [28] A lack of fit of 3.73 with a corresponding p-value of 0.087 implies that it is not significant. The value of R2 of 0.9131 indicates that 91.31% of the variation of surface roughness can be explained by the empirical model and only 8.69% of total variation cannot be described by the developed model. When predicted R2 and adjusted R2 values are in reasonable agreement, it indicates strong relationship in observed and predicted values. Here, adjusted R2 is 0.9842 and the predicted value of R2 is 0.9473. Table 10 shows the model summary for all the responses. The closeness of the values depicts a strong correlation between them. Figure 7 shows the normal probability plot of residuals for SR. The plots highlight that that developed regression model fits well with the observed values.
Figure 8 shows the impact of machining parameters on SR. The effect of Ton is directly proportional to SR and has a constantly increasing trend because the higher the Ton, the higher the energy of discharge and spark intensity, resulting in poor surface texture and vice versa [55]. This indicates that with the increase in Ton, the SR value also increases. The effect of Toff on SR shows a reducing trend with growth in Toff (inversely proportional) [50,55]. Thus, with the increase in Toff, SR is decreasing as discharge energy falls with the rise in Toff, leading to the plunging of the crater dimensions and hence reducing SR. The effect of current on SR shows a growing tendency with a rise in current value and vice versa. Higher SR with the high current can be explained by the fact that crater sizes are also dependent on the ionization of the dielectric fluid and ionization takes place at faster rate with a higher current [50]. Thus, crater size increases with an increase in current value, resulting in higher SR. The positive dependency of SR on Ton and negative dependency on Toff was also reported by Arikatla et al. [28].
The ANOVA analysis of all three responses showed that developed mathematical models are significant for predicting the responses. It indicates that the experimental error is very minimal and collected output data can be used for multi-objective optimization.
The main effect plots of the considered responses (Figure 4, Figure 6 and Figure 8) can be used to identify optimum factor-level combinations from a single-objective optimization perspective. In the present study, VC and MRR are of the “higher the better” category and SR falls under the “lower the better category”. Thus, the optimum parameter settings obtained considering a single objective are shown in Table 11. It can be observed that the optimum settings are conflicting in nature when all three responses are considered together. The industry demands process parameter setting, which can result in higher productivity (higher VC and MRR) with good quality (low SR). This requires an improved means of optimization which can take care of such conflicting situations.

3.5. Optimization Using Integrated Fuzzy AHP and Fuzzy TOPSIS

3.5.1. Fuzzy AHP

The hierarchical structure of the objective, criterion, and alternatives of the problem is constructed as per Step 1 of fuzzy AHP. The linguistic data of the pairwise comparison matrix were collected from domain experts and converted to crisp values using Table 4. The pairwise comparison matrix with crisp values was prepared using Equation (3) and is shown in Table 12.
The consistency ratio (CR) was calculated as 0.028. A value of CR less than 0.1 confirms good consistency in the judgments made by domain experts while assigning values in a pairwise comparison matrix [56]. Triangular fuzzy numbers from Table 4 were used for converting crisp values in the pairwise comparison matrix to fuzzy numbers. The obtained fuzzified comparison matrix is shown in Table 13.
The fuzzy mean geometric values were calculated using Equation (8) and were used to find fuzzy weights. The fuzzy weights of the criterion were calculated using Equation (9) and are shown in Table 14. These weights represent the lower, modal, and upper values of the fuzzy numbers, respectively.

3.5.2. Fuzzy TOPSIS

In fuzzy TOPSIS, the alternative selection criteria were decided from the CCD design of the response surface methodology. Initially, the response data obtained were normalized using Equations (10) and (11) and are shown in Table 6. Then, a fuzzified normalized decision matrix was obtained by assigning weights to each sub-criteria grade using Equation (12). The fuzzified normalized decision matrix obtained is shown in Table 15.
The fuzzy weights obtained from the fuzzy AHP technique (Table 13) were then multiplied to each performance rating using Equation (13), resulting in a weighted normalized fuzzy decision matrix as shown in Table 16.
FPIS and FNIS were calculated using Equations (14) and (15). In the present study, VC and MRR are benefit criteria and SR is cost criteria. Thus, we assumed FPIS V j + as (1, 1, 1) for VC, MRR and (0, 0, 0) for SR. Additionally, we assigned FNIS V j as (0, 0, 0) for VC, MRR and (1, 1, 1) for SR.
Equations (17) and (18) were used to find FPIS and FNIS separation measures for each alternative, and the obtained results of separation measures are shown in Table 17. Finally, the closeness coefficient (CCi) was found for each alternative using Equation (19). The value of CCi indicates whether the alternative is nearest to theoretical FPIS and furthest from the theoretical FNIS or vice versa [57]. The highest rank is given to the alternative with the highest value of the closeness coefficient, as shown in Table 17.
After calculating the closeness coefficient using fuzzy TOPSIS, the optimized process parameters were determined using a single to noise (S/N) ratio of CCi values. High CCi values are always preferred. Hence, the (S/N) ratio was calculated using the “higher the better” strategy in Minitab. The obtained main effect plot of the (S/N) ratio is shown in Figure 9. The graph indicates the effect of factors’ concerning level. The optimum levels of process parameters were picked based on a higher value of the S/N ratio for the levels of factors. The optimized process parameters based on the highest values S/N ratio were a Ton of 40 µs, Toff of 15 µs, and current of 2A. The determined optimized process parameters correspond to alternative 12.
Confirmation experiments were performed to validate the optimal settings. Table 18 shows that the experimental results correlate well with the predicted results.

4. Conclusions

In this study, multi-criteria decision-making techniques, fuzzy AHP and fuzzy TOPSIS, were integrated for parameter optimization problems in Ti6Al4V alloy machining. The conclusions of the present work are summarized as follows:
  • Response surface methodology is effective for systematically designing the experiments. The mathematical relations developed between dependent and independent parameters are significant for predicting the responses at a 95% confidence interval.
  • ANOVA analysis confirmed that the input parameters Ton, Toff, and current significantly affect cutting speed, material removal rate, and surface roughness.
  • Fuzzy AHP can be incorporated to prioritize the responses using data collected from experts. The use of the fuzzy approach eliminates the aleatory uncertainty present in the natural language. The weights calculated using fuzzy AHP can be incorporated in fuzzy TOPSIS without bias.
  • For the considered range of process parameters, the optimal process parameters for WEDM are Ton = 40 µs, Toff = 15 µs, and current = 2A.
  • The confirmatory experiments proved that fuzzy logic is an effective and efficient solution for the optimization of WEDM process parameters. The proposed integrated approach of RSM, fuzzy AHP, and fuzzy TOPSIS can be further extended for different machining processes.

Author Contributions

Conceptualization, K.F.; methodology, D.M.; software, A.D. and D.M.; validation, J.V., R.C. and P.P.; formal analysis, K.F., D.Y.P. and S.W.; investigation, K.F., R.C., J.V. and S.W.; resources, D.Y.P.; data curation, K.F., R.C., K.G. and D.Y.P.; writing—original draft preparation, K.F., J.V., R.C., D.Y.P., P.P. and K.G.; writing—review and editing, J.V., R.C. and S.W.; visualization, K.F.; supervision, R.C. and S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data presented in this study are available in this article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AHPAnalytical hierarchy process
ANOVAAnalysis of variance
CCiCloseness coefficient index
CCDCentral composite design
CRConsistency ratio
FMGVFuzzy mean geometric value
FPISFuzzy positive ideal solutions
FNISFuzzy negative ideal solutions
GRAGray relational analysis
HTSHeat transfer search
MCDMMulti-criteria decision making
MRRMaterial removal rate
RSMResponse surface methodology
S/NSingle to noise
SRSurface roughness
TOPSISTechnique for Order Preference by Similarity to Ideal Solution
TonPulse-on time
ToffPulse-off time
VCCutting speed
WFWire feed rate
WEDMWire electrical discharge machining process

References

  1. Ezugwu, E.; Wang, Z. Titanium alloys and their machinability—A review. J. Mater. Process. Technol. 1997, 68, 262–274. [Google Scholar] [CrossRef]
  2. Chaudhari, R.; Vora, J.J.; Prabu, S.M.; Palani, I.; Patel, V.K.; Parikh, D. Pareto optimization of WEDM process parameters for machining a NiTi shape memory alloy using a combined approach of RSM and heat transfer search algorithm. Adv. Manuf. 2019, 9, 64–80. [Google Scholar] [CrossRef]
  3. Baltatu, M.S.; Vizureanu, P.; Sandu, A.V.; Florido-Suarez, N.; Saceleanu, M.V.; Mirza-Rosca, J.C.J.M. New Titanium Alloys. Promis. Mater. Med Devices 2021, 14, 5934. [Google Scholar]
  4. Nicholson, J.J.P.W. Titanium alloys for dental implants: A review. Prosthesis 2020, 2, 100–116. [Google Scholar] [CrossRef]
  5. Saravanan, R.; Rani, M.P. Metal and Alloy Bonding—An Experimental Analysis: Charge Density in Metals and Alloys; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  6. Vora, J.; Chaudhari, R.; Patel, C.; Pimenov, D.Y.; Patel, V.K.; Giasin, K.; Sharma, S. Experimental Investigations and Pareto Optimization of Fiber Laser Cutting Process of Ti6Al4V. Metals 2021, 11, 1461. [Google Scholar] [CrossRef]
  7. Hourmand, M.; Sarhan, A.A.; Sayuti, M.; Hamdi, M.J.A. A Comprehensive Review on Machining of Titanium Alloys. Engineering 2021, 1–37. [Google Scholar] [CrossRef]
  8. Chaudhari, R.; Vora, J.; Parikh, D.; Wankhede, V.; Khanna, S. Multi-response Optimization of WEDM Parameters Using an Integrated Approach of RSM–GRA Analysis for Pure Titanium. J. Inst. Eng. Ser. D 2020, 101, 117–126. [Google Scholar] [CrossRef]
  9. Gupta, V.; Singh, B.; Mishra, R.K. Machining of titanium and titanium alloys by electric discharge machining process: A review. Int. J. Mach. Mach. Mater. 2020, 22, 99–121. [Google Scholar]
  10. Chaudhari, R.; Vora, J.J.; Pramanik, A.; Parikh, D. Optimization of Parameters of Spark Erosion Based Processes. Spark Erosion Machining; CRC Press: Boca Raton, FL, USA, 2020; pp. 190–216. [Google Scholar]
  11. Sheth, M.; Gajjar, K.; Jain, A.; Shah, V.; Patel, H.; Chaudhari, R.; Vora, J. Multi-objective optimization of inconel 718 using Combined approach of taguchi—Grey relational analysis. In Advances in Mechanical Engineering; Springer: Singapore, 2021; pp. 229–235. [Google Scholar]
  12. Joshi, R.; Zinzala, G.; Nirmal, N.; Fuse, K. Multi-Response Optimization of EDM for Ti-6Al-4V Using Taguchi-Grey Relational Analysis. Solid State Phenomena; Trans Tech Publications Ltd.: Bäch, Switzerland, 2017. [Google Scholar]
  13. Aggarwal, V.; Pruncu, C.I.; Singh, J.; Sharma, S.; Pimenov, D.Y.J.M. Empirical investigations during WEDM of Ni-27Cu-3.15 Al-2Fe-1.5 Mn based superalloy for high temperature corrosion resistance applications. Materials 2020, 13, 3470. [Google Scholar] [CrossRef]
  14. Sen, B.; Hussain, S.A.I.; Gupta, A.D.; Gupta, M.K.; Pimenov, D.Y.; Mikołajczyk, T.J.M. Application of Type-2 Fuzzy AHP-ARAS for Selecting Optimal WEDM Parameters. Metals 2021, 11, 42. [Google Scholar] [CrossRef]
  15. Rathi, P.; Ghiya, R.; Shah, H.; Srivastava, P.; Patel, S.; Chaudhari, R.; Vora, J. Multi-Response Optimization of Ni55. 8Ti Shape Memory Alloy Using Taguchi–Grey Relational Analysis Approach. Recent Advances in Mechanical Infrastructure; Springer: Singapore, 2020; pp. 13–23. [Google Scholar]
  16. Goyal, A. Investigation of material removal rate and surface roughness during wire electrical discharge machining (WEDM) of Inconel 625 super alloy by cryogenic treated tool electrode. J. King Saud Univ.-Sci. 2017, 29, 528–535. [Google Scholar] [CrossRef]
  17. Lenin, N.; Sivakumar, M.; Selvakumar, G.; Rajamani, D.; Sivalingam, V.; Gupta, M.K.; Mikolajczyk, T.; Pimenov, D.Y. Optimization of Process Control Parameters for WEDM of Al-LM25/Fly Ash/B4C Hybrid Composites Using Evolutionary Algorithms: A Comparative Study. Metals 2021, 11, 1105. [Google Scholar] [CrossRef]
  18. Patel, S.; Fuse, K.; Gangvekar, K.; Badheka, V. Multi-response optimization of dissimilar Al-Ti alloy FSW using Taguchi-Grey relational analysis. In Key Engineering Materials; Trans Tech Publications Ltd.: Bäch, Switzerland, 2020. [Google Scholar]
  19. Vora, J.; Patel, V.K.; Srinivasan, S.; Chaudhari, R.; Pimenov, D.Y.; Giasin, K.; Sharma, S. Optimization of Activated Tungsten Inert Gas Welding Process Parameters Using Heat Transfer Search Algorithm: With Experimental Validation Using Case Studies. Metals 2021, 11, 981. [Google Scholar] [CrossRef]
  20. Chaudhari, R.; Khanna, S.; Vora, J.; Patel, V.K.; Paneliya, S.; Pimenov, D.Y.; Wojciechowski, S. Experimental investigations and optimization of MWCNTs-mixed WEDM process parameters of nitinol shape memory alloy. J. Mater. Res. Technol. 2021, 15, 2152–2169. [Google Scholar] [CrossRef]
  21. Nain, S.S.; Garg, D.; Kumar, S. Investigation for obtaining the optimal solution for improving the performance of WEDM of super alloy Udimet-L605 using particle swarm optimization. Eng. Sci. Technol. Int. J. 2018, 21, 261–273. [Google Scholar] [CrossRef]
  22. Sharma, N.; Khanna, R.; Gupta, R.D. WEDM process variables investigation for HSLA by response surface methodology and genetic algorithm. Eng. Sci. Technol. Int. J. 2015, 18, 171–177. [Google Scholar] [CrossRef] [Green Version]
  23. Phate, M.R.; Toney, S.B. Modeling and prediction of WEDM performance parameters for Al/SiCp MMC using dimensional analysis and artificial neural network. Eng. Sci. Technol. Int. J. 2019, 22, 468–476. [Google Scholar] [CrossRef]
  24. Sharma, N.; Khanna, R.; Gupta, R.D.; Sharma, R. Modeling and multiresponse optimization on WEDM for HSLA by RSM. Int. J. Adv. Manuf. Technol. 2013, 67, 2269–2281. [Google Scholar] [CrossRef]
  25. Kavimani, V.; Prakash, K.S.; Thankachan, T. Multi-objective optimization in WEDM process of graphene–SiC-magnesium composite through hybrid techniques. Measurement 2019, 145, 335–349. [Google Scholar] [CrossRef]
  26. Bose, S.; Nandi, T. A novel optimization algorithm on surface roughness of WEDM on titanium hybrid composite. Sādhanā 2020, 45, 236. [Google Scholar] [CrossRef]
  27. Daniel, G.; Kumar, P.M. Parametric study and optimization of regression model in WEDM using genetic algorithm. IOP Conf. Ser. Mater. Sci. Eng. 2020, 764, 012047. [Google Scholar] [CrossRef]
  28. Prasad Arikatla, S.; Tamil Mannan, K.; Krishnaiah, A. Parametric Optimization in Wire Electrical Discharge Machining of Titanium Alloy Using Response Surface Methodology. Mater. Today Proc. 2017, 4, 1434–1441. [Google Scholar] [CrossRef]
  29. Chaudhari, R.; Vora, J.J.; Mani Prabu, S.S.; Palani, I.A.; Patel, V.K.; Parikh, D.M.; de Lacalle, L.N.L. Multi-Response Optimization of WEDM Process Parameters for Machining of Superelastic Nitinol Shape-Memory Alloy Using a Heat-Transfer Search Algorithm. Materials 2019, 12, 1277. [Google Scholar] [CrossRef] [Green Version]
  30. Saedon, J.; Jaafar, N.; Yahaya, M.A.; Nor, N.M.; Husain, H. A study on kerf and material removal rate in wire electricaldischarge machining of Ti-6Al-4V: Multi-objectives optimization. In Proceedings of the 2014 2nd International Conference on Technology, Informatics, Management, Engineering & Environment, Bandung, Indonesia, 19–21 August 2014. [Google Scholar]
  31. Payal, H.; Maheshwari, S.; Bharti, P.S.; Sharma, S.K. Multi-objective optimisation of electrical discharge machining for Inconel 825 using Taguchi-fuzzy approach. Int. J. Inf. Technol. 2019, 11, 97–105. [Google Scholar] [CrossRef]
  32. Modi, D.; Dalsaniya, A.; Fuse, K.; Kanakhara, M. Optimizing the Design of Brake Disc Using Multi-criteria Decision-making Method AHP-TOPSIS for All-terrain Vehicle. In Proceedings of the 2020 9th International Conference on Industrial Technology and Management (ICITM), Oxford, UK, 1–13 February 2020. [Google Scholar]
  33. Ananthakumar, K.; Rajamani, D.; Balasubramanian, E.; Davim, J.P. Measurement and optimization of multi-response characteristics in plasma arc cutting of Monel 400™ using RSM and TOPSIS. Measurement 2019, 135, 725–737. [Google Scholar] [CrossRef]
  34. Prabhu, S.R.; Shettigar, A.; Herbert, M.; Rao, S. Multi response optimization of friction stir welding process variables using TOPSIS approach. IOP Conf. Ser. Mater. Sci. Eng. 2018, 376, 012134. [Google Scholar] [CrossRef] [Green Version]
  35. Tamjidy, M.; Baharudin, B.; Paslar, S.; Matori, K.; Sulaiman, S.; Fadaeifard, F. Multi-objective optimization of friction stir welding process parameters of AA6061-T6 and AA7075-T6 using a biogeography based optimization algorithm. Materials 2017, 10, 533. [Google Scholar] [CrossRef] [Green Version]
  36. Sudhagar, S.; Sakthivel, M.; Mathew, P.J.; Daniel, S.A.A. A multi criteria decision making approach for process improvement in friction stir welding of aluminium alloy. Measurement 2017, 108, 1–8. [Google Scholar] [CrossRef]
  37. Gaidhani, Y.; Kalamani, V. Abrasive water jet review and parameter selection by AHP method. IOSR J. Mech. Civil. Eng. 2013, 8, 1–6. [Google Scholar] [CrossRef]
  38. Babu, K.A.; Venkataramaiah, P. Multi-response optimization in wire electrical discharge machining (WEDM) of Al6061/SiCp composite using hybrid approach. J. Manuf. Sci. Prod. 2015, 15, 327–338. [Google Scholar]
  39. Nayak, B.; Mahapatra, S. Multi-response optimization of WEDM process parameters using the AHP and TOPSIS method. Int. J. Theor. Appl. Res. Mech. Eng. 2013, 2, 109–215. [Google Scholar]
  40. Chou, Y.-C.; Yen, H.-Y.; Dang, V.T.; Sun, C.-C. Assessing the Human Resource in Science and Technology for Asian Countries: Application of Fuzzy AHP and Fuzzy TOPSIS. Symmetry 2019, 11, 251. [Google Scholar] [CrossRef] [Green Version]
  41. Sirisawat, P.; Kiatcharoenpol, T. Fuzzy AHP-TOPSIS approaches to prioritizing solutions for reverse logistics barriers. Comput. Ind. Eng. 2018, 117, 303–318. [Google Scholar] [CrossRef]
  42. Roy, T.; Dutta, R.K. Integrated fuzzy AHP and fuzzy TOPSIS methods for multi-objective optimization of electro discharge machining process. Soft Comput. 2019, 23, 5053–5063. [Google Scholar] [CrossRef]
  43. Saaty, T.L. A scaling method for priorities in hierarchical structures. J. Math. Psychol. 1977, 15, 234–281. [Google Scholar] [CrossRef]
  44. Wijitkosum, S.; Sriburi, T. Fuzzy AHP integrated with GIS analyses for drought risk assessment: A case study from upper Phetchaburi River basin, Thailand. Water 2019, 11, 939. [Google Scholar] [CrossRef] [Green Version]
  45. Leśniak, A.; Kubek, D.; Plebankiewicz, E.; Zima, K.; Belniak, S. Fuzzy AHP application for supporting contractors’ bidding decision. Symmetry 2018, 10, 642. [Google Scholar] [CrossRef] [Green Version]
  46. Chen, S.-J.; Hwang, C.-L. Fuzzy Multiple Attribute Decision Making Methods. Fuzzy Multiple Attribute Decision Making; Springer: Berlin, Germany, 1992; pp. 289–486. [Google Scholar]
  47. Nădăban, S.; Dzitac, S.; Dzitac, I. Fuzzy TOPSIS: A General View. Procedia Comput. Sci. 2016, 91, 823–831. [Google Scholar] [CrossRef] [Green Version]
  48. Chen, C.-T.; Lin, C.-T.; Huang, S.-F. A fuzzy approach for supplier evaluation and selection in supply chain management. Int. J. Prod. Econ. 2006, 102, 289–301. [Google Scholar] [CrossRef]
  49. Torfi, F.; Farahani, R.Z.; Rezapour, S. Fuzzy AHP to determine the relative weights of evaluation criteria and Fuzzy TOPSIS to rank the alternatives. Appl. Soft Comput. 2010, 10, 520–528. [Google Scholar] [CrossRef]
  50. Chaudhari, R.; Vora, J.J.; Patel, V.; López de Lacalle, L.; Parikh, D.J.M. Surface analysis of wire-electrical-discharge-machining-processed shape-memory alloys. Materials 2020, 13, 530. [Google Scholar] [CrossRef] [Green Version]
  51. Pannerselvam, R. Design and Analysis of Experiments; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2012. [Google Scholar]
  52. Wankhede, V.; Jagetiya, D.; Joshi, A.; Chaudhari, R. Experimental investigation of FDM process parameters using Taguchi analysis. Mater. Today Proc. 2020, 27, 2117–2120. [Google Scholar] [CrossRef]
  53. Chaurasia, A.; Wankhede, V.; Chaudhari, R. Experimental Investigation of High-Speed Turning of INCONEL 718 Using PVD-Coated Carbide Tool under Wet Condition. Innovations in Infrastructure; Springer: Singapore, 2019; pp. 367–374. [Google Scholar]
  54. Chaudhari, R.; Vora, J.; Lacalle, L.; Khanna, S.; Patel, V.K.; Ayesta, I. Parametric Optimization and Effect of Nano-Graphene Mixed Dielectric Fluid on Performance of Wire Electrical Discharge Machining Process of Ni55. 8Ti Shape Memory Alloy. Materials 2021, 14, 2533. [Google Scholar] [CrossRef] [PubMed]
  55. Chaudhari, R.; Vora, J.J.; Patel, V.; Lacalle, L.; Parikh, D.J.M. Effect of WEDM process parameters on surface morphology of nitinol shape memory alloy. Materials 2020, 13, 4943. [Google Scholar] [CrossRef]
  56. Thirumalai, R.; Seenivasan, M.; Panneerselvam, K. Experimental investigation and multi response optimization of turning process parameters for Inconel 718 using TOPSIS approach. Mater. Today Proc. 2021, 45, 467–472. [Google Scholar] [CrossRef]
  57. Gegovska, T.; Koker, R.; Cakar, T. Green Supplier Selection Using Fuzzy Multiple-Criteria Decision-Making Methods and Artificial Neural Networks. Comput. Intell. Neurosci. 2020, 26, 8811834. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Membership function fuzzy scale for AHP along with linguistic variables [44].
Figure 1. Membership function fuzzy scale for AHP along with linguistic variables [44].
Materials 14 07408 g001
Figure 2. Triangular fuzzy number [45].
Figure 2. Triangular fuzzy number [45].
Materials 14 07408 g002
Figure 3. Residual plot for cutting speed.
Figure 3. Residual plot for cutting speed.
Materials 14 07408 g003
Figure 4. Main effect plot of cutting speed.
Figure 4. Main effect plot of cutting speed.
Materials 14 07408 g004
Figure 5. Residual plot for MRR.
Figure 5. Residual plot for MRR.
Materials 14 07408 g005
Figure 6. Main effect plot for MRR.
Figure 6. Main effect plot for MRR.
Materials 14 07408 g006
Figure 7. Residual plot for SR.
Figure 7. Residual plot for SR.
Materials 14 07408 g007
Figure 8. Main effect plot for SR.
Figure 8. Main effect plot for SR.
Materials 14 07408 g008
Figure 9. Main effect plot of S/N ratio of CCi.
Figure 9. Main effect plot of S/N ratio of CCi.
Materials 14 07408 g009
Table 1. Chemical composition (wt.%) of Ti6Al4V.
Table 1. Chemical composition (wt.%) of Ti6Al4V.
CFeAlN2CuVTi
0.050.206.200.040.0014.0Balanced
Table 2. Input parameters with working range and their levels.
Table 2. Input parameters with working range and their levels.
ParameterSymbolUnitLevel 1Level 2Level 3
Pulse-on time (Ton)Aµs4070100
Pulse-off time (Toff)Bµs152025
CurrentCA234
Table 3. Input parameters using a central composite design of RSM.
Table 3. Input parameters using a central composite design of RSM.
Std. OrderRun OrderTonToffCurrent
141223
32131
93122
154222
115212
186222
67313
138221
179222
1210232
1611222
512113
113111
1914222
1015322
716133
417331
818333
219311
2020222
Table 4. Membership function [44].
Table 4. Membership function [44].
Fuzzy NumberLinguistic ScaleFuzzy Number
9Perfect8910
8Absolute789
7Very good678
6Fairly good567
5Good456
4Preferable345
3Not bad234
2Weak advantage123
1Equal111
Table 5. Transformation of fuzzy membership function [49].
Table 5. Transformation of fuzzy membership function [49].
RankSub-Criteria GradeMembership Function
Very Low (VL)01(0.00, 0.10, 0.25)
Low (L)02(0.15, 0.30, 0.45)
Medium (M)03(0.35, 0.50, 0.65)
High (H)04(0.55, 0.70, 0.85)
Very High (VH)05(0.75, 090, 1.00)
Table 6. Experimental results and normalized values of the output responses.
Table 6. Experimental results and normalized values of the output responses.
Run OrderTon
(µs)
Toff
(µs)
Current
(A)
Experimental ValuesNormalized Values
VC
(mm/min)
MRR
(mm3/min)
SR
(µm)
VCMRRSR
1702042.7153.7305.800.66320.75970.5498
2402521.2341.5803.270.00000.00000.0249
3402032.0122.5603.150.34840.34630.0000
4702032.3603.0005.550.50430.50180.4979
5701532.9173.7105.890.75370.75270.5685
6702032.4413.1095.200.54050.54030.4253
71001543.4674.4107.971.00001.00001.0000
8702021.7792.2604.010.24410.24030.1784
9702032.4443.1105.330.54190.54060.4523
10702532.0332.5805.220.35780.35340.4295
11702032.4773.1505.830.55670.55480.5560
12401543.0133.8303.960.79670.79510.1680
13401522.1142.6902.980.39410.39220.1037
14702032.4863.1705.000.56070.56180.3838
151002032.7313.5006.100.67040.67840.6120
16402541.8902.4004.200.29380.28980.2178
171002521.6732.1004.830.19660.18370.3485
181002542.4903.1705.700.56250.56180.5290
191001522.3813.0804.710.51370.53000.3237
20702032.4773.2105.600.55670.57600.5083
Table 7. ANOVA for cutting speed.
Table 7. ANOVA for cutting speed.
SourceSum of
Squares
DfMean Sum of
Square
F Valuep-ValueContributionSignificance
Model4.8387590.53764112.300.00099.02%significant
Ton0.6145410.61454128.370.00012.58%significant
Toff2.0903212.09032436.630.00042.78%significant
Current1.9307211.93072403.290.00039.51%significant
Ton × Toff0.0126410.012642.640.1350.26%
Ton × Current0.0151410.015143.160.1060.31%
Toff × Current 0.0327710.032776.840.0260.67%significant
Ton × Ton0.0056610.005661.180.3020.12%
Toff × Toff0.0092910.009291.940.1940.19%
Current × Current0.0793510.0793516.570.0021.62%significant
Residual0.04787100.00479 0.98%
Lack of Fit0.0369450.007393.380.1040.76%Insignificant
Pure Error0.0109350.00219 0.22%
Total4.8866219 100.00%
Table 8. ANOVA for MRR.
Table 8. ANOVA for MRR.
SourceSum of
Squares
Degree of FreedomAdjusted Mean Sum of
Square
F Valuep-ValueContributionSignificance
Model8.1708490.9078763.960.00098.29%significant
Ton1.0240011.0240072.140.00012.32%significant
Toff3.4692113.46921244.390.00041.73%significant
Current3.3988913.39889239.440.00040.89%significant
Ton × Toff0.0128010.012800.900.3650.15%
Ton × Current0.0242010.024201.700.2210.29%
Toff × Current 0.0420510.042052.960.1160.51%
Ton × Ton0.0272310.027231.920.1960.33%
Toff × Toff0.0006610.000660.050.8340.01%
Current × Current0.0497510.049753.500.0910.60%
Residual0.14195100.01420 1.71%
Lack of Fit0.1159750.023194.460.06261.40%Insignificant
Pure Error0.0259850.00520 0.31%
Total8.3127919 100.00%
Table 9. ANOVA for SR.
Table 9. ANOVA for SR.
SourceSum of
Squares
Degree of FreedomMean Sum of
Square
F Valuep-ValueContributionSignificance
Model22.377692.486411.670.00091.31%significant
Ton12.2766112.276657.650.00050.09%significant
Toff0.876210.87624.110.0703.58%Not significant
Current5.126615.126624.070.00120.92%significant
Ton × Toff0.505010.50502.370.1552.06%
Ton × Current1.044011.04404.900.0514.26%
Toff × Current 0.391610.39161.840.2051.60%significant
Ton × Ton0.982510.98254.610.0574.01%
Toff × Toff0.303610.30361.430.2601.24%
Current × Current0.277610.27761.300.2801.13%significant
Residual2.1297100.2130 8.69%
Lack of Fit1.679450.33593.730.0876.85%Insignificant
Pure Error0.450350.0901 1.84%
Total24.507319 100.00%
Table 10. Model summary for VC, MRR, and SR.
Table 10. Model summary for VC, MRR, and SR.
ResponseUnitStandard DeviationR-sqR-sq (adj)
VCmm/min0.069191299.02%98.14%
MRRmm3/min0.11914498.29%96.76%
SRµm0.46148591.31%83.49%
Table 11. Optimum parameter setting considering single objective optimization.
Table 11. Optimum parameter setting considering single objective optimization.
ResponseUnitOptimum Parameter Setting Considering Single Objective Optimization
VCmm/minA3B1C3
MRRmm3/minA3B1C3
SRµmA1B3C1
Table 12. Comparison matrix.
Table 12. Comparison matrix.
VCMRRSR
VC11/31/7
MRR311/4
SR741
Table 13. Fuzzified comparison matrix.
Table 13. Fuzzified comparison matrix.
VCMRRSR
VC(1, 1, 1)(0.25, 0.33, 0.50)(0.13, 0.14, 0.17)
MRR(2, 3, 4)(1, 1, 1)(0.2, 0.25, 0.33)
SR(6, 7, 8)(3, 4, 5)(1, 1, 1)
Table 14. Fuzzy weights.
Table 14. Fuzzy weights.
Weights
VC(0.5286, 0.7049, 0.9312)
MRR(0.1486, 0.2109, 0.2996)
SR(0.0635, 0.0841, 0.1189)
Table 15. Fuzzified normalized data of the output responses.
Table 15. Fuzzified normalized data of the output responses.
AlternativesVC (mm/min)MRR (mm3/min)SR (µm)
1(0.55, 0.70, 0.85)(0.55, 0.70, 0.85)(0.35, 0.50, 0.65)
2(0.00, 0.10, 0.25)(0.00, 0.10, 0.25)(0.00, 0.10, 0.25)
3(0.15, 0.30, 0.45)(0.15, 0.30, 0.45)(0.00, 0.10, 0.25)
4(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)
5(0.55, 0.70, 0.85)(0.55, 0.70, 0.85)(0.35, 0.50, 0.65)
6(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)
7(0.75, 0.90, 1.0)(0.75, 0.90, 1.00)(0.75, 0.90, 1.00)
8(0.15, 0.30, 0.45)(0.15, 0.30, 0.45)(0.00, 0.10, 0.25)
9(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)
10(0.15, 0.30, 0.45)(0.15, 0.30, 0.45)(0.35, 0.50, 0.65)
11(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)
12(0.55, 0.70, 0.85)(0.55, 0.70, 00.85)(0.00, 0.10, 0.25)
13(0.15, 0.30, 0.45)(0.15, 0.30, 0.45)(0.00, 0.10, 0.25)
14(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)(0.15, 0.30, 0.45)
15(0.55, 0.70, 0.85)(0.55, 0.70, 0.85)(0.55, 0.70, 0.85)
16(0.15, 0.30, 0.45)(0.15, 0.30, 0.45)(0.15, 0.30, 0.45)
17(0.00, 0.10, 0.25)(0.00, 0.10, 0.25)(0.15, 0.30, 0.45)
18(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)
19(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)(0.15, 0.30, 0.45)
20(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)(0.35, 0.50, 0.65)
Table 16. Weighted normalized data.
Table 16. Weighted normalized data.
AlternativesVC (mm/min)MRR (mm3/min)SR (µm)
1(0.034, 0.059, 0.101)(0.082, 0.148, 0.255)(0.185, 0.352, 0.605)
2(0.000, 0.008, 0.030)(0.000, 0.021, 0.075)(0.000, 0.070, 0.233)
3(0.010, 0.025, 0.054)(0.022, 0.063, 0.135)(0.000, 0.070, 0.233)
4(0.022, 0.042, 0.077)(0.052, 0.105, 0.195)(0.185, 0.352, 0.605)
5(0.035, 0.059, 0.101)(0.082, 0.148, 0.255)(0.185, 0.352, 0.605)
6(0.022, 0.042, 0.077)(0.052, 0.105, 0.195)(0.185, 0.352, 0.605)
7(0.048, 0.076, 0.119)(0.111, 0.190, 0.300)(0.396, 0.634, 0.931)
8(0.010, 0.025, 0.054)(0.022, 0.063, 0.135)(0.000, 0.070, 0.233)
9(0.022, 0.042, 0.077)(0.052, 0.105, 0.195)(0.185, 0.352, 0.605)
10(0.010, 0.025, 0.054)(0.022, 0.063, 0.135)(0.185, 0.352, 0.605)
11(0.022, 0.042, 0.077)(0.052, 0.105, 0.195)(0.185, 0.352, 0.605)
12(0.035, 0.059, 0.101)(0.082, 0.148, 0.255)(0.000, 0.070, 0.233)
13(0.010, 0.025, 0.054)(0.022, 0.063, 0.135)(0.000, 0.070, 0.233)
14(0.022, 0.042, 0.077)(0.052, 0.105, 0.195)(0.079, 0.211, 0.419)
15(0.035, 0.059, 0.101)(0.022, 0.063, 0.134)(0.291, 0.493, 0.792)
16(0.010, 0.025, 0.054)(0.022, 0.063, 0.135)(0.079, 0.211, 0.419)
17(0.000, 0.008, 0.030)(0.000, 0.021, 0.075)(0.079, 0.211, 0.419)
18(0.022, 0.042, 0.077)(0.052, 0.105, 0.195)(0.185, 0.352, 0.605)
19(0.022, 0.042, 0.077)(0.052, 0.105, 0.195)(0.079, 0.211, 0.419)
20(0.022, 0.042, 0.077)(0.052, 0.105, 0.195)(0.185, 0.352, 0.605)
Table 17. Closeness coefficient index.
Table 17. Closeness coefficient index.
AlternativesD+D–CCi
12.1950.8900.456
22.0960.9670.477
32.0391.0260.490
42.2560.8270.443
52.1950.8900.456
62.2560.8270.443
72.4130.7100.414
82.0391.0260.490
92.2560.8270.443
102.3170.7640.432
112.2560.8270.443
121.9181.1510.522
132.0391.0260.490
142.1120.9600.473
152.3410.7640.427
162.1730.8980.460
172.2310.8390.448
182.2560.8270.443
192.1120.9600.473
202.2560.8270.443
Table 18. Results of the confirmation experiment.
Table 18. Results of the confirmation experiment.
Performance ResponseOptimal SettingPredicted ValuesExperimental Values% Error
VC (mm/min)Ton 40 µs, Toff 15 µs, Current 2A2.0672.1142.22
MRR (mm3/min)2.6162.6902.75
SR (µm)3.1172.984.39
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14237408

AMA Style

Fuse K, Dalsaniya A, Modi D, Vora J, Pimenov DY, 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(23):7408. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14237408

Chicago/Turabian Style

Fuse, Kishan, Arrown Dalsaniya, Dhananj Modi, Jay Vora, Danil Yurievich Pimenov, Khaled Giasin, Parth Prajapati, Rakesh Chaudhari, and Szymon Wojciechowski. 2021. "Integration of Fuzzy AHP and Fuzzy TOPSIS Methods for Wire Electric Discharge Machining of Titanium (Ti6Al4V) Alloy Using RSM" Materials 14, no. 23: 7408. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14237408

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