In this work, some solutions such as Mori-Tanaka approximation (MTA), Differential approximations (DA), and Fast Fourier transformation method (FFT) were applied to estimate the elastic bulk shear modulus of three-phase composites in 2D. which two different sizes circular inclusions are arranged randomly non-overlapping a continuous matrix. The numerical using FFT analysis compared with DA, MTA...