Hi Jean-Paul and Wolfgang,

Exactly, I would like to use what my DataPostprocessor computes in other 
parts of the programme. 

Here my problem:

In my postprocessor, I compute stresses and principal stresses for all 
quadrature points (to give you the context, I am using quad elements with 4 
quadrature points. The stress is constant for each cell/element as I am 
using quadratic functions for the displacement. I am solving a plane stress 
problem, isotropic elastic material). Principal stresses can be also 
computed. 

So, for each element, I get the same principal stress at the 4 quadrature 
points (they are constant for each element). If I plot them, I get the 
following in Paraview (see screenshot.PNG)


Shared vertices with different principal stress components of different 
elements. To solve this, I would like to get the value of the principal 
stress of each element and plot it in the centre of the element, so that I 
can draw the arrow with the principal stress direction. For this I need the 
coordinates of the centre of the element and the value of stress. I do not 
know how to export these variables in other parts of the programme, so I 
can build the paraview file to plot that. 

I still struggle to understand some of the c++ structures and how they can 
be accessed in the virtual functions of the dealii library. 

I hope my problem is a little clearer now, thank you!!!

Raúl

El martes, 13 de diciembre de 2022 a las 1:05:13 UTC+1, Wolfgang Bangerth 
escribió:

> On 12/12/22 08:01, Raúl Aparicio Yuste wrote:
> > These variables are saved on the data_out object, so I use 
> > “data_out.write_vtk(output);” and I get them on the vtk file. However, I 
> do 
> > not know how to access to these new variables (computed_quantities) from 
> > another file different from the one where I declare the functions.
> > 
> > The function evaluate_vector_field is a virtual function, and it is 
> declared 
> > as const override. I do not know how to export the variables 
> > computed_quantities or modify them, because it is only allowed to read 
> them 
> > since the function is declared as constant. Is there any way to access 
> the 
> > structure postprocessor? Something like:
> > 
> > postprocessor.getVariables(xxxxxxxx)
> > 
> > Or an example where I can see how to implement it. Thank you very much 
> in advance
> > 
>
> Raul -- in addition to J-P's answer, start by explaining *what* it is you 
> want 
> to do, not *how* you want to do it. Do I understand right that you want to 
> use 
> what the DataPostprocessor computes in other parts of the program as well? 
> If 
> so, why? Are you trying to solve a nonlinear system and you want to reuse 
> the 
> computed stress for other computations?
>
> Best
> W.
>
> -- 
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: [email protected]
> 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 [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/5bd4d9ef-ec02-4469-8fdd-bf5c586540edn%40googlegroups.com.

Reply via email to