An effective numerical stress intensity factor calculation with no crack discretization

An effective numerical procedure for calculating stress intensity factors (SIF) in plane problems based on a modified boundary element technique not requiring any crack discretization was proposed by Snyder 1]. Instead of the usual fundamental solution, he used Green's function for the problem of a traction-free central crack in an anisotropic plate.In the first part of the present paper, the corresponding Green's function for the isotropic problem, not explicitly included in 1], is presented. In addition to the central crack, a semi-infinite edge crack is considered. Both Green's functions are given for the case of the anti-plane state of strain as well. In the first step of the proposed procedure, the tractions and displacements along the outer boundary are calculated. In the second step, the SIF for modes I, II and III are derived in terms of simple boundary integrals over quantities known from the previous step. Contrary to Snyder's derivation, the determination of the SIF is based on the asymptotic displacement field at the crack tip. The method can easily be extended to multiple crack problems by using the subregion technique. Some illustrative examples demonstrate the effectiveness of the method.

---

Résumé Snyder a proposé une procédure numérique pour le calcul des facteurs et intensité de contraintes dans les problèmes plans, en se basant sur une modification de la technique des éléments aux limites ne requérant pas de discrétisation de la fissure. Au lieu d'une solution fondamentale habituelle, il a utilisé une fonction de Green pour traiter le problème d'une fissure centrale libre de contraintes dans un plaque anisotrope.Dans la première partie de la présente étude, on présente la fonction de Green qui correspond à un problème isotrope, qui n'était pas explicitement couvert par le travail de Snyder. Outre le cas de la fissure centrale, on considère le cas de la fissure de bord dans un milieu semi-infini. Les deux fonctions de Green relatives au cas d'un état de déformation antiplanaire sont également communiquées. Dans une première étape de la procédure proposée, on calcule les sollicitations et déplacements le long du contour extérieur; dans une deuxième étape, on établit les facteurs d'intensité de contraintes relatifs aux modes I, II et III en terms d'intégrals sur un contour simple défini par les valeurs résultant de l'étape précédente. Contraitement à l'approche de Snyder, on détermine le facteur d'intensité des contraintes sur base du champs de déplacement asymptotique à l'extrémité de la fissure. Par la technique des sous-régions, on peut aisément étendre la méthode à des problèmes de fissuration multiples. L'efficacité de la méthode est illustrée par divers exemples.