Automated simulation of fatigue crack propagation for two-dimensional linear elastic fracture mechanics problems by boundary element method

In this paper, automated simulation of multiple crack fatigue propagation for two-dimensional (2D) linear elastic fracture mechanics (LEFM) problems is developed by using boundary element method (BEM). The boundary element method is the displacement discontinuity method with crack-tip elements proposed by the author. Because of an intrinsic feature of the boundary element method, a general growth problem of multiple cracks can be solved in a single-region formulation. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not necessary. Local discretization on the incremental crack extension is performed easily. Further the new adding elements and the existing elements on the existing boundaries are employed to construct easily the total structural mesh representation. Here, the mixed-mode stress intensity factors are calculated by using the formulas based on the displacement fields around crack tip. The maximum circumferential stress theory is used to predict crack stability and direction of propagation at each step. The well-known Paris’ equation is extended to multiple crack case under mixed-mode loadings. Also, the user does not need to provide a desired crack length increment at the beginning of each simulation. The numerical examples are included to illustrate the validation of the numerical approach for fatigue growth simulation of multiple cracks for 2D LEFM problems.