Abstract: Natural coordinates are used to define multibody systems. Dynamic equations of multibody systems in natural coordinates are first established by means of the nullspace of constraint Jacobian matrix. They are then derived in the form of generalized momenta. The process of the unmerical dynamic analysis is given. One illustrative example is finally presented.