CN101451920B - All-band characteristic prewarning method for mechanical vibration signal - Google Patents

All-band characteristic prewarning method for mechanical vibration signal Download PDF

Info

Publication number
CN101451920B
CN101451920B CN2007101719820A CN200710171982A CN101451920B CN 101451920 B CN101451920 B CN 101451920B CN 2007101719820 A CN2007101719820 A CN 2007101719820A CN 200710171982 A CN200710171982 A CN 200710171982A CN 101451920 B CN101451920 B CN 101451920B
Authority
CN
China
Prior art keywords
vibration
amplitude
centerdot
value
band characteristic
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.)
Expired - Fee Related
Application number
CN2007101719820A
Other languages
Chinese (zh)
Other versions
CN101451920A (en
Inventor
刘晗
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Baosteel Industry Technological Service Co Ltd
Shanghai Baosteel Industry Inspection Corp
Original Assignee
Shanghai Baosteel Industry Inspection Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Baosteel Industry Inspection Corp filed Critical Shanghai Baosteel Industry Inspection Corp
Priority to CN2007101719820A priority Critical patent/CN101451920B/en
Publication of CN101451920A publication Critical patent/CN101451920A/en
Application granted granted Critical
Publication of CN101451920B publication Critical patent/CN101451920B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses an early-warning method for mechanical vibration signal full frequency band characteristic, including acquisition multiple sets of the mechanical vibration signals, calculation of components of each weight of the vibration signals, and obtainment of a plurality of sets of vibration amplitudes and frequency value sequence sets, namely, for specific fi frequency, a plurality of sets of vibration amplitudes can be obtained, an amplitude early-warning value for the specific fi frequency can be obtained by calculating average value, deviation value of the plurality of sets of vibration amplitudes; for the plurality of sets of fi frequencies, an early-warning line for the mechanical vibration signal full frequency band characteristic is drawn on an original vibration signal spectrogram by manufacturing the spectrogram, in order to determine early warning for the mechanical vibration signal full frequency band characteristic. Hence, the method can realize early warning for the mechanical vibration signal full frequency band characteristic even when certain frequency component vibration amplitude exceeds standard to early warn.

Description

All-band characteristic prewarning method for mechanical vibration signal
Technical field
The present invention relates to the method for machine performance monitoring and diagnosis, relate in particular to a kind of method of mechanical oscillation signal all-band characteristic prewarning.
Background technology
At mechanical field, 60%~70% mechanical fault is arranged according to statistics all is by vibration and have the noise of vibration radiation to reflect, therefore, for other plant equipment Condition Monitoring Technology, whether in allowed limits and be in which rank vibration monitoring is simple and easy to do with diagnosis, normal adopt the actual measurement vibration level and differentiate this vibration.
Judge whether plant equipment exists the degree of fault and fault, must measured value and the specified standard value that characterize machine state be compared.A large amount of example explanations exists special association between the mechanical fault of vibration value and machine, and vibration values does not show that greatly fault must occur.Will inevitably in vibration, reflect to some extent when on the one hand, machine breaks down; On the other hand, vibration may not necessarily be broken down again to the machine of certain value greatly.In view of the situation, the former becomes the benchmark of vibration standard, and the latter becomes the difficult point of standard formulation again.
In ISO mechanical vibration standard, adopt vibration signal entire spectrum global vibration value OA to characterize the vibration of machine mostly,
OA = Σ i = 1 n A i 2 N BF = A 1 2 + A 2 2 + · · · + A n 2 1.5
OA=vibration signal full range global vibration value n=FFT resolution spectral line number wherein
A i=every FFT spectral line amplitude N BF=Hanning window function grass parameter
And stipulated the vibration early warning threshold value of all kinds of machines, its absolute vibration, axle of only having stipulated to estimate axle relatively and the vibration of structure member and the vibration of bearing seat and is only assessed the wide-band vibration value, does not have the evaluation criterion of frequency component and phase place; Its interpretational criteria is the amplitude of wide-band vibration, and another is the variation of considering amplitude, increases or reduces no matter be.
In engineering practice was used, the method for early warning of above-mentioned standard code had been obtained big effect, but has also been had certain deficiency in machine diagnostic, mainly contain:
Although certain frequency content vibration amplitude may exceed tolerance interval in 1 vibration signal, the global vibration value does not exceed standard, and early warning in time causes and fails to report, and is insensitive to initial failure.
Even 2 global vibrations exceed standard, can't know that also those concrete frequency content component amplitudes exceed standard, the disappearance diagnostic message causes accurately fault location position and reason.
3, former early warning belongs to the early warning of entire spectrum global vibration monodrome, can't realize vibration signal full range band arbitrfary point or the numerous points early warning is set.
Summary of the invention
Technical matters to be solved by this invention provides a kind of all-band characteristic prewarning method for mechanical vibration signal, and it can carry out the timely early warning that the global vibration value exceeds standard, and also can realize the timely early warning that certain frequency content vibration amplitude exceeds standard.
For solving the problems of the technologies described above, all-band characteristic prewarning method for mechanical vibration signal of the present invention comprises the steps,
Step 1, vibration signal pretreatment, the mechanical oscillation signal of collecting device, the mechanical vibration discrete signal of establishing collection is: x (n Δ t), n=1 wherein, 2....N, Δ t are sampling time interval, by the Fast Fourier Transform (FFT) formula S ( f ) = Δt Σ n = - ∞ ∞ x ( nΔt ) e - j 2 πfΔt
Calculate each frequency amplitude of vibration signal, and obtain
Figure S2007101719820D00022
Group (A i, f i) vibration amplitude and frequency values;
Step 2, rumble spectrum data processing to an equipment, are gathered the n group vibration signal of its running status, form the n group length by said method and are
Figure S2007101719820D00023
(A i, f i) vibration amplitude and frequency values sequence sets, promptly for specific f iFrequency can obtain n group A I1, A I2, A I3... .A InVibration amplitude;
Step 3, employing statistical method draw f iThe amplitude early warning value of frequency,
Calculate vibration amplitude mean value
Figure S2007101719820D00024
M ‾ i = A il + A i 2 + A i 3 + A i 4 + · · · · · · + A in n
Calculate vibration amplitude and mean value Standard deviation δ i:
δ i = ( Ai i - M ‾ i ) 2 + ( A i 2 - M ‾ i ) 2 + · · · + ( A in - M ‾ i ) 2 n - 1
Determine to judge f iThe amplitude early warning value M of frequency Ia:
M ia = M ‾ i + ( 1 ~ 2 ) δ i
Step 4, mechanical oscillation signal all-band characteristic prewarning,
Method according to step 3 is right
Figure S2007101719820D00032
Group (f i) all just calculate M separately Ia, and make spectrogram, on the original vibration signal spectrogram, draw out mechanical oscillation signal all-band characteristic prewarning line;
To the machine vibration signal of new collection, calculate separately frequency values and respective magnitudes by step 1, promptly Group (A i, f i) vibration amplitude and frequency values, the early warning value (M that its value calculates with step 3 Ia, f i) relatively, exceed the all-band characteristic prewarning line, realize early warning.
Because all-band characteristic prewarning method for mechanical vibration signal of the present invention has adopted technique scheme, it is the mechanical oscillation signal of collecting device, and this mechanical oscillation signal carried out particular processing, calculate each frequency amplitude of vibration signal, form the frequency values sequence sets, calculate for specific f iThe amplitude early warning value time, also can draw the mechanical oscillation signal all-band characteristic prewarning under the frequency by the making of spectrogram.
Description of drawings
The present invention is described in further detail below in conjunction with drawings and embodiments:
Fig. 1 is the spectrogram of each frequency amplitude of vibration signal of the present invention,
Fig. 2 is an all-band characteristic prewarning spectrogram of the present invention.
Embodiment
All-band characteristic prewarning method for mechanical vibration signal of the present invention comprises the steps,
Step 1, vibration signal pretreatment, the mechanical oscillation signal of collecting device, the mechanical vibration discrete signal of establishing collection is: x (n Δ t), n=1 wherein, 2....N, Δ t are sampling time interval, by the Fast Fourier Transform (FFT) formula S ( f ) = Δt Σ n = - ∞ ∞ x ( nΔt ) e - j 2 πfΔt
As shown in Figure 1, calculate each frequency amplitude of vibration signal, and obtain
Figure S2007101719820D00035
Group (A i, f i) vibration amplitude and frequency values;
Step 2, rumble spectrum data processing to an equipment, are gathered the n group vibration signal of its running status, form the n group length by said method and are
Figure S2007101719820D00036
(A i, f i) vibration amplitude and frequency values sequence sets, promptly for specific f iFrequency can obtain n group A I1, A I2, A I3... .A InVibration amplitude;
Step 3, employing statistical method draw f iThe amplitude early warning value of frequency,
Calculate vibration amplitude mean value
Figure S2007101719820D00041
M ‾ i = A i 1 + A i 2 + A i 3 + A i 4 + · · · · · · + A in n
Calculate vibration amplitude and mean value
Figure S2007101719820D00043
Standard deviation δ i:
δ i = ( Ai 1 - M ‾ i ) 2 + ( A i 2 - M ‾ i ) 2 + · · · + ( A in - M ‾ i ) 2 n - 1
Determine to judge f iThe amplitude early warning value M of frequency Ia:
M ia = M ‾ i + ( 1 ~ 2 ) δ i
Step 4, mechanical oscillation signal all-band characteristic prewarning,
Method according to step 3 is right
Figure S2007101719820D00046
Group (f i) all just calculate M separately Ia, as shown in Figure 2, make spectrogram, on the original vibration signal spectrogram, draw out mechanical oscillation signal all-band characteristic prewarning line.
To the machine vibration signal of new collection, calculate separately frequency values and respective magnitudes by step 1, promptly
Figure S2007101719820D00047
Group (A i, f i) vibration amplitude and frequency values, the early warning value (M that its value calculates with step 3 Ia, f i) relatively, exceed the all-band characteristic prewarning line, realize early warning.
In accompanying drawing 1 and the accompanying drawing 2, X-axis is frequency (Hz) or rotating speed (RPM), and its span is [0,1/ (2 Δ t)], and Δ t is a sampling interval; Y-axis is displacement (um) or speed (mm/s) or acceleration (mm/s 2), also can be its db decibel value.

Claims (1)

1. all-band characteristic prewarning method for mechanical vibration signal is characterized in that: this method comprises the steps,
Step 1, vibration signal pretreatment, the mechanical oscillation signal of collecting device, the mechanical vibration discrete signal of establishing collection is: x (n Δ t), n=1 wherein, 2....N, Δ t is a sampling time interval, selects suitable window function, by the Fast Fourier Transform (FFT) formula
Figure FSB00000325415800011
Calculate each frequency amplitude of vibration signal, and obtain
Figure FSB00000325415800012
Group vibration amplitude A iWith frequency values f i
Step 2, rumble spectrum data processing to an equipment, are gathered the n group vibration signal of its running status, form the n group length by first step and are
Figure FSB00000325415800013
Vibration amplitude A iWith frequency values f iSequence sets is promptly for specific frequency values f i, can obtain n group A I1, A I2, A I3... .A InVibration amplitude;
Step 3, employing statistical method draw frequency values f iThe amplitude early warning value,
Calculate vibration amplitude mean value
M i ‾ = A i 1 + A i 2 + A i 3 + A i 4 + · · · · · · + A in n
Calculate vibration amplitude and vibration amplitude mean value
Figure FSB00000325415800016
Standard deviation δ i:
δ i = ( Ai 1 - M i ‾ ) 2 + ( A i 2 - M i ‾ ) 2 + · · · + ( A in - M i ‾ ) 2 n - 1
Determine determination frequency value f iAmplitude early warning value M Ia:
M ia = M i ‾ + ( 1 ~ 2 ) δ i
Step 4, mechanical oscillation signal all-band characteristic prewarning,
Method according to step 3 is right
Figure FSB00000325415800019
Class frequency value f iAll just calculate amplitude early warning value M separately Ia, and make spectrogram, on the original vibration signal spectrogram, draw out mechanical oscillation signal all-band characteristic prewarning line;
To the machine vibration signal of new collection, calculate separately frequency values and corresponding vibration amplitude by step 1, promptly
Figure FSB000003254158000110
Group vibration amplitude A iWith frequency values f i, the amplitude early warning value M that its value calculates with step 3 IaWith frequency values f iRelatively, exceed the all-band characteristic prewarning line, realize early warning.
CN2007101719820A 2007-12-07 2007-12-07 All-band characteristic prewarning method for mechanical vibration signal Expired - Fee Related CN101451920B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2007101719820A CN101451920B (en) 2007-12-07 2007-12-07 All-band characteristic prewarning method for mechanical vibration signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2007101719820A CN101451920B (en) 2007-12-07 2007-12-07 All-band characteristic prewarning method for mechanical vibration signal

Publications (2)

Publication Number Publication Date
CN101451920A CN101451920A (en) 2009-06-10
CN101451920B true CN101451920B (en) 2011-05-11

Family

ID=40734287

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007101719820A Expired - Fee Related CN101451920B (en) 2007-12-07 2007-12-07 All-band characteristic prewarning method for mechanical vibration signal

Country Status (1)

Country Link
CN (1) CN101451920B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5421669B2 (en) * 2009-06-23 2014-02-19 アズビル株式会社 Amplitude calculation apparatus and amplitude calculation method
CN102032942B (en) * 2009-09-25 2012-07-25 上海宝钢工业检测公司 Equipment vibration signal full frequency band characteristic extraction and alarm method based on Gaussian smoothing
EP2547991A4 (en) * 2010-03-19 2017-10-04 ABB Oy Health monitoring method and system for drives
CN102890750B (en) * 2012-09-11 2015-06-10 中国包装科研测试中心 Data analysis method for transportation safety recorder
CN105651376B (en) * 2014-11-10 2019-08-06 上海宝钢工业技术服务有限公司 The analysis of mechanical equipment off-line checking system vibration signals spectrograph and alarm method
CN106249146B (en) * 2016-08-30 2019-03-22 河南中烟工业有限责任公司洛阳卷烟厂 A kind of analysis of cigar mill's electric motor operation state and method for early warning
CN107843334A (en) * 2017-12-07 2018-03-27 广东省水利水电科学研究院 The intelligent monitor system that a kind of soft base sluice plate comes to nothing

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5305235A (en) * 1991-07-10 1994-04-19 Mitsubishi Denki Kabushiki Kaisha Monitoring diagnosis device for electrical appliance
CN1370983A (en) * 2001-02-14 2002-09-25 深圳市创为实测控技术有限公司 Sensitive rotating machines monitoring and data storing method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5305235A (en) * 1991-07-10 1994-04-19 Mitsubishi Denki Kabushiki Kaisha Monitoring diagnosis device for electrical appliance
CN1370983A (en) * 2001-02-14 2002-09-25 深圳市创为实测控技术有限公司 Sensitive rotating machines monitoring and data storing method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JP特开2003-343338A 2003.12.03

Also Published As

Publication number Publication date
CN101451920A (en) 2009-06-10

Similar Documents

Publication Publication Date Title
CN101451920B (en) All-band characteristic prewarning method for mechanical vibration signal
CN102788695B (en) Identification method of rolling bearing abrasion
EP3043062B1 (en) Condition monitoring apparatus for wind turbine
CN104075795A (en) Method and system for monitoring vibration state of impeller of wind generating set
CN101363420B (en) Event monitoring via combination of signals
CN101592590B (en) Fault indirect diagnosis technique of rotating blade
US7409854B2 (en) Method and apparatus for determining an operating status of a turbine engine
CN101265878A (en) Tonal emission control for wind turbines
CN104865400B (en) A kind of detection recognition method and system of Wind turbines rotating speed
CN104535323A (en) Locomotive wheelset bearing fault diagnosis method based on angular domain-time domain-frequency domain
US20120010852A1 (en) Method for monitoring wind turbines
CN101430247A (en) Real-time diagnosis method for random vibration fault of steam generator set
US9476860B2 (en) Component failure detection system
CN105756864A (en) Blade imbalance fault diagnosis method based on stator current of double-fed wind generating set
CN109708884B (en) Universal shaft fault detection method and equipment
CN102242695B (en) Mutative-peak-index-based wind generating set abnormal condition early warning method
CN111594395B (en) Wind turbine generator tower resonance identification method and device and monitoring alarm system
CN102966526B (en) Based on the compressor low frequency flow induced vibrations class method for diagnosing faults that axial vibration is analyzed
CN102778354A (en) Method for diagnosing failure of wind-powered rotary support based on wavelet analysis
US20230272783A1 (en) Device for identifying an accumulation of ice on rotor blades of a wind turbine and method for teaching such a device
Shan et al. A novel experimental research on vibration characteristics of the running high-speed motorized spindles
CN111795824A (en) Fault diagnosis method and system for mechanical power system
CN110082101A (en) Planetary gear system fault monitoring method based on input and output torque power characteristic
CN115111181A (en) Blade health degree intelligent monitoring, analyzing and evaluating method, system device and application of axial flow fan
CN114278518A (en) Tower drum bolt state monitoring method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C56 Change in the name or address of the patentee

Owner name: SHANGHAI BAOSTEEL INDUSTRIAL INSPECTION CO., LTD.

Free format text: FORMER NAME: SHANGHAI BAOSTEEL INDUSTRY INSPECTION CORP.

Owner name: SHANGHAI BAOSTEEL INDUSTRY TECHNOLOGICAL SERVICE C

Free format text: FORMER NAME: SHANGHAI BAOSTEEL INDUSTRIAL INSPECTION CO., LTD.

CP01 Change in the name or title of a patent holder

Address after: 201900 Shanghai city Baoshan District Meipu Road No. 335

Patentee after: Shanghai Baosteel Industrial Inspection Co.,Ltd.

Address before: 201900 Shanghai city Baoshan District Meipu Road No. 335

Patentee before: Shanghai Baosteel Industry Inspection Corp.

Address after: 201900 Shanghai city Baoshan District Meipu Road No. 335

Patentee after: SHANGHAI BAOSTEEL INDUSTRY TECHNOLOGICAL SERVICE Co.,Ltd.

Address before: 201900 Shanghai city Baoshan District Meipu Road No. 335

Patentee before: Shanghai Baosteel Industrial Inspection Co.,Ltd.

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110511

Termination date: 20211207