This invention relates to analysis of signals. More particularly, this invention relates to a system for separating periodic amplitude peaks from nonperiodic amplitude peaks in machine vibration data.
Periodicity is a term for quantifying how much energy in a signal is periodic. The autocorrelation coefficient function is a mathematical process that can be used to determine how much of the energy in a signal waveform is periodic. The pattern of the periodic peaks—or lack thereof—can be very helpful in identifying fault types. The term “periodic peaks” is a general characterization describing one or more family of harmonic peaks. Periodic peaks related to running speeds of a machine are referred to as synchronous periodic peaks. A family of periodic peaks that are not harmonically related to the running speed (such as bearing and belt frequencies) is referred to as non-synchronous periodic peaks.
Traditionally, recognizing these patterns and how to apply them requires an experienced vibration analyst. The desire is to first calculate an indicator that is representative of general periodic patterns. This value, referred to herein as a “periodic signal parameter” (PSP), is calculated based on statistical measures derived from an autocorrelation waveform of the originating waveform. Based on the peaks in the originating spectrum in combination with the peaks from the spectrum taken from autocorrelation waveform and the PSP value, a Periodic Information Plot (PIP) can be produced. This plot contains only the periodic peaks (synchronous and non-synchronous periodic peaks) from the originating spectra.
An estimate of how much periodicity is present in the original waveform can be determined from the largest amplitude found in the autocorrelated waveform. The square root of the largest amplitude in the autocorrelated waveform is used to estimate Percent Periodic Energy (defined as a percentage of energy in the original spectra that is related to periodic signals). In general, if the value of the PSP is greater than 0.1, or the Percent Periodic Energy is greater than 50%, a periodic signal is present.
By knowing the estimated energy in the original spectrum related to periodic peaks, a “basic” Periodic Information Plot (PIP) can be generated. This “basic” PIP is created by plotting the periodic peaks one peak at a time—starting with the largest peak and continuing in descending order—until the energy in the resulting plot is equal to the Percent Periodic Energy. However, since the Percent Periodic Energy value is an estimate, the resulting “basic” PIP often contains excessive noise.
Described herein is a statistical method to more precisely determine which peaks in the autocorrelation spectrum are periodic. Once the precise periodic peaks are obtained, a “cleaner” (less noisy) PIP can be constructed. The PIP algorithm is described in more detail hereinafter.
The statistical method for separating periodic from non-periodic data in autocorrelation spectra has been developed based on an understanding of the distribution associated with each set of peaks. As a whole, the autocorrelation spectrum is not normally distributed because the periodic peaks are significantly large and random in amplitude relative to the noise floor. However, the noise floor is generally normally distributed (assuming Gaussian noise). Based on this, a method has been devised to separate the periodic peaks from the noise (non-periodic peaks).
A statistical parameter that best characterizes the noise floor is kurtosis. For a normally distributed signal, the kurtosis should have a value of 3. Therefore, after sorting peaks from largest to smallest amplitudes, the kurtosis of the set of peaks is calculated, and the largest peaks are sequentially removed from the set until the kurtosis is equal to a kurtosis threshold. In a preferred embodiment, the kurtosis threshold is 3. In other embodiments, the kurtosis threshold may be 2.9, or 3.1, or some other value near 3. As each large peak is removed from the peak set, the cumulative energy of the remaining peaks in the set is calculated. All the peaks that were culled out (forming a periodic peak set) are considered to be periodic. The energy of the periodic peak set is considered to be the Periodic Energy of the original spectrum. As the process of culling peaks proceeds, if the total energy of the “growing” periodic peak set becomes greater than or equal to the Percent Periodic Energy before the peak set kurtosis reaches a kurtosis threshold of 3 or less, then the process stops, and the periodic peak set is thus defined.
The kurtosis algorithm described herein can be applied to any FFT, not just autocorrelated FFT's. By applying this algorithm to a velocity FFT for instance, the major located peaks can be culled out and used for determining the most influential peaks in an FFT. (These peaks will not necessarily be all the periodic peaks such as can be found when applied to the current disclosure.) These influential peaks can be used as a basis for improving speed detection as well as major harmonic families. The harmonic families can be related to mechanical components generating gear mesh, gear teeth, eccentricity (static and dynamic) and bearings anomalies to name a few.
Implementing the periodicity technology described herein is particularly useful when applied to PeakVue™ data. “PeakVue Plus” refers to a process that combines the PeakVue™ process with the periodicity described herein. Any periodicity that is not a harmonic of a running speed (referred to herein as “non-synchronous periodicity”) is typically associated with a bearing fault (such as inner or outer race faults or ball spin or cage faults). The severity of a bearing fault may be determined by the peak amplitude of the associated PeakVue™ waveform. This severity is proportional to the fault levels dictated by the turning speed of the bearing. When a gearbox is under surveillance, any synchronous periodicity is related to the gear teeth health. The severity of gear teeth faults is related to the peak amplitude of the PeakVue™ waveform and is proportional to the fault levels dictated by the associated gear speed. When large peak amplitude values are present in the PeakVue™ waveform and PSP≤0.1, then lubrication issues are suspected in roller element bearings and/or gearboxes.
By determining the periodicity of a waveform, the PIP can be produced to allow an analyst to easily visualize the condition of the machinery being monitored. From this simplified plot, the analysist can predict faults to be acted upon or investigated if desired. Preferred embodiments of a process for generating the PIP is described hereinafter.
Some embodiments described herein provide a computer-implemented method for analyzing periodic information in digital vibration data to determine a rotational speed of a shaft in a machine. A preferred embodiment of the method includes:
In some embodiments, the kurtosis threshold is 3.
In some embodiments, the vibration spectrum is a velocity spectrum.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee. Further advantages of the invention are apparent by reference to the detailed description in conjunction with the figures, wherein elements are not to scale so as to more clearly show the details, wherein like reference numbers indicate like elements throughout the several views, and wherein:
In an alternative embodiment depicted in
With regard to sensor placement for bearing and gear diagnosis, the sensor 104 is typically mounted orthogonal to the shaft. It is preferably mounted on a rigid and massive piece of metal that is near the source of the signal (i.e., bearing or gear). The large mass of metal on which the sensor is mounted helps prevent resonances entering the signal due to the surface of the machine as opposed to what is happening internal to the machine. The sensor 104 should be mounted so as to minimize loss of signal integrity during transmission. This requires a rigid connection—typically by stud mounting the sensor 104. In some circumstances, such as where the mounting surface of the machine is rough or covered with many layers of paint, the surface will need to be sanded.
Periodic Signal Parameter
If MaxPeak is greater than or equal to 0.3 (step 20) and
If MaxPeak is greater than or equal to 0.3 (step 20) and
then Y=0 (step 25).
If MaxPeak is less than 0.3 (step 20) and CF1 less than 4 and σ is less than or equal to 0.1 (step 26), then Z=0.025 (step 28). If MaxPeak is less than 0.3 (step 20) and CF1 is not less than 4 or σ is greater than 0.1 (step 26), then Z=0 (step 30).
If CF2 is greater than or equal to 4 and the number of discarded peaks is greater than 2 (step 36), then W=0.025 (step 38). If CF2 is less than 4 or the number of discarded peaks is not greater than 2 (step 36), then W=0 (step 40).
If
and σ is between 0.1 and 0.9 (step 44), then X=0.1 (step 46). If
or σ is not between 0.1 and 0.9 (step 44), then X=σ (step 48).
The PSP is the sum of the values of X, W, Y and Z (step 50).
In general, smaller PSP values are indicative of more non-periodic signals and less distinctive frequencies, while larger PSP values are symptomatic of more periodic signals relating to large single frequencies. As shown in
Following are some advantages of generating a PSP.
A further embodiment of the present invention employs a programmable central processing unit, such as the processor 114, programmed with program logic to assist a user with an interpretation of waveform information. The program logic compares the Periodic Signal Parameter and Stress Wave analysis information with expected or historical or empirically-derived experiential values to discern a relative ranking from low to high. Then discrete or graduated outputs, such as those portrayed in Table 2 above, are employed to select logically arrayed observations, findings, and recommendations. In addition to evaluating PSP and Stress Wave Analysis information, program logic sometimes prompts a user to supply additional information or obtains additional information from another source such as from a knowledge base, to enable the logic to distinguish between two or more possible logical results. For example, program logic that returns a high PSP and a high Stress Wave Analysis finding may select a rolling element defect finding rather than other possible findings within that category because a similarity is calculated when program logic compares a periodic frequency finding and a bearing fault frequency for a machine component identified in a knowledge base.
Another technique to differentiate between lubrication and pump cavitation is to look at the trend of the impacting as indicated by Stress Wave analysis. If it increases slowly, then insufficient lubrication should be suspected. If it increases suddenly on a pump, then it is likely pump cavitation. If combined with logic or inputs on a control system, then the logic could look for process configuration changes that occurred at the same time as the increase in impacting—along with a low PSP—to confirm pump cavitation. In some embodiments, the system suggests to the operator what action caused the cavitation, so that the operator can remove the cause and stop the machine from wearing excessively and failing prematurely.
Periodic Information Plot
A preferred embodiment of the invention creates a new type of vibration spectrum, referred to herein as a Periodic Information Plot (PIP). The PIP provides the user an easily viewed summary of the predominate periodic peaks from the originating spectrum, which would be a PeakVue spectrum in a preferred embodiment.
In a first embodiment, a signal is collected from plant equipment (e.g. rotating or reciprocating equipment) and is processed using two different sets of analysis techniques as depicted in
First, a waveform is acquired (step 60 of
The waveform from step 60 is also autocorrelated (step 66) to generate a waveform referred to herein as the autocorrelation waveform 68, having time on the X-axis and the correlation factor on the Y-axis. The autocorrelation process accentuates periodic components of the vibration waveform, while diminishing the presence of random events in the original signal. As a result of the autocorrelation calculations, the autocorrelation waveform 68 has half the x-axis (time) values as that of the original vibration waveform 60. Therefore, the timespan of the autocorrelation waveform 68 will be half of that of the original vibration waveform 60. An optional step (70) takes the square root of the autocorrelation waveform (Y-axis values) to provide better differentiation between lower amplitude values.
An FFT of the autocorrelation waveform 68 is taken (step 72), resulting in an autocorrelation spectrum (AS) 74. Since random events have largely been removed from the autocorrelation waveform 68, the remaining signal in the autocorrelation spectrum 74 is strongly related to periodic events. As shown in
In the first embodiment, the vibration spectrum 64 and the autocorrelation spectrum 74 are processed to derive a graph referred to herein as the Periodic Information Plot (PIP) (step 76). Several methods for processing the vibration spectrum 64 and the autocorrelation spectrum 74 may be used according to the first embodiment, three of which are described below.
Because the vibration spectrum is twice the resolution of the autocorrelation spectrum, a point-to-point comparison for values on the x-axis (frequency) between the two spectra is not possible. However, a point-to-point comparison can be made by mathematically combining the amplitude values of two x-axis values in the vibration spectrum (step 65) for each associated x-axis value in the autocorrelation spectrum. Each XAS(n) value of the autocorrelation spectrum (where n=1 . . . N, and N is the number of lines of resolution for the autocorrelation spectrum) is mapped to the XVS(2n) value on the vibration spectrum. The mathematically combined x-axis value is defined such that XMCVS(n)=XVS(2n). The mathematically combined amplitude values YVS(2n) and YVS(2n−1) (herein termed YMCVS(n)) associated with the XMCVS(n) value from the vibration spectrum are calculated from the amplitudes of both the XVS(2n) and XVS(2n−1) frequencies from the x-axis. The calculation for deriving the mathematically combined amplitude value associated with the XMCVS(n) value from the vibration spectrum is:
Y
MCVS(n)=√{square root over ((YVS(2n−1))2+(YVS(2n))2)} Eq. (0)
where n=1 . . . N and N is the number of lines of resolution found in the autocorrelation spectrum.
In a first method (step 76a), for each X-value in the PIP (XPIP1), the Y-value in the PIP (YPIP1) is determined by multiplying the mathematically combined Y-value in the vibration spectrum (YMCVS) by the corresponding Y-value in the autocorrelation spectrum (YAS), according to:
Y
PIP1(n)=YMCVS(n)×YAS(n) Eq. (1)
for n=1 to N, where N is the number of X-values (frequency values) in the autocorrelation spectrum. Since amplitudes of periodic signals in the autocorrelation spectrum are higher than the amplitudes of random signals, the multiplication process will accentuate the periodic peaks while decreasing non-periodic peaks. An example of a PIP formed by the first method is depicted in
In a second method (step 76b), for each X-value in the PIP (XPIP2), the Y-value in the PIP (YPIP2) is determined by comparing the corresponding Y-value in the autocorrelation spectrum (YAS) to a predetermined threshold value (YTHR). For each autocorrelation spectrum amplitude greater than this threshold value, the associated amplitude for PIP (YPIP2(n)) will be set to the corresponding mathematically combined value from the vibration spectrum (YMCVS(n)). YAS values above the predetermined threshold indicate data that is largely periodic. Thus, the YPIP2 values are determined according to:
If YAS(n)>YTHR,YPIP2(n)=YMCVS(n) Eq. (2a)
If YAS(n)≤YTHR,YPIP2(n)=0 (or some other default level) Eq. (2b)
for n=1 to N.
In one preferred embodiment of the second method, YTHR is set to only include a percentage of the largest peaks from the autocorrelation spectrum. The percentage may be calculated based on the percent periodic signal in the autocorrelation waveform. The percent periodic signal is calculated based on the autocorrelation coefficient, which is the square root of the Y-value of the largest peak in the autocorrelation waveform. For this method, only the percent periodic signal of the total number of autocorrelation spectrum peaks will be evaluated. An example of a PIP formed by this method, with YTHR set to 59%, is depicted in
In another preferred embodiment of the second method, YTHR is set to include only peaks with values that are within the “percent periodic signal” of the largest peak value in the autocorrelation spectrum. These peaks, along with their harmonics that appear in the autocorrelation spectrum, will be utilized as the group of peaks to be intersected with those in the vibration spectrum to form the PIP. An example of a PIP formed by this method, with YTHR set to 59%, is depicted in
In a third method (step 76c), the PIP is determined according to the first method described above, and then the threshold of the second method is applied to the PIP according to:
If YPIP1(n)>YTHR,YPIP3(n)=YPIP1(n) Eq. (3a)
If YPIP1(n)≤YTHR,YPIP3(n)=0 (or some other default level) Eq. (3b)
for n=1 to N. An example of a PIP formed by this method is depicted in
Some embodiments also derive a Non-periodic Information Plot (NPIP) that consists of only the Y-values of the autocorrelation spectrum that are less than a predetermined threshold (step 78). Thus, the NPIP includes only non-periodic components. An example of an NPIP formed by this method is depicted in
Some embodiments also derive a Periodicity Map from the vibration spectrum and the autocorrelation spectrum (step 82). The Periodicity Map is created by pairing the mathematically combined Y-values from the vibration spectrum and the autocorrelation spectrum corresponding to any given X-value of the autocorrelation spectrum. These pairs are plotted with the mathematically combined Y-value from the vibration spectrum YMCVS(n) as the X-value of the point on the map XPM(n), and the Y-value from the autocorrelation spectrum YAS(n) as the corresponding Y-value on the map YPM(n), according to:
X
PM(n)=YMCVS(n) Eq. (4a)
Y
PM(n)=YAS(n) Eq. (4b)
for n=1 to N. As shown in
Some embodiments also derive a Circular Information Plot from any of the Periodic Information Plots described above (step 80). Once a linear PIP is calculated, an inverse FFT can be applied to generate an “information waveform.” A Circular Information Plot can then be generated from this information waveform. An example of a Circular Information Plot formed by this method is depicted in
Although preferred embodiments of the invention operate on vibration signals, the invention is not limited to only vibration signals. Periodic Signal Parameters and Periodic Information Plots may be derived from any signal containing periodic components.
In a second embodiment, a signal is collected from plant equipment (i.e. rotating or reciprocating equipment) and is processed using the method 300 depicted in
First, a waveform is generated (step 302 of
The waveform from step 302 is autocorrelated (step 314) to generate an autocorrelation waveform 316, having time on the X-axis and the correlation factor on the Y-axis. An FFT of the autocorrelation waveform 316 is calculated using the same Fmax as was used in the calculation of the FFT of the vibration waveform (step 318), resulting in an autocorrelation spectrum 320. Using the same Fmax forces the lines of resolution (LOR) of the autocorrelation spectrum 320 to be half of the LOR used in calculating the vibration spectrum 306. Since random events have largely been removed from the autocorrelation waveform 316, the remaining signal in the autocorrelation spectrum 320 is strongly related to periodic events. As shown in
Percent Periodic Energy (% Periodic Energy) is the percentage of energy in the vibration spectrum 306 that is related to periodic signals. It is calculated at step 322 based on the autocorrelation waveform 316 according to:
% Periodic Energy=√{square root over (MaxPeak (after 3% of autocorrelation waveform))}.
In a preferred embodiment, the total energy of the vibration spectrum 306 is calculated as the square root of the sum of the squares of each bin value in the vibration spectrum 306 ranging from zero to Fmax. For purposes of finding bearing and/or gear teeth faults, the vibration spectrum 306 is the PeakVue spectrum.
The percent energy of the vibration spectrum 306 is calculated at step 308 according to:
% Energy of Original=Total energy of vibration spectrum×% Periodic Energy.
A list of peaks from the vibration spectrum 306 is generated, wherein each listed peak is a located peak having a located frequency and an associated located amplitude (step 310). A list of peaks from the autocorrelation spectrum 320 is also generated, wherein each listed peak is a located peak having a located frequency and an associated located amplitude (step 324). In both lists, the peaks are arranged in order of descending amplitude, such that the peak having the largest amplitude is first in the list and the peak having the smallest amplitude is last (steps 312 and 326).
For the frequency value of each peak in the peak list generated for the autocorrelation spectrum, an associated matching peak is found in the peak list generated for the vibration spectrum (step 328). For a peak to “match,” the frequency value of the peak from the vibration spectrum 306 must be within N×ΔFrequency of the frequency value of the peak from the autocorrelation spectrum 320, where in a preferred embodiment N=4 and ΔFrequency is expressed as:
Thus, a match exists when
|original peak frequency−autocorrelation peak frequency|≤N×ΔFrequency.
For each matching peak from the vibration spectrum 306 found in step 328, the values of the located frequency and located amplitude is added to a PIP peak list (step 330). As each matching peak is added to the PIP peak list, a running Total Peak Energy value of all peaks in the PIP peak list is calculated (step 332). Because a Hanning window is used in the FFT calculation for this embodiment, the energy of a located peak is the result of energy from three bin values used in the creation of the located peak.
For each Total Peak Energy≤% Energy of Original, discard the associated peak in step 330 from the Autocorrelation Spectrum peak list before returning to step 328 (step 335). This process of matching peaks and adding matched peaks to the PIP peak list continues until
Total Peak Energy>% Energy of Original (step 334).
The Periodic Information Plot (PIP) is created by plotting the three points associated with each peak in the PIP peak list (step 336). In the preferred embodiment, the three points correspond to three bins associated with each located peak, assuming a Hanning window is used for FFT calculations. Examples of PIP's created using the method 300 of
Periodic Peaks
Periodic peaks in a spectrum are classified as either synchronous or non-synchronous peaks. Synchronous peaks are peaks that occur at the running speed of a shaft and its harmonic frequencies. For a gearbox having multiple shafts, there are also multiple families of synchronous peaks, wherein each family is associated with the speed of a particular shaft in the gearbox. In addition to running speed peaks, synchronous peaks associated with a gearbox also occur at all hunting tooth fundamental frequencies and their harmonics. Non-synchronous peaks are periodic families of harmonic peaks that are not members of a synchronous family. A family of non-synchronous, periodic peaks is most likely related to a bearing defect.
Because there may be many families of peaks related to either synchronous or non-synchronous peaks, a preferred embodiment provides a display color scheme to separate the different families of peaks. By color coding the different families in a spectrum, it is easy to distinguish between frequencies related to bearings (non-synchronous) and those related to running speed. In a gearbox, analysis of these running speed harmonic families (synchronous) can lead to the discovery of gear teeth problems. Using colors to designate the different families of peaks in a spectrum display or in the Periodic Information Plot simplifies the analysis for both the novice and experienced analyst.
Methods for Sorting and Discarding Statistically Outlying Peaks in the Autocorrelation Waveform (Step 34 in
The following routine takes an array of data values, such as values of positive peaks in the autocorrelation waveform, and discards values outside the statistically calculated boundaries. In a preferred embodiment, there are four methods or criteria for setting the boundaries.
Method 1: Non-Conservative, Using Minimum and Maximum Statistical Boundaries
Consider an array of P values (or elements) where P0 represents the number of values in the present array under evaluation. Now let P−1 represent the number of values in the array evaluated a single step before P0, let P−2 represent the number of values in the array evaluated a single step before P−1, and let P−3 represent the number of values in the array evaluated a single step before P−2.
Step 1:
Step 2:
Step 3:
Method 2: Non-Conservative, Using Maximum Statistical Boundary Only (No Minimum Boundary)
Use the same procedure as in Method 1 except only values exceeding the upper statistical boundaries are discarded. The minimum boundary is set to zero.
Method 3: Conservative, Using Minimum and Maximum Statistical Boundaries
Discard values based on Method 1, Step 1 only.
Method 4: Conservative, Using Maximum Statistical Boundary Only (No Minimum Boundary)
Discard values based on Method 1, Step 1 only and based on values exceeding the upper statistical boundaries. The minimum boundary is set to zero.
As an example of the sorting Method 1, consider an original set of values, P0, containing the twenty-one values listed below in Table 3 below, with n=1.
The mean (μ) of this original set, P0, is 0.54955 and standard deviation (σ) is 0.13982. Therefore, in Step 1 of Method 1,
Since 0.25442 is greater than 0.1, calculate
μ−nσ=0.54955−1*0.13982=0.409735
and
μ+nσ=0.54955+1*0.13982=0.689373.
Next, define the set P−1=P0 and define a new set P0, the values of which are all the values of P−1 that are between the values μ+σ=0.689343 and μ−σ=0.409735. The set P0 now contains the values listed below in Table 4, wherein three outlier values have been eliminated.
Since P0≠P−1, Step 1 is repeated, where for the set P0:
μ=0.50234,
σ=0.06946,
σ/μ=0.138263,
μ+σ=0.571797, and
μ−σ=0.432887.
Now define the set P−2=P−1, and P−1=P0 and define a new set P0, the values of which are all the values of P−1 that are between the values μ+σ=0.571797 and μ−σ=0.432887. The set P0 now contains the values listed below in Table 5, wherein four more outlier values have been eliminated.
Since P0≠P−1, Step 1 is repeated, where for the set P0:
μ=0.481311,
σ=0.037568, and
σ/μ=0.078053.
Since
σ/μ=0.078053≤1,
all the members of the array P0 are statistically close in value and need no more sorting.
If at any point in the calculations P0=P−1 and P−1≠P−2, then Step 2 would be executed instead of Step 1. In the example above, since P0≠P−1 for every iteration, only Step 1 was necessary for the calculations.
Predicting Bearing Faults Based on Periodic Signal Parameter (PSP)
In a preferred embodiment, alert amplitude limit levels (in g's) are determined based on the nominal turning speed according to the relationship depicted in
Before calculations of severity values can be made, Percent Periodic Energy must be calculated. Percent Periodic Energy (step 414) is calculated from the autocorrelation waveform according to:
% Periodic Energy=√{square root over (MaxPeak(after first 3%))}
wherein the maximum peak in the autocorrelation waveform does not include the first 3% of the waveform. Generally, the Percent Periodic Energy calculation is not as accurate for values less than 50%. Accordingly, as indicated in
In a preferred embodiment, the severity value is normalized by multiplying the result of step 416 by a desired maximum gauge value x according to:
Normalized General Severity=General Severity×x (step 418).
For the gauges shown in
Normalized General Severity=General Severity×10.
If the PSP is greater than 0.1 (step 419), a bearing fault is possibly present. Bearing Fault Severity (BFS) may be calculated according to:
BFS=Normalized Severity×% Periodic Energy (step 430).
If the resulting answer is greater than x (10 in this example), then the answer is truncated to be x.
In some embodiments, knowledge of the turning speed improves confidence that the periodicity is related to bearing faults and not turning speed incidences. When the turning speed is known, periodic peaks from the periodic information plot (PIP) can be classified as synchronous and non-synchronous. If only synchronous peaks are present, no bearing fault is indicated. If significant non-synchronous peaks are present, a possible bearing issue is confirmed, as indicated by:
If PSP≤0.1 and MaxPeak is <alert level, no fault is indicated by the measurement, meaning the asset is in good condition.
If PSP is less than or equal to 0.1 and MaxPeak is greater than the alert amplitude limit level (step 420), a deficiency in bearing lubrication is indicated. In addition, there may be lubrication issues when a bearing fault is present. (This is shown in
As shown in
The Lubrication Severity (LS) value is determined according to:
where x is the normalization value (step 426). For the Lubrication Severity gauge shown in
In an alternative embodiment, instead of determining whether PSP is greater than 0.1 in step 114, it is determined whether % Periodic Energy is greater than Y, where in most cases Y is 50%.
While the preferred embodiment of the algorithm described above and depicted in
Following are four examples that demonstrate use of the algorithm of
Predicting Gearbox Faults Based on Periodic Signal Parameter (PSP)
The rotational speed of at least one of the shafts in the gearbox is measured, such as using a tachometer (step 212), and the speed of each of the other shafts in the gearbox is calculated based on the speed measured in step 212 and knowledge of the gear ratios for the other shafts (step 214). In addition, based on shaft running speeds, hunting tooth frequencies are calculated based on techniques known to those of ordinary skill in the art. In a preferred embodiment, alert amplitude limit levels (in g's) are determined based on the nominal turning speed according to the relationship depicted in
Before calculations of specific severity values can be made, Percent Periodic Energy must be calculated. In a preferred embodiment, Percent Periodic Energy is calculated from the autocorrelation waveform according to:
% Periodic Energy=√{square root over (MaxPeak (after first 3%))}
wherein the MaxPeak of the autocorrelation waveform does not include the first 3% of the waveform (step 218). Generally, the Percent Periodic Energy calculation is not as accurate for values less than 50%. Accordingly, as indicated in
In order to calculate severity values for different faults, a general severity value is determined. General Severity may be calculated according to:
The severity value is normalized by multiplying the result of step 220 by a desired maximum gauge value x according to:
Normalized General Severity=General Severity×x (step 222).
For the gauge shown in
Normalized General Severity=General Severity×10.
The PIP is generated using the procedure described herein with reference to
If the PSP is greater than 0.1 (step 225), periodic frequencies related to the gearbox and/or bearings are present.
Based on knowledge of the turning speed, periodic peaks from the periodic information plot (PIP) can be classified as synchronous and non-synchronous. If non-synchronous peaks are present in the PIP (step 226), a bearing fault severity (BFS) value may be calculated (step 228) and displayed (step 234) according to:
If synchronous peaks are present (step 230) and fault limits are exceeded, gear teeth degradation is indicated. A gearbox fault severity (GFS) value may be calculated (step 232) and displayed (step 234) according to:
If the resulting answer is greater than x (10 in this example), then the answer is truncated to be x.
If PSP≤0.1 and MaxPeak is <alert level, no fault is indicated by the measurement, meaning the asset is in good condition.
If PSP is less than or equal to 0.1 and MaxPeak is greater than the alert amplitude limit level (step 234), a deficiency in bearing and/or gearbox lubrication is indicated. In addition, there may be lubrication issues along with mechanical faults present. (This is shown in
As discussed above, Percent Non-periodic energy (% NPE) is a function of Percent Periodic Energy and can be determined using the plot of
The bearing or gearbox lubrication severity value is determined and displayed according to:
where x is the normalization value (steps 240 and 242). For the Lubrication Severity gauge shown in
In an alternative embodiment, instead of determining whether PSP is greater than 0.1 in step 218, it is determined whether % Periodic Energy is greater than Y, where in most cases Y is 50%.
In a third embodiment, a signal is collected from plant equipment (i.e. rotating or reciprocating equipment) and is processed using the method 500 depicted in
First, a waveform is generated (step 502 of
The waveform from step 502 is autocorrelated (step 514) to generate an autocorrelation waveform 516, having time on the X-axis and the correlation factor on the Y-axis. An FFT of the autocorrelation waveform 516 is calculated using the same Fmax as was used in the calculation of the FFT of the original waveform (step 518), resulting in an autocorrelation spectrum 320. Using the same Fmax forces the lines of resolution (LOR) of the autocorrelation spectrum 520 to be half of the LOR used in calculating the vibration spectrum 506. Since random events have largely been removed from the autocorrelation waveform 516, the remaining signal in the autocorrelation spectrum 520 is strongly related to periodic events. As shown in
Percent Periodic Energy (% Periodic Energy) is calculated at step 522 based on the autocorrelation waveform 516 according to:
% Periodic Energy=√{square root over (MaxPeak (after 3% of autocorrelation waveform))}.
In a preferred embodiment, the total energy of the vibration spectrum 506 is calculated as the square root of the sum of the squares of each bin value in the vibration spectrum 506 ranging from zero to Fmax (step 507). For purposes of finding bearing and/or gear teeth faults, the original spectrum 506 is the PeakVue spectrum.
The percent energy of the original spectrum 506 is calculated at step 508 according to:
% Energy of Original=Total energy of original spectrum×% Periodic Energy.
A list of the largest peaks from the vibration spectrum 506 is generated, wherein each listed peak is a located peak having a located frequency and an associated located amplitude (step 510). A list of the largest peaks from the autocorrelation spectrum 520 is also generated, wherein each listed peak is a located peak having a located frequency and an associated located amplitude (step 524). In both lists, there may be N number of peaks (such as N=300) that are preferably arranged in order of descending amplitude, such that the peak having the largest amplitude is first in the list and the peak having the smallest amplitude is last (steps 512 and 526).
Further details of step 526 are depicted in
For the frequency value of each of the K number of peaks in the periodic peak list, an associated matching peak is found in the peak list generated for the vibration spectrum 506 (step 528). For a peak to “match,” the frequency value of the peak from the vibration spectrum 506 must be within N×ΔFrequency of the frequency value of the peak from the periodic peak list, where in a preferred embodiment N=3 and ΔFrequency is expressed as:
Thus, a match exists when
|original peak frequency−periodic peak list peak frequency|≤N×ΔFrequency.
The total number of matches may be designated as M.
For each matching peak from the vibration spectrum 506 found in step 528, the values of the located frequency and located amplitude is added to the periodic peak list (step 530). As each matching peak is added to the periodic peak list, a running Total Peak Energy value of all peaks in the periodic peak list is calculated (step 532). Because a Hanning window is used in the FFT calculation for this embodiment, the energy of a located peak is the result of energy from three bin values used in the creation of the located peak.
For each Total Peak Energy value that is less than or equal to the % Energy of the original vibration spectrum (calculated at step 508), the associated largest peak (from step 350) is discarded from the periodic peak list before returning to step 528 (step 535). This process of matching peaks and adding matched peaks to the periodic peak list continues until the number of iterations equal K, or the Total Peak Energy is greater than or equal to the % Energy of the original vibration spectrum, whichever event occurs first (step 534).
The Periodic Information Plot (PIP) is created by plotting the three points associated with each peak in the periodic peak list (step 536). The three points correspond to three bins associated with each located peak, assuming a Hanning window is used for FFT calculations. In a preferred embodiment, P number of PIP peaks are plotted starting with the first peak in the periodic peak list, wherein the number P is the maximum desired number of peaks to be plotted, and P≤M. For example, in the case of a CSI 9420 Wireless Vibration Transmitter, only twenty peaks (P=20) may be displayed for a “thumbnail” plot. An example of a PIP created using the method 500 of
Periodic Peaks
The periodic peaks in a spectrum are classified as either synchronous or non-synchronous peaks, as defined below.
Synchronous periodic peaks are all running speed peaks and their harmonics. For a gearbox, there are several families of synchronous periodic peaks—one family associated with the speed of each shaft in a gearbox. In addition to running speed peaks, synchronous periodic peaks are also defined for a gearbox as all hunting tooth fundamental frequencies and their harmonics.
Non-synchronous periodic peaks are any periodic family of harmonic peaks that are not members of a synchronous family. Most non-synchronous periodic peaks are related to bearing frequencies.
It is very important to obtain the correct turning speed (and associated turning speeds of other shafts for gearboxes) of the machine being evaluated. As for any analysis, incorrect speed will produce incorrect locations of the synchronous and non-synchronous periodic peaks. Precise speeds are typically obtained from tachometer/strobe outputs or manual analysis of spectra. If the precise speed is unknown, a Speed Detection Algorithm should be employed.
Because there may be many families of peaks related to either synchronous or non-synchronous periodic peaks, a color scheme may be devised to separate the different family of peaks. By color coding the different families of peaks, it is easy to distinguish frequencies related to bearings (non-synchronous) from those related to running speed. In a gearbox, these running speed harmonic families (synchronous) can lead to the discovery of gear teeth problems. Using colors to designate the different families of periodic peaks in the Periodic Information Plot simplifies analysis for both the novice and experienced analyst.
The use of colors to designate the different families of periodic peaks is depicted in
An alternative to using color is to display each family of peaks with a different line thickness and/or pattern, an example of which is depicted in
Although preferred embodiments described herein are applied to PeakVue™ vibration data, the techniques can be applied to any analysis in which random (noise) data is to be separated from desired informative data. Examples include analysis of acceleration and/or velocity vibration data.
Embodiments described herein also have application in analysis of vibration of pumps due to cavitation. Cavitation in pumps causes vibration signals that are random in nature and therefore appear as noise in an FFT spectrum. This randomness implies that vibration signals caused by cavitation are not periodic in nature, and therefore can be separated from the periodic signals using the processes described herein.
Cavitation appears in the FFT as noise (random signals) about the blade pass frequency of a pump, which is the number of fan blades in the pump multiplied by the running speed. Cavitation may also show up as high frequencies that are detectable as non-periodic signals using PeakVue™ waveform analysis. The basic periodicity algorithms described in U.S. Pat. No. 9,791,422 and US Patent Publication No. 2018/0011065A1 may be used to isolate these indicators of cavitation. Autocorrelation can be applied to both velocity, acceleration, and PeakVue™ waveforms to isolate the randomness generated by cavitation.
Generally, the detection of cavitation differs from the detection of bearing or gear defects in that auto-correlation must be applied to velocity waveforms and the “non-periodic” indicator must be trended to see changes in the non-periodicity.
Generating List of Vibration Peaks Associated with Gear Mesh Frequency
One preferred embodiment provides a method to separate desired periodic vibration peaks, such as peaks associated with the gear mesh frequency (GMF), from noise within a certain frequency region of interest in a vibration spectrum. Generally, GMF peaks have amplitudes that are much higher than the noise floor. However, there are times when it may be difficult to easily separate such peaks from the noise using known methods, such as by simply picking the 300 largest peaks in a spectrum. Described herein with reference to
The GMF is calculated as
GMF=n×StartFreq,
where n is the number of teeth of the gear on the shaft associated with the GMF, and StartFreq is the rotational speed of the shaft (step 602). The ROI is calculated using the GMF as the FOI (step 604). In the example of
FreqDiff=PctRange×Fmax/100=5×Fmax/100,
wherein Fmax is the maximum frequency of the vibration spectrum. If the FOI is greater than Fmax, then the FOI is set to be equal to Fmax. The minimum frequency in the ROI is determined as
MinRangeFreq=FOI−FreqDiff,
and the maximum frequency in the ROI is determined as
MaxRangeFreq=FOI+FreqDiff.
If MinRangeFreq is less than or equal to zero, then MinRangeFreq is set equal to zero, and if MaxRangeFreq is greater than or equal to Fmax, then MaxRangeFreq is set equal to Fmax.
In step 606, a located peak subset comprising the M number of highest amplitude peaks having frequency values between MinRangeFreq and MaxRangeFreq is created. This subset of highest amplitude peaks, referred to herein as GB_LocPeak(m), contains an amplitude value (GB_AmplLocPeak(m)) and a frequency value (GB_FreqLocPeak(m)) for each of the located peaks, where m=1 to M. In a preferred embodiment, M is 300.
In step 608, the located peaks in the GB_LocPeak(m) list are then sorted from highest to lowest amplitude to create a sorted peak list, with GB_LocPeak(1) associated with the highest amplitude peak, GB_LocPeak(2) associated with the next highest amplitude peak, and GB_LocPeak(M) associated with the lowest amplitude peak.
In step 610, the Kurtosis value for the M number of peaks in the sorted list is calculated. If the Kurtosis value is greater than 3, the highest amplitude peak (GB_LocPeak(1))) is removed from the sorted list, and its amplitude and frequency values are stored in a largest peak list as GB_AmplPeak(1) and GB_FreqPeak(1). The Kurtosis value is again calculated for the remaining number (M−1) of peaks in the sorted list, and if the value is still greater than 3, the next highest amplitude peak (GB_LocPeak(2)) is removed from the sorted list, and its amplitude and frequency values are stored in the largest peak list as GB_AmplPeak(2) and GB_FreqPeak(2). This process of moving the highest amplitude peaks from the sorted list to the largest peak list continues until the Kurtosis value for the remaining peaks in the sorted list is less than or equal to 3.
In step 612, the J number of gear mesh candidate frequencies are determined to be the frequencies of the peaks in the largest peak list for which the following condition is met:
LowerFreqLimit<GB_FreqPeak(i)<UpperFreqLimit
wherein the frequencies and amplitudes of the candidate peaks are
FreqCandidate(j)=GB_FreqPeak(i)
AmplCandidate(j)=GB_AmplPeak(i).
In step 614, the candidate rotational speeds of the gear shaft are determined as:
SpeedCandidate(j)=FreqCandidate(j)/n.
In a preferred embodiment, the frequency limits—LowerFreqLimit and UpperFreqLimit—are based on the nameplate speed information, such as may be used in a speed detection algorithm. The lower limit is preferably the “Service Factor Speed” (Synchronous Speed*(1−1.2*% slip) which is 120% of full load speed. The upper limit is preferably the synchronous speed.
The process of
The foregoing description of preferred embodiments for this invention has been presented for purposes of illustration and description. They are not intended to be exhaustive or to limit the invention to the precise form disclosed. Obvious modifications or variations are possible in light of the above teachings. The embodiments are chosen and described in an effort to provide the best illustrations of the principles of the invention and its practical application, and to thereby enable one of ordinary skill in the art to utilize the invention in various embodiments and with various modifications as are suited to the particular use contemplated. All such modifications and variations are within the scope of the invention as determined by the appended claims when interpreted in accordance with the breadth to which they are fairly, legally, and equitably entitled.
This application claims priority as a continuation-in-part of U.S. non-provisional patent application Ser. No. 16/515,793 filed Jul. 18, 2019, titled “System for Separating Periodic Amplitude Peaks from Non-periodic Amplitude Peaks in Machine Vibration Data,” which claims priority as a continuation-in-part of U.S. non-provisional patent application Ser. No. 15/697,911 filed Sep. 7, 2017, titled “Machine Fault Prediction Based on Analysis of Periodic Information in a Signal,” which claims priority as a continuation-in-part of U.S. non-provisional patent application Ser. No. 14/316,883 filed Jun. 27, 2014, titled “Analysis of Periodic Information in a Signal,” which claims priority to U.S. provisional patent application No. 61/842,035 filed Jul. 2, 2013, titled “Periodic Signal Parameter.” All of the subject matter described in the above referenced applications is incorporated herein by reference.
Number | Date | Country | |
---|---|---|---|
62813863 | Mar 2019 | US | |
62410935 | Oct 2016 | US | |
61842035 | Jul 2013 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 16515793 | Jul 2019 | US |
Child | 17227271 | US | |
Parent | 15697911 | Sep 2017 | US |
Child | 16515793 | US | |
Parent | 14316883 | Jun 2014 | US |
Child | 15697911 | US |