ABSTRACT
This study puts forward a multiattribute decision model that determines good production and inventory parameter settings in maketostock (MTS) environments controlled by the Conwip (constant workinprocess) method. This model uses a discrete event simulation to evaluate the performance of a system in relation to workinprocess and a finished goods inventory. Based on MultiAttribute Value Theory (MAVT), a compromise solution is found by taking into account the decisionmaker’s preferences and tradeoff judgments as to the attributes of cycle time, throughput rate, holding cost and stockout cost. A numerical application is given to illustrate the use of the proposed multiattribute model which includes a sensitivity analysis.
Keywords:
MultiAttribute Value Theory (MAVT); tradeoff elicitation; simulation; Conwip
1 INTRODUCTION
This paper focuses on MTS production systems controlled by the CONWIP (Constant WorkinProcess) order release rule (Spearman et al., 199018. SPEARMAN ML, WOODRUFF DL & HOPP WJ. 1990. Conwip: a pull alternative to Kanban. Production Research, 28(5): 87994.; Hopp & Spearman, 20138. HOPP WJ & SPEARMAN ML. 2013. A Ciência na Fábrica (Factory Physics). 3rd ed. Porto Alegre: Bookman. 720 pages.; Gong et al., 20147. GONG Q, YANG Y & WANG S. 2014. Information and decisionmaking delays in MRP, KANBAN, and CONWIP. Production Economics, 156: 208213.). As demonstrated by Little’s Law (Little, 196112. LITTLE JDC. 1961. A Proof for the Queuing Formula: L = λW. Operations Research, 9: 83387.); Hopp & Spearman (2013), the performance of such systems can be highly influenced by different factors, including the WIP and finished goods inventory (FGI). This raises difficulties when the planning and control of production seeks to maximize performance in more than one attribute. Since cycle time (CT), throughput rate (TH), holding cost (HC) and stockout cost (SC) are the most common performance metrics used in MTS/CONWIP systems, the decision problem addressed in this paper consists of selecting the most appropriate level of WIP and FGI which best satisfies the preferences of the decisionmaker (DM) in relation to these attributes.
Although there is a multiple objective nature in this problem, little research has been undertaken that uses simulation together with multicriteria decision aiding (MCDM/A). Pergher & de Almeida (201715. PERGHER I & DE ALMEIDA AT. 2017. A multiattribute decision model for setting production planning parameters. Manufacturing Systems, 42: 224232.) present a multiattribute, expected utility model to define planning parameters in MTS/CONWIP production environments. For this kind of production system, Pergher & Vaccaro (201316. PERGHER I & VACCARO GLR. 2013. Work in process level definition: a method based on computer simulation and Electre tri. Produção, 24(3): 536547.) proposed a noncompensatory approach based on a simulation and the Electre tri method to determine the level of WIP. Another study using the Electre tri method to deal with stock management policies can be found in Fontana & Cavalcante (20136. FONTANA ME & CAVALCANTE CAV. 2013. Electre tri method used to storage location assignment into categories. Pesquisa Operacional, 33(2): 283303.). Lu et al. (201113. LU J, YANG T & WANG C. 2011. A lean pull system design analyzed by value stream mapping and multiple criteria decisionmaking method under demand uncertainty. Computer Integrated, 24(3): 211228.) proposed a lean pull system implementation procedure based on simulation, the Taguchi method and the TOPSIS multicriteria method. Xu et al. (201119. XU T, MOON DH & BAEK SG. 2011. Simulation study integrated with analytic hierarchy process (AHP) in an automotive manufacturing system. DOI: 10.1177/0037549711407781.
https://doi.org/10.1177/0037549711407781...
) applied simulation and the Analytic Hierarchy Process (AHP) to the design of a transmission case line in a Korean automotive factory. Another study using AHP to model the service and manufacturing activities in a supply chain can be found in Rabelo et al. (200717. RABELO L, ESKANDARI H, SHAALAN T & HELAL M. 2007. Value chain analysis using hybrid simulation and AHP. Production Economics, 105(2): 536547.).
The brief literature review given above indicates that there is a lack of studies on considering how the DM’s tradeoff assessments amongst criteria within the axiomatic structure of MAVT can be regarded. Consequently, the solution recommended may be inadequate when the DM wishes to compensate the low average performance of one of the decision attributes as a result of the high performance of another attribute. For this kind of decision context, a compensatory approach based on MAVT (Keeney & Raiffa, 19769. KEENEY RL & RAIFFA H. 1976. Decision Making with Multiple Objectives, Preferences, and Value Tradeoffs. New York: John Wiley & Sons. 592 pages.; Belton & Stewart, 20025. BELTON V & STEWART TJ. 2002. Multiple criteria decision analysis. Kluwer Academic Publishers. 372 pages.; de Almeida et al., 2015a1. DE ALMEIDA AT, CAVALCANTE CAV, ALENCAR MH, FERREIRA RJP, DE ALMEIDAFILHO AT & GARCEZ TV. 2015a. Multicriteria and Multiobjective Models for Risk, Reliability and Maintenance Decision Analysis. International Series in Operations Research & Management Science. New York: Springer.) is proposed. The solution for this MTS/CONWIP is found by taking into account the DM’s assessments among the multiple consequences of underestimating, and the potential consequences of overestimating, the level of WIP and FGI. The MAVT model is implemented with a Decision Support System (DSS) described in the application section.
The paper is structured as follows. In Section 2, the multiattribute decision model is proposed. Using realistic data, a numerical application is reported and discussed in Section 3. Finally, some conclusions are drawn and further lines of research are given in Section 4.
2 BUILDING THE DECISION MODEL
This section describes a framework (de Almeida, 2013a2. DE ALMEIDA AT, LUGO SDR, DOWSLEY BS, SANTOS VAPA & CLEMENTE TRN. 2015b. Elicitation for Tradeoff Additive Model with Linear or NonLinear Value Function and Sensitivity Analysis  Webbased  Código TUT2OMOWT1. Computational program developed by CDSID (Center for Decision Systems and Information Development). Registro INPI: BR5120160013214.; de Almeida et al., 2015a3. DE ALMEIDA AT. 2013a. Processo de Decisão nas Organizações. São Paulo: Atlas.) that is used to build multicriteria decision models on the basis of constructs obtained from the context of the decision problem itself. This framework allows the main characteristics of the method for developing simulation experiments proposed by Law & Kelton (200011. LAW AM & KELTON WD. 2000. Simulation modeling and analysis. 3rd ed. Boston: McGrawHill. 760 pages.) and the MAVT compensatory approach (Keeney & Raiffa, 19769. KEENEY RL & RAIFFA H. 1976. Decision Making with Multiple Objectives, Preferences, and Value Tradeoffs. New York: John Wiley & Sons. 592 pages.; de Almeida et al., 2015a4. DE ALMEIDA AT. 2013b. Additiveveto models for choice and ranking multicriteria decision problems. AsiaPacific Journal of Operational Research, 30(6): 120.) to be incorporated into a multiattribute decision approach. An adaptation of that framework, with details of the current application, is presented in Figure 1. The proposed model assumes that there is a single DM with a mutual preference independence condition amongst the criteria (Keeney & Raiffa, 1976).
In the preliminary phase, it is relevant to identify the actors who are involved in the decisionmaking process. The DM provides value judgments on multiple consequences in order to support the choice of an alternative that reflects the company’s objectives. Other relevant actors for this process are experts and a decision analyst. The experts provide technical information about the MTS/CONWIP production system required by the proposed model, whereas the decision analyst gives methodological support so that the proposed model can be applied. Furthermore, the actors must reach a consensus on the set of decision attributes that will be applied to evaluate the alternatives. Due to the characteristics of the problem, we consider TH, CT, HC and SC as the most relevant decision attributes for our investigation.
The first step is to collect the production and information flow data required for the simulation experiment. In Step 2, the production data are input into the MicroSaint Sharp simulation: a commercial, discrete event simulator adopted for this study. Afterwards, the simulation model should be validated in accordance with a structured technique (see Law & Kelton, 200011. LAW AM & KELTON WD. 2000. Simulation modeling and analysis. 3rd ed. Boston: McGrawHill. 760 pages.; MurraySmith, 201514. MURRAYSMITH DJ. 2015. Testing and Validation of Computer Simulation Models, Principles, Methods and Applications. Series: Simulation foundations, methods and applications. Springer: Cham. 252 pages.).
At Step 4, the set of alternatives A is developed based on WIP and FGI. This step also includes the simulation runs with set A considering the experimental parameters: number of replications R and warmup conditions. Thereafter, the consequence matrix is established to represent the performance of each alternative a _{i} ∈ A = {a _{1} , a _{2} , ..., a _{i} ... , a _{n} } over the attributes TH, CT, HC and SC. The procedure consists of calculating the average outcome x _{ij} for each combination of one alternative a _{i} and one attribute j, taking into account that the set of simulation outputs x is produced in R replications.
Step 6 consists of obtaining the subjective information from the DM within the axiomatic structure of MAVT. This preference modeling is assisted by a webbased DSS (de Almeida et al., 2015b) that performs the classical tradeoff elicitation protocol proposed by Keeney & Raiffa (19769. KEENEY RL & RAIFFA H. 1976. Decision Making with Multiple Objectives, Preferences, and Value Tradeoffs. New York: John Wiley & Sons. 592 pages.). This DSS was developed by CDSID (Center for Decision Systems and Information Development), which CDSID can made available on request (www.cdsid.org.br).
After inputting the required data set into the DSS, a sequence of questions is put to the DM, in order to obtain the indifference values between pairs of consequences associated to ordered attributes. The first group of questions is designed to rank the attribute weights. Then, other questions allow the DM to understand the consequence space. Finally, the DM should express the indifference value between pairs of consequences related to the ranking of attributes. As a result, the scaling constants k _{j} are obtained and used to aggregate the DM’s tradeoff judgments into a multiattribute value function (1).
where v(a _{i} ) is the overall value of alternative a _{i} and V _{j} (a _{i} ) represents the normalized value of a simulated consequence for attribute j . The following property holds (2):
The use of the additive form in (1) relies on the assumption that DM’s preferences over some attribute are not influenced by the level of another attribute (the mutual preference independence condition).
At step 7, the sensitivity analysis verifies to what extent the ranking of alternatives obtained by applying (1) is sensitive to variations in the scaling constants from the multiattribute function in (1). For this purpose, a Monte Carlo simulation and a Kendal t rank correlation (Kendall, 197010. KENDALL MG. 1970. Rank Correlation Methods. 4 ed. London: Griffin.) are applied to test the null hypothesis (h _{0}) that there is no correlation between each new ranking and the original one, h _{0}: τ = 0.
Finally, step 8 consists of analyzing the results from the sensitivity analysis and producing the report for the DM, with a detailed description of the final recommendations. In addition, a discussion on the main assumptions and simplifications of the model and its impact on the solution proposed should be included in the report.
3 DECISION MODEL APPLICATION
The case under investigation is an MTS/CONWIP multiproduct manufacturing line of a metalwork company. This system consists of five workstations W _{h} (h = 1, ... , 5) with a stable bottleneck and a storage area. In the storage area, finished products await orders from customers. If there are finished goods that match orders, the goods are dispatched to customers and this information is sent to W _{1} so as to update the planned order release. When the WIP falls below the target level, a product type is released to the system in accordance with the FIFO (firstin, firstout) dispatching rule. All product types (A, B, C and D) are manufactured in batch quantities. The statistical parameters of processing times and the interarrival time of customer orders are summarized in Table 1.
The plant works three 7h shifts a day, five days per week and is idle during the weekend unless onrequest. In relation to the simulation experiment, the following assumptions were made:

There is an infinite supply of raw materials.

Workstation W _{4} is the only bottleneck in the entire line.

Transfer time between workstations is negligible.

The scheduling of orders in the system is based on the FIFO rule.

The holding cost per producttype is determined by multiplying the average of inventoryonhand by an associated unit holding cost. Then, HC is determined as the sum of these partial values.

The stockout cost is determined by multiplying the number of canceled orders by the invariant gross profit yielded by each customer order.

Due to market characteristics, if the product is not available in sufficient quantity to meet an order in full, the customer cancels the order.
To validate the model, the comparison between the CT in the simulated period was considered. Adhesion between the real and the simulated system represented 86.1%. Considering the interarrival time of orders and the capacity of the storage unit, the simulation model was used to run different combinations of FGI and WIP. Trial runs reveal that 32 alternatives are necessary to explore the tradeoffs amongst criteria. Each treatment consists of 30 runs and the results were collected over 1680 hours (70 days of operations). The data from an initial warmup period of 24 hours were discarded. After making pilot runs, the consequence matrix was established to represent the average performance value of each combination of alternative a _{i} and attribute j. Table 2 shows the alternatives, their evaluation, and the consequence matrix.
As can be observed in the consequence matrix (Table 2), increasing the WIP level from 2 to 4 lots results in a better performance at TH. However, the CT performance is reduced. Likewise, increasing WIP from 4 to 5 lots does not yield a better consequence for TH and the performance of CT is the worst.
This performance behavior of MTS/CONWIP systems can be explained by Little’s Law (Little, 196112. LITTLE JDC. 1961. A Proof for the Queuing Formula: L = λW. Operations Research, 9: 83387.). In addition, alternatives that propose a WIP level of 2 lots yield the highest SC value, due to the low level of the use of capacity. However, these alternatives provide the lowest HC value. Alternatives a _{3}, a _{4}, a _{7}, a _{8} and a _{12} provide the most satisfactory performance for the SC attribute, while producing an extra holding cost in comparison to the alternatives associated with lower levels of FGI.
Assuming a DM with the mutual independence condition amongst criteria, the MAVT approach is applied in order to deal simultaneously with the conflicting nature of the decision attributes. After MAVT concepts have been explained to the DM, an assessment is made of his preferences by using the webbased DSS. This process starts by including the basic inputs to webDSS, namely: problem description, the set of attributes, preference direction to each attribute j (maximization or minimization) and consequence matrix. Subsequently, the rankorder of the weights of the attributes is elicited, using questions like: “consider a hypothetical alternative, with the worst performance in all criteria, and suppose you have to choose it. Now suppose that you can improve the performance of this alternative in only one of the criteria to the maximum value, which criterion would you choose?” As a result, the order obtained was: SC, HC, TH and CT.
After the ranking of criteria weights, the DSS provides a graphical visualization of the range of the consequences space $\left[{x}_{j}^{0},{x}_{j}^{*}\right]$, where ${x}_{j}^{0}$ and ${x}_{j}^{*}$ are the worst and the best outcomes respectively. At this point, other questions are put to the DM in order to provide a better understanding of the consequence space. Thereafter, the elicitation of k _{j} begins by asking the DM for indifferences between two consequences
that have been represented graphically. Then, the consequence ${x}_{SC}^{*}$ is reduced to the level of ${x}_{1}^{\text{'}}$ such that the DM is indifferent between ${x}_{1}^{\text{'}}$ and x _{2}. Based on these kinds of questions, the order three indifference values were elicited from the DM. In accordance with the DM’s preferential information, the DSS calculates the scale constants by using a simple system of equations related to the classical tradeoff method (see Keeney & Raiffa, 19769. KEENEY RL & RAIFFA H. 1976. Decision Making with Multiple Objectives, Preferences, and Value Tradeoffs. New York: John Wiley & Sons. 592 pages.; de Almeida et al., 2015a). Following this elicitation procedure, the values of k j obtained were:
By applying (1), this leads to the ranking shown in Table 3.
As shown in Table 3, a _{2} provides the highest value for the multiattribute value function (1) among the set of 32 alternatives examined. The proposed solution for this problem offers the most preferred tradeoff between the multidimensional consequences in the performance of the MTS/CONWIP system which is influenced by different levels of WIP and FGI. In comparison with the a _{6} , v(a _{6}) = 0.762, the compromise solution enables a better result in terms of HC and SC. From a single attribute perspective, the alternatives a _{3} and a _{9} provide the best results for SC and HC attributes, respectively. However, a _{3} increases HC to $3,973.73, while a _{9} increments the cost associated with lost sales in $3,388.91 in relation to a _{2}. Since there is no alternative providing the best results in all attributes modeled, the use of MAVT allows a compensating disadvantage in one of the performance attributes because of an advantage in other attributes, taking into account the goals that the DM wants to achieve.
A sensitivity analysis was conducted in order to examine the robustness of the compromise solution as to variations in the scaling constants. For this purpose, 10,000 independent ranks of alternatives were randomly generated from a uniform distribution with parameters (1 − 20%)k _{j} and (1+20%)k _{j} , where k _{j} is the scaling value elicited previously. Simulation results are summarized as follows: τ (max) = 1, τ (min) = 0.95, τ (mean) = 0.98, τ (median) = 0.99, τ (mode) = 0.99, τ (standard deviation) = 0.013. The pvalue of 0.078 was not statistically significant at a 5% level of significance. Thus, we can infer that the compromise solution a _{2} remains in first position in the ranking, since the simulated and the original ranks are independent. Although the compromise solution remains invariant within the (±20%) range, this does not mean that k _{j} can be elicited without careful evaluation of the mutual preference independence condition. In addition, the proposed MAVT model should be applied to MTS/CONWIP decision problems when the DM wishes to make compensations among multiple conflicting attributes. These can be obtained by the model given in (1). For situations in which the DM seeks to maximize the performance of the MTS/CONWIP system for a single objective, classical optimization approaches can be considered.
Due to the compensatory nature of MAVT, the values of the scale constants cannot be simply interpreted as a relative degree of importance of the attributes. The determination of k _{j} is based on subjective information provided by the DM, considering the range of the consequences space.
That is, when the range $\left[{x}_{j}^{0},{x}_{j}^{*}\right]$ changes, it is necessary to adjust the value of k _{j} . For instance, by varying the values of the consequence matrix shown in Table 2 within the range (1 − 25%)x _{ij} and (1 + 25%)x _{ij} without adjusting k _{j} , the Kendall tau test rejects h _{0}: t = 0 ( pvalue of 0.034), which indicates significant differences between ranks evaluated at a 5% level of significance. This feature is a distinctive feature of this multiattribute model in relation to the models proposed hitherto. However, the proposed model assumes that the DM is willing to answer all questions required by the tradeoff method (Keeney & Raiffa, 19769. KEENEY RL & RAIFFA H. 1976. Decision Making with Multiple Objectives, Preferences, and Value Tradeoffs. New York: John Wiley & Sons. 592 pages.).
To address the uncertainties of the performance of the MTS/CONWIP system, Pergher and de Almeida (201715. PERGHER I & DE ALMEIDA AT. 2017. A multiattribute decision model for setting production planning parameters. Manufacturing Systems, 42: 224232.) integrated discrete event simulation with Multiattribute Utility Theory  MAUT (Keeney & Raiffa, 19769. KEENEY RL & RAIFFA H. 1976. Decision Making with Multiple Objectives, Preferences, and Value Tradeoffs. New York: John Wiley & Sons. 592 pages.). Nonetheless, this paper assumes that the average value of consequences provides a satisfactory representation of the performance of the system and hence the use of MAVT is suitable in this context.
4 CONCLUDING REMARKS AND FUTURE RESEARCH
In this paper, a decision model is presented that integrates discrete event simulation with an MAVT compensatory approach so as to identify good settings for production/inventory parameters, namely cycle time (CT), throughput rate (TH), holding cost (HC) and stockout cost (SC).
We use a webbased DSS developed by the CDSID to implement the MAVT compensatory approach. The solution recommended for this multiattribute problem is a _{4}, which reflects the DM’s tradeoff judgment among multiple conflicting attributes.
The use of a discrete event simulation and MAVT allows not only the economic results to be evaluated but also those intangible consequences yielded by the level of WIP and FGI, without stopping the real system.
The proposed model can be further improved, for example by considering the additiveveto model (de Almeida, 2013b). For decision problems where the DM does not wish to select one alternative with a very low performance in one attribute, even though it is compensated by the high level of one or more of the other attributes, additiveveto models can be employed.
ACKNOWLEDGMENTS
This work had partial support from CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico), the Brazilian Research Council. The CNPq was not involved in the study or in writing this paper. The authors would like to thank the Editor and the anonymous reviewers for their insightful critique of an earlier version of this paper and for their valuable suggestions for improving it.
REFERENCES

^{1}DE ALMEIDA AT, CAVALCANTE CAV, ALENCAR MH, FERREIRA RJP, DE ALMEIDAFILHO AT & GARCEZ TV. 2015a. Multicriteria and Multiobjective Models for Risk, Reliability and Maintenance Decision Analysis. International Series in Operations Research & Management Science. New York: Springer.

^{2}DE ALMEIDA AT, LUGO SDR, DOWSLEY BS, SANTOS VAPA & CLEMENTE TRN. 2015b. Elicitation for Tradeoff Additive Model with Linear or NonLinear Value Function and Sensitivity Analysis  Webbased  Código TUT2OMOWT1. Computational program developed by CDSID (Center for Decision Systems and Information Development). Registro INPI: BR5120160013214.

^{3}DE ALMEIDA AT. 2013a. Processo de Decisão nas Organizações. São Paulo: Atlas.

^{4}DE ALMEIDA AT. 2013b. Additiveveto models for choice and ranking multicriteria decision problems. AsiaPacific Journal of Operational Research, 30(6): 120.

^{5}BELTON V & STEWART TJ. 2002. Multiple criteria decision analysis. Kluwer Academic Publishers. 372 pages.

^{6}FONTANA ME & CAVALCANTE CAV. 2013. Electre tri method used to storage location assignment into categories. Pesquisa Operacional, 33(2): 283303.

^{7}GONG Q, YANG Y & WANG S. 2014. Information and decisionmaking delays in MRP, KANBAN, and CONWIP. Production Economics, 156: 208213.

^{8}HOPP WJ & SPEARMAN ML. 2013. A Ciência na Fábrica (Factory Physics). 3rd ed. Porto Alegre: Bookman. 720 pages.

^{9}KEENEY RL & RAIFFA H. 1976. Decision Making with Multiple Objectives, Preferences, and Value Tradeoffs. New York: John Wiley & Sons. 592 pages.

^{10}KENDALL MG. 1970. Rank Correlation Methods. 4 ed. London: Griffin.

^{11}LAW AM & KELTON WD. 2000. Simulation modeling and analysis. 3rd ed. Boston: McGrawHill. 760 pages.

^{12}LITTLE JDC. 1961. A Proof for the Queuing Formula: L = λW. Operations Research, 9: 83387.

^{13}LU J, YANG T & WANG C. 2011. A lean pull system design analyzed by value stream mapping and multiple criteria decisionmaking method under demand uncertainty. Computer Integrated, 24(3): 211228.

^{14}MURRAYSMITH DJ. 2015. Testing and Validation of Computer Simulation Models, Principles, Methods and Applications. Series: Simulation foundations, methods and applications. Springer: Cham. 252 pages.

^{15}PERGHER I & DE ALMEIDA AT. 2017. A multiattribute decision model for setting production planning parameters. Manufacturing Systems, 42: 224232.

^{16}PERGHER I & VACCARO GLR. 2013. Work in process level definition: a method based on computer simulation and Electre tri. Produção, 24(3): 536547.

^{17}RABELO L, ESKANDARI H, SHAALAN T & HELAL M. 2007. Value chain analysis using hybrid simulation and AHP. Production Economics, 105(2): 536547.

^{18}SPEARMAN ML, WOODRUFF DL & HOPP WJ. 1990. Conwip: a pull alternative to Kanban. Production Research, 28(5): 87994.

^{19}XU T, MOON DH & BAEK SG. 2011. Simulation study integrated with analytic hierarchy process (AHP) in an automotive manufacturing system. DOI: 10.1177/0037549711407781.
» https://doi.org/10.1177/0037549711407781
Publication Dates

Publication in this collection
JanApr 2018
History

Received
24 May 2017 
Accepted
16 Dec 2017