Microsoft Word - numero_38_art_18 A. Znaidi et alii, Frattura ed Integrità Strutturale, 38 (2016) 135-140; DOI: 10.3221/IGF-ESIS.38.18 135 Focussed on Multiaxial Fatigue and Fracture Theorical study on mechanical properties of AZ31B Magnesium alloy Sheets under multiaxial loading A. Znaidi, O.Daghfas, S.Guellouz, R.Nasri LR-MAI-ENIT BP37, National School of Engineers of Tunis, Université de Tunis Manar 1002, Tunisie amna.znaidi@laposte.net, daghfasolfa@yahoo.fr, Sameh.guellouz@laposte.net, Rachid.nasri@enit.rnu.tn ABSTRACT. Numerical simulation by plastic deformation of the shaping processes currently has a large industrial interest. It allows you to shorten the time of design and construction related products and tools to analyze and to optimize processes. An essential part of simulation tools is the constitutive law used to describe the material used. The activity of characterization and modeling of material behavior of the plastic deformation shaping remains a very important research field of activity; the objective of proposing laws of behavior used in computer codes, essentially based on finite element is sufficiently to represent the real behavior of materials. Considering the nature of the materials used and the stresses they experience the behavior laws account for several requirements which make them increasingly complicated. Among these requirements, we cite in particular plastic anisotropy, the great transformations, the complexity and diversity of loads, etc. The complexity of these laws makes them more difficult to implement and in particular to identify: the classic tests are no longer sufficient for identification. The objective of this work is based on two essential points: Suggest a construction strategy, particularly of identifying laws elastoplastic behavior anisotropic operational for the numerical simulation of plastic deformation shaping processes with particular attention to sheet metal magnesium. Magnesium sheet metal manufacturing process involves rolling operation. In a cost-cutting goal, this operation now takes place cold, implying a very marked anisotropy of the material at the output of the mill. KEYWORDS. Mechanical behavior; Tensile test; Loading direction; Lankford coefficient; Identification. Citation: Znaidi, A., Daghfas, O., Guellouz, S., Nasri, R., Theorical study on mechanical properties of AZ31B Magnesium alloy sheets under multiaxial loading, Frattura ed Integrità Strutturale, 38 (2016) 135-140. Received: 30.05.2016 Accepted: 25.06.2016 Published: 01.10.2016 Copyright: © 2016 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. A. Znaidi et alii, Frattura ed Integrità Strutturale, 38 (2016) 135-140; DOI: 10.3221/IGF-ESIS.38.18 136 INTRODUCTION large majority of the metal parts are obtained by forming during which the material is plastically deformed. These forming processes are optimized to reduce cost; this requires manufacturers to increasingly use numerical simulation and therefore need to describe the material behavior. These simulations are often flawed by a simplified description of the plastic behavior of the material; particularly the anisotropy of rolled sheets [1-3]. Therefore, it is important to accurately model the plastic behavior of metals in large deformation in order to better predict the behavior of the part during the forming. The formulation of the anisotropic elasto-plastic behavior in large deformations is well understood now: using the formalism of rotating frame ensures the objectivity of the behavior law regardless of the constitutive model functions, [2-5]. To describe the plastic behavior of the material, it is necessary to clarify two concepts: (i) a load surface related to a plasticity criterion [6] that indicates the conditions of plastic flow, (ii) an anisotropy evolution. The experimental determination of these areas through various mechanical testing and mathematical modeling has been the subject of many research efforts such as using the Von Mises criterion: due to its implementation in most commercial finite element analysis codes. This criterion is called energy criterion in which the elastic deformation energy of the material must not exceed a limit value in order to remain within the elastic range. In the case sheet metal, the material is treated as having orthotropic plasticity where retained three preferred directions used in the expression of the Hill criterion. Also, in order to describe the asymmetric behavior in tension and compression such as the anisotropy of a structure of a sheet metal Cazacu, Barlat and al proposed in a new orthotropic criterion [7-11]. The objective of this work is to provide a model for the numerical simulation of forming processes by plastic deformation of thin metal sheets. Hence, the importance of developing a general framework for elasto-plastic orthotropic models (initial orthotropic and isotropic hardening) based on the choice of an equivalent stress, a hardening law and a plastic potential [12-14] and a model identification strategy using and experimental database [15]. This database consists of various hardening curves from various tests interpreted as homogeneous [2] and their Lankford coefficients. These plates are obtained from a hot-rolling process. At this point, the identification of constitutive parameters of the material behavior laws is an important step. A new identification strategy with its validation using Barlat's criterion will be presented. ANISOTROPIC ELASTO-PLASTIC CONSTITUTIVE LAWS n this work, we are interested in a plastic hardening behavior. The materials are treated as incompressible with negligible elastic deformations. The plastic hardening constitutive laws that we have to study fall within the following framework:  f ( , ) 0,    σ Q (1) with Q the transformation tensor between the initial time t0 and the current time t. α Represents the internal hardening variable. h( , )   σD (2) . 1( , )    σ (3) With λ plastic multiplier that can be determined from the consistency condition . f and D is the plastic strain rate tensor. This work is limited to plastic orthotropic behavior. Models are formulated for standard generalized materials with an isotropic hardening described by an internal hardening variable, a law of evolution and an equivalent deformation. The material is initially orthotropic and remains orthotropic; isotropic hardening is assumed to be captured by a single scalar internal hardening variable denoted by . In particular, we will assume that the elastic range evolves homothetically, the yield criterion is then written as follows: A I A. Znaidi et alii, Frattura ed Integrità Strutturale, 38 (2016) 135-140; DOI: 10.3221/IGF-ESIS.38.18 137 D D c sf ( , ) ( ) ( )    σ σ (4) Dσ is the deviator of the Cauchy stress tensor (incompressible plasticity). Using the special setup of the space deviators, the general form of the equivalent orthotropic plan stress, is thus: D D f1 2 3c c( ) ( X , X , X ) / ( , 2 )    σ σ (5) With D D D 1 2 3X cos ; X sin cos2 ; X sin sin2      σ σ σ Any type of criterion (4) can be written in the form: D sf ( , 2 ) / ( )    σ (6) Where  is the angle that defines the test and  the off-axis angle [12-14]. Test Expansions Equibiaxes (E.E) Simple Traction (S.T) Large Traction (L.T)     Table 1: Values of θ relative to the various tests IDENTIFICATION PROCEDURES n this section we focus on the phenomenology of plastic behavior; especially modeling plasticity and hardening based on experimental data represented as families of hardening curves, and Lankford coefficient data. In order to simplify our identification process, the following assumptions are adopted: Identification through “small perturbations” process, the tests used are treated as homogeneous tests, we neglect the elastic deformation; the behavior is considered rigid plastic incompressible, the plasticity surface evolves homothetically (isotropic hardening)and all tests are performed in the plane of the sheet resulting in a plane stress condition. The identified model is defined by an equivalent stress σc (A: σD) when σc is an isotropic function; it is assumed that the shape is defined by coefficient of the form m, A is the 4th order orthotropic tensor defined by anisotropy coefficients f, g, h, n’ and the hardening curve s (). Knowing that hardening’s curve ( )   is determined from experimental tests as r ( ). We can also begin our identification procedures using the Lankford coefficient; as determined from the tensile test by: . . . . 2 3 1 2r / -1/ (1 / )       (7) We can notice that the Lankford coefficient is independent of . This coefficient is equal to one in the case of isotropy, and remains constant in the case of transverse isotropy. However, in the case of orthotropy, it varies depending on the off-axis angle . This coefficient completely characterizes the anisotropy of the sheet when loaded in its plane. RESULTS AND DISCUSSIONS his identification strategy requires:  Experimental database.  Criterion for anisotropic plasticity. I T A. Znaidi et alii, Frattura ed Integrità Strutturale, 38 (2016) 135-140; DOI: 10.3221/IGF-ESIS.38.18 138  Validation strategy In the particular case of Magnesium sheets where anisotropy is present, the identification of this constitutive law requires the identification of the hardening function, the anisotropy coefficients, the form factor m and r ( ) the Lankford coefficient. Using as a hardening function a Hollomon law: n s k   (8) And as a plasticity criterion using the Barlat model [6]: m mmm c q q q q q q1 2 2 3 1 3        (9) Where q1, q2 and q3 are the eigenvalues of the tensor q c s ai i(q) - ( ) 0. Where (a ) : thus (a ) :      q A q Aσ (10) Using the simplex algorithm and using the plastic Barlat model and respecting the assumptions, the identification of the thin Magnesium sheet is equivalent to choosing the model coefficients while minimizing the squared difference between the theoretical and experimental results. ψ k n 00° 408.5 0.102 22.5° 410.6 0.096 45° 425.6 0.101 67.5° 432.5 0.0901 90° 437.1 0.0873 Table 2: Identification of the constants of hardening law for different tractions tests. Knowing that the coefficient n is the same for all tests as demonstrated at the beginning of this work, by convention we choose n for traction in direction ψ = 00° as reference. For n=0.102, we present different values of k (see Tab. 3). ψ k 00° 418.1 22.5° 418.1 45° 427.1 67.5° 447.3 90° 455 Table 3: Identification of the constant hardening law for fixed n. In Fig.1, the experimental hardening curve (exp) and the curve identified from a model (iden) using an average value of n are represented. For tensile tests, the models (iden) give a clear fit between the theoretical and experimental results. Our second identification step amounts to determining the coefficients of anisotropy (f, g, h, n’) and the shape coefficient m (Tab. 4), considering the non-quadratic Barlat's criterion (9). A. Znaidi et alii, Frattura ed Integrità Strutturale, 38 (2016) 135-140; DOI: 10.3221/IGF-ESIS.38.18 139 Figure 1: Identification of the hardening curve:  =00°. f g h n’ m 0.137 0.0992 1.4889 3.104 1.1999 Table 4: Identification of anisotropic coefficients and a shape coefficient m In Fig. 2 the model (ident) allows us to study the load surface on each test. We note that this material is resistant to simple shear like simple traction. Figure 2: Evolution of the load surface in the deviatory plan  x x2 3, . Figure 3: Evolution of the anisotropy. 0,5 0,7 0,9 1,1 0 0,4 0,8 1,2 1,6       A. Znaidi et alii, Frattura ed Integrità Strutturale, 38 (2016) 135-140; DOI: 10.3221/IGF-ESIS.38.18 140 Using the identified anisotropic coefficients, we represent in Fig. 3 the evolution of the anisotropy based on off-axis angles; using the Barlat's criterion this material is very anisotropic at  =45°. CONCLUSION In this work we show that identification strategy results can be extracted. This identification has focused both on plastic material parameters of the constitutive law and Lankford coefficient. Thus, the plastic behavior model: Hollomon Law and Barlat criterion with 5 parameters are identified. A validation by comparing model / to experiment data was performed.The model using the Barlat's criterion is in good agreement with experimental results relating to Lankford coefficients. Following this strategy, we observed very pronounced anisotropy of AZ31BMagnesium and the load surface for different tests at the end of this identification. With this strategy, we can study in a more precise way the anisotropy of the Magnesium model by integrating Barlat in integrating the model Hill. REFERENCES [1] Ghouati, O., Gelin, J.C., A finite element based identification method for complex metallic material behavior, Computational Materials Science, 68 (2001) 2157. [2] Znaidi, A., Plasticity orthotrope in large deformation, PhD thesis at the Faculty of Sciences of Tunis (2004). [3] Znaidi, A., Gahbiche A., Identification of the plastic behavior of orthotropic thin plates based on non-quadratic yield criterion. 6th National Conference of Physics Research Hammamet (1999). [4] Boubakar, L., Boisse, P., Anisotropic elastoplastic behavior for numerical analysis of thin shells in great transformations. Européene review Mechanics, 7(6) (1998). [5] Manget, B., Perre, P., A large displacement formulation for anisotropic constitutive laws, European journal of Mecanique, 18(5) (1999). [6] Barlat, F., Brem, D. L. J., A six components yield function for anisotropic materials, Int. J. Plasticity, 7 (1991) 693- 712. [7] Gronostajsk, Z., The constitutive equations for fem analysis. Journal of material processing technology, 106 (2000) 40-44. [8] Hill, R., A theory of the yielding and plastic flow of anisotropic metals, Proc. Roy. Soc of London, A93 (1948) 281- 297. [9] Cazacu, O., Ionescu, I. R., Yoon, J. W., Orthotropic strain rate potential for the description of anisotropy in tension and compression of metals, Int. J. Plasticity, (2009). [10] Cazacu, O., Plunkett, B., Barlat, F., Orthotropic yield criterion description of anisotropy in tension and compression of sheet metals, Int. J. Plasticity, 24 (2008) 847-866. [11] Plunkett, B., Cazacu, O., Barlat, F., Orthotropic yield criterion for hexagonal closed packed metals, Int. J. Plasticity, 22 (2006) 1171-1194. [12] Daghfas, O., Znaidi, A., Nasri, R., Anisotropic behavior of mild steel subjected to isotropic and kinematic hardening, In: 6th International Conference on Advances in Mechanical Engineering and Mechanics (ICAMEM2015), Hammamet, (Tunisia) (2015) 128. [13] Baganna, M., Znaïdi, A., Kharroubi, H., Nasri, R., Identification of anisotropic plastic behavior laws for aluminum 2024 after heat treatment materials from off-axis testing, Toulouse, (2010). [14] Znaidi, A., Daghfas, O., Toussaint, F., Nasri, R., Strategy for the identification of anisotropic behavior laws for Ti40 sheets, AMPT2015 Madrid, (2015). [15] Yueqian, J., Xiang, L., Yuanli, B., Theorical Study On Mechanical Properties of AZ31B Magnesium Alloy Sheets Under Multiaxial Loading, J. Automotive Safety and Energy, 4 (2012). << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice