Abstract

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.

1. 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 [24].

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 [69]. PSO has been widely used in identification in many research fields [1015]. 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 important 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 [2027]. 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.

2. 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 (𝑥 and 𝑦).

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]̈𝜂+2ΩΞ̇𝜂+Ω2=Φ𝑇𝑟𝑐𝑓𝑐+Φ𝑇𝑟𝑐𝑓ex,(1)𝑢𝑟𝑐=Φ𝑟𝑐𝜂,(2) where 𝜂 is the modal coordinate vector, 𝑢𝑟𝑐 is the displacement in 𝑥-direction and 𝑦-direction, Φ𝑟𝑐 is the modal matrix, Ω is a diagonal matrix, and Ξ is the modal damping matrix, 𝑓ex denotes the sinusoidal disturbance causing rotor vibration disturbance, and 𝑓𝑐 the electro mechanical force acting on the rotor. The model structure is the same as in [32]𝑑𝑑𝑡̇𝜂𝜂𝑖=2ΩΞΦ𝑇𝑟𝑐𝑃𝑒𝑚Φ𝑟𝑐Ω2Φ𝑇𝑟𝑐𝐶𝑒𝑚𝐼00𝑆𝑒𝑚Φ𝑟𝑐𝑄𝑒𝑚Φ𝑟𝑐𝐴𝑒𝑚̇𝜂𝜂𝑖+0Φ𝑇𝑟𝑐00𝐵𝑒𝑚0𝑣𝑓ex,𝑢𝑟𝑐=0Φ𝑟𝑐0̇𝜂𝜂𝑖.(3) The input 𝑣=(Re(̂𝑢𝑐,0),Im(̂𝑢𝑐,0)), is the control voltage, the first block row is the mechanical model, and 𝑖=(Re(̂𝑖𝑐,0),Re(𝑀𝑟,𝑐,2̂𝑖𝑟,2,0)T,Im(̂𝑖𝑐,0),Im(𝑀𝑟,𝑐,2̂𝑖𝑟,2,0))T, where ̂𝑖𝑐,0 and ̂𝑖𝑟,2,0 are current space vectors and 𝑀𝑟,𝑐,2 included in 𝑖 is mutual inductance of control winding and the rotor cage. In variables the subscript 𝑐 denotes a relation to the control winding, the subscript 𝑟 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 coordinateŝ𝑢𝑐=𝑅𝑐̂𝑖𝑐+𝑑𝜓𝑐,2𝑑𝑡,0=𝑅𝑟,2̂𝑖𝑟,2+𝑑𝜓𝑟,2𝑑𝑡𝑗2𝜔𝑚𝜓𝑟,2.(4) The four-pole harmonic flux leakages are𝜓𝑐,2=𝐿𝑐̂𝑖𝑐+𝑀𝑟,𝑐,2̂𝑖𝑟,2+𝑋𝑐,𝜖,2𝐵1𝑧𝑟,𝜓𝑟,2=𝐿𝑟,2̂𝑖𝑟,2+𝑀𝑐,𝑟,2̂𝑖𝑐+𝑋𝑟,𝜖,2𝐵1𝑧𝑟.(5) By substitution ̂𝑢𝑐,0=𝑅𝑐+𝑗𝜔1𝐿𝑐̂𝑖𝑐,0+𝐿𝑐𝑑̂𝑖𝑐,0𝑑𝑡+𝑗𝜔1𝑀𝑟,𝑐,2̂𝑖𝑟,2,0+𝑑𝑑𝑡𝑀𝑟,𝑐,2̂𝑖𝑟,2,0+𝑋𝑐,𝜖,2𝐵1̇𝑧𝑟+𝑗𝜔1𝑧𝑟,0=𝛾𝑟,2+𝑗𝜔2𝑀𝑟,𝑐,2̂𝑖𝑟,2,0+𝑑𝑑𝑡𝑀𝑟,𝑐,2̂𝑖𝑟,2,0+𝑗𝜔2𝑀𝑟,𝑐,2𝑀𝑐,𝑟,2̂𝑖𝑐,0+𝑀𝑟,𝑐,2𝑀𝑐,𝑟,2𝑑̂𝑖𝑐,0𝑑𝑡+𝑀𝑟,𝑐,2𝑋𝑟,𝜖,2𝐵1̇𝑧𝑟+𝑗𝜔2𝑧𝑟,(6) where 𝑧 is a complex coordinate representation of the rotor eccentricity. In addition, the control force 𝑓𝑐 in terms of the parameters in Table 1 is given by𝑓𝑐=𝜋𝑑𝑟𝑙𝑟4𝜇0𝛿0𝛼𝑧𝑟+𝛽𝑒𝑗2𝜔1𝑡𝑧𝑟𝑃𝑒𝑚(𝑡)+𝜋𝑑𝑟𝑙𝑟4𝛿0𝐵1𝑘𝑟,2𝑀𝑟,𝑐,2𝑀𝑟,𝑐,2̂𝑖𝑟,2,0+𝑘𝑐̂𝑖𝑐,0.𝐶𝑒𝑚(7) These equations include all the parameters to identify. Since the complex coordinate system of 𝑧 is decomposed in real and imaginary parts in the final model, matrices 𝑃𝑒𝑚(𝑡) and 𝐶𝑒𝑚, are given by𝑃𝑒𝑚(𝑡)=𝜋𝑑𝑟𝑙𝑟4𝜇0𝛿0Re{𝛼}Im{𝛼}Im{𝛼}Re{𝛼}+Re𝛽𝑒𝑗2𝜔1𝑡Im𝛽𝑒𝑗2𝜔1𝑡Im𝛽𝑒𝑗2𝜔1𝑡Re𝛽𝑒𝑗2𝜔1𝑡,(8)𝐶𝑒𝑚𝐶=𝜋𝑑𝑟𝑙𝑟4𝛿0𝑘𝑐𝐵1𝑘𝑟,2𝑀𝑟,𝑐,2𝐵1,(9)𝐶𝑒𝑚=Re𝐶𝑒𝑚𝐶Im𝐶𝑒𝑚𝐶Im𝐶𝑒𝑚𝐶Re𝐶𝑒𝑚𝐶.(10) It should be noted that in the time-invariant case 𝑃𝑒𝑚(𝑡) is constant 𝑃𝑒𝑚.

Equation (6) can be written in matrix form𝐿𝑐1𝑀𝑟,𝑐,2𝑀𝑐,𝑟,21×𝑑𝑑𝑡̂𝑖𝑐,0𝑀𝑟,𝑐,2̂𝑖𝑟,2,0+𝑅𝑐+𝑗𝜔1𝐿𝑐𝑗𝜔1𝑗𝜔2𝑀𝑟,𝑐,2𝑀𝑐,𝑟,2𝛾𝑟,2+𝑗𝜔2̂𝑖𝑐,0𝑀𝑟,𝑐,2̂𝑖𝑟,2,0+𝑋𝑐,𝜖,2𝐵1𝑀𝑟,𝑐,2𝑋𝑟,𝜖,2𝐵1̇𝑧𝑟+𝑗𝜔1𝑋𝑐,𝜖,2𝐵1𝑗𝜔1𝑀𝑟,𝑐,2𝑋𝑟,𝜖,2𝐵1𝑧𝑟=10̂𝑢𝑐,0.(11)

Rearranging this will give the rest of the matrices, which are needed for the combined model 𝑑𝑑𝑡̂𝑖𝑐,0𝑀𝑟,𝑐,2̂𝑖𝑟,2,0=𝐿𝑐1𝑀𝑟,𝑐,2𝑀𝑐,𝑟,211𝑋𝑐,𝜖,2𝐵1𝑀𝑟,𝑐,2𝑋𝑟,𝜖,2𝐵1𝑆𝑒𝑚𝐶̇𝑧𝑟𝐿𝑐1𝑀𝑟,𝑐,2𝑀𝑐,𝑟,211𝑗𝜔1𝑋𝑐,𝜖,2𝐵1𝑗𝜔1𝑀𝑟,𝑐,2𝑋𝑟,𝜖,2𝐵1𝑄𝑒𝑚𝐶𝑧𝑟𝐿𝑐1𝑀𝑟,𝑐,2𝑀𝑐,𝑟,211𝑅𝑐+𝑗𝜔1𝐿𝑐𝑗𝜔1𝑗𝜔2𝑀𝑟,𝑐,2𝑀𝑐,𝑟,2𝛾𝑟,2+𝑗𝜔2𝐴𝑒𝑚𝐶×̂𝑖𝑐,0𝑀𝑟,𝑐,2̂𝑖𝑟,2,0+𝐿𝑐1𝑀𝑟,𝑐,2𝑀𝑐,𝑟,21110𝐵𝑒𝑚𝐶̂𝑢𝑐,0.(12) Matrices 𝐴𝑒𝑚, 𝐵𝑒𝑚, 𝑆𝑒𝑚, and 𝑄𝑒𝑚 are finally obtained like 𝐶𝑒𝑚 in  (10)𝐴𝑒𝑚=Re𝐴𝑒𝑚𝐶Im𝐴𝑒𝑚𝐶Im𝐴𝑒𝑚𝐶Re𝐴𝑒𝑚𝐶,𝐵𝑒𝑚=Re𝐵𝑒𝑚𝐶Im𝐵𝑒𝑚𝐶Im𝐵𝑒𝑚𝐶Re𝐵𝑒𝑚𝐶,𝑆𝑒𝑚=Re𝑆𝑒𝑚𝐶Im𝑆𝑒𝑚𝐶Im𝑆𝑒𝑚𝐶Re𝑆𝑒𝑚𝐶,𝑄𝑒𝑚=Re𝑄𝑒𝑚𝐶Im𝑄𝑒𝑚𝐶Im𝑆𝑒𝑚𝐶Re𝑄𝑒𝑚𝐶.(13)

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 𝑢(𝑡) consists of the voltage input 𝑣 to the control winding and the disturand the outputs are the displacements 𝑢𝑟𝑐 of the rotor center. The system matrices 𝐴, 𝐵, and 𝐶 are the functions of the unknown parameters in the model which are reformed in a vector as𝑃=𝐿𝑐,𝑘𝑐𝐵1,𝑀𝑟,𝑐,2𝑀𝑐,𝑟,2,𝛾𝑟,2,𝑋𝑐,𝜖,2𝐵1,𝛼,𝑘𝑟,2𝐵1𝑀𝑟,𝑐,2,𝑋𝑐,𝜖,2𝐵1𝑀𝑟,𝑐,2.(14) The corresponding significances of the unknown parameters, and the detailed derivation of the model can be found in [32].

2.2. The Artificial Fish-Swarm Algorithm

Suppose that the problem under consideration has 𝐷-dimensions. Initialize the swarm with 𝑁 artificial fish. The state of one artificial fish can be formulated as 𝑋𝑖=(𝑥𝑖1,𝑥𝑖2,,𝑥𝑖𝐷), 𝑖=1,,𝑁, where 𝑋𝑖 represents the target variable for the problem under consideration. 𝑦=𝑓(𝑋𝑖) 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 (𝑑𝑖𝑗=𝑋𝑗𝑋𝑖 stands for the Euclidean distance between 𝑋𝑖 and 𝑋𝑗, “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 𝑋𝑖, the artificial fish selects a state 𝑋𝑗 randomly within the visual distance, such as 𝑋𝑗=𝑋𝑖+rand(0,1)×vd. If 𝑓(𝑋𝑗)<𝑓(𝑋𝑖) the artificial fish moves from 𝑋𝑖 towards𝑋𝑗, meaning that 𝑋𝑡𝑖𝑋𝑡+1𝑖. The formulation can be depicted as follows: 𝑋𝑡+1𝑖=𝑋𝑡𝑖+rand(0,1)×𝑠×𝑋𝑡𝑗𝑋𝑡𝑖𝑋𝑡𝑗𝑋𝑡𝑖.(15) If 𝑓(𝑋𝑗)>𝑓(𝑋𝑖), 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 𝑋𝑡+1𝑖=𝑋𝑡𝑖+rand(0,1)×𝑠.(16)

(2) Swarming
The current state of the artificial fish is 𝑋𝑖,  nf is the number of its fellows within the visual distance, equal to the number of elements in the set of 𝐵={𝑋𝑗𝑑𝑖𝑗Visual}. If nf0 that is the set 𝐵 is not empty, let 𝑋center=nf𝑗=1𝑋𝑗/nf and then 𝑦center=𝑓(𝑋center) stands for the fitness of the centre position. If nf×𝑦center<𝛿×𝑦𝑖, meaning that this area is not crowded, then if 𝑦center<𝑦𝑖, the artificial fish moves one step towards the centre position: 𝑋𝑡+1𝑖=𝑋𝑡𝑖+rand(0,1)×𝑠×𝑋𝑡𝑐𝑋𝑡𝑖𝑋𝑡𝑐𝑋𝑡𝑖.(17) Otherwise it executes the behaviour of preying.

(3) Chasing
The current state of the artificial fish is 𝑋𝑖, 𝑋minstands for the best artificial fish individual within 𝑋𝑖’s visual distance. nf is the number of 𝑋min’s fellows within the visual distance. 𝑦min=𝑓(𝑋min), if 𝑦min<𝑦𝑖 and nf×𝑦min<𝛿×𝑦𝑖, the artificial fish moves one step to 𝑋min: 𝑋𝑡+1𝑖=𝑋𝑡𝑖+rand(0,1)×𝑠×𝑋𝑡min𝑋𝑡𝑖𝑋𝑡min𝑋𝑡𝑖.(18) Otherwise it executes the behaviour of preying.

3. 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.

3.1. 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, 𝑆=𝑆𝑡𝑆𝑡={𝑠𝑡1,𝑠𝑡2,,𝑠𝑡𝑛}, 𝑆𝑡 stands for the best fish in the swarm at iteration 𝑡. In other words, it can be initialized with the best fish in the initial fish swarm and updated by the following function: 𝑠𝑡+1𝑗=𝑋𝑡+1𝑔best,𝑗if𝑓𝑋𝑡+1𝑔best<𝑓𝑠𝑡,𝑠𝑡𝑗otherwise,(19) where 𝑋𝑡+1𝑔best denotes the best artificial fish individual in the fish swarm at generation 𝑡+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 𝑁=𝐼,𝑈,𝐿,𝐷,(20) where 𝑈, 𝐿 and 𝐷 are 𝑛-dimensional vectors, and 𝐼={𝑥𝑙𝑥𝑢}, 𝑛 is the number of the variables, 𝑙𝑗 and 𝑢𝑗 are the lower and upper bounds for the 𝑗th variable, respectively, 𝐿𝑗 and 𝑈𝑗are the values of the fitness function associated with the bound 𝑙𝑗 and 𝑢𝑗. Generally, 𝑙𝑗 and 𝑢𝑗 are initialized with the lower and upper bounds of individuals. 𝐿𝑗 and 𝑈𝑗 are usually initialized with positive infinity. The formulation for normative knowledge updating can be depicted in Table 2.

3.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]: 𝑓𝑎(𝑁,𝑡)=𝑁𝛽+𝑁𝛽𝑡,(21) where 𝑁 is the size of the swarm, 𝑡 is the iteration number, and 𝛽 is a constant (usually chosen as 0.2).

3.3. 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(𝐼𝑘)=𝑢𝑘𝑙𝑘 is the size of the belief interval which is decided by the normative knowledge for the 𝑘th 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.

3.4. Crossover Operator

A criterion is set up to judge whether the algorithm falls into local optimum:|||||𝑓𝑋𝑡𝑖𝑓𝑋𝑡1𝑖𝑓𝑋𝑡1𝑖|||||<0.1.(22) When the criterion is satisfied, the crossover operator will be applied to the 𝑖th artificial fish 𝑋𝑖(𝑖=1,,𝑁):𝑥𝑖=𝑥𝑟1+𝛼×𝑥𝑟2𝑥𝑟1,(23) where 𝑥𝑟2, 𝑥𝑟1 are two individuals selected randomly, 𝑟1, 𝑟2 are integers, which are generated randomly in the interval of [1,𝑁] and 𝑟1𝑟2𝑖.  𝛼 is random number uniformly distributed in interval of [𝑑,1+𝑑] and 𝑑 is constant (0.25). Evaluate the child 𝑥𝑖 and replace the individual 𝑥𝑖 with the child if 𝑥𝑖 performs better [35].

4. Application for Parameter Identification of Induction Motor

4.1. Identification Using CAFAC and PEM

In the identification process we disturb the system by voltage 𝑣 and force 𝑓ex excitations. In order to decrease the dimension of the problem we remove the disturbance from the data. So the input 𝑓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 𝑣. The basic idea of parameter estimation is to find a vector 𝑃 to minimize the following cost function:𝐽(𝑃)=𝑁𝑛=1|||𝑢𝑟𝑐(𝑛)𝑢𝑟𝑐(𝑛)|||𝑁,(24) where 𝑢𝑟𝑐(𝑛) are the measurement of output for the real system under certain inputs; 𝑢𝑟𝑐(𝑛) 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 𝑃 and given inputs 𝑣(𝑡) and 𝑓ex(𝑡). The fit value is computed as [30]𝑚𝑑=100×1𝐽(𝑃)𝑁𝑛=1|||𝑢𝑟𝑐(𝑛)𝑢𝑟𝑐|||21/2,(25) where 𝑢𝑟𝑐 is the average value of the measurement data {𝑢𝑟𝑐(𝑛)}. 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 𝑁 artificial fish in the search scope with random positions.(2)Evaluate all the artificial fishes using the fitness function 𝑦 and initialize the belief space.(3)For each 𝑖th artificial fish, simulate the preying pattern, swarming, and chasing patterns separately, and select the best child fish. If the child is better, replace the 𝑖th artificial fish with the child.(4)Update the belief space.(5)If the crossover criterion is satisfied, apply the crossover operator to the 𝑖th artificial fish got from Step (3).(6)Switch to PSO until the termination criterion is satisfied.(7)Switch to (3) 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𝑉𝑁𝜃,𝑍𝑁=1𝑁𝑁𝑡=0𝜀𝐹(𝑡,𝜃),(26) where 𝜀𝐹(𝑡,𝜃) is the prediction error filtered through a linear filter and () is a scalar function [36]. 𝑍𝑁 is the data consisting of 𝑁 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.

4.2. 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.

5. 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 and the normative knowledge were employed to guide the evolution of the Artificial Fish-Swarm optimization. The crossover operation can help the artificial fish jump out of the local optimum without losing the characteristics of the previous generation. The proposed knowledge-based Artificial Fish-Swarm optimization can improve the performance of the original Artificial Fish-Swarm optimization and can be applied, for example, to find parameter values for a model of an actuator used for vibration control of rotor in an induction motor. Realistic values for the components of a structural first-principles electromechanical model were obtained, which improved the earlier identification results 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 nature-inspired 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 time-varying 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.

Acknowledgment

This research work was funded by the Academy of Finland under Grants 135225 and 127299 and the NSFC under Grant no. 60874084.