CN105510010A - Characteristic parameter model for rotating machinery in misaligned structural anomaly state - Google Patents

Characteristic parameter model for rotating machinery in misaligned structural anomaly state Download PDF

Info

Publication number
CN105510010A
CN105510010A CN201510835820.7A CN201510835820A CN105510010A CN 105510010 A CN105510010 A CN 105510010A CN 201510835820 A CN201510835820 A CN 201510835820A CN 105510010 A CN105510010 A CN 105510010A
Authority
CN
China
Prior art keywords
characteristic parameter
parameter model
ratio
rotating machinery
rotary machine
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.)
Pending
Application number
CN201510835820.7A
Other languages
Chinese (zh)
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.)
Jiangnan University
Original Assignee
Jiangnan 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 Jiangnan University filed Critical Jiangnan University
Priority to CN201510835820.7A priority Critical patent/CN105510010A/en
Publication of CN105510010A publication Critical patent/CN105510010A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a characteristic parameter model for the rotating machinery in the misaligned structural anomaly state and belongs to the diagnostic field of mechanical faults, which mainly solves the fault detecting and fault type distinguishing problem during the fault occurrence early stage of the rotating machinery. Through analyzing the low-frequency spectrum characteristics of the rotating machinery in the misaligned structural anomaly state, associated characteristic parameters are constructed. In this way, the vibration characteristics of the rotating machinery in the above state can be reflected. Based on a resolution index (DI), the suitability of the characteristic parameter model and the resolution sensitivity of the characteristic parameter model are evaluated. Meanwhile, the problem that the fault type cannot be determined only based on conventional characteristic parameters can be solved.

Description

A kind of for the abnormal characteristic parameter model for misaligning of rotary machine configuration
Technical field
The present invention relates to a kind of for the abnormal characteristic parameter model for misaligning of rotary machine configuration, belonging to mechanical fault diagnosis field.
Background technology
Structural anormaly is also referred to as structural failure, and namely rotary machine breaks down, and at the low frequency region place indicating characteristic frequency spectrum of excited frequency.Structural anormaly causes the axle portion of rotating machinery overworked, is also the main cause causing bearing, gear etc. to break down.That is, structural anormaly may cause mechanical system to break down, and may cause huge production loss.Therefore, detection and specification configuration fault are very important to guarantee production efficiency and factory safety.In some sense, feature extraction can be described as the bottleneck problem in the research of Current mechanical fault diagnosis, and it is directly connected to the accuracy of fault diagnosis and the reliability of fault early prediction.In actual applications, often there is not simple relation one to one between fault and sign, if the corresponding multiple sign of a kind of fault possibility, otherwise a kind of sign also may caused by various faults, and these extract to fault signature and bring difficulty.
For a certain concrete fault type, its feature is showed by some physical parameter usually, and there is certain contacting with the strong and weak situation of each physical parameter, as long as the state of mechanical system changes, just inherently have influence on each dynamic physical parameter be associated with it, concern is comparatively wide, and these sensitive to fault, reliable and stable physical parameters are then used to build corresponding characteristic parameter model.Although the generation of certain fault type can cause the change of multiple physical parameter, the parameter that can be used as fault signature is limited, and Failure Characteristic Parameter model must possess high susceptibility, reliability and practicality.
Misaligning is the structural failure that the axle portion of rotating machinery often occurs.What is called misaligns, and refers to that two shaft centre lines with shaft coupling links up exist deviation, and as produced axis being parallel skew, axis produces cheap or both combination in angle.Though dimensionless group can the stage carries out fault diagnosis in early days for traditional characteristic parameter model such as margin index, flexure index, kurtosis index etc., be difficult to determine fault type.
Summary of the invention
Based on the technical matters of characteristic parameter model existence in the past, the present invention proposes a kind of novel characteristic parameter model for rotary machine configuration extremely for misaligning.Being the spectrum signature at low frequency place under condition of misalignment by analyzing rotary machine configuration abnormal, building 2 higher harmonics rates, direction of principal axis vibration index, vertical and horizontal amplitude ratio three parametric models.
Wherein, 2 higher harmonics rates are the spectrum value of 2 times of gyro frequencys that horizontal or vertical direction detects and gyro frequency
Wherein, the ratio of the spectrum value sum on the odds ratio turning the spectrum value sum under the spectrum value and gyro frequency 2 ~ 20 times of order harmonic frequencies that frequently component ratio is the gyro frequency that horizontal or vertical direction detects under the spectrum value of normal condition gyro frequency and gyro frequency 2 ~ 20 times of order harmonic frequencies, namely
Higher harmonics ratio be the spectrum value of 3 times of gyro frequencys and the spectrum value of gyro frequency that horizontal or vertical direction detects odds ratio on the ratio of the spectrum value of normal condition 3 times of gyro frequencys and the spectrum value of gyro frequency, namely N 2 = T d ( 3 f r ) / T d ( f r ) T n ( 3 f r ) / T n ( f r ) .
The present invention has the following advantages:
This characteristic parameter model resolution is high, practical, reliability is high.
Accompanying drawing explanation
Fig. 1 is the spectrogram of horizontal direction under centrifugal fan normal condition.
Fig. 2 is the spectrogram of vertical direction under centrifugal fan normal condition.
Fig. 3 is the spectrogram in centrifugal fan normal condition lower axis direction.
Fig. 4 is the spectrogram of horizontal direction under centrifugal fan non-equilibrium state.
Fig. 5 is the spectrogram of vertical direction under centrifugal fan non-equilibrium state.
Fig. 6 is the spectrogram in centrifugal fan non-equilibrium state lower axis direction.
Specific embodiments
Below in conjunction with accompanying drawing and concrete process of the test, the invention will be further described.
The diagnostic test of centrifugal fan: three accelerometers (PCBMA352A60) export from 5HZ to 60KHZ with the bandwidth measurement level of 10mV/g, vertical and axial vibration signal respectively under normal and non-equilibrium state.Vibration signal measured by accelerometer is transformed into signal recorder (scope scrambler DL750) after being amplified by transducer signal conditioner (PCBICP model 480C02).Former vibration signal in a frequency domain respectively as shown in figs. 1 to 6.These signals are all measured with constant speed (600 turns).The sample frequency of signal measurement is 50KHZ, and the sampling time is 20s.
Each parameter value is calculated as following table according to measurement result:
Utilize resolving index (DI) to carry out the resolution sensitivity of evaluating characteristic parameter as quality index, concrete proof procedure is as follows:
Suppose x 1and x 2be the value of the characteristic parameter that in measuring state 1 and state 2, signal calculates respectively, and meet normal distribution N (μ respectively 1, σ 1) and N (μ 2, σ 2).Here, μ and σ is mean value and the standard deviation of characteristic parameter.| x 2-x 1| value larger, distinguish the sensitivity of two states by characteristic parameter higher.Due to z=x 2-x 1also normal distribution N (μ is met 21, σ 1+ σ 2), there is the following density function about Z:
f ( z ) = 1 2 π ( σ 1 2 + σ 2 2 ) exp { { z - ( μ 2 - μ 1 ) } 2 2 ( σ 1 2 + σ 2 2 ) } - - - ( 1 )
Here, μ 2>=μ 1(in like manner, work as μ 1>=μ 2), probability can be calculated as follows:
P 0 = ∫ - ∞ 0 f ( z ) d z - - - ( 2 )
Here, 1-P 0be called as " resolution (DR) ".By following formula:
μ = 2 - ( μ 2 - μ 1 ) σ 1 2 + σ 2 2 - - - ( 3 )
In substitution formula (1) and (2), P 0obtained by following formula:
P 0 = 1 2 π ∫ - ∞ - D I exp ( - μ 2 2 ) d μ - - - ( 4 )
Here, being calculated as follows of DI (resolving index):
D I = μ 2 - μ 1 σ 1 2 + σ 2 2 Or D I = x 2 ‾ - x 1 ‾ σ 1 2 + σ 2 2 - - - ( 5 )
Clearly, DI value is larger, " resolution (DR=1-P 0) " value also will be larger, therefore, characteristic parameter also will be better.
According to above-mentioned proof procedure, the DI value of the character pair parameter calculated is as following table:
M 1 M 2 M 3
DI N-M 2.17765 10.2071 7.48199
Show the DI value of three characteristic parameters is all greater than 2.17 from above, according to area and the ordinate table of standardized normal distribution curve, as DI=2.17, resolution reaches 98.499%, therefore, three characteristic parameter models all reach more than 98.499% to rotary machine configuration is abnormal for the diagnosis resolution misaligning fault.

Claims (4)

1. the characteristic parameter model for rotary machine configuration extremely for misaligning, it is characterized in that: by analyzing the abnormal spectrum signature for low frequency place under condition of misalignment of rotary machine configuration, specifically comprising the characteristic parameter model of 2 higher harmonics rates, direction of principal axis vibration index, vertical and horizontal amplitude ratio.
2. the characteristic parameter model for rotary machine configuration extremely for misaligning according to right 1, it is characterized in that 2 times described higher harmonics rates be the spectrum value of 2 times of gyro frequencys and the spectrum value of gyro frequency that horizontal or vertical direction detects odds ratio on the ratio of the lower spectrum value of 2 times of gyro frequencys of normal condition and the spectrum value of gyro frequency, namely
3. the characteristic parameter model for rotary machine configuration extremely for misaligning according to right 1, it is characterized in that described direction of principal axis vibration index be the root-mean-square value of vibration signal on the axis direction and horizontal direction that detect odds ratio on the ratio of root-mean-square value of vibration signal in normal condition lower axis direction and horizontal direction, namely
4. the characteristic parameter model for rotary machine configuration extremely for misaligning according to right 1, it is characterized in that described vertical and horizontal amplitude ratio be the root-mean-square value of vibration signal in the vertical direction and horizontal direction that detect odds ratio on the ratio of root-mean-square value of vibration signal under normal condition in vertical direction and horizontal direction, namely
CN201510835820.7A 2015-11-26 2015-11-26 Characteristic parameter model for rotating machinery in misaligned structural anomaly state Pending CN105510010A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510835820.7A CN105510010A (en) 2015-11-26 2015-11-26 Characteristic parameter model for rotating machinery in misaligned structural anomaly state

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510835820.7A CN105510010A (en) 2015-11-26 2015-11-26 Characteristic parameter model for rotating machinery in misaligned structural anomaly state

Publications (1)

Publication Number Publication Date
CN105510010A true CN105510010A (en) 2016-04-20

Family

ID=55718153

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510835820.7A Pending CN105510010A (en) 2015-11-26 2015-11-26 Characteristic parameter model for rotating machinery in misaligned structural anomaly state

Country Status (1)

Country Link
CN (1) CN105510010A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109488630A (en) * 2018-11-13 2019-03-19 上海金艺检测技术有限公司 Centrifugal blower rotor misalignment method for diagnosing faults based on harmonic wave relative indicatrix
CN110553844A (en) * 2019-07-24 2019-12-10 西安交通大学 Method and system for detecting misalignment fault of rotary machine

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102721545A (en) * 2012-05-25 2012-10-10 北京交通大学 Rolling bearing failure diagnostic method based on multi-characteristic parameter
CN102778354A (en) * 2012-05-08 2012-11-14 南京工业大学 Method for diagnosing failure of wind-powered rotary support based on wavelet analysis
CN104215395A (en) * 2014-09-09 2014-12-17 中国石油大学(北京) Method and device for detecting imbalance fault of rotor
CN104655423A (en) * 2013-11-19 2015-05-27 北京交通大学 Rolling bearing fault diagnosis method based on time-frequency domain multidimensional vibration feature fusion

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102778354A (en) * 2012-05-08 2012-11-14 南京工业大学 Method for diagnosing failure of wind-powered rotary support based on wavelet analysis
CN102721545A (en) * 2012-05-25 2012-10-10 北京交通大学 Rolling bearing failure diagnostic method based on multi-characteristic parameter
CN104655423A (en) * 2013-11-19 2015-05-27 北京交通大学 Rolling bearing fault diagnosis method based on time-frequency domain multidimensional vibration feature fusion
CN104215395A (en) * 2014-09-09 2014-12-17 中国石油大学(北京) Method and device for detecting imbalance fault of rotor

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HONGTAOXUE,ET AL.: "Automatic diagnosis method for structural fault of rotating machinery based on distinctive frequency components and support vector machines under varied operating conditions", 《NEUROCOMPUTING》 *
赵利华: "旋转机械故障机理与轴心轨迹识别方法研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109488630A (en) * 2018-11-13 2019-03-19 上海金艺检测技术有限公司 Centrifugal blower rotor misalignment method for diagnosing faults based on harmonic wave relative indicatrix
CN110553844A (en) * 2019-07-24 2019-12-10 西安交通大学 Method and system for detecting misalignment fault of rotary machine

Similar Documents

Publication Publication Date Title
CN103575523B (en) The rotary machinery fault diagnosis method of kurtosis-envelope spectrum analysis is composed based on FastICA-
CN109488630B (en) Centrifugal fan rotor misalignment fault diagnosis method based on harmonic relative index
CN104165925B (en) The centrifugal compressor half-opened impeller crack fault detection method of accidental resonance
CN104343709B (en) A kind of device and method detecting fan trouble
Saidi et al. The use of SESK as a trend parameter for localized bearing fault diagnosis in induction machines
CN107313816B (en) A kind of T Steam Turbine SI measuring device and measuring method
CN106441896A (en) Characteristic vector extraction method for rolling bearing fault mode identification and state monitoring
CN107727338A (en) A kind of bridge damnification diagnostic method based on Vehicle-Bridge Coupling System
CN107704695A (en) A kind of full-scale quantitative Diagnosis method of housing washer defect
CN104048595A (en) Angle misalignment fault quantitative detection system and method for rotating machine
WO2015072130A1 (en) Leakage determination system and leakage determination method
CN104297003A (en) Fault monitoring method of bogie rotating part based on dynamic alarm threshold values
CN105675274A (en) Time-domain parameter and D-S evidence theory-based rotor running state monitoring method
CN104792519A (en) Offline automobile gearbox NVH (noise, vibration and harshness) detecting method
CN109211390B (en) Safety test and evaluation method for vibration and strength of power transmission tower
CN104215395A (en) Method and device for detecting imbalance fault of rotor
CN104483118A (en) Rotor dynamic and static rub impact fault diagnosis method based on instantaneous frequency shaft centerline orbit
CN104764622A (en) Bridge state detection device and detection method
CN104008294A (en) System and method for detecting abnormality of bearing
CN107478729A (en) The acoustic emission detection method of fluid machinery blade multiple crackses
CN103149869A (en) Systems and methods for use in monitoring an industrial facility
CN105510010A (en) Characteristic parameter model for rotating machinery in misaligned structural anomaly state
Jiang et al. Fault detection of rolling element bearing based on principal component analysis
CN105509965A (en) Characteristic parameter model for rotating machinery in unbalanced structural anomaly state
CN103335842A (en) Bearing fault recognition system based on probability density curve and method thereof

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160420

WD01 Invention patent application deemed withdrawn after publication