Fatigue crack growth simulation in heterogeneous material using s-version FEM

•A fatigue crack growth simulation system in a heterogeneous material is developed.•The system is fully automated and employs the s version Finite Element Method (s-FEM).•The system gives an adequate accurate estimation of the stress intensity factors.•We showed the difference in Young’s modulus has a large effect on the crack path.•s-FEM successfully simulates the crack growth process in a heterogeneous material.

A fully automatic fatigue crack growth simulation system is developed using the s-version Finite Element Method (s-FEM). This system is extended to fractures in heterogeneous materials. In a heterogeneous material, the crack tip stress field has a mixed-mode condition, and the crack growth path is affected by inhomogeneous materials and mixed-mode conditions. Stress intensity factors (SIFs) in the mixed-mode condition are evaluated using the virtual crack closure method (VCCM). The criteria for the crack growth amount and crack growth path are based on these SIFs, and the growing crack configurations are obtained. At first, the basic problem is solved, and the results are compared with some results available in the literature. It is shown that this system gives an adequate accurate estimation of the SIFs. Then, two-dimensional fatigue crack growth problems are simulated using this system. The first example is a plate with an interface between hard and soft materials. The cracks tend to grow in soft materials through the interface. A second example is a plate with distributed hard inclusions. The crack takes a zig-zag path by propagating around the hard inclusions. In each case, the crack growth path changes in a complicated manner. Changes of the SIFs values are also shown and discussed.