No presente trabalho é proposto um modelo constitutivo e um método numérico implícito para a análise do comportamento dúctil de rochas salinas sob deformações finitas. Para modelar a fase primária e secundária da fluência emprega-se o modelo de Multimecanismos de Deformação proposto por Munson, Fossum e Senseny, o qual incorpora os fenômenos de endurecimento e recuperação das rochas salinas quando sujeitas a um processo de deformação. No modelo de viscoplasticidade supõe-se uma decomposição multiplicativa do gradiente de deformação em uma parte elástica e uma parte viscoplástica, e considera-se um modelo não associativo de Druker-Prager. A formulação do modelo segue com uma descrição Lagrangeana Total e considera as equações constitutivas escritas em termos da medida de deformação logarítmica de Hencky e da tensão rotacionada de Kirchoff. Além disso, propôs-se um método numérico implícito com o objetivo de determinar soluções numéricas aproximadas utilizando o método de elementos finitos de Galerkin. Adicionalmente, um operador tangente consistente contínuo, associada ao método implícito, é derivado por uma linearização adequada da forma fraca da equação de equilíbrio. Alguns resultado numéricos são apresentados, sob condições axissimétricas e de estado plano de deformações, para verificar a eficiência do modelo constitutivo e o desempenho do método numérico proposto.