On 10/25/19 11:06 AM, 'Andreas Hegendörfer' via deal.II User Group wrote:
>
> for (unsigned int i=0; i<cell->get_fe().dofs_per_cell; ++i)
> {
> strain=fe_values[displacement].symmetric_gradient (i, q_point)
> *solution [local_dof_indices[i]]+strain;
You can do the first part simpler by using
fe_values[displacement].get_function_symmetric_gradients(...)
> }
> strain_solution[local_dof_indices[q_point]]=strain;
So here you put 'strain' into a solution vector. This solution vector
must be associated with a finite element. And if it is associated...
> In the programm I do have only one hp::DoFHandler<dim> dof_handler
> overall. I am using the two elements: substructure_fe (FE_Q<dim>(2),
> dim, FE_Nothing<dim>(), 1 ) and piezo_fe (FE_Q<dim>(2), dim,
> FE_Q<dim>(2), 1).
...with this element, then you will necessarily have a continuous space
into which you put the strains. If you don't like that, you'll have to
create a different (discontinuous) FE, DoFHandler, and associated vector.
> I save the result of the last calculation at one support point in
> strain_solution[local_dof_indices[q_point]]. So the ohter calculation
> results are lost but it should be no effort to save them, too.
Correct.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/1e523ccf-c490-a493-0599-d5b90e5c476d%40colostate.edu.
