Voronoi tessellation has been employed to characterise material features in Finite Element Method (FEM) simulation, however, a poor mesh quality of the voronoi tessellations causes problems in explicit dynamic simulation of forming processes. Although centroidal voronoi tessellation can partly improve the mesh quality by homogenisation of voronoi tessellations, small features, such as short edges and small facets, lead to an inferior mesh quality. Further, centroidal voronoi tessellation cannot represent all real micro structures of materials because of the almost equal tessellation shape and size. In this paper, a density function is applied to control the size and distribution of voronoi tessellations and then a Laplacian operator is employed to optimise the centroidal voronoi tessellations. After optimisation, the small features can be eliminated and the elements are quadrilateral in 2D and hexahedral in 3D cases. Moreover, the mesh quality is significantly higher than that of the mesh generated on the original voronoi or centroidal voronoi tessellation. This work is beneficial for explicit dynamic simulation of forming processes, such as micro deep drawing processes.