A análise de concentração de tensões é um tema de extrema importância em diversas áreas da engenharia, porém de díficil solução analítica. Com o desenvolvimento dos microcomputadores, a análise numérica destes problemas, por meio do Método dos Elementos Finitos, tornou-se comumente difundida. Neste trabalho, investiga-se a evolução dos campos de tensão em uma malha poliédrica por meio de técnicas MEF, visando obter numericamente o valor das concentrações de tensão induzidas pela geometria do domíno escolhido. Primeiro, são definidos os domínios de interesse de forma implícita. Elementos poliedrais são gerados, na sequência, pela triangulação de Delaunay e diagrama de Voronoi, sendo estes analisados de pela função de distância sinalizada, para verificação do pertencimento ao domínio. Após a etapa inicial de geração de elementos, é feita simulação por elementos finitos, e então obtem-se as tensões médias em cada elemento, pela média das tensões de cada nó do mesmo, a partir do critério de Von Mises. Os elemento com maiores tensão, são automaticamente selecionado entrando numa rotina de refino da malha, até que se obtenha uma melhor percepção dos campos de tensão. Finalmente, apresentam-se diversos exemplos que demonstram os recursos e o potencial da abordagem proposta e validam os os resultados obtidos.