First by using the first kind of Lagrange's equations,this paper presented the dynamical equations of multibody system with friction constraints,and the generalized forces of friction forces were described in the matrix form.Then by using the Baumgarte's method and the augmentation approach,the differential-algebraic equations were converted into the matrix form of differential equation,which can improve the calculating efficiency.Finally the numerical results proved the method to be effective.