A new crack-tip element for the logarithmic stress-singularity of Mode-III cracks in spring interfaces

A new crack-tip finite element able to improve the accuracy of Finite Element Method (FEM) solutions for cracks growing along the Winkler-type spring interfaces between linear elastic adherents is proposed. The spring model for interface fracture, sometimes called Linear-Elastic (perfectly) Brittle Interface Model (LEBIM), can be used, e.g., to analyse fracture of adhesive joints with a thin adhesive layer. Recently an analytical expression for the asymptotic elastic solution with logarithmic stress-singularity at the interface crack tip considering spring-like interface behaviour under fracture Mode III was deduced by some of the authors. Based on this asymptotic solution, a special 5-node triangular crack-tip finite element is developed. The generated special singular shape functions reproduce the radial behaviour of the first main term and shadow terms of the asymptotic solution. This special element implemented in a FEM code written in Matlab has successfully passed various patch tests with spring boundary conditions. The new element allows to model cracks in spring interfaces without the need of using excessively refined FEM meshes, which is one of the current disadvantages in the use of LEBIM when stiff spring interfaces are considered. Numerical tests carried out by h-refinement of uniform meshes show that the new singular element consistently provides significantly more accurate results than the standard finite elements, especially for stiff interfaces, which could be relevant for practical applications minimizing computational costs. The new element can also be used to solve other problems with logarithmic stress-singularities.