The present work considers the use of the two-fluid (Euler-Euler) CFD approach for the continuum description of vibrated fluidized beds as a less computationally demanding alternative to the discrete description given by Lagrangian-Eulerian methods such as DEM. In particular, a novel simulation strategy consisting on solving the two-fluid model equations in a coordinate reference system that moves with the vibrating walls of a gas-solid fluidized bed is proposed. By this way, vibration is transformed into simple alternating acceleration terms that are introduced through body forces in both the gas and the particle phase equations. The results of a series of two-fluid model simulations compare well with discrete particle simulations as well as with experimental data reported for beds containing Geldart group B particles. In general, the results of a series of two-fluid model simulations show similar trends to those seen in discrete particle simulations as well as in experimental data reported for beds containing Geldart group B particles. Exception of that is the velocity of bubbles, for which the two-fluid simulations compare less satisfactorily with the available experimental data. The two-fluid model simulations are also able to reproduce expected phenomena like the bubble growth with the vibration amplitude and the dependence of the pressure drop fluctuation on the vibration strength. In view of these promising results, the proposed two-fluid model formulation opens the possibility of increasing the scale of the vibrated fluidized beds currently simulated. (C) 2012 Elsevier B.V. All rights reserved.