Research on Green Vehicle Path Planning of AGVs with Simultaneous Pickup and Delivery in Intelligent Workshop

: In this study, we present and discuss a variant of the classic vehicle routing problem (VRP), the green automated guided vehicle (AGV) routing problem, which involves simultaneous pickup and delivery with time windows (GVRPSPDTW) in an intelligent workshop. The research object is AGV energy consumption. First, we conduct a comprehensive analysis of the mechanical forces present during AGV transportation and evaluate the overall operational eﬃciency of the workshop. Then, we construct a green vehicle path planning model to minimize the energy consumption during AGV transportation and the standby period. Hence, the problems considered in this study are modeled based on asymmetry, making the problem solving more complex. We also design a hybrid diﬀerential evolution algorithm based on large neighborhood search (DE-LNS) to increase the local search ability of the algorithm. To enhance the optimal quality of solutions, we design an adaptive scaling factor and use the squirrel migration operator to optimize the population. Last, extensive computational experiments, which are generated from the VRPSPDTW instances set and a real case of an intelligent workshop, are designed to evaluate and demonstrate the eﬃciency and eﬀectiveness of the proposed model and algorithm. The experimental results show that DE-LNS yields competitive results, compared to advanced heuristic algorithms. The eﬀectiveness and applicability of the proposed algorithm are veriﬁed. Additionally, the proposed model demonstrates signiﬁcant energy-saving potential in workshop logistics. It can optimize energy consumption by 15.3% compared with the traditional VRPSPDTW model. Consequently, the model proposed in this research carries substantial implications for minimizing total energy consumption costs and exhibits promising prospects for real-world application in intelligent workshops.


Introduction
With the rapid development of Industry 4.0, the traditional forms of industrial work are facing unprecedented challenges. The construction of an intelligent workshop becomes an important starting point and a technical way for the manufacturing industry to promote digital transformation. Its core content is to improve the overall level of product production. As an important auxiliary link in the workshop, material distribution occupies a large proportion of the operation cost [1,2]. As one of the important carriers of workshop logistics, the material handling system composed of AGVs can greatly improve the production flexibility and competitiveness of enterprises.
According to a report by the International Energy Agency (IEA), the carbon dioxide emissions of the EU power industry increased by 28 million tons in 2022. However, fossil fuels, especially coal resources, support nearly half of the electricity demand. The production process of green manufacturing, energy saving, and emission reduction has become the main development trend of the modern manufacturing industry [3]. Its purpose is to achieve the goal of net zero emissions in 2050. When the AGV performs handling tasks, AGV is powered by its own ba ery, which needs to be supplemented after 8-10 h of operation [4]. Therefore, AGV execution tasks may be interrupted due to insufficient power, which will cause the production plan to be disrupted. Therefore, it is very important for the intelligent workshop to reasonably arrange and plan the path of reducing energy consumption for AGV. It can not only save costs for enterprises to achieve the development goal of energy conservation and emission reduction, but also improve the stability of material handling systems in intelligent workshops.
There are periodic logistics transportation tasks in the intelligent workshop. When the task information is released, AGVs start from the production scheduling center. They perform simultaneous pickup and delivery logistics transportation tasks on the workstation of the workshop along the planned route and return to the production scheduling center after completing the tasks. VRPSPD of AGV in intelligent workshops can be regarded as a variant of the typical vehicle routing problem (VRP). In many studies of the current green vehicle routing problem, travel distance is no longer the only influencing factor. Energy consumption or fuel consumption is closely related to the real-time payload weight of the vehicle body [5][6][7]. The AGV load is always changing during simultaneous pickup and delivery. The traction force of AGV will change with the change of load, in the process of driving. Therefore, we add the consideration of an AGV load change to the research of AGV green path planning in the intelligent workshop.
The production characteristics of an intelligent workshop are many batches and high output. If the AGV arrives at the workstation too early, there will be idle waiting causing energy waste. Moreover, the operator and the production machine work is suspended, resulting in production losses and delays in subsequent processing. We comprehensively consider the overall operation efficiency of the workshop and the development goal of energy conservation and further add time window constraints based on GVRPSPD. For this reason, G-VRPSPDTW (GVRPSPD With Time Windows) is the focus of this research.
Currently, the AGV path planning problem within intelligent workshops has garnered significant research a ention as a subfield of the VRP problem. Numerous successful methods have been employed to address AGV path planning challenges [8][9][10]. However, in the majority of these studies, researchers tend to prioritize cost optimizations such as travel distance, travel time, number of AGVs, and penalty costs [11]. Regre ably, energy-saving path planning for AGVs during simultaneous pickup and delivery processes is seldom explored. Leveraging insights from contemporary scholarly research, we formulate a green vehicle path planning model with the optimization objective of minimizing energy consumption during AGV transportation and standby periods. At the same time, we devise a hybrid differential evolution algorithm and subsequently analyze the experimental outcomes.
This research focuses on the sustainable development of the intelligent manufacturing industry and the greater a ention to the energy-saving and environmental protection application environment of the intelligent factory workshop. The energy-saving path planning of AGV simultaneous pickup and delivery is studied. We hope to minimize the energy consumption of the AGV fleet to complete the pickup and delivery tasks. In this research, our main contributions are summarized as follows:


We established a mathematical model by analyzing the mechanical force of AGV and the overall efficiency of the workshop. The model takes into account the real-time load of AGV, the rated power of AGV and the time window of workstation, so it is more in line with the actual application environment.  We propose a hybrid differential evolution algorithm. The algorithm combines the differential evolution algorithm with the large neighborhood search algorithm. The purpose is to increase the local search ability of the differential evolution algorithm, so that it can find the optimal solution faster, avoid falling into the local optimum, and improve the optimized solution.
 Our experiments show that the hybrid algorithm proposed in this paper is effective and applicable in solving the GVRPSPDTW problem of intelligent workshop. The model established in this paper can effectively reduce energy consumption compared with the traditional VRPSPDTW model.
The remainder of this article is structured as follows: Section 2 provides a review of green vehicle path planning and vehicle routing problem with simultaneous pickup and delivery with time windows. Section 3 introduces the problem, outlining its definition and mathematical formulation. In Section 4, we present the design of our differential evolution algorithm, incorporating novel enhancements such as squirrel search and a large neighborhood search. Section 5 outlines the experimental methodology utilized to validate the effectiveness and applicability of our algorithm, while also investigating the key factors impacting model optimization. Lastly, Section 6 offers a comprehensive conclusion and outlines potential avenues for future research endeavors.

Literature Review
The AGV path planning problem has been the subject of extensive research over the past few decades, and it remains a topic of great interest for scholars and manufacturing companies. However, it is worth noting that the GVRPSPDTW of AGVs in intelligent workshops has been li le explored in the existing literature. Therefore, we conducted a thorough review of the literature related to the green vehicle routing problem and the simultaneous pickup and delivery vehicle routing problem with time window constraints.

Green Vehicle Route Problem
The transportation and transportation industry accounts for a large proportion of global carbon emissions. Therefore, it is very important to study energy conservation and emission reduction and environmental pollution reduction in the transportation and transportation industry. According to a report from the International Energy Agency (IEA), the transport sector has generated about 25% of global CO2 emissions in recent years [12]. The concept of GVRP is proposed by considering the reduction of carbon dioxide emissions and the realization of green and energy-efficient logistics and transportation [13]. It is a branch of the classic VRP and belongs to the NP problem. The GVRP in the workshop is to find the best route set by scientifically planning the number of departures and transportation routes of the fleet to meet various constraints. the aim is to minimize energy consumption [14].
At present, the research on the green vehicle routing problem is mainly divided into two aspects: traditional fuel and new energy (electric energy). The research on GVRP of traditional fuel (fossil fuel)-driven vehicles are mostly concentrated in intercity logistics. The construction of the GVRP model is based on vehicle fuel consumption and carbon emissions as a measure, its goal is to reduce carbon dioxide emissions and improve environmental pollution and reduce energy consumption [15]. For example, Bektas [16] further extended the goal of traditional VRP, and he considered the problem of vehicle pollution emissions. The optimization objectives in this paper include vehicle cost, energy consumption, pollution gas emission and labor cost. Based on this, Demir [17] proposed an adaptive large neighborhood search algorithm to be er solve the GVRP problem. In this paper, the actual geographical data is used to analyze the algorithm, in order to have greater applicability and authenticity for GVRP research. Antonio Giallanza [18] constructed a multi-objective fuzzy chance-constrained programming model to minimize total cost and carbon emissions and confirmed the robustness of the model by applying it to the case study of Hoh Xil. In recent years, power resources instead of traditional fuels have gradually become the focus of GVRP research and applications. Amiri Afsane [19] studied the GVRP problem of heterogeneous vehicles (traditional trucks and electric trucks). Transportation costs and greenhouse gases are considered combinatorial optimization objectives. Surendra Reddy Kancharla [20] considered the influence of load on ba ery consumption to calculate energy demand. He proposed an adaptive large neighborhood search algorithm with special operators to solve energy consumption and verified the effectiveness of the algorithm by testing on benchmark instances.
The GVRP problem in the workshop is to study the energy consumption of AGVs. Most AGVs used in intelligent workshops are also powered by ba eries. At present, the research direction of scholars on AGV energy consumption is mainly divided into two aspects: One is to take the power demand of the AGV motor and its motion state as the research object [21,22]. With the increase in production demand in the intelligent workshop, AGVs need to complete deliveries of the workstations at the same time. The realtime change in AGV's load weight caused by AGV loading and unloading of workstations has become another research focus. Some professional institutions and academic researchers have obtained a large number of statistical data on AGV energy consumption and its loading weight. These statistics show that these two variables are proportional [23]. Zhou [24] analyzed AGV energy consumption from AGV loading weight and travel distance and then constructed a multi-objective mathematical model to minimize total energy consumption and total weighted delay. He proposed a tabu-enhanced particle swarm optimization (TEPSO) algorithm to solve multi-objective problems and verified the effectiveness of the algorithm. Wang [25] took the body load and AGV driving speed during AGV transportation as the optimization objectives and established the AGV driving energy consumption model. He evaluated and weighed the two optimization objectives of distribution cost and energy consumption by a weighted index.
According to the above research, it can be concluded that the vehicle energy consumption index is closely related to the travel distance and the real-time load weight of the body. The reason is that they are the main parameters of the mechanical force of AGVs. In addition, the current research on GVRP mainly focuses on the path planning and vehicle scheduling of intercity logistics. Although some experts and scholars have studied the energy consumption of AGVs, the energy consumption characteristics of AGVs' actual logistics transportation process and the overall operation efficiency of the workshop are rarely fully considered. For example, [21] considered the AGV motion state and time window constraints, but did not consider the impact of AGV load on energy consumption. Similarly, [24] considered the influence of AGV load and driving distance on energy consumption but did not consider the influence of workstation time window constraints on energy consumption. The environment of GVRP in intelligent workshops is different from that of intercity logistics. At present, the VRPSPDTW problem in the intelligent workshop has not been studied. Therefore, we propose a model, which considers two major factors of AGV energy consumption. namely, the main mechanical force and the overall operating efficiency of the workshop. In terms of mechanical force, we study its influencing factors including travel distance, real-time body load and AGV rated power. In terms of operating efficiency, we increase the workstation service time window constraint. The model can directly calculate the energy consumption of AGV performing the entire logistics transportation tasks in the intelligent workshop.

Vehicle Routing Problem with Simultaneous Pickup and Delivery with Time Windows
VRPSPD with a time window, namely VRPSPDTW, can be regarded as a combination of the VRPTW problem and VRPSPD problem. VRPSPDTW is widely used in modern logistics [26], therefore, it may be the most studied problem in the variants of VRPSPD. The production scheduling system usually limits the loading and unloading time of each workstation at a specific time interval in the intelligent workshop. Due to the actual correlation and inherent complexity of VRPSPDTW, it is a hot research direction of algorithm innovation and practical path planning in the field of logistics.
In the existing research, many path search algorithms are proposed based on mathematical models and environmental representations. They can be divided into two categories: exact algorithms and heuristic algorithms. Angelelli [27] was the first and so far the only scholar to solve the VRPSPDTW with an exact algorithm. The exact algorithm he designed was mainly the classical branch pricing method, but due to the limitations of the exact algorithm, the algorithm was only tested on 20 customers. Cao and Lai [28] first a empted to solve VRPSPDTW by implementing heuristics. They proposed an improved genetic algorithm, which contains an encoding for the customer order problem. The experimental results show that their proposed algorithm can effectively find the optimal or approximate optimal solution when the number of customers is eight. Later, in order to further solve the VRPSPDTW problem, Lai and Cao [29] designed an IDE. In the experiment, they used a case with eight customers and two vehicles. Wang and Chen [30] studied VRPSPDTW using co-evolutionary genetic algorithm. They designed a benchmark instance for solving VRPSPDTW based on the famous Solomon [31] to evaluate the performance of their proposed algorithm. Recently, Hof and Schneider [32] designed an adaptive large neighborhood search algorithm (ALNS-PR) to solve VRPSPD problems, including VRPSPDTW problems. In their published results, it can be seen that ALNS-PR can obtain a be er solution than the previous benchmark results, and the effectiveness of the algorithm is verified. Wu Hongguang [33] proposed an ant colony optimization algorithm with a damage and repair strategy (ACO-DR) to solve the VRPSPDTW problem and verified the effectiveness of the algorithm.
In the past few decades, VRPSPDTW has been widely studied by scholars. However, as we know, the current research on VRPSPDTW in intelligent workshops is scant in the literature. Therefore, we reviewed the literature on AGV pickup and delivery in the workshop. Sungbum Jun [34] took minimizing the total delay of transportation requests as the optimization goal in the pickup and delivery problem of intelligent workshops. He proposed two constructive heuristic algorithms with high computational speed. The proposed algorithm can find a near-optimal solution in a reasonable time. Qiu [35] studied the routing problem of the AGV pickup and delivery process in the depot. He took minimizing energy consumption as the optimal goal and adopted the particle swarm optimization (PSO) algorithm based on the rapidity of the algorithm. The simulation results showed that the model has good performance and can effectively reduce the energy consumption of the system. Zou [11] established a multi-objective mixed integer linear programming model. Customer satisfaction and total driving cost were the optimization objectives of AGV pickup and delivery problems. He developed an effective multi-objective evolutionary algorithm to solve this problem.
From the brief review above, we can observe that many methods have been successfully applied to the pickup and delivery problems of AGVs in intelligent workshops. Most experts and scholars use heuristic algorithms to solve such problems. Due to the characteristics of the VRPSPDTW problem, we usually abstract the logistics transportation network in the intelligent workshop into a topological graph to plan the path and find the optimal path. With the increase of the difficulty of the path planning problem and the increase of the constraints in the corresponding model, it is a very difficult task to find the optimal path to get the optimal solution. For this reason, On the basis of the DE algorithm, we propose an effective hybrid differential evolution algorithm (DE-LNS) to solve the VRPSPDTW problem in intelligent workshops.
The comparison of the literature is shown in Table 1.

Mathematical Model
A typical intelligent workshop is mainly composed of some workstations, material handling equipment (AGV), and a production scheduling center. In general, the actual layout of the workshop is asymmetric. This also increases the difficulty of problem solving. The workstation is mainly composed of two parts, one is the buffer area, and the other is the computer numerical control (CNC). The buffer is used to store raw materials and finished products. The continuous work of CNC machine tools will reduce the inventory of raw materials in the buffer area and increase the inventory of finished products, the production scheduling center formulates the AGV scheduling plan, then AGVs perform logistics transportation tasks. Each AGV starts from the production scheduling center along the route to the designated workstation for pickup and delivery as shown in Figure  1. To ensure the overall operation efficiency in the workshop, an AGV is required to not arrive at the workstation at the latest time beyond the specified time window and return to the production scheduling center after each transportation cycle. The problem is how to assign workstations to AGVs and how to determine the order of workstations in each AGV. The objective is to minimize the energy consumption of AGVs in performing logistics transportation tasks, including AGV driving energy consumption and standby energy consumption. To effectively illustrate the problems studied, we made the following basic assumptions:


The production inside the workshop is stable, and there will be no shutdown, failure, or other phenomena.  The location of the workstation and the production scheduling center in the workshop is unchanged.  The type of AGV in the workshop is the same, and its loading capacity is consistent.  Before the dispatch command is issued, all AGVs are parked at the production scheduling center. When AGV performs logistics transportation tasks, there may be some traffic problems. In this study, the infrared obstacle avoidance module is installed on the AGV body, which can prevent the occurrence of traffic problems in the workshop.  The speed of AGVs in turning and going straight is a uniform constant.

Notation
The summary of the notations used in this paper is shown in Table 2.

Energy Consumption Analysis
With the above description, we define the workshop GVRPSPDTW as an undirected The node 0 in the set is represented as the production scheduling center, which is the starting point of AGV. 0 {1,2, } V n   is the workstation loading and unloading point. A represents the set of all edges connecting each node, Suppose that an AGV travels from node i to node j the traveling distance and traveling time are respectively calculated by using the following formula: If there is an obstacle between i and j , the distance between two points is calculated by using the Dijkstra algorithm.
The arrival time of node j is calculated by using the following formula: When AGVs perform logistics transportation tasks in the intelligent workshop, they roll on a flat, concrete surface, such as that commonly encountered in factories. As discussed in Garcia [36], The energy consumption of such vehicles is mainly affected by three mechanical forces, namely traction ( traction F ), friction ( friction F ) and air resistance ( airdrag F ), in Newton (N). Because the vehicle itself has a higher quality and the vehicle's speed is lower, the air resistance is negligible, which means that the vehicle's aerodynamics will not play a significant role in this AGV's logistics transportation tasks. There is traction when there is movement between the object and the tangential surface. AGV traction for uniform motion depends on its rated power and travel speed, as shown in Formula (4). The friction F is expressed as the reaction force that hinders the movement of the object when it rolls on the surface. According to Formula (5), friction F varies with the weight load, rolling coefficient, and gravity acceleration of an AGV.
The research in [37] shows that any plan to minimize energy consumption should take into account the three parameters of load, number of stops, and travel distance. Because of the characteristics of the problem studied in this paper, an AGV performs pickup and delivery on all workstations in the workshop, and the number of stops is fixed. But the AGV load is constantly changing due to the nature of the VRPSPD problem. In order to solve the energy consumption change caused by load, they use a fuel consumption function (FCR) where the fuel consumption rate per unit volume is proportional to the driving distance and vehicle weight.
In this formula, 0 Q is the weight of the vehicle without load, 1 Q is the load carried by the vehicle, and b is the automobile fuel consumption coefficient.
Let Q be the capacity of the vehicle. The FCR when fully loaded ( * p ), and the FCR without load ( 0 p ) are calculated by Equations (7) and (8), respectively.
the mixed integer programming Formula (6) of GVRPSPD is constructed in [5].
In the formula, 0 c is the unit fuel cost.  is a parameter.
V are the total demand for delivery and pickup on an arc.
The model in [5] takes into account the real-time change in load when calculating energy consumption. However, the environment of intercity logistics is still different from that of an intelligent workshop. The author only models the path problem of one vehicle, and multiple AGVs perform tasks in the intelligent workshop.
For this reason, we combine the energy consumption calculation method proposed in the literature [5,37] with the GVRPSPDTW in the workshop and define the energy consumption p E of the AGV transportation process as: In the case of AGVs working normally, the sum of speed and power can be assumed to be known and constant because they are automatically managed by the production scheduling center. The rolling coefficient r C is also constant, as plants generally have a concrete floor. When AGV performs logistics transportation tasks, the parameters affecting AGV energy consumption are real-time changing load and driving distance.
Due to the limited ba ery capacity of the AGV body, if the AGV arrives at the workstation in advance, there will be a certain standby window period. Excessive standby time will not only cause the punctuality rate of AGV to performing the loading task to decrease, but also may lead to excessive consumption of AGV ba ery power and inability to support the completion of the remaining tasks. Conversely, the workstation cannot work because the cache is overloaded. Based on the above analysis, we define the service time window of the workstation as a semi-soft time window to ensure the comprehensiveness of the AGV energy consumption calculation. We allow an AGV to arrive at the workstation earlier but do not allow an AGV to reach the workstation beyond the established limit. The standby energy consumption t E is defined as follows:

Mathematical Model
The problem solved in this paper is how to solve the service order of each AGV arriving at the specified workstation to perform logistics transportation tasks. The output is also the order of AGV service to the workstation, and the variable is an integer. We hope to achieve the purpose of reducing energy consumption by changing the order of AGV access. For this reason, the objective of this study is to lower the energy consumption. The mathematical model was developed as follows: The constraints are as follows: In the above model, goal (12) shows the solution goal of the problem, which is to minimize the total energy consumption. It is separated into two sections: AGV transportation energy consumption and standby power consumption. Constraint (14) ensures that each workstation is only accessed by one AGV and only accessed once. Constraint (15) states that backward paths and possible loops are not allowed. Constraint (16) imposes that each AGV starts from the production scheduling center and returns back to it. Constraint (17) limits the number of AGVs used. Constraint (18) ensures that the service time of AGVs to workstation i cannot exceed the i LT of i point. i.e., the time window constraints. Constraint (19) ensures that the load weight of any workstation of the AGV transportation route does not exceed the maximum load of the body. Constraints (20)-(23) determines the upper and lower limits of pickup and delivery on arc i j  . Constraints (24)-(25) deal with the balance of pickup and delivery demands on the route and the amount of load carried on arcs. Constraint (26) is a decision variable constraint.

Algorithm Design
We design a differential evolution algorithm with a hybrid large neighborhood search (DE-LNS) for the complexity of GVRPSPDTW. The overall framework of DE-LNS is shown in Algorithm 1. return final external archive

Algorithm 1 General Framework of DE-LNS
The algorithm uses a differential evolution algorithm as a global search algorithm. Firstly, we design a mutation strategy with an adaptive scaling factor and an adaptive crossover operator. Then we refer to the foraging process of the squirrel search algorithm [38] to optimize and screen the population. Finally, we introduce a large neighborhood search as a search mechanism to perform local searches on excellent chromosomes. The large neighborhood search algorithm has a stronger ability in the local search stage and can overcome the shortcomings of the differential evolution algorithm. Among them, P is the population size, MaxGen is the maximum number of iterations, D is the optimal number of individuals in the population, MaxN is the number of iterations of large neighborhood search, and f is the fitness function.

Encoding and Decoding Design
According to the characteristics of GVRPSDPTW, the chromosome coding in this algorithm adopts the natural number coding mechanism. The total number of genes on one chromosome in the population is n m  . The workstation number is {1,2, , 1, } n n   , and the AGV number is   1, 2, , n n n m     . We make the chromosome gene composition multiple gene fragments by inserting the AGV number in the workstation number. Each gene fragment represents a path of the AGV fleet in the logistics transportation task.
We give a simple example to solve the problem under consideration to express our method clearly. When

Initialization Population
We combine the initial population with random traversal to generate the initial path according to the nature of the simultaneous pickup and delivery problem. We constrain each route according to the service time window and AGV loading status. The specific steps of Algorithm 2 are shown. Among them, N is the number of workstations, 0 V is the workstation set, Q is the maximum load of an AGV, and the left time window of ET workstation.

Mutation Operator
The mutation is to obtain a new solution different from the current population. The value of the scaling factor F is crucial to the algorithm in the mutation strategy. If the value of F is large, it is easy to obtain the optimal solution, but the convergence speed of the algorithm will decrease. Conversely, the convergence speed of the algorithm becomes faster, but it is prone to stagnation. Therefore, the scaling factor F should maintain a large value at the beginning and then change to a smaller value. We refer to the scaling factor designed by Zhou [39] and design an adaptive scaling factor based on the logistic function. The specific formula of scaling factor F is: start F in the formula is the initial value, In the early stage of the algorithm, the scaling factor number is large, and the purpose is to maintain the diversity of the population. The scaling factor is gradually reduced, and the convergence speed and local optimization ability are improved in the iterative process, as shown in Figure 3. The value of F is in a reasonable range during the iteration process. When the value of start F is different, the value range of F is different. The larger the value range is, the more conducive it is to enhancing the diversity of individuals. In order to ensure the comprehensiveness of F in the iterative process, the initial value of start F is 0.9.
We use Formula (28) to mutate the current individual according to the scaling factor obtained from each generation.
In the formula, First, if the gene g in the mutated individual is greater than P , then g g P   ; if g is less than 0, then g g P   . Then, the same gene is deleted. Finally, the missing gene fragment is completed. A possible mutation process is shown in the Figure 4.

Crossover Operator
The crossover operator is similar to the mutation operator. We also use an adaptive crossover operator in the algorithm. However, the crossover operator is opposite to the scaling factor, and the crossover operator needs to increase with iteration. The smaller crossover operator can avoid destroying the mutant individual. A larger crossover operator can avoid the algorithm from falling into the local optimum in the later stage of evolution and can increase the population diversity. as shown in the Figure 5. The value of CR is in a reasonable range during the iteration process. When the value of start CR is different, the value range of CR is different. The larger the value range is, the more conducive it is to enhancing the diversity of individuals. In order to ensure the comprehensiveness of CR in the iterative process, the initial value of start CR is 0.9. Therefore, the self-adaptive crossover operator is shown as follows:

1
( 1) The start CR in the formula is the initial value,  We use the OX crossover operator to expand the search range of the solution space. The purpose is to generate different AGV task allocation schemes. As shown in Figure 6, taking the generation of offspring 2 as an example, we randomly selected the crossover position of two chromosomes to perform the crossover operation. The sequence 4-5-6-7-8 in parent 1 is the first segment of offspring 2. The sequence of parent 2 after screening points 4, 5, 6, 7 and 8 was used as the second segment of offspring 2. Similarly, offspring 1 is generated.

Squirrel Migration Operator
The squirrel migration operator designed in this paper is based on the original squirrel search algorithm (SSA). We divide all individuals in the population into three different categories of squirrels according to their fitness function. The best squirrel remains unchanged; it is the best individual in the current population. The number of sub-optimal squirrels accounted for 10% of the total. Sub-optimal squirrels use one of the best squirrels as a model, and the OX crossover operation is performed between them. The remaining number of individuals are common squirrels. In the second case, common squirrels used one of the sub-optimal squirrels as a model [40], and the OX crossover operation is performed between them. We perform crossovers on different species of squirrels. Algorithm 3 shows the migration process. The squirrel migration operator of the SSA algorithm is used to optimize the population solved by the DE algorithm. At the same time, the migration operation can filter the population and retain some of the optimal solution sets for large neighborhood search operations. return The top 25% individuals

Large Neighborhood Search
The differential evolution algorithm easily falls into local optimal solutions in the process of solving. We use the large neighborhood search algorithm to perform local search operations on the filtered excellent chromosomes. The purpose is to improve the quality of the solution. For the consideration of the characteristics of the problem, we use the method of removing and inserting to meet all constraints. The point is to destroy and recreate some parts of the solution set to search for new feasible solutions. This process is to provide a be er mechanism to escape the local optimal solution created by DE.

Removal Operator
For the AGV path solution set of the GVRPSPDTW problem in the intelligent workshop, we use the Shaw removal (SR) method. The general idea of the SR operator is to remove multiple sets of points. The similarity between two workstations i and j is defined by the correlation weight ( , ) w i j . This includes two items: distance item ij C and relation item ij R . We use the distance ij c between workstations i and j to define their correlation ij C . The lower the ij C , the greater the relationship between the two workstations. The determination of ij R is the same as that of ij C , because the relationship between the two workstations i and j on the same transport path is greater than the relationship on different paths.
The definition of ( , ) w i j is as follows: i and j are on the same path, the value is 1, otherwise it is 0.

Insertion Operator
The greedy repair operator is the most commonly used insertion operator method. However, the greedy insert operators only consider improvements resulting from inserting a contract into its best position. This is its obvious disadvantage, and it does not consider the impact of this insertion on subsequent decisions. The regret insertion operator tries to circumvent this problem [41,42]. Let  . We consider a regret Regret-2 operator [43] in our study.
The first inserted workstation i satisfies: We insert workstation i at the best position, and the process iterates until the required number of workstations is inserted. The flow-chart of the DE-LNS algorithm structure can be seen in Figure 7.

Experimental Studies
The algorithm proposed in this paper is coded by MATLAB language and was tested on a PC with Intel Core (TM) i7-8700 3.20 GHz CPU, 16 GB RAM, and Windows 10 OS. GVRPSPDTW in the intelligent workshop is an important extension of VRPSPDTW. Therefore, the computational experiment part of this paper is divided into two stages, respectively in terms of algorithm performance effect and model energy saving effect. In the analysis of algorithm effectiveness, we use the international benchmark sets of VRP-SPDTW proposed by Wang and Chen [30] to test the efficiency and performance of the DE-LNS algorithm. At the same time, we use the data of the real workshop to conduct experiments to measure the applicability of the algorithm in solving the GVRPSPDTW of the intelligent workshop. In the energy saving effect, we conduct simulation experiments with the traditional VRPSPDTW model, and compare the calculated energy consumption to verify the energy-saving performance of the model.

Instances and Results Comparison
The international benchmark sets of VRPSPDTW contain 9 small-scale instances and 56 medium-scale instances. It can be divided into categories 1 and 2 based on the width of the time window. Class 1 indicates that the time window required by the customer is narrow, while class 2 indicates that the time window is wide. It can also be divided into R, C, and RC categories based on the distribution of customer points. R class indicates that the geographical distribution of customer points is random and sca ered, the C class indicates that the geographical distribution of customer points is aggregated, and the RC class is a mixture of the above two types. All of the small-scale instances are RC category. The number of customer points is 10, 25, and 50, respectively. The medium-sized instances include 23 R categories, 17 C categories and 16 RC categories. The set classification is shown in Table 3: In this experiment, all these instances are defined in a two-dimensional Euclidean space. This means that all nodes (customer or depot) has their own unique coordinate ( , ) i i x y , and the distance between two nodes i and j is Euclidean distance The states in the VRPSPDTW example set can be further divided into six types (subsets) according to the distribution of customer locations and the strength of time window constraints and capacity constraints. Table 3. VRPSPDTW instance set classification.

1
The customer distribution is concentrated, the time window is narrow and the vehicle load is small.
The customer distribution is random, the time window is narrow and the vehicle load is small.
The customer distribution is mixed, the time window is narrow and the vehicle load is small.

2
The customer distribution is distributed, the time window is wide and the vehicle load is heavy.
The customer distribution is random, the time window is wide and the vehicle load is heavy.
The customer distribution is mixed, the time window is wide and the vehicle load is heavy.
In this study, we solve the above six sets of instances under three different customer scales. The purpose is to illustrate the effectiveness of DE-LNS in solving VRPSPDTW. To ensure the comprehensiveness and completeness of numerical experiments, we consider three state-of-the-art algorithms, including Co-GA [30], p-SA [44], and ALNS-PR [32].
Each of these algorithms is tested on some or all instances in the Wang and Chen set. We analyze the experimental mathematical results obtained by DE-LNS, as shown in Tables 4  and 5. We take the number of vehicles (NV) and total travel distance (TD) as the main optimization objectives. The gap in NVs is expressed as an integer, and the gap in TD is expressed as a percentage. Minimizing NV has a higher priority than minimizing TD in Wang and Chen's instances. For this reason, we refer to the objective function of Shengcai Liu [45] to solve the example, as shown in Formula (32).
The total cost of S is recorded as In this experiment, we use the value of TC as the evaluation index of the effectiveness of the algorithm. If TC of one solution is less than that of another, then the solution is be er than the other, and the algorithm is be er. Table 4 lists the results of small-scale instances. The table includes the commercial mathematical software CPLEX (which was used in [30] to find the optimal solution for small-scale cases), the optimal experimental results of CO-GA, P-SA, and ALNS-PR, and the average and optimal values of the DE-LNS algorithm running 20 times. Since ALNS-PR does not report results from small-scale customers, it is not shown in the table. It can be seen from the data in the table that CPLEX successfully solved five examples, and DE-LNS and CO-GA also found the optimal solution, while P-SA performed slightly worse on rcdp2501. For the other four examples, CPLEX terminates early due to the 'memory overflow' condition, while the other three heuristic algorithms achieve be er solutions. The solution of the DE-LNS algorithm is be er than that of P-SA on rcdp5001 and rcdp5007. It is worth mentioning that DE-LNS performed very stably in these instances, achieving a standard deviation of 0 across 20 independent runs.
In general, the algorithms' performance on instances of larger scales is of more interest. The medium-scale instances were first solved by Co-GA, but later p-SA and ALNS-PR algorithms refreshed most of the results. Table 5 gives the results of the DE-LNS algorithm and the best results of the three comparison algorithms obtained from the literature on medium-scale instances. We evaluate the performance of DE-LNS mainly from two perspectives, namely, the best and average performance achieved in 20 independent runs. Gap1, Gap2, and Gap3 in the table represent the deviation of the optimal solution between DE-LNS and CO-GA, P-SA, and ALNS-PR, respectively.
DE-LNS obtains be er results on 41 instances compared with Co-GA and obtains approximately equal solutions to Co-GA on 15 instances, and there is no poor solution. More precisely, there are currently 16 instances that have been optimized in terms of vehicle number (the first optimization objective) in 56 instances, and the number of vehicles in the remaining 40 instances is the same as that of the Co-GA algorithm. For instances with the same number of vehicles, the DE-LNS algorithm improves the travel distance of 25 instances, and the other 15 instances obtain the same or approximately equal travel distance as CO-GA. From the data in the table, it can be found that the travel distance in some cases increases to a certain extent when the number of vehicles is optimized. This is a normal phenomenon because the number of vehicles is the first optimization goal, and its complete reduction may sacrifice optimization in terms of travel distance.
The DE-LNS algorithm achieves be er solutions in 40 instances compared with the p-SA algorithm in 56 instances, and the same or approximate results are obtained in 16 instances. Specifically, the DE-LNS algorithm optimizes the number of vehicles in 6 instances and the travel distance in 35 instances. The number of vehicles in the remaining 14 instances is the same as that of the p-SA algorithm, and only 2 instances have one more vehicle. A total of 10 solutions are the same as the results of the p-SA algorithm in 16 identical or approximate solutions, and the deviation of the remaining 6 solutions from the p-SA algorithm is less than 1%. Compared with the results obtained by the ALNS-PR algorithm, DE-LNS performs well, and it optimizes five instances. On the remaining 51 instances, the DE-LNS algorithm achieves the same or a close solution as the ALNS-PR algorithm, and there is no worse solution.
From the data in the table, it can be found that DE-LNS obtained the best solution in 21 of the 56 instances. Specifically, these instances are evenly distributed in all instances, including 6 in Rdp, 12 in Cdp, and 3 in Rcdp. These results show that DE-LNS has strong robustness and stability. It is worth mentioning that the average performance of 20 independent runs of DE-LNS is the same as its optimal performance in 37 instances out of 56 instances. The average performance of DE-LNS is still significantly be er than that of Co-GA and p-SA. The results obtained by DE-LNS still have great competitiveness compared with the ALNS-PR algorithm. Such results demonstrate that DE-LNS performs very stably, although it is a randomized algorithm in nature. In conclusion, all of the above observations prove that the DE-LNS algorithm has strong robustness and optimization ability and provides a high-quality solution for the VRPSPDTW problem.

Real Case Application
We take the production data in the actual production process of an intelligent workshop as an experimental case and then simulate the process of AGVs performing logistics transportation tasks. The purpose is to verify the benefits of the proposed green vehicle path planning model and to evaluate the practical applicability of DE-LNS. The workshop is mainly composed of automation workstations (material buffer and CNC), material handling equipment (AGV), a production scheduling center, and a cargo storage center (depot). The workshop topology map we established is shown in Figure 8. The map is asymmetric, and it consists of 27 nodes and 32 edges. The number 0 is the production scheduling center, and numbers 1-26 are the workstation. There are 10 AGVs in the workshop. The AGV specific parameters required for this experiment include the AGV's own weight and the maximum load of the body, as well as the rated power and speed. These parameters are derived to calculate the energy consumption of AGV driving. The parameter data comes from the AGV working in the actual workshop. The parameters are shown in Table 6. Each edge of the workshop presents a transport path. The topology map of the workshop is displayed as 26 points. We use the production data of the workshop to verify the energy saving effect of the model and the applicability of the algorithm. Points are represented in the following Table 7. ID is the number of the workstation and the production scheduling center, coordinate is their coordinates in the topology network, and pickups and deliveries are the weight of the goods taken and delivered respectively. We extract the specific data of a logistics transportation task and abstract the task execution time as a time window. Because of the asymmetry of the graph, the shortest distance between any two nodes i and j in the map was obtained by using the Dijkstra algorithm in advance. There is no production priority relationship between workstations, and each AGV corresponds to a specified transportation path. We analyze the calculation results of DE-LNS and current advanced heuristic algorithms under the case conditions designed in this paper. The purpose is to verify the effectiveness of DE-LNS in the context of real cases and the performance of De-LNS on GVRPSPDTW in the workshop.
The experimental conditions and parameters of the algorithm are the same. Each algorithm repeats the calculation 20 times, and then statistical analysis is performed. We take the optimal results from 20 calculations as shown in Table 8. The optimized curve is shown in Figure 9.
According to the experimental data, the calculation result of the DE-LNS algorithm is the best. The energy consumption of the transportation process obtained by the DE-LNS algorithm in the actual AGV logistics transportation task is 3.8% lower than that of GA and 2.1% lower than that of the DE algorithm. The standby energy consumption is optimized by 2.6% compared with GA and 36.2% compared with the DE algorithm. Although the results of the DE-LNS algorithm are larger than those of the VNS algorithm, LNS algorithm, and DE-VNS algorithm in terms of standby energy consumption, there are more optimizations in total energy consumption and transportation energy consumption. The total energy consumption calculated by the DE-LNS algorithm is 0.06% lower than that of VNS, 0.08% lower than that of LNS, and 9.34% lower than that of DE-VNS. Note that the VNS algorithm calculates more AGVs. The reason is that the algorithm designed in this paper uses more optimization methods. We have improved the defects of the DE algorithm, such as premature convergence and search stagnation, and made up for the defects of DE through the LNS algorithm. At the same time, we optimized and screened the population by the squirrel migration strategy. Based on the above reasons, the DE-LNS algorithm optimizes the service order of AGV to workstations in the total transportation route, and the route can control the service time of AGVs to more workstations within its corresponding time window. Therefore, the DE-LNS algorithm designed in this paper has higher solution quality in optimizing the total energy consumption of AGVs. The solution also shows good optimization ability for actual production cases. For this reason, the effectiveness and applicability of the algorithm are verified.
From the optimization curve in Figure 9, it can be seen that the total energy consumption of an AGV in performing logistics transportation tasks gradually decreases with the increase in the number of iterations. The optimal solutions of GA, DE, and DE-VNS are obtained in 434, 336 and 202 generations, respectively. The convergence speed of DE-LNS is 58.1%, 45.8%, and 9.9% higher than that of GA, DE, and DE-VNS, respectively. The LNS algorithm obtains the optimal solution in 181 generations, and the convergence speed is almost the same as that of DE-LNS. Although the VNS algorithm obtains the optimal solution in 174 generations, the calculation result is worse because the VNS algorithm falls into the local optimal solution. In the experimental environment of this actual case, the energy consumption value and convergence speed calculated by the algorithm are used as evaluation indexes to evaluate the effectiveness and applicability of the algorithm. The reason is that DE-LNS designed in this paper designs a variety of strategies to improve the convergence speed of the algorithm. We design an adaptive scaling factor based on the logistic function, which gradually decreases during the iteration process to improve the convergence speed and local optimization ability. The squirrel migration strategy was used to optimize the population on a large scale and screen out some elite individuals. The powerful local search ability of LNS is used to destroy and repair these individuals,and finally, a series of new sets are obtained. This process can help the algorithm quickly find the optimal solution set and accelerate the convergence of the algorithm.

Energy Consumption Cost Comparison of the GVRPSPDTW and the VRPSPDTW
In this section, we study the impact of the green vehicle path planning model designed in this paper on energy consumption costs. Therefore, we compare VRPSPDTW and GVRSPDTW in a real case. The energy consumption value and driving distance calculated by the same algorithm are used as the evaluation indexes of the model. When the objective function (8) is replaced by the minimum total travel distance, GVRPSPDTW is reduced to VRPSPDTW. As shown in Formula (33) In the formula, i CT is the unit loss cost when an AGV arrives earlier than the time window, and i CL is the unit time penalty cost when an AGV arrives later than the time window. The constraints of the two models are the same. We hope to reduce how much energy consumption can be reduced by observing the increase of AGV driving distance. In order to achieve this goal, we calculate the optimal total travel distance and the optimal the total energy consumption obtained by GVRPSPDTW and VRPSPDTW respectively. We design 5 experimental cases in the logistics transportation tasks of different time periods in the workshop and use the DE-LNS algorithm to perform 20 experiments on the case to obtain the best solution. The results are shown in Tables 9 and 10.
According to the experimental data, although the total travel distance calculated by the green vehicle path planning model increased by 10.03% on average, the total energy consumption decreased by 8.01% on average. The reason is that the green vehicle path planning model not only considers the influence of the AGV driving path on energy consumption, but also deeply analyzes the mechanical force affecting energy consumption during AGV driving. We add the consideration of real-time changing load and AGV power. The results show that the green vehicle path planning model designed in this paper is very effective in terms of energy consumption. In terms of total transportation time, the green mathematical model is 15.3% lower than the VRPSPDTW model. This result shows that the green vehicle path planning model not only achieves the purpose of reducing energy consumption but also improves the overall efficiency of AGV logistics transportation tasks in the workshop. For this reason, the green vehicle path planning model can have great application value in an actual intelligent workshop.

Conclusions and Future Research Work
In this paper, the GVRPSPDTW in the intelligent workshop is studied, and the conclusions are as follows.
(1) We build a green AGV path planning model. The combined optimization goal is to minimize the energy consumption of the AGV transportation process and standby energy consumption. The model can objectively and comprehensively reflect the total energy consumption of AGVs in an intelligent manufacturing workshop. (2) We design a differential evolution algorithm with large neighborhood search to solve the problem studied in this paper. The algorithm can improve the quality of the initial solution by improving the neighborhood searchability. (3) DE-LNS adopts adaptive crossover and mutation strategies in evolutionary operations. We design an adaptive scaling factor based on a logistic function and adaptive crossover operator to ensure effective convergence and robustness. We verify the effectiveness of the algorithm through the VRPSPDTW international standard dataset. (4) The experiments in this study show that the designed algorithm still has strong applicability to the AGV routing problem in an actual intelligent workshop. The solution quality of DE-LNS is be er than that of advanced heuristic algorithms. The green mathematical model we established is very effective for the energy consumption of AGV's logistics transportation tasks. There is no doubt that AGVs, as the main transportation equipment for intelligent workshops, have great energy-saving potential.
Therefore, this paper provides a potential basis for future research. At the same time, in order to be er solve the problem, this paper does not consider the speed change caused by AGV turning. To further realize the development goal of energy saving and emission reduction in the workshop, we will consider the impact of AGV speed changes on energy consumption in future work. Of course, whether different types of AGV can save more energy is also one of the focuses of our research. We will conduct in-depth research on the integrated scheduling of processing resources and transportation resources from the perspective of energy saving. The GVRPSPDTW problem of multi-process processing in intelligent workshops is also one of our important research contents. In the future, the solution algorithm will be further optimized.
Author Contributions: Conceptualization, X.Z. and F.G.; methodology, X.Z. and F.G.; investigation, X.Z. and F.G.; writing-original draft preparation, X.Z. and F.G.; writing-review and editing, F.G.; supervision, F.G. and X.T.; and funding acquisition, X.Z. and X.T.; All authors have read and agreed to the published version of the manuscript.