An Integrated Production-Distribution Planning Problem under Demand and Production Capacity Uncertainties: New Formulation and Case Study

In this study, we propose to solve a biobjective tactical integrated production-distribution planning problem for a multisite, multiperiod, multiproduct, sea-air intermodal supply chain network under uncertainties. Two random parameters are considered simultaneously: product replenishment orders and production capacity, which are modelled via a finite set of scenarios, using a two-stage stochastic approach. A corresponding mathematical model is developed, coded, and solved using the LINGO 18.0 software optimisation tool. ,is model aims to simultaneously minimise the total costs of production in both regular and overtime, inventory, distribution, and backordering activities and maximise the customer satisfaction level over the tactical planning horizon. ,e AUGMECON technique is applied to handle with the multiobjective optimisation. ,e applicability and the performance of the proposedmodel are tested through a real-life case study inspired from amedium-sized Tunisian textile and apparel company. Sensitivity analysis on stochastic parameters and managerial insights for the studied supply chain network are argued based on the empirical findings.


Introduction
e tremendous competition in productivity, quality and time response, introduction of products with short life cycles, and intensive customers' expectations are the main concerns of today's global markets.
ese circumstances have pushed industries to aggressively seeking better ways to establish a closer interaction and collaboration between the different activities of the supply chain [1]. Handling decisions of different functions of the supply chain such as production planning and distribution planning are vital for companies to obtain an optimal strategy that improves efficiency and enhances the global performance. Due to the complex nature and dynamics of the relationships amongst the different actors and the hidden customer preference, firms have faced the threat of many sources of uncertainties such as customer demand, prices, and resource capacities [2]. Underestimating uncertainties and failing to account for their significant fluctuations in supply chain planning could lead to unsatisfied customer demand translating to loss of market share. It ultimately leads to undesirable and unexpected impacts on the ability of reaching supply chain goals. us, a major challenge for the decision-maker is to determine how to continuously manage and monitor these uncertainties.
In this context, we propose to develop a biobjective integrated production-distribution planning model under both demand and production capacity uncertainties. A twostage stochastic programming formulation is proposed to help the decision-makers to simultaneously minimise total costs and maximise the customer satisfaction level in terms of on-time delivery. A case study of a manufacturing company producing and exporting fashion products to the European market is tackled in this paper. e efficiency of the proposed approach is illustrated throughout a comparison with the current planning of the company. e remaining of this paper is structured as follows: A literature review of related works is highlighted in Section 2. Section 3 is devoted to the problem statement considering a two-stage stochastic model for an integrated productiondistribution planning problem under uncertainties. e proposed mathematical model is developed in Section 4. e resolution method of the multiobjective optimisation is described in Section 5. In Section 6, a computational study inspired from the textile and apparel industry is presented to assess the performance of the proposed two-stage stochastic programming method. Section 7 displays the related results as well as sensitivity analysis on the variability of stochastic parameters and the main managerial insights. Finally, conclusions and future works are highlighted in Section 8.

Related Literature
e coordinated production and distribution planning problems have drawn more and more attention for many years. As these problems are mutually related, it is critical to optimise these two functions and plan them simultaneously in an integrated manner instead of optimising in a sequential way with little or no integration [3]. ere are a number of review papers on supply chain management tackling this issue [4,5]. Traditionally, the field of supply chain design has focused on mono-objective optimisation quantifying economic objectives, such as cost minimisation [6,7] and profit maximisation [8,9]. Nevertheless, real-world problems often require the consideration of multiple conflicting criteria to measure the customer service level and supply chain responsiveness. e multiobjective optimisation mainly includes goals related to maximising products' quality level [10], minimising the average tardiness, and minimising delivery [11]. Owing to the difficulties of collecting data and knowledge about production and distribution planning problems, some works examined real-life case studies such as automotive manufacturing [11], apparel industry [12], pharmaceutical industry [13], and plastic industry [14]. Many researches suppose the availability of deterministic data when developing a production-distribution planning model. However, the real-world planning problems imply some level of uncertain and incomplete data in one or more of the input parameters. e literature reveals that optimisation problems with uncertainty have been widely tackled and developed. Many survey articles on such models have appeared in [1,15,16].
ere are diverse factors that generate uncertainties in supply chains. e most commonly and extensively stated sources of uncertainty are connected to the customers' demand, process, and supply that have huge effects on the decisionmaking process. To incorporate uncertainty in an optimisation model, several approaches can be addressed such as stochastic programming methods [17], fuzzy programming [10], and robust optimisation techniques [18]. One of the most widely used formulations for decisionmaking under uncertainty is stochastic programming with recourse. e basic idea of this modelling approach is to formulate the problem in a two-stage setting.
is formulation was introduced for the first time in [19]. Another main issue in scenario-based (two-and multistage) stochastic programs consists of generating scenarios associated with their probabilities. Some studies of the production-distribution modelling efforts in the literature that are more relevant to our paper are summarised in Table 1.
In the context of stochastic optimisation, many authors state that the most obvious and extensively studied source of uncertainty has been demand. For instance, Guillén et al. [25] investigated a two-stage stochastic model to incorporate the uncertainty associated with the demand and aimed for the maximisation of profit, demand satisfaction, and financial risk. In addition to the demand uncertainty, several other uncertain parameters have been considered in SC network designs including prices, cost, capacities, and lead time. Rahimi et al. [27] used two-stage stochastic programming considering demand and cost uncertainties. Mohammadi and Mohd [23] proposed a two-stage stochastic programme under uncertainties of costs, demand, and capacity of facilities. Azaron et al. [26] focused on a multiobjective two-stage stochastic programming approach for the supply chain network design in which demands, supplies, processing time, shortage, transportation, and capacity costs were assumed as uncertain parameters. Mirzapour Al-e-hashem et al. [24] addressed a multisite, multiperiod, multiproduct three-echelon supply chain under uncertainties of cost parameters and demand fluctuations. Nemati and Alavidoost [29] studied a multiobjective model through the integration of sales, production, distribution, and procurement planning considering sea-air intermodal transportation modes. e aims were to minimise the total supply chain cost and to minimise the customer service level in terms of maximising total sales. Pasandideh et al. [28] proposed a biobjective optimisation of the supply chain network under uncertainties of costs, customer demand, unit processing time, available production time, and time of setup.
To the best of our knowledge, we notice that it is the first time that the uncertainty of the production capacity in terms of workforce availability is considered in a two-stage stochastic programming approach to model a biobjective tactical integrated production-distribution planning considering sea-air intermodal transportation. Finally, thanks to a real-world case study of the textile-apparel industry, we contemplate a more realistic situation than those examined in [20,21] by considering two sources of uncertainty and multiobjective scheme in our production-distribution model in order to minimise total costs along with maximising the customer satisfaction level in terms of minimising late deliveries.

Problem Statement
e supply chain network, illustrated in Figure 1, consists of multiple production stages composed of several production plants that receive and transform raw materials to    semifinished and end products and deliver them to the customer zones. e tactical planning horizon is often divided into several periods. Note that each plant is characterized by its own production processing time, costs, and capacities. e transportation time between factories from the last production stage and customers' zones depends on the used transportation modes: marine or air transportation. e shortage of not meeting customers' demand on time is acceptable with penalty cost. Nevertheless, all the orders should be fully met at the end of the tactical planning horizon since lost sales are not permitted. Both stochastic demand and stochastic production capacity are considered in this problem, and we assume that these uncertainties can be modelled by defining several different scenarios with a given probability of occurrence. We assume that the company has two types of customers' orders: firm orders that are supposed to be known with certainty at the beginning of the planning horizon and replenishment orders with short due dates that can be processed from a given period. e latter is a part of demand that cannot be initially predicted. So, it will be considered a stochastic order. Moreover, there are a variety of internal factors, which cannot be easily identified and cause huge variations in workers' productivity. We note that the most important ones are the resource availability, workforce skills, and breakdown machines. In our study, we would focus on the workforce absenteeism factor, which leads to a significant fluctuation of the production capacity.
For the considered tactical planning horizon, the optimisation goals are firstly minimising the overall costs, which mainly include production costs in normal working hours and in overtime, inventory-holding costs, shipping costs between production plants, distribution costs to the final customers, and backorder costs. Secondly, we aim to maximise the customer service level by ensuring the delivery of demand on time.
We develop a two-stage scenario-based stochastic programming planning problem, which is formally defined in terms of the following assumptions: (i) Customer demand and production capacity are considered random parameters. (ii) e initial values of inventory and backlog are null. (iii) Production plants from the last stage deliver products directly to overseas customer zones bypassing the distribution network (airline and waterway modes). (iv) It is assumed that all the products produced and transported have no defects. (v) Backorder and overtime production decisions are allowed to avoid volatility in demand and production capacity. (vi) Demand needs not be fully fulfilled all the time since backorders are allowed. However, the lost sales are not permitted.

Mathematical Formulation
is section describes the two-stage stochastic programming approach for an integrated production-distribution planning problem and formulates the mathematical model. We examine a biobjective stochastic integrated productiondistribution planning problem using sea-air intermodal transportation. e stochastic two-stage programming approach represents the random event through a finite number S � c 1, c 2 , . . . , c s of possible realisations, called scenarios with known probabilities δ(s), and basically, the decisions are made in two different stages. In this approach, first-stage decisions are made "here and now" at the beginning of the planning horizon, prior to the realisation of the random events. After the first-stage decisions are taken and the Mathematical Problems in Engineering random events are realised, the second-stage or recourse decisions always referred to as "wait-and-see" decisions are taken as a corrective action after the realisation of the uncertainty [30]. e presence of uncertainty is correlated with the stochastic nature of the demand and production capacity associated with the second-stage decisions. erefore, the objective function consists of the sum of the first-stage decision costs and the expected second-stage recourse costs.
To formulate the problem mathematically, the sets, indices, both deterministic and random parameters, and decision variables are introduced in Table 2. Let ST, I, P, M, and K be the sets of production stages, products, plants, transportation modes, and customers, respectively. e representation set for the probability and the number of scenarios are expressed as S and δ(s), respectively (Table 2). Based on the abovementioned assumptions and description, the two-stage stochastic programming model for the tactical production and distribution planning problem could be formulated as a multiobjective integer linear programming model. e objectives correspond to minimisation of the overall costs of the supply chain and maximisation of customer services that can be rendered to customers in terms of the rate of on-time deliveries.

Constraints
6 Mathematical Problems in Engineering e first objective function (1) aims to minimise the total expected costs of the whole supply chain including, respectively, production costs in normal working hours and overtime, inventory-holding costs for both semifinished and finished products, transportation costs between plants, shortage costs, and distribution costs to customers. e second objective (2) seeks to maximise the customer satisfaction level by minimising the amount of products delivered with lateness. Constraints (3) and (4) are the flow balance constraint for the plants in production stages ST j , j � 1. . .N − 1, where production amounts are either stored in the inventory or further transported to its successor plants. Constraints (5) and (6) correspond to the inventory balance in plants from the last production stage.
ey are equal to the inventory of the previous period and the production amounts minus the shipped products to customers. Constraints (7) and (8) ensure that the quantity shipped from one production plant will be received by its successor. Equation (9) represents that the customer demand over the period must be satisfied or partially backordered. e backorder level of the product for the customer k in the period t is equal to the maximum difference between his/her demand and the quantity shipped from the plants from the last stage to customer zones taking into account distribution lead time. Constraint (10) ensures that the demands of customers for every product under each scenario and for the entire planning horizon are satisfied. In fact, the backorder level at the last period should be zero under the scenario s in order to avoid lost sales. Equation (11) states that the shipping amount of products from the plant for its successor is received in the next period. e limit of backorder's rate is illustrated in constraint (12). Equation (13) limits the production amounts to the random production capacity in regular time under the scenario s. e limitation of production capacity in overtime is written in equation (14). Constraint (15) limits the product inventory levels of production plants to their related inventory storage capacities. e shipping quantity of products between plants is delimited in equation (16) and the distribution capacity to customers in constraint (17). e constraint set (18) enforces the nonnegativity of decision variables.

Resolution Approach
In this paper, we applied one of the efficient multiobjective optimising solvers based on the well-known ε-constraint method as it is widely used in the multiobjective productiondistribution planning model. is approach is firstly defined in [31]. As shown in equations (19) and (20), the biobjective optimisation problem becomes a mono-objective optimisation problem by converting the second objective function into a constraint using the value ofε 2 .
e augmented ε-constraint method (AUGMECON) introduced in [32] is employed to generate the Pareto solutions of the multiobjective optimisation problem. e generated Pareto optimal solutions by the ε-constraint method may be dominated or inefficient solutions. e augmented ε-constraint technique is suggested to remedy this deficiency. In this method, the inequality constraint (20) of the ε-constraint approach is transformed into equality equation (21) by introducing a positive surplus variable S 2 . min subject to where r 2 is the range of the second objective function, c is a small positive number generally between 10 −3 and 10 −6 , and S 2 is the introduced surplus variable for the constraint.
To generate the payoff matrix, we should first optimise the two objective functions individually to obtain F * 1 and F * 2 . en, the feasible solution range of the considered objective 8 Mathematical Problems in Engineering functions would be determined with the help of the payoff table and equation (23). Next, we calculate ε k 2 by dividing the feasible solution r 2 into q distributed intervals, as depicted in equation (24). Finally, we optimise the selected objective function with the ε constraint and then repeat this step considering the different values of k to generate the set of Pareto optimal solutions.

Description of the Case Study.
e company considered in this study is called TE INTER, which is an export-oriented company located in Tunisia. e planning concerns a 3month time horizon, with each discrete period of one week. e supply chain network design is composed of 4 production stages, 6 manufacturers, and three customers, as illustrated in Figure 2.
e production process includes cutting of the fabric to obtain pattern pieces of the garment, embroidery/serigraphy of these pieces, assembly activity through clothes making process, and finally the packaging operations. TE INTER has the option of subcontracting some orders to local partners for the aim of covering the lack of technical competence and resources (I2) as well as extending the production capacities (I4 and I5). e finished products are packaged, batched, and distributed to 3 customers located in France. Different transportation modes could be used for delivery operations: trucks are used to tranship the products between plants and the seaport and/or airport for the exportation activities. Each transportation mode is characterized by its own cost, travel time, and capacities. e travel time is considered three periods (three weeks) for a ship and delivered at the same period for the airplane. For the sake of simplicity, we suppose that raw materials (fabrics and accessories) are fully available and with no defects at the beginning of the planning horizon.
Trendy apparel products are characterized by high demand uncertainty, quick response times, and a high risk of obsolescence. Forecasting the demand within the apparel and textile industry is a very difficult issue as demand is often highly volatile and cannot obviously be predicted. ese reasons make the decision-makers consider the demand uncertainty in the planning activities. Based on the historical data and experts of the company, we should mention that TE INTER acquires two types of orders from the potential customers: preseason orders and replenishment orders (also called the in-season orders). Preseason orders are that the initial customers' orders comprise the major part of the demand ordered prior to the start of the season. ey are firm orders and known with certainty at the beginning of the horizon plan. Accordingly, the due date of the production and the distribution planning of this deterministic demand is 3 months. However, it would be a difficult and hard task to predict the in-season orders because of the emergency of producing and delivering these orders, which serve to prevent shortages resulting from possible initial underestimation of demand or to replace unsold products during the selling season.
is is why retailers usually resort to inseason order opportunities after updating their forecasts based on the current sales and trends over the weeks of the current season. e quantity of the in-season orders would be much smaller than that of the preseason orders, and the associated lead time would not exceed few weeks. e second type of uncertainty considered in this industrial field is associated with the internal production capacity. is uncertainty depends on numerous factors and conditions. We should notice that manufacturing processes in the textile and apparel industry are labour-intensive as they are directly connected with the availability of workforce of each production site. In these circumstances, the absenteeism remains the most critical random indicator that would have a huge effect on the actual production planning. erefore, the actual production time capacity can be obtained by subtracting the absenteeism ratio based on the historical data of one year from total production time resources.

Data.
e considered data correspond to production collection established through 6 months of the summer season.
Overtime production and inventory-holding costs are proportional to regular production cost, following a respective factor of 1.4 and 0.2. Overtime production capacity is limited to 20% of the regular production according to the government policy. e distribution costs and capacity depend on the used transportation mode. e preseason demand, which is considered known with certainty at the beginning of the planning horizon, is given in Table 3. In addition, the production capacity without taking into account the absenteeism indicator is reported in Table 4. Based on the historical data of TE INTER, we consider that the customers' replenishment orders or the in-season orders and the production capacity parameters can be either very high, high, medium, or poor, i.e., four scenarios for each random parameter with associated probabilities based on the historical data of the summer season. e probabilities assigned to scenarios represent the importance of each scenario in an uncertain environment. e in-season demand could be estimated as a percentage of the preseason demand. We simply estimate the uncertain rate of workers' availability to define the production capacity by subtracting the ratio of workers' absenteeism and the random threshold of the in-season orders calculated as a percentage of the preseason orders (deterministic demand). e scenarios and the probabilities of the rate of workforce's availability and the ratio of the replenishment orders are depicted in Table 5.
Ratio of workforce's availability for each scenario � (1 − ratio of workforce's absenteeism for each scenario).
CAPP it,s � ratio of workforce's availability for each scenario * deterministic production capacity. RO pkt,s � ratio of the in-season orders for each scenario * preseason orders.

Computational Results
is section reports the results of the application of the proposed formulation to the industrial case study as well as    T1-T8  T9  T10  T11  T12   P1 T  T1  T2  T3  T4  T5  T6  T7  T8  T10  T11   some useful managerial insights that can be highly applicable to the fashion industry.

Performance of Stochastic Programming Solutions.
We use the augmented epsilon constraint (AUGMECON) approach to solve this multiobjective problem. e mathematical model has 1162 continuous decision variables and 822 constraints. e two-stage stochastic programming model and the sampling technique were implemented using LINGO 18 on PC Intel Core i5 with a 2.71 GHz processor and 8 GB memory.
We notice that combining these demand and production capacity uncertainty states resulted in the development of 16 scenarios shown in Table 6.
To assess the importance of the uncertainty effect, we analysed both the EVPI (expected value of perfect information) and the VSS (value of stochastic solution). e VSS is determined by the difference between the optimal solution of the recourse problem (RP) and the expected result of using the expected value solution (EEV) (equation (25). e latter is obtained by solving the RP with the first-stage decision variables fixed until (t 1 − 1) at the optimal values achieved by solving the average-scenario deterministic model using the mean values of the uncertain parameters. e VSS is a measure of how much can be saved by implementing the "here and now" solution as opposed to the deterministic expected value solution [33].
We evaluate the EVPI, as shown in equation: where RP is the objective value of the two-stage stochastic programming problem and WS = δ(s) * WS(s) is the expected wait-and-see solution, which expresses the objective value of the deterministic problem for the scenario s. Table 7 depicts the performance of stochastic programming solutions considering the implementation of the problem of the two-stage stochastic programming approach considering both random replenishment orders and production capacity. For each value of epsilon, we obtain relatively high values of the EVPI and VSS, which imply that it is important to consider uncertainty. e EVPI ranges from 24666.98 to 30861.44, which is approximately 6.691% to 8.091% of the two-stage stochastic objective value (RP), and e VSS ranges from 20851.97 to 35982.28, which is approximately 5.62% to 9.11% of the RP. Based on equations (25) and (26), diverse conclusions can be withdrawn by inspecting Table 8. In fact, the great values of the EVPI mean that incomplete information about the future may prove costly. Consequently, neglecting uncertainty in these cases deteriorates significantly the demand fulfillment policies and the capacity planning strategies. Besides, it is clear from the interpretation of the VSS that the solutions of the proposed two-stage stochastic programming model exceed those of the deterministic model. Accordingly, we point out that ignoring uncertainty via replacing the random parameters by their mean value leads to nonoptimised planning during making decisions. Furthermore, VSS/RP (%) has its highest value in case of the total customer satisfaction level (F2 = 1) which proves that considering uncertain factors has become more critical. ese results imply the potential worth of obtaining more accurate forecast as well as efficient stochastic solutions. Figure 3 presents the Pareto optimal solution of the RP and EEV and the solution provided by TE INTER. It clearly indicates that the considered objective functions are conflicting. It seems that an increase in one causes a decline in the other one because the on-time delivery loading effect will entail a strong percentage augmentation of distribution costs and overtime production costs, which strongly leads to decreasing profit. e total cost of the company's current practice is 421247.31 (monetary unit), and the customer satisfaction level is 56%. Our numerical results show that our stochastic approach not only provides better stochastic solutions (RP) than the expected value solution (EEV) but also is better than those proposed by the company's planner, as clearly depicted in Figure 3. In fact, we notice that the improvements stated between stochastic solutions (RP) and solution provided by the company are evaluated as 12.48% in terms of total costs and 28.77% in terms of the service level. To enhance the practitioners' understanding of these gains, it is important to consider the integration of production and distribution planning when demand and production capacity parameters are uncertain. Hence, the proposed model uses historical data and mathematical model that hold the ability to manage uncertain information. Anticipating absenteeism and demand helps firms achieve their targets, customer satisfaction, and increase productivity as well as flexibility and responsiveness of the supply chain. Figure 4 gives information about the total costs, regular production costs, overtime production costs, backorder costs, and transportation costs of each mode for the different scenarios (ε � 0.72). Based on the numerical results in Figure 4, we note that the total costs (F1) increase dramatically, while the customers' orders increase or the production capacity decreases. In fact, fixing, for example, of all "very high" in-season order scenarios (S1, S2, S3, and S4), we conclude that F1 increases when the production capacity is becoming lower. Taking into account the same production capacity scenario (for example, for medium scenario of production capacity: S3, S7, S11, and S15), it is obvious that overall costs (F1) increase in case of the increasing replenishment orders. e results reveal that the main contributions of the cost increasing in case of increasing quantity of the random orders urge the decision maker to react and plan more overtime production and use the airplane as a tool to speed up the production and the distribution activities. In conclusion, in order to hedge the huge effect of uncertainties, the practitioner would choose to procure either more overtime production or delivery planning using the combined sea-air transportation to fulfill the order delivery on time.

Sensitivity Analysis.
Here, we focus on a single-factor sensitivity analysis to study how the variation of the scenario probabilities of uncertain parameters affects the performance of the model.
We generate randomly 7 instances conserving the same supply chain structure of the industrial case study as we can  S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 S13 S14 S15 S16   see in Table 7. We firstly vary in instances 1-4 the probabilities of the in-season demand and consider the same probability of the production capacity of the case study. As for instance 5-6, we assume the same in-season demand probability and we consider a variation in production capacity's probability. In each instance, a scenario is considered to be far more likely than other economic scenarios. As shown in Figure 5, the total cost is more sensitive to the probability of uncertain demand compared with the probability of production capacity for instances from 1 to 6.
is observation could be explained by the impact of demand on the supply chain network performance. However, the cost augments slightly with the distribution of production capacity in these numerical examples. In instance 7, the cost increases dramatically assuming that the 4 th scenario of lowest production capacity has the highest probability.
is might be explained by planning more overtime production and air transportation.

Managerial Insights.
Demand uncertainty is the major challenge for any organisation in the current competitive world. In our case study, a part of the demand known as the in-season order is stochastic. In the real practice of the company, managers receive these unexpected orders during the selling season of the garments. e uncertain amount of these orders should be produced and delivered within a shorter time lapse to satisfy the customers and end users. On the other side, employee absenteeism is also identified as a crucial random factor that affects labour productivity and the quick response of the orders. Generally, the term "absenteeism" refers to unscheduled and unapproved absences of workers from the workplace. Based on our empirical findings, we provide decisionmakers with a pathway to follow for enhancing the level of the planning. e improvement area is achieved by a strong trade-off of the production in normal working hours and overtime, subcontracting, combined sea-air distribution, and backordering costs, whereas the current strategy based on the sequential planning method focuses solely on production. In terms of the customer service level, distribution network restructuring provides greater benefit to the apparel industry. In fact, the proposed model takes advantage of the airway transportation option to save time on trips, and this consequently increases on-time deliveries. To solve the problem of the absenteeism fluctuation, decision-makers can boost productivity by planning the annual leave taking into account the low and peak production season.
us, the managers would consider changing schedules when appropriate to accommodate different employee needs since the planning becomes more viable especially in case of high demand and low productivity. erefore, it is believed that these findings can guide decision-makers in making improvements with regard to costs and service level performance. is paper not only provides insights for policy-makers in terms of policy support but also contributes to the literature particularly within a competitive environment such as the textile-apparel industry.

Conclusion and Perspectives
is paper presents a new formulation of an integrated production-distribution planning problem in a stochastic multisite, multiproduct, multiperiod, and intermodal supply chain context. Two objectives are considered: minimising the total costs and maximising the customer satisfaction level in terms of orders delivered on time. A two-stage stochastic programming approach is developed to deal with demand and production capacity uncertainties. Results not only provide optimal production and distribution planning but also show the importance of considering uncertainties in the stochastic environment. A case study inspired from the manufacturing company producing trading clothes has been managed to illustrate and validate the stochastic model. Significant values of EVPI and VSS are obtained showing, respectively, the expected gain of perfect information and the outperformance of the stochastic solutions (RP) over the expected value solution (EEV). ese results justify the effectiveness of the proposed two-stage programming approach. Compared to real practice of the company, the observed improvements are evaluated as 12.48% in total costs and 28.77% in the customer service level. It is worth mentioning that the application of the proposed model may give some interesting managerial insights for decisionmakers regarding managing the fluctuation of the uncertainties of the apparel supply chain network. For instance, managers can settle a planning of annual leaves taking into account the peak and low production season in order to foster the production capacity and the flexibility of their plants.
For the future work, other sources of randomness affecting the manufacturing systems such as quality and supply uncertainties can be integrated. Another interesting issue consists of developing a risk management approach to ensure the robustness of the stochastic model. Finally, solving substantially larger production-distribution planning problems,  compared to those we addressed in our real case study, may call for the development of metaheuristics that can reach optimal solutions within reasonable lengths of time.
Data Availability e data provided by the company TE INTER used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.