CN102759397B - Efficient extraction method of friction fault characteristics in vibration signal of rotating shaft - Google Patents
Efficient extraction method of friction fault characteristics in vibration signal of rotating shaft Download PDFInfo
- Publication number
- CN102759397B CN102759397B CN201210237365.7A CN201210237365A CN102759397B CN 102759397 B CN102759397 B CN 102759397B CN 201210237365 A CN201210237365 A CN 201210237365A CN 102759397 B CN102759397 B CN 102759397B
- Authority
- CN
- China
- Prior art keywords
- signal
- centerdot
- friction
- rotating shaft
- degree coefficient
- 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
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention discloses an efficient extraction method of friction fault characteristics in a vibration signal of a rotating shaft, comprising the following steps of: selecting a plurality of monitoring sections along the direction of the rotating shaft; mounting an eddy-current transducer to monitor the vibration signal (t) of the rotating shaft; carrying out two times of differential continuously on time by virtue of the actually-monitored vibration signal and carrying out burr point enhancement operation on the signal processed by the two times of differential to obtain Z(t); calculating statistic characteristics of the Z(t) signal to obtain a sample average, a sample variance and a signal Z'(t) without an average value; calculating a signal distortion amplitude fault degree coefficient and a distortion frequency rate fault degree coefficient; calculating average values of the signal distortion amplitude fault degree coefficients and the distortion frequency rate fault degree coefficients in a plurality of periods, and taking the maximum value umax in the two average values as friction fault degree coefficients to judge the friction fault degree; and judging a friction position according to the friction fault degree coefficient of each measuring point of the rotating shaft. The method disclosed by the invention can effectively amplify the early-period friction fault characteristics and improve the diagnosis accuracy of early-period friction faults.
Description
Technical field
The present invention relates to tribological failure feature extracting method in a kind of shaft vibration signal.
Background technology
In order to prevent that steam, pressurized air, lubricated wet goods actuating medium from leaking outside, all kinds of rotating machineries such as steam turbine, generator, gas turbine, compressor, pump, blower fan have all been installed all kinds of sealings in rotating shaft, blade tip and bearing both sides, gap between sound parts is less, is in servicely easy to occur the friction between sound parts.After tribological failure occurs, unit vibration can increase, and affect unit safety stable operation, will cause the generation of the serious accidents such as permanent bending of rotor, bearing shell be cracked in the time of seriously.According to incompletely statistics, the more than 70% cambered axle accident of domestic electric system Turbo-generator Set is because tribological failure is caused.Efficient monitoring, analysis and the diagnostic method of tribological failure are most important for the safe operation of rotating machinery.
Monitoring and the analytical approach of tribological failure mainly contain at present: (1) thermal effect method.After friction occurs, the one week friction light and heavy degree difference in rotating shaft surface, can form thus the temperature difference and cause rotor to occur thermal deformation.The vibration variable quantity that monitoring thermal deformation causes can help to analyze and judge tribological failure.This method is used widely in Turbo-generator Set.The degree of functioning of this method is very relevant with friction light and heavy degree, and fault is difficult to accurately judgement in early days.The rotating machineries such as Turbo-generator Set also exist due to the caused thermal deformation fault of other factors, and the feature between friction and these factors is very similar, are difficult to make a distinction, thereby have affected the accuracy of method.(2) time domain method.According to waveform and the judgement of orbit of shaft center feature.Under friction impact force action, on waveform, likely there will be the phenomenons such as burr, it is disorderly that orbit of shaft center likely can become.When friction force is very large, above-mentioned phenomenon likely occurs, but now friction has entered middle and advanced stage, has lost the best opportunity of tribological failure diagnosis.(3) Spectrum Method.Under large impact force action, system becomes non-linear, there will be the high fdrequency components such as two frequencys multiplication and frequency tripling in vibration signal, according to the variation of vibration signal frequency spectrum medium-high frequency component, can judge tribological failure.When friction force is larger, said frequencies feature likely occurs.Under early stage Frotteurism, frequecy characteristic is not clearly.(4) power sensor method.At the back of friction member installing force sensor, during friction, power sensor can be experienced impulsive force signal, judges thus tribological failure.This method is mainly used in laboratory, cannot be applied to engineering reality.
In tribological failure, late period, friction characteristic is obvious, the large multipotency of said method is judged fault.While rubbing in early days or slightly, friction force is smaller, and thermal deformation is less, and the impact of countershaft vibrational waveform is also smaller, and the features such as wave form distortion and frequency spectrum are also not obvious.
Summary of the invention
The present invention proposes tribological failure feature highly effective extraction method in a kind of shaft vibration signal, and the method can be amplified early stage tribological failure feature effectively, improves the accuracy rate of diagnosis of early stage tribological failure.
In order to realize above object, the present invention adopts following technical scheme: rotating machinery sound tribological failure feature highly effective extraction method, comprises the steps:
1) in rotating shaft, fix reflecting strips, aim at the fixing photoelectric sensor of reflecting strips and be used for measuring rotational speed pulse signal; In axial direction select multiple monitorings cross section, eddy current sensor is installed respectively, eddy current sensor detects shaft vibration signal, and shaft vibration signal is designated as: y (t); Rotational speed pulse signal is as the analysis benchmark of the shaft vibration signal of multiple measuring points;
Wherein, Δ t is sampling time interval;
3) signal after second differential is made to burr point enhance operation:
4) statistical nature of calculating Z (t) signal, remove mean value:
Sample average
Sample variance
Remove the signal after average
5) in definition sampling period Z ' (t) in signal the point of amplitude > 3 σ be distortional point, maximal value Z '
max=max[Z ' is (t)], calculate the distortion n that counts and account for the ratio lambda of sample number N:
6) calculate signal distortion amplitude fault degree coefficient u in arbitrary sampling period
1j:
This formula is known technology, is not described further;
7) calculate signal distortion frequency fault degree coefficient u in arbitrary sampling period
2j:
This formula is known technology, is not described further;
8) mean value of distortion amplitude fault degree coefficient and aberration frequency fault degree coefficient in L the sampling period of calculating:
10) according to measuring point tribological failure degree coefficient, judge tribological failure degree, tribological failure degree coefficient value is larger, represents that tribological failure degree is more serious; According to the tribological failure degree coefficient of each measuring point in rotating shaft, judge frictional position, this measuring point tribological failure degree coefficient value is larger, represents the closer to friction point place.
Tribological failure feature highly effective extraction method in shaft vibration signal, wherein signal analysis frequency f
a>=10KHz, signals collecting frequency f
s>=2.56f
a, the sampling period is counted L>=8.
Beneficial effect: compared with prior art, tool of the present invention has the following advantages:
(1) original signal is done after second differential and burr enhance operation, the small distortional point in signal will be exaggerated, and can effectively amplify early stage tribological failure feature, improves the accuracy rate of diagnosis of early stage tribological failure;
(2) the tribological failure degree coefficient value calculating according to axial multiple measuring points distributes can judge frictional position;
(3) when tribological failure judges, considered the impact of distortional point quantity and distortion amplitude simultaneously;
(4) provide tribological failure and diagnose quantitative determination methods, make friction diagnosis quantitative from qualitative trend.
Sound friction detection method proposed by the invention can be monitored tribological failure and friction position effectively; the method can be embedded in all kinds of rotating machinery vibrating state protection TDM systems such as Turbo-generator Set, can effectively help technician to carry out the condition monitoring for rotating machinery such as Turbo-generator Set and fault diagnosis work.
Accompanying drawing explanation
Fig. 1 is tribological failure feature highly effective extraction method schematic flow sheet in shaft vibration signal of the present invention.
Fig. 2 is the rotor experiment table friction monitoring schematic diagram of the embodiment of the present invention.
Wherein, the 1st, motor; The 2nd, shaft coupling; The 3rd, photoelectric sensor; The 4th, bearing seat; The 5th, eddy current sensor; The 6th, friction screw and support; The 7th, rotor; The 8th, base.
Fig. 3 is without original signal waveform under Frotteurism.
Fig. 4 is original signal waveform under Frotteurism.
Fig. 5 is without under Frotteurism, original vibration signal being carried out to the waveform after second differential computing.
Fig. 6 is without under Frotteurism, the waveform after second differential being made to the waveform after burr enhance operation again.
Fig. 7 carries out the waveform after second differential computing to original vibration signal under Frotteurism.
Fig. 8 makes the waveform after burr enhance operation again to the waveform after second differential under Frotteurism.
Embodiment
The present embodiment is take the rotor experiment table shown in Fig. 1 as example, the early stage tribological failure highly effective extraction method analysis of carrying out.
The axis friction fault high efficiency extraction analytical approach that the present embodiment relates to is:
(1) on the shaft part exposing, paste reflecting strips, photoelectric sensor is aimed to reflecting strips, measure rotational speed pulse signal, the benchmark using this as the synchronous integer-period sampled analysis of multi-measuring point.
(2) at rotating shaft different cross section place, arrange that sensitivity is the non-contact turbulent flow sensor of 7.87mV/ μ m, monitoring shaft vibration signal y (t).(3) signalization collection, analytical parameters:
Analysis frequency f
a=20KHz
Frequency acquisition f
s=5.12KHz
Sampling period is counted L=20
(4) original signal is once differentiated:
(5) signal after a subdifferential is remake to a subdifferential:
(6) signal after second differential is made to burr enhance operation:
(7) statistical nature of calculating signal Z (t), remove mean value:
(8) definition distortional point is the point of amplitude > 3 σ in one-period.The n that counts by distorting in one-period accounts for the ratio lambda of sample number N and calculates distortional point frequency fault degree coefficient:
(9) by maximum distortion amplitude Z ' in one-period
max=max[Z ' is (t)] calculating distortion amplitude fault degree coefficient:
(10) calculate u in L cycle
1jand u
2jmean value:
(12) tribological failure criterion is as follows:
1) tribological failure degree
The larger expression tribological failure of the tribological failure degree coefficient value degree of measuring point is more serious.
2) friction position
In axial each measuring point, the cross section, some place of tribological failure degree coefficient value maximum is the most close tribological failure cross section.
Claims (2)
1. rotating machinery sound tribological failure feature highly effective extraction method, is characterized in that, comprises the steps:
1) in rotating shaft, fix reflecting strips, on the shaft part exposing, paste reflecting strips, aim at the fixing photoelectric sensor of reflecting strips and be used for measuring rotational speed pulse signal, described reflecting strips is to paste in rotating shaft along the axial direction of rotating shaft; In axial direction select multiple monitorings cross section, eddy current sensor is installed respectively, eddy current sensor detects shaft vibration signal, and shaft vibration signal is designated as: y (t); Rotational speed pulse signal is as the analysis benchmark of the shaft vibration signal of multiple measuring points:
2) will survey shaft vibration signal y (t) time will be done to two subdifferentials, obtain
Wherein, Δ t is sampling time interval;
3) signal after second differential is made to burr point enhance operation:
4) calculate the statistical nature of Z (t) signal, and remove mean value:
Sample average
Sample variance
Wherein N is sample number;
5) in definition sampling period Z ' (t) in signal the point of amplitude > 3 σ be distortional point, maximal value
Z '
max=max[Z ' is (t)], calculate the distortion n that counts and account for the ratio lambda of sample number N:
6) calculate signal distortion amplitude fault degree coefficient u in arbitrary sampling period
1j:
7) calculate signal distortion frequency fault degree coefficient u in arbitrary sampling period
2j:
8) mean value of distortion amplitude fault degree coefficient and aberration frequency fault degree coefficient in L the sampling period of calculating:
10) first, according to measuring point tribological failure degree coefficient, judge tribological failure degree, tribological failure degree coefficient value is larger, represents that tribological failure degree is more serious; Then, according to the tribological failure degree coefficient of each measuring point in rotating shaft, judge frictional position, this measuring point tribological failure degree coefficient value is larger, represents the closer to friction point place.
2. rotating machinery sound tribological failure feature highly effective extraction method according to claim 1, is characterized in that signal analysis frequency f
a>=10KHz, signals collecting frequency f
s>=2.56.f
a, the sampling period is counted L>=8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210237365.7A CN102759397B (en) | 2012-07-10 | 2012-07-10 | Efficient extraction method of friction fault characteristics in vibration signal of rotating shaft |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210237365.7A CN102759397B (en) | 2012-07-10 | 2012-07-10 | Efficient extraction method of friction fault characteristics in vibration signal of rotating shaft |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102759397A CN102759397A (en) | 2012-10-31 |
CN102759397B true CN102759397B (en) | 2014-04-16 |
Family
ID=47053937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210237365.7A Expired - Fee Related CN102759397B (en) | 2012-07-10 | 2012-07-10 | Efficient extraction method of friction fault characteristics in vibration signal of rotating shaft |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102759397B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106441641B (en) * | 2016-10-14 | 2019-01-29 | 广西科技大学 | A kind of measurement method of flexible axle frictional force in the bent state |
CN109813544B (en) * | 2019-03-26 | 2019-11-05 | 武汉众犇慧通科技有限公司 | A kind of rotating machinery Incipient Fault Diagnosis method and system based on on-line monitoring |
CN114757226A (en) * | 2022-04-03 | 2022-07-15 | 昆明理工大学 | Bearing fault characteristic enhancement method of parameter self-adaptive decomposition structure |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101556210A (en) * | 2002-05-16 | 2009-10-14 | Lg电子株式会社 | Device and method for detecting wrong operation of motor in washing machine |
CN102063547A (en) * | 2011-01-06 | 2011-05-18 | 广东电网公司电力科学研究院 | Modeling process of friction energy model of dynamic and static crushing and grinding faults of generator as well as fault detection method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3609982B2 (en) * | 2000-04-20 | 2005-01-12 | リオン株式会社 | Fault diagnosis method and apparatus |
JP2008268187A (en) * | 2007-03-26 | 2008-11-06 | Nippon Steel Corp | Method and device for diagnosing abnormality of extremely low speed rotary machine |
-
2012
- 2012-07-10 CN CN201210237365.7A patent/CN102759397B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101556210A (en) * | 2002-05-16 | 2009-10-14 | Lg电子株式会社 | Device and method for detecting wrong operation of motor in washing machine |
CN102063547A (en) * | 2011-01-06 | 2011-05-18 | 广东电网公司电力科学研究院 | Modeling process of friction energy model of dynamic and static crushing and grinding faults of generator as well as fault detection method |
Non-Patent Citations (1)
Title |
---|
李洁等.发电机转子同相振动分析与诊断.《汽轮机技术》.2012,第54卷(第3期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN102759397A (en) | 2012-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gubran et al. | Shaft instantaneous angular speed for blade vibration in rotating machine | |
CN103645052A (en) | Wind turbine set gearbox remote online state monitoring and life assessment method | |
CN104089693B (en) | Rotation shaft torsional vibration modality online measurement method for turbo generator set | |
CN102519579A (en) | Detection method of turbine collision and abrasion faults | |
CN102759397B (en) | Efficient extraction method of friction fault characteristics in vibration signal of rotating shaft | |
CN102966526B (en) | Based on the compressor low frequency flow induced vibrations class method for diagnosing faults that axial vibration is analyzed | |
CN102072764B (en) | Low-frequency vibration single-peak real-time analysis method for steam turbine generator set | |
CN102095491B (en) | Method for analyzing low-frequency vibration mutability of steam turboset in real time | |
Ali et al. | Observations of changes in acoustic emission parameters for varying corrosion defect in reciprocating compressor valves | |
CN104111108A (en) | Torsional vibration impact signal characteristic extracting method for rotating mechanism | |
CN111458016A (en) | Method for rapidly detecting vibration of drive axle assembly | |
Leite et al. | Bearing fault detection in induction machine using squared envelope analysis of stator current | |
CN104121178A (en) | Fault diagnosis system and method for fuel delivery pump | |
CN105987809A (en) | Centrifugal-compressor semi-open-type impeller crack detection method based on random resonance | |
Vercoutter et al. | Tip timing spectral estimation method for aeroelastic vibrations of turbomachinery blades | |
CN102175307B (en) | Real-time quantitative analysis method for low-frequency vibration spectrum array of steam turbine generator set | |
Lim et al. | Improved blade fault diagnosis using discrete Blade Passing Energy Packet and rotor dynamics wavelet analysis | |
Ma et al. | Compound fault diagnosis of wind turbine bearing under ultra-low speed operations using generalized sparse spectral coherence | |
CN102692303B (en) | High-efficiency identification method of steam excited vibration fault for steam turbine generator unit | |
CN102879084B (en) | On-line low-frequency vibration unstable state early warning method for steam turbine generator unit | |
CN102095564B (en) | Method for identifying fluctuation rub-impact fault of turbo generator set in real time | |
Thanagasundram et al. | Autoregressive based diagnostics scheme for detection of bearing faults | |
Ompusunggu et al. | Low-cost vibration sensor with low frequency resonance for condition monitoring of low speed bearings: a feasibility study | |
Azeem et al. | An octave analysis approach based on vibrational data for early detection of multiple faults in rolling bearing | |
CN102680243B (en) | Online judgment method for steam flow shock excitation fault of steam turbine generator unit |
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: 20140416 Termination date: 20170710 |
|
CF01 | Termination of patent right due to non-payment of annual fee |