The influence of anisotropic properties and layering of saturated soils on dynamic response of a circular foundation under time-harmonic vertical loading is studied in this paper. The foundation is modeled as a rigid disk, with either permeable or impermeable and smooth contact surface, whereas a multilayered transversely isotropic poroelastic half-space is defined for the supporting soil. To investigate the dynamic interaction problem, a displacement and flow boundary condition at the contact surface is constituted by employing a discretization technique and the displacement influence functions, which are obtained from the exact stiffness matrix scheme, with the stiffness matrices being presented explicitly for the first time. A set of numerical results are presented to portray the effects of anisotropic and poroelastic properties, hydraulic boundary conditions, soil heterogeneity, and frequency of excitation on vertical dynamic response of rigid circular foundations. In addition, vertical vibrations of a ring foundation and a massive circular foundation are also investigated.