A short Matlab implementation of the P1 finite element method
is provided for the
numerical solution of viscoplastic and elastoplastic evolution problems
in 2D and 3D for von-Mises yield functions and Prandtl-Reuss
The material behaviour includes perfect plasticity as well as isotropic
and kinematic hardening with or without a
viscoplastic penalisation in a dual model, i.e. with displacements
and the stresses as the main variables. The numerical realisation,
however, eliminates the internal variables and becomes
displacement-oriented in the end. Any adaption from the given
examples to more complex applications can easily be performed
because of the shortness of the program and the given
documentation, which is available here.
In the numerical 2D and 3D examples an efficient error
estimator is realized to monitor the stress error.