CN104614161B - A kind of rotating machinery rotatable parts come off weight and location recognition method - Google Patents
A kind of rotating machinery rotatable parts come off weight and location recognition method Download PDFInfo
- Publication number
- CN104614161B CN104614161B CN201510013858.6A CN201510013858A CN104614161B CN 104614161 B CN104614161 B CN 104614161B CN 201510013858 A CN201510013858 A CN 201510013858A CN 104614161 B CN104614161 B CN 104614161B
- Authority
- CN
- China
- Prior art keywords
- measuring point
- rotor
- impeller
- oil
- 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
- 238000000034 method Methods 0.000 title claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 230000035945 sensitivity Effects 0.000 claims abstract description 19
- 206010064930 age-related macular degeneration Diseases 0.000 claims abstract description 5
- 239000003921 oil Substances 0.000 claims description 21
- 238000013016 damping Methods 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 239000010687 lubricating oil Substances 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Structures Of Non-Positive Displacement Pumps (AREA)
- Testing Of Balance (AREA)
Abstract
Come off weight and location recognition method the invention discloses a kind of rotating machinery rotatable parts, oil film bearingses rigidity and damped coefficient are calculated by ARMD softwares, rotary machine rotor bearing arrangement kinetic equation is set up using finite element method.Apply mass unbalance power respectively at rotor impellers at different levels, obtain being vibrated at two ends of rotor bearing the vibration variable quantity of measuring point, each vibration measuring point is obtained to impeller unbalance response sensitivity at different levels using vector calculus rule.Variable quantity is vibrated according to measured rotor bear vibration measuring point, linkage disequilibrium response sensitivity matrix is that can obtain weight and the position that rotary part comes off.
Description
Technical field
Come off position and weight recognition methods the present invention relates to a kind of rotating machinery rotatable parts, help technical staff's analysis
Rotatable parts come off on rotor position and the size of the weight that comes off.Main application includes：Power, metallurgy, petrochemical industry, aviation etc. are big
Type rotating machinery, such as steam turbine, compressor, gas turbine, blower fan.
Background technology
Vibration is the important of the large high-speed rotating machinery safe operations such as influence steam turbine, compressor, gas turbine, blower fan
Factor.In work process, movable vane piece is rotated with high speed rotor, and its condition of work is extremely complex, except because of high-speed rotation
Acted on steam flow and born outside static stress and dynamic stress higher, being also in high temperature super heated steam area, two-phase respectively because of it flows through
Cross work in area and wet-steam region and bear high temperature, corrosion and erosion effect.This is resulted in running, and movable vane piece is fragile
Fracture.
The content of the invention
The technical problems to be solved by the invention be to provide a kind of operating rotatable parts of rotating machinery come off weight and
Location recognition method, the present invention can quickly, accurately, easily analyze the weight of releasing part and position on rotor, improve event
Barrier recognition efficiency.
Weight and location recognition side in order to solve the above technical problems, a kind of rotating machinery rotatable parts of present invention offer come off
Method, it is characterised in that comprise the following steps：
1) according to rotating machinery sliding bearing geometric parameter and lubricating oil physical parameter, sliding axle is calculated using ARMD softwares
Hold the stiffness coefficient and damped coefficient of oil film；
2) it is rotary machine rotor is discrete for impeller units and shaft part unit, set up rotor-support-foundation system using finite element method and transport
Dynamic equation is as follows：
Formula (1)
Wherein [M1]、[J1]、[K1] and Ω difference rotors lumped mass matrix, revolution matrix, stiffness matrix and anglec of rotation speed
Degree, U1、U2It is modal displacement vector, Q1、Q2It is broad sense Unbalanced force vector；It is U1And U2Second dervative,It is U1And U2First derivative；
3) using oil film rigidity coefficient and damped coefficient as the generalized force acted at axle journal node, rotor-support-foundation system fortune is substituted into
Dynamic equation, obtains the sliding bearing support rotor-support-foundation system equation of motion：
Formula (2)
In formula,
kxx、kyyIt is oil film bearingses Horizontal Stiffness Coefficient and oil film bearingses vertical stiffness coefficient, kxyRepresent axle journal
The horizontal direction oil-film force increment that vertical direction unit displacement causes, kyxRepresent that axle journal horizontal direction unit displacement causes vertical
Direction oil-film force increment., cxx、cyyIt is oil film bearingses horizontal damping coefficient and oil film bearingses vertical damping coefficient,
cxyRepresent the horizontal direction oil-film force increment that axle journal vertical direction unit speed causes, cyxRepresent axle journal horizontal direction unit speed
The vertical direction oil-film force increment that degree causes；
4) assume to have n vibration measuring point in rotating shaft, out-of-balance force is applied first at the 1st grade of impellerSolution formula (2)
The middle rotor-support-foundation system equation of motion, draw vibrated on rotor measuring point 1, measuring point 2 ..., the unbalance response at measuring point n is respectivelyThen unbalance response sensitivity of the n measuring point to out-of-balance force at first order impeller on rotor
Vector is：
Similarly, out-of-balance force is applied at i-stage impellerI represents any one-level, solves rotor-support-foundation system fortune in formula (2)
Dynamic equation, draw vibrated on rotor measuring point 1, measuring point 2 ..., the unbalance response at measuring point n is respectivelyThen unbalance response sensitivity of n, the two ends measuring point to out-of-balance force at i-stage impeller is vectorial
For：
Assuming that m grades of impeller is had, the unbalance response sensitivity coefficient combination by n measuring point to out-of-balance force at each impeller
It is following matrix form,
Formula (3)
(5) after rotating machinery rotatable parts come off, the vibration variable quantity of n measuring point of two ends of rotor is respectively：Δvib1、
Δvib2、……、Δvibn, the matrix form (3) of the unbalance response sensitivity obtained according to step (4), then releasing part weight
Vector expression with position is：
Formula (4).
Δ m is a vector, there is direction and a size, and the absolute value of Δ m is the size of quality of coming off, and the direction of vector is exactly
The position for coming off.
The beneficial effect that the present invention is reached：
Compared with prior art, the beneficial effect that is reached of the present invention is：
1) existing rotator model is utilized, equation of rotor motion can be drawn with FInite Element, surveyed so as to obtain each vibration
Point saves as data file to the unbalance response sensitivity matrix of impellers at different levels, is called at any time with facilitating.
2) according to vibration measuring point vibration amplitude and SPA sudden phase anomalies amount, rotatable parts can conveniently be obtained using sensitivity matrix
Weight and position are flown off, fault diagnosis efficiency is improved.
3) fixed algorithm can be formed, it is convenient and swift accurately to judge that rotating machinery rotatable parts fly off weight and position
Put.
Brief description of the drawings
Fig. 1 is model schematic of the invention.
Specific embodiment
The present invention is described further below in conjunction with the accompanying drawings.
Rotating machinery rotatable parts of the invention come off weight and location recognition method, it is characterised in that including following step
Suddenly：
1) according to rotating machinery sliding bearing geometric parameter and lubricating oil physical parameter, sliding axle is calculated using ARMD softwares
Hold the stiffness coefficient and damped coefficient of oil film；
2) it is rotary machine rotor is discrete for impeller units and shaft part unit, set up rotor-support-foundation system using finite element method and transport
Dynamic equation is as follows：
Formula (1)
Wherein [M1]、[J1]、[K1] and Ω difference rotors lumped mass matrix, revolution matrix, stiffness matrix and anglec of rotation speed
Degree, U1、U2It is modal displacement vector, Q1、Q2It is broad sense Unbalanced force vector；It is U1And U2Second dervative,It is U1And U2First derivative；
3) using oil film rigidity coefficient and damped coefficient as the generalized force acted at axle journal node, rotor-support-foundation system fortune is substituted into
Dynamic equation, obtains the sliding bearing support rotor-support-foundation system equation of motion：
Formula (2)
In formula,
kxx、kyyIt is oil film bearingses Horizontal Stiffness Coefficient and oil film bearingses vertical stiffness coefficient, kxy、kyxIt is oil
Film intersects stiffness coefficient, cxx、cyyIt is oil film bearingses horizontal damping coefficient and oil film bearingses vertical damping coefficient,
cxy、cyxIt is oil film cross damping coefficient；
4) assume to have n vibration measuring point in rotating shaft, out-of-balance force is applied first at the 1st grade of impellerSolution formula (2)
The middle rotor-support-foundation system equation of motion, draw vibrated on rotor measuring point 1, measuring point 2 ..., the unbalance response at measuring point n is respectivelyThen unbalance response sensitivity of the n measuring point to out-of-balance force at first order impeller on rotor
Vector is：
Similarly, out-of-balance force is applied at i-stage impellerI represents any one-level, solves rotor-support-foundation system fortune in formula (2)
Dynamic equation, draw vibrated on rotor measuring point 1, measuring point 2 ..., the unbalance response at measuring point n is respectivelyThen unbalance response sensitivity of n, the two ends measuring point to out-of-balance force at i-stage impeller is vectorial
For：
Assuming that m grades of impeller is had, the unbalance response sensitivity coefficient combination by n measuring point to out-of-balance force at each impeller
It is following matrix form,
Formula (3)
(5) after rotating machinery rotatable parts come off, the vibration variable quantity of n measuring point of two ends of rotor is respectively：Δvib1、
Δvib2、……、Δvibn, the matrix form (3) of the unbalance response sensitivity obtained according to step (4), then releasing part weight
Vector expression with position is：
Formula (4).
Δ m is a vector, there is direction and a size, and the absolute value of Δ m is the size of quality of coming off, and the direction of vector is exactly
The position for coming off.
As a example by model as shown in Figure 1, the rotor includes level Four impeller.
Rotating machinery rotatable parts of the invention come off weight and location recognition method, comprise the following steps：
(1) bearing diameter 300mm, 250mm long, bearing clearance 0.33mm, calculate the bearing and exist using ARMD softwares
Oil film rigidity and damping parameter in the case of 3000 revs/min；
(2) it is 25 shaft parts and impeller pattern by the rotor modelling；
(3) rotor-bearing system motion finite element equation is set up；
(4) according to a large amount of engineering experiences, in actual motion, rotary part comes off and mostly occurs on blade, shroud, so only
Demand goes out the sensitivity that unit mass at impeller section comes off to bearings at both ends horizontal and vertical measuring point vibratory response.
(5) when amount of unbalance is 0 ° of 1kg ∠ at each impeller, the imbalance at the unit bearings at both ends horizontal and vertical measuring point
Response sensitivity matrix α is：
(6) assume that the unit third level impeller comes off quality 2.7kg, the two end axles for now causing in circumferential 190 degree of positions
Both horizontally and vertically vibration variable quantity is held to be respectively：
The bearing measuring point of table 1 vibrates variable quantity (amplitude μm ∠ angles °)
| Measuring point | 1 level | 1 is vertical | 2 levels | 2 is vertical |
| Fundamental vibration variation delta vib | 113.4∠337 | 126.9∠350 | 173∠39 | 164.3∠57 |
Each measuring point vibration values variable quantity is expressed as plural form：
Then, the come off complex expression Δ m of quality and position of the unit is：
By above-mentioned complex representation for the form of amplitude and angle is as follows：
As can be seen that identify that third level impeller comes off 2.7kg 190 degree of positions using sensitivity matrix, with the weight that comes off
Amount is consistent with circumferential position.Below the present invention is disclosed with preferred embodiment, so it is not intended to limiting the invention, all uses
The technical scheme that equivalent or equivalent transformation mode are obtained, is within the scope of the present invention.
Claims (1)
1. a kind of rotating machinery rotatable parts come off weight and location recognition method, it is characterised in that comprise the following steps：
1) according to rotating machinery sliding bearing geometric parameter and lubricating oil physical parameter, sliding bearing oil is calculated using ARMD softwares
The stiffness coefficient and damped coefficient of film；
2) it is rotary machine rotor is discrete for impeller units and shaft part unit, set up rotor-support-foundation system motion side using finite element method
Journey is as follows：
Wherein [M1]、[J1]、[K1] and Ω be respectively rotor lumped mass matrix, revolution matrix, stiffness matrix and the anglec of rotation speed
Degree, U1、U2It is modal displacement vector, Q1、Q2It is broad sense Unbalanced force vector；It is U1And U2Second dervative,It is U1And U2First derivative；
3) oil film rigidity coefficient and damped coefficient are substituted into rotor-support-foundation system motion side as the generalized force acted at axle journal node
Journey, obtains the sliding bearing support rotor-support-foundation system equation of motion：
In formula,
kxx、kyyIt is oil film bearingses Horizontal Stiffness Coefficient and oil film bearingses vertical stiffness coefficient, kxyRepresent that axle journal is vertical
The horizontal direction oil-film force increment that direction unit displacement causes, kyxRepresent the vertical direction that axle journal horizontal direction unit displacement causes
Oil-film force increment；cxx、cyyIt is oil film bearingses horizontal damping coefficient and oil film bearingses vertical damping coefficient, cxyRepresent
The horizontal direction oil-film force increment that axle journal vertical direction unit speed causes, cyxRepresent what axle journal horizontal direction unit speed caused
Vertical direction oil-film force increment；
4) assume to have n vibration measuring point in rotating shaft, out-of-balance force is applied first at the 1st grade of impellerSolution formula (2) transfer
The subsystem equation of motion, draw vibrated on rotor measuring point 1, measuring point 2 ..., the unbalance response at measuring point n is respectivelyThen unbalance response sensitivity of the n measuring point to out-of-balance force at first order impeller on rotor
Vector is：
Similarly, out-of-balance force is applied at i-stage impellerI represents any one-level, solves rotor-support-foundation system motion side in formula (2)
Journey, draw vibrated on rotor measuring point 1, measuring point 2 ..., the unbalance response at measuring point n is respectivelyThen unbalance response sensitivity of n, the two ends measuring point to out-of-balance force at i-stage impeller is vectorial
For：
Assuming that have m grades of impeller, by n measuring point the unbalance response sensitivity coefficient of out-of-balance force at each impeller is combined as
Lower matrix form,
(5) after rotating machinery rotatable parts come off, the vibration variable quantity of n measuring point of two ends of rotor is respectively：△vib1、△
vib2、……、△vibn, the matrix form (3) of the unbalance response sensitivity obtained according to step (4), then releasing part weight and
The vector expression of position is：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510013858.6A CN104614161B (en) | 2015-01-12 | 2015-01-12 | A kind of rotating machinery rotatable parts come off weight and location recognition method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510013858.6A CN104614161B (en) | 2015-01-12 | 2015-01-12 | A kind of rotating machinery rotatable parts come off weight and location recognition method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104614161A CN104614161A (en) | 2015-05-13 |
| CN104614161B true CN104614161B (en) | 2017-06-23 |
Family
ID=53148708
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510013858.6A Expired - Fee Related CN104614161B (en) | 2015-01-12 | 2015-01-12 | A kind of rotating machinery rotatable parts come off weight and location recognition method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104614161B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110929419B (en) * | 2018-12-29 | 2021-08-13 | 山东大学 | Method for quickly predicting instability limit of steam turbine rotor system based on shroud zero damping |
| CN110261112B (en) * | 2019-07-11 | 2021-03-23 | 中国航发哈尔滨东安发动机有限公司 | Dynamic balance evaluation method for ultra-precise journal abrasion condition |
| CN114323608A (en) * | 2021-12-20 | 2022-04-12 | 东方电气集团东方汽轮机有限公司 | Method for judging axial position of shedding part of steam turbine rotor |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4193305A (en) * | 1978-08-16 | 1980-03-18 | Lee Hunter | Vehicle wheel and tire balance testing machine |
| CN101059386A (en) * | 2007-06-04 | 2007-10-24 | 西安交通大学 | Precession vector based holographic on-spot dynamic balance method |
| CN102494847A (en) * | 2011-12-12 | 2012-06-13 | 山东电力研究院 | Site high-speed dynamic balance method for large steam turbine generator unit shafting |
| CN103592081A (en) * | 2013-11-01 | 2014-02-19 | 西安电子科技大学 | Rotor non-trial-weight dynamic balancing method suitable for distributed unbalance |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8056411B2 (en) * | 2009-03-02 | 2011-11-15 | GM Global Technology Operations LLC | Shaft balancing system and methods |
-
2015
- 2015-01-12 CN CN201510013858.6A patent/CN104614161B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4193305A (en) * | 1978-08-16 | 1980-03-18 | Lee Hunter | Vehicle wheel and tire balance testing machine |
| CN101059386A (en) * | 2007-06-04 | 2007-10-24 | 西安交通大学 | Precession vector based holographic on-spot dynamic balance method |
| CN102494847A (en) * | 2011-12-12 | 2012-06-13 | 山东电力研究院 | Site high-speed dynamic balance method for large steam turbine generator unit shafting |
| CN103592081A (en) * | 2013-11-01 | 2014-02-19 | 西安电子科技大学 | Rotor non-trial-weight dynamic balancing method suitable for distributed unbalance |
Non-Patent Citations (3)
| Title |
|---|
| 改进传递矩阵法计算转子***不平衡响应和灵敏度;杨建刚等;《机械工程学报》;20010630(第06期);第109-130页 * |
| 汽轮机叶片脱落故障诊断方法的研究;晋风华;《热力透平》;20070931(第03期);第178-180、190页 * |
| 转子***不平衡响应的灵敏度和估算方法;高建民等;《振动工程学报》;19910930(第03期);第83-87页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104614161A (en) | 2015-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tian et al. | Dynamic behaviours of a full floating ring bearing supported turbocharger rotor with engine excitation | |
| Ma et al. | Effects of eccentric phase difference between two discs on oil-film instability in a rotor–bearing system | |
| Xu et al. | Active magnetic bearings dynamic parameters identification from experimental rotor unbalance response | |
| CN104614161B (en) | A kind of rotating machinery rotatable parts come off weight and location recognition method | |
| Rashidi et al. | Bifurcation and nonlinear dynamic analysis of a rigid rotor supported by two-lobe noncircular gas-lubricated journal bearing system | |
| Ma et al. | Numerical and experimental analysis of the first-and second-mode instability in a rotor-bearing system | |
| Dupont | Robust rotor dynamics for high-speed air bearing spindles | |
| Jung et al. | Rotordynamic modelling and analysis of a radial inflow turbine rotor-bearing system | |
| CN108287952B (en) | Low-temperature liquid expander rotor critical rotation speed prediction method considering radial fluid force influence | |
| Fegade et al. | Unbalanced response of rotor using ansys parametric design for different bearings | |
| CN105488286B (en) | Vertical multi-stage centrifugal pump rotor dynamics modeling method under a kind of band liquid operating mode | |
| Alekseeva et al. | Evaluation of effect of oil film of rotor bearing | |
| Childs et al. | Rotordynamic performance of a negative-swirl brake for a tooth-on-stator labyrinth seal | |
| Basumatary et al. | A numerical study on effect of electromagnetic actuator on rigid rotor supported on gas foil bearing | |
| He et al. | Rotor dynamics of multistage centrifugal pump | |
| Miranda et al. | Finite element method applied to the eigenvalue analysis of flexible rotors supported by journal bearings | |
| CN113326647A (en) | Water turbine shafting rotor dynamic modal calculation method | |
| Jamil et al. | Study the dynamic behavior of rotor supported on a worn journal bearings | |
| Ogbonnaya et al. | Analysis of gas turbine blade vibration due to random excitation | |
| Wang et al. | Nonlinear vibration analysis for a Jeffcott rotor with seal and air-film bearing excitations. | |
| EP3006729A1 (en) | Rotating machine and installation for converting energy comprising such a machine | |
| Liu et al. | Research of vibration characteristics based on a crystal format model of rotor | |
| Tessarzik et al. | Experimental evaluation of multiplane-multispeed rotor balancing through multiple critical speeds | |
| Wang et al. | Dynamic analysis for a two-shaft rotor system under a sudden loss of mass excitation | |
| Pandya et al. | Vibration Analysis of Pump Shaft Using Finite Element Analysis Software: A Review |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190529 Address after: 211100 Tianyuan East Road 1009, Jiangning District High-tech Park, Nanjing City, Jiangsu Province Patentee after: Nanjing Yunqi Resonant Power Technology Co.,Ltd. Address before: 211167 1 Hong Kong Avenue, Jiangning Science Park, Nanjing, Jiangsu Patentee before: Nanjing Institute of Technology |
|
| TR01 | Transfer of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170623 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |