Dynamic Inventory Routing and Pricing Problem with a Mixed Fleet of Electric and Conventional Urban Freight Vehicles

: Urban freight transport is essential for supporting our society regarding providing the daily needs of consumers and local businesses. In addition, it allows for the movement of goods, is distributed within urban environments, provides thousands of jobs, and supports economic growth. However, a number of issues are associated with urban freight transport, including environmental impacts, road congestion, and land use of freight facilities that conﬂicts with residential land use. Electric freight vehicles create zero emissions and provide a sustainable delivery system in comparison with conventional freight vehicles. In this study, a novel dynamic inventory routing and pricing problem under a mixed ﬂeet of electric and conventional vehicles was formulated to minimize the total travel and charging costs. The proposed model is capable of deciding on replenishment times and amounts and vehicle routes. We aimed to determine the maximum social welfare (SW) capable of providing an optimal trade-off between the supplier cost and customer delay that uses a mixed ﬂeet of vehicles. Our computational study was conducted on real data generated from a delivery dataset in Tehran. Under the proposed policy with a ﬂeet of only electric vehicles, the SW increased by 3% while the average customer delay reduced by 15% compared with a ﬂeet of conventional vehicles. The results show that the number of served customers and customer delay would be affected by transitioning conventional urban freight vehicles to electric vehicles. Therefore, the proposed delivery system has a signiﬁcant impact on energy savings and emissions.


Introduction
The demand for the transport of goods within cities is rising while significant progress is being made in increasing the number of freight vehicles in daily city logistics operations. Many companies are striving to improve their supply chains by making them more transparent although conventional urban freight transport is a major contributor to greenhouse gases (GHGs) and is associated with substantially increased health risks. Hence, delivery companies in urban logistics must switch to zero-emission-capable road freight in urban areas using electric or hybrid trucks. Electric freight vehicles (EFVs) offer a potential solution for these issues. More precisely, they allow for the movement of goods and are distributed within an urban environment that reduces air and noise pollutions compared with conventional vehicles. However, the actual application of EFVs in urban logistics operations remains restricted.
EFVs become more competitive for logistic companies in order to save operational costs, leading to significantly increasing environmental benefits. However, the implementation of EFVs in urban logistic operations has several operational challenges, including charging, load capacity, and the limited range of energy. Therefore, vehicle routes should be precisely planned as a result of the limited range of vehicles. Our contributions include the following features: (1) We propose a non-myopic dynamic inventory routing and pricing problem under a mixed fleet of electric and conventional vehicles. (2) We introduce a heterogeneous fleet of vehicles and use a realistic energy consumption that considers vehicle speed, cargo load and gradients, vehicle miles traveled (VMT), vehicle emissions, and charging costs. (3) An empirical study is conducted on real data generated from a delivery dataset in Tehran that investigates both fleet sizing (conventional and electric vehicles) and allocation aiming at maximizing social welfare. The results show that the number of served customers and customer delay would be affected by transitioning conventional vehicles to electric vehicles. Figure 1 presents an overview of the design of the inventory and delivery problem. The supply-chain optimization problems (e.g., the inventory routing problem) aiming at optimal joint decisions of inventory management, vehicle routing, and dispatching of urban freight vehicles have recently received considerable attention. In general, there is a set of pick-up nodes (suppliers) where a set of delivery nodes (customers) should be served by a mixed fleet of electric and conventional vehicles. The products are transported using a fleet of vehicles and customers gradually consuming the products. Pallets and boxes can be used to deliver the products from supply points to demand points. The items can be delivered by each vehicle that has a limited capacity. The type of truck route that must be taken depends on the origins and destinations of deliveries. Due to the activities of private logistics operators, three main classifications for road trucks based on weight are heavy, medium, and light trucks. Heavy-, medium-, and light-duty trucks are vehicles with a gross vehicle weight rating (GVWR) over 26,001 lbs, between 10,001 and 26,000 lbs, and under 10,000 lbs, respectively. This study focuses on a set of trucks to transport goods to a set of geographically scattered retailers. EFVs become more competitive for logistic companies in order to save operational costs, leading to significantly increasing environmental benefits. However, the implementation of EFVs in urban logistic operations has several operational challenges, including charging, load capacity, and the limited range of energy. Therefore, vehicle routes should be precisely planned as a result of the limited range of vehicles. Our contributions include the following features: (1) We propose a non-myopic dynamic inventory routing and pricing problem under a mixed fleet of electric and conventional vehicles. (2) We introduce a heterogeneous fleet of vehicles and use a realistic energy consumption that considers vehicle speed, cargo load and gradients, vehicle miles traveled (VMT), vehicle emissions, and charging costs. (3) An empirical study is conducted on real data generated from a delivery dataset in Tehran that investigates both fleet sizing (conventional and electric vehicles) and allocation aiming at maximizing social welfare. The results show that the number of served customers and customer delay would be affected by transitioning conventional vehicles to electric vehicles. Figure 1 presents an overview of the design of the inventory and delivery problem. The supply-chain optimization problems (e.g., the inventory routing problem) aiming at optimal joint decisions of inventory management, vehicle routing, and dispatching of urban freight vehicles have recently received considerable attention. In general, there is a set of pick-up nodes (suppliers) where a set of delivery nodes (customers) should be served by a mixed fleet of electric and conventional vehicles. The products are transported using a fleet of vehicles and customers gradually consuming the products. Pallets and boxes can be used to deliver the products from supply points to demand points. The items can be delivered by each vehicle that has a limited capacity. The type of truck route that must be taken depends on the origins and destinations of deliveries. Due to the activities of private logistics operators, three main classifications for road trucks based on weight are heavy, medium, and light trucks. Heavy-, medium-, and light-duty trucks are vehicles with a gross vehicle weight rating (GVWR) over 26,001 lbs, between 10,001 and 26,000 lbs, and under 10,000 lbs, respectively. This study focuses on a set of trucks to transport goods to a set of geographically scattered retailers. In addition, the present study investigates a mixed fleet of conventional and EFV freight vehicles that must deliver orders to a set of geographically scattered customers over a long-term horizon look-ahead. The supplier takes on the responsibility of product deliveries and inventory management to customers. Further, a dynamic inventory routing In addition, the present study investigates a mixed fleet of conventional and EFV freight vehicles that must deliver orders to a set of geographically scattered customers over a long-term horizon look-ahead. The supplier takes on the responsibility of product deliveries and inventory management to customers. Further, a dynamic inventory routing problem is developed between supply and retailer nodes while considering an elastic demand (i.e., a demand that depends on inventory, routing, and pricing). A non-myopic approximation model for the specific case of a closed-loop supply chain problem is studied that considers future states. The policy aims to determine the optimal prices and the customers' arrival rate in order to maximize social welfare (SW), which is computed as the Sustainability 2021, 13, 6703 3 of 16 difference between the revenue and the supplier and customer costs. The advantage of this method is its "self-regulating" nature since it encompasses different relationships between the effective arrival rate, replenishment times and amounts, vehicle route, dynamic price, and control service capacity that needs automation. These automated systems can be used in fully automated freight vehicles for a delivery system [1]. The remaining sections of this study are organized as follows: Section 2 presents a relevant literature review on the topic to reveal the contribution of the research. A novel design of dynamic routing and the delivery of items under a mixed fleet of trucks is presented in Section 3, followed by an analysis of numerical examples and simulations of a real case in Section 4. Finally, Section 5 provides the concluding remarks.

Literature Review
Electric freight vehicles are much more energy efficient compared with conventional freight vehicles, and, therefore, are expected to bring considerable energy savings to urban freight. For example, [2,3] studied several approaches to solve the EVR under partial recharge and multiple recharging technologies for using EFVs. Similarly, [4] proposed a vehicle routing problem related to the use of EFVs under a limited capacity of the battery. Likewise, [5,6] investigated a delivery problem for urban trucks that considers power consumption for managing a limited battery capacity. The authors of [7] proposed a static VRP model under a mixed fleet of electric and conventional trucks while considering linear energy consumption. The authors of [8] proposed a mixed-integer linear programming method for a delivery routing problem with a mixed fleet of electric and conventional vehicles. The authors of [9] also introduced an electric vehicle routing problem (EVRP) model that evaluates slow and fast recharging stations by a nonlinear charging approach. Furthermore, [10] studied an EVRP model with recharging stations under a mixed fleet of vehicles. Moreover, [11] extended the model that focuses on a combined fleet of vehicles, such as conventional, plug-in hybrid, and electric vehicles.
Many existing studies have also considered the design of recharging infrastructures in the EVRP. For example, [12] studied a routing-location problem under stochastic demands aiming at minimizing the total cost including the construction cost of the recharging station and the expected traveling cost of the route. In their study, the authors of [13] investigated a two-stage mixed-integer linear program for an EVRP with a stochastic waiting time at recharging stations. Additionally, [14] proposed a Markov decision process (MDP) for an EVRP with a stochastic waiting time at public recharging stations. In another study, [15], the authors introduced a robust optimization location-routing framework that considers the problem of simultaneously routing vehicles and locating charging stations.
To the best of our knowledge, most studies in this area are related to the EVRP in the context of stochastic decision-making problems. For instance, [16] presented a stochastic EVRP without en-route recharging under a stochastic power consumption function. Strategies that consider future states are defined as non-myopic policies. However, a limited body of research has addressed the dynamic routing of electric vehicles under a non-myopic policy. For example, [17] evaluated online EVs with battery swaps as an MDP using an approximated dynamic programming method with linear temporal differencing. The authors of [18] introduced a new social-based optimization model for the distribution of batteries, where the approach comprises decisions related to loaded and unloaded batteries, routing, and scheduling of vehicles between the battery charging and battery swapping stations. A summary of recent contributions to routing and delivery systems is described in Table 1.
This study proposes a new routing and scheduling model that considers replenishment times and amounts under a mixed fleet for urban freight transportation problems, stemming from the work of [18]. Currently, there are no dynamic models with a look-ahead policy that incorporate heterogeneous electric freight vehicles and use energy consumption for the problem analysis. In addition, no model has so far investigated both fleet sizing (conventional and electric vehicles) and allocation aiming at maximizing SW or considered a dynamic dispatch system by taking into account the future state that minimizes the total travel and charging costs. Using the real-time urban freight dataset, we show that both customers and system providers are satisfied by the proposed methodology. Table 1. Summary of routing and delivery studies.

Studies
Model Feature(s) [19] Incorporated energy consumption in a myopic electric routing problem while not considering vehicle loads.
[20] Introduced mixed-integer linear programming for the probabilistic inventory routing problem with a homogenous vehicle that considers emissions and energy consumption.
[16] Minimized the total costs of a myopic electric routing problem. [21] Introduced mixed-integer linear programming for a production inventory and delivery problem with a heterogeneous vehicle while not considering energy consumption and environmental concerns. [18] Maximized social welfare for the non-myopic distribution of the electric battery problem with homogeneous vehicles while not considering energy consumption, vehicle loads, and environmental concerns.

Current paper
Integrated the non-myopic dynamic optimization of routing, transportation, and inventory with a heterogeneous fleet of vehicles that considers energy consumption and environmental impacts.

The Mathematical Model
The proposed algorithm starts by determining the load size, the effective arrival rates, the reward upon completion of service, and the degree of look-ahead for each new customer, followed by updating the locations and service states of the tours determined for the previous analyzed customer. Then, an optimization formulation Problem 1 is applied for each truck to obtain a potential tour, and the optimal arrival rate of the current customer is derived, along with the system revenue. Finally, the truck maximizing SW is chosen, and the potential tour is updated with the selected one while keeping all the others the same as before. In summary, there are three simultaneous decisions to make, including how much to deliver when serving a customer under a mixed fleet of conventional and electric freight trucks, how to route the vehicle among the customers to be served, and when to serve a customer and the amounts. Table 2 presents the input and outputs of Problem 1. Table 2. Input parameters and outputs of the proposed model.

Symbol Description
Indices and sets i, j vertex the fuel-to-air mass ratio p supply nodes χ the heating value of a typical diesel fuel t time periods ς the conversion factor N the set of all points in the network = 0.5C d ρA S the set of the supply nodes, S ⊂ N C d expected air drag coefficient D the set of the demand nodes, D ⊂ N ρ expected air density τ the current location of the trucks, τ ∈ N A frontal area e the dummy end node that all trucks will go to at the end of their tour, e ∈ N γ γ = (1/1000Ψ) T the set of needed steps to finish the delivery efficiency parameter for diesel engines O i the number of orders of customer nodes i ∈ D π π = g sin(α) + gC r cos(α) S i maximum supply nodes i ∈ S C r rolling friction coefficient K the capacity of the truck α road angle (rad) l 0 an initial load of the truck g gravitational constant d ij the distance between points i and j ∈ N ϑ curb-weight (kg) ∆ ij truck speed (km/h) w the weight of a full load of items η η = LH Γ 1 the energy consumption cost of a conventional vehicle L the engine displacement Γ 2 the energy consumption cost of an electric vehicle type of fuel f = {1, 2}, where 1 refers to conventional freight vehicles and 2 refers to electric freight vehicles, v f = 1, . . . , V f Decision variables x ijt if the truck goes from i ∈ N to j ∈ N at step t 1; otherwise 0 y ijt the number of loaded items on the truck going from i ∈ N to j ∈ N at step t l it the number of loaded items in i ∈ S by the truck at time period t u it the number of unloaded items in i ∈ D by the truck at time period t

Energy Consumption
The total amount of fuel for the conventional freight vehicle C is calculated by the engine module, which is linear in the travel cost and investigates vehicle speed and cargo load, and gradients [22][23][24][25]. Table 3 presents the notations and parameters.
Likewise, the total amount of the energy consumption of the road segment for electric freight vehicle E is computed by [16,23,26]

Incorporating Uncertainty
In this section, a value function in terms of SW under the infinite-horizon look-ahead is determined using the demand function that considers a one-to-one connection between non-myopic pricing and effective customers joining the system, the supplier costs (in terms of the length of tours under a mixed fleet of urban trucks), and the expected customer delay. In this study, an approximation method is proposed to investigate future states by examining the trade-off between "exploitation" and "exploration" policies [27,28]. This study applies an approximation approach according to [18]: where SW λ n ;ξ * n ,l v * ,û v * ,ŷ v * ,p * n is the SW function and v f represents the set of vehicles by sizes v and the type of fuels f = {1, 2}, where 1 refers to conventional freight vehicles and 2 refers to electric freight vehicles, v f = 1, . . . , V f . Let χ indicate the set of decision variables that contains an optimal decision on replenishment times and amounts wherê l v * denotes the optimal number of loaded items at supply nodes andû v * is the optimal number of unloaded items delivered to customer nodes. In addition,ŷ v is the number of loaded items on the vehicle. Further,ξ * n is a new tour upon the arrival of a customer and places a new request for new items, and, finally,p * n represents the optimal price [18].
where T ξ n ,l v ,û v ,ŷ v is the length of a tour by vehicle v, which is determined by solving time, which is obtained from the simulation module regarding the potential toursξ n . A weight differentiating the objectives is selected by θ, 0 ≤ θ ≤ 1, and β is the policy parameters that differentiate among non-myopic policies and myopic policies.

Mathematical Formulation
An optimization formulation is run to obtain potential tours. It aims to minimize the traveled distance by the mixed fleet of vehicles to visit customers in a distribution and transportation network. The objective function (8) of this study is to minimize the total charging costs, composed of two terms. The first and the second terms correspond to the sum of the total charging costs under the conventional freight vehicle and the electric freight vehicle, respectively.

Problem 1: Integrated optimization of routing and inventory with a mixed fleet of electric and conventional vehicles
Subject to: Sustainability 2021, 13, 6703 The fuel rate of electric vehicles for a path is presented by constraint (9). The flow equivalence in the nodes is guaranteed by constraints (10)- (15). Constraints (17) and (20) cover the case with the loading and unloading of items in the demand/supply nodes. The dummy node is defined to remove the necessity of returning a truck to the initial point. Constraint (18) guarantees that the truck will not pick more items than what is needed for demand points in its scheduled route. Constraints (16)- (21) guarantee that the load of the truck and all picked items from a supply node will not violate the capacity of the truck limitation at the supply node, respectively. The provision of the needed demand for each customer node is stated in constraints (22) and (23). Finally, the type of variable is defined by constraint (24).
The following algorithm, called Algorithm 1, is applied in the online operational routing and delivery problem (and for simulated runs). The status of each truck location regarding the supply and demand nodes is updated upon the arrival of a new request. The remaining tour (ξ) shows the status of a mixed fleet of trucks v ∈ {V 1 , V 2 }, where V 1 denotes a set of conventional trucks and V 2 refers to a set of electrical trucks. The initially loaded items in the truck and the currently loaded items are determined by l 0 and l, respectively. The set of all nodes in network N is determined by supply points in the route of the truck (dem(v.ξ)), the new demand points (n), and the dummy end node (e). The distance between points d ij is determined where the distance between points and the dummy nodes is zero. Then, the inventory and routing with the mixed fleet model are solved and are shown by constraints (8), along with constraint (24) for the given parameters of the truck v, the nodes N, the capacity of supply nodes C, and the demand of customers O. The SW function is maximized to find the best decisions, such as the best truck (v * ) and the optimal delivery items and routes between supply and demand nodes. Eventually, the outputs include the effective optimal arrival time, fare prices, and the sequence of customers visited by each truck and the tour length.

Algorithm 1. Integrated optimization of vehicle routing, transportation, and inventory with the mixed fleet of vehicles
Step 0. Initialization: Input parameters. For i ∈ n For v ∈ {V1, V 2 }: Step 1. Update

Results
This section describes our tests of the proposed dynamic model. Section 4.1 explains the generation of our test data. In Section 4.2, the efficiency of our proposed approach is demonstrated in a real-life case study, and scenarios are created in a distribution company in Tehran. The main point of the proposed model is an online delivery of items from supply points to demand points under a mixed fleet of trucks. The branch and bound (B&B) algorithm was implemented using Matlab, running on a desktop with an Intel ® Core™ i5-8550U processor with 16GB of RAM and a 64-bit platform using a Windows 10 operating system.

Data Collection
For testing, our computational study was run on real data generated from the daily activity of a distribution company in Tehran. The selected area includes three suppliers (S 1 , S 2 , S 3 ) that serve 13 demand nodes (D 1 . . . , D 13 ). The distribution company is presently one of the largest delivery companies in Iran. It blends, packs, sells, and deliveries teas as well as a wide variety of other types of food products, including rice, spices, pistachios, saffron, and so forth. The effects of other types of urban freight vehicles, such as refrigerated trucks, on the charging cost and energy consumption will be modeled in a future study. Figure 2 shows the distribution of items between the supply and demand networks.
The travel cost between road intersections was determined by the Euclidean distance method. The customer arrival rateλ of items at each demand node appears as the Poisson process and was computed by the historical dataset of a distribution company. The service rateμ of items at each supply node was assumed to follow an exponential distribution and was calculated by the service times from the historical dataset. Table 3 presents the applied input parameters in our simulations.

Experimental Results
In the following section, the routing and delivery system is examined under the fleet of only electric vehicles or the fleet of only conventional vehicles, followed by evaluating the impact of feet size and demand rates on the performance of the distribution model.

Data Collection
For testing, our computational study was run on real data generated from the daily activity of a distribution company in Tehran. The selected area includes three suppliers ( , , ) that serve 13 demand nodes ( … , ). The distribution company is presently one of the largest delivery companies in Iran. It blends, packs, sells, and deliveries teas as well as a wide variety of other types of food products, including rice, spices, pistachios, saffron, and so forth. The effects of other types of urban freight vehicles, such as refrigerated trucks, on the charging cost and energy consumption will be modeled in a future study. Figure 2 shows the distribution of items between the supply and demand networks. The travel cost between road intersections was determined by the Euclidean distance method. The customer arrival rate of items at each demand node appears as the Poisson process and was computed by the historical dataset of a distribution company. The service rate ̂ of items at each supply node was assumed to follow an exponential distribution and was calculated by the service times from the historical dataset. Table 3 presents the Then, the impact of the feet size on the performance of the delivery system was studied. Fleet sizes V are considered in the range 9-15. Next, the VMT and the amount of emissions were determined for each fleet size. Table 5 presents the results of vehicle routes, the loading/unloading of items, and times under the fleet of only conventional vehicles. The VMT for the only conventional vehicles scenario is 175.962, 89.392, 155.970, 34.396,  55.590, 66.512, 94.751, 70.937, 51.937, and 173.317 km, respectively. Thus, the proposed model has a significant impact on energy savings and the environment. Table 4. The results of only electric vehicles. Route: S3-D3-S3-D4-S3-D8-S1-D8-S1-D8-S1-D8-S1-D8-S1-  Route: S3-D11-S2-S2-D10-S1-D10-S1-D10-S1-D10-S1-D10-S1-D10-S1-D10  Route: S1-D9-S2-D9-S2-S2-D10-S1-D10-S1-D10-S1-D10 Time   Route: S3-D4-S3-D3-S3-S1-D10-S1-D10-S1-D10-S1-D10-S1-S10-S1-D10-S1-D10-S1-D10-S1-D10-S1-D10 Time

V 5
Route: S2-D13-S1-D8-S1-D8-S1-D8-S1-D9-S2-D9 Time The application of electric trucks is extremely slow since the operational cost and user cost remain extremely high in comparison with diesel or petrol cars. Using the proposed model, the customer delays would not be affected by transitioning conventional urban freight vehicles to electric vehicles. We compared the time that a customer is waiting to receive their demands with electric and conventional vehicles. The average customer delay for the two types of fleets for the delivery of items is shown in Figure 3. The results show that the average customer delay for electric vehicles and conventional vehicles is 52.54 min and 62.39 min, respectively. The optimal price that a customer is willing to pay is presented in Figure 4. The results show that the socially optimal prices are sensitive to customer delays due to service demands with electric and conventional vehicles. Therefore, the proposed model proposes a competitive price to the customer that considers the customer waiting time and the value of time.   The optimal price that a customer is willing to pay is presented in Figure 4. The results show that the socially optimal prices are sensitive to customer delays due to service demands with electric and conventional vehicles. Therefore, the proposed model proposes a competitive price to the customer that considers the customer waiting time and the value of time. The study also focused on the type of fleets used to distribute items between supply and demand nodes and evaluated the solutions considering different fleet sizes. Then, the impact of the two types of fleets on vehicle miles traveled (VMT) was determined, which is an interesting result for delivery companies. Accordingly, guidelines are provided to make the right fleet sizing decisions according to different fleet size policies regarding the minimization of SW. The levels of the emissions of GHGs due to the delivery of items between supply and demand nodes were computed using the results of our simulated runs: Emissions = VMT × Emission factor (g/km) (25) NOx emissions average 0.225 g/km while CO2 emissions average 450.49 g/km. Table  6 summarizes the performance differences of the two types of fleets for the delivery of items. More precisely, it presents the results of electric and conventional fleets, as expressed by the SW, the VMT, the amount of emissions in terms of NOx and CO2 emissions, and the customer delay. Under the results of only electric vehicles when the fleet size is 10 (Prob 2), the SW increases by 0.28%, and the average customer delay reduces by 15% compared with a delivery system with a fleet of conventional vehicles. Therefore, it is proven that the number of served customers and customer delays would not be affected by transitioning conventional urban freight vehicles to electric vehicles.
In order to consider the pollution from the power charging system, we have to study different types of recharging methods. For instance, battery swapping is an efficient and fast recharging method enabling drivers to go to a battery swapping station (BSS) and replace their empty batteries with full ones. The battery charging stations (BCS) are responsible for distributing full batteries to the BSS nodes. However, batteries can be charged at BCSs with renewable energy sources [18,31]. The study also focused on the type of fleets used to distribute items between supply and demand nodes and evaluated the solutions considering different fleet sizes. Then, the impact of the two types of fleets on vehicle miles traveled (VMT) was determined, which is an interesting result for delivery companies. Accordingly, guidelines are provided to make the right fleet sizing decisions according to different fleet size policies regarding the minimization of SW. The levels of the emissions of GHGs due to the delivery of items between supply and demand nodes were computed using the results of our simulated runs: NOx emissions average 0.225 g/km while CO 2 emissions average 450.49 g/km. Table 6 summarizes the performance differences of the two types of fleets for the delivery of items. More precisely, it presents the results of electric and conventional fleets, as expressed by the SW, the VMT, the amount of emissions in terms of NOx and CO 2 emissions, and the customer delay. Under the results of only electric vehicles when the fleet size is 10 (Prob 2), the SW increases by 0.28%, and the average customer delay reduces by 15% compared with a delivery system with a fleet of conventional vehicles. Therefore, it is proven that the number of served customers and customer delays would not be affected by transitioning conventional urban freight vehicles to electric vehicles.
In order to consider the pollution from the power charging system, we have to study different types of recharging methods. For instance, battery swapping is an efficient and fast recharging method enabling drivers to go to a battery swapping station (BSS) and replace their empty batteries with full ones. The battery charging stations (BCS) are responsible for distributing full batteries to the BSS nodes. However, batteries can be charged at BCSs with renewable energy sources [18,31].
In our experiment, we analyzed the performance of different density rates of customer arrivals at supply nodes for the base case with respect to four key performance indicators with a fleet of 10 vehicles, including VMT, customer delay, SW, and the amount of emissions. Table 7 provides the results of different density rates with only electric and only conventional fleets. Under a high density (0.1), the VMT and customer delay decrease by 8 and 7%, respectively, while the SW increases by 3% in comparison with a fleet of only conventional vehicles. Therefore, the results reveal that the proposed model provides better performance under the high density of demands, because each vehicle can use its capacity. However, the conventional vehicles under a lower density have better performance.

Conclusions
Urban freight vehicles in cities increase traffic-air and noise-pollution issues. Gases emitted by combustion engines cause illnesses and deaths and have a significant impact on economics and public health. However, the application of zero-emissions traffic is extremely slow (i.e., less than 1% of the overall automotive market) since costs remain extremely high in comparison with diesel or petrol cars. Electronic urban vehicles play a key role in reducing the greenhouse effect and saving the driver's energy expenditure. This study proposed a novel dynamic non-myopic model for deciding about the allocation of a mixed fleet for the delivery of items between supply and demand nodes. For this purpose, it was shown how to assign incoming customer requests to fleet types, derive an optimized schedule for the delivery service, and reroute the vehicles as a dynamic non-myopic dynamic distribution system. Accordingly, a dynamic system was provided for integrating routing, transportation, and inventory problems using a mixed fleet of electric and conventional trucks. Additionally, a look-ahead policy was used in the SW function to improve the opportunity costs by considering future states. These components are embedded within a simulation system to evaluate the performance of various fleet for serving customer demands. To investigate the proposed strategy, an experimental example was generated from a delivery system in Tehran. Our experiments indicated substantial improvements in energy and the environment over simpler methods.
For future research, a dynamic model will be proposed to guide trucks by finding an optimal route that considers traffic congestion in the road network. In addition, an energy consumption function for electric vehicles can be developed to manage the limited battery capacity. Many cities around the world are starting to discuss what automated vehicles will mean to their land use, traffic management, and transit planning. A dynamic model can be considered for fully automated freight vehicles in future research that considers joint decisions on charging location, the scheduling of charging, routing, and inventory for a delivery system. Finding the optimal hub location and parking locations for urban freight network design can be considered in future research [32,33]. The impacts of other types of vehicles, such as refrigerated trucks, on energy consumption and charging costs can be obtained and modeled in a future study [34]. Eventually, some customer nodes will be able to receive items indirectly using a nearby node drop, so a flexible routing and delivery system with elastic logistics can be considered in future research [35]. Funding: This research did not receive any grant from funding agencies in the public, commercial, or not-for-profit sectors.
Institutional Review Board Statement: Not applicable.