In view of the shortcomings of the existing calculation methods of electromagnetic forces, this paper proposes an analytical method (AM) combining the magnetomotive force (MMF)-permeance method and the improved subdomain method. The amplitude of each force density component considering stator teeth saturation and current harmonics is calculated, and the specific order and frequency components in the radial and tangential electromagnetic force densities are analyzed and sorted out. In modal analysis, the stator core and the casing are regarded as cylindrical shells with axial ribs, the natural frequencies of the stator are obtained based on the energy method, and the accuracy is well verified by the modal test. Finally, the frequency response functions of low-order modes of the stator only related to the stator structure parameters are determined, and the vibration accelerations under different working conditions are obtained based on the linear superposition principle by only updating the electromagnetic forces. Compared with the multi-physical coupling finite element analysis, the total calculation time of the proposed AM is reduced by more than 92%, which facilitates the fast prediction and effective reduction of electromagnetic vibrations in the motor design stage.