Alternative Computational Methods For Electromagnetic Circuits of Disk-Type Induction Motors - Comparison and Verification
Alternative Computational Methods For Electromagnetic Circuits of Disk-Type Induction Motors - Comparison and Verification
net/publication/351760722
CITATIONS READS
0 30
1 author:
Tomasz Wolnik
Łukasiewicz Research Network - Upper Silesian Institute of Technology
23 PUBLICATIONS 42 CITATIONS
SEE PROFILE
All content following this page was uploaded by Tomasz Wolnik on 06 July 2022.
Tomasz Wolnik,Ph.D.
Łukasiewicz Research Network – Institute of Electrical Drives and Machines KOMEL
Katowice, Poland
[email protected], ORCID ID 0000-0002-2636-8423
Abstract—Due to specific shape of electromagnetic circuit practise this may disqualify a particular design, since the
in disk-type motors, electromagnetic calculations should make time framework of the project is usually restricted. In
use of three-dimensional models. However, such calculations addition, time limitations occur, for example, at the stage of
are very time-consuming and this poses practical problems. submitting an offer to the customer. Alternative equivalent
The answer might be to use other alternative methods (not methods (other than FEM 3D) help to solve this problem,
based on 3D models), which would make it possible to obtain since they make it possible to obtain calculation results in a
appropriately accurate results. Comparison of calculation much shorter time.
results obtained with elaborated alternative computational Besides, in practise it’s also important that the input data of
methods for electromagnetic circuits of disk-type induction
the computational algorithm is based on the geometrical
motors is given in this paper. The results have been compared
with FEM 3D results and verified by lab test results (tests run
dimensions of the electromagnetic circuit. Then it is possible
with a physical motor model). When definite simplifications to analyze and monitoring the impact of design changes on
are adopted, the alternative computational methods make it the motor's operating parameters.
possible to reduce computational time in a significant way. The In [10,11] authors presents alternative calculation
achieved results are highly convergent with results acquired methods, but for axial flux permanent magnet motors
for real model. (AFPM), not axial flux induction motors (AFIM).
Keywords—disk-type motors, induction motors, In [6,12,13,14,15] authors presents the analytical method
computational methods. of calculating the electromagnetic circuit of the AFIM (Axial
Flux Induction Motor), but the input data are not the
I. INTRODUCTION geometric and winding data of the electromagnetic circuit.
In many applications disk-type motors present an Many physical quantities, such as electromagnetic force,
attractive alternative to cylindrical motors. Examples of efficiency, power factor etc. are assumed, not calculated on
specific disk-motor applications are lift drives, hoisting the basis of input data. This makes it impossible to analyze
machine drives as well as traction drives including vehicle the influence of design changes on motor parameters.
low-speed wheel drives etc. [1-6]. In most cases, brushless
In paper [13], the infinite permeability of the stator and
PM (Permanent Magnet) motors are used on account of their
rotor core was assumed, so the saturation of the magnetic
good parameters and resultant advantages [3,7]. The
circuit is not taken into account.
development trends in PM motors are strongly justified by
technical factors such as high efficiency, high torque An interesting approach is presented in [16]. For the
overload capacity, high power density coefficients. Apart calculations of the AFIM, the author used the
from all technical advantages, the PM motors possess two electromagnetic circuit of a classic radial flux induction
significant shortcomings: they are relatively expensive and motor (RFIM). However, the paper does not present the
the supply of main component (i.e. neodymium magnets) is method of converting the three-dimensional electromagnetic
strongly dependent on one monopolistic world source. In circuit to 2-dimensional and the obtained results were not
accordance with [8]c. 80% of neodymium magnet production compared with the test results on the physical model.
is centred in China. This is why even though PM disk motors
are being steadily developed, it’s important to research Comparison of the calculation results for single
simultaneously the construction of disk-type induction stator/single rotor disk-type induction motor (AFIM-11)
motors. obtained with alternative methods, i.e. the elaborated
analytical method and FEM (Finite Element Method) method
Due to specific flow of main magnetic flux in disk-type with 2D (2-dimensional) cylindrical model and standard
motors, where the length of magnetic path is subject to FEM 3D (3-dimensional) model is presented in this paper.
change when stator/rotor disk radii change, it is assumed that The developed alternative methods are based on geometric
electromagnetic circuits of these machines should be quantities of electromagnetic circuit and all winding data as
analysed with the help of three-dimensional models input quantities. These computational results have been
[1,5,9,10]. Then, the whole spectrum of different effects may verified by lab tests results (tests run with a real motor
be considered. However, the main drawback of such analysis model).
is long computational time. This is particularly true when
there is no possibility of limiting the computational model on
the basis of repeatable magnetic symmetry. Frequently in
𝑈 = 𝑈 + 𝐼 ∙ 𝑅 + 𝑗𝐼 ∙ 𝑋 , (1)
𝑈 = 𝑗𝐼 ∙ 𝑋 + 𝐼 ∙ , (2)
𝐼 =𝐼 +𝐼 , (3)
a)
a) b)
Fig. 5. Calculation model of AFIM11 motor, FEM 3D algorithm a) Motor
with stator windings in place, b) stator core and rotor core including end
ring
FEM calculations
alculations make it possible to determine
distribution of flux density in different constructional
elements of electromagnetic circuit. On the basis b of this
information, the designer is able to assess the saturation rate
in the machine elements for a given working point. In disk-
type motor the value of flux density in different elements is
closely associated
ted with diameter, so application of 3D model
mod
in calculations is very helpful. Distribution of flux density in
b) selected elements of electromagnetic circuit of AFIM-11 is
shown in Fig.7.
The final results of FEM 3D calculations are waveforms
of such physical quantities as electromagnetic torque,
current, supply voltage, flux linkage and stator iron core
losses. Example of calculated electromagnetic torque
waveform for the discussed motor is
i shown in Fig.6.
c)
Fig. 4. Stages of converting electromagnetic circuit of disk model to 2D
plane of cylindrical motor: (a) core of disk motor, (b) linear development of Fig. 6. Electromagnetic torque time waveform for the model motor; speed
the stator/rotor cross-section at Dsr diameter, (c) 2D model at a plane of n=950 rpm. [17]
cylindrical motor. [18]
The main comparison criterion adopted for assessment of
elaborated computational equivalent methods is accuracy of
calculations (in relation to real values obtained for the
physical model) as welll as time necessary for working out
the computational model together with actual calculation
time.
Computational models have been worked out using
computer equipped with Intel Core i7-4770K
i7 processor,
3.5 GHz, 32 GB RAM, licence restricted to a single core.
In order to compare the expenditure of time associated
with analysed methods, average calculation time for rated
working point of the motor for different FEM methods is
shown in Table II. Computational time for analytical method
is counted
ted in seconds and does not cover transient-type
transient time
analysis.
FEM 2d FEM 3D -
Fig. 7. Flux density distribution in selected elements of electromagnetic -RFIM AFIM
circuit of AFIM11 motor
Average calculation time 2.5 h 80 h
VI. COMPARISON AND VERIFICATION
CATION OF COMPUTATIO
COMPUTATIONAL
Number of calculation periods 6 6
METHODS
Final calculation time 120 ms 120 ms
In order to verify the presented computational algorithms
and to evaluate the equivalentvalent computational methods, Time step 0.2 ms 0.4 ms
a physical model of the electromagnetic circuit of AFIM11
motor has been worked out and executed. Lab tests have
been run. This physical model is shown in Fig.
Fig.8. Comparison Data on number of finite elements used in 2dRFIM
of the methods has been performed for rated parameters of model of the equivalent method and 3D model of the
the motor and in the operational range. discussed disk-type
type motor AFIM-11
AFIM is presented in Table III.
Distribution of finite element mesh for both models
mode is shown
graphically in Fig.9.
FEM 2d FEM 3D -
-RFIM AFIM
Band (rotor surroundings) 1 120 32 900
b)
Fig. 8. Physical model of electromagnetic circuit of AFIM11 motor,
a) b)
a) stator with windings, b) rotor core with aluminium cage
cage. Fig. 9. Finite element mesh of computational model a) FEM 2dRFIM,
b) FEM 3D AFIM-11 [8]
When we compare data set out in Tables II and III and As far as the rated parameters are concerned, we may argue
images shown in Fig.9, we may observe that even though that equivalent computational methods and analytical
total number of finite elements in FEM 3d AFIM model is method in particular, are characterized by very good
about 3 times higher, the density of finite element mesh is convergence with results obtained from physical model
significantly lower. Computational time for 3D model is tests. For this method relative error does not exceed 5% and
about 30 times longer even though computational time step this may be accepted as a very good result. In case of FEM
is twice as long as in 2D model (i.e. computational time step 2dRFIM method maximum relative error is 13.5% and this
is halved); we may state that it is about 60 times as long as relates to calculated load current I1. Still, at this error value
that of 2D model developed in equivalent plane of cylindrical
we may state that accuracy of calculations is good.
motor.
The remaining operational parameters of the motor are
Comparison of obtained results for different motor characterized by greater divergence from measured
parameters is shown in Table IV for all applied methods. The parameters. For FEM 2dRFIM method, the results in the
values of relative errors (calculated in accordance with operational range should be treated as rough estimates, since
Eq.(7) are given in Table V. relative errors for nearly all parameters attain 30 to 40%. In
| | case of analytical method this error is equal to 25%
δ[%] = ∙ 100 (7) maximum (this relates to slip corresponding to maximum
current). These differences are due to (among other things)
TABLE IV. COMPARISION OF RATED AND OPERATING PARAMETERS
FOR DIFFERENT MODELS
differences in calculated values of stator/rotor leakage
reactances and rotor resistance used in equivalent scheme
Rated and parameters.
Analytical FEM 2d FEM 3D
operational Tests
algorithm -RFIM -AFIM
parameters
T [Nm] 15 15 15 15.1 Apart from rated parameters and specific operational
parameters, operation of induction motor may also be
P1 [W] 1868 1832 1839 1882
described with selected curves. Comparison of
P2 [W] 1508 1514 1486 1520 computational methods for the most important curves is
I1 [A] 3.54 3.21 3.34 3.71 shown in Fig.10 [17].
n [rpm] 960 964 946 958
30
cosφ - 0.77 0.805 0.78 0.73
η % 80.7% 82.7% 80.8% 80.8%
25
15
sk - 0.163 0.17 - 0.13
Mr [Nm] 13.8 18.4 11.7 13.2 Lab tests
10
Ir [A] 13.7 15.1 11.6 11.5 Analytically
5 MES 2dRFIM
MES 3dAFIM
TABLE V. RELATIVE CALCULATION ERRORS FOR DIFFERENT MODELS 0
Rated and 900 920 940 960 980 1000
Analytical FEM 2d FEM 3D
operational n [ rpm]
algorithm -RFIM -AFIM
parameters
a)
T [Nm] 0.7% 0.7% 0.7%
P1 [W] 0.8% 2.7% 2.3% 6
P2 [W] 0.8% 0.4% 2.2%
I1 [A] 4.6% 13.5% 10.0% 5
ACKNOWLEDGMENT
The project has been financed by the Narodowe Centrum
Nauki (National Science Centre) within the framework of the
research grant No. UMO 2012/07/B/ST8/04099.REFERENCES