Pieter,

I would definitely second Wolfgang's suggestion here. I struggled for a long time to implement a finite-deformation elastoplasticity model by moving vertices, but this approach introduces mapping inconsistencies when using higher-order elements. The MappingQEulerian class provides an easy fix for this problem.

Josh




Pieter,

I have a problem similar to the step-18 tutorial program, but would like
to use higher order (biquadratic) elements in combination with a
MappingQ of degree 2.

I am trying to adapt the move_mesh method from example program 18 to not
only move the triangulation vertices, but to move all support points.
However, changing the support points does not seem to move the mesh.

Does it make sense to try to move these support points? Is there an easy
way to do so?

I think this is the wrong way. Support points are mapped from reference points by the Mapping classes, which is independent of the mapping you do with
vertex locations.

Have you looked at the MappingEulerianQ class?

W.

-------------------------------------------------------------------------
Wolfgang Bangerth                email:            [email protected]
www: http://*www.*math.tamu.edu/ ~bangerth/

_______________________________________________
dealii mailing list http://*poisson.dealii.org/mailman/listinfo/dealii


_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to