On 12/13/22 00:00, Deepika Kushwah wrote:
1. I have 2 materials in a single domain for which I already defined
material_id in triangulation. Now how to call these material_id while
calculating strain values? (as material_ids are only related with lambda and mu).
In my previous email, I already pointed at the link to the function with which
you can obtain the cell iterator from the DataPostprocessor class you implement.
You then only need to do something along the lines of this:
DoFHandler<dim>::cell_iterator cell = ...;
double lambda;
if (cell->material_id == ...)
lambda = ...; // first material
else
lambda = ...; // second material
2. I wanted to calculate stress for the same problem which I defined above. I
have the following relation but how to write this constitutive matrix (C)
while using datapostprocessor class?
In the formula you show, you are using Voigt notation. If you want to use
that, you need to copy the solution gradients into a Voigt vector, multiply it
with the matrix you show (probably using the FullMatrix class), and the result
will be a Voigt vector of the stress.
We tend not to use Voigt notation in the tutorials and instead represent
strain and stress via SymmetricTensor<2,dim> objects and the C operator as a
SymmetricTensor<4,dim>. You can find numerous examples of using these in the
tutorials.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/62c185b4-96c9-cc41-83c8-b62c276883ad%40colostate.edu.
