In this paper, the dynamic analysis of shear deformable shells was conducted using the hybrid displacement
boundary element method (HD-BEM). The approach involved utilizing the fundamental solution as a trial
function to derive the 8-node element stiffness matrix, while approximating the boundary variables with
quadratic shape functions. To ensure efficient computations and avoid singular integrals, the source points were
strategically positioned outside the domain. For accurate calculations of the proposed 8-node element mass
matrix, only four Gauss points were relied upon. It is important to note that both the element mass and stiffness
matrices demonstrate symmetry. Throughout the study, free, forced, and harmonic analyses were considered.
Encouragingly, the results obtained from the study exhibited exceptional accuracy levels, even when employing a
coarse discretization approach.