Aspects of three‐dimensional constrained Delaunay meshing

Presented in this paper are the theoretical aspects of node addition to a non-convex, multiboundary mesh of tetrahedral elements as used in finite element modelling. The method used is derived from Watson1 and Shenton and Cendes2 and is extended to deal with node addition on inter-material boundaries. Several situations are identified that result in an illegal insertion polyhedron (IP), these could be caused by the ‘constrained’ nature of the mesh, adjacent objects with different material properties, or degenerate node configurations. A new Delaunay algorithm is described that checks for illegal cases of the IP and then corrects them, this checking relies on the consistent ordering of the element nodes. It is shown that a particular type of illegal IP can easily be identified and corrected using this technique. The Delaunay algorithm is then applied to automatic mesh generation, and modification to the basic Delaunay algorithm is described so that previously meshed edges and faces of the current object being meshed are not deleted during the addition of subsequent nodes. This ‘protection’ method only becomes viable by recognizing the node ordering sense of the IP faces.