Dear everyone,
I am now interested in solving the elasticity problem. I have followed
step 8, and adapted it a little using step 20 in order to use
BlockSparseMatrix and BlockVector which I found very elegant.
Now, after solving the linear system, I would like to compute the strain
and stress tensors on each integration point and output them using
DataPostprocessor, exactly as I did it for the heat flux following steps
21 & 29. DataPostprocessor has compute_derived_quantities_scalar and
compute_derived_quantities_vector but no
compute_derived_quantities_tensor, so I tried to output the tensors in
"engineering notations" (e.g. for the strain tensor in 2d, [\epsilon_11
\epsilon_22 \epsilon_12]). However, when running the code I have the
following message:
Furthermore, vectors must always consist of exactly <dim> components. However, the
vector component at position 2 with name <strain> does not satisfy these conditions.
It seems that I can not output vector of size dim*(dim+1)/2 using
DataPostprocessor. Is there another way do it, or do I have to split the
strain and stress tensors into many components?
Thanks in advance,
Martin.
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
