We developed a functional autoencoder framework for cluster analysis of multi-dimensional functional data, data that are random realizations of vector-valued random functions.