Esse trabalho apresenta uma abordagem baseada em programação matemática para a solução de problemas de valor inicial de contorno constitutivo elastoplástico. Considerando a plasticidade associativa, as equações constitutivas locais, em sua forma discreta, são formuladas como problemas de otimização cônica. Especificamente, é demonstrado que métodos implícitos de retorno a superfície para os critérios mais conhecidos da literatura, como o de Rankine, von Mises, Tresca, Drucker-Prager e Mohr Coulomb, podem ser expressos como problemas de otimização cônica de segunda ordem e semidefinida. Além disso, um novo método numérico para a determinação do operador elastoplástico consistente, baseado na derivada paramétrica de primeira ordem das soluções ótimas, é proposto.