The fault diagnosis of the gearbox is primarily based on the analysis of the harmonic response. However, the generation mechanism of the harmonic response of the gear pair is still not distinctly explained in a healthy gearbox. Therefore, this paper is dedicated to establishing a dynamic model of the gearbox for this issue. In this model, the contact force of the bearing and the meshing stiffness of the gear pair are determined by the Hertz theory and the potential energy method, then the shaft motion is used as the connection to establish the coupling relationship between the gear pair and the bearing. After this dynamic model is confirmed by experiments and existing models, the coupling relationship inside the gearbox and the harmonic response of the gear pair are studied using this model. The result demonstrates that a closed control loop is formed due to the vibration transmission between the gear pair, the bearing and the shaft, and there is a strong coupled vibration between the gear pair and the bearing in the gearbox. Due to the generation of the center distance error, misalignment angle and disorganized meshing process of the gear pair, the meshing stiffness is altered, which excites the harmonic response of the gear pair. When the speed of the gearbox is increased, although the vibration amplitude of the gearbox is increased, the system predominantly vibrates at the fundamental frequency, and the harmonic of the gear pair is also shifted to the low frequency band.