Comparison of Various Criteria Determining the Direction of Crack Propagation Using the UDMGINI User Procedure Implemented in Abaqus.

Abstract

This paper describes a method of predicting the direction of crack propagation implemented by user subroutines in the Simulia-Abaqus FEA system with the use of the extended finite element method (X-FEM). This method is based on displacements and stresses according to Westergaard’s solution of Griffith’s crack problem. During the calculations, in each crack increment, the algorithm reads the stresses and displacements in the model around the crack tip, calculates the criterion values at the read points, reduces them to a unit distance from the crack tip, fits a polynomial to these points, and finds the minimum of the function closest to the last propagation angle. The algorithm also decides when the crack grows, depending on a chosen criterion. Four criteria have been implemented to predict the direction of failure propagation: the maximum principal stress criterion, the Ottosen–Podgórski criterion, the new criterion described here based on the minimum component values of the displacement vector, and the maximum circumferential tensile stress (MTS). These criteria were verified in two tests: the three-point bending test of the notched beam and the anchor pull-out test. For these tests, the criterion built into Simulia Abaqus does not correctly define the crack path, which causes the crack propagation direction to “rotate” when simulating the fracture. The criteria developed here, in most cases, determine the crack path and the maximum force very well compared to real laboratory tests.