A coupled method between the Scaled Boundary Finite Element Method (SBFEM) and Finite Element Method (FEM) for evaluating the Stress Intensity Factors (SIFs) is presented and achieved on the platform of the commercial finite element software ABAQUS by using Python as the programming language. Automatic transformation of the finite elements around a singular point to a scaled boundary finite element subdomain is realized. This method combines the high accuracy of the SBFEM in computing the SIFs with the ability to handle material nonlinearity as well as powerful mesh generation and post processing ability of commercial FEM software. The validity and accuracy of the method is verified by analysis of several benchmark problems. The coupled algorithm shows a good converging performance, and with minimum additional treatment can be able to handle more problems that cannot be solved by either SBFEM or FEM itself. For fracture problems, it proposes an efficient way to represent stress singularity for problems with complex geometry, loading condition or certain nonlinearity.
An open-source ABAQUS implementation of the scaled boundary finite element method to study interfacial problems using polyhedral meshes
