3D predictions of the local effective stress intensity factor as the fatigue crack propagation driving force

Abstract The objective of this study is to develop a numerical tool using the commercial software, Abaqus and Python, to predict the fatigue crack front shape while taking into account the influence of plasticity induced crack closure on crack propagation in three Dimensional (3D) structures. In this aim, a 3D model of a compact tension specimen made out of stainless steel 304L, and subjected to a constant loading scheme, is proposed. The crack propagation is considered to be driven by the stress fields developed in the vicinity of the crack tip and thus by the stress intensity factor K. Two parallel simulations are used: an elastic simulation intends to calculate the local maximum stress intensity factor while the other, an elasto-plastic one, aims at obtaining the plastic wake and the resulting crack closure load. The results of both simulations are combined in order to constitute the effective stress intensity factor range, which is in turn used, along with Paris law, to calculate the crack propagation along the thickness. The local crack advancements obtained allow to construct the new crack front. Finally, a node release technique is used with geometry remeshing to issue new iterations with new boundary conditions that respond to the changes in the crack front. The procedure is repeated until the stabilization of the effective stress intensity factor values all along the specimen thickness is reached. The results obtained are compared with previously issued experimental results, showing very good results in small scale yielding and beyond that a large dependency on the plastic zone size developed in the neighborhood of the crack front.