Next Article in Journal
Digital Design, Analysis and 3D Printing of Prosthesis Scaffolds for Mandibular Reconstruction
Next Article in Special Issue
A Macroscopic Strength Criterion for Isotropic Metals Based on the Concept of Fracture Plane
Previous Article in Journal
Metal Oxide Nanoparticle-Based Coating as a Catalyzer for A-TIG Welding: Critical Raw Material Perspective
Previous Article in Special Issue
Modeling of Ti6Al4V Alloy Orthogonal Cutting with Smooth Particle Hydrodynamics: A Parametric Analysis on Formulation and Particle Density
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Extended Iterative Identification Method for the GISSMO Model

1
State Key Laboratory of Mechanical Transmissions, Chongqing University, Chongqing 400044, China
2
School of Automobile Engineering, Chongqing University, Chongqing 400044, China
*
Author to whom correspondence should be addressed.
Metals 2019, 9(5), 568; https://doi.org/10.3390/met9050568
Submission received: 22 April 2019 / Revised: 8 May 2019 / Accepted: 13 May 2019 / Published: 15 May 2019
(This article belongs to the Special Issue Failure Mechanisms in Alloys)

Abstract

:
This study examines an extended method to obtain the parameters in the Generalized Incremental Stress State Dependent Damage (GISSMO) model. This method is based on an iterative Finite Element Method (FEM) method aiming at predicting the fracture behavior considering softening and failure. A large number of experimental tests have been conducted on four different alloys (7003 aluminum alloy, ADC12 aluminum alloy, ZK60 magnesium alloy and 20CrMnTiH Steel), here considering tests that span a wide range of stress triaxiality. The proposed method is compared with the two existing methods. Results show that the new extended Iterative FEM method gives the good estimate of the fracture behaviors for all four alloys considered.

1. Introduction

The research into lightweight and high-durability materials used in automobiles has greatly increased in recent years, due to an increased focus on preserving natural resources and reducing air pollution. One effective way to achieve weight reduction is through material development, but with this approach comes a demand for accurately predicting the flow and damage behaviors of the new material. [1,2] Flow behavior and damage evolution is, however, a challenging task.
A large amount of research has been focused on evaluating the flow and fracture properties through micromechanics-based consideration [3,4,5,6,7,8,9,10,11]. Among these are the porous plasticity models, the best-known being the classical model by Gurson [3] that includes the growth of voids in a metallic matrix and a yield condition that relates to porosity. The Gurson model and many of its subsequent extensions agree well with experiments; however, it fails to predict the material behavior under shear dominated loading conditions. To remedy this issue, Nahshon and Hutchinson [12] distinguished shear dominated states by the third invariant of stress. Dæhli et al. [13] extended the Gurson model by Lode-dependent void evolution, both successfully predicting fracture under low triaxiality. However, the Gurson model is computationally inefficient as the element size must scale with the dominant void sparking. Thus, the phenomenological fracture models outweigh the porous models in industrial applications, as these models often require fewer material constants to be determined through data fitting. The typical fracture models [14,15,16,17] are formalized on experimental observation, which directly related to the stress state. Identified through the stress, the triaxiality and the lode angle are proposed. Compared to the Gurson model, the Johnson-Cook (JC) model [17,18,19] predicts the fracture locus accounting for the strain, the strain rate and the temperature, and the model parameters in the JC model is easily to calibrate. Other fracture models in Ref. [14,15,20] predicts the fracture locus in 3D space with more accuracy. Moreover, it is worth noticing that these fracture models are uncoupled from flow behavior, which is described by the standard material model. The Generalized Incremental Stress State Dependent Damage (GISSMO) Model which is developed by Neukamm [21,22], fully describes the ductile damage accounting for material softening and fracture. The fracture behavior is described by the damage variable and the damage softening behavior is represented by an instability measure. As the GISSMO model meets the industrial requirements, it finds wide use in crashworthiness and forming process simulations [23,24,25,26,27].
The parameter calibration strongly influences the accuracy of a numerical prediction. The constitutive constants are classically obtained by graphical analysis [28] and linear regression methods [29] to fit experimental data. For example, the force-displacement curves obtained from experimental data is transformed into equivalent stress-strain curves, from which constitutive parameters are obtained [30]. However, these calculation methods are conducted under assumptions of relatively small strain variation, failing to predict fracture behavior in several industrial applications. Recently, an iterative FEM method [31,32,33,34] has been developed for parameter calibration that, based on trial and error, can determine a set of parameters, which leads to good agreements between experimental data and numerical results. Xiao [31] obtained the fracture parameters through the iterative FEM method to predict the fracture locus in the 3D stress state. However, fracture parameters by the iterative FEM method are difficult to predict the fracture behavior considering softening.
In this paper, an extended iterative FEM method is proposed to obtain the GISSMO parameters. The accuracy of the inverse iterative FEM method is then validated by considering four different types of alloy through an extensive experimental program that subject the materials in different stress states. FEM procedures are conducted for three sets of validated parameters obtained by three methods. The accuracy of the three methods is evaluated. These stress-strain curves and fracture modes observed throughout the tests are compared to the model predictions. Finally, a discussion on the iterative FEM method is conducted.

2. Materials and Methods

2.1. Experiment Method

Four kinds of alloy, listing 7003-T6 aluminum (Test # 1–7), ADC12 aluminum (Test # 8–11), ZK60 magnesium (Test # 12–18) and 20CrMnTiH steel (Test # 19–25), are designated as the testing material with the chemical composition shown in Table 1. These alloys are extensively used in aerospace and automobile industries. Figure 1 presents 25 different specimens for a wide range of stress states. The dimensions of tensile specimen are recommended by the GB/T228.1-2010 standard. All tests are performed in the ETM504C electronic universal testing machine (Suzhou Dymeter Automotive testing technology co., LTD, Suzhou, China) with a load range of (5–50) kN and a speed range of (0.001–500) mm/min. The deforming results are measured by the Digital Image Correlation (DIC) method, and the loading results are measured by force cell. For DIC method, signature is marked in undeformed specimen, and displacement is measured by searching the signature in deformed specimen. The crosshead speed of the testing machine is kept constant at 1.8 mm/min for smooth tensile tests and 0.2 mm/min for notched and shear tests. These tests are performed in the ambient air at the room temperature. Following experimental data is recorded: the stress-strain relation, fracture mode, and elongation.

2.2. Effective Stress and Strain Curves

Experimental obtained engineering curves cannot be directly used in parameter calibration. In small strain, engineering curves is transformed into true stress-strain curves by volume constancy law,
ε true = ln ( 1 + ε eng ) σ true = σ eng ( 1 + ε eng )
where σ true and ε true is true stress and true strain. σ eng and ε eng is engineering stress and strain. For some metals Equation (1) cannot compensate for the large necking effects, constitutive models are used for extrapolating the after-necking stress-strain curves. Taken smooth tension for example, the engineering curves are transformed into true stress-strain curves by volume constancy law. Then, the soften part in the true stress-strain curve is extrapolated by the power law of Ludwik [35]. The extrapolated curve is shown in Figure 2 as the yellow dash one. The effective stress-strain curve is obtained by subtracting off the elastic strain and inputting the true stress. The effective stress-strain curve is used in material model calibration. The effective strain at fracture point is input as initial value in iterative FEM method.

2.3. GISSMO Model

The GISSMO, developed by Neukamm et al. [21,22] and Basaran et al. [36], ensures flexibility for a wide range of metal as it can correctly predict damage regardless the details of the material model formulation. Besides, the GISSMO targets damage prediction that accounts for material instability, localization, and failure. Damage accumulation, which is illustrated by Johnson [17] and Xue [37], is based on an incremental formulation found from
Δ D = n ε f D ( 1 1 n ) ε ˙ p ,
This equation considers the nonlinear relation between plastic strain and internal damage [38,39], and the damage exponent n allows for a nonlinear accumulation of damage until failure. The equivalent plastic strain increment is denominated as ε ˙ p . ε f represents the equivalent plastic failure strain which allows for an arbitrary definition of triaxiality dependent failure strains ε f ( σ * ) by inputting a tabulated curve. Triaxiality is defined as p / σ m , in which p is the average main stress and σ m is the von Mises stress. D is the damage value, which accumulates in each element during deformation. When D = 1, the element is deleted.
Similar to damage accumulation, material instability is determined as
Δ F = n ε c F ( 1 1 n ) ε ˙ p ,
where ε c is the triaxiality dependent critical strain ε c ( σ * ) which acts as an activation for coupling damage and stress. F is the instability measure. If F = 0, the material is undeformed. F = 1 corresponds to the onset of localization.
As soon as F reaches unity, the element stress is reduced by
σ = σ ˜ [ 1 ( D D crit 1 D crit ) m ] ,
where σ is the modified stress and σ ˜ is the current stress. Dcrit is the damage only comes to some value when F reaches unity. m is a fading exponent intended to better depict the rate of material softening.

2.4. Extended Iterative Finite Element Method

The extended Iterative FEM (EIFEM) method is a trial-and-error method, aiming at obtaining stress and strain at large strain. The EIFEM method is inspired by the method in Ref. [31], which can obtain parameters in the JC model. This paper extended the model to calibrate parameters in the GISSMO model. Five parameters need to be calibrated, including stress triaxiality σ * , fracture strain ε f , critical strain ε c , fading exponent m, and damage exponent n. The proposed procedure of the extended Iterative FEM method (EIFEM) is,
  • Initial value of ε 0 f , ε 0 c , m0 and n0. ε 0 f is the fracture strain and ε 0 c is the necking strain at obtained effective stress-strain curve. m0 and n0 can be set to arbitrary value. in this paper, m0 = 1, n0 = 3.
  • Iteration till the shape of numerical force and displacement curve after necking coincide with experimental one. 3D FEM simulation by LS-DYNA is conducted to calculate elongation Δ L FEM . In simulation, the deforming process is predicted by the JC material model and the fracture model is the GISSMO model. If numerical shape differs experimental shape, n and m are modified by ni = ni – 1 1 and mi = mi − 1 + 0.5 till convergence.
  • Iteration till the experimental elongation coincides with the numerical elongation Δ L EXP = Δ L FEM . If Δ L EXP Δ L FEM , ε f would be modified by ε i   +   1 f × Δ L EXP = ε i f × Δ L i FEM . If ε i f = ε i c but Δ L FEM still unequal to Δ L EXP , ε c is modified by ε i + 1 c × Δ L n EXP = ε i c × Δ L n ( i ) FEM . Δ L n is the displacement at necking point. A satisfied pair of ε c and ε f is obtain when Δ L EXP = Δ L FEM .
  • Check if the standard deviation below 3%. The standard deviation (Std) between experimental and numerical curves is defined as
    Std = 0 l f | f Exp ( l ) f FEM ( l ) | d l Max { f Exp ( l ) } l f ,
    where l is displacement. The iteration is continued until the Std below 3%. Otherwise return to Step 2 and renew the number of m, n, ε c and ε f till convergence. In the last iterative FEM simulation, σ * is obtained at the fracture element right before deleting element.

2.5. Finite Element Method

The finite element (FE) analysis is conducted by the explicit mechanical solver of the commercial finite element code LS-DYNA. FE models for tensile tests are meshed through 8-node hexahedron elements by one-point integration, with single-point-constrain aligned to one side and the prescribed motion to the other side. The average mesh size is 0.2 for each specimen in gauge length. As the FEM mesh are similar between materials, the representative FEM mesh for ZK60 alloy are shown in Figure 3. Tensile specimens all move in constants speeds according to experimental speeds. Due to the strain rate insensitivity, the velocity is set at a relatively larger value compared to experiment to reduce the calculating time. The velocity is set at 10 mm/s in this paper. The plastic flow and damage are considered separately. The GISSMO model describes the fractured and damage behavior. In LS-DYNA, the GISSMO model is embedded in the ADD_EROSION which provides the failure and erosion option. The JC material model describes the flow behavior. The JC material model, which is one of the most widely used empirical models, provides an accurate prediction of the flow behavior considering the effects of stress state and strain rate. As all tests are performed at room temperature, the temperature effect is ignored. The strain rate part is also ignored as all experiments are in quasi-static state. The JC material model [40] is simplified as follows:
σ e = A + B ε p c ,
where A, B, and c are material constants and σ e is the flow stress. By achieving the good fit of smooth tensile experimental data, four parameter sets are calibrated and listed in Table 3. Please note that the elastic modulus is obtained from the elastic region ranging from start point to yielding point (the value stress value is A). This elastic modulus would make the numerical predicted results more accurate.
Fracture curves ε f   ( σ * ) , critical strain curves ε c   ( σ * ) , the fading exponent m and the damage exponent n are identified through the numerical simulation. Although different GISSMO parameters combinations may lead to a similar numerical result, the effect could be minimized by increasing the number of tests for GISSMO model calibration. Recall the method in Section 2.4, the parameters m and n can be calibrated. The fading exponent m and the damage exponent n are listed in Table 3. ε f   ( σ * ) and ε c   ( σ * ) are cubic spline interpolations of fracture strain and critical strain.

3. Results and Discussion

3.1. Experimental Results

The original experimental results are presented in Figure 4. The engineering stress is calculated by load, which is directly obtained by force cell. The engineering strain is obtained by DIC method. Recall the method in Section 2.2, the effective stress and strain curves are obtained. Then, the triaxiality σ * , critical strain ε c and fracture strain ε f are obtained based on the EIFEM method in Section 2.4. Then σ * , ε c and ε f is used in numerical simulation.

3.2. Numerical Results and Validation

Table 4 presents the experimental and numerical elongation. The relative errors are below 5% for all tests. Figure 5 shows the force and displacement curves for experiment and simulation. Figure 5a presents curves for tensile tests. Numerical results agree with the experimental results well. For 20CrMnTiH steel and 7003 aluminum alloy, the after-necking region in numerical force and displacement curves coincides with the experimental curves. The good agreements are also observed in notched and shear tests shown in Figure 5b–e.
Visual comparisons between experimental and numerical predicted fracture modes are shown in Figure 6, Figure 7, Figure 8 and Figure 9. Good agreements on the fracture modes are observed in those figures, indicating that the EIFEM method can reproduce the fracture behavior for different metals over a wide stress triaxiality with certain accuracy.

3.3. Comparison

To further study the application range of the proposed method, the numerical results predicted by the GISSMO model are compared with other two methods. The first method is widely used in industry for rough calculation. To get the fracture strain in the first method, the experimental obtained force and displacement curves are transposed into effective stress and strain based on method in Section 2.2. Then the fracture strain is obtained at the end point in the effective stress and strain curve. Usually this fracture strain is smaller than the EIFEM obtained fracture strain. To verify which fracture strain is more accurate, the numerical simulation for first method also carried out in LS-DYNA. For better clarification, the notation in first method adds A, so the fracture strain is written as ε A f . The second method is a commonly used iterative FEM method as in Ref. [31]. Different value of fracture strain can be obtained by this method, and the corresponding notation adds a B, so ε B f . The stepwise procedure for the second method is shown in Appendix A. The numerical simulation for validating the second method is also conducted.
Numerical simulations for the aforementioned two methods carry out in LS-DYNA. The material model is John-Cook model and the fracture model is the GISSMO model, both are consistent with the EIFEM method. The method A and the method B cannot obtain the critical strain ε c , the ε c in the two method are set to a relatively large constant. Then the GISSMO model cannot display the soften region. The value of ε c influence the value of ε f . In this paper, the ε c is set to 1 in method A and method B. m and n are set according to Table 3. The fracture strain vs. stress triaxiality for the three methods are shown in Figure 10.
Force and displacement relationships for method A and B is presented in Figure 11. The elongation is discussed firstly. Elongation calculated by method A is much smaller than experimental elongation. Elongations calculated by method B are closer to the experimental data. However, the relative errors for Test No. 9, 10, 17, 20 and 22 (numbered in Table 2) are beyond 5%, indicating the method B cannot accurately predict the failure behavior. Then the force and displacement curves are compared. For method A, the numerical force and displacement curves overlap the elastic region, failing to predict the harden and soften region. For method B, numerical force and displacement curves in Test No. 2, 3, 8–11,16–18, 23–25 match the experimental data. However, for other tests, the numerical curves cannot capture the soften region. The experiments well predicted by methods B include tests of ADC12 alloy and all shear tests. The experiments failing to be predicted by method B include the uniaxial tensile tests for 7003 alloy and 20CrMnTiH steel and notched tests for 7003 alloy, ADC12 alloy, and 20CrMnTiH steel. The well predicted experiments are brittle fracture (ADC12 alloy) and shear fracture, with relatively small strain and no soften effects. The bad estimated force and displacement curves are featured by obvious soften effects. So, the method B cannot predict soften effect. The reason lies in the constant value of critical strain ε c . ε c will influence the value of fracture strain ε f . Usually, ε f in condition ε c <   ε f is larger than condition ε c   ε f . So, the ε f obtain by the EIFEM method is larger than ε f by method B. Physically, the soften effect leads to the inaccurate estimate of fracture stain. Usually, ε EIFEM f ε B f ε A f . In conclusion, method A is not recommending unless for small strain conditions or very rough calculation. Method B can be safely used in deformation with little or no soften effects. The EIFEM method can accurately predict the large deformation and obvious soften conditions.

4. Conclusions

This paper validates the accuracy for an extended parameter identification method by experiments and simulations, followed by comparison with two existing methods. The EIFEM method obtains the parameters in the GISSMO fracture model. The numerical obtained elongation, force, and displacement curves and the fractured position agree well with experiment. The EIFEM method achieves the best estimate compared with other two methods. Method A is recommended to use in small strain conditions or rough calculation. Method B can be used in metal with little or no soften effects. The EIFEM method can be used in predicting the large deformation and obvious soften effects.

Author Contributions

Conceptualization, Y.X.; methodology, Y.X.; software, Y.X.; validation, Y.X.; formal analysis, Y.X.; data curation, Y.X.; writing—original draft preparation, Y.X.; writing—review and editing, Y.X. and Y.H.; visualization, Y.X.; supervision, Y.H.; project administration, Y.H.; funding acquisition, Y.H.

Funding

This research was funded by Fundamental Research Funds for the Central Universities, grant number 2018CDXYTW0031 and Graduate Scientific Research and Innovation Foundation of Chongqing, grant number CYB18060; Fundamental Research Funds for the Central Universities, grant number 0009005202001002.

Acknowledgments

Special thanks to Kim Lau Nielsen in the Technical University of Denmark for his comments and improvements to this paper. This work is also supported by State Key Laboratory of Mechanical Transmission 2017 Open Fund, grant number SKLMT-ZZKT-2017M15

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The iterated FEM method is used to obtain stress triaxiality σ * and fracture strain ε f . In stepwise Step L is the gauge length and Δ L is elongation.
Step 1: Obtaining the initial value of ε 0 f . ε 0 f =   Δ L / L .
Step 2: Numerical simulation for tensile tests. In this step, the FEM method is performed to simulate the tensile tests. Through the first attempt of tensile simulation, σ * can be obtained in the first fractured element at the fractured moment.
Step 3: Iteration of ε f . The elongation Δ L ( FEM ) is not identical to the Δ L ( EXP ) under ε 0 f . Therefore, the fracture strain is modified by ε i   +   1 f × Δ L ( EXP ) = ε i f × Δ L ( FEM ) to reduce the difference between Δ L ( FEM ) and Δ L ( EXP ) . This iteration is continued until the relative error between Δ L ( FEM ) and Δ L ( EXP ) below 3%. Through the FEM procedures, the precise values of stress triaxiality σ * and fracture strain ε f are identified.

References

  1. Hu, M.Q.; Zhang, J.J.; Sun, B.Z.; Gu, B.H. Finite element modeling of multiple transverse impact damage behaviors of 3-D braided composite beams at microstructure level. Int. J. Mech. Sci. 2018, 148, 730–744. [Google Scholar] [CrossRef]
  2. Spagnoli, A.; Terzano, M.; Brighenti, R.; Artoni, F.; Stahle, P. The fracture mechanics in cutting: A comparative study on hard and soft polymeric materials. Int. J. Mech. Sci. 2018, 148, 554–564. [Google Scholar] [CrossRef]
  3. Gurson, A.L. Continuum theory of ductile rupture by void nucleation and growth: Part I Yield criteria and flow rules for porous ductile media. J. Eng. Mater-T ASME 1977, 99, 2–15. [Google Scholar] [CrossRef]
  4. McClintock, F.A. A criterion for ductile fracture by the growth of holes. Int. J. Appl. Mech. 1968, 35, 363–371. [Google Scholar] [CrossRef]
  5. Nguyen, N.T.; Kim, D.Y.; Kim, H.Y. A continuous damage fracture model to predict formability of sheet metal. Fatigue Fract. Eng. Mater. Struct. 2013, 36, 202–216. [Google Scholar] [CrossRef]
  6. Nielsen, K.L.; Tvergaard, V. Ductile shear failure or plug failure of spot welds modelled by modified Gurson model. Eng. Fract. Mech. 2010, 77, 1031–1047. [Google Scholar] [CrossRef]
  7. Rice, J.R.; Tracey, D.M. On the ductile enlargement of voids in triaxial stress fields. J. Mech. Phys. Solids 1969, 17, 201–217. [Google Scholar] [CrossRef]
  8. Tvergaard, V.; Needleman, A. Analysis of the cup-cone fracture in a round tensile bar. Acta Metall. 1984, 32, 157–169. [Google Scholar] [CrossRef]
  9. Xue, L. Constitutive modeling of void shearing effect in ductile fracture of porous materials. Eng. Fract. Mech. 2008, 75, 3343–3366. [Google Scholar] [CrossRef]
  10. Cricrì, G. A consistent use of the Gurson-Tvergaard-Needleman damage model for the R-curve calculation. Fract. Struct. Integrity 2013, 7, 161–174. [Google Scholar] [CrossRef]
  11. Sepe, R.; Lamanna, G.; Caputo, F. A robust approach for the determination of Gurson model parameters. Fract. Struct. Integrity 2016, 10, 369–381. [Google Scholar] [CrossRef]
  12. Nahshon, K.; Hutchinson, J.W. Modification of the Gurson model for shear failure. Eur. J. Mech. A Solids 2008, 27, 1–17. [Google Scholar] [CrossRef]
  13. Dæhli, L.E.; Morin, D.; Børvik, T.; Hopperstad, O.S. A Lode-dependent Gurson model motivated by unit cell analyses. Eng. Fract. Mech. 2018, 190, 299–318. [Google Scholar] [CrossRef]
  14. Bao, Y.; Wierzbicki, T. On fracture locus in the equivalent strain and stress triaxiality space. Int. J. Mech. Sci. 2004, 46, 81–98. [Google Scholar] [CrossRef]
  15. Bao, Y.; Wierzbicki, T. A comparative study on various ductile crack formation criteria. J. Eng. Mater-T ASME 2004, 126, 314–324. [Google Scholar] [CrossRef]
  16. Bao, Y.; Wierzbicki, T. On the cut-off value of negative triaxiality for fracture. Eng. Fract. Mech. 2005, 72, 1049–1069. [Google Scholar] [CrossRef]
  17. Johnson, G.R.; Cook, W.H. Fracture characteristics of three metals subjected to various strains, strain rates, temperatures and pressures. Eng. Fract. Mech. 1985, 21, 31–48. [Google Scholar] [CrossRef]
  18. Senthil, K.; Iqbal, M.A. Effect of projectile diameter on ballistic resistance and failure mechanism of single and layered aluminum plates. Theor. Appl. Fract. Mech. 2013, 67–68, 53–64. [Google Scholar] [CrossRef]
  19. Xue, L.; Ling, X.; Yang, S.S. Mechanical behaviour and strain rate sensitivity analysis of TA2 by the small punch test. Theor. Appl. Fract. Mech. 2019, 99, 9–17. [Google Scholar] [CrossRef]
  20. Benzerga, A.A.; Surovik, D.A.; Keralavarma, S.M. On the path-dependence of the fracture locus in ductile materials–analysis. Int. J. Plasticity 2012, 37, 157–170. [Google Scholar] [CrossRef]
  21. Neukamm, F.; Feucht, M.; Haufe, A. Considering damage history in crashworthiness simulations. In Proceedings of the 7th European LS-DYNA Conference, Salzburg, Austria, 14–15 May 2009. [Google Scholar]
  22. Neukamm, F.; Feucht, M.; Haufe, A.; Roll, K. On closing the constitutive gap between forming and crash simulation. In Proceedings of the 10th Internation LS-DYNA Users Conference, Dearborn, MI, USA, 8–10 June 2008. [Google Scholar]
  23. Anderson, D.; Butcher, C.; Pathak, N.; Worswick, M.J. Failure parameter identification and validation for a dual-phase 780 steel sheet. Int. J. Solids. Struct. 2017, 124, 89–107. [Google Scholar] [CrossRef]
  24. Andrade, F.X.C.; Feucht, M.; Haufe, A.; Neukamm, F. An incremental stress state dependent damage model for ductile failure prediction. Int. J. Fract. 2016, 200, 127–150. [Google Scholar] [CrossRef]
  25. Effelsberg, J.; Haufe, A.; Feucht, M.; Neukamm, F.; Du Bois, P. On parameter identification for the GISSMO damage model. In Proceedings of the 12th International LS-DYNA® Users Conference, Dearborn, MI, USA, 3–5 June 2012. [Google Scholar]
  26. Haufe, A.; Feucht, M.; Neukamm, F. The Challenge to Predict Material Failure in Crashworthiness Applications: Simulation of Producibility to Serviceability. In Predictive Modeling of Dynamic Processes: A Tribute to Professor Klaus Thoma; Hiermaier, S., Ed.; Springer: Boston, MA, USA, 2009. [Google Scholar]
  27. Heibel, S.; Nester, W.; Clausmeyer, T.; Tekkaya, A.E. Influence of different yield loci on failure prediction with damage models. In 36th Iddrg Conference—Materials Modelling and Testing for Sheet Metal Forming; Volk, W., Ed.; IOP Publishing: Munich, Germany, 2017; Volume 896. [Google Scholar]
  28. Cingara, A.M.; McQueen, H.J. New formula for calculating flow curves from high temperature constitutive data for 300 austenitic steels. J. Mater. Process. Technol. 1992, 36, 31–42. [Google Scholar] [CrossRef]
  29. Khoddam, S.; Lam, Y.C.; Thomson, P.F. The effect of specimen geometry on the accuracy of the constitutive equation derived from the hot torsion test. Steel Res. 1995, 66, 45–49. [Google Scholar] [CrossRef]
  30. Hu, Y.; Xiao, Y.; Jin, X.; Zheng, H.; Zhou, Y.; Shao, J. Experiments and FEM simulations of fracture behaviors for ADC12 aluminum alloy under impact load. WITH Mater. Int. 2016, 22, 1015–1025. [Google Scholar] [CrossRef]
  31. Xiao, Y.; Tang, Q.; Hu, Y.; Peng, J.; Luo, W. Flow and fracture study for ZK60 alloy at dynamic strain rates and different loading states. Mater. Sci. Eng. A Struct 2018, 724, 208–219. [Google Scholar] [CrossRef]
  32. Gavrus, A.; Massoni, E.; Chenot, J.L. An inverse analysis using a finite element model for identification of rheological parameters. J. Mater. Process. Technol. 1996, 60, 447–454. [Google Scholar] [CrossRef]
  33. Gelin, J.C.; Ghouati, O. An inverse method for determining viscoplastic properties of aluminium alloys. J. Mater. Process. Technol. 1994, 45, 435–440. [Google Scholar] [CrossRef]
  34. Maniatty, A.; Zabaras, N. Method for solving inverse elastoviscoplastic problems. J. Eng. Mech. 1989, 115, 2216–2231. [Google Scholar] [CrossRef]
  35. Ludwik, P. Elemente der technologischen Mechanik; Springer: Berlin, Germany, 2013. [Google Scholar]
  36. Basaran, M.; Wölkerling, S.D.; Feucht, M.; Neukamm, F.; Weichert, D.; AG, D. An extension of the GISSMO damage model based on lode angle dependence. In Proceedings of the LS-DYNA Anwenderforum, Bamberg, Germany, 30 September–1 October 2010. [Google Scholar]
  37. Xue, L. Ductile fracture modeling: theory, experimental investigation and numerical verification. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, UK, 2007. [Google Scholar]
  38. Tasan, C.C. Micro-mechanical characterization of ductile damage in sheet metal. Ph.D. Thesis, Eindhoven University of Technology, Eindhoven, The Netherlands, 2010. [Google Scholar]
  39. Weck, A.G. The role of coalescence on ductile fracture. Ph.D. Thesis, McMaster University, Hamilton, ON, Canada, 2007. [Google Scholar]
  40. Johnson, G.R. A Constitutive Model and Data for Metals Subjected to Large Strains, High Strain Rates and High Temperatures. In Proceedings of the International Symposium on Ballistics, Hague, The Netherlands, 19–21 April 1983; pp. 541–548. [Google Scholar]
Figure 1. Schematic sketch of the specimen. specified value of R, H and θ can be found in Table 2. All dimensions are in mm.
Figure 1. Schematic sketch of the specimen. specified value of R, H and θ can be found in Table 2. All dimensions are in mm.
Metals 09 00568 g001
Figure 2. Stress-strain curves for experimental data processing. The red solid curve is the engineering stress-strain curves. The blue solid curve is transformed through the red one by plastic incompressibility condition. To compensate for the non-uniform deformation beyond necking, the blue solid curve is extrapolated by the power law of Ludwik. Yellow dash curve represents modified curve which could be used in calibrating constitutive model.
Figure 2. Stress-strain curves for experimental data processing. The red solid curve is the engineering stress-strain curves. The blue solid curve is transformed through the red one by plastic incompressibility condition. To compensate for the non-uniform deformation beyond necking, the blue solid curve is extrapolated by the power law of Ludwik. Yellow dash curve represents modified curve which could be used in calibrating constitutive model.
Metals 09 00568 g002
Figure 3. The representative FEM model for ZK60 alloy, (a) smooth tensile test, (b) pure shear test, (c) notched tensile test and (d) tensile and notched tests.
Figure 3. The representative FEM model for ZK60 alloy, (a) smooth tensile test, (b) pure shear test, (c) notched tensile test and (d) tensile and notched tests.
Metals 09 00568 g003
Figure 4. Engineering stress and strain curves for (a) 7003 aluminum alloy, (b) ADC12 aluminum alloy, (c) ZK60 magnesium alloy and (d) 20CrMnTiH Steel.
Figure 4. Engineering stress and strain curves for (a) 7003 aluminum alloy, (b) ADC12 aluminum alloy, (c) ZK60 magnesium alloy and (d) 20CrMnTiH Steel.
Metals 09 00568 g004
Figure 5. Numerical and experimental force and displacement curves. (a) Tensile tests, (b) notched and shear tests for 7003 aluminum alloy, (c) notched and shear tests for ADC12 aluminum alloy, (d) notched and shear tests for ZK60 magnesium alloy and (e) notched and shear tests for 20CrMnTiH steel. Scatters represent experimental data and solid lines represent numerical results hereinafter.
Figure 5. Numerical and experimental force and displacement curves. (a) Tensile tests, (b) notched and shear tests for 7003 aluminum alloy, (c) notched and shear tests for ADC12 aluminum alloy, (d) notched and shear tests for ZK60 magnesium alloy and (e) notched and shear tests for 20CrMnTiH steel. Scatters represent experimental data and solid lines represent numerical results hereinafter.
Metals 09 00568 g005
Figure 6. Fracture mode for 7003Al-T6 aluminum alloy for (a) smooth tensile test, (b) pure shear test, (c) tensile and shear test, (d) notched test for notched radii 5 mm, (e) notched test for notched radii 10 mm, (f) notched test for notched radii 15 mm and (g) notched test for notched radii 20 mm. Hereinafter, the numerical fracture mode is in plastic strain contour.
Figure 6. Fracture mode for 7003Al-T6 aluminum alloy for (a) smooth tensile test, (b) pure shear test, (c) tensile and shear test, (d) notched test for notched radii 5 mm, (e) notched test for notched radii 10 mm, (f) notched test for notched radii 15 mm and (g) notched test for notched radii 20 mm. Hereinafter, the numerical fracture mode is in plastic strain contour.
Metals 09 00568 g006
Figure 7. Fracture mode for ADC12 aluminum alloy for (a) smooth tensile test, (b) notched test for notched radii 4 mm, (c) notched test for notched radii 8 mm and (d) pure shear test.
Figure 7. Fracture mode for ADC12 aluminum alloy for (a) smooth tensile test, (b) notched test for notched radii 4 mm, (c) notched test for notched radii 8 mm and (d) pure shear test.
Metals 09 00568 g007
Figure 8. Fracture mode for ZK60 magnesium alloy for (a) smooth tensile test, (b) notched test for notched radii 1mm, (c) notched test for notched radii 5mm, (d) notched test for notched radii 10mm, (e) pure shear test, (f) shear and tensile test at shear angle 30° and (g) shear and tensile test at shear angle 60°.
Figure 8. Fracture mode for ZK60 magnesium alloy for (a) smooth tensile test, (b) notched test for notched radii 1mm, (c) notched test for notched radii 5mm, (d) notched test for notched radii 10mm, (e) pure shear test, (f) shear and tensile test at shear angle 30° and (g) shear and tensile test at shear angle 60°.
Metals 09 00568 g008
Figure 9. Fracture mode for 20CrMnTiH steel for (a) smooth tensile test, (b) notched test for notched radii 2 mm, (c) notched test for notched radii 0.8 mm, (d) notched test for notched radii 0.4 mm, (e) pure shear test, (f) shear and tensile test at shear angle 30° and (g) shear and tensile test at shear angle 60°.
Figure 9. Fracture mode for 20CrMnTiH steel for (a) smooth tensile test, (b) notched test for notched radii 2 mm, (c) notched test for notched radii 0.8 mm, (d) notched test for notched radii 0.4 mm, (e) pure shear test, (f) shear and tensile test at shear angle 30° and (g) shear and tensile test at shear angle 60°.
Metals 09 00568 g009
Figure 10. Strain and triaxiality curves for (a) 7003 aluminum alloy, (b) ADC12 aluminum alloy, (c) ZK60 magnesium alloy and (d) 20CrMnTiH steel.
Figure 10. Strain and triaxiality curves for (a) 7003 aluminum alloy, (b) ADC12 aluminum alloy, (c) ZK60 magnesium alloy and (d) 20CrMnTiH steel.
Metals 09 00568 g010
Figure 11. Force and displacement curves for (a) uniaxial tensile tests based on method A, (b) uniaxial tensile tests based on method B, (c) 7003Al alloy based on method A, (d) 7003Al alloy based on method B, (e) ZK60 Mg alloy based on method A, (f) ZK60 Mg alloy based on method B, (g) 20CrMnTiH based on method A and (h) 20CrMnTiH based on method B.
Figure 11. Force and displacement curves for (a) uniaxial tensile tests based on method A, (b) uniaxial tensile tests based on method B, (c) 7003Al alloy based on method A, (d) 7003Al alloy based on method B, (e) ZK60 Mg alloy based on method A, (f) ZK60 Mg alloy based on method B, (g) 20CrMnTiH based on method A and (h) 20CrMnTiH based on method B.
Metals 09 00568 g011
Table 1. Chemical composition of the materials used in the study (wt.%). The single value means the maximum value of element’s concentration.
Table 1. Chemical composition of the materials used in the study (wt.%). The single value means the maximum value of element’s concentration.
Materials AlCCrCuFeMgMnNiPPbSSnSiTiZrZn
7003AlBal.--0.10.350.5–10.3-----0.30.20.05–0.255.5–6.5
ADC12Bal.--1.5–3.51.20.30.50.5-0.1-0.19.6–12.0--1
ZK60-----Bal.--------0.45–0.94.8–6.2
CrMnTiH-0.17–0.231–1.350.3Bal.-0.8–1.150.30.04-0.04-0.17–0.370.01–0.1--
Table 2. Summary of experiments.
Table 2. Summary of experiments.
Mat.No.DescriptionGauge Length/mm
7003-T61Smooth tensile50
2Pure Shear19
3Tensile Shear, 45°, θ = 45°17.5
4Notched, R5, R = 5 mm, H = 10 mm20
5Notched, R10, R = 10 mm, H = 10 mm20
6Notched, R15, R = 15 mm, H = 8 mm40
7Notched, R20, R = 20 mm, H = 8 mm50
ADC128Smooth Tensile35
9Notched, R4, R = 4mm, H = 5mm14
10Notched, R8, R = 8mm, H = 5mm18
11Pure shear, flat40.6
ZK6012Smooth Tensile25
13Notched, R112
14Notched, R520
15Notched, R1030
16Pure Shear38
17Tensile Shear, 30°, θ = 30°21
18Tensile Shear, 60°, θ = 60°21
20CrMnTiH19Round smooth tensile20
20Round notch, R0.4, R = 0.4 mm, H = 9.2 mm10.4
21Round notch, R0.8, R = 0.8 mm, H = 8.4 mm10.8
22Round notch, R2, R = 2 mm, H = 6 mm12
23Pure shear38
24Tensile shear, 30°, θ = 30°21
25Tensile shear, 60°, θ = 60°21
Table 3. Parameters for numerical simulation.
Table 3. Parameters for numerical simulation.
Mat.E/GPaABcnm
7003-T6663482520.4432.5
ADC123311519380.67210
ZK60312213160.4321.5
20CrMnTiH1259447540.2833.5
Table 4. Summary of experimental and numerical results.
Table 4. Summary of experimental and numerical results.
Mat.No.Description σ * ε c ε f Elongation/mm
Exp.Num.R.E. (%) 1
7003-T61Smooth tensile0.330.30.798.168.301.80
2Pure Shear−0.010.30.831.941.844.87
3Tensile Shear, 45°−0.080.050.732.552.492.33
4Notched, R50.550.020.470.960.970.72
5Notched, R100.460.10.441.191.171.32
6Notched, R150.40.010.501.321.320.24
7Notched, R200.380.010.501.421.452.12
ADC128Smooth Tensile0.340.0050.0190.700.691.50
9Notched, R40.360.0050.0330.120.133.52
10Notched, R80.350.0050.0350.230.232.50
11Pure shear, flat−0.030.0050.0250.100.101.52
ZK6012Smooth Tensile0.330.080.6974.114.070.84
13Notched, R10.410.100.2620.360.351.51
14Notched, R50.460.050.1070.650.650.28
15Notched, R100.480.030.2891.131.093.39
16Pure Shear0.070.260.3430.270.270.77
17Tensile Shear, 30°0.220.050.4200.170.171.50
18Tensile Shear, 60°0.300.090.5870.310.311.10
20CrMnTiH19Round smooth tensile0.330.030.171.951.940.51
20Round notch, R0.40.500.030.080.130.134.83
21Round notch, R0.80.640.050.150.080.081.15
22Round notch, R20.690.030.090.130.124.78
23Pure shear−0.070.10.450.340.340.41
24Tensile shear, 30°0.20.10.420.210.222.57
25Tensile shear, 60°0.270.10.340.340.332.08
1 R.E. denotes the relative error.

Share and Cite

MDPI and ACS Style

Xiao, Y.; Hu, Y. An Extended Iterative Identification Method for the GISSMO Model. Metals 2019, 9, 568. https://doi.org/10.3390/met9050568

AMA Style

Xiao Y, Hu Y. An Extended Iterative Identification Method for the GISSMO Model. Metals. 2019; 9(5):568. https://doi.org/10.3390/met9050568

Chicago/Turabian Style

Xiao, Yue, and Yumei Hu. 2019. "An Extended Iterative Identification Method for the GISSMO Model" Metals 9, no. 5: 568. https://doi.org/10.3390/met9050568

APA Style

Xiao, Y., & Hu, Y. (2019). An Extended Iterative Identification Method for the GISSMO Model. Metals, 9(5), 568. https://doi.org/10.3390/met9050568

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop