Experimental and Numerical Investigation of Al Properties Fabricated by CGP Process
Experimental and Numerical Investigation of Al Properties Fabricated by CGP Process
net/publication/257757672
CITATION READS
1 175
3 authors:
F. Djavanroodi
Prince Mohammad University
107 PUBLICATIONS 2,430 CITATIONS
SEE PROFILE
All content following this page was uploaded by F. Djavanroodi on 30 March 2019.
a r t i c l e i n f o a b s t r a c t
Article history: Constrained groove pressing (CGP) process is one of the favorable severe plastic deformation methods to
Received 24 January 2012 fabricate ultra-fine grain sheet materials. In this paper, commercial pure aluminum sheet was subjected to
Received in revised form 31 March 2012 CGP process using conventional and covered sheet casing (CSC) methods. Three dies with different groove
Accepted 30 April 2012
lengths and angles (t = 1 mm and = 45◦ ; t = 1.8 mm and = 45◦ ; t = 1 mm and = 53◦ ) were designed
Available online 15 May 2012
and manufactured. Mechanical properties including strength, elongation to failure and hardness are
obtained and discussed. The influences of these two methods and die groove angle on specimen strain
Keywords:
behavior and required force value are analyzed by the finite element method in plain strain condition.
CGP
CSC
The results show that the CSC method leads to better hardness distribution, uniform strain dispersal and
FEM finer surface quality as compared to conventional method. Moreover, the new method is a good candidate
Mechanical properties for substituting the conventional CGP process.
Strain distribution © 2012 Elsevier B.V. All rights reserved.
1. Introduction and lower grooved dies as shown in Fig. 1a. In the first cycle, the
sheet sample is pressed by grooved die (Fig. 1b). Then in the sec-
In the last decade, the importance of severe plastic deformation ond cycle, the grooved work-piece is straightened with the flat die,
(SPD) techniques is extensively recognized for all researchers in the see Fig. 1c. The theoretical relationships for estimating shear and
materials science field [1]. Various industries are now interested effective strains after two cycles are shown in Eqs. (1) and (2),
to use ultra-fine grained (UFG) or nano-structured (NS) materi- respectively [15]. After second cycle, the specimen is rotated by
als because of their high physical, mechanical and superplastic 180◦ around the axis perpendicular to the sheet plane; see Fig. 1d.
properties [2–4]. In general, SPD methods can be divided into This allows the non-deformed regions to be deformed by further
two categories: first group such as equal channel angular pressing pressing due to the asymmetry of the groove die (Fig. 1(e)) and
(ECAP) [5], accumulative back extrusion (ABE) [6], tubular chan- finally the flat die (Fig. 1(f)) results in a relatively uniform distribu-
nel angular pressing (TCAP) [7] and high pressure torsion (HPT) tion of plastic strain throughout the specimen.
[8] has been used for bulk-form materials. Second group includ-
t
ing accumulative roll bonding (ARB) [9], repetitive corrugation and xy = =1 (1)
t
straightening (RCS) [10], groove rolling (GR) [11], and constrained
groove pressing (CGP) [12] has been applied for sheet-form metals 2
4(xy /2) 1
and alloys. In all techniques, an intensive shear strain is imposed εeff = = √ = 0.557 (2)
3 3
to the specimen and deformed work-piece is obtained without any
change in the initial specimen dimensions. This characteristic is the So far, several experimental and numerical works have been
distinctive property of the SPD processes [2,13,14]. carried out on the CGP process. The experimental and simu-
CGP process firstly developed by Shin et al. [12] has the poten- lated evaluation of CGP and CGR (constrained grooved rolling)
tial for fabrication of sheet-form metals and alloys. The schematic techniques for the pure Al had been explored by Lee and Park
representation of CGP process and the sequences of pressing have [16]. They have reported that a non-uniform strain distribution is
been displayed in Fig. 1. Each pass is consisting of four cycles, two observed in the thickness direction and the effective strain mag-
for asymmetrically grooved die and the other two for the flat die. nitude at the center portion of sheet is greater than that at the
At first, a sheet specimen is placed in the gap between the upper surface. In two separate researches, microstructure and mechani-
cal properties of CGPed pure copper and pure aluminum had been
surveyed by Krishnaiah et al. [17,18] at two ambient and cryo-
∗ Corresponding author. Tel.: +98 21 77240203; fax: +98 21 77240203. genic temperatures. No significant effect was observed for these
E-mail address: [email protected] (F. Djavanroodi). two temperatures. Also the results show that the grain size is
0921-5093/$ – see front matter © 2012 Elsevier B.V. All rights reserved.
http://dx.doi.org/10.1016/j.msea.2012.04.121
Author's personal copy
reduced from 78 m to a cell block structure approximately 0.5 m 2D stress and 2D strain conditions. In the bio-field research, the
for pure Cu and from 38 m to 1 m for pure Al. Peng et al. [19] investigation of Thirugnanam et al. [26] had displayed that good
showed that a wider groove reduces the stress at the most severely wettability and higher surface energy are obtained with CGPed
loaded portion of the sheet and allows more cycles of pressings to specimens as compared to un-CGPed ones.
attain better grain refinement and improved mechanical behavior. In this research, commercial purity aluminum is subjected to
The strain distribution in the CGPed sample had been investigated CGP process up to failure in three different die conditions as shown
by Yoon et al. [20] using the commercial Elasto-Plastic FEM code in Table 1. Scanning electron microscopy (SEM) and XRD patterns
ABAQUS/Standard. The main finding is that the strain localization using classical Williamson–Hall technique are applied to mea-
is more intensified with the number of passes, although average sure grain size magnitudes to ensure that the material after CGP
effective strain level increases. Introduction of cross-orientation process have the UFG structure. Mechanical properties including
CGP process on the 38 Zn alloy showed the permission of more pass strength, elongation to failure and hardness are obtained for all
numbers and thus higher magnitude of effective strain, finer grain above three conditions. Finally finite element analysis, after veri-
size and higher mechanical properties without any crack formation fication by experimental work, is utilized to investigate the strain
in the sample are obtained as compared to single-orientation CGP behavior and required force value for all above conditions.
process [21]. Hosseini and Kazeminez [22] had showed that CGP
process decreases work-hardening rate and increases strain rate
sensitivity of samples. The effect of intermediate and post anneal- 2. Experimental procedure
ing on the mechanical and microstructure of CGPed Cu sheet had
been studied by Rafizadeh et al. [23]. The result indicated that the Commercial purity aluminum with the chemical composition
influence of first pass on the strength of material is much more than of (in wt.%) 0.262% Fe, 0.0696% Si, 0.0214% Ti, 0.0201% Cu, 0.0134%
the subsequent passes. CGP of low carbon steel were investigated Zn, 0.0076% Ni, 0.0072% V, 0.006% Mg and Al as the balance were
by Khodabakhshi et al. [24]. They showed that the inhomogeneity prepared in 50 mm × 50 mm. Before pressing, all sheet-formed
factor value is decreased with increasing the pass number indicat- samples were annealed at 370 ◦ C for 1 h and then cooled slowly
ing increase of strain homogeneity. Niranjan and Chakkingal [15] at furnace. The initial average grain size, obtained by SEM, is about
had studied the deep drawability of pure aluminum before and after 15 m as shown in Fig. 2.
CGP process. The initial sheets for pressing are produced by 0◦ , 45◦
and 90◦ related to rolling direction. The LDR (limit drawability ratio)
enhances to a maximum value for the CGPed work-piece at 45◦ to Table 1
rolling direction orientation as compared to 0◦ and 90◦ orienta- The dies and specimen conditions used in this research.
tions. Shirdel et al. [25] investigated the semi-constrained groove Die condition Die parameters Sheet dimensions
pressing (SCGP) on pure Al. Increasing the number of passes leads
(◦ ) t (mm) Cross-section (mm2 ) Thickness (mm)
to accumulation of strain in the material and raises the level of
flow stress and sheet hardness. Also, comparison of force diagrams Die 1 45 1 50 × 50 1
obtained from 2D plane stress, 2D plane strain, 3D and experimen- Die 2 45 1.8 50 × 50 0.4Cu + 1Al + 0.4Cu
Die 3 53 1 50 × 50 1
tal work indicated that the real state of SCGP is a combination of
Author's personal copy
For this research, three CGP dies were designed and manu-
factured. The die parameters and sheet dimensions are listed in
Table 1. Fig. 3 represents the three CGP dies and corresponding
deformed samples after one cycle. A container used to prevent
the material flow in the direction parallel to groove (not shown
in Fig. 3). For die 2 condition, 1 mm thick aluminum sheet is shel-
tered by the 0.4 mm thick copper sheet (CSC method) was used to
prevent micro-cracking in the surface of Al sheet. After one pass,
the sample was rotated 90◦ about its normal axis (cross-orientation Fig. 5. The die 1 required force curves to press Al sheet in plane strain, plane stress
and 3D simulations.
[21]) and subsequent pass was performed up to failure. All press-
ing were carried out with a hydraulic press with a ram speed of
0.1 mm s−1 . Molybdenum disulfide (MoS2 ) was used as a lubricant 0.1 mm s−1 same as the experimental procedure and friction coef-
to decrease frictional effects between specimens and dies [26]. ficient was assumed to be 0.1 [16]. Fig. 4 displays true stress-true
After pressing, the mechanical properties including tensile test strain curve of annealed aluminum, obtained by tensile test, which
and Vickers hardness according to ASTM E8M and ASTM E92 and is imported to the software. Adaptive meshing, mass scaling and
grain size using XRD pattern are obtained for all three constructed automatic re-meshing were applied and justified for all of the sim-
die work-pieces after each pass. The hardness tests were performed ulations to prevent failure of the mesh during large deformation
on the centerline of the CGPed sheet by the increment of about and to reduce computation time. Mesh sensitivity was performed
5 mm. and the optimum element size was found. The numbers of mesh
for aluminum sheet and cover copper sheet were 3200 and 1600,
2.1. Finite element analysis and verification respectively.
Verification of FEM was based on required force value to
The commercial finite element code, ABAQUS/Explicit software, press aluminum sheet for the die 1 condition. Plane stress, plane
is applied in order to investigate strain behavior, homogeneity of strain and 3D situations are simulated and the extracted required
strain distribution and required force magnitude in the CGPed spec- force curves are presented in Fig. 5. Comparison of simulated
imens. Both flat and grooved dies were assumed to be discrete and experimental results indicates that the plane strain and 3D
rigid (no deformation). The constant upper die speed was equal to situations are in good agreement with experimental work and due
Fig. 3. The CGP dies and corresponding deformed sheets used in this study.
Author's personal copy
100 A. Sajadi et al. / Materials Science and Engineering A 552 (2012) 97–103
Fig. 6. Die 1 first cycle corner teeth filling for plane stress, plane strain and 3D conditions.
to computation running time consideration, plane strain situation pass 3. The same trends had been reported in previous studies [22].
is chosen to simulate FEM work. Finally, the teeth corner filling Also, the magnitude of elongation is decreased with increasing pass
at the first cycle of CGPed Al for the die 1 situation are considered number and this tendency is started with high incline in the first
for plane stress, plane strain and 3D conditions, as shown in Fig. 6. pass. The magnitudes of YS and UTS at final pass in die 1 situa-
It can be observed that for plane stress situation, the corner teeth tion are 98 MPa and 110 MPa, respectively. Die 2 configuration is a
of die 1 configuration is not fully filled with material; but, both new method (CSC) based on sandwiching the Al work-piece with
plane strain and 3D situations have the same behavior in filling copper sheets; so that surface micro-cracking can be prevented
die corners which is seen in the experimental CGPed sample. or postponed as the number of passes increases. With this new
The experimental load magnitudes of 157 kN, 188 kN and 175 kN method, no decreasing in the YS and UTS values at the fourth pass
are achieved for dies 1, 2 and 3 configurations at first cycle, is observed. High magnitudes of YS and UTS (119 MPa and 141 MPa)
respectively, which are used to verify the numerical works. The are achieved with die 2 configuration. In die 3, the grooved angle is
results show that the percentage difference between FEM and 53◦ and the failure of sheet occurred after three passes. Lower mag-
experimental achievement is lower than 10% for all three die nitudes of YS and UTS are obtained with this die condition (80 MPa
conditions which is acceptable for engineering works. and 87 MPa). In the other words, increasing grooved angle leads to
After validation, the three experimental die situations (die 1, increase stresses in the specimen at the corner of the teeth and thus,
die 2 and die 3) were simulated up to 1 pass (4 cycles) and the enhance the cracking tendency. Hence, higher grooved angle does
magnitude of effective strain, homogeneity of effective strain and not permit more cycles of pressing and enhancement of mechan-
required force curves to press sheet are obtained and discussed. ical properties. In general micro-cracking and flow softening are
two mechanisms which can influence the strength reduction of
3. Results and discussion deformed sample before failure [22]. The UTS and YS values after
third pass for the aluminum specimen which was covered with
3.1. Microstructure evaluation copper sheet (CSC method – die 2) shows increase unlike other
dies 1 and 3 at this pass number. Also, for the new method a bet-
The annealed aluminum sheet exhibits equiaxed grain size ter surface quality than the ones in the conventional CGP process
structure. To confirm the reduction of grain size, final pass of have been observed. Surface condition is an important parameter
CGPed Al for all three die conditions were examined by Philips X- in micro-cracking mechanism, using copper sheet which has a sac-
ray diffraction (XRD) and the results were analyzed using classical rificing role for the Al sample can promote the surface condition
Williamson–Hall method. The details of this method can be studied and cause to delay strength dropping due to micro cracking. Thus,
elsewhere [24]. The average grain size of as-received and CGPed the strength dropping does not occur in 3rd pass as observed in dies
aluminum sheets at final pass for dies 1, 2 and 3 conditions are 1 and 3.
15 m, 870 nm, 865 nm and 890 nm, respectively. Table 2 compares the percentage changes in yield and ultimate
tensile strengths for Al sheets produced in each pass number with
respect to initial condition for three die geometries. In all die con-
3.2. Mechanical properties
ditions, increases of about 100% for YS and 50% for UTS have been
observed in the first pass. The increasing of strengths continues
After pressing the aluminum sheet samples up to failure, ten-
sile and hardness tests have been performed to obtain mechanical
properties of sheets for each pass in the three die conditions. It is Table 2
noted that the Al sheets in dies 1 and 2 with grooved angle of 45◦ The percentage changes in yield and ultimate tensile strengths for Al sheets pro-
were pressed up to 4 passes and in the die 3, where the grooved duced in each pass number with respect to initial condition for three die geometries.
angle is 53◦ , the specimen was pressed up to 3 passes. Fig. 7 dis- YS UTS
plays the magnitudes of yield strength, ultimate tensile strength
and elongation to failure versus the pass numbers in all die situa- Die 1 Die 2 Die 3 Die 1 Die 2 Die 3
tions. For die 1 condition, it can be observed that increasing number 1st pass 105 125 112 44 58 49
of passes results in high magnitudes of yield and ultimate tensile 2nd pass 154 150 151 60 81 53
3rd pass 177 198 99 67 88 16
strengths up to 3 passes. YS and UTS are increased with high slope
4th pass 148 198 – 46 87 –
in the first pass. For pass 4, strength is decreased as compared to
Author's personal copy
A. Sajadi et al. / Materials Science and Engineering A 552 (2012) 97–103 101
with lower percentage magnitude in subsequent passes and then where Hi , H̄ and n are the magnitude of hardness for ith location,
decreases. average value of hardness measurements and the number of hard-
Considering that the strength and hardness are proportional, ness tests for each sheet, respectively. In general, less value for
the same trend is anticipated and the results of hardness measure- I.F. leads to better homogeneity of hardness distribution. Fig. 10
ment confirm that. Fig. 8 shows the average magnitude of CGPed Al displays the I.F. values of hardness measurements at each pass
sheet hardness for each pass in the three die conditions. Hardness numbers for the three die configurations.
Fig. 9. The hardness distribution for all dies conditions at each pass number.
Fig. 10. The inhomogeneity factor (I.F.) for all three die conditions versus pass numbers.
Author's personal copy
102 A. Sajadi et al. / Materials Science and Engineering A 552 (2012) 97–103
Fig. 12. The magnitudes of average effective strain and standard deviation for three
die conditions after one pass CGP.
Fig. 11. The effective strain contours for three die conditions after 1 pass CGP: (a)
die 1 condition, (b) die 2 condition and (c) die 3 condition.
A. Sajadi et al. / Materials Science and Engineering A 552 (2012) 97–103 103
4. Conclusion [2] R.Z. Valiev, T.G. Langdon, Prog. Mater. Sci. 51 (2006) 881–981,
http://dx.doi.org/10.1016/j.pmatsci.2006.02.003.
[3] F. Djavanroodi, M. Ebrahimi, Mater. Sci. Eng. A 527 (2010) 1230–1235,
The plastic deformation behavior of commercial purity alu- http://dx.doi.org/10.1016/j.msea.2009.09.052.
minum processed by constrained groove pressing was investigated [4] R.B. Figueiredoa, T.G. Langdon, Mater. Sci. Eng. A 503 (2009) 141–144,
using 2D FEM and experimental method. The process was per- http://dx.doi.org/10.1016/j.msea.2008.04.080.
[5] S.D. Wu, Z.G. Wang, C.B. Jiang, G.Y. Li, I.V. Alexandrov, R.Z. Valiev, Mater. Sci. Eng.
formed in the three configurations utilizing three dies with A 387–389 (2004) 560–564, http://dx.doi.org/10.1016/j.msea.2003.12.087.
different conditions. Three CGP dies (die 1 with the = 45◦ and [6] S.M. Fatemi-Varzaneh, A. Zarei-Hanzaki, Mater. Sci. Eng. A 504 (2009) 104–106,
t = 1 mm; die 2 with = 45◦ and t = 1.8 mm; die 3 with = 53◦ and http://dx.doi.org/10.1016/j.msea.2008.10.027.
[7] G. Faraji, M.M. Mashhadi, H.S. Kim, Mater. Lett. 65 (2011) 3009–3012,
t = 1 mm) were designed and manufactured. For each condition, CGP
http://dx.doi.org/10.1016/j.matlet.2011.06.039.
process has been carried out up to failure (4 passes for dies 1 and [8] T. Hebesberger, H.P. Stüwe, A. Vorhauer, F. Wetscher, R. Pippan, Acta Mater. 53
2; 3 passes for die 3). The results of conventional and covered sheet (2005) 393–402, http://dx.doi.org/10.1016/j.actamat.2004.09.043.
[9] M. Shaarbaf, M.R. Toroghinejad, Mater. Sci. Eng. A 473 (2008) 28–33,
casing (CSC) methods were compared. The following conclusions
http://dx.doi.org/10.1016/j.msea.2007.03.065.
were drawn: [10] J. Huang, Y.T. Zhu, D.J. Alexander, X. Liao, T.C. Lowe, R.J. Asaro, Mater. Sci. Eng.
A 371 (2004) 35–39, http://dx.doi.org/10.1016/S0921-5093(03)00114-X.
1. Using of XRD patterns and classical Williamson–Hall method, it [11] S.M. Byon, D.H. Na, Y. Lee, J. Mater. Process. Technol. 209 (2009) 4465–4470,
http://dx.doi.org/10.1016/j.jmatprotec.2008.10.045.
is confirmed that the pure aluminum sheets subjected to the CGP [12] D.H. Shin, J.-J. Park, Y.-S. Kim, K.-T. Park, Mater. Sci. Eng. A 328 (2002) 98–103,
process have a UFG structure with grain sizes of 15 m, 870 nm, http://dx.doi.org/10.1016/S0921-5093(01)01665-3.
865 nm and 890 nm for as-received, final pass for dies 1, 2 and 3 [13] S.C. Baik, Y. Estrin, H.S. Kim, R.J. Hellmig, Mater. Sci. Eng. A 351 (2003) 86–97,
http://dx.doi.org/10.1016/S0921-5093(02)00847-X.
conditions, respectively. [14] S.M. Fatemi-Varzaneh, A. Zarei-Hanzaki, M. Naderi, A.A. Roostaei, J. Alloys
2. The ultimate tensile strength increased to 110 MPa after the 4th Compd. 507 (2010) 207–214, http://dx.doi.org/10.1016/j.jallcom.2010.07.157.
pass for die 1 situation. Highest magnitudes of YS (119 MPa) [15] G.G. Niranjan, U. Chakkingal, J. Mater. Process. Technol. 210 (2010) 1511–1516,
http://dx.doi.org/10.1016/j.jmatprotec.2010.04.009.
and UTS (141 MPa) are achieved with die 2 condition (CSC–CGP [16] J.W. Lee, J.J. Park, J. Mater. Process. Technol. 130–131 (2002) 208–213,
method). With this method, no decreasing in the YS and UTS http://dx.doi.org/10.1016/S0924-0136(02)00722-7.
values at the 4th pass is observed. Lower magnitudes of UTS and [17] A. Krishnaiah, U. Chakkingal, P. Venugopal, Mater. Sci. Eng. A 410–411 (2005)
337–340, http://dx.doi.org/10.1016/j.msea.2005.08.101.
number of passes are obtained with die 3 condition (87 MPa) at
[18] A. Krishnaiah, U. Chakkingal, P. Venugopal, Scripta Mater. 52 (2005) 1229–1233,
the 3rd pass. Higher grooved angle does not promote more cycles http://dx.doi.org/10.1016/j.scriptamat.2005.03.001.
of pressing and improvement of mechanical properties. [19] K. Peng, L. Su, L.L. Shaw, K.W. Qian, Scripta Mater. 56 (2007) 987–990,
http://dx.doi.org/10.1016/j.scriptamat.2007.01.043.
3. The experimental results showed that by increasing the num-
[20] S.C. Yoon, A. Krishnaiah, U. Chakkingal, H.S. Kim, Comput. Mater. Sci. 43 (2008)
ber of passes, the specimen hardness increased from 24 VHN to 641–645, http://dx.doi.org/10.1016/j.commatsci.2008.01.007.
about 50 VHN for all the dies conditions. [21] K. Peng, Y. Zhang, L.L. Shaw, K.W. Qian, Acta Mater. 57 (2009) 5543–5553,
4. Better homogeneity of mechanical properties is obtained with http://dx.doi.org/10.1016/j.actamat.2009.07.049.
[22] E. Hosseini, M. Kazeminezhad, Mater. Sci. Eng. A 526 (2009) 219–224,
increasing number of passes. On the other hand, CSC–CGP http://dx.doi.org/10.1016/j.msea.2009.07.028.
method gave the best uniformity in mechanical properties. [23] E. Rafizadeh, A. Mani, M. Kazeminezhad, Mater. Sci. Eng. A 515 (2009) 162–168,
5. As a final result, the CSC method is a good candidate for substi- http://dx.doi.org/10.1016/j.msea.2009.03.081.
[24] F. Khodabakhshi, M. Kazeminezhad, A.H. Kokabi, Mater. Sci. Eng. A 527 (2010)
tuting the conventional CGP process. 4043–4049, http://dx.doi.org/10.1016/j.msea.2010.03.005.
[25] A. Shirdel, A. Khajeh, M.M. Moshksar, Mater. Des. 31 (2010) 946–950,
References http://dx.doi.org/10.1016/j.matdes.2009.07.035.
[26] A. Thirugnanam, T.S. Sampath Kumar, U. Chakkingal, Mater. Sci. Eng. C 30 (2010)
203–208, http://dx.doi.org/10.1016/j.msec.2009.10.002.
[1] Y. Iwahashi, Z. Horita, M. Nemoto, T.G. Langdon, Acta Mater. 46 (9) (1998)
3317–3331, http://dx.doi.org/10.1016/S1359-6454(97)00494-1.