This work presentes the theory and the implementation techniques for na automatic h-adaptive procedure of finite elemento mesh refinements. Themetodology is applied to plane problems inlinear elasticity using triangular and quadrilateral shape elements with linear and quadratic lagrangian interpolation functions. Solution errors are a-posteriori evaluated using the Zienkiewicz- Zhu error estimator based on numerical solutions recovered from the finite element method solution – superconvergent patch recovery technique. For the mesh adaptivity of triangular and quadrilateral elements the remeshing concepts proposed by. Peraire and Zhu employed, respectively. The presence of solution singularities and their intensity is identified throughout the finite cloment mesh from the evaluation of the strain energy concentration, considering the numerical solution obtained. This scheme combined with adaptivity methodology allows for fully automatic procedure of mesh refinement/derefinement at the cycle of solution searching. Three examples are considered to demonstrate the procedure ability in handling elasticity problems: a cantiliver short beam under transverse distributed loading, a square plate with a middle square hole under inplace traction distributed loading and a double notch square plate under in plane distributed loading.