An automatic adaptive refinement procedure for finite element analysis for two-dimensional stress analysis problems is presented. Through the combined use of the new mesh generator developed by the authors (to appear) for adaptive mesh generation and the Zienkiewicz-Zhu [Int. J. numer. Meth. Engng 31, 1331 1382(1992)] error estimator based on the superconvergent patch recovery technique, an adaptive refinement procedure can be formulated which can achieve the aimed accuracy very economically in one or two refinement steps. A simple method is also proposed to locate the existence and the position of singularities in the problem domain. Hence, little or no a priori knowledge about the location and strength of the singularities is required. The entire adaptive refinement procedure has been made fully automatic and no user intervention during successive cycles of mesh refinements is needed. The robustness and reliability of the refinement procedure have been tested by solving difficult practical problems involving complex domain geometry with many singularities. We found that in all the examples studied, regardless of the types of meshes employed, triangular and quadrilateral meshes, nearly optimal overall convergence rate is always achieved.