A new procedure is suggested for quantization of the system potential energy variation for elastic isotopic bodies as applied to the construction of resolving relations of the finite elements methods. Divergence and vortex of the displacement vector are taken as basic characteristics approximated by the deformation characteristics element. The procedure for constructing the rigidity matrices is illustrated for the convex polygonal elements.