Hi,
I'm thinking about implementing membrane in triaxial test, in order
to be able to see localization due to internal shearing.
questions:
1) was it done already by someone, and I missed this?
2) it is maybe not necessary, because periodic boundaries also allow
to see localization due to shearing? (I doubt that, because the shift
on one side is on different height than the shift on right side).
I would do this in following way: detect which spheres are on the
outer boundary (so they are touching the membrane), and push them
inside with some force.
To detect spheres that touch the membrane I will do a Delaunay
triangulation of whole sample, with extra four points added. Each of
those four points is on the exterior of the sample and is a
"membrane point". Any sphere that gets connected through Delaunay
with this point must be lying on the boundary and is in contact with
membrane.
do you have any comments on that?
I suppose that eigen has Delaunay triangulation. Previously I was
using wm3 for that.