UDK 669.715:620.178	ISSN 1580-2949
Original scientific article/Izvirni znanstveni članek	MTAEC9, 44(3)129(2010)
USE OF GREY BASED TAGUCHI METHOD IN BALL BURNISHING PROCESS FOR THE OPTIMIZATION OF SURFACE ROUGHNESS AND MICROHARDNESS OF AA 7075 ALUMINUM ALLOY
UPORABA GREY-TAGUCHIJEVE METODE PRI PROCESU GLAJENJA ZA OPTIMIZACIJO POVRŠINSKE HRAPAVOSTI IN MIKROTRDOTE ALUMINIJEVE ZLITINE AA 7075
Ugur Esme
Mersin University Tarsus Technical Education Faculty, Department of Mechanical Education, 33140, Tarsus-Mersin/Turkey
uguresme@gmail.com
Prejem rokopisa - received: 2009-11-09; sprejem za objavo - accepted for publication: 2010-03-20
This study investigated the multi-response optimization of burnishing process for an optimal parametric combination to yield favorable surface roughness and microhardness using the Grey relational analysis and Taguchi method. Sixteen experimental runs based on an orthogonal array of Taguchi method were performed to derive objective functions to be optimized within experimental domain. The objective functions have been selected in relation of burnishing parameters; burnishing force, number of passes, feed rate and burnishing speed. The Taguchi approach followed by Grey relational analysis was applyed to solve the multi-response optimization problem. The significance of the factors on overall quality characteristics of the burnishing process has also been evaluated quantitatively with the variance method (ANOVA). Optimal results were verified through confirmation experiments. This shows application feasibility of the Grey relation analysis in combination with Taguchi technique for continuous improvement in product quality in manufacturing industry.
Keywords: ball burnishing process, Grey relation analysis, Taguchi method
V tej študiji je raziskana večodgovorna optimizacija procesa glajenja z dosego optimalnih kombinacij parametrov za ugodno površinsko hrapavost in mikrotrdoto z uporabo Greyjeve analize odvisnosti in Taguchijeve metode. Šestnajst eksperimentov v ortogonalni porazdelitvi po metodi Taguchi je bilo uporabljenih za razvoj objektivnih funkcij za optimizacijo v eksperimentalnem polju. Objektivne funkcije so bile izbrane v odvisnosti od parametrov glajenja; sila glajenja, število prehodov, hitrost podajanja in hitrost glajenja. Taguchijev približek in Greyjeva analiza odvisnosti sta bila uporabljena za rešitev večodgovornega problema. Kvantitativno je bil ocenjen tudi pomen dejavnikov kakovosti procesa glajenja z metodo variance (ANOVA). Optimalni rezultati so bili potrjeni s preizkusi. Delo dokazuje uporabnost Greyjeve analize odvisnosti in Taguchijeve tehnike za stalno izboljšanje kakovosti proizvodov v predelovalni industriji.
Ključne besede: krogelno glajenje, Greyjeva analiza odvisnosti, Taguchijeva metoda
1 INTRODUCTION
surface layers4. Accordingly, burnishing distinguishes itself from chip-forming finishing processes such as
The function performance of a machined component grinding, honing, lapping and super-finishing which
such as fatigue strength, load bearing capacity, friction,	induce residual tensile stresses at the machined surface
etc. depends to a large extent on the surface as	layers5,6. Also, burnishing is economically desirable, topography, hardness, nature of stress and strain induced because it is a simple and cheap process, requiring less
on the surface region. Nowadays, about 50% of the	time and skill to obtain a high-quality surface finish4,5.
energy supplied is lost in the friction of elements in	Beside producing a good surface finish, the burnish-
relative motion1,2. Roughness values less than 0.1 mm	ing process has additional advantages over other machi-are required for good aesthetic appearance, easy mould ning processes, such as securing increased hardness,
release, good corrosion resistance, and high fatigue	corrosion resistance and fatigue life as a result of
strength. During recent years, however, considerable	producing compressive residual stress. Residual stresses
attention has been paid to the post-machining metal	are probably the most important aspect in assessing
finishing operations such as burnishing which improves	integrity because of their direct influence on perfor-
the surface characteristics by plastic deformation of the	mance in service. Thus, control of the burnishing process
surface layers2,3.	(burnishing conditions) in such a way as to produce
Burnishing is considered as a cold-working finishing	compressive residual stresses in the surface region could
process, differing from other cold-working, surface	lead to considerable improvement in component life. A treatment processes such as shot peening and sand comprehensive classification of burnishing tools and
blasting, etc. in that it produces a good surface finish and	their application has been given by Shneider7. A lite-
also induces residual compressive stresses at the metallic	rature survey shows that work on the burnishing process
has been conducted by many researchers and the process improves also the properties of the parts, e.g. higher wear resistance2,8,9 increased hardness10-12, surface quality2,3,14 and increased maximum residual stress in compression11. The parameters affecting the surface finish are: burnishing force, feed rate, ball material, number of passes, workpiece material, and lubrication2,3. It is necessary to find an optimal process condition capable of producing desired surface quality and hardness. However, this optimization should be performed in such a way that all the objectives should fulfill simultaneously. Such an optimization technique is called multi-response optimization15.
The majority of the research existing in literature on the effect of burnishing parameters on the burnished surface is of experimental nature and very few analytical models are available in the literature.
The Taguchi method is very popular for solving optimization problems in the field of production engineering.16,17 The method utilizes a well-balanced experimental design (allows a limited number of experimental runs) called orthogonal array design, and signal-to-noise ratio (S/N ratio), which serve as objective function to be optimized (maximized) within the experimental domain. However, traditional Taguchi method cannot solve multi-objective optimization problem. To overcome this, the Taguchi method coupled with Grey relational analysis has a wide area of application in manufacturing processes. This approach can solve multi-response optimization problem simultaneously15,18.
Planning the experiments through the Taguchi orthogonal array has been used quite successfully in process optimization19-24. Therefore, in this study the Taguchi L16(44) orthogonal array was applied to plan the experiments on burnishing process.
Four controlling factors including burnishing force (F), number of passes (N), feed rate (f) and burnishing speed (V) on the surface roughness (Ra) and micro-hardness (HV) with four levels for each factor were selected. The Grey relational analysis was then applied to examine how the burnishing parameters influenced the surface roughness and microhardness and an optimal parameter combination was then obtained. Through analyzing the Grey relational grade matrix, the most influential factors for individual quality targets of burnishing process can be identified. Additionally, the analysis of variance (ANOVA) was also utilized to examine the most significant factors for the surface roughness and microhardness in burnishing process.
2 GREY RELATIONAL ANALYSIS
In Grey relational analysis, experimental data i.e.,
measured features of quality characteristics are first normalized ranging from zero to one. This process is known as Grey relational generation. Next, based on normalized experimental data, Grey relational coefficient
is calculated to represent the correlation between the desired and actual experimental data15. Then overall Grey relational grade is determined by averaging the Grey relational coefficient corresponding to selected responses. The overall performance characteristic of the multiple response process depends on the calculated Grey relational grade. This approach converts a multiple response process optimization problem into a single response optimization situation with the objective function which is the overall Grey relational grade. The optimal parametric combination is then evaluated which would result in the highest Grey relational grade. The optimal factor setting for maximizing overall Grey relational grade can be obtained by Taguchi method15.
In Grey relational generation, the normalized Ra values corresponding to the smaller-the-better (SB) criterion which can be expressed as:
(k) =
max y. (k) - y. (k) max y^ (k) - min y^ (k)
(1)
HV100 should follow the larger-the-better (LB) criterion, which can be expressed as:
y^ (k) - min y^ (k)
X, (k) =
max y. (k) - min y. (k)
(2)
where x,(k) is the value after the Grey relational generation, min yi(k) is the smallest value of y,(k) for the klh response, and max yi(k) is the largest value of y,(k) for the klh response15. An ideal sequence is [xo(k) (k=1,
2, 3......, 16)] for the responses. The definition of Grey
relational grade in the course of Grey relational analysis is to reveal the degree of relation between the 16
sequences [xo(k) and x,(k), i=1, 2, 3,.......,16]. The Grey
relational coefficient ^,(k) can be calculated as:
^i(k)=
A min - „
A 0 i (k)+^A „
(3)
where =
Xo (k)-Xi (k)
the absolute value of the dif-
ference of Xo(k) and Xi(k); ^ is the distinguishing coefficient 0 <xp< 1; Ami„ =V;mi" e iVkmi"||x0(k)-Xj (k) the smallest value of A0i; and
Amax =V;ma" e iVk|x0(k)-x^. (k)\\ is the largest value
of A0i. After averaging the Grey relational coefficients, the Grey relational grade yi can be computed as:
1 n
y i = -1 ^ i (k)
n k=1
(4)
where n is the number of process responses. The higher value of Grey relational grade corresponds to intense relational degree between the reference sequence X0(k) and the given sequence x,■(k). The reference sequence X0(k) represents the best process sequence. Therefore, higher Grey relational grade means that the corresponding parameter combination is closer to the optimal15. The mean response for the Grey relational grade with its
grand mean and the main effect plot of Grey relational grade are very important because optimal process condition can be evaluated from this plot15.
3 EXPERIMENTAL DETAILS AND TEST RESULTS
3.1 Workpiece Material
In this study, high strength precipitation hardening 7XXX series wrought aluminum alloy AA 7075 was used. The strength and good mechanical properties make the AA 7075 aluminum alloy appropriate for the use in aerospace industry. The chemical composition and mechanical properties of the workpiece material is given in Table 1.
Table 1: Chemical composition and mechanical properties of AA7075 aluminum alloy Tabela 1: Kemi~na s AA 7075
Chemical composition (%)
Mechanical properties
Al
90.0
Tensile strength (MPa)
220
Parameters	Notation	Unit	Levels of factors			
			1	2	3	4
Burnishing force	F	N	58 . 86*	117.72	176.58	235.44
Number of passes	N	-	1*	2	3	4
Feed rate	f	mm/ min	0.1*	0.2	0.3	0.4
Burnishing speed	V	rpm	200*	300	500	700
Yield strength (MPa)
9.5
Shear strength (MPa)
150
Fatique strength (MPa)
160
Zn
Hardness (HV100)
150
The workpiece material, as shown in Figure 1, was prepared with the diameter of 30 mm and 70 mm in length as a three part each having 20 mm length.
Figure 1: Dimensions of workpiece material25 Slika 1: Mere preizkusanca
3.2	Machines and Equipments
An 18 mm diameter ball was used for burnishing. The detailed and drawing is shown in Figure 2. When the ball or roller is pressed against the surface of the metallic specimen, a pre-calibrated spring was compressed. This spring is being used mainly to reduce the possible sticking of the tool onto the surface25.
The experiments were performed on a FANUC GT-250B CNC lathe. The burnishing tool was mounted on the CNC turret. Dry turning and burnishing were used in all the experimental work, but alcohol was used to clean the specimens before burnishing. Cleaning of the ball was carried out continuously in order to prevent hard particles from entering on the contact surface between the tool and the specimen, such hard particles usually leaving deep scratches, which may damage the burnished surface of the specimen25.
Phynix TR-100 model surface roughness tester was used to measure the surface roughness of the burnished samples. Cut off length was chosen as 0.3 for each roughness measurement. Vickers microhardness tester with 100 g load (HV100) was used for microhardness measurements25. Six measurements of surface roughness and microhardness were taken from the samples and average of the values were used in the multi-criteria optimization.
3.3	Process Parameters and Test Results
In full factorial design, the number of experimental runs exponentially increases as the number of factors as well as their level increases. This requires huge experimentation cost and considerable time. So, in order to compromise these two adverse factors and to search the optimal process condition through a limited number of experimental runs Taguchi's L16(44) orthogonal array consisting of 16 sets of data was selected to optimize the multiple performance characteristics of surface roughness. The burnishing parameters used in this study are shown in Table 2.
Table 2: Process parameters and their limits Tabela 2: Parametri in omejitve procesa
Parameters
Burnishing force
Number of
passes
FpprI rntp
Burnishing speed
Notation
N
V
Unit
N
mm/
min
rpm
Levels of factors
58 86'
0 1
900'
117 72
0.2
300
3
176.58
0.3
500
235 44
0 4
700
Figure 2: Detailed drawing of the ball burnishing tool: (1) casing; (2) adapter cover; (3) spring25
Slika 2: Na~rt gladilnega orodja: (1) ohišje; (2) prilagoditveni pokrov, (3) vzmet25
'Initial factor settings
Table 3 shows the selected design matrix based on Taguchi L16(44) orthogonal array consisting of 16 sets of coded conditions and the experimental results for the
P
i
responses of Ra and HV100. All these data were utilized for the analysis and evaluation of optimal parameter combination required to achieve desired surface quality within the experimental domain.
Table 3: Orthogonal array L16(44) of the experimental runs and results Tabela 3: Ortogonalna porazdelitev L16(44) preizkusov in rezultati
Run no	Process parameters				Experimental results	
	F	N	f	V	Ra (^m)	HV100
1		1	1	1	0.61	160
2		2	2	2	0.57	167
3		3	3	3	0.53	173
4		4	4	4	0.59	187
5		1	2	3	0.36	165
6		2	1	4	0.18	178
7		3	4	1	0.25	181
8	2	4	3	2	0.20	193
9	3	1	3	4	0.20	172
10	3	2	4	3	0.22	187
11	3	3	1	2	0.08	196
12	3	4	2	1	0.23	210
13	4	1	4	2	0.30	191
14	4	2	3	1	0.15	199
15	4	3	2	4	0.14	205
16	4	4	1	3	0.19	212
4 PARAMETRIC OPTIMIZATION OF BURNISHING PROCESS
4.1 Evaluation of Optimal Process Condition
First, by using Eqs. (1) and (2), experimental data were normalized to obtain Grey relational generation15. The normalized data and A0i for each of the responses are given in Table 4 and Table 5 respectively. For Ra smaller-the-better (SB) and for HV larger-the-better (LB) criterion has been selected.
Table 4: Grey relational generation of each performance characteristics
Tabela 4: Generacija Greyjeve odvisnosti za karakteristike vsake performance
Run no	Ra	HV
	Smaller-the-better	Larger-the-better
Ideal sequence	1	1
1	0.000	0.000
2	0.043	0.075
3	0.250	0.151
4	0.519	0.038
5	0.096	0.472
6	0.346	0.811
7	0.404	0.679
8	0.635	0.774
9	0.231	0.774
10	0.519	0.736
11	0.692	1.000
12	0.365	0.717
13	0.596	0.585
14	0.750	0.868
15	0.865	0.887
16	1.000	0.792
Table 5: Evaluation of A0i for each of the responses		
Tabela 5: Ocena A0i za vsak odgovor		
Run no	Ra	HV
Ideal sequence	1	1
1	1.000	1.000
2	0.957	0.925
3	0.750	0.849
4	0.481	0.962
5	0.904	0.528
6	0.654	0.189
7	0.596	0.321
8	0.365	0.226
9	0.769	0.226
10	0.481	0.264
11	0.308	0.000
12	0.635	0.283
13	0.404	0.415
14	0.250	0.132
15	0.135	0.113
16	0.000	0.208
Table 6 shows the calculated Grey relational coefficients (with ^Ra = 0.67, ^hv = 0.33) of each performance characteristic using Eq. (3).
Table 6: Grey relational coefficient of each performance characteristics (^Ra = 0.67, ^HV = 0.33)
Tabela 6: Greyjev odvisnostni koeficient za karakteristike vsake performance (ipRa = 0.67, ■^HV = 0.33)
Run no	Ra	HV
Ideal sequence	1	1
1	0.333	0.333
2	0.343	0.351
3	0.400	0.371
4	0.510	0.342
5	0.356	0.486
6	0.433	0.726
7	0.456	0.609
8	0.578	0.688
9	0.394	0.688
10	0.510	0.654
11	0.619	1.000
12	0.441	0.639
13	0.553	0.546
14	0.667	0.791
15	0.788	0.815
16	1.000	0.707
The Grey relational coefficients, given in Table 7, for each response was accumulated by using Eq. (4) to evaluate Grey relational grade, which is the overall representative of all the features of burnishing quality. Thus, the multi-criteria optimization problem was transformed into a single equivalent objective function optimization problem using the combination of Taguchi
Run no	Grey relational grade	Rank
1	0.333	16
2	0.348	15
3	0.380	14
4	0.397	13
5	0.443	12
6	0.629	6
7	0.558	10
8	0.651	5
9	0.591	8
10	0.606	7
11	0.874	1
12	0.573	9
13	0.548	11
14	0.750	4
15	0.806	2
16	0.803	3
10
Run no	S/N
1	-9.54
2	-9.16
3	-8.40
4	-8.02
5	-7.07
6	-4.02
7	-5.06
8	-3.72
9	-4.57
10	-4.34
11	-1.17
12	-4.83
13	-5.21
Grey relational grade
0.333
0.348
0.380
0.397
0.443
Factors	Grey relational grade				
	Level 1	Level 2	Level 3	Level 4	Max-Min
F	0.36	0.57	0.66	0.73	0.37
N	0.48	0.58	0.65	0.61	0.17
f	0.66	0.54	0.59	0.53	0.13
V	0.55	0.61	0.56	0.61	0.06
0.606
0.874
0.573
0.548
0.750
0.806
0.803
Rank
16
15
14
13
12
11
a function would be offered (efficiency evaluation) about one design factor as the standard to evaluate the efficiency and to understand the experiment efficiency. While Taguchi experimental design took the quality loss as the base, it would design one statistic to evaluate efficiency, which is called the S/N ratio. In this section we had one quality characteristic that is grey relational grade for Taguchi analysis26. Table 8 shows the S/N ratio based on the larger-the-better criterion for overall Grey relational grade calculated by using Eq. (5).
r 1 11
s / N=-10lg-y
n i=1 v:
(5)
where n is the number of measurements, and yi is the measured characteristic value.
Table 8: S/N ratio for overall Grey relational grade
Tabela 8: Razmerie S/N za splošno Greyievo odvisnostno stopnio
Run no
1
10
11
12
13
S/N
-9.54
-9.16
-8.40
-8 02
-7.07
-4.02
-5.06
-3.72
-4.57
-4.34
-1.17
-4.83
-5.21
F1 F2 F3 F4 N1 N2 N3 N4 f1 f2 f3 f4 VI V2 V3 V4
Parameter level
Figure 3: S/N ratio plot for the overall Grey relational grade Slika 3: S/N razmerje za splošno Greyjevo stopnjo
14
15
16
-2.50
-1.87
-1.90
Graphical representation of S/N ratio for overall Grey relational grade is shown in Figure 3. The dashed line is the value of the total mean of the S/N ratio.
As indicated in Figure 3, the optimal condition for the burnishing of AA7075 aluminum alloy becomes F4N3fiV4. Table 9 shows the mean Grey relational grade ratio for each level of the process parameters.
Table 9: Response table for the mean Grey relational grade
Tabela 9: Tabela odgovorov za povprečno Greyjevo odvisnostno
stopnjo
Factors
N
V
Level 1 Level 2 Level 3 Level 4 Max-Min
0.36
0.48
0.66
0.55
Grey relational grade
0.57
0.58
0.54
0.61
0.66
0.65
0.59
0.56
0.73
0.61
0.53
0.61
0.37
0.17
0.13
0.06
Total mean Grey relational grade= 0.58
4.2 Analysis of Variance (ANOVA)
The purpose of the analysis of variance (ANOVA) is to investigate which burnishing parameters significantly affect the performance characteristices. This is accomplished by separating the total variability of the grey relational grades, which is measured by the sum of the squared deviations from the total mean of the grey relational grade, into contributions by each burnishing parameters and the error. Thus;
SS T = SS F + SS e	(6)
where
SS
^ =1t (y. -
j=1
T. )
(7)
1
2
4
5
6
7
9
7
1
9
4
F
f
2
a
6
1
Q
and
SSt	Total sum of squared deviations about the mean
yj	Mean response for experiment
Ym	Grand mean of the response
p	Number of experiments in the orthogonal array
SSf	Sum of squared deviations due to each factor
SSe	Sum of squared deviations due to error
In addition, the F test was used to determine which burnishing parameters have a significant effect on the performance characteristic. Usually, the change of the burnishing parameter has a significant effect on the performance characteristics when the F value is large. ANOVA for overall Grey relational grade is shown in Table 10.
Table 10: ANOVA results Tabela 10: ANOVA rezultati
Parameter
N
V
Error
Total
Degree of Freedom
15
Sum of Square
0.300
0.066
0.043
0.009
0.001
0.419
Mean Square
0.100
0.022
0.014
0.003
0.005
19.00
4.20
2.72
0.63
Contribution (%)
71.59
15.75
10.26
2.14
0.23
100
According to this analysis, the most effective parameters with respect to surface roughness and microhardness are burnishing force, number of passes, burnishing feed and burnishing speed. Percent contribution indicates the relative power of a factor to reduce variation. For a factor with a high percent contribution, has a great influence on the performance. The percent contributions of the burnishing parameters on the surface roughness and microhardness are shown in Table 10 and Figure 4. Burnishing force (71.59%) was found to be the major factor affecting surface roughness and microhardness, whereas number of passes (15.75%), burnishing feed (10.26%) and burnishing speed (2.14%) were found to be the second, third and fourth ranking factor respectively.
Figure 4: Contribution percentage of the burnishing parameters Slika 4: Odstotna porazdelitev parametrov glajenja
4.3 Confirmation Test
After evaluating the optimal parameter settings, the next step is to predict and verify the enhancement of quality characteristics using the optimal parametric combination. The estimated Grey relational grade yY using the optimal level of the design parameters can be calculated as:
y/ = y
+£ (Y j -
i=1
■y m )
(8)
where ym is the total mean Grey relational grade, yi is the mean Grey relational grade at the optimal level, and o is the number of the main design parameters that affect the quality characteristics. Table 11 indicates the comparison of the predicted surface roughness and microhardness with that of actual by using the optimal burnishing conditions; good agreement between the actual and predicted results was obtained. Also, improvement in overall Grey relational grade was found to be as 0.65.
Table 11: Results of confirmation test Tabela 11: Rezultati potrditvenih preizkusov
Factor levels
Ra
HV100
S/N ratio of overall Grey relational grade
Overall Grey relational grade
Initial factor settings
F1N1f1V1
0.61
160
-9.54
0.33
Optimal process condition
Prediction
F4N3f1V4
-0.51
0.91
Experiment
F4N3f1V4
0.12
200
-0.80
0.98
Improvement in Grey relational grade=0.65
In Taguchi method, the only performance feature is the overall Grey relational grade and the aim should be to search a parameter setting that can achieve highest overall Grey relational grade. The Grey relational grade is the representative of all individual performance characteristics. In this study, objective functions have been selected in relation to parameters of surface roughness and microhardness. The weight calculations were done by using Analytic Hierarchy Process (AHP) and the weights were found to be as 0.67 and 0.33 for the responses of surface roughness and microhardness respectively. The results showed that using optimal parameter setting (F4N3f1V4) caused lower microhardness together with higher microhardness.
5 CONCLUSIONS
Taguchi method is a very effective tool for process optimization under limited number of experimental runs. Essential requirements for all types burnishing processes are smoother surface with higher surface hardness. This study has concentrated on the application of Taguchi method coupled with Grey relation analysis for solving
F
F
3
f
3
multi criteria optimization problem in the field of burnishing process. Experimental results have shown that surface roughness and microhardness of burnished aluminum alloy are greatly improved by using Grey based Taguchi method.
6 REFERENCES
1L. Kukielka, Journal of Mechanical Technology, 19 (1989), 319-356
2M. H. El-Axir, International Journal of Machine Tools & Manufacture, 40 (2000), 1603-1617
3	N. H. Loh, S. C. Tam, S. Miyazawa, Journal of Mechanical Working Technology, 18 (1989), 53-61
4	A. M. Hassan, Journal of Materials Processing Technology, 72 (1997), 385-391
5	C. Wick, R. F. Veilleux, Soc. Manuf. Eng., 3 (1985), 16-38
6	T. Siva Prasad, B. Kotiveerachari, J. Inst. Eng., 69 (1988), 55-58
7	G. Y. Shneider, Machines and Tooling, 1 (1989), 34-35
8	A. N. Niberg, Soviet Engineering Research, 7 (1987), 67
9	P. C. Michael, N. Saka, E. Rabinowicz, Wear, 132 (1989), 265
10	N. H. Loh, S. C. Tam, S. Miyazawa, Wear, 129 (1989), 235
11	I. Yashcheritsyn, E. I. Pyatosin, V. V. Votchuga, Soviet Journal of Friction and Wear, 8 (1987), 87
12	M. Fattouh, M. H. El-Axir, S. M. Serage, Wear, 100 (1988), 127
13	A. S. H. Al-Bsharat, Investigation into the effect of burnishing process on some properties of cylindrical metallic materials, Jordan University of Science and Technology, MSc Thesis, (1994)
14	O. M. Abd Al-Wahhab, The effects of roller-burnishing on some properties of non-ferrous metals, Jordan University of Science and Technology, MSc Thesis, (1996)
15	S. Datta, A. Bandyopadhyay, P. K. Pal, International Journal of Advanced Manufacturing Technology, 39 (2008), 1136-1143
16	W. H., Yang, Y. S. Tarng, Journal of Materials Processing Technology, 84 (1998), 122-129
17	H. Rowlands, J. Antony, The TQM Magazine, 12 (2000), 78-83
18	S. H. Lim, C. M. Lee, W. J. Chung, International Journal of Precision Engineering and Manufacturing, 7 (2006), 18-23
19	D. C. Chen, C. F. Chen, Journal of Materials Processing Technology, 190 (2007), 130-137
20	C. P. Fung, P. C. Kang, Journal of Materials Processing Technology, 170 (2005), 602-610
21	S. H. Tang, V. J. Tan, S. M. Sapuan, S. Sulaiman, N. Ismail, R. Samin, Journal of Materials Processing Technology, 182 (2007), 418-426
22	P. Vijian, V. P. Arunachalam, Journal of Materials Processing Technology, 180 (2006), 161-166
23	L. J. Yang, Journal of Materials Processing Technology, 185 (2007), 113-119
24	J. Z. Zhang, J. C. Chen, E. D. Kirby, Journal of Materials Processing Technology, 184 (2007), 233-239
25	U. Esme, A. Sagbas, F. Kahraman, M. K. Kulekci, Materials and Technology 42 (2008), 215-219
26F. C. Kung, Y. J. Weng, Polymer-Plastics Technology and Engineering, 47 (2008), 1154-1161