CN115014637B - Modal dynamic balance method based on low-rotation-speed measurement - Google Patents
Modal dynamic balance method based on low-rotation-speed measurement Download PDFInfo
- Publication number
- CN115014637B CN115014637B CN202210411949.5A CN202210411949A CN115014637B CN 115014637 B CN115014637 B CN 115014637B CN 202210411949 A CN202210411949 A CN 202210411949A CN 115014637 B CN115014637 B CN 115014637B
- Authority
- CN
- China
- Prior art keywords
- rotor
- order
- modal
- balance
- speed
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005259 measurement Methods 0.000 title claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 15
- 230000001360 synchronised effect Effects 0.000 claims abstract description 5
- 238000013016 damping Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000013459 approach Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
- G01M1/14—Determining imbalance
- G01M1/16—Determining imbalance by oscillating or rotating the body to be tested
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Balance (AREA)
Abstract
A modal dynamic balance method based on low-rotation-speed measurement relates to the field of rotary machinery and vibration test. In order to eliminate unbalanced response of the former N-order modes, test weights are added on N planes, unbalanced calibration quantities to be added are solved, N balance planes are taken for balancing the former N-order modes, synchronous vibration response of measuring points is obtained when a rotor runs at any rotating speed, the former N+M-order modes are obtained at a lower rotating speed, different rotating speeds are measured in a range lower than a first-order critical rotating speed of the rotor, known test weights are added on the first balance planes, influence coefficients are obtained after the measurement of different rotating speeds are also carried out, and unbalanced calibration quantities are solved repeatedly for the rest balance planes. The rotor can be measured at a lower rotating speed, and the balance parameters are fitted by combining with the rotor modal parameter information, so that modal dynamic balance at a low speed is realized, resonance caused by the supercritical rotating speed in the testing process is avoided, and the rotor is ensured to have a more reliable balance state in a wider rotating speed range.
Description
Technical Field
The invention relates to the field of rotary machinery and vibration testing, in particular to a modal dynamic balance method based on low-rotation-speed measurement, which can be applied to dynamic balance of a flexible rotor.
Background
Industrial rotating machinery such as aeroengines and large gas turbines are high in running speed, rotors of the industrial rotating machinery often need to work in a supercritical speed state, and in the process of passing through the supercritical speed, the rotors of the industrial rotating machinery often cause remarkable vibration and deformation due to structural resonance due to unbalanced rotor mass, so that the running performance and safety of equipment are seriously affected.
In order to ensure that the rotor system can pass through the critical rotation speed more stably, the conventional modal dynamic balance method generally requires dynamic balance test to be carried out near the critical rotation speed, but has the problem that vibration is too large when approaching or passing through the critical rotation speed in the measurement process, so that the dynamic balance test cost is higher, and a larger potential safety hazard exists. Therefore, a new dynamic balance method which can be safer and more reliable is urgently needed.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a modal dynamic balance method based on low-rotation-speed measurement, so that a rotor can be measured at a lower rotation speed, resonance is avoided in the test process as much as possible, modal balance is realized, and a rotor is ensured to have a more reliable balance state in a wider rotation speed range.
The specific steps of the invention are as follows:
1) For a flexible rotor system, it is in the axial direction x s The synchronous vibration response at the point due to mass imbalance is as follows:
wherein ,mr ,ω r ,ξ r ,ψ r (x) Respectively representing the r-th order modal mass, modal frequency, modal damping and modal shape, wherein omega is the shaft rotation speed, U (x) is the unbalanced distribution of the rotor, and l is the length of the rotor; typically, as the rotational speed approaches the r-th order modal frequency, the order modal response dominates, written as:
to eliminate the unbalanced response of the first N-order modes, the test weight lambda needs to be added on N planes k K=1, 2,..n, such thatThe method comprises the following steps:
in the formula ,representing the integral of the modal unbalance component of the r-th order over the entire axis, which may be referred to as the modal unbalance factor of the r-th order, then equation (3) is written in the form of a matrix as: />
The required added unbalance calibration amount can be solved according to the above formula:
Λ=-[Ψ] -1 u (5)
wherein Λ= [ λ ] 1 ,λ 2 ,...λ N ] T Represents the set of unbalance calibration quantities that need to be added, [ ψ ]] r =[ψ r (x 1 ),ψ r (x 2 ),...,ψ r (x N )]Represents the r-order vibration mode vector, u= [ u ] 1 ,u 2 ,...u N ] T Representing a vector of the first N-order modal unbalance factors.
2) Taking N balance surfaces to balance the front N-order modes, and setting a measurement point to be positioned at x s Where each order mode shape is ψ i (x s ) Each row [ ψ ] of the vibration mode matrix] i Multiplying byThen write as:
when the rotor runs at any rotation speed, the rotor is at a measuring point x s The synchronous vibration response of (2) is:
wherein ,defined as a complex parameter term related to rotational speed, when considering the balanced pre-N-order mode, consider that the mode after the n+m-th order contributes negligible to the pre-N-order, at lower rotational speeds Ω (Ω < ω) 1 <ω N <<ω N+M+1 ) Then, the mode is truncated to the former N+M order mode to obtain:
when the different rotational speeds Ω= { Ω are performed in the first-order critical rotational speed range below the rotor 1 ,Ω 2 ,...Ω p Measurement of } (p.gtoreq.N+M), then the above formula is written as:
for simplicity and convenience, recordZ=[Z 1 (Ω) Z 2 (Ω)...Z N+M (Ω)]When the critical rotation speed omega of each order of the rotor is known i (i=1, 2,) n+m and critical damping ratio ζ i (i=1, 2,..n+m), then solved by the formula:
in the formula ,(·)+ Representing the pseudo-inverse of the matrix, y= [ Y (x s ,Ω 1 ),Y(x s ,Ω 2 ),...,Y(x s ,Ω p )] T 。
Adding a known test weight Q to the first equilibrium surface 1 The same applies to different rotational speeds Ω= { Ω 1 ,Ω 2 ,...Ω p The influence coefficient is obtained after the measurement of the (p is more than or equal to N+M):
writing the above in the form of a matrix:
H 1 =ZX 1 (12)
repeating the steps for the rest balance surfaces to obtain:
Compared with the prior art, the invention has the beneficial effects that:
the method has the advantages that the mode shape is known on the premise that the mode balance is required to be carried out by the rotor when the rotor is operated to be close to each order of critical rotation speed, but the method is limited by practical conditions, such as complex rotor structure, the mode shape cannot be obtained accurately, and the operation of the rotor to the critical rotation speed is dangerous and difficult to realize.
Drawings
FIG. 1 is a schematic diagram of a rotor imbalance distribution.
Fig. 2 is a front 5-order mode shape curve of the rotor.
Fig. 3 is a vibration response curve of the rotor due to unbalance.
Fig. 4 is a graph of vibration response curve measurements at initial imbalance.
Fig. 5 is a comparison of vibration response amplitudes before and after balancing.
Detailed Description
The invention is further illustrated with reference to the following examples.
The embodiment adopted by the invention is a numerical simulation case, and is specifically described as follows: a rotor of uniform cross section is shown in fig. 1, the rotor length being 1m, there being unbalanced masses at different cross section positions, the specific positions and sizes being shown in table 1. The maximum working speed of the rotor is set between the third order critical speed and the fourth order critical speed, and the contribution of the modes of the sixth order and above to the vibration response of the rotor in the whole working speed is considered to be negligible, so that the superposition result of the mode components of the first 5 orders is taken to represent the real vibration response of the rotor, the first 5 orders critical speed, the damping ratio and the mode quality of the rotor are set as shown in the table 2, and the first 5 orders mode shape curve of the rotor is shown in the figure 2. From the given parameters, an imbalance response curve of the rotor over the first 5 th order critical speed range can be calculated as shown in fig. 3.
Table 1 unbalance amount size and distribution
TABLE 2 critical speed to damping ratio of the first 5 th order of the rotor
In this embodiment, the target is to balance the modal unbalance response of the first 3 rd order, and by estimating the first 5 th order critical rotation speed and damping ratio information of the rotor through simulation or experiment, the remaining information (unbalance distribution, modal mass and modal shape are all unknown amounts) of the rotor is implemented as follows:
1) Selected U 1 、U 2 and U3 The plane is used as a balance surface, and the sensor positions are arranged at U 2 The plane is located;
2) Imbalance response Y to rotor in low speed range (Ω=200-600 RPM) 0 (Ω,) The amplitude and phase of (a) are measured to obtain an unbalanced vibration response curve as shown in figure 4;
3) Let us estimate the first 5 th order critical rotation speed and damping ratio information of the rotor through simulation or experiment (see table 2), usingCalculating Z in the measured rotation speed range omega=200-600 RPM r (Ω), r=1, 2,3,4,5, further solved:
4) In U 1 The plane is added with a known test block Q 1 (1 g.mm in magnitude, 0 ° in phase), measured at Ω=200-600 RPM, gives an imbalance responseCalculating H 1 (Ω), and solve: />
5) Take off Q 1 In U 2 The plane is added with a known test block Q 2 (magnitude 1 g.mm, phase 0 °) and likewise measured at Ω=200-600 RPM, gives an imbalance responseCalculating H 2 (Ω), solving:
6) Take off Q 2 In U 2 The plane is added with a known test block Q 3 (magnitude 1 g.mm, phase 0 °) and likewise measured at Ω=200-600 RPM, gives an imbalance responseCalculating H 3 (Ω), solving:
7) Γ, X obtained in the steps 3) to 6) 1 、X 2 and X3 Substituting the values into the table 3, the balance calibration amounts to be added to the three balance surfaces were obtained:
TABLE 3 Table 3
Calibration quantity | λ 1 | λ 2 | λ 3 |
Size (g, mm) | 0.44 | 1.69 | 0.36 |
Phase (°) | 14.961 | -89.86 | -27.28 |
8) Balancing according to the balancing calibration quantity obtained in the step 7), and finishing balancing of the rotor. As shown in fig. 5, it can be seen that the imbalance response curve (broken line in the figure) after balancing by the method according to the present invention is significantly lower in the entire operating rotation speed range than before balancing (solid line in the figure), and the vibration response due to modal imbalance in the rotation speed range can be largely eliminated.
Claims (1)
1. The modal dynamic balance method based on low-rotation-speed measurement is characterized by comprising the following specific steps of:
1) For a flexible rotor system, it is in the axial direction x s The synchronous vibration response at the point due to mass imbalance is as follows:
wherein ,mr ,ω r ,ξ r ,ψ r (x) Respectively representing the r-th order modal mass, modal frequency, modal damping and modal shape, wherein omega is the shaft rotation speed, U (x) is the unbalanced distribution of the rotor, and l is the length of the rotor; typically, as the rotational speed approaches the r-th order modal frequency, the order modal response dominates, written as:
to eliminate the unbalanced response of the first N-order modes, the test weight lambda needs to be added on N planes k K=1, 2,..n, such thatThe method comprises the following steps:
solving the unbalance calibration quantity required to be added according to the above formula:
Λ=-[Ψ] -1 u
wherein Λ= [ λ ] 1 ,λ 2 ,...λ N ] T ,[Ψ] r =[ψ r (x 1 ),ψ r (x 2 ),...,ψ r (x N )],u=[u 1 ,u 2 ,...u N ] T ;
2) Taking N balance surfaces to balance the front N-order modes, and setting a measurement point to be positioned at x s Where each order mode shape is ψ i (x s ) Each row [ ψ ] of the vibration mode matrix] i Multiplying byThen write as:
when the rotor runs at any rotation speed, the rotor is at a measuring point x s The synchronous vibration response of (2) is:
wherein ,when the N-th order mode is considered before balancing, the mode after the n+m-th order is considered to contribute negligible to the N-th order before, at a lower rotation speed Ω (Ω < ω) 1 <ω N <<ω N+M+1 ) Then, the mode is truncated to the former N+M order mode to obtain:
when the different rotational speeds Ω= { Ω are performed in the first-order critical rotational speed range below the rotor 1 ,Ω 2 ,...Ω p Measurement of } (p.gtoreq.N+M), then the above formula is written as:
recording deviceZ=[Z 1 (Ω) Z 2 (Ω)...Z N+M (Ω)]When the critical rotation speed omega of each order of the rotor is known i (i=1, 2,) n+m and critical damping ratio ζ i (i=1, 2,..n+m), then solved by the formula:
in the formula ,(·)+ Representing the pseudo-inverse of the matrix, y= [ Y (x s ,Ω 1 ),Y(x s ,Ω 2 ),...,Y(x s ,Ω p )] T ;
Adding a known test weight Q to the first equilibrium surface 1 The same applies to different rotational speeds Ω= { Ω 1 ,Ω 2 ,...Ω p The influence coefficient is obtained after the measurement of the (p is more than or equal to N+M):
writing the above in the form of a matrix:
H 1 =ZX 1
repeating the steps for the rest balance surfaces to obtain:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210411949.5A CN115014637B (en) | 2022-04-19 | 2022-04-19 | Modal dynamic balance method based on low-rotation-speed measurement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210411949.5A CN115014637B (en) | 2022-04-19 | 2022-04-19 | Modal dynamic balance method based on low-rotation-speed measurement |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115014637A CN115014637A (en) | 2022-09-06 |
CN115014637B true CN115014637B (en) | 2023-05-26 |
Family
ID=83067770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210411949.5A Active CN115014637B (en) | 2022-04-19 | 2022-04-19 | Modal dynamic balance method based on low-rotation-speed measurement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115014637B (en) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110702314A (en) * | 2019-10-10 | 2020-01-17 | 西北工业大学 | Non-trial-weight transient high-speed dynamic balancing method based on flexible rotor acceleration response information |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5214585A (en) * | 1989-06-30 | 1993-05-25 | General Electric Company | Balancing method and product |
US5544073A (en) * | 1994-06-02 | 1996-08-06 | Computational Systems, Inc. | Rotor balancing calculator |
CN1204385C (en) * | 2003-07-10 | 2005-06-01 | 西安交通大学 | Flexible rotor low speed holographic dynamic balancing method |
RU2453818C1 (en) * | 2011-01-11 | 2012-06-20 | Федеральное государственное унитарное предприятие "Российский федеральный ядерный центр-Всероссийский научно-исследовательский институт технической физики имени академика Е.И. Забабахина" | Method of adjusting balancing stand for determining rotor mass-inertia asymmetry parameters |
CN102507120B (en) * | 2011-11-18 | 2014-02-26 | 上海交通大学 | Generator rotor unbalance magnetic tension testing bench supported by sliding bearings |
RU2499985C1 (en) * | 2012-04-11 | 2013-11-27 | Федеральное государственное унитарное предприятие "Российский федеральный ядерный центр - Всероссийский научно-исследовательский институт технической физики имени академика Е.И. Забабахина" | Rotor balancing method in one correction plane |
US9181804B1 (en) * | 2015-02-17 | 2015-11-10 | Borgwarner Inc. | Ball bearing turbocharger balancer |
EP3179611B1 (en) * | 2015-12-10 | 2018-06-27 | Skf Magnetic Mechatronics | Balancing method for balancing at high speed a rotor of a rotary machine |
CN112729681A (en) * | 2020-12-22 | 2021-04-30 | 厦门大学 | Detection method and detection system for rotor field dynamic balance |
-
2022
- 2022-04-19 CN CN202210411949.5A patent/CN115014637B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110702314A (en) * | 2019-10-10 | 2020-01-17 | 西北工业大学 | Non-trial-weight transient high-speed dynamic balancing method based on flexible rotor acceleration response information |
Non-Patent Citations (1)
Title |
---|
李志炜等.双转子耦合系统低速动平衡实验研究.《化学工程与装备》.2019,(第第6期期),第1-5页. * |
Also Published As
Publication number | Publication date |
---|---|
CN115014637A (en) | 2022-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5214585A (en) | Balancing method and product | |
CN104568313A (en) | Influence coefficient dynamic balance method on rotating machine with multiple plane, multiple-points and multiple revolving speed shafting | |
CN111475903A (en) | Large-scale high-speed rotation equipment multistage part dynamic characteristic step-by-step measuring, adjusting and distributing method based on multi-bias error synchronous compensation | |
CN110501107B (en) | Spacecraft rotation load dynamic balance measurement method based on six-dimensional force tester | |
CN109115408B (en) | Dynamic balance test method of large hydroelectric generating set based on centrifugal force equation | |
CN104101464B (en) | Multi-wheel-disc rotor dynamic balancing test method based on rotating coordinate system | |
US20220252478A1 (en) | Method and drivetrain test bench for detecting an imbalance and/or a misalignment | |
CN108225783B (en) | Method and device for balancing fan rotor of aviation turbofan engine | |
CN109406053B (en) | dynamic balancing method for rotor without removing weight | |
CN115014637B (en) | Modal dynamic balance method based on low-rotation-speed measurement | |
CN113190786B (en) | Vibration prediction method for large-scale rotating equipment by utilizing multidimensional assembly parameters | |
CN110702314B (en) | Non-trial-weight transient high-speed dynamic balancing method based on flexible rotor acceleration response information | |
CN104165729B (en) | A kind of dynamic balance method of high speed rotor | |
Pesek et al. | Numerical analysis of dry friction damping effect of tie-boss couplings on three blade bundle | |
CN111473859B (en) | Method for formulating vibration limiting value of complete machine | |
JP2012002615A (en) | Unbalance amount measuring method and device therefor | |
CN113945230A (en) | Identification method for high-order error coefficient of inertial device | |
Tiwari et al. | An extended field balancing procedure for flexible rotors fully levitated by active magnetic bearings | |
Pennacchi et al. | A model-based prediction of balancing behavior of rotors above the speed range in available balancing systems | |
CN110702313A (en) | Method for high-precision identification of unbalanced excitation force of flexible rotor based on variable-speed starting | |
CN104483067A (en) | Measuring method of dynamic unbalance amount of fan | |
Shorr et al. | Numerical and experimental estimation of the turbine blade damper efficiency | |
Grządziela et al. | The impact estimation of damping foundations in dynamics of the rotor system | |
CN113358282B (en) | Low-speed orthogonal fusion dynamic balancing method for composite material tail shaft on dynamic balancing machine | |
Li et al. | Effect analysis of rotor balancing under conditions of external forces |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |