From GridLAB-D Wiki
Jump to: navigation, search

Three Phase Induction Motors

Implement three phase induction motor dynamic phasor model into GridLAB-D.

Synopsis

This work will incorporate the ability to model three phase induction motor and other more advanced power flow capabilities into GridLAB-D.

Classes

Variables

Variable Definition
[math] V_{p} [/math] Positive sequence voltage
[math] V_{n}^{*} [/math] Complex conjugate of negative sequence voltage
[math] r_{s} [/math] Stator resistance
[math] L_{s}[/math] Stator inductance
[math] w_{s}[/math] Angular speed
[math] L_{m}[/math] Magnetizing inductance
[math] L_{r}[/math] Rotor inductance
[math] \Omega _{r,0}[/math] Mechanical speed dc component
[math] \Omega _{r,2} [/math] Mechanical speed second harmonic component
[math] I_{P}^{s} [/math] Positive sequence of stator current
[math] I_{n,s}^{*}[/math] Complex conjugate of negative sequence stator current
[math] I_{n,r}^{*}[/math] Complex conjugate of negative sequence rotor current
[math] J [/math] Inertia
[math] B [/math] Damping constant
[math] \psi _{f}^{R} + j. \psi _{f}^{I}[/math] Forward rotating flux
[math] \psi _{f}^{R} + j. \psi _{f}^{I}[/math] Backward rotating flux
[math] \phi [/math] Voltage phasor angle

Dynamic Phasor Equations

[math] V_{p} = [r_{s}+j.w_{s}.L_{s} +L_{s} \frac{\mathrm{d} }{\mathrm{d} t} ]I_{P}^{s}+ [j.w_{s}.L_{m} +L_{m} \frac{\mathrm{d} }{\mathrm{d} t}]I_{P}^{r} [/math]

[math]0 = [r_{r}+j.w_{s}.L_{r} +L_{r} \frac{\mathrm{d} }{\mathrm{d} t} ]I_{P}^{r}+ [j.w_{s}.L_{m} +L_{m} \frac{\mathrm{d} }{\mathrm{d} t}]I_{P}^{s}- j.\Omega_{r,0} \frac{\mathrm{P} }{\mathrm{2}} [L_{m} I_{P}^{s}+ L_{r}.I_{P}^{r}] - j.\Omega_{r,2} \frac{\mathrm{P} }{\mathrm{2}} [L_{m} I_{n,s}^{*}+ L_{r}.I_{n,r}^{*}][/math]

[math]V_{n}^{*} = [r_{s} - j.w_{s}.L_{s} + L_{s} \frac{\mathrm{d} }{\mathrm{d} t} ]I_{n,s}^{*}+ [-j.w_{s}.L_{m} +L_{m} \frac{\mathrm{d} }{\mathrm{d} t}]I_{n,r}^{*}[/math]

[math] 0 = [- j.w_{s}.L_{m} +L_{m} \frac{\mathrm{d} }{\mathrm{d} t}]I_{n,s}^{*} + [r_{r} - j.w_{s}.L_{r} + L_{r} \frac{\mathrm{d} }{\mathrm{d} t} ]I_{n,r}^{*} - j.\Omega_{r,0} \frac{\mathrm{P} }{\mathrm{2}} [L_{m} I_{n,s}^{*} + L_{r}.I_{n,r}^{*}] - j.\Omega_{r,2}^{*} \frac{\mathrm{P} }{\mathrm{2}} [L_{m} I_{P}^{s}+ L_{r}.I_{P}^{r}] [/math]

[math] J.\frac{\mathrm{d} }{\mathrm{d} t}.\Omega_{r,0} = \frac{\mathrm{P} }{\mathrm{2}}.L_{m}(I_{P}^{s}.I_{P}^{r} + I_{n,s}^{*}.I_{n,r}^{*})- B.\Omega_{r,0} - T_{L} [/math]

[math] J.\frac{\mathrm{d} }{\mathrm{d} t}.\Omega_{r,2} = \frac{\mathrm{P} }{\mathrm{4}}.L_{m}(I_{P}^{s}.I_{n}^{r} + I_{n,s}^{*}.I_{p,r})-( B + j. 2.J.w_{s}).\Omega_{r,2} [/math]


Simulation study

Case 1

Case 1 simulates a step change in mechanical torque at constant grid voltage. A torque step is applied at time t=10 sec. Fig 1. demonstrates constant quadrature axis voltage. Fig.2 indicates step change in torque. Following are the simulation results. Fig. 3 shows the impact of a step change in torque on direct axis current. Fig. 4 shows the impact of a step change in torque on quadrature axis current. Fig. 5 and Fig.6 show dc and second harmonic component of speed respectively.


Figure 2. Q axis voltage
Figure 2. Step change in torque
Figure 3. Direct axis current
Figure 4. Quadrature axis current
Figure 5. Speed
Figure 6. Second harmonic component of Speed

Case 2

Case 2 simulates a step change in grid voltage. A constant load torque is considered. Following are the simulation results. Fig. 9 shows impact of step change in voltage on direct axis current. Fig. 10 shows impact of step change in voltage on quadrature axis current. Fig.11 and Fig.12 shows dc and second harmonic component of speed respectively.


Figure 7. Q axis voltage
Figure 8. Step change in torque
Figure 9. Direct axis current
Figure 10. Quadrature axis current
Figure 11. Speed
Figure 12. Second harmonic component of Speed

References

P. Krause et al. Analysis of electric machinery and drive systems. Vol. 75. John Wiley & Sons, 2013.