I'm trying to compute and output nodal stress and strain. Although there
are already some helpful discussions, for example this
<https://groups.google.com/forum/#!searchin/dealii/project$20stress$20and$20strain|sort:relevance/dealii/lvwwrT6zOFQ/C4r3JVAQAQAJ>
and this
<https://groups.google.com/forum/#!searchin/dealii/output|sort:relevance/dealii/uME5kwaJIEU/s5EN7hutAAAJ>,
none of them addresses my problem.
I computed the stress and strain on the quadrature points in each cell
and used FETools::compute_projection_from_quadrature_points_matrix
> [...]
This is one way of doing it, but an easier way -- at least for
visualization purposes -- is to use the DataPostprocessor framework to
directly output the gradient (or the stress, if you want) to the viz files.
Rather than describing this here, I thought I'd take the time to
actually write this up properly since your question is not the first one
to go in this direction. Take a look here:
https://github.com/dealii/dealii/pull/4991
This is not quite what you want: I only show how to do it for a scalar
solution where the gradient is a vector. In your case, the solution is a
vector and the gradient/stress a tensor. We don't currently support
tensor-valued output, so you'll have to output dim^2 components via the
DataPostprocessor class, but maybe you get the idea.
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.
For more options, visit https://groups.google.com/d/optout.
