A representação do comportamento de materiais elastoplásticos a partir de modelos com múltiplas superfícies de plastificação é uma alternativa para representar o comportamento de materiais como concreto, rochas e solos, que apresentam diferentes tipos de resposta não linear, a depender do estado de tensão atuante. No entanto, o emprego desses modelos requer a definição de muitos parâmetros que, por vezes, não possuem significado físico. Além disso, a implementação de modelos elastoplásticos com multiplas envoltórias traz complexidades adicionais. O emprego desse tipo de modelo requer um esquema robusto de integração das equações de evolução das variáveis plásticas. Nessa pespectiva, é apresentado um algoritmo de mapeamento de retorno baseado em um método de otimização sem restrições, o método de Newton-Raphson com busca unidimensional. Propõe-se uma expressão para o tensor constitutivo elastoplástico consistente para modelos com múltiplas superfícies de plastificação. Adicionalmente, apresenta-se uma metodologia de calibração dos parâmetros de tais modelos a partir da solução de um problema de otimização, solucionado via algoritmo genético. Para melhor compreender os parâmetros envolvidos nesse algoritmo, desenvolve-se um estudo paramétrico, solucionando uma série de problemas de otimização global. A robustez e eficácia dos algoritmos são avaliadas por meio de aplicações, dentre elas algumas disponíveis na literatura, num modelo constitutivo idealizado para concreto, rochas e solos: Cap Model. Por fim, calibrase tal modelo, considerando dados experimentais disponíveis na literatura. Assim, este trabalho tem como objetivo contribuir para viabilizar o emprego de modelos elastoplásticos complexos em problemas de engenharia.