Next Article in Journal
A Parameter Self-Calibration Method for GNSS/INS Deeply Coupled Navigation Systems in Highly Dynamic Environments
Next Article in Special Issue
Software Sensor for Activity-Time Monitoring and Fault Detection in Production Lines
Previous Article in Journal
Extending QGroundControl for Automated Mission Planning of UAVs
Previous Article in Special Issue
Fault Detection and Isolation via the Interacting Multiple Model Approach Applied to Drive-By-Wire Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Partial Inductance Model of Induction Machines for Fault Diagnosis

by
Manuel Pineda-Sanchez
,
Ruben Puche-Panadero
*,
Javier Martinez-Roman
,
Angel Sapena-Bano
,
Martin Riera-Guasp
and
Juan Perez-Cruz
Institute for Energy Engineering, Universitat Politècnica de València, Cmno. de Vera s/n, 46022 Valencia, Spain
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(7), 2340; https://doi.org/10.3390/s18072340
Submission received: 31 May 2018 / Revised: 12 July 2018 / Accepted: 16 July 2018 / Published: 18 July 2018
(This article belongs to the Special Issue Sensors for Fault Detection)

Abstract

:
The development of advanced fault diagnostic systems for induction machines through the stator current requires accurate and fast models that can simulate the machine under faulty conditions, both in steady-state and in transient regime. These models are far more complex than the models used for healthy machines, because one of the effect of the faults is to change the winding configurations (broken bar faults, rotor asymmetries, and inter-turn short circuits) or the magnetic circuit (eccentricity and bearing faults). This produces a change of the self and mutual phase inductances, which induces in the stator currents the characteristic fault harmonics used to detect and to quantify the fault. The development of a machine model that can reflect these changes is a challenging task, which is addressed in this work with a novel approach, based on the concept of partial inductances. Instead of developing the machine model based on the phases’ coils, it is developed using the partial inductance of a single conductor, obtained through the magnetic vector potential, and combining the partial inductances of all the conductors with a fast Fourier transform for obtaining the phases’ inductances. The proposed method is validated using a commercial induction motor with forced broken bars.

1. Introduction

Induction machines (IMs) are key components of modern industrial processes, and their unexpected failures can generate unscheduled stops and high economic losses. To avoid this risk, on-line condition monitoring of IMs has become an integral part of industrial maintenance systems [1,2]. A great research effort has been devoted in the last decades to the development of fault diagnostic systems that can detect multiple machine faults in an early stage, giving a reliable information about the machine condition, and avoiding false alarms [3,4]. Different signals have been proposed to build fault diagnosis systems, such as vibrations [5], stator currents [6], thermal images [7], acoustic signals [8,9], etc. In particular, diagnostic algorithms that can detect the faults through their characteristic signature in the stator currents, known as motor current signal analysis (MCSA), have gain a great relevance [6,10], because they can be applied on-line using a simple hardware and software, just current clamps to measure the stator currents and the fast Fourier transform (FFT) to analyze them in the frequency domain [11,12]. Multiple types of faults can be detected by identifying the characteristic harmonics that they generate in the stator current, with frequencies that have been established theoretically.
For example, the characteristic frequencies of the fault harmonic components are, for bar breakages in squirrel cage induction machines [12,13,14] and rotor asymmetries in wound-rotor induction machines f b b [15]
f b b = ( 1 + 2 k s ) f 1 k = ± 1 , ± 2 , ± 3
for mixed eccentricity faults f e c c [16,17,18]
f e c c = f 1 1 + k 1 s p k = ± 1 , ± 2 , ± 3
for inter-turn short-circuits in the stator winding f s c h [19,20,21,22]
f s c h = f 1 m 1 s p ± k m = 0 , 1 , 2 , , k = 0 , 1 , 3 ,
and for mechanical looseness f l o o s e [23]
f l o o s e = f 1 1 + k 1 s p · n k = 1 , 2 , , n = 2 , 3 ,
where s is the slip, f 1 is the frequency of the power supply and p is the pole pairs number. Other types of faults, such as bearing faults [24,25,26], or gearbox faults [27] generate components in the stator current with frequency expressions similar to Equations (1)–(4). These equations are valid both in steady-state regime and in transient regime, but in this case the frequencies given by Equations (1)–(4) change along the time [15], and advanced time–frequency transforms must be used instead of the FFT [25,28].
Commonly, fault diagnostic systems that are based on the analysis of the current signal, such as MCSA, try to measure the amplitude of the harmonic components with frequencies given by Equations (1)–(4) in the spectrum of the stator current. These values are then compared with predefined threshold alarms, in data based diagnostic systems [5], or with the output of machine models, in model based diagnostic systems [5,20,29,30,31]. In both cases, accurate machine models are needed, for fixing the values of the threshold alarms, for running the model-based fault diagnostic system [32], or for developing new diagnostic algorithms [33]. The most accurate machine models are built using the finite elements analysis (FEA) [34], because it can take into account the actual geometry of the machine, including the asymmetries induced by the fault (eccentric machines), the characteristic of the magnetic materials, which depend on their operating point (saturation, hysteresis), and the configuration of the windings, which can be highly asymmetrical in case of rotor asymmetries or inter-turn short-circuits faults [35]. However, FEA models demand a great amount of computational resources and computing time [17]. This is a serious drawback for using FEA models on embedded on-line diagnostic hardware [24], such as DSPs or FPGAs [36], which have limited computing resources, or in systems which require a fast response, such as real-time model-based diagnostic systems [37]. A growing trend is the use of the IM model in hardware-in-the-loop (HIL) systems, which are used for developing and testing real-time diagnostic hardware, [38], or for training advanced artificial intelligence (AI) based diagnostic systems, such as the set membership identification (SMI) approach [14], neural networks [39,40,41], support vector machines [42] or fuzzy inference systems [43].
To alleviate these drawbacks, several solutions have been proposed in the technical literature, such the use of coupled finite elements and circuit equations [44], improved conformal mappings [45], or field reconstruction methods (FRM), based on linear superposition of FEA solutions [46]. An alternative approach to FEA models are fast and simpler analytical models of the induction machine [47], such as the magnetic equivalent circuit (MEC) model [36,48], which represents the induction machine as a mesh of reluctance elements, or the permeance network model (PNM) [49]. Recent works also suggest the use of generalized two-phase model of an induction machine of fifth order [50], space phasor models [27], complex-vector models [51], or equivalent circuit models [52].
Among the analytical models, the multiple coupled circuit model (MCCM) [53,54], also called the winding function approach [55], has been extensively used for modeling the IM in faulty conditions [56,57,58]. As Ojaghi and Faiz [55] pointed out, it is the most detailed and complete model used to analyze the performance of cage IMs, especially under various faults. It relies on the time-varying self and mutual inductances of the machine phases and rotor-loops, and their derivatives, which are calculated by integral expressions representing the placement of the winding turns along the air-gap periphery. However, these integral expressions must be computed for every rotor position and phase winding layout, which can be altered due to the asymmetries generated by the machine’s fault, resulting in the need of heavy computing resources to build the machine inductances matrices under faulty conditions. As it is stated in [54], this task typically consumes a high amount of time, so that only discrete curves of inductance versus rotor position are calculated and linear interpolation is applied at intermediate rotor positions. To solve this problem, in this paper, a novel, fast approach is proposed for obtaining the matrices of self and mutual inductances, making use of the concept of partial inductance of a single conductor [59,60]. Instead of solving the integral expression for every different winding layout, the partial inductance of a single conductor [61] is computed, using the magnetic vector potential (MVP), and the total self and mutual inductances of the windings are obtained by simply combining the partial inductances of their conductors. This combination is made in a very fast way using the fast Fourier transform. The partial inductance approach has been used in the partial element equivalent circuit (PEEC) method [62] for the electromagnetic modeling of power electronic modules [63], antennas [64], high voltage transmission lines [65] or printed circuit boards (PCBs) [66], but it is the first time, to the authors’ best knowledge, that it is proposed for obtaining the phase inductances of IMs. Compared with the WFA, this novel method is characterized by the following main points:
  • The conductor, instead of the coil, is used as the basic winding unit, which simplifies the modeling of arbitrarily complex winding layouts.
  • The partial inductance of a single conductor, computed using the MVP, is taken as the characteristic function of the machine’s windings, instead of the winding functions based on the magnetomotive force (MMF) of the coils.
  • A fast circular convolution, computed through the FFT, is used to obtain the mutual and self inductances of the phases, instead of solving the integral expressions of the winding functions for every rotor position.
The structure of this paper is as follows. In Section 2, the electro-mechanical equations used for modeling the IM are briefly presented. In Section 3, the partial inductance of a single conductor is obtained through its MVP, and in Section 4 the FFT is used for assembling the matrices of self and mutual inductances of the machine, combining the partial inductances of the phases’s conductors. Section 5 presents the experimental validation of the proposed approach, using a commercial IM with broken bar faults. Finally, Section 6 presents the conclusions of this work.

2. Electro-Mechanical System of Equations of an Induction Machine

The following system of equations can be written for an induction machine with m stator and n rotor phases with arbitrary layout, that is, even under fault conditions, [67]
[ U s ] = [ R s ] [ I s ] + d [ Ψ s ] / d t
[ 0 ] = [ R r ] [ I r ] + d [ Ψ r ] / d t
[ Ψ s ] = [ L s s ] [ I s ] + [ L s r ] [ I r ]
[ Ψ r ] = [ L s r ] T [ I s ] + [ L r r ] [ I r ]
[ U s ] = [ u s 1 , u s 2 , , u s m ] T
[ I s ] = [ i s 1 , i s 2 , , i s m ] T , [ I r ] = [ i r 1 , i r 2 , , i r n ] T
where [ U ] is the phase voltages matrix, [ I ] is the phase currents matrix, [ R ] is the resistances matrix, [ Ψ ] is the flux linkages matrix and [ L ] is the inductances matrix. Subscripts s and r are used for the stator and for the rotor, respectively. The mechanical equations are
T e = [ I s ] T [ L s r ] φ [ I r ]
T e T L = J d Ω d t = J d 2 θ d t 2
where T e is the electromechanical torque of the machine, T L is the load torque, J is the total system inertia (rotor plus load), Ω is the mechanical speed, θ is the angular position of the rotor, and φ is the angular coordinate. Two additional constraints, namely the sum of all the stator currents equals 0, and the sum of all the rotor currents equals 0, must be added to this set of equations.
To solve the system of Equations (5)–(12), the inductance matrices and their derivatives must be calculated very accurately, in order to be able to reproduce the fault harmonics given by Equations (1)–(4). The inductances in Equations (7), (8) and (11) correspond to the self and mutual inductances of the machine’s phases, and of the rotor loops in the case of squirrel cage machines (two consecutive bars and the sections of the end-rings between them). The inductances of these windings, such as the ones represented in Figure 1, can be computed summing up the coil inductances for all the phases’ coils, as in WFA. However, this approach must deal with the wide variety of coil types that can be used for building a phase winding. Besides, these inductances may change with the rotor position, as in the case of rotor-stator mutual inductances, and also in the case of winding faults, such as inter-turn short-circuits, or bar breakages. This approach results in complex and lengthy computations that typically consume a high amount of time [54]. The use of the partial inductance concept can greatly simplify this process, because the partial inductance of the conductor, instead of the coil inductance, is chosen as the building block for the computation of the inductance matrices, which is performed using a very simple and efficient FFT-based procedure, even in the case of complex or faulty winding layouts.
In this work, the following assumptions will be made for calculating the phase inductances, common to many analytical IM models presented in the technical literature:
  • There is negligible saturation.
  • There is a uniform air-gap.
  • The eddy currents are neglected. The windings are considered to be made of filamentary conductors placed on the external rotor surface or on the internal stator surface.
  • The end-coil leakage inductances need to be pre-calculated, and are treated as constants in Equations (7), (8) and (11), using explicit expressions that can be found in [52,68].
Some of these assumptions can be removed by using a modified inverse air-gap function which takes into account the saturation [55], the rotor eccentricity [30], or the slotting effects [69]. Nevertheless, as the focus of this paper is the introduction of the partial inductance model of the IM, a constant air-gap has been used to present the proposed approach, applied to the simulation and experimental test of an IM with a rotor asymmetry (broken bars fault), the same fault analyzed in [11,12,14,15,56], among many others.

3. Phase Inductance as a Sum of the Conductors’ Partial Inductances

Let us consider a single coil of the machine windings, represented as a rectangular current loop in Figure 2, where conductors c 1 and c 2 are parallel to the machine’s axis, and conductors c 3 and c 4 belong to the end windings.
The self inductance of this single coil, under the assumptions given in Section 2, is given by Equation
L c o i l = Ψ I
where Ψ and I are the flux linked by the coil and its current, respectively, assuming that I is the only current in the machine. The same equation holds for the mutual inductance between two coils, just using the current of the second coil in Equation (13). However, the coil inductance given by Equation (13) is a global coil characteristic which can’t be assigned to any segment of the loop represented in Figure 2. On the contrary, a partial inductance can be assigned individually to each of the conductors that forms the coil loop. In effect, Equation (13) can be expressed as Equation
L c o i l = S B · d S I
where B is the magnetic flux density through the loop area S . Using B = × A , where A is the magnetic vector potential, and applying Stokes’ theorem, gives Equation
L c o i l = S × A · d S I = C A · d l I
where C is the contour of the current loop. However, the contour integral in Equation (15) can be broken into the four segments of the loop of Figure 2 as
L c o i l = C 1 A · d l I + C 2 A · d l I + C 3 A · d l I + C 4 A · d l I = L c o i l 1 + L c o i l 2 + L c o i l 3 + L c o i l 4
where C i is the contour of each segment of the loop. The net inductance of each segment C i , L i in Equation (16), represents the contribution of the segment C i to the total inductance of the current loop.
Using the concept of partial inductance, the inductance of each segment i of the coil in Equation (16) L c o i l i ( i = 1 4 ) can be expressed as the sum of the partial inductances of segment i in Equation
L c o i l i = C i A · d l I = j = 1 4 L p i j
where the partial inductance L p i j is defined [59,60] as the ratio between the magnetic flux crossing the rectangular surface between the segment i and infinity, and the current I j that produces that flux (which, in this case, is the same for all the loop conductors). That is,
L p i j = C i A ij · d l I j
Using Equations (16) and (17), the total inductance of the rectangular loop in Equation (13) can be expressed as the sum of all the partial inductances of its segments, as
L c o i l = i = 1 4 j = 1 4 L p i j
As stated in the previous section, the end-coil leakage inductance is considered in this work as a constant value, computed a priori, so Equation (19) must be applied only to the axial conductors of the coil, C 1 and C 2 in Figure 2. With this assumption, the sums in Equation (19) must be extended only to these two conductors, as
L c o i l = i = 1 2 j = 1 2 L p i j + L e n d _ c o i l
In a similar way, the self-inductance of a phase A with N A conductors can be computed extending Equation (20) to account for all the phase’s conductors, as
L p h a s e A = i = 1 N A j = 1 N A L p i j + L e n d _ p h a s e
where L e n d _ p h a s e represents the total end-winding leakage inductance of the phase.
The expression for the mutual inductance between phase A and other phase B, with N B conductors, is analogous to Equation (21),
L A B = i = 1 N A j = 1 N B L p i j
with the additional assumption of neglecting the mutual inductance between the end-coils of both phases.
Comparing Equations (21) and (22), it is clear that L p h a s e A in Equation (21) can be obtained using Equation (22) as
L p h a s e A = L A A + L e n d _ p h a s e
so this work will focus in the effective computation of Equation (22).
The practical implementation of the proposed method seems to be cumbersome. First, a partial differential equation (PDE) must be solved for every conductor to obtain the MVP in Equation (18). Second, two double summations must be made in Equations (21) and (22), for all the stator and rotor phases. This process must be repeated for every position of the rotor, for obtaining the mutual inductance between a rotor phase and a stator phase. Nevertheless, in the next sections, an analytical expression for Equation (18) is derived, and a simple procedure is given for obtaining Equation (22), for any possible displacement between phases A and B, using a single FFT.

3.1. Partial Inductance between Axial Conductors

As only the partial inductances between the axial conductors of the phases are considered in Equations (21) and (22), only the axial component of the MVP must be taken into account. Neglecting end effects, a central cross section, perpendicular to the machine axis, is analyzed, thus using a 2D approximation for computing the MVP. Under these conditions, the MVP has a single, axial component, A z . If the length of the axial conductors is equal to the effective length of the magnetic core of the machine l m a q , the value of L p i j in Equation (18) is simply
L p i j = A z i j · l m a q I j
In an induction machine, the boundary condition at the infinity in Equation (18) can be replaced with a zero Dirichlet boundary condition at the external surface of the stator, where the MVP is assumed to be zero (no flux crossing the external surface of the stator). With this boundary condition, A z i j in Equation (18) is just the MVP at the conductor’s position, which is equal to the yoke flux per unit length at this position.
The method used in this work for obtaining the partial inductance L p i j in Equation (24) between two axial conductors i and j consists in feeding conductor j with a unit current, being this the only current in the machine, and computing the MVP at the position of conductor i, so that Equation (24) becomes
L p i j = l m a q · A z i j
For applying the proposed method, it is necessary to obtain the MVP generated by a single conductor, fed with a unit current. The MVP generated by a linear current in an air-gap with concentric circular cylindrical boundaries, as shown in Figure 3, must satisfy the Poisson’s equation
2 A z ( r , φ ) r 2 + 1 r A z ( r , φ ) r + 1 r 2 2 A z ( r , φ ) φ 2 = μ 0 J ( r , φ )
Equation (26) can be solved using the method of separation of variables, as in [70,71]. Placing the conductor at the origin of the angular coordinate ( φ = 0 in Figure 3), and assuming that the iron of the cylindrical boundaries is infinitely permeable, the solution of Equation (26) for the conductor shown in Figure 3 is
A z 0 ( r , φ ) = μ 0 J π × n = 1 1 n · c n c 2 n + b 2 n b 2 n a 2 n r 2 n + a 2 n r n cos ( n φ ) ln ( c ) r c n = 1 1 n · c n c 2 n + a 2 n b 2 n a 2 n r 2 n + b 2 n r n cos ( n φ ) ln ( r ) r > c
The constant, cyclic part of Equation (27) ( ln ( c ) , ln ( r ) ) can be discarded because, in any induction machine, the net sum of the linear currents through a transversal section is always zero. Figure 4 plots the harmonic part of Equation (27) for an example machine with a = 1 , c = 1.4 and b = 1.5 .
If the conductors are assumed to lie on the outer rotor surface or in the inner stator surface (Figure 5), then the linear currents are restricted to c = a or c = b in Figure 5. Besides, for obtaining the partial inductances of the conductors, the MVP must be calculated only on these same surfaces, that is, r = a or r = b in Figure 5. These four possibilities ( c = a , c = b , r = a , r = b ) reduce Equation (27) to only two different expressions: the MVP generated by an axial conductor on the same surface where it lies, or on the opposite one:
A z 0 ( φ ) = μ 0 J π × n = 1 1 n b 2 n + a 2 n b 2 n a 2 n cos ( n φ ) same surface n = 1 1 n 2 a n b n b 2 n a 2 n cos ( n φ ) opposite surface
The partial inductance of Equation (25) between two axial conductors i and j can now be computed using the values given by Equation (28). If conductor i is placed at the origin ( φ = 0 in Figure 3), and conductor j is placed at an arbitrary angular position φ , then the partial inductance between these two axial conductors is
L p 0 ( φ ) = μ 0 · l m a q π × n = 1 1 n b 2 n + a 2 n b 2 n a 2 n cos ( n φ ) same surface n = 1 1 n 2 a n b n b 2 n a 2 n cos ( n φ ) opposite surface
If now conductor i is assumed to be at an angular position φ = α , instead of φ = 0 , then the partial inductance with conductor j will be simply
L p α ( φ ) = L p 0 ( φ α )
The derivative of the partial inductance with respect to the angle, which is used for computing the torque in Equation (11), can be obtained directly from Equation (29) as
L p 0 ( φ ) φ = μ 0 · l m a q π × n = 1 2 a n b n b 2 n a 2 n sin ( n φ )
In Equation (31), only the derivative of the mutual inductance between conductors in opposite surfaces is needed, due to the uniformity of the air-gap.

3.2. Discretization of the Expression of the Partial Inductance between Axial Conductors

Equation (29) must be discretized in order to be used in an induction machine’s model. Therefore, the air-gap circumference is divided into N equal-length intervals, with an angular width equal to Δ φ = 2 π / N . In this case, the maximum number of harmonics of the infinite series in Equation (29) that are needed to represent the partial inductance is, by Shannon’s theorem, N / 2 . Thus, the sequence that gives the partial inductance between a conductor placed at the origin, and other conductor j placed at an angular position j · Δ φ , with j = 0 N 1 , is
L p 0 [ j ] = μ 0 · l m a q π × n = 1 n = N / 2 1 1 n b 2 n + a 2 n b 2 n a 2 n cos ( n · j Δ φ ) same surface n = 1 n = N / 2 1 1 n 2 a n b n b 2 n a 2 n cos ( n · j Δ φ ) opposite surface
Analogously to Equation (30), if conductor i is not placed at the origin, but at an angular position i · Δ φ , then Equation (32) becomes
L p i [ j ] = L p 0 [ ( ( j i ) ) N ]
where ( ( j i ) ) N = | j i | mod N is understood as a modulo N operation. In this work, to simplify the mathematical expressions, all the matrix indexes are considered as modulo N, without using the full notation of Equation (33). Following this convention, Equation (33) is expressed as
L p i [ j ] = L p 0 [ j i ]
The derivative with respect to the angular coordinate of the partial inductance between two conductors placed on opposite surfaces, given by Equation (31), can be also discretized with the same procedure, giving a sequence defined by
d L p 0 [ j ] = μ 0 · l m a q π × n = 1 n = N / 2 1 2 a n b n b 2 n a 2 n sin ( n · j Δ φ )

4. Computation of the Mutual Inductances between Two Phases Based on the Partial Inductance between Two Axial Conductors

4.1. Vectors Containing the Distribution of Conductors of Phases A and B and the Partial Inductances between Two Axial Conductors

The process for computing the mutual inductance between two phases A and B, such as those represented in Figure 5, is based on the substitution in Equation (22) of Equation (33). To perform this process, using discrete arithmetic, three column vectors of N elements are needed.
  • The vector L p 0 , whose element j, L p 0 [ j ] , contains the partial inductance between a conductor placed at the origin and other conductor placed at an angular position j · Δ φ , using the corresponding Equation (32) (depending if both conductors are in the same surface or in opposite surfaces).
    L p 0 = { L p 0 [ 0 ] , L p 0 [ 1 ] , , L p 0 [ N 1 ] }
  • The vector d L p 0 , whose element j, d L p 0 [ j ] , contains the derivative with respect to the angular coordinate of the partial inductance between a conductor placed at the origin and other conductor placed at an angular position j · Δ φ , using the corresponding Equation (35).
    L p 0 φ = d L p 0 = { d L p 0 [ 0 ] , d L p 0 [ 1 ] , , d L p 0 [ N 1 ] }
  • The vector Z A 0 , whose element j, n A 0 [ j ] , contains the number of conductors of phase A located in the angular interval of length Δ φ centered at position j · Δ φ . The direction of the currents at each position is represented using positive and negative values, depending on the current direction
    Z A 0 = { Z A 0 [ 0 ] , Z A 0 [ 1 ] , , Z A 0 [ N 1 ] }
    In the definition of Z A 0 , it is assumed that the axis of the distribution of conductors of the phase coincides with the origin of coordinates φ = 0 . If the phase is shifted by an angle φ A = k A · Δ φ , then all the elements of Z A 0 are shifted by k A , that is
    Z A k A [ j ] = Z A 0 [ j k A ]
  • The vector Z B 0 , which contains the distribution in the airgap of the conductors of phase B, which is defined in the same way as vector Z A 0
    Z B 0 = { Z B 0 [ 0 ] , Z B 0 [ 1 ] , , Z B 0 [ N 1 ] }
    In a similar way to the definition of Z A 0 , it is assumed that the axis of the distribution of conductors Z B 0 coincides with the origin of coordinates φ = 0 . If the phase is shifted by an angle φ B = k B · Δ φ , then all the elements of Z B 0 are shifted by k B , that is
    Z B k B [ j ] = Z B 0 [ j k B ]
It is worth mentioning that no restriction has been put on the definition of vectors Z A 0 and Z B 0 , so that arbitrarily complex winding layouts can be used in the proposed model, including the uniform distribution of the conductors along the slot width, or along skewed slots.

4.2. Vector L A B Containing the Mutual Inductances between Phases A and B for All Their Relative Positions

With the vectors defined in Section 4.1, the expression of the mutual inductance between phase A, with its axis located at an angular position φ A = k A · Δ φ , and phase B, placed at other angular position φ B = k B · Δ φ , L k A k B , can be formulated, using Equation (22), as
L k A k B = i = 0 N 1 j = 0 N 1 Z B k B [ i ] · L p i [ j ] · Z A k A [ j ]
Equation (42) can be formulated with vectors L p 0 , Z A 0 , Z B 0 , using Equations (34), (39), and (41) respectively, as
L k A k B = i = 0 N 1 j = 0 N 1 Z B 0 [ i k B ] · L p 0 [ j i ] · Z A 0 [ j k A ]
Using the change of variables j = j k A , Equation (43) becomes
L k A k B = i = 0 N 1 j = 0 N 1 Z B 0 [ i k B ] · L p 0 [ j ( i k A ) ] · Z A 0 [ j ]
and, using the change of variables i = i k A , Equation (44) becomes
L k A k B = i = 0 N 1 j = 0 N 1 Z B 0 [ i ( k B k A ) ] · L p 0 [ j i ] · Z A 0 [ j ]
That is, the mutual inductance between phases A and B depends only on their relative position ( k B k A ) · Δ φ . In this way, it can be defined a vector L A B whose element k, L A B [ k ] , contains the mutual inductance between phases A and B for an angular separation between them k · Δ φ
L A B [ k ] = i = 0 N 1 j = 0 N 1 Z B 0 [ i k ] · L p 0 [ j i ] · Z A 0 [ j ]
Equation (46) is simple to setup, because it contains only the expression of the partial inductance between two axial conductors (Equation (32)), and the distribution of the conductors (not the coils) of phases A and B, all of them obtained using a common origin of coordinates. However, at the same time, it is cumbersome to calculate, because a double summation is needed for obtaining the mutual inductance L A B [ k ] , for each value of k at N different positions. Nevertheless, this process can be performed in an extremely efficient way using the FFT, as proposed in the present work.

4.3. Matrix Formulation of the Expression of the Mutual Inductance between Phases A and B

Equation (45) can be expressed in matrix form as
L A B = TcZ B 0 T · TcL p 0 · Z A 0
where TcL p 0 is a square N × N matrix defined as
TcL p 0 = L p 0 [ 0 ] L p 0 [ N 1 ] L p 0 [ 1 ] L p 0 [ 1 ] L p 0 [ 0 ] L p 0 [ 2 ] L p 0 [ N 1 ] L p 0 [ N 2 ] L p 0 [ 0 ]
and TcZ B 0 T is the transpose of a square N × N matrix defined as
TcZ B 0 = Z B 0 [ 0 ] Z B 0 [ N 1 ] Z B 0 [ 1 ] Z B 0 [ 1 ] Z B 0 [ 0 ] Z B 0 [ 2 ] Z B 0 [ N 1 ] Z B 0 [ N 2 ] Z B 0 [ 0 ]
Matrix TcL p 0 is built using the column vector L p 0 in Equation (36) as its first column, and rotating each successive column by one position. That is, TcL p 0 is a Toeplitz circulant matrix. Matrix TcZ B 0 has been built using the same procedure, thus it is also a Toeplitz circulant matrix.

4.4. Computation of the Mutual Inductance between Phases A and B Using the FFT

A nice property of the Toeplitz circulant matrices is that they become diagonal by a transformation into the frequency domain. Therefore, Equation (48) can be reduced to a simple element-wise multiplication of three vectors in the frequency domain, using the FFT (and its inverse, the IFFT), as
L A B = IFFT FFT ( Z B 0 ) * FFT ( L p 0 ) FFT ( Z A 0 )
where the superscript * stands for the conjugate of each vector element, and the operator ∘ stands for the element-wise or Hadamard product of two vectors (written as “*” in MatLab environment). This equation is similar to the one presented in [72], where the convolution theorem has been applied to obtain the mutual inductances between two phases, based on the magnetomotive force generated by a single conductor, instead of the MVP.
Equation (50) gives, with a single element-wise multiplication of three vectors in the frequency domain, the values of the mutual inductance between phases A and B for every angular displacement between them (with a precision Δ φ ). It relies only in the distribution of the phases’ conductors, and on the analytical expression of the partial inductance between two axial conductors, Equation (32), which makes it very simple to implement. This approach differs from the approaches followed in [70,71], where the equation of the MVP in cylindrical coordinates (Equation (26)) must be solved again for each distribution of phase conductors, which limits its practical application to simple distributions, such as the sinusoidal ones.
The derivative with respect to the angular coordinate of the mutual inductance between phases A and B, located on opposite surfaces, has the same expression as Equation (50), just replacing the sequence L p 0 in Equation (32) by its angular derivative, the sequence d L p 0 in Equation (35), giving
L A B φ = d L A B = IFFT FFT ( Z B 0 ) * FFT ( d L p 0 ) FFT ( Z A 0 )
The procedure for obtaining the mutual inductance between two phases A and B, L A B in Equation (50), has been depicted graphically in the block diagram shown in Figure 6.
The same block diagram in Figure 6 is valid for computing the derivative of the mutual inductance between phases A and B, d L A B in Equation (51), just replacing in Figure 6 the vector L p 0 in Equation (36) by its derivative, the vector d L p 0 in Equation (37).

5. Experimental Validation

The proposed method has been applied to the fault diagnosis of a commercial squirrel-cage induction motor, whose characteristics are given in Appendix A. The type of fault that has been used for the experimental validation of the proposed model is a broken bars fault. The use of the proposed partial inductance approach is particularly well suited for modeling a squirrel cage induction motor with broken bars. In this case, instead of using rotor meshes formed by two consecutive bars, as usually done, the faulty rotor winding can be modeled using directly the bar currents in Equation (10), and the partial inductances of the rotor bars in Equations (7), (8) and (11). Following this approach, the presence of multiple broken bars, consecutive or not, can be introduced in the model just by eliminating the broken bars entries in the current and voltage vectors, Equations (9) and (10) , instead of modifying the inductances of the rotor meshes affected by the broken bars. This approach avoids the cumbersome process of recomputing the rotor inductances matrix to take into account the variations in some rotor meshes generated by the broken bars fault.
It is worth mentioning that the election of a commercial squirrel cage motor for performing the experimental validation of the proposed method does not exclude its validity for other types of induction machines, such as wound rotor induction machines. In fact, Equations (50) and (51) do not impose any constraint on the winding distribution, and are equally valid for both single rotor bars and wounded rotor phases.

5.1. Experimental Setup

The test equipment, displayed in Figure 7, consists of a current clamp, whose characteristics are given in Appendix B, a 200 pulse/revolution incremental encoder, a Yokogawa DL750 Oscilloscope and a Personal Computer connected to it via an intranet network. As shown in Figure 7, the broken bar fault has been artificially produced by drilling a hole in the selected bars.
The experimental validation of the proposed approach has been carried out using a set of artificially damaged rotors with one or two broken bars, beginning with two consecutive broken bars (Positions 1–2), and increasing gradually their separation (Positions 1–3, 1–4, etc.), up to a separation of seven slots between the two broken bars (Positions 1–8), as shown in Figure 8. An additional healthy rotor, with no broken bars, has also been used for comparison purposes. This gives a total number of nine experimental tests: a healthy motor, a motor with one broken bar, and seven motors with two broken bars at different positions.
The same stator has been used in all the experimental tests (Figure 9), which allows comparing the effects of the broken bar fault in the stator currents in a controlled environment. To this end, the motors were disassembled, and one of the stators was mounted in the test bed. All tests have been carried out mounting the various rotors (see Figure 9) in this stator. The induction motor under test (Appendix A) is connected via a belt coupling to a DC generator, which feds a resistive load (see Figure 10). Both the resistive load and the field excitation of the generator can be controlled, so that the induction machine works at rated speed 1410 r/min (s = 0.06). The current of a stator phase has been measured using a sampling frequency of 5000 Hz, during an acquisition time of 50 s, and its spectrum has been obtained for identifying the fault harmonics given by Equation (1).

5.2. Model Setup

The partial inductance between two conductors of the simulated machine in Equation (32) has been represented in Figure 11, for the cases of the two conductors lying in the same surface (Figure 11, top), or in opposite surfaces (Figure 11, middle). The derivative of the mutual inductance of a stator conductor and a rotor one in Equation (35) has been also represented in Figure 11, bottom. These values depend only on the geometrical dimensions of the machine.
The distribution of the conductors of a stator phase and the one that corresponds to a skewed rotor bar of the simulated machine are shown in Figure 12.
A direct substitution in Equations (50) and (51) of the values presented in Figure 11 and Figure 12 gives the mutual inductance between two phases, for every relative angular position, and its derivative respect to the angular coordinate. The mutual inductance between a stator phase and a rotor bar of the simulated machine, as a function of the rotor position, is given in Figure 13, top, and its angular derivative is given in Figure 13, bottom.
The simulations with the motor model have been carried using a constant load torque T L in Equation (12), with a value equal to the motor rated torque, that is, T L = 9550 × 1.1 / 1410 = 7.45 Nm.

5.3. Diagnosis of a Single Broken Bar Fault

The procedure for the diagnosis of a single broken bar fault is based on the MCSA technique. In the case of a broken bar fault, the magnitude of the characteristic fault harmonics given by Equation (1) increases significantly. For the tested and simulated motors, the rotor speed is the rated one (1410 r/min , f 1 = 50 Hz). The main fault harmonics used for the diagnosis are those corresponding to a value k = ± 1 in Equation (1): the lower sideband harmonic (LSH), with a frequency f L S H = ( 1 2 s ) f 1 , and the upper sideband harmonic (USH), with a frequency f U S H = ( 1 + 2 s ) f 1 . In the case of the tested and simulated motors, s = 0.06 , so that f L S H = ( 1 2 × 0.06 ) × 50 = 44 Hz, and f U S H = ( 1 + 2 × 0.06 ) × 50 = 56 Hz. Figure 14 shows the spectrum of the currents measured in the experimental tests, with a healthy and a faulty motor with a single broken bar. Figure 15 shows the same spectra obtained from the simulated motor.
The magnitude of the fault harmonics are presented in Table 1, showing a good agreement between the simulated and the experimental data. It is worth mentioning that, in the case of the experimental motor, two additional harmonics appear at frequencies of 43.5 Hz and 56.5 Hz, due to the belt used for coupling the load to the test bed. In fact, when the motor is tested unloaded, with the belt removed, these harmonic do not appear, so they are probably generated by an axial eccentricity induced by the asymmetric load coupling to the motor shaft. These harmonics do not appear in the case of the simulated motor, because the partial inductances model presented in this work does not take into account the effect of axial eccentricity.

5.4. Diagnosis of a Double Breakage Fault with Non-Consecutive Broken Bars

Major motor manufacturers have reported cases where the damaged bars appear randomly distributed around the rotor perimeter, indicating that the failure of non-adjacent bars is fairly common in large cage induction motors. Figure 16 shows the rotor of a 3.15 MW, 6 kV induction motor with a double breakage fault, affecting to non-consecutive rotor bars.
In the case of a double breakage, the magnitude of the LSH is a function of the relative position of the two broken bars. In [73] it has been shown that the ratio between the LSH in the case of double and single breakages depends on the angle between the broken bars as Equation
L S H p u = | L S H d o u b l e L S H s i n g l e | = | 2 cos ( p α b b ) |
where p is the number of pole pairs and α b b is the angle between the two broken bars. From Equation (52), it can be deduced that if α b b approximates π / 2 p , that is, half a pole pitch, then the second breakage reduces the magnitude of the LSH to a value lower than its magnitude in the case of a single breakage. Therefore, in this case, a motor with two broken bars could be erroneously diagnosed as a healthy motor. This behavior is more challenging to simulate than the single broken bar fault, and it has been selected to verify the validity of the proposed model, based on partial inductances, for fault diagnosis.
The following procedure has been followed to perform the experimental tests:
  • Seven rotors in healthy condition have been successively mounted and tested, to verify that there are no significant differences between the tested rotors.
  • A bar has been drilled in each of the rotors, denoted as b 1 bar. The motors have been tested at rated speed, and, using the spectrum of one of the stator phase currents, the magnitude of the LSH ( L S H s i n g l e in Equation (52)) has been recorded.
  • A second bar has been drilled at different positions from the first one (see Figure 8). The motors with a double breakage have been tested again at rated speed, and, using the spectrum of one of the stator phase currents, the magnitude of the LSH ( L S H d o u b l e in Equation (52)) has been recorded.
  • Finally, the ratio L S H p u has been computed for every rotor, and compared with the same value obtained with the simulated motor. The results are presented in Figure 17 and in Table 2. The time spent in each simulation (50 s in steady state) was 12 s, using a computer whose characteristic are given in Appendix C.
The results shown in Figure 17 and in Table 2 indicate that the results obtained with the proposed model clearly follow the experimental trend, which confirms its validity as a tool for fault diagnosis. It is worth mentioning that, with the proposed approach, all the experimental tests, using the rotors presented in Figure 8, have been simulated using a single inductances matrix, which contains the partial inductances between stator phases and rotor bars. Each type of the considered faults has been simulated just by eliminating the faulty bars from the set of Equations (7), (8) and (11), without any modification to the rest of the terms in these equations.

6. Conclusions

The development, improvement and implementation of fault diagnostics techniques for induction machines requires the use of fast and accurate dynamic models of the machine, which can take into account the asymmetries generated by the machine faults. In this paper, a novel approach has been proposed, using the concept of partial inductance. Instead of using the coil as the basic winding unit, the partial inductance between conductors has been proposed for building the matrices of phase inductances and their derivatives. The partial inductance of a single conductor has been obtained analytically using the magnetic vector potential, and the combination of the partial inductances of all phases has been solved using the FFT, which makes the proposed approach very fast to compute and very easy to implement. The proposed method has been theoretically presented and experimentally validated using the diagnosis of double breakage faults in the squirrel cage of a commercial induction motor for different, non-consecutive positions of the broken bars.
The application of this novel approach is not limited to the specific fault that has been used to perform the experimental validation, namely a broken bars fault. Other types of fault that can be simulated with the procedure presented in this paper include rotor asymmetries of wound-rotor inductions machines, stator inter-turn faults, bearing faults, oscillating loads and supply imbalances. Future work will extend this approach to induction machines with non-uniform air gap, which will allow the simulation of eccentricity faults (static, dynamic and mixed), and the effect of core saturation. In these cases, the analytical computation of phase inductances is a highly complex issue. The proposed method can alleviate this complexity by computing analytically only the MVP generated by a single conductor, which, inserted in the formulae presented in this work, Equations (50) and (51), can provide seamlessly the phases mutual inductances, regardless the complexity of the windings layout.

Author Contributions

This work was performed in collaboration among the authors: M.P.-S. directed the research; J.M.-R. contributed to the theory; M.R.-G. and A.S.-B. designed and validated the main methods and experiments; J.P.-C. and R.P.-P. analyzed the data.

Funding

This work was supported by the Spanish “Ministerio de Economía y Competitividad” in the framework of the “Programa Estatal de Investigación, Desarrollo e Innovación Orientada a los Retos de la Sociedad” (project reference DPI2014-60881-R).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Commercial IM

Three-phase induction machine. Rated characteristics: P = 1.1 kW , f = 50 Hz , U = 230 / 400 V , I = 2.7 / 4.6 A , n = 1410 r / min , cos φ = 0.8 .
Machine dimensions: Effective length of the magnetic core = 70.2 mm, radius at the middle of the air gap = 41.1 mm, air gap length = 1.2 mm.
Stator: Three-phase winding, 36 slots, 78 wires/slot, winding pitch = 7/9, slot opening width = 2.1 mm, phase resistance 7.68 Ω , end winding leakage = 2.3 mH.
Rotor: Squirrel-cage winding, 28 bars, slot opening width = 1.4 mm , skew = one slot pitch, bar resistance = 0.00202 m Ω , end winding leakage = 2.45 × 10 5 mH.

Appendix B. Current Clamp

Chauvin Arnoux MN60, Nominal measuring scope: 100 mA–20A, ratio input/output: 1 A/100 mV, intrinsic error: ≤ 2% + 50 mV, frequency use: 400 Hz–10 kHz.

Appendix C. Computer Features

CPU: Intel Core i7-2600K CPU @ 3.40 GHZ RAM memory: 16 GB, Matlab Version: 9.4.0.813654 (R2018a).

References

  1. Merizalde, Y.; Hernández-Callejo, L.; Duque-Perez, O. State of the Art and Trends in the Monitoring, Detection and Diagnosis of Failures in Electric Induction Motors. Energies 2017, 10, 1056. [Google Scholar] [CrossRef]
  2. Hossain, M.L.; Abu-Siada, A.; Muyeen, S.M. Methods for Advanced Wind Turbine Condition Monitoring and Early Diagnosis: A Literature Review. Energies 2018, 11, 1309. [Google Scholar] [CrossRef]
  3. Medina-García, J.; Sánchez-Rodríguez, T.; Galán, J.A.G.; Delgado, A.; Gómez-Bravo, F.; Jiménez, R. A Wireless Sensor System for Real-Time Monitoring and Fault Detection of Motor Arrays. Sensors 2017, 17, 469. [Google Scholar] [CrossRef] [PubMed]
  4. Pang, B.; Tang, G.; Tian, T.; Zhou, C. Rolling Bearing Fault Diagnosis Based on an Improved HTT Transform. Sensors 2018, 18. [Google Scholar]
  5. Liu, Y.; Bazzi, A.M. A review and comparison of fault detection and diagnosis methods for squirrel-cage induction motors: State of the art. ISA Trans. 2017, 70, 400–409. [Google Scholar] [CrossRef] [PubMed]
  6. Artigao, E.; Koukoura, S.; Honrubia-Escribano, A.; Carroll, J.; McDonald, A.; Gómez-Lázaro, E. Current Signature and Vibration Analyses to Diagnose an In-Service Wind Turbine Drive Train. Energies 2018, 11, 960. [Google Scholar] [CrossRef]
  7. Glowacz, A.; Glowacz, Z. Diagnosis of the three-phase induction motor using thermal imaging. Infrared Phys. Technol. 2017, 81, 7–16. [Google Scholar] [CrossRef]
  8. Glowacz, A.; Glowacz, W.; Glowacz, Z.; Kozik, J. Early fault diagnosis of bearing and stator faults of the single-phase induction motor using acoustic signals. Measurement 2018, 113, 1–9. [Google Scholar] [CrossRef]
  9. Glowacz, A.; Glowacz, Z. Diagnosis of stator faults of the single-phase induction motor using acoustic signals. Appl. Acoust. 2017, 117, 20–27. [Google Scholar] [CrossRef]
  10. Cablea, G.; Granjon, P.; Bèrenguer, C. Three-phase electrical signals analysis for mechanical faults monitoring in rotating machine systems. Mech. Syst. Signal Process. 2017, 92, 278–292. [Google Scholar] [CrossRef]
  11. Li, D.Z.; Wang, W.; Ismail, F. A Spectrum Synch Technique for Induction Motor Health Condition Monitoring. IEEE Trans. Energy Convers. 2015, 30, 1348–1355. [Google Scholar] [CrossRef]
  12. Samanta, A.K.; Naha, A.; Routray, A.; Deb, A.K. Fast and accurate spectral estimation for online detection of partial broken bar in induction motors. Mech. Syst. Signal Process. 2018, 98, 63–77. [Google Scholar] [CrossRef]
  13. Sapena-Bano, A.; Pineda-Sanchez, M.; Puche-Panadero, R.; Martinez-Roman, J.; Kanovic, Z. Low-Cost Diagnosis of Rotor Asymmetries in Induction Machines Working at a Very Low Slip Using the Reduced Envelope of the Stator Current. IEEE Trans. Energy Convers. 2015, 30, 1409–1419. [Google Scholar] [CrossRef]
  14. Mustafa, M.O.; Nikolakopoulos, G.; Gustafsson, T.; Kominiak, D. A fault detection scheme based on minimum identified uncertainty bounds violation for broken rotor bars in induction motors. Control Eng. Pract. 2016, 48, 63–77. [Google Scholar] [CrossRef]
  15. Gritli, Y.; Rossi, C.; Casadei, D.; Filippetti, F.; Capolino, G.A. A Diagnostic Space Vector-Based Index for Rotor Electrical Fault Detection in Wound-Rotor Induction Machines Under Speed Transient. IEEE Trans. Ind. Electron. 2017, 64, 3892–3902. [Google Scholar] [CrossRef]
  16. Ojaghi, M.; Yazdandoost, N. Oil-Whirl Fault Modeling, Simulation, and Detection in Sleeve Bearings of Squirrel Cage Induction Motors. IEEE Trans. Energy Convers. 2015, 30, 1537–1545. [Google Scholar] [CrossRef]
  17. Faiz, J.; Moosavi, S. Eccentricity fault detection - From induction machines to DFIG—A review. Renew. Sustain. Energy Rev. 2016, 55, 169–179. [Google Scholar] [CrossRef]
  18. Bayram, D.; Seker, S. Redundancy-Based Predictive Fault Detection on Electric Motors by Stationary Wavelet Transform. IEEE Trans. Ind. Appl. 2017, 53, 2997–3004. [Google Scholar] [CrossRef]
  19. Maraaba, L.; Al-Hamouz, Z.; Abido, M. An Efficient Stator Inter-Turn Fault Diagnosis Tool for Induction Motors. Energies 2018, 11. [Google Scholar] [CrossRef]
  20. Nguyen, V.; Seshadrinath, J.; Wang, D.; Nadarajan, S.; Vaiyapuri, V. Model-Based Diagnosis and RUL Estimation of Induction Machines Under Interturn Fault. IEEE Trans. Ind. Appl. 2017, 53, 2690–2701. [Google Scholar] [CrossRef]
  21. Barzegaran, M.; Mazloomzadeh, A.; Mohammed, O.A. Fault Diagnosis of the Asynchronous Machines Through Magnetic Signature Analysis Using Finite-Element Method and Neural Networks. IEEE Trans. Energy Convers. 2013, 28, 1064–1071. [Google Scholar] [CrossRef]
  22. Surya, G.N.; Khan, Z.J.; Ballal, M.S.; Suryawanshi, H.M. A Simplified Frequency-Domain Detection of Stator Turn Fault in Squirrel-Cage Induction Motors Using an Observer Coil Technique. IEEE Trans. Ind. Electron. 2017, 64, 1495–1506. [Google Scholar] [CrossRef]
  23. Jung, J.; Lee, S.B.; Lim, C.; Cho, C.H.; Kim, K. Electrical Monitoring of Mechanical Looseness for Induction Motors With Sleeve Bearings. IEEE Trans. Energy Convers. 2016, 31, 1377–1386. [Google Scholar] [CrossRef]
  24. Lu, S.; He, Q.; Yuan, T.; Kong, F. Online Fault Diagnosis of Motor Bearing via Stochastic-Resonance-Based Adaptive Filter in an Embedded System. IEEE Trans. Syst. Man Cybern. Syst. 2017, 47, 1111–1122. [Google Scholar] [CrossRef]
  25. Elbouchikhi, E.; Choqueuse, V.; Amirat, Y.; Benbouzid, M.E.H.; Turri, S. An Efficient Hilbert-Huang Transform-Based Bearing Faults Detection in Induction Machines. IEEE Trans. Energy Convers. 2017, 32, 401–413. [Google Scholar] [CrossRef]
  26. Elbouchikhi, E.; Choqueuse, V.; Benbouzid, M. Induction machine bearing faults detection based on a multi-dimensional MUSIC algorithm and maximum likelihood estimation. ISA Trans. 2016, 63, 413–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Ottewill, J.; Ruszczyk, A.; Broda, D. Monitoring tooth profile faults in epicyclic gearboxes using synchronously averaged motor currents: Mathematical modeling and experimental validation. Mech. Syst. Signal Process. 2017, 84, 78–99. [Google Scholar] [CrossRef]
  28. Burriel-Valencia, J.; Puche-Panadero, R.; Martinez-Roman, J.; Sapena-Bano, A.; Pineda-Sanchez, M. Short-Frequency Fourier Transform for Fault Diagnosis of Induction Machines Working in Transient Regime. IEEE Trans. Instrum. Meas. 2017, 66, 432–440. [Google Scholar] [CrossRef] [Green Version]
  29. Mustafa, M.O.; Varagnolo, D.; Nikolakopoulos, G.; Gustafsson, T. Detecting broken rotor bars in induction motors with model-based support vector classifiers. Control Eng. Pract. 2016, 52, 15–23. [Google Scholar] [CrossRef]
  30. Ojaghi, M.; Aghmasheh, R.; Sabouri, M. Model-based exact technique to identify type and degree of eccentricity faults in induction motors. IET Electr. Power Appl. 2016, 10, 706–713. [Google Scholar] [CrossRef]
  31. Mojallal, A.; Lotfifard, S. Multi-Physics Graphical Model based Fault Detection and Isolation in Wind Turbines. IEEE Trans. Smart Grid 2017, 1. [Google Scholar] [CrossRef]
  32. Singh, A.; Grant, B.; DeFour, R.; Sharma, C.; Bahadoorsingh, S. A review of induction motor fault modeling. Electr. Power Syst. Res. 2016, 133, 191–197. [Google Scholar] [CrossRef]
  33. Go, Y.; Song, M.H.; Kim, J.Y.; Choi, W.; Lee, B.; Kim, K.M. A New Algorithm for Analyzing Method of Electrical Faults of Three-Phase Induction Motors Using Duty Ratios of Half-Period Frequencies According to Phase Angle Changes. In Mechatronics and Robotics Engineering for Advanced and Intelligent Manufacturing; Zhang, D., Wei, B., Eds.; Springer: Cham, Switzerland, 2017; pp. 303–317. [Google Scholar]
  34. Strous, T.D.; Wang, X.; Polinder, H.; Ferreira, J.A. Evaluating Harmonic Distortions in Brushless Doubly Fed Induction Machines. IEEE Trans. Magn. 2017, 53, 1–10. [Google Scholar] [CrossRef]
  35. Berzoy, A.; Mohamed, A.A.S.; Mohammed, O. Impact of Inter-Turn Short-Circuit Location on Induction Machines Parameters Through FE Computations. IEEE Trans. Magn. 2017, 53, 1–4. [Google Scholar] [CrossRef]
  36. Tavana, N.R.; Dinavahi, V. Real-Time Nonlinear Magnetic Equivalent Circuit Model of Induction Machine on FPGA for Hardware-in-the-Loop Simulation. IEEE Trans. Energy Convers. 2016, 31, 520–530. [Google Scholar] [CrossRef]
  37. Guezmil, A.; Berriri, H.; Pusca, R.; Sakly, A.; Romary, R.; Mimouni, M.F. Detecting Inter-Turn Short-Circuit Fault in Induction Machine Using High-Order Sliding Mode Observer: Simulation and Experimental Verification. J. Control Autom. Electr. Syst. 2017, 28, 532–540. [Google Scholar] [CrossRef]
  38. Jandaghi, B.; Dinavahi, V. Hardware-in-the-Loop Emulation of Linear Induction Motor Drive for MagLev Application. IEEE Trans. Plasma Sci. 2016, 44, 679–686. [Google Scholar] [CrossRef]
  39. Sun, W.; Zhao, R.; Yan, R.; Shao, S.; Chen, X. Convolutional Discriminative Feature Learning for Induction Motor Fault Diagnosis. IEEE Trans. Ind. Inform. 2017, 13, 1350–1359. [Google Scholar] [CrossRef]
  40. Sun, W.; Shao, S.; Zhao, R.; Yan, R.; Zhang, X.; Chen, X. A sparse auto-encoder-based deep neural network approach for induction motor faults classification. Measurement 2016, 89, 171–178. [Google Scholar] [CrossRef]
  41. Pietrowski, W. Detection of time-varying inter-turn short-circuit in a squirrel cage induction machine by means of generalized regression neural network. COMPEL Int. J. Comput. Math. Electr. Electron. Eng. 2017, 36, 289–297. [Google Scholar] [CrossRef]
  42. Yin, Z.; Hou, J. Recent advances on SVM based fault diagnosis and process monitoring in complicated industrial processes. Neurocomputing 2016, 174, 643–650. [Google Scholar] [CrossRef]
  43. Chudasama, K.; Shah, V.; Shah, S. Induction Motor Relaying Scheme for External Faults Detection and Classification Using Subtractive Clustering Based Sugeno Fuzzy Inference System. Electr. Power Compon. Syst. 2016, 44, 1149–1162. [Google Scholar] [CrossRef]
  44. Mezani, S.; Hamiti, T.; Belguerras, L.; Lubin, T.; Gerada, C. Computation of Wound Rotor Induction Machines Based on Coupled Finite Elements and Circuit Equation Under a First Space Harmonic Approximation. IEEE Trans. Magn. 2016, 52, 1–4. [Google Scholar] [CrossRef]
  45. Rezaee-Alam, F.; Rezaeealam, B.; Faiz, J. Unbalanced Magnetic Force Analysis in Eccentric Surface Permanent-Magnet Motors Using an Improved Conformal Mapping Method. IEEE Trans. Energy Convers. 2017, 32, 146–154. [Google Scholar] [CrossRef]
  46. Ajily, E.; Ardebili, M.; Abbaszadeh, K. Magnet Defect and Rotor Eccentricity Modeling in Axial-Flux Permanent-Magnet Machines via 3-D Field Reconstruction Method. IEEE Trans. Energy Convers. 2016, 31, 486–495. [Google Scholar] [CrossRef]
  47. Ojaghi, M.; Daliri, S. Analytic Model for Performance Study and Computer-Aided Design of Single-Phase Shaded-Pole Induction Motors. IEEE Trans. Energy Convers. 2017, 32, 649–657. [Google Scholar] [CrossRef]
  48. Amrhein, M.; Krein, P.T. Induction Machine Modeling Approach Based on 3-D Magnetic Equivalent Circuit Framework. IEEE Trans. Energy Convers. 2010, 25, 339–347. [Google Scholar] [CrossRef]
  49. Fourati, A.; Bourdon, A.; Feki, N.; Rémond, D.; Chaari, F.; Haddar, M. Angular-based modeling of induction motors for monitoring. J. Sound Vib. 2017, 395, 371–392. [Google Scholar] [CrossRef]
  50. Kiselychnyk, O.; Bodson, M.; Wang, J. Comparison of Two Magnetic Saturation Models of Induction Machines and Experimental Validation. IEEE Trans. Ind. Electron. 2017, 64, 81–90. [Google Scholar] [CrossRef]
  51. Berzoy, A.; Mohamed, A.A.S.; Mohammed, O. Complex-Vector Model of Interturn Failure in Induction Machines for Fault Detection and Identification. IEEE Trans. Ind. Appl. 2017, 53, 2667–2678. [Google Scholar] [CrossRef]
  52. Hafner, M.; Popescu, M.; Boglietti, A.; Cavagnino, A. Analytic Modeling of Inverter-Fed Induction Machines—A Practical Approach for Matching Measurement and Simulation Data. IEEE Trans. Ind. Appl. 2016, 52, 4710–4718. [Google Scholar] [CrossRef]
  53. Oumaamar, M.E.K.; Maouche, Y.; Boucherma, M.; Khezzar, A. Static air-gap eccentricity fault diagnosis using rotor slot harmonics in line neutral voltage of three-phase squirrel cage induction motor. Mech. Syst. Signal Process. 2017, 84, 584–597. [Google Scholar] [CrossRef]
  54. Toliyat, H.A.; Lipo, T.A.; White, J.C. Analysis of a concentrated winding induction machine for adjustable speed drive applications. I. Motor analysis. IEEE Trans. Energy Convers. 1991, 6, 679–683. [Google Scholar] [CrossRef]
  55. Ojaghi, M.; Faiz, J. Extension to Multiple Coupled Circuit Modeling of Induction Machines to Include Variable Degrees of Saturation Effects. IEEE Trans. Magn. 2008, 44, 4053–4056. [Google Scholar] [CrossRef]
  56. El Malek, M.A.; Abdelsalam, A.K.; Hassan, O.E. Induction motor broken rotor bar fault location detection through envelope analysis of start-up current using Hilbert transform. Mech. Syst. Signal Process. 2017, 93, 332–350. [Google Scholar] [CrossRef]
  57. Cheng, M.; Han, P.; Hua, W. General Airgap Field Modulation Theory for Electrical Machines. IEEE Trans. Ind. Electron. 2017, 64, 6063–6074. [Google Scholar] [CrossRef]
  58. Joksimovic, G. Parameterized dynamic model of cage induction machine. In Proceedings of the 2016 XXII International Conference on Electrical Machines (ICEM), Lausanne, Switzerland, 4–7 September 2016; pp. 513–518. [Google Scholar]
  59. Paul, C.R. The Concept of Partial Inductance. In Inductance: Loop and Partial; Wiley-IEEE Press: Hoboken, NJ, USA, 2010; pp. 195–245. [Google Scholar]
  60. Holloway, C.L.; Kuester, E.F.; Ruehli, A.E.; Antonini, G. Partial and Internal Inductance: Two of Clayton R. Paul’s Many Passions. IEEE Trans. Electromagn. Compat. 2013, 55, 600–613. [Google Scholar] [CrossRef]
  61. Ni, C.; Zhao, Z.; Cui, X. Partial Inductance of Conductor Segments With Coulomb Gauge in Quasi-Static Field. IEEE Trans. Electromagn. Compat. 2017, 59, 1125–1132. [Google Scholar] [CrossRef]
  62. Romano, D.; Antonini, G.; Ruehli, A.E. Time-Domain Partial Element Equivalent Circuit Solver Including Non-Linear Magnetic Materials. IEEE Trans. Magn. 2016, 52, 1–11. [Google Scholar] [CrossRef]
  63. Nguyen, T.S.; Duc, T.L.; Tran, T.S.; Guichon, J.M.; Chadebec, O.; Meunier, G. Adaptive Multipoint Model Order Reduction Scheme for Large-Scale Inductive PEEC Circuits. IEEE Trans. Electromagn. Compat. 2017, 59, 1143–1151. [Google Scholar] [CrossRef]
  64. Bauernfeind, T.; Baumgartner, P.; Biro, O.; Magele, C.; Preis, K.; Torchio, R. PEEC-Based Multi-Objective Synthesis of Non-Uniformly Spaced Linear Antenna Arrays. IEEE Trans. Magn. 2017, 53, 1–4. [Google Scholar] [CrossRef]
  65. Martinho, L.B.; Siau, J.; Bannwarth, B.; Guichon, J.M.; Chadebec, O.; Meunier, G.; Silva, V.C. Generalized PEEC Analysis of Inductive Coupling Phenomena in a Transmission Line Right-of-Way. IEEE Trans. Magn. 2016, 52, 1–4. [Google Scholar] [CrossRef]
  66. Cao, Y.S.; Makharashvili, T.; Cho, J.; Bai, S.; Connor, S.; Archambeault, B.; Jiang, L.; Ruehli, A.E.; Fan, J.; Drewniak, J.L. Inductance Extraction for PCB Prelayout Power Integrity Using PMSR Method. IEEE Trans. Electromagn. Compat. 2017, 59, 1339–1346. [Google Scholar] [CrossRef]
  67. Naderi, P. Modified magnetic-equivalent-circuit approach for various faults studying in saturable double-cage-induction machines. IET Electr. Power Appl. 2017, 11, 1224–1234. [Google Scholar] [CrossRef]
  68. Bortolozzi, M.; Tessarolo, A.; Bruzzese, C. Analytical Computation of End-Coil Leakage Inductance of Round-Rotor Synchronous Machines Field Winding. IEEE Trans. Magn. 2016, 52, 1–10. [Google Scholar] [CrossRef] [Green Version]
  69. Ojaghi, M.; Nasiri, S. Modeling Eccentric Squirrel-Cage Induction Motors With Slotting Effect and Saturable Teeth Reluctances. IEEE Trans. Energy Convers. 2014, 29, 619–627. [Google Scholar] [CrossRef]
  70. Hague, B. The Principles of Electromagnetism Applied to Electrical Machines; Dover: New York, NY, USA, 1962. [Google Scholar]
  71. Hoeijmakers, M.J. An Induction Machine Model Based on Analytic Two-Dimensional Field Computations. In Proceedings of the Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM 2004), Anacapri, Italy, 16–18 June 2004; pp. 43–48. [Google Scholar]
  72. Sapena-Bano, A.; Martinez-Roman, J.; Puche-Panadero, R.; Pineda-Sanchez, M.; Perez-Cruz, J.; Riera-Guasp, M. Induction machine model with space harmonics for fault diagnosis based on the convolution theorem. Int. J. Electr. Power Energy Syst. 2018, 100, 463–481. [Google Scholar] [CrossRef]
  73. Riera-Guasp, M.; Cabanas, M.F.; Antonino-Daviu, J.A.; Pineda-Sanchez, M.; García, C.H.R. Influence of Nonconsecutive Bar Breakages in Motor Current Signature Analysis for the Diagnosis of Rotor Faults in Induction Motors. IEEE Trans. Energy Convers. 2010, 25, 80–89. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Stator windings and rotor bars of a three-phase squirrel cage induction motor.
Figure 1. Stator windings and rotor bars of a three-phase squirrel cage induction motor.
Sensors 18 02340 g001
Figure 2. Rectangular current loop representing a winding coil.
Figure 2. Rectangular current loop representing a winding coil.
Sensors 18 02340 g002
Figure 3. Linear current in an air-gap with concentric circular cylindrical boundaries.
Figure 3. Linear current in an air-gap with concentric circular cylindrical boundaries.
Sensors 18 02340 g003
Figure 4. (Left) MVP that produces in the air-gap an axial conductor located at φ = 0 and r = c when fed with a 1A current; and (right) induction lines in the air-gap.
Figure 4. (Left) MVP that produces in the air-gap an axial conductor located at φ = 0 and r = c when fed with a 1A current; and (right) induction lines in the air-gap.
Sensors 18 02340 g004
Figure 5. Location of the winding conductors in the air-gap on the outer surface of the rotor (phase A) and on the inner surface of the stator (phase B), for two phases separated an angular distance φ .
Figure 5. Location of the winding conductors in the air-gap on the outer surface of the rotor (phase A) and on the inner surface of the stator (phase B), for two phases separated an angular distance φ .
Sensors 18 02340 g005
Figure 6. Block diagram for the computation of the mutual inductance between phases A and B with (50). (a) Construction of the vectors with the layout of the conductors of phases A, Z A , and B, Z B , and the partial inductance between a conductor placed at the origin and other conductor at successive positions, L p 0 ; (b) Element-wise multiplication in the frequency domain of these three vectors; (c) Conversion into the spatial domain of the result of the multiplication, which gives the vector L A B . The element k of this vector, L A B [ k ] , is the mutual inductance between phases A and B for a relative angular position between them of k × 2 π / N .
Figure 6. Block diagram for the computation of the mutual inductance between phases A and B with (50). (a) Construction of the vectors with the layout of the conductors of phases A, Z A , and B, Z B , and the partial inductance between a conductor placed at the origin and other conductor at successive positions, L p 0 ; (b) Element-wise multiplication in the frequency domain of these three vectors; (c) Conversion into the spatial domain of the result of the multiplication, which gives the vector L A B . The element k of this vector, L A B [ k ] , is the mutual inductance between phases A and B for a relative angular position between them of k × 2 π / N .
Sensors 18 02340 g006
Figure 7. Experimental setup: measurement equipment (left); and rotor of the motor whose characteristics are given in Appendix A, with a broken bar (right).
Figure 7. Experimental setup: measurement equipment (left); and rotor of the motor whose characteristics are given in Appendix A, with a broken bar (right).
Sensors 18 02340 g007
Figure 8. Tested rotors with two broken bars in different relative positions.
Figure 8. Tested rotors with two broken bars in different relative positions.
Sensors 18 02340 g008
Figure 9. Tested rotors with two broken bars in different relative positions, with the stator used for performing all the experimental tests.
Figure 9. Tested rotors with two broken bars in different relative positions, with the stator used for performing all the experimental tests.
Sensors 18 02340 g009
Figure 10. Schema of the loading of the experimental machine (left) and experimental setup (right). The induction motor under test (Appendix A) is connected to a DC generator via a belt coupling. The DC machine feeds a resistive load. Both the resistive load and the field excitation can be controlled so that the induction machine works at rated speed.
Figure 10. Schema of the loading of the experimental machine (left) and experimental setup (right). The induction motor under test (Appendix A) is connected to a DC generator via a belt coupling. The DC machine feeds a resistive load. Both the resistive load and the field excitation can be controlled so that the induction machine works at rated speed.
Sensors 18 02340 g010
Figure 11. Partial inductance between two conductors of the simulated machine (32), if the conductors are in the same surface (top), or in opposite surfaces (middle). Derivative of the partial inductance with the angle (35), when the conductors are in opposite surfaces (bottom).
Figure 11. Partial inductance between two conductors of the simulated machine (32), if the conductors are in the same surface (top), or in opposite surfaces (middle). Derivative of the partial inductance with the angle (35), when the conductors are in opposite surfaces (bottom).
Sensors 18 02340 g011
Figure 12. Distribution of the conductors of a stator phase (top) and a rotor bar (bottom).
Figure 12. Distribution of the conductors of a stator phase (top) and a rotor bar (bottom).
Sensors 18 02340 g012
Figure 13. Mutual inductance between a stator phase and rotor bar of the simulated machine as a function of the rotor position (top), and its angular derivative (bottom).
Figure 13. Mutual inductance between a stator phase and rotor bar of the simulated machine as a function of the rotor position (top), and its angular derivative (bottom).
Sensors 18 02340 g013
Figure 14. Spectrum of the experimental motor in healthy conditions (top); and with a single broken bar (bottom). The position and magnitude of the fault harmonics (LSH and USH) are indicated with arrows.
Figure 14. Spectrum of the experimental motor in healthy conditions (top); and with a single broken bar (bottom). The position and magnitude of the fault harmonics (LSH and USH) are indicated with arrows.
Sensors 18 02340 g014
Figure 15. Spectrum of the simulated motor in healthy conditions (top); and with a single broken bar (bottom). The position and magnitude of the fault harmonics (LSH and USH) are indicated with arrows.
Figure 15. Spectrum of the simulated motor in healthy conditions (top); and with a single broken bar (bottom). The position and magnitude of the fault harmonics (LSH and USH) are indicated with arrows.
Sensors 18 02340 g015
Figure 16. Rotor of an induction machine of 3.15 MW, 6 kV, with two non-consecutive broken bars.
Figure 16. Rotor of an induction machine of 3.15 MW, 6 kV, with two non-consecutive broken bars.
Sensors 18 02340 g016
Figure 17. Comparison between the simulated and the experimental data obtained for the motor with one broken bar (first column), and with two broken bars separated a variable number of healthy bars.
Figure 17. Comparison between the simulated and the experimental data obtained for the motor with one broken bar (first column), and with two broken bars separated a variable number of healthy bars.
Sensors 18 02340 g017
Table 1. Magnitude of the main fault harmonics generated by a broken bar fault in the tested and simulated motors.
Table 1. Magnitude of the main fault harmonics generated by a broken bar fault in the tested and simulated motors.
Healthy MotorFaulty Motor
LSHUSHLSHUSH
Experimental−63.03 dB−69.43 dB−32.69 dB−52.39 dB
Simulated−160 dB−166.5 dB−37.93 dB−55.52 dB
Table 2. Comparison between the ratio L S H p u of Equation (52) computed with the experimental tests and with the simulated motors.
Table 2. Comparison between the ratio L S H p u of Equation (52) computed with the experimental tests and with the simulated motors.
L S H p u
ExperimentalSimulated
Single broken bar ( b 1 )11
Broken bars b 1 & b 2 1.521.472
Broken bars b 1 & b 3 1.0981.151
Broken bars b 1 & b 4 0.75270.5965
Broken bars b 1 & b 5 0.43580.4089
Broken bars b 1 & b 6 0.98270 1.07
Broken bars b 1 & b 7 1.4251.467
Broken bars b 1 & b 8 1.7371.604

Share and Cite

MDPI and ACS Style

Pineda-Sanchez, M.; Puche-Panadero, R.; Martinez-Roman, J.; Sapena-Bano, A.; Riera-Guasp, M.; Perez-Cruz, J. Partial Inductance Model of Induction Machines for Fault Diagnosis. Sensors 2018, 18, 2340. https://doi.org/10.3390/s18072340

AMA Style

Pineda-Sanchez M, Puche-Panadero R, Martinez-Roman J, Sapena-Bano A, Riera-Guasp M, Perez-Cruz J. Partial Inductance Model of Induction Machines for Fault Diagnosis. Sensors. 2018; 18(7):2340. https://doi.org/10.3390/s18072340

Chicago/Turabian Style

Pineda-Sanchez, Manuel, Ruben Puche-Panadero, Javier Martinez-Roman, Angel Sapena-Bano, Martin Riera-Guasp, and Juan Perez-Cruz. 2018. "Partial Inductance Model of Induction Machines for Fault Diagnosis" Sensors 18, no. 7: 2340. https://doi.org/10.3390/s18072340

APA Style

Pineda-Sanchez, M., Puche-Panadero, R., Martinez-Roman, J., Sapena-Bano, A., Riera-Guasp, M., & Perez-Cruz, J. (2018). Partial Inductance Model of Induction Machines for Fault Diagnosis. Sensors, 18(7), 2340. https://doi.org/10.3390/s18072340

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop