I am trying to implement mixed elasticity with Neumann conditions using deal.ii. The formulation seeks for stress, displacement and rotation as the unknowns, and the Neumann conditions are to set stress normal components to be equal to some function g.
I know that for mixed Darcy it is possible to do this using project_boundary_values as from what I understood from the implementation of this method, it does not do anything to pressure DoFs as they are not on the boundary. However, in my setting I am using BDM1 - Q0 - Q1_cts for elasticity, and while this approach does not constrain displacements, it actually restricts rotations, which I would want to avoid. I looked at interpolate_boundary_values which takes ComponentMask as an argument, but unfortunately it does not work with the elements I am using. Could you please point me to the possible way of solving this issue? Another question, for mixed Darcy, is there a difference in using project_boundary_values or project_boundary_values_div_conforming dealing with Neumann conditions? Thank you. -- 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.
