Next Article in Journal
A Separation Method of Superimposed Gratings in Double-Projector Fringe Projection Profilometry Using a Color Camera
Next Article in Special Issue
On the Development of a Surrogate Modelling Toolbox for Virtual Assembly
Previous Article in Journal
The Moderating Role of IT in Process of Shaping Organizational Performance by Dynamic Capabilities of Controlling
Previous Article in Special Issue
Analysis and Experimental Verification of Mechanical Errors in Nine-Link Type Double-Toggle Mold/Die Clamping Mechanisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Unbalance Optimization Method for a Multi-Stage Rotor Based on an Assembly Error Propagation Model

1
Centre of Ultra-Precision Optoelectronic Instrument Engineering, Harbin Institute of Technology, Harbin 150080, China
2
Key Lab of Ultra-Precision Intelligent Instrumentation, Harbin Institute of Technology, Ministry of Industry and Information Technology, Harbin 150080, China
*
Author to whom correspondence should be addressed.
Submission received: 1 December 2020 / Revised: 9 January 2021 / Accepted: 15 January 2021 / Published: 19 January 2021
(This article belongs to the Special Issue Mechanical Tolerance Analysis in the Era of Industry 4.0)

Abstract

:
For the assembly of a multi-stage rotor, such as an aero-engine or gas turbine, the parts need to be assembled optimally to avoid excessive unbalance. We propose a method to optimize the unbalance of a multi-stage rotor during assembly. First, we developed an assembly error propagation model for a multi-stage rotor. The alignment process and distribution of the screw holes of the adjacent rotors was considered for the first time. Secondly, we propose a new assembly datum for unbalance optimization to ensure consistency with the actual conditions of a dynamic balance test. Finally, the unbalance optimization of a multi-stage rotor was achieved using a genetic algorithm, and the corresponding optimal assembly orientations of rotors at different stages were also identified. The results of the simulations showed that the assembly error propagation model had high accuracy and that the genetic optimization process had good convergence. The effect of unbalance optimization was also proven with experiments.

1. Introduction

The core components of rotating machinery, such as the high-pressure compressor of an aero-engine, are typically assembled by several single-stage rotors [1,2,3]. The machining errors of a single-stage rotor will propagate and accumulate during the assembly, which makes the actual assembly location of the rotor deviate from its ideal design location and the center of mass of the rotor deviate from the central inertia principal axis of the assembly, resulting in unbalance [4,5]. The unbalance caused by improper assembly is the main reason for vibration faults of rotating machinery.
The traditional way to eliminate unbalance is to measure the unbalance of the rotor with a dynamic balancing machine and then grind the rotor on a pre-set balanced surface [6,7]. This often requires repeated test runs and repairs to ensure the rotor meets the vibration requirements. This method not only costs time but also changes the initial unbalance of the single-stage rotor with the balanced surfaces. If the multi-stage rotor is reassembled, the unbalance of the single-stage rotor with balanced surfaces should be measured again. Therefore, there is an urgent need for an assembly optimization method that can not only predict the cumulative error of a multi-stage rotor after assembly but also achieve optimal matching of the assembly orientations of rotors at different stages to improve the assembly efficiency and one-time assembly qualified rate of multi-stage rotors.
The assembly optimization method of multi-stage rotors has been derived in many studies seeking to minimize the assembly cumulative errors. Hussain et al. [8] proposed an assembly error propagation model for 2-D rotational parts. For the first time, they came up with a strategy of straight-build assembly. To be specific, the line of centers of rotors at diverse stages can be as straight as possible by altering the assembly orientations of these rotors. In [9], five assembly optimization methods were evaluated to reduce the error build-up in a straight-build assembly of 2D rectangular components.
Yang et al. [10] proposed another strategy of a parallelism-build assembly based on the above 2-D error propagation model. In line with such a strategy, the angular error between the actual and nominal mounting surfaces was minimized for the last-stage rotor. Such a 2-D error propagation model was further modified into a 3-D model by Yang et al. [11]. Then, they assumed that the key dimensions of rotors at various stages obeyed a normal distribution within a specified tolerance zone, and the eccentricity error of assembled rotors were predicted by tolerance analysis.
Hussain et al. [12] derived a progressive eccentricity error optimization model, and the corresponding physical significance still conformed to the concept of a straight-build assembly. Yang et al. [13] also probed into how the number of the circumferential assembly orientations affects the cumulative error of assembly. By constructing a probability density function of cumulative eccentric error of the last-stage rotor, Yang et al. [14] compared the global optimization results obtained by probability and the Monte Carlo method.
Jin et al. [15] proposed a solution to partial parallel chains for error analysis in aero-engine assembly. Sun et al. [16] proposed a neural network model for predicting the concentricity and perpendicularity of a multi-stage rotor. Depending on the location and orientation tolerance of rotor parts, an assembly error propagation model of a multi-stage rotor was constructed by Wang et al. [17] to achieve the optimal assembly orientations of rotors at different stages. Sun et al. [18] also derived the same error propagation model according to the coordinate propagation principle, which can provide guidance for the maintenance of aero-engines.
In all the above studies, the optimization objectives and known conditions all focused on the geometric attributes of rotors. The optimization of a rotor’s mass attributes should be especially highlighted to prevent the vibrations of a multi-stage rotor from exceeding its upper limits. Piskin et al. [19] proposed a novel turbine balancing method using an Ant Colony algorithm. To optimize rotor unbalance, an assembly method was put forward by Liu et al. [20]. However, the center-of-mass coordinate in this study was not obtained by actual measurement but estimated by assuming a conversion ratio and based on the centroid coordinates. To overcome this defect, Sun et al. [21] made use of a vertical dynamic balancing machine to measure the mass attributes of rotors at different stages; on this basis, an assembly approach was proposed for the purpose of optimizing both the coaxiality and unbalance of rotors.
It can be seen from the current studies that an accurate error propagation model is very important for the prediction of assembly cumulative error and the guidance of the assembly process. The existing error propagation models were derived based on the axisymmetric characteristics of the rotating body and combined with coordinate propagation theories; however, the alignment processes and distribution of the screw holes around the rotors at different stages were ignored, and the selection of assembly orientations should be discrete.
The above problems are not clearly specified in the existing assembly optimization methods. When optimizing the rotor unbalance, the selection of the assembly datum should not be the same as that of the coaxiality optimization but should be consistent with the actual rotating axis of the dynamic balance test.
Targeted at the above problems, an unbalance optimization method of a multi-stage rotor is proposed in this study. In Section 2.1, an assembly error propagation model is developed, and the alignment processes and distribution of screw holes of adjacent rotors are introduced into the model. In Section 2.2, a new assembly datum for unbalance optimization of multi-stage rotor is calculated to ensure that the optimization datum and the measurement datum is truly realized. In Section 2.3, the unbalance optimization of a multi-stage rotor is achieved using a genetic algorithm (GA). In Section 3.1, the calculated accuracy of the assembly error propagation model is validated, and the optimization performance of GA is presented in Section 3.2. Finally, in Section 4, the effect of the unbalance optimization is proven in experiments.

2. Methods

2.1. Assembly Error Propagation Model

Limited to the current industrial level, the geometric deviations of the single-stage rotor are inevitable in the process of machining. The concentricity and parallelism errors of the assembly mounting surfaces of rotors at different stages are propagated stage by stage continuously, forming assembly cumulative errors. In Figure 1, the initial state of a 3-D two-stage simulation rotor before assembly is portrayed. Clearly, when the bottom mounting surface of each single-stage rotor is selected as the measurement datum, P1 refers to a parallelism error of the top mounting surface of Rotor 1, P2 to that of the top mounting surface of Rotor 2, and C1 and C2 are the concentricity errors of Rotor 1 and Rotor 2, respectively.
A pair of screw holes is randomly selected in the circumferential direction of the rotor for calibration, and the sampling angle of the highest point of the top mounting surface is recorded with the center of the calibrated screw hole as the sampling starting point. The red and blue dots in Figure 1 are the calibrated screw holes of Rotor 1 and Rotor 2, respectively. M1 and M2 are the highest points of the top mounting surface of Rotor 1 and Rotor 2, respectively. N1 and N2 are the lowest points of the top mounting surface of Rotor 1 and Rotor 2, respectively.
In Figure 1, the bottom mounting surface of Rotor 1 is used as an XY plane, and the normal axis passing through it is taken as the Z axis. The X axis is formed by a projection axis on the bottom surface of a line connecting the highest and the lowest points of the top mounting surface of Rotor 1. δ refers to the phase difference between the calibrated screw holes of Rotor 1 and Rotor 2. After that, the assembly process of the two-stage rotor can be divided into the following four steps (see Figure 2):
Step-1: Alignment of the calibrated screw holes.
In actual assembly, the screw holes of Rotor 1 and Rotor 2 must be aligned first (i.e., Rotor 2 rotates an angle of δ around the axis Z relative to Rotor 1), and then the assembly orientation can be selected with this state as the initial state. Although the angle of δ cannot be determined, the phase difference (φ2φ1) between the highest points of Rotor 1 and Rotor 2 after assembly can be obtained by measuring the sampling angle between the calibrated screw hole and its highest point of each stage rotor.
Step-2: Selection of the optimal assembly orientation.
Assuming that Rotor 1 is fixed, the calibrated screw holes of Rotor 2 and Rotor 1 must be aligned first. This is the initial state to select the appropriate assembly orientation of Rotor 2, and the rotation matrix of Rotor 2 relative to Rotor 1 about the Z axis can be expressed as follows:
R z 2 = [ cos [ θ z 2 + ( φ 2 φ 1 ) ] sin [ θ z 2 + ( φ 2 φ 1 ) ] 0 sin [ θ z 2 + ( φ 2 φ 1 ) ] cos [ θ z 2 + ( φ 2 φ 1 ) ] 0 0 0 1 ]
where θz2 refers to an angle by which Rotor 2 rotates around the Z axis relative to Rotor 1, φ2 refers to an angle between the center of the calibrated screw hole of Rotor 2 and its highest point, and φ1 to that of Rotor 1. (φ2φ1) refers to the phase difference between the highest point of Rotor 2 and that of Rotor 1 after alignment of the calibrated screw holes of Rotor 2 and Rotor 1. Considering the discrete and uniform distribution of the screw holes, the selection range of θz2 was 0−180°, and the optional angle was q(180/k), where k is the number of circumferential screw holes on the mounting surface, and q represents how many distribution angles the calibrated screw hole needs to rotate.
Step-3: Rotation.
Rotor 2 is rotated by an angle of β2 around the Y axis so that the bottom mounting surface of Rotor 2 is parallel to the top mounting surface of Rotor 1. Then the rotation matrix of Rotor 2 relative to Rotor 1 about the Y axis can be expressed as follows:
R y 2 = [ cos ( β 2 ) 0 sin ( β 2 ) 0 1 0 sin ( β 2 ) 0 cos ( β 2 ) ]
where β2 is an angle formed by the top and bottom mounting surfaces of Rotor 2 and can be calculated as follows:
β 2 = arctan ( P 2 2 D 2 )
where D2 is the measured radius of the top mounting surface of Rotor 2.
Step-4: Translation.
Rotor 2 is translated so that the center of the bottom mounting surface of Rotor 2 coincides with the top mounting surface of Rotor 1, which is equivalent to directly accumulating the center coordinates of the top mounting surface of Rotor 1 on the basis of Step-2.
Through the above steps, the assembly error propagation model of a 2-stage rotor can be expressed as follows:
X 2 = I 2 R z 2 R y 2 + S 1
where X2 refers to a coordinate vector of any point in Rotor 2 after assembly, and I2 to that in Rotor 2 before assembly. Rotor 1 does not move by default, and S1 is the coordinate vector of the center of the top mounting surface of Rotor 1.
The assembly process of an n-stage rotor can be regarded as the assembly processes of (n − 1) two-stage rotor. Therefore, the assembly error propagation model of an n-stage rotor can be expressed as follows:
X n = I n [ n : 1 : 1 ( R z n R y n ) ] + S 1 ( n N , n > 1 )
where Xn refers to a coordinate vector of any point in rotor n after assembly, In to that in rotor n before assembly, Rzn to a rotation matrix of rotor n relative to rotor (n − 1) about the Z axis, and Ryn to a rotation matrix of rotor n relative to rotor (n − 1) about the Y axis.
For a rigid rotor, the coordinate vector of any point in rotors at different stages before and after assembly can be accurately predicted by Equation (5) on the premise that the geometric errors of rotors at different stages are known. The parameters needed to solve Equation (5) include concentricity (Cn), parallelism (Pn), measured radius (Dn), the angle between the calibrated screw hole and the highest point (φn), and the height of mounting surface (Hn), which can be obtained using a coordinate measuring instrument (CMM); only θzn is unknown and is to be optimized.

2.2. Assembly Datum for Unbalance Optimization

The mathematical definition of unbalance is the product of the mass of a rotor and the eccentricity error of the center of mass and can also be expressed as the product of an unbalanced mass block and its effect radius. During the dynamic balance test of a rotor, journals of the front and back shafts of the rotor are respectively positioned on the support frames on left and right sides of a dynamic balancing machine.
The axis of rotation should be approximately the central line between the journals of the front and back shafts, and the measured vector radius of the mass eccentricity should be perpendicular to such an axis of rotation (i.e., the orientation-varying axis in Figure 3). This axis of rotation will vary along with variations in the assembly orientations of rotors at different stages. For the optimization of the rotor unbalance, this axis of rotation should be used as the assembly datum. The unbalance of rotors at different stages after assembly can be obtained as follows:
U n = m n r n ,
where Un is the unbalance of rotor n after assembly, that is a product of mn (unbalanced mass) and rn′ (i.e., the vector of the effect radius of rotor n). As for rn′, it can be solved by the following steps:
Step-1: A linear equation expressing the orientation-varying axis can be written as follows:
x x 0 e n x x 0 = y y 0 e n y y 0 = z z 0 e n z z 0 = t n ,
where (x0, y0, z0) is a coordinate vector of the bottom center of the bottom rotor. Under the circumstance that the table axis serves as the datum for measurement of the rotor’s geometric parameters, such a point is the origin of the coordinates by default.
Step-2: As for the plane Qn that is extended across an unbalanced mass point of rotor n and perpendicular to the orientation-varying axis, it is expressed in the equation below:
{ ( e n x x 0 ) ( x r n x ) + ( e n y y 0 ) ( y r n y ) + ( e n z z 0 ) ( z r n z ) = 0 ( R n = [ r n x r n y r n z ] , R n = [ r n x r n y r n z ] ) ,
where Rn′ refers to a coordinate vector of the unbalanced mass point of rotor n after assembly. Rn refers to the initial coordinate vector of the unbalanced mass point of rotor n before assembly, and it can be measured by the dynamic balancing machine or artificially added. By substituting Rn into Equation (5), Rn′ is solved.
Step-3: Coordinates of an intersecting point between the orientation-varying axis and the plane Qn are expressed as follows:
First, Equations (7) and (8) form a simultaneous equation. Through such a simultaneous equation, parameter tn in the linear equation expressing the orientation-varying axis can be obtained:
t n = ( e n x x 0 ) ( r n x x 0 ) + ( e n y y 0 ) ( r n y y 0 ) + ( e n z z 0 ) ( r n z z 0 ) ( e n x x 0 ) 2 + ( e n y y 0 ) 2 + ( e n z z 0 ) 2 .
By substituting tn into Equation (7), the coordinates of the intersection point Jn can be acquired.
Step-4: The effect radius of the unbalanced mass relative to the orientation-varying axis, that is the modulus of, rn′, is expressed in the following equation:
r n = ( j n x r n x ) 2 + ( j n y r n y ) 2 + ( j n z r n z ) 2 ( J n = [ j n x j n y j n z ] ) .
By substituting Equation (10) into Equation (6), a value of Un is acquired. After that, the unbalance of rotors at diverse stages should be decomposed to the balanced surfaces A and B through Equation (11):
{ U A = j = 1 n l B l j l B l A U j U B = j = 1 n l j l A l B l A U j ,
where UA and UB turn out to be the unbalance of an n-stage rotor on the balanced Surfaces A and B, respectively. lA and lB are the distances between the balanced Surfaces A and B from the support point of the journal, respectively. lj is the distance between the unbalanced measuring surface of each rotor and the balanced Surface A. In this context, the maximum unbalance of the n-stage rotor can be obtained by the following equation:
U = max ( U A , U B ) .

2.3. Genetic Optimization Algorithm

As a global optimization algorithm, genetic algorithms (GA) are aimed at mapping the solution space to a genetic space, encoding each set of possible solutions into a chromosome. First, a series of candidate solutions are randomly generated, and they form an initial population. Based on a fitness function designed in advance according to the objective function, the fitness of each individual to the problem solving environment is calculated. Then, corresponding selections are made in consistency with the obtained fitness, suppressing chromosomes with low fitness, but promoting those with high fitness.
After that, genetic manipulations, such as crossovers and mutations are performed for the purpose of evolving the next-generation population. Through repeated operations as described above, constant evolution toward optimized solutions is achieved. At last, except for a population that meets the corresponding convergence conditions and is the most adaptive to the problem solving environment, the optimal solution is also obtained. Depending on the above thoughts of genetic optimization, a GA fitness function of the unbalance of an n-stage rotor can be expressed as follows:
{ f i t n e s s _ U ( x ) = U ( x ) s . t .   x = ( θ z 1 , θ z 2 , , θ z n ) , 0 θ z n 180 ° ,
where design variable x represents a vector formed by assembly orientations of the n-stage rotor. In terms of GA, its steps are detailed below:
Step-1 (initial population): the n elements in an assembly orientation sequence x are regarded as n genes; and each gene varies from 0° to 180°. Each chromosome is composed of n genes. In this way, 1000 chromosomes are randomly generated and act as an initial population.
Step-2 (individual evaluation): the fitness of each chromosome is determined for individual evaluation.
Step-3 (selection): chromosomes with the optimal fitness are selected and take the place of a chromosome of the worst fitness.
Step-4 (mutation): the purpose of mutation is to improve the current chromosome. According to Equation (14), the chromosome with higher fitness has a smaller variation range. As for chromosomes with lower fitness, their variation ranges are larger accordingly.
x n e w = x [ 1 ± γ ( 1 f f b e s t ) 2 ] ,
where xnew and x are chromosomes after/before mutation respectively. γ is a random number between 0 and 1, and f stands for the fitness of x, while fbest stands for the optimal fitness of the present population.
Step-5 (crossover): some genes are randomly selected from two chromosomes for replacement and recombination. In this way, a new chromosome is generated. By means of crossover, the search capability of the GA can be enhanced.
Step-6 (termination criteria): the algorithm is ended when the fitness of optimal individuals and the population ceases to rise and the number of iterations reaches the designed value.

3. Simulation

3.1. Validation of Assembly Error Propagation Model

In Section 2.1, an assembly error propagation model of a multi-stage rotor was developed to predict the changes of coordinate values of any point in each rotor before and after assembly on the premise that the geometric parameters of rotors at different stages are known. To verify the calculated accuracy of the assembly error propagation model, the geometric parameters of a three-stage rotor presented in [17] (see Table 1) are used for reference for a simulation. In addition, the angle between the calibrated screw holes and the highest points is 0° by default (i.e., φ1 = φ2 = φ3 = 0°). The assembly orientations of the first rotor, second rotor, and third rotor were set to 0°, 30°, and 60° respectively (i.e., θz1 = 0°; θz2 = 30°; and θz3 = 60°).
The new coordinates of the center of the top mounting surface of the second and the third stage rotors after assembly were calculated using the model. SolidWorks, a 3-D software for drawing, was utilized to reconstruct the corresponding three-stage rotor model. Through a measurement function of the software, the coordinates of the center of the top mounting surface of the assembled rotors at different stages were measured (see Figure 4). As shown in Table 2, the calculated results achieved by the model are in good agreement with those measured by the software.

3.2. GA-Based Unbalance Optimization

We assumed that the three-stage simulation rotor described in Section 3.1 is made of steel. The initial coordinates of the center of mass of the rotors at different stages were measured using software (see Table 3). The initial unbalance of each rotor is the product of its mass and the vertical distance from the center of mass to the central inertia principal axis, which can be obtained using Equation (6). As shown in Figure 4, the balanced Surfaces A and B of the three-stage rotor are perpendicular to the central principal axis of inertia (i.e., the line connecting the center of the top mounting surface and the bottom surface) and, respectively, pass through the center of the top mounting surface and the bottom surface.
lA is set as 0, and lB is equal to the distance between the center of the top and bottom mounting surfaces of the three-stage rotor. The unbalance on the balanced Surfaces A and B can be obtained by Equation (11), and the integral unbalance of the three-stage rotor can be obtained by Equation (12). The minimum and maximum unbalance are taken as the optimization objectives, and the genetic algorithm proposed in Section 2.3 is applied to calculate the optimal assembly orientations of the second and third stage rotors. In terms of the first stage rotor, its assembly orientation was designed at 0° by default; each chromosome was designed to hold two genes, that is (θz2, θz3); and, the number of initial populations was set at 1000.
In addition, the selection principle of the internal parameters of the GA referred to the literature [22], and the mutation rate, crossover rate, and iterative times were set as 0.01, 0.9, and 500, respectively. The convergence progress of the GA-based unbalance optimization of the three-stage rotor is shown in Figure 5. When θz2 and θz3 were 0° and 180°, respectively, the corresponding minimum unbalance reached 32.2568 g·mm (minimum unbalance), as shown in Table 4. When θz2 and θz3 were 180° and 180°, respectively, the corresponding maximum unbalance reached 129.6123 g·mm (maximum unbalance).
Further considering the distribution of the assembly screw holes, the optional assembly orientations of the rotors at different stage are discrete, and the optimal corresponding relationship of the screw holes of the rotors should be obtained by a secondary calculation according to the above optimal solution. As shown in Figure 6, assuming the optimal orientation obtained by the optimization is θ, θ1 and θ2 are the two optional orientations adjacent the optimal orientation. The fitness in both cases must be calculated again, and the optimal solution can be selected. The angles of the screw holes adjacent to the optimal solution were determined by calculating all possible alignment relations (number: 2n−1) for n-stage rotors.
As assumed, there were 24 circumferential screw holes on the top and bottom mounting surfaces of each rotor of the above three-stage simulation rotor; and the distributed angle of the screw holes was set at 15°. Since the optimal solution is exactly a multiple of the distribution angle (15°), the minimum unbalance can be obtained when Rotor 2 rotates zero distributed angles (0 × 15 = 0°) relative to Rotor 1, and Rotor 3 rotates 12 distributed angles (12 × 15 = 180°) relative to Rotor 2. When Rotor 2 rotates 12 distributed angles (12 × 15 = 180°) relative to Rotor 1 and Rotor 3 rotates 12 distributed angles (12 × 15 = 180°) relative to Rotor 2, the unbalance of the assembly reaches the maximum value.

4. Experiments

An aero-engine was scaled in a certain proportion and simplified to a four-stage simulation rotor. A section view of the four-stage rotor is shown in Figure 7. As observed, it was composed of a front shaft, a low-pressure compressor (LPC), a high-pressure compressor (HPC), and a back shaft. Each single-stage rotor was assembled step by step through its own axial and radial front edge, and connected by a certain number of screws. The distances between the two pre-set balanced Surfaces A and B from the support point were lA and lB, respectively. To further prove the validity of the proposed assembly optimization method in this paper, experiments were performed, and the corresponding experimental steps are detailed below:
Step-1: The geometric parameters of the rotors at different stages were measured by a 3-D CMM.
Figure 8 shows the measured actual scene and schematic diagram of the HPC. First, the center of a pair of screw holes on the mounting surface of HPC was calibrated. Taking this center as the starting point of sampling, the concentricity (Cn), the parallelism (Pn), the measured radius (Dn), an angle between this center and the highest point (φn), and the high point (H) of the mounting surface of the HPC were measured successively. As shown in Figure 8b, the circle of diameter t1 is the tolerance zone defined by the concentricity specification, and the distance t2 is the tolerance zone defined by the parallelism specification. The measuring principle of the LPC and front and back shaft was the same as that of the HPC. The geometric parameters of component of the four-stage rotor were measured by a 3-D CMM for six times, and presented in Table 5, Table 6, Table 7 and Table 8, respectively. As shown in Table 9, the standard uncertainties of different measured parameters were calculated by the equation below:
u ( x ) = i = 1 n ( x i x a ) 2 n ( n 1 ) ,
where n refers to the number of measurements, and xa refers to the mean value of the six measurement results. The distributed numbers of the screw holes of rotors at different stages are shown in Table 10. The optimal assembly orientations of each rotor can only be selected from these discrete angles.
Step-2: The initial unbalance of the four-stage rotor was artificially added to each rotor.
To better reflect the experimental effect of the unbalance optimization, the unbalanced mass was artificially added to the rotors at different stages, which is far greater than the unbalance caused by uneven material distribution. Hexagon-socket-head-screws were used for the connection of the four-stage rotor, as shown in Figure 9. The method was to replace an original M4 × 12-screw with a M4 × 20-screw at the calibrated screw hole of rotors at different stages to produce an unbalanced mass with a mass difference of about 2 g, whose effect radius and phase are known. The initial unbalance parameters of the four-stage rotor are shown in Table 11.
Step-3: The unbalance of the four-stage rotor was measured using a dynamic balancing machine.
According to the geometric and unbalance parameters obtained in Step-1 and Step-2, the genetic algorithm derived in Section 3.2 was used to calculate the optimal assembly orientations of the rotors at different stages with three assembly strategies as follows:
Strategy-1: the minimum unbalance was set as the optimization objective;
Strategy-2: the maximum unbalance was set as the optimization objective; and
Strategy-3: direct assembly (i.e., the rotors at each stage remained in their initial state).
The optimal assembly orientations of each component and the corresponding unbalance of the four-stage rotor were calculated according to Strategy-1, Strategy-2, and Strategy-3, respectively, as shown in Table 12, Table 13 and Table 14. θz1, θz2, θz3, and θz4 were the optimal assembly orientations of the front shaft, the LPC, the HPC, and the back shaft, respectively, without considering the distribution of the screw holes of each component. The mean value and standard deviations of the optimal results were provided in Table 15. Furthermore, the optimal assembly orientations where the screw hole exists were obtained by a secondary calculation. As shown in Table 16, the optimal assembly orientations achieved during six iterations were exactly same. This indicates that the measurement errors of each parameter had little influence on the optimal results, especially when considering the existence of the distributed screw holes of each component.
The four-stage rotor was assembled according to the above three optimization results, and the unbalance of the assemblies was measured by a dynamic balancing machine further. The dynamic balancing machine was manufactured by Shanghai Shenzhong Co., Ltd. (Shanghai, China), and its minimum residual unbalance reaches 0.5 g·mm/kg. The experimental device was shown in Figure 10, and the corresponding results of the experiments were shown in Table 17.

5. Results and Discussion

Comparing the calculated unbalance in Table 16 and the measured results in Table 17, the simulation results were in good agreement with the experimental results. As shown in Table 17, when the θz2, θz3, and θz4 were 120°, 105°, and 120°, respectively (i.e., the calibrated screw holes of the LPC, HPC, and back shaft should rotate four, seven, and four distributed angles), the unbalance of the four-stage rotor reached the minimum value (184.6 g·mm). When the θz2, θz3, and θz4 were 150°, 60°, and 60°, respectively (i.e., the calibrated screw holes of the LPC, HPC, and back shaft should rotate five, four, and two distributed angles), the unbalance of the four-stage rotor reached the maximum value (404.1 g·mm). When the θz2, θz3, and θz4 were all 0° (i.e., the initial state), the unbalance of the four-stage rotor was 404.1 g·mm. Compared with the maximum unbalance and the initial unbalance, the minimum unbalance was notably reduced by 54.3% and 43.3%, respectively. These phenomena revealed that the proposed method of unbalance optimization for multi-stage rotor in this paper was effective.
The contents of the unbalance optimization method were mainly reflected in three aspects. First, an accurate assembly error propagation model is important to calculate the unbalance of a multi-stage rotor after assembly. The relevant simulation results in Section 3.1 showed that the assembly error propagation model developed in Section 2.1 was able to accurately predict the coordinates of any point in the assembled rotor. Moreover, the distribution of the screw holes of each single-stage rotor was not considered in the existing assembly error propagation models. In the actual assembly process of the multi-stage rotor, the assembly orientations of each single-stage rotor are finite and discrete. This is not simply to turn the continuous variables (0−180°) into the discrete variables (i.e., 180° divided by the number of the screw holes). In order to calculate the assembly cumulative error under any corresponding relationship of the screw holes of rotors at different stage, the sampling angles (φn) between the calibrated screw hole and the highest point of each single-stage rotor must be known. Therefore, the sampling angles were introduced into the error propagation model in our study to calculate the optimal corresponding relationship of each rotor. With comparison to the existing studies, the distribution of the screw holes of rotors at different stages was directly introduced into the model for the first time, which increases the suitability for guidance on actual assemblies.
The second aspect is concerned with the selection and calculation of the assembly datum for unbalance optimization. Considering the actual operating conditions of dynamic balancing for a rotor, two-plane unbalance measured by the dynamic balancing machine must be perpendicular to the axis of rotation of the measured rotor. Such an axis of rotation is simply a line connecting the midpoints of the journals of the front and back shaft. In current studies, the assembly datum for calculating unbalance of a multi-stage rotor was selected as the same as the measured datum for geometric parameter. It is clearly unreasonable that the simulated assembly axis is inconsistent with the axis of rotation of the dynamic balancing test. In addition, the results of the simulation are not comparable with the results of the experiment. In this study, we solved the problem by selecting the “orientation-varying axis” in Figure 3 as the assembly datum for unbalance optimization. The proposed assembly datum can ensure the constraint conditions of the unbalance optimization have good consistency with those of the dynamic balance test.
In terms of algorithm optimization, the bionics algorithm is now seldom applied in assembly optimization. Here, a GA was employed to achieve the optimal matching of assembly orientations for rotors at different stages. Such a strategy is applicable not only to the optimization of continuous angles but also to the optimization of discrete angles in the presence of screw holes for assembly.

6. Conclusions

Unbalance caused by improper assembly is the main reason for vibration faults in multi-stage rotors. A method is urgently needed to not only predict the state of a multi-stage rotor after assembly, but also to find the optimal matching relationship of a rotor at different stages to minimize the unbalance of the rotor. In this study, we developed an unbalance optimization method of a multi-stage rotor based on an assembly error propagation model. The main contributions of this study can be summarized as follows:
  • An accurate assembly error propagation model was developed to calculate the coordinates of any point in a rotor at different stages after assembly. The alignment process and distribution of the screw holes of adjacent rotors were considered in this model for the first time.
  • A new assembly datum for unbalance optimization was proposed to ensure consistency with the actual conditions of the dynamic balance test.
  • Unbalance optimization of a multi-stage rotor was achieved using a GA. Additionally, the optimal assembly orientations of rotors at different stages were also identified.

Author Contributions

Conceptualization, Y.C. and J.C.; methodology, Y.C.; software, Y.C.; validation, Y.C.; formal analysis, Y.C.; investigation, Y.C.; resources, J.C.; data curation, Y.C.; writing—original draft preparation, Y.C.; writing—review and editing, Y.C.; visualization, Y.C.; supervision, X.S.; project administration, J.C.; funding acquisition, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Outstanding Youth Project of Natural Science Foundation of Heilongjiang Province, Grant JQ2019E002. The authors are grateful for the financial support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge Cui, S. in the College of Mechanical Engineering, Tianjin University of Science and Technology for the active support of the research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ding, S.; Zheng, X.; Bao, J.; Zhang, J. A comprehensive study of three dimensional deviation analysis methods for aero-engine rotors assembly. IOP Conf. Ser. Mater. Sci. Eng. 2019, 688, 033039. [Google Scholar] [CrossRef]
  2. Qin, Z.; Han, Q.; Chu, F. Bolt loosening at rotating joint interface and its influence on rotor dynamics. Eng. Fail. Anal. 2016, 59, 456–466. [Google Scholar] [CrossRef]
  3. Wang, H.; Ge, Q. Aero structural assembly deformation prediction considering drilling-induced stresses. J. Aerosp. Eng. 2015, 28, 04015007. [Google Scholar] [CrossRef]
  4. Wang, M.; Han, Q.; Wen, B. Modal characteristics and unbalance responses of fan rotor system with flexible support structures in aero-engine. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2017, 231, 1686–1705. [Google Scholar] [CrossRef]
  5. Hu, X.; Gao, F.; Cui, C. Analysis of mass unbalance torque on a spinning superconducting rotor. IEEE Trans. Appl. Supercond. 2013, 24, 3600204. [Google Scholar] [CrossRef]
  6. Ercoli, L.; La, M.S. The development of a dynamic balancing machine. Int. J. Mech. Eng. Educ. 1998, 26, 149–155. [Google Scholar] [CrossRef]
  7. Khazali, A.; Hisham, A.H.; Askari, M.R. The experimental analysis of vibration monitoring in system rotor dynamic with validate results using simulation data. Isrn Mech. Eng. 2014, 2012, 1–17. [Google Scholar] [CrossRef] [Green Version]
  8. Hussain, T.; Yang, Z.; Popov, A.A. Straight-build assembly optimization: A method to minimize stage-by-stage eccentricity error in the assembly of axisymmetric rigid components (two-dimensional case study). J. Manuf. Sci. Eng. 2011, 133, 031014. [Google Scholar] [CrossRef]
  9. Hussain, T.; Yasinshaikh, G.; Shaikh, S.A. Variation propagation control in straight-build assemblies: 2d case study. Mehran Univ. Res. J. Eng. Technol. 2013, 32, 71–80. [Google Scholar]
  10. Yang, Z.; Hussian, T.; Popov, A.A. A comparison of different optimization techniques for variation propagation control in mechanical assembly. IOP Conf. 2011, 26, 012017. [Google Scholar] [CrossRef] [Green Version]
  11. Yang, Z.; Hussian, T.; Popov, A.A. Novel optimization technique for variation propagation control in an aero-engine assembly. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2010, 225, 100–111. [Google Scholar] [CrossRef]
  12. Hussain, T.; Mcwilliam, S.; Popov, A.A. Geometric error reduction in the assembly of axis-symmetric rigid components: A two-dimensional case study. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2012, 226, 1259–1274. [Google Scholar] [CrossRef]
  13. Yang, Z.; Popov, A.A.; Mcwilliam, S. Variation propagation control in mechanical assembly of cylindrical components. J. Manuf. Syst. 2012, 31, 162–176. [Google Scholar] [CrossRef] [Green Version]
  14. Yang, Z.; Popov, A.A.; Mcwilliam, S.; Hussian, T. Dimensional variation propagation analysis in straight-build mechanical assemblies using a probabilistic approach. J. Manuf. Syst. 2013, 32, 348–356. [Google Scholar] [CrossRef]
  15. Jin, S.; Ding, S.; Li, Z. Point-based solution using Jacobian-Torsor theory into partial parallel chains for revolving components assembly. J. Manuf. Syst. 2017, 46, 46–58. [Google Scholar] [CrossRef]
  16. Sun, C.; Li, C.; Liu, Y. Prediction method of concentricity and perpendicularity of aero engine multistage rotors based on PSO-BP neural network. IEEE Access 2019, 7, 132271–132278. [Google Scholar] [CrossRef]
  17. Wang, L.; Sun, C.; Tan, J. Improvement of location and orientation tolerances propagation control in cylindrical components assembly using stack-build assembly technique. Assem. Autom. 2015, 34, 358–366. [Google Scholar] [CrossRef]
  18. Sun, W.; Li, T.; Yang, D. Dynamic investigation of aeroengine high pressure rotor system considering assembly characteristics of bolted joints. Eng. Fail. Anal. 2020, 112, 104510. [Google Scholar] [CrossRef]
  19. Piskin, A.; Aktas, H.E.; Topal, A. Rotor balancing with turbine blade assembly using ant colony optimization for aero-engine applications. Int. J. Turbo Jet-Engines 2017, 8, 0060. [Google Scholar] [CrossRef]
  20. Liu, Y.; Zhang, M.; Sun, C. A method to minimize stage-by-stage initial unbalance in the aero engine assembly of multistage rotors. Aerosp. Sci. Technol. 2018, 85, 270–276. [Google Scholar] [CrossRef]
  21. Sun, C.; Liu, Z.; Liu, Y. An adjustment method of geometry and mass centers for precision rotors assembly. IEEE Access 2019, 7, 169992–170002. [Google Scholar] [CrossRef]
  22. Yang, X.S. Nature-Inspired Meta-Heuristic Algorithms; Luniver Press: Beckington, UK, 2008; pp. 314–322. [Google Scholar]
Figure 1. The initial state of a 3-D two-stage simulation rotor before assembly.
Figure 1. The initial state of a 3-D two-stage simulation rotor before assembly.
Applsci 11 00887 g001
Figure 2. The assembly processes of a 3-D two-stage simulation rotor.
Figure 2. The assembly processes of a 3-D two-stage simulation rotor.
Applsci 11 00887 g002
Figure 3. Assembly datum for unbalance optimization of four-stage rotor.
Figure 3. Assembly datum for unbalance optimization of four-stage rotor.
Applsci 11 00887 g003
Figure 4. The center coordinates of the top mounting surface of the assembled three-stage rotor were measured using software.
Figure 4. The center coordinates of the top mounting surface of the assembled three-stage rotor were measured using software.
Applsci 11 00887 g004
Figure 5. The convergence progress of the genetic algorithm (GA)-based unbalance optimization. (a) Taking the minimum unbalance as the optimization objective. (b) Taking the maximum unbalance as the optimization objective.
Figure 5. The convergence progress of the genetic algorithm (GA)-based unbalance optimization. (a) Taking the minimum unbalance as the optimization objective. (b) Taking the maximum unbalance as the optimization objective.
Applsci 11 00887 g005
Figure 6. Selection of the optimal orientations between two optional orientations.
Figure 6. Selection of the optimal orientations between two optional orientations.
Applsci 11 00887 g006
Figure 7. The section view of the four-stage rotor.
Figure 7. The section view of the four-stage rotor.
Applsci 11 00887 g007
Figure 8. Measurement of the geometric parameters of a high-pressure compressor (HPC). (a) Measured scene. (b) Schematic diagram.
Figure 8. Measurement of the geometric parameters of a high-pressure compressor (HPC). (a) Measured scene. (b) Schematic diagram.
Applsci 11 00887 g008
Figure 9. The setting of the initial unbalance of the four-stage rotor.
Figure 9. The setting of the initial unbalance of the four-stage rotor.
Applsci 11 00887 g009
Figure 10. Measurement of the unbalance of the four-stage rotor.
Figure 10. Measurement of the unbalance of the four-stage rotor.
Applsci 11 00887 g010
Table 1. The geometric parameters of the three-stage rotor [17].
Table 1. The geometric parameters of the three-stage rotor [17].
StageHeight (Hn) [mm]Measured Radius (Dn)
[mm]
Concentricity (Cn) [mm]Parallelism (Pn) [mm]Sampling Angle (φn)
[°]
1701000.0050.0050
2701000.0050.0050
3701000.0050.0050
Table 2. The center coordinates of the top mounting surface of the assembled three-stage rotor.
Table 2. The center coordinates of the top mounting surface of the assembled three-stage rotor.
StageCalculated Results [mm]Measured Results [mm]
1(0.005, 0, 70)(0.005, 0, 70)
2(0.0076, 0.0025, 140)(0.0076, 0.0025, 140)
3(0.0043, 0.0066, 210)(0.0043, 0.0066, 210)
Table 3. The mass parameters of the three-stage rotor measured using software.
Table 3. The mass parameters of the three-stage rotor measured using software.
StageMass [kg]Initial Coordinates of the Center of Mass [mm]
117.2630(0.0036, 0, 35)
217.2630(0.0036, 0, 35)
317.2630(0.0036, 0, 35)
Table 4. Results of the GA-based unbalance optimization.
Table 4. Results of the GA-based unbalance optimization.
Optimization Objectiveθz1 [°]θz2 [°]θz3 [°]Unbalance [g·mm]
Minimum unbalance0018032.2568
Maximum unbalance0180180129.6123
Table 5. The measured geometric parameters of the front shaft.
Table 5. The measured geometric parameters of the front shaft.
SequenceCn [mm]Pn [mm]Dn [mm]φn [°]Hn [mm]
10.01530.013050.01307996.5410
20.01610.014750.01368196.5425
30.01580.011050.01227996.5480
40.01390.012350.01457996.5438
50.01430.012850.01718096.5450
60.01520.013150.01258196.5422
Table 6. The measured geometric parameters of the low pressure compressor (LPC).
Table 6. The measured geometric parameters of the low pressure compressor (LPC).
SchemeCn [mm]Pn [mm]Dn [mm]φn [°]Hn [mm]
10.09440.0235101.0910174105.1940
20.09120.0194101.0921175105.1915
30.09350.0200101.0850174105.1970
40.09900.0223101.0915174105.1984
50.09250.0219101.0946175105.1960
60.09500.0231101.0972176105.1921
Table 7. The measured geometric parameters of the high pressure compressor (HPC).
Table 7. The measured geometric parameters of the high pressure compressor (HPC).
SequenceCn [mm]Pn [mm]Dn [mm]φn [°]Hn [mm]
10.09630.015136.019013114.7130
20.09270.014336.018812114.7189
30.09300.014836.020412114.7205
40.09790.016736.017512114.7150
50.09110.015636.017013114.7137
60.09540.015236.019513114.7118
Table 8. The measured geometric parameters of the back shaft.
Table 8. The measured geometric parameters of the back shaft.
SequenceCn [mm]Pn [mm]Dn [mm]φn [°]Hn [mm]
10.04960.011150.120021080.2250
20.04850.012550.122521080.2245
30.05010.010550.123821180.2210
40.05120.013450.115021180.2265
50.04950.011550.122021180.2193
60.04900.012850.123121280.2200
Table 9. The standard uncertainties of different measured parameters.
Table 9. The standard uncertainties of different measured parameters.
Componentu(Cn) [μm]u(Pn) [μm]u(Dn) [μm]u(φn) [°]u(Hn) [μm]
Front shaft0.350.490.740.401.02
LPC1.100.681.670.331.13
HPC1.040.330.520.221.41
Back shaft0.380.450.320.311.22
Table 10. The distributed numbers of the screw holes of each component.
Table 10. The distributed numbers of the screw holes of each component.
ComponentDistributed Numbers
(Top Mounting Surface)
Distributed Numbers
(Bottom Mounting Surface)
Front shaft1212
LPC2412
HPC1224
Back shaft1212
Table 11. The initial unbalance parameters of the four-stage rotor.
Table 11. The initial unbalance parameters of the four-stage rotor.
Component Unbalanced Mass [G]Phase Relative to Calibrated Screw Hole [°]Effect Radius [Mm]Axial Coordinate of The Unbalanced Mass Point Relative to Measured Datum [Mm]
Front shaft205896.5
LPC2064105
HPC20900
Back shaft20510
Table 12. The optimal results were calculated according to Strategy-1.
Table 12. The optimal results were calculated according to Strategy-1.
Sequenceθz1 [°]θz2 [°]θz3 [°]θz4 [°]Unbalance [g·mm]
10120.01116.03123.00181.15
20120.00118.04122.00180.90
30120.04117.04121.00182.46
40120.00117.04121.00182.50
50120.01117.03122.00181.53
60120.01118.04121.00181.60
Table 13. The optimal results were calculated according to Strategy-2.
Table 13. The optimal results were calculated according to Strategy-2.
Sequenceθz1 [°]θz2 [°]θz3 [°]θz4 [°]Unbalance [g·mm]
10140.8045.2030.17402.15
20143.9245.0730.01401.90
30141.9945.0030.00402.46
40141.9245.0730.04402.50
50141.9745.0230.01401.53
60142.9545.0430.00404.60
Table 14. The optimal results were calculated according to Strategy-3.
Table 14. The optimal results were calculated according to Strategy-3.
Sequenceθz1 [°]θz2 [°]θz3 [°]θz4 [°]Unbalance [g·mm]
10000324.60
20000328.13
30000324.92
40000324.90
50000327.14
60000329.72
Table 15. The mean value and standard deviation of the optimal results.
Table 15. The mean value and standard deviation of the optimal results.
Optimal ResultsStrategy-1Strategy-2Strategy-3
Mean ValueStandard DeviationMean ValueStandard DeviationMean ValueStandard Deviation
θz1 [°]000000
θz2 [°]120.010.0147142.261.061700
θz3 [°]117.200.756345.070.070900
θz4 [°]121.670.816530.040.066200
Unbalance [g·mm]181.690.6631402.521.0800326.571.1011
Table 16. The optimal results with considering the distributed screw holes of each component.
Table 16. The optimal results with considering the distributed screw holes of each component.
Strategyθz1 [°]θz2 [°]θz3 [°]θz4 [°]Unbalance [g·mm]
10120(4 × 30)105(7 × 15)120(4 × 30)189.68
20150(5 × 30)60(4 × 15)60(2 × 30)403.07
30000320.65
Table 17. The results of three experiments.
Table 17. The results of three experiments.
Strategyθz1 [°]θz2 [°]θz3 [°]θz4 [°]Unbalance [g·mm]
10120(4 × 30)105(7 × 15)120(4 × 30)184.6
20150(5 × 30)60(4 × 15)60(2 × 30)404.1
30000325.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Y.; Cui, J.; Sun, X. An Unbalance Optimization Method for a Multi-Stage Rotor Based on an Assembly Error Propagation Model. Appl. Sci. 2021, 11, 887. https://0-doi-org.brum.beds.ac.uk/10.3390/app11020887

AMA Style

Chen Y, Cui J, Sun X. An Unbalance Optimization Method for a Multi-Stage Rotor Based on an Assembly Error Propagation Model. Applied Sciences. 2021; 11(2):887. https://0-doi-org.brum.beds.ac.uk/10.3390/app11020887

Chicago/Turabian Style

Chen, Yue, Jiwen Cui, and Xun Sun. 2021. "An Unbalance Optimization Method for a Multi-Stage Rotor Based on an Assembly Error Propagation Model" Applied Sciences 11, no. 2: 887. https://0-doi-org.brum.beds.ac.uk/10.3390/app11020887

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