This study focuses on the examination of the effect of cutting parameters on surface roughness when drilling medium-density fiberboard (MDF) with a parallel robot. Taguchi technique was applied to find the optimum drilling parameters and, later, the drilling processing. Experimental layout was established using the Taguchi L18 orthogonal array, and experimental data were examined via a statistical analysis of variance (ANOVA). Experimental results were performed by multiple regression analysis (linear and quadratic). Correlation coefficient (R2) was found 99.46% for surface roughness with the quadratic regression model. By the Taguchi analysis, the optimum values for the surface roughness were found to be a point angle of 118°, a cutting speed of 47.1 m/min, and a feed rate of 0.01 mm/rev. The optimization outcomes presented that the Taguchi technique had been successfully performed to decide the optimal surface roughness of the MDF in the drilling.
1. Introduction
Medium-density fiberboard (MDF) consists of wood fibres and a synthetic binder. Diverse machining processes such as drilling, milling, and turning are applied to obtain the required size and form of the MDF [1–4]. The most common technique for examining machining efficiency is based on the surface finishing and covers visual evaluation of surface quality measurement [5, 6]. Nowadays, especially, algorithms based on experimental design are used to optimize various engineering problems [7, 8]. The Taguchi method is an empirical design procedure that greatly reduces the number of experiments using orthogonal arrays and, at the same time, tries to minimize the effects of uncontrollable factors. Optimum process parameters are obtained by a signal-to-noise (S/N) ratio and analysis of variance (ANOVA) [9–12]. Gaitonde [11] researched the impact of cutting parameters on multiple performance characteristics throughout drilling of the MDF material via the Taguchi method. Both, the feed rate and cutting speed were found to be significant drilling parameters to control delamination using ANOVA. Palanikumar [13] examined the delamination of the MDF board using orthogonal array, S/N, and ANOVA. The results showed that the feed rate was the main factor affecting the delamination of the MDF board throughout the drilling process. The authors recommended high spindle speed, low feed rate, and low drill diameter when drilling the MDF boards. Prakash et al. [14] examined the Taguchi method and response surface methodology (RSM) for surface roughness analysis in drilling of medium-density fiberboard panels. Experimental results showed that the most important drilling parameter for surface roughness is feed rate and then cutting speed. Prakash et al. [15] studied the estimation of surface roughness in machining of the MDF material by drilling parameters such as spindle speed, feed rate, and drill diameter using the Taguchi method. Palanikumar and Valarmathi [16] used the Taguchi method, RSM, and ANOVA to estimate the effect of machining parameters on thrust force in drilling of MDF materials. Szwajka et al. [17] examined the influence of the coating type on thrust force, cutting torque, cutting tool temperature, and surface roughness of the hole on the drill cutting blades in drilling MDF panels. The researchers analyzed the results of the experiments statistically via a multifactor analysis of variance (ANOVA).
Many papers have been published on optimization of cutting parameters and estimation of shear forces, surface roughness, and damage factor in machining metals [11 and 12] and composites [18–20]. Köklü et al. [18], drilling performance of functionally graded composite (FGC), depending on delamination, thrust force, and vibration, investigated experimentally by using different cutting parameters, which are feed rate, spindle speed, and material directions (carbon/epoxy and glass/epoxy). They used the Taguchi method, ANOVA, and ANN to optimize parameters. Sur and Erkan [19] drilled carbon fibre reinforced plastic (CFRP) composite plates with standard geometry drill bits with different point angles using different cutting parameters and optimized the hole damages with the Taguchi method. Ay et al. [20] examined the drillability performance of the carbon fibre reinforced polymeric composite material using different shear parameters under dry and cryogenic conditions. They used the Gray relational analysis method to determine optimum drilling conditions. In machining, vibrations occur due to changes in parameters such as cutting speed, feed rate, and depth cut applied to the workpiece. Nonlinear interaction between the tool and workpiece leads to the formation of chip formation and consequently to irregularities on the workpiece surface and directly affects the surface roughness quality [21]. Many properties of the MDF have been investigated in numerous studies in the past, but no work has been done on drilling with a parallel robot. The paper combines the Taguchi method and multiple regressions to optimize process parameters in MDF drilling for minimizing surface roughness. Nonetheless, this study was inspired by very limited or no work on the application of the Taguchi method and multiple regressions in this field. For this reason, this work aims to implement a Taguchi L18 orthogonal array to make changes as follows: (1) a systematic research of the impact of drilling parameters such as point angle, cutting speed, and feed rate on the MDF surface roughness; (2) optimize these process variables considered to find the lowest drilling surface roughness.
2. Materials and Methods
The experiments were realized on a parallel robot with a 0.75 kW motor spindle and a maximum spindle speed of 24,000 rpm [22]. A MDF material with thickness of 8 mm was used as the test sample. The MDF is supplied and produced by Turkey. The MDF is commercially convenient and is operated for producing panel doors, flush doors, tables, cabinets, etc. Mechanical and physical properties of the MDF material are demonstrated in Table 1.
Mechanical and physical properties of the MDF material.
Properties
Value
Tensile strength
0.8 MPa
Modulus of rupture
28 MPa
Elasticity modulus
2800 MPa
Humidity
5–8%
Density
600–900 kg/m3
The experiments realized in MDF materials using HSS (high-speed steel) drills with 5 mm of diameter and point angles of 118° and 135°. The MDF material in the form of 50 × 50 × 8 mm was used as the workpiece. The experimental setup is shown in Figure 1.
The experimental setup.
The surface roughness was measured with the Mahr Marsurf PS 10. The surface roughness (Ra) value was assigned as the mean of three measurements taken from the drilled surfaces. The drilling tests were implemented at three different cutting speeds (15.7, 31.4, and 47.1 m/min), three feed rates (0.01, 0.03, and 0.05 mm/rev), and two point angles (118° and 135°).
3. Experimental Design3.1. Taguchi Method
The Taguchi method improved by Taguchi has been used in engineering analysis to optimize performance characteristics by adjusting design parameters [9]. Also, the Taguchi method offers a simple and effective systematic approach to optimize designs in terms of quality, performance, and cost. By using the Taguchi method based on orthogonal arrays, time and cost can be reduced. Taguchi recommends using signal-to-noise (S/N) ratio (η) to determine the quality characteristics applied to engineering design problems. There are generally three types of quality characteristics; the lower the better, the nominal the best, and the higher the better [23]. In this work, the lower the better quality characteristic was used to minimize surface roughness:(1)η=SN=−10log1N∑i=1nyi2,where n is the number of data of the experiment and yi is the measured data at the ith experiment [24]. In addition to the signal-to-noise ratio, a statistical analysis of variance (ANOVA) was employed to demonstrate the impact of process parameters on surface roughness. Drilling parameters and their levels considered in this work are demonstrated in Table 2. The empirical layout for the three cutting parameters via full factorial design by orthogonal array of Taguchi L18 (21 × 32) is demonstrated in Table 3.
Drilling parameters and their levels.
Parameters
Symbol
Level 1
Level 2
Level 3
Point angle
A
118
135
—
Cutting speed (m/min)
B
15.7
31.4
47.1
Feed rate (mm/rev)
C
0.01
0.03
0.05
Full factorial design with orthogonal array of Taguchi L18 (21 × 32).
Experiment no.
Factor A
Factor B
Factor C
1
1
1
1
2
1
1
2
3
1
1
3
4
1
2
1
5
1
2
2
6
1
2
3
7
1
3
1
8
1
3
2
9
1
3
3
10
2
1
1
11
2
1
2
12
2
1
3
13
2
2
1
14
2
2
2
15
2
2
3
16
2
3
1
17
2
3
2
18
2
3
3
3.2. Multiple Regression Method
Regression analysis proposed by Galton examines the statistical relationship between one or more independent variables and dependent variables [25]. In addition, multiple regressions are a collection of mathematical and statistical techniques beneficial for modeling and analyzing problems in which the response of interest is affected by several variables [26]. In this study, the empirical data was used to construct a linear and quadratic model by the regression method. The independent variables were point angle (α), cutting speed (V), and feed rate (f), whereas the dependent variable was surface roughness (Ra). Multiple regression analysis was operated to obtain the predict equation (the linear and quadratic model) for Ra, as seen in equations (2) and (3), respectively [27]:(2)yl=β0+β1α+β2V+β3f,(3)yq=β0+β1V+β2f+β3α+β4V2+β5f2+β6α2+β7Vf+β8Vα+β9fα,where β0, β1, …, β9 are coefficients used to predict in the multiple regression. The experimental design and statistical analysis were applied by Minitab 18 software.
4. Results and Discussion4.1. Analysis of the S/N Ratio
The S/N ratios were calculated using equation (1) for each of the 18 trials, and the values are presented in Table 4 along with their experimentally measured values. Valarmathi et al. [28] reported that when drilling MDF panels with drills with different point angles, the surface roughness value decreases by selecting the low tip angle, low feed rate, and high cutting speed. Similarly, when Table 4 is studied, it can be seen that the drill bits using a point angle of 118° generally exhibited lower surface roughness values than the drill bits using a point angle of 135°. The results demonstrate that the surface roughness increases with the increase of feed rate, whereas the increase of cutting speed reduces the surface roughness.
The results of experiments and S/N ratio values.
Experiment no.
Control factors
Surface roughness, Ra (μm)
S/N ratio for Ra (dB)
A
B
C
Point angle (α)
Cutting speed (V)
Feed rate (f)
1
118
15.7
0.01
3.08
−-9.7710
2
118
15.7
0.03
4.072
−12.1962
3
118
15.7
0.05
4.862
−13.7363
4
118
31.4
0.01
2.479
−7.8855
5
118
31.4
0.03
3.269
−10.2883
6
118
31.4
0.05
4.186
−12.436
7
118
47.1
0.01
2.122
−6.5349
8
118
47.1
0.03
2.883
−9.1968
9
118
47.1
0.05
3.486
−10.8465
10
135
15.7
0.01
3.837
−11.6798
11
135
15.7
0.03
4.886
−13.7791
12
135
15.7
0.05
5.864
−15.3639
13
135
31.4
0.01
3.249
−10.235
14
135
31.4
0.03
3.975
−11.9867
15
135
31.4
0.05
4.758
−13.5485
16
135
47.1
0.01
3.029
−9.6259
17
135
47.1
0.03
3.43
−10.7059
18
135
47.1
0.05
4.212
−12.4898
The average value of the surface roughness was found to be 3.759 μm at the end of the drilling tests. The same, average value of S/N ratio for surface roughness was found to be −11.2392 dB. Analysis of the effect of each control factor (Ct, V, and f) on the surface roughness has been implemented with the S/N ratio response table in Table 5. The level values of control factors for Ra depicted in Table 5 are demonstrated in graph forms in Figure 2. These figures demonstrate the level to be chosen for the optimum drilling parameters (the level with the highest point on the graphs) and, at the same time, the relative influence each parameter has on the S/N ratio. The best surface roughness values were specified as factor A1 (S/N = −10.321), factor B3 (S/N = −9.900), and factor C1 (S/N = −9.289). The optimal values for the surface roughness assigned for the point angle, cutting speed, and feed rate were 118°, 47.1 m/min, and 0.01 mm/rev, respectively (Figure 2). As seen in the S/N ratio impact graphs, the feed rate is demonstrated to have a potent influence on surface roughness and its S/N ratios.
The S/N response table for Ra factor.
Level
Control factors
Surface roughness (Ra)
A
B
C
Level 1
−10.321
−12.754
−9.289
Level 2
−12.157
−11.063
−11.359
Level 3
—
−9.900
−13.070
Delta
1.836
2.854
3.781
The effect of process parameters on average S/N ratio for Ra.
4.2. Analysis of Variance
The analysis of variance (ANOVA) is used to determine the effects of design parameters on performance characteristics. In this work, the application of ANOVA, the percent contribution ratio (PCR) of surface roughness is investigated, and the ANOVA conclusions for surface roughness are demonstrated in Table 6. This analysis was performed at a 5% (P<0.05) significance level and a 95% confidence level to identify the drilling parameters that affect the surface roughness. For the control factors to be statistically significant on the output, the P value should be less than 0.05, and the control factor with the highest F value is considered to be the most effective factor [29]. When Table 6 is examined, the PCR of the A, B, and C factors on the surface roughness reached to be 16.89%, 30.99%, and 50.18%, respectively. The feed rate (factor C, 50.18%) was determined as the most significant factor affecting surface roughness.
Results of ANOVA for surface roughness.
Variance source
Degree of freedom (DoF)
Sum of squares (SS)
PCR (%)
Mean of squares (MS)
F ratio
P value
Ra
A
1
2.5696
16.89
2.56964
104.21
0.001
B
2
4.7151
30.99
2.35757
95.61
0.001
C
2
7.6358
50.18
3.81788
154.83
0.001
Error
12
0.2959
1.94
0.02466
Total
17
15.2165
100.00
Bold values show the most influential factor.
4.3. Multiple Regression Analysis of Surface Roughness
In this study, multiple regression analysis for both linear and quadratic regression models was used for the purpose of estimating surface roughness (Ra) under the effect of parameters such as point angle (α), cutting speed (V), and feed rate (f). The linear regression model equation obtained for surface roughness is given below:(4)Ral=−1.820+0.04445α−0.03949V+39.88f.
The comparison of empirical and estimated values obtained with the linear regression model is specified in Figure 3. R2 value of the equation obtained by the linear regression model for Ra was found to be 97.37%. The quadratic regression model equation obtained for surface roughness is given below:(5)Raq=−2.96−0.0344V+62.3f+0.05413α+0.000653V2+28f2−0.502Vf−0.000245Vα–0.066fα.
Comparison of surface roughness between predicted and experimental values by the linear regression model.
The comparison of empirical and estimated values obtained with the quadratic regression model is specified in Figure 4. The R2 value of the equation obtained by the quadratic regression model for Ra was reached to be 99.46%. As a result, the quadratic regression model gave better results than the linear regression model and proved to be successful in predicting surface roughness.
Comparison of surface roughness between predicted and experimental values by the quadratic regression model.
4.4. Optimization of Surface Roughness
Optimized parameters must be verified to estimate optimal surface roughness [29]. For this purpose, a confirmation experiment by using the Taguchi optimization technique was performed. Equation (6) was used in order to find the optimum surface roughness [30]:(6)Raopt=A1−TRa+B3−TRa+C1−TRa+TRa,where (A1, B3, and C1) show the optimal-level mean values of surface roughness (Raopt) in Table 2 and TRa declares the mean of whole of the surface roughness value obtained from Table 7. According to these equations, the minimum surface roughness can be found as 2.022 μm.
The mean response table for Ra factor.
Level
Control factors
Surface roughness (Ra)
A
B
C
Level 1
3.382
4.433
2.966
Level 2
4.138
3.653
3.752
Level 3
—
3.194
4.561
Delta
0.756
1.240
1.595
The results of the confirmation experiment performed with Table 7 are evaluated, taking into account the confidence interval (CI) value calculated by equations (7) and (8) [30]:(7)CIRa,a=Fα,1,feVe1neff+1R,(8)neff=N1+Tdof,where Tdof is the total main factor degrees of freedom, N is the total number of experiments, neff is the effective number of replications, R is the number of replications for confirmation experiments, Ve is error variance, fe is the degrees of freedom of error, α is the significance level, Fα,1,fe is the F ratio at a 95 percent CI. F0.05, 1, 12 = 4.747 (from F test table), VeRa = 0.02466 (Table 6), R = 3, N = 18, Tdof = 5, and neff = 3. The confidence intervals were calculated as CIRa = ±0.279. The mean optimal Ra (CI 95 percent) is given as follows:(9)Raopt−CIRa<Raexp<Raopt−CIRa=2.022−0.279<2.122<2.022+0.279=1.743<2.122<2.301.
The empirical value of Raexp was within admissible CI limits. As a result, the system for surface roughness was successfully optimized by the Taguchi technique (P=0.05). Estimated values and confirmation test conclusions by Taguchi technique and regression values are specified in Table 8. The estimated values and empirical values are very close to each other. When the confirmation test results are studied, it can be said that the conclusions obtained are sufficient and Taguchi optimization is successful.
Estimated values and confirmation test results by Taguchi technique and regression values.
Level
Taguchi method
Linear regression method
Quadratic regression method
Exp.
Pred.
Error (%)
Exp.
Pred.
Error (%)
Exp.
Pred.
Error (%)
Ra (μm)
A1B3C1 (optimum)
2.122
2.012
5.18
2.122
1.964
7.44
2.122
2.201
3.72
A2B2C3 (random)
4.758
4.834
1.59
4.758
4.935
3.72
4.758
4.820
1.30
5. Conclusions
In this paper, the MDF material is drilled via different cutting parameters and drill bits. The experimental layout was established according to the Taguchi L18 orthogonal array. The Taguchi technique is applied to obtain the optimal drilling parameters and, later, the drilling and experimental conclusions were examined via ANOVA. The conclusions of the experiment are given below:
By using S/N, the optimal values for the surface roughness assigned for the point angle, cutting speed, and feed rate were 118°, 47.1 m/min, and 0.01 mm/rev, respectively (A1B3C1).
According to the variance analysis, the PCR of the A, B, and C factors on the surface roughness were found to be 16.89%, 30.99%, and 50.18%, respectively. The feed rate (factor C, 50.18%) was determined as the most important factor affecting surface roughness.
The R2 values of the equations attained via the quadratic regression model for Ra was found to be 99.46%. As a result, the quadratic regression model gave better results than the linear regression model and proved to be successful in predicting surface roughness.
The optimization outcomes showed that the Taguchi technique had been successfully performed to decide the optimal surface roughness of the MDF material in the drilling operation.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request and within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
PrakashS.MercyJ. L.SaluguM. K.VineethK. S. M.Optimization of drilling characteristics using grey relational analysis (GRA) in medium density fiber board (MDF)201524-51541155110.1016/j.matpr.2015.07.0802-s2.0-84947334086BenlatrecheY.NouveauC.MarchalR.MartinsJ. P. F.AknoucheH.Applications of CrAlN ternary system in wood machining of medium density fibreboard (MDF)20092675–81056106110.1016/j.wear.2009.01.0442-s2.0-65949103207KocK. H.ErdinlerE. S.HazirE.ÖztürkE.Effect of CNC application parameters on wooden surface quality2017107121810.1016/j.measurement.2017.05.0012-s2.0-85019071930DipponJ.RenH.AmaraF. B.AltintasY.Orthogonal cutting mechanics of medium density fiberboards2000507-82530SouzaA. J.JeremiasT. D.GonzalezA. R.AmorimH. J.Assessment of melamine-coated MDF surface finish after peripheral milling under different cutting conditions201977455956810.1007/s00107-019-01413-y2-s2.0-85064568268AyyildizM.Modeling for prediction of surface roughness in milling medium density fiberboard with a parallel robot201939571672310.1108/sr-02-2019-00512-s2.0-85071309662SofuogluS. D.Determination of optimal machining parameters of massive wooden edge glued panels which is made of Scots pine Pinus sylvestris L.) using Taguchi design method2017751334210.1007/s00107-016-1028-z2-s2.0-84961205800HazirE.KocK. H.Optimization of wood machining parameters in CNC routers: Taguchi orthogonal array based simulated angling algorithm2019214493510TaguchiG.ElsayedE.HsiangT.1989New York, NY, USAMcGraw-HillKaraF.Taguchi optimization of surface roughness and flank wear during the turning of DIN 1.2344 tool steel2017591090390810.3139/120.1110852-s2.0-85043597947KurtM.KöklüU.Minimization of the shape error in the interrupted grinding process by using Taguchi method2012186677682KöklüU.Optimisation of machining parameters in interrupted cylindrical grinding using the Grey-based Taguchi method2013268696702GaitondeV. N.KarnikS. R.DavimJ. P.Taguchi multiple-performance characteristics optimization in drilling of medium density fibreboard (MDF) to minimize delamination using utility concept20081961–3737810.1016/j.jmatprotec.2007.05.0032-s2.0-36549015127PalanikumarK.PrakashS.ManoharanN.Experimental investigation and analysis on delamination in drilling of wood composite medium density fiber boards2019241213411348PrakashS.MercyJ. L.GoswamiK.A systemic approach for evaluating surface roughness parameters during drilling of medium density fiberboard using Taguchi method20147111888189410.17485/ijst/2014/v7i11.21PalanikumarK.ValarmathiT. N.Experimental investigation and analysis on thrust force in drilling of wood composite medium density fiberboard panels201640139140010.1007/s40799-016-0044-62-s2.0-84965149071SzwajkaK.Zielińska-SzwajkaJ.TrzepiecinskiT.Experimental study on drilling MDF with tools coated with TiAlN and ZrN201912338610.3390/ma120303862-s2.0-85060619287KöklüU.MaydaM.MorkavukS.AvcıA.DemirO.Optimization and prediction of thrust force, vibration and delamination in drilling of functionally graded composite using Taguchi20196808533510.1088/2053-1591/ab26172-s2.0-85069680549SurG.ErkanÖ.Cutting tool geometry in the drilling of CFRP composite plates and Taguchi optimisation of the cutting parameters affecting delamination2018363619628AyM.BasmaciG.Said YorukA.EtyemezA.Optimization of drilling of carbon fiber reinforced polymer with the grey-based method2020137453553810.12693/aphyspola.137.535LitakG.Chaotic vibrations in a regenerative cutting process20021371531153510.1016/s0960-0779(01)00176-x2-s2.0-0036604641Aşkar AyyıldızE.AyyıldızM.Development of a 3-axis parallel kinematic machine for milling wood material–Part 1: design201712493269337PhadkeM. S.1989Englewoodcliffs, NJ, USAPrentice-HallKaraF.KöklüU.KabasakaloğluU.Taguchi optimization of surface roughness in grinding of cryogenically treated AISI 5140 steel202062101041104710.3139/120.111583GaltonF.Regression towards mediocrity in hereditary stature18861524626310.2307/2841583MontgomeryD. C.1984New York, NY, USAJohn Wiley2ndAsiltürkI.ÇunkaşM.Modeling and prediction of surface roughness in turning operations using artificial neural network and multiple regression method201138558265832ValarmathiT. N.PalanikumarK.SekarS.Modeling of surface roughness in drilling of MDF panels2015766-76783183610.4028/www.scientific.net/amm.766-767.831KaraF.ÖztürkB.Comparison and optimization of PVD and CVD method on surface roughness and flank wear in hard-machining of DIN 1.2738 mold steel2019391243310.1108/sr-12-2017-02662-s2.0-85043754297GünayM.YücelE.Application of Taguchi method for determining optimum surface roughness in turning of high-alloy white cast iron2013462913919