International Journal of Electrical and Computer Engineering (IJECE) Vol.
No.
August 2017, pp.
ISSN: 2088-8708.
DOI: 10.
11591/ijece.
One-dimensional Lumped-Circuit for Transient Thermal Study of an Induction Electric Motor Radhia Jebahi.
Helmi Aloui.
Moez Ayadi Laboratory of Advanced Electronic Systems and Sustainable Energy.
University of Sfax.
ENETAocom.
Sfax.
Tunisia
Article Info
ABSTRACT
Article history:
Electrical machines lifetime and performances could be improved when along the design process both electromagnetic and thermal behaviors are taken into account.
Moreover, real time information about the device thermal state is necessary to an appropriate control with minimized losses.
Models based on lumped parameter thermal circuits are: generic, rapid, accurate and qualified as a convenient solution for power systems.
The purpose of the present paper is to validate a simulation platform intended for the prediction of the thermal state of an induction motor covering all operation regimes.
do so, in steady state, the proposed model is validated using finite element calculation and experimental records.
Then, in an overload situation, obtained temperatures are compared to finite elementAos ones.
It has been found that, in both regimes, simulation results are with closed proximity to finite elementAos ones and experimental records.
Received Aug 31, 2016 Revised Nov 2, 2016 Accepted Nov 16, 2016 Keyword:
Finite element Induction motor Lumped circuit Thermal analysis Thermal model Copyright A 2017 Institute of Advanced Engineering and Science.
All rights reserved.
Corresponding Author:
Helmi Aloui.
Laboratory of Advanced Electronic Systems and Sustainable Energy.
National School of Electronics and Telecommunications (ENETAoCo.
University of Sfax.
ENETAoCom.
BP 1163, 3018 Sfax.
Tunisia.
Email: helmi.
aloui@enig.
INTRODUCTION
The energy conversion inside an electric machine leads inherently to the heat occurrence, essentially due to different proportions of iron and copper losses.
In addition, an overheating of a motor is directly linked to the degradation of its performances .
Thus, in order to reduce such a limitation and avoid the damage of the device, it is necessary to gather the thermal analysis to the electromagnetic study during design process .
Indeed, the prediction of the temperature distribution in the different parts of the machine allows us to evaluate if the machine will reach the thermal class for which it is being designed and also to check if the supplied air flow of the cooling system is sufficient to a normal operation at rated conditions.
Thereby, a satisfactory life span of the machine gathered to an improvement of its performances could be guaranteed.
That explains the existence of several prediction methods of temperatures distribution inside electric motors, like: Computational Fluid Dynamics.
Lumped Parameter Thermal Circuits.
Finite Elements (FE) A, .
Thermal models based on lumped parameter were developed since a long time.
In such models, each node represents a part of the machine and the thermal circuit, in the steady state, consists of thermal resistances and heat sources connected between motor components.
For transient analysis, the heat capacitances are additionally used to take into account the ongoing changes.
Despite the accuracy of this method, it needs several improvement.
Indeed, in .
, the developed thermal model focused only stator Whereas, when the permanent magnets (PM) are placed in rotor, it is necessary to expand the thermal model and take into account this part.
In fact.
PMs are very sensitive to the temperature rise and the choice of their polarization is directly related to their thermal behaviour.
Moreover, in .
, the temperatures obtained by lumped parameter method were a little higher than the measured ones due to the hypothesis considered to elaborate the thermal model.
Journal homepage: http://iaesjournal.
com/online/index.
php/IJECE IJECE
ISSN: 2088-8708
However, the originality of this paper consists of the development of a generic refined model giving thermal state of all parts of the machine, along all operation regimes, in an attempt to a future exploitation to develop a control strategy with minimum losses rates and to size appropriate cooling system.
Presented works are applied on an induction motor and are organized as follows: the first part is devoted to the presentation of the developed lumped circuit for the thermal analysis of the considered device.
Then, the second part is intended to the validation of the proposed model through a finite element calculation using a 2D model and also an experimental test, considering a steady state operation regime.
Finally, the lumped circuit is exploited to predict the thermal state of the motor, when an overload regime is considered, and to compare obtained results to FE ones.
THERMAL MODELING OF THE INDUCTION MOTOR
The present section is devoted to the description on the induction motor going to be studied and to the presentation of two thermal models.
The first is based on lumped-circuit method and the second is intended to validate this latter through a finite element calculation.
Motor description The induction motor is the ultimate industrial motor due to its robustness and affordable price, compared to the rest.
Due to the hard constraints in industry, it is necessary to have an efficient tool telling about the thermal state of the machine, especially during overload or short-circuits situation, in order to guarantee equipment security and optimal control with reduced losses.
Indeed, losses in electrical machines lead to the rise of the components temperatures which conducts to the degradation of their performances.
Figure 1 shows the layout of the motor going to be modeled.
It is a cage asynchronous motor with: pole pair number .
, voltage .
V), power .
37 kW), power factor .
, frequency .
, weight .
, protection (IP.
and insulation class (F).
In order to prepare thermal modeling, the motor geometry is divided on layers as shown in Figure 2 and corresponding dimensions are measured, calculated and summarized in Table 1.
Figure 1.
Layout of the studied motor Figure 2.
Layers distribution for thermal modeling of the Table 1.
Measured motor dimensions Name Outer diameter rotor (R.
Inner diameter of rotor (R.
Air-gap thickness .
Inner diameter stator (R.
Insulator height .
Stator yoke height .
Outer diameter of stator (R.
Carter height .
Active length (L.
Value .
Lumped-circuit of the motor The 1D thermal modeling was used to integrate thermal aspects in the design of active and passive components in electronic circuits .
Thanks to the simplicity and efficiency of this technique, it becomes increasingly used for the analysis of thermal behavior of electric machine .
, .
, .
One-dimensional Lumped-Circuit for Transient Thermal Study of an Induction .
(Radhia Jebah.
A ISSN: 2088-8708 Thermal models based on lumped-circuit method are built using thermal-electrical analogy.
Indeed, this technique represents the thermal problems by the development of an electrical equivalent circuit using resistances and thermal capacitances .
Depending on the kind of thermal exchange, we can distinguish two types of thermal resistance: the conduction thermal resistance and the convection thermal resistance.
steady state, the thermal circuit is a network consisting of thermal resistances and heat sources connected between nodes representing motor components.
Besides, for transient analysis, thermal capacitances are additionally inserted to take into account the ongoing internal energy changes.
The general expressions of the conduction thermal resistance (Rcon.
and one of the convection thermal resistance (Rcon.
are respectively given by Equations .
Where: e is the distance between the point masses.
S is the interface area, is the heat conductivity.
Sc is the cooling cross section between the two regions and h is the convection coefficient calculated from proven empirical dimensionless analysis algorithms.
In addition, the thermal capacitance represents the ability of a component to absorb heat and is generally expressed as in Equation .
Where: V is the volume.
A is the density and c is the heat capacity of the material.
In different works .
, .
, the thermal studies were based on limited models focused only on stator components.
However, we propose in this paper a 1D thermal model.
Figure 3, allowing us to have real time information about happening temperatureAos changes in different motor components.
In this model, each component is represented by thermal resistances Ri and thermal capacitances Ci whereas heat sources correspond to copper and iron losses.
Since the structure of the induction motor is cylindrical, the thermal resistances of the different parts of motor are calculated after the resolution of the heat transfer Equation .
in the radial direction.
Taking into account the boundary conditions, the heat sources, and motor dimensions given in Figure 2 and Table 1, the thermal resistances and capacitances of the proposed model are expressed as in the following Equations.
Calculated values of those resistances and capacities are summarized in Table 2.
C Rotor C Coil C Insulator .
IJECE Vol.
No.
August 2017 : 1714 Ae 1724
IJECE
ISSN: 2088-8708
C Stator-yoke .
C Carter Table 2.
Calculated values of resistances and thermal capacity Motor components Rotor Coil Insulator Stator-yoke Carter Thermal resistance (W/K) Thermal capacity (J/K) Figure 3.
1D Thermal model of the motor implemented on Matlab-Simulink Finite element model of the motor The finite element method is used for numerical analysis of engineering problems (Mechanical.
Electromagnetic.
Thermal.
Biomechanics, etc.
) to obtain approximate solutions with very satisfying accuracy and closed proximity to the reality.
This technique is increasingly being employed in the thermal analysis of electrical machines .
, .
, .
Several software platforms are devoted for this kind of calculations as long as the study domain is defined around specific model and boundary conditions.
In the present case, the motor is supposed to be thermally isotropic along active length direction.
Thus, a 2D finite element model is sufficient to carry out the motor thermal state for any operation regime.
Figure 4 shows the 2D study domain, the cooling system is modeled as a forced convection through the boundary condition .
and the temperature rise in air-gap caused by the copper losses in rotor is taken into account as a boundary condition .
In addition.
Table 3 shows the thermal proprieties of motor components Simulation results and discussion At this stage, we are interested in the comparison of temperature evolutions inside the motor obtained respectively by the lumped circuit and the FE model for both a static study and a transient one.
Static analysis The static study is conducted in an attempt to focus the maximal temperature that could be reached in each element of the motor.
Such a temperature is obtained when steady-state is reached, after a calculation time of 4200 s.
In fact, we illustrated in Figure 5 the evolution of temperature inside the motor armature.
this figure we can see that, at the starting point the temperature is uniform and equal to the ambient One-dimensional Lumped-Circuit for Transient Thermal Study of an Induction .
(Radhia Jebah.
A ISSN: 2088-8708 temperature .
Then, after few operating seconds, the obtained distribution is illustrated in .
These temperatures rise gradually according to the nodeAos position to reach finally the steady-state .
, where the temperature is maximal and uniform.
As an illustration.
Figure 6 shows the temperature evolution of the stator coil .
and the stator yoke .
versus the stator copper losses carried out using lumped-circuit .
and FE model .
Shown results are the extrapolation of the three marked points in the mentioned figure, corresponding to an average ambient temperature of 29 AC and the calculation conditions itemized in Table 4.
We can remark the total harmony between lumped-circuit values and FE ones.
Table 3.
Thermal properties of the motorAos materials Material Copper (Coi.
Insulator Aluminum (Carte.
Iron (Statoryok.
Conductivities (W.
K-.
Mass heat (J.
K-.
Density (Kg.
Figure 4.
Meshed 2D finite element thermal model Figure 5.
Temperature distribution in stator coil Figure 6.
Evolution of the stator coil temperature .
and the stator yoke temperature .
versus the stator copper losses.
Legend: .
: Lumped-circuit results and .
: FE results IJECE Vol.
No.
August 2017 : 1714 Ae 1724
IJECE
ISSN: 2088-8708
Table 4.
Simulation conditions of the static study Conditions Iron losses .
Stator copper losses .
Rotor copper losses .
Point 1 Point 2 Point 3 Transient analysis In order to prove the accuracy the developed lumped-circuit, the thermal state of the motor is carried out in transient regime as shown in Figure 7 for the motorAos armature.
Results are obtained for when the following conditions are considered: ambient temperature: 28.
2 AC, iron losses: 26,5W, stator copper losses:
45W and rotor copper losses: 25.
79 W.
A good accuracy is noticed between the lumped-circuit values and the FE ones where the introduced error does not exceed 5% in the worst case.
Moreover, it is to be signaled that, from the point of view of simulation cost, the 1D thermal model is extremely fast compared to the FE one: A CPU-cost ratio of .
is registered between the two simulations when they are performed with a Pentium processor.
Thus, the lumped-circuit could be convenient solution to satisfy control requirement consisting of real time information about happening temperatureAos changes in different motor components.
EXPERIMENTAL VALIDATION
As mentioned previously, our purpose is to establish a simulation platform prediction of the thermal state of an induction motor covering all operation regimes.
To do so and before expanding our works to other rygimes, we intend in the present section to validate experimentally previous presented models .
umpedcircuit and FE mode.
Two tests are performed: at no-load then and under load operation regimes, considering the conditions given in Table 5 and using the built test bench of Figure 8 composed of:
The induction motor: on which thermal measurement are performed An autotransformer: gradual star of the induction motor The motor load: A coupled DC generator feeding a resistive load A digital power meter: measurement of electrical magnitudes .
ower, current and voltag.
An infrared thermometer: to pick up the temperature in each component of the motor, where: Xr is the measurement point on the rotor.
Xco is the measurement point on the armature.
Xsy is the measurement point on the stator yoke and (Xca1 and Xca.
are two symmetric measurement points on the carter.
addition, in order to have access to the internal parts, two apertures are made in the motor carter.
The first one, .
, serves for the temperatures of coil and stator-yoke, while the second, .
, is used for the rotor's temperature.
Figure 8.
At this stage, it is to be mentioned that each experimental lasts for 4200 seconds and the temperature of each element is measured every minute.
Temperature(AC) Table 5.
Energetic conditions of experimental tests I(A) Pf(W) Pjs(W Pjr(W No-load Load operation 100(W) 400(W) Time.
Figure 7.
Evolutions of the stator coil temperature No-load operation study The energetic state of the motor corresponding to the no-load operation test is given in Table 5.
One-dimensional Lumped-Circuit for Transient Thermal Study of an Induction .
(Radhia Jebah.
A ISSN: 2088-8708 The evolutions of the picked temperatures in the carter symmetric points (Xca1 and Xca.
are illustrated in Figure 9.
We can observe that the temperature is uniform both point which permits to verify the symmetry hypothesis used to elaborate the 1D thermal model and the FE one.
Furthermore.
Figure 10 and Table 6 show a good accuracy between the experimental measurements of the coil temperature and its values obtained by the 1D thermal and the finite element models with an error less than 2.
Figure 8.
Experimental test bench Temperature(AC) Xca1 Xca2 Time.
Figure 9.
Evolutions of the temperatures in carter under no-load operation Temperature(AC) Time.
Figure 10.
Evolutions of the coil temperature during no-load operation IJECE Vol.
No.
August 2017 : 1714 Ae 1724
IJECE
ISSN: 2088-8708
Load operation study The energetic state of the motor corresponding to the present test is given in Table 5.
Figure 11 shows the evolutions of the armature temperature picked experimentally and evaluated using the lumped-circuit and the FE model considering two values for the load connected to the DC generator .
W and 400 W).
Referring to such figure and also to Table 7, we can notice a high accuracy between the experimental measurements of the coil temperature and its values obtained by the 1D thermal and the finite element models with an error less than 2.
Table 6.
Comparison results (No-load operatio.
Time .
Measured Temperature(AC) Error % 400(W) 100(W) Time.
Figure 11.
Evolutions of the coil temperature during load operation Table 7.
Comparison results (Load operatio.
Tim e .
Measured 100(W) Error % Measured 400(W) Error % One-dimensional Lumped-Circuit for Transient Thermal Study of an Induction .
(Radhia Jebah.
A ISSN: 2088-8708 All performed tests, at no-load and under load regimes, prove a high satisfactory accuracy level between thermal states of the induction motor obtained through experimental measurement and carried out using the established models .
umped-circuit and FE mode.
Consequently, we can qualify both models as valid tools to be used for the exploration of the thermal state of the motor, especially when experimentation is either difficult to implement or harmful as in overload operation regimes like short-circuits under nominal voltage supply.
OVERLOAD OPERATION
The motor components are subjected to the load and temperatures variations according to the operating cycle and ambient conditions.
They are stressed heavily due to the effects of temperature cycle and face to the overheating risk that can cause the degradation of the motor performances and a raise of losses which must be kept below critical values during all operation regime.
In fact, a raise of temperature level caused by an overload operation changes coils parameters .
esistance and inductanc.
which leads to wrong control instructions in numerous strategies where such parameters are considered fix.
Figure 12, shows the evolution of the armatureAos resistance versus the temperature of the studied induction motor, using Equation .
Where:
R (T): the resistance at the temperature T R (T.
: the armature resistance at the ambient temperature T0.
In the present case.
R(T.
has been measured at T0=29.
9AC and was equal to =14.
: the temperature coefficient (=3.
9*10-3 C-.
Inconvenient control orders are the origin of several phenomena limiting the device lifetime like vibration and mechanical fatigue.
Thus, a predictive rapid model providing real-time information about the thermal state of the motor with acceptable accuracy could be a good solution permitting to perform advanced and complete control schemes taking into account the evolution of the machine parameters and the prevention of damages.
For this reason, it is necessary to consider a coupled problem based on an electromagnetic model, for losses quantification, and an electro-thermal one in order to carry out real-time information about the thermal state of the motor with acceptable accuracy leading to adjusted d-q control currents.
Figure 13.
Temperature (AC) Figure 12.
Coil resistance versus motor temperature Figure 13.
Control algorithm The present work focuses the electro-thermal part of the coupled problem by the establishment of a refined lumped circuit giving accurate results, as proved previously through the comparison with FE values and experimental records.
Now, we are going to test the ability of such model to detect sudden variation of IJECE Vol.
No.
August 2017 : 1714 Ae 1724
IJECE
ISSN: 2088-8708
the coil temperature when an overload situation occurs.
On another hand, the effect of the subdivision of the coil block into multiple layers, as shown in Figure 14 and Figure 15, on values accuracy is investigated.
Figure 14.
CoilAos block Refined Coil Figure 15.
1D thermal model with refined coil block Using the model of Figure 15, the evolution of coilAos temperature is carried out in the case of an overload: a short-circuit under nominal voltage with an armature current equal to ten times the nominal one and corresponding to: iron losses .
5 W), stator copper losses .
0 W) and rotor copper losses .
W).
Moreover, simulations are conducted for a number coilAos blocks equal to: 1, 4, 8, 16 and 24, respectively, as illustrated in Figure 16 and summarized in Table 8.
Furthermore, a FE calculation is performed in order to validate lumped-circuit results.
Figure 17.
Analyzing obtained results, we can notice that:
The accuracy of the lumped-circuit follows the increase of the number of coil blocks, despite an increase of simulation time.
For that, a compromise should be always found between calculation time and tolerated error.
For all considered coil blocks numbers, the error does not exceed 7%.
Thus, we judge not necessary to refine the model of Figure 3 and to assign only one block to the coil as the rest of motor components.
Figure 16.
Evolutions of the coil temperatures under overload .
D thermal result.
Figure 17.
Temperature distribution in the coil at t=0.
Table 8.
Coil temperature at 0.
Number of divisions Coil temperature (AC) Percentage error % CONCLUSION Electrical machines performances are strongly affected by the thermal states of their components.
Thus, the development of tool giving real time information about the temperatures of the deviceAos elements is One-dimensional Lumped-Circuit for Transient Thermal Study of an Induction .
(Radhia Jebah.
A ISSN: 2088-8708 a recommended solution, permitting an appropriate control .
n normal regime.
and also to predict and avoid damages .
n overload case.
Thermal models based on lumped parameters are known by their accuracy, but several of them focus only stator parts.
Besides, some model didn't provide the optimal accuracy and rapidity.
The originality of this paper consists of the development of a generic refined model giving thermal state of all parts of the machine, along all operation regimes.
A lumped-circuit was proposed and validated, through a first study, at no-load and under load regimes.
Indeed.
All performed tests proved a high satisfactory accuracy level between thermal states of the induction motor obtained through experimental measurement and carried out using the established models .
umped-circuit and FE mode.
, with an error less than 2.
Consequently, both models .
umped and FE) were adopted to predict thermal state of the motor in overload operation regime.
A short-circuit situation, under nominal voltage and ten times the nominal current, was It has been found that, for all considered coil blocks numbers, the error does not exceed 7%.
Thus, in order to avoid the increase of simulation time with no significant improvement of the accuracy, it is not necessary to refine the model and only one block will be assigned to the coil as the rest of motor
To exploit the developed model, and based on registered progress in software tools, authors are actually working on the development of an advanced platform intended for a control of an electrical machine control integrating a real-time thermal module telling about the thermal state of the machine-converter REFERENCES