METHOD OF MONITORING A TEMPERATURE RISE OF A SQUIRREL CAGE INDUCTION MOTOR

Abstract: The invention relates to a method of monitoring a temperature rise of a stator winding of a squirrel cage induction motor, in which method a stator current (I) of the motor is measured and the temperature rise of the motor at each moment is determined on the basis of this current (I), when a rated temperature rise ( theta N) of the motor at a rated current (IN), the rated current (IN) and a temperature rise time constant ( tau ) of the motor are known. The method comprises steps of determining the stator current (I) of the motor as discrete moments separated by DELTA t long time periods, determining the temperature rise ( theta n(1)) of the motor during one time period (n) from the equation: theta n(1) = I<2>/I<2>N. theta N(1-e<- DELTA t/ tau >) determining a cooling ( theta n(2)) of the motor during the same time period (n) from the equation: theta n(2) = theta n-1e<- DELTA t/ tau ?, where theta >n-1< is the temperature rise of the motor at the end of the preceding time period, and summing the temperature rise ( theta >n(1)<) and the cooling ( theta >n(2)<) occurred during the same time period to obtain the temperature rise ( theta >N<) of the motor at the end of the time period (n).

Description: Method of monitoring a temperature rise of a squirrel cage induction motor The present invention relates to a method of monitoring a temperature rise of a stator winding of a squirrel cage induction motor supplied by a frequency converter, in which method a stator current I of the motor is measured and the temperature rise of the motor at each moment is determined on the basis of this cur- rent I, when a rated temperature rise Θ N of the motor at a rated current I N , the rated current I N and a tempera¬ ture rise time constant τ of the motor are known, the method comprising steps of determining the stator current I of the motor as discrete moments separated by Δt long time periods, determining the temperature rise θ n(1) of the motor during one time period n from the equation determining a cooling θ n(2) of the motor during the same time period n from the equation β _ a p -Δt/t β n {2) "n-1 where Q n _ 1 is the temperature rise of the motor at the end of the preceding time period, and summing the temperature rise θ n(X) and the cooling θ n(2) occurred during the same time period to obtain the temperature rise Θ N of the motor at the end of the time period n. Similar methods applicable to monitoring tem¬ perature rises, though not for frequency converter sup¬ ply, are known from GB Patent Application 2151862 and EP Patent Application 350 507, for instance. Electromechanical transmission generates power losses in a squirrel cage induction motor and the losses generated appear as temperature rises of the motor. A magnetizing current to be fed to the motor already causes a temperature rise and upon adding the load the motor warms up more. Part of the heat generated can be caused to flow out of the motor by means of cooling. The outflow of the heat depends on the cooling technique of the motor. The net temperature rise of the motor (tempe¬ rature rise caused by power losses minus heat outflow) determines the maximum allowed loading capacity. The more the motor is cooled, the more power losses may be allowed for it. The motors are provided with temperature limits, which should not be exceeded. The limits are generally informed both as maximum temperature rise of the insulation of the winding and as maximum temperature rise of the whole motor, whereby the endurance of bear¬ ings, for instance, is considered. If the motor exceeds the allowed temperature, the result is either a shorten¬ ed service life of the motor or, in the worst case, a damaged motor. The limits of the temperature rise of the motor thus restrict the loading capacity. This invention relates to an examination of a temperature rise at the generally most sensible point of the motor, the stator winding. Motor identification plates disclose how much the motor can be loaded by sine wave feed at rated frequency, but the situation becomes more complicated when the supply frequency is changed. Frequency converter supply causes more power losses in a squirrel cage induction motor than sine wave feed, due to which a motor supplied by frequency con- verter has a lower loading capacity than a motor having sine wave feed. Additionally, the cooling of the motor causes problems with loading capacity, when a frequency converter is used. The cooling capacity of air-cooled rib-cooled motors changes as a function of frequency, which also causes a change in the maximum loading capa¬ city of the motor as a function of frequency. The max¬ imum loading capacity of especially self-ventilated motors decreases considerably at low speeds when the ventilation capacity is weakening. Loading capacity is generally informed in rela¬ tion to rated moment. Figure 1 of the attached drawing shows a typical moment ratio curve of the loading capa¬ city of a self-ventilated motor when a frequency conver- ter is used, i.e. moment/rated moment in relation to the frequency. The first part (0...45 Hz) of the curve con¬ siders the effect caused by weakening ventilation capa¬ city and the latter part (50...100 Hz) of the curve con¬ siders the effect caused by a field weakening on load moment. The curve of Figure 1 is a simplification, which does not illustrate the actual loading capacity accurately. A more accurate idea of the effects of one frequency converter on the loading capacity is obtained from practical measurement results shown in Figure 2 of the drawing. When measuring the curve of Figure 2, the oper¬ ation of a modulator of the frequency converter used is based on a star modulated PWM (Pulse Width Modulation) technique. The winding behaviour of the loading capacity according to Figure 2 is caused, except by changes in cooling, by the modulator generating an output voltage. The loading capacity of the motor is changed at points where modulation is changed, being points at which the number of slices/pulses is changed. The slice number signifies the number of slices during a 60° sec¬ tor. The number of pulses depends, except on the number of slices, also thereon how many indicators are used per slice. The object of the present invention is to set forth an advantageous and simple thermal model for cur¬ rent, based on a measurement of the motor current and simulating an actual temperature rise of a stator winding of a squirrel cage induction motor reliably enough, by means of which model the temperature rise of the stator winding of the motor then can be reliably monitored. There are many kinds of protective methods for a squirrel cage induction motor on the market, but there are drives for which no one of the normally used methods is fully suitable. There are drives the motors of which can be reliably protected, but the protective methods are unreasonably expensive, while methods having an advantageous price are unreliable. The reliable methods are generally based on a measurement of motor winding temperatures. One object of this invention is to provide a method of monitoring a temperature rise, which method can be used sensibly also for motors in the protection of which it is not desirable (worth while) to invest much money, but the protective method must nevertheless be moderately reliable. An objective has been set that an equally good technical protection level shall be achieved within the whole supply frequency range of the frequency converter as is achieved at con¬ stant frequency by protective relays based on bimetal technique. When supplying a motor by a frequency con¬ verter, the above objects are achieved by means of the method according to the invention, which is character¬ ized in that the method comprises, except for the steps mentioned in the preamble, an additional step of multiplying the temperature rise Θ N of the motor by a multiplier K F , which is determined from the equation where Z I0 = maximum allowed continuous load current/ rated current of the motor at zero frequency, Z Iπιt = maximum allowed continuous load current/ rated current of the motor at a frequency f^, f act = actual frequency value of the stator current of the motor, f mt = frequency of the stator current of the motor, below which frequency the cooling capacity of the motor decreases and S yp = constant dependent on slice number/ modulation range of the frequency converter, this range corresponding to the actual frequency value f act . Further, for the purpose of considering the different temperature rise time constants of squirrel cage induction motors of different sizes, the method additionally comprises a step of determining the temperature rise time constant τ of the motor, used for the determination of tempera¬ ture rise and cooling and corresponding to the actual frequency value f act of the stator current of the motor, from the equation τ = τ N x f 1 + (Jc-1) ( a. ct f N f N = rated frequency of the stator current, τ N = temperature rise time constant at rated frequency and k = constant dependent on motor size. In the following, the method according to the invention and its theoretical background are described in more detail with reference to the attached drawing, in which Figure 1 shows a typical continuous loading capacity curve of a self-ventilated squirrel cage induc¬ tion motor supplied by a frequency converter. Figure 2 shows a loading capacity curve of a self-ventilated squirrel cage induction motor based on practical measurements of a frequency converter drive, Figure 3 shows a moment ratio curve of the loading capacity of a squirrel cage induction motor supplied by frequency converter and a current ratio curve formed on the basis thereof, and Figure 4 shows a loading capacity current curve of a self-ventilated squirrel cage induction motor sup¬ plied by frequency converter, which curve has been determined by measurements. On the basis of the requirements and objects constituting the starting point of the method according to the invention, an overload protection is implemented by a model calculating the temperature rise by means of mathematical, equations. In practice, the calculation is performed by a microprocessor, preferably by a micropro- cessor of the frequency converter or a drive control system. The method of the invention is based on simple basic equations of temperature rise, the number of input parameters being small and the values of the parameters generally known by the planners. The method of the in¬ vention concentrates on protecting the critical part of the motor, the stator winding. The basic starting point is the current to be fed to the motor. As known, the losses and the temperature rise of the motor are pro¬ portional to the second power of the current. The basic starting point of the method of the invention is con- stituted by the basic equations of temperature rise and cooling of a body: Temperature rise: θ = θ. d - e-^) ( 1 ) and cooling : θ = θ 0 e- c/τ ( 2 ) where θ = temperature rise "oo = temperature rise achieved after an infinite time θ = initial temperature rise τ = temperature rise time constant Interest is directed especially to the propor- tion of the temperature rise to the rated temperature rise Θ„ of the motor, corresponding to a temperature rise to be obtained at a rated point by a rated load, i.e. a rated current I N . Considering the fact that the temperature rise is proportional to the second power of the load current, the following equations are obtained for the temperature rise and the cooling: The load current going to the motor always warms up the motor, no matter how small it is. The cool- ing may nevertheless be stronger than the temperature rise caused by the load current, whereby the combined effect is cooling. The motor warms up and cools conti¬ nuously according to changes in the load current. Conse¬ quently, it is not possible to choose only one of the preceding equations 3 and 4. Both of them, cooling and temperature rise, have to be considered continuously. Temperature rise and cooling can be thought to occur as a non-continuous, i.e. discrete process. Accordingly, the temperature rise is divided as per time into succes- sive short parts and each discrete Δt long part of time (period) is examined separately. The current acting during this period warms up the motor during the whole period, while the temperature rise occurred before that period and caused by the currents of the preceding periods is cooling during the whole period to be examin¬ ed. In a way, this corresponds to a real physical situa¬ tion, in which both temperature rise and cooling take place in the motor all the time. The temperature rise achieved during one period n (the length of which thus is Δt) is obtained from equation 3 when time is indicated by Δt, and the cooling occurred during the same period n from the temperature rise of the preceding period n-1 is obtained from equa¬ tion 4 when time t is indicated by Δt. By combining these two equations, the temperature rise caused by the combined effects thereof is obtained at the end of period n: ~n ~ -n (l) + Q n {2 ) = - N (l-e- Δt/τ ) +6^ e "Δt/τ (5) IN Primarily, it depends on the length of the period Δt how big the temperature rise/cooling is during a period at a predetermined current. Secondly, the size of the temperature rise depends on the time during which the final temperature is achieved, i.e. on the tempera¬ ture rise time constant of the motor. When the tempera- ture rise is calculated at such short, i.e. discrete, time intervals, the actual value of the temperature rise is continuously known. Equation 5 does not consider the effects caused by a frequency converter on the loading capacity., for which reason it is suitable for being used as a thermal model of the stator winding of a motor fed by sine wave. As mentioned above, the loading capacity of a squirrel cage induction motor supplied by a frequency converter is not constant within the whole frequency range. Especially, as far as self-ventilated squirrel cage induction motors are concerned, the loading capa¬ city decreases abruptly because of weakened cooling capacity at frequencies lower than the rated frequency. In equation 5 is assumed that the rated temperature rise is achieved at the rated current. This is not the case when a frequency converter is used, but the rated tempe¬ rature rise of the motor is achieved even at the rated frequency at a current lower than the rated current. The changes in the loading capacity caused by frequency con- verter supply shall be considered in the thermal model of such drives. It is most preferable to process the changes in the loading capacity as a function of frequency into an equation form, which can be inserted as a coefficient in the thermal model equation 5. The loading capacity of a motor supplied by frequency converter is traditio¬ nally presented as a moment ratio according to a curve Z τ of Figure 3. The moment ratio curve of the loading capacity is there even simpler than the curve of Figure 1, due to which it suits well for a starting point of a coefficient equation. Magnitude Z τo of the curve Z τ of Figure 3 is the maximum allowed moment/rated moment at zero frequency in continuous operation that the motor endures thermally, Z Tmt is the maximum allowed moment/ rated moment at a threshold frequency f mt in continuous operation that the motor endures thermally, and f mt is the frequency below which the cooling capacity of the motor begins to go down. Curve Z τ of Figure 3 informs the loading capa- city as a moment ratio. However, the thermal model of current utilizes the measurement data of the current, due to which it is important in view of the current model to form a current ratio curve for the maximum allowed continuous loading capacity. Connection between current and moment varies greatly depending on the type of motor. At the same load moment, the currents taken by different motors differ from each other. This is due to the magnetizing current. The higher the pole number of the motor is, the lower is the power factor cosφ and the higher the magnetizing current of the motor. First- ly, the rated current of the motor shall be divided into a magnetizing current component 1^ and a current compo¬ nent I causing a moment reaction. The division is made by means of equations: T^ x [5 x ( sinφ ) - 1] J» - ' 5 - - s -i- nφ - "he < 6 > After this, a connection is obtained between the moment curve Z τ and the current curve Z z : Z j = - x yJ ( Z T x I qN ) 2 + ∑ϊw (9) -~N When the loading capacity decreases, the power factor cosφ decreases as well. A decreasing power factor is due to a decreasing current component causing a mo¬ ment reaction when the load is decreased. When the mo¬ ment ratio curve of the loading capacity is changed into a current ratio curve by means of equation 9, the changes in the current component, which influence the very moment, are considered and through this the changes in the power factor. When points of vertical intersec¬ tion of the moment ratio curve Z τ shown in Figure 3 are marked in the equation, the respective points of verti¬ cal intersection of the current ratio curve can be calculated. The load current curve Z x provided in this way is also shown in Figure 3. The points of vertical intersection of the curve are marked therein by indi¬ cations Z I0 and Z Imt corresponding to the points Z τo and z τmt °f *-he moπ^en ratio curve. In the case shown in Figure 3, motor identification plate values are used in equation 9, the values being: I„ = 242 A, cosφ = 0,85. The points of vertical intersection of the moment load curve Z τ are: Z^ _ 40 %/ z ^ ___ 90 % and ^ = 45 Hz# From Figure 3 is seen that the loading capacity decreases linearily at frequencies lower than the thre¬ shold frequency f mt , whereby the equation considering the changes in the loading capacity is an equation of a straight line: Zlact ~ ZlO = ∑ " ' 0 ( act ~ 0) (10) J -π.t u <S> f Ziac t = ~~ — x ( Z Imt - Z I0 + Z I0 ) , ( 11 ) Λ7t where z ιact = maximum allowed continuous load current/ rated current of the motor at actual frequency value ^ac ' Z I0 = maximum allowed continuous load current/ rated current of the motor at zero frequency, Z Iπ)t = maximum allowed continuous load current/ rated current of the motor at a frequency f mt , f act = actual frequency value of the stator current of the motor, f mt = stator current frequency of the motor, below which the cooling capacity of the motor decreases. Figure 2 showed that the loading capacity does not depend on the cooling only, but it is also affected by the quality of the voltage supplied by the frequency converter. The curve in Figure 4 shows how the modulator of the frequency converter influences the loading capa-. city, which curve is a loading capacity current curve of a self-ventilated squirrel cage induction motor determined by measurements and supplied by frequency converter. A significant change occurs in the loading capacity always at a change of the slice number of the modulator/the pulse number of the voltage. Figure 4 also shows the basic current curve shown in Figure 3 and based on equation 11. From Figure 4 is seen that, within each slice number area, the loading capacity follows the direction of the basic straight line rather accurately, but not the level of the basic straight line. In Figure 4 are drawn straight lines indicating the loading capa- city of each slice number area and being parallel with the above-formed basic straight line of the loading capacity. Accordingly, the actual loading capacity is approached very accurately by making at each slice number area a parallel transmission of the basic straight line solved above, which transmission corre¬ sponds to the loading capacity of the area. From Figure 4 is seen that the straight line must be lifted only at three slice areas of an upper modulation when the motor flow is kept at the maximum. In the other slice number areas, both in upper and lower modulation, the original straight line illustrates well the loading capacity. The areas marked by dotted lines in the figure represent areas within which the slice number is changed. The change is affected by hysteresis, due to which, when the reference frequency value is in¬ creased, the change of slice number occurs at a point slightly different from the respective changing point when the reference frequency value is decreased. How¬ ever, the frequency converter knows all the time the present slice number/modulation area, due to which the model can choose for its use a parallel transmission of correct size on the basis of this status data. A term S yp causing the parallel transmission is added to equa¬ tion 11 of the straight line, whereby is obtained: f a. ct 'lact x (Z f mt - z J I τn 0)> + '1 τn 0 + S yp (12) m. t The equation (12) can be made to the coeffi¬ cient of the thermal model equation 5 of current by raising it to the second power and by taking an inverse from it: In addition to the weakening ventilation capa¬ city of the motor, the coefficient considers the whole decrease of the loading capacity of the motor caused by the use of frequency converter compared to sine wave feed. This coefficient guarantees that the same tempe¬ rature rise is achieved at all points of the loading capacity current curve as by means of a sine wave fed motor, the load being a hundred per cent. Variables are the points of vertical intersection of the loading capacity curve and the term of the slice number area. The temperature rise time constant is one basic factor of the thermal model. So it is decisively, impor¬ tant for a reliable operation of the model to provide the thermal model with correct time constant. Measure¬ ments have shown that the temperature rise time con¬ stant of the motor is not constant in the whole fre¬ quency range. In the same way as the coefficient equa¬ tion was made of the changes in the loading capacity for the basic thermal model 5, it is also preferable to make a coefficient equation for the temperature rise time constant changing as a function of frequency. The temperature rise time constant of a 132 kW motor used for the measurements changes as a function of frequency according to equation 14 in such a way that the time constant is about double the rated temperature rise time constant at zero frequency. τ = • N K. = 'N 1 + l -act \ (14) f N The time constant of a second measured small (30 kW) squirrel cage induction motor follows equation 15 in such a way that the time constant is about four times the rated temperature rise time constant at zero frequency. τ = τ N x K τ = τ N x " 1 + 3 ( ^f- ) (15) ^N It is thus conceivable that the temperature rise time constant follows the equation f N = rated frequency of stator current τ N = temperature rise time constant at rated frequency and k = constant depending on motor size. On the basis of the above practical measure¬ ment results, the size of the constant k is chosen to be for instance 2 or 4 depending on the size category of the motor. If the thermal model is provided with a too low time constant, an error occurs always in a safer direction in view of protection. The model may decouple the drive too quickly. On the other hand, a too high time constant of the thermal model may cause a damage of the motor. By combining the above calculated coefficients caused by the loading capacity and the temperature rise time constant with the mathematical thermal model of current (equation 5) obtained for sine wave feed, a thermal model equation of current is obtained, which is suitable for frequency converter use: θ n = θ π(l) + -n (2) ^ e„ [i - e - 4t w'] + e n . i e - 4t " JA> (17) N As appeared earlier, the method of monitoring a temperature rise according to the invention can be arranged either to be calculated by the microprocessor of the frequency converter or by the microprocessor of the drive control system. Only a few parameter data are needed for the use of the method. As parameters are fed identification plate values of the drive and the point of vertical intersection of the loading capacity moment curve. In direct and group drives, software for the thermal model of current shall be positioned in the frequency converter. A direct drive is a frequency converter drive, which has a separate supply and one converter unit supplying one or several motors. The structure of group and line drives differs from the direct drive, but their basic structure is mutually similar. Inverter units of the drive groups have a common supply unit and a DC bar. The line drives differ from the group drives in such a way that the drive groups and thus also the motors therein are controlled by a common drive control. The control of the whole line drive is attended to by a process computer. Different drive applications also cause differences in the mea¬ surement of the motor current required by the model. When the inverter feeds only one squirrel cage induction motor, the separate current measurement system of the inverter can be utilized. In a so-called parallel motor drive, i.e. when more than one parallel-connected motors are fed by inversion, the separate current measurement of the inverter cannot be utilized. Then the current of (at least) one phase of each motor shall be measured by a separate current measurement system.

Claims: Claims : 1. A method of monitoring a temperature rise of a stator winding of a squirrel cage induction motor supplied by a frequency converter, in which method a stator current I of the motor is measured and the tem¬ perature rise of the motor at each moment is determined on the basis of this current I, when a rated temperature rise θ„ of the motor at a rated current I N , the rated current I N and a temperature rise time constant τ of the motor are known, the method comprising steps of determining the stator current I of the motor as discrete moments separated by Δt long time periods, determining a temperature rise θ n(1) of the motor during one time period n from the equation determining a cooling θ n(2) of the motor during the same time period n from the equation ft σ __(2) = ft lp e~Δt τ ' where θ n _.. is the temperature rise of the motor at the end of the preceding time period, and summing the temperature rise θ n(1) and the cooling θ n(2) occurred during the same time period to obtain the temperature rise Θ N of the motor at the end of the time period n, c h a r a c t e r i z e d in that the method additionally comprises a step of multiplying the temperature rise Θ N by a multi- plier K F , which is determined from the equation Kp - f. '10 \ Z 'I T m mt t- Z τn ) x act + S £ VP m, t where Z I0 = maximum allowed continuous load current/ rated current of the motor at zero frequency, Z τ _ = maximum allowed continuous load current/ rated current of the motor at a frequency f mt , f act = actual frequency value of the stator current of the motor, f mt = frequency of the stator current of the motor, below which frequency the cooling capacity of the motor decreases and S,^ = constant dependent on slice number/modu¬ lation range of the frequency converter, this range cor¬ responding to the actual frequency value f act . 2. A method according to claim 1, c h a r ¬ a c t e r i z e d in that it additionally comprises a step of determining the temperature rise time constant τ of the motor, used for the determination of tempera¬ ture rise and cooling and corresponding to the actual frequency value f act of the stator current of the motor, from the equation τ = τ N f. 1 + (Jc-1) ( act f N f H = rated frequency of the stator current, τ N = temperature rise time constant at rated frequency and k = constant dependent on motor size.

Inventors: ERKKILAE ILKKA, KAUHANEN MATTI, VEIJALAINEN VESA, ERKKILAE, ILKKA, KAUHANEN, MATTI, VEIJALAINEN, VESA
Applicants: ABB INDUSTRY OY, ERKKILAE ILKKA, KAUHANEN MATTI, VEIJALAINEN VESA, ERKKILAE, ILKKA, KAUHANEN, MATTI, VEIJALAINEN, VESA
Application year: 1994
Main classification: H02H6
IPC codes: H02H6, H02H7