David,
I'm trying to compute the effective elastic properties of a heterogeneous,
linear and a bi-periodic system (i.e., left-right and top-bottom periodic
displacement fields). To this system, I would like to apply a global shearing
by prescribing the displacement field of the surfaces in the form of Dirichlet
BCs. This seems slightly contradictory since a bi-periodic system doesn't have
surfaces.
However, we can still think of the global shearing as an average surface
displacement around which periodic fluctuations occur (the origin of such
fluctuations is due to heterogeneous elastic properties). This is illustrated
in the picture below.
bitmap.png
I would like to know which is the best way to do this in deal.II (I have tried
with make_periodicity_constraints and interpolate_boundary_values, but the
problem is that, as I explain before, we set apparently contradictory
constraints).
It seems to me like the correct boundary values are of the form
u(left) = u(right) + offset
and similarly for the bottom/top. The point is that you encode the shearing in
the 'offset', so you have a variation of periodic boundary conditions that
includes this nonzero offset.
I believe that the make_periodicity_constraints() function takes such an
offset argument. Have you tried that?
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/55564a33-65b9-f0f6-eccb-d9e41b4c1f61%40colostate.edu.