We propose an isogeometric boundary element (IGABEM) parallel fast direct solver based on OpenMP to solve three-dimensional multiple crack problems. In this paper, the orderly and randomly distributed cracks are studied to observe the influence of the multiple crack interaction on the effective elastic properties. This paper designs the parallel implementation of the matrix low rank approximation and the matrix assembly for the hierarchical off-diagonal low rank (HODLR) matrix. In the matrix low rank approximation, we propose the m-n tree to improve the accelerated hybrid algorithm to balance the CPU performance and memory usage. The numerical results show that the effective elastic properties can respond to the strong interaction to a certain extent, but the maximum stress intensity factors (SIFs) of multiple cracks still need to be calculated to diagnose the possible crack initiation points. On the other hand, the numerical results also show that the IGABEM parallel fast direct solver can greatly improve the computational efficiency on the premise of ensuring the computational accuracy.
Read full abstract