Analysis of Void Growth During Superplastic Deformation of Commercial AL5083 Alloy

Superplastic alloys and metals possess the ability to undergo large uniform strains prior to failure. A number of materials are subject to the cavitation during superplastic deformation. Cavitation usually leads to either the undesirable post-forming characteristics or to the premature tensile failure. It is also apparent that the cavities can preexist in the form of cracks and decohered interfaces, which develop during thermo-mechanical processing necessary to produce the superplastic microstructures. Evidently, extensive cavitation imposes significant limitations on their commercial application. The material constitutive equation constants of commercial AL5083 alloy contain strength coefficient, and strain rate sensitivity index is determined by superplastic bulge forming tests for 400, 450, 500, 550 °C. By comparing the results of a deformed sample, good accordance between experimental and FEM results is observed. Using the calculated values for C and m parameters, the effect of material properties such as the cavity growth rate, strain rate sensitivity index, strain hardening exponent and number of intentionally preexisted voids on specimens voids growth subject to the tensile deformation and to the biaxial deformation has been determined numerically. The tensile tests have been simulated by the FEM software ABAQUS v6.9 using commercial aluminum 5083 alloy that presents superplastic properties at temperatures in the range 400–550 °C. The simulations are implemented at 400, 450, 500 and 550 °C. The results of the numerical prediction obtained are in good agreement with the results of the experiments done by some other authors.


Introduction
The reason for the spread of superplastic materials is due to their exceptional ductility in particular conditions. Indeed, in some cases tensile elongations above 7000% were achieved, even if, in industrial applications typical values are comprised between 200 and 1000% (Langdon 1995). It is important to note that such values are greater by one or two orders of magnitude compared to those of conventional metallic materials. As regards isothermal superplasticity, which is the one considered in industrial cases, the conditions that must be respected for a metal or a metallic alloy to exhibit superplastic behavior regard both material microstructure and the operating parameters adopted during the forming process. In particular, the material must have a grain size of less than 10 lm, and forming must take place at a temperature that is constant and equal to at least half the absolute melting point, and at a very low strain rate, in the order of 10 À5 to 10 À3 s À1 (Langdon 1995). It is observed that the value of the strain rate sensitivity index (m) has a strong effect on the ductility of superplastic materials. In general, the higher the m value, the greater the elongation to failure (Kim et al. 1996;Carrino et al. 2001Carrino et al. , 2003. Shehata et al. (1978) examined the formability of several Al-Mg alloys at temperatures from room temperature to 300°C over a wide range of strain rates by performing uniaxial and biaxial stretch forming tests. Naka et al. (2001) investigated the effect of forming speed and temperature on the formability of AA5083 alloy sheet by stretch forming tests with a flat head cylindrical punch at various forming speeds and temperatures from room temperature to 300°C. Hosseinipour (2009) studied the hot deformation behavior of AA5083 with tensile tests at various temperatures and strain rates. The results showed that the formability increased with decreasing speed for any strain paths at high temperature, while at room temperature it was not as sensitive to speed. Hosseinipour (2010) investigated the strain rate sensitivity and cavitation in a commercial 5083 aluminum alloy. The results showed that with increasing temperature the maximum strain rate sensitivity decreases and shifts to the lower strain rates. The failure surface is wide and failure occurs by cavitation. Chung and Cheng (2002) showed that, in some superplastic materials, the fracture mode is dominated by unstable plastic flow. The instability of superplastic deformation was studied by several investigators by analytical approaches. Pearce (1989) showed that in the tensile test the shrinkage rate is inversely proportional to the cross section of the specimen and highly sensitive to the strain rate sensitivity index.
Strain rate control assumes a fundamental role. So that the strain rate sensitivity index, m, attains a maximum value. This condition permits high material formability, m, being directly proportional to the material elongation capability.
This high ductility can be easily exploited to form complex shapes making superplastic forming technology of notable interest for the metal forming industry, particularly for the aerospace applications. This ability to reach large strains is limited by the fact that such materials are susceptible to internal voids formation, which may eventually lead to failure in the specimen. The cavitation process is comprised of three distinct stages, which in most cases occur simultaneously: cavity nucleation, cavity growth and cavity coalescence.
Cavities, which preferentially nucleate at irregularities of the grain boundary, triple point or second-phase particles, grow by either plastic deformation (Hancock 1976;Cocks and Ashby 1980) or diffusion-controlled mechanism (Chokshi and Langdon 1990), or a combination of the two. The cavity coalescence occurs during the last stages of the deformation process by which time a large cavity volume fraction has developed (Pilling and Ridley 1989). The plasticity-controlled cavity growth is described by the following equation (Pilling and Ridley 1989): where r is the cavity radius, e the strain, c the surface energy of cavity, r e the effective stress and g ¼ d ln m=de is the cavity growth rate factor, where m is the volume of a single cavity. This volumetric growth rate factor, g, described by Cocks and Ashby (1980), Stowell et al. (1984) and Pilling and Ridley (1989), in terms of strain rate sensitivity (m), the applied stress state and stress concentration factor (k s ), is expressed as follows: where P is the imposed pressure, k s depends on the geometry of deformation and the extent of grain boundary sliding (Cocks and Ashby 1980;Pilling and Ridley 1989) and: where r m is the mean stress. Stowell (1980) by assuming a certain level of preexisting voids within the material gives: where C v is the cavity volume fraction and C 0 the initial volume fraction of cavities. To reduce the level of cavitation for a given strain, it is necessary to reduce the magnitude of g, the cavity growth parameter, which depends on m and the ratio r m =r e (Eqs. 2 and 3). Diffusional cavity growth occurs by the stress-directed diffusion of vacancies in the cavities, usually along grain boundaries. This model predicts the following growth rate: where r is the cavity radius, dr=de the cavity growth rate per unit strain, X the atomic volume, d the grain boundary width, D gb the coefficient for grain boundary diffusion, k the Boltzmann's constant, T the absolute temperature and c the surface energy (Chokshi 1986). In the early stages, the voids shape of deformation is elliptical, the hole increases in length along the tensile axis but it decreases along the perpendicular direction to the tensile axis (transverse contraction), subsequently, the hole grows both along and perpendicular to the tensile axis (transverse growth), and finally, a crack nucleates on either side of the hole and it propagates to cause failure (crack propagation) (Chokshi and Langdon 1996). A special aspect of the growth problem is the effect of voids interaction on the macroscopic mechanical behavior (ductility and flow stress) of the investigated material. The numerical and experimental studies have revealed that the parameters affecting voids growth in the superplastic deformation are the strain rate sensitivity (m), the voids size and the number of the existing voids (Khraishi et al. 2001). Carrino et al. (2004) investigated the effects of some material properties such as void growth rate, void numbers, strain rate sensitivity index and strain hardening exponent on void growth and material ductility of fine-grained Pb-Sn alloy that presents superplastic behavior at room temperature by experimental procedure and FEM at tensile deformation and biaxial deformation. They stated that increasing the void numbers and decreasing of m and n leads to decrease in void growth and increase in material ductility. Lin et al. (2005) collected contents for different damage mechanisms of metallic materials at various deformation conditions of high-temperature creep, cold metal forming, superplastic forming and hot metal forming. They reviewed and discussed the ruling equations and damage calibration techniques for these deformation conditions. Mulholland et al. (2006) conducted a set of parametric experiments of simple superplastic tensile test on a pre-machined superplastic eutectic tin lead alloy with one or more holes. Their results indicated an increase in ductility by increasing hole numbers up to 10 holes, and decrease in ductility by increasing the hole numbers more than 10 holes. Hosokawa et al. (2013) experimentally investigated the influences of material work-hardening behavior of some materials on fracture mechanism and void growth and coalescence, using in situ X-ray computed tomography (XCT) along with tensile test on specimens including artificial void array. Also, they used Thomason, Pardeon and Hutchinson models for void coalescence.
In this paper, at first the material constitutive equation constants including strength coefficient and strain rate sensitivity index are calculated by gas blow forming bulge tests at different superplastic hot temperatures 400, 450, 500, 550°C for commercial AL5083 alloy. Then, by using these values, a comparison of a sample dome heights between experimental and FEM results has been made. A good agreement between two sets of results has been obtained. Similar to work done by Carrino et al. (2004), the effect of material parameters (m and n), g and number of voids on the voids growth is analyzed with finite element method. In the case of uniaxial tension 1, 3 and 5, preexisted voids have been considered and an only void in the case of biaxial tension has been considered. Finally, the true stress-strain curves have been investigated for mentioned hot temperatures.

Experimental Material and Procedures
The simple tensile test samples are clamped at the bottom side and at the top side are stretched with a velocity in accordance with the strain rate defined. This velocity can be calculated using Eq. (6): where L is the instantaneous length of the sample and _ e is the constant strain rate. A power law form of the constitutive relationship is assumed: where r is the effective stress, _ e is the effective strain rate, C is a constant, and m is the strain rate sensitivity index.
The material characteristics in the hot condition were determined by bulge tests with the gas blow forming process at various temperatures and pressures using the following equation (Giuliano and Franchitti 2008;Sorgente and Scintilla 2010): where t 1 and t 2 are the forming times necessary to obtain the same dome height at constant pressures of p 1 and p 2 , respectively. The dome height during the bulge test was measured by an ultrasonic telemetry sensor. The equivalent stress, strain and strain rate in the dome apex were calculated using the following equations (Liu et al. 2008;Koç and Billur 2011): where: and e ¼ ln s s 0 and where r is the dome radius, s is the final thickness of the dome apex, s 0 is the initial thickness of the blank, h is the dome height, d is the die diameter, and t is the forming time.
In order to validate the finite element simulations, a comparison between experimental test and FEM simulation has been made. The experiments were conducted using a commercial AL5083 alloy, which is known to exhibit superplastic properties at temperatures 400-550°C. The bulge tests were carried out using a set of equipment ( Fig. 1) at constant pressures 0.2 MPa and 0.3 MPa to establish the constitutive equation of the material. If we assume these two pressures as p 1 and p 2 , for a constant value of bulge height, the times to reach this bulge height are recorded as t 1 and t 2 . The strain rate sensitivity, m, is calculated by Eq. (8) (Giuliano and Franchitti 2008). The material constitutive equation is assumed as Eq. (7).
These tests have been implemented for temperatures 400, 450, 500 and 550°C, and strain rate sensitivity and strength coefficient were calculated (Table 1). Figure 2 shows the comparison between two samples bulged at constant pressure 2 bar and at temperature 500°C, FEM is simulated on top and experimental sample at the bottom. Figure 3 shows the comparison of dome height of sample during sheet forming, between experimental and simulation results. The figure shows the good agreement between experimental and numerical simulation results.

Numerical Analysis
In this work, the numerical studies of the growth of the voids and the parameter that affect it were carried out. The finite element calculations were performed using the ABAQUS v6.9 software package. Initially, each finite element model was used to simulate a thin plate (plane stress conditions) with preexisting holes subjected to a constant pressure. Figure 4 shows the meshes of tensile specimens with three different hole configurations: a central hole, three equally spaced holes and five equally spaced holes. The 2D modeling, using four-noded continuum elements, was examined. The initial holes have the form of a circle and a radius that is equal to 640 lm.
Due to the symmetry of the loading and the geometry, only the first quadrant was used in the calculations. The single-hole mesh had 235 quadrilateral elements with 277 nodes. The three-hole mesh had 505 elements with 568 nodes, and finally, the five-hole mesh had 825 elements with   In the case of biaxial deformation, a single hole was realized in the center of a sheet with a width of 150 mm and, due to the symmetry, the first quadrant was examined. The mesh had 105 quadrilateral elements for axisymmetric applications with 212 nodes. The FEM mesh and the die are shown in Fig. 5. The nodes along the die entry radius are fixed in the y-direction except for one, which is fixed in the xdirection in order to simulate the presence of a blank-holder. The element type CPS4R was used in the FEM simulations. This is a 4-node bilinear plane stress quadrilateral, reduced integration, hourglass control element. The element shape is quad, and the technique of meshing is structured.
A constitutive equation, which defined the relationship between the flow stress, r, the strain, e and the strain rate, _ e, was employed as: where k is the strength coefficient, m the strain rate sensitivity index and n the strain hardening index. The critical variables in the present model are the material parameters (m and n) and the strain. The m value has been varied between 0.25 and 0.35 (0.25, 0.32, 0.342 and 0.35), and n had the following values: 0, 0.0334, 0.0667 and 0.1. The voids growth was valued, considering the parameters: normalized minor radius r min =r 0 , normalized major radius r max =r 0 , normalized area A=A 0 and the void growth rate g, depending on elongation(%), m and n. If we assume C v =C 0 ¼ A=A 0 , according to Eq. (4), it can be shown that: 4 Results and Discussions Figures 6 and 7 show the material constants at elevated temperatures and various strain rates which were obtained by bulge tests using the gas blow forming process. The material constants were calculated at different dome heights, and the mean values were reported. As can be seen, with increasing temperature the m value increased, and the maximum m value was obtained at 500°C, but the C constant decreased with increasing temperature. As shown in Fig. 8, the results of numerical analysis indicate that, in the case of a single hole preexisted, the growth of the voids increases with increasing sample elongation and depends on m and n parameters. The void growth (measured by A=A 0 ) is relatively small for high m values, fixed n = 0, and always becomes important for small m values (m = 0.25).
Analogously, as shown in Fig. 9, fixed m = 0.25, the void growth evolves with decreased n values as also reported in other works (Khraishi et al. 2001;Carrino et al. 2003). By increasing the values of true strain, the void growth rate parameter, g, at first decreases and then begins to increase. The variation tilt of void growth rate, at bigger n values is low, but at low n values the g parameter varies with higher intensity. The g values decrease by increase in   The effect of multiple holes was analyzed, considering three and five holes aligned along the major axis. The tensile test output is represented by the variation of normalized minor and major radii of voids with sample elongation, normalized void area with elongation and void growth rate parameter, g, with true strain. For brevity, only the curves relative to five holes, fixed m = 0.25, are reported in Fig. 9.
The performance of the multiple-hole curves is the same as the single-hole curves with varying m and n parameters, but by increasing the number of holes, along the major axis the effect of increasing the metal ductility is observed. This   Fig. 9 The effects of strain hardening index on the growth of the five voids m = 0.25) occurs when the extent of necking is reduced and its onset delayed through dispersing the plastic deformation at the different hole sites and also reducing the stress concentration. Besides, the elongation of the major radius is restricted by neighboring holes and cannot achieve values attained by a single-hole specimen for the same given strain. Figure 10 shows a decrease in a major radius and a neck-in, in accordance with the holes, by increasing the number of holes. Those results are in good qualitative agreement with the results of Khraishi et al. (2001). The effect of multiaxial modes of deformation was investigated, varying m and n parameters and valuing the single void growth (measured by normalized area) depending on true strain (measured by the height of free forming specimens). In the case of the biaxial deformation, the form of the hole does not change during the superplastic forming. Fixing n, the normalized area increases by decreasing the m parameter, which is obviously shown in Fig. 11 which shows that at high m values ðm ! 0:32Þ, the void growth procedure happens gradually, but at m ¼ 0:25 which refers to the temperature of 400°C, by increasing the true strain value, void growth occurs with high tilt and intensity. The results of numerical analysis are shown in Fig. 11. Figure 12 shows the deformed mesh, fixed h = 17 mm and n = 0, for the m values of 0.25 and 0.35, and the major void growth is observed at m = 0.25. Besides, the numerical results indicate that the radius of the hole is increased by increasing the height (true strain) of the specimen. High ductility increasing can be seen by the increase in m value.
In Fig. 13, the evolution of the void growth is reported in four different levels of deformation. Comparing uniaxial and biaxial deformation (Fig. 14) in the case of a single void, the void growth (measured by the normalized area) is more important for the biaxial deformation by increasing the true strain. Figure 15 shows the parameter void growth rate, g, for different m values. For all m values, the void growth parameter g sharply increased at low strains, but after reaching a maximum value, g began to decrease by plastic deformation evolvement. g values rose with decrease in m parameter. Figure 16 shows the evolution of void growth in simple tensile test. At first, the void growth velocity is low, but suddenly begins to grow fast.
The stress-strain curves at various temperatures are illustrated in Fig. 17. It appears that at lower strain after high initial stress, flow stress begins to gradually fall with strain, and at higher strain, the flow stress reaches a steady-state stress. This kind of strain softening is normally observed when the material exhibits continuous recrystallization during hot deformation where the deformed grains are replaced by strain-free grain.

Conclusion
The material constitutive equation constants, strength coefficient (c) and strain rate sensitivity index (m), were obtained using gas blow forming bulge tests at different superplastic hot temperatures 400, 450, 500, 550°C for commercial AL5083 alloy. Using these values and by comparing sample dome heights between experimental and FEM results, good agreement between these two sets of results was observed. Finite element software ABAQUS v6.9 was used to investigate the parameters affecting the growth of the preexisted voids in the commercial AL5083 alloy by the numerical tests in two stress states, uniaxial superplastic tensile test and biaxial superplastic bulge test. The stress-strain curves at various temperatures were illustrated and discussed. The m values increased by temperature evolving at 500°C; the maximum m values were obtained, but the strength coefficient decreased by increase in temperature. The voids growth increased with elongation increasing. The void growth is independent from the strength coefficient value and the applied strain rate, but it is dependent on the strain rate sensitivity and the strain hardening index. It is also small for high m values and important for small m values. The void growth rate    The stress-strain curves at various temperatures parameter, g, decreases with an increase in n. The effect of strain hardening is unimportant for high m values. The multiple holes aligned along the major tensile axis had the effect of decreasing the growth of the voids and, therefore, the ductility of the materials. The normalized area of voids increased by decreasing the m parameter. Generally by increasing m value, high ductility increase was observed. The biaxial deformation involved an increase in the void growth, compared with the uniaxial deformation, by increasing the true strain. Besides, in the biaxial deformation, the void growth increases by decreasing the m and n values. The void growth rate parameter values, g, shifted upward by decrease in m parameter. For all stress-strain curves at different temperatures, at lower strain after high initial stress, the stress values started to fall gradually with increase in strains, and at higher strains, the stress value achieved a steady-state stress.