Transient analysis of the dynamic stress intensity factors using SGBEM for frequency-domain elastodynamics

In this paper, a two-dimensional symmetric-Galerkin boundary integral formulation for elastodynamic fracture analysis in the frequency domain is described. The numerical implementation is carried out with quadratic elements, allowing the use of an improved quarter-point element for accurately determining frequency responses of the dynamic stress intensity factors (DSIFs). To deal with singular and hypersingular integrals, the formulation is decomposed into two parts: the first part is identical to that for elastostatics while the second part contains at most logarithmic singularities. The treatment of the elastostatic singular and hypersingular singular integrals employs an exterior limit to the boundary, while the weakly singular integrals in the second part are handled by Gauss quadrature. Time histories (transient responses) of the DSIFs can be obtained in a post-processing step by applying the standard fast Fourier transform (FFT) and algorithm to the frequency responses of these DSIFs. Several test examples are presented for the calculation of the DSIFs due to two types of impact loading: Heaviside step loading and blast loading. The results suggest that the combination of the symmetric-Galerkin boundary element method and standard FFT algorithms in determining transient responses of the DSIFs is a robust and effective technique.