3D constrained triangulations