Optimization of EDM process parameters using Response Surface Methodology for AISI D3 Steel

Similar documents
Associate Professor, Department of Mechanical Engineering, Birla Institute of Technology, Mesra, Ranchi, India

RESPONSE SURFACE ANALYSIS OF EDMED SURFACES OF AISI D2 STEEL

Modeling and Optimization of WEDM Process Parameters on Machining of AISI D2 steel using Response Surface Methodology (RSM)

Parameter Optimization of EDM on En36 Alloy Steel For MRR and EWR Using Taguchi Method

Optimization of Machining Parameters in Wire Cut EDM of Stainless Steel 304 Using Taguchi Techniques

Modeling of Wire Electrical Discharge Machining of AISI D3 Steel using Response Surface Methodology

International Journal of Scientific & Engineering Research, Volume 5, Issue 1, January ISSN

Modeling of Wire Electrical Discharge Machining Parameters Using Titanium Alloy (Ti-6AL-4V)

An investigation of material removal rate and kerf on WEDM through grey relational analysis

Optimization of WEDM Parameters for Super Ni-718 using Neutrosophic Sets and TOPSIS Method

Optimization of Cutting Parameter of (SS302) on EDM using Taguchi Method Chintan A. Prajapati 1 Prof. Dr. Prashant Sharma 2 Prof.

Study of the effect of machining parameters on material removal rate and electrode wear during Electric Discharge Machining of mild steel

ELECTRIC DISCHARGE MACHINING AND MATHEMATICAL MODELING OF Al-ALLOY-20 % SiC p COMPOSITES USING COPPER ELECTRODE

Impact of Microchannel Geometrical Parameters in W-EDM Using RSM

A Parametric Optimization of Electric Discharge Drill Machine Using Taguchi Approach

Neuro-fuzzy model and Regression model a comparison study of MRR in Electrical discharge machining of D2 tool steel

THere is a heavy demand of the advanced materials with

Effect and Optimization of EDM Process Parameters on Surface Roughness for En41 Steel

Optimization of MRR and SR by employing Taguchis and ANOVA method in EDM

Influence of Input Parameters on Characteristics of Electro Chemical Machining Process

MULTI-RESPONSE OPTIMIZATION OF EDM PERFORMANCE CHARACTERISTICS USING RESPONSE SURFACE METHODOLOGY AND DESIRABILITY FUNCTION

OPTIMIZATION OF MATERIAL REMOVAL RATE AND SURFACE ROUGHNESSIN WED-MACHINING OF TiNi SMA USING GREY RELATION ANALYSIS

Study of EDM Parameters on Mild Steel Using Brass Electrode

Experimental Study on Parametric Optimization of Titanium based Alloy (Ti-6Al-4V) in Electrochemical Machining Process

Modelling And Optimization Of Electro Discharge Machining In C45 Material

EXPERIMENTAL INVESTIGATION OF MRR, SURFACE ROUGHNESS AND OVERCUT OF AISI 304 STAINLESS STEEL IN EDM A THESIS SUBMITTED IN PARTIAL FULFILMENT

AnalysisoMRRandSRwithDifferentElectrodeforSS316onDi-SinkingEDMusingTaguchiTechnique

Optimization of Machining Parameters in ECM of Al/B4C Composites Using Taguchi Method

MODELING OF SURFACE ROUGHNESS IN WIRE ELECTRICAL DISCHARGE MACHINING USING ARTIFICIAL NEURAL NETWORKS

EXPERIMENTAL INVESTIGATIONS ON ORBITAL ELECTRO DISCHARGE MACHINING OF INCONEL 718 USING TAGUCHI TECHNIQUE

Statistical and regression analysis of Material Removal Rate for wire cut Electro Discharge Machining of SS 304L using design of experiments

Study of water assisted dry wire-cut electrical discharge machining

Post Graduate Scholar, Department of Mechanical Engineering, Jalpaiguri Govt. Engineering College, India. 2

Parametric Study and Optimization of WEDM Parameters for Titanium diboride TiB2

FEM ANALYSIS WITH EXPERIMENTAL RESULTS TO STUDY EFFECT OF EDM PARAMETERS ON MRR OF AISI 1040 STEEL

Optimization of Wire EDM Process Parameters of Al 6061/Al 2 O 3 /3%red mud MMC

The Influence of EDM Parameters in Finishing Stage on Surface Quality, MRR and EWR

Optimization of Radial Force in Turning Process Using Taguchi s Approach

Percentage of harmful discharges for surface current density monitoring in electrical discharge machining process

*Corresponding author: Received November 06, 2013; Revised November 20, 2013; Accepted November 29, 2013

Puttur, Andhra Pradesh, India , Andhra Pradesh, India ,

Journal of Chemical and Pharmaceutical Research, 2017, 9(10): Research Article

Experimental study of electrical discharge drilling of stainless steel UNS S30400

Mr. Harshit K. Dave 1, Dr. Keyur P. Desai 2, Dr. Harit K. Raval 3

Experimental Investigation of Micro-EDM Process on Brass using Taguchi Technique

A study to achieve a fine surface finish in Wire-EDM

EXPERIMENTAL INVESTIGATION OF MRR AND SURFACE ROUGHNESS OF EN-18 STEEL IN ECM

2. To compare results with a previous published work[6].

International Journal of Scientific & Engineering Research, Volume 5, Issue 12, December ISSN IJSER

A CASE STUDY ON APPLICATION OF FUZZY LOGIC IN ELECTRICAL DISCHARGE MACHINING (EDM)

Optimization of Process Parameters in CNC Drilling of EN 36

Investigation of effect of process parameters in micro hole drilling

Chapter 4 - Mathematical model

Analysis and Mathematical Modeling of Wire EDM Process Parameters for Commercial Graphite by Using Response Surface Methodology

MODELING AND OPTIMIZATION FOR DRILLING OF HIGH ASPECT RATIO BLIND MICRO HOLES IN MICRO EDM

Optimization of process parameter in electrochemical machining. Of Inconel 718 by Taguchi analysis

Multi-Objective Optimization of Electrochemical machining of EN31 steel by Grey Relational Analysis

Study on Erosion Mechanism of Magnetic-field-assisted Micro-EDM

Investigation on Electric Discharge Machine for optimization of Surface Roughness using Response Surface Methodology approach

Research Article Study of Tool Wear and Overcut in EDM Process with Rotary Tool and Magnetic Field

Drilling Mathematical Models Using the Response Surface Methodology

Taguchi-grey relational based multi response optimization of electrical process parameters in electrical discharge machining

Optimization of machining parameters of Wire-EDM based on Grey relational and statistical analyses

PREDICTIVE MODELLING AND ANALYSIS OF SURFACE ROUGHNESS IN ELECTRO-DISCHARGE MACHINING OF D2 TOOL STEEL USING REGRESSION AND NEURAL NETWORKS APPROACH

A Review Paper on Rotary Electro-Discharge Machining

APPLICATION OF GREY RELATIONAL ANALYSIS TO MACHINING PARAMETERS DETERMINATION OF WIRE ELECTRICAL DISCHARGE MACHINING

MICRO ELECTRO DISCHARGE MACHINING OF Γ -TITANIUM ALUMINIDE ALLOY

Decision Science Letters

MATHEMATICAL MODELING OF MATERIAL REMOVAL RATE OF T90Mn2W50Cr45 TOOL STEEL IN WIRE ELECTRICAL DISCHARGE MACHINE

Performance Evolution and Selection of Controllable Process variables in ECM for Al, B4C Metal Matrix Composites

Response surface Methodology and Desirability Approach to Optimize EDM Parameters

S. S.Mahapatra National Institute of Technology Department of Mechanical Engineering Rourkela. INDIA. Amar Patnaik.

CHAPTER -6 ANALYSIS AND DISCUSSION OF RESULTS

Analysis and Multi-objective Optimisation of Surface

ANALYSIS OF PARAMETRIC INFLUENCE ON DRILLING USING CAD BASED SIMULATION AND DESIGN OF EXPERIMENTS

Institutionen för matematik och matematisk statistik Umeå universitet November 7, Inlämningsuppgift 3. Mariam Shirdel

A Comparative Modeling and Multi- Objective Optimization in Wire EDM Process on H21 Tool Steel Using Intelligent Hybrid Approach

ADVANCES in NATURAL and APPLIED SCIENCES

Experimental Investigation of Machining Parameter in Electrochemical Machining

Optimization of Machining Process Parameters in Drilling of

Electrode set-up for EDM-drilling of large aspect-ratio microholes

Algorithm for Modeling Wire Cut Electrical Discharge Machine Parameters using Artificial Neural Network

All about sparks in EDM

Prediction and analysis of radial overcut in holes drilled by electrochemical machining process

Influence of Electrode Material and Process Parameters on Surface Quality and MRR in EDM of AISI H13 using ANN

CHAPTER 4 EXPERIMENTAL DESIGN. 4.1 Introduction. Experimentation plays an important role in new product design, manufacturing

Material removal characteristics of microslot (kerf) geometry

AS INGPEED Development Journal in Integrated Engineering (ISSN:requested) 2013, Vol. 1

Application of Taguchi method in optimization of control parameters of grinding process for cycle time reduction Snehil A. Umredkar 1, Yash Parikh 2

CHAPTER 6 MACHINABILITY MODELS WITH THREE INDEPENDENT VARIABLES

Design of arrow-head electrode in electropolishing of cylindrical part

Chapter 6 The 2 k Factorial Design Solutions

Drilling Microholes in Hot Tool Steel by Using Micro-Electro Discharge Machining

ECM PROCESS CHARACTERISTICS

Performance analysis of µed-milling process using various statistical techniques

MODELING, ANALYSIS, EVALUATION, SELECTION AND EXPERIMENTAL INVESTIGATION OF EDM PROCESS

Statistical and Experimental Study on the Influence of Input Parameters on the Dimensional Accuracy of Workpiece in EDM

Proceedings of the International Conference on Advances in Production and Industrial Engineering

THERMAL-ELECTRICAL MODELLING OF ELECTRICAL DISCHARGE MACHINING PROCESS

Transcription:

Optimization of EDM process parameters using Response Surface Methodology for AISI D3 Steel Mr.B.Gangadhar 1, Mr.N. Mahesh Kumar 2 1 Department of Mechanical Engineering, Sri Venkateswara College of Engineering and Technology,India 2 Department of Mechanical Engineering, Sri Venkateswara College of Engineering and Technology,India ABSTRACT The present work demonstrates the optimization process of material removal rate (MRR) of electrical discharge machining (EDM) by RSM (Response Surface Methodology). The work piece material was EN31 tool steel. The pulse on time, pulse off time, pulse current and voltage were the control parameters of EDM. RSM method was used to design the experiment using rotatable central composite design as this is the most widely used experimental design for modeling a second order response surface. The process has been successfully modeled using response surface methodology (RSM) and model adequacy checking is also carried out using Minitab software. The second-order response models have been validated with analysis of variance. Finally, an attempt has been made to estimate the optimum machining conditions to produce the best possible responses within the experimental constraints. KEYWORDS: Methodology CNC Machining, EDM, Material Removal Rate(MRR), Response Surface I.INTRODUCTION Technologically advanced industries like aeronautics, automobiles, nuclear reactors, missiles, turbines etc., requires materials like high strength temperature resistant alloys which have higher strength, corrosion resistance, toughness, and other diverse properties. With rapid development in the field of materials it has become essential to develop cutting tool materials and processes which can safely and conveniently machine such new materials for sustained productivity, high accuracy and versatility at automation. Consequently, nontraditional techniques of machining are providing effective solutions to the problem imposed by the increasing demand for high strength temperature resistant alloys, the requirement of parts with intricate and compacted shapes and materials so hard as to defy machining by conventional methods. Electrical Discharge Machining is a non-traditional machining technique, which is widely used to produce finish parts through the action of an electrical discharge of short duration and high current density between the tool and work piece. The tool and the work piece are free from the physical contact with each other. Generally, the EDM is used for machining of electrical conductive materials in the presence of a dielectric fluid. These are submersed in a dielectric liquid 337 P a g e

such as kerosene or deionized water. Its unique feature of using thermal energy to machine electrically conductive parts regardless of hardness has been its distinctive advantage. The electrical discharge machining process is widely used in the aerospace, automobile, die manufacturing and plastic mould industries to machine hard metals and its alloy [1]. The basic principle in EDM is the conversion of electrical energy into thermal energy through a series of discrete electrical discharges occurring between the electrode and work piece immersed in the dielectric fluid. The insulating effect of the dielectric is important in avoiding electrolysis of the electrodes during the EDM process. A spark is produced at the point of smallest inter-electrode gap by a high voltage, overcoming the strength dielectric breakdown strength of the small gap between the cathode and anode at a temperature in the range of 8000 to 12,000 C. Erosion of metal from both electrodes takes place there. The numerical control monitors the gap conditions (voltage and current) and synchronously controls the different axes and the pulse generator. The dielectric liquid is filtrated to remove debris particles and decomposition products [2].The prediction of optimal machining conditions for good surface finish plays a very important role in process planning. Speeding and Wang [3] have attempted to optimize the process parametric combinations by modeling the process using ANN and characterize the surface in wire electrical discharge machining (WEDM) on AISI 420 through time series techniques. Zhang et al. [4] have investigated the effects on material removal rate, surface roughness and diameter of discharge points in electro-discharge machining (EDM) on ceramics. From the experimental results, they have shown that the material removal rate, surface roughness and the diameter of discharge point all increase with increasing pulse-on time and discharge current. Tsai and Wang [5] have established a semi-empirical model of surface finish on work for various materials (three different grades of steel) in electrical discharge machining and the parameters of the model viz. peak current, pulse duration, electric polarity and properties of materials have been fitted based on the experimental data using Taguchi method. It is seen that the developed model is dependent on work and tool materials. Singh et al. [6] have developed a model for multi-response optimization of process parameters viz. metal removal rate, tool wear rate, taper, radial overcut and surface roughness on electrical discharge machining of Al-10%SiCp composites. Yih-fong and Fu-chen [7] have presented an approach for optimizing high-speed electrical discharge machining (EDM) using Taguchi methods. They have concluded that the most important factors affecting the EDM process robustness have been identified as pulse-on time, duty cycle, and pulse peak current. Fig. 1 shows Fishbone diagram showing parameters affecting MRR. 338 P a g e

Figure 1 Fishbone diagram showing parameters affecting MRR II.EXPERIMENTATION The design of experiments technique is a very powerful tool, which permits to carry out the modeling and analysis of the influence of process variables on the response variables. Improving the MRR and surface quality are still challenging problem that restrict the expanded application of the technology. Semi-empirical models of MRR for various work piece and tool electrode combinations have been presented by various researchers. The influence of pulse current, pulse time, duty cycle, open circuit voltage and dielectric flushing pressure over the MRR and surface roughness on EN 31 tool steel have also been studied. The optimum processing parameters are very much essential to boost up the production rate to a large extent and shrink the machining time, since these materials, which are processed by EDM are costly and the process is very expensive too. The rotatable central composite design is the most widely used experimental design for modeling a second order response surface. A design is called rotatable when the variance of the predicted response at any point depends only on the distance of the point from the center point of design. Table I shows the components of central composite second order rotatable design. Table 1 Components of central composite second order Rotatable Design Variable(K) Fractional Start point 2K Center Point Total (N) Value of α Point(2 k ) 3 8 6 6 20 1.682 4 16 8 7 31 2 5 16 10 6 32 2 6 32 12 9 53 2.378 339 P a g e

As the number of variables is 4, a total of 31 experiments were planned for this investigation. Experiments were carried out using CNC EDM (EMT 43) Electronica die sinking machine. Table II shows the specification of die sinking EDM machine. Table 2 Specifications of die sinking EDM machine Machining conditions Machine Used Electrode Electrode polarity Workpiece Dielectric CNC EDM (EMT 43) (Electronica) Electrolytic Copper ( 99.9% Purity) Positive Oil Hardened Non Shrinking Steel ( 48 50 RC) EDM Oil The composition of AISI D3 steel work-piece material used for experimentation in this work is given in Table 1. Table 3 Chemical composition of AISI D3 steel (wt %) Material C Cr Mn Mo V Si Ni AISI D3 2.05 11.10 0.589 0.042 0.055 0.498 0.065 An electrolytic pure copper with 25 mm X 25 mm is used as a tool electrode (positive polarity).the machining parameters and their levels are shown in Table 4. Table 4 Different variables used in the experiment and their levels coding Level Variable 1 2 3 Pulse On (T on ) in A 200 300 400 μs Pulse Off (T off ) in B 1800 1700 1600 ìs Discharge Current C 8 12 16 (Ip) in A Voltage (V) in V D 40 60 80 The parameter MRR is selected as response variable, which refers to the machining efficiency of the EDM process and defined as follows: MRR (gm/min) = Wi Wf t : Where, Wi = Initial weight of work piece material (gms), Wf= Final weight of work piece material (gms),t = Time period of trials in minutes The work piece is weighed before and after each experiment using an electric balance to determine the value of MRR. For efficient evolution of the EDM process, the larger MRR is regarded as the best machining performance. 340 P a g e

III.METHODOLOGY In statistics, Response surface methodology (RSM) investigates the interaction between several illustrative variables and one or more response variables. Box and Draper [9] were introducing RSM in 1951.The most important purpose of RSM is to use a series of designed experiments to attain an optimal response. A seconddegree polynomial model is use in RSM. These models are only an approximation, but used because such a model is easy to estimate and apply, even when little is known about the process. The process of RSM includes designing of a series of experiments for sufficient and reliable measurement of the response and developing a mathematical model of the second order response surface with the best fittings. Obtaining the optimal set of experimental parameters, thus produce a maximum or minimum value of the response. The Minitab Software was used to analyze the data [10] IV.RESULT AND DISCUSSION Using the experimental results for MRR (Table IV), response surface model is developed and analysis of variance (ANOVA) for the adequacy of the model is then performed in the subsequent step. The F ratio is calculated for 95% level of confidence. Table 5 ANOVA table for MRR (before elimination) Estimated Regression Coefficients for MRR. Term Coefficient SE Coefficient T p Constant 0.22752 0.00508 44.787 0 A 0.000681 0.002744 0.248 0.807 B 0.053207 0.002744 19.393 0 C 0.096158 0.002744 35.049 0 D -0.0431 0.002744-15.71 0 A*A -0.00201 0.002513-0.799 0.436 B*B -0.00335 0.002513-1.331 0.202 C*C 0.012858 0.002513 5.116 0 D*D 0.01777 0.002513 7.07 0 A*B -0.00214 0.00336-0.636 0.534 A*C 0.00251 0.00336 0.747 0.466 A*D 0.000303 0.00336 0.09 0.929 B*C 0.021314 0.00336 6.343 0 B*D -0.013 0.00336-3.868 0.001 C*D -0.01826 0.00336-5.433 0 S = 0.0134406 PRESS = 0.0166486 R-Sq = 99.21% R-Sq(pred) = 95.46% R-Sq(adj) = 98.52% After eliminating the non-significant terms, the final response equation for MRR is given as follows: MRR = 0.222579 + ( 0.000681 X T on ) + (0.053207 X T off )+( 0.096158 X I p )- (0.043101 X V ) 341 P a g e

+ (0.013372 X T on 2 )+ (0.018285 X V 2 )+ (0.021314 X T off X I p )- (0.012999 X T off )- (0.018288 X I p X V) The final model is tested for variance analysis (F-test) and indicates that the adequacy of the test is established which are justified with F-values. For the analysis the data, the checking of goodness of fit of the model is very much required. The model adequacy checking includes the test for significance of the regression model, test for significance on model coefficients, and test for lack of fit. For this purpose, analysis of variance (ANOVA) is performed. Table V shows the ANOVA for the response, MRR. The fit summary recommends that the quadratic model is statistically significant for analysis of MRR. Table 6 Analysis of Variance for MRR Source DF Seq. SS Adj. SS Adj. MS F P Regression 9 0.363588 0.363588 0.040399 244.87 0 Linear 4 0.334449 0.334449 0.083612 506.8 0 Square 2 0.013834 0.013834 0.006917 41.93 0 Interaction 3 0.015304 0.015304 0.005101 30.92 0 Residual 21 0.003465 0.003465 0.000165 Error Lack-of-Fit 15 0.003465 0.003465 0.000231 Pure Error 6 0 0 0 Total 30 0.367052 The normal probability plot of residuals for MRR is illustrated in Fig 2.It is expected that data from experiments form a normal distribution. It reveals that the residual fall on a straight line, implying that the errors are spread in a normal distribution. Here a residual means difference in the observed value (obtained from the experiment) and the predicted value or fitted value. This is also, confirmed by the variations between the experimental results and model predicted values analyzed through residual graphs, and are presented in Fig 3. Figure 2 Normal Probability Figure 3 Residual Plot 342 P a g e

Figure 4 Main Effects Plot From this main effects plot, Figure 4 it is clear that the parameters pulse off and current have highest inclination, so these are most significant but pulse on and voltage are nearly horizontal, so these are nonsignificant. The parametric analysis has been carried out to study the influences of the input process parameter such as Ton, Toff, Ip and V on the process response, MRR during die-sinking EDM process. Contour plots and threedimensional response surface plots are formed based on the quadratic model to evaluate the variation of response. The plots are shown in Figures 5-10.These plots can also give further assessment of the correlation between the process parameters and response as under: 1. MRR increases with increase in Ip and Ton. This is due to higher spark energy from high temperature. (Figure 2. MRR decreases with increase in Toff. 3. Increase in applied voltage also increase MRR. Figure 5 Variation of MRR according to change of Ip and change of Ton Hold value: T off =1700(µs), V= 60 =60 Figure 6.Variation of MRR according to Hold value: Ton=300(µs), V 343 P a g e

Fig 7.Variation of MRR according to change of voltage and current Fig 8. Variation of MRR change of voltage and Hold on Values: Ton=300(μs), Toff= 1700(μs) Hold value: Ip= 12A, Toff=1700(μs) Fig 9.Variation of MRR according to change of Toff according to change of Hold value: Ip= 12A Ton=300(ìs) voltage value: V=60, Ip= 12A Fig 10.Variation of MRR Toff and Hold To check the developed model one confirmation test is carried out at the mid-levels of the process parameters. Table 5and 6 shows the result of the confirmation run for MRR. It is observed that the calculated error is small (about 2%) This confirms the reproducibility of experimental conclusion. Table 7Conformation test result and comparison with predicted result as per model Ton(μs) Toff(ìs) Ip (A) V (Volt) MRR(gm/min) Experimental Model error (%) Predicted 300 1700 12 60 0.22752 0.222579 2.17 Finally an optimum condition is obtained from RSM with an objective of maximum MRR within the experimental range and the levels of the process parameters are Pulse on 500 μs, Pulse off 1500 μs, Current 20 A and Voltage 60 V. 344 P a g e

V. CONCLUSION Experimental investigation on electrical discharge machining of EN 31 tool steel is performed with a view to correlate the process parameters with the responses for MRR. The process has been successfully modeled using response surface methodology (RSM) and model adequacy checking is also carried out. The second-order response models have been validated with analysis of variance. Finally, an attempt has been made to estimate the optimum machining conditions to produce the best possible response within the experimental constraints. This study can help researchers and industries for developing a robust, reliable knowledge base and early prediction of MRR without experimenting with EDM process for EN 31 tool steel. REFERENCES [1] Marafona J.D., Jo A.A. (2009), Influence of work piece hardness on EDM performance, International Journal of Machine Tools & Manufacture, Vol. 49, pp. 744 748. [2] Kunieda M., Lauwers B., Rajurkar K. P., Schumacher B. M. (2005), Advancing EDM through Fundamental Insight into the Process, Journal of Materials Processing Technology, Annals of CIRP, Vol. 54(2), pp. 599-622. [3] Spedding, T.A. and Wang, Z. Q., (1997), Parametric optimization and surface characterization of wire electrical discharge machining process, Precision Engineering, Vol. 20, pp.5-15. [4] Zhang, J.H., Lee, T.C. and Lau, W.S., (1997), Study on the electro-discharge machining of a hot pressed aluminum oxide based ceramic, Journal of Materials Processing Technology, Vol. 63, pp.908-912. [5] Tsai, K.M. and Wang, P.J., (2001), Predictions on surface finish in electrical discharge machining based upon neural network models, International Journal of Machine Tools and Manufacture, Vol. 41, pp.1385 1403. [6] Singh, P. N., Raghukandan, K. and Pai, B.C., (2004), Optimization by Grey relational analysis of EDM parameters on machining Al-10%SiCp composites, Journal of Materials Processing Technology, Vol.155 156, pp.1658 1661. [7] Yih-fong, T. and Fu-chen, C., (2003), A simple approach for robust design of high-speed electrical discharge. Machining technology, International Journal of Machine Tools and Manufacture, Vol. 43, pp.217 227. [8] Cochran, G., and Cox, G.M. (1962), Experimental design, Asia Publishing House, New Delhi. [9] Box, G. E. P. and N.R. Draper, 1987. Empirical Model- Building and Response Surfaces, Jon Wiley & Sons, New York. [10] Minitab14 (2003). Minitab User Manual Release 14. State College, PA, USA. 345 P a g e