A FAST MUTIPOLE DUAL BEM FOR THE ANALYSIS OF 2-D CRACKS

Z. H. Yao,P. B. Wang
2005-01-01
Abstract:Fast mutipole method (FMM) is applied to Dual BEM for 2-D elastic crack analysis. Dual BEM overcomes the mathematical degeneration of the displacement boundary integral equation when the two surfaces of the same crack are co-planar, by introduction of the traction boundary integral equation. The concept of finite-part integral is applied to deal with the hyper-singular integrals in traction boundary integral equation. Crack surfaces are discretized using discontinuous elements to satisfy continuity requirements. In order to achieve more run-time and space efficiency, fast mutipole method (FMM) is applied to Dual BEM. Due to the realization of mutipole expansion and local expansion, both computational complexity and memory requirement is reduced to O(N), where N is the number of DOF . In mutipole expansion, a new form of complex Taylor series reduces the number of mutipole moments and local moments. Generalized minimum residual method (GMRES) is adopted as iterative solver. Sparse approximate inverse type is selected as the left preconditioner. An effective sparsity pattern is developed to reduce the cost of preconditioni ng and improve convergence rate. Two numerical examples are given to verify the numerical method and corresponding program. The numerical results of normalized SIF for a center crack in square plate show good agreement with that in literature. The results of COD for four cracks in square plate show good agreement with FEM results of MSC/MARC. Finally a large -scale numerical example of 3500 irregularly distributed cracks in square plate, with 1,054,824 DOF, is given. The numerical results show that the presented fast mutipole Dual BEM scheme is very efficient for large-scale crack problems.
What problem does this paper attempt to address?