Next Article in Journal
Technological Solutions and Tools for Circular Bioeconomy in Low-Carbon Transition: Simulation Modeling of Rice Husks Gasification for CHP by Aspen PLUS V9 and Feasibility Study by Aspen Process Economic Analyzer
Previous Article in Journal
Influence of Thermal Aging on the Combustion Characteristics of Cables in Nuclear Power Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Efficient Hybrid Algorithm for Energy Expenditure Estimation for Electric Vehicles in Urban Service Enterprises

Faculty of Transport, Warsaw University of Technology, 00-662 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Energies 2021, 14(7), 2004; https://doi.org/10.3390/en14072004
Submission received: 21 February 2021 / Revised: 19 March 2021 / Accepted: 26 March 2021 / Published: 5 April 2021

Abstract

:
The article deals with the decision problems of estimating the energy expenditure of low-emission fleets in urban service companies due to environmental safety. One of the most important problems of today’s transport policy of many city authorities is the ecological safety of its inhabitants. The basic measures are aimed at banning high-emission vehicles from city centers and promoting the introduction of zero-emission vehicles, such as electric or hybrid cars. The authors proposed an original approach to the decision model, in which the energy expenditure from the use of electric vehicles was defined as a criterion function. The boundary conditions took into account limitations typical of an electric vehicle, e.g., maximum range or battery charging time. To solve the problem, the authors proposed an efficient hybrid algorithm based on ant colony algorithm and genetic algorithm. The verification was made for the example of a utility company serving a medium-sized city in the eastern part of Poland.

1. Introduction

The environmental security of urban life is closely related to sustainability as the ability to meet current needs without compromising future generations [1,2]. Particular attention is paid to the development of ecological means of transport, which are characterized by low emissions of harmful fumes such as electric and hybrid cars. Green mobility or the so-called electric mobility seems to be one of the best solutions to achieve both sustainability goals and mobility needs [3,4]. Many studies point to the effectiveness of electric cars [5] both in terms of the so-called zero local emissions, which is a great advantage for their use in densely populated cities, and in terms of consumers’ willingness to pay for environmental safety [6]. Nevertheless, research is also being conducted on the production costs, autonomy and environmental impact of electric vehicle (EV) production, traffic, and recycling [7].
Due to the restrictions on conventional vehicles entering city centers, a number of analyses are being carried out in service companies on the viability of introducing electric vehicles, including an energy balance analysis using zero or low-emission fleets. Unfortunately, compared to conventional vehicles, electric vehicles require recharging at charging stations due to their limited battery capacity. This implies planning sustainable transport infrastructure in many aspects, including environmental, social, and economic aspects [8]. Moreover, the use of electric vehicles determines the development of modern infrastructure, including charging stations [9]. Ongoing research indicates that a modern and efficient transport system infrastructure both at a macro [10] and micro [11] scale favors the development of both the region of an area, urban centers, or entire countries. Many publications provide methods and tools to support the assessment of pollution and noise emissions [12], assessment of the impact of infrastructure on the quality of life [13], and assessment of the economic viability of a given investment project [14]. The introduction of innovative solutions in the field of transport means also determines the need for research in the field of infrastructural potential at the suppliers providing services to the inhabitants [15]. Moreover, in urban networks, with dense traffic streams and high risk of interference, the probability of collisions is very high. To increase safety, solutions are being introduced to improve traffic flow, e.g., using different ITS (Intelligent Transportation System) service packages [16]. Intelligent transportation systems are information and communication systems that are designed to provide services related to transport and traffic management.
Analysis of the development of urban traffic organization and the introduction of instruments promoting green mobility indicate that electric vehicles are gaining more and more attention as an important step toward alleviating the emission problem [17]. In this context, service companies with tasks in metropolitan areas must adapt their fleets to modern requirements for low-emission vehicles. Due to new regulations related to greenhouse gas emissions in the transport sector, service companies are paying higher fines for every gram of emissions/km. The inclusion of electric vehicles in the service of the recipients of services in the city (customers) affects the reduction of, among others, greenhouse gas emissions or noise emissions.
Urban electric public transport fleets play an important role in e-mobility [18,19] due to a large share of transport service delivery and the obsolete fleets in many companies. The use of electric vehicles, including e-buses, contributes significantly to reducing pollutants in the air [20,21]. Another area of service delivery in urban agglomerations is small parcel delivery services or waste collection services provided by municipal service companies. These companies are a specific type of transport company in which the cargo flow takes place between a collective cargo pick-up or delivery point and the customers [22]. Customers are located in different areas of the city, so tasks in such an enterprise are interpreted as vehicle routes between individual customers and a collective point. Typical service companies operating in a given metropolitan area include municipal companies collecting waste from inhabitants [22], or courier companies [23,24,25]. Transportation by urban service companies is described in terms of the problems of VRP (Vehicle Routing Problem) with Pick-Up and Delivery [26,27,28] in which vehicles visit individual customers, pick up or deliver cargo, and return to a collective point. Vehicle routing problems are problems that consist of determining optimal transport routes taking into account various kinds of restrictions.
This paper proposes a new approach for electric vehicle routing in urban service enterprises in the context of the VRP with Pick-Up and Delivery problem. The authors focused on determining the energy expenditure of fleets in urban service companies served by electric vehicles by applying a hybrid heuristic algorithm. The analysis of the literature highlighted that the application of electric vehicles in urban service enterprises, e.g., courier companies or municipal enterprises, is still a great challenge for researchers.
The whole work is divided into four main parts. The first one is a critical analysis of the literature on methods and tools supporting decision-making on the benefits of using electric cars in an enterprise. The second part deals with the formalization of the original decision model records including the task parameters, boundary conditions, and the criterion for evaluating the quality of the solution. As an evaluation criterion, the authors proposed an energy balance for the analyzed region. The next section presents an original efficient hybrid algorithm based on genetic algorithm [29,30] and ant colony algorithm [31,32,33], which simultaneously determines tasks and assigns them to vehicles. The hybrid combination of algorithms consists of determining by the ant algorithm initial solutions (routes for vehicles), which will be further improved by the genetic algorithm. The ant colony algorithm was used to generate the initial population for the genetic algorithm. The initial population is not possible to be generated in a traditional random way because of the way the chromosomes in this population are represented. Chromosomes contain a sequence of sequentially aligned vertices of a transport network with vehicle route interpretations, while the random selection of these vertices would generate invalid routes. In the light of the literature, these two algorithms are efficient optimization algorithms used to determine NP-hard problems; thus, they have been implemented in the problem of energy expenditure estimation for electric vehicles in urban service enterprises. In addition, the introduction of two optimization algorithms is a mechanism to verify the correctness of the results generated by each algorithm. An important part of the article is the verification of the correctness of the proposed algorithm at the example of a municipal service company dealing with the collection of waste from a given area of the city.

2. Research Status in the Literature

2.1. Electric Vehicles—Development, Benefits, and Weaknesses

There is no doubt that electric vehicles are experiencing a renaissance. This is related to the concept of ecological safety, which is seen through the prism of the development of technologies and means of transport in terms of zero emissions. In 2015, the United Nations adopted the 2030 Agenda for Sustainable Development and presented 17 Sustainable Development Goals (SDGs), reaffirming the World Community’s commitment to sustainable development [34]. Among the 17 objectives, two relate to transportation. One goal points to the need to build resilient infrastructure, promote sustainable industrialization and foster innovation, while the other goal points to the need to make cities socially friendly, safe, and sustainable environments.
Studies on the estimation of exhaust emissions, their extent and magnitude, have been indicated in many publications [35,36]. The European Commission, setting certain limits in the transportation sector in order to significantly reduce emissions, has paid particular attention to zero-emission urban mobility by 2050 [37] and zero-emission urban freight transport by 2030 [38]. In order to reduce CO2 emissions, the need for instruments to encourage the introduction of clean and energy-efficient vehicles for light commercial vehicles by 2025 and 2030 has been identified [39].
Many concerns of organizations for sustainable development indicate that one of the most important issues of modern city managers is electric mobility. This determines the pro-electricity policy in many cities or regions [40]. Electric car purchase incentives are proposed to induce car manufacturers to develop and produce electric vehicles with higher performance standards [41,42]. The problem of eco-mobility, often referred to as electric mobility, is addressed in many studies on car parameters [43], the applicability of quantitative methods and procedures and cost–benefit analysis [44], multi-criteria analysis [45,46,47], and cost-effectiveness analysis [48].

2.2. Routing Problem Including Electric Vehicles

It is a challenge for courier service companies or municipal transport companies and public transport to carry out analyses in order to assess the efficiency of routes and the feasibility of introducing electric vehicles. In this type of enterprise, the task completion problem is a complex decision problem. The task is to route the vehicle from the pick-up point (start of the journey) to the delivery point (end of the journey). First, travel routes (tasks) are determined, and then, in general, vehicles are assigned to these routes. The travel routing problem, combined with the vehicle-to-task allocation problem, is treated as an NP-hard problem (non-deterministic polynomial-time hard problem). NP-hard problems are computational problems for which finding solutions is not possible with polynomial computational complexity. The travel routing problems are related to the traveling salesman problem. In the light of the literature, there is no known effective algorithm that guarantees finding the optimal solution to the traveling salesman problem.
In the paper [49], the authors presented a very in-depth analysis of the state of the art of the Multi-Depot Vehicle Routing Problem (MDVRP), considering models including time windows, split delivery, heterogeneous fleet, periodic delivery, as well as pick-up and delivery. A similar analysis was conducted in the works [50,51]. A rather interesting approach of extending the classical VRP with Hybrid Vehicle Routing (HVRP) is presented in work [52]. The authors mapped the routes using a fleet in which both electric and conventional vehicles were distinguished. The authors assumed that the vehicle could change the drive mode at any time. The solution of Vehicle Routing Problems is carried out based on heuristic or metaheuristic optimization algorithms [53,54].
However, the problem of allocation of vehicles to tasks in service companies is presented in the context of determining the work schedules of vehicles and drivers. In this context, the scheduling process can be interpreted as the process of determining driver work schedules and vehicle schedules [55] or assigning drivers to work shifts. In the literature, these problems are referred to as Bus Driver Scheduling Problems and Crew Scheduling Problems [56,57,58]. In contrast, the allocation of vehicles to transport routes is referred to as the Vehicle Scheduling Problem [59]. In this problem, each vehicle must be assigned routes in the correct order so that the total cost of completing all routes is minimal. The routes interpreted as transportation tasks are known. The problem also lies in determining the minimum number of resources (vehicles, aircraft) to perform all routes.
In the paper [60], the authors presented the problem of routing and scheduling for electric vehicles. They developed an EVRP decision model with constraints relating to capacity, time windows, and pre-defined vehicle charge levels. The problem of determining the routes for electric vehicles taking into account the constraints on battery service life and battery replacement stations was presented by the authors of [61]. To measure the energy consumption and carbon dioxide emissions of electric vehicles, the authors proposed a decision model considering speeds, loads, and distances. They developed an adaptive genetic algorithm based on hill-climbing optimization and neighborhood search to solve the problem. Crossover and mutation probabilities are designed to adaptively adjust to changes in population fitness.
A rather interesting approach to the management of energy consumption was presented by the authors of the paper [62]. The authors presented the problem of vehicle routing in order to find the shortest path as well as the best location for building a storage center and charging station for electric vehicles. The estimated energy consumption of electric buses on the shortest circular route was calculated and compared for electric buses with different energy consumption rates. This makes it possible to plan delivery bus routes in such a way as to minimize the number of battery recharges, whereas energy management strategy for hybrid electric buses (HEBs) based on reinforcement learning (RL) was presented by the authors of the papers [63,64].
An optimization model for electric bus charging schedules that determines both planning and operational decisions while minimizing the total annual cost is presented in [65]. The authors implemented the model using an example of a real transit network in Davis, California. The authors showed that the vehicle operating range problem can be eliminated by adopting certain charging strategies. Comparative analyses have shown that the use of electric buses is more economical and environmentally friendly than diesel buses. On the other hand, the authors of the paper [66] presented the problem of electric vehicle routing with time windows and mixed fleet (E-VRPTWMF). They performed route optimization of a mixed fleet of electric commercial vehicles (ECVs) and conventional internal combustion commercial vehicles (ICCVs) based on an energy consumption model that accounts for speed, slope, and load distribution. To solve the decision problem, they developed an Adaptive Large Neighborhood Search algorithm enhanced with local search.
In the paper [67], the authors pointed out that in order to effectively manage the fleet of electric vehicles, it is necessary to develop algorithms that take into account the charging of vehicles at charging stations. They presented the problem of routing electric vehicles with time windows (E-VRPTW) along with the strategy of multiple or single charging during the route. In doing so, the authors considered a homogeneous fleet of EVs with limited charge and battery capacity, customer time windows, and full linear recharging at charging stations. They used a metaheuristic algorithm based on the ruin–recreate principle to solve the problem.
A rather important problem concerning electric vehicles was pointed out by the authors in paper [68]. In the article, they presented the problem of determining routes for electric vehicles taking into account the decisions concerning routing of vehicles and location of charging stations to support strategic decisions of logistics fleet operators. The authors have considered different charging options due to the constraints present in different decision situations of task execution. They presented alternative objective functions that minimize not only the distance traveled but also the number of vehicles needed and the number of charging stations deployed, as well as the total cost.
Similarly, the authors of the work [69] point to the need for the effective management of the delivery fleet in order to achieve high quality and timely deliveries in logistic distribution processes. The authors described the challenges that have arisen with the integration of electric vehicles into delivery processes, the characteristics of electric vehicles, and the latest energy consumption models. The paper indicates that to address these challenges, it was necessary to adapt effective heuristics and meta-heuristics relevant to the VRP problem. In the field of fleet management for urban public transport, the authors of the paper [70] indicated that the proper design of an electric bus network and fleet planning is important for a bus operator because this determines the purchase of the right number of electric buses at the right time. The authors analyzed the feasibility of operating electric buses as a replacement for conventional buses using the example of Malaysia. They investigated different scenarios of an electric bus operating system using traffic modeling software. The results obtained showed that the developed system of electric bus operation is more efficient with respect to conventional buses not only in terms of generating a higher profit margin for the bus operator but also in terms of better meeting the needs of passengers. The study shows that by using electric vehicles, the operator can carry more passengers per bus unit with lower energy consumption.
Similar conclusions were also reached by the authors of the study [71], who found that diesel buses significantly contribute to air pollution. They proposed a new life additional benefit–cost (NLABC) approach to solve the mixed bus fleet management (MBFM) problem. The NLABC analysis was the basis for developing a total-count program that maximizes the total net benefit of early replacement. The authors show that using this approach, both the optimal fleet size and its composition under budget constraints can be determined. In the comparative analysis, they presented four types of buses, including electric bus, compressed natural gas bus, hybrid bus, and diesel bus.
Having analyzed the literature, it may be stated that most of the publications concerning the routing of electric vehicles refer to enterprises providing passenger transportation services in urban or non-urban areas. The specificity of functioning of municipal service companies carrying out cargo transportation, e.g., courier transport or municipal services [72,73] and the complexity of the NP-hard problem impose the necessity of applying heuristic algorithms for determining the routes of electric vehicles in these enterprises. The routes in enterprises providing passenger transportation services are known and determined by the timetable, whereas the routes in municipal service companies carrying out cargo transportation must be determined. It is advisable to develop new algorithms suitable for to these companies.

3. A Decision for Estimating Energy Expenditure Using Electric Vehicles

3.1. General Assumptions for the Model

Analysis of the following assumptions was introduced to develop a decision model for estimating energy expenditure for electric vehicles in urban service enterprises:
  • A transportation network refers to an urban network where the point elements of the network are intersections and the linear elements are the roads between intersections.
  • A transportation task is interpreted as a travel route of an electric vehicle where the vehicle visits the individual sections of the transport network at which the cargo pick-up or delivery points (collective points) are located. The task starts with the departure of the vehicle from a collective pick-up or delivery point, e.g., a cargo distribution point, and it ends at the same point.
  • Tasks are assigned to vehicles for a given working day.
  • The access to a collective point is possible by various routes; hence, in the model, intermediate points should be determined in the close vicinity of that point indicating the beginnings of various routes of entry and exit to that point.
  • In the case of companies collecting cargo from suppliers, it was assumed that the cargo is located on a given section of the transport network. The volume of cargo in a given section is not greater than the capacity of the vehicle designed to pick it up.
  • A battery charging station can be located either at a collective point or on a particular section of the transport network.
  • The criterion function for determining the energy consumption of vehicles takes into account the energy required to overcome drag and rolling resistance forces. The forces generating the acceleration of the vehicle were not taken into account, as it was assumed that the vehicle moves with uniform motion on the sections. Accelerations due to the start of the vehicle travel have not been taken into account. The vehicle was assumed to be moving on a flat surface, so the uphill drag forces were also neglected. The energy expenditure due to vehicle lighting, air-conditioning etc. is not included in the criterion function

3.2. Model Parameters

The input data of the decision model are shown in Table 1.

3.3. Quantities Sought

For the problem posed in this way, decision variables were defined to determine the travel routes of individual electric vehicles, with the interpretation described in Table 2.

3.4. Restrictions

The model constraints take the following form:
v e h V E H t y p T Y P t a s k T A S K x v e h , t y p , t a s k t w o r k , v e h , t y p v e h V E H t y p T Y P t a s k T A S K p P [ s S i 1 I 1 x 1 s , i 1 , v e h , t y p , t a s k , p d 1 s , i 1 v v e h , t y p , p + i 2 I 2 i 2 I 2 x 2 i 2 , i 2 , v e h , t y p , t a s k , p ( d 2 i 2 , i 2 v v e h , t y p , p + t 1 i 2 , i 2 , t y p ) + i 1 I 1 s S x 3 i 1 , s , v e h , t y p , , t a s k , p ( d 3 i 1 , s v v e h , t y p , p + t 2 s , t y p ) ] + + v e h V E H t y p T Y P t a s k T A S K b B y v e h , t y p , b , t a s k t c h , v e h , t y p , b
v e h V E H ,   t y p T Y P p P t a s k T A S [ s S i 1 I 1 x 1 s , i 1 , v e h , t y p , t a s k , p d 1 s , i 1 v v e h , t y p , p + i 2 I 2 i 2 I 2 x 2 i 2 , i 2 , v e h , t y p , t a s k , p ( d 2 i 2 , i 2 v v e h , t y p , p + t 1 i 2 , i 2 , t y p ) + + i 1 I 1 s S x 3 i 1 , s , v e h , t y p , t a s k , p ( d 3 i 1 , s v v e h , t y p , p + t 2 s , t y p ) + b B y v e h , t y p , b , t a s k t c h , v e h , t y p , b ] t w o r k , v e h , t y p
t y p T Y P v e h V E H t a s k T A S K x v e h , t y p , t a s k n t y p
t a s k T A S K v e h V E H t y p T Y P x v e h , t y p , t a s k = 1
v e h V E H   t y p T Y P ,   t a s k T A S K ,   p P i 2 I 2 i 2 I 2 x 2 i 2 , i 2 , v e h , t y p , t a s k , p q i 2 , i 2 c v e h , t y p
v e h V E H ,   t y p T Y P ,   t a s k T A S K ,   p P , x 1 s , i 1 , v e h , t y p , t a s k , p d 1 s , i 1 + i 2 I 2 i 2 I 2 x 2 i 2 , i 2 , v e h , t y p , t a s k , p d 2 i 2 , i 2 + x 3 i 1 , s , v e h , t y p , t a s k , p d 3 i 1 , s d p e m , v e h , t y p
x 1 s , i 1 , v e h , t y p , t a s k , p d 1 s , i 1 + i 2 I 2 i 2 I 2 x 2 i 2 , i 2 , v e h , t y p , t a s k , p d 2 i 2 , i 2 d p e m , v e h , t y p
Constraint (1) indicates that the disposable time capacity of the number of vehicles used must be greater than the sum of the completion times of all tasks including loading, unloading, and battery charging times of the vehicles. Whereas (2) refers to the completion time of all tasks for a single vehicle, which may not exceed the vehicle’s working time resulting from the driver’s work constraints. The next constraint (3) is the limit on the number of vehicles of a given type, while constraint (4) informs that the current task is performed by exactly one vehicle. The vehicle load capacity constraint (5) is very important, which describes that the load capacity of the vehicle during the task must not be exceeded. This constraint occurs in service companies where the cargo is brought to a collective point. The last two constraints are due to the range of the electric vehicle. The permissible driving distance of the electric vehicle shall not be exceeded, in case the charging station is located at a collective point, Formula (6), or on a section of the transport network, Formula (7) applies.

3.5. Energy Expenditure as a Criterion

The energy expenditure resulting from the transportation tasks was determined as a criterion. Energy expenditure considers the energy used to overcome aerodynamic forces of drag and rolling resistance:
F ( X 1 , X 2 , X 3 ) = = v e h V E H t y p T Y P t a s k T A S K p P [ s S i 1 I 1 x 1 s , i 1 , v e h , t y p , t a s k , p d 1 s , i 1 ( F air + F 1 t ) + i 2 I 2 i 2 I 2 x 2 i 2 , i 2 , v e h , t y p , t a s k , p d 2 i 2 , i 2 ( F air + F 2 t ) + i 1 I 1 s S x 3 i 1 , s , v e h , t y p , t a s k , p d 3 i 1 , s ( F air + F 3 t ) ] min
where:
  • drag force:
    F air = 1 2 ρ c x A ( v v e h , t y p , p ) 2
  • rolling resistance force:
    F 1 t = μ m 1 s , i 1 , v e h g c o s ( α )
    F 2 t = μ m 2 i 2 , i 2 , v e h g c o s ( α )
    F 3 t = μ m 3 i 1 , s , v e h g c o s ( α ) .

4. An Efficient Hybrid Algorithm for Solving the Problem

4.1. Preliminary Assumptions

The hybrid algorithm was developed to determine the energy expenditure of electric vehicles performing transportation tasks in urban service companies. In order to determine the energy expenditure, it is necessary to determine the vehicle travel routes. Two heuristic algorithms have been used to achieve this goal: i.e., ant colony algorithm [74,75] and genetic algorithm [76]. The hybrid combination of these two algorithms involves the initial determination of an initial population of admissible solutions by the ant colony algorithm, which is processed by the genetic algorithm in further steps. The hybrid combination of two heuristic algorithms also serves as a mechanism to verify the correctness of each algorithm and thus the quality of generated solutions. The mechanisms of these algorithms are different, so this increases the plane of finding the final solution.
Ant colony algorithms are algorithms whose working principle is to mimic the existence of ants in the natural environment. The deciding factor in each ant’s route choice is pheromone, which is a chemical secreted by ants. Ants, given a choice of many routes, head for the one with the strongest concentration of pheromone possible.
Genetic algorithms use mechanisms based on the natural phenomena of selection, crossover, and mutation. An adaptation function is introduced to evaluate each individual (chromosome). The selection process involves selecting the best individuals for the next generation (iteration of the algorithm). The crossover process involves the exchange of genetic material between individuals of a population. Mutation involves a random rearrangement of genes, swapping their values or transposition. The crossover and mutation operators are designed to increase the search plane for the optimal solution. The random population generation is selected randomly or by other heuristic algorithms.
Ant colony algorithm and genetic algorithm are iterative algorithms; in successive iterations, the algorithms improve their solution. The algorithm runs until the stop condition is reached. The stop condition is a certain number of iterations. The number of ants in the population as well as individuals (chromosomes) is determined at the beginning of the implementation of the algorithms.

4.2. Operation of the Ant Colony Algorithm

According to the theory, in order to implement the ant colony algorithm, the task to be performed for the ant must first be specified. Its task is to determine the transportation tasks (travel routes) for individual electric vehicles. In order to determine the decision variables described in the mathematical model, the ant creates a route consisting of points defining the type of vehicles selected for the tasks and the number of vehicles of a given type, which are placed in layer 1, and points defining the transportation tasks with the interpretation of connections between intersections, which are placed in layer 2. Graphically, the elements of the ant’s route are shown in Figure 1. The number of elements in layer I depends on the number of available vehicles in the service company, while the number of points in layer II depends on the accuracy of mapping point elements (intersections) in the urban network
The ant starts the route by selecting the type of vehicle and the number of vehicles of the indicated type that will carry out the designated transportation tasks. It is assumed that the ant is a representative of selected vehicles from layer I. When the ant determines the number of vehicles to be used for transportation tasks, the engagement time of the ant is determined, which is the sum of the allowable working times of all selected vehicles. The reselection by the ant of vehicles from layer I occurs when the designated engagement time of the ant is less than the total time of all transportation tasks (constraint of model (1)).
In the next steps, the ant visits individual customers (buyers, suppliers) to deliver the cargo or to pick it up and thus determines the transportation tasks (determination of the decision variables of the model X1, X2, X3, and Y). The start and end point of each transportation task is a collective point of cargo delivery or pick-up from which the ant starts assigning tasks.
An ant goes to the collective point from a route in several situations:
  • When it is fully loaded with cargo collected from customers (constraint of model (5)) or empty when all cargo is delivered to customers,
  • When all the customers have been served,
  • When the battery needs to be recharged at a charging station, if the station is located at a collective delivery point or pick-up point.
The ant also has the ability to recharge batteries at charging stations assigned to particular sections of the transport network. After the ant has mapped the entire route and served all customers (Figure 1), the route is divided into transportation tasks. The task as defined in the mathematical model starts and ends at the collective point of cargo delivery or pick-up, so the tasks in the ant’s route will be sequences of transport network sections defining the ant’s route between the ant’s departure from and return to the collective point.
Next, the transportation tasks are combined into sets of tasks and assigned to a vehicle of a certain type (determination of the decision variable X) in such a way that the permissible operating time of a given vehicle is not exceeded (constraint (2)). The number of task sets created corresponds to the actual number of vehicles used to carry out the transportation tasks.
The starting point of each ant’s route is the collective pick-up or delivery point. The further route of the ant, and thus the selection of subsequent points on the route, occurs with a certain probability:
P R m r y z ( t ) = { [ τ y z ( t ) ] α [ η y z ( t ) ] β l Ω m r [ τ y l ( t ) ] α [ η y l ( t ) ] β , z Ω m r 0 , z Ω m r
where:
  • τ y z ( t ) —intensity of the pheromone trace between y-th point of the ant’s route (crossing) and z-th point (the next crossing) in t-iteration,
  • η y z ( t ) —heuristic information:
η y z ( t ) = 1 d ( y , z ) ,
where:
  • d ( y , z ) —distance between y-th point of ant’s route and z-th point in t-iteration;
  • α , β —effects of pheromones and heuristic data on ant behavior,
  • Ω m r —set of all point elements of the transportation network, l—potential ant route points taken into account at the moment of choosing the next ant route point.
The transition of an ant from the y-th point of the route to the z-th point of the route can be presented in the following steps: step 1—computing the probability of transition from point y to potential points l, step 2—computing the distribution for each connection (y, l) and drawing a number from the interval [0, 1], step 3—choosing the transition route tr with the value of the distribution qtr satisfying the relation qtr − 1 < r qtr, where tr is the number of the route between the y-th point of the route and the z-th point of the route.
Once all ants have finished building routes, the pheromone trail is updated. A cyclic update was used in the process. Initially, the trace is assumed to be equally strong on the links between route points. In subsequent iterations, the pheromone trace is calculated using the formula:
τ y z ( t + 1 ) = ( 1 ρ ) τ y z ( t ) + m r = 1 M R Δ τ y z m r ( t )
where:
  • m r —another ant in the anthill m r M R ,
  • ρ —pheromone volatilization rate ( 0 < ρ 1 ) (0 < ρ 1),
  • τ y z ( t + 1 ) —pheromone reinforcement, for the first iteration takes a value at each connection equal to τ 0 .
The first component of Formula (15) determines the pheromone volatilization rate, while the second component determines the pheromone reinforcement and takes the value (when the route (y, z) was used by the ant mr—(1); otherwise, 0—(2)):
Δ τ y z m r ( t ) = { ( 1 ) 1 D m r ( t ) K m r ( t ) ( 2 ) 0
where:
  • D m r ( t ) —total route traveled in t -iteration by m r -th ant performing all tasks, the route is the basis for the calculation of energy expenditure, criterion function (8);
  • K m r ( t ) —penalty for exceeding the admissible working time by an individual vehicle (model constraint (2)) in the route created by m r -th ant in t -th iteration, the algorithm assumes that the penalty for exceeding is half of the accumulated pheromone on the route;
According to Equation (16), the intensity of the pheromone increases when the denominator of the quotient tends to the minimum value. Wrong routes set by ants are not removed. The introduced penalty function gradually weakens the pheromone at the connection and thus reduces its attractiveness to ants. Removing the entire route could destroy those routes that would generate the best solution in further iterations.
The steps of the ant algorithm can be described as follows:
  • Step 1. Selection of vehicle type and their number by an ant. Based on the number of vehicles, the engagement time of the ant in the transport tasks is determined. When this time is exceeded and all customers (loading or unloading points) have not been served, the ant again indicates the type of vehicle and their number to complete the tasks.
  • Step 2. The ant determines transportation tasks and selects individual points of the route according to the calculated probability (13). In order to avoid searching the whole transportation network for loading or unloading points, at the beginning of each ant’s route, a random order of visiting particular points is generated.
  • Step 3. Before each ant passes to the next point of the route, the distance resulting from the range of the electric vehicle constraint is controlled (6) and (7). After crossing the safety threshold, the ant heads to the charging station and then continues its route in accordance with the drawn order of visiting individual points.
  • Step 4. The capacity constraint (constraint (5)) is controlled before each ant’s transition to the next point of the route. If the constraint is exceeded, the ant heads to a collective pick-up or delivery point.
  • Step 5. Each time the ant returns to the cargo pick-up or delivery point, constraint (1) is checked for the time of the ant’s engagement in the transportation tasks. If the constraint is met, the ant starts creating a new route; if the constraint is not met, the last task (ant’s route) is eliminated, and the ant performs step 1 and chooses again the type and number of vehicles.
  • Step 6. Steps 1–5 are repeated until the ant has created the entire route.
  • Step 7. The next ant in the population starts creating the route. Steps 1–6 are repeated until all routes have been completed by the ants in the population.
  • Step 8. Pheromone update according to Equation (15).
  • Step 9. Steps 1–8 are repeated until all iterations are completed.
  • Step 10. Selection of the ant route with the highest pheromone intensity from all routes generated by ants from the population.
  • Step 11. Allocation of vehicles to a set of tasks within a vehicle type.
When allocating vehicles to tasks, constraint (2) resulting from the completion time of all tasks assigned to a single vehicle should be taken into consideration. The steps of the algorithm designed in this way also guarantee that the constraints arising from the execution of one task by one vehicle (4) are kept and the number of available vehicles (3) is preserved.

4.3. Operation of the Genetic Algorithm

The genetic algorithm is based on the initial population generated by the ant colony algorithm. The genetic algorithm acts as a tool to verify the results generated by the ant colony algorithm. The stages of the genetic algorithm can be represented as follows:
  • Stage 1. Determining the structure processed by the algorithm. A vector structure based on the tasks determined by the ant colony algorithm was used to represent the chromosome. The initial population consists of a finite number of these structures and forms a matrix structure. The number of individuals (vector structures) in the population is fixed at the beginning of the algorithm operation. Genes in the chromosome represent the section numbers of the transport network. Exemplary vector structures processed by the algorithm are shown in Figure 2.
Each row of the structure represents the routes of the ant defining the tasks (exemplary transportation tasks in structure one: (0,3,4,1,2,0), (0,1,3,0), 0-collective point) along with a selection of the type and number of vehicles. The number of columns defining the Type part of the structure depends on the number of vehicle types in the company. The binary value is used to indicate the types of vehicles in the company. As a consequence of the designated vehicle type, the number of these vehicles in the Number section is also indicated. Ant routes are shown in the part of the structure called Tasks.
  • Stage 2. Identification of the adaptation function that evaluates a given structure. Adaptation function for the k-th vector structure M(k, t) can be represented as follows (K = {1, …, k, …, K}—set of structures in the population, t—iteration of the algorithm):
    F ( k , t ) = 1 F ( k , t ) K ( k , t ) m a x
    where:
    • F ( k , t ) —the function defining the length of the route in all transportation tasks; based on this function, the energy expenditure in all transportation tasks is determined;
    • K ( k , t ) —penalty for exceeding the constraint on the completion time of all tasks for a single vehicle (2); the algorithm assumes that the penalty for exceeding halves the generated value of the adaptation function for a given vector structure;
The adaptation function reaches maximum values when the denominator of the quotient tends to minimum values.
  • Stage 3. Selection. The selection process selects the best individuals (chromosomes) from populations, and they are then introduced into the next generation (iteration). The developed algorithm uses the known roulette method, in which linear scaling was applied to counteract premature algorithm convergence in the initial iterations. Scaling coefficient C = 2.0 has been adopted.
  • Stage 4. Crossover. According to the crossover probability parameter, the number of chromosomes to be crossed is selected. The crossover process involves randomly selecting a transportation task from the first chromosome (Figure 3, task: (0,1,2,3,4,0)), identifying genes in the task to which customers are assigned (cargo pick-up and delivery points) e.g., 2 and 3, identifying a task with the same customers in the other chromosome (0,2,6,3,1,0), and exchanging transportation tasks between chromosomes. A prerequisite for the crossover process to take place is that the types of vehicles in both chromosomes are compatible. There is an automatic exchange of the number of vehicles between chromosomes. After the crossover process, it is necessary to check constraint (1) on the time of the ant’s engagement in the tasks. If this constraint is not met, the crossover process is stopped. A graphic interpretation of the crossover process is shown in Figure 3.
  • Stage 5. Mutation. The mutation process involves the exchange of genes in a given chromosome structure. Mutation is infeasible because it generates an unacceptable solution.
The principle of the algorithm can be presented in the following steps:
  • Step 1. Generation of the initial population with the ant colony algorithm,
  • Step 2. Chromosome update according to Equation (17).
  • Step 3. Selection of the best chromosomes for subsequent iterations using the roulette method.
  • Step 4. Crossover.
  • Step 5. Repeating steps 2–4 until the stop condition is reached. The stop condition is the number of iterations specified at the beginning of the algorithm.
  • Step 6. Selection of a chromosome from the population characterized by the maximum value of the adaptation function. This chromosome contains transportation routes characterized by the minimum energy expenditure of electric vehicles.

4.4. The Computational Complexity of the Algorithm

In order to determine the time efficiency of both algorithms, their computational complexity was estimated. The number of operations for each algorithm was determined based on the number of all its loops. The input data for both algorithms include the size of the population (MR), the number of iterations (N), and the number of genes in chromosomes (G). The computational complexity of the ant (AA) and genetic (GA) algorithm was presented in Table 3. In order to estimate the computational complexity, the notation O was determined [77]. The estimation consists in finding a function that limits the function that determines the number of loops in the algorithm. The main loops of algorithms that increase the computational complexity are presented in Table 4. Based on Table 3, it may state that the ant algorithm has the higher computational complexity than the genetic algorithm, which affects the longer operation time of the algorithm. Considering the fact that both algorithms have the polynomial computational complexity and the analyzed problem belongs to NP-hard problems, it may be stated that the generated solutions will be suboptimal solutions.

5. Verification of Algorithms at the Example of an Urban Utility Company

5.1. Model Input Data

The transport network proposed to test the performance of the proposed optimization algorithms is shown in Figure 4. The analyzed network is an urban network, so the point elements of this network are intersections, and additionally, a collective point was introduced to define the waste dump located in the non-urban area, and two vehicle charging points were defined. Municipal waste was located along road sections of the transportation network.
The utility company has nine DAF CF Electric 6 × 2 vehicles, which are characterized by a maximum total weight of 28 tons and load capacity of 7 tons, a battery capacity of 170 kWh, a maximum distance of driving without charging of 100 km in urban conditions, a battery charging time of 1.5 h; and 5 Renault Trucks electric vehicles with a battery capacity of 200 kWh, a range of 120 km of continuous driving in urban conditions, a maximum total weight of 26 tons, a load capacity of 6 tons, and a battery charging time of 2 h. The vehicles work one shift of 8 h. Municipal waste collection begins at 7:00 am. Two traffic peaks were assumed, i.e., the morning peak from 7:00 to 11:00 a.m. and the afternoon peak from 11:00 a.m. to 3:00 p.m., with the following average travel speeds: morning peak—25 km/h, midday peak—35 km/h. The amount of waste in individual sections of the transportation network is shown in Table 5. Waste loading times for a given section are shown in Table 6. Waste discharge time at the landfill was assumed to be 5 min. The distances between intersections in the transportation network are shown in Table 7; the opposite directions are characterized by the same values. The distance between the landfill and the starting points of the route is 0.4 km and 0.5 km.

5.2. Discussion of the Results for the Sensitivity Analysis of the Ant Colony Algorithm

The sensitivity analysis of the ant colony algorithm was carried out for the two variants described in Section 5.1. Both algorithms have been implemented in programming language C#. The following parameter values were assumed to test the ant colony algorithm α : 1; 3; 5; 10; 20, parameter β : 0.5; 1; 5 and parameter ρ : 0.2, 0.4, 0.6, 0.8. The test settings take into account typical values of parameters adopted in the process of algorithm implementation, e.g., α = 1, β = 0.5, ρ = 0.5 [74], and boundary values have been introduced in order to broaden the area of search for the best solution.
A total of 60 possible test settings were generated i.e., (5 (parameter α ) × 3 (parameter β ) × 4 (parameter ρ ) = 60), which were presented in Table 8. Table 9 presents a summary of the results of the algorithm operation.
The energy expenditure was determined according to the criterion function described in the mathematical model (8) based on the pheromone value indicating the minimum travel route of all vehicles (15). Experimentally, the number of iterations of the algorithm was set to 200 repetitions, and the population size was set to 50 ants.
Table 9 shows that the minimum value of the criterion function (8) is F(X1,X2,X3) = 204 [kWh] for test number 24 and algorithm parameters α = 1 , β = 1 ,   ρ = 0.8 . The energy expenditure was determined for three DAF vehicles and one Renault vehicle for a total route of 122 km.
The ant colony algorithm is a probabilistic algorithm, so it can aim for different local optima in each run. For each test (Table 4), the algorithm was run 30 times, and then, the best value for that test was selected. Exemplary results of the ant colony algorithm for the parameters α = 1 , β = 1 ,   ρ = 0.8 are presented in Table 10.
Based on the results presented in Table 9, it may be concluded that the ant colony algorithm is sensitive to parameter changes and directs its search to different regions of the space of admissible solutions. The task of the algorithm calibration process is to indicate such a region in the space of admissible solutions for which the algorithm will generate local optima or will set a global optimum. Graphic interpretation of the ant colony algorithm operation for exemplary settings (test 24, test 25) are shown in Figure 5. The convergence to the optimum starts after more than 100 iterations. Premature convergence of the ant colony algorithm in the initial iterations is minimized by updating the pheromone in each iteration of the algorithm. This action results in the determination of new attractive routes for ants and thus the search for different local optima. Beyond 100 iterations, the pheromone stagnates, and thus, the ants fail to set new solutions.

5.3. Discussion of Results for the Sensitivity Analysis of the Genetic Algorithm

For the sensitivity analysis of the genetic algorithm, the following values of pcross = 0.2; 0.4; 0.6; 0.8; 1.0 were adopted. The mutation parameter was not defined due to the generation of unacceptable solutions after the implementation of this process. The test parameters of the genetic algorithm settings along with the generated energy expenditure are shown in Table 11. The initial population is the final ant population generated in the last iteration of the ant colony algorithm. The minimum energy expenditure in the initial population is 211 [kWh]. Experimentally, the number of iterations of the algorithm was set to 200 repetitions, and the population size is set to 50 chromosomes.
Table 7 shows that the minimum energy expenditure of 199 [kWh] was generated for a crossover parameter equal to pcross = 0.8. The energy expenditure was generated with three Renault vehicles and one DAF vehicle for a total route of 103 km. Results from Table 9 and 11 are the results generated in the last iteration of the algorithm and computed in accordance with criterion function (8) for the chromosome characterized by the maximum adaptation function (17). For small values of the crossover parameter (tests 1, 2), the genetic algorithm replicated the results of the ant colony algorithm because the population was not dominated by better individuals resulting from the crossover process. The genetic algorithm is also a probabilistic algorithm, so it can generate different solutions in each run.
For each test setting, the algorithm was also run 30 times, and then, the best value of the adaptation function for that test was selected and the energy expenditure was determined based on the minimum travel route of all vehicles performing the transportation tasks. Exemplary results of the genetic algorithm for the parameter pcross = 0.8 are presented in Table 12. The convergence to the optimum starts after more than 150 iterations, which is due to the high competition of individuals in the population and thus the emergence of new solutions in subsequent iterations. The consequence of strong competition between individuals in the population is that the algorithm achieves several local optima, Figure 6b.
The calculation time of the ant algorithm in the analyzed case was on average 7.9 min, while for the genetic algorithm, it was 4.3 min. The algorithms were tested on a Dell core i7, 2.20 GHz, 8.00 GB. The genetic algorithm and the ant algorithm can be used to plan vehicle driving routes and determine vehicle work schedules due to the relatively favorable time of generating results.

6. Conclusions

The calibration process of the hybrid algorithm is a key process that affects the correctness and efficiency of the algorithm. This process directs the algorithm into different regions of the space of admissible solutions in order to search for potential local optima or global optima. The parameters for both algorithms were selected experimentally on the basis of the adopted test settings.
For complex decision problems, which is the problem of estimating the energy expenditure of electric vehicles in service companies, heuristic algorithms generate suboptimal solutions, i.e., near-optimal solutions. The probabilistic nature of the algorithms and the search of a limited space of acceptable solutions makes it impossible to determine the global optimum for complex decision problems.
Taking into account the fact that the analyzed problem is an NP-hard problem with many local minimums, it is necessary to use a set of mechanisms that prevent premature convergence of both algorithms to the first encountered local minimum. Premature convergence in the genetic algorithm is blocked by using linear scaling in the roulette method, whereas in the ant algorithm, it is blocked by using a cyclic updating pheromone trail at the beginning of each iteration.
In this paper, 60 test settings for the ant colony algorithm and five settings for the genetic algorithm were examined. The ant colony algorithm reached a local optimum of 211 [kWh]. The genetic algorithm further improved the result; the value of energy expenditure 199 [kWh] appeared three times, so it can be assumed that the solution generated by the genetic algorithm is a suboptimal solution. Considering the results generated by both algorithms, it can be concluded that the combination of these two heuristic algorithms with different operation gives an effective tool for determining the energy expenditure of vehicles in service companies.
For further research on the problem of determining the energy expenditure of electric vehicles in service enterprises, other variations of these algorithms should be applied, e.g., the use of anthill systems. In addition, the problem should be examined in stochastic terms by introducing random variables into the mathematical description.
The authors emphasize the fact that the important aspect of the carried out research is the comparison of the results of the proposed heuristics and their computational times with the result of the exact algorithm [78,79,80,81]. Due to the fact that we have focused strictly on the problem of determining the routes of delivery collection by electric vehicles, the analysis of the algorithm operation time was not presented in detail and can be considered in a future publication.
The results obtained in this paper can serve as the base results for testing other optimization algorithms.

Author Contributions

Conceptualization, M.I.; methodology, M.I., M.J.; software, M.I.; validation, M.I.; formal analysis, M.I. and M.J.; investigation, M.I. and M.J.; resources, M.I. and M.J.; data curation, M.I.; writing—original draft preparation, M.I. and M.J.; writing—review and editing, M.I. and M.J.; visualization, M.I.; supervision, M.J.; project administration, M.J.; funding acquisition, M.I. and M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Warsaw University of Technology within the IDUB-Open Science Program.

Acknowledgments

The authors would like to gratefully acknowledge the reviewers provided helpful comments and insightful suggestions on draft of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. What Is Meant by the Term “Sustainability”? Available online: http://www.fao.org/3/ai388e/AI388E05.htm (accessed on 16 January 2021).
  2. UN Secretary-General. World Commission on Environment and Development: Our Common Future; Oxford University Press: Oxford, UK, 1987. [Google Scholar]
  3. Cieśla, M.; Sobota, A.; Jacyna, M. Multi-Criteria Decision Making Process in Metropolitan Transport Means Selection Based on the Sharing Mobility Idea. Sustainability 2020, 12, 7231. [Google Scholar] [CrossRef]
  4. Westbrook, M.H.; Westbrook, M. The Electric Car: Development and Future of Battery, Hybrid and Fuel-Cell Cars; IET: London, UK, 2001. [Google Scholar]
  5. Pielecha, J.; Skobiej, K.; Kurtyka, K. Exhaust Emissions and Energy Consumption Analysis of Conventional, Hybrid, and Electric Vehicles in Real Driving Cycles. Energies 2020, 13, 6423. [Google Scholar] [CrossRef]
  6. Zhou, G.; Hu, W.; Huang, W. Are Consumers Willing to Pay More for Sustainable Products? A Study of Eco-Labeled Tuna Steak. Sustainability 2016, 8, 494. [Google Scholar] [CrossRef] [Green Version]
  7. Carteni, A.; Henke, I.; Molitierno, C.; Errico, A. Towards E-Mobility: Strengths and Weaknesses of Electric Vehicles. In Proceedings of the Workshops of the International Conference on Advanced Information Networking and Applications, Caserta, Italy, 15–17 April 2020; Springer: Cham, Switzerland, 2020; pp. 1383–1393. [Google Scholar]
  8. Henke, I.; Biggiero, L.; Pagliara, F. The Environmental Risks Related to Visitors’ Trips to Festivals: Transport Planning for Sustainability. In Proceedings of the 20th International Conference on Environment and Electrical Engineering, Madrid, Spain, 9–12 June 2020. [Google Scholar]
  9. Erickson, L.; Ma, S. Solar-Powered Charging Networks for Electric Vehicles. Energies 2021, 14, 966. [Google Scholar] [CrossRef]
  10. Jacyna, M.; Wasiak, M.; Lewczuk, K.; Kłodawski, M. Simulation model of transport system of Poland as a tool for developing sustainable transport. Arch. Transp. 2014, 31, 23–35. [Google Scholar] [CrossRef]
  11. Karoń, G.; Żochowska, R. Problems of Quality of Public Transportation Systems in Smart Cities—Smoothness and Disruptions in Urban Traffic. In Modelling of the Interaction of the Different Vehicles and Various Transport Modes, Lecture Notes in Intelligent Transportation and Infrastructure; Sładkowski, A., Ed.; Springer: New York, NY, USA, 2020. [Google Scholar]
  12. Wasiak, M.; Niculescu, A.I.; Kowalski, M. A generalized method for assessing emissions from road and air transport on the example of Warsaw Chopin Airport. Arch. Civ. Eng. 2020, 66, 399–419. [Google Scholar] [CrossRef]
  13. Cascetta, E.; Carteni, A.; Henke, I. Acceptance and Equity in Advanced Path-Related Road Pricing Schemes. In Proceedings of the 2017 5th IEEE International Conference on Models and Technologies for Intelligent Transportation Systems [MT-ITS], Naples, Italy, 26–28 June 2017; pp. 492–496. [Google Scholar]
  14. Bagloee, S.A.; Asadi, M.; Bozic, C. A Sustainability Approach in Road Project Evaluation, Case-Study: Pollutant Emission and Accident Costs in Cost Benefit Analysis. In Sustainable Automotive Technologie; Springer: Berlin/Heidelberg, Germany, 2012; pp. 295–303. [Google Scholar]
  15. Izdebski, M.; Jacyna, M. The organization of the municipal waste collection: The decision model. Rocz. Ochr. Sr. 2018, 20, 919–933. [Google Scholar]
  16. Karoń, G.; Żochowska, R. Modelling of Expected Traffic Smoothness in Urban Transportation Systems for ITS Solutions. Arch. Transp. 2015, 33, 33–45. [Google Scholar] [CrossRef]
  17. Jacyna, M.; Kotylak, P. Decision-making problems of collective transport development in terms of sustainable urban mobility. J. KONBIN 2020, 50. [Google Scholar] [CrossRef]
  18. Carteni, A. A Plug-in Hybrid Electric Bus Fleet as a Rational and Sustainable Urban Transport Policy: A Real Casa Application in Italy. Int. J. Transp. Syst. 2017, 2, 46–53. [Google Scholar]
  19. Kühne, R. Electric Buses–An Energy Efficient Urban Transportation Means. Energy 2010, 35, 4510–4513. [Google Scholar] [CrossRef]
  20. Bakker, S.; Konings, R. The Transition to Zero-Emission Buses in Public Transport—The Need for Institutional Innovation. Transp. Res. Part D Transp. Environ. 2018, 64, 204–215. [Google Scholar] [CrossRef]
  21. Guida, U.; Leonard, S. ZeEUS: Zero Emission Urban Bus System. In Proceedings of the 2014 IEEE International Electric Vehicle Conference [IEVC], Florence, Italy, 17–19 December 2014; pp. 1–7. [Google Scholar]
  22. Quak, H.; Nesterova, N.; van Rooijen, T.; Dong, Y. Zero Emission City Logistics: Current Practices in Freight Electromobility and Feasibility in the near Future. Transp. Res. Procedia 2016, 14, 1506–1515. [Google Scholar] [CrossRef] [Green Version]
  23. Han, H.; Ponce-Cueto, E. Waste Collection Vehicle Routing Problem: A Literature Review. Promet Traffic Transp. 2015, 27, 345–356. [Google Scholar] [CrossRef] [Green Version]
  24. Kummer, S.; Hribernik, M.; Herold, D.M.; Mikl, J.; Dobrovnik, M.; Schoenfelder, S. The impact of courier-, express- and parcel (CEP) service providers on urban road traffic: The case of Vienna. Transp. Res. Interdiscip. Perspect. 2021, 9, 100278. [Google Scholar] [CrossRef]
  25. Ducret, R. Parcel deliveries and urban logistics: Changes and challenges in the courier express and parcel sector in Europe—The French case. Res. Transp. Bus. Manag. 2014, 11, 15–22. [Google Scholar] [CrossRef]
  26. Nowak, M.; Ergun, O.; White, C.C., III. Pickup and delivery with split loads. Transp. Sci. 2008, 42, 32–43. [Google Scholar] [CrossRef]
  27. Psaraftis, H. A multi-commodity, capacitated pickup and delivery problem: The single and two-vehicle cases. Eur. J. Oper. Res. 2011, 215, 572–580. [Google Scholar] [CrossRef]
  28. Liu, R.; Xie, X.; Augusto, V.; Rodriguez, C. Heuristic algorithms for a vehicle routing problem with simultaneous delivery and pickup and time windows in home health care. Eur. J. Oper. Res. 2013, 230, 475–486. [Google Scholar] [CrossRef] [Green Version]
  29. Izdebski, M.; Jacyna-Gołda, I.; Gołębiowski, P.; Plandor, J. The Optimization Tool Supporting Supply Chain Management in the Multi-Criteria Approach. Arch. Civ. Eng. 2020, 66, 505–524. [Google Scholar] [CrossRef]
  30. Hubert, J.; Wilck, J.; Cavalier, T.M. A Genetic Algorithm for the Split Delivery Vehicle Routing Problem. Am. J. Oper. Res. 2012, 2, 207–216. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, C.Y.; Yu, J.J. Multiple depots vehicle routing based on the ant colony with the genetic algorithm. J. Ind. Eng. Manag. 2012, 6, 1013–1026. [Google Scholar] [CrossRef] [Green Version]
  32. Yu, B.; Yang, Z.Z. An ant colony optimization model: The period vehicle routing problem with time windows. Transp. Res. Part E 2011, 47, 166–181. [Google Scholar] [CrossRef]
  33. Montemanni, R.; Gambardella, L.M.; Rizzoli, A.E.; Donati, A.V. Ant Colony System for a Dynamic Vehicle Routing Problem. J. Comb. Optim. 2005, 10, 327–343. [Google Scholar] [CrossRef]
  34. United Nations. Transforming Our World: The 2030 Agenda for Sustainable Development; Division for Sustainable Development Goals; UN: New York, NY, USA, 2015. [Google Scholar]
  35. Jacyna, M.; Merkisz, J. Proecological approach to modelling traffic organization in national transport system. Arch. Transp. 2014, 2, 43–56. [Google Scholar] [CrossRef] [Green Version]
  36. Jacyna-Gołda, I.; Żak, J.; Gołębiowski, P. Models of traffic flow distribution for various scenarios of the development of proecological transport system. Arch. Transp. 2014, 4, 17–28. [Google Scholar] [CrossRef]
  37. UNFCCC. Kyoto Protocol to the United Nations Framework Convention on Climate Change Adopted at COP3 in Kyoto, Japan, on 11 December 1997. Available online: https://www.eea.europa.eu/data-and-maps/indicators/primary-energy-consumption-by-fuel/unfccc-1997-kyoto-protocol-to (accessed on 12 January 2021).
  38. Niestadt, M.; Bjornavold, A. Electric Road Vehicles in the European Union; European Parliament: Brussels, Belgium, 2019. [Google Scholar]
  39. Ulrich, P.; Lehr, U. Economic Effects of an E-Mobility Scenario–Input Structure and Energy Consumption. Econ. Syst. Res. 2020, 32, 84–97. [Google Scholar] [CrossRef]
  40. Jakovljevic, B.; Paunovic, K.; Belojevic, G. Road-traffic noise and factors influencing noise annoyance in an urban population. Environ. Int. 2009, 35, 552–556. [Google Scholar] [CrossRef]
  41. Dijk, M.; Iversen, E.; Klitkou, A.; Kemp, R.; Bolwig, S.; Borup, M.; Møllgaard, P. Forks in the Road to E-Mobility: An Evaluation of Instrument Interaction in National Policy Mixes in Northwest Europe. Energies 2020, 13, 475. [Google Scholar] [CrossRef] [Green Version]
  42. Cartenì, A.; De Guglielmo, M.L.; Henke, I. Design of Sustainable Urban Transport Infrastructures: A Real Case Application in Italy. Int. J. Civ. Eng. Technol. 2018, 9, 2131–2147. [Google Scholar]
  43. Merkisz, J.; Jacyna, M.; Merkisz-Guranowska, A.; Pielecha, J. The parameters of passenger cars engine in terms of real drive emission test. Arch. Transp. 2014, 32, 43–50. [Google Scholar] [CrossRef] [Green Version]
  44. Carteni, A. A Cost-Benefit Analysis Based on the Carbon Footprint Derived from Plug-in Hybrid Electric Buses for Urban Public Transport Services. WSEAS Trans. Environ. Dev. 2018, 14, 125–135. [Google Scholar]
  45. Tzeng, G.-H.; Lin, C.-W.; Opricovic, S. Multi-Criteria Analysis of Alternative-Fuel Buses for Public Transportation. Energy Policy 2005, 33, 1373–1383. [Google Scholar] [CrossRef]
  46. Vögele, S.; Ball, C.; Kuckshinrichs, W. Multi-Criteria Approaches to Ancillary Effects: The Example of e-Mobility. In Ancillary Benefits of Climate Policy; Springer: Cham, Switzerland, 2020; pp. 157–178. [Google Scholar]
  47. Watróbski, J.; Malecki, K.; Kijewska, K.; Iwan, S.; Karczmarczyk, A.; Thompson, R.G. Multi-Criteria Analysis of Electric Vans for City Logistics. Sustainability 2017, 9, 1453. [Google Scholar] [CrossRef] [Green Version]
  48. Ernst, C.-S.; Hackbarth, A.; Madlener, R.; Lunz, B.; Sauer, D.U.; Eckstein, L. Battery Sizing for Serial Plug-in Hybrid Vehicles: A Model-Based Economic Analysis for Germany. Energy Policy 2011, 39, 5871–5882. [Google Scholar] [CrossRef]
  49. Montoya-Torresa, J.R.; Franco, J.L.; Isaza, S.N.; Jiménez, H.F.; Herazo-Padilla, N. A literature review on the vehicle routing problem with multiple depots. Comput. Ind. Eng. 2015, 79, 115–129. [Google Scholar] [CrossRef]
  50. Braekers, K.; Ramaekers, K.; Van Nieuwenhuyse, I. The vehicle routing problem: State of the art classification and review. Comput. Ind. Eng. 2016, 99, 300–313. [Google Scholar] [CrossRef]
  51. Laporte, G. Fifty Years of Vehicle Routing. Transp. Sci. 2009, 43, 408–416. [Google Scholar] [CrossRef]
  52. Mancini, S. The hybrid vehicle routing problem. Transp. Res. Part C 2017, 78, 1–12. [Google Scholar] [CrossRef]
  53. Jacyna-Gołda, I.; Izdebski, M.; Podviezko, A. Assessment of efficiency of assignment of vehicles to tasks in supply chains: A case study of a municipal company. Transport 2017, 32, 243–251. [Google Scholar] [CrossRef] [Green Version]
  54. Bianchessi, N.; Righini, G. Heuristic algorithms for the vehicle routing problem with simultaneous pick-up and delivery. Comput. Oper. Res. 2007, 34, 578–594. [Google Scholar] [CrossRef]
  55. Valouxis, C.; Housos, E. Combined bus and driver scheduling. Comput. Oper. Res. 2002, 29, 243–259. [Google Scholar] [CrossRef]
  56. Desrochers, M.; Soumis, F. A column generation approach to the urban transit crew scheduling problem. Transp. Sci. 1989, 23, 1–13. [Google Scholar] [CrossRef] [Green Version]
  57. Haase, K.; Desaulniers, G.; Desrosiers, J. Simultaneous vehicle and crew scheduling in urban mass transit systems. Transp. Sci. 2001, 35, 286–303. [Google Scholar] [CrossRef]
  58. Huisman, D.; Freling, R.; Wagelmans, A.P.M. Multiple-depot integrated vehicle and crew scheduling. Transp. Sci. 2005, 39, 491–502. [Google Scholar] [CrossRef] [Green Version]
  59. Burkhard, R.; Dell’Amico, M.; Marttelo, S. Assignment Problems; Society for Industrial and Applied Mathematics; SIAM: Philadelphia, PA, USA, 2009. [Google Scholar]
  60. Anagnostopoulou, A.; Boile, M.; Theofanis, S.; Sdoukopoulos, E.; Margaritis, D. Electric Vehicle Routing Problem with Industry Constraints: Trends and Insights for Future Research. Transp. Res. Procedia 2014, 3, 452–459. [Google Scholar] [CrossRef] [Green Version]
  61. Li, J.; Wang, F.; He, Y. Electric Vehicle Routing Problem with Battery Swapping Considering Energy Consumption and Carbon Emissions. Transp. Res. Part D Transp. Environ. 2018, 60, 104–118. [Google Scholar] [CrossRef]
  62. Taweepworadej, W.; Buasri, P. Vehicle Routing Problem for Electric Bus Energy Consumption and Planning. Int. J. Adv. Agric. Environ. Eng. 2016, 3, 224–226. [Google Scholar] [CrossRef]
  63. Wu, J.; Wei, Z.; Li, W.; Wang, Y.; Li, Y.; Sauer, D. Battery Thermal- and Health-Constrained Energy Management for Hybrid Electric Bus based on Soft Actor-Critic DRL Algorithm. IEEE Trans. Ind. Inform. 2020, 17, 3751–3761. [Google Scholar] [CrossRef]
  64. Wu, J.; Wei, Z.; Liu, K.; Quan, Z.; Li, Y. Battery-involved Energy Management for Hybrid Electric Bus Based on Expert-assistance Deep Deterministic Policy Gradient Algorithm. IEEE Trans. Veh. Technol. 2020, 69, 12786–12796. [Google Scholar] [CrossRef]
  65. Wang, Y.; Huang, Y.; Xu, J.; Barclay, N. Optimal recharging scheduling for urban electric buses: A case study in Davis. Transp. Res. Part E Logist. Transp. Rev. 2017, 100, 115–132. [Google Scholar] [CrossRef]
  66. Goeke, D.; Schneider, M. Routing a mixed fleet of electric and conventional vehicles. Eur. J. Oper. Res. 2015, 245, 81–99. [Google Scholar] [CrossRef]
  67. Erdelić, T.; Carić, T.; Erdelić, M.; Tišljarić, L. Electric vehicle routing problem with single or multiple recharges. Transp. Res. Procedia 2019, 40, 217–224. [Google Scholar] [CrossRef]
  68. Schiffer, M.; Walther, G. The electric location routing problem with time windows and partial recharging. Eur. J. Oper. Res. 2017, 260, 995–1013. [Google Scholar] [CrossRef]
  69. Erdelić, T.; Carić, T. A Survey on the Electric Vehicle Routing Problem: Variants and Solution Approaches. J. Adv. Transp. 2019, 54, 1–48. [Google Scholar] [CrossRef]
  70. Teoh, L.E.; Khoo, H.L.; Goh, S.Y.; Chong, L.M. Scenario-based electric bus operation: A case study of Putrajaya, Malaysia. Int. J. Transp. Sci. Technol. 2018, 7, 10–25. [Google Scholar] [CrossRef]
  71. Lu, L.; Hong, K.L.; Feng, X.; Xuekai, C. Mixed bus fleet management strategy for minimizing overall and emissions external costs. Transp. Res. Part D Transp. Environ. 2018, 60, 104–118. [Google Scholar] [CrossRef]
  72. Beliën, J.; De Boeck, L.; Van Ackere, J. Municipal solid waste collection and management problems: A literature review. Transp. Sci. 2012, 48, 78–102. [Google Scholar] [CrossRef]
  73. Bautista, J.; Fernández, E.; Pereira, J. Solving an urban waste collection problem using ants heuristics. Comput. Oper. Res. 2008, 35, 3020–3033. [Google Scholar] [CrossRef] [Green Version]
  74. Dorigo, M.; Stutzle, T. Ant Colony Optimization; Bradford Books: Cambridge, MA, USA, 2004. [Google Scholar]
  75. Dorigo, M.; Blum, C. Ant Colony Optimization Theory: A survey. Theor. Comput. Sci. 2005, 344, 243–278. [Google Scholar] [CrossRef]
  76. Goldberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning, 1st ed.; Addison-Wesley Professional: Boston, MA, USA, 1989. [Google Scholar]
  77. Arora, S.; Barak, B. Computational Complexity: A Modern Approach; Cambridge University Press: New York, NY, USA, 2009. [Google Scholar]
  78. Oelschlägel, T.; Knust, S. Solution approaches for storage loading problems with stacking constraints. Comput. Oper. Res. 2021, 127, 105–142. [Google Scholar] [CrossRef]
  79. Dean, B.C. A simple expected running time analysis for randomized “divide and conquer” algorithms. Discret. Appl. Math. 2006, 154, 1–5. [Google Scholar] [CrossRef] [Green Version]
  80. Early, D.; Schellekens, M. Running time of the Treapsort algorithm. Theor. Comput. Sci. 2013, 487, 65–73. [Google Scholar] [CrossRef]
  81. Sudholt, D.; Thyssen, C. Running time analysis of Ant Colony Optimization for shortest path problems. J. Discret. Algorithms 2012, 10, 165–180. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Elements of the ant’s route.
Figure 1. Elements of the ant’s route.
Energies 14 02004 g001
Figure 2. Structure of the genetic algorithm.
Figure 2. Structure of the genetic algorithm.
Energies 14 02004 g002
Figure 3. Crossover process: (a) chromosomes before crossover, (b) chromosomes after crossover.
Figure 3. Crossover process: (a) chromosomes before crossover, (b) chromosomes after crossover.
Energies 14 02004 g003
Figure 4. Transportation network for municipal waste collection.
Figure 4. Transportation network for municipal waste collection.
Energies 14 02004 g004
Figure 5. Operation of the ant colony algorithm, (a) test settings 24, (b) test settings 25.
Figure 5. Operation of the ant colony algorithm, (a) test settings 24, (b) test settings 25.
Energies 14 02004 g005
Figure 6. Operation of the genetic algorithm, (a) test settings 4, (b) test settings 3.
Figure 6. Operation of the genetic algorithm, (a) test settings 4, (b) test settings 3.
Energies 14 02004 g006
Table 1. Symbols used in the decision model.
Table 1. Symbols used in the decision model.
SymbolMeaning
SThe set of collective starting points of a vehicle route, e.g., cargo distribution points, dumps where s is an element of the S set.
I1The set of intermediate points in the vicinity of a collective point, i.e., intersections defining different directions of vehicle travel routes, where i1 is an element of the I1 set.
I2The set of all intermediate points, i.e., intersections in close and distant proximity to the collective point, where i2′, i2″ are elements of the I2 set.
BA set of electric vehicle battery charging stations, where b is an element of the B set.
VEHThe set of vehicles in the service enterprise, where veh is an element of the VEH set.
TYPEThe set of vehicle types in a service company that differ in cargo capacity, where type is an element of the TYPE set.
TASKA set of transportation tasks to be completed on a given working day, where task is an element of the TASK set.
PThe set of periods with interpretations of the morning, midday, and afternoon peaks, where p is an element of the P set.
vveh,typ,pSpeed of a vehicle of a given type during a given period of a working day, where V is the speed matrix.
ntypThe number of vehicles of a particular type in a service company, where N is a vector of the number of vehicles of a particular type.
twork,veh,typThe vehicle operating time resulting from driver work restrictions, where Twork is a matrix of vehicle operating time of a specific type.
qi2′,i2″The amount of cargo ordered to be taken from a given route segment or ordered to be delivered, where Q is the cargo size matrix.
cveh,typThe load capacity of a vehicle of a given type, where C is the vehicle load capacity matrix.
t1i2′,i2″,typThe time of loading or unloading in a given route section, depending on the type of vehicle, where T1 is a matrix of loading or unloading times depending on the services provided by the company.
t2s,typUnloading or loading time at a collective point dependent on vehicle type, where T2 is a matrix of unloading or loading times.
tch,veh,typ,bThe battery charging time of a vehicle of a given type at a given charging station, where Tch is a matrix of battery charging times of vehicles of a given type at a given charging station.
d1s,i1The distance between the collective point and the initial intermediate point of the vehicle route, where D1 is the distance matrix between the collective and the initial intermediate point.
d2i2′,i2″The distance between intermediate points of the transport network, where D2 is the distance matrix between intermediate points of the transport network.
d3i1,sThe distance between the intermediate points of the transport network and the collective point, where D3 is the distance matrix between the intermediate points of the transport network and the collective point.
dpem,veh,typThe allowable travel distance for each electric vehicle, where Dpem is the matrix of allowable travel distances for each vehicle.
m1s,i1,vehThe weight of the vehicle together with load between the collective point and the intermediate points of the route, where M1 is a matrix defining the weight of the vehicle on the indicated route section.
m2i2′,i2″,vehThe weight of the vehicle including load between intermediate points of the route, where M2 is a matrix defining the weight of the vehicle on the indicated route section.
m3i1,s,vehThe weight of the vehicle including load between intermediate points of the route and the collective point, where M3 is a matrix representing the weight of the vehicle on the indicated route section.
gEarth acceleration, where g = 9.81 m/s2.
μThe coefficient of road friction was assumed to be identical on each section of the route, μ = 2.02 × 10−5.
cxThe drag coefficient, where cx = 0.2.
AVehicle frontal area, for each type of vehicle an identical frontal area is assumed, A = 20 m2.
Table 2. Decision variables of the model.
Table 2. Decision variables of the model.
SymbolMeaning
xveh,typ,taskThe assignment of a vehicle of a certain type to a task, where X is a matrix defining the assignment of vehicles to tasks. If xveh,typ,task = 1, there is an assignment of a given vehicle to a task; otherwise, xveh,typ,task = 0.
x1s,i1,veh,typ,task,pA connection between a collective point and a route intermediate point made by a vehicle of a given type during a given period of the day in the task in progress, where X1 is a matrix specifying all these connections. If x1s,i1,veh,typ,task,p = 1, the vehicle of the type performs the route between the collective point and the intermediate point; otherwise, x1s,i1,veh,typ,task,p = 0.
x2i2′,i2″,veh,typ,task,pA connection between intermediate points made by a vehicle of a given type during a given period of the day in the task in progress, where X2 is a matrix specifying all these connections. If x2i2′,i2″,veh,typ,task,p = 1, the vehicle of a given type performs the route between intermediate points; otherwise, x2i2′,i2″,veh,typ,task,p = 0.
x3i1,s,veh,typ,task,pA connection between a route intermediate point and a collective point made by a vehicle of a given type during a given period of the day in the task in progress, where X3 is a matrix defining all these connections. If x3i1,s,veh,typ,task,p = 1, the vehicle of the type performs the route between the intermediate point and the collective point; otherwise, x3i1,s,veh,typ,task,p = 0.
yveh,typ,b,taskMaking a decision to charge a vehicle at a given station in a task in progress, where Y is a matrix specifying all decisions made to charge vehicles at given stations. If yveh,typ,b,task = 1, the vehicle of the given type charges the batteries at the given station; otherwise, yveh,typ,b,task = 0.
Table 3. The computational complexity of the algorithms.
Table 3. The computational complexity of the algorithms.
AlgorithmsInput DataLoopsEstimationComplexity
AAMR, N, G2∙G2∙MR∙N + G∙MR∙N2∙G2∙MR∙N + G∙MR∙N ≤ (MR + N + G)4O((MR + N + G)4)
polynomial
GAMR, N, G3∙G∙MR∙N3∙G∙MR∙N ≤ (MR + N + G)3O((MR + N + G)3)
polynomial
Table 4. The main loops in the algorithms.
Table 4. The main loops in the algorithms.
AALoopsGALoops
Determination of the probabilityG2∙MR∙NSelectionG∙MR∙N
Determination of the distributionG2∙MR∙NCrossoverG∙MR∙N
A cyclic updateG∙MR∙NPopulation assessmentG∙MR∙N
Table 5. Amount of municipal waste [tons].
Table 5. Amount of municipal waste [tons].
q2,3 = 0.4q3,2 = 0.6q4,5 = 0.2q5,4 = 0.6q7,8 = 0.2q8,7 = 0.8
q15,16 = 1.2q16,15 = 0.4q21,22 = 0.4q22,21 = 0.2q24,25 = 0.6q14,23 = 0.4
q19,28 = 0.6q29,30 = 0.2q32,31 = 0.4q26,35 = 0.4q44,45 = 0.6q45,44 = 0.2
q42,41 = 1.2q55,54 = 0.4q52,51 = 0.6q51,50 = 0.2q49,48 = 1.2q59,58 = 0.6
q64,65 = 0.8q65,64 = 0.6q67,58 = 0.2---
Table 6. Loading time [minutes].
Table 6. Loading time [minutes].
t12,3,(1)(2) = 8t13,2,(1)(2) = 10t114,5,(1)(2) = 6t15,4,(1)(2) = 10t17,8,(1)(2) = 6t18,7,(1)(2) = 16
t115,16,(1)(2) = 20t116,15,(1)(2) = 8t121,22,(1)(2) = 8t122,21,(1)(2) = 6t124,25,(1)(2) = 10t114,23,(1)(2) = 8
t119,28,(1)(2) = 10t129,30,(1)(2) = 6t132,31,(1)(2) = 8t126,35,(1)(2) = 8t144,45,(1)(2) = 10t145,44,(1)(2) = 8
t142,41,(1)(2) = 20t155,54,(1)(2) = 8t152,51,(1)(2) = 10t151,50,(1)(2) = 6t149,48,(1)(2) = 20t159,58,(1)(2) = 10
t164,65,(1)(2) = 16t165,64,(1)(2) = 10t167,58,(1)(2) = 6---
Table 7. Distances between intersections km.
Table 7. Distances between intersections km.
d21,2 = 1.0d22,3 = 0.8d23,4 = 1.3d24,5 = 0.7d25,6 = 1.4d26,7 = 0.8
d27,8 = 1.6d28,9 = 2.0d21,10 = 1.2d22,11 = 1.2d23,12 = 1.2d24,13 = 1.2
d25,14 = 1.2d26,15 = 1.2d27,16 = 1.2d28,17 = 1.2d29,18 = 1.2d21,11 = 1.5
d22,10 = 1.5d22,12 = 1.0d23,11 = 1.0d23,13 = 1.8d24,12 = 1.8d24,14 = 1.0
d25,13 = 1.0d25,15 = 1.8d26,14 = 1.8d26,16 = 1.0d27,15 = 1.0d27,17 = 2.0
d28,16 = 2.0d28,18 = 2.3d29,17 = 2.3d210,11 = 1.0d211,12 = 0.8d212,13 = 1.3
d213,14 = 0.7d214,15 = 1.4d215,16 = 0.8d216,17 = 1.6d217,18 = 2.0d210,19 = 1.0
d211,20 = 1.0d212,21 = 1.0d213,22 = 1.0d214,23 = 1.0d215,24 = 1.0d216,25 = 1.0
d217,26 = 1.0d218,27 = 1.0d210,20 = 1.4d211,19 = 1.4d211,21 = 1.3d212,20 = 1.3
d212,22 = 1.8d213,21 = 1.8d213,23 = 1.2d214,22 = 1.2d214,24 = 1.7d215,23 = 1.7
d215,25 = 1.3d216,24 = 1.3d216,26 = 2.3d217,25 = 2.3d217,27 = 2.3d218,26 = 2.3
d219,20 = 1.2d220,21 = 0.8d221,22 = 1.3d222,23 = 0.7d223,24 = 1.4d224,25 = 0.8
d225,26 = 1.6d226,27 = 2.0d219,28 = 1.5d220,29 = 1.5d221,30 = 1.5d222,31 = 1.5
d223,32 = 1.5d224,33 = 1.5d225,34 = 1.5d226,35 = 1.5d227,36 = 1.5d219,29 = 2.1
d220,28 = 2.1d220,30 = 1.7d221,29 = 1.7d221,31 = 2.0d222,30 = 2.0d222,32 = 1.6
d223,31 = 1.6d223,33 = 2.0d224,32 = 2.0d224,34 = 1.7d225,33 = 1.7d225,35 = 2.1
d226,34 = 2.1d226,36 = 2.5d227,35 = 2.5d228,29 = 1.0d229,30 = 0.8d230,31 = 1.3
d231,32 = 0.7d232,33 = 1.4d233,34 = 0.8d234,35 = 1.6d235,36 = 2.0d228,37 = 1.2
d229,38 = 1.2d230,39 = 1.2d231,40 = 1.2d232,41 = 1.2d233,42 = 1.2d234,43 = 1.2
d235,44 = 1.2d236,45 = 1.2d228,38 = 1.6d229,37 = 1.6d229,39 = 1.4d230,38 = 1.4
d230,40 = 1.8d231,39 = 1.8d231,41 = 1.4d232,40 = 1.4d232,42 = 1.8d233,41 = 1.8
d233,44 = 1.4d234,42 = 1.4d234,44 = 2.0d235,43 = 2.0d235,45 = 2.3d236,44 = 2.3
d237,38 = 1.0d238,39 = 0.8d239,40 = 1.3d240,41 = 0.7d241,42 = 1.4d242,43 = 0.8
d243,44 = 1.6d244,45 = 2.0d237,48 = 2.0d238,49 = 2.0d239,50 = 2.0d240,51 = 2.0
d241,52 = 2.0d242,53 = 2.0d243,54 = 2.0d244,55 = 2.0d245,56 = 2.0d237,49 = 2.2
d238,48 = 2.2d238,50 = 2.1d239,49 = 2.1d239,51 = 2.4d240,50 = 2.4d240,52 = 2.1
d241,51 = 2.1d241,53 = 2.4d242,52 = 2.4d242,54 = 2.1d243,53 = 2.1d243,55 = 2.6
d244,54 = 2.6d244,56 = 2.8d245,55 = 2.8d248,49 = 1.0d249,50 = 0.8d250,51 = 1.3
d251,52 = 0.7d252,53 = 1.4d253,54 = 0.8d254,55 = 1.6d255,56 = 2.0d248,57 = 1.4
d249,58 = 1.4d250,59 = 1.4d251,60 = 1.4d252,61 = 1.4d253,62 = 1.4d254,63 = 1.4
d255,64 = 1.4d256,65 = 1.4d248,58 = 1.7d249,57 = 1.7d249,59 = 1.6d250,58 = 1.6
d250,60 = 1.9d251,59 = 1.9d251,61 = 1.6d252,60 = 1.6d252,62 = 2.0d253,61 = 2.0
d253,52 = 1.6d254,62 = 1.6d254,64 = 2.1d255,63 = 2.1d255,65 = 2.4d256,64 = 2.4
d257,58 = 1.0d258,59 = 0.8d259,60 = 1.3d260,61 = 0.7d261,62 = 1.4d262,63 = 0.8
d263,64 = 1.6d264,65 = 2.0d257,66 = 1.0d258,67 = 1.0d259,68 = 1.0d260,69 = 1.0
d261,70 = 1.0d262,71 = 1.0d263,72 = 1.0d264,73 = 1.0d265,74 = 1.0d257,67 = 1.4
d258,66 = 1.4d258,68 = 1.3d259,67 = 1.3d259,69 = 1.6d260,68 = 1.6d260,70 = 1.2
d261,69 = 1.2d261,71 = 1.7d262,70 = 1.7d262,72 = 1.3d263,71 = 1.3d263,73 = 1.9
d264,72 = 1.9d264,74 = 2.2d265,73 = 2.2d266,67 = 1.0d267,68 = 0.8d268,69 = 1.3
d269,70 = 0.7d270,71 = 1.4d271,72 = 0.8d272,73 = 1.6d273,74 = 2.0-
Table 8. Test settings of ant colony algorithm parameters.
Table 8. Test settings of ant colony algorithm parameters.
Test α β ρ Test α β ρ Test α β ρ
110.50.221110.241150.2
210.50.422110.442150.4
310.50.623110.643150.6
410.50.824110.844150.8
530.50.225310.245350.2
630.50.426310.446350.4
730.50.627310.647350.6
830.50.828310.848350.8
950.50.229510.249550.2
1050.50.430510.450550.4
1150.50.631510.651550.6
1250.50.832510.852550.8
13100.50.2331010.2531050.2
14100.50.4341010.4541050.4
15100.50.6351010.6551050.6
16100.50.8361010.8561050.8
17200.50.2372010.2572050.2
18200.50.4382010.4582050.4
19200.50.6392010.6592050.6
20200.50.8402010.8602050.8
Table 9. Summary of the results of the ant colony algorithm kWh.
Table 9. Summary of the results of the ant colony algorithm kWh.
TestEnergy ExpenditureTestEnergy ExpenditureTestEnergy ExpenditureTestEnergy ExpenditureTestEnergy Expenditure
123213290252613732049282
224114293262593831350273
324515288272633934051291
423416295282664033552274
524817323292714124553291
625618333302824222354289
726419343312774324855302
825520355322804425056313
926721223332884526257320
1026522212342904627158333
1126123216353104726159312
1227324211363114826760315
Table 10. Results of the ant colony algorithm [kWh].
Table 10. Results of the ant colony algorithm [kWh].
Energy ExpenditureEnergy ExpenditureEnergy ExpenditureEnergy ExpenditureEnergy ExpenditureEnergy Expenditure
243220211221237234
221232236223211218
211227332227233220
233211218230231235
218214235228222223
Table 11. Test input parameters of the genetic algorithm [kWh].
Table 11. Test input parameters of the genetic algorithm [kWh].
TestpcrossEnergy ExpenditureTestpcrossEnergy ExpenditureTestpcrossEnergy Expenditure
10.221120.423130.6331
40.819951225---
Table 12. Exemplary results of the genetic algorithm [kWh].
Table 12. Exemplary results of the genetic algorithm [kWh].
Energy ExpenditureEnergy ExpenditureEnergy ExpenditureEnergy ExpenditureEnergy ExpenditureEnergy Expenditure
220231220201216241
215220224213220223
199223215221212215
229234223203230221
218218237206232218
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Izdebski, M.; Jacyna, M. An Efficient Hybrid Algorithm for Energy Expenditure Estimation for Electric Vehicles in Urban Service Enterprises. Energies 2021, 14, 2004. https://doi.org/10.3390/en14072004

AMA Style

Izdebski M, Jacyna M. An Efficient Hybrid Algorithm for Energy Expenditure Estimation for Electric Vehicles in Urban Service Enterprises. Energies. 2021; 14(7):2004. https://doi.org/10.3390/en14072004

Chicago/Turabian Style

Izdebski, Mariusz, and Marianna Jacyna. 2021. "An Efficient Hybrid Algorithm for Energy Expenditure Estimation for Electric Vehicles in Urban Service Enterprises" Energies 14, no. 7: 2004. https://doi.org/10.3390/en14072004

APA Style

Izdebski, M., & Jacyna, M. (2021). An Efficient Hybrid Algorithm for Energy Expenditure Estimation for Electric Vehicles in Urban Service Enterprises. Energies, 14(7), 2004. https://doi.org/10.3390/en14072004

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop