CN101929917A - Fault diagnosis method for rotary machine - Google Patents

Fault diagnosis method for rotary machine Download PDF

Info

Publication number
CN101929917A
CN101929917A CN 201010130450 CN201010130450A CN101929917A CN 101929917 A CN101929917 A CN 101929917A CN 201010130450 CN201010130450 CN 201010130450 CN 201010130450 A CN201010130450 A CN 201010130450A CN 101929917 A CN101929917 A CN 101929917A
Authority
CN
China
Prior art keywords
vibration
vector
rotor
vector area
frequency
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
Application number
CN 201010130450
Other languages
Chinese (zh)
Other versions
CN101929917B (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN2010101304504A priority Critical patent/CN101929917B/en
Publication of CN101929917A publication Critical patent/CN101929917A/en
Application granted granted Critical
Publication of CN101929917B publication Critical patent/CN101929917B/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 relates to a fault diagnosis method for a rotary machine, comprising the following steps: acquiring original vibration signals of two perpendicular channels on a vibration measurement section, and carrying out anti-aliasing filtering, low pass filtering and the like on collected signals; removing unnecessary interference noise signals, adopting complex FFT conversion to calculate a primary vibration vector a and a secondary vibration vector b; further calculating relevant numerical values of a vector area and a vector vibration ratio and obtaining graphs of the vector area and the vector vibration ratio relative to various variables; and utilizing the vector area to express the vibration strength and direction of a measuring point under single frequency, and utilizing the vector vibration ratio to express eccentricity size, flat degree and forward and reverse precessional motion of a vector ellipse. The method has the advantage of good effect in diagnosing the steady-state faults and instantaneous faults of the rotary machine; and the vector area graph and the vector vibration ratio are combined to diagnose the faults, thus the erroneous judgment generated by only utilizing the vector area is avoided, the diagnosis accuracy is enhanced, and the diagnosis efficiency and precision are improved.

Description

A kind of method for diagnosing faults of rotating machinery
Technical field
Patent of the present invention belongs to the mechanical fault diagnosis field, is specifically related to rotary machinery fault diagnosis and control field.
Background technology
General orthogonal two probes in same cross section that adopt of the detection of rotating machinery are finished, but because the whirling motion characteristic of rotor, the detected information of arbitrary probe all can not reflect the actual state of unit operation on spectrum value and structure, difference is very big under many circumstances, can directly influence diagnostic result, cause erroneous judgement.In order to reflect the vibrational state of rotor objectively, the erroneous judgement of avoiding so causing, many scholars and researchist have done a large amount of research and application work, comprising domestic holographic spectral analysis technology, vow spectral analysis technology and external full spectral analysis technology entirely, these analytical approachs and technology have obtained great success in practical engineering application, but also find the place of some shortcomings, Quan shown in Fig. 1-a~Fig. 1-c, be followed successively by holographic spectrum, vow spectrum and compose the graph of a relation that each component changes with octave entirely.Holographic spectral analysis technology be under each frequency multiplication with the oval state that vibrates represented, but, directly adopt ellipse to represent to be difficult to expressed intact, so its figure spectral resolution is low owing to differ more for a long time at each harmonic amplitude.The full spectral analysis technology of vowing is to represent the intensity vibrated with oval major semi-axis size under each frequency multiplication, figure spectral resolution height, but be difficult to distinguish the equal and eccentricity anisochrouous vibrational state of major semi-axis.Full spectral analysis technology is to show on positive and negative frequency axis in oval two the positive radius of circle sizes of decomposing of following orbit of shaft center of each frequency multiplication, judge the precession direction by the amplitude size on the positive and negative frequency axis under the two same harmonic wave of comparison, figure spectral resolution height, but also need determine with the amplitude phase Calais on the positive and negative frequency axis of same harmonic wave owing to judge each frequency lower rotor part oscillation intensity size, so the expression of intensity is carried out indirectly, is difficult to provide the presentation of bulk strength structure.
Summary of the invention
The object of the present invention is to provide a kind of method for diagnosing faults of rotating machinery, can remedy the incomplete and not high shortcoming of resolution of existing analytical approach information representation, improve the accuracy of judging.
For achieving the above object, the technical scheme of the method for diagnosing faults of rotating machinery provided by the invention is as follows: a kind of stable state method for diagnosing faults of rotating machinery, and the step of this method is as follows:
(1) arranges two mutually perpendicular displacement transducers in the rotor cross section, gather the rotor oscillation signal: in vertical plane, the X transducer arrangements in the horizontal direction, Y transducer arrangements in the vertical direction, consistent with the rotor speed direction to the direction of Y by X, employing waits the tach signal collection and carries out integer-period sampled synchronously with key signal;
(2) signal of gathering is carried out anti-aliasing filtering, low-pass filtering, remove unnecessary interfering noise signal;
(3) adopt multiple FFT conversion, lead the arrow that shakes
Figure 2010101304504100002DEST_PATH_IMAGE002
, the pair arrow that shakes
Figure 2010101304504100002DEST_PATH_IMAGE004
Calculate;
(4) calculate the vector area of vowing ellipse
Figure 2010101304504100002DEST_PATH_IMAGE006
And vibration sagittal ratio
Figure 2010101304504100002DEST_PATH_IMAGE008
: vector area is that the orbit of shaft center that contains the precession direction is vowed oval area, with the measuring point oscillation intensity size under the size expression single-frequency of vector area value, the positive and negative precession of the positive and negative expression rotor of vector area value: promptly
Figure 2010101304504100002DEST_PATH_IMAGE010
Size represent oscillation intensity state under this frequency;
Figure 2010101304504100002DEST_PATH_IMAGE012
The time, rotor is positive precession,
Figure 2010101304504100002DEST_PATH_IMAGE014
The time, rotor is anti-precession;
Figure 560284DEST_PATH_IMAGE008
Be used to explain X axis vibration absolute value and the preceding magnitude relationship of Y axis vibration absolute value fusion, the relatively suffered exciting force size of rotor both direction:
Figure 2010101304504100002DEST_PATH_IMAGE016
=1 expression X-axis and Y axis vibration absolute value are equal, =
Figure 2010101304504100002DEST_PATH_IMAGE018
It is infinitely great that expression X axis vibration absolute value and Y axis vibration absolute value differ;
(5) vector area and vibration sagittal ratio the graph of a relation that draws with respect to frequency multiplication or other variable;
(6) in the stable state fault diagnosis, to vector area
Figure 70211DEST_PATH_IMAGE006
, vibration sagittal ratio
Figure 403103DEST_PATH_IMAGE008
And graph of a relation is analyzed: the characteristics of variables under the first read fortune fault characteristic frequency together, analyze possible fault type, the size of vibration sagittal ratio and axial vibration and amplitude are differentiated under the corresponding fault characteristic frequency of basis again, and be common in conjunction with the particular type of determining fault; When handling transient process information, in conjunction with the method for expressing of Bode diagram, the instantaneous graph of a relation of the vector area of drawing, pairing frequency of vector area value maximum and tachometer value are the resonant frequency and the critical rotary speed value of rotor-support-foundation system in the analysis collection of illustrative plates.
The shake computing method of arrow, vibration sagittal ratio and vector area of the main arrow that shakes, pair are as follows:
Figure 2010101304504100002DEST_PATH_IMAGE020
(1)
Figure 418469DEST_PATH_IMAGE002
In the formula,
Figure 342109DEST_PATH_IMAGE002
,
Figure 164572DEST_PATH_IMAGE004
The main arrow that shakes of expression, the pair arrow that shakes respectively, ε represents vibration sagittal ratio,
Figure 108302DEST_PATH_IMAGE006
The expression vector area.
Patent of the present invention has following function and characteristics:
1. this method has been used for reference the advantage of holographic spectrum with the oval vibrational state of describing, and has been integrated into ellipse and has quantized; Absorb in the full arrow spectrum major semi-axis and represented the advantage that oscillation intensity quantizes, be integrated into the vibration performance of different eccentricities; Utilize full spectrum to quantize and described the advantage that two positive radius of circles are judged the precession direction, be integrated into the feature that to distinguish vibrational state.
2. utilize vector area spectrogram and vibration sagittal ratio to combine fault is judged, the misjudgment phenomenon that has produced when having avoided independent use vector area has strengthened the accuracy of diagnosis.
3. than additive method, because the characteristic that vector area itself is had makes this method have characteristics highly sensitive, good stability.
Description of drawings
Fig. 1-a is the synoptic diagram of holographic spectrum;
Fig. 1-b is full synoptic diagram of vowing spectrum;
Fig. 1-c is the synoptic diagram of full spectrum;
Fig. 2-a is that vector area is vowed the curve map that changes with main shaking;
To be vector area with pair shake Fig. 2-b vows the curve map that changes;
Fig. 2-c is the curve map that vector area changes with octave;
Full Fig. 3-a vows in the stable state fault diagnosis that the spectral method master shakes to vow the curve map that changes with octave;
Full Fig. 3-b vows in the stable state fault diagnosis that the spectral method pair is shaken to vow the curve map that changes with octave;
Full Fig. 3-c vows the curve map that the spectral method phasing degree changes with octave in the stable state fault diagnosis;
Fig. 4-a is the curve map that vector area method vector area changes with octave in the stable state fault diagnosis;
Fig. 4-b is the curve map that vector area method vibration sagittal ratio changes with octave in the stable state fault diagnosis;
Fig. 5-a be in the transient fault diagnosis X-axis with the Bode diagram of rotation speed change;
Fig. 5-b be in the transient fault diagnosis Y-axis with the Bode diagram of rotation speed change;
Fig. 5-c is the full Bode diagram of vowing in the transient fault diagnosis;
Fig. 5-d be in the transient fault diagnosis vector area with the Bode diagram of rotation speed change.
Embodiment
Rotating machinery is along with the response of the change procedure of some technological parameter and operational factor is one group of complete process signal, and it carries out fault diagnosis and parameter recognition etc. and have important value for the stability of analytical equipment.By the analysis of one group of process signal, can analyze the vibration characteristics of rotor, the reason of inquiry amplitude and phase change, obtain the key parameters such as critical rotary speed of unit actual motion, determine the unbalance response of rotor, the resonance of research structure and element etc. are one of important evidence of fault diagnosis.
Be that existing holography is being composed among the present invention, the full spectrum of vowing, on the basis of full spectral analysis technology, systematically analyze each parameter and rotor-support-foundation system parameter, the relation of exciting force parameter, and the advantage of the existing rotary machinery fault diagnosis analytical approach of absorption, proposed to be used for a kind of new method of rotary machinery fault diagnosis---vector area method, promptly under each frequency multiplication, represent vibrational state with the area of the rotor axis movement locus ellipse that contains the precession direction, successively combine with vibration sagittal ratio and to be applied in the actual diagnosis, the result shows that it has unique diagnosis effect, be a kind of effective Analysis on Fault Diagnosis new method, be to the innovation of full arrow spectral analysis method and replenish.
Based on the rotor dynamics theory, under each response frequency, the axes movement locus of rotor is an ellipse.The vector vibration signal is one group of Dynamic Signal describing the rotor space vibration, is called for short vector signal; Vow that ellipse is to contain the orbit of shaft center ellipse that the vector signal of precession direction is described; The orbit of shaft center that contains the precession direction is vowed that oval area is referred to as vector area, and algebraic expression itself contains the vector meaning, the positive and negative precession of positive and negative representative, and numerical value is oval area value.
When using the vector area method to carry out fault diagnosis, at first will arrange the sensor acquisition signal, concrete operations are as follows: arrange two mutually perpendicular displacement transducers in the rotor cross section, gather the rotor oscillation signal.In vertical plane, the X transducer arrangements in the horizontal direction, Y transducer arrangements in the vertical direction, consistent with the rotor speed direction to the direction of Y by X, employing waits the tach signal collection also to carry out integer-period sampled synchronously with key signal.
According to full arrow spectral technology, the signal that collects is carried out anti-aliasing filtering, low-pass filtering etc., remove unnecessary interfering noise signal after, adopt multiple FFT conversion, leading shakes vows that a, pair are shaken and vows the calculating of b:
(1)
On the resultant pair of shaking according to the full master who vows that spectral technology is obtained is shaken the basis of vowing, with vowing ellipse area statement vibrational state feature:
Promptly
Figure 578634DEST_PATH_IMAGE004
Size represent oscillation intensity state under this frequency. The time, rotor is positive precession, The time, rotor is anti-precession.
Vibration sagittal ratio
Figure DEST_PATH_IMAGE028
, can calculate in conjunction with formula (1), obtain following result:
Figure DEST_PATH_IMAGE030
Figure 474506DEST_PATH_IMAGE016
Size is represented the size and the oval flat degree of the oval excentricity of track, eccentricity,
Figure 149201DEST_PATH_IMAGE016
Big more, eccentricity big more (ellipse is flat more),
Figure 859537DEST_PATH_IMAGE008
The positive and negative precession of positive and negative representative:
Figure 289381DEST_PATH_IMAGE008
For just representing the positive precession of rotor,
Figure 468690DEST_PATH_IMAGE008
Be the anti-precession of negative indication rotor.
Figure 743814DEST_PATH_IMAGE016
Can be used for understanding XThe shaft vibration absolute value and YMagnitude relationship before the shaft vibration absolute value merges, the relatively suffered exciting force size of rotor both direction.Can prove, if XThe axle and YThe shaft vibration absolute value is equal, then vibration sagittal ratio
Figure 75700DEST_PATH_IMAGE016
Equal 1; If XThe axle and YThe shaft vibration absolute value does not wait, then vibration sagittal ratio
Figure 410866DEST_PATH_IMAGE016
Greater than 1; If XThe shaft vibration absolute value is far longer than YThe shaft vibration absolute value (or XThe shaft vibration absolute value is far smaller than YThe shaft vibration absolute value), vibration sagittal ratio then
Figure 77471DEST_PATH_IMAGE016
Be tending towards infinitely great.
Vibration sagittal ratio can remedy the defective that new method produces erroneous judgement.As when minor axis is far smaller than major axis, null value appears in the vector area value that new method is drawn in collection of illustrative plates, can produce misjudgment phenomenon, but in conjunction with vibration sagittal ratio, then can successfully differentiate: because this moment, the vibration sagittal ratio absolute value occurred infinitely great in collection of illustrative plates, this frequency lower rotor part axle center vibration shape approximation is a straight line, and promptly phase differential is
Figure DEST_PATH_IMAGE032
Integral multiple, corresponding in such cases a kind of specific fault mode misaligns as shaft coupling.
As Fig. 2 (a~c) is depicted as the curve map that the vector area of new method of the present invention is drawn with different parameters, comparison diagram 1, it is complete that we can clearly relatively draw vector area method and holography spectrum, vow spectrum and full spectrometry to compare the advantage that is had as follows:
(1) used for reference the advantage of holographic spectrum, be integrated into the oval advantage that quantizes with the oval vibrational state of describing;
(2) use for reference in the full arrow spectrum major semi-axis and represented the advantage that oscillation intensity quantizes, be integrated into the vibration performance of different eccentricities;
(3) use for reference full spectrum and quantized and describe the advantage that two positive radius of circles are judged the precession directions, be integrated into the feature that to distinguish vibrational state;
(4) this method is compared with diagnostic method in the past, have at rotary machinery fault diagnosis and controlling party mask highly sensitive, and stable outstanding feature.
New method of the present invention both can be applied in the stable state fault diagnosis of rotating machinery, simultaneously when handling the transient process information of rotating machinery, equally also show outstanding characteristic, the implementation procedure that is applied as routine process in detail during the imbalance fault in stable state fault diagnosis diagnosis and transient fault are diagnosed respectively with regard to this method below.
Certain company double-suction centrifugal pump formula induced draft fan, since 1993 put into operation, the fan platform vibration was strong, had a strong impact on the ordinary production of system.Though, produce little effect through repeatedly testing and keeping in repair.Analyze with full arrow spectral technology, (shown in a~c), find that 1 frequency multiplication place oscillation intensity is big as Fig. 3.Big according to oscillation intensity under fault characteristic frequency 1 frequency multiplication, judge that this fault may be due to basis resonance or the imbalance, which kind of fault can't specifically differentiate is.With this paper method field data is analyzed, shown in Fig. 4 (a), found that area value is big under 1 frequency multiplication, the vertical radial vibration of Y is greater than the vibration on X transverse axis and the Z axle both direction, and these features are all similar to basis resonance or the uneven phenomenon of the failure that causes.
Suppose the vibration that basis resonance causes, so Ji Chu natural frequency is very approaching a little less than the frequency of operation of blower fan or the two, will produce the excessive vibration of some directions.But learn under 1 frequency multiplication by this paper method diagnostic analysis
Figure 343236DEST_PATH_IMAGE016
Be approximately 1, shown in Fig. 4 (b), by analysis before as can be known,
Figure 778897DEST_PATH_IMAGE016
=1 expression XThe axle and YThe shaft vibration absolute value equates, does not have excessive vibration in a certain direction in this explanation physical fault, is caused by basis resonance so fix a breakdown.
Suppose the vibration that imbalance causes, though there is phase differential in the vibration of each point, vertical diameter is to all bigger (promptly with the horizontal radial vibration Less) and 1 frequency multiplication with this paper methods analyst under Be approximately 1 and be consistent, for this reason the fault occurrence cause.Carry out balance test through dismounting and find that fault is just by due to the rotor unbalance.
Vector area method collection of illustrative plates and additive method single channel analytic approach, dual-channel analysis method (as holography spectrum, full spectrum, the full spectrum etc. of vowing) are compared, and the vector area value under more obvious, the single or fault characteristic frequency of fault characteristic frequency is bigger and single.By the further investigation failure mechanism, and, further sum up vibration sagittal ratio features different under the power frequency fault signature with new method of the present invention through a large amount of tests.Under identical fault characteristic frequency, can distinguish different faults effectively in conjunction with vibration sagittal ratio, may be used in the actual diagnosis, improve the accuracy of diagnosis.
The power frequency failure symptom that table 1 is summed up by the vector area method
? Characteristic frequency | vibration sagittal ratio | Axial vibration Orbit of shaft center Amplitude is with rotation speed change
Uneven 1X Near 1 Less Near circle Change little
Base flexible 1X Bigger ? Ellipse is more flat May reduce
Rotor bow
1X Near 1 Bigger Near circle Amplitude can drop to minimum value in the raising speed
The sound friction 1X, companion's higher hamonic wave Bigger ? " wedge angle " arranged Change little
In last table, axial vibration be dewater in the three dimensions square to and vertical direction (this both direction information fusion) outside, another channel axis is taken out to value and is assisted the characteristic quantity analyzed separately.Orbit of shaft center is the trace graphics after binary channels merges, i.e. the calculating figure of vector area, and the size of vector area value is the size of axle center trace graphics area value under each frequency just; And amplitude refers to the size of vector area here with " amplitude " in the rotation speed change, refers to herein under the natural frequency that the vector area value is with change in rotational speed.
In the start-stop process, rotor has experienced various rotating speeds, and its vibration signal is the response of rotor-support-foundation system to rotation speed change, it is the external reflection of rotor dynamics and failure symptom, comprised and be difficult to the abundant information that obtains at ordinary times, therefore, the start-stop process analysis procedure analysis is an important process of rotor monitoring.The same with full arrow spectral technology, the vector area method is when handling rotating machinery transient process information, equally also show outstanding characteristic---can accurately find the resonant frequency and the critical rotary speed value of rotor-support-foundation system, it is the method for expressing by Bode diagram, thereby accurately judge the critical rotary speed of rotor-support-foundation system, guidance is provided with the working speed of machine and avoids critical rotary speed, in order to avoid machine resonates.
(be followed successively by Bode diagram that X-axis draws, Bode diagram that Y-axis is drawn shown in a~d), vow the Bode diagram that Bode diagram and vector area method are drawn entirely, the Bode diagram that the vector area method is drawn is that the rotor axis movement locus area value under each rotating speed is plotted in the curve map that rectangular coordinate is fastened as Fig. 5.
Shown in Fig. 5-a, the vibration under the 3150r/s rotating speed is 110 to the maximum
Figure DEST_PATH_IMAGE034
, and in Fig. 5-c, can obtain, the vibration under the 3250r/s rotating speed is 148 to the maximum
Figure 9786DEST_PATH_IMAGE034
This shows that this traditional single channel analytic approach changes with the selection of X, Y passage, can't specifically determine the critical rotary speed value size of this equipment.Can analyze vibration maximum under the 3200r/s by full arrow Bode diagram 5-c, and the Bode diagram 5-d that draws with the vector area method equally accurately intuitive analysis go out vibration maximum under the 3200r/s, and vow that more full Bode diagram is more clear and intuitive.This figure ordinate is represented is not rotor value of shaking of a certain passage under the different rotating speeds in rising the car process, but the area value of rotor axis movement locus, thereby make the transient signal feature performance that rises the car process obvious especially, and this remarkable performance can not change with the passage variation as traditional analysis.

Claims (2)

1. the stable state method for diagnosing faults of a rotating machinery, it is characterized in that: the step of this method is as follows:
(1) arranges two mutually perpendicular displacement transducers in the rotor cross section, gather the rotor oscillation signal: in vertical plane, the X transducer arrangements in the horizontal direction, Y transducer arrangements in the vertical direction, consistent with the rotor speed direction to the direction of Y by X, employing waits the tach signal collection and carries out integer-period sampled synchronously with key signal;
(2) signal of gathering is carried out anti-aliasing filtering, low-pass filtering, remove unnecessary interfering noise signal;
(3) adopt multiple FFT conversion, lead the arrow that shakes
Figure DEST_PATH_IMAGE002
, the pair arrow that shakes
Figure DEST_PATH_IMAGE004
Calculate;
(4) calculate the vector area of vowing ellipse
Figure DEST_PATH_IMAGE006
And vibration sagittal ratio
Figure DEST_PATH_IMAGE008
: vector area is that the orbit of shaft center that contains the precession direction is vowed oval area, with the measuring point oscillation intensity size under the size expression single-frequency of vector area value, the positive and negative precession of the positive and negative expression rotor of vector area value: promptly
Figure DEST_PATH_IMAGE010
Size represent oscillation intensity state under this frequency;
Figure DEST_PATH_IMAGE012
The time, rotor is positive precession,
Figure DEST_PATH_IMAGE014
The time, rotor is anti-precession;
Figure 512676DEST_PATH_IMAGE008
Be used to explain X axis vibration absolute value and the preceding magnitude relationship of Y axis vibration absolute value fusion, the relatively suffered exciting force size of rotor both direction:
Figure DEST_PATH_IMAGE016
=1 expression X-axis and Y axis vibration absolute value are equal, = It is infinitely great that expression X axis vibration absolute value and Y axis vibration absolute value differ;
(5) vector area and vibration sagittal ratio the graph of a relation that draws with respect to frequency multiplication or other variable;
(6) in the stable state fault diagnosis, to vector area , vibration sagittal ratio
Figure 376093DEST_PATH_IMAGE008
And graph of a relation is analyzed: the characteristics of variables under the first read fortune fault characteristic frequency together, analyze possible fault type, the size of vibration sagittal ratio and axial vibration and amplitude are differentiated under the corresponding fault characteristic frequency of basis again, and be common in conjunction with the particular type of determining fault; When handling transient process information, in conjunction with the method for expressing of Bode diagram, the instantaneous graph of a relation of the vector area of drawing, pairing frequency of vector area value maximum and tachometer value are the resonant frequency and the critical rotary speed value of rotor-support-foundation system in the analysis collection of illustrative plates.
2. the method for diagnosing faults of rotating machinery according to claim 1 is characterized in that: the shake computing method of arrow, vibration sagittal ratio and vector area of the main arrow that shakes, pair are as follows:
Figure DEST_PATH_IMAGE020
(1)
Figure 917883DEST_PATH_IMAGE002
In the formula,
Figure 474761DEST_PATH_IMAGE002
,
Figure 588210DEST_PATH_IMAGE004
The main arrow that shakes of expression, the pair arrow that shakes respectively, ε represents vibration sagittal ratio,
Figure 777690DEST_PATH_IMAGE002
The expression vector area.
CN2010101304504A 2010-03-24 2010-03-24 Fault diagnosis method for rotary machine Expired - Fee Related CN101929917B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101304504A CN101929917B (en) 2010-03-24 2010-03-24 Fault diagnosis method for rotary machine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101304504A CN101929917B (en) 2010-03-24 2010-03-24 Fault diagnosis method for rotary machine

Publications (2)

Publication Number Publication Date
CN101929917A true CN101929917A (en) 2010-12-29
CN101929917B CN101929917B (en) 2012-08-29

Family

ID=43369217

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101304504A Expired - Fee Related CN101929917B (en) 2010-03-24 2010-03-24 Fault diagnosis method for rotary machine

Country Status (1)

Country Link
CN (1) CN101929917B (en)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102322942A (en) * 2011-09-13 2012-01-18 大连圣力来监测技术有限公司 Vibration signal vector ellipse spectrum analyzing method of rotating machinery
CN102998110A (en) * 2012-11-29 2013-03-27 西安交通大学 Rotary machine fault characteristic extraction method based on order-holospectrum principle
CN103323274A (en) * 2013-05-24 2013-09-25 上海交通大学 Rotating machinery condition monitoring and fault diagnosing system and method
CN103398831A (en) * 2013-07-09 2013-11-20 西安瑞特快速制造工程研究有限公司 Method for identifying supporting difference fault of rotary machine by procession trajectories of cross sections
CN103471841A (en) * 2013-09-30 2013-12-25 国家电网公司 Method for diagnosing vibration faults of rotary machine
CN103674545A (en) * 2013-11-26 2014-03-26 成都阜特科技股份有限公司 Mechanical fault detecting method
CN104374557A (en) * 2014-11-10 2015-02-25 西安陕鼓动力股份有限公司 Rotary machine second harmonic generation fault detecting method and system based on shaft center orbit
CN104677619A (en) * 2015-02-11 2015-06-03 中煤科工集团重庆研究院有限公司 Rotary machine fault feature extracting method on basis of complex signal bilateral spectrum
CN105067105A (en) * 2015-05-04 2015-11-18 西安交通大学 Kinetic parameter identification method utilizing rotary machine start and stop transient signal feature
CN105806622A (en) * 2016-05-26 2016-07-27 沈阳工程学院 Experimental method for measuring critical speed influencing factors for steam turbine rotor
CN106017956A (en) * 2016-05-18 2016-10-12 重庆大学 New method for fault diagnosis of rotating machine based on precession energy difference density spectrum analysis
CN106441840A (en) * 2016-09-09 2017-02-22 西安交通大学 Rotary machine pseudo-subsynchronous fault holographic diagnosis method based on model
CN106596025A (en) * 2016-12-06 2017-04-26 招商局重庆交通科研设计院有限公司 Highway tunnel hanging fan base stability detection method and system based on impulse response
CN108051189A (en) * 2017-11-20 2018-05-18 郑州工程技术学院 A kind of rotary machine fault characteristic extraction method and device
CN109488630A (en) * 2018-11-13 2019-03-19 上海金艺检测技术有限公司 Centrifugal blower rotor misalignment method for diagnosing faults based on harmonic wave relative indicatrix
CN112345248A (en) * 2019-08-09 2021-02-09 郑州工程技术学院 Fault diagnosis method and device for rolling bearing
CN112816200A (en) * 2021-01-15 2021-05-18 东方电气集团科学技术研究院有限公司 Rotary machine axis track identification method based on deep learning
CN112987613A (en) * 2021-02-24 2021-06-18 深圳中广核工程设计有限公司 Distributed rotating machinery state online monitoring and fault analysis system
CN113358282A (en) * 2021-06-04 2021-09-07 重庆大学 Composite material tail shaft low-speed orthogonal fusion dynamic balancing method on dynamic balancing machine
CN113785180A (en) * 2019-03-18 2021-12-10 蒙斯大学 Resonance frequency vibration test
CN115376301A (en) * 2022-07-25 2022-11-22 广东粤电博贺能源有限公司 Blade fracture early warning method and system of steam turbine, electronic equipment and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0585623A2 (en) * 1992-09-04 1994-03-09 Carl Schenck Ag Method for early detection of cracks in a rotating shaft
DE19855745A1 (en) * 1998-12-03 2000-08-03 Medav Digitale Signalverarbeit Fault diagnosis and quality grading system for production test of rotating machines uses sensing of rotation noise and oscillation with training of processor.
CN101266197A (en) * 2007-03-12 2008-09-17 斯奈克玛 Method of detecting damage to an antifriction bearing of a motor
CN101430239A (en) * 2008-11-28 2009-05-13 华北电力大学 Real-time diagnosis method for oil film whirl fault of large steam turbine-generator

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0585623A2 (en) * 1992-09-04 1994-03-09 Carl Schenck Ag Method for early detection of cracks in a rotating shaft
DE19855745A1 (en) * 1998-12-03 2000-08-03 Medav Digitale Signalverarbeit Fault diagnosis and quality grading system for production test of rotating machines uses sensing of rotation noise and oscillation with training of processor.
CN101266197A (en) * 2007-03-12 2008-09-17 斯奈克玛 Method of detecting damage to an antifriction bearing of a motor
CN101430239A (en) * 2008-11-28 2009-05-13 华北电力大学 Real-time diagnosis method for oil film whirl fault of large steam turbine-generator

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《机械科学与技术》 20090831 李利娟等 基于矢椭圆面积的双截面信息融合及其应用 第1023-1026页 第28卷, 第8期 2 *
《现代制造技术与装备》 20070831 陈先利等 全矢谱及其新方法在旋转机械故障中的应用 第11-12,61页 , 第4期 2 *
《设备管理与维修》 20080831 陈先利等 基于全矢谱的矢振比在旋转机械故障诊断中的应用 第40-41页 , 第8期 2 *

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102322942A (en) * 2011-09-13 2012-01-18 大连圣力来监测技术有限公司 Vibration signal vector ellipse spectrum analyzing method of rotating machinery
CN102998110A (en) * 2012-11-29 2013-03-27 西安交通大学 Rotary machine fault characteristic extraction method based on order-holospectrum principle
CN102998110B (en) * 2012-11-29 2015-07-01 西安交通大学 Rotary machine fault characteristic extraction method based on order-holospectrum principle
CN103323274B (en) * 2013-05-24 2015-10-14 上海交通大学 Condition monitoring for rotating machinery and fault diagnosis system and method
CN103323274A (en) * 2013-05-24 2013-09-25 上海交通大学 Rotating machinery condition monitoring and fault diagnosing system and method
CN103398831A (en) * 2013-07-09 2013-11-20 西安瑞特快速制造工程研究有限公司 Method for identifying supporting difference fault of rotary machine by procession trajectories of cross sections
CN103398831B (en) * 2013-07-09 2015-08-05 西安瑞特快速制造工程研究有限公司 The method of cross section whirl orbit identification rotating machinery supporting difference fault
CN103471841A (en) * 2013-09-30 2013-12-25 国家电网公司 Method for diagnosing vibration faults of rotary machine
CN103471841B (en) * 2013-09-30 2015-08-12 国家电网公司 A kind of rotating machinery vibrating failure diagnosis method
CN103674545A (en) * 2013-11-26 2014-03-26 成都阜特科技股份有限公司 Mechanical fault detecting method
CN103674545B (en) * 2013-11-26 2016-01-13 成都阜特科技股份有限公司 A kind of mechanical fault method for detecting
CN104374557B (en) * 2014-11-10 2017-02-15 西安陕鼓动力股份有限公司 Rotary machine second harmonic generation fault detecting method and system based on shaft center orbit
CN104374557A (en) * 2014-11-10 2015-02-25 西安陕鼓动力股份有限公司 Rotary machine second harmonic generation fault detecting method and system based on shaft center orbit
CN104677619A (en) * 2015-02-11 2015-06-03 中煤科工集团重庆研究院有限公司 Rotary machine fault feature extracting method on basis of complex signal bilateral spectrum
CN105067105B (en) * 2015-05-04 2018-04-17 西安交通大学 Utilize the kinetic parameter identification method of rotating machinery start and stop car transient signal feature
CN105067105A (en) * 2015-05-04 2015-11-18 西安交通大学 Kinetic parameter identification method utilizing rotary machine start and stop transient signal feature
CN106017956A (en) * 2016-05-18 2016-10-12 重庆大学 New method for fault diagnosis of rotating machine based on precession energy difference density spectrum analysis
CN105806622A (en) * 2016-05-26 2016-07-27 沈阳工程学院 Experimental method for measuring critical speed influencing factors for steam turbine rotor
CN106441840A (en) * 2016-09-09 2017-02-22 西安交通大学 Rotary machine pseudo-subsynchronous fault holographic diagnosis method based on model
CN106441840B (en) * 2016-09-09 2018-07-20 西安交通大学 A kind of pseudo- subsynchronous failure holographic diagnostics method of the rotating machinery based on model
CN106596025B (en) * 2016-12-06 2019-07-02 招商局重庆交通科研设计院有限公司 Stability of foundation of fan suspended in road tunnel detection method and system based on impulse response
CN106596025A (en) * 2016-12-06 2017-04-26 招商局重庆交通科研设计院有限公司 Highway tunnel hanging fan base stability detection method and system based on impulse response
CN108051189A (en) * 2017-11-20 2018-05-18 郑州工程技术学院 A kind of rotary machine fault characteristic extraction method and device
CN109488630A (en) * 2018-11-13 2019-03-19 上海金艺检测技术有限公司 Centrifugal blower rotor misalignment method for diagnosing faults based on harmonic wave relative indicatrix
CN113785180A (en) * 2019-03-18 2021-12-10 蒙斯大学 Resonance frequency vibration test
CN112345248A (en) * 2019-08-09 2021-02-09 郑州工程技术学院 Fault diagnosis method and device for rolling bearing
CN112345248B (en) * 2019-08-09 2022-11-25 郑州工程技术学院 Fault diagnosis method and device for rolling bearing
CN112816200A (en) * 2021-01-15 2021-05-18 东方电气集团科学技术研究院有限公司 Rotary machine axis track identification method based on deep learning
CN112987613A (en) * 2021-02-24 2021-06-18 深圳中广核工程设计有限公司 Distributed rotating machinery state online monitoring and fault analysis system
CN113358282A (en) * 2021-06-04 2021-09-07 重庆大学 Composite material tail shaft low-speed orthogonal fusion dynamic balancing method on dynamic balancing machine
CN113358282B (en) * 2021-06-04 2023-04-21 重庆大学 Low-speed orthogonal fusion dynamic balancing method for composite material tail shaft on dynamic balancing machine
CN115376301A (en) * 2022-07-25 2022-11-22 广东粤电博贺能源有限公司 Blade fracture early warning method and system of steam turbine, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN101929917B (en) 2012-08-29

Similar Documents

Publication Publication Date Title
CN101929917B (en) Fault diagnosis method for rotary machine
CN101639395B (en) Improved holographic dynamic balancing method of high-speed main shaft
CN100561158C (en) Holographic spot dynamic balance method based on precession vector
CN101556200B (en) Vector spectrum based dynamic balance method for flexible rotor
CN103364071B (en) Thin-wall cylinder modal test system and method for single-point laser continuous scanning vibration test
CN101907089B (en) Fault diagnosis method of compressor shafting based on three-dimensional space axle center orbit
CN100434890C (en) Flexible rotor holographic dynamic balancing method based on empirical mode decomposition
CN105300513A (en) Method and device for identifying resonance frequency multiplication number of blade vibration of multiple tip timing sensors
CN111337250A (en) Machine tool state fault diagnosis system and method based on virtual instrument
CN102072797B (en) Method for measuring unbalance of spindle in measurement of dynamic balance of tire and tire balancing machine
CN102928222B (en) Method for testing and identifying dynamic characteristic coefficients of sliding bearing
CN102721462B (en) Method for quickly computing Bode plot and Nyquist plot of rotary mechanical vehicle starting and parking processes
CN111397877B (en) Rotary machine beat vibration fault detection and diagnosis method
CN102175412A (en) Rotor torsion oscillation testing device and testing method thereof
CN105067106B (en) A kind of intershaft bearing vibration signals collecting method
CN103884482A (en) Compressor-based vibration test method and system
CN108827536A (en) Online dynamic balance detection system and detection method are assembled by a kind of combination revolution system
CN105092255A (en) Turbofan engine fan complete machine balancing method and system
CN103292958B (en) A kind of rotor based on model is without test mass unbalance parameter discrimination method
CN104101464A (en) Multi-wheel-disc rotor dynamic balancing test method based on rotating coordinate system
CN205426330U (en) Device is discerned to many apex timing sensor blade vibration resonance frequencies
CN103698088A (en) Testing method for asymmetric rigidity of turbonator rotating shaft
CN211121975U (en) Measuring device for dynamic impact load of wheel edge motor bearing
CN109282991A (en) A kind of DCT vehicle power drive system torsional vibration measurement and system
CN101451882B (en) Short time amplitude frequency spectrum array for single section shaft vibration analysis for mechanical rotor

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120829

Termination date: 20130324