Tensor Network B-splines for high-dimensional function approximation