Next Article in Journal
Optical Properties of Au-Based and Pt-Based Alloys for Infrared Device Applications: A Combined First Principle and Electromagnetic Simulation Study
Next Article in Special Issue
Charged Satellite Drop Avoidance in Electrohydrodynamic Dripping
Previous Article in Journal
A Novel Fabricating Process of Catalytic Gas Sensor Based on Droplet Generating Technology
Previous Article in Special Issue
Investigation of Product and Process Fingerprints for Fast Quality Assurance in Injection Molding of Micro-Structured Components
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Response Optimization of Electrical Discharge Machining Using the Desirability Function †

by
Rafał Świercz
*,
Dorota Oniszczuk-Świercz
and
Tomasz Chmielewski
Institute of Manufacturing Technology, Warsaw University of Technology, 00-661 Warsaw, Poland
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in the XIII Internation al Conference Electromachining 2018 that was held in Bydgoszcz, Poland, 9–11 May 2018.
Micromachines 2019, 10(1), 72; https://doi.org/10.3390/mi10010072
Submission received: 16 November 2018 / Revised: 9 January 2019 / Accepted: 14 January 2019 / Published: 20 January 2019
(This article belongs to the Special Issue Product/Process Fingerprint in Micro Manufacturing)

Abstract

:
Electrical discharge machining (EDM) is a modern technology that is widely used in the production of difficult to cut conductive materials. The basic problem of EDM is the stochastic nature of electrical discharges. The optimal selection of machining parameters to achieve micron surface roughness and the recast layer with the maximal possible value of the material removal rate (MRR) is quite challenging. In this paper, we performed an analytical and experimental investigation of the influence of the EDM parameters: Surface integrity and MRR. Response surface methodology (RSM) was used to build empirical models on the influence of the discharge current I, pulse time ton, and the time interval toff, on the surface roughness (Sa), the thickness of the white layer (WL), and the MRR, during the machining of tool steel 55NiCrMoV7. The surface and subsurface integrity were evaluated using an optical microscope and a scanning profilometer. Analysis of variance (ANOVA) was used to establish the statistical significance parameters. The calculated contribution indicated that the discharge current had the most influence (over the 50%) on the Sa, WL, and MRR, followed by the discharge time. The multi-response optimization was carried out using the desirability function for the three cases of EDM: Finishing, semi-finishing, and roughing. The confirmation test showed that maximal errors between the predicted and the obtained values did not exceed 6%.

1. Introduction

Electrical discharge machining is a precision method of manufacturing hard, complex shaped, conductive materials. The removal mechanism of the material in EDM is the result of the electrical discharge, which causes melting and evaporation in the local surface layers of both the workpiece and the working electrode. Owing to the impact of the thermal and chemical processes, the electrical discharge properties of the surface layer of the material are changed [1,2,3]. Craters form a specific surface texture. Surface roughness parameters directly depend on the discharge current and the time pulse [4]. Heat flux changes the surface integrity. Metallographic images of the machined samples of tool steel show new layers, i.e., the external melting layer (white layer), the heat-affected zones, and the tempered layer. The white layer is non-homogeneous and discontinuous and may vary in thickness on the analyzed sample. Discontinuity of the melted layer is caused by the random occurrence of electrical discharge and the machining conditions. This layer is characterized by high variations in thickness and microstructure defects of the material, such as micro-cracks. Microcracks are an undesirable effect, resulting in reduced fatigue resistance and corrosion resistance. It is important to choose the appropriate parameters and processing conditions to obtain the smallest thickness of the white layer or its complete elimination. The quality of the surface after the EDM process does not always meet expectations [5,6]. Therefore, additional technological operations are used to change the surface integrity. Some of the most applicable operations are electrochemical machining [7,8], laser surface modification [9,10], applying coatings [11,12], or the use of hybrid machining [13,14] or non-conventional finishing [15,16,17,18]. However, the use of additional technological operations significantly increases the production costs. Therefore, this work analyzed the impact of the following parameters of the EDM process: discharge current, pulse time, and pulse interval. Furthermore, the optimization of EDM using the desirability function will allow for the selection of the most favorable processing conditions, which reduce the use of additive treatments to the bare minimum.
Electrical discharge machining belongs to a group of non-conventional manufacturing techniques. The material is removed from the workpiece using electrical discharges occurring between the working electrode and the workpiece. Owing to the conducted thermal energy from the electrical discharge, the local temperature increases (in the range of 8000–12,000 °C), leading to the melting and evaporation of a small volume of the surface workpiece and the working electrode. Then, a collapsing plasma channel at the end of the discharge induces high-pressure waves that rinse the molten and evaporated metal [19]. The physics of the material removal phenomenon in EDM is complex. A model of the EDM process proposed by Izquierdo et al. [20], shows that by using the superposition of multiple discharges and calculating the temperature fields inside the workpiece, it is possible to predict the surface roughness parameters and the material removal rate. One of the main problems in modeling is the appropriate determination of the influence of EDM parameters like discharge voltage, discharge current, and pulse time on ionization and growth of the plasma channel. Information about the percentage of discharge energy devoted to the heat flux, the mechanism of plasma channel growth, and the temperature which facilitates material removal, is used to model not only surface roughness, but also to build models of structural changes in the surface layers and their thickness [21]. Ming et al. [22] analyzed the distributions of the energy of the workpiece for different materials like Al 6061, Inconel 718, and SKD11. The authors indicated that the analyses of energy efficiency simultaneously with MRR, could be used in existing thermal–physical models to improve the technical performance of the models. Theoretical models of EDM are significant because the models analyze the processes and their implication on the manufactured material. However, the application of the theoretical model in the manufacturing process is not always possible. Given the complexity in describing the physical phenomena of EDM, such as the random distribution of the electric field and the temperature field, the formation of the plasma channel, and changes in the properties of the dielectric, mathematical models have been built based on the empirical studies. Several process variables promote the application of the optimization method to achieve high productivity.
Optimization of the EDM machining process can be carried out using various methods like response surface methodology [23,24], Taguchi analyses [25], artificial networks [26], grey-based response surface methodology [27,28], and the Deringer desirable [29] or hybrid methods [30,31]. One of the most common manufacturing processes is the response surface methodology (RSM). RSM is extensively used in an analytical and industrial application like turning [32,33], milling [34], welding [35], grinding [36], and erosion machining [37]. Ghodsiyeh et al. [38] used RSM to optimize the wire electrical discharge machining of the titanium alloy Ti-6Al-4 V. Presented results indicated that the pulse time had the highest impact on surface roughness, whilst discharge current had the same role for the wire wear ratio and white layer thickness. Alavi et al. [39] analyzed the influence of the EDM parameters on tool wear, crater size, and microhardness on the titanium alloy Ti-6Al-4V. The presented research indicates that the main effect on the crater size was the discharge voltage, whilst the capacitance was the most important for tool wear and surface micro-hardness. Increasing the capacitance caused a reduction in tool wear and an increase in the micro-hardness of the machined surface. Selvarajan et al. [40] investigated the possibilities of the EDM composite Si3N4-TiN. The authors used RSM to determine the optimal parameters for the machining of ceramic composites. The conducted research showed that the MRR and surface roughness of the manufactured parts depended on the value of the discharge current and the pulse time, and the results were consistent with the results for the processing of tool steel [41]. Kumaran et al. [42] used a grey fuzzy logic approach to optimize the EDM parameters during the machining of carbon fiber reinforced plastic composite. Their research showed that established optimal parameters in ultrasonic-assisted EDM allowed improvements in the deburring rate, with a simultaneous improvement in the tool wear rate (TWR). Gu et al. [43] indicated that the machining of new alloys, which have a high melting point and good thermal conductivity, like titanium–zirconium–molybdenum, require optimizations of the EDM in connection with the analysis. The presented results showed that the crater diameter was much smaller than the plasma-affected zone. To improve the machining performance (Ra, MRR), the response surface methodology was used. Dang [44] proposed the Kriging regression model and particle swarm in the optimization of EDM of P20 steel. The authors indicated that the Kriging model could capture the nonlinear characteristics and was better able to obtain the optimum parameters for the MRR, tool wear, and surface roughness. Mohanty et al. [45] pointed out that the choice of the electrode material should be considered in the optimization of the EDM process. The authors found that in the machining of Inconel 718, the material removal rate and tool wear could be improved by using a graphite electrode and to improve surface integrity, the better choice was a brass and copper electrode. Using the utility concept and the quantum particle swarm optimization (QPSO) algorithm, the optimal parametric setting was developed with the objectives to maximize the MRR and minimize tool wear, surface roughness, and radial overcut. Research presented by Maity et al. [46] on the influence of EDM parameters on the thickness of the recast layer, material removal rate, and overcut on the machining of Inconel 718, showed that the optimization parameters using the RSM and Artificial Bee Colony algorithm gave an average prediction accuracy of about 3.5% in relation to confirmation tests. The predictive efficiency of neural networks may be affected by different factors like noise corruption, spatial distribution, and the size of the data used to construct the artificial neural network (ANN) model. Tripathy et al. [47], in order to optimize the powder mixed (SiC) electro-discharge machining of H-11 die steel, which seeks to maximize the MRR and minimize electrode wear and surface roughness, used a different method of optimization. The authors used the grey relational analysis and the technique for the order of preference by similarity, where the Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) solution achieved a similar effect of an improved performance of the quality characteristics. Nguyen et al. [48] investigated the influence of powder mixed (Ti) electro-discharge machining of SKD61, SKD11, SKT4 steel on the surface roughness (SR), MRR, and microhardness. The presented results showed that the addition of Ti powder in dielectric resulted in reduced SR and increased microhardness. The authors indicated that optimization with the Taguchi–TOPSIS made it difficult to select the optimal parameters. The presented research showed that the measured distance could lead to confusion in selecting the best alternative. A fuzzy analytic hierarchy process (AHP) and fuzzy TOPSIS method were used by Roy et al. [49] to optimize multiple responses of the material removal rate, tool wear, and tool overcut in EDM based on various process parameters. Kandpal et al. [50] investigated the influence of EDM parameters on the MRR, tool wear, and overcut of aluminum matrix composites. The authors showed optimization with the utility concept, which provided the collective optimization of both responses for improving the mean of the process. D’Urso et al. [51] proposed the optimization of EDM micro drilling using a cost index which combined two opposite effects of the material removal rate and tool wear. The minimization of the cost index enabled optimal working conditions. Parsana et al. [52] indicated that in the optimization process of EDM drilling, an important point to consider was the roundness of the holes. Using the RSM and passing vehicle search algorithm, normalized weights proved to be useful in obtaining the Pareto fronts for a combination of different objectives at a time. Research carried out by Hadad et al. [53] showed that the analysis of the optimal parameters of EDM machining, in addition to the electrical parameters, should also include the initial roughness of the working electrodes, which has significant effects on the machining performance during the finishing, semi-rough, and rough EDM processes.
The published literature indicates that few studies have reported on the optimization of the EDM, which considers the surface 3D roughness parameter, white layer thickness, and the MRR in the three stages of machining: roughing, semi-finishing, and finishing, with the desirable function. Therefore, in this paper, a multi-response optimization of the EDM of tool steel 55NiCrMoV7 was conducted. This material has a wide range of industry applications on die matrices, matrix inserts, and hydraulic and mechanical press dies. The surface roughness was investigated using the 3D roughness parameter. The parameter, Sa, gives more information about the surface properties. In electrical discharge machining, the surface roughness is obtained by overlapping the craters of individual discharges. EDM with parameters corresponding to the roughing and semi-finishing operations results in surfaces which have different profiles on the cross-sections of samples. The calculated Ra parameters on the profile cross-sections may provide inaccurate results.
The optimization of parameters was performed using the desirable function. Optimization of electrical discharge machining was divided into three cases: finishing, semi-finishing, and roughing. In each case, different goals were set. For finishing, the goal was to minimize the surface roughness (Sa) simultaneously whilst minimizing the thickness of the white layer (WL), with a possible maximized MRR. In the case of semi-finishing, the EDM goal was to obtain a specific value of the MRR whilst possibly minimizing the roughness (Sa), simultaneously minimizing the thickness of the white layer (WL). In the last case of roughing, the goal was to maximize the material removal rate with the possibility of minimizing the roughness (Sa) and the thickness of the white layer (WL).
The topic of the article focuses on describing the changes occurring in the material as a result of the local thermal processes due to electric discharges. Experimental studies allowed for a better understanding of the relationship between the changes in surface integrity of tool steel 55NiCrMoV7 and how to optimize the process to achieve a micron roughness and recast layer, with the possibility of achieving the maximal value of the material removal rate.

2. Materials and Methods

Industry applications of electrical discharge machining are limited by the obtained specific surface integrity and low material removal rate. The purpose of the research was to develop the multi-response optimization of the EDM process of the tool steel 55NiCrMoV7 for three cases: finishing, semi-finishing, and roughing. Tool steel 55NiCrMoV7 was chosen because of its wide industry applications on die matrices, matrix inserts, and hydraulic and mechanical press dies. This material is characterized by high dimensional stability and crack resistance, with dynamically changing pressures and rapid heating and cooling during operation. Heat-treated samples of the tool steel (55 HRC) had the dimensions of 12 mm × 12 mm × 3 mm. Experimental studies were conducted on the electrical discharge machine, Charmilles Form 2LC ZNC (GF Solutions, Geneva, Switzerland). The electrode used was graphite (EDM-3 POCO), and the EDM fluid 108 MP-SE 60 was used as the dielectric. The present paper was focused on the selection of optimal parameters for EDM, which led to minimum surface roughness, as well as the thickness of the white layer and maximum productivity.
The main object of the study was the optimization of the EDM process using statistical models on the influence of EDM parameters on surface roughness (Sa), the thickness of the white layer, and the material removal rate. To achieve this goal, experimental research was carried out using a completely orthogonal design of the experiment, three-level three parameters full factorial design. Choice of this type of experiment design allowed the reduction in the number of experimental runs required to generate sufficient information for a statistically adequate result. A schematic diagram of the experimental set-up is shown in Figure 1. Investigation of the surface roughness parameters after the EDM was carried out on a Taylor–Hobson FORM TALYSURF Series 2 scan profilometer (Taylor Hobson, Leicester, United Kingdom). The roughness parameter (Sa), the arithmetic mean of the deviations from the mean, was measured on a surface area of 2 mm × 4 mm with a discretization step (10 μm) in the X-axis and Y-axis. The Sa (average value of the absolute heights over the entire surface) parameter responded to the 2D roughness profile parameters Ra. This may be obtained by adding the individual height values, without regard to sign, and dividing the sum by the number of the data matrix, where M is the number of points per profile, N is the number of profiles, and z, x, y are the heights of the profile at a specific point.
S a = 1 N M x = 0 N = 1 y = 0 M = 1 | z x , y |
Metallographic surface structure studies were performed using a Nikon Eclipse LV 150 optical microscope (Nikon, Tokyo, Japan), coupled to an NIS-Elements BR 3.0 image analyzer (Nikon). Specimens were included in the resin, and were then machined with grinding and polishing. Micro-etching was performed with nital (5%) to reveal the microstructure of the material. The maximum thickness of the white layer in sections was measured for each sample.

2.1. Uncertainty Evaluation Procedure

To verify the quality of the measurements, an uncertainty evaluation was carried out. The measurements were carried out inside a metrological laboratory with a 20 °C ± 0.5 °C controlled temperature. The thermal deformation of the samples was neglected. The surface topography measurements were carried out using a Taylor–Hobson FORM TALYSURF Series 2 scan profilometer (Taylor Hobson, Leicester, United Kingdom). The raw surface acquisitions were post-processed with the dedicated image metrology software TalyMap (Taylor Hobson, Leicester, United Kingdom). The calibration of the scanning profilometer was performed with a calibrated roughness artifact (nominal value: Ra = 810 nm). According to Reference [54], it is possible to establish uncertainties for surface roughness Sa measurements following ISO 15530-3 [55].
The uncertainties for surface roughness measurements when using the induction scanning profilometer UPROF was calculated according to the following equation:
U PROF = u 2 cal + u 2 p   + u 2 res ,   PROF  
where ucal is the standard calibration uncertainty of the roughness standard; up is the standard uncertainty related to the measurement procedure and is calculated as a standard deviation of ten repeated measurements on the calibrated standard; and ures,PROF is the resolution standard uncertainty related to the declared 3 nm vertical resolution of the scanning profilometer for measuring range 0.2 mm.
The expanded uncertainty of the surface roughness Sa measurement was calculated as follows:
U 95 , Sa = k   U 2 PROF + U 2 S a , EDM
where k is the coverage factor, equal to 2 for a 95% confidence interval, uSa,EDM is calculated using the standard deviation of repeated Sa measurements on the electrical discharge machining sample. Table 1 presents the uncertainty budget.
The measurements of the thickness of the white layer were carried out using a Nikon Eclipse LV 150 optical microscope, coupled to an NIS-Elements BR 3.0 image analyzer (Nikon). The uncertainty was calculated using the above method. The calibrated slide with 10 µm division was selected as the calibrated artifact.
The uncertainties for the thickness of the white layer when using optical microscope UOM was calculated according to the following equation:
U OM = u 2 cal + u 2 p   + u 2 res ,   OM
where ucal is the standard calibration uncertainty of the calibration slide, up is the standard uncertainty related to the measurement procedure and is calculated as the standard deviation of ten repeated measurements on the calibrated slide; and ures,OM is the resolution standard uncertainty related to the objective magnification (50x). The expanded uncertainty of the thickness of the white layer measurement was calculated as follows:
U 95 , WL = k   U OM  
where k is the coverage factor, equal to 2 for a 95% confidence interval. Table 2 presents the uncertainty budget.
The material removal rate (MRR) was calculated based on the volume of material removed from the workpiece divided by the machining time:
MRR = m 1 m 2 ρ   Δ t   [ m m 3 m i n ]
where m 1 is the sample weight before processing, m 2 is the sample weight after processing, ρ is specific material density, Δ t is a time of manufacturing.
Each sample was weighed before manufacturing on a precision electronic balance (Radwag, Radom, Poland). The samples after the EDM process were cleansed with the compressed air and then weighed again. The measurement uncertainties for the weight measurements were calculated according to the following equation:
U B = u 2 m 1 + u 2 res + u 2 i   + u 2 ie  
where um1 is the standard uncertainty related to the measurement procedure and is calculated as the standard deviation of the ten repeated measurements of the sample, ures, is the resolution uncertainty related to the declared resolution of the balance with a readability 0.01 mg for the measuring samples of a maximum capacity 50 g, ui is the uncertainty related to a balance indication error, and uie is the uncertainty on a determining indication error.
The expanded uncertainty of the weight measurement was calculated as follows:
U 95 , W = k   U B  
where k is the coverage factor, equal to 2 for a 95% confidence interval. Table 3 reports the uncertainty budget.

2.2. Analyses of Current and Voltage Waveforms.

The complexity of the physical phenomena of the EDM process and its conditions caused considerable difficulties in describing and identifying the impact of individual machining parameters on the surface roughness, integrity, and the MRR. In the first stage of the conducted research, the measurement circuit was developed to determine the current–voltage characteristics of the generator machines. A primary test was conducted to investigate a range of stability discharges for different values of the discharge current, pulse time, and time interval. The measurement of the current and voltage waveforms in the EDM process conditions was done using a National Instruments NI5133 oscilloscope card (National Instruments, Austin, TX, USA). An application was developed in the LabView environment, which enabled the control of the work of the oscilloscope card. The current measurement was done using the indirect method as the voltage drop on the non-inductive current sensor. The maximum value of the voltage drop for the set current values did not exceed 3V, so the signal was fed directly to the oscilloscope card. The measurement of the voltage during the electric discharge was done with the Tektronix probe (Tektronix UK Ltd., Berkshire, UK). The sampling rate was 100 MS/s, 2-Channel registration. Analyses of the obtained data were performed in DIAdem (National Instruments).
Exemplary current and voltage waveforms registered for the investigated EDM are shown in Figure 2. The workpiece was machined at the moment when the supply voltage Uz dropped to the discharge voltage Uc, with the increase of the discharge current I, during the pulse ton. Then, during the toff interval, the conditions in the gap stabilized, and the process was cyclically repeated.
The following parameters can characterize the voltage–current waveforms (Figure 2):
  • I = the height of the peak current during discharging,
  • Uz = open circuit voltage, this is the system voltage when the EDM circuit is in the open state, and the energy has been built up for discharge,
  • Uc = discharge voltage,
  • ton = pulse time, the time required for the current to rise and fall during discharging,
  • toff = time interval, this is the time from the end of one pulse to the beginning of the next pulse with the current.
Analysis of the obtained voltage and current waveforms showed that at the highest adjustable currents, pulse duration, and minimum values of the break times, in most cases, arc discharges or short circuits occurred. For short time intervals, toff, the plasma channel may not be completely deionized, which increases the probability of another discharge being in the same place. Furthermore, the ineffective removal of the products of erosion from the gap causes a reduction in the dielectric resistance and destabilization of the conditions. There is a high probability of a short circuit. The machine's control system resists the phenomena described above by increasing the gap (temporarily raising the electrode), whilst at the same time extending the break time (Figure 3).
The proper operation of the generator control system ensures the energy repeatability of discharges. Therefore, the presented disturbances will not have a significant impact on the quality of the treated surfaces. Nevertheless, unfavorably selected ranges of the set parameters will significantly affect the efficiency of the process [56]. The following are examples of stable U(t), I(t) waveforms, which enabled the selection of the range of variability of parameters used in the experimental research (Figure 4).
Analysis of the recorded voltage and current waveforms enabled the selection of stable parameters in the EDM process for roughing, semi-finishing, and finishing machining. Experimental studies were conducted using the orthogonal full factorial design DOE (design of experiments) methodology: Orthogonal full factorial design. Preliminary experiments were conducted to obtain the stable discharges in all ranges of the design matrix. After analysis of the results, the following machining conditions were selected: discharge current in the range I = 3–14 A, pulse time in the range ton = 10–400 µs, and time interval toff = 10–150 μs, with the following constants :open voltage U0 = 225 V, discharge voltage Uc = 25 V. Table 4 shows the levels of machining parameters carried out in the experimental design.

3. Results and Discussion

3.1. Analysis of Surface Integrity

The primary factor that affected the properties of the machined parts was the surface texture. Experimental investigations showed that parameters such as surface roughness (Sa), directly depended on the applied machining parameters. Surface topography after EDM was the result of the overlapping of craters from single discharges and it had an isotropic structure (Figure 5).
Owing to rapid local thermal processes during electrical discharge machining, phase changes occurred on the surface layer of the workpiece. The analysis of images of the metallographic structure of tool steel 55NiCrMoV7 showed the occurrence of three characteristic sublayers for the whole range of the investigated machining parameters (Figure 6).
An external molten layer (commonly referred to as a white layer), was formed by melting and rapidly solidifying a thin layer of metal not removed from the surface of the crater during an electric discharge. The white layer in its structure may have chemical decompositions from both the core material and the working electrode. The heat-affected zone (HAZ), which is located directly under the melted layer (visible as a light structure), was characterized by an increased hardness around the core material. The last observed layer was the tempered layer, which was visible as a dark streak immediately below the heat-affected zone layer. The thickness of each observed layer depended on the investigated EDM parameters. The white layer, unlike the heat-affected zone and the tempered layer, was characterized by local discontinuities and thickness changes (Figure 7). The industrial application of the results of the conducted experiments requires information on the largest thickness of the white layer. This information allows for the correct selection of machining allowances for semi-finishing and finishing manufacturing.
Electrical discharges caused the local melting and evaporation of material. The thermal processes of removing material and the rapid re-solidification of the molten metal which was not removed from the discharge crater generated thermal stress. Exceeding the maximum tensile strength of the material caused the generation of micro-cracks. Micro-cracks are an undesirable effect, resulting in reduced fatigue resistance and corrosion resistance. In most cases, micro-cracks propagate to the end of the white layer (Figure 8). In rare cases, the propagation of a crack has been observed to penetrate the core of the material. Micro-cracks can be observed directly on the machining surface (Figure 9).

3.2. Response Surface Methodology

Experimental investigation of the influence of the EDM parameters on the surface roughness (Sa), the thickness of the white layer, and the MRR was carried out using response surface methodology. In RSM, the dependence between the desired response and the independent variables can be represented by the following:
Y = f (I,ton, toff) ± ε
where Y is the response; f is the response function; ε is the experimental error. I the discharge current, ton (μs) the pulse time, and toff (μs) the time interval are independent parameters. In the study, the polynomial regression model was chosen to fit the response function to the experimental results.
The experimental investigation was carried out based on the full factor orthogonal experiment design: three-level three-parameter. The study of the influence of the input factors on three equidistant levels of variation allows for the determination of regression equations with a high degree of correlation and a small spread of values. According to the full factor orthogonal design plan, twenty-eight samples, with one additional replication in the center point, were manufactured and measured. Based on the experimental data, an empirical model of the influence of the discharge current I, pulse time ton, and time interval toff was built. The results of the experimental studies are presented in Table 5. The surface roughness (Sa) was in the range of 1.88 µm to 12.7 µm. The maximal thickness of the white layer was in the range of 5.5 µm to 33.5 µm. The material removal rate was in the range 0.1 mm3/min to 29.19 mm3/min. The obtained value of roughness (Sa), the maximal thickness of the white layer, and the MRR corresponded to the finishing and roughing machining.
Analysis of variance (ANOVA) was used to check the significance of each independent variable in the response function. The ANOVA test was conducted at a 5% significant level. The F-value corresponded to a continuous probability distribution. If this probability (Prob > f) value for each factor was less than 0.05, this indicated that the model factor was significant (i.e., at a 95% confidence level). Values of Prob > f higher than 0.05 indicated that a model factor was non-significant.
The ANOVA results for the Sa, white layer thickness, and the MRR are shown in Table 6, Table 7 and Table 8, respectively. Table 6 shows the ANOVA results for surface roughness (Sa). The calculated contribution indicated that the discharge current had the most influence on the surface roughness (Sa) (57.6%). Second, the affecting variable was pulse time (15.7%) and the interaction of the discharge current with pulse time (14.5%). Other variables and their interactions had a significant influence on the surface roughness (Sa), but each of the contributions did not exceed 5%. The ANOVA results presented in Table 7 indicated that the most significant influence on the maximal thickness of the white layer was the discharge current (47.5%), followed by the pulse time (27.8%) and the squared pulse time (9.2%). The contribution of other variables on the Wl was significant but less important. Table 8 presents the ANOVA results for the MRR. The calculated contributions indicated that the discharge current (55.9%) had the most influence on the MRR, followed by the interaction of the discharge current with the pulse time (12.7%) and also the pulse time (11.6%). Other variables and their interactions were significant, but their contributions were smaller and contained in the range (1.6%–5.6%). From the presented ANOVA Table 6, Table 7 and Table 8, calculated Fisher coefficients for the models Sa, WL, and MRR were 150.95, 141.26, and 208.65, respectively. The results implied that all the developed models were significant at a 95% confidence level.
Regression analysis with a backward elimination process was performed. For each equation, we calculated the coefficient of determination, R-squared, and the adjusted coefficient of determination, R-Adj. The coefficients represented the percentage of variance explained by the model. When the value of the R-sqr and R-Adj approaches unity, a more accurate fit of the regression equation to the research results would be obtained.
After eliminating the non-significant factors in the response equations for the surface roughness (Sa), maximal white layer thickness (WL), and MRR, this was described by the following polynomial function:
Sa = 0.38 + 0.54 I − 0.027 I2 − 0.0004 ton + 0.00002 ton2 + 0.00006 I ton
− 0.000007 I ton2 + 0.0003 I2 ton
WL = 1.422 + 1.697 I − 0.075 I2 − 0.101 ton + 0.0003 ton2 − 0.002 toff + 0,035 I ton − 0.0001 I ton2
− 0.0014 I2 ton + 0.000005 I2 ton2 + 0.00027 ton toff − 0.000001 ton2 toff
MRR = −1.2087 + 0.342 I + 0.02967 I2 − 0.00817 ton + 0.00004 ton2 + 0.02287 toff
+ 0.00096 I ton − 0.00001 I ton2 + 0.00057 I2 ton − 0.00668 I toff
Analyses of the results of the MRR showed that the values of the R-squared for surface roughness (Sa), the maximal thickness of the white layer, and the MRR were over 98%, 99%, and 99%, respectively. This result indicated that the regression models provided an excellent explanation of the relationship between the independent variables and the response Sa, WL, and MRR. Differences between the R-squared and the R-adjustable were smaller than 0.2, which indicated that the established model was adequate in representing the process. The developed models can be used to predict the values of surface roughness (Sa), the maximal value of white layer thickness (WL), and the material removal rate (MRR). Comparisons between the results of experimental studies and the values calculated based on the developed models for Sa, WL, and MRR are shown in Figure 10. The results indicated that the predicted values were very close to the experimental data.
Developed models for the Sa, WL, and MRR were checked using additional statistical tests, which confirmed the basic assumptions using the ANOVA. Residuals had a normal distribution, constant variance, and were independent of an order of data. The assumption of constant variance was checked by plotting the residuals versus the predicted values. The normality assumption and independence of residuals were checked by plotting the expected normal value versus the residuals and the residuals versus the order of data, respectively. Analysis of the residual normal probability plots (Figure 11a, Figure 12a and Figure 13a) showed that the residuals had normal distributions. Plots of the residuals versus the predicted values (Figure 11b, Figure 12b and Figure 13b); and the residuals versus the case number values (Figure 11c, Figure 12c and Figure 13c) showed that the residuals had a stochastic nature. The analysis of the plotted residuals versus the case values indicated that the error terms were independent of one another. The analyses of the residuals confirmed that the developed models were adequate.
To better understand the influence of the EDM parameters on surface roughness (Sa), the maximal thickness of white layer, and the material removal rate, the response surface plots were estimated. Based on the regression models (Equations (10)–(12)), the influence of the discharge current I, pulse time ton, and time interval toff on the Sa, WL, and MRR is shown in Figure 14, Figure 15 and Figure 16, respectively.
The results of the experimental studies indicated that the main parameters that influenced surface roughness (Sa), were the discharge current and the pulse time (Figure 14). Time interval, in the case of stability discharges, does not have a significant impact on the surface texture properties. The surface roughness (Sa) increases with the growth of the discharge current and the pulse time. These two parameters, with constant voltage, determine the amount of energy of the electrical discharge. At the lowest value of the discharge current, the changing of the pulse time does not generate a crater with greater depth. The surface roughness (Sa) does not change significantly. With the increase of the discharge current, the amount of energy delivered to the workpiece causes the melting and evaporation of a higher volume of material, which generates a crater with a larger depth. The presented dependence also has effects on the material removal rate. The MRR, similar to surface roughness, is influenced by the volume of material which is removed in single discharges, and it mainly depends on the discharge current (Figure 16). The time interval is responsible for the stabilization of the conditions in the gap after discharges. In a stable EDM process, increases of the time interval result in the decrease of the material removal rate. Discharge energy devoted to the heat flux, the mechanism of growing the plasma channel, and the removal process of the material implies the thickness of the white layer. Figure 15 shows the estimated response surface plot for the white layer thickness in relation to the EDM parameters. The increase in the pulse time and current resulted in an increase in the amount of melted and evaporated material from the discharge zone. However, more material which was melted in the single crater was not removed from the surface of the workpiece and it re-solidified on the core.
In the industrial application of the developed models, splitting the electrical discharge machining into several steps should be considered. In the first stage, the material will be removed from the workpiece using the highest discharge energy (roughing technology). In the roughing step, the applied parameters should ensure the maximum removal rate. In the next step, the parameters of the process should be changed to achieve proper surface layer properties and a low roughness for the machined surface of the manufactured parts. The EDM process will be conducted with the semi-finishing and finishing step, with respectively lower discharge energies. In the last step—the finishing treatment—it is vital to obtain the appropriate surface roughness and thickness of the white layer. However, the finishing can take more time than the roughing. The result is that a significant increase in the cost of production may be observed. In the finishing machining, a combination of minimum surface roughness with the minimum value of white layer thickness, and with a possibly maximum MRR is desirable. In the case of EDM, the simultaneous achievement of these three goals is conflicting. For that reason, in considering the properties of the EDM optimization, it should be based on the desirability technique. This method uses the Derringer's [57] desirability function, which in the case of the same importance of each response can be described by the equation:
D = ( d 1 × d 2 × ..... × d n ) 1 / n = ( i = 1 n d i ) 1 n
where n is the number of responses in the measure.
The desired function is established for each investigated response, di( y ^ i ), and it has a range from zero to one (one being the most desirable). Different desirable functions can be built, depending on the adopted optimization criteria which determine the desirable value, maximal (upper-Ui) or minimal (lower-Li). If the response for the investigated parameter should be minimized, then di( y ^ i ) can be calculated according to the following equation:
d i ( y ^ i ) = { 1 y ^ i < L i ( U i y ^ i U i L i ) t , L i y ^ i U i 0 y ^ i > U i
If the desirable function should be maximized, then it can be expressed by the following equation:
d i ( y ^ i ) = { 0 y ^ i < L i ( y ^ i L i U i L i ) s , L i y ^ i U i 1 y ^ i > U i
Calculations of the desirable function consider the extent to which the estimated values ( y ^ i ) are close to the minimum or maximum. Figure 17 presents a graphical interpretation of the desirability functions with the “importance” levels s and t. When considering the case when the “importance” t (for minimum) and s (for maximum) is large, the desirability is low unless the response moves close to the target. For low-value parameters, t and s desirability has a high value for a wide range of responses. It means that it is possible to achieve satisfactory desirability not only in the target value (minimum Li or maximum Ui).
The multi-response optimization was divided into three cases. In the first case, we examined the optimal parameters of the finishing EDM. For this task, the goal was defined as minimizing the surface roughness and white layer thickness, whilst maximizing the maximum material removal rate. The second case was to find the optimal parameters for semi-finishing EDM. In this task, we aimed to achieve an average value of the MRR (about 14.5 mm3/min), with the possibility of minimizing the surface roughness and white layer thickness. In the last case, the optimal parameters for roughing EDM were considered. In this task, the goal of optimization was to achieve the maximum MRR possible, with the possibility of minimizing the surface roughness and white layer thickness. These three cases of the optimization of electrical discharge machining were carried out using the desirability function, based on the regression Equations (10)–(12). It should be maintained that the success of the optimization with the desirability function mainly depends on the quality of the regressions models. In this study, each established model had a coefficient of determination, R-squared, that was over 98%, and the differences between the R-squared and the R-adjustable were smaller than 0.2, which indicated that the models were adequate in representing the process. For each EDM parameter (discharge current, pulse time, time interval) there was a simultaneous analysis of every combination of the factors for each of the nine responses (Figure 14, Figure 15 and Figure 16). A multi-response optimization procedure was performed for the global desirability function. The ranges for the constraints and factors for optimization are shown in Table 9.
The results of the multi-response optimization procedure of the global desirability function for the finishing, semi-finishing, and roughing operations are shown in the contour plots in Figure 18, Figure 19, Figure 20 and Figure 21, and in Table 10. The desirable function in the first case (finishing EDM) reached 0.95 (Figure 19). The semi-finishing and roughing operations reached 0.98 and 0.99, respectively. The optimal EDM parameters for finishing electrical discharge machining were a discharge current I = 3 A, pulse time ton = 176 µs, and pulse interval toff = 10 µs. The predicted surface roughness (1.7 µm) and the white layer thickness (6 µm) after optimization were close to the results obtained in the experimental studies for sample number four. Nevertheless, the material removal rate grew almost seven times and reached an MRR = 1.1 mm3/min. The increase of the MRR was achieved, along with the minimization of the surface roughness and the white layer thickness, which has a significant effect on productivity.
Figure 20 and Table 10 present the results of the optimization for semi-finishing EDM. Values of the optimal EDM parameters were as follows: discharge current I = 14 A, pulse time ton = 52 µs, and pulse interval toff = 24 µs. In this case, if the material removal rate reached 14.5 mm3/min (i.e., the average value of the MRR from experimental studies), the optimized surface roughness and the white layer thickness were as follows: Sa = 5.2 µm and WL = 15 µm. Figure 21 and Table 10 present the results of the optimization for roughing EDM. Values of the optimal EDM parameters as follows: I = 14 A, ton = 361 µs, and toff = 24 µs. The optimized material removal rate, surface roughness, and white layer thickness were as follows: MRR = 29.2 mm3/min, Sa = 12.1 µm, WL = 28.8 µm.
In the last stage of the experimental investigations for the optimal EDM parameters, a confirmation test was conducted on the EDM machine Charmilles Form 2LC ZNC (Bern, Switzerland). Table 10 presents the results of the validations of multi-response optimizations. The maximal errors between the predicted and the obtained values were 6%, which could be considered a very good result. The calculated measurement uncertainty (at 95% confidence) of the Sa and MRR is much lower than the measurement uncertainty of the thickness of the white layer. Moreover, prediction errors are in the range of maximal measurement uncertainty.

4. Conclusions

In this study, the EDM of tool steel 55NiCrMoV7 was analyzed and optimized for three cases: roughing, semi-finishing, and finishing machining. Based on the theoretical analyses and experimental research, the following conclusions were obtained:
  • Experimental research on the influence of discharge current, pulse time, and pulse interval on the surface roughness (Sa), white layer thickness, and the MRR showed that the discharge current had the main effect on Sa, WL, and the MRR. With an increase in the discharge current and pulse time, the amount of energy delivered to the workpiece caused the melting and evaporation of a higher volume of material, which generated craters with a larger depth and diameter. However, more material which melted in the single crater was not removed from the surface of the workpiece and it re-solidified on the core. The time interval between pulses did not significantly affect the change in surface integrity and the MRR, but it played an important role in the stability of the process.
  • The desirability function was used in the multi-response optimization of three functions: Sa, WL, and MRR. For the three cases of EDM—finishing, semi-finishing, and roughing operations—the optimal parameters were established. The confirmation tests for the established optimal parameters showed that the maximal errors between the predicted and the obtained values did not exceed 6%, which could be considered as a very good result.
  • The developed regression equations could be used in electrical discharge machining as a guideline for the selection of EDM parameters.

Author Contributions

Conceptualization, R.Ś.; Methodology, R.Ś., D.O.-Ś. and T.C.; Software, D.O.-Ś.; Validation, R.Ś., D.O.-Ś. and T.C.; Formal Analysis, R.Ś. and T.C.; Investigation, R.Ś. and D.O.-Ś.; Resources, R.Ś.; Data Curation, R.Ś. and D.O.-Ś.; Writing-Original Draft Preparation, R.Ś.; Writing-Review & Editing, RŚ, D.O.-Ś., and T.C.; Visualization, D.O.-Ś.; Supervision, R.Ś.; Project Administration, R.Ś.; Funding Acquisition, R.Ś.

Funding

This research was funded by using the statutory subsidy of the Faculty of Production Engineering of the Warsaw University of Technology in 2018.

Acknowledgments

Thanks are given to the assistance from the technical specialist Eng. Adrian Kopytowski and Eng. Rafał Nowicki.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hsu, W.-H.; Chien, W.-T. Effect of electrical discharge machining on stress concentration in titanium alloy holes. Materials 2016, 9, 957. [Google Scholar] [CrossRef] [PubMed]
  2. Ramulu, M.; Spaulding, M. Drilling of hybrid titanium composite laminate (HTCL) with electrical discharge machining. Materials 2016, 9, 746. [Google Scholar] [CrossRef] [PubMed]
  3. Górka, J. Assessment of steel subjected to the thermomechanical control process with respect to weldability. Metals 2018, 8, 169. [Google Scholar] [CrossRef]
  4. Świercz, R.; Oniszczuk-Świercz, D.; Dąbrowski, L.; Zawora, J. Optimization of machining parameters of electrical discharge machining tool steel 1.2713. AIP Conf. Proc. 2018, 2017, 020032. [Google Scholar] [CrossRef]
  5. Liang, J.F.; Liao, Y.S.; Kao, J.Y.; Huang, C.H.; Hsu, C.Y. Study of the EDM performance to produce a stable process and surface modification. Int. J. Adv. Manuf. Technol. 2018, 95, 1743–1750. [Google Scholar] [CrossRef]
  6. Valíček, J.; Držík, M.; Hryniewicz, T.; Harničárová, M.; Rokosz, K.; Kušnerová, M.; Barčová, K.; Bražina, D. Non-contact method for surface roughness measurement after machining. Meas. Sci. Rev. 2012, 12, 184–188. [Google Scholar] [CrossRef]
  7. Klocke, F.; Zeis, M.; Klink, A. Interdisciplinary modelling of the electrochemical machining process for engine blades. CIRP Ann. 2015, 64, 217–220. [Google Scholar] [CrossRef]
  8. Kozak, J.; Zybura-Skrabalak, M. Some problems of surface roughness in electrochemical machining (ECM). Procedia CIRP 2016, 42, 101–106. [Google Scholar] [CrossRef]
  9. Gusarov, A.V.; Grigoriev, S.N.; Volosova, M.A.; Melnik, Y.A.; Laskin, A.; Kotoban, D.V.; Okunkova, A.A. On productivity of laser additive manufacturing. J. Mater. Process. Technol. 2018, 261, 213–232. [Google Scholar] [CrossRef]
  10. Wyszyński, D.; Ostrowski, R.; Zwolak, M.; Bryk, W. Laser beam machining of polycrystalline diamond for cutting tool manufacturing. AIP Conf. Proc. 2017, 1896, 180007. [Google Scholar] [CrossRef]
  11. Chmielewski, T.; Golański, D.; Włosiński, W.; Zimmerman, J. Utilizing the energy of kinetic friction for the metallization of ceramics. Bull. Pol. Acad. Sci. Tech. Sci. 2015, 63, 201–207. [Google Scholar] [CrossRef] [Green Version]
  12. Rokosz, K.; Hryniewicz, T.; Matýsek, D.; Raaen, S.; Valíček, J.; Dudek, Ł.; Harničárová, M. SEM, EDS and XPS analysis of the coatings obtained on titanium after plasma electrolytic oxidation in electrolytes containing copper nitrate. Materials 2016, 9, 318. [Google Scholar] [CrossRef] [PubMed]
  13. Skoczypiec, S.; Ruszaj, A. A sequential electrochemical–electrodischarge process for micropart manufacturing. Precis. Eng. 2014, 38, 680–690. [Google Scholar] [CrossRef]
  14. Gołąbczak, M.; Święcik, R.; Gołąbczak, A.; Nouveau, C.; Jacquet, P.; Blanc, C. Investigations of surface layer temperature and morphology of hard machinable materials used in aircraft industry during abrasive electrodischarge grinding process. Materialwissenschaft und Werkstofftechnik 2018, 49, 568–576. [Google Scholar] [CrossRef]
  15. Kelemesh, A.; Gorbenko, O.; Dudnikov, A.; Dudnikov, I. Research of wear resistance of bronze bushings during plastic vibration deformation. East.-Eur. J. Enterp. Technol. 2017, 2, 16–21. [Google Scholar] [CrossRef] [Green Version]
  16. Golański, D.; Dymny, G.; Kujawińska, M.; Chmielewski, T. Experimental investigation of displacement/strain fields in metal coatings deposited on ceramic substrates by thermal spraying. Solid State Phenom. 2016, 240, 174–182. [Google Scholar] [CrossRef]
  17. Salacinski, T.; Winiarski, M.; Chmielewski, T.; Świercz, R. Surface finishing using ceramic fibre brush tools. In Proceedings of the 26th International Conference on Metallurgy and Materials, Brno, Czech Republi, 24–26 May 2017; pp. 1220–1226. [Google Scholar]
  18. Guo, J.; Wang, H.; Goh, M.H.; Liu, K. Investigation on surface integrity of rapidly solidified aluminum RSA 905 by magnetic field-assisted finishing. Micromachines 2018, 9, 146. [Google Scholar] [CrossRef]
  19. Kunieda, M.; Lauwers, B.; Rajurkar, K.P.; Schumacher, B.M. Advancing EDM through fundamental insight into the process. CIRP Ann. Manuf. Technol. 2005, 54, 64–87. [Google Scholar] [CrossRef]
  20. Izquierdo, B.; Sánchez, J.A.; Plaza, S.; Pombo, I.; Ortega, N. A numerical model of the EDM process considering the effect of multiple discharges. Int. J. Mach. Tools Manuf. 2009, 49, 220–229. [Google Scholar] [CrossRef]
  21. Izquierdo, B.; Sánchez, J.A.; Ortega, N.; Plaza, S.; Pombo, I. Insight into fundamental aspects of the EDM process using multidischarge numerical simulation. Int. J. Adv. Manuf. Technol. 2011, 52, 195–206. [Google Scholar] [CrossRef]
  22. Ming, W.; Zhang, Z.; Wang, S.; Huang, H.; Zhang, Y.; Zhang, Y.; Shen, D. Investigating the energy distribution of workpiece and optimizing process parameters during the EDM of Al6061, Inconel718, and SKD11. Int. J. Adv. Manuf. Technol. 2017, 92, 4039–4056. [Google Scholar] [CrossRef]
  23. Gulbinowicz, Z.; Świercz, R.; Oniszczuk-Świercz, D. Influence of electrical parameters in electro discharge machining of tungsten heavy alloys on surface texture properties. AIP Conf. Proc. 2018, 2017, 020007. [Google Scholar] [CrossRef]
  24. Salcedo, A.T.; Arbizu, I.P.; Pérez, C.J.L. Analytical modelling of energy density and optimization of the EDM machining parameters of Inconel 600. Metals 2017, 7, 166. [Google Scholar] [CrossRef]
  25. Rahang, M.; Patowari, P.K. Parametric optimization for selective surface modification in EDM using taguchi analysis. Mater. Manuf. Process. 2016, 31, 422–431. [Google Scholar] [CrossRef]
  26. Vagaská, A.; Gombár, M. Comparison of usage of different neural structures to predict AAO layer thickness. Tehnički Vjesnik 2017, 24, 333–339. [Google Scholar] [CrossRef]
  27. Wojciechowski, S.; Maruda, R.W.; Królczyk, G.M. The application of response surface method to optimization of precision ball end milling. MATEC Web Conf. 2017, 112, 01004. [Google Scholar] [CrossRef] [Green Version]
  28. Adalarasan, R.; Santhanakumar, M.; Rajmohan, M. Optimization of laser cutting parameters for Al6061/SiCp/Al2O3 composite using grey based response surface methodology (GRSM). Measurement 2015, 73, 596–606. [Google Scholar] [CrossRef]
  29. Vera Candioti, L.; De Zan, M.M.; Cámara, M.S.; Goicoechea, H.C. Experimental design and multiple response optimization. Using the desirability function in analytical methods development. Talanta 2014, 124, 123–138. [Google Scholar] [CrossRef]
  30. Ayesta, I.; Izquierdo, B.; Sanchez, J.A.; Ramos, J.M.; Plaza, S.; Pombo, I.; Ortega, N. Optimum electrode path generation for EDM manufacturing of aerospace components. Robot. Comput.-Integr. Manuf. 2016, 37, 273–281. [Google Scholar] [CrossRef]
  31. Abidi, M.H.; Al-Ahmari, A.M.; Siddiquee, A.N.; Mian, S.H.; Mohammed, M.K.; Rasheed, M.S. An investigation of the micro-electrical discharge machining of nickel-titanium shape memory alloy using grey relations coupled with principal component analysis. Metals 2017, 7, 486. [Google Scholar] [CrossRef]
  32. Rubio, E.M.; Villeta, M.; Valencia, J.L.; Sáenz de Pipaón, J.M. Experimental study for improving the repair of magnesium–aluminium hybrid parts by turning processes. Metals 2018, 8, 59. [Google Scholar] [CrossRef]
  33. 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]
  34. Kilickap, E.; Yardimeden, A.; Çelik, Y.H. Mathematical modelling and optimization of cutting force, tool wear and surface roughness by using artificial neural network and response surface methodology in milling of Ti-6242S. Appl. Sci. 2017, 7, 1064. [Google Scholar] [CrossRef]
  35. Rogalski, G.; Fydrych, D.; Łabanowski, J. Underwater wet repair welding of API 5L X65M pipeline steel. Pol. Marit. Res. 2017, 24, 188–194. [Google Scholar] [CrossRef]
  36. Unune, D.R.; Mali, H.S. Parametric modeling and optimization for abrasive mixed surface electro discharge diamond grinding of Inconel 718 using response surface methodology. Int. J. Adv. Manuf. Technol. 2017, 93, 3859–3872. [Google Scholar] [CrossRef]
  37. Hlaváč, L.M.; Krajcarz, D.; Hlaváčová, I.M.; Spadło, S. Precision comparison of analytical and statistical-regression models for AWJ cutting. Precis. Eng. 2017, 50, 148–159. [Google Scholar] [CrossRef]
  38. Ghodsiyeh, D.; Golshan, A.; Izman, S. Multi-objective process optimization of wire electrical discharge machining based on response surface methodology. J. Braz. Soc. Mech. Sci. Eng. 2014, 36, 301–313. [Google Scholar] [CrossRef]
  39. Alavi, F.; Jahan, M.P. Optimization of process parameters in micro-EDM of Ti-6Al-4V based on full factorial design. Int. J. Adv. Manuf. Technol. 2017, 92, 167–187. [Google Scholar] [CrossRef]
  40. Selvarajan, L.; Manohar, M.; Kumar, A.U.; Dhinakaran, P. Modelling and experimental investigation of process parameters in EDM of Si3N4-TiN composites using GRA-RSM. J. Mech. Sci. Technol. 2017, 31, 111–122. [Google Scholar] [CrossRef]
  41. Świercz, R.; Oniszczuk-Świercz, D. Influence of electrical discharge pulse energy on the surface integrity of tool steel 1.2713. Proceedings of 26th International Conference on Metallurgy and Materials, Brno, Czech Republic, 24–26 May 2017; pp. 1452–1457. [Google Scholar]
  42. Kumaran, S.T.; Ko, T.J.; Kurniawan, R. Grey fuzzy optimization of ultrasonic-assisted EDM process parameters for deburring CFRP composites. Measurement 2018, 123, 203–212. [Google Scholar] [CrossRef]
  43. Gu, L.; Zhu, Y.; Zhang, F.; Farhadi, A.; Zhao, W. Mechanism analysation and parameter optimisation of electro discharge machining of titanium-zirconium-molybdenum alloy. J. Manuf. Process. 2018, 32, 773–781. [Google Scholar] [CrossRef]
  44. Dang, X.-P. Constrained multi-objective optimization of EDM process parameters using kriging model and particle swarm algorithm. Mater. Manuf. Process. 2018, 33, 397–404. [Google Scholar] [CrossRef]
  45. Mohanty, C.P.; Mahapatra, S.S.; Singh, M.R. An intelligent approach to optimize the EDM process parameters using utility concept and QPSO algorithm. Eng. Sci. Technol. Int. J. 2017, 20, 552–562. [Google Scholar] [CrossRef] [Green Version]
  46. Maity, K.; Mishra, H. ANN modelling and Elitist teaching learning approach for multi-objective optimization of μ-EDM. J. Intell. Manuf. 2018, 29, 1599–1616. [Google Scholar] [CrossRef]
  47. Tripathy, S.; Tripathy, D.K. Multi-response optimization of machining process parameters for powder mixed electro-discharge machining of H-11 die steel using grey relational analysis and topsis. Mach. Sci. Technol. 2017, 21, 362–384. [Google Scholar] [CrossRef]
  48. Nguyen, H.-P.; Pham, V.-D.; Ngo, N.-V. Application of TOPSIS to Taguchi method for multi-characteristic optimization of electrical discharge machining with titanium powder mixed into dielectric fluid. Int. J. Adv. Manuf. Technol. 2018, 98, 1179–1198. [Google Scholar] [CrossRef]
  49. Roy, T.; Dutta, R.K. Integrated fuzzy AHP and fuzzy TOPSIS methods for multi-objective optimization of electro discharge machining process. Soft Comput. 2018. [Google Scholar] [CrossRef]
  50. Kandpal, B.C.; Kumar, J.; Singh, H. Optimization and characterization of EDM of AA 6061/10%Al2O3 AMMC using Taguchi’s approach and utility concept. Prod. Manuf. Res. 2017, 5, 351–370. [Google Scholar] [CrossRef]
  51. D’Urso, G.; Giardini, C.; Quarto, M.; Maccarini, G. Cost index model for the process performance optimization of micro-EDM drilling on tungsten carbide. Micromachines 2017, 8, 251. [Google Scholar] [CrossRef]
  52. Parsana, S.; Radadia, N.; Sheth, M.; Sheth, N.; Savsani, V.; Prasad, N.E.; Ramprabhu, T. Machining parameter optimization for EDM machining of Mg–RE–Zn–Zr alloy using multi-objective Passing Vehicle Search algorithm. Arch. Civil Mech. Eng. 2018, 18, 799–817. [Google Scholar] [CrossRef]
  53. Hadad, M.; Bui, L.Q.; Nguyen, C.T. Experimental investigation of the effects of tool initial surface roughness on the electrical discharge machining (EDM) performance. Int. J. Adv. Manuf. Technol. 2018, 95, 2093–2104. [Google Scholar] [CrossRef]
  54. Baruffi, F.; Parenti, P.; Cacciatore, F.; Annoni, M.; Tosello, G. On the application of replica molding technology for the indirect measurement of surface and geometry of micromilled components. Micromachines 2017, 8, 195. [Google Scholar] [CrossRef]
  55. ISO. IOS 15530-3: Geometrical Product Specifications (GPS)—Coordinate Measuring Machines (CMM): Technique for Determining the Uncertainty of Measurement; IOS: Geneva, Switzerland, 2011. [Google Scholar]
  56. Xin, B.; Li, S.; Yin, X.; Lu, X. Dynamic observer modeling and minimum-variance self-tuning control of EDM interelectrode gap. Appl. Sci. 2018, 8, 1443. [Google Scholar] [CrossRef]
  57. 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]
Figure 1. The schematic diagram of the experimental set-up.
Figure 1. The schematic diagram of the experimental set-up.
Micromachines 10 00072 g001
Figure 2. The recorded current and voltage waveforms.
Figure 2. The recorded current and voltage waveforms.
Micromachines 10 00072 g002
Figure 3. The recorded current and voltage waveforms.
Figure 3. The recorded current and voltage waveforms.
Micromachines 10 00072 g003
Figure 4. The recorded voltage and current waveforms for the following parameters: (a) Uc = 25 V, Uz = 230 V, I = 3 A, ton = 400 µs, toff = 100 µs; (b) Uc = 25 V, Uz = 230 V, I = 3 A, ton = 13 μs, toff = 13 μs.
Figure 4. The recorded voltage and current waveforms for the following parameters: (a) Uc = 25 V, Uz = 230 V, I = 3 A, ton = 400 µs, toff = 100 µs; (b) Uc = 25 V, Uz = 230 V, I = 3 A, ton = 13 μs, toff = 13 μs.
Micromachines 10 00072 g004
Figure 5. The surface texture of the tool steel 55NiCrMoV7 after electrical discharge machining (EDM): (a) Uc = 25 V, I = 8.5 A, ton = 400 µs, toff = 150 µs; (b) Uc = 25 V, I = 3 A, ton = 400 µs, toff = 80 µs.
Figure 5. The surface texture of the tool steel 55NiCrMoV7 after electrical discharge machining (EDM): (a) Uc = 25 V, I = 8.5 A, ton = 400 µs, toff = 150 µs; (b) Uc = 25 V, I = 3 A, ton = 400 µs, toff = 80 µs.
Micromachines 10 00072 g005
Figure 6. The metallographic structure of tool steel 55NiCrMoV7 after (EDM).
Figure 6. The metallographic structure of tool steel 55NiCrMoV7 after (EDM).
Micromachines 10 00072 g006
Figure 7. The metallographic structure of tool steel 55NiCrMoV7 after EMD: Uc = 25 V, I = 14 A, ton = 13 µs, toff = 10 µs.
Figure 7. The metallographic structure of tool steel 55NiCrMoV7 after EMD: Uc = 25 V, I = 14 A, ton = 13 µs, toff = 10 µs.
Micromachines 10 00072 g007
Figure 8. The metallographic structure of tool steel 55NiCrMoV7 after EMD: (a) Uc = 25 V, I = 3 A, ton = 206 µs, toff = 80 µs; (b) Uc = 25 V, I = 14 A, ton = 400 µs, toff = 80 µs.
Figure 8. The metallographic structure of tool steel 55NiCrMoV7 after EMD: (a) Uc = 25 V, I = 3 A, ton = 206 µs, toff = 80 µs; (b) Uc = 25 V, I = 14 A, ton = 400 µs, toff = 80 µs.
Micromachines 10 00072 g008
Figure 9. The structure of the surface after EMD: Uc = 25 V, I = 14 A, ton = 400 µs, toff = 150 µs.
Figure 9. The structure of the surface after EMD: Uc = 25 V, I = 14 A, ton = 400 µs, toff = 150 µs.
Micromachines 10 00072 g009
Figure 10. The comparison between the results of experimental studies to the values calculated based on the developed models for (a) surface roughness (Sa); (b) maximal thickness of white layer (WL); (c) material removal rate (MRR).
Figure 10. The comparison between the results of experimental studies to the values calculated based on the developed models for (a) surface roughness (Sa); (b) maximal thickness of white layer (WL); (c) material removal rate (MRR).
Micromachines 10 00072 g010
Figure 11. The plots to check the model for surface roughness (Sa): (a) the normal plot of residuals; (b) the residuals versus the predicted values; and (c) the residuals versus the case values.
Figure 11. The plots to check the model for surface roughness (Sa): (a) the normal plot of residuals; (b) the residuals versus the predicted values; and (c) the residuals versus the case values.
Micromachines 10 00072 g011
Figure 12. The plots to check the model for maximal white layer thickness: (a) the normal plot of residuals; (b) the residuals versus the predicted values; and (c) the residuals versus the case values.
Figure 12. The plots to check the model for maximal white layer thickness: (a) the normal plot of residuals; (b) the residuals versus the predicted values; and (c) the residuals versus the case values.
Micromachines 10 00072 g012
Figure 13. The plots to check the model MRR: (a) the normal plot of residuals; (b) the residuals versus the predicted values; and (c) the residuals versus the case values.
Figure 13. The plots to check the model MRR: (a) the normal plot of residuals; (b) the residuals versus the predicted values; and (c) the residuals versus the case values.
Micromachines 10 00072 g013
Figure 14. The estimated response surface plot for roughness (Sa): (a) constant toff = 80 µs; (b) constant ton = 206 µs; and (c) constant I = 8.5 A.
Figure 14. The estimated response surface plot for roughness (Sa): (a) constant toff = 80 µs; (b) constant ton = 206 µs; and (c) constant I = 8.5 A.
Micromachines 10 00072 g014
Figure 15. The estimated response surface plot for the WL: (a) constant toff = 80 µs; (b) constant ton = 206 µs; and (c) constant I = 8.5 A.
Figure 15. The estimated response surface plot for the WL: (a) constant toff = 80 µs; (b) constant ton = 206 µs; and (c) constant I = 8.5 A.
Micromachines 10 00072 g015
Figure 16. The estimated response surface plot for the MRR: (a) constant toff = 80 µs; (b) constant ton = 206 µs; and(c) constant I = 8.5 A.
Figure 16. The estimated response surface plot for the MRR: (a) constant toff = 80 µs; (b) constant ton = 206 µs; and(c) constant I = 8.5 A.
Micromachines 10 00072 g016
Figure 17. The desirability functions for the target value: (a) maximized; (b) minimized criteria.
Figure 17. The desirability functions for the target value: (a) maximized; (b) minimized criteria.
Micromachines 10 00072 g017
Figure 18. The contour plots of the desirability function for the finishing EDM optimization.
Figure 18. The contour plots of the desirability function for the finishing EDM optimization.
Micromachines 10 00072 g018
Figure 19. The example of profiles for the predicted values and desirability for the finishing EDM.
Figure 19. The example of profiles for the predicted values and desirability for the finishing EDM.
Micromachines 10 00072 g019
Figure 20. The contour plots of the desirability function for semi-finishing EDM optimization.
Figure 20. The contour plots of the desirability function for semi-finishing EDM optimization.
Micromachines 10 00072 g020
Figure 21. The contour plots of the desirability function for roughing EDM optimization.
Figure 21. The contour plots of the desirability function for roughing EDM optimization.
Micromachines 10 00072 g021
Table 1. The uncertainty contributions for the Sa roughness measurements on the scanning profilometer.
Table 1. The uncertainty contributions for the Sa roughness measurements on the scanning profilometer.
Uncertainty Contributions (nm)
ucalupures,PROFUPROFUSa, EDMU95,Sa
202320.5542
Table 2. The uncertainty contributions for the thickness of the white layer measurements on microscope.
Table 2. The uncertainty contributions for the thickness of the white layer measurements on microscope.
Uncertainty Contributions (μm)
ucalupures,OMUOMU95,WL
0.0600.0480.3120.3210.6
Table 3. The uncertainty contributions for the weight measurements on the precision balance.
Table 3. The uncertainty contributions for the weight measurements on the precision balance.
Uncertainty Contributions (mg)
um1uresuiuieUBU95,W
0.020.00290.00580.010.0230.046
Table 4. The levels of machining parameters carried out in the experimental design.
Table 4. The levels of machining parameters carried out in the experimental design.
EDM ParametersLevel 1Level 2Level 3
discharge current I (A)38.514
pulse time ton (μs)13206400
time interval toff (μs)980150
Table 5. The design of the experimental matrix.
Table 5. The design of the experimental matrix.
Exp.
no.
EDM ParametersObserved Values
Discharge Current I (A)Pulse Time ton (μs)Time Interval
toff (μs)
Surface Roughness Sa (μm)Maximal Thickness of the White Layer (μm)MRR
(mm3/min)
1313102.05.50.54
28.513103.111.53.47
31413103.81211.06
4313801.960.17
58.513803.0121.18
61413803.411.53.21
73131501.960.10
88.5131503.011.50.55
914131503.3121.31
103206101.970.51
118.5206106.2228.09
1214206109.325.428.46
133206801.9100.36
148.5206806.0245.77
15142068010.52819.23
1632061501.8100.29
178.52061505.4254.68
181420615011.73215.48
193400102.4120.37
208.5400103.9176.58
21144001012.32829.19
223400802.413.50.34
238.5400804.0205.61
24144008012.72924.84
2534001502.5140.28
268.54001504.918.42.56
271440015011.533.520.31
288.5206806.124.55.88
Table 6. The analysis of variance (ANOVA) table for the Sa (after elimination).
Table 6. The analysis of variance (ANOVA) table for the Sa (after elimination).
SourceSum of SquaresDegrees of FreedomMean SquareF-ValueProb > fContribution %
Model344.7600749.027150.95<0.0001-
I198.54681198.5468611.29<0.000157.6
I25.809715.809717.880.00041.7
ton54.1840154.1840166.82<0.000115.7
ton216.1085116.108549.59<0.00014.7
I ton50.0208150.0208154.01<0.000114.5
I ton28.923518.923527.47<0.00012.6
I2 ton11.1696111.169634.39<0.00013.2
Error6.4953200.32479---
Total SS351.256027R-sqr = 0.98R-Adj = 0.97
Table 7. The ANOVA table for the maximal thickness of the white layer (after elimination).
Table 7. The ANOVA table for the maximal thickness of the white layer (after elimination).
SourceSum of SquaresDegrees of FreedomMean SquareF-ValueProb > fContribution %
Model1886.36611171.48141.26<0.0001-
I896.6561896.656738.590.002247.5
I216.041116.04113.21<0.00010.8
ton524.8801524.880432.35<0.000127.8
ton2174.3661174.366143.620.00029.2
toff27.406127.40622.57<0.00011.4
I ton90.750190.75074.750.00044.8
I ton261.584161.58350.720.00313.3
I2ton37.210137.21030.650.00712.0
I2ton244.018144.01736.25<0.00012.3
tontoff5.60315.6034.61<0.00010.3
ton2toff7.86017.8596.470.04730.4
Error19.424161.2140738.590.0216-
Total SS1905.79927R-sqr = 0.99 R-Adj = 0.98
Table 8. The ANOVA table for the material removal rate (MRR) (after elimination).
Table 8. The ANOVA table for the material removal rate (MRR) (after elimination).
SourceSum of SquaresDegrees of FreedomMean SquareF-ValueProb > fContribution %
Model2243.499247.881208.65<0.0001-
I1253.28711253.2871055.08<0.000155.9
I2126.2431126.243106.27<0.00015.6
ton260.6551260.655219.43<0.000111.6
ton256.489156.48947.55<0.00012.5
toff101.3811101.38185.34<0.00014.5
I ton285.9481285.948240.72<0.000112.7
I ton236.030136.03030.33<0.00011.6
I2ton44.106144.10637.13<0.00012.0
I toff79.350179.35066.80<0.00013.5
Error21.381181.188---
Total SS2264.8727R-sqr = 0.99R-Adj = 0.99
Table 9. The goals and factor range for optimization.
Table 9. The goals and factor range for optimization.
FactorsGoalLower LimitUpper LimitWeightImportance
Finishing EDMSemi-FinishingRoughing
I (A)In range3141---
ton (µs)In range134001---
toff (µs)In range101501---
Sa (µm)Minimize1.8512.71t = 5t = 3t = 0.3
WL (µm)Minimize5.533.51t = 5t = 3t = 0.3
MRR (mm3/min)Maximize0.0129.191s = 0.3s = 3s = 5
Table 10. The experimental validations of the multi-response optimizations.
Table 10. The experimental validations of the multi-response optimizations.
Optimal EDM ParametersSummary of Values Obtained in Optimization
ResponsePredictedExperimental VerificationError%
Finishing I = 3 A
ton = 176 µs
toff = 10 µs
Sa (µm)1.71.86
WL (µm)66.35
MRR (mm3/min)1.131.066
Semi- finishing I = 14 A
ton = 52 µs
toff = 24 µs
Sa (µm)5.2 5.44
WL (µm)1515.85
MRR (mm3/min)14.5153
RoughingI = 14 A
ton = 361 µs
toff = 24 µs
Sa (µm)12.1 12.75
WL (µm)28.830.56
MRR (mm3/min)29.228.14

Share and Cite

MDPI and ACS Style

Świercz, R.; Oniszczuk-Świercz, D.; Chmielewski, T. Multi-Response Optimization of Electrical Discharge Machining Using the Desirability Function. Micromachines 2019, 10, 72. https://doi.org/10.3390/mi10010072

AMA Style

Świercz R, Oniszczuk-Świercz D, Chmielewski T. Multi-Response Optimization of Electrical Discharge Machining Using the Desirability Function. Micromachines. 2019; 10(1):72. https://doi.org/10.3390/mi10010072

Chicago/Turabian Style

Świercz, Rafał, Dorota Oniszczuk-Świercz, and Tomasz Chmielewski. 2019. "Multi-Response Optimization of Electrical Discharge Machining Using the Desirability Function" Micromachines 10, no. 1: 72. https://doi.org/10.3390/mi10010072

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