key: cord-1041059-2ht80hax
authors: Momenitabar, Mohsen; Dehdari Ebrahimi, Zhila; Arani, Mohammad; Mattson, Jeremy; Ghasemi, Peiman
title: Designing a sustainable closed-loop supply chain network considering lateral resupply and backup suppliers using fuzzy inference system
date: 2022-04-30
journal: Environ Dev Sustain
DOI: 10.1007/s10668-022-02332-4
sha: 30d6b553be26c6340c0acf72ada479fae884e520
doc_id: 1041059
cord_uid: 2ht80hax

Sustainability is key factor for transforming traditional supply chain networks into modern ones. This study, for the first time, considers the impacts of the backup suppliers and lateral transshipment/resupply simultaneously on designing a Sustainable Closed-Loop Supply Chain Network (SCLSCN) to decrease the shortage that may occur during the transmission of produced goods in the network. In this manner, the fuzzy multi-objective mixed-integer linear programming model is proposed to design an efficient SCLSCN resiliently. Moreover, the concept of circular economy has been studied in this paper to reduce environmental effects. This study aims to optimize total and environmental costs, including energy consumption and pollution emissions, while increasing job opportunities. A demand uncertainty component is considered to represent reality more closely. Due to the importance of demand, this parameter is estimated using the Fuzzy Inference System (FIS) as an input into the proposed mathematical model. Then, the fuzzy robust optimization approach is applied in a fuzzy set’s environment. The model is tackled by a Multi-Choice Goal Programming Approach with Utility Function (MCGP-UF) to be solved in a timely manner, and the equivalent auxiliary crisp model is employed to convert the multi-objective function to a single objective. The proposed model is tested on the case study of the tire industry in terms of costs, environmental impacts, and social effects. The result confirmed that considering the concept of lateral resupply and backup supplier could considerably decrease the total costs and reduce shortages on the designed SCLSCN. Finally, sensitivity analysis on some crucial parameters is conducted, and future research directions are discussed.

During the last two decades, the concept of Sustainability in Closed-Loop Supply Chain Network (CLSCN) has been widely expanded by multi-national corporations that have committed to reducing the undesirable effects of industrialization (Ahmed et al., 2020; Ebrahimi & Momenitabar, 2017) . Many businesses are still improving CLSCN by adopting sustainability factors into their network through considering real challenges, such as climate change, social legislation/regulation, environmental pollution, etc. Therefore, they must find ways to address these challenges that have unprecedented adverse effects on the entire network.

On the other hand, the concept of Circular Economy (CE) is increasingly considered in CLSCN, which empowers companies and organizations to collaborate, innovate in ways that assume resource scarcity and climate risk, and respond to consumer and community pressure to reduce waste. Indeed, CE aims to maintain components, products, and materials in CLSCN at the highest level of utility and its value always. This concept is a new approach for production and distribution companies to deal with sustainability challenges and resource consumption efficiency.

The uncertainty concept in designing a Sustainable Closed-Loop Supply Chain Network (SCLSCN) makes it a strategic necessity for companies to identify products' future demand and market share. In this way, the Fuzzy Inference System (FIS) is one of the efficient and widely used approaches in demand estimation utilized in this study (Sridharan, 2021) . The features of this approach not only make it possible to analyze and predict the state of complex markets but also lead to a comprehensive understanding of the fundamental market mechanism. Therefore, in this study, the criteria to forecast demand are identified using FIS to be imported into the mathematical model as a parameter discussed later in detail.

Our model was inspired by a real-world difficulty faced by manufacturing; however, because many businesses experience similar problems, the model was built to be general. The problem here is that customers order products based on their needs, and the manufacturers have to make the products customers want. On the other hand, manufacturers must obtain their raw material from primary or secondary suppliers, whichever opened/selected (output of proposed mathematical model). By obtaining raw materials, manufacturers build customers' products based on the Bill of Material (BOM). The produced products are sent to the distribution center to be delivered to the customers. During this time, customers may be unsatisfied with their products, so they might want to return them to the manufacturers. In this way, the returned items will be gathered at collection centers, and they will be checked in terms of suitability. After inspection, if the returned items/goods are detected to be usable again, they will send to the manufacturing centers to be repaired for reusing. If they are seen as not functional, they will send to the disposal center to be disposed. Based on that, in order to address this problem, in reality, we selected the tire industry in Iran as the case study of this research. The tire industry is a growing industry with recyclable goods that are expected to increase at a rate of 3.8 percent each year through 2025, according to predictions (Markets, 2025) . This growth in demand shows how the tire industry's actions can have an impact on a country's economy, environment, and society. Figure 1 shows a framework of this study for designing an efficient SCLSCN under uncertainty. This framework intends first to estimate demand using FIS and then optimize SCLSCN by minimizing the total costs of the network, the total environmental effects, and maximizing the social effects through applying the Fuzzy Robust Optimization (FRO) approach. Moreover, we consider both planning decisions, including the inventory level of products at manufacturing centers and distribution centers, the number of the finished product at manufacturing centers, and the number of product transferred between different centers of the network and strategic decisions, including the selection of primary and secondary suppliers and facility locations into the model. Indeed, this study considered the secondary suppliers as a backup supplier to support primary suppliers in the first level and concept of lateral transshipment/resupply through substituting products between manufacturing centers as well as distribution centers to avoid shortage during the transmission of produced goods in the network.

We develop a novel fuzzy multi-objective mixed-integer linear programming (FMO-MILP) model to design an efficient SCLSCN by minimizing the total costs of the network, the total environmental effects, and the social effects. This study, for the first time, considered the backup suppliers and lateral transshipment/resupply simultaneously to decrease the shortage that may occur during the transmission of produced goods in the network. Because of the essence of demand parameters, we estimate this parameter using FIS that will be used as input into the proposed mathematical model. Then, we utilized the FRO approach to present the mathematical model in a fuzzy sets environment and tackled it by Multi-Choice Goal Programming Approach with Utility Function (MCGP-UF) approach in a timely manner. Furthermore, before solving the proposed model, we employed the Equivalent Auxiliary Crisp (EAC) model to convert the multiobjective function to a single objective. The model's computational requirements are analyzed as well. Regarding this, we are looking for answers to the following questions as follows: Fig. 1 The framework of this research • Which centers of the SCLSCN should be opened to transfer the finished products to the customer and return it to the collection centers? • How many products should be transformed to avoid shortages that may occur at manufacturing centers and distribution centers? • What would be the total costs, total environmental effects, and social effects of various centers when they are opened?

The remaining parts of the paper are organized as follows. In Sect. 2, the literature review on SCLSCN are discussed. The main contribution of this work has been explained in Sect. 3. Also, the mathematical modeling and solution methodology are proposed and discussed in Sects. 4 and 5. The tire case study, results, and the sensitivity of the model are investigated in Sect. 6. Finally, Sect. 7 provides concluding remarks and future research directions of this study.

In this section, we briefly present some works done in the past by authors and some literature reviews on SCLSCN design. Some comprehensive literature reviews have been conducted in recent years (Battini et al., 2017; Govindan et al., 2015; ıldızbaşı, 2021; Khalili Nasr et al., 2020; Mehrjerdi & Shafiee, 2020; Rajak, 2021; Souza, 2013; Zhalechian et al., 2016) . In the following paragraphs, we will discuss different aspects of SCLSCN comprehensively.

The SCLSCN consists of two parts, including forward and reverse flow (Cilacı Tombuş et al., 2017) . A lack of concern for uncertainty in the real world can prevent companies from meeting the needs of their customers. Uncertainty in today's world is an essential part of every SCLSCN. Many researchers have considered the concept of uncertainty in both demand and supply when designing a SCLSCN (Ahmed et al., 2020; Arani et al., 2021; Zhalechian et al., 2016) . Zhalechian et al. (2016) formulated a stochastic-possibilistic programming model to handle the uncertainty of SCLSCN for large-sized problems by using a hybrid meta-heuristic algorithm and lower bounds. Mehrjerdi et al. (2020) included the interaction of sustainability and resilience simultaneously to design a SCLSCN. They used the fuzzy TOPSIS method to rank supply chain strategies, such as holding additional inventory, backup capacity, multiple sourcing, and so on (ldızbaşı, 2021). Khalili Nasr et al. (2020) utilized the fuzzy best-worst method to select the best suppliers. By proposing a multi-objective mixed-integer linear programming model, they could design a multi-product, multi-period, SCLSCN. Uncertain parameters such as demand and return products were considered by Amin and Zhang (2013) . They employed a scenario-based analysis to face uncertainty in their proposed model by using stochastic programming formulation and applying the weighted sums method as well as the Epsilon constraint method to solve their proposed model. Burton et al. (2005) analyzed the cost effects of two lateral (intra-echelon) transshipment approaches in a two-echelon supply chain network, with a single supply source at the higher echelon and multiple retail locations at the lower, using a series of simulation experiments under different operating conditions. Ghomi-Avili et al.

(2019) proposed a multi-objective model for designing the CLSCN with a price-dependent demand while considering random disruptions and shortage. They have considered various resilience strategies, including lateral transshipment among production centers and multisource allocation to model the resilient supply chain network design problem. Ahmed et al. (2020) designed a new green supply chain design approach for maintaining the financial virtue accompanying the environmental factors, proposing a multi-objective mathematical model to minimize the total costs and CO 2 emissions. They solved the model by applying two optimization models, ε-constraint method, and genetic algorithm (GA). Goodarzian et al. (2021) used a production-distribution-inventory-allocation-location problem to design a sustainable medical supply chain network by considering three hybrid meta-heuristic algorithms, including ant colony optimization, fish swarm algorithm, and firefly algorithm to solve their model. Soleimani et al. (2017) studied a SCLSCN with multiple levels, multiple products, and multiple periods by applying the GA to solve their model (Safaei et al., 2022) . Sahebjamnia et al. (2018) considered the SCLSCN to minimize total costs, as well as the environmental impacts of the facilities, and to maximize the social effects by proposing four new hybrid meta-heuristic algorithms to decrease the computational time and improve the quality of solutions. Rezaei and Kheirkhah (2018) proposed a mathematical model for designing the SCLSCN based on economic, environmental, and social requirements by applying cross-docking operations. Indeed, they utilized a cuckoo optimization algorithm for the first time to tackle their model.

According to Hillier and Lieberman (2012) , the lateral transshipment concept is a variant of the minimum cost flow problem in which products move through intermediary transfer locations on their way from one source to another. Reduced total network costs, increased service levels (Tracht et al., 2011) , and reduced interruption risk are all advantages of transshipment (Jabbarzadeh et al., 2018) . Only a few scholars have incorporated transshipment into their models, despite the fact that it delivers the aforesaid benefits and simplifies CLSC design. Arani et al., (2021) presented a mixed-integer mathematical programming model to design a sustainable supply chain network design under uncertainties through considering the concept of lateral transshipment. They showed that considering this concept could allow hospitals to satisfy their demand from more than one source. Jabbarzadeh et al. (2018) considered an SCLSCN subject to random disruption risks. They found that including lateral transshipment was useful to overcome disruption risks. Another study by Samuel et al. (2021) integrated a model to design a multi-component multi-product SCLSCN with economic scale to locate the processing centers, transshipment of inspected products between centers, and reliability of returned products. Lastly, Tracht et al. (2011) applied a simulation model to design a multi-echelon spare-part closedloop supply chain, in which spare parts are repaired and led back to stock after removal from a broken machine. Also, they found that utilizing the concept of lateral transshipments leads to reducing the total cost of the model. It was also found that higher loan costs resulted in longer stock-outs. For more information, readers are referred to Ahmed et al. (2020) , Ebrahimi & Momenitabar (2017) , Safaei et al. (2022) , Momenitabar et al. (2013) .

Some studies considered the concept of CE in SCLSCN Nasir et al., 2017) . Abdul Nasir et al. (2017) applied a case study in the construction industry to show how the environmental effects by including carbon emissions can be achieved by CE principles. They could benefit from integrating the concept of CE into a sustainable supply chain network in terms of environmental gains. Genovesa et al. (2017) considered two studies to compare the performance of traditional and circular production systems. They concluded that the integration of CE with SCLSCN provided clear advantages in terms of environmental effects. Rashedul Kabir et al. (2021) proposed an improved network design of an open-loop reverse supply chain of a water treatment plant using circular economy that had benefits such as reducing environmental burdens, enhancing the protection of raw material supply, increasing competition, promoting productivity, boosting sustainable economic growth, and creating employment. Computational experiments showed that when the number of the returned product increases, the total profit increases, which leads to achieving the circular economy with the most efficiency.

Some studies designed a decision support system to predict the demand for medical equipment in the event of a COVID-19 outbreak pandemic (Govindan et al., 2019) . A FIS is used to predict the amount of demand. Criteria for predicting demand include patient age, fatigue, number of coughing hours, fever, and fatigue. The results indicated the proper performance of the proposed system (Govindan et al., 2019) . Nabavi et al. (2019) predicted white rice demand by considering some criteria, including transportation, machinery, electricity, and human labor. The neural network approach has been used to calculate the weight of the criteria. The results showed an accurate prediction of the proposed model. Souza et al. (2020) predicted the demand of the beauty industry by four FIS-based methods. The result indicated that the proposed predicting models for short-term and long-term performance had performed well.

To the authors' knowledge, a few studies considered the concept of lateral transshipment/ resupply by substituting the final products at manufacturing and distribution centers to prevent shortage as well as concepts of the backup supplier in their studies. Indeed, lateral transshipment/resupply is a type of inventory system that permits manufacturing center to satisfy its demand by other manufacturing inventories in the absence of other required products. Also, considering backup suppliers help the manufacturing centers to meet their needs. More interestingly, demand estimation is an essential task for reaching more accurate results which plays as an input to the model by solving the mathematical model of this study. Regarding this, the main contribution of this study that makes a sharp distinction between this study and related studies are as follows:

• Proposing a new formulation of designing a SCLSCN by applying Fuzzy Robust Optimization model in the fuzzy set environment, • Utilizing the Multi-Choice Goal Programming with Utility Function (MCGP-UF) to solve the proposed mathematical model in a deterministic environment, • Considering the backup supplier as supporter role for a primary supplier in the first level of the network as well as the concept of lateral transshipment/resupply in manufacturing and distribution centers to prevent shortages, and • Utilizing the customized Fuzzy Inference Systems method to estimate the required demandas an input to the designed network.

This section firstly plans to propose the problem statement and then presents the mathematical model of this study. Before proposing the mathematical model, we will explain how we estimate the demand parameter using FIS as an input to the mathematical model .

This paper addresses the design of multi-echelon, multi-products, and multi-periods SCLSCN, which comprises various layers, including primary suppliers and secondary suppliers at the first level, manufacturing centers at the second level, distribution centers at the third level, customers at the fourth level, and collection and disposal centers at the fifth and sixth levels, respectively. The SCLSCN studied in this paper is constituted by forward and reverse flows. Primary suppliers, secondary suppliers, manufacturing centers, distribution centers, and customers are placed in forwarding flow, while collectors and disposal centers are placed in reverse flows of the proposed network. In the first level, primary suppliers and secondary suppliers are responsible for preparing the raw material for manufacturing centers at the second level by sending them to the manufacturing centers to produce the products according to BOM at manufacturing centers. Then, the made products are shipped to the customers through distribution centers. In the reverse flow, the returned products are gathered by collection centers. The purpose of collecting the returned products by customers is to fulfill the needs of customers through recovering or repairing the returned products. The returned products are sent to the collection centers to check if they can be used again or not. After careful inspection, the returned products are either returned to the manufacturing centers or sent to the disposal centers as scrapped products.

The demand for products happens in the fourth layer, where the customer is located. As we can see from Fig. 2 , the demand can be fulfilled by sending generated products from the distribution centers to the customers. On the other hand, due to the nature of the case study discussed in Sect. 6, the demand for this study is considered as an uncertain parameter as estimated using FIS that will be addressed in the following sub-section.

The FIS procedure applied in this research is shown in Fig. 3 . First, the criteria for determining the amount of demand are selected by experts and converted to a fuzzy number. Then, fuzzified criteria are entered into the fuzzy inference engine. In this process, linguistic membership functions are first defined, and then the desired operators are customized. In the next step, the rules set by the experts are evaluated and aggregated by the relevant operators. In the last step, the required amount of demand is estimated.

Mamdani Inference System was first introduced by Mamdani et al. (1975) . Due to its visual nature, this system can be widely used in decision support systems. The steps of Mamdani FIS are described in the following sections.

In this study, an expert is a person who has at least five years of experience in relevant work and at least a bachelor's degree. According to the statistical population, 30 people have the conditions as an expert to be selected. Based on experts' opinions, five criteria have been chosen by using the Fuzzy Delphi Approach to predict demand. Selected criteria include the product price, product quality, environmental impact, advertising, and product availability. In Fig. 4 , we have shown these criteria selected by the experts. The terms Moderate Importance (MI), Weak Importance (WI), and Strong Importance (SI) have also been used to weight the selected criteria. Table 1 shows how these three criteria are weighted.

In this step, the membership function, inputs, and outputs (demand) are determined. The numbers considered in this function are fuzzy triangular. Also, the linguistic terms considered in this research include High (H), Medium (M), and Low (L). Figure 5 shows the fuzzy structure that has been considered in this study.

Mamdani Inference System has been used as fuzzy operators. Therefore, the operator's min and multiplication are used to process each rule. Equations (1) and (2) show how these operators are used in this study. 

Rules evaluation is the central part of the FIS, and it is determined as "if-then" based on expert opinion. A fuzzy rule can be "if x 1 is a 1 AND x 2 is b 1 THEN y is c 1 " so x 1 , x 2 , y are considered variables, and a , b, c are the linguistic variables. As an example, Table 2 shows the rules produced by Decision Maker 1. Rule 1, for example, states that demand is high if price, quality, environmental impact, advertising, and availability are low, medium, low, medium, and high, respectively.

This sub-section details the mathematical model of this study, as illustrated in Sect. 4.1. Also, before proposing the formulation of this study, we need to define some assumptions to build the model in a proper way. The assumptions include but are not limited to:

• The concept of "Lateral Resupply" is considered to avoid the shortage that may occur at manufacturing centers and distribution centers. • The length of each period is one day. • The candidate locations of primary and secondary suppliers are fixed.

• Manufacturing centers are fed by both primary and secondary suppliers. Indeed, the secondary suppliers play a role as a supporter for primary suppliers to avoid shortages in different periods. • Capacity is limited at each level of the proposed SCLSCN.

• The demand and supply are considered uncertain. Also, the demand parameter is estimated by the FIS procedure, as is discussed in Sect. 4.2 in detail. Figure 6 displays the framework of the proposed model, considering main inputs or parameters, objective functions that include two minimizing objective functions and one maximizing objective function, constraints that build the feasible solution of the proposed model, and main outputs that state the number of decisions variables in a deterministic space as they will be reported in Sect. 6.

The sets, parameters, and variables of this study are presented in detail in Tables 3, 4 , 5, 6, and 7. Fig. 6 The framework of the proposed model After defining sets, parameters, and variables, the mathematical model of this study is presented as follows:

(3) 

The set of secondary suppliers

The set of product type ( e = 1, 2 … , E) a

The set of raw material ( a = 1, 2 … , A) t

The set of the period ( t = 1, 2 … , T)

The first objective function, which is shown by Eq. (8) consists of five sub-functions that aim to minimize the total costs of the network. The first term is Eq. (3) that minimizes the total fixed costs of SCLSCN, while Eq. (4) aims at minimizing the total transportation (8) 

The demand of product type e in customer zone c in period t costs of the network. Equation (5) minimizes the total production cost at manufacturing centers as well as purchasing prices of products at all the centers except primary suppliers, secondary suppliers, and customer's centers. Equation (6) is designed to minimize the number of inventories at manufacturing centers and distribution centers, while Eq. (7) minimizes the disposal costs of disposal centers for scrapped products returned from the collection centers after inspection.

(9) Holding cost of product type e in collection center o in period t

Disposal cost of product type e in disposal center r in period t a The quantity of raw material a to produce one finished product 1 ec

The percentage of scrapped products type e returned from customer zone c 2 eo

The percentage of products type e scrapped from collection center o The second objective function, shown by Eq. (11), plans to minimize total emissions and resource consumption of the network, which is embodied by two parts. The first part is related to the water consumption, electricity consumption, and pollution created by these utilities, which is shown by Eq. (9). While, the second part, Eq. (10), tries to minimize the emissions created by transporting finished products and raw materials between different centers of the network.

(10) The amount of transferring raw material a from secondary suppliers b to manufacturing center m in period t The number of product type e that is transshipped from manufacturing center m to manufacturing center m′ in period t

The number of product type e that is transshipped from distribution center d to distribution center d′ in period t Q emt

The production quantity of product e at manufacturing center m in period t

The third objective function proposed by Eq. (14) aims to maximize the social effects comprised of the total number of workers employed when some centers are opened. This objective function is structured into two parts, including fixed-and variable-job opportunities which are shown by Eqs. (12) and (13).

After proposing the three objective functions above, we need to build a feasible solution by providing some constraints.

Assignment constraints:

The first group of constraints is assignment constraints which are stated by Eq. (15)-(18). Indeed, the aims of constructing these equations are to set up some centers in the proposed network Capacity constraints:

Equations (19)-(24) ensure that the number of products or raw materials that have been transferred from one center to another center should not exceed the capacity of the original center. This statement is written mathematically by Eqs. (19)-(24) for different layers of the network. Additionally, two Eqs. (25) and (26) are written for the amount of supply from first layers that primary and secondary suppliers are located to second layers (manufacturing centers) that should be at least a specific number to be used by manufacturing centers to produce goods.

Supply constraints:

Equations (27)-(32) are related to balanced constraints that show the flows of products between different layers of the network. For instance, Eq. (27) indicates that the number of raw materials sent from primary and secondary suppliers to manufacturing centers should be equal to the number of products produced in manufacturing centers. These statements are acceptable for other constraints such as Eq. (32) as well.

Balanced constraints:

The inventory equations are written by Eqs. (33)-(35), which display the amount of inventory in each period in manufacturing centers, distribution centers, and collection centers as well. For each center, including manufacturing centers, distribution centers, and collection centers, we have written the inventory equations. For example, Eq. (33) indicates that inventory level for manufacturing center should come up with inventory level of each product in the previous period, number of produced products in the manufacturing center, number of returned products from collection center, number of products sent out from the

manufacturing center and number of substituted products between manufacturing centers (whichever have inventory level to meet other centers to avoid shortage). Inventory constraints:

Finally, the nature of each defined variable is shown by Eqs. (36) and (37).

This section presents the solution methodology to solve the proposed mathematical model in the previous section. Also, the flowchart of the proposed solution procedure is depicted in Fig.7 . What is clear from Fig. 7 is that firstly, the demand was estimated using FIS, and then it was used as a parameter to be imported to a mathematical model. After that, the FMO-MILP model is proposed in Sect. 3.3 in a fuzzy environment. Then, it is converted to a deterministic environment by using an Equivalent Auxiliary Crisp (EAC) model (Jiménez et al., 2007) . Finally, the Deterministic Multi-objective Mixed-Integer Linear Programming (DMO-MILP) model is solved by utilizing the Multi-Choice Goal Programming Approach with Utility Function (MCGP-UF) to obtain the results of the proposed model.

As we discussed in the previous section, some parameters in the proposed model were uncertain. Uncertainty can be classified into different types, including randomness, epistemic, and deep uncertainty. Figure 8 shows the classification of uncertain data based on the essence of data availability. Because of the nature of our research problem, the type of uncertainty that we use is epistemic (Bairamzadeh et al., 2018) .

The Equivalent Auxiliary Crisp (EAC) model was first introduced by Jiménez et al. (2007) . The EAC converts optimization problems that have triangular/trapezoidal fuzzy numbers (33) 1}, ∀m, d, o, r, p, b (37) ∀e, a, p, b, m, d, c, o, r, t to deterministic problems. The first term of the first objective function, which is shown by Eqs. (3), (12), and (13), are converted to Eqs. (38), (39), and (40) as follows: Fig. 7 The flowchart of the proposed solution Also, constrains (19)- (24) and (30) are fuzzy constrains that need to be converted to deterministic constraints. Then, we have: 

In this sub-section, an efficient solution methodology called Multi-Choice Goal Programming Approach with Utility Function (MCGP-UF) is proposed to solve the deterministic multi-objective MILP model in Sect. 5.1. Since the proposed FMO-MILP model has multiple decision variables, it needs to be solved within a feasible solution space in a faster time. While considering a case study, the proposed model is coded on a computer with an intel CPU @2.9 GHz processor and 64 GB of RAM using CPLEX solver.

The Goal Programming (GP) method is an efficient tool that is still widely used to solve multi-objective and multi-criteria decision-making problems. In the goal programming method, an aspiration level is defined for each objective function and is calculated based on an expert's opinion. This paper applies the MCGP-UF approach proposed by Chang (2011) . Chang (2011) presented a methodology to solve the model by substituting the aspiration level with the scalar value for multiple objective problems, including classic GP and MCGP in terms of the level achieved for the Utility Function (UF). Indeed, the reasons for utilizing this technique are to aid decision-makers in developing the best strategy for their goals and to achieve the highest level of usefulness Nayeri et al., 2018) . The MCGP with UF applied in this study is defined as follows: y k : continuous variable, Uk, max k,min : rank of kth aspiration level. d + k , d − k : negative and positive deviations, − k normalized deviation of y k from U k,min ,w k : weight of − k , k : utility value. According to Change (2011), the normalized MCGP with UF can be shown as follows: − k does not need to be normalized because it is between 0 and 1. On the other hand, f + k = min f k (X) and f − k = max f k (X) . To calculate the upper and lower level of each objecti function (aspiration level), we need to calculate Uk, max k,min .

In this section, we considered the tire industry as a case study of this study to validate the proposed mathematical model in Sect. 4. Then, we report the result of the proposed model by running the model using the CPLEX solver.

The tire industry in Iran is interdisciplinary. Figure 9 shows the share of different tire factories production in Iran. As we can see, the Barez factory is located at the top with 36.40%, and the Kian factory and Pars factory are among the lowest level with 4.90%. In this study, the Barez Tire factory was used to validate the proposed model (Shabani, 2018) . The Barez tire factory was established in Kerman/Iran, and the construction of the factory began with a nominal capacity of 25,500 tons (Shabani, 2018) . The company's operating capacity is currently 27,900 tons of all kinds of tires and rims (flaps) per year. The factory also provides direct and full-time employment for 1200 people and indirect employment for 2000 people. The main goal of this factory is to minimize supply chain costs, minimize pollution emissions and resource consumption, and maximize social impacts. The considered period includes three periods. In this research, two products, Truck Tire (product-type 1) and Bus Tire (product-type 2), have been considered for planning. Table 8 shows the amount of demand estimated by the FIS in each period. For example, the demand for product 1 and first customer in period 1 is equal to (1000, 1200, 1400) . Also, Tables 9, 10, and 11 have shown the number of input parameters for the proposed model. As discussed in Sect. 4.3, some parameters are certain numbers, and others are fuzzy numbers. So, these parameters are defined and numbered through uniform distribution and fuzzy triangular numbers, as shown in Tables 9, 10, and 11.

As discussed in Sect. 3.4, some parameters are certain, and others are fuzzy numbers. So, these defined parameters are numbered through uniform distribution and fuzzy triangular numbers, as shown in Tables 10 and 11.

This section provides the results of the proposed mathematical model in the tire industry by applying the MCGP-UF approach. Firstly, the weight of each objective function is considered 0.3, 0.4, 0.3, and the parameters of MCGP-UF are given in Table 12 . Secondly, the deterministic model is solved under six scenarios-two optimum values are sought for each objective function-to obtain the required results discussed in Table 13 .

As we stated previously, the strategic decisions of the proposed model are pertinent to the selection of primary and secondary suppliers and facilities' locations. The result of the proposed model is shown in Table 14 . Table 14 displays which suppliers and facility centers are opened in the network for all periods. As we can see, secondary supplier number 1 is selected as the supportive supplier to help primary suppliers to avoid any interruptions between levels 1 and 2 of the networks. Two manufacturing centers, as well as two distribution centers, are selected to produce the products and deliver them to the customers. Additionally, one collection center and one disposal center are opened to collect the returned products from customers and for scrapped items that are not able to be reused. 

This sub-section shows the value of planning variables such as the number of finished products transferred from distribution centers to customers and the substitution between distribution centers. For instance, Table 15 shows that in period one, 648 units of product type 1 are transferred from distribution center 2 to customer 1. Moreover, Table 16 shows the number of finished products substituted between distribution centers. For instance, in the third period, 1139 units of products type 1 are replaced between distribution centers.

In this sub-section, sensitivity analyses on demand and lateral resupply are presented. By doing the sensitivity analysis, we aim to understand the variations on demand as well as supply how much can impact the first, second, and third objective functions.

The model was solved considering different levels of demand (Base case, − 10%, − 20%, + 10%, and + 20%). Results are shown in Figs. 10, 11, and 12. Figure 10 shows that as demand increases from − 20 to + 20%, the value of the first objective function, which is total costs, increases from approximately $1,200,000 to $1,700,000. More interestingly, a positive relationship is also found between demand and the value of the second objective function where demand is increased (see Fig. 11 ). Lastly, Fig. 12 shows that by increasing demand, the value of the third objective function, which is related to the maximization of social effects, also increases. In other words, when the demand increased, the number of direct and indirect employees also increased.

Another critical parameter for SCLSCN is to consider the supply amount. Sensitivity analysis was also conducted for the impacts of supply variations on the three objective functions (see Figs. 13, 14, and 15 ). The sensitivity analysis of supply variations is periodic for the second objective function (both sides of Fig. 14) and the right sides of Fig. 13 and 15. As we can see, these figures are repeated intermittently by changing the supply amount in each period. On the other hand, the decreasing orientation occurs for the left side of Fig. 13 . The same happens for the left side of Fig. 15 , except for a 20% increase in demand causes the third objective function to increase. 

This sub-section is dedicated to the examination of how lateral resupply could have effects on the total costs of the SCLSCN. The analysis first considers the network without lateral resupply and then analyzes the same network while allowing for lateral resupply. Results are shown in Fig. 16 .

Results show that lateral resupply could considerably decrease the total cost of the SCLSCN. In other words, this strategy could remarkably help managers to lower the total costs of the network. The lateral resupply is considered in Eqs. (33) and (34) the shortage that may occur in the network by substituting the products between manufacturing centers at the second level and distribution centers at the third level.

This study addressed the design of multi-echelon, multi-products, and multi-periods SCLSCN, which comprises primary suppliers, secondary suppliers, manufacturing centers, distribution centers, customers, collection centers, and disposal centers. The proposed mathematical model aimed at minimizing the total costs of the network, minimizing the environmental effects, and maximizing social outcomes. This study, for the first time, considered the backup suppliers and lateral transshipment/resupply simultaneously to decrease the shortage that may occur during the transmission of produced goods in the network by proposing an FMO-MILP model. Also, the concept of circular economy is included in this study to reduce environmental effects. Given the importance of the demand, this parameter was estimated by the FIS procedure. After screening, the criteria including price, quality, environmental impact, advertising, and availability were selected as the final criteria. After fuzzification and determining the membership functions, the criteria are weighted, and then the chosen criteria and rules enter the FIS procedure. Finally, the output of the proposed FIS was the amount of product demand that was imported to the mathematical model as a parameter . This study utilizes the MCGP-UF approach to solve the FMO-MILP in Sect. 4. Before using this approach, we employed the EAC model to convert the multi-objective function to a single objective one. Location and allocation of centers, along with the study of the flow of raw materials and products, have been among the decisions made in this study. In addition, in this research, the level of inventory of the collection, production, and distribution centers and the amount of production of products in each period has been determined. The solution confirmed that the proposed model could reach the target level as the decision-makers defined. Also, we examined each objective function's response using sensitivity analysis based on demand variations. By conducting this study, we found that increasing demand led to an increase in socialeconomic impacts such as creating direct jobs and indirect jobs at the tire manufacturing companies while leading to a decrease in pollution emission and energy consumption (see Figs. 10, 11, and 12) . Also, by conducting sensitivity analysis on supply variation, we understood that the supply variations are periodic from one period to another period. (see Figs. 13, 14, and 15) . By considering the concept of lateral resupply and secondary suppliers in the designed SCLSCN, we have shown in Fig. 16 that the proposed methodology could decrease the total costs as considered in the first objective function.

The limitation of this research should be addressed as well. Disruption can be considered in various centers of designed SCLSCN, including primary suppliers, secondary suppliers, manufacturing centers, disposal centers, and distribution centers . Indeed, disruption can have negative impacts on the whole network as it disturbs the supply chain of transporting goods. Considering other decision variables such as vehicle routing and scheduling can be other ideas for future research as well. Also, another study can be prioritizing backup suppliers on sustainable criteria. There are many approaches for handling uncertainty for various types of data availability (Fig. 8) , such as two-stage stochastic programming, fuzzy robust programming, and robust convex programming. Lastly, the proposed model can be solved in large size by considering several metaheuristic algorithms and using criteria reported in Zhang and Gen (2014) to compare the algorithms.

An environmentally sustainable closed-loop supply chain network design under uncertainty: Application of optimization

A multi-objective facility location model for closed-loop supply chain network under uncertain demand and return

Optimizing the total production and maintenance cost of an integrated multi-product process and maintenance planning (IPPMP) model

A lateral resupply blood supply chain network design under uncertainties

Modelling different types of uncertainty in biofuel supply network design and planning: A robust optimization approach

Closed loop supply chain (CLSC): Economics, modelling, management and control

Cost-parametric analysis of lateral transshipment policies in two-echelon supply chains

Multi-choice goal programming with utility functions

Designing distribution systems with reverse flows

Design of mathematical modeling in a green supply chain network by collection centers in the environment

Sustainable supply chain management and the transition towards a circular economy: Evidence and some applications

A multi-objective model for the closed-loop supply chain network design with a price-dependent demand, shortage and disruption

An integrated sustainable medical supply chain network during COVID-19

A decision support system for demand management in healthcare supply chains considering the epidemic outbreaks: A case study of coronavirus disease 2019 (COVID-19)

Closed-loop supply chain network design under disruption risks: A robust approach with real world application

Linear programming with fuzzy parameters: An interactive method resolution

A novel fuzzy multi-objective circular supplier selection and order allocation model for sustainable closed-loop supply chains

An experiment in linguistic synthesis with a fuzzy logic controller

Automotive tire market to 2025 -global analysis and forecasts by tire type; vehicle type; and distribution channel

A resilient and sustainable closed-loop supply chain using multiplesourcing and information sharing strategies

A proposed lean distribution system for solar power plants using mathematical modeling and simulation technique

Fuzzy Mathematical modeling of distribution network through location allocation model in a three-level supply chain design

Comprehensive model of energy, environmental impacts and economic in rice milling factories by coupling adaptive neuro-fuzzy inference system and life cycle assessment

Comparing linear and circular supply chains: A case study from the construction industry

Goal programming-based post-disaster decision making for allocation and scheduling the rescue units in natural disaster with time-window

A comprehensive approach in designing a sustainable closed-loop supply chain network using cross-docking operations

Multi-objective mixed-integer linear optimization model for sustainable closed-loop supply chain network: a case study on remanufacturing steering column

Application of mamdani fuzzy inference system in predicting the termal performnce of solar distillation still

Designing a new multi-echelon multiperiod closed-loop supply chain network by forecasting demand using time series model: A genetic algorithm

Sustainable tire closed-loop supply chain network design: Hybrid metaheuristic algorithms for large-scale networks

Multicomponent multiproduct closedloop supply chain design with transshipment and economies of scale considerations

Tire industry analysis in Iran

Fuzzy multi-objective sustainable and green closed-loop supply chain network design

Closed-loop supply chains: A critical review, and future research

Demand forecasting in the beauty industry using fuzzy inference systems

Effects of lateral transshipments in multi-echelon closedloop supply chains. Glocalized solutions for sustainability in manufacturing

Assessing the social sustainable supply chain indicators using an integrated fuzzy multi-criteria decision-making methods: A case study of Turkey

Exploring the knowledge base and trajectories of knowledge dissemination in closed loop supply chain

Sustainable design of a closed-loop location-routing-inventory supply chain network under mixed uncertainty

Hybrid sampling strategy-based multiobjective evolutionary algorithm for process planning and scheduling problem

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations