# Effects of the Slot Harmonics on the Stator Current in an Induction Motor with Bearing Fault

Authors: Xiangjin Song; Jingtao Hu; Hongyu Zhu; Jilong Zhang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Bearing faults are the most frequent faults of induction motors. The current spectrum analysis is an easy and popular method for the monitoring and detection of bearing faults. After a presentation of the existing fault models, this paper illustrates an analytical approach to evaluate the effects of the slot harmonics on the stator current in an induction motor with bearing fault. Simple and clear theoretical analysis results in some new current characteristic frequencies. Experimental tests with artificial bearing outer raceway fault are carried out by measuring stator current signals. The experimental results by spectral analysis of the stator current agree well with the theoretical inference.


---

## Body

## 1. Introduction

Induction motors (IMs) play a critical role in agricultural, industrial, and commercial sectors. Such popularity is increasing mainly because of their simple organizational structure, easy maintenance, premium power efficiency, and high reliability. Even if induction motors are reliable, unexpected failures may still occur under environmental stresses which could lead to production shutdown and extra maintenance cost. Thus, the online monitoring and fault diagnosis of induction motors have been extensively studied and many techniques have been published about this issue, as described in [1]. The IMs faults can be divided into two major categories: electrical defects (rotor faults, stator winding failures, etc.) and mechanical faults (bearing faults, eccentricities, and shaft misalignment). Among them, bearing-related faults account for over 40% (large machines) to 90% (small machines) of total faults [2]. Therefore, it is important to research on rolling-element bearing faults diagnosis.Depending on the type of the measured variables, methods that are used for diagnosing IM bearing faults can be categorized as acoustic emission, temperature measurements, vibration monitoring, motor current signal analysis (MCSA), and so forth. Among them, vibration signals have been widely applied in IM bearing faults detection [3–6]. However, this method is invasive and needs to install special sensor on the motor. Not only does this increase the investment and maintenance workload, but also the vibration sensors are sensitive to mounting position [7, 8]. Therefore, MCSA has received considerable attention due to its low cost and noninvasive nature in recent years [7, 9–14]. Moreover, stator current signals are often previously measured for control purposes in the drive system and can be analyzed remotely and online [9, 11].The relationship between bearing vibration characteristics and current spectrum effects was first proposed by Schoen et al. based on the generation of a rotating eccentricity, but the author did not give a detailed theoretical derivation [14]. Based on the studies mentioned above, Blödt et al. put forward a new mathematical model by analyzing the combinations of permeance wave and magneto-motive force (MMF) harmonics, taking into consideration the radial rotor movement and bearing faulted-related load torque variations [13]. The rotor was assumed to be perfectly centered when there is no contact between a ball and the defect, and the electromagnetic torque was supposed to be stable in Blödt et al. [13]. A theoretical analysis of the physical link between faults which was modeled as a torque disturbance and current components was obtained by accounting that the electromagnetic torque was subjected to minor variations [6]. A rectangular pulse model was proposed by Cheng et al. viewing the effect of bearing banded failures as a rectangular pulse [15]. A new current eigenfrequencies expression was introduced by Ma et al. considering the IM carried with bearing fault and inherent eccentricity [16]. The aforementioned models, however, neglected rotor and stator slotting effects which modify the valid airgap permeance. These techniques may lead to unreliable IM fault diagnosis.In view of the advantages of the aforementioned models, together with the effects of rotor and stator slotting, a more comprehensive mathematical model is proposed in this paper based on rotating wave approach where airgap flux density is determined by the product of permeance and magneto-motive force waves. At last, the theoretical analysis is verified by bearing outer raceway fault detection and the experimental results are well in agreement with the theoretical investigation.The rest of the paper is constructed as follows. Section2 gives a brief description of the structure and fault types of a rolling-element bearing. The proposed current characteristics analysis theory is described in Section 3. Section 4 reports the experiment setup and the obtained results for testing the effectiveness of the proposed theory on a three-phase induction motor with one healthy and one artificially damaged bearing. Finally, some conclusions of this paper are summarized in Section 5.


## 2. The Structure and Fault Types of a Rolling-Element Bearing

The typical geometric structure of a rolling-element bearing is shown in Figure1. Rolling-element bearings are mainly composed of the outer and inner raceway, the ball elements, and the cage, which ensures the distances between the ball elements uniformly. On ordinary conditions, bearing faults can be separated into two types, that is, single-point (also called localized or cyclic) and generalized-roughness (also called distributed or noncyclic) faults [17]. Localized faults could be imagined as a small hole, a pit, or a missing piece of material and they affect a localized region [13, 17]. When a localized defective bearing runs at a constant speed, a periodic impact between rolling elements and the raceway is produced [6, 13]. The existence of such impulsive forces gives rise to an increase in the vibration level [18, 19]. The characteristic frequencies of these vibrations rely on fault location, bearing geometry, and operating speed [20, 21]. Not alike localized faults, generalized-roughness faults greatly worsen the entire area of a bearing element because of the absence of lubrication, erosion, or bearing pollution and are hard to represent by specific frequencies [22, 23]. In this paper, only localized fault is studied.Figure 1
Typical structure of a rolling-element bearing with main parameters.Localized faults can be usually divided into four categories depending on the affected element, namely, outer raceway fault, inner raceway fault, ball fault, and cage fault. The vibration characteristic frequencies of these faults, when the outer raceway is stationary, are expressed as follows [20]:(1)fof=Nb2fr1-DbDccos⁡β,fif=Nb2fr1+DbDccos⁡β,fbf=DcDbfr1-DbDccos⁡β2,fcf=12fr1-DbDccos⁡β,where fof is the outer raceway fault frequency, fif is the inner raceway fault frequency, fbf is the ball fault frequency, fcf is the cage fault frequency, fr is the shaft rotational frequency, Nb is the number of balls, Db is the ball diameter, Dc is the ball pitch diameter, and β is the contact angle between a ball and the raceway.If the bearings are with a number of rolling elements between 6 and 12, the outer and inner raceway fault frequencies can be computed as follows [12]:(2)fof=0.4Nbfr,fif=0.6Nbfr.In rolling-element bearing, the stationary ring (outer raceway in this case) fault happens first due to the stationary ring material in load area receiving more dynamic load cycles than other bearing components (e.g., the rotating ring and rolling elements) [7]. Hence, the outer raceway faults are within the scope of investigations in this paper.


## 3. Current Characteristics Analysis Theory

In a typical induction motor, the dominant components of stator current are the fundamental supply frequency, the eccentricity harmonics, the slot harmonics, and other components including environmental noise [24, 25]. The models for bearing fault detection proposed by Blödt et al. are limited that the rotor is perfectly centered when there is no contact between a ball and the defect [13]. The new current eigenfrequencies expression considering the IM carried with bearing fault and inherent eccentricity was introduced by Ma et al. [16]. However, the other prime components of stator current, namely, the slot harmonics, are ignored. In this section, a more comprehensive mathematical model is proposed, accounting for stator and rotor slot harmonics, airgap eccentricity harmonics, and bearing outer raceway fault harmonics. Some new current characteristic frequencies are discovered around the rotor slot harmonics and can be used for bearing fault detection.The airgap permeance, accounting for stator and rotor slotting, can be approximated, expressed as follows [26]:(3)Pθ,t=δsθ,tδrθ,tgθ,t,where g(θ,t) is the radial airgap length taking into account bearing fault and inherent eccentricity and(4)δsθ,t=∑ns=0+∞δnscos⁡nsSθis the permeance of airgap with restriction of a slotted stator and a smooth rotor and(5)δrθ,t=∑nr=0+∞δnrcosnrRθ-nrRwrtis the permeance of airgap with restriction of a slotted rotor and a smooth stator.S and R are the number of stator and rotor slots, wr is the rotor mechanical rotational speed, θ is the angular position referring to stator, δns and δnr are the amplitudes of permeance related to the stator and rotor slots effects, ns and nr are any integer, and t is the time. Hereafter, n represents any integer.In the following text, after the explicit mathematical derivation process of the airgap permeance and the airgap flux density, some new characteristic frequencies which can be used for bearing outer raceway fault detection around the rotor slot harmonics are obtained.Under ideal condition, the airgap length in a machine is constant. However, due to unsuitable assembling and manufacturing process, an inherent level of static and dynamic eccentricity tends to coexist. In the case of static eccentricity, by assuming a smooth stator and rotor surface, the length of the airgap can be expressed as follows [27]: (6)gseθ=g01-δsecos⁡θ,where g0 is the airgap length in an induction motor under ideal condition and δse is the static eccentricity degree.In the case of dynamic eccentricity, the length of the airgap can be expressed as follows [27]:(7)gdeθ,t=g01-δdecos⁡θ-wrt,where δde is the dynamic eccentricity degree.If Fourier series are developed, then,(8)11-δsecos⁡θ=∑nse=0+∞δnsecos⁡nseθ,11-δdecos⁡θ-wrt=∑nde=0+∞δndecos⁡ndeθ-ndewrt.The airgap length of an ideal motor with bearing outer raceway fault can be described as follows [13]:(9)gofθ,t=g01-ecos⁡θ∑nof=-∞+∞δt-noffof,where e is the eccentric degree because of the outer raceway fault and δ(t-nof/fof) is Dirac function.When the Dirac function is expressed as a complex Fourier series, (9) can be represented as follows:(10)gofθ,t=g01-∑nof=0+∞ecnofcos⁡θ±nofwoft,where cnof=fof∀nof and wof=2πfof.In order to simplify the airgap permeance coefficient, the fraction1/(1-x)=1+x+x2+x3+⋯≈1+x (for x<1) is used. The condition x<1 is constantly met because the degree of eccentricity confirmed e<1 to prevent any contact between the rotating rotor and nonrotating stator. Then the airgap permeance with bearing outer raceway fault can be obtained(11)δofθ,t=1gofθ,t=δ0+∑nof=0+∞δnofcos⁡θ±nofwoft,where δ0=1/g0 and δnof=eδ0cnof∀nof.The airgap length for the induction motor with bearing outer raceway fault and mixed eccentricity can thus be represented as follows [16]:(12)gθ,t=g01-δsecos⁡θ1-δdecos⁡θ-wrt1-∑nof=0+∞ecnofcos⁡θ±nofwoft.A good approximation for the airgap permeance which takes into account the mixed eccentric motor with bearing outer raceway fault is given by(13)δθ,t=1gθ,t=δ0∑nse=0+∞∑nde=0+∞δnse,ndecos⁡nse±ndeθ-ndewrt+∑nse=0+∞∑nde=0+∞∑nof=0+∞δnse,nde,nofcos⁡nse±nde±1θ-ndewr±nofwoft.Then taking into account the interaction of stator and rotor slot harmonics, airgap eccentricity harmonics, and bearing outer raceway fault harmonics, the permeance of the airgap could be written as follows:(14)Pθ,t=δsθ,tδrθ,tδθ,t=δ0∑ns=0+∞∑nr=0+∞∑nse=0+∞∑nde=0+∞δns,nr,nse,ndecos⁡nsS±nrR±nse±ndeθ-nrR±ndewrt+∑ns=0+∞∑nr=0+∞∑nse=0+∞∑nde=0+∞∑nof=0+∞δns,nr,nse,nde,nofcos⁡nsS±nrR±nse±nde±1θ-nrR±ndewr±nofwoft.Only the stator magneto-motive force (MMF) is taken into account and, neglecting skew, it can be represented by(15)Ftotθ,t=∑nθf=1+∞∑nwf=-∞+∞Fnθf,nwfcos⁡nθfpθ-nwfw1t,where F is the amplitude of stator MMF and w1 is the electrical supply frequency.Then the airgap flux density, deduced from the product of the permeance (14) and MMF (15), can be given by(16)Bθ,t=∑me,ΩeBme,Ωecos⁡meθ-Ωet+∑mf,ΩfBmf,Ωfcos⁡mfθ-Ωft,where(17)me=nsS±nrR±nse±nde±nθfp,Ωe=nrR±ndewr±nwfw1,mf=nsS±nrR±nse±nde±nθfp±1,Ωf=nrR±ndewr±nofwof±nwfw1.Then, the change in the flux density distribution results in the current having harmonics as follows:(18)fe=nrR±ndefr±nwff1,(19)ff=nrR±ndefr±noffof±nwff1.The models for bearing fault detection proposed by Blödt et al. show that bearing outer raceway characteristic frequenciesf1±noffof can be obtained in the stator current spectrum. The current eigenfrequencies expression f1±ndefr±noffof was introduced by Ma et al. considering the IM carried with bearing fault and inherent eccentricity. From (19), we can see that the new bearing outer raceway characteristic frequencies |f1±noffof±(nrR±nde)fr| (with nr≠0) can also arise around the rotor slot harmonics when considering not only the slot harmonics but also mixed eccentricity in an induction motor with bearing fault. These new bearing outer raceway characteristic frequencies have not been found in the aforementioned papers. It should be possible to detect these characteristic frequencies by analyzing the stator current spectrum.


## 4. Experimental Results

To assess the effectiveness of the proposed techniques for IM bearing fault detection, a series of comparisons have been conducted between the spectrum of the stator current in healthy and faulty cases. In a rolling-element bearing, artificially produced defect is introduced by simply drilling a hole in the outer raceway to simulate a localized fault, which is similar to other literatures [9–12].This section is divided into two parts: first, an overview of the experiment test rig is presented, and then the proposed techniques are verified offline using stator current signals for bearing outer race fault detection.


### 4.1. Test Rig

Schematic diagram of the experiment setup used in this paper is shown in Figure2. The rated data of the used induction motor are given in Table 1. The induction motor is mechanically coupled to a DC machine, which feeds an external resistor bank. The motor can be identically loaded at different speed from zero to rated load by adjusting the field voltage of the DC machine. Both machines are fixed on a steel plate test bench resting on a cement floor. A photoelectric speed sensor (SZGB-6) is used to measure the shaft speed. As the used National Instrument Data Acquisition (NI DAQ 6216) card accepts only voltage signal with maximum amplitudes of ±10 V, JLB-21 Hall current sensor is used to gather the current signal first, and then the acquired current signal is converted to proportional voltage signals by conditioning circuit. A PC equipped with DAQ card through a USB cable is used to store them in the memory. A picture of the experimental setup is shown in Figure 3.Table 1
Test induction motor rated parameters.
Nominal power
3 kW
Nominal rated voltage
380 V
Nominal rated current
6.8 A
Frequency
50 Hz
Pole pairs
2
Stator slots
36
Rotor slots
32
Stator winding connection
YFigure 2
Schematic diagram of the experiment setup.Figure 3
Photograph of the experimental setup used for bearing fault detection.The induction machine has two CU6206RZ type bearings (single row and deep groove ball bearings) with a width ofw = 16 mm, an outside diameter of D = 62 mm, an inside diameter of d = 30 mm, and a pitch diameter of Dc = 46.0 mm. Each bearing has nine balls (Nb = 9) with an approximate diameter of Db = 9.52 mm. The contact angle of the ball with the race (β) is assumed to be 0°. Two identical motors are prepared: one motor carries healthy bearings and the other carries one defected bearing which are mounted on the load side of the mechanics. Here, a 6 mm diameter hole is drilled through the outer race of the bearing to simulate a localized fault. Figure 4 shows the healthy and defected bearings used in the experimental setup.Figure 4
Tested bearings: (a) healthy, (b) outer raceway fault.
(a)
(b)In the experiments, two induction motors have been considered: a healthy motor and a faulty one affected by bearing outer raceway fault. Firstly, tests have been performed with faulty motor under several different loads by randomly adjusting the field voltage of the DC machine. Then, supplementary test is performed with the healthy motor for the sake of spectra contrastive analysis. All the tests are done in steady-state conditions. The motors under study are fed by a 50 Hz power supply. The stator current and speed signals are sampled at 50 kHz with a recording length of 25 s.


### 4.2. Experimental Results

Linear plots of the current spectrum when the motor runs at the shaft speed of 1431 r/min (23.85 Hz) are given in Figures5 and 6. Figure 5 shows the current spectrum around the fundamental region while Figure 6 depicts the PSHs (principle slot harmonics, PSHs) region. In each figure, the spectrums from a healthy case and a faulty bearing are provided for comparisons. In the frequency range (0–200 Hz) of the spectrum in the healthy machine (Figure 5(a)), one can see that the frequency components at |f1±fr|, f1+2fr, f1+3fr, and f1+5fr are presented due to the inherent level of eccentricity. Comparisons of Figures 5(a) and 5(b) show that a great number of eccentricity harmonics appeared in Figure 5(b). The magnitudes of eccentricity harmonics are also increased and the frequency component at 3f1+fr can also be observed with the faulty bearing due to the fault-related eccentricity according to (18) and (19) (with nr=nof=0). The fault characteristic frequency at f1+fo is visible and the other fault-related component atf1+fr+focan also be noticed in the faulty spectrum (Figure 5(b)). These results match with existing studies.Figure 5
Comparison of stator current spectrums with 50 Hz supply frequency around fundamental frequency region between healthy IM and IM with outer raceway defect: (a) healthy IM and (b) IM with a bearing outer raceway defect.
(a)
(b)Figure 6
Comparison of stator current spectrums with 50 Hz supply frequency around PSHs region between healthy IM and IM with outer raceway defect: (a) healthy IM and (b) IM with a bearing outer raceway defect.
(a)
(b)A comparison of the stator current spectrums from a healthy case and a faulty bearing around the PSHs region is shown in Figure6, considering the effects of the slot harmonics on the current characteristic frequency in an induction motor with inherent eccentricity and bearing fault. The frequency components at |f1±Rfr| and |3f1-Rfr| are already visible in the spectrum of the healthy machine due to the interaction of stator MMF and rotor slot permeance (Figure 6(a)). With an increase in the 3rd slot time harmonic, the frequency components at |3f1-R+1fr| and |3f1-R+2fr| are also observed with the faulty bearing, because of PSHs interaction with eccentricity harmonics (Figure 6(b)). This is true also of the frequency component at f1+R-1fr. The components at f1-Rfr+fo and f1+Rfr-fo do emerge in the situation of the outer raceway fault, but not in the health spectrum, as they are the result of an intermodulation between PSHs and the characteristic frequencies for bearing fault. Another two fault-related components at f1-R-1fr-fo and f1+R-1fr+fo can also be observed in the case of the fault, owing to the effects of PSHs with eccentricity harmonics and bearing fault characteristic frequencies. The results attained by outer raceway fault prove, hence, earlier theoretical analysis.In order to further prove the effectiveness of the proposed model for bearing fault detection, the current signals under several different load conditions have been processed using MATLAB. The offline experimental results are summarized in Tables2 and 3. Table 2 highlights the presence of the fault frequency components and the magnitudes of these frequencies for healthy and faulty motors under these load conditions. Moreover, some average and standard deviation values of fault frequency components with a number of runs are also provided in Table 3 for the purpose of showing repeatability. In Table 3, there are six runs altogether for the same load conditions. From these experimental results we can clearly see that the proposed approach is reliable to be used for bearing fault detection. The estimated fundamental supply frequency is f1=50.0679Hz in this paper.Table 2
The experimental results of healthy and faulty motors under several different load conditions.
State
f
r
f
o
f
1
+
f
o
f
1
+
f
r
+
f
o
f
1
-
R
f
r
+
f
o
f
1
+
R
f
r
-
f
o
f
1
-
(
R
-
1
)
f
r
-
f
o
f
1
+
R
-
1
f
r
+
f
o
Fault
23.8500
85.8600
Calculated/Hz
135.9279
159.7779
627.2721
727.4079
775.1421
875.2779
Estimated/Hz
135.9463
159.6451
626.4210
727.1767
774.5743
874.6624
Amplitude/dB
−66.9323
−69.9651
−72.6639
−74.7578
−77.3310
−73.2667
Health
23.8500
—
Amplitude/dB
−76.2459
−77.6033
−81.5884
−82.4756
−82.0215
−82.0451
Fault
23.9655
86.2758
Calculated/Hz
136.3437
160.3092
630.5523
730.6881
779.1384
879.2742
Estimated/Hz
136.3277
160.1696
629.2343
730.6576
777.1015
878.1433
Amplitude/dB
−67.8734
−75.5752
−77.3139
−78.9183
−81.8125
−77.7823
Health
23.9655
—
Amplitude/dB
−79.1248
−84.4931
−85.4721
−85.9430
−85.1713
−85.4535
Fault
24.0907
86.7265
Calculated/Hz
136.7944
160.8851
634.1080
734.2438
783.4703
883.6061
Estimated/Hz
136.7569
160.6941
633.7643
733.6617
783.1097
882.5779
Amplitude/dB
−66.9565
−80.2746
−77.7492
−78.7337
−81.7378
−77.8168
Health
24.0907
—
Amplitude/dB
−84.4098
−85.5045
−85.3693
−85.1180
−85.2143
−85.2862
Fault
24.1980
87.1128
Calculated/Hz
137.1807
161.3787
637.1553
737.2911
787.1829
887.3187
Estimated/Hz
137.1384
161.2186
637.2929
736.4750
788.3072
888.9675
Amplitude/dB
−70.0095
−77.2724
−77.8660
−79.1447
−81.7262
−78.4580
Health
24.1980
—
Amplitude/dB
−78.0160
−82.9799
−85.0438
−85.4245
−85.2521
−85.6223Table 3
The average and standard deviation values of fault frequency components under three different load conditions (/Hz).
State
f
r
f
o
f
1
+
f
o
f
1
+
f
r
+
f
o
f
1
-
R
f
r
+
f
o
f
1
+
R
f
r
-
f
o
f
1
-
(
R
-
1
)
f
r
-
f
o
f
1
+
R
-
1
f
r
+
f
o
Fault
23.8500
85.8600
Calculated
135.9279
159.7779
627.2721
727.4079
775.1421
875.2779
Average value
135.9384
159.6292
627.1203
727.4466
774.6935
875.4969
Standard deviation
0.0359
0.0718
0.4097
0.7135
0.3501
0.9974
Fault
23.9655
86.2758
Calculated
136.3437
160.3092
630.5523
730.6881
779.1384
879.2742
Average value
136.3744
160.2650
630.0769
729.7913
779.9864
880.0586
Standard deviation
0.1626
0.3358
0.5213
0.5148
0.6060
0.8409
Fault
24.0907
86.7265
Calculated
136.7944
160.8851
634.1080
734.2438
783.4703
883.6061
Average value
136.7648
160.7179
634.2332
734.4882
782.8474
883.2693
Standard deviation
0.0925
0.0261
0.7259
0.4272
0.7084
0.4389


## 4.1. Test Rig

Schematic diagram of the experiment setup used in this paper is shown in Figure2. The rated data of the used induction motor are given in Table 1. The induction motor is mechanically coupled to a DC machine, which feeds an external resistor bank. The motor can be identically loaded at different speed from zero to rated load by adjusting the field voltage of the DC machine. Both machines are fixed on a steel plate test bench resting on a cement floor. A photoelectric speed sensor (SZGB-6) is used to measure the shaft speed. As the used National Instrument Data Acquisition (NI DAQ 6216) card accepts only voltage signal with maximum amplitudes of ±10 V, JLB-21 Hall current sensor is used to gather the current signal first, and then the acquired current signal is converted to proportional voltage signals by conditioning circuit. A PC equipped with DAQ card through a USB cable is used to store them in the memory. A picture of the experimental setup is shown in Figure 3.Table 1
Test induction motor rated parameters.
Nominal power
3 kW
Nominal rated voltage
380 V
Nominal rated current
6.8 A
Frequency
50 Hz
Pole pairs
2
Stator slots
36
Rotor slots
32
Stator winding connection
YFigure 2
Schematic diagram of the experiment setup.Figure 3
Photograph of the experimental setup used for bearing fault detection.The induction machine has two CU6206RZ type bearings (single row and deep groove ball bearings) with a width ofw = 16 mm, an outside diameter of D = 62 mm, an inside diameter of d = 30 mm, and a pitch diameter of Dc = 46.0 mm. Each bearing has nine balls (Nb = 9) with an approximate diameter of Db = 9.52 mm. The contact angle of the ball with the race (β) is assumed to be 0°. Two identical motors are prepared: one motor carries healthy bearings and the other carries one defected bearing which are mounted on the load side of the mechanics. Here, a 6 mm diameter hole is drilled through the outer race of the bearing to simulate a localized fault. Figure 4 shows the healthy and defected bearings used in the experimental setup.Figure 4
Tested bearings: (a) healthy, (b) outer raceway fault.
(a)
(b)In the experiments, two induction motors have been considered: a healthy motor and a faulty one affected by bearing outer raceway fault. Firstly, tests have been performed with faulty motor under several different loads by randomly adjusting the field voltage of the DC machine. Then, supplementary test is performed with the healthy motor for the sake of spectra contrastive analysis. All the tests are done in steady-state conditions. The motors under study are fed by a 50 Hz power supply. The stator current and speed signals are sampled at 50 kHz with a recording length of 25 s.


## 4.2. Experimental Results

Linear plots of the current spectrum when the motor runs at the shaft speed of 1431 r/min (23.85 Hz) are given in Figures5 and 6. Figure 5 shows the current spectrum around the fundamental region while Figure 6 depicts the PSHs (principle slot harmonics, PSHs) region. In each figure, the spectrums from a healthy case and a faulty bearing are provided for comparisons. In the frequency range (0–200 Hz) of the spectrum in the healthy machine (Figure 5(a)), one can see that the frequency components at |f1±fr|, f1+2fr, f1+3fr, and f1+5fr are presented due to the inherent level of eccentricity. Comparisons of Figures 5(a) and 5(b) show that a great number of eccentricity harmonics appeared in Figure 5(b). The magnitudes of eccentricity harmonics are also increased and the frequency component at 3f1+fr can also be observed with the faulty bearing due to the fault-related eccentricity according to (18) and (19) (with nr=nof=0). The fault characteristic frequency at f1+fo is visible and the other fault-related component atf1+fr+focan also be noticed in the faulty spectrum (Figure 5(b)). These results match with existing studies.Figure 5
Comparison of stator current spectrums with 50 Hz supply frequency around fundamental frequency region between healthy IM and IM with outer raceway defect: (a) healthy IM and (b) IM with a bearing outer raceway defect.
(a)
(b)Figure 6
Comparison of stator current spectrums with 50 Hz supply frequency around PSHs region between healthy IM and IM with outer raceway defect: (a) healthy IM and (b) IM with a bearing outer raceway defect.
(a)
(b)A comparison of the stator current spectrums from a healthy case and a faulty bearing around the PSHs region is shown in Figure6, considering the effects of the slot harmonics on the current characteristic frequency in an induction motor with inherent eccentricity and bearing fault. The frequency components at |f1±Rfr| and |3f1-Rfr| are already visible in the spectrum of the healthy machine due to the interaction of stator MMF and rotor slot permeance (Figure 6(a)). With an increase in the 3rd slot time harmonic, the frequency components at |3f1-R+1fr| and |3f1-R+2fr| are also observed with the faulty bearing, because of PSHs interaction with eccentricity harmonics (Figure 6(b)). This is true also of the frequency component at f1+R-1fr. The components at f1-Rfr+fo and f1+Rfr-fo do emerge in the situation of the outer raceway fault, but not in the health spectrum, as they are the result of an intermodulation between PSHs and the characteristic frequencies for bearing fault. Another two fault-related components at f1-R-1fr-fo and f1+R-1fr+fo can also be observed in the case of the fault, owing to the effects of PSHs with eccentricity harmonics and bearing fault characteristic frequencies. The results attained by outer raceway fault prove, hence, earlier theoretical analysis.In order to further prove the effectiveness of the proposed model for bearing fault detection, the current signals under several different load conditions have been processed using MATLAB. The offline experimental results are summarized in Tables2 and 3. Table 2 highlights the presence of the fault frequency components and the magnitudes of these frequencies for healthy and faulty motors under these load conditions. Moreover, some average and standard deviation values of fault frequency components with a number of runs are also provided in Table 3 for the purpose of showing repeatability. In Table 3, there are six runs altogether for the same load conditions. From these experimental results we can clearly see that the proposed approach is reliable to be used for bearing fault detection. The estimated fundamental supply frequency is f1=50.0679Hz in this paper.Table 2
The experimental results of healthy and faulty motors under several different load conditions.
State
f
r
f
o
f
1
+
f
o
f
1
+
f
r
+
f
o
f
1
-
R
f
r
+
f
o
f
1
+
R
f
r
-
f
o
f
1
-
(
R
-
1
)
f
r
-
f
o
f
1
+
R
-
1
f
r
+
f
o
Fault
23.8500
85.8600
Calculated/Hz
135.9279
159.7779
627.2721
727.4079
775.1421
875.2779
Estimated/Hz
135.9463
159.6451
626.4210
727.1767
774.5743
874.6624
Amplitude/dB
−66.9323
−69.9651
−72.6639
−74.7578
−77.3310
−73.2667
Health
23.8500
—
Amplitude/dB
−76.2459
−77.6033
−81.5884
−82.4756
−82.0215
−82.0451
Fault
23.9655
86.2758
Calculated/Hz
136.3437
160.3092
630.5523
730.6881
779.1384
879.2742
Estimated/Hz
136.3277
160.1696
629.2343
730.6576
777.1015
878.1433
Amplitude/dB
−67.8734
−75.5752
−77.3139
−78.9183
−81.8125
−77.7823
Health
23.9655
—
Amplitude/dB
−79.1248
−84.4931
−85.4721
−85.9430
−85.1713
−85.4535
Fault
24.0907
86.7265
Calculated/Hz
136.7944
160.8851
634.1080
734.2438
783.4703
883.6061
Estimated/Hz
136.7569
160.6941
633.7643
733.6617
783.1097
882.5779
Amplitude/dB
−66.9565
−80.2746
−77.7492
−78.7337
−81.7378
−77.8168
Health
24.0907
—
Amplitude/dB
−84.4098
−85.5045
−85.3693
−85.1180
−85.2143
−85.2862
Fault
24.1980
87.1128
Calculated/Hz
137.1807
161.3787
637.1553
737.2911
787.1829
887.3187
Estimated/Hz
137.1384
161.2186
637.2929
736.4750
788.3072
888.9675
Amplitude/dB
−70.0095
−77.2724
−77.8660
−79.1447
−81.7262
−78.4580
Health
24.1980
—
Amplitude/dB
−78.0160
−82.9799
−85.0438
−85.4245
−85.2521
−85.6223Table 3
The average and standard deviation values of fault frequency components under three different load conditions (/Hz).
State
f
r
f
o
f
1
+
f
o
f
1
+
f
r
+
f
o
f
1
-
R
f
r
+
f
o
f
1
+
R
f
r
-
f
o
f
1
-
(
R
-
1
)
f
r
-
f
o
f
1
+
R
-
1
f
r
+
f
o
Fault
23.8500
85.8600
Calculated
135.9279
159.7779
627.2721
727.4079
775.1421
875.2779
Average value
135.9384
159.6292
627.1203
727.4466
774.6935
875.4969
Standard deviation
0.0359
0.0718
0.4097
0.7135
0.3501
0.9974
Fault
23.9655
86.2758
Calculated
136.3437
160.3092
630.5523
730.6881
779.1384
879.2742
Average value
136.3744
160.2650
630.0769
729.7913
779.9864
880.0586
Standard deviation
0.1626
0.3358
0.5213
0.5148
0.6060
0.8409
Fault
24.0907
86.7265
Calculated
136.7944
160.8851
634.1080
734.2438
783.4703
883.6061
Average value
136.7648
160.7179
634.2332
734.4882
782.8474
883.2693
Standard deviation
0.0925
0.0261
0.7259
0.4272
0.7084
0.4389


## 5. Conclusions

This paper has studied the effects of the slot harmonics on the stator current signals in an induction motor with bearing fault. Based on the airgap flux density calculation by the product of the permeance and the MMF, a more comprehensive mathematical model has been obtained when considering not only inherent eccentricity but also the slot harmonics in an induction motor with bearing fault. Some experimental investigations in this paper have proved that some new characteristic fault frequencies can be found around the rotor slot harmonics in the stator current spectrum. These experimental results agree well with the theoretical inference and the new characteristic frequencies expression can be used as effective assessment indicator for bearing fault detection.



---

*Source: 2640796-2017-04-05.xml*