Efficient optimization strategy of multibody systems is developed in this paper. Aug- mented Lagrange method is used to transform constrained optimal problem into unconstrained form firstly. Then methods based on second order sensitivity are used to solve the unconstrained problem, where the sensitivity is solved by hybrid method. Generalized-α method and generalized-α projection method for the differential-algebraic equation, which shows more efficient properties with the lager time step, are presented to get state variables and adjoint variables during the optimization procedure. Numerical results validate the accuracy and efficiency of the methods is presented.