yade-users team mailing list archive

[Jan Stránský <question690521@xxxxxxxxxxxxxxxxxxxxx>]

Question #690521 on Yade changed:
https://answers.launchpad.net/yade/+question/690521

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hello,

just follow standard plasticity concepts:
- you know current total displacement dTot and plastic displacement dPl
- compute elastic displacement dEl = dTot - dPl
- compute trial stress sTrial = k * dEl
- do stress return to admissible value according to F (possibly also involving change of damage) => sActual
- compute "plastic stress" sPlastic = sTrial - sActual
- plastic strain increment is simply sPlastic / k

cheers
Jan

PS: there is a bug in the current implementation [1]...
[1] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/dem/ConcretePM.cpp#L436

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.