Two-yield elastoplasticity solver
Provided Matlab software demonstrates two-yield elastoplasticity quasi-static calculation. The two-yield models is the simplest version of a multi-yield model, which generalizes the classical model of linear kinematic hardening in elastoplasticity. Theory and numerics can be found in Jan Valdman's PhD thesis:
Mathematical and Numerical Analysis of Elastoplastic Material with Multi-Surface Stress-Strain Relation
submited in 2001 (under supervision of Martin Brokate and Carsten Carstensen) at CAU Kiel and published at Dissertation.de Verlag in Berlin in 2002 under ISBN 3-89825-501-8. An improved pdf version of the thesis can be downloaded here.
Two journal papers (-> file-link, -> file-link) directly based on the thesis followed and other two papers extended original results (-> file-link, -> file-link).
Two animations from the solver (colors meaning: black - elastic zone, brown - first plastic zone, light yellow - second plastic zone):
Animation 1: Time evolution of 2D two-yield beam Animation 2: ZZ space adaptivity within one time-step
Software development:
28.08.2002 Original Matlab Solver from the PhD thesis provided.
08.11.2002 This html page added.
14.02.2003 rewriting Solver in C++ (as a part of NGSolve), some pictures and animations here.
14.04.2010 Matlab Solver extended and upgraded here.