Dear Prof. Wolfgang, Thank you for the comprehensive and informative response. As my goal is to access the original Gauss point values and I do not want to involve any interpolation or approximation scheme so I would like to opt for the first option among the two you mentioned:
*"If you want to output the stress/strain information, you have two options: * ** You actually do show them as values defined only at individual points, rather than as fields. In this scheme, the stress/strain really only exists at the quadrature points, but not anywhere in between. You can't create surface plots, you can't create isocontour plots, etc. * You create a field that somehow approximates these point values. There are different ways of doing that, and the approach you are currently using is one way for this."* So far it would be enough if I have the Gauss point values only at the points rather than having complete field. I like to access these point values of stress and strain tensors in .vtk file or .pvtu file through ParaView (same as I am accessing the results currently). I would be grateful if you or any other user has experience in this regard or know about relevant deal.ii tutorial dealing with same feature. Thank you! On Thursday, July 16, 2020 at 4:59:31 AM UTC+2 Wolfgang Bangerth wrote: > On 7/15/20 11:05 AM, Muhammad Mashhood wrote: > > > > My question is that along with this, can I also export the stress and > strain > > tensors data of quadrature points i.e. > > "local_history_strain_values_at_qpoints[i][j]" directly to the output > file > > (.pvtu or .vtu etc.) i.e. without mapping or averaging etc. on the node. > In > > this way I want to directly visualize the quadrature point values of > stress > > and strain in Paraview. > > The question is more a philosophical one than one of how you can achieve > this. > Typically, when we output information in finite element contexts, we > output > them as "fields", i.e., functions of space. This allows us to show them as > surfaces, with color gradients, etc. The strategy you have found of taking > values defined in (quadrature) points and converting them to fields is a > way > to make that happen. > > If you want to output the stress/strain information, you have two options: > > * You actually do show them as values defined only at individual points, > rather than as fields. In this scheme, the stress/strain really only exists > at the quadrature points, but not anywhere in between. You can't create > surface plots, you can't create isocontour plots, etc. > * You create a field that somehow approximates these point values. There > are > different ways of doing that, and the approach you are currently using is > one way for this. > > So the question is mostly: What's your goal with this? Do you want to show > these quantities as fields, or as data to be visualized at individual > points? > > Best > W. > > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/bc35a6d4-f3d2-4c57-a193-b26afb6aaa5dn%40googlegroups.com.