Published: 09 September 2021

Development of a mathematical model of frequency controlled induction electric drive

Nazira Darimbaeva1
Marat Mustafin2
Nurgul Almuratova3
Algazy Zhauyt4
1, 2, 3Department of Electrical Machines and Electric Drive, Almaty University of Power Engineering and Telecommunications, Almaty, Kazakhstan
4Department of Electronics and Robotics, Almaty University of Power Engineering and Telecommunications, Almaty, Kazakhstan
Corresponding Author:
Marat Mustafin
Views 277
Reads 112
Downloads 634

Abstract

A mathematical model of a frequency-controlled induction (asynchronous) electric drive is proposed, which makes it possible to study electromagnetic and electromechanical processes in an electric drive with any drive mechanism and takes into account the non-sinusoidality of the output voltage of the frequency converter. The block diagram of the algorithms of calculations performed in the “Mathcad” computer application, the results of modeling and their analysis are presented.

Development of a mathematical model of frequency controlled induction electric drive

1. Introduction

At present, the most massive consumers of electrical energy are electric drives based on Induction (asynchronous) electric motors (IM). At the design stages of new and modernization of existing electric drives with IM, an important place is occupied by mathematical modeling of the object [1]. The mathematical description of the electromagnetic and electromechanical processes of an induction electric drive has been known for a long time [2]. Different authors and at different times have used different approaches, methods of analysis, forms of mathematical description of stationary and dynamic modes of operation of induction drive systems [3]. However, the existing models of the electric drive with IM either do not take into account the features of the driven mechanism or designed for one type of mechanisms [5]. So, in, a model of an induction drive for centrifugal mechanisms with a quadratic dependence of the torque resistance on speed is proposed [6]. The authors propose the universalization of the mathematical model of the electric drive in order to use it for mechanisms with any torque resistance. In addition, in the known models, one of the assumptions was the sinusoidality of the supply voltage of the IM. In reality, the output voltage of frequency converters contains higher harmonic components [7-9], leading to distortion of the shape of the motor currents and the electromagnetic torque. The article proposes to supplement the asynchronous motors (AM) model in order to take into account the real form of voltages and currents in a frequency-controlled electric drive [10]. The model, created on the basis of the computer mathematical application “Mathcad”, and this model is quite universal, since it covers all options for IM control. By setting the necessary initial conditions [11], all types of transient processes (start, stop, reverse, load shedding and surge) are modeled, various forms of setting influences for their formation. Calculation of instantaneous values of currents, torques, and powers is accompanied by the calculation of integral characteristics - powers during the transient process, power factors and efficiency [12]. The complete system of initial differential equations describing the processes in the induction electric drive of centrifugal pumps consists of three groups: the equations of electromagnetic equilibrium (Kirchhoff), which describe the electromagnetic processes in the engine, the equations of electromechanical energy conversion, and the equations characterizing the features of the mechanical load on the motor shaft [13-15].

2. Materials and methods

Fig. 1 shows a block diagram of the mathematical model of an electric drive based on IM. The added blocks for calculating the mechanical load and the parameters of the frequency converter, which are discussed in this article, are highlighted here [16].

Fig. 1Block diagram of the mathematical model

Block diagram of the mathematical model

The model described in includes a mandatory block of electromechanical conversions in the drive [1]. In general, the electromechanical relations in the program are calculated by the components of the equation of motion of the electric drive:

1
M-MC=Jdωdt.

In Eq. (1), the electromagnetic torque developed by the engine is determined quite simply:

2
M=-32pnlmi1xi2xi1yl2y.

To determine the dependence Mcω, it is necessary to carry out a calculation and analysis of mechanical characteristics for each investigated mechanism. For example, under centrifugal loading, the moment of resistance of the MC mechanism is a complex function of speed [17]. In addition, for the numerical solution of the system of differential equations of the drive, it is necessary to write in an analytical form (calculate) the first derivative (dMC/dt value). In this case, various types of simplifications and assumptions are used that reduce the accuracy of the simulation results.

In [1], for instance, dMC/dt calculated as:

dMCdt=dMCdsdsdt,
dsdt=-1ω0dωdt=-1Jω03pnlm2i1yi2x-i1xi2y-MC.

In the proposed model, it has been proposed to unify this calculation by making it applicable for any dependence Mcω. The mechanical characteristic of the most mechanisms is a continuous “smooth” function without jumps and discontinuities and, therefore, this can be approximated by a polynomial function [18]. The computer mathematical application “Mathcad” has a number of built-in functions that simplify this solution. One of the is used in the developed model. Let MC change during speed regulation according to a nonlinear dependence, uncertain analytically. It is determined experimentally or according to the passport characteristics presented by the manufacturer in the form of graphs or Tables. Fig. 2 shows two such characteristics corresponding to a mechanism with a constant power on the motor shaft (Mc1) and with a “fan” load (Mc2).

Fig. 2Changes in the moment of resistance

Changes in the moment of resistance

In the “Calculation of mechanical load” block, the mechanical characteristic presented in the form of a table is converted into a polynomial regression of the form:

3
Mcω=b0 +b1ω+b2ω2+..+bnωn.

Next, the value of dMp/dω is calculated, which in our case has the form:

4
dMcω/dω=b1 +2b2ω+3b3ω2+..+nbnωn-1

Fig. 3 shows a fragment of the program and the results of the calculation in the “Calculation of mechanical load” block for an object with a mechanical characteristic Mc1.

Fig. 3Fragment of the program “Calculation of mechanical load”

Fragment of the program “Calculation of mechanical load”

In the program, the dependence Mcω is presented in the form of the ММ1 vector. To solve the problem, the program “regressX,Y,k” is used which built into Mathcad, which returns the vector “coeffs” containing the coefficients of the k degree polynomial, which best approximates the “cloud” of points with coordinates stored in vectors X and Y.

Graphical interpretation of the result is shown in Fig. 4. Tabular values of the mechanical characteristics (X, Y, indicated by crosses) are approximated with high accuracy by the polynomial Mpω. The calculated derivative dMpω is also shown here [19].

To obtain the characteristics of any other mechanism in the initial data of the program, it is enough to set the tabular data of its mechanical characteristics. In the figure below (see Fig. 5) the calculation results for the mechanism with the “fan” load Mc2 are shown [20].

Fig. 4Exact approximation of the polynomials Mpω and dMp(ω)

Exact approximation of the polynomials Mpω and dMp(ω)

Fig. 5Fan load moment of resistance

Fan load moment of resistance

In accordance with the proposed changes, the block for solving the system of differential equations of the mathematical model of an electric drive with IM takes the form:

5
di1xdt=1l1l2-lm2(-r1l2i1x+ωl1l2i1y-ωslm2iy1+r2lmi2x +ωl2lmi2y-ωl2lmsi2y+(l2u1x-lmu2x)),di1ydt=1l1l2-lm2(-ωl1l2i1x+ωslm2i1x-r1l2i1y-ωl2lmi2x +ωl2lmsi2x+r2lmi2y+(l2u1y-lmu2y)),di2xdt=1l1l2-lm2(r1lmi1x-ωl1lmi1y+ωsl1lmiy1-r2l1i2x -ωlm2i2y+ωl1l2si2y+(l1u2x-lmu1x)),di2ydt=1l1l2-lm2(ωl1lmi1x-ωsl1lmi1x+r1lmi1y-ωl1l2si2x +ωlm2i2x-r2l1i2y+(l1u2y-lmu1y)),dsdt=-3pnlm2Jωi1yi2x+3pnlm2Jωi1xi2y+MCJω,dMCdt=bnsn3pnlm2Jωi1yi2xMn+3pnlm2Jωi1xi2yMns-MCJωMn .

Fig. 6 and Fig. 7 show the dependences of the torque resistance Mc1t and Mc2t at the start of the IM, obtained on this model.

The curves of the IM electromagnetic torque at start-up with various types of loads, obtained on the proposed model, are presented in Fig. 8 and Fig. 9.

Fig. 6Dependence of the moment of resistances Mc1t at the state of the IM

Dependence of the moment of resistances Mc1t at the state of the IM

Fig. 7Dependence of the moment of resistances Mc2t at the start of the IM

Dependence of the moment of resistances Mc2t at the start of the IM

Fig. 8Electromagnetic moment of IM at start with MC1

Electromagnetic moment of IM at start with MC1

Fig. 9Electromagnetic moment of IM at start with MC2

Electromagnetic moment of IM at start with MC2

3. Results and discussion

The obtained curves correspond to the concept of transient processes in an electric drive at the torque resistance of the mechanism, depending on the speed of the motor shaft, as well as to the mechanisms considered in the article (see Fig. 2). The results of simulating a drive with a centrifugal load (Fig. 7-9) coincide with the results obtained in on a model of an induction electric drive of a centrifugal pump [1]. In the previous model, one of the assumptions was the sinusoidality of the supply voltage of the IM. In reality, the output voltage of frequency converters contains higher harmonic components, leading to distortion of the shape of the motor currents and the electromagnetic torque [2]. The article proposes to supplement the IM model in order to take into account the real form of voltages and currents in a frequency-controlled electric drive. For this, the system is supplemented with equations that reflect the output parameters and properties of power supplies [5]. The block diagram adds two blocks (“Frequency converter parameters” and “Calculation of harmonics of the inverter voltage”), in which the harmonic analysis of the inverter output voltage is carried out and the obtained values are converted into synchronous coordinates x, y. In the calculations, we use the basic concepts of the vector-harmonic method for analyzing electromagnetic processes in IM. It is known that the higher harmonic components of voltages and currents can, like the first harmonics, be represented by vectors rotating with a frequency that is a multiple of the order of harmonics, and, therefore, they can be transformed in different coordinate systems. This makes it possible to use the same algorithms and formulas for calculating the distorting components of motor currents as for calculating electromagnetic processes in IM by the first harmonics. For ν-х harmonic components of IM currents and voltages, the equation of electrical equilibrium Eq. (1) retains the structure:

6
U1xυ U1yυ 00=r1+p(l1+lm)-ωυ(l1+lm)plm-ωυlmωυ(l1+lm)r1+pn(l1+lm)ωυlmpnlmpnlm-ωυsνlmr2pn(l2+lm)-ωυsν(l2+lm)ωυsνlmpnlmωυsν(l2+lm)r2(l2+lm)I1xυ I1yυ I2xυ I2yυ ,

where ν – is the serial number of the harmonic component of current and voltage; U1xν – is the projection of the ν-th harmonic of stator voltage on the x-axis; U1yν – is the projection of the ν-th harmonic of the stator voltage onto the y-axis; I1xν, I1yν – is the projections of the ν-th harmonic of the stator current, respectively, on the x and y axes; I2xν, I2yν – is the projections of the ν-th harmonic of the rotor current, respectively, on the x and y axes; l1, l2 – is the phase leakage inductance of the stator and rotor; lm – is the inductance of the magnetizing circuit; ω – is the angular frequency of rotation of the rotor; r1, r2 – is the active resistances of the stator and rotor phases; ων – is the rotation frequency ν-th harmonic component.

The difference from the first harmonic equation is as follows [1]. In Eq. (6), the value of sν appears – slip relative to higher harmonic ones:

7
Sν=ων±ωων.

In Eq. (7), the minus sign refers to the harmonic components of the positive sequence (including the first harmonic), the plus sign refers to the harmonic negative sequence, the vector of which rotates in the opposite direction.

In the equations of electromechanical energy conversion, the time derivative of sν:

8
dsνdt=-1ωνdωdt=-1Jων3pnlm2i1yνi2xν-i1xνi2yν-MC.

In the program for calculating the dynamic characteristics of the centrifugal mechanism drive, the quadratic dependence of the moment of resistance on the rotation frequency was used.

One more equation in the system (dMC/dt value) is obtained by differentiating the equation of motion of the electric drive with a quadratic dependence of the moment of resistance on the rotation frequency [12]:

9
dMCdt=dMCdsνdsνdt=-2Mn+2Mnsν-1Jων3pnlm2i1yνi2xν-i1xνi2yν-MC.

As a result, taking into account expressions Eqs. (6-9), the system of equations of an asynchronous electric drive takes the form:

10
di1xνdt=1l1l2-lm2(-r1l2i1xν+ωνl1l2i1yν-ωνsνlm2iy1ν+r2lmi2xν +ωνl2lmi2yν-ωνl2lmsνi2yν+(l2u1xν-lmu2xν)),di1yνdt=1l1l2-lm2(-ωνl1l2i1xν+ωνsνlm2i1xν-r1l2i1yν-ωνl2lmi2xν +ωνl2lmsνi2xν+r2lmi2yν+(l2u1yν-lmu2yν)),di2xνdt=1l1l2-lm2(r1lmi1xν-ωνl1lmi1yν+ωνsνl1lmiy1ν-r2l1i2xν -ωνlm2i2yν+ωνl1l2sνi2yν+(l1u2xν-lmu1xν)),di2yνdt=1l1l2-lm2(ωνl1lmi1xν-ωνsνl1lmi1xν+r1lmi1yν-ωνl1l2sνi2xν +ωνlm2i2xν-r2l1i2yν+(l1u2yν-lmu1yν)),dsνdt=-3pnlm2Jωνi1yνi2xν+3pnlm2Jωνi1xνi2yν+MCJων,dMCdt=3pnlmJωνi1yνi2xνMn-3pnlmJωνi1yνi2xνMnsν-3pnlmJωνi1xνi2yνMn +3pnlmJωνi1xνi2yνMnsν-2MCJωνMn+2MCJωνMnsν.

Figs. 10-14 show the results of calculations of transient processes taking into account the nonsinusoidality of the output voltage of the inverter by the example of starting an asynchronous motor with a squirrel-cage rotor.

Fig. 10Active and reactive components of the IM stator current (first harmonics)

Active and reactive components of the IM stator current (first harmonics)

Fig. 11a) Speed and b) torque IM (first harmonics)

a) Speed and b) torque IM (first harmonics)

a)

a) Speed and b) torque IM (first harmonics)

b)

The curves of the transient process for the first harmonics (see Figs. 2-3) are similar to the results obtained in [1], which confirms the adequacy of the model. For calculations for higher harmonic voltages and currents, a frequency converter with a controlled rectifier, made according to a 6-pulse (bridge) circuit, is adopted [13-15]. The output voltage of such an inverter contains higher harmonic components of the order of: ν=6n±1, where n= 1, 2,….

Figs. 4-6 show the results of modeling of transient processes for the seventh harmonics. The frequency and attenuation of the components of the higher harmonics of the stator current of the IM correspond to the well-known concepts of electrical engineering.

Fig. 12Active and reactive components of the AM stator current (seventh harmonic)

Active and reactive components of the AM stator current (seventh harmonic)

Fig. 13 shows the curve of the change in the electromagnetic torque of the IM from the interaction of the seventh harmonics of currents, calculated by the formula:

11
M7=32pnLmi1y7i2x7-i1x7i2y7.

Fig. 13Electromagnetic moment of IM from the interaction of seventh harmonics

Electromagnetic moment of IM from the interaction of seventh harmonics

Obviously, the value of such a torque (1.1 Nm) is negligible (the nominal torque of the IM in the model is 1200 Nm). But in the engine, as a result of the interaction of harmonic flux and currents of different orders, pulsating moments are formed [16]. Their main source is the interaction of the main flow with the higher harmonic components of the rotor current:

M1,ν=Ф-1i-2ν,

or for the interaction of 1st and 7th harmonics:

12
M1,7=32pnLmi1yi2x7-i1xi2y7.

The curve of the IM electromagnetic moment calculated according to Eq. (12) is shown in Fig. 6.

Fig. 14Electromagnetic moment of IM from the interaction of the first and seventh harmonics

Electromagnetic moment of IM from the interaction of the first and seventh harmonics

The amplitude of oscillations of the electromagnetic moment is 115 Nm, i.e. up to 9 % of the nominal value [17, 18]. The total pulsating moment from the interaction of the main flow with all harmonic currents of the rotor can be determined as conclusions:

M1,ν=32pnLmν=1i1yi2xν-i1xi2yν.

As a result of the further development of the mathematical model of an adjustable induction electric drive, it becomes possible to use the model to study electric drives of mechanisms with any moment of resistance with minimal preparation of data on the load [19]. At the same time, the adequacy of the model is preserved with a high accuracy of representation of the real mechanical characteristics of the driven mechanism, obtained experimentally or from the manufacturer, by a polynomial [20]. The model allows, in addition to options for IM control, types of transient processes, the nature of the load, to take into account the features of the circuit and the principle of operation of the frequency converter.

4. Conclusions

The mathematical model proposed in the article, in addition to the blocks describing electromagnetic and electromechanical processes in an asynchronous motor, includes blocks for calculating the mechanical load and parameters of the frequency converter. The frequency converter (FC) parameter block calculates the harmonic composition of the FC output voltage and prepares data for calculating the highest harmonic components of the IM currents. As a result, the currents and the electromagnetic moment of the IM are calculated taking into account the real (non-sinusoidal) form of the FC voltage. In this case, the calculation is possible for any circuit solution and method of forming the FC voltage. The proposed upgrade allows not only to refine the calculation of the dynamic and static characteristics of the electric drive, but also to assess the impact of the properties of the FC on its energy performance. Taking into account the change in the mechanical load on the motor shaft when adjusting its speed presents significant difficulties. In the proposed model, it is proposed to unify the mathematical description of the mechanical characteristics of the working mechanism. To do this, in the mechanical load calculation block, the entered tabular load data is converted into the coefficients of a polynomial of a given degree, which most accurately approximates the mechanical characteristic of any complexity. The obtained result is further used in solving the equation of motion of the electric drive. As a result, it becomes possible to use the model to study electric drives of mechanisms with any moment of resistance with minimal preparation of load data. Thus, a mathematical model of an adjustable asynchronous electric drive is obtained, which allows, in addition to options for controlling the IM, types of dynamic and stationary processes, to adequately take into account the features of the circuit and the principle of operation of the frequency converter, as well as the features of the mechanical load of any complexity.

References

  • Mustafin M. and Almuratova N., “Calculation of transient processes of electric drives of centrifugal mechanisms,” International Journal of Pharmacy and Technology, Vol. 8, No. 3, pp. 17745–17756, 2016.
  • L. S. Bocii and V. Muller, “Mathematic model of three-phase induction machine connected to advanced inverter for traction system for electric trolley,” Anais da Academia Brasileira de Ciências, Vol. 85, No. 2, pp. 849–858, May 2013, https://doi.org/10.1590/s0001-37652013005000028
  • A. P. Mark, R. Vairamani, and G. C. R. Irudayaraj, “Mathematical modeling and analysis of different vector controlled CSI fed 3-phase induction motor drive,” Journal of Applied Mathematics, Vol. 2014, pp. 1–13, 2014, https://doi.org/10.1155/2014/623982
  • K. P. Phillips, “Current-source converter for AC motor drives,” IEEE Transactions on Industry Applications, Vol. IA-8, No. 6, pp. 679–683, Nov. 1972, https://doi.org/10.1109/tia.1972.349851
  • M. Veerachary, “Optimal control strategy for a current source inverter fed induction motor,” Computers and Electrical Engineering, Vol. 28, No. 4, pp. 255–267, Jul. 2002, https://doi.org/10.1016/s0045-7906(00)00051-3
  • M. Salo and H. Tuusa, “A vector-controlled PWM current-source-inverter-fed induction motor drive with a new stator current control method,” IEEE Transactions on Industrial Electronics, Vol. 52, No. 2, pp. 523–531, Apr. 2005, https://doi.org/10.1109/tie.2005.844243
  • M. Babaei and H. Heydari, “Direct torque control of pulse width modulation current source inverter-fed induction motor by novel switching method,” Electric Power Components and Systems, Vol. 38, No. 5, pp. 514–532, Mar. 2010, https://doi.org/10.1080/15325000903376933
  • D. G. Holmes, B. P. Mcgrath, and S. G. Parker, “Current regulation strategies for vector-controlled induction motor drives,” IEEE Transactions on Industrial Electronics, Vol. 59, No. 10, pp. 3680–3689, Oct. 2012, https://doi.org/10.1109/tie.2011.2165455
  • O. Kuznyetsov, “Mathematical model of a three-phase induction machine in a natural abc reference frame utilizing the method of numerical integration of average voltages at the integration step and its application to the analysis of electromechanical systems,” Mathematical Problems in Engineering, Vol. 2019, pp. 1–13, Dec. 2019, https://doi.org/10.1155/2019/4581769
  • A. Dineva et al., “Review of soft computing models in design and control of rotating electrical machines,” Energies, Vol. 12, No. 6, p. 1049, Mar. 2019, https://doi.org/10.3390/en12061049
  • D. Wang, X. Yuan, and M. Zhang, “Power-balancing based induction machine model for power system dynamic analysis in electromechanical timescale,” Energies, Vol. 11, No. 2, p. 438, Feb. 2018, https://doi.org/10.3390/en11020438
  • A. Parizad, S. Mohamadian, M. E. Iranian, and J. M. Guerrero, “Power system real-time emulation: a practical virtual instrumentation to complete electric power system modeling,” IEEE Transactions on Industrial Informatics, Vol. 15, No. 2, pp. 889–900, Feb. 2019, https://doi.org/10.1109/tii.2018.2837079
  • M. Dagbagi, A. Hemdani, L. Idkhajine, M. W. Naouar, E. Monmasson, and I. Slama-Belkhodja, “ADC-based embedded real-time simulator of a power converter implemented in a low-cost FPGA: application to a fault-tolerant control of a grid-connected voltage-source rectifier,” IEEE Transactions on Industrial Electronics, Vol. 63, No. 2, pp. 1179–1190, Feb. 2016, https://doi.org/10.1109/tie.2015.2491883
  • M. D. Omar Faruque et al., “Real-time simulation technologies for power systems design, testing, and analysis,” IEEE Power and Energy Technology Systems Journal, Vol. 2, No. 2, pp. 63–73, Jun. 2015, https://doi.org/10.1109/jpets.2015.2427370
  • L. Wang et al., “Methods of interfacing rotating machine models in transient simulation programs,” IEEE Transactions on Power Delivery, Vol. 25, No. 2, pp. 891–903, Apr. 2010, https://doi.org/10.1109/tpwrd.2009.2039809
  • M. Y. Pustovetov, “Approach to computer implementation of mathematical model of 3-phase induction motor,” IOP Conference Series: Materials Science and Engineering, Vol. 327, p. 022085, Mar. 2018, https://doi.org/10.1088/1757-899x/327/2/022085
  • M. Sokola and E. Levi, “A novel induction machine model and its application in the development of an advanced vector control scheme,” The International Journal of Electrical Engineering and Education, Vol. 37, No. 3, pp. 233–248, Jul. 2000, https://doi.org/10.7227/ijeee.37.3.3
  • S. E. Zirka, Y. I. Moroz, and C. M. Arturi, “Accounting for the influence of the tank walls in the zero-sequence topological model of a three-phase, three-limb transformer,” IEEE Transactions on Power Delivery, Vol. 29, No. 5, pp. 2172–2179, Oct. 2014, https://doi.org/10.1109/tpwrd.2014.2307117
  • Maksimkina J., “Direct start-up of the squirrel-cage induction motor at variable parameters of the rotor,” Latvian Journal of Physics and Technical Sciences, Vol. 50, No. 1, pp. 10–21, 2013.
  • Benecke M., Doebbelin R., Griepentrog G., and Lindemann A., “Skin effect in squirrel cage rotor bars and its consideration in simulation of non-steady-state operation of induction machines,” in Progress in Electromagnetics Research Symposium Proceedings, pp. 1451–1455, 2011.

About this article

Received
12 May 2021
Accepted
18 August 2021
Published
09 September 2021
SUBJECTS
Dynamics and oscillations in electrical and electronics engineering
Keywords
mathematical modeling
induction motor
characteristic
frequency
energy