CC BY
13
-□ □-
Проведено теоретичнг дослгдження електромагттних про-цесгв при випробуваннях силових трансформаторе в режимг дослгдного короткого замикання на основi чисельног реатзаци тривимгрног моделг магнгтного поля у частотному формулю-ваннг. Шляхом верифгкацп данихрозрахункгв частотних i змгн-них у часг моделей магттного поля обгрунтовано достовiрнiсть i точнгсть визначення параметрiв дослгдного короткого замикання силового трансформатора у частотному формулюван-нг. Визначено основнг закономгрностг розподглу магттного поля в об'емг активног частини трансформатора. В зонах локалгза-цг магнгтного поля 30 розподгл напруженостг е рiвномiрним i визначаеться 20розподглом у горизонтальному перерЫ актив-ног частини на середин висоти фазних обмоток. Значення осьо-вог компоненты напруженостг магнгтного поля наближають-ся до 96-97 % вiд модуля вектора напруженостг. Реализовано ефективний пiдхiд до польового моделювання на основi деком-позицг розрахунковог областг на просторовi зони. Кожнгй роз-рахунковгй зонi поставлено у вiдповiднiсть електричний контур схеми замщення. Розподгл електричних потенцгалгв у горизон-тальних перергзах провiдникiв мгж котушками або мгж витками обмоток прийнято рiвномiрним. Суперпозицию магнгтних полгв у просторових зонах здгйснено засобами динамгчного синтезу за критергями мгнгмальног струмовог похибки для електричних контурiв схеми замщення. Декомпозиция 30 областг польового моделювання на центральну i торцевi зони здгйсню-еться на вгдстанг 10-15 % висоти фазних обмоток, що забез-печуе високу точтсть розрахунку напруженостг магнгтного поля гз похибкою не больше 1,62 %. Витрати часу для польового моделювання електромагттних процесгв у режимг дослгдно-го короткого замикання зменшено у 5 раз, а вимоги до потуж-ностг обчислювальних апаратних ресурсгв знижено у 4 рази. Високу точтсть гдентифжацп параметргв дослгдного короткого замикання трифазних трансформаторiв пгдтвердже-но поргвнянням даних розрахункгв гз результатами випробу-вань в умовах приватного тдприемства «Елтз» (Запорiжжя, Украгна). Похибки розрахункгв не перевищують 1,42 % для активних втрат i 1,39 % для напруги короткого замикання. Запропонований шдхд гз використанням методiв декомпозицп та динамгчного синтезу дозволяе значно шдвищити ефектив-нгсть попереднього етапу конструкторськог пгдготовки вироб-ництва г може бути використаний при реалгзацгг задач оптимг-зацгг конструктивних ргшень
Ключовг слова: електромагнгтне поле, трифазний трансформатор, дослгдне коротке замикання, декомпозицгя, дина-мгчний синтез
-□ □-
UDC 621.314.21
[DOI: 10.15587/1729-4061.2018.1402361
ENHANCING THE EFFECTIVENESS OF CALCULATION OF PARAMETERS FOR SHORT CIRCUIT OF THREE-PHASE TRANSFORMERS USING FIELD SIMULATION METHODS
D. Yarym bash
Doctor of Technical Sciences, Associate Professor* E-mail: yarymbash@gmail.com S. Yarym bash PhD, Associate Professor* E-mail: kstj06@gmail.com M . K o t s u r PhD, Associate Professor Department of electrical and electronic apparatuses** E-mail: kotsur8@gmail.com T. D i vc h u k Senior Lecturer Department of electrical machines** E-mail: DIV2009@i.ua *Department of electrical machines** **Zaporizhzhya National Technical University Zhukovskoho str., 64, Zaporizhzhya, Ukraine, 69063
1. Introduction
Power transformers are some of the most common electrical devices that play a major role in the transportation of electricity from the place of its generation to its users. The established capacity of power transformers can by order of magnitude exceed the capacity of electricity generation, the transformation multiplicity of which can be equal to 7-8 [1, 2]. That is why there are high requirements to technical and economic indicators and reliability of power transformers that must be provided even at the stage of design preparation of production [3].
Under the steady operating modes of power systems and distribution networks, the main effect of transformer equipment is determined by the parameters of the replacement scheme for the test short circuit (SC) [3-5]. Resistances of SC are calculated according by specifications of the transformer for voltage and short circuit losses [1, 2].
Engineering design techniques based on the methods for researching linear systems fragmentally take into consideration design features of the active part of power transformer and its other structural elements [1, 2]. This can lead to a distortion of real processes of electromagnetic conversion of AC energy in power transformers [4, 5]. Therefore, the accuracy
©
of the design calculations of specifications of short circuit of the transformer is low, and permissible errors can approximate 10-15 % [3]. This causes the relevance of development of high-precision and high-performance methods for calculation of the transformer parameters in the mode of the studied SC.
_2. Literature review and problem statement_
To determine the parameters and characteristics of power transformers, the schematic models based on the theory of electric and magnetic circuits are normally used [6, 7]. Comparative simplicity of the algorithms of computer implementation caused the use of schematic models in engineering calculations at the stage of design preparation of production [4, 5]. However, a number of assumptions and simplifications, which laid the basis of schematic models, lead to deviations from the actual electromagnetic processes of energy conversion of AC during its transformations [8, 9]. Simplified descriptions of surface effects, nonlinear properties of active and structural materials [10], disregarding the peculiarities of a multi-component spatial structure of the active part of transformers [11] can lead to significant calculation errors [12-14]. It imposes a limitation to the application of these methods for calculation of currents, voltage and electrical resistance under the dynamic conditions of short circuit with the shock electro-dynamic efforts and thermal loads [15, 16].
It is possible to take into consideration the influence of structural factors on parameters of a short circuit of power transformers more completely, using field simulation [17-19]. 3D field simulation can ensure high accuracy in determining the parameters of power transformers. However, its computer implementation, based on the methods of finite elements in multi-component areas of the active part with non-linear electro-physical properties, is complicated by great consumption of time and requirements for computational resources [17]. Paper [20, 21] proposed the 3D model of electromagnetic processes to determine electrical and electro-power parameters of electro-technical systems of the AC transformation. However, its application is limited only to systems that have linear magnetic properties of materials. In papers [22-24], it is proposed to use the differentiation of dimensions of finite elements and approximation by Lagrange first order polynomials for effective numerous implementations of the field models. The application of this approach to power transformers is complicated by the significant nonlinearity of magnetic properties and the difference of geometric dimensions of separate elements [25-27]. Some authors try to introduce simplifications into the geometric model [28, 29], neglect nonlinearity of the magnetic properties of materials [24] or to apply 2D models [30, 31]. These assumptions and simplifications decrease the accuracy of results of simulation and in some cases can distort a real picture of electromagnetic processes. The field calculation model of the power transformer in the mode of the studied short circuit should in detail display the multi-component structure of the active part and take into consideration the nonlinearity of electro-physical and magnetic properties of the active materials. The effectiveness of the numerical implementation of the field model must meet the requirements of the problems of design parameters optimization. This causes the need for development of specific approaches that improve the computation accuracy and efficiency and ensure adap a-tion of 3D-field simulation to the tasks of optimum design of power transformers.
_3. The aim and objectives of the study_
The aim of this study is to develop an effective approach to field simulation of the electromagnetic processes of the three-phase transformer through the decomposition of the 3D area of the active parts into sub -areas for the synthesis of parameters of the corresponding electric circuits in SC parameters with high precision.
To accomplish the aim, the following tasks have been set:
- to determine mathematical statement of the AC field model that reflects the features of electromagnetic energy conversion in the mode of the studied SC of the power transformer with the multicomponent active part of a complex structure;
- to determine parameters of decomposition of the calculation region of the active parts of the transformer in the spatial zones of the reduced volume in order to reduce the number of finite elements, to ensure the accuracy and effectiveness of the computation processes;
- to develop the method of synthesis of field zones into the electric circuit by the criterion of current error minimization for the refined determining of SC parameters.
4. Statement of the mathematical model of electromagnetic processes in the transformer for the mode of the studied short circuit
The 3D area of the field simulation of electromagnetic processes in the power transformer includes the active part, which is placed in the tank - 6 (Fig. 1, a). The active part consists of the systems of windings, the magnetic system, the systems of longitudinal and main electric insulation. The system of windings consists of the main and adjusted windings on the high voltage side (HV) - 1 and the winding on the low voltage side (LV) - 2 with coils, turns, layers, etc. [33]. In the charged magnetic system of the transformer - 3, one distinguishes limbs and yokes that are brought together by yoke beams. The system of electric insulation of the transformer consists of solid insulation from electric cardboard - 4, paper electric insulation and transformation oil - 5 (Fig. 1, a).
Fig. 1. Three-phase power transformer: a — calculation; b, c, d — upper, central and lower decomposition zones
b
a
According to data from [17], electro-physical processes in dielectric materials can be considered similar to the processes of electrical conductivity in materials with relatively low electrical conductivity. That is why electromagnetic fields in the calculation region fl (Fig. 1, a) can be described by Maxwell equations in frequency statements [20, 21]:
vxi(msm ) vxAj
=-o, K )( ,»A j+vv, );
v[°j (e, \vv, ] = (e, )A, )>
with potentiality conditions [20]: {E, = -grad(V); B, = rot (A, ), relationships of form [20]:
{H j =[mom (| ^j | )]-1 »j ; J =°j (e )E j
(1)
(2)
(3)
where A, V are the vector and magnetic potentials, B is the induction of the magnetic field, E, H are the voltages of electric and magnetic fields, J is the density of current, a(0) is the temperature dependence of specific electric conductivity, 0 is the temperature, |o is the magnetic permeability of vacuum, |(|B|) is the equivalent relative magnetic permeability, ro is the angular frequency.
The system of the field equations (1) is supplemented by conditions of combination of magnetic and electric fields on the boundaries of the cells with different electro-physical properties [21]:
{n x(H+-H-) = 0; n ■( J+-J-) = 0,
(4)
H=1 J№ j-B *)/m>m, (IB ,1 )
v i
pl,=JJJo,rl ( J,-'J', )■ AxdyAz;
J, = -o,-(grad (v ) + iraAl),
(5)
as well as currents in the phase windings of HV and LL.
Phase high-voltage windings are attached to the three-phase power supply source by scheme Y, and low-voltage windings, connected by scheme A, are shunted by branches with small active resistance. The restrictions for the power source of test voltage are taken into account using linear resistances between this source and the windings of the
transformer. This will not lead to any significant increase in current error provided that these resistances are by several orders of magnitude lower than resistance of SC, determined by specifications. Under the mode of the studied SC, currents in the windings are equal to the rated values, in accordance with the phase voltage at the HV side {DU k }, k = 1,2, 3:
DUk = -,»4'k + Ik ■ Rk, k = 1,2,3; 0 + îk ■ Rk, k = 4,5,6.
(6)
and homogeneous boundary conditions at the external surface of the tank [26] or impedance conditions [32] at its internal side.
For calculation region fl, the temperature mode of the power transformer is accepted as stationary. The temperatures of the elements of the active parts are determined by the project or by the data, which are admissible for heating resistance A [1]. The anisotropy of magnetic properties of cold rolled electro-technical steel was taken into account according to [25].
For the elements of the power transformer in the 3D field fl, densities of currents, active losses (for electrically conductive materials) and energy of the magnetic field are calculated [20]:
Computer implementation of the field model (1) to (4) is carried out in the structure of the specialized software of field simulation COMSOL Multiphysics [33]. 3D area is divided into three-dimensional finite elements - tetrahedrons with edges, which are approximated by Lagrange second order polynomials. To increase effectiveness of the numerical implementation of the model, differentiation of dimensions of finite elements is performed. In sub-regions with distinctly pronounced nonlinearity of magnetic properties and localization of magnetic field, dimensions of finite element decrease and increase near the outer boundary of the calculation region.
5. Results of the study of magnetic fields of the transformer under the mode of test short circuit. Statement of criteria for the decomposition of active part and the synthesis of electric circuits of field zones
Field simulation of electromagnetic fields in the area of the active part was carried out for the three-phase power transformer of voltage class of 35 kV and power of 1,600 kWA. The windings are connected by Y/Y0-0 circuit. Rated values of linear voltages on the side of HV and LV were 20 kV and 0.4 kV. For the mode of the test CS, the amplitudes of phase voltages were determined by specified value of SC voltage, which amounted to 6.5 % of the rated phase voltage from the HV side of the transformer. The symmetric three-phase system of electrical voltages with the phase shift vector [0; 2 n/3;
4 n/3] for frequency of 50 Hz of AC was used.
The results of 3D simulation are represented by distributions of voltage of the magnetic field and its axial components in vertical and horizontal cross sections of the active part of the power transformer (Fig. 2, a).
According to field calculations, it was determined that for the mode of the test SC the localization of the magnetic field occurs in the phase windings systems. In the determined localization areas, directions and module of vectors of magnetic field intensity almost never change. The biggest component is the projection of Hz, the value of which is 96-97 % of the module of magnetic intensity vector. This allows us to consider that the 3D distribution of the magnetic field intensity in the area of its localization is determined by the 2D distribution in the horizontal plane of cross section of the active part. The change in electrical potentials on the surface of horizontal cross-sections of the conductors between the coils or between the turns of the windings can be neglected and we can accept the assumption of the uniform distribution of electrical potentials.
The features of the magnetic fields distribution allow applying for calculation region fl the principle of superposition of electromagnetic fields, calculated separately in the central (Fig. 1, c) and end zones (Fig. 1, b, d).
iuuneu u ¡usics
To reduce the number of finite elements and dimensionality of equations of the finite element model, it is advisable to apply the distribution of the active part structure into the subsystems with separate operating elements or components [33]. This will significantly reduce the number of functional relations between the parameters of each subsystem and the number of independently varied parameters to minimize current errors for phase currents of the studied SC.
a b
Fig. 2. Magnetic field of the active part of the power transformer: a — amplitudes of magnetic field intensity; b — amplitudes of the projection of magnetic field intensity onto the vertical axis of the transformer
Decomposition of the calculation region Q implies the division into several calculation areas, each of which can be put in conformity with the circuits of the replacement scheme. These areas are separated from each other by horizontal cross sections that are perpendicular to axis 0z and pass through the axial insulating channels between the coils or between the turns of the windings. The distance of horizontal cross-sections from the upper and lower ends is 10-15 % of the height of the windings.
The central area includes the central part of phase windings of HV, LV and limbs, main insulation, transformer oil and the central part of the tank (Fig. 1, c). The end zones include the upper or the lower part of the tank, the upper or the lower yokes, relevant parts of HV, LV windings and limbs of different phases, solid insulation and transformer oil (Fig. 1, b, d). Division of the field of simulation Q into the zones ensures a reduction of time of implementation of the field models due to the fewer number of finite elements in the zones of smaller volume.
It is possible to display the volume of the region of electromagnetic energy conversion Q by the sum of the volumes of zones and to determine electrical losses and magnetic energy by the sum of values in the respective areas.
For the synthesis of electric circuits on the boundaries of connection of the areas of phase windings, we accept the conditions of the equality of electric potentials:
where {AU} is the voltage drop in the sections of phase windings; i is the number of sequential branches in the circuit of the phase winding, indices k e(1,2,3) belong to the HV windings of A, B, C phases, indices k e(4,5,6) belong to LV windings of A, B, C phases, respectively.
In conditions (7), each value of const is "free". That is why
a voltage drop in the section of phase windings {AU^ }k
cannot be assigned in advance and should be determined from the solution of the problem of current errors minimization. For the synthesis of electric circuits and superposition of magnetic fields, in the calculation region, it is possible to state objective function with the components that are determined by the sum of the squares of current errors for currents in the sections of windings:
min (M2 (D,R)) = I min (M2 (D,R), ) =
=I
min
I i * - i *
(9)
In the general case, component of the vector of objective function M2(D,R)k is the function of vectors of independently varying D and dependent R parameters for each winding of the transformer. A voltage drop on the section of the circuit of the winding is considered as projections of the vector, irrespective varying parameters D. This makes it possible to convert the components of objective function M2(D,R)k into components of function M 2(D)k.
To implement conditions (9), it is advisable to apply the method of dynamic programming [34]. The vector of optimization parameters with components
D = {AU.}
I i ie(1,3),ke(1,6)
can be represented as the sum of vectors of initial approximation
D,
= {DU »..Il k
lke(1.6)
and its specification
8D={5AÛ M U,,«)
That is why optimality conditions (9) can be put in conformity with the system of recurrent equations of Bellman [9] for each winding of the transformer:
min AM2 =o(D0.k ) =
= mm... min
5D1 5D,
[(AM i2 (Do, , 5DU ) +... AM 2 (D m_i,, 5D M , ))]
.(10)
ke(1.6)
V = Vm|r = constj,i,AUhJ = Vi. -Vm|r ,
[_ 11 j ,i+1 11 j ,i+1
for which system of equations (6) will take the form:
(7)
AÛ =I(AÛ y }, = I(MJ.},
. ,=1.2.3 .
AÛk, = -jœ*,. + i,. • R,..k e(1,6).
¿=4.5.6
(8)
Relations (9). (10) determine the conditions for dynamic connections of "free" parameters in equations (7) for the central and end 3D zones of region Q. The iteration calculation process of the method of dynamic programming is implemented in the assigned limits for current e/<0.1 % and angular errors e9<0.5 %.
For the mode of the studied SC. we applied the conditions that the sum of instantaneous values of voltage drop in the circuits of the parts of phase windings of each decomposition zone must be equal to instantaneous values of SC voltage for
£
each phase. Spatial distribution of magnetic field intensity in sub-areas of the upper - I, central - II and lower - III zone of the decomposition was compared with the distribution of intensity in the 3D area of the active part (Fig. 3). For the decomposition method, the largest deviations of calculations of magnetic field intensity are characteristic of horizontal surfaces of the areas that cross windings (Fig. 3). However, the maximum value did not exceed 1.62 %.
Fig. 3. Distribution of normal component of the magnetic field intensity in the calculation area and zones
connection of electric circuits of phase windings sections ensures effective implementation of the 3D modeling of magnetic fields in the regions of a complex spatial form. High reliability and accuracy of the data of modeling and determining parameters of a short circuit are achieved thanks to taking into account of finite dimensions of the magnetic system, the winding system and features of reciprocal location. Time consumption decreases by 5 times, and requirements for capacity of computation resources for the refined calculation of the SC parameters decrease by 4 times.
The application area of the results of the conducted research is limited to the most common types of three-phase power transformers of voltage classes of up to 35 kV including the circuits of connection of windings Y/Yo-0, Y/D-11, D/Yo-11 [1]. For the transformers of voltage class of 110 kV and above, the linear input inside HV windings is applied, which is performed by an electric circuit with two parallel branches [1]. In addition, high voltage classes require a substantial increase in dimensions of axial insulation gaps in the middle of LV windings, which leads to a multiple increase in the number of decomposition zones for calculation of the SC of the transformer. That is why effective approaches to high-precision determining the SC parameters for transformers with voltage class of 110 kV and above require further development.
6. Discussion of results of simulation of electromagnetic processes in the mode of the test short circuit of the three-phase transformer
Table 1
Errors of the method of dynamic synthesis of field models for the criterion of minimum current error
Calculation regions (zones)
Magnetic field energy
Active power losses
Voltage drop (module/phase)
Verification of SC parameters by the data of field simulation of the parameters of SC of the transformer was performed for SC voltage, the value of which was 740.2 V (6.41 %) for the frequency-domain model and 741.7 V (6.42 %) for the time-dependent model. Corresponding calculation values of SC losses were equal to 17.745 kW and 17.815 kW.
Thus, relative discrepancies in frequency calculations did not exceed 0.4 % for losses and 0.2 % for voltage of short circuit. The relative error in calculation of losses and short-circuit voltage did not exceed 1.5 %.
Errors of calculation of active power losses, magnetic field energy and voltage drops using the method of decomposition of the 3D area of the active part of the power transformer are shown in Table 1. The relative errors in the decomposition areas do not exceed: 0.47 % for magnetic field energy, 0.53 % for electric losses, 0.51 % for the voltage drop module, 0.29 % for the voltage drop phase. The relative errors of calculation of SC parameters make up 1.42 % for losses and 1.39 % for SC voltage. Validation of calculation parameters of the short circuit was carried out by comparing the calculation data using technique [26] with the results of tests, obtained at private enterprise "Eltiz" (Zaporizhye, Ukraine).
Reliability and accuracy of the results of simulation using the methods of decomposition and of dynamic synthesis of the windings sections circuits are ensured provided that the distribution of electrical potential in the connection boundaries approaches the uniform one. In this case, the sum of active losses and magnetic field energy in the calculation regions will correspond with high accuracy to active losses and magnetic energy for the 3D area of the active part.
The proposed combination of the methods of decomposition and of dynamic synthesis with "free" conditions of
Upper zone
0.47
0.53 %
0.51 %/0.29 %
Central zone
0.19
0.31 %
0.38 %/0.18 %
Lower zone
0.25 %
0.48 %
0.46 %/0.28 %
Area of active part
0.76 %
1.42 %
1.39 %/0.78 %
In the future, the proposed methods of decomposition and of dynamic synthesis can be adapted for studying high-voltage reactors and rotating electric machinery and determining their parameters in test modes.
7. Conclusions
1. It was found that mathematical statements for the model of the magnetic field reflect the features of the electromagnetic energy conversion in the mode of the test SC of the power transformer with the multicomponent active part of a complex structure. Using the verification of the data of calculation of frequency-domain and time dependent models of the magnetic field, we substantiated validity and accuracy of results when using frequency statements for determining the parameters of the studied SC of the power transformer. The relative discrepancies of frequency calculations did not exceed 0.4 % for losses and 0.2 % for short circuit voltage. Comparison of calculation data with test results from the private enterprise "Eltiz" (Zaporizhye, Ukraine) is proved by high accuracy of identification of the parameters of the test SC. Relative error of calculation of losses and voltage of short circuit does not exceed 1.5 %.
2. It was established that in the test SC, the magnetic field is located in phase windings. In the localization areas,
directions and the module of the vector of magnetic field almost never change. The axial component, the value of which is 96-97 % of the module of the magnetic tension vector, is the largest. This makes it possible to determine the 3D distribution of the magnetic field intensity in the localization area according to the data of 2D distribution in the horizontal plane of the cross section of the active part in the middle of the phase windings height. The change in electrical potential on the surface of horizontal cross-sections of the conductors between the coils or between the turns of the windings can be neglected and the assumption of uniform distribution can be accepted. The calculation area is divided into the central and end zones. Each of the zones was set in compliance with the electrical circuit in the substitution scheme. Decomposition of the 3D area of field simulation
into the zones at the distance of 10...15 % of the height of the phase windings from their ends ensures high precision of calculation of the spatial distribution of the magnetic field intensity with the error of not more than 1.62 %.
3. It was found that the dynamic synthesis of field areas with electric circuits in the substitution scheme by the criterion of a minimum current error significantly increases computational efficiency of determining SC parameters. In the set restrictions for current s/<0.1 % and angular errors s<p<0.5 %, relative calculation errors do not exceed 1.42 % for active losses and 1.39 % for SC voltage. Application of the methods of decomposition and of dynamic synthesis of electric circuits make it possible to decrease time consumption by 5 times and decrease the requirements for computing resources by 4 times.
References
1. Biki M. A. Proektirovanie silovyh transformatorov. Raschet osnovnyh parametrov. Moscow: Znak, 2013. 612 p.
2. Khaparde S., Kulkarni S. Transformer Engineering: Design and Practice. CRC Press, 2004. 496 p. doi: https://doi.org/ 10.1201/9780203970591
3. C57.12.90-2006 - IEEE Standard Test Code for Liquid-Immersed Distribution, Power, and Regulating Transformers. doi: https:// doi.org/10.1109/ieeestd.2006.320496
4. Lul'e A. I. Process vklyucheniya silovogo transofrmatora na holostoy hod i korotkoe zamykanie // Elektrotekhnika. 2008. Issue 2. P. 2-18.
5. Leytes L. V. Elektromagnitnye raschety trasnfomatorov i reaktorov. Moscow: Energiya, 1981. 365 p.
6. Novash I. V., Rumiantsev Y. V. Three-phase transformer parameters calculation considering the core saturation for the matlab-sim-ulink transformer model // ENERGETIKA. Proceedings of CIS higher education institutions and power engineering associations. 2015. Issue 1. P. 12-24.
7. Schiop A., Popescu V. Pspice simulation of power electronics circuit and induction motor drives // Revue Roumaine des Sciences Techniques-Serie Electrotechnique et Energetique. 2007. Vol. 52, Issue 1. P. 33-42.
8. Speed synchronization methods of the energy-efficient electric drive system for induction motors / Kotsur M., Yarymbash D., Kotsur I., Bezverkhnia Y., Andrienko D. // 2018 14th International Conference on Advanced Trends in Radioelecrtronics, Telecommunications and Computer Engineering (TCSET). 2018. doi: https://doi.org/10.1109/tcset.2018.8336208
9. Jamali M., Mirzaie M., Asghar-Gholamian S. Calculation and Analysis of Transformer Inrush Current Based on Parameters of Transformer and Operating Conditions // Electronics and Electrical Engineering. 2011. Vol. 109, Issue 3. doi: https://doi.org/10.5755/ j01.eee.109.3.162
10. Singh A. K., Patel S. Mitigation of Inrush Current For Single Phase Transformer by Control Switching Method // International Journal of Electronics, Electrical and Computational System. 2015. Vol. 4. P. 146-150.
11. Taghikhani M. A., Sheikholeslami A., Taghikhani Z. Harmonic Modeling of Inrush Current in Core Type Power Transformers Using Hartley Transform // IJEEE. 2015. Vol. 11, Issue 2. P. 174-183 doi: https://doi.org/10.22068/IJEEE.1U.174
12. Chiesa N., Mork B. A., H idalen H. K. Transformer Model for Inrush Current Calculations: Simulations, Measurements and Sensitivity Analysis // IEEE Transactions on Power Delivery. 2010. Vol. 25, Issue 4. P. 2599-2608. doi: https://doi.org/10.1109/ tpwrd.2010.2045518
13. Khederzadeh M. Mitigation of the impact of transformer inrush current on voltage sag by TCSC // Electric Power Systems Research. 2010. Vol. 80, Issue 9. P. 1049-1055. doi: https://doi.org/10.1016/j.epsr.2010.01.011
14. Tykhovod S. M. Transients modeling in transformers on the basis of magnetoelectric equivalent circuits // Electrical Engineering and Power Engineering. 2014. Issue 2. P. 59-68. doi: https://doi.org/10.15588/1607-6761-2014-2-8
15. Increasing of thermal reliability of a regulated induction motor in non-standard cycle time conditions / Kotsur M., Kotsur I., Bezverkhnia Y., Andrienko D. // 2017 International Conference on Modern Electrical and Energy Systems (MEES). 2017. doi: https://doi.org/10.1109/mees.2017.8248960
16. Lazarev N. S., Shul'ga R. N., Shul'ga A. R. Toki vklyucheniya silovyh transformatorov // Elektrotekhnika. 2010. Issue 11. P. 11-17.
17. Podol'cev A. D., Kontorovich L. N. Chislenniy raschet elektricheskih tokov, magnitnogo polya i elektrodinamicheskih sil v silovom transformatore v avariynyh rezhimah s ispol'zovaniem MATLAB/SIMULINK i COMSOL // Tekhnichna elektrodynamika. 2011. Issue 6. P. 3-10.
18. Majumder R., Ghosh S., Mukherjee R. Transient Analysis of Single Phase Transformer Using State Model // International Journal of Innovative Research in Science, Engineering and Technology. 2016. Vol. 5, Issue 3. P. 3300-3306.
19. Computer simulation of electromagnetic field with application the frequency adaptation method / Yarymbash D. S., Yarym-bash S. T., Kotsur M. I., Litvinov D. O. // Radio Electronics, Computer Science, Control. 2018. Issue 1. P. 65-74. doi: https:// doi.org/10.15588/1607-3274-2018-1-8
20. Yarymbash D. S., Oleinikov A. M. On specific features of modeling electromagnetic field in the connection area of side busbar packages to graphitization furnace current leads // Russian Electrical Engineering. 2015. Vol. 86, Issue 2. P. 86-92. doi: https:// doi.org/10.3103/s1068371215020121
21. Yarymbash D. S. The research of electromagnetic and thermoelectric processes in the AC and DC graphitization furnaces // Nau-kovyi Visnyk Natsionalnoho Hirnychoho Universytetu. 2015. Issue 3. P. 95-102.
22. A New Approach of the Induction Motor Parameters Determination in Short-Circuit Mode by 3D Electromagnetic Field Simulation / Kotsur M., Yarymbash D., Yarymbash S., Kotsur I. // 2017 IEEE International Young Scientists Forum on Applied Physics and Engineering (YSF). 2017. P. 207-210. doi: https://doi.org/10.1109/ysf.2017.8126620
23. A New Simulation Approach of the Electromagnetic Fields in Electrical Machines / Yarymbash D., Kotsur M., Subbotin S., Oliinyk A. // 2017 International Conference on Information and Digital Technologies (IDT). 2017. P. 452-457. doi: https:// doi.org/10.1109/dt.2017.8024332
24. Milykh V. I., Polyakova N. V. An analysis of harmonic composition the AC magnetic field associated with a rotating rotor turbine generator, at idle speed and short circuit modes // Electrical Engineering And Power Engineering. 2013. Issue 2. P. 5-13.
25. Features of Defining Three-Phase Transformer No-Load Parameters by 3D Modeling Methods / Yarymbash D., Yarymbash S., Kylymnyk I., Divchuk T., Litvinov D. // 2017 International Conference on Modern Electrical and Energy Systems (MEES). 2017. P. 132-135. doi: https://doi.org/10.1109/mees.2017.8248870
26. An Application of Scheme and Field Models for Simulation of Electromagnetic Processes of Power Transformers / Yarymbash D., Kotsur M., Yarymbash S., Kylymnyk I., Divchuk T. // 2018 14th International Conference on Advanced Trends in Radioelecrtron-ics, Telecommunications and Computer Engineering (TCSET). 2018. P. 308-313. doi: https://doi.org/10.1109/tcset.2018.8336209
27. Dual Reversible Transformer Model for the Calculation of Low-Frequency Transients / Jazebi S., de Leon F., Farazmand A., Deswal D. // IEEE Transactions on Power Delivery. 2013. Vol. 28, Issue 4. P. 2509-2517. doi: https://doi.org/10.1109/tpwrd. 2013.2268857
28. Tang Q., Guo S., Wang Z. Magnetic flux distribution in power transformer core with mitred joints // Journal of Applied Physics. 2015. Vol. 117, Issue 17. P. 17D522. doi: https://doi.org/10.1063/1.4919119
29. Cundeva S. A transformer model based on the Jiles-Atherton theory of ferromagnetic hysteresis // Serbian Journal of Electrical Engineering. 2008. Vol. 5, Issue 1. P. 21-30. doi: https://doi.org/10.2298/sjee0801021c
30. Rashtchi V., Rahimpour E., Fotoohabadi H. Parameter identification of transformer detailed model based on chaos optimisation algorithm // IET Electric Power Applications. 2011. Vol. 5, Issue 2. P. 238. doi: https://doi.org/10.1049/iet-epa.2010.0147
31. Paikov I. A., Tikhonov A. I. Analysis of power transformer electromagnetic calculation models // Vestnik IGEU. 2015. Issue 3. P. 38-43. doi: https://doi.org/10.17588/2072-2672.2015.3.038-043
32. Ostrenko M. V., Tykhovod S. M. Calculation of losses in elements of construction of power transformers and reactors by finite element method with surface impedance boundary conditions // Electrical Engineering and Power Engineering. 2016. Issue 2. P. 33-42. doi: https://doi.org/10.15588/1607-6761-2016-2-4
33. Analysis of inrush currents of the unloaded transformer using the circuitfield modelling methods / Yarymbash D., Yarymbash S., Kotsur M., Divchuk T. // Eastern-European Journal of Enterprise Technologies. 2018. Vol. 3, Issue 5 (93). P. 6-11. doi: https:// doi.org/10.15587/1729-4061.2018.134248
34. Avetisyan D. A., Sokolov V. S., Han V. H. Optimal'noe proektirovanie elektricheskih mashin na EVM. Moscow: Energiya, 1976. 215 p.
