The dynamics of k coincident D-branes in string theory is described effectively by U(k) Yang-Mills theory at low energy. While these configurations appear as classical solutions in matrix models, it was not clear whether it is possible to realize the k =/= 1 case as the true vacuum. We consider massive Yang-Mills-Chern-Simons matrix model and investigate the generation of higher gauge groups U(2), U(3) in this model.