Neste trabalho analisa-se o problema de contato entre dois corpos bidimensionais ou tridimensionais com comportamento elasto-plástico submetidos a grandes deformações através de duas metodologias de solução implementadas. O procedimento de solução numérica inclui o tratamento conjunto das não-linearidades geométricas e físicas como também da não-linearidade das condições de contato baseando- se numa estratégia incremental-iterativa. Na modelagem dos corpos utiliza-se o método dos elementos finitos em conjunto com uma formulação Lagrangeana Atualizada e relações constitutivas para grandes deformações elasto- plásticas. Os corpos são considerados ambos deformáveis, supondo a superfície de contato sem atrito. Os métodos implementados foram: o método da Penalidade onde as restrições de contato são satisfeitas de forma aproximada através do emprego do parâmetro de penalidade e o método da Programação Linear Complementar que formula um problema de programação matemática para cada configuração de equilíbrio incremental onde as condições de Kuhn-Tucker são resolvidas pelo esquema de pivoteamento de Lemke [43]. Foi realizada uma análise comparativa entre os métodos e também um estudo paramétrico verificando a influência de parâmetros relevantes para o problema. Foi implementado um algoritmo para consideração da geometria de contato onde são supostas condições genéricas de contato: nó-superfície e onde as relações cinemáticas são dadas em termos de uma função diferenciável da distância.Foram utilizados elementos lineares formulados em deslocamentos, hexaédricos de 8 nós e elementos com uma formulação híbrida - Enhanced Assumed Strain - EAS, que possuem, além de oito nós, três parâmetros internos de deformação, com um campo de deformações trilinear completo, minorando os efeitos de -locking- presentes nos elementos lineares formulados em deslocamentos. O emprego desses dois tipos de elementos teve como finalidade comparar as suas performances em problemas de contato.