RESUMO: A ampla aplicação de vibrações na engenharia tem motivado pesquisas recentes nessa área. No entanto, para a grande maioria dos problemas, soluções exatas não estão disponíveis, tornando imprescindível o uso de métodos numéricos, sendo o mais usual, o Método dos Elementos Finitos (MEF). Diante disto, o presente trabalho possui o objetivo de investigar como os parâmetros da discretização do MEF influenciam na aproximação das frequências naturais em um problema de vibração livre não amortecida. Os parâmetros considerados incluem: tipo, tamanho e ordem polinomial dos elementos finitos. Além disso, pretende-se explorar a relação entre a exatidão da solução aproximada com as diferentes configurações das deflexões dinâmicas dos modos de vibração. O problema de engenharia adotado é de uma viga bidimensional engastada na extremidade esquerda, elástica e linear e a simulação é realizada pelo software Ansys®, considerando um estado plano de tensão. Os resultados obtidos indicam uma melhor aproximação para os modos predominantemente axiais. Foi observado também que, em comparação aos elementos de ordem linear, os elementos de ordem quadrática apresentam maior exatidão na aproximação. Por fim, constatou-se que os elementos quadrangulares são mais eficientes na aproximação do que os elementos triangulares.