Next Article in Journal
Densities of Liquid Tm2O3, Yb2O3, and Lu2O3 Measured by an Electrostatic Levitation Furnace Onboard the International Space Station
Previous Article in Journal
Estimation Method of Relative Slip in Fretting Fatigue Contact by Digital Image Correlation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Ductility Failure of Advanced High Strength Dual Phase Steel DP590 during Warm Forming Based on Extended GTN Model

School of Traffic and Vehicle Engineering, Shandong University of Technology, Zibo 255000, China
*
Author to whom correspondence should be addressed.
Metals 2022, 12(7), 1125; https://doi.org/10.3390/met12071125
Submission received: 9 June 2022 / Revised: 25 June 2022 / Accepted: 27 June 2022 / Published: 30 June 2022

Abstract

:
Under warm forming, the damage parameters of extended GTN model are easily affected by temperature, and the failure of materials under warm forming is less studied by using this model. In this paper, based on the extended GTN model, the damage parameters of DP590 at different temperatures are determined by experiment and numerical simulation, studying the trend of damage parameters changing with temperature. Through the analysis of shear specimens at different angles and temperatures, we studied the changes in shear damage and void damage under different conditions, discussing the influence of shear damage and void damage on competitive fracture failure under warm forming. We modify the damage by using a function based on stress triaxiality, and present a competitive failure equation considering temperature and stress triaxiality. It is found that the extended GTN model can be applied to the failure study of DP590 steel under warm forming.

1. Introduction

Continuum damage (CDM) and micromechanical damage are two representative methods to study the ductile damage of materials. For the CDM, Lemaitre [1] puts forward Lemaitre model by considering the influence of stress. Mashayekh [2] proposed an improved Lemaitre model and studied the forming process of materials through this model. This approach is derived on a thermodynamic framework and it takes no account of the mechanical behavior of micro defects in the material.
Gurson et al. [3,4,5] proposed that the failure of materials was mainly due to the nucleation, growth, and aggregation of microvoids in materials by studying the nucleation of spherical cavities. In order to more accurately explain the hardening of materials and the interaction between cavities, Tvergaard and Needleman imported [6] three extended parameters and the GTN model is obtained. GTN model has been widely used in the study of material failure and damage, but there are still some limitations. The ordinary GTN model is used to predict the triaxiality of low stress, that is, when the material with shear stress is destroyed, it cannot accurately describe the evolution law of void volume. Many experts and scholars have studied this. Nahshon and Hutchinson [7] introduce a new parameter to determine the development rate of shear damage without changing the mechanical behavior of materials under axisymmetric stress state, the model with this parameter can simulate the fracture form under shear stress under low stress triaxiality. Wei J et al. [8] extended the formula of shear deformation on damage evolution, so that the model can be applied to material failure under extensive stress triaxiality. Liang X [9] proposed a new shear damage mechanism based on McClintock considering the accumulation of void volumes in the shear band. Based on this model, Ying L [10] further calibrated the new parameters in the tensile process of 22MnB5, and established the corresponding numerical implementation method. He W [11] evaluated the evolution of void volume fraction and shear damage. The improved GTN model was used to predict the ductile fracture during piperotation, and good prediction results were obtained. Zao H A, et al. [12] proposed an improved shear propagation GTN model containing two independent damage mechanisms to predict the ductile fracture of materials under different stress states. Gatea S [13] and others used the extended GTN in the research of progressive forming and found that the extended GTN model improved the modeling accuracy of cracks compared with the original GTN model. Liu J [14] combined the GTN model with the maximum shear stress criterion and established a shear extended GTN model to predict the void volume fracture process of Q345. It is found that its prediction ability is better than the original GTN model. Through research, Sun Q [15] found that the identification parameters of shear modified GTN damage model can effectively characterize the mechanical properties, damage evolution and ductile failure of materials during small punch test. Li X [16] et al., established a coupling model by shear extended GTN damage model, and found that the shear size can promote the shear damage of materials. Wu P [17] proved that the main mode of failure changed from microvoid growth to shear damage under the condition of shear stress by superimposing hydrostatic pressure in numerical simulation. Zhou et al. [18] combined the damage mechanics concept of Lemaitre and Jean [19] with the GTN cavity growth model and proposed an extended GTN model. The model combines the concept of damage mechanics with the void growth model, and couples the void volume parameters and shear damage parameters into the yield function. It is proposed that the void volume damage and shear damage adopt a separate critical damage state. When the damage parameters are unified, the material will fail. In order to study the material failure and fracture more accurately, the finite element method has been widely used. Yao D [20] proposed a finite element aided test (fat) method to obtain the uniaxial full range constitutive relationship of A5083 steel and SS316L before failure. The fracture stress and strain were obtained by finite element analysis. Churyumova T A [21] using the experimental values of the relative crosssection reduction and fifinite element calculation of the stress triaxiality got the critical values of the modified Rice and Tracy fracture criteria. Aldakheel F [22] develop a new theoretical and computational framework for the phase fifield modeling of ductile fracture in conventional elastic–plastic solids under fifinite strain deformation, and verify its accuracy by numerical simulation.
To our knowledge, it is common to use the extended GTN model to study the failure of materials at room temperature. However, under the condition of warm forming, the damage parameters will be affected by temperature, it is necessary to recalibrate and understand the change trend of damage parameters with temperature. On the other hand, the evolution of shear damage and void damage under different temperatures and different stress states needs to be further studied. If we can establish influence of shear damage and void damage on competitive fracture failure, it will be very helpful to verify whether the extended GTN model can be applied to the failure study of DP590 steel under warm forming.
In this paper, the damage parameters of materials at different temperatures can be determined by employing experiments and numerical simulation. Through the analysis of shear specimens at different angles and temperatures, we will study the changes of shear damage and void damage under different conditions. The effects of shear damage and void damage on material failure under warm forming conditions will be also discussed. The damage is corrected by using a function based on stress triaxiality, and a failure equation considering temperature and stress triaxiality is given.

2. Experiments

2.1. Material

Advanced high strength dual phase steel is mainly composed of ferrite and martensite. In this study, we selected DP590 dual phase steel which produced by Shanghai Baosteel company. Its chemical composition shows in Table 1.
According to the Chinese GB/T 4338–2006 high temperature tensile test method for metallic materials. Uniaxial tensile tests were carried out on DP590 specimens with a thickness of 1 mm at 293K, 373K, 473K, 573K, 673K, and 773K. The tensile specimen is cut from the DP590 plate in the rolling direction (RD), 45 direction and transverse direction (TD). The test was carried out on WDW-20D material performance testing machine with the tensile speed 0.2 mm/min at different temperatures. The dimensions of uniaxial tensile specimens are shown in Figure 1.
Before the tensile test, heating the DP590 test specimen to the test temperature, and the temperature of the test specimen remains stable until the test specimen is stretched to fracture. Uniaxial tensile tests were carried out on the same three groups of samples at different temperatures. In order to ensure the accuracy of the results, the average value of each group of tests was taken. Test results are shown in Figure 2.
According to experimental results, when the temperature is higher than 673 K, the plasticity of DP590 will obviously improve. Because when the temperature is high, the warm forming of DP590 can reduce the obstruction of grain boundary to the internal dislocation movement of the material, improve the ability of dynamic recovery and recrystallization of the material, and promote the elimination of material defects and internal stress, which is reflected in the improvement of plasticity and the reduction of strength. When the temperature is 473 K–573 K, the strength of DP590 is high and the phenomenon of “blue brittleness” will appear. This is because C and N atoms anchor the dislocation rapidly, which makes it difficult for the dislocation to slip, forming the blue embrittlement phenomenon.
We use Equations (1) and (2) to convert the engineering stress–strain curve obtained from the test into the real stress–strain curve.
ε = ln ( 1 + ε nom )
σ = σ nom ( 1 + ε nom )
ε , ε nom stand for real strain and engineering strain respectively, σ , σ nom stand for real stress and engineering stress respectively.
Using the famous swift model to describe the mechanical properties of DP590 steel at different temperatures, taking 293 K as an example. As shown in Figure 3.
σ = K ( ε 0 + ε ¯ p ) n
where K is the hardening coefficient and n is the material hardening index.
The determined elastoplastic parameters are shown in Table 2.

2.2. Tensile Shear Test

In order to study whether the extended GTN model is applicable to the failure of DP590 dual phase steel at different temperatures and low stress triaxiality, we designed a series of shear specimens (as shown in Figure 4), and the shear tests at different temperatures were carried out by tensile testing machine.
Before the tensile shear test, heating the DP590 shear specimen to the test temperature, and the temperature of the specimen remains stable until the specimen is stretched to fracture. Figure 5 shows the partial shear specimen. Figure 6 shows the load–displacement curve of different angle shear specimens at different temperatures.

3. Constitutive Model

3.1. GTN Damage Model

Gurson model considers the factors of void volume damage. On the basis of Gurson model, tvegard and niedman modified it into GTN damage model by introducing void volume interaction parameters. GTN damage model is a commonly used constitutive model in meso damage analysis of materials. The yield function of GTN damage model is expressed as:
Φ = ( σ e q σ m ) 2 + 2 q 1 f * cos h ( 3 q 2 σ h 2 σ m ) [ 1 + q 3 f * 2 ] = 0
where σ e q is the Mises equivalent stress, σ h is the hydrostatic pressure, σ m is the equivalent stress of the material, q 1 , q 2 , q 3 is the strengthening parameter, f * is the void volume fraction function:
f * = { f f c + δ ( f f c ) f f c f > f c
δ = f u * f c f f f c
where f c is the critical void volume fraction when the voids begin to polymerize, f f is the void volume fraction when the material ruptures, δ is the void growth acceleration factor, f u * is the f * value when the stress in the yield equation is 0.

3.2. Extended GTN Model

GTN model is more accurate in predicting the necking fracture phenomenon of materials, but the effect of the model is not accurate in predicting the material failure dominated by shear damage. In the case of low stress triaxiality, the growth of void volumes is restrained, only a small number of void volumes can nucleate, the void volume fraction changes little, and the ordinary GTN model cannot accurately predict the failure of materials. In order to overcome these defects, Zhou et al. [18] introduced the concept of continuous damage into the GTN model, the shear damage variable is defined as DS. Shear damage only affects the deviatoric stress state, and the change of plastic volume is affected by porosity. The yield function of the extended GTN model can be expressed as:
Φ = ( σ e q σ m ) 2 + 2 q 1 f * cos h ( 3 q 2 σ h 2 σ m ) [ 1 + ( q 1 f * + D s ) 2 2 D s ] = 0
D = q 1 f * + D S
where D is the total damage, DS is shear damage.
Shear damage is regarded as a function of plastic strain and stress state. The definition formula is as follows:
D s = ( ε m p l ε f s ) n
where ε f s is the fracture strain under pure shear state; n is the weakening index of the material and the ε m p l plastic strain of the matrix. In the above formula, stress triaxiality and Lode parameter are introduced and extended to any stress state.
D s · = ψ ( θ , T ) n D s n 1 n ε f s ε m p l ·
where ψ ( θ , T ) is the weight function [23].
ψ ( θ , T ) = { g ( θ ) g ( θ ) ( 1 k ) + k T > 0 T 0
where k is the weight factor when the triaxiality of stress is negative, which can be calibrated by the data of tensile test.
g ( θ ) = 1 6 | θ | π
where θ is Lode angle. At this time, the evolution of plastic strain is expressed material equivalent plastic power and macro equivalent power:
ε m p l · = σ : d ε p ( 1 D / q 1 ) σ m
where ( 1 D / q 1 ) can ensure that here is no shear damage.

4. Parameters Calibration

4.1. The Yield Surface Coefficients

Before applying the extended GTN damage model to predict the damage of dual phase steel at different temperatures, 11 uncertain parameters in the model need to be calibrated. Among them, the extended GTN model is used to calibrate n and ε f parameter, and the GTN model is used to calibrate other parameters. By consulting a number of literature [24,25,26,27,28,29,30,31,32], most researchers use the strengthening coefficient q1 = 1.5, q2 = 1, q3 = 2.25, proposed by Tvergaard when using GTN model to study the failure of different materials. In this paper, the coefficients are selected according to the results of literature review.

4.2. Nucleation Rate Function Parameters and Void Nucleation Parameters

For the standard deviation of nucleation strain Sn, according to the above literature research, it can be found that the standard deviation of nucleation strain is generally taken as Sn = 0.1.
For the normal temperature forming of DP590, the average ε N strain void volume nucleation is generally taken as 0.3. According to Liu Wenquan’s research [33], when the forming temperature of highstrength steel is considered, the average strain of void volume nucleation should consider the influence of temperature. Therefore, ε N is the strain value when strain local deformation occurs in the tensile sample at different temperatures, as shown in Table 3.
In this paper, the parameter inverse method is used to identify the damage parameters at different temperatures. In the process of parameter calibration, central composite design and response surface analysis method are adopted [34,35,36], and the peak point and fracture point of the material tensile curve are used as the key positions of damage parameter calibration. The load R1 at the peak point, the displacement R2 at the peak point, the load R3 at the fracture point, and the displacement R4 at the fracture point are taken as the optimization objectives. When the error between the prediction curve and the test curve is small, the appropriate damage parameter value can be found taking the damage parameter calibration at 573 K as an example. Table 4 shows the damage parameter coefficient levels of the central composite design.
Table 5 lists 28 groups of damage parameter combinations and response values obtained from numerical simulation.
Analysis of ANVOA is an important method to evaluate the significance of response surface model. According to the analysis of ANOVA, we can find out whether the regression model is applicable.
Through ANVOA analysis (as shown in Table 6), it can be seen that the model is significant and effective. Therefore, the regression model of damage parameters at 573 K temperature is:
R 1 = 6.791 + 17.96 f c 374.03 f 0 68.49 f n 1.348 f f + 10986.04 f c f 0 + 876.54 f 0 f f 814.77 f c 2 + 1240.94 f n 2
R 2 = 7024.08 + 91.04 f c 13767.94 f 0 + 397.68 f n 60.98 f f 8062.16 f c f n + 13513.88 f 0 f f + 802566 f 0 2 3023.72 f n 2 + 76.531 f f 2
R 3 = 11.628 152.948 f c 487.46 f 0 + 55.61 f n 4.66 f f + 5365.07 f c f 0 1143.65 f c f n + 7493 f 0 f f + 1816.03 f c 2 + 14.13 f f 2
R 4 = 4472.325 + 48981.1 f c + 546198 f 0 24768.43 f n + 1690.19 f f 1258571 f c f 0 + 117353 f c f n 1703910 f 0 f n + 341183 f 0 f f 81471 f c 2 48665100 f 0 2 41569.49 f f 2
Finally, the damage parameter values under different temperatures are obtained by genetic algorithm (Table 7). The flow chart of genetic algorithm is shown in the Figure 7.
According to the above damage parameters, we write VUMAT subroutine and use ABAQUS software to simulate the damage parameters at different temperatures. We compare the load displacement curve obtained by numerical simulation with the test load displacement curve, and the results are shown in Figure 8.
By comparing the curves, we found that the maximum error of the curve is 4.72%. Therefore, it can be considered that the GTN damage parameters obtained by numerical simulation combined with the test and RSM experimental design are accurate and can be used to predict the failure of unidirectional tensile specimens at different temperatures. Figure 9 shows the variation law of void volumerelated damage parameters with temperature.
It can be found that the initial void volume fraction increases slightly with the increase of temperature. When the temperature is high, the material will recover and recrystallize, so as to reduce the internal defects, reduce the volume fraction f N of cavity nucleation, and increase the critical void volume fraction f c and fracture void volume fraction f f .

4.3. Shear Damage Parameters

The shear damage introduced into the extended model is an important state variable. The evolution of shear damage is affected by the material shear damage softening parameter n and effective failure strain. The shear parameters were calibrated by shear test and numerical simulation. Taking 293K as an example, the load–displacement curve under two parameter sets (Figure 10) and other model parameters take the above given values, and finally determine the value of n and ε f at different temperatures.
The load–displacement curve test and numerical simulation curve of 0°shear parts at different temperatures are compared as shown in Figure 10. It can be seen that the calibrated shear damage parameters are appropriate.

5. Verification and Discussion

In this section, the damage parameters of the above test and calibration are numerically simulated by using VUMAT subroutine in the previous section and using ABAQUS software to study the ductile failure of DP590 at different temperatures.

5.1. Finite Element Model

According to the shear tests and the damage parameters calibrated above, the numerical simulation of uniaxial tensile shear test is carried out for shear specimens from other angles. In order to better capture local deformation and damage, refined meshes and coarse meshes are used inside and outside the central region respectively. The size of the refined mesh is 0.2 mm × 0.2 mm, and the size of the coarsened mesh is 0.8 mm × 0.8 mm. Updated Lagrange formulations were used to handle the geometrical and material nonlinearities, and mass scaling factor was set to 20 so as to reduce the total simulation time. Fixed constraint was imposed at the bottom of the samples and the displacement loading was implemented on the top boundary. The finite element model of the test piece is shown in Figure 11.

5.2. Numerical Simulation Results

The numerical simulation results such as equivalent plastic strain, vvf, etc., are shown in Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16 below (take the shear specimen at 673 K as an example).
It can be seen from the above Figure 12, Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17 that in the 0° tensile test, the shear damage starts at both ends of the sample notch and propagates to the center, and the crack propagation direction is parallel to the loading direction. With the increasing of specimen, the influence of shear damage becomes smaller and smaller, the initial position of shear damage gradually shifts to the center, and the crack propagation direction also develops vertical to the loading direction. When any variable of shear damage or void volume damage reaches the critical value, the material will fail.
Through numerical simulation, changes in void volume damage and shear damage at different temperatures are obtained. It can be seen from the above Figure 18 and Figure 19 that for samples with the same temperature and different angles, the influence of void volume damage on material damage increases with the increase of shear angle. This is consistent with the gradual development of fracture form from shear fracture to tensile fracture. For the shear specimen with the same angle, when the temperature changes from 293 K to 573 K, the void volume damage value is small under the same equivalent plastic strain. When the temperature changes from 573 K to 773 K, the void volume damage value will gradually increase. This shows that when the formability of DP590 becomes better, the void volume damage has a great influence on the failure of DP590.
In order to more accurately measure the influence of shear damage and void damage on competitive fracture failure, introduce a weight coefficient to describe the evolution process of shear damage and void damage.
D = ( T T L ) 2 × q 1 f * + 1 ( ( T T L ) 2 ) × D S ( T < T L )
where T L is the stress triaxiality in pure tension, which can be calibrated by uniaxial tensile test specimen and T is the stress triaxiality. With the continuous increase of stress triaxiality, the influence of shear damage decreases gradually. When the stress triaxiality of the specimen exceeds the critical stress triaxiality, the total damage is caused by void damage. In this way, the shear damage and void damage in the case of low stress triaxiality can be corrected, which can better conform to the evolution of total damage.
In order to find the influence of shear damage and void damage on competitive fracture failure and triaxiality of temperature and stress, taking 673 K as an example to extract the shear damage and stress triaxiality before the critical damage of materials as shown in Table 8 and we used a multiple linear regression equation to describe the shear damage and stress triaxiality at different temperatures.
D S = 0.937 2.46 T 3.85 e 7 × t
It can be seen that temperature has little effect on shear damage, so we choose quadratic polynomial to more accurately describe the relationship between shear damage and stress triaxiality.
D s = 4.49 × T 2 4.17 × T + 1.02
Because f c has little change under low stress triaxiality, only look for the relationship between f c and temperature. By fitting the relationship between f c and temperature, we get the expression of critical void volumes at different temperatures.
f c = 91.36 91 × exp ( 0.5 × ( t 495 ) 2 12816 )
When the stress triaxiality is less than the critical stress triaxiality, shear damage is the dominant factor causing material damage. With the continuous increase of stress triaxiality, void damage gradually dominates. Therefore, establish a failure expression based on temperature and stress triaxiality with T L as the threshold value.
{ D s > 4.49 × T 2 4.17 × T + 1.02 f c > 91.36 91 × exp ( 0.5 × ( t 495 ) 2 12816 ) T T L T > T L
when the stress triaxiality is greater than T L , the specimen will undergo tensile fracture. When the stress triaxiality is less than T L , the specimen will shear. At this time, the critical value can be determined by D S and f c respectively.

5.3. Result Analysis

According to the damage parameters determined above and the introduced correction function, we numerically simulate the shear specimens with different temperatures and angles, and extract the load–displacement curve (Figure 20).
Comparing the load–displacement curve in the Figure 20, the maximum error is within 5%. It can be considered that the extended GTN damage parameters obtained by combining test and numerical simulation are accurate, the introduced correction function can better predict the damage and the extended GTN model can be applicable to the shear damage prediction of DP590 dual phase steel under warm forming conditions. In order to verify the accuracy of the competitive failure of the extended GTN model, the tensile bending test and numerical simulation of DP590 are carried out. The results are shown in Figure 21.
By analyzing the numerical simulation results of different fillet radii, it is found that the fillet radius of 1 mm is shear fracture, the fillet radius of 5 mm is mixed fracture, and the fillet radius of 15 mm is tensile fracture. The fracture form and fracture position are consistent with the test results. It is verified that the competitive failure can predict the fracture of DP590 under different stress states.

6. Conclusions

In this paper, the damage parameters of materials at different temperatures are determined by experiment and numerical simulation. Through the analysis of shear specimens at different angles and temperatures, the changes of shear damage and void damage under different conditions are studied, and the effects of shear damage and void damage on material failure under warm forming conditions are discussed. A function based on stress triaxiality is used to modify the total damage D, and a competitive failure equation considering temperature and stress triaxiality is given.
  • Uniaxial tensile tests are carried out on plane specimens of DP590 dual phase steel at different temperatures and the central composite design (CCD)—response surface method (RSM)—genetic algorithm (GA) method was used to back correct the elasticplastic mechanical parameters and extended GTN damage parameters at different temperatures. By comparing the load–displacement curves of test and numerical simulation, the calibrated damage parameters are found to be accurate, and was found that the damage parameter f0 of the extended GTN model changes little with temperature at different temperatures. At higher temperature, the critical void volume integral fc and fracture void volume fraction ff increase due to the good recovery and recrystallization ability of the material. At the same time, due to the reduction of internal defects, the nucleation volume fraction fn of the cavity reduced.
  • Combined with test results, we simulated and analyzed the shear parts with different angles at different temperatures. For specimens with the same temperature and different angles, with the increase of stress triaxiality, the influence of void volume damage on material damage is greater. Moreover, under the same stress triaxiality, when the formability of DP590 increases, the void volume damage will gradually dominate.
  • Using a function based on stress triaxiality to modify the total damage D, by extracting the load–displacement curves of tensile shear numerical simulation at different temperatures and angles, the maximum error is 4.75%, which shows that the correction of this function is accurate. A competitive failure equation was presented to judge the critical damage of materials. Through tensile bending test and numerical simulation, we find that the fracture form and fracture position are consistent, which shows that the competitive failure equation is effective. Experiment and simulation show that the extended GTN model can be applied to the study of ductile fracture of DP590 under warm forming, and can obtain accurate ductile fracture prediction.

Author Contributions

Conceptualization, X.L. and D.L.; methodology, X.L.; software, Z.L. and H.S.; validation, X.L., D.L. and H.C.; formal analysis, X.L.; investigation, X.L.; resources, J.X.; data curation, N.J.; writing—original draft preparation, X.L.; writing—review and editing, D.L.; visualization, X.L.; supervision, N.J.; project administration, D.L.; All authors have read and agreed to the published version of the manuscript.

Funding

The Project supported by Natural Science Fund of Shandong Province (ZR2020QA040) and China SinoTruck Group industryuniversityresearch cooperation program project (No. 20200810).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lemaitre, J. A Continuous Damage Mechanics Model for Ductile Fracture. J. Eng. Mater. Technol. 1985, 107, 83–89. [Google Scholar] [CrossRef]
  2. MashayeKhi, M.; Torabian, N.; Poursina, M. Continuum damage mechanics analysis of strip tearing in a tandem cold rolling process. Simul. Model. Pract. Theory 2011, 19, 612–625. [Google Scholar] [CrossRef]
  3. Mcclintock, F.A. Erratum: A Criterion for Ductile Fracture by the Growth of Holes. J. Appl. Mech. 1968, 35, 363–371. [Google Scholar] [CrossRef]
  4. 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] [Green Version]
  5. Gurson, A.L. Continuum theory of ductile rupture by void nucleation and growth: Part I dash yield criteria and flow rules for porous ductile media. J. Eng. Mater. Technol. Trans. ASME 1977, 99, 2–15. [Google Scholar] [CrossRef]
  6. Tvergaard, V.; NEe Dleman, A. Analysis of the cupcone fracture in a round tensile bar. Acta Metall. 1984, 32, 157–169. [Google Scholar] [CrossRef]
  7. Nahshon, K.; Hutchinson, J.W. Modification of the Gurson Model for shear failure. Eur. J. Mech. 2008, 27, 1–17. [Google Scholar] [CrossRef]
  8. Wei, J.; Li, Y.; Jie, S. Modified GTN model for a broad range of stress states and application to ductile fracture. Eur. J. Mech.—A/Solids 2016, 57, 132–148. [Google Scholar]
  9. Liang, X. Damage accumulation and fracture initiation in uncracked ductile solids subject to triaxial loading. Int. J. Solids Struct. 2007, 44, 5163–5181. [Google Scholar]
  10. Ying, L. On the numerical implementation of a shear modified GTN damage model and its application to small punch test. Int. J. Mater. Form. Off. J. Eur. Sci. Assoc. Mater. Form. Esaform 2018, 11, 527–539. [Google Scholar] [CrossRef]
  11. He, W.; Xu, W.; Shan, D. An extended GTN model for low stress triaxiality and application in spinning forming. J. Mater. Process. Technol. 2018, 263, 112–128. [Google Scholar]
  12. Zao, H.A.; Hao, Z.; Yh, A. An improved shear modified GTN model for ductile fracture of aluminium alloys under different stress states and its parameters identification—ScienceDirect. Int. J. Mech. Sci. 2018, 192, 106081. [Google Scholar]
  13. Gatea, S.; Xu, D.; Ou, H. Evaluation of formability and fracture of pure titanium in incremental sheet forming. Int. J. Adv. Manuf. Technol. 2018, 95, 625–641. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, J.; Shen, Y.; Zhao, X. Simulation of fracture of a tubular Xjoint using a shearmodified GursonTvergaardNeedleman model. Thin-Walled Struct. 2018, 132, 120–135. [Google Scholar]
  15. Sun, Q.; Lu, Y.; Chen, J. Identification of material parameters of a shear modified GTN damage model by small punch test. Int. J. Fract. 2020, 222, 25–35. [Google Scholar] [CrossRef]
  16. Li, X.; Chen, Z.; Dong, C. Size effect on the damage evolution of a modified GTN model under high/low stress triaxiality in mesoscaled plastic deformation. Mater. Today Commun. 2020, 26, 101782. [Google Scholar] [CrossRef]
  17. Wu, P. Numerical Study about the Influence of Superimposed Hydrostatic Pressure on Shear Damage Mechanism in Sheet Metals. Metals 2021, 11, 1193. [Google Scholar]
  18. Zhou, J.; Gao, X.; SobotKa, J.C. On the extension of the Gursontype porous plasticity models for prediction of ductile fracture under sheardominated conditions. Int. J. Solids Struct. 2014, 51, 3273–3291. [Google Scholar] [CrossRef] [Green Version]
  19. Lemaitre, J. A Course on Damage Mechanics; Springer: Berlin/Heidelberg, Germany, 1996. [Google Scholar]
  20. Yao, D.; Cai, L.; Bao, C. A new fracture criterion for ductile materials based on a finite element aided testing method. Mater. Sci. Eng. A 2016, 673, 633–647. [Google Scholar] [CrossRef]
  21. Churyumova, T.A. United Approach to Modelling of the Hot Deformation Behavior, Fracture, and Microstructure Evolution of Austenitic Stainless AISI 316Ti Steel. Appl. Sci. 2021, 11, 3204. [Google Scholar] [CrossRef]
  22. Aldakheel, F.; Wriggers, P.; Miehe, C. A modified Gursontype plasticity model at finite strains: Formulation, numerical analysis and phasefield coupling. Comput. Mech. 2018, 62, 815–833. [Google Scholar] [CrossRef]
  23. Sun, C.; Li, D.; Zhao, L. Meso damage parameter calibration method based on power exponent hardening model. J. Henan Univ. Sci. Technol. (Nat. Sci. Ed.) 2016, 37, 7–10. (In Chinese) [Google Scholar]
  24. Schmitt, W.; Sun, D.Z.; Blauel, J.G. Damage mechanics analysis (gurson model) and experimental verification of the behaviour of a crack in a weldcladded component. Nucl. Eng. Des. 1997, 174, 237–246. [Google Scholar] [CrossRef]
  25. Hambli, R. Comparison between lemaitre and gurson damage models in crack growth simulation during blanKing process. Int. J. Mech. Sci. 2001, 43, 2769–2790. [Google Scholar] [CrossRef]
  26. Achouri, M.; Germain, G.; Santo, P.D.; Saidane, D. Numerical integration of an advanced gurson model for shear loading: Application to the blanKing process. Comput. Mater. Sci. 2013, 72, 62–67. [Google Scholar] [CrossRef] [Green Version]
  27. Springmann, M.; Kuna, M. Identification of material parameters of the gursontvergaardneedleman model by combined experimental and numerical techniques. Comput. Mater. Sci. 2005, 33, 500. [Google Scholar] [CrossRef]
  28. Kiran, R.; Khandelwal, K. Gurson model parameters for ductile fracture simulation in ASTM A992 steels. Fatigue Fract. Eng. Mater. Struct. 2014, 37, 171–183. [Google Scholar] [CrossRef]
  29. Kami, A.; Dariani, B.M.; Vanini, A.S.; Comsa, D.S.; Banabic, D. Numerical determination of the forming limit curves of anisotropic sheet metals using GTN damage model. J. Mater. Process. Technol. 2014, 216, 472–483. [Google Scholar] [CrossRef]
  30. Lemiale, V.; Chambert, J.; Picart, P. Description of numerical techniques with the aim of predicting the sheet metal blanking process by FEM simulation. J. Mater. Process. Technol. 2009, 209, 2723–2734. [Google Scholar] [CrossRef]
  31. Achouri, M.; Germain, G.; Dal Santo, P.; Saidane, D. Experimental and numerical analysis of micromechanical damage in the punching process for High-Strength Low-Alloy steels. Mater. Des. 2014, 56, 657–670. [Google Scholar] [CrossRef] [Green Version]
  32. Zhao, P.J.; Chen, Z.H.; Dong, C.F. Experimental and numerical analysis of micromechanical damage for DP600 steel in fineblanking process. J. Mater. Process. Technol. 2016, 236, 16–25. [Google Scholar] [CrossRef]
  33. Liu, W. Research on Damage Evolution, Mechanical Behavior and Formability Prediction of Warm Stamping. Ph.D. Thesis, Dalian University of Technology, Dalian, China, 2018. (In Chinese). [Google Scholar]
  34. Ying, L.; Liu, W.; Wang, D.; Hu, P. Parameter Calibration of GTN Damage Model and Formability Analysis of 22MnB5 in Hot Forming Process. J. Mater. Eng. Perform. 2017, 26, 5155–5165. [Google Scholar] [CrossRef]
  35. Wang, L.Y.; Li, L. Parameter Identification of GTN Model Using Response Surface Methodology for HighStrength Steel BR1500HS. J. Mater. Eng. Perform. 2017, 26, 3831–3838. [Google Scholar] [CrossRef]
  36. More, S.; Kumar, A.; Narasimhan, K. Parameter identification of GTN damage model using response surface methodology for single point incremental sheet forming of IF steel. Adv. Mater. Process. Technol. 2021, 12, 1–16. [Google Scholar] [CrossRef]
Figure 1. Dimension diagram of the tensile specimen.
Figure 1. Dimension diagram of the tensile specimen.
Metals 12 01125 g001
Figure 2. (a) Specimen test diagram and (b) engineering stress–strain curve of DP590 steel at different temperatures.
Figure 2. (a) Specimen test diagram and (b) engineering stress–strain curve of DP590 steel at different temperatures.
Metals 12 01125 g002
Figure 3. (a) Engineering and true stress–strain plots and (b) true plastic stress–strain approximate curve for realizing hardening parameters.
Figure 3. (a) Engineering and true stress–strain plots and (b) true plastic stress–strain approximate curve for realizing hardening parameters.
Metals 12 01125 g003
Figure 4. 0° shear specimen and other angles shear specimen ( α = 0 ° , 15 ° , 30 ° , 45 ° , 60 ° ) .
Figure 4. 0° shear specimen and other angles shear specimen ( α = 0 ° , 15 ° , 30 ° , 45 ° , 60 ° ) .
Metals 12 01125 g004
Figure 5. Partial shear specimen test diagram.
Figure 5. Partial shear specimen test diagram.
Metals 12 01125 g005
Figure 6. Load–displacement curves of shear specimens at different temperatures (af).
Figure 6. Load–displacement curves of shear specimens at different temperatures (af).
Metals 12 01125 g006aMetals 12 01125 g006b
Figure 7. Flow chart of genetic algorithm.
Figure 7. Flow chart of genetic algorithm.
Metals 12 01125 g007
Figure 8. Comparison curve of uniaxial tensile load and displacement.
Figure 8. Comparison curve of uniaxial tensile load and displacement.
Metals 12 01125 g008
Figure 9. Temperature variation of damage parameters f 0 f c f N and f f .
Figure 9. Temperature variation of damage parameters f 0 f c f N and f f .
Metals 12 01125 g009
Figure 10. (a) Numerical simulation diagram of shear damage calibrations and (b) comparison of load–displacement curves of 0°shear specimen at different temperatures.
Figure 10. (a) Numerical simulation diagram of shear damage calibrations and (b) comparison of load–displacement curves of 0°shear specimen at different temperatures.
Metals 12 01125 g010
Figure 11. (ae) Finite element grid diagram of shear specimen.
Figure 11. (ae) Finite element grid diagram of shear specimen.
Metals 12 01125 g011
Figure 12. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 0° specimen, and the result diagram of the specimen.
Figure 12. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 0° specimen, and the result diagram of the specimen.
Metals 12 01125 g012
Figure 13. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 15° specimen, and the result diagram of the specimen.
Figure 13. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 15° specimen, and the result diagram of the specimen.
Metals 12 01125 g013
Figure 14. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 30°specimen, and the result diagram of specimen.
Figure 14. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 30°specimen, and the result diagram of specimen.
Metals 12 01125 g014
Figure 15. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 45° specimen, and the result diagram of the specimen.
Figure 15. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 45° specimen, and the result diagram of the specimen.
Metals 12 01125 g015
Figure 16. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 60°specimen and the result diagram of specimen.
Figure 16. Equivalent plastic strain, VVF, shear damage at the beginning of microcrack in 60°specimen and the result diagram of specimen.
Metals 12 01125 g016
Figure 17. Test diagram before fracture of 0°, 30°, and 60° specimens.
Figure 17. Test diagram before fracture of 0°, 30°, and 60° specimens.
Metals 12 01125 g017
Figure 18. VVF and equivalent plastic strain diagram at different temperatures (af).
Figure 18. VVF and equivalent plastic strain diagram at different temperatures (af).
Metals 12 01125 g018
Figure 19. Shear damage and equivalent plastic strain diagram at different temperatures (af).
Figure 19. Shear damage and equivalent plastic strain diagram at different temperatures (af).
Metals 12 01125 g019aMetals 12 01125 g019b
Figure 20. Comparison between numerical simulation and test load–displacement curves at different temperatures (af).
Figure 20. Comparison between numerical simulation and test load–displacement curves at different temperatures (af).
Metals 12 01125 g020
Figure 21. Results and numerical simulation diagram of tensile bending specimens with different fillets (the rounded corners from left to right are 15 mm, 5 mm, and 1 mm, respectively).
Figure 21. Results and numerical simulation diagram of tensile bending specimens with different fillets (the rounded corners from left to right are 15 mm, 5 mm, and 1 mm, respectively).
Metals 12 01125 g021
Table 1. The chemical composition of DP5900 and its mass fraction%.
Table 1. The chemical composition of DP5900 and its mass fraction%.
CSiMnPSAl
0.080.461.750.0140.0040.032
Table 2. Elastoplastic mechanical parameters.
Table 2. Elastoplastic mechanical parameters.
Temperature (K)E (Mpa)vK (Mpa)n
293215,0000.311450.198
373212,0000.39710.137
473198,0000.310640.129
573182,0000.311120.149
673176,0000.36540.06
773165,0000.33660.063
Table 3. Values of ε N at different temperatures.
Table 3. Values of ε N at different temperatures.
Temperature (K) ε N
2930.3
3730.289
4730.275
5730.271
6730.36
7730.45
Table 4. Range of values of damage parameters.
Table 4. Range of values of damage parameters.
Damage ParametersLow LevelHigh LevelIntermediate Level
f 0 0.0010.00280.0019
f c 0.0180.0390.00285
f N 0.0230.0410.032
f f 0.120.220.17
Table 5. Design scheme of damage parameters and response values.
Table 5. Design scheme of damage parameters and response values.
f c f 0 f N f f R 1 R 2 R 3 R 4
0.0180.0010.0410.125.857015.3810.534853.26
0.0180.00280.0410.125.856999.5110.264500.2
0.00750.00190.0320.175.407003.711.434561
0.0390.0010.0230.125.3867013.548.195957.3
0.0390.00280.0230.225.856994.298.15956.5
0.0390.0010.0410.225.4907008.048.7295692
0.0390.0010.0230.225.21987007.128.1895928.15
0.02850.00190.0320.175.6707001.698.815817
0.0180.0010.0230.125.75987011.89.715089.2
0.04950.00190.0320.175.2206998.548.115762.78
0.0180.0010.0230.225.75987011.89.84472.6
0.0180.00280.0230.225.58016996.779.185595.14
0.02850.00190.0320.175.6707001.698.815817
0.0180.00280.0230.125.58016996.779.265387.09
0.0390.00280.0230.125.4906992.418.15895.69
0.02850.00190.0320.075.6707001.698.9995374.74
0.02850.00010.0320.175.4007017.389.185303.39
0.02850.00190.0320.175.6707001.698.95652.21
0.02850.00190.0140.175.8227001.228.185869.5
0.0390.00280.0410.225.82669938.375885.26
0.02850.00190.0320.175.6707001.699.0255152.32
0.02850.00190.0320.275.7367001.259.2224836.23
0.02850.00190.050.176.3216998.239.724936.23
0.02850.00370.0320.175.8256989.238.925423.62
0.0180.00280.0410.225.88026996.7710.344422.26
0.0390.0010.0410.125.4907008.048.45986125
0.0390.00280.0410.125.77626992.418.645296.56
0.0180.0010.0410.225.8257011.810.1694255.11
Table 6. ANVOA analysis.
Table 6. ANVOA analysis.
Response R 1 R 2 R 3 R 4
Model p value<0.0001<0.0001<0.00010.0002
R 2 %93.2796.8598.2790.52
Table 7. Damage parameter values at different temperatures.
Table 7. Damage parameter values at different temperatures.
Temperature f N f 0 f c f f
2930.0160.00020.0450.136
3730.01950.00040.03950.128
4730.0210.000920.0370.126
5730.0320.00120.030.123
6730.0290.001320.0430.158
7730.0120.001380.0560.172
Table 8. Stress triaxiality and shear damage values at 673 K.
Table 8. Stress triaxiality and shear damage values at 673 K.
AngleStress TriaxialityShear Damage
00.00250.947
150.1330.586
300.2340.492
450.2820.350
600.3060.0786
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, X.; Li, D.; Song, H.; Lu, Z.; Cui, H.; Jiang, N.; Xu, J. Study on Ductility Failure of Advanced High Strength Dual Phase Steel DP590 during Warm Forming Based on Extended GTN Model. Metals 2022, 12, 1125. https://doi.org/10.3390/met12071125

AMA Style

Liu X, Li D, Song H, Lu Z, Cui H, Jiang N, Xu J. Study on Ductility Failure of Advanced High Strength Dual Phase Steel DP590 during Warm Forming Based on Extended GTN Model. Metals. 2022; 12(7):1125. https://doi.org/10.3390/met12071125

Chicago/Turabian Style

Liu, Xingfeng, Di Li, Hui Song, Zipeng Lu, Hongjian Cui, Ning Jiang, and Jiachuan Xu. 2022. "Study on Ductility Failure of Advanced High Strength Dual Phase Steel DP590 during Warm Forming Based on Extended GTN Model" Metals 12, no. 7: 1125. https://doi.org/10.3390/met12071125

APA Style

Liu, X., Li, D., Song, H., Lu, Z., Cui, H., Jiang, N., & Xu, J. (2022). Study on Ductility Failure of Advanced High Strength Dual Phase Steel DP590 during Warm Forming Based on Extended GTN Model. Metals, 12(7), 1125. https://doi.org/10.3390/met12071125

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