Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 55, 3, pp. 923-935, Warsaw 2017 DOI: 10.15632/jtam-pl.55.3.923 STUDY OF THE INFLUENCE OF COLD WORKING ON MECHANICAL BEHAVIOR AND DUCTILE FRACTURE OF 5754 ALUMINUM ALLOY: EXPERIMENTAL AND NUMERICAL SIMULATIONS Wafa Taktak Laboratoire des Systemes Electro-mecaniques (LASEM), National School of Engineers of Sfax, Tunisia e-mail: wafa.taktak@yahoo.fr Rim Taktak, Nadder Haddar Laboratoire Genie des Materiaux et Environnement (LGME), National School of Engineers of Sfax, Tunisia Riadh Elleuch Laboratoire des Systemes Electro-mecaniques (LASEM), National School of Engineers of Sfax, Tunisia The ductile damage of automotive aluminum sheet alloy AA5754-H111 is investigated by experiments and numerical simulation using the Gurson-Tvergaard-Needleman (GTN)mo- del. The GTN parameters were determined by a uni-axial tensile test and the inverse finite elementmethod. The same parameters were employed to provide the ductile damage beha- vior of central crackedpanel (CCP) specimens. A good prediction can be established among the numerical simulation and experimental data in from of the force opening displacement. As an application, the identified GTN model is used to predict the influence of cold wor- king on deformation and ductile damage. The numerical simulation results obtained are assimilated with experimental data. Keywords: GTNmodel, ductile fracture, AA5754-H111, cold working 1. Introduction Inautomotive applications, thealuminumalloys are extensively used for obtaining lightmassand high strength structures. Aluminum-Magnesium (Al-Mg) aluminum alloys, indicated by 5xxx series, have a very good formability but a relativity low strength. This series of alloys strengthen only bywork hardening (Burger et al., 1995). The results of investigation of mechanical damage of 5754-H111 aluminum alloy has indicated that damage and fracture aremostly results of nuc- leating, growing and coalescing of micro cavities or micro voids. In order to predict the ductile fracture process, several theoretical models using local approaches have been presented in the literature (Betegon et al., 1997; Corigliano et al., 2000; Imad et al., 2003; Achouri et al., 2012). So, the selection of an adaptedmicromechanical model allows understanding of the fractureme- chanismof 5457-H111 aluminumalloy. Thebasic researchwas started byMcClintock (1968) and afterwards by Rice and Tracey (1969) who investigated the growth of cylindrical and spherical voids in ductile solids. Founded on theses analysis, Gurson (1977) proposed a micromechanical approachmodel. Later, theGursonmodel wasmodified byTvergaard andNeedleman (1984) by introducing parameters q1 and q2. They founded the void fusion equation f ∗ to describe ductile failure by nucleation, growth and coalescence of spherical micro voids. Generally, the modified Gursonmodel is named theGTNmodel.Many researchers have used theGTNmodel to provide ductile porousmaterials. Benseddiq and Imad (2007) used theGTNdamagemodel to investiga- te ductile tearing of 2024-T351 aluminum alloy. Oh et al. (2007) proposed a phenomenological model of ductile fracture for API X65 using the GTN model. Yan et al. (2013) employed the GTN damage model to study the initiation and propagation of the crack near the edge under 924 W. Taktak et al. rolling condition.Guo et al. (2013), Zhang et al. (2000), Huang et al. (2007) andHu et al. (2014) used numerical results and experimental data to obtain parts of theGTNparameters. Unknown parameters of the model were obtained by using the inverse finite element method. In thepresent research, ductile tearingof 5754aluminumalloyhasbeenanalyzedbyusing the GTNmodel.TheGTNparameters are identifiedbycombininguni-axial tensile tests anddetailed finite element analyses. The identified parameters have been employed to predict ductile failure of central cracked panel specimens.Agood agreements can be established between the numerical and experimental results in form of the force versus displacement curves. As an application, the calibrated model is used to measure the cold working influence on deformation and fracture of 5754 aluminum alloy. The theoretical results are compared with experimental ones in the case of cold worked tensile and central cracked panels (CCP). 2. The GTN model Gurson (1977) proposed a model to describe damage of ductile materials based on a micro- mechanical approach to porous plastic solids. This model takes into account degradation of the load carrying capacity due to the presence of porosity in isotropicmaterials. The originalGurson model wasmodified byTvergaard andNeedleman (1984) who relied on the coalescence of voids. The yield surface is presented by Φ(σy,σeq,f)= σ2eq σ2y +2f∗q1cosh (3 2 q2 σm σy ) − [1+q3(f∗)2] = 0 (2.1) where q1, q2 and q3 = (q1) 2 are the constitutive parameters incorporated by Tvergaard (1981) to amplify the hydrostatic stress effect for all strain levels. σeq is the conventional von Mises equivalent stress definedby:σeq = √ 3 2 SijSij (Sij is the stress deviator) andσy is the yield stress of the undamagedmatrix material. For good prediction of the effect of void coalescence, Tvergaard and Needlemen (1984) in- troduced the damage function f∗ defined as follows f∗(f)=        f for f σy (3.1) with the hardening exponent n and ductility coefficient k. Themechanical proprieties of 5457-H111 aluminum alloy are summarized in Table 2. As shown in Fig. 3, the true stress-strain results of 5457-H111 aluminum alloy are similar in the three directions. 926 W. Taktak et al. Fig. 2. (a) Geometry of the tensile specimen (in mm), (b) dimensions of the CCP specimens (in mm) Table 2.Mechanical proprieties of 5457-H111 aluminum alloy in three directions:E – Young’s modulus, σy – yield stress, σu – ultimate stress,A% – elongation, n – hardening exponent and k – ductility coefficient E [MPa] σy [MPa] σu [MPa] δ [%] n K [MPa] RD 70612 100 270 15.28 0.586 575.554 TD 70134 99 265 14.89 0.589 599.659 DD 69978 98 259 15.13 0.580 590.034 Fig. 3. Strain-stress curves of 5754-H111 aluminum alloy in the rolling direction (RD), transverse direction (TD) and diagonal direction (DD) Ductile tearing tests were carried out on central cracked panels (CCP) at a constant displa- cement rate of 1mm/min. The specimen geometry is presented in Fig. 2b. The ductile tearing specimens were manufactured in the rolling direction (RD) and transverse direction (TD) of the plates. In order to obtain a normalized crack length ratio a/w equal to 0.36, the specimens Study of the influence of cold working on mechanical behavior and... 927 were pre-cracked by the fatigue test. These tests were achieved using the single specimen me- thod (Taktak et al., 2009). Three specimens were operated at a constant cross-head speed of 1mm/min. The experimental load versus displacement curves of 5754-H111 aluminumalloy are exposed in Fig. 4. It is noted that these curves show a low experimental dispersion. Fig. 4. Load versus displacement curves of 5754 aluminum alloy in the rolling direction (RD) and transverse direction (TD) Also the fractographic examinations of 5754 aluminum alloy were performed on the fracture surfaces of the broken specimens from tensile testing by employing a scanning electron micro- scope (SEM). The SEM fractographs for 5754 aluminum alloy indicate the significant presence of ductile dimples (large voids beside smaller voids), which show the characteristic micro-void coalescence mechanism of ductile fracture (Fig. 5a-5d). Fig. 5. SEM fractographs of the facture surface for a tensile specimen abstracted from 5754 aluminum alloys under different magnifications: a overviewmorphology; b extended image of indicated region of the picture (a), c extended image of the indicated region of the picture (b), and d extracted image of the indicated region of the picture (c) 928 W. Taktak et al. 4. Numerical analysis 4.1. The meshing and boundary conditions To obtain the theoretical force versus displacement curves, three-dimensional computation was performed using the finite element programANSYS. In this study, thematerial is supposed to be elastoplastic, homogeneous and exhibit isotropic hardening. Two forms have been studied: the tensile specimen (Fig. 1a) and the CCP specimen (Fig. 1b). Themeshing corresponding to each specimen using SOLID 185 finite element of ANSYS is presented in Fig. 6. The meshing near the crack tip consists of squaremesh with 0.2mm size (Taktak et al., 2008). All the nodes on the ligament in front of the crack tip are fully constrained to have a zero displacement in the y-direction normal to the plane of the crack. Due to symmetry, only a quarter of the tensile specimen (Fig. 6a) and CCP specimen (Fig. 6b) are modeled. Fig. 6. (a)Meshing of a quarter of tensile specimen, (b) meshing of a quarter of the CCP specimen 4.2. Identification of damage parameters The (GTN) ductile damage model has 9 necessary parameters: Three constitutive parame- ters related to the refined yield locus, q1, q2 and q3, three void nucleation parameters εn, Sn and fn and three parameters for void growth and coalescence f0, fc and fF . In the present inve- stigation, seven parameters are determined by typical values proposed in literature (Tvergaard andNeedleman, 1984; Lievers et al., 2004): εn =0.65,Sn =0.03, fn =0.00035, q1 =1.5, q2 =1, q3 = q 2 1 =2.25 and ff =0.25. In accordance with the discussion in (Zhang et al., 2000; Rousse- lier, 2001), the critical void volume fraction fc can be chosen as 0.15 for aluminum alloys. Thus, only the initial volume fraction f0 needs to be determined. Generally, in aluminum alloys, the void is made up of a fragile intermetallic phase (Ghahremaninezhad and Ravi-Chandar, 2012). The evaluation of the initial void volume fraction f0 is obtained by metallorgraphic examina- tion on a polished surface of undamagedmaterials. Also, f0 is very little for this kinds of alloy. The initial void volume fraction f0 can be estimated by the inverse finite element method using the finite element modeling and experiments. Guo et al. (2013) identified this parameter using the original Rousselier model to obtain three values of f0 0.0001, 0.001, 0.005, respectively. By comparisonwith experimental results, he found a good agreement with f0 =0.0001 forAA5052. Hu et al. (2014) identified the parameter f0 by an inverse method using experimental results. The simulations were done by setting the value of f0 from 0.001 to 0.03. The initial void volume fraction f0 for AA 6016 equal 0.001 improved the good agreement between the simulation and experimental results. In thiswork, to estimate theparameter f0, an inversemethod is usedbasedon theFEmethod presented in this Section and the experimental results presented in Section 2. Six analyses have been made using the GTN model and the following initial void volume fraction f0: 0.00001, 0.0001, 0.001, 0.01, 0.05 and 0.1, respectively. Study of the influence of cold working on mechanical behavior and... 929 Fig. 7. The true stress-strain curves with different initial void volume fractions f0 obtained by FE using the GTNmodel The true stress-strain curvesmadeby simulationswithdifferent f0 are shown inFig. 7. It can also be shown that, for f0 =0.001, the difference between the simulation and experiment results is admissible. So, the initial void volume fraction f0 was identified as 0.001 for this material. The whole GTN parameters of the material are now presented in Table 3. Table 3.Calibrated parameters of the GTNmodel for the 5754-H111 aluminum alloy q1 q2 q3 f0 fn εn Sn fc ff 1.5 1 2.25 0.001 0.00035 0.65 0.03 0.15 0.25 To validate the identified parameters listed in Table 3, the experimental load versus displa- cement curve issued from tearing tests on the CCP specimen are compared with the simulated ones using these parameters in Fig. 8. A good agreement is observed between the experimental results and those prediced by the finite element method. This shows that the GTN parameters identified in this work are acceptable. 5. Application to cold working effects In the preceding Sections, the ductile fracture model based on the GTN model has been found for 5754-H111 aluminum alloy. This model has many potential application domains: it can be used for predicting not only the failure behavior of ductile tearing tests specimens but also for predicting size effects of ductile tearing tests specimens.Theapplication of 5754-H111 aluminum alloy in automotive industry needs a high strength/weight ratio (Burger et al., 1995). So, it is possible to carry out the strain hardening in the rolling direction. In the literature, many experimental investigations have beenpublishedonquantification of the strain hardening impact onmechanical coldworking fractureproperties for aluminumalloys.Our essential interest in this Section is the application and validation of the identified GTNmodel to measure cold working effects on numerical simulation. Tests were carried tomeasure the coldworking effects on tensile andductile tearingproprieties of 5754-H111 aluminumalloy. In these tests, theusedmaterialwas 930 W. Taktak et al. Fig. 8. Load versus displacement. Comparison between the experiment and FE simulation for the GTNmodel similar condition to the one used in the previous Sections. A plate of the alloy (with thickness of 2.45mm,width 150mmand length of 150mm)was cold rolled by a laboratory rollingmill to the reduction of 25% and 50% in area (CW=25% and CW=50%). Both the tensile and CCP specimens were abstracted from the plate in the longitudinal direction. The experimental conditions and geometry of the tensioned specimen and ductile te- aring tests were same as those in the previous Sections. The tensile properties of the material subjected to different cold working rates are summarized in Table 4. It shows that the yield and tensile strength increase with the increasing percent of cold working, but the ductility decre- ases (Cosham, 2001; Mansourinejad and Mirzakhani, 2012). The strengthening of the material can be described by the increase of dislocation density with plastic deformation. The average distance between dislocations decreases and the dislocations start blockingmotion of each other (Hajizadeh et al., 2014). It is widely acknowledged that the cold working effect on tensile curves can be quantified simply by shifting the true strain-stress curve by an increase of cold working (Ainsworth, 1986; Cosham, 2001), which is carried out by analyzing the present experimental results, as shown in Fig. 9. Table 4.Mainmechanical characteristics for different coldworking rates:E –Young’smodulus, σy – yield stress, σu – ultimate stress,A%– elongation CW E [MPa] σy [MPa] σu [MPa] A [%] 0% 70612 100 270 15.28 25% 71334 168 278 7.50 50% 69789 190 307 5.29 As shown in Fig. 10, the values of the hardening exponent n and ductility coefficient k are similar in the two directions. For this reason, we can consider that the behavior of these work-hardenedmaterials is isotropic. Study of the influence of cold working on mechanical behavior and... 931 Fig. 9. Experimental results of the cold working effect on tensile curves Fig. 10. Values of the hardening exponent n and ductility coefficient k in the rolling direction (RD), transverse direction (TD) The SEMpictures of the cold worked specimens indicated a fine network of dimples (elonga- ted voids with the fibrous structure) and little quasi-cleavage parts corresponding to the ductile fragile failuremechanismwhen the percentage of cold working increases (shown in Fig. 11). Ge- nerally, inFCCmetals like 5457-H111 aluminumalloys, even at low temperature, the dislocation of leavings is important and the material rests ductile enough. Thus, the morphology of fracture ought to be fundamentally fibrous plus some cleavage, reflecting the fragile failure behavior of the alloys when the percentage of cold working high. The experimental load-displacement curves are summarized in Fig. 12. It indicates that the maximum loads for cold working are lower than those for without cold working. Using the overhead information, finite element damage analyses based on the GTN model are used for simulating tensile tests and ductile tearing using ANSYS. To incorporate cold working into finite element damage analyses, twomodifications are added. First, the true stress- -strain curve ismodified according to coldworking. Secondly, the value of the initial void volume fraction f0 is changed (Oh et al., 2007). As noted, the value of the initial void volume fraction f0 increases when the cold working rate increases (Oh et al., 2007). 932 W. Taktak et al. Fig. 11. SEM fractographs of 5457 aluminum alloy at different percentage of cold working: (a) 0%, (b) 25% and (c) 50% Fig. 12. Experimental results of cold working effects on load versus displacement curves Figure 13 compares the experimental true stress-strain results with the numerical ones from finite element damage analyses based on the GTNmodel for different choices of f0 for selected cases (CW=25% and 50%). From Figs. 13a and 13b, it can be seen that the choice of f0 =0.006 and f0 =0.009 for the cold working rate 25% and 50%, respectively, enables achieving the best agreement between the (FE) damage analyses and experimental results. For validation, the results from finite element damage analyses for cold working rates 25% and 50% are compared with the experimental load-displacement curves, and reasonably good comparisons are found. It is shown in Fig. 14. Study of the influence of cold working on mechanical behavior and... 933 Fig. 13. True stress-strain curves with different initial void volume fractions f0 obtained by FE using the GTNmodel. (a) The initial void volume fraction f0 is calibrated as 0.006 for CW=25%. (b) The initial void volume fraction f0 is calibrated as 0.009 for CW=50% Fig. 14. Comparisons of FE simulations with experiments for the load versus displacement curves (a) CW=25%, (b) CW=50% The results show that the experimental tensile tests and load versus displacement have been successfully predicted by theGTNmodel, and can take into account the cold working effect not only on plastic deformation but also on ductile fracture. 6. Summary and conclusions • The continuum damage mechanics model (GTN) has been used to simulate the ductile tearing behavior of 5754-H111 aluminum alloy sheet metal. • The parameters of the GTN model have been identified by an experimental tensile test (true stress versus true strain) and the inverse finite element method. • The validity of the proposedparameters has been investigated by comparing the simulated results with the experimental ones from the tensile and ductile tearing tests (load versus displacement). • The determined GTN model has been applied to predict the cold working influence on deformation and fracture. Comparison of experimental data of cold working, tensile tests andductile tearing testswithfinite element damage analyses have showngood agreements. 934 W. Taktak et al. References 1. Achouri M., Germain G., Dal Santo P., Saidane D., 2012, Implementation and validation of a Gurson damagemodel modified for shear loading: effect of void growth rate andmesh size on the predicted behavior,Key Engineering Materials, 504, 691-696 2. AinsworthR.A., 1986,Anassessmentof the effects of prestrainonupper shelf fracture toughness, The Journal of Strain Analysis for Engineering Design, 21, 219-224 3. Benseddiq N., Imad A., 2007, A ductile fracture analysis using a local damagemodel, Interna- tional Journal of Pressure Vessels and Piping, 85, 219-227 4. Betegon C., Rodriguez C., Belzunce F.J., 1997, Analysis and modelisation of short crack growth by ductile fracture micromechanisms, Fatigue and Fracture of Engineering Material and Structures, 20, 633-644 5. Burger G.B., Gupta A.K., Jeffrey P.W., Lloyd D.J., 1995, Microstructural control of aluminum sheet used in automotive applications,Materials Characterization, 35, 23-39 6. Cosham A., 2001, A model of pre-strain effects on fracture toughness, Journal of Offshore Me- chanics and Arctic Engineering, 123, 182-190 7. Corigliano A., Mariani S., Orsatti B., 2000, Identification of Gurson-Tvergaard material model parameters via Kalman filtering technique. I. Theory, International Journal of Fracture, 104, 349-373 8. GhahremaninezhadA., Ravi-ChandarK., 2012,Ductile failure behavior of polycrystallineAl 6061-T6, International Journal of Fracture, 174, 177-202 9. Guo J., Zhao S.,Murakami R.I., Zang S., 2013,Experimental and numerical investigation for ductile fractureofAl-alloy 5052usingmodifiedRousseliermodel,ComputationalMaterials Science, 71, 115-123 10. Gurson A.L., 1977, Continuum theory of ductile rupture by void nucleation and growth: Part I – Yield criteria and flow rules for porous ductile media, Journal of Engineering Materials and Technology, 99, 2-15 11. Hajizadeh K., Tajall M., Emadoddin E., Borhani E., 2014, Study of texture, anisotropy and formability of cartridge brass sheets, Journal of Alloys and Compounds, 588, 690-696 12. Hu Y.M., Chen M.Z., Xiao Y., Xiao J., Tan X., Tang Q., Cui T.S., 2014, Parameters determination of GTN model and damage analysis of aluminum alloy 6016 sheet, International Conference on Material Science and Applications (ICMSA-15), Atlantis Press 13. Huang X.-C., Chen Y.-Z., Chen Y.-M., Niu W., 2007, Identification of mesodamage parame- ters of 2169 steel, Journal of Materials Engineering, 4, 50-52 14. Imad A., Wilsius J., Abdelaziz M.N., Mesmacque G., 2003, Experiments and numerical approaches toductile tearing in an2024-T351aluminiumalloy, International Journal ofMechanical Sciences, 45, 1849-1861 15. Lievers W.B., Pilkey A.K., Lloyd D.J., 2004, Using incremental forming to calibrate a void nucleationmodel for automotive aluminum sheet alloys,Acta Materialia, 52, 3001-3007 16. Mansourinejad M., Mirzakhani B., 2012, Influence of sequence of cold working and aging treatment on mechanical behaviour of 6061 aluminum alloy, Transactions of Nonferrous Metals Society of China, 22, 2072-2079 17. McClintockF.A., 1968,Acriterion for ductile fractureby thegrowthofholes,Journal of Applied Mechanics, 35, 363-371 18. Oh C.K., Kim Y.J., Baek J.H., Kim Y.P., 2007, A phenomenologicalmodel of ductile fracture for API X65 steel, International Journal of Mechanical Sciences, 49, 1399-1412 19. Rice J.R., Tracey D.M., 1969, On the ductile enlargement of voids in triaxial stress fields, Journal of the Mechanics and Physics of Solids, 17, 201-217 Study of the influence of cold working on mechanical behavior and... 935 20. Rousselier G., 2001, The Rousselier model for porous metal plasticity and ductile fracture, Handbook of Materials Behavior Models, 2, 436-445 21. Taktak R., Benseddiq N., ImadA., 2009,Analysis of ductile tearing using a local approach to fracture,Fatigue and Fracture of Engineering Materials and Structures, 32, 525-530 22. Tvergaard V., 1981, Influence of voids on shear band instabilities under plane strain conditions, International Journal of Fracture, 17, 389-407 23. Tvergaard V., Needleman A., 1984, Analysis of the cup-cone fracture in a round tensile bar, Acta Metallurgica, 32, 157-169 24. YanY., SunQ., Chen J., PanH., 2013,The initiation and propagation of edge cracks of silicon steel during tandemcold rollingprocess basedon theGurson-Tvergaard-Needlemandamagemodel, Journal of Materials Processing Technology, 213, 598-605 25. Zhang Z.L., Thaulow C., Odegard, 2000, A complete Gurson model approach for ductile fracture,Engineering Fracture Mechanics, 67, 155-168 Manuscript received October 14, 2016; accepted for print March 31, 2017