Hi all,
I have a hyperelastic code similar to Step-44. I am able to compute the
Cauchy stress \sigma at the quadrature points. Based on this, how do I
further compute the divergence (w.r.t. current configuration) of Cauchy
stress as a vector?
The following is my thought:
1. Define the stress field associated to the nodes using either DG or FE_Q.
Compute the stress at quadrature points and then extropolate to the nodes.
Depending on how the stress field is represented, cell averaging may or may
not be needed, which is essentially the same as the task of "outputing
stress" that has been discussed in the mailing list.
2. Given a stress field, FEValues::get_function_gradients can be used to
compute the stress gradient from which the divergence can be obtained.
My questions are,
a. is the gradient computed from get_function_gradients w.r.t. reference
configuration or current configuration? In case of large deformation, do I
need to multiply the computed gradient with a F^{-1}?
b. currently I define every component in the nodal stress as a scalar
field, which is cumbersome. I've only seen scalar field and vector field in
the tutorials, but is it possible to have a tensor-valued field?
Thank you
Riku
--
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.
