Dear all,
the fourth-order referential deviatoric tensor as returned by Physics::Elasticity::StandardTensors <https://www.dealii.org/current/doxygen/deal.II/classPhysics_1_1Elasticity_1_1StandardTensors.html>< dim >::Dev_P includes the fourth-order referential/spatial unit *symmetric* tensor Physics::Elasticity::StandardTensors <https://www.dealii.org/current/doxygen/deal.II/classPhysics_1_1Elasticity_1_1StandardTensors.html>< dim >::S = identity_tensor <https://www.dealii.org/current/doxygen/deal.II/symmetric__tensor_8h.html#ab3e890348aa219805e84f7d367e098c3> <dim>(). In the literature, for instance G. A. Holzapfel: "Nonlinear solid mechanics. A Continuum Approach for Engineering" (2007), however, the general fourth-order unit tensor I_{ijkl} = delta_{ik} delta_{jl} is used to compute Dev_P. I implemented a hyperelastic material model with a volumetric / isochoric split of the strain energy function and it only converges when using 'S' - as dealii does it. Using the general fourth-order unit tensor to define Dev_P, my solver does not converge at all. Is it neccessary to use 'S' due to the way dealii stores and accesses the elements of symmetric tensors? Best Simon -- 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/cc73ad68-f056-421d-b717-1bfcf1be9d4cn%40googlegroups.com.
