In my question I meant whether is it possible to evaluated plastic strain
component from currently implemented plasticity algorithm as a further
development of this code?
I am pretty sure that it is, by just computing the difference between the
elastic stress (C eps(u)) and the actual stress computed. In fact, for the
current situation, the actual stress computed equal to the elastic stress
where it is less than the yield stress (and so the plastic strain is zero),
and it is simply a fraction of the elastic stress where it exceeds the yield
stress. Once you have the plastic stress, you can compute the plastic strain
by multiplying it by C^{-1}.
Then as a step further I would be trying to store this plastic strain in cells
or Gauss points along with the modified yield strength (due to isotropic
hardening) so that history of loading is stored too in the domain.
I imagine that that, too, can be done. I'm not an expert in plasticity, but I
see no fundamental reasons why what you want to do should not be possible.
There are also classes CellDataStorage and TransferableQuadraturePointData and
parallel::distributed::ContinuousQuadratureDataTransfer that can help you with
storing information at quadrature points.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/833e007c-7f5a-31d4-37ed-eafd021ce0ff%40colostate.edu.
