Stress-state dependency of ductile fracture in an extruded magnesium alloy and its underlying mechanisms
Introduction
Exhibiting a combination of high specific strength, good machinability, and being environmentally friendly, magnesium alloys are promising to serve as structural products in industries where lightweight design is vital (Nie et al., 2020). However, magnesium alloys have relatively poor formability at room temperature due to their hexagonal close-packed (HCP) crystal structure and the consequently limited number of active slip systems. The ductile fracture prone to occur during deformation has become a big challenge for engineering applications of magnesium alloy. In this context, a comprehensive study on the deformation and the ductile fracture of magnesium alloys at various stress states is indispensable. This demands both elaborately tailor-designed experiments and effective modeling strategies.
It is well established that fracture behaviors of metallic materials are dictated by the stress state, which can be uniquely determined by stress triaxiality and Lode parameter. Investigations on the stress-state dependency of ductile fracture in magnesium alloy sheets were reported, e.g., Abedini et al. (2018) studied the fracture behaviors of ZEK100 sheet in proportional and non-proportional cases; Habib et al. (2019) investigated the fracture behaviors of ZEK100 sheet at different stress states and strain rates; Jia and Bai (2016) conducted a research on the fracture behaviors of AZ31B sheet under various loading conditions; Lee et al. (2018) examined the fracture limits of AZ31 and ZE10 sheets under various stress states. In these reports, notched plate specimens were frequently used for tensile fracture tests. As the notch acuity changes, local stress states at the specimen center vary between uniaxial tension and plane strain tension. Although widely used, these tests are restricted to plane stress states.
To study the fracture behaviors at triaxial stress states, notched round specimens were almost exclusively used. Specific to magnesium alloys, tension experiments of the smooth and differently notched round specimens were conducted to explore the triaxiality effect on room-temperature ductility of hot-rolled AZ31 (Kondori and Benzerga, 2014) and WE43 (Kondori and Benzerga, 2015) plates. It was shown that as stress triaxiality increases, fracture strains of these two magnesium alloys first increase to a maximum and then decrease in varying degrees. The authors (Kondori and Benzerga, 2014, 2015) attributed this trend to fracture mode transition associated with the activating deformation mechanisms, second phase particles and grain-boundary precipitates, etc. Motivated by the experiments of Kondori and Benzerga (2014), Selvarajou et al. (2016, 2017) quantified the distributions of stress triaxiality and the evolutions of relative activities of key deformation mechanisms in smooth and notched round specimens through crystal plasticity FE simulations. Although modeling of the mechanism-based failure was not attempted in their work, some fundamental understanding of ductile fracture in magnesium alloy was inspired. One advantage of using the notched round specimens for tension tests is that stress triaxialities generated at the notch area increase as the notch gets sharper. This renders notched round specimen ideal for the investigation of ductile fracture at high stress triaxialities. Another advantage is that plastic deformations are localized at the notch area, and ductile fracture can be expected to initiate there. Nevertheless, the exact site of fracture initiation is somewhat ambiguous. Typically, ductile fracture initiates from the inner center of the notched round specimen in tension tests of aluminum alloys (Bao and Wierzbicki, 2004) and steels (Benzerga et al., 2004a, b), while in tension tests of AZ31 magnesium alloy, fracture initiation was reported to occur near the notch root rather than the inner center when the notch is severe (Kondori, 2015, 2018).
The study reviewed above mainly concerned the ductile fracture at positive stress triaxialities. By comparison, ductile fracture at negative stress triaxialities has received fewer research efforts. One important reason might be that ductile fracture in metals is generally the result of the nucleation, growth and coalescence of micro-voids (Shang et al., 2020), and a compressive state would probably impede such processes. Whereas, there is growing experimental evidence that ductile fracture can occur even at very low stress triaxialities due to the evolution of micro-cracks (Khan and Liu, 2012; Brünig et al., 2018). The fact is that microscopic fracture mechanisms vary at different stress states, and this leads to diverse macroscopic fracture behaviors (Pineau et al., 2016). For the magnesium alloy with a limited ductility at room temperature, ductile fracture at both positive and negative stress triaxialities merit serious attention.
Cylindrical specimens are conventionally used for compression tests. These tests are easy to implement but have the following limitations in the characterization of ductile fracture. On the one hand, compression results are highly sensitive to friction conditions between the specimen and platen, which is unfavorable for obtaining reliable experimental data. On the other hand, ductile fracture generally initiates at the equatorial bulging area of the cylindrical specimen, where the deformation has experienced a transition from compression to shear. This means that using cylindrical specimen is insufficient to study the ductile fracture at low negative stress triaxialities.
By modification of the cylindrical specimen, some novel specimens were proposed. For example, the cylindrical specimen with a spherical recess at the equatorial area (Kubík et al., 2016) and that with a non-prismatic notch along the length (Kubík et al., 2018) were adopted. Such modifications provide the exact fracture initiation site with an average stress triaxiality lower than −1/3. However, due to the ‘eccentric’ geometry of the specimen, obvious distortion appears during compressive loading and stress triaxiality at the fracture initiation site changes enormously as well. Additionally, specimens with only one certain geometry were adopted and no method to change or control the generated stress triaxiality was introduced. To achieve various negative stress triaxialities, some biaxially loaded specimens were proposed. Brünig et al. (2018) designed new cruciform specimen with a notch at its center and loaded it in vertical and horizontal directions concurrently. The vertical load leads to shear deformation and the horizontal load superimposes tensile or compressive deformation in the notched part. By applying a vertical and horizontal load ratio of 1: −8, the cruciform specimen fractured at a stress triaxiality low to −0.56. In the same vein, Gerke et al. (2017) designed double symmetric H-specimen with four notched regions and loaded it in two dimensions. With different loading ratios, stress triaxialities in the range of (−0.6, 0.8) were obtained. Although such biaxially loaded specimens appear to be versatile, the specimen processing is cumbersome and a special biaxial testing machine is needed. Actually, negative stress triaxialities can be realized in a notched round specimen under the compressive condition according to the Bridgman formula (Bridgman, 1952). However, unlike notched tensile specimens, notched compressive specimens are seldom reported. Bao and Wierzbicki (2004) used the notched round specimen for compression tests to study the fracture behaviors of 2024-T351 aluminum alloy. In their study, only one type of notched round specimen was adopted, and it was designed mainly to eliminate the friction effect during compression. Moreover, stress triaxiality at the fracture initiation site of the designed specimen gets closer to zero as deformation increases, which resembles the situation in compression tests of the cylindrical specimen. It seems that to date there have not been satisfactory specimens for fracture tests in the negative stress triaxiality regime.
Experiments aside, ductile fracture modeling is another important subject, since fracture tests at a lot of stress states are not achievable. In this respect, some advanced phenomenological fracture models, e.g. MMC model by Bai and Wierzbicki (2010), DF2014 model by Lou et al. (2014), and HC model by Mohr and Marcadet (2015), work well by taking both accuracy and efficiency into consideration. Previously, we proposed a new framework for establishing multi-component fracture models, and the established two-component DF2014 model was shown to have higher accuracy than the single-component counterpart when applied to describe the fracture strains in a wide range of stress triaxialities (Li and Fang, 2018). In phenomenological models, a damage index is defined as the integral of a stress-state function along strain path. Ductile fracture is assumed to occur when the damage index reaches a critical value. To obtain the accurate stress and strain fields during deformation, a hybrid experimental-numerical procedure is always exploited (Dunand and Mohr, 2010; Papasidero et al., 2014, 2015). In this procedure, the plasticity model employed in numerical simulations is of great importance for successful implementation of the fracture model.
In the present work, ductile fracture of an extruded magnesium alloy at both positive and negative stress triaxialities is investigated experimentally and numerically. Differently notched round specimens are used for tensile fracture tests, spanning a wide range of positive stress triaxialities. To realize negative stress triaxialities, new types of notched round specimens are proposed for compressive fracture tests. The new specimens are designed to have different notch radii, hence the triaxiality effect in the range of compressive states can be characterized quantitatively. What is worth mentioning is that during all fracture tests of notched round specimens, the loading process is successfully stopped at the onset of macroscopic crack, which facilitates the correct location of the fracture initiation site. In conjunction with the plasticity model based on Yoon2014 function (Yoon et al., 2014), FE simulations of all fracture tests are performed. The obtained experimental and numerical data are utilized to calibrate both the single- and two-component DF2014 fracture models. After that, a comparison of the description results between these two models is performed. Finally, the fractured surfaces of various specimens are detected and analyzed by scanning electron microscopy (SEM), to shed some light on the fracture mechanisms of magnesium alloy at different stress states.
Section snippets
Material and experimental procedure
The material to study is an extruded Mg-Al-Zn-RE alloy bar (Bai et al., 2019), with a rectangular cross-section of 50 mm × 30 mm. Plastic deformation characteristics of the magnesium alloy in three principal directions, i.e. extrusion direction (ED), transverse direction (TD), and normal direction (ND), as well as three off-axis orientations, were studied in detail previously (Li and Fang, 2020). To learn about its ductile fracture behaviors with respect to ED, various types of fracture tests
Plasticity model based on Yoon2014 function
In a hybrid experimental-numerical procedure aimed to obtain accurate stress and strain fields at any stage of deformation, establishing an appropriate constitutive model is essential. In our previously published work (Li and Fang, 2020), it has been shown that the plasticity model comprised with Yoon2014 yield function (Yoon et al., 2014), non-associated flow rule (Stoughton, 2002; Stoughton and Yoon, 2009), and yield surface interpolation method in hardening description, hereinafter referred
Load-displacement responses
The load-displacement responses of the magnesium alloy in all fracture tests obtained by experiments and simulations are compared in Fig. 4. Note that for a certain test, only one representative experimental result is presented. The displacements refer to the relative vertical movements of two reference points B and B’ tracked by DIC, instead of those recorded by the position sensor of the testing machine. The fracture initiation points are denoted by blue diamond symbols on load-displacement
Conclusion
Ductile fracture behaviors and the underlying fracture mechanisms of an extruded magnesium alloy are investigated by experiments and simulations. In particular, new notched compressive specimens are proposed for fracture tests at negative stress triaxialities. FE simulations of all fracture tests using the plasticity model based on Yoon2014 function are validated by global load-displacement responses and local surface strain evolutions. On this basis, fracture strains and corresponding loading
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
The present research is financially supported by the National Natural Science Foundation of China (No. 52075288), and the authors would like to express acknowledgement for it.
References (43)
- et al.
Experimental fracture characterisation of an anisotropic magnesium alloy sheet in proportional and non-proportional loading conditions
Int. J. Solids Struct.
(2018) - et al.
The role of deformation twinning in the fracture behavior and mechanism of basal textured magnesium alloys
Mater. Sci. Eng. A
(2014) - et al.
On fracture locus in the equivalent strain and stress triaxiality space
Int. J. Mech. Sci.
(2004) - et al.
Anisotropic ductile fracture: part I: experiments
Acta Mater
(2004) - et al.
Anisotropic ductile fracture: part II: theory
Acta Mater
(2004) - et al.
Damage and failure at negative stress triaxialities: experiments, modeling and numerical simulations
Int. J. Plast.
(2018) - et al.
Hybrid experimental–numerical analysis of basic ductile fracture experiments for sheet metals
Int. J. Solids Struct.
(2010) - et al.
New biaxially loaded specimens for the analysis of damage and fracture in sheet metals
Int. J. Solids Struct.
(2017) - et al.
Fracture of an anisotropic rare-earth-containing magnesium alloy (ZEK100) at different stress states and strain rates: experiments and modeling
Int. J. Plast.
(2019) - et al.
A new approach for ductile fracture prediction on Al 2024-T351 alloy
Int. J. Plast.
(2012)
On the notch ductility of a magnesium-rare earth alloy
Mater. Sci. Eng. A
Void growth and coalescence in a magnesium alloy studied by synchrotron radiation laminography
Acta Mater
Calibration of ductile fracture criteria at negative stress triaxiality
Int. J. Mech. Sci.
Notched specimen under compression for ductile failure criteria
Mech. Mater.
Fracture prediction based on a two-surface plasticity law for the anisotropic magnesium alloys AZ31 and ZE10
Int. J. Plast.
Establishment and verification of two-component ductile fracture criteria with high estimation accuracy
Eng. Fract. Mech.
Modeling of 3D plastic anisotropy and asymmetry of extruded magnesium alloy and its applications in three-point bending
Int. J. Plast.
Modeling of shear ductile fracture considering a changeable cut-off value for stress triaxiality
Int. J. Plast.
An improved curvature method of detecting the onset of localized necking in Marciniak tests and its extension to Nakazima tests
Int. J. Mech. Sci.
Micromechanically-motivated phenomenological Hosford–Coulomb model for predicting ductile fracture initiation at low stress triaxialities
Int. J. Solids Struct.
Ductile fracture of aluminum 2024-T351 under proportional and non-proportional multi-axial loading: Bao–Wierzbicki results revisited
Int. J. Solids Struct.
Cited by (30)
Anisotropic elastoplasticity and fracture of SCFR-PEEK composites in complex biaxial loading: Experiments and modelling
2024, Composites Science and TechnologyGurson-Cohesive modeling (GCM) for 3D ductile fracture simulation
2024, International Journal of PlasticityA forming limit framework accounting for various failure mechanisms: Localization, ductile and cleavage fracture
2024, International Journal of PlasticityNew geometry-inspired numerical convex analysis method for yield functions under isotropic and anisotropic hardenings
2024, International Journal of Solids and StructuresAn uncoupled ductile fracture model considering void shape change and necking coalescence
2023, Engineering Fracture MechanicsInvestigation into the fracture behavior of ZK60 Mg alloy rolling sheet under different stress triaxiality
2023, Journal of Materials Research and Technology