Este trabalho tem por objetivo o desenvolvimento de uma ferramenta numérica para avaliação do fluxo saturado-não saturado em encostas de grandes dimensões, com aplicação na análise de estabilidade dessas áreas. Emprega-se o método dos elementos finitos na solução da equação de Richards, considerando a carga de pressão como variável primária e utilizando formulação adequada para minimização dos problemas de conservação de massa, freqüentemente, associados a esse fato. O modelo constitutivo utilizado para a curva característica e função de condutividade hidráulica é o proposto por Van Genuchten (1980). Para solução da não-linearidade, emprega-se um método quasi- Newton (BFGS). Com o objetivo de minimizar os requisitos de memória computacional, utiliza-se esquema de armazenamento de matriz esparsa, associado ao método de gradiente biconjugado, na solução do sistema de equações. Paralelamente, é apresentado algoritmo de geração de malha tridimensional de elementos finitos, a partir de uma malha superficial de triângulos, representativa da topografia. Análises numéricas são executadas com a finalidade de validação do código gerado, comparando-se os resultados obtidos com aqueles gerados por outros programas já consagrados na literatura técnica. É proposta metodologia para geração de mapas de susceptibilidade a escorregamentos translacionais rasos, empregando-se o método do talude infinito, associado à estrutura da malha de elementos finitos e aos resultados do problema de fluxo, incorporando- se, assim, o efeito do estado de não saturação na resistência do material.