Next Article in Journal
Design Optimization Method of Feet-Lock Steel Pipe for Soft-Rock Tunnel Based on Load-Deformation Coordination
Previous Article in Journal
Numerical Simulation Study of Brittle Rock Materials from Micro to Macro Scales Using Digital Image Processing and Parallel Computing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Continuous Rotor Dynamics of Multi-Disc and Multi-Span Rotors: A Theoretical and Numerical Investigation of the Identification of Rotor Unbalance from Unbalance Responses

1
School of Mechanical Electronic & Information Engineering, China University of Mining & Technology, Beijing 100083, China
2
School of Mechanics and Civil Engineering, China University of Mining & Technology, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(8), 3865; https://doi.org/10.3390/app12083865
Submission received: 17 March 2022 / Revised: 7 April 2022 / Accepted: 8 April 2022 / Published: 11 April 2022
(This article belongs to the Topic Dynamical Systems: Theory and Applications)

Abstract

:
Rotor unbalance identification plays a critical role in balancing rotors. In this paper, concerned with multi-disc and multi-span rotor-bearing systems, two novel algorithms called the Single Direction Algorithm (SDA) and the Two Orthogonal Direction Algorithm (TODA) are proposed for identifying rotor unbalance from unbalance responses. A matrix method is proposed to solve the problem of the equations being non-linear transcendental, there being too many unknown variables in the equations, and rotor unbalances and bearing coefficients being coupled together. The unbalance responses at all the eccentric discs are necessary for identifying their unbalances. Numerical simulations are conducted to validate the proposed methods. Moreover, an adjustment point is found, and a proper sensor resolution is suggested to achieve high identification accuracy by means of numerical studies. In addition, the identification accuracy of SDA is better than TODA, and SDA is more practical and suitable for medium-speed and high-speed rotors. The proposed algorithms have the flexibility to incorporate any number of bearings and discs and provide a technique for monitoring rotor unbalance without test runs or external exciters.

1. Introduction

1.1. Background and Formulation of the Problem

Rotating machines, such as steam turbine generators, turbine compressor units, pumps, etc., play a critical role in industry and are widely used. The rotor-bearing system (rotor) is the main element in rotating machines. Rotor unbalance is a typical fault of rotors. It is inevitable due to errors in manufacturing, installation, and operation. Vibrations and even sudden breakdowns can be caused because of rotor unbalance. Identifying rotor unbalance (amplitude and phase angle) is necessary for balancing rotors to reduce vibrations and ensure the safety of rotating machines. In view of these issues, rotor unbalance estimation, which is one inverse problem of rotor dynamics, has been an active area of research.
Although rotor unbalance can be estimated so that a rotor can be balanced off-site using a dynamic balancing machine in the laboratory, complex processes, such as uncovering, dismantling, and transporting, are required. The whole process is time-consuming and expensive and is difficult to implement, especially for large rotors. Hence, monitoring rotor unbalance on-line, by which means balancing time and cost can be reduced, is very important.

1.2. Literature Survey

Rotor unbalance identification is an old problem. The modal method proposed in 1959 [1,2] and the influence coefficient method proposed in 1964 [3] are the two classical techniques used to solve the problem. Although they have been developed by many other researchers [4,5,6,7,8,9,10,11,12,13], test runs are still demanding. For large field rotors, starting and stopping the machine several times is costly and can reduce its service life [14].
Identifying rotor unbalance without test runs is the current trend [15]. This approach can be divided into two types: a method with external excitation and a method without external excitation. For the external excitation method, Bently and Muszynska [16] used frequency excitation to estimate rotor unbalance as well as bearing coefficients. Hiroshi Iida [17] applied impulse excitation on a double-disc and single-span rotor-bearing system to identify rotor unbalance as well as the stiffness and damping coefficients of two bearings. Lou [18,19] estimated rotor unbalance and bearing coefficients using an active magnetic executor to generate external excitation. However, high-power exciters are necessary for large rotors and the excitation may damage the rotors.
Consequently, methods that do not require external excitation have been developed. Aiming at a single-span rotor-bearing system, Shrivastava and Mohanty [20] identified single-disc unbalance parameters from unbalance responses using an unbalance force estimation technique. A rotor unbalance estimation method using the joint-input state estimation technique [21] and least-squares technique was proposed [22] to solve the problem that the required response measurements at different locations on the shaft may not always be accessible. Yao [23] identified single-disc and double-disc unbalance based on modal expansion combined with optimisation algorithms. The errors caused by modal expansion were reduced by an integrated modal expansion inverse problem methodology combined with an optimisation procedure. Zou [24] developed a double-disc unbalance identification method using the finite element model combined with augmented Kalman filter algorithms. While based on the continuous dynamic model and analytical solution methods [25,26], Wang [27] provided an analytical model-based algorithm to identify single-disc unbalance. Moreover, some scholars have even tried to estimate rotor unbalance and bearing coefficients simultaneously without external excitation. Tiwari [28,29] formulated an estimation algorithm using unbalance responses from three different unbalance configurations for both clockwise and anti-clockwise rotor rotations. Wang [30] proposed a simultaneous estimation of the rotor unbalance and bearing coefficients of a continuous single-disc and single-span rotor-bearing system using the Rayleigh beam model. Jamadar [31] developed a mathematical model of an unbalanced rotor using dimensional analysis and a rigid rotor approach. The factorial regression analysis method is used to solve the model. Based on it, a numerical technique for the detection of unbalance magnitude has been proposed. Ambur [32] presented an estimation method for unbalance magnitude and phase from the vibration in frequency domain. Sanches [33] proposed an identification method of unbalance for a rotor with residual shaft bent based on the finite element method and correlation analysis. Zhang [34] proposed an unbalance identification method for a high-speed rotor without trial weights based on modal analysis and the modal equivalent principle. The shortcomings of the above-mentioned methods are that they cannot incorporate any number of bearings and discs.
Therefore, Bin [14] proposed an approach based on the multi-plane influence coefficient balancing method for muti-disc and muti-span rotors. Tiwari [35,36] developed an algorithm to simultaneously estimate the rotor unbalances, four stiffness coefficients, and four damping coefficients of bearings from impulse responses and run-down responses. However, test runs or external exciters are required in these studies.

1.3. Scope and Contribution of This Study

Focusing on multi-span and multi-disc rotors, two novel algorithms, called SDA and TODA, are proposed based on the continuous rotor dynamic analysis method (CRDAM) in this paper to realize rotor unbalance identification without excitation or test runs. These algorithms have the flexibility to incorporate any number of bearings and discs. The equations of the inverse problem, which are developed based on CRDAM using unbalance responses as inputs, cannot be directly solved to obtain rotor unbalances as the equations are non-linear transcendental and there are too many unknown variables. Moreover, rotor unbalances and bearing coefficients are coupled together. A matrix method is proposed to solve the problem and identify rotor unbalance from unbalance responses. Four kinds of numerical simulations considering sensor resolutions and measurement errors are conducted for the validation. The adjustment point, which can greatly improve identification accuracy, is found by means of simulations. Proper sensor resolutions are recommended for engineering applications.
Compared with the methods described in references [14,35,36], SDA and TODA do not require test runs or external exciters. Compared with other existing methods in the literature, SDA and TODA can be applied to rotors with any number of bearings and any number of discs. Moreover, SDA and TODA can not only be used for rotors with rolling bearings but also for rotors supported by oil journal bearings. The proposed methods provide a technique for on-line monitoring of the rotor unbalances of multi-bearing and multi-disc rotors without using test runs or external exciters.

1.4. Organization of the Paper

Section 1 describes the background, the formulation of the problem, the literature, and the scope and contribution of this study. Section 2 discusses the derivation process of the proposed algorithms. Section 3 describes the numerical investigations for examining the algorithms and presents the discussion of the simulation results. Section 4 summarizes the conclusions of the study.

2. Theory

2.1. Revisting the CRDAM

There are m discs and n bearings in the rotor shown in Figure 1. The unbalance response of any position on the rotor shaft can be expressed as a function of position, rotor unbalance, and bearing stiffness and damping coefficients, according to the CRDAM. Equations (1) and (2) represent the dimensionless form of the unbalance response in the frequency domain.
U ( q ) = { j = 1 m [ π m j u w 2 e j ( sin α j j cos α j ) + w 2 m j d L U j d ] G u ( q , q j d ) j j = 1 n L ( k j j y x V j j b + k j j y y U j j b + i w c j j y x V j j b + i w c j j y y U j j b ) G u ( q , q j j b ) } L 2 E I
V ( q ) = { j = 1 m [ π m j u w 2 e j ( cos α j + i sin α j ) + w 2 m j d L V j d ] G v ( q , q j d ) j j = 1 n L ( k j j x x V j j b + k j j x y U j j b + i w c j j x x V j j b + i w c j j x y U j j b ) G v ( q , q j j b ) } L 2 E I
where U ( q ) and V ( q ) are the dimensionless unbalance response in the frequency domain along the y-axis and x-axis, respectively; m j u are the eccentric masses of # j disc; e j is the eccentric distance of # j disc; α j are the eccentric angles, which are defined as the angles between the x-axis and the eccentric position in the direction of rotation; m j d are the masses of # j disc; w is the rotation frequency; L is the length of the shaft; E is the elastic modulus of the shaft; I is the diametric shaft cross-sectional geometric moment of inertia; k j j x x , k j j x y , k j j y x , and k j j y y are the stiffness coefficients of # j j bearings respectively; c j j x x , c j j x y , c j j y x , and c j j y y are the damping coefficients of # j j bearing, respectively; z is the axial position of the shaft; z j d are the z coordinate positions of each disc; z j j b is the z coordinate position of each bearing; q , q j d , and q j j b are their dimensionless values; q = z / L , q j d = z j d / L , q j j b = z j j b / L ; G u ( q , q j d ) , G u ( q , q j j b ) , G v ( q , q j d ) , and G v ( q , q j j b ) are Green’s coefficients, which can be calculated using Green’s functions G u ( q , q i ) and G v ( q , q i ) ; and Green’s functions G u ( q , q i ) and G v ( q , q i ) can be found in reference [30]. Further, U j d represents the dimensionless unbalance response of each disc in the frequency domain in the y direction; V j d represents the dimensionless unbalance response of each disc in the frequency domain in the x direction; U j j b represents the dimensionless unbalance response of each bearing in the frequency domain in the y direction; V j j b represents the dimensionless unbalance response of each bearing in the frequency domain in the x direction.

2.2. Single Direction Algorithm

Identification of rotor unbalance amplitude m j u e j and angle α j is the inverse problem of the CRDAM. According to Equations (1) and (2), the relationships between unbalance response and rotor unbalance and the bearing coefficients are non-linear transcendental, although the eight stiffness and damping coefficients of the bearing are linearized. As the number of bearings and discs are unknown, the numbers of equations and unknown variables are also unknown. Moreover, rotor unbalance and bearing coefficients are coupled together according to Equations (1) and (2). Although enough equations can be obtained by using the unbalance responses as inputs and considering rotor unbalances and bearings coefficients as unknown variables, they cannot be solved directly to obtain the unknown variables (amplitude and angle of rotor unbalance) due to the non-linear transcendental equations, the unfixed number of equations, and the unknown variables. Hence, it is proposed that Equation (1) is written in matrix form as follows:
U ( q ) E I L 2 = [ G u ( q , q 1 d ) G u ( q , q m d ) G u ( q , q 1 b ) G u ( q , q n b ) ] T [ π w 2 m 1 u e 1 ( sin α 1 i cos α 1 ) + L w 2 m 1 d U 1 d π w 2 m m u e m ( sin α m i cos α m ) + L w 2 m m d U m d L k 1 s y x V 1 b L k 1 s y y U 1 b L k n s y x V n b L k n s y y U n b ] m + n
m + n equations can be obtained based on Equation (3) by using m + n unbalance responses as inputs. There are two unknown variables (amplitude and angle) for a disc’s unbalance. Moreover, the eight coefficients of a bearing are also unknown. Hence, there are 2 m + 8 n unknown variables for a rotor with m discs and n bearings, while only m + n equations exist. The greatest difficulty is due to the fact that the relationship among these variables is non-linear transcendental. A direct solution is not feasible. Therefore, by writing the m + n equations in matrix form, Equation (4) can be obtained as the following:
E I L 2 [ U ( q 1 s ) U ( q m s ) U ( q ( m + 1 ) s ) U ( q ( m + n ) s ) ] = H 1 [ π w 2 m 1 u e 1 ( sin α 1 i cos α 1 ) + L w 2 m 1 d U 1 d π w 2 m m u e m ( sin α m i cos α m ) + L w 2 m m d U m d L k 1 s y x V 1 b L k 1 s y y U 1 b L k n s y x V n b L k n s y y U n b ] m + n
where q 1 s , …, q m s and q ( m + 1 ) s , …, q ( m + n ) s are the dimensionless values of locations on the shaft excluding the locations of all discs and bearings; U ( q 1 s ) , …, U ( q m s ) and U ( q ( m + 1 ) s ) , …, U ( q ( m + n ) s ) are the measured dimensionless unbalance responses in the frequency domain in the y direction of locations, excluding all discs and bearings; the number of U ( q 1 s ) , …, U ( q m s ) is m and the number of U ( q ( m + 1 ) s ) , …, U ( q ( m + n ) s ) is n ; H 1 is a ( m + n ) × ( m + n ) matrix and can be calculated using Green’s functions G u ( q , q i ) as long as the locations of the measuring points, the bearings, and the discs are known; H 1 is as follows:
H 1 = [ G u ( q 1 s , q 1 d ) G u ( q 1 s , q m d ) G u ( q 1 s , q 1 b ) G u ( q 1 s , q n b ) G u ( q m s , q 1 d ) G u ( q m s , q m d ) G u ( q m s , q 1 b ) G u ( q m s , q n b ) G u ( q ( m + 1 ) s , q 1 d ) G u ( q ( m + 1 ) s , q m d ) G u ( q ( m + 1 ) s , q 1 b ) G u ( q ( m + 1 ) s , q n b ) G u ( q ( m + n ) s , q 1 d ) G u ( q ( m + n ) s , q m d ) G u ( q ( m + n ) s , q 1 b ) G u ( q ( m + n ) s , q n b ) ] ( m + n ) × ( m + n ) .
Equation (5) can be obtained according to Equation (4).
[ π w 2 m 1 u e 1 ( sin α 1 i cos α 1 ) + L w 2 m 1 d U 1 d π w 2 m m u e m ( sin α m i cos α m ) + L w 2 m m d U m d L k 1 s y x V 1 b L k 1 s y y U 1 b L k n s y x V n b L k n s y y U n b ] m + n = E I L 2 H 1 1 [ U ( q 1 s ) U ( q m s ) U ( q ( m + 1 ) s ) U ( q ( m + n ) s ) ]
The right side of Equation (5) is known. Define H 2 as following.
H 2 = E I L 2 H 1 1 [ U ( q 1 s ) U ( q m s ) U ( q ( m + 1 ) s ) U ( q ( m + n ) s ) ]
Equation (6) is obtained according to Equation (5).
{ π w 2 m 1 u e 1 ( sin α 1 i cos α 1 ) + L w 2 m 1 d U 1 d = H 2 ( 1 , 1 ) π w 2 m m u e m ( sin α m i cos α m ) + L w 2 m m d U m d = H 2 ( m , 1 )
According to Equations (6) and (7) is obtained.
{ m 1 u e 1 ( sin α 1 i cos α 1 ) = [ H 2 ( 1 , 1 ) L w 2 m 1 d U 1 d ] π w 2 m m u e m ( sin α m i cos α m ) = [ H 2 ( m , 1 ) L w 2 m m d U m d ] π w 2
According to Equations (7)–(9) are obtained.
{ m 1 u e 1 = | [ H 2 ( 1 , 1 ) L w 2 m 1 d U 1 d ] π w 2 | m m u e m = | [ H 2 ( m , 1 ) L w 2 m m d U m d ] π w 2 |
{ sin α 1 i cos α 1 = [ H 2 ( 1 , 1 ) L w 2 m 1 d U 1 d ] | H 2 ( 1 , 1 ) L w 2 m 1 d U 1 d | sin α m i cos α m = [ H 2 ( m , 1 ) L w 2 m m d U m d ] | H 2 ( m , 1 ) L w 2 m m d U m d |
Equation (8) is the amplitude of rotor unbalance for each eccentric disc. Moreover, according to Equations (9) and (10) can be obtained, by means of which the angle of rotor unbalance for each eccentric disc can be calculated.
{ α 1 = angle ( imag ( [ H 2 ( 1 , 1 ) L w 2 m 1 d U 1 d ] π w 2 ) + i real ( [ H 2 ( 1 , 1 ) L w 2 m 1 d U 1 d ] π w 2 ) ) α m = angle ( imag ( [ H 2 ( m , 1 ) L w 2 m m d U m d ] π w 2 ) + i real ( [ H 2 ( m , 1 ) L w 2 m m d U m d ] π w 2 ) )
where angle ( ) is the function used to obtain the angle of a complex number, real ( ) is the real part of a complex number, and imag ( ) is the imaginary part of a complex number.
Hence, the amplitude and angle of each eccentric disc’s unbalance can be calculated based on Equations (8) and (10), in which the unbalance responses in the y direction are required as inputs. According to Equations (4) and (5), the matrix H 2 in Equations (8) and (10) can be calculated using m + n unbalance responses as inputs. The unbalance response of the eccentric disc whose unbalance is to be identified must be measured according to Equation (7). Hence, the total number of input unbalance responses is equal to m + n + 1 . To identify all the eccentric discs’ unbalances simultaneously, the unbalance responses of all the discs should be included.
Similarly, Equation (11) can be obtained according to Equation (2) in the x direction.
V ( q ) E I L 2 = [ G v ( q , q 1 d ) G v ( q , q m d ) G v ( q , q 1 b ) G v ( q , q ( m + n ) b ) ] T [ π w 2 m u 1 e 1 ( cos α 1 + i sin α 1 ) + w 2 m 1 d L V 1 d π w 2 m m u e m ( cos α m + i sin α m ) + w 2 m m d L V m d L k 1 s x y U 1 b L k 1 s x x V 1 b L k n s x y U n b L k n s x x V n b ] m + n
m + n equations can be obtained based on Equation (11) using the unbalance responses as inputs and considering rotor unbalances and bearing coefficients as unknown variables. Equation (12) can be obtained by writing the equations in matrix form.
E I L 2 [ V ( q 1 s ) V ( q m s ) V ( q ( m + 1 ) s ) V ( q n s ) ] = H 3 [ π w 2 m u 1 e 1 ( cos α 1 + i sin α 1 ) + w 2 m 1 d L V 1 d π w 2 m m u e m ( cos α m + i sin α m ) + w 2 m m d L V m d L k 1 s x y U 1 b L k 1 s x x V 1 b L k n s x y U n b L k n s x x V n b ] m + n
where V ( q 1 s ) , …, V ( q m s ) and V ( q ( m + 1 ) s ) , …, V ( q ( m + n ) s ) are the measured dimensionless unbalance responses in the frequency domain in the y direction of locations, excluding the discs and bearings; the number of V ( q 1 s ) , …, V ( q m s ) is m and the number of V ( q ( m + 1 ) s ) , …, V ( q n s ) is n ; H 3 is a ( m + n ) × ( m + n ) matrix and can be calculated using Green’s functions G v ( q , q i ) as long as the locations of the measuring points, the bearings, and the discs are known; H 3 is as follows:
H 3 = [ G v ( q 1 s , q 1 d ) G v ( q 1 s , q m d ) G v ( q 1 s , q 1 b ) G v ( q 1 s , q n b ) G v ( q m s , q 1 d ) G v ( q m s , q m d ) G v ( q m s , q 1 b ) G v ( q m s , q n b ) G v ( q ( m + 1 ) s , q 1 d ) G v ( q ( m + 1 ) s , q m d ) G v ( q ( m + 1 ) s , q 1 b ) G v ( q ( m + 1 ) s , q n b ) G v ( q n s , q 1 d ) G v ( q n s , q m d ) G v ( q n s , q 1 b ) G v ( q n s , q n b ) ] .
Hence, Equation (13) is obtained.
[ π w 2 m u 1 e 1 ( cos α 1 + i sin α 1 ) + w 2 m 1 d L V 1 d π w 2 m m u e m ( cos α m + i sin α m ) + w 2 m m d L V m d L k 1 s x y U 1 b L k 1 s x x V 1 b L k n s x y U n b L k n s x x V n b ] m + n = E I L 2 H 3 1 [ V ( q 1 s ) V ( q m s ) V ( q ( m + 1 ) s ) V ( q n s ) ]
Define H 4 as following.
H 4 = E I L 2 H 3 1 [ V ( q 1 s ) V ( q m s ) V ( q ( m + 1 ) s ) V ( q n s ) ] .
According to Equations (13) and (14) can be obtained.
{ π w 2 m u 1 e 1 ( cos α 1 + i sin α 1 ) + w 2 m 1 d L V 1 d = H 4 ( 1 , 1 ) π w 2 m m u e m ( cos α m + i sin α m ) + w 2 m m d L V m d = H 4 ( m , 1 )
According to Equations (14) and (15) can be obtained.
{ m u 1 e 1 ( cos α 1 + i sin α 1 ) = [ H 4 ( 1 , 1 ) w 2 m 1 d L V 1 d ] π w 2 m m u e m ( cos α m + i sin α m ) = [ H 4 ( m , 1 ) w 2 m m d L V m d ] π w 2
According to Equations (15)–(17) are obtained, by means of which the amplitude and angle of the rotor unbalance can be calculated. The unbalance responses in the x direction are used as inputs. To calculate the matrix H 4 , the number of measured unbalance responses should be m + n according to Equations (12) and (13). To identify an eccentric disc’s unbalance, its unbalance response should also be measured. Therefore, the total amount of the measured unbalance responses is m + n + 1 and the unbalance responses of all the discs should be measured in order to estimate all the eccentric discs’ unbalances.
{ m u 1 e 1 = | [ H 4 ( 1 , 1 ) w 2 m 1 d L V 1 d ] π w 2 | m m u e m = | [ H 4 ( m , 1 ) w 2 m m d L V m d ] π w 2 |
{ α 1 = angle ( [ H 4 ( 1 , 1 ) w 2 m 1 d L V 1 d ] π w 2 ) α m = angle ( [ H 4 ( m , 1 ) w 2 m m d L V m d ] π w 2 )

2.3. Two Orthogonal Direction Algorithms

According to Equations (6) and (14), Equation (18) can be obtained.
{ H 2 ( 1 , 1 ) π w 2 m 1 u e 1 ( sin α 1 i cos α 1 ) H 4 ( 1 , 1 ) π w 2 m u 1 e 1 ( cos α 1 + i sin α 1 ) = U 1 d V 1 d = C 1 d H 2 ( m , 1 ) π w 2 m m u e m ( sin α m i cos α m ) H 4 ( m , 1 ) π w 2 m m u e m ( cos α m + i sin α m ) = U m d V m d = C m d
Let { c x 1 = m 1 u e 1 sin α 1 c y 1 = m 1 u e 1 cos α 1 { c x m = m m u e m sin α m c y m = m m u e m cos α m . According to Equations (18) and (19) is obtained.
{ H 2 ( 1 , 1 ) π w 2 ( c x 1 i c y 1 ) = C 1 d H 4 ( 1 , 1 ) C 1 d π w 2 m u 1 e 1 ( c y 1 + i c x 1 ) H 2 ( m , 1 ) π w 2 ( c x m i c y m ) = C m d H 4 ( m , 1 ) C m d π w 2 ( c y m + i c x m )
In Equation (19), H 2 ( 1 , 1 ) H 2 ( m , 1 ) , H 4 ( 1 , 1 ) H 4 ( m , 1 ) , C 1 d C m d are complex numbers and the others are real numbers. According to Equations (19) and (20) can be obtained.
{ π w 2 [ ( 1 + i C 1 d ) r C 1 d r C 1 d ( 1 + i C 1 d ) ] [ c x 1 c y 1 ] = [ r ( C 1 d H 4 ( 1 , 1 ) ) r H 2 ( 1 , 1 ) i ( C 1 d H 4 ( 1 , 1 ) ) i H 2 ( 1 , 1 ) ] π w 2 [ ( 1 + i C m d ) r C m d r C m d ( 1 + i C m d ) ] [ c x m c y m ] = [ r ( C m d H 4 ( m , 1 ) ) r H 2 ( m , 1 ) i ( C m d H 4 ( m , 1 ) ) i H 2 ( m , 1 ) ]
where r H 2 ( 1 , 1 ) r H 2 ( m , 1 ) , r H 4 ( 1 , 1 ) r H 4 ( m , 1 ) , r C 1 d r C m d are the real parts of the above-mentioned complex number; and i H 2 ( 1 , 1 ) i H 2 ( m , 1 ) , i H 4 ( 1 , 1 ) i H 4 ( m , 1 ) , i C 1 d i C m d are the imaginary parts of the above-mentioned complex number.
Hence, Equation (21) is obtained according to Equation (20).
{ [ c x 1 c y 1 ] = 1 π w 2 [ ( 1 + i C 1 d ) r C 1 d r C 1 d ( 1 + i C 1 d ) ] 1 [ r ( C 1 d H 4 ( 1 , 1 ) ) r H 2 ( 1 , 1 ) i ( C 1 d H 4 ( 1 , 1 ) ) i H 2 ( 1 , 1 ) ] [ c x m c y m ] = 1 π w 2 [ ( 1 + i C m d ) r C m d r C m d ( 1 + i C m d ) ] 1 [ r ( C m d H 4 ( m , 1 ) ) r H 2 ( m , 1 ) i ( C m d H 4 ( m , 1 ) ) i H 2 ( m , 1 ) ]
The amplitude and angle of the rotor unbalance are:
{ { m 1 u e 1 = abs ( c y 1 + i c x 1 ) α 1 = angle ( c y 1 + i c x 1 ) { m m u e m = abs ( c y m + i c x m ) α m = angle ( c y m + i c x m )
where abs ( ) is the function used to obtain the module of a complex number and angle ( ) is the function used to obtain the angle of a complex number.
Therefore, the amplitude and angle of each eccentric disc’s unbalance can be calculated based on Equations (21) and (22). There should be m + n unbalance responses in both the x and y directions to obtain the matrix H 2 and H 4 . If an eccentric disc’s unbalance is to be identified, its unbalance responses in both the x and y directions are required according to Equation (18). Hence, there should be m + n + 1 measured unbalance responses in the two orthogonal directions. To identify all the eccentric discs’ unbalances, the unbalance responses of the measured eccentric discs must be included.

2.4. Identification Procedures of the Two Algorithms

The identification procedures used in SDA and TODA to estimate the rotor unbalances ( m u j e j , α i ) of all discs are defined in Figure 2.
Firstly, each disc’s unbalance response and the other n + 1 unbalance responses should be measured and changed to dimensionless unbalance responses in the frequency domain according to Equation (23). Meanwhile, the rotating speed should also be measured. The inherent parameters, which are the length of the shaft, the mass per unit length of the rotor shaft, the elastic modulus of the shaft, and the diameter of the shaft, should be known, as prior knowledge and the location of the selected measured points on the shaft should also be used as inputs.
Secondly, the matrix H 1 and H 3 could be calculated according to Equations (4) and (12), respectively. Then, H 2 and H 4 can be calculated according to Equations (5) and Equation (13), respectively.
Thirdly, using H 2 and the dimensionless unbalance responses in the y direction in frequency domain obtained in the first step, each disc’s unbalance can be calculated according to Equations (8) and (10). Or, using H 4 and the dimensionless unbalance responses in the x direction in the frequency domain obtained in the first step, each disc’s unbalance can be calculated according to Equations (16) and (17). This is the identification procedure used in SDA. While for TODA, C1d, ……, Cmd should be calculated using the dimensionless unbalance responses in the frequency domain in both the x and y directions according to Equation (18). Then, using H 2 and H 4 , all rotor unbalances can be calculated based on Equations (20) and (21).
U D = A π L [ cos ( α ) + i sin ( α ) ]
where U D is the dimensionless unbalance response in the frequency domain; and A and α are the amplitude and phase of the unbalance responses in the time domain, respectively.

3. Numerical Simulations and Discussion

3.1. Methodology of the Numerical Simulations

The numerical simulations were conducted to validate the proposed methods by comparing the identified rotor unbalances with the set value of rotor unbalances.
As shown in Figure 3 and Figure 4, six computational examples, which represent single-span and single-disc rotors (g1.1, h1.1), single-span and four-disc rotors (g1.4, h1.4), and four-span and four-disc rotors (g4.4, h4.4) are used in the simulation. The rotors g1.1, g1.4, and g4.4 are supported by rolling bearings, and h1.1, h1.4, and h4.4 are sustained by oil journal bearings. The parameters of the rotors are summarized in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6. The positions of each bearing and disc on the shaft are shown in Figure 3 and Figure 4. There are 121 nodes in the rotors in Figure 3 by dividing the shaft into 120 segments equally. There are 61 nodes in the rotors in Figure 4 by dividing the shaft into 60 segments equally. All the nodes, where the bearings and the discs located, are chosen for calculating simulated unbalance response by the CRDAM. Moreover, besides the above points, any other point on the shaft is also chosen in the simulation.
Four kinds of numerical simulations were conducted. Firstly, the simulated unbalance responses calculated by the CRDAM were directly fed into SDA and TODA to estimate rotor unbalances. Secondly, a similar identification exercise was performed by contaminating simulated unbalance responses by the set measured error. The relative error for the unbalance response amplitude was 5% and the absolute error for the unbalance response angle was 5°. Thirdly, the resolution of the vibration displacement sensor should be considered when measuring the unbalance response. Hence, identification exercises were performed by contaminating simulated unbalance responses at a resolution of 0.1 nm, which is the highest used at present. Finally, identification exercises were performed considering four kinds of typical sensor resolution (0.1 nm, 1 nm, 0.1 um, and 1 um). By limiting the number of digits after the decimal point in the unbalance responses, the resolutions of unbalance response measurement systems are considered.
In the above simulations, the calculation frequency was from 1 to 2001 Hz and the interval was 2 Hz. The relative error for the identified rotor unbalance amplitude to the set rotor unbalance amplitude and the absolute error for the identified rotor unbalance angle to the set rotor unbalance angle were used for the analysis of the simulation results. Moreover, the statistical results of the allowable frequency points (AFPs), at which the identification error for the rotor unbalance amplitude was less than 20% and the rotor unbalance angle was less than 10°, were obtained.
In addition, as shown in Figure 2, SDA includes two algorithms which are represented by Equations (8), (10), (16), and (17), respectively. The difference between them is that one requires unbalance responses in the y direction as input and the other uses unbalance responses in the x direction. They can be considered the same algorithm. Hence, for SDA, only Equations (8) and (10), which require unbalance responses in the y direction, are validated by the simulations.

3.2. Accuracy of SDA and TODA

3.2.1. Results

(1)
Results of the first kind of simulation
Figure 5 and Figure 6 were obtained by using the computational example h1.4 based on SDA and TODA, respectively. They represent the maximum identification errors of each rotor unbalance. In the simulation, #10 point was used as one of the required measuring points and the other measuring points were #3, #47, #11, #21, #31, and #41 points at the location of the two bearings and the four discs. The form α%∠β°, in which α% is the identification error for the rotor unbalance amplitude and β° represents the identification error for the rotor unbalance angle, is used to express the identification error for rotor unbalance.
According to Figure 5, for SDA, the maximum relative errors of the amplitude of rotor unbalance for #1-#4 discs are almost equal to zero and are 1.70545 × 10−8%, 1.4489 × 10−6%, 3.22562 × 10−5%, and 2.8886 × 10−5%, respectively. For the rotor unbalance angle, the maximum absolute errors, which are 5.505.9 × 10−9°, 5.99272 × 10−7°, 5.83537 × 10−5°, and 1.93098 × 10−5°, respectively, also equal zero.
According to Figure 6, for TODA, the maximum relative errors of the rotor unbalance amplitude of each disc are 1.14223 × 10−6%, 2.66709 × 10−5%, 0.00102%, and 1.6262 × 10−4%, respectively. The maximum absolute errors of the rotor unbalance angle are 4.68235 × 10−7°, 1.59148 × 10−5°, 4.41457 × 10−4°, and 5.47386 × 10−5°, respectively. Although they are bigger than the maximum errors obtained by SDA, they can be also considered equal to zero.
Similar results can be found for the first kind of calculation simulation example for g1.1, h1.1, h4.4, and g1.4, according to Table A1 and Figure A1, Figure A2, Figure A3 and Figure A4 in Appendix A.
However, according to Figure 7 and Figure 8, which show the maximum identification errors of each rotor unbalance for g4.4, for TODA, the first two errors are very small, but the maximum identification errors of rotor unbalance for #3 and #4 discs are too big, according to Figure 8. The first two errors are 7.52405 × 10−5%∠1.81997 × 10−4°, 0.17108%∠0.03405°, while for the errors for #3 and #4 discs, they are 123.71543%∠21.8563° and 1281.45771%∠21.69524°, respectively. Whereas, for SDA, the maximum identification errors of rotor unbalance for #1-#4 discs are 4.70411 × 10−8%∠1.49579 × 10−7°, 5.00642 × 10−4%∠1.41912 × 10−4°, 0.28694%∠0.12378°, and 0.16497%∠1.57779°, respectively, according to Figure 8.
Moreover, Figure 9 represents the rotor unbalance identification errors for #3 and #4 discs from 1 to 2001 Hz obtained. The maximum identification error for rotor unbalance for #3 disc occurs at 1 Hz. For #4 disc, the maximum identification errors of the rotor unbalance amplitude and angle occur at 1 Hz and 2 Hz, respectively. However, when the rotating frequency becomes high, the identification errors are almost equal to zero.
(2)
Results of the second kind of simulation
Figure 10 and Figure 11, which represent the maximum identification errors of each rotor unbalance for h1.4 are obtained based on SDA and TODA, respectively. In this simulation, #10 point was used as one of the required measuring points and the other measuring points were #3, #47, #11, #21, #31, and #41 at the locations of the two bearings and the four discs.
According to Figure 10, for SDA, the maximum relative errors of the rotor unbalance for #1–#4 discs are 5%∠5.000001199°, 5%∠5.000020523°, 5.00003%∠5.00022366°, and 5%∠5.000170751°, respectively, which are almost equal to the setting errors (5%∠5°) of the inputted unbalance responses. According to Figure 11, for TODA, the maximum relative errors of rotor unbalance for each disc also almost equal the setting error. They are 5%∠5°, 5%∠5.00002°, 5.00006%∠5.00008°, and 5%∠5.00005°, respectively. Similar results can be found in the second kind of calculation simulation example for g1.1, h1.1, h4.4, and g1.4, according to Table A2 and Figure A5, Figure A6, Figure A7 and Figure A8 in Appendix B.
However, according to Figure 12 and Figure 13, which show the maximum identification errors of each rotor unbalance for g4.4, for TODA, the maximum identification errors of rotor unbalance for #3 and #4 discs are 134.9012142%∠16.8563° and 1350.530732%∠26.69523°, respectively, although the errors for #1 and #2 discs are almost equal the setting error 5%∠5°. Whereas, for SDA, the maximum identification errors for the rotor unbalance for #1–#4 discs are 5%∠5°, 5.00001%∠5.00014°, 5.0039%∠5.00404°, and 5.17322%∠6.57779°, respectively. Only the rotor unbalance angle of #4 disc is a bit bigger than the setting error of 5%. The others are almost equal to the input setting error of 5%∠5°.
Hence, the rotor unbalance identification errors for #3 and #4 discs from 1 to 2001 Hz are as shown in Figure 14. The maximum error occurs at the beginning. The maximum identification error for rotor unbalance for #3 disc is at 1 Hz. For #4 disc, the maximum identification errors of rotor unbalance amplitude and angle are at 1 Hz and 2 Hz, respectively. However, when the rotating frequency becomes high, the identification errors are almost equal the input setting error of 5%∠5°.

3.2.2. Discussion

In the first kind of numerical simulation, the simulated unbalance responses calculated by the CRDAM were fed into SDA and TODA to estimate the rotor unbalance. For SDA, the identified rotor unbalance was equal to the set value of rotor unbalance. For TODA, the identified rotor unbalance was equal to the set value of rotor unbalance except for the simulation of g4.4. The error for the identified rotor unbalance for #3 and #4 discs in g4.4 is big at the low rotating frequency (1 Hz, 2 Hz) but very small at other rotating frequencies.
In the second kind of numerical simulation, the simulated unbalance responses, which were calculated by the CRDAM and contaminated the set measured error (5%, 5°), were fed into SDA and TODA. For SDA, the identified error was equal to the set error. For TODA, the identified rotor unbalance was almost equal to the set error except for the simulation of g4.4. The error for the identified rotor unbalance for #3 and #4 discs in g4.4 was bigger than the set errors at the low rotating frequency (1 Hz, 2 Hz) but almost equal to the set error at other frequencies.
The unbalance responses in both the x and y directions were needed when using TODA. At a low frequency, the value of the unbalance response is very small. The errors in computer calculations, such as rounding errors, have a considerable influence on small values. Accordingly, there were errors with regard to the unbalance responses in both the x and y directions at low frequencies. If the errors in the x and y  directions are quite different, a big identification error may occur at a low frequency for TODA because the errors in the x direction are divided by the errors in the y direction, while, for SDA, the unbalance response in only one direction is needed. Hence, the identification error can be equal to the set measured error.
Therefore, the error for the identified rotor unbalance response will be equal to the set error if the errors of all the measured unbalance responses are the same. It is indicated that the repeatability precision of each measuring channel of the unbalance response measurement system plays a critical role when it comes to using the proposed algorithms. In engineering, the vibration caused by rotor unbalance can be very small for low-speed rotors and the measuring system cannot accurately detect tiny vibrations. Therefore, the proposed method, especially TODA, cannot be applied to low-speed rotors.

3.3. Adjustment Point

3.3.1. Results

In the third kind of simulation of g1.4, the matrix of the maximum identified error for rotor unbalance for each disc was obtained with Equations (24)–(27). The numbers of AFPs were also obtained. They are shown in Figure 15, Figure 16, Figure 17 and Figure 18. The simulations were conducted by using simulated unbalance responses as input data to SDA and TODA and by changing one measuring point which was close to the identified disc. The #20, #40, #60, and #80 points were used for the respective changes. The other m + n measuring points were at the discs and bearings. The simulated unbalance responses were calculated by means of the CRDAM contaminating the sensor resolution at 0.1 nm.
In Equations (24)–(27), the four element values in a row of the matrix are the maximum identification errors of unbalance for #1-#4 discs, respectively, when using the same measuring point. The four elements in a column are the maximum identification errors of unbalance for the same disc when the measuring points #20, #40, #60, and #80 were applied, respectively. By comparing the biggest identification error obtained under different adjustment point conditions when the rotating frequency is changed from 1 to 2001 Hz, it is easy to find that the adjustment point plays a critical role in improving the identification accuracy for SDA and TODA. The details are as follows:
(1)
According to Equations (24) and (25), for SDA, when #20 point, which is near #1 disc, is used as one of the required measuring points, the maximum identification error for rotor unbalance for #1 disc is the smallest among that of the four discs. The maximum identification error for rotor unbalance is 10.84%∠2.12°, while for #2–4 discs, their maximum errors are much bigger. When #40 point, which is close to #2 disc is used, the maximum identification error for rotor unbalance for #2 disc becomes the smallest. It is 42.54%∠8.60°. When #60 point, which is beside #3 disc, is used, #3 disc’s rotor unbalance identification error, which is 257.24%∠57.41°, is the second smallest and is close to the smallest error of 254.72%∠51.14°. When #80 point, which is beside #4 disc, is used, #4 disc’s rotor unbalance identification error, which is 9.82%∠16.25°, is much smaller than the others. From the perspective of the column in the matrix, when the measuring point is changed from #20 to #40, the maximum identification error for rotor unbalance for #1 disc becomes bigger, while the maximum identification error for rotor unbalance for #4 disc becomes smaller. The maximum identification errors for rotor unbalance for #2 and #3 discs are also changed apparently. As #20 point is applied, #1 disc acquires the best recognition accuracy, while #2 disc has the best recognition accuracy as #40 point is applied, #3 disc has the best recognition accuracy as #60 point is applied, and #4 disc has the best recognition accuracy as #80 point is applied.
(2)
As for TODA, similar results can be obtained, although the values of maximum identification errors are much bigger than the maximum identification errors for SDA, according to Equations (26) and (27). When #20 point, which is near #1 disc, is used as one of the required measuring points, the maximum identification error for rotor unbalance for #1 disc is the smallest among that of the four discs. The maximum identification error for rotor unbalance is 4.22 × 102%∠170.38°, while for #2–4 disc, their maximum errors are much bigger. When #40 point, which is close to #2 disc, is used, the maximum identification error for rotor unbalance for #2 disc becomes the smallest. It is 1.64 × 104%∠140.73°. When #60 point, which is beside #3 disc, is used, #3 disc’s rotor unbalance identification error, which is 6.03 × 104%∠179.57°, is the second smallest. When #80 point, which is beside #4 disc, is used, #4 disc’s rotor unbalance identification error is 2.68 × 1015%∠179.95°. Although it is the biggest among the obtained maximum identification errors of the four discs, it is the smallest error for #4 disc, which can be obtained by changing the measuring point from #20 to #40. From the perspective of the column in the matrix, when the measuring point is changed from #20 to #40, the maximum identification error for rotor unbalance for #1 disc becomes bigger. The maximum identification errors of rotor unbalance for #2, #3, and #4 discs are also changed apparently. As #20 point is applied, #1 disc acquires the best recognition accuracy, while #2 disc has the best recognition accuracy as #40 point is applied, #3 disc has the best recognition accuracy of unbalance amplitude as #60 point is applied, and #4 disc has the best recognition accuracy as #80 point is applied.
S D A _ M 1 g 14 = [ 10.84 376.18 14934.67 22570.13 400.78 42.54 7251.41 4823.78 382.89 254.72 257.54 3203.79 1153.90 475.79 1148.01 9.82 ]
S D A _ M 2 g 14 = [ 2.12   152.05 179.19 170.46 83.28 8.60 178.20 173.14 132.99 51.14 57.41 173.82 121.58 141.59 139.37 16.25 ]  
T O D A _ M 1 g 14 = [ 4.22 × 10 2 1.54 × 10 5 1.04 × 10 7 7.39 × 10 18 4.97 × 10 4 1.46 × 10 4 3.45 × 10 6 9.50 × 10 17 3.66 × 10 4 6.58 × 10 4 6.03 × 10 4 1.37 × 10 17 9.14 × 10 5 1.94 × 10 6 3.89 × 10 6 2.68 × 10 15 ]
T O D A _ M 2 g 14 = [ 170.38   179.83 178.53 179.98 178.23 140.73 179.01 180.00 179.27 179.76 179.57 179.99 179.75 179.67 179.67 179.95 ]  
where S D A _ M 1 g 14 and S D A _ M 2 g 14 are the matrix of the maximum identification error for rotor unbalance amplitude and the angle of #1 to #4 discs under different adjustment point conditions using SDA; T O D A _ M 1 g 14 and T O D A _ M 2 g 14 are the matrix of the maximum identification error for rotor unbalance amplitude and the angle of #1 to #4 discs under different adjustment point conditions using TODA.
According to the number of AFPs of rotor unbalance obtained by SDA in Figure 15, Figure 16, Figure 17 and Figure 18, the results are as follows:
(1)
The numbers of AFPs of rotor unbalance amplitude for #1–#4 discs are 1001, 979, 744, and 762, respectively, according to Figure 15a. For the rotor unbalance angle, there are 1001, 990, 629, and 612 AFPs, respectively, according to Figure 15b. This indicates that the identification results for rotor unbalance for #1 disc are best when #20 point is used as one of the required measuring points.
(2)
According to Figure 16, there are 994, 999, 824, and 929 AFPs, respectively, for the rotor unbalance amplitude of #1 disc, #2 disc, #3 disc, and #4 disc, and 988, 1001, 835, and 787 AFPs, respectively, for the rotor unbalance angle. This indicates that the identification results for rotor unbalance for #2 disc are best when #40 point is used as one of the required measuring points.
(3)
According to Figure 17, the numbers of AFPs of rotor unbalance amplitude are 980, 987, 995, and 988, respectively, and the numbers of rotor unbalance angles are 974, 995, 995, and 960 respectively. This indicates that the identification results for rotor unbalance for #3 disc are best when #60 point is used as one of the required measuring points.
(4)
According to Figure 18, for the rotor unbalance amplitude of #1 disc, #2 disc, #3 disc, and #4 disc, there are 917, 929, 897, and 1001 AFPs, respectively. For the angle, there are 861, 966, 881, and 999 AFPs, respectively. This indicates that the identification results for rotor unbalance for #4 disc are best when #80 point is used as one of the required measuring points.
For TODA, the same conclusion can be obtained according to Figure 15, Figure 16, Figure 17 and Figure 18, although the number of AFPs obtained by TODA is obviously lesser than that of SDA. When #20 point is used as one of the required measuring points, the numbers of AFPs of rotor unbalance amplitude for #1-#4 discs are 867, 163, 29, and 12, respectively. For the rotor unbalance angle, there are 896, 117, 57, and 186 AFPs, respectively. Apparently, #1 disc has the most AFPs. When #40 point is used, the numbers of AFPs of rotor unbalance amplitude for #1-#4 discs are 587, 662, 107, and 38, respectively. For the rotor unbalance angle, there are 541, 720, 137, and 158 AFPs, respectively. Apparently, #2 disc has the most AFPs. When #60 point is used, the number of AFPs of rotor unbalance amplitude of #1-#4 discs are 408, 297, 461, and 234, respectively. For the rotor unbalance angle, there are 352, 271, 421, and 201 AFPs, respectively. Obviously, #3 disc has the most AFPs. When #80 point is used, the number of AFPs of rotor unbalance amplitude of #1-#4 discs are 78, 48, 39, and 707, respectively. For the rotor unbalance angle, there are 166, 128, 97, and 675 AFPs, respectively. Obviously, #4 disc has the most AFPs.
Similar results can be found in the third kind of calculation simulation example for g4.4, h1.4, and h4.4, according to Equations (A1)–(A12) and Figure A9, Figure A10, Figure A11, Figure A12, Figure A13, Figure A14, Figure A15, Figure A16, Figure A17, Figure A18, Figure A19 and Figure A20 in Appendix C. However, for h4.4, no result can be obtained based on SDA and TODA from 1 Hz to 20 Hz.

3.3.2. Discussion

In the third kind of numerical simulation, the highest sensor resolution (0.1 nm) at present is considered. In the simulation of g1.4 and g4.4, when #20 point, which is near #1 disc, is used as one of the required measuring points, the identification accuracy for the rotor unbalance for #1 disc (at location of #21 point) is high. Similarly, #2 disc at #41 point, #3 disc at #61 point, and #4 disc at #81 point can obtain high identification accuracy for rotor unbalance as #40, #60, and #80 points are used respectively. In the simulation of h1.4 and h4.4, #1disc at #11 point, #2 disc at #21 point, #3 disc at #31point, and #4 disc at #41 point can obtain a high identification accuracy for rotor unbalance as #10, #20, #30, and #40 points are used, respectively.
From the above, to obtain good identification results, there should be a measuring point near the disc, for which the rotor unbalance should be identified. The measuring point near the disc determines the identification error for the rotor unbalance and it should be near the disc the unbalance of which is to be identified. Hence, the measuring point is called an adjustment point.
By applying the proposed adjustment point, the identification errors for rotor unbalance for #3 and #4 discs in the first and second kind of simulations are reduced and Figure 19 and Figure 20 are obtained. According to Figure 19, the maximum identification error is almost equal to zero. It is 0.00358%∠0.00193° for #3 disc and 2.04046 × 10−4%∠5.84276 × 10−5° for #4 disc, while they are 123.71543%∠21.8563° and 1281.45771%∠21.69524°, respectively, before the application of the adjustment point. According to Figure 20, the ultimate error is 5.0006%∠5.00193° for #3 disc and 4.99979%∠5.00377° for #4 disc. They are almost equal to the setting error of 5%∠5°, while they are 134.9012142%∠16.8563° and 1350.530732%∠26.69523°, respectively, before the application of the adjustment point. Therefore, there should be a sensor mounted near each disc if the unbalances of all the discs are to be identified.
It was established that SDA has a much better identification accuracy than TODA. The reason is that the unbalance responses in both the x and y directions are involved in TODA and the unbalance response in only one direction is needed for SDA. This means that the measured errors of unbalance responses in both the x and y directions will contribute to the identification error of TODA. Moreover, they are combined by division in TODA. When the measured errors in the two directions are not the same and the error, whose value is smaller, happens to be the denominator, a high identification error occurs. Based on the simulation results, the value of the denominator is smaller in most cases, but the value of the denominator can be bigger, in which case the identification error of SDA is bigger than that of TODA, according to Figure A9, Figure A17 and Figure A20. Moreover, according to Figure A9, Figure A17 and Figure A20, similar results for SDA and TODA can be obtained when the comprehensive effect of the response measuring errors in the x and y directions is equal to that in the x or y direction. It is indicated that the identification result of TODA is greatly affected by the measured error for the unbalance response due to the unbalance responses in both the x and y directions being required. In practical engineering, it is difficult to obtain completely consistent measurement errors for each unbalance response. The measured errors for unbalance responses in the x direction are different from those in the y direction. This may lead to big identification errors for TODA. Hence, SDA is better than TODA from the perspective of engineering applications.

3.4. Affect of Sensor Resolution

3.4.1. Results

According to the result in the third kind of simulation, the maximum identification error of TODA was too big. Hence, the fourth kind of simulation was conducted based on SDA. The simulated unbalance responses calculated by CRDAM containing different sensor resolutions (0.1 nm, 1 nm, 0.1 um, and 1 um) were used as input data to the method. The computational example g1.4 is used. In the simulation, the measuring points for the two discs and the four bearings are utilized and #20, #40, #60, and #80 points are used as the respective adjustment points in order to obtain good identification results.
The curves for the identification errors, which were obtained when the sensor resolutions were 0.1 nm, 1 nm, 0.1 um, and 1 um, are shown in Figure 21, Figure 22, Figure 23 and Figure 24. According to the figures, the identified rotor unbalance changes with the rotating frequencies and there are several peak values for the identified rotor unbalance. The results are as follows:
(1)
Sensor resolution of 0.1 nm
According to Figure 21, for #1 disc, the peak values of the amplitude of the identified rotor unbalance are 10.84%, 2.02%, 0.31%, and 0.42% at frequencies of 1 Hz, 271 Hz, 643 Hz, and 1267 Hz, respectively. The peak values of the absolute error for the angle are 2.12°,1.47°,0.32°, and 0.32° at 1 Hz, 271 Hz, 645 Hz, and 1267 Hz, respectively.
For #2 disc, the peak values of the amplitude of the identified rotor unbalance are 3.46%, 42.54%, 4.33%, and 1.42% at frequencies of 1 Hz, 269 Hz, 645 Hz, and 1263 Hz, respectively. The peak values of the absolute error for the angle are 2.29°, 8.6°, 0.31°, and 0.8° at 1 Hz, 267 Hz, 645 Hz, and 1265 Hz, respectively.
For #3 disc, the peak values of the amplitude of identified rotor unbalance are 26.02%, 257.24%, 11.64%, and 20.64% at frequencies of 1 Hz, 269 Hz, 645 Hz, and 1265 Hz, respectively. The peak values of the absolute error for the angle are 19.49°, 57.41°, 11.72°, and 7.76° at 1 Hz, 271 Hz, 645 Hz, and 1267 Hz, respectively.
For #4 disc, the peak values of the amplitude of identified rotor unbalance are 9.82%, 6.46%, 0.82%, 0.56%, and 7.61% at frequencies of 1 Hz, 267 Hz, 645 Hz, 809 Hz, and 1267 Hz, respectively. The peak values of the absolute error for the angle are 16.259°, 1.88°, 1.39°, and 4.15° at 269 Hz, 645 Hz, 809 Hz, and 1265 Hz, respectively.
(2)
Sensor resolution of 1 nm
According to Figure 22, for #1 disc, the peak values of the amplitude of identified rotor unbalance are 83.36%, 1.92%, 0.31%, and 0.42% at frequencies of 1 Hz, 269 Hz, 643 Hz, and 1267 Hz, respectively. The peak values of the absolute error for the angle are 8.77°, 1.45°, 0.32°, and 0.32°at 1 Hz, 271 Hz, 645 Hz, and 1267 Hz, respectively.
For #2 disc, the peak values of the amplitude of identified rotor unbalance are 95.6%, 42.54%, 4.33%, and 1.42% at frequencies of 1 Hz, 269 Hz, 645 Hz, and 1263 Hz, respectively. The peak values of the absolute error for the angle are 32.5°, 8.6°, and 0.8° at 1 Hz, 267 Hz, and 1265 Hz, respectively.
For #3 disc, the peak values of the amplitude of identified rotor unbalance are 242.24%, 257.23%, 11.64%, and 20.64% at frequencies of 1 Hz, 269 Hz, 645 Hz, and 1265 Hz, respectively. The peak values of the absolute error for the angle are 108.42°, 57.41°, 11.72°, and 7.76°at 1 Hz, 271 Hz, 645 Hz, and 1267 Hz, respectively.
For #4 disc, the peak values of the amplitudes of identified rotor unbalance are 17.48%, 6.46%, 0.82%, 0.66%, and 7.61% at frequencies of 1 Hz, 267 Hz, 645 Hz, 809 Hz, and 1267 Hz, respectively. The peak values of the absolute error for the angle are 165.08°, 16.25°, and 4.15°at 1 Hz, 269 Hz, and 1265 Hz, respectively.
(3)
Sensor resolution of 0.1 um
With the sensor resolution being reduced to 0.1 um, the minimum starting frequency of the identification was 4 Hz (SDA does not work at a low frequency). According to Figure 23, the results are as follows:
For #1 disc, the peak values of the amplitude of identified rotor unbalance are 695.26%, 19.45%, 0.65%, and 0.62% at frequencies of 4 Hz, 268 Hz, 646 Hz, and 1266 Hz, respectively. The peak values of the absolute error for the angle are 140.04°, 10.86°, 0.55°and 0.45° at 4 Hz, 268 Hz, 646 Hz and 1262 Hz, respectively.
For #2 disc, the peak values of the amplitude of identified rotor unbalance are 597.53%, 22.7%, 5.32% and 1.63% at frequencies of 4 Hz, 270 Hz, 644 Hz, and 1268 Hz, respectively. The peak values of the absolute error for the angle are 94.44°, 5.72°1.9°, and 0.71° at 4 Hz, 266 Hz, 644 Hz, and 1268 Hz, respectively.
For #3 disc, the peak values of the amplitude of identified rotor unbalance are 1722.04%, 262.49%, 39.32%, 9.24%, and 11.76% at frequencies of 8 Hz, 268 Hz, 644 Hz, 814 Hz, and 1264 Hz, respectively. The peak values of the absolute error for the angle are 120.31°, 20.6°, 8.75°, 3.49°, and 11.3° at 16 Hz, 270 Hz, 644 Hz, 814 Hz, and 1268 Hz, respectively.
For #4 disc, the peak values of the amplitude of identified rotor unbalance are 853.26%, 4.43%, 13.87%, and 3.3% at frequency 6 Hz, 268 Hz, 810 Hz, and 1266 Hz, respectively. The peak values of the absolute error for the angle are 165.1°,10.04°, 0.79°, 4.34°, and 1.5° at 10 Hz, 270 Hz, 644 Hz, 810 Hz, and 1264 Hz, respectively.
(4)
Sensor resolution of 1 um
As the sensor resolution was reduced to 1 um, the starting frequency of the identification was 856 Hz (SDA does not work at a low frequency). According to Figure 24, for #1 disc, the peak value of the amplitude of identified rotor unbalance is 0.63% at a frequency of 1266 Hz. The peak value of the absolute error for the angle is 0.45° at 1262 Hz. For #2 disc, the peak value of the amplitude of identified rotor unbalance is 1.58% at a frequency of 1268 Hz. The peak value of the absolute error for the angle is 0.75° at 1268 Hz. For #3 disc, the peak value of the amplitude of identified rotor unbalance is 11.73% at a frequency of 1264 Hz. The peak value of the absolute error for the angle is 11.56° at 1268 Hz. For #4 disc, the peak value of the amplitude of identified rotor unbalance is 3.47% at a frequency of 1268 Hz. The peak value of the absolute error for the angle is 1.49° at 1264 Hz.
(5)
According to the above results, the peak values that were outside the assumed allowable range—the amplitudes bigger than 20% and the angles bigger than 10°—are all counted in Table 7, Table 8, Table 9 and Table 10, with the sensor resolutions set at 0.1 nm, 1 nm, 0.1 um, and 1 um, respectively. From Table 7, Table 8, Table 9 and Table 10, the frequencies at which these peak values appear are at or around the critical frequencies listed in Table 11 and the low frequencies such as 1 Hz, 4 Hz, and so on. This indicates that the rotor unbalance identification error is big when the rotor works at a low speed or near the critical speeds when using SDA.
Moreover, similar results can also be obtained in simulations of g4.4, h1.4, and h4.4.

3.4.2. Discussion

In the fourth kind of numerical simulation, the proposed adjustment point is applied and four kinds of typical sensor resolutions are considered for the validation of SDA. At or around the critical frequencies which are obtained by CRDAM and low frequencies, the identification error for rotor unbalance is bigger than at other frequencies. The resolution of sensors is applied by limiting the number of digits after the decimal point in the unbalance responses. At a low frequency, the unbalance responses are very small. The rounding digits have a considerable influence on small values and causes big input errors for unbalance responses. At or near the critical frequency, the unbalance responses vary greatly with frequencies. It is indicated that the rounding digits may lead to big differences among the measured errors of unbalance responses. The repeatability precision of the input responses is poor. Hence, big identification error occurs at low frequencies or at critical frequencies. Similarly, the peak values that are outside the assumed allowable range are at or around the critical frequencies and low frequencies. Far from the critical frequencies, the identification errors are small. The rotating machines’ operational frequency (speed) is designed to be far from the critical frequencies to avoid resonance. This indicates that SDA can be applied to rotors not working under very low rotating speeds.
The statistical results for peak values of the identification error for #1 disc of rotor g1.4 are listed in Table 12. According to Table 12, it is indicated as follows:
(1)
At a low frequency, the peak values of identification errors become bigger as sensor resolution is reduced and rotor unbalance cannot be identified when the sensor resolution is 1 um.
(2)
At or around the first order critical frequency, the identification error, which is obtained using a sensor resolution of 0.1 nm, is little different from the identification error obtained using a sensor resolution of 1 nm. The identification error becomes much bigger when the sensor resolution is 0.1 um and the rotor unbalance cannot be identified when the sensor resolution is close to 1 um.
(3)
At or around the second order critical frequency, the identification error obtained when the sensor resolution is 0.1 nm is equal to the identification error obtained when the sensor resolution is 1 nm. However, when the sensor resolution is 0.1 um, the identification error becomes a little bigger. Moreover, when the sensor resolution comes close to 1 um, the rotor unbalance cannot be identified.
(4)
At or around the third order critical frequency, the identification error obtained when the sensor resolution is 0.1 nm is equal to the identification error obtained when the sensor resolution is 1 nm. However, when the sensor resolution comes close to 0.1 um and 1 um, the identification errors become a little bigger.
The above indicates that sensors with a resolution of 1 um should be avoided and sensors with a resolution of 1 nm and 0.1 nm are recommended for practical application. Similar results can also be obtained from the identification errors for rotor unbalance for #2~#4 discs. These indicate that SDA requires high sensor resolution. The higher the sensor resolution is, the smaller the identification error. In engineering, the resolution of most vibration displacement sensors is 1 um and sensors with resolutions of 0.1 nm and 1 nm are very expensive. Hence, the cost of SDA can be very high when there are multiple discs and bearings in a rotor.

4. Conclusions

In this paper, two kinds of algorithms, SDA and TODA, are proposed to identify rotor unbalances in multi-disc and multi-span rotors from unbalance responses. Rotors are modeled based on the continuous rotor dynamic analysis method. By using the unbalance responses as inputs and considering the unbalance responses at the bearings and the eccentric discs as unknown variables, an inverse problem model is developed based on the matrix method to eliminate the coupling between bearing coefficients and rotor unbalances. Four kinds of numerical simulations considering sensor resolutions and measurement errors were conducted to study the algorithms. Compared with existing methods, the proposed algorithms have the flexibility to incorporate any number of bearings and discs. Moreover, test runs and external exciters are not required. The conclusions are summarized as follows.
(1)
The proposed algorithms provide a technique with which to monitor all the rotor unbalances on-line under operational conditions. For SDA, the unbalance responses in only one direction are needed. For a rotor with m discs and n bearings, the number of required unbalance responses is m + n + 1. While for TODA, the unbalance responses in both the x and y directions are required. The necessary measuring position of the two methods should be at the disc whose unbalance is to be monitored. For monitoring the unbalance of all discs, there should be a sensor mounted at each disc, while the other measuring positions can be at or around the bearings or discs. Moreover, numerical simulations indicate that there should be one measuring point called the adjustment point to achieve a high identification accuracy. The proposed adjustment point should be near the disc whose unbalance is to be monitored. In order to identify all the discs’ unbalances accurately and simultaneously, there should be m + n + n measuring points, among which n adjustment points near each disc are necessary.
(2)
The identification accuracy of the proposed algorithms requires a high performance of the unbalance response measurement system. Numerical simulations indicate that if the measuring errors of all the required unbalance responses are zero, the identification error will be zero, too. When the measuring errors are the same, the identification error will be equal to the measuring errors. It is indicated that the consistency of each channel’s measurement errors plays a critical role in identifying rotor unbalance when using the proposed algorithms. In addition, SDA has a better identification accuracy than TODA when considering sensor resolution from the perspective of engineering applications. The identification error of SDA is only high at a very low frequency and the critical frequencies when sensor resolutions are considered. Hence, SDA is suitable for medium-speed and high-speed rotors. Moreover, identification accuracy is strongly related to sensor resolution. Sensors with resolutions of 1 um should be avoided and sensors with 1 nm and 0.1 nm resolutions are recommended.
For future research, there are various possible directions. Experiments should be conducted to validate the proposed methods. The limitation of the proposed algorithms is that they are not suitable for low-speed rotors and a high accuracy of measurement of unbalance responses is strongly demanded. Therefore, sensors and instruments with high accuracy and good resolutions for measuring unbalance responses should be developed. Moreover, further investigations could focus on decreasing the high requirements on measurement accuracy. Sensors with a resolution of 1 um can be used for the identification. The Timoshenko model, in which gyroscopic moments are considered, could be taken into consideration based on CRDAM when modelling a continuous shaft. The research method of this paper can be regarded as a tool for future study.

Author Contributions

Conceptualization, A.W.; methodology, A.W.; software, A.W., Y.B. and Y.F.; formal analysis, A.W.; investigation, Y.F., X.C. and J.Y.; writing—original draft preparation, A.W. and Y.B.; writing—review and editing, A.W. and Y.B.; supervision, G.M.; funding acquisition, G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities, grant number 8000150A084, and the Yue Qi Young Scholar Project, China University of Mining & Technology, Beijing.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

q Dimensionless value of position z due on the z-axis
U ( q ) Dimensionless unbalance response in the frequency domain due on the y-axis
V ( q ) Dimensionless unbalance response in the frequency domain due on the x-axis
# For instance, #1 disc means the first disc
m j u Eccentric masses of # j disc
e j Eccentric distance of # j disc
α j Eccentric angle of # j disc
m j d Masses of # j disc
w Rotation frequency rotors
L Length of rotor shaft
E Elastic modulus of rotor shaft
I Diametric shaft cross-sectional geometric moment of inertia
k j j x x # j j bearing’s main stiffness coefficient in the x direction
k j j x y # j j bearing’s cross-coupled stiffness coefficient in the x direction
k j j y x # j j bearing’s cross-coupled stiffness coefficient in the y direction
k j j y y # j j bearing’s main stiffness coefficient in the y direction
c j j x x # j j bearing’s main damping coefficient in the x direction
c j j x y # j j bearing’s cross-coupled damping coefficient in the x direction
c j j y x # j j bearing’s cross-coupled damping coefficient in the y direction
c j j y y # j j bearing’s main damping coefficient in the y direction
k s x x Bearing’s main complex coefficient in the x direction
k s x y Bearing’s cross-coupled complex coefficient in the x direction
k s y x Bearing’s cross-coupled complex coefficient in the y direction
k s y y Bearing’s main complex coefficient in the y direction
z j d z coordinate position of # j disc
z j j b z coordinate position of # j j bearing
q j d Dimensionless value of z j d
q j j b Dimensionless value of z j j b
G u ( q , q i ) Green’s functions in the y direction
G v ( q , q i ) Green’s functions in the x direction
G u ( q , q j d ) Green’s coefficients for # j disc in the y direction
G u ( q , q j j b ) Green’s coefficients for # j j bearing in the y direction
G v ( q , q j d ) Green’s coefficients for # j disc in the x direction
G v ( q , q j j b ) Green’s coefficients for # j j bearing in the x direction
U j d # j disc’s dimensionless unbalance response in frequency domain in the y direction
V j d # j disc’s dimensionless unbalance response in frequency domain in the x direction
U j j b # j bearing’s dimensionless unbalance response in frequency domain in the y direction
V j j b # j bearing’s dimensionless unbalance response in frequency domain in the x direction
m Number of discs
n Number of bearings

Appendix A

(1)
Simulation results for g1.1 and h1.1
In this simulation, #60 and #30 are used as the required measuring points for g1.1 and h1.1, respectively, and the other measuring points are at the locations of the two bearings and the disc.
Table A1. The maximum identification errors of the disc calculated by SDA and TODA.
Table A1. The maximum identification errors of the disc calculated by SDA and TODA.
MethodsRotor Unbalanceg1.1h1.1
SDAAmplitude6.77977 × 10−91.23934 × 10−8
Angle1.6455 × 10−97.73824 × 10−9
TODAAmplitude4.21247 × 10−52.79355 × 10−8
Angle3.32871 × 10−62.49249 × 10−8
(2)
Simulation results for rotor h4.4
In this simulation, #10 is used as one of the required measuring points and the other measuring points are at the locations of the two bearings and the four discs.
Figure A1. The maximum identification error for each rotor unbalance in h4.4 based on SDA using #10 point as one of the required measuring points: (a) relative error for rotor unbalance amplitude; (b) absolute error for rotor unbalance angle.
Figure A1. The maximum identification error for each rotor unbalance in h4.4 based on SDA using #10 point as one of the required measuring points: (a) relative error for rotor unbalance amplitude; (b) absolute error for rotor unbalance angle.
Applsci 12 03865 g0a1
Figure A2. The maximum identification error for each rotor unbalance in h4.4 based on TODA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure A2. The maximum identification error for each rotor unbalance in h4.4 based on TODA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g0a2
(3)
Simulation results for rotor g1.4
In this simulation, #20 is used as one of the required measuring points and the other measuring points are #15, #90, #21, #41, #61, and #81, which are the locations of the two bearings and the four discs.
Figure A3. The maximum identification error for each rotor unbalance in g1.4 based on SDA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure A3. The maximum identification error for each rotor unbalance in g1.4 based on SDA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g0a3
Figure A4. The maximum identification error for each rotor unbalance in g1.4 based on TODA using #20 point as one of the required measuring points: (a) relative error for rotor unbalance amplitude; (b) absolute error for rotor unbalance angle.
Figure A4. The maximum identification error for each rotor unbalance in g1.4 based on TODA using #20 point as one of the required measuring points: (a) relative error for rotor unbalance amplitude; (b) absolute error for rotor unbalance angle.
Applsci 12 03865 g0a4

Appendix B

(1)
Simulation results for g1.1 and h1.1
In this simulation, #60 and #30 are used as one of the required measuring points for g1.1 and h1.1, respectively, and the other measuring points are at the locations of the two bearings and the disc.
Table A2. The maximum identification error for the disc calculated by SDA and TODA.
Table A2. The maximum identification error for the disc calculated by SDA and TODA.
MethodsRotor UnbalanceG1.1H1.1
SDAAmplitude5.0000000075.000000013
Angle55.000000008
TODAAmplitude5.0000001085.000000022
Angle5.0000033295.000000025
(2)
Simulation results for rotor g1.4
In this simulation, #20 is used as one of the required measuring points and the other measuring points are #15, #90, #21, #41, #61, and #81, where the two bearings and the four discs are.
Figure A5. The maximum identification error for each rotor unbalance in g1.4 based on SDA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure A5. The maximum identification error for each rotor unbalance in g1.4 based on SDA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g0a5
Figure A6. The maximum identification error for each rotor unbalance in g1.4 based on TODA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure A6. The maximum identification error for each rotor unbalance in g1.4 based on TODA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g0a6
(3)
Simulation results for rotor h4.4
In this simulation, #10 is used as one of the required measuring points and the other measuring points are at the locations of the two bearings and the four discs.
Figure A7. The maximum identification error for each rotor unbalance in h4.4 based on SDA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure A7. The maximum identification error for each rotor unbalance in h4.4 based on SDA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g0a7
Figure A8. The maximum identification error for each rotor unbalance in h4.4 based on TODA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure A8. The maximum identification error for each rotor unbalance in h4.4 based on TODA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g0a8

Appendix C

(1)
Simulation results for g4.4
In this simulation, #20, #40, #60, and #80 points are used, respectively, and the other m + n measuring points are at the discs and bearings.
S D A _ M 1 g 4 4 = [ 7.74 7263.22 4675809.49 38337218.12 290.65 4.75 15819.92 131153.82 47436.53 2396.85 270.45 2287.39 2659853.86 134047.34 11889.74 4.97 ]  
S D A _ M 2 g 44 = [ 0.70   176.85 179.33 178.63 94.64 8.20 177.39 179.98 179.58 170.81 19.58 104.65 171.64 177.59 138.72 1.63 ]
T O D A _ M 1 g 44 = [ 2.72 × 10 3 2.20 × 10 5 6.15 × 10 8 5.97 × 10 9 4.66 × 10 4 5.71 × 10 3 7.17 × 10 6 6.97 × 10 7 7.48 × 10 6 3.64 × 10 4 4.97 × 10 4 1.78 × 10 5 3.49 × 10 9 1.60 × 10 7 6.20 × 10 6 1.27 × 10 4 ]
T O D A _ M 2 g 44 = [ 172.85   177.23 173.68 175.19 179.57 115.55 179.88 179.76 179.98 179.66 178.09 179.98 179.92 179.92 179.91 177.42 ]
where S D A _ M 1 g 44 and S D A _ M 2 g 44 are the matrix of the maximum identification error for rotor unbalance amplitude and the angle of #1 to #4 discs under different adjustment point conditions using SDA; T O D A _ M 1 g 44 and T O D A _ M 2 g 44 are the matrix of the maximum identification error for rotor unbalance amplitude and the angle of #1 to #4 discs under different adjustment point conditions using TODA.
Figure A9. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of g4.4 using #20 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A9. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of g4.4 using #20 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a9
Figure A10. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of g4.4 using #40 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A10. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of g4.4 using #40 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a10
Figure A11. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of g4.4 using #60 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A11. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of g4.4 using #60 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a11
Figure A12. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of g4.4 using #80 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A12. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of g4.4 using #80 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a12
(2)
Simulation results for h1.4
In this simulation, #10, #20, #30, and #40 points are used, respectively, and the other m + n measuring points are at the discs and bearings.
S D A _ M 1 h 14 = [ 62.06 377.46 15340.98 8311.13 495.48 9.29 5473.69 3128.49 1050.54 302.84 146.67 515.27 1741.08 422.80 650.49 42.56 ]    
S D A _ M 2 h 14 = [ 21.67   176.86 178.86 179.80 119.47 7.69 172.06 179.22 142.30 139.97 136.16 163.97 140.05 83.30 170.82 12.68 ]
T O D A _ M 1 h 14 = [ 678.59 11343.28 126989.97 68381.74 9068.03 611.79 96905.97 20941.13 41233.01 8545.67 4524.67 5890.04 91145.20 24188.79 18490.77 294.55 ]  
T O D A _ M 2 h 14 = [ 151.57   171.11 179.62 179.20 175.38 134.00 179.86 177.95 176.74 165.76 173.57 179.56 179.20 179.39 179.86 107.60 ]
where S D A _ M 1 h 14 and S D A _ M 2 h 14 are the matrix of the maximum identification error for rotor unbalance amplitude and the angle of #1 to #4 discs under different adjustment point conditions using SDA; T O D A _ M 1 h 14 and T O D A _ M 2 h 14 are the matrix of the maximum identification error for rotor unbalance amplitude and the angle of #1 to #4 discs under different adjustment point conditions using TODA.
Figure A13. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h1.4 using #10 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A13. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h1.4 using #10 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a13
Figure A14. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h1.4 using #20 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A14. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h1.4 using #20 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a14
Figure A15. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h1.4 using #30 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A15. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h1.4 using #30 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a15
Figure A16. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h1.4 using #40 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A16. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h1.4 using #40 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a16
(3)
Simulation results for h4.4
In this simulation, #10, #20, #30, and #40 points are used, respectively, and the other m + n measuring points are at the discs and bearings.
S D A _ M 1 h 44 = [ 0.10 40.21 648.96 330916.69 32.56 1.24 5473.69 5238.11 438.92 44.95 4.46 61.05 132071.11 3609.63 402.43 1.49 ]
S D A _ M 2 h 44 = [ 0.03 45.47 179.93 178.81 25.58 0.22 114.52 179.55 179.59 12.10 6.72 47.90 179.90 178.59 71.69 0.15 ]
T O D A _ M 1 h 44 = [ 0.73 764.47 28245.31 862557.51 3091.55 10.72 469.71 34447.58 43906.69 387.83 16.43 624.62 9631294.06 30595.13 267.13 6.50 ]
T O D A _ M 2 h 44 = [ 0.56 179.40 179.89 179.22 172.28 6.54 179.54 179.97 178.12 143.85 5.10 166.32 179.77 174.86 167.72 5.17 ]
where S D A _ M 1 h 44 and S D A _ M 2 h 44 are the matrix of the maximum identification error for rotor unbalance amplitude and the angle of #1 to #4 discs under different adjustment point conditions using SDA; T O D A _ M 1 h 44 and T O D A _ M 2 h 44 are the matrix of the maximum identification error for rotor unbalance amplitude and the angle of #1 to #4 discs under different adjustment point conditions using TODA.
Figure A17. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h4.4 using #10 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A17. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h4.4 using #10 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a17
Figure A18. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h4.4 using #20 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A18. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h4.4 using #20 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a18
Figure A19. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h4.4 using #40 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A19. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h4.4 using #40 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a19
Figure A20. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h4.4 using #40 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Figure A20. Statistical results of the number of AFPs of rotor unbalance of each disc in the simulation of h4.4 using #40 point as one of the required measuring points and a 0.1 nm resolution: (a) results of unbalance amplitude; (b) results of unbalance angle.
Applsci 12 03865 g0a20

References

  1. Bishop, R.E.D. The vibration of rotating shafts. J. Mech. Eng. Sci. 1959, 1, 50–65. [Google Scholar] [CrossRef]
  2. Gladwell, G.M.L.; Bishop, R.E.D. The vibration of rotating shafts supported in flexible bearings. J. Mech. Eng. Sci. 1959, 1, 195–206. [Google Scholar] [CrossRef]
  3. Goodman, T.P. A least-squares method for computing balance corrections. Trans. ASME J. Eng. Ind. 1964, 864, 273–279. [Google Scholar] [CrossRef]
  4. Bishop, R.E.D.; Gladwell, G.M.L. The vibration and balancing of an unbalanced flexible rotor. J. Mech. Eng. Sci. 1959, 1, 66–67. [Google Scholar] [CrossRef]
  5. Bishop, R.E.D.; Parkinson, A.G. On the isolation of modes in the balancing of flexible shafts. Proc. Inst. Mech. Eng. 1963, 177, 407–426. [Google Scholar]
  6. Lindley, A.L.G.; Bishop, R.E.D. Some recent research of the balancing of large flexible rotors. Proc. Inst. Mech. Eng. 1963, 177, 881–897. [Google Scholar] [CrossRef]
  7. Bishop, R.E.D.; Parkinson, A.G. On the use of balancing machines for flexible rotors. ASME J. Eng. Ind. 1972, 94, 561–576. [Google Scholar] [CrossRef]
  8. Kellenberger, W. Balancing flexible rotors on two generally flexible bearings. Brown Boveri Rev. 1970, 54, 603–619. [Google Scholar]
  9. Kellenberger, W. Should a flexible rotor be balanced in N or N + 2 planes? Trans. ASME J. Eng. Ind. 1972, 94, 548–560. [Google Scholar] [CrossRef]
  10. Lund, J.W.; Tonnesen, J. Analysis and experiments in multi-plane balancing of flexible rotors. Trans. ASME J. Eng. Ind. 1972, 194, 70–77. [Google Scholar] [CrossRef]
  11. Lund, J.W.; Tonnesen, J. Experimental and Analytic Investigation of High-Speed Rotor Balancing-Phase I; Research Report No. FR-8, Project No. FP-4; Department of Mechine Design, Technical University of Denmark: Copenhagen, Denmark, 1970. [Google Scholar]
  12. Shamsah, S.M.I.; Sinha, J.K.; Mandal, P. Estimating rotor unbalance from a single run-up and using reduced sensors. Measurement 2019, 136, 11–24. [Google Scholar] [CrossRef]
  13. Liu, Y.; Zhang, M.; Sun, C.; Hu, M.; Chen, D.; Liu, Z.; Tan, J. A method to minimize stage-by-stage initial unbalance in the aero engine assembly of multistage rotors. Aerosp. Sci. Technol. 2019, 85, 270–276. [Google Scholar] [CrossRef]
  14. Bin, G.; Li, X.; Shen, Y.; Wang, W. Development of whole-machine high speed balance approach for turbomachinery shaft system with N + 1 supports. Measurement 2018, 122, 368–379. [Google Scholar] [CrossRef]
  15. Xia, Y.; Ren, X.; Qin, W.; Yang, Y.; Lu, K.; Fu, C. Investigation on the transient response of a speed-varying rotor with sudden unbalance and its application in the unbalance identification. J. Low Freq. Noise Vib. Act. Control 2019, 39, 1065–1086. [Google Scholar] [CrossRef] [Green Version]
  16. Bently, D.E.; Muszynska, A. Modal testing and parameter identification of rotating shaft fluid lubricated bearing system. In Proceedings of the 4th International Conference on Modal Analysis, Orlando, FL, USA, April 1986; pp. 1393–1402. Available online: https://xueshu.baidu.com/usercenter/paper/show?paperid=635ff95826102054142b1957b8cde2db&site=xueshu_se&hitarticle=1 (accessed on 16 March 2022).
  17. Iida, H. Application of the Experimental Determination of Character Matrics to the Balancing of a Flexible Rotor. In Proceedings of the IFToMM 4th International Conference on Rotor Dynamics, Chicago, IL, USA, September 1994; pp. 111–115. Available online: https://xueshu.baidu.com/usercenter/paper/show?paperid=f42b876c8d2000aa6dd8baf7b534a8e8&site=xueshu_se&hitarticle=1 (accessed on 16 March 2022).
  18. Lou, X.; Zheng, S.; Wang, X. Some factors of leading to nonlinear electromagnetic force in active magnetic actuator. J. Mech. Electr. Eng. 1999, 2, 40–42. [Google Scholar]
  19. Liu, S.; Zheng, S. On-line unbalance identification of flexible rotor system. J. Zhejiang Univ. 2004, 38, 257–261. [Google Scholar]
  20. Shrivastava, A.; Mohanty, A.R. Estimation of single plane unbalance parameters of a rotor-bearing system using Kalman filtering based force estimation technique. J. Sound Vib. 2018, 418, 184–199. [Google Scholar] [CrossRef]
  21. Gillijns, S.; De Moor, B. Technical communique: Unbiased minimum-variance input and state estimation for linear discrete-time systems with direct feedthrough. Automatica 2007, 43, 934–937. [Google Scholar] [CrossRef]
  22. Shrivastava, A.; Mohanty, A.R. Identification of unbalance in a rotor system using a joint input-state estimation technique. J. Sound Vib. 2019, 442, 414–427. [Google Scholar] [CrossRef]
  23. Yao, J.; Liu, L.; Yang, F.; Scarpa, F.; Gao, J. Identification and optimization of unbalance parameters in rotor-bearing systems. J. Sound Vib. 2018, 431, 54–69. [Google Scholar] [CrossRef] [Green Version]
  24. Zou, D.; Zhao, H.; Liu, G.; Ta, N.; Rao, Z. Application of augmented Kalman filter to identify unbalance load of rotor-bearing system: Theory and experiment. J. Sound Vib. 2019, 463, 114972. [Google Scholar] [CrossRef]
  25. He, C.H.; Amer, T.S.; Tian, D.; Abolila, A.F.; Galal, A.A. Controlling the kinematics of a spring-pendulum system using an energy harvesting device. J. Low Freq. Noise Vib. Act. Control 2022. [Google Scholar] [CrossRef]
  26. He, J.H.; Amer, T.S.; Abolila, A.F.; Galal, A.A. Stability of three degrees-of-freedom auto-parametric system. Alex. Eng. J. 2022, 61, 8393–8415. [Google Scholar] [CrossRef]
  27. Wang, A.; Cheng, X.; Meng, G.; Xia, Y.; Wo, L.; Wang, Z. Dynamic analysis and numerical experiments for balancing of the continuous single-disc and single-span rotor-bearing system. Mech. Syst. Signal Process. 2017, 86, 151–176. [Google Scholar] [CrossRef]
  28. Tiwari, R.; Chakravarthy, V. Simultaneous estimation of the residual unbalance and bearing dynamic parameters from the experimental data in a rotor-bearing system. Mech. Mach. Theory 2009, 44, 792–812. [Google Scholar] [CrossRef]
  29. Tiwari, R. Conditioning of regression matrices for simultaneous estimation of the residual unbalance and bearing dynamic parameters. Mech. Syst. Signal Process. 2005, 19, 1082–1095. [Google Scholar] [CrossRef]
  30. Wang, A.; Yao, W.; He, K.; Meng, G.; Cheng, X.; Yang, J. Analytical modelling and numerical experiment for simultaneous identification of unbalance and rolling-bearing coefficients of the continuous single-disc and single-span rotor-bearing system with Rayleigh beam model. Mech. Syst. Signal Process. 2019, 116, 322–346. [Google Scholar] [CrossRef]
  31. Jamadar, I.M. A Model to Estimate Synchronous Vibration Amplitude for Detection of Unbalance in Rotor-Bearing System. Int. J. Acoust. Vib. 2021, 26, 161–169. [Google Scholar] [CrossRef]
  32. Ambur, R.; Rinderknecht, S. Unbalance detection in rotor systems with active bearings using self-sensing piezoelectric actuators. Mech. Syst. Signal Process. 2018, 102, 72–86. [Google Scholar] [CrossRef]
  33. Sanches, F.D.; Pederiva, R. Simultaneous identification of unbalance and shaft bow in a two-disk rotor based on correlation analysis and the SEREP model order reduction method. J. Sound Vib. 2018, 433, 230–247. [Google Scholar] [CrossRef]
  34. Zhang, Y.; Li, M.; Yao, H.; Gou, Y.; Wang, X. A modal-based balancing method for a high-speed rotor without trial weights. Mech. Sci. 2021, 12, 85–96. [Google Scholar] [CrossRef]
  35. Tiwari, R.; Chakravarthy, V. Simultaneous identification of residual unbalances and bearing dynamic parameters from impulse responses of rotor–bearing systems. Mech. Syst. Signal Process. 2006, 20, 1590–1614. [Google Scholar] [CrossRef]
  36. Tiwari, R.; Chakravarthy, V. Identification of the bearing and unbalance parameters from rundown data of rotors. In IUTAM Symposium on Emerging Trends in Rotor Dynamics; Springer: Dordrecht, The Netherlands, 2011; pp. 479–489. [Google Scholar]
Figure 1. A multi-disc and multi-span rotor-bearing system.
Figure 1. A multi-disc and multi-span rotor-bearing system.
Applsci 12 03865 g001
Figure 2. Identification procedures of SDA and TODA.
Figure 2. Identification procedures of SDA and TODA.
Applsci 12 03865 g002
Figure 3. Rotors supported by rolling bearings: (a) single-span single-disc rotor (g1.1); (b) single-span four-disc rotor (g1.4); (c) four-span four-disc rotor (g4.4).
Figure 3. Rotors supported by rolling bearings: (a) single-span single-disc rotor (g1.1); (b) single-span four-disc rotor (g1.4); (c) four-span four-disc rotor (g4.4).
Applsci 12 03865 g003
Figure 4. Rotors supported by oil journal bearings: (a) single-span single-disc rotor (h1.1); (b) single-span four-disc rotor (h1.4); (c) four-span four-disc rotor (h4.4).
Figure 4. Rotors supported by oil journal bearings: (a) single-span single-disc rotor (h1.1); (b) single-span four-disc rotor (h1.4); (c) four-span four-disc rotor (h4.4).
Applsci 12 03865 g004
Figure 5. The maximum identification error for each rotor unbalance in h1.4 based on SDA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure 5. The maximum identification error for each rotor unbalance in h1.4 based on SDA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g005
Figure 6. The maximum identification error for each rotor unbalance in h1.4 based on TODA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure 6. The maximum identification error for each rotor unbalance in h1.4 based on TODA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g006
Figure 7. The maximum identification error for each rotor unbalance in g4.4 based on SDA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure 7. The maximum identification error for each rotor unbalance in g4.4 based on SDA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g007
Figure 8. The maximum identification error for each rotor unbalance in g4.4 based on TODA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure 8. The maximum identification error for each rotor unbalance in g4.4 based on TODA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g008
Figure 9. Identification error changed with rotating frequency based on TODA using #20 point as one of the required measuring points in the simulation of g.4.4: (a) rotor unbalance amplitude of #3 disc; (b) rotor unbalance angle of #3 disc; (c) rotor unbalance amplitude of #4 disc; (d) rotor unbalance angle of #4 disc.
Figure 9. Identification error changed with rotating frequency based on TODA using #20 point as one of the required measuring points in the simulation of g.4.4: (a) rotor unbalance amplitude of #3 disc; (b) rotor unbalance angle of #3 disc; (c) rotor unbalance amplitude of #4 disc; (d) rotor unbalance angle of #4 disc.
Applsci 12 03865 g009
Figure 10. The maximum identification errors of each rotor unbalance in h1.4 based on SDA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure 10. The maximum identification errors of each rotor unbalance in h1.4 based on SDA using #10 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g010
Figure 11. The maximum identification error for each rotor unbalance in h1.4 based on TODA using #10 point as one of the required measuring points (the other measuring points are at the location of the two bearings and the four discs): (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure 11. The maximum identification error for each rotor unbalance in h1.4 based on TODA using #10 point as one of the required measuring points (the other measuring points are at the location of the two bearings and the four discs): (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g011
Figure 12. The maximum identification error for each rotor unbalance in g4.4 based on SDA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure 12. The maximum identification error for each rotor unbalance in g4.4 based on SDA using #20 point as one of the required measuring points: (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g012
Figure 13. The maximum identification error for each rotor unbalance in g4.4 based on TODA using #20 point as one of the required measuring points (the other measuring points are at the location of the two bearings and the four discs): (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Figure 13. The maximum identification error for each rotor unbalance in g4.4 based on TODA using #20 point as one of the required measuring points (the other measuring points are at the location of the two bearings and the four discs): (a) relative error for the rotor unbalance amplitude; (b) absolute error for the rotor unbalance angle.
Applsci 12 03865 g013
Figure 14. Identification error changed with rotating frequency based on TODA using #20 point as one of the required measuring points in the simulation for g.4.4: (a) rotor unbalance amplitude of #3 disc; (b) rotor unbalance angle of #3 disc; (c) rotor unbalance amplitude of #4 disc; (d) rotor unbalance angle of #4 disc.
Figure 14. Identification error changed with rotating frequency based on TODA using #20 point as one of the required measuring points in the simulation for g.4.4: (a) rotor unbalance amplitude of #3 disc; (b) rotor unbalance angle of #3 disc; (c) rotor unbalance amplitude of #4 disc; (d) rotor unbalance angle of #4 disc.
Applsci 12 03865 g014
Figure 15. Statistical results of the number of AFPs of rotor unbalance for each disc in the simulation of g1.4 using #20 point as one of the required measuring points and using a 0.1 nm resolution: (a) results for the unbalance amplitude; (b) results for the unbalance angle.
Figure 15. Statistical results of the number of AFPs of rotor unbalance for each disc in the simulation of g1.4 using #20 point as one of the required measuring points and using a 0.1 nm resolution: (a) results for the unbalance amplitude; (b) results for the unbalance angle.
Applsci 12 03865 g015
Figure 16. Statistical results of the number of AFPs of rotor unbalance for each disc in the simulation of g1.4 using #40 point as one of the required measuring points and using a 0.1 nm resolution: (a) results for the unbalance amplitude; (b) results for the unbalance angle.
Figure 16. Statistical results of the number of AFPs of rotor unbalance for each disc in the simulation of g1.4 using #40 point as one of the required measuring points and using a 0.1 nm resolution: (a) results for the unbalance amplitude; (b) results for the unbalance angle.
Applsci 12 03865 g016
Figure 17. Statistical results of the number of AFPs of rotor unbalance for each disc in the simulation of g1.4 using #60 point as one of the required measuring points and using a 0.1 nm resolution: (a) results for the unbalance amplitude; (b) results for the rotor unbalance angle.
Figure 17. Statistical results of the number of AFPs of rotor unbalance for each disc in the simulation of g1.4 using #60 point as one of the required measuring points and using a 0.1 nm resolution: (a) results for the unbalance amplitude; (b) results for the rotor unbalance angle.
Applsci 12 03865 g017
Figure 18. Statistical results of the number of AFPs of rotor unbalance for each disc in the simulation of g1.4 using #80 point as one of the required measuring points and using a 0.1 nm resolution: (a) results for the unbalance amplitude; (b) results for the unbalance angle.
Figure 18. Statistical results of the number of AFPs of rotor unbalance for each disc in the simulation of g1.4 using #80 point as one of the required measuring points and using a 0.1 nm resolution: (a) results for the unbalance amplitude; (b) results for the unbalance angle.
Applsci 12 03865 g018
Figure 19. Identification error changed with frequency based on TODA for the first kind of g.4.4 simulation in which the proposed adjustment points (#60 and #80 points) are applied: (a) rotor unbalance amplitude of #3 disc; (b) rotor unbalance angle of #3 disc; (c) rotor unbalance amplitude of #4 disc; (d) rotor unbalance angle of #4 disc.
Figure 19. Identification error changed with frequency based on TODA for the first kind of g.4.4 simulation in which the proposed adjustment points (#60 and #80 points) are applied: (a) rotor unbalance amplitude of #3 disc; (b) rotor unbalance angle of #3 disc; (c) rotor unbalance amplitude of #4 disc; (d) rotor unbalance angle of #4 disc.
Applsci 12 03865 g019
Figure 20. Identification error changed with frequency based on TODA for the second kind of g.4.4 simulation in which the proposed adjustment points (#60 and #80 points) are applied: (a) rotor unbalance amplitude of #3 disc; (b) rotor unbalance angle of #3 disc; (c) rotor unbalance amplitude of #4 disc; (d) rotor unbalance angle of #4 disc.
Figure 20. Identification error changed with frequency based on TODA for the second kind of g.4.4 simulation in which the proposed adjustment points (#60 and #80 points) are applied: (a) rotor unbalance amplitude of #3 disc; (b) rotor unbalance angle of #3 disc; (c) rotor unbalance amplitude of #4 disc; (d) rotor unbalance angle of #4 disc.
Applsci 12 03865 g020
Figure 21. Identified rotor unbalance changing with frequency in g1.4 using a 0.1 nm resolution: (a) obtained rotor unbalance amplitude; (b) obtained rotor unbalance angle.
Figure 21. Identified rotor unbalance changing with frequency in g1.4 using a 0.1 nm resolution: (a) obtained rotor unbalance amplitude; (b) obtained rotor unbalance angle.
Applsci 12 03865 g021
Figure 22. Identified rotor unbalance changing with frequency in g1.4 using a 1 nm resolution: (a) obtained rotor unbalance amplitude; (b) obtained rotor unbalance angle.
Figure 22. Identified rotor unbalance changing with frequency in g1.4 using a 1 nm resolution: (a) obtained rotor unbalance amplitude; (b) obtained rotor unbalance angle.
Applsci 12 03865 g022
Figure 23. Identified rotor unbalance changing with frequency in g1.4 using a 0.1 um resolution: (a) obtained rotor unbalance amplitude; (b) obtained rotor unbalance angle.
Figure 23. Identified rotor unbalance changing with frequency in g1.4 using a 0.1 um resolution: (a) obtained rotor unbalance amplitude; (b) obtained rotor unbalance angle.
Applsci 12 03865 g023
Figure 24. Identified rotor unbalance with frequency in g1.4 using a 1 um resolution: (a) obtained rotor unbalance amplitude; (b) obtained rotor unbalance angle.
Figure 24. Identified rotor unbalance with frequency in g1.4 using a 1 um resolution: (a) obtained rotor unbalance amplitude; (b) obtained rotor unbalance angle.
Applsci 12 03865 g024
Table 1. Meanings of symbols in the computational example of the rotor.
Table 1. Meanings of symbols in the computational example of the rotor.
ParameterMeaning
r_shaftRadius of the shaft
p_shaftDensity of the shaft
E_shaftElastic modulus of the shaft
L_shaftLength of the shaft
r_discRadius of the disc
p_discDensity of the disc
E_discElastic modulus of the disc
L_discWidth of the disc
Table 2. Parameters of the rotor shafts in the rotors.
Table 2. Parameters of the rotor shafts in the rotors.
ParameterValueParameterValue
r_shaft of rotor g1.1, g1.4 and g4.410 × 10−3 mL_shaft of the rotor g1.41200 × 10−3 m
r_shaft of rotor h1.1, h1.4 and h4.415 × 10−3 mL_shaft of the rotor g4.41600 × 10−3 m
p_shaft7800 kgm−3L_shaft of the rotor h1.11400 × 10−3 m
E_shaft2.1 × 1011 PaL_shaft of the rotor h1.41400 × 10−3 m
L_shaft of the rotor g1.1800 × 10−3 mL_shaft of the rotor h4.43600 × 10−3 m
Table 3. Parameters of the eccentric discs in the rotors g1.1, g1.4, and g4.4.
Table 3. Parameters of the eccentric discs in the rotors g1.1, g1.4, and g4.4.
ParameterValueParameterValueParameterValue
m 1 u 0.12031 kg m 2 u 0.15 kg m 3 u 0.01 kg
e 1 50 × 10−3 m e 2 10 × 10−3 m e 3 20 × 10−3 m
α 1 225° α 2 120° α 3 −120°
m 4 u 0.10 kgr_disc of #1 disc60 × 10−3 mp_disc of of #1~#4 disc7800 kgm−3
e 4 15 × 10−3 mr_disc of #2~#450 × 10−3 mE_disc of of #1~#4 disc2.1 × 1011 Pa
α 4 −170°L_disc of of #1~#4 disc10 × 10−3 m
Table 4. Parameters of the eccentric discs in the rotors h1.1, h1.4, and h4.4.
Table 4. Parameters of the eccentric discs in the rotors h1.1, h1.4, and h4.4.
ParameterValueParameterValueParameterValue
m 1 u 0.05 kg m 2 u 0.15 kg m 3 u 0.01 kg
e 1 30 × 10−3 m e 2 10 ×10−3 m e 3 20 × 10−3 m
α 1 45° α 2 90° α 3 170°
m 4 u 0.10 kgr_disc of #1–4 disc50 × 10−3 mp_disc of of #1~#4 disc7800 kgm−3
e 4 15 × 10−3 mL_disc of of #1~#4 disc10 × 10−3 m
α 4 −170°E_disc of of #1~#4 disc2.1 × 1011 Pa
Table 5. Parameters of the eight oil journal bearings in the rotors.
Table 5. Parameters of the eight oil journal bearings in the rotors.
ParameterValueParameterValue
[ k 1 x x k 1 x y k 1 y x k 1 y y ] [ 9 × 10 6 1 × 10 6 6 × 10 6 4 × 10 6 ] N/m [ k 2 x x k 2 x y k 2 y x k 2 y y ] [ 3 × 10 6 0.6 × 10 6 1.5 × 10 6 1 × 10 6 ] N/m
[ c 1 x x c 1 x y c 1 y x c 1 y y ] [ 8 × 10 5 1 × 10 5 1 × 10 5 10 × 10 5 ] N·s/m [ c 2 x x c 2 x y c 2 y x c 2 y y ] [ 1 × 10 5 1.4 × 10 5 1.4 × 10 5 1.5 × 10 5 ] N·s/m
[ k 3 x x k 3 x y k 3 y x k 3 y y ] [ 4 × 10 6 7 × 10 6 1.6 × 10 6 1.1 × 10 6 ] N/m [ k 4 x x k 4 x y k 4 y x k 4 y y ] [ 4.5 × 10 6 0.9 × 10 6 2 × 10 6 1.4 × 10 6 ] N/m
[ c 3 x x c 3 x y c 3 y x c 3 y y ] [ 2 × 10 5 1.5 × 10 5 1.5 × 10 5 1.6 × 10 5 ] N·s/m [ c 4 x x c 4 x y c 4 y x c 4 y y ] [ 1 × 10 5 2 × 10 5 2 × 10 5 2 × 10 5 ] N·s/m
[ k 5 x x k 5 x y k 5 y x k 5 y y ] [ 5 × 10 6 1.5 × 10 6 2.5 × 10 6 2 × 10 6 ] N/m [ k 6 x x k 6 x y k 6 y x k 6 y y ] [ 3 × 10 6 1 × 10 6 1 × 10 6 4 × 10 6 ] N/m
[ c 5 x x c 5 x y c 5 y x c 5 y y ] [ 2 × 10 5 3 × 10 5 3 × 10 5 3 × 10 5 ] N·s/m [ c 6 x x c 6 x y c 6 y x c 6 y y ] [ 1 × 10 5 0.6 × 10 5 0.6 × 10 5 2 × 10 5 ] N·s/m
[ k 7 x x k 7 x y k 7 y x k 7 y y ] [ 3 × 10 6 6 × 10 6 1.5 × 10 6 1 × 10 6 ] N/m [ k 8 x x k 8 x y k 8 y x k 8 y y ] [ 4 × 10 6 7 × 10 6 1.6 × 10 6 1.1 × 10 6 ] N/m
[ c 7 x x c 7 x y c 7 y x c 7 y y ] [ 4 × 10 5 1.4 × 10 5 1.4 × 10 5 1.5 × 10 5 ] N·s/m [ c 8 x x c 8 x y c 8 y x c 8 y y ] [ 7 × 10 5 1.5 × 10 5 1.5 × 10 5 16 × 10 5 ] N·s/m
Table 6. Parameters of the eight rolling bearings in the rotors g1.1, g1.4, and g4.4.
Table 6. Parameters of the eight rolling bearings in the rotors g1.1, g1.4, and g4.4.
ParameterValueParameterValue
[ k 1 x x k 1 x y k 1 y x k 1 y y ] [ 1.8 × 10 7 0 0 1.84 × 10 7 ] N/m [ k j j x x k j j x y k j j y x k j j y y ] ,   j j = 2 to 8 [ 9 × 10 6 0 0 8.9 × 10 6 ] N/m
[ c 1 x x c 1 x y c 1 y x c 1 y y ] , [ 75 0 0 70 ] N·s/m [ c j j x x c j j x y c j j y x c j j y y ] ,   j j = 2 to 8 [ 80 0 0 75 ] N·s/m
Table 7. Statistical results of peak values outside the assumed allowable range when the resolution is 0.1 nm.
Table 7. Statistical results of peak values outside the assumed allowable range when the resolution is 0.1 nm.
DiscUnbalanceValueValueValueValue
#1 Amplitude(10.84%, 1 Hz)---
Angle----
#2 Amplitude--(42.54%, 645 Hz)-
Angle----
#3 Amplitude(26.02%, 1 Hz)(257.24%, 269 Hz)-(20.64%, 1265 Hz)
Angle(19.49°, 1 Hz)-(11.72°, 645 Hz)-
#4Amplitude----
Angle-(16.259°, 269 Hz)--
Table 8. Statistical results of peak values outside the assumed allowable range when the resolution is 1 nm.
Table 8. Statistical results of peak values outside the assumed allowable range when the resolution is 1 nm.
DiscUnbalanceValueValueValueValue
#1 Amplitude(83.36%, 1 Hz)---
Angle----
#2 Amplitude(95.6%, 1 Hz)(42.54%, 269 Hz)--
Angle(32.5°, 1 Hz)---
#3 Amplitude(242.24%, 1 Hz)(257.24%, 269 Hz)(11.64%, 645 Hz)(20.64%, 1265 Hz)
Angle(108.42°, 1 Hz)(57.41°, 271 Hz)(11.72°, 645 Hz)-
#4Amplitude----
Angle(165.08°, 1 Hz)(16.259°, 269 Hz)--
Table 9. Statistical results of peak values outside the assumed allowable range when the resolution is 0.1 um.
Table 9. Statistical results of peak values outside the assumed allowable range when the resolution is 0.1 um.
DiscUnbalanceValueValueValueValue
#1 Amplitude(695.26%, 4 Hz)(19.45%, 268 Hz)--
Angle(140.04°, 4 Hz)(10.86°, 268 Hz)--
#2 Amplitude(597.53%, 4 Hz)(22.7%, 270 Hz)--
Angle(94.44°, 4 Hz)---
#3 Amplitude(1722.04%, 8 Hz)(262.49%, 268 Hz)-(11.76%, 1268 Hz)
Angle(120.31°, 16 Hz)(20.6°, 270 Hz)-(11.3°, 1268 Hz)
#4Amplitude(853.26%, 6 Hz)---
Angle(165.1°, 10 Hz)(10.04°, 270 Hz)--
Table 10. Statistical results of peak values outside the assumed allowable range when the resolution is 1 um.
Table 10. Statistical results of peak values outside the assumed allowable range when the resolution is 1 um.
DiscUnbalanceValue
#1 discAmplitude-
Angle-
#2 discAmplitude-
Angle-
#3 discAmplitude(11.73%, 1264 Hz)
Angle(11.56°, 1268 Hz)
#4discAmplitude-
Angle-
Table 11. First three critical frequencies of g1.4 from 1 Hz to 2000 Hz.
Table 11. First three critical frequencies of g1.4 from 1 Hz to 2000 Hz.
First OrderSecond OrderThird Order
269 Hz645 Hz1267 Hz
Table 12. Statistical results of peak values of identification errors of rotor unbalance for #1 disc in g1.4.
Table 12. Statistical results of peak values of identification errors of rotor unbalance for #1 disc in g1.4.
Sensor ResolutionLow Frequency (1 Hz)First Order (269 Hz)Second Order (645 Hz)Third Order (1267 Hz)
0.1 nm(10.84%, 1 Hz)(2.02%, 271 Hz)(0.31%, 643 Hz)(0.42%, 1267 Hz)
(2.12°, 1 Hz)(1.47°, 271 Hz)(0.32°, 645 Hz)(0.32°, 1267 Hz)
1 nm(83.36%, 1 Hz)(1.92%, 269 Hz)(0.31%, 643 Hz)(0.42%, 1267 Hz)
(8.77°, 1 Hz)(1.45°, 271 Hz)(0.32°, 645 Hz)(0.32°, 1267 Hz)
0.1 um(695.26%, 4 Hz)(19.45%, 268 Hz)(0.65%, 646 Hz)(0.62%, 1266 Hz)
(140.04°, 4 Hz)(10.86°, 268 Hz)(0.55°, 646 Hz)(0.45°, 1262 Hz)
1 umCannot be identifiedCannot be identifiedCannot be identified(0.63%, 1266 Hz)
Cannot be identifiedCannot be identifiedCannot be identified(0.45°, 1262 Hz)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, A.; Bi, Y.; Feng, Y.; Yang, J.; Cheng, X.; Meng, G. Continuous Rotor Dynamics of Multi-Disc and Multi-Span Rotors: A Theoretical and Numerical Investigation of the Identification of Rotor Unbalance from Unbalance Responses. Appl. Sci. 2022, 12, 3865. https://doi.org/10.3390/app12083865

AMA Style

Wang A, Bi Y, Feng Y, Yang J, Cheng X, Meng G. Continuous Rotor Dynamics of Multi-Disc and Multi-Span Rotors: A Theoretical and Numerical Investigation of the Identification of Rotor Unbalance from Unbalance Responses. Applied Sciences. 2022; 12(8):3865. https://doi.org/10.3390/app12083865

Chicago/Turabian Style

Wang, Aiming, Yujie Bi, Yu Feng, Jie Yang, Xiaohan Cheng, and Guoying Meng. 2022. "Continuous Rotor Dynamics of Multi-Disc and Multi-Span Rotors: A Theoretical and Numerical Investigation of the Identification of Rotor Unbalance from Unbalance Responses" Applied Sciences 12, no. 8: 3865. https://doi.org/10.3390/app12083865

APA Style

Wang, A., Bi, Y., Feng, Y., Yang, J., Cheng, X., & Meng, G. (2022). Continuous Rotor Dynamics of Multi-Disc and Multi-Span Rotors: A Theoretical and Numerical Investigation of the Identification of Rotor Unbalance from Unbalance Responses. Applied Sciences, 12(8), 3865. https://doi.org/10.3390/app12083865

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop