Metamaterials, i.e. artificial structures with unconventional properties, have shown to be highly potential lightweight and compact solutions for the attenuation of noise and vibrations in targeted frequency ranges, called stop bands. In order to analyze the performance of these metamaterials, their stop band behavior is typically predicted by means of dispersion curves, which describe the wave propagation in the corresponding infinite periodic structure. The input for these calculations is usually a finite element model of the corresponding unit cell. Most common in literature are 2D plane metamaterials, which often consist of a plate host structure with periodically added masses or resonators. In recent literature, however, full 3D metamaterials are encountered which are periodic in all three directions and which enable complete, omnidirectional stop bands. Although these 3D metamaterials have favorable vibro-acoustic characteristics, the computational cost to analyze them quickly increases with unit cell model size. Model order reduction techniques are important enablers to overcome this problem. In this work, the Bloch Mode Synthesis (BMS) and generalized BMS (GBMS) reduction techniques are extended from 2D to 3D periodic structures. Through several verifications, it is demonstrated that dispersion curve calculation times can be strongly reduced, while accurate stop band predictions are maintained.