SummaryIn this paper, new spherical Hankel shape functions are used to reformulate boundary element method for 2‐dimensional elastostatic and elastodynamic problems. To this end, the dual reciprocity boundary element method is reconsidered by using new spherical Hankel shape functions to approximate the state variables (displacements and tractions) of Navier's differential equation. Using enrichment of a class of radial basis functions (RBFs), called spherical Hankel RBFs hereafter, the interpolation functions of a Hankel boundary element framework has been derived. For this purpose, polynomial terms are added to the functional expansion that only uses spherical Hankel RBF in the approximation. In addition to polynomial function fields, the participation of spherical Bessel function fields has also increased robustness and efficiency in the interpolation. It is very interesting that there is no Runge phenomenon in equispaced Hankel macroelements, unlike equispaced classic Lagrange ones. Several numerical examples are provided to demonstrate the effectiveness, robustness and accuracy of the proposed Hankel shape functions and in comparison with the classic Lagrange ones, they show much more accurate and stable results.