Open access

Performances of the PCA Method in Electrical Machines Diagnosis Using Matlab

Written By

Jacques Fanjason Ramahaleomiarantsoa, Eric Jean Roy Sambatra, Nicolas Héraud and Jean Marie Razafimahenina

Submitted: 26 March 2012 Published: 26 September 2012

DOI: 10.5772/46449

From the Edited Volume

MATLAB - A Fundamental Tool for Scientific Computing and Engineering Applications - Volume 1

Edited by Vasilios N. Katsikis

Chapter metrics overview

5,408 Chapter Downloads

View Full Metrics

1. Introduction

Nowadays, faults diagnosis is almost an inevitable step to be maintained in the optimal safety operating of every physical system. Electrical machines, main elements of every electromechanical system, are among the research topics of many academic and industrial laboratories because of the importance of their roles in the industrial process. Lots of technologies of these machines are old and well controlled. However, they remain the seat of several electrical and mechanical faults [1-4]. Thus, this article deals with faults detection of a wound rotor synchronous machine (WRIM) by the principal component analysis (PCA) method.

Several diagnostic methods have been proposed and used in the literature for the electrical machines diagnosis [1-4]. The PCA method, which showed his effectiveness in the fault detection and isolation (FDI), was implemented recently for the system diagnosis [5-8].

This work is then to prove the strength of PCA method in faults diagnosis of systems using WRIM as application device.

To proceed with, in the first, we propose an accurate analytical model of the WRIM without or in the presence of faults [1, 9]. This model provides the matrix data of several characteristic quantities of the machine. These data will be included as input variables of the PCA method.

Then, we present a complete approach of PCA method based on the study of residues [10]. Special attention has been paid for the choice of the number of principal components to be maintained [11, 12].

These models are then implemented in the Matlab software. Simulation results of several variables (stator and rotor currents, shaft rotational speed, electrical power, electromagnetic torque and other variables issued from mathematical transformations) of healthy and faulted WRIM are analyzed. Comparisons of simulation results with those of other diagnostic methods are performed to show the effectiveness and importance of the PCA method in fault diagnosis systems [9].

The following are the different steps of the approach:

Figure 1.

Synoptic diagram of the different steps of the data treatment

The Figure 1 shows that the proposed approach is divided in four blocs:

  • WRIM modeling: mathematical equations calculation and simulation.

  • Simulation :graph showing the output states of the system (healthy and faulted operation)

  • Results Analysis: system diagnosis.

  • PCA: data treatment.

Advertisement

2. WRIM modeling

In the process of faults survey and diagnosis, an accurate modeling of the machine is necessary. In this paper, three phases model based on magnetically coupled electrical circuits were chosen.

The aim of the modeling is to highlight the electrical faults influences on the different state variables of the WRIM. For that, some modeling assumptions given in the following section are necessary.

2.1. Modeling assumptions

In the proposed approach, we assumed that:

  • the magnetic circuit is linear, and the relative permeability of iron is very large compared to the vacuum,

  • the skin effect is neglected,

  • hysteresis and eddy currents are neglected,

  • the airgap thickness is uniform,

  • magnetomotive force created by the stator and the rotor windings is sinusoidal distribution along the airgap,

  • the stator and the rotor have the same number of turns in series per phase,

  • the coils have the same properties,

  • the WRIM stator and rotor coils are coupled in star configuration and connected to the considered balanced state grid.

2.2. Differential equation system of the WRIM

We note that the voltage vectors ([VS], [VR]), the current vectors ([IS], [IR]) and the flux vectors ([ΦS], [ΦR]) of the stator and rotor are respectively: [VS]=[VAVBVC];[IS]=[IAIBIC]; [ϕS]=[ϕAϕBϕC][VR]=[VaVbVc]; [IR]=[IaIbIc]; [ϕR]=[ϕaϕbϕc]

Figure 2.

Equivalent electrical circuit of the WRIM

Vj, Ij and Φj (j : A, B, C for the stator phases and a, b, c, for the rotor phases) are respectively the voltages, the electrical currents and the magnetic flux of the stator and the rotor phases, θ is the angular position of the rotor relative to the stator.

The Figure 2 shows the equivalent electrical circuit of the WRIM. Each coil, for both stator and rotor, is modelised with a resistance and an inductance connected in series configuration (Figure 3).

Figure 3.

Equivalent electrical circuit of the WRIM coils

[VS]=[RS] [IS]+d[ϕS]dtE1
[VR]=[RR] [IR]+d[ϕR]dtE2
[ϕS]=[LS] [IS]+[MSR] [IR]E3
[ϕR]=[LR] [IR]+[MRS] [IS]E4

[RS] and [RR] are the resistance matrices, [LS] and [LR] the self inductance matrices, and [MSR] and [MRS] the mutual inductances matrix between the stator and the rotor coils.

With (3) and (4), (1) and (2) become:

[VS]=[RS]  [IS]+d{[LS] [IS]}dt+d{[MSR] [IR]}dtE5
[VR]=[RR]  [IR]+d{[LR] [IR]}dt+d{[MRS] [IS]}dtE6

By applying the fundamental principle of dynamics to the rotor, the mechanical motion equation is [13]:

JtdΩdt+fvΩ=CemCrE7
Ω=dθdtE8

with:

Cem=12[I]t*d( [L] )dθ*[I]E9

is the total inertia brought to the rotor shaft, the shaft rotational speed, [I]=[IA IB IC Ia Ib Ic]t the current vectors, the viscous friction torque, the electromagnetic torque, the load torque, the angular position of the rotor relative to the stator and [L] the inductance matrix of the machine.

Introducing the cyclic inductances of the stator and the rotorLSC=32LS and LRC=32LR(LS is the self inductance of each phase of the stator and LR is the self inductance of each phase of the rotor), the mutual inductances between the stator and the rotor coils MSR and pole pair number p, the inductance matrix of the WRIM can be written as follow:

[L]=[LSC00MSRf1MSRf2MSRf30LSC0MSRf3MSRf1MSRf200LSCMSRf2MSRf3MSRf1MSRf1MSRf3MSRf2LRC00MSRf2MSRf1MSRf30LRC0MSRf3MSRf2MSRf100LRC]E10
f1=cos(pθ)E11
f2=cos(pθ+2π3)E12
f3=cos(pθ2π3)E13

By choosing the stator and rotor currents, the shaft rotational speed and the angular position of the rotor relative to the stator as state variables, the differential equation system modeling the WRIM is given by:

[X˙]=[A]1([U][B][X])E14

with: [X]=[IA  IB  IC  Ia  Ib  Ic  Ω  θ]t[A]=[[L]000Jt0001];[U]=[[V]Cr0]; [V]=[VA  VB  VC  Va  Vb  Vc]t;

[B]=[[R]+Ωd[L]dθ0012[I]td[L]dθfv0010]E15

This model of the WRIM will be used to simulate both the healthy and the faulted configuration of the stator and the rotor.

2.3. WRIM faults

Despite the constant improvements on technical design of reliable machine, different types of faults still exist. The faults can be resulted by normal wear, poor design, poor assembly (misalignment), improper use or combination of these different causes.

Figure 4 and Figure 5 show the faults distribution carried out by a German company on industrial systems. The Figure 4 shows the faults of the low and medium power machines (50KW to 200KW), and the Figure 5 those of the high power machines (from 200KW) [1-2].

Figure 4.

Low and medium power induction machines faults [1-2]

Figure 5.

High power induction machines faults [1-2]

Figure 4 shows that the most encountered faults of the low and medium power on the induction machines are the stator faults and the Figure 5 shows that the faults due to mechanical defects give the highest percentage. The induction machine faults can be classified into four categories [2]:

The stator faults can be found on the coils or the breech. In most cases, the winding failure is caused by the inter-turns faults. These last grow and cause different faults between coils, between several phases or between phase and earth point before the deterioration of the machine [3]. The breech of electrical machines is built with insulated thin steel sheets in order to minimize the eddy currents for a greater operational efficiency. For the medium and great power machines, the core is compressed before the steel sheets emplacement to minimize the rolling sheets vibrations and to maximize the thermal conduction. The core problems are very low, only 1% if compared to winding problems [4].

The rotor faults can be bar breaks, coils faults or rotor eccentricities.

The bearings faults can be caused by a poor choice of materials during the manufacturing steps, the problems of rotation within the breech caused by damaged, chipped or cracked bearing and can create disturbance within the machines.

The other defaults might be caused by flange or shaft defaults. The faults created by the machine flange are generally caused during the manufacturing step.

2.4. Considered faults

The considered faults are on the resistance values which increase due to the rise of their temperature. In normal operation, a resistance value variation compared to its nominal value (in ambient temperature, 25°C) is a faulted machine due to machine overload or coils fault [1,9]. The resistance versus the temperature is expressed as:

R=R0(1+αΔT)E16

R0 is the resistance value at T0 = 25°C, α the temperature coefficient of the resistance and ΔT the temperature variation.

Advertisement

3. PCA methodology

The PCA method is based on simple linear algebra. It can be used as exploring tool, analyzing data and models design. The PCA method is based on the transformation of the data space representation. The new space dimension is smaller than that the original space dimension. It is classified as without model method, [5] and it can be considered as full identification method of physical systems [6]. The PCA method allows providing directly the redundancy relations between the variables without identifying the state representation matrix of the system. This task is often difficult to achieve.

3.1. PCA method formulation

We note by xi(j) = [x1 x2 x3 …xm] the measurements vector. « i » represents the measurement variables that must be monitored (i = 1 to m) and « j » the number of measurements for each variable « m », j = 1 to N.

The measurements data matrix (XdRN*m) can be written as follows:

Xd=(x1(1)...xm(1).........x1(N)...xm(N))E17

The data matrix is described by a smallest new matrix, that is an orthogonal linear projection of a subspace of m dimension on a less dimension subspace l (l<m). The method consists in identifying the PCA model and it is based on two steps [10]:

  • Determination on the eigenvalues and the eigenvectors of the covariance matrix R.

  • Determination of the structure of the model, which consists in calculating the component number « l » to be retained in the PCA model.

3.2. Eigenvalues and eigenvectors determination

The first step is the data normalization. The variables must be centered and reduced. Then, the obtained normalized matrix is:

X=[X1...Xm]E18

And the covariance matrix R is given by:

R=1N1XXTE19

By decomposing R, (17) can be expressed as:

R=PΛPTE20

With

PPT=PTP=ImE21

Λis the diagonal matrix of the eigenvectors of R and their eigenvalues are ordered in descending order with respect to magnitude values (λ1λ2...λm).

The eigenvectors matrix P is expressed as:

P=[p1,p2,...,pm]E22

piis the orthogonal eigenvectors corresponding to the eigenvalueλi. Then, the principal components matrix can be calculated using:

T=XP
TN*mE23

3.3. PCA model construction

To obtain the structure of the model, the components number « l » to be retained must be determined. This step is very important for PCA construction. The component number can be determined by using the following:

(i=1lλik=1mλk)*100thc,l<mE24

Where “thc” is an user defined threshold expressed as percentage. Now, user should retain only the components number « l » which was associated in the first term of (23). By reordering the eigenvalues, the minimum numbers of components are retained while still reaching the minimum variance threshold [14].

By taking into account the number of components to be retained and by partitioning the principal component matrix T, the eigenvectors matrix P and the eigenvalues matrix Λ[12], the constructed PCA model is given by:

T=[TpN*lTrN*(ml)]E25
P=[PpN*lPrN*(ml)]E26
Λ=[Λl*l00Λ(ml)(ml)]E27

Tpand Tr are respectively the principal and residual parts of T, Ppand Pr are respectively the principal and residual parts of P.

With this PCA model, the centered and reduced matrix X can be written as:

X=PpTpT+PrTrTE28

By considering:

Xp=PpTpT=i=1lPiTiTE29
E=PrTrT=i=l+1mPiTiTE30

The centered and reduced data matrix is given by:

X=Xp+EE31

Xpis the principal estimated matrix and E the residues matrix which represents information losses due to data matrix X reduction. It represents the difference between the exact and the approached representations of X. This matrix is associated with the lowest eigenvaluesλl+1,...,λm. Therefore, in this case, the data compression preserves all the best information that it conveys.

Advertisement

4. PCA method application on WRIM

4.1. Simulation conditions

Nine state variables (m=9) have been chosen to be monitored and 10000 measures (N=10000) during 4s are considered. The WRIM faults are introduced from the initial time (t=0s) to the final time (t=4s) of the different simulations. The machine is coupled to a mechanical load torque (10Nm) at t=2s. The considered faults are respectively, increases from 10% to 40% of the resistance value of both stator and rotor coils.

4.2. Choice of the number of principal components

The Figure 6 and the Figure 7 represent the residues variation of the WRIM stator current versus time and show impact of the « l » number in the diagnosis approach.

Figure 6.

Stator current residue for l = 5

Figure 6 show that the chosen number of components is too high then the residual space dimension is reduced. Some faults are projected in the principal space and the stator current residues can not be detectable.

However, with the Figure 6, the number of components is well chosen. Faults can be detected and localized and the PCA model is well reconstructed.

Generally, the detection approach in the case of diagnosis based on analytical model is linked with the residues generation step. From these residues analysis, the decision making step must indicate if faults exist or not. The residues generation approach can be the state estimation approach or the parameter estimation approach.

The residue indicates the information losses given by the matrix dimension reduction of the state variables matrix data to be monitored. Indeed, a small residue means that the estimated value tends to approach the exact value in healthy operation case.

In our case, the eigenvalues corresponding to the number of the retained principal components represent 93% of the total sum of eigenvalues. 0nly 7% of the total represent the residues subspace. One can conclude that the PCA model has been well constructed.

Figure 7.

Stator current residue for l = 6

Advertisement

5. Simulation results

The different simulation results have been performed with respect to the simulation conditions mentioned earlier.

Figure 8 to Figure 13 and Figure 16 represent the real variations without PCA method, and Figure 14, Figure 15 and Figure 17 represent the residue variations with PCA application of the faulted WRIM state variables in considering the stator defaults.

With the WRIM state variables, other quantities obtained by their transformations have been calculated:

  • quadrature axis and direct axis currents with Park transformation,

  • αβ

Figure 8.

Real variations versus time of the stator current of the healthy and faulted WRIM

Figure 9.

Real variations versus time of the rotor current of the healthy and faulted WRIM

Figure 10.

Real variations versus time of the shaft rotational speed of the healthy and faulted WRIM

Figure 11.

Real variations versus time of the electromagnetic torque of the healthy and faulted WRIM

Figure 12.

Real variations of β axis current versus the phase α axis current of the stator phase

Figure 13.

Real variations of the quadrature axis current versus the phase direct axis current of the stator phase

Figure 14.

Variations of the stator phase “A” current residues versus time of the healthy and faulted WRIM

Figure 15.

Variations of the rotor phase “a” current residues versus time of the healthy and faulted WRIM

Figure 16.

Real variations of electromagnetic torque versus the shaft rotational speed of the WRIM

Figure 17.

Variations of electromagnetic torque residues versus the shaft rotational speed residues of the WRIM

Advertisement

6. Discussion

Several types of representations are used in the signals processing domain, especially for the electrical machines diagnosis. We can mention the temporal representation (Figure 8 to Figure11, Figure 14 and Figure 15) and the signal frequency analysis. Although they have proved their efficiency, the state variable representations between them also show their advantages. They can be performed without mathematical transformation (Figure 16) and with mathematical transformation (Figure 12 and Figure 13).

The latter representation type and the temporal representation are confronted with the PCA method application results (Figure 14, Figure 15 and Figure 17). Only the simulation results with stator faults are presented because the global behavior of the state variables in both rotor and stator faults are almost similar.

For the temporal variations case, the rotor currents (Figure 9) and the shaft rotational speed (Figure 10) are the variables which produce the most information in presence of defaults. The defaults occur on the rotor current frequency and the shaft rotational speed magnitude.

The electromagnetic torque variations versus the shaft rotational speed also show clearly the WRIM operation zone in the presence of defaults (Figure 16). On the opposite, the representations with mathematical transformations (Figure 12 and Figure 13) do not provide significant information due to the fact that the stator currents remain almost unchanged in the presence of defaults (Figure 8).

With PCA method application, every representation type shows precisely the differences between healthy and faulted WRIM (Figure 14, Figure 15 and Figure 17). In the healthy case, residues are zero. When defaults appear, the residue representations have an effective value with an absolute value superior to zero.

In the figure 17, the healthy case is represented by a point situated on the coordinate origins. Therefore, one can show several right lines corresponding to the faulted cases. This behavior is due to the proportional characteristic of the considered faults.

PCA method proved to be very effective in electrical machines faults detection. This requires a good choice of the number of the principal components to be retained so that information contained in residues is relevant.

Advertisement

7. Implementation of the WRIM and the PCA models in the Matlab software

The differential equations system governing the WRIM is composed of linear differential equation which has the following form:

dydt=f(y)E32

Pre-programmed solvers are available in the Matlab software to solve easily this type of equation. These pre-programmed functions (ode45, ode113 …) proposed by the software helped to solve correctly with scalable computation time by the number of data to be processed. We adopted “ode45”, solver based on the Runge-Kutta 4, 5 numerical resolution method. After creating a function detailing the differential equation system, we have to use it in the chosen solver to calculate numerically the equations governing the WRIM.

The following extract lines of code illustrate the use of the pre-programmed function “ode45” (Figure 18):

Figure 18.

Use of ode45 solver to solve the differential equations system of the WRIM

The differential equations system is established in the “Diff_Equa_WRIM” function.

One of the major strengths of Matlab is the matrix manipulation. With the amount of data in matrix form that we consider in this paper, this feature of the software allows us to treat easily and without complexity these data. Thus, for the ACP method, matrix manipulations are done by simple operations because all variables in Matlab are intrinsically represented by matrix forms. In addition, pre-programmed functions are available to perform some precise operations such as the descendant sorting with “descend” function.

And finally, Matlab offers a multitude of possibilities for graphic representations. At the end of the PCA process, the original data and those from the treatment are represented graphically. This allowed more comparative studies as well as quantitative and qualitative analysis of the entire device. A function was reserved to the automatic superposition of curves of the same variables for the different considered defaults.

To summarize, three major functions have been developed to carry out the approach:

  • resolution of de differential equations system governing the WRIM,

  • resolution of the PCA method,

  • visualization of comparative curves.

We would like to note that each approach has been developed as a separate function, but our program runs automatically. These functions are executed automatically, one after the other, in another function.

Advertisement

8. Conclusion

PCA method based on residues analysis has been established and applied on WRIM diagnosis.

An accurate analytical model of the machine has been proposed and simulated to perform the healthy and faulted data for PCA approach need.

Several representations of nine state variables of the machine have been analyzed. For the temporal variation without PCA, the rotor current and the shaft rotational speed are the more affected by the considered fault type. The representations of the electromagnetic torque versus the shaft rotational speed in both with and without PCA approach show clearly the presence of defaults. Indeed, PCA method is interesting for all types of representation compared to some other signal processing types.

Simulation results show the efficiency of the detection but require a good choice of the number of principal components.

References

  1. 1. Benbouzid M. E. H2000A Review of Induction Motors Signature Analysis as a Medium for Faults Detection, IEEE Transactions On Industrial Electronics, 47N° 5, 984993
  2. 2. Chia-ChouY.al2008A reconfigurable motor for experimental emulation of stator winding interturn and broken bar faults in polyphase induction machines, 01, IEEE transactions on energy conversion, 23n°4. 10051014
  3. 3. SinM. L.SoongW. L.ErtugrulN.2003Induction machine on-line condition monitoring and fault diagnosis- a survey, University of Adelaide.
  4. 4. Negrea M.D2006Electromagnetic flux monitoring for detecting faults in electrical machines, PhD, Helsinki University of Technology, Department of Electrical and Communications Engineering, Laboratory of Electromechanics, TKK Dissertations 51, Espoo Finland.
  5. 5. LiuL.2006Robust fault detection and diagnosis for permanent magnet synchronous motors, PhD dissertation, College of Engineering, The Florida State University, USA.
  6. 6. BenaïchaA.GuerfeL. M.BouguilaN.BenothmanK.2010New PCA-based methodology for sensor fault detection and localization, 8th International Conference of Modeling and Simulation, MOSIM’10, “Evaluation and optimization of innovative production systems of goods and services”, Hammamet, Tunisia.
  7. 7. KuW.StorerR. H.GeorgakisC.1995Disturbance detection and isolation by dynamic principal component analysis. Chemometrics and Intelligent Laboratory Systems, 30. 179196
  8. 8. HuangB.2001Process identification based on last principal component analysis, Journal of Process Control, 11. 1933
  9. 9. Ramahaleomiarantsoa J. F, Héraud N, Sambatra E. J. R, Razafimahenina J. M2011Principal Components Analysis Method Application in Electrical Machines Diagnosis, ICINCO’08, Noordwykerhout, Pays Bas.
  10. 10. LiW.QinS. J.2001Consistent dynamic pca based on errors-in-variables Subspace identification, Journal of Process Control, 11(6). 661678
  11. 11. ValleS.WeihuaL.QinS. J.1999Selection of the number of principal components: The variance of the reconstruction error criterion with a comparison to other methods, Industrial and Engineering Chemistry Research, 3843894401
  12. 12. BenaïchaA.MourotG.GuerfelM.BenothmanK.RagotJ.2010A new method for determining PCA models for system diagnosis, 18th Mediterranean Conference on Control & Automation Congress Palace Hotel, Marrakech, Morocco June 23-25. 862867
  13. 13. WieczorekM.RosołowskiE.2010Modelling of induction motor for simulation of internal faults, Modern Electric Power Systems, MEPS’10, paper 29Wroclaw, Poland.
  14. 14. HalliganR. G.JagannathanS.2009PCA Based Fault isolation and prognosis with application to Water Pump, paper1, thesis, Fault detection and prediction with application to rotating machinery, Missouri University of Science and Technology.

Written By

Jacques Fanjason Ramahaleomiarantsoa, Eric Jean Roy Sambatra, Nicolas Héraud and Jean Marie Razafimahenina

Submitted: 26 March 2012 Published: 26 September 2012