Descrevem-se as considerações teóricas que serviram de base à elaboração dum programa de càlculo automático de análise da estabilidade global de aterros sobre solos moles reforçados com geossintéticos. O programa, formulado segundo pressupostos de equilíbrio limite, utiliza os resultados de aplicação dum modelo numérico bidimensional baseado no método dos elementos finitos. Apresentam-se e analisam-se os resultados de aplicação desse programa a uma obra de aterro sobre solos moles, comparando-se os resultados correspondentes à mesma obra reforçada e não reforçada.