Vehicle Routing Optimization with Cross-Docking Based on an Artiﬁcial Immune System in Logistics Management

: Background: Manufacturing companies optimize logistics network routing to reduce transportation costs and operational costs in order to make proﬁts in an extremely competitive environment. Therefore, the efﬁciency of logistics management in the supply chain and the quick response to customers’ demands are treated as an additional source of proﬁt. One of the warehouse operations for intelligent logistics network design, called cross-docking (CD) operations, is used to reduce inventory levels and improve responsiveness to meet customers’ requirements. Accordingly, the optimization of a vehicle dispatch schedule is imperative in order to produce a routing plan with the minimum transport cost while meeting demand allocation. Methods: This paper developed a two-phase algorithm, called sAIS, to solve the vehicle routing problem (VRP) with the CD facilities and systems in the logistics operations. The sAIS algorithm is based on a clustering-ﬁrst and routing-later approach. The sweep method is used to cluster trucks as the initial solution for the second phase: optimizing routing by the Artiﬁcial Immune System. Results: In order to examine the performance of the proposed sAIS approach, we compared the proposed model with the Genetic Algorithm (GA) on the VRP with pickup and delivery benchmark problems, showing average improvements of 7.26%. Conclusions: In this study, we proposed a novel sAIS algorithm for solving VRP with CD problems by simulating human body immune reactions. The experimental results showed that the proposed sAIS algorithm is robustly competitive with the GA on the criterion of average solution quality as measured by the two-sample t-test.


Introduction
The concept of a smart city includes a high degree of information technology (IT) integration and communication, which is the same concept as supply chain management today, which relies on the use of various ITs and techniques.One of the important elements of smart cities is smart logistics, aiming to help manufacturers gain performance from reusable transport packaging, such as pallets, racks, and bins, as well as tracking packages.By building Internet of Things (IoT) monitoring technology, smart logistics IT systems can be integrated into virtually any transportation asset to track its location in real-time, optimize inventory planning, and monitor environmental conditions, especially when the world enters the fifth-generation mobile communication technology or even the six-generation mobile communication technology in the future [1-3].However, even with advanced IT today, logistics costs are still an important part of a global company's worldwide operation management to construct an enterprise's competitive edge over competitors, especially when dealing with post-COVID-19 supply chain disruptions.Intelligent logistics network functions allow companies to quickly respond to customer requirements and create their own competitive advantage over competitors in the design of smart logistics by using IT, such as radio frequency identification, to track shipments in real-time.
gistics by using IT, such as radio frequency identification, to track shipments in real-time.Moreover, logistics costs weigh significantly on a companyʹs total costs.The authors of [4] identified the logistics costs from three different aspects of a case company: (1) the share of procurement costs reduced from approximately 0.65 (2013) to 0.45% (2015); (2) the share of production costs increased from approximately 50 (2013) to 65% (2015); and (3) the share of sales costs increased from approximately 7 (2013) to 8.5% (2015).As a result, in order to increase overall profits, lower operational costs, and improve a company's service level, a well-designed and highly cost-efficient logistics mechanism becomes essential.The paper focuses on controlling transport costs by optimizing routes and scheduling of vehicles, which is known as the vehicle routing problem (VRP) [5].
VRPs are a fundamental activity in the fields of transportation systems, distribution channel design, and logistics networks.The study of the VRPs is indeed vital in optimizing the physical flow of goods in the logistics operation in order to reduce transport costs and increase supply chain network performance.Among numerous studies of the VRPs, the problem with optimal routes and scheduling of vehicles, considering both pickup and delivery processes simultaneously, is called the VRP with cross-docking (VRPCD).Kulwiec (2004) pointed out that the cross-docking (CD) facility is an important supply chain strategy.They classified the CD facility into six different types.One of the CD operations is "Truck/Rail Consolidation." [6].Both suppliers' and retailers' sides of the supply chain have to be considered at the same time in the process of transporting goods.Therefore, CD facilities are an important element in a synchronized supply chain as well as in sustainable logistics management.Generally, there would be no interruption between upstream and downstream operations as long as the physical flow of goods in the pickup process is simultaneously delivered at the CD depot and then delivered to customers after the consolidation process.Therefore, no inventory was stocked, and no delay in customer orders occurred at the facility.As a result, the construction of a CD system in the logistics network makes companies able to: (1) facilitate the efficiency and effectiveness of supply chain management; (2) lower space requirements and reduce transportation costs; and (3) better control the distribution process [7]. Figure 1 shows a simple layout of the CD depot.Studying an efficient heuristics methodology is essential to acquiring an optimal or near-optimal solution within a reasonable amount of computation time because the VRP is a well-known non-deterministic polynomial-time-hard (NP-hard) problem.As the global pandemic of the new coronavirus (COVID-19) has shaken the world since 2020, the Studying an efficient heuristics methodology is essential to acquiring an optimal or near-optimal solution within a reasonable amount of computation time because the VRP is a well-known non-deterministic polynomial-time-hard (NP-hard) problem.As the global pandemic of the new coronavirus (COVID-19) has shaken the world since 2020, the study of the human immune system has piqued the interest of many researchers.The Artificial Immune Systems (AIS) is the simulation algorithm of human body defense systems and is applied to solve many research fields, such as clustering/classification, fault detection, combinatorial optimization, and learning problems [8][9][10][11][12].
The AIS has presented numerous studies on solving optimization problems.The results showed that their algorithm was a feasible and effective method for the VRP.For example, the authors of [13] applied clonal selection to tackle the VRP by using clonal selection operators, super mutation operators, and clonal proliferation to improve global convergence speed.The results indicate that their algorithm has a remarkable reliability of global convergence and avoids prematurity when solving the VRP effectively.
This paper aims to propose a two-phase optimization approach based on an artificial immune system combined with the sweep method, called sAIS, to solve the vehicle routing problem with the CD facilities in the logistics network as an important element of smart city infrastructure.The two-phase approach begins with the grouping method to fulfill vehicle capacity constraints and is followed by the optimization engine to find near-optimal solutions for each truck routing sequence while utilizing a minimum number of trucks.The experimental results showed that the proposed sAIS algorithm is robustly competitive with the GA on the criterion of average solution quality.
The remainder of this paper is organized as follows.Section 2 reviews the related literature, and Section 3 defines the VRPCD formulation.Section 4 details the methodology and procedure of our proposed sAIS heuristic algorithm, following a number of experimentation examples presented in Section 5. Section 6 concludes the research with a summary based on the computational results.

Related Works
The CD system is a lean supply chain model of transporting raw materials or products from pickup to delivery without ever storing them in the warehouse.It can significantly reduce inventory levels, required space, handling costs, and lead times, as well as customer response times.Packages are unloaded from inbound trucks immediately after arriving at the depot, followed by a sorting, repacking, and dispatching process as shown in Figure 2, then loaded onto outbound trucks for delivery to retailers in a distribution channel [14][15][16][17][18].The primary objective is to avoid a high inventory level and reduce handling costs so that there will be no inventory being stored in the depot.This is the same concept as the Toyota Production System or lean operations.A well-designed CD operation can provide companies with significant benefits, including decreased inventory levels, low storage space requirements, low transportation costs, a fast response to customer requirements, and better control of the distribution process.Both receiving and shipping processes must be considered simultaneously, to enable the CD facility to be integrated into the logistics network effectively.Lee, Jung, and Lee (2006) first proposed that the most important process for CD operations is the pickup phase, in which all trucks must arrive at the depot at the same time in order to control the start time of warehouse operations [17].There are four types of immune mutations running inside the human body: IgM, IgG, IgE, and IgA [47][48][49][50].In this paper, we simulate each mode having a different mathematical function.First, the calculation of affinity is the total correlations between the antigenantibodies and antibodies-antibodies in our mathematical programming model.If the affinity value of a new string in IgM mode is smaller than that of the old one, a total of four Before the CD system served as a mathematical constraint, Dantzig and Ramaser (1959) initially introduced the VRP concept as a solution to the "Truck Dispatching Problem."In their study, a linear programming model was designed to acquire a near-best solution for the truck scheduling problem, which was concerned with the optimum routing of a fleet of delivery trucks supplied by the terminal [19].Following that, numerous research proposals were made to address the developing VRPs in their study.The VRPs have been studied for many decades.The VRP is a set of customers with known demands who are serviced by a fleet of trucks from one or more depots to a number of geographically dispersed locations and customers based on optimally designed routes [20][21][22][23][24].
There are many solving methods for variants of the VRP proposed in the literature.The authors of [25] dealt with capacitated VRP and distance restrictions by using an integer programming method that used a constraint relaxation approach and sub-tour elimination.The author of [26] proposed tour-partitioning heuristics to solve the pickup and delivery VRPs, whereas the authors of [27] developed a hybrid heuristic method based on the Genetic Algorithm (GA) with neighborhood search to solve the basic VRPs.They showed that the hybrid GA had a significant improvement over the pure GA and was competitive with the simulated annealing approach [28] and the Tabu search approach [29][30][31] in their experiment results.The authors of [32,33] proposed hybrid ant colony optimizations (ACOs) to solve the VRPs with time windows and found that they were applicable and effective in practical problems.The authors of [34] used an ACO approach for the multiple VRPs with pickup and delivery along with time windows and heterogeneous fleets and applied it to large-scale problems.The authors of [18] developed a matheuristic approach consisting of two phases: adaptive large neighborhood search (ALNS) and setting partitioning to solve VRPCD.The authors of [35] proposed a particle swarm optimization approach to solve VRPCD and carbon emissions reduction.Moreover, sustainable logistics management is a popular research topic nowadays to follow the United Nations' sustainable development goals [36].
Following the same concept to simulate an ant's behavior as ACO, Jerne (1974) proposed the first AIS model to simulate the immune system as a mathematical formulation to solve optimization problems, which had an interaction network of lymphocytes and molecules that had variable regions [37].Following Jerne's research, Farmer, Packard, and Perelson (1986) proposed the dynamic immune network, which could simulate and solve classification problems, showing the AIS can be extended to solve big data prediction problems nowadays as machine learning pioneers [38].Kephart (1994) 2022) proposed an immune system programming with local search (ISPLS) algorithm with a tree data structure to be used in the meta-heuristics programming approach to develop new practical machine learning tools [42].As a result of many years of development, the AIS has become a well-known meta-heuristic that is widely applied to solve combinatorial optimization and abnormal detection problems.
The AIS mainly simulates the relationship between antigens and antibodies, and the core idea of immune reactions includes antibodies' reproduction, clonal expansion, and immune memory properties in the biological immune system.Organisms have two kinds of immunity: the innate immune system and the adaptive immune system.The innate immune system is capable of recognizing molecular patterns in pathogens and signaling other immune cells to start fighting against the pathogens.Adaptive immune systems can maintain a stable memory of known patterns.Living organisms use the immune system to defend their bodies from invasion by outside substances.Lymphocytes are parts of the immune system, which includes T and B cells.Both types of lymphocytes have surface receptors capable of recognizing molecular patterns present on antigens (binding with epitopes).When the receptors bind to epitopes and exceed a threshold, a lymphocyte becomes activated.
Activation triggers a series of reactions that can lead to the elimination of pathogens.During the infection response, the immune system's B cells produced antibodies.Pathogens are bound by antibodies called antigens.Antigens and antibodies can bind together when complementary shapes exist.After binding, the antibody disables the pathogens so that the immune system can easily destroy them.Figure 2 depicts how to select immune system antibodies [43][44][45][46].
There are four types of immune mutations running inside the human body: IgM, IgG, IgE, and IgA [47][48][49][50].In this paper, we simulate each mode having a different mathematical function.First, the calculation of affinity is the total correlations between the antigenantibodies and antibodies-antibodies in our mathematical programming model.If the affinity value of a new string in IgM mode is smaller than that of the old one, a total of four kinds of immune mutations, including IgG mode, IgE mode, and IgA mode, are randomly selected to operate in the next step of mutation in the same generation for optimization purposes.

1.
Somatic Hypermutation Simulation Function: • IgM mode: Inverse mutation is used in the IgM mode.In the sequence s, randomly selected two positions, i and j.Inverse the sequence of cells between the i and j positions in the neighbor of s.It should be noted if |i − j| < 2 it cannot be mutated; • IgG mode: Pairwise mutation is used in the IgG mode.In the sequence s, randomly selected two positions, i and j.Swap these two cells in the neighbor of s; • IgA mode: Insertion mutation is used in the IgA mode.In the sequence s, randomly selected two positions, i and j.Insert the cell i to the position j in the neighbor of s; • IgE mode: Both the swap and insertion mutations are used in the IgE mode.In the sequence s, randomly selected two positions, i and j.A new sequence s' is provided by swapping i and j.Then, in the sequence, s' randomly selected i' and j' to be a cell and a position, respectively.Inserting cell i at position j in the neighbor of s'.

Affinity Maturation Simulation Function
Affinity is a positive or negative correlation between an antigen and an antibody.When the affinity is higher, it can generate a better fit between the interacting surfaces of the antibody and antigen.After somatic hypermutation, some variant antibodies with higher affinities may be produced.Then, in the next response, the cells will have a greater affinity.This phenomenon is called affinity maturation;

Elimination Simulation Function
The human body's bone marrow generates billions of B cells into circulation every day.Those B cells are saved in a catalog, where there is limited storage space for the B cells.The B cells are dead within a certain period if they meet the specific antigen.Hence, the B-cell catalog is not static, and newly produced antibodies are continually tested against antigen infection.Therefore, new antibodies are generated and old ones are deleted.

VRPCD Formulation
The VRP is traditionally illustrated as a graph network with vertices denoting terminal points and arcs as vehicle routes, as shown in Figure 3.The basic notations for all VRPs are shown as follows: where V = {v 0 , . . ., v n } is a vertex set; v 0 represents the central warehouse from which deliveries are made; where V = {v0, …,vn} is a vertex set; v0 represents the central warehouse from which deliveries are made; A = {(vi, vj): i ≠ j, vi, vj ∈ V} represents the directed arc set; E = {(vi, vj): i < j, vi, vj ∈ V} is a set of undirected edges.Lee, Jung, and Lee (2006) addressed the constraints of vehicle routing with the one CD warehouse problem.The CD facility plays a key role in synchronizing the distribution process on both sides of the supply chain [17].As simulated, the CD facility can be treated as the home depot as in the traditional VRPs, except that the model of CD specifies the simultaneous arrival of each vehicle from the receiving trip.Therefore, several assumptions are made.First, we have n nodes, denoting a total number of suppliers and retailers, who are serviced by m vehicles.Every truck must depart and come back to the depot (i = 0), with the simultaneous arrival of trucks from pickup routes.Second, for each customer, only one truck is assigned and associated with a cost amount of Cij.Every customer has a homogeneous demand d, which is related to the capacity limit qk for truck k.Additionally, a constraint of the planning time horizon T specifies that the total distance traveled by vehicles cannot exceed it.Two types of costs are considered: transportation costs and operational costs.The objective of this research in the mathematical programming model is to find optimized routing solutions by using a minimum number of trucks to service and complete the assignment.The following presents the basic notations of the VRPCD model.

Decision Variables:
Xijk: a binary variable representing the route from i to j is serviced by vehicle k, where

Notation of Variables:
Yijk: loaded quantity of vehicle k from pickup trip i to j; Zijk: unloaded quantity of vehicle k from delivery trip i to j; tcijk: the transportation cost of vehicle k from customer i to j; etijk: time for vehicle k to move from i to j; δik: service time required by vehicle k to load/unload the quantity demanded at i; Lee, Jung, and Lee (2006) addressed the constraints of vehicle routing with the one CD warehouse problem.The CD facility plays a key role in synchronizing the distribution process on both sides of the supply chain [17].As simulated, the CD facility can be treated as the home depot as in the traditional VRPs, except that the model of CD specifies the simultaneous arrival of each vehicle from the receiving trip.Therefore, several assumptions are made.First, we have n nodes, denoting a total number of suppliers and retailers, who are serviced by m vehicles.Every truck must depart and come back to the depot (i = 0), with the simultaneous arrival of trucks from pickup routes.Second, for each customer, only one truck is assigned and associated with a cost amount of C ij .Every customer has a homogeneous demand d, which is related to the capacity limit q k for truck k.Additionally, a constraint of the planning time horizon T specifies that the total distance traveled by vehicles cannot exceed it.Two types of costs are considered: transportation costs and operational costs.The objective of this research in the mathematical programming model is to find optimized routing solutions by using a minimum number of trucks to service and complete the assignment.The following presents the basic notations of the VRPCD model.
Decision Variables: X ijk : a binary variable representing the route from i to j is serviced by vehicle k, where

Notation of Variables:
Y ijk : loaded quantity of vehicle k from pickup trip i to j; Z ijk : unloaded quantity of vehicle k from delivery trip i to j; tc ijk : the transportation cost of vehicle k from customer i to j; et ijk : time for vehicle k to move from i to j; δ ik : service time required by vehicle k to load/unload the quantity demanded at i; m: number of vehicles; n: number of customers; c k : fixed cost of vehicle k; q k : maximum capacity for each vehicle k; T: planning horizon; P: set of unit demand from each pickup stop; D: set of unit demand from each delivery stop; DT jk : departure time of vehicle k at node j; AT k : arrival time at the depot of vehicle k.Objective Function: (1) Subject to: Equation ( 1) is the objective function whose goal is to minimize both transportation costs and fixed costs.The constraints of each customer being serviced by only one vehicle are described in Equations ( 2) and (3), while Equation (4) represents that each vehicle arriving at the node must also leave from the same node.Equations ( 5) and ( 6) state the constraints of each vehicle being only permitted to start from and return to the crossdocking depot and being used to serve at most one route, respectively.The loaded and uploaded demand from pickup and delivery processes cannot exceed the vehicle quantity limit detailed by Equation (7), while Equations ( 8) and (9) each represent the quantity limit for pickup and delivery en route.The total distance visited and time traveled cannot exceed the planning horizon, which is limited by Equation (10).The departure and arrival times are functioned by Equations ( 11) and (12), respectively.The vehicles' simultaneous arrival at the CD hub is stated by Equation (13).

Clustering Method
In order to minimize total cost, each dispatched vehicle must serve the maximum number of customers within its capacity limit while traveling the shortest distance possible.In this research, we simulate this practical strategy to solve the problem of generating vehicle route sequences with one CD operation in the supply chain.Moreover, for faster convergence, we use a two-phase algorithm to solve VRPCD.
In the initial route-generating phase, we use the sweep method as the first phase to provide the initial solution for the second phase of sAIS.Usually, the sweep method is applied to a polar coordinate system, and the center of the coordinate system serves as the depot of the VRP.In this case, for each traveling distance of route O i (i = 1, 2, . . ., n, where n denotes the number of customers), the ith customer's C i served as the first customer in the cluster.Next, search for the closest customer from C i + 1 to C i by increasing the angle, and add it to the same cluster without exceeding vehicle capacity.When a customer cannot be added to the current cluster due to a vehicle capacity limit, this customer becomes the first customer in the next cluster.After all customers are clustered, we calculate the objective value of this route.After O 1 , O 2 , . . ., O n , were calculated, we chose the minimum distance of n routes as the initial solution of the AIS algorithm.Figure 4 illustrates an example of the clustering process for 12 nodes.
times are functioned by Equations ( 11) and (12), respectively.The vehicles' simultaneous arrival at the CD hub is stated by Equation (13).

Clustering Method
In order to minimize total cost, each dispatched vehicle must serve the maximum number of customers within its capacity limit while traveling the shortest distance possible.In this research, we simulate this practical strategy to solve the problem of generating vehicle route sequences with one CD operation in the supply chain.Moreover, for faster convergence, we use a two-phase algorithm to solve VRPCD.
In the initial route-generating phase, we use the sweep method as the first phase to provide the initial solution for the second phase of sAIS.Usually, the sweep method is applied to a polar coordinate system, and the center of the coordinate system serves as the depot of the VRP.In this case, for each traveling distance of route Oi (i = 1, 2, …, n, where n denotes the number of customers), the ith customer's Ci served as the first customer in the cluster.Next, search for the closest customer from Ci + 1 to Ci by increasing the angle, and add it to the same cluster without exceeding vehicle capacity.When a customer cannot be added to the current cluster due to a vehicle capacity limit, this customer becomes the first customer in the next cluster.After all customers are clustered, we calculate the objective value of this route.After O1, O2, …, On, were calculated, we chose the minimum distance of n routes as the initial solution of the AIS algorithm.Figure 4 illustrates an example of the clustering process for 12 nodes.The total cost is accumulated by the fixed cost of vehicles used (operational cost) and transportation costs as vehicles travel from one customer to another.In some cases, the aggregated service time consumed by each customer is also considered.Overall, the shorter the distance traveled by vehicles, the lower the total cost incurred.Figure 5 exemplifies the hierarchies of the proposed two-phase algorithm for 10 customers with homogeneous demand served by 2 vehicles as well as the chromosome encoding format in this study.
aggregated service time consumed by each customer is also considered.Overall, the shorter the distance traveled by vehicles, the lower the total cost incurred.Figure 5 exemplifies the hierarchies of the proposed two-phase algorithm for 10 customers with homogeneous demand served by 2 vehicles as well as the chromosome encoding format in this study.

AIS Procedure
As mentioned previously, the AIS was successfully implemented to solve combinatorial optimization and abnormal detection problems.For data representation, we used chromosomes to denote routing sequences.All suppliers and retailers were assigned a unique number.We made the variation field of each cell, and its dimensions (length of chromosome) equivalent to the suppliers' (for pickup routes) or customers' (for delivery routes) numbers to be served, as shown in Figure 5. Therefore, the cells' chromosomes represent the nodes of the positions to be visited.In the meantime, the suppliers' or customers' sequences also mean the routing of orders.First, somatic hypermutation will be selected to make cells evolve, which can search for any feasible cells in the solution field.Next, we calculate the affinities between the antigen-antibodies and antibodies-antibodies according to the objective function.After iterations, the relationship between antigens and antibodies has a greater affinity.Moreover, elimination is a mechanism that will delete worse cells with each iteration.With affinity maturation and elimination, the cells will mutate toward the better-quality field of solution.In this case, each iteration individually will preserve the best cell, which means a higher affinity value.
In the proposed sAIS optimization algorithm, the principles of each antibody mutation have four parts

AIS Procedure
As mentioned previously, the AIS was successfully implemented to solve combinatorial optimization and abnormal detection problems.For data representation, we used chromosomes to denote routing sequences.All suppliers and retailers were assigned a unique number.We made the variation field of each cell, and its dimensions (length of chromosome) equivalent to the suppliers' (for pickup routes) or customers' (for delivery routes) numbers to be served, as shown in Figure 5. Therefore, the cells' chromosomes represent the nodes of the positions to be visited.In the meantime, the suppliers' or customers' sequences also mean the routing of orders.First, somatic hypermutation will be selected to make cells evolve, which can search for any feasible cells in the solution field.Next, we calculate the affinities between the antigen-antibodies and antibodies-antibodies according to the objective function.After iterations, the relationship between antigens and antibodies has a greater affinity.Moreover, elimination is a mechanism that will delete worse cells with each iteration.With affinity maturation and elimination, the cells will mutate toward the better-quality field of solution.In this case, each iteration individually will preserve the best cell, which means a higher affinity value.
In the proposed sAIS optimization algorithm, the principles of each antibody mutation have four parts.The basic components of the sAIS applied in the VRPCD are shown as follows: • Antigen: the antigen mimics the equations of the VRPCD model, which represent the formulas for the objection function and constraint functions; Stopping criterion: the sAIS optimization procedure will not be terminated until the maximum number of iterations is reached.
The proposed sAIS optimization approach has a mechanism with the inside operation processes: 1.

Somatic Recombination
Somatic recombination is one of the gene rearrangement processes that involves cutting out small regions of DNA and then putting the remaining pieces of DNA back together in an error-prone way in the adaptive immune system.For every iteration of our proposed algorithm, we randomly create a new set of antibodies by using the somatic recombination process in order to search for optimized solutions;

Somatic Hypermutation
There are five types of immune mutations that we use in our proposed sAIS algorithm.Each type has a different function.If the cost of objection function of antibodies in IgM is larger than the cost of the old ones, we randomly selected one of four kinds of immune mutations to operate in the next step of mutation in the same iteration.In addition to the traditional four types of immune mutations, we proposed a new immune mutation, called IgG 2 , to further improve the performance of the sAIS optimizations:

•
IgM mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Inverse the sequence of suppliers or customers between i and j as shown in Figure 6.It is noted that if |i − j| < 2 it is not allowed to be mutated; • IgG mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Swap these two suppliers or customers in the neighbor of s; • IgA mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Insert the suppliers or customers i to the sequence's jth position in the neighbor of s as shown in Figure 7; • IgG 2 mode: In the vehicle's route of sequence s, let i 1 , i 2 , j 1 , and j 2 be randomly selected as four suppliers or customers in the sequence.Swap i 1 to j 1 and i 2 to j 2 in the neighbor of s as shown in Figure 8; • IgE mode: Both the IgG mode and the IgA mode are applied in this mutation.In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.A new sequence s' is provided by swapping i and j.Then, in the sequence s' randomly selected i' and j' will be suppliers and customers, respectively, and have a position in the sequence.Inserting the suppliers or customers i to the sequence position j in the neighbor of s'; FOR PEER REVIEW 10 of 19


Stopping criterion: the sAIS optimization procedure will not be terminated until the maximum number of iterations is reached.
The proposed sAIS optimization approach has a mechanism with the inside operation processes: 1. Somatic Recombination Somatic recombination is one of the gene rearrangement processes that involves cutting out small regions of DNA and then putting the remaining pieces of DNA back together in an error-prone way in the adaptive immune system.For every iteration of our proposed algorithm, we randomly create a new set of antibodies by using the somatic recombination process in order to search for optimized solutions;

Somatic Hypermutation
There are five types of immune mutations that we use in our proposed sAIS algorithm.Each type has a different function.If the cost of objection function of antibodies in IgM is larger than the cost of the old ones, we randomly selected one of four kinds of immune mutations to operate in the next step of mutation in the same iteration.In addition to the traditional four types of immune mutations, we proposed a new immune mutation, called IgG2, to further improve the performance of the sAIS optimizations:  IgM mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Inverse the sequence of suppliers or customers between i and j as shown in Figure 6.It is noted that if |i − j| < 2 it is not allowed to be mutated;


IgG2 mode: In the vehicle's route of sequence s, let i1, i2, j1, and j2 be randomly selected as four suppliers or customers in the sequence.Swap i1 to j1 and i2 to j2 in the neighbor of s as shown in Figure 8;


Stopping criterion: the sAIS optimization procedure will not be terminated until the maximum number of iterations is reached.
The proposed sAIS optimization approach has a mechanism with the inside operation processes: 1. Somatic Recombination Somatic recombination is one of the gene rearrangement processes that involves cutting out small regions of DNA and then putting the remaining pieces of DNA back together in an error-prone way in the adaptive immune system.For every iteration of our proposed algorithm, we randomly create a new set of antibodies by using the somatic recombination process in order to search for optimized solutions;

Somatic Hypermutation
There are five types of immune mutations that we use in our proposed sAIS algorithm.Each type has a different function.If the cost of objection function of antibodies in IgM is larger than the cost of the old ones, we randomly selected one of four kinds of immune mutations to operate in the next step of mutation in the same iteration.In addition to the traditional four types of immune mutations, we proposed a new immune mutation, called IgG2, to further improve the performance of the sAIS optimizations:  IgM mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Inverse the sequence of suppliers or customers between i and j as shown in Figure 6.It is noted that if |i − j| < 2 it is not allowed to be mutated;


IgG mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Swap these two suppliers or customers in the neighbor of s;  IgA mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Insert the suppliers or customers i to the sequence's jth position in the neighbor of s as shown in Figure 7;


IgG2 mode: In the vehicle's route of sequence s, let i1, i2, j1, and j2 be randomly selected as four suppliers or customers in the sequence.Swap i1 to j1 and i2 to j2 in the neighbor of s as shown in Figure 8;  IgE mode: Both the IgG mode and the IgA mode are applied in this mutation.In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.A new sequence s' is provided by swapping i and j.Then, in the sequence s' randomly selected i' and j' will be suppliers and customers, respectively, and tion processes: 1. Somatic Recombination Somatic recombination is one of the gene rearrangement processes that involves cutting out small regions of DNA and then putting the remaining pieces of DNA back together in an error-prone way in the adaptive immune system.For every iteration of our proposed algorithm, we randomly create a new set of antibodies by using the somatic recombination process in order to search for optimized solutions;

Somatic Hypermutation
There are five types of immune mutations that we use in our proposed sAIS algorithm.Each type has a different function.If the cost of objection function of antibodies in IgM is larger than the cost of the old ones, we randomly selected one of four kinds of immune mutations to operate in the next step of mutation in the same iteration.In addition to the traditional four types of immune mutations, we proposed a new immune mutation, called IgG2, to further improve the performance of the sAIS optimizations:  IgM mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Inverse the sequence of suppliers or customers between i and j as shown in Figure 6.It is noted that if |i − j| < 2 it is not allowed to be mutated;  IgG mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Swap these two suppliers or customers in the neighbor of s;  IgA mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Insert the suppliers or customers i to the sequence's jth position in the neighbor of s as shown in Figure 7;


IgG2 mode: In the vehicle's route of sequence s, let i1, i2, j1, and j2 be randomly selected as four suppliers or customers in the sequence.Swap i1 to j1 and i2 to j2 in the neighbor of s as shown in Figure 8;  IgE mode: Both the IgG mode and the IgA mode are applied in this mutation.In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.A new sequence s' is provided by swapping i and j.Then, in the sequence s' randomly selected i' and j' will be suppliers and customers, respectively, and

Affinity Maturation
Affinity means the total cost of the vehicles routed from the objection function of the VRPCD model.Some mutant antibodies with lower costs may be generated after several iterations of somatic hypermutation.Then, from iteration to iteration, the cost of producing antibodies will steadily decrease.This phenomenon is called affinity maturation;

Elimination Process
The antibodies will be stored in a limited list per iteration.Except for the best antibody, the antibodies that were produced and mutated per iteration were deleted from previous iterations.As a result, antibodies will search for alternative solution spaces for the VRPCD;

Escape Criterion
We designed an escape criterion for searching the antibody in the global optimization to prevent the search from entering the same solution space and remaining in the local optimum.
Figure 9 and Algorithm 1 show the proposed procedure of the sAIS algorithm.
Algorithm 1: The proposed sAIS algorithm.
1. Initialization: Create a population of p antibodies randomly (p is the size of the antibody population) 2.
For each iteration do: 3.
For each antibody do: 4.

Somatic recombination 5.
Calculate affinity of the antibody 6.
Decode the new string 9.
Decode the new string 14.
If the objective function (new) < the objective function (old): 15.
Eliminate antibodies except the best antibody in the population 19.
Reserve the best antibody in the population to affect velocity 20.
If count attain list length: 23.
Generate the new antibodies randomly 24.
End For Loop 25.

Escape Criterion
We designed an escape criterion for searching the antibody in the global optim to prevent the search from entering the same solution space and remaining in th optimum.
Figure 9 and Algorithm 1 show the proposed procedure of the sAIS algorithm

Computational Experiments
First, a preliminary investigation of the optimal parameter settings for the algorithm is presented here.Next, 60 VRPPD benchmark problems were retrieved from the VRP Web (http://www.bernabe.dorronsoro.es/vrp/accessed on 3 February 2023) as the criteria to evaluate the performance of the proposed heuristics.The comparison set is based on the GA algorithm, which is functioned by the customized software Palisade Evolver (https://www.palisade.com/evolver/)industrial edition 5.5.1, the Genetic Algorithm Solver plug-in for Microsoft Excel 2016.

Preliminary Tests
In the preliminary experiment, the design of experiments (DOE) was conducted to decide the optimal parameter settings for the sAIS optimization.The size of the population and escape list were acquired through the experimentation methodology.One randomly selected instance from the 60 testing problems was used to execute the experiment methodology.The parameter of population size, after an initial trial-and-error process, was located in the interval of [100, 500], and we used the set of {100, 200, 300, 400, and 500} to run the single factor factorial design experiment at 5 levels.Figure 10 shows the mean plot of the population size factor from the results of the experiments.The variable y is the objective function's minimized value (fitness value), as is the y-axis in Figures 10 and 11.

FOR PEER REVIEW 13 of 19
to run the single factor factorial design experiment at 5 levels.Figure 10 shows the mean plot of the population size factor from the results of the experiments.The variable y is the objective functionʹs minimized value (fitness value), as is the y-axis in Figures 10 and 11.When the length of the escape list is attained, the AIS algorithm will jump out of the present field.The DOE also determined the length from large numbers in the interval [300 700], and we chose {300, 400, 500, 600, and 700} as the 5 levels for the single factor factoria design experiment.Figure 8 shows the mean plot of the escape list factor from the results of the experiment.to run the single factor factorial design experiment at 5 levels.Figure 10 shows the mean plot of the population size factor from the results of the experiments.The variable y is the objective functionʹs minimized value (fitness value), as is the y-axis in Figures 10 and 11  When the length of the escape list is attained, the AIS algorithm will jump out of the present field.The DOE also determined the length from large numbers in the interval [300 700], and we chose {300, 400, 500, 600, and 700} as the 5 levels for the single factor factoria design experiment.Figure 8 shows the mean plot of the escape list factor from the results of the experiment.It is found that the sAIS method can generate a better-quality solution when the population size and escape list factors are assigned to 300 and 500.According to the literature review, the maximum iteration was set at 100,000 trials as the termination criterion in our When the length of the escape list is attained, the AIS algorithm will jump out of the present field.The DOE also determined the length from large numbers in the interval [300, 700], and we chose {300, 400, 500, 600, and 700} as the 5 levels for the single factor factorial design experiment.Figure 8 shows the mean plot of the escape list factor from the results of the experiment.
It is found that the sAIS method can generate a better-quality solution when the population size and escape list factors are assigned to 300 and 500.According to the literature review, the maximum iteration was set at 100,000 trials as the termination criterion in our experiments.These settings were used in the following experiment for the sAIS algorithm.On the other hand, the population size of genes in the GA method was set at 50; the mutation rate is 0.06; and the crossover rate is 0.15.The maximum iteration number was 5000 trials with 30 replications for each instance.

Experiment Results
A total of 60 benchmark datasets with homogeneous demands were used in our experiment.Each instance is replicated for 30 runs for both the sAIS algorithm and the GA method.The solution from the sweep method, which generates proper clusters for shipment grouping of vehicles, is then picked up as the initial solution for both the sAIS and the GA methods.The comparison of performance between the two methods is computed based on the average improvement rate (AIR), as shown in Equation (14) As listed, the 60 benchmark problems were sorted into 20 subgroups with three instances in each, based on the node coordinates assigned.Instances within one group are given identical node coordinates, except for the position where the cross-docking facility is located.Each set of the 60 instances contains 100 nodes, comprising both pickup and delivery customers.Each node is associated with a known demand of 10 units, while each vehicle dispatched is constrained by a capacity limit of 100 units.Examples of optimized routing sequences are shown in the Appendix A.
The computational results showed that the solution quality of the sAIS method is superior to the GA, outperforming each with an average improvement rate and lowest improvement rate of 11.98% and 9.47%, respectively, on the basis of average solution quality.In addition, the sAIS optimization, which has the lowest cost, was able to discover new, better solutions than the GA did in all 60 benchmarks, even though the computation time of the sAIS algorithm method takes longer; however, it is still a reasonable amount of time.Moreover, the maximum rate of improvement on average solution quality is 30.03%,whereas the minimum rate is acquired at 1.28%, indicating that the performance of the sAIS method is robust and competitive with the GA method.Nevertheless, all problems' average improvement rate was better able to discover new solutions than the GA did.
Finally, a one-sided, two-sample t-test was conducted to verify the performance of the two methods.The hypothesis test is: H 0 : µ GA − µ sAIS = 0 H A : µ GA − µ sAIS > 0 the results showed that the t-value = 2.09 and the p-value = 0.019.Because the p-value is less than the 0.05 significant level, we reject the null hypothesis and conclude that µ GA is significantly greater than µ sAIS .That is, the total cost, including transportation costs and operational costs, generated by the sAIS approach is smaller than the total cost generated by the GA in our experiments.

Conclusions and Future Research
In this research, a novel sAIS algorithm is proposed to approach the combinatorial optimal solution of the VRPCD.The primary objectives of this work include the integration of the operation of cross-docking and the optimal vehicle routing schedule into the design of supply chain optimization.A significant development lies in the synchronization between upstream suppliers and downstream retailers, where both sides of the supply chain are simultaneously considered to collaborate on the physical flow of goods in each inbound and outbound process.Manufacturers can reduce logistics costs by building IoT monitoring technology into smart logistics IT systems to track its location in real-time, reduce inventory levels in warehouses, monitor environmental conditions, and optimize the routing sequence for trucks for smart city infrastructure.
The computational results show that the sAIS model is effective for solving the VRPCD.The effectiveness of the method comes from the two-phase mechanism.In the initial route generation phase, the initial solution was generated by the sweep method before being input into the route's optimization phase with the AIS algorithm.The combination of the two-phase approach ensures that the sAIS method yields quality solutions.
A total of 60 benchmark problems were deployed to investigate the applicability of the proposed sAIS method.The experimental results showed that the sAIS method was able to produce significant improvements over the GA, surpassing each testing problem.Additionally, the sAIS method was able to discover better solutions than the GA method for all 60 benchmark problems, and the sAIS's search time is faster than the GA.It found sAIS to be a useful methodology.As artificial intelligence research rapidly increases, more simulations from human body systems could be used to improve the solution quality of the AIS optimization, and the AIS optimization could be applied to other research problems in the operations research field, such as multiple-criteria decision-making problems or supplier selection problems, etc.

Appendix A
Examples of optimized routing sequences are shown in the appendix.

Figure 1 .
Figure 1.A layout of a typical CD facility.

Figure 1 .
Figure 1.A layout of a typical CD facility.

Figure 2 .
Figure 2. The selection process of antibodies.

Figure 3 .
Figure 3.A simple illustration of the VRP routing network.

Figure 3 .
Figure 3.A simple illustration of the VRP routing network.

Figure 4 .
Figure 4.An illustration of clustering by the sweep method.

Figure 4 .
Figure 4.An illustration of clustering by the sweep method.

Figure 5 .
Figure 5.An illustration of 10 nodes with homogeneous demand served by 2 vehicles.
. The basic components of the sAIS applied in the VRPCD are shown as follows:  Antigen: the antigen mimics the equations of the VRPCD model, which represent the formulas for the objection function and constraint functions;  Antibody: The antibody is the candidate solution for the VRPCD model.The strings of solutions refer to the visiting sequence of suppliers or customers.The antibodies are feasible sets of models for the problem;  Cell: the cell represents the ID of suppliers or customers;  Population size: the number of antibodies (feasible solutions of the VRPCD model) that are used per iteration;  Number of clones for each cell: number of times antibodies are reproduced per iter-

Figure 5 .
Figure 5.An illustration of 10 nodes with homogeneous demand served by 2 vehicles.

Figure 6 .Figure 7 .
Figure 6.The operation of IgM mode.IgGmode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Swap these two suppliers or customers in the neighbor of s;  IgA mode: In the vehicle's route of sequence s, randomly select two suppliers or customers, i and j.Insert the suppliers or customers i to the sequence's jth position in the neighbor of s as shown in Figure7;

Figure 6 .
Figure 6.The operation of IgM mode.

Figure 6 .
Figure 6.The operation of IgM mode.

Figure 7 .
Figure 7.The operation of IgA mode.

Figure 7 .
Figure 7.The operation of IgA mode.

Figure 6 .
Figure 6.The operation of IgM mode.

Figure 7 .
Figure 7.The operation of IgA mode.

Figure 9 .
Figure 9.The processes of the proposed sAIS flowchart.

Figure 9 .
Figure 9.The processes of the proposed sAIS flowchart.

Figure 10 .
Figure 10.The mean plot of population size for a one-factor factorial design.

Figure 10 .
Figure 10.The mean plot of population size for a one-factor factorial design. .

Figure 10 .
Figure 10.The mean plot of population size for a one-factor factorial design.

Figure 11 .
Figure 11.The mean plot of the escape list for a one-factor factorial design.

Figure 11 .
Figure 11.The mean plot of the escape list for a one-factor factorial design.

Figure A2 .
Figure A2.Optimized delivery routes from sAIS for instance 53P1.

Figure A3 .
Figure A3.Optimized pickup and delivery routes from sAIS for instance 53P1.

Figure A2 .
Figure A2.Optimized delivery routes from sAIS for instance 53P1.

Figure A3 .
Figure A3.Optimized pickup and delivery routes from sAIS for instance 53P1.

Figure A2 .
Figure A2.Optimized delivery routes from sAIS for instance 53P1.

Figure A2 .
Figure A2.Optimized delivery routes from sAIS for instance 53P1.

Figure A3 .
Figure A3.Optimized pickup and delivery routes from sAIS for instance 53P1.

Figure A3 .
Figure A3.Optimized pickup and delivery routes from sAIS for instance 53P1.

Table 1 .
. The experimentation results are given in Table1.Performance comparisons between sAIS and GA.