Based on the high order variational integrators for dynamic equations of multibody systems, gradient-based optimization methods and genetic algorithm are developed for the general model of design optimization of multibody dynamics.Augmented Lagrange method and unconstrained methods based on second order sensitivity are compared with the modified genetic algorithm with adaptive crossover and mutation rates.