# Motor Parameter Estimation from Steady-State Models

Estimation of an induction motor's equivalent circuit parameters from steady-state models and manufacturer performance data is the most common approach for offline power system studies, e.g. dynamic motor starting, transient stability, etc [5].

## Contents

## Parameter Estimation Problem

The characteristics of an induction motor are normally provided by manufacturers in the form of a standard set of performance parameters, with the following parameters being the most common:

- Nominal voltage, [math]U_n \, [/math] (V)
- Nominal frequency, [math]f \, [/math] (Hz)
- Rated asynchronous speed, [math]n_{fl} \, [/math] (rpm)
- Rated (stator) current, [math]I_{s,fl} \, [/math] (A)
- Rated mechanical power, [math]P_{m,fl} \, [/math] (kW)
- Rated torque, [math]T_{n} \, [/math] (Nm)
- Full load power factor, [math]\cos{\phi_{fl}} \, [/math] (pu)
- Full load efficiency, [math]\eta_{fl} \, [/math] (pu)
- Breakdown torque, [math]T_{b} / T_{n} \, [/math] (normalised)
- Locked rotor torque, [math]T_{lr} / T_{n} \, [/math] (normalised)
- Locked rotor current, [math]I_{lr} / I_{s,fl} \, [/math] (pu)

We know that a set of equivalent circuit parameters can yield specific torque-speed and current-speed curves. So given a set of performance parameters that contain features on the torque-speed and current-speed curves (e.g. breakdown torque, locked-rotor current, etc), is it possible to determine the corresponding equivalent circuit parameters that yield these features? This is the crux of the parameter estimation problem and can be posed as follows - "How can the motor performance parameters be converted into equivalent circuit parameters?".

While all of the performance parameters in the above set can be used in an estimation procedure, there are actually only six indpendent magnitudes that can be formed from them: [math]P_{m,fl} \, [/math], [math]Q_{fl} \, [/math], [math]T_{b} \, [/math], [math]T_{lr} \, [/math], [math]I_{lr} \, [/math] and [math]\eta_{fl} \, [/math] [1]. These independent magnitudes will thus form the basis of the problem formulation, where the independent magnitudes calculated from the equivalent circuit are matched with the performance parameters supplied by the manufacturer.

The basic double cage model is used to illustrate how these six independent magnitudes can be calculated from the equivalent circuit model. Stator and rotor currents at slip "s" can be readily calculated from the equivalent circuit.

Quantities for per-unit active power [math]P \, [/math], reactive power [math]Q \, [/math] and power factor [math]\cos{\phi} \, [/math] at slip "s" can be calculated as follows:

- [math] S(s) = U_{n} I_{s}(s)^{*} \, [/math]
- [math] P(s) = T(s) (1 - s) \, [/math]
- [math] Q(s) = \Im \{ S(s) \} \, [/math]
- [math] \cos{\phi} (s) = \frac{\Re \{ S(s) \}}{|| S(s) ||} \, [/math]

Nominal speed [math]n_{s} \, [/math] and full load slip [math]s_{f} \, [/math] is calculated as follows:

- [math] n_{s} = \frac{120f}{p} \, [/math]
- [math] s_{f} = 1 - \frac{n_{fl}}{n_s} \, [/math]

where [math]p \, [/math] is the number of motor poles

- [math]f \, [/math] is the nominal frequency (Hz)
- [math]n_{fl} \, [/math] is the asynchronous speed at full load (rpm)

Calculating the slip at maximum torque [math]s_{max} \, [/math] is found by solving the equation:

- [math]\frac{dT}{ds} = 0 \, [/math]

(Under the condition that the second derivative [math] \frac{d^{2}}{ds^{2}} \lt 0 \, [/math])

In the double cage model, the solution to this equation is not trivial and it is more convenient to use an estimate, e.g. based on an interval search between s=0 and s=0.5.

## Problem Formulation Ignoring Core Losses

### Single Cage Model (Ignoring Core Losses)

In the single cage model, the locked rotor torque [math]T_{lr} \, [/math] and locked rotor current [math]I_{lr} \, [/math] are not used because the single cage model does not have enough degrees of freedom to capture both the starting and breakdown torque characteristics without introducing significant errors [1]. As a result, it is more commonplace to only consider the breakdown torque [math]T_{b} \, [/math] in the single cage model and simply ignore the torque and current characteristics at locked rotor. For wound-rotor motors, this yields sufficiently accurate results (i.e. in terms of the resulting torque-speed curve). However, a single-cage model is unable to accurately model the torque-speed characteristics of squirrel cage motors, especially those with deep bars.

Without taking into account core losses, the full load motor efficiency [math]\eta_{fl} \, [/math] also cannot be used. Therefore, there are only three independent parameters that can be used in the problem formulation: [math]P_{m,fl} \, [/math], [math]Q_{fl} \, [/math] and [math]T_{b} \, [/math].

These independent parameters can be used to formulate the parameter estimation in terms of a non-linear least squares problem, with a set of non-linear equations of the form [math]\boldsymbol{F}(\boldsymbol{x}) = \boldsymbol{0} \, [/math]:

- [math] f_{1} (\boldsymbol{x}) = P_{m,fl} - P(s_{f}) = 0 \, [/math]
- [math] f_{2} (\boldsymbol{x}) = \sin{\phi} - Q(s_{f}) = 0 \, [/math]
- [math] f_{3} (\boldsymbol{x}) = T_{b} - T(s_{max}) = 0 \, [/math]

where [math]\boldsymbol{F} = ( f_1, f_2, f_3 ) \, [/math] [math]\boldsymbol{x} = ( R_s, X_s, X_m, R_{r}, X_{r} ) \, [/math] are the equivalent circuit parameters of the single cage model

### Double Cage Model (Ignoring Core Losses)

In the double cage model, the locked rotor torque [math]T_{lr} \, [/math] and locked rotor current [math]I_{lr} \, [/math] are included as independent parameters. As in the single cage model, the full load motor efficiency [math]\eta_{fl} \, [/math] cannot be used without taking into account core losses. Therefore, there are five independent parameters and the following non-linear least squares problem:

- [math] f_{1} (\boldsymbol{x}) = P_{m,fl} - P(s_{f}) = 0 \, [/math]
- [math] f_{2} (\boldsymbol{x}) = \sin{\phi} - Q(s_{f}) = 0 \, [/math]
- [math] f_{3} (\boldsymbol{x}) = T_{b} - T(s_{max}) = 0 \, [/math]
- [math] f_{4} (\boldsymbol{x}) = T_{lr} - T(s=1) = 0 \, [/math]
- [math] f_{5} (\boldsymbol{x}) = I_{lr} - I(s=1) = 0 \, [/math]

where [math]\boldsymbol{F} = ( f_1, f_2, f_3, f_4, f_5 ) \, [/math] [math]\boldsymbol{x} = ( R_s, X_s, X_m, R_{r1}, X_{r1}, R_{r2}, X_{r2} ) \, [/math] are the equivalent circuit parameters of the double cage model

## Problem Formulation Considering Core Losses

Without taking into account the core (and mechanical) losses, the motor full load efficiency [math]\eta_{fl} \, [/math] cannot be used as an independent parameter in the problem formulation. This is because efficiency is calculated based on the ratio of output mechanical power to input electrical power. If the heat losses through the core and rotor frictional losses are not taken into account, then the equivalent circuit is not suitable to accurately estimate motor efficiency. It follows that attempting to use the motor full load efficiency in the estimation of the equivalent circuit without a core loss component would cause errors in the parameter estimates (e.g. the stator resistance would be overestimated).

When core losses are included in the model, then the motor full load efficiency [math]\eta_{fl} \, [/math] can also be used as an independent parameter. The problem formulations are restated below for the single cage and double cage models with core losses taken into account.

### Single Cage Model (With Core Losses)

The non-linear least squares problem for the single cage model with core losses is as follows:

- [math] f_{1} (\boldsymbol{x}) = P_{m,fl} - P(s_{f}) = 0 \, [/math]
- [math] f_{2} (\boldsymbol{x}) = \sin{\phi} - Q(s_{f}) = 0 \, [/math]
- [math] f_{3} (\boldsymbol{x}) = T_{b} - T(s_{max}) = 0 \, [/math]
- [math] f_{4} (\boldsymbol{x}) = \eta{fl} - \eta(s_{f}) = 0 \, [/math]

where [math]\boldsymbol{F} = ( f_1, f_2, f_3 f_4 ) \, [/math] [math]\boldsymbol{x} = ( R_s, X_s, X_m, R_{r}, X_{r}, R_{c} ) \, [/math] are the equivalent circuit parameters of the single cage model (with core losses)

### Double Cage Model (With Core Losses)

The non-linear least squares problem for the double cage model with core losses is as follows:

- [math] f_{1} (\boldsymbol{x}) = P_{m,fl} - P(s_{f}) = 0 \, [/math]
- [math] f_{2} (\boldsymbol{x}) = \sin{\phi} - Q(s_{f}) = 0 \, [/math]
- [math] f_{3} (\boldsymbol{x}) = T_{b} - T(s_{max}) = 0 \, [/math]
- [math] f_{4} (\boldsymbol{x}) = T_{lr} - T(s=1) = 0 \, [/math]
- [math] f_{5} (\boldsymbol{x}) = I_{lr} - I(s=1) = 0 \, [/math]
- [math] f_{6} (\boldsymbol{x}) = \eta{fl} - \eta(s_{f}) = 0 \, [/math]

where [math]\boldsymbol{F} = ( f_1, f_2, f_3, f_4, f_5 ) \, [/math] [math]\boldsymbol{x} = ( R_s, X_s, X_m, R_{r1}, X_{r1}, R_{r2}, X_{r2}, R_{c} ) \, [/math] are the equivalent circuit parameters of the double cage model (with core losses)

## Parameter Estimation Algorithms

The most common algorithms used to solve the non-linear least squares problems for motor parameter estimation are as follows:

## Software

Many commercial software packages for power systems analysis such as ETAP and DIgSILENT PowerFactory contain integrated routines for estimating motor parameters from manufacturer data. The commercial packages typically use some variant of the Newton-Raphson algorithm.

Moto is a free standalone tool for estimating induction motor parameters based on commonly available manufacturer data (i.e. breakdown torque, locked rotor torque, full load power factor, etc). The program supports a number of estimation algorithms including Newton-Raphson, Levenberg-Marquardt and genetic algorithms.

## References

- [1] Pedra, J., “Estimation of typical squirrel-cage induction motor parameters for dynamic performance simulation”, IEEE Proceedings on Generation, Transmission and Distribution, Vol. 153, No. 2, 2006

- [2] Filho, E., Lima, A., and Jacobina, C., “Parameter Estimation for Induction Machines via Non-Linear Least Squares Method”, Conference proceedings of IECON, 1991

- [3] Rogers, G., Shirmohammadi, D., “Induction Machine Modelling for Electromagnetic Transient Program”, IEEE Transactions on Energy Conversion, Vol. EC-2, No. 4, 1987

- [4] Waters, S. S., Willoughby, R. D., "Modeling Induction Motors for System Studies", IEEE Transactions on Industry Applications, Vol 1A-19, No. 5, 1983

- [5] Lindenmeyer, D., Dommel, H. W., Moshref, A., and Kundur, P., "An induction motor parameter estimation method", Electrical Power and Energy Systems, Vol. 23, 2001