# Prediction of Crack for Drilling Process on Alumina Using Neural Network and Taguchi Method

Authors: Kingsun Lee
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

This study analyzes a variety of significant drilling conditions on aluminum oxide (withL18 orthogonal array) using a diamond drill. The drilling parameters evaluated are spindle speed, feed rate, depth of cut, and diamond abrasive size. An orthogonal array, signal-to-noise (S/N) ratio, and analysis of variance (ANOVA) are employed to analyze the effects of these drilling parameters. The results were confirmed by experiments, which indicated that the selected drilling parameters effectively reduce the crack. The neural network is applied to establish a model based on the relationship between input parameters (spindle speed, feed rate, depth of cut, and diamond abrasive size) and output parameter (cracking area percentage). The neural network can predict individual crack in terms of input parameters, which provides faster and more automated model synthesis. Accurate prediction of crack ensures that poor drilling parameters are not suitable for machining products, avoiding the fabrication of poor-quality products. Confirmation experiments showed that neural network precisely predicted the cracking area percentage in drilling of alumina.


---

## Body

## 1. Introduction

Alumina is the more common name of aluminum oxide (Al2O3) and combination of hardness, high temperature operation, and good electrical insulation makes it useful for a wide range of applications, that is, semiconductor, electronics, and medical and automotive applications. Ceramic drilling is usually achieved by laser, but the use of laser drilling often creates splashes, recast layers, and a tapered bore [1].Grinding is the most common abrasive process. It is a material removal process accomplished by the application of a bonded grinding wheel—comprised of abrasive particles that are contained in a bonded grinding wheel—applied to a work piece surface while rotating at very high speeds. The typical ceramic is characterized by high hardness and low toughness, along with lack of plasticity and ductility. Following contact with abrasive particles, cracks will be produced in the ceramic. The cutter shape, depth of cut, material property, cutting speed, and abrasive particle are the key points of the grinding process [2].Funkenbusch et al. [3] used a diamond wheel to grind single crystal alumina (sapphire). They found that increased frequency of use of a diamond grinding wheel resulted in larger grinding force, but grinding efficiency was reduced. Finer abrasive particles resulted in good surface roughness, but the better grinding efficiency of a diamond grinding wheel will reduce surface roughness.Li et al. [4] used the diamond tool grinding to get the material removed efficiently, the diamond should be tightly integrated on the substrate, and metal bonded wheels are measured at higher bonding strength than resin bonded wheels. The optimal depth of cut can achieve the maximum rate of material removal efficiency, where the cut depth is higher or lower than the optimal depth; however, material removal efficiency will be gradually decreased.Mohan et al. [5] employed a signal-to-noise ratio to analyze the influence of various parameters on peel-up and push-down delamination factors in drilling of glass fiber-reinforced plastic (GFRP). Using Taguchi and response surface methodology, the main objective was to determine the factors that influence delamination and to achieve optimal machining conditions.Gill and Singh [6] demonstrated the use of an Adaptive Neurofuzzy Inference System technique of fuzzy-based systems to simulate the material removal rate in ultrasonic drilling of alumina ceramic. The performance of the model was validated by comparing the predicted results with practical results obtained from corroboration experiments. Davim [7] presented a study of the influence of cutting parameters and cutting time on drilling metal-matrix composites. The objective was to establish a correlation between cutting parameters and evaluators.Gaitonde et al. [8] presented the methodology of Taguchi optimization method for simultaneous minimization of delamination at drilling entry and exit holes while drilling a melamine coating layered (MDF) panel. It was revealed that the delamination can be effectively reduced in drilling of MDF materials by employing a higher cutting speed and lower feed rate. Janardhan and Krishna [9] applied CNC cylindrical grinding machine with L9 orthogonal array to investigate an optimal metal removal rate (MRR) and surface roughness (Ra). The results revealed that feed rate and depth of cut influence predominantly the output responses of metal removal rate (MRR) and surface roughness (Ra). The results are further confirmed by conducting confirmation experiments.The optimal control factors for the hole quality were determined by using Taguchi technique [10]. Two cutting tools, cutting speeds, and feed rates were considered as control factors, and L8 orthogonal array was determined for experimental trials. Multiple regression analysis was employed to derive the predictive equations of the surface roughness and roundness error achieved via experimental design. Confirmation experiments showed that Taguchi method precisely optimized the drilling parameters in drilling of stainless steel.Abductive network draws on statistical and multiple regression analysis methods as well as neural networks, resulting in faster and more automated model synthesis. The development of this approach for network learning through self-organization has followed the track of the group method of data handling (GMDH) algorithm and the closely related adaptive learning network (ALN) technique [11]. In this study, abductive modeling technique is able to represent the complex and uncertain relationships between crack and process parameters. It shows that the crack can be predicted with reasonable accuracy based on the developed abductive network.


## 2. Taguchi Method and Experimental Design

Taguchi methods are statistical methods developed by Genichi Taguchi to improve the quality of manufactured products. Robust design method helps ensure customer satisfaction. Robust design focuses on improving the fundamental function of the product or process, thus facilitating flexible designs and concurrent engineering. Indeed, it is the most powerful method available to reduce product cost, improve quality, and simultaneously reduce development interval. Design of experiment (DOE) parameter design techniques are used to determine which controllable factors and which noise factors are the significant variables. The aim is to set the controllable factors at those levels that will result in a product or process being robust with respect to the noise factors [12, 13].Taguchi’s approach is based on statistical design of experiments, which fulfills the requirements to solve engineering problems and enhance process optimization. The Taguchi method uses a practical design of orthogonal arrays and a relatively small number of experiments to find one set of optimization. This study involved four factors: spindle speed, feed rate, depth of cut, and the grain size of the abrasive diamond. Each factor was used at three levels. Table1 lists the experimental factors and levels. L18 orthogonal array is selected (Table 2) where 18 experiments were carried out to study the effect of drilling parameters. Each experiment was repeated twice in order to reduce experimental errors.Table 1
Experimental levels and factors.
Factor
Levels
1
2
3
A, spindle speed (rev/min)
30000
25000
20000
B, feed rate (mm/min)
10
20
30
C, depth of cut (mm)
0.02
0.04
0.06
D, abrasive grain size (#)
100
140
180Table 2
L
1
8 orthogonal array.
Trial number
A
B
C
D
1
1
1
1
1
2
1
2
2
2
3
1
3
3
3
4
2
1
1
2
5
2
2
2
3
6
2
3
3
1
7
3
1
2
1
8
3
2
3
2
9
3
3
1
3
10
1
1
3
3
11
1
2
1
1
12
1
3
2
2
13
2
1
2
3
14
2
2
3
1
15
2
3
1
2
16
3
1
3
2
17
3
2
1
3
18
3
3
2
1The Taguchi method experiments were used to investigate the relationship between process parameters and crack area in this study. These experiments provide a proper dataset to train the associated abductive network model. Once the abductive network model was constructed, the relationships between input and output variables became available.


## 3. Abductive Network Synthesis and Evaluation

The abductive neural network analysis method has been used for simulation with the aid of a program written in C-language. It has been shown that the prediction accuracy in abductive neural networks [14] is much higher than other networks. Abductive neural network analysis based on the abductive modeling technique is able to represent the complex and uncertain relationships between input and output parameters. In an abductive network, a complex system can be decomposed into smaller, simpler subsystems grouped into several layers using polynomial function nodes. The polynomial network proposed by Ivakhnenko [15] is a group of methods of data handling (GMDH) techniques. These nodes evaluate the limited number of inputs by a polynomial function and generate an output to serve as an input to subsequent nodes of the next layer. The general polynomial function in a polynomial functional node can be expressed as follows:(1)y0=B0+∑i=1nBixi+∑i=1n∑j=1nBijxixj+∑i=1n∑j=1n∑k=1nBijkxixjxk+⋯,where xi, xj, xk are the inputs, y0 is the output, and B0, Bi, Bij, Bijk are the coefficients of the polynomial functional nodes.In this paper, several specific types of polynomial functional nodes were used in a polynomial network to predict the crack of the alumina. These polynomial functional nodes were called normalizer (N), unitizer (U), white node (W), single node (S), double node (D), and triple node (T). These are defined as follows.(1) Normalizer. A normalizer transforms the original input variables into a relatively common region:(2)a1=q0+q1x1,where a1 is the normalized input, q0 and q1 are the coefficients of the normalizer, and x1 is the original input.(2) White Node. A white node consists of the linear weight sums of all the outputs of the previous layer:(3)b1=r0+r1y1+r2y2+r3y3+⋯+rnyn,where y1,y2,y3,…,yn are the inputs of the node, b1 is the output of the node, and r0,r1,r2,r3,…,rn are the coefficients of the white node.(3) Single, Double, and Triple Node. These are named according to the numbers of the input variables. The algebraic form of each of these nodes is shown in the following equation:(4)single:c1=s0+s1z1+s2z12+s3z13,double:d1=t0+t1n1+t2n12+t3n13+t4n2+t5n22+t6n23+t7n1n2,triple:e1=u0+u1o1+u2o12+u3o13+u4o2+u5o22+u6o23+u7o3+u8o32+u9o33+u10o1o2+u11o2o3+u12o1o3+u13o1o2o3,where z1,z2,z3,…,zn, n1,n2,n3,…,nn, and o1,o2,o3,…,on are the inputs of the node; c1, d1, and e1 are the outputs of the node; and s0,s1,s2,s3,…,sn, t0,t1,t2,t3,…,tn, and u0,u1,u2,u3,…,un are the coefficients of the single, double, and triple nodes.These nodes are third-degree polynomial equations. The double node and triple node have cross-terms, allowing interaction among the node input variables.(4) Unitizer. A unitizer coverts the output to a real output:(5)f1=v0+v1i1,where i1 is the output of the network, f1 is the real output, and v0 and v1 are the coefficients of the unitizer.To build a complete abductive network, the first requirement is to train the database. During training, the information given by the input and output parameters must be sufficient. A predicted square error (PSE) criterion is then used to automatically determine an optimal structure. The principle of the PSE criterion is to select the least complex but still accurate network as possible. The PSE is composed of two terms:(6)PSE=FSE+KP,where FSE is the average square error of the network for fitting the training data and KP is the complex penalty of the network, as shown in the following equation: (7)KP=CPM2σp2KN,where CPM is the complex penalty multiplier, K is the number of coefficients in the network, N is the number of training datasets to be used, and σp2 is a prior estimate of the model error variance.


## 4. Experiment and Analysis

A high-speed CNC milling machine was used to drill alumina. Crack areas were measured by tool microscope and crack images were taken by a digital microscope.The experimental specimens used in this study were alumina (Al2O3) ceramic. The specifications are listed in Table 3. Each specimen, with thickness of 3 mm, was placed on a moveable platform of fixture, as shown in Figure 1.Table 3
Ceramic specifications.
Item
Specification
Material
99.7%Al2O3
Size (mm3)
26 × 33 × 3.2
Hardness (HV)
1,700
Bending strength (Kgf/mm2)
3,500
Compression strength (Kgf/mm2)
30,000Figure 1
Drilling platform.4 mm diameter diamond grinding rods were used in the drilling experiment. The substrate of the diamond grinding rod is SCM440. The diamond particles were electroformed on the grinding rod.The crack generated during drilling was measured via microscope. The image was amplified with 150x, and Micro CAM measurement software was applied to evaluate cracking area. The drilling cracking percentage around the holes is determined by the ratio of the crack area (Ac) at the exit side to the ideal whole area (A0). Figure 2 showed the scheme of crack area after drilling. The value of the cracking percentage Pc can be expressed as follows:(8)Cracking  percentagePc=AcA0100%.Figure 2
Cracking percentage scheme.


## 5. Results and Discussion

Ceramics offer many advantages compared to other materials. They are harder and stiffer than steel. Small cracking area at the exit side of ceramics denotes better quality. A set of experiments using the Taguchi method were conducted to investigate the relationship between process parameters and cracking percentage. The “smaller the better” quality characteristic for minimum crack should be used to obtain optimal drilling performance. Optimal process design is achieved when theS/N ratio is maximized. Table 4 gives the calculated S/N ratio for the ceramic crack area according to L18 orthogonal array. Figure 3 shows the cracking condition of trial number 1.Table 4
Experimental results for cracking percentage.
Trial number
A
B
C
D
Cracking area percentage (%)
S
/
N ratio
y
1
y
2
1
1
1
1
1
11.496
11.109
-
21.065
2
1
2
2
2
17.925
18.721
-
25.262
3
1
3
3
3
102.952
104.514
-
40.319
4
2
1
1
2
16.219
15.279
-
23.949
5
2
2
2
3
18.135
16.563
-
24.794
6
2
3
3
1
41.316
39.838
-
32.167
7
3
1
2
1
14.158
14.687
-
23.182
8
3
2
3
2
38.682
37.350
-
31.601
9
3
3
1
3
25.149
24.122
-
27.833
10
1
1
3
3
88.434
84.236
-
38.726
11
1
2
1
1
7.983
8.494
-
18.321
12
1
3
2
2
33.928
32.161
-
30.385
13
2
1
2
3
25.284
22.371
-
27.558
14
2
2
3
1
9.261
10.583
-
19.951
15
2
3
1
2
15.823
17.460
-
24.434
16
3
1
3
2
39.655
43.122
-
32.345
17
3
2
1
3
19.118
20.347
-
25.908
18
3
3
2
1
30.436
28.187
-
29.347Figure 3
Cracking condition of ceramic.The Taguchi method usesS/N ratio to measure the variations of the experimental design. The “smaller the better” quality characteristic was selected for the calculation of S/N ratio since the lowest values of crack were the desired result in terms of good product quality. A greater S/N ratio corresponds to a small cracking percentage.A bigger control factor level difference results in a greater effect on the crack area at the exit side of alumina. The effect of all control factors and level differences can be calculated, as listed in Table5. S/N ratio response value at level 1 of factor A in Table 5 is the summation of S/N ratio corresponding to factor A at level 1 divided by the number of repetitions of the factor level 1 in Table 4. Factor C (depth of cut) is the most effective process factor, because it has the highest level of difference among all control factors.Table 5
S
/
N ratio response for cracking percentage.
Level
A (spindle speed)
B (feed rate)
C (depth of cut)
D (abrasive grain)
1
-
2
9.01
-
27.80
-
23.59
-
24.01
2
-
25.48
-
24.31
-
26.75
-
28.00
3
-
2
8.37
-
30.75
-
32.52
-
30.86
Effect
3.54
6.44
8.93
6.85
Rank
4
3
1
2The response diagram is more clearly indicating the effect of each control factor, as shown in Figure4. Selected factors with higher S/N ratio values will yield minimum value of crack. It is clear from the figure that the optimum combination of drilling parameters for a minimum crack area is A2, B2, C1, and D1 (i.e., spindle speed = 25000 rpm; feed rate = 20 mm/min; depth of cut = 0.02 mm; and abrasive grain = 100 mesh).Figure 4
Response diagram for cracking percentage.The results of the experiments were evaluated by the analysis of variance (ANOVA). The main objective of the analysis was to determine the influence of every parameter on the variance of the results, regarding the total variance of all the parameters. In this study, ANOVA established the relative significance of factors in terms ofF-ratio and P value to the response. ANOVA for S/N ratio can be calculated, as shown in Table 6.Table 6
ANOVA for cracking percentage.
Sources of variance
Degrees of freedom (DOF)
Sum of squares (SS)
Variance (Var.)
F-ratio
P value
A
2
42.602
21.301
4.12
0.054
B
2
124.78
62.39
12.05
0.003
C
2
246.13
123.065
23.77
0
D
2
142.075
71.038
13.72
0.002
Error
9
46.588
5.176
Total
17
602.175The calculation of ANOVA was made on the basis of the equations in [16]. The total sum of the squared deviations SST is decomposed into two components: the sum of the squared deviations SSm due to each factor and the sum of the squared errors SSe.In ANOVA, total sum of squares (SST) can be calculated as(9)SST=∑i=1n∑j=1rYij2-n∗r∗Y¯2,where Yij is the experimental data, n is the number of experiments, each experiment has r experimental data, and Y¯ is the grand mean.The sum of squares (SS) for each factorm can be expressed as(10)SSm=n∗rL∑k=1LY¯k-Y¯2,where Y¯k is the response values at level k of factor and L is the number of levels:(11)Total  degrees  of  freedom  (DOF)=Total  number  of  experiments-1,DOF  for  each  factor=Number  of  levelsL-1,VarianceVar=SSDOF,F-ratio=VarVar  of  error  (Ve),Pvalue=FDISTF-ratio,DOF  of  factor,  DOF  of  error,where FDIST is Microsoft Excel function.From Table6, it can be observed that speed (Factor A) is insignificant factor due to P value > 0.05. The depth of cut (Factor C) is the most significant drilling factor, indicated by its P value = 0. The P value was used as a tool to check the significance of each of the factors. The smaller value of P was more significant to the response. The response model is significant at the considered confidence level (95%) since the response has P value < 0.05.To check the accuracy of analysis, a confirmation test needs to be carried out. The estimatedS/N ratio for the optimal level of the drilling parameters can be calculated as(12)ηopt=ηavg+∑i=1nηi-ηavg,where ηavg is the mean value of S/N ratio of all control factors, ηi is the mean S/N ratio at the optimal level i, and n is the number of parameters.Two trials were made in the confirmation test; the results are presented in Table7. The table shows the comparison of the estimated S/N ratio with the actual S/N ratio using the optimal parameters. It is seen that there is good agreement between the estimated and actual S/N ratios observed. The error is found to be less than 7.501 per cent. Furthermore, the confirmation test verifies that the obtained optimal drilling parameters provide high experimental reliability.Table 7
Confirmation test results.
Trial number
Cracking percentage (%)
Confirmation testS/N ratio
Estimated valueS/N ratio
Error %
1
5.892
-
15.62
-
14.53
7.501
2
6.183


## 6. Prediction of Abductive Network

The abductive network established the relationship between input parameter and the output parameter during drilling process. To build a complete abductive network, the first requirement is to train the database. Table4 illustrates the data of process parameters and the cracking percentage of product obtained from the drilling process.The three-layer abductive networks created were based on the development of the optimal drilling system model. These were comprised of the drilling system parameters and the drilling results synthesized automatically and are capable of predicting the product crack under various drilling modeling techniques.In the abductive network (appendix), each layer of output had three layers of input (e.g., T2 (Triple2) had T1 (Triple1), T4 (Triple4) and T5 (Triple5) had three layers, T1 had NX2 (NormalizerX2), and NX3 (NormalizerX3) and NX4 (NormalizerX4) had three layers). The input parameters were spindle speed, feed rate, depth of cut, and the grain size of the abrasive diamond. The output result was the cracking percentage. Each node was related to a specific polynomial expression. These polynomial function nodes constituted the complex relationship between the input parameters and the individual output result. If the input parameters are changed, the whole model of this abductive network will change completely. Once the abductive network model was set up, any engineer would have been able to find the cracking area under drilling parameters without any complicated computation required. All polynomial equations used in this network are shown in Figure5, where(13)NX1=-5.95+0.000238X1,NX2=-2.38+0.119X1,NX3=-2.38+59.5X1,NX4=-4.17+0.0298X1,T1=-0.903-0.078NX2+0.687NX3+0.718NX4+0.29NX22+0.443NX32+0.258NX42+0.0439NX2NX3-0.00296NX2NX4+0.324NX3NX4+0.115NX2NX3NX4+0.219NX23-0.25NX43,T4=-0.87+0.449NX1-0.43NX2+0.702NX3+0.296NX12+0.213NX22+0.399NX32-0.0146NX1NX2+0.441NX1NX3+0.0167NX2NX3+0.044NX1NX2NX3-0.125NX13+0.312NX23,T5=-0.682+0.342NX1+0.0866NX2+0.737NX4+0.35NX12+0.0843NX22+0.289NX42+0.108NX1NX2+0.504NX1NX4-0.0892NX2NX4+0.0774NX1NX2NX4-0.0625NX13+0.0312NX23-0.25NX43,D1=-0.703-0.161NX2+0.701NX3+0.336NX22+0.408NX32+0.0187NX2NX3+0.219NX23,T2=-0.284+0.859T1+0.404T4-0.407T5+0.062T12-0.0897T42+1.31T52-0.271T1T4+0.383T1T5-0.527T4T5-0.185T1T4T5-0.833T13-0.288T43+2.1T53,T3=0.00215+0.944T2+0.0195D1+0.016NX3+0.0477T22+0.0845D12-0.0514NX32-0.254T2D1+0.0628T2NX3+0.0684D1NX3+0.00345T2D1NX3+0.0198T23-0.0565D13,Ucracking  percentage=Output=32.1+25.3T3.Figure 5Table8 compares the errors predicted by the abductive model with the experimental value. These three experimental cases are excluded from the 18 sets (Table 4) of experimental data for establishing the model. This set of data is used to test the appropriateness of the abductive network model. Table 8 shows that the maximum error of cracking percentage is approximately 8.611 per cent, indicating that the proposed abductive model is suitable for drilling process on alumina ceramic.Table 8
Comparison of cracking percentage between neural network and actual experiment.
Trial number
Spindle speed (rpm)
Feed rate (mm/min)
Depth of cut (mm)
Abrasive grain (#)
Cracking percentage measured by actual experiment
Cracking percentage predicted by neural network
Error %
1
27100
11
0.02
140
17.349
16.760
3.395
2
23000
18
0.06
180
45.236
42.448
6.163
3
20500
15
0.05
100
13.158
12.025
8.611


## 7. Conclusion

A short hole is defined as a hole with a small ratio of depth to diameter. In this study, a 4 mm diameter short hole is perforated on 3.2 mm thickness ceramic, and the crack area obviously appeared when bad drill parameters are selected. This study used the Taguchi method and neural network on alumina (Al2O3) ceramics to optimize drilling conditions with the aim of reducing the crack area at the exit side of alumina. Neural network is used for predicting crack area percentage.The experimental results may be summarized as follows:(1)
ANOVA verified that factorsC, D, and B are very significant for drilling alumina. Factor A had less influence.(2)
The optimal combination of drilling parameters was obtained by Taguchi method on alumina ceramic, including spindle speed (factorA) of 25000 rpm, feed rate (factor B) of 20 mm/min, depth of cut (factor C) of 0.02 mm, and abrasive grain size (factor D) of 100 mesh.(3)
The optimal combination of drilling parameters was confirmed through confirmation experiments. The results indicate that the error between confirmation test and predicted values was within 7.501%.(4)
A comparison was made between the actual drill experiments and the model prediction value of the neural network. The comparison showed that the abductive network model prediction was very close to the actual experimental value, with an error rate of less than 8.611%. The paper demonstrated that engineers can use this abductive model to simulate drilling conditions and predict the crack under any drilling parameters. In addition, the technique can reduce the time needed to experiment.



---

*Source: 304691-2015-09-10.xml*