Numerical Estimation and Experimental Verification of Optimal Parameter Identification Based on Modern Optimization of a Three Phase Induction Motor

: The parameters of electric machines play a substantial role in the control system which, in turn, has a great impact on machine performance. In this paper, a proposed optimal estimation method for the electrical parameters of induction motors is presented. The proposed method uses the particle swarm optimization (PSO) technique. Further, it also considers the influence of temperature on the stator resistance. A complete experimental setup was constructed to validate the proposed method. The estimated electrical parameters of a 3.8-hp induction motor are compared with the measured values. A heat run test was performed to compare the effect of temperature on the stator resistance based on the proposed estimation method and the experimental measurements at the same conditions. It is shown that acceptable accuracy between the simulated results and the experimental measurements has been achieved.


Introduction
Induction motors (IMs) are employed extensively in several industrial applications because of their simple and rugged structure, easy operation and low cost maintenance. The induction motors can work direct online or via a controlled system. The main aim of using a control system for induction motors is to improve their dynamic performance. Several control techniques exist and are employed on induction motors, for example, vector control, direct torque control and sensorless control [1]. The electrical parameters of the motor are very essential to design, performance assessment, and feasibility of the control techniques. Therefore, it is very important to estimate the parameters of the electric motor accurately. Any difference between the actual motor parameters and those implemented in the drive system will significantly affect the system performance as well as the operation characteristics. Thereby, it is necessary to monitor the variation of the parameters during motor operation [2][3][4].
Most conventional techniques that are used to estimate the induction motor parameters require locked-rotor and no-load tests. However, such techniques cannot easily be applied due to mechanical rotor locking. In addition, in the locked rotor test, the rotor frequency is the same as the supply frequency; however, during operating conditions, the value of rotor frequency is lower than the supply. Therefore, such improper rotor frequency has an effect on locked rotor test results. In addition, in motors with high power, this test is impractical. These problems have encouraged researchers to investigate different procedures for motor parameter estimation. An analytical noniterative method is presented in Reference [2] to estimate the parameters of the squirrel cage induction motor representing core losses from catalog or nameplate data. In Reference [3], the electrical parameters of squirrel cage and double-cage induction motors, considering and neglecting iron losses for both models by measuring current, voltage, rotor speed and electromagnetic torque, are identified. A review of parameter values estimation of three phase IMs based on information provided by manufacturers catalogs, in particular by those dedicated to efficiency estimation, has been presented in Reference [4].
Modern optimization has been used in different engineering applications like in References [5][6][7][8]. The idea of parameters estimation of electric machines is examined by several techniques like Particle Swarm Optimization (PSO), Extended Kaman filter (EKF) and Genetic algorithm (GA). For example, Reference [9] presents the application of PSO for estimating and tracking the changes of IM parameters using measured three phase currents, voltages and rotor speed. But this method applied 80 particle numbers and takes about 40 iterations to converge the actual values to the estimated ones. In Reference [10] an online estimation of IM parameters with Kalman filter is suggested, but the modelling errors cause state noise in the system. Rotor resistance and the time constant of IM are also identified using a fourth-order sliding-mode flux observer as presented in Reference [11] but high errors are obtained in their values when operating the motor at low torques. Moreover, rotor and stator resistances online estimation, using adaptive observer while motor speed and one phase current are available and high observer gains, are applied to achieve observer stability [12]. A GA is employed to estimate the values of the parameters of squirrel cage IM through some experimental data sets in Reference [13]. GA is also used in References [14][15][16][17] for determining the parameters of IM based on the manufacturer's published data. The Particle Swarm Optimization (PSO) algorithm is employed to determine the parameters of inductions motors in several publications, for example, References [18][19][20][21].
In this paper, a PSO algorithm for estimating the electrical parameters of IM is presented. In addition, a comparison between the performance of the IM based on the measured and estimated parameters using the PSO algorithm is provided. Moreover, a heat run test is performed to study the effect of temperature on the stator circuit and compare it with that obtained using the proposed PSO at the same conditions. The manuscript has been organized in two main sections; Section 2 presents induction motor parameters estimation based on PSO algorithm; Section 3 shows the results and discussion using both simulations and experimental systems.

IM Parameter Determination Based on PSO Algorithm
The electrical equivalent circuit of the IM at steady-state is shown in Figure 1. In this circuit, five electrical parameters are required to simulate the performance of the machine. The five parameters are the stator and rotor resistances ( and ), stator and rotor leakage inductances ( , ) and magnetizing inductance ( ) [17][18][19][20][21][22][23].
In order to estimate these five parameters, a modified PSO technique is proposed hereafter. The PSO algorithm depends on the behavior of a bird's motion searching for food to trace its evolution. It has a figure of particles to generate a swarm. For a PSO system, the particles try to adjust their position depending on its own knowledge and the knowledge of neighboring particles and use the best position obtained from themselves and their neighbors [24].
Assume X and V are the particle's position and its consistent velocity in the search space respectively. At iteration k, the particle's position and velocity could be estimated from the following equations: where k is the iteration step, i is the particle, W denotes the weight factor, C1 and C2 are constants used to scale the contribution of the cognitive and social components, respectively, β1 and β2 denote random values [0-1] by which stochastic elements are introduced in the algorithm and these parameters do not depend on time; Pi and Pg are the best position of the individual's and of the group's respectively.

Parametric Characterization
Considering a system has an input and an output . We can use the following state space form to describe the induction motor model.
where: x is states vector, is input vector, and is output vector as follows: where and are the stator and rotor resistances, and are stator and rotor leakage inductances, is the magnetizing inductance, , , are the q-axis and d-axis stator and rotor voltages respectively, , , are the q-axis and d-axis stator and rotor currents respectively, , , , are the q-axis and d-axis stator and rotor flux linkages and is the rotor speed.
The parameters of the actual induction motor are represented in a vector called "θ" which is as follows: Since is the unknown parameters vector and its value must be identified exactly, a modeled system similar to the actual one is constructed with the same initial condition From this description, the input is provided for the actual and the modeled systems as Equations (3) and (4). The states and outputs of the actual system are ascertained based upon the information of , and the estimated parameters. x ∧ and y ∧ denote the states and the output of the model. To assess the parameters to be distinguished, the fitness function for this identification problem will be contrasted.
To serve this purpose, the objective function can be formulated by way of [25]: where id k and iq k are suitable gains and T is a suitable time.
The two values of ids and iqs are used as inputs to the optimizer. The cost function represented by the square of the difference values between the real and estimated output. This cost function is minimized during the optimization process in order to find out the optimal parameters of IM. The cost function will be zero or minimized only if the estimated IM parameters nearly equal their actual ones, that is, θ θ ∧ = , to simplify the identification process, the IM system model in Equation (3) must be discretized.
The configuration of the parameter identification methodology in this study is demonstrated in Figure 2. The steps of the optimization process are explained in Figure 3. First, the system input u is known to the actual and modeled systems. After that, the outputs of them are compared and fed to a fitness function. Finally, the calculated value of this function is used as an input to the PSO identifier for the identification process.
In order to apply the proposed optimization method to the IM model, a proper cost function is essential to be well-defined. From Equation (9) The motor parameters recieves the current particle position The currents i ds and i qs are calculated by (4) The current error is calculated by (9), (10) and (11) The local best parameters update criterion are applied The particles positions and velocities are updated by (1) and (2) yes No iteration index = iteration index+ 1

Convergence criterion = True? No
The global best parameters update criterion are applied End Particle index = maximum ?
yes Figure 3. Flowchart of PSO algorithm.

Non-Parametric Characterization
For achieving steady state performance, studying these characteristics will be very essential and they are the stator current ( ), power factor ( ), efficiency ( ), input power ( ) and output power ( ). During this test, the IM is mechanically loaded with a DC generator of almost equal rating load and the stator current, rotor speed, the input and the output powers are measured. The overview of the virtual operation of the induction motor is illustrated in Figure 4.

Characteristic of Induction Motor
To display the performance of the proposed PSO algorithm, the complete model presented before is simulated using the Matlab environment. The simulation is performed at low speed of 100 rpm and loading of 10 N.m. The estimated motor parameters are plotted versus iterations number as shown in Figures 5-8. The PSO algorithm runs with five particles in the swarm and 200 iterations. The stator resistance versus iterations number is plotted in Figure 5 and shows that the parameter tries to be identified at iteration 38. Further, when the stator resistance is increased by 50% from the nominal value (due to increasing the temperature of stator windings), the estimated value corresponds with the actual value only after 10 iterations. Similarly, the rotor resistance is reported in Figure 6, which shows that the algorithm has a reduced sensitivity and begins to precisely track the actual value within 20 iterations. The prediction of magnetizing and stator inductances ( and = ) are confirmed in Figures 7 and 8. It is clear that the PSO system converges faster with less fluctuation and greatly outperformed the others.
In this part, a comparison between three methods including performance depending on calculation based on standard data specification of the motor (see Appendix A), on data obtained from PSO system and finally on a direct loading test of the induction motor. To confirm the effectiveness of the proposed method in parameter identification, computer simulations which have been carried out using MATLAB software results will be compared with others from experimental tests. The motor data specifications are identical to those used in the simulation. Figures 9-12 show the results obtained from the parameters and data specifications of the manufacturer, PSO algorithm and experimental measurements of the induction motor. Figures 9 and  10 show standard, estimated and experimental value of stator current and input power factor. These figures indicate that the estimated values are well matched and agree with their measured ones. On the other hand, the large deviation between the standard values of stator current and input power factor and their experimental ones is owed to the consideration that the parameters of the machine are assumed to be constant at all operating conditions.      Figure 12 shows the standard, estimated and experimental values of the motor input and output power versus motor speed, respectively. From this figure, it is clear that the estimated values of both motor input and output power and measured values have small deviation. The discrepancy between estimated and experimental values, especially at high values of motor speed, is attributed to the mechanical loss and magnetization saturation, which is disregarded from the equivalent circuit.

Temperature Influence on Stator Resistance
To investigate the temperature effect on induction motor stator resistance at different loading conditions using the proposed estimation method, a heat run test was performed. The induction motor was run on 3 phase supply and loaded by a DC generator coupled to it as shown in Figure 4. Different loading conditions were performed by adjusting load resistance across the DC generator. The motor was run for four hours at each load and after that time, and at motor standstill, the stator windings resistance were measured using the DC test but this measured value was multiplied by a factor to convert it from its DC value to its AC value. This is done to account for the skin effect.
The stator circuit resistance was estimated via PSO at different loading conditions and then compared with the measured values acquired from the heat run test. Figure 13 demonstrates that the values of stator circuit resistance, which have been obtained experimentally from the direct loading test, are in agreement with the estimated ones. The discrepancy between estimated and experimental values is attributed to the mechanical loss and magnetization saturation which have been neglected.

Conclusions
An estimation method of induction motor electrical parameters based on the Particle Swarm Optimization (PSO) algorithm has been presented in this paper. Besides, comparison between the motor performance indicators, that is, power factor, efficiency and output power, obtained from experimental measurements and simulated based on the proposed parameters estimation method is provided. Moreover, the effect of temperature on the stator resistance is also shown. It is proved that the proposed estimation method of the motor parameters is capable of predicting an accurate motor performance, with an acceptable error. The main reason of the error is due to neglecting some effects on the implemented mathematical model of the induction motor such as mechanical losses and saturation of the magnetization inductance.
Finally, experimental measurements of the steady state performance and heat run test of the motor validate the theoretical results obtained using the proposed method.