Motor Thermal Model
Motor Thermal Model
1. Abstract
This paper discusses the fundamentals of a motor thermal 12%, as in the Pulp and Paper industry. Downtime in a factory
model and its mathematical interpretation and physics for the can be very expensive and, in some instances, may exceed
different stages of motor operation. (overload, locked rotor, too the cost of motor replacement. Proper machine protection is
frequent or prolonged acceleration, duty cycling applications). required to minimize the motor failure rate, prevent damage
It explains Thermal Model Time Constants and other technical to associated equipment and to ensure both personnel safety
parameters that cause the biasing of the thermal model and production targets.
algorithm. Other topics covered in this paper show that (a) The document “Report of Large Motor Reliability Survey of
detailed motor data sheet information, and (b) coordination Industrial and Commercial Installations”, published by the IEEE
between the protection engineer and the motor supplier, can Motor Reliability Working Group [3] contains the results of IEEE
lead to proper selection of motor thermal protection parameters. and EPRI surveys on motor reliability and major causes of motor
This paper presents a closer look at motor stall, acceleration failure. The summary of these results is shown in Table I.
and running thermal limit curves. It also explains the concept
of thermal capacity and elaborates on how thermal capacity is In spite of different approaches and criteria (IEEE failure groups
evaluated in motor protection devices. The following points are are formed according to “cause of failure” and EPRI according to
also covered in this paper: “failed component”) both studies indicate a very similar failure
percentage associated with mechanical- and electrical-related
• Discusses some additional methods, such as voltage- machine problems.
dependant and slip-dependant motor overload curves, Analyzing the data from this table we can conclude that
employed to evaluate thermal capacity in many failures are directly or indirectly related to, or caused
nonstandard motor applications,
IEEE Study EPRI Study Average
• Presents the concept of matching thermal Failure Contributor % Failed Component % %
time constants for motor cyclic loads cases.
Persistent Overload 4.2% Stator Ground Insulation 23.00 Electrical
In addition, the response of a thermal Related
model algorithm in practical applications is Normal Deterioration 26.40% Turn Insulation 4.00
Failures
demonstrated. Bracing 3.00
Core 1.00
• Describes a real case example showing how Cage 5.00 33%
to apply and fine-tune the thermal model in Electrical Related Total 30.60% Electrical Related Total 36.00%
high-inertia load application. High Vibration 15.50% Sleeve Bearings 16.00 Mechanical
Related
Poor Lubrication 15.20% Antifriction Bearings 8.00
• Explores in this context, some of the key Failures
Trust Bearings 5.00
topics that will ensure safe operation of the
Rotar Shaft 2.00
motor while promoting satisfactory motor
Rotor Core 1.00 31%
design characteristics.
Mechanical Related Total 30.70% Mechanical Related Total 32.00%
High Ambient Temp. 3 Bearing Seals 6.00 Environmental
2. Introduction Abnormal Moisture 5.8 Oil Leakege 3.00
Maintenance
& Other
Abnormal Voltage 1.5 Frame 1.00 Reasons
Induction motors are the workhorses of any Related
industrial plant. Typical motor applications Abnormal Frequency 0.6 Wedges 1.00
Failures
include pumps, fans, compressors, mills, Abrasive Chemicals 4.2
shredders, extruders, de-barkers, refiners, Poor Ventilation Cooling 3.9
cranes, conveyors, chillers, crushers, and Other Reasons 19.7 Other Components 21.00
blowers. Statistics have shown that despite their Environmental Reasons & 38.70% Maintenance Related & 32.00% 35%
reliability and simplicity of construction, annual Other Reasons Total Other Parts Total
motor failure rate is conservatively estimated Table 1.
at 3-5% per year, and in extreme cases, up to Summary of IEEE and EPRI Motor Reliability Surveys.
ln(1 x ) | ln(1) x | x x 1
W (12)
? t max I | I 2 !! 1
I2
dT t I 2 t 1 (13)
dt W CM
Fig 2.
Single Constant Thermal Model vs Relay Overload Curve Comparison To gain insights into what the response is to a cycling load,
we will consider a simple cycling load in which the current
time-constant implied by the curve multiplier of the overload alternates between no load and an overload value:
curve must be set equal to the time-constant of the single time- Ilow 0 | current during the low cycle
constant model. In Figure 2, the ratio of the time divided by the
time-constant is plotted against per unit current. It can be seen I high current during the high cycle
that although equation (6) is not exactly the same as equation (14)
tlow time interval for the low cycle
(8), the approximation is very close, particularly for large values
of current. For lower values of current, the standard overload thigh time interval for the high cycle
curves are a better approximation to typical motor overload
curves than a single time-constant model. That is because Motor heating is proportional to the square of the current, so
there are at least two time-constants in the thermal response the effective current for heating over the cycle is:
of a motor. Over short time intervals, the thermal response
of a motor is dominated by heat transfer from the stator and
rotor conductors to iron. Over longer time intervals, the thermal
response is dominated by heat transfer from the iron to cooling (15)
air. A single time-constant model cannot be accurate over the
full range of operation and tends to overprotect a motor when
it is operated near its rated load. A standard overload curve
provides protection that is a closer match to a motor’s thermal
limit. The close proximity of the two curves for large values of Equation (15) can also be expressed in terms of a duty cycle
current is not a coincidence because both models are equivalent ratio:
to an adiabatic model for large values of current. This can be H effective D I high 2 1 D I low 2
shown mathematically by finding the asymptotic behavior of
the two curves. First, equation (8) is given approximately by: t high (16)
D duty cycle ratio
W CM t low t high
t max I
I 1
2
If the current and heating are expressed in per unit and the low
(10)
W CM cycle current is approximately equal to zero, the steady state
t max I | I !! 1
2
I2 boundary condition for tripping the motor becomes:
'Thigh |
1
W CM
I high
2
1 t high (20) • Motor Starting: State is declared when previous state was
“Stopped” and current greater than 2% of the motor full
load amps has been detected. The motor current must
The approximate temperature drop estimated by the cooling
increase to the level of overload pickup (service factor
model during the unloaded portion of the duty cycle is
times full load amps) within 1 second otherwise motor will
computed by multiplying equation (19) by the appropriate time,
transfer into the next state: “Running”
with per unit temperature equal to 1, because that is what it will
be approximately equal to during a limit cycle that approaches • Motor Running: State is declared when previous state was
tripping: “Starting” and motor current drops below overload pickup
1 level.
'Tlow | tlow (21)
W cool • Motor Overloaded: State is declared when previous state
The overload detection boundary is determined by setting the was “Running” and motor current raises above the overload
net temperature change equal to zero. This implies that the pickup level.
total of the right hand sides of equations (20) and (21) taken
together is equal to zero:
I high 1 t high
1 1
'Thigh 'Tlow 2
t low 0 (22)
W CM W cool
Equation (22) can be rearranged to show that standard overload
curves respond correctly to cycling loads. Equation (22) also
reveals how to properly select parameters for a load cycling
applications:
W cool
1 D I high
2
(23)
W CM
Equation (23) expresses the actual overload detection boundary
of an overload curve model in terms of its settings, the duty
cycle, and the amount of overload. Except for the factor
of ττ cool
CM
, equation (23) is exactly the same as the ideal overload
detection boundary, specified by equation (17). Equation (23)
and equation (17) will be identical, provided that ττ coolCM
is set Fig 3.
Motor Thermal Limit Curves
Motor Thermal Model Protection Applications 47
Thermal Capacity Used (TCU) evaluates the thermal condition Thus HCR is 8 sec / 10 sec = 0.8 and the level of stabilized TCU
of the motor. TCU is expressed as percentage of the thermal featuring the hot motor is equal to 20%, or in other words the
limit used during motor operation. Per IEEE Std 620-1996 (10) allowed motor thermal withstand time at overload conditions
the motor thermal limit is presented in the form of a time- will effectively decrease by 20%. If the motor load is lower then
current curve for 3 possible motor overload conditions: locked 100% the TCU level corresponding to the hot motor condition
rotor, acceleration and running overload. Every point on this is proportionally lower: 75% load – 15% TCU, 50% load – 10%
curve represents the maximum allowable save time at a stator TCU and so on.
current above normal load.
The unbalanced stator phase current will cause additional
rotor heating due to the developed negative sequence current
TCU is incrementally updated every 100 milliseconds and the
and flux rotating in the opposite direction to rotor rotation with
integrated value of TCU is stored in the thermal memory register
approximately double the power system frequency. The skin
of MPD according to the following equation.
effect in the rotor bars at this frequency will cause a substantial
TIME INTERVAL increase in rotor resistance and hence increased heating, which
TCU @ T TCU @ T 1 u 100% (25)
TIME TO TRIP is not accounted for by the regular thermal model. In order
to account for this additional heating factor the Equivalent
Current concept is introduced. The idea is that the current
The following example can be a good illustration of TCU
input into the thermal model is biased to reflect the additional
accumulation during the cold motor start; initial TCU is equal to
heating caused by the negative sequence component of the
0%. Motor starting pattern (1) and relay overload curve (2) are
load current.
shown at Figure 3.
For simplicity assume that the time interval for TCU update is 1 I EQ 2
IM u (1 K u ( I 2 I 1 ) 2 ) (26)
second. Every point of motor current on this plot corresponds
where:
to the number of seconds that motor can withstand before
I EQ - equivalent motor heating current
tripping on overload. The numerical values showing the progress
of TCU accumulation during 17 seconds of motor acceleration I M - real motor current
are presented in table 2. We can observe that by the end of a I1 - positive sequence component of real motor current
successful starting the thermal memory of the motor protection I 2 - negative sequence component of real motor current
device (MPD) accumulates 46.7% of TCU.
K - unbalance bias factor
The Unbalance Bias K factor reflects the degree of extra heating
caused by the negative sequence component of the load
current and can be defined as the ratio of Positive Sequence
Rotor Resistance to Negative Sequence Rotor Resistance. It
is practical and quite accurate to use the estimate method
to define the K factor. Equations for typical and conservative
estimates are presented below.
(27)
I eq
TCU END u 1 HCR u 100% (30)
SF u FLA
Where:
The RTD-TCU Curve is constructed based on the 3 key points. Let us assume that the following motor information is available
See Figure 4 for details. to us.
• Motor thermal limit curves are as presented at Figure 3.
1. RTD bias minimum
• Motor Cold and Hot Locked Rotor Times at 100% of the
Set to 40°C or another value of ambient temperature, if the system voltage are 34 and 26 seconds respectively. At 80%
appropriate RTD is available. TCU is equal to 0%. of the system voltage Cold and Hot Locked Rotor Times are
50 and 38 seconds respectively.
2. RTD bias mid point • Motor Acceleration at 100% of the system voltage is 17
The mid-point temperature is set according to the motor’s hot seconds. Maximum locked rotor current is 6 times that of
running temperature and is calculated as follows: full-load amperes (FLA). The MPD overload curve that we
employ as a limit to calculate TCU, is shown in Figure 3.
Rated Temperature Rise + Ambient Temperature Please note that the location of this curve is between the
For example: The temperature rise for NEMA Class B motors hot and cold thermal limit curves supplied by the motor
with 1.15 Service Factor, is 90ºC. Thus the temperature value for manufacturer. The time-current relation in this curve is per
this point is 130ºC. The TCU quantity for this point is the value of following equation:
a steady-state running condition @ rated motor load, and can
be found as: 87.4 u 12
Tto trip (sec) (32)
2
I EQ 1
TCU CENTER = (1 − HCR ) × 100% (31)
Fig 7.
Stall Trip. 80% Voltage
Fig 9.
Thermal model response to cyclic load
6. Application Description
This case study examines the Induced Draft (ID) fan application
Fig 10. on the A. B. Brown Unit 2 Selective Ccatalytic Reduction (SCR)
Voltage Dependent Thermal Limit Curves Project, located in Evansville, Indiana. Unit 2 is owned by Vectren
Corporation, and the role of Black & Veatch (B&V) on this project
For example, for a 100% voltage start (Figure 10, curve 1) the
was to construct an SCR facility in this plant.
locked rotor thermal limit is calculated based on a LRC of 6 times
full load current (FLC) and 8 seconds of the allowed locked rotor The SCR Project Scope of Work included modifying both ID fans
safe time, and I²T is equal to 288. After 14 seconds the motor for catalyst draft losses. The motors were powered from 13.8
accelerates to approximately 80% of the rated speed and the kV switchgear.
current drops to the level of 4.8 times that of FLC.
The allowed time to withstand 4.8 FLA for this stage of
acceleration is 40 seconds; I²T=922. Now let us consider the
The A. B. Brown case study presents a unique situation where the [13] GE Multilin, “Setting The 469 Motor Management Relay for
close coordination between the parties involved (the end-user, A Cycling Load Application”, GE Publication GET-8478.
the application engineer, the relay manufacturer and the motor
manufacturer) allowed the refinement of application-related
information so that proper thermal protection and coordination
during commissioning could be provided.