CN101354312B - Bearing failure diagnosis system - Google Patents

Bearing failure diagnosis system Download PDF

Info

Publication number
CN101354312B
CN101354312B CN2008100702352A CN200810070235A CN101354312B CN 101354312 B CN101354312 B CN 101354312B CN 2008100702352 A CN2008100702352 A CN 2008100702352A CN 200810070235 A CN200810070235 A CN 200810070235A CN 101354312 B CN101354312 B CN 101354312B
Authority
CN
China
Prior art keywords
signal
value
bearing
vibration acceleration
enter
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
CN2008100702352A
Other languages
Chinese (zh)
Other versions
CN101354312A (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN2008100702352A priority Critical patent/CN101354312B/en
Publication of CN101354312A publication Critical patent/CN101354312A/en
Application granted granted Critical
Publication of CN101354312B publication Critical patent/CN101354312B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention provides a bearing fault diagnosis system which comprises an intelligent bearing provided with six output ports with a first output port and a second output port outputting vibration acceleration signals, a third output port and a fourth output port speed signals and a fifth output port and a sixth output temperature signals. The system is characterized in that: the six groups of signals output by the intelligent bearing are sent to the input port of an A/D conversion module, six groups of digital signals output by the A/D conversion module are sent to a processor internally provided with a status monitor and a fault recognizer, wherein, the status monitor monitors the status of the bearing according to the vibration acceleration digital signals, the speed digital signals and the temperature digital signals, and the fault recognizer utilizes a signal processing tool to get a fault judgment result. The system of the invention has the advantages of being capable of monitoring the status of the bearing in an on-line manner according to the three types of signals of the intelligent bearing on the basis of the current intelligent bearing provided with a compound sensor soas to judge which part a fault takes place in.

Description

The intelligent bearing fault diagnosis system
Technical field
The present invention relates to sensor, specifically, relate to the wheel shaft, wheel box of a kind of wheel shaft that is applied in high-speed railway vehicle (bullet train etc.) and sedan limousine, large truck and some other adopts the rotating machinery on-line monitoring and the fault diagnosis system of intelligent bearing supporting way.
Background technology
Bearing, as one of vital part in the plant equipment, the quality of its running status will directly influence equipment running status and usability.To bearing, the monitoring of the bearing running status on particularly large-scale dynamoelectric equipment axis system and the hot-short and fault detect have become the important research content and the application of fault diagnosis technology.Through facts have proved, the quality of bearing diagnosis system performance depends on its detectability to the early stage fault in service of bearing, promptly detects abnormal conditions before catastrophic failure takes place.A host of facts prove that many major accidents are all owing to mechanical disorder caused with regard to appearance before expected life.These accidents have illustrated the potential danger that bearing fault brought that is difficult to predict, development intelligent bearing technology, and creating online bearing diagnosis and early warning system is current bearing failure diagnosis Developing Trend in Technology.
The intelligent bearing technology also can be described as under the background of technology develop rapidlys such as machinery, electronics, communication and produces.Because the continuous development of science and technology, improving constantly of the measure of precision of machinery and automaticity, early stage bearing fault detection method has fallen behind, and melts Modern Transducer Technology, signal transmission and treatment technology and computer technology is arisen at the historic moment in " intellectuality " diagnostic techniques of one.The performance monitoring and the method for diagnosing faults of therefore, the primary element that rotates freely as providing of widespread use in the mechanical system---bearing also develop towards the direction of " intellectuality ".Definition for " intelligent bearing ", current a kind of representative saying is, the sensing device of integrated different purposes on the basis of traditional bearing, it is combined into one and forms the particular structure unit, carry out information processing by microcomputer again, reach the purpose of real time on-line monitoring.
Compare with traditional bearing, a remarkable advantage that is integrated with the intelligent bearing of microsensor is exactly the serviceability that can monitor oneself in real time, and the signal to noise ratio (S/N ratio) height.Secondly because intelligent bearing centre bearer and sensor are combined into an one-piece construction, thereby the product that can be used as an independent completion produce and buy, make things convenient for user's practical application.
The combination of bearing and sensor at present comprises external hanging type and embedded two kinds of versions.External hanging type is meant that sensor is not in inner ring, outer ring or the rolling body that is embedded into bearing, but is attached on the bearing.This structure can not destroyed the integrality (being that maximum stress and the distortion that bearing allows do not change) of bearing.Traditional external hanging type is attached to sensory package around the bearing, and the object of monitoring and scope are very limited.More difficult in practice widespread use.And do not change bearing contour dimension, and the embedded intelligence bearing that can carry out multiple parameter testing is a better developing direction of intelligent bearing.The embedded intelligence bearing unit then is that sensor is embedded in the bearing body, the advantage of this method be sensing device can be very near the generation source of measured signal.The intermediate interface of signal transmission reduces, and the signal of collection more can reflect the practical working situation of bearing, signal to noise ratio (S/N ratio) height truly.
But because where sensing element is embedded into bearing actually, and sensing device these key issues such as reasonable Arrangement on compound sensor, make Embedded application be restricted.Currently used damascene structures is that sensing device is concentrated on the bearing, and this can produce stress and concentrate, and the bad layout of lead, brings difficulty to assembling.In addition, the sensing device of existing intelligent bearing only adopts a kind of mostly, and modal is the intelligent bearing unit that has the intelligent bearing unit of acceleration (vibration) sensing device and have speed sensing device.At high speed vehicles such as bullet train, bullet trains, because the speed of a motor vehicle is too high, the bearing vibration aggravation must detect vibration signal, with the huge person and the economic loss of preventing to cause.For preventing the generation of hot-short locking phenomenon, need to detect in real time tach signal.In addition, in order to prevent since bearing fault cause cut an accident, and owing to the wearing and tearing that the bearing self-heating causes, burn and need detect the temperature signal of bearing inner race (axle), outer ring simultaneously.So just must design vibration acceleration sensing device, temperature sensing device, the integrated intelligent bearing unit of revolution speed sensing device.
The inventor has developed a kind of (application number: 2007100925746) of intelligent bearing with compound sensor, this intelligent bearing is provided with six road delivery outlets, wherein first and second road delivery outlet AI0, AI1 respectively export one group of vibration acceleration signal, third and fourth road delivery outlet AI2, AI3 respectively export one group of tach signal, and the 5th, six road delivery outlet AI4, AI5 respectively export one group of temperature signal.Vibration acceleration sensing device, temperature sensing device, revolution speed sensing device and bearing is integrated, can extract the fault-signal of bearing effectively.
But also lack how the state of bearing to be carried out on-line monitoring, the concrete technical method that fault is discerned according to vibration acceleration signal, temperature signal, tach signal.
Summary of the invention
The purpose of this invention is to provide a kind of intelligent bearing fault diagnosis system, can be on the basis of the intelligent bearing that has the band compound sensor now, vibration acceleration signal, temperature signal, tach signal according to intelligent bearing carry out on-line monitoring to the state of bearing, and fault is discerned.
For achieving the above object, a kind of intelligent bearing fault diagnosis system of the present invention, comprise intelligent bearing, this intelligent bearing is provided with six road delivery outlets, wherein first, two road delivery outlets are respectively exported one group of vibration acceleration signal, the 3rd, four road delivery outlets are respectively exported one group of tach signal, the 5th, six road delivery outlets are respectively exported one group of temperature signal, its key is: six groups of signals of described intelligent bearing output send to the input port of A/D modular converter, this A/D modular converter is a digital signal with six groups of conversion of signals, six groups of digital signals of A/D modular converter output send to processor, be provided with state monitor and Fault Identification device in this processor, wherein state monitor is according to described first, two groups of vibration acceleration digital signals, the 3rd, four groups of rotating speed digital signals and the 5th, six groups of temperature digital signals are monitored the state of bearing, and generation failure alarm signal, state monitor is exported this alerting signal and is driven described Fault Identification device, this Fault Identification device utilizes the signal Processing instrument, comprehensive described vibration acceleration digital signal and rotating speed digital signal draw the fault judgement result.
This diagnostic system utilizes the intelligent bearing technology, gather vibration acceleration, speed and the temperature signal of bullet train hub bearing, treatment circuit by intelligent bearing inside amplifies vibration acceleration signal, tach signal, and temperature signal amplifies and compensates.Conversion enters monitoring of software analysis to pretreated 6 tunnel simulating signals through A/D.Rely on temperature value, velocity amplitude and vibration acceleration statistic to carry out status monitoring, when monitoring variable surpasses preset threshold, the proof bearing breaks down, at this moment state monitor fault alarm, start the Fault Identification function of Fault Identification device simultaneously, failure judgement occurs on which parts of bearing.The Fault Identification function relies on method for processing signals that digital signal is analyzed and handled and realizes.
Need to prove: the collection of intelligent bearing data is a collection of finishing, and the number of the sampled point that each batch data held is relevant with the cache size of A/D modular converter.After a batch data reads in by the A/D modular converter, to carry out multiple operation,, carry out digital demonstration, graphic presentation, also will realize multiple functions such as status monitoring and Fault Identification by calculating such as to comprising a batch data of 8192 sampling point informations.Only finished all operations, the A/D modular converter just can read in the next group data does repetitive operation.
Be provided with in the described state monitor:
Be used for the device that commencing signal is gathered;
Be used to obtain the time domain statistical parameter threshold value S of vibration acceleration signal R, S KDevice;
Be used for device with channel number CH, the first variable N1, the second variable N2 zero setting;
Be used to calculate root-mean-square value RMS, the kurtosis COEFFICIENT K of AI0 passage vibration acceleration FDevice with temperature data value W;
Owing to can't directly obtain vibration performance, use vibration acceleration time domain statistic to reflect the size of vibrational energy.The RMS value has reflected the body vibration energy, and it can estimate the health status of bearing and more stable, but this parameter is responsive inadequately to the initial failure of bearing.The kurtosis COEFFICIENT K FBe dimensionless group, reflected the size of impact shock, can be used for the initial failure diagnosis, but it can reduce gradually along with the intensification of fault.This diagnostic system adopts the method for two time domain parameter comprehensive evaluations to carry out the monitoring of vibrational state.
Be used to judge that temperature data value W is whether greater than 120 device;
If W is greater than 120 for the temperature data value, then be introduced into the device that is used to generate heat alarm;
Enter again and be used to judge whether RMS surpasses its threshold value S RDevice;
If temperature data value W is not more than 120, then directly enters and describedly be used to judge whether RMS surpasses the device of its threshold value SR;
If RMS does not surpass its threshold value S R, then directly enter and be used to judge the kurtosis COEFFICIENT K FWhether surpass its threshold value S KDevice;
If RMS surpasses its threshold value S R, then enter and be used for the first variable N1 is added 1 device;
Be used to judge that the first variable N1 is whether greater than 5 device;
If greater than 5, then entering, the first variable N1 is used to generate first alerting signal, and to the device of the first variable N1 zero clearing;
Enter the described kurtosis COEFFICIENT K that is used to judge again FWhether surpass its threshold value S KDevice;
If the first variable N1 is not more than 5, then directly enter the described kurtosis COEFFICIENT K that is used to judge FWhether surpass its threshold value S KDevice;
If K FDo not surpass its threshold value S K, then directly enter and be used for channel number CH added 1 device;
If K FSurpass its threshold value S K, then enter and be used for the second variable N2 is added 1 device;
Be used to judge that the second variable N2 is whether greater than 5 device;
If greater than 5, then entering, the second variable N2 is used to generate second alerting signal, and to the device of the second variable N2 zero clearing;
Enter described 1 the device that is used for channel number CH added again;
If the second variable N2 is not more than 5, then directly enter described 1 the device that is used for channel number CH added;
Enter again and be used to judge that channel number CH is whether greater than 1 device;
If CH is not more than 1, then directly enter the root-mean-square value RMS, the kurtosis COEFFICIENT K that are used to calculate AI1 passage vibration acceleration FAnd the device of temperature data value W;
If CH greater than 1, then is introduced into the device that is used for channel number CH zero setting;
Return described root-mean-square value RMS, the kurtosis COEFFICIENT K that is used to calculate AI0 passage vibration acceleration again FAnd the device of temperature data value W;
Wherein said
RMS = ( 1 N Σ i = 1 N x i 2 ) 1 / 2 ,
K F = 1 N Σ i = 1 N x i 4 / ( 1 N Σ i = 1 N x i 2 ) 2
The vibration acceleration signal of described intelligent bearing output is that the form with aanalogvoltage exists, and unit is mv; Through becoming digital voltage signal after the A/D conversion, unit is mv; In processor digital voltage signal is carried out conversion again, obtain the vibration acceleration value xi of sampled point, unit is g, and its conversion formula is:
X i = X i sm
x i': digital voltage signal, unit are mv;
S: the sensitivity coefficient of vibrating sensing device, unit are mv/g;
M: the enlargement factor of intelligent bearing vibration acceleration signal;
N: bearing inner race rotates 2 and encloses counting of being collected, and is to be that basic calculation draws with the tach signal:
N = 120 f nk
F: the sample frequency of every batch data, unit is Hz;
K: port number;
N: rotating speed, unit are r/min;
Described temperature data value W is: W = ( k N Z Σ i = 1 N z k y i t ) / S W
N Z: total sampling number of every batch data;
K: port number;
y i: the value of the temperature digital signal sampling point after amplification and compensation, unit is mv;
T: the enlargement factor of intelligent bearing temperature signal;
S W: the temperature-sensitivity coefficient of intelligent bearing, unit are mv/ ℃.
Be provided with in the described Fault Identification device:
Be used to begin the device of Fault Identification;
Be used to obtain rolling body and on outer ring raceway, pass through frequency f OpDevice;
Be used to obtain rolling body and on inner ring raceway, pass through frequency f IpDevice;
Be used to obtain the gyrofrequency f of the relative retainer of rolling body BcDevice;
Be used to read vibration acceleration value x iDevice;
Be used to utilize Fourier transformer, draw vibration acceleration signal value x iThe device of frequency spectrum;
The analytic function of data comprises multiple signal processing methods such as cepstrum, filtering, relevant, envelope, can choose suitable analytical approach according to the characteristic of signal, and wherein the frequency spectrum of vibration acceleration is that method by Fourier transform obtains.
Be used in the frequency range less than f/2k, choosing the device of the 1st resonance peak P=1 greater than 3000;
Be used for determining the device of filter bandwidht B=500;
Be used for bandpass filtering, obtain the device of filtering after vibration acceleration time domain signal;
Be used to obtain the device of signal envelope spectrum after the filtering;
Be used for finding out envelope spectrum maximum amplitude Y mDevice;
Be used to judge Y m-f OpThe value device between-5~5 whether;
If Y m-f OpValue between-5~5, then enter the device that is used to generate the outer ring fault-signal;
If Y m-f OpValue outside-5~5, then enter and be used to judge Y m-f IpThe value device between-5~5 whether;
If Y m-f IpValue between-5~5, then enter the device that is used to generate the inner ring fault-signal;
If Y m-f IpValue outside-5~5, then enter and be used to judge Y m-f BcThe value device between-5~5 whether;
If Y m-f BcValue between-5~5, then enter the device that is used to generate the rolling body fault-signal;
If Y m-f BcValue outside-5~5, then enter and be used for filter bandwidht B and add 100 device;
Vibration acceleration signal is carried out the arrowband bandpass filtering treatment, do the Hilbert conversion then, the analytic signal Q (t) that signal after the conversion and filtered vibration acceleration signal constitute is:
Q(t)=q(t)+jq 1(t)
Q (t): filtered signal;
q 1(t): the signal after the Hilbert conversion.
Right | Q (t) | remake Fourier transform and obtain envelope spectrum, find out the maximum amplitude Y of envelope spectrum m, at first the fault characteristic frequency with bearing outer ring contrasts, if satisfy-5<Y m-f Op<5 condition is diagnosed as bearing outer ring and breaks down, and does not carry out the fault judgement of inner ring, rolling body more successively if do not satisfy condition, and Rule of judgment is respectively :-5<Y m-f Ip<5 ,-5<Y m-f Bc<5.If three conditions do not satisfy, illustrate the selection of filtering parameter improperly, the program redirect is so that reselect relevant filtering parameter, and same resonance peak is being increased progressively with the filtering bandwidth of 100Hz under as the prerequisite of analytic target, and promptly B adds 100.
Enter again and be used to judge that B is whether greater than 1000 device;
If B is less than or equal to 1000, then return the described bandpass filtering that is used for, obtain the device of filtering after vibration acceleration time domain signal;
If B greater than 1000, then enters the device that is used for choosing at frequency spectrum following 1 resonance peak P=P+1;
If still can't discern the generation of fault when bandwidth surpasses 1000Hz, next so that frequency is higher resonance peak judges again as analytic target, i.e. P=P+1.
Be used to judge that P is whether greater than the frequency sequence number P of high resonance peak nDevice;
If P is less than or equal to P n, then return the described device that is used for determining filter bandwidht B=500;
If P is greater than P n, then enter the device that is used for the output safety signal;
If to all having carried out analyzing the generation that also can't discern fault, i.e. P less than all resonance peaks in the frequency range of f/2k greater than 3000〉P n, judge the bearing non-fault so.
Wherein:
f op = n 2 f r ( 1 - d D cos α ) ,
f ip = n 2 f r ( 1 + d D cos α ) ,
f bc = D 2 d f r ( 1 + d 2 D 2 cos α 2 )
In the formula, D is the pitch diameter of bearing, i.e. the circle at place, rolling body center, and d is the diameter of rolling body, α is a contact angle, f r=f i-f 0Be the frequency that relatively rotates of Internal and external cycle, when outer ring fixedly the time, f rBe the rotational frequency of axle.
Native system adopts the method for resonance and demodulation to carry out the Fault Identification of bearing, thereby which parts is failure judgement occur on.Resonance demodulation technique is low frequency to be impacted the high-frequency resonance ripple evoked carry out envelope detection, obtains one that impact corresponding to low frequency and amplify broadening also the resonance and demodulation ripple.By this resonance and demodulation ripple is carried out spectrum analysis, judge the value and the type of fault.
Remarkable result of the present invention is: utilize the intelligent bearing technology, gather bearing vibration acceleration, rotating speed and temperature signal, the treatment circuit by intelligent bearing inside amplifies vibration acceleration signal, tach signal, and temperature signal amplifies and compensates.Conversion enters monitoring of software analysis to pretreated 6 tunnel simulating signals through A/D.Rely on temperature value, velocity amplitude and vibration acceleration statistic to carry out status monitoring, when monitoring variable surpasses preset threshold, the proof bearing breaks down, at this moment state monitor fault alarm, start the Fault Identification function of Fault Identification device simultaneously, the Fault Identification device adopts the method for resonance and demodulation to carry out the Fault Identification of bearing, thereby which parts is failure judgement occur on.
Description of drawings
Accompanying drawing 1 is a connection block diagram of the present invention;
Accompanying drawing 2 is the workflow diagram of state monitor;
Accompanying drawing 3 is the workflow diagram of Fault Identification device.
Embodiment
Below in conjunction with the drawings and specific embodiments the present invention is described in further detail.
As shown in Figure 1, a kind of intelligent bearing fault diagnosis system, by intelligent bearing 1, A/D modular converter 2, processor, state monitor and Fault Identification device are formed, wherein intelligent bearing 1 is provided with six road delivery outlets, wherein first, two road delivery outlet AI0, AI1 respectively exports one group of vibration acceleration signal, the 3rd, four road delivery outlet AI2, AI3 respectively exports one group of tach signal, the 5th, six road delivery outlet AI4, AI5 respectively exports one group of temperature signal, six groups of signals of described intelligent bearing 1 output send to the input port of A/D modular converter 2, this A/D modular converter 2 is a digital signal with six groups of conversion of signals, six groups of digital signals of A/D modular converter 2 outputs send to processor, be provided with state monitor and Fault Identification device in this processor, wherein state monitor is according to described first, two groups of vibration acceleration digital signals, the 3rd, four groups of rotating speed digital signals and the 5th, six groups of temperature digital signals are monitored the state of bearing, and generation failure alarm signal, state monitor is exported this alerting signal and is driven described Fault Identification device, this Fault Identification device utilizes the signal Processing instrument, comprehensive described vibration acceleration digital signal and rotating speed digital signal draw the fault judgement result.
This diagnostic system utilizes the intelligent bearing technology, gather vibration acceleration, speed and the temperature signal of bullet train hub bearing, treatment circuit by intelligent bearing inside amplifies vibration acceleration signal, tach signal, and temperature signal amplifies and compensates.Conversion enters monitoring of software analysis to pretreated 6 tunnel simulating signals through A/D.Rely on temperature value, velocity amplitude and vibration acceleration statistic to carry out status monitoring, when monitoring variable surpasses preset threshold, the proof bearing breaks down, at this moment state monitor fault alarm, start the Fault Identification function of Fault Identification device simultaneously, failure judgement occurs on which parts of bearing.The Fault Identification function relies on method for processing signals that digital signal is analyzed and handled and realizes.
Need to prove: the collection of intelligent bearing data is a collection of finishing, and the number of the sampled point that each batch data held is relevant with the cache size of A/D modular converter 2.After a batch data reads in by A/D modular converter 2, to carry out multiple operation,, carry out digital demonstration, graphic presentation, also will realize multiple functions such as status monitoring and Fault Identification by calculating such as to comprising a batch data of 8192 sampling point informations.Only finished all operations, A/D modular converter 2 just can read in the next group data does repetitive operation.
As shown in Figure 2: be provided with in the described state monitor:
Be used for the device that commencing signal is gathered;
Be used to obtain the time domain statistical parameter threshold value S of vibration acceleration signal R, S KDevice; S R, S KAs predetermined value, be stored in the storer in the processor.
Be used for device with channel number CH, the first variable N1, the second variable N2 zero setting;
Be used to calculate root-mean-square value RMS, the kurtosis COEFFICIENT K of AI0 passage vibration acceleration FDevice with temperature data value W;
Owing to can't directly obtain vibration performance, use vibration acceleration time domain statistic to reflect the size of vibrational energy.The RMS value has reflected the body vibration energy, and it can estimate the health status of bearing and more stable, but this parameter is responsive inadequately to the initial failure of bearing.The kurtosis COEFFICIENT K FBe dimensionless group, reflected the size of impact shock, can be used for the initial failure diagnosis, but it can reduce gradually along with the intensification of fault.This diagnostic system adopts the method for two time domain parameter comprehensive evaluations to carry out the monitoring of vibrational state.
Wherein said RMS = ( 1 N Σ i = 1 N x i 2 ) 1 / 2 , K F = 1 N Σ i = 1 N x i 4 / ( 1 N Σ i = 1 N x i 2 ) 2
The vibration acceleration signal of described intelligent bearing output is that the form with aanalogvoltage exists, and unit is mv; Through becoming digital voltage signal after the A/D conversion, unit is mv; In processor, digital voltage signal is carried out conversion again, obtain the vibration acceleration value x of sampled point i, unit is g, its conversion formula is:
X i = X i ′ sm
x i': digital voltage signal, unit are mv;
S: the sensitivity coefficient of vibrating sensing device, unit are mv/g;
M: the enlargement factor of intelligent bearing 1 vibration acceleration signal;
N: bearing inner race rotates 2 and encloses counting of being collected, and is to be that basic calculation draws with the tach signal:
N = 120 f nk
F: the sample frequency of every batch data, unit is Hz;
K: port number;
N: rotating speed, unit are r/min;
Be used to judge that temperature data value W is whether greater than 120 device;
If W is greater than 120 for the temperature data value, then be introduced into the device that is used to generate heat alarm;
Whether 120 ℃ be in an important indicator of safe work state as bearing, becomes the threshold value of temperature alarming.
Enter again and be used to judge whether RMS surpasses its threshold value S RDevice;
If temperature data value W is not more than 120, then directly enters and describedly be used to judge that whether RMS surpasses its threshold value S RDevice;
Described temperature data value W is: W = ( k N z Σ i = 1 N z k y i t ) / S W
N z: total sampling number of every batch data;
K: port number;
y i: the value of the temperature digital signal sampling point after amplification and compensation, unit is mv;
T: the enlargement factor of intelligent bearing 1 temperature signal;
s w: the temperature-sensitivity coefficient of intelligent bearing 1, unit are mv/ ℃.
If RMS does not surpass its threshold value S R, then directly enter and be used to judge the kurtosis COEFFICIENT K FWhether surpass its threshold value S KDevice;
If RMS surpasses its threshold value S R, then enter and be used for the first variable N1 is added 1 device;
Be used to judge that the first variable N1 is whether greater than 5 device;
If greater than 5, then entering, the first variable N1 is used to generate first alerting signal, and to the device of the first variable N1 zero clearing;
Enter the described kurtosis COEFFICIENT K that is used to judge again FWhether surpass its threshold value S KDevice;
If the first variable N1 is not more than 5, then directly enter the described kurtosis COEFFICIENT K that is used to judge FDo not surpass its threshold value S KDevice;
If K FDo not surpass its threshold value S K, then directly enter and be used for channel number CH added 1 device;
If K FSurpass its threshold value S K, then enter and be used for the second variable N2 is added 1 device;
Be used to judge that the second variable N2 is whether greater than 5 device;
If greater than 5, then entering, the second variable N2 is used to generate second alerting signal, and to the device of the second variable N2 zero clearing;
Enter described 1 the device that is used for channel number CH added again;
If the second variable N2 is not more than 5, then directly enter described 1 the device that is used for channel number CH added;
Enter again and be used to judge that channel number CH is whether greater than 1 device;
If CH is not more than 1, then directly enter the root-mean-square value RMS, the kurtosis COEFFICIENT K that are used to calculate AI1 passage vibration acceleration FAnd the device of temperature data value W;
If CH greater than 1, then is introduced into the device that is used for channel number CH zero setting;
Return described root-mean-square value RMS, the kurtosis COEFFICIENT K that is used to calculate AI0 passage vibration acceleration again FAnd the device of temperature data value W;
As shown in Figure 3: be provided with in the described Fault Identification device:
Be used to begin the device of Fault Identification;
Be used to obtain rolling body and on outer ring raceway, pass through frequency f OpDevice;
Be used to obtain rolling body and on inner ring raceway, pass through frequency f IpDevice;
Be used to obtain the gyrofrequency f of the relative retainer of rolling body BcDevice;
Rolling bearing has n rolling body, and then rolling body passes through frequency f on outer ring and inner ring raceway OpAnd f Ip, and the gyrofrequency f of the relative retainer of rolling body BcBe expressed as respectively:
f op = n 2 f r ( 1 - d D cos α ) , f i p = n 2 f r ( 1 + d D cos α ) , f bc = D 2 d f r ( 1 + d 2 d 2 cos α 2 )
D is the pitch diameter of bearing, i.e. the circle at place, rolling body center, and d is the diameter of rolling body, α is a contact angle, f r=f i-f 0Be the frequency that relatively rotates of Internal and external cycle, when outer ring fixedly the time, f rBe the rotational frequency of axle;
Be used to read vibration acceleration value x iDevice;
Be used to utilize Fourier transformer, draw vibration acceleration value x iThe device of frequency spectrum;
The analytic function of data comprises multiple signal processing methods such as cepstrum, filtering, relevant, envelope, can choose suitable analytical approach according to the characteristic of signal, and wherein the frequency spectrum of vibration acceleration is that method by Fourier transform obtains.
Be used in the frequency range less than f/2k, choosing the device of the 1st resonance peak P=1 greater than 3000;
Be used for determining the device of filter bandwidht B=500;
Be used for bandpass filtering, obtain the device of filtering after vibration acceleration time domain signal;
Be used to obtain the device of signal envelope spectrum after the filtering;
Be used for finding out envelope spectrum maximum amplitude Y mDevice;
Be used to judge Y m-f OpThe value device between-5~5 whether;
If Y m-f OpValue between-5~5, then enter the device that is used to generate the outer ring fault-signal;
If Y m-f OpValue outside-5~5, then enter and be used to judge Y m-f IpThe value device between-5~5 whether;
If Y m-f IpValue between-5~5, then enter the device that is used to generate the inner ring fault-signal;
If Y m-f IpValue outside-5~5, then enter and be used to judge Y m-f BcThe value device between-5~5 whether;
If Y m-f BcValue between-5~5, then enter the device that is used to generate the rolling body fault-signal;
If Y m-f BcValue outside-5~5, then enter and be used to filter bandwidht B to add 100 device;
Vibration acceleration signal is carried out the arrowband bandpass filtering treatment, do the Hilbert conversion then, the analytic signal Q (t) that signal after the conversion and filtered vibration acceleration signal constitute is:
Q(t)=q(t)+jq 1(t)
Q (t): filtered signal;
q 1(t): the signal after the Hilbert conversion.
Right | Q (t) | remake Fourier transform and obtain envelope spectrum, find out the maximum amplitude Y of envelope spectrum m, at first the fault characteristic frequency with bearing outer ring contrasts, if satisfy-5<Y m-f Op<5 condition is diagnosed as bearing outer ring and breaks down, and does not carry out the fault judgement of inner ring, rolling body more successively if do not satisfy condition, and Rule of judgment is respectively :-5<Y m-f Ip<5 ,-5<Y m-f Bc<5.If three conditions do not satisfy, illustrate the selection of filtering parameter improperly, the program redirect is so that reselect relevant filtering parameter, and same resonance peak is being increased progressively with the filtering bandwidth of 100Hz under as the prerequisite of analytic target, and promptly B adds 100.
Enter again and be used to judge that B is whether greater than 1000 device;
If B is less than or equal to 1000, then return the described bandpass filtering that is used for, obtain the device of filtering after vibration acceleration time domain signal;
If B greater than 1000, then enters the device that is used for choosing at frequency spectrum following 1 resonance peak P=P+1;
If still can't discern the generation of fault when bandwidth surpasses 1000Hz, next so that frequency is higher resonance peak judges again as analytic target, i.e. P=P+1.
Be used to judge that P is whether greater than the frequency sequence number P of high resonance peak nDevice;
If P is less than or equal to P n, then return the described device that is used for determining filter bandwidht B=500;
If P is greater than P n, then enter the device that is used for the output safety signal;
If to all having carried out analyzing the generation that also can't discern fault, i.e. P less than all resonance peaks in the frequency range of f/2k greater than 3000〉P n, judge bearing safety so, non-fault.

Claims (2)

1. intelligent bearing fault diagnosis system, comprise intelligent bearing (1), this intelligent bearing (1) is provided with six road delivery outlets, wherein first, two road delivery outlet (AI0, AI1) respectively export one group of vibration acceleration signal, the 3rd, four road delivery outlet (AI2, AI3) respectively export one group of tach signal, the 5th, six road delivery outlet (AI4, AI5) respectively export one group of temperature signal, it is characterized in that: six groups of signals of described intelligent bearing (1) output send to the input port of A/D modular converter (2), this A/D modular converter (2) is a digital signal with six groups of conversion of signals, six groups of digital signals of A/D modular converter (2) output send to processor (3), be provided with state monitor (4) and Fault Identification device (5) in this processor (3), wherein state monitor (4) is according to described first, two groups of vibration acceleration digital signals, the 3rd, four groups of rotating speed digital signals and the 5th, six groups of temperature digital signals are monitored the state of bearing, and generation failure alarm signal, state monitor (4) is exported this alerting signal and is driven described Fault Identification device (5), this Fault Identification device (5) utilizes the signal Processing instrument, comprehensive described vibration acceleration digital signal and rotating speed digital signal draw the fault judgement result;
Described state monitor is provided with in (4):
Be used for the device that commencing signal is gathered;
Be used to obtain the time domain statistical parameter threshold value S of vibration acceleration signal R, S KDevice;
Be used for device with channel number CH, the first variable N1, the second variable N2 zero setting;
Be used to calculate root-mean-square value RMS, the kurtosis COEFFICIENT K of AI0 passage vibration acceleration FDevice with temperature data value W;
Be used to judge that temperature data value W is whether greater than 120 device;
If W is greater than 120 for the temperature data value, then be introduced into the device that is used to generate heat alarm;
Enter again and be used to judge whether RMS surpasses its threshold value S RDevice;
If temperature data value W is not more than 120, then directly enters and describedly be used to judge that whether RMS surpasses its threshold value S RDevice;
If RMS does not surpass its threshold value S R, then directly enter and be used to judge the kurtosis COEFFICIENT K FWhether surpass its threshold value S KDevice;
If RMS surpasses its threshold value S R, then enter and be used for the first variable N1 is added 1 device;
Be used to judge that the first variable N1 is whether greater than 5 device;
If greater than 5, then entering, the first variable N1 is used to generate first alerting signal, and to the device of the first variable N1 zero clearing;
Enter the described kurtosis COEFFICIENT K that is used to judge again FThe device that whether surpasses its threshold value SK;
If the first variable N1 is not more than 5, then directly enter the described kurtosis COEFFICIENT K that is used to judge FWhether surpass its threshold value S KDevice;
If K FDo not surpass its threshold value S K, then directly enter and be used for channel number CH added 1 device;
If K FSurpass its threshold value S K, then enter and be used for the second variable N2 is added 1 device;
Be used to judge that the second variable N2 is whether greater than 5 device;
If greater than 5, then entering, the second variable N2 is used to generate second alerting signal, and to the device of the second variable N2 zero clearing;
Enter described 1 the device that is used for channel number CH added again;
If the second variable N2 is not more than 5, then directly enter described 1 the device that is used for channel number CH added;
Enter again and be used to judge that channel number CH is whether greater than 1 device;
If CH is not more than 1, then directly enter the root-mean-square value RMS, the kurtosis COEFFICIENT K that are used to calculate AI1 passage vibration acceleration FAnd the device of temperature data value W;
If CH greater than 1, then is introduced into the device that is used for channel number CH zero setting;
Return described root-mean-square value RMS, the kurtosis COEFFICIENT K that is used to calculate AI0 passage vibration acceleration again FAnd the device of temperature data value W;
Wherein said
Figure 610556DEST_PATH_RE-F2008100702352C00031
Figure 877589DEST_PATH_RE-F2008100702352C00032
The vibration acceleration signal of described intelligent bearing output is that the form with aanalogvoltage exists, and unit is mv; Through becoming digital voltage signal after the A/D conversion, unit is mv; In processor (3), digital voltage signal is carried out conversion again, obtain the vibration acceleration value x of sampled point i, unit is g, its conversion formula is:
Figure 76489DEST_PATH_RE-F2008100702352C00033
x i': digital voltage signal, unit are mv;
S: the sensitivity coefficient of vibrating sensing device, unit are mv/g;
M: the enlargement factor of intelligent bearing (1) vibration acceleration signal;
N: bearing inner race rotates 2 and encloses counting of being collected, and is to be that basic calculation draws with the tach signal:
Figure 745368DEST_PATH_RE-F2008100702352C00034
F: the sample frequency of every batch data, unit is Hz;
K: port number;
N: rotating speed, unit are r/min;
Described temperature data value W is:
Figure 738732DEST_PATH_RE-F2008100702352C00035
N z: total sampling number of every batch data;
K: port number;
y i: the value of the temperature digital signal sampling point after amplification and compensation, unit is mv;
T: the enlargement factor of intelligent bearing (1) temperature signal;
S W: the temperature-sensitivity coefficient of intelligent bearing (1), unit are mv/ ℃.
2. intelligent bearing fault diagnosis system according to claim 1 is characterized in that:
Be provided with in the described Fault Identification device (5):
Be used to begin the device of Fault Identification;
Be used to obtain rolling body and on outer ring raceway, pass through frequency f OpDevice;
Be used to obtain rolling body and on inner ring raceway, pass through frequency f IpDevice;
Be used to obtain the gyrofrequency f of the relative retainer of rolling body BcDevice;
Be used to read vibration acceleration value x iDevice;
Be used to utilize Fourier transformer, draw vibration acceleration signal value x iThe device of frequency spectrum;
Be used in the frequency range less than f/2k, choosing the device of the 1st resonance peak P=1 greater than 3000;
Be used for determining the device of filter bandwidht B=500;
Be used for bandpass filtering, obtain the device of filtering after vibration acceleration time domain signal;
Be used to obtain the device of signal envelope spectrum after the filtering;
Be used for finding out envelope spectrum maximum amplitude Y mDevice;
Be used to judge Y m-f OpThe value device between-5~5 whether;
If Y m-f OpValue between-5~5, then enter the device that is used to generate the outer ring fault-signal;
If Y m-f OpValue outside-5~5, then enter and be used to judge Y m-f IpThe value device between-5~5 whether;
If Y m-f IpValue between-5~5, then enter the device that is used to generate the inner ring fault-signal;
If Y m-f IpValue outside-5~5, then enter and be used to judge Y m-f BcThe value device between-5~5 whether;
If Y m-f BcValue between-5~5, then enter the device that is used to generate the rolling body fault-signal;
If Y m-f BcValue outside-5~5, then enter and be used for filter bandwidht B and add 100 device;
Enter again and be used to judge that B is whether greater than 1000 device;
If B is less than or equal to 1000, then return the described bandpass filtering that is used for, obtain the device of filtering after vibration acceleration time domain signal;
If B greater than 1000, then enters the device that is used for choosing at frequency spectrum following 1 resonance peak P=P+1;
Be used to judge that P is whether greater than the frequency sequence number P of high resonance peak nDevice;
If P is less than or equal to P n, then return the described device that is used for determining filter bandwidht B=500;
If P is greater than P n, then enter the device that is used for the output safety signal;
Wherein:
Figure DEST_PATH_FA20184056200810070235201C00031
Figure DEST_PATH_FA20184056200810070235201C00033
In the formula, D is the pitch diameter of bearing, i.e. the circle at place, rolling body center, and d is the diameter of rolling body, α is a contact angle, f r=f i-f 0Be the frequency that relatively rotates of Internal and external cycle, when outer ring fixedly the time, f rBe the rotational frequency of axle;
Vibration acceleration signal is carried out the arrowband bandpass filtering treatment, do the Hilbert conversion then, the analytic signal Q (t) that signal after the conversion and filtered vibration acceleration signal constitute is:
Q(t)=q(t)+jq 1(t)
Q (t): filtered signal;
q 1(t): the signal after the Hilbert conversion;
Wherein envelope spectrum is that to remake Fourier transform behind the mould resulting by analytic signal Q (t) each point is asked.
CN2008100702352A 2008-09-05 2008-09-05 Bearing failure diagnosis system Expired - Fee Related CN101354312B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008100702352A CN101354312B (en) 2008-09-05 2008-09-05 Bearing failure diagnosis system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008100702352A CN101354312B (en) 2008-09-05 2008-09-05 Bearing failure diagnosis system

Publications (2)

Publication Number Publication Date
CN101354312A CN101354312A (en) 2009-01-28
CN101354312B true CN101354312B (en) 2010-09-22

Family

ID=40307228

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008100702352A Expired - Fee Related CN101354312B (en) 2008-09-05 2008-09-05 Bearing failure diagnosis system

Country Status (1)

Country Link
CN (1) CN101354312B (en)

Families Citing this family (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9008997B2 (en) * 2009-10-26 2015-04-14 Fluke Corporation System and method for vibration analysis and phase analysis of vibration waveforms using dynamic statistical averaging of tachometer data to accurately calculate rotational speed
DE102010002294A1 (en) * 2010-02-24 2011-08-25 Siemens Aktiengesellschaft, 80333 System or method for determining a storage condition
CN101819092B (en) * 2010-03-25 2012-02-01 重庆大学 Coupling type intelligent bearing monitoring device arranged on bearing
CN102269644B (en) * 2010-06-07 2014-06-04 北京化工大学 Diagnosis method for impact type failure between rolling bearing and gear based on optimal self-adaptive wavelet filter
JP5442553B2 (en) * 2010-07-21 2014-03-12 日立建機株式会社 Bearing damage detection method
CN102175453A (en) * 2011-01-21 2011-09-07 太原重工股份有限公司 Bearing on-line monitoring method of casting crane reducer
CN102435436B (en) * 2011-11-24 2015-05-27 电子科技大学 Method for evaluating degradation of state of fan bearing
EP2795280A4 (en) * 2011-12-21 2015-12-16 Skf Ab Method of monitoring a health status of a bearing with a warning device in a percentage mode
CN102564759B (en) * 2011-12-26 2014-10-29 重庆大学 Service life predicting system for automobile rear axle based on intelligent bearing
CN102809423A (en) * 2012-08-10 2012-12-05 上海卫星工程研究所 On-orbit satellite micro-vibration measurement system
CN103076168B (en) * 2013-01-09 2015-01-21 西安交通大学 Diagnosis method for mechanical faults of circuit breaker
CN104061178A (en) * 2013-03-21 2014-09-24 北京奇宏科技研发中心有限公司 Monitor device for service life of fan bearing
CN103335842A (en) * 2013-06-14 2013-10-02 江苏永钢集团有限公司 Bearing fault recognition system based on probability density curve and method thereof
CN103674538A (en) * 2013-12-18 2014-03-26 北京航天测控技术有限公司 Multi-fault mode identification method and device of swashplate of helicopter
CN103630361B (en) * 2013-12-19 2016-06-08 重庆大学 Intelligent bearing monitoring system and method thereof based on transition speed
CN103884506B (en) * 2014-03-25 2019-07-23 福建省工业设备安装有限公司 Production equipment on-line monitoring and fault diagnosis system and method based on wireless network
CN103868694B (en) * 2014-03-26 2016-01-20 东南大学 A kind of embedded variable speed bearing fault diagnosing apparatus
CN105502115A (en) * 2014-09-26 2016-04-20 刘一 Maintenance quality online assessment method, device and system based on elevator
CN104535097A (en) * 2014-11-05 2015-04-22 中国人民解放军第二炮兵工程大学 Small state monitoring and recording, fault diagnosing and alarming device
CN105784338B (en) * 2014-12-19 2018-06-29 安徽容知日新科技股份有限公司 The higher hamonic wave localization method of slewing base frequency
CN104913913B (en) * 2015-05-25 2019-10-25 广西科技大学 A kind of steel ring failure self-identifying intelligence system
CN105065639B (en) * 2015-06-19 2019-01-25 北京唐智科技发展有限公司 Transmission system resonance rupture is automatic to scout identification and the precaution and controlling method
CN104990709B (en) * 2015-08-07 2017-10-13 杨玉娇 Method for detecting Locomotive Bearing failure
US9978248B2 (en) * 2015-10-02 2018-05-22 Skf Canada Limited Low cost online condition monitoring
CN107121268A (en) * 2016-02-24 2017-09-01 王智中 Smart machine detection method
CN106017921B (en) * 2016-07-13 2018-07-06 广东石油化工学院 A kind of rotary machinery fault diagnosis method based on peak index
CN108225769B (en) * 2016-12-15 2020-02-14 唐智科技湖南发展有限公司 Mechanism diagnosis method for preventing holder characteristic spectrum from being misdiagnosed as fault spectrum
CN106872168A (en) * 2017-03-09 2017-06-20 武汉科技大学 A kind of method for diagnosing faults for rolling bearing
CN108507788A (en) * 2018-01-22 2018-09-07 内蒙古久和能源装备有限公司 A kind of rolling bearing fault degree judgment method
CN109856538B (en) * 2019-03-04 2021-07-09 中国矿业大学 Induction motor broken bar fault detection method based on short-time correction FFT
CN110315139B (en) * 2019-07-25 2020-11-10 马鞍山钢铁股份有限公司 Monitoring system for bearing degradation state of sliding seat bearing wheel on section steel saw machine
CN111060317A (en) * 2020-01-03 2020-04-24 上海电器科学研究所(集团)有限公司 Method for judging fault signal of rolling bearing of mining fan motor
CN114646467B (en) * 2020-12-21 2023-08-18 宇通客车股份有限公司 Driving motor bearing detection method suitable for whole vehicle environment
CN112857800A (en) * 2021-03-03 2021-05-28 嘉兴市特种设备检验检测院 Detection device and method for detecting reliability of high-speed elevator main shaft assembly
CN113295415B (en) * 2021-05-07 2022-06-10 湖北精瑞通流体控制技术有限公司 Bearing fault detection method based on multi-scale spectrum sensing technology
CN114580460A (en) * 2022-01-17 2022-06-03 西南交通大学 Railway vehicle wheel rail fault diagnosis method based on morphological filtering and HHT conversion

Also Published As

Publication number Publication date
CN101354312A (en) 2009-01-28

Similar Documents

Publication Publication Date Title
CN101354312B (en) Bearing failure diagnosis system
CN104535323B (en) A kind of train wheel Method for Bearing Fault Diagnosis based on angular domain time-domain and frequency-domain
CN104990709B (en) Method for detecting Locomotive Bearing failure
US7640139B2 (en) Abnormality diagnosing system for mechanical equipment
EP2522977B1 (en) Abnormality diagnosis device for rolling bearing, wind power generator, and abnormality diagnosis system
CN102483368B (en) Method for detecting structural defect in mechanical assembly including rotary member
CN108515984B (en) Wheel damage detection method and device
Wang et al. A smart sensing unit for vibration measurement and monitoring
CN106769033A (en) Variable speed rolling bearing fault recognition methods based on order envelope time-frequency energy spectrum
CN103471848A (en) Rolling bearing fault feature extraction method based on independent component analysis and cepstrum theory
WO2006120550A2 (en) Vehicle vibration analyzer
CN110553844B (en) Method and system for detecting misalignment fault of rotary machine
CN104931262A (en) Bearing fault diagnosis device used in cooperation with wheel set bearing running-in machine and diagnosis method thereof
Barbieri et al. Analysis of automotive gearbox faults using vibration signal
CN109253244A (en) A kind of multiple tooth wheel system big machinery gearbox fault detection method
KR20110131411A (en) Hub bearing test method
CN111595515B (en) Indirect detection method for unbalance fault of universal shaft
CN107101827A (en) A kind of low-speed heavy-loaded gear crack fault online test method
CN111458016A (en) Method for rapidly detecting vibration of drive axle assembly
CN102798413B (en) A kind of railway dynamic detection system
CN112594356B (en) High-speed railway and subway gear box operation monitoring and fault diagnosis system
CN111060302A (en) Multi-parameter equipment state comprehensive monitoring and diagnosing device and method
CN112284720B (en) Acoustic test-based fault diagnosis method for central transmission bevel gear of aircraft engine
CN105651512B (en) A kind of automobile drive axle master subtracts gear enclosed mass detection method
CN202885861U (en) Dynamic detection system for railway

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

Granted publication date: 20100922