On the large D expansion of Hermitian multi-matrix models