Abstract: Dynamic optimization based on multibody system dynamics is a complex process, which composed of traditional analysis and simulation of multibody system dynamics, design sensitivity analysis and iterative processes of optimization. The mixed penalty function method of optimization is presented, which is based on general kinematical, dynamical analysis formulation, general objective functions and adjoint variable method. Finally, an example of simple spring/damper-mass system is given to test the method presented in this paper.