A New Fuzzy TOPSIS-TODIM Hybrid Method for Green Supplier Selection Using Fuzzy Time Function

Today green supply chain is considered all around the world and supplier selection has been changed regarding these green and carbon emission criteria, so green supplier selection has been a major problem in this area. In this study we use fuzzy time function to assist managers in green supplier selection under uncertainty and ambiguity. This function will consider derivation from the goal during the time and by using it, and we will be able to have the best supplier in every period after having some modification in legal limitations for green supplier selection criteria. We use a fuzzy TOPSIS to have better initial weighting in TODIM, a discrete multicriteria method based on prospect theory in uncertainty (known as TODIM in Portuguese) decision making method. The results indicated that our proposed approach can easily and effectively accommodate criteria with gains and loss functions during time and also by using this method we will have a more reasonable predict of our suppliers ranking in future and that will help us in future investment in these suppliers. Finally it has been shown in car industries in Iran.


Introduction
In recent years, the European Union (EU) has established various environmental policies, including the RoHS (Restricted Use of Hazardous Substances in Electronics and Electrical Equipment) as well as WEEE (Waste Electronics and Electrical Equipment) directives.So far, environmental management has evolved to include boundary-spanning activities in the upstream and the downstream supply chains.Sirvastava defined green supply chain management (GSCM) as a combination of environmental and supply chain management activities, including product design, material selection, manufacturing process, final product delivery, and end-of-life product management.With GSCM, firms can select from a wide variety of suppliers and leverage resources throughout the firm to eliminate the environmental impact of supply chain activities, Tseng [1].
Firms typically expect their suppliers to surpass environmental compliance and to develop efficient and green product design.In addition, suppliers are expected to assess the life cycle of a product.Although the qualitative criteria are littered with subjective perception because the GSCM evolution criteria tend to be subjective, qualitative, or described with linguistic information.Thus, it is extremely difficult for the decision makers to express their preference using exact numerical values, Zhang et al. [2], so it is important to use linguistic number to distinguish between these criteria.Nevertheless, a firm's suppliers must satisfy green criteria under the constrains of subjective human preferences, which include uncertainty, in addition to uncertainties of time and planning for future; in this study we consider time in all decision making processes and this will be more effective in combination with fuzzy number for linguistic terms.Most of studies assess GSCM based on its alignment with the firm's identified objectives and fulfillment of a set of assessment criteria.In general, the evaluation of criteria is highly subjective and unstructured, as it relies significantly on manager's experience, knowledge, and intuition.However, managers cannot consider all relevant criteria due to bounded rationality and limited capacity for information processing; see Tseng [1].The fuzzy TODIM approach can easily and effectively accommodate criteria with gains and loss functions; see Tseng et al. [3], but with all advantages of using fuzzy TODIM two problems have been remained: how can we predict our supplier behavior and how much time this supplier can be effective in our system and when we must change our supplier.Previous studies also used MADM tools for selecting suppliers according to criteria which was made it easier with using fuzzy approach.Izadikhah [4] used TOPSIS method with interval-valued fuzzy number in supplier selection and selected best supplier according to criteria; in our study, we used TODIM fuzzy because of its sensitivity to evaluate data and show how much we can improve our system with changing supplier, also we used fuzzy time function to show changes in other suppliers, and finally we can improve our supply chain with respect to this data and estimation of suppliers action.Hence, the evolution approach often implemented ineffectively because management does not effectively use appropriate knowledge and experience on GSCM.Due to the impact of this lack, management is not confident that supplier selection has been studied and applied to a set of criteria to maximize firms' green supply chain benefits.How could management apply knowledge of previously (successful and unsuccessful) supplier selections to support future decision making?The objective of this study is to create a mechanism that could assist managers in analyzing and selecting green suppliers during time horizons and this will help us to predict a picture for future of the organization activities; for these targets we use a hybrid fuzzy TOPSIS-TODIM decision making process considering fuzzy time function.We will continue this paper by literature review on fuzzy time function and fuzzy TOPSIS and also TODIM in Section 2 and also have some necessary explanations about the above techniques in Section 3; we show our proposed method in this paper by using an illustrative example for green supplier selection.Finally in Section 4 conclusions will be presented.

Fuzzy Logic.
A fuzzy set is a class of objects with grades of membership.A membership function is between zero and one, Zadeh [5].Fuzzy logic is derived from fuzzy set theory to deal with reasoning that is approximate rather than precise.It allows the model to easily incorporate various subject experts' advice in developing critical parameter estimates, Zimmermann [6].In other words, fuzzy logic enables us to handle uncertainty.There are some kinds of fuzzy numbers.Among the various shapes of fuzzy number, triangular fuzzy number (TFN) is the most popular one.It is represented with three points as  = ( 1 ,  2 ,  3 ).The membership function is illustrated in (1).Let  and  be defined as  = ( 1 ,  2 ,  3 ),  = ( 1 ,  2 ,  3 ).
is the subtraction of them: (1)

Fuzzy Time Function (FTF).
Fuzzy time function is an approach to considering time as an important factor in uncertainty.In many situations we have to change our approach because of the uncertainty and changes in criteria during time.Sometimes these kinds of changes cause extra cost for organization.With fuzzy time function we could predict when some new factors like technology will be used in our system and how they could affect our decisions.de Figueiredo and Perkusich, [7] used fuzzy time function in fault and timing analysis in real time, where fuzzy time function used to verify the possibility of the input tokens in places.Yoneyama [8] used fuzzy set theory for time-delay systems that considered the uncertain Takgai-Sugeno fuzzy model with time delay.In this paper we used a new aspect of fuzzy time function which represents a combination of time and fuzzy triangular numbers.
A fuzzy time function consists of three sections, optimistic, normal, and pessimistic related to time and these lines do not cross each other's (Figure 1).
In fact, this function in each time represents a fuzzy triangular number (3), which considers all three possibilities.In many organizations some methods or some criteria will be out of date after some years or new equipment or method will be better than the older ones and these changes will effect in supplier selection methods and factors.By using FTF we could consider them in our decisions and can also upgrade our systems based on expert's viewpoints.Considering the above descriptions we will have FTF for every criterion as 1 for optimistic  2 for normal  3 for pessimistic. ( According to (2) FTF represents the triangular fuzzy number during the time: For computing the FTF in this paper we introduce two ways.
First.In this case we use average of slopes when the deviation between slopes of criteria is low and is not very important and after this we use (2) for these criteria.
To get precise decisions, it is recommended to use the separate functions for each period of time.
Second.In this method we do not use mean of the slopes because of more deviation in numbers and in these cases we use separate functions in calculations that will be more reasonable.
In some cases which fuzzy time function has more complicated function, we must use the combination of these two methods.
2.3.Fuzzy TOPSIS.TOPSIS, one of the classical multicriteria decision making methods, was developed by Hwang and Yoon [9].It is based on the concept that the selected alternative should have the shortest distance from the positive ideal solution (PIS) and the farthest from the negative ideal solution (NIS).TOPSIS also provides an easily understandable and programmable calculation procedure.It has the ability of taking various criteria with different units into account simultaneously; see Ekmekc ¸ioĝlu et al. [10].A number of fuzzy TOPSIS methods have been developed in recent years.Chen and Hwang [11] first applied fuzzy numbers to establish fuzzy TOPSIS method in which relative closeness for each alternative is evaluated based on fuzzy arithmetic operations.Chen [12] extends the TOPSIS method to fuzzy group decision making situations by considering triangular fuzzy numbers and defining crisp Euclidean between two fuzzy numbers.Chu [13] and Chu and Lin [14] further improved the methodology proposed by Chen [12].Jahanshahloo et al. [15] and Chu and Lin [16] extended the fuzzy TOPSIS method based on alpha level sets with interval arithmetic.Fuzzy TOPSIS has been introduced for various multiattribute decision-making problems.Yong [17] used fuzzy TOPSIS for planet location selection, Chen et al [18] used fuzzy TOPSIS for supplier selection, and Kahraman et al. [19] used fuzzy TOPSIS to select municipal solid waste disposal method and site.Kutlu and Ekmekc ¸ioĝlu [20] used a modified fuzzy TOPSIS integrated with fuzzy AHP to propose new FMEA failure modes and effects analysis which overcomes the shortcomings of traditional FMEA.Kaya and Kahraman [21] proposed a modified fuzzy TOPSIS for the best energy technology selection.
In the following, Chen's fuzzy TOPSIS method is explained.Chen [12] extends the TOPSIS method to fuzzy group decision making situations by considering triangular fuzzy numbers and defining crisp Euclidean distance between two fuzzy numbers.In Chen's fuzzy TOPSIS, linguistic preferences can easily be converted to fuzzy numbers which
In this study the linguistic variables used for each criterion and we will make FTF for each criterion during times.At the beginning weights of the criteria and fuzzy ratings of alternatives with respect to each criterion have been calculated, and the fuzzy multicriteria of decision-making problem can be expressed in matrix format as where   is the rating of the alternative   according to criterion  (i.e.,   ) and   denotes the importance weight of   .These linguistic variables can be described by triangular fuzzy numbers:   = (  ,   ,   ).To avoid the complicated normalization formula we used in classical TOPSIS, the linear scale transformation is used here to transform the various criteria scales into a comparable scale.Therefore, we can obtain the normalized fuzzy decision matrix denoted by : where  and  are the set of benefit criteria and cost criteria, respectively, and The above normalization method preserves the ranges of normalized fuzzy decision matrix as where According to the weighted normalized fuzzy decision matrix (  ), we know that the element V  is positive normalized triangular fuzzy numbers and will be in [0, 1].Then, we can define the fuzzy positive-ideal solution (FPIS,  * ) and fuzzy negative-ideal solution (FPIS;  − ) as where The distance of each alternative from  * and  − can be currently calculated as where (⋅) is the distance between two fuzzy numbers calculated as follows: where  = ( 1 ,  2 ,  3 ) and  = ( 1 ,  2 ,  3 ) are two triangular fuzzy numbers.The closeness coefficient of each alternative is calculated as Obviously, when   is near to 1 alternative,   is closer to the (FPIS,  * ) and farther from (FPIS,  − ).Therefore, according to the closeness coefficient, we can determine the ranking order of all alternatives and select the best one among a set of feasible alternatives.

TODIM
2.4.1.TODIM Method.TODIM is a discrete multicriteria method founded on prospect theory.The TODIM method has been successfully used and empirically validated in different applications.This is an experimental method based on how people make effective decisions in risky conditions.The shape of the value function of TODIM is identical to prospect theory's gain and loss function.The global multicriteria value function of TODIM aggregates all measures of gains and losses by considering all criteria.Gomes and Rangel [23] apply TODIM to investigate and recommend options for upstream projects for the natural gas reserves recently discovered in the Mexilho field in the Santos Basin, Brazil.In addition, Gomes and Rangel presented an evaluation of residential properties with real estate agents in Brazil and define a reference value for the rents of these properties' characteristics using the TODIM method for multicriteria decisions.This approach can assist professionals in the real estate market to evaluate alternatives clearly using the criteria defined by specialists.In general, TODIM can be used for qualitative and quantitative criteria.Verbal scales of qualitative criteria are converted into cardinal scales, and both types of scales are normalized.The relative measure of dominance of one alternative over another alternative is determined for each pair of alternatives.This measure is computed as the sum of all criteria of relative gain and loss values for these alternatives.This sum will be a gain, a loss, or zero, depending on the performance of each alternative with respect to each criterion.Tseng et al. [3] apply TODIM and TFNs to select green supplier chain.They used TODIM to find the best supplier, but the first weighting for TODIM with using fuzzy set theory does not consider the relationship between the criteria more obvious than fuzzy TOPSIS, with combination of these two approaches the TODIM input will be more accurate and the result will be more correct consequently.We use a new combination of TODIM and fuzzy TOPSIS and FTF in our proposed method.
In previous methods used in TODIM or other MCDM methods, data have been collected based on alternative and criteria comparisons, but we use interval-valued triangular fuzzy numbers which consider the criteria with respect to alternatives in a deterministic time and these data collection will be continued in other times which are important for decision makers, or when we have some changes in criteria or alternatives.Table 2 shows transformation of linguistic criteria to fuzzy-interval triangular numbers.
Let  and  be defined as TFN:  = ( 1 ,  2 ,  3 ),  = ( 1 ,  2 ,  3 ).The distance between  and  is The TFN is based on a three-value judgment: the minimum possible value  1 , the mean value  2 , and the maximum possible value  3 .But the interval-TFN values that have been used in this paper have 5 parameters; on the other hand these are two TFN numbers which have the same middle point (Figure 2).
is expressed as an interval-value TFN where Given   = [( 1 ,   1 );  2 ; (  3 ,  3 )], the normalized performance rating is  The distance between the reference value and each comparison value can be calculated by using definition (2) as follows: These calculations are used to determine the distance between the reference value and the comparison value in the interval; after calculation we have a new interval TFN for every TFN as   = [ (1)   ,  (2)   ].The weight vector of the criteria is calculated according to Zhnag et al. [2]   = ∑  =1 ( (1)   +  (2)   ) ∑  =1 ∑  =1 ( (1)   +  (2)  ) .

Prospect Theory.
The value function used in the prospect theory is described in the form of a power law expressed as where  and  are parameters related to gains and losses, respectively.The parameter  is a risk factor that is considered in model and must be greater than one.Figure 3 shows prospect value function that must be concave and has an  shape form.

TODIM Formulation.
The TODIM method uses pairwise comparisons between the criteria by using technically simple resources to eliminate occasional inconsistencies resulting from these comparisons.TODIM allows value judgments to be performed on a verbal scale using hierarchy of criteria, fuzzy value judgments, and interdependence relationships among the alternatives.The decision matrix consists of alternatives and criteria.
TODIM then calculates the partial dominance matrices and the final dominance matrix.The first calculation that the decision makers must define is a reference criterion (typically the criterion with the greatest importance weight).Therefore, w  indicates the weight of the criterion  by the reference criterion .TODIM is expressed by the following equations.

Advances in Fuzzy Systems
The dominance of an alternative over the other is as follows: where The term   (  ,   ) represents the contribution of criterion  ( = 1, . . ., ) to the function  (  ,   ) while comparing alternative  with alternative .The parameter  represents the attention factor of the losses, whose mitigation depends on the specific problem.A positive (  −   ) represents a gain, whereas a nil or a negative (  −  ) represents a loss.The final matrix of dominance is obtained by summing the partial matrices of dominance for each criterion; see Tseng et al. [3].The global value of the alternative  is determined by normalizing the final matrix of dominance according to the following expression: Ordering the values   provides the rank of each alternative, and better alternatives have higher values of   .The use of numerical values in rating alternatives may be limited in their capacity to address uncertainties.Therefore, an extension of TODIM is proposed to solve problems with uncertain data resulting in fuzzy TODIM.In practical applications, the triangular shape of the membership function is often used to represent fuzzy numbers.Fuzzy models that used TFNs proved highly effective for solving decision making problems when information is imprecise.Hence, this study provides some basic definitions of fuzzy set theory; see Tseng,[24].
In this study we first use Fuzzy TOPSIS for weighting the criteria considering time factor and then we use this data to combine with data from TODIM which calculate from comparisons between alternatives regarding to criteria, and finally we will able to determine which alternative will be more reliable and more effective in any duration of time.

Illustrative Method and Example
3.1.Proposed Approach.We use the abovementioned method to find out the best supplier in green supply chain management considering time variations as follows (Figure 4).
Step 1.A group of decision makers identified the criteria in GSCM which are important and also will be changed during time horizons.
Step 2. Collect the opinion of decision makers with linguistic variables (Table 1) and define FTF for each important criterion.
Step 3. Use fuzzy TOPSIS to evaluate the criteria during time and initial weight of each criterion.
Step 4. Collect the opinion of decision makers on alternatives, respectively, with linguistic variables (Table 2).
Step 5. Use TODIM for evaluating the final weight of each criterion against alternatives and the relationship between them.
Step 6. Combine the results of TODIM and fuzzy TOPSIS to find out which supplier will be more effective from our company's imagination now and in the future regarding the condition changes.

An Illustrative Example
Step 1.In this section we study on green supplier selection problem based on time factor in a Tier company in Iran.
In this company regarding the expert researches we have 6 important criteria.The data have been collected from three  expert decision makers who have more than 10 years of experience in this area and also have the ability to predict the market and its requirements in future.Table 3 represents the criteria of GSCM which is considered in this study for our three suppliers.Table 4 represents data and then uses them to produce three functions for each criterion, but because of lack of data we consider just three points of time: 10, 20, and 30 months after data collection.Step 2. The data in Table 4 were prepared based on linguistic fuzzy variable from Table 1, and data are collected from three experts to find out their opinion on criteria with respect to time.After collecting data with linguistic variable, we can change them to fuzzy triangular number as represented in Table 1 and then FTF must be represented.There are two methods to calculate the FTF.
First.In this case we use average of slopes, when the deviation between slopes is low and is not very important; in our example for criterion 1 we have the following.Criterion 1 (0.33,1. The FTF from (2) for criterion 1 is To get precise decisions, it is recommended to use the separate functions for each period of time.For example for criterion 1, using the function between 10 and 20 months is more exact than considering the FTF for any time.Calculations for finding the best green supplier in our study for the 15th and 25th months are presented in Table 5; for example for Criterion   By using FTF, we will have the variance for calculations that is useful in some other analysis which could not be achieved by ordinary triangular fuzzy numbers.
Second.In this method we did not use mean of the slopes because of more deviation in numbers and in these cases we use separate functions in calculations that will be more reasonable.For example for criterion 4 we have We use the first method for C1 and second method for other criteria.
Step 3. Table 5 shows that calculated numbers come from FTF for criteria; these data will be the input for fuzzy TOPSIS in (12) and (13).The results of calculations are represented in Table 6.This data will be the initial weight of TODIM approach.
By using the fuzzy TOPSIS method we have Table 6.As shown in Table 6, we will have different ranking for criteria with respect to expert opinion changes during the time that will be measurable by FTFs.This evaluation helps us to know about the importance of each criterion in these times and gives us the perspective for more effective actions in our company in future and Table 7 presents the ranking of criteria with respect to time.
Step 4. Table 8 shows the experts' opinion on alternatives in each criteria; these data come from the linguistic variables and then change to interval valued triangular fuzzy numbers by Table 2. Step 5. Use TODIM methods; see Tseng et al. [3].
By applying TODIM approach, first  − ,  * have been calculated from (17) which is shown in Table 9 after that the weight of criteria has been calculated from (18) which is shown in Table 10.
After calculating weight of each criterion, the data must be normalized for calculation dominance weight of criteria.The normalized data is shown in Table 11.
Equation (23) calculates the overall value of alternative by normalizing the corresponding dominance measurements.The rank of each alternative derives from ordering the alternatives values.The global measures computed the complete rank ordering of all alternatives.In addition, a sensitivity analysis should then be applied to verify the stability of the results based on the decision makers' preferences.The results are presented in Table 12.
Considering results in Table 12, we will find out which criterion is the most important one, we can also find out ranking of alternatives regarding each criterion.
Step 6.Now the final conclusion will be conducted from Tables 12 and 6 and we can find out in next 15 months and 25 months which criteria will be more effective and finally which alternative will be the first.By combination of Fuzzy TOPSIS which represents the weight of criteria of green supply chain with respect to time and TODIM method which shows the weight of alternatives according to criteria of supply chain we will have weights of criteria and alternatives in Tables 13 and 14.

Conclusions
This study represents a new way to select supplier during time periods by using a hybrid MCDM.Using fuzzy TOPSIS to have more accurate weighting method for TODIM and then  combining it with fuzzy time function (FTF) helped us to rank the criteria and alternatives in several time periods.The important aim of this study was to use fuzzy time function with a new approach to consider time for each triangular fuzzy number that helped us to find out which supplier in the future according to criteria is suitable for our green supply chain and when we should change our supplier.According to the results shown in Table 14, supplier A3 has a decreasing rate from month 15 to 25, but its decreasing rate is less than supplier 2's.Regarding these results we can assign our future improvement programs and orders to these suppliers considering criteria changes in green supplier selections and also have some repairing programs due to our vision and managerial decisions in future for these companies.Also we can assign our orders to these companies and we proposed an order plan for these suppliers regarding their portions.

Figure 4 :
Figure 4: Flowchart of the method.

Table 1 :
Linguistic variable representing triangular fuzzy numbers.

Table 3 :
Criteria for green supply chain.Annual growth in green products ( 1 ) Cost of revenue: extent that it remains flat to decreases each year ( 2 ) Industry leadership: green market share ( 3 ) Customer retention/percentage of growth with existing customers ( 4 ) Customer acquisition: the number of new green customers/total revenue to new green customers ( 5 ) Life cycle assessment ( 6 )

Table 4 :
Experts data collection.

Table 5 :
Green supply chain criteria considering time.

Table 6 :
Criteria weights using fuzzy TOPSIS considering time periods.

Table 7 :
Criteria ranking for green supply chain considering time.

Table 9 :
Matrix of alternative scores with respect to criteria.

Table 12 :
Final weight of criteria by using TODIM.
Negative numbers explain losses from goal.

Table 13 :
Criteria weight's changes during time.

Table 14 :
Alternative weights during time.