Based on the 2D BEM software THBEM2 developed by the authors' group, which can be applied to simulate the elastic body with randomly distributed circular holes. A scheme of BEM for the simulation of elastic bodies with randomly distributed identical circular inclusions is presented. The workload of the computation is nearly the same as for the corresponding case of circular holes, the numerical examples given show that the scheme is applicable and effective.