Dual boundary element method and finite element method for mixed-mode crack propagation simulations in a cracked hollow shaft

Three-dimensional mixed-mode crack propagation simulations were performed by means of the dual boundary element method code BEASY and 2 finite element method-based crack propagation codes: ZENCRACK (ZC) and CRACKTRACER3D (CT3D). The stress intensity factors (SIFs) along the front of an initial semielliptical crack, initiated from the external surface of a shaft, were calculated for 4 different load cases: bending, press fit, shear, and torsion. The methods used for the SIF assessment along the crack front were the J-integral for BEASY and ZC and the quarter point element stress method for CT3D. Subsequently, crack propagation simulations were performed, with the crack growth rate evaluated by using Paris' law, calibrated for the material at stake (American Society for Testing and Materials A469 steel). The kink angles were evaluated by using the minimum strain energy density and maximum tangential stress criteria for BEASY, the maximum energy release rate and maximum tangential stress for ZC, and the maximum principal asymptotic stress for CT3D. The results obtained in terms of SIFs and crack propagation life show very good agreement among the 3 codes. Also, the shape of the propagated crack, which is significantly out-of-plane for the shear and torsion loading, matched very well.