In this paper, a multiphysics simulation workflow based on the use of reduced-order models (ROMs) to decrease the computation time is presented. Using FE pre-computed flux-linkage and lumped force lookuptables (LUTs) in function of dq0 currents and rotor position, an electromagnetic model of a permanent-magnet synchronous machine (PMSM) is developed. Coupled to a mechanical state-space representation, the complete workflow allows to predict the vibrations of the motor. Using our model, the influence of the 0-sequence current on the vibration spectrum is investigated. A detailed description on the parameter identification for homogenized stator core of the considered 12-slot 10-pole delta-connected fractional-slot concentrated-winding (FSCW) spoke-type PMSM is also given. The numerical results are validated with experimental data for the machine at hand. Because of the large amount of FE data required, the workflow is best suited for the optimization of motor control schemes (instead of the initial design optimization) with the purpose of reducing the noise, vibration, and harshness (NVH).