Computationally Efficient Problem Reformulations for Capacitated Lot Sizing Problem

Abstract

In this article, we propose novel reformulations for capacitated lot sizing problem. These reformulations are the result of reducing the number of variables (by eliminating the backorder variable) or increasing the number of constraints (time capacity constraints) in the standard problem formulation. These reformulations are expected to reduce the computational time complexity of the problem. Their computational efficiency is evaluated later in this article through numerical analysis on randomly generated problems.

Share and Cite:

Sharma, R. , Sinha, P. and Verma, M. (2018) Computationally Efficient Problem Reformulations for Capacitated Lot Sizing Problem. American Journal of Operations Research, 8, 312-322. doi: 10.4236/ajor.2018.84018.

1. Introduction

Lot sizing problem aims to optimally utilize the available production resources while meeting the demand targets. It is classified as medium-term planning in production planning taxonomy. Lot sizing problem formulation depends upon the layout and the operating constraints in the production system. In the manufacturing industry, we come across many types of production systems. These production systems further give rise to different types of lot sizing problems (with different constraints and operating conditions) and their solution methodologies. Hence there is a rich literature on lot sizing problem and their solution methods. In this article, we restrict our discussion to general dynamic multi-level capacitated lot sizing problem.

This problem was first proposed by Billington et al. [1] . It addresses the following scenario: a finite planning horizon is given and is divided into discrete time periods. There is a dynamic demand for items which needs to be satisfied for each time period while honoring the production capacity constraints. Problem aims to develop a production plan over all time periods while minimizing the total cost comprising of setup cost, inventory cost, and backordering cost.

A capacitated lot sizing problem is a well-known NP-hard problem. If the capacity constraints of the problem are relaxed, then the problem can be solved in polynomial time [2] . Many model formulations have been proposed to develop an efficient numerical solution. These formulations differ with each other due to variables and constraints used. Formulations given in this article belongs to inventory and lot-size (I & L) formulations category, which are among the most popular in the literature due to their computational efficiency. These formulations use production quantity and inventory levels as the variables. Relaxation of each constraint in the standard I&L model affects the problem structure and hence its numerical complexity. Relaxation of capacity constraints decomposes CLSP into single-level lot sizing problem popularly known as Wagner-Whitin problem [2] . Further, some popular extensions to the standard problem have been suggested in the literature, to address certain practical issues. For example, Dillenberger et al. [3] have extended the problem to incorporate setup carryover. A setup cost and setup time are not incurred if the same item is being produced in the next time period, carried forward by the previous period. Our formulation incorporates binary variable for setup to address these issues. Binary setup variable to address the carryover of setup to the next period was earlier used by Hasse [4] , and Surie and Stadtler [5] .

We next discuss certain problem reformulations from the literature which are computationally efficient. CLSP can be formulated to assign each production quantity to a demand in a specific time period while minimizing the production cost. Shortest route formulation was proposed by Eppen and Martin [6] for a single level case. It was later extended by Tempelemier and Helber [7] for multi-level CLSP. Stadtler proposed an improvement in this formulation ( [8] [9] ), which decreases the number of non-negative coefficients. Rosling [10] introduced a formulation based on Plant location problem analogy. This formulation was further extended by Maes et al. [11] for the capacitated case of a lot-sizing problem. Capacity constraints were included in the original formulation for this purpose. Equivalence of Shortest route and SPL formulation in terms of objective function was shown by Denizel et al. [12] .

Apart from reformulation, additional inequalities can be added to the problem formulation to tighten the bound while reducing the search space. Important researches in this category are discussed next. Barany et al. [13] included lot sizing and inventory variables for the single level uncapacitated lot-sizing problem. Additional valid constraints are included in the formulation to tighten the convex bound of the uncapacitated lot-sizing problem. Pochet and Wolsey [14] ; and Clark and Armentano [15] extended the work of Barany [13] for the multi-level case. Miller et al. [16] have proposed additional valid inequalities for the capacitated case of a lot-sizing problem. Further Surie and Stadtler [5] have proposed valid inequalities for multi-level capacitated lot sizing problem with set up carry over. Setup carryover constraints are redefined to achieve a computational advantage in this formulation.

Research in this article is based on the appropriate reformulation of the standard capacitated lot sizing problem. We state the standard problem formulation and then derive three reformulations of the problem by eliminating the backordering variable or/and adding two capacity constraints. Efficacies of these formulations in terms of reduced computational complexity are demonstrated through numerical analysis of random problems on GAMS.

2. Research Methodology

As stated in the previous section, we intend to evaluate the improvement in computational efficiency of the model when the number of decision variables are decreased, or constraints are added to tighten the bound of the solution space. Model A1 is the reference standard model, which is tinkered to develop model A2, A3, and A4 accordingly. In model A2 (proposed later), we have eliminated the backordering variable; hence it is expected to be computationally efficient when compared with model A1. Similarly, we have added two extra constraints (Equation (27), Equation (28)) in our standard model (A1) which is referred to as model A3. Further, we eliminate backordering variable while adding two constraints in model A1 and refer it to model A4. Hence model A4 is expected to perform best among all. Efficacy of each model is evaluated by performing paired t-test of the computational time of random problem instances on A1, A2, A3, and A4. Branch and bound method is used in GAMS for solving random problem optimally by these models. Finally it is concluded in section 6 that most computationally efficient formulation should be used solving capacitated lot sizing problem.

3. Problem Formulation/Reformulation

Table 1. Notations used in the model.

3.1. Model A1

Minimize Z = i = 1 I t = 1 T [ C P i t * X P i t + C S i t * Y S i t + C I N V i t * X I N V i t + C B O i t * X B O i t ] (1)

Subject to:

X P i t + X I N V i , t 1 + X B O i t = D i t + X I N V i t + X B O i , t 1 i I , t T (2)

i = 1 I ( P T i t * X P i t + S T i * Y S i t ) C A P T t t T (3)

i = 1 I ( P T i t * D i t + S T i * Y S i t ) C A P T t t T (4)

X P i t C A P i t Y S i t i I , t T (5)

t = 1 T X P i t t = 1 T D i t i I (6)

X I N V i 0 = 0 i I (7)

X I N V i T = 0 i I (8)

X B O i 0 = 0 i I (9)

X B O i T = 0 i I (10)

Y S i t { 0 , 1 } i I , t T (11)

X I N V i t , X P i t , X B O i t 0 i I , t T (12)

3.2. Model A2

Minimize Z = i = 1 I t = 1 T [ C P i t * X P i t + C S i t * Y S i t + C I N V i t * X I N V i t ] + t 1 = 1 T i = 1 I C B O i t 1 * ( t = 1 t 1 D i t + X I N V i t 1 t = 1 t 1 X P i t ) (13)

Subject to:

t = 1 t 1 D i t + X I N V i t 1 t = 1 t 1 X P i t 0 i I , t 1 T (14)

i = 1 I ( P T i t * X P i t + S T i * Y S i t ) C A P T t t T (15)

X P i t C A P i t Y S i t i I , t T (16)

t = 1 T X P i t t = 1 T D i t i I (17)

X I N V i 0 = 0 i I (18)

X I N V i T = 0 i I (19)

Y S i t { 0 , 1 } i I , t T (20)

X I N V i t , X P i t 0 i I , t T (21)

3.3. Model A3

Minimize Z = i = 1 I t = 1 T [ C P i t * X P i t + C S i t * Y S i t + C I N V i t * X I N V i t + C B O i t * X B O i t ] (22)

Subject to:

X P i t + X I N V i , t 1 + X B O i t = D i t + X I N V i t + X B O i , t 1 i I , t T (23)

i = 1 I ( P T i t * X P i t + S T i * Y S i t ) C A P T t t T (24)

i = 1 I ( P T i t * D i t + S T i * Y S i t ) C A P T t t T (25)

X P i t C A P i t Y S i t i I , t T (26)

i = 1 I t = 1 T ( P T i t * D i t + S T i * Y S i t ) t = 1 T C A P T t (27)

i = 1 I t = 1 T ( P T i t * X P i t + S T i * Y S i t ) t = 1 T C A P T t (28)

t = 1 T X P i t t = 1 T D i t i I (29)

X I N V i 0 = 0 i I (30)

X I N V i T = 0 i I (31)

X B O i 0 = 0 i I (32)

X B O i T = 0 i I (33)

Y S i t { 0 , 1 } i I , t T (34)

X I N V i t , X P i t , X B O i t 0 i I , t T (35)

3.4. Model A4

Minimize Z = i = 1 I t = 1 T [ C P i t * X P i t + C S i t * Y S i t + C I N V i t * X I N V i t ] + t 1 = 1 T i = 1 I C B O i t 1 * ( t = 1 t 1 D i t + X I N V i t 1 t = 1 t 1 X P i t ) (36)

Subject to:

t = 1 t 1 D i t + X I N V i t 1 t = 1 t 1 X P i t 0 i I , t 1 T (37)

i = 1 I ( P T i t * X P i t + S T i * Y S i t ) C A P T t t T (38)

X P i t C A P i t Y S i t i I , t T (39)

i = 1 I t = 1 T ( P T i t * D i t + S T i * Y S i t ) t = 1 T C A P T t (40)

i = 1 I t = 1 T ( P T i t * X P i t + S T i * Y S i t ) t = 1 T C A P T t (41)

t = 1 T X P i t t = 1 T D i t i I (42)

X I N V i 0 = 0 i I (43)

X I N V i T = 0 i I (44)

Y S i t { 0 , 1 } i I , t T (45)

X I N V i t , X P i t 0 i I , t T (46)

Problem notations are tabulated in Table 1. Equation (1), Equation (13), Equation (22), and Equation (36) minimize the total production cost of the system in model A1, A2, A3, and A4 respectively. Equation (2), Equation (14), Equation (23), and Equation (37) are the state equations as they ensure that the total quantities in a particular time period is a function of total quantities carried forward from the preceding time period while satisfying demand. It must be noted that in model A2 and A4, the backorder variable has been eliminated from the model by appropriate substitution (backorder variable X B O i t , has been eliminated by substituting it in terms of other decision variables from Equation (2). Its value has been substituted in the objective function Equation (13) in model A2 and objective function Equation (36) in model A4). These changes are reflected in Equation (23) and Equation (37). Reduction in number of variables improves the time complexity of the model. Equation (3), Equation (4), Equation (15), Equation (24), Equation (25), Equation (38), Equation (40), and Equation (41) are the time capacity constraints. They ensure that production of the items, and demand satisfied (through production) in a particular time period does not violate the production time available in that time period. Similarly, Equation (5), Equation (16), Equation (26), and Equation (39) are the production capacity constraints. These constraints ensures that capacity constraints in terms of production resources are not violated in any time period. Equation (27), Equation (28), Equation (40), and Equation (41) are the additional capacity constraints added in model A3, and A4 for tightening the bound and subsequently achieving computational advantage as discussed earlier. These constraints are derived from constraint Equation (24) and Equation (25) by extending the time capacity constraints over the entire time horizon T. It must be noted that Equation (24) and Equation (25) ensures the time capacity constraints are honored only for individual time period. Equation (6), Equation (17), Equation (29), and Equation (42) ensures that demand is satisfied in each period. Equations (7)-(10), Equation (18), Equation (19), Equations (30)-(33), Equation (43), Equation (44) sets the initial and final conditions (boundary conditions) over the production horizon. Equation (11), Equation (12), Equation (20), Equation (21), Equation (34), Equation (35), Equation (45), Equation (46) are the binary and non-negativity constraints for decision variables.

4. Numerical Experiments

40 problems each of size 5 × 5 and 6 × 6 are randomly generated in GAMS. 5 × 5, 6 × 6 problems denotes the lot sizing problem to find an optimum production plan of 5 items over 5 time periods, and 6 items over 6 time periods respectively. Only feasible problems are retained for data analysis (6 × 6―29 problems, 5 × 5―31 problems). Value of constants in these problems is randomly generated according to normal distribution (Table 2) and uniform distribution (Table 3).

Table 2. Random data generation (Normal distribution).

Table 3. Random data generation (Uniform Distribution).

5. Data Analysis

All the problems are implemented in GAMS. Solution to these sample problems are tabulated in Appendix. According to t-test performed on data, problem A3 is computationally efficient to problem A1 with a statistical significance of 0.009317 (p-value). Similarly A4 is better than A2 with a statistical significance of 0.003071 (p-value). Model A2 is computationally more efficient than model A1 with a statistical significance of 0.000695 (p-value). Model A4 is computationally more efficient than model A3 with a statistical significance of 0.00473 (p-value).

6. Conclusion

In this article, we have demonstrated the effect of reducing the number of variables, increasing the number of constraints on the computational time of lot sizing problem through 4 models. We infer from our data analysis that model A4 is the most computationally efficient model, and hence is recommended to be used for solving capacitated lot sizing problem.

Appendix: Data Analysis Results

(a) (b)

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Billington, P.J., McClain, J.O. and Joseph Thomas, L. (1983) Mathematical Programming Approaches to Capacity-Constrained MRP Systems: Review, Formulation and Problem Reduction. Management Science, 29, 1126-1141. https://doi.org/10.1287/mnsc.29.10.1126
[2] Wagner, H.M. and Whitin, T.M. (1958) Dynamic Version of the Economic Lot Size Model. Management Science, 5, 89-96. https://doi.org/10.1287/mnsc.5.1.89
[3] Dillenberger, C., Escudero, L.F., Wollensak, A. and Zhang, W. (1993) On Solving a Large-Scale Resource Allocation Problem in Production Planning. Operations Research in Production Planning and Control, Springer, Berlin, 105-119. https://doi.org/10.1007/978-3-642-78063-9_7
[4] Haase, K. (2012) Lotsizing and Scheduling for Production Planning. Vol. 408, Springer Science & Business Media.
https://books.google.co.in/books?hl=en&lr=&id=QVL-CAAAQBAJ&oi=fnd&pg=PA1&dq=Lot sizing+and+scheduling+for+production+planning&ots=HSTSWCFLRT&sig=nDlJCwQpjGkD0DPW ZUGI2RMa3Kg#v=onepage&q=Lotsizing%20and%20scheduling%20for%20production%20pl anning&f=false
[5] Suerie, C. and Stadtler, H. (2003) The Capacitated Lot-Sizing Problem with Linked Lot Sizes. Management Science, 49, 1039-1054. https://doi.org/10.1287/mnsc.49.8.1039.16406
[6] Eppen, G.D. and Kipp Martin, R. (1987) Solving Multi-Item Capacitated Lot-Sizing Problems Using Variable Redefinition. Operations Research, 35, 832-848.
https://doi.org/10.1287/opre.35.6.832
[7] Tempelmeier, H. and Helber, S. (1994) A Heuristic for Dynamic Multi-Item Multi-Level Capacitated Lotsizing for General Product Structures. European Journal of Operational Research, 75, 296-311. https://doi.org/10.1016/0377-2217(94)90076-0
[8] Stadtler, H. (1996) Mixed Integer Programming Model Formulations for Dynamic Multi-Item Multi-Level Capacitated Lotsizing. European Journal of Operational Research, 94, 561-581.
https://doi.org/10.1016/0377-2217(95)00094-1
[9] Stadtler, H. (1997) Reformulations of the Shortest Route Model for Dynamic Multi-Item Multi-Level Capacitated Lotsizing. Operations-Research-Spektrum, 19, 87-96.
https://doi.org/10.1007/BF01545506
[10] Rosling, K. (1986) Optimal Lot-Sizing for Dynamic Assembly Systems. Multi-Stage Production Planning and Inventory Control, Springer, Berlin, 119-131.
https://doi.org/10.1007/978-3-642-51693-1_7
[11] Maes, J., McClain, J.O. and Van Wassenhove, L.N. (1991) Multilevel Ca-pacitated Lotsizing Complexity and LP-Based Heuristics. European Journal of Operational Research, 53, 131-148.
https://doi.org/10.1016/0377-2217(91)90130-N
[12] Denizel, M., Tevhide Altekin, F., Süral, H. and Stadtler, H. (2008) Equivalence of the LP Relaxations of Two Strong Formulations for the Capacitated Lot-Sizing Problem with Setup Times. OR Spectrum, 30, 773-785. https://doi.org/10.1007/s00291-007-0094-3
[13] Barany, I., Van Roy, T.J. and Wolsey, L.A. (1984) Strong Formulations for Multi-Item Capacitated Lot Sizing. Management Science, 30, 1255-1261.
https://doi.org/10.1287/mnsc.30.10.1255
[14] Pochet, Y. and Wolsey, L.A. (1991) Solving Multi-Item Lot-Sizing Problems Using Strong Cutting Planes. Management Science, 37, 53-67. https://doi.org/10.1287/mnsc.37.1.53
[15] Clark, P.W. and Armentano, L.E. (1997) Replacement of Alfalfa Neutral Detergent Fiber with a Combination of Nonforage Fiber Sources. Journal of Dairy Science, 80, 675-680.
https://doi.org/10.3168/jds.S0022-0302(97)75986-1
[16] Miller, A.J., Nemhauser, G.L. and Savelsbergh, M.W.P. (2000) On the Capacitated Lot-Sizing and Continuous 0-1 Knapsack Polyhedra. European Journal of Operational Research, 125, 298-315.
https://doi.org/10.1016/S0377-2217(99)00461-0

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.