Three-phase permanent magnet synchronous machine with sinusoidal or trapezoidal back electromotive force, or five-phase permanent magnet synchronous machine with sinusoidal back electromotive force
Simscape / Electrical / Specialized Power Systems / Electrical Machines
The Permanent Magnet Synchronous Machine block implements a three-phase or a five-phase permanent magnet synchronous machine. The stator windings are connected in wye to an internal neutral point.
The three-phase machine can have a sinusoidal or trapezoidal back EMF waveform. The rotor can be round or salient-pole for the sinusoidal machine. The rotor is round when the machine is trapezoidal. Preset models are available for the sinusoidal back EMF machine.
The five-phase machine has a sinusoidal back EMF waveform and round rotor.
The Permanent Magnet Synchronous Machine block operates in either generator or motor mode. The mode of operation is dictated by the sign of the mechanical torque (positive for motor mode, negative for generator mode). The electrical and mechanical parts of the machine are each represented by a second-order state-space model.
The sinusoidal model assumes that the flux established by the permanent magnets in the stator is sinusoidal, which implies that the electromotive forces are sinusoidal.
The trapezoidal model assumes that the winding distribution and flux established by the permanent magnets produce three trapezoidal back EMF waveforms.
These equations are expressed in the rotor reference frame (qd frame). All quantities in the rotor reference frame are referred to the stator.
$$\frac{d}{dt}{i}_{d}=\frac{1}{{L}_{d}}{v}_{d}-\frac{R}{{L}_{d}}{i}_{d}+\frac{{L}_{q}}{{L}_{d}}p{\omega}_{m}{i}_{q}$$
$$\frac{d}{dt}{i}_{q}=\frac{1}{{L}_{q}}{v}_{q}-\frac{R}{{L}_{q}}{i}_{q}-\frac{{L}_{d}}{{L}_{q}}p{\omega}_{m}{i}_{d}-\frac{\lambda p{\omega}_{m}}{{L}_{q}}$$
$${T}_{e}=1.5p[\lambda {i}_{q}+({L}_{d}-{L}_{q}){i}_{d}{i}_{q}]$$
L_{q}, L_{d} | q-axis and d-axis inductances |
R | Resistance of the stator windings |
i_{q}, i_{d} | q-axis and d-axis currents |
v_{q}, v_{d} | q-axis and d-axis voltages |
ω_{m} | Angular velocity of the rotor |
λ | Amplitude of the flux induced by the permanent magnets of the rotor in the stator phases |
p | Number of pole pairs |
T_{e} | Electromagnetic torque |
The L_{q} and L_{d} inductances represent the relation between the phase inductance and the rotor position due to the saliency of the rotor. For example, the inductance measured between phase A and B (when phase C is left open) is given by:
$${L}_{ab}={L}_{d}+{L}_{q}+\left({L}_{q}-{L}_{d}\right)\mathrm{cos}\left(2{\theta}_{e}+\frac{\pi}{3}\right),$$
where Θ_{e} represents the electrical angle.
The next figure shows the variation of the phase-to-phase inductance in function of the electrical angle of the rotor.
For a round rotor, there is no variation in the phase inductance:
$${L}_{d}={L}_{q}=\frac{{L}_{ab}}{2}.$$
For a salient round rotor, the dq inductances are given by:
$${L}_{d}=\frac{\mathrm{max}({L}_{ab})}{2}$$
$${L}_{q}=\frac{\mathrm{min}({L}_{ab})}{2}$$
These equations are expressed in the rotor reference frame using an extended Park transformation (q1d1 and q2d2 frame). All quantities in the rotor reference frame are referred to the stator.
$$\frac{d}{dt}{i}_{d1}=\frac{1}{L}{v}_{d1}-\frac{R}{L}{i}_{d1}+\frac{{L}_{q}}{L}p{\omega}_{m}{i}_{q1}$$
$$\frac{d}{dt}{i}_{q1}=\frac{1}{L}{v}_{q1}-\frac{R}{L}{i}_{q1}-\frac{{L}_{d}}{L}p{\omega}_{m}{i}_{d1}-\frac{\lambda p{\omega}_{m}}{L}$$
$$\frac{d}{dt}{i}_{d2}=\frac{1}{L}{v}_{d2}-\frac{R}{L}{i}_{d2}$$
$$\frac{d}{dt}{i}_{q2}=\frac{1}{L}{v}_{q2}-\frac{R}{L}{i}_{q2}$$
$${T}_{e}=2.5p\lambda {i}_{q1}$$
L | Armature inductance |
R | Resistance of the stator windings |
i_{q1}, i_{d1} | q1-axis and d1-axis currents |
v_{q1}, v_{d1} | q1-axis and d1-axis voltages |
i_{q2}, i_{d2} | q2-axis and d2-axis currents |
v_{q2}, v_{d2} | q2-axis and d2-axis voltages |
ω_{m} | Angular velocity of the rotor |
λ | Amplitude of the flux induced by the permanent magnets of the rotor in the stator phases |
p | Number of pole pairs |
T_{e} | Electromagnetic torque |
These equations are expressed in the phase reference frame (abc frame). Note that the phase inductance L_{s} is assumed to be constant and does not vary with the rotor position.
$$\begin{array}{c}\frac{d}{dt}{i}_{a}=\frac{1}{3{L}_{s}}\left(2{v}_{ab}+{v}_{bc}-3{R}_{s}{i}_{a}+\lambda p{\omega}_{m}\left(-2{\Phi}_{a}^{\prime}+{\Phi}_{b}^{\prime}+{\Phi}_{c}^{\prime}\right)\right)\\ \frac{d}{dt}{i}_{b}=\frac{1}{3{L}_{s}}\left(-{v}_{ab}+{v}_{bc}-3{R}_{s}{i}_{b}+\lambda p{\omega}_{m}\left({\Phi}_{a}^{\prime}-2{\Phi}_{b}^{\prime}+{\Phi}_{c}^{\prime}\right)\right)\\ \frac{d}{dt}{i}_{c}=-\left(\frac{d}{dt}{i}_{a}+\frac{d}{dt}{i}_{b}\right)\\ {T}_{e}=p\lambda \left({\Phi}_{a}^{\prime}\cdot {i}_{a}+{\Phi}_{b}^{\prime}\cdot {i}_{b}+{\Phi}_{c}^{\prime}\cdot {i}_{c}\right)\end{array}$$
L_{s} | Inductance of the stator windings |
R | Resistance of the stator windings |
i_{a}, i_{b}, i_{c} | a, b and c phase currents |
Φ_{a}^{'}, Φ_{b}^{'}, Φ_{c}^{'} | a, b and c phase electromotive forces, in per-unit value to the amplitude of the flux λ |
v_{ab}, v_{bc} | ab and bc phase to phase voltages |
ω_{m} | Angular velocity of the rotor |
λ | Amplitude of the flux induced by the permanent magnets of the rotor in the stator phases |
p | Number of pole pairs |
T_{e} | Electromagnetic torque |
The electromotive force Φ^{'} is represented by:
$$\begin{array}{c}\frac{d}{dt}{\omega}_{m}=\frac{1}{J}\left({T}_{e}-{T}_{f}-F{\omega}_{m}-{T}_{m}\right)\\ \frac{d\theta}{dt}={\omega}_{m}\end{array}$$
J | Combined inertia of rotor and load |
F | Combined viscous friction of rotor and load |
θ | Rotor angular position |
T_{m} | Shaft mechanical torque |
T_{f} | Shaft static friction torque |
ω_{m} | Angular velocity of the rotor (mechanical speed) |
When you use Permanent Magnet Synchronous Machine blocks in discrete systems, you might have to use a small parasitic resistive load, connected at the machine terminals, to avoid numerical oscillations. Large sample times require larger loads. The minimum resistive load is proportional to the sample time. Remember that with a 25 μs time step on a 60 Hz system, the minimum load is approximately 2.5% of the machine nominal power. For example, a 200 MVA permanent magnet synchronous machine in a power system discretized with a 50 μs sample time requires approximately 5% of resistive load, or 10 MW. If the sample time is reduced to 20 μs, a resistive load of 4 MW is sufficient.
The Permanent Magnet Synchronous Machine block assumes a linear magnetic circuit with no saturation of the stator and rotor iron. This assumption can be made because of the large air gap usually found in permanent magnet synchronous machines.
The power_brushlessDCmotor
example illustrates the use of the
Permanent Magnet Synchronous Machine block.
[1] Grenier, D., L.-A. Dessaint, O. Akhrif, Y. Bonnassieux, and B. LePioufle. “Experimental Nonlinear Torque Control of a Permanent Magnet Synchronous Motor Using Saliency.” IEEE^{®} Transactions on Industrial Electronics, Vol. 44, No. 5, October 1997, pp. 680-687.
[2] Toliyat, H.A. “Analysis and Simulation of Multi-Phase Variable Speed Induction Motor Drives Under Asymmetrical Connections.” Applied Power Electronics Conference and Exposition, Vol. 2, March 1996, pp. 586-592.
[3] Beaudart, F., F. Labrique, E. Matagne, D. Telteux, and P. Alexandre. “Control under normal and fault tolerant operation of multiphase SMPM synchronous machines with mechanically and magnetically decoupled phases.” International Conference on Power Engineering, Energy and Electrical Drives, March 2009, pp. 461-466.
[4] Krause, P.C., O. Wasynczuk, and S.D. Sudhoff. Analysis of Electric Machinery and Drive Systems. IEEE Press, 2002.