This paper presents simulation of 3D composite materials containing a large number of particles using a new version of adaptive fast multipole boundary element method (BEM). A scheme of similar sub-domain approach is also applied for the case of identical sphere particles. Generalized Minimum Residual Method (GMRES) is selected as an iterative solver for the equation systems implicitly arising from fast multipole BEM. The number of particles reaches 100 by the combination of above algorithms on an personal computer. Numerical results show that fast multipole BEM is applicable for large scale simulation of certain composite materials.