This study presents a new nonlinear dynamic model for a gear-shaft-disk-bearing system. A nonlinear dynamic model of a spur gear pair is coupled with linear finite element models of shafts carrying them, and with discrete models of bearings and disks. The nonlinear elasticity term resulting from backlash is expressed by a describing function, and a method developed in previous studies to determine multi harmonic responses of nonlinear multi-degree-of-freedom systems is employed for the solution. The excitations considered in the model are external static torque and internal excitation caused by mesh stiffness variation, gear errors and gear tooth profile modifications. The model suggested and the solution method presented combine the versatility of modeling a shaft-bearing-disk system that can have any configuration without a limitation to the total degree of freedom, with the accuracy of a nonlinear gear mesh interface model that allows to predict jumps and double solutions in frequency response. Thus any single stage gear mesh configuration can be modeled easily and accurately. With the model developed it is possible to calculate dynamic gear loads, dynamic bearing forces, dynamic transmission error and bearing displacements. Theoretical results obtained by using the method suggested are compared with the experimental data available in literature, as well as with the theoretical values calculated by employing a previously developed nonlinear single degree of freedom model.