Dear David, > I studied the step- 44 Tutorial Program and read already a lot about how to > apply a Neumann contribution. But I have still two questions: > In my case, I like to apply a force (not a force per unit reference area) as > a Neumann contribution on the RHS. But for the Integral, Ni * traction* JxW > I need the Piola- Kirchhoff traction. > My idea would be the following: F= PN dA = T dA <=> F/dA= T > I would simply divide the force by the area of the cellsurface. > So my first question is: Can I do it like this or is there a simpler way of > doing it?
It depends. The approach you’ve described here (that is, to specify the total force applied to a surface) is the same as is done in this code-gallery example (see line 1690): https://github.com/dealii/code-gallery/tree/master/Quasi_static_Finite_strain_Compressible_Elasticity <https://github.com/dealii/code-gallery/tree/master/Quasi_static_Finite_strain_Compressible_Elasticity> But this only makes sense if you know the area over which the traction is integrated and if the traction force is uniformly distributed. For a geometrically complex shapes, this may no longer be the best solution, and you should probably consider precomputing the total surface area before defining how the traction is distributed. This point relates to your second question... > And the second question is: How can I access this sort of Information? The > area might for example be computed using the norm of a cross product. But I > am not sure how to get the two vectors for this operation. You can just use numerical integration for this (which can also take account of the mapping used on curved boundaries). Look at the compute_vol_current() function in step-44, which can be easily adopted to compute surface areas. I hope this helps. Best, 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
