D-optimal designs for multi-response linear mixed models