Research Article  Open Access
Application of TwoPhase Fuzzy Optimization Approach to Multiproduct Multistage Integrated Production Planning with Linguistic Preference under Uncertainty
Abstract
This paper tackles the challenges for a production planning problem with linguistic preference on the objectives in an uncertain multiproduct multistage manufacturing environment. The uncertain sources are modelled by fuzzy sets and involve those induced by both the epistemic factors of process and external factors from customers and suppliers. A fuzzy multiobjective mixed integer programming model with different objective priorities is proposed to address the problem which attempts to simultaneously minimize the relevant operations cost and maximize the average safety stock holding level and the average service level. The epistemic and external uncertainty is simultaneously considered and formulated as flexible constraints. By defining the priority levels, a twophase fuzzy optimization approach is used to manage the preference extent and convert the original model into an auxiliary crisp one. Then a novel interactive solution approach is proposed to solve this problem. An industrial case originating from a steel rolling plant is applied to implement the proposed approach. The numerical results demonstrate the efficiency and feasibility to handle the linguistic preference and provide a compromised solution in an uncertain environment.
1. Introduction
Production planning in a multistage multiproduct production system concerns key production activities of each interconnected element. Material flows involved start from suppliers to the manufacturing network and finally to the customers. Coordination of lot sizes schedule for each element contributes to enhancing the utilization of activities that produce more values and thus attracts substantial interests. Typically, production planning can be categorized into three classes in terms of time horizon for decision making: longterm, mediumterm, and shortterm. Longterm planning focuses on identifying strategic decisions over a relatively long time, such as facility location and extent of additional investments in processing networks. Mediumterm planning involves tactical decisions for the optimal use of various resources and anticipated lot sizing for production stages. Models for mediumterm planning comparatively account for the presence of production system topology and concern controlling the material flow, inventory level, transports, suppliers, and so forth. Shortterm planning is related to detailed scheduling operations of a job, such as hourtohour job sequencing. The production planning problem in this work focuses on a mediumterm level which integrates activities going across the interconnection from suppliers to customers.
Normally, for traditional manufacturers, the material flows within the production system present in two common forms: parallel and spread. Parallel flows usually apply to singleproduct structure, and switching product families are rather expensive. In contrast, spread flows are related to diversification in product mix and can provide services for customization. However, the dynamics in market and customers challenge the manufacturers to response fast and accurately. It is usually difficult for the above two production forms to handle overachievement or underachievement for the forecasted demand of specific products, since the production stages within the system are either isolate or less correlate. When the demand for a specific product exceeds the production capacity of its dedicated production stages, it may result in backorders or lost sales, and on the contrary, it may also lead to a high level on inventory or excess capacities. Therefore, many manufacturers nowadays prefer to the production system in which products could be manufactured in a flexible way.
In this regard, we consider a multistage production system in which material and information flows may be interconnected or interwoven with each other. Other than the conventional parallel or spread topology structure, the addressed production system allows each production stage to schedule and allocate the items to be produced independently upon the information received from the customers. In other words, each product is produced through a designated route consisting of a series of production stages, and each stage may process items coming from one or more upstream stages.
To deal with demand fluctuation from the customers, several operations strategies have been developed and incorporated in planning decisions, such as backlogging, lost sales, or inventory holding [1–3]. These techniques show several advantages on obtaining feasible solutions in a more flexible way. However, in practical, the dynamic nature of production environment and markets induces a high degree of uncertainty, thus increasing the risk for decisionmakers. Peidro et al. [4] identified several metrics for distinct uncertain sources in supply chain planning and typically addressed the uncertainty for supply, process, and demand. In this paper, we consider three types of practical problems encountered for decisionmaking: epistemic uncertainty that originates from internal factors due to a lack of knowledge in production parameters or coefficients, external uncertainty concerned with customers and suppliers, and tradeoff among the conflicting objectives with decisionmakers’ preference which is derived from the complicated environment.
Several techniques and researches attempted to model the uncertainty in the production planning, among which four modelling approaches are identified: conceptual, analytical, artificial intelligence, and simulation models [5]. Among the researches, stochastic programming and probability distribution show the applicability in major production categories [6–8]. Stochastic programming formulates the uncertain sources based on the random characteristics and applies probability distributions that can be expressed in different forms such as normal, Beta, and Poisson [9]. Kira et al. [10] proposed a stochastic linear programming model to formulate the hierarchical production planning under random demand. Fleten and Kristoffersen [11] presented a multistage mixedinteger linear stochastic programming model to address a shortterm production plan for a hydropower plant. Zanjani et al. [12] developed a hybrid scenario tree with stationary probability distributions to integrate the uncertainty in the quality of raw materials and demand. However, as pointed out by Inuiguchi and Ramík [13], a stochastic programming problem may not be solved easily, since the computational efforts would be enhanced dramatically with the problem scale. Besides, when there is a lack of evidence on historical or statistical data, the probabilistic reasoning modelling methods are not always reliable and appropriate [14].
Alternatively, fuzzy set theory and possibility theory provide an attractive tool to account for the uncertainty [15]. There are several practical advantages of the fuzzy set theory with respect to handling the uncertainty for production decisions. It is able to deal with the uncertainty that is lack of data or evidence and is thus applicable to illdefined scenarios. It allows the decisionmakers to incorporate judgements for improving the solution interactively. The fuzzy model allows for development of more flexible and reliable decision tools that could be expressed linguistically based on human perception. The distinct fuzzy membership functions provide broader alternatives that are of high computational efficiency. Bellman and Zadeh [16] originally introduced the concept of fuzzy aggregation operators, upon which Zimmermann [17] presented a fuzzy linear programming approach to aggregating the fuzzy goals and constraints.
Following the development of fuzzy optimization approaches, the fuzzy set theory shows its effectiveness and superiority to handle the uncertainty in a production planning problem. Fung et al. [18] discussed a multiproduct aggregate production planning with fuzzy demand and capacity. They formulated the dynamic balance constraints as fuzzy/soft equations and solved the model using parametric programming. Wang and Liang [19] developed a fuzzy multiobjective linear programming model for a multiproduct aggregate production planning problem in a fuzzy environment. The fuzzy objectives were formulated using piecewise linear membership functions and were aggregated by a max operator. Wang and Liang [20] then addressed this problem by incorporating the imprecise forecasted demand, operating cost, and capacity. They adopted the triangular fuzzy number to represent the imprecise data and transform the fuzzy objectives by minimizing the three prominent points. Vasant [21] proposed a fuzzy linear programming model with a modified scurve membership function to solve a production planning problem with vague parameters and objective coefficients. Mula et al. [22] studied a material resource planning with flexible constraints in a multiproduct multilevel manufacturing environment. They applied a fuzzy linear programming approach with three kinds of aggregation schemes to decompose the original model into three fuzzy models. Torabi and Hassini [23] developed a novel multiobjective possibilistic mixed integer linear programming model for integrating procurement, production, and distribution planning. The model aimed to simultaneously minimize the total cost of logistics and maximize the total value of purchasing. Torabi et al. [24] further studied a fuzzy hierarchical production planning problem. The imprecise parameters along with the soft constraints are introduced to provide required consistency between decisions of the adjacent levels. Peidro et al. [4] proposed a fuzzy model for a tactical supply chain planning which contemplates the different uncertain sources from supply, demand, and process. They adopted a linear ranking function and triangular fuzzy numbers to transform the fuzzy model into a crisp equivalent problem. Taghizadeh et al. [25] presented a fuzzy multiobjective linear programming model for a multiperiod multiproduct production planning problem that simultaneously minimizes production cost and maximizes machine utilization. The model is solved using piecewise linear membership functions through a nonsymmetric decision for the fuzzy objectives and constraints. Aliev et al. [26] introduced a fuzzygeneric approach for solving an aggregate productiondistribution planning problem in a fuzzy environment. Generally, the above researches take advantage of various tools based on the concepts of fuzzy set theory that are applied in goal programming and soft/flexible mathematical programming. Other research studies relevant to fuzzy programming models include Hsu and Wang [27], Lan et al. [28], FigueroaGarcía et al. [29], Yaghin et al. [30], and Su and Lin [31].
In realworld applications, the production plan usually concerns conflicting objectives regarding the performance evaluated by various factors. Preference on the objectives imposed by the decisionmakers which is a critical issue is rather difficult to be represented in the conflicting objectives properly. Fuzzy set theory provides an alternative to express the preference; that is, the more important the objective, the higher the satisfying degree [32]. With respect to the production planning in an uncertain environment, however, few attentions are paid on the decisionmakers’ preference in the existing contributions in the literature. Conventionally, the weighted additive methods are used by assigning different weights for the objectives [33, 34]. However, the decisionmakers still encounter the problems of weighing the relative importance of each objective since the weights are quite dependent on the preference and are not easy to be specified upon the subjective judgment. With regard to the objective priorities, Tiwari et al. [35] constructed priority levels by specifying the membership grades and solved the subproblem for each level in sequence. Chen and Tsai [36] used the concept of desirable achievement degree to reflect the relative importance of goals explicitly; that is, the more important the goal, the higher the desirable achievement degree. However, these models are rigid in nature and may encounter the difficulties in obtaining a feasible solution when acquiring high overall satisfying degrees or desirable preference extent between the objectives.
In this paper, a novel fuzzy multiobjective mixed integer linear programming model (FMOMILP) is proposed for the multiproduct multistage production planning problem with linguistic preference in a fuzzy environment. The epistemic and external uncertainty is simultaneously formulated as flexible constraints and is handled by using a weighted average method and a fuzzy ranking method. The model defines two key performance indicators as the objectives, along with which the economic objectives are assessed in terms of relative importance. The model attempts to minimize total cost of production, overtime, raw materials, and inventory and in the meanwhile maximize the average safety stock holding level and the average service level. The production constraints include supply and production capacity, warehouse space, production route and product family, inventory balance, and backorders. Linguistic preference on the objectives imposed by the decisionmakers is considered and constructed by priority levels. To acquire desirable solutions that fully reflect the preference, a twophase interactive satisfying optimization method is applied through relaxing the overall satisfying degree. The FMOMILP model is transformed into an equivalent crisp model by treating the objectives and constraints separately. Then we developed a novel interactive solution approach to solving the FMOMILP model until a compromised solution is obtained. Finally, the proposed approach is implemented on realworld production planning in a steel rolling plant.
In summary, the main contributions of this work mainly lie in presenting a novel multiproduct multistage production planning model with linguistic preference and various practical constraints, simultaneously taking into account the epistemic and external uncertainty; definition and formulation of average safety stock holding level and average service level as fuzzy objectives for key indicators, together with fuzzy objectives for the operations cost, where priority levels are constructed; a twostage interactive solution approach to solving the fuzzy model which is able to treat not only the fuzzy objectives with linguistic preference but also the epistemic and external uncertainty in a flexible way.
The remainder of this work is organized as follows. Section 2 analyses the problem, details the assumptions, and formulates the FMOMILP model for the multiproduct multistage production planning problem. Subsequently, Section 3 presents the approach to handling the uncertainty and the objectives with linguistic preference and develops the twophase interactive solution approach. Next, an industrial case is used to illustrate the applicability and potential in Section 4. Finally, conclusions and remarks for further researches are given in Section 5.
2. Model Formulation
2.1. Problem Descriptions and Assumptions
The multiproduct multistage production planning problem addresses a small or medium sized plant that manufactures various types of products for the customers over a given planning horizon. Each product is processed stagebystage through a designated production route. As the flexible nature of the production topology, the production routes specified for different products may be intersected, thus complicating the modelling. To model the production topology, each production stage is formulated regarding its successor and predecessor stages. Since the operations parameters or coefficients are often incomplete and/or unobtainable, epistemic uncertainty from the productioninventory system involves the available production capacity, planned maintenance time, production efficiency, warehouse space, and safety stock level. These uncertain sources are modelled by fuzzy numbers using possibility theory. External uncertainty from both the suppliers and customers resulting volatile material and information flows is formulated as flexible constraints by fuzzy sets. The objectives evaluated by the decisionmakers can be divided into two different categories: operations cost and key performance indicators concerned. The proposed model aims to simultaneously minimize total cost related to production and inventory operations and maximize the average safety stock holding level and the average service level. In practice, such complex and conflicting objectives along with different constraints imposed by suppliers, production topology, and customers really challenge the efforts on decisionmaking. The assessment on each objective is vague when considering preference from the decisionmakers. The preference decisions are made in linguistic terms to distinguish the relative importance for the objectives, such that “average service level is the most important during this planning horizon” or “average safety stock holding level should be approximately higher than certain value and is relatively more important than inventory cost.” Therefore, this work also defines and incorporates the relative importance of the objectives in a fuzzy environment.
The proposed FMOMILP model is formulated based on the following assumptions:(1)Items produced can be stored to be used in certain periods, but within the planning horizon. The items are shipped to either the successor production stages or customers for delivery.(2)Backorders are allowed for multiple periods but should be fulfilled within the planning horizon. This strategy provides more flexibility for the production system to schedule the resources in a feasible way.(3)Supply capacity for each raw material is estimated and limited by the suppliers. The total procurement capacity for all the raw materials is limited by purchasing department of the manufacturer.(4)Suppliers, different production stages, and warehouses are assumed located in close proximity, so the lead time among them is ignored.(5)Raw materials, workinprocess, and enditems are stored in separated warehouses and thus formulated independently.(6)Some items can be acted as both workinprocess and enditems according to their usage with regard to the inventory level.(7)Dynamic safety stock strategy is considered by providing a buffer to meet fluctuating customer demand.(8)Priority levels are constructed for all the objectives to reflect decisionmakers’ judgement and preference on the relative importance of each objective. The priority structure may not be constant and is subject to the preference.
Indices, sets, parameters, and decision variables for the proposed FMOMILP model are defined as follows.
Indices are as follows: : index for production stages, : index for the immediate successor production stages of , : index for items including raw materials, workinprocess, and end items, : index for the immediate downstream items of , : index for item families, : index for time period, : index for raw material suppliers, : index for the periods after , : index for the periods before .
Sets are as follows: : set for all the items, : set for all the production stages, : set for time horizon, : subset of denoting raw materials, : subset of denoting workinprocess, : subset of denoting end items, : set for raw materials suppliers , : set for the production stages that produce , : set for the first production stage of , : set for the suppliers that supply raw material , : set for the terminal production stage of , : set for item families produced by production stage , : set for the items produced by production stage , : set for items involved in family , : set for the immediate predecessor production stages of for producing , : set for the items produced from by production stage .
Deterministic parameters are as follows: : hours in period , : cost of supplying one unit of raw material , : cost of backlogging one unit of end item for one time period, : cost of holding one unit of raw material in warehouse for one time period, : cost of holding one unit of workinprocess in warehouse for one time period, : cost of holding one unit of end item in warehouse for one time period, : overtime cost per hour in period , : cost of producing one unit of at stage , : setup cost of stage for producing , : large positive number, : planned regular production capacity of stage in period (in hours), : planned maintenance time of stage in period (in hours), : setup time for family at stage , : penalty coefficient for backorders that are fulfilled in the next planning horizon, : warehouse capacity occupied by one unit quantity of the item .
Fuzzy parameters are as follows: : production efficiency of stage for producing , : demand for end item in period , : supply capacity of raw material for supplier in period , : raw material supply capacity for supplier in period , : raw material warehouse capacity in period , : end item warehouse capacity in period , : workinprocess warehouse capacity in period , : end item safety stock in period .
Continuous variables are as follows: : inventory level for raw material in period , : inventory level for end item in period , : end item inventory deviation below the safety stock level in period , : inventory level for workinprocess , : minimum production quantity of family at stage in period , : backlogging quantity of end item in period that will be fulfilled in period , : quantity of end item that fails to meet the order demand in period , : supply quantity of raw material by supplier in period , : quantity of end item shipped to the warehouse in period , : production quantity of item at stage in period , : production quantity of family at stage in period , : overtime production capacity at stage in period , : quantity of raw material shipped to the warehouse in period and which will be stored to be used in period , : quantity of end item shipped to the warehouse in period and which will be stored to be used in period , : quantity of workinprocess shipped to the warehouse in period and which will be stored to be used in period , : maximum production quantity of family at stage in period .
Integer variables are as follows: : binary variable that indicates whether the end items shipped in period will be stored to be used in the future, : binary variable that indicates whether the end item demand in period will be backlogged, : binary variable that indicates whether item family is produced at stage in period .
2.2. Fuzzy Programming Model
2.2.1. Objective Functions
Two conflicting objective categories are considered simultaneously in the FMOMILP model: the operations cost and the key indicators concerned.
(1) Objectives for the Operations Cost. The operations cost includes those of production, overtime, raw materials, and inventory holding. Accordingly, the following describes these objective functions that are expected to be minimized:
The symbol “” represents the fuzzified version of “” and reads “essentially equal to.” The production cost in (1) is generated by processing and setup activities. The processing cost is calculated specific to each product and production stage independently. Setup cost is only related to changeovers between products of a family where similar products are grouped. Overtime production is allowed but relatively will incur high cost in labour, as defined in (2). Raw materials are purchased from diverse suppliers with different prices, as specified in (3). Additionally, inventory holding cost in (4) involves warehouses for raw materials, workinprocess, and enditems.
(2) Objectives for Key Indicators. Decisions on operations assessed by cost are those objectives that could be measured by economic performance easily and immediately. However, some objectives having indirect effects on the cost are not suitable to evaluate them by economic performance, such as customer satisfaction and planning nervousness [37]. Conventional researches attempted to estimate the parameters by cost and unify the measuring unit so that the multiobjective optimization problem can be managed through a simple additive method, interactive methods, or fuzzy methods [38–40]. Typical examples include estimation on backlogging penalty and lost sale penalty [41]. These estimations are usually empirical and may not satisfy the dynamic and diversified requirements since the production environment and market are rather volatile. With regard to these problems, this work formulates the objectives independently and carries out two key indicators as noneconomic measured objectives that are frequently encountered in the real application: average safety stock holding level and average service level. To incorporate the above indicators into the model, the following definitions are given:(i)Safety stock holding level (SHL) is as follows:(ii) Average safety stock holding level (ASHL) is as follows:
SHL evaluates deviation below the safety stock level for each product. Generally, ASHL measures the average level of safety stock holding for all the products, where indicates the number of elements in set .(iii) Service level (SL) is as follows:(iv) Average service level (ASL) is as follows:
SL concerns customer satisfactory in terms of backlogging level for each product. It should be noted that the backorders to be fulfilled in the next planning horizon have a greater impact on those backorders that are fulfilled during the current planning horizon. Thus, a weight penalty coefficient is added for this scenario. Similarly, ASL defines the average level on customer satisfactory for all the products.
Following the concepts of the above indicators, the decisionmakers seek to enhance the service and keep the safety stock level, which are presented by the objective functions:
2.2.2. Constraints
(1) Material Balance for Raw Material. The raw materials from various suppliers are stored in the warehouse to be used for further processing. In other words, material flow of raw material is controlled by both suppliers and initial production stages. The following constraints are applied to these two aspects for raw material balance:
Equations (11) and (12) ensure the material balance at raw material warehouse. Due to the operations strategy, raw materials are allowed to be stored for multiple periods within the planning horizon. Therefore, raw materials supplied during the last period of planning horizon would only be used for processing, as shown in (12). Equation (13) states that the raw materials stored in the warehouse are only used at initial production stages:
Equation (14) represents the procurement capacity limit from the perspective of manufactures. This restriction may be related to the shipment capacity or budget. Equation (15) implies that the supply capacity for each supplier is also limited by an upper bound. These two bounds are often lack of fully knowledge and thus are expressed as fuzzy coefficients.
(2) Material Balance for WorkinProcess. Similarly, material flow of workinprocess is controlled by its upstream and downstream flows, especially by its successor and predecessor production stages:
Equation (16) presents material balance of the workinprocess warehouse for period , in which two adjacent production stages are modelled through constructing sets with respect to their logical relationships. Equation (17) shows the material balance for the last period , where the items shipped to the warehouse should be used during the current period. On account of the production topology, items produced by a certain production stage may be either stored for further processing or shipped for delivery in terms of requirements. In this regard, (18) specifies that those items only for further processing are not allowed to be shipped to the end item warehouse.
(3) Material Balance for End Items. Consider
Equations (19) and (20) formulate the material balance at the end item warehouse for periods and , respectively. Safety stock is incorporated in the equations to provide a buffer against uncertainty, where the inventory level is allowed to fluctuate around the safety stock level. Hence, deviation below the safety stock is utilized to relax this restriction on the end item warehouse. Equation (21) indicates the relation of material flow into the end item warehouse with the items produced by terminal production stages. In realworld scenarios, product demand along with the safety stock is usually forecasted and estimated upon the historical statistics and cannot reflect the real market precisely. Therefore, these two variables are treated as fuzziness.
(4) Inventory Level Constraints. Due to the fact that raw materials, workinprocess, and end items are stored separately under different inventory management strategies, inventory constraints should be formulated independently upon the dynamic characteristics:
Since the model depicts exact usage of raw materials from one period to another, lower boundary of the raw material inventory level for periods and can be obtained by (22) and (23), respectively. The inventory level should be essentially less than a limited storage capacity as specified by (24):
Similarly, (25)–(27) imply the upper limit and lower limit for the workinprocess warehouse:
The inventory level for the end items should be appropriately controlled to satisfy the fluctuating demand by utilizing the flexible strategy of safety stock. Equation (28) ensures that the end item inventory level for period should cover both surplus parts of the safety stock and those end items stored for future use, and (29) states the same restriction for period . However, the buffer capacity of safety stock is also limited by the setting level itself, thus the lower deviation is restricted by (30). Equation (31) represents the storage capacity for the end item warehouse.
(5) Production Capacity Constraints. Consider
Various items are grouped in families by product and processing specification to reduce capacity consumption caused by frequent setup. Equation (32) defines attribution set and quantity relation to the families with the items. Equations (33) and (34) specify the logical relationship between binary setup and production quantity of the families:
The production capacity is consumed by both regular processing and setup. Equation (35) ensures that the total capacity occupied is not allowed to exceed the maximum level, which involves planned regular capacity and overtime capacity as well. Since the capacity is measured by time, the upper bound given by (35) should be further restricted by physical time and maintenance plans in a period, as shown by (36). Typically, (35) and (36) are modelled as flexible constraints due to the uncertainty existing in capacity allocation and efficiency.
(6) Backlogging Constraints. Consider
Although inventory holding and backorders may be sustained for several periods, for each product, they are not allowed to concurrent during a period. Equations (37)–(39) follow “the best I can do” policy; that is, the decisionmakers should make the most of inventory to satisfy the customers before considering backorders [42]. In particular, (37) and (38) convert the inventory holding and backorders into two binary variables, and then (39) can simply restrict the concurrence. It should be noted that, according to the operations strategy, the end items shipped to the warehouse during the last period of the planning horizon can be only used within the same period; however, backorders may be fulfilled during the next planning horizon. In other words, the concurrence would not exist during period and can be neglected.
(7) Decision Variable Constraints. Consider
Equations (40) and (41) indicate the integral constraints. Other decision variables default to nonnegativity as shown in (42).
2.3. Model Comparisons
The proposed FMOMILP model for the multiproduct multistage production planning problem is novel which considers meaningful technical features. The following Table 1 shows the qualitative comparisons of the proposed approach with some representative techniques for the multiproduct multistage production planning, including Torabi and Hassini [23], Peidro et al. [43], and Jamalnia and Soukhakian [44]. Several significant findings and advantages of the proposed approach can be summarized as follows. First, the proposed model involves several practical techniques to formulate the meaningful characteristics in the multiproduct multistage production planning problem, such as constructing the production route using set theory, allowing backorders for multiple periods, and holding the safety stock level. These technical factors provide more flexibility in dealing with the variants on production and market demand. Second, the proposed model defines two significant indicators ASHL and ASL as maximized objectives, together with the operations cost, which are treated to be fuzzy. The construction of ASHL and ASL presents the requirements from both manufactures and customers directly, instead of quantifying them by cost. Third, the multiple objectives have different priorities described by linguistic terms, which can fully express the preference from decisionmakers. The twostage interactive solution approach provides a flexible tool to deal with the objective priority where the preference extent can be quantified. Finally, the proposed model simultaneously considers the epistemic and external uncertainty which affects such a multiproduct multistage production environment, reflecting the uncertainty from productioninventory system, suppliers, and customers. Different strategies by fuzzy set theory and possibility theory are jointly developed to deal with the fuzzy constraints.

3. Solution Strategy
In the proposed problem, the uncertainty is expressed as fuzzy numbers to construct the FMOMILP model. In this section, the original FMOMILP model can be converted into an equivalent crisp multiobjective mixed integer linear model by treating the constraints in different forms independently. Then we apply a twophase fuzzy multiobjective optimization approach to handling the objectives taking into account linguistic preference.
3.1. Handling the Uncertain Sources
In the FMOMILP model, the fuzzy numbers represent the uncertainty on either one side or both sides of the constraints. The constraints with fuzzy numbers on one side hold a strict equality since the converted crisp constraints by defuzzing the uncertainty are still rigid. By contrast, those constraints with fuzzy numbers on both sides may be rather flexible due to the fact that the fuzzy numbers come from different sources.
In the literature, fuzzy numbers can be modelled as various membership functions in terms of application scenarios, such as triangular, trapezoidal, monotonic linear, piecewise linear, and Sshape [45]. Here, we adopt triangular membership function to construct the fuzzy numbers modelled in the constraints. As stated by Liang [46] and Li et al. [47], the triangular fuzzy number is commonly adopted due to its ease in simplistic data acquisition and flexibility of fuzzy arithmetic operations. Decisionmakers are easy to estimate the maximum and minimum deviations from the most likely value. Generally, a triangular fuzzy number can be characterized based on the following three prominent values [45]:(i)The most likely value that defines the most belonging member to the set of available values (membership degree = 1 if normalized).(ii)The most pessimistic value that defines the least belonging member to the set of available values on the minimum accepted value.(iii)The most optimistic value that defines the least belonging member to the set of available values on the maximum accepted value.
Considering the above definitions on prominent values, the membership functions for the triangular fuzzy numbers in the th fuzzy constraint can be represented as follows:where , , and refer to the most likely, most pessimistic, and most optimistic value for the th fuzzy constraint, respectively, and means the membership function for fuzzy number . If inducing the minimal acceptable degree , the prominent values can be reformulated as follows:where , , and denote the most likely, most pessimistic, and most optimistic value with the minimal acceptable degree, respectively. Then the membership functions can be reconstructed with respect to these three values.
Recalling the original model, (14), (15), (24), (27), (28), (29), (30), and (31) are those constraints with fuzzy numbers on one side. Tanaka and Asai [48] handled the triangular fuzzy number by using the weighted average of the most pessimistic value and most optimistic value. Alternatively, the weighted average method is utilized to treat the fuzzy numbers and construct crisp constraints [25, 49]. Then the above constraints can be converted into the following equations:where , , and represent the corresponding weights of the most likely, most pessimistic, and most optimistic value for the th fuzzy number, respectively. The most likely value with the highest membership should therefore be assigned to more weights than the others. When the weights for the above three prominent values are given, the solution can be further compromised by adjusting the minimal acceptable degree .
On the other hand, (19), (20), (35), and (36) represent the fuzzy numbers on both sides of the constraints. In this context, we treat the constraints based on comparison of their corresponding fuzzy numbers, which is referred to as fuzzy ranking. In the literature, several fuzzy ranking methods for comparing fuzzy numbers are given in terms of various classification schemes [50, 51]. We consider the concept of setinclusion along with the fuzzy ranking method proposed by Ramík and ímánek [52]. Then (19), (20), (35), and (36) can be converted into crisp constraints by ordering the lefthand and righthand side fuzzy numbers. The equivalent auxiliary constraints are expressed as
3.2. Treatment of the Objectives with Linguistic Preference
The objective functions in the FMOMILP model concern both the operations cost and the key indicators which are conflicting and noncommensurable. In such a fuzzy environment, decisionmakers are even more difficult to obtain feasible solutions as desired among the multiple objectives. Membership approaches and fuzzy set theory are efficient tools to reformulate the objective functions as commensurable units. Conventionally, the relative importance of the objectives is either specified with a weight vector or/and given satisfying degrees [23, 53] or compared by the corresponding satisfying degrees [44]. The above methods may lead to low computational efficiency, and additional constraints may be too strict to obtain satisfying even feasible solution. Thus, we apply the interactive satisfying method by Li and Hu [32] to overcome these shortcomings, where the satisfying degree is regarded as optimization variables and relaxed using varyingdomain optimization method [54]. The tradeoff among the conflicting objectives with linguistic preference can then be achieved by relaxing the overall satisfying degree through an interactive procedure. To do so, we implement a twophase fuzzy optimization approach to solving the proposed problem with linguistic preference.
3.2.1. Constructing the Priority Levels and Objective Membership Functions
Decisionmakers usually use linguistic terms to express the preference, such that objective A is more important than B. Considering the fact that the objective with higher priority has higher desirable achievement, the objectives involved can be grouped into different priority levels according to the preference. To quantify as priority levels, the FMOMILP problem can be constructed aswhere is a set for priority levels in sequence of relative importance; represents a set of objectives in the th priority level; and refers to system constraints.
By means of fuzzy set theory, the objective functions concerned can be generally expressed in the form of or , which means the objective is approximately less than or greater than the aspiration level , respectively. To describe the conflict in the objectives, a monotonically decreasing linear membership function is used to normalize the fuzzy objectives in (1)–(4):where for denotes the highest possible value with the membership degree of 1 and for is the smallest possible value with the membership degree of 0. Accordingly, the fuzzy objectives modelled by (9) are identified by fuzzy sets having an increasing linear membership function:where for denotes the highest possible value with the membership degree of 0 and for is the smallest possible value with the membership degree of 1.
3.2.2. Phase I Model
Phase I is to obtain the maximum overall satisfying degree for all the objectives without linguistic preference. Introduced by Bellman and Zadeh [16], the min/max operator can aggregate the objective membership functions using a new variable . According to the treatment on the fuzzy constraints and the objective membership functions aforementioned, the auxiliary formulation for the FMOMILP model regardless of linguistic preference can be given by
However, the optimal solution solved by the logical min operator in model (50) provides no preference judged by the decisionmakers. When considering the relative importance between the objectives, the subjective constraints induced may lead to a decline in the overall satisfying degree.
3.2.3. Phase II Model
In phase II, the maximum overall satisfying degree yielded in phase I will be relaxed to certain extent by the varyingdomain method, so that the linguistic preference could be involved. By inducing a relaxation parameter , the overall satisfying degree can be relaxed as shown by the following for the objectives :
The extent of relaxation corresponds to the priority imposed by . To compare the relative importance between the objectives in terms of preference, the varyingdomain method is applied by using a relaxation parameter to extend differences among the objective priorities. Concerning the priority levels defined in (47), the satisfying degree comparison can be stated bywhere means the satisfying degrees for the objectives in the th priority level. Equation (52) ensures that the objectives located in the higher priority level will hold higher satisfying degrees. In this regard, phase II also uses min operator to quantify the preference extent. Then the FMOMILP model can be further formulated by incorporating the priority levels and represented as follows:
In the above model (53), is considered to quantify the extent of the relative importance defined by the priority levels, which is relatively dependent on relaxation parameter . When the relaxation on is not enough, the preference expressed by the objective priorities may not be achieved, and in this case . Therefore, relaxation should be altered incrementally until the preference extent is achieved. This criterion relies on decisionmaker’s judgement, since the more the preference extent, the lower the overall satisfying degree.
3.3. Interactive Solution Procedure
The proposed approach concerns treating both the fuzzy constraints in different forms and incorporating the linguistic preference in the objectives. Decisionmakers can take fully advantage of the flexible parameters to achieve a preferred tradeoff between the overall satisfying degree and objective priorities. In summary, the original fuzzy model with linguistic preference is solved through a novel interactive method which incorporates decisionmakers’ judgement.
Step 1. Model the fuzzy numbers in the constraints using triangular possibility distributions by specifying the three prominent values.
Step 2. Given the value of minimal acceptable degree , convert the fuzzy constraints (14), (15), (24), (27), (28), (29), (30), and (31) into crisp ones using the weighted average method as shown in (45) and simultaneously, concert the fuzzy constraints (19), (20), (35), and (36) into crisp ones using the fuzzy ranking method.
Step 3. Model the fuzzy objectives and using the monotonically decreasing linear membership function in (48) and the monotonically increasing linear membership function in (49), respectively.
Step 4. Regardless of linguistic preference on the objectives, aggregate the fuzzy model using a min operator into a crisp phase I model. Then the corresponding overall satisfying degree can be obtained by solving the model as shown in model (50).
Step 5. Construct the priority levels for the objective functions in terms of linguistic preference.
Step 6. Relax the overall satisfying degree by (51) and formulate the priority comparison by (52).
Step 7. Assume that the relaxation parameter and aggregate the FMOMILP model and objective priorities into a crisp phase II model.
Step 8. Solve the crisp model (53), and the minimal value of preference extent can be obtained.
Step 9. Check whether the preference extent obtained in Step 7 satisfies the relative importance by the decisionmakers. If which means all the objectives are with the same importance and dissatisfy the linguistic preference, then go to Step 10. If but the preference extent is not satisfied by the decisionmakers, then go to Step 10 as well. Otherwise, stop the interactive procedure and the result obtained provides a compromise solution.
Step 10. Relax the overall satisfying degree obtained in Step 4 by augmenting the value of parameter , and then go back to Step 8.
4. Case Study
4.1. Case Description
This study investigates a steel rolling plant as a realworld case to illustrate the applicability of the proposed fuzzy programming model and evaluate the performance of the solution approach. The plant manufactures small and midsize steel plates classified in 10 types. These products are generally processed from 3 types of steel slabs as the raw materials, going through different production routes. The production system consists of 6 major production stages, and the production route for each product is constructed by a collection of these production stages according to the process. Table 2 presents the raw material and product matching matrix in which pair indicates that the product is manufactured by the raw material , and otherwise . Table 3 depicts the production stage and product matching matrix in which pair implies that the production of product goes through stage , and otherwise . It should be noted that the product is processed in sequence of , constituting a production route. Further, we can describe the workinprocess for each product upon its specific production route. In addition, the output items processed by production stage are indicated in terms of each product. Then in Table 3, for example, we can observe that product is manufactured through a production route 12456 with corresponding output workinprocess by each production stage. Concerning the production families modelled in (32)(33), Table 4 gives the family structure with regard to each production stage and its output items.



The case addresses a planning horizon of 7 periods. The uncertainty in the constraints modelled by the triangular fuzzy numbers is defined according to the most likely values with their deviations. Table 5 gives the relevant fuzzy parameters within a certain range of intervals in which the three prominent values by (43) are generated. The operational strategy implies that the production planning is made in a maketoorder and maketostock environment. It means that the products are manufactured based on the orders received and in the meantime to keep the safety stock. The dynamic demand is forecasted for each period prior to the planning horizon and is given in Table 6.

