Procedure for Detection of Stator Inter-Turn Short Circuit in AC Machines Measuring the External Magnetic Field †

: This paper presents a non-invasive procedure to detect inter-turn short circuit faults in the stator windings of AC electrical machines. It proposes the use of the stray external magnetic ﬁeld measured in the vicinity of the machine to determine stator faults. The originality introduced by this procedure is the analysis method presented in the paper, which when compared to usual diagnosis methods, does not require any data on the healthy state of the machine. The procedure uses the magnetic unbalance created by the rotor poles and the load variation in faulty cases. The presented method can be applied to induction and synchronous machines used as a motor or generator. It is based on the variation of sensitive spectral lines obtained from the external magnetic ﬁeld when the load changes. Analytical relationships are developed in the paper to justify the proposed method and to explain the physical phenomenon. To illustrate these theoretical considerations, practical experiments are also presented.


Introduction
In recent decades the reliability and the operational safety of electrical machines become an essential issue, so many research studies are focused on their monitoring, which is a crucial phase to prevent severe unexpected failure [1].For that, the development of acquisition, analysis, and decision techniques are necessary to ensure the detection and diagnosis of electrical machines.The success of these techniques requires a good knowledge of the machine and its behavior in the presence of an internal fault.Among all the diagnosis methods used for rotating electrical machines, one can find: For detection of electrical and mechanical faults in electrical AC machines, different techniques have been tested as those based on spectral analysis of the machine's vibrations [5] or current signature analysis [6][7][8][9][10].However, the interpretation of results requires an expertise level even more advanced than the method hereto described, making it difficult to ensure real democratization of those techniques.Therefore, monitoring systems are applied only to systems that require high operational safety (for example, in power generation plants).A reliable diagnosis technique which can detect a failure and avoid total damage of motors or generators with a simple and non-invasive monitoring system is of great importance.For this reason, the technology in this field is still in permanent evolution to develop advanced methods [11][12][13][14][15][16].
In the 1970s, a new technique using the analysis of external magnetic field was developed by Penman [17].It is a non-invasive technique and easy to implement.The drawback of the latter is the modeling of the magnetic field which depends on the motor housings with an important shielding effect or on the stator yoke.The determination of the external magnetic field requires the modeling of the internal sources and the ferromagnetic influence and the machine conducting materials.The computation of such a problem can be made using finite element software.However, the accurate modeling requires a large computational effort [18,19], especially when 3D modeling is performed.Another approach consists of adapting analytical solutions existing for simple geometries [20] but these methods, based on simplified geometry and under particular hypotheses, can be hardly exploited for electrical machines.In [21], a method based on the definition of attenuation coefficients can be easily combined with an analytical model of the machine.
Fault detection methods using the external magnetic field analysis are based on the property that any fault changes the magnetic field in the near vicinity of the machine.Difficulties for modeling and in the interpretation of this variable lead to exploit only qualitative features of the spectrum, like the appearance of sensitive spectral lines [22].More usually, studies on the exploitation of the external magnetic field for fault detection are generally limited to model internal consequences of the fault such as: changes in the m.m.f.distribution [23], interaction with the slotting effect [24], magnetic of electrical unbalance [25].Other researchers prove that the axial field can provide additional information [26][27][28].In [29], it is shown that the analysis of the stray flux after supply disconnection can be useful for diagnosis.Advanced exploitation has been developed to provide deeper information: In [30], a sophisticated inverse problem is used for fault detection.In [31], it is shown that the external magnetic field can give information concerning the location of the fault.To improve the diagnostic process, the use of two flux sensors is proposed in [32].Furthermore, all the diagnosis methods usually require the knowledge of the healthy state of the machine regardless of the physical variable considered.The fault detection is then based on the comparison of the signature for a given state with that of the presumed healthy state by considering an indicator determined from a measurement that is known to be sensitive to a fault.On the other hand, the machine load can be a disturbing factor for diagnosis, because it induces several healthy states.A further difficulty lies in the fact that the healthy state is practically never known until the failure occurs because the user never records the healthy signature beforehand.
This paper proposes a new solution that exploits the information of the external magnetic field measured around the machine.Generally, the load is a disturbing factor for diagnosis methods.However, in the proposed solution, the load variation is used to improve the diagnosis with the advantage that it does not require the knowledge of the machine's healthy state.The detection of a stator fault is based on a comparison between no-load and load operating conditions.Initially, the analytical modeling of the stray flux analysis in the presence of the stator inter-turn short circuit is proposed.Then, the experimental validation for an induction machine (IM) working in a motor and generator case and a wound smooth rotor synchronous generator are presented.The reliability of the method has been tested in a self-excited induction generator (SEIG) which makes it possible to test an unbalanced load case with large frequency slip.The experimental safety measures that have to be taken to ensure the reliable diagnosis are also presented.

Analytical Approach for Healthy Machine
The analytical developments concern first a three-phase, p pole pair induction machine with N s stator slots and N r rotor bars per pole pair.In this machine with n s turns per pole pair, the stator coils are series connected and supplied by three-phase balanced sinusoidal voltage with frequency f (angular frequency ω).For the modeling, it is assumed that the air-gap flux density content b is constant during the load variation and that the considered harmonic components are those generated by i s q sinusoidal no load currents of I s 0 rms, where q is the phase number.In the proposed model, the ε s air-gap magnetomotive force (m.m.f) generated by the stator winding is multiplied by the λ per area unit air-gap permeance (which takes into account the slotting effect) to obtain b.

Air-Gap Permeance
To determine the expression of the air-gap permeance, a specific model is considered as presented by Figure 1.In this model defined in [33], the slot shape is considered rectangular and the field lines crossing the air-gap are radial.

Analytical Approach for Healthy Machine
The analytical developments concern first a three-phase, p pole pair induction machine with N s stator slots and N r rotor bars per pole pair.In this machine with n s turns per pole pair, the stator coils are series connected and supplied by three-phase balanced sinusoidal voltage with frequency f (angular frequency ω).For the modeling, it is assumed that the air-gap flux density content b is constant during the load variation and that the considered harmonic components are those generated by s q i sinusoidal no load currents of s 0 I rms, where q is the phase number.In the proposed model, the ε s air-gap magnetomotive force (m.m.f) generated by the stator winding is multiplied by the  per area unit air-gap permeance (which takes into account the slotting effect) to obtain b.

Air-Gap Permeance
To determine the expression of the air-gap permeance, a specific model is considered as presented by Figure 1.In this model defined in [33], the slot shape is considered rectangular and the field lines crossing the air-gap are radial.In this model, kskr is a permeance coefficient that depends on the slot geometry: As the field lines never join the bottom of the slots, practically, the air-gap can be modeled considering fictitious slots with a depth equal to the fifth of their opening.With that assumption, can be defined as: where θ represents the angular position of the rotor tooth 1 axis relatively to the stator referential d s tied to the phase 1 axis, ks and kr are negative, null, or positive integers, α s is the angular abscissa of considered point M in the air-gap related to d s , and kskr is a permeance coefficient that depends on the slot geometry.Considering the slip s which appears in IM, θ is defined as: In this model, Λ kskr is a permeance coefficient that depends on the slot geometry: where µ 0 = 4π10 −7 is the permeability of the vacuum approximately equal to that of air:r s d = l s d /(l s e + l s d ), p s = l s e /5, r r d = l r d /(l r e + l r d ).In (1) the component, A sr is defined as: where e s = e + p s , e r = e + p r , e M = e + p s + p r .As the field lines never join the bottom of the slots, practically, the air-gap can be modeled considering fictitious slots with a depth equal to the fifth of their opening.With that assumption, λ can be defined as: where θ represents the angular position of the rotor tooth 1 axis relatively to the stator referential d s tied to the phase 1 axis, ks and kr are negative, null, or positive integers, α s is the angular abscissa of considered point M in the air-gap related to d s , and Λ kskr is a permeance coefficient that depends on the slot geometry.Considering the slip s which appears in IM, θ is defined as: Energies 2021, 14, 1132 4 of 22

Healthy Machine m.m.f
In the healthy machine, the m.m.f ε calculated relatively to the stator d s can be expressed as: where A s h s is a function which depends on the winding coefficient corresponding to the rank h s defined by: h s = 6k + 1, where k varies between −∞ to +∞.

Air-Gap Flux Density
The calculus developments lead to define b = ε s λ in the reference frame related to d s as the following: After regrouping the components of the same frequency and same polarity, it comes: where b K,H is an elementary component of K frequency rank, and H pole pair number defined as: and: The frequency rank K only depends on the rotor rank permeance kr.kr = 0 leads to define the fundamental (K = 1), kr = ±1 leads to the first slotting harmonics.

Transverse Field
The stray external magnetic field is obtained from leakage flux created by the different elements of the machine.It is the result of its combination of transverse and axial components.The transverse field measured in the plane perpendicular to the machine axis is an image of the b attenuated by the stator yoke and the axial field located in a plane containing the axis of the machine and created by the overhang effects of the winding.The attenuation caused by the external machine frame will not be taken into account.According to the position of the sensor, the measured field does not come from the same source, and does not result from the same physical phenomenon.Therefore, the measured signal will be sensitive to different field components and it will depend on the position of the sensor around the machine periphery.So, in practice, it is possible to measure mainly the transverse field by choosing an adequate position of the sensor so that the effects of the axial field are minimal.This position corresponds roughly to the middle of the sheet package.
Figure 2 presents a simplified geometry of the machine with the main geometrical parameters which appear in the analytical model with a smooth air gap [34].In the analytical model used to determine the flux density, only the normal component of the traverse field is considered.It takes into consideration the attenuation coefficient CH related to the stator yoke that influences the magnitude of bK,H component.It depends on the magnetic permeability μr influenced by the stator lamination, H pole pair number, and the geometrical parameters inner s R int and outer s ext R radius and can be defined as follows [35]: The evolution of CH versus H pole pair number with  A physical justification of the evolution given in Figure 3 can be done by considering the line field distribution.Actually, flux density components of high polarity have a line field concentrated at the level of the air-gap whereas those of low polarity have a line field that spreads on the whole stator core.This coefficient influences the amplitude of the traverse flux field.In the analytical model used to determine the flux density, only the normal component of the traverse field is considered.It takes into consideration the attenuation coefficient C H related to the stator yoke that influences the magnitude of b K,H component.It depends on the magnetic permeability µr influenced by the stator lamination, H pole pair number, and the geometrical parameters inner R s int and outer R s ext radius and can be defined as follows [35]:

Measurement of the Transverse Field
The evolution of C H versus H pole pair number with R s int = 71.6 mm, R s ext = 121 mm and µr = 1000 is presented in Figure 3.One can notice that the value of the attenuation coefficient decreases with the increase of H.In the analytical model used to determine the flux density, only the normal component of the traverse field is considered.It takes into consideration the attenuation coefficient CH related to the stator yoke that influences the magnitude of bK,H component.It depends on the magnetic permeability μr influenced by the stator lamination, H pole pair number, and the geometrical parameters inner s R int and outer s ext R radius and can be defined as follows [35]: The evolution of CH versus H pole pair number with int 71.
and μr = 1000 is presented in Figure 3.One can notice that the value of the attenuation coefficient decreases with the increase of H.A physical justification of the evolution given in Figure 3 can be done by considering the line field distribution.Actually, flux density components of high polarity have a line field concentrated at the level of the air-gap whereas those of low polarity have a line field that spreads on the whole stator core.This coefficient influences the amplitude of the traverse flux field.

Measurement of the Transverse Field
For the proposed diagnosis method, the normal component of the traverse field is measured with flux sensors.The sensors are placed in the closed vicinity of the stator A physical justification of the evolution given in Figure 3 can be done by considering the line field distribution.Actually, flux density components of high polarity have a line field concentrated at the level of the air-gap whereas those of low polarity have a line field that spreads on the whole stator core.This coefficient influences the amplitude of the traverse flux field.

Measurement of the Transverse Field
For the proposed diagnosis method, the normal component of the traverse field is measured with flux sensors.The sensors are placed in the closed vicinity of the stator yoke, allowing to consider in the model only the C H coefficient.Supposing the sensor placed at radius x = R s ext , the normal traverse flux density b x is defined as: Let us introduce b x K the harmonic of K rank of b x at the given point M' in the closed vicinity of the stator (x = R s ext , α s = α s 0 ), corresponding to the center of the wound flux sensor.b x K can be defined by: b x K can be computed by introducing complex quantities: At given α s 0 , the resulting flux density harmonics at Kω angular frequency is composed of several elementary components of different polarity H.

Structure of the Faulty Machine
In order to determine the influence of the faulty turns in the change of the flux density, a model considering a three-phase stator winding was developed.In this model, it is supposed that "y" turns from the n' s turns of an elementary section belonging to the phase q are short-circuited and that y is small compared with pn s , the total number of turns per phase.Therefore, it can be assumed that the current remains unchanged and has the same values in each phase in the faulty case.This hypothesis can therefore characterize the short circuit thanks to a model that preserves the original structure of the machine.This model assumes that the stator winding in presence of the fault is equivalent to the healthy winding, associated to "y" independent turns in which the short-circuit current circulates.It will be assumed that these two circuits have independent running.The healthy part of the winding generates therefore the same flux density components without fault.
The model of a faulty winding is presented in Figure 4 where the whole phase winding is composed of an elementary healthy section and one with short circuit turns.For both structures, it is assumed that the magnetic reaction of the rotor is such that only the fundamental of the stator currents for a running at no load will be considered to characterize the air-gap flux density.This way, the resulting air-gap flux density b * is equal to the initial one, b, to which the flux density b sc generated by the "y" turns flowing through by the current i s qsc : b * = b + b sc is added.
Energies 2021, 14, x FOR PEER REVIEW 6 of 22 yoke, allowing to consider in the model only the CH coefficient.Supposing the sensor placed at radius x = s ext R , the normal traverse flux density b x is defined as: x K b ˆ can be computed by introducing complex quantities: At given s 0  , the resulting flux density harmonics at Kω angular frequency is com- posed of several elementary components of different polarity H.

Structure of the Faulty Machine
In order to determine the influence of the faulty turns in the change of the flux density, a model considering a three-phase stator winding was developed.In this model, it is supposed that "y" turns from the n' s turns of an elementary section belonging to the phase q are short-circuited and that y is small compared with pn s , the total number of turns per phase.Therefore, it can be assumed that the current remains unchanged and has the same values in each phase in the faulty case.This hypothesis can therefore characterize the short circuit thanks to a model that preserves the original structure of the machine.This model assumes that the stator winding in presence of the fault is equivalent to the healthy winding, associated to "y" independent turns in which the short-circuit current circulates.It will be assumed that these two circuits have independent running.The healthy part of the winding generates therefore the same flux density components without fault.
The model of a faulty winding is presented in Figure 4 where the whole phase winding is composed of an elementary healthy section and one with short circuit turns.For both structures, it is assumed that the magnetic reaction of the rotor is such that only the fundamental of the stator currents for a running at no load will be considered to characterize the air-gap flux density.This way, the resulting air-gap flux density b * is equal to the initial one, b, to which the flux density bsc generated by the "y" turns flowing through by the current s qsc i : The short circuit current is defined as follows: The short circuit current is defined as follows: ϕ sc is the phase lag between the short circuit current and the phase 1 current as shown in Figure 5.This phase actually depends on several parameters such as the impedance Energies 2021, 14, 1132 7 of 22 that limits the short circuit current, the short circuit winding, and the position of the fundamental air-gap flux density relative to the phase q current (depending on the load).
Energies 2021, 14, x FOR PEER REVIEW 7 of 22 sc is the phase lag between the short circuit current and the phase 1 current as shown in Figure 5.This phase actually depends on several parameters such as the impedance that limits the short circuit current, the short circuit winding, and the position of the fundamental air-gap flux density relative to the phase q current (depending on the load).

Faulty Turns m.m.f
The magnetomotive force s qsc  generated by the "y" short circuit turns, shifted of  from ds, is shown on Figure 6    and h is a non-null relative integer, which can take consequently all the values of h s .Here, h  is defined as: As , the calculus developments lead to define this quantity in the reference frame related to d s .A grouping of parameters with same polarity and same frequency provides: with:

Faulty Turns m.m.f
The magnetomotive force ε s qsc generated by the "y" short circuit turns, shifted of α s q from ds, is shown on Figure 6 in the case of a 4 poles machine.It also shows the m.m.f ε s qel generated by the healthy elementary winding.ε s qsc is an unidirectional m.m.f and can be decomposed in rotating fields that rotate in opposite directions.In a stator referential, ε s qsc can be written as follows: Energies 2021, 14, x FOR PEER REVIEW 7 of 22 sc is the phase lag between the short circuit current and the phase 1 current as shown in Figure 5.This phase actually depends on several parameters such as the impedance that limits the short circuit current, the short circuit winding, and the position of the fundamental air-gap flux density relative to the phase q current (depending on the load).
with: A s h is a function that can be determined from the Fourier series of ε s qsc and h is a non-null relative integer, which can take consequently all the values of h s .Here, φ h is defined as: φ h = hα s q + φ sc .As b sc = λε s qsc , the calculus developments lead to define this quantity in the reference frame related to d s .A grouping of parameters with same polarity and same frequency provides: with: where ks and kr are the permeance rank of the rotor and stator which vary from −∞ to + ∞, equivalent to ks and kr defined in (10).b x sc,K sc calculated at the point M' is the harmonic of K rank, of magnitude bx sc,K sc : By comparing the values of the frequency rank K sc given by ( 18), which appears in the case of short circuit with frequency value taken by K in healthy case defined by (10), we can notice that, there is no bring new frequency in the signal spectrum.Therefore, to detect the failure, the classic diagnosis method which analyzes the increase (or decrease) of already existing lines needs to know the amplitude value in the healthy case as a reference, thus limiting their practice application.Concerning the polarities H and H sc , one can observe that H sc can take all positive and negative integers whereas H is a multiple of p. H sc can especially be equal to ±1 corresponding to components that are weakly attenuated by the stator iron.In the following, the properties relating to the dissymmetry generated by such components will be exploited.

Searching for Sensitive Components
The harmonic at Kω angular frequency is composed of several components of various polarities.The ones that have the highest contribution can be found by searching the lowest values of h s , h, ks, kr, k's, k'r, leading to the smallest possible H (or H sc ) corresponding to the less attenuated components.
Let us consider an induction machine with p = 2, N s = 24, N r = 16, supplied by 50 Hz power source.Tables 1 and 2 give the magnitude of the components relative to the magnitude of the fundamental for the healthy and faulty machines respectively.The short circuit current rms value is three times that of the line current: I s sc = 3I s 0 .In Tables 1 and 2, it can be observed that the highest components in the external magnetic field are mainly issued from air-gap components of the lowest polarity.In healthy conditions, the first slotting harmonics are obtained for: In the following, one will consider the harmonics of rank K = K sc = 17 (f = 850 Hz) because it is generated by components of the lowest polarity in healthy and faulty conditions (H = −2, H sc = −1).These components are highlighted in grey in Tables 1 and 2 and they will be considered as the sensitive components for the inter-turn short circuit fault.The sensitivity has been shown in previous works [24].It can also be noticed that for the healthy machine, only one predominant component can be associated to one given harmonic whereas for the faulty machine, one given harmonic is generated by several components of different polarity and of similar magnitude.In this case, the components with the highest magnitudes has the polarity Hsc = +1 and Hsc = −1.Moreover, these components, weakly attenuated by the stator frame, have a magnitude similar to that of the component related to the healthy machine as those highlighted in yellow in Table 2. Therefore, the fault does lead to a significant change in the magnitude of the spectral line at 850 Hz in the external magnetic field.In order to improve the detection, the analysis will be focused on the load-induced variation of sensitive spectral lines measured at two points in the transverse external magnetic field.

Influence of the Load
Only the sensitive components defined in the previous section are now considered: b x K relative to the healthy machine and b x sc,K sc relative to the faulty machine.They both merge together to generate the resulting harmonic flux density b * x K .According to the hypothesis related to the rotor magnetic reaction, it will be assumed that bx K does not change when the load increases.The same hypothesis will be considered for bx sc,K sc .Actually, only the phase lag ϕ sc and consequently φ h change with the load.Consequently, for position 1 (Pos.1)defined for α s = 0, and position 2 (Pos.2) defined for α s = π, the resulting harmonic b * x K can be expressed as follows: Pos. 2 : Considering resulting harmonics, b * x K obtained from ( 22) and ( 23), their difference consists in the sign of the second term b x sc,K sc which changes in the faulty case as a consequence of the polarity H sc = 1 in the expression of cosines ((cos(γ + π) = −cos(γ)), contrary to polarity H = 2 of the healthy term when it does not change (cos(γ + 2π) = cos(γ)).This dissymmetry, generated by the fault in (22) and ( 23) is, from a physical point of view, the base of the proposed detection procedure.Figure 7a   b ˆ).The theoretical approach neglects the phase variation, which leads to a partial justification of the physical phenomenon.The different directions of variation between the magnitudes of the in faulty case allow one to propose the presented procedure:

•
In the healthy case: The air-gap flux density does not change when the machine is loaded compared to the no load case and the term x K sc sc b , is null.In this case, the amplitudes of the terms

•
In the faulty case: In loaded conditions, the magnitude of the component at Kω angular frequency measured in Pos. 1 and Pos. 2 will evolve in opposed directions and this particularity is proposed to be exploited in the analysis of the harmonic Kω as a fault indicator.Compared to another existing proposed methods, one does not require the knowledge of the healthy state to be applied.

Application to a Salient Pole Synchronous Machine
Different kinds of synchronous machines exist, they differ only by the structure of their rotor: wound salient, or wound smooth rotor (turbo alternator), surface mounted, or The different directions of variation between the magnitudes of the b * x K in faulty case allow one to propose the presented procedure:

•
In the healthy case: The air-gap flux density does not change when the machine is loaded compared to the no load case and the term b x sc,K sc is null.In this case, the amplitudes of the terms b * x1

•
In the faulty case: In loaded conditions, the magnitude of the component at Kω angular frequency measured in Pos. 1 and Pos. 2 will evolve in opposed directions and this particularity is proposed to be exploited in the analysis of the harmonic Kω as a fault indicator.Compared to another existing proposed methods, one does not require the knowledge of the healthy state to be applied.

Application to a Salient Pole Synchronous Machine
Different kinds of synchronous machines exist, they differ only by the structure of their rotor: wound salient, or wound smooth rotor (turbo alternator), surface mounted, or interior permanent magnet rotor.The modeling of the intact synchronous machine and the faulty one is deduced from the studied induction machine where the slotting effect created by the dampers is neglected.Since the stator has the same structure for both machines, in the analytical model, it is considered that the effect created by the saliency rotor and wound smooth rotor of the synchronous machine is similar to the rotor slots of the induction machine.Consequently, if the rotor is a wound smooth rotor, the components generated by a healthy and faulty synchronous machine are given by ( 10) and (18) with particularity s = 0.For a salient rotor, as the pole number is equal to the rotor saliencies of the machine, the frequency ranks and pole pair numbers given by ( 10) can be likewise considered in the case of a synchronous machine, with N r = 2 and s = 0: For a faulty machine, (18) becomes: Tables 3 and 4 give the magnitude of the components relative to the magnitude of the fundamental respectively for the healthy and for the faulty synchronous machines with 4 poles, and N s = 18.The short circuit current rms value is three-time that of the line current: I sc s = 3I s 0 .The results are similar to those of induction machine.It appears that each harmonic at given kr mainly originates from the components of lowest polarity.It also appears that the sensitive harmonics are: kr = 7, ks = 0, h s = −5 lead to K = 15, H = 18 (f = 750 Hz); kr = 7, ks = 0, h s = 1 lead to K = 15, H = 30 (f = 750 Hz); and kr = 8, ks = 0, h s = −5 lead to K = 17, H = 22 (f = 850 Hz).In faulty conditions, the same harmonics are obtained for: kr = 7, ks = 0, h s = −5 lead to K = 15, H = 23 (f = 750 Hz); kr = 7, ks = 0, h s = 1 lead to K = 15, H = 29 (f = 750 Hz), and kr = 8, ks = 0, h s = −5 lead to K = 17, H = 27 (f = 850 Hz).In the proposed method, which is the same as in the induction machine, the harmonics obtained for kr' = ±7 and ±8 can be analyzed.The choice between two harmonics must be made taking into account as practical measurements of harmonics amplitudes

Presentation of the IM Test Bench
The tests are performed using a three-phase squirrel-cage induction machine with 4 poles, 50 Hz, 11 kW, 380/660 V, 22.3/13A, 1450 rpm, 48 stator slots, and 32 rotor bars (N s = 24, N r = 16).This machine presented in Figure 8a is supplied by the grid and is specially modified to enable inter-turn short circuits.For measurement of the external magnetic field, the sensors are 180 • spatially shifted around the frame of the machine (Figure 8b).The output connections to the terminal box can short-circuit an elementary section housed in a slot, which corresponds to 12.5% of the full phase winding for the induction machine (Figure 8c).
The equipment above the machine makes it possible to simulate a fault by shortcircuiting the coils.The machine can operate under no-load or various loading conditions.For the induction machine, the short circuit current is limited to 15.3 A using an external rheostat.Figure 8b displays locations 1 and 2 of the flux sensors required for the method.The tests consist of measuring and analyzing the magnetic field outside the machine in order to validate the proposed diagnosis method.The measurements are performed using two identical manufactures wound flux sensors, placed against the machine.It has been checked that the wire current from the machine to the terminal block does not disturb the measurement.As announced in the previous section, the position of sensors must be on the one hand α s 1 = 0 (Pos. 1) and, on the other hand, as accurately as possible in α s 2 = π with the input (Pos.2) so that the measurements could be interpreted in the context of the proposed method.For both measurements, the sensor signal measured by the PULSE in the input Signal 1 and Signal 2 is placed in the middle of the machine to reduce the influence of end windings.Tests have been performed on the two induction machines for three operating conditions (motor, generator connected to the power system, and self-excited induction generator).In the proposed method, which is the same as in the induction machine, the harmonics obtained for kr' = ±7and ±8 can be analyzed.The choice between two harmonics must be made taking into account as practical measurements of harmonics amplitudes

Presentation of the IM Test Bench
The tests are performed using a three-phase squirrel-cage induction machine with 4 poles, 50 Hz, 11 kW, 380/660 V, 22.3/13A, 1450 rpm, 48 stator slots, and 32 rotor bars (N s = 24, N r = 16).This machine presented in Figure 8a is supplied by the grid and is specially modified to enable inter-turn short circuits.For measurement of the external magnetic field, the sensors are 180° spatially shifted around the frame of the machine (Figure 8b).The output connections to the terminal box can short-circuit an elementary section housed in a slot, which corresponds to 12.5% of the full phase winding for the induction machine (Figure 8c).The equipment above the machine makes it possible to simulate a fault by short-circuiting the coils.The machine can operate under no-load or various loading conditions.For the induction machine, the short circuit current is limited to 15.3 A using an external rheostat.Figure 8b displays locations 1 and 2 of the flux sensors required for the method.The tests consist of measuring and analyzing the magnetic field outside the machine in order to validate the proposed diagnosis method.The measurements are performed using two identical manufactures wound flux sensors, placed against the machine.It has been checked that the wire current from the machine to the terminal block does not disturb the measurement.As announced in the previous section, the position of sensors must be on the one hand 1 0 s   (Pos. 1) and, on the other hand, as accurately as possible in 2 s   with the input (Pos.2) so that the measurements could be interpreted in the context of the proposed method.For both measurements, the sensor signal measured by the PULSE in the input Signal 1 and Signal 2 is placed in the middle of the machine to reduce the influence of end windings.Tests have been performed on the two induction machines for three

Induction Machine in Motor Case
The time variation of the electromotive force (emf) and the spectrum of the signal delivered by the sensor in the case of the healthy induction machine, measured in the near vicinity of the machine are shown in Figure 9a where the magnitudes are presented in dB.The first tooth harmonics at 750 and 850 Hz can be clearly observed in Figure 9b.The study is focused on the line at 850 Hz, where the theoretical analysis shows a polarity lower than that of the line at 750 Hz.This is reflected by a greater magnitude of the line at 850 Hz, a property that is measured experimentally.operating conditions (motor, generator connected to the power system, and self-excited induction generator).

Induction Machine in Motor Case
The time variation of the electromotive force (emf) and the spectrum of the signal delivered by the sensor in the case of the healthy induction machine, measured in the near vicinity of the machine are shown in Figure 9a where the magnitudes are presented in dB.The first tooth harmonics at 750 and 850 Hz can be clearly observed in Figure 9b.The study is focused on the line at 850 Hz, where the theoretical analysis shows a polarity lower than that of the line at 750 Hz.This is reflected by a greater magnitude of the line at 850 Hz, a property that is measured experimentally.It should be kept in mind that in practice, the line studied can move with the slip, but will be identified as the line at 850 Hz.The results obtained for this machine operating as a motor in healthy and short-circuit fault conditions are presented, taking into account the variation of the spectral line at 850 Hz under the load influence.Figure 10a   It should be kept in mind that in practice, the line studied can move with the slip, but will be identified as the line at 850 Hz.The results obtained for this machine operating as a motor in healthy and short-circuit fault conditions are presented, taking into account the variation of the spectral line at 850 Hz under the load influence.Figure 10a gives the variation of the spectral line at 850 Hz on the signal delivered by the sensor in Pos.It should be kept in mind that in practice, the line studied can move with the slip, b will be identified as the line at 850 Hz.The results obtained for this machine operating a motor in healthy and short-circuit fault conditions are presented, taking into account t variation of the spectral line at 850 Hz under the load influence.Figure 10a  Figure 10 shows that a difference appears between the amplitudes of significant lin at 850 Hz for the healthy machine in positions 1 and 2. Actually, due to the real geome of the machine, the elements that lead to the attenuation are not exactly the same around the machine.However, as aforementioned in the analytical development, o finds the same direction of variation (positive) of this component for both positions wh the machine is loaded.Figure 10 shows that a difference appears between the amplitudes of significant lines at 850 Hz for the healthy machine in positions 1 and 2. Actually, due to the real geometry of the machine, the elements that lead to the attenuation are not exactly the same all around the machine.However, as aforementioned in the analytical development, one finds the same direction of variation (positive) of this component for both positions when the machine is loaded.
For the faulty machine (Figure 11), the presence of the shorted turns introduces twopole flux density components that will be combined with those presented for the healthy machine.Analyzing the results presented above, it can be seen that:

•
in healthy conditions, the lines at 850 Hz evolve in the same direction with the change of the load, • in faulty conditions, the lines at 850 Hz in position 1, compared with those in position 2 vary in opposite directions.

•
the value of the short circuit current Isc can influence the evolution of the lines

Induction Machine in Generator Case
A second test of the proposed method is carried out in the case where the machine is used as a generator connected to the grid.The results are shown in Figure 11.
For this kind of operation, the angle ϕ sc of short-circuit undergoes a phase change of approximately π compared to motor operating mode.As in the case of motor operation, the obtained results show that for a healthy machine, the magnitudes of lines at 850 Hz increase with the load (Figure 12a).In both measurement positions, they vary in the same direction.For the faulty machine with the short-circuit current Isc limited to 8.2A (Figure 12b) harmonics magnitudes change the variation, its decrease in position 1 for 3.3 kW power load and increase in position 2. Figure 12c,d gives a zoom in the evolution of flux harmonics at frequency 850 Hz measured in Pos. 1 and Pos. 2. The experimental results confirm that the proposed method can be applied also as an induction machine operating as a generator.
For the faulty machine (Figure 11), the presence of the shorted turns introduces twopole flux density components that will be combined with those presented for the healthy machine.Analyzing the results presented above, it can be seen that:

Induction Machine in Generator Case
A second test of the proposed method is carried out in the case where the machine is used as a generator connected to the grid.The results are shown in Figure 11.
For this kind of operation, the angle sc  of short-circuit undergoes a phase change of approximately  compared to motor operating mode.As in the case of motor operation, the obtained results show that for a healthy machine, the magnitudes of lines at 850 Hz increase with the load (Figure 12a).In both measurement positions, they vary in the same direction.For the faulty machine with the short-circuit current Isc limited to 8.2A (Figure 12b) harmonics magnitudes change the variation, its decrease in position 1 for 3.3 kW power load and increase in position 2. Figure 12c,d gives a zoom in the evolution of flux harmonics at frequency 850 Hz measured in Pos. 1 and Pos. 2. The experimental results confirm that the proposed method can be applied also as an induction machine operating as a generator.

Synchronous Machine
In order to test the efficiency of the proposed method, the following tests are performed on a synchronous machine operating as a generator.This machine is a 4 poles wound smooth rotor (cylindrical rotor) with concentric rotor windings presented in Figure 13

Synchronous Machine
In order to test the efficiency of the proposed method, the following tests are performed on a synchronous machine operating as a generator.This machine is a 4 poles wound smooth rotor (cylindrical rotor) with concentric rotor windings presented in Figure 13

Synchronous Machine
In order to test the efficiency of the proposed method, the following tests are performed on a synchronous machine operating as a generator.This machine is a 4 poles wound smooth rotor (cylindrical rotor) with concentric rotor windings presented in Figure 13   A picture of the test bench is shown in Figure 14 where a DC machine is used as a prime mover and external flux sensors are placed around the synchronous generator frame.A picture of the test bench is shown in Figure 14 where a DC machine is used as a prime mover and external flux sensors are placed around the synchronous generator frame.For diagnosis, we use data obtained from the signals delivered by two flux sensors which are 180° spatially shifted around the frame of the machine as shown in Figure 8b.Here, the load induces a variation of the spectral line delivered by both sensors at 750 Hz.For presented results, a single turn is short-circuited (0.26%) in phase A of the stator.Figure 15a gives the variations of the line in healthy conditions and Figure 15b, in faulty ones.Figure 15c,d shows the zoom of harmonics variation in the faulty case for signals measured in Pos. 1 and Pos. 2. One can notice in a healthy case an identical evolution of the spectral line with load variation.In the faulty case, it is measured an increase in position 2 and a wave variation in position 1.This asymmetry evolution is in agreement with ( 22), (23), and phasor diagram presented in Figure 7.For the experimental tests, the synchronous machine has been modified, permitting a small short-circuit fault.For diagnosis, we use data obtained from the signals delivered by two flux sensors which are 180 • spatially shifted around the frame of the machine as shown in Figure 8b.Here, the load induces a variation of the spectral line delivered by both sensors at 750 Hz.For presented results, a single turn is short-circuited (0.26%) in phase A of the stator.Figure 15a gives the variations of the line in healthy conditions and Figure 15b, in faulty ones.Figure 15c,d shows the zoom of harmonics variation in the faulty case for signals measured in Pos. 1 and Pos. 2. One can notice in a healthy case an identical evolution of the spectral line with load variation.In the faulty case, it is measured an increase in position 2 and a wave variation in position 1.This asymmetry evolution is in agreement with (22), (23), and phasor diagram presented in Figure 7.For the experimental tests, the synchronous machine has been modified, permitting a small short-circuit fault.
For presented results, a single turn is short-circuited (0.26%) in phase A of the stator.Figure 15a gives the variations of the line in healthy conditions and Figure 15b, in faulty ones.Figure 15c,d shows the zoom of harmonics variation in the faulty case for signals measured in Pos. 1 and Pos. 2. One can notice in a healthy case an identical evolution of the spectral line with load variation.In the faulty case, it is measured an increase in position 2 and a wave variation in position 1.This asymmetry evolution is in agreement with ( 22), (23), and phasor diagram presented in Figure 7.For the experimental tests, the synchronous machine has been modified, permitting a small short-circuit fault.

Self-Excited Induction Generator
To confirm the applicability of the proposed method, another test has been carried out with an induction machine operating as a self-excited induction generator in a balanced and unbalanced load as it is indicated in Figure 16.This machine is a three-phase squirrel cage induction generator with 4 poles, 50 Hz, 3 kW, 380/660 V, 7.

Self-Excited Induction Generator
To confirm the applicability of the proposed method, another test has been carried out with an induction machine operating as a self-excited induction generator in a balanced and unbalanced load as it is indicated in Figure 16.This machine is a three-phase squirrel cage induction generator with 4 poles, 50 Hz, 3 kW, 380/660 V, 7.This machine does not offer the possibility to create artificial faults, but tests with variable frequency and the unbalanced load was carried out in order to estimate the reliability of the method in case of an unbalanced load.Actually, the aim is to check that the method provides a healthy response in unbalanced operating conditions.For IM operating as SEIG, the frequency changes with load and reactive power.If an unbalanced load appears, the symmetry of the magnetic field around the machine will be disturbed.However, this dissymmetry should not be confused with that created by a stator short-circuit fault.
For the tests, a resistive load star connected is used and the unbalanced case is obtained by disconnecting phase A of the SEIG from the load.The reactive power is kept constant and the flux sensors are placed 180° spatially shifted around the machine frame.
In this case, the frequency changes with the slip which is negative in generator operating condition (s < 0). Figure 17a shows the harmonic evolution for the load variation in This machine does not offer the possibility to create artificial faults, but tests with variable frequency and the unbalanced load was carried out in order to estimate the reliability of the method in case of an unbalanced load.Actually, the aim is to check that the method provides a healthy response in unbalanced operating conditions.For IM operating as SEIG, the frequency changes with load and reactive power.If an unbalanced load appears, the symmetry of the magnetic field around the machine will be disturbed.However, this dissymmetry should not be confused with that created by a stator shortcircuit fault.
For the tests, a resistive load star connected is used and the unbalanced case is obtained by disconnecting phase A of the SEIG from the load.The reactive power is kept constant and the flux sensors are placed 180 • spatially shifted around the machine frame.
In this case, the frequency changes with the slip which is negative in generator operating condition (s < 0). Figure 17a shows the harmonic evolution for the load variation in balanced mode.Here, the line at 564 Hz (corresponding to s = −3.5%,kr = −1 and fundamental harmonic frequency 32 Hz) evolves in the same direction when the load increases.This machine does not offer the possibility to create artificial faults, but tests with variable frequency and the unbalanced load was carried out in order to estimate the reliability of the method in case of an unbalanced load.Actually, the aim is to check that the method provides a healthy response in unbalanced operating conditions.For IM operating as SEIG, the frequency changes with load and reactive power.If an unbalanced load appears, the symmetry of the magnetic field around the machine will be disturbed.However, this dissymmetry should not be confused with that created by a stator short-circuit fault.
For the tests, a resistive load star connected is used and the unbalanced case is obtained by disconnecting phase A of the SEIG from the load.The reactive power is kept constant and the flux sensors are placed 180° spatially shifted around the machine frame.
In this case, the frequency changes with the slip which is negative in generator operating condition (s < 0). Figure 17a shows the harmonic evolution for the load variation in balanced mode.Here, the line at 564 Hz (corresponding to s = −3.5%,kr = −1 and fundamental harmonic frequency 32 Hz) evolves in the same direction when the load increases.In practice, the studied line moves with the slip, but here it will be identified as the line at 564 Hz.The time variation of the electromotive force measured by both sensors and the spectrum of the signal delivered by sensor 1 in a balanced case are presented in Figure 18.We can notice identical emf curves measured by flux sensors in Pos. 1 and Pos. 2. As In practice, the studied line moves with the slip, but here it will be identified as the line at 564 Hz.The time variation of the electromotive force measured by both sensors and the spectrum of the signal delivered by sensor 1 in a balanced case are presented in Figure 18.We can notice identical emf curves measured by flux sensors in Pos. 1 and Pos. 2. As in 11 kW IM, the analyzed harmonics is chosen to take into account (10) and the highest measured value.Figure 17b shows the magnitude of the line at 617 Hz (36 Hz for the fundamental harmonic frequency) in the case of unbalanced resistive load connected at the SEIG output.In this case (which is not a faulty case), the lines show the same variation in the two measurement positions although the emf curves shown in Figure 19 are not identical.Figure 17b shows the magnitude of the line at 617 Hz (36 Hz for the fundamental harmonic frequency) in the case of unbalanced resistive load connected at the SEIG output.
In this case (which is not a faulty case), the lines show the same variation in the two measurement positions although the emf curves shown in Figure 19 are not identical.
(c) Figure 17b shows the magnitude of the line at 617 Hz (36 Hz for the fundamental harmonic frequency) in the case of unbalanced resistive load connected at the SEIG output.In this case (which is not a faulty case), the lines show the same variation in the two measurement positions although the emf curves shown in Figure 19  This test confirms that the proposed method is sensitive to field disturbance generated by the short-circuit current and is not sensitive to frequency variation or unbalanced load.This test confirms that the proposed method is sensitive to field disturbance generated by the short-circuit current and is not sensitive to frequency variation or unbalanced load.

Practical Precautions
The method is applicable for the diagnosis of all AC machines (motors and generators) uses in an industrial process, except 2 pole machines because it does not generate a dissymmetry in the phase winding.It uses the stray external magnetic field which is an image of the air-gap flux density attenuated by the stator magnetic circuit.By measuring the field, a difference can be observed which appears between the specific amplitudes (850 Hz, 750 Hz) of the magnetic spectrum in two specific positions around the machine.As the difference appears in the faulty machine but also in the healthy case, the load (which generally is a disturbing factor) is used here for discrimination.The method consists in comparing the signals delivered by each sensor during of the load variation.As shown in Figure 12a, the sensitive harmonic increases in the healthy case for Pos. 1 and Pos. 2 when the load increases.Whereas, in the faulty case, the increase of the load leads to a decrease followed by an increase (Figure 12c).This asymmetric amplitude variation between Pos. 1 and Pos. 2 is used for diagnosis and it does not require a high level of expertise to analyze the data.It concerns:

•
the kind of sensor: the sensor has to be small enough so that one can consider that the measured value of the flux density is the same as the flux density at the center of the sensor.In order to get an acceptable value of the output voltage, the small size has to be compensated by the number of turns.But as the increase of the number of turns decreases the resonance frequency of the sensor, this latter has to be checked with an impedance analyzer to get sure that the frequency resonance is far from the analyzed frequency.The used sensors are 3.2 cm diameter and 1200 turns, which leads to a resonance frequency at 93 KHz.

•
the location of the sensors: both sensors have to be placed symmetrically from the machine axis, approximately in the middle of the stator, between the end bells.Concerning the accuracy of the symmetry of the sensors, it increases when the pole pair number p of the machine is high.Actually, an error in the sensor alignment will be multiplied by p and it could not be possible to consider that the first terms of ( 22) and ( 23) are the same, leading to a disturbance in the diagnosis method.

•
the number of measurements: practical tests have shown that when the sensors are far from the axis of the damaged coil, the results give the impression that the machine is healthy because the magnitudes at sensitive spectral lines always vary in the same direction with the load variation.Theoretically, this can be taken into account considering all the elementary components of Tables 2 and 4 instead of the predominant components grey highlighted.To avoid this confusion, it is necessary to perform several measurements such as at least one is closer to the faulty winding than the healthy winding, among the p windings that constitute a whole phase.Consequently, the minimal number of measurements should be 2p, but to increase the probability of fault detection, we recommend 3p measurement.Of course, in a practical point of view, the locations of the sensors have to take into account the nonaccessible areas, e.g., the areas of the mounting base or the terminal box (Figure 20).
number p of the machine is high.Actually, an error in the sensor alignment will be multiplied by p and it could not be possible to consider that the first terms of ( 22) and ( 23) are the same, leading to a disturbance in the diagnosis method.• the number of measurements: practical tests have shown that when the sensors are far from the axis of the damaged coil, the results give the impression that the machine is healthy because the magnitudes at sensitive spectral lines always vary in the same direction with the load variation.Theoretically, this can be taken into account considering all the elementary components of Tables 2 and 4 instead of the predominant components grey highlighted.To avoid this confusion, it is necessary to perform several measurements such as at least one is closer to the faulty winding than the healthy winding, among the p windings that constitute a whole phase.Consequently, the minimal number of measurements should be 2p, but to increase the probability of fault detection, we recommend 3p measurement.Of course, in a practical point of view, the locations of the sensors have to take into account the non-accessible areas, e.g., the areas of the mounting base or the terminal box (Figure 20).Then, when the obtained results give at least one position with opposed magnitude variations with the load, then a fault can be suspected.AI methods can be also used to limit the number of measurements [36].
• case of inverter fed: it should be checked that the chopping frequency is higher than the sensitive one.

Conclusions
This paper presents the application and the reliability of a noninvasive procedure for the detection of inter-turn short-circuit faults in the stator winding of induction and synchronous machines.This procedure uses the stray external magnetic field measured in the vicinity of the AC machines by two sensors placed at a specific angular shift to extract information about the faulty state of the machines.The analysis is carried out considering the change of the transverse field obtained from the air-gap flux density which brings up the properties related to the characteristics of the flux density components being found in the transverse external field.The new contribution of the proposed diagnosis method is the analysis strategy which eliminates the main disadvantage presented by other diagnosis methods using the comparison with a healthy state assumed to be known before the fault appears.Moreover, the described procedure is non-invasive and inexpensive, characterized by the exploitation of the spatial dissymmetry of the field created by a short circuit current in the area where the fault appears.This dissymmetry modifies the rotor slotting harmonics and induces specific signatures in the external magnetic field spectrum which are analytically demonstrated, identified, and presented in this paper.In order to eliminate the constructive dissymmetry which exists practically, the presented diagnosis method proposes a procedure where the external field is measured in the vicinity of the electrical machine in the two cases: loaded and not-loaded machine.In order to test the reliability of the method, different Isc levels and an unbalanced load are tested.The paper shows that the proposed procedure can be applied to induction and synchronous machines with wound smooth rotor, considering specific spectral lines for each machine, with frequencies depending on the number of rotor slots.These specific frequencies facilitate the use of this procedure in the presence of other magnetic fields as is the case in an industrial environment.It is also shown that unbalanced operating conditions does disturb the procedure, especially in the healthy case.

•
methods of appreciation, including techniques that use artificial intelligence [2,3]; • methods for the identification and estimation of physical parameters of the machine [4]; • methods based on modeling of signals that analyze the time variation and the spectral content of different physical quantities.The work presented in this paper concerns that kind of analysis.

Figure 1 .
Figure 1.Geometrical parameters of the machine used in the analytical model.
is the permeability of the vacuum approximately equal to that of air:

Figure 1 .
Figure 1.Geometrical parameters of the machine used in the analytical model.

Figure 2 .
Figure 2. Geometrical parameters of machine used in the analytical model.
is presented in Figure3.One can notice that the value of the attenuation coefficient decreases with the increase of H.

Figure 3 .
Figure 3. Attenuation coefficient evolution CH versus pole pair number H.

ForFigure 2 .
Figure 2. Geometrical parameters of machine used in the analytical model.

Figure 2 .
Figure 2. Geometrical parameters of machine used in the analytical model.

Figure 3 .
Figure 3. Attenuation coefficient evolution CH versus pole pair number H.

Figure 3 .
Figure 3. Attenuation coefficient evolution C H versus pole pair number H.

)
Let us introduce x K b the harmonic of K rank of b x at the given point M' in the closed vicinity of the stator (x = , corresponding to the center of the wound flux sensor.x K b can be defined by:

Figure 4 .
Figure 4. Proposed model for a faulty winding.

Figure 4 .
Figure 4. Proposed model for a faulty winding.

Figure 5 .
Figure 5. Current diagram in faulty case.
in the case of a 4 poles machine.It also shows the m.m.f s qel  generated by the healthy elementary winding.s qsc  is an unidirectional m.m.f and can be decomposed in rotating fields that rotate in opposite directions.In a stator referential, s qsc  can be written as follows:

Figure 6 .
Figure 6.m.m.f generated by the faulty turns.
function that can be determined from the Fourier series of s qsc 

Figure 5 .
Figure 5. Current diagram in faulty case.

Figure 5 .
Figure 5. Current diagram in faulty case.

3. 2 .
Faulty Turns m.m.fThe magnetomotive force s qsc  generated by the "y" short circuit turns, shifted of s q  from ds, is shown on Figure 6 in the case of a 4 poles machine.It also shows the m.m.f s qel  generated by the healthy elementary winding.s qsc  is an unidirectional m.m.f and can be decomposed in rotating fields that rotate in opposite directions.In a stator referential, s qsc  can be written as follows:

Figure 6 .
Figure 6.m.m.f generated by the faulty turns.
a function that can be determined from the Fourier series of s qsc  and h is a non-null relative integer, which can take consequently all the values of h s .Here, developments lead to define this quantity in the reference frame related to d s .A grouping of parameters with same polarity and same frequency provides:

Figure 6 .
Figure 6.m.m.f generated by the faulty turns.
gives the associated time phasor diagram for ϕ x K = 0.It clearly shows that the magnitudes of the complex quantities b * x1 K and b * x2 Kare different.If the load changes, then a variation of ϕ cc leads to a variation of ϕ x sc,K sc .Figure7bgives the time phasor diagram after a variation of ϕ x sc,K sc .It can be observed that the resulting magnitude has increased in Pos. 1 ( b * x1 K ) and has decreased in Pos. 2 ( b * x2 K ).The theoretical approach neglects the phase variation, which leads to a partial justification of the physical phenomenon.

22 K..
are different.If the load changes, then a variation of cc leads to a variation of Figure7bgives the time phasor diagram after a variation of x K sc sc , It can be observed that the resulting magnitude has increased in Pos. 1 (
in Pos. 1 and Pos. 2 keep similar values or at least evolve in the same way if the load of the machine change.

K
and b * x2K measured in Pos. 1 and Pos. 2 keep similar values or at least evolve in the same way if the load of the machine change.

Figure 8 .
Figure 8. Experimental test bench with special induction machine modified to enable inter-turn short circuits (a) test bench, (b) flux sensors spatially shifted at 180°, (c) stator with a short circuit box and output connections.

Figure 8 .
Figure 8. Experimental test bench with special induction machine modified to enable inter-turn short circuits (a) test bench, (b) flux sensors spatially shifted at 180 • , (c) stator with a short circuit box and output connections.

Figure 9 .
Figure 9. Signal delivered by PULSE acquisition system; (a) emf delivered by flux coil sensor, (b) signal spectrum delivered for the healthy machine.
gives the variation of the spectral line at 850 Hz on the signal delivered by the sensor in Pos. 1 and Figure 10b by the sensors in Pos. 2.

Figure 9 .
Figure 9. Signal delivered by PULSE acquisition system; (a) emf delivered by flux coil sensor, (b) signal spectrum delivered for the healthy machine.

Figure 9 .
Figure 9. Signal delivered by PULSE acquisition system; (a) emf delivered by flux coil sensor, (b) signal spectrum delivered for the healthy machine.

Figure 10 .
Figure 10.Harmonic components at 850 Hz under both measurement positions for healthy motor and load variation (a) measurement in position 1, (b) measurement in position 2.

Figure 10 .
Figure 10.Harmonic components at 850 Hz under both measurement positions for healthy motor and load variation (a) measurement in position 1, (b) measurement in position 2.

•
in healthy conditions, the lines at 850 Hz evolve in the same direction with the change of the load, • in faulty conditions, the lines at 850 Hz in position 1, compared with those in position 2 vary in opposite directions.• the value of the short circuit current Isc can influence the evolution of the lines (a) (b)

Figure 11 .
Figure 11.Harmonic components at 850 Hz under both measurement positions for faulty motor and load variation (a) short-circuit current Isc limited at 8.2 A, (b) short-circuit current limited at 15.3 A, (c) load variation for a sensor in position 1 and Isc = 15.3A, (d) load variation for a sensor in position 2 and Isc = 15.3 A.

Figure 11 .Figure 12 .
Figure 11.Harmonic components at 850 Hz under both measurement positions for faulty motor and load variation (a) short-circuit current Isc limited at 8.2 A, (b) short-circuit current limited at 15.3 A, (c) load variation for a sensor in position 1 and Isc = 15.3A, (d) load variation for a sensor in position 2 and Isc = 15.3 A. Energies 2021, 14, x FOR PEER REVIEW 15 of 22 and the following characteristics: 10 kVA, 230/400 V, 25/15 A, 50 Hz, 54 stator slots, (N s = 27), N r = 16,

Figure 12 .
Figure 12.Harmonic components at 850 Hz under both measurement positions for generator operating at load variation in healthy and faulty case: (a) healthy case, (b) faulty case, (c) variation in faulty case for measurement in position 1, (d) variation in the faulty case for measurement in position 2.

Figure 12 .
Figure 12.Harmonic components at 850 Hz under both measurement positions for generator operating at load variation in healthy and faulty case: (a) healthy case, (b) faulty case, (c) variation in faulty case for measurement in position 1, (d) variation in the faulty case for measurement in position 2.

Energies 2021 , 22 Figure 14 .
Figure 14.Experimental test bench with special synchronous generator modified to enable interturn short circuits.

Figure 14 .
Figure 14.Experimental test bench with special synchronous generator modified to enable inter-turn short circuits.

Figure 15 .
Figure 15.Harmonic components at 750 Hz under both measurement positions for synchronous generator operating at variable load: (a) healthy case, (b) faulty case, (c) faulty case in position 1, (d) faulty case in position 2.

Figure 15 .
Figure 15.Harmonic components at 750 Hz under both measurement positions for synchronous generator operating at variable load: (a) healthy case, (b) faulty case, (c) faulty case in position 1, (d) faulty case in position 2.

Figure 16 .
Figure 16.Self-excited induction generator connected to a resistive balanced and unbalanced load.

Figure 16 .
Figure 16.Self-excited induction generator connected to a resistive balanced and unbalanced load.

Figure 16 .
Figure 16.Self-excited induction generator connected to a resistive balanced and unbalanced load.

Figure 17 .
Figure 17.Amplitude variation for SEIG harmonic components at load variation: (a) healthy machine with balanced resistive load, (b) healthy machine with unbalanced resistive load.

Figure 17 .
Figure 17.Amplitude variation for SEIG harmonic components at load variation: (a) healthy machine with balanced resistive load, (b) healthy machine with unbalanced resistive load.

Energies 2021 ,Figure 18 .
Figure 18.Signals measured by the flux sensors for balanced load: (a) emf measured in position 1, (b) emf measured in position 2, (c) signal spectrum delivered from flux sensor in position 1.

Figure 18 .
Figure 18.Signals measured by the flux sensors for balanced load: (a) emf measured in position 1, (b) emf measured in position 2, (c) signal spectrum delivered from flux sensor in position 1.

Figure 18 .
Figure 18.Signals measured by the flux sensors for balanced load: (a) emf measured in position 1, (b) emf measured in position 2, (c) signal spectrum delivered from flux sensor in position 1.

Figure 19 .
Figure17bshows the magnitude of the line at 617 Hz (36 Hz for the fundamental harmonic frequency) in the case of unbalanced resistive load connected at the SEIG output.In this case (which is not a faulty case), the lines show the same variation in the two measurement positions although the emf curves shown in Figure19are not identical.

Figure 19 .
Figure 19.Signals measured by the flux sensors for unbalanced load: (a) emf measured in position 1, (b) emf measured in position 2.

Figure 20 .
Figure 20.Illustration of the 3 measured positions of the sensors around the machine.Figure 20.Illustration of the 3 measured positions of the sensors around the machine.

Figure 20 .
Figure 20.Illustration of the 3 measured positions of the sensors around the machine.Figure 20.Illustration of the 3 measured positions of the sensors around the machine.

Table 1 .
Components generated by the healthy induction machine.

Table 2 .
Components generated by the faulty induction machine.

Table 3 .
Components generated by the healthy synchronous machine.

Table 4 .
Components generated by the faulty synchronous machine.