Stress analysis of inclusion problems of various shapes in an infinite anisotropic elastic medium

A boundary integral equation approach is used to solve an infinite anisotropic elastic inclusion problem subjected to remote loading. Continuous and discontinuous quadratic isoparametric boundary elements are employed to model the interfaces between the inclusions and the matrix. The inclusion–matrix interfaces are assumed to be perfectly bonded. Inclusions of various shapes and their interaction are investigated. Numerical examples are compared with existing analytical solutions. Relative to the finite element method and the volume integral method, the present method is more efficient and accurate in the analysis of elastic inclusion problems.