Linear mixed models have become popular in many statistical applications during recent years. However design issues for multi-response linear mixed models are rarely discussed. The main purpose of this paper is to investigate D-optimal designs for multi-response linear mixed models. We provide two equivalence theorems to characterize the optimal designs for the estimation of the fixed effects and the prediction of random effects, respectively. Two examples of the D-optimal designs for multi-response linear mixed models are given for illustration.