An Improved Shuffled Frog Leaping Algorithm for Multiload AGV Dispatching in Automated Container Terminals

Multiload AGVs, which can carry more than one container at a time, are widely used in automated container terminals. ,e dispatching decisions for multiload AGVs serving in automated container terminals on the target of minimum travel distance are significant in the process of container transportation in terms of improving operating efficiency. Previous work usually focused on AGVs working in a single-carrier mode, which was not only inconsistent with actual circumstances but also a waste of resources. In this paper, we establish a new mathematical model to describe multiload AGVs operating in automated container terminals, which is closer to the actual situation in real terminals. Based on this improved model, we propose a priority rule-based algorithm, termed as shuffled frog leaping algorithm with a mutant process (SFLAMUT), which can increase the diversity of the population and improve the convergence rate. Experiments were carried out based on data generated randomly according to the working properties of container terminals, and it is observed that the proposed SFLAMUT presents an effective and efficient exploration process and yields promising results in solving the proposed mathematical model.


Introduction
Automatic guided vehicles (AGVs), driverless and batterypowered vehicles, are applied extensively for container handling operations in container terminals. e use of multiload AGVs, which can carry up to two 20 ft containers or one 40 ft container in container terminals, is not fully investigated. In this paper, we try to investigate multiload AGV dispatching at highly automated container terminals, which can be divided into three main areas, namely, the berthing area, transport area, and yard area. Figure 1 illustrates the layout of automated container terminals. In the berthing area, containers are loaded or unloaded from a vessel by several quay cranes. e yard area is divided into several blocks, where containers are loaded and unloaded by stacking cranes. e sequence of the loading and unloading process is predetermined before a vessel arrives. e function of AGVs is to transport containers between the yard and the berth within the transport area according to the predetermined sequence.
In most container terminals, AGVs are used to carry a 20 ft or 40 ft container in a single-load carrier mode, which means an AGV cannot carry another container until it has completed the previous delivery task. In this way, there is no need to consider the capacity of AGVs and the size of containers, and thus the complexity of AGV dispatching is reduced. However, when delivering a 20 ft container, the half capacity of AGVs is not fully utilized, leading to a waste of resources and imposing a great burden on the traffic, especially when there are a large number of AGVs serving in the terminals. Given the high automation level of quay cranes and stacking cranes these days, the unloading and loading process of quay cranes and stacking cranes is of high efficiency. AGV dispatching in the transport area is becoming the bottleneck in automated container terminals.
Specifically, the AGV dispatching problem is composed of three parts, respectively, assigning tasks, routing AGVs, and traffic control [1]. Since algorithms for routing and traffic control are already integrated into the electronic control systems provided by AGV manufacturers, we only focus on algorithms for assigning tasks in automated container terminals. In the case of single-load AGVs, which can only carry one container at a time, the task assignment problem can be seen as the classical m : n assignment problem without considering the sequence of the pickup and delivery operations. While referring to multiload AGVs, the task assignment problem is more complicated, as we may not only take the capacity of AGVs and the size of containers (20 ft or 40 ft) into consideration but also develop a pickup-delivery rule to determine the sequence of carrying out pickup tasks and delivery tasks, which dramatically increases the space complexity and time complexity of this problem.
In order to improve the operation efficiency further, we establish a new mathematical model to describe the multiload AGV dispatching problem in automated container terminals. Based on the new model, we propose an improved shuffled frog leaping algorithm, which can increase the diversity of the population and exploit the potential to find the globally optimal solution to the dispatching problem.
e rest of this article is presented as follows: in Section 2, we provide a comprehensive view of existing works. In Section 3, we describe the multiload AGV dispatching problem we try to address in detail and establish a mathematical model based on the dispatching process we developed. e improved shuffled frog leaping algorithm used to solve the model is presented in Section 4. Finally, the computational experimental results of different algorithms are compared in Section 5, and conclusions are made in Section 6.

Research on AGV Dispatching
Problems. AGV dispatching has been studied extensively in the past decades. However, most of the studies focused on single-load AGV dispatching, assuming that there is only one type of container in the terminal, which is not the case in real practice. Bish et al. [2] divided the AGV problem into three subproblems, including assignment, scheduling, and routing. Based on the subproblems Bish et al. raised, the following researchers had made a further study. Zaghdoud et al. [3] proposed a hybrid approach that combined the Dijkstra algorithm, genetic algorithm, and a heuristic method to solve the assignment of containers to AGVs problem. Angeloudis and Bell [4] proposed a flexible dispatching algorithm along with a new AGV dispatching approach for job assignments problems in container terminal settings  under various conditions of uncertainty. Rashidi and Tsang [5] addressed a scheduling problem for automated guided vehicles in container terminals by designing an extended network simplex algorithm (NSA+); however, it was proved to be time-consuming when solving problems with a large number of tasks. Miyamoto and Inoue [6] formulated the dispatch and conflict-free routing problem of capacitated AGV systems (DCFRPC) as an integer program and proposed local/random search methods. A two-stage ant colony algorithm (ACA) was proposed by Saidi-Mehrabad et al. [7] to address the combinatorial problem of job shop scheduling and conflict-free routing.

Research on Multiload AGV Dispatching Problems.
Recently, as the latest trend in transfer equipment, the multiload AGVs show enormous potential in improving the operation efficiency of automated container terminals. Grunow et al. [1] developed a flexible priority rule-based approach for both single and dual-load vehicles and analyzed the MILP model in different scenarios concerning the total lateness of the AGVs. By combining different pickupdispatching rules, Azimi et al. [8] generated several control strategies and determined the best control strategy in terms of different criteria such as system throughput (ST), mean flow time of parts (MFTP), mean tardiness of parts (MTP), and AGV idle time (AGVIT). Fazlollahtabar et al. [9] proposed a mathematical program and an optimization method in two stages, namely, searching the solution space and finding optimal solutions, by minimizing the penalized earliness and tardiness. Chawla et al. [10] presented the modified memetic particle swarm optimization algorithm (MMPSO), which is a combination of particle swarm optimization (PSO) for global search and memetic algorithm (MA) for local search, for scheduling of multiload AGVs in the flexible manufacturing systems. Ho and Liu [11] addressed the problem of load selection and pickup-dispatching of multiple-load AGVs by proposing various loadselection rules and pickup-dispatching rules, aiming at understanding not only the performance of load-selection rules and pickup-dispatching rules in two performance criteria (throughput and tardiness) but also the effects that pickup-dispatching rules and load-selection rules have on each other's performance. Ho and Liu [12] also focused on the pickup-dispatching problem, for which nine pickup rules are proposed and studied in different performance measures, e.g., the system's throughput, the mean flow time of parts, and the mean tardiness of parts. Huo and Hu [13] established a mixed-integer linear programming model for multiload AGV dispatching and solved it by applying the genetic algorithm. Huo et al. [14] proposed GUROBI and the genetic algorithm to solve the single-load and multiload AGV dispatching problem and tested and verified the superiority of multiload AGV.

Intelligent Optimization Algorithms.
In order to address the multiload AGV dispatching problem efficiently and effectively, developing an appropriate algorithm is necessary. Previous work usually concentrated on the genetic algorithm, which can fall into a local optimum easily, especially when solving complicated problems. Some research shows that intelligent optimization algorithms can address many complicated computing problems efficiently. Among them, shuffled frog leaping algorithm (SFLA) is one of the most commonly used algorithms, which can improve the searching ability theoretically. However, the iterative times will increase as the problem becomes more complicated when SFLA converges to the global optimum solution [15]. Later, Niknam et al. [16] integrated a mutant process into SFLA to increase the diversity of the frog population, which can improve the searching ability of the global optimum solution. However, improving the convergence rate of the problem with large scale and high complexity is still needed. Quantum computation has many advantages, such as diversity and parallelism. erefore, the intelligent optimization algorithm and quantum computation are combined to generate the quantum-inspired algorithm which attracts extensive attention in the field. e quantum-inspired algorithm is mainly based on the advantages of quantum diversity, which has a relatively large searching space in a small-scale range.
e authors of References [17,18] proposed a quantum-inspired shuffled frog leaping algorithm (QSFLA) to improve the diversity further. However, it increases the randomness of the solution which leads to the slowness of convergence speed. To solve the randomness of the solution, Sun et al. [19] once proposed quantum-behaved particle swarm optimization (QPSO) based on delta potential well, the solution space of which is transformed from quantum space into real space. Subsequently, Li et al. [20] further proposed QPSO for continuous space optimization. Quantum bit's phase was used to describe the solution updated by the quantum rotation gate, which is used as a mutant process to avoid falling into the local optimum solution and improve the convergence speed.
From the study of previous works, it is learned that multiload AGV dispatching in automated container terminals is not fully investigated. Considering the research gap, we establish an improved mathematical model to describe the multiload AGV dispatching problem and propose a new algorithm to find a better solution to this model.

Problem Statement
In the case we consider, the multiload AGVs serving in container terminals can carry one 40 ft container or two 20 ft containers at one time, so we mainly focus on the dispatching problem of multiload AGVs which can carry up to two containers at a time and describe this kind of problem in details in this section.
A new assignment of tasks to available AGVs has to be determined, every moment a dispatching request is initiated. e definition of the multiload AGVs' availability is sparklingly important in this case. e concept of the availability of multiload AGVs has been introduced by Grunow et al. [1]. An AGV is considered unavailable if its capacity is fully utilized by the tasks assigned to it, otherwise available if an AGV can provide capacity partially or fully.

Mathematical Problems in Engineering
Based on the availability status of vehicles Grunow et al. proposed, we investigate the dispatching process and availability status of multiload AGVs further. An AGV is fully available after it unloads a 40 ft container or the second 20 ft container or two 20 ft containers simultaneously, and an AGV is partially available after loading the first 20 ft container. On the contrary, an AGV is unavailable after loading the second 20 ft container or a 40 ft container or two 20 ft containers simultaneously. As presented in Figure 2, different kinds of lines represent different kinds of status. e solid line means that the AGV is not available both during the trip to the next stop and at the next stop, while the dotted line and dashed line represent that the AGV is fully available or partially available, respectively. For an AGV with a different status, tasks at the next stop vary. A fully available AGV can carry a 40 ft container or a 20 ft container or two 20 ft containers simultaneously, and a partially available AGV is allowed to load a second 20 ft container or unload the first 20 ft container it carries. However, an unavailable AGV can only carry out the delivery task to unload the container after finishing the previous pickup task. In order to determine which type of container should be carried and whether the next task is a pickup task or a delivery task, we develop a priority rule presented in the subsequent sections.
Based on the multiload AGV dispatching process stated above, we formulate the mathematical model of the dispatching problem as follows, together with corresponding assumptions. e formulated dispatching problem is subjected to several constraints and optimized for minimum travel distance.

Basic Assumptions.
e following assumptions are considered in the present study: (1) At the initial stage, all of the AGVs are empty and are in the waiting place (2) e distance between two points is known in advance, and AGVs travel from one point to another point through the predetermined guide path (3) e location and number of delivery points and pickup points are fixed during the whole dispatching process (4) e number of AGVs is also fixed (5) Quay cranes, yard cranes, and multiload AGVs are reliable (6) ere is no task preemption (7) ere is no traffic congestion (8) After completing all the tasks, the AGVs will travel to the waiting place

Model Formulation.
Before introducing the newly proposed mathematical model, the following notation has to be defined. Parameters include the number of tasks, the number of AGVs, the node-set of all task points for AGV k , all available routes for AGV k , the maximum capacity of AGV k , and the traveling distance between different task points, which are known in advance according to the operation plan of automated container terminals. e detailed descriptions are shown below:

Parameters
n: number of tasks u: number of task points (including pickup points, delivery points, starting point, and termination point) K: number of AGVs P: node-set of pickup points D: node-set of delivery points T � P ∪ D: set of nodes for all task points P k : node-set of pickup points for AGV k D k : node-set of delivery points for AGV k T k � P k ∪ D k : node-set of pickup and delivery points for AGV k s k : starting point of AGV k d k : termination point of AGV k V k � T k ∪ s k , d k : node-set of all task points for AGV k A k � V k × V k : all available routes for AGV k C k : the maximum capacity of AGV k d i,j,k : the traveling distance between points i and j Nodes i and n + i: the pickup and delivery point of the ith task, respectively n i : the number of containers transported from point i to point n + i Decision variables include the routes AGV k visits, the number of containers AGV k carries, and the type of container loaded by AGV k , which are determined according to the shortest distance rule proposed in Section 4. e detailed descriptions are presented below:  l i,k : the number of containers AGV k carries, after completing the task at point i. Initially, l i,k � 0. s i,k : the type of container loaded at point i by AGV k . If it is 40 ft, s i,k � 1. If it is 20 ft, s i,k � 0.

Decision Variables
Applying the above notation, multiload AGV dispatching problem in container terminals can be formulated as an integer programming model. e objective function and constraints of the problem are presented as follows:

Constraints
j∈T k j∈P k Equation (1) is the objective function that minimizes total travel distances. Equations (2) and (3) are the constraints on the beginning and the end of the stage, which indicates that multiload AGVs must start traveling from the starting point before performing any tasks and move to the termination point to wait for the new pickup request after completing all tasks. Equation (4) illustrates that every task can only be executed once by one AGV. Equation (5) ensures that the pickup and delivery task of a container must be carried out by one AGV. Equation (6) means that if an AGV arrives at a task point, it will exit that point. Equations (7) and (8) are constraints on the sequence of delivery tasks and pickup tasks, which indicates that an AGV must go to the associated delivery point to unload the container if it has just loaded a 40 ft container, and an AGV can visit the associated delivery point or any pickup point where 20 ft containers are waiting to be loaded after loading a 20 ft container. Equations (9)∼(13) are constraints on the capacity of an AGV, which illustrate that the number of containers carried by the AGV will never exceed its capacity and the number of containers transported by the AGV from one point to another point will never exceed its load.

Shuffled Frog Leaping Algorithm with a Mutant Process (SFLAMUT)
Multiload AGVs are assigned several tasks in the very beginning, and then each AGV carries out the tasks assigned according to a predetermined AGV control process, which works as a rule giving multiload AGV guidance about task determination (whether a delivery task or a pickup task is performed). To solve the model established above, we propose a priority rule-based algorithm, called shuffled frog leaping algorithm with a mutant process (SFLAMUT), which conforms to the shortest distance rule and an improved shuffled frog leaping scheme. e shortest distance rule, under which the multiload AGV will give priority to the containers in the nearest pickup point or delivery point, has been proved to have the best throughput performance in AGV dispatching problem [12]. erefore, we adopt the shortest distance rule as our dispatching rule and design the multiload AGV control process based on this rule. When the operation begins, every AGV performs the pickup task or delivery task by following the control process illustrated in Figure 3.
SFLAMUT, based on the shuffled frog leaping scheme, is designed to address the integer programming problem as Section 3.2 stated, which has fast convergence speed and powerful searching ability to find the global optimum solution.

Population Initialization.
Assuming that the number of tasks is l, the number of AGVs is k, the frog population size is F, and the number of subgroups is m. Set the state of the ith frog as N i � [n i1 , n i2 , . . . , n il ], and the fitness value of its position as fitness i � f(N i ). e state of every frog in the frog population is initialized randomly. All frogs are sorted Mathematical Problems in Engineering in the ascending order according to the fitness value, before being divided into m subgroups containing s frogs, where m × s � F. e grouping process is stated as follows: the first frog is assigned to the first subgroup, the second frog is assigned to the second subgroup, . . . , the mth frog is assigned to the mthsubgroup, the (m + 1)th frog is assigned to the first subpopulation, . . . , and so on, until every frog individual in the frog population is assigned to a subgroup. e whole frog population can be expressed as follows: P � n 11 n 12 · · · n 1j · · · n 1l n 21 n 22 · · · n 2j · · · n 2l ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ n i1 n i2 · · · n ij · · · n il where n ij ∈ [1, k] is the serial number of an AGV, representing a specific AGV.

Fitness Function Calculation Based on the Shortest Distance Rule.
e fitness function is mainly used to evaluate the frog position, which reflects the extent that the frog is near the best solution. e fitness value of a frog is derived from the objective function shown in equation (1), which is calculated by following the multiload AGV control process based on the short distance rule (illustrated in Figure 3). e relationship between the fitness function and objective function is shown below. F GA represents the best solution the genetic algorithm generates: e pseudocode of the shortest distance rule and the calculation process of the objective function is presented in Pseudocode 1. Following parameters of this problem is predetermined: node-set of pickup points P � [p 1 , p 2 , . . . , p j , . . . , p l ]; node-set of delivery points D � [d 1 , d 2 , . . . , d j , . . . , d l ]; container size S � [s 1 , s 2 , . . . , s j , . . . , s l ]; node-set of tasks T � [t 1 , t 2 , . . . , t j , . . . , t l ], where t j � [p j , d j , s j ], which means the jth task is to transport the container of size s j from pickup point p j to delivery point d j ; the traveling distance between different task points distance � , where u represents the number of different task points and d ij represents the distance between the ith task point and the jth task point.

Updating with a Mutant
where randint is a random integer ranging from 0 to Δ(n), or Δ(n) to 0, where Δ(n) � ch lb u (t) − ch lw u (t), and t is the iterative time.
In each subgroup, ch lw u is updated by equation (17). During the updating process, a temporary variable new ch lw u � [n new u1 , n new u2 , . . . , n new uj , . . . , n new ul ] is used to store the updated solution. If new ch lw u is better than ch lw u , which means that the fitness value of new ch lw u is lower than the fitness value of ch lw u , new ch lw u will replace ch lw u . Otherwise, G will take the place of ch lb u in equation (16), and new ch lw u is recalculated. After completing the above process, if ch lw u is still not updated, a randomly generated new ch lw u will replace ch lw u . e previous process is repeated before reaching the local maximum iterative times. After the local search, all frogs are mixed and regrouped into m subgroups, and the above process of local search is repeated until reaching the boundary limits. e above update strategy may lead to convergence to the local optimal solution. To avoid falling into the local optimum and improve the diversity of the frog population, introducing a mutant process to this algorithm is necessary. e mutant process is developed according to the fitness value after all frogs being regrouped into subgroups, which can guarantee the convergence speed and avoid falling into the local optimum. e process is expressed as follows: rand ch l u � r 1 , r 2 , . . . , r j , . . . , r l , Initialize:/ * initialize tasks and AGVs T � InitializeTasks();/ * initialize the pickup point, delivery point, and size of container of tasks AGVs � InitializeAGVs();/ * initialize the status of AGV and assign tasks to AGVs randomly GetRoutes:/ * AGV route planning according to tasks assigned While there are tasks waiting or AGV is not in the waiting place:/ * for every AGV in AGVs IF AGV is in the waiting place or has just completed a pickup task or a delivery task:

Mathematical Problems in Engineering
If the fitness values of rand ch l u is better than ch lb u , ch lb u will be substitute by rand ch l u ; Otherwise, if mut ch l u is better than ch lb u , ch lb u will be substitute by mut ch l u ; otherwise, ch lb u will remain the same, and ch lw u will be substitute by rand ch l u .

Algorithm Process.
In summary, the algorithm process of SFLAMUT is stated as follows: (1) Initialize population P, and set iteration times t � 0 (2) Calculate the fitness value of every frog by equation (15), keep the best frog (3) Sort the frogs in population according to their fitness value in the ascending order and group the frogs into several subgroups (4) Update the worst frog in every subgroup by equations (16) and (17)  (5) Divide the updated population into new subgroups (6) Apply the mutant process by equations (18)∼(21) (7) Judge whether the termination criteria is met; if so, the current best frog is the optimum result yielded; otherwise, set t � t + 1 and turn to Step (2)

Local Convergence Analysis.
To facilitate the convergence analysis, we assume that the number of subgroup m � 1, indicating that the global best frog G equals the best frog ch lb u in the subgroup at the beginning. Describe equations (16) and (17) as follows. e tth iteration: e (t + 1)th iteration: where Round(x) is the integer number of x through rounding. Let e i be equal to the difference between the integer number of r(ch lb u − ch lw u (i − 1)) through rounding and r(ch lb u − ch lw u (i − 1)), equations (24) and (25) can be derived as follows: Substitute equations (24) and (25) into equation (26): where ]. If ch lb u − ch lw u (t) ≥ 0, we can yield that Δch lw u (t + 1) ≥ 0 and Δch lw u (t) ≥ 0. In this case, − Δch lw u (t) ≤ ΔRound ≤ 0. If Δch lw u (t) ≠ 0, equation (28) can be derived: Let q(t) be equal to 1 + (ΔRound/Δch lw u (t)). From equation (28), we can identify that 0 < q(t) < 1. In this case, lim t⟶∞ Δch lw u (t) � 0 and lim t⟶∞ ch lw u (t) � ch lb u , indicating that the result converges to the local optimal solution.
If ch lb u − ch lw u (t) ≤ 0, we can yield the same result according to the above derivation similarly.

Global Convergence Analysis.
As SFLAMUT is a kind of random search techniques, we determine whether SFLAMUT can get the global optimization solution through the convergence criteria of the stochastic algorithm proposed by Solis and Wets [21] as follows.
For optimization problem 〈A, f〉, we seek a point x in A which minimizes f on A. D is a random search algorithm, X k is the solution of the kth iteration, and the solution of the next iteration is X k+1 � D(X k , α), where A is the feasible solution space, f is the fitness function, α is a solution algorithm D once found.
Condition H1: f (D(x, α) Convergence eory (Global Search) [21]: suppose that f is a measurable function, A is a measurable subset of R n 8 Mathematical Problems in Engineering and Conditions H1 and H2 are satisfied. Let x k ∞ k�0 be a sequence generated by algorithm D. en, where R ε is the set of global best solutions and P(x k ∈ R ε ) is the probability that the point x k generated by the algorithm is in R ε at step k.
To verify whether SFLAMUT can get the global optimization solution or not, we just need to verify whether SFLAMUT satisfies Conditions H1 and H2. For SFLAMUT, algorithm D can be defined as follows: It is obvious that algorithm D satisfies Condition H1. If algorithm satisfies Condition H2, the union of all solutions algorithm D once found can consist of the feasible solution space, that is, where R i,k is the solution algorithm D found at iteration times k and A is the feasible solution space. From equation (22), we have that From local convergence analysis, we have proved that lim t⟶∞ ch lw u (t) � ch lb u , so when t ⟶ ∞, , which means that the algorithm cannot ensure the global convergence if just updated through equations (16) and (17). erefore, we introduce the mutant process (equations (18)∼(21)) to the updating process. X i (t) will be generated randomly through the mutant process after every iteration. In this case, if t ⟶ ∞, we have that ∪ F i�1 R i,t � A, indicating that equation (31) is satisfied. In this case, for any subset B of A, we have that 0 ) k � 0, which means the probability that algorithm D cannot find the solution in B is zero, and therefore, Condition H2 is satisfied.
In conclusion, the updating process of equations (16) and (17) can ensure that the algorithm converges to the current optimal position; however, it cannot guarantee the global convergence. When equations (16) and (17) cannot generate a better solution, equations (18)∼(21) randomly generate a solution to substitute a current solution, which can expand the sample space. In this case, as iteration times k ⟶ ∞, the probability that SFLAMUT converges to the global optimum solution is 1, which indicates the global convergence of SFLAMUT.

Description of the Experiments.
e proposed mathematical model and SFLAMUT were tested by computational experiments on hypothetical automated container terminals, in which multiload AGVs are employed. e computational experiment was performed by the Python program. We choose algorithms including genetic algorithm (GA), shuffled frog leaping algorithm (SFLA), and quantum-inspired shuffled frog leaping algorithm (QSFLA) as our benchmarks. e difference between the given algorithms is presented as follows: (1) GA: traditional genetic algorithm with roulette wheel selection process, single-point crossover process, and a mutant process (2) SFLA: traditional shuffled frog leaping algorithm with improved updating process suitable for solving integer programming problems (3) QSFLA: shuffled frog leaping algorithm and quantum computation are combined to generate QSFLA, the solution space of which is transformed from quantum space into real space (4) SFLAMUT: an improved shuffled frog leaping algorithm with a mutant process as stated in Section 4

Parameter Settings.
e parameters in the mathematical model we proposed were set as follows: e number of pickup points was set to 5, the number of delivery points to 5, and the number of multiload AGVs to 5. e parameters of different algorithms are presented in Table 1. e sign "-" represents that the parameter of the algorithm does not exist.
We used GA, QSFLA, SFLA, and SFLAMUT to solve the mathematical model, respectively, and evaluated the performance of the algorithms by comparing the convergence speed, optimal results, and stability at a fixed hybrid iteration number of 500.
For tasks of different numbers ranging from 10 to 80, each algorithm was executed 50 times independently in order to reduce the impact of random factors. We calculated the best value, worst value, and mean value of the optimal results generated by 50 experimental runs of each algorithm. Details are shown in Table 2. Table 2, we can figure out that when it comes to searching for a better optimum solution, intelligent optimization algorithms (SFLA, QSFLA, and SFLAMUT) perform better than the traditional heuristic algorithm (GA) in terms of the best value, worst value, and mean value in most cases. is is because GA with selection and crossover process will retain the genes of the best individuals to the next generation, leading to falling into the local optimum especially when the population size is small. However, SFLA, QSFLA, and SFLAMUT with the updating process of global information exchange can avoid falling into the local optimum to some extent. Although SFLA, QSFLA, and SFLAMUT perform better with respect to optimal solutions, they have longer running time than GA due to a more complex searching process. is limitation can be accepted in the case we consider because the task assignment process should be completed in advance, not in operation, which means there is plenty of time to execute each algorithm. Among intelligent optimization algorithms (SFLA, QSFLA, and SFLAMUT), SFLAMUT we proposed yielded the minimum best value, the minimum worst value, and the minimum mean value, indicating far better optimal results than the other three algorithms in all experiments we carried out. at means SFLAMUT has the potential to find the global optimum solution, despite a slightly longer running time. Besides, SFLAMUT has a relatively low value with respect to the standard deviation, which indicates better search accuracy and credibility. e optimization process of these four algorithms is compared in Figures 4-11. As depicted in each figure, the vertical axis is the objective function, whose value is obtained from the mean value of 50 experimental runs of each algorithm, and the horizontal axis is the evolutionary generation (iteration times).

Results and Discussion. From
As Figures 4-11 show, intelligent optimization algorithms (SFLA, QSFLA, and SFLAMUT) have a faster convergence speed than the traditional heuristic algorithm (GA) in most cases. However, as the complexity of the problem increases, QSFLA and GA have almost similar performance concerning convergence speed and optimum results, which indicates QSFLA performs worse when solving complex problems. is is reasonable because QSFLA transforms the solution space from real space into quantum space, which improves the diversity but increases the randomness of the solution, leading to the slowness of convergence speed, especially when solving complicated problems. erefore, when the iteration times are fixed, the results QSFLA yielded become worse and worse as problems become more and more complex (as shown in Figures 9-11). SFLA and SFLAMUT always perform better than QSFLA and GA because they have similar improved updating process as equations (16) and (17) show. Such an updating process maintains the advantages (avoiding falling into the local optimum and a faster convergence speed) of traditional SFLA and is suitable for solving the integer programming problem in the meantime. Compared with SFLA, SFLAMUT has the best performance among the four algorithms. We add the mutant process (equations (18)∼(21)) into the updating process of SFLAMUT, which enables SFLAMUT to randomly generate a solution compared with the local best solution, leading to more chances to find the global optimum solution when iteration times are fixed.

Conclusions
In this paper, we establish an improved mathematical model for multiload AGV dispatching in automated container terminals, aiming at fully utilizing the capacity of multiload AGVs and improving its operation efficiency further. To find the global optimal solution to this model, we propose the shortest distance rule-based shuffled frog leaping algorithm with a mutant process (SFLAMUT). e experimental results of SFLAMUT outperform the other three algorithms (GA, SFLA, and QSFLA). Future research work can be extended to analyze the proposed model by considering other objective functions with respect to maximizing capacity utilization or minimizing total travel time. More in-depth research of other efficient algorithms to address the dispatching problem, such as machine learning algorithms, can be conducted based on the mathematical model proposed in this paper.
Data Availability e data used in the experiments are generated randomly based on the actual flow of automated container terminals. All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest
e authors declare no conflicts of interest.