Simulation of ductile crack propagation in dual-phase steel

The modified Mohr–Coulomb and the extended Cockcroft–Latham fracture criteria are used in explicit finite-element (FE) simulations of ductile crack propagation in a dual-phase steel sheet. The sheet is discretized using tri-linear solid elements and the element erosion technique is used to model the crack propagation. The numerical results are compared to quasi-static experiments conducted with five types of specimens (uniaxial tension, plane-strain tension, in-plane shear, 45° and 90° modified Arcan) made from a 2 mm thick sheet of the dual-phase steel Docol 600DL. The rate-dependent J 2 flow theory with isotropic hardening was used in the simulations. The predicted crack paths and the force–displacement curves were quite similar in the simulations with the different fracture criteria. Except for the 45° modified Arcan test, the predicted crack paths were in good agreement with the experimental findings. The effect of using a high-exponent yield function in the prediction of the crack path was also investigated, and it was found that this improved the crack path prediction for the 45° modified Arcan test. In simulations carried out on FE models with a denser spatial discretization, the prediction of localized necking and crack propagation was in better accordance with the experimental observations. In four out of five specimen geometries, a through-thickness shear fracture was observed in the experiments. By introducing strain softening in the material model and applying a dense spatial discretization, the slant fracture mode was captured in the numerical models. This did not give a significant change in the global behaviour as represented by the force–displacement curves.