The thesis studies nonlinear solution algorithms for problems in
infinitesimal elastoplasticity and their numerical realization within
a parallel computing framework. New algorithms like Active Set and
Augmented Lagrangian methods are proposed and analyzed within a
semismooth Newton setting. The analysis is often carried out in
function space which results in stable algorithms. Large scale
computer experiments demonstrate the efficiency of the new algorithms.