Optimization of Production-Distribution Problem in Supply Chain Management under Stochastic and Fuzzy Uncertainties

Production-Distribution Problem (PDP) in Supply Chain Management (SCM) is an important tactical decision. One of the challenges in this decision is the size and complexity of supply chain system (SCS). On the other side, a tactical operation is a midterm plan for 6–12 months; therefore, it includes different types of uncertainties, which is the second challenge. In the literature, the uncertain parameters were modeled as stochastic or fuzzy. However, there are a few studies in the literature that handle stochastic and fuzzy uncertainties simultaneously in PDP. In this paper, themodeling and solution approaches of PDPwhich contain stochastic and fuzzy uncertainties simultaneously are investigated for a SCS that includes multiple suppliers, multiple products, multiple plants, multiple warehouses, multiple retailers, multiple transport paths, andmultiple time periods, which, to the best of the author’s knowledge, is not handled in the literature.The PDP contains deterministic, fuzzy, fuzzy random, and random fuzzy parameters. To the best of the author’s knowledge, there is no study in the literature which considers all of them simultaneously in PDP. An analytic solution approach has been developed by using possibilistic programming and chance-constrained programming approaches. The proposed modeling and solution approaches are implemented in a numerical example. The solution has shown that the proposed approaches successfully handled uncertainties and produce robust solutions for PDP.


Introduction
The global competition enforces the firms to manage their facilities more effectively and to make right decisions in the market.Supply Chain Management (SCM), which is defined as the integration of key business processes from end user through original suppliers which provides products, services, and information that add value for customers and other stakeholders by the Global Supply Chain Forum (GSCF) [1], is a useful management approach to survive in the global market.SCM includes several processes such as supplier relation management, product development and commercialization, procurement, order fulfillment, manufacturing flow management, demand management, customer relationship management, returns management, and information management [1].Production-Distribution Problem (PDP) in SCM is an important planning operation that affects several processes such as procurement, order fulfillment, and manufacturing flow management.PDP starts to plan by determining raw materials provided by suppliers and makes decisions about the production planning and the distribution of final products to customers.The researchers and practitioners have been interested in PDP over the past years.Fahimnia et al. [2] indicated that there might be two main reasons increasing the number of studies on PDP: (1) affecting the profitability and (2) responding to the market changes quickly.The studies on the PDP can be classified into the different clusters according to the different criteria such as complexity of the supply chain system (SCS), decision levels, solution approaches, and structure of parameters.
PDP can be handled at different decision levels such as operational, tactical, and strategical levels.The strategic decisions are long-term plans that have vital effects on surviving in the market.The papers in this cluster focus on supply chain network design.They also consider opening plants, warehouses, and so forth [3][4][5].In tactical perspective, PDP can be used to determine the production and transportation quantities for aggregate production planning and distribution planning.Besides, it is useful for capacity and resources planning decisions [6][7][8][9].The PDP in operational level seeks to optimize the SCS by adding operational decisions to the aggregate models such as scheduling problem and routing problem [10][11][12][13].
There are several differences between operational, tactical, and strategical decisions such as time period, detail of information, responsibility, and the cost of a wrong decision.One of them is uncertainty which depends on the length of time period.The precision and exactness of information about problem parameters decrease when the time period of decision increases.Therefore, uncertainty is a challenge in PDP.PDP can be classified into four groups according to the structure of parameters.The first group is deterministic parameters: these models do not include any uncertainty in their parameters.All of the parameters are exact and are known at the beginning of solution process [12,14].The second group is stochastic parameters: the parameters include stochastic uncertainty.The probability theory models these parameters [3,[15][16][17].The third group is fuzzy parameters: the fuzzy set theory is an effective modeling approach when the information on parameters is imprecise or inexact.It enables reflecting the decision maker's judgements into the problem [18][19][20][21].The fourth group is fuzzy and stochastic parameters: in some situations, both fuzzy and stochastic uncertainties can occur in parameters simultaneously such as fuzzy random or random fuzzy parameters [22][23][24].
The PDP requires using various techniques for solving this problem because of the properties of the PDP which are discussed above.Fahimnia et al. [2] classified these techniques into four clusters: analytic techniques, heuristic techniques, simulation, and genetic algorithms.For analytic techniques, the studies in this cluster use mathematical programming to solve PDP, that is, linear programming, nonlinear programming, mixed integer programming, and Lagrangian relaxation [44][45][46].For heuristic techniques, since analytic techniques have a limitation on solving large-scale PDP, the researchers developed heuristic techniques that obtain feasible solution close to an optimal solution [16,35,47].For simulation modeling, simulation is a very useful tool to analyze the system's behavior and performance criteria when the considered system is very complex to solve analytically [30,48].For genetic algorithms (GA), they are effective algorithms that use direct and stochastic search methods to solve large-scale problems [49,50].
In this paper, the PDP has been handled from a tactical perspective for a SCS.The SCS includes multiple suppliers, multiple products, multiple plants, multiple warehouses, multiple retailers, multiple transport paths, and multiple time periods, which, to the best of the author's knowledge, is not handled in the literature.A, 0-1 mixed-integer programming model has been developed for the PDP which includes deterministic, fuzzy, fuzzy random, and random fuzzy parameters.To the best of the author's knowledge, there is no study in the literature which considers deterministic, fuzzy, fuzzy random, and random fuzzy parameters simultaneously in PDP.An analytic solution approach has been developed for 0-1 mixed-integer programming model by using possibilistic programming and chance-constrained programming approaches.
The paper is organized as follows: the modeling uncertainty is given in Section 2. Section 3 represents mathematical model and uncertain parameters for PDP.The proposed solution approach is given in Section 4. The implementation of the proposed solution approach for a real-life industry case is presented in Section 5.The paper is finalized with concluding remarks in Section 6.

Modeling Uncertainty
Let us give the definitions of some uncertainty types such as random, fuzzy, random fuzzy, and fuzzy random variables.Definition 1.If  is an experiment having sample space ℘ and  is a function that assigns a real number () to every outcome  ∈ ℘, then () is called a random variable [51].Definition 2. Let Ω be a set of all outcomes of a random experiment.A (nonempty) collection  of subsets (called events) of Ω is assumed to have the following properties: (a) Ω ∈ ; (b) if A ∈ , then A  ∈ ; and (c) if A  ∈  is a countable sequence of events, then ⋃    ∈ .Such a collection  is called a -algebra.For each random event , there is a nonnegative number Pr{}, called its probability, such that (i) Pr{⌀} = 0 and Pr{Ω} = 1 and (ii) Pr{⋃    } = ∑  Pr{  } for every countable sequence of mutually disjoint events   .The triplet (Ω, , Pr) is called a probability space and the function Pr is referred to as a probability measure.A random variable on the probability space (Ω, , Pr) is a function  from Ω to the real line R for any Borel set  of R [52].
Normal distribution is very important in both theory and application of statistics.The notation  ∼ (,  2 ) is often used to indicate that the random variable  is normally distributed with mean  and variance  2 .
After random variable definition, now we can get fuzzy variable definition and properties.Fuzzy set theory was proposed by L. Zadeh and applications of his theory can be found, for example, in artificial intelligent, computer science, control engineering, operation research, and decision theory [53].Definition 3. Let  denote a universal set.Then a fuzzy subset Ã of  is defined by its membership function  Ã :  → [0, 1] which assigns to each element  ∈  a real number  Ã() in the interval [0, 1], where the value of  Ã() at  represents the grade of membership of  in Ã.A fuzzy variable is defined as a function from the possibility space (Θ, (Θ), Pos) to the real line R [52].
Triangular fuzzy variable is the most known and used fuzzy variable which is denoted by the triplet (  , ,   ) and has the shape of a triangle.
The concept of the random fuzzy variable was initialized by Liu and defined as a fuzzy variable taking "random values."Definition 4. A random fuzzy variable is defined as a function from the possibility space (Θ, (Θ), Pos) to the set of random variables [54].
Assume that  1 ,  2 , . . .,   are random variables and that A random fuzzy variable is defined as a function that assigns a random value to each fuzzy subset.On the other hand, a fuzzy random variable is defined as a function that assigns a fuzzy subset to each possible output of a random experiment.

Mathematical Model and Uncertain
Parameters for PDP

Mathematical Model.
The SCS includes multiple suppliers, multiple products, multiple plants, multiple warehouses, multiple retailers, multiple transport paths, and multiple time periods.In the PDP, multiple raw materials are supplied from multiple suppliers and transported to the multiple plants by using multiple transport paths.In plants, multiple products are manufactured by using regular time and overtime and the final products are transported to the multiple warehouses by using multiple transport paths.Multiple warehouses deliver multiple products to the multiple retailers by using multiple transport paths.The customers pick up their products from multiple retailers.
Several assumptions have been made to construct a 0-1 mixed-integer programming model which are given as follows: (i) The quantities of raw materials in suppliers are restricted.
(ii) The number and capacity of transport paths between all the components in SCS are restricted.
(iii) The starting and ending inventories of product and raw materials in plants, warehouses, and retailers are zero.
(iv) The plants have ability to produce several products.
(v) The plants have ability to store raw materials and products.
(vi) The storage capacities of plants are restricted.
(vii) The plants have regular-time and overtime production.
(viii) The warehouses have ability to store products.
(ix) The storage capacities of warehouses are restricted.
(x) The retailers have ability to store products.
(xi) The storage capacities of retailers are restricted.
(xii) The unsatisfied demands are lost.
In this paper, the PDP in SCM has been considered in tactical level.Since it is a mid-term plan, it includes a lot of uncertainties.The uncertainties make the problem more complex compared to the deterministic ones because there are challenges in modeling parameters and obtaining robust solutions.In the literature, the researchers have tried to overcome these challenges by using fuzzy set theory or probability theory.The fuzzy set theory provides a highly effective means of handling imprecise data.It enables incorporating the decision-maker's expertise and judgements into the problem.However, it is not a powerful theory like probability theory for modeling and solution.On the other side, the probability theory is an effective tool for modeling uncertainties in the stochastic process.It acts with the past data analysis for the forecasting of future events and does not include decisionmakers into the decision-making process.

Mathematical Problems in Engineering
However, the decision-makers have an impact on the future events by the way of their decisions.In PDP, the unit cost of raw materials may change based on the purchased quantity.The unit production cost directly depends on lot size.Producing in overtime or producing and holding in regular time at previous periods are based on planning manager's decision.The unit transportation cost is related to path type and transported quantity.The unit price of a product may be changed by making discount, giving an advertisement.The capacity of raw materials supplied from the market is based on the contracts made by the decisionmakers.The decision-maker can change the workforce level by hiring and firing; therefore, the production capacity can be changed.The parameters in the model related to the above discussion can be modeled by using triangular fuzzy numbers.The triangular fuzzy numbers are well known and are commonly used in many applications because the decision-maker has opinions about pessimistic, optimistic, and most possible values by using his/her expertise and expectation.
The transportation capacities of all echelons in the SCS are related to the number of the transporters in the portfolio of the decision-makers, transportation quantities, and vehicle routing decisions that make the transportation capacity uncertain.Therefore, transportation capacities can be modeled as triangular fuzzy numbers by using the decisionmaker's expertise and judgements.However, the available transportation capacity can occur in different situations based on the suitability of the transporter in the market which are defined as discrete events.These discrete events are determined as high, medium, and low capacities.It is possible to increase the number of situations; however, it will cause confusion in the categorization process.By analyzing the past data, probability levels can be determined for occurrences of each of the situations.Therefore, the transportation capacities can be modeled as fuzzy random parameters.
The demand of product can be modeled as a probability distribution by analyzing the past data.Since the PDP is a mid-term plan, a sum of identically distributed independent demand variables has a normal distribution according to the central limit theorem.However, the decision-maker can affect the demand quantity by making discount, advertisement, or other strategies.These marketing strategies are based on management decisions.Therefore, the demand quantity is modeled as a random fuzzy variable.
The mathematical model is given in Notations.
Mathematical Problems in Engineering TRQ  , SRP  , RPQ  , OPQ  , SLP  , TPQP  ≥ 0, The objective function, given in (2), maximizes the total profit.Total profit is obtained by total revenue, which is gained from total sales, plus total cost.Total cost includes raw material purchasing cost, fixed costs of using path  for transportations between all components of SCS, unit transportation costs between all components of SCS, unit production costs in regular time and overtime, holding costs in plants, warehouses, and retailers, and backorder cost.Equation ( 3) is a capacity constraint for the supplier that ensures that the total transported quantity from supplier  for material  at each period will be less than or equal to total capacity.Equations ( 4) and ( 5) are constructed to select the transportation path from supplier to plant and not to exceed its capacity where  is a big number.Equation ( 6) is an inventory balance constraint for the raw material in a plant.Equations ( 7), (8), and ( 9) are capacity constraints for regular production, overtime production, and inventory level in plants, respectively.Equation ( 10) is an inventory balance constraint for the product in a plant.Equations ( 11) and ( 12) are designed to select the transportation path from a plant to the warehouse and not to exceed its capacity.Equation ( 13) is an inventory balance constraint in a warehouse.Equation ( 14) is an inventory capacity constraint for a warehouse.Equations ( 15) and ( 16) are designed to select the transportation path from a warehouse to the retailer and not to exceed its capacity.Equation ( 17) is a balance constraint for inventory and backorder level.Equation ( 18) is an inventory capacity constraint for a retailer.Equation ( 19) ensures meeting customer demand.Equation (20) gives the definitions of the decisions variables.

Uncertain Parameters in PDP.
where Pr  , Pr  , and Pr  represent the probabilities of transportation capacity situations such as high, medium, and low capacities.On the other hand, triangular fuzzy numbers (  , ,   ), (  , ,   ), and (  , ,   ) represent the amounts of each of the transportation capacities for probability levels.

Random Fuzzy Parameters.
There is one random fuzzy parameter, CDP  , which represents the customer demand.The customer demand includes two main parameters: the probability and quantity.Therefore, the demand can be calculated as the sum of multiplication of probability value and quantity which can be referred to as expected value of discrete random variable.The probability and quantity of demand are random fuzzy variables.
The probability of demand is modeled as follows: there are three states about the demand; it may be high, medium, or low.Assume that these three probabilities are represented as Pr() high , Pr() medium , and Pr() low .The probability of demand state is affected by three indicators which are (1) political developments, (2) competitors' strategies, and (3) sectoral expectation.For example, if the competitors perform a strong strategy in the market, the demand quantity will be affected by this situation; most likely it will decrease.These indicators are related to the expertise and expectations.Therefore, it is a very difficult task to model the demand states in deterministic or stochastic case.However, random fuzzy variables enable modeling these situations easier than the remaining ones.It is possible to reflect the decisionmaker's judgements and expectations into the demand state by using random fuzzy variables.The modeling of demand states by using random fuzzy variables can be explained with an example for situations (1) and (2) in Table 1.
The first situation assumes that political development will be good, competitors' strategy will be medium, and sectoral expectation will be good.According to the decision-maker's judgements, expertise, and expectation, the possibility of occurrence of situation (1) is one.It means that situation (1) is an event that can absolutely occur.The analysis of historical data shows that when situation (1) occurs, the probabilities of demand which may be high, medium, and low are 0.8, 0.15, and 0.05, respectively.Situation (2) can be interpreted like situation (1).Consequently, the only way to model Pr() high , Pr() medium , and Pr() low is using random fuzzy variables.According to Definition 4, Pr() high can be modeled as follows:

Solution Approach
The idea of uncertain programming is to convert the uncertain nature of a model into an equivalent deterministic one [56].Therefore, the uncertain parameters in PDP will be transformed into some equivalent deterministic ones by using properties of fuzzy, fuzzy random, and random fuzzy variables.

Transforming Uncertain Parameters into Deterministic
Equivalents.The uncertain parameters have occurred in both objective function and constraints.Therefore, transforming operations of uncertain parameters are considered based on the location of uncertain parameters in the mathematical model.

Uncertain Parameters in Constraints.
Let transformation operation of fuzzy parameter start and that operation is called "defuzzification" in the literature [57].
Definition 8.The -cut of a fuzzy set Ã is a crisp subset of  and is denoted by Ã = { |  Ã() ≥  and  ∈ }.
The -cut of the triangular fuzzy variable Definition 9.The multiplication of a fuzzy variable Ã by a real number  > 0 can be defined [58]: A real number  can be defined as a triangular fuzzy number by the triplet (  , ,   ), where   = ,  = , and   = .Definition 10.Assume that X and Ỹ are two fuzzy numbers.The result  of the addition of the fuzzy numbers X and Ỹ can be defined by the -cut sets [59].That is, All of the fuzzy parameters in the right-hand sides of (3), (7), and ( 8) can be transformed into deterministic close interval by using -cut approach.Now let consider fuzzy random parameters.Definition 11.Let  be a discrete random variable taking values  1 ,  2 , . . .with probabilities Pr 1 , Pr 2 , . .., respectively.Then the expected value of this random variable is the infinite sum All of the fuzzy random parameters in the right-hand sides of ( 4), (11), and ( 15) are transformed into deterministic close interval according to Corollary 12. Now let consider random fuzzy parameters.As defined in Section 3.2, the customer demand has three discrete events; it may be high, medium, or low with probability values Pr() high , Pr() medium , and Pr() low , respectively.On the other side, the demand quantity for each event follows normal distribution with different fuzzy mean parameters and different variances which are ( X( high/medium/low ), ( high/medium/low )).According to Definition 11, total customer demand can be written as follows: Total customer demand = Pr () high *  ( X ( high ) ,  ( high )) + Pr () medium *  ( X ( medium ) ,  ( medium )) + Pr () low *  ( X ( low ) ,  ( low )) , (25) where Pr() high/medium/low and ( X( high/medium/low ), ( high/medium/low )) are random fuzzy parameters; therefore, total customer demand is a random fuzzy parameter.
In order to transform total customer demand into its deterministic equivalent, it is required to transform the probabilities and quantities into deterministic cases.
The following definition and corollary have been made for transforming the probabilities.Definition 13.Let  be a normalized discrete fuzzy variable whose possibility distribution function is defined by The expected value of  is as follows: where the weights are given by ( 0 = 0,  +1 = 0) for  = 1, 2, . . .,  and satisfy the following constraints:   ≥ 0 and ∑  =1   = max 1≤≤   = 1, since  is a normalized fuzzy variable [60].

Corollary 14. If the demand state is a discrete random fuzzy variable, the probabilities of demand states are fuzzy variables according to Definition 5 and then by using expected value of the fuzzy variable (Definition 13), crisp expected probability values can be calculated for high, medium, and low demand states.
According to Corollary 14, the expected probability values for high, medium, and low demand states, which are represented by  high ,  medium , and  low to prevent confusions in next formulations, are written as follows:  (25), can be transformed into a random variable by using Definitions ( 8), ( 9), (10), and (11) According to Corollary 16, the right-hand side of ( 19) is transformed into a random parameter; however, it is still uncertain.The chance-constrained programming can be used to obtain its deterministic equivalent.
The structure of a chance-constraint is as follows [56]: It means that the constraint is realized with a minimum probability of 1 − .If  is normally distributed parameter,  ∼ (  ,  2  ), the constraint is converted as follows: where (−  )/√ 2  represents a standard normal variate with a mean of zero and a variance of one.Then, the stochastic chance-constraint is transformed into the following inequality: where Φ( 1− ) = 1 −  and Φ( ) represents the standard normal cumulative distribution function.This yields the following linear deterministic constraint: 4.1.2.Uncertain Parameters in Objective Function.The uncertainties in the constraints are converted into their deterministic equivalents.However, objective function still includes fuzzy parameters.Therefore, Lai and Hwang's [59] approach has been used to obtain deterministic equivalent of the objective function.
Lai and Hwang [59] had handled a mathematical model as given in the following equation: max c, where Ã, b, and c are triangular fuzzy numbers.
The fuzzy objective function is fully defined by three corner points (  , 0), (  , 1), and (  , 0) geometrically.Lai and Hwang [59] suggested that maximizing the fuzzy objective can be obtained by pushing these three critical points in the direction of the right-hand side.The vertical coordinates of the critical points are fixed at 1 or 0. The only considerations then are the three horizontal coordinates.Therefore, the objective function is translated to the form given in the following equation: Instead of maximizing these three objectives simultaneously, Lai and Hwang [59] proposed maximizing   , minimizing [   −   ], and maximizing [   −   ].The proposed approach involves maximizing the most possible value of the profit, minimizing the risk of obtaining lower profit, and maximizing the possibility of obtaining higher profit.The last two objectives actually are relative measures from   .This leads us to the auxiliary multiobjective linear programming model given in the following equation: Lai and Hwang suggested using Zimmermann's [61] fuzzy programming method to convert the auxiliary multiobjective linear programming model into an equivalent single-goal LP problem.First, the positive ideal solutions (PIS) and negative ideal solutions (NIS) of the objective functions can be specified as follows [59]: The linear membership function of each objective function is defined as follows: Lai and Hwang used fuzzy ranking concepts for the constraints and combined them with their strategy for imprecise objective function.The constraints can be modeled by using -cut approach as follows: If only the right-hand sides include fuzzy parameters, Lai and Hwang propose the weighted average method to obtain crisp right-hand side values.Assume that only the right-hand side of the constraint in (35) ( b) is fuzzy.For a given minimum acceptable possibility, , the crisp equality constraints can be constructed as follows: where The objective function which is fully fuzzy has been handled by using Zimmermann's [60] fuzzy programming method.Therefore, there is no different technique in the proposed approach to convert the objective function.However, different techniques are used in constraints.
The goals of determining positive and negative ideal solutions are to calculate the minimum and maximum values of objective functions.Therefore, the positive and negative ideal solutions are determined according to the pessimistic and optimistic scenarios in uncertain models to obtain robust solutions.However, Lai and Hwang proposed a weighted average method in constraints that only includes fuzziness on right-hand side for obtaining positive and negative ideal solutions.In weighted average method,  1 =  3 = 1/6 and  2 = 4/6.This method produces a crisp value that is very close to the most possible value.Therefore, the weighted average method prevents obtaining lower and higher ideal solutions.Naturally, the weighted average method may produce unfeasible solution.
The algorithm of the solution methodology for practical PDP decisions is as follows.
Step 1. Formulate the PDP model.
Step 2. Model the fuzzy parameters as triangular fuzzy numbers, model the discrete fuzzy random parameters (the capacities of transportation paths) as triangular fuzzy numbers with probability values, and model the random fuzzy parameter (demand) as normal distributions with triangular fuzzy mean parameters.
Step 3. Develop three new crisp objective functions of the auxiliary MOMILP problem from the fully fuzzy objective function which are equivalent simultaneously maximizing the most possible total profit, minimizing the risk of obtaining lower profit, and maximizing the of possibility of obtaining higher profit.
Step 4. Determine an  value and transform the fuzzy parameters in the right-hand side of the constraints into deterministic close intervals by using -cut approach.Produce two separate constraints from these constraints that one of them uses lower bound of close interval and the other one uses upper bound of close interval.
Step 5. Transform the discrete fuzzy random parameters in the right-hand side of the constraints into deterministic close intervals according to Corollary 12. Produce two separate constraints from these constraints that one of them uses lower bound of close interval and the other one uses upper bound of close interval.
Step 6. Transform the random fuzzy parameters in the right-hand side of the constraints into normally distributed random parameters with deterministic close interval mean parameters according to Corollaries 14 and 16.
Step 7. Determine an acceptable probability value () and model the constraint obtained in Step 6 as deterministic linear constraint by using chance-constraint approach.Produce two separate constraints from this constraint that one of them uses lower bound of close interval and the other one uses upper bound of close interval.
Step 8. Solve  1 PIS ,  1 NIS ,  2 PIS ,  2 NIS ,  3 PIS , and  3 NIS .Obtain maximum and minimum values for  1 ,  2 , and  3 .Specify the linear membership functions for each of them, and then convert the auxiliary MOMILP problem into a singleobjective 0-1 mixed-integer programming model.
Step 9. Solve the single-objective 0-1 mixed-integer programming model and obtain the solution.
Step 10.If the DM is not satisfied with the initial solution, return to Step 4 and modify  and  values and repeat the remaining steps until a satisfactory solution is found.

Implementation of the Proposed Solution Approach
The proposed modeling and solution approaches have been implemented for furniture manufacturer in Turkey.The firm has two plants ( = 1, 2) in Ankara and Bursa.It produces three different types of products ( = 1, 2, 3) by using six different types of raw materials ( = 1, 2, . . ., 6) which are supplied from three suppliers ( = 1, 2, 3).The firm uses three warehouses ( = 1, 2, 3) and five retailers ( = 1, 2, . . ., 5) to deliver the products to the costumers by using three different types of transportation paths ( = 1, 2, 3) which are trucks, trains, and planes.The demands of the products significantly increase in summer; therefore, time period is three months, that is, June, July, and August ( = 1, 2, 3), for planning horizon.The costumers were not categorized ( = 1) because of tactical planning decision.The proposed modeling and solution approaches were performed at February 2017 for summer season in 2017.The assumptions given in Section 3.1 are held for this real-life application.The implementation of the algorithm for real-life application is given as follows.
Step 2. Unit prices of raw materials ( RUP  ) which are given in Table 8, capacities of raw materials ( RC  ) which are given in Table 9, unit transportation costs of raw material from suppliers to plants ( ṼTCS  ) which are given in  22, unit transportation costs of products from plants to warehouses ( ṼTCP  ) which are given in Table 26, unit transportation costs of products from warehouses to retailers ( ṼTCW  ) which are given in Table 30, unit holding costs of products in warehouses ( WHC  ) which are given in Table 31, unit holding costs of products in retailers ( HCR  ) which are given in Table 33, backorder costs of products ( BCR  ) which are given in Table 35, and price of products ( POP  ) which is given in Table 37 have been modeled as triangular fuzzy numbers by making a discussion and analysis with a group of staff which contains procurement, production planning, marketing, and warehouse and retailers managers.
It has been observed that the capacities of transportation paths between echelons of SCS can occur in three different situations, high, medium, and low capacities, when the historical transportation data was analyzed.The probabilities of obtaining high, medium, and low capacities have been determined as 0.5, 0.35, and 0.15, respectively.However, for each situation, the quantities of capacities can change based on the availability of transporters.Therefore, total capacities of transportation from suppliers to plants (TCSP  ) which are given in Table 10, total capacities of transportation from plants to warehouses (TCPW  ) which are given in Table 24, and total capacities of transportation from warehouses to retailers (TCWR  ) which are given in Table 28 have been modeled as triangular fuzzy numbers by making a discussion and analysis with a group of staff.
It has been observed that the demands of products can occur in three different states, high, medium, and low demands, by analyzing the historical demand data with bubble charts and the quantity of demand follows a probability distribution for each state.Three indicators have been identified which affect the demand state, that is, political development, competitors' strategy, and sectoral expectation.Six alternative situations have been generated by using these three indicators according to the expertise of management and a possibility value of occurrence has been assigned to each situation by the managers.The probability values of demand that can be high, medium, or low at each situation have been calculated according to the frequency analysis.Therefore, the probability of demand has been modeled as random fuzzy number which is given in Table 1.On the other side, the quantities of the demand for high, medium, and low states follow a normal distribution with a mean and variance parameter according to the results of Anderson-Darling test.However, the managers mentioned that they can affect the demand by using advertisements and discounts.Therefore the mean parameters of the normal distributions have been modeled as triangular fuzzy numbers which are given in Table 36.Fixed costs of using transport paths from suppliers to plants, required transportation capacities of raw materials, required amounts of raw materials for unit products, unit production times of products in plants at each period, inventory capacities in plants, required capacities to store products, fixed cost of using transport paths from plants to warehouses, required capacities to transport unit products, Mathematical Problems in Engineering fixed cost of using paths from warehouses to retailers, and inventory capacities in warehouses are given in Tables 11,13,14,16,21,23,25,27,29,and 32, respectively.
Step 3. Three new crisp objective functions of the auxiliary MOMILP problem have been developed from the fully fuzzy objective function.
Step 4. The  value has been determined as 0.4 and the constraints that contain fuzzy parameters in the right-hand side have been transformed into deterministic close intervals by using -cut approach.Two separate constraints have been produced from these constraints that one of them uses lower bound of close interval and the other one uses upper bound of close interval.
Step 5.The constraints that contain the discrete fuzzy random parameters in the right-hand side have been transformed into deterministic close intervals according to Corollary 12. Two separate constraints have been produced from these constraints that one of them uses lower bound of close interval and the other one uses upper bound of close interval.
Step 6.The constraints which contain the random fuzzy parameters in the right-hand sides have been transformed into normally distributed random parameters with deterministic close interval mean parameters according to Corollaries 14 and 16.
Step 7.An acceptable probability value () has been determined as 0.95 and the constraint obtained in Step 6 has been modeled as deterministic linear constraint by using chance-constraint approach.Two separate constraints have been produced from this constraint that one of them uses lower bound of close interval and the other one uses upper bound of close interval.
Step 8.    The linear membership functions for each of them have been specified and then the auxiliary MOMILP problem has been converted into a single-objective 0-1 mixed-integer programming model.
Step 9. When the equivalent single-objective 0-1 mixedinteger programming model has been solved, the total profit is obtained as a triangular fuzzy variable with 7313900, 9021600, and 10941200, which is given in Figure 1.The overall degree of DM satisfaction with multiple goal values  is achieved at 0.703.The solutions of the model are represented in Tables 2-6.
Step 10.DM has been satisfied with the initial solution and the algorithm has been terminated.
Summary of the results during the three time periods is given in Tables 2-6.Total sales quantities are represented in Table 2.According to these solutions, capacities of the retailers, given in Table 34, are adequate to meet the demands.Therefore, the managers do not focus on the retailers.
Total quantities of products transported from warehouses to retailers are given in Table 3.It is seen that transportation capacities between warehouses and retailers, given in Table 28, are adequate when compared with the results in Table 3.Therefore, the managers do not seek additional transportation capacities.
In Table 4, it is observed that plant 2 has not sent any products to warehouse 3 based on high transformation costs.Therefore, the managers search decreasing transformation costs or any other transformation alternatives from plant 2 to warehouse 3.
On the other hand, plant 1 does not manufacture product 3 according to Table 5.The high production cost may be the reason.The managers seek the way of decreasing production costs.
Table 6 shows that the SCS has not satisfied the demands exactly; therefore, there is a bottleneck in the SCS system.When the SCS system is analyzed in detail, it is identified that the reason of the bottleneck in the system is the limited transportation capacity of products from plants to warehouses.Therefore, there is abundance in the capacities of the remaining components of SCS.
It is obvious that the decision-makers can decrease the backorders costs and consequently increase the total profit by increasing the transportation capacity from plants to warehouses supplied by new transporters.However, they predict that the cooperation with new transporters may cause harmonization problems and extra costs.Therefore, the decision-makers plan to find new transporters and cooperate with them for next planning periods after performing a transporter evaluation system.The real-life application has been performed at different  and  parameter values.Total profits in triangular fuzzy number and defuzzified form for different parameter values are presented in Table 7.The weighted average method was used in defuzzification with  1 =  3 = 1/6 and  2 = 4/6.In the first column of Table 7, the effects of changes in  parameter are examined when  is fixed.The change in  directly affects the demand quantity.Total unsatisfied demand quantity increases when  increases; therefore, total profit decreases for  = 0.98 and  = 0.99.On the other side, it is expected to obtain lower profit when  decreases because of lower demand.However, it is observed that total profit is increased for  = 0.90 based on decreasing unsatisfied demand quantity.
The change in  parameter directly/only affects the SCS capacity (raw material, transportation, and production).tactical level.A modeling approach and a solution approach have been proposed for PDP.
In modeling approach, a 0-1 mixed-integer mathematical model has been developed for PDP.The uncertain parameters in the model have been handled as fuzzy, random, fuzzy random, and random fuzzy variables.A solution approach is developed for 0-1 mixed-integer mathematical model by integrating possibilistic programming and chance-constrained programming techniques.
The proposed modeling and solution approaches have been implemented on a PDP of furniture firm.The PDP has been solved globally optimally by using GAMS optimization package.The solutions have been represented to the firm managers and managers have been satisfied with these solutions.The main problem in SCS is the lower transportation capacity between plants and warehouses.Therefore, managers decide to seek new transporters for the next planning periods.
The proposed modeling and solution approaches can be compared with the other PDP approaches such as deterministic models, stochastic models, or fully fuzzy models.The solution of the deterministic PDP can be easily infeasible in dynamic structure of the system at tactical level because of ignoring uncertainty.It is an alternative way to perform the deterministic model several times with different parameter values to analyze the changes of the parameters.However, hundreds of solutions are obtained for a SCS which includes a lot of uncertain parameters.Consequently, it causes another uncertainty called abundance of the solutions.Therefore, it is not an effective way to make decision at tactical level.
On the other hand, stochastic models only contain probabilistic uncertainties in SCS and do not consider the other uncertainties that are based on decision-maker's expertise and judgements.Probability theory successfully handles the variation in problem nature by assuming that the conditions of the experiment will not change.However, the conditions such as political developments, competitors' strategy, or manager's decision may change in a dynamic system and that creates a situation which has not been observed in the past.Therefore, the solutions of the stochastic models can be infeasible.
Contrary to stochastic models, fully fuzzy models for PDP do not consider the variations in the problem nature.The fully fuzzy models produce rough and subjective solutions.
The main advantage of the proposed modeling and solution approaches is including the decision-maker into the problem formulation and solving processes with probabilistic uncertainties.Decision-makers directly affect the problem parameters by their decisions.Fuzzy, fuzzy random, and random fuzzy variables enable including the decision-maker's judgements and expertise into the probabilistic model.In this way, the mathematical model produces robust solutions.
There are some limitations of the proposed modeling and solution approaches.One of them is the determination of the  value.The  value represents the minimum acceptable possibility degree, in other words, satisfaction level of decision-maker.It may not be meaningful in decisionmaker's mind like probability level.Therefore, it should be explained clearly to the decision-maker in implementation process.Another one is the limitation of the optimization packages.The optimization packages cannot solve the PDP globally optimally for big-size SCS because of the number of the binary variables; therefore, it is required to develop a metaheuristic algorithm.
In this study, single-objective function, maximization of the total profit, has been considered for the PDP in uncertain environment.In future studies, it can be modeled as biobjective functions by considering maximization of the customer satisfaction level or minimization of the total transportation time.The proposed modeling approach can be used in different optimization problems that have fuzzy and random parameters in their nature.
Figure 1 Definition 7. A fuzzy random variable is a function  from a probability space (Ω, , Pr) to the set of fuzzy variables such that Pos(() ∈ ) is a measurable function of  for any Borel set  of R.
[55]with possibility  1 ,  2 , with possibility  2 ,   , with possibility   (1)is a clearly discrete random fuzzy variable[54].Definition 5. Assume that  is a random fuzzy variable.Then the probability Pr{() ∈ } is a fuzzy variable for any Borel set  of R[54].Definition 6.A random fuzzy variable  is said to be normal if, for each , () is a normally distributed random variable; that is, () ∼ ((), ()), with  and  being fuzzy variables defined on the space Θ such that  > 0. A normally distributed random fuzzy variable is usually denoted as () ∼ (, ), and the fuzziness of random fuzzy variable  is said to be characterized by fuzzy vector (, )[55].

Table 12 ,
unit holding costs of raw materials in plants ( SRC  ) which are given inTable 15, unit production costs of products in regular time ( RPC  ) which are given in Table 17, available regular time capacities of plants ( ÃRC  ) which are given in Table 18, unit production costs of products in overtime ( ÕPC  ) which are given in Table 19, available overtime capacities of plants ( ÃOC  ) which are given in Table 20, unit holding costs of products in plants ( PHC  ) which are given in Table 1 PIS ,  1 NIS ,  2 PIS ,  2 NIS ,  3 PIS , and  3 NIS have been solved globally optimally by using GAMS/CPLEX solver.Maximum and minimum values for  1 ,  2 , and  3 have been obtained as follows:

Table 2 :
Total sale quantities of product  in retailer .

Table 3 :
Total quantities of product  transported from warehouse  to retailer .

Table 4 :
Total quantities of product  transported from plant  to warehouse .

Table 5 :
Total quantities of product  manufactured in plant .

Table 6 :
Total backorder quantities of product  in retailer .

Table 7 :
Total profits at different parameter levels.

Table 8 :
Unit price of raw material  in supplier  at period .

Table 9 :
Capacity of raw material i supplied by s at period .

Table 10 :
Capacity of transport path  from supplier s to plant  at period .

Table 19 :
Unit production cost of product  at overtime in plant  at period .

Table 20 :
Available overtime capacity (time) in plant  at period .

Table 21 :
Inventory capacity in plant .

Table 22 :
Unit holding cost of product  in plant  at period .

Table 23 :
Required capacity to store unit product .

Table 24 :
Transportation capacity of path  from plant  to warehouse  at period .

Table 25 :
Fixed cost of using path  from plant  to warehouse  at period .

Table 26 :
Unit transportation cost of product  from plant  to warehouse  by using path  at period .

Table 27 :
Required capacity to transport unit product .
: Transportation path : Customers : Time period.Parameters RUP  : Unitpriceofrawmaterial in supplier  at period  RC  : Capacityofrawmaterial supplied by  at period  TCSP  : Total capacity of transport path  from supplier  to plant  at period