The dynamic response of rigid and flexible three-dimensional (3D) foundations of arbitrary shape, placed on a homogeneous or layered halfspace, which may contain a finite, longitudinally invariant structure or inhomogeneity, is determined numerically. A 2.5D coupled Integral Transform Method (ITM) - Finite Element Method (FEM) approach is used to compute the dynamic compliance at the surface of the stratified soil with inclusion, whereas 3D ITM fundamental solutions are applied in case of a homogeneous or layered subgrade. The foundation is modelled by 3D finite elements in both cases and coupled to the underlying ground by enforcing the compatibility conditions at the interface. Compliance functions at the soil foundation interface are pre...