This work deals with the elastostatic identification of industrial manipulators. By reviewing the basics of the physical elastic properties of both links and joints in the framework of the lumped stiffness modeling techniques, the Gramian nature of the stiffness matrices has been found out adequate to do so. Then, a novel optimization method has been developed, which incorporates the Gramian matrix formulation along a non-linear optimization process, acting as an intrinsic constraint for the conservativeness of the elastostatic modeling. Numerical and experimental analyses evince the effectiveness of the proposed method, as the elastostatic models obtained by means of the proposed technique predict more than 93.7% of the compliance deviations of a real industrial robot. The proposed method is simple enough to be jointly applicable to the most recent elastostatic model reduction techniques.