Methods of Modeling Permanent Magnet Synchronous Motor
Methods of Modeling Permanent Magnet Synchronous Motor
In dynamic modeling and simulation of permanent magnet synchronous motor for improved
performance, MATLAB software will be used to model and observe the dynamic behavior of the
Surface-Mounted and Interior PMSM under various operating condition. A generalize method of
solution and analysis by computer simulation using machine variable (voltage, current, or flux
linkage) by describing them either in rotor reference frame, stationary frame or synchronous
frame. This will enable the motor current, voltage, speed, power, torque characteristics to be
viewed under different operating condition and the results of relationship between them and time
will be showed.
Detailed modeling of Permanent Magnet motor drive system is required for proper simulation of
the system. The d-q model has been developed on rotor reference frame as shown in figure 3.1.
At any time t, the rotating rotor d-axis makes and angle θr with the fixed stator phase axis and
rotating stator MMF makes an angle α with the rotor d-axis. Stator MMF rotates at the same
The model of PMSM without damper winding has been developed on rotor reference frame
i. Saturation is neglected.
The equivalent circuit shown in figure 3.2 express about the dynamic equation of PMSM without
considering core loss. It should be notified that, all lower case symbols in introduce
V q =Rs i q + ρ ( λq ) +ωr λ d 1
V d =Rs i d + ρ ( λ d ) −ω r λq 2
Where the q and d-axes flux linkages in the rotor reference frame is
λ q=Lq i q 3
λ d=Ld i d + λm 4
V q =Rs i q + ρ ( Lq i q ) +ωr ( Ld i d + λm ) 5
V d =Rs i d + ρ ( Ld i d + λm ) −ωr ( Lq i q ) 6
Resolving further,
V q =(R ¿ ¿ s+ ρ Lq )i q+ ωr Ld i d + ωr λ m 7 ¿
V d =−ω r Lq i q + ( R s + ρ Ld ) i d 8
[ ][
V q (R ¿ ¿ s + ρ Lq)¿ ω r L d
Vd
=
¿
i
][
ω λ
]
( R s + ρ Ld ) ¿ q +[ r m ]9
id 0
λ m=λ af =Lm i fr 10
ρ ( λ m )=0
d
ρ :Operator
dt
λ m :The flux linkage due ¿ the rotor magnets linking the stator .
3 P 3 P
T e = ( ) ( λ d i q −λq i d ) = ( )( λm i q + ( Ld −Lq ) i q i d )11
2 2 2 2
3
Pem =ω rm T e = ωr ( λd i q− λq i d ) 12
2
P
ω r= ω 13
2 rm
Where:
P: Number of poles
T e =T l +T d + B ω rm+ J ρωrm 14
Solving for the rotor mechanical speed from equation 9, assuming the dry friction is equals zero,
T e −T l−B ω rm
ω rm=∫ (¿ )dt 15 ¿
J
and
2
ω rm= ω 16
P r
Where:
Here, the effect of core loss on the various performance characteristics equations is considered.
(Rs + R c ) (R s+ Rc )
V q =Rs i q + ρ ( λq)+ ω λ d 17
Rc Rc r
(R s + Rc ) (R s + Rc )
V d =Rs i d + ρ ( λd )− ωr λ q 18
Rc Rc
The expression of Te is the same as without core resistance, this implies that the effect of this
resistance does not appear directly in this expression but as its effect in l d also lq and
consequently on Te.