Summary: We investigate a finite element formulation of the exponentiated Hencky-logarithmic model whose strain energy function is given by \[ W_{\mathrm{eH}}(F) = \frac{\mu}{k}\, e^{k ||\mathrm{dev}_n \log U| |^2} + \frac{\kappa}{2 \hat{k}}\, e^{\hat k [\mathrm{tr} (\log U)]^2}, \] where \(\mu >0\) is the (infinitesimal) shear modulus, \(\kappa >0\) is the (infinitesimal) bulk modulus, \(k\) and \(\hat{k}\) are additional dimensionless material parameters, \(U=\sqrt{{F^T}{F}}\) is the right stretch tensor corresponding to the deformation gradient \(F\), \(\log \) denotes the principal matrix logarithm on the set of positive definite symmetric matrices, \(\mathrm{dev}_n X = X-\frac{\mathrm{tr}X}{n}1\) and \(|| X | | = \sqrt{\mathrm{ tr}X^TX}\) are the deviatoric part and the Frobenius matrix norm of an \(n\times n\)-matrix \(X\), respectively, and \(\mathrm{tr}\) denotes the trace operator. To do so, the equivalent different forms of the constitutive equation are recast in terms of the principal logarithmic stretches by use of the spectral decomposition together with the undergoing properties. We show the capability of our approach with a number of relevant examples, including the challenging “eversion of elastic tubes” problem.