A Multilayer Genetic Algorithm for Automated Guided Vehicles and Dual Automated Yard Cranes Coordinated Scheduling

At present, a lot of studies on automatic terminal scheduling are aimed at the shortest operating time. An effective way to reduce the operating time is to increase the amount of operating equipment. However, people often ignore the additional costs and energy consumption caused by increasing the amount of equipment. .is paper comprehensively considers the two aspects of the equipment operation time and equipment quantity matching. With the minimum total energy consumption of the operating equipment as the objective function, a cooperative scheduling model of Automated Guided Vehicles (AGVs) and dual Automated Yard Cranes (AYCs) is established. In the modelling process, we also considered the interference problem between dual Automated Yard Cranes (AYCs). In order to solve this complex model, this paper designs an improved multilayer genetic algorithm. Finally, the calculation results from CPLEX and a multilayer genetic algorithm are compared, and the effectiveness of the model and algorithm is proved by experiments. In addition, at the same time, it is proved that it is necessary to consider the interference problem of dual Automated Yard Cranes (AYCs), and the optimal quantity matching scheme for the equipment and the optimal temporary storage location is given.


Introduction and Literature Review
With the deepening of economic globalization, the status of automated container terminals is becoming increasingly prominent. Reasonable scheduling of loading and unloading machinery and equipment has become the key to improving the efficiency of terminal operations. In addition, with the increasingly serious problem of climate change, low carbon and energy-saving have become urgent issues in container terminal production. e main operating equipment includes dual Automated Yard Cranes (AYCs), and Automated Guided Vehicles (AGVs) are also energy-intensive equipment in the automated container terminal. Sim [1] pointed out that the total carbon emissions of the container terminal were comprised of 37.34% from the container loading and unloading process, 1.04% from the container transportation process, and 9.92% from the container receiving and delivery process. erefore, it is of great significance for container terminals to realize green and low-carbon development by jointly optimizing the configuration and scheduling of double AYCs and AGVs to reduce the energy consumption of terminal loading and unloading operations. Figure 1 shows the operating system of an automated container terminal. Taking the unloading operation as an example, the containers are transported to the container yard by Automated Guided Vehicles (AGVs) from vessel. e dual Automated Yard Cranes (AYCs) are responsible for stocking and retrieving containers. e loading process is reversed. e dual AYCs often interfere with each other when completing their loading and unloading tasks. erefore, the key point of this article is to solve the scheduling problems in the environment of double AYC interference: (1) AGV resource configuration; (2) the energy consumption of AGVs; (3) the operation path of AYCs; and (4) the energy consumption of AYCs.
At present, a great deal of research has been done on the energy consumption of container terminal loading and unloading operations and the configuration and scheduling of AYCs and AGVs. He et al. [2] pointed out that the YC scheduling problem was firstly converted into a vehicle routing problem with soft time windows (VRPSTW). is problem was formulated as a mixed-integer programming (MIP) model, whose two objectives minimize the total completion delay of all task groups and the total energy consumption of all YCs. Liu and Ge [3] proposed a convex mathematical programming model for the QC assignment problem, in which the queuing theory is used to model the queuing behavior of automatic guided vehicles (AGVs). e objective of the proposed model was to minimize CO 2 emission during an unloading process of containers from QCs to AGVs by optimizing the number of QCs. Huang Xiaobo et al. [4] considered three sources of carbon emissions during the moving process, the loading and unloading process, and the preparation process, based on the feature of RTGs that they cannot cross each other. A mathematical route programming model for RTGs is developed to minimize the carbon emissions. A path strategy for RTGs is designed to address the computational complexity of a mixed-integer programming model. A simulated annealing algorithm is applied to find the near-optimal solution. Xin et al. [5] provided a methodology for determining the trajectory of interacting machines that transport containers between the quayside area and the stacking area in an automated container terminal. Simulation studies illustrate that energy consumption of container handling can indeed be reduced by the proposed methodology. Yang et al. [6] proposed a mixed-integer programming (MIP) model for the integrated scheduling issue of AGVs and Rail-Mounted Gantry Cranes (RMGs) so as to minimize the makespan of unloading operations with the task allocation constraints of AGVs and RMGs. S. Hu and Z. Hu [7] considered the cooperative scheduling mechanisms among the three types of devices, a full freedom optimization problem for an integrated quay crane, yard crane, and yard truck is studied and a mixed-integer programming model is built. Via simulation, the sequencing and operating times of tasks with different combinations of quay cranes, yard cranes, and yard trucks are analysed. e study provides a basic model to coordinate the allocation and dispatch of critical operating resources in container terminals. Chang and Zhu [8] developed an integrated scheduling model to improve the coordination of different types of equipment in a container terminal.
is model considers not only loading and unloading simultaneously but also the gantry crane interference and safety margin and the gantry crane and yard crane travel times. Moreover, the buffer area and congestion between inner trucks are also considered. en, an improved multilayer genetic algorithm is proposed to solve the problem. Le and Bo [9] considered the synchronization between YCs and YTs. Based on this information, the constraints that exist in actual operations are first considered, such as the noncrossing constraint when several YCs share a bidirectional lane. ere are other constraints, such as fixed YC separation distances and job-precedence constraints. A mathematical model is formulated to describe the problem, and the objective is to minimize the makespan. e interference in among multiple cranes is also widely studied. Liang et al. [10] used two operation modes, which are the relay mode and the mixed mode, to study the effects of different operation modes on the efficiency of dual-ARMG operations. Considering the relaying problem of the dual-ARMG task, the safety distance between the temporary buffer area and the dual-ARMG was set. Huang and Han [11] proposed a mixed-integer programming model for the collaborative scheduling problem of dual automatic stacking cranes (AYCs) with both of the operations (storage and retrieval) on a single container block of an automated container terminals. A genetic algorithm based on job sequence coding was designed and multiple sets of practical examples were solved by using CPLEX and a genetic algorithm and by considering the interference between dual AYCs. Zhan et al. [12] discussed the load scheduling problem of multiple yard cranes. A mathematical model, which considers the interference between adjacent yard cranes, is provided that use a time-space network to formulate the problem and a two-stage hybrid algorithm composed of a greedy algorithm and dynamic programming is developed to solve the proposed model. Liang et al. [13] studied the rail-mounted gantry cranes scheduling problem (RMGCSP) that considers the interference and safety distance between these cranes. A firefly algorithm (FA) is proposed. Park et al. [14] proposed heuristic-based and local-search-based real-time scheduling methods for dual rail-mounted gantry (RMG) cranes working in a block at an automated container terminal. Yang et al. [15] proposed a multiagent model to solve the problem in which moving the yard cranes is always hampered by each other.
ere is a great deal of research on the improvement of the algorithm when studying the scheduling problem. Gharehgozli et al. [16] studied an operational problem arising at a container terminal, consisting of scheduling a yard crane to carry out a set of container storage and retrieval requests in a single container block. e problem is modelled as a continuous-time integer programming model and the complexity is proven. ey used the intrinsic properties of the problem to propose a two-phase solution method to optimally solve the problem. Lee et al. [17] proposed a novel approach that integrates these two problems as one. e objective is to minimize the weighted sum of the total delay of requests and the total travel time of yard trucks. Due to the intractability of the proposed problem, a hybrid insertion algorithm is designed to find effective problem solutions. Ng [18] examined the problem of scheduling multiple yard cranes to perform a given set of jobs with different ready times in a yard zone with only one bidirectional traveling lane. ey developed a dynamic programming-based heuristic to solve the scheduling problem and an algorithm to find the lower bounds for benchmarking the schedules found by the heuristic. Li et al. [19] developed an efficient model for YC scheduling by taking into account realistic operational constraints such as intercrane interference, fixed YC separation distances, and simultaneous container storage/retrievals. ey show how the model can be solved quickly using heuristics and a rolling-horizon algorithm, yielding close-to-optimal solutions in seconds. Ulrich and Schneider [20] described an approach for scheduling triple crossover stacking cranes in an automated container storage block with asynchronous handover at the transfer areas at both block front ends. Wang and Xiao [21] applied ant colony optimization for efficient crane scheduling and reducing some of the terminal costs. In 2011, they applied differential evolution optimization to achieve the same goal. Jin et al. [22] presented a scheduling model for the subsystem to minimize the overall operational time of a gantry crane. A hybrid algorithm called the NSGA based on the nearest neighbour strategy (NS) and genetic algorithm (GA) was developed. Chen et al. [23] proposed a mixed-integer programming model to minimize the total operating time and determined the crane movement and work status (moving or handling) by considering the time and space synchronization constraints of dual crane parallel operations. A genetic algorithm was developed to quickly solve the near-optimal solutions for large-scale problems. Pei and Chang [24] proposed a mixed-integer programming model to minimize the total operating time and the distance between dual ARMGs on the basis of considering the coordination of the ARMGs and the capacity of the seaside handover point. In view of the complexity of solving the problem, an innovative strategy was used that combines the genetic algorithm with simulation in which the simulation model is designed to optimize the problem, and a genetic algorithm was employed to optimize the initial generated solutions. Zeng et al. [25] developed a scheduling optimization model to consider the characteristics of no waiting, block and batch processing. e lower bound of the model was formulated and algorithms based on the taboo search and heuristics dispatching rule were designed. Feng et al. [26] developed an integrated optimization model of a handling and retrieving sequence for inbound containers to optimize the retrieving sequence yard crane assignment and handling scheme. To solve the model, a dynamic programming-based heuristic is designed.
Compared with other literatures, the innovations of this paper are shown in Table 1. e rest of this paper is organized as follows. In Section 2, we describe the technical aspects of the problem and present the mathematical model. In Section 3, the solution method is developed. Section 4 presents the computational experiments, and Section 5 contains the conclusions.

Problem Description.
e automated yard crane used in this paper is a relay-type double-yard crane. Each block is equipped with two cranes, which are, respectively, recorded as AYC 1 and AYC 2 .
e two AYCs cooperate with each other to complete tasks. e specific process is shown in Figure 2. First, the AGV transports a container to the handover area of the yard from the vessel. en, AYC 1 unloads the container from the AGV and places it in the temporary storage area of the block. At the same time, AYC 2 unloads the container in the temporary storage area and moves it to the designated position until all unloading tasks are completed.

Mathematical Problems in Engineering
In addition to considering the cooperative scheduling problem of the AGVs and AYCs, the interference between the two AYCs should also be considered in this paper. In order to effectively solve the problem of coscheduling interference in dual AYC access blocks, the following task splitting rules are adopted. e temporary storage area is set in the block, and the yard is divided into two blocks. Each AYC can only operate in its own half-block, and the tasks across the half-blocks need to be completed by the double AYC using the relay area. e interference between the two AYCs only occurs in the relay area. e import blocks and export blocks designed in this paper are separated. ere are only two situations in which interference may occur. Assume that the import container (the container is transported from vessels to yard by AGV) is represented by i and the export container (the container is transported from yard to vessels by AGV) is represented by j.
ese two situations include the following: (1) import container i + 1 is unloaded by AYC 1 and import container i is loaded by AYC 2 , and the times of the two tasks are the same; (2) export container j + 1 is unloaded by AYC 2 and export container j is loaded by AYC 1 , and the times of the two tasks are the same. e interference process from the two AYCs is shown in Table 2.
In this paper, the cooperative scheduling of AGVs and dual AYCs is studied in two stages. In the first stage, AGV is configured according to the number and location of containers to be loaded and unloaded, and the initial scheduling scheme for AGV is obtained. In the second stage, the mutual interference in double AYCs is considered to obtain the walking track of double AYCs. And the initial scheduling scheme for AGV is adjusted to obtain the optimal scheduling scheme. It can realize the minimum energy consumption of AGVs and AYCs.

Mathematical Model.
We make the following assumptions based on problems. (1) e AGVs can transport all containers (import containers and export containers) and the AGVs can serve any AYC. (2) In the process of import and export operations, the AYCs choose the container group with small turnovers and short empty driving times as the priority. (3) e dual AYCs work at the same rate and energy consumption. (4) All AGVs have the same performance such as driving speed and energy consumption. (5) Uncertain factors such as path conflict in AGV transportation are not considered. e first stage is to assign loading and unloading tasks to each AGV and arrange the sequence of operations, so as to ensure that the loading and unloading from each AYC 1 is not delayed. e optimization goal is to minimize the total energy consumption of the AGV. Table 3 shows the relevant symbols.
{ }, 1 if the AGV v executes container j ′ after container j, 0 otherwise: Equation (1) shows that the energy consumption of the AGV transportation process is minimized. It specifically refers to the minimization of the total energy consumption including the load and no-load energy consumption of the AGV and the waiting energy consumption of the AGV: Equations (2) and (3) indicate that a container i(j) is transported by only one AGV v: ET kj � T vj , Equations (4) and (5) represent the time window constraint that the AGV v should meet in the import (export) container handover area: Table 3: e description of symbols.

Symbol
Description i e number of import containers (i � 1, 2, . . . , n) j e number of export containers (j � 1, 2, . . . , m) v e number of AGVs (v � 1, 2, . . . , V) k e number of AYC 1 s in the import (export) container yard (k � 1, 2, . . . , K) T vi e moment when the AGV v unloads the container i in the import container handover area T vj e moment when the AGV v loads the container j in the export container handover area ET ki e earliest moment when the AYC k unloads the container i from the AGV v in the import container handover area LT ki e latest moment when the AYC k unloads the container i from the AGV v in the import container handover area ET kj e earliest moment when the AYC k loads the container j from the AGV v in the export container handover area LT kj e latest moment when the AYC k loads the container j from the AGV v in the export container handover area t k e average time required for the AYC k to complete a container load/unload operation t vi e time for the AGV v to transport the import container i from the vessel to handover area t vj e time for the AGV v to transport the export container j from the handover area to vessel t ii′ e time for the AGV v to transport empty to the next container i′ after delivery of the container i t jj′ e time for the AGV v to transport empty to the next container j′ after delivery of the container j t ij e time for the AGV v to transport empty to the container j after delivery of the container i w vi e waiting time of the AGV v for the AYC k to unload the container i w vj e waiting time of the AGV v for the AYC k to load the container j C 1 e energy consumption for each AGV to transport containers per unit time C 2 e energy consumption for each AGV to transport empty per unit time C 3 Each AGV's waiting energy consumption per unit time at the handover area Mathematical Problems in Engineering Equations (6) and (7) represent the waiting time of the AGV v in the import and export container handover area: T vi ≥ 0, Equation (8) represents the utilization rate of the AGV and equation (9) represents the value range of the parameters.
In the second stage, AYC 1 and AYC 2 alternately complete the task of loading and unloading containers. Considering the interference constraints on the dual AYCs operation in the temporary storage area, the goal is to minimize the total energy consumption. Table 4 shows parameters and variables. e decision variables are as follows: x ki ∈ 0, 1 { }, 1 if the container i is assigned dual AYC k, 0 otherwise. y kj ∈ 0, 1 { }, 1 if the container j is assigned dual AYC k, 0 otherwise. z kii′ ∈ 0, 1 { }, 1 if the container i is assigned dual AYC k after container i ′ , 0 otherwise. u kjj′ ∈ 0, 1 { }, 1 if the container j is assigned dual AYC k after container j ′ , 0 otherwise: Equation (10) shows that the energy consumption of the AYC operation process is minimized. It specifically refers to the minimization of the total energy consumption including the load and unload energy consumption of the AYC, the transportation energy consumption, and the waiting energy consumption of the AYC: K k�1 y kj � 1, j � 1, 2, . . . , m.
Equations (11) and (12) indicate that each container can only be assigned to one dual AYC: Equations (13) and (14) indicate that there is a task in the front and at back of each container: Since each container needs to be operated alternately by AYC 1 and AYC 2 , equations (15) and (16) mean that each container is split into a container pair: T fi + t ki′ ≤ T fi′ , Equations (17)- (20) represent the interference between the dual AYCs.

Solution Algorithms
According to the characteristics of the two-way cooperation model of the AGVs and the AYCs, and considering the interference constraint on the dual AYCs, a multilayer genetic hybrid algorithm is designed to solve the problem. e genetic algorithm used in this paper is divided into two layers. e first layer is used to determine the order for AGV transporting containers, and the second layer is used to determine the optimal AYCs completion order according to the order of the first layer. erefore, the generation of chromosomes in the second layer is limited by the chromosomes in the first layer.

Encoding and Decoding.
In order to deal with the cooperative scheduling problem of AGVs, AYC 1 and AYC 2 , the multilayer chromosome and integer code methods are adopted. First, the import and export containers' blocks of the automated containers yard are coded. Each gene is composed of several subgenes, and the import and export containers corresponded to subgenes. It is assumed that each container block comes equipped with dual AYCs. e number of AGVs is 5. Second, it is assumed that the initial number of tasks of each import container block is m and the initial number of tasks of each export container block is n. Import container tasks i and i + 1 represent that the subtasks are split from the same initial task i, and i is the precedence task. Export container tasks j and j + 1 represent that the subtasks are split from the same initial task j, and j is the precedence task. Suppose that the number of the tasks in the import container blocks is 10. According to the task splitting rules, 20 subtasks are obtained, which are jointly completed by AYC 1 and AYC 2 . Figure 2 is the schematic diagram of the chromosome coding. e individual structure of the first layer represents an alternative solution of AGV operation sequence. Correspondingly, the individual structure of the second layer represents an alternative solution of AYC. e initial solution of the first layer individual is generated under the constraint of priority order, and the second layer individual   (j � 1, 2, . . . , m) T si e moment of the dual AYC starts to operate the import container i T fi e moment of the dual AYC finish the import container i T sj e moment of the dual AYC starts to operate the export container j T fj e moment of the dual AYC finish the export container j t ki e time for the dual AYC k to transport the import container i t kj e time for the dual AYC k to transport the export container j w ki e waiting time of the dual AYC k for the AGV v to transport the container i w kj e waiting time of the dual AYC k for the AGV v to transport the container j w kk′ e time that AYC 1 and AYC 2 wait for each other C 4 e energy consumption for each AYC to transport containers per unit time C 5 e energy consumption for each AYC to load (unload) containers per unit time C 6 Each AYC's waiting energy consumption per unit time h A safe distance between the dual AYCs is generated randomly under the constraint of the first layer individual. Take 5 AGVs, 1 dual AYC, and 10 container tasks as an example. As shown in Figure 3, the first group is the code of containers (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20). A container task was split into a pair (12, 34, 56, . . .). e second group is the first level individual. 0 is used to separate the operating sequence groups of different AGVs. e operation sequence of the first AGV is (1,11), the second AGV is (3,19), the third AGV is (5), the fourth AGV is (7,13,17), and the fifth AGV is (9,15). e third group is the second layer individual. e assignment sequence of AYC 1 is (1,11,3,19,5,7,13,17,9,15), and the assignment sequence of AYC 2 is (2,12,4,20,6,8,14,18,10,16).

Crossover and Mutation.
e first level individuals adopted order crossover (OX). e order crossover is used as the crossover operator for the container chromosome (the first vector in a chromosome). Since the aforementioned initialization procedure and order crossover ensure that the first gene of each offspring is the first container of a route, we can determine the specified AGV for the first gene. After the crossover operation of the container chromosome, the AGV chromosome (the second vector in a chromosome) is determined by the following procedure. First, we randomly select a substring in the parent chromosome. e second step is to generate an offspring and ensure that the selected genes in the offspring are in the same positions as the parents. e third step is to find the positions of the genes selected in the first step in the other parent and then place the remaining genes into the offspring generated in the previous step in order. e order crossover is shown in Figure 4.
For the second level of individuals, the structure is more complicated because it is related to the structure of the first level. e ordinary order crossover method is not suitable for the first level individuals. e crossover method of the second level requires to merge two-parent generations to produce an offspring, as shown in Figure 5.
In addition, two genetic positions are randomly selected as crossing points in the two-parent chromosomes of the AYC segments, and the corresponding matching segments are exchanged to obtain the offspring chromosomes. e crossover process is shown in Figure 6. e main task of the mutation is to promote the diversity of population and avoid the GA's premature convergence to local optimal solutions. In this paper, a swap mutationoperator is used for the container chromosome. e mutation method used in this article randomly selects two elements and then swaps their positions. e operation sequence of AGV aftermutation is (1,11), (3,17), (5), (7,13,19), (9,15). e assignment sequence of AYC 1 after mutation is (1,11,3,17,5,7,13,19,9,15), and the assignment sequence of AYC 2 after mutation is (2,12,4,18,6,8,14,20,10,16). e mutation results of the first level are shown in Figure 7, and the mutation results of the second level are shown in Figure 8.

Computational Experiments
We take the Qingdao port automated container terminal as an example. It is assumed that there are four import container blocks and four export container blocks. e length of the container blocks is 20 bays, and the bay 0 is the handover area between the AGVs and the AYCs. e number of containers is 8-400. e number of AGVs is 4-20. Each container blocks is equipped with dual AYCs. Tables 5 and 6 represent the basic parameters of the AGVs and the AYCs, respectively.
In order to make the experimental results have better convergence, this paper debugs the basic parameters of the genetic algorithm through many experiments and finally obtains the parameters shown in Table 7.
At first, an example is given to illustrate the result of the algorithm. Figure 9 is the distribution diagram of containers to be loaded and unloaded in a container block. e number of AGVs is 4, the number of dual AYCs is 1, and the number of containers is 10.
As shown in Figure 10, it can be seen which bay of the dual AYCs at any time. e operation track of AYC 1 and AYC 2 intersects in Figure 10(b). It indicates that there is a crossover between dual AYCs, and it is also not allowed. is fully proves the importance of considering interference issues.
In addition, in order to verify the effectiveness of the algorithm, small-scale and large-scale examples are used to estimate the results. e number of containers for smallscale problems is 8-40. e CPLEX and GA are used to solve 8 Mathematical Problems in Engineering         the problem. In order to reduce the random error in the GA solution, the average operation time T − and the average optimal value C − are recorded. e difference in the optimal value DOV indicates the difference between the optimal value obtained by the GA and the optimal value obtained by CPLEX. e calculation results are shown in Table 8.
According to examples 1-3, when the number of containers is small, CPLEX can quickly calculate the results. However, as the number of containers increases, the solution time of CPLEX becomes increasingly longer, as shown in examples 6-8. When the number of containers reaches 50, CPLEX cannot obtain the optimal solution in an acceptable time, as shown in examples 9-10. As the number of containers increases, the solution time of the GA does not change dramatically and is almost stable within 4-10 s. ere is little difference between the optimal solution obtained by the GA and the optimal solution obtained by CPLEX. e biggest difference is in example 8, which is only 4.8%. is shows that the GA is very effective at solving the dual scheduling problem.
When the example size is large, it is difficult for CPLEX to obtain the optimal solution in an acceptable time. However, the GA has unique advantages in solving largescale examples, and the results are shown in Table 9.
According to the results of the experiments, the GA can obtain the optimal solution in a short time. It can be seen that when the number of containers increases and the numbers of AGVs and AYCs remain unchanged, the energy consumption will inevitably increase, as shown in examples 11 and 13. With the number of containers and AYCs remaining unchanged, the energy consumption will increase as the number of AGVs increases, as shown in examples 13,14,19, and 20. is shows that the impact of the number of AGVs on the energy consumption is greater than that of the operation time. According to examples 15-18, as the number of AGVs increases, the number of AYCs decreases, and the energy consumption increases significantly. To sum up, the impact of the number of AGVs on the energy consumption is greater than that of the AYCs. erefore, the scheduling of the AGVs should be considered first, and then the AYCs must be coordinated to improve the overall operating efficiency and reduce the energy consumption. e interference in the dual AYCs will increase the waiting times of AYCs and AGVs, thus increasing the energy consumption. erefore, the following section will analyze the interference problem of dual AYCs.
A temporary storage area is set in the container block, as shown in Figure 3, so that the AYCs can complete the  loading and unloading from containers through the temporary storage area, which can reduce the frequency of interference among the dual AYCs. erefore, determining the location of the temporary storage area is the key to optimizing the interference problem in the dual AYCs. We take experiment 5 with the smallest difference in Table 7 as an example and use the CPLEX and the GA to solve the interference problem. e experimental results for when the location of the temporary storage area is changed are as shown in Figure 11. e locations of the temporary storage area obtained by the two algorithms are the same, and there are 23 bays. is result further proves the effectiveness of the GA in solving such problems.

Conclusion and Further Research
is paper studies the cooperative scheduling problem of AGV and dual Automated Yard Cranes (AYCs) and considers the influence of the interference in the dual AYCs on the overall scheduling efficiency. is article considers both the amount of operating equipment and the operating time. Aiming at minimizing the total energy consumption of the operating equipment, a mixed-integer programming model is established. Due to the complexity of the problem model, this paper designs a multilayer genetic algorithm to solve the problem. In addition, large-and small-scale numerical examples are analysed, and the effectiveness of the GA is verified by comparing CPLEX with the GA. e main conclusions of this article are as follows. (1) By evaluating the performance of the algorithm for solving AGV and AYC collaborative scheduling problems under different scale tasks, it is concluded that the multilayer genetic algorithm is very effective at solving complex scheduling problems. When the container's loading and unloading tasks are determined, the collaborative scheduling optimization scheme of AGVs and AYCs can be obtained by using the multilayer genetic algorithm. (2) By changing the location of the temporary storage area, CPLEX and the GA are used to determine the optimal location of the temporary storage area, and the walking paths of the dual AYCs are also obtained. is fully proves the effectiveness of the proposed model and algorithm to solve the interference problem between double AYCs.
is article only considers the coordinated scheduling between AGVs and AYCs. However, the scheduling of automated container terminals requires the integration of automated shore cranes, automated guided vehicles, automated yard cranes, trucks, trains, etc. In future research, we will focus on the cooperative scheduling of multiple operating equipment while considering uncertain environments,  Data Availability e data in the Qingdao port automation terminal is from Qingdao port official website (http://www.qdport.com/).

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