Inspection Allocation Optimization with Resource Constraints Based on Modified NSGA-II in Flexible Manufacturing Systems

With the development of smart manufacturing, quality has become an indispensable issue in themanufacturing process. Although there is increasing publication about inspection allocation problems, inspection allocation optimization research considering resource capability is scarce. *is paper focuses on the inspection allocation problem with resource constraints in the flexible manufacturing system. Combined with the inspection resource capability model, a bi-objective model is developed to minimize the cost and balance loads of the inspection station. Amodified NSGA-II algorithm with adaptive mutation operators is suggested to deal with the proposed model. Finally, a simulation experiment is conducted to test the performance of the modified algorithm and the results demonstrate that modified NSGA-II can obtain acceptable inspection solutions.


Introduction
In manufacturing systems, the quality characteristics (QC) of the manufactured product may inevitably fail to meet the requirement of the design because of machines, human error, etc. erefore, inspection is a necessary part of the product manufacturing process [1]. Inspection is an organized assessment process, and it can be arranged after the processing procedure or even after the completion of the product. e quality and the condition of determining material of products will be inspected to know whether the product is suitable for use and the product meets the regulated requirement or not [2]. In reality, if there are too many rounds of inspection in the process of product manufacturing, it will raise the inspection cost and prolong the production cycle. us, many scholars conducted inspection process planning and inspection allocation research to arrange the inspection flow properly and dispatch the inspection resource, thus optimizing the cost of product inspection [3,4]. Lindsay and Bishop [5] firstly studied the inspection allocation (IA) problem and proposed a perfect inspection model with the lowest cost for serial production. A dynamic planning model was used to get inspection solutions. Hurst et al. [6] built an imperfect inspection model which, for the first time, added the first and the second type of inspection errors into the inspection allocation model. Shiau et al. [7] built a model which integrated processing planning and inspection planning problems under limited resource. e model took a series of manufacturing systematic factors like manufacturing capability, inspection capability, and the cost of product quality into consideration to minimize the total cost. Korytkowski [8] proposed a multiproduct product inspection model, in which the inspection station was regarded as a G/G/C queuing model. e genetic algorithm was adapted to deal with the model to shorten the waiting time of inspection. Mohammadi et al. [1] brought up an optimizing framework for the inspection process and the framework aimed at two types of inspection, conformity inspection and monitoring inspection, balancing the cost and the quality.
Even many scholars researched the inspection allocation problem. However, the problems that current research mainly focused on are the inspection location and inspection process planning. Many scholars regarded the inspection stations as the same inspection unit, which is an ideal processing way to some extent. In the actual condition, different inspection stations own specific inspection capabilities, such as precision, process time. erefore, in this paper, two subproblems are taken into consideration. e first one is the inspection location problem and the second one is the inspection resource allocation problem based on the inspection capability model of Shiau [9]. Following the inspection process planning definition of Narahari et al. [10], the paper summarized the inspection allocation problem with resource constraints as follows: (i) e location of inspection in the product manufacturing cycle. (ii) Which inspection station will be used to inspect the QC?
For the type of manufacturing system, people paid more attention to the research on the IA problem in the series manufacturing system. ere are few papers concerning the IA problem in the Flexible Manufacturing System (FMS). In the FMS, the product flow and inspection flow of a production line have many types of parts. e inspection stations shall inspect different types of components that contain separate product QCs, and each of the components owns a specific inspection plan. erefore, the allocation of inspection resource becomes very complicated in a complete view.
On the other hand, plenty of scholars took the inspection cost, quality improvement, and uncertain factors into consideration in the inspection allocation problem. e inspection unit is a crucial resource, and its working load is worthy of concern. However, it seldom was mentioned by scholars. When the inspection station is put into use and each working unit gets a relatively balanced working quantity, the overall inspection efficiency will be improved. erefore, in the paper, a multi-objective inspection allocation with resource constraints (IARC) model is built in the FMS to trade-off manufacturing cost and inspection station load. is research looks forward to reducing the planning difficulty in the inspection allocation and giving support to decision-making. e methodology to solve the inspection planning problem varies, such as dynamic programming, integrated programming [11]. With the increase of system size, the amount of computation of previous methods increased and these methods became impractical [12]. In recent years, more and more scholars have paid attention to heuristic algorithms and metaheuristic algorithms to solve the inspection planning problem, which belongs to an NP-hard problem [13]. Especially, genetic algorithm (GA) is the most suitable method to deal with the inspection allocation problem as the codes express the inspection policy directly. Hsin Rau et al. [12] proposed a GA approach for solving the inspection allocation problem, and they also analysed the similarity of the codes used in the chromosome of the GA approach and the representation of the inspection allocation policy for workstations in the production system. Due to the complexity of inspection planning problems, some scholars treated them as multi-objective optimization problems and proposed corresponding solutions.
Mohammadi et al. [14] presented a bi-objective, mixed-integer linear programming (BOMILP) model to tackle the inherent uncertainty of the concerned parameters. A differential evolution algorithm combined with the Taguchi and Monte Carlo methods was proposed to solve the model. Besides, simulation methods were also used in this problem [15]. However, simulation just investigated a limited number of production scenarios for optimization.
In this paper, the proposed IARC model belongs to a multi-object optimization problem (MOP), in which objects contract with each other. us, finding the optimal solutions for every object is impossible. Multi-objective heuristic algorithms are practical to solve the MOP problem; many methods were developed in recent years, such as NSGA-II, NASG-III, MOEA/D, Monarch butterfly optimization, and differential evolution algorithm, quantum-inspired differential evolution [16][17][18][19][20]. For the inspection allocation problem, the methods for a single-objective problem are abundant, such as the dynamic programming algorithm, genetic algorithm [12,[21][22][23][24]. For multi-objective problems, there are a fewer number of methods than methods in the single-objective problem. NSGA-II is one of the most popular methods, and it is used to solve the inspection planning problem and obtain an acceptable result [11]. In this paper, the modified NSGA-II algorithm is used for solving the IARC model. e following are the main contributions of this paper: (1) Inspection-capability-based IARC is the first model which allocates inspection location and inspection station in the FMS. Besides, the IARC is a bi-objective optimization model which tries to trade-off the cost and machine load. (2) A modified algorithm is proposed with three adaptive mutation operators for the IARC model. e simulation experiment validates the performance of the proposed algorithm. e remainder of the paper is organized as follows. In the next Section, problem description and mathematical modeling are presented. In Section 3, a modified NSGA-II with three kinds of mutation operators is described. In Section 4, simulation and analysis are conducted to verify the effectiveness of the model and algorithm. Finally, conclusions are stated in Section 5.

Problem Description.
e FMS production line includes many types of parts. Here, Ω � {1, . . ., P} is used to express parts. A part p∈Ω has K p quality characteristics QC p,q , and each QC p,q owns the target value M p,q , upper deviation USL p,q , and the lower deviation LSL p,q . In the system, there are machine tools M � {1, . . ., N} and inspection stations Λ � {1, . . ., I}. According to the production technique, the work in process (WIP) gets QC p,q after being processed in the machine n. If the next procedure is not an inspection for part p, it will be transferred to the next working procedure. If the component should be inspected, then it will be conveyed to the inspection station i∈Λ for inspection. e products in the system share the same manufacturing resource M and inspection resources Λ. e inspection strategy for manufactured products is as follows: qualified (to the next procedure), reworked, and discarded. When the inspection value LSL p,q < v p,q < USL p,q , the product is qualified and shall be put into the next procedure. When v p,q > USL p,q , the product shall be reworked. When v p,q < LSL p,q , the product is not qualified and will be discarded. e workstation n and inspection station i have their specific resource capability. erefore, the first and second types of errors exist. e actual value x of processing QC p,q of processing station n is assumed to be subjected to the normal distribution f(n) ∼ N(u n , σ n ). Similarly, the measured value y of inspection QC p,q of the inspection station is assumed to be subjected to the normal distribution φ i ∼ N(u i ,σ i ). e inspection results after processing of QC could be divided into 7 types. For specific content, please refer to Table 1, and for the calculation regulation, please refer to Shiau [25].
Pursuing lower costs is a vital issue in manufacturing [25]. Costs in the production process include manufacturing cost, inspection cost, scrap cost, and loss cost of defective product sold to customers. Inspection plans have an indirect effect on costs as inspection will affect the products' flow in manufacturing. If the inspection procedure increases, the inspection cost will increase, but the subsequent manufacturing cost will be reduced as the defective goods will be detected out in advance in manufacturing [14]. At the same time, defective goods will be prevented from flowing into the market. e decrease in the number of defective goods flowing into the market will reduce the damage caused by defective goods. erefore, the total cost is the first objective function obj1 of the IARC model. In the FMS, the inspection unit shall process many kinds of inspection tasks. e method of assigning tasks to inspection stations is worth studying. When each inspection station is assigned a relatively balanced inspection task, the overall inspection efficiency will be improved. To prevent the large task gap between inspection stations, the maximum load of the inspection unit is the second objective function obj 2 .

Model Hypothesis
(1) In manufacturing, each machine tool processes only one product QC p,q (2) In an inspection procedure, each inspection station inspects only one QC p,q (3) Suppose that the whole system is working correctly, and the processing flow is known After inspecting QC p,q , the processing strategy includes rework, discard, and entering the next step. If the QC p,q does not need an inspection, the product can enter the next procedure directly. While the QC p,q needs an inspection, the qualified one can enter the next working procedure. Even the second type of error occurs, namely, the defective goods will also enter the next step by mistake. Table 1: e inspection probability based on the resource capability.

Parameters of the IARC Model
e probability of product p in the q working procedure could flow into the next working procedure: (1) e probability of product p in the q working procedure being discarded: (2) e probability of product p in the q working procedure being reworked: e probability of the product p entering the q working procedure:

Cost Model.
Following the possible results of the product after each working procedure, we can calculate the inspection cost, manufacturing cost, reworking cost, discarding cost, and damage cost caused by selling the defective goods to the customer. Based on the above probabilities, we can get the related costs as follows: e manufacturing cost of product p: e inspection cost of product p: e reworking cost of product p: e damage caused by discarding the product p: e damage caused by selling the defective goods p to the customer: By integrating all the cost factors, and we can get the total cost of product p.

Time Model.
Following the distribution condition of the inspection station, the total working time of inspection station i is:

Mathematic Model.
Aiming at the optimization problem of inspection allocation under the FMS, the minimum total cost and minimum biggest inspection station load were taken as targets to build a multi-objective MIP model: Objective function: TMC p , TIC p , TRC p , TSC p , TCSC p , t i,p,q > 0, ∀p ∈ Ω, ∀q � 1 . . . K p , ∀i ∈∧.
Formula (12) expresses how to minimize the total cost of all the products. Formula (13) expresses how to minimize the maximum inspection time among all inspection stations. Formula (14) shows the relation between the maximum load and the working time of inspection stations. Formula (15) represents the relation between the total cost and related costs. Formula (16) expresses that each QC shall be arranged to one inspection station. In formula (17), w p,q,i , h p,q are the variates of 0-1. Formula (18) means that the value of costs and time are positive numbers.

The Modified NSGA-II Algorithm
e IARC model mentioned in this paper includes two fitness function formulas, (12) and (13). To address this model, NSGA-II is adopted to calculate the Pareto optimal solution of targets and balance these two optimizing targets. Besides, some necessary modification is done with the algorithm to adapt the IARC model and enhance the performance. A four-layer coding scheme for chromosomes is put forward, and it can express the inspection locations and inspection resource assignment scheme directly. In the mutation scheme, three types of mutation operators are raised to enhance neighborhood searching capability. e multipoint crossover method is adopted to the NSGA-II algorithm. e flow chart of the modified NSGA-II is shown in Figure 1.
e crowding distance calculation and nondominated sorting method are adopted in this paper referring to Deb [26].

Coding.
e coding scheme adopted four-layer chromosomes to express the inspection location and inspection resource assignment scheme, and it is shown in Figure 2. e first layer is the production layer, and the value of each unit in the layer is the serial number of the products. e second layer is the QC layer, which shows the serial number of QCs. e third layer shows the inspection scheme. When the QC needs inspection according to the plan, the value is 1, otherwise it is 0. e fourth layer is the resource assignment layer, which shows the serial number of inspection stations. When the value is 0, the associated QC does not need an inspection. Each layer of chromosomes is a task unit O, which contains information on inspection strategy and inspection resource. For instance: O (1,1,1,1) means that the QC 1,1 shall be inspected in inspection station 1 after processing. It is worth noting that layer one and layer two are static in the evolution process.

Population Initialization.
is paper proposed three ways for initial solution generation: (1) 20% Random Generation Method (RGM), (2) 40% Minimum Load Selection Method (MLSM), and (3) 40% Earliest Stage Assignment Method (ESAM). e ESAM was proposed by Shiau [25]. In ESAM, the earlier the manufacturing procedure is, the higher weight the QC will be allocated. ere is a higher probability that this QC will be allocated an inspection station. e MLSM selects the inspection station with the minimum weight among all the available resources to inspect the QC. e detailed steps of the initial solution generation are as follows: Step 1: All QCs are traversed in order. When the QC exists, skip to step 2. When all the QCs have been traversed, end.
Step 2: e 0-1 variable v 1 is generated randomly. When v 1 � 1, the QC needs an inspection, and then skip to step 3 to allocate the inspection station. When v 1 � 0, the QC does not need an inspection and returns to step 1.

Multipoint Crossover Method.
e multipoint crossover method is adopted in this algorithm. P1 and p2 are paternal chromosomes, and c1 and c2 are, respectively, the offspring chromosomes generating from p1 and p2. Two cross positions are selected randomly and gene fragments of p1 and p2 are exchanged. Finally, the new generations c1 and c2 are created. e crossover method is shown in Figure 3.

Adaptive Mutation Operators.
e most common mutation method is single point variation, which has a shortage.
e generated population has a low diversity, and it further makes the solutions fall into the local optimal solution prematurely. Besides, multilayer encoding is characteristic of deep dependency between different layers of genes. If the method of single point variation is adopted, the infeasible offspring may be created, as the other layer of genes may not change with the mutated gene.
In recent years, neighborhood search strategy has been successfully applied to combinatorial problems with good performance [27]. In this section, three different variation strategies are proposed by combining the neighborhood search strategy with variation strategy to improve the diversity of the population. e mutation operator selects the task to be optimized based on the current population. Once the task subset is selected, the corresponding variable is deemed to be optimized. Mutation can be considered as a neighborhood search operation on an existing solution structure. e specific contents include: (1) Random inspection task change operator (RIC) : RIC randomly selects a QC and changes the inspection policy of the QC. If QC does not need inspection before, this QC is reassigned as an inspection procedure. en, an available inspection station will be allocated to the QC. If the QC needs an inspection, the inspection policy of the QC is cancelled. (2) Random equipment replacement operator (RER) : is operator can randomly select an inspection procedure and randomly reassign a piece of available equipment for this inspection procedure. e mutation operators are designed based on human experience on the inspection allocation problem. e three mutation strategies differ from each other but do not have the same importance. In the process of population evolution,

Problem input and parameters initialization
Reach the temination condition?
Multi-point crossover Crossover  an adaptive selection strategy is proposed. For each generation of chromosome mutation, roulette is used to select the mutation operator according to the weight of the mutation operator w θ , θ∈Θ. e initial value w θ is 1. With the increase of generation, the mutation operator weights w θ are updated according to the improvement of chromosome fitness S gen after mutation. e computation rules are as follows: S gen � obj1 gen − obj1 gen−1 obj1 gen−1 + obj2 gen − obj2 gen−1 obj2 gen−1 .

Modified Elite Selection Strategy.
Selection is an essential operation, which transmits excellent and adaptable parents to the next generation and maintains the quantity of the population unchanged. e individuals in the population are divided into different frontiers of inferior quality according to govern relations by using the nondominated sorting method. e elite selection strategy always retains the nondominant solution in the lower frontier to the next generation. In this case, the majority of chromosomes come from level one of the nondominant frontier of the last generation. us, the optimal result always entraps into the local optimum solution. e elite selection strategy is improved by adjusting the probability of the lower-level solution to the next generation. e specific steps are as follows: Step 1: e chromosomes of the parent and offspring are combined P ai � Ch i ∪F ai . e total quantity of chromosomes is |F ai |. e individuals are divided into several frontiers, which are expressed as F t � {F 1 , . . ., F i }. e number of optimal solutions at inferior t is |F t |.
Step 2: e next generation is F ai + 1 � ∅, and the initial Pareto frontier is t � 1.
Step 3: e quantity of individuals that are selected to the next generation at inferior t is n i � (1-0.1 t)|F t |. If n i + |F ai + 1| ≤ |Ch i |, execute step 4. Otherwise, execute step 5.
Step 4: e individual at the present level |F t | and the parent generation P ai + 1 are combined, P ai � Ch i ∪F ai + 1 t � t+1. Skip to the last step.

Simulation and Analysis
In this Section, a simulation experiment is conducted to verify the correctness of the proposed inspection planning model and the feasibility of the modified NSGA-II algorithm. e experiment bases on a flexible production line of an aircraft manufacturer and the parameters are estimated following Shiau [25].
e proposed inspection planning model belongs to a multi-objective programming problem whose outcomes are normally difficult to evaluate. us, a quantified indicator called "hyper-volume indicator" is adapted to evaluate the performance of NSGA-II, modified NSGA-II. Besides, the outcome of a specific example is listed individually to analyze the inspection solutions and related cost.
To verify the effectiveness of the proposed algorithm, NSGA-II and modified NSGA-II are used to calculate the examples. Matlab R2020a was used to implement the program, which was run on the Intel (R) Core (TM) i5-7200 2.50 GHz computer under the Windows 7 operating system. e parameters of both algorithms were set as follows: popsize � 50, max-_gen � 100, p c � 0.75. e mutation probability for NSGA-II is only 0.25. To evaluate the performance of the two algorithms, the following evaluations of the algorithm are adopted.

Metric and Objective Normalization.
e hyper-volume (HV) [28,29] is a qualitative metric that can evaluate the convergence and diversity of multi-objective problems synchronously [30,31]. is metric is used by calculating the volume of the objective space that is covered by individuals of a nondominated solution set. e value of HV depends on the number of objective functions, the number of individuals in the nondominated solutions set, and reference points. e equation is where ND is a nondominated solution set D 1 , D 2 , . . . , D Normalization of the objective function values is required before calculating the hyper-volume to obtain a nondimensionalised optimization problem [32]. A simple normalization method is adapted to replace the objective [33]. rough normalization, the values of objective functions range from [0, 1].
where z * i is the max point in the objective space, z nad i is the nadir point in the objective space.

Result Analysis.
is section consists of two parts: the outcome analysis and performance analysis. e outcome analysis aims to analyse the costs and the load of inspection tasks by comparing the result of a specific inspection task leveraging by the IARC model and the NSGA-II algorithm. e objective of performance analysis is to compare the performance of NSGA-II and modified NSGA-II in the inspection planning problem by using the HV metric.

Outcome Analysis.
e inspection instance I4_A was selected as an example to carry out the outcome analysis. ree methods, namely, Modified NSGA-II, NSGA-II, and all inspection strategy, were adapted to deal with this inspection task. All inspection strategy means that all the QCs of 4 parts from class A will be inspected by 3 inspection stations, and the inspection stations were allocated to these inspection tasks randomly. Finally, the result is obtained, and the selected scheme, the costs, the objective value, and the HV are shown in Table 4. e obj1 value is the sum of MC, IC, DC, RC, and CNC, and it means that all factors affect the result collectively. It is clear that the proportion of MC, DC, and CNC is relatively large among all factors by three methods, but the value of MC changes slightly among the three plans. us, it has a slight influence on obj1 as the manufacturing operations are inevitable. For the comparison of the outcome of NSGA-II and modified NSGA-II, there are more nondominated solutions by using modified NSGA-II. By comparing the outcome of all inspection strategies with heuristic algorithms, the results from the optimization algorithm are distinctly superior to the solutions from all inspection strategies. Especially, the cost of inspection and discard is visibly higher than the other two methods, because more inspection operations were set in the plan and more defective products were taken out along the process. e value of CNC in all inspection strategy is lower than others, as the inspection is strict and defective products can seldom flow into the market.
For the comparison of the two optimal algorithms, please refer to Figure 4, which shows the distribution of nondominated solutions in the objective space of two algorithms. e horizontal axis represents the total cost, and the vertical axis represents the load balance. e modified NSGA-II can generate more nondominated solutions compared with NSGA-II, and the solutions of modified NSGA-II are dominated to the nonmodified one. To further compare the performance of the two algorithms, performance analysis is done in the next section.

Performance Analysis.
Performance analysis is conducted by dealing with the 17 instances provided in Table 3 ten times using two algorithms, NSGA-II and the modified NSGA-II. e HV is used to assess the results, and statistical analysis is carried out according to the results. e statistic variables are shown in Table 5.
For the mean value of 10 times runs, both the Modified NSGA-II and NSGA-II can see a slight growth in mean value  Table 3: e instances of the experiment.   I1_A I1_B NA I2_A I2_B I2_AB I4_A I4_B I4_AB  I6_A I6_B I6_AB I8_A I8_B I8_AB I10_A I10_B I10_AB   8 Mathematical Problems in Engineering as the number of products involved in the experiment increased. When there is only one product in the experiment, the difference between the mean values of the two algorithms is tiny. e gap between the two results becomes large when the number of products involved in the experiment increased. e largest difference between the two algorithms is 7.73% in I10_AB and the average difference of the mean value of all instances in 10 runs is 3.36%. It means that the modified NSGA-II is excellent when dealing with complex inspection allocation task and the modified NSGA-II outperforms NSGA-II in dealing with the IARC model. e variance of the modified NSGA-II is smaller than the corresponding value of NSGA-II in most cases, which means that the solutions of the previous algorithm converge evidently than the latter one. e average number of solutions in the nondominated solutions set generated by modified NSGA-II is far more than the other one and it means that the previous one performs well in terms of variety. For the performance of algorithms in different tasks, the modified NSGA-II performs better when dealing with hybrid class AB than class A and class B. Figure 5 shows a box plot of the HV metric in I4_A. e computation complexity of the algorithms is compared, and the average time of ten runs is listed in Table 5. e average runtime of modified NSGA-II is higher than NSGA-II by around 8%. e reason is that the adaptive mutation operators need more extra time, which is in an acceptable range. Figure 6 shows a box plot of runtime in example I4_A, and there is no significant difference between the two algorithms in terms of runtime.
It is worth mentioning that the performance of different mutation operators is different when dealing with the IARC model. e performance is evaluated by using the average changes of weight value for mutation operators, which is shown in Table 6. e initial values of all operators are the same.
us, the larger the value of average changes, the greater the contribution of the method to objective optimization. e MLR is the best mutation operator in this experiment. RIC is the second one, and the RER has the worst performance. e MLR focuses on the load balance to optimize the obj 2 . In the simulation, the value of obj 2 improvement is large, which can be concluded by observing the obj 2 in Table 4. erefore, the MLR has the best performance in this experiment.

Discussion
rough the outcome analysis, the wrong idea that all inspection strategy is always the most economically beneficial solution was corrected, but all inspection strategy is widely adopted in some factories. By using the quantitative methods, the lowest-cost solution could be found. However, the formulation of the method often considers more than a single factor in the real application. In this study, the inspection resource and the cost are both important for stakeholders. e IARC model can generate compromise solutions to make a trade-off between the cost and the load balance. Decision-makers can choose one among the multiple solutions generated by modified NSGA-II algorithms.
e experiment results show that the modified NSGA-II algorithm has a good performance in terms of diversity and convergence compared with NSGA-II in the inspection allocation problem. e modified algorithm can obtain more solutions that perform better in an objective space by using the HV metrics. It proves the effectiveness of adaptive mutation operators, which can improve the neighborhood search capability of the modified NSGA-II.
Regarding future work, more influence factors will be included in the IARC model as the objectives, such as energy efficiency, and a multi-objective optimization model will be   built to find the optimum inspection solutions. Besides, a production system normally has many activities, such as production, inspection, and maintenance, and these activities need to be scheduled to improve the efficiency of the whole system. us, joint scheduling of production, inspection, and maintenance activities need further research. In terms of algorithm, the proposed modified algorithm could be employed in similar problems, and it is interesting to investigate the performance in other MOP problems. e experiment result shows that the proposed IARC model has a certain guidance meaning for the realization of the intelligent, automation inspection planning systems. Although the model still has a long way to go from the practical application as it has some hypotheses, it is an interesting direction that integrates the IARC model into the IPS, which is commercially used in factories. By combining the inspection resources in the industrial field, the model can   provide real-time inspection planning solutions, which can support the decision-making process.

Conclusions
is paper researched inspection allocation problems considering inspection resources. e main contributions are as follows: (1) Based on the characteristics of the FMS, an IARC model is proposed to analyse the possibility and probability value of the product in the product flow by using the inspection capability model. e IARC mixed-integer programming model is built to tradeoff the total cost, the load of the inspection station. (2) To solve the multi-objective optimization model, the NSGA-II algorithm is improved. A four-layer coding scheme, the improved elite strategy, and mutation operators are integrated into NSGA-II. e neighborhood research capability of the methods is enhanced so that the solutions do not fall into the first layer of the Pareto solution prematurely. Finally, the simulation experiment was conducted. e result verified the effectiveness of the IARC model, the performance of the modified NSGA-II in terms of diversity and convergence.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest
e authors declare no conflicts of interest.