Hi Tomás,
tomás méndez <tmsmendez(a)gmail.com> writes:
> My name is Tomas Mendez im an architecture PhD candidate in Turin. Sorry to
> contact you directly but i found no other way.
*cough*
http://lists.tiker.net/listinfo/meshpy
*cough*
:)
> I'm trying to generate a periodic(able) mesh with a series of holes using
> meshpy for FEM structural calculations.
>
> I use the meshpy.triangle.subdivide_facets() function and it works just
> fine, all of the edges have X ammount of subdivisions on them.
>
> Because the mesh is forFEM purposes i need to use the minimum angle
> function. This also works great.
>
> The problem comes when there is a hole that gets too close to the first or
> last edge of the mesh. Triangle generates new points to satisfy all of the
> requirements, and the points on the first and last edges no longer
> coincide, therefore the mesh is no longer periadicable.
>
> What do you suggest that i do?
> Is there a way to forbid new points and still satisfy the rest of my
> requirements?
Yes, there is a way to prohibit new points on the boundary, and that's
what I'd recommend doing:
http://documen.tician.de/meshpy/reference.html#meshpy.triangle.build
What you want is 'allow_boundary_steiner=False', I think. Of course then
you need to put enough points on the boundary in your initial mesh.
HTH,
Andreas
