1. Introduction
The optimization of aero-engine rotor unbalance is an important part of the assembly process [
1,
2,
3]. Blades need to work under a high-load, high-speed, and high-vibration environment, which directly affects the start-stop performance, working reliability, efficiency, and cost of the aero-engine. The blades have the characteristics of complex shapes and high precision requirements. Therefore, in the production process, in order to meet the technological requirements, the blades need to be manually polished, so that the mass moment of each blade after assembly is different, resulting in the overall rotor having a large residual unbalance. Therefore, a reasonable blade arrangement can reduce the unbalance of the rotor, and a large number of scholars have done research in this area [
4]. Piskin et al. introduced the grouping methods of blade assembly. These methods optimize the imbalance of the rotor and provide a basis for the application of intelligent algorithms [
5,
6]. Lavagnoli developed a numerical strategy to select the best blade arrangement around the rotor, taking into account the measured blade weight distribution and any residual disk unbalance, which makes the algorithm calculate the best blade sort in the fan-constrained rotor row in a few seconds [
7]. Amiouny developed a heuristic of a problem that models the static balance of a turbofan: load point mass at regularly spaced locations on the circumference so that the residual imbalance around the center corresponding to the fan’s axis of rotation is as small as possible [
8]. Pitsoulis proposed a heuristic algorithm to solve an NP-Hard combinatorial optimization problem in turbine engine manufacturing and maintenance [
9].
Li et al. used the ant colony algorithm to find the blade arrangement, and successfully reduced the impeller unbalance over one year to 1% of the standard value, and the impeller unbalance over 99% was reduced to 1% of the official value [
10,
11]. Mason used a neighborhood search algorithm to determine the best position of the turbine blades, so that the total time required for balancing could be significantly reduced. It also showed that using the starting points obtained from the Lagrangian dual scheme can greatly improve the results of the problem [
12]. Thompson et al. used the simulated annealing algorithm to solve the large-scale combinatorial optimization problem of seeking the best layout through blade exchange [
13,
14,
15]. Dai et al. applied it to the optimal arrangement of blades of impeller machinery on the basis of the genetic algorithm, which solved the problem of a long solution time and an inability to obtain optimization results that were encountered in the previous optimization of a large number of blades by the exhaustive method [
16,
17,
18]. Pan et al. improved the genetic algorithm based on the genetic algorithm, and proposed an improved genetic algorithm to solve the assembly sort optimization problem [
19,
20,
21]. Wang improved the non-dominated sorting genetic algorithm (NSGA), introduced the control elite and dynamic crowding distance, and showed good performance when dealing with multi-objective optimization problems of wind turbines [
22]. In this paper, we use a cloud adaptive genetic algorithm (CAGA) to solve the blade sorting problem for the rotor unbalance situation, and its convergence speed and sorting results are explained.
The structure of this paper is as follows. Firstly, the vector sum model of blade mass moment and eccentric moment is established to obtain the unbalanced optimization function. Secondly, the cloud adaptive genetic algorithm is used to optimize the rotor unbalance, and the blade sort corresponding to the minimum unbalance is given. Finally, combined with the blade mass moment weighed by the torque, the effectiveness of the method is verified by comparing six commonly used blade arrangement methods.
3. Cloud Adaptive Genetic Algorithm
The genetic algorithm (GA) was first proposed by the American professor Holland in 1975 [
23]. It is a kind of random search algorithm that draws lessons from natural selection and the natural genetic mechanism in the biological world. The genetic algorithm simulates the selection, crossover, and mutation phenomena of natural selection and the genetic process. Each iteration of blades with a set of candidate solutions, the optimal chromosomes are selected according to the indicators, and these chromosomes are combined with genetic operators to generate a new generation of candidate groups. The process is repeated until the convergence index is satisfied. The specific flow chart is shown in
Figure 7.
The initial population is the size of the search area, which determines the search time and search accuracy. Fitness is the only measure of population evolution. The greater the population fitness, the higher the accuracy of the search results. The number of iterations is used as the number of population updates. The selection operation selects individuals with high fitness from the parent population as the offspring. The crossover operation randomly selects the two chromosomes left by the selection operation for gene exchange, thereby forming a new individual and placing it in the offspring. The mutation operation randomly selects the chromosomes in the selection operation for gene mutation, to form a new individual. Both crossover and mutation operations have crossover probability and mutation probability to determine whether to perform the operation.
Genetic algorithm has the characteristics of a low order, short defining moment, and high average fitness. If the probability of crossover and mutation is small and the fitness is low, it is difficult to produce excellent new chromosomes. If a large crossover and mutation probability are adopted, the population will easily fall into local optimization. Therefore, the adaptive genetic algorithm (AGA) appeared, while AGA only considers the trend of the evolutionary process and ignores the randomness of environmental evolution.
The cloud adaptive genetic algorithm improves the above defects. It adopts the cloud model theory [
24]. The cloud model is a set of random numbers with stable regularity following the regular distribution law, represented by the expected value (
Ex), entropy (
En), and hyper entropy (
He). The cloud model makes use of the uncertain transformation model between qualitative and quantitative concepts of linguistic values, which is both random and fuzzy, and provides an effective means for the combination of qualitative and quantitative information. It is mainly reflected in the crossover and mutation operators of the genetic algorithm.
3.1. The Coding and Operation of the Genetic Algorithm for Optimal Blade Sorting
In this paper,
N sets of blades were randomly generated as the initial population. Each blade is equivalent to a chromosome, which has a certain code. Sequential coding is adopted in blade sorting, that is, there is a one-to-one correspondence between the blade and feasible position. Since the wheel disk is a circle, it is stipulated that the position that coincides with the positive direction of the X-axis is the first position, and the counterclockwise direction is the positive direction of blade installation. The sorting scheme given in
Figure 2 is [6, 8, 4, 2, 5, 7, 3, 1].
3.1.1. Selection
The selection reflects the survival principle of the fittest. It retains large fitness and eliminates low fitness chromosomes. Its role is to avoid gene deletion, and improve the global convergence and computational efficiency. In this paper, the operator is selected using roulette, and chromosomes with large fitness are highly likely to be selected. For each chromosome, if its fitness value is
f(
xi), then the relative value of its fitness is:
where
pop is the size of the population, and we take
A as the selection probability of selecting this chromosome. The specific operation of the roulette blade selection is: Accumulate the chromosome’s selection probability(
pi) one by one to get the cumulative probability(
pj), and divide the interval [0–1] into pop intervals. Then, generate a random number belonging to [0–1]. If the number falls in the interval [
pj−1,
pj], select the chromosome corresponding to
pj. For example, if the selection probabilities of chromosomes
A,
B,
C, and
D are 0.1, 0.2, 0.3, and 0.4, then their cumulative probabilities set the interval to [0–0.1], [0.1–0.3], [0.3–0.6], and [0.6–1]. If the random number is 0.5, then the chromosome
C is selected.
3.1.2. Crossover
Crossover allows the exchange of genetic material between chromosomes to produce better chromosomes. In this paper, a blade is regarded as a gene in a chromosome, so two blades of the same code cannot appear in a chromosome. General single-point crossover operators and multi-point crossover operators are obviously not eligible, so we use the recombination crossover operator. This method first builds an edge list for the blades in the two parent chromosomes, indicating the blades connected to this blade and the number of occurrences. Assume two parental chromosomes are:
The two chromosomes are crossed, and the resulting edge list is shown in
Table 1. If an edge appears twice in the parents, a “−” sign is added to the vertex of the edge in the list. The edge recombination crossover operator starts constructing offspring by selecting an initial point.
We can choose the first position of the parent (P
1) as the starting point of the descendants. Then, the offspring (C
1) is:
The principle of selecting a chromosome from the parents is to select the blade with the least number of blades adjacent to it. If the number of blades connected to a certain two are equal, the blade with “−” is preferred. If the above conditions are the same for both blades, then one of them is selected at random. As shown by P
1 and P
2, there are two, six, and eight blades connected to blade 1. They have the same number of blades, but there is a “−” sign before eight, so eight is chosen. There are one, three, and seven blades connected to eight, and one has been selected before, so we chose between three and seven; because there are more adjacent blades to seven, we chose three. Repeatedly, we can get:
The crossover probability(
pc) is given by the cloud adaptive genetic algorithm:
3.1.3. Mutation
The mutation can restore the lost or untapped genetic material of the chromosome to prevent the chromosome from converging prematurely during the formation of the optimal solution. For the problem of blade arrangement, this paper uses a two-element optimization mutation operator. That is, the chromosomes selected according to the mutation probability are randomly selected at a certain position of the chromosome to exchange blades, and a new blade order is obtained. All the chromosomes of the parents are mutated according to the mutation probability until the next generation population is generated. The mutation probability(
pm) is given by the cloud adaptive genetic algorithm:
where
RANDN(
En,
He) generates a normal random number with an expected value of
En and a standard deviation of
He.
fmax is the maximum fitness of the population,
is the average fitness,
f′ is the fitness of the mutant chromosomes, and
f is the larger value of the fitness of the crossed chromosomes,
k1~
k4 ε[0~1]. In this paper,
k1 =
k3 = 1.0,
k2 =
k4 = 0.5,
c1 = 2.9,
c3 = 3.0,
c2 =
c4 = 10.
3.2. Fitness Function
According to the principle of biological evolution, the larger the fitness requirement, the better, and it is non-negative. Combined with the total unbalance given above, the fitness function given in this paper is as follows:
4. Data Validation
In order to view the actual optimization effect and calculation speed of the algorithm, the blades of the first-stage rotor of an aeroengine were selected for cloud adaptive genetic algorithm assembly optimization. Due to special requirements, we cannot disclose the test pieces of blades and disk, so we used the rotor in
Figure 8 as the experimental schematic diagram. There are 32 blades in this group. The mass moment of this group of blades was measured on the mass moment scale (MW0). The value is shown in
Table 2.
According to the mass moments given in
Table 2, it was assumed that the disk unbalance is 0 g.mm. The cloud adaptive genetic algorithm was used to sort the blades. The set algorithm parameters were as follows: the population size
N is 30; the maximum number of iterations is 500. The C sharp language was used to program the algorithm, and the algorithm was run on an ordinary PC with I7 2.4 GHz 4 GB running memory. It was found that the minimum unbalance of the population obtained varies with the number of iterations as shown in
Figure 9.
According to
Figure 9, we can see that when the number of iterations reaches 235, the maximum fitness result was achieved. For general or unreasonable sorting, the rotor unbalance amount can reach 6000 g.mm or more. As the number of iterations of the algorithm increases, the unbalance becomes smaller and smaller, until the minimum unbalance is 99.63 g.mm when the number of iterations reaches 235. The corresponding sorting result is shown in
Table 3.
In order to compare the optimization effect of the cloud adaptive genetic algorithm, the blades were sorted according to several methods (weight sorted, sorting type 2, sorting-1/4 skip) mentioned in the paper [
5], and the obtained unbalance sorting results are shown in
Figure 10a–c. According to sorting type 2, this paper designed sorting type 4, that is, sorting type 2 was used to re-distribute the blades in four sectors, and the resulting unbalance sorting result is shown in
Figure 10d. The unbalance sorting results obtained by the genetic algorithm (GA) and cloud adaptive genetic algorithm (CAGA) are shown in
Figure 10e,f.
Figure 10 describes the assembly position of each blade and the arrangement of the mass moment under the six blade arrangement methods. These methods can obtain the unbalance value under several blade arrangements. Assuming that a mass block needs to be added to the circumference of a disk with a radius of 500 mm for the balancing process, the required balance mass is recorded in
Table 4.
It can be seen from
Table 4 that the weight sorted method causes the largest unbalance in the assembly process. The sorting type 2, sorting type 4, and sorting-1/4 skip method reduce the unbalance after arrangement to a certain extent relative to the weight sorted. GA also optimized the unbalance to a certain extent, but it still cannot meet the requirement that the unbalance cannot be brought into the tolerance range from the process, and the mass of the corresponding mass that needs to be balanced does not meet the requirements. From the perspective of the balance method, the intelligent algorithm optimization effect is better than the empirical method. The unbalance of GA is at least 57.5% optimized relative to the first four methods, and the unbalance of the CAGA algorithm is optimized by 95.7% relative to GA.
In order to check the accuracy and calculation time of the cloud adaptive genetic algorithm and genetic algorithm, the results of the two methods when the number of iterations are 50, 100, 200, and 500, respectively, are shown in
Figure 11 and
Figure 12.
Figure 11 shows that CAGA is much better than GA in terms of the algorithm accuracy, and the imbalance of the two methods decreases with the increase of the number of times.
Figure 12 illustrates that CAGA is slower than GA in terms of the calculation time.
For the general blade sorting process, the initial unbalance of the disk is not considered. This approach reduces the unbalance of the blade sorting result, but due to the existence of the initial unbalance of the disk, the assembly unbalance after the sorting is not reduced. Therefore, the calculation results obtained by setting the disk unbalance in the X direction as 100, 200, and 500 g.mm in this paper are shown in
Table 5.
It can be seen from
Table 5 that at the same time, different initial unbalance can still get better results. When the unbalance amount is 200 g.mm, the unbalance amount after optimization even reaches 5.75 g.mm, which shows that the method is still effective for the optimization of blade sorting with the initial unbalance amount. The second and third columns of the table compare the final imbalance of the same blade sequence when the imbalance of
Table 5 is included, and the imbalance of the blade is only considered when the optimal rotor imbalance is reached. It can be seen from the table that the model that considers the imbalance of the leaf disc is much smaller than the imbalance that only considers the residual error of the blade, which proves the effectiveness of the model.