Esta dissertação tem por objetivo a implementação de uma ferramenta numérica para avaliação do fluxo transiente 3D saturado-não saturado em barragens de terra e enrocamento, baseado no método dos elementos finitos e no programa GEOFLUX implementado por Machado Jr. (2000) para análise de problemas 2D. Nesta nova versão, foram incluídos elementos triangulares de 3 nós para análises 2D e elementos tetraédricos de 4 nós para análises 3D. Implementam-se também subrotinas que oferecem a possibilidade de variação das condições de contorno com o tempo. A equação de Richards é solucionada considerando a formulação mista e o método iterativo de Picard Modificado para solução do sistema de equações não- lineares. Para a solução do sistema de equações utiliza-se um armazenamento especial para matrizes esparsas associado com o método do gradiente bi-conjugado, tornando o processo muito rápido, mesmo em sistemas de grande porte. Utilizam- se dois modelos para representar as curvas características: o modelo exponencial proposto por Srivastava e Yeh (1991) e o modelo proposto por van Genuchten (1980). O programa computacional desenvolvido (GEOFLUX3D) foi aplicado na análise de fluxo na barragem de enrocamento de Gouhou, China, e na barragem de terra Macusani, Peru. Os resultados numéricos indicam a necessidade de análises numéricas 3D em barragens situadas em vales estreitos, onde os efeitos de geometria nas condições de fluxo são significativos e não podem ser ignorados.