FFT-based solvers are increasingly used by many researcher groups interested in modelling the mechanical behavior associated to a heterogeneous microstructure. A development is reported here that concerns the viscoelastic behavior of composite structures generally studied experimentally through Dynamic Mechanical Analysis (DMA). A parallelized computation code developed with complex-valued quantities provides virtual DMA experiments directly in the frequency domain on a heterogeneous system described by a voxel grid of mechanical properties. The achieved precision and computation times are very good. An effort has been made to illustrate the application of such a virtual DMA tool through two examples from the literature: the modelling of glassy/amorphous systems at a small scale and the modelling of experimental data obtained in temperature sweeping mode by DMA on a particulate composite made of glass beads and a polystyrene matrix, at a larger scale. Both examples show how virtual DMA can contribute to question, analyze, and understand relaxation phenomena on either theoretical or experimental points of view.