Research Article  Open Access
Solving Fuzzy Multiproduct Aggregate Production Planning Problems Based on Extension Principle
Abstract
Aggregate production planning (APP) plays a critical role in supply chain management (SCM). This paper investigates multiproduct, multiperiod APP problems with several distinct types of fuzzy uncertainties. In contrast to the existing studies, the modelling in this work conserves the fuzziness such that the obtained APP is more effective. Based on Zadeh’s extension principle, the results obtained are fuzzy solutions described by membership functions, in contrast to results from previous studies. A pair of twolevel parametric mathematical programs is formulated to calculate the lower and upper bounds of the optimum fuzzy performance measure. The membership function of the fuzzy total cost is constructed by enumerating various possibility levels. A case studied in previous research is investigated to demonstrate the validity of the proposed model and solution procedure. Because the optimal objective value and associated decision variables are expressed using fuzzy numbers rather than crisp values, the proposed approach is able to represent APP systems more accurately, and therefore, the results obtained can provide decision makers with more effective and informative APPs and more chance to achieve the optimal disaggregate plan.
1. Introduction
Current trends in the highly competitive and dynamic business environment are driving companies across the globe to move towards aggregate production planning (APP) with the intent of finding optimum balance among capacity, forecasted demand, and fluctuating customer orders over the midterm, often from 3 to 18 months ahead. An ample body of research has indicated that this optimum is achieved through management that includes the primary objective of adjusting the controllable factors (e.g., regular and overtime production rates, inventory levels, and labour levels) under a given set of production resources and constraints during each phase of the planning horizon [1–5].
The APP problem has been extensively investigated over the last few decades, and a significant amount of effort has been expended in the development of new decision models. By reviewing previous studies, Saad [6] examined the conventional decision models and classified them into the following six categories: the linear decision rule (LDR) [7], the transportation method [8], linear programming (LP) [9], the management coefficient approach [10], simulation approaches [11], and the search decision rule [12]. In the following decades, studies were focused on developing a variety of optimisation models [1], and more recent studies include those of Baykasoglu [13], Jain and Palekar [12], Gomes da Silva et al. [14], Leung et al. [15], and Leung and Chan [16].
The development of the abovementioned models took place under crisp environments, and the resulting models belong to deterministic or stochastic optimisation methods [1, 17, 18]. However, Mula et al. [18] also noted that the APP problem is associated with elements that raise the uncertainty of the nonstochastic elements (e.g., customer demand, lead times, production fluctuation, and their associated costs). Many attempts have been made to categorise problems of this type. In practical applications, in place of probability distributions, such linguistic terms as “approximately 1.6 meters,” “often,” “moderate,” and “rare” may be embedded in APP, and this information is obtained subjectively from management [19, 20].
Fuzzy set theory was first introduced by Zadeh [21] to handle uncertainties in the nonstochastic sense. Because the APP problem in crisp environments can be modelled as an LP model, several fuzzy optimisation methods for solving fuzzy APP problems have been developed based on the fuzzy LP (FLP) [22–28]. For example, Lee [29] investigated singleproduct APP problems with fuzzy objectives, workforce levels, and demands. Other studies include the works by Lai and Hwang [30], Dai et al. [31], Wang and Fang [28, 32], and Wang and Liang [19, 33].
In practice, most companies produce multiple products, and thus, many fuzzy parameters of different types exist. Therefore, practical multiproduct APP problems are more complicated than the singleproduct APP often presented in modelling and solution procedures, and these more complicated problems deserve further investigation. For multiproduct APP problems, Tang et al. [34] developed a model characterised by fuzzy resource righthandside constraints (i.e., demands and capacities). Soon after, Wang and Fang [32] discussed a multiobjective APP problem in which the objective function coefficients (i.e., production prices and subcontracting costs) and resource righthandside constraints (i.e., manpower levels, production capacities, and market demands) were characterised by fuzziness values that were in turn included in the model as trapezoidal fuzzy numbers. Wang and Liang [33] and Liang [35] applied a possibilistic linear programming (PLP) [30] approach to solving multiproduct APP problems associated with multiple imprecise objectives and cost coefficients, and the demands were described by triangular possibility distributions in uncertain environments. Aliev et al. [36] proposed an approach based on the genetic algorithm (GA) for the fuzzy multiproduct aggregate productiondistribution planning (APDP) problem in supply chain management (SCM).
Tang et al. [34] emphasised that any single optimal solution should not be taken as a guarantee of reaching a viable disaggregate plan. Furthermore, Chen and Huang [37] reviewed the models and solution methods for APP problems and noted that most of the existing methods may fail to maintain the fuzziness of the input information, which can be used to represent the fuzzy APP more accurately. It thus becomes apparent that the performance measures are characterised by fuzziness because multiple parameters in an APP model are also fuzzy. Furthermore, conserving the fuzziness of the input information is possible if the model offers the opportunity to derive descriptive membership functions of certain performance measures. Better representation of the fuzziness will thus facilitate the user in attaining a performance measure that is a more logical and more accurate representation of the fuzzy APP. The abovementioned approaches offer crisp solutions that are prone to either under or overestimates, thus leading users to make faulty decisions.
To overcome this shortcoming, Chen and Huang [37] proposed a solution procedure that is able to find the fuzzy objective value of the fuzzy APP model. However, they only investigated the singleproduct APP problem with two sets of fuzzy parameters, the maximum workforce available and the forecasted demand. In fact, their model is a FLP model with crisp objective coefficients, crisp lefthandside constraints, and fuzzy righthand sides with two sets of constraints. As stated previously, multiproduct APP problems are more complicated than their singleproduct counterparts. This paper investigates the multiproduct, multiperiod APP problem with several fuzzy parameters of different types. More importantly, the solutions obtained by solving the proposed model are also fuzzy and thus conserve fuzziness in the general nature, indicating that the proposed model can provide an effective APP and a feasible disaggregate plan can be obtained.
This paper proposes a costbased multiproduct/multiperiod APP model with several sets of fuzzy parameters, including the unit production costs excluding labour costs, overtime labour costs, labour costs, inventory carrying costs, costs to hire one worker, costs to lay off one worker, unit backorder costs, conversion factors in hours of labour, forecasted demands, and maximum workforce. A solution procedure that is able to calculate the fuzzy objective value of the fuzzy APP model is also developed. The main concept is based on the application of αcuts and Zadeh’s extension principle [38, 39] to transform the fuzzy APP model into a family of crisp APP models. A pair of twolevel mathematical programs [40, 41] is formulated to calculate the lower and upper bounds of the αcut of the fuzzy minimum total cost. Consequently, the membership function of the fuzzy minimum total cost is derived analytically or numerically by enumerating different values of α. Because the proposed approach is based on Zadeh’s extension principle, it is significantly different from those of several related studies (e.g., [30, 31]), which may fail to compute the sets of possible values of the fuzzy minimum total cost.
This paper is organised as follows. First, an extensive fuzzy multiproduct APP model with multiple fuzzy parameters is constructed in Section 2 based on the work of Lai and Hwang [30]. Based on Zadeh’s extension principle, Section 3 presents a pair of mathematical programs for calculating the αcuts of the fuzzy minimum total cost. Next, a modified twoproduct sixperiod APP example referenced by Lai and Hwang [30] is successfully investigated in Section 4 to demonstrate the validity of the proposed approach and is compared with one of the classical methods, that is, Chanas’ approach [30], in Section 5. Finally, conclusions are presented in Section 6.
2. Modelling the Fuzzy Multiproduct APP
This paper constructs the multiproduct APP model with many fuzzy parameters of different types based on the LP formulation. In this section, the crisp LP and general fuzzy LP models are briefly introduced.
2.1. Fuzzy Linear Programming
The LP method has been demonstrated as one of the most frequently used operational research/management science (OR/MS) techniques for practical problems [42]. The LP method focuses on the effective and efficient allocation of limited resources to interrelated activities with the objective of achieving a prespecified goal, such as profit maximisation [43]. The general maximisation LP model can be described in the canonical form using matrix notation [44]: where denotes the transpose of the vector of the profit coefficients of the objective function, A is the constraint matrix, x is the vector of decision variables (or activity levels), and b is the vector of total resources available. Note that in this model, all coefficients of A, b, and c are crisp numbers and each constraint must be strictly satisfied. In practical situations, the input data for this model are typically imprecise due to incomplete information [30]; in addition, the vagueness in these input data may not be probabilistic. As Rommelfanger [27] noted, “LP requires much welldefined and precise data which involves high information costs,” and its assumptions [42] are often too strict to limit the LP applications.
To quantitatively address the imprecise information in the decisionmaking process, Zadeh [21] introduced the concept of fuzzy set theory. Since then, additional research has developed specific FLP models and proposed the associated solution methods. In fact, so many different types of FLP models exist that Lai and Hwang [30] provided a detailed and systematic classification. In general, the basic types of FLP can be grouped into the following six categories: fuzzy objective and precise constraints, precise objective and fuzzy constraint matrix, precise objective and fuzzy total resources available, fuzzy objective and fuzzy constraint matrix, fuzzy objective and fuzzy total resources available, and fuzzy objective and fuzzy constraints, including the constraint matrix and total resources available [26]. In other words, the possible combinations of the fuzziness of c, A, and b are fuzzy c only, fuzzy A only, fuzzy b only, fuzzy c and A, fuzzy c and b, and fuzzy c, A, and b. Numerous studies have been devoted to developing different types of solution methods for these six types of FLP models.
Instead of describing all of the types of FLP models, the general FLP model with fuzzy c, A, and b is exemplified as follows: Clearly, the assumption of certainty for the crisp LP is relaxed to fit realworld situations. Moreover, the operations of addition and multiplication are operations of fuzzy arithmetic, and “” denotes the ordering of fuzzy numbers [26]. Note that the proposed model in the next subsection is of this type.
Lai and Hwang [30] also recognised that fuzzy linear programming (FLP) problems differ from possibilistic linear programming (PLP) problems because FLP connects fuzzy input data through a subjective preferencebased membership function, while PLP correlates imprecise data by means of possibility distributions in which the latter are an analogue of probability distributions that are either objective or subjective by nature.
2.2. Fuzzy Multiproduct APP
The investigated APP problem in this paper originates from the study previously conducted by Lai and Hwang [30]. We consider a multiproduct manufacturer that aims to satisfy the market demand over a medium planning horizon . The management aims to efficiently satisfy the forecasted demand and thus adopts APP in an attempt to determine the optimum of inventory levels, production rates, overtime work, subcontracting rates, and other factors. In particular, the challenge is to determine a feasible method of adjusting the above factors appropriately such that the total costs are minimised and the market demand is satisfied. The forecasted demand may either be met or backordered [30]. At least one of the backorders and inventory levels in a period must be zero because it is crucial that any backorder be met within the following period. Furthermore, during any period, either net hiring or net firing may occur but not both. A fuzzy decision model is thus constructed to solve this problem following a subset of Lai and Hwang’s assumptions [30].(1)The product quantities include those of regulartime production, overtime production, production due to hiring additional employees, and those required to meet the demands.(2)The initial inventory refers to the level of inventory that is necessary for initiating an order.(3)The initial level of workforce for the production is known to the decision maker.(4)The processing of successive orders does not induce setup costs.(5)The finished goods in process can be stored in the inventory storage space.(6)The new employees are equally productive as the old employees, indicating the presence of a reliable workforce pool. The new employees are assumed to be as fully productive as the old employees when they begin work.
Tables 1 and 2 list the decision variables and parameters of the proposed model. Notably, the first 10 sets of parameters listed in Table 2 are assumed to be fuzzy without loss of generality.


2.2.1. Objective Function
The objective is to minimise the total incurred cost including the following:(i)the total production cost: ;(ii)the total labour cost: , where is a conversion factor for transforming the unit of to manhours (refer to Table 2);(iii)the total inventory carrying cost: ;(iv)the total backorder cost: ;(v)the total costs of changes in labour levels, including the costs to hire and layoff workers: , where , indicating that either net hiring or net firing of labour occurs during a period but not both [30].
2.2.2. General Constraints
The following constraints are considered for each product and time period.
(1) Labour Level Constraints.(i)The workforce level should not be greater than the maximum available workforce level during any period: .(ii)The workforce level in this period should be equal to the workforce level in the preceding period plus the new hires minus the layoffs : , where .(iii)The regulartime production capacity should not be larger than the available labour capacity: , where and are conversion factors for transforming a unit of and to manhours, respectively (refer to Table 2).(iv)The variation in a workforce should not exceed the permitted level of the company policy during any period: , where and are the same as those described above.
(2) Capacity Constraints. At the end of each time period, either a net inventory or backorder will occur but not both. Therefore, the statement that the total available product quantity should be equal to or greater than the minimum demand in period can be written as where , indicating that the demand over a particular period can be either met or backordered but not both.
(3) Inventory Level Constraints. In a similar manner, either a net inventory or backorder will occur during a period but not both. Therefore, the fundamental material quantity balance [30, 31, 45] can be written as where . This equality indicates that the inventory level or backorder level in this period is equal to that of the preceding period plus the regulartime and overtime production and minus the forecasted demand. This equation ensures that the amount of each product sold in a period plus the inventory (or backorder) at the end of the period equals the total supply consisting of inventory (or backorder) from the previous period plus the regular and overtime production in the current period [37, 46].
2.2.3. The Proposed Model
By incorporating the above objective function, general constraints, and nonnegative constraints , the proposed fuzzy multiproduct APP model can be constructed as follows [30]:
It is clear that Model (5) is intrinsically a linear programming (LP) model. Therefore, no requirement exists to enforce two sets of constraints, and , on Model (5) because according to the complementary slackness theory of LP [44], the basic solution to Model (5) has a special structure that does not allow and to exist as basic variables at the same time, as is also true for and .
3. The Solution Procedure
Because Model (5) is a FLP model, it can be solved using the existing FLP solution methods. However, as is generally known, most of the results obtained from these methods are crisp rather than fuzzy, an outcome that will result in the loss of certain important fuzzy information and undermine the quality of the APP. To overcome this shortcoming and conserve the fuzzy information, Chen and Huang [37] proposed a solution method for singleproduct APP problems with fuzzy parameters on the righthand side of the general constraints. However, Model (5), an FLP with several parameters in the objective function and the coefficients of the general constraints, is more complicated than the model conducted by Chen and Huang [37]; therefore, it deserves further investigation.
In Model (5), many parameters are approximately known and can be described by fuzzy numbers [26]. For example, the forecasted demands () are stated as where is the membership function and is the support of the fuzzy number that denotes the universal set of the forecasted demands. Other fuzzy parameters can be described by the following fuzzy numbers in the same manner: According to Zadeh’s extension principle, if certain parameters in Model (5) are fuzzy numbers, the minimum of the objective function is also a fuzzy number, and its membership function is defined as follows: where is the crisp objective value of Model (5) with the fuzzy parameters degenerated simultaneously to their crisp values. Equation (8) is theoretically correct, but it is not easy to use in practice.
Instead of directly using (8), this present paper uses the concept of parametric mathematical programming to find the membership function . A pair of mathematical programs is developed to find the cuts of based on Zadeh’s extension principle. First, the definition of the cut of a fuzzy number is briefly introduced. The cut of the forecasted demand is defined as follows [26]: Using , can be represented by different possibility levels of the confidence intervals [26]. Therefore, the cut of the fuzzy number defined in (9) represent a crisp interval that can be expressed in the following form: This interval indicates where the value of lies at possibility level . Clearly, from Definition (8), the membership function is also parameterised by . Thus, one approach to constructing is the cut approach.
Next, the of other fuzzy parameters can be, respectively, described as follows:
The cut approach is intended to derive the cuts of or, in other words, to derive the crisp interval for all . This derivation will be introduced later. According to Zadeh’s extension principle [38, 39], the defined in (8) is a fuzzy number that possesses convexity [26]. Consequently, according to the attractive feature that all cuts form a nested structure with respect to [47], the set of intervals builds the approximated membership function of by enumerating different values of .
According to Definition (8), is the minimum of these 10 sets of membership functions of fuzzy parameters (e.g., and ), which complicates the situation. Thus, to derive the crisp interval of and address the membership value, all of the membership values of these 10 sets of fuzzy parameters that are not smaller than are required, and at least one of them must be equal to such that to satisfy . It is clear that the membership function can be defined by determining the leftshape function and the rightshape function of . Based on this idea, the lower bound and upper bound of can be derived using the solution of the following pair of twolevel mathematical programs [40, 41]:
To satisfy as required by Definition (8), at least one of , and must meet the boundary of their αcuts. For Model (12a), can be directly set to its upper bound , , because it defines the largest feasible region as . Furthermore, because this model aims to find the minimum of all minimum objective values, the constraints of the first level can be inserted into the second level to transform the twolevel mathematical program of Model (12a) into the following classical onelevel mathematical program:
This model is a classical linear program that can be solved easily using LP solvers. Because all , have been set to the upper bound of their αcut in this model, is guaranteed, as required by Definition (8) based on Zadeh’s extension principle.
For Model (12b), however, because the first level is a maximisation problem that is not consistent with the minimisation operation of the second level, the constraints of the first level cannot be directly inserted into the second level. To address this situation, the dual of the secondlevel problem is formulated as consistent with the maximisation operation of the first level. According to the duality theorem of LP, it is well known that the primal and dual models share the same objective value [42]. Consequently, Model (12b) can be reformulated as follows:where , and , are the sets of dual variables defined for the first to sixth sets of constraints in Model (12b), respectively.
Similar to the lower bound case, in Model (14), the upper bound of the objective value can be derived by setting , , which yields the largest feasible region. Furthermore, because both the first and second levels perform the same maximisation operation, their constraints can be combined to form a classical mathematical program. Consequently, Model (14) can be rewritten as which is a linearly constrained nonlinear programming (NLP) model that can be solved using many efficient and effective solution methods, including the generalised reduced gradient (GRG) method and successive quadratic programming (SQP) [48].
The solution of Models (13) and (15) yields the crisp interval . By enumerating different values of α, the approximated membership function of , or the fuzzy minimum total cost of the APP model, can be derived by numerically identifying the set of intervals . Meanwhile, the corresponding optimal APP is also completed.
Note that the minimum total cost obtained by adopting the proposed approach is a fuzzy performance measure. The significance of such a fuzzy performance measure is that it conserves the fuzziness of the input parameters; thus, it is able to more accurately represent a system with nonstochastic factors. Consequently, the proposed approach can obtain more realistic performance measures when certain parameters in the APP are fuzzy. These features are important and useful to the decision maker in practice.
4. Numerical Application
The following example inspired by previous studies illustrates the validity of the proposed model and solution procedure in this paper [30]. The example is a sixperiod (monthly base) APP problem involving two types of products with several parameters characterised by triangular fuzzy numbers as follows: the unit production cost excluding labour cost, overtime labour cost, labour cost, inventory carrying cost, cost to hire one worker, cost to lay off one worker, unit backorder cost, conversion factor in hours of labour, forecasted demands, and maximum workforce, as shown in Tables 3, 4, 5, 6, 7, and 8. The assumed trend for future demand is considered optimistic in periods 1 and 2, pessimistic in periods 5 and 6, and neutral in periods 3 and 4. The example is further characterised by the following.(i)The regular labour force works eight hours a day; that is, .(ii)The initial workforce consists of 100 workers for product 1 and 150 workers for product 2; that is, and .(iii)Overtime production is limited to no more than 30% of regulartime production; that is, .(iv)The initial inventories of products 1 and 2 are 0; that is, and .(v)Backorders must not be carried over for more than one period.



