My problem is not the theory. I am just not a deal.II expert and therefore it is still quite difficult for me to implement my approach. Going in to theoretical detail would take maybe much longer. Maybe it is better instead, if I reformulate my question like this:
There is a vector called cell_cf for instance. It contains all integrated material forces (summed up material forces from integration points) for each cell. Each cell has 4 vertices and 8 DoFs in 2D. Hence, it has to be a 8x1 vector. For instance: CONFIGURATIONAL FORCES WITHIN ELEMENT 0.010126 0.007502 0.011509 -0.012249 -0.008940 -0.009762 -0.012695 0.014508 Now imagine you have this vector which represents forces for each vertex, e.g. vertex 1 has 0.010126 in X-direction and 0.007502 in Y-direction, vertex 2 has 0.011509 in X-direction and -0.012249 in Y-direction and so on... How is one able now to take these vectors for each cell and assemble them correspondent to the DoF and summarize the values for coinciding nodes. Remember the values are computed within Gauss points and are just summarized over integration points to one single vector for each cell. It means there should be nothing stored on the nodes yet, but we have to take the information to the global structure and node/vertices level. I already tried to take the approach we implemented within our C++ code, but it is quite difficult to store data in deal.II since there is no processinfo or such a thing which tracks all data created and stored over the computation. I cannot save my configurational force vector and access it anywhere except I store it by using the setup_quadrature_point_history() and update_quadrature_point_history() functions from step-18. But this means I have to store my elastic energy then update and compute my material forces which are based on the energy and update again. This is not very nice to call a function several times for every little thing you want to save. Isn't there a better way to store data efficiently during runtime? Sorry for causing headaches ^^' Kind regards, S. A. Mohseni -- 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. For more options, visit https://groups.google.com/d/optout.
