In this article, we present a free-vertex tetrahedral finite-element representation of irregularly shaped small bodies, which provides an alternative solution for estimating asteroid density distribution. We derived the transformations between gravitational potentials expressed by the free-vertex tetrahedral finite elements and the spherical harmonic functions. Inversely, the density of each free-vertex tetrahedral finite element can be estimated via the least-squares method, assuming a spherical harmonic gravitational function is present. The proposed solution is illustrated by modeling gravitational potential and estimating the density distribution of the simulated asteroid 216 Kleopatra.