CN103076174A - Bearing fault and fault degree diagnostic method on basis of phase-locked amplification algorithm - Google Patents
Bearing fault and fault degree diagnostic method on basis of phase-locked amplification algorithm Download PDFInfo
- Publication number
- CN103076174A CN103076174A CN2013100003895A CN201310000389A CN103076174A CN 103076174 A CN103076174 A CN 103076174A CN 2013100003895 A CN2013100003895 A CN 2013100003895A CN 201310000389 A CN201310000389 A CN 201310000389A CN 103076174 A CN103076174 A CN 103076174A
- Authority
- CN
- China
- Prior art keywords
- fault
- signal
- bearing
- max
- phase
- 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
- 230000003321 amplification Effects 0.000 title claims abstract description 35
- 238000003199 nucleic acid amplification method Methods 0.000 title claims abstract description 35
- 238000002405 diagnostic procedure Methods 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 17
- 239000000203 mixture Substances 0.000 claims description 16
- 238000013016 damping Methods 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 8
- 230000035939 shock Effects 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 7
- 101100129500 Caenorhabditis elegans max-2 gene Proteins 0.000 claims description 6
- 230000002950 deficient Effects 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 4
- 238000005096 rolling process Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 239000000470 constituent Substances 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims description 2
- 238000003745 diagnosis Methods 0.000 abstract description 23
- 238000000926 separation method Methods 0.000 abstract description 4
- 230000005540 biological transmission Effects 0.000 abstract 1
- 239000004615 ingredient Substances 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000000975 dye Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 240000004859 Gamochaeta purpurea Species 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a bearing fault and fault degree diagnostic method on the basis of a phase-locked amplification algorithm. According to the invention, due to the advantages of stable center frequency, narrow transmission band, high quality factor and the like of a phase-locked amplifier, a phase-locked amplification principle is applied to diagnosis on a bearing fault and a fault degree by utilizing the phase-locked amplifier. A reference signal on the basis of the bearing fault signal characteristics is established and the phase-locked amplification principle is utilized to remain and amplify information related to the fault and reduce noise signals which are not related to fault ingredients by multiplication and integral operation, so that rapid separation of signal noise and positioning on the fault are implemented. Meanwhile, a reference signal added with fault degree parameters is established and diagnosis on the bearing fault degree is implemented by combining a kurtosis value.
Description
Technical field
The present invention relates to a kind of bearing fault and fault degree diagnostic method, particularly a kind of bearing fault and fault degree diagnostic method based on phase-locked interpolator arithmetic.
Background technology
The serviceability of industry kinematic train directly affects the work efficiency of whole production line.Bearing is as the important component part of industrial kinematic train, and its effect is very crucial.The inner topmost part of bearing comprises rolling bearing and gear, moves under high speed, the fully loaded transportation condition because these two kinds of parts are everlasting, and rate of breakdown is higher, will cause serious casualties and economic loss in case break down.Therefore bearing is carried out condition monitoring and fault diagnosis to guaranteeing production safety, prevent major accident, important in inhibiting reduces production costs.
At present, adopting the vibration signal processing technology based on high-speed computer is the main flow of bearing failure diagnosis technology.The key of bearing failure diagnosis is how to extract fault signature from the fault vibration signal.But because the environmental facies of the operation of the actual centre bearer of engineering are worked as badly, its vibration signal is very complicated, contains much noise and labile factor, is a kind of typical non-stationary signal, and when particularly early defect appearred in bearing inside, the signal fault feature was very faint.And the equipment to diagnosing malfunction in engineering reality often only possesses some simple signal analysis abilities, for the obvious situation of fault certain effect is arranged still, feeble signal or the larger signal of signal to noise ratio (S/N ratio) often be can not get effectively analysis result, analytical effect to fault degree also has much room for improvement, and has affected the efficient of signal analysis.
Therefore how adopting effective analysis tool and algorithm, the fault of bearing and fault degree are analyzed and diagnosed, extract fault signature and realize fault progression status monitoring and diagnosis, is a large difficult point of bearing being carried out fault detection and diagnosis.
Lock-in amplifier is owing to have CFS center frequency stabilization, and passband is narrow, and the quality factor advantages of higher is having good performance aspect feeble signal amplification and the noise separation.
Yet current phaselock technique great majority are applied to the amplification of sine, cosine signal, and the contrast signal kind of its input is also very single.By to the analysis and modeling of bearing fault signal as can be known, the bearing fault signal has the feature of periodic signal and can utilize the principle of phase locking unit to carry out the signal processing, from its mathematic(al) representation, can find out, as long as seek based on the contrast signal of bearing fault characteristics and process accordingly, can realize theoretically the amplification of fault-signal and noise signal to zero.
Therefore the principle design algorithm that utilizes lock-in amplifier based on bearing fault and the fault degree diagnostic method of phase-locked interpolator arithmetic described herein causes the state of the bearing fault of significant impact to study mainly for meeting to production.After the feature of fully having analyzed the bearing fault signal, determine the input signal of lock-in amplifier, after multiplying and integral operation, relevant with fault information is stayed and amplify, those and the incoherent noise signal of trouble unit are extremely zero, realize the quick separation of noise and the location of fault.In input signal, add simultaneously the fault degree parameter, realized the diagnosis of bearing fault degree.
Summary of the invention
Above-mentioned technical matters in bearing failure diagnosis and the fault degree diagnosis the invention provides a kind of bearing fault based on phase-locked interpolator arithmetic and fault degree diagnostic method.
The technical scheme that the present invention solves the problems of the technologies described above comprises sets up the input reference signal model, carries out phase-locked amplification and calculates, and calculates fault degree.
Because the bearing fault characteristics signal also has periodic characteristics, therefore can consider to select suitable reference signal to utilize phase-locking and amplification principle to weaken noise signal in the bearing burst, amplify and extract wherein fault characteristic signals.
By analyzing bearing fault characteristics as can be known, bearing fault signal v
s(t) can be expressed as
That is:
k
1For the bearing fault signal amplitude can directly be measured, p wherein
1For the damping coefficient of bearing fault signal shock response can be according to bearing parameter measurements and calculations, f
1Corresponding to the damped natural frequency of measured bearing be each impact of bearing fault signal frequency can the corresponding frequency of spectrum peak position directly obtains T in the bearing fault signal spectrum figure resonant belt by reading
0For the bearing fault characteristics cycle can calculate by bearing parameter, all t are the time and can directly measure in the literary composition, all δ are the impulse function identifier in the literary composition, to be the maximum occurrences of time t be that the length of measured signal s (t) can directly be measured to all max (t) in the literary composition, in the literary composition all e be exponential function represent the symbol, all n are bearing rotary week number in the literary composition, represent with nonnegative integer.
Measured signal can be expressed as:
s(t)=v
s(t)+n(t)
Wherein n (t) is random noise, v
s(t) be the bearing fault signal
According to phase-locking and amplification principle, if select the signal close with tested composition and choose suitable frequency and can amplify and de-noising measured signal.Traditional lock-in amplifier is all selected sinusoidal signal and is adjusted cosine composition in the suitable frequency abstraction analyzed signal.For the bearing fault signal, its most basic effective constituent is the impact-attenuating signal, therefore can select the impact-attenuating signal as reference signal extraction and amplify effective trouble unit.
Establishing reference signal at this is:
Wherein y (t) is reference signal, k
2For the reference signal amplitude can be set as arbitrarily number greater than 1, the larger enlargement factor of value is larger, p
2Can get the arbitrary value value for the damping coefficient of reference signal shock response in 1000-1200, the larger kurtosis of value is larger, f
2For the measured bearing damped natural frequency that dopes is that reference signal frequency needs according to f
1Value be defined as equaling f
1Value.
Measured signal and contrast signal then through the output signal of multiplier are:
v(t)=s(t)+y(t)
Therefore signal v (t) can be expressed as through the output signal V of integrator again:
Following formula is phase-locked amplification expression formula, and wherein u is the initial position of reference signal, and max (t) is that the length of measured signal s (t) can directly be measured for the maximum occurrences of time t.
Following formula launches to be expressed as:
Order
Then V can be expressed as:
V=V
s+V
n
It is the addition of fault-signal item and noise signal item.V wherein
sBe integration fault-signal item later, V
nBe integration noise signal item later.
Wherein
Because f among the reference signal y (t)
2For the measured bearing damped natural frequency that dopes is that reference signal frequency needs according to f
1Value be defined as equaling f
1Value, i.e. f
1=f
2, so noise item V
nGet zero, i.e. V
n=0 V=V
s, so the signal V of integrator output equals:
Compare with measured signal s (t), dispelled noise contribution, and amplified the amplification of trouble unit.Namely realized the function of the amplification of denoising and fault-signal, wherein K is enlargement factor K=k
1k
2
Above-mentioned multiply each other and the process of integration is that phase-locked amplification process can be expressed as:
Be the sequence of impacts signal of a string not Noise, and fault signature has obtained K and has doubly amplified.Hence one can see that, when obtaining reference signal according to y (t) expression formula, the frequency values f in the reference signal
2The centre frequency of analyzed signal s (t) resonant frequency is f
2=f
1The time, can to utilize phase-locking and amplification principle be phase-locked amplification expression formula amplifies and carries out effective de-noising to each failure impact signal of bearing.
In order to diagnose the bearing fault size, can utilize bearing fault signal accurate model as the reference signal, introduced fault size parameter d in the model
xAt first need calculate the linear velocity of rolling body motion and the pulse width that different faults causes in order to obtain this model
Bearing roller linear velocity s:
S=πdf
r
Pulse width p
x:
The pulse x (t) that can obtain thus the defective generation can be expressed as:
The reference signal that the impact that is produced by defective namely contains the fault degree parameter can be expressed as:
φ′
imp(p
2,u,f
2,dx,d,fr)=conv(x(t),φ
imp(p
2,u,f
2))
φ wherein
Imp(p
2, u, f
2) for its expression formula of decaying exponential function be:
Be φ '
Imp(p
2, u, f
2, dx, d, fr) and can be expressed as φ
Imp(p
2, u, f
2) with the convolution of x (t), (conv is the convolution algorithm symbol).
Wherein d is that the bearing path can be determined f according to bearing designation
rFrequently can record by the sensor of special measurement rotating speed for turning.d
xFor fault diameter (unit: mm), p
2Can get the arbitrary value value for the damping coefficient of reference signal shock response in 1000-1200, the larger kurtosis of value is larger, f
2For the measured bearing damped natural frequency that dopes is that reference signal frequency needs according to f
1Value be defined as equaling f
1Value.
Parameter p in the model
2, u, f
2Value is identical with the value of using y (t) as with reference to signal the time.Make d in the model
xSpan is [0.5, max (d
x)], get a value every 0.5, get altogether
Individual d
xValue, wherein max (d
x) be maximum fault diameter, this
Individual amount is correspondence 1,2 from small to large,,
This
Individual integer as the fault level label, is NO with letter representation.I.e. 0.5 corresponding numeral 1,1 corresponding numeral 2,, 5 corresponding numerals 10,, max (d
x) corresponding numeral
Be 0.5 each standards change d of 0.5 that increases according to initial value
xValue, d of every variation
xValue is utilized the phase-locked simplification expression formula of putting that signal is carried out single treatment, and is asked for the kurtosis index after the processing, and the record kurtosis is with d
xChange curve.
Kurtosis ks can be expressed as:
Wherein V is the analyzed signal after phase-locked amplification is processed, and mean is mean value.
What try to achieve
Choose three maximum kurtosis values among the individual kurtosis, record their corresponding integer fault level labels, multiply by corresponding failure grade label with kurtosis value after sum obtain the fault degree eigenwert divided by three kurtosis sums.
Technique effect of the present invention is: utilize lock-in amplifier owing to have CFS center frequency stabilization, passband is narrow, and the quality factor advantages of higher is applied to phase-locking and amplification principle the diagnosis of bearing fault and fault degree.The contrast signal of having set up based on bearing fault characteristics utilizes phase-locking and amplification principle, after multiplying and integral operation, relevant with fault information is stayed and amplify, those and the incoherent noise signal of trouble unit to zero, are realized the quick separation of noise and the location of fault.Set up simultaneously a kind of reference signal that adds the fault degree parameter, realized the diagnosis of bearing fault degree in conjunction with the kurtosis index.
Description of drawings
The invention will be further described below in conjunction with the drawings and specific embodiments.
Fig. 1 is the chief component figure based on phase-locked interpolator arithmetic bearing failure diagnosis of the present invention.
Fig. 2 is of the present invention based on phase-locked interpolator arithmetic bearing failure diagnosis process flow diagram.
Fig. 3 is the emulation time domain waveform of housing washer fault.
Result after Fig. 4 processes through phase-locked interpolator arithmetic.
Fig. 5 is the process flow diagram of diagnosing based on phase-locked interpolator arithmetic bearing fault degree among the present invention.
Fig. 6 utilizes fault degree and the fault degree eigenwert graph of a relation that obtains based on phase-locked interpolator arithmetic bearing fault degree diagnosis among the present invention.
Embodiment
Fig. 1 is the chief component figure based on phase-locked interpolator arithmetic bearing failure diagnosis of the present invention.Formed by bearing failure diagnosis and bearing fault degree diagnosis two parts based on phase-locked interpolator arithmetic bearing failure diagnosis
Fig. 2 is the process flow diagram based on phase-locked interpolator arithmetic bearing failure diagnosis of the present invention.Below in conjunction with process flow diagram to being elaborated based on phase-locked interpolator arithmetic Method for Bearing Fault Diagnosis.
(1) utilize the acceleration vibration transducer that the vibration signal in the bearing operational process is measured, obtain vibration acceleration signal as signal s to be analyzed (t), sampling length is decided to be 2 integer power, sets sample frequency according to bearing rotating speed and model;
(2) set up reference signal
P wherein
2Be the vibration damping coefficient, according to the p as can be known of the reasoning in the summary of the invention
2Value do not affect amplification effect, any value in theory is for so that the larger p of the kurtosis of phase-locked amplified signal
2Can be taken as the arbitrary value among the 1000-1200, the larger kurtosis of value is larger.f
2For being reference signal frequency for the measured bearing damped natural frequency that dopes, can from the frequency spectrum of analyzed signal, read, namely get in the frequency spectrum resonant belt the corresponding frequency values in spectrum peak and namely get f
2=f
1k
2The value that is taken as greater than 1 is the reference signal amplitude, as long as can obtain amplification effect, k
2The larger enlargement factor of value larger, amplification effect is more obvious.
(3) utilize the phase-locked amplification expression formula after simplifying to carry out phase-locked amplification calculating, expression formula is:
Wherein u is the initial position of reference signal, and max (t) can directly measure for the maximum occurrences of time t is the length of measured signal s (t).
Fig. 3 is the emulation time domain waveform of the housing washer fault of adding noise.Signal length is 8192 points, and the Signal-to-Noise SNR that dyes after making an uproar is-9.87dB.As can be seen from the figure, the impact composition that dyes the rear signal of making an uproar has been submerged substantially.
Fig. 4 is the result after Fig. 3 signal is processed through phase-locked interpolator arithmetic, and as seen from the figure, fault signature composition has after treatment obtained amplification, and noise contribution has obtained inhibition.
Fig. 5 is the process flow diagram based on phase-locked interpolator arithmetic bearing fault degree diagnosis of the present invention.Below in conjunction with process flow diagram to being elaborated based on phase-locked interpolator arithmetic bearing fault degree diagnostic method.
(1) utilize the acceleration vibration transducer that the vibration signal in the bearing operational process is measured, obtain vibration acceleration signal as signal s to be analyzed (t), sampling length is decided to be 2 integer power, sets sample frequency (use based on the signal that collects in the phase-locked interpolator arithmetic bearing failure diagnosis and get final product) here according to bearing rotating speed and model;
(2) set up the reference signal that adds the fault size parameter, concrete grammar is as described below:
At first need calculate the linear velocity of rolling body motion and the pulse width that different faults causes in order to obtain this model
Bearing roller linear velocity s:
S=πdf
r
Pulse width p
x:
The pulse x (t) that can obtain thus the defective generation can be expressed as:
The reference signal that the impact that is produced by defective namely contains the fault degree parameter can be expressed as:
φ′
imp(p
2,u,f
2,dx,d,fr)=conv(x(t),φ
imp(p
2,u,f
2))
φ wherein
Imp(p
2, u, f
2) for its expression formula of decaying exponential function be:
Be φ '
Imp(p
2, u, f
2, dx, d, fr) and can be expressed as φ
Imp(p
2, u, f
2) with the convolution of x (t), (conv is the convolution algorithm symbol).
Wherein d is that the bearing path can be determined f according to bearing designation
rFrequently can record by the sensor of special measurement rotating speed for turning.d
x(unit: mm), u is the initial position of reference signal, p for fault diameter
2Can get the arbitrary value value for the damping coefficient of reference signal shock response in 1000-1200, the larger kurtosis of value is larger, f
2For the measured bearing damped natural frequency that dopes is that reference signal frequency needs according to f
1Value be defined as equaling f
1Value.
Parameter p in the model
2, u, f
2Value is identical with the value of using y (t) as with reference to signal the time.
(3) utilize phase-locked interpolator arithmetic that analyzed signal is processed, detailed process is:
Make d in the model
xSpan is [0.5, max (d
x)], get a value every 0.5, get altogether
Individual d
xValue, wherein max (d
x) be maximum fault diameter, this
Individual amount is correspondence 1,2 from small to large,,
This
Individual integer as the fault level label, is NO with letter representation.I.e. 0.5 corresponding numeral 1,1 corresponding numeral 2,, 5 corresponding numerals 10,, max (d
x) corresponding numeral
Be 0.5 each standards change d of 0.5 that increases according to initial value
xValue, d of every variation
xValue, utilize the phase-locked simplification expression formula of putting that signal is carried out single treatment namely:
The kurtosis index of signal V after every processing is once got and processed, the record kurtosis is with d
xChanging value.
Kurtosis ks can be expressed as:
Wherein V is the analyzed signal after phase-locked amplification is processed, and mean is mean value.
(4) calculate fault degree eigenwert A:
What try to achieve
Choose three maximum kurtosis values among the individual kurtosis
,
,
, record their corresponding integer fault level label NO (max
1(d
x)), NO (max
2(d
x)), NO (max
3(d
x)) (maximum kurtosis value
Corresponding fault level label is NO (max
1(d
x)), second largest kurtosis value
Corresponding fault level label is NO (max
2(d
x)), the third-largest kurtosis value
Corresponding fault level label is NO (max
3(d
x))).Sum obtains the fault degree eigenwert divided by three kurtosis sums after multiply by corresponding failure grade label with kurtosis value, can be expressed as:
Fault degree and fault degree eigenwert graph of a relation that Fig. 6 obtains for test, wherein horizontal ordinate is fault degree value d
x, ordinate is the fault degree eigenwert A that tries to achieve.Along with the increase of fault, the fault degree eigenwert A that tries to achieve increases gradually, has proved the validity of method as seen from the figure.
Claims (6)
1. bearing fault and fault degree diagnostic method based on a phase-locked interpolator arithmetic may further comprise the steps:
(1) gathers bearing vibration signal as measured signal;
(2) set up reference signal;
(3) utilizing phase-locked interpolator arithmetic to carry out fault signature extracts and denoising;
(4) set up the reference signal that contains the fault degree parameter;
(5) utilize phase-locked interpolator arithmetic that analyzed signal is processed;
(6) calculating fault features value.
2. bearing fault and fault degree diagnostic method based on phase-locked interpolator arithmetic according to claim 1, the reference signal method for building up y (t) in the described step (2) is as follows:
Wherein y (t) is reference signal, k
2Be amplitude, wherein p
2Be the damping vibration attenuation characteristic of shock response, f
2Damped natural frequency corresponding to system.
3. bearing fault and fault degree diagnostic method based on phase-locked interpolator arithmetic according to claim 1 utilize phase-locked interpolator arithmetic to carry out fault signature and extract and denoising in the described step (3), according to the principle device process prescription of phase-locked amplification be:
Measured signal and reference signal then through the output signal v of multiplier (t) are:
v(t)=s(t)y(t)
Wherein s (t) is the bearing fault based on phase-locked interpolator arithmetic according to claim 1 and fault degree diagnostic method, and the measured signal that described step (1) collects with equation expression is:
s(t)=v
s(t)+n(t)
Wherein n (t) is random noise, v
s(t) for the bearing fault signal indication be:
That is:
k
1For the bearing fault signal amplitude is directly measured, p wherein
1For the damping coefficient of bearing fault signal shock response according to bearing parameter measurements and calculations, f
1Be that the corresponding frequency of spectrum peak position directly obtains T in the bearing fault signal spectrum figure resonant belt by reading for the frequency of each impact of bearing fault signal corresponding to the damped natural frequency of measured bearing
0For the bearing fault characteristics cycle calculates by bearing parameter, all t directly measure the time in the literary composition, all δ are the impulse function identifier in the literary composition, to be the maximum occurrences of time t be that the length of measured signal s (t) is directly measured to all max (t) in the literary composition, in the literary composition all e be exponential function represent the symbol, all n are bearing rotary week number in the literary composition, represent with nonnegative integer;
Measured signal is expressed as:
s(t)=v
s(t)+n(t)
Wherein n (t) is random noise, v
s(t) be the bearing fault signal;
According to phase-locking and amplification principle, if select the signal close with tested composition and choose suitable frequency and namely measured signal is amplified and de-noising; Traditional lock-in amplifier is all selected sinusoidal signal and is adjusted cosine composition in the suitable frequency abstraction analyzed signal; For the bearing fault signal, its most basic effective constituent is the impact-attenuating signal, therefore selects the impact-attenuating signal as reference signal extraction and amplifies effective trouble unit;
Establishing reference signal at this is:
Wherein y (t) is reference signal, k
2For the reference signal amplitude is set as arbitrarily number greater than 1, the larger enlargement factor of value is larger, p
2Get the arbitrary value value for the damping coefficient of reference signal shock response in 1000-1200, the larger kurtosis of value is larger, f
2For the measured bearing damped natural frequency that dopes is that reference signal frequency is defined as equaling f
1Value;
Measured signal and contrast signal then through the output signal of multiplier are:
v(t)=s(t)y(t)
Therefore signal v (t) is expressed as through the output signal V of integrator again:
Following formula is phase-locked amplification expression formula, and wherein u is the initial position of reference signal, and max (t) is that the length of measured signal s (t) is directly measured for the maximum occurrences of time t;
Following formula launches to be expressed as:
Following formula launches to be expressed as:
Order
Then V is expressed as:
V=V
s+V
n
It is the addition of fault-signal item and noise signal item; V wherein
sBe integration fault-signal item later, V
nBe integration noise signal item later;
Wherein
Work as f
1=f
2The time, noise item is zero, i.e. V=V
s, so the signal of integrator output:
Compare with measured signal s (t), dispelled noise contribution, and amplified the amplification of trouble unit; Namely realized the function of the amplification of denoising and fault-signal, wherein K is enlargement factor K=k
1k
2
Above-mentioned multiply each other and the process of integration is that phase-locked amplification process is expressed as:
Be the sequence of impacts signal of a string not Noise, and fault signature has obtained K and has doubly amplified; Know thus, when obtaining reference signal according to y (t) expression formula, the frequency values f in the reference signal
2The centre frequency of analyzed signal s (t) resonant frequency is f
2=f
1The time, namely to utilize phase-locking and amplification principle be phase-locked amplification expression formula amplifies and carries out effective de-noising to each failure impact signal of bearing.
4. bearing fault and fault degree diagnostic method based on phase-locked interpolator arithmetic according to claim 1 are set up the reference signal that contains the fault degree parameter and are comprised the steps: in the described step (4)
In order to diagnose the bearing fault size, utilize bearing fault signal accurate model as the reference signal, introduced fault size parameter d in the model
xAt first need calculate the linear velocity of rolling body motion and the pulse width that different faults causes in order to obtain this model
Bearing roller linear velocity s:
S=πdf
r
Pulse width p
x:
The pulse x (t) that obtains thus the defective generation is expressed as:
The reference signal that the impact that is produced by defective namely contains the fault degree parameter is expressed as:
Φ′
imp(p
2,u,f
2,dx,d,fr)=conv(x(t),φ
imp(p
2,u,f
2))
φ wherein
Imp(p
2, u, f
2) for its expression formula of decaying exponential function be:
Be Φ '
Imp(p
2, u, f
2, dx, d, fr) and be expressed as φ
Imp(p
2, u, f
2) with the convolution of x (t), conv is the convolution algorithm symbol;
Wherein d is that the bearing path is determined f according to bearing designation
rFor the sensor that turns frequently by special measurement rotating speed records; d
xBe fault diameter, unit: mm, u are the initial position of reference signal, get the arbitrary value value for the damping coefficient of reference signal shock response in 1000-1200, and the larger kurtosis of value is larger, f
2For the measured bearing damped natural frequency that dopes is that reference signal frequency needs according to f
1Value be defined as equaling f
1Value.
5. bearing fault and fault degree diagnostic method based on phase-locked interpolator arithmetic according to claim 1, utilize phase-locked interpolator arithmetic that analyzed signal is processed in the described step (5) and be described as:
Make d in the model
xSpan is [0.5, max (d
x)], get a value every 0.5, get altogether
Individual d
xValue, wherein max (d
x) be maximum fault diameter, this
Individual amount is correspondence 1,2 from small to large,, 1-
This
Individual integer as the fault level label, is NO with letter representation; I.e. 0.5 corresponding numeral 1,1 corresponding numeral 2,, 5 corresponding numerals 10,, max (d
x) corresponding numeral
Be 0.5 each standards change d of 0.5 that increases according to initial value
xValue, d of every variation
xValue, utilize the phase-locked simplification expression formula of putting that signal is carried out single treatment namely:
The kurtosis index of signal V after every processing is once got and processed, the record kurtosis is with d
xChanging value;
Kurtosis ks is expressed as:
Wherein V is the analyzed signal after phase-locked amplification is processed, || be absolute value sign, mean is mean value.
6. the process prescription that calculates fault eigenvalue in bearing fault and the fault degree diagnostic method based on phase-locked interpolator arithmetic according to claim 1, described step (6) is:
What try to achieve
Choose three maximum kurtosis values among the individual kurtosis
,
,
, record their corresponding integer fault level label NO (max
1(d
x)), NO (max
2(d
x)), NO (max
3(d
x)), maximum kurtosis value
Corresponding fault level label is NO (max
1(d
x)), second largest kurtosis value
Corresponding fault level label is NO (max
2(d
x)), the third-largest kurtosis value
Corresponding fault level label is NO (max
3(d
x))); Sum obtains the fault degree eigenwert divided by three kurtosis sums after multiply by corresponding failure grade label with kurtosis value, is expressed as:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310000389.5A CN103076174B (en) | 2013-01-01 | 2013-01-01 | A kind of bearing fault based on lock-in amplify algorithm and fault degree diagnostic method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310000389.5A CN103076174B (en) | 2013-01-01 | 2013-01-01 | A kind of bearing fault based on lock-in amplify algorithm and fault degree diagnostic method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103076174A true CN103076174A (en) | 2013-05-01 |
CN103076174B CN103076174B (en) | 2015-10-28 |
Family
ID=48152779
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310000389.5A Expired - Fee Related CN103076174B (en) | 2013-01-01 | 2013-01-01 | A kind of bearing fault based on lock-in amplify algorithm and fault degree diagnostic method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103076174B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103712796A (en) * | 2013-12-13 | 2014-04-09 | 中国航空工业集团公司北京航空精密机械研究所 | Bearing fatigue signal detection system |
CN104156585A (en) * | 2014-08-04 | 2014-11-19 | 北京工业大学 | Double-dictionary matching pursuit bearing fault degree evaluation method based on multiple index properties |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03207284A (en) * | 1989-12-29 | 1991-09-10 | Ntn Corp | Motor load detector |
CN1321241A (en) * | 1999-08-25 | 2001-11-07 | 通用电气公司 | Vibration phasor monitoring system for rotating members |
JP2006153855A (en) * | 2004-10-26 | 2006-06-15 | Nsk Ltd | Anomaly diagnostic apparatus |
CN102455239A (en) * | 2010-10-29 | 2012-05-16 | 中国科学院沈阳自动化研究所 | Method utilizing phase-locking amplifier to detect weak signals of dynamic balancer and device for method |
-
2013
- 2013-01-01 CN CN201310000389.5A patent/CN103076174B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03207284A (en) * | 1989-12-29 | 1991-09-10 | Ntn Corp | Motor load detector |
CN1321241A (en) * | 1999-08-25 | 2001-11-07 | 通用电气公司 | Vibration phasor monitoring system for rotating members |
JP2006153855A (en) * | 2004-10-26 | 2006-06-15 | Nsk Ltd | Anomaly diagnostic apparatus |
CN102455239A (en) * | 2010-10-29 | 2012-05-16 | 中国科学院沈阳自动化研究所 | Method utilizing phase-locking amplifier to detect weak signals of dynamic balancer and device for method |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103712796A (en) * | 2013-12-13 | 2014-04-09 | 中国航空工业集团公司北京航空精密机械研究所 | Bearing fatigue signal detection system |
CN104156585A (en) * | 2014-08-04 | 2014-11-19 | 北京工业大学 | Double-dictionary matching pursuit bearing fault degree evaluation method based on multiple index properties |
CN104156585B (en) * | 2014-08-04 | 2018-06-01 | 北京工业大学 | A kind of doubledictionary match tracing bearing fault degree evaluation method based on many index feature |
Also Published As
Publication number | Publication date |
---|---|
CN103076174B (en) | 2015-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Detecting of transient vibration signatures using an improved fast spatial–spectral ensemble kurtosis kurtogram and its applications to mechanical signature analysis of short duration data from rotating machinery | |
CN103575523B (en) | The rotary machinery fault diagnosis method of kurtosis-envelope spectrum analysis is composed based on FastICA- | |
Zhang et al. | Adaptive maximum second-order cyclostationarity blind deconvolution and its application for locomotive bearing fault diagnosis | |
Singh et al. | Rolling element bearing fault diagnosis based on over-complete rational dilation wavelet transform and auto-correlation of analytic energy operator | |
Zhao et al. | Detection and recovery of fault impulses via improved harmonic product spectrum and its application in defect size estimation of train bearings | |
Patel et al. | Defect detection in deep groove ball bearing in presence of external vibration using envelope analysis and Duffing oscillator | |
Chen et al. | IGIgram: An improved Gini index-based envelope analysis for rolling bearing fault diagnosis | |
Antoni et al. | The spectral kurtosis: application to the vibratory surveillance and diagnostics of rotating machines | |
CN101221066B (en) | Engineering nonlinear vibration detecting method | |
Chen et al. | Optimal frequency band selection using blind and targeted features for spectral coherence-based bearing diagnostics: A comparative study | |
Chen et al. | Compound fault identification of rolling element bearing based on adaptive resonant frequency band extraction | |
CN101886977B (en) | Method for self-adaptively detecting periodic instant component in signal | |
CN104165925B (en) | The centrifugal compressor half-opened impeller crack fault detection method of accidental resonance | |
Schmidt et al. | A methodology for identifying information rich frequency bands for diagnostics of mechanical components-of-interest under time-varying operating conditions | |
CN103471848A (en) | Rolling bearing fault feature extraction method based on independent component analysis and cepstrum theory | |
CN102759448B (en) | Gearbox fault detection method based on flexible time-domain averaging | |
Xu et al. | The Enfigram: A robust method for extracting repetitive transients in rolling bearing fault diagnosis | |
CN104215456B (en) | Plane clustering and frequency-domain compressed sensing reconstruction based mechanical fault diagnosis method | |
Dong et al. | A frequency-shifted bispectrum for rolling element bearing diagnosis | |
CN101403684A (en) | Diagnosis method for combined failure of rotating machinery based on blind source analysis technology | |
CN103792000A (en) | Method and device for detecting transient components in signal based on sparse representation | |
CN105699082B (en) | A kind of maximum humorous make an uproar of rarefaction compares deconvolution method | |
CN103018044A (en) | Bearing combined failure diagnosis method based on improved impact dictionary matching pursuit | |
Luo et al. | Multiple discolored cyclic harmonic ratio diagram based on meyer wavelet filters for rotating machine fault diagnosis | |
US11422024B2 (en) | Method for detecting a defect in a vibration sensor, associated device and computer program |
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 | ||
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: 20151028 |