Hallöchen!
Robert Cimrman writes:
[...]
The shape is, in general, (n_elements, n_quadrature_points, n_row, n_col). In "el_avg" mode, you get an average of all quadrature points values per element, so the shape in your case is (n_elements, 1, n_row, n_col). So each group, as you say, is the value of the stress tensor in the corresponding element. It has n_row = 3, n_col = 1, as it is symmetric 2x2 and we store only the items 11, 22, 12 (in this order).
Is this clear enough?
Absolutely, thank you. It enabled me to calculate the "equivalent tensile stress" (see http://en.wikipedia.org/wiki/Von_Mises_yield_criterion) for each element according to sqrt(a11**2 - a11 * a22 + a22**2 + 3*a12**2).
If this exceeds a certain threshold, the element is deleted. This way, I plan to let a crack extend through the material.
Tschö, Torsten.
-- Torsten Bronger Jabber ID: torsten...@jabber.rwth-aachen.de or http://bronger-jmp.appspot.com