Failure Modeling for QP980 Steel by a Shear Ductile Fracture Criterion

: The wide application of advanced high strength steels with high speciﬁc strength in the automotive industry can signiﬁcantly reduce energy consumption and contribute to carbon neutrality. Accurate prediction of the ductile fracture behavior of advanced high strength steels under complex stress states is of great signiﬁcance for its application in automobile industry. In this study, the ductile fracture behavior of QP980 under complex stress states, covering shear, uniaxial tension, and plane strain tension, is investigated by conducting the hybrid experiment and simulation. The pressure-coupled Drucker yield function is chosen to characterize the effect of stress states on yielding for QP980, considering its high accuracy compared with the von Mises yield function. Failure limit of the stress states is modelled by ﬁve uncoupled ductile fracture criteria (Brozzo, Oh, Rice-Tracey, Ko-Huh, and DF2012). To improve the numerical prediction accuracy, the parameters of the constitutive model are optimized by using the inverse engineering approach. The numerical predicted results are compared with the experimental load-stroke curves with the onset of fracture. The comparison indicates that the prediction error of the DF2012 criterion is signiﬁcantly lower than those of the other four criteria. In addition, the prediction accuracy is greatly improved with the parameters of the constitutive model optimized by the inverse engineering.


Introduction
Energy conservation is crucial to the automotive industry, which can effectively reduce the costs of automobiles and greenhouse gas emissions to promote carbon-neutral. Research shows that approximately 75% of fuel consumption is positively correlated with the weight of automobile. A 10% reduction in weigh could contribute to 6-8% reduction in energy consumption and 13% reduction in carbon dioxide emissions [1]. Therefore, the lightweight automobile can greatly reduce energy consumption. There are currently three main ways to achieve lightweight in the automotive industry, including topology optimization design, new lightweight materials, and new manufacturing technology. The application of new lightweight materials is beneficial to the weight loss of automobiles. Among them, QP980 of advanced high strength steel (AHSS) has high specific strength, not only reducing the weight but also increasing the safety of automobiles. The microstructure of QP980 at room temperature is a mixture of the martensite, ferrite, and retained austenite. The hard martensite improves the strength of QP980, while the soft ferrite enhances its ductility [2]. Consequently, QP980 is widely utilized in the automobile structures, such as B-pillar reinforcement plate. Nevertheless, ductile fracture is the main failure mode during the forming processes of AHSS. With the development of computer technology, the numerical simulation method can effectively predict the plastic deformation behavior of sheet metals in forming processes, thus avoiding the ductile fracture under wide loading conditions. Metals usually fail in virtue of the nucleation, growth, and coalescence of microscopic voids [3]. Many ductile fracture criteria have been proposed to predict the plastic deformation of sheet metals. Those ductile fracture criteria fall into two categories, namely coupled ductile fracture criteria considering the damage accumulation in constitutive model and uncoupled ones [4,5]. Bai and Wierzbicki [6] modified the Mohr-Coulomb criterion (MMC) and then successfully characterized the ductile fracture behavior of Al2024-T351 in a wide range of stress triaxiality. Zhang et al. [7] investigated the strain hardening behavior of AA5182-O aluminum alloy of under the stress states varying from shear to equibiaxial tension. Luo et al. [8] predicted the strain path changing effect on forming limits of AA 6111-T4 based on a shear ductile fracture criterion. Ghadikolaee et al. [9] studied the U-bending of AA6061-T6 aluminum alloy by using Ayada, Rice-Tracey, and normalized Cockroft-Latham fracture criteria. Chow and Jie [10] accurately predicted the forming limit of Al6022 aluminum alloy by using the Hill's quadratic anisotropic yield criterion based on the continuum damage mechanics. Luo and Wierzbicki [11] presented the failure of Dual Phase steel during stretch-bending operations by using the modified MMC. Xu et al. [12] proposed a new ductile criterion based on two typical fracture mechanisms, tension fracture and shear fracture, to predict the ductile fractures with stress triaxiality less than −1/3 for Al 6061-T6 and Al 2024-T351. Lou et al. [13,14] proposed ductile fracture criteria based on micromechanisms of ductile fracture: strain-controlled void nucleation, triaxiality-governed void growth, and shear coalescence of voids. Lou and Yoon [15] extended a stress-invariant-based function to model fracture limits of sheet metals. Mu et al. [16] developed a mathematical model of ductile fracture behavior by considering two major void deformation modes and calibrated the model for DP780 using a hybrid experimental-numerical method. On the one hand, all the above research are to investigate the ductile fracture behavior of sheet metals under different stress triaxiality, which can accurately characterize the deformation behavior of sheet metal under complex stress state. On the other hand, a single ductile fracture criterion is mainly used to simulate the ductile fracture behavior of sheet metals. Therefore, various ductile fracture criteria were used to characterize the ductile fracture behavior of QP980 under complex stress states.
In this paper, the plastic deformation behavior of QP980 sheet metal under various stress states of shear (in-plane shear specimen), uniaxial tension (specimen with a central hole), and plane strain tension (notched specimen) was investigated by conducting experiments and simulations. The material strength was subsequently predicted by using the von-Mises and pressure-coupled Drucker yield criteria. On this basis, four traditional uncoupled ductile fracture criteria and DF2012 criterion were used to simulate the ductile fracture of QP980. In addition, the parameters of the constitutive models and ductile fracture criteria were optimized by using an inverse engineering approach to improve the prediction accuracy.

Material
QP980 with the thickness of 1.0 mm was used to evaluate the ductile fracture under complex stress state during forming processes, and its chemical compositions are shown in Table 1. Four types of specimens shown in Figure 1 were fabricated to investigate the deformation behavior under uniaxial tension, plane strain tension, and shear stress states. The dog-bone specimen is used to characterize plastic behavior under uniaxial tension. The specimen with a central hole characterizes fracture at uniaxial tension, the notched specimen for fracture at plane strain tension, and the in-plane shear specimen for fracture at shear. It should be pointed out that prior fracture would occur at circular edge for the notched specimen. Accordingly, a parallel groove with 3 mm was embedded for the notched specimens to avoid stress concentration and ensure the fracture initiating from the notched specimen central during deformation. Specimens along the rolling direction (RD), diagonal direction (DD) and transverse direction (TD) respectively were prepared for each of the four types to study the effect of anisotropy on the deformation behavior. the notched specimen central during deformation. Specimens along the rolling direction (RD), diagonal direction (DD) and transverse direction (TD) respectively were prepared for each of the four types to study the effect of anisotropy on the deformation behavior.

Tensile Tests
Each specimen was stretched to ductile fracture by using a universal tensile testing machine with the maximum load capacity of 50 kN. In order to achieve a strain rate of 0.001/s, particular tensile velocities were defined for different specimens, as shown in Table 2 [17]. The initial gauge lengths are illustrated in Figure 1 for different specimens. The deformation processes of the specimens painted with stochastic patterns (Figure 1e) were recorded by using a XTOP digital image correlation (DIC) system at a framing rate of five per second.

Experimental Results
Figure 2 depicts experimental load-stroke response for dog-bone specimens along RD, DD, and TD. The curves with good repeatability demonstrate that the anisotropic strength presents weak in the uniaxial tension. The true longitudinal strain and true width strain of the dog-bone specimens are calculated by DIC and GOM software with the initial gauge length denoted in Figure 1a. The relationships of the longitudinal and width strains of the dog-bone specimens along different loading directions are demonstrated in Figure  3. It shows that the anisotropy in plastic deformation is insignificant for QP980. Accordingly, the plastic behavior is approximately isotropic for QP980.

Tensile Tests
Each specimen was stretched to ductile fracture by using a universal tensile testing machine with the maximum load capacity of 50 kN. In order to achieve a strain rate of 0.001/s, particular tensile velocities were defined for different specimens, as shown in Table 2 [17]. The initial gauge lengths are illustrated in Figure 1 for different specimens. The deformation processes of the specimens painted with stochastic patterns (Figure 1) were recorded by using a XTOP digital image correlation (DIC) system at a framing rate of five per second.

Specimens with a Central Hole
Notched Specimens

Experimental Results
Figure 2 depicts experimental load-stroke response for dog-bone specimens along RD, DD, and TD. The curves with good repeatability demonstrate that the anisotropic strength presents weak in the uniaxial tension. The true longitudinal strain and true width strain of the dog-bone specimens are calculated by DIC and GOM software with the initial gauge length denoted in Figure 1a. The relationships of the longitudinal and width strains of the dog-bone specimens along different loading directions are demonstrated in Figure 3. It shows that the anisotropy in plastic deformation is insignificant for QP980. Accordingly, the plastic behavior is approximately isotropic for QP980.   By analytical calculation of the load-stroke curve along RD in Figu hardening curve before necking, where the elastic strain is removed from is depicted in Figure 4. The maximum plastic strain before necking is appr The strain hardening curve is fitted by the Swift and Voce hardening law (1)-(3) respectively, illustrating the hardening behavior with good con necking. Nevertheless, the predicted results by using those two hardening l different after the necking strain. The Swift model predicts a higher load the sheet metal than that obtained from the Voce model. Hence, the combi model is used to modulate the predicted load capacity characterizing the st after necking. The parameters of the combined Swift-Voce model are show where K is the strengthen coefficient, n is the hardening exponent, and A, B   By analytical calculation of the load-stroke curve along RD in Figur hardening curve before necking, where the elastic strain is removed from t is depicted in Figure 4. The maximum plastic strain before necking is appro The strain hardening curve is fitted by the Swift and Voce hardening law (1)-(3) respectively, illustrating the hardening behavior with good cons necking. Nevertheless, the predicted results by using those two hardening la different after the necking strain. The Swift model predicts a higher loadi the sheet metal than that obtained from the Voce model. Hence, the combin model is used to modulate the predicted load capacity characterizing the str after necking. The parameters of the combined Swift-Voce model are show where K is the strengthen coefficient, n is the hardening exponent, and A, B material constants,  represents pre-strain,  denotes the equivalent p By analytical calculation of the load-stroke curve along RD in Figure 2, the strain hardening curve before necking, where the elastic strain is removed from the total strain, is depicted in Figure 4. The maximum plastic strain before necking is approximately 0.17. The strain hardening curve is fitted by the Swift and Voce hardening laws as Equations (1)-(3) respectively, illustrating the hardening behavior with good consistency before necking. Nevertheless, the predicted results by using those two hardening laws are largely different after the necking strain. The Swift model predicts a higher loading capacity of the sheet metal than that obtained from the Voce model. Hence, the combined Swift-Voce model is used to modulate the predicted load capacity characterizing the strain hardening after necking. The parameters of the combined Swift-Voce model are shown in Table 3.
where K is the strengthen coefficient, n is the hardening exponent, and A, B, and C are the material constants, ε 0 represents pre-strain, ε p denotes the equivalent plastic strain.

Drucker Yield Function
It is shown that the strength of materials is strongly affected by stress states [7]. A proper yield function should be selected to predict the effect of stress states on strength Considering the isotropic strength of QP980, the pressure-coupled Drucker (P_Drucker yield function [18] is used to model the plastic behavior under various loading conditions including shear, uniaxial tension, and plane strain tension. Additionally, the P_Drucker yield function couples the pressure effect into the Drucker yield function, which is ex pressed in terms of three stress invariants as: Where I1 is the first stress tensor invariant, J2 and J3 are the second and third deviatoric stress tensor invariant, respectively, a , b , and c are material constants. The applica tion of the P_Drucker function is to precisely describe the strain hardening behavior under various loading conditions. Equation (8) is substituted into Equation (6) to describe the principal stresses as be low: The stress triaxiality  and Lode parameter L are defined as follows:

Drucker Yield Function
It is shown that the strength of materials is strongly affected by stress states [7]. A proper yield function should be selected to predict the effect of stress states on strength.
Considering the isotropic strength of QP980, the pressure-coupled Drucker (P_Drucker) yield function [18] is used to model the plastic behavior under various loading conditions, including shear, uniaxial tension, and plane strain tension. Additionally, the P_Drucker yield function couples the pressure effect into the Drucker yield function, which is expressed in terms of three stress invariants as: where I 1 is the first stress tensor invariant, J 2 and J 3 are the second and third deviatoric stress tensor invariant, respectively, a, b, and c are material constants. The application of the P_Drucker function is to precisely describe the strain hardening behavior under various loading conditions.
Equation (8) is substituted into Equation (6) to describe the principal stresses as below: The stress triaxiality η and Lode parameter L are defined as follows: Based on Equations (8), (10), and (11), the principal stresses can be presented by stress triaxiality η, Lode parameter L, and von Mises effective stress σ as follows: According to Equations (12)- (14), Equation (7) can be transformed as follows: Equation (10) is substituted into Equation (5) in the form of principal stress Based on Equations (9), (15), and (16), Equation (4) can be transformed as a function of η, L, and σ as below: Based on the stress-strain curves obtained by uniaxial tensile tests in Figure 5, the material constant a is determined by Equation (18)

Element Size Sensitivity
The Swift-Voce hardening model and P_Drucker function are programmed into VU-MAT subroutine to predict the deformation processes of QP980 under complex stress states. In order to improve the calculation efficiency, only one-eighth of a dog-bone specimen, a specimen with a central hole and a notched specimen, and a half of an in-plane shear specimen are built as the simulation model. The symmetric boundary conditions are then forced in the simulations with respect to the x-y, y-z and x-z planes for the dog-bone specimen, notched specimen, and the specimen with a central hole. Symmetric boundary condition is applied for the in-plane shear specimen since a half model along the thickness is used for simulation of this test as shown in Figure 6d. Then constant velocity is applied to the fixed end of the specimens according to the tensile velocity in the tests. The mesh type is C3D8R, and the mesh refinement is defined on the deformation concentrated area of the specimens. The element size sensitivity is studied by comparing the experimental and simulated loadstroke curves with six different element sizes for all the specimens, of which the notched specimen is shown in Figure 5. The number of element layers along thickness increases from one to five from set #1 to set #5. The number of elements for each layer is identical to that in set #1. The finite element model of set #6 has five layers along the thickness, but the number of elements in each layer is 1.5 times of set #1. The comparison of the load-stroke curves in Figure 5b indicates that the force response is insensitive to the element size. The local strain evolution at the center of the notched specimen is compared with experimental results. It is observed that element size does not affect the local deformation significantly. Because the parameter optimization by the inverse engineering approach needs to conduct the simulation for dozens of times, the reduction of simulation time is very important for the inverse engineering approach. Consequently, Set #1 is selected as the element model for the notched specimens to improve the numerical simulation efficiency. By using the similar element size analysis, the finite element models for other shaped specimens are determined, as shown in Figure 6. to the element size. The local strain evolution at the center of the notched specimen is compared with experimental results. It is observed that element size does not affect the local deformation significantly. Because the parameter optimization by the inverse engineering approach needs to conduct the simulation for dozens of times, the reduction of simulation time is very important for the inverse engineering approach. Consequently, Set #1 is selected as the element model for the notched specimens to improve the numerical simulation efficiency. By using the similar element size analysis, the finite element models for other shaped specimens are determined, as shown in Figure 6.

Calibration of the Constitutive Model Parameters by Using the Inverse Engineering
The finite element models in Section 4.2 are used to simulate the load-stroke response of the four shaped specimens along RD, DD, and TD, respectively. The combined Swift-Voce hardening law as shown in Equation (3) is utilized to characterize the strain hardening of the material. The von Mises and P_Drucker yield function are used to model the yielding phenomena under various arbitrary stress states. The experimental and simulated load-stroke curves of the four shaped specimens along RD are compared in Figure  7. The load errors predicted by the two yield functions are calculated for simulation along different loading directions according to Equations (19) and (20) as shown in Figure 8. The prediction errors of the two yield functions are approximately similar for the dog-bone specimen and the specimen with a central hole along RD, DD, and TD. For the notched specimens, the results by the P_Drucker yield function is slightly lower than that predicted by the von Mises yield function along DD and TD. The error values of in-plane shear specimens of the P_Drucker function is observed to be much smaller than the von Mises function. Comparing the total error values, it can be seen that the total error values of

Calibration of the Constitutive Model Parameters by Using the Inverse Engineering
The finite element models in Section 4.2 are used to simulate the load-stroke response of the four shaped specimens along RD, DD, and TD, respectively. The combined Swift-Voce hardening law as shown in Equation (3) is utilized to characterize the strain hardening of the material. The von Mises and P_Drucker yield function are used to model the yielding phenomena under various arbitrary stress states. The experimental and simulated load-stroke curves of the four shaped specimens along RD are compared in Figure  7. The load errors predicted by the two yield functions are calculated for simulation along different loading directions according to Equations (19) and (20) as shown in Figure 8. The prediction errors of the two yield functions are approximately similar for the dog-bone specimen and the specimen with a central hole along RD, DD, and TD. For the notched specimens, the results by the P_Drucker yield function is slightly lower than that predicted by the von Mises yield function along DD and TD. The error values of in-plane shear specimens of the P_Drucker function is observed to be much smaller than the von Mises function. Comparing the total error values, it can be seen that the total error values of

Calibration of the Constitutive Model Parameters by Using the Inverse Engineering
The finite element models in Section 4.2 are used to simulate the load-stroke response of the four shaped specimens along RD, DD, and TD, respectively. The combined Swift-Voce hardening law as shown in Equation (3) is utilized to characterize the strain hardening of the material. The von Mises and P_Drucker yield function are used to model the yielding phenomena under various arbitrary stress states. The experimental and simulated loadstroke curves of the four shaped specimens along RD are compared in Figure 7. The load errors predicted by the two yield functions are calculated for simulation along different loading directions according to Equations (19) and (20) as shown in Figure 8. The prediction errors of the two yield functions are approximately similar for the dog-bone specimen and the specimen with a central hole along RD, DD, and TD. For the notched specimens, the results by the P_Drucker yield function is slightly lower than that predicted by the von Mises yield function along DD and TD. The error values of in-plane shear specimens of the P_Drucker function is observed to be much smaller than the von Mises function. Comparing the total error values, it can be seen that the total error values of P_Drucker yield function are all lower than that of von Mises yield function, and the prediction accuracy of P_Drucker yield function is about 15% higher than that of von Mises yield function.
Total err = ∑ err m (20) where n is the number of data point, F Although P_Drucker yield function increases the prediction accuracy of the plastic deformation of QP980 compared to von Mises yield function, the predicted load values by P_Drucker yield function for the specimens with a center hole are higher than the experimental values. However, the predicted load values by P_Drucker yield function for notched specimens are lower than the experimental values. To accurately clarify the material strength of QP980 under complex stress states, the inverse engineering (IE) is used to calibrate the constitutive model parameters. As the model parameter correction method combining experiment and numerical simulation, IE has been widely used in many fields. By constantly optimizing the parameters of the constitutive model by IE, the error between the simulated and the experimental results is minimized according to Equation (19), thus obtaining the optimal parameters of the constitutive model. The optimization algorithm used in IE is Newton downhill method. The optimization process of IE is illustrated in Figure 9. P_Drucker yield function are all lower than that of von Mises yield function, and the prediction accuracy of P_Drucker yield function is about 15% higher than that of von Mises yield function.     is the average value of experimental load. Although P_Drucker yield function increases the prediction accuracy of the deformation of QP980 compared to von Mises yield function, the predicted load by P_Drucker yield function for the specimens with a center hole are higher than perimental values. However, the predicted load values by P_Drucker yield funct notched specimens are lower than the experimental values. To accurately clarify t terial strength of QP980 under complex stress states, the inverse engineering (IE) i to calibrate the constitutive model parameters. As the model parameter correction m combining experiment and numerical simulation, IE has been widely used in many By constantly optimizing the parameters of the constitutive model by IE, the er tween the simulated and the experimental results is minimized according to Eq (19), thus obtaining the optimal parameters of the constitutive model. The optim algorithm used in IE is Newton downhill method. The optimization process of IE i trated in Figure 9. The P_Drucker yield function has better prediction accuracy than the von Mise function. Hence, the parameters of the combined Swift-Voce model and the P_D yield function are calibrated by IE. Considering dog-bone specimen along RD as th mization model of IE, the optimized constitutive model parameters are summari Table 4. The optimized constitutive model parameters are used to simulate the load curves for the specimen with a central hole, the notched specimen, and the in-plane specimen. Figure 10 presents the predicted load capability after optimization with fitting to the experimental values, decreasing the deviation of the P_Drucker yield fu in the numerical simulation of the specimen with a center hole and the notched spe The error comparison in Figure 11 demonstrates that the constitutive model optimi IE can accurately characterize the strain hardening behavior of QP980.  The P_Drucker yield function has better prediction accuracy than the von Mises yield function. Hence, the parameters of the combined Swift-Voce model and the P_Drucker yield function are calibrated by IE. Considering dog-bone specimen along RD as the optimization model of IE, the optimized constitutive model parameters are summarized in Table 4. The optimized constitutive model parameters are used to simulate the load-stroke curves for the specimen with a central hole, the notched specimen, and the in-plane shear specimen. Figure 10 presents the predicted load capability after optimization with better fitting to the experimental values, decreasing the deviation of the P_Drucker yield function in the numerical simulation of the specimen with a center hole and the notched specimen. The error comparison in Figure 11 demonstrates that the constitutive model optimized by IE can accurately characterize the strain hardening behavior of QP980. The maximum value and distribution of the equivalent plastic strain obtained from the simulation by using the optimized constitutive model are close to those measured by DIC, as shown in Figure 12. On the basis of the equivalent plastic strain distribution, the equivalent plastic strain of Node A in numerical simulation is extracted and then compared with DIC measurement, as shown in Figure 13. The error between the simulated and DIC measured equivalent plastic strains of the notched specimen, as shown in Figure 13b, is found to be less than 5%. The comparison further reveals that the strain hardening behavior of those specimens can be accurately characterized by using the combined Swift-Voce hardening model and P_Drucker yield function, and the IE can effectively enhance the prediction accuracy.  The maximum value and distribution of the equivalent plastic strain obtained from the simulation by using the optimized constitutive model are close to those measured by DIC, as shown in Figure 12. On the basis of the equivalent plastic strain distribution, the equivalent plastic strain of Node A in numerical simulation is extracted and then compared with DIC measurement, as shown in Figure 13. The error between the simulated and DIC measured equivalent plastic strains of the notched specimen, as shown in Figure  13b, is found to be less than 5%. The comparison further reveals that the strain hardening behavior of those specimens can be accurately characterized by using the combined Swift-Voce hardening model and P_Drucker yield function, and the IE can effectively enhance the prediction accuracy.   The maximum value and distribution of the equivalent plastic strain the simulation by using the optimized constitutive model are close to those DIC, as shown in Figure 12. On the basis of the equivalent plastic strain di equivalent plastic strain of Node A in numerical simulation is extracted a pared with DIC measurement, as shown in Figure 13. The error between Figure 11. Comparison of the prediction error value before and after parameter optimization.
According to the experimental load-stroke curves shown in Figure 13a-c, the stroke that drop of load capability is clearly observed is defined as the fracture stroke, and the equivalent strain at the fracture stroke is defined as the fracture strain that is marked by blue pentagram. The ductile fracture strains of the three shaped specimens are summarized in Table 5.  blue pentagram. The ductile fracture strains of the three shaped specimens are summarized in Table 5.  According to the experimental load-stroke curves shown in Figure 13a-c, the stroke that drop of load capability is clearly observed is defined as the fracture stroke, and the equivalent strain at the fracture stroke is defined as the fracture strain that is marked by blue pentagram. The ductile fracture strains of the three shaped specimens are summarized in Table 5.

Fracture Prediction
The accurate prediction of the fracture in sheet metal forming requires an appropriate ductile fracture criterion as well as a proper parameter calibration. The uncoupled ductile fracture criterion is widely used to predict the ductile fracture of sheet metals during the forming processes. However, the traditional uncoupled ductile fracture criteria generally do not take the effect of the stress states into account, resulting in a consequence that the prediction accuracy is not promising in a wide range of stress triaxiality [19]. In this research, the fracture is modeled by a shear-controlled ductile fracture criterion which is proposed based on the micro-mechanism of the ductile fracture for nucleation, growth, and coalescence of voids [3]. The ductile fracture criterion is referred to as DF2012 as below: where η is stress triaxiality and τ max is maximum shear stress. The traditional uncoupled ductile fracture criteria used in this section are summarized as below: Ko-Huh criterion [20]: Brozzo criterion [21]: Oh criterion [22]: Rice-Tracey criterion [23]: DF2012 and the traditional uncoupled ductile fracture criteria, namely Ko-Huh criterion, Brozzo criterion, Oh criterion, and Rice-Tracey criterion, are respectively used to predict the ductile fracture of QP980 under the complex stress states. The stress triaxiality, Lode parameters, and equivalent plastic strain from the initiation of plastic strain to fracture in the numerical simulation are shown in Figure 14a-c, the extraction positions of these parameters are depicted in Node A of Figure 13. It is observed that the stress state is approximately proportional for the specimen with a central hole. There is a slight rise in the Lode parameter. The stress state evolution is very strong for the in-plane shear specimen. It is a big challenge to characterize the plasticity and fracture behavior under shear. Different experimental methods were developed to characterize mechanical behavior under shear for sheet metals, such as torsion tests. Zhang et al. [7] combined the in-plane torsion test with an inverse engineering approach to characterize shear flow curves up to large strain. Even though the torsion tests perform better than in-plane shear specimen to characterize mechanical properties under shear, the specimen geometry and test machine are complicated for the torsion test. Therefore, the in-plane shear specimen is the simplest method to characterize the shear properties of sheet metals by experiments. The average stress triaxiality and Lode parameter of three shaped specimens are calculated by Equations (26) and (27) and then summarized in Table 5. The results in Table 5 are used to calibrate the material constants in the ductile fracture criterion, denoted as Table 6. Based on the optimized constitutive model of Section 4.3, those ductile fracture criteria are conducted to predict the ductile fracture of QP980 under complex stress states. The minimum difference between the experimental and the simulated fracture stroke are employed to calibrate the ductile fracture criterion parameters by using IE, as shown in Table 6. The comparison of the stroke-load curves between the experiments and simulations is plotted in Figure 15a-c. The prediction accuracy of the DF2102 criterion is significantly higher than that of the traditional uncoupled ductile fracture criteria. The prediction errors of the different fracture criteria are expressed in Figure  16. The prediction error of DF2012 by IE is the lowest, compared to the traditional uncou- The results in Table 5 are used to calibrate the material constants in the ductile fracture criterion, denoted as Table 6. Based on the optimized constitutive model of Section 4.3, those ductile fracture criteria are conducted to predict the ductile fracture of QP980 under complex stress states. The minimum difference between the experimental and the simulated fracture stroke are employed to calibrate the ductile fracture criterion parameters by using IE, as shown in Table 6. The comparison of the stroke-load curves between the experiments and simulations is plotted in Figure 15a-c. The prediction accuracy of the DF2102 criterion is significantly higher than that of the traditional uncoupled ductile fracture criteria. The prediction errors of the different fracture criteria are expressed in Figure 16. The prediction error of DF2012 by IE is the lowest, compared to the traditional uncoupled ductile fracture criteria. This is because the influences of the hydrostatic stress, equivalent stress and maximum principal stress are considered in the traditional uncoupled ductile fracture criteria, and those mechanical variables are closely related to stress triaxiality. Bao et al. [24] found that the fracture mechanisms depends on the stress triaxiality ranges. The influence of the stress triaxiality and normalized maximum shear stress (τ max /σ) are considered in DF2012. Hence, DF2012 can effectively predict the ductile fracture of QP980 during the forming processes. The fracture locus of QP980 predicted by DF2012 and the traditional uncoupled ductile fracture criteria is depicted in Figure 17. The comparison in Figure 17b shows that DF2012 can describe the fracture behavior of QP980 in a wide range of stress triaxiality, and the ductile fracture of sheet metal needs to consider the influence of stress triaxiality and Lode parameters. DF2012 and the traditional uncoupled ductile fracture criteria is depicted in Figure 17. T comparison in Figure 17b shows that DF2012 can describe the fracture behavior of QP in a wide range of stress triaxiality, and the ductile fracture of sheet metal needs to c sider the influence of stress triaxiality and Lode parameters.

Conclusions
This research characterizes the failure of QP980 sheet metals under wide stress states from shear to plane strain tension by experiments, analytical modeling, and numerical simulation. It is found that the four specimens used in this study are proper to measure the fracture strain under different stress states of shear, uniaxial tension, and plane strain tension. The P_Drucker yield function with the Swift-Voce hardening law is shown to be capable of precise modeling of plastic response of QP980 from the onset of plastic deformation to fracture at different stress states. The DF2012 criterion is observed to predict the onset of ductile fracture with good agreement compared with experiments. The predicting accuracy of the constitutive models is strongly dependent on the calibration of the material constants in the constitutive models. The strain hardening law, yield function, and fracture criterion are recommended to be calibrated by the inverse engineering approach to improve the performance of the constitutive models and the reliability of numerical simulations.
Author Contributions: Writing-original draft preparation, formal analysis, S.W.; methodology, Figure 17. The fracture locus of QP980: (a) in the space of (η, L, ε f ) predicted by DF2012; (b) in the space of (η, ε f ) predicted by the traditional uncoupled ductile fracture criteria and DF2012.

Conclusions
This research characterizes the failure of QP980 sheet metals under wide stress states from shear to plane strain tension by experiments, analytical modeling, and numerical simulation. It is found that the four specimens used in this study are proper to measure the fracture strain under different stress states of shear, uniaxial tension, and plane strain tension. The P_Drucker yield function with the Swift-Voce hardening law is shown to be capable of precise modeling of plastic response of QP980 from the onset of plastic deformation to fracture at different stress states. The DF2012 criterion is observed to predict the onset of ductile fracture with good agreement compared with experiments. The predicting accuracy of the constitutive models is strongly dependent on the calibration of the material constants in the constitutive models. The strain hardening law, yield function, and fracture criterion are recommended to be calibrated by the inverse engineering approach to improve the performance of the constitutive models and the reliability of numerical simulations.
Author Contributions: Writing-original draft preparation, formal analysis, S.W.; methodology, software, validation, data curation, Y.W.; writing-review and editing, L.Y.; writing-review and editing, K.J.; writing-review and editing, X.L.; conceptualization, supervision, writing-review and editing, project administration, Y.L. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.