Next Article in Journal
Potential of Jatropha curcas L. as Biodiesel Feedstock in Malaysia: A Concise Review
Next Article in Special Issue
A Process-Based Modeling Method for Describing Production Processes of Ship Block Assembly Planning
Previous Article in Journal
A Novel Bearing Fault Diagnosis Method Based on GL-mRMR-SVM
Previous Article in Special Issue
Introducing Risk Considerations into the Supply Chain Network Design
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Capacitated Lot-Sizing Problem with Sequence-Dependent Setup, Setup Carryover and Setup Crossover

Department of Industrial Engineering, Chosun University, 309 Pilmun-Daero, Dong-Gu, Gwangju 61452, Korea
Processes 2020, 8(7), 785; https://doi.org/10.3390/pr8070785
Submission received: 30 May 2020 / Revised: 3 July 2020 / Accepted: 3 July 2020 / Published: 5 July 2020
(This article belongs to the Special Issue Advances in Sustainable Supply Chains)

Abstract

:
Since setup operations have significant impacts on production environments, the capacitated lot-sizing problem considering arbitrary length of setup times helps to develop flexible and efficient production plans. This study discusses a capacitated lot-sizing problem with sequence-dependent setup, setup carryover and setup crossover. A new mixed integer programming formulation is proposed. The formulation is based on three building blocks: the facility location extended formulation; the setup variables with indices for the starting and the completion time periods; and exponential number of generalized subtour elimination constraints (GSECs). A separation routine is adopted to generate the violated GSECs. Computational experiments show that the proposed formulation outperforms models from the literature.

1. Introduction

1.1. Background

Master production plans of a main production process of a supply chain provide all the participants with the visibility regarding the operation. Production planning is referred to as “lot-sizing” in a discrete production environment [1]. The capacitated lot-sizing problem (CLSP) has been considered in various industries to generate production plans reflecting production restriction [1,2].
Generally, lot-sizing problem assumes that planning horizon is finite and can be divided into multiple time periods, which are terms as “time buckets”. The time buckets are classified into small- and big-time buckets, when the comparison of the relative length of the time period is made with respect to the relative length of the production lot [3]. Models considering small-time buckets allow at most one or two products to be produced in a period and models considering big-time buckets allow more than two products to be produced in a period [1].
Setup operations prepare processing units (e.g., materials and machines) to manufacture production lots. In process industries, the setup time required for replacing production models accounts for a large proportion of capacity utilization. Moreover, variations in the setup time may be highly depending on the product models and their operating sequence. The setup time between similar products is relatively short, allowing engineers to setup multiple times in a day or a shift. However, setup operations requiring more than one day or several days may be necessary between products when the processing material or equipment needs to be changed. For solving the big-time bucket CLSP model, several modeling techniques have been developed to reflect these characteristics of setup time. Sequence dependent setup time model is used to reflect setup time that depends on product order.
Several studies [4,5] have introduced CLSP models, in which the setup state can be preserved across several periods and have used various terms such as “linked lot sizes” and “setup carryover”. In this paper, the concept is called as setup carryover. The setup carryover allows a setup state to be maintained from one period to the following one—namely, a setup carryover model enables a product to be produced in the following buckets without additional setup time. In this concept, production planers in process industries often prefer to produce only a single product without model change during several days or shifts.
In some manufacturing environments, the setup started in period t can be split between the periods t and t + 1 [6,7]. In this paper, the concept is called as “setup-splitting”.
In most continuous manufacturing industries such as process industries, due to the presence of long setup times, manufacturers allow the setup to be carried across more than one period [1], the concept is called as “setup crossover”. Setup crossover enables the incomplete setup operation to cross over the boundaries of time periods. In the case that setup times are longer than a period length, the setup operation may be performed in more than two periods. However, if all the setup times are less than a period, the setup crossover is exactly same as the setup-splitting.

1.2. Literature Review

There are several studies published on big-bucket CLSP with setup times. However, studies on CLSPs with setup carryover are limited. Sox and Gao [4] presented two formulations for the CLSP with setup carryover. Suerie and Stadtler [5] presented a formulation for the CLSP with setup carryover and their computational tests have shown that their formulation is better than the formulation provided by [4]. Further, Menezes et al. [6] proposed CLSPs with setup-splitting considering sequence-dependent and non-triangular setups. In the formulation, the classical subtour elimination constraints, known as the Miller–Tucker–Zemlin (MTZ) constraints, introduced by Dantzig et al. [8] are used. Kopanos et al. [9] introduced a CLSPs with setup carryover and setup-splitting problem with items classified into product families. They considered sequence-dependent setups between products in different families and sequence-independent setups between products in the same family. Mohan et al. [10] extended the formulation of Suerie [5] to address setup-splitting. Camargo et al. [11] presented three formulations for the two-stage lot-sizing and scheduling problem. With computational tests, they showed that the formulation with setup-splitting is the most flexible to incorporate setup-related features. Ramya et al. [7] proposed a formulation for the CLSP with setup carryover and setup-splitting and no backorders or lost sales. Fiorotto et al. [12] compared two formulations proposed by Mohan [10] and by Menezes [6] and proposed symmetry braking constraints. With computational experiments, they showed that the proposed constraints are effective for the formulation proposed by Mohan.
Sung and Maravelias [13], Belo-Filho et al. [14] considered the CLSPs with setup carryover and setup crossover with the setup operations can be split into more than two periods. Sung and Maravelias [13] proposed a big-bucket formulation. Belo-Filho et al. [14] proposed two formulations for the case with backlog. One of the formulations involves a time index disaggregation, defining the setup variables with indices for the starting and the completion time periods of the setup operation.
For the case of sequence dependent setup, two groups of studies have been published. First group of studies consider the problems where at most one lot could be generated per a product for each time bucket. The other group of studies allow multiple lots per a product for each time bucket. The former problems usually consider triangular setup times and costs, while the latter problems usually consider non-triangular setup times or minimum lot sizes. Almada-Lobo et al. [15] presented two models for the CLSP with sequence-dependent and triangular setup times and costs using the MTZ subtour prohibition constraints. Clark et al. [16] formulated a sequencing and lot-sizing model with non-triangular setup times based on the asymmetric traveling salesman problem. Ramya et al. [1] presented two models based on MTZ subtour prohibition constraints for the CLSP with sequence-dependent setup times and setup cost and with setup carryover and setup crossover. In the model, no idle times are allowed during setup operation and production of a product. In the models of [1,15,16] at most one lot per product can be produced in a time period.
Menezes et al. [6] allowed production of multiple lots per period. Clark et al. [17] presented a stronger formulation than Menezes et al. [6] for modeling the production of multiple many a product per period by using a polynomial number of multicommodity-flow-type constraints [18]. Mahdieh et al. [19] presented formulations based on Clark’s formulation for lot-sizing and scheduling with minimum lot-size, non-triangular sequence-dependent setup times and costs with setup carryover and setup-splitting.

1.3. Contribution and Organization

In this paper, a CLSP with sequence-dependent and triangular setup times and with setup carryover and setup crossover is considered. For the problem, a new mathematical formulation is proposed. The formulation is based on the facility location extended formulation [2]. Moreover, it uses the setup variables with indices for the starting and the completion time periods as Belo-Filho’s formulation [14] and exponential number of generalized subtour elimination constraints (GSECs) to prevent subtours in a time bucket.
The mixed integer programming (MIP) formulations are usually compared in terms tightness or compactness [20]. The tightness of an MIP formulation is defined as the difference between the optimal values for the continuous relaxation and the original MIP problem. Tightening an MIP formulation, usually by adding cutting planes or valid inequalities, can reduce the search space that the solver needs to explore [21]. Meanwhile, the compactness refers to the quantity of data that must be processed when solving the problem [20]. Because there are many relaxations that must be solved when solving an MIP problem, the compactness of such a problem always refers to the size of these relaxations. A more compact formulation can speed up the search for the optimal solutions. The tightness and the compactness of the formulation are shown by comparison with one of Mahdieh’s formulation (deferred to as MLOV-SM) [19] and one of Ramya’s formulation (deferred to as MM2) [1].
The remainder of the paper is organized as follows. In Section 2, a formal description of the problem is given. In Section 3, the formulation for the problem is presented. The formulation has exponential numbers of the GSECs. The separation algorithms for the GSECs are explained in Section 4. In addition to comparison with other models, computational experiments and their results are reported in Section 5. Finally, Section 6 concludes this study and suggests some directions for further research.
During the computational experiments, the author found that the setup-splitting is not implemented properly in MLOV-SM. The modified and fixed version of MLOV-SM is reported in Appendix A.

2. Problem Statement

In the following, a formal description of the problem is presented with assumptions for clarifying and simplifying the mathematical model.
A planning horizon consists of finite number of time buckets. For each bucket, an index ranging from 1 to n is assigned; that is, the number of buckets is denoted by n . For each bucket t = 1 , , n , a positive amount of capacity denoted by C t is assigned. In the problem, the production lot sizes and schedules for a set of items, denoted by I , are planned. For each item i in I and bucket t , a non-negative demand quantity denoted by d t i is given. The demand could be satisfied through production before the due date with inventory holding cost h t i for each bucket or after the due date with backlogging cost b t i . Production of one unit of items i in bucket t consumes a positive amount of capacity denoted by p t i , namely the process time per unit. All the demands must be met until the end of the planning horizon. Let T i denote the set of buckets where a positive amount of demands is defined for item i . Considering sequence-dependent setup time and cost, production sequences were decided for each bucket. Production change from item i to j requires setup time τ i j and causes cost c i j . If capacity allows, multiple setup could occur in a bucket. With long setup time or inadequate capacities, setup operation could be performed during consecutive multiple buckets. Production of the last lot could resume in the following bucket without any burden of setup operations. The initial and final setup states are not given. The objective is to minimize the overall cost, which includes backlogging, inventory holding and setup costs.
Under the above descriptions, we can introduce the following two assumptions without loss of generality.
Assumption 1.
No idle time is allowed during setup operations.
Bucket-independent setup times and costs were considered, therefore, any idle during setup operations could be moved to the buckets before setup started or after setup finished without increasing the cost. Therefore, the Assumption 1 is acceptable.
Assumption 2.
In a bucket, at most one lot can be generated for an item.
In [14], the setup variables with indices for the starting and completion time periods of the setup operation were introduced to solve CLSP with setup carryover and setup crossover. According to Assumption 1, for each type of setup (for each item in [14]) with setup time τ , a setup variable can be defined for each ordered bucket pair ( s , t ) satisfying the following constraint:
1 s t n   and u = s + 1 t 1 C u < τ u = s t C u
Proposition 1.
The number of ordered bucket pairs satisfying constraints (1) are limited by 2 n 1 . (Proof: refer to [14]).
For a sequence-dependent setup operation from item i to j , we can consider the setup variable indexed by bucket pairs in the following set:
T i j { ( s , t ) | 1 s t n , u = s + 1 t 1 C u < τ i j u = s t C u }

3. Formulation

The problem formulation is based on three building blocks: the facility location extended formulation [2]; the setup variables with indices for the starting and the completion time buckets [14]; and exponential number of generalized subtour elimination constraints (GSECs) [18].
The following two variables come from the facility location extended formulation: setup state variable, y t i , for item i and bucket t ; and production portion variable, x t u i , for item i and bucket t and bucket u T i (having positive demand d u i ). y t i has value one if setup states for item i in bucket t is on and item i can be produced and 0 otherwise. Moreover, x t u i denotes the proportion of the demand d u i produced in bucket t . For example, in the case ( d 2 i , d 3 i ) = ( 50 , 100 ) , ( y 1 i , y 2 i , y 3 i , x 12 i , x 13 i , x 23 i ) = ( 1 , 1 , 0 , 1 , 0.5 , 0.5 ) means that total demand of d 2 i and half of demand d 3 i are produced in bucket 1 and the remaining half of demand d 3 i is produced in bucket 2.
In the formulation, model change or setup operation is indicated by the following setup-indicating variable: e s t i j , for a pair of items i j I and an ordered pair of buckets ( s , t ) T i j . The variable has a value of one if setup operation from item i to j starts in bucket s and finishes in bucket t and the value is zero otherwise. If s < t , e s t i j indicates setup crossover. When s + 1 < t , e s t i j indicates setup crossover longer than two buckets. By Proposition 1, the number of the variables ( e s t i j ) is bounded by ( 2 n 1 ) × | I | × ( | I | 1 ) . R t i denotes the remaining setup times for item i at the beginning of bucket t in the case of setup crossover. In the formulation, a setup carryover indicating variable, z t i , is used for each item i and bucket t = 1 , , n + 1 . z t i is has value of one if the setup state for item i is maintained at the boundary between bucket t 1 and t and zero otherwise. As mentioned in the problem statement, the initial and finial setup states are not given. Hence, z 1 i ’s decide initial setup state and z n + 1 i ’s decide final setup state.
Lot sequencing with GSECs needs exponential number of constraints, but do not need any additional variables. The process of deriving the constrains is introduced at the end of this section. Table 1 presents the parameters and the variables used in the formulation.
The formulation M1 is as follows:
Min     i I t T i ( s = 1 t 1 u = s t 1 h u i d t i x s t i + s = t + 1 n u = t s 1 b u i d t i x s t i ) + i , j I | i j ( s , t ) T i j c i j e s t i j
Subject to
s = 1   n x s t i = 1 ,   for   all   i I , t T i
x s t i y s i ,   for   all   i I , t T i , s = 1 , , n
u T i p t i d u i x t u i C t y t i ,   for   all   i I , t = 1 , , n
i , j I | i j ( t , u ) T i j τ i j e t u i j + i I u T i p t i d u i x t u i + i I R t i + L t = C t + i I R t + 1 i ,   for   all   t = 1 , , n
i I z t i + i , j I | i j ( s , u ) T i j | s < t , u t e s u i j   1 ,   for   all   t = 1 , , n + 1
z t j + i I | i j ( s , t ) T i j e s t i j = y t j ,   for   all   j I , t = 1 , , n
z t + 1 i + j I | i j ( t , u ) T i j e t u i j = y t i ,   for   all   i I , t = 1 , , n
R t j i I | i j ( s , u ) T i j | s < t , t u τ i j e s u i j   ,   for   all   j I , t = 2 , , n
i , j I | i j ( t , t ) T i j τ i j e t t i j + i I u T i p t i d u i x t u i + L t         C t ( 1 i , j I | i j ( s , u ) T i j | s < t , t < u e s u i j ) ,   for   all   t = 2 , , n 1
i j S | ( t , t ) T i j e t t i j i S \ { k } y t i ,   for   all   S I , k S , t = 1 , , n
e s t i j { 0 , 1 }   for   all   i j I ,   ( s , t ) T i j
z t i { 0 , 1 }   for   all   i I , t = 1 , , n + 1
x s t i 0   for   all   i I , t T i , s T
y t i { 0 , 1 }   for   all   i I , t T
L t 0   for   all   t T
R t i 0   for   all   i I , t = 2 , , n ,   R 1 i = R n + 1 i = 0   for   all   i I .
Here, the objective function (3) minimizes the total cost including backlogging, inventory holding and setup cost. Constraints (4) ensure that each demand is met within the planning horizon. Constraints (5) and (6) ensure that production occurs only in the bucket where a lot is generated. Constraints (7) ensure that capacity consumption during production and setup change plus the idle time equal to the capacity of each bucket, considering the remaining setup times during setup crossover. Constraints (8) ensure that setup state is dedicated to at most one item at the beginning of a bucket. Constraints (9) ensure that a lot comprising one item is generated in a bucket when the setup of the item is carried over or a setup operation to the item ends in the bucket. Constraints (10) ensure that one lot is followed by a setup carryover operation to the next bucket or a setup change operation to another item. Constraints (11) ensure that remaining setup times are allowed only during setup crossover. From Assumption 1, idle time is not allowed during setup operation, this is reflected in constraints (12). First two terms of the left-hand side are added to lift the constraint. Constraints (13) prevent the subtours in each bucket. Constraints (14)–(19) describe the domains of the variables. From Assumption 2, the number of lots is limited to one for each item and bucket and this is reflected in constraints (17). In the formulation, ( 2 n 1 ) × | I | 2 + 2 × | I | binary integer variables, and n × ( i I | T i | + | I | + 1 ) fractional variables are defined. Moreover, total number of constraints except constraints (13) is ( n + 1 ) × i I | T i | + 4 n × | I | + ( n 1 ) × | I | + 3 n 1 .
In the remainder of this section, the process of deriving the constrains (13) is introduced. Constraints (9), (10) and the variables in the constraints can be converted into a network as Figure 1; here, the constraints and variables are converted to nodes and arcs, respectively. Without constraints (13), a sequence-dependent setup usually causes subtours, indicated by the circle with arcs y t i , e t t i j , y t j , and e t t j i in the figure. To discard such a solution, subtour elimination constraints are used.
There are many studies on subtour elimination [18]. To obtain the required results, let us consider a converted network G ˜ t ( V t , A t ) for each bucket t , as depicted in Figure 2. We generated the network using the following procedure:
  • contracting of the arcs relating to y ’s and merging the end nodes into one, denoted by y t i ;
  • merging the tail nodes of arcs crossing the border between bucket t 1 and t into one node, denoted by “source”;
  • and merging the head nodes of arcs crossing the border between bucket t and t + 1 into one node, denoted by “sink”.
Therefore, we obtained the followings:
V t { y t i , i I } { s o u r c e ,   s i n k }
A t { z t i , i I } { z t + 1 i , i I } { e t t i j , i j I , ( t , t ) T i j } { e s t i j , i j I , ( s , t ) T i j | s < t } { e t u i j , i j I , ( t , u ) T i j | t < u } { e s u i j , i j I , ( s , u ) T i j | s < t < u }
Balas [22] introduced the prize-collecting traveling salesman problem to derive a model for scheduling the daily operation of a steel rolling mill. For the problem, Balas proposed the following generalized subtour elimination constraint:
( i , j ) δ + ( S ) e i j ( k , j ) δ + ( k ) e k j for   all   k S , S V \ { s o u r c e , s i n k } ,
where e i j = 1 if setup changes from item i to j , 0 otherwise; δ + ( S ) is the set of outgoing arcs from S . Constraints (22) ensure that the number of selected outgoing arcs from S is at least equal to the number of selected arcs going out from any node k in S . Consequently, subtours can be prevented.
Tacaari [18] rewrote constraints (22) as follows:
( i , j ) A ( S ) e i j i S \ { k } ( i , j ) δ + ( i ) e i j for   all   k S , S V \ { s o u r c e , s i n k } ,
where A ( S ) is the set of arcs between two nodes in S . In G ˜ t ( V t , A t ) , S { y t i , i I } and A ( S ) { e t t i j , i j I | ( t , t ) T i j } ; therefore, by replacing the ( i , j ) δ + ( i ) e i j with y t i (this relation comes from constraints (10)), we can rewrite constraint (23) to fit into our formulation as follows:
i j S | ( t , t ) T i j e t t i j i S \ { k } y t i ,   for   all   S I , k S
Through generating constraints (24) for each bucket t , we can derive constraints (13), the subtour elimination constraint.

4. Branch-and-Cut Algorithm

There are exponential numbers of subsets S I ; therefore, there are exponential numbers of constraints (13) in the formulation. To solve a formulation with exponential number of constraints, we usually use the branch-and-cut algorithm. In the algorithm, most of the constraints are outside of the running formulation and constraints violated by the current solution of linear programming (LP) relaxation are identified dynamically and added into the running formulation. Therefore, we can control the size of the running formulation. The procedure of identifying the violated constraints is referred to as separation (details regarding separation and branch-and-cut algorithm are available in [21]).

4.1. Separation

After obtaining an optimal solution to the LP relaxation (let e ¯ s t i j , y ¯ t i denote the values of the variables e s t i j , y t i , respectively, in the optimal solution) of the current running formulation, for each bucket t = 1 , , n , we use two separation routines to find the violated subtour elimination constraints. One routine is based on the minimum cut, and the other is based on the strongly connected components. In the routine based on the minimum cut, we generate a network G ¯ t ( V ¯ t , A ¯ t ) , illustrated in Figure 3, where V ¯ t { y t i , i I | y ¯ t i > 0 } { s i n k } and A ¯ t A t 1 A t 2   . Here, A t 1 , A t 2 are defined as follows:
A t 1 { e t t i j , i j I | ( t , t ) T i j , e ¯ t t i j > 0 } ,
A t 2 { ( y t i , s i n k ) | i I , y ¯ t i > 0 } .
For each arc e t t i j A t 1 , we assign capacity e ¯ t t i j and for each arc ( y t i , s i n k ) A t 2 , we assign capacity y ¯ t i j | e t t i j A t 1 e ¯ t t i j .
For each vertex y t i V ¯ t \ { s i n k } , the minimum cut between y t i and s i n k yields a node partition ( S , V ¯ t \ S ) , where y t i S and s i n k V ¯ t \ S . If the total capacity of the arcs from S to V ¯ t \ S is less than y ¯ t i , the constraint (24) with k = y t i is violated by the current solution. Thus, we add the constraint into the running formulation. For each time bucket t = 1 , , n , we determine the minimum cut | V ¯ t | 1 times. The time complexity of finding a minimum cut (or maxflow) is O ( | V ¯ t | 2 | A ¯ t | ) by using Goldberg and Tarjan algorithm [23]. Therefore, the overall time complexity is O ( | V ¯ t | 3 | A ¯ t | ) . This routine needs a considerably long computational time but always succeeds in finding the violated constraint.
In the separation routine based on the strongly connected components, we generate a network G ^ t ( V ^ t , A ^ t ) , where V ^ t { y t i , i I | y ¯ t i > 0 } and A ^ t { e t t i j , i j I | ( t , t ) T i j , e ¯ t t i j > 0 } . In the network, we could find some strongly connected components with more than two vertices. In each component S , we select the vertex y t i S with the largest value y ¯ t i ; then, for S and k = y t i , we check if constraint (22) is violated by the current solution. If the constraint is variolated, we add it into current running formulation. The strongly connected components are found in a O ( | V ^ t | + | A ^ t | ) operations by Tarjan’s algorithm [24], and O ( | V ^ t | + | A ^ t | ) operations are needed to check the violation. Therefore, the overall time complexity is O ( | A ^ t | ) . This routine is significantly fast. However, it successfully determines the violated constraint only in the case of integer solution. More information about the separation routine based on the strongly connected components is available in [25].

4.2. Implementation

The formulation and the branch-and-cut procedure are implemented using the Gurobi LP/MIP solver 9.0.2 with Python 3.7. The separation routines are executed by the embedded callback function. To find the minimum cut (or maxflow) and the strongly connected components, the efficient implementations in the Python library, Python-igraph 0.8.2, were used.
The following separation scenario was used in the computational experiments:
  • For integral solutions, the separation routine based on the strongly connected components generates the violated inequalities for each time bucket;
  • For fractional solution founded in i -th branch-and-bound nodes;
    If i 1000 , the separation routine based on the strongly connected components generates the violated inequalities for each time bucket. In the absence of a violated inequality, the separation routine based on the minimum cut is performed for each time bucket.
    If 1000 < i 10 , 000 , the separation routine based on the strongly connected components generates the violated inequalities for each time bucket.
  • If i > 10 , 000 , no separation routine is performed.

5. Computational Experiments

This section describes three sets of computational experiments. The first is for comparison with MLOV-SM in [19], the CLSP with sequence dependent setup, setup carryover and “setup-splitting”. The second is for the analysis of performance of M1 with test instances having long setup times. The third is for comparison with MM2 in [1], the CLSP with sequence dependent setup, setup carryover and “setup crossover”.
Experiments were conducted on a Windows 2010 PC with 16 GB main memory, which was equipped with an Intel(R) Core(TM) i5-9500 @ 3.00 GHz CPU. Default parameter settings for Gurobi LP/MIP optimizer are used, except 600 s of time limit.

5.1. Test Instances

For the first and second tests, three groups of artificial test instances (“T”, “L”, “LA”) were generated. The instances in “T” have short setup times and were designed to be used in the first test. The instances in “L”, “LA” have long setup times and were designed to be used in the second test. They have the following features:
  • Problem dimension: the problem dimension is set by ( m , n ) , where m is the number of items and n is the number of time buckets. Therefore, each item is indexed as i = 1 , , m . The following six problem dimensions were used: (10,10), (10,20), (10,30), (15,10), (15,20), (15,30).
  • Setup time: for i , j = 1 , , m
    For group “T”,
    τ i j = 240 m × | i j | + 30 m ,
    For m = 10 , the minimum and maximum values of τ i j are 44 and 236, respectively.
    For m = 15 , the minimum and maximum values of τ i j are 31 and 239.
    For group “L”,
    τ i j = { 24 × | i j | + 210 if   i 5 = j 5 24 × | i j | + 10 otherwise . ,
    For m = 10 , the minimum and maximum values of τ i j are 34 and 426, respectively. Furthermore, items are divided into three groups: { 1 , 2 , 3 , 4 } ; { 5 , 6 , 7 , 8 , 9 } ; and { 10 } , and most inter-group model changes have a longer setup change times than a bucket capacity of 240. For m = 15 , the minimum and maximum values of τ i j are 34 and 546, respectively. Moreover, items are divided by four groups: { 1 , 2 , 3 , 4 } ; { 5 , 6 , 7 , 8 , 9 } ; { 10 , 11 , 12 , 13 , 14 } ; and { 15 } , and most inter-group model changes have setup change times longer than 240.
    For group “LA”,
    τ i j = { 24 × | i j | + 490 if   i 5 = j 5 24 × | i j | + 10 otherwise .
    For m = 10 , the minimum and maximum values of τ i j are 34 and 706, respectively. Moreover, items are divided into three groups: { 1 , 2 , 3 , 4 } ; { 5 , 6 , 7 , 8 , 9 } ; and { 10 } , and an inter-group model change has setup change time longer than the sum of the capacity of two buckets, 480. For m = 15 , the minimum and maximum values of τ i j are 34 and 826, respectively. Moreover, items are divided into 4 groups: { 1 , 2 , 3 , 4 } ; { 5 , 6 , 7 , 8 , 9 } ; { 10 , 11 , 12 , 13 , 14 } ; and { 15 } , and an inter-group model changes has setup change time longer than 480.
  • Process time per unit:
    p t i = 1 , for all i = 1 , , m , t = 1 , , n .
  • Capacity:
    C t = { 240 , for   t = 1 , , n 1 , 2400 , for   t = n .
    Enough capacity was given to the last bucket to prevent loss of sales.
  • Costs:
    h t i = 3 , for all i = 1 , , m , t = 1 , , n .
    b t i = { 30 , for   all   i = 1 , , m , t = 1 , , n 2 , 300 , for   all   i = 1 , , m , t { n 1 , n } .
    Large backorder cost was given to bucket n 1 to minimize production in the last bucket.
    c i j = τ i j 10 , for i , j = 1 , , m .
  • Demand: for problem dimension ( m , n )
    For each group and problem dimension, initial demand set, named as (1), was generated manually according to the following criteria:
    i I | T i | = 2 n
    | T i | 2 n m , for i I
    All production can be finished before the last bucket.
To increase the reliability of the experiments, for each initial demand set, the due dates were randomly perturbed to generate four additional demand sets, named as (2), (3), (4), (5). The demands of an item having the same due date are merged into one.
For the last test, a test instance introduced in [1] was used. All the above test instances are available on [26].

5.2. First Test: Comparison with the Setup-Splitting Model (MLOV-SM)

As mentioned earlier, MLOV-SM is the model for the CLSP with sequence dependent setup, setup carryover and “setup-splitting”. In addition, it considers situations in which multiple lots can be generated for an item within a time bucket. For direct comparison with our model (M1), MLOV-SM was modified and represented in Appendix A. The modified features are as follows:
  • Setup operation from i to i is not allowed;
  • At most one lot can be generated for an item within a time bucket;
  • The initial setup state is not given but determined by solving the model.
M1 has ( 2 n 1 ) × m 2 + 2 m binary integer variables and n × ( i I | T i | + m + 1 ) fractional variables, whereas MLOV-SM has n m 3 + 2 n m 2 + 2 n m binary variables and 3 n m + 2 n real variables, much more than M1. Because the setup time of the data instance in the group “T” is less than the capacities of the buckets, we can directly compare the lower bounds of the two formulation. We compare the two formulations with respect to the size, tightness of the LP relaxation bound and solution quality. Table 2 shows the tightness of the LP relaxation bound of the formulations. The column headers of the table are given as follows:
  • n : number of buckets;
  • m : number of items;
  • LB0: the objective value of LP relaxation without any additional cuts except GSECs;
  • LB: the objective value of LP relaxation with cuts generated by Gurobi;
  • BEST OBJ: the best integer solution value found within the time limit. Bold means optimal;
  • Gap0: relative gap between LB0 and the best integer solution value (BEST OBJ) at the end.
    G A P 0 = 100 × B E S T   O B J L B 0   B E S T   O B J
  • GAP0 improved: the difference in Gap0 between MLOV-SM and M1.
In Table 2, we can find that M1 give much tight lower bound for the test instances.
Table 3 and Table 4 represent the size and the test result of the formulations. The column headers of table are given as follows:
  • #Vars: number of variables;
  • #Cons: number of constraints (no subtour elimination constraints are counted here);
  • #BVars: number of binary variables;
  • GAP: relative gap between the best lower bound ( Z L B ) and the best integer solution value ( Z I ) at the end;
    Gap = 100 × Z I Z L B   Z I
  • #B node: number of nodes visited in the branch-and-bound tree;
  • Runtime: total running time;
  • GSECs: number of generated subtour elimination constraints.
Table 3 and Table 4 show that M1 has much compact sizes of LP relaxation and gives much better solutions then MLOV-SM.

5.3. Second Test: For the Instance of Long Setup Times (L, LA)

The second test is to see if M1 is stable in solving problems with long setup time. Table 5 and Table 6 show test results for L and LA, respectively. The gap value ‘-’ in Table 5 indicates that no feasible solution found in the time limit.
Table 4, Table 5 and Table 6 show that M1 reliably optimizes the problem instances with ten buckets and give relatively good feasible solution for the problem instances with twenty buckets.

5.4. Third Test: Comparison with the Setup Crossover Model (MM2)

Ramya, et al. [1] introduced the only research result on CLSP with sequence-dependent setup, setup carryover, and setup crossover. In the model, no idle times are allowed during setup operation and setup carryover. MM2 is a formulation introduced in [1] for the problem. But MM2 is too big to solve even problems of moderate size. MM2 has 8 m n 3 + 4 m 2 n 2 + 6 m n 2 + 4 m 2 n + 6 m n + n binary integer variables and more than 150 types of constraints. Ramya, et al. [1] reported that it took more than 26,000 s to solve a test instance with 15 items and 10 buckets and that the solver failed with a file error after 8.5 h due to insufficient space in the hard disk while solving a test instance with 15 items and 15 buckets.
For the comparison with MM2, a modified version of M1 was developed and used for the test. The modified formulation, referred to as M2, needs two types of new variables described in Table 7.
With the variables, M2 is defined as follows:
Min (3) subject to (4)–(7), (9) and (11)–(19) and
i I ( z t i + z ¯ t i ) + i , j I | i j ( s , u ) T i j | s < t , u t e s u i j 1 ,   for   all   t = 1 , , n + 1
z t + 1 i + w t i = y t i ,   for   all   i I , t = 1 , , n
z ¯ t + 1 i + j I | i j ( t , u ) T i j e t u i j = w t i + z ¯ t i ,   for   all   i I , t = 1 , , n
L t C t i I ( w t i + z ¯ t i ) ,   for   all   t = 1 , , n
z ¯ t i { 0 , 1 }   for   all   i I , t = 1 , , n + 1
y t i { 0 , 1 }   for   all   i I , t T
z ¯ 1 0 = 1 .
Here, 0 in z ¯ 1 0 means the initial setup item. M2 has ( 2 n 1 ) × m 2 + 2 n m + 3 m binary integer variables, whereas MM2 has 8 m n 3 + 4 m 2 n 2 + 6 m n 2 + 4 m 2 n + 6 m n + n binary variables. We can say that M2 is a much compact formulation than MM2.
Ramya et al. [1] used six problem instances for the test, but just one instance with 10 items and 15 buckets was reported, so that we can use only that. They reported that MM2 found an optimal solution of the problem in 425 s. M2 was tested with the instance and get the result reported in Table 8. M2 solved the problem just in 2 s to the optimality. In addition, the optimal solution provided the same production plan as [1], except when idle time occurred. Even considering the performance differences in the computers used and limited test instance, it is judged that these differences show the superiority of M2.

6. Conclusions

In this paper, the CLSP with sequence dependent setup, setup carryover, and setup crossover was considered. A new mixed integer programming formulation was introduced. The formulation is based on three building blocks: the facility location extended formulation [2]; the setup variables with indices for the starting and the completion time periods [14]; and exponential number of generalized subtour elimination constraints (GSECs) [18]. A modified version of the separation routine of [25] was adopted to generate the violated GSECs.
Three groups of artificial test instances and an instance from [1] were used for computational experiments. The computational results showed that the proposed formulation outperformed the models from the literature, but the formulation has still large number of binary setup variables, O ( n × | I | 2 ) , to solve the problems with more than twenty buckets in 10 min.
To overcome this drawback, studies to reduce the number of setup variables are needed. These studies could involve using the column generation technique or heuristics based on variable fixing.

Funding

This study was supported by research fund from Chosun University, 2017.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Modified MLOV-SM Model

The formulation for the modified MLOV-SM model [19] is follows:
Min     i , j I t = 1   n c i j y i j t + i I t = 1   n h t i I i t + i I t = 1   n b t i B i t
Subject to
I j t 1 B j t 1 + x j i I j t + B j t = d t j ,   for   all   j I , t = 1 , , n
i I p t i x i t + i , j I τ i j y i j t + S t 1 S t + s l k t = C t ,   for   all   t = 1 , , n
x j t C t p t i × z j t ,   for   all   j I , t = 1 , , n
i I α i t = 1 ,   for   all   t = 1 , , n + 1
α i t + j I ( y j i t O L S j i t ) = z i t ,   for   all   i I , t = 1 , , n
α i t + 1 + j I y i j t j I O L S j i t = z i t ,   for   all   i I , t = 1 , , n
O L S i j t + a i j t k y i j t ,   for   all   i , j , k I , t = 1 , , n
α k t + i I a i k t k = z i t ,   for   all   k I , t = 1 , , n
α i t + j I a j i t k j I a i j t k ,   for   all   i k I , t = 1 , , n
S t i , j I τ i j O L S i j t ,   for   all   t = 1 , , n
j I O L S j i t α i , t + 1 ,   for   all   i I , t = 1 , , n
O L S i j t y i j t ,   for   all   i , j I , t = 1 , , n
x j t , I j t , B j t 0 ,   for   all   j I , t = 1 , , n
S t , s l k t 0 ,   for   all   t = 1 , , n
z i t , y i j t , O L S i j t , a i j t k { 0 , 1 } ,   for   all   i , j , k I , t = 1 , , n
α i t { 0 , 1 } ,   for   all   i I , t = 1 , , n + 1
I i n , B i n ,   I i 0 , B i 0 = 0 ,   for   all   i I
a k j t k = 0 ,   for   all   j , k I , t = 1 , , n
S 0 = 0
y j j t = 0 ,   for   all   j I , t = 1 , , n

References

  1. Ramya, R.; Rajendran, C.; Ziegler, H.; Mohapatra, S.; Ganesh, K. Capacitated Lot Sizing Problems in Process. Industries; Springer: Cham, Switzerland, 2019. [Google Scholar]
  2. Pochet, Y.; Wolsey, L.A. Production Planning by Mixed Integer Programming; Wiley: New York, NY, USA, 2006. [Google Scholar]
  3. Salomon, M.; Kroon, L.G.; Kuik, R.; Wassenhove, L.N. Some extensions of the discrete lot sizing and scheduling problem. Manag. Sci. 1991, 37, 801–812. [Google Scholar] [CrossRef] [Green Version]
  4. Sox, S.R.; Gao, Y. The capacitated lot sizing with setup carry-over. IIE Trans. 1999, 31, 173–181. [Google Scholar] [CrossRef]
  5. Suerie, C.; Stadtler, H. The capacitated lot-sizing problem with linked lot sizes. Manag. Sci. 2003, 49, 1039–1054. [Google Scholar] [CrossRef]
  6. Menezes, A.A.; Clark, A.; Almada-Lobo, B. Capacitated lot-sizing and scheduling with sequence-dependent, period-overlapping and non-triangular setups. J. Sched. 2011, 14, 209–219. [Google Scholar] [CrossRef] [Green Version]
  7. Ramya, R.; Rajendran, C.; Ziegler, H. Capacitated lot-sizing problem with production carry-over and set-up splitting: Mathematical models. Int. J. Prod. Res. 2016, 54, 2332–2344. [Google Scholar] [CrossRef]
  8. Dantzig, G.; Fulkerson, R.; Johnson, S. Solution of a large-scale traveling-salesman problem. Oper. Res. 1954, 2, 393–410. [Google Scholar] [CrossRef]
  9. Kopanos, G.M.; Puigjaner, L.; Maravelias, C.T. Production planning and scheduling of parallel continuous processes with product families. Ind. Eng. Chem. Res. 2011, 50, 1369–1378. [Google Scholar] [CrossRef] [Green Version]
  10. Mohan, S.; Gopalakrishnan, M.; Marathe, R.; Rajan, A. A note on modelling the capacitated lot-sizing problem with set-up carryover and set-up splitting. Int. J. Prod. Res. 2012, 50, 5538–5543. [Google Scholar] [CrossRef]
  11. Camargo, V.C.B.; Toledo, F.M.B.; Almada-Lobo, B. Three time-based scale formulations for the tow-stage lot sizing and scheduling in process industries. J. Oper. Res. Soc. 2012, 63, 1613–1630. [Google Scholar] [CrossRef]
  12. Fiorotto, D.J.; Jans, R.; De Araujo, S.A. An analysis of formulations for the capacitated lot sizing problem with setup crossover. Comput. Ind. Eng. 2017, 106, 338–350. [Google Scholar] [CrossRef] [Green Version]
  13. Sung, C.; Maravelias, C.T. A mixed-integer programming formulation for the general capacitated lot-sizing problem. Comput. Chem. Eng. 2008, 32, 244–259. [Google Scholar] [CrossRef]
  14. Belo-Filho, M.A.F.; Almada-Lobo, B.; Toledo, F.M.B. Models for capacitated lot-sizing problem with backlogging, Setup carryover and crossover. J. Oper. Res. Soc. 2014, 65, 1735–1747. [Google Scholar] [CrossRef]
  15. Alamda-Lobo, B.; Klabjan, D.; Carravilla, M.A.; Oliveira, J.F. Single machine multi-product capacitated lot sizing with sequence-dependent setups. Int. J. Prod. Res. 2007, 45, 4873–4894. [Google Scholar] [CrossRef] [Green Version]
  16. Clark, A.; Morabito, R.; Toso, E. Production setup-sequencing and lot-sizing at an animal nutrition plant through ATSP subtour elimination and patching. J. Sched. 2010, 13, 111–121. [Google Scholar] [CrossRef]
  17. Clark, A.; Mahdieh, M.; Rangel, S. Production lot sizing and scheduling with non-triangular sequence-dependent setup times. Int. J. Prod. Res. 2014, 52, 2490–2503. [Google Scholar] [CrossRef] [Green Version]
  18. Taccari, L. Integer programming formulations for the elementary shortest path problem. Eur. J. Oper. Res. 2016, 252, 122–130. [Google Scholar] [CrossRef]
  19. Mahdieh, M.; Clark, A.; Bijari, M. A novel flexible model for lot sizing and scheduling with non-triangular, period overlapping and carryover setups in different machine configurations. Flex. Serv. Manuf. J. 2018, 30, 884–923. [Google Scholar] [CrossRef]
  20. Yang, L.; Zhang, C.; Jian, J.; Meng, K.; Xu, Y.; Dong, Z. A novel projected two-binary-variable formulation for unit commitment in power systems. Appl. Energy 2017, 187, 732–745. [Google Scholar] [CrossRef]
  21. Nemhauser, G.L.; Wolsey, L.A. Integer and Combinatorial Optimization; Wiley: New York, NY, USA, 1988. [Google Scholar]
  22. Balas, E. The prize collecting traveling salesman problem. Networks 1989, 19, 621–636. [Google Scholar] [CrossRef]
  23. Goldberg, A.V.; Tarjan, R.E. A new approach to the maximum-flow problem. J. ACM 1988, 35, 921–940. [Google Scholar] [CrossRef]
  24. Tarjan, R. Depth-first search and linear graph algorithms. SIAM J. Comput. 1972, 1, 146–160. [Google Scholar] [CrossRef]
  25. Drexl, M. A note on the separation of subtour elimination constraints in elementary shortest path problems. Eur. J. Oper. Res. 2013, 229, 595–598. [Google Scholar] [CrossRef]
  26. CLSP Data Set. Available online: http://ie1.chosun.ac.kr/~jhkang/data/CLSP.zip (accessed on 26 June 2020).
Figure 1. Illustration of a network with setup related constrains and variables and a subtour. (10) Green nodes denote constraints and (9) blue nodes denote constraints. Variables corresponding to thick arcs have values of one, and the others have values of zero. Circular setup change between ( i , j ) generates a subtour.
Figure 1. Illustration of a network with setup related constrains and variables and a subtour. (10) Green nodes denote constraints and (9) blue nodes denote constraints. Variables corresponding to thick arcs have values of one, and the others have values of zero. Circular setup change between ( i , j ) generates a subtour.
Processes 08 00785 g001
Figure 2. Illustration of the converted network G ˜ t ( V t , A t ) for bucket t. Variables corresponding to thick arcs have values of one and the others have values of zero. Circular setup change between ( i , j ) generates a subtour.
Figure 2. Illustration of the converted network G ˜ t ( V t , A t ) for bucket t. Variables corresponding to thick arcs have values of one and the others have values of zero. Circular setup change between ( i , j ) generates a subtour.
Processes 08 00785 g002
Figure 3. Illustration of a network G ¯ t ( V ¯ t , A ¯ t ) for bucket t for the use of cut-based separation.
Figure 3. Illustration of a network G ¯ t ( V ¯ t , A ¯ t ) for bucket t for the use of cut-based separation.
Processes 08 00785 g003
Table 1. Notation.
Table 1. Notation.
TypesNotationMeaning
Set
ISet of items, index i ,   j ,   k
TSet of buckets, { 1 , 2 , , n }
T i Set of buckets where a positive amount of demands is defined for item i
T i j Set of ordered pair of buckets where the first means the starting bucket of the setup operation from i to j and the second means the ending bucket of the setup
Parameters
nNumber of buckets, index s , t , u
C t Capacity of bucket t = 1 , , n
d t i Demand of item i in bucket t
p t i Process time per unit of item i in bucket t
h t i Inventory holding cost of item i in bucket t
b t i Backlogging cost of item i in bucket t
τ i j Setup time from item i to j
c i j Setup cost from item i to j
Variables
x t u i Proportion of the demand for item i in bucket u produced in bucket t , for i I , u T i
y t i Variable indicating setup state: 1 if setup states for item i in bucket t is on and item i can be produced, and 0 otherwise.
z t i Variable indicating setup carryover: 1 if the setup state is for item i at the beginning of bucket t , and 0 otherwise
e s t i j Variable indicating setup change: 1 if setup operation from item i to item j starts in bucket s and finishes in bucket t , and 0 otherwise, for i j I ,
R t i Remaining setup times for item i at the beginning of bucket t in the case of setup crossover
L t Idle time in bucket t
Table 2. The objective values of linear programming (LP)-relaxation of analysis of performance of (M1) and setup-splitting model (MLOV-SM).
Table 2. The objective values of linear programming (LP)-relaxation of analysis of performance of (M1) and setup-splitting model (MLOV-SM).
m nSetMLOV-SMM1GAP0 ImprovedBEST OBJ
LB0LBGAP0LB0LBGAP0
1010(1)144285100%3612482881%19%19,431
1010(2)493529100%1655347991%9%17,970
1010(3)253852100%2390343182%18%13,144
1010(4)96810,69796%904110,72262%33%24,102
1010(5)876369989%3174353261%28%8146
1020(1)406409496%2493389877%19%11,070
1020(2)55714,48298%13,397593960%38%33,840
1020(3)1255535990%4555526563%27%12,259
1020(4)1630673788%5830678957%31%13,463
1020(5)933478794%4197477573%21%15,343
1030(1)2606813484%7387850155%29%16,510
1030(2)1102570193%4704590469%23%15,388
1030(3)417423,46693%20,61526,52964%29%57,589
1030(4)12,67847,23388%45,94853,93356%32%104,245
1030(5)1040743093%6153731861%32%15,975
1510(1)233431100%1628305489%11%15,053
1510(2)139586299%4047776776%23%17,099
1510(3)271882100%1144187087%13%8647
1510(4)1157394484%3736399349%35%7305
1510(5)887749100%6350754368%31%20,002
1520(1)334363097%2715389580%18%13,321
1520(2)107266399%1716279788%11%14,863
1520(3)314569799%3415652188%11%28,001
1520(4)552888298%6761921571%27%23,141
1520(5)438426,92692%24,37331,46758%35%57,563
1530(1)957600194%5173666067%27%15,492
1530(2)1384759993%6331847866%27%18,469
1530(3)1276656192%4505707072%20%16,245
1530(4)520904498%553110,63175%22%22,412
1530(5)938553193%4212586568%25%13,060
Table 3. Test results of MLOV-SM for instance group T.
Table 3. Test results of MLOV-SM for instance group T.
m nSet#Vars#Cons#BVarsGap#B NodeRuntime
1010(1)12,53013,63112,2100%11,54454
1010(2)12,53013,63112,2100%36,390255
1010(3)12,53013,63112,2100%26,189154
1010(4)12,53013,63112,2100%10,295100
1010(5)12,53013,63112,2100%554152
1020(1)25,05027,26124,41023%31,564600
1020(2)25,05027,26124,4100%26,580399
1020(3)25,05027,26124,41027%25,209600
1020(4)25,05027,26124,41031%16,111600
1020(5)25,05027,26124,41042%11,777600
1030(1)37,57040,89136,61027%14,949600
1030(2)37,57040,89136,61043%13,386600
1030(3)37,57040,89136,61042%6217600
1030(4)37,57040,89136,61029%10,364600
1030(5)37,57040,89136,61041%11,469600
1510(1)39,03541,43138,56538%6169600
1510(2)39,03541,43138,56512%7393600
1510(3)39,03541,43138,56529%6431600
1510(4)39,03541,43138,56518%10,003600
1510(5)39,03541,43138,56527%6125600
1520(1)78,05582,86177,11576%2429600
1520(2)78,05582,86177,11583%1996600
1520(3)78,05582,86177,11581%1508600
1520(4)78,05582,86177,11533%4051600
1520(5)78,05582,86177,11550%1814600
1530(1)117,075124,291115,66551%2801600
1530(2)117,075124,291115,66576%1778600
1530(3)117,075124,291115,66564%2267600
1530(4)117,075124,291115,66555%1752600
1530(5)117,075124,291115,66563%2194600
Table 4. Test results of M1 for instance group T.
Table 4. Test results of M1 for instance group T.
m nSet#Vars#Cons#BVarsGap#B NodeRuntimeGSECs
1010(1)237063920700%22825807
1010(2)235061720700%10,117261675
1010(3)236062820700%22,125362007
1010(4)236062820700%8217221468
1010(5)235061720700%7327201731
1020(1)5240164742700%58,8131552323
1020(2)5240164742700%8761484003
1020(3)52001605427012%213,4546004400
1020(4)52801689427029%70,0786005317
1020(5)52601668427018%132,3856003867
1030(1)84402984647017%66,8596005603
1030(2)85303077647030%58,1346006837
1030(3)85003046647014%102,7906009691
1030(4)8470301564707%123,3806005100
1030(5)85303077647028%62,7836008550
1510(1)487583445300%45,0972354429
1510(2)487583445300%7003462876
1510(3)487583445300%28,266933358
1510(4)487583445300%38,9341453944
1510(5)487583445300%21,151923326
1520(1)10,4352084933057%29,7706009788
1520(2)10,3952042933065%30,0176009787
1520(3)10,3952042933040%20,8286009864
1520(4)10,375202193300%46,6903058245
1520(5)10,3952042933013%75,1916006768
1530(1)16,305364114,13036%29,30860018,372
1530(2)16,275361014,13030%25,51160011,348
1530(3)16,365370314,13035%24,84160017,170
1530(4)16,335367214,13015%26,17760011,645
1530(5)16,275361014,13036%42,64060016,000
Table 5. Test results of M1 for instance group L.
Table 5. Test results of M1 for instance group L.
m nSet#Vars#Cons#BVarsGap#B NodeRuntimeGSECs
1010(1)215662818660%30,233441052
1010(2)214661718660%22,589281038
1010(3)215662818660%17,661271005
1010(4)216663918660%19,602321024
1010(5)216663918660%25,84531954
1020(1)4836164738660%59,0502662108
1020(2)48161626386625%174,2866002203
1020(3)48161626386612%187,7866002106
1020(4)47761584386632%108,1656002741
1020(5)4816162638666%176,8936001650
1030(1)78362984586633%40,9776003491
1030(2)78362984586630%48,6416002755
1030(3)78362984586643%29,9506004432
1030(4)77462891586636%63,4326003510
1030(5)78062953586652%30,3856004655
1510(1)447682341410%40,8421332120
1510(2)448683441410%24,924882090
1510(3)448683441410%17,389761865
1510(4)446681241410%36,0681242408
1510(5)447682341410%16,542792336
1520(1)97262063864168%15,1876005991
1520(2)96862021864170%19,4656005356
1520(3)96662000864149%28,4316004613
1520(4)97262063864164%22,6796003930
1520(5)96862021864157%26,4966003846
1530(1)15,316364113,14198%89686019291
1530(2)15,316364113,14184426007097
1530(3)15,286361013,14197276015861
1530(4)15,376370313,14180%10,5406004914
1530(5)15,316364113,14185%74736006544
Table 6. Test results of M1 for instance group LA.
Table 6. Test results of M1 for instance group LA.
m nSet#Vars#Cons#BVarsGap#B NodeRuntimeGSECs
1010(1)209462818040%10,21215745
1010(2)208461718040%577913929
1010(3)209462818040%25,70931953
1010(4)210463918040%3980131105
1010(5)210463918040%879717833
1020(1)47741647380431%135,0296001888
1020(2)47541626380423%125,3426001960
1020(3)4754162638040%48,7641651965
1020(4)47141584380430%120,8606001990
1020(5)47541626380429%118,2986001980
1030(1)77742984580454%23,6086003127
1030(2)77742984580451%51,1176002873
1030(3)77742984580464%33,0126003085
1030(4)76842891580464%33,1786003777
1030(5)77442953580455%33,4846003640
1510(1)429482339590%81,4981602164
1510(2)430483439590%55,2731652116
1510(3)430483439590%106,1692492128
1510(4)428481239590%127,9832562905
1510(5)429482339590%220,7373441769
1520(1)95442063845966%35,1276003679
1520(2)95042021845962%27,6436003878
1520(3)94842000845958%35,2916003297
1520(4)95442063845972%21,7496003921
1520(5)95042021845970%34,1656004398
1530(1)15,134364112,95988%98766015509
1530(2)15,134364112,95991%94166005512
1530(3)15,104361012,95980%10,5036005484
1530(4)15,194370312,95990%10,9286005329
1530(5)15,134364112,95992%77556005304
Table 7. New variables for M2.
Table 7. New variables for M2.
w t i Variable indicating cleanup state: 1 if production of item i is completed and model change is prepared in bucket t , and 0 otherwise.
z ¯ t i Variable indicating carryover of completed job: 1 if the setup state is for completed item i at the beginning of bucket t , and 0 otherwise
Table 8. Test results of M2.
Table 8. Test results of M2.
m n#Vars#Cons#BVarsGap#B NodeRuntimeGSECs
1015375290032380%11.5344

Share and Cite

MDPI and ACS Style

Kang, J. Capacitated Lot-Sizing Problem with Sequence-Dependent Setup, Setup Carryover and Setup Crossover. Processes 2020, 8, 785. https://doi.org/10.3390/pr8070785

AMA Style

Kang J. Capacitated Lot-Sizing Problem with Sequence-Dependent Setup, Setup Carryover and Setup Crossover. Processes. 2020; 8(7):785. https://doi.org/10.3390/pr8070785

Chicago/Turabian Style

Kang, Jangha. 2020. "Capacitated Lot-Sizing Problem with Sequence-Dependent Setup, Setup Carryover and Setup Crossover" Processes 8, no. 7: 785. https://doi.org/10.3390/pr8070785

APA Style

Kang, J. (2020). Capacitated Lot-Sizing Problem with Sequence-Dependent Setup, Setup Carryover and Setup Crossover. Processes, 8(7), 785. https://doi.org/10.3390/pr8070785

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop