Distribution Route Optimization for Electric Vehicles in Urban Cold Chain Logistics for Fresh Products under Time-Varying Traffic Conditions

Electric vehicles (EVs) have been widely used in urban cold chain logistic distribution and transportation of fresh products. In this paper, an electric vehicle routing problem (EVRP) model under time-varying traﬃc conditions is designed for planning the itinerary for fresh products in the urban cold chain. The object of the EVRP model is to minimize the total cost of logistic distribution that includes economic cost and fresh value loss cost. To reﬂect the real situation, the EVRP model considers several inﬂuencing factors, including time-varying road network traﬃc, road type, client’s time-window requirement, freshness of fresh products, and en route queuing for charging. Furthermore, to address the EVRP, an improved adaptive ant colony algorithm is designed. Simulation test results show that the proposed method can allow EVs to eﬀectively avoid traﬃc congestion during the distribution process, reduce the total distribution cost, and improve the performance of the cold chain logistic distribution process for fresh products.


Introduction
With the increase of greenhouse effect, clean energy is the energy of the future for traffic development. With government support, electric vehicles (EVs) are becoming an industry trend as people seek ways to reduce emissions to protect the environment. EVs, an ideal means of transportation to replace conventional fuel vehicles, are energy-efficient, green, and environmentally friendly. erefore, EVs have received extensive attention to combat the deteriorating air pollution and global climate change. Many countries have set policies and developed several projects to improve the penetration of EVs in their daily lives. By 2017, EVs (including pure electric vehicles and plug-in hybrid vehicles) exceeded 3 million around the world, an increase of 57% compared with 2016. Data show that the number of pure electric vehicles in China rose to 3.1 million in 2019, which is the highest number of EVs compared with the EV numbers for a few years ago, accounting for about 40% of the global EV ownership. It is foreseen that the number of EVs will break through 200 million in 2030 [1]. For instance, Beijing released the "Beijing New Energy Logistics and Distribution Vehicle Priority Implementation Plan" in 2019: by the end of 2020, light electric logistic vehicles under 4.5 tons will account for 90 percent of the total freight vehicles. While they have advantages (e.g., low energy consumption and zero pollution), EVs also have disadvantages, namely, long charging time, limited battery capacity, and relatively low mileage on a single charge [2,3]. With the increasing popularization of EVs, how to scientifically dispatch EVs and optimize logistic distribution routes under the constraints of EV technology to achieve economic and environmental protection goals has continually been a focus of researchers. e EVRP is an extension of vehicles routing problem (VRP), but there are obvious differences between them. (1) Different transport means result in different energy consumption functions. Compared with traditional vehicles, EVs generally have a range constraint of about 100-200 km. Due to the characteristics of urban time-varying networks, road congestion, and customer point service requirements, the EV range may be reduced. To meet distribution needs, EVs may need to be recharged in transit. (2) e EVRP needs to consider the charging strategy of EVs in the routing planning process [4]. (3) If EVs need en route fast charging in the EVRP, the selection of social charging stations, queuing, charging time, and charging efficiency have an important impact on path planning. erefore, the EVRP is more complex compared with traditional VRP.
Currently, the EVRP is a focal area of research. Beltran et al. [5] proposed using EVs in urban transportation and established an optimization model with the minimization of the transportation distance as an objective. Erdoan and Miller-Hooks [6] were the first to put forward the EVRP. In view of the problems of EVs (i.e., low mileage on a single charge and a lack of charging facilities), they studied an EVRP with time windows that accounted for charging strategies. In addition, they established a mixed integer model with the minimization of the travel distance as an objective and solved it using an improved saving algorithm and a clustering algorithm. Schneider et al. [7] studied the last-kilometer logistics of EVs and took into account the effects of client's time window, vehicle load and single-charge mileage constraints, and en route charging. ey solved the problem using a variable neighborhood search and tabu search-combined heuristic algorithm. Goeke and Schneider [8] studied an EV power consumption model based on relevant factors, such as vehicle speed, vehicle load, and road gradient, and established an EVRP model for driving and energy consumption costs by considering charging strategies. Yang et al. [9] considered the time-of-use (ToU) electricity price when optimizing and solving the EVRP and designed an EV route planning and charging strategy selection model under the ToU electricity price. ey also designed a genetic algorithm for solving this model. Desaulniers et al. [10] took into consideration four charging strategies when optimizing an EVRP with time windows. ey obtained an optimum charging method by optimizing the four charging strategies under a client timewindow constraint. Jie et al. [11] and Hiermann et al. [12] examined a combination of multiple vehicle types and the logistic distribution route problem of EVs. Jie et al. [11] distinguished various vehicle types based on factors such as battery capacity, power consumption rate per unit, and vehicle load and designed a branch-and-price algorithm for finding an optimal solution to the problem. Hiermann et al. [12] rationally configured the fleet structure for distribution based on the differences between EV types in terms of energy consumption and load. To better conform to reality, Li et al. [13] planned mixed fleets of conventional vehicles and EVs for enterprises and thereby provided decision-making support to the government for choosing locations of charging facilities and to enterprises for logistic operation management.
Cold chain logistics for fresh food products is also a focus of research in China and elsewhere. Tiwari and Chang [14] optimized the route of vehicles in cold chain logistics for fresh products with the minimization of the travel distance and carbon emissions as objectives. Kumar et al. [15] established a multiobjective (total cost and total emissions) cold chain transportation route problem. Rabbani et al. [16], Amorim and Almada-Lobo [17], and Wu et al. [18] introduced the freshness of fresh products into cold chain distribution route optimization for fresh products. Zhang and Liang [19] established a cold chain vehicle route optimization model with the minimization of the total cost of cold chain logistic network construction and transportation as an objective. Devapriya et al. [20] established a coordinated production and distribution optimization model for perishable products with the minimization of the total distribution cost as an objective. ey also designed a mixed heuristic algorithm to solve the model. Li et al. [21] investigated a new integrated planning problem for intelligent food logistic systems. ey formulated a biobjective mixed integer linear programming model that considers the two objectives: minimizing total production, inventory, and transportation cost and maximizing average food quality. e VRP is a typical NP-(nondeterministic polynomial-) hard problem. As the number of nodes increases rapidly, an exact algorithm is rather difficult to employ to solve largescale VRP. Currently, there are many VRP optimization algorithms, such as genetic algorithm [22], local search algorithm [23], tabu search algorithm (TSA) [24], simulated annealing search algorithm [25], and ant colony algorithm (ACA). According to the characteristics of the EVRP model in this paper, the ant colony algorithm is designed to solve the problem. Ant colony algorithm is a bionic random search algorithm originated from nature. Many scholars at home and abroad have applied it to various fields of life. It has made more outstanding achievements in combinatorial optimization, such as traveling salesman problem [26], scheduling problem [27], robot control system problem [28], and personalized trip planning problem [29]. ACA has the characteristics of distributed computing, no central control, and indirect communication between individuals, so it is easy to combine with other algorithms. For example, ACA combines with BP neural networks [30], genetic algorithm [31], etc. VRP is a complex discrete-time system, and many literature studies show that ACA has a broad development prospect in the field of VRP. Chen et al. [32] studied the optimization of the transportation routing problem for fresh food by an improved ant colony algorithm based on tabu search. Tang et al. [33] studied the low-carbon vehicle routing problem model; they designed an improved ant colony algorithm for solving this model, and chaotic disturbance mechanism was introduced to update the ant pheromones on the path. Fang and Ai [34] studied the hybrid ant colony algorithm to solve the cold chain logistic distribution optimization model.
A review of the literature published by researchers in China and elsewhere finds that in-depth research can still be conducted in the following several areas. (1) e available studies have mostly independently investigated the EVRP and the cold chain distribution route optimization problem but less frequently considered the optimization of urban cold chain distribution with EVs. (2) Studies on the effects of the external environment on the distribution with EVs mostly have considered that the vehicle speed is constant and there is only one road network type, overlooking the impact of dynamic traffic in a road network on the actual distribution process with EVs. (3) e available studies have mostly considered the biggest continue voyage course, but in the process of electric vehicle distribution, when and where to quickly charge is also a key issue. (4) When examining en route EV charging, most studies have only considered the charging time but overlooked the queuing time at the charging station. All of the above are the primary motivations of this research. e main contributions of this paper can be summarized as follows: (1) According to the characteristics of fast charging of EVs in social charging stations, the queuing and charging process model of EVs is employed. (2) An improved adaptive ant colony algorithm is designed to solve the distribution optimization model with time windows for EVs in cold chain logistics for fresh products. (3) e charging station placement strategy based on the minimum power consumption principle is designed to solve the key problem of when and where to fast charge in the distribution process. (4) Based on the perspective of considering the economic cost of distribution and the cost of fresh value loss, the EVRP for urban fresh distribution under the time-varying road network is studied. e remainder of this paper is organized as follows. In Section 2, we describe the mathematical model and the cost analysis of EVs in cold chain logistics. In Section 3, the improved ant colony algorithm is proposed with the entire design process of the algorithm. In Section 4, experiments with the Solomon benchmark are performed to analyze the impact of different types of cases and optimization objectives on the performance of the algorithm and compared with other excellent algorithms. Finally, Section 5 gives an overall conclusion.

Description of the Problem and Basic Assumptions.
e problem of distribution route optimization for EVs in cold chain logistics for fresh products under time-varying urban traffic conditions can be described as follows. A certain distribution center distributes fresh agricultural products to end clients via EVs. e EVs leave the distribution center, complete all the distribution tasks, and then return to the distribution center. e location, time window, and demand at each client point are known. e time-varying urban road network traffic information can be obtained from traffic control authorities. During the distribution process, the EVs can be charged at public charging stations as needed. Under these conditions, a distribution plan is formulated with the minimization of the total distribution cost as an objective. To clarify the scope of application of the study, the following assumptions are made: (1) ere is only one distribution center in the distribution network. is center has an ample supply of EVs of the same type. (2) Each EV leaves the distribution center on a full charge and immediately returns after completing the tasks. (3) e demand at each client point is less than the capacity of each EV.
ere is also a service time-window requirement. An EV can reach a client point to provide service at an earlier or later time but then will pay a certain penalty. (4) Vehicle speeds in all the sections of the road network are updated every 10 min. Considering the time-varying nature of network traffic, road sections are classified based on the areas where two arbitrary nodes are located. e road sections differ in time-varying characteristics. (5) When its remaining battery capacity cannot meet the distribution need, an EV will go to a charging station for a fast charge. (6) Each EV has a fixed maximum load capacity. In addition, one and only one EV serves each client point. (7) e locations of the public charging stations are known. e public charging stations have the same relevant charging parameters. e queuing time for charging is uncertain and depends on the time of day and the area.

Symbolic
Parameters. N is the set of customer points; 0 represents the distribution center; F is the set of charging stations; V represents the set of all nodes in a logistics network, V � N ∪ F ∪ 0 { }. T is the set of the time period throughout the day, T � T 1 , T 2 , . . . , T P ; P is the total number of time periods. R is the set of the road section in the road network, R � R 1 , R 2 , . . . , R H ; H is the number of road section types. A is the set of area, A � A 1 , A 2 , . . . , A G ; G is the number of the area types. Other descriptions about symbols are explained in Table 1.

Analysis of Charging Needs and Power Consumption
Performance. Currently, there are three modes of charging an EV, namely, battery replacement, slow charging, and fast charging. Considering the high timeliness requirements of urban cold chain distribution of fresh products, fast charging in the delivery process is considered in this study. Time requirements are not considered during the slow-charging process. To overcome its low single-charge mileage problem during the distribution process and improve its distribution efficiency, an EV needs to undergo a fast charge at a public charging station. In addition, the relatively long queuing and charging process is an important component of distribution with EVs. See elsewhere [35] for a discussion of the queuing process. e queuing and charging process of an EV at a public charging station conforms to the M/G/s standard multiserver queuing system model (M means that the time interval for an EV to reach a charging station follows the exponential distribution; G means that the charging time of an EV follows the normal distribution; there are a total of s servers). Let ρ i be the utilization ratio of charging equipment in the charging station i (ρ i � (λ i /s i μ i )). According to the classical little equation, the average queuing time for the k th EV at charging station i is During the fast-charging process, the charging time for the k th EV at charging station i is as follows: e energy consumption of electric vehicles is not only related to the nature of the vehicle but also related to the load and speed. In addition, part of the output power of the electric refrigerated vehicle is converted into mechanical power under certain transmission efficiency, and the parameters are shown in Table 2. When an EV whose load is Q v Mathematical Problems in Engineering runs on a flat road at a speed of v, the running power P(Q V , v) is expressed as follows [36]:  [37,38], the real-time speed of a vehicle within a sufficiently short time is used as its fixed speed within the time period. To reflect the variation in road sections where vehicles travel, it is also necessary to take into consideration the types of road sections. On this basis, the vehicle speed in the current road section is determined at the current moment. e following shows a step function for the speed at various times on an arbitrary road section: where the velocity value does not change in the same period of time in the road (i, j). A number of factors, including road e fixed cost of per EV (yuan) P 2 e labor and rent cost of EV per unit time (yuan/min) P 3 e unit electricity consumption (yuan/kWh) P 4 e unit electricity price of fast charging (yuan/kWh) P 5 e refrigeration costs which generate during transportation process of unit time (yuan/h) P 6 e refrigeration costs which generate during unloading process of unit time (yuan/h) P 7 e price of the unit commodity (yuan/kg) P 8 e cost of waiting for the unit time if refrigerated EV arrives at customer node in advance (yuan/min) P 9 e cost of punishing for the unit time if refrigerated EV is late to the customer node (yuan/min) λ i e number of EVs that reach the charging station i (vehicles/h) μ i e number of EVs that have completed charging per unit time in the charging station i (vehicles/h) s i e number of charging piles in the charging station i z 1 e freshness decay coefficient during the transportation process z 2 e freshness decay coefficient during the unloading process t Q ik e average queuing time for EV k th at charging station i (min) ω A i e area factor for charging demand which is introduced to differentiate the charging demand between different areas, A i ∈ A t C ik e charging time for k th EV at charging station i (min) r c e charging speed (kWh/min) of the public charging station CN max e maximum number of charging cycles CE e cost of a battery (yuan) x ijk e binary variable that is 1 when the k th EV is driving on road (i, j) and 0 otherwise y ik e binary variable that is 1 when the k th EV is servicing customer point i and 0 otherwise Z ik e binary variable that is 1 when the k th EV is charging at point i and 0 otherwise type, travel time period, and time-varying vehicle speed, are analyzed. e whole-day distribution time is evenly divided into multiple time periods. In a time-varying network, vehicle travel speeds are updated once every time period. Here, the traveling of a vehicle on a specific road section is used as an example. e travel time calculation process for the vehicle in the time-varying network is analyzed as follows: (1) Determining the road type r based on the areas where the points i and j are located, and, on this basis, the vehicle speed during the corresponding time period is determined. Go to (2). (2) e distance traveled by the EV when it travels at the corresponding speed in the remaining time that is in the current period, plus the distance that the EV has already traveled, is calculated, and whether the total distance obtained exceeds the distance of road (i, j) is determined. (3) If the total distance exceeds the distance of road (i, j), the travel time is the sum of the travel time experienced by the EV before entering the current time period and in the current time period. Where the driving time of the EV in the current period is obtained as follows: the distance of road (i, j) minus the distance previously experienced divided by the driving speed corresponding to the current period. (4) If the total distance does not exceed the distance of road (i, j), the EV will continue to drive in the next period. Go to (2) and repeat.
e pseudocode for calculating the traveling time t ijk in the time-varying traffic network is detailed in Algorithm 1, and the corresponding parameters and descriptions are shown in Table 3.

Cost Analysis of EVS in Cold Chain Logistics.
Based on the actual overall profit of fresh-food cold chain distribution enterprises, the minimization of the total distribution cost (including economic cost of distribution and fresh value loss cost, where the economic cost of distribution consists of the fixed cost, time cost, energy consumption cost, charging cost, refrigeration cost, and penalty cost) is set as the optimization objective. On this basis, a route optimization model for EVs in cold chain logistic distribution of fresh products is established.
(1). Fixed Cost of EVS. e fixed cost of EVs (C 1 ) is primarily their front-end cost. is cost is only related to the number of refrigerated EVs that are used: (2). Travel Cost of EVS.
is cost (C 2 ) is primarily the product of the cost per unit time and the total running time.
e running time mainly includes the road travel time and the client-point service time as well as the queuing time and charging time at public charging stations. us, the time cost of EVs can be calculated using the following equation: (3). Energy Consumption Cost. Energy consumption cost refers to the energy consumption cost of the vehicle in the whole process of driving. Based on the calculation of power consumption above, the total energy consumption cost (C3) is (4). Charging Cost. e charging cost (C 4 ) is the total cost of fast charging at public charging stations during the distribution process. us, the fast charging cost of EVs is (5). Refrigeration Cost. e refrigeration cost of a refrigerated vehicle (C 5 ) is generally the cost of refrigerating agents and energy consumption to maintain the temperature within the vehicle. According to Fang et al. [34], the refrigeration cost consists of the refrigeration cost during the transportation process and the refrigeration cost during the loading and unloading of goods:  Mathematical Problems in Engineering 5 (6). Penalty Cost. is cost (C 6 ) is incurred by the vehicle as a result of violating the client's time window during the distribution process. Client i accepts service outside the time window, but there is a penalty cost for failing to provide service within the appointed time, which is equal to (7). Battery Loss Cost. Fast charging will significantly reduce the service life of batteries. When the number of times a battery has fast charged exceeds the maximum number of times that it can withstand fast charging, the battery will be discarded. us, in this study, the battery loss cost (C 7 ) is represented by a function positively correlated with the number of fast charges an EV has undergone: (8). Fresh Value Loss Cost. ere are two primary causes of damaged goods, namely, the passage of time during the distribution and transportation process and changes in the surrounding environment during the loading/unloading and service processes. According to Fang et al. [34] and Li and Fan [39], a fresh-product freshness decay function is introduced to calculate the cost of damaged goods (C 8 ): Equation (12) shows the fresh value loss cost during the transportation process. Equation (13) shows the fresh value loss cost when loading goods onto the vehicle and unloading goods from the vehicle.

2.3.4.
Modeling. Based on the above description, a distribution route optimization model for EVs in cold chain logistics for fresh products in a time-varying road network is established, as shown below, by considering the difference in en route charging time and queuing time for charging and with the minimization of the total cost of logistic distribution that includes economic cost (consisting of the fixed EV cost, time cost, energy consumption cost, refrigeration cost, charging cost, time penalty cost, and battery loss cost) and fresh value loss cost as the objective function Z: Constraint (16) indicates that each customer point can only be serviced by one EV and all customer points are serviced; constraint (17) indicates that each vehicle's starting point and returning point are the distribution center; constraint (18) indicates that the delivery vehicle's cargo weight cannot exceed the vehicle weight; constraint (19) indicates that the vehicle is fully charged when departing from the distribution center; constraint (20) indicates that the  (21) indicates the vehicle's power level limits when arriving at any node, which cannot be lower than the minimal power level for normal driving or higher than the battery capacity constraint; constraint (22) indicates the power levels when the vehicle leaves node i and arrives at node j, and M0 is a very large positive constant number; constraints (23) to (24) indicate the time constraints: constraint (23) indicates the time constraints on a vehicle arriving at and departing from node i, constraint (24) indicates the time relationship between the electric vehicle leaving node i and arriving node j, and M0 is a very large positive constant number; constraints (25) to (27) indicate variable value constraints.

Improved Ant Colony Algorithm (IACA) Design
Because the ant colony algorithm has a positive feedback mechanism, strong solving ability and robustness, and can be integrated with other optimization algorithms, in this study, we modified the basic ant colony algorithm to propose an improved ant colony algorithm for solving the above problem. e corresponding parameters and descriptions of IACA are shown in Table 4.

State Transition Strategy.
e ants determine the subsequent service client point j through the state transition strategy based on client point i. e following goes over the procedure in detail.
First, the transition probability P m ij for ant m to travel from the node i to the node j is calculated: where α, β, and θ are the weights of the pheromone concentration, the visibility, and the time-window span, respectively. According to the method proposed by Li et al. [40] for obtaining the adaptive and expected pheromone heuristic factors, en, to prevent the algorithm from falling into a local optimum and avoid premature stagnation, the state transition strategy is improved based on the Roulette wheel method using the following equation: randomly selected point j, j ∈ C T m , based on Roulette wheel method, ∀R 1 > R 0 .
An algorithm parameter R 0 in the interval of [0, 1] is introduced. When ant m selects the subsequent client point, a random number in the interval [0, 1] will be generated. e next client point to which the ant transitions is determined by comparing R 0 and R 1 . When R 1 ≤ R 0 , the point with the highest transition probability is selected as the next service client point; otherwise, the ant selects the next client point according to the procedure of the Roulette wheel method [34].

Pheromone Update Strategy.
To allow the search process to be more instructive, the pheromone of the road sections composed of the routes established by all the ants upon completion of traversing is updated. e update rules are as follows: Δτ m ij , where τ old ij is the pheromone of the road section (i, j) before the current iteration and τ new ij is the updated pheromone after the current iteration.

Charging Station Placement Strategy Based on the Minimum Power Consumption Principle.
During the distribution process, an EV needs to satisfy not only the vehicle load constraint but also the power constraint. at is, the en route charging behavior of an EV needs to be taken into consideration. erefore, when constructing a feasible solution for a single ant, it is necessary to first construct a feasible solution that satisfies the vehicle load constraint and then determine the power constraint for the EV at each client point. If the solution does not satisfy the vehicle load constraint, a charging station is placed based on the minimum power consumption principle. e following shows the procedure in detail: (1) First, an initial feasible solution that satisfies the vehicle load constraint is generated. en, the remaining power E A ik is calculated when the ants reach each client point. If E A ik > E low at all the client points, then the EV does not need to be charged en route. e operation is terminated. Otherwise, the first client point i is determined where the power constraint is not satisfied and go to operation (2). e EV goes on to fulfill the service at client point i. If the constraint is unsatisfied, another client point is found by forward recursion and operation (2) is repeated until finding the appropriate charging station, which is then placed. e corresponding process is shown in Figure 1.

Algorithm Process.
By the integration of all the steps, described above, we propose an improved ant colony algorithm for solving the distribution route optimization model for EVs in urban cold chain logistics for fresh products under time-varying traffic conditions. e optimization objective of the model is the total distribution cost. e general structure of the improved ant colony algorithm is shown in Algorithm 2, the corresponding process is shown in Figure 2, and the relevant parameters are shown in Table 4. In order to further analyze the relationship between cost elements in the model, we take economic cost and fresh value loss cost as two targets, respectively, to solve the problem. e biobject optimization question is dealt with by the main object method [41]. e economic cost is chosen as the main objective, and the fresh value loss is taken as the constraint condition to solve the biobject optimization problem. e selection of the main objective function depends on the strategic planning of fresh e-commerce enterprises. e specific transformation is as follows: e objective function (32) minimizes the economic cost; the objective function minimizing the fresh value loss cost is converted to constraint (33) by the transformation. Constraint (33) which indicates the ratio of the fresh value loss to the total fresh value of each customer point should be lower than the minimum threshold ε. According to the characteristics of the model, we further improve the ant colony algorithm. For each ant route selection, we added the fresh value loss ratio constraint. e specific process is shown in Algorithm 3.

Test Settings.
Because of the diversity of customer distribution in fresh cold chain distribution, the data of the test example were from the international standard EV route databases [42]. Considering various factors of the EV distribution, we used the data of the type C (concentrated distribution), type R (random distribution), and type RC (mixed distribution) examples of the international standard of the EVRPTW database. In each example, there are 20 charging stations (No.1 to No.20), 100 customer points (No.21 to No.120), and 1 distribution center (No.0). Considering that there are multiple scenarios for the coordinates of fresh-product cold chain distribution clients, (1) the starting and ending times of the model are set to 7:00 (corresponding to time 0, i.e., 0 min) and 24:00 (corresponding to 1,020 min), respectively. e length of each time period is 10 min. According to urban traffic patterns, the period 7:30-9:00 (i.e., 4 th -12 th periods) and the period 17: 30-19:00 (i.e., 64 th -72 nd periods) are set as traffic congestion periods (vehicle speeds vary with the road-section type), and the other periods are set as normal periods. According to Xiao and Konak [37], let ϖ � 0.1. ree vehicle speeds [VR(1 − ϖ), VR(1 + 2ϖ), VR(1 − 3ϖ)] are allowed. Based on the T value of the time period, the remainder function  while C T m ≠ ∅ (7) Select customer point j according to state transition strategy (8) if L T m + q j ≤ Q (9) if E A m (j) < E low (10) Insert the charging station in the route according to the minimum power consumption principle and update the R T m information (11) end if (12) Insert the customer point j to the end of the route R T m , update tabu list, delete j from C T m . If j is the first service node, calculate the departure time. (13) else L T m + q j > Q (14) Ant m returns to the distribution center, insert 0 to the end of R T m , VN m � VN m + 1, L T m � 0 (15) end if (16) end while (17) end for (18) Calculate the fitness of each ant and select the optimal value for this iteration f NC B (19) Update the segment pheromones according to each ant's path and fitness (20)  (1) For m � 1: Select customer point j according to state transition strategy (4) if Insert the charging station in the path according to the minimum power consumption principle and update the R T m information (7) end if (8) if Insert the customer point j to the end of the path R T m and delete j from C T m (10) if j is the first customer point that is served by EV (11) Calculate the EV departure time (12) end if (13) else (F A m (j)/P 7 · d j ) > ε (14) Ant m returns to the distribution center, insert 0 to the end of R T m , VN m � VN m + 1, L T m � 0 (15) end if (16) else L T m + q j > Q π � mod(T, 3) is used. e π values of 0, 1, and 2 correspond to three vehicle speeds during the normal time periods. e vehicle speeds vary with the road type. (2) In the actual distribution, the nodes are not connected by straight roads but by multiple roads; therefore, using the straight line distance instead of the distance between two points will lead to quite different from the actual situation. is paper introduces the circuitous coefficient δ(δ � 1.5) [43], and d ij is the product of the straight-line distance between two points and the circuitous coefficient. (3) To be more consistent with real life, the areas where any nodes are located within the road network are classified into three types with urban inner-ring highways as boundaries ((40, 40) is the center of each inner ring; the radii of the second and third rings of the city are 20 and 30, respectively). e area within the second ring is denoted by A 1 . e area between the second and third rings is denoted by A 2 . e area outside the third ring is denoted by A 3 . (4) Based on the area where any two arbitrary nodes are located, road sections are classified into six types. Table 5 summarizes the attributes of the road sections. e vehicle had an unloaded weight of 4000 kg, with a capacity of 200 weight units; assuming that one weight unit is 10 kg, the vehicle's load capacity is 2 tons. e model parameters are shown in Table 6, and the algorithm parameters are shown in Table 7.
e calculation was performed using MATLAB R2016 on a PC with a 2.0 GHz CPU and 8G memory.

Analysis of the Stability of the Algorithm, Charging Planning, and Departure Time Based on a Large-Scale Dataset.
To examine its stability, the proposed algorithm was tested on the R204 dataset. e program was run 10 times. e average running time was 378.491 s. is demonstrates that the proposed algorithm can produce an optimum EV route planning scheme within a relatively short time. Figure 3 shows the optimum EV route planning based on the R204 dataset (the cyan and brown rings are the second and third rings of the city, respectively). e R204 dataset was used in the example, and the detailed route optimization scheme is shown in Table 8. In Table 8, "VN" represents the vehicle number; "EVR" represents the vehicle route; 0 stands for the distribution center, 1-20 for the charging stations, and 21-120 for the customer points; "ST" represents the EV departure time/(min); "ET" represents the end of the EV distribution time; "CN" represents the number of times EV is fast charging; "CT" represents the EV charging time/(min); "QT" represents the EV queuing time for charging/(min); "TT" represents the total EV travel time/(min). e test calculation results show the following. e total distribution cost is 5,579.2 yuan. Eight EVs are used. e travel distance is 1,164.6 km. e fixed, travel, refrigeration, fresh value loss, energy consumption, charging, penalty, and battery loss costs are 1,600, 1,429.3, 397.4, 741.3, 238.2, 187.2, 760.8, and 225 yuan, respectively. Table 8 summarizes the routes and times in detail. e following can be derived from Table 8. (1) A comprehensive analysis of the EVR and the CN finds that there is a significant difference in the distribution route between EVs due to limiting factors (e.g., time-varying road network conditions, load, and battery capacity). e 5 th and 7 th EVs have the most distribution tasks, serving 15 clients. e 8 th EV has the fewest, only 9 clients. e 5 th , 6 th , 7 th , and 8 th EVs need to be charged at public charging stations during the distribution process so that they can continue to complete their distribution tasks, whereas the other EVs do not need to be charged during the distribution process.
is suggests that to complete their distribution tasks, EVs need to go to public charging stations for charging to increase their single-charge mileage and expand their distribution area within a short time. (2) An analysis of the ST finds that the 1 st , 2 nd , 3 rd , 4 th , 6 th , 7 th , and 8 th EVs depart at time 0 from the distribution center to distribute goods, while the 5 th EV departs at the other time. If departing at time 0, an EV needs to bear a relatively high time-window penalty cost. is suggests that a logistic enterprise needs to scientifically plan routes based on the actual conditions (e.g., road network conditions and time windows at client points). (3) An analysis of the QT and CT finds that the 5 th , 6 th 7 th , and 8 th EVs extend their passage range by fast charging, but the travel time of EV is increased due to charging time and queuing time. It is obvious that in all charged EVs, the highest proportion of queuing time to the travel time is 6.6% (5 th EV), and the smallest proportion is 4.7% (7 th EV). e highest proportion of charging time to the travel time is 7.9% (6 th EV), and the smallest proportion is 6.3% (7 th EV). In addition, the highest proportion of the sum of charging time and queuing time to the travel time is 13.9%. erefore, in the distribution process, charging time and queuing time cannot be ignored. (4) An analysis of the ST and ET finds that the 1 st , 2 nd , 3 rd , 4 th , 6 th , 7 th , and 8 th EVs operate during the morning peak congestion time period only, and the 5 th EV operates during the evening peak congestion time period only.
is demonstrates that the proposed algorithm can relatively effectively avoid congestion time periods and improve the distribution efficiency of vehicles.

Analysis of Vehicle Path Planning for Various Types of
Cases. Fresh-product electronic commerce enterprises have various types of clients distributed in various patterns. In view of this, solutions are found based on various Solomon datasets. In the type-C dataset, the clients are distributed in a concentrated area, the time windows are less than 1 h, and the client service time is relatively long (90 min). In the type-R dataset, the clients are distributed randomly, the time windows are greater than 1 h, and the service time is relatively short (10 min). In the type-RC dataset, the clients are distributed in a mixed manner, most of the time windows are 30 min, and the service time is 10 min. Figure 4 shows the optimal route map based on various datasets. e experimental results are shown in Table 9. In Table 9, "EX" is the dataset type, "TC" is the total distribution cost (yuan), "VD" is the total travel distance (km), "ECC" is the energy consumption and charging cost (yuan), "TMC" is the time management cost (yuan) (including the time cost and the penalty cost), "RC" is the refrigeration cost (yuan), "VL" is the cost of fresh value loss (yuan), "VN" is the number of vehicles, and "CPUT" is the running time (s). Table 9 shows the following. (1) e longest running time is only 395.1 s.
is suggests that the proposed algorithm can produce an effective optimization result for each dataset within a short time.
(2) e TC, TMC, RC, VL, and VN for the type-C dataset are the highest. is is primarily because the service time for all the type-C clients is 90 min (compared to the service time of 10 min for the R-and RC-type clients), and the total distribution time is relatively long. Consequently, the TMC, VL, and RC, which are related to the total   (3) e VD and ECC for the type-C dataset are lower than those for the type-R and -RC datasets. is is mainly because in the type-C dataset, the clients are distributed in a relatively concentrated area, so the travel distance and time are relatively short during the distribution process, resulting in relatively low power consumption. (4) In these datasets, the ECC accounts for a very small proportion of the TC, with the highest proportion being 7.8% (in the R203-type dataset). In comparison, the TMC, RC, and VL, which are related to the delivery time, each account for a relatively large proportion of the TC. e highest proportion of the sum of these three costs to the TC is 80.0% (in the C103-type dataset), and the smallest is 60.0% (in the R204-type dataset). is suggests that reducing transportation time is the key to effectively reducing the cold chain transportation cost of EVs. erefore, how to improve the single-charge mileage of EVs and reduce their charging time and frequency is the key to reducing the distribution cost of EVs in cold chain logistics for fresh products.

Analysis of Vehicle Route Planning for Various Optimization Objectives.
e route planning produced by the proposed algorithm is compared with that produced with the minimization of the total travel distance of vehicles as an objective. Table 10 summarizes the test results. In Table 10, "min TC" is the minimization of the total distribution cost (yuan), "min VD" is the minimization of the total travel distance (km), and "min PC" is the minimization of energy consumption (kWh). e meanings of "EX," "TC," "VD," and "ECC" are the same as in Table 9. Table 10 shows the following. (1) e proposed model can produce a lower total cost than models with other objectives. (2) In the type-C dataset, the TC values obtained with the three objectives differ relatively insignificantly. e TC obtained with the minimization of the VD as the objective is 4.2% higher than that produced by the proposed algorithm. e TC obtained with the minimization of the PC as the objective is only 3.4% higher than that produced by the proposed algorithm. In comparison, in the type-R and type-RC datasets, the TC values obtained with the objectives differ   "VN" represents the vehicle number; "EVR" represents the vehicle route; 0 stands for the distribution center, 1-20 for the charging stations, and 21-120 for the customer points; "ST" represents the EV departure time; "ET" represents the end of the EV distribution time; "CN" represents the number of times EV is fast charging; "CT" represents the EV charging time; "QT" represents the EV queuing time for charging; "TT" represents the total EV travel time.    "EX" is the dataset type, "TC" is the total distribution cost (yuan), "VD" is the total travel distance (km), "ECC" is the energy consumption and charging cost (yuan), "TMC" is the time management cost (including the travel cost and the penalty cost) (yuan), "RC" is the refrigeration cost (yuan), "VL" is the cost of fresh value loss (yuan), "VN" is the number of vehicles, and "CPUT" is the running time (s). relatively significantly. For example, in the R204-type dataset, the TC obtained with the minimization of the PC as the objective is 18.9% higher than that produced by the proposed algorithm, and in the RC204-type dataset, the TC obtained with the minimization of the VD as the objective is 14.6% higher than that produced by the proposed algorithm.

Analysis on the Economic Cost Analysis under the Different Constraints of Fresh Value Loss.
To analyze the relationship between cost elements in the model, we take economic cost and fresh value loss cost as two targets, respectively, to solve the problem. We use the main object method to deal with it. e economic cost is the main object, and the fresh value loss cost is a secondary goal that is converted to the constraint (the ratio of fresh value loss to total value of each customer point should be lower than the minimum threshold). e experiment was tested on the R204 dataset. Table 11 summarizes the test results. In Table 11, "ε" is the minimum threshold for the ratio of the fresh loss value to the fresh value, "TC" is the total distribution cost (yuan), "EC" is the economic cost (yuan), "VL" is the cost of fresh value loss (yuan), "VN" is the number of vehicles, and "CN" represents the number of times EV is fast charging. Table 11 shows the following. (1) As the proportion of fresh value loss decreases, the fresh value loss cost decreases, but the total cost and economic cost increase. Because of the improvement of customers' requirements on the quality of fresh goods, the travel time and the number of customers serviced of each EV are relatively reduced, so more electric vehicles are needed for distribution. (2) When the proportion of fresh value loss is reduced, the number of customers served, travel time, and the travel distance of each EV are reduced; therefore, the probability of en routing charging decreases. In addition, when the proportion is less than 5%, EVs do not need en route charging.

Comparative Analysis of Various Algorithms.
To examine its effectiveness, the proposed algorithm (IACA) is compared with the classical ACA and the tabu search algorithm (TSA). In the classical ACA designed by Devapriya et al. [20] and the TSA designed by Li and Zhang [44], the departure time for all the used vehicles is 0. In the classical ACA, there is no time-window span factor or an adaptive operator in the state transition rule. e R204 and RC204 datasets are used in the test. Table 12 summarizes the results. "EX," "TC," "VD," and "CPUT" are the same as in Table 9. e algorithms take the total distribution cost as the target to solve the EVRP. e following can be derived from Table 12. (1) e proposed algorithm (IACA) outperforms the ACA and the TSA in terms of TC and VD. (2) In regard to computational time, while the proposed algorithm takes into consideration the calculation of the departure time and the improvement of the state transition probability, which are not considered in the ACA and the TAS, there is no significant difference in computation time of the three algorithms. e IACA can produce an optimal result within a relatively short time. (3) Introducing the adaptive heuristic operator and the expected heuristic factor effectively improves the local search and global convergence performance of the ACA.

Conclusion
e characteristics of EVs used in urban cold chain distribution of fresh products are analyzed. On this basis and in view of the time-varying nature of urban road network traffic and the overall power deficiency in EVs, a distribution route optimization model for EVs in urban cold chain logistics for fresh products in time-varying urban traffic networks is established with the minimization of the total cost (including economic cost and fresh value loss cost) as an objective. is model accounts for the differences between urban areas, road types, and queuing times at charging stations. is "min TC" is the minimization of the total distribution cost (yuan), "min VD" is the minimization of the total travel distance (km), and "min PC" is the energy consumption (kWh). e meanings of "EX," "TC," "VD," and "ECC" are the same as in Table 9.  "TC" is the total distribution cost (yuan), "VD" is the total travel distance (km), and "CPUT" is the running time (s).
model also satisfies load, power, and client time-window constraints. Moreover, a corresponding ACA is designed to solve the model. Furthermore, the state transition strategy and pheromone update are improved. Results obtained from multiple tests show the following. (1) When using EVs in cold chain distribution, it is necessary to reasonably plan their distribution routes based on the actual conditions (e.g., urban traffic network, client's time window, power consumption and battery capacity of EVs, and public charging stations) to rationally avoid congestion time periods and reduce the distribution cost. (2) Reducing the distribution time is the key to reducing the cost of cold chain transportation of fresh products. When using EVs in urban cold chain distribution of fresh products, it is important to reduce charging time and frequency by means such as improving battery capacity and increasing singlecharge mileage, to effectively reduce cost. (3) With the increase of customers' requirement for fresh value loss cost reduction, the total cost will increase accordingly. erefore, it is necessary to scientifically plan the distribution route according to the specific requirements of customers to improve customer satisfaction. (4) e model and improved ACA presented in this study can be used to rationally plan distribution time, effectively avoid traffic congestion periods, and reduce the distribution cost.
However, in the actual operation of urban cold chain distribution of fresh products, logistic enterprises face a more complex environment, such as the fresh commodity diversity, the uncertainty about demand for fresh goods, the dynamic nature of traffic information, multiple distribution centers, multiple EV models, and time-of-use electricity prices.
us, the following several aspects need to be addressed in future systematic in-depth studies: first, in the actual traffic environment, the EVRP problem should consider the dynamics and the speed randomness and study its specific representation. Second, further consideration should be given to different fresh goods damage and refrigeration temperature control.
ird, more practical constraints can be introduced, e.g., conducting a comprehensive analysis of the urban cold chain distribution of the EVRP problem under different situations including multiple distribution centers, multiple vehicle types, and multimodal transport.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare no conflicts of interest.