Os recentes colapsos catastróficos de barragens de mineração no Brasil intensificaram as exigências regulatórias, impondo o desafio de promover a descaracterização de um grande número de estruturas com o objetivo de mitigar riscos. Entretanto, ainda são incipientes as abordagens que integram o fluxo em regime parcialmente saturado às análises de estabilidade de estruturas descaracterizadas. Este estudo apresenta a implementação de um código computacional em MATLAB para avaliar a estabilidade de barragens descaracterizadas sem a remoção total do maciço, integrando simulações de fluxo em meio parcialmente saturado, fundamentadas na equação de Richards, e análises de estabilidade por meio da Análise Limite Numérica (ALN). A metodologia foi aplicada a uma barragem real já descaracterizada e validada com base em dados de instrumentação de campo. O estudo permitiu determinar a evolução temporal das poropressões e dos fatores de segurança (FS), que apresentaram variação anual da ordem de 2%. Os resultados demonstram o potencial do modelo para subsidiar o monitoramento de longo prazo e aprimorar os protocolos de avaliação de risco aplicáveis a essas estruturas geotécnicas.