Venugopal2003.pdf Selection of Optimum Conditions For Maximum Material Removal
Venugopal2003.pdf Selection of Optimum Conditions For Maximum Material Removal
Received 4 March 2003; received in revised form 3 June 2003; accepted 10 June 2003
Abstract
Efficient grinding of structural ceramics requires judicious selection of operating parameters to maximize removal rate while
controlling surface integrity. Grinding of silicon carbide is difficult because of its low fracture toughness, making it very sensitive
to cracking. In the present work, experiments were carried out to study the effect of wheel parameters; grain size and grain density
and grinding parameters; depth of cut and feed on the surface roughness and surface damage. The significance of the grinding
parameters on the selected responses was evaluated using analysis of variance. Mathematical models were developed using the
experimental data considering only the significant parameters. A genetic algorithm (GA) code has been developed to optimize the
grinding conditions for maximum material removal, using a multi-objective function model, by imposing surface roughness and
surface damage constraints. The choice of including manufacturer’s constraints on the basis of functional requirements of the
component for maximizing the production rate was also embedded in the GA code.
2003 Elsevier Ltd. All rights reserved.
1. Introduction which are visible, and the later damage is due to median
and lateral cracks that are formed below the affected
There has been an increased interest in the use of grinding zone which are not visible [2].
advanced ceramic materials in the recent past due to their Low thermal coefficient of expansion, low density and
unique physical and mechanical properties. The advan- relatively high thermal conductivity are the special features
tages of ceramics over other materials include high hard- of silicon carbide ceramics. In view of these properties,
ness and strength at elevated temperatures, chemical SiC is expected to be used increasingly for heat resistant
stability, attractive high temperature wear resistance and parts [3]. Grinding is often the method of choice for mach-
low density [1]. Structural ceramics such as silicon ining ceramics in large-scale production and automation.
nitride and silicon carbide are now being increasingly Despite various research efforts in ceramic grinding over
used in bearings, valves, rotors and other applications the past two decades, much needs to be established to stan-
where a close dimensional tolerance is required. dardize models for process optimization for improving
Ceramic materials possess very low fracture tough- product quality and increasing productivity to reduce
ness compared with metals and alloys. This means that machining cost. The effective use of ceramics in industrial
they are very sensitive to the forces introduced due to applications demands the machining of ceramics with
machining. The major form of damage usually occurs as good surface finish and low surface damage. Therefore, in
surface and subsurface damage. The first type of damage the present work, an attempt has been made to examine
is due to radial cracks formed on the ground surface the effect of various process parameters and wheel proper-
ties on these responses during grinding of silicon carbide
with diamond abrasives. A genetic algorithm (GA) based
∗
Corresponding author. Tel.: +91-11-659-1443; fax: +91-11-868- optimization procedure has been developed to optimize
2037. grinding conditions, viz. depth of cut, work speed, grit size
E-mail address: [email protected] (P. Venkateswara Rao). and density, using a multi-objective function model.
0890-6955/$ - see front matter 2003 Elsevier Ltd. All rights reserved.
doi:10.1016/S0890-6955(03)00165-2
1328 A. Venu Gopal, P. Venkateswara Rao / International Journal of Machine Tools & Manufacture 43 (2003) 1327–1336
Nomenclature
C constant in mathematical model
C1 constant in surface roughness model
C2 constant in percentage area of surface damage model
d depth of cut, µm
f table feed rate, m/min
M grit size (mesh)
R grain density
Y machining response
%D percentage area of surface damage
Ra surface roughness measured, µm
d.o.f. degrees of freedom
MRR material removal rate, mm3/mm width/min
a depth of cut exponent in mathematical model
a1 depth of cut exponent in surface roughness model
a2 depth of cut exponent in percentage area of surface damage model
b feed rate exponent in mathematical model
b1 feed rate exponent in surface roughness model
b2 feed rate exponent in percentage area of surface damage model
g grit size exponent in mathematical model
g1 grit size exponent in surface roughness model
g2 grit size exponent in percentage area of surface damage model
h grain density exponent in mathematical model
h1 grain density exponent in surface roughness model
h2 grain density exponent in percentage area of surface damage model
lyzed the mechanism of material removal in ceramic roughness prediction model for turning mild steel using
grinding with the help of indentation fracture mechanics response surface methodology to give the factor effects
approach and the machining approach. It was shown of the individual process parameters. They have also
with the first approach that median/radial cracks are usu- attempted to optimize the surface roughness as an objec-
ally associated with strength degradation, and lateral tive function using GA.
cracks with material removal. Xu et al. [8] have demon- Most of the research has been conducted with silicon
strated that the mechanism of material removal and the nitride as workpiece material to investigate material
effect of machining induced damage on strength of removal rate, effect of residual stresses on fracture
advanced ceramics can be controlled by approximately strength and surface integrity in grinding. This paper
tailoring the microstructure. This not only promotes easy focuses on the effect of grit size, grit density, depth of
and well controlled material removal by grain dislodge- cut and work feed on surface finish and damage pro-
ment during machining, but also suppresses the forma- duced during grinding of silicon carbide, as this is one
tion of strength degrading cracks. of the materials being used increasingly in structural
Detailed knowledge on the effect of the grinding pro- applications. The results of the experimental work are
cess on surface integrity gives the opportunity for a bet- used to find the significance of the grinding variables
ter exploitation of ceramic materials by improved pro- on the machining responses by developing mathematical
cess conditions. Pfeiffer and Hollstein [9] have used the models. The mathematical models thus developed are
X-ray diffraction technique for determining the damage further utilized to find the optimum grinding parameters
induced in ground silicon nitride and alumina and ther- using genetic algorithms employing a multi-objective
eby established correlations between micro-plastic defor- function model.
mation and amount of damage. These investigations
show that in the case of lapped and ground alumina and
of ground silicon nitride, bending strength is dominated 3. Methodology
by machining induced damage. In the case of lapped sili-
con nitride, the effect of damage is compensated by To obtain applicable and practical predictive quanti-
machining induced compressive residual stresses. Dani- tative relationships, it is necessary to model the grinding
els [10] has investigated the influence of surface grinding responses and the process variables. These models would
parameters such as diamond abrasive type, wheel speed be of great use during optimization of the ceramic grind-
and down feed on the rupture strength of silicon carbide. ing process. The experimental results were used to model
It was found that more severe grinding conditions, such the various responses using the multiple regression
as higher normal forces and power consumption, did not method by using a non-linear fit among the response and
significantly reduce the mean rupture strength of the the corresponding significant parameters. The signifi-
material. The most encouraging aspect inferred from cance of the parameters on the responses could be found
these results was that grinding conditions could be using ANOVA. In this work, a commercially available
changed in order to optimize the process without sig- mathematical software package MATLAB was used for
nificant structural damage to the work material. the computation of regression constants and parameters.
Models contribute significantly to the comprehension The mathematical modeling of surface roughness and
of the process itself, and form the basis for the simul- surface damage in the grinding of ceramic material
ation of the grinding processes. They thus create a pre- involved lots of other factors, such as wheel speed, type
condition for increased efficiency while ensuring a high of abrasive grit, grit density, etc. However, to facilitate
product quality at the same time. Tonshoff et al. [11] experimental data collection, only four dominant factors
have described the state of the art in the modeling and were considered in the planning of the experimentation.
simulation of grinding processes, comparing the differ- A high precision hydraulic surface grinding machine was
ent approaches to modeling. Furthermore, the benefits as used to grind silicon carbide work material with diamond
well as the limits of model application and simulation grinding wheels at a speed of 36 m/s. The silicon carbide
were discussed in their paper. Liao and Chen [12] have material has been made from 3 µm particle sized powder
used back-propagation neural networks for modeling and at a sintering temperature of 2100 °C, which gives a
optimizing the creep-feed grinding of alumina with dia- density of 3.17 g/cm3 after processing. The factors con-
mond wheels. The authors have shown that it is possible sidered were depth of cut, feed rate, grit size, and grit
to obtain a global optimal solution by the proper use of density. The experiments were planned using a complete
the Boltzman factor. Jain and Jain [13] have also used 34 factorial design. Based on this, a total of 81 experi-
a similar approach for optimum selection of machining ments, each having a combination of different levels of
conditions in abrasive flow machining using neural net- factors as shown in Table 1, were carried out.
works. The results have been validated by comparing the The responses measured were tangential grinding
optimized machining conditions obtained using genetic force, surface roughness and percentage area of surface
algorithms. Suresh et al. [14] have developed a surface damaged. The grinding force was measured with an
1330 A. Venu Gopal, P. Venkateswara Rao / International Journal of Machine Tools & Manufacture 43 (2003) 1327–1336
Fig.(a) d:5µm; f:15m/min; M:500; R:100 and %D:1.38, (b) d:5µm; f:5m/min; M:120; R:50 and %D:3.07, (c) d:15µm; f:15m/min; M:120; R:75 and %D:3.21,
(d) d:5µm; f:15m/min; M:240; R:50 and %D:1.61, (e) d:30µm; f:15m/min; M:500; R:100 and %D:2.34, (f) d:30µm; f:15m/min; M:120; R:75 and %D:3.92
(magnification is same for all the figures shown in Fig.(a))
nation could not be carried out due to the limitation of Due to limitation of resources, testing could be carried
experimental resources. Accordingly, it has been out on only 27 specimens. Since all the factors had three
assumed that the four factor interaction is not present levels, an L27 orthogonal array (OA) has been used to
and the corresponding sum of squares and degrees of select the specimens for testing. The results of percent-
freedom have been taken as residual/error to conduct age area of surface damage are shown in Table 3. How-
the ANOVA. ever, since the experiments involved four factors (d, f,
The results of ANOVA for surface roughness, shown M and R) only, these factors are assigned to columns 1,
in Table 2, indicate that the variables d, f, M and R are 5, 8 and 11, respectively, of L27 OA. This assignment
significant at the 99% confidence level and their interac- makes it possible to test the significance of not only the
tions are not significant even at the 95% confidence factors but also the interactions between “d and f”, “d
level. Thus, the parameters d, f, M and R are considered and M” and “d and R” for the levels of factors con-
while developing the model for roughness, and it is sidered.
given by: The results of ANOVA for the percentage area of sur-
face damage, shown in Table 4, indicate that the vari-
Ra ⫽ 0.86d0.1843f 0.5253M⫺0.2866R⫺0.2444 (6)
ables d, f and M are only significant at the 99% confi-
dence level, whereas the parameter R and the interactions
4.2. Surface damage model
are not significant even at the 95% confidence level.
Thus, only the parameters d, f and M have been con-
The measurement of machining induced damage in
sidered in developing the percentage area of surface
ceramics is of interest for two reasons. First, the damage
damage model, and it is given by:
induced will influence the service behavior of the mech-
anical parts. The effect of damage on strength is parti- %D ⫽ 24.44d0.2857f⫺0.3M⫺0.4140 (7)
cularly important in brittle materials, like ceramics.
Second, the damage induced is a consequence of forces It can be observed from the above models that the
occurring during machining. In order to understand the surface roughness decreases with a decrease in depth of
influence of various factors on the induced grinding cut and feed and with an increase in grit number as well
damage, the development of analytical or empirical as grit density. It is also observed that decrease in depth
models for reliable prediction of machining performance of cut or an increase in feed or grit number would
becomes another key issue. decrease the percentage area of surface damage. The
A. Venu Gopal, P. Venkateswara Rao / International Journal of Machine Tools & Manufacture 43 (2003) 1327–1336 1333
Table 2
ANOVA for surface roughness
Total 80 0.9469
Table 4
ANOVA for percentage area of surface damage
Raⱕ(Ra)max,
%Dⱕ(%D)max and
xliⱕxiⱕxui
MRR ⫽ fd
Ra ⫽ 0.86(d)0.1843(f)0.5253(M)⫺0.2866(R)⫺0.2444
%D ⫽ 24.44(d)0.2857(f)⫺0.3(M)⫺0.4140
5 µmⱕdⱕ30 µm
5 m / minⱕfⱕ15 m / min
120ⱕDⱕ500
50ⱕRⱕ100
where xli and xui are the lower and upper bounds on pro-
cess variables xi, and (Ra)max and (%D)max are the
maximum allowable values of surface roughness and
number of flaws.
Fig. 4. General procedure of GA. The GA program had been written in the MATLAB
environment. The maximizing function was written faci-
litating the user to set the constraints. The optimization
process can be described as maximizing the MRR sub- was carried out for different values of constraints on sur-
ject to a criteria set of constraints on surface roughness, face roughness (0.16–0.4 µm) and percentage area of
percentage area of surface damage and input variables. surface damage (1.6–4%). Some typical results, obtained
In order to optimize the present problem using GA: by GA for different surface roughness values imposing
a constraint on percentage area of surface damage as 2,
(a) the following parameters were specified by practice, are presented in Table 5. In a similar way, another set
to get optimal solutions with low computational of typical results, obtained by GA for different percent-
effort: age area of surface damage values placing a constraint
on surface roughness as 0.25 µm, are presented in
앫 maximum number of generations = 300, Table 6.
앫 total string length = 30, It can be observed from the results of optimization
앫 mutation probability = 0.1 and that the MRR increases by almost 10 times by increasing
앫 crossover probability = 0.03; the roughness constraint from 0.16 to 0.4 µm at 2% area
(b) the constrained optimization problem is stated as fol- of surface damage, whereas it increases only 2.5 times
lows: by increasing the damage constraint from 1.6% to 3.0%
at a roughness value of 0.25 µm and remains constant
앫 maximize MRR subject to thereafter. This shows that no improvement in MRR can
A. Venu Gopal, P. Venkateswara Rao / International Journal of Machine Tools & Manufacture 43 (2003) 1327–1336 1335
Table 5
Optimization results of GA at various values of Ra and at 2% area of surface damage
Table 6
Optimization results of GA at various percentage area of surface damages and at Ra = 0.25 µm
also obtained for maximization of material removal formation mechanisms in grinding silicon nitride, Journal of
using surface roughness and percentage damage, com- Materials Research 11 (7) (1996) 1717–1724.
[5] J.E. Mayer Jr., G.P. Fang, Effect of grinding parameters on sur-
puted by the models developed, as constraints. The face finish of ground ceramics, Annals of the CIRP 44 (1) (1995)
results of optimization conclude that the material 279–282.
removal and the cost of grinding are influenced more [6] R.L. Allor, T.J. Whalen, J.R. Baer, K.V. Kumar, Machining of
by the constraint on surface roughness than by surface silicon nitride: experimental determination of process/property
damage. The optimized model developed offers a sol- relationships, in: Proceedings of International Conference on
Machining Advanced Materials, NIST SP847, 1993, pp. 223–234.
ution to reduce the cost of machining, thereby making [7] S. Malkin, T.W. Hwang, Grinding mechanisms for ceramics,
silicon carbide a more commercially viable material for Annals of the CIRP 45 (2) (1996) 569–580.
industrial applications. [8] H.H.K. Xu, N.P. Padutre, S. Jahanmir, Effect of microstructure
on material removal mechanisms and damage tolerance in abras-
ive machining of silicon carbide, Journal of the American Cer-
amic Society 78 (9) (1995) 2443–2448.
Acknowledgements [9] W. Pfeiffer, T. Hollstein, Damage determination and strength pre-
diction of machined ceramics by X-ray diffraction techniques, in:
Proceedings of Conference on Experimental Mechanics, Copen-
The authors would like to express their sincere thanks hagen, Denmark, 1990, pp. 235–245.
to the Department of Science and Technology, New [10] W.H. Daniels, Superabrasives for ceramic grinding and finishing,
Delhi, India, for financially supporting this research pro- SME Technical Paper EM89-125, 1989.
ject and also to the International Advanced Research [11] H.K. Tonshoff, J. Peters, I. Inasaki, T. Paul, Modelling and simul-
Center for Powder Metallurgy and New Materials, Hyd- ation of grinding processes, Annals of the CIRP 41 (2) (1992)
677–688.
erabad, India, for supplying the work material. [12] T.W. Liao, L.J. Chen, A neural network approach for grinding
processes: modelling and optimization, International Journal of
Machine Tools and Manufacture 34 (7) (1994) 919–937.
[13] R.K. Jain, V.K. Jain, Optimum selection of machining conditions
References in abrasive flow using neural networks, Journal of Materials Pro-
cessing Technology 108 (2000) 62–67.
[1] B.P. Bandyopadhyay, The effects of grinding parameters on the [14] P.V.S. Suresh, P. Venkateswara Rao, S.G. Deshmukh, A genetic
strength and surface finish of two silicon nitride ceramics, Journal algorithmic approach for optimization of surface roughness pre-
of Materials Processing Technology 53 (1995) 533–543. diction model, International Journal of Machine Tools and Manu-
[2] T.M.A. Maksoud, A.A. Maksoud, A.A. Mokbel, J.E. Morgan, facture 42 (2002) 675–680.
Evaluation of surface and sub-surface cracks of ground ceramic, [15] E.J.A. Armarego, R.H. Brown, in: The Machining of Metals,
Journal of Materials Processing and Technology 88 (1999) Prentice Hall, New Jersey, 1984, pp. 254–291.
222–243. [16] D.E. Goldberg, Genetic Algorithms, Addison Wesley Longman
[3] I. Inasaki, Grinding of hard and brittle materials, Annals of the Inc, India, 1999.
CIRP 36 (2) (1987) 463–471. [17] D.C. Montgomery, Design and Analysis of Experiments, John
[4] H.H.K. Xu, S. Jahanmir, L.K. Ives, Material removal and damage Wiley and Sons (Asia) Pvt. Ltd, Singapore, 2001.