A comparative study on modelling material removal rate by ANFIS and polynomial methods in electrical discharge machining process

Chia sẻ: Đăng Xuân Phương | Ngày: | Loại File: PDF | Số trang:15

lượt xem
  Download Vui lòng tải xuống để xem tài liệu đầy đủ

(BQ) The study sheds light on the powerful learning capability of ANFIS models and its superiority over the conventional polynomial models in terms of modelling complex non-linear machining processes

Chủ đề:

Nội dung Text: A comparative study on modelling material removal rate by ANFIS and polynomial methods in electrical discharge machining process

Computers & Industrial Engineering 79 (2015) 27–41<br /> <br /> Contents lists available at ScienceDirect<br /> <br /> Computers & Industrial Engineering<br /> journal homepage:<br /> <br /> A comparative study on modelling material removal rate by ANFIS and<br /> polynomial methods in electrical discharge machining process<br /> Khalid Al-Ghamdi, Osman Taylan ⇑<br /> King Abdulaziz University, Faculty of Engineering, Department of Industrial Engineering, P.O. Box 80204, Jeddah 21589, Saudi Arabia<br /> <br /> a r t i c l e<br /> <br /> i n f o<br /> <br /> Article history:<br /> Received 27 April 2013<br /> Received in revised form 26 October 2014<br /> Accepted 29 October 2014<br /> Available online 7 November 2014<br /> Keywords:<br /> EDM<br /> MRR<br /> Polynomial model<br /> Neuro-fuzzy model<br /> Non-conventional machining<br /> <br /> a b s t r a c t<br /> Due to the controversy associated with modelling Electrical Discharge Machining (EDM) processes based<br /> on physical laws; this task is predominantly accomplished using empirical modelling methods. The<br /> modelling studies reported in the literature deal predominantly with quantitative parameters i.e. ones<br /> with numerical levels. In fact, modelling categorical parameters has been devoted a scant attention. This<br /> study reports the results of an EDM experiment conducted on the Ti–6Al–4V alloy. Its aim was to model<br /> the relationship between the Material Removal Rate (MRR) and the parameters of the process, namely,<br /> current, pulse on-time and pulse off-time along with a categorical factor (electrode material). The modelling process was accomplished using adaptive neuro-fuzzy inference system (ANFIS) and polynomial<br /> modelling approaches. In fact, one purpose of this study was to compare the performance of these modelling approaches as no study was found contrasting their prediction capability in the literature. Regarding the polynomial model, two numerical parameters (current and pulse on-time) were declared<br /> significant in the ANOVA together with the electrode material and its interaction with pulse on-time.<br /> Thus, they were all incorporated in the developed polynomial model. Furthermore, five ANFIS models<br /> with 6, 9, 19, 21 and 51 rules were developed utilizing the first order Sugeno fuzzy approach by backpropagation neural networks training algorithm. Of these, the ANFIS model with 21 rules was the best.<br /> This model also outperformed the polynomial model remarkably in terms of predicting error, residuals<br /> range and the correlation coefficient between the experimental and predicted MRR values. The study<br /> sheds light on the powerful learning capability of ANFIS models and its superiority over the conventional<br /> polynomial models in terms of modelling complex non-linear machining processes.<br /> Ó 2014 Elsevier Ltd. All rights reserved.<br /> <br /> 1. Introduction<br /> Recent years have witnessed remarkable increase in the<br /> demand for the titanium based alloys owing to their high<br /> strength-weight ratio, excellent fracture and corrosion resistance<br /> and superior mechanical properties at elevated temperature. An<br /> extensive utilization of these alloys has been observed in the<br /> aerospace industry (Armendia, Garay, Iriarte, & Arrazola, 2010;<br /> Moiseev, 2005). Indeed, many of the aero-engine components are<br /> produced using these alloys due, in no small part, to their relatively<br /> lightweight which helps to reduce the aircraft’s overall weight<br /> thereby moderating fuel consumption (Singh & Khamba, 2007).<br /> Moreover, titanium alloys are noticeably replacing the traditionally<br /> used materials in aerospace industry such as aluminium (Lütjering<br /> & Williams, 2007). They are also used increasingly in other indus⇑ Corresponding author. Tel.: +966 500031056; fax: +966 26952486.<br /> E-mail addresses: (K. Al-Ghamdi),<br /> (O. Taylan).<br /><br /> 0360-8352/Ó 2014 Elsevier Ltd. All rights reserved.<br /> <br /> trial and commercial applications such as, chemical processing,<br /> pollution control, surgical implantation, nuclear waste storage,<br /> food processing and petroleum refining. A major hindrance for<br /> the wider use of titanium alloys is the difficulties associated with<br /> machining them conventionally. Being chemically reactive to<br /> almost all tool materials, they tend to weld the cutting tools while<br /> machining thereby accelerating tool failure (Hartung, Kramer, &<br /> Von Turkovich, 1982). The damage to cutting tools is compounded<br /> by the fact that the heat generated at the tool-work-piece interface<br /> while machining is concentrated at the tool’s cutting edges. This<br /> heat cannot be transferred through the work-piece due to the poor<br /> thermal conductivity of titanium alloys. Moreover, a very high<br /> mechanical stress is accumulated on the immediate vicinity of<br /> the cutting edge as a result of high strength and hardness that<br /> these alloys retain at elevated temperature poses a further impediment to their machine-ability (Che-Haron & Jawaid, 2005; Trent &<br /> Wright, 2000). Ezugwu and Wang (1997) asserted that it was due<br /> to the limitations associated with machining these alloys that large<br /> companies such as Rolls-Royce and General Electric invested large<br /> <br /> 28<br /> <br /> K. Al-Ghamdi, O. Taylan / Computers & Industrial Engineering 79 (2015) 27–41<br /> <br /> sums of money in developing techniques to deal with them.<br /> Ezugwu, Bonney, Da Silva, and Cakir (2007) argued that, despite<br /> these limitations, machining will remain the favoured option for<br /> manufacturing complex titanium alloys at a competitive cost in<br /> the foreseeable future.<br /> Non-conventional machining processes provide effective alternatives to the conventional ones in dealing with the technical difficulties associated with machining titanium alloys. Of these,<br /> Electrical Discharge Machining (EDM) is the most widely known<br /> and used process for the manufacture of engineering components<br /> (Aspinwall, Soo, Berrisford, & Walder, 2008). In this process, the<br /> material removal takes place as a result of the discharge of energy<br /> between a tool and a work-piece, which are separated by a small<br /> gap filled with a dielectric fluid. The process involves a succession<br /> of discrete discharge pulses initiated using a DC pulse generator for<br /> a certain duration followed by a similar period during which deionisation of the dielectric occurs and the gap is flushed of debris<br /> (McGeough, 1988). A very high temperature is generated by the<br /> discharge energy at the point of the spark on the work-piece and<br /> the material is removed by melting and vaporization. The stability<br /> of the process is maintained by a servo-controlled mechanism<br /> whereby a constant gap between the tool and the work-piece is<br /> sustained whilst impressing the former pre-selected shape<br /> (Singh, Maheshwari, & Pandey, 2004). The EDM theory was postulated by two Soviet scientists in the middle of 1940s (Ho &<br /> Newman, 2003). Its process is very effective in machining materials of any hardness given that they conduct electricity (Abbas,<br /> Solomon, & Bahari, 2007; Yan, Chung, & Yuan Huang, 2005). A<br /> major advantage of the EDM process is its electro-thermal nature<br /> which eliminates the need for any direct contact between the electrode and work-piece during machining thereby eradicating the<br /> mechanical stress, chatter, tool deformation and vibration problems (Tsai & Masuzawa, 2004). A recent review of titanium alloys<br /> EDM is reported in (Kumar, Singh, Batish, & Singh, 2012). The scope<br /> of the current work is centred on the EDM of the Ti–6Al–4V alloy,<br /> which accounts for 50–60% of total titanium alloy production<br /> (Hood, Lechner, Aspinwall, & Voice, 2007).<br /> Material Removal Rate (MRR) is one of the key performance<br /> measures in EDM. The ability to rigorously understand its mechanism and accurately predict its values based on the well-known<br /> physical laws are still debatable (Markopoulos, Manolakos, &<br /> Vaxevanidis, 2008; Tsai & Wang, 2001a). In fact, the large number<br /> of parameters involved and their interdependence together with<br /> the stochastic mechanism of the EDM process appear to significantly impair the ability to establish a scientifically admissible<br /> theory correlating the process variables with the main performance measures such as the MRR. Consequently, this task is<br /> predominantly accomplished by means of empirical or semiempirical modelling methods. Broadly, two approaches are<br /> adopted; namely, the use of conventional statistical models and<br /> the application of soft computing algorithm such as ANFIS. For<br /> instance, the second order (quadratic) polynomials, conventionally<br /> employed in response surface method, were used to relate the<br /> EDM process parameters with its key performance measures in<br /> the case of machining cobalt-bonded tungsten carbide (Puertas,<br /> Luis, & Alvarez, 2004), silicon carbide (Luis, Puertas, & Villa,<br /> 2005), Al–4Cu–6Si alloy–10 wt.% SiCP composites (Dhar, Purohit,<br /> Saini, Sharma, & Kumar, 2007) and DIN 1.2714 steel (Zarepour,<br /> Tehrani, Karimi, & Amini, 2007). Artificial neural network (ANN)<br /> techniques, on the other hand, were used to model and predict<br /> the main EDM performance characteristics in the case of machining aluminium and iron (Tsai & Wang, 2001b) and various steel<br /> grades (Panda & Bhoi, 2005; Pradhan & Biswas, 2010). The performance of ANN was compared with that of the ANFIS in two studies.<br /> In one of these (Puertas et al., 2004), the ANFIS method outper-<br /> <br /> formed the examined ANN techniques in predicting the performance measures whereas comparable performance of the two<br /> approaches was reported by (Pradhan & Biswas, 2010).<br /> Some gaps were noted in the current literature on modelling<br /> die-sinking EDM process. First, while some of the studies compared the performance of ANN with that of ANFIS, none of them<br /> contrasted the prediction capability of polynomial models with<br /> that of ANFIS. This paper attempts to fill this void by comparing<br /> the performance of a polynomial model with a Neuro-Fuzzy<br /> model derived based on different number of rules and fuzzy<br /> inference systems. Second, the parameters examined in published modelling studies were, by and large, quantitative i.e.<br /> their levels were numerical. In fact, scant attention has been<br /> accorded to the case of correlating qualitative (categorical)<br /> parameters with the key performance measures of EDM, particularly when polynomial modelling was adopted. In this study,<br /> the relationship between three numerical parameters along with<br /> a categorical one and MRR was modelled using a polynomial<br /> function as well as a Neuro-Fuzzy system. The implications of<br /> incorporating a qualitative factor and the way in which it was<br /> dealt with are elaborated on. Third, in the context of machining<br /> Ti–6Al–4V, Caydas and Hascalik (2008) appears to be the only<br /> study that addressed the modelling of the EDM process of this<br /> alloy. The present work extends that of Caydas and Hascalik<br /> (2008) in two aspects. First, while Caydas and Hascalik (2008)<br /> examined the influence of three parameters; namely, pulse ontime, pulse off-time and current on the electrode wear and<br /> recast thickness, the current research assessed the impact of<br /> the same three parameters (at different levels) along with the<br /> effect of electrode material on the MRR, which was not examined<br /> in the above mentioned study. Second, although the response surface method was the only modelling approach adopted in Caydas<br /> and Hascalik (2008), this paper utilizes and compares the performance of two methods; namely, polynomial modelling and ANFIS,<br /> bearing in mind that, to date, the latter has not been applied to<br /> model the die-sinking EDM of Ti–6Al–4V.<br /> The methods adopted in the current study are firstly described.<br /> A discussion of the experimental design and procedures follows.<br /> The obtained results are then presented and discussed. After comparing the performance of the investigated modelling approaches,<br /> the paper culminates with the main conclusions.<br /> <br /> 2. Methodology<br /> Design of Experiments (DoE) is an effective and efficient<br /> approach for exploring, understanding and empirically modelling<br /> the cause-and-effect relationship between the engineering processes’ parameters and their performance measures. An experimental investigation involves purposeful changes of the input<br /> variables of a process or system to explore and justify the changes<br /> that may occur in the output response (Montgomery, 2010).<br /> Generally a DOE study is performed in three stages: plan ning,<br /> conducting and analysis and interpretation (Antony, 2003). The<br /> planning stage involves recognizing the problem or the improvement opportunity, stating the objectives, selecting the performance measure(s) and the measurement system(s), determining<br /> the factors that may influence the chosen performance measure(s),<br /> choosing the levels for the factors, finding an appropriate design to<br /> vary the factor levels in accordance with and assigning the factors<br /> and interactions to the selected design. Having accomplished these<br /> tasks, the next stage is to carry out the experiment as planned.<br /> Finally, the obtained results should be analyzed and interpreted.<br /> In so doing, two approaches were adopted in this study; namely,<br /> the polynomial modelling approach and the ANFIS.<br /> <br /> 29<br /> <br /> K. Al-Ghamdi, O. Taylan / Computers & Industrial Engineering 79 (2015) 27–41<br /> <br /> 2.1. Polynomial modelling approach<br /> Polynomial modelling is a major phase of the statistical analysis<br /> of experimental data which comprises four stages: estimation,<br /> hypothesis testing, modelling and identifying the best parameters<br /> setting. The estimation involves obtaining a number to quantify<br /> the impact of each of the studied parameters and their interactions<br /> on the investigated performance measure(s). To accomplish this<br /> task in the current study, a measure that was initially proposed<br /> by Kelley (1935) and labelled the Percentage Contribution (PC)<br /> by Ross (1996) was adopted. This was obtained using Eq. (1).<br /> <br /> PC ¼<br /> <br /> SSP À ðdfpÞMSE<br /> Â 100<br /> SStotal<br /> <br /> ð1Þ<br /> <br /> where SSp is the parameter’s Sum of Squares, dfp is its degrees of<br /> freedom and the MSE is the mean square error of the experiment.<br /> The PC gives the percentage each parameter contributes to the<br /> total performance measure’s variation, and is a measure of how<br /> much the performance could be improved if the parameter was<br /> controlled accurately. The quantities required for estimating the<br /> PC can be obtained from the Analysis of Variance (ANOVA) which<br /> is the main technique for accomplishing the hypothesis testing<br /> task (the process through which the statistical significance of each<br /> of the examined parameter and interaction is tested). The ANOVA<br /> involves decomposing the observed total variance into components due to different sources of variation (factors, interactions<br /> and experimental error). Dividing the estimated variance for each<br /> parameter by that of error lends a random variable that follows<br /> an F-distribution and can be used to test the statistical significance<br /> of the concerned parameter. However, when the experiment is not<br /> replicated, there will be no degrees of freedom to estimate the<br /> error variance. In such cases, it is necessary to identify the factors<br /> or interactions that can reasonably be pooled to provide an estimate of the error variance. A powerful tool for accomplishing this<br /> task is the Half Normal probability Plot (HNP) (Daniel, 1959) of<br /> which the use is conventionally confined to analyzing unreplicated<br /> two-level experiments. Essentially, HNP can be constructed for any<br /> design provided that the plotted effect estimates are uncorrelated<br /> and have the same variance. For three-level designs, Wu and<br /> Hamada (2011) showed how HNP can be used in the case of<br /> decomposing the studied factors into single degree of freedom<br /> components with equal variance. Whitcomb and Oehlert (2007)<br /> proposed an alternative strategy for using the HNP that relaxed<br /> the condition of splitting the effects into single degree of freedom<br /> elements. In this strategy, the sums of squares are used as alternative measures of the factorial effects which are commonly estimated using the contrasts among the averages of response values<br /> associated with the levels of the studied factors. Furthermore, an<br /> iterative method that bears considerable resemblance with the<br /> stepwise variable selection procedure is adopted to estimate a provisional error variance. Dividing the sum of squares of a selected<br /> effect by this estimate renders a chi-square statistic for which a<br /> provisional p-value can be computed and translated in terms of a<br /> provisional z-score of a Half Normal distribution.<br /> The factors and interactions that pronounced significant in the<br /> ANOVA are used to develop an empirical model correlating them<br /> with the response under study (Wu & Hamada, 2011).<br /> Polynomials are very powerful for empirically modelling the<br /> cause-and-effect relationship between the process parameters<br /> (X1, X2..Xm) and the performance measure of interest (Yi). Mathematically, this relationship is developed using Taylor series expansion of terms involving the studied factors and interactions.<br /> Generally, a polynomial model is said to be of order ‘N’ if this is<br /> the power of one or more of its factor terms or if ‘N’ is the sum<br /> of the powers of the factors involved in one of its interaction terms.<br /> Polynomials of second order are commonly used in modelling EDM<br /> <br /> (Puertas et al., 2004; Luis et al., 2005; Dhar et al., 2007; Zarepour<br /> et al., 2007). For ‘m’ factors, a model of this order comprises<br /> (m + 1)(m + 2)/2 terms and takes the form that is given in Eq. (2):<br /> <br /> ^ ^<br /> y ¼ b0 þ<br /> <br /> m<br /> m<br /> m<br /> X<br /> X<br /> X<br /> ^<br /> ^<br /> ^<br /> bi xi þ<br /> bii x2 þ<br /> bij xi xj<br /> i<br /> i¼1<br /> <br /> i¼1<br /> <br /> ð2Þ<br /> <br /> i150 ls), arcing was observed with all of the electrode materials. Consequently, the discharge efficiency was compromised<br /> rendering a noticeable decrease in MRR.<br /> A polynomial model was used to represent the relationship<br /> between the parameters and the interaction that were declared<br /> significant in the ANOVA and the MRR. As a categorical factor,<br /> one way of dealing with the electrode material is to develop a separate model consisting only of the quantitative factors for each of<br /> its levels rendering three separate models. This approach was<br /> adopted by (Sahoo, Routara, & Bandyopadhyay, 2009) in examining<br /> the impact of different work-piece materials. Helpful though in<br /> terms of avoiding explicit modelling of categorical factors, this<br /> approach becomes quite cumbersome as the number of these<br /> factors and their levels increases. In fact, explicit modelling of<br /> categorical factors has the advantage of being more practical as it<br /> yields a single predication equation. Moreover, it renders one estimate of the experimental error variance and more residual degrees<br /> of freedom than fitting several regression models. Therefore, this<br /> approach was adopted in this study.<br /> To model the process under study, a model matrix was constructed from the design array given in Table 2. This is a coded<br /> matrix that has a column for each term in the model to be developed. As three main effects and one interaction of three-level factors were found significant, ten terms were needed. Of these four<br /> corresponded to the linear and quadratic components of the two<br /> quantitative factors: pulse on-time and current. The qualitative<br /> factor (electrode material) has two terms associated with its two<br /> degrees of freedom. Finally, the interaction between the electrode<br /> material and pulse on-time has four terms corresponding to the<br /> four possible combinations of the linear and quadratic elements<br /> of pulse on-time and the two components of the electrode material. To construct the model matrix, it was necessary to code the<br /> levels of the studied parameters so that they become dimensionless and can be used to estimate comparable model coefficients.<br /> Coding also improves the precision with which the coefficients<br /> are estimated (Montgomery, 2010). Linear main effects of the<br /> quantitative parameters were coded so that the highest and lowest<br /> levels of each factors become 1 and À1 respectively. Any value in<br /> between can be calculated so that its difference from 1 and À1 is<br /> proportional to the difference between the corresponding actual<br /> value and the actual highest and lowest values. The codes for the<br /> quadratic main effect can then be obtained by squaring the coded<br /> linear main effect levels. Using this system, a least squares linear<br /> main effect coefficient is interpreted as the average change in the<br /> response variable per unit change in its factor. The quadratic main<br /> effect coefficient, on the other hand, is a measure of how its factor<br /> changes when its own values are changed. The components of the<br /> electrode material factor can be coded by assigning one dummy<br /> variable for each of the three electrode types. For each variable, a<br /> value of one can be assigned to the runs where the appropriate<br /> electrode material was used while a value of zero should be written for the remaining runs. This entails adding three dummy variable columns of which the row-wise sum is one, rendering a<br /> column identical to the one that should be used to estimate the<br /> model intercept. Consequently, it will not be possible to estimate<br /> any of the model coefficients as the model matrix will be singular.<br /> One way out of this problem is to drop one of the three dummy<br /> variables, in which case, each of the coefficients of the remaining<br /> <br />



Đồng bộ tài khoản