Collective Losses of Low Power Cage Induction Motors—A New Approach

: Energy efﬁciency of systems of water pumping is a complex problem since efﬁciency of two distinct interacting systems needs to be combined: water and power supply. This paper introduces a non-intrusive method of calculating the so-called “collective losses” of a cage induction motor. The term “collective losses”, which the authors deﬁne, allows for accurate estimation of motor efﬁciency. Control system of a pump determines operating point of a pumping station, and thus its efﬁciency. General estimated performance characteristics of a motor, components of a control system, are assumed to serve selection of a range of pumping speed variations. Rotational speed has a direct effect on motor load torque, pump power and head, and thus on motor performance. Hellwig’s statistical method was used to specify characteristics of estimated collective losses on the basis of experimental studies of 21 motors rated at up to 2.2 kW. The results of simulations and experiments are used to verify validity and efﬁciency of the suggested method. The method is non-intrusive, simple to use, and requires minimum data.


Introduction
Analysis of energy efficiency of a water plant system is complex. Variations of plant output pressure and efficiency of pump and motor must be addressed. Energy efficiency is commonly decided by a control system applied to pump operation. In parallel, the control system is designed to provide adequate water pressure for each consumer. A pump and motor should operate in the region of maximum efficiencies, and thus of minimum losses. Minimum losses are a criterion of a control system's operation [1]. The losses are an indication of rotational speed and motor load torque. They are directly correlated with a plant's efficiency and output pressure. For relatively low pump capacity, rotational speed reduces relatively much, which affects the losses and motor efficiency [2].
Estimation of losses of a motor in operation (in situ) helps to develop an energy efficient control system of a pumping plant. Most traditional methods of efficiency assessment require measurements of rotor speed and shaft torque to compute output power. In most cases, torque and speed transducers cannot be installed because motors may be embedded in a machine where space is insufficient for attachment of such transducers between the motor and the load. The so-called non-intrusive motor efficiency regulator is the part of a control system that helps to improve energy efficiency without using a transducer of rotor speed and shaft torque. Therefore, indirect methods for estimating induction motor efficiency are preferable, since they do not require direct measurements of torque or speed [3,4].
The model of main losses of an induction machine (IM) consists of stator and rotor losses, mechanical losses, and core losses. Additional losses in effect of switching and Efficiency optimization can be realized by load observer control [25]. A motor equivalent model is used with online speed measurement. The load observer accurately estimates an unknown load disturbance. However, total motor loss consists only of stator and rotor copper loss and stator iron loss. Knowledge of parameters of the motor equivalent model is also required.
Efficiency is determined relying on equivalent circuit parameters, which are updated as a motor operates. At the same time, symmetry of these parameters is assumed for three phases of the motor [26]. In addition, power supply voltage, stator, and rotor resistances influence accuracy of efficiency calculations; what is more, this method does not cover stray load losses or mechanical losses. [27]. In effect, it is not applicable to industrial variable speed drives. An improved method of estimating parameters that serves identification of energy efficiency for an equivalent circuit is presented by [28,29]. Stray loss and friction and windage loss are derived from IEEE 112. A detailed method for determination of stray load losses is described in [30,31]. Estimation of iron losses in a motor under sinusoidal voltage supply plays an important part in the process of determining efficiency characteristics. Iron losses account for a substantial part of all the losses during frequency variations [32,33]. A new method of start-up helps to calculate motor efficiency. The test is very easy, practical, and quick to perform. It requires neither statistical nor nameplate data information. The results for the rated frequency are only described, since they are of no use to control systems [34].
The method of interpolation of total losses in induction motors in the power range 0.12-1000 kW is presented in the International Standard IEC 61800-9-2 [35]. Losses are calculated taking into account regression coefficients in terms of function of motor frequency and torque. The torque measurement for loss calculations is high intrusive. DTC or FOC electric drive control also enables torque determination. The high accuracy of the torque calculation (2%) according to the field-oriented methods makes them interesting for the determination of motor losses. The scalar control method cannot be applied. The efficiency interpolation data of typical motors two-and four-poles and efficiency class IE2 is provided.
A group of methods compute motor efficiency by means of slip measurement. They employ a linear dependence of output power on slip [36]. This method allows to estimate efficiency for rated frequency as a function of output power. High errors of efficiency estimation for other frequencies make it useless [37,38].
Another group of methods for determining efficiency of a squirrel cage motor relies on measurements of stator current [39,40]. They produce satisfactory results only for a limited range (0.7-1) of rated torque T N . Accuracy of efficiency calculation using this estimation method is not relevant.
There are many (over 30) methods for determining efficiency of a squirrel cage induction motor. A general division is based on: These methods generate results which estimate efficiency of induction motors only after performing additional investigations accurately or effectively. They all need no-load losses and stray losses to minimize the estimated efficiency error. This involves some intrusive steps. Table 1 shows a comparison of the methods analyzed above and what testing needs to be conducted for a given method, whether a method is highly or weakly intrusive, and what its accuracy level is. The authors review of control systems of existing and new pump and fan drives demonstrates they often do not apply energy optimization. In order to change this situation, control systems should be modernized without substantial spending or changes to design of electric drives. A modernized control system should address the following considerations:

•
A control system should operate in-situ, that is, without additional intervention in drive design. • A control system should not require additional preliminary motor testing, e.g., no-load test, or inductance measurement. • Determination of mechanical losses and stray load losses should not require additional measurements.

•
The above assumptions are adopted for estimation of IM losses and are intended to simplify implementation procedures of control systems to new and existing electric drives with any control techniques.
In order to calculate motor efficiency, a new non-intrusive method that refers to the so-called "collective losses" is postulated. The "collective losses", introduced by the authors, are not new losses but a group of total losses in a squirrel cage induction motor. The study is designed to determine a characteristic for their estimation as a function of frequency, voltage, and current supplied to a motor. Arguments of the estimation function are easily measurable, do not require high financial costs, and are non-intrusive. The main goal of the paper is to compare estimated and measured collective losses using the experimental model.
The newly proposed method of calculating the coefficients of explained variable, i.e., collective losses, is expressed as follows:

1.
Collective losses are defined on the basis of a Sankey diagram, 2.
Explaining variables necessary for a multiple regression model are determined on the basis of relative collective losses, 3.
Hellwig's method is adapted to eliminate those explaining variables that have less significant impact on the estimated collective loss function-the estimated model is simplified, 4.
Coefficients, components of the predicted variable, are determined.

Collective Losses-Definition
Induction motor losses for in-service intervals are evaluated further on for motors of rated power below 2.2 kW. Power loss components as described by the classical theory of electric machines are shown in Figure 1. Depending on their location in a motor, the losses are classified as Figure 1. Distribution of power loss components: P 1 -input power, P 2 -output power, P s -stator losses, P r -rotor losses, P m -mechanical losses, P mb -bearings losses, P mf -windage and fan losses, P Ψ -electromagnetic power, P js -stator Joule loss, P jr -rotor Joule loss, P ins -insulation losses, P cors , P corr -stator and rotor core losses, P strs , P strr -stator and rotor stray load losses. P s -stator losses, equal to P els resulting from adverse electromagnetic effects across the stator, P r -rotor losses, equalling P elr and resulting from negative electromagnetic effects across the rotor and mechanical losses P m (bearing P mb and windage and fan P mf losses).
A detailed analysis of losses in a motor is highly complex and irrelevant to their in-service evaluation. The authors suggest a new, simplified distribution of power losses across a squirrel-cage induction motor to calculate output power P 2 . The collective losses P coll are introduced as a new term whose interpretation is the purpose of this paper. They comprise parts of traditional induction motor losses. Core losses P cor , stray load losses P str , mechanical losses P m , and insulation losses P ins , calculated in complicated ways, are assumed to form a single grouping, i.e., collective losses P coll . The remaining losses, determined by distribution of power loss components in Figure 1, are stator winding (Joule) loss P js and rotor winding (Joule) loss P jr . These assumptions simplify Sankey diagram ( Figure 2) [41][42][43][44]. A sinusoidal supply voltage is assumed in this discussion.
In line with the modified Sankey diagram (Figure 2), collective losses P coll are The suggested method for estimation of collective losses relies on some assumptions: a. Accurate analysis of the individual losses is not the main goal. b.
We assume "a priori" that an estimation characteristic of collective losses can be produced. c.
The excess eddy-current losses are ignored. Finally, core losses P cor consist of stator and rotor hysteresis losses P hyss , P hysr and stator and rotor eddy-current losses P eddys , P eddyr . Analysis of collective losses in this paper is a new approach, therefore, simplifi-cation of this method is important. If the errors are not acceptable, a new assumption is intoduced. d. The insulation losses P ins are omitted as they are negligible for 2.2 kW motors. e.
Stator and rotor stray load losses P strs , P strr and stator and rotor stray no-load losses P str0s , P str0r are presented independently. This is connected with a broad range of load torque variations as defined by affinity laws for pumps and fans. f.
The motor is supplied with a sine wave and has a variable frequency voltage source. g.
The following relationship is realized by the voltage source: U f = const. A component of (1), rotor Joule loss P jr , equals slip s times electromagnetic power P Ψ : The following was produced on the basis of distribution of power loss components ( Figure 1): The following equation is obtained on the basis of (1), (2), and (3): The component including slip s at (4) becomes negligible small relative to P coll , which ultimately become: For nominal operating conditions of an induction motor, rated collective losses P collN are defined as where P 1N -rated input power, P jsN -rated stator Joule loss, n N -rated rotational speed, P 2N -rated output power. A comparison between distribution of power losses components in Figures 1 and 2 shows collective losses P coll consist of: P ins -insulation losses, P hyss , P hysr -stator and rotor hysteresis losses, P eddys , P eddyr stator and rotor eddy-current losses, P strs , P strr -stator and rotor stray load losses, P str0s , P str0r -stator and rotor stray no-load losses, and P mmechanical losses. P coll are expressed by: P coll = P ins + P hyss + P eddys + P hysr + P eddyr + P strs + P strr + P str0s + P str0r + P m Equation (7) can be simplified as the following assumptions are taken into account: • The rotor hysteresis losses P hysr are ignored because: f s f r . Such an assumption is correct for pump and fan load torque that reduces in line with the squared rotational speed in accordance with affinity laws.

•
The rotor eddy-current losses P eddyr are omitted because: The rotor stray load losses P str0r are ignored. They include rotor current I r , which is difficult to measure. However, I r is strongly correlated with the stator current I s . • Linearity of the magnetic circuit is assumed.

•
Magnetic flux densities B zs and B zr of the stator and rotor teeth are considered proportional to supply voltage U s and to the rotational speed n. • A heating power P heat component which defines temperature increment of the motor ∆T is added, representing effect of heat on P coll .
According to (7) and the above assumptions, the collective losses P coll can be shown in Figure 3. Figure 3. A block diagram showing parts of collective losses P coll , P cor -core losses, P hyss -stator hysteresis losses, P eddys stator eddy-current losses, P strs -stator stray load losses, P str0s -stator stray no-load losses, P m -mechanical losses, P heat -heating power. and can be represented as follows: where c hyss , c eddys , c strs , c str0s , c m , c τ -constant loss coefficients, c τ -temperature constant. ∆T = T f − T i ; T f -final stator winding temperature, T i -initial stator winding temperature. Constant loss coefficients have different values for each electric machine. The precise value of a constant helps to calculate selected loss correctly. Collective losses P coll represent a group of motor losses and can be estimated at an assumed accuracy. It means that other common constants can be computed for a given group of squirrel cage induction motors.
The above assumptions and set of Equation (8) are employed to convert Equation (7) into: In order to standardize the analysis and considering (9) and (6), relative collective losses P * coll are defined as follows: where relative collective losses P * coll are and relative constant loss coefficients c 1 , ..., c 6 are

Test Stand
Collective losses P * coll can be computed for each individual motor according to (10). It is assumed for the purposes of this paper there is a single characteristic of estimated collective loses for squirrel-cage induction motors rated below 2.2 kW. Such a characteristic can be generated on the basis of our testing of squirrel-cage induction motors. The test stand and the process of testing are described in the following part of this article.
Twenty-one squirrel-cage induction motors of both old and new motors and motors of various IE efficiency classes by six manufacturers were used in the experiments. Their rated powers ranged 0.75-2.2 kW and rated rotational speeds were 620-2900 rpm for frequency 50 Hz. Table 2 shows the rated data of some motors selected for the testing. A maximum rated power of 2.2 kW was adopted as the sole limitation. This assumption was expected to conform with testing conditions of pump energy efficiency as set out in [2]. Those authors analyzed pumps from four manufacturers with input power of up to 2.5 kW including squirrel-cage induction motors by unknown manufacturers. The sole assumption, of the maximum rated power of the motors, is reasonable. The testing was conducted in a laboratory stand shown in Figure 4. A synchronous generator GS was driven by a squirrel-cage induction motor M1 powered from a frequency converter ACS800. Change of frequency f p of U p at ACS800 output enabled variations of the M1 rotational speed and of voltage frequency f s across GS terminals. The excitation wiring was powered with voltage U extS from an SE source to arrive at a set voltage U g . Separation Delta-Wye transformer TR was installed at the output of the synchronous generator in order to eliminate the third harmonic of U g . Voltage U s of TR at f s was supplied to the motors tested, M2. In order to produce a desirable torque T across the shaft of M2, a DC separately excited machine DCM was coupled with M2. Winding of DCM armature was supplied with voltage U A and excitation winding was supplied with U extD from a thyristor converter DML [45]. Signal from a tachometer generator PT was provided to the DML. In effect, T across the motor shaft depended on rotational speed n of M2. A temperature sensor PT100 was mounted directly on the winding of each tested motor. To begin with, the load torque value is as desired. The motor temperature rises during the testing. In the initial time interval, close to the instant t 1 , the motor temperature becomes stable (temperature cannot change more than 1 • C in 1 h). This allows to measure all test values necessary to calculate collective losses. After these variable values are read, the load torque increases, then the motor temperature rises again till the next time interval close to the instant t 2 . With a further increase in torque, measurements are taken at subsequent instants t 3 , t 4 and so on. An example of temperature distribution and instants t 1 -t 8 of the measurements for the supply voltage frequency of 40 Hz are illustrated in Figure 5. Every instant t 1 -t 8 corresponds to a different load torque.  The tests were performed using the following measuring devices:

1.
Phase currents were measured with E3N current probes. Their accuracy was δ I = 0.03. A probe signal was converted into digital signal by a 16-bit converter in a NI-USB 6361 measurement card.

2.
Differential probes MTX9030-Z served to measure inter-phase voltages. The voltage measurements had a relative error of δU = 0.03 and signal damping 1/200. The probe output signal was sent to the analogue input of NI-USB 6361 16-bit measurement card.

3.
Rotational speed and motor torque were measured with MT-20 and MT-10 load cells with the accuracy 0.2% of the measurement range. 4.
MPS41 XX METROL mains parameter metre measured electric quantities supplied to the motor. The device analyzes the mains parameters with a satisfactory accuracy. The measurement errors are: current measurement error δ I = ±0.5%, voltages measurement δU = ±0.5%, active power measurement δ P = ±1%.
Over 1600 measurements were conducted for 21 tested motors. Measurement errors were below 10%. They were maximum for T lower than 50% of the rating, a consequence of measurement ranges of the metres which had been adopted.

Adapting Hellwig's Method
Test results of induction motors rated P N < 2.2 kW were the starting point for one estimated characteristic of collective losses P * coll . Selection of variables for econometric models was applied to this end. It was necessary to determine a computational method of finding a best-fitting subset of orthogonal linear combinations of all available predictor variables. Linear models or those that can be reduced to such models by single-equation methods are the best developed [46][47][48][49]. When there are two or more predictor variables, the model is called a "multiple regression model" with the general form: where y-the variable to be forecast (explained variable). X 1 , . . . ,X k -k predictor variables (explaining variable). Each of the explaining variables must be numerical.
t-time. β 1 , . . . ,β k -the coefficients measure the effect of each explaining variable after taking into account the effects of all the other explaining variables in the model.
The basic postulate concerning the variables' strong correlation between each explaining and explained variable and a weak correlation between the explaining variables causes most variable selection methods to employ some versions of the linear correlation. A solution found to be optimal depends on the lowest number of mutually uncorrelated explaining variables. The explaining variables show a maximum possible correlation with the explained variable [50].
The method of optimum predicate selection, adapting Hellwig's method [47], served to limit the number of elements of set X of originally selected explaining variables. Hellwig's method is most popular with empirical econometric research. Its idea is to select as explaining variables such variables that are strongly correlated with an explained variable and weakly correlated with one another. Estimating a matrix R (13) of correlation coefficients between potential explaining variables and vector R (13) of correlation coefficients between potential explaining variables and vector R 0 (14) of correlation coefficients between an explained variable and potential explaining variables is the starting point.
Correlation coefficients r x l x k between the particular explaining variables (X l , X k ) form the following matrix of correlation coefficients R: Correlation coefficients r yx k between an explained variable Y and the particular explaining variables X k form the following matrix of correlation coefficients R 0 : The correlation coefficients r yx k and r x l x k are measured by means of Pearson's coefficient of linear correlation r yx k between the explained variable Y and the explaining variable X k , expressed with [49,51,52]: The correlation coefficient r x l x k between two explaining variables X l and X k is defined as where N-number of measurements, x ik -i th value of the explaining variable X k , x il -i th value of the explaining variable X l , y i -i th value of the explained variable Y, and x k -mean value of the explaining variable X k , defined as x ik (17) y-mean value of the explained variable Y, defined as As part of Hellwig's method, all possible combinations of potential explaining variables are generated and then the so-called integral individual ratio of information capacity is studied, produced by application of each possible combination of these variables. Selection of explaining variables for the model according to Hellwig's method follows the procedure below.
In the first step, all possible L combinations are formed for all m explaining variables X 1 , ..., X m . The number of combinations L equals the number of all possible subsets of the m-element set, that is, L = 2 m − 1. For instance, if m = 2, 3 variable combinations can be produced: As part of step two, for each l th combination (l = 1, ..., L), an individual ratio of information capacity h lj (19) is calculated for j th explaining variable belonging to the l th combination of variables, as formulated below: where I l = {i : X i ∈ K l }-set of variable indices as part of l th combination K l ; h lj -individual capacity of j th variable in l th combination K l ; r yx k -correlation coefficient of k th explaining variable with the endogenous (explained) variable y; ∑ i=I l r ij -total sum of absolute correlation coefficients of j th explaining variable with the remaining explaining variables included in the combination K l .
The third step involves calculation of individual capacity of potential explaining variables (H l ) as the total sum of individual capacities as part of each combination, as formulated below: Explaining variables X belonging to that subset of l th combination to which a maximum value of H l corresponds are selected to the multiple regression model (12), i.e.,

Statistical Model of Collective Losses P * coll
The prediction of relative collective losses P * coll (10) is performed in this chapter. The linear multiple regression model in line with (12) was accepted for the purposes of analysis. Explaining variables X k (where k∈{1,2 . . . m}) are expressed with the quantities, parts of vector X, which are tested in the real model as a function of time [2]: Coefficients β 1 , . . . ,β k (12) correspond to the constants c 1 , ..., c 6 from (10); the definition is offered: where vector C consists of relative constant loss coefficients c 1 , ..., c 6 : Finally, the predicted variable Y as a linear multiple regression model is defined as The linear multiple regression model required optimization of the number of explaining variables and determination of coefficients c 0 ...c 6 . To this end, we measured the following electric and mechanical quantities as part of our experimentation: torque, rotational speed, voltages, currents, motor winding temperature, and frequency. About 1600 test samples were produced for 21 motors. Using the measurement data, the vector X explaining variables (22) was generated on the basis of (10) that defines relative collective losses and the multiple regression model (12).
Explaining variables X k for the predicted variable Y was selected according to the following formal and statistical criteria: • the explaining variables X k in the predicted variable Y should have coefficients of correlation ν k > 0.2-the method of eliminating quasi-constant variables, • the coefficients of correlation between the explaining variables X k , and the explained variable Y should tend towards one, • the coefficients of correlation between the explaining variables X k should tend towards one zero, • the degree of matching of the predicted variable Y to results of real measurements should be maximum possible, which is expressed as maximization of the determination coefficient R 2 .
The analyses showed the above criteria are met for all the explaining variables X k .
Hellwig's method was applied to the next stage of selection of the explaining variables X k . Based on the correlation coefficients r yx k (15) and r x l x k (16), Hellwig's ratio of integral capacity H l is analyzed for all combinations of l without repeating the explaining variables remaining in the model. L = 2 m − 1 of the combinations for m = 6 explaining variables equals 63. The individual ratio of information capacity for i th combination, K l , is calculated as per (19).
Out of all the combinations of explaining variables X k one l is selected for which Hellwig's ratio of integral capacity becomes maximum H lmax = max (H l ) (21). Following Hellwig's method, the number of explaining variables X k is reduced. Two explaining variables, X 1 and X 3 , remain in the model of relative collective losses for which the explained variable Y e is defined as follows: where c e0 , c e1 , c e2 -coefficients of explained variable Y e of relative collective losses P * coll . The simplified model of the predicted variable Y was subject to statistical verification to check if the selected variables X 1 and X 3 have significant impact on the explained variable Y and whether the criteria of applicability of the least-squares method to determination of the coefficients c 0 , c 1 , c 3 of the model are fulfilled. The model's statistical verification proceeded by:

1.
Testing significance of the model's coefficients with Student's t-test.

2.
Testing of the model's matching to empirical data: • calculation of the determination coefficient, • calculation of the convergence coefficient, • calculation of the random variability coefficient.

3.
Verification of the random structure's properties: • testing randomness of a random component with the runs test, • testing normalcy of a random component's distribution, • testing homogeneity of a random component's variance with White's test.
The analyses demonstrated the selected variables X 1 and X 3 have significant impact on the explained variable Y e . Finally, the explained variable Y e (26) depends on two variables, X 1 and X 3 , and is described with three coefficients: c 0 , c 1 and c 3 . The mean square approximation, also named the minimum chi-square method, where a minimum sum ∆ min is sought, served to compute these coefficients [48,49]: where i-number of a measurement, y mi -measured value of relative collective losses, y ei -value of predicted variable, estimated by the model (26) for the i th measurement. The minimum chi-square method reaches function factors by means of matrix calculus. The method can be applied if the following assumptions are met: a.
The model is linear or can be reduced to a linear form, b.
A linear dependence holds between no interpretative variables, c.
The number of measurements serving to estimate factors exceeds the number of factors estimated in this manner, d. The expected value of a random component equals zero E(ε) = 0.
In order to estimate the model coefficients (c e0 , c e1 , c e2 ), a coefficient vector C e is arrived at according to [51,52]: where Y m -vector of measured values y mi of relative collective losses, defined as C e -vector of the model coefficients is given as M-observation matrix of the explaining variables is defined as c e0 , c e1 , c e2 , computed as per (28), were substituted to (26) and values of the predicted variable Y were calculated as The predicted variable Y e considering (10) enables to define the model for estimation of relative collective losses P * coll : Relative estimation collective losses P * coll (33) were compared to measured experimental values.

Simulation and Testing Analyses
Variations of estimated relative collective losses P * ecoll were analyzed on the basis of calculation results as per (33). Based on (5) and (6), variations of tested relative collective losses P * tcoll expressed as follows could be analyzed: Calculations of estimated P * ecoll (33) and tested P * tcoll (34) relative collective losses required measurement of the following quantities: f -frequency of supply voltage, U ssupply voltage to induction motor, I s -stator current, R-resistance of stator winding, n-rotational speed, P 1 -input power, T L -torque across motor shaft, T-temperature of stator winding. Temperature is not directly present in the model of relative collective losses, nevertheless, it affects winding resistance.
To compare estimated P * ecoll and tested P * tcoll collective losses, losses of 21 motors operating at different rotational speeds and loads were tested. Rotational speed was in the range (0.2-1) n N . For selected motors, maximum rotational speed reached 1.2n N . A thyristor converter DML and a separately excited machine DCM helped to maintain a constant load torque T L for the tested motor M2 during a single test. The load torque ranged (0.1-1.1) T LN (T LN -rated load torque). The real range of load torque variations T L depended on the frequency of voltage supplied to M2. For lower frequencies, maximum load torque reduces if the motor current is below the rated value I N . The calculation results are shown in Figure 6. Tested P * tcoll collective losses are represented with black dots, while estimated P * ecoll collective losses are represented with red dots. The results form sets of dots around selected values of rotational speed n * . For instance, the set of dots around the relative rotational speed n * = 1 determines values of relative collective losses P * tcoll and P * ecoll for f = 50 Hz selected values of load torque T L in the range (0.1-1.1) T N . As f declines, corresponding to lower rotational speeds n * , maximum load torque reduces. The resultant black and red dot sets for a selected frequency largely overlap in Figure 6a,b,d. The results illustrated in Figure 6c are evidence of relatively high divergences.  Table 2).
The coefficient of determination R 2 is an indication of how correct the results are. R 2 for the entire test group of 21 motors is calculated as 80%. This means the selection of explaining variables is good. Higher values of R 2 represent quality of matching of an estimated model to a real motor. R 2 values for individual motors range 70-97%. This means the selection of explaining variables for estimation of collective loss characteristic curve P * ecoll is satisfactory, good, or very good. Examples of R 2 are presented in Table 2. The results can be seen as satisfactory. They reaffirm the prior assumption an estimated characteristic of collective losses exists in a squirrel cage induction motor is correct.
The statistical group encompasses 21 new and in-service motors by different manufacturers. This ensures universality of the results. The suggested approach is non-intrusive and can be used by industry to control on-line motor losses. To estimate collective losses, only measured values of frequency, supply voltage, and motor current must be taken into consideration. The instrumentation in this method is cost-effective, which justifies its particular applicability to low power pump drives.
The results suggest a need for further testing of 2.2-7.5 kW motors. They are part of many electric drives of pump systems. Existence of an estimated characteristic for collective losses shows motor efficiency needs to be verified. Such an estimated efficiency would be especially useful in control of water pump stations in order to improve energy efficiency.

Conclusions
In this paper, a non-intrusive method of estimating the so called "collective losses" of squirrel-cage induction motors of rated power 0.75-2.2 kW has been utilized. The term "collective losses" has been introduced. In order to determine coefficients of a collective loss curve, Hellwig's information capacity method has been used. The results have demonstrated the method's accuracy, with a good agreement between measured and estimated collective losses in an induction motor fed with a sinusoidal voltage source in the range of 10-60 Hz.
Determination of estimated collective losses is non-intrusive. It is necessary only to measure (33) frequency, voltage, and current of an induction motor. Additionally, temperature of the motor windings should be monitored. A comparison of this method of collective loss calculation with known techniques of motor loss calculation leads to the following conclusions: In the new method of collective loss calculation:

1.
No load test is required.

2.
IEEE 112 standard or other test methods are not used to determine stray load losses, 3.
A torque metre does not need to be mounted on the motor shaft-the method is non-intrusive. 4.
Any control system of electric drive can be used-V/f, DTC, FOC, VVC, or others, 5.
The estimated characteristic of the collective losses includes a minimum of explaining variables-Hellwig's optimization. 6.
The estimated characteristic of the collective losses addresses both new motors and motors after years of service. 7.
The estimated characteristic of the collective losses addresses motors of various efficiency classes. 8.
Collective loss calculation requires only measurement of supply voltage, motor current, and frequency.
These results are part of our research oriented towards improvement of pumping efficiency of water or sewage. If a satisfactory convergence is confirmed between estimated and measured characteristics of motor efficiency, the efficiency of an entire pumping unit can be controlled. Practical application of this solution to PLCs will improve energy effectiveness of new and in-service water and sewage pumping stations.

Data Availability Statement:
There is no report of any data.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
f frequency (Hz) f N rated frequency (Hz) I s current of stator winding (A) I N rated current (A) n rotational speed (rpm) P power (kW) P N rated power (kW) U s voltage across stator winding (V) P 1 input motor power (kW) P 2 output motor power (kW) P Ψ electromagnetic power (kW) P js stator Joule loss (kW) P jr rotor Joule loss (kW) P coll collective losses (kW) P collN rated collective losses (kW) P corr rotor core losses (kW) P cors stator core losses (kW) P strr rotor stray load losses (kW) P strs stator stray load losses (kW) P str0 stray no-load losses (kW) P m mechanical losses (kW) P ins insulation losses (kW) s slip T N rated motor torque (Nm) T L load torque (Nm) T LN rated load torque (Nm) U s voltage across stator winding (V) η motor efficiency P * coll = P coll