Re: [deal.II] Output of Gauss point stress tensor

2020-09-29 Thread Muhammad Mashhood
Thanks a lot Prof. Wolfgang. Your guidance was very fruitful and right on point. Solved my problem! :) On Thursday, August 6, 2020 at 1:53:48 AM UTC+2 Wolfgang Bangerth wrote: > > > Thanks for the guidance. I tried > > replacing the "

Re: [deal.II] Output of Gauss point stress tensor

2020-08-05 Thread Wolfgang Bangerth
Thanks for the guidance. I tried replacing the " source/particles/particle_handler.cc

Re: [deal.II] Output of Gauss point stress tensor

2020-08-05 Thread Muhammad Mashhood
Dear Prof. Wolfgang, Thanks for the guidance. I tried replacing the " source/particles/particle_handler.cc " file with the one present in

Re: [deal.II] Output of Gauss point stress tensor

2020-07-31 Thread Wolfgang Bangerth
On 7/30/20 5:26 PM, Muhammad Mashhood wrote: _But it gives the following error on running:_ /An error occurred in line <298> of file in function     void dealii::Particles::Particle::set_properties(const dealii::ArrayView&) [with int dim = 3; int spacedim = 3] The violated condition was:

Re: [deal.II] Output of Gauss point stress tensor

2020-07-30 Thread Muhammad Mashhood
Dear Prof. Wolfganag, Thank you for the suggestion. I ran debugger with the ".prm" file and was successful to backtrace the error and resolved it. After this, using the step-19 particles tutorial refered by you, so far I am also successful in creating and

Re: [deal.II] Output of Gauss point stress tensor

2020-07-29 Thread Wolfgang Bangerth
On 7/28/20 11:05 AM, Muhammad Mashhood wrote: /.~/working_dir/step-42$ /step-42 p1_chinese.prm Segmentation fault (core dumped)/ I tried to comment out the line "  , particle_handler(triangulation, mapping, /*n_properties=*/dim) " and by doing so the code starts running normally. I have

Re: [deal.II] Output of Gauss point stress tensor

2020-07-28 Thread Muhammad Mashhood
Dear Prof. Wolfgang, Thank you very much for the suggestions. To have a quick solution, I am opting the first suggestion i.e. to make particles on quadrature points and assigning them the tensor quantities. And if I gain some reasonable results specific to

Re: [deal.II] Output of Gauss point stress tensor

2020-07-16 Thread Wolfgang Bangerth
On 7/16/20 9:15 AM, Muhammad Mashhood wrote: 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

Re: [deal.II] Output of Gauss point stress tensor

2020-07-16 Thread Muhammad Mashhood
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

Re: [deal.II] Output of Gauss point stress tensor

2020-07-15 Thread Wolfgang Bangerth
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

[deal.II] Output of Gauss point stress tensor

2020-07-15 Thread Muhammad Mashhood
Dear Deal.ii users, Hi! I am working on the deal.ii code. I have data of stress and strain tensors at quadrature points. Currently using the results_output function, I am transferring the data at nodal points of my triangulation using "