Induction Machine On-Line Parameter Identiﬁcation for Resource-Constrained Microcontrollers Based on Steady-State Voltage Model

: This paper presents a new, computationally modest on-line identiﬁcation method for the simultaneous estimation of the rotor resistance and magnetizing inductance of an induction machine suitable for electric drives that use an indirect ﬁeld-oriented control strategy (IFOC), and their control hardware is equipped with a resource-constrained microcontroller. Such drives can be found both in the manufacturing industry and railway traction vehicles in the thousands, having either older control hardware that cannot cope with computationally excessive identiﬁcation methods or being in cost-sensitive applications, thus being equipped with a low-cost microcontroller. IFOC is a very common control strategy for such drives due to its good dynamic properties and comparatively simple implementation. However, it is sensitive to inaccuracies of rotor resistance and magnetizing inductance. These two parameters change during the operation of the drive, being inﬂuenced by the temperature, frequency, and saturation of the magnetic circuit. Improper values of parameters in the controller can degrade the performance of IFOC, resulting in slower acceleration or unnecessary oversaturation of the machine. Respecting these changes can therefore bring signiﬁcant beneﬁts such as the good dynamic properties of the drive, which can shorten operations in the manufacturing industry or travel times of vehicles. A number of on-line identiﬁcation methods for monitoring the parameter changes have been published so far, but the majority of them are demanding on microcontroller time or its memory. The proposed method, on the contrary, is comparatively simple and thus easy for implementation with low requirements to the microcontroller. Therefore, it is suitable for both upgrades of existing drives or new low-cost applications. Derivation of the method from the mathematical model and the ﬁnal algorithm for the microcontroller are presented. The performance of the proposed method is validated with experimental results obtained with a 3.5 kW induction machine drive with an industrial microcontroller during a warming test and under various loads and frequencies.


Introduction
Induction machines have many favorable properties such as robustness, reliability, and low maintenance cost, which have caused them to become one of the most widespread types of electromechanical converters in industrial and railway applications. Probably the most significant difficulty of induction machines is the complexity of their control when they are used in a variable speed drive. Among various modern methods of control of induction machines, Field-Oriented Control (FOC), Model Predictive Control (MPC), and Direct Torque Control (DTC) are the most frequently used ones [1][2][3][4][5][6][7][8][9]. All these methods of control have a common feature in that they rely on some form of mathematical model of the induction machine.
resources. It can be so demanding that for its implementation a microcontroller with a real-time control accelerator may be needed [6]. GAs are demanding on the storage of historical data, while OAs, such as Particle Swarm Optimization, need a large number of calculations [23].
Different types of observers are another vast group of methods used. The Luenberger observer was used for flux estimation with an ensuing stator and rotor resistance calculation in [24] and for a stator resistance estimation in [25,26]. Sliding mode observer usage for stator resistance estimation has been proposed in [27]. An adaptive observer for stator and rotor resistance was utilized in [28]. These observer-based methods need matrix calculations in their algorithms, which can be a limiting point for their deployment in a microcontroller with resource-constrained hardware.
The Extended Kalman Filter (EKF) has also been widely studied for parameter identification [29][30][31], but this method is also characterized by its excessive computational demands [5].
Recursive Least Square (RLS) [10,32,33] is easy to implement and is not so computationally intensive, but it is very sensitive to noise [34], which implies the implementation of more sophisticated filters. Additionally, the tuning of the forgetting factor of RLS can be tricky [5] and can lead to error accumulation. To overcome this, a time-varying forgetting factor has been proposed [35], which, however, makes the RLS more complex.
Model Reference Adaptive System (MRAS)-based methods can be also counted among computationally undemanding techniques with a simple structure and implementation [21,22,36,37]. Usually, they are used to identify one or two parameters because it is difficult to find a stable adaptive law [5]. The adaptive law is often implemented by PI or I controllers, which can be seen as a burden from the point of view of limited-resources hardware because it needs an additional controller in the structure.
Signal injection is generally not demanding on the microcontroller, but it brings other related difficulties with its implementation. The use of a high-frequency signal is not so popular because of its demands on the speed of the switching devices of the converter [5]. DC signal injection is thus more favorable; however, it can cause unwanted additional torque pulsations. These could be avoided, for instance, by introducing second-order harmonics [38], which can lead to more complicated structures in the microcontroller. Further, it should be used only in specified time spans to avoid unnecessary additional power losses caused by the DC test signal [39].
Among the methods that do not require excessive microcontroller resources are those based on voltage model calculation. In [40], the authors used voltage-model-derived differential equations to calculate rotor resistance and magnetizing inductance. If at least a part of the working cycle of the drive can be considered steady-state, a steady-state model of the machine can be used, which enables the usage of algebraic equations which are far less complex to solve [22].
Another area of computationally modest approaches is based on compensations of machine resistances based on temperature measurements. In [19], the authors used a temperature sensor mounted to the stator winding of the machine. Based on the temperature change, they calculated the values of both stator and rotor resistances assuming that the temperature of the whole machine is the same. In [41], the authors also equipped the rotor of a machine with a wireless temperature sensor besides the wired one attached to the stator winding. However, a temperature sensor is usually not standard equipment of a machine and has to be mounted additionally, which can be seen as a burden, especially in the case of a rotor temperature sensor.
In general, additional sensors are not favored due to cost and reliability reasons, which has led to a great advent of so-called sensor-less drives in recent years. However, a number of drives are still equipped with a speed sensor, which can be used to help with parameter identification. Such examples can be found in the manufacturing industry [42][43][44] or railway traction vehicles (locomotives, EMUs) [45], where induction machines are still the dominant type of electromechanical converter. Other sensors used in industrial or traction electric drives are current sensors on machine terminals and voltage sensors in the DC-links of power converters [44,46]. A method suitable for parameter identification of the drive with a microcontroller with limited resources should only use the quantities that are already measurable within the standard drive configuration.
The combination of a steady-state voltage model and information from a mechanical speed sensor was used in [12]. Stator voltages and currents, stator frequency, and mechanical frequency measured with a speed sensor were used to calculate rotor resistance to obtain the temperature change of the rotor of the induction machine in a battery electric vehicle drive. In [22], the same combination was used for the estimation of rotor magnetic flux and its angle.
This paper proposes the new method of the simultaneous on-line identification of rotor resistance and magnetizing inductance suitable for IFOC. This method has significantly lower computational demands than those mentioned in the introduction as it does not need excessive calculations, nor does it need to store large data in memory. It does not need any additional PI or I controllers in its structure and does not need to solve any differential equations. It is focused on drives with induction machines that are equipped with a mechanical speed sensor, a microcontroller with limited resources and where a part of their working cycle is in a steady-state mode. These could be existing drives with an older microcontroller system whose hardware cannot cope with computational or memoryintensive algorithms or new drives that need a speed sensor due to technical standards, or for high precision speed or position measurement, but otherwise are equipped with low-cost microcontrollers.

Drive with Field-Oriented Control
A block scheme of a variable-speed industrial drive with an induction machine is shown in Figure 1. As the voltage measurement of the switched voltage at the output of the inverter is complicated due to its rectangular nature with steep rising and falling edges, voltage is usually reconstructed from the switching pattern and the measured DC-link voltage that normally has far slower changes [17]. The nonlinearity of the inverter is also taken into account in most cases [12,46]. With the help of this procedure, stator voltage and frequency (or angular speed of the stator magnetic field) are obtained. Current sensors for measuring the machine stator current are placed on the inverter output, normally on two phases, while the third current is calculated from the other two [44]. A speed sensor is mounted on the shaft of the machine or on the gearbox, giving the information about the mechanical speed of the machine [8,45].  [42][43][44] or railway traction vehicles (locomotives, EMUs) [45], where induction machines are still the dominant type of electromechanical converter. Other sensors used in industrial or traction electric drives are current sensors on machine terminals and voltage sensors in the DClinks of power converters [44,46]. A method suitable for parameter identification of the drive with a microcontroller with limited resources should only use the quantities that are already measurable within the standard drive configuration. The combination of a steady-state voltage model and information from a mechanical speed sensor was used in [12]. Stator voltages and currents, stator frequency, and mechanical frequency measured with a speed sensor were used to calculate rotor resistance to obtain the temperature change of the rotor of the induction machine in a battery electric vehicle drive. In [22], the same combination was used for the estimation of rotor magnetic flux and its angle.
This paper proposes the new method of the simultaneous on-line identification of rotor resistance and magnetizing inductance suitable for IFOC. This method has significantly lower computational demands than those mentioned in the introduction as it does not need excessive calculations, nor does it need to store large data in memory. It does not need any additional PI or I controllers in its structure and does not need to solve any differential equations. It is focused on drives with induction machines that are equipped with a mechanical speed sensor, a microcontroller with limited resources and where a part of their working cycle is in a steady-state mode. These could be existing drives with an older microcontroller system whose hardware cannot cope with computational or memory-intensive algorithms or new drives that need a speed sensor due to technical standards, or for high precision speed or position measurement, but otherwise are equipped with low-cost microcontrollers.

Drive with Field-Oriented Control
A block scheme of a variable-speed industrial drive with an induction machine is shown in Figure 1. As the voltage measurement of the switched voltage at the output of the inverter is complicated due to its rectangular nature with steep rising and falling edges, voltage is usually reconstructed from the switching pattern and the measured DClink voltage that normally has far slower changes [17]. The nonlinearity of the inverter is also taken into account in most cases [12,46]. With the help of this procedure, stator voltage and frequency (or angular speed of the stator magnetic field) are obtained. Current sensors for measuring the machine stator current are placed on the inverter output, normally on two phases, while the third current is calculated from the other two [44]. A speed sensor is mounted on the shaft of the machine or on the gearbox, giving the information about the mechanical speed of the machine [8,45]. The controller of the drive is thus enabled to acquire stator current components Isd and Isq, stator voltage components Vsd and Vsq, stator angular frequency ωs and rotor angular frequency ωm. The proposed identification method requires no other quantities or sensors than these. The controller of the drive is thus enabled to acquire stator current components I sd and I sq , stator voltage components V sd and V sq , stator angular frequency ω s and rotor angular frequency ω m . The proposed identification method requires no other quantities or sensors than these.

Mathematical Model of the Machine
An induction machine mathematical model is usually represented by an equivalent circuit scheme. For the parameters mentioned in the introduction (stator resistance R s , rotor resistance R r , stator leakage inductance L σs , rotor leakage inductance L σr , and magnetizing inductance L m ), the so-called T-equivalent circuit is used ( Figure 2).

Mathematical Model of the Machine
An induction machine mathematical model is usually represented by an equivalent circuit scheme. For the parameters mentioned in the introduction (stator resistance Rs, rotor resistance Rr, stator leakage inductance Lσs, rotor leakage inductance Lσr, and magnetizing inductance Lm), the so-called T-equivalent circuit is used ( Figure 2). This circuit is derived based on the following vector equations that describe the induction machine in steady-state with harmonic supply: with the meaning of particular designations shown in Nomenclature. These vector equations can also be expressed in the form of their components in a dq reference frame, which is a common form of representation of vectors in microcontrollers. Capital letters are used for respective vector components to denote that they are steadystate values. The component notation of Equations (1) and md sd rd mq sq rq This circuit is derived based on the following vector equations that describe the induction machine in steady-state with harmonic supply: with the meaning of particular designations shown in Nomenclature. These vector equations can also be expressed in the form of their components in a dq reference frame, which is a common form of representation of vectors in microcontrollers. Capital letters are used for respective vector components to denote that they are steady-state values. The component notation of Equations (1)-(4) is: and

Solution of the Mathematical Model for Calculating Rotor Resistance R r and Magnetizing Inductance L m
Values of I sd , I sq , V sd , V sq , ω s , and ω r are supposed to be known in the microcontroller, as described in Section 2. Further, the values of parameters L σs , L σr , and R s are also assumed to be known, either provided by the manufacturer or by standard (laboratory) tests [13] or some off-line method. Aspects of using these three parameters are further discussed in Section 5. Assuming the steady-state of the machine, all those quantities are assumed to be constant during the calculation.
Back EMF components V id and V iq can be calculated from Equation (5) as all the other quantities are known: One of the vectors, e.g., the stator voltage, can be aligned with an axis of the dq reference frame, resulting in the particular component of such a vector being equal to zero. This zero component could possibly be omitted from the notation of the equations, but for generality it is assumed in this paper that the vectors can have any position.
For the clarity and simplification of the notation of the calculation, equivalent rotor resistance R req is defined as: Rotor current components I rd and I rq can be calculated by merging and rearranging Equation (7): From the equivalent circuit ( Figure 2), it is obvious that magnetizing current components I md and I mq can be expressed not only as stated in Equation (8), but also as: Substituting I md and I mq from Equation (12) to Equation (8) and I rd and I rq from Equation (11) to Equation (8), Equation (8) can be rewritten as: Converting I sd and I sq to a common denominator with the rest of the right side of Equation (13) forms: Now, the first equation in Equation (14) is multiplied with V id and the second equation in Equation (14) is multiplied with V iq , and the values are then added together to eliminate the left sides and obtain Equation (15): Rearranging the matching constituents of Equation (15) forms: As the left side of Equation (16) is equal to zero, the denominator can be removed. The expression in the nominator of Equation (16) can be rearranged as: which is a reduced quadratic equation of the type x 2 + px + q = 0 that can be solved by the common formula for finding its roots. For clarity, the denominator of p is defined as an inner power of the machine P i and the nominator of p as a squared magnitude of a back EMF vector |V i | 2 : For the motoric mode of the machine, the inner power P i is positive and therefore the positive sign of the square root in the quadratic formula is used: If the machine acts as a generator (regenerative braking), then the inner power is negative and therefore the negative sign is used, which is shown in Equation (20). The resulting value of R req will then be negative. Despite this, the resulting value of R r will be positive, as the value of slip s is also negative in the regenerative brake mode of the machine.
The final value of R r is then easily calculated from R req by multiplying it with slip s. With a known value of R req (or R r ), magnetizing inductance L m can be calculated from the second or preferably the first equation of Equation (12) as it does not contain a negative sign. The calculation is performed only with half of the components of the vectors V i and I m , which can have different signs based on actual vector positions. To prevent the negative resulting value of L m , the absolute value of the right side is added, forming the final equation for calculating L m as: As deriving the calculation is based on a vector description that can be in the scale of both maximum values (amplitudes) or RMS values, the calculation gives the same results both for maximum or RMS input values.

Algorithm for Calculation in Microcontroller
To summarize the previous deduction, a list of equations (Algorithm 1) that need to be performed by the microcontroller is presented in this chapter. Maximum or RMS values of input quantities can be used, considering particular implementation. First, back EMF components are calculated from stator current and voltage (Equations (a) and (b)). Then, the inner power (Equation (c)) and squared magnitude of the back EMF vector (Equation (d)) values are needed for the coefficient p (Equation (e)) of the reduced quadratic equation. The q coefficient is calculated in Equation (f) and the discriminant in Equation (g). Following the equivalent rotor resistance R req calculation (Equation (h)), its value is used to calculate the rotor current component (Equation (i)). In the denominator of Equation (i), the value of the q component is used once again as this part has the same value (see Equation (11)). The rotor current component is needed for the magnetizing current component calculation (Equation (j)). With the value of the magnetizing current component, L m can finally be calculated (Equation (k)). Rotor resistance R r could be calculated already after Equation (h), but is left till the end of the calculation in order to have the resulting values of parameters in the last two equations. So, in Equation (l), R req is multiplied by the value of the actual slip s (see Equation (10)) to obtain the R r . Algorithm 1. Proposed method for simultaneous estimation of R r and L m .

Experimental Results
The proposed method was verified experimentally on a laboratory drive with a squirrel cage induction machine with a rated power of 3.5 kW coupled to a 9.4 kW DC machine. The structure of this test bench is depicted in Figure 3.

Experimental Results
The proposed method was verified experimentally on a laboratory drive with a squirrel cage induction machine with a rated power of 3.5 kW coupled to a 9.4 kW DC machine. The structure of this test bench is depicted in Figure 3.  Table 1 and parameters obtained by standard tests (laboratory measurements) according to [13] are in Table 2.  Table 2. Values of parameters of tested machine used for calculations obtained by standard tests specified in [13]. Rated values are in Table 1 and parameters obtained by standard tests (laboratory measurements) according to [13] are in Table 2.  Table 2. Values of parameters of tested machine used for calculations obtained by standard tests specified in [13].

Parameter Value
Stator resistance R s 1.11 Ω Stator leakage inductance L σs 8.25 mH Rotor leakage inductance L σr 8.25 mH Other components of the test bench (see Figure 4) were a two-level three-phase voltage type DC/AC IGBT inverter, voltage transducer LEM LV25-P for DC-link voltage measurement and current transducers LEM LF205-S for stator current measurement. Stator voltage and its angular frequency were obtained as described in Section 2. Rotor mechanical speed was measured with a LARM IRC 300/1024 optical incremental encoder with 1024 impulses per revolution. The heart of the control system was a TI TM4C123G-H6PM microcontroller with 12-bit successive approximation AD converters. For torque measurements, a torque-measuring shaft KTR DATAFLEX 22/100 was used. As

Warming Test
To observe how the method can track the two parameters in the course of time, a 60 min long warming test with a constant load (6 N·m) was performed. The change in rotor resistance Rr and frame warming is depicted in Figure 5. Data for each point were averaged from 12 sets of samples of input quantities, and the time span between these samples was 10 s. The temperature of the machine frame was measured with a PT100 sensor. The initial temperature was 32 °C and the final temperature was 80 °C, resulting in a warming of 48 °C. Rr rose from 1.05 Ω to 1.28 Ω, which is an increase of 22%. Considering the well-known equation for the temperature dependence The heart of the control system was a TI TM4C123G-H6PM microcontroller with 12-bit successive approximation AD converters. For torque measurements, a torque-measuring shaft KTR DATAFLEX 22/100 was used.
As stated at the end of Section 3.1, maximum (amplitudes) or RMS input values can be chosen for input. The following experiments use maximum values in all cases. The voltage vector is aligned with the q axis, so the d component (V sd ) is zero in all presented cases.

Warming Test
To observe how the method can track the two parameters in the course of time, a 60 min long warming test with a constant load (6 N·m) was performed. The change in rotor resistance R r and frame warming is depicted in Figure 5.

Warming Test
To observe how the method can track the two parameters in the course of time, a 60 min long warming test with a constant load (6 N·m) was performed. The change in rotor resistance Rr and frame warming is depicted in Figure 5. Data for each point were averaged from 12 sets of samples of input quantities, and the time span between these samples was 10 s. The temperature of the machine frame was measured with a PT100 sensor. The initial temperature was 32 °C and the final temperature was 80 °C, resulting in a warming of 48 °C. Rr rose from 1.05 Ω to 1.28 Ω, which is an increase of 22%. Considering the well-known equation for the temperature dependence of metals (Equation (22)), and assuming the linear behavior of this dependence, this change in resistance estimated by the proposed method is in good accordance with it. 2 20 20 20 20 (1 α Δ ) Data for each point were averaged from 12 sets of samples of input quantities, and the time span between these samples was 10 s. The temperature of the machine frame was measured with a PT100 sensor. The initial temperature was 32 • C and the final temperature was 80 • C, resulting in a warming of 48 • C. R r rose from 1.05 Ω to 1.28 Ω, which is an increase of 22%. Considering the well-known equation for the temperature dependence of metals (Equation (22)), and assuming the linear behavior of this dependence, this change in resistance estimated by the proposed method is in good accordance with it.
where ∆ϑ = 48 • C (warming during our test) and α 20 = 0.0038 K −1 is the temperature coefficient for copper (material of the cage of the tested machine). The result of Equation (22) means a rise of 18% from the initial value R 20 compared to 22% obtained by the proposed method. It is correct that the value obtained by this method is higher than that obtained by the temperature change of the stator. The temperature of the rotor is normally higher than that of the stator and this difference increases the longer the machine runs [12]. This comparison is, however, only rough as, for example, the absolute temperature of the rotor is unknown to us. In Figure 6, there is no visible influence of temperature upon magnetizing inductance L m , which corresponds with theoretical assumptions (see Section 1). Both stator voltage and current were kept almost constant during this warming test (see Figure 7)-stator current was supplied from a converter with a current loop. Because of this, there is a small variance in stator current. We used this variance to show the dependence of L m on the magnetizing component of the stator current (I sd ) in Figure 8, where it is depicted in better scale. This dependence under different loads and frequencies is also shown later in Figure 11. Such a dependence is presented in [11,47], so it can be regarded as another proof of the results obtained by the proposed method. and current were kept almost constant during this warming test (see Figure 7)-stator current was supplied from a converter with a current loop. Because of this, there is a small variance in stator current. We used this variance to show the dependence of Lm on the magnetizing component of the stator current (Isd) in Figure 8, where it is depicted in better scale. This dependence under different loads and frequencies is also shown later in Figure  11. Such a dependence is presented in [11,47], so it can be regarded as another proof of the results obtained by the proposed method.   current was supplied from a converter with a current loop. Because of this, there is a small variance in stator current. We used this variance to show the dependence of Lm on the magnetizing component of the stator current (Isd) in Figure 8, where it is depicted in better scale. This dependence under different loads and frequencies is also shown later in Figure  11. Such a dependence is presented in [11,47], so it can be regarded as another proof of the results obtained by the proposed method.

Measurements with Different Loads and Stator Frequencies
Another set of measurements was conducted with different load torques and four different stator frequencies, namely 20, 30, 40, and 50 Hz. Each measurement is represented by five working points with different loads. The temperature difference among measurement points did not exceed 5 °C, so its influence should not be significant in this case. Each measurement point was averaged from five samples with a time span of 10 s.
Stator voltage was set according to the V/f (voltage to frequency) ratio. For measurements at stator frequencies 20, 30, and 40 Hz, this ratio was 6.5. At 50 Hz, the field weakening area was reached; therefore, the V/f ratio was only 5.6. The voltage vector is again aligned with the q axis, so the d component (Vsd) is zero in all presented cases.
The measured data and estimated values of Rr and Lm are presented in the tables

Measurements with Different Loads and Stator Frequencies
Another set of measurements was conducted with different load torques and four different stator frequencies, namely 20, 30, 40, and 50 Hz. Each measurement is represented by five working points with different loads. The temperature difference among measurement points did not exceed 5 • C, so its influence should not be significant in this case. Each measurement point was averaged from five samples with a time span of 10 s.
Stator voltage was set according to the V/f (voltage to frequency) ratio. For measurements at stator frequencies 20, 30, and 40 Hz, this ratio was 6.5. At 50 Hz, the field weakening area was reached; therefore, the V/f ratio was only 5.6. The voltage vector is again aligned with the q axis, so the d component (V sd ) is zero in all presented cases.
The measured data and estimated values of R r and L m are presented in the tables below (Tables 3-6). The tables show stator current components I sd and I sq , measured load torque T load , and rotor angular frequency ω m . The representation of speed by angular frequency was chosen because angular frequency is used in the equations in Section 3. For the calculation of rotations per minute of the shaft, note that the machine has six poles (see Table 1). Table 3. Measurement at f s = 20 Hz (ω s = 125.66 rad·s −1 ) and V sq = 130 V.  Table 4. Measurement at f s = 30 Hz (ω s = 188.50 rad·s −1 ) and V sq = 195 V.  Table 5. Measurement at f s = 40 Hz (ω s = 251.33 rad·s −1 ) and V sq = 260 V.  Table 6. Measurement at f s = 50 Hz (ω s = 314.16 rad·s −1 ) and V sq = 280 V. It is generally difficult to compare the estimated parameter values with some references. Laboratory standard tests [13], for example, use the line frequency (50 Hz in our case) for locked rotor tests, meaning that slip frequency is equal to 50 Hz. On the other hand, during the operation of the drive with scalar control or IFOC, slip frequency is low (in order of hertz). Magnetizing inductance has a similar problem with the comparison of its estimated value to the value obtained by standard tests: the proposed identification method determines its value for a particular working point, while the standard test value is obtained from no-load measurements. So far, no on-line method has been generally or officially standardized to become a source of reference measurement that other methods could be compared to. As mentioned in the introduction, only R s is directly measurable in the case of an induction machine with a squirrel cage rotor. Therefore, we compare our results with trends of parameter changes published elsewhere [11,19,28,47] and by a verification scheme described below in this chapter.

Measured Input Quantities
To demonstrate that the derivation of the equations of the proposed method is correct, the verification scheme depicted in Figure 9 was proposed. The idea behind this verification scheme is to calculate some measurable quantities from the mathematical model of the machine with the help of estimated parameters and then compare measured and calculated values. In the presented case, stator current components were selected. With the proposed on-line identification method, Rr and Lm were estimated and then their values were used in a mathematical model of the machine. In this model, represented again by a T-equivalent circuit (Figure 2), the knowledge of all the equivalent circuit parameters and some of the measurable quantities, namely stator voltage components Vsd and Vsq, stator angular frequency ωs, and rotor angular frequency ωm, was used. From these, the stator current components Isdc and Isqc were calculated, with the The idea behind this verification scheme is to calculate some measurable quantities from the mathematical model of the machine with the help of estimated parameters and then compare measured and calculated values. In the presented case, stator current components were selected. With the proposed on-line identification method, R r and L m were estimated and then their values were used in a mathematical model of the machine. In this model, represented again by a T-equivalent circuit (Figure 2), the knowledge of all the equivalent circuit parameters and some of the measurable quantities, namely stator voltage components V sd and V sq , stator angular frequency ω s , and rotor angular frequency ω m , was used. From these, the stator current components I sdc and I sqc were calculated, with the letter "c" in their subscript standing for "calculated".
The calculated values of I sdc and I sqc are exactly the same as those that can be measured (I sd , I sq ) because this T-equivalent circuit is not capable of introducing some other parasitic effects. Despite this, it demonstrates that the values of R r and L m have values appropriate for the T-equivalent circuit representation of the machine.
A number of parameter characteristics can be seen from the presented measurements. Figure 10 shows the dependence of R r on slip frequency f r , which corresponds with observations published elsewhere [11,19,28,47]. With the rising load torque T load , the mechanical speed decreased, which resulted in an increase in slip and therefore the slip frequency f r . There is also a visible influence of stator frequency f s . R r is higher mainly due to the skin effect in rotor bars, which is more significant the higher the frequencies are.

Discussion
Due to the character of the results, their interpretation has been carried out in the previous section (Section 4) along with their presentation.
An estimation of the proper value of Rr by the proposed method depends on the accuracy of the speed measurement. If the drive is not equipped with a sensor with a suitable number of pulses per revolution, approaches such as longer measurement or the measurement of the period of the pulse can be deployed. Prolonging the speed measurement should not be a problem as the drive should operate in steady-state mode during measurement. Parameters Lσs, Lσr, and Rs are assumed to be constant. This simplification can be There is again a visible dependence of L m on the magnetizing component of stator current component I sd (described in Section 4.1), as shown in Figure 11. In this figure, data from all the measurements with different loads and different frequencies are grouped and the trend corresponds with the results published elsewhere, e.g., in [11,47].

Discussion
Due to the character of the results, their interpretation has been carried out in the previous section (Section 4) along with their presentation.
An estimation of the proper value of Rr by the proposed method depends on the accuracy of the speed measurement. If the drive is not equipped with a sensor with a suitable Figure 11. L m as a function of stator current component I sd from measurements with different loads and different frequencies.

Discussion
Due to the character of the results, their interpretation has been carried out in the previous section (Section 4) along with their presentation.
An estimation of the proper value of R r by the proposed method depends on the accuracy of the speed measurement. If the drive is not equipped with a sensor with a suitable number of pulses per revolution, approaches such as longer measurement or the measurement of the period of the pulse can be deployed. Prolonging the speed measurement should not be a problem as the drive should operate in steady-state mode during measurement. Parameters L σs , L σr , and R s are assumed to be constant. This simplification can be performed without a significant loss of accuracy in the case of L σs and L σr as they represent the part of magnetic flux that flows mainly through air. Thus, it is not influenced by a change in temperature and the effect of saturation can also be considered negligible. R s is, however, influenced by the temperature, but its influence on the calculation is not so significant in most of the operation range. The voltage drop on L σs is higher than that on R s , and this difference rises with the rising stator frequency and torque-producing stator current component I sq . Moreover, the temperature of the rotor is normally higher than that of the stator [12], which means that the change in R s is smaller than the change in R r . Frequency should not have any significant effect on R s as the stator winding is normally manufactured as parallel filaments to prevent the skin effect.
The proposed method is designated for machines that can be described by a model in the form of a T-equivalent circuit, as presented in Figure 2, which is the majority of induction machines used in industry or traction. An example of a machine type that could not be described with such a model can be a machine with a double rotor cage that would need an additional parallel rotor branch in its equivalent circuit. When the machine conforms to the stated criterion, then the results should be correct regardless of the construction details of the machine. Different construction details would be, however, visible on the results. For example, the change in rotor resistance due to the skin effect caused by slip frequency will depend on the shape of the cross section of the rotor bars. The deeper the bars, the more visible the effective rotor resistance rise due to the skin effect. On the other hand, the wound rotor will not show any significant influence of the skin effect on the rotor resistance as the cross section of the particular winding filaments is small compared to the bars of the squirrel cage.
The shape of the slots (both on the stator and rotor) will have a further effect on the course of magnetizing inductance. Slots whose openings widen towards the air gap (e.g., made for being closed with wedges) will cause the machine to be more sensitive to saturation as this makes the magnetic circuit path narrower compared to closed slots.
The proposed method is thus comparatively universal regardless of the construction features of a particular machine, unless it can be described by a T-equivalent circuit.
Both R r and L m can be regarded as slowly varying during steady-state operation. In this mode, R r is mostly influenced by temperature change, which rises slowly, as can be seen from Figure 5. The flux, which can influence the value of L m by saturating the magnetic circuit, also does not change quickly during the normal operation of IFOC. Therefore, the method is enough to be executed within minutes during steady-state operation. In the case of a change in a working point of the drive, it should be executed after this change as magnetizing current or slip change can result in a sharp change in R r or L m . The averaging of multiple measurements of input quantities is advised to prevent the variation of results.

Conclusions
This paper presents a new, computationally modest on-line method for the simultaneous identification of rotor resistance R r and magnetizing inductance L m of induction machines. These two parameters are crucial for IFOC, which is a popular control strategy in industrial or traction drives. However, their values change during operation due to the temperature change, frequency, or saturation of the magnetic circuit. Tracking their values can thus improve the performance of the control strategy.
The proposed method utilizes a steady-state voltage model of the machine, which enables a significant simplification of necessary equations. The required quantities are obtained only through equipment that is already present in common electric drives, such as stator current measurement, DC-link voltage measurement, and a speed sensor. It is suitable for drives with resource-constrained microcontrollers that have at least a part of their working cycle in a steady-state mode. The proposed method is designated for machines that can be described by a model in the form of a T-equivalent circuit, which is the majority of induction machines used in industry or traction. When the machine conforms to the stated criterion, then the results should be correct regardless of the construction details of the machine. From the point of view of the overall construction of the drive, a part that can have a significant influence on the results is the speed sensor. This is because the slip of the machine is used for the calculation of R r . If a sensor with a low number of pulses per revolution is deployed, the measurement of speed can be improved by uncomplicated techniques such as longer measurement or the measurement of the period of the pulse, which will not significantly increase the complexity of the method.
Deriving the method is presented in detail, including a final algorithm of equations for calculation in the microcontroller. Experimental validation was carried out by measurements with different loads at different frequencies and a warming test. Difficulties of the comparison of the resulting values of R r and L m to some reference values are presented. To validate the results, two other approaches were therefore used. Firstly, the resulting values of parameters were used in a mathematical model of the machine to calculate stator current components. These calculated stator current components were compared with the measured ones, and both of these values corresponded with each other. This confirms that the resulting values of R r and L m are suitable for a working mathematical model of the machine and also confirms that the algorithm of the method is derived correctly. Secondly, the results were further discussed and compared with other publications to prove that the proposed method gives relevant results.
Author Contributions: This article presents the collective work of both authors (T.K. and P.K.). All authors have read and agreed to the published version of the manuscript.
Funding: This research was funded by the Czech CRRC Science and Technology Development Ltd. in Prague, Czech Republic.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to company restrictions.

Conflicts of Interest:
The authors declare no conflict of interest. Nomenclature R s stator resistance R r rotor resistance L σs stator leakage inductance L σr rotor leakage inductance L m magnetizing inductance V s , V sd , V sq stator voltage vector and its components in dq reference frame I s , I sd , I sq stator current vector and its components in dq reference frame I sdc , I sqc calculated values of stator current components in dq reference frame V i , V id , V iq back EMF vector and its components in dq reference frame I m , I md , I mq magnetizing current vector and its components in dq reference frame I r , I rd , I rq rotor current vector and its components in dq reference frame ω s stator angular frequency ω m rotor angular frequency j imaginary unit s slip defined as s = ω s −ω m ω s R req equivalent rotor resistance defined as R req = R r s P i inner machine power p, q coefficients of the quadratic equation R ϑ2 value of resistance at elevated temperature R 20 value of resistance at 20 • C α 20 temperature coefficient ∆ϑ temperature difference ϑ frame machine frame temperature f s stator frequency f r slip frequency T load load torque