Apresenta-se neste trabalho a implementação de um método adaptativo de refinamento automático de malhas de Elementos de Finitos. O método aplica-se a problemas planos de elasticidade linear. Elementos triangulares e quadrilaterais (lineares e quadráticos) são utilizados. Estimativas de erro a-posteriori são obtidas através do estimador de erro proposto por Zienkiewicz-Zhu, empregando-se uma técnica de recuperação da solução numérica fornecida pelo MEF (superconvergent patch recovery technique). A metodologia da adaptação da malha emprega o conceito de remalhamento do domínio na forma proposta por Peraire e Zhu para malhas de elementos triangulares e quadrilaterais, respectivamente. Implementa-se um esquema de localização de singularidades presentes no domínio baseado no conceito de concentração da energia de deformação. O uso deste esquema combinado com a metodologia de adaptação permite um processo automático sem intervenção do usuário durante os ciclos de refinamento da malha. Para demonstrar-se a performance do método são analisados três problemas: uma viga curta em balanço sob carregamento distribuído, uma placa quadrada com furo quadrado e uma placa quadrada com dupla trinca.