An efficient numerical model for contact-induced crack propagation analysis

This paper presents an efficient boundary element formulation for the simulation of contact-induced crack propagation problems. The formulation combines a load-scaling boundary element technique for frictional contact analysis and a dual boundary element method for axisymmetric crack propagation. These previously developed contact and fracture mechanics algorithms are revisited and adapted for the proposed coupled formulation. The novelty of the present work is the development of an efficient algorithm capable of simulating contact-induced crack propagation, while keeping proper equilibrium conditions in the frictional contact area. An efficient combination of the two separate formulations is not straightforward, and possible approaches are explained. The algorithm is validated by comparing the results of several indentation tests with available experimental data and analytical expressions. Its efficiency is also verified in frictional, conforming and non-conforming, indentation simulations.