The aim of this investigation is the proposal of 3D mesh generation method based on the Delaunay triangulation. The method is valid for the finite element modelling of any convex 3D domain into tetrahedra with optimum geometrical configuration. This paper includes the mathematical background of the mesh generation method, its detail, proposal of some efficient tools for faster and more rigorous computations, and some examples of the mesh generation.