An Efficient Hybrid Algorithm for Energy Expenditure Estimation for Electric Vehicles in Urban Service Enterprises
Abstract
:1. Introduction
2. Research Status in the Literature
2.1. Electric Vehicles—Development, Benefits, and Weaknesses
2.2. Routing Problem Including Electric Vehicles
3. A Decision for Estimating Energy Expenditure Using Electric Vehicles
3.1. General Assumptions for the Model
- 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
3.3. Quantities Sought
3.4. Restrictions
3.5. Energy Expenditure as a Criterion
- drag force:
- rolling resistance force:
4. An Efficient Hybrid Algorithm for Solving the Problem
4.1. Preliminary Assumptions
4.2. Operation of the Ant Colony Algorithm
- 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.
- —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,
- —heuristic information:
- —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,
- —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.
- —another ant in the anthill ,
- —pheromone volatilization rate (0 < ρ 1),
- —pheromone reinforcement, for the first iteration takes a value at each connection equal to .
- —total route traveled in -iteration by -th ant performing all tasks, the route is the basis for the calculation of energy expenditure, criterion function (8);
- —penalty for exceeding the admissible working time by an individual vehicle (model constraint (2)) in the route created by -th ant in -th iteration, the algorithm assumes that the penalty for exceeding is half of the accumulated pheromone on the route;
- 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.
4.3. Operation of the Genetic Algorithm
- 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.
- 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):
- —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;
- —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;
- 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.
- 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
5. Verification of Algorithms at the Example of an Urban Utility Company
5.1. Model Input Data
5.2. Discussion of the Results for the Sensitivity Analysis of the Ant Colony Algorithm
5.3. Discussion of Results for the Sensitivity Analysis of the Genetic Algorithm
6. Conclusions
Author Contributions
Funding
Acknowledgments
Conflicts of Interest
References
- What Is Meant by the Term “Sustainability”? Available online: http://www.fao.org/3/ai388e/AI388E05.htm (accessed on 16 January 2021).
- UN Secretary-General. World Commission on Environment and Development: Our Common Future; Oxford University Press: Oxford, UK, 1987. [Google Scholar]
- 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]
- Westbrook, M.H.; Westbrook, M. The Electric Car: Development and Future of Battery, Hybrid and Fuel-Cell Cars; IET: London, UK, 2001. [Google Scholar]
- 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]
- 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]
- 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]
- 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]
- Erickson, L.; Ma, S. Solar-Powered Charging Networks for Electric Vehicles. Energies 2021, 14, 966. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- Izdebski, M.; Jacyna, M. The organization of the municipal waste collection: The decision model. Rocz. Ochr. Sr. 2018, 20, 919–933. [Google Scholar]
- 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]
- Jacyna, M.; Kotylak, P. Decision-making problems of collective transport development in terms of sustainable urban mobility. J. KONBIN 2020, 50. [Google Scholar] [CrossRef]
- 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]
- Kühne, R. Electric Buses–An Energy Efficient Urban Transportation Means. Energy 2010, 35, 4510–4513. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Nowak, M.; Ergun, O.; White, C.C., III. Pickup and delivery with split loads. Transp. Sci. 2008, 42, 32–43. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- United Nations. Transforming Our World: The 2030 Agenda for Sustainable Development; Division for Sustainable Development Goals; UN: New York, NY, USA, 2015. [Google Scholar]
- 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]
- 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]
- 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).
- Niestadt, M.; Bjornavold, A. Electric Road Vehicles in the European Union; European Parliament: Brussels, Belgium, 2019. [Google Scholar]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Laporte, G. Fifty Years of Vehicle Routing. Transp. Sci. 2009, 43, 408–416. [Google Scholar] [CrossRef]
- Mancini, S. The hybrid vehicle routing problem. Transp. Res. Part C 2017, 78, 1–12. [Google Scholar] [CrossRef]
- 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]
- 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]
- Valouxis, C.; Housos, E. Combined bus and driver scheduling. Comput. Oper. Res. 2002, 29, 243–259. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- Burkhard, R.; Dell’Amico, M.; Marttelo, S. Assignment Problems; Society for Industrial and Applied Mathematics; SIAM: Philadelphia, PA, USA, 2009. [Google Scholar]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Goeke, D.; Schneider, M. Routing a mixed fleet of electric and conventional vehicles. Eur. J. Oper. Res. 2015, 245, 81–99. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Dorigo, M.; Stutzle, T. Ant Colony Optimization; Bradford Books: Cambridge, MA, USA, 2004. [Google Scholar]
- Dorigo, M.; Blum, C. Ant Colony Optimization Theory: A survey. Theor. Comput. Sci. 2005, 344, 243–278. [Google Scholar] [CrossRef]
- Goldberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning, 1st ed.; Addison-Wesley Professional: Boston, MA, USA, 1989. [Google Scholar]
- Arora, S.; Barak, B. Computational Complexity: A Modern Approach; Cambridge University Press: New York, NY, USA, 2009. [Google Scholar]
- Oelschlägel, T.; Knust, S. Solution approaches for storage loading problems with stacking constraints. Comput. Oper. Res. 2021, 127, 105–142. [Google Scholar] [CrossRef]
- 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]
- Early, D.; Schellekens, M. Running time of the Treapsort algorithm. Theor. Comput. Sci. 2013, 487, 65–73. [Google Scholar] [CrossRef]
- 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]
Symbol | Meaning |
---|---|
S | The set of collective starting points of a vehicle route, e.g., cargo distribution points, dumps where s is an element of the S set. |
I1 | The 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. |
I2 | The 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. |
B | A set of electric vehicle battery charging stations, where b is an element of the B set. |
VEH | The set of vehicles in the service enterprise, where veh is an element of the VEH set. |
TYPE | The set of vehicle types in a service company that differ in cargo capacity, where type is an element of the TYPE set. |
TASK | A set of transportation tasks to be completed on a given working day, where task is an element of the TASK set. |
P | The set of periods with interpretations of the morning, midday, and afternoon peaks, where p is an element of the P set. |
vveh,typ,p | Speed of a vehicle of a given type during a given period of a working day, where V is the speed matrix. |
ntyp | The 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,typ | The 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,typ | The load capacity of a vehicle of a given type, where C is the vehicle load capacity matrix. |
t1i2′,i2″,typ | The 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,typ | Unloading or loading time at a collective point dependent on vehicle type, where T2 is a matrix of unloading or loading times. |
tch,veh,typ,b | The 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,i1 | The 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,s | The 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,typ | The allowable travel distance for each electric vehicle, where Dpem is the matrix of allowable travel distances for each vehicle. |
m1s,i1,veh | The 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″,veh | The 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,veh | The 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. |
g | Earth 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. |
cx | The drag coefficient, where cx = 0.2. |
A | Vehicle frontal area, for each type of vehicle an identical frontal area is assumed, A = 20 m2. |
Symbol | Meaning |
---|---|
xveh,typ,task | The 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,p | A 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,p | A 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,p | A 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,task | Making 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. |
Algorithms | Input Data | Loops | Estimation | Complexity |
---|---|---|---|---|
AA | MR, N, G | 2∙G2∙MR∙N + G∙MR∙N | 2∙G2∙MR∙N + G∙MR∙N ≤ (MR + N + G)4 | O((MR + N + G)4) polynomial |
GA | MR, N, G | 3∙G∙MR∙N | 3∙G∙MR∙N ≤ (MR + N + G)3 | O((MR + N + G)3) polynomial |
AA | Loops | GA | Loops |
---|---|---|---|
Determination of the probability | G2∙MR∙N | Selection | G∙MR∙N |
Determination of the distribution | G2∙MR∙N | Crossover | G∙MR∙N |
A cyclic update | G∙MR∙N | Population assessment | G∙MR∙N |
q2,3 = 0.4 | q3,2 = 0.6 | q4,5 = 0.2 | q5,4 = 0.6 | q7,8 = 0.2 | q8,7 = 0.8 |
q15,16 = 1.2 | q16,15 = 0.4 | q21,22 = 0.4 | q22,21 = 0.2 | q24,25 = 0.6 | q14,23 = 0.4 |
q19,28 = 0.6 | q29,30 = 0.2 | q32,31 = 0.4 | q26,35 = 0.4 | q44,45 = 0.6 | q45,44 = 0.2 |
q42,41 = 1.2 | q55,54 = 0.4 | q52,51 = 0.6 | q51,50 = 0.2 | q49,48 = 1.2 | q59,58 = 0.6 |
q64,65 = 0.8 | q65,64 = 0.6 | q67,58 = 0.2 | - | - | - |
t12,3,(1)(2) = 8 | t13,2,(1)(2) = 10 | t114,5,(1)(2) = 6 | t15,4,(1)(2) = 10 | t17,8,(1)(2) = 6 | t18,7,(1)(2) = 16 |
t115,16,(1)(2) = 20 | t116,15,(1)(2) = 8 | t121,22,(1)(2) = 8 | t122,21,(1)(2) = 6 | t124,25,(1)(2) = 10 | t114,23,(1)(2) = 8 |
t119,28,(1)(2) = 10 | t129,30,(1)(2) = 6 | t132,31,(1)(2) = 8 | t126,35,(1)(2) = 8 | t144,45,(1)(2) = 10 | t145,44,(1)(2) = 8 |
t142,41,(1)(2) = 20 | t155,54,(1)(2) = 8 | t152,51,(1)(2) = 10 | t151,50,(1)(2) = 6 | t149,48,(1)(2) = 20 | t159,58,(1)(2) = 10 |
t164,65,(1)(2) = 16 | t165,64,(1)(2) = 10 | t167,58,(1)(2) = 6 | - | - | - |
d21,2 = 1.0 | d22,3 = 0.8 | d23,4 = 1.3 | d24,5 = 0.7 | d25,6 = 1.4 | d26,7 = 0.8 |
d27,8 = 1.6 | d28,9 = 2.0 | d21,10 = 1.2 | d22,11 = 1.2 | d23,12 = 1.2 | d24,13 = 1.2 |
d25,14 = 1.2 | d26,15 = 1.2 | d27,16 = 1.2 | d28,17 = 1.2 | d29,18 = 1.2 | d21,11 = 1.5 |
d22,10 = 1.5 | d22,12 = 1.0 | d23,11 = 1.0 | d23,13 = 1.8 | d24,12 = 1.8 | d24,14 = 1.0 |
d25,13 = 1.0 | d25,15 = 1.8 | d26,14 = 1.8 | d26,16 = 1.0 | d27,15 = 1.0 | d27,17 = 2.0 |
d28,16 = 2.0 | d28,18 = 2.3 | d29,17 = 2.3 | d210,11 = 1.0 | d211,12 = 0.8 | d212,13 = 1.3 |
d213,14 = 0.7 | d214,15 = 1.4 | d215,16 = 0.8 | d216,17 = 1.6 | d217,18 = 2.0 | d210,19 = 1.0 |
d211,20 = 1.0 | d212,21 = 1.0 | d213,22 = 1.0 | d214,23 = 1.0 | d215,24 = 1.0 | d216,25 = 1.0 |
d217,26 = 1.0 | d218,27 = 1.0 | d210,20 = 1.4 | d211,19 = 1.4 | d211,21 = 1.3 | d212,20 = 1.3 |
d212,22 = 1.8 | d213,21 = 1.8 | d213,23 = 1.2 | d214,22 = 1.2 | d214,24 = 1.7 | d215,23 = 1.7 |
d215,25 = 1.3 | d216,24 = 1.3 | d216,26 = 2.3 | d217,25 = 2.3 | d217,27 = 2.3 | d218,26 = 2.3 |
d219,20 = 1.2 | d220,21 = 0.8 | d221,22 = 1.3 | d222,23 = 0.7 | d223,24 = 1.4 | d224,25 = 0.8 |
d225,26 = 1.6 | d226,27 = 2.0 | d219,28 = 1.5 | d220,29 = 1.5 | d221,30 = 1.5 | d222,31 = 1.5 |
d223,32 = 1.5 | d224,33 = 1.5 | d225,34 = 1.5 | d226,35 = 1.5 | d227,36 = 1.5 | d219,29 = 2.1 |
d220,28 = 2.1 | d220,30 = 1.7 | d221,29 = 1.7 | d221,31 = 2.0 | d222,30 = 2.0 | d222,32 = 1.6 |
d223,31 = 1.6 | d223,33 = 2.0 | d224,32 = 2.0 | d224,34 = 1.7 | d225,33 = 1.7 | d225,35 = 2.1 |
d226,34 = 2.1 | d226,36 = 2.5 | d227,35 = 2.5 | d228,29 = 1.0 | d229,30 = 0.8 | d230,31 = 1.3 |
d231,32 = 0.7 | d232,33 = 1.4 | d233,34 = 0.8 | d234,35 = 1.6 | d235,36 = 2.0 | d228,37 = 1.2 |
d229,38 = 1.2 | d230,39 = 1.2 | d231,40 = 1.2 | d232,41 = 1.2 | d233,42 = 1.2 | d234,43 = 1.2 |
d235,44 = 1.2 | d236,45 = 1.2 | d228,38 = 1.6 | d229,37 = 1.6 | d229,39 = 1.4 | d230,38 = 1.4 |
d230,40 = 1.8 | d231,39 = 1.8 | d231,41 = 1.4 | d232,40 = 1.4 | d232,42 = 1.8 | d233,41 = 1.8 |
d233,44 = 1.4 | d234,42 = 1.4 | d234,44 = 2.0 | d235,43 = 2.0 | d235,45 = 2.3 | d236,44 = 2.3 |
d237,38 = 1.0 | d238,39 = 0.8 | d239,40 = 1.3 | d240,41 = 0.7 | d241,42 = 1.4 | d242,43 = 0.8 |
d243,44 = 1.6 | d244,45 = 2.0 | d237,48 = 2.0 | d238,49 = 2.0 | d239,50 = 2.0 | d240,51 = 2.0 |
d241,52 = 2.0 | d242,53 = 2.0 | d243,54 = 2.0 | d244,55 = 2.0 | d245,56 = 2.0 | d237,49 = 2.2 |
d238,48 = 2.2 | d238,50 = 2.1 | d239,49 = 2.1 | d239,51 = 2.4 | d240,50 = 2.4 | d240,52 = 2.1 |
d241,51 = 2.1 | d241,53 = 2.4 | d242,52 = 2.4 | d242,54 = 2.1 | d243,53 = 2.1 | d243,55 = 2.6 |
d244,54 = 2.6 | d244,56 = 2.8 | d245,55 = 2.8 | d248,49 = 1.0 | d249,50 = 0.8 | d250,51 = 1.3 |
d251,52 = 0.7 | d252,53 = 1.4 | d253,54 = 0.8 | d254,55 = 1.6 | d255,56 = 2.0 | d248,57 = 1.4 |
d249,58 = 1.4 | d250,59 = 1.4 | d251,60 = 1.4 | d252,61 = 1.4 | d253,62 = 1.4 | d254,63 = 1.4 |
d255,64 = 1.4 | d256,65 = 1.4 | d248,58 = 1.7 | d249,57 = 1.7 | d249,59 = 1.6 | d250,58 = 1.6 |
d250,60 = 1.9 | d251,59 = 1.9 | d251,61 = 1.6 | d252,60 = 1.6 | d252,62 = 2.0 | d253,61 = 2.0 |
d253,52 = 1.6 | d254,62 = 1.6 | d254,64 = 2.1 | d255,63 = 2.1 | d255,65 = 2.4 | d256,64 = 2.4 |
d257,58 = 1.0 | d258,59 = 0.8 | d259,60 = 1.3 | d260,61 = 0.7 | d261,62 = 1.4 | d262,63 = 0.8 |
d263,64 = 1.6 | d264,65 = 2.0 | d257,66 = 1.0 | d258,67 = 1.0 | d259,68 = 1.0 | d260,69 = 1.0 |
d261,70 = 1.0 | d262,71 = 1.0 | d263,72 = 1.0 | d264,73 = 1.0 | d265,74 = 1.0 | d257,67 = 1.4 |
d258,66 = 1.4 | d258,68 = 1.3 | d259,67 = 1.3 | d259,69 = 1.6 | d260,68 = 1.6 | d260,70 = 1.2 |
d261,69 = 1.2 | d261,71 = 1.7 | d262,70 = 1.7 | d262,72 = 1.3 | d263,71 = 1.3 | d263,73 = 1.9 |
d264,72 = 1.9 | d264,74 = 2.2 | d265,73 = 2.2 | d266,67 = 1.0 | d267,68 = 0.8 | d268,69 = 1.3 |
d269,70 = 0.7 | d270,71 = 1.4 | d271,72 = 0.8 | d272,73 = 1.6 | d273,74 = 2.0 | - |
Test | Test | Test | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
1 | 1 | 0.5 | 0.2 | 21 | 1 | 1 | 0.2 | 41 | 1 | 5 | 0.2 |
2 | 1 | 0.5 | 0.4 | 22 | 1 | 1 | 0.4 | 42 | 1 | 5 | 0.4 |
3 | 1 | 0.5 | 0.6 | 23 | 1 | 1 | 0.6 | 43 | 1 | 5 | 0.6 |
4 | 1 | 0.5 | 0.8 | 24 | 1 | 1 | 0.8 | 44 | 1 | 5 | 0.8 |
5 | 3 | 0.5 | 0.2 | 25 | 3 | 1 | 0.2 | 45 | 3 | 5 | 0.2 |
6 | 3 | 0.5 | 0.4 | 26 | 3 | 1 | 0.4 | 46 | 3 | 5 | 0.4 |
7 | 3 | 0.5 | 0.6 | 27 | 3 | 1 | 0.6 | 47 | 3 | 5 | 0.6 |
8 | 3 | 0.5 | 0.8 | 28 | 3 | 1 | 0.8 | 48 | 3 | 5 | 0.8 |
9 | 5 | 0.5 | 0.2 | 29 | 5 | 1 | 0.2 | 49 | 5 | 5 | 0.2 |
10 | 5 | 0.5 | 0.4 | 30 | 5 | 1 | 0.4 | 50 | 5 | 5 | 0.4 |
11 | 5 | 0.5 | 0.6 | 31 | 5 | 1 | 0.6 | 51 | 5 | 5 | 0.6 |
12 | 5 | 0.5 | 0.8 | 32 | 5 | 1 | 0.8 | 52 | 5 | 5 | 0.8 |
13 | 10 | 0.5 | 0.2 | 33 | 10 | 1 | 0.2 | 53 | 10 | 5 | 0.2 |
14 | 10 | 0.5 | 0.4 | 34 | 10 | 1 | 0.4 | 54 | 10 | 5 | 0.4 |
15 | 10 | 0.5 | 0.6 | 35 | 10 | 1 | 0.6 | 55 | 10 | 5 | 0.6 |
16 | 10 | 0.5 | 0.8 | 36 | 10 | 1 | 0.8 | 56 | 10 | 5 | 0.8 |
17 | 20 | 0.5 | 0.2 | 37 | 20 | 1 | 0.2 | 57 | 20 | 5 | 0.2 |
18 | 20 | 0.5 | 0.4 | 38 | 20 | 1 | 0.4 | 58 | 20 | 5 | 0.4 |
19 | 20 | 0.5 | 0.6 | 39 | 20 | 1 | 0.6 | 59 | 20 | 5 | 0.6 |
20 | 20 | 0.5 | 0.8 | 40 | 20 | 1 | 0.8 | 60 | 20 | 5 | 0.8 |
Test | Energy Expenditure | Test | Energy Expenditure | Test | Energy Expenditure | Test | Energy Expenditure | Test | Energy Expenditure |
---|---|---|---|---|---|---|---|---|---|
1 | 232 | 13 | 290 | 25 | 261 | 37 | 320 | 49 | 282 |
2 | 241 | 14 | 293 | 26 | 259 | 38 | 313 | 50 | 273 |
3 | 245 | 15 | 288 | 27 | 263 | 39 | 340 | 51 | 291 |
4 | 234 | 16 | 295 | 28 | 266 | 40 | 335 | 52 | 274 |
5 | 248 | 17 | 323 | 29 | 271 | 41 | 245 | 53 | 291 |
6 | 256 | 18 | 333 | 30 | 282 | 42 | 223 | 54 | 289 |
7 | 264 | 19 | 343 | 31 | 277 | 43 | 248 | 55 | 302 |
8 | 255 | 20 | 355 | 32 | 280 | 44 | 250 | 56 | 313 |
9 | 267 | 21 | 223 | 33 | 288 | 45 | 262 | 57 | 320 |
10 | 265 | 22 | 212 | 34 | 290 | 46 | 271 | 58 | 333 |
11 | 261 | 23 | 216 | 35 | 310 | 47 | 261 | 59 | 312 |
12 | 273 | 24 | 211 | 36 | 311 | 48 | 267 | 60 | 315 |
Energy Expenditure | Energy Expenditure | Energy Expenditure | Energy Expenditure | Energy Expenditure | Energy Expenditure |
---|---|---|---|---|---|
243 | 220 | 211 | 221 | 237 | 234 |
221 | 232 | 236 | 223 | 211 | 218 |
211 | 227 | 332 | 227 | 233 | 220 |
233 | 211 | 218 | 230 | 231 | 235 |
218 | 214 | 235 | 228 | 222 | 223 |
Test | pcross | Energy Expenditure | Test | pcross | Energy Expenditure | Test | pcross | Energy Expenditure |
---|---|---|---|---|---|---|---|---|
1 | 0.2 | 211 | 2 | 0.4 | 231 | 3 | 0.6 | 331 |
4 | 0.8 | 199 | 5 | 1 | 225 | - | - | - |
Energy Expenditure | Energy Expenditure | Energy Expenditure | Energy Expenditure | Energy Expenditure | Energy Expenditure |
---|---|---|---|---|---|
220 | 231 | 220 | 201 | 216 | 241 |
215 | 220 | 224 | 213 | 220 | 223 |
199 | 223 | 215 | 221 | 212 | 215 |
229 | 234 | 223 | 203 | 230 | 221 |
218 | 218 | 237 | 206 | 232 | 218 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
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
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 StyleIzdebski, 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 StyleIzdebski, 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