Next Article in Journal
Investigating the Obstacle Climbing Ability of a Coal Mine Search-and-Rescue Robot with a Hydraulic Mechanism
Next Article in Special Issue
Research and Evaluation Method for Echelon Usage of Double-Life Equipment Considering Calendar Life
Previous Article in Journal
Lipreading Using Liquid State Machine with STDP-Tuning
Previous Article in Special Issue
The Characteristics and Mechanisms of High-Intensity Sound in a High-Speed Multistage Compressor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and I-NSGA-III-VLC Solution of Aircraft Equipment Rotation and Echelon Usage under Uncertainty

School of Equipment Management and UAV Engineering, Air Force Engineering University, Xi’an 710051, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(20), 10482; https://doi.org/10.3390/app122010482
Submission received: 22 September 2022 / Revised: 9 October 2022 / Accepted: 11 October 2022 / Published: 17 October 2022
(This article belongs to the Special Issue Aircrafts Reliability and Health Management)

Abstract

:

Featured Application

In this study, an effective multi-objective optimization method (I-NSGA-III-VLC) is proposed and applied to the rotation and echelon usage problem of aircraft equipment under uncertainty, proving the effectiveness of the method.

Abstract

Optimizing the aircraft equipment usage scheme of different units according to their task intensity has great significance in improving aircraft reliability and health management. This paper studied the modeling and solving methods of the rotation and echelon usage problems of aircraft equipment measured by dual-life indexes, one of which cannot be controlled. In order to maximize the waste rate of the rotation quantity, echelon uniformity index, life matching index and life utilization index, a decision-making model of the equipment rotation and echelon usage problem under uncertainty was constructed, and an improved NSGA-III with a variable length chromosome was proposed. An improved segmented coding method and operators were proposed, and the repeated individual control mechanism was used to improve the population diversity. When the scale of the problem was large, this method could search a wider range in a short time and obtain more feasible solutions, which verified the feasibility of this method.

1. Introduction

The life of many pieces of aircraft equipment is measured by two indexes, namely calendar life (CL) and working life (WL), and this type of equipment is called dual-life equipment (hereinafter referred to as equipment). Similar to aviation-guided munitions, the “power on time” is the working life of aviation-guided munitions.
The specified available calendar life for maintenance is called the stage calendar life, and the specified available working life for maintenance is called the stage working life. The calendar life is gradually consumed over time regardless of whether the task is executed or not, and it is uncontrollable, while the working life is only consumed by the execution of the task, which is controllable. In the actual usage process, the number of tasks for each unit and the number of pieces of equipment rotated between units are uncertain due to the lack of a rotation and echelon usage scheme formulation method. The scheme selection method is random and subjective [1], resulting in the fact that when the calendar life of the equipment stage is exhausted, its working life only consumes 40~60% of its stage working life [2], and the working life is seriously wasted and the echelon state is difficult to control. In the case of no rotation, due to the different equipment states and working intensities of different units, the life utilization index and life matching index of equipment are greatly different. Therefore, an effective rotation and echelon usage control method has an important role in solving the problems associated with echelon uniformity, life utilization and the life-matching degree of each unit.
Liu [3] established a set of comprehensive evaluation indexes, including working life combat-readiness reserves and echelon-uniformity indexes. However, these indexes cannot adapt to the situation of considering calendar life, and cannot control the echelon status. Zhang [4] proposed an improved echelon usage diagram method, which is concise and intuitive. Zhou et al. [5] proposed a formation and control method of motorcycle working life echelon reserves using SAA to formulate a usage and repair plan of a motorcycle. For assignment problems, Dorterler [6] proposed a new genetic algorithm with agent-based crossover for the generalized assignment problem, but it is based on the idea that each task should be assigned to one agent only. Fu et al. [7] studied the uncertain multi-objective assignment problem, in which the number of tasks each person undertakes is uncertain, but each task can only be undertaken by one person. Xiong [8] solved the assignment problem of multiple people participating in a task with a new generalized nonlinear integer programming model and verified its effectiveness through an example. Ding [9,10] proposed an α-optimistic uncertain assignment model, but this model has the limitations of low efficiency and large storage space requirements.

2. Problem Description and Modeling

2.1. Problem Description

In the case of rotation, the echelon usage problem of equipment belongs to the uncertain assignment problem. The characteristics of this kind of assignment problem are that the amount of equipment in each rotation is uncertain, the quantity of tasks that each equipment can undertake is uncertain and the objective function result is uncertain before the decision matrix is determined.
According to the remaining-life data of equipment and task quantity in different units, the arrangement of the equipment to be rotated is calculated first, followed by the performance of a task of a months, where a is a positive integer, aiming at maximizing the average working-life echelon-uniformity index, the average calendar-life- and working-life-matching index, the average life-utilization-rate index, and the waste rate-of-rotation-quantity index. The result is the number of pieces of equipment rotated by each unit and the quantity of tasks completed by each piece equipment each month.
The average working-life echelon uniformity refers to the average value of the working-life echelon uniformity of each unit, hereinafter referred to as the average echelon uniformity. The average calendar-life- and working-life-matching-degree index refer to the average value of the calendar-life- and working-life matching-degree index of each unit, hereinafter referred to as the average matching-degree index. The average life-utilization rate refers to the average value of the life-utilization rate of each unit, hereinafter referred to as the average utilization-rate index. The waste rate-of-rotation-quantity index refers to the ratio of the unused rotation equipment quantity to the maximum equipment rotation quantity of all units.

2.2. Problem Assumptions

  • Assumption 1: There is complete interchangeability between the equipment;
  • Assumption 2: The quantity of tasks per month of each unit is fixed;
  • Assumption 3: The quantity of equipment required for each unit is fixed;
  • Assumption 4: The working-life consumption per task is fixed;
  • Assumption 5: The calendar life of each piece of equipment is an integer, with the unit of months;
  • Assumption 6: There is no limit to the number of times each piece of equipment can be used every month;
  • Assumption 7: Equipment failure is not considered during rotation and usage;
  • Assumption 8: The upper limit of the amount of equipment in each rotation is determined;
  • Assumption 9: The calendar life consumption caused by rotation is negligible.

2.3. Problem Inputs

  • Task quantity sequence every month of each unit U ;
  • Quantity of equipment required for each task v ;
  • Working-life consumption per task w ;
  • Stage working life of equipment M W L ;
  • Stage calendar life of equipment M C L ;
  • Task time range a ;
  • Maximum quantity of equipment for each rotation g ;
  • Equipment quantity sequence of each unit N ;
  • Unit quantity c .

2.4. Problem Ouputs

  • Equipment number of each rotation;
  • Number of pieces of equipment used for each task.

2.5. Mathematical Model

Taking unit u 1 as an example, when rotating from u 1 to other units, the rotation number is c 1 , and the rotation process is shown in Figure 1.
The decision matrix is divided into two parts.
Let us introduce the 0-1 variable x i , j ( i = 1 , 2 , , c 1 , j = S 1 , 1 , S 1 , 2 , , S 1 , N 1 , c 2 ) , where S 1 is the equipment number sequence of u 1 , N 1 is the equipment quantity of u 1 . By analogy, the 0-1 variable is also introduced for other units to build a rotation decision matrix. The selected equipment is moved to the corresponding unit, and the equipment numbers of the current unit and the corresponding unit are updated. S is the updated equipment number sequence of each unit and N is the updated equipment quantity sequence of each unit.
x i , j = 0   S i   i s   n o t   s e l e c t e d   i n   j t h   r o t a t i o n 1   S i   i s   s e l e c t e d   i n   j t h   r o t a t i o n
Let us introduce the 0-1 variable y i , j = ( i = 1 , 2 , , a U 1 ; j = 1 , 2 , , n ) , where U 1 is the task quantity of u 1 . By analogy, the 0-1 variable is also used for other units to construct the task decision matrix and update the working life of the selected equipment.
y i , j = 0   S i   i s   n o t   s e l e c t e d   i n   j t h   t a s k 1   S i   i s   s e l e c t e d   i n   j t h   t a s k
Next, let us establish the rotation decision matrix, X , and task decision matrix, Y , of u 1 , as shown in Table 1 and Table 2.
X = x 1 , 1 x 1 , N 1 x c 1 , 1 x c 1 , N 1
Y = y 1 , 1 y 1 , N 1 y a U 1 , 1 y a U 1 , N 1
Constraints:
s . t . j = 1 N 1 x i , j g 0 i = 1 a U 1 y i , j T 1 , j = F l o o r ( W L 1 , j w ) F l o o r ( W L 1 , j w ) U 1 C L 1 , j U 1 C L 1 , j F l o o r ( W L 1 , j w ) > U 1 C L 1 , j i = 1 c 1 x i , j 1 i = 1 N 1 y i , j = v
Update data:
W L 1 , j = W L 1 , j w i = 1 a U 1 x i , j C L 1 , j = C L 1 , j a
T i , j are the remaining usable times of S i , j , where W L and C L are the updated remaining working life sequence and calendar life sequence. The fitness functions of this problem are the average echelon-uniformity function of each unit, average life-matching function of each unit, average utilization function and rotation-quantity waste function.
1.
The average echelon-uniformity function, F 1 .
The improved echelon-uniformity function, Q 1 , of u 1 is shown in Formula (1)
Q 1 = 1 i = 1 N 1 1 U L 1 , i + 1 U L 1 , i e 1 N 1 2 2 + U L 1 , 1 + e 1 U L 1 , N 1 e 1 N 1 2 2 N 1 1 e 1 N 1 2 2 + e 1 e 1 N 1 2 2 × 100 %
This is continued for other units to acquire F 1 = i = 1 c Q i c .
U L 1 is the sequence of sorted W L 1 in ascending order according to the C L 1 of u 1 . Its purpose is to increase the working life along the direction of increasing calendar life. e 1 is the maximum value in U L 1 , and N 1 is the quantity of equipment of u 1 .
2.
The average life-matching function, F 2 .
The life-matching function, Z 1 , of u 1 is shown in Formula (2)
Z 1 = 1 i = 1 N 1 H 1 , i N 1 V 1 H 1 , i = M W L W L 1 , i M C L C L 1 , i M W L M C L                               C L i M C L ( M W L W L 1 , i ) M W L M C L         C L i = M C L
This is continued for other units to acquire F 2 = i = 1 c Z i c .
H 1 , i is the absolute value of the difference between the updated used working life and the used calendar life ratio and the ideal ratio of each equipment. When C L 1 , i = M C L , the used working life of the equipment is calculated as a ratio.
3.
The average utilization function, F 3 .
It is defined that the standard for equipment to enter into consumption is that the working life must be less than 3 h or the calendar life must be less than 1 month. When the equipment reaches the consumption standard, it can be considered that the utilization rate is 100%. The process is as follows:
For u 1 , define k 1 , 1 : = 0 , k 1 , 2 : = 0 , R W L 1 : = 0 , R C L 1 : = 0 .
Traverse the updated data.
When C L 1 , i : = 0 , k 1 , 1 : = k 1 , 1 + 1 ;
When W L 1 , i w , R W L 1 : = R W L 1 + W L 1 , i ;
When W L 1 , i < w , R W L 1 : = R W L 1 + 0 ;
When W L 1 , i < w , k 1 , 2 : = k 1 , 2 + 1 , R C L 1 : = R C L 1 + C L 1 , i .
The utilization function of u 1 is shown in Formula (3).
R 1 = k 1 , 1 M W L R W L 1 k 1 , 1 M W L + k 1 , 2 M C L R C L 1 k 1 , 2 M C L 2 100 %       k 1 , 1 > 0 & k 1 , 2 > 0 k 1 , 1 M W L R W L 1 k 1 , 1 M W L 100 %     k 1 , 1 > 0 & k 1 , 2 = 0 k 1 , 2 M C L R C L 1 k 1 , 2 M C L 100 %       k 1 , 2 > 0 & k 1 , 1 = 0
k 1 , 1 is the quantity of equipment with an exhausted calendar life of u 1 , k 1 , 2 is quantity of equipment with an exhausted calendar life of u 1 , R W L 1 is the sum of the remaining working life of the equipment with an exhausted calendar life after updating u 1 , and R C L 1 is the sum of the remaining calendar life of the equipment with an exhausted working life after updating u 1 .
This is continued for other units to acquire F 3 = i = 1 c R i c .
4.
The rotation-quantity waste function, F 4 .
The transportation cost and risk level caused by rotation can be replaced by the equipment rotation quantity. In order to reduce the equipment transportation cost and risk caused by rotation, the number of rotating pieces of equipment should be as small as possible.
Let us assume that the number of pieces of equipment rotated from a certain unit at each time cannot exceed g , the rotation phase can be divided into from u 1 to u 2 , from u 1 to u 3 ,…, from u 1 to u c ,…, and from u c to u c 1 , where the total time of rotations is c ( c 1 ) , and the maximum quantity of rotations is c ( c 1 ) g . Let us assume that the rotation quantity represented by an individual chromosome is P , the rotation-quantity waste rate function of the individual F 4 is shown in Formula (4).
F 4 = c ( c 1 ) g P c ( c 1 ) g
Aiming at the multi-objective assignment problem with an uncertain rotation quantity, with the task quantity completed by each equipment under the conditions of rotation, this paper constructed a decision-making model for the rotation and echelon usage of dual-life equipment and proposed an improved NSGA-III method with a variable length chromosome (I-NSGA-III-VLC) to solve it. This method can obtain a variety of optimization schemes in a relatively short time. Through scheme screening, the problem of poor working-life echelon status, low matching between calendar life and working life and the low utilization of the working life of equipment in different units can be solved to a certain extent with as few rotations as possible.

3. Principle of I-NSGA-III-VLC Method

Optimizing the equipment echelon uniformity, life-matching degree, life utilization and rotation-quantity waste rate at the same time is a high-dimensional multi-objective optimization problem. Deb proposed a fast, non-dominated-sorting genetic algorithm based on reference points (NSGA-III) [11], which replaced the congestion by associating the reference points, and solved the problem of high-dimensional objective optimization. Many scholars have solved high-dimensional multi-objective optimization problems based on the above algorithms, including the multi-objective classification problem [12], reservoir flood-control-operation problem [13], resource-allocation problem [14,15,16], location-routing problem [17,18,19,20], high-dimensional target power-flow-optimization problem of a power system [21,22], etc.
The feasible region of many problems is uncertain in real life. Many scholars [23,24,25,26,27,28,29] have improved the NSGA-III-VLC according to the characteristics of the problem. In the formulation of the equipment usage plan considering rotation, because the quantity of rotation will affect the rotation cost and risk, the length of the chromosome also affects the decision-making. In this paper, the chromosome-coding method, crossover operator, mutation operator and repetitive individual control operator were improved according to the characteristics of the problem. The specific improvement contents are shown from Section 3.1, Section 3.2, Section 3.3, Section 3.4 and Section 3.5.

3.1. Variable Length Cegment Hybrid Coding Method

Due to the need for considering the rotation between different units, the chromosome-coding method needs to be changed to distinguish different units. The rotation stage is placed before the task stage. This paper adopted the variable length and hybrid segmented-coding methods. Let us add a number “prefix” before each rotation, such as 1, 2, …, c ( c 1 ) . Therefore, the shortest chromosome length at the rotation stage is c ( c 1 ) , and the maximum quantity of equipment for each rotation is g . Therefore, the longest chromosome length is c ( c 1 ) ( g + 1 ) .
After the rotation phase, according to the task intensity sequence U and V of each unit, according to the unit quantity, let us add a number “prefix” before each unit task, such as c ( c 1 ) + 1 , c ( c 1 ) + 2 ,…, c 2 . Therefore, the total length of chromosomes in the task phase is i = 1 c N i a + c .
Each coding segment is divided by a digital “prefix”. With this coding method, we can not only distinguish each unit but also cross and mutate under the same “prefix”. The range of crossover and mutation is reduced, and it is easy to control the chromosome length. The chromosome initialization process is shown below.
Step 1: Arrange all equipment numbers S i and remaining usable times T in each unit to form a sequence in S i in the format of “Number: Remaining Usable Times”. Like all equipment numbers, S 1 of u 1 is S 1 = [‘A’: 4, ‘C’: 5, ‘E’: 5, ‘F’: 6, ‘B’: 6, ‘D’: 7], and like all equipment numbers, S 2 of u 2 is S 2 = [‘AA’: 3, ‘CC’: 6, ‘EE’: 7, ‘BB’: 7, ‘FF’: 8, ‘DD’: 9].
Step 2: Adopt a segmental coding method where there are c ( c 1 ) coding segments in the rotation phase, and the shortest length of each coding segment is 1 and the longest is g + 1 . Assume c = 2 , g = 3 . Randomly select l 1 equipment from u 1 and randomly select l 2 equipment from u 2 , 0 l 1 , l 2 g , where the equipment number and “prefix” together constitute the code segment, i.e., [1, ‘A’, ‘B’, 2, ‘CC’, ‘FF’]. This indicates that the first rotation is to rotate the equipment numbered ‘A’ and ‘B’ in u 1 to u 2 , and the second rotation is to rotate the equipment numbered ‘CC’ and ‘FF’ in u 2 to u 1 .
Add ‘A’ and ‘B’ to u 1 according to the remaining usable times, delete them in u 2 , add ‘CC’ and ‘FF’ to u 2 according to the remaining usable times, delete them in u 1 , update S 1 and S 2 , and obtain S 1 = [‘C’: 5, ‘E’: 5, ‘F’: 6, ‘CC’: 6, ‘D’: 7, ‘FF’: 8] and S 2 = [‘AA’: 3, ‘A’: 4, ‘B’: 6, ‘EE’: 7, ‘BB’: 7, ‘DD’: 9].
Step 3: In the task phase, the order of executing the task of the equipment is not considered in the same coding segment, but it cannot be reused. Gene j can be reused in different coding segments, but the number of repetitions cannot exceed the corresponding remaining usable times. The total length of code in a task stage is U i v . Assume U 1 = 2 , U 2 = 1 , v = 4 , initialize the chromosome, select v equipment from u 1 as a group, repeat U 1 times and select v equipment from u 2 as a group, repeat U 2 times, such as [3, [‘CC’, ‘C’, ‘E’, ‘F’], [‘C’, ‘E’, ‘FF’, ‘D’], 4, [‘AA’, ‘B’, ‘EE’, ‘DD’]]. This coding method can be directly interpreted as selecting the equipment to complete the task. The segmented real-number coding conforms to the practical significance, and it is easy to control the upper limit of usage times.
Step 4: Combine the chromosomes of the two stages to acquire C 1 = [1, ‘A’, ‘B’, 2, ‘CC’, ‘FF’, 3, [‘CC’, ‘C’, ‘E’, ‘F’], [‘C’, ‘E’, ‘FF’, ‘D’], 4, [‘AA’, ‘B’, ‘EE’, ‘DD’]]. By the same token, C 2 = [1, ‘C’, 2, ‘DD’, ‘AA’, 3, [‘AA’, ‘A’, ‘E’, ‘F’], [‘B’, ‘E’, ‘DD’, ‘D’], 4, [‘C’, ‘FF’, ‘EE’, ‘BB’]]. C 1 is shown in Figure 2, and C 1 is shown in Figure 3.

3.2. Population Initialization

Step 1: Population size M , p o p = [   ] , s : = 0 , equipment numbers S , such as S 1 = [‘A’: 4, ‘C’: 5, ‘E’: 5, ‘F’: 6, ‘B’: 6, ‘D’: 7] and S 2 = [‘AA’: 3, ‘CC’: 6, ‘EE’: 7, ‘BB’: 7, ‘FF’: 8, ‘DD’: 9];
Step 2: Initialize chromosome C , acquire S 1 and S 2 and calculate the usage time p of each equipment number in the chromosome in the format of “Number: Usage Times”. For example, p 1 of C 1 is [‘C’: 2, ‘E’: 2, ‘F’: 1, ‘CC’: 1, ‘D’: 1, ‘FF’: 1] and p 2 of C 1 is [‘AA’: 1, ‘B’: 1, ‘EE’: 1, ‘DD’: 1];
Step 3: Traverse the number in p i . If the usage time of the number is greater than the remaining usable time of the corresponding number in S 1 , return to Step 2. Otherwise, enter Step 4.
Step 4: Add C to the p o p , where s : = s + 1 . If s < M return to Step 2. Otherwise, the population initialization is completed.

3.3. Crossover Operator

For different coding stages, the ectopic single point crossover method was adopted in the rotation stage, and the two-point crossover method was adopted in the task stage. The crossover operator steps are as follows.
Step 1: Given the crossover rate p c , 0 p c 1 , the new population size is M , where n e w _ p o p = [   ] , and s : = 0 .
Step 2: Select two parent individuals and define c o u n t 1 : = 0 and c o u n t 2 : = 0 . The paternal chromosomes are segmented according to the number “prefix”, as shown in Figure 4.
Step 3: Combine the split parent individuals according to the same “prefix”. Because the first c ( c 1 ) combinations belong to the rotation stage, the length of coding segments may be different, so the heterotopic single point crossover method is adopted.
For the code segment combination in the rotation stage:
Step 4: Randomly select two points n u m 1 , j and n u m 2 , j , where n u m i , j is the point selected from the coding segment with “prefix” j of the parent i . In addition, 0 n u m 1 , j l e n 1 , j g , 0 n u m 2 , j l e n 2 , j g , 0 l e n 1 , j n u m 1 , j + n u m 2 , j g , and 0 l e n 2 , j n u m 2 , j + n u m 1 , j g , where l e n i , j is the length of the coding segment of the parent i with “prefix” j .
Step 5: Generate a random number r c between 0 and 1. If r c p c and c o u n t 1 < 10 , crossover the first n u m 1 , j genes of the coding segment with j as the “prefix” in parent1 and the first n u m 2 , j genes of the coding segment with j as the “prefix” in parent2. c o u n t 1 : = c o u n t 1 + 1 , generate crossover results, as shown in Figure 5, and enter Step 6. If p c < r c , do not crossover, return Step 4. If c o u n t 1 = 10 , the parent individuals will be preserved in the new population.
The specific steps of modifying operator 1 are as follows.
Step 6: Calculate the repetitions g i of each number after crossover in the format “Number: Repetitions”. Iterate through the numbers in g i . If the repetitions of the number are 1 and the set of numbers is a subset of the number S i , acquire S i and enter Step 7. Otherwise, return to Step 10.
For the code segment combination in the task stage:
Step 7: Adopt the method of two-point crossing and randomly select two points, p o t 1 and p o t 2 , where p o t 1 < p o t 2 , p o t 1 % v = 0 and p o t 2 % v = 0 .
Step 8: Generate a random number r c between 0 and 1. If r c p c and c o u n t 2 < 10 , exchange genes between two points, and if c o u n t 2 : = c o u n t 2 + 1 , generate crossover results, as shown in Figure 6, and enter Step 9. If p c < r c and c o u n t 2 < 10 do not cross, return to Step 7. If c o u n t 2 = 10 , the parent individuals are added to the new population, with s : = s + 2 .
The specific steps of modifying operator 2 are as follows.
Step 9: Calculate the number of repetitions b i , j of each character in each unit after crossover, where i represents the individual number of the offspring after crossover, and j represents the unit number in the format “Number: Quantity of Repetitions”. According to the different units, traverse the number in b i , j . If the number of repetitions is less than the number of repetitions of the corresponding number in S i of the corresponding unit and the number of b i , j is a subset of S i , enter Step 10. Otherwise, return to Step 7.
Step 10: Merge the cross results according to the “prefix” and add the results to the new population n e w _ p o p = [   ] , where s : = s + 1 . If s < M , return to Step 2. If s = M , the crossover operator ends.

3.4. Mutation Operator

For different coding stages, the methods of gene mutation and gene length mutation are used, respectively. The steps of the mutation operator are as follows.
Step 1: Define the mutation rate p m , 0 p m 1 .
Step 2: Select the parent individuals, define c o u n t 1 : = 0 , c o u n t 2 : = 0 , and divide the chromosome according to the number “prefix”.
For the combination of rotation stage:
Step 3: Randomly select two points n u m 1 , j and n u m 2 , j , where 1 n u m 1 l e n 1 g and 1 n u m 2 l e n 2 g .
Step 4: Generate a random number r m between 0 and 1. If r m p m , delete the number except n u m 1 in S 1 to get S S 1 , and randomly select 0 ~ g n u m 1 + 1 numbers from S S 1 to replace it. Similarly, delete the number except n u m 2 in S 2 to get S S 2 , and randomly select 0 ~ g n u m 2 + 1 numbers from S S 2 to replace it. Generate the mutation results, as shown in Figure 7, and enter Step 5. If p m < r m , do not generate the mutation, and return to Step 3.
Step 5: After obtaining the mutation results in the rotation stage, S i can be obtained according to the mutation results, such as S 1 = [‘AA’: 3, ‘E’: 5, ‘F’: 6, ‘B’: 6, ‘D’: 7] and S 2 = [‘A’: 4, ‘C’: 5, ‘CC’: 6, ‘EE’: 7, ‘BB’: 7, ‘FF’: 8, ‘DD’: 9].
For the code segment in the task phase:
Step 6: Traverse the code segment. If there is a number that does not belong to S i , it will be replaced with the number in S i according to the number of the remaining usage time.
Step 7: Randomly select the code segment of the parent to generate a random number r m between 0 and 1.
Step 8: If r m p m , enter Step 9, otherwise, return to Step 7.
Step 9: Calculate the number of gene repetitions excluding coding segment k i , and subtract k i from S i to acquire q i .
Step 10: This coding segment is mutated into a gene segment generated by the set q i , and no repeat genes are allowed.
This mutation operator ensures that the number repetitions of the chromosome gene after mutation is lower than the limit of the number repetitions, without using a new modifying operator, and improves the efficiency of the algorithm.

3.5. Population-Diversity-Preserving Operator

Step 1: Initialize the population, with the iteration number e p o c h : = 0 and the diversity-keeping parameter d .
Step 2: Conduct the crossover, mutation and selection operators to generate a new generation of population, where e p o c h : = e p o c h + 1 .
Step 3: Judge whether e p o c h % d is equal to 0. If the result is 0, enter Step 4; if the result is not 0, enter Step 5.
Step 4: Traverse all the individuals in the population, sequence the genes in the rotation phase of the individuals and count the number of gene repeats in the task phase. If there are individuals with the same rotation gene and the same number of repeats, the repeats will be replaced by a new, randomly generated individual.

3.6. Selection Operator

The selection operator of I-NSGA-III-VLC is used to select the population of the i iteration from the mixed population consisting of the population of the i iteration (the number is M ) and the new offspring obtained by the genetic operator (the number is M ). The individuals are sorted from largest to the smallest according to the non-dominated level [30], and the individuals in the same non-dominated set are selected according to the reference point [31,32,33,34,35,36,37]. A total of M individuals are obtained, which will be used as the next generation for the next iteration.

4. Examples and Results

In this example, there were two units, among which there were twenty pieces of equipment, respectively. Each task required four pieces of equipment and consumed their 3 h working life. The upper limit g of the quantity of equipment for each rotation was three. The four fitness functions were optimized under different circumstances by using enumeration method (①), improved genetic algorithm (②) and I-NSGA-III-VLC (③). A total of five experiments were carried out. The comparison of the time and the average value of the optimal solution of each fitness function is shown in Table 3, The population number of ② and ③ was set to 600, the crossover probability was set to 0.4, the mutation probability was set to 0.8, the diversity maintenance parameter was set to 6 and the number of iterations was set to 100 generations. M W L and M C L were set to 180 and 36, respectively. When the enumeration method had the problem of combination explosion, the best combination was selected from the first 100,000 combinations.
The two units u 1 and u 2 and the equipment data of u 1 and u 2 are shown in Table A1 and Table A2. The average monthly task intensities were U 1 = 6 and U 2 = 4. The iteration results of the 3-month equipment-echelon-usage problem using I-NSGA-III-VLC are shown in Figure 8.
By the weight formula β = α 1 F 1 + α 2 F 2 + α 3 F 3 + α 4 F 4 , where α 1 = 0.3 , α 2 = 0.1 , α 3 = 0.58 , α 4 = 0.02 , the acquired maximum value of β was 94.7%. The corresponding equipment number selected in the rotation phase is shown in Table 4. The equipment information of the two units and the number of tasks completed in 3 months are shown in Table 5 and Table 6.
The remaining working life and remaining calendar life of each equipment after the task was completed are shown in Figure 9 and Figure 10.
Equipment U and F in u 1 were rotated in u 2 , u 1 and formed a good echelon status.
By rotation, the equipment U reduced the life waste of u 2 , and u 2 had a good echelon status.

5. Discussion

It can be seen from Table 3 that the I-NSGA-III-VLC could obtain solutions equal to or even better than the enumeration method within a certain error and in a shorter time than other methods, with a higher efficiency, which means it had a larger search range. It was proven that the method was feasible in solving the problem of equipment echelon usage under the conditions of rotation.
It can be seen from Figure 8 that the I-NSGA-III-VCL could search for four fitness functions at the same time and obtain a Pareto optimal solution set. In addition to the three fitness functions represented by the coordinate system in Figure 8, there also was the number of rotations represented by the color, where rotation x representing the equipment quantity of rotations is x .
Using the usage scheme that maximizes the weight formula β , it can be seen from Figure 9 and Figure 10 that the remaining working life of the equipment of each unit formed a good echelon, and the average echelon uniformity index F 1 reached 86.4%, of which Q 1 was 86.4%, Q 2 was 86.4% and the average life matching degree F 2 reached 94.5%. The average life utilization rate of the remaining life of equipment F 3 was 100%, the life waste was the lowest and the total quantity of rotations was 2, which is acceptable.

6. Conclusions

In this paper, the echelon use problem of dual-life equipment considering rotation was modeled, and I-NSGA-III-VLC was proposed to solve this problem. The results of the example showed that this method could find multiple effective schemes. After scheme screening, the problem of low average echelon uniformity, average life-matching degree and average life-utilization rate caused by the extensive selection of dual-life equipment under rotation could be solved to a certain extent, and it provided a feasible and efficient method to build and solve the echelon usage problem model of dual-life equipment considering rotation. This problem was a method study. It may be impossible to guarantee that the usage of aircraft equipment will meet all the assumptions over a long time in real life, but if the time range is narrowed, the usage of aircraft equipment will meet the assumptions over a short time. In the future, research will continue under the conditions of uncertain task quantity and consumption, etc., and it is necessary to consider the difference in the quantity of equipment between units to prevent the situation in which the number of equipment is too small due to rotation and the task cannot be completed.

Author Contributions

Software, Q.M.; validation, Q.M., X.Z. and P.Q.; writing—original draft preparation, Q.M.; writing—review and editing, Q.M.; formal analysis, X.L.; investigation, Q.M., X.L.; resources, Q.M., X.Z.; data curation, Q.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Equipment data of u 1 .
Table A1. Equipment data of u 1 .
NumberRemaining Working LifeRemaining Calendar LifeNumberRemaining Working LifeRemaining Calendar Life
A15024K14524
B459L10018
C6513M356
D16026N17032
E9016O9017
F11519P16530
G7013Q11018
H6512R12020
I9518S8015
J253T13022
U203
Table A2. Equipment data of u 2 .
Table A2. Equipment data of u 2 .
NumberRemaining Working LifeRemaining Calendar LifeNumberRemaining Working LifeRemaining Calendar Life
AA14525KK11518
BB459LL356
CC7515MM16530
DD15027NN9517
EE8016OO15529
FF12520PP11018
GG6512QQ12520
HH7013RR10017
II9016SS13022
JJ14024

References

  1. Cheng, X.; Zhang, Q. Research on the aircraft usage decision based on BP neural network. Mach. Des. Manuf. Eng. 2014, 43, 76–78. [Google Scholar]
  2. Zhang, X.; Liu, X.; Zhang, T. Calendar life of aircraft based on annual flight frequency. J. Nanjing Univ. Aeronaut. Astronaut. 2017, 49, 56–59. [Google Scholar]
  3. Liu, Q.; Li, L.; Su, T. Evaluation indexes for echelon usage control of military aircraft. Ordnance Ind. Autom. 2016, 35, 24–27. [Google Scholar]
  4. Zhang, Q.; Wang, Y. Research in echelon usage of aircraft considering the calendar life. Aircr. Des. 2014, 34, 33–36. [Google Scholar]
  5. Wang, J.P.; Zhao, C.; Zhang, Z.W. Formation and control method of motor hour echelon storage for a kind of armored equipment. J. Acad. Armored Force Eng. 2012, 26, 24–28. [Google Scholar]
  6. Dorterler, M. A new genetic algorithm with agent-based crossover for the generalized assignment problem. Inf. Technol. Control 2019, 48, 389–400. [Google Scholar] [CrossRef] [Green Version]
  7. Fu, X.; Guo, Q.; Ma, Q. A multiple objective uncertainty assignment problem and the study of its algorithm. Oper. Res. Manag. Sci. 2013, 22, 34–38. [Google Scholar]
  8. Xiong, S.; Sun, Z.; Wang, X.; Ni, M. Generalized assignment model of talk allocation problem and its method. J. Mil. Commun. Technol. 2016, 37, 55–59. [Google Scholar] [CrossRef]
  9. Ding, S.; Zeng, X.J.; Zhang, H. Uncertain random simulation algorithm with application to bottleneck assignment problem. Soft Comput. 2019, 23, 10977–10982. [Google Scholar] [CrossRef] [Green Version]
  10. Ding, S.; Zeng, X.J. Uncertain random assignment problem. Appl. Math. Model. 2018, 56, 96–104. [Google Scholar] [CrossRef] [Green Version]
  11. Deb, K.; Jain, H. An evolutionary many-objective optimization algorithm using reference-point-based nondominated sorting approach, part I: Solving problems with box constraints. IEEE Trans. Evol. Comput. 2013, 18, 577–601. [Google Scholar] [CrossRef]
  12. Liu, Q.; Liu, X.; Wu, J.; Li, Y. An improved NSGA-III algorithm using genetic K-means clustering algorithm. IEEE Access 2019, 7, 185239–185249. [Google Scholar] [CrossRef]
  13. Chen, C.; Yuan, Y.; Yuan, X. An improved NSGA-III algorithm for reservoir flood control operation. Water Resour. Manag. 2017, 31, 4469–4483. [Google Scholar] [CrossRef]
  14. Bai, H.; Fan, T.; Niu, Y.; Cui, Z. Multi-UAV Cooperative Trajectory Planning Based on Many-Objective Evolutionary Algorithm. Complex Syst. Model. Simul. 2022, 2, 130–141. [Google Scholar] [CrossRef]
  15. Hou, Y.; Wu, N.Q.; Li, Z.W.; Zhang, Y.; Qu, T.; Zhu, Q.H. Many-objective optimization for scheduling of crude oil operations based on NSGA-III with consideration of energy efficiency. Swarm Evol. Comput. 2020, 57, 100714. [Google Scholar] [CrossRef]
  16. Ai, Y.; Du, M.; Pan, Z.; Li, G. The optimization of reactive power for distribution network with PV generation based on NSGA-III. CPSS Trans. Power Electron. Appl. 2021, 6, 193–200. [Google Scholar]
  17. Wang, S.; Tao, F.; Shi, Y. Optimization of location–routing problem for cold chain logistics considering carbon footprint. Int. J. Environ. Res. Public Health 2018, 15, 86. [Google Scholar] [CrossRef]
  18. Wang, X.; Wei, J.; Zhou, X.; Xia, Z.; Gu, X. Dual-Objective Collision-Free Path Optimization of Arc Welding Robot. IEEE Robot. Autom. Lett. 2021, 6, 6353–6360. [Google Scholar] [CrossRef]
  19. Wang, D.; Zhu, J.; Wei, X.; Cheng, T.C.E.; Yin, Y.; Wang, Y. Integrated production and multiple trips vehicle routing with time windows and uncertain travel times. Comput. Oper. Res. 2019, 103, 1–12. [Google Scholar] [CrossRef]
  20. Shao, C.; Wang, H.; Yu, M. Multi-Objective Optimization of Customer-Centered Intermodal Freight Routing Problem Based on the Combination of DRSA and NSGA-III. Sustainability 2022, 14, 2985. [Google Scholar] [CrossRef]
  21. Tang, Q. An Improved NSGA-III Algorithm to the Many-Objective Optimal Power Flow of Power System. Ph.D. Thesis, Xiamen University, Xiamen, China, 2017. [Google Scholar]
  22. Wang, S. I-NSGA-III Based Many-Objective Optimal Power Flow for Power System and Its Decision Support. Ph.D. Thesis, Xiamen University, Xiamen, China, 2019. [Google Scholar]
  23. Ghannami, A.; Li, J.; Hawbani, A.; Al-Dubai, A. Stratified opposition-based initialization for variable-length chromosome shortest path problem evolutionary algorithms. Expert Syst. Appl. 2021, 170, 114525. [Google Scholar] [CrossRef]
  24. Yu, B.J. A Genetic Algorithm Framework using Variable Length Chromosomes for Vehicle Maneuver Planning. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 2022. [Google Scholar]
  25. Asim, M.; Abd El-Latif, A.A. Intelligent computational method for multi-unmanned aerial vehicle-enabled autonomous mobile edge computing systems. ISA Trans. 2021, in press. [Google Scholar] [CrossRef] [PubMed]
  26. Rodríguez-Molina, A.; Solís-Romero, J.; Villarreal-Cervantes, M.G.; Serrano-Pérez, O.; Flores-Caballero, G. Path-planning for mobile robots using a novel variable-length differential evolution variant. Mathematics 2021, 9, 357. [Google Scholar] [CrossRef]
  27. Ni, J.; Wang, K.; Cao, Q.; Khan, Z. A memetic algorithm with variable length chromosome for robot path planning under dynamic environments. Int. J. Robot. Autom. 2017, 32. [Google Scholar] [CrossRef] [Green Version]
  28. Nazarahari, M.; Khanmirza, E.; Doostie, S. Multi-objective multi-robot path planning in continuous environment using an enhanced genetic algorithm. Expert Syst. Appl. 2019, 115, 106–120. [Google Scholar] [CrossRef]
  29. Cruz-Piris, L.; Marsa-Maestre, I.; Lopez-Carmona, M.A. A variable-length chromosome genetic algorithm to solve a road traffic coordination multipath problem. IEEE Access 2019, 7, 111968–111981. [Google Scholar] [CrossRef]
  30. Jiang, N.; Fan, W.; Xie, X.D.; Guo, F.; Li, E.; Zhao, S. Comparative study of NSGA-II and NSGA-III on multi-objective optimization of heat exchanger network. Chem. Ind. Eng. Prog. 2020, 7, 2534–2547. [Google Scholar]
  31. Qian, J.; Long, H.; Long, Y.; Zhao, C. Improved NSGA-Algorithm and BP Fuel-cost Prediction Network for Many-objective Optimal Power Flow Problems. IAENG Int. J. Appl. Math. 2021, 51, 1–14. [Google Scholar]
  32. Tong, Z.; Chen, Y.; Tong, S.; Yu, Y.; Li, J.; Hao, G. Multi-objective optimization design of low specific speed centrifugal pumps based on NSGA-III algorithm. China Mech. Eng. 2020, 31, 2239. [Google Scholar]
  33. Cui, Z.; Chang, Y.; Zhang, J.; Cai, X.; Zhang, W. Improved NSGA-III with selection-and-elimination operator. Swarm Evol. Comput. 2019, 49, 23–33. [Google Scholar] [CrossRef]
  34. Li, H.; Deb, K.; Zhang, Q.; Suganthan, P.N.; Chen, L. Comparison between MOEA/D and NSGA-III on a set of novel many and multi-objective benchmark problems with challenging difficulties. Swarm Evol. Comput. 2019, 46, 104–117. [Google Scholar] [CrossRef]
  35. Zhu, Y.; Liang, J.; Chen, J.; Zhong, M. An improved NSGA-III algorithm for feature selection used in intrusion detection. Knowl. Based Syst. 2017, 116, 74–85. [Google Scholar] [CrossRef]
  36. Vesikar, Y.; Deb, K.; Blank, J. Reference point based NSGA-III for preferred solutions. In Proceedings of the 2018 IEEE Symposium Series on Computational Intelligence (SSCI), Bangalore, India, 18–21 November 2018; pp. 1587–1594. [Google Scholar]
  37. Chaudhari, P.; Thakur, A.K.; Kumar, R.; Banerjee, N.; Kumar, A. Comparison of NSGA-III with NSGA-II for multi objective optimization of adiabatic styrene reactor. Mater. Today Proc. 2022, 57, 1509–1514. [Google Scholar] [CrossRef]
Figure 1. Rotation process.
Figure 1. Rotation process.
Applsci 12 10482 g001
Figure 2. Chromosome C 1 .
Figure 2. Chromosome C 1 .
Applsci 12 10482 g002
Figure 3. Chromosome C 2 .
Figure 3. Chromosome C 2 .
Applsci 12 10482 g003
Figure 4. Segmented chromosome C 1 and C 2 .
Figure 4. Segmented chromosome C 1 and C 2 .
Applsci 12 10482 g004
Figure 5. Heterotopic single-point crossover result.
Figure 5. Heterotopic single-point crossover result.
Applsci 12 10482 g005
Figure 6. Two-point crossover result.
Figure 6. Two-point crossover result.
Applsci 12 10482 g006
Figure 7. Mutation result.
Figure 7. Mutation result.
Applsci 12 10482 g007
Figure 8. Iterative results of I-NSGA-III-VLC.
Figure 8. Iterative results of I-NSGA-III-VLC.
Applsci 12 10482 g008
Figure 9. R C L of u 1 .
Figure 9. R C L of u 1 .
Applsci 12 10482 g009
Figure 10. R C L of u 2 .
Figure 10. R C L of u 2 .
Applsci 12 10482 g010
Table 1. Construction rotation decision matrix X of u 1 .
Table 1. Construction rotation decision matrix X of u 1 .
Equipment Number S 1 , 1 S 1 , 2 S 1 , N 1 1 S 1 , N 1
Rotation Time
1 x 1 , 1 x 1 , 2 x 1 , N 1 1 x 1 , N 1
2 x 2 , 1 x 2 , 2 x 2 , N 1 1 x 2 , N 1
c 2 x c 2 , 1 x c 2 , 2 x c 2 , N 1 1 x c 2 , N 1
c 1 x c 1 , 1 x c 1 , 2 x c 1 , N 1 1 x c 1 , N 1
Table 2. Construction task decision matrix Y of u 1 .
Table 2. Construction task decision matrix Y of u 1 .
Equipment Number S 1 , 1 S 1 , 2 S 1 , N 1 1 S 1 , N 1
Task
1 y 1 , 1 y 1 , 2 y 1 , N 1 1 y 1 , N 1
2 y 2 , 1 y 2 , 2 y 2 , N 1 1 y 2 , N 1
a U 1 1 y a U 1 1 , 1 y a U 1 1 , 2 y a U 1 1 , N 1 1 y a U 1 1 , N 1
a U 1 y a U 1 , 1 y a U 1 , 2 y a U 1 , N 1 1 y a U 1 , N 1
Table 3. Comparison of time and results of different methods.
Table 3. Comparison of time and results of different methods.
U 1 U 2 a MethodTime F 1 F 2 F 3 F 4
1231690.52 s84.57%90.34%91.11%100%
37.20 s84.66%90.34%91.11%100%
513.30 s88.65%91.14%92.50%100%
1631894.00 s87.16%90.45%91.11%100%
40.22 s87.16%90.34%92.13%100%
558.33 s88.58%93.23%93.47%100%
223385.45 s88.06%77.09%96.11%100%
50.83 s88.20%79.39%96.11%100%
516.19 s88.22%91.44%98.06%100%
2631511.76 s83.99%81.37%96.11%100%
102.30 s84.97%82.31%96.11%100%
531.87 s87.46%93.45%98.06%100%
6231426.69 s84.33%92.41%100%100%
89.17 s85.33%92.19%100%100%
531.22 s87.99%92.75%100%100%
6631473.19 s82.46%92.24%100%100%
99.73 s83.12%92.99%100%100%
560.79 s88.96%94.67%100%100%
Table 4. Equipment number of rotation in 3 months.
Table 4. Equipment number of rotation in 3 months.
Rotation StageEquipment Number
1‘F’, ‘U’
Table 5. The equipment number performing the task of u 1 .
Table 5. The equipment number performing the task of u 1 .
Task NumberEquipment NumberTask NumberEquipment Number
1‘J’, ‘M’, ‘B’, ‘S’10‘J’, ‘H’, ‘K’, ‘T’
2‘C’, ‘Q’, ‘O’, ‘R’11‘D’, ‘H’, ‘G’, ‘T’
3‘L’, ‘J’, ‘H’, ‘D’12‘J’, ‘S’, ‘N’, ‘Q’
4‘H’, ‘A’, ‘E’, ‘P’13‘A’, ‘P’, ‘O’, ‘G’
5‘A’, ‘K’, ‘P’, ‘D’14‘A’, ‘K’, ‘P’, ‘D’
6‘E’, ‘O’, ‘N’, ‘M’15‘I’, ‘B’, ‘N’, ‘K’
7‘P’, ‘D’, ‘I’, ‘L’16‘E’, ‘J’, ‘B’, ‘S’
8‘M’, ‘N’, ‘J’, ‘B’17‘H’, ‘N’, ‘A’, ‘K’
9‘T’, ‘J’, ‘C’, ‘M’18‘D’, ‘J’, ‘M’, ‘C’
Table 6. The equipment number performing the task of u 2 .
Table 6. The equipment number performing the task of u 2 .
Task NumberEquipment NumberTask NumberEquipment Number
1‘U’, ‘SS’, ‘JJ’, ‘KK’7‘DD’, ‘GG’, ‘CC’, ‘U’
2‘DD’, ‘FF’, ‘II’, ‘LL’8‘QQ’, ‘PP’, ‘CC’, ‘U’
3‘QQ’, ‘BB’, ‘AA’, ‘RR’9‘U’, ‘OO’, ‘HH’, ‘BB’
4‘LL’, ‘GG’, ‘KK’, ‘HH’10‘U’, ‘SS’, ‘JJ’, ‘KK’
5‘LL’, ‘GG’, ‘KK’, ‘HH’11‘GG’, ‘BB’, ‘AA’, ‘QQ’
6‘LL’, ‘MM’, ‘GG’, ‘U’12‘MM’, ‘RR’, ‘JJ’, ‘GG’
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Miao, Q.; Zhang, X.; Qin, P.; Liu, X. Modeling and I-NSGA-III-VLC Solution of Aircraft Equipment Rotation and Echelon Usage under Uncertainty. Appl. Sci. 2022, 12, 10482. https://doi.org/10.3390/app122010482

AMA Style

Miao Q, Zhang X, Qin P, Liu X. Modeling and I-NSGA-III-VLC Solution of Aircraft Equipment Rotation and Echelon Usage under Uncertainty. Applied Sciences. 2022; 12(20):10482. https://doi.org/10.3390/app122010482

Chicago/Turabian Style

Miao, Qinglin, Xiaofeng Zhang, Pisheng Qin, and Xianguang Liu. 2022. "Modeling and I-NSGA-III-VLC Solution of Aircraft Equipment Rotation and Echelon Usage under Uncertainty" Applied Sciences 12, no. 20: 10482. https://doi.org/10.3390/app122010482

APA Style

Miao, Q., Zhang, X., Qin, P., & Liu, X. (2022). Modeling and I-NSGA-III-VLC Solution of Aircraft Equipment Rotation and Echelon Usage under Uncertainty. Applied Sciences, 12(20), 10482. https://doi.org/10.3390/app122010482

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