NZ554049A - Method and system for vibration signal processing - Google Patents
Method and system for vibration signal processingInfo
- Publication number
- NZ554049A NZ554049A NZ55404907A NZ55404907A NZ554049A NZ 554049 A NZ554049 A NZ 554049A NZ 55404907 A NZ55404907 A NZ 55404907A NZ 55404907 A NZ55404907 A NZ 55404907A NZ 554049 A NZ554049 A NZ 554049A
- Authority
- NZ
- New Zealand
- Prior art keywords
- vibration
- samples
- digital
- signals
- electrical
- Prior art date
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
A method of vibration signal processing is disclosed. The method comprising the steps of: (a) obtaining electrical vibration signals having amplitudes and frequencies that represent mechanical vibrations of a machine; (b) converting the electrical vibration signals into digital vibration samples; (c) dividing the digital vibration samples into equal time intervals; (d) determining the average absolute amplitude of the digital vibration samples for each time interval; (e) generating a time waveform comprising the determined average amplitudes; and (f) processing the time waveform as if it were an independently-detected signal.
Description
*10055395177
NEW ZEALAND PATENTS ACT, 1953
No: 554049 Date: 22 March 2007
COMPLETE SPECIFICATION
METHOD AND SYSTEM FOR VIBRATION SIGNAL PROCESSING
We, COMMTEST INSTRUMENTS LIMITED, a New Zealand company of 28B Moorhouse Avenue, Christchurch 8001, New Zealand, do hereby declare the invention for which we pray that a patent may be granted to us, and the method by which it is to be performed, to be particularly described in and by the following statement:
IPONZ
2 9 JAr, 2U08
2
FIELD OF THE INVENTION
The present invention relates to a method and system for vibration signal processing. In particular, but not exclusively, the present invention relates to a method and system for 5 machine fault detection using a demodulation process to account for situations where impact or impulse events are hidden in a normal vibration reading.
BACKGROUND TO THE INVENTION
Demodulation is a vibration analysis process used to account for situations where certain machine problems remain hidden in a normal vibration reading, preventing early corrective action. Usually, these situations arise when the vibration is associated with an impact or impulse event. These events have a high frequency but a short duration. Unfortunately, when taking vibration readings, these vibration events are overwhelmed by the prominent 15 long duration of low-frequency vibrations.
In use, demodulation first removes the low-frequency components of the vibration signals using a high-pass filter. To make the impulse events more prominent and easily identifiable later, the high-pass filtered signal is envelope-detected. The envelope-detected signal is 20 then processed as if it were an independently-detected vibration signal. One disadvantage in demodulation is that, due to the number of filter stages, its digital implementation requires intensive signal processing.
Instead of envelope-detecting the high-pass filtered vibration signal, a technique known as 25 PeakVue (disclosed in US Patent No. 5,895,857 to Robinson et al.) exists to extract the impulse events by signal processing the peaks of a detected vibration signal. The reference to the Robinson et al. patent and any other reference to patent specifications, external documents, or sources of information in this specification has been provided for the purpose of providing a context for discussing the features of the invention. Unless 30 specifically stated otherwise, reference to such documents or sources of information is not to be construed as an admission that such documents or sources of information are prior art or form part of the common general knowledge in the art in any jurisdiction.
1310770
It is an object of the present invention to either provide a method and system that implements a demodulation type process, with a high signal-to-noise ratio (SNR), without intensive signal processing, or at least provide the public with a useful choice.
SUMMARY OF THE INVENTION
In one aspect, the present invention relates to a method of vibration signal processing comprising the steps of: obtaining electrical vibration signals having amplitudes and 10 frequencies that represent mechanical vibrations of a machine; converting the electrical vibration signals into digital vibration samples; dividing the digital vibration samples into equal time intervals; determining the average absolute amplitude of the digital vibration samples for each time interval; generating a time waveform comprising the determined average amplitudes; and processing the time waveform as if it were an independently-15 detected signal.
In another aspect, the present invention relates to a method of vibration signal processing comprising the steps of: measuring the vibration level of a machine using one or more sensors to obtain electrical vibration signals; conditioning the electrical vibration signals; 20 high pass filtering the electrical vibration signals; converting the high pass filtered electrical vibration signals into digital vibration samples; dividing the digital vibration samples into equal time intervals and, for each time interval, rectifying the digital vibration samples and determining the arithmetic mean or the root mean square (RMS) of the rectified digital vibration samples; buffering the mean digital vibration samples to create a time waveform 25 representing the mean digital vibration samples; and processing the time waveform as if it were an independently-detected signal.
Preferably, the step of conditioning the electrical vibration signals comprises: buffering the electrical vibration signals; adjusting the gain of the buffered electrical vibration signals; and 30 low-pass filtering the gain-adjusted electrical vibration signals.
Preferably, the step of high-pass filtering includes the step of selecting a lowest frequency of interest, Dmm.
1310770
Preferably, the step of converting the electrical vibration signals into digital vibration samples includes the steps of:
selecting the highest frequency of interest, Dmax, and converting the filtered electrical vibration signals into digital vibration samples at a sampling rate of typically Dmax x 2.56.
Preferably, the step of dividing the digital vibration samples into equal time intervals includes the steps of:
selecting a maximum fault frequency of interest, DFmax;
dividing the digital vibration samples into equal time intervals, each having approximately Dmax / DFmax number of samples.
Preferably, the step of processing the time waveform comprises the steps of:
windowing the time waveform;
transforming the windowed time waveform into a frequency spectrum; and determining an average of the frequency spectrum.
In another aspect, the present invention relates to a system for vibration signal processing comprising: one or more vibration sensors configured to sense mechanical vibrations and to represent the sensed mechanical vibrations as electrical vibration signals; a conditioning circuit to condition the electrical vibration signals; a high-pass filter circuit; an analogue-to-digital converter (ADC) to convert the electrical vibration signals from the high-pass filter circuit into digital vibration samples; an averaging circuit configured to divide the digital vibration samples into equal time intervals, the averaging circuit including a rectifier to rectify the digital vibration samples in each time interval, and an averager to determine the arithmetic mean or the root mean square (RMS) of the rectified digital vibration samples in each time interval; and a buffer to receive the mean of the rectified digital vibration samples of each time interval, and to output the samples as a time waveform for further analysis.
Preferably, the averaging circuit and buffer are implemented in a processor.
Preferably, the processor is also configured to: window the time waveform; transform the windowed time waveform into a frequency spectrum; and determine an average of the frequency spectrum. These steps could also be carried out as part of a second processing chain on another processor.
Preferably, the conditioning circuit comprises a buffer, a gain-adjuster, and a low-pass filter.
Preferably, the system includes an input means for a user to select a lowest frequency of interest for the high-pass filter.
Preferably, the system includes an input means for a user to select the highest frequency of interest for use in setting the sampling rate of the ADC at typically D/m. X 2.56.
Preferably, the system includes an input means for a user to select a maximum fault 15 frequency of interest (DF//M.) for use in setting each of the time intervals with a length of approximately Daax / DFmax.
The term 'comprising' as used in this specification means 'consisting at least in part of, that is to say when interpreting statements in this specification which include that term, the 20 features, prefaced by that term in each statement, all need to be present but other features can also be present. Related terms such as 'comprise' and 'comprised' are to be interpreted in similar manner.
The term "windowing" is well known to those skilled in the art. As used in this 25 specification the term means the application of a window function to a signal. In signal processing, a window function is a function that is zero-valued outside of some chosen interval. When another function or a signal is multiplied by a window function, the product is also zero-valued outside the interval.
As used herein the term "(s)" following a noun means the plural and/or singular form of that noun.
As used herein the term "and/or" means "and" or "or", or where the context allows both.
1310770
To those skilled in the art to which the invention relates, many changes in construction and widely differing embodiments and applications of the invention will suggest themselves without departing from the scope of the invention as defined in the appended claims. The disclosures and the descriptions herein are purely illustrative and are not intended to be in any sense limiting.
The invention consists in the foregoing and also envisages constructions of which the following gives examples only.
BRIEF DESCRIPTION OF THE FIGURES
Preferred forms of the method and system of the present invention will now be described with reference to the accompanying figures in which:
Figure 1 shows a simplified processing chain of the system of the present invention,
Figure 2A shows a worst-case type of test signal used for testing some of the prior techniques and that of the present invention,
Figure 2B shows the results of a prior art demodulation process carried out on the test signal of Figure 2A,
Figure 2C shows the results of another prior art demodulation process carried out on the test signal of Figure 2A, and
Figure 2D shows the results of applying the technique of the present invention to the test signal of Figure 2A.
7
DETAILED DESCRIPTION OF THE PREFERRED FORMS
The present invention relates to vibration signal analysis that can be used to account for situations where impact or impulse events are hidden in a normal vibration reading by low-5 frequency vibrations and high-frequency noise. By applying the techniques of the present invention to a noisy vibration signal, an improved signal-to-noise ratio (SNR) can be obtained. This improved SNR makes spectral peaks of interest more easily distinguished from the noise floor, allowing earlier detection of a developing machine fault.
The Preferred Form System
Referring to Figure 1, a simplified processing chain of the preferred form system is shown generally as 100. The system 100, which is used for vibration signal processing, includes one or more vibration sensors 102 to generate electrical vibration signals representing 15 sensed mechanical vibrations. The sensor 102 is typically an accelerometer, but other sensors capable of sensing mechanical vibrations and representing them as electrical vibration signals could also be used.
Once the vibration signals are obtained, they are conditioned in a conditioning circuit 104. 20 In the preferred form, the conditioning circuit 104 includes a buffer, amplifier and anti-alias filter (i.e. low-pass filter). After being conditioned, the signals are high-pass filtered at 106. This results in only those signals above a user-selected lowest frequency (/)„,„) being allowed to pass through the filter 106. The value for D/mi/ is preferably selected via an input means, which may be one or more buttons, switches, knobs and the like, whether in 25 hardware form or software form.
Skilled persons will appreciate that it is not essential to provide the high-pass filter 106 separately from the conditioning circuit 104. Where necessary or desired, the high-pass filter 106 may form part of the conditioning circuit, for instance, by combining the high-30 pass filter 106 with the low-pass filter of the conditioning circuit 104 such that band-pass filter is provided.
1310770
8
The signals that are output by the high-pass filter 106 are then sampled by an analogue-to-digital converter (ADC) 108 to produce digital vibration samples. In the preferred form, the user selects via an input means the highest frequency of interest (DMflJ, which is used to select the appropriate Nyquist sampling rate of typically Dmax X 2.56 for the ADC.
Once digital vibration samples are obtained, the samples are sent to an averaging circuit 110. The averaging circuit 110 divides or breaks the samples into equal time intervals. The number of samples in the interval will be approximately the ratio of the user-selected highest frequency of interest and a user-selected maximum fault frequency of
interest (DFmJ:
Samples in Interval ~ D,mx / DFtmx (1)
Although not shown, the averaging circuit 110 includes a rectifier to rectify the samples in 15 each time interval. There is also an averager to determine the arithmetic mean or root mean square (RMS) of the rectified samples in each time interval. The determination of the arithmetic mean of the samples helps distinguish between signal peaks of interest and unwanted noise peaks. Once computed, the arithmetic mean values are saved into a waveform buffer 112. The buffer 112 collects the digital samples over a period and 20 outputs the collected samples as a time waveform 114. The time waveform 114, which is made up of only the averaged values, is then processed as if it were an independendy-detected vibration signal. For instance, the time waveform may be sent through a windowing and Fast Fourier Transform (FFT) process at 116 and an averaging process at 118 to produce a frequency spectrum of the vibration at 120.
Typically, the averaging circuit 110, waveform buffer 112, the processes of windowing and a Fast Fourier Transform (FFT) at 116 and averaging at 118 are performed using a microprocessor. To relieve the microprocessor of computations, the high-pass filter 106 is ideally implemented in analogue hardware. However, skilled persons will appreciate that 30 the high-pass filter could also be implemented as a digital filter by the microprocessor.
1310770
9
The Preferred Form Method
In one preferred form, the method of the present invention begins with the step of obtaining electrical vibration signals having amplitudes and frequencies that represent 5 mechanical vibrations of a machine. The signals are typically obtained directly from one or more sensors. In certain cases, however, the signals may be obtained from a storage device that stores signals that were sensed previously.
The signals are high pass filtered and then divided into time intervals, each with an interval 10 length set according to equation (1) noted earlier. For each interval, the method then determines the average absolute (rectified) amplitude of the digital vibration samples. The average amplitude of each interval is then used to generate a time waveform such that the waveform comprises the determined average amplitudes. This time waveform can then be processed as if it were an independently-detected signal.
In another form, the above preferred method may include the step of conditioning the electrical vibration signals that are obtained using vibration sensors. This step may comprise: (i) buffering the electrical vibration signals, (ii) adjusting the gain of the signals, and (iii) low-pass filtering the signals for anti-aliasing.
Example
The above method of the present invention was simulated and compared against a classic demodulation process and the prior art PeakVue process. The simulation was based on a 25 worst-case type of test signal, shown in Figure 2A. Specifically, the signal included a IV low-frequency noise at 500Hz, ldB white noise and a simulated 60Hz impact signal. The simulated impact signal was a IV, 5 kHz sinusoid passed through a Blackman-Harris window with a duration of 5ms.
The result of classic demodulation is shown in Figure 2B, the result of PeakVue processing is shown in Figure 2C and the result of the method of the present invention is shown in Figure 2D. Although all processes picked up the 60Hz impacting signal, the method of the
1310770
present invention provided the best SNR of ~15dB, compared to PeakVue's 12dB and classic demodulation's lOdB SNR.
The foregoing describes the invention including preferred forms thereof. Alterations and 5 modifications as will be obvious to those skilled in the art are intended to be incorporated within the scope hereof, as defined by the accompanying claims.
1310770
11
Claims (9)
1. A method of vibration signal processing comprising the steps of: obtaining electrical vibration signals having amplitudes and frequencies that represent mechanical vibrations of a machine; converting the electrical vibration signals into digital vibration samples; dividing the digital vibration samples into equal time intervals; determining the average absolute amplitude of the digital vibration samples for each time interval; generating a time waveform comprising the determined average amplitudes; and processing the time waveform as if it were an independently-detected signal.
2. A method of vibration signal processing comprising the steps of: measuring the vibration of a machine using one or more sensors to obtain electrical vibration signals; conditioning the electrical vibration signals; high-pass filtering the conditioned signals; converting the high-pass filtered electrical vibration signals into digital vibration samples; dividing the digital vibration samples into equal time intervals and, for each time interval, rectifying the digital vibration samples and determining the arithmetic mean of the rectified digital vibration samples; buffering the mean digital vibration samples to create a time waveform representing the mean digital vibration samples; and processing the time waveform as if it were an independently-detected signal.
3. The method as claimed in claim 2, wherein the step of conditioning the electrical vibration signals comprises: buffering the electrical vibration signals; adjusting the gain of the buffered electrical vibration signals; and low-pass filtering the gain-adjusted electrical vibration signals. 12
4. The method as claimed in claim 2 or 3, wherein the step of high-pass filtering includes the step of selecting a lowest frequency of interest, Dm/I.
5. The method as claimed in any one of the preceding claims, wherein the step of converting the electrical vibration signals into digital vibration samples includes the steps of: selecting the highest frequency of interest, Dmax, and converting the electrical vibration signals into digital vibration samples at a sampling rate of typically Dmi, x 2.56.
6. The method as claimed in any one of the preceding claims, wherein the step of dividing the digital vibration samples into equal time intervals includes the steps of: selecting a maximum fault frequency of interest, DFmax; dividing the digital vibration samples into equal time intervals, each having approximately D/mx / DFmax number of samples.
7. The method as claimed in any one of the preceding claims, wherein the step of processing the time waveform comprises the steps of: windowing the time waveform; transforming the windowed time waveform into a frequency spectrum; and determining an average of the frequency spectrum.
8. A system for vibration signal processing comprising: one or more vibration sensors configured to sense mechanical vibrations and to represent the sensed mechanical vibrations as electrical vibration signals; a conditioning circuit to condition the electrical vibration signals; a high-pass filter circuit; an analogue-to-digital converter (ADC) to convert the high-pass filtered electrical vibration signals into digital vibration samples; an averaging circuit configured to divide the digital vibration samples into equal time intervals, the averaging circuit including a rectifier to rectify the digital vibration samples in each time interval, and an averager to determine the arithmetic mean or root mean square (RMS) of the rectified digital vibration samples in each time interval; and 13 a buffer to receive the mean of the rectified digital vibration samples of each time interval, and to output the samples as a time waveform for further analysis. 5 implemented in a processor. 10. The system as claimed in claim 8 or 9, wherein the processor is also configured to: window the time waveform; transforming the windowed time waveform into a frequency spectrum; and 11. The system as claimed in any one of claims 8 to 10, wherein the conditioning circuit comprises a buffer, an amplifier, a low-pass filter and a high-pass filter. 15 12. The system as claimed in claim 11, wherein the system includes an input means for a user to select a lowest frequency of interest (Daif) for the high-pass filter. 13. The system as claimed in any one of claims 8 to 12, wherein the system includes an input means for a user to select the highest frequency of interest for use in setting 20 the sampling rate of the ADC at typically Dmtlx x 2.56. 14. The system as claimed in claim 13, wherein the system includes an input means for a user to select a maximum fault frequency of interest (DFmaJ for use in setting each of the time intervals with a length in samples of approximately Dmax / DFmax. 25 15. A method of vibration signal processing substantially as herein described with reference to Figure 1. 16. A system for vibration signal processing substantially as herein described with 30 reference to Figure 1.
9. The system as claimed in claim 8, wherein the averaging circuit and buffer are 10 determining an average of the frequency spectrum. COMMTEST INSTRUMENTS LIMITED By the authorized agents A]- Per: IPONZ 2 9 JAN 21108 1310770
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
NZ55404907A NZ554049A (en) | 2007-03-22 | 2007-03-22 | Method and system for vibration signal processing |
CN2008100041612A CN101271018B (en) | 2007-03-22 | 2008-01-23 | Method and system for vibration signal processing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
NZ55404907A NZ554049A (en) | 2007-03-22 | 2007-03-22 | Method and system for vibration signal processing |
Publications (1)
Publication Number | Publication Date |
---|---|
NZ554049A true NZ554049A (en) | 2008-03-28 |
Family
ID=39203985
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NZ55404907A NZ554049A (en) | 2007-03-22 | 2007-03-22 | Method and system for vibration signal processing |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN101271018B (en) |
NZ (1) | NZ554049A (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5421669B2 (en) * | 2009-06-23 | 2014-02-19 | アズビル株式会社 | Amplitude calculation apparatus and amplitude calculation method |
CN103398846B (en) * | 2013-08-16 | 2016-06-22 | 大连美恒时代科技有限公司 | A kind of decelerator health analysis method and analysis platform system thereof |
CN104568133B (en) * | 2014-12-23 | 2018-01-02 | 潍柴动力股份有限公司 | A kind of determination methods and system of motor vehicle abnormal vibrations |
DE102015121621B4 (en) * | 2015-12-11 | 2018-03-01 | Endress+Hauser Gmbh+Co. Kg | Device for the safe determination and / or monitoring of a process variable |
JP2017173076A (en) * | 2016-03-23 | 2017-09-28 | Ntn株式会社 | State monitoring system and wind generator system with the same |
US10788395B2 (en) * | 2017-02-10 | 2020-09-29 | Aktiebolaget Skf | Method and device of processing of vibration sensor signals |
JP7010031B2 (en) * | 2018-02-01 | 2022-01-26 | 株式会社島津製作所 | Amplitude detection method and material tester |
CN109032961B (en) * | 2018-07-11 | 2019-10-01 | 中国科学院地质与地球物理研究所 | A kind of underground vibrating impact data record method |
CN108918982B (en) * | 2018-07-26 | 2021-08-17 | 汉宇集团股份有限公司 | Method for acquiring signal peak value |
DE102021121071A1 (en) * | 2020-08-19 | 2022-02-24 | Computational Systems, Inc. | Eliminating the effects of an asymptotically decaying DC bias from a vibration waveform |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6321602B1 (en) * | 1999-09-28 | 2001-11-27 | Rockwell Science Center, Llc | Condition based monitoring by vibrational analysis |
US6668234B2 (en) * | 2002-03-22 | 2003-12-23 | Abb Inc. | Method and apparatus for calculating the amplitude of a complex waveform associated with a rotating machine shaft after removing the running speed frequency |
CN2824176Y (en) * | 2005-05-23 | 2006-10-04 | 湖南科技大学 | High-integrated portable vibration signal treater |
-
2007
- 2007-03-22 NZ NZ55404907A patent/NZ554049A/en not_active IP Right Cessation
-
2008
- 2008-01-23 CN CN2008100041612A patent/CN101271018B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101271018B (en) | 2010-06-02 |
CN101271018A (en) | 2008-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7493220B2 (en) | Method and system for vibration signal processing | |
NZ554049A (en) | Method and system for vibration signal processing | |
CA2530344C (en) | Improvements in or relating to vibration analysis | |
JP2008525808A (en) | Method for identifying impact location | |
EP2495996A3 (en) | Method for measuring critical gain on a hearing aid | |
US20120157870A1 (en) | Noise reduction of breathing signals | |
US7713207B2 (en) | Ultrasonic apparatus and method for generating audio signals using doppler | |
US6684700B1 (en) | Stress wave sensor | |
CA2582560C (en) | Method and system for vibration signal processing | |
JP2013014286A (en) | Vehicle approach notification device and program | |
KR20170124279A (en) | Method and Apparatus for DEMON Processing in order that Removal of External Target Noise When measuring underwater radiated noise | |
CN118243214A (en) | Test method for measuring audible tone in environmental noise | |
JPWO2019049667A1 (en) | Heart rate detector, heart rate detection method and program | |
Guski et al. | Impulsive noise detection in sweep measurements | |
US20130096855A1 (en) | Detection of an Abnormal Signal in a Compound Sampled Signal | |
EP1769267B1 (en) | Movement detection system and method | |
WO2016207951A1 (en) | Shunt murmur analysis device, shunt murmur analysis method, computer program, and recording medium | |
CN104034970B (en) | A kind of sound level measurement method based on computer sound card | |
CN109682958B (en) | Acceleration sensor signal compensation method for thrombelastogram instrument | |
JPWO2020090763A1 (en) | Processing equipment, systems, processing methods, and programs | |
JP2019017406A (en) | Noise removal device of kinetic recall output signal of time series data, noise removal method, program, analysis device of respiration, heart beat and voice by cluster analysis, analysis method and program | |
JP7122225B2 (en) | Processing device, system, processing method, and program | |
JP2849711B2 (en) | Accelerometer | |
WO2016207950A1 (en) | Shunt murmur analysis device, shunt murmur analysis method, computer program, and recording medium | |
JP2019111347A (en) | Shunt sound analysis device, shunt sound analysis method, computer program and recording medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PSEA | Patent sealed | ||
RENW | Renewal (renewal fees accepted) | ||
ASS | Change of ownership |
Owner name: GENERAL ELECTRIC COMPANY, US Effective date: 20130627 |
|
RENW | Renewal (renewal fees accepted) |
Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 29 JAN 2016 BY COMPUTER PACKAGES INC Effective date: 20141231 |
|
LAPS | Patent lapsed |