This paper extends the scatterer polymerization methodology to three-dimensional multiple scattering of elastic waves by spherical inclusions. The methodology was originally developed for analyzing multiple scattering of elastic antiplane shear waves in two-dimensional spaces. The analytically exact solution of multiple scattering is reformulated by using this methodology, which is verified by using different ways, with or without scatterer polymerization, to solve physically the same multiple scattering problem. As an application example, the band gap formation for elastic wave propagating in a cubic lattice of spherical scatterers is observed through a series of numerical simulations. These simulations also demonstrate the capability of the present computational system for simulating three-dimensional multiple scattering of elastic waves.