Problemas de engenharia geotécnica onde a resistência ao cisalhamento é fator determinante para o projeto são geralmente avaliados com base no critério de ruptura de Mohr-Coulomb; no entanto, modelos constitutivos avançados são necessários para abordar grande parte de problemas geotécnicos onde estimativas de deformação são necessárias. Lade e Kim propuseram um modelo constitutivo elastoplástico avançado para investigar o comportamento tensão vs. deformação de materiais com atrito, sob condição drenada ou não drenada. Nesta dissertação, este modelo foi implementado no programa computacional de diferenças finitas FLAC3D com o objetivo de estudar o comportamento elasto-plástico de aterros. As equações governantes do modelo constitutivo de Lade-Kim são apresentadas e um algoritmo explícito de integração de tensões foi implementado na presente pesquisa. Os resultados numéricos obtidos com o modelo de Lade e Kim foram comparados com aqueles publicados na literatura, envolvendo ensaios triaxiais de laboratório realizados em areia e análises numéricas 2D do comportamento de barragem. Finalmente, após validação da implementação do modelo, análises do comportamento 2D e 3D de aterros foram executadas, verificando a influência destes tipos de análise nos valores calculados, bem como resultados numéricos foram também comparados com aqueles determinados pelo tradicional modelo elasto-plástico de Mohr-Coulomb.