CC BY
9
UDC 621.313.33
doi: 10.20998/2074-272X.2020.3.02
V.S. Malyar, O.Ye. Hamola, V.S. Maday
MODELLING OF DYNAMIC MODES OF AN INDUCTION ELECTRIC DRIVE AT PERIODIC LOAD
Goal. Development of methods and mathematical models, based on them, for the calculation of transients and steady-state modes of induction electric drives operating in periodic load mode. Methodology. The developed algorithms are based on a mathematical model of an induction motor, which takes into account the saturation of the magnetic core and the displacement of current in the rotor bars. The processes are described by a system of nonlinear differential equations in the orthogonal axes x, y, which enables the results to be obtained with the smallest amount of calculations. The magnetization characteristics by the main magnetic flux and the leakage fluxes are used to calculate the electromagnetic parameters of the motor. To account for the current displacement in the rotor bars, the short-circuited winding is considered as a multilayer structure formed by dividing the bars in height by several elements. Results. Due to the variable load on the motor shaft, electromagnetic processes in both transient and steady state modes of the electric drive in any coordinate system are described by a system of nonlinear differential equations. The result of the calculation of the transients is obtained as a result of their integration time dependencies of coordinates (currents, electromagnetic torque, etc.) at a given law of change of the moment of loading. The proposed method of calculating steady-state mode is based on algebraization of differential equations on the mesh of nodes of the process cyclicity period and allows to obtain periodic dependencies in the time domain. Originality. The problem of calculating a steady-state periodic mode is solved as a boundary problem for a system of first-order differential equations with periodic boundary conditions, which allows to obtain instantaneous dependences during the period of currents, electromagnetic torque, capacities and other coordinates. Practical significance. Using the developed algorithm, it is possible to calculate the static characteristics of periodic processes as dependencies on different parameters of the cycle of periodic load or other coordinates, which is the basis for the choice of the motor for overload, power, heating, etc., as well as to detect the possibility of resonance. References 9, figures 4. Key words: induction motor, periodic load, mathematical model, steady-state dynamic mode, transient, static characteristics, saturation of the magnetic core, displacement of current.
Розроблено математичт модел i алгоритми, з використанням яких складен програми розрахунку nepexidHux процеав i усталених рeжuмiв асинхронних eлeктропрuводiв, як працюють в рeжuмi перюдичноТ змши навантаження. В Тх основу покладено математичну модель асинхронного двигуна, розроблену на основi теори кт i зображувальних вeкторiв електричних координат, в якш враховуеться насичення магнтопроводу i виткнення струму в стержнях короткозамкненого ротора. Внаслдок змнного навантаження на валу двигуна eлeктромагнiтнi процеси як в пeрexiднux, так i усталених режимах в будь-якш систем координат описуються системою немншних дuфeрeнцiальнux рiвнянь. В роботi використано систему ортогональних координатних осей x, y, яка обертаеться з довтьною швидшстю. Для обчислення електромаг^тних парамeтрiв двигуна використовуються характеристики намаг^чування основним магнiтнuм потоком, а також потоками розсювання статора i ротора. Для урахування виткнення струму в стержнях ротора короткозамкнена обмотка подаеться у вuглядi багатошаровоТ структури, утвореноТ розбиттям стержтв по вuсотi на ктька eлeмeнтiв. Усталений перюдичний режим розраховуеться методом розв'язування крайовоТзадачi, розробленим на основi апроксимаци координат кубiчнuмu сплайнами, що дае змогу отримати Тх перюдичн залежностг в позачасовш областi i розраховувати статuчнi характеристики як залeжностi в'д парамeтрiв циклу перюдично-змжного навантаження або нших координат. Бiбл. 9, рис. 4. Ключовi слова: асинхронний двигун, перюдичне навантаження, математична модель, усталений динамiчний режим, перехвдний процес, крайова задача, резонанс, статичш характеристики, насичення магштопроводу, виткнення струму.
Разработаны математические модели и алгоритмы, с использованием которых составлены программы расчета переходных процессов и установившихся режимов асинхронных электроприводов, которые работают в режиме периодического изменения нагрузки. В их основу положено математическую модель асинхронного двигателя, разработанную на основе теории цепей и изображающих векторов электрических координат, в которой учитывается насыщение магнитопровода и вытеснение тока в стержнях ротора. Вследствие переменной нагрузки на валу двигателя электромагнитные процессы как в переходных, так установившихся режимах в любой системе координат описываются системой нелинейных дифференциальных уравнений. В работе используется система ортогональных осей x, y, которая вращается с произвольной скоростью. Для вычисления электромагнитных параметров двигателя используются характеристики намагничивания основным магнитным потоком, а также потоками рассеивания статора и ротора. Для учета вытеснения тока в стержнях ротора короткозамкнутая обмотка представляется в виде многослойной структуры, образованной разделением стержней по высоте на несколько элементов. Установившийся периодический режим рассчитывается методом решения краевой задачи, разработанным на основе аппроксимации координат кубическими сплайнами, что дает возможность получить периодические зависимости во вневременной области и рассчитать статические характеристики как зависимости от параметров цикла периодически изменяющейся нагрузки или других координат. Библ. 9, рис. 4. Ключевые слова: асинхронный двигатель, периодическая нагрузка, математическая модель, установившийся динамический режим, переходный процесс, статические характеристики, резонанс, насыщение магнитопровода, вытеснение тока.
Introduction. In modern conditions of development of science and technology, the problem of the development of induction electric drives requires new
approaches to their practical implementation, which can be realized only on the basis of the development of
© V.S. Malyar, O.Ye. Hamola, V.S. Maday
adequate mathematical models of electric drive systems that are adapted to their operating conditions. Their use allows not only to correctly select the necessary induction motor (IM), but also to develop a control system under which the motor, operating under these conditions, would ensure the maximum possible efficiency of the electric drive system as a whole.
Modern factory methods make it possible to design an IM that is highly likely to meet the technical conditions of operation in a steady nominal mode with constant load. Such calculations are usually performed using classical substitution circuits [1, 2], but classical substitution circuits are not suitable for the calculation of dynamic modes, and their various adaptations need to be checked on a case-by-case basis.
In the practice of operation, IMs are used not only for actuating mechanisms that operate with unchanged mechanical torque of loading, but also for drives with periodic repeated short-term loading [3, 4]. The duration of the cycle of periodic re-alternating load T consists of two parts: the duration of the action of the load pulse and the pause. In particular, for repeated short-term operation (S3), the duration of the load pulse is expressed as a percentage to the duration of the full cycle. Switch-on time (ST) = 15; 25; 40; 60 % (e.g. S3 - 25 %; S3 - 40 %) is considered standard, with a cycle time of 10 minutes [5]. The industry produces IMs to operate in different standard-defined S3 modes. The selection of the motor power for the repeated short-term operation mode S3 can be made for equivalent power or torque for a given load schedule. Knowing the power of the catalog IM for motors designed to operate in a particular S3 mode, it is possible to select a motor to check for starting torque, overload capacity and heating [4].
Both standard long-run motors and motors specially designed for repeated short-term mode can operate in the repeated short-term mode. After all, often the values of the duration of switching on the IM are not standard. There is a need for a comprehensive study of the operation of the motor under the conditions specified by the working mechanism of the periodic load torque, which can be accomplished through mathematical modelling.
The goal of the work is to develop mathematical models for the analysis of the dynamic modes of induction motors operating under conditions of periodic-variable loading.
A mathematical model for the calculation of transients. For the analysis of the operation of electric drives operating in dynamic modes, mathematical models of IM, built on the basis of substituting circuits or linear differential equations (DEs) can be used only for approximate calculations. Because the electromagnetic torque is determined by the flux linkages and currents of the circuits of the motor, the inaccuracy of their determination leads to the inaccuracy of the calculation of the mechanical characteristic [1, 2]. In particular, the value of the inductive resistances of the windings is significantly influenced by the saturation of the magnetic
core, the change in the active resistances of the rotor winding due to displacement of current. Taking them into account in dynamic modes with the help of corresponding coefficients [2] does not guarantee the accuracy of the calculation results, especially for deep-slot motors.
The object of study is an IM with short-circuited rotor winding, which is powered by a three-phase network with a symmetric voltage system. For the analysis of electromagnetic processes in IM, we use a mathematical model, created using orthogonal coordinate axes, which allows to consider processes by computer simulation taking into account both saturation and displacement of current in bars of the short-circuit rotor windings with minimal computation. The magnetization characteristics of the main magnetic flux and the scattering fluxes are used to take into account saturation, and to take into account displacement of current the bars are separated by a height into n layers (2 < n < 5) which results in the n windings being covered by different magnetic scattering fluxes. The calculation algorithms are based on a mathematical model of the IM in the x, y axes, developed on the basis of the theory of imaging vectors [7], which allows to consider processes in the IM based on the theory of circuits.
Dynamics of motion of the rotor of the IM, operating in the mode of periodically variable loading, is described by the system of the DEs of electromechanical equilibrium which in the system of orthogonal axes x, y taking into account the division of each bar in height into n elementary ones, as well as subjecting the image vector of the supply voltage along x axis that is commonly practiced looks like
dVsx = . dt dySy
Rsisx ^ U m
dt
~m0Wsy = -moysx - Rsisy ;
dVix dt
dViy dt
d¥n
dm dt
dt dVny dt
' J
■ = (mo -mViy- Riiix ;
= -(®o - Riiiy ;
■ = (m0 -m)ny - Rninx ;
(i)
= -(m0 -m)nx - Rni
n'ny ■
Po
(y/sxisy sx )~Mc(t )!
where the indices sx, sy denote that the flux linkages ( y), currents (i), and active resistances (r) belong to the corresponding stator circuits; and ix,...,nx, iy,...,ny to the rotor ones; Um, m0 are the amplitude value and angular frequency of the stator winding phase voltage; m is the angular velocity of rotation of the rotor; J is the moment of inertia of the moving parts of the electric drive reduced to the shaft of the IM; p0 is the number of pole pairs.
The loading diagram of the mechanism must be known for modelling. Taking into account that the time dependence of the load torque is periodic, it is necessary to represent it in the form of a law of change, which corresponds to a complete cycle in the form Mc(t) = Mc(t + T), where T is the period.
An algorithm of calculation of characteristics. If the IM operates in one of the standard modes (full cycle is 10 minutes), then the transient is almost complete, and for a complete analysis of the motor operation it is enough to calculate the transient during the period. This can be done by integrating the DE system (1) using the numerical method [6].
DE system (1) includes 2 + 2n electrical equilibrium equations and one rotor dynamics equation. Therefore, when calculating the transient, it is necessary to rotate a matrix of the same order at each step (sub-step). In order to reduce the amount of calculations, we reduce the DE system (1), based on the following considerations.
The flux linkages of each IM circuit according to accepted assumptions consists of the sum
Wj _WjJ + Wcj of working flux linkage y/j which is nonlinearly dependent on the currents of all circuits, and scattering flux linkage which has a linear dependence, respectively, only on the stator currents or only on the rotor one. In addition, the flux linkage caused by the main work flow and the flux linkage of the slit scattering for all rotor circuits along the x axis are equal. The same applies to similar circuits along the y axis. The above makes it possible to divide the equilibrium equations of the DE system (1) into two parts by distinguishing a linear part in it. For this purpose it is necessary to replace the 5th equation by the difference of the 5th and the 3rd ones, the 6th equation by the difference of the 6th and 4th ones, etc. The first one is of the fourth order
dWsx y
- aoWsy
Write these two systems in the form
dt
dySy
Rslsx + U m
dt
d(lx dt d( dt
= -a0(sx - Rslsy ;
= (®o -®Viy- Rihx ;
= -(a0 -a)lx - Rllly
and is nonlinear, and the second of the 2(n-l) order is linear
d (lXdt(2 x )=(o- - ( y )- i'ix+r2'2 x ;
d (fly "(2y )
dt
■ = -(o -a)((lx -(2x )- rllly + r2l2y ;
d((lx -(nx ^
dt
:(0 -a)((
(ly -(ny)- rlllx + rnlnx ;
d (ly-(ny) ( V ) ---= -(®o - a Wlx -(nx ) - Il y + rnlny •
All Al
Al2
A22.
dii / dt din I dt
Bl h.
Determine the derivative of equation (2)
d--(All -Al2A2-2A2lI"1 (l
- Al2 A22lB2
in which only the elements of the matrices A11 and A12 depend on saturation. This allows to calculate the
elements of the matrices A221 and A21 once and use them
to determine at each step the integration the derivative
diii = A-l HT=A22
f
B2 - A
2l'
dii dt
Therefore, it is enough to rotate once the matrix of the 2(n-1) order and rotate the 4th order matrix at each integration step. The obtained formulas make it possible by numerical methods to reduce to Cauchy form the system (2) of the DEs of the electric equilibrium of the circuits.
Flux linkages of circuits are determined based on the use of magnetization curves by the main magnetic flux y and the scattering fluxes of the stator y and rotor windings
У/M_УM('M), Wcs _Wjs ('s ), Wcr _Wjr ('r ) ,
where
lM =
^sx + irx ) +iy + iry )
■ = ï~2-- • ■ = f
ls = y lsx + lsy ; lr = y '
■ = r 2 + - 2
lr = \ lrx + lry •
The currents of the rotor circuits are defined as the sum of the currents of n bar elements
irx _I j ; _Zj .
j_1 j_1 A mathematical model for calculating steady-state dynamic mode. In order to reduce the presentation of the computation of the steady-state dynamic calculation algorithm, we write the DE system (1) as a vector equation of the form
dx = (£T ?(y,x,3,f),
dt lsx I v '
where — =
dt
xy 0
is the matrix in which Lxy =
(3)
d( dl
is the complete matrix of differential inductances of the IM in coordinate axes x, y [7];
y =
(sx Um ' " 0 " lsx
(sy 0 0 lsy
(lx 0 0 llx
(ly ; 3 = 0 ; f = 0 ; x = lly
(nx 0 0 lnx
(ny 0 0 lny
a 0 Mc (( )_ a
x
In steady state mode of the electric drive system with periodic change of the load torque M(t) = M(t + T) flux linkages, currents, rotor speed, electromagnetic torque, etc. are varied by periodic laws. The task of calculating a periodic mode is to determine these dependencies. The solution of the system of equations (3) is the periodic dependencies of the components of the vector x(t) = x(t + T). The calculation by the stable method is inefficient for many reasons. In particular, CPU time is wasted, and if the process is becomes stable too slowly, then the oscillations at time t are little different from those for time t + T, so there is a problem of determining the time when the transient ends. Finally, the stable method is practically unsuitable for optimization calculations.
The most effective approach to calculating a steady periodic mode is to consider the problem as a boundary one [7], which allows to obtain periodic dependencies of coordinates in the timeless domain, that is, without resorting to the calculation of the transient. To do this, the system of continuous DEs (1) must be reduced to discrete ones, which are a point mapping of the dependencies of the coordinates during the process repetition period. There are many methods of algebraization in the literature that have both positive and negative sides: difference, collocation, including trigonometric, differential transformations, etc. The method, based on spline approximations of coordinates, developed in [8], makes it possible to formalize the algebraization process and is also numerically stable. It allows to obtain continuous dependencies of coordinates on a period on the basis of the obtained by calculation their discrete values in nodes of a mesh on a period. Note that the mesh of nodes can be taken uniformly. In the system of algebraic equations obtained by the approximation of variables, the values of the coordinates in m nodes of the period are unknown. As a result, taking into account periodic boundary conditions
Y t) = Y t+T), X t) = X t+T), we obtain the system of m* (2n+3) nonlinear algebraic equations, which can be represented as a vector equation
Y tX )= H _1Z (Y, X), (4)
in which H is the square matrix of size m(3+2n) of transition from continuous change of coordinates to their nodal values, whose elements are determined only by the
mesh step [8]; Y = (,..., ym), Z = (,..., ?m), X = ((,..., xm) are the vectors made up of values of vectors y , x, z in m nodes of the period.
By defining from vector (4) the vector X , it is possible to construct periodic dependencies of all coordinates, including electromagnetic of torque, power, etc.
Direct application of the iterative method to the solution of system (4) is practically impossible due to the divergence of the iterative process. A reliable method of solving the problem is the method of continuation by parameter [9]. However, in the system of nonlinear
algebraic equations there are two disturbing actions: applied voltage - vector U = ((,...,um)and vector of nodal values of load torque - F = (/1,.,fm). It is impossible to increase them at the same time, so the problem is solved in two stages, the essence of which is to increase them alternately in proportion to a certain parameter. First, we increase the applied voltage, and then, taking it unchanged, we increase the nodal values of the applied torque. This makes it possible to determine the time dependencies of the coordinates in the steady-state periodic mode of operation of the IM at a given law of change of the applied torque.
The steady-state mode calculation algorithm is the basis for the calculation of static characteristics, which can be obtained as a sequence of steady-state modes calculated with a set of coordinate values, which is taken as an independent variable, which can be any value: moment of inertia, pulse density of the load torque; ratio between pulse duration and pause, pulse rate, maximum and minimum torque, period duration, etc. In addition, under cyclic loading mechanical resonance is possible, which can be detected by mathematical modelling.
The problem of calculating static characteristics can be solved by a differential method, the essence of which is the differentiation of algebraic equation (4) on an independent variable, for example s, as a parameter. As a result of differentiation, we obtain a nonlinear system of DE in the form
^ = (5)
ds ds
The static multidimensional characteristic as a dependence of periodic curves on the independent variable s is obtained by integrating system (5) with parameter s. The initial conditions should be those obtained as a result of the first stage of the calculation at a given supply voltage. At each integration step, the result can be refined by the Newton method. During integration as well as iterative refinement, it is necessary to determine the differential inductance of circuits as nonlinear functions of currents.
Results of investigations. Below are examples of calculation results performed using the above algorithms on the example of the IM with short-circuited rotor 4AP160S4Y3 (P = 15 kW, U = 220 V, I = 29.9 A, p0 = 2).
Figure 1 shows the time dependences of the relative values of the electromagnetic torque (Fig. 1,a) and the current value (Fig. 1,b) in the transient during start up of the IM with cyclic loading, at which the torque of load varies with the period T = 0.16 s in the range from idling to nominal value, the moment of inertia J = 0.5 kg-m2, and the density is 60 %, and Fig. 2 presents the same dependencies, but with less moment of inertia J = 0.1 kg-m2.
1.6 "'"0
Fig. i
0.4
a
0.4
b
Fig. 2
Figures 3, 4 present an example of calculated by the presented in this work which correspond dependencies in method of solving the boundary value problem of periodic steady state shown in Fig. 2. curves of current, electromagnetic torque and load,
2.5 1.8 1.1 0.4 0.3 _ 1
ht- ht-K'K
m; j
A
j \ /
j \\M' /
j \ /
j /
1 V /
2
i.6
1.2
0.8
0.4 t 0
A
f \ i j \
1 \
1 \ 11 K \
/ \ A ] \y\
i
0
772
3T/2
2T
772
3T/2
t 2 T
b
Fig. 3. Periodic dependencies (two periods shown) of the relative values of the load torque (Mc ), electromagnetic torque (Me ) and current (I*) calculated at the moment of inertia J = 0.i kg-m2 by the method of solving the boundary value problem
a
Fig. 4. Dependence of the electromagnetic torque of the motor on the relative value of the duration of the period of change of load (at the point TJTe = 0.16 there is a mechanical resonance)
Conclusions.
1. The developed calculation methods and corresponding algorithms make it possible to use mathematical modelling to analyze the operation of induction motors with a short-circuited rotor taking into account the saturation and displacement of currents in the rotor bars under different laws of change of periodic load.
2. The algorithm of calculation of the steady-state periodic modes at cyclic loading allows to obtain the periodic dependencies of coordinates in the timeless domain, which ensures high speed.
3. The mathematical models developed can be used to design and analyze the operation of electric drives with periodic load.
REFERENCES
1. Voldek A.I., Popov V.V. Elektricheskiye mashiny. Mashiny peremennogo toka [Electric machines. AC machines]. Saint Petersburg, Piter Publ., 2010. 350 p. (Rus).
2. Safaryan V.S., Gevorgyan S.G. Ascertainment of the equivalent circuit parameters of the asynchronous machine. Energetika. Proceedings of CIS higher education institutions and power engineering associations, 2015, no. 6, pp. 20-34. (Rus).
3. Rogal V.V., Kapshtik V.S. Reactive power compensation in intermittent duties. Electronics and Communication. Thematic issue «Electronics and Nanotechnology», 2011, no. 3, pp. 101108. (Ukr).
4. Petrushin V.S., Plotkin J.R., Yenoktaiev R.N., Bendahmane Boukhalfa. Development of energy-efficient asynchronous electric drive for intermittent operation. Bulletin of the National Technical University «KhPI». Series: Problems of automated electrodrive. Theory and practice, 2019, no. 16 (1341), pp. 7079. doi: 10.20998/2079-8024.2019.16.13.
5. Petuhov S.V., Krishyanis M.V. Elektroprivod promyishlennyih ustanovok [Electric driver industrial-scale plants]. Arkhangelsk, S(A)FU Publ., 2015. 303 p. (Rus).
6. Hrisanov V.I. Transient process analysis at various methods of starting asynchronous machines. Technical electrodynamics. Thematic issue «Electric drive», 2000, pp. 24-27. (Rus). ZFil'ts R.V. Matematicheskie osnovy teorii elektromekhanicheskikh preobrazovatelei [Mathematical foundations of the theory of electromechanical transducers]. Kyiv, Naukova dumka Publ., 1979. 208 p. (Rus).
8. Malyar V.S., Malyar A.V. Mathematical modeling of periodic modes of operation of electrical devices. Electronic Modeling, 2005, vol.27, no.3, pp. 39-53. (Rus).
9. Zhulin S.S. The method of continuation by parameter and its application to the tasks of optimal control. Numerical methods and programming, 2007, vol. 8, no. 2, pp. 205-217. (Rus).
Received 11.02.2020
V.S. Malyar', Doctor of Technical Science, Professor, O.Ye. Hamola1, Candidate of Technical Science, Associate Professor,
V.S. Maday1, Candidate of Technical Science, Associate Professor,
1 Lviv Polytechnic National University, 12, S. Bandera Str., Lviv, 79013, Ukraine, e-mail: [email protected], [email protected], [email protected]
How to cite this article:
Malyar V.S., Hamola O.Ye., Maday V.S. Modelling of dynamic modes of an induction electric drive at periodic load.
Electrical engineering & electromechanics, 2020, no. 3, pp. 9-14. doi: 10.20998/2074-272X.2020.3.02.
