CN110617964A - Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing - Google Patents
Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing Download PDFInfo
- Publication number
- CN110617964A CN110617964A CN201910687690.5A CN201910687690A CN110617964A CN 110617964 A CN110617964 A CN 110617964A CN 201910687690 A CN201910687690 A CN 201910687690A CN 110617964 A CN110617964 A CN 110617964A
- Authority
- CN
- China
- Prior art keywords
- frequency
- instantaneous frequency
- order
- ratio analysis
- analysis method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
Abstract
The invention discloses a synchronous compression transformation order ratio analysis method for fault diagnosis of a rolling bearing, which comprises the steps of firstly analyzing variable-speed vibration signals of the bearing by utilizing SST to obtain time-frequency distribution of the variable-speed vibration signals, then extracting frequency-converted instantaneous frequency ridge lines and carrying out high-order polynomial fitting on the frequency-converted instantaneous frequency ridge lines to obtain a fitted instantaneous frequency curve, solving a phase discrimination time scale according to the fitted instantaneous frequency curve, then carrying out equal-angle resampling on original vibration signals to obtain angle domain signals, then carrying out Hilbert envelope demodulation on the angle domain signals to obtain an envelope order spectrum of the angle domain signals, judging whether the rolling bearing has faults and fault types of the rolling bearing by analyzing the envelope order spectrum, giving diagnosis results and providing maintenance suggestions, thereby ensuring normal operation of subway vehicles. The synchronous compression transformation order ratio analysis method can accurately extract instantaneous frequency conversion information contained in the vibration signal, realizes order ratio tracking, does not need a specific hardware device, greatly simplifies the installation process and reduces the fault monitoring cost.
Description
Technical Field
The invention relates to the technical field of metro vehicles, in particular to a synchronous compression transformation order ratio analysis method for fault diagnosis of a rolling bearing.
Background
The bearing is used as an important part of the subway axle box and plays an important role in the stable operation of the subway. The vibration signal of the rolling bearing during variable-speed operation contains more abundant characteristic information than constant-speed. In recent years, order ratio analysis is increasingly used in variable-speed mechanical fault diagnosis, and the basic principle is to convert an unstable time domain signal into a quasi-stable angular domain signal through equiangular resampling.
At present, the scale analysis technology at home and abroad is divided into three types: hardware-based step ratio tracking using phase detection means, computational step ratio tracking (COT) and estimation of rotating mechanical step ratio tracking based on instantaneous frequency.
Hardware step ratio tracking requires constant angle incremental sampling of vibration signals by phase detection equipment, but since these equipment are not only complex in result and inconvenient to install, but also expensive, certain production cost is increased. The step ratio tracking is calculated by independently acquiring rotating speed signals by using a tachometer, obtaining a phase discrimination time scale through a corresponding interpolation algorithm and realizing equal-angle interval sampling on the basis of software, but for some equipment, the tachometer is inconvenient to install, and when the rotating speed of rotating equipment changes too fast, the obtained rotating speed change information is inaccurate, so that although the step ratio tracking is more convenient than hardware type step ratio tracking, a specific device is not required, and the limitation of the tachometer can not be got rid of. The method for tracking the order ratio of the rotating machinery based on instantaneous frequency estimation comprises the following key steps of extracting instantaneous frequency, wherein the commonly used methods for estimating the instantaneous frequency comprise the following steps: short-time Fourier transform (STFT), Wigner Distribution (WD), Wavelet Transform (WT), and the like. The time-frequency resolution of the STFT is not high, and is susceptible to noise interference, which may affect the extraction of the instantaneous frequency. WD has difficulty extracting instantaneous frequency curves due to the presence of cross terms. WT cannot improve the aggregation of time and frequency, and is easily interfered, and the extracted instantaneous frequency curve has low precision.
Disclosure of Invention
In order to solve the above problems, an object of the present invention is to provide a synchronous compression transformation order ratio analysis method for diagnosing a fault of a rolling bearing, which can effectively extract a frequency conversion in a vibration signal of the bearing, realize order ratio analysis based on instantaneous frequency estimation, and provide a reliable theoretical basis for diagnosing a fault of a metro vehicle during a speed change process.
The invention provides a synchronous compression transformation order ratio analysis method for fault diagnosis of a rolling bearing, which comprises the following steps:
step 1, carrying out low-pass filtering and down-sampling treatment on the collected vibration signals of the variable rotating speed of the bearing;
step 2, SST time frequency analysis processing is carried out on the processed vibration signal to obtain the time frequency spectrum thereof, the instantaneous frequency in the vibration signal is estimated, and the instantaneous frequency ridge line f of the frequency conversion is extractedcs(k) Namely the instantaneous frequency of the rotating speed, wherein k is the serial number of a sampling point;
step 3, extracting the ridge line f of the instantaneous frequencycs(k) Performing high-order polynomial fitting to obtain a fitted instantaneous frequency curve fi(t);
Step 4, utilizing the obtained fitting instantaneous frequency curve fi(T) carrying out interpolation filtering, and calculating to obtain a phase discrimination time scale sequence Tn;
Step 5, according to the phase discrimination time scale sequence TnCarrying out equal-angle resampling on the original vibration signal to obtain a quasi-steady-state angular domain signal sequence x (T)n);
Step 6, diagonal domain signal sequence x (T)n) Further carrying out Hilbert envelope demodulation, and solving an envelope order spectrum through FFT (fast Fourier transform), so as to realize order ratio analysis on the vibration signal;
and 7, judging whether the rolling bearing has faults and fault types thereof by analyzing the envelope order spectrum, giving a diagnosis result, proposing a maintenance suggestion, and completing the fault diagnosis of the bearing in the subway speed change process.
As a further improvement of the present invention, in step 1, the cut-off frequency of the low-pass filtering cannot be lower than the maximum value of the conversion frequency.
As a further improvement of the invention, in the step 2, the peak search is used on the time frequency spectrum to obtain the ridge line f of the instantaneous frequency of the frequency conversioncs(k)。
As a further improvement of the invention, in the step 3, the ridge line f of the instantaneous frequencycs(k) When fitting is carried out, high-order polynomial fitting and spline fitting are adopted in a small range.
As a further improvement of the present invention, in step 3, a second-order polynomial fitting is adopted, and a fitting equation is:
fi(t)=at2+bt+c (1)
assuming ε is the error of the fitting equation, the sum of the squares of ε is:
if so:andfitting coefficients a, b and c can be obtained.
As a further improvement of the invention, in the step 4, the phase discrimination time scale sequence T is calculatednThe formula used is as follows:
in the formula, N is 1,2,3, N is a sampling time sequence number, and N is a sampling sequence length; delta theta pi/Omax, OmaxIs the maximum theoretical order;
integrating equation (3) yields:
get T0When the value is equal to 0, goOne-step calculation is carried out to obtain a phase discrimination time scale sequence T required by equal-angle resamplingn。
As a further improvement of the invention, in the step 5, a phase discrimination time scale sequence T is utilizednAccording to the interpolation principle, original data are subjected to equal-angle resampling, and unsteady-state time domain signals are converted into quasi-steady-state angular domain signal sequences x (T)n):
In the formula,. DELTA.tsIs a time domain sampling interval; m Δ tsIs at TnA value in the vicinity; h issAnd (t) is an interpolation filter.
The invention has the beneficial effects that:
the SST order ratio analysis method has good capability centralization, can concentrate the energy in time-frequency distribution on a main frequency curve when performing time-frequency analysis, can accurately extract instantaneous frequency conversion information contained in a vibration signal, and realizes order ratio tracking;
the SST order ratio analysis method of the invention does not need a specific hardware device, breaks the limitation of hardware, greatly simplifies the installation process, reduces the fault monitoring cost, can realize order ratio analysis of the variable-speed vibration signal of the bearing only by analyzing the vibration signal, and completes the fault diagnosis of the bearing.
Drawings
Fig. 1 is a schematic flow chart of a synchronous compression transformation order ratio analysis method for fault diagnosis of a rolling bearing according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a waveform of a constructed simulated acceleration signal and its frequency spectrum;
FIG. 3 is a schematic diagram of an ideal instantaneous frequency ridge compared to an instantaneous frequency ridge extracted using SST;
FIG. 4 is a schematic diagram comparing a fitted instantaneous frequency curve to an ideal instantaneous frequency curve fitted to the instantaneous frequency ridge extracted from SST of FIG. 3;
FIG. 5 is a graphical representation of the results of SST order ratio analysis of the fitted instantaneous frequency curve of FIG. 4;
FIG. 6 is a schematic diagram of an SQI rotor test bed used in the experiment;
FIG. 7 is a schematic view of the arrangement of the measuring points of the SQI rotor test bed;
FIG. 8 is a schematic diagram of a fault vibration signal of an ER-12K inner ring of a bearing collected when an SQI rotor test bed is started;
FIG. 9 is a schematic diagram illustrating the STFT and SST comparing the instantaneous frequency extraction of the filtered and down-sampled inner ring fault vibration signal;
FIG. 10 is a fitted instantaneous frequency curve obtained by second order polynomial fitting of the instantaneous frequency ridge extracted from SST of FIG. 9;
FIG. 11 is a schematic diagram of the angular domain signal obtained after SST order ratio analysis of the fitted instantaneous frequency curve of FIG. 10;
fig. 12 is a diagram illustrating an envelope of the angular domain signal and an envelope order spectrum thereof obtained by performing Hilbert envelope demodulation on the angular domain signal of fig. 11;
FIG. 13 is a schematic diagram of a fault vibration signal of an ER-12K outer ring of a bearing collected when an SQI rotor test bed is started;
FIG. 14 is a schematic diagram showing the comparison between the instantaneous frequency ridge line extracted by SST of the outer ring fault vibration signal and the actually calculated rotation speed;
FIG. 15 is a fitted instantaneous frequency curve obtained by second order polynomial fitting of the instantaneous frequency ridge extracted from SST of FIG. 14;
FIG. 16 is a schematic diagram of the angular domain signal obtained after SST order ratio analysis of the fitted instantaneous frequency curve of FIG. 15;
fig. 17 is a diagram illustrating an envelope of the angular domain signal and an envelope order spectrum thereof obtained by performing Hilbert envelope demodulation on the angular domain signal of fig. 16.
Detailed Description
The present invention will be described in further detail below with reference to specific embodiments and with reference to the attached drawings.
The synchronous compression transformation order ratio analysis method for fault diagnosis of the rolling bearing comprises the steps of firstly analyzing variable-rotating-speed vibration signals of the bearing by utilizing SST to obtain time-frequency distribution of the variable-rotating-speed vibration signals, then extracting instantaneous frequency ridges of frequency conversion, carrying out high-order polynomial fitting on the instantaneous frequency ridges to obtain a fitted instantaneous frequency curve, obtaining a phase discrimination time scale according to the fitted instantaneous frequency curve, carrying out equal-angle resampling on original vibration signals to obtain angle domain signals, further carrying out Hilbert envelope demodulation on the angle domain signals to obtain an envelope order spectrum of the angle domain signals, judging whether the rolling bearing has faults and fault types of the rolling bearing by analyzing the envelope order spectrum, giving diagnosis results, and giving maintenance suggestions, so that normal operation of subway vehicles is guaranteed.
As shown in fig. 1, the synchronous compression transformation (SST) order ratio analysis method of the present invention specifically includes the following steps:
step 1, low-pass filtering and down-sampling processing are carried out on the collected vibration signals of the variable rotating speed of the bearing.
Due to the interference of high-frequency noise, the rotating speed information is not easy to extract, so that the signal is firstly subjected to low-pass filtering and down-sampling processing. After low-pass filtering processing, high-frequency components in the vibration signals are removed, energy is mainly concentrated in the low-frequency components, the down-sampling also has a noise reduction effect, the analyzed data volume is reduced, and time is saved for extracting the instantaneous frequency ridge line.
Further, in order to ensure the integrity of the speed curve in this interval, the low-pass filter cut-off frequency cannot be lower than the maximum value of the rotational frequency.
Step 2, SST time frequency analysis processing is carried out on the processed vibration signal to obtain a time frequency spectrum thereof, peak value searching is utilized on the time frequency spectrum to estimate the instantaneous frequency in the vibration signal, and the instantaneous frequency ridge line f of frequency conversion is extractedcs(k) I.e. the instantaneous frequency of the rotation speed, where k is the sample point number.
The order ratio analysis method is based on a synchronous compression transform (SST) method, is based on wavelet transform, utilizes a synchronous compression operator to improve the resolution of a time-frequency ridge line on a time-frequency spectrum, and realizes extraction and reconstruction of instantaneous frequency. With ψ (b) as the mother wavelet function, the continuous wavelet transform of signal x (t) is:
wherein a is a scale factor, b is a translation factor,is the conjugate of ψ (t).
Assuming that the input signal is a single-frequency signal, x (t) ═ a · cos (ab), its fourier transform is x (f), and the fourier transform defining its wavelet function is Ψ (ξ), the wavelet transform and its reduction are:
the derivation of the wavelet transform in the translation direction can result in:
the instantaneous frequency information for the (a, b) position in the wavelet domain is therefore:
on the scale a, W (a, b), regardless of the value, the oscillation characteristic on b points to the initial frequency Ω, so:
according to the defined synchronous compression transform, the wavelet inverse transform:
coefficient of phase difference
Omega containing instantaneous frequency information for (a, b) locations in the wavelet domainx(a, b) are integrated along the direction of the scale a, and are classified into the frequency domain, wherein omega is omegaxAt the position of (a, b), the synchronous compression transform is defined as:
wherein A (b) { a, Wx(a,b)≠0}。
By the formula, the amplitude of the signal is classified to the position of the time frequency domain, and finally, the time frequency spectrum with high resolution can be obtained.
Step 3, extracting the ridge line f of the instantaneous frequencycs(k) Performing high-order polynomial fitting to obtain a fitted instantaneous frequency curve fi(t)。
Further, a method of high-order polynomial fitting and spline fitting may be performed in a small range to achieve more accurate approximation.
A second order polynomial fit is used here. Fitting by using a second-order polynomial, wherein the fitting equation is as follows:
fi(t)=at2+bt+c (1)
assuming ε is the error of the fitting equation, the sum of the squares of ε is:
if so:andfitting coefficients a, b and c can be obtained.
Step 4, utilizing the obtained fitting instantaneous frequency curve fi(T) carrying out interpolation filtering, and calculating to obtain a phase discrimination time scale sequence Tn。
In the calculation of phase discrimination time scale sequence TnThe formula used is as follows:
in the formula, N is 1,2,3, N is a sampling time sequence number, and N is a sampling sequence length; delta theta pi/Omax, OmaxIs the maximum theoretical order;
integrating equation (3) yields:
taking T in general0And (5) further calculating to obtain a phase discrimination time scale sequence T required by equal-angle resamplingn。
Step 5, utilizing the phase discrimination time mark sequence TnAccording to the interpolation principle, original data are subjected to equal-angle resampling, and unsteady-state time domain signals are converted into quasi-steady-state angular domain signal sequences x (T)n):
In the formula,. DELTA.tsIs a time domain sampling interval; m Δ tsIs at TnA value in the vicinity; h issAnd (t) is an interpolation filter.
Step 6, diagonal domain signal sequence x (T)n) And further carrying out Hilbert envelope demodulation, and solving an envelope order spectrum through FFT (fast Fourier transform), thereby realizing order ratio analysis on the vibration signal.
And 7, judging whether the rolling bearing has faults and fault types thereof by analyzing the envelope order spectrum, giving a diagnosis result, proposing a maintenance suggestion, and completing the fault diagnosis of the bearing in the subway speed change process.
The simulation verification of the SST order ratio analysis method proposed by the present invention will be described in detail below.
Since the subway vehicle is accelerated to a certain set speed at a constant acceleration after being out of the station, the constant speed operation is maintained, and the acceleration process is usually linear change, an acceleration signal with the linear change of the speed is constructed during simulation verification, and the following formula is shown below:
during simulation: the sampling frequency 2000HZ and the number of sampling points 2000 are set, and white gaussian noise with SNR-8 (signal to noise ratio) is added. As can be seen from the equation (6), there are three components in x (t), and assuming that the first component is the rotation speed variation component, the instantaneous frequency is f0200t, the instantaneous frequencies of the second and third components are f1=2.5f0、 f2=3.5f0. The order of the three components is 1, 2.5 and 3.5 respectively according to the relation between the order and the frequency conversion.
Fig. 2 is a waveform of a constructed simulated acceleration signal and a frequency spectrum thereof. As can be seen from fig. 2, the amplitude of the signal waveform gradually increases while becoming denser; the frequency in the frequency spectrum has a gradual change phenomenon, and fault characteristic information is difficult to find.
Fig. 3 is a graph showing the comparison of an ideal instantaneous frequency ridge with an instantaneous frequency ridge extracted using SST. As can be seen from fig. 3, at the beginning, due to the fact that x (t) is relatively small in amplitude and is buried by noise, the time-frequency distribution of the instantaneous frequency obtained by SST is not clear, but can reflect the linear variation trend of the frequency. With the increase of the amplitude, the time-frequency ridge line obtained by the SST is clearer and has strong concentration. The instantaneous frequency ridge of the SST extraction speed substantially coincides with the ideal instantaneous frequency ridge.
And fitting the instantaneous frequency ridge extracted by the SST in the graph 3 to obtain a fitted instantaneous frequency curve. FIG. 4 is a graph showing a comparison of a fitted instantaneous frequency curve with an ideal instantaneous frequency curve. As can be seen from fig. 4, the fitted instantaneous frequency curve almost coincides with the ideal instantaneous frequency curve, which shows that the SST can accurately extract the frequency conversion information in the vibration signal.
And performing interpolation filtering on the fitted instantaneous frequency curve in the graph 4 to obtain a phase discrimination time scale sequence, and performing equal-angle resampling on the original vibration signal according to the phase discrimination time scale sequence. The maximum order O needs to be set before equal-angle samplingmaxMaximum order follows Nyqu at 20ist theorem of sampling.
The order of the order analysis corresponds to the frequency, and for a rotating machine, the order refers to the relationship between all frequencies in the vibration signal and the ratio of the frequency to the rotating frequency, and the order of 1 indicates that the frequency is consistent with the rotating speed frequency.
Wherein f is frequency (HZ), l is order, and R is rotational speed (rpm).
Figure 5 is a graph showing the results of SST order ratio analysis performed on the fitted instantaneous frequency curve of figure 4. As can be seen from fig. 5, the angular domain signals after equal angle sampling have relatively uniform distribution on the time axis, and the frequency does not change with time. The order 1.002 in the envelope order spectrum of fig. 5 appears very close to the ideal order 1 and can therefore be considered as frequency-converted information, while the orders 2.516 and 3.518 correspond exactly to the second and third component components in the simulated signal. The SST order ratio analysis method provided by the invention can accurately and effectively convert the unsteady variable-speed signal into the steady angular domain signal and obtain the characteristic order, thereby realizing order ratio analysis.
Experimental verification of the SST scale ratio analysis proposed by the present invention will be detailed below.
As shown in FIG. 6, the SQI rotor test bench is used as an experimental platform during the experiment. As shown in fig. 7, the measuring points of the SQI rotor test stand are arranged as follows: the measuring point 1 is a tachometer, the measuring point 2 is a motor extending end bearing, the measuring point 3 is a bearing near the motor end, and the measuring point 4 is a bearing far away from the motor end. During the experiment, the position of the pre-buried fault bearing is located at the end close to a motor (a measuring point 3 in the drawing), and the model of the bearing is ER-12K. The failure characteristic frequency coefficients of the bearing are shown in table 1, for example: the inner ring fault characteristic frequency is 4.95 × revolution frequency. When the SQI rotor test bed is started, the collected fault vibration signal of the bearing inner ring is shown in figure 8. The characteristic order after the order ratio analysis is 4.95. The speed variation range of the rotor test bed is set to be 0-24.5 Hz, the sampling frequency is 6000Hz, and in order to obtain a good acceleration process, the length of the collected signal is 30000.
TABLE 1 ER-12K bearing failure coefficient
Because of the interference of high-frequency noise, the rotating speed information is not easy to extract, and therefore the inner ring fault vibration signal is firstly subjected to pass filtering and down-sampling treatment. It is verified that the extracted instantaneous frequency is the best when the cut-off frequency of the low-pass filtering is 40Hz and the down-sampling interval is 20 s. Fig. 9 is a comparison diagram of the filtered down-sampled instantaneous frequency extraction of the signal for STFT and SST. As can be seen from the time-frequency distribution diagram of fig. 9, compared with STFT, SST has good concentration, and the transient frequency ridge extracted by SST can better reflect the change details of the transient frequency.
Second order polynomial fitting is performed on the instantaneous frequency ridge extracted by SST, and the obtained fitted instantaneous frequency curve is shown in fig. 10. In fig. 10, the solid line is an actual rotating speed curve, the dotted line is a fitted instantaneous frequency curve, and after comparison, the fitted instantaneous frequency curve can basically reflect the variation trend of the rotating speed, and the fitting effect is good.
And solving a phase discrimination time scale through the fitted instantaneous frequency curve obtained in the step 10, and performing equal-angle resampling on the inner ring fault vibration signal in the step 8 by using an interpolation method. Fig. 11 is an angular domain signal obtained by performing SST order ratio analysis on the fitted instantaneous frequency curve of fig. 10, wherein the abscissa is the number of revolutions, and the ordinate is the amplitude, and the angular domain waveform is observed to find an obvious periodic impact phenomenon. FIG. 12 shows that Hilbert envelope demodulation is performed on the angular domain signal of FIG. 11 to obtain the envelope of the angular domain signal and the envelope order spectrum thereof, and it is found that the orders 1.017, 2.024 and 2.996 in the envelope order spectrum are quite obvious, and the signals are analyzed as the frequency conversion lr1-3 times of frequency. In addition, orders 4.979, 9.969, 14.95 are also more prominent than inner ring fault signature order liAnd the second and third harmonic order thereof are very closeWhile l is also foundr-liAnd lr+liThe orders are very close, and obvious bearing inner ring fault characteristics can be determined.
And collecting bearing outer ring fault vibration signals with the same length to continue experimental verification, wherein the collected bearing outer ring fault vibration signals are shown in fig. 13 when the SQI rotor test bed is started. As can be seen from fig. 13, the signal waveform diagram has a relatively obvious impact, but the impact does not appear periodically, so the SST processing is performed on the vibration signal in the acceleration process to extract the rotation speed change information.
Fig. 14 is a schematic diagram showing the comparison between the instantaneous frequency ridge line extracted by SST of the outer ring fault vibration signal and the actually calculated rotating speed. As shown in fig. 14, some fluctuation of the instantaneous frequency ridge of SST extraction was found due to noise interference, and the change of the frequency with time was smoother, but the trend of the frequency change could be reflected.
Fitting the instantaneous frequency ridge extracted from the SST in FIG. 14 to obtain a fitted instantaneous frequency curve as shown in FIG. 15, wherein the solid line is the actual rotating speed variation curve, the dotted line is the fitted instantaneous frequency curve, and the fitting instantaneous frequency curve is basically overlapped with the actual rotating speed variation curve after comparison, so that the fitting effect is proved to be good.
The phase discrimination time scale is obtained by fitting the instantaneous frequency curve obtained in fig. 15, and the outer ring fault vibration signal in fig. 13 is subjected to equal-angle resampling by using an interpolation method. Fig. 16 is an angular domain signal obtained by performing SST order ratio analysis on the fitted instantaneous frequency curve of fig. 15, wherein the abscissa is the number of revolutions, and the ordinate is the amplitude, and it can be seen from observing the angular domain waveform that there is an obvious periodic impact in the angular domain signal, and the angular domain signal substantially conforms to the outer ring fault vibration signal characteristics collected at a constant rotational speed in a general analysis.
Fig. 17 shows that the angular domain signal of fig. 16 is subjected to Hilbert envelope demodulation to obtain an envelope of the angular domain signal and an envelope order spectrum thereof, and an order of 3.036 and 2-5 harmonic orders thereof are found to be substantially consistent with an outer ring fault order of 3.05, so that the angular domain signal can be determined as a bearing outer ring fault feature.
As shown by simulation verification and experimental verification, the SST order ratio analysis method can effectively and accurately extract instantaneous frequency conversion information contained in the vibration signals, and converts unsteady vibration signals into quasi-steady angular domain signals by using the extracted rotating speed change ridge line, so that order ratio analysis is realized, and rolling bearing faults are successfully diagnosed.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (7)
1. A synchronous compression transformation order ratio analysis method for fault diagnosis of a rolling bearing is characterized by comprising the following steps:
step 1, carrying out low-pass filtering and down-sampling treatment on the collected vibration signals of the variable rotating speed of the bearing;
step 2, SST time frequency analysis processing is carried out on the processed vibration signal to obtain the time frequency spectrum thereof, the instantaneous frequency in the vibration signal is estimated, and the instantaneous frequency ridge line f of the frequency conversion is extractedcs(k) Namely the instantaneous frequency of the rotating speed, wherein k is the serial number of a sampling point;
step 3, extracting the ridge line f of the instantaneous frequencycs(k) Performing high-order polynomial fitting to obtain a fitted instantaneous frequency curve fi(t);
Step 4, utilizing the obtained fitting instantaneous frequency curve fi(T) carrying out interpolation filtering, and calculating to obtain a phase discrimination time scale sequence Tn;
Step 5, according to the phase discrimination time scale sequence TnCarrying out equal-angle resampling on the original vibration signal to obtain a quasi-steady-state angular domain signal sequence x (T)n);
Step 6, diagonal domain signal sequence x (T)n) Further carrying out Hilbert envelope demodulation, and solving an envelope order spectrum through FFT (fast Fourier transform), so as to realize order ratio analysis on the vibration signal;
and 7, judging whether the rolling bearing has faults and fault types thereof by analyzing the envelope order spectrum, giving a diagnosis result, proposing a maintenance suggestion, and completing the fault diagnosis of the bearing in the subway speed change process.
2. The synchronous compression transform order ratio analysis method of claim 1, wherein in step 1, the cut-off frequency of the low-pass filtering cannot be lower than the maximum value of the conversion frequency.
3. The synchronous compression transformation order ratio analysis method as claimed in claim 1, wherein in step 2, the peak search is used on the time frequency spectrum to obtain the instantaneous frequency ridge line f of the frequency conversioncs(k)。
4. The synchronous compression transformation order ratio analysis method as claimed in claim 1, wherein in step 3, the instantaneous frequency ridge line f is measuredcs(k) When fitting is carried out, high-order polynomial fitting and spline fitting are adopted in a small range.
5. The synchronous compressive transformation order ratio analysis method of claim 1, wherein in the step 3, a second-order polynomial fitting is adopted, and the fitting equation is as follows:
fi(t)=at2+bt+c (1)
assuming ε is the error of the fitting equation, the sum of the squares of ε is:
if so:andfitting coefficients a, b and c can be obtained.
6. The synchronous compressive transformation order ratio analysis of claim 1Method, characterized in that in step 4, a phase discrimination time scale sequence T is calculatednThe formula used is as follows:
in the formula, N is 1,2,3, N is a sampling time sequence number, and N is a sampling sequence length; delta theta pi/Omax,OmaxIs the maximum theoretical order;
integrating equation (3) yields:
get T0And (5) further calculating to obtain a phase discrimination time scale sequence T required by equal-angle resamplingn。
7. The synchronous compression transformation order ratio analysis method as claimed in claim 1, wherein in step 5, a phase-detecting time scale sequence T is usednAccording to the interpolation principle, original data are subjected to equal-angle resampling, and unsteady-state time domain signals are converted into quasi-steady-state angular domain signal sequences x (T)n):
In the formula,. DELTA.tsIs a time domain sampling interval; m Δ tsIs at TnA value in the vicinity; h issAnd (t) is an interpolation filter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910687690.5A CN110617964A (en) | 2019-07-29 | 2019-07-29 | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910687690.5A CN110617964A (en) | 2019-07-29 | 2019-07-29 | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110617964A true CN110617964A (en) | 2019-12-27 |
Family
ID=68921523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910687690.5A Pending CN110617964A (en) | 2019-07-29 | 2019-07-29 | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110617964A (en) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111307460A (en) * | 2020-03-14 | 2020-06-19 | 中国石化销售股份有限公司华南分公司 | Rolling bearing fault diagnosis method based on calculation order tracking and spectral kurtosis |
CN111307452A (en) * | 2020-03-05 | 2020-06-19 | 江苏天沃重工科技有限公司 | Intelligent fault diagnosis method for rotating machinery at time-varying rotating speed |
CN111368713A (en) * | 2020-03-02 | 2020-07-03 | 西南交通大学 | Vehicle network system harmonic time-frequency analysis method based on synchronous compression wavelet transform |
CN111665051A (en) * | 2020-07-01 | 2020-09-15 | 天津大学 | Bearing fault diagnosis method under strong noise variable-speed condition based on energy weight method |
CN111779573A (en) * | 2020-06-28 | 2020-10-16 | 河南柴油机重工有限责任公司 | Diesel engine online fault detection method and device |
CN111797789A (en) * | 2020-07-10 | 2020-10-20 | 合肥工业大学 | Time-frequency ridge line extraction method and device for fault feature interpretation |
CN111855208A (en) * | 2020-07-02 | 2020-10-30 | 南京航空航天大学 | Rotary machine weak fault signal extraction method based on order analysis and sparse coding |
CN111975784A (en) * | 2020-09-03 | 2020-11-24 | 昆明理工大学 | Joint robot fault diagnosis method based on current and vibration signals |
CN112484998A (en) * | 2020-11-16 | 2021-03-12 | 苏州大学文正学院 | Wind turbine generator bearing fault diagnosis method based on synchronous modal spectrum |
CN112577746A (en) * | 2020-12-07 | 2021-03-30 | 东南大学 | Method for extracting rolling bearing envelope order spectrum fault features under rotation speed fluctuation |
CN112665851A (en) * | 2020-10-28 | 2021-04-16 | 西安交通大学 | Key-free phase change rotating speed gearbox fault diagnosis method |
CN112710488A (en) * | 2020-12-08 | 2021-04-27 | 重庆川仪自动化股份有限公司 | Stress wave signal acquisition and analysis method for mechanical power equipment |
CN112781709A (en) * | 2020-12-24 | 2021-05-11 | 上海核工程研究设计院有限公司 | Method for analyzing early failure and extracting characteristics of equipment vibration signal under variable speed working condition |
CN113358212A (en) * | 2021-06-21 | 2021-09-07 | 重庆理工大学 | Electromechanical fault diagnosis method and system based on relative harmonic order and modeling method |
CN113702046A (en) * | 2021-09-13 | 2021-11-26 | 长沙理工大学 | Bearing fault diagnosis method based on mobile equipment under variable rotating speed working condition |
CN113866631A (en) * | 2021-09-02 | 2021-12-31 | 浙江大学 | Cage type asynchronous motor rotor broken bar fault diagnosis method |
CN114383718A (en) * | 2021-12-13 | 2022-04-22 | 北京化工大学 | High-frequency blade passing frequency extraction method based on vibration signals of external casing of gas turbine |
CN114486304A (en) * | 2022-01-30 | 2022-05-13 | 中国铁道科学研究院集团有限公司 | Method and device for tracking rotating part based on rolling stock |
CN114972350A (en) * | 2022-08-01 | 2022-08-30 | 深圳市信润富联数字科技有限公司 | Method, device and equipment for detecting abnormality of mold and storage medium |
CN115808236A (en) * | 2023-02-02 | 2023-03-17 | 武汉理工大学 | Fault on-line monitoring and diagnosing method and device for marine turbocharger and storage medium |
CN116257739A (en) * | 2023-05-16 | 2023-06-13 | 成都飞机工业(集团)有限责任公司 | Rapid visual diagnosis method for high-speed motorized spindle |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103884502A (en) * | 2014-04-02 | 2014-06-25 | 清华大学 | Method for diagnosing faults of planetary gear system of wind driven generator under variable rotating speed |
CN104266747A (en) * | 2014-06-09 | 2015-01-07 | 中能电力科技开发有限公司 | Fault diagnosis method based on vibration signal order analysis |
CN105446940A (en) * | 2015-11-23 | 2016-03-30 | 西安交通大学 | Amplitude correction method based on synchronous compression transformation and reconstruction |
CN105547698A (en) * | 2015-12-31 | 2016-05-04 | 新疆金风科技股份有限公司 | Fault diagnosis method and apparatus for rolling bearing |
CN105572501A (en) * | 2015-12-17 | 2016-05-11 | 西安理工大学 | Power quality disturbance identification method based on SST conversion and LS-SVM |
CN107525674A (en) * | 2017-05-27 | 2017-12-29 | 苏州大学 | Frequency method of estimation and detection means are turned based on crestal line probability distribution and localised waving |
CN108120597A (en) * | 2017-12-15 | 2018-06-05 | 天津金岸重工有限公司 | Seat type crane hoisting mechanism fault signature extracting method under variable speed |
CN108225764A (en) * | 2017-12-05 | 2018-06-29 | 昆明理工大学 | It is a kind of based on the high-precision of envelope extraction without key signal Order Tracking and system |
CN108760294A (en) * | 2018-06-13 | 2018-11-06 | 苏州大学 | Rotating machinery auto-check system and method |
CN108776031A (en) * | 2018-03-21 | 2018-11-09 | 南京航空航天大学 | A kind of rotary machinery fault diagnosis method based on improved synchronous extruding transformation |
CN109239780A (en) * | 2017-07-10 | 2019-01-18 | 中国石油化工股份有限公司 | Based on the synchronous method for squeezing wavelet transformation removal surface wave |
-
2019
- 2019-07-29 CN CN201910687690.5A patent/CN110617964A/en active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103884502A (en) * | 2014-04-02 | 2014-06-25 | 清华大学 | Method for diagnosing faults of planetary gear system of wind driven generator under variable rotating speed |
CN104266747A (en) * | 2014-06-09 | 2015-01-07 | 中能电力科技开发有限公司 | Fault diagnosis method based on vibration signal order analysis |
CN105446940A (en) * | 2015-11-23 | 2016-03-30 | 西安交通大学 | Amplitude correction method based on synchronous compression transformation and reconstruction |
CN105572501A (en) * | 2015-12-17 | 2016-05-11 | 西安理工大学 | Power quality disturbance identification method based on SST conversion and LS-SVM |
CN105547698A (en) * | 2015-12-31 | 2016-05-04 | 新疆金风科技股份有限公司 | Fault diagnosis method and apparatus for rolling bearing |
CN107525674A (en) * | 2017-05-27 | 2017-12-29 | 苏州大学 | Frequency method of estimation and detection means are turned based on crestal line probability distribution and localised waving |
CN109239780A (en) * | 2017-07-10 | 2019-01-18 | 中国石油化工股份有限公司 | Based on the synchronous method for squeezing wavelet transformation removal surface wave |
CN108225764A (en) * | 2017-12-05 | 2018-06-29 | 昆明理工大学 | It is a kind of based on the high-precision of envelope extraction without key signal Order Tracking and system |
CN108120597A (en) * | 2017-12-15 | 2018-06-05 | 天津金岸重工有限公司 | Seat type crane hoisting mechanism fault signature extracting method under variable speed |
CN108776031A (en) * | 2018-03-21 | 2018-11-09 | 南京航空航天大学 | A kind of rotary machinery fault diagnosis method based on improved synchronous extruding transformation |
CN108760294A (en) * | 2018-06-13 | 2018-11-06 | 苏州大学 | Rotating machinery auto-check system and method |
Non-Patent Citations (2)
Title |
---|
冯志鹏等: "《机械系统复杂非平稳信号分析方法原理及故障诊断》", 31 March 2018 * |
高威 等: ""同步压缩变换及其在机械振动信号处理中的应用"", 《时代农机》 * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111368713A (en) * | 2020-03-02 | 2020-07-03 | 西南交通大学 | Vehicle network system harmonic time-frequency analysis method based on synchronous compression wavelet transform |
CN111368713B (en) * | 2020-03-02 | 2022-06-28 | 西南交通大学 | Vehicle network system harmonic time-frequency analysis method based on synchronous compression wavelet transform |
CN111307452A (en) * | 2020-03-05 | 2020-06-19 | 江苏天沃重工科技有限公司 | Intelligent fault diagnosis method for rotating machinery at time-varying rotating speed |
CN111307460A (en) * | 2020-03-14 | 2020-06-19 | 中国石化销售股份有限公司华南分公司 | Rolling bearing fault diagnosis method based on calculation order tracking and spectral kurtosis |
CN111779573B (en) * | 2020-06-28 | 2022-02-11 | 河南柴油机重工有限责任公司 | Diesel engine online fault detection method and device |
CN111779573A (en) * | 2020-06-28 | 2020-10-16 | 河南柴油机重工有限责任公司 | Diesel engine online fault detection method and device |
CN111665051A (en) * | 2020-07-01 | 2020-09-15 | 天津大学 | Bearing fault diagnosis method under strong noise variable-speed condition based on energy weight method |
CN111855208A (en) * | 2020-07-02 | 2020-10-30 | 南京航空航天大学 | Rotary machine weak fault signal extraction method based on order analysis and sparse coding |
CN111855208B (en) * | 2020-07-02 | 2021-09-21 | 南京航空航天大学 | Rotary machine weak fault signal extraction method based on order analysis and sparse coding |
CN111797789B (en) * | 2020-07-10 | 2024-02-06 | 合肥工业大学 | Time-frequency ridge line extraction method and device for fault feature interpretation |
CN111797789A (en) * | 2020-07-10 | 2020-10-20 | 合肥工业大学 | Time-frequency ridge line extraction method and device for fault feature interpretation |
CN111975784A (en) * | 2020-09-03 | 2020-11-24 | 昆明理工大学 | Joint robot fault diagnosis method based on current and vibration signals |
CN112665851A (en) * | 2020-10-28 | 2021-04-16 | 西安交通大学 | Key-free phase change rotating speed gearbox fault diagnosis method |
CN112484998B (en) * | 2020-11-16 | 2022-12-27 | 苏州大学文正学院 | Wind turbine generator bearing fault diagnosis method based on synchronous modal spectrum |
CN112484998A (en) * | 2020-11-16 | 2021-03-12 | 苏州大学文正学院 | Wind turbine generator bearing fault diagnosis method based on synchronous modal spectrum |
CN112577746B (en) * | 2020-12-07 | 2022-07-22 | 东南大学 | Method for extracting rolling bearing envelope order spectrum fault features under rotation speed fluctuation |
CN112577746A (en) * | 2020-12-07 | 2021-03-30 | 东南大学 | Method for extracting rolling bearing envelope order spectrum fault features under rotation speed fluctuation |
CN112710488A (en) * | 2020-12-08 | 2021-04-27 | 重庆川仪自动化股份有限公司 | Stress wave signal acquisition and analysis method for mechanical power equipment |
CN112710488B (en) * | 2020-12-08 | 2023-01-20 | 重庆川仪自动化股份有限公司 | Stress wave signal acquisition and analysis method for mechanical power equipment |
CN112781709A (en) * | 2020-12-24 | 2021-05-11 | 上海核工程研究设计院有限公司 | Method for analyzing early failure and extracting characteristics of equipment vibration signal under variable speed working condition |
CN113358212A (en) * | 2021-06-21 | 2021-09-07 | 重庆理工大学 | Electromechanical fault diagnosis method and system based on relative harmonic order and modeling method |
CN113866631A (en) * | 2021-09-02 | 2021-12-31 | 浙江大学 | Cage type asynchronous motor rotor broken bar fault diagnosis method |
CN113702046A (en) * | 2021-09-13 | 2021-11-26 | 长沙理工大学 | Bearing fault diagnosis method based on mobile equipment under variable rotating speed working condition |
CN114383718A (en) * | 2021-12-13 | 2022-04-22 | 北京化工大学 | High-frequency blade passing frequency extraction method based on vibration signals of external casing of gas turbine |
CN114486304A (en) * | 2022-01-30 | 2022-05-13 | 中国铁道科学研究院集团有限公司 | Method and device for tracking rotating part based on rolling stock |
CN114486304B (en) * | 2022-01-30 | 2024-03-22 | 中国铁道科学研究院集团有限公司 | Rolling stock-based rotating part tracking method and device |
CN114972350A (en) * | 2022-08-01 | 2022-08-30 | 深圳市信润富联数字科技有限公司 | Method, device and equipment for detecting abnormality of mold and storage medium |
CN114972350B (en) * | 2022-08-01 | 2022-11-15 | 深圳市信润富联数字科技有限公司 | Method, device and equipment for detecting abnormality of mold and storage medium |
CN115808236B (en) * | 2023-02-02 | 2023-05-05 | 武汉理工大学 | Marine turbocharger fault on-line monitoring and diagnosing method and device and storage medium |
CN115808236A (en) * | 2023-02-02 | 2023-03-17 | 武汉理工大学 | Fault on-line monitoring and diagnosing method and device for marine turbocharger and storage medium |
CN116257739A (en) * | 2023-05-16 | 2023-06-13 | 成都飞机工业(集团)有限责任公司 | Rapid visual diagnosis method for high-speed motorized spindle |
CN116257739B (en) * | 2023-05-16 | 2023-08-04 | 成都飞机工业(集团)有限责任公司 | Rapid visual diagnosis method for high-speed motorized spindle |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110617964A (en) | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing | |
CN110763462B (en) | Time-varying vibration signal fault diagnosis method based on synchronous compression operator | |
Wang et al. | Rolling element bearing fault diagnosis via fault characteristic order (FCO) analysis | |
CN111665051A (en) | Bearing fault diagnosis method under strong noise variable-speed condition based on energy weight method | |
Urbanek et al. | Comparison of amplitude-based and phase-based method for speed tracking in application to wind turbines | |
CN111307460A (en) | Rolling bearing fault diagnosis method based on calculation order tracking and spectral kurtosis | |
CN105092241B (en) | A kind of gear local fault diagnosis method and system | |
CN110044610A (en) | Gear failure diagnosing method | |
CN112105907B (en) | Method and apparatus for monitoring a gear system | |
CN113125179A (en) | Keyless phase order tracking method for rotating speed fluctuation of rotary machine | |
CN110119764B (en) | Purification method of axis track under variable rotating speed working condition | |
CN108120597A (en) | Seat type crane hoisting mechanism fault signature extracting method under variable speed | |
JPH1026580A (en) | Method and device for diagnosing speed-change-type rotary mechanical equipment | |
CN108398260B (en) | Method for quickly evaluating instantaneous angular speed of gearbox based on mixed probability method | |
CN108278184B (en) | Wind turbine generator impeller imbalance monitoring method based on empirical mode decomposition | |
Lin et al. | A review and strategy for the diagnosis of speed-varying machinery | |
CN111307426A (en) | Rotating machinery fault feature extraction method based on FrFT-EWT principle | |
CN112465068A (en) | Rotating equipment fault feature extraction method based on multi-sensor data fusion | |
Zhu et al. | Fault diagnosis method for rolling element bearings under variable speed based on TKEO and fast-SC | |
CN113984282B (en) | Rotor dynamic balance fault feature extraction method under keyless phase-change rotating speed | |
CN117836599A (en) | Method for detecting bearing defects in a rotating system and monitoring system for implementing said method | |
CN112781723B (en) | Harmonic component detection method based on frequency spectrum variance | |
CN116256158A (en) | Rotary machine instantaneous phase self-adaptive extraction method based on depth signal separation | |
CN114112396B (en) | Bearing fault diagnosis method under rotating speed fluctuation quasi-stationary working condition by adopting current signal analysis | |
CN112784702A (en) | Signal processing method based on self-adaptive rotation synchronous extraction frequency modulation transformation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191227 |
|
RJ01 | Rejection of invention patent application after publication |