# An Optimized Modal Balancing Approach for a Flexible Rotor Using a Vibration Response While the Rotor Is Speeding Up

Authors: Liqing Li; Yuanhang Hou; Shuqian Cao
Journal: Shock and Vibration
 (2022)

---

## Abstract

The modal balancing method (MBM) is an effective method for reducing the vibration caused by the unbalancing of a rotor. A rotor’s modal parameters, especially its modal shapes, need to be accurately calculated in this method. This paper proposes an optimized modal balancing approach for flexible rotors. The vibration modes of the rotor are first obtained with experimental modal analysis based on the rotor’s response signals while the rotor is speeding up. The rotor balancing strategy is subsequently optimized by the sensitivity analysis of the mode shapes. Orthogonal trial masses are obtained based on the orthogonality of each vibration mode, and the correction masses are finally calculated by using the influence coefficients of the trial masses. An experimental result is shown to demonstrate and validate that the proposed approach is able to achieve superior accuracy compared to the conventional MBM.


---

## Body

## 1. Introduction

A rotor is the core component of rotating machinery. The main vibration of rotating equipment is caused by the unbalance of a rotor. Rigid rotors were used in early rotating machinery. The working speed of these rotors was far lower than their critical speed, and hence, the rotors did not produce large bending deformations [1]. With the growing demand for power, flexible rotors have been increasingly used in rotating machineries, such as in gas turbines and aerospace engines. The working speed of these rotors far exceeds the first critical speed and is even higher than the second and third critical speeds [2]. When passing through different critical speeds, a rotor will show different forms of bending deformation. The deformation at the critical speed will change the initial unbalanced distribution of a rotor and even cause additional unbalance. Therefore, flexible rotors must be accurately balanced to safely pass through each critical speed to the designed speed [3, 4].Generally, rotor balancing methods include the influence coefficient method (ICM), modal balancing method (MBM), and comprehensive balancing method [5–10]. The ICM obtains the influence coefficient by adding trial masses to different correction planes of the rotor in trial runs [11–13]. Because each influence coefficient is measured via experiments, a significant number of revolutions are required to obtain the correction data at the balancing speed. Turbomachinery users are constantly concerned about the increased cost and possible downtime associated with the ICM due to its limitations and difficulties in flexible rotor balancing [14, 15]. In addition, the selection of balancing planes is heavily dependent on the operator’s experience. The lack of rotor dynamic characteristics may lead to balance failure [16–18]. Therefore, the application of the ICM in flexible rotor balancing is seriously restricted.The MBM is a method based on flexural deformation for an unbalanced rotor. The unbalance can be expressed in the form of each vibration mode, and the rotor deflection of each mode can be excited by the corresponding vibration unbalance. The method generally assumes that the rotor system has a plane vibration mode. Balancing one mode will not affect any other modes and the effects of higher modes are not considered. Grobel [19] proposed a balancing technique similar to MBM for turbine rotors. Kellenberger [20] performed an in-depth study of the MBM to continuously improve the theory of balancing flexible rotors with the MBM and proposed an N-plane method and an N + 2 planes method for flexible rotor balancing. Determining the mode shapes of a rotor is a key step for modal balancing, and these shapes are generally calculated with the finite element (FE) method. Bertoneri [21] developed a dynamical balancing system by using the FE model to calculate the rotor vibration mode for the correction masses in balancing. The FE model was relatively simple in the balancing system. It was only applicable to simply supported rotors. Prior to balancing a rotor, the structural parameters of the rotor need to be imported into the FE model. Saldarriaga [22] regarded the balancing of a flexible rotor as a process for solving two types of inverse problems. The first inverse problem is used to obtain dynamical characteristics (modal parameters) of a rotor system and the second inverse problem is used to determine the unbalance. The FE method is also used to obtain the rotor modal parameters, while the parameters of the pedestals are identified using the frequency response function (FRF). Khulief [23] adopted a balancing method and developed a hybrid experiment and analysis technique to balance a high-speed flexible rotor. The technology combines FE modeling, experimental modal analysis, vibration measurement, and mathematical identification. Li [24] proposed a trialless rebalancing method based on the FE model. This method uses the vibration response and the accurate FE model of a rotor system to obtain the correction masses of a rotor. These methods without trial masses eliminate trial runs in rotor balancing to save time [25]. In most MBMs, a highly accurate FE model of a rotor is required to obtain the modal information. Sinha [26] found that the maximum error in the estimated unbalance amplitude could reach up to 45% for simultaneous random errors of 5% in the FE model of a rotor. It is always difficult to obtain accurate FE models for complex rotors. This difficulty will be obvious for the operator when balancing various rotors. Additionally, it is very dangerous in most MBMs to run a rotor stably at a critical speed.In this paper, an optimized MBM is proposed for balancing flexible rotors. The procedure of the approach is shown in Figure1. Firstly, the modal parameters of a rotor, such as the mode shapes and critical speeds, are calculated by using the vibration response of the rotor while the rotor is speeding up. Secondly, the modal parameters are used to calculate the orthogonal trial masses and calculate the sensitivity. Thirdly, the mode shape sensitivity analysis for the local mass changes of the rotor is used to optimize the balancing strategy in selecting correction planes, and the correction masses are calculated by using the MBM. Finally, experimental validation is conducted to assess the achievable accuracy of the proposed method and the conventional MBM.Figure 1
The schematic illustration of the approach.


## 2. Modal Parameter Estimation for the Rotor

An isotropic rotor is considered to be a multidegree of freedom linear vibration system. The generalized complex coordinate expression of its motion equation is shown in equation (1) [27]. The gyroscopic effect of the rotor disc has a great influence on the dynamics of the rotor system, especially for large rotors. For most rotors, due to the existence of the gyroscopic effect, the Campbell diagram of the system is not as “flat” as the Jeffcott rotor.(1)Mq¨t+ΩG+Cq˙t+Kqt=ft,where M, C, and K are the mass, stiffness, and damping matrices of the rotor system, respectively, and they are independent of the rotational speed of the rotor; G is the gyroscopic matrix of the rotor, which is dependent on the rotational speed of the rotor; qt is the dimensional generalized coordinate vector; ft is the corresponding generalized force vector. The FRF matrix of the rotor system described by (1) results from the inversion of the impedance matrix. It can be expressed as follows:(2)Hω=Zs−1|s=jω=K+jΩG+Cω−Mω2−1,where Zs is the impedance matrix of the rotor system, which can be obtained by Laplace transform. The FRF for the two physical coordinates k and s can be expressed using the residue as follows:(3)Hksω=∑r=1NRksrjω−λr+Rksr∗jω−λr∗=∑r=1Narφkrφsrjω−λr+ar∗φkr∗φsr∗jω−λr∗=∑r=12Narφkrφsrjω−λr,where ∗ represents the conjugation and λr is the pole of the system. λr=−σr+jωr and λr∗=−σr−jωr, where σr is the attenuation coefficient, ωr is the modal frequency, N is the number of modes, ar is the modal scale factor for mode r, and φkr and φsr represent the components of r-th mode shape vectors of the rotor at the coordinates k and s, respectively.Generally, the initial unbalanced mass distributions of a rotor are generated during rotor manufacturing. The initial distributions are always unknown. Since the unbalance distributions can be expressed in the form of vibration modes,n modal unbalances can be obtained, and these can be written as the column vector u=u1u2⋯unT. The excitation force generated by the unbalancing of the rotor is as follows:(4)ft=uΩ2ejΩt.For the case ofω=Ω, the Fourier transform of the unbalanced force of the rotor system can be expanded as follows:(5)fω=1T∫−T/2T/2uΩ2ejΩte−jωtdt=uω2.The vibration of thek-th point on the rotor is a superposition response by S excitation forces. It can be written as follows:(6)xkω=∑s=1SHksωfsω=∑s=1SHksωusω2.The ratio of the response at testing pointk to the response at testing point p on a structure is defined as the transfer rate (TR) from testing point p to testing point k, where the testing point p is the reference point. If the modes of the rotor system can be effectively separated, it means that the coupling between each mode is small, and the response of the system will be dominated by the r-th modal vibration at the modal frequency ωr. The contribution of the other modes at this frequency can be ignored. The TR at the critical speed ωr can be written as follows:(7)xkωxpω≈urω2∑s=1SHksωurω2∑s=1SHpsω=∑s=1S∑r=12Narφkrφsr/jω−λr∑s=1S∑r=12Narφprφsr/jω−λr≈∑s=1Sarφkrφsr/jω−λr∑s=1Sarφkrφsr/jω−λr=φkrφpr.The rotational speed can be measured accurately. The vibration response at each test point on the rotor can be measured in the integer period (T) of the rotational speed. For the finite length signal xt|t∈0,T, the finite Fourier transform is xω=∫0Txte−jωtdt. With this definition, the autopower spectrum and the cross-power spectrum of the two signals xkt and xpt are defined as follows:(8)Gppωr=1Txp∗ωrxpωr,Gpkωr=1Txp∗ωrxkωr.Combining (7) and (8), we arrive at the following:(9)φkrφpr=xkωrxpωr=GpkωrGppωr.It is assumed that the rotor system is linear. In the power spectrum amplitude curve, the frequency corresponding to every peak value is the critical speed of a rotor. The power spectrum of each testing position on the rotor should contain the same peak frequency unless the testing position is the modal node of a rotor. According to (9), the ratio of the r-th mode shape component of testing point k to testing point p is equal to the ratio of the cross-power spectrum Gpkω to the autopower spectrum Gppω at the critical speed. φkr/φpr is complex, and it contains the magnitude and the phase angle. This can be regarded as the normalized mode shape with reference point p.Since the excitation forces are unmeasurable, only the response spectrum can be used to determine the modal parameters of the rotor. It is necessary to synthesize the response spectrum of all the measuring points and the structural characteristics to comprehensively determine the modal parameters.


## 3. Correction Plane Optimization with Sensitivity Analysis

Assuming(10)P=ΩGMM0,Q=K00−M,q′t=qtq˙t.The state space equation of the rotor system described by (1) can be written as follows [28]:(11)Pq˙′t+Qq′t=0.The eigenvalues and eigenvectors in the state space are set toλi, Ψi, and P and Q are assumed to be the coefficient matrices of the state equation. From the orthogonality of the eigenvectors, the following expression is obtained:(12)ΨkTPΨi=0,i≠k,1,i=k,(13)ΨkTQΨi=0,i≠k,−λi,i=k,in whichΨi has been normalized according to certain conditions.From the generalized eigenvalue problem of the state spaceλP+QΨ=0, the following expression is obtained:(14)λiPΨi+QΨi=0.Calculating the partial derivative of the local massm gives the following:(15)∂λi∂mPΨi+λi∂P∂mΨi+λiP∂Ψi∂m+∂Q∂mΨi+Q∂Ψi∂m=0.MultiplyingΨkT by the left side of (15) gives the following:(16)∂λi∂mΨkTPΨi+ΨkTλi∂P∂m+∂Q∂mΨi+λiΨkTP+ΨkTQ∂Ψi∂m=0.From (14), the following expression is obtained:(17)λiΨiTP+ΨiTQ=0.Wheni≠k, from equations (12), (16), and (17), the following expression is obtained:(18)λi−λkΨkTP∂Ψi∂m+ΨkTλi∂P∂m+∂Q∂mΨi=0.For eigenvectorsΨs constituting a set of bases of N-dimensional vector space, ∂Ψi/∂m is a linear combination of the eigenvectors:(19)∂Ψi∂m=∑s=1nasΨs.Substituting (19) into (18) and considering (12) yields(20)ak=1λk−λiΨkTλi∂P∂m+∂Q∂mΨi,k≠i.Substituting (10) into (20), the following expression is obtained:(21)ak=1λk−λiΦkTλi2∂M∂m+λi∂ΩG∂m+∂K∂mΦi=λi2λk−λiΦkTΦi,k≠i.Wheni=k, the partial derivative of (12) is calculated directly, and considering the symmetry of P, the following expression is obtained:(22)2ΨkTP∂Ψi∂m+ΨkT∂P∂mΨi=0.Substituting (19) into (22) and considering (12), the following expression is obtained:(23)ak=ai=−12ΨiT∂P∂mΨi,i=k.In the same way as (21), by substituting (10) into (2) 3, the following expression is obtained:(24)ak=ai=−12ΦiT2λi∂M∂m+∂ΩG∂mΦi=−λΦi2,k=i.The mode shape sensitivity of the rotor for mass changes is defined as follows:(25)SΦ,m=∂Φi∂m=∑k=12nakΦk.In whichak can be obtained from (21) and (24). It can be written in more detail as follows:(26)whenk=iak=−2λiφeiφfi,e≠f,−λiφei2,e=f.(27)whenk≠i:ak=λi2λk−λiφekφfi+φeiφfk,e≠f,λi2λk−λiφekφei,e=f.In (26) and (27), λi=2πfi is the i-th eigenvalue of the rotor. φ is the component of the eigenvector. i and k represent the modal numbers. e and f represent the testing points.From the foregoing analysis, the mode shape sensitivity of the rotor for the mass changes can be expressed as a linear combination of the mode shapes, whereak is the weight of the k-th mode. By analyzing the mode shape sensitivities for the local mass, the balancing strategy can be optimized by selecting the appropriate correction planes.


## 4. The Optimized Modal Balancing Approach

According to theN-plane rotor balancing theory, N correction planes are required to balance N modal unbalances [20]. The balancing condition is as follows:(28)φx11φx21⋯φxN1φx12φx22⋯φxN2⋮⋮⋮φx1Nφx2N⋯φxNNU1x1U2x1⋯UNx1U1x2U2x2⋯UNx2⋮⋮⋮U1xNU2xN⋯UNxN=−u10⋯00u2⋯0⋮⋮⋮00⋯uN,where φxNN is the component of the mode shapes of the rotor, the subscript xN is the correction plane, and N is the number of modes, as well as the number of critical speeds, to be balanced. UNxN is the component of the corrections matrix. uN and UNxN are complex and include the correction mass and the correction phase mass. uN is the N-th component of the modal unbalance, which is equivalent to the initial mass unbalance distribution of the rotor.To ensure that the correction of each mode only affects the current mode shape, the correction value vector must be orthogonal. Each column of the correction matrix is a set of mass corrections installed on correction planes. As can be seen from (28), each column of the correction matrix is a set of mass corrections installed on the correction planes that only balance the corresponding modal unbalance and do not affect the others.It is assumed that the trail massesWi=wix1wix2⋯wixNT of each mode on the correction plane satisfy modal orthogonality:(29)φx11φx21⋯φxN1φx12φx22⋯φxN2⋮⋮⋮φx1Nφx2N⋯φxNNw1x1w2x1⋯wNx1w1x2w2x2⋯wNx2⋮⋮⋮w1xNw2xN⋯wNxN=I.where all the mode shapes have been obtained with the method mentioned in Section 2. The trial masses matrix can be obtained via matrix inversion:(30)w1x1w2x1⋯wNx1w1x2w2x2⋯wNx2⋮⋮⋮w1xNw2xN⋯wNxN=φx11φx21⋯φxN1φx12φx22⋯φxN2⋮⋮⋮φx1Nφx2N⋯φxNN−1.Comparing (28) and (29), the following expression is obtained:(31)Uix1Uix2⋮UxN=−kiwix1wix2⋮wixN,Ui=−kiWi.The trial masses are proportional to the correction masses. The multipleski are obtained with trial runs that are equal to the modal unbalance ui. From (28), the vibration caused by the initial unbalances and the correction masses of the rotor cancel each other out. In other words, the original vibration must be equal to the vibration caused by the correction masses in the rotor balancing. ui is unknown, but the proportionality constant ki can be obtained with an experiment. The total process of this approach is as follows:Step 1.
Speeding up the rotor to a speed that exceeds all critical speeds to be balanced. In this approach, due to the huge amplitude at the critical speed, the rotor should pass each critical speed quickly. The vibration data of the rotor during this process must be tested. The initial vibration of the unbalanced rotor isV0ixn, where i is the modal order of the rotor and n the test point number of the rotor. Because the rotor is rotating at the critical speed, its response can be considered purely caused by a modal unbalance with the same order.Step 2.
Stopping the rotor. The trail masses corresponding to each critical speed should be installed on the rotor.Step 3.
Repeating Step1. The vibration of the rotor V1ixj can be measured. At this time, the vibration is caused by both the initial unbalance and the trail masses. Then the following relationship can be established:(32)V0iV1i−V0i=Uiwi=−ki.Step 4.
The correction masses of each plane can be calculated with equation (32).
Since the trial masses are calculated based on the modal shapes of the rotor, theoretically, the ratioki should be equivalent at any correction planes or corresponding testing points. In practice, the data for any of the vibration sensors can be chosen to calculate the correction masses unless the sensor is only arranged at the modal node. Additionally, the least square method can be used to optimize the data in the calculation.


## 5. Experimental Validation

### 5.1. The Rotor Experimental System

The rotor balancing experiment is carried out on a DHRMT rotor bench. The schematic diagram of the rotor system is shown in Figure2, and a picture of the bench is shown in Figure 3. The bench is composed of a V-shaped pedestal, a speed-regulating motor, two bearings, an elastic coupling, a rotor shaft, three rotor discs, three eddy current sensors, and a photoelectric sensor. The bench is 810 mm long, 335 mm wide, and 133 mm high. The bench can be affixed to the ground with six feet. The rotor shaft has a length of 560 mm and a diameter of 10 mm. In order to make the rotor work flexibly, two different dimensions of the rotor discs are designed and manufactured. Disc 1 and disc 3 have a diameter of 78 mm, a thickness of 25 mm, a weight of 800 g, and a correction radius of 50 mm. The diameter of disc 2 is 78 mm, the thickness is 12.5 mm, the weight is 400 g, and the correction radius is 50 mm. With the rotor discs installed, the span of the two bearings is 500 mm. The three rotor discs divide the shaft into four segments, and the length of each segment is 125 mm. The bench uses a speed-regulating motor to drive the shaft through an elastic coupling to rotate at different speeds. The maximum output power of the motor is 148 W. The motor driver rectifies the 220 V AC power supply and outputs a PWM signal to drive the motor. The motor speed range is 0–6000 RPM.Figure 2
Schematic diagram of the rotor experimental system.Figure 3
The setup of the rotor experimental system.A photoelectric sensor is installed on the side of the coupling with a bracket. A piece of a reflective paper is attached to the coupling, and the photoelectric sensor generates a pulse in a rotating circle. This pulse signal serves as a reference for the rotor speed and vibration phase. Three eddy current sensors are installed near the discs with the bracket to test the vertical vibration response of the shaft at the three discs. The technical parameters of the three eddy current sensors after calibration are shown in Table1. A DH5922 N dynamic signal collector is used to collect the vibration signal of these sensors. The maximum sampling frequency of DH5922 N is 128 kHz, and the A/D module is 24 bits.Table 1
Eddy current sensor parameters.
Parameters#1#2#3Measuring range (mm)111Sensitivity (V/mm)9.989.9810.03Amplitude linearity0.21%0.60%0.40%Amplitude repeatability0.33%0.64%0.37%Zero error0.01%0.01%0.02%Supply voltage (V)±15±15±15The sensor number corresponds to the disc number. That is, sensor #1 tests the vibration of disc 1, sensor #2 tests the vibration of disc 2, and sensor #3 tests the vibration of disc 3. The three discs are also the correction planes. Threaded holes are set in the discs to mount the correction masses during balancing. These are shown in Figure4 and 5.Figure 4
Disc with holes for mounting correction masses.Figure 5
Electronic balance and small screws for correction masses.


### 5.2. Modal Parameters of the Rotor Experimental System

The motor speed of the bench shown in Figure5 can only reach 6000 rpm, and it is difficult to reach the third critical speed of the rotor. Therefore, vibration testing and balancing are mainly carried out at the first and the second critical speeds. The rotor is dynamically balanced at a low speed to reduce the original imbalance of the rotor beforehand so that the rotor can pass the first two critical speeds.As is generally known, the vibration of a structure is very large at resonance, and this is dangerous to a user. When the stress in the vibration exceeds the elastic stress limit of the structural material, it will cause permanent damage to the structure. The same is true for rotors. In a traditional MBM, a rotor runs at a certain critical speed constantly. The vibration response of the rotor will be too large to cause damage to the rotor shaft, bearing, or other parts. In addition, this is dangerous for operators. In the process of speeding up, the rotor can pass through the critical speed quickly without staying. The calculation of the modal parameters and correction masses using the vibration signals of a rotor while the rotor is speeding up is important not only for safety but also for efficiency and precision.The vibration displacements of the three discs of the rotor before flexible balancing are shown in Figure6. It can be seen from these curves that although the rotor can pass the two first critical speeds after low-speed balancing, the vibration displacements at the two critical speeds are still relatively large for the vibration mode excited by the corresponding unbalance. The maximum amplitude is 0.18 mm.Figure 6
Vibration displacement of three discs while the rotor is speeding up.Taking the vibration displacement on the first disc as the reference, the correlation functions of the vibration displacement for the three discs are as shown in Figures7 and 8.Figure 7
Magnitude curve of the power spectrum.Figure 8
Phase angle curve of the power spectrum.Three power spectrums are calculated as shown in Figures7 and 8, and these spectrums can be used to identify the modal frequencies and the mode shapes. It can be determined that the first critical speed of the rotor is 23.44 Hz, which corresponds to 1406.4 RPM. The second critical speed is 77.15 Hz, which corresponds to 4629 RPM. From Figures 7 and 8, the mode shapes of the rotor listed in Table 2 can be obtained. These mode shapes have been normalized by the maximum. Figure 9 shows the mode shape of the rotor.Table 2
Rotor mode shapes normalized by the maximum.
1st magnitude1st angle2nd magnitude2nd anglePlane 10.850.0°1.000.0°Plane 21.00−1.3°0.09−151.1°Plane 30.94−6.4°0.86179.9°Figure 9
Mode shapes of the rotor.


### 5.3. Sensitivity Analysis

The mode shape sensitivity data for the mass changes are calculated by using the method described in Section3. The mode shapes of the rotor are shown in Figure 9. The calculation results are shown in Table 3. In the table, mi is the local mass in the i-th plane and φir is the r-th mode shape component of the i-th plane.Table 3
Sensitivity matrix for mass changes in the three planes.
Sensitivitym1m2m3φ11−0.0521−0.1342−0.0481φ21−0.1185−0.1562−0.1189φ31−0.0480−0.1342−0.0436φ12−0.9497−0.0395−1.0158φ22−0.6605−0.0430−0.6663φ32−0.9823−0.0398−1.0536Surveying the data horizontally for the first mode in Table3, the absolute values of the sensitivity φ21 for m1,m2 and m3 in the first mode are larger. This implies that changing the three local masses has an obvious effect on the second component in the first mode. In these three sensitivities, the absolute value of the sensitivity of φ21 for m2 is about 30% larger than that for m1 and m3. Surveying the data vertically in Table 4, the absolute values of the sensitivity φ11, φ21, and φ31 for m2 are much larger than those for m1 and m3. The absolute values of sensitivity of φ11 and φ31 for m1 and m3 are much smaller than the other sensitivities. This means that changing the local mass of the first and third correction planes (discs) has little effect on the first mode while changing the local mass of the second correction plane (disc) has a significant effect on the first mode. This result can be clearly shown in Figure 10.Table 4
Comparison of different methods for the rotor balancing at the first critical speed.
Conventional MBMOptimization approachCorrection planePlane1Plane 2Plane 3Plane 1Plane 2Plane 3Initial amplitude (mm)0.14080.16760.15250.14080.16760.1525Initial phase (°)−31.4−30.8−31.1−31.4−30.8−31.1Correction mass (g)0.450.370.44\1.1\Correction phase (°)150.0150.0150.0\150.0\Actual installation phase (°)157.5157.5157.5\157.5\Residual amplitude (mm)0.10110.11420.13110.08210.08720.0574Balancing rate28%32%14%42%48%63%Average balancing rate25%51%Figure 10
Histogram of the first mode sensitivity matrix.Surveying the data about the second mode horizontally in Table3, the changes of m1. m2 and m3 have a relatively large influence on the components φ12 and φ32 in the second mode. The sensitivities of φ12 and φ32 for m1 and m3 are almost the same. Taking the data in the first row as an example, the absolute values of the sensitivity φ12 for m1 and m3 are about 25 times larger than the sensitivity φ12 for m2. Surveying the data vertically in Table 5, the absolute sensitivity values of the three mode components for m2 are much smaller than those for m1 and m3. The sensitivity of φ22 for m2 has the largest absolute value in line 2, which is 0.0430. The sensitivity of φ22 for m1 has the smallest absolute value in line 1 and line 3, which is 0.6605. The latter is about 16 times larger than the former. This means that changing the local mass of the first and third correction planes has a subequal influence on the second mode shape, which is more considerable than changing the local mass of the second plane. This result can be clearly seen in Figure 11.Table 5
Comparison of different methods for the rotor balancing at the second critical speed.
Conventional MBMOptimization approachCorrection planePlane 1Plane 2Plane 3Plane 1Plane 2Plane 3Initial amplitude (mm)0.11380.04310.13630.11380.04310.1363Initial phase (°)152.1−42.3−41.5152.1−42.3−41.5Correction mass (g)0.881.060.840.7\0.8Correction phase (°)−28.0138.0138.0−28.0\138.0Actual installation phase (°)337.5135.0135.0337.5\135.0Residual amplitude (mm)0.08340.04200.09800.03410.04270.0780Balancing rate27%3%28%70%3%43%Average balancing rate19%39%Figure 11
Histogram of the second mode sensitivity matrix.


### 5.4. Experimental Results of the Rotor Balancing

After low-speed balancing, the rotor can pass two critical speeds safely. However, the vibration of the rotor is still high at critical speeds. Therefore, the rotor needs high-speed balancing at critical speeds next. The initial vibration amplitude and the phase angle of the rotor at these critical speeds are extracted from the time domain signals while the rotor is speeding up, as shown in Figure6. The phase angles of the vibrations for the three discs are calculated based on the pulse signal of the photoelectric sensor.The optimization analysis shows that balancing the vibration at the first critical speed only requires operating disc 2; balancing the vibration at the second critical speed requires operating disc 1 and disc 3 concurrently. Since the first mode shape sensitivities of the mass changes on disc 1 and disc 3 are relatively low, the second mode shape has an extremely low sensitivity for the mass changes on disc 2. That is, disc 1 and disc 3 have less effect than disc 2 in the first mode. Disc 2 has less effect than disc 1 and disc 3 on the second mode. According to this strategy, the optimized balancing method mentioned above is used to rebalance the rotor.The data for balancing by using the conventional MBM and the optimization approach are shown in Tables4 and 5. Table 4 shows the comparison of the two methods at the first critical speed of the rotor. Table 5 shows the second critical speed.After balancing the flexible rotor with the conventional MBM, the vibration amplitude of each plane has been reduced. Unfortunately, the reduction of the vibration amplitude in the three discs is not obvious. The highest balancing rate of the first critical speed is only 32%. The average balancing rate in the first critical speed is only 25%. After the correction planes have been optimized, the balancing rate has improved. The maximum balancing rate at the first critical speed is 63%. The average balancing rate is 51%, which is twice as much as that before optimization. From Table5, it can also be seen that the optimization approach at the second critical speed of the rotor also has a higher balance accuracy than the traditional method.The vibration response plots of the three discs while the rotors are speeding up are presented in Figure12. In this figure, the balancing results before balancing, for the conventional MBM and for the optimization approach for the rotor are shown. It is worth mentioning that the vibration response of disc 2 at the second critical speed changes slightly before and after balancing. This is because the second disk is at the node of the second vibration mode of the rotor. It is obvious from the figure that the optimization approach has a better balance effect.Figure 12
Vibration response of three discs during speeding up.


## 5.1. The Rotor Experimental System

The rotor balancing experiment is carried out on a DHRMT rotor bench. The schematic diagram of the rotor system is shown in Figure2, and a picture of the bench is shown in Figure 3. The bench is composed of a V-shaped pedestal, a speed-regulating motor, two bearings, an elastic coupling, a rotor shaft, three rotor discs, three eddy current sensors, and a photoelectric sensor. The bench is 810 mm long, 335 mm wide, and 133 mm high. The bench can be affixed to the ground with six feet. The rotor shaft has a length of 560 mm and a diameter of 10 mm. In order to make the rotor work flexibly, two different dimensions of the rotor discs are designed and manufactured. Disc 1 and disc 3 have a diameter of 78 mm, a thickness of 25 mm, a weight of 800 g, and a correction radius of 50 mm. The diameter of disc 2 is 78 mm, the thickness is 12.5 mm, the weight is 400 g, and the correction radius is 50 mm. With the rotor discs installed, the span of the two bearings is 500 mm. The three rotor discs divide the shaft into four segments, and the length of each segment is 125 mm. The bench uses a speed-regulating motor to drive the shaft through an elastic coupling to rotate at different speeds. The maximum output power of the motor is 148 W. The motor driver rectifies the 220 V AC power supply and outputs a PWM signal to drive the motor. The motor speed range is 0–6000 RPM.Figure 2
Schematic diagram of the rotor experimental system.Figure 3
The setup of the rotor experimental system.A photoelectric sensor is installed on the side of the coupling with a bracket. A piece of a reflective paper is attached to the coupling, and the photoelectric sensor generates a pulse in a rotating circle. This pulse signal serves as a reference for the rotor speed and vibration phase. Three eddy current sensors are installed near the discs with the bracket to test the vertical vibration response of the shaft at the three discs. The technical parameters of the three eddy current sensors after calibration are shown in Table1. A DH5922 N dynamic signal collector is used to collect the vibration signal of these sensors. The maximum sampling frequency of DH5922 N is 128 kHz, and the A/D module is 24 bits.Table 1
Eddy current sensor parameters.
Parameters#1#2#3Measuring range (mm)111Sensitivity (V/mm)9.989.9810.03Amplitude linearity0.21%0.60%0.40%Amplitude repeatability0.33%0.64%0.37%Zero error0.01%0.01%0.02%Supply voltage (V)±15±15±15The sensor number corresponds to the disc number. That is, sensor #1 tests the vibration of disc 1, sensor #2 tests the vibration of disc 2, and sensor #3 tests the vibration of disc 3. The three discs are also the correction planes. Threaded holes are set in the discs to mount the correction masses during balancing. These are shown in Figure4 and 5.Figure 4
Disc with holes for mounting correction masses.Figure 5
Electronic balance and small screws for correction masses.


## 5.2. Modal Parameters of the Rotor Experimental System

The motor speed of the bench shown in Figure5 can only reach 6000 rpm, and it is difficult to reach the third critical speed of the rotor. Therefore, vibration testing and balancing are mainly carried out at the first and the second critical speeds. The rotor is dynamically balanced at a low speed to reduce the original imbalance of the rotor beforehand so that the rotor can pass the first two critical speeds.As is generally known, the vibration of a structure is very large at resonance, and this is dangerous to a user. When the stress in the vibration exceeds the elastic stress limit of the structural material, it will cause permanent damage to the structure. The same is true for rotors. In a traditional MBM, a rotor runs at a certain critical speed constantly. The vibration response of the rotor will be too large to cause damage to the rotor shaft, bearing, or other parts. In addition, this is dangerous for operators. In the process of speeding up, the rotor can pass through the critical speed quickly without staying. The calculation of the modal parameters and correction masses using the vibration signals of a rotor while the rotor is speeding up is important not only for safety but also for efficiency and precision.The vibration displacements of the three discs of the rotor before flexible balancing are shown in Figure6. It can be seen from these curves that although the rotor can pass the two first critical speeds after low-speed balancing, the vibration displacements at the two critical speeds are still relatively large for the vibration mode excited by the corresponding unbalance. The maximum amplitude is 0.18 mm.Figure 6
Vibration displacement of three discs while the rotor is speeding up.Taking the vibration displacement on the first disc as the reference, the correlation functions of the vibration displacement for the three discs are as shown in Figures7 and 8.Figure 7
Magnitude curve of the power spectrum.Figure 8
Phase angle curve of the power spectrum.Three power spectrums are calculated as shown in Figures7 and 8, and these spectrums can be used to identify the modal frequencies and the mode shapes. It can be determined that the first critical speed of the rotor is 23.44 Hz, which corresponds to 1406.4 RPM. The second critical speed is 77.15 Hz, which corresponds to 4629 RPM. From Figures 7 and 8, the mode shapes of the rotor listed in Table 2 can be obtained. These mode shapes have been normalized by the maximum. Figure 9 shows the mode shape of the rotor.Table 2
Rotor mode shapes normalized by the maximum.
1st magnitude1st angle2nd magnitude2nd anglePlane 10.850.0°1.000.0°Plane 21.00−1.3°0.09−151.1°Plane 30.94−6.4°0.86179.9°Figure 9
Mode shapes of the rotor.


## 5.3. Sensitivity Analysis

The mode shape sensitivity data for the mass changes are calculated by using the method described in Section3. The mode shapes of the rotor are shown in Figure 9. The calculation results are shown in Table 3. In the table, mi is the local mass in the i-th plane and φir is the r-th mode shape component of the i-th plane.Table 3
Sensitivity matrix for mass changes in the three planes.
Sensitivitym1m2m3φ11−0.0521−0.1342−0.0481φ21−0.1185−0.1562−0.1189φ31−0.0480−0.1342−0.0436φ12−0.9497−0.0395−1.0158φ22−0.6605−0.0430−0.6663φ32−0.9823−0.0398−1.0536Surveying the data horizontally for the first mode in Table3, the absolute values of the sensitivity φ21 for m1,m2 and m3 in the first mode are larger. This implies that changing the three local masses has an obvious effect on the second component in the first mode. In these three sensitivities, the absolute value of the sensitivity of φ21 for m2 is about 30% larger than that for m1 and m3. Surveying the data vertically in Table 4, the absolute values of the sensitivity φ11, φ21, and φ31 for m2 are much larger than those for m1 and m3. The absolute values of sensitivity of φ11 and φ31 for m1 and m3 are much smaller than the other sensitivities. This means that changing the local mass of the first and third correction planes (discs) has little effect on the first mode while changing the local mass of the second correction plane (disc) has a significant effect on the first mode. This result can be clearly shown in Figure 10.Table 4
Comparison of different methods for the rotor balancing at the first critical speed.
Conventional MBMOptimization approachCorrection planePlane1Plane 2Plane 3Plane 1Plane 2Plane 3Initial amplitude (mm)0.14080.16760.15250.14080.16760.1525Initial phase (°)−31.4−30.8−31.1−31.4−30.8−31.1Correction mass (g)0.450.370.44\1.1\Correction phase (°)150.0150.0150.0\150.0\Actual installation phase (°)157.5157.5157.5\157.5\Residual amplitude (mm)0.10110.11420.13110.08210.08720.0574Balancing rate28%32%14%42%48%63%Average balancing rate25%51%Figure 10
Histogram of the first mode sensitivity matrix.Surveying the data about the second mode horizontally in Table3, the changes of m1. m2 and m3 have a relatively large influence on the components φ12 and φ32 in the second mode. The sensitivities of φ12 and φ32 for m1 and m3 are almost the same. Taking the data in the first row as an example, the absolute values of the sensitivity φ12 for m1 and m3 are about 25 times larger than the sensitivity φ12 for m2. Surveying the data vertically in Table 5, the absolute sensitivity values of the three mode components for m2 are much smaller than those for m1 and m3. The sensitivity of φ22 for m2 has the largest absolute value in line 2, which is 0.0430. The sensitivity of φ22 for m1 has the smallest absolute value in line 1 and line 3, which is 0.6605. The latter is about 16 times larger than the former. This means that changing the local mass of the first and third correction planes has a subequal influence on the second mode shape, which is more considerable than changing the local mass of the second plane. This result can be clearly seen in Figure 11.Table 5
Comparison of different methods for the rotor balancing at the second critical speed.
Conventional MBMOptimization approachCorrection planePlane 1Plane 2Plane 3Plane 1Plane 2Plane 3Initial amplitude (mm)0.11380.04310.13630.11380.04310.1363Initial phase (°)152.1−42.3−41.5152.1−42.3−41.5Correction mass (g)0.881.060.840.7\0.8Correction phase (°)−28.0138.0138.0−28.0\138.0Actual installation phase (°)337.5135.0135.0337.5\135.0Residual amplitude (mm)0.08340.04200.09800.03410.04270.0780Balancing rate27%3%28%70%3%43%Average balancing rate19%39%Figure 11
Histogram of the second mode sensitivity matrix.


## 5.4. Experimental Results of the Rotor Balancing

After low-speed balancing, the rotor can pass two critical speeds safely. However, the vibration of the rotor is still high at critical speeds. Therefore, the rotor needs high-speed balancing at critical speeds next. The initial vibration amplitude and the phase angle of the rotor at these critical speeds are extracted from the time domain signals while the rotor is speeding up, as shown in Figure6. The phase angles of the vibrations for the three discs are calculated based on the pulse signal of the photoelectric sensor.The optimization analysis shows that balancing the vibration at the first critical speed only requires operating disc 2; balancing the vibration at the second critical speed requires operating disc 1 and disc 3 concurrently. Since the first mode shape sensitivities of the mass changes on disc 1 and disc 3 are relatively low, the second mode shape has an extremely low sensitivity for the mass changes on disc 2. That is, disc 1 and disc 3 have less effect than disc 2 in the first mode. Disc 2 has less effect than disc 1 and disc 3 on the second mode. According to this strategy, the optimized balancing method mentioned above is used to rebalance the rotor.The data for balancing by using the conventional MBM and the optimization approach are shown in Tables4 and 5. Table 4 shows the comparison of the two methods at the first critical speed of the rotor. Table 5 shows the second critical speed.After balancing the flexible rotor with the conventional MBM, the vibration amplitude of each plane has been reduced. Unfortunately, the reduction of the vibration amplitude in the three discs is not obvious. The highest balancing rate of the first critical speed is only 32%. The average balancing rate in the first critical speed is only 25%. After the correction planes have been optimized, the balancing rate has improved. The maximum balancing rate at the first critical speed is 63%. The average balancing rate is 51%, which is twice as much as that before optimization. From Table5, it can also be seen that the optimization approach at the second critical speed of the rotor also has a higher balance accuracy than the traditional method.The vibration response plots of the three discs while the rotors are speeding up are presented in Figure12. In this figure, the balancing results before balancing, for the conventional MBM and for the optimization approach for the rotor are shown. It is worth mentioning that the vibration response of disc 2 at the second critical speed changes slightly before and after balancing. This is because the second disk is at the node of the second vibration mode of the rotor. It is obvious from the figure that the optimization approach has a better balance effect.Figure 12
Vibration response of three discs during speeding up.


## 6. Conclusions

This paper proposes an optimized balancing approach to balance a flexible rotor at critical speeds. The critical speeds and mode shapes of the rotor are obtained with the experimental modal analysis method using the vibration response while the rotor is speeding up. Furthermore, the optimized selection of the rotor correction planes is obtained in the rotor balancing strategy by means of the vibration mode shape sensitivity analysis for the mass changes. The mode shapes of the rotor are used to calculate the orthogonal trial masses in balancing. Then the influence coefficients of the orthogonal trial weights in each mode are calculated with experimentation. The correction masses are obtained to balance the rotor. In every operation of the rotor, it is safe for the operator to use the vibration response signal of the rotor while the rotor is speeding up. Compared to the conventional MBM, the balancing approach proposed in this work is able to achieve better accuracy. This approach is expected to be used in the accurate balancing variety of flexible rotors, especially for rotor balancing in the field.



---

*Source: 5261279-2022-03-07.xml*