A two-level nesting smoothed extended meshfree method for static and dynamic fracture mechanics analysis of orthotropic materials

This paper presents a two-level nesting smoothed extended mesh-free method (S-XMM-N) for the static and dynamic fracture mechanics analysis in orthotropic materials. To capture the jump of displacement fields across the crack surface, as well as singularities of stress fields in the vicinity of the crack tip, both Heaviside function and asymptotic solutions from the orthotropic fracture modeling are incorporated into the radial point interpolation shape functions associated with the partition of unity approach. To improve the accuracy and accelerate the computation of the extended mesh-free method, the two-level nesting smoothing integration scheme based on the first and second-level triangular smoothing sub-domains is developed to perform the stiffness integration. Compared to the extended mesh-free method using the conventional Gaussian quadrature, the presented S-XMM-N uses even less integration sampling points. At the same time, the complex and time-consuming derivative calculation of extended shape functions is completely avoided. This leads to the dramatically improvement of efficiency. Moreover, the presented method gives very exactly numerical solutions by optimally combining the contributions from the two-level nesting smoothing sub-domains based on the Richardson extrapolation method. We use the interaction integral method in conjunction with the asymptotic near crack tip fields of orthotropic materials to extract the static and dynamic stress intensity factors. The numerical solutions obtained from the S-XMM-N are further compared with reference solutions from the literature to verify the accuracy of the proposed method.