For general nonholonomically constrained systems, variation identity is used to define three kinds of unfree variations, i.e., nonholonomic variations: the vakonomic, the Hlder and the Suslov by means of vector fields on the extended configuration manifold. The relations among the three kinds variations are discussed and a necessary and sufficient condition for the variations to become free ones is obtained. The nonholonomic variations and the corresponding integral variational principles are utilized to derive the two kinds of dynamical equations: vakonomic equations and Routh's equations or Chaplygins equations. By comparing vakonomic equations with Rouths equations and Chaplygins equations respectively, two necessary and sufficient conditions for the two kinds of equations to have common solutions are obtained, which are not integrable conditions of the constraints.