[deal.II] Issues with single element boundary conditions

Hello, it is the first time i am using the GridGenerator and i am trying to create just a single element (3D) and add dirichlet BC on one side. I did the following: repetitions[0] = 1; repetitions[1] = 1; // Only allow one element through the thickness // (modelling a plane str

[deal.II] Re: Changing time step

Alright guys i was able to solve this. Thanks BBah schrieb am Montag, 18. Juli 2022 um 11:00:56 UTC+2: > Hello everybody, > > i am trying to change my timestep after few simulations steps but i am > getting a kink into my solution once i reduce my timestep. > I did this act

[deal.II] Changing time step

Hello everybody, i am trying to change my timestep after few simulations steps but i am getting a kink into my solution once i reduce my timestep. I did this actually simple by adjusting my timestep once a certain time is reached when i apply the dirichlet boundary conditions: if(ti

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

/20/22 01:20, BBah wrote: > > Hello Mr.Bangerth. Thank you very much for your input and your help. I > did > > implement the same loop actually in my code but i am getting a > segmentation > > fault but i guess this is a minor issue that i cannot find. Maybe > because

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

quot;);* Since i am using here the CellDataStorage i was not able to use the same loop as in step 18. I hope this helps somebody too. BBah schrieb am Freitag, 20. Mai 2022 um 09:20:36 UTC+2: > Hello Mr.Bangerth. Thank you very much for your input and your help. I did > implement the same lo

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

MappingQEulerian q_mapping(degree, dof_handler_ref, soln);* *data_out.build_patches(q_mapping, degree);* *std::ostringstream filename;* *filename << "solution-" << time.get_timestep() << ".vtk";* *std::ofstream output(filename.str().c_str());* *data_o

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

soln(i) = solution_n(i);* *MappingQEulerian q_mapping(degree, dof_handler_ref, soln);* *data_out.build_patches(q_mapping, degree);* *std::ostringstream filename;* *filename << "solution-" << time.get_timestep() << ".vtk";* *std::ofstream

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

Like u mentioned i should be able to use the same loop without the workstream::run function but i do not know why this is not working. I guess it is a really small issue but i cannot find it. BBah schrieb am Freitag, 20. Mai 2022 um 09:03:00 UTC+2: > Hello Mr.Bangerth. I already implemen

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

:cout< *>(cell-> user_pointer())[q].old_stress; } } } I am able to compile this but running leads to a segmentation fault. Wolfgang Bangerth schrieb am Freitag, 20. Mai 2022 um 02:02:33 UTC+2: > On 5/19/22 00:47, BBah wrote: > > > > i am just tr

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

One more hint maybe. I have just modified the code " Quasi-Static Finite-Strain Compressible Elasticity " for my purposes. If anybody know how to put out stresses from this code i would be able to implement it into my code too. BBah schrieb am Donnerstag, 19. Mai 2022 um 12:17:07 U

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

I just can't find any sample or example how to put out my stress tensor. It should be pretty easy because my stress is stored for each cell in my PointHistory class but how to access these stresses and put them out ? I hope somebody can help me BBah schrieb am Donnerstag, 19. Mai 2022 um

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

: > > BBah, > I have to admit that I don't understand what precisely your question is -- > your first email does not contain enough information what you have already > tried, or what exactly you want to do. > > But regarding what you have below: segmentation faults are th

[deal.II] Re: Outputting stress tensor at each quadrature point

atics::F(Grad_u_n); const NumberType det_F = determinant(F); const SymmetricTensor<2,dim,NumberType> C = Physics::Elasticity:: Kinematics::C(F); old_stress = material->get_s_piola(det_F,C); } } } ... BBah schrieb am Mittwoch, 18. Mai 2022 um 14:52:56 UTC+2

[deal.II] Outputting stress tensor at each quadrature point

Hello everybody, i am trying to output my stresses by using a similar approach like in Step 18: *Vector norm_of_stress(triangulation