A four-degree-of-freedom non-linear model with time varying mesh stiffness has been developed for the dynamic analysis of spur gear trains. The model includes a spur gear pair, two shafts, two inertias representing load and prime mover. In the model, developed several factors such as time varying mesh stiffness and damping, separation of teeth, teeth collision, various gear errors and profile modifications have been considered. Two computer programs are developed to calculate stiffness of a gear pair and transmission error and the dynamic analysis of modeled system using time integration method. Dynamic tooth and mesh forces, dynamic factors are calculated. Numerical examples have been given, which shows the time varying mesh stiffness ha a significant effect upon the dynamic tooth force and torsional vibrations.