Open access peer-reviewed chapter

# Induction Electrical Machine Simulation at Three-Phase Stator Reference Frame: Approach and Results

By Mikhail Pustovetov

Submitted: July 31st 2018Reviewed: July 30th 2019Published: September 23rd 2019

DOI: 10.5772/intechopen.88906

Downloaded: 352

## Abstract

This chapter provides the equations of mathematical model of three-phase induction electrical machine recorded in the three-phase stator reference frame. The possibilities of the proposed mathematical model are described in detail. This chapter also discusses the development of the computer model of an induction motor based on the abovementioned mathematical model. It uses an approach that allows combining during preparation the computer model dual methods: means of visual programming circuitry (in the form of electrical schematics) and logical one (in the form of block diagrams). The approach enables easy integration of the model of an induction motor as part of more complex models of electrical complexes and systems. The developed computer model gives the user access to the beginning and the end of a winding of each of the three phases of the stator and rotor. This property is particularly important when considering the asymmetric modes of operation or when powered by the special circuitry of semiconductor converters. Simulation results show the adequacy of the proposed mathematical model of induction electrical machine.

### Keywords

• induction electrical machine
• mathematical model
• stator reference frame
• electromagnetic torque
• instantaneous power
• simulation

## 1. Introduction

The induction motor (IM) plays a very important role in industrial sectors and transport, primarily due to its robustness and low cost. There are some authors, who devote their publications to problems of mathematical modeling of IM [1, 2, 3, 4, 5, 6]. The author wishes to offer its own version of a mathematical model (MM) of the IM, suitable for unbalanced modes of simulation purposes.

Writing the equations of the MM of a three-phase IM in a three-phase stator reference frame (SRF) is useful in the analysis, comparing the calculated and actual curves of currents and phase voltages, suitable without additional transformation of equations for the consideration of modes of IM operation at asymmetrical characteristics of feed or parameters of IM. Simulation results for the coordinatesα, β, γcorrespond to the actual processes in phases A, B, and C of the stator, in the case of the short-circuited rotor, i.e., in most cases, only and can be experimentally measured. IM simulation in three-phase coordinates is useful in detecting and diagnosing defects in the stator  and implementing algorithms for direct torque control.

It can be argued—subject to review only, the fundamental spatial harmonics of the magnetic field in the air gap is symmetrical in a construction three-phase IM, powered by non-sinusoidal voltages asymmetrical system—that simulation results will be correct, but the same one in two-phase orthogonal coordinates will not.

Advertisement

## 2. Mathematical model equations

The suggested MM of a three-phase IM is based on a three-phase electrical machine MM at SRF coordinate system axis α, β, γ, which are aligned with the stator phase axes A, B, and C. This basic MM supplemented by iron losses resistances rμincluded in magnetization circuits in each phase of IM in parallel with magnetizing inductance Lμ (analogy with the T-shaped equivalent circuit of induction electric machine). The equations of the electromagnetic processes in IM are given in Eq. (1), adopted by the usual generalized electrical machine assumptions: each of the phase stator windings creates in a smooth air gap in the sinusoidal-allocated magnetomotive force; magnetic saturation coefficient is constant. In the expressions (Eq. (1)), the following notations are further adopted: v, voltage; i, current; t, time; r, resistance; Ψ, magnetic flux linkage; ωr, mechanical rotor speed; and p, the number of pole pairs. Lower indexes α, β, γindicate the affiliation to the appropriate phase. The subscript sindicates the affiliation to the stator, the index rbelongs to the rotor, and index μbelongs to the magnetization branch. Lσsαis the leakage inductance of stator phase; Lσrαis the leakage inductance of the rotor winding phase.

Traditionally, the rotor parameters are given to the stator winding. Detailed components of the system (Eq. (1)) are described in Eqs. (2)(7). The form of Eq. (1) as much as possible is unified with published equations (Eq. (1)) of a three-phase transformer MM :

vsα=rsαisα+dLσsαisαdt+v0α;iμα_active=v0αrμα;vsβ=rsβisβ+dLσsβisβdt+v0β;iμβ_active=v0βrμβ;vsγ=rsγisγ+dLσsγisγdt+v0γ;iμγ_active=v0γrμγ;vrα=e0αerotαdLσrαirαdtrrαirα;vrβ=e0βerotβdLσrβirβdtrrβirβ;vrγ=e0γerotγdLσrγirγdtrrγirγ;E1

Voltages at the terminals of phase magnetization branches (derived from flux linkage of mutual induction)

v0α=dΨμαdt=rμαisα+irα12isβ+irβ+isγ+irγΨμαM;
v0β=dΨμβdt=rμβisβ+irβ12isα+irα+isγ+irγΨμβM;E2
v0γ=dΨμγdt=rμγisγ+irγ12isα+irα+isβ+irβΨμγM.

Mutual inductance of IM windings defined as

M=2/3LμE3

M—mutual inductance of the phase windings of the rotor and stator of IM in the case of coincidence of their axes. In other words main inductance of the winding of the stator when magnetic flux calculated in the absence of currents in other phases of the stator and rotor windings (i.e., from the part of main magnetic flux created by the stator phase itself).

Lμ—complete phase inductance of the stator winding from the main magnetic flux, which takes into account the presence of currents in other phases. In other words the inductance of the main magnetic flux part is created by the winding itself M, and the inductance of the portion of the main magnetic flux is created by two other stator windings M/2.

Per-phase magnetizing currents:

iμα=isα+irα=23isα12isβ+isγ+irα12irβ+irγ=iμα_active+iμα_reactive;
iμβ=isβ+irβ=23isβ12isα+isγ+irβ12irα+irγ=iμβ_active+iμβ_reactive;E4
iμγ=isγ+irγ=23isγ12isα+isβ+irγ12irα+irβ=iμγ_active+iμγ_reactive.

In the expressions (Eq. (4)), iμ_activeand iμ_reactiveare active and reactive (inductive) components of per-phase magnetizing current, respectively.

e0α=v0α;e0β=v0β;e0γ=v0γ.E5
erotα=ΨrβΨrγpωr3;erotβ=ΨrγΨrαpωr3;erotγ=ΨrαΨrβpωr3.E6

Flux linkages of IM rotor phases:

Ψrα=Lσrαirα+Ψμα;Ψrβ=Lσrβirβ+Ψμβ;Ψrγ=Lσrγirγ+Ψμγ.E7

The electromagnetic torque of IM equation is

Tem=p32Misαirγ+isβirα+isγirβisαirβ+isβirγ+isγirα.E8

The equation of motion for the IM shaft in the case of single-mass mechanical part:

dωrdt=1JTemTload,E9

where Jis the moment of inertia of masses coupling with the rotor shaft and Tloadis the static torque of mechanical load coupling with the IM shaft.

In the general case, Eq. (9) can be written as

dωrdt=1JTemTloadtωrΘr=TdynamicJ,E10

Tdynamic—dynamic torque on the IM shaft;

ωr=ωr0+0tTdynamicJdt,E11

ωr0—initial angular speed of the rotor shaft of IM.

The angle of rotation of the IM rotor shaft, radian

Θr=Θr0+0tωrdt,E12

Θr0—initial angular position of the IM rotor shaft, radian.

It is interesting to know the energy performance of the electric machine in the transition process. For IM in the absence of power from the rotor instantaneous value of the active power consumption is calculated as

P1=vsαisα+vsβisβ+vsγisγ.E13

The instantaneous reactive power consumption

Q1=13vsαisβisγ+vsβisγisα+vsγisαisβ.E14

The instantaneous value of useful shaft power of IM

P2=ωrTload.E15

The abovementioned equations (Eqs. (1)(15)) of IM MM may be supplemented by the expressions (Eq. (16)), which allow to go to the description of the rotor phase currents at frequency of f2=sf1, where sis the sleep of the rotor and f1is the frequency of stator voltage.

ira=irα13+23cospΘr+irβ13+23cospΘr2π3+ irγ13+23cospΘr+2π3;
irb=irα13+23cospΘr+2π3+irβ13+23cospΘr+ irγ13+23cospΘr2π3;E16
irc= irα13+23cospΘr2π3+ irβ13+23cospΘr+2π3+irγ13+23cospΘr.

Such a representation of the rotor currents demonstrates the nature of changes in the transition processes (see Figure 1 , which shows the results of the rotor phase current simulation in the case of balanced power supply sinusoidal voltages for IM type AE92-4, 4-pole, 40 kW). In Figure 1 , curve 1 is irγcurrent, curve 2 is irccurrent, and curve 3 is rotor speed. According to the frequency f2, you can check the value of the rotor speed of IM. Figure 1.The results of the rotor phase current simulation for IM type AE92-4.

Figure 2 shows the applicability of the developed MM for the simulation of electromagnetic and electromechanical processes in IM with its power supply from the power semiconductor converter—autonomous voltage inverter. Figure 2.The results of simulation for traction IM type DTA-350M (6-pole, 350 kW) under power supply from autonomous voltage inverter condition.

## 3. Requirements for computer model

In the process of the computer model of the IM according to the equations (Eqs. (1)(12)), development using OrCAD—computer-aided design system—intended primarily for the design and simulation of electronic and electrical devices , there was the task of ensuring the embedded model of IM in an electrical power supply circuit, including converters, by simple connection of virtual terminals. That is, relative to the model of the electrical circuit that is attached to this circuit, the IM model should also have the properties of electric circuits: one can apply a potential difference to the terminals, providing a bidirectional electric current in the connecting circuits, including a pass-through current between them. Another task is the formation of such structure of the computer model which will be used perhaps as a universal template, which records the values of variables (parameters of IM), asked simultaneously for all equations. The resulting computer model is suitable for the description of both the squirrel-cage and wound rotor IM in any of four quadrants. The winding phases can be connected in a triangle scheme, wye, independently joined with each of its voltage to have any other wiring of each other or of power source. The developed computer model gives the user access to the beginning and the end of a winding of each of the three phases of the stator and the rotor. This property is particularly important when considering the asymmetric modes of operation or when powered by the special circuitry of semiconductor converters .

## 4. The implementation of a computer model

The tasks are solved by combining two approaches in the IM computer model development: circuit for electrical parts and the method of block diagrams for the magnetic and mechanical parts. The computer model of IM is packed into a hierarchical block with a specified list of variables (parameters of IM). Within a single OrCAD project, using the copy operation, one can get the required number of hierarchical blocks (models of IM), for each of which it is possible to set unique parameter values. Sensors and signal inputs in the IM model used controlled sources of currents and voltages. Such voltage source is controlled by current (VSCC) and by voltage (VSCV), and current source is controlled by voltage (CSCV). A graphical image of the computer model of phase Awindings of the stator and the rotor of IM, composed of Eqs. (1)(7), is shown in Figure 3 . Figure 3.The graphical image of the computer model of phaseAwindings of the stator and the rotor of IM.

VSCC1 performs a function of the phase current of the stator isαsensor; similarly, VSCC2 is a sensor of rotor phase current irαgiven to the stator winding. CSCV1 is used for entering irαinto the scheme. VSCV1 is introduced into the phase magnetizing circuit v0αvoltage. A similar function is performed by VSCV2 entering into the circuit of the rotor winding; the difference is between the electromotive force of a branch of magnetization and rotation electromotive force, i.e., the voltage drop value. Active resistances rserv(see Figure 3 ) have a large value, for example, 10 MΩ. These resistors are introduced for service purposes. Without affecting numerical calculation results, they stabilize the solution (simulation) by maintaining the current circuit (the physical sense—the way for leakage current), which is especially important when a discrete change of the resistance of a IM circuit occurs, for example, when powered from the frequency converter or in the case of the phase circuit breaking. A similar solution is described in .

The rest of the equations of a three-phase IM mathematical model is implemented by the author in the computer model in the form of block diagrams . Figure 4 shows a part of the model, where one obtains instantaneous value of Lμ—complete phase inductance of the stator winding from the main magnetic flux in relative units. ΨμΣmratedis the rated value of ΨμΣm—the instantaneous value of the amplitude of the representing vector of the flux linkage of mutual induction. Figure 4.The block diagram for obtaining instantaneous values of complete phase inductance of the stator winding from the main magnetic flux in relative unitsLμ∗.

Signal ΨμΣmcan be calculated as

ΨμΣm=Ψμx2+Ψμy2,E17

where Ψμxand Ψμyare the projections of the representing vector of the flux linkage of mutual induction on orthogonal coordinate axes X and Y. In the case of three-phase stator reference frame of the coordinate system axis, these projections are derived from mutual induction flux linkages of each phase.

Ψμx=23Ψμα+Ψμβcos2π3+Ψμγcos2π3;E18
Ψμy=23Ψμβsin2π3+Ψμγsin2π3.E19

Figure 5 shows the piecewise linear approximation by five points of LμΨμΣfor IM type AGV250 (2-pole, 110 kW). Such a way successfully tested in IM simulator is developed by means of OrCAD PSpice : the one standard component table is used for programming the abovementioned approximated dependence. Figure 5.The piecewise linear approximation by five points ofLμ∗ΨμΣm∗for IM type AGV250.

In Figure 6 there is the model part corresponding to the abovementioned (Eqs. (1)(3)) equations. Figure 7 shows a part of the IM model, designed to determine the rotation electromotive force in each phase. Model rotor speed ωr_modelin ptimes is higher than the real mechanical speed of rotor ωr. Figure 8 shows a part of the IM model corresponding to equations (Eq. (2)). Figure 9 shows a part of the model, which forms the signals of electromagnetic torque Tem, angular speed ωr, and the angle of rotation Θrof the rotor. In Figure 9 , nris mechanical angular rotor speed, rotation per minute. Figure 6.The block diagram that implements the calculation of the instantaneous value of the amplitude of the representing vector of the flux linkage of mutual induction. Figure 7.Part of the IM model, designed to determine the rotation electromotive force in each phase. Figure 8.Block diagrams corresponding to equations (Eq. (2)). Figure 9.Block diagrams, which form the signals of electromagnetic torqueTem, angular speedωr, and the angle of rotation of rotorΘr.

## 5. An example of using a mathematical model of the motor when the power supply voltage is unbalanced

A good example of functionality for previously provided MM under unbalanced supply voltage maybe modeling of electromechanical phase splitter, made on the basis of a three-phase IM.

Electromechanical phase splitters are used in Russia and India  onboard of electric trains and electric locomotives for transform AC single-phase voltage in three-phase voltage to feed auxiliary electric drives with IM loaded by fans and air compressors. Figure 10 is a schematic diagram of the rotary phase splitter, where С1is operating and С2is start-up capacity. Figure 10.Schematic diagram of the rotary phase splitter.

Phase splitters are IM with a symmetrical or nonsymmetrical stator winding and no load (or low load) on the shaft. Phase splitter can be considered as combined single-phase IM and three-phase synchronous generator. In accordance with the terminology, adopted in India, the phase splitter is called the ARNO converter .

On contemporary freight AC electric locomotives of family “Ermak” in Russia, IM of NVA-55 type (4-pole, 55 kW) is used as a phase splitter. To drive the fans and compressor, the same type of IM is used. Let us make simulation of the phase splitter start without connecting electrical loads to check for phase splitting effect. Define С1=968uF and С2=2904uF. The torque of mechanical losses on the shaft (of the load) will take 28 Nm at the rotor speed of 1500 rotation per minute.

The simulated results of phase splitter start are shown in Figures 11 and 12 , where the graphs are 1, 2, 3—line voltages between phases A-B, B-C, C-A— and 4, rotor speed. Start-up capacity С2is switched off when the current value of the voltage between one of the wires of a single-phase supply and wire of the phase splitter stator winding phase, which is not connected to single-phase supply (see Figure 10 ), will be over the value of 300 volts—such a situation means the end of the rotor acceleration and ending the formation of a three-phase voltage system (signal 5 in Figure 11 ). Figure 11.Results of phase splitter start simulation. Figure 12.Results of phase splitter simulation at steady-state mode.

Figure 12 shows the results of simulation of phase splitting in steady state (at average rotor speed of phase splitter 1500 rotation per minute). In Figure 12 curves 6–8 are the phase currents A, B, and C. Table 1 presents the numerical results at steady-state mode. Obviously, the phase splitting effect is clearly demonstrated by the results of simulation.

Phase-to-phase voltages (V)Phase currents (A)
A-BB-CC-AABC
444380408131.9113.922.1
SplittedObtained from a single-phase transformer secondary windingSplittedWithout any connected electrical loads

### Table 1.

Results of phase splitter simulation at steady-state mode.

The proposed mathematical model of a three-phase IM and the method of its computer implementation confirmed their adequacy and effectiveness to address some of the problems of analysis and design of electrotechnical complexes and systems, having IM in its composition.

Advertisement

## 6. Conclusion

In conclusion, the author would like to identify ways of improving the suggested mathematical model of a three-phase IM [13, 14]. Priority for criticism might be the lack in the mathematical model accounting for the skin effect in the conductors of the rotor winding. In fact, this is a notable disadvantage since in some cases, the current displacement in the rotor is able to accelerate the transitional process of IM start-up in several times [15, 16].

There are two types of problems whose solution requires a mathematical model accounting for the skin effect in the conductors of the IM rotor winding: the correct description of transient processes at direct start-up and correct account of the losses in the rotor winding in the case of feeding from the source of non-sinusoidal voltage or current. When deciding tasks of the first type, the author sometimes puts active resistance and leakage inductive resistance of the rotor, which change depending on the slip. The solution of tasks of the second type requires the coordination of the complex resistance components of the rotor values to current frequency in it. The second approach is more versatile as it can be used for solution of both types of the abovementioned tasks. According to the author, to implement the second approach, the best way is to create the computer implementation of the mathematical model of a three-phase IM in the original coordinates.

## Download for free

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Mikhail Pustovetov (September 23rd 2019). Induction Electrical Machine Simulation at Three-Phase Stator Reference Frame: Approach and Results, Applied Electromechanical Devices and Machines for Electric Mobility Solutions, Adel El-Shahat and Mircea Ruba, IntechOpen, DOI: 10.5772/intechopen.88906. Available from:

### chapter statistics

352total chapter downloads

1Crossref citations

### More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

### Related Content

#### Applied Electromechanical Devices and Machines for Electric Mobility Solutions

Edited by Adel El-Shahat

Next chapter

#### CMOS Active Inductor and Its Applications

By Dhara Pinkesh Patel

#### Advanced Applications for Artificial Neural Networks

Edited by Adel El-Shahat

First chapter

#### Introductory Chapter: Artificial Neural Networks

By Adel El-Shahat

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.