A CEGA-Based Optimization Approach for Integrated Designing of a Unidirectional Guide-Path Network and Scheduling of AGVs

In the current industrial fields, automatic guided vehicles (AGVs) are widely employed to constitute the flexible manufacturing system (FMS), owing to their great advantages of routing flexibility and high efficiency. However, one main challenge lies in the coupling process of the design problem of the unidirectional guide-path network (UGN) and the task scheduling problem of AGVs. To reduce the complexity, most pertinent literatures only handle these problems one by one, based on the stepwise design methods, thereby neglecting the constraint conditions and the optimization objectives caused by the FMS environment. )e motivation of the paper is to bring the coupling factors into the integrated design and solution process. Firstly, an integrated design model of designing UGN and scheduling AGVs simultaneously is proposed, with the objective of minimizing the makespan (i.e., the maximum completion time of all handling tasks), in the consideration of the practical constraints, e.g., the job handling and processing sequence constraints and the AGV number constraint. Secondly, a dual-population collaborative evolutionary genetic algorithm (CEGA) is developed to solve the problems of designing and scheduling in a parallel way. )e solutions of the integrated model, i.e., the potential strongly connected UGN and the feasible processing and handling sequence, are, respectively, coded as two different subpopulations with independent and concurrent evolution processes. )e neighbourhood search operation, the niche technique, and the elitism strategy are combined to improve the convergence speed and maintain the population diversity. )e experimental results show that the integrated design model can formulate the problem more accurately, and the CEGA algorithm is computationally efficient with high solution quality.


Introduction
As a flexible and high-efficiency logistics system, an automatic-guided vehicle system (AGVS) has broad application prospects in fields such as port transportation, warehousing and distribution, and manufacturing [1]. e unidirectional guide-path network (UGN) design and the AGVS task scheduling [2][3][4] are the two main issues that should be addressed when the AGVS is adopted. Many design models and algorithms for these two problems have been proposed in the literature. For example, aiming at the UGN design problem, Kaspi et al. [5] proposed a branch and bound method to minimize the travel distance of loaded vehicles, which can obtain the optimal solution more stably. Since the time complexity of the algorithm increases exponentially with the scale of the network, the method is only suitable for small-scale design problems. To decrease the design time, some heuristic algorithms have been developed. For example, Ko and Egbelu [6] proposed a multistage design method that determines the direction of the path sequentially, according to the traffic flow, with the objective of minimizing the travel distance of loaded vehicles. However, this algorithm does not take into account the empty travel distance and the result may substantially differ from the actual situation. Seo et al. [7] solved the UGN design problem by using a tabu search (TS) algorithm but without the consideration of the empty travel distance. Besides, the algorithm has the comparatively strong dependence on the initial status, thereby being prone to fall into the local optimum. Guan et al. [8] proposed an improved electromagnetism-like mechanism algorithm to minimize both loaded and empty travel distance. Functionally similar approaches also include the genetic algorithm [9,10] and the met heuristics algorithm [11].
In recent years, many mathematical models and algorithms have been developed for the AGVs task scheduling problem in the manufacturing environment. For example, aiming at the logistics dynamic scheduling problem in an intelligent manufacturing workshop (IMW), Xu et al. [12] proposed an intelligent logistics scheduling model and response method with Automated-Guided Vehicles (AGVs) based on the mode of "requestscheduling-response." Umar et al. [13] investigated the simultaneous scheduling of machines and AGVs in an FMS environment. A mixed-integer nonlinear programming model was formulated firstly, and a hybrid multiobjective genetic algorithm was then proposed to optimize the makespan, AGV travel time, and penalty cost. Kumar et al. [14] addressed simultaneous scheduling of both machines and material handling systems with alternative machines for the makespan minimizing objective. A machine selection heuristic and a vehicle assignment heuristic were incorporated in the differential evolution approach to assign tasks to appropriate machines and vehicles and to minimize cycle time. Similar research works can be found in Nouri et al. [15], Ahmadi-Javid and Hooshangi-Tabrizi [16], Lin et al. [17], Lacomme et al. [18], and Zhang and Li et al. [19]. Ahmadi-Javid and Hooshangi-Tabrizi [20] addressed a ternary integration scheduling problem, which incorporated employee timetabling into the schedules of machines and transportation systems in a job shop environment. A timed coloured Petri net (TCPN) approach with two performance objectives was devised in [21] for the simultaneous scheduling problem of machines and AGVs. Moreover, two heuristics (A * and ALS [22]) were developed to partially explore the graph. Lyu et al. [23] focused on the machine and AGV scheduling problem in a flexible manufacturing system by simultaneously considering the optimal number of AGVs, the shortest transportation time, a path planning problem, and a conflict-free routing problem (CFRP). To study these problems simultaneously, they proposed a genetic algorithm combined with the Dijkstra algorithm that is based on a time window. Li et al. [24] proposed a mathematical model and an improved harmony search algorithm for material transfer in a realworld manufacturing system.
In fact, the design process of the UGN is potentially correlated with the task scheduling problem of AGVs in an FMS environment [25]. e mutual coupling relationship is mainly reflected in the following two aspects. On the one hand, the UGN determines the job transportation time among different workstations [26], while the transportation time is an important factor that must be considered in solving the AGVs task scheduling problem [27][28][29]. On the other hand, the AGVs task scheduling scheme influences the traffic flow [30] between each two workstations, which provides one main reference for the UGN design process. In this sense, these two problems are highly interacted with each other.
However, the mutual coupling relationship is often neglected in the pertinent literature, in order to decrease the complexity of the problem. e stepwise design method, which separates the UGN design from the AGV task scheduling, is mostly employed in the current design process. Firstly, the loaded traffic flow matrix among all workstations is generated on the basis of the job technological route information. Secondly, an optimal UGN is designed by optimizing only one objective, i.e., the loaded and empty travel distance. Other objectives and constraints in the FMS environment, e.g., the job handling and processing sequence constraints introduced by the job technological route, are neglected in this design stage. Finally, the simultaneous scheduling of machines and AGVs is addressed by taking into account the practical constraints and other objectives in the industrial field, e.g., minimizing the makespan. It is obvious that the stepwise design method reduces the complexity of the problem significantly. However, the solution quality of UGN cannot be always guaranteed due to the model infidelity without the coupling relationship.
e motivation behind this article lies in developing an integrated design model and a parallel optimization approach of the UGN design and AGV task scheduling, by considering the practical constraints imposed by the FMS environments. Our main contributions are twofold. On the one hand, an integrated design model of the UGN design and AGV scheduling in an FMS environment is constructed with the makespan minimization objective. e practical constraints in the industrial field, e.g., the job handling and processing sequence constraints, and the AGV number constraint, are considered in the model. On the other hand, a dual-population collaborative evolutionary genetic algorithm (CEGA) is proposed to solve the problems of designing and scheduling in a parallel way. Two kinds of chromosome are used, respectively, for coding the strongly connected UGN and the feasible processing and handling sequence as two different subpopulations. ey evolve independently and concurrently, only having an information exchange when an individual need to be evaluated. Moreover, the neighbourhood search operation, the niche technique, and the elite retention strategy are combined to improve the convergence speed and maintain the population diversity. e rest of this paper is organised as follows. e problem statement is described in detail in Section 2 and a new integrated model of simultaneously designing a UGN and scheduling AGVs in an FMS is formulated. In Section 3, a CEGA approach with specific evolutionary operators is proposed to solve the integrated model. In Section 4, computational experiments are carried out and experimental results are analysed and compared. Finally, conclusions and recommendations for future research are given in Section 5. Figure 1 shows a typical FMS layout and an undirected guided-path network of AGVS. In the figure, M 0 and M 9 represent the input workstation and output workstation, respectively. M 1 to M 8 are the processing workstations, and each workstation M m is composed of an input buffer, a processing machine, and an output buffer. e input workstation only has an output buffer, and the output workstation has only an input buffer. ere is an undirected guided-path network G � (V, E) in the middle, where V is the node set, which is composed of the loading points, unloading points, and other path intersection points. e edge (path segment) set is E, and the length of the edge e i,j is d i,j .

Problem Formulation
ere are N J kinds of jobs J � J j | 1 ≤ j ≤ N J to be processed in this workshop. As shown in Figure 2, each job J j is composed of n j processing tasks and n j + 1 handling tasks. All handling tasks are performed by N G AGVs. e processing workstation and the time of each processing task are known. e questions needed to be settled are as follows: the direction of each guided-path section, the execution AGV of each handling task, and the execution sequence of all processing and handling tasks. First, some assumptions are given as follows: (1) A processing machine can only perform one processing task at a time. Once a processing task is started, it cannot be interrupted. Each processing machine processes the jobs in the input buffer based on the first-come-first-served rule. (2) When the FMS starts to run, each job is located at the output buffer of the input workstation and is waiting to be transported to the first processing workstation according to its technological route. e job will enter the input buffer, when it arrives at the unloading point of the workstation. After the job completes a processing task by the processing machine, it will enter the output buffer of the corresponding workstation and wait to be transported to the next workstation. After completing all processing tasks, the job will leave the system via the output workstation.
(3) When the FMS starts to run, all AGVs are located at the loading point of the input workstation. ey will return to the loading point of the input workstation after completing all handling tasks. (4) e capacity of all buffers in the system is sufficiently large, so deadlock problems triggered by the insufficient buffer capacity are not considered. (5) An AGV can only carry one job at a time. e loading and unloading time are not taken into account. (6) An idle AGV will stay at the unloading points of its last handling task but does not affect the operation of other AGVs.
(7) Routing deadlocks caused by path conflicts between AGVs are not considered. N J is the number of job types in the system. Ω J � J j | j ∈ N and 1 ≤ j ≤ N J , the set of jobs. n j is the total number of processing tasks for the job J j .

Symbol Definitions
the set of all processing tasks for the job J j . Ω G � AGV k | k ∈ N and 1 ≤ k ≤ N G , the set of AGVs, and N G represents the number of AGVs in the system. b i,j is the direction of the path segment e i,j . If the value is 1, the direction is from node v i to v j ; the reverse direction occurs, if the value is 0. N P is the number of loading points. N D is the number of unloading points. N T is the number of handling tasks. X T i,j,k is the assignment relationship between the handling tasks T ij and the vehicle AGV k . e value of 1 means task T ij will be executed by AGV k ; otherwise, the value is 0. [D DP ] N D ×N P is the matrix of the shortest logistic distance among all unloading and loading points. d E w,u is the element located at the w − th row and u − th column in [D DP ] N D ×N P , representing the shortest logistic distance from the unloading point D w to the loading point P u . C k is the completion time for the AGV k as it returns to the loading point of the input workstation after finishing all handling tasks. S T i,j,k,m is the starting time of handling task T i,j , which is the m − th handling task of vehicle AGV k. E T i,(j− 1),k,a is the completion time of handling task T i,(j− 1) , which is the a − th handling task of the vehicle AGV k. S O i,j,w,n is the starting time of processing task O i,j , which is the n − th processing task of workstation M w. E O i,j,w,n is the completion time of processing task O i,j , which is the n − th processing task of workstation M w . v is the speed of the AGV.

Decision-Making Variables
(a) b i,j is the direction of all the path segments in the guide-path network. is variable is the decision variable of the traditional UGN design problem.
is the assignment matrix between all the AGVs and the handling tasks. is variable is the first decision variable of the AGVs task scheduling problem. (c) e feasible execution sequence of all the processing and handling tasks. is variable is the second decision variable of the AGVs task scheduling problem.

Objective Function.
e objective is to minimize the makespan, which is the maximum completion time of all AGVs return to the loading point of the input workstation after finishing all handling tasks. So, the objective function is

Decision Variable Constraints.
e decision variable is a 0-1 integer, that is, e direction of each path segment is unique:

Auxiliary Decision Variable
Constraints. e auxiliary decision variable is a 0-1 integer, that is, Each handling task can only be assigned to one AGV:

e UGN Must Be Strongly Connected.
Each node contains at least one outdegree (i.e., the number of edges starting from this node) and one indegree (i.e., the number of edges ending at this node), which is a necessary condition for a strongly connected UGN: ere is a route to link any pair of unloading point and loading point, so it can be written as follows: ere is a route to link any pair of loading point and unloading point, which can be expressed as follows: Handling task and its identifier T NJ,1 Processing task and its identifier O NJ,1

Processing and Handling Sequence Constraints.
Handling task T i,(j+1) can be executed only after the processing task O i,j is completed: Handling task T i,j can be started only after the AGV k has completed the previous handling task T c,d and arrived at the loading point P u of the task T i,j : Only after the job arrives at the current workstation, the current processing task O i,j can start, which can be expressed by Each processing machine can only perform one processing task at a time. So, only after the previous processing task O g,h is finished, the next processing task O i,j can start, which can be expressed by As said before, the UGN design problem is potentially correlated with the AGVs task scheduling problem. e mutual coupling relationship is mainly reflected in the following two aspects. On the one hand, it can be seen from constraints (10) and (11) that the shortest distances among all loading or unloading points (d E w,u and d L u,w ) are the important factors to determine the beginning time of all processing and handling tasks, while the shortest distance between all loading points depends on the decision variable of the traditional UGN design problem. On the other hand, only when the decision variables of the AGVs scheduling problem are known, the completion time C k can be determined by constraints (9) to (12), so as to evaluate the quality of the designed UGN.
As both the UGN design and the AGVs scheduling problem in the FMS environment are nondeterministic polynomial time-(NP-) hard problems [31,32], the integration of these two problems will become a more complex NP-hard problem. As a high-performance multipopulation heuristic intelligent search algorithm, the collaborative evolutionary genetic algorithm (CEGA) [33] divides the system decision variables into several variable subgroups. Each variable subgroup can select appropriate encoding scheme to obtain a subpopulation. Moreover, appropriate genetic operators can be designed for each subpopulation according to its encoding characteristics. e CEGA is very suitable for solving the problem of multivariable combinatorial optimization [34,35]. Based on the characteristics of the abovementioned mathematical model, a dual-population CEGA involving the collaborative evolution of the UGN subpopulation and the task sequence subpopulation is proposed.

Optimization Algorithm
Drawing upon the idea of collaborative evolution, two subpopulations are adopted to represent the UGN and the feasible task sequence, respectively. e UGN subpopulation is labelled as P1, in which each chromosome represents a strongly connected UGN. e task sequence subpopulation is labelled as P2, in which each chromosome corresponds to a feasible handling and processing task sequence. In the evolutionary process, the two subpopulations independently and concurrently evolve. Information exchange between the subpopulations occurs only during the process of fitness calculation. e algorithm mainly includes chromosome encoding, fitness calculation, subpopulations initialization, and the genetic evolutionary operations. e detailed methods are listed as follows.

Encoding of the UGN.
To reduce the complexity of the problem, the guide-path network is preprocessed before encoding. e node with only two path segments connected has special significance. Taking the path network shown in Figure 1, for example, edge (1, 2) should have the same direction as edge (2,3), in order to make node 2 satisfy the constraint shown in formula (6). erefore, in the design process, edges (24, 1), (1, 2), (2, 3), (3,4), (4,5), and (5, 6) can be merged into one combined edge (24,6). According to this preprocessing principle, the guide-path network shown in Figure 1 has nine combined edges. All combined edges and serial numbers are as shown in Table 1 Figure 2, a corresponding handling task occurs before each processing task. erefore, the task sequence is coded based on the processing task sequence. Each gene in the task sequence chromosome denotes a processing task, which is expressed by the job type number. e length of code is equal to the total number of processing tasks of all jobs put together. For instance, when there are three jobs need to be processed and each job requires three processing tasks. If a feasible processing task sequence is

Fitness Calculation.
Since the individual in each subpopulation represents only partial optimization problem, it cannot be evaluated independently. erefore, if an individual in a subpopulation need to be evaluated, it has to combine with an individual in the other subpopulation to form a complete chromosome. e fitness calculation method for a complete chromosome is described as follows.

Mathematical Problems in Engineering 5
A complete chromosome is composed of two parts. Take the complete chromosome of [111111100][ 4 1 3 4 2 1 3 2  2 1 2 3 4 1 4 1 3 2 4 3 1 2 3 4 ] as an example. e first half is the UGN code. Referring to Table 1, the corresponding UGN is constructed, as shown in Figure 3. And the shortest logistic distance between all the loading and unloading points can be obtained by using the Floyd algorithm [36], as shown in Tables 2 and 3. If the AGV travelling speed is given, the logistics transportation times between all the loading and unloading points can be calculated. erefore, the UGN determines the logistic time between workstations. e last half part of the complete chromosome is the task sequence code. First, the task sequence code is converted to a processing task sequence. Second, according to the fitness calculation flow shown in Figure 4, the machine tool and the AGV are scheduled according to the task sequence to calculate the fitness. Assuming that the technological route information of the job set is provided, as shown in Table 4, the number of AGVs in the system is 4, the AGV speed is 1 m/s. e Gantt chart obtained by the fitness calculation flow is shown in Figure 5, and the makespan is 1,260 seconds.

Subpopulation Initialization.
e methods of the initial subpopulation generation for the two subpopulations are discussed as follows.

UGN Subpopulation.
e size of the UGN subpopulation is N P1 . e procedure for generating the initial UGN subpopulation is as follows. First, generate a random UGN chromosome code and assign each bit of the code a random value of 0 or 1. en, decode this code and check whether the corresponding UGN is strongly connected; if yes, keep the code; otherwise, generate a new code. Repeat the above steps until N P1 feasible codes are obtained. ese feasible codes constitute the initial UGN subpopulation G 0 P1 .

Task Sequence
Subpopulation. e size of the task sequence subpopulation is N P2 . N P2 distinct chromosome individuals are generated randomly, and these individuals constitute the initial task sequence subpopulation G 0 P2 .

Evolutionary Operations for the UGN Subpopulation.
e evolutionary algorithm flow for the UGN subpopulation is shown in Figure 6.
e key steps include selection, crossover, mutation, individual evaluation, neighbourhood search, and niche elimination operations. e specific operations are described as follows.

UGN Selection.
Firstly, all individuals in the parent subpopulation G i P1 are sorted from large to small according to their fitness. en, two different individuals are selected to enter the crossover operation by using the roulette selection method. e probability of an individual being selected is proportional to its fitness.

UGN Crossover.
is operator is implemented using the double-point crossover method. A crossover operation example is shown in Figure 7. Firstly, two crossover positions are randomly generated. en, the genes between the two positions of the parent individual 1 are replaced by the genes at the corresponding positions of the parent individual 2. e changed parent individual 1 becomes a new child individual.

UGN Mutation Operation.
e two-point mutation is employed, and a mutation operation example is shown in Figure 8. Firstly, two mutation positions are randomly generated. en, the genes at the two positions are reversed. Finally, the strongly connectivity of the UGN that corresponds to the code is verified. If the UGN is strongly connected, evaluate the individual; otherwise, repeat the steps of selection, crossover, and mutation operations again.

UGN Individual Evaluation.
Since the UGN individual represents only partial optimization problem, it cannot be evaluated directly. erefore, in order to evaluate an UGN individual, several individuals are randomly selected from the task sequence subpopulation by using the roulette method. Multiple complete chromosomes are formed by combining these task sequence individuals with the UGN chromosome code. e fitness of each complete chromosome is separately calculated, and the greatest fitness is selected as the individual fitness for the UGN chromosome. To improve the convergence speed of the algorithm, the complete chromosome with the highest fitness will be used by the neighbourhood search operation.

Neighborhood Search.
As already mentioned, a complete chromosome is composed of the UGN code and the task sequence code. A neighbourhood solution will be obtained by reversing a certain bit of the UGN code, while other bits are unchanged. For an original complete chromosome with a UGN code length of N B , there are N B neighbourhood solutions. First, the strong connectivity of the UGN that corresponds to each neighbourhood solution is verified. If the UGN is strongly connected, calculate its fitness value; otherwise, set the fitness value to 0. en, find  (18,12) the neighbourhood solution with the highest fitness. If this neighbourhood solution is better than the original individual, use it to replace the original one.

Niche Technique and Elitism Strategy.
To measure the genetic differences between two individuals, the Hamming distance between two individual chromosomes is defined: In the formula, b i k is the value of the k − th position in the UGN individual code B i .
To avoid the early maturity in the evolution process of the subpopulation, it is necessary to maintain genetic diversity among the individuals in the subpopulation. Firstly, the N P1 offspring individuals are merged with all N P1 parent individuals to obtain 2N P1 individuals as the candidate subpopulation.
en, calculate the Hamming distance H(B i , B j ) between any two chromosomes. If the calculated H(B i , B j ) value is smaller than H min (which means they are similar chromosomes), multiply the fitness of the poor individual by an elimination factor η < 1. Finally, these 2N P1 individuals in the candidate subpopulation are sorted according to the adjusted fitness, and the optimal N P1 individuals are selected to the next generation.

Evolutionary Operations for the Task Sequence
Subpopulation.
e flowchart of the evolutionary algorithm for the task sequence subpopulation is shown in Figure 9.
e key steps include the selection, crossover, mutation, and individual evaluation operations. e details are provided as follows.

Task Sequence Selection.
Using the roulette selection method, two different individuals from G i P2 are selected to enter the crossover operation step.

Task Sequence Crossover.
is operator is implemented by two-point crossover. An example of crossover operation is shown in Figure 10. Assuming that the two parent individuals are Parent1 and Parent2, the crossover operation is described as follows: (i) Step 1: randomly select two crossover positions Pos1 and Pos2 (ii) Step 2: translate the parent individual code into the processing task sequence     Step 3: rearrange the processing tasks of Parent1 between the two crossover positions according to their order in Parent2 (iv) Step 4: copy the rearranged Parent1 to the child individual (v) Step 5: convert the child processing task sequence to the task sequence code.

Task Sequence
Mutation. An example of the task sequence mutation operation is shown in Figure 11. Randomly select two crossover positions Pos1 and Pos2. If the job types are different at these two positions, then the codes of the two positions are exchanged. Whether the mutated individuals have appeared in the parent subpopulation needs to be verified. If this individual has appeared in the parent population, then the selection, crossover, and mutation operations will be repeated until a different individual is generated.

Task Sequence Evaluation.
Firstly, the roulette method is used to select multiple individuals from the UGN subpopulation. Secondly, the task sequence individual is merged with these UGN individuals to form multiple complete chromosomes. Finally, calculate the fitness of each complete chromosome and select the largest fitness as the evaluated individual fitness.

Task Sequence Elitism Strategy.
New N P2 offspring individuals are merged with N P2 parent individuals. ese 2N P2 individuals constitute the candidate population. N P2 optimal individuals of the candidate population are selected to constitute a new generation task sequence subpopulation. Figure 1. e AGV speed is set at 1 m/s. e technological information of the job set is shown in Table 4. e integrated design method proposed in this paper is compared with the traditional stepwise design methods, which can be implemented as follows: firstly, the technological route information is converted into the loaded flow matrix among all loading and unloading points (Table 5). Secondly, with the objective of minimizing both loaded and empty travel distance, the method proposed by GUAN and DAI [8] is used to design the UGN. e results obtained by the method is the optimal UGN code [111111100], which corresponds to the UGN shown in Figure 3. e best empty flow matrix among all unloading and loading points is shown in Table 6. e total loaded and empty travel distance End Fetch the next processing task of the task sequence chromosome.

UGN Design Experiment. e layout of an example FMS is shown in
Initialize system status, input the task sequence chromosome.
Determine the workstation where the job is located at the present moment.
Find which AGV can reach the loading point of the workstation early and dispatch this AGV to the loading point.
AGV waits for the completion of the previous processing task of the job and loads the job as soon as the previous processing task is completed.
Transport the job to the next workstation according to the technological route. The processing task will start once the machine is idle.
Update the status of the job, corresponding workstation, and AGV.
Dispatch all AGVs to the loading point of the input workstation.
Output the objective function (makespan). is J � J L + J E � 1920 + 800 � 2720, which is calculated by formula (1) in reference [8]. Finally, with the makespan minimization objective, the simultaneous scheduling of machines and AGVs are addressed under this UGN by taking into account the constraints of the FMS application environment, e.g., the job handling and processing sequence constraint, as shown by formula (9)- (12), and the number of AGVs constraint, as shown in formula (5). e optimal results of the traditional methods and our approach for an AGVS with different numbers of AGVs and jobs are compared in Table 7.

Start
As shown in Table 7, the number of AGVs and the task sequence constraints introduced by the job technological route are all neglected by the traditional design methods based on the loaded flow matrix. Since the loaded flow matrix that corresponds to the technological routes is unchanged, the optimal UGN results obtained by the traditional method remain the same in all cases. However, when taking into account the constraints of the FMS application environment and the makespan minimization objective, the integrated design model proposed in this paper can adjust the optimal UGN according to the increase number of AGVs, and it can obtain a better solution than the traditional design model. In addition, as shown in Table 7, when both the number of jobs and AGVs is 4, the traditional stepwise design method gets a solution with the makespan of 1260 and the actual travel distance of 2720. e corresponding Gantt chart is shown in Figure 5. On the other hand, the integrated design method can get a better solution with the actual travel distance of 3490. is indicates that the total loaded and empty travel distance is not directly related to the makespan. However, the traditional UGN design model can only optimize the UGN by minimizing the total loaded and empty travel distance and reducing the UGN quality of the FMS application environment inevitably. is experiment proves the superiority of the integrated design model.  e CEGA proposed in this paper. e parameters of the CEGA are set as follows. e population size of the UGN subpopulation N P1 is set to be 20. e population size of task sequence subpopulation N P2 is set at 100. Because the CEGA adopts elitist retention strategy to ensure the stability of the algorithm, so, we set the crossover rate P c � 1 and the mutation rate P m � 1. (iii) TGA: among the CEGA, if without incorporating the niche technology and the neighbourhood search strategy, it becomes the traditional collaborative evolution genetic algorithm (TGA). (iv) NGA: adding the niche strategy to the TGA will transform it into the niche collaborative evolution genetic algorithm (NGA). e parameters of the TGA and the NGA are consistent with that of the CEGA. All programs are written in Plant Simulation 15.0 from Siemens Ltd. e program interface is shown in Figure 12. e running environment is Windows 10 with Intel Xeon E3-1535M V6 3.1 GHz, 32 GB RAM. Each algorithm is executed ten times using the same calculation example (both the number of jobs and AGVs is 4). e optimization process obtained by each algorithm is recorded every 0.4 minutes. e average makespan and the standard deviation performance of all the algorithms are shown in Figures 13 and 14.

Optimization Algorithm
As shown in Figure 13, the TS algorithm has a fast convergence speed in the early stage, when the search time is less than 0.8 minutes, it can get the local optimal solution as the best solution, but due to the lack of global search ability, the optimal result of the TS mainly depends on the initial

Replace old population with new population
The The (i + 1)-th generation population G P1 i+1 Figure 6: Collaborative evolutionary algorithm for the UGN subpopulation.   Generate initial task sequence population G 0 p2 i+1 Figure 9: Collaborative evolutionary algorithm for the task sequence subpopulation.  random solution, which has greater randomness. As shown in Figure 14, the standard deviation of the TS is always greater than 180, which indicates that the stability of the algorithm is poor. Different from TS algorithm, the standard deviation of the three genetic algorithms decreases with the progress of the search process. In the later stage of the search, the standard deviation of the three genetic algorithms can be reduced to less than 75, which shows that the stability of genetic algorithm is much higher than the TS. Among the three genetic algorithms, the TGA has the good initial convergence effect. As the diversity among populations cannot be guaranteed, the TGA is easily trapped into the   local optimum. Although the addition of niche technology causes a relatively poor initial convergence of the NGA, it can maintain the diversity of the populations. Hence, the later search performance of the NGA is better than that of the TGA. e CEGA combines both benefits of the niche technique and the neighbourhood search strategy, which is computationally efficient with the high solution quality. As shown in Table 7, when both the number of jobs and AGVs is four, the optimal makespan is 1,200. e proposed CEGA can stably obtain the global optimal solution after 3.6 minutes during each of ten runs, exhibiting the excellent overall performance.

Conclusions
A mutual coupling relationship exists between the UGN design process and the AGV task scheduling problem in the FMS environments. However, the coupling factors are often   neglected by using the stepwise design method in the pertinent literature, in order to reduce the complexity of the problem at the cost of the solution quality of UGN. e motivation of the paper is to bring the coupling factors into the integrated design process. On the one hand, an integrated design model of designing a UGN and scheduling AGVs simultaneously is proposed, with the objective of minimizing the makespan, in the consideration of the practical constraints. On the other hand, a dual-population collaborative evolutionary genetic algorithm (CEGA) is developed to solve the problems of designing and scheduling in a parallel way. e experimental results show the satisfactory efficiency and performances of our model and approach. It is noteworthy that the AGVs in this research work still pertain to the ordinary single-load type. However, multiple-load AGVs are becoming increasingly popular in modern factories, owing to their less traffic congestion, smaller fleet size, and increased system throughput [37,38]. Hence, the integrated design and solution process oriented to multipleload AGVs will attract great interest in the future research work.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.