Várias formulações numéricas da análise limite têm sido propostas para a solução de problemas geotécnicos de estabilidade. Nestas formulações, o problema de estabilidade é colocado como um problema da Programação Matemálica. Este problema pode ser resolvido diretamente por métodos de programação não linear (PNL) ou aproximado por um problema da programação Iinear (PL). Este trabalho apresenta uma formulação mista de análise limite para solucionar problemas de estabilidade da geotecnia em meios contínuos, meios fraturados ou solos reforçados, sob condição de deformação plana, pelo método dos elementos finitos. Para a descrição do comportamento do maciço rochoso, o critério de Hoek e Brown pode ser usado como função de escoamento. Para o solo é empregado o critério de escoamento de Mohr-Coulomb na forma linearizada e a condição de escoamento de interfaces é descrita pelo critério de Coulomb. O solo, as interfaces e o reforço são idealizados com comportamento rígido-plásticos perfeitos, na hipótese de plasticidade associada. Os exemplos de problemas de estabilidade em rocha foram resolvidos pela programação não linear. Os problemas de estabilidade em solos e solos reforçados foram resolvidos pela técnica de PL. Os resultados obtidos nos exemplos estudados pelas implementações numéricas indicam que a metodologia a presentada t em considerável potencial para ser aplicada a problemas de engenharia geotécnica.