Extending Mohr Coulomb state variables #742
Open
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Description
This PR includes new state variables for the Mohr-Coulomb material:
4 (binary) flags specifying the type of failure (if any). It can now be distinguished whether the material is currently (or has ever) failed in shear and/or tension
stored elastic energy. Note that it is denoted
E_el
, should it be rathereEl
?cumulated plastic work, i.e. dissipation, given as a positive number. (well, let's not argue here whether this really is a state variable from a theoretical perspective ;-) ). Which is currently
W_pl
(should it be ratherwPl
?)Additional context
I was about to modify mpm-doc as well to document this change, but have not yet fully decided about the most appropriate place to do so. Probably in source corresponding to https://mpm.cb-geo.com/#/theory/material/mohr-coulomb, in a new (final) section "Output state variables" ?
This code comes from @SchDvr work during his PhD thesis. There is quite a number of other updates that could be proposed to be applied following his work, but this should probably be further discussed elsewhere.