Modeling layered systems with dislocations is very challenging; yet, it is important since most smart structures are made of multilayers to make best use of the combined effective property. As such, during the manufactures, defects, such as dislocations, could be introduced in the multilayers. In this article, we analytically find, for the first time, the response of three-dimensional multilayered magneto-electro-elastic systems due to time-harmonic dislocations. The dislocations are the most general, containing the elastic dislocations and discontinuity of the electric potential and/or magnetic potential over a circular region in any layer in the medium. The fully coupled partial differential equations of motion and the Gauss law for the magneto-electro-elastic materials are solved in terms of cylindrical system of vector functions, and the dual variable and position method is further introduced to treat the multilayers. Numerical examples are carried out based on the derived analytical solution to demonstrate the effects of the time-harmonic dislocations on the induced magneto-electro-elastic fields. This analytical solution is important in both electrodynamics and elastodynamics, with possible applications in material sciences and physics. The numerical results are useful in design process of smart devices made of magneto-electro-elastic solids applicable to other engineering fields like renewable energy.