This paper discusses a fast multipole method for three dimensional anisotropic elastodynamics in time domain. A PWTD (plane wave time domain) approach is presented. The proposed approach is verified via a simple numerical example. The accuracy of the approach is shown to be satisfactory for engineering purposes, although its efficiency still needs to be improved.