Method and apparatus for detecting stator faults in rotary dynamoelectric machines

Abstract: In order to determine the occurrence of a stator fault in an induction motor during operation of the motor, the axial leakage flux is monitored during operation of the motor. The occurrence of a stator fault is detected by observing a change in selected harmonic components of the axial leakage. The angular position of the fault with respect to a datum position is determined by noting and recording changes in the axial leakage flux at a plurality of positions distributed symmetrically around the axis of the motor and computing the angular position by triangulation. The axial leakage flux is detected by individual search coils located at said plurality of positions.

Description: This invention relates to a method of detecting and locating stator faults in dynamoelectric machines during operation thereof, and to an apparatus for implementing the method. Specifically, the invention is concerned with the detection and location of a shorted turn of the stator winding by monitoring the axial leakage flux. Studies on reliability and failure mechanisms of dynamoelectric machines with multiturn stator coils indicate that stator insulation breakdown is one of the major causes of failure. Additionally there is evidence that the majority of stator winding failures in such machines result from turn insulation breakdown. In some cases where there has been minimal or no core damage due to the fault, it is possible to isolate the damaged coil from the rest of the winding and restart the machine. Such an operation may necessitate isolation of other coils in order to retain balanced operation of the machine. In this way the machine can be returned to service quickly, albeit with degraded performance until a more complete repair or winding replacement can be scheduled. Unfortunately location of the failed coil is difficult and time consuming. IEEE Proc. B. El. Power Appl., Vol. 133, no. 3, pages 142 to 148 discloses a method of monitoring electric drives in order to pre-empt occurrences of failures and to give an early indication of the onset of possible failures. Fault signatures are used to identify the type of fault which may occur and to give an early indication of the onset of a possible malfunction. Specifically, a possible malfunction is detected by measuring the harmonic components in axial leakage flux and comparing the measured harmonics with known fault signatures. The object of the present invention is to provide a method of, and a device for, detecting the occurrence of a shorted turn and indicating its location in the stator winding of the machine while the machine is still running. According to a first aspect of the present invention, there is provided a method of detecting and locating stator faults in a dynamoelectric machine having a rotor and a stator with multiturn stator coils by monitoring changes in axial leakage flux during operation of the machine, comprising: identifying the frequencies of selected harmonic components of axial leakage flux that are to be significantly affected by the occurrence of a stator fault, detecting the axial leakage flux at each of at least four positions, including a datum position, distributed symmetrically around the axis of the rotor adjacent to an end winding of the stator, deriving from the detected flux at each of said positions a first signal having a value corresponding to the magnitude of the detected flux, storing the values of said first signals derived prior to the occurrence of a stator fault, deriving from the detected axial leakage flux a further signal having a value corresponding to the magnitude of said selected harmonic components, monitoring said further signal and detecting a change therein thereby detecting the occurrence of a stator fault, storing the values of said first signals derived after the occurrence of the stator fault, and comparing the values of said first signals stored before and after the occurrence of a stator fault and storing the differences therebetween, characterized in that the method further comprises: determining from said comparison a value representing the angular position of the stator fault with respect to the datum position by: (i) storing a value R 1 representing the radial distance of said four positions from the axis of the stator, (ii) storing a value R 2 representing the mean radius of the stator end winding, (iii) computing values B 1 and B 2 corresponding respectively to the moduli of said differences between the stored values of said first signals derived from the detected flux at one diametrically opposed pair of said flux detecting positions, (iv) computing values B 3 and B 4 corresponding respectively to the moduli of said differences between the stored values of said first signals derived from the detected flux at the other diametrically opposed pair of said flux detecting positions, (v) computing values ψ 1 and ψ 2 , where ± ψ 1 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 1 -B 2 )/(B 1 +B 2 )] ± ψ 2 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 3 -B 4 )/(B 3 +B 4 )], (vi) determining therefrom the angular position of the stator fault ψ o with respect to said datum position from the relationship ψ o = ψ 1 = ψ 2 + 90° and (vii) displaying the value ψ o . The occurrence of a stator fault representing a shorted turn is thus detected by monitoring selected axial flux harmonics which are sensitive to the occurrence. The angular position of the fault is computed from stored values, namely machine parameters which are initially known and measured values of the flux after detection of the fault. According to a second aspect of the present invention, there is provided an apparatus for detecting and locating stator faults in an induction motor by monitoring changes in axial leakage flux during operation of the motor, comprising: means for detecting the axial leakage flux during operation of the motor, said means including means for deriving a plurality of first signals each having a value corresponding to the detected flux at a respective one of at least four positions, including a datum position, distributed around the axis of the motor, means for deriving from the detected flux a further signal corresponding to the magnitude of the detected flux at selected harmonic frequencies, first data processing means for monitoring said further signal and for detecting the occurrence of a stator fault by a change in the value thereof, and second data processing means responsive to said first data processing means, said second data processing means including: (a) means for storing said first signals prior to the occurrence of a stator fault, and (b) means responsive to the detection of a stator fault for storing said first signals after the occurrence of the fault, characterized in that said second data processing means further includes: (c) means for computing from said stored values of said first signals the angular position of the stator fault with respect to said datum position, (d) said means for computing including: (d1) means for storing a value R 1 representing the radial distance of said four positions from the axis of the stator, (d2) means for storing a value R 2 representing the mean radius of the stator end winding, (d3) means for computing values B 1 and B 2 corresponding respectively to the moduli of said differences between the stored values of said first signals derived from the detected flux at one diametrically opposed pair of said flux detecting positions, (d4) means for computing values B 3 and B 4 corresponding respectively to the moduli of said differences between the stored values of said first signals derived from the detected flux at the other diametrically opposed pair of said flux detecting positions, (d5) means for computing values ψ 1 and ψ 2 , where ± ψ 1 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 1 -B 2) /(B 1 +B 2 )] ± ψ 2 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 3 -B 4 )/(B 3 +B 4 )], (d6) means for determining therefrom the angular position of the stator fault ψ o with respect to said datum position from the relationship ± ψ o = ψ 1 = ψ 2 + 90° and (d7) displaying the value ψ o . The preferred method according to the invention, and an apparatus therefor, will now be described by way of example with reference to the accompanying drawings. In the drawings: Figure 1 is a graph showing the waveform of detected magnetomotive force (MMF) due to a shorted stator turn; Figures 2 and 3 are diagrams illustrating the geometric considerations required in magnetic field calculation; Figure 4 is a diagram illustrating the geometry of sensors used to detect axial leakage flux; Figure 5 is a partially diagrammatic illustration of an apparatus used to monitor axial leakage flux; Figure 6 is a diagram illustrating the general algorithm used to detect and locate a stator fault; and Figure 7 is a diagram illustrating the detailed algorithm used to locate the stator fault. Basic Principles In the case of an ideal electrical machine there is zero net axial magnetic leakage flux because the stator and rotor currents should exactly cancel one another. However, in practice, due to nonuniformities in materials and construction methods, such a situation does not occur. Hence associated with all rotating machines there is a small, but measurable, axial leakage flux. This magnetic flux can be measured by the simple expedient of placing a search coil, or other appropriate magnetic flux sensing device, in close proximity to the end region of the machine. The output from such a device, when analyzed in the frequency domain, is manifest as a series of harmonics of varying magnitude. From basic electromagnetic considerations of electrical machines it is possible to derive expressions which permit the harmonic components inherent in the axial leakage flux to be predicted. The equations, which are simple algebraic expressions, only require knowledge of the number of poles and rotational speed of the machine, in order to calculate all possible harmonics for a particular machine. The occurrence of a stator fault in an induction motor, for example, will inevitably result in a change in the air gap space harmonic distribution. These space harmonics cannot be detected directly by a search coil. However, the search coil can detect the time harmonics of the axial flux. Thus it is necessary to derive the relationship between the space and time harmonics in order to correctly interpret the frequency spectrum obtained from the search coil. For the purpose of this work, considering only the supply fundamental and the third harmonic component due to saturation, the space harmonic distribution of MMF due to a balanced, full pitched, three phase winding fed from a balanced supply frequency, ω, is given by [1], m = 0.955N 1 [k w1 cos(ωt-p)+0.2k w5 cos(ωt+5p) -0.14k w7 cos(ωt-7p)+0.09k w11 cos(ωt+11p)-...] where k wn is the n th winding factor p is the number of pole pairs  is the angular displacement from the stator datum This represents a rotating set of harmonics of order 6n±1 which can be simplified to the corresponding air gap fluxes, B s = B 1 cos(ωt-p)+B 5 cos(ωt+5p)-B 7 cos(ωt-7p) +B 11 cos(ωt+11p) where B n are the spatial harmonic fluxes. This expression is in the stator frame of reference. Consequently, because the shaft flux of the rotor is of interest, it is necessary to refer equation (2) to the rotor reference frame. Consider the situation in which β is the angular displacement between the rotor and stator datum positions, and a is defined to be the angular displacement from the rotor datum. Then  = α + β. If the angular rotor speed is ω r , then,  = α+ω r t now using the normal expression for the slip of the motor, i.e. s = (ω s -ω r )/ω s , where ω s , the synchronous speed, = ω/p, ω r = ω(1-s)/p Now, the general term of equation (2) is, B ns = B n cos(ωt±np) substituting equations (3) and (4) into (5) leads to, B ns = B n cos[(1±n(1-s))ωt±npα] Expanding the expression for the first few terms gives, B s = B 1 cos(sωt-pα)+B 5 cos[(6-5s)ωt+5pα] -B 7 cos[7s-6)ωt-7pα]+B 11 cos[(12-11s) + 11 pα]- Equation (7) gives the frequency components of the currents that are induced in the rotor due to the air gap space harmonics of a balanced winding and supply. In addition to these harmonics, the fundamental of the supply frequency will also appear in the axial flux spectrum. The presence of additional, higher order, harmonics can be accounted for by using the term nω rather than ω. Having derived the harmonic components of axial flux due to asymmetries in a "healthy" machine, the situation in which an interturn short circuit is introduced is now considered. The effect of the interturn fault is to remove a turn from the stator winding. This will have a small, but finite, effect on the main air gap flux distribution. In addition, an emf will be induced in the shorted turn which will result in a current flow limited only by the self impedance of the fault. This impedance essentially determines the transition time between turn and groundwall insulation failure. The fault current due to the shorted turn is the source of an additional MMF pulse. This MMF pulse also has a space harmonic distribution which is superposed on the main field distribution. From previous considerations, this will lead to a change in the time harmonics observed in the leakage field. The changes to be expected can also be predicted mathematically, as follows: Simple consideration of the MMF distribution due to an interturn short circuit leads to the characteristic illustrated in Figure 1. This is the case for a four pole machine. The analysis to be outlined below is for the general case of a 2p pole machine. Fourier analysis of the MMF waveform, illustrated in Figure 1, shows that it contains all harmonics except the fourths, i.e., B s = 0.5ΣB n cos(ωt±n),n ≠ 4m, all m For the general case, the corresponding waveform would have a mark-space ratio of 1:(2p-1) causing every 2p th harmonic to be absent. Hence the time harmonics produced by the rotor are given by, B s = 0.5ΣB n cos[(1±n(1-s)/p)ωt±nα], n ≠ 2pm, Adding in supply time harmonics of higher order k, leads to the completely general expression, B s = 0.5ΣΣB n cos[(k±n(1-s)/p)ωt±nα], n ≠ 2pm, Although this is a large series, only the lower order harmonics are significant. The key element of this expression is, [k±n(1-s)/p]f 1 for k = 1,3 and n = 1,2,3,.........,(2p-1) The nα term in the argument of (10) causes the components defined above to beat at the slip frequency of the rotor current. Thus the following procedure can be employed to identify the stator fault: 1. Calculate air gap space harmonics due to balanced supply and winding conditions. 2. Derive associated time harmonic currents (fluxes) in rotor and stator reference frames. 3. Refer all time harmonic components to the stator frame, since the search coil is stationary. These are the components to be expected in the "healthy" machine. 4. Calculate additional space harmonics injected by the occurrence of an interturn short circuit. 5. Relate these to additional time harmonics in the spectrum of the axial field. 6. Look for changes in these components, which can be predicted for a given machine, to indicate the presence of an interturn fault. Having determined that a stator fault condition exists, the next task is to locate the position of the fault within the stator winding. The introduction of a shorted turn in the stator winding produces an asymmetry in the magnetic field in the endwinding of the machine. This is because the short circuit current flowing in the fault is not balanced by an equal current flowing in the corresponding phase belt which is situated diametrically opposite the faulted turn. Thus, by measuring the asymmetry in the endwinding magnetic field the position of the faulted coil within the stator winding can be determined. The technique relies upon the use of an array of at least four search coils, or other appropriate sensors, distributed axisymmetrically about the drive shaft of the motor. These coils provide a local measure of the magnetic field in the end region of the machine. By triangulating the outputs from the search coils, it is possible to determine the location of the stator coil containing the interturn short circuit. In order to accomplish this function it is necessary to develop an expression for the field at any point on the circumference of the circle passing through the centres of the search coils due to the current flowing in the end winding of an arbitrarily positioned turn of the stator winding. By invoking the Biot-Savart law and with reference to Figure 2, one can deduce d B = (µ o Id l x a z )/4πr 2 where B is the magnetic field strength vector µ o is the permeability of free space a z is the unit vector in the z direction From Figure 2, and rewriting equation (12) dB = (µ o Idlsin)/4πr 2 but since dlsin = rd, dB = (µ o Id)/4πr Figure 3 illustrates the general case in which the shorted turn is displaced an angular amount ψ, from an arbitrary datum position. In this case R 2 is the mean radius of the end winding, and R 1 is the radius of the circle on which the field is to be calculated. The circles are assumed to be coplanar and concentric. The length of arc of the coil is L and the field measurement point is displaced an angle δ from the datum. Using the cosine rule, r 2 = R 1 2 +R 2 2 -2R 1 R 2 cos(ψ+δ) If R 2 is very much greater than R 1 , then L ≈ r, thus  = L/[R 1 2 +R 2 2 -2R 1 R 2 cos(ψ+δ)] 0.5 From equations (14) and (16), the value of B at the field point is given by, B = µ o IL/[4π(R 1 2 +R 2 2 -2R 1 R 2 cos(ψ+δ))] By fixing the position of a number of coils with respect to the datum, and arbitrarily assigning a reference coil at δ = 0, the value of B at each coil can be found as a function of ψ. From the measured values of B, ψ can be found and consequently the position of the faulted turn. In order to optimize the detection scheme, it is necessary to determine the minimum number of coils required for simple and reliable fault location. Consider the situation in which two coils are spaced 180° apart. Each coil will have an induced emf, proportional to the flux linking it. With reference to Figure 4, and taking the coil at position 1 to be the datum, from equation (17), B 1 = k/(a-bcos ψ) and B 2 = k/(a-bcos(ψ+180°)) thus B 2 = k/(a+bcosψ) where k = µ o IL/4π a = (R 1 2 +R 2 2 ) b = 2R 1 R 2 From equations (18) and (20), and some rearrangement of the terms, cosψ = [a(B 1 -B 2 )]/[b(B 1 +B 2 )] This expressions locates the fault to within ±ψ because, ±ψ = cos -1 [ 0.5(R 1 /R 2 +R 2 /R 1 ){(B 1 -B 2 )/(B 1 +B 2 )}] In practice the measured values of the emf induced in the coils can be used because the B values are ratioed in equation (22). Consequently the values used when attempting fault location should be the modulus of the difference between the "healthy" and "faulted" conditions. In order to locate uniquely the fault, a second set of coils, 3 and 4, is required, as shown in Figure 4. For convenience they are displaced at right angles to the axis of symmetry of coils 1 and 2. However, the angular displacement of the coils is not critical. Repeating the location procedure in the same fashion as above with the second set of coils, will yield ±ψ p2 , which differentiates it from ±ψ p1 , found from coils 1 and 2. If the higher magnitude coil signal is chosen as the point from which to measure, then ψ will always fall within the range of 0° to 90°. In addition, only one quadrant will contain an angular position identified using both sets of coils. This is the approximate location of the fault. In order to clarify the procedure for locating the position of a faulted coil within the stator winding, reference should be made to Figure 4. For the purpose of the location algorithm, the coils are grouped in pairs. In this case coils 1 and 2 form one pair, and coils 3 and 4 form the other pair. Initially, at a point in time at which the motor is considered to be healthy, the induced emf value on each of the coils is measured and recorded. In the event of an interturn short circuit occurring, the emf values on each of the coils are again recorded. These latter values are subtracted from the initial values and the modulus of the resultant is taken. Subsequently, the four positive numbers thus derived, i.e. B 1 , B 2 , B 3 and B 4 , are input into the fault location algorithm. This operation results in two expressions, ± ψ p1 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )((B 1 -B 2 )/(B 1 +B 2 ))] and, ±ψ p2 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )((B 3 -B 4 )/(B 3 +B 4 ))] The resultants from these expressions are angles. The reference points are taken as the coils, in each pair, upon which the largest emf value was recorded. Thus, if for example coil 1 had the largest reading between coils 1 and 2, one would say that the position of the fault was ±ψp 1 from coil 1. Similarly, if the reading on coil 4 was the larger of the readings on coils 3 and 4, the fault would be located ±ψ p2 from coil 4. Hence, it is a simple step to draw the conclusion that, in this case, the the fault must lie in the quadrant between coils 1 and 4 at the angular position ψ p1 -ψ p2 /2. To summarize,the detection and location scheme is as follows: from machine parameters calculate the axial leakage flux components under normal conditions, calculate which of the harmonics will change as the result of an interturn short circuit search for these components in the axial leakage field of the machine, and upon fault indication apply the fault location technique described above. Implementation of the Method Figure 5 illustrates a preferred apparatus used to detect and locate the occurrence of a shorted turn of the stator winding 11 of an induction motor 10 by monitoring changes in the axial leakage flux during operation of the motor. The apparatus comprises two subsystems, namely means for sensing the axial leakage flux and means for processing the signals derived therefrom. A further subsystem, which may be necessary in some cases, is some form of speed measurement in order to gain a measure of the slip of the motor when loaded. The flux sensors shown in Figure 5 are air cored coils 13, which are inexpensive and easily fabricated. The preferred location of the search coils is axisymmetric with the drive shaft of the motor, either inside or outside the motor casing. In the illustrated embodiment there are four search coils 13 distributed symmetrically around the axis of the motor adjacent to an end winding of the stator 11. In principle there could be more than four search coils, but the use of four search coils arranged as diametrically opposed pairs simplifies the signal processing procedures and circuitry. Speed measurement, if required, can be accomplished in a number of ways. In some cases, speed information is available at the motor control centre in the plant. However, it is possible that some type of transduction system may be necessary. A measure of the rotational speed of the rotor is needed in order to calculate the slip of the motor. The slip is required because the harmonics of the axial leakage flux are dependent on the slip of the motor. Rotational speed can easily be obtained by measuring the time required for a reference point on the drive shaft to complete one revolution. In practice this can be done by painting a stripe on the drive shaft and using a device such as a photo transistor to generate an electrical pulse each time the stripe passes under the photo transistor. It is necessary to condition the voltage signals derived from the search coils, which are normally of only a few millivolts. In order to bring these signals to more manageable levels they are amplified by preamplifiers 14. Furthermore, in order to remove the power frequency component, which is relatively high and likely to saturate subsequent amplifier stages, the signals are filtered by power frequency notch filters 15. Alternatively, very wide dynamic range amplifiers may be used. The four analogue signals, after conditioning, are digitized to produce respective "first" signals having magnitudes corresponding to the harmonic flux detected by the four search coils 13. These "first" signals are subsequently processed by a computer 16 as hereinafter described, to determine the angular position of a detected stator fault. The analogue signals from the notch filters 15 are also summed by a summing amplifier 17, the summed signal then being digitized. This summed signal must now be decomposed into its harmonic components. Two practical approaches to this task are possible. Given that the components of the axial leakage flux are mathematically predictable, a set of filters, tuned to the frequencies of interest, can be used. The outputs of these filters can be monitored for changes which potentially indicate the onset of a fault condition. In practice, in view of the variety of motor designs which exist, construction of discrete filters tailored to the motor parameters is not a practical approach. However, by using digitally controlled analogue filters, it is possible to set up an array of filters which can be readily customized to the particular motor. For example, a single controllable filter may be sequentially stepped through the predetermined frequencies of interest. The preferred method, according to the present invention, of decomposing the signal into its harmonic components is to digitize the signal from the summing amplifier 17 and to transform the digital signal to the frequency domain by means of a fast Fourier transform (FFT). In this way, a further signal having a value corresponding to the magnitude of selected harmonic components of axial leakage flux is derived. The harmonic components are selected as being those which are significantly affected by the occurrence of a stator fault. The preferred method will now be described with reference to Figures 6 and 7. In order to detect the occurrence of a stator fault by monitoring changes in the axial leakage flux, it is first necessary to identify the frequencies of the harmonic components which are significantly affected by the occurrence. These can be calculated from the machine parameters and the conditions under which the motor is operating. Referring to Figure 6, the first thing to do is to input into the computer the relevant input parameters, specifically the number of poles of the motor and its nominal speed (BLOCK 21). From these parameters the harmonic components of axial leakage flux under faulted and no-fault conditions are computed (BLOCK 22). If the motor is on-load (BLOCK 23), the rotor velocity has to be taken into account (BLOCK 24); in either case the results of the calculation are input into the initialization procedure (BLOCK 25). The summed outputs from the search coils are recorded at a first period t(n) and converted to the frequency domain, → A(f,t(n)) (BLOCK 26) as previously described. The results are tested (BLOCK 27) and should confirm the calculations of (BLOCK 23). The next step in the procedure is to measure and record the summed output of induced emf on each of the search coils at a second period t(n+1), and converted to the frequency domain, → A(f,t(n+1)) (BLOCK 29).L If a possible fault is detected by applying the test A(f,t(n+1)) >> A(f,t(n)) then the summed output from the search coils is measured and recorded at a third period and converted to the frequency domain, →A(f,t(n + 2)) (BLOCK 31). The possible fault is confirmed or ignored by testing whether A(f,t(n + 2)) >> A(f,t(n)) If the result is positive, the fault is confirmed; otherwise it can be taken that the apparent fault was due to a temporary or spurious event. If the fault is confirmed, the rms value of induced emf in each of the search coils in the faulted condition of the motor is measured and recorded, (BLOCK 33). The emf value for the faulted condition is subtracted from the stored emf value for the healthy condition, for each search coil, and the modulus of the difference is taken (BLOCK 34). From the moduli value so obtained the angular position of the fault with respect to a datum position can be calculated, (BLOCK 35), as will be described with reference to Figure 7. The occurrence of the stator fault, and its position in the stator winding, are then displayed (BLOCK 36). Upon confirmation of the occurrence of a stator fault (BLOCK 32), the angular position of the fault is determined by using the algorithm illustrated in Figure 7. This algorithm involves the use of values derived from signals obtained after the occurrence of the fault, as will now be described, as well as the following values: R 1 = radius of the circle on which the axes of the search coils lie, this value having been entered and stored initially; R 2 = mean radius of the end winding of the stator, this value having been entered and stored initially; EMF (healthy) = the values of said "first" signals derived from the search coils, these values having been stored prior to the occurrence of the fault. Referring to Figure 7, the number of search coils n was previously entered into the computer and stored (BLOCK 41). In the present example n = 4, and in the algorithm the individual search coils are identified by the letter x, where x = 1, x = 2, x = 3 and x = 4, respectively. The signals from the search coils, previously referred to as the "first" signals, are examined in turn, (BLOCKS 42, 43, 44). The "first" signal from each of the search coils corresponds to the axial leakage flux detected by the search coil, being derived from the voltage signal which is a measure of the rms value of the emf induced in the coil. The value of this signal, denoted by EMF (faulted, position x), is subtracted from the previously stored value of the signal EMF (healthy, position x), (BLOCK 45), and the modulus |B(x)| is taken (BLOCK 46). The computer now has in store all the necessary values and parameters from which the angular position of the stator fault can be computed, namely R 1 , R 2 , B 1 , B 2 , B 3 , B 4 . From these values the following relationships are derived, (see BLOCK 47): ± ψ 1 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 1 -B 2 )/(B 1 +B 2 )] ± ψ 2 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 3 -B 4 )/(B 3 +B 4 )] Since ψ = ψ 2 + 90°, the required angular position of the fault ψ 0 can be determined. In other words the quadrant in which the fault lies can be readily determined. Thus, considering the first pair of diametrically opposed search coils 1, 2, one must determine whether EMF(fault, 1) is greater or less than EMF(fault, 2), (BLOCK 48). If it is greater, then the fault is at one of the positions ±ψ 1 from coil 1,(BLOCK 49). If it is less, then the fault is at one of the positions ±ψ 1 from coil 2, (BLOCK 50). In order to resolve the ambiguity, one must now determine whether EMF(fault, 3) is greater or less than EMF(fault, 4), (BLOCK 51). If it is greater, then the fault is at a position ±ψ 2 from coil 3, (BLOCK 52), while if it is less the fault is at one of the positions ±ψ 2 from coil 4, (BLOCK 53). The fault obviously lies within the quadrant defined by ±ψ 1 and ±ψ 2 , (BLOCK 54); that is to say, the angular position ψ 0 of the fault with respect to the datum position is given by ψ o = ψ 1 = ψ 2 + 90°. This position is displayed along with indication that a fault has occurred.

Claims: A method of detecting and locating stator faults in a dynamoelectric machine having a rotor and a stator with multiturn stator coils by monitoring changes in axial leakage flux during operation of the machine, comprising: identifying the frequencies of selected harmonic components of axial leakage flux that are to be significantly affected by the occurrence of a stator fault, detecting the axial leakage flux at each of at least four positions, including a datum position, distributed symmetrically around the axis of the rotor adjacent to an end winding of the stator, deriving from the detected flux at each of said positions a first signal having a value corresponding to the magnitude of the detected flux, storing the values of said first signals derived prior to the occurrence of a stator fault, deriving from the detected axial leakage flux a further signal having a value corresponding to the magnitude of said selected harmonic components, monitoring said further signal and detecting a change therein thereby detecting the occurrence of a stator fault, storing the values of said first signals derived after the occurrence of the stator fault, and comparing the values of said first signals stored before and after the occurrence of a stator fault and storing the differences therebetween, characterized in that the method further comprises: determining from said comparison a value representing the angular position of the stator fault with respect to the datum position by: (i) storing a value R 1 representing the radial distance of said four positions from the axis of the stator, (ii) storing a value R 2 representing the mean radius of the stator end winding, (iii) computing values B 1 and B 2 corresponding respectively to the moduli of said differences between the stored values of said first signals derived from the detected flux at one diametrically opposed pair of said flux detecting positions, (iv) computing values B 3 and B 4 corresponding respectively to the moduli of said differences between the stored values of said first signals derived from the detected flux at the other diametrically opposed pair of said flux detecting positions, (v) computing values ψ 1 and ψ 2 , where ± ψ 1 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 1 -B 2 )/(B 1 +B 2 )] ± ψ 2 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 3 -B 4 )/(B 3 +B 4 )], (vi) determining therefrom the angular position of the stator fault ψ o with respect to said datum position from the relationship ψ o = ψ 1 = ψ 2 + 90° and (vii) displaying the value ψ o . A method according to claim 1, wherein the storage of said first signals derived after the occurrence of a stator fault is initiated by the detection of the fault. A method according to claim 1 or 2, wherein the axial leakage flux is detected at each of four positions, including said datum position, distributed symmetrically around said axis. A method according to any one of claims 1 to 3, wherein the spectrum of axial leakage flux is analysed by deriving a respective voltage signal corresponding to the detected leakage flux at each of said plurality of positions, summing the derived voltage signals, digitising the sum of said voltage signals, and transforming the digitised signal to the frequency domain by means of a fast Fourier transform. Apparatus for detecting and locating stator faults in an induction motor (10) by monitoring changes in axial leakage flux during operation of the motor, comprising: means (13) for detecting the axial leakage flux during operation of the motor, said means including means for deriving a plurality of first signals each having a value corresponding to the detected flux at a respective one of at least four positions, including a datum position, distributed around the axis of the motor, means for deriving from the detected flux a further signal corresponding to the magnitude of the detected flux at selected harmonic frequencies, first data processing means (16) for monitoring said further signal and for detecting the occurrence of a stator fault by a change in the value thereof, and second data processing means (16) responsive to said first data processing means, said second data processing means including: (a) means for storing said first signals prior to the occurrence of a stator fault, and (b) means responsive to the detection of a stator fault for storing said first signals after the occurrence of the fault, characterized in that said second data processing means further includes: (c) means for computing from said stored values of said first signals the angular position of the stator fault with respect to said datum position, (d) said means for computing including: (d1) means for storing a value R 1 representing the radial distance of said four positions from the axis of the stator, (d2) means for storing a value R 2 representing the mean radius of the stator end winding, (d3) means for computing values B 1 and B 2 corresponding respectively to the moduli of said differences between the stored values of said first signals derived from the detected flux at one diametrically opposed pair of said flux detecting positions, (d4) means for computing values B 3 and B 4 corresponding respectively to the moduli of said differences between the stored values of said first signals derived from the detected flux at the other diametrically opposed pair of said flux detecting positions, (d5) means for computing values ψ 1 and ψ 2 , where ± ψ 1 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 1 -B 2 )/(B 1 +B 2 )] ± ψ 2 = cos -1 [0.5(R 1 /R 2 +R 2 /R 1 )•(B 3 -B 4 )/(B 3 +B 4 )], (d6) means for determining therefrom the angular position of the stator fault ψ o with respect to said datum position from the relationship ± ψ o = ψ 1 = ψ 2 + 90° and (d7) displaying the value ψ o . Apparatus according to claim 5, wherein said means for detecting the axial leakage flux comprises four flux sensing coils (13) mounted symmetrically around the axis of the stator (11) adjacent to an end winding thereof.

Inventors: LLOYD BLAKE A, PENMAN JAMES, SEDDING HOWARD G, LLOYD, BLAKE A., PENMAN, JAMES, SEDDING, HOWARD G.
Applicants: ONTARIO HYDRO
Application year: 1993
Main classification: G01R31
IPC codes: G01R31