CN114061729A - Rotor sudden-change vibration fault positioning method and system - Google Patents
Rotor sudden-change vibration fault positioning method and system Download PDFInfo
- Publication number
- CN114061729A CN114061729A CN202111232415.8A CN202111232415A CN114061729A CN 114061729 A CN114061729 A CN 114061729A CN 202111232415 A CN202111232415 A CN 202111232415A CN 114061729 A CN114061729 A CN 114061729A
- Authority
- CN
- China
- Prior art keywords
- vibration
- rotor
- fault
- beta
- alpha
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000008859 change Effects 0.000 claims abstract description 22
- 239000013598 vector Substances 0.000 claims abstract description 20
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 3
- 238000003745 diagnosis Methods 0.000 description 5
- 238000012549 training Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000012706 support-vector machine Methods 0.000 description 3
- 238000005299 abrasion Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
- G01H1/006—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines of the rotor of turbo machines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a method and a system for positioning sudden change vibration fault of a rotor, which comprises the following steps: step 1, recording vibration vectors before sudden change of vibration at two ends of a rotor in the same directionAndabrupt change of vibration vectorAndstep 2, calculating vibration variationStep 3, the vibration variation quantity is adjustedAndharmonic component decompositionWhereinThe phase angle of (a) is alpha,has a phase angle of beta; step 4, if alpha-beta is approximately equal to 90 degrees or alpha-beta is approximately equal to-90 degrees, the vibration fault point is close to the side A; if it isIf the value is far greater than 1, the vibration variable is mainly the same direction component, and the corresponding fault point is located in the middle of the rotor. The invention adopts the harmonic component decomposition of the vibration variation before and after the rotor fault, thereby obtainingThe phase angle of (a) is alpha,the phase angle is beta, and the fault point is judged according to the phase angle, so that the whole calculation process is simple and the positioning is accurate.
Description
Technical Field
The invention relates to the technical field of mechanical vibration fault analysis of a steam turbine, in particular to a method and a system for positioning sudden-change vibration fault of a rotor.
Background
When the steam turbine takes place the vibration trouble, often can lead to some harmful results to produce, make the normal operation of steam turbine and production receive great influence, can lead to the steam turbine overall to take place the condition of destroying even. When some vibration faults occur in the unit, not only the nature of the fault needs to be judged, but also the fault point is located at the axial and circumferential positions of the rotor.
For diagnosis and positioning of some sudden-change vibrations, the method can not only provide auxiliary criteria for judging the nature of the fault, but also guide the direction of unit maintenance and avoid meaningless expansion of the maintenance range for positioning hidden faults such as turn-to-turn short circuit of the generator rotor. The prior art can detect that the rotor has faults and can not accurately judge fault points. The wavelet-gray level co-occurrence matrix-based engine rotor collision and abrasion fault diagnosis method disclosed in application number 202110668493.6 solves the problems of low fault identification precision, high requirements of a neural network method on samples and hardware and tedious training of the traditional method. The method comprises the following specific steps: 1) collecting engine friction data; 2) performing wavelet transformation on the engine rub-impact data to obtain a vibration signal wavelet transformation diagram; 3) converting the wavelet transform map of the vibration signal into a gray scale map, extracting the gray scale co-occurrence matrix image characteristic quantity of the gray scale map, and constructing a characteristic vector; 4) setting labels for the feature vectors, and dividing the feature vectors into a training set and a test set; 5) inputting the training set into a machine learning algorithm support vector machine for training to obtain a trained support vector machine; 6) and inputting the test set into a trained support vector machine to realize the collision and abrasion fault diagnosis of the engine rotor. Although this method can solve the problem of failure diagnosis, it cannot locate the failure position.
Disclosure of Invention
The invention aims to solve the technical problem of how to quickly locate the fault position of the rotor.
The invention solves the technical problems through the following technical means:
the method for positioning the sudden change vibration fault of the rotor comprises the steps that the symmetrical rotor comprises the rotor, a plurality of level blades symmetrically fixed on the rotor and supporting points A, B; the support points A, B respectively provide support for two ends of the rotor; assuming that a certain level of blade falls off, the method for judging the position of the fault point comprises the following steps:
step 1, recording two ends of the rotor to be the sameVibration vector before sudden change of vibration in one directionAndabrupt change of vibration vectorAnd
Step 3, the vibration variation quantity is adjustedAndharmonic component decompositionWhereinThe phase angle of (a) is alpha,has a phase angle of beta;
The invention carries out harmonic component on the vibration variation before and after the rotor faultIs decomposed to obtainThe phase angle of (a) is alpha,the phase angle is beta, and the fault point is judged according to the phase angle, so that the whole calculation process is simple and the positioning is accurate.
Further, the vibration in the other direction on the same rotor is calculated and judged according to the steps 1 to 4.
Further, in the step 4, when the angle α - β ≈ 90 °, the fault position isWherein the distance between the falling point and the A side weighted plane is LaThe axial distance of the weighted planes on both sides of the supporting point A, B is L,for the co-directional component influence coefficient,is an inverse component influence coefficient.
Further, in the step 4, when the angle alpha-beta is approximately equal to-90 degrees, the fault position is
Corresponding to the method, the invention also provides a rotor sudden-change vibration fault positioning system, wherein a symmetrical rotor comprises a rotor, a plurality of level blades symmetrically fixed on the rotor, and supporting points A, B; the support points A, B respectively provide support for two ends of the rotor; assuming that a blade of a certain level is detached, the system for determining the position of the fault point includes:
the rotor sudden change front and back vibration vector recording module is used for recording the vibration vectors of the two ends of the rotor before sudden change in the same directionAndabrupt change of vibration vectorAnd
Harmonic component decomposition module for the vibration variationAndharmonic component decomposition
if alpha-beta is approximately equal to 90 degrees or alpha-beta is approximately equal to-90 degrees, the fault point judgment module indicates that the vibration fault point is close to the side A; if it isIf the vibration variable is far more than 1, the vibration variable is mainly the same direction component, and the corresponding fault point is positioned in the middle part of the rotor。
Further, the vibration in the other direction on the same rotor is calculated and judged according to the steps 1 to 4.
Further, in the fault point judgment module, when the angle alpha-beta is approximately equal to 90 degrees, the fault position isWherein the distance between the falling point and the A side weighted plane is LaThe axial distance of the weighted planes on both sides of the supporting point A, B is L,for the co-directional component influence coefficient,is an inverse component influence coefficient.
Furthermore, in the fault point judgment module, when the angle alpha-beta is approximately equal to-90 degrees, the fault position is
The invention has the advantages that:
the invention belongs to a rotor vibration fault diagnosis and fault location method, which is suitable for rotary equipment provided with an online vibration monitoring systemThe phase angle of (a) is alpha,the phase angle is beta, the relative position of a fault point in the axial direction is judged according to the phase angle, the whole calculation process is simple, the structural fault is accurately positioned, and the position of a component with the fault of the rotor can be accurately judged by combining the structural characteristics of rotors of different devices.
Drawings
FIG. 1 is a schematic structural diagram of a symmetrical rotor applied to a method for locating a sudden rotor vibration fault in an embodiment of the present invention;
fig. 2 is a flow chart of a rotor sudden-change vibration fault locating method in the embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the embodiments of the present invention, and it is obvious that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment provides a method for positioning sudden change vibration fault of a rotor, as shown in fig. 1, wherein fig. 1 is a schematic structural diagram of a symmetrical rotor, and the symmetrical rotor comprises a rotor, blades and supporting bearing bushes a and B. For such a rotor and supporting point A, B system, the rotor has mechanism failure (e.g. some stage of blade has part falling off, the falling mass is m, the distance from the falling point to the A side weighted plane is LaA, B, the axial distance of the weighted planes on both sides is L), vibration changes, and as shown in FIG. 2, the method for determining the position of the fault point mainly comprises the following steps:
step 1, recording vibration vectors of two ends of a rotor in the same direction (X direction or Y direction) before sudden change of vibrationAndabrupt change of vibration vectorAnd
as shown in FIG. 1, the mass m causes an additional centrifugal force F in the direction opposite to the drop-off position, which force F is equivalent to that at the weight plane a according to the force translation principleAnd couple FLaIt can be considered that there are forces in the same magnitude and opposite directions in the plane bThus, in-phase component force is formed in the planes a and bAnd couple of opposite phase components(couple FLa withIn the opposite direction).
According to the harmonic component method, the rotor vibration at the working rotation speed is decomposed into in-phase (symmetrical) component vibration and anti-phase (anti-symmetrical) component vibration, the in-phase component is caused by first-order unbalance of the rotor, the anti-phase component vibration is caused by second-order unbalance (couple unbalance) of the rotor and accords with the orthogonal relation, and therefore in-phase weighting mass and anti-phase weighting mass are obtained.
Reverse weighting qualityWherein,for the co-directional component influence coefficient,is an inverse component influence coefficient;
according to the couple balance mode, if the couple at the weighting surfaces a and b and the drop-off position m is 0, the couple is equal to
If alpha-beta is approximately equal to-90 degrees, the vibration fault point is close to the side A, and the corresponding fault position is
If it isIf the value is far greater than 1, the vibration variable is mainly the same direction component, and the corresponding fault point is located in the middle of the rotor.
And calculating and judging the vibration in the other direction on the same bearing according to the steps 1 to 4.
Corresponding to the method, the embodiment also provides a rotor sudden-change vibration fault positioning system, which comprises
The rotor sudden change front and back vibration vector recording module is used for recording the vibration vectors of the two ends of the rotor before sudden change in the same directionAndabrupt change of vibration vectorAnd
Harmonic component decomposition module for the vibration variationAndharmonic component decomposition
as shown in FIG. 1, the mass m causes an additional centrifugal force F in the direction opposite to the drop-off position, which force F is equivalent to that at the weight plane a according to the force translation principleAnd couple FLaIt can be considered that there are forces in the same magnitude and opposite directions in the plane bThus, in-phase component force is formed in the planes a and bAnd couple of opposite phase components(couple FLa withIn the opposite direction).
According to the harmonic component method, the rotor vibration at the working rotation speed is decomposed into in-phase (symmetrical) component vibration and anti-phase (anti-symmetrical) component vibration, the in-phase component is caused by first-order unbalance of the rotor, the anti-phase component vibration is caused by second-order unbalance (couple unbalance) of the rotor and accords with the orthogonal relation, and therefore in-phase weighting mass and anti-phase weighting mass are obtained.
Reverse weighting qualityWherein,for the co-directional component influence coefficient,is an inverse component influence coefficient;
according to the couple balance mode, if the couple at the weighting surfaces a and b and the drop-off position m is 0, the couple is equal to
If alpha-beta is approximately equal to 90 degrees or alpha-beta is approximately equal to-90 degrees, the fault point judgment module indicates that the vibration fault point is close to the side A; if it isIf the value is far greater than 1, the vibration variable is mainly the same direction component, and the corresponding fault point is located in the middle of the rotor.
The above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (8)
1. The method for positioning the sudden change vibration fault of the rotor is characterized in that the symmetrical rotor comprises the rotor, a plurality of grade blades symmetrically fixed on the rotor and supporting points A, B; the support points A, B respectively provide support for two ends of the rotor; assuming that a certain level of blade falls off, the method for judging the position of the fault point comprises the following steps:
step 1, recording vibration vectors before sudden change of vibration at two ends of a rotor in the same directionAndabrupt change of vibration vectorAnd
Step 3, the vibration variation quantity is adjustedAndharmonic component decompositionWhereinThe phase angle of (a) is alpha,has a phase angle of beta;
step 4, if alpha-beta is approximately equal to 90 degrees or alpha-beta is approximately equal to-90 degrees, the vibration fault point is close to the side A; if it isIf the value is far greater than 1, the vibration variable is mainly the same direction component, and the corresponding fault point is located in the middle of the rotor.
2. The method for locating the sudden rotor vibration fault according to claim 1, wherein the vibration in the other direction on the same rotor is calculated and judged according to steps 1 to 4.
3. The method for locating the sudden rotor vibration fault according to claim 1 or 2, wherein in the step 4, when the angle α - β ≈ 90 °, the fault position isWherein the distance between the falling point and the A side weighted plane is LaThe axial distance of the weighted planes on both sides of the supporting point A, B is L,for the co-directional component influence coefficient,is an inverse component influence coefficient.
5. The rotor sudden change vibration fault location system of claim 1, wherein a symmetrical rotor comprises a rotor, a plurality of stage blades symmetrically fixed on the rotor, a support point A, B; the support points A, B respectively provide support for two ends of the rotor; assuming that a blade of a certain level is detached, the system for determining the position of the fault point includes:
the rotor sudden change front and back vibration vector recording module is used for recording the vibration vectors of the two ends of the rotor before sudden change in the same directionAndabrupt change of vibration vectorAnd
Harmonic component decomposition module for the vibration variationAndharmonic component decomposition
if alpha-beta is approximately equal to 90 degrees or alpha-beta is approximately equal to-90 degrees, the fault point judgment module indicates that the vibration fault point is close to the side A; if it isIf the value is far greater than 1, the vibration variable is mainly the same direction component, and the corresponding fault point is located in the middle of the rotor.
6. The system for locating sudden rotor vibration fault according to claim 5, wherein vibration in another direction on the same rotor is calculated and judged according to steps 1-4.
7. The rotor sudden-change vibration fault location system according to claim 5 or 6, wherein in the fault point judgment module, when α - β ≈ 90 °, the fault position isWherein the distance between the falling point and the A side weighted plane is LaThe axial distance of the weighted planes on both sides of the supporting point A, B is L,for the co-directional component influence coefficient,is an inverse component influence coefficient.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111232415.8A CN114061729A (en) | 2021-10-22 | 2021-10-22 | Rotor sudden-change vibration fault positioning method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111232415.8A CN114061729A (en) | 2021-10-22 | 2021-10-22 | Rotor sudden-change vibration fault positioning method and system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114061729A true CN114061729A (en) | 2022-02-18 |
Family
ID=80235328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111232415.8A Withdrawn CN114061729A (en) | 2021-10-22 | 2021-10-22 | Rotor sudden-change vibration fault positioning method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114061729A (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095561A (en) * | 2010-12-01 | 2011-06-15 | 浙江省电力试验研究院 | Falling fault positioning method for rotating part of large-size steam turbine |
US20140230555A1 (en) * | 2012-12-20 | 2014-08-21 | Zapadoceska Univerzita V Plzni | Method of detecting and localizing partial rotor-stator rubbing during the operation of a turbine |
CN110231161A (en) * | 2019-07-12 | 2019-09-13 | 中国大唐集团科学技术研究院有限公司华东电力试验研究院 | Single-span rotor fault position diagnostic method and system based on three-dimensional system of coordinate |
-
2021
- 2021-10-22 CN CN202111232415.8A patent/CN114061729A/en not_active Withdrawn
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095561A (en) * | 2010-12-01 | 2011-06-15 | 浙江省电力试验研究院 | Falling fault positioning method for rotating part of large-size steam turbine |
US20140230555A1 (en) * | 2012-12-20 | 2014-08-21 | Zapadoceska Univerzita V Plzni | Method of detecting and localizing partial rotor-stator rubbing during the operation of a turbine |
CN110231161A (en) * | 2019-07-12 | 2019-09-13 | 中国大唐集团科学技术研究院有限公司华东电力试验研究院 | Single-span rotor fault position diagnostic method and system based on three-dimensional system of coordinate |
Non-Patent Citations (2)
Title |
---|
刘宝富 等: "基于TDM的汽轮发电机组轴承振动异常分析", 《沈阳工程学院学报(自然科学版)》 * |
应光耀 等: "基于谐分量法的汽轮机叶片飞脱故障定位方法研究及应用", 《动力工程学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Vibration signal analysis using parameterized time–frequency method for features extraction of varying-speed rotary machinery | |
CN109253882A (en) | A kind of rotor crack fault diagnostic method based on variation mode decomposition and gray level co-occurrence matrixes | |
Pinheiro et al. | Vibration analysis in turbomachines using machine learning techniques | |
CN110905732B (en) | Method and system for identifying unbalance of wind wheel of wind turbine generator and storage medium | |
CN110146279B (en) | Ship shafting unbalance fault diagnosis method based on vector evidence reasoning | |
CN109211546A (en) | Rotary machinery fault diagnosis method based on noise reduction autocoder and incremental learning | |
CN109115408B (en) | Dynamic balance test method of large hydroelectric generating set based on centrifugal force equation | |
CN106762452A (en) | Fan master control system fault diagnosis and on-line monitoring method based on data-driven | |
CN117077064B (en) | Fault detection method, device and equipment for wind power equipment and storage medium | |
CN105352726B (en) | A kind of method for diagnosing faults of gear | |
CN110646138B (en) | Rotary machine keyless phase non-test weight dynamic balance method and analysis device | |
CN117113200B (en) | Rotor fault diagnosis method, device, electronic equipment and medium | |
CN114061729A (en) | Rotor sudden-change vibration fault positioning method and system | |
WO2022064038A1 (en) | Method and system for wind speed determination using vibration data | |
JP2020064024A (en) | Rotary machine abnormality factor specifying device | |
van Vondelen et al. | Modal analysis of an operational offshore wind turbine using enhanced Kalman filter‐based subspace identification | |
Barszcz et al. | Concept of automated malfunction detection of large turbomachinery using machine learning on transient data | |
Michalski et al. | Applying Kalman filtering to unbalance estimation in rotating machinery | |
CN110231161B (en) | Single-span rotor fault position diagnosis method and system based on three-dimensional coordinate system | |
KR101507248B1 (en) | Method for determine eigen value and apparatus for using the same | |
CN113029620A (en) | Method for predicting vibration response of shaft-disk-blade non-axisymmetric rotating machine | |
Perišic et al. | Gearbox fatigue load estimation for condition monitoring of wind turbines | |
Craggs | A component mode method for modelling the dynamics of turbo-generator sets | |
Wang et al. | Conductor galloping prediction on imbalanced datasets: Svm with smart sampling | |
Pennacchi et al. | Unbalance identification in large steam turbo-generator unit using a model-based method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220218 |