Next Article in Journal
Analysis of the Activities That Make Up the Reverse Logistics Processes and Their Importance for the Future of Logistics Networks: An Exploratory Study Using the TOPSIS Technique
Previous Article in Journal
Industrial Packaging Performance Indicator Using a Group Multicriteria Approach: An Automaker Reverse Operations Case
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Green Logistics Technique for Planning Merchandise Deliveries: A Case Study

by
Carlos Hernández-Mejía
1,*,
Delia Torres-Muñoz
2,
Everardo Inzunza-González
3,
Carlos Sánchez-López
4 and
Enrique Efrén García-Guerrero
3
1
Subdirección de Posgrados e Investigación, Instituto Tecnológico Superior de Misantla, Km. 1.8 Carretera a Loma del Cojolite, Misantla 93821, Mexico
2
Instituto Tecnológico Superior de San Martín Texmelucan, Camino a Barranca de Pesos S/N, San Martín Texmelucan 74120, Mexico
3
Facultad de Ingeniería, Arquitectura y Diseño, Universidad Autónoma de Baja California, Carretera Transpeninsular Ensenada-Tijuana No. 3917, Ensenada 22860, Mexico
4
Department of Electronics, Autonomous University of Tlaxcala, Av. Universidad No. 1, Tlaxcala 90070, Mexico
*
Author to whom correspondence should be addressed.
Logistics 2022, 6(3), 59; https://doi.org/10.3390/logistics6030059
Submission received: 24 July 2022 / Revised: 4 August 2022 / Accepted: 9 August 2022 / Published: 15 August 2022

Abstract

:
Background: The health crisis due to COVID-19 has changed the habits of social coexistence and therefore has had a significant impact on several economic sectors, including logistics. Nowadays, this sector faces one of the most difficult challenges in history and hence has carried out innovative strategies to adjust to the new normal and guarantee the permanence of the supply chain. In this paper, a novel green logistics technique for planning merchandise deliveries is introduced herein. Methods: This technique is based on a modified version of the Resistive Grid Path Planning Methodology (RGPPM) to establish a path through georeferential locations for delivering merchandise to customers. To do this, multi-connected resistive grids are based at customer locations for searching the path connecting the positions through electronic circuit analysis techniques. Results: Experimental results reveal that the proposed method can find a minimum spanning tree that connects all hand-over points by a continuous path. This finding has represented a reduction of around 45% in the length of the path with respect to the longest path. Conclusions: Finally, this technique is capable of addressing different optimization strategies, locating a merchandise distribution center and exploring environmental standards to reduce fuel consumption.

1. Introduction

At the end of 2019, the outbreak of coronavirus disease (COVID-19) emerged in China. SARS-COV-2 virus spreading has significantly disrupted global production chains causing serious disturbances in the supply chain [1]. According to International Labor Organization (ILO), the health security measures such as social distancing and the stoppage of non-essential activities have affected around 2.7 billion workers; i.e., 81% of the world’s workforce. Therefore, it is estimated that 1.25 billion workers in the world remain at risk of losing their job because they belong to the sectors where the impact of COVID-19 has been most significant, particularly in accommodation services, the manufacturing industry, commerce, and food services.
Health contingency has devastated the world of work causing immense human suffering and exposing the extreme vulnerability of millions of people and companies. According to the report [2] of an Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES) workshop on diversity and pandemics, the pandemics in the near future will emerge more frequently, and they will have a greater impact on the global economy. Therefore, it is important that the most vulnerable companies such as: micro, small, and medium-sized companies explore strategies and technological techniques that help them to go through these pandemic events. This paper explores an emerging technique in order to deliver merchandise during confinement events.
Food deliveries have increased dramatically in many countries since 2020 due to mobility restrictions and restaurant closures as a result of the coronavirus pandemic event [3]. Companies in the online food distribution industry estimate that many of their first-time online shoppers have experienced the benefits and intend not to revert to their pre-crisis habits [4]. This opens up the possibility for exponential growth in activities related to the distribution and transportation of meals and food. However, the movement of vehicles to customers in order to deliver merchandise generates a significant amount of CO2 emissions that affect the quality of life of citizens and the climate [5,6]. Therefore, it is important to establish a green strategy for the logistics and distribution of merchandise [7].
Logistics [8,9] is a fundamental stage in the global economy and therefore it must be strategically transformed with regard to carbon emissions. Green logistics [10,11] explore the measurement and minimization of environmental impact during logistics activities.
Green strategies [12,13] are a series of actions aimed at the interrelation between social awareness and business problems. Companies with large volumes of transport and distribution pursue green climate protection strategies to avoid CO2 emissions. On the other hand, for the transport industry, green strategies mean not only a reduction in CO2 emissions but also significant economic savings due to optimization methods for planning the distribution of merchandise, for example, path planning is a critical point for optimizing the transfer of merchandise because, depending on how well the green strategy is carried out, it it is possible to reduce from 10% to 30% of the distance in kilometers and, consequently, the carbon emissions into the atmosphere and fuel consumption, which represents social and economic benefits. This paper proposes a green strategy focused on path planning for delivering merchandise with the purpose of reducing CO2 emissions and therefore cause social and economic benefits.
Lecheria y Quesos Artesanales “Los Ángeles” is a Mexican startup that markets products derived from milk. Since the beginning of 2020, the enterprise has maintained the distribution of merchandise in four principal cities: Puebla, Xalapa, Veracruz, and Boca del Rio. As a result of the health emergency and the restrictions for opening a commercial space in order to serve customers, the startup has been forced to migrate towards a strategy for delivering merchandise at private homes. This paper exposes the experimental results of having carried out a new green technique for planning merchandise deliveries using the database of customers coming from the business Lecheria y Quesos Artesanales “Los Ángeles”.
Due to the increase in e-commerce and the high demand in the delivery of merchandise to the particular locations of customers, the vast majority of businesses and companies need to carry out innovative strategies to remove problems in last mile logistics [14]. According to the literature [15], there is a famous problem that has been explored in practical real-life situations for the purpose of minimizing last mile costs, it is known as the traveling salesman problem. This problem has been considered a fundamental piece for vehicle routing in the logistics area [16]. Basically, it is possible to relate this problem to the need for customers to be visited with the intention of delivering merchandise or packages starting from a base position on the map. This paper explores a novel technique for delivering merchandise in the last mile logistics.
Although the best technique for this type of problem is unknown, the following research proposes an original strategy related to the behavior of the physical phenomenon of electric current, from the perspective of circuit theory, in order to deliver merchandise in the last mile. Moreover, the exploration of this approach opens the possibility of laying the foundations to face computational problems as the number of customers increases.
The manuscript is organized as follows: Section 2 introduces the proposed RGPPM algorithm for planning merchandise deliveries. In Section 3, the experimental results are presented. Section 4 presents a brief discussion about the previous work and the experimental results. Finally, in Section 5, some conclusions are drawn.

2. RGPPM for Delivering Merchandise

RGPPM is a novel emerging alternative for planning trajectories from a start point to an end point. This methodology has been successfully evaluated for the movement of terrestrial robotic systems [17,18]. However, until now, configuration spaces with regular geometrical patterns have been exclusively explored, as shown in Figure 1.
In order to expand RGPPM for planning merchandise deliveries, it is necessary to establish configuration spaces with multi-connected irregular geometrical patterns, as shown in Figure 2. This configuration represents a novel contribution in research with the purpose of merging path planning and merchandise deliveries using geometrical patterns.
As it can be seen, the nodes have been fully connected to each other and therefore the edges represent all possible connections between nodes. If it is assumed that each node represents a position to deliver merchandise, then the edges describe the paths between the different destinations. Subsequently, the modeling of the configuration space using resistive grids is carried out, as shown in Figure 3. In the same way that it has been shown previously, this configuration represents a novel contribution in research with the purpose of merging path planning and merchandise deliveries using resistive grids.
According to the irregular geometrical configuration of Figure 2, the values of the edges are related to the length between the nodes and therefore the values of the resistors in Figure 3, R 1 , R 2 , R 3 , R 4 , R 5 , R 6 , R 7 , R 8 , R 9 , and R 10 are related to the length of the paths between locations of handing over 1 , 2 , 3 , 4 , 5 .
With the purpose of carrying out a simple example, the resistive grid in Figure 3 has been characterized by the resistor values shown in Table 1. These values have been calculated by locating the nodes on a Cartesian plane and measuring the distance between each one.
After establishing values for resistors in Figure 3, it is necessary to connect power supplies on all nodes in the resistive grid. In order to do this, E 1 , E 2 , E 3 , E 4 , and E 5 power supplies have been connected to nodes 1 , 2 , 3 , 4 , 5 , respectively. It is important to determine which is the node of the starting delivery for assigning to the power supply located in that position a value different from 0 V. In this case, E 1 = 10 V, and from E 2 to E 5 , it is equal to 0 V. It is important to mention that the above numerical value of E 1 has been established for the purpose of carrying out a simple example since all positive values are considered feasible.
When the characterization of the resistive grid has been achieved, the Modified Nodal Analysis (MNA) [19] procedure is carried out in order to formulate the equilibrium equations. Therefore, conductance and independent voltage source stamps [20,21] are carried out as follows:
v i v j RHS i G G 0 j G G 0
v i v j i E RHS i 0 0 + 1 0 j 0 0 1 0 BE + 1 1 0 E
where v i and v j represent the nodal voltage with respect to ground, i and j are the nodes of the elements, E is the supplied voltage by the independent voltage source (IVS), G = 1 R is the equivalent conductance, RHS means the right-hand side of the equilibrium equation, and i E is the unknown current that crosses by the IVS.
According to Equations (1) and (2), the formulation of the linear system of equations in Figure 3 is based on the form A x = b as follows:
A = 1.04 0.17 0.13 0.44 0.28 1 0 0 0 0 0.17 0.65 0.12 0.13 0.22 0 1 0 0 0 0.13 0.12 0.66 0.16 0.23 0 0 1 0 0 0.44 0.13 0.16 1.03 0.28 0 0 0 1 0 0.28 0.22 0.23 0.28 1.02 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0
x = v 1 v 2 v 3 v 4 v 5 i E 1 i E 2 i E 3 i E 4 i E 5 b = 0 0 0 0 0 10 0 0 0 0
where A is a sparse matrix (MNA matrix), and x and b are vectors (unknowns and stimulus vector, respectively). It is important to mention that RGPPM uses a direct method for solving systems of linear equations called a multifrontal method from Unsymmetric MultiFrontal PACK (UMFPACK) [22].
After carrying out the multifrontal method in the linear system of equations, the solutions are shown in Table 2.
With the values of the nodal voltages, the next step consists of finding the path from the starting delivery to the new position of the next delivery. In order to do this, RGPPM uses a Local Current Comparison (LCC) [23] algorithm for calculating the branch currents and therefore it establishes the largest branch current.
According to Figure 4, the LCC algorithm starts at node one. Firstly, it determines the neighboring nodes to the actual node 2 , 3 , 4 , 5 . Then, it calculates the branch currents to neighboring nodes using the following expression:
I j = V V n R j
where R j , V n , V, and I j represent the resistance value in the j-branch, the voltage of the neighboring node in the j-branch, the actual nodal voltage, and the current in the j-branch, respectively.
Table 3 shows the currents of the branches connected to node one. According to this list of values, the highest current flows through resistor R 2 and therefore this branch must be considered as the path to the new actual node 4 .
After that, the previous node and the branches connected to it should be removed from the multi-connected resistive grid configuration before continuing with the new actual node. However, the previous node and the branch of the highest current connected to it must be stored because these elements belong to the complete path for merchandise deliveries.
Figure 5 shows the modified multi-connected resistive grid configuration after removing the previous node and the branches connected to it. Finally, with this new multi-connected resistive grid configuration, the stages of assigning power supplies to the nodes, the formulation of the equilibrium equations, the solution of the system of equations, the path search procedure and the modification of the resistive grid configuration must be carried out again and repeatedly until all positions for delivering merchandise have been passed.
As can be seen in Figure 6, the resistors in red R 2 , R 3 , R 4 , R 9 form a connection tree with all the nodes of the resistive grid and therefore it is possible to establish the complete path for delivering merchandise. In addition, the sequence of nodes in the complete path 1 , 4 , 5 , 3 , 2 represents the sequence of positions in order to deliver merchandise. In addition, the sum of the resistor values represents the total length of the path for delivering merchandise; for this example, it is 18.12.
It is important to mention that, for this example, the starting point to deliver merchandise has been established at node one and therefore the route sequence is 1 , 4 , 5 , 3 , 2 . However, the technique is capable of modifying the starting node to a different point and thus it generates a new tree connection with all nodes.
Figure 7 shows two alternate full path resistive grid configurations using different starting nodes. In the case of the resistive grid on the left side, the starting node is the number three which generates a connection tree formed by the resistors in red R 4 , R 3 , R 2 , R 1 and a path sequence 3 , 5 , 4 , 1 , 2 . On the other hand, the resistive grid on the right side starts at node five generating a connection tree compound of the resistors in red R 3 , R 2 , R 1 , R 9 and a path sequence 5 , 4 , 1 , 2 , 3 . In addition, the total lengths of each path shown in left-side and right-side Figure 7 are 15.89 and 19.65, respectively.
With the change of starting node, the possibility of exploring the different starting points in the resistive grid is opened, and therefore it is possible to determine which is the connection tree that reduces the total distance for the path between all the nodes of the resitive grid. Particularly, for this example, five different connection trees can be generated by changing the starting node.
Table 4 shows the five possible starting points for delivering merchandise. Based on the results for the length of the path that passes through all the nodes, the choice of node three as the starting point is the most convenient due to its smaller value.

3. Experimental Results

With the purpose of evaluating the feasibility of the RGPPM extension in order to deliver merchandise, Lecheria y Quesos Artesanales “Los Ángeles” has shared three maps of the locations for the deliveries of its products from three different cities: Puebla, Xalapa, and Veracruz-Boca del Rio. In these maps, the locations are transformed into nodes belonging to multi-connected resistive grid configurations in order to carry out the modified RGPPM algorithm and therefore establish the complete path for delivering merchandise.
For each experimental case, a partial distribution zone has been established in order to carry out the technique based on resistive grids. In addition, the search of two alternative full path resistive grid configurations have been explored using different starting nodes. Finally, the behavior of the computation time as the number of georeferential locations increases has been evaluated.

3.1. Puebla Logistical Map

In the city of Puebla, Lecheria y Quesos Artesanales “Los Ángeles” distributes its merchandise to the locations of fifty-four customers, as shown in Figure 8a. As it can be seen, the distribution of the merchandise delivery points contains disperse and concentrate areas.
With the purpose of showing the transformation from the distribution map to the resistive grid, a partial distribution zone has been selected as shown in Figure 8b. This set of eleven georeferential locations has been converted into a multi-connected resistive grid, as shown in Figure 9.
As it can be seen, the eleven geopositions have been interconnected by fifty-five resistors. In order to assign numerical values to resistors, it is necessary to calculate the distance between pairs of georeferential coordinates. In this part of the experimental results, the distance between two points on the earth’s surface has been estimated using the mathematical approximation of the Harvesine formula [24]
d = 2 r arcsin sin 2 ϕ 2 ϕ 1 2 + cos ( ϕ 1 ) cos ( ϕ 2 ) sin 2 λ 2 λ 1 2
where ϕ 1 , ϕ 2 , λ 1 , λ 2 are the latitude and longitude, expressed in radians, of the first and second georeferential points, respectively, and r is the terrestrial radius (6371 km approx.)
It is important to mention that the Haversine expression is more accurate for calculating longer distances between two points on the surface than two-dimensional Euclidean distances due to consideration of the curvature of the earth. Therefore, this addition opens the possibility to explore multi-connected irregular geometrical grid configurations with the greatest possible dispersion between customer locations.
After carrying out Equation (4), the distances between all the pairs of locations of the merchandise deliveries have been calculated and therefore the numerical values (in Ω ) of the fifty-five resistors for the multi-connected resistive grid from Puebla map have been established as shown in Table 5.
Hereinafter, the modified version of RGPPM begins iteratively searching for the path that goes through all the locations for delivering merchandise.
According to Figure 10, the path through the eleven locations for delivering merchandise contains the resistors in red R 1 , R 13 , R 28 , R 28 , R 41 , R 46 , R 51 , R 55 , R 54 , R 25 and the sequence of nodes in the complete path is 1 , 2 , 5 , 4 , 6 , 7 , 8 , 10 , 11 , 9 , 3 . Moreover, the total length of the path is 16.88 km.
Figure 11 shows two alternate full path resistive grid configurations using different starting nodes. In the case of the resistive grid on the left side, the starting node is the number seven, which generates a connection tree formed by the resistors in red R 46 , R 51 , R 44 , R 29 , R 28 , R 13 , R 1 , R 2 , R 25 , R 54 and a path sequence 7 , 8 , 10 , 6 , 4 , 5 , 2 , 1 , 3 , 9 , 11 . On the other hand, the resistive grid on the right side starts at node four generating a connection tree compound of the resistors in red R 29 , R 35 , R 13 , R 1 , R 2 , R 23 , R 46 , R 51 , R 55 , R 54 and a path sequence 4 , 6 , 5 , 2 , 1 , 3 , 7 , 8 , 10 , 11 , 9 . In addition, the total lengths of each path shown on the left-side and the right-side Figure 11 are 17.14 km and 16.21 km, respectively.
As it can be seen in Figure 12, the behavior of the CPU-time with respect to the number of merchandise delivery points has been successfully estimated. Start from the partial distribution zone (seven customers), the computation time, in order to transform the georeferential locations into a resistive grid, formulate the equilibrium equations, solve the system of equations, and find the path that runs through all the nodes that has been established according to increases in the number of customers by one until the maximum amount (fifty-four customers) has been achieved.
It is important to mention that, as the number of customers increases, the CPU-time increases, with the maximum computation time reached around 360 ms.

3.2. Xalapa Logistical Map

In the city of Xalapa, Lecheria y Quesos Artesanales “Los Ángeles” distributes its merchandise to the locations of twenty-five customers, as shown in Figure 13a. In the same way as in the city of Puebla, the location of the customers displays concentrated and disperse areas.
Figure 13b shows the partial distribution zone in order to be transformed from the distribution map to the resistive grid. In order to do this, the six georeferential locations have been converted into a multi-connected resistive grid, as shown in Figure 14.
After having converted the georeference locations into a multi-connected resistive grid, it is necessary to carry out the mathematical expression shown in Equation (4) in order to estimate the numerical values (in Ω ) of the fifteen resistors as shown in Table 6.
From this moment, the formulation of the equilibrium equations, the solution of the system of equations, and the path search procedure are carried out for establishing the path that goes through all the locations for delivering merchandise.
According to Figure 15, the path through the six locations for delivering merchandise contains the resistors in red R 1 , R 6 , R 28 , R 10 , R 13 , R 15 and a path sequence 1 , 2 , 3 , 4 , 5 , 6 . Therefore, it reaches a total distance of 6.84 km.
Figure 16 shows two alternate full path resistive grid configurations using different starting nodes. In the case of the resistive grid on the left side, the starting node is the number three, which generates a connection tree formed by the resistors in red R 10 , R 7 , R 1 , R 4 , R 15 and a path sequence 3 , 4 , 2 , 1 , 5 , 6 . On the other hand, the resistive grid on the right side starts at node five generating a connection tree compound of the resistors in red R 15 , R 12 , R 10 , R 7 , R 1 and a path sequence 5 , 6 , 3 , 4 , 2 , 1 . Moreover, the total lengths of each path shown on the left-side and right-side (Figure 16) are 8.62 km and 7.17 km, respectively.
With the same purpose as in the city of Puebla, the information in Table 7 shows that, as the number of customers increases, the computation time remains constant in most of the range of variation. The maximum computation time reached is 46.87 ms.

3.3. Veracruz-Boca del Rio Logistical Map

In the Veracruz-Boca del Rio metropolitan area, Lecheria y Quesos Artesanales “Los Ángeles” distributes its merchandise to the locations of sixty-four customers, as shown in Figure 17a.
In the same way that it has been carried out previously, a partial distribution zone is transformed into a resistive grid by converting five georeferential locations to a multi-connected resistive grid as shown in Figure 18.
Table 8 shows the numerical values (in Ω ) of the ten resistors calculated with the mathematical expression shown in Equation (4).
Next, the extension of RGPPM for green logistics is carried out in order to establish the path that passes through all the locations for delivering merchandise. Figure 19 shows the path through the five locations for delivering merchandise. It contains the resistors in red R 1 , R 5 , R 8 , R 10 , the path sequence 1 , 2 , 3 , 4 , 5 , and therefore reaches a total distance of 5.19 km.
According to the steps previously carried out in the maps of the cities of Puebla and Xalapa, the exploration of the alternate full path resistive grid configurations using different starting nodes and the behavior of the computation time as the number of clients increases must be accomplished. However, these results have not been shown because they are similar to the previous cases.
Finally, the experimental results confirm, through the exploration of three logistics maps for the delivery of merchandise from Lecheria y Quesos Artesanales “Los Ángeles” business, that the novel green logistic technique based on RGPPM extension is capable of establishing path trees connected to all georeferential points. Additionally, it has been shown that the technique is capable of changing the starting point and therefore re-generating new alternatives in order to deliver merchandise. Moreover, it has been possible to observe the behavior of the computing time as the number of clients increases.

4. Discussion

Although the research does not intend to carry out a numerical comparison between the traditional methods or techniques of easy computational implementation for the delivery of merchandise and the novel technique based on RGPPM, it is convenient to discuss the similarities and differences with respect to the principal characteristics of each technique. According to the experience of the authors and as far as the knowledge about the area reaches, it is possible to consider two fundamental methods: Dijkstra [25] and Kruskal [26] algorithms.
The classical Dijkstra algorithm [27] solves the single origin shortest path problem on a weighted graph. Because the problem of delivering merchandise to all hand-over points during the shortest continuous path is related to the connection of a weighted graph, Dijkstra’s algorithm can be considered as an alternative to attack this problem. However, the shortest path between two locations is not the same as the minimum spanning tree that connects all hand-over points in a continuous path. Therefore, Dijkstra’s original conception must be modified to address the search for the minimum weight spanning tree. At this point, it is important to mention that RGPPM has previously shown to be able to estimate the shortest path from one point to another for weighted graphs from resistive grids and as a consequence of the present investigation RGPPM can be expanded towards the problem of the search for the minimum spanning tree. Due to this, RGPPM represents an advantage since it is capable of tackling two different problems with the same original conception of the algorithm.
The Kruskal algorithm [28] generates a spanning tree with a minimal cost for the sum of all edge weights, i.e., a minimum spanning tree. However, this algorithm has been considered short-sighted for its estimates since it makes decisions using the available information without worrying about the effect that these decisions may cause in the future and therefore many problems may not have been solved correctly. According to the methodological structure of RGPPM, it is possible to establish an advantage compared to the Kruskal method because RGPPM uses the MNA formulation to concentrate all the information of the weighted graph in the form of a resistive grid and therefore be able to estimate the changes in the grid using theorems from circuit theory. This feature causes a focus-sighted in order to carry out decisions.
According to Figure 3, the multi-connected resistive grid configuration is capable of containing twenty-four different paths in order to traverse the locations for delivering merchandise. Table 9 shows the paths, the sequence of nodes, and the length of the path.
As it can be seen in the previous table, the path established by RGPPM contains the sequence of nodes 1 , 4 , 5 , 3 , 2 and therefore the length of the path is 18.12. Because the length of the path is related to the displacement of a terrestrial vehicle for delivering merchandise, it is possible to assert that RGPPM reduces the total distance of the path by around 30% compared to the longest path. Moreover, it is important to mention that reducing the length of the path has a considerable impact on different variables such as: amount of fuel, travel time, and emission of pollutants.
Based on simulation results, it is possible to determine that the location of the starting node on the multi-connected resistive grid is a fundamental piece for the length of the path in order to deliver merchandise. Table 10 shows the length of the paths and the percentage of reduction with respect to the longest path as the starting node changes in Figure 3.
As it can be seen in the table above, the reduction in the length of the path for delivering merchandise undergoes significant changes for different starting nodes. For example, the difference in percentage between the shortest path and the longest path during changes in the starting node has been estimated from 25% to 45%. Therefore, it is possible to reduce or increase the displacement path of a transport system for delivering merchandise by changing the initial location.
A strategical and decisive factor for the success or failure of a business is the location of a merchandise distribution center. According to experimental results from the RGPPM-based technique, it is possible to glimpse the feasible positions for locating a merchandise distribution center, as shown in Figure 20.
Figure 20 shows the complete distance for delivering merchandise by changing the start node. According to these results, it is possible to establish that the maximum distance is around 84 km, and the minimum distance is around 72 km. However, there are thirteen different starting nodes where the complete distance to deliver merchandise remains close to the minimum distance. It is important because it opens up the possibility of various alternatives for the location of a merchandise distribution center.
For the purpose of evidencing the feasibility of this new technique, in order to recalculate the tree connection with all the nodes based on the estimation of the best place to locate a merchandise distribution center, Figure 21 shows the multiconnected resistive grid (black lines), the actual location of the merchandise distribution center (orange marker), and the shortest path for delivering merchandise (red lines).
Figure 12 and Table 7 have shown the evolution of the computation time as the number of merchandise delivery points increases. Despite of observing that the CPU-time tends to increase as the number of georeferential points increases in the logistic map of the cities, it is possible to agree that the behavior of the curve and the data reflect a smooth nonlinear trend. It has been feasible because RGPPM uses the multifrontal method to solve sparse matrices with large dimensions.
Future research should be focused on the following topics:
  • Incorporating real information from the road route such as distance across transportation paths, traffic speed, waiting time at traffic lights, and others with the purpose of estimating statistical indicators so that the algorithm carries out decisions on the planning of merchandise deliveries.
  • Exploring the fusion of the RGPPM algorithm with Dijkstra and Bellman–Ford algorithms for assembling a hybrid and competitive strategy with the purpose of finding the path with the shortest distance that connects all locations in a continuous path.

5. Conclusions

A novel green logistics approach for organizing item delivery during pandemic events has been proposed. It is based on a modified version of RGPPM. This approach translates georeferential locations to a multi-connected resistive grid in order to search a path with a shorter distance than the possible paths in the logistical map. Additionally, the location of a merchandise distribution center has been estimated in a city using iterative exploration of all the starting points of the resistive grid. In addition, the technique maintains an acceptable performance with respect to computation time as the number of delivery points increases.
Due to the fact that the core of this novel technique is based on RGPPM, it has the uniqueness of relating the fundamental principles of circuit theory with planning and logistics problems and therefore contributes significantly to a new vision of the optimization process. The ability to process large amounts of georeferencing locations, adaptability in changing the initial conditions for the search of the minimun spanning tree, not high computation times for the exploration of logistic indicators and growth prospects as the technique matures can be considered as strengths. In this phase of the research, the limitations of the technique are centered on the inability to use multi-agent systems for the distribution of merchandise and the consideration of the distance of the trajectory as the unique variable in the delivery process.
Additionally, the most significant findings using this technique for the planning of merchandise deliveries can be summarized as follows:
Spanning trees that reduce around 45% of the total displacement with respect to the longest connecting tree have been achieved for the different scenarios;
Computation time for finding the shortest spanning tree is no more than half a second as the number of customers increases around the maximum number of locations for delivering merchandise;
Change of starting point has made it possible to explore the impact of the different spanning trees with respect to the total distance in order to deliver the merchandise;
Numerical measurement of spanning trees contributes to strategical location of distribution centers based on customer grids.
Finally, two recommendations for industry managers have been established as follows:
  • Experimental results of this research recommend those in charge of locating the merchandise distribution centers in the cities to methodologically explore the different positions of the clients for strategically estimating the shortest spanning tree in order to deliver merchandise;
  • Nowadays, companies must assume high environmental commitments by implementing best practices in order to operate without going over the environment. Managers of last mile deliveries can meet these environmental standards by using this planning technique since it decreases fuel consumption, and therefore the costs associated with delivery activities also decrease.

Author Contributions

Conceptualization, C.H.-M. and D.T.-M.; Data curation, E.E.G.-G.; Methodology, C.H.-M., D.T.-M., E.I.-G. and C.S.-L.; Formal analysis, C.H.-M., D.T.-M., E.I.-G. and C.S.-L.; Investigation, C.H.-M., D.T.-M., E.I.-G. and C.S.-L.; Resources, E.I.-G. and E.E.G.-G.; Visualization, C.H.-M., E.E.G.-G. and E.I.-G.; Writing—original draft preparation, C.H.-M. and D.T.-M.; Writing—review and editing, C.H.-M., D.T.-M. and C.S.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nagao, T.; Ijuin, H.; Yamada, T.; Nagasawa, K.; Zhou, L. COVID-19 Disruption Strategy for Redesigning Global Supply Chain Network across TPP Countries. Logistics 2022, 6, 2. [Google Scholar] [CrossRef]
  2. IPBES. Workshop Report on Biodiversity and Pandemics of the Intergovernmental Platform on Biodiversity and Ecosystem Services (IPBES); IPBES: Bonn, Germany, 2020. [Google Scholar]
  3. Figliozzi, M.; Unnikrishnan, A. Home-deliveries before-during COVID-19 lockdown: Accessibility, environmental justice, equity, and policy implications. Transp. Res. D Transp. Environ. 2021, 93, 102760. [Google Scholar] [CrossRef]
  4. Wang, X.C.; Kim, W.; Holguín-Veras, J.; Schmid, J. Adoption of delivery services in light of the COVID pandemic: Who and how long? Transp. Res. Part Policy Pract. 2021, 154, 270–286. [Google Scholar] [CrossRef] [PubMed]
  5. Li, C.; Mirosa, M.; Bremer, P. Review of Online Food Delivery Platforms and their Impacts on Sustainability. Sustainability 2020, 12, 5528. [Google Scholar] [CrossRef]
  6. Allen, J.; Piecyk, M.; Cherrett, T.; Juhari, M.N.; McLeod, F.; Piotrowska, M.; Bates, O.; Bektas, T.; Cheliotis, K.; Friday, A.; et al. Understanding the transport and CO2 impacts of on-demand meal deliveries: A London case study. Cities 2021, 108, 102973. [Google Scholar] [CrossRef]
  7. Wakeland, W.; Cholette, S.; Venkat, K. Food transportation issues and reducing carbon footprint. In Green Technologies in Food Production and Processing; Boye, J.I., Arcand, Y., Eds.; Springer US: Boston, MA, USA, 2012; pp. 211–236. [Google Scholar]
  8. Gleissner, H.; Femerling, J. Logistics: Basics-Exercises-Case Studies; Springer International Publishing: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  9. Rushton, A.; Croucher, P.; Baker, P. The Handbook of Logistics and Distribution Management: Understanding the Supply Chain; Kogan Page: London, UK, 2014. [Google Scholar]
  10. Psaraftis, H.N. Green Transportation Logistics; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  11. Fahimnia, B.; Bell, M.G.; Hensher, D.A.; Sarkis, J. Green Logistics and Transportation; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  12. Derbel, H.; Jarboui, B.; Siarry, P. Modeling and Optimization in Green Logistics; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  13. Zijm, H.; Klumpp, M.; Regattieri, A.; Heragu, S. Operations, Logistics and Supply Chain Management; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  14. Ranieri, L.; Digiesi, S.; Silvestri, B.; Roccotelli, M. A Review of Last Mile Logistics Innovations in an Externalities Cost Reduction Vision. Sustainability 2018, 10, 782. [Google Scholar] [CrossRef]
  15. Gutin, G.; Punnen, A.P. The Traveling Salesman Problem and Its Variations; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  16. Palhares, R.A.; AraÚjo, M.C.B. Vehicle Routing: Application of Travelling Salesman Problem in a Dairy. In Proceedings of the 2018 IEEE International Conference on Industrial Engineering and Engineering Management (IEEM), Bangkok, Thailand, 16–19 December 2018; pp. 1421–1425. [Google Scholar]
  17. Hernández-Mejía, C.; Vázquez-Leal, H.; Sánchez-González, A.; Corona-Avelizapa, Á. A Novel and Reduced CPU Time Modeling and Simulation Methodology for Path Planning Based on Resistive Grids. Arab. J. Sci. Eng. 2019, 44, 2321–2333. [Google Scholar] [CrossRef]
  18. Hernández-Mejía, C.; Hernández-Méndez, S.; Garrido-Grajales, M.; Corona-Avelizapa, A. Robotical Implementations of Collision-Free Path Planning for Terrestrial Systems using RGPPM. In Proceedings of the 2019 IEEE International Conference on Engineering Veracruz (ICEV), Boca del Rio, Mexico, 14–17 October 2019; Volume I, p. 16. [Google Scholar]
  19. Ho, C.W.; Ruehli, A.; Brennan, P. The modified nodal approach to network analysis. IEEE Trans. Circuits Syst. 1975, 22, 504–509. [Google Scholar]
  20. Ogrodzki, J. Circuit Simulation Methods and Algorithms; Electronic Engineering Systems; Taylor & Francis: Abingdon, UK, 1994. [Google Scholar]
  21. Schwarz, A. Computer-Aided Design of Microelectronic Circuits and Systems: General Introduction and Analog-Circuit Aspects; Academic Press: Cambridge, MA, USA, 1987. [Google Scholar]
  22. Davis, T.A. Algorithm 832: UMFPACK V4.3—An Unsymmetric-pattern Multifrontal Method. ACM Trans. Math. Softw. 2004, 30, 196–199. [Google Scholar] [CrossRef]
  23. Kanaya, M.; Cheng, G.X.; Watanabe, K.; Tanaka, M. Shortest path searching for robot walking using an analog resistive network. In Proceedings of the IEEE International Symposium on Circuits and Systems-ISCAS’94, London, UK, 30 May–2 June 1994; Volume 6, pp. 311–314. [Google Scholar]
  24. Lawhead, J. Learning Geospatial Analysis with Python; Packt Publishing: Birmingham, UK, 2019. [Google Scholar]
  25. Dijkstra, E.W. A note on two problems in connexion with graphs. Numer. Math. 1959, 1, 269–271. [Google Scholar] [CrossRef]
  26. Kruskal, J.B. On the shortest spanning subtree of a graph and the traveling salesman problem. Proc. Am. Math. Soc. 1956, 7, 48–50. [Google Scholar] [CrossRef]
  27. Edelkamp, S.; Schrödl, S. Heuristic Search-Theory and Applications; Academic Press: Cambridge, MA, USA, 2012. [Google Scholar]
  28. Graham, R.; Hell, P. On the History of the Minimum Spanning Tree Problem. Ann. Hist. Comput. 1985, 7, 43–57. [Google Scholar] [CrossRef]
Figure 1. Regular geometrical grid configurations from [17].
Figure 1. Regular geometrical grid configurations from [17].
Logistics 06 00059 g001
Figure 2. Multi-connected irregular geometrical grid configuration.
Figure 2. Multi-connected irregular geometrical grid configuration.
Logistics 06 00059 g002
Figure 3. Multi-connected resistive grid configuration.
Figure 3. Multi-connected resistive grid configuration.
Logistics 06 00059 g003
Figure 4. A segment of the grid configuration from Figure 3.
Figure 4. A segment of the grid configuration from Figure 3.
Logistics 06 00059 g004
Figure 5. A modified multi-connected resistive grid configuration.
Figure 5. A modified multi-connected resistive grid configuration.
Logistics 06 00059 g005
Figure 6. Resistive grid configuration with a complete path for delivering merchandise.
Figure 6. Resistive grid configuration with a complete path for delivering merchandise.
Logistics 06 00059 g006
Figure 7. Resistive grid configurations with alternative complete paths using a different starting node.
Figure 7. Resistive grid configurations with alternative complete paths using a different starting node.
Logistics 06 00059 g007
Figure 8. Distribution of deliveries in the city of Puebla.
Figure 8. Distribution of deliveries in the city of Puebla.
Logistics 06 00059 g008
Figure 9. A multi-connected resistive grid from Figure 8b.
Figure 9. A multi-connected resistive grid from Figure 8b.
Logistics 06 00059 g009
Figure 10. Path for delivering merchandise on the Puebla map.
Figure 10. Path for delivering merchandise on the Puebla map.
Logistics 06 00059 g010
Figure 11. Path for delivering merchandise on the Puebla map using a different starting node.
Figure 11. Path for delivering merchandise on the Puebla map using a different starting node.
Logistics 06 00059 g011
Figure 12. Performance of computing time for different amount of customers on the Puebla map.
Figure 12. Performance of computing time for different amount of customers on the Puebla map.
Logistics 06 00059 g012
Figure 13. Distribution of deliveries in the city of Xalapa.
Figure 13. Distribution of deliveries in the city of Xalapa.
Logistics 06 00059 g013
Figure 14. A multi-connected resistive grid from Figure 13b.
Figure 14. A multi-connected resistive grid from Figure 13b.
Logistics 06 00059 g014
Figure 15. Path for delivering merchandise on the Xalapa map.
Figure 15. Path for delivering merchandise on the Xalapa map.
Logistics 06 00059 g015
Figure 16. Path for delivering merchandise on the Xalapa map using different starting nodes.
Figure 16. Path for delivering merchandise on the Xalapa map using different starting nodes.
Logistics 06 00059 g016
Figure 17. Distribution of deliveries in the city of Veracruz-Boca del Rio.
Figure 17. Distribution of deliveries in the city of Veracruz-Boca del Rio.
Logistics 06 00059 g017
Figure 18. A multi-connected resistive grid from Figure 17b.
Figure 18. A multi-connected resistive grid from Figure 17b.
Logistics 06 00059 g018
Figure 19. Path for delivering merchandise on the Veracruz-Boca del Rio map.
Figure 19. Path for delivering merchandise on the Veracruz-Boca del Rio map.
Logistics 06 00059 g019
Figure 20. Path distances changing start node from the Veracruz-Boca del Rio map.
Figure 20. Path distances changing start node from the Veracruz-Boca del Rio map.
Logistics 06 00059 g020
Figure 21. Exploration of a strategic location for a merchandise distribution center in the Veracruz-Boca del Rio map.
Figure 21. Exploration of a strategic location for a merchandise distribution center in the Veracruz-Boca del Rio map.
Logistics 06 00059 g021
Table 1. Resistor values for the multi-connected grid configuration.
Table 1. Resistor values for the multi-connected grid configuration.
ResistorValue Ω
R 1 5.83
R 2 2.23
R 3 3.53
R 4 4.30
R 5 3.53
R 6 7.28
R 7 7.28
R 8 4.52
R 9 8.06
R 10 6.00
Table 2. Summary of voltages and currents.
Table 2. Summary of voltages and currents.
VariableValue
v 1 10.0 V
v 2 0 V
v 3 0 V
v 4 0 V
v 5 0 V
i E 1 −10.4 A
i E 2 1.7 A
i E 3 1.3 A
i E 4 4.4 A
i E 5 2.8 A
Table 3. Branch currents from Figure 4.
Table 3. Branch currents from Figure 4.
CurrentValue
i R 1 1.71 A
i R 5 2.83 A
i R 6 1.37 A
i R 2 4.48 A
Table 4. Performance of the paths from Figure 3 by changing the start node.
Table 4. Performance of the paths from Figure 3 by changing the start node.
Starting NodeSequence of Nodes Using RGPPMLength
1 1 , 4 , 5 , 3 , 2 18.12
2 2 , 5 , 1 , 4 , 3 16.28
3 3 , 5 , 4 , 1 , 2 15.89
4 4 , 1 , 5 , 3 , 2 18.12
5 5 , 4 , 1 , 2 , 3 19.65
Table 5. Resistor values for the Puebla map.
Table 5. Resistor values for the Puebla map.
ResistorValueResistorValueResistorValueResistorValueResistorValue
R 1 0660 R 12 1460 R 23 2840 R 34 4950 R 45 4030
R 2 2280 R 13 1380 R 24 3630 R 35 1650 R 46 0810
R 3 1720 R 14 2470 R 25 4280 R 36 3360 R 47 2730
R 4 2010 R 15 3960 R 26 4220 R 37 3680 R 48 1590
R 5 2850 R 16 4490 R 27 5110 R 38 6010 R 49 2470
R 6 4150 R 17 6310 R 28 1360 R 39 3560 R 50 2880
R 7 4770 R 18 4590 R 29 1170 R 40 5650 R 51 0900
R 8 6310 R 19 6400 R 30 2490 R 41 1720 R 52 1970
R 9 4980 R 20 1820 R 31 3060 R 42 2070 R 53 3710
R 10 6620 R 21 3140 R 32 4900 R 43 4410 R 54 2160
R 11 2570 R 22 2580 R 33 3280 R 44 2140 R 55 2360
Table 6. Resistor values for Xalapa map.
Table 6. Resistor values for Xalapa map.
ResistorValueResistorValueResistorValue
R 1 1620 R 6 0870 R 11 1730
R 2 2370 R 7 1510 R 12 1890
R 3 2830 R 8 1840 R 13 2200
R 4 3340 R 9 2170 R 14 2240
R 5 3720 R 10 0680 R 15 1470
Table 7. Performance of CPU-time for different amounts of customers in the Xalapa map.
Table 7. Performance of CPU-time for different amounts of customers in the Xalapa map.
Number of CustomersCPU-Time (ms)Number of CustomersCPU-Time (ms)
615.621615.62
715.621731.25
815.621831.25
915.621931.25
1015.622031.25
1115.622131.25
1215.622231.25
1315.622331.25
1415.622446.87
1515.622546.87
Table 8. Resistor values for the Veracruz-Boca del Rio map.
Table 8. Resistor values for the Veracruz-Boca del Rio map.
ResistorValueResistorValue
R 1 1060 R 6 2010
R 2 1830 R 7 2250
R 3 1590 R 8 2730
R 4 1530 R 9 3020
R 5 0800 R 10 0600
Table 9. Performance of the paths from Figure 3.
Table 9. Performance of the paths from Figure 3.
PathSequence of NodesLength
1 1 , 2 , 3 , 4 , 5 23.42
2 1 , 2 , 3 , 5 , 4 21.72
3 1 , 2 , 4 , 3 , 5 23.41
4 1 , 2 , 4 , 5 , 3 20.94
5 1 , 2 , 5 , 3 , 4 20.65
6 1 , 2 , 5 , 4 , 3 19.88
7 1 , 3 , 2 , 4 , 5 26.15
8 1 , 3 , 2 , 5 , 4 23.39
9 1 , 3 , 4 , 2 , 5 25.08
10 1 , 3 , 4 , 5 , 2 21.33
11 1 , 3 , 5 , 2 , 4 23.38
12 1 , 3 , 5 , 4 , 2 22.39
13 1 , 4 , 2 , 3 , 5 21.87
14 1 , 4 , 2 , 5 , 3 18.33
15 1 , 4 , 3 , 2 , 5 20.81
16 1 , 4 , 3 , 5 , 2 17.05
17 1 , 4 , 5 , 2 , 3 18.34
18 1 , 4 , 5 , 3 , 2 18.12
19 1 , 5 , 2 , 3 , 4 22.11
20 1 , 5 , 2 , 4 , 3 21.33
21 1 , 5 , 3 , 2 , 4 23.17
22 1 , 5 , 3 , 4 , 2 21.11
23 1 , 5 , 4 , 2 , 3 22.40
24 1 , 5 , 4 , 3 , 2 21.12
Table 10. Performance of the paths from Figure 3 by changing the start node.
Table 10. Performance of the paths from Figure 3 by changing the start node.
Starting NodeSequence of Nodes Using RGPPMPorcentage of Reduction
1 1 , 4 , 5 , 3 , 2 30.70%
2 2 , 5 , 1 , 4 , 3 32.42%
3 3 , 5 , 4 , 1 , 2 44.47%
4 4 , 1 , 5 , 3 , 2 30.70%
5 5 , 4 , 1 , 2 , 3 24.85%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hernández-Mejía, C.; Torres-Muñoz, D.; Inzunza-González, E.; Sánchez-López, C.; García-Guerrero, E.E. A Novel Green Logistics Technique for Planning Merchandise Deliveries: A Case Study. Logistics 2022, 6, 59. https://doi.org/10.3390/logistics6030059

AMA Style

Hernández-Mejía C, Torres-Muñoz D, Inzunza-González E, Sánchez-López C, García-Guerrero EE. A Novel Green Logistics Technique for Planning Merchandise Deliveries: A Case Study. Logistics. 2022; 6(3):59. https://doi.org/10.3390/logistics6030059

Chicago/Turabian Style

Hernández-Mejía, Carlos, Delia Torres-Muñoz, Everardo Inzunza-González, Carlos Sánchez-López, and Enrique Efrén García-Guerrero. 2022. "A Novel Green Logistics Technique for Planning Merchandise Deliveries: A Case Study" Logistics 6, no. 3: 59. https://doi.org/10.3390/logistics6030059

APA Style

Hernández-Mejía, C., Torres-Muñoz, D., Inzunza-González, E., Sánchez-López, C., & García-Guerrero, E. E. (2022). A Novel Green Logistics Technique for Planning Merchandise Deliveries: A Case Study. Logistics, 6(3), 59. https://doi.org/10.3390/logistics6030059

Article Metrics

Back to TopTop