Revolve in PyGmsh leaving point of revolution in mesh

mesh
1

I am trying to run code to perform some analysis of an elastic arch, in doing this I have to build a circular arch using PyGmsh. I have done this with the following code:


with pygmsh.geo.Geometry() as geom:

    poly = geom.add_polygon(
        [
            [-100*np.sin(Theta),100*np.cos(Theta), -1/2],
            [-100*np.sin(Theta),100*np.cos(Theta), 1/2],
            [-101*np.sin(Theta),101*np.cos(Theta), 1/2],
            [-101*np.sin(Theta),101*np.cos(Theta), -1/2],
        ],
        mesh_size=1,
    )

    geom.revolve(poly,[0.0,0.0,-1/2],[0.0,0.0,1/2],angle=2*Theta)

    mesh = geom.generate_mesh(dim=3)

This generates the shape I want, but with the two points at [0,0,1/2] and [0,0,-1/2] still there. This makes it very difficult to view in paraview, and I’m not sure if it will effect the later calculations?

2

This is a PyGMSH question, please post it at: Issues · nschloe/pygmsh · GitHub

Please note that PyGMSH is no longer actively maintained, so I would strongly suggest using the GMSH Python API.
I’ve got a tutorial for this at: Using the GMSH Python API to generate complex meshes | Jørgen S. Dokken