CN113702046A - Bearing fault diagnosis method based on mobile equipment under variable rotating speed working condition - Google Patents
Bearing fault diagnosis method based on mobile equipment under variable rotating speed working condition Download PDFInfo
- Publication number
- CN113702046A CN113702046A CN202111068927.5A CN202111068927A CN113702046A CN 113702046 A CN113702046 A CN 113702046A CN 202111068927 A CN202111068927 A CN 202111068927A CN 113702046 A CN113702046 A CN 113702046A
- Authority
- CN
- China
- Prior art keywords
- signal
- bearing
- fault
- envelope
- mobile equipment
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000003745 diagnosis Methods 0.000 title claims abstract description 12
- 238000001228 spectrum Methods 0.000 claims abstract description 14
- 230000005236 sound signal Effects 0.000 claims abstract description 9
- 230000003595 spectral effect Effects 0.000 claims abstract description 5
- 238000001914 filtration Methods 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 10
- 238000012952 Resampling Methods 0.000 claims description 4
- 230000010363 phase shift Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 244000062793 Sorghum vulgare Species 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 235000019713 millet Nutrition 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
Images
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
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention relates to a bearing fault diagnosis method under the variable rotating speed working condition based on mobile equipment, which comprises the steps of sending a carrier wave through a loudspeaker of the mobile equipment, demodulating the amplitude of a sound signal collected by a microphone, searching and restoring an instantaneous phase signal of a rotating shaft where a bearing is located through a maximum value and fitting a least square polynomial, extracting a fault impact envelope through a spectral kurtosis method to carry out order tracking, and finally analyzing and judging the fault of the bearing through an envelope order spectrum.
Description
Technical Field
The invention relates to the technical field of mechanical fault diagnosis and signal processing analysis, in particular to a bearing fault diagnosis method under a variable-rotating-speed working condition based on mobile equipment.
Background
The bearing is a widely used mechanical part, and the health monitoring of the bearing is very important for the safe and stable operation of rotary mechanical equipment. Vibration and noise analysis is a common means for bearing fault diagnosis, generally needs a professional data acquisition card and sensor hardware, and is too expensive for non-key equipment objects. In addition, a large number of bearings work under the working condition of variable rotating speed, and the bearing fault diagnosis under the working condition of variable rotating speed also depends on speed measuring hardware such as a photoelectric probe and an encoder, so that the test cost is increased.
In the document [1], an author acquires radiation noise of a bearing through a built-in microphone of a mobile phone and diagnoses bearing faults through analysis and processing of acoustic signals, the technology is favorable for popularization and application of a bearing fault diagnosis technology to a certain extent, but the method is only suitable for a constant rotating speed working condition and needs to know rotating speed parameters of the bearing in advance.
The patent with patent application number 201910416881.8 realizes the measurement of the rotating speed of the rotating shaft by using the mobile device, and the invention innovatively utilizes the built-in loudspeaker of the mobile device to actively sound, and simultaneously performs signal processing on the sound signal collected by the built-in microphone to realize the estimation of the speed, but the method is only suitable for the rotating speed estimation under the working condition of fixed rotating speed.
[1]Orman M,Rzeszucinski P,Tkaczyk A,et al.Bearing fault detection with the use of acoustic signals recorded by a hand-held mobile phone[C]//2015International Conference on Condition Assessment Techniques in Electrical Systems(CATCON).IEEE,2015:252-256.
Disclosure of Invention
The invention aims to enable field equipment maintainers to use mobile equipment (smart phones and tablet computers) to conveniently diagnose the bearing fault by the proposed invention method without additional special test hardware. The bearing fault under the variable-speed working condition is accurately diagnosed by utilizing the built-in loudspeaker and microphone hardware of the mobile equipment and the calculation display capacity through a signal processing method of envelope demodulation, least square fitting, spectral kurtosis, order tracking and envelope order spectrum.
In order to achieve the purpose, the invention adopts the technical scheme that:
the bearing fault diagnosis method based on the variable rotating speed working condition of the mobile equipment comprises the following steps:
the method comprises the following steps: closely aligning a loudspeaker of the mobile equipment to a rotating shaft where a test bearing is positioned, and sending out a set carrier frequency f through the loudspeakercThe microphone records the received sound signal, denoted as x (t);
step two: constructing a zero phase shift band-pass filter, wherein the center frequency of the band-pass filter is the carrier frequency fcThe bandwidth is 2 times of the limit frequency of the measured bearing, and the signal x (t) is subjected to band-pass filtering to obtain a filtering signal xF(t), filtering the signal xF(t) obtaining an amplitude modulation signal x by Hilbert transformE(t);
Step three: searching for modulated signal xEThe maximum point in (t) is taken as the zero crossing { x ] in the key phase signalE(t1),xE(t2),…,xE(tn) And (4) when the rotating shaft rotates for one circle every time of zero crossing, obtaining an instantaneous phase signal theta of the rotating shaft where the bearing is located by least square polynomial fitting to eliminate zero point estimation errorsS(t) converting the signal θS(t) differentiating to obtain an instantaneous rotating speed curve S (t) of the shaft;
step four: determining the frequency band of the fault impact component in the signal x (t) by a spectral Kurtogram (Kurtogram) method, and performing band-pass filtering and Hilbert transform on the frequency band to extract the fault impact envelope xI(t);
Step five: envelope the fault impact xI(t) combining the instantaneous phase signal θS(t) carrying out order tracking to obtain the fault impact envelope x after angle domain resamplingI(θS) And further obtaining an envelope order spectrum through Fourier transform, and judging whether a fault occurs and the type according to whether the bearing fault characteristic order occurs in the envelope order spectrum.
The invention has the following beneficial effects:
a) the method of the invention uses the mobile device to accurately diagnose the bearing fault under the variable rotating speed working condition, does not need to use extra vibration noise test equipment and speed measurement hardware, has low test cost, and is beneficial to technical popularization and use of field maintainers in the form of application of mobile phone programs.
b) The bearing fault characteristic order is constant and is irrelevant to the working rotating speed of the bearing, and compared with the traditional envelope spectrum, the existence and the type of the fault can be easily determined through the envelope order spectrum.
c) For bearing fault diagnosis under the working condition of constant rotating speed, the method is also effective, and the working rotating speed of the bearing is not required to be known in advance.
Drawings
FIG. 1 is a schematic test diagram of an embodiment of the present invention.
FIG. 2 is a flow chart of the present invention.
Fig. 3 shows a sound signal x (t) recorded by a microphone according to an embodiment of the invention.
FIG. 4 shows an amplitude modulation signal x according to an embodiment of the present inventionE(t)。
FIG. 5 is a diagram showing an instantaneous phase signal θ of a rotating shaft on which a bearing according to an embodiment of the present invention is locatedS(t)。
FIG. 6 is a graph S (t) showing the instantaneous rotational speed of the rotating shaft on which the bearing of the embodiment of the present invention is located.
FIG. 7 shows a fault strike envelope signal x according to an embodiment of the inventionI(t)。
FIG. 8 shows a fault strike envelope signal x according to an embodiment of the inventionI(t) spectrum.
FIG. 9 shows a fault impact envelope signal x after angle resampling according to an embodiment of the present inventionI(θS)。
FIG. 10 is a graph of an envelope order spectrum according to an embodiment of the present invention.
FIG. 11 is a photograph of a bearing failure in accordance with an embodiment of the present invention.
Detailed Description
The invention is described in detail below with reference to the figures and examples.
Fig. 1 is a schematic test diagram of an embodiment, where 1 is a mobile device, 2 is a rotation axis where a faulty bearing is located, 3 is the faulty bearing, 4 is a carrier component and a modulation echo transmitted by a mobile phone, and 5 is a fault impact radiated by the faulty bearing.
The specific parameters are as follows: 1) the mobile equipment 1 is a mobile phone (millet 11), and the distance between the mobile equipment and the surface of the rotating shaft 2 is 5 cm; 2) the rotating shaft 2 is driven to rotate by the servo motor, and the rotating speed is increased from about 340RPM to about 730RPM in 10 seconds; 3) structural parameters of the failed bearing 3: a) roller diameter: 12.3mm, b) pitch diameter: 52mm, c) contact angle: 0 °, d) number of rollers: 8; 4) the frequency of the carrier wave sent by the mobile phone is 5 kHz; 5) the mobile phone carries out high-frequency sampling and data processing on the sound signal, wherein the sampling frequency of the sound signal is 44100Hz, and the sampling time duration is 10 seconds.
The method for analyzing the sound signal and diagnosing the bearing fault is applied, and as shown in figure 2, the method for diagnosing the bearing fault under the variable rotating speed working condition based on the mobile equipment comprises the following steps:
the method comprises the following steps: aligning a loudspeaker of the mobile equipment to the surface of a rotating shaft where a test bearing is located at a distance of 5cm, and emitting a set carrier frequency f through the loudspeakercA 5000Hz sine wave, and a microphone records a received sound signal, which is denoted as x (t), and x (t) in this embodiment is shown in fig. 3;
step two: constructing a zero phase shift band-pass filter, wherein the center frequency of the band-pass filter is the carrier frequency fcThe bandwidth is 5000Hz and is 2 times of the limit frequency of the measured bearing, the bandwidth is 60Hz in the embodiment, and the signal x (t) is subjected to band-pass filtering to obtain a filtering signal xF(t), filtering the signal xF(t) obtaining an amplitude modulation signal x by Hilbert transformE(t), this example xE(t) as shown in FIG. 4;
step three: searching for modulated signal xEThe maximum point in (t) is taken as the zero crossing { x ] in the key phase signalE(t1),xE(t2),…,xE(tn) And the zero-crossing points are marked by asterisks in figure 4, each time a zero-crossing point is passed, the rotating shaft rotates for one circle, and in order to eliminate zero-point estimation errors, an instantaneous phase signal theta of the rotating shaft where the bearing is located is obtained through least square polynomial fittingS(t), example θS(t) As shown in FIG. 5, the signal θS(t) differentiating to obtain an instantaneous speed curve s (t) of the shaft, where the instantaneous speed curve s (t) is consistent with a speed regulation curve of the set motor as shown in fig. 6 in this embodiment;
step (ii) ofFourthly, the method comprises the following steps: determining the frequency band of the fault impact component in the signal x (t) by a spectral Kurtogram (Kurtogram) method, and performing band-pass filtering and Hilbert transform on the frequency band to extract the fault impact envelope xI(t), this example xI(t) impact envelope x, as shown in FIG. 7I(t) the frequency spectrum of the rotation speed tends to become dense in the time domain along with the increase of the rotation speed, as shown in fig. 8, obvious frequency blurring can be seen, and the bearing fault cannot be judged according to the traditional envelope analysis method;
step five: envelope the fault impact xI(t) combining the instantaneous phase signal θS(t) carrying out order tracking to obtain the fault impact envelope x after angle domain resamplingI(θS) Example xI(θS) As shown in fig. 9, the resampled impact envelopes are distributed at equal intervals in an angle domain, and then an envelope order spectrum is obtained through fourier transform, the envelope order spectrum of this embodiment is shown in fig. 10, harmonic waves of the bearing outer ring fault order appear in the envelope order spectrum, and therefore it can be judged that the bearing has an outer ring fault and is consistent with the actual fault of the bearing, and a fault bearing photograph is shown in fig. 11.
The bearing fault diagnosis method based on the variable-speed working condition of the mobile equipment does not need special test instrument equipment, fully exerts the loudspeaker and microphone of the mobile equipment and the calculation and display capacity, extracts fault characteristic information, and effectively diagnoses the bearing fault under the variable-speed working condition.
Claims (1)
1. The bearing fault diagnosis method based on the variable rotating speed working condition of the mobile equipment comprises the following steps:
the method comprises the following steps: closely aligning a loudspeaker of the mobile equipment to a rotating shaft where a test bearing is positioned, and sending out a set carrier frequency f through the loudspeakercThe microphone records the received sound signal, denoted as x (t);
step two: constructing a zero phase shift band-pass filter, wherein the center frequency of the band-pass filter is the carrier frequency fcThe bandwidth is 2 times of the limit frequency of the measured bearing, and the signal x (t) is subjected to band-pass filtering to obtain a filtering signal xF(t), filtering the signal xF(t) obtaining an amplitude modulation signal x by Hilbert transformE(t);
Step three: searching for modulated signal xEThe maximum point in (t) is taken as the zero crossing { x ] in the key phase signalE(t1),xE(t2),…,xE(tn) And (4) when the rotating shaft rotates for one circle every time of zero crossing, obtaining an instantaneous phase signal theta of the rotating shaft where the bearing is located by least square polynomial fitting to eliminate zero point estimation errorsS(t) converting the signal θS(t) differentiating to obtain an instantaneous rotating speed curve S (t) of the shaft;
step four: determining the frequency band of the fault impact component in the signal x (t) by a spectral Kurtogram (Kurtogram) method, and performing band-pass filtering and Hilbert transform on the frequency band to extract the fault impact envelope xI(t);
Step five: envelope the fault impact xI(t) combining the instantaneous phase signal θS(t) carrying out order tracking to obtain the fault impact envelope x after angle domain resamplingI(θS) And further obtaining an envelope order spectrum through Fourier transform, and judging whether a fault occurs and the type according to whether the bearing fault characteristic order occurs in the envelope order spectrum.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111068927.5A CN113702046B (en) | 2021-09-13 | 2021-09-13 | Bearing fault diagnosis method under variable rotation speed working condition based on mobile equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111068927.5A CN113702046B (en) | 2021-09-13 | 2021-09-13 | Bearing fault diagnosis method under variable rotation speed working condition based on mobile equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113702046A true CN113702046A (en) | 2021-11-26 |
CN113702046B CN113702046B (en) | 2024-06-11 |
Family
ID=78660117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111068927.5A Active CN113702046B (en) | 2021-09-13 | 2021-09-13 | Bearing fault diagnosis method under variable rotation speed working condition based on mobile equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113702046B (en) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103048137A (en) * | 2012-12-20 | 2013-04-17 | 北京航空航天大学 | Fault diagnosis method of rolling bearing under variable working conditions |
EP2810027A1 (en) * | 2012-01-30 | 2014-12-10 | SPM Instrument AB | Apparatus and method for analysing the condition of a machine having a rotating part |
CN105784366A (en) * | 2016-03-30 | 2016-07-20 | 华北电力大学(保定) | Wind turbine generator bearing fault diagnosis method under variable speed |
CN106053871A (en) * | 2016-07-25 | 2016-10-26 | 昆明理工大学 | Method for rotation speed extraction through peeling off fault corresponding impact by employing rolling ball track |
CN110108467A (en) * | 2019-05-20 | 2019-08-09 | 长沙理工大学 | Active sounding speed-measuring method based on portable mobile apparatus |
CN110617964A (en) * | 2019-07-29 | 2019-12-27 | 中国铁道科学研究院集团有限公司城市轨道交通中心 | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing |
CN110887663A (en) * | 2019-10-30 | 2020-03-17 | 中国石油化工股份有限公司 | Bearing fault diagnosis method combining variable working condition calculation order tracking and spectral kurtosis |
CN111307460A (en) * | 2020-03-14 | 2020-06-19 | 中国石化销售股份有限公司华南分公司 | Rolling bearing fault diagnosis method based on calculation order tracking and spectral kurtosis |
CN111553178A (en) * | 2019-02-12 | 2020-08-18 | 斯凯孚公司 | Intelligent identification method for vibration characteristics of rotating machinery |
CN111665051A (en) * | 2020-07-01 | 2020-09-15 | 天津大学 | Bearing fault diagnosis method under strong noise variable-speed condition based on energy weight method |
CN112577746A (en) * | 2020-12-07 | 2021-03-30 | 东南大学 | Method for extracting rolling bearing envelope order spectrum fault features under rotation speed fluctuation |
CN113125179A (en) * | 2021-03-11 | 2021-07-16 | 同济大学 | Keyless phase order tracking method for rotating speed fluctuation of rotary machine |
-
2021
- 2021-09-13 CN CN202111068927.5A patent/CN113702046B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2810027A1 (en) * | 2012-01-30 | 2014-12-10 | SPM Instrument AB | Apparatus and method for analysing the condition of a machine having a rotating part |
CN103048137A (en) * | 2012-12-20 | 2013-04-17 | 北京航空航天大学 | Fault diagnosis method of rolling bearing under variable working conditions |
CN105784366A (en) * | 2016-03-30 | 2016-07-20 | 华北电力大学(保定) | Wind turbine generator bearing fault diagnosis method under variable speed |
CN106053871A (en) * | 2016-07-25 | 2016-10-26 | 昆明理工大学 | Method for rotation speed extraction through peeling off fault corresponding impact by employing rolling ball track |
CN111553178A (en) * | 2019-02-12 | 2020-08-18 | 斯凯孚公司 | Intelligent identification method for vibration characteristics of rotating machinery |
CN110108467A (en) * | 2019-05-20 | 2019-08-09 | 长沙理工大学 | Active sounding speed-measuring method based on portable mobile apparatus |
CN110617964A (en) * | 2019-07-29 | 2019-12-27 | 中国铁道科学研究院集团有限公司城市轨道交通中心 | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing |
CN110887663A (en) * | 2019-10-30 | 2020-03-17 | 中国石油化工股份有限公司 | Bearing fault diagnosis method combining variable working condition calculation order tracking and spectral kurtosis |
CN111307460A (en) * | 2020-03-14 | 2020-06-19 | 中国石化销售股份有限公司华南分公司 | Rolling bearing fault diagnosis method based on calculation order tracking and spectral kurtosis |
CN111665051A (en) * | 2020-07-01 | 2020-09-15 | 天津大学 | Bearing fault diagnosis method under strong noise variable-speed condition based on energy weight method |
CN112577746A (en) * | 2020-12-07 | 2021-03-30 | 东南大学 | Method for extracting rolling bearing envelope order spectrum fault features under rotation speed fluctuation |
CN113125179A (en) * | 2021-03-11 | 2021-07-16 | 同济大学 | Keyless phase order tracking method for rotating speed fluctuation of rotary machine |
Non-Patent Citations (3)
Title |
---|
张亢;程军圣;: "基于LMD和阶次跟踪分析的滚动轴承故障诊断", 振动.测试与诊断, no. 03, 15 June 2016 (2016-06-15) * |
武英杰;辛红伟;王建国;王晓龙;: "基于VMD滤波和极值点包络阶次的滚动轴承故障诊断", 振动与冲击, no. 14, 28 July 2018 (2018-07-28) * |
赵德尊: "变转速下滚动轴承时变非平稳故障特征提取方法研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》, no. 11, 15 November 2018 (2018-11-15), pages 029 - 1 * |
Also Published As
Publication number | Publication date |
---|---|
CN113702046B (en) | 2024-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104655380B (en) | A kind of rotating machinery fault signature extracting method | |
CN110108467B (en) | Active sounding speed measurement method based on portable mobile equipment | |
Lei et al. | EEMD method and WNN for fault diagnosis of locomotive roller bearings | |
CN106226078B (en) | A kind of bearing calibration of Doppler's distortion acoustic signal based on microphone array | |
CN110617964A (en) | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing | |
CN108572075A (en) | Fault Diagnosis of Roller Bearings based on Wavelet Packet Energy Spectrum and modulation double-spectrum analysis | |
Liu et al. | Rolling bearing fault diagnosis via STFT and improved instantaneous frequency estimation method | |
CN109946081B (en) | Fault diagnosis method for rolling bearing slipping under variable rotating speed | |
CN104215456B (en) | Plane clustering and frequency-domain compressed sensing reconstruction based mechanical fault diagnosis method | |
CN106769041B (en) | Permanent magnet synchronous motor bearing online fault diagnosis device and method under variable rotating speed working condition | |
CN105509876A (en) | Undersampled leaf apex timing vibration signal reconstruction method and device thereof | |
CN110057583A (en) | A kind of bearing fault recognition methods, device and computer equipment | |
CN103308305A (en) | Gear fault diagnosis device and method | |
CN114739667B (en) | Multi-mode information fusion bearing lubrication state monitoring device and method | |
JP2009115537A (en) | Vibration measuring method | |
JP5765881B2 (en) | Equipment diagnostic equipment | |
Peter et al. | The sparsogram: A new and effective method for extracting bearing fault features | |
CN206300859U (en) | A kind of test system that material friction performance is tested and analyzed based on acoustic emission signal | |
CN113702046A (en) | Bearing fault diagnosis method based on mobile equipment under variable rotating speed working condition | |
CN104090126A (en) | Accelerometer bandwidth testing method | |
CN116358873B (en) | Rolling bearing fault diagnosis method, device, equipment and medium | |
CN203011692U (en) | Portable fault diagnosing device for rolling bearing | |
CN1049285C (en) | Trouble testing technique with multiple sensors, resonance and demodulation | |
Xu et al. | Mobile device-based bearing diagnostics with varying speeds | |
CN107271184B (en) | A kind of the kernel regression decomposition method and system of rolling bearing fault diagnosis |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |