In this work we present a general error estimator for the finite element solution of solid mechanics problems based on the Variational Multiscale method. The main idea is to consider a rich model for the subgrid scales as an error estimator. The subscales are considered to belong to a space orthogonal to the finite element space (Orthogonal Subgrid Scales) and we take into account their contribution both in the element interiors and on the element boundaries (Subscales on the Element Boundaries). A simple analysis shows that the upper bound for the obtained error estimator is sharper than in other error estimators based on the Variational Multiscale Method. Numerical examples show that the proposed error estimator is an accurate approximation for the energy norm error and can be used both in simple linear constitutive models and in more complex non-linear cases.