Dear Seyed, On top of what Wolfgang has already said, I have a few more direct comments about what you posted a few days ago.
> The computation of configurational forces is accomplished by the following > formula: > > > > <https://lh3.googleusercontent.com/-B38RUw135-c/WIjDUOvqsfI/AAAAAAAAAEI/4rD24zbgMeoI9O-2MFpXt77UTbGAHjmLQCLcB/s1600/config_forces.png> > > If one reinterprets the Eshelby (energy momentum) stress tensor as any other "regular" stress tensor, then the computation of $\mathbf{g}$ is exactly the same as the assembly of the RHS vector in nonlinear elasticity (without body forces and tractions). You can compare your implementation to that of step-42 and step-44. > > > > How can I use these new features? Is there a step example or tutorial? Is > this TransferableQuadraturePointData the same like in step-18? > CellStorageData is now used in step-44 (development version); for TransferableQuadraturePointData and ContinuousQuadratureDataTransfer you'll have to look at the unit tests; see this post <https://groups.google.com/forum/#!topic/dealii/GRZMUTLIm2I> for details on how to do that. > I just did this: > > std::vector<types::global_dof_index> local_dof_indices(dofs_per_cell); > > typename DoFHandler<dim>::active_cell_iterator cell = dof_handler. > begin_active(), endc = dof_handler.end(); > for (; cell != endc; ++cell) > { > if ( cell->is_locally_owned() ) > { > fe_values.reinit(cell); > > ... configurational forces computation ... > > } > > configurational_forces = cell_cf; > > cell->get_dof_indices(local_dof_indices); > > constraints.distribute_local_to_global(cell_cf, local_dof_indices, > configurational_forces); > } > > > There are a few potential points that you can check here, and a few to be concerned about: - Is the Eshelby stress tensor non-trivial on all (locally owned) cells? - Is cell_cf zero on cells that aren't locally owned? Why do you still assemble on non-locally-owned cells? - Is cell_cf non-trivial on cells that are locally owned? - Why do you set "configurational_forces = cell_cf;" ? Isn't the one a global vector and the other one with local contributions? This is suggested by the call to distribute_local_to_global. Does this even run in debug mode? So what exactly do you mean by repopulating each cell? > I queried whether you called cell->get_dof_indices(local_dof_indices); for each cell before distribution. I agree with Wolfgang that its probably most helpful to you if you become more familiar with the "deal.II way" by studying the tutorials (even some not directly related to your topic of interest) more thoroughly. Best regards, Jean-Paul -- 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.
