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

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

Thanks for the guidance. I tried replacing the " source/particles/particle_handler.cc

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

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

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

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

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

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

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

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

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 "