Next Article in Journal
Designing Our Own Board Games in the Playful Space: Improving High School Student’s Citizenship Competencies and Creativity through Game-Based Learning
Next Article in Special Issue
A Hybrid Multi-Criteria Decision-Making Model Combining DANP with VIKOR for Sustainable Supplier Selection in Electronics Industry
Previous Article in Journal
Mechanical and Microstructural Properties of a Stabilized Sand Using Geopolymer Made of Wastes and a Natural Pozzolan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Incorporating Vehicle-Routing Problems into a Closed-Loop Supply Chain Network Using a Mixed-Integer Linear-Programming Model

1
Independent Researcher, North Vancouver, BC V7P 1T7, Canada
2
Department of Business Administration, University of Gothenburg, SE-405 30 Gothenburg, Sweden
3
Prime School of Logistics, Saito University College, Petaling Jaya 46200, Selangor, Malaysia
4
Apadana Institute of Higher Education, Shiraz 7187985443, Iran
5
Independent Researcher, Addis Ababa P.O. Box 1164, Ethiopia
6
Department of Industrial Engineering, Semnan Branch, Islamic Azad University, Semnan 3381774895, Iran
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(4), 2967; https://doi.org/10.3390/su15042967
Submission received: 18 October 2022 / Revised: 2 January 2023 / Accepted: 2 February 2023 / Published: 6 February 2023

Abstract

:
In recent years, there has been a tremendous increase in environmental awareness, due to concerns about sustainability. Designing an efficient supply chain network that fulfills the expectation of both business owners and customers and, at the same time, pays attention to environmental protection is becoming a trend in the commercial world. This study proposes a theoretical model incorporating vehicle routing problems (VRPs) into the typical CLSC (closed-loop supply chain) network architecture. This combination assists all operators to act more efficiently in terms of environmental protection and profitability. A mixed-integer-linear-programming model for CLSC network design with fuzzy and random uncertain data is developed to achieve the goals. The parameters of the CLSC network are also programmed using hybrid fuzzy-stochastic mathematical programming. The model is for a single product and a single timeframe. Several numerical examples are provided to demonstrate the validity of the proposed mixed-integer-linear-programming (MILP) model. This study also investigated probabilistic possibilities for recourse variables with a trapezoidal fuzzy number using a problem size of four cases. The result indicates that the model performed well in the numerical test, suggesting it can help the operation to be more profitable if this model is implemented in their daily routines.

1. Introduction

The Paris Agreement has highlighted the significance of preserving natural resources and reducing fossil fuel consumption [1]. Although the forum’s main objective was reducing global warming, global warming is the direct product of excessive consumption of fossil fuel, which is caused by the unnecessary presence of vehicles in the supply chain network. Another side of the inattentive consumption of fossil fuels that accompanies the destruction of forests and wildlife is transmitting viruses hidden in the natural environment to our living places, which may cause several pandemics with shorter time spans in the future [2] .Moreover, destroying biodiversity can put human health at risk and even put our food security in danger [3], and consequently, the human race is on the precipice of extinction. The global economy is also facing many challenges, such as global imbalances, economic exclusion, prospects for growth, deregulated markets, inflation, energy and the environment, inequality, labor issues, emerging markets, and the impact of new technologies. Policymakers set common sustainability goals, which explain what should be developed and what is to be sustained, for how long, and for the benefit of whom [4]. However, implementing sustainability concepts in the global society is challenging, since it will confront man’s current lifestyle, which has been exercised for centuries [5]. With all the challenges, sustainability has been institutionalized to develop rules and regulations to re-structure social intervention, create awareness, and shape behaviors globally. One of the main challenges of sustained economic growth is using a linear production model [6].
In the linear production model, resources are fed to the production line and transported to manufacturers to be processed into some parts and products. After some time, the finalized products will be discarded. This routine will be repeated frequently, increasing resource prices because of resource depletion. Therefore, the goal of sustainable development goal is the long-term stability of natural resources, the environment, societal concerns, and the economy, acknowledging sustainable development through the decision-making process of policymakers. In this regard, the key policy adopted by companies is to take more responsibility in collecting used products at the end of their life cycle and provide recovery alternatives such as remanufacturing, recycling, and the disposal of used products in an environmentally friendly manner. Switching from the linear model, which mainly focuses on the concept of take-make-use-and-dispose, to a circular economy necessitates changes in which companies generate value, understand, and do business. This transformation requires a new business model for sustainability and circularity. Moreover, these changes may lead managers to ambiguous situations in designing and managing the distribution and collection channels for new and used products. Therefore, developing an efficient network that supports the management of both forward and reverse flows is a challenge. In this situation, a closed-loop supply chain network is a sustainable approach.
A closed-loop supply chain is the integration of a forward supply chain (the processes of converting raw material to finished product) and a reverse supply chain (the collection and recovery of the products returned in the supply chain). Due to interdependent decisions in forward and reverse supply chain networks, considering the networks separately could lead to sub-optimal results [7], and therefore it is important to take into account both forward and reverse networks simultaneously or to create a network configuration in the decision-making process.
In designing a CLSC, it is crucial to integrate the operational, strategic, and tactical decisions, ensuring a better business perspective concerning reducing costs and improving customer service levels [8]. A recent literature review by Oliveira and Machado [9] showed that 48% of studies are involved in strategic decision making. Moreover, the research focusing on the tactical decision amounted to 72 papers, representing 23.5% of the reviewed papers. Only three studies have considered strategic, tactical, and operational decision making in their model. One of the critical operational decisions is vehicle-routing management, which finds an efficient route for transferring products through the network [8]. Although identifying an optimal vehicle route has the potential to reduce costs and improve customer-service levels significantly, research studies in the field of CLSC have mostly overlooked the importance of this operational decision [10]. Another limitation of the CLSC model that has yet to be addressed in the literature of CLSC is the simultaneous delivery of products and pick-up of end-of-life products; almost all operators plan and perform different fleets for picking up malfunctioned or end-of-life products. Moreover, some manufacturers prefer to outsource their reverse supply chain or assign them to retailers or third parties [11,12,13], since they think recycling the used products is costly and is not profitable for their firms. This study develops a new CLSC model with simultaneous delivery and pick-up integrated with location-allocation decisions implying uncertainty. Different types of vehicles, such as trucks, planes, and ships, typically carry out the movement of materials through the supply chain network. They all operate based on precise and empirical programs, which are typically costly for operators and business owners. With this consideration, logistics planning comes forward to facilitate the management of the operation. Vehicle routing programs (VRPs) are part of these procedures that assign each vehicle to each customer on the nodes on the graphs usually designed for logistics planning [14].
Moreover, uncertainties such as the growth in the complexity of structures, market fluctuation, and service facilities make the planning more intricate. Several approaches have been suggested, for instance, stochastic programming [15], robust optimization [16], and fuzzy programming [17]. Researchers have attempted to combine these methods to address uncertainties more efficiently [18]. In order to align this study with the above research stream, this paper proposes a fuzzy-stochastic programming approach to overcome the uncertainties associated with the problem. The general form of an optimization problem is to find the maximum (or minimum) of a particular “objective function” under some “constraints”; the solution to this problem is known as the optimization method. In classical optimizations, the objective function and constraints are deterministic (i.e., not fuzzy). Nevertheless, in the real world, the values are hardly found in deterministic form and typically fluctuate; it needs to be decided which criteria should be implemented in a system or eliminated, so solving uncertain optimization problems is essential for both theory and application. Fuzzy mathematical programming is considered, based on the concept of a fuzzy variable that is one of the complete patterns of fuzziness. This study proposes a new mixed-integer linear mathematical-programming model to minimize the total costs, including transportation, investment, and operational costs.

2. Literature Review

The concept of CLSC has emerged to reconcile the environmental and economic objectives of the supply chain. As a result, many studies can be found in the literature that have attempted to integrate these two aspects into a research frame [19]. In this regard, Govindan and Soleimani [20,21], as well as Govindan, Soleimani [8] conducted a literature review that summarized the growing interest in the closed-loop supply chain and also analyzed the environmental, legal and social, as well as the economic factors involved in CLSC. Quattrociocchi et al. [22] provided a comprehensive overview of the tourism supply chain, which is multi-billion-dollar business; any flaw in this era brings about dire consequences for business owners, consumers, and the environment. Moreover, identifying risk and calculating its effects in an operation can be very helpful for mitigating the side effects of any business practices. Mital et al. [23] formulated a model by implementing the analytic hierarchy process for risk identification and assessing the impacts on business performance. Since this research paper focuses on the latest modeling of CLSC, and the integration of CLSC with the vehicle routing problem, the literature related to this area is briefly reviewed in two streams, namely the CLSC network and VRP, in order to make a comparison with the earlier research works in this field.
CLSC mainly addresses economic-performance improvement through value recovery from end-of-life products [24]. Ghomi-Avili et al. [25] formulated a closed-loop supply chain recruiting a fuzzy bi-objective bi-level function to maintain the supply network at an optimum level in times of disruption when the whole network comes to arbitrary stagnation. Jabbarzadeh et al. [26] fashioned a model utilizing stochastic-robust-optimization methods to guarantee the resiliency and performance of a closed-loop supply chain when an interruption occurs. Their model was mainly focused on locating the facilities and enhancing transportation to minimize the cost of the supply network. Wu and Zhou [27] developed a modeled underlying game theory that analyzed the renewal of manufacturing equipment, conducted by third-party manufacturers, and its economic effects on different players throughout the supply chain. Zheng et al. [12] developed game models and numerical studies analyzing how to share profits among players in CLSC equally and fairly, and maintain their overall satisfaction. Zhen et al. [28] proposed a green closed-loop supply chain with a balance between the normal operations of the supply chain and the operation of environmentally friendly procedures. The authors recruited the Lagrangean relaxation method to solve the problems. Wang et al. [11] formulated a scenario-based model to scrutinize different problems in a CLSC environment where competition among players, from manufacturers to retailers, prevailed in an uncontrolled manner. The authors provided mathematical models to examine different scenarios between manufacturers and recycling firms, to study what happens if a company totally or partially outsources its recycling operations. Liu et al. [29] examined the impact of product design on the profitability of companies. The model revealed that performing CLSC was not profitable for a corporation, but redesigning the products intensified and guaranteed a firm’s profitability. Mohtashami et al. [30] presented a comprehensive CLSC model focusing on the logistic portion of the supply network to analyze how fleet-network optimization could provide golden opportunities to minimize the negative environmental impacts typically created by vehicles. The authors provided numerical examples using a queuing system to solve the problems. Govindan et al. [31] provided a model to inspect supply networks from various dimensions, from supplier selection to routing problems. They utilized many mathematical methods for solving problems, such as the fuzzy-analytic-network process, fuzzy-decision-making trial and evaluation laboratory, and multi-objective mixed-integer linear programming. Yavari and Zaker [32] formulated a green closed-loop supply chain mainly focusing on perishable products and the disruption of the supply network, probably instigated by the power outage. The authors also offered several strategies for eliminating the impacts of the problems. Fathollahi-Fard et al. [33] provided a multi-objective stochastic optimization model to tackle the problem of potable water supply and sewerage collection within an uncertain environment. The authors utilized a multi-objective stochastic optimization model to solve the problems. The chief goals of the research were to offer a methodology to curb water-wastage problems, save valuable material as much as possible, and offer some solutions for implementing the strategy in the real world. Chen et al. [13] developed a model to facilitate and rationalize the procedures of selecting providers who conduct reverse supply chains for manufacturers. The critical point for conducting this research was assisting people willing to outsource their circular-economy agenda. Szmelter-Jarosz et al. [34] formulated a convenient model using the neutrosophic-fuzzy method to handle the supply chain in the era of the COVID-19 pandemic, in order to eliminate vehicle congestion and deliver cargo to consumers in an efficient manner. The research aimed to maintain the service level and follow the environmentally friendly operating principles in an era when the public stress level was soaring. Shabbir et al. [35] integrated resiliency, sustainability, and reliability into the supply chain network, following CLSC principles. The authors employed the Lp-metric and Lagrangian relaxation methods to solve the problems. Liao et al. [36] focused on agricultural products, considered to be among those perishable materials, to develop an appropriate supply chain with regard to environmentally friendly principles. Harvesting fresh products is a business that produces excessive pollution if business owners do not follow agricultural laws and regulations. This study provided a comprehensive study on organizing and operating in this era. The researchers utilized the genetic algorithm and simulated annealing to solve the problems.
Amin and Zhang [37] developed a single-product, single-period model for a CLSC network that helps remanufactured products to be sent to a secondary market. Shi, Nie [38] studied a CLSC which utilized third-party-logistics service providers for the forward supply chain. Along the same line of thought, Sasikumar and Haq [39] developed a CLSC network in which a third-party reverse logistics provider (3PRLP) was responsible for the reverse supply chain. Das and Rao Posinasetti [40] addressed environmental concerns arising from harmful emissions and energy consumption in the CLSC-network design. In addition to the economic objective, several researchers incorporate social and environmental objectives, as interest in environmental protection, customer satisfaction, and sustainable development is increasing [24]. Garg, Kannan [41] developed a bi-objective mixed-integer linear-programming model for CLSC-network design to maximize total profits and minimize the number of hired vehicles in the forward supply chain. Taleizadeh, Haghighi [42] developed a multi-objective CLSC to maximize profits, minimize environmental effects, and minimize social objectives. Hasanov, Jaber [43] addressed the effects of inventory policy on CLSC costs by considering emissions from production and transportation. The results showed a higher rate for the collection of used items, improving the environmental performance and reducing its cost. Mawandiya, Jha [44] developed a two-echelon closed-loop supply chain to find the optimal-lot sizing and shipment policies. Giri and Masanta [45] developed a closed-loop supply chain model where the production process is subject to learning and forgetting.
In CLSS-network development, many parameters must be defined precisely, due to the associated uncertainties [46]. In the literature, different approaches have been specified by researchers in handling uncertainty, namely fuzzy optimization [47], stochastic optimization [48], and robust optimization [49]. Soleimani, Seyyed-Esfahani [50] studied a location-allocation problem in a CLSC-network design. They integrated three types of risk measures: mean absolute deviation, value at risk, and conditional value at risk (CVaR) into the two-stage stochastic- programming model. Their profit-oriented model revealed that risk-neutral approaches are not efficient. Later, Subulan, Baykasoğlu [51] developed a mixed-integer linear-programming model by assuming financial and collection risk with almost identical risk measurements as that of Soleimani, Seyyed-Esfahani [50]. Next, Dai and Zheng [52] developed a multi-period, multi-product CLSC-network design addressing uncertainties in the disposing of rates, demand, and capacities. They used a Monte Carlo simulation along with hybrid GA/fuzzy programming and chance-constrained programming. Khatami, Mahootchi [53] developed a mixed-integer linear-programming model for CLSC. The main concerns of their model were the uncertainties about product demand and return. Bender’s decomposition was applied to overcome uncertainties. Accordingly, Radhi and Zhang [54] developed a mixed-integer non-linear-programming model by taking into account the uncertainty of quality and demand. Keyvanshokooh, Ryan [55] proposed a hybrid robust-stochastic programming method to address uncertainties about the demand, return, and transportation cost. In another study, Ma, Yao [48] developed a robust multi-objective mixed-integer nonlinear-programming model for the CLSC-network design. This model focused on the environmental impact of the network as well as its cost. Saedinia, Vahdani [56] developed a bi-objective in the oil and gas industry. Jeihoonian, Kazemi Zanjani [57] developed a two-stage stochastic for the CLSC model, where the quality of return products was uncertain. The authors of [58] developed a two-stage stochastic non-linear programming model to minimize total costs. The model was solved by a heuristic-tabu search algorithm.
Physical distribution is one of the most prominent and costly functions of any logistics system, and requires the transportation of products from the manufacturer through the distribution center to the customers. The vehicle routing problem (VRP) is a generic name that refers to optimization problems, in that customers are served using a number of vehicles [59]. Nagy and Salhi [60] extended the basic VRP model by including customers who may receive and send goods. Ahmadi Javid and Azad [61] developed a stochastic supply chain model with location-allocation, an optimum amount of capacity and inventory, and an optimum decision for selecting the best route. Nekooghadirli, Tavakkoli-Moghaddam [62] established a bi-objective location-routing-inventory (LRI) model that considered a multi-period and multi-product system which is used to make strategic and tactical decisions for location and inventory-routing. Bae and Moon [63] developed a multi-depot vehicle-routing problem with time windows to study the delivery and installation of electronic products. They showed the possibility for cost minimization of the depot, delivery, and installation of the vehicle, as well as travel distance and labor. Archetti, Desaulniers [64] developed an inventory-routing problem that minimizes a logistic ratio, which is the ratio of the total routing cost to the total quantity distributed. Iassinovskaia, Limbourg [65] developed an inventory-routing problem where the returnable transport items can be collected at the customer’s location. Wang, Shao [66] developed a two-echelon capacitated-vehicle-routing problem with environmental considerations. Ahmadi-Javid, Amiri [67] developed a model for location-routing pricing problems, which considers profit maximization. Madankumar and Rajendran [68] considered a special case of VRP that addresses the routing problem in a semiconductor supply chain. They developed an MILP (mixed-integer linear-programming) model for solving the green VRP with pickup and deliveries in a semiconductor supply chain. Abu Al Hla, Othman [69] solved a VRP model that improves driver satisfaction, customers’ perceived quality, and the company’s financial objective. They considered three levels of strategic, tactical, and operational decisions. The result shows that giving autonomy to drivers will not involve a significant cost, and the model can be optimum. Madankumar and Rajendran [70] developed an MILP model for VRP with simultaneously delivery and pickup, and considered time windows to improve the performance and responsiveness of the model. Similarly, Zhou, Qin [71] developed VRP with time windows and simultaneous pickup and delivery. Their network consisted of two echelons. They formulated the model and solved it with a neighborhood tabu-search algorithm. The result shows that the algorithm that they used can save 19% computational time, on average.
To better recognize the novelty of the current article compared with the previously published papers, a summary of the previous works is illustrated in Table 1.
The discussion mentioned above highlighted the research gaps which lack CLSC and VRP integration and are inadequate in addressing uncertainty in the CLSC and the simultaneous delivery and taking of end-of-life materials from customers, which has already been discussed in the introduction. Two researchers, Garg, Kannan [41] and Zhalechian, Tavakkoli-Moghaddam [10] tried to integrate CLSC and VRP; however, Garg, Kannan [41] developed a CLSC model that optimized total profit and the number of hired vehicles in a forward supply chain only, and did not include the reverse supply chain. Zhalechian, Tavakkoli-Moghaddam [10] developed a model to optimize the inventory model of CLSC. None addressed the CLSC-network design and VRP concurrently, a research gap that this paper attempts to cover. Integration of CLSC-network development with VRP supports the company’s decision-making by encompassing all three types of strategic, tactical, and operational decisions in parallel. The strategic decision is addressed in selecting the location, the tactical decision is applied by determining the flow of materials and products among the facilities, and the operational decision is conveyed by defining a route to distribute products to the customers.
To address the gaps mentioned above, adopting an integrated fuzzy-mathematical and stochastic-programming approach, such as a hybrid fuzzy-stochastic one, is a suitable choice for addressing uncertainties in CLSC-network design. Therefore, this study developed a mixed-integer linear-programming model for CLSC-network design with routing decisions under fuzzy and random uncertain data.

3. Model Description and Formulation

As stated previously, CLSC is complex by nature, and thus our proposed network provides hybrid facilities to reduce the total costs and complexity of the network. The hybrid facility (HF) at the proposed CLSC network (Figure 1) coordinates a distribution and collection center. The application of HF is of economic benefit to the network because, in general, the co-location of facilities reduces considerable investment in human resources, equipment, and infrastructure [80].
As shown in Figure 1, OEMs and customers are the existing sites in the network, while HFs are potential new facilities. The OEMs purchase raw materials and components from suppliers (S) and recycling centers (RC). The OEM has special vehicles assigned to deliver products to HFs. These vehicles are used to distribute and collect end-of-life products returned from customers within assigned routes. The challenge in this phase is selecting a route that minimizes total costs while ensuring that each customer receives his demand. Each HF, in addition to its collection activity, is responsible for inspecting, testing, and sorting the used products. For simplicity, inspection, testing, and sorting will be subsequently referred to as sorting. The end-of-life products collected from the HF are sent to the recycling centers to be recycled and reused as materials for production. Customers are grouped in zones called aggregate delivery points [81].
There are two main goals for developing a CLSC model:
  • To select the optimum CLSC-network configuration.
  • Finding the routes for the limited number of vehicles to serve a group of customers with their demand.
Every customer has an uncertain delivery demand, D ˜ r , and a pick-up, P ˜ r , for new and end-of-life products simultaneously. Moreover, the following assumptions are postulated:
  • Facility locations are known beforehand.
  • The following of products between HFs is not allowed.
  • The set-up cost of facilities is considered as fixed and predefined.

3.1. Objective Function

The proposed mixed-integer linear-programming model aims to minimize the total cost, including the establishment cost of HF, purchasing cost from the supplier, production cost, purchasing cost from the recycling center, and transportation costs.
The fixed opening cost includes a decision on the establishment of HFs that can be written as:
h H F 1 h U 1 h
The processing cost involves both purchasing and production costs. The purchasing cost includes costs related to the purchasing cost of raw materials from recycling centers and suppliers, and the production cost includes the sum of setup cost, labor, and all costs of production, which can be written as follows:
s S i I m M P S s i X 0 s i m
m M h H v i V I P M m X 1 m h v i
c C i I m M P C c i X 4 c i m
The transportation cost includes the traveling of a vehicle along its defined route between each node. It involves traveling from HSP to HF, from HF to customers, and from HF to the recycling center, which can be written as:
m M h H v i V I d 1 m h T C 1 m h v i X 1 m h v i
j b j b k K d 2 j j T C 2 j j k X 2 j j k
h H c C v j V J d 3 h c T C 3 h c v j X 3 h c v j
Although the main objective is to minimize cost, revenue could be achieved from selling the end-of-life products to the recycling center, which can be written as:
h H c C v j V J S P X 3 h c v j
Therefore:
M i n Z = h H F 1 h U 1 h + s S i I m M P S s i X 0 s i m + m M h H v i V I P M m X 1 m h v i h H c C v j V J S P X 3 h c v j + c C i I m M P C c i X 4 c i m + m M h H v i V I d 1 m h T C 1 m h v i X 1 m h v i + j b j b k K d 2 j j T C 2 j j k X 2 j j k + h H c C v j V J d 3 h c T C 3 h c v j X 3 h c v j

3.2. Constraints

The constraints of the model are structured as follows:
h H v i V I X 1 m h v i C M m m M (10)
R 2 h C D h U 1 h h H (11)
h H v j V J X 3 h c v j C C c c C (12)
m M X 4 c i m C C c U i c C ; i I (13)
m M h H X 1 m h v i C 1 v i v i V I (14)
L 2 k C 2 k k K (15)
L t j j k C 2 k j , j b ; k K (16)
L t j j k C 2 k j , j b ; k K (17)
h H c C X 3 h c v j C 3 v j v j V J (18)
h H X 4 c i m + s S X 0 s i m = h H v i V I U i X 1 m h v i m M ; i I (19)
h H v i V I U i X 3 h c v j = m M X 4 c i m c C ; i I (20)
h H v j V J X 3 h c v j = R 2 h h h (21)
k K j b X 2 j j k = 1 j R (22)
j b X 2 j j k = j b X 2 j j k j R ; k K (23)
j R X 2 j j k = j R X 2 j j k k K ; j H (24)
j H j R X 2 j j k 1 k K (25)
r R Y 2 h r U 1 h h H (26)
M 2 r k M 2 r k + B 2 X 2 j j k B 2 1 r , r R ; k K (27)
X 2 j j k = 0 k K ; j , j H (28)
j b j b d 2 j j X 2 j j k M D 2 k k K (29)
j b X 2 j j k + j b X 2 i j k Y 2 h r 1 j R ; i H ; k K (30)
h H Y 2 h r = 1 r R (31)
m M v i V I X 1 m h v i R 2 h h H (32)
L 2 k j b j R D ˜ r X 2 j j k k K (33)
R 2 h r R D ˜ r y 2 h r h H (34)
L t j j k X 2 j j k r R D r j , j b ; k K (35)
k K b B L t b r k k K b B L t r b k + D r r R (36)
L 2 k j b j R P ˜ r X 2 j j k k K (37)
R 2 h r R P ˜ r Y 2 h r j H (38)
L t j j k X 2 j j k r R P r j , j b ; k K (39)
k K b B L t b r k k K b B L t r b k P r r R (40)
Constraints (10)–(18) are capacity constraints that ensure the vehicle and facilities meet their capacities. Constraints (10)–(13) ensure production centers, hybrid facilities, and recycling centers cannot exceed their capacity. Constraint (14) is capacity requirements for a vehicle that guarantee the number of products carried with vehicle vj meets the capacity. Constraints (15)–(17) are related to capacity vehicle k, and guarantee that vehicle k cannot load more than its capacity. Constraint (18) is associated with vehicle vj, and ensures vehicle vj meets its capacity. Constraints (19)–(21) are flow-balance constraints at production centers, hybrid facilities, and recycling centers, which specify the relationship between input and output at respective facilities. Constraints (22)–(31) are associated with routing. Constraint (22) ensures every customer belongs to one, and only one, route. Constraint (23) ensures every customer belongs to one, and only one, HF. Constraint (24) indicates that each customer arrives and leaves with the same vehicle. Constraints (25) and (26) secure one HF in each route. Constraint (27) is a sub-tour elimination. Constraint (28) guarantees that there is no route between HFs. Constraint (29) limits the maximum distance of vehicle type k. Constraint (30) is the linkage and allocation constraints that ensure that if vehicle k starts its trip from HF j and serves customer r during its trip, then customer r should be assigned to HF. Constraint (31) assures each customer is assigned to one HF. Constraint (32) assures vehicle load after visiting all HFs. Constraint (33) indicates vehicle load for customers’ demand. Constraint (34) satisfies the vehicle load after visiting all customers. Constraints (35) and (36) ensure that the vehicle load is not exceeding the demand quantity. Constraint (37) indicates expected return products for HFs. Constraint (38) assures the vehicle load for the pick-up product after visiting all HFs. Constraints (39) and (40) ensure that the vehicle pick-up load does not exceed the pick-up quantity.

3.3. Proposed Solution

In the stochastic-programming model, uncertainty is characterized by the random nature of the parameters. The objective function and the constraints are developed using fuzzy numbers: this transforms the problem from deterministic to stochastic. The fuzzy-stochastic type of uncertainty can be formulated using two-stage stochastic programming (TSP) with a recourse model. In this method, the decisions are made in two phases:
  • Parameters need to be determined before the value of a random variable.
  • Parameters need to be determined after the random event has happened.
The second-stage decision can be made for the minimization of the penalty that occurs because of infeasibility [82]. Therefore, the TSP model can be formulated as follows:
M a x     f = c x E [ Q ( x , ω ) ]
s . t .   A x b
x 0
where x is the first-stage anticipated decisions made before random variables are observed and Q ( x , ξ ) is the optimum value, for any given Ω , of the following nonlinear program:
                                  m i n                 q ( y , ω )
                                                                            s . t .   W ( ω ) y = h ( ω ) T ( ω ) x
y 0
where y is the second-stage decision variable that depends on the realization of the first-stage random vector; q ( y , ω ) denotes the second-stage cost functions; T ( ω ) , W ( ω ) , h ( ω ) | ω Ω are model parameters with reasonable dimensions, and are a function of the random vector, ( ω ) . For given values of the first-stage variables, ( x ) , the second-stage problem can be decomposed into independent, linear sub-problems, with one sub-problem for each realization of the uncertain parameters. Then, model (41)–(43) can be reformulated as:
M a x       f = c s E [ min q y 0 ( y , ω ) | T ( ω ) x + W ( ω ) y = h ( ω ) ]
  s . t .   A x b
x 0
The above TSP problem is nonlinear [83]; however, the problem can be converted to a linear form by assuming discrete distributions for the uncertain parameters. Thus, the expected value of the nonlinear term E [ Q ( x , ω ) ] can be linearized as follows [83]:
E [ Q ( x , ω ) ] = h = 1 s p h Q ( x , ω h )
where ω possesses a discrete and finite distribution, with support Ω = ω 1 , ω 2 , , ω s ; p h represents the probability of realization of scenario ω h , with p h 0 and h = 1 s p h = 1 .
Therefore, the model (41)–(43) can be converted to a linear form as in the following:
M a x   f = c x h = 1 s p h q ( y h , ω h )
s . t .   A x b
T ( ω h ) x + W ( ω h ) y ω h = h ( ω h ) ,   ω h Ω
x 0
y h 0
The model (51)–(55) can be formulated with chance-constrained programming (CCP), as follows:
M a x   f = j = 1 n 1 c j x j j = 1 n 2 h = 1 s p h d j y j h
s . t .                 j = 1 n 1 a r j x j b r                         r = 1 , 2 , , m 1
Pr j = 1 n 1 a i j x j + j = 1 n 2 a i j y j h w i h γ i                     i = 1 , 2 , , m 2 h = 1 , 2 , , s
x j 0                   j = 1 , 2 , , n 1
y j h 0                     j = 1 , 2 , , n 2 h = 1 , 2 , , s
where γ i ( γ i [ 0 , 1 ] ) is the probability level that constraint i should be satisfied with at least.
In addition to the lack of historical data or the high cost of data acquisition, errors in obtained information, variations in a spatial and temporal unit, and incomplete or imprecise observed information [84,85], make it challenging to determine the probability distribution of uncertain parameters. Fuzzy mathematical programming provides a tool to describe uncertain parameters mathematically. Fuzzy possibilistic programming (FPP), among others, is a suitable method that enables decision-makers to model uncertainty, particularly in cases where the uncertainty lies either in the coefficients of the left-hand side and right-hand side of constraints or coefficients of the objective function [83]. The FPP represents the parameter in a fuzzy form that contains possibility distributions. Therefore, a hybrid fuzzy-stochastic programming model can be formulated as follows [83]:
M a x   f = j = 1 n 1 c ˜ j x j j = 1 n 2 h = 1 s p h d ˜ j y j h
s . t .               j = 1 n 1 a ˜ r j x j b ˜ r               r = 1 , 2 , , n 1
j = 1 n 1 a ˜ i j x j + j = 1 n 2 a ˜ i j y j h w ˜ i h               i = 1 , 2 , , m 2 h = 1 , 2 , , s
x j 0               j = 1 , 2 , , n 1
y j h 0               j = 1 , 2 , , n 2 h = 1 , 2 , , s
where x j and y j h are first and second-stage decision variables, respectively; c j and d j are fuzzy coefficients in the objective function; a ˜ i j and a ˜ r j are fuzzy left-hand side coefficients, and w ˜ i h is the independent random variables with a known probability distribution.
This hybrid fuzzy-stochastic model can be solved by using chance constraint programming (CCP). CCP can be applied when some right-hand-side parameters are uncertain, with a known probability distribution. The CCP approach converts the model to a deterministic form by applying a fixed, certain level of probability, q ˜ i [ 0 , 1 ] for every uncertain constraint, i, and imposing the condition that the constraint is satisfied with at least a probability level of γ ˜ i = 1 q ˜ i with γ ˜ i [ 0 , 1 ] . Thus, the fuzzy-stochastic constraint can be described as [83]:
Pr { j = 1 n 1 a ˜ i j x j + j = 1 n 2 a ˜ i j y j h w ˜ i h } γ ˜ i               i = 1 , 2 , , m 2 h = 1 , 2 , , s
Generally, constraint (66) is nonlinear, and the set of feasible constraints is convex only for some particular distribution and certain level of q ˜ i . Therefore, this constraint can be formulated as the equivalent fuzzy deterministic form [86], as follows:
j = 1 n 1 a ˜ i j x j + j = 1 n 2 a ˜ i j y j h w ˜ i h q ˜ i               i = 1 , 2 , , m 2 h = 1 , 2 , , s
where w ˜ q ˜ i i h = F i 1 ( q ˜ i ) , given the cumulative distribution function of w ˜ i h and the probability of violating constraint i ( q ˜ i ) . Therefore, the model can be written as follows:
M a x   f = j = 1 n 1 c ˜ j x j j = 1 n 2 h = 1 s p h d ˜ j y j h
s . t .                 j = 1 n 1 a ˜ r j x j b ˜ r                 r = 1 , 2 , , m 1
j = 1 n 1 a ˜ i j x j + j = 1 n 2 a ˜ i j y j h w ˜ i h q ˜ i                 i = 1 , 2 , , m 2 h = 1 , 2 , , s
x j 0                 j = 1 , 2 , , n 1
y j h 0                 j = 1 , 2 , , n 2 h = 1 , 2 , , s
The model can be solved with the help of a fuzzy set and α c u t approaches. The possibility distribution of fuzzy parameters b ˜ r = ( b ¯     r , b r 1 , b r 2 , b ¯ r ) can be characterized as a trapezoidal fuzzy set when b r 1 b r 2 , and can be specified as a triangular fuzzy set when b r 1 = b r 2 . Parameter a ˜ i j under each α c u t level can be formulated as a closed crisp interval: [ ( 1 α ) b ¯     r + α b r 1 , ( 1 α ) b r 2 + α b ¯ ) ] . Kaufmann, Gil Aluja [87] established 11 levels for feasibility degree starting from an unacceptable solution (α = 0) up to a completely acceptable solution (α = 1). The model can be written as a form of two deterministic forms with the range of lower and upper bounds of the objective-function value.
The lower bound is as follows:
f l = j = 1 n 1 [ ( 1 α ) c ¯     j + α c j 1 ] x j j = 1 n 2 h = 1 s p h [ ( 1 α ) d ¯ j + α d j 2 ] y j h
s . t .           j = 1 n 1 [ ( 1 α ) a ¯ r j + α a r j 2 ] x j [ ( 1 α ) b ¯     r + α b r 1 ]                 r = 1 , 2 , , m 1
j = 1 n 1 [ ( 1 α ) a ¯ i j + α a i j 2 ] x j + j = 1 n 2 [ ( 1 α ) a _   i j + α a i j 1 ] y j h w i h [ ( 1 α ) q ¯ i + α q i 1 ]                 i = 1 , 2 , , m 2 h = 1 , 2 , , s
x j 0                 j = 1 , 2 , , n 1
y j h 0                 j = 1 , 2 , , n 2 h = 1 , 2 , , s
The upper bound is as follows:
f u = j = 1 n 1 [ ( 1 α ) c ¯ j + α c j 2 ] x j j = 1 n 2 h = 1 s p h [ ( 1 α ) d ¯ j + α d j 1 ] y j h
j = 1 n 1 [ ( 1 α ) a ¯ r j + α a r j 1 ] x j [ ( 1 α ) b ¯     r + α b r 2 ]                 r = 1 , 2 , , m 1
j = 1 n 1 [ ( 1 α ) a ¯ i j + α a i j 1 ] x j + j = 1 n 2 [ ( 1 α ) a ¯ i j + α a i j 2 ] y j h w i h [ ( 1 α ) q ¯ i + α q i 2 ]                 i = 1 , 2 , , m 2 h = 1 , 2 , , s
x j 0                 j = 1 , 2 , , n 1
y j h 0                 j = 1 , 2 , , n 2 h = 1 , 2 , , s
Then, after solving the lower and upper form of the model with various α c u t , a set of optimum solutions related to the hybrid fuzzy-stochastic model can be obtained as the following:
x j o p t = [ x j o p t l , x j o p t u ]                 j
y j h o p t = [ y j h o p t l , y j h o p t u ]                 j , h
f o p t = [ f o p t l , f o p t u ]

3.4. Numerical Study

Solving the model with the proposed method requires identifying first-stage and second-stage variables. First-stage variables related to the opening of the facilities are the first step, and variables related to allocation are recourse decision variables. In addition, it requires identifying fuzzy parameters and writing models according to their relevant sub-models’ lower-bound model (63)–(67) and upper-bound model (68)–(72). To illustrate the validity of the proposed model and the usefulness of the proposed solution approach, several numerical experiments are implemented. This also includes numerical examples that demonstrate how an increase in the size of the problem, which results in model complexity, affects the model and solving time. Table 2 shows the size of the examples with four different problem sizes. In this paper, we considered two scenarios for recourse variables with probability. A value of parameters of the model is assumed to have a trapezoid fuzzy number, and is shown in Table 3. The distance between the manufacturer and the hybrid facility is a uniform distribution between 60 and 450. The distance between hybrid facilities and demand zones and between demand zones and hybrid facilities are uniform distributions between 60 and 500. The utilization value is a uniform distribution between 1 and 4. Six α-cut levels (0, 0.2, 0.4, 0.6, 0.8, 1) were examined for the objective function. The MILP model was written in GAMS optimization software and solved with ILOG CPLEX 15.5.1 It was discovered that the size of the model has a direct relationship with its running time; as the size of the model increases, so does the time required to solve the mode. For example, problem number 2 for α = 0.4 took 0.42 s (CPU time) to solve and problem number 1 for α = 0.4 took 0.25 s (CPU time) to solve. The value of the objective functions for all the problem numbers are presented in Table 4.
Figure 2 shows the network structure of the model when α- cut is 0.4 for the lower bound for problem number one, and the model took 0.25 s (CPU time) to be solved. It shows that the model selected suppliers 2 and 3 to send raw materials to the manufacturer. It also shows that hybrid facilities 1, 5, and 6 were selected to support the demand zones. As it can be seen, hybrid facilities 1 supports demand zone 6, 4, and 8, consecutively.
Figure 3 depicts the model’s network structure when the cut is set to 0.6 for the lower bound on problem number one. Figure 2 and Figure 3 show that retailer number 5 is assigned to hybrid facility number 6, and each supplier is dedicated to only one manufacturer. Furthermore, retailer 7 has been added to the route of retailers 1 and 2 for hybrid facility 5. When these two figures are compared, the number of arrows decreases, indicating that the model is attempting to optimize routes and assign hybrid functions to retailers as efficiently as possible.
Figure 4 depicts a network structure of the model when the α- cut is 0.4 for the lower bound for problem number two and the model took 0.42 s (CPU time) to be solved. It shows suppliers 2,3 and 5 were selected from among five suppliers. It also shows that hybrid facilities 1, 4, 7, 9, and 11 were selected to support demand zones. Figure 3 and Figure 4 also demonstrate how changes in the size of the model will affect the complexity of the model and the network.
Figure 5 shows the value of the objective function for problem number 2 for the upper and lower bounds. As can be seen, the objective-function value changes with the change in α-cut levels. It shows that with the increase in α-cut levels, the value of the objective function decreases for the upper bound. On the other hand, with the increase in α-cut levels, the value of the objective function increases for the lower bound of the objective function. Moreover, it shows that the acceptable range of the objective function has a reverse relationship with the α-cut levels: an increase in the α-cut-level objective value has a narrower range, and vice versa. This can be explained with the Kaufmann, Gil Aluja [87] rules of feasibility degree. As previously stated, they established 11 levels of feasibility degree, ranging from completely unacceptable (=0) to completely acceptable (=1). When α = 0.1 the level of confidence is very low; therefore, the model generates a high range of value for the objective function. As the level of confidence grows, the model becomes more certain, and the objective function value has a narrower range.
Figure 6 shows the value of objective functions when the number of return products increased by 20%. It shows that with the increase in α- cut levels, the range of acceptable objective-function values decreases.

4. Conclusions

Because of the growing concern about sustainability and the future, there has been a rise in recent years in the industry-wide awareness of environmental protection [88,89,90]; despite its price, CLSC is a solution for sustainability [91]. This paper explained that it is possible to achieve significant cost savings and revenue growth while also protecting the environment by designing an efficient supply chain network that considers both forward and reverse networks during the planning process. Moreover, to accomplish this, the authors developed a model that incorporates the vehicle-routing problem into the conventional CLSC-network design and a strategy for dealing with uncertain parameters that arise during the development of the CLSC network. In order to achieve the objectives, a mixed-integer linear-programming model for CLSC-network design with routing decisions under fuzzy and random uncertain data was developed. A hybrid fuzzy-stochastic mathematical-programming approach was employed to overcome the complexity and uncertainty in the parameters of the CLSC network.
Nevertheless, Kabak and D. Ruan [92] proposed a mathematical model with resource allocation and outsourcing decisions. They were only two decision variables in fuzzy form; in this study, all the variables are in fuzzy numbers to assist decision makers to make accurate decisions regarding cost reduction, increasing the revenue of a company as well as environmental protection. The objective function and the constraints are developed using fuzzy numbers, which transforms the problem from deterministic to stochastic. The fuzzy-stochastic type of uncertainty can be formulated using two-stage stochastic programming (TSP) with a recourse model. In this method, the decisions are made in two phases; phase 1 parameters need to be determined before the value of a random variable, and phase 2 parameters need to be determined after the random event has happened. The second-stage decision can be made for the minimization of penalties that occur because of infeasibility. Fallahpour et al. [93] proposed a hybrid fuzzy-programming approach for the modeling of a sustainable, resilient supply chain network. A case study of the palm oil industry in Malaysia was proposed. Their study focuses on the forward supply chain only. The advantage of this current research is that in this research the authors developed a CLSC-network design by integrating VRP. Integrating CLSC with VRP furnishes companies with a useful decision-support system encompassing all three types of strategic, tactical, and operational decisions in parallel. The strategic decision is addressed in selecting the location, the tactical decision is applied by determining the flow of materials and products among the facilities, and the operational decision is conveyed by defining a route to distribute products to the customers.
A series of different examples were formed to demonstrate the validity of the proposed MILP model. In this paper, we considered a variety of scenarios. The model performed well, and it can be used with a single product and for a single period of time, which shows great promise for contributing to the advancement of both knowledge and practice. Every study has limits, regardless of how carefully it is performed, and this study is no exception. When developing the models, there were several limitations, most of which were connected to the models’ assumptions. As with most mathematical models which contribute to basic research outcomes, the results of this study are highly dependent on the assumptions that were made in formulating the problem. Future studies can broaden the scope of the study to include a broader range of products and time periods; another area of investigation in the near term is the recognition of a model using real-world data.

Author Contributions

Conceptualization, A.P.; methodology, A.P.; writing—original draft preparation, A.P., F.M.; writing—review and editing A.A.; S.S. 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

The data will be availabe per contacting th efirst author, A.P.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

The following notations have been used in this paper:
sset of suppliers s S
iset of raw materials i I
mset of plants m M
hset of hybrid facilities ( H F s ) h H
r&jset of customers
baggregate set of hybrid facilities and customers ( h r )
cset of recycling centers c C
viset of vehicles v i V I at node M
kset of vehicles k K
vjset of vehicles v j V J at node RC
Parameters:
d 1 m h distance between plants m M and hybrid facility h H
d 2 j j distance between customer in node j J and j J
d 3 h c distance between hybrid facility h H and recycling center c C
D ˜ r demand of customer r R
P ˜ r pickup customer r R
M D 2 k maximum distance which vehicle k K covers in a tour
C 1 v i capacity of vehicle v i V I
C 2 k capacity of vehicle k K
C 3 v j capacity of vehicle v j V J
C S s i capacity of supplier s S to supply raw material i I
C M m production capacity of plant m M
C D h distribution capacity of hybrid facility h H
C C c recycling capacity of recycling center c C
P S s i purchasing cost per unit of material i I from supplier s S
P M m production cost per unit of product at plant m M
P R c recycling cost per unit of product at recycling center c C
P C c i purchasing cost of unit of material i I from recycling center c C
B 2 total number of customers
F 1 h fixed opening cost of hybrid facility h H
T C 1 m h v i transportation cost per unit of product of vehicle v i V I between plant m M and hybrid facility h H
T C 2 j j k transportation cost per unit of product of vehicle k K between customer j J and customer j J
T C 3 h c v j transportation cost per unit of product of vehicle v j V J between hybrid facility h H and recycling center c C
u i utilization rate of material i I
S P selling price of end-of-life product
Decision variables:
X 2 j j k binary variable indicating whether vehicle k K travels directly from node   j J to node j J
Y 2 h r binary variable if hybrid facility h H is assigned to customer r R
U 1 h binary   variable   if   hybrid   facility   h H is open
L 2 k load of vehicle k K when leaving hybrid facility
L 2 k load of vehicle k K after having serviced all assigned customers
M 2 r k sub-tour elimination variable for customer r R   in   k K
R 2 h distribution quantity of hybrid facility h H
R 2 h pick-up quantity for return product of hybrid facility h H
X 0 s i m shipment quantity of raw material i I between supplier s S and plant m M
X 1 m h v i shipment quantity between plant m M and hybrid facility h H with vehicle v i V I
X 3 h c v j shipment quantity between hybrid facility h H and recycling center c C with vehicle v j V J
X 4 c i m shipment quantity of raw material i I between recycling center c C and plant m M
L t j j k unload of demand for vehicle k K from node j J to node j J
L t j j k load of vehicle k K from node j J to node j J

References

  1. United Nations Climate Change. 2015. Available online: https://unfccc.int/process-and-meetings/the-paris-agreement/the-paris-agreement (accessed on 1 December 2022).
  2. Carlson, C.J.; Albery, G.F.; Merow, C.; Trisos, C.H.; Zipfel, C.M.; Eskew, E.A.; Olival, K.J.; Ross, N.; Bansal, S. Climate change increases cross-species viral transmission risk. Nature 2022, 607, 555–562. [Google Scholar] [CrossRef]
  3. World Health Organization. Biodiversity and Health. 2015. Available online: https://www.who.int/news-room/fact-sheets/detail/biodiversity-and-health (accessed on 1 December 2022).
  4. Acero, L.; Savaget, P. Plural understandings of sociotechnical progress within the OECD. In Proceedings of the 12th Globelics International Conference, Addis Ababa, Ethiopia, 29–31 October 2014. [Google Scholar]
  5. Hák, T.; Janoušková, S.; Moldan, B. Sustainable Development Goals: A need for relevant indicators. Ecol. Indic. 2016, 60, 565–573. [Google Scholar] [CrossRef]
  6. Bonviu, F. The European economy: From a linear to a circular economy. Rom. J. Eur. Aff. 2014, 14, 78. [Google Scholar]
  7. Pishvaee, M.S.; Farahani, R.Z.; Dullaert, W. A memetic algorithm for bi-objective integrated forward/reverse logistics network design. Comput. Oper. Res. 2010, 37, 1100–1112. [Google Scholar] [CrossRef]
  8. Govindan, K.; Soleimani, H.; Kannan, D. Reverse logistics and closed-loop supply chain: A comprehensive review to explore the future. Eur. J. Oper. Res. 2015, 240, 603–626. [Google Scholar] [CrossRef]
  9. Oliveira, L.S.; Machado, R.L. Application of optimization methods in the closed-loop supply chain: A literature review. J. Comb. Optim. 2021, 41, 357–400. [Google Scholar] [CrossRef]
  10. Zhalechian, M.; Tavakkoli-Moghaddam, R.; Zahiri, B.; Mohammadi, M. Sustainable design of a closed-loop loca-tion-routing-inventory supply chain network under mixed uncertainty. Transp. Res. Part E Logist. Transp. Rev. 2016, 89, 182–214. [Google Scholar] [CrossRef]
  11. Wang, N.; He, Q.; Jiang, B. Hybrid closed-loop supply chains with competition in recycling and product markets. Int. J. Prod. Econ. 2019, 217, 246–258. [Google Scholar] [CrossRef]
  12. Zheng, X.-X.; Li, D.-F.; Liu, Z.; Jia, F.; Sheu, J.-B. Coordinating a closed-loop supply chain with fairness concerns through varia-ble-weighted Shapley values. Transp. Res. Part E Logist. Transp. Rev. 2019, 126, 227–253. [Google Scholar] [CrossRef]
  13. Chen, Z.-S.; Zhang, X.; Govindan, K.; Wang, X.-J.; Chin, K.-S. Third-party reverse logistics provider selection: A computa-tional semantic analysis-based multi-perspective multi-attribute decision-making approach. Expert Syst. Appl. 2021, 166, 114051. [Google Scholar] [CrossRef]
  14. Roboredo, M.; Sadykov, R.; Uchoa, E. Solving vehicle routing problems with intermediate stops using VRPSolver models. Networks 2022. [Google Scholar] [CrossRef]
  15. Zeballos, L.J.; Mendez, C.; Barbosa-Póvoa, A.P.; Novais, A.Q. Multi-stage stochastic optimization of the design and planning of a Closed-Loop Supply Chain. Comput. Aided Chem. Eng. 2013, 32, 691–696. [Google Scholar] [CrossRef]
  16. Pishvaee, M.S.; Rabbani, M.; Torabi, S. A robust optimization approach to closed-loop supply chain network design under un-certainty. Appl. Math. Model. 2011, 35, 637–649. [Google Scholar] [CrossRef]
  17. Khalili-Damghani, K.; Tavana, M.; Amirkhan, M. A fuzzy bi-objective mixed-integer programming method for solving supply chain network design problems under ambiguous and vague conditions. Int. J. Adv. Manuf. Technol. 2014, 73, 1567–1595. [Google Scholar] [CrossRef]
  18. Talaei, M.; Moghaddam, B.F.; Pishvaee, M.S.; Bozorgi-Amiri, A.; Gholamnejad, S. A robust fuzzy optimization model for carbon-efficient closed-loop supply chain network design problem: A numerical illustration in electronics industry. J. Clean. Prod. 2016, 113, 662–673. [Google Scholar] [CrossRef]
  19. Kazemi, N.; Modak, N.M.; Govindan, K. A review of reverse logistics and closed loop supply chain management studies published in IJPR: A bibliometric and content analysis. Int. J. Prod. Res. 2018, 57, 4937–4960. [Google Scholar] [CrossRef]
  20. Govindan, K.; Soleimani, H. A review of reverse logistics and closed-loop supply chains: A Journal of Cleaner Production focus. J. Clean. Prod. 2017, 142, 371–384. [Google Scholar] [CrossRef]
  21. Cammarano, A.; Perano, M.; Michelino, F.; Del Regno, C.; Caputo, M. SDG-Oriented Supply Chains: Business Practices for Procurement and Distribution. Sustainability 2022, 14, 1325. [Google Scholar] [CrossRef]
  22. Quattrociocchi, B.; Mercuri, F.; Perano, M.; Calabrese, M. Tourism supply chain & strategic partnerships for managing the com-plexity in tourism industry. Enl. Tour. A Pathmaking J. 2017, 7, 62–93. [Google Scholar]
  23. Mital, M.; Del Giudice, M.; Papa, A. Comparing supply chain risks for multiple product categories with cognitive mapping and Analytic Hierarchy Process. Technol. Forecast. Soc. Chang. 2018, 131, 159–170. [Google Scholar] [CrossRef]
  24. Yu, H.; Solvang, W.D. A fuzzy-stochastic multi-objective model for sustainable planning of a closed-loop supply chain considering mixed uncertainty and network flexibility. J. Clean. Prod. 2020, 266, 121702. [Google Scholar] [CrossRef]
  25. Ghomi-Avili, M.; Naeini, S.G.J.; Tavakkoli-Moghaddam, R.; Jabbarzadeh, A. A fuzzy pricing model for a green competitive closed-loop supply chain network design in the presence of disruptions. J. Clean. Prod. 2018, 188, 425–442. [Google Scholar] [CrossRef]
  26. Jabbarzadeh, A.; Haughton, M.; Khosrojerdi, A. Closed-loop supply chain network design under disruption risks: A robust approach with real world application. Comput. Ind. Eng. 2018, 116, 178–191. [Google Scholar] [CrossRef]
  27. Wu, X.; Zhou, Y. Buyer-specific versus uniform pricing in a closed-loop supply chain with third-party remanufacturing. Eur. J. Oper. Res. 2018, 273, 548–560. [Google Scholar] [CrossRef]
  28. Zhen, L.; Huang, L.; Wang, W. Green and sustainable closed-loop supply chain network design under uncertainty. J. Clean. Prod. 2019, 227, 1195–1209. [Google Scholar] [CrossRef]
  29. Liu, Z.; Li, K.W.; Li, B.-Y.; Huang, J.; Tang, J. Impact of product-design strategies on the operations of a closed-loop supply chain. Transp. Res. Part E Logist. Transp. Rev. 2019, 124, 75–91. [Google Scholar] [CrossRef]
  30. Mohtashami, Z.; Aghsami, A.; Jolai, F. A green closed loop supply chain design using queuing system for reducing en-vironmental impact and energy consumption. J. Clean. Prod. 2020, 242, 118452. [Google Scholar] [CrossRef]
  31. Govindan, K.; Mina, H.; Esmaeili, A.; Gholami-Zanjani, S.M. An Integrated Hybrid Approach for Circular supplier selection and Closed loop Supply Chain Network Design under Uncertainty. J. Clean. Prod. 2019, 242, 118317. [Google Scholar] [CrossRef]
  32. Yavari, M.; Zaker, H. Designing a resilient-green closed loop supply chain network for perishable products by considering disruption in both supply chain and power networks. Comput. Chem. Eng. 2019, 134, 106680. [Google Scholar] [CrossRef]
  33. Fathollahi-Fard, A.M.; Ahmadi, A.; Al-E-Hashem, S.M. Sustainable closed-loop supply chain network for an integrated water supply and wastewater collection system under uncertainty. J. Environ. Manag. 2020, 275, 111277. [Google Scholar] [CrossRef]
  34. Szmelter-Jarosz, A.; Ghahremani-Nahr, J.; Nozari, H. A Neutrosophic Fuzzy Optimisation Model for Optimal Sustainable Closed-Loop Supply Chain Network during COVID-19. J. Risk Financial Manag. 2021, 14, 519. [Google Scholar] [CrossRef]
  35. Shabbir, M.S.; Mahmood, A.; Setiawan, R.; Nasirin, C.; Rusdiyanto, R.; Gazali, G.; Arshad, M.A.; Khan, S.; Batool, F. Closed-loop Supply Chain Network Design with Sustainability and Resiliency Criteria. Ph.D. Thesis, Petra Christian University, Surabaya, Indonesia, 2021. [Google Scholar]
  36. Liao, Y.; Kaviyani-Charati, M.; Hajiaghaei-Keshteli, M.; Diabat, A. Designing a closed-loop supply chain network for citrus fruits crates considering environmental and economic issues. J. Manuf. Syst. 2020, 55, 199–220. [Google Scholar] [CrossRef]
  37. Amin, S.H.; Zhang, G. A proposed mathematical model for closed-loop network configuration based on product life cycle. Int. J. Adv. Manuf. Technol. 2011, 58, 791–801. [Google Scholar] [CrossRef]
  38. Shi, Y.; Nie, J.; Qu, T.; Chu, L.-K.; Sculli, D. Choosing reverse channels under collection responsibility sharing in a closed-loop supply chain with re-manufacturing. J. Intell. Manuf. 2013, 26, 387–402. [Google Scholar] [CrossRef]
  39. Sasikumar, P.; Haq, A.N. Integration of closed loop distribution supply chain network and 3PRLP selection for the case of battery recycling. Int. J. Prod. Res. 2011, 49, 3363–3385. [Google Scholar] [CrossRef]
  40. Das, K.; Posinasetti, N.R. Addressing environmental concerns in closed loop supply chain design and planning. Int. J. Prod. Econ. 2015, 163, 34–47. [Google Scholar] [CrossRef]
  41. Garg, K.; Kannan, D.; Diabat, A.; Jha, P. A multi-criteria optimization approach to manage environmental issues in closed loop supply chain network design. J. Clean. Prod. 2015, 100, 297–314. [Google Scholar] [CrossRef]
  42. Taleizadeh, A.A.; Haghighi, F.; Niaki, S.T.A. Modeling and solving a sustainable closed loop supply chain problem with pricing decisions and discounts on returned products. J. Clean. Prod. 2018, 207, 163–181. [Google Scholar] [CrossRef]
  43. Hasanov, P.; Jaber, M.; Tahirov, N. Four-level closed loop supply chain with remanufacturing. Appl. Math. Model. 2018, 66, 141–155. [Google Scholar] [CrossRef]
  44. Mawandiya, B.K.; Jha, J.; Thakkar, J. Optimal production-inventory policy for closed-loop supply chain with reman-ufacturing under random demand and return. Oper. Res. 2018, 20, 1623–1664. [Google Scholar]
  45. Giri, B.C.; Masanta, M. A closed-loop supply chain model with uncertain return and learning-forgetting effect in pro-duction under consignment stock policy. Oper. Res. 2022, 22, 947–975. [Google Scholar]
  46. Wang, H.-F.; Hsu, H.-W. Resolution of an uncertain closed-loop logistics model: An application to fuzzy linear programs with risk analysis. J. Environ. Manag. 2010, 91, 2148–2162. [Google Scholar] [CrossRef]
  47. Zceylan, E. Simultaneous optimization of closed- and open-loop supply chain networks with common components. J. Manuf. Syst. 2016, 41, 143–156. [Google Scholar] [CrossRef]
  48. Ma, R.; Yao, L.; Jin, M.; Ren, P.; Lv, Z. Robust environmental closed-loop supply chain design under uncertainty. Chaos Solitons Fractals 2016, 89, 195–202. [Google Scholar] [CrossRef]
  49. Kisomi, M.S.; Solimanpur, M.; Doniavi, A. An integrated supply chain configuration model and procurement management under uncertainty: A set-based robust optimization methodology. Appl. Math. Model. 2016, 40, 7928–7947. [Google Scholar] [CrossRef]
  50. Soleimani, H.; Seyyed-Esfahani, M.; Kannan, G. Incorporating risk measures in closed-loop supply chain network design. Int. J. Prod. Res. 2013, 52, 1843–1867. [Google Scholar] [CrossRef]
  51. Subulan, K.; Baykasoğlu, A.; Özsoydan, F.; Taşan, A.; Selim, H. A case-oriented approach to a lead/acid battery closed-loop supply chain network design under risk and uncertainty. J. Manuf. Syst. 2015, 37, 340–361. [Google Scholar] [CrossRef]
  52. Dai, Z.; Zheng, X. Design of close-loop supply chain network under uncertainty using hybrid genetic algorithm: A fuzzy and chance-constrained programming model. Comput. Ind. Eng. 2015, 88, 444–457. [Google Scholar] [CrossRef]
  53. Khatami, M.; Mahootchi, M.; Farahani, R.Z. Benders’ decomposition for concurrent redesign of forward and closed-loop supply chain network with demand and return uncertainties. Transp. Res. Part E Logist. Transp. Rev. 2015, 79, 1–21. [Google Scholar] [CrossRef]
  54. Radhi, M.; Zhang, G. Optimal configuration of remanufacturing supply network with return quality decision. Int. J. Prod. Res. 2015, 54, 1487–1502. [Google Scholar] [CrossRef]
  55. Keyvanshokooh, E.; Ryan, S.M.; Kabir, E. Hybrid robust and stochastic optimization for closed-loop supply chain network design using accelerated Benders decomposition. Eur. J. Oper. Res. 2016, 249, 76–92. [Google Scholar] [CrossRef]
  56. Saedinia, R.; Vahdani, B.; Etebari, F.; Nadjafi, B.A. Robust gasoline closed loop supply chain design with redistricting, service sharing and intra-district service transfer. Transp. Res. Part E Logist. Transp. Rev. 2019, 123, 121–141. [Google Scholar] [CrossRef]
  57. Jeihoonian, M.; Zanjani, M.K.; Gendreau, M. Closed-loop supply chain network design under uncertain quality status: Case of durable products. Int. J. Prod. Econ. 2017, 183, 470–486. [Google Scholar] [CrossRef]
  58. Zhen, L.; Sun, Q.; Wang, K.; Zhang, X. Facility location and scale optimisation in closed-loop supply chain. Int. J. Prod. Res. 2019, 57, 7567–7585. [Google Scholar] [CrossRef]
  59. Montoya-Torres, J.R.; Franco, J.L.; Isaza, S.N.; Jiménez, H.F.; Herazo-Padilla, N. A literature review on the vehicle routing problem with multiple depots. Comput. Ind. Eng. 2015, 79, 115–129. [Google Scholar] [CrossRef]
  60. Nagy, G.; Salhi, S.D. Heuristic algorithms for single and multiple depot vehicle routing problems with pickups and de-liveries. Eur. J. Oper. Res. 2005, 162, 126–141. [Google Scholar] [CrossRef]
  61. Ahmadi Javid, A.; Azad, N. Incorporating location, routing and inventory decisions in supply chain network design. Transp. Res. Part E Logist. Transp. Rev. 2010, 46, 582–597. [Google Scholar] [CrossRef]
  62. Nekooghadirli, N.; Tavakkoli-Moghaddam, R.; Ghezavati, V.; Javanmard, S. Solving a new bi-objective loca-tion-routing-inventory problem in a distribution network by meta-heuristics. Comput. Ind. Eng. 2014, 76, 204–221. [Google Scholar] [CrossRef]
  63. Bae, H.; Moon, I. Multi-depot vehicle routing problem with time windows considering delivery and installation vehicles. Appl. Math. Model. 2016, 40, 6536–6549. [Google Scholar] [CrossRef]
  64. Archetti, C.; Desaulniers, G.; Speranza, M.G. Minimizing the logistic ratio in the inventory routing problem. EURO J. Transp. Logist. 2017, 6, 289–306. [Google Scholar] [CrossRef]
  65. Iassinovskaia, G.; Limbourg, S.; Riane, F. The inventory-routing problem of returnable transport items with time windows and simultaneous pickup and delivery in closed-loop supply chains. Int. J. Prod. Econ. 2017, 183, 570–582. [Google Scholar] [CrossRef]
  66. Wang, K.; Shao, Y.; Zhou, W. Matheuristic for a two-echelon capacitated vehicle routing problem with environmental considerations in city logistics service. Transp. Res. Part D Transp. Environ. 2017, 57, 262–276. [Google Scholar] [CrossRef]
  67. Ahmadi-Javid, A.; Amiri, E.; Meskar, M. A Profit-Maximization Location-Routing-Pricing Problem: A Branch-and-Price Algorithm. Eur. J. Oper. Res. 2018, 271, 866–881. [Google Scholar] [CrossRef]
  68. Madankumar, S.; Rajendran, C. Mathematical models for green vehicle routing problems with pickup and delivery: A case of semiconductor supply chain. Comput. Oper. Res. 2018, 89, 183–192. [Google Scholar] [CrossRef]
  69. Abu Al Hla, Y.; Othman, M.; Saleh, Y. Optimising an eco-friendly vehicle routing problem model using regular and occasional drivers integrated with driver behaviour control. J. Clean. Prod. 2019, 234, 984–1001. [Google Scholar] [CrossRef]
  70. Madankumar, S.; Rajendran, C. A mixed integer linear programming model for the vehicle routing problem with sim-ultaneous delivery and pickup by heterogeneous vehicles, and constrained by time windows. Sādhanā 2019, 44, 39. [Google Scholar] [CrossRef]
  71. Zhou, H.; Qin, H.; Zhang, Z.; Li, J. Two-echelon vehicle routing problem with time windows and simultaneous pickup and delivery. Soft Comput. 2022, 26, 3345–3360. [Google Scholar] [CrossRef]
  72. Özceylan, E.; Demirel, N.; Çetinkaya, C.; Demirel, E. A closed-loop supply chain network design for automotive industry in Turkey. Comput. Ind. Eng. 2017, 113, 727–745. [Google Scholar] [CrossRef]
  73. Farrokh, M.; Azar, A.; Jandaghi, G.; Ahmadi, E. A novel robust fuzzy stochastic programming for closed loop supply chain network design under hybrid uncertainty. Fuzzy Sets Syst. 2018, 341, 69–91. [Google Scholar] [CrossRef]
  74. Jerbia, R.; Boujelben, M.K.; Sehli, M.A.; Jemai, Z. A stochastic closed-loop supply chain network design problem with multiple recovery options. Comput. Ind. Eng. 2018, 118, 23–32. [Google Scholar] [CrossRef]
  75. Almaraj, I.I.; Trafalis, T.B. An integrated multi-echelon robust closed- loop supply chain under imperfect quality production. Int. J. Prod. Econ. 2019, 218, 212–227. [Google Scholar] [CrossRef]
  76. Zhou, J.; Xia, W.; Wang, K.; Li, H.; Zhang, Q. Fuzzy Bi-Objective Closed-Loop Supply Chain Network Design Problem with Multiple Recovery Options. Sustainability 2020, 12, 6770. [Google Scholar] [CrossRef]
  77. Diabat, A.; Jebali, A. Multi-product and multi-period closed loop supply chain network design under take-back legislation. Int. J. Prod. Econ. 2020, 231, 107879. [Google Scholar] [CrossRef]
  78. Chouhan, V.K.; Khan, S.; Hajiaghaei-Keshteli, M. Metaheuristic approaches to design and address multi-echelon sug-arcane closed-loop supply chain network. Soft Comput. 2021, 25, 11377–11404. [Google Scholar] [CrossRef]
  79. Chiu, C.-Y.; Cheng, C.-Y.; Wu, T.-Y. Integrated Operational Model of Green Closed-Loop Supply Chain. Sustainability 2021, 13, 6041. [Google Scholar] [CrossRef]
  80. Easwaran, G.; Uster, H. A closed-loop supply chain network design problem with integrated forward and reverse channel decisions. IIE Trans. 2010, 42, 779–792. [Google Scholar] [CrossRef]
  81. Pedram, A.; Bin Yusoff, N.; Udoncy, O.E.; Mahat, A.B.; Pedram, P.; Babalola, A. Integrated forward and reverse supply chain: A tire case study. Waste Manag. 2017, 60, 460–470. [Google Scholar] [CrossRef]
  82. Guo, P.; Huang, G. Inexact fuzzy-stochastic mixed-integer programming approach for long-term planning of waste management—Part A: Methodology. J. Environ. Manag. 2009, 91, 461–470. [Google Scholar] [CrossRef]
  83. Li, Y.; Liu, J.; Huang, G. A hybrid fuzzy-stochastic programming method for water trading within an agricultural system. Agric. Syst. 2013, 123, 71–83. [Google Scholar] [CrossRef]
  84. Freeze, R.A.; Massmann, J.; Smith, L.; Sperling, T.; James, B. Hydrogeological Decision Analysis: 1. A Framework. Groundwater 1990, 28, 738–766. [Google Scholar] [CrossRef]
  85. Mendoza, G.A.; Bare, B.B.; Zhou, Z. A fuzzy multiple objective linear programming approach to forest planning under uncertainty. Agric. Syst. 1993, 41, 257–274. [Google Scholar] [CrossRef]
  86. Iskander, M.G. A suggested approach for possibility and necessity dominance indices in stochastic fuzzy linear programming. Appl. Math. Lett. 2005, 18, 395–399. [Google Scholar] [CrossRef]
  87. Kaufmann, A.; Aluja, J.G.; Kaufmann, J.A. Técnicas de Gestión de Empresa: Previsiones, Decisiones y Estrategias; FisicalBook: Madrid, Spain, 1992. [Google Scholar]
  88. Aghelie, A.; Sorooshian, S.; Azizan, N. Research gap in sustainopreneurship. Indian J. Sci. Technol. 2016, 9, 1–6. [Google Scholar] [CrossRef]
  89. Sorooshian, S.; Qi, L.C.; Fei, L.L. Characterization of ISO 14001 implementation. Environ. Qual. Manag. 2018, 27, 97–105. [Google Scholar] [CrossRef]
  90. Sorooshian, S.; Ting, K.C. Reasons for implementing ISO 14001 in Malaysia. Environ. Qual. Manag. 2018, 27, 125–133. [Google Scholar] [CrossRef]
  91. Pedram, A.; Pedram, P.; Yusoff, N.; Sorooshian, S. Development of closed–loop supply chain network in terms of cor-porate social responsibility. PLoS ONE 2017, 12, e0174951. [Google Scholar]
  92. Kabak, Ö.; Ruan, D. A comparison study of fuzzy MADM methods in nuclear safeguards evaluation. J. Glob. Optim. 2010, 51, 209–226. [Google Scholar] [CrossRef]
  93. Fallahpour, A.; Nayeri, S.; Sheikhalishahi, M.; Wong, K.; Tian, G.; Fathollahi-Fard, A. A hyper-hybrid fuzzy deci-sion-making framework for the sustainable-resilient supplier selection problem: A case study of Malaysian Palm oil industry. Environ. Sci. Pollut. Res. 2021, 1–21. [Google Scholar] [CrossRef]
Figure 1. A proposed closed-loop supply chain network.
Figure 1. A proposed closed-loop supply chain network.
Sustainability 15 02967 g001
Figure 2. Network configuration when α-cut level is 0.4 for the lower bound for problem number one.
Figure 2. Network configuration when α-cut level is 0.4 for the lower bound for problem number one.
Sustainability 15 02967 g002
Figure 3. Network configuration when α-cut level is 0.6 for the lower bound for problem number one.
Figure 3. Network configuration when α-cut level is 0.6 for the lower bound for problem number one.
Sustainability 15 02967 g003
Figure 4. Network configuration when α- cut level is 0.4 for the lower bound for problem number two.
Figure 4. Network configuration when α- cut level is 0.4 for the lower bound for problem number two.
Sustainability 15 02967 g004
Figure 5. Acceptable level of the objective function.
Figure 5. Acceptable level of the objective function.
Sustainability 15 02967 g005
Figure 6. Acceptable level of the objective function with a 20% increase in return products.
Figure 6. Acceptable level of the objective function with a 20% increase in return products.
Sustainability 15 02967 g006
Table 1. Literature comparison.
Table 1. Literature comparison.
Model CharacteristicsÖzceylan, Demirel [72]Farrokh, Azar [73]Jabbarzadeh, Haughton [26]Jerbia, Kchaou Boujelben [74]Almaraj and Trafalis [75]Zhou, Xia [76]Diabat and Jebali [77]Chouhan, Khan [78]Chiu, Cheng [79]This Paper
Product
Single X X XXX
MultipleXXX X X
Period
Single X X X
MultipleXXX X XXX
Modeling Approach
Deterministic XXX
Stochastic-robust optimization
Fuzzy programmingX X
Scenario-based robust optimization X
Mixed-integer linear programmingX XXXXXX X
Two-stage stochastic program X
Robust optimization X
Fuzzy-stochastic programming X X
Solution approach
Optimization software packageXX XX X XX
Lagrangian relaxation X
Metaheuristics algorithms X
Decomposition method X
Uncertain parameters
DemandXX XX X
CapacityXX X X
Cost X XXX X
Return XX X
Recovery rate X
Revenue X X
Error type X
Delivery time X
Disposal and repair X
Objective-function components
Transportation cost (min)XXXX XXXXX
Inventory cost (min) XX XXXX
Facility fixed-cost opening (min)XX X XX X
Penalty cost (min) XX
Disposal (min) X X X
Purchasing (min)XX X XX
Collection (min)XXX XX
Manufacturing (min) XXXXX XX
Disassembling (min)X
Recycling cost (min) X X X
Lost cost (min) X
Profit (max) X X
Remanufacturing cost (min) X X
Repair cost (min) X
Distribution cost (min) X
Procurement cost (min) X
Labor cost (min) X
Allocation cost (min) X
Processing cost (min) XX X
Table 2. The size of numerical examples.
Table 2. The size of numerical examples.
Problem No.No. of Potential SuppliersNo. of Potential ManufacturersNo. of Hybrid FacilitiesNo. of Existing CustomersNo. of Potential Recycling Centers
142683
26411174
35310156
45515205
Table 3. Trapezoidal fuzzy number for parameters.
Table 3. Trapezoidal fuzzy number for parameters.
ParameterValue
D ˜ r ~Unif (20,000, 30,000)~Unif (30,000, 40,000)~Unif (40,000, 50,000)~Unif (50,000, 60,000)
P ˜ r ~Unif (5000, 10,000)~Unif (10,000,15,000)~Unif (15,000,20,000)~Unif (20,000, 25,000)
C ˜ 1 v i ~Unif (3,000,000, 4,000,000)~Unif (4,000,000, 5,000,000)~Unif (5,000,000, 6,000,000)~Unif (6,000,000, 7,000,000)
C ˜ 2 k ~Unif (3,000,000, 4,000,000)~Unif (4,000,000, 5,000,000)~Unif (5,000,000, 6,000,000)~Unif (6,000,000, 7,000,000)
C ˜ 3 v j ~Unif (3,000,000, 4,000,000)~Unif (4,000,000, 5,000,000)~Unif (5,000,000, 6,000,000)~Unif (6,000,000, 7,000,000)
C ˜ S s i ~Unif (2,000,000, 3,000,000)~Unif (3,000,000, 4,000,000)~Unif (4,000,000, 5,000,000)~Unif (5,000,000, 6,000,000)
C ˜ M m ~Unif (2,000,000, 2,500,000)~Unif (2,500,000, 3,500,000)~Unif (3,500,000, 4,500,000)~Unif (4,500,000, 5,500,000)
C ˜ D h ~Unif (900,000, 1,000,000)~Unif (1,000,000, 1,500,000)~Unif (1,500,000, 2,000,000)~Unif (2,000,000, 2,500,000)
C ˜ C c ~Unif (900,000, 1,000,000)~Unif (1,000,000, 1,500,000)~Unif (1,500,000, 2,000,000)~Unif (2,000,000, 2,500,000)
F ˜ 1 h ~Unif (200,000, 250,000)~Unif (250,000, 300,000)~Unif (300,000, 350,000)~Unif (350,000, 400,000)
P ˜ S s i ~Unif (1, 3)~Unif (3, 6)~Unif (6, 9)~Unif (9, 12)
P ˜ M m ~Unif (2, 4)~Unif (4, 6)~Unif (6, 8)~Unif (8, 10)
P ˜ R c ~Unif (2, 4)~Unif (4, 6)~Unif (6, 8)~Unif (8, 10)
P ˜ C c i ~Unif (2, 4)~Unif (4, 6)~Unif (6, 8)~Unif (8, 10)
T ˜ C 1 m h v i ~Unif (10, 12)~Unif (12, 14)~Unif (14, 16)~Unif (14, 18)
T ˜ C 2 j j k ~Unif (10, 12)~Unif (12, 14)~Unif (14, 16)~Unif (14, 18)
T ˜ C 3 h c v j ~Unif (10, 12)~Unif (12, 14)~Unif (14, 16)~Unif (14, 18)
S P 20304050
Table 4. Optimum value for different sample problems for different levels.
Table 4. Optimum value for different sample problems for different levels.
Problem No.
α c u t   l e v e l 1234
α = 0 [190,528,300, 773,505,400][474,328,100, 1,719,479,000][590,739,600, 2,150,690,000][306,343,100, 1,191,767,000]
α = 0.2 [217,807,900, 726,333,400][527,852,300, 1,663,605,000][662,951,500, 2,150,690,000][342,264,400, 1,138,206,000]
α = 0.4 [246,618,800, 680,223,700][584,420,000, 1,494,726,000][727,251,800, 2,032,501,000][421,803,700, 1,076,690,000]
α = 0.6 [276,812,200, 6,302,216,00][631,762,400, 1,387,619,000][775,098,100, 1,748,506,000][464,338,500, 995,462,700]
α = 0.8 [308,375,900, 577,293,400][679,946,900, 128,397,3000][845,901,700, 1,599,929,000][500,146,300, 912,892,800]
α = 1 [341,418,500, 526,630,300][748,914,300, 1,183,985,000][986,029,000, 1,457,822,000][554,242,300, 833,767,800]
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

Pedram, A.; Sorooshian, S.; Mulubrhan, F.; Abbaspour, A. Incorporating Vehicle-Routing Problems into a Closed-Loop Supply Chain Network Using a Mixed-Integer Linear-Programming Model. Sustainability 2023, 15, 2967. https://doi.org/10.3390/su15042967

AMA Style

Pedram A, Sorooshian S, Mulubrhan F, Abbaspour A. Incorporating Vehicle-Routing Problems into a Closed-Loop Supply Chain Network Using a Mixed-Integer Linear-Programming Model. Sustainability. 2023; 15(4):2967. https://doi.org/10.3390/su15042967

Chicago/Turabian Style

Pedram, Ali, Shahryar Sorooshian, Freselam Mulubrhan, and Afshin Abbaspour. 2023. "Incorporating Vehicle-Routing Problems into a Closed-Loop Supply Chain Network Using a Mixed-Integer Linear-Programming Model" Sustainability 15, no. 4: 2967. https://doi.org/10.3390/su15042967

APA Style

Pedram, A., Sorooshian, S., Mulubrhan, F., & Abbaspour, A. (2023). Incorporating Vehicle-Routing Problems into a Closed-Loop Supply Chain Network Using a Mixed-Integer Linear-Programming Model. Sustainability, 15(4), 2967. https://doi.org/10.3390/su15042967

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