Abstract-An interior-permanent-magnet motor is modeled by a combined analytical-numerical approach, in which the relationships between the stator currents and flux linkages are identified with static finite-element (FE) analysis. In addition to the previous approaches using the current space vector as the state variable, new models are also developed using the fluxlinkage space vector, which leads to more convenient timeintegration of the voltage equations. In order to account for the zero-sequence effects in delta connection, the models also include either the zero-sequence flux or current as an additional state variable. Finally, the possibilities of deriving the required quantities as partial derivatives of the magnetic field energy are discussed. The energy-based approaches avoid inaccuracies related to torque computation and thus allow better satisfying the power balance in the state-space model. We show the ability of the developed state-space models to predict the currents and torque equally to a nonlinear time-stepping FE model with much less computational burden. The results are validated by means of measurements for a prototype machine in both star and delta connections. In addition, we also demonstrate the effect of the zero-sequence current on the torque ripple in case of a deltaconnected stator winding.Index Terms-Field energy, finite element methods, magnetic saturation, permanent magnet machines, reluctance machines, state-space methods, torque ripple, variable speed drives.
I. INTRODUCTIONERMANENT-MAGNET (PM) synchronous machines with interior magnets have become popular in variablespeed applications as both motors and generators. Their main advantages are high power densities, relatively simple construction, small rotor losses as well as the possibility of taking advantage of the reluctance torque owing to the magnetic saliency [1]. However, the spatial permeance harmonics due to the saliency and the interaction of the permanent magnet flux with the stator slotting also cause unwanted ripple in the electromagnetic torque of the machine [2]. The torque ripple typically deteriorates the properties of the application and may excite mechanical resonances. In addition, the harmonic effects complicate the modeling of the machine. The voltage equations for a PM synchronous machine can be written asin which the column vectors u, i and ψ include either the phase-domain or two-axis components (in the stator frame of reference) of the voltage, current and flux linkage, respectively, and R s is the resistance of a stator phase. When solving this system, the relationship between the current and flux linkage is usually expressed using analytical inductance functions, the simplest ones of which reduce to constant inductances when (1) is transformed into the rotor directquadrature (d-q) frame of reference. This traditional d-q model only accounts for the fundamental spatial permeance variation and neglects magnetic saturation. Although more complicated functions for both higher-order permeance harmonics [3]- [6...