The Havriliak-Negami model for dynamic viscoelastic material behavior and Biot’s theory of poroelasticity are employed to develop an exact solution for three-dimensional scattering effect of harmonic plane P–SV waves from a circular cavity lined with a multilayered fluid-filled shell of infinite length containing viscoelastic damping materials and embedded within a fluid-saturated permeable surrounding soil medium. The analytical results are illustrated with numerical examples where the effects of liner/coating structural arrangement, viscoelastic material properties, liner-soil interface bonding condition, seismic excitation frequency, and angle of incidence on the induced dynamic stress concentrations are evaluated and discussed to obtain representative values of the parameters that characterize the system. It is demonstrated that incorporating viscoelastic damping materials with a low shear modulus in the constrained layer configuration is an efficient means of enhancing the overall seismic isolation performance, especially for near-normally incident seismic shear waves where the amplitudes of induced dynamic stresses may be reduced by up to one-third of those without isolation in a relatively wide frequency range. Some additional cases are considered and good agreements with solutions available in the literature are obtained. The Havriliak-Negami model for dynamic viscoelastic material behavior and Biot’s theory of poroelasticity are employed to develop an exact solution for three-dimensional scattering effect of harmonic plane P-SV waves from a circular cavity lined with a multilayered fluid-filled shell of infinite length containing viscoelastic damping materials and embedded within a fluid-saturated permeable surrounding soil medium. The analytical results are illustrated with numerical examples where the effects of liner / coating structural arrangements, viscoelastic material properties, liner-soil interface bonding conditions, seismic excitation frequency, and angle of incidence on the induced dynamic stress concentrations are evaluated and discussed to obtain the representative values ​​of the parameters that characterize the system. It is characterized that the incorporating viscoelastic damping materials with a low shear modulus in the constrained layer configuration is an efficient means of enhancing the ove rall seismic isolation performance, especially for near-normal incident seismic shear waves where the amplitudes of induced dynamic stresses may be reduced by up to one-third of those without isolation in a relatively wide frequency range. Some additional cases are considered and good agreements with solutions available in the literature are obtained.