Most of our finite element codes are built using the Warp & Blend Lagrange elements. The node distribution of these elements are constructed using an edge warping and interior blending construction. The interior blending is tuned to optimize the node's interpolation property for up to degree 15 polynomial elements.
T. Warburton, An Explicit Construction for Interpolation Nodes on the Simplex, Journal of Engineering Mathematics, Volume 56, Number 3, pp. 247-262, 2006, (paper).