사각형입니다.

https://doi.org/10.6113/JPE.2018.18.2.492

ISSN(Print): 1598-2092 / ISSN(Online): 2093-4718



A Novel Method for the Identification of the Rotor Resistance and Mutual Inductance of Induction Motors Based on MRAC and RLS Estimation


Gwon-Jae Jo* and Jong-Woo Choi


†,*Department of Electrical Engineering, Kyungpook National University, Daegu, Korea



Abstract

In the rotor-flux oriented control used in induction motors, the electrical parameters of the motors should be identified. Among these parameters, the mutual inductance and rotor resistance should be accurately tuned for better operations. However, they are more difficult to identify than the stator resistance and stator transient inductance. The rotor resistance and mutual inductance can change in operations due to flux saturation and heat generation. When detuning of these parameters occurs, the performance of the control is degenerated. In this paper, a novel method for the concurrent identification of the two parameters is proposed based on recursive least square estimation and model reference adaptive control.


Key words: Induction motor, Model reference adaptive control, Mutual inductance, Parameter identification, Recursive least square estimation, Rotor resistance


Manuscript received Jun. 16, 2017; accepted Oct. 6, 2017

Recommended for publication by Associate Editor Kwang-Woon Lee.

Corresponding Author: cjw@knu.ac.kr Tel: +82-53-950-5515, Kyungpook National University

*Dept. of Electrical Engineering, Kyungpook National University, Korea



NOMENCLATURE

Rs

Stator resistance

Rr

Rotor resistance

Lm

Mutual inductance

Ls

Stator self-inductance

Lr

Rotor self-inductance

Lls

Stator leakage inductance

Llr

Rotor leakage inductance

σLs

Stator transient inductance

Tr

Rotor time constant

s, r

Superscripts: denote the stationary and rotor reference frames

s, r

Subscripts: denote the stator and rotor

^

Denotes estimated quantities

s

Derivative operator

Bold form x

Vector xdjxq



Ⅰ. INTRODUCTION

Induction motors are widely used in industrial and home appliances. For the speed control of induction motors, rotor-flux oriented control is usually utilized due to its simplicity and its ability to control flux and torque components separately. In the rotor-flux control algorithm, the precise tuning of Lm and Rr is very important for efficient control of the currents of both the flux and torque components. Moreover, these parameters are essential in flux observers. However, identifying them is more difficult than identifying stator parameters such as Rs and σLs. This is attributed to the fact that Lm varies due to nonlinear characteristics in the B-H curve. It is also due to the fact that Rr is more affected by thermal temperature in operation than Rs since the rotor is positioned inside the motor.

Several identification methods for Rr and Tr have been reported. A method by which Rr can be estimated using the relationship between the flux and the developed torque was proposed in [1]. A method, which helps identify Tr by solving simultaneous equations for the rotor voltage and the torque, was proposed in [2]. Methods where Rr is identified by adopting adaptive observers were reported in [3], [4]. A method using a particle swarm optimization algorithm was reported in [5]. Methods for the identification of Rr based on model reference adaptive control (MRAC) were proposed in [6], [7], where the variables observed in the MRAC are the voltage [6] and the instantaneous reactive power [7]. A method that compares fluxes estimated by flux observers, was reported in [8]. The rotor time constant Tr can be identified from a model reference adaptive system speed estimator in sensorless systems [9]. Using the dynamics of a synchronous reference frame current controller, the identification of Tr was studied in [10].

In addition to the identification of single rotor parameters, techniques for the identification of various other parameters were reported in [11]-[18]. Parameters have been identified offline based on fitting the experimental data and numerical analysis [11]. Voltage equations of the stator and rotor have been solved directly by integration to identify Lm and Rr [12], [13]. Based on the MRAC technique, currents have been compared for identification in [14]. Furthermore, parameters have been identified using recursive least square (RLS) estimation [15], [16]. Methods using sliding mode observers and artificial neural networks have been proposed [17], [18] for parameter identification.

In this paper, a novel method for the identification of Lm and Rr based on MRAC is proposed. In the MRAC algorithm, the mechanism is to compare the flux estimated by the flux observers designed from current and voltage models. Furthermore, the differences in fluxes estimated by the two flux observers are analyzed. The proposed method can operate at a sufficiently high frequency in order to obtain the complete reference model in the MRAC system. Thus, it is difficult to apply the proposed method at a low frequency including standstill. However, most of the papers published so far have reported the identification of single parameter, whereas this paper identifies Lm and Rr concurrently by separating the two parameters with respect to the equation. Thus, this paper has advantage since there are no dependencies on each other for Lm and Rr in the proposed method. Additionally, the proposed method can guarantee effective behavior in both the steady-state and transient-state because any assumption for the steady state in the algorithm is not involved. RLS estimation is utilized for the realization of the parameters identification. The feasibility of the proposed method has been verified by simulation and experimental results under diverse conditions.



Ⅱ. GOVERNING EQUATIONS

The flux observers utilized in this paper are based on the governing equations for the model of an induction motor. The governing equations on an arbitrary rotating axis, given by ω in squirrel cage-type induction motors consist of the voltage equations and linkage flux equations in the stator and rotor. The voltage equations are expressed as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c127e.bmp
원본 그림의 크기: 가로 624pixel, 세로 86pixel   (1)

그림입니다.
원본 그림의 이름: CLP00001b0c0001.bmp
원본 그림의 크기: 가로 798pixel, 세로 77pixel     (2)

The linkage flux equations are expressed as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c0002.bmp
원본 그림의 크기: 가로 436pixel, 세로 81pixel           (3)

그림입니다.
원본 그림의 이름: CLP00001b0c0003.bmp
원본 그림의 크기: 가로 439pixel, 세로 77pixel           (4)



Ⅲ. FLUX OBSERVERS


A. Voltage Model

The voltage model flux observer is designed from the stator voltage equation in [19]. The equation for the estimation of the rotor flux in the voltage model can be expressed using (1), (3), and (4) by defining these equations on the stationary reference frame (ω = 0). Thus, the rotor flux vector can be estimated from the stator voltage and current as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c0004.bmp
원본 그림의 크기: 가로 802pixel, 세로 182pixel     (5)

However, the voltage model flux observer cannot be implemented because the pure integration in (5) can diverge due to a dc offset voltage. Thus, the closed-loop Gopinath style flux observer in [20] has been employed. This flux observer model operates the same as the voltage model flux observer in a sufficiently high frequency region with solving the problem with the dc offset voltage.


B. Current Model

The current model is a typical flux observer in the rotor-flux oriented control [21]. This model is based on the rotor equation. A differential equation on the rotor reference frame (ω = ωr) can be expressed from (2) and (4) as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c0005.bmp
원본 그림의 크기: 가로 358pixel, 세로 188pixel     (6)

Equation (6) yields the low-pass filter (LPF) form constructed from only the rotor flux and stator current. The stator current only can be used to estimate the rotor flux. The rotor flux vector on the stationary reference frame can be calculated from 그림입니다.
원본 그림의 이름: CLP00001b0c0006.bmp
원본 그림의 크기: 가로 62pixel, 세로 79pixel in (6) as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c0007.bmp
원본 그림의 크기: 가로 289pixel, 세로 92pixel        (7)

where 그림입니다.
원본 그림의 이름: CLP00001b0c0008.bmp
원본 그림의 크기: 가로 238pixel, 세로 150pixel



Ⅳ. MRAC

In this paper, parameter identification is based on the principle of MRAC. Fig. 1 shows a MRAC system that consists of a reference model, an adjustable model, and an adaptation mechanism. The reference model is independent of the parameters, whereas the adjustable model is dependent on the parameters. The adaptation mechanism generates the estimated parameters by comparing the variables (xref and xadj) from the reference and adjustable models.


그림입니다.
원본 그림의 이름: image9.png
원본 그림의 크기: 가로 507pixel, 세로 320pixel

Fig. 1. MRAC system.


The current model flux observer is dependent on Lm and Rr according to (6). On the other hand, the voltage model flux observer is independent of these parameters. Firstly, the estimation equation of the voltage model (5) does not include Rr. Secondly, the ratio 그림입니다.
원본 그림의 이름: CLP00001b0c0009.bmp
원본 그림의 크기: 가로 167pixel, 세로 80pixel is nearly unity due to a small leakage ratio. This is shown by the following equation:

그림입니다.
원본 그림의 이름: CLP00001b0c000a.bmp
원본 그림의 크기: 가로 612pixel, 세로 259pixel    (8)

Therefore, the detuning of Lm barely affects the estimated flux in the voltage model. From the principle of MRAC, the voltage and current models can be chosen as reference and adjustable models in the identification algorithm for Lm and Rr. In the detuned case of Lm and Rr, the rotor flux estimated by the current model has errors. On the other hand, the rotor flux can be obtained quite accurately from the voltage model. Therefore, correct values of Lm and Rr can be obtained from the adaptation mechanism wherein the fluxes estimated by the two flux observers can be compared. As the parameters are calibrated by the adaptation mechanism, the adjustable model will follow the reference model.



Ⅴ. PARAMETER IDENTIFICATION


A. Expression of the Flux Difference

The estimation equation given by (6) can be rearranged if the current model is used as an adjustable model, as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c000b.bmp
원본 그림의 크기: 가로 682pixel, 세로 190pixel          (9)

where 그림입니다.
원본 그림의 이름: CLP00001b0c000c.bmp
원본 그림의 크기: 가로 115pixel, 세로 77pixel denotes 그림입니다.
원본 그림의 이름: CLP00001b0c000d.bmp
원본 그림의 크기: 가로 64pixel, 세로 80pixel from the current model. Equivalently, if the voltage model is used as a reference model, (9) is rewritten as (10) given by:

그림입니다.
원본 그림의 이름: CLP00001b0c000e.bmp
원본 그림의 크기: 가로 687pixel, 세로 167pixel          (10)

where 그림입니다.
원본 그림의 이름: CLP00001b0c000f.bmp
원본 그림의 크기: 가로 115pixel, 세로 72pixel denotes 그림입니다.
원본 그림의 이름: CLP00001b0c000d.bmp
원본 그림의 크기: 가로 64pixel, 세로 80pixel from the voltage model. Equation (10) consist of actual parameters because this model is robust to the detuning of Lm and Rr. Fig. 2 shows the proposed MRAC system for parameter identification. When (9) is subtracted from (10), an equation with a flux difference 그림입니다.
원본 그림의 이름: CLP000025440026.bmp
원본 그림의 크기: 가로 557pixel, 세로 85pixel is obtained as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c0010.bmp
원본 그림의 크기: 가로 1244pixel, 세로 201pixel      (11)


그림입니다.
원본 그림의 이름: image13.png
원본 그림의 크기: 가로 721pixel, 세로 305pixel

Fig. 2. MRAC system for the proposed parameter identification.


Equation (11) can be rearranged to obtain the following equation:

그림입니다.
원본 그림의 이름: CLP00001b0c0011.bmp
원본 그림의 크기: 가로 1288pixel, 세로 204pixel     (12)

The ratio 그림입니다.
원본 그림의 이름: CLP00001b0c0015.bmp
원본 그림의 크기: 가로 409pixel, 세로 81pixel in (12) is approximately unity due to a small leakage ratio as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c0012.bmp
원본 그림의 크기: 가로 861pixel, 세로 270pixel         (13)

Equations (12) and (13) lead to the following rearranged equation:

그림입니다.
원본 그림의 이름: CLP00001b0c0013.bmp
원본 그림의 크기: 가로 1182pixel, 세로 201pixel        (14)

Equation (14) implies that the flux difference can be regarded as a linear equation with the two filtered signals, and the two coefficients given by:

그림입니다.
원본 그림의 이름: CLP00001b0c0014.bmp
원본 그림의 크기: 가로 764pixel, 세로 190pixel       (15)


B. RLS Estimation

For solving (14), RLS estimation, which is an algorithm for estimating unknown coefficients in solving regression problems, is employed. In this paper, RLS estimation with a forgetting factor is utilized in order to enhance the rate of convergence [22]. A regression equation is expressed by (16), where y is the output, X is the regression vector, and A is the coefficient vector:

그림입니다.
원본 그림의 이름: CLP00001b0c0016.bmp
원본 그림의 크기: 가로 240pixel, 세로 81pixel          (16)

In this system, the unknown coefficients are estimated by (17) as follows, where K is the gain vector and e is the error:

그림입니다.
원본 그림의 이름: CLP00001b0c0017.bmp
원본 그림의 크기: 가로 703pixel, 세로 91pixel         (17)

Here, e is defined by (18) and K is calculated from (19), where P is the covariance matrix and λ is the forgetting factor.

그림입니다.
원본 그림의 이름: CLP00001b0c0018.bmp
원본 그림의 크기: 가로 680pixel, 세로 96pixel           (18)

그림입니다.
원본 그림의 이름: CLP00001b0c0019.bmp
원본 그림의 크기: 가로 761pixel, 세로 160pixel       (19)

Equation (17) implies that the estimated coefficient is compensated by the gain times the error so that the error becomes zero. λ plays the role of putting more weight on current state than the previous state, which promotes the convergence of 그림입니다.
원본 그림의 이름: CLP00001b0c001a.bmp
원본 그림의 크기: 가로 47pixel, 세로 67pixel. Lastly, P is renewed as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c001b.bmp
원본 그림의 크기: 가로 988pixel, 세로 141pixel       (20)

On application of the RLS algorithm to (14), the output and regression vector in the equation become:

그림입니다.
원본 그림의 이름: CLP00001b0c001c.bmp
원본 그림의 크기: 가로 203pixel, 세로 71pixel   (21)

그림입니다.
원본 그림의 이름: CLP00001b0c001d.bmp
원본 그림의 크기: 가로 809pixel, 세로 206pixel     (22)

The output and regression vector can be calculated from the fluxes estimated by the flux observers and the filters which have a known time constant 그림입니다.
원본 그림의 이름: CLP00001b0c001e.bmp
원본 그림의 크기: 가로 60pixel, 세로 80pixel. The process for the calculation of the output and regression vector is depicted in Fig. 3. In this way, the coefficients in (15) are estimated by RLS. Note that the parameters (Lm and Rr) to be identified are separate in terms of the equation. Therefore, the detuning of Rr does not affect the identification of Lm and vice versa.


그림입니다.
원본 그림의 이름: image22.png
원본 그림의 크기: 가로 428pixel, 세로 303pixel

Fig. 3. Calculation of the regression vector and output.


C. Coefficient Controller

The compensation components of the respective parameters are calculated from the coefficients in (15) for tuning them. These compensation components are obtained by regulating the estimated coefficients (â1 and â2). When â1 and â2 are reduced to zero, 그림입니다.
원본 그림의 이름: CLP00001b0c001f.bmp
원본 그림의 크기: 가로 239pixel, 세로 80pixel and 그림입니다.
원본 그림의 이름: CLP00001b0c0020.bmp
원본 그림의 크기: 가로 219pixel, 세로 81pixel are implicated from (15). Based on MRAC, the flux difference becomes zero according to (14). A proportional-integral (PI) controller is employed for the regulation of â1 and â2 as shown Fig. 4. The outputs of the controller (ΔLm and ΔRr) are used as compensation components of the parameters as follows:

그림입니다.
원본 그림의 이름: CLP00001b0c0021.bmp
원본 그림의 크기: 가로 430pixel, 세로 92pixel   (23)

그림입니다.
원본 그림의 이름: CLP00001b0c0022.bmp
원본 그림의 크기: 가로 407pixel, 세로 84pixel   (24)


그림입니다.
원본 그림의 이름: image31.png
원본 그림의 크기: 가로 661pixel, 세로 211pixel

Fig. 4. Coefficient controller.


where ΔLm and ΔRr are the compensation components of Lm and Rr, and 그림입니다.
원본 그림의 이름: CLP00001b0c0023.bmp
원본 그림의 크기: 가로 155pixel, 세로 79pixel and 그림입니다.
원본 그림의 이름: CLP00001b0c0024.bmp
원본 그림의 크기: 가로 147pixel, 세로 78pixel are their initial setting values. In addition, the final outputs ΔLm and ΔRr are filtered by a LPF to prevent noise.

The process for the proposed parameter identification is shown in Fig. 5. A block diagram of a rotor-flux oriented control system with the proposed parameter identification algorithm in an induction motor is shown in Fig. 6.


그림입니다.
원본 그림의 이름: image37.png
원본 그림의 크기: 가로 1468pixel, 세로 165pixel

Fig. 5. Process of the proposed parameter identification.


그림입니다.
원본 그림의 이름: image38.png
원본 그림의 크기: 가로 1264pixel, 세로 506pixel

Fig. 6. Rotor-flux oriented control with the proposed parameter identification in an induction motor.



Ⅵ. SIMULATION RESULTS

The specifications and parameters of the utilized induction motor are specified in Table I. All of the controllers in the simulation are executed digitally, where the control periods of the current controller and speed controller are 100 μs and 400 μs, respectively. The control period of the parameter identification is synchronized with the speed controller. The test speeds are 300 (low-speed), 600 (mid-speed), and 1200 rpm (high-speed). Furthermore, various loads (30%, 50% and 80% of the rating) are applied to the induction motor. The forgetting factor in the RLS is set to 0.99. The initial value of Lm is 1.5 times (206 mH) the actual value, and Rr is detuned to be 0.5 times (0.37 Ω) the actual value. First, a situation with perfectly tuned stator parameters (σLs and Rs) is tested. Second, the detuned cases for the stator parameters are simulated, where the voltage model flux observer generated an erroneous flux due to the detuning. The parameter identification algorithm starts in the steady state when the speed becomes consistent with the applied load torque. Furthermore, Lm and Rr, estimated by the algorithm, are updated to all of the controllers and the algorithm in real time, as shown in Fig. 6.


TABLE I SPECIFICATIONS AND PARAMETERS OF AN INDUCTION MOTOR

Specifications and Parameters

Quantities

Rated power

1.5 kW

Pole

4 pole

Rated speed

1,455 rpm

Rated torque

9.2 N·m

Rs

1.67 Ω

Rr

0.73 Ω

Lm

137 mH

Lls,Llr

6.5 mH


In fact, the vector form (14), which involves both the d-axis and q-axis, can be analyzed on one of the two axes. First, the identification with respect to only the d-axis is tested. Fig. 7(a) shows the responses of the actual coefficients (a1) and the estimated coefficients (â1) for the d-axis. â1 followed a1 from its initial value since the start of identification. Although the coefficient converges to zero asymptotically in the steady state, â1 has oscillation in the transient state due to the fact that λref is sinusoidal signal with a slip frequency. For generating the parameter compensation components, it is appropriate that â1 and â2 to be regulated by the PI controller should follow a1 and a2 smoothly without oscillations in the transient state. To achieve this, the average of two coefficients on the d-axis and q-axis is considered. It is confirmed from Fig. 7(b) that the oscillation of the coefficient in the transient state is removed after averaging the outputs on the two axes.


Fig. 7. Responses of coefficients. (a) On the d-axis. (b) Upon averaging.

그림입니다.
원본 그림의 이름: CLP00001b0c002f.bmp
원본 그림의 크기: 가로 1745pixel, 세로 864pixel

(a)

 

그림입니다.
원본 그림의 이름: CLP00001b0c0030.bmp
원본 그림의 크기: 가로 1750pixel, 세로 863pixel

(b)


Fig. 8(a) shows the result of identification for Lm and Rr under given conditions. Upon the implementation of the proposed algorithm, both of the parameters converge to their actual values within 5 s. Although the approximation in (13) is applied, Rr is accurately estimated since Lm is simultaneously calibrated to yield the actual value. When the parameters are calibrated, both of d-axis and q-axis currents change slowly. Fig. 8(b) is the flux vector difference between the reference and adjustable models. The flux vector difference is defined as the distance between the flux positions estimated by the voltage model and current model, which is given by:

그림입니다.
원본 그림의 이름: CLP00001b0c003a.bmp
원본 그림의 크기: 가로 1045pixel, 세로 106pixel     (25)


Fig. 8. Well-tuned case for all of the stator parameters (600 rpm, 50 % load). (a) 그림입니다.
원본 그림의 이름: CLP00001b0c0033.bmp
원본 그림의 크기: 가로 63pixel, 세로 77pixel, 그림입니다.
원본 그림의 이름: CLP00001b0c0034.bmp
원본 그림의 크기: 가로 79pixel, 세로 80pixel, and currents. (b) 그림입니다.
원본 그림의 이름: CLP00001b0c0035.bmp
원본 그림의 크기: 가로 101pixel, 세로 64pixel.

그림입니다.
원본 그림의 이름: CLP00001b0c0031.bmp
원본 그림의 크기: 가로 1390pixel, 세로 917pixel

(a)

 

그림입니다.
원본 그림의 이름: CLP00001b0c0032.bmp
원본 그림의 크기: 가로 1778pixel, 세로 793pixel

(b)


The flux vector difference converges to zero when the parameters are well tuned. This implies that the adjustable model follows the reference model both on the d-axis and the q-axis.

Fig. 9(a) shows results of the detuned case of Rs. Rs is set to 1.15 times the actual value. In this situation, the relative error of Lm is -6.7 % and that of Rr is -0.7 % in the steady state. Finally, the detuned case of σLs  is shown in Fig. 9(b). It is usually harder to tune σLs than Rs. Thus, σLs is set to -25 %, to have a higher error than the Rs case. In this case, the relative errors of estimation are 0.6 % and 2.5 % for Lm and Rr, respectively. In the detuned cases of the stator parameters, the errors of the estimated Lm and Rr at various speeds and loads are enumerated in Table II and Table III. In terms of their speed and load, their errors are different. However, these errors are not as obvious as these are low. For the detuned case of σLs, equal errors appear as different speeds. In conclusion, the numerical results in Table II and Table III demonstrate the good robustness of the proposed algorithm for detuning the stator parameters.


Fig. 9. Detuned case of the stator parameters. (a) 그림입니다.
원본 그림의 이름: CLP00001b0c0038.bmp
원본 그림의 크기: 가로 317pixel, 세로 79pixel case (300 rpm, 80 % load). (b) 그림입니다.
원본 그림의 이름: CLP00001b0c0039.bmp
원본 그림의 크기: 가로 378pixel, 세로 80pixel case (1200 rpm, 30 % load).

그림입니다.
원본 그림의 이름: CLP00001b0c0036.bmp
원본 그림의 크기: 가로 1403pixel, 세로 920pixel

(a)

 

그림입니다.
원본 그림의 이름: CLP00001b0c0037.bmp
원본 그림의 크기: 가로 1398pixel, 세로 915pixel

(b)


TABLE II ERRORS IN THE ESTIMATED VALUES OF 그림입니다.
원본 그림의 이름: CLP0000242c0002.bmp
원본 그림의 크기: 가로 60pixel, 세로 56pixel AND 그림입니다.
원본 그림의 이름: CLP0000242c0003.bmp
원본 그림의 크기: 가로 68pixel, 세로 58pixel IN 그림입니다.
원본 그림의 이름: CLP00001b0c0046.bmp
원본 그림의 크기: 가로 315pixel, 세로 76pixel IN VARIOUS SITUATIONS (UNIT: %)

Speeds (rpm)

30 % load

50 % load

80 % load

그림입니다.
원본 그림의 이름: CLP0000242c0002.bmp
원본 그림의 크기: 가로 60pixel, 세로 56pixel

그림입니다.
원본 그림의 이름: CLP0000242c0003.bmp
원본 그림의 크기: 가로 68pixel, 세로 58pixel

그림입니다.
원본 그림의 이름: CLP0000242c0002.bmp
원본 그림의 크기: 가로 60pixel, 세로 56pixel

그림입니다.
원본 그림의 이름: CLP0000242c0003.bmp
원본 그림의 크기: 가로 68pixel, 세로 58pixel

그림입니다.
원본 그림의 이름: CLP0000242c0002.bmp
원본 그림의 크기: 가로 60pixel, 세로 56pixel

그림입니다.
원본 그림의 이름: CLP0000242c0003.bmp
원본 그림의 크기: 가로 68pixel, 세로 58pixel

300

5.0

-2.7

1.7

-4.4

-0.7

-6.7

600

2.4

-1.5

0.8

-2.4

-0.5

-3.7

1200

1.2

-0.8

0.4

-1.3

-0.3

-2.0


TABLE III ERRORS IN THE ESTIMATED VALUES OF 그림입니다.
원본 그림의 이름: CLP0000242c0002.bmp
원본 그림의 크기: 가로 60pixel, 세로 56pixel AND 그림입니다.
원본 그림의 이름: CLP0000242c0003.bmp
원본 그림의 크기: 가로 68pixel, 세로 58pixel IN 그림입니다.
원본 그림의 이름: CLP00001b0c0047.bmp
원본 그림의 크기: 가로 377pixel, 세로 78pixel IN VARIOUS SITUATIONS (UNIT: %)

Speeds (rpm)

30 % load

50 % load

80 % load

그림입니다.
원본 그림의 이름: CLP0000242c0002.bmp
원본 그림의 크기: 가로 60pixel, 세로 56pixel

그림입니다.
원본 그림의 이름: CLP0000242c0003.bmp
원본 그림의 크기: 가로 68pixel, 세로 58pixel

그림입니다.
원본 그림의 이름: CLP0000242c0002.bmp
원본 그림의 크기: 가로 60pixel, 세로 56pixel

그림입니다.
원본 그림의 이름: CLP0000242c0003.bmp
원본 그림의 크기: 가로 68pixel, 세로 58pixel

그림입니다.
원본 그림의 이름: CLP0000242c0002.bmp
원본 그림의 크기: 가로 60pixel, 세로 56pixel

그림입니다.
원본 그림의 이름: CLP0000242c0003.bmp
원본 그림의 크기: 가로 68pixel, 세로 58pixel

300

2.5

0.6

2.6

-0.6

2.8

-3.2

600

1200



Ⅶ. EXPERIMENTAL RESULTS

The information on the induction motor used in the experiment has been specified in Table I. The utilized DSP in the experiment is a TMS320F28335. The switching frequency of the PWM inverter is 5 kHz, and the sampling frequency of the line currents is 10 kHz. All of the settings for the control period are same as those for the aforementioned simulations. The utilized load motor is a 2.26 kW synchronous motor and it is set to run in the torque control mode to generate a specific torque. Before the drive implementation, in order to obtain accurate output voltages, the methods for dead time compensation in [23] were applied, where the switch compensation time is calibrated considering the dead time, the voltage drops of the power devices, etc. Rs was identified before the speed drive. Rs was 1.67 Ω at the beginning of the experiment. The initial settings of Lm and Rr are 202.5 mH and 0.35 Ω which have errors of approximate 50 %.

In the well-tuned case of Rs and σLs, results of parameter identification are shown in Fig. 10(a). The two parameters are estimated correctly. These results correspond to Fig. 8(a) from the simulation results. The currents change a bit with the identification. The voltages and speed are plotted in Fig. 10(b). The speed remains constantly at 600 rpm during both the steady-state and transient-state under the identification algorithm. The voltages decrease as the parameters are being tuned. Fig. 10(c) shows the flux vector difference between the reference and the adjustable models. 그림입니다.
원본 그림의 이름: CLP00001b0c004d.bmp
원본 그림의 크기: 가로 109pixel, 세로 66pixel converges to zero asymptotically as Lm and Rr are being tuned.


Fig. 10. Well-tuned case for all of the stator parameters (600 rpm, 50 % load). (a) 그림입니다.
원본 그림의 이름: CLP00001b0c003e.bmp
원본 그림의 크기: 가로 66pixel, 세로 79pixel, 그림입니다.
원본 그림의 이름: CLP00001b0c003f.bmp
원본 그림의 크기: 가로 75pixel, 세로 75pixel and currents. (b) Voltages and speed. (c) 그림입니다.
원본 그림의 이름: CLP00001b0c0040.bmp
원본 그림의 크기: 가로 101pixel, 세로 63pixel.

그림입니다.
원본 그림의 이름: CLP00001b0c003b.bmp
원본 그림의 크기: 가로 1385pixel, 세로 911pixel

(a)

 

그림입니다.
원본 그림의 이름: CLP00001b0c003c.bmp
원본 그림의 크기: 가로 1381pixel, 세로 918pixel

(b)

 

그림입니다.
원본 그림의 이름: CLP00001b0c003d.bmp
원본 그림의 크기: 가로 1380pixel, 세로 635pixel

(c)


Fig. 11 and Fig. 12 describe the results of the detuned case of Rs and σLs, respectively. The two detuned cases show the similar results to Fig. 10, although errors of the stator parameters exist. Since the two parameters are tuned simultaneously, the currents and voltages change. Meanwhile, the speeds stay at the command speed. Data collected from 27 situations (speeds: 300/600/1200 rpm, loads: 30/50/80 %, and the tuned-states of Rs or σLs) are shown in Fig. 13(a). This data shows consistent results under various conditions. As a result, these experimental results proves that the proposed algorithm is robust to the detuning of the stator parameters.


Fig. 11. 그림입니다.
원본 그림의 이름: CLP00001b0c0043.bmp
원본 그림의 크기: 가로 313pixel, 세로 72pixel case (300 rpm, 80 % load). (a) 그림입니다.
원본 그림의 이름: CLP00001b0c0044.bmp
원본 그림의 크기: 가로 66pixel, 세로 74pixel, 그림입니다.
원본 그림의 이름: CLP00001b0c0045.bmp
원본 그림의 크기: 가로 76pixel, 세로 76pixel and currents. (b) Voltages and speed.

그림입니다.
원본 그림의 이름: CLP00001b0c0041.bmp
원본 그림의 크기: 가로 1315pixel, 세로 867pixel

(a)

 

그림입니다.
원본 그림의 이름: CLP00001b0c0042.bmp
원본 그림의 크기: 가로 1315pixel, 세로 863pixel

(b)


Fig. 12. 그림입니다.
원본 그림의 이름: CLP00001b0c004a.bmp
원본 그림의 크기: 가로 382pixel, 세로 76pixel case (1200 rpm, 30 % load). (a) 그림입니다.
원본 그림의 이름: CLP00001b0c004b.bmp
원본 그림의 크기: 가로 67pixel, 세로 78pixel, 그림입니다.
원본 그림의 이름: CLP00001b0c004c.bmp
원본 그림의 크기: 가로 78pixel, 세로 75pixel and currents. (b) Voltages and speed.

그림입니다.
원본 그림의 이름: CLP00001b0c0048.bmp
원본 그림의 크기: 가로 1344pixel, 세로 893pixel

(a)

 

그림입니다.
원본 그림의 이름: CLP00001b0c0049.bmp
원본 그림의 크기: 가로 1343pixel, 세로 860pixel

(b)


Fig. 13. Sets of estimated parameters under various conditions. (a) Beginning of experimental activity. (b) Two hours after experimental activity.

그림입니다.
원본 그림의 이름: CLP00001b0c004e.bmp
원본 그림의 크기: 가로 1377pixel, 세로 877pixel

(a)

 

그림입니다.
원본 그림의 이름: CLP00001b0c004f.bmp
원본 그림의 크기: 가로 1364pixel, 세로 795pixel

(b)


Two hours after the speed drive operation, Rs was identified again. Rs increased a bit than the beginning, which was 1.78 Ω. In this situation, experiments in the selected condition were carried out repeatedly. Fig. 13(b) shows these results. The results for Rr have 6~7 % higher values than those for the previous experiments. Lm has values that are equal to those of the previous experiment.

In practical speed drive systems, coupled loads fluctuate for any reason. Thus, it is necessary to prove that the proposed method operates properly under such conditions. It is expected that the proposed algorithm can operate smoothly because it does not have any assumptions for the steady state, as can be seen in chapter V. In other words, proper operation of the algorithm can be guaranteed in states where the speed or current is not constant. In making these situations, a fluctuating load is generated, which has a constant component 50 % plus a sinusoidal component 20 % for the rated torque of the induction motor. Here, the frequency of the sinusoidal component is synchronized with the speed. Fig. 14 shows the experimental results at 300, 600 and 1200 rpm. The Q-axis current fluctuates to compensate the fluctuating load. The speed also fluctuates since load torque is not perfectly compensated. Despite this condition, the parameters are well identified without fluctuations. Thus, the proposed identification method may be used in practical speed drive systems with fluctuating loads.


Fig. 14. In cases of speed operation with fluctuating load. (a) At 300 rpm. (b) At 600 rpm. (c) At 1200 rpm.

그림입니다.
원본 그림의 이름: CLP00001b0c0050.bmp
원본 그림의 크기: 가로 1362pixel, 세로 833pixel

(a)

 

그림입니다.
원본 그림의 이름: CLP00001b0c0051.bmp
원본 그림의 크기: 가로 1357pixel, 세로 834pixel

(b)

 

그림입니다.
원본 그림의 이름: CLP00001b0c0052.bmp
원본 그림의 크기: 가로 1363pixel, 세로 836pixel

(c)



Ⅷ. CONCLUSIONS

In this paper, a novel method for the identification of Lm and Rr has been proposed. In order to realize the proposed method, the MRAC principle is applied, and RLS estimation is utilized for estimating these parameters. Because a closed- loop Gopinath style flux observer is employed as the reference model in the adaptation mechanism, the proposed algorithm can operate at a sufficiently high frequency to get the complete reference model. However, the proposed algorithm has the following advantages:

- Simple structure and well-known algorithms: flux estimations, filtering calculations, and RLS estimation.

- The two parameters can be identified simultaneously.

- Independence of the two parameters: the detuning of Lm does not affect the estimation accuracy of Rr and vice versa.

- It guarantees effective operation in both the steady- state and the transient-state.

The proposed method was verified by simulation and experimental results where the parameters were estimated accurately through a speed control test with a load. Additionally, in detuned cases of the stator parameters and cases with a fluctuating load, experimental results demonstrated the excellence of the proposed algorithm in diverse situations.



ACKNOWLEDGMENT

This work was supported by the Korea Institute of Energy Technology Evaluation and Planning (KETEP) and the Ministry of Trade, Industry & Energy (MOTIE) of the Republic of Korea (No. 20174030201490).



REFERENCES

[1] S. Wade, M. W. Dunnigan, and B. W. Williams, “A new method of rotor resistance estimation for vector-controlled induction machines,” IEEE Trans. Ind. Electron., Vol. 44, No. 2, pp. 247-257, Apr. 1997.

[2] Ba-Razzouk, A. Chérity, and P. Sicard, “Implementation of a DSP based real-time estimator of induction motors rotor time constant,” IEEE Trans. Pow. Electron., Vol. 17, No. 4, pp. 534-542, Jul. 2002.

[3] S. -H. Jeon, K. -K. Oh, and J. -Y. Choi, “Flux observer with online tuning of stator and rotor resistances for induction motors,” IEEE Trans. Ind. Elec., Vol. 49, No. 3, pp. 653-664, Jun. 2002.

[4] F. R. Salmasi, T. A. Najafabadi, and P. J. Maralani, “An adaptive flux observer with online estimation of DC-link voltage and rotor resistance for VSI-based induction motors,” IEEE Trans. Pow. Electron., Vol. 25, No. 5, pp. 1310-1319, May 2010.

[5] L. Zhao, J. Huang, J. Chen, and M. Ye, “A parallel speed and rotor time constant identification scheme for indirect field oriented induction motor drives,” IEEE Trans. Power Electron., Vol. 31, No. 9, pp. 6494-6503, Sep. 2016.

[6] X. Yu, M. W. Dunnigan, and B. W. Williams, “A novel rotor resistance identification method for an indirect rotor flux-orientated controlled induction machine system,” IEEE Trans. Power Electron., Vol. 17, No. 3, pp. 353-364, May 2002.

[7] S. Maiti, C. Chakraborty, Y. Hori, and M. C. Ta, “Model reference adaptive controller-based rotor resistance and speed estimation techniques for vector controlled induction motor drive utilizing reactive power,” IEEE Trans. Ind. Electron., Vol. 55, No. 2, pp. 594-601, Feb. 2008.

[8] K. Wang, B. Chen, G. Shen, W. Yao, K. Lee, and Z. Lu, “Online updating of rotor time constant based on combined voltage and current mode flux observer for speed- sensorless AC drives,” IEEE Trans. Ind. Elec., Vol. 61, No. 9, pp. 4583-4593, Sep. 2014.

[9] D. P. Marćetić and S. N. Vukosavić, “Speed-sensorless AC drives with the rotor time constant parameter update,” IEEE Trans. Ind. Electron., Vol. 54, No. 5, pp. 2618-2625, Oct. 2007.

[10] A. Yoo, “Robust on-line rotor time constant estimation for induction machines,” J. Power Electron., Vol. 14, No. 5, pp. 1000-1007, 2014.

[11] D. M. Reed, H. F. Hofmann, and J. Sun, “Offline identification of induction machine parameters with core loss estimation using the stator current locus,” IEEE Trans. Energy Convers., Vol. 31, No. 4, Dec. 2016.

[12] S.-K. Sul, “A novel technique of rotor resistance estimation considering variation of mutual inductance,” IEEE Trans. Ind. Appl., Vol. 25, No. 4, pp. 578-587, Jul./Aug. 1989.

[13] S.-H. Lee, A. Yoo, H.-J. Lee, Y.-D. Yoon, and B.-M. Han, “Identification of induction motor parameters at standstill based on integral calculation,” IEEE Tran. Ind. App., Vol. 53, No. 3, May/Jun. 2017.

[14] C. Attaianese, A. Damiano, G. Gtto, I. Marongiu, and A. Perfetto, “Induction motor drive parameters identification,” IEEE Trans. Power Electron., Vol. 13, No. 6, pp. 1112-1122, Nov. 1998.

[15] D. Telford, M. W. Dunnigan, and B. W. Williams, “Online identification of induction machine electrical parameters for vector control loop tuning,” IEEE Trans. Ind. Electron., Vol. 50, No. 2, pp. 253-261, Apr. 2003.

[16] Y. He, Y. Wang, Y. Feng, and Z. Wang, “Parameter identification of an induction machine at standstill using the vector constructing method,” IEEE Trans. Power Electron., Vol. 27, No. 2, pp. 905-915, Feb. 2012.

[17] B. Proca and A. Keyhani, “Sliding-mode flux observer with online rotor parameter estimation for induction motors,” IEEE Trans. Ind. Electron., Vol. 54, No. 2, pp. 716-723, Apr. 2007.

[18] M. Wlas, Z. Krzemiński, and H. A. Toliyat, “Neural- network-based parameter estimations of induction motors,” IEEE Trans. Ind. Electron., Vol. 55, No. 4, pp. 1783-1794, Apr. 2008.

[19] X. Xu and D. W. Novotny, “Implementation of direct stator flux orientation control on a versatile DSP based system,” IEEE Trans. Ind. Appl., Vol. 27, No. 4, pp. 694-700, Jul./Aug. 1991.

[20] P. L. Jansen and R. D. Lorenz, “A physically insightful approach to the design and accuracy assessment of flux observer for field oriented induction machine drives,” IEEE Trans. Ind. Appl., Vol. 30, No. 1, pp. 101-110, Jan./Feb. 1994.

[21] I. Takahasi and T. Noguchi, “A new quick-response and high-efficiency control strategy of an induction motor,” IEEE Trans. Ind. Appl., Vol. IA-22, No. 5, pp. 820-827, Sep./Oct. 1986.

[22] M. H. Hayes, “Recursive Least Squares,” in Statistical Digital Signal Processing and Modeling, John Wiley & Sons, Inc., Chap. 9.4, pp. 541-553, 1996.

[23] J.-W. Choi and S.-K. Sul, “Inverter output voltage synthesis using novel dead time compensation,” IEEE Trans. Power Electron., Vol. 11, No. 2, pp. 221-227, Mar. 1996.



그림입니다.
원본 그림의 이름: image62.jpeg
원본 그림의 크기: 가로 190pixel, 세로 225pixel

Gwon-Jae Jo was born in Daegu, Korea. He received his B.S. and M.S. degrees in Electrical Engineering from Kyungpook National University, Daegu, Korea, in 2013 and 2016, respectively, where he is presently working towards his Ph.D. degree in Electrical Engineering. His current research interests include AC motor drives and custom power devices.


그림입니다.
원본 그림의 이름: image61.jpeg
원본 그림의 크기: 가로 169pixel, 세로 225pixel

Jong-Woo Choi was born in Daegu, Korea. He received his B.S., M.S. and Ph.D. degrees in Electrical Engineering from Seoul National University, Seoul, Korea, in 1991, 1993 and 1996, respectively. From 1996 to 2000, he worked as a Research Engineer at LG Industrial Systems Co., Korea. Since 2001, he has been a Professor in the Department of Electrical Engineering at Kyungpook National University, Daegu, Korea. His current research interests include static power conversion and electric machine drives.