Hybrid Swarm Algorithms for Parameter Identification of an ActuatorModel in an Electrical Machine

Efficient identification and control algorithms are needed, when active vibration suppression techniques are developed for industrial machines. In the paper a new actuator for reducing rotor vibrations in electrical machines is investigated. Model-based control is needed in designing the algorithm for voltage input, and therefore proper models for the actuator must be available. In addition to the traditional prediction error method a new knowledge-based Artificial Fish-Swarm optimization algorithm (AFA) with crossover, CAFAC, is proposed to identify the parameters in the new model. Then, in order to obtain a fast convergence of the algorithm in the case of a 30 kW two-pole squirrel cage induction motor, we combine the CAFAC and Particle Swarm Optimization (PSO) to identify parameters of the machine to construct a linear time-invariant(LTI) state-space model. Besides that, the prediction error method (PEM) is also employed to identify the induction motor to produce a black box model with correspondence to input-output measurements.


Introduction
One of the fundamental processes of control is the parameter identification.In control engineering a lot of effort has been done to develop methods to identify the system model and its parameters.A wide range of techniques such as the least squares method, the maximum likelihood method, and the cross correlation method, exist for system identification [1].However, a drawback of traditional identification methods is their dependence on unrealistic assumptions such as unimodal performance landscapes and differentiability of the cost function.Consequently, some problems are oversimplified to fulfil such assumptions.Stochastic search algorithms such as evolutionary-based algorithms seem to be a promising alternative to these conventional methods.Swarm Intelligence is one of the most important branches for Evolutionary Algorithms, which has been widely used in identification problems [2][3][4].
The Artificial Fish-Swarm Algorithm (AFA) and the Particle Swarm Optimization (PSO) are two kinds of typical Swarm Intelligence methods [5].The AFA is an animal behaviour-based optimization method and, like the PSO, which was firstly developed by Li and Eberhart and Kennedy, respectively [6][7][8][9].PSO has been widely used in identification in many research fields [10][11][12][13][14][15].However, there are few researches to apply the AFA in such kind of problems.The AFA can search for the global optimum effectively and has a certain adaptive ability for searching space.But the AFA individual behaviour is to hunt for local optimum.Therefore, avoiding individual premature becomes difficult.In this case, artificial fish will be stuck into local optima when dealing with multimodal optimization problems.To improve the performance of AFA is an im-portant challenge for applying the AFA in real problems [16].
In order to guide the evolutionary-based stochastic algorithms, a novel optimization method, the Cultural Algorithm (CA) proposed by Reynolds in 1994 was developed to demanding problems, due to its flexibility and efficiency [17].The CA is a class of computational models derived from the principles of the culture evolution in nature, composed by the population space, the belief space and a communication protocol.In the CA, the evolution takes place in the population space.The knowledge generated from the experiences of the individuals is stored in the belief space.The communication protocol can be depicted as an acceptance function and an influence function, which is utilized to make interaction between the other two components [18,19].Various evolutionary algorithms have been utilized in the population space of the CA leading to good results [20][21][22][23][24][25][26][27].The framework of the culture algorithm can be depicted in Figure 1.
The motivation behind this work is twofold.Firstly, in the literature there seems to be no results on hybrid AFA and CA.The knowledge stored in the belief space do act as a significant role in the process of evolution, and we try to find a proper framework for updating the belief space in combining the CA and AFA.Secondly, the new hybrid algorithm is applied to identify the parameters of a new kind of actuator, which is used to suppress rotor vibrations of an electronic machine.
As part of our investigation, two typical kinds of knowledge in CA, the situation knowledge and the normative knowledge, are stored in the belief space to update the population space and to establish the relationship between the two spaces in the CA.The performance of the CAFAC is explored using offline parameter identification of the actuator-rotor system in an electrical machine.The identification is performed based on a lower-order physical linear time invariant (LTI) parametric state-space model of the actuator-rotor system.The resulting model can be utilized to design model-based active control algorithms for vibrations reduction.
The rest of this paper is organized as follows.Section 2 briefly introduces the background knowledge of the system under control and the AFA.Section 3 proposes and discusses the underlying principle of the CAFAC.In Section 4 the identification of an actuator-rotor system is considered using CAFAC and PEM and simulation results showing the control performance is established.A few conclusions and remarks are given in Section 5.

Problem Statement and Preliminaries
2.1.The Parametric Model of the Actuator-Rotor System.We examine a two-pole cage induction motor equipped with a built-in force actuator, which actively generates force on the rotor (Figure 2(a)).The actuator is a four-pole extra winding inside the stator slots (Figure 2(b)).This design allows producing additional forces without harming the normal operation of the motor.The actuator is controlled by two voltage signals that by design induce forces on the rotor in horizontal and vertical directions (x and y).
The most important measurements for the identification are obtained using eddy current sensor.The sensors, conforming to the control signals, measure the rotor position also in horizontal and vertical directions.With them it is possible to record the rotor movement in any conditions accurately.Another set of sensors has been added on the right side of the motor in order to improve accuracy.Another important measurement devise is the encoder that provides the rotational angle and frequency of the rotor.
In the identification measurement the motor was operating at 32.085 Hz and the excitation input (control voltages in horizontal and vertical directions) was a uniform random number signal with frequency content up to 500 Hz.The output data was then processed so that the effects of vibrations were removed and only the response to the excitation signals remained [28].The processed measurement is displayed in Figure 3.
The data from this setup is used to obtain a mathematical actuator model that can be used for control design purposes, when the objective is compensating rotor vibrations.The motivation of obtaining a new parametric physical model is increased understanding of the model compared to a black box model, which has been used before successfully [29].Parametric models have been identified before with data from a finite element model, but they do not describe the test motor and some of the identified parameters had complex values, although for a realistic machine the parameters such as resistances must be realvalued.
The linear time-invariant parametric (LTI) model of the system will be discussed.The model for the induction motor according to [30] is based on the mechanical model in modal coordinates [31] where η is the modal coordinate vector, u rc is the displacement in x-direction and y-direction, Φ rc is the modal matrix, Ω is a diagonal matrix, and Ξ is the modal damping matrix, f ex denotes the sinusoidal disturbance causing rotor vibration disturbance, and f c the electro mechanical force  acting on the rotor.The model structure is the same as in [32] d dt ( The input v = (Re( u c,0 ), Im( u c,0 )), is the control voltage, the first block row is the mechanical model, and i = (Re( i c,0 ), Re(M r,c,2 i r,2,0 ) T , Im( i c,0 ), Im(M r,c,2 i r,2,0 )) T , where i c,0 and i r,2,0 are current space vectors and M r,c,2 included in i is mutual inductance of control winding and the rotor cage.In variables the subscript c denotes a relation to the control winding, the subscript r denotes a relation to the rotor cage inside the motor, and the number in subscript denotes harmonic component [33].All the parameters for the model are based on measured data and their explanations are listed in Table 1 [30].
Following [30,32] a combined electromechanical model can be derived starting from the voltage-flux equation for the control winding and the equation for rotor cage four-pole harmonic in stator coordinates The four-pole harmonic flux leakages are ( By substitution The eccentricity coupling factor of the control winding multiplied with the complex conjugate of the first air-gap field harmonic [34].
M r,c,2 M c,r,2 mH Mutual inductance of rotor cage and control winding multiplied by itself and divided by rotor cage inductance of four-pole harmonic (M c,r,2 , where R r,2 is the rotor cage resistant of 4-pole harmonic.
The eccentricity coupling factor of control winding multiplied with first air-gap field harmonic.
A coefficient related to the unbalanced magnetic pull towards the shortest air-gap (gap between stator and rotor) [33].
The coupling factor of rotor cage multiplied with the complex conjugate of the first air-gap field harmonic and divided by rotor cage inductance of four-pole harmonic.
The eccentricity coupling factor of rotor cage four-pole harmonic multiplied with the first air-gap field harmonic and divided by rotor cage inductance of four-pole harmonic.
where z is a complex coordinate representation of the rotor eccentricity.In addition, the control force f c in terms of the parameters in Table 1 is given by These equations include all the parameters to identify.Since the complex coordinate system of z is decomposed in real and imaginary parts in the final model, matrices P em (t) and C em , are given by It should be noted that in the time-invariant case P em (t) is constant P em .Equation ( 6) can be written in matrix form Rearranging this will give the rest of the matrices, which are needed for the combined model Matrices A em , B em , S em , and Q em are finally obtained like C em in ( 10) These matrices and the mechanical model ( 1) can be combined in one state-space model (3).
The parameters of the electromechanical model are listed in Table 1.Certain practical constraints were set for the parameters for fulfilling physical quantities.Here the resistance is constant and known (14.5), so it does not need to be optimized.When the effect of the unbalanced magnetic pull is ignored we can get the LTI model.The input u(t) consists of the voltage input v to the control winding and the disturand the outputs are the displacements u rc of the rotor center.The system matrices A, B, and C are the functions of the unknown parameters in the model which are reformed in a vector as The corresponding significances of the unknown parameters, and the detailed derivation of the model can be found in [32].

The Artificial Fish-Swarm Algorithm.
Suppose that the problem under consideration has D-dimensions.Initialize the swarm with N artificial fish.The state of one artificial fish can be formulated as X i = (x i1 , x i2 , . . ., x iD ), i = 1, . . ., N , where X i represents the target variable for the problem under consideration.y = f (X i ) stands for the food concentration of the artificial fish currently, where it is the objective function.The basic behaviours of artificial fish can be expressed as follows (d i j = X j − X i stands for the Euclidean distance between X i and X j , "vd" is visual distance of the artificial fish individual, "s" is the size of the movement of artificial fish, δ is the crowd factor of the artificial fish) [6].
(1) Preying.The current state of the artificial fish is X i , the artificial fish selects a state X j randomly within the visual distance, such as X j = X i + rand(0, 1) × vd.If f (X j ) < f (X i ) the artificial fish moves from X i towards X j , meaning that X t i → X t+1 i .The formulation can be depicted as follows: If f (X j ) > f (X i ), the artificial fish selects another state randomly again.If the artificial fish cannot meet the requirement in a given time, then it moves one step randomly as (2) Swarming.The current state of the artificial fish is X i , nf is the number of its fellows within the visual distance, equal to the number of elements in the set of nf j=1 X j /nf and then y center = f (X center ) stands for the fitness of the centre position.If nf × y center < δ × y i , meaning that this area is not crowded, then if y center < y i , the artificial fish moves one step towards the centre position: Otherwise it executes the behaviour of preying.
(3) Chasing.The current state of the artificial fish is X i , X min stands for the best artificial fish individual within X i 's visual distance.nf is the number of X min 's fellows within the visual distance.y min = f (X min ), if y min < y i and nf × y min < δ × y i , the artificial fish moves one step to X min : Otherwise it executes the behaviour of preying.

Hybrid Optimization Methods of Articificial Fish Algorithm and Culture Algorithm
Here, the fish swarm is regarded as the population space, where the domain knowledge is extracted from.Then the domain knowledge is formed and stored in belief space so as to model and impact the evolution of the population at iteration.In four versions of CAFAC, we use the situation knowledge and the normative knowledge to guide the direction and the step size of the evolution.Both of them can be depicted as follows.

Structures of Belief Space in CAFAC.
The situational knowledge provides a set of best individuals available for the interpretation of specific individual experience [19].Here, the situational exemplar set consists of only the best fish found so far, S = S t | S t = {s t 1 , s t 2 , . . ., s t n } , S t stands for the best fish in the swarm at iteration t.In other words, it can be initialized with the best fish in the initial fish swarm and updated by the following function: where X t+1 gbest denotes the best artificial fish individual in the fish swarm at generation t + 1.
The normative knowledge can give the feasible solution space of the optimization problems under consideration [19].It is a set of information for each variable and is given as where U, L and D are n-dimensional vectors, and I = {x | l ≤ x ≤ u}, n is the number of the variables, l j and u j are the lower and upper bounds for the jth variable, respectively, L j and U j are the values of the fitness function associated with the bound l j and u j .Generally, l j and u j are initialized with the lower and upper bounds of individuals.L j and U j are usually initialized with positive infinity.The formulation for normative knowledge updating can be depicted in Table 2.

Acceptance Function in CAFAC.
The acceptance function determines which individuals and their performances can have impact on the knowledge in the belief space.The number of the individuals accepted for the update of the belief space is obtained according to the following function [1]: where N is the size of the swarm, t is the iteration number, and β is a constant (usually chosen as 0.2).

Influence Functions in CAFAC.
The belief space can influence the evolution in the population space in three ways: (i) determining the step size of the evolution, (ii) determining the direction of the evolution, (iii) determining the visual distance of AFA.
More precisely, if the normative knowledge is used to determine the step size of the evolution and visual distance in AFA, our knowledge-based AFA is named as CAFAC (Ns).In four versions of CAFAC, all of the behaviours, preying, swarming, and chasing, are modelled by the knowledge.The influence function for the CAFAC is defined as in Tables 3, 4 and 5.
If the situational knowledge is used to guide the direction of the evolution, our knowledge-based AFA is named as CAFAC (Sd).
If the normative knowledge guides the step size and the visual distance meanwhile the situational knowledge is used to determine the direction of the evolution, respectively, our knowledge-based AFA is named as CAFAC (Ns+Sd).
If the normative knowledge is used to determine the step size and direction of the evolution and the visual distance, our knowledge-based AFA is named as CAFAC (Ns+Nd).
In Tables 3, 4, and 5, size(I k ) = u k − l k is the size of the belief interval which is decided by the normative knowledge for the kth variable.The rand(0, 1) is a random number uniformly distributed in the interval (0,1).The other parameters are depicted in Section 2.2.

Crossover Operator.
A criterion is set up to judge whether the algorithm falls into local optimum: When the criterion is satisfied, the crossover operator will be applied to the ith artificial fish X i (i = 1, . . ., N ): where x r2 , x r1 are two individuals selected randomly, r 1 , r 2 are integers, which are generated randomly in the interval of [1, N ] and Evaluate the child x i and replace the individual x i with the child if x i performs better [35].

Application for Parameter Identification of Induction Motor
4.1.Identification Using CAFAC and PEM.In the identification process we disturb the system by voltage v and force f ex excitations.In order to decrease the dimension of the problem we remove the disturbance from the data.So the input f ex is set to zero.Then the problem to be considered is to match the model output with processed data using the recorded band limited white noise signal as input v.The basic idea of parameter estimation is to find a vector P to minimize the following cost function: where u rc (n) are the measurement of output for the real system under certain inputs; u rc (n) is the estimated output by exciting the model ( 1)-( 2) using the same inputs.The experiments are carried out using a simulated reference signal generated by the real system's parameters, which have been verified by simulations for a small 30 kW two-pole cage induction motor.
Here we use an indicator of fit value to evaluate the accuracy of an identified model for a specific parameters vector P and given inputs v(t) and f ex (t).The fit value is computed as [30] Table 2: Formulation for normative knowledge updating.
where the ith individual affects the lower bound for variable j, and the kth individual affects the upper bound for variable j.Note that t denotes the current generation of the belief space.
Table 3: Influence function for swarming.

Swarming
Ns where u rc is the average value of the measurement data {u rc (n)}.This indicator was chosen because it has been used in the research of optimal control of an induction motor using the prediction error method (PEM) for identification.
In our previous work, we found that the NsSd version has better performance than the three others we mentioned in Section 3 [26].Even if the CAFAC has better optimization function than PSO, the speed of the PSO is an obvious advantage.Therefore we combine the CAFAC and the basic PSO so as to speed up the optimization process.The optimization process is implemented using the software MATAB.The basic procedure of CAFAC algorithm can be described as follows.
(1) Set all the values for the parameters and initialize the N artificial fish in the search scope with random positions.
(2) Evaluate all the artificial fishes using the fitness function y and initialize the belief space.
(3) For each ith artificial fish, simulate the preying pattern, swarming, and chasing patterns separately, and select the best child fish.If the child is better, replace the ith artificial fish with the child.
(5) If the crossover criterion is satisfied, apply the crossover operator to the ith artificial fish got from Step (3).
(6) Switch to PSO until the termination criterion is satisfied.
(8) End the program if the final termination criterion is satisfied.
The termination criterion here is the same as (22) and the final termination criterion is the maximum of the whole iteration.
In PEM identification the general idea is to produce a model that minimizes a norm, such as where ε F (t, θ) is the prediction error filtered through a linear filter and (•) is a scalar function [36].Z N is the data consisting of N samples and θ is the model parameterization (i.e., polynomials).The problem is usually solved iteratively using gradient-based methods.
With the identification data, the PEM identification will result in fit values of 79.1 and 80.6 in horizontal and vertical directions.The CAFAC identification produced a model that had fit values of 73.9 and 77.7. Figure 4 displays the filtered measurement, and the model outputs from PEM and the parametric model.It can be stated that the models estimate the data reasonably well.The main differences occur near zero where the response is modest.

Simulation
Results.The obtained model was tested in simulations by implementing a state observer and a linear quadratic controller (LQ) for the model.The same was done in [29] with a model, which was identified using PEM.A real measurement of rotor vibrations (similar to  the identification data, but including the disturbance and not input signal) in the 2-pole induction motor was used as the input in the simulations.The expected result was that the observer would follow this signal and the LQ controller, fed through the model of the actuator and would compensate the vibrations.Figure 5 shows the estimation error of the observer.It can be concluded that the observer works decently since the amplitude of the estimation error is about 5 μm and the measured amplitude varies between 130 and 240 μm.The control result is displayed in Figure 6, where the controller is switched on after three seconds.During the 17-second period when control is on, the reductions of the vibrations are up to 67.0% and 70.3% and the control also reduces the fluctuation of the amplitude, which depends on the slip angle.
It is also necessary to consider the control signal because the test equipment can produce voltages in the range of ±100 V, and naturally the voltage levels should not be very high in a properly working controller.Figure 7 shows that there is at first a peak in the control signal, which reaches the saturation level, but it settles close to ±20 V range after 0.02 seconds.
The same simulations were also done with the black box model, which was identified using PEM, and the controller was tuned in the same manner.In this case the controller reduced vibrations in the simulations by 70.0% and 64.2%.The result is displayed in Figure 8.The control signals with the PEM model stay between ±50 V after that (see Figure 9).Strictly based on this data, the parametric model is slightly better than the PEM identification.However, the parametric model has not yet been successfully used in the real environment, whereas the PEM identification has been proven to work and to provide a reduction of the vibrations up to 90% depending on how the controller is tuned.

Conclusions
In this paper, we used a knowledge-based Artificial Fish-Swarm optimization algorithm to identify the parameters of an actuator model in an electrical machine.The culture framework was invested to direct the crossover operation in the AFA.In the culture framework, the situation knowledge Time (s) x-direction y-direction considerably.The results compared well with those obtained by PEM identification, but because the latter does not provide physical parameter values, the method presented in the paper can be considered an improvement.Also, the natureinspired evolutionary algorithms considered are particularly interesting, because they are computationally effective and reasonably easy to understand.Also, they do not need computation of gradients.For the future works, a timevarying model must be considered.In order to do that, more knowledge in the culture framework should be invested and more swarm intelligence must be tested to succeed in parameter identification.

Figure 2 :
Figure 2: (a) The test motor [28].Eddy current sensors are positioned on top of the rotor and behind it on the left.Encoder is located on the left side of the machine in the end of the rotor.(b) Control windings are built in the stator slots inside the motor [28].

Figure 3 :
Figure 3: The response to control voltages in xand y-directions.

Figure 4 :
Figure 4: Measured data and modelled outputs using PEM and CAFAC identification methods.

Figure 5 :Figure 6 :
Figure 5: The estimation error of the observer in xand y-directions.

Figure 7 :Figure 8 :Figure 9 :
Figure 7: Control voltage of the LQ controller with CAFAC model when control is switched off after 3 seconds.

Table 1 :
Parameters in the electromechanical model.

Table 4 :
Influence function for preying.

Table 5 :
Influence function for chasing.