Abstract:Based on the basic idea solving multi-domain combinations with virtual boundary element least square method, that each crack can be treated as a pair of sub-domains, and learning from the interpolation of the compactly supported radial basis function used in boundary-type meshless methods to approximately construct the virtual source function on the virtual boundary corresponding to each sub-domain, the computational scheme with the virtual boundary meshless least squares analyzing two-dimensional multi-crack problems is established. According to the definition about sub-domain in this paper, the added extra sub-domains on the boundary extended along the crack surface as “conventional sub-domain method” in the direct boundary element method do not have to be considered, thereby reducing the computational, especially avoiding this calculation error caused due to inadequate number of the elements or with the collocation points configured on the boundary of the additional sub-domains and its improper configuration. In order to verify feasibility and accuracy of the numerical algorithm proposed in the article and discuss the interaction between multiple cracks with arbitrary distribution, some examples, such as a single center crack, third-adjacent spacing of different length collinear cracks and one horizontal crack and one different inclined angle crack by one-way tension in infinite plate, are given. The results show that this method leads to higher accuracy in comparison with the other methods considered in this study.