Next Article in Journal
Efficient Numerical Solutions for Fuzzy Time Fractional Convection Diffusion Equations Using Two Explicit Finite Difference Methods
Next Article in Special Issue
Generalization of Fermatean Fuzzy Set and Implementation of Fermatean Fuzzy PROMETHEE II Method for Decision Making via PROMETHEE GAIA
Previous Article in Journal
Bidiagonal Factorizations of Filbert and Lilbert Matrices
Previous Article in Special Issue
Multi-Objective Non-Linear Programming Problems in Linear Diophantine Fuzzy Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving a Multimodal Routing Problem with Pickup and Delivery Time Windows under LR Triangular Fuzzy Capacity Constraints

1
SDU-ANU Joint Science College, Shandong University, Weihai 264209, China
2
School of Management Science and Engineering, Shandong University of Finance and Economics, Jinan 250014, China
*
Author to whom correspondence should be addressed.
Submission received: 28 February 2024 / Revised: 21 March 2024 / Accepted: 22 March 2024 / Published: 26 March 2024
(This article belongs to the Special Issue Recent Developments in Fuzzy Control Systems and Their Applications)

Abstract

:
This study models a container routing problem using multimodal transportation to improve its economy, timeliness, and reliability. Pickup and delivery time windows are simultaneously formulated in optimization to provide the shipper and the receiver with time-efficient services, in which early pickup and delayed delivery can be avoided, and nonlinear storage periods at the origin and the destination can be minimized. Furthermore, the capacity uncertainty of the multimodal network is incorporated into the advanced routing to enhance its reliability in practical transportation. The LR triangular fuzzy number is adopted to model the capacity uncertainty, in which its spread ratio is defined to measure the uncertainty level of the fuzzy capacity. Due to the nonlinearity introduced by the time windows and the fuzziness from the network capacity, this study establishes a fuzzy nonlinear optimization model for optimization problem. A chance-constrained linear reformulation equivalent to the proposed model is then generated based on the credibility measure, which makes the global optimum solution attainable by using Lingo software. A numerical case verification demonstrates that the proposed model can effectively solve the problem. The case analysis points out that the formulation of pickup and delivery time windows can improve the timeliness of the entire transportation process and help to achieve on-time transportation. Furthermore, improving the confidence level and the uncertainty level increases the total costs of the optimal route. Therefore, the shipper and the receiver must prepare more transportation budget to improve reliability and address the increasing uncertainty level. Further analysis draws some insights to help the shipper, receiver, and multimodal transport operator to organize a reliable and cost-efficient multimodal transportation under capacity uncertainty through confidence level balance and transportation service and transfer service selection.

1. Introduction

The operational-level multimodal routing plays a crucial role in improving the operations and management of multimodal transportation systems [1] and is a promising approach to fully strengthen the advantages of multimodal transportation over unimodal transportation. In China, by enhancing the service level to attract more freight sources, multimodal routing enables multimodal transportation to replace road-dominated transportation that is low efficiency and energy-consuming to optimize the transportation structure [2]. Due to its great significance in practice, the multimodal routing problem (MRP) has been widely discussed in the research field of transportation planning [3] and is one of the six research topics in the decision support of multimodal transportation (policy support, terminal network design, multimodal service network design, multimodal routing, drayage operations, and ICT innovations) [4].
Modeling MRP provides quantitative support for the decision-making of the transportation actors and stakeholders. MRP with the lead time constraint can avoid delayed delivery of containers at the destination that would result in economic loss. Shahedi et al. [5] build a multi-objective mixed integer linear programming model to formulate the lead time-based multimodal freight routing problem considering the Physical Internet. Verma and Verter [6] explore a rail–truck multimodal routing problem of dangerous goods with lead time constraints, in which a bi-objective optimization model is proposed. Although the delayed delivery of containers can be avoided, the lead time constraint might result in early delivery that causes post-transportation storage. Therefore, increasing multimodal routing studies focus on the MRP with time window(s), in which the time window uses its lower bound and upper bound separately to reduce or avoid both lateness and earliness of the container delivery at the destination. Therefore, compared to the lead time constraint, the concept of a time window is more suitable for achieving time-efficient transportation, especially in the context that a just-in-time strategy is resorted by companies, and transportation at the right time is required by them to minimize storage [7,8].
Both hard and soft time windows are established by the MRP literature. Those with the hard time window force that the container delivery by the multimodal route at the destination must be accomplished within the time window and propose the upper and lower bound constraints. Currently, there are a few studies using such a kind of time window. Zhang et al. [9] consider a green multimodal service network design problem with hard time window constraints, construct an integer linear optimization model, and develop genetic and heuristic algorithms. However, a hard time window is too strict and cannot reflect the customers’ attitudes that early and delayed transportation is acceptable to a certain degree. Consequently, the hard time window reduces the flexibility of the MRP. As a result, most MRP studies are interested in soft time windows that make the earliness and lateness allowable on condition that penalty or storage costs must be paid when the time window is violated, e.g., Zhang et al. [10], Fazayeli et al. [11], Yuan et al. [12], and Li et al. [13]. In this case, the penalty or storage period is calculated by the continuous piecewise linear function that introduces nonlinearity into the modeling of MRP. In transportation practice, delayed transportation could disrupt the further processing of the goods and could cause a long-term economic loss that is difficult to determine [14]. Therefore, the penalty costs in the soft time window cannot be given accurately. As a result, both hard and soft time windows yield disadvantages that reduce their feasibility when applied in practical transportation.
Furthermore, multimodal transportation starts with the container pickup from the shipper and ends with container delivery to the receiver. Therefore, both pickup and delivery determine the time efficiency of the entire multimodal transportation process [15]. However, the majority of relevant studies focus more on the delivery than the pickup when modeling the time window and assume that containers are released at a fixed time [13,16] or there is the earliest release time [17]. Currently, quite a few works simultaneously model the two types of time windows, in which Sun et al. [2] and Zhang et al. [18] use the soft time window. Although these studies discuss the pickup and delivery time windows, their modeling still has the disadvantage of the soft time windows indicated above.
In addition to the time window, uncertainty is a key direction of the MRP that seeks to improve the reliability of the multimodal route in practice [19] and has been highlighted by supply chain management [20]. The constantly changing status of the dynamic transportation network makes it difficult to accurately predict the real-time state of the parameters of the multimodal network in advanced transportation planning, which causes uncertainty. Formulating the uncertainty makes advanced multimodal routing more reliable in practical transportation where the real-time state of the multimodal network is eventually known. Currently, transportation planning studies pay attention to the capacity uncertainty of the transportation network [21]. Regarding MRP, the optimal route is infeasible if the capacities are overestimated and might not be the best solution if the capacities are underestimated. Compared to the MRP articles with deterministic capacity constraints, there are limited studies considering capacity uncertainty. Huang [22] assumes that the capacities of water and rail transportation obeys the normal distribution and designs a two-stage stochastic programming model for the multimodal routing problem in the stochastic environment. Uddin et al. [23,24] utilize the scenario-based approach to model the uncertain capacities of rail transportation and terminals and build stochastic programming models for the MRP. Contrary to stochastic programming, based on the fuzzy set theory [25], Sun et al. [8], Li et al. [13], and Lu et al. [26] use triangular fuzzy numbers to model capacity uncertainty and establish the fuzzy chance-constrained programming (FCCP) models to deal with their MRPs. A group of studies by Sun et al. [2,16,27] model trapezoidal fuzzy capacities and propose FCCP models for specific MRPs.
All these papers cited above state the significant influence of capacity uncertainty on the optimization of the MRP and promote the research progress on MRP under capacity uncertainty. However, the stochastic models taken by References [22,23,24] need to use large numbers of data to fit the probability distribution of the uncertain capacity, while there are not enough data in most cases [28]. The scenario-based approach to model the uncertain capacity as large numbers of scenarios challenges the computational efficiency of problem-solving [11]. Moreover, assuming that the capacity uncertainty obeys a certain distribution might not match practice. Under the above situation, fuzzy programming based on the use of fuzzy numbers to model uncertainty [29,30] is an alternative for MRP under capacity uncertainty. However, References [2,8,13,16,26,27] use regular fuzzy numbers when modeling uncertain capacity and only analyze the influence of the confidence level introduced by fuzzy programming approaches on the MRP. In these studies, the influence of the uncertainty level of the fuzzy capacity on the MRP is not considered. Additionally, apart from Li et al. [13], the studies cited above only formulate the uncertainty of the travel capacity of rail transportation, while taking the travel capacities of other transportation modes to be deterministic, together with the node transfer capacities between different transportation modes.
Consequently, there is a need to design time windows that can avoid inestimable economic loss and simultaneously maintain a certain flexibility. It is also important to use such time windows to improve the quality of both pickup and delivery services of multimodal transportation. Secondly, modeling the capacity uncertainty in MRP should cover the entire multimodal network, including the travel process on arcs and the transfer process at the nodes, which makes the modeling of capacity uncertainty complete and enhances the feasibility of the MRP under capacity uncertainty. Finally, modeling the capacity uncertainty should enable its uncertainty level to be measurable, so that we can analyze the influence of the uncertainty level on the MRP and draw insights to deal with uncertainty to organize reliable multimodal transportation. All of these three aspects are not considered comprehensively in the relevant MRP studies. Faced with the research gaps mentioned above, this study continues to investigate the multimodal routing problem with pickup and delivery time windows under capacity uncertainty (MRPPDTWCU) by making the following contributions. Specifically, this study is an extension of the study by Li et al. [13] in the following three aspects:
(1)
Both the pickup and delivery time windows are incorporated into the MRP to improve the time efficiency of the entire transportation process, in which early pickup and delayed delivery are forbidden, and the storage caused by delayed pickup and early delivery is minimized. Through this setting, the weaknesses of both soft and hard time windows are fixed.
(2)
The uncertainty of both travel capacities of transportation modes on the arcs and transfer capacities between different transportation modes at the nodes of the multimodal network is modeled by LR triangular fuzzy numbers (LRTFNs), in which we define the uncertainty level of the LR triangular fuzzy capacity (LRTFC).
(3)
Based on the chance-constrained linear programming (CCLP) model using the credibility measure, we analyze the influence of the confidence level and the uncertainty level on the MRPPDTWCU and summarize the insights that help multimodal transportation to deal with the uncertain environment.
The remaining sections of this study are organized as follows. In Section 2, we introduce the problem scenario, the formulation of pickup and delivery time windows, and the modeling of uncertain capacity in the MRPPDTWCU. In Section 3, a fuzzy optimization model is initially constructed based on the specific problem description. In Section 4, we process the model in Section 3 to remove its fuzziness and nonlinearity to make its global optimum solution attainable for the shipper, receiver, and multimodal transport operator (MTO), which supports their decision-making. In Section 5, a numerical case based on a previous study is presented to verify the feasibility of the model, and further reveal some insights that help to organize a cost-efficient and reliable multimodal transportation. Finally, we summarize the conclusions in Section 6.

2. Problem Description

In multimodal transportation, the MTO is responsible for the shipper and the receiver (that is, the customers), and coordinates the carriers of rail, road, and water transportation to carry out the origin-to-destination transportation based on the customer demands [31]. Therefore, it is the work of the MTO that plans the multimodal route, and therefore the MRP is oriented on the MTO’s viewpoint. The MTO needs to satisfy the customer demands by presenting an optimal route with high economy, timeliness, and reliability or a balanced state on these objectives. Among these objectives, improving the transportation economy by minimizing transportation-related costs is the most important consideration for customers to improve the profits of their products and maintain competition in the market. As a result, the objective of the MRP explored in this study is to minimize the total costs for multimodal transportation.
In this study, both the shipper and the receiver have a time window. For the container shipper, the lower bound of the pickup time window is the earliest release time of the containers at the origin. The container pickup time should not be earlier than that time. When the container pickup is delayed and later than the upper bound of the pickup time window, containers should be stored first before pickup. In this case, the shipper should pay for the storage. The receiver also has a fixed delivery time window. The upper bound is the latest delivery time, and the delivery of containers at the destination must be no later than that time to avoid disruption in the further processing of the containers and the resulting economic loss that is difficult to determine. However, container delivery is allowed to be accomplished earlier than the lower bound of the delivery time window, under the condition that the receiver needs to pay for the storage before the containers receive further processing. Under the above time window setting, improving the on-time pickup and delivery to optimize the timeliness of the entire multimodal transportation process can be achieved by reducing the storage costs covered in the optimization objective and following the time window constraints of the MRPPDTWCU. Under the above consideration, the MRPPDTWCU not only aims to give a physical path that connects the origin and the destination but also needs to plan the container pickup time at the origin.
In the MRPPDTWCU that should be conducted before the beginning of practical transportation, this study fully formulates the capacity uncertainty of the entire multimodal network. This study uses fuzzy set theory to build the uncertain capacity and constructs the fuzzy optimization model for the MRPPDTWCU. We adopt LR fuzzy numbers to model the fuzzy capacity. A fuzzy number g ~ is an LR fuzzy number if it has membership functions L for left and R for right, and left hand and right hand spreads α and β that are both positive numbers [32,33]. The membership function of LR fuzzy number g ~ is as Equation (1).
μ g ~ x = L g x α if   x g R x g β if   x g
In Equation (1), g is the mean value of g ~ . Symbolically, g ~ can be represented by g , α , β L R . L g x α and R x g β are non-increasing functions, and there are L 0 = 1 and R 0 = 1 . In this study, we consider the above two membership functions to be linear. Then, g ~ = g , α , β L R is an LRTFN [34,35] and can be used to express fuzzy capacity. In this expression, g is the mean value of the uncertain capacity, and α and β reflect the uncertainty levels that the multimodal network suffers capacity insufficiency (the worst condition) and is under capacity sufficiency (the best condition). Consequently, the utilization of the LRTFC enables the advanced MRPPDTWCU to consider all the possible capacity conditions that the multimodal network can face in practical transportation. Furthermore, we define α g and β g as the left hand and right hand spread ratios of the LRTFN. It is obvious that bigger α g and β g increase the uncertainty level of LRTFN g ~ . Compared to the regular TFN, with the help of the LRTFN to model the fuzzy capacity, we can analyze the influence of the uncertainty level of the fuzzy capacity on the MRPPDTWCU.

3. Fuzzy Nonlinear Optimization Model for MRPPDTWCU

3.1. Symbols in the Model

(1)
Sets, indices, and parameters:
N : set of nodes;
A : set of arcs;
S : set of transportation modes;
i , j , k : indices of transportation nodes, and i , j , k N ;
Φ j : set of predecessor nodes to node j , and Φ j N ;
Ψ j : set of successor nodes to node j , and Ψ j N ;
i , j : arc from node i to node j , and i , j A ;
S i j : set of transportation modes on arc i , j , and S i j S ;
S j : set of transportation modes linking node j , and S j S ;
m , n : indices of transportation modes, and m , n S ;
l i j m : travel distance in km of transportation mode m on arc i , j ;
t i j m : travel time in h of transportation mode m on arc i , j ;
t j m n : transfer time in h/TUE from transportation mode m to transportation mode n at node j ;
g ~ i j m = g i j m , α i j m , β i j m : LRTFC in TEU of transportation mode m on arc i , j ;
g ~ j m n = g j m n , α j m n , β j m n : LRTFC in TEU of node j to transfer containers from transportation mode m to transportation mode n ;
c i j m : travel costs in CNY/TEU of transportation mode m on arc i , j ;
c j m n : transfer costs in CNY/TEU from transportation mode m to transportation mode n at node j ;
o : origin of container transportation, and o N ;
d : destination of container transportation, and d N ;
q : transportation volume in TEU of containers;
c o s t o r e : storage costs in CNY/TEU·h at the origin;
c d s t o r e : storage costs in CNY/TEU·h at the destination;
t w o , t w o * : pickup time window of the shipper;
t w d , t w d * : delivery time window of the receiver;
(2)
Variables
x i j m : 0–1 variable. If transportation mode m moves containers on arc i , j , x i j m = 1 ; otherwise, x i j m = 0 ;
y j m n : 0–1 variable. If there is a transfer of containers from transportation mode m to transportation mode n at node j , y j m n = 1 ; otherwise, y j m n = 0 ;
z o : non-negative variable meaning the container pickup time at the origin;
z d : non-negative variable meaning the container delivery time at the destination;
w o : non-negative variable in h meaning the storage period of the containers at the origin;
w d : non-negative variable in h meaning the storage period of the containers at the destination.

3.2. Fuzzy Nonlinear Optimization Model

min i ,   j A m S i j c i j m · q · x i j m + j N m S j n S j c j m n · q · y j m n + c o s t o r e · q · w o + c d s t o r e · q · w d
s.t.
i Φ j m S i j x i j m k Ψ j n S j k x j k n = 1 j = o 0 j N \ o ,   d 1 j = d
m S i j x i j m 1 ( i , j ) A
m S j n S j y j m n 1 j N \ o ,   d
i Φ j x i j m = n S j y j m n j N \ o ,   d m S j
m S j y j m n = k Ψ j x j k n j N \ o ,   d n S j
z o t w o
w o = z o t w o * i f   z o > t o * 0 i f   z o t o *
z d = z o + i ,   j A m S i j t i j m · x i j m + j N m S j n S j t j m n · q · y j m n
z d t w d *
w d = t w d z d if   t w d > z d 0 if   t w d z d
x i j m · q g ~ i j m ( i , j ) A m S i j
y j m n · q g ~ j m n j N \ o ,   d m S j n S j
x i j m 0,1 ( i , j ) A m S i j
y j m n 0,1 j N \ o ,   d m S j n S j
z o 0
z d 0
w o 0
w d 0
  • Equation (2) is the optimization objective of the MRPPDTWCU and aims to minimize the total costs of transportation that include travel costs, transfer costs, and storage costs at the origin and the destination.
  • Equation (3) is the container flow equilibrium constraint.
  • Equations (4) and (5) ensure that containers are unsplittable in the transportation process from the shipper to the receiver.
  • Equations (6) and (7) ensure that the optimal route yields a smooth connection between the travel process on the selected arcs and the transfer process at the selected nodes. (Equations (3)–(7) are the general constraints of the MRP considering unsplittable flow [13]).
  • Equation (8) ensures that container pickup time at the origin should be no earlier than the lower bound of the shipper’s pickup time window.
  • Equation (9) calculates the storage period of the containers at the origin using a continuous piecewise linear function.
  • Equation (10) determines the container delivery time.
  • Equation (11) ensures that the container delivery time at the destination should be no later than the upper bound of the receiver’s delivery time window.
  • Equation (12) uses the same function as Equation (9) to present the storage period of the containers at the destination.
  • Equations (13) and (14) ensure that the container volume does not exceed the LRTFCs of the optimal route.
  • Equations (15)–(20) are the variable domain constraints.

4. Model Defuzzification and Linearization

The fuzzy nonlinear optimization model proposed in Section 3.2 is not solvable since its capacity constraints contain imprecise information. Therefore, we should remove the fuzziness of the model to make its solution attainable. As a widely used defuzzification approach, FCCP constructs the chance constraints of its fuzzy constraints, in which the confidence level is introduced into the fuzzy optimization model [36]. The value of the confidence level represents the attitudes of decision-makers on the constraints. It should be predetermined before solving the FCCP model.
The establishment of an FCCP model should first determine the type of fuzzy measure. Compared to the possibility and necessity measures separately representing the extremely optimistic and pessimistic attitudes, the credibility measure is a compromised attitude. It shows the attitudes of decision-makers that are neither too optimistic nor too pessimistic and matches the practical decision-making situation. Additionally, the credibility measure is self-dual and can guarantee that the fuzzy event must hold when its credibility reaches 1.0, and must fail when the credibility is equal to 0 [37]. Under the above consideration, we utilize the credibility measure to construct the chance constraints of Equations (13) and (14) as follows:
C r x i j m · q g ~ i j m δ ( i , j ) A m S i j
C r y j m n · q g ~ j m n δ j N \ o ,   d m S j n S j
As indicated by the above equations, the credibility (Cr) that the transportation volume of containers does not exceed the LRTFCs of the route should not be less than a minimum confidence level δ . Increasing the value of δ makes the optimal route more credible in practical transportation where the actual capacities emerge, as it improves the credibility that the optimal route satisfies the constraints of the actual capacities. As a result, with the increase in δ , the optimal route yields higher reliability in practical transportation. Although δ falls into a closed interval 0,1 , in practical routing decision-making, decision-makers prefer a δ that belongs to 0.5,1 [38]. As an interactive parameter, δ should be set by the decision-makers, including the shipper, receiver, and MTO.
Suppose a positive LR triangular fuzzy g ~ = g , α , β L R where g α > 0 and a positive deterministic number τ . The distribution of the credibility that g ~ τ is as Equation (23) [39].
C r g ~ τ = 1 if   g α τ g + α τ 2 α if   g α τ g g + β τ 2 β if   g τ g + β 0 if   τ g + β
Therefore, given a confidence level δ 0.5,1 , the crisp reformulation of C r g ~ τ δ can be realized by Equation (24).
C r g ~ τ δ g + 1 2 δ α τ
Accordingly, Equations (21) and (22) can be reformulated as Equations (25) and (26), respectively.
g i j m + 1 2 δ α i j m x i j m · q ( i , j ) A m S i j
g j m n + 1 2 δ α j m n y j m n · q j N \ o ,   d m S j n S j
Above all, we can obtain an FCCP model by replacing Equations (13) and (14) with Equations (25) and (26). This model is solvable since it removes the fuzziness. However, it is still a nonlinear optimization model due to the nonlinearity of Equations (9) and (10). To further improve the efficiency of finding the global optimum solution for the MRPPDTWCU, this study linearizes the two nonlinear constraints that use continuous piecewise linear functions as follows [13,16,27]:
w o z o t w o *
w d t w d z d
Although the two equivalent linear equations are not straightforwardly based on the definitions of the variables, they can significantly improve the computational efficiency of problem-solving [40]. Finally, after the model processing, we present a CCLP model that takes Equation (2) as the objective and covers Equations (3)–(8), (10), (11), (15)–(20) and (25)–(28) as the constraints. The global optimum solution to the problem can be obtained by exact solution algorithms.

5. Numerical Case Verification

This section presents a numerical case in the Chinese scenario to verify the proposed model. The multimodal network presented in Sun and Lang’s work [41] has been commonly modified and further used to carry out case verification (e.g., Li et al. [13] and Xiong and Wang [42]) and provide a foundation for us to design the case in this study. The multimodal network used to serve the container transportation from the shipper to the receiver is shown in Figure 1. The travel distances of the rail, road, and water transportation on the arcs are set as those given by Sun and Lang [41].
The travel costs and time of the three transportation modes are listed in Table 1, and the transfer time and costs between these transportation modes are presented in Table 2.
According to the specific MRPPDTWCU, we assume that the capacities in Sun and Lang [41] are the mean values of the LRTFCs of the multimodal network in this numerical case, which also makes the difference in the numerical case design between this study and Li et al. [13] more significant. The MTO serves a transportation order whose shipper is in Node 1 and receiver is in Node 35. The transportation volume of containers is 40 TEU. The pickup time window of the shipper is from 8:00 to 12:00 on day 1, and the delivery time window of the receiver is from 21:00 on day 2 to 3:00 on day 3. The storage costs at the origin and the destination are CNY 10/TEU·h and CNY 20 /TEU·h, respectively.
We assume that the left hand spread ratio is 20%, i.e., the spread is 20% of the mean values of the LRTFCs of the multimodal network. The shipper and the receiver prefer high reliability of the transportation of their containers and determine that the confidence level should be 0.9. In this case, we use Lingo software to run the branch-and-bound algorithm to solve the CCLP model. We obtain the optimal route “1-rail-3-road-6-rail-11-road-15-rail-16-road-21-rail-27-rail-28-road-35” and the container pickup time at Node 1 that is 11:33 on day 1. The total costs of using this route are CNY 289,341. The total costs consist of the travel costs of CNY 287,941 and transfer costs of CNY 1400, which means that the pickup and delivery of the containers are both within the corresponding time windows, and the optimal route can achieve on-time transportation. In this case, the pickup and delivery time windows have the same function as hard time windows. When the left hand spread ratio is 5% and the confidence level is 1.0, the optimal route with pickup and delivery time windows formulated in this study yields total costs of CNY 269,820, including the storage costs of CNY 6374. If the pickup and delivery time windows are the hard ones, the total costs of the optimal route are CNY 280,537 without any storage costs. Compared to the time windows designed in this study, the hard time windows increase the total costs of the optimal route by approximately 4%. Therefore, under this situation, the pickup and delivery time windows in this study can achieve a balance between the travel and transfer costs and the storage costs to minimize the total costs consisting of the above three parts.
In the numerical case, we further analyze how the spread ratio and confidence level influence the MRPPDTWCU. First of all, we set the left hand spread ratio as 25% and obtained the optimization results of the MRPPDTWCU under the values of the confidence level from 0.5 to 1.0 with a step of 0.1. The results are indicated in Figure 2.
As we can see in Figure 2, with the enhancement of the confidence level, the total costs of the optimal route increase, because the increase in confidence level makes the chance constraints (i.e., Equations (25) and (26)) tighter, and further suppresses the feasible solution space of the MRPPDTWCU. Consequently, the best solution given under a smaller confidence level might be found to be infeasible by the model under a larger confidence level. Improving the confidence level helps to make the optimal route more reliable in practical transportation. Therefore, the economy and reliability of the MRPPDTWCU are conflicting. The decision-makers in multimodal transportation including the MTO, shipper, and receiver should make tradeoffs between the two objectives. The sensitivity shown in Figure 2 can provide a reference for them to make a decision.
The majority of the relevant studies consider multimodal routing under capacity certainty. These studies prefer to use the mean values of the LRTFNs that represent the most likely conditions to formulate the deterministic capacities [8]. In this case, the optimal route under capacity certainty in the numerical case yields total costs of CNY 239,668, which equals the one under fuzzy capacity constraints when the confidence level is 0.5. However, the modeling of MRP under capacity uncertainty using the methods proposed by this study can provide diverse solutions that depend on the shipper, receiver, and MTO’s selection of the confidence level based on their attitudes. The solutions also include the one under capacity certainty that yields a low credibility, which means a low reliability in practical transportation. Consequently, compared to the deterministic modeling, modeling the MRP under capacity uncertainty can provide the shipper, receiver, and MTO with more optimal route schemes and thereby enable their decision-making to be more flexible, in which it is possible to make tradeoffs between lowering the costs and improving the reliability.
Then we analyze the influence of the left hand spread ratio on the MRPPDTWCU. In this case, we set the confidence level as 0.9 and present the optimization results of the MRPPDTWCU under different left hand spread ratios from 5% to 30% with a step of 5%. The results are shown in Figure 3.
Figure 3 clarifies that improving the spread ratio increases the total costs of the optimal route, which shows a similar change to the variation of the total costs caused by the increase in the confidence level. The reason resulting in the sensitivity shown in Figure 3 is the same as that leading to the variation illustrated by Figure 2, i.e., improving the left hand spread ratio increases α i j m and α j m n , thereby making the chance constraints tighter and the feasible solution space smaller, and finally increasing the objective values of the optimization model. Since the spread ratio indicates the uncertainty level of the LRTFCs of the multimodal network, Figure 3 demonstrates that the shipper and the receiver need to prepare more budget for a more costly multimodal route to deal with the increasing capacity uncertainty of the multimodal network, so that their containers can be transported successfully from origin to destination.
Additionally, we change the left hand spread ratio from 5% to 30% with a step of 5% and increase the confidence level from 0.5 to 1.0 with a step of 0.1. Then, we use the proposed model to generate the optimal routes under different combinations of the spread ratio and the confidence level. The results are illustrated in Figure 4. It should be noted that there is no feasible solution to the MRPPDTWCU when the left hand ratio is 30% and the confidence level is 1.0.
Figure 4 verifies the sensitivities shown in Figure 2 and Figure 3. Furthermore, Figure 4 shows that the degree of increase in total costs of the optimal route caused by improving the confidence level is reduced when the left hand spread ratio becomes smaller, which is clearly illustrated by Figure 5. In Figure 5, the gap is calculated by Equation (29) where c o s t 1.0 and c o s t 0.5 are the total costs of the optimal route when the confidence level is 0.5 and 1.0, respectively.
g a p = c o s t 1.0 c o s t 0.5 c o s t 0.5 100 %
Figure 4 and Figure 5 indicate that if a more reliable multimodal route is needed and meanwhile the shipper and the receiver would like to reduce the budget prepared for the multimodal transportation, the MTO needs to select the transportation services and transfer services with smaller uncertainty levels of fuzzy capacities, i.e., the transportation services and transfer services with stable capacities.

6. Conclusions

This study explores a MRPPDTWCU. The modeling of pickup and delivery time windows aims to improve the timeliness of the entire multimodal transportation process from its beginning to its end. The capacity uncertainty is modeled by LRTFNs, which enables the uncertainty level of the fuzzy capacity of the multimodal network to be measurable. Under the above consideration, we propose a fuzzy nonlinear optimization model to address the MRPPDTWCU and design an equivalent CCLP reformulation based on the FCCP method and linearization processing. A numerical case verifies that the proposed model yields high efficiency in obtaining the global optimum solution to the MRPPDTWCU. The findings of the numerical case using the proposed model are:
(1)
The modeling of pickup and delivery time windows is able to provide both the shipper and the receiver with on-time transportation services. Such a consideration also enables the shipper, receiver, and MTO to make a balance between the travel and transfer costs and the storage costs to realize the cost minimization.
(2)
Compared to the deterministic modeling, the MRP under capacity uncertainty enables the shipper, receiver, and MTO to make more flexible decisions, in which they can make tradeoffs between the economy and reliability of transportation.
(3)
The capacity uncertainty shows a significant influence on the MRPPDTWCU from two aspects. The first is the confidence level that is introduced into the MRPPDTWCU by the FCCP and reflects the reliability of transportation, and the second is the uncertainty level of the fuzzy capacities.
(4)
Improving the confidence level to achieve a reliable multimodal route scarifies the transportation economy. Therefore, the shipper and the receiver should make tradeoffs between the reliability and economic objectives by determining a suitable confidence level, in which the sensitivity shown in Figure 2 can provide a solid reference. Then, the MTO can plan the optimal route based on their demand using the proposed model.
(5)
To address the higher uncertainty level of the fuzzy capacity, the shipper and the receiver need to increase their transportation budget. To help them reduce the budget and meanwhile maintain a high confidence level to ensure reliable transportation, the MTO needs to use transportation services and transfer services with stable capacities.
In practical transportation, the shipper and the receiver first negotiate to determine their respective time windows and the transportation volume of the containers, and then propose their transportation order to the MTO. After receiving the transportation order, the MTO establishes the multimodal network for container transportation preferentially using the transportation services and transfer services with relatively stable capacities. The sensitivity analysis can help the shipper and the receiver to determine the confidence level based on their attitudes toward the reliability of transportation and the budget that they can afford for multimodal transportation. After setting the above issues, the MTO can use the proposed CCLP model to find the optimal route that yields the minimum costs to move the containers from the shipper to the receiver, in which the combination of economy and timeliness of multimodal transportation gets optimized.

Author Contributions

Conceptualization, J.G. and Y.S.; data curation, J.G. and Y.S.; formal analysis, J.G. and Y.S.; funding acquisition, Y.S.; investigation, Y.S.; methodology, J.G. and Y.S.; resources, Y.S.; software, J.G. and Y.S.; supervision, Y.S.; validation, J.G. and Y.S.; writing—original draft, J.G. and Y.S.; writing—review and editing, Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Shandong Provincial Natural Science Foundation of China under Grant No. ZR2023MG020 and the University–Industry Collaborative Education Program of the Ministry of Education of China under Grant Nos. 230711255407320 and 220601369102256.

Data Availability Statement

All data used in this study are included in this paper and can be obtained from the cited reference [41].

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wang, Q.Z.; Chen, J.M.; Tseng, M.L.; Luan, H.M.; Ali, M.H. Modelling green multimodal transport route performance with witness simulation software. J. Clean. Prod. 2020, 248, 119245. [Google Scholar] [CrossRef]
  2. Sun, Y.; Sun, G.; Huang, B.; Ge, J. Modeling a carbon-efficient road–rail intermodal routing problem with soft time windows in a time-dependent and fuzzy environment by chance-constrained programming. Systems 2023, 11, 403. [Google Scholar] [CrossRef]
  3. Macharis, C.; Bontekoning, Y.M. Opportunities for OR in intermodal freight transport research: A review. Eur. J. Oper. Res. 2004, 153, 400–416. [Google Scholar] [CrossRef]
  4. Caris, A.; Macharis, C.; Janssens, G.K. Decision support in intermodal transport: A new research agenda. Comput. Ind. 2013, 64, 105–112. [Google Scholar] [CrossRef]
  5. Shahedi, A.; Gallo, F.; Saeednia, M.; Sacco, N. Lead-time-based freight routing in multi-modal networks considering the Physical Internet. J. Supply Chain Manag. Sci. 2023, 4, 61–80. [Google Scholar] [CrossRef]
  6. Verma, M.; Verter, V. A lead-time based approach for planning rail-truck intermodal transportation of dangerous goods. Eur. J. Oper. Res. 2010, 202, 696–706. [Google Scholar] [CrossRef]
  7. Dua, A.; Sinha, D. Quality of multimodal freight transportation: A systematic literature review. World Rev. Intermodal Transp. Res. 2019, 8, 167–194. [Google Scholar] [CrossRef]
  8. Sun, Y.; Hrušovský, M.; Zhang, C.; Lang, M. A time-dependent fuzzy programming approach for the green multimodal routing problem with rail service capacity uncertainty and road traffic congestion. Complexity 2018, 2018, 8645793. [Google Scholar] [CrossRef]
  9. Zhang, D.; He, R.; Li, S.; Wang, Z. A multimodal logistics service network design with time windows and environmental concerns. PLoS ONE 2017, 12, e0185001. [Google Scholar] [CrossRef] [PubMed]
  10. Zhang, X.; Jin, F.Y.; Yuan, X.M.; Zhang, H.Y. Low-carbon multimodal transportation path optimization under dual uncertainty of demand and time. Sustainability 2021, 13, 8180. [Google Scholar] [CrossRef]
  11. Fazayeli, S.; Eydi, A.; Kamalabadi, I.N. Location-routing problem in multimodal transportation network with time windows and fuzzy demands: Presenting a two-part genetic algorithm. Comput. Ind. Eng. 2018, 119, 233–246. [Google Scholar] [CrossRef]
  12. Yuan, X.; Jiang, Y.; Zhang, X. Research on robust optimization of interval-based fuzzy intermodal transport paths under low-carbon policies. Ind. Eng. Manag. 2021, 26, 134–141. [Google Scholar]
  13. Li, X.; Sun, Y.; Qi, J.; Wang, D. Chance-constrained optimization for a green multimodal routing problem with soft time window under twofold uncertainty. Axioms 2024, 13, 200. [Google Scholar] [CrossRef]
  14. Tang, J.; Pan, Z.; Fung, R.Y.; Lau, H. Vehicle routing problem with fuzzy time windows. Fuzzy Sets Syst. 2009, 160, 683–695. [Google Scholar] [CrossRef]
  15. Dragomir, A.G.; Doerner, K.F. Solution techniques for the inter-modal pickup and delivery problem in two regions. Comput. Oper. Res. 2020, 113, 104808. [Google Scholar] [CrossRef]
  16. Sun, Y. Fuzzy approaches and simulation-based reliability modeling to solve a road–rail intermodal routing problem with soft delivery time windows when demand and capacity are uncertain. Int. J. Fuzzy Syst. 2020, 22, 2119–2148. [Google Scholar] [CrossRef]
  17. Qi, Y.; Harrod, S.; Psaraftis, H.N.; Lang, M. Transport service selection and routing with carbon emissions and inventory costs consideration in the context of the Belt and Road Initiative. Transp. Res. Part E Logist. Transp. Rev. 2022, 159, 102630. [Google Scholar] [CrossRef]
  18. Zhang, Y.; Guo, W.; Negenborn, R.R.; Atasoy, B. Synchromodal transport planning with flexible services: Mathematical model and heuristic algorithm. Transp. Res. Part C Emerg. Technol. 2022, 140, 103711. [Google Scholar] [CrossRef]
  19. Delbart, T.; Molenbruch, Y.; Braekers, K.; Caris, A. Uncertainty in intermodal and synchromodal transport: Review and future research directions. Sustainability 2021, 13, 3980. [Google Scholar] [CrossRef]
  20. Lu, K.; Liao, H.; Zavadskas, E.K. An overview of fuzzy techniques in supply chain management: Bibliometrics, methodologies, applications and future directions. Technol. Econ. Dev. Econ. 2021, 27, 402–458. [Google Scholar] [CrossRef]
  21. Hosseini, A.; Pishvaee, M.S. Capacity reliability under uncertainty in transportation networks: An optimization framework and stability assessment methodology. Fuzzy Optim. Decis. Mak. 2022, 21, 479–512. [Google Scholar] [CrossRef]
  22. Huang, Y. A two-stage stochastic model for intermodal transportation operational planning problem under capacity shortage. In Proceedings of the 2021 IEEE International Conference on Recent Advances in Systems Science and Engineering (RASSE), Shanghai, China, 12–14 December 2021; pp. 1–9. [Google Scholar]
  23. Uddin, M.; Huynh, N. Reliable routing of road-rail intermodal freight under uncertainty. Netw. Spat. Econ. 2019, 19, 929–952. [Google Scholar] [CrossRef]
  24. Uddin, M.; Huynh, N.N.; Ahmed, F. Assignment of freight traffic in a large-scale intermodal network under uncertainty. Highlights Sustain. 2024, 3, 1–15. [Google Scholar] [CrossRef]
  25. Zimmermann, H.J. Applications of fuzzy set theory to mathematical programming. In Readings in Fuzzy Sets for Intelligent Systems; Elsevier: Amsterdam, The Netherlands, 1993; pp. 795–809. [Google Scholar]
  26. Lu, Y.; Lang, M.; Sun, Y.; Li, S. A fuzzy intercontinental road-rail multimodal routing model with time and train capacity uncertainty and fuzzy programming approaches. IEEE Access 2020, 8, 27532–27548. [Google Scholar] [CrossRef]
  27. Sun, Y.; Yu, N.; Huang, B. Green road–rail intermodal routing problem with improved pickup and delivery services integrating truck departure time planning under uncertainty: An interactive fuzzy programming approach. Complex Intell. Syst. 2022, 8, 1459–1486. [Google Scholar] [CrossRef]
  28. Zarandi, M.H.F.; Hemmati, A.; Davari, S. The multi-depot capacitated location-routing problem with fuzzy travel times. Expert Syst. Appl. 2011, 38, 10075–10084. [Google Scholar] [CrossRef]
  29. Sakawa, M.; Nishizaki, I.; Uemura, Y. Interactive fuzzy programming for multi-level linear programming problems with fuzzy parameters. Fuzzy Sets Syst. 2000, 109, 3–19. [Google Scholar] [CrossRef]
  30. Kundu, P.; Kar, S.; Maiti, M. Multi-objective multi-item solid transportation problem in fuzzy environment. Appl. Math. Model. 2013, 37, 2028–2038. [Google Scholar] [CrossRef]
  31. Bierwirth, C.; Kirschstein, T.; Meisel, F. On transport service selection in intermodal rail/road distribution networks. Bus. Res. 2012, 5, 198–219. [Google Scholar] [CrossRef]
  32. Wang, Z.X.; Liu, Y.J.; Fan, Z.P.; Feng, B. Ranking L–R fuzzy number based on deviation degree. Inf. Sci. 2009, 179, 2070–2077. [Google Scholar] [CrossRef]
  33. Zhou, J.; Yang, F.; Wang, K. Fuzzy arithmetic on LR fuzzy numbers with applications to fuzzy programming. J. Intell. Fuzzy Syst. 2016, 30, 71–87. [Google Scholar] [CrossRef]
  34. Compaoré, A.; Somé, K.; Poda, J.; Somé, B. Efficiency of MOMA-Plus method to solve some fully fuzzy LR triangular multiobjective linear programs. J. Math. Res. 2018, 10, 77–87. [Google Scholar] [CrossRef]
  35. Gani, A.N.; Duraisamy, C.; Veeramani, C. A note on fuzzy linear programming problem using LR fuzzy number. Int. J. Algorithms Comput. Math. 2009, 2, 93–106. [Google Scholar]
  36. Xu, J.; Zhou, X. Approximation based fuzzy multi-objective models with expected objectives and chance constraints: Application to earth-rock work allocation. Inf. Sci. 2013, 238, 75–95. [Google Scholar] [CrossRef]
  37. Zheng, Y.; Liu, B. Fuzzy vehicle routing model with credibility measure and its hybrid intelligent algorithm. Appl. Math. Comput. 2006, 176, 673–683. [Google Scholar] [CrossRef]
  38. Zhu, H.; Zhang, J. A credibility-based fuzzy programming model for APP problem. In Proceedings of the 2009 International Conference on Artificial Intelligence and Computational Intelligence, Shanghai, China, 7–8 November 2009; Volume 1, pp. 455–459. [Google Scholar]
  39. Peykani, P.; Hosseinzadeh Lotfi, F.; Sadjadi, S.J.; Ebrahimnejad, A.; Mohammadi, E. Fuzzy chance-constrained data envelopment analysis: A structured literature review, current trends, and future directions. Fuzzy Optim. Decis. Mak. 2022, 21, 197–261. [Google Scholar] [CrossRef]
  40. Roshanaei, V.; Luong, C.; Aleman, D.M.; Urbach, D.R. Reformulation, linearization, and decomposition techniques for balanced distributed operating room scheduling. Omega 2020, 93, 102043. [Google Scholar] [CrossRef]
  41. Sun, Y.; Lang, M. Bi-objective optimization for multi-modal transportation routing planning problem based on Pareto optimality. J. Ind. Eng. Manag. 2015, 8, 1195–1217. [Google Scholar] [CrossRef]
  42. Xiong, G.; Wang, Y. Best routes selection in multimodal networks using multi-objective genetic algorithm. J. Comb. Optim. 2014, 28, 655–673. [Google Scholar] [CrossRef]
Figure 1. Multimodal network consisting of 35 nodes and 136 arcs [13,41,42].
Figure 1. Multimodal network consisting of 35 nodes and 136 arcs [13,41,42].
Axioms 13 00220 g001
Figure 2. Sensitivity of the MRPPDTWCU with respect to the confidence level.
Figure 2. Sensitivity of the MRPPDTWCU with respect to the confidence level.
Axioms 13 00220 g002
Figure 3. Sensitivity of the MRPPDTWCU with respect to the left hand spread ratio.
Figure 3. Sensitivity of the MRPPDTWCU with respect to the left hand spread ratio.
Axioms 13 00220 g003
Figure 4. Sensitivity of the MRPPDTWCU with respect to the combination of the left hand spread ratio and the confidence level.
Figure 4. Sensitivity of the MRPPDTWCU with respect to the combination of the left hand spread ratio and the confidence level.
Axioms 13 00220 g004
Figure 5. Gaps between the total costs of the optimal routes under confidence levels of 0.5 and 1.0.
Figure 5. Gaps between the total costs of the optimal routes under confidence levels of 0.5 and 1.0.
Axioms 13 00220 g005
Table 1. Travel costs and time of the transportation modes [13].
Table 1. Travel costs and time of the transportation modes [13].
Transportation ModesTravel Costs (CNY/TEU)Travel Time (h)
Rail c i j m = 500 + 2.03 d i j m t i j m = d i j m / 60
Road c i j m = 15 + 8 d i j m t i j m = d i j m / 80
Water c i j m = 950 t i j m = d i j m / 30
Table 2. Transfer time and costs between different transportation modes [13,41].
Table 2. Transfer time and costs between different transportation modes [13,41].
Transfer TypesTransfer Time (h/TEU)Transfer Costs (CNY/TEU)
Rail~Road0.0675
Rail~Water0.1337
Road~Water0.10010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ge, J.; Sun, Y. Solving a Multimodal Routing Problem with Pickup and Delivery Time Windows under LR Triangular Fuzzy Capacity Constraints. Axioms 2024, 13, 220. https://doi.org/10.3390/axioms13040220

AMA Style

Ge J, Sun Y. Solving a Multimodal Routing Problem with Pickup and Delivery Time Windows under LR Triangular Fuzzy Capacity Constraints. Axioms. 2024; 13(4):220. https://doi.org/10.3390/axioms13040220

Chicago/Turabian Style

Ge, Jie, and Yan Sun. 2024. "Solving a Multimodal Routing Problem with Pickup and Delivery Time Windows under LR Triangular Fuzzy Capacity Constraints" Axioms 13, no. 4: 220. https://doi.org/10.3390/axioms13040220

APA Style

Ge, J., & Sun, Y. (2024). Solving a Multimodal Routing Problem with Pickup and Delivery Time Windows under LR Triangular Fuzzy Capacity Constraints. Axioms, 13(4), 220. https://doi.org/10.3390/axioms13040220

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