Interactive Fuzzy Multi Criteria Decision Making Approach for Supplier Selection and Order Allocation in a Resilient Supply Chain

Abstract

Modern supply chains are vulnerable to high impact, low probability disruption risks. A supply chain usually operates in such a network of entities where the resilience of one supplier is critical to overall supply chain resilience. Therefore, resilient planning is a key strategic requirement in supplier selection decisions for a competitive supply chain. The aim of this research is to develop quantitative resilient criteria for supplier selection and order allocation in a fuzzy environment. To serve the purpose, a possibilistic fuzzy multi-objective approach was proposed and an interactive fuzzy optimization solution methodology was developed. Using the proposed approach, organizations can tradeoff between cost and resilience in supply networks. The approach is illustrated using a supply chain case from a garments manufacturing company.

1. Introduction

Outsourcing is a competitive strategy in the global supply chain. Evaluating and selecting the best set of suppliers is a challenging decision in outsourcing and it plays a significant role in supply chain performance [1,2]. Traditionally the supplier selection and order allocation decision is made based on cost and quality criteria. However, modern supply chains are more prone to unexpected High Impact Low Probability (HILP) and Low Impact High Probability (LIHP) disruption events [3]. HILP disruption events are commonly known as random disruptions risks such as man-made and natural disasters, whereas LIHP disruptions are targeted disruptions such as day-to-day operational risks. Tang and Tomlin [4] proposed six disruption sources in the supply chain and among them, supplier performance is most frequent. The role of the supplier selection decision in these supply chain risk has only been partially explored in the literature [5]. Multi-sourcing strategies and are now common to many supply chains in order to minimize the supplier’s disruption risks [6]. For example, during a fire in a plant of Philips Electronics in 2010 disrupt two of its major customers: Ericsson and Nokia. Ericsson lost about a month of production and suffered $200 million while Nokia recovered due to its multi-sourcing ability [7]. Toyota Motor Corp. lost billions of dollars in 2010 during product recall due to its part sourcing from one supplier for many car models. These examples show that multi-sourcing strategies work well. On the contrary, multi-sourcing strategies failed during some HILP disruption events. For example, Japan earthquake disrupted many semiconductor supply chains. Chinese Firm ZTE Corp. faced shortages of batteries and LCD screens due to all of its suppliers in the affected region. Ford Motor Co. and General Motor Co. faced shortages of auto parts and stop production due to the shutdown of two Hitachi Ltd.’s plants. These historical events suggest that the supplier selection criteria should be extended to new resilience capabilities [8]. Thus it crucial to provide a reliable level of resilience to the supply side to protect such shortages especially during HILP events [9].

Several studies have been conducted to consider resilience in the supply chain [8,10,11]. The concept of resilience in specific to the supplier selection problem has also been discussed by several authors [12,13,14,15,16,17,18,19,20,21,22,23,24,25]. Most of these studies focused on multiple sourcing and operational performance of suppliers. However, to the best of authors knowledge, this is the first study which focuses on supply network by considers supply density, resilience score of supplier’s locations, and transit time in addition to other operational criteria. This paper aims to develop a supplier selection and order allocation model to build a resilient supply chain in response to HILP disruptions. To do so, a possibilistic multi-objective fuzzy optimization-based model with a new resilience objective is proposed which consists of supply density, resilience score, and transit time. Furthermore, the proposed model is solved using Tiwari, et al. [26] weighted additive approach and Werners [27] fuzzy and operator methods. Fuzzy based multi-objective approaches are widely used in supplier selection problem to deal with uncertain information [28,29]. This research answers the following questions: (i) Which supplier is selected based on the importance given to each objective? (ii) How much to purchase from each selected supplier?

The remainder of the paper is organized as follows. Section 2 provides related literature. Section 3 includes problem description and mathematical model for supplier selection and order allocation with a new resilience objective. Section 4 comprises of the proposed possibilistic fuzzy based solution methodology. Section 5 presents a numerical example to show the application of proposed supplier selection and order allocation model. Section 6 and Section 7 discussed the results of the proposed mathematical model and solution methodology. Finally, Section 8 presents some conclusions and future directions drawn from the study.

2. Literature Review

The word resilience first coined by Holling [30] in the context of ecology. According to Holling [30], the resilience is the ability of a system to absorb changes in state variables, driving variable and parameters, and still persist. Due to the increase in complexity and uncertainty in the business environment, several studies shown interest in the concept of resilience in a managerial perspective. Hamel and Valikangas [31] defined resilience as a capacity for continuous reconstruction. Sheffi [32] defined resilience in terms of enterprise resilience as the ability of an organization to successfully confront the unforeseen. Sutcliffe and Vogus [33] stated that resilience is the (1) ability to absorb strain and improve the functionality of organization despite the presence of difficulty or (2) ability to bounce back after disturbances. More recently, Woods [34] defined resilience in simple terms as system’s ability to bounce back after disruptions and to bounce forward through learning from those disruption events and increase the system’s adaptive capacity for handling uncertain events.

The concept of supply chain resilience gained prominent importance during recent years in supply chain risk management research [35]. Supply chain resilience is a relatively new concept to mitigate risks that can be defined as the ability to reduce the probability of a disruption, to reduce the impact of disruption, and to reduce the recovery time to normal performance [36]. Supply chain resilience has been defined by several authors in simpler and broader terms. Most of these studies described supply chain resilience as the ability to withstand disruptions and converge to the original state or to a new desirable state. Despite the increasing number of publication in supply chain resilience area, most of the researches provided qualitative insights and there is a limited number of quantitative modelling techniques available [37]. These qualitative models used different performance measures for designing resilient supply chains. Priya Datta, et al. [38] developed the framework to improve operational resilience. Falasca, Zobel and Cook [36] proposed three determinants (density, complexity, and node criticality) of supply chain resilience for supply chain design. Azevedo, et al. [39] proposed GResilient index to assess supply chain resilience using the Delphi technique. Miller-Hooks, et al. [40] proposed a transportation network resilience model using stochastic programming. They presented the expected fraction of demand fulfilment after disruption as resilience metric.

Supply chain usually functions in the system of parties where the resilience of one party (e.g. a supplier) is critical for overall supply chain resilience. As discuss earlier, many major disruptions break down supply networks and it takes a long time to recover. Whereas, the probability of disruption may be reduced by developing a resilient network and it takes considerably less recovery time [41]. Suppliers constitute the most important role in the performance of the supply chain, therefore, the resilience of the supply network is expected to contribute and increase overall supply chain resilience [42]. Despite its importance, there is very limited research conducted which consider resilience of supply network. The literature on resilient supplier selection is summarized in

Table 1. The classification of literature in resilient supplier selection.

Authors	Resilience Criteria	Types of Risks		Solution Methodology
		Network Risks	Operational Risks	Possibilistic Programming	Stochastic Programming	Fuzzy Programming	Grey Programming	Other
Wang, Herty and Zhao [6]	Multi-supplier, production capacity, product quality, production cost	-	✓	-	-	-	-	Fluid-dynamic models
Torabi, Baghersad and Mansouri [9]	Multiple sourcing, fortifying supplier, pre-positioned inventories, backup supplier, and supplier’s business continuity	-	✓	✓	✓	-	-	-
Parkouhi and Ghadikolaei [12]	Benefits, Opportunities, costs, and risks	-	✓	-	-	✓	✓	-
Hosseini and Al Khaled [13]	Absorptive capacity, adaptive capacity, and restorative capacity	✓	✓	-	-	-	-	Predictive analytics models
Sahu, Datta and Mahapatra [14]	Investment capacity, Responsiveness, and Inventory capacity	-	✓	-	-	✓	-	-
López and Ishizaka [15]	Flexibility, Visibility, Anticipation, Recovery, Security, Adaptability, Financial strength, Market position, and Collaboration	-	✓	-	-	✓	-	-
Sabouhi, Pishvaee and Jabalameli [16]	multi-sourcing, supplier fortification, and emergency inventory	-	✓	-	✓	✓	-	-
Jabbarzadeh, Fahimnia and Sabouhi [17]	Extra production capacities, Multiple sourcing, and Backup suppliers	-	✓	-	✓	✓	-	-
Hosseini and Barker [19]	Absorptive capacity, adaptive capacity, and restorative capacity	✓	✓	-	-	-	-	Bayesian network
Haldar, Ray, Banerjee and Ghosh [21]	Investment capacity, Responsiveness, and Emergency inventory holding capacity	-	✓	-	-	✓	-	-
Rajesh and Ravi [25]	Responsiveness, risk reduction, and Technical support	-	✓	-	-	-	✓	-
Parkouhi, Ghadikolaei and Lajimi [46]	Safety, Visibility, Environmental Controls, Trust, Flexibility, Support Services, Future Manufacturing Capabilities, and others	-	✓	-	-	-	✓	-
This paper	Supply density, Transit time, Resilience score of supplier’s locations	✓	✓	✓	-	✓	-	-

. Haldar, Ray, Banerjee and Ghosh [20] developed a resilient supplier selection model using AHP-QFD method to rate the suppliers. They used five resilience criteria (density, complexity, node criticality, responsiveness and re-engineering) for the supplier selection process. More recently, Haldar, Ray, Banerjee and Ghosh [21] proposed a fuzzy group decision-making approach for resilient supplier selection using triangular and trapezoidal fuzzy numbers. They considered investment, responsiveness, and capacity of holding inventory as resilient criteria. Sawik [22] proposed resilience strategies for supplier selection and order allocation problem under disruption risks. Recently, Sawik [23] proposed a stochastic mixed integer model for supplier selection and customer order scheduling using disruption risks for single and dual sourcing. Yilmaz-Börekçi, İşeri Say and Rofcanin [24] proposed a scale for measuring supplier resilience within supply networks. Mari, et al. [43] proposed a resilient and sustainable supply chain network model. They used expected disruption cost as a resilient metric. Torabi, Baghersad and Mansouri [9] developed a bi-objective stochastic model to trade-off the resilience level of the supply network and system cost. Memon, Lee and Mari [28] proposed supplier selection model using grey and uncertainty theories. Rajesh and Ravi [25] proposed a resilient supplier selection method using grey relational analysis. They consider supplier’s performance, responsiveness, risk reduction, technical capabilities and sustainability as selection criteria. Sahu, Datta and Mahapatra [14] proposed fuzzy-VIKOR based supplier selection framework by combining the general selection strategy and resilience strategy. They considered investment capacity, responsiveness, and inventory capacity as resilience criteria. Hosseini and Al Khaled [13] developed a resilience score of suppliers based on eight criteria namely backup supplier contracting, surplus inventory, location separation, robustness, reliability, reorganization, rerouting, and restoration. They proposed a hybrid ensemble and AHP approach to select a suitable supplier. Pramanik, et al. [44] developed a resilient supplier selection model using AHP-TOPSIS and Sen, et al. [45] proposed g-resilient supplier selection model using dominance-based fuzzy decision making. Wang, Zhang, Chong and Wang [18] developed a resilient supplier selection model for construction supply chain by integrating building information modelling and geographical information system. Parkouhi and Ghadikolaei [12] proposed FANP and grey-VIKOR based resilient supplier selection model. They considered benefits, opportunities, costs, and risks as resilience criteria. López and Ishizaka [15] proposed a coupled method of FCM and AHP for supplier selection for a resilient supply chain. Recently, Jabbarzadeh, Fahimnia and Sabouhi [17] proposed stochastic bi-objective optimization model considering fuzzy c-means to developed sustainable supply chain network that performs resiliently. Sabouhi, Pishvaee and Jabalameli [16] presented integrated stochastic and fuzzy DEA based model for supplier selection considering resilience. They considered multi-sourcing, supplier fortification, and emergency inventory as resilience criteria. More recently, Parkouhi, et al. [46] proposed a resilient supplier selection framework using grey-DEMATEL approach.

Most of the above studies considered operational risks for resilient supplier selection problem. This is the first time that a possibilistic fuzzy multi-objective model is proposed for resilient supplier selection and order allocation problem. Furthermore, for the first time in the literature supply density, resilience index score, and transit time are considered as supply selection criteria for resilient supply network.

Falasca, Zobel and Cook [36] proposed three characteristics (node criticality, supply chain complexity, and supply chain density) for building a resilient supply chain. Among them, supply chain density is the most important resilient criteria when designing supply networks under HILP disruptions. This is due to fact that denser supply networks are vulnerable to HILP disruption risks [36]. For example, 1999s Taiwan earthquake ended up having a significant effect on the entire global PC supply chain, because of the high concentration of computer component manufacturers in Hsinchu, Taiwan [47]. This example shows that the selection of a large number of suppliers from each region is a vulnerable multi-sourcing strategy, hence, this paper proposed supply density-based approach to tackle this problem. Furthermore, every country or territory has different resilient capabilities FMGlobal [48] and it affects the performance of the supply chain. Therefore, this study also proposed resilience index score-based criteria to supplier selection. The resilience index score is proposed by FMGlobal [48] is a data-driven tool to rank the countries to supply chain disruption risks. Nine key drivers of supply chain risks are considered and grouped into three categories namely: economic, risk quality, and supply chain factors. These nine drivers include local supplier quality, quality of fire risk management, GDP per capita, oil intensity, quality of hazard risk management, exposure to natural hazards, corruption control, infrastructure, and political risks [49]. Transit time is the last resilient criteria for supplier selection considered in this study. Transit time is an important indicator of supply chain flexibility [50]. Transit time reduction is one of the wildly used criteria to mitigate supply risks [51].

3. Problem Formulation

In this study, a garment manufacturer is assumed which want to select a suitable set of suppliers for the required material. All the model parameters are considered as fuzzy parameters. Five objectives considered in this study, namely: a cost which includes purchase and transportation costs, the rejection rate of suppliers, transit time from suppliers, supply density, and supplier resilience score based on their locations.

3.1. Mathematical Model Notations

• Indices
s existing suppliers	s = 1, 2, … S
l number of objective functions	l = 1, 2, … L
• Parameters
$\tilde{d}$	Total demand of required material
${\tilde{u}}_s$	Purchase cost of required material from supplier s
${\tilde{c}}_s$	Transportation cost of material from supplier s
${\tilde{\rho }}_s$	Percentage of the rejected material delivered by the supplier s
${\tilde{\nu }}_s$	Capacity of supplier s
${\tilde{k}}_s$	Minimum acceptable purchase quantity from supplier s
mxs	Maximum number of suppliers selected for required material
${\tilde{t}}_s$	Transit time from supplier s
dsab	Distance between selected supplier a ∈ s and selected supplier b ∈ s (a ≠ b)
rss	Resilience index score of supplier s
• Decision Variables
$q_s$	Purchase quantity from supplier s
${\mathsf{\Phi }}_s=\{\begin{array}{c}1\\ 0\end{array}$	1, If supplier s is selected, otherwise 0.
$\begin{array}{c}{\omega }_{ab}\\ {\omega }_{ab}^{\prime }\end{array}=\{\begin{array}{l}1, \mathrm{If} \mathrm{supplier} a \mathrm{and} \mathrm{supplier} b \mathrm{are} \mathrm{selected}.\\ 0,\mathrm{otherwise}\end{array}$	$\forall (a,b)\in s \mathrm{and} a\ne b$

3.2. Model Objectives

The Equation (1) represents the objective function for the cost. It is the sum of procurement cost and transportation cost from selected suppliers. The objective function of the total rate of rejection is estimated in Equation (2). The objective function for transit time from all selected supplier is calculated as shown in Equation (3). The Equation (4) shows the supply density for all selected suppliers. Total resilience index score objective is estimated in Equation (5). Where rs_s represents the resilience index score of supplier location obtained from FMGlobal [48].

$$\mathrm{Minimize} f_{cost}={\displaystyle \sum _s({\tilde{u}}_s+{\tilde{c}}_s )q_s}$$

(1)

$$\mathrm{Minimizes} f_{rej}={\displaystyle \sum _s{\tilde{\rho }}_s{q}_s}$$

(2)

$$\mathrm{Minimizes} f_{\mathit{time}}={\displaystyle \sum _s{\tilde{t}}_s{\mathsf{\Phi }}_s}$$

(3)

$$\mathrm{Maximize} f_{\mathit{den}}=\frac{1}{\tilde{d}}\left({\displaystyle \sum _{a\in s}{\displaystyle \sum _{\begin{array}{l}b\in s\\ a\ne b\end{array}}d{s}^{ab}{\mathsf{\Phi }}_s}}\right)$$

(4)

$$\mathrm{Maximize} f_{\mathit{res}}={\displaystyle \sum _s\frac{r{s}_s{q}_s}{\tilde{d}}}$$

(5)

3.3. Model Constraints

Constraint (6) ensures that total procured material should satisfy its demand.

$${\displaystyle \sum _s{q}_s}=\tilde{d}$$

(6)

Constraint (7) is capacity restrictions on the supplier. Also, it controls the flow between the supplier and the buyer through a binary variable. Constraint (8) ensures that purchase quantity from the selected supplier will be more than its acceptable order quantity limit. Constraint (9) restricts the maximum allowable supplier to be selected for the required material.

$$q_s\le {\tilde{\nu }}_s{\mathsf{\Phi }}_s$$

(7)

$$q_s\ge {\tilde{k}}_s{\mathsf{\Phi }}_s$$

(8)

$${\displaystyle \sum _s{\mathsf{\Phi }}_s}\le m{x}_s$$

(9)

Constraints (10) and (11) determine the intra-stage flow between suppliers and buyer. If buyer received material from both supplier a and supplier b then ${\omega }_{ab}={\mathsf{\Phi }}_{aϵs}={\mathsf{\Phi }}_{bϵs}=1$ and ${\omega }_{ab}^{\prime }=0$. On the contrary, if buyer not received material from both supplier a and supplier b then ${\omega }_{ab}^{\prime }=1$, ω_ab = 0 and ${\mathsf{\Phi }}_{aϵs}\ne {\mathsf{\Phi }}_{bϵs}$.

$$\begin{gathered} 2{\omega }_{ab}+{\omega }_{ab}^{\prime }={\mathsf{\Phi }}_{a\in s}+{\mathsf{\Phi }}_{b\in s} \\ \forall (a,b)\in s, \mathrm{and} a\ne b \end{gathered}$$

(10)

$$\begin{gathered} {\omega }_{ab}+{\omega }_{ab}^{\prime }\le 1 \\ \forall (a,b)\in s, \mathrm{and} a\ne b \end{gathered}$$

(11)

4. Fuzzy Based Solution Methodology

Fuzzy based programming methods are highly used for multi-objective optimization because of their capability in measuring and adjusting the decision maker’s satisfaction level of each objective function explicitly. In addition, the fuzzy theory is helpful to tackle the uncertain parameters related to supply chain optimization problem [52]. The main advantage of interactive fuzzy based approaches is that decision maker can efficiently achieve his/her preferences by controlling the search direction. The proposed solution methodology consists of the following steps.

Step 1: Convert uncertain mathematical model to equivalent auxiliary crisp

The proposed mathematical model is converted to an equivalent auxiliary crisp model. In this study, Jiménez, Arenas, Bilbao and Rodrı [51] approach is used which is based on an expected interval (EI) and expected value (EV) of fuzzy numbers. According to Jiménez, Arenas, Bilbao and Rodrı [51], the EI and EV of triangular fuzzy number (TFN) can be defined as in equation (12) and (13) respectively. Where ϑ^pes is the pessimistic value, ϑ^mos is the most likely value, and ϑ^opt is the optimimum value of triangular fuzzy number (ϑ). This research considered TFN because it is frequently used for a practical purpose [53].

$$EI(\tilde{\vartheta })=[E_1^{\vartheta },E_2^{\vartheta }]=\left[{\displaystyle \underset{0}{\overset{1}{\int }}{f}_{\vartheta }^{-1}(x)dx,{\displaystyle \underset{0}{\overset{1}{\int }}{g}_{\vartheta }^{-1}(x)dx,}}\right]=\left[\frac{1}{2}({\vartheta }^{pes}+{\vartheta }^{mos}),\frac{1}{2}({\vartheta }^{mos}+{\vartheta }^{opt})\right]$$

(12)

$$EV(\tilde{\vartheta })=\frac{E_1^{\vartheta }+E_1^{\vartheta }}{2}=\frac{{\vartheta }^{pes}+2{\vartheta }^{mos}+{\vartheta }^{opt}}{4}$$

(13)

Using the above Equations (12) and (13), the equivalent auxiliary crisp model can be formulated as follows.

$$\mathrm{Minimize} f_{\mathit{cost}}={\displaystyle \sum _s\left(\frac{u_s^{pes}+2u_s^{mos}+u_s^{opt}+c_s^{pes}+2c_s^{mos}+c_s^{opt}}{4}\right)q_s}$$

(14)

$$\mathrm{Minimizes} f_{\mathit{rej}}={\displaystyle \sum _s\left(\frac{{\rho }_s^{pes}+2{\rho }_s^{mos}+{\rho }_s^{opt}}{4}\right)q_s}$$

(15)

$$\mathrm{Minimizes} f_{\mathit{time}}={\displaystyle \sum _s\left(\frac{t_s^{pes}+2t_s^{mos}+t_s^{opt}}{4}\right){\mathsf{\Phi }}_s}$$

(16)

$$\mathrm{Maximize} f_{\mathit{den}}=\frac{1}{\frac{d^{pes}+2d^{mos}+d^{opt}}{4}}\left({\displaystyle \sum _{a\in s}{\displaystyle \sum _{\begin{array}{l}b\in s\\ a\ne b\end{array}}d{s}^{ab}{\mathsf{\Phi }}_s}}\right)$$

(17)

$$\mathrm{Maximize} f_{\mathit{res}}={\displaystyle \sum _s\frac{r{s}_s{q}_s}{\frac{d^{pes}+2d^{mos}+d^{opt}}{4}}}$$

(18)

Subject to

$${\displaystyle \sum _s{q}_s}\ge \left[\frac{\alpha }{2}\left(\frac{d^{mos}+d^{opt}}{2}\right)+(1-\frac{\alpha }{2})\left(\frac{d^{pes}+d^{mos}}{2}\right)\right]$$

(19)

$${\displaystyle \sum _s{q}_s}\le \left[\frac{\alpha }{2}\left(\frac{d^{pes}+d^{mos}}{2}\right)+(1-\frac{\alpha }{2})\left(\frac{d^{mos}+d^{opt}}{2}\right)\right]$$

(20)

$$q_s\le {\mathsf{\Phi }}_s\left[\alpha \left(\frac{{\nu }_s^{pes}+{\nu }_s^{mos}}{2}\right)+(1-\alpha )\left(\frac{{\nu }_s^{mos}+{\nu }_s^{opt}}{2}\right)\right]$$

(21)

$$q_s\ge {\mathsf{\Phi }}_s\left[\alpha \left(\frac{k^{mos}+k^{opt}}{2}\right)+(1-\alpha )\left(\frac{k^{pes}+k^{mos}}{2}\right)\right]$$

(22)

$${\displaystyle \sum _s{\mathsf{\Phi }}_s\le m{x}_s}$$

(23)

$$2{\omega }_{ab}+{\omega }_{ab}^{\prime }={\mathsf{\Phi }}_{a\in s}+{\mathsf{\Phi }}_{b\in s}$$

(24)

$${\omega }_{ab}+{\omega }_{ab}^{\prime }\le 1$$

(25)

Step 2: Determine α- extreme solutions

To estimate the upper (α-UB) and lower (α-LB) bounds to each objective, the crisp model developed in Step 1 is solved for each objective along with its constraint.

Step 3: Determine fuzzy membership function

Develop the fuzzy membership function for each objective using lower (α-LB) and upper (α-UB) bound values. The linear memberships for fuzzy goals are given as follows. It is assumed that membership functions are linear based on preferences and satisfaction level.

$${\mu }_{\mathit{cost}}(x)=\{\begin{array}{ll}1,& f_{\mathit{cost}}(x)\le f_{\mathit{cost}}^{\alpha -LB}\\ \frac{f_{\mathit{cost}}^{\alpha -UB}-f_{\mathit{cost}}(x)}{f_{\mathit{cost}}^{\alpha -UB}-f_{\mathit{cost}}^{\alpha -LB}},& f_{\mathit{cost}}^{\alpha -LB}< f_{\mathit{cost}}(x)\le f_{\mathit{cost}}^{\alpha -UB} \\ 0,& f_{\mathit{cost}}(x)\ge f_{\mathit{cost}}^{\alpha -UB}\end{array}$$

(26)

$${\mu }_{\mathit{rej}}(x)=\{\begin{array}{ll}1,& f_{rej}(x)\le f_{rej}^{\alpha -LB}\\ \frac{f_{rej}^{\alpha -UB}-f_{rej}(x)}{f_{rej}^{\alpha -UB}-f_{rej}^{\alpha -LB}},& f_{rej}^{\alpha -LB}\le f_{rej}(x)\le f_{rej}^{\alpha -UB} \\ 0,& f_{rej}(x)\ge f_{rej}^{\alpha -UB}\end{array}$$

(27)

$${\mu }_{time}(x)=\{\begin{array}{ll}1,& f_{time}(x)\le f_{time}^{\alpha -LB}\\ \frac{f_{time}^{\alpha -UB}-f_{time}(x)}{f_{time}^{\alpha -UB}-f_{time}^{\alpha -LB}},& f_{time}^{\alpha -LB}\le f_{time}(x)\le f_{time}^{\alpha -UB} \\ 0,& f_{time}(x)\ge f_{time}^{\alpha -UB}\end{array}$$

(28)

$${\mu }_{den}(x)=\{\begin{array}{ll}1,& f_{den}(x)\ge f_{den}^{\alpha -UB}\\ \frac{f_{den}(x)-f_{den}^{\alpha -LB}}{f_{den}^{\alpha -UB}-f_{den}^{\alpha -LB}},& f_{den}^{\alpha -UB}\le f_{den}(x)\le f_{den}^{\alpha -LB}\\ 0,& f_{den}(x)\le f_{den}^{\alpha -LB}\end{array}$$

(29)

$${\mu }_{res}(x)=\{\begin{array}{ll}1,& f_{res}(x)\ge f_{res}^{\alpha -UB}\\ \frac{f_{res}(x)-f_{res}^{\alpha -LB}}{f_{res}^{\alpha -UB}-f_{res}^{\alpha -LB}},& f_{res}^{\alpha -UB}\le f_{res}(x)\le f_{res}^{\alpha -LB} \\ 0,& f_{res}(x)\le f_{res}^{\alpha -LB}\end{array}$$

(30)

where $f_l^{\alpha -LB}$ is a minimum value of $f_l\left(x\right)$ and $f_l^{\alpha -UB}$ is a maximum value of $f_l\left(x\right)$ with predefined α. These values the l^th objective depends on its nature. $f_l^{\alpha -LB}$ are set as the aspiration level of cost, rejection rate, and transit time objective. Whereas $f_l^{\alpha -UB}$ are set as the aspiration level of supply density and resilience index.

Step 4: Convert the multi-objective model into a single objective

The proposed model is converted into the single objective in this stage. In this paper, two most popular fuzzy based approaches i.e., Tiwari, Dharmar and Rao [26] weighted additive approach and Werners [27] fuzzy and operator are implemented.

The weighted additive approach allows the buyer to assign different weights to objectives in the simple additive fuzzy achievement function. Mathematical formulation of the weighted additive method is as follows.

$$\begin{array}{lll}\mathrm{maximize}& {\displaystyle \sum _l{w}_l}{\mu }_{zl}(x)& \\ \mathrm{subject} \mathrm{to}& {\mu }_{zl}(x)\in [0,1],& \forall l\\ & x\ge 0& \end{array}\}$$

(31)

where w_l represents the weight of i^th objective. Selection of weights are subjected choice of decision makers and some good techniques can be used to determine the weights such as FAHP and structural equation modelling.

Werners’ fuzzy and operator method are widely used interactive method. The advantage of this method is that it is positively related to the compensation rate due to its strong monotonicity. Additionally, it is easy to handle and has generated reasonable consistent results in applications [54]. By adopting the Werner’s’ method following a single objective model can be formed.

$$\begin{array}{lll}\mathrm{maximize}& \mathsf{\gamma }{\mathsf{\zeta }}_0+(1-\mathsf{\gamma }){\displaystyle \sum _l{\mathsf{\zeta }}_1}& \\ \mathrm{subject} \mathrm{to}& {\mathsf{\mu }}_1(x)\ge {\mathsf{\zeta }}_0+{\mathsf{\zeta }}_1,& \forall l\\ & {\mathsf{\zeta }}_0,{\mathsf{\zeta }}_1,\mathsf{\gamma }\in [0,1]& \end{array}\}$$

(32)

where, ζ₁ is the difference between satisfaction level of objectives their minimum satisfaction level ζ₀, That is, ζ₁ = μ₁ − ζ₀. γ denotes the coefficient of compensation.

Step 5: Determine the solution method parameter

Determine the values of the relative importance of objectives (w_l) and coefficient of compensation (γ) to solve the mathematical model using both weighted additive approach and Werners’ fuzzy and operator methods.

Step 6: Solve the model

In the last step, solve the model by using parameters of the mathematical model and solution methods. This process continues by varying the solution method parameters (i.e., γ and w_l) until decision makers are satisfied with the final solution. If decision makers want to modify the value of α, then restart the process from step 2.

5. An Illustration

The effectiveness of the proposed resilient supplier selection model and solution methodology is demonstrated in this section. The data relates to a realistic situation of a garment manufacturing sector as shown in Figure 1. The adopted situation can easily be extended to any other industry. Initially, the data are estimated as most likely values, these most likely values of the fuzzy parameter (fz^mos) are estimated using available information or hypothetically set based on realistic assumption. The pessimistic and optimistic values are estimated using $f{z}^{pes}=(1-{\eta }_1)f{z}^{mos}$ and, where two random numbers η₁ and η₂ are assumed between 0.2 and 0.8 to estimate.

Table 2. Model input data.

Potential Supplier Location	Purchase Cost of Material ($/unit)	Transportation Cost of Material ($/unit)	Percentage of the Rejected Material	Capacity of Suppliers (1000 units)	Resilience Score
Korea	(6,8,10)	(0.08,0.13,0.21)	(0.01,0.01,0.02)	(3.08,5,8.08)	42.1
China	(1,2,4)	(0.11,0.19,0.30)	(0.04,0.06,0.10)	(3.7,6,9.7)	45.3
Thailand	(4,6,8)	(0.05,0.09,0.14)	(0.02,0.03,0.05)	(3.08,5,8.08)	39
Bangladesh	(1,3,5)	(0.11,0.17,0.28)	(0.02,0.03,0.05)	(3.7,6,9.7)	29
India-(Calcutta)	(3,5,7)	(0.10,0.16,0.26)	(0.01,0.02,0.03)	(2.46,4,6.46)	27.1
India-(Hyderabad)	(3,4,5)	(0.01,0.02,0.03)	(0.01,0.02,0.03)	(4.62,7.5,12.12)	27.1
Pakistan	(3,5,7)	(0.06,0.10,0.16)	(0.01,0.02,0.03)	(3.08,5,8.08)	22.2
Turkey	(8,10,12)	(0.08,0.12,0.20)	(0.01,0.01,0.02)	(3.08,5,8.08)	38.4
Brazil	(6,8,10)	(0.09,0.14,0.22)	(0.02,0.04,0.06)	(3.7,6,9.7)	47.8
Mexico	(7,9,11)	(0.11,0.17,0.28)	(0.02,0.03,0.05)	(4.31,7,11.31)	44.8

shows the data set of unit purchase cost, transportation cost, rejection rate, capacity, and resilience index score of potential suppliers. The unit purchase costs from each supplier are hypothetical set based on labour cost, land value, resource availability at supplier locations. Transportation cost and transit time (see

Table 3. Transit time from suppliers.

Supplier Location	Transit Time (days)
Korea	(7.39,12,19.39)
China	(9.24,15,24.24)
Thailand	(1.23,2,3.23)
Bangladesh	(1.85,3,4.85)
India-(Calcutta)	(1.23,2,3.23)
India-(Hyderabad)	(0.62,1,1.62)
Pakistan	(3.08,5,8.08)
Turkey	(8.01,13,21.01)
Brazil	(14.78,24,38.78)
Mexico	(19.71,32,51.71)

) from each supplier to manufacturer are estimated from sea rates (https://www.searates.com/). The resilience index score of each supplier is estimated from FM, Global resilience index data-driven tool (https://www.fmglobal.com/) based on the location of suppliers. The distance (as the crow flies) between suppliers (see

Table 4. The distance between suppliers (kilometers).

Suppliers’ location.	Korea	China	Thailand	Bangladesh	India (Calcutta)	India (Hyderabad)	Pakistan	Turkey	Brazil	Mexico
Korea	N/A	1028.16	3723.83	3783.61	4037.39	5212.7	5768.78	7641.12	16,739.92	12,029.52
China	1028.16	N/A	2743.18	3059.41	3304.68	4478.2	5273.64	7587.05	17,589.43	13,051.91
Thailand	3723.83	2743.18	N/A	1526.77	1616.16	2388.26	3709.84	6926.22	17,348.59	15,721.92
Bangladesh	3783.61	3059.41	1526.77	N/A	250.06	1432.18	2374.24	5408.55	16,067.89	15,092.77
India (Calcutta)	4037.39	3304.68	1616.16	250.06	N/A	1180.87	2186.91	5327.76	15,902.02	15,299.76
India (Hyderabad)	5212.7	4478.2	2388.26	1432.18	1180.87	N/A	1461.88	4788.23	14,979.22	15,878.77
Pakistan	5768.78	5273.64	3709.84	2374.24	2186.91	1461.88	N/A	3334.57	13,710.04	14,838.79
Turkey	7641.12	7587.05	6926.22	5408.55	5327.76	4788.23	3334.57	N/A	10,723.38	11,983.21
Brazil	16,739.92	17,589.43	17,348.59	16,067.89	15,902.02	14,979.22	13,710.04	10,723.38	N/A	5692.39
Mexico	12,029.52	13,051.91	15,721.92	15,092.77	15,299.76	15878.77	14,838.79	11,983.21	5692.39	N/A

) are calculated using the Google maps (https://www.google.com/maps). The demand for raw materials is assumed as most likely 8000 units, minimum acceptable order quantity is assumed as 1000 units, and maximum three suppliers can be selected for required material.

6. Model Solution and Result Analysis

The proposed fuzzy multi-objective model is solved in Lingo 14.0. Lingo optimization software has been widely used in supply chain optimization problems [55,56]. According to the steps of the proposed methodology, the payoff values are estimated by solving the single objective model.

Table 5. Payoff values.

Objective	Cost ($)	Rejection (%)	Transit Time (days)	Supply Density	Resilience Score
Minimize Cost	22,514.93	0.097	19.04	0.022	38.2
Minimize Rejection	75,796.12	0.025	26.45	0.056	39.22
Minimize Transit time	38,143.47	0.052	3.17	0.017	29.79
Maximize Supply density	70,159.31	0.1	61.11	0.286	43.41
Maximize Resilience score	50,560.14	0.1	41	0.12	47.89

shows payoff values estimated from the single objective model. The aspiration level of cost, rejection, transit time, supply density, and resilience score are estimated as $22,514.93, 0.025, 3.17 days, 0.286, and 47.89 respectively. Once the payoff values are estimated, the fuzzy membership functions of objectives are estimated as shown below.

$${\mu }_{\mathit{cost}}(q_s)=\{\begin{array}{l}1,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{\mathit{cost}}(q_s)\le 22514.93\\ \frac{75796.12-f_{\mathit{cost}}(q_s)}{75796.12-22514.93},\hspace{1em}22514.93< f_{\mathit{cost}}(q_s)\le 75796.12 \\ 0,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{\mathit{cost}}(q_s)\ge 75796.12\end{array}$$

$${\mu }_{\mathit{rej}}(q_s)=\{\begin{array}{l}1,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{rej}(q_s)\le 0.025\\ \frac{0.10-f_{rej}(q_s)}{0.10-0.025},\hspace{1em}0.025\le f_{rej}(q_s)\le 0.10 \\ 0,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{rej}(q_s)\ge 0.10\end{array}$$

$${\mu }_{time}(q_s)=\{\begin{array}{l}1,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{time}(q_s)\le 3.17\\ \frac{61.11-f_{time}(q_s)}{61.11-3.17},\hspace{1em}3.17\le f_{time}(q_s)\le 61.11\\ 0,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{time}(q_s)\ge 61.11\end{array}$$

$${\mu }_{den}(q_s)=\{\begin{array}{l}1,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{den}(q_s)\ge 0.28\\ \frac{f_{den}(q_s)-0.017}{0.28-0.017},\hspace{1em}0.017\le f_{den}(q_s)\le 0.28\\ 0,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{den}(q_s)\le 0.017\end{array}$$

$${\mu }_{res}(q_s)=\{\begin{array}{l}1,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{res}(q_s)\ge 47.89\\ \frac{f_{res}(q_s)-29.79}{47.89-29.79},\hspace{1em}29.79\le f_{res}(q_s)\le 47.89\\ 0,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{f}_{res}(q_s)\le 29.79\end{array}$$

The proposed mathematical model is solved using both methods, that is, a weighted additive approach and Werners’ 'fuzzy and' operator methods. A solution of illustrated case example using both methods is shown in

Table 6. Solution of a case example.

Model Objectives	Weighted Additive Approach a		Werners’ 'Fuzzy and' Operator b
	μl (x)	fl (x)	μl (x)	fl (x)
Cost	0.346	78,238.03	0.195	65,374.99
Rejection	0.051	0.12	0.366	0.072
Transit time	0.457	42.0	0.383	38.89
Supply density	0.975	0.27	0.972	0.272
Resilience score	0.786	75.5	0.803	44.33

α = 0.9, ^a w₁ = 0.2, w₂ = 0.1, w₃ = 0.2, w₄ = 0.2, w₅ = 0.3, ^b γ = 0.7.

. The result shows that both methods produce a comprehensive optimal solution. However, the weighted additive approach considers the importance of objectives based on the weight given to each objective. For example, when more importance is given to resilience score (i.e., w₅ = 0.3), it results in 78.6% achievement of aspiration level. On the other hand, when all objectives and equal importance (Werners’ method) than resilience score is given, the objective results in 34% achievement of aspiration level. Comparative analysis of both methods is shown in Figure 2.

7. Sensitivity Analysis

Sensitivity analysis of both methods is carried out by varying methodology parameters (i.e., α, γ and w_l) as discussed in the final step of the proposed solution methodology.

Table 7. Sensitivity analysis of the weighted additive approach.

Alpha	Objective Weights	μcost	μrej	μtime	μden	μres	fcost	frej	ftime	fden	fres
0.9	w1 = 0.2, w2 = 0.1, w3 = 0.2, w4 = 0.2, w5 = 0.3	0.184	0.200	0.365	0.97	0.881	65,941.28	0.085	39.9	0.27	45.7
	w1 = 0.3, w2 = 0.2, w3 = 0.2, w4 = 0.1, w5 = 0.2	0.920	0.200	0.762	0.061	0.498	26,730.90	0.085	16.9	0.03	38.8
0.6	w1 = 0.2, w2 = 0.1, w3 = 0.2, w4 = 0.2, w5 = 0.3	0.111	0.100	0.566	0.910	0.935	66,577.34	0.092	28.3	0.25	47.4
	w1 = 0.3, w2 = 0.2, w3 = 0.2, w4 = 0.1, w5 = 0.2	0.851	0.200	0.762	0.061	0.634	27,727.65	0.085	16.9	0.03	42.6

and

Table 8. Sensitivity analysis of Werners’ ‘fuzzy and’ operator.

α	γ	μcost	μrej	μtime	μden	μres	fcost	frej	ftime	fden	fres
0.6	0.0–0.8	0.107	0.366	0.383	0.972	0.793	66,774.85	02.072	38.8	0.27	46.4
	0.9	0.366	0.366	0.383	0.972	0.451	53,167.88	0.072	38.89	0.27	38.3
	1.0	0.407	0.566	0.543	0.407	0.406	51,054.51	0.057	29.62	0.12	37.3
0.9	0.0–0.8	0.195	0.366	0.383	0.972	0.803	65,374.99	0.072	38.89	0.27	44.33
	0.9	0.366	0.366	0.383	0.972	0.541	56,259.06	0.072	38.89	0.27	39.59
	1.0	0.392	0.566	0.543	0.407	0.392	54,899.34	0.057	29.6	0.12	36.88

show the analysis result of the weighted additive approach and Werners’ ‘fuzzy and’ operator method, respectively. Two cases are solved with the weighted additive approach: (1) more importance is given to resilience objective and (2) more importance given to cost objective.

The result shows that there is a tradeoff in economic objective and resilience objective. Hence, it can be said that the increase is supply chain resilience will tend to increase the total cost of the supply network. This is possible because economic supply networks are more dense networks in order to minimize transportation cost, hence any HILP disruption event such as earthquake or tsunami may disrupt more than one supplier. Therefore, it is suggested that companies should avoid denser supply networks to minimize risks from high impact low probability disruption event. Analysis result shows that 92% achievement of cost aspiration level will bring achievement of about only 6.1% supply density and 49.8% resilience index score aspiration levels respectively. On the other hand, 18.4% achievement of cost aspiration level will bring achievement of about 97% supply density and 88.1% resilience index score respectively.

Assume that decision makers required balance results and choose the best outcome at α = 0.9, and γ = 1.0 (highlighted row in Table 8). Hence, the final decision of the illustrated example is shown in Figure 3. The optimal quantity of material purchase from suppliers are q (Korea) = 3751 units, q (Bangladesh) = 3655 units, and q (Turkey) = 1258 units.

8. Conclusion and Future Suggestions

This paper introduced a novel possibilistic fuzzy environment for resilient supplier selection and order allocation considering new resilience criteria (supply density, transit time and resilience index score) for the selection of suppliers. For this purpose, an interactive fuzzy multi-objective programming approach is introduced to reduce uncertainties inherent in the supplier selection decision. A six-step solution methodology is designed to solve the proposed uncertain multi-objective model and a numerical case example is provided to show the applicability of the proposed model in a real situation.

This study significantly helps the practitioners who are trying to consider resilience in their supply network. The research results show the importance of supplier’s location in order to minimize disruption risks. Furthermore, the proposed model will help the managers to effectively evaluate and select a suitable set of suppliers while considering cost and resilience simultaneously. Moreover, the proposed possibilistic fuzzy based solution methodology will be helpful for practitioners and academicians to tackle cognitive and stochastic uncertainties related to supplier evaluation and selection problem. This study also helps the academicians to analyze the importance of resilient supply networks under disruption risks. Additionally, the proposed multi-objective possibilistic fuzzy-based approach can be useful in another area of supply chain optimization.

Although this research gained important insights from the implementation of proposed resilience criteria and solution methodology, there are some limitations which may be considered in future research. This paper assumed triangular fuzzy numbers for model parameters, it will be interesting to compare the research results with other fuzzy numbers such as trapezoidal in future studies. Furthermore, this research only considered a disruption of the location where suppliers are located. However, it will be important to consider the disruption of transport links between suppliers and buyers. Another interesting possible direction for future research is to extend the proposed model to multi-commodity and multi-period planning horizon. Lastly, the analysis result indicates that economic networks are denser network and they may be vulnerable to disruption risks, it will be valuable to further investigate the relationship between the economic network (denser network) and disruption risks.