This paper proposes a high-fidelity finite element model of a permanent synchronous motor (PMSM) to predict electromagnetic responses. The proposed method aims to generate electromagnetic responses from the PMSM under various operational conditions-including normal and faulty conditions-by coupling several partial differential equations governing the electromagnetics of a PMSM. The rotor eccentricity is considered to be a representative fault of a PMSM, which has electromagnetic characteristics that differ from the healthy state of a PMSM. Note that eccentricity is the most frequent fault during PMSM operation. Therefore, the proposed model could replicate the defected torque responses of an actual motor system. The effectiveness of the proposed model is validated using measurements from a PMSM test bench. Quantitative comparison reveals that the proposed model could replicate both the transient- and steady-state torque responses of the PMSM of interest at a variety of operational conditions, including a faulty status. The proposed model could be used to generate virtual electromagnetic responses of a PMSM, which could be used for data-driven fault detection methods of electric motor systems.