Hello everyone, I am trying to follow Wolfgang's suggestion here (https://groups.google.com/forum/#!topic/dealii/KrSp6k0FSx8) to implement a mixed formulation for elasticity. I'm experimenting with the Hu-Washizu 3-field variational form.
The most straightforward approach appears to be to declare dim basis functions for displacements, and dim*(dim+1)/2 basis functions for strains and stresses respectively. I am not exactly sure how the FEValuesExtractors::SymmetricTensor orders the elements it extracts into the symmetric tensor object. For instance, for a rank 2 tensor, does it assume that the first dim entries are the diagonal entries, and the next dim*(dim-1)/2 entries are the off diagonals? As an example for the stresses: [sigma_11, sigma_22, sigma_33, sigma_23, sigma_13, sigma_12] I believe this is the typical Voigt notation rules, but I'm not sure if dealii does it this way. I did try to look into the documentation but couldn't quite put my finger on it. Thanks for any inputs! Arjun -- 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]. For more options, visit https://groups.google.com/d/optout.
