Abstract

In order to meet the lead time that the customers require, work-in-process inventory (WIPI) is necessary at almost every station in most make-to-order manufacturing. Depending on the station network configuration and lead time at each station, some of the WIPI do not contribute to reducing the manufacturing lead time of the final product at all. Therefore, it is important to identify the optimal set of stations to hold WIPI such that the total inventory holding cost is minimized, while the required due date for the final product is met. The authors have presented a model to determine the optimal position and quantity of WIPI for a given simple bill of material (S-BOM), in which any part in the BOM has only one immediate parent node. In this paper, we extend the previous study to the general BOM (G-BOM) in which parts in the BOM can have more than one immediate parent and present a new solution procedure using genetic algorithm.

1. Introduction

Since 1964, extensive investigations have been conducted on the materials requirement planning (MRP), which ensures materials to be available to meet manufacturing and market demands. Utilizing master production schedule (MPS) and bill of material (BOM) as input data, MRP determines when and how many items to produce or purchase to meet the MPS. MRP has been widely used in most manufacturing industry for generating a detailed production schedule and purchase order for two decades since introduced.

However, due to the rapid development of technologies, wider choice of customers, and lower forecast accuracy, the performance of MRP became unsatisfactory. In 2011, Ptak and Smith [1] developed a new type of MRP, named as demand-driven MRP (DDMRP), which replaces the previous convention of safety stock with strategically replenished positions. As an innovative multiechelon pull methodology, DDMRP can plan inventory of materials, enable a company to build plans more closely to actual market demands, and promote better and quicker decisions and actions in the planning and execution. Traditional MRP focuses on answering for how much inventory to hold, and when to release order, while DDMRP focuses on answering for where to position the work-in-process inventory (WIPI).

The make-to-order (MTO) strategy makes the end product only when the customer places the order. The lead time for MTO is relatively longer compared to the make-to-stock (MTS) for the customer to get the product. In order to shorten the overall response time to the customer in MTO manufacturing, WIPI needs to be held in the station network. However, holding WIPI in some stations does not contribute to shortening the manufacturing lead time of the final product at all, only to increase the total inventory holding cost. So it is important to identify the optimal set of stations to place WIPI in the station network. We named this problem as the strategic inventory positioning (SIP) problem in our previous work (Rim et al. [2]). In the previous work, we address the SIP problem for the simple BOM (S-BOM) in which any part in the BOM has only one immediate parent part. The current study, however, is an extension or generalization of the previous work in that we address the SIP problem for the general BOM (G-BOM) in which some parts in the BOM have more than one immediate parent.

Only the SIP problem with BOM having a small number of parts can obtain the optimal solution by enumeration within a reasonable computation time. For S-BOM with items, SIP problem is to choose a set of zero-one decisions out of alternatives, of which one must be optimal. However, for G-BOM, if one part has parent parts, the number of alternatives needs to be multiplied by !. As more items have multiple parents in the BOM, it is impractical to enumerate all solutions due to the excessively long computation time.

Therefore the objective of this paper is to develop a new mathematical model and a new genetic algorithm to solve the SIP problem for G-BOM, in which some parts in the BOM have more than one immediate parent. Our goal is to find a solution method that can balance the solution quality with the computational feasibility. We will present genetic operators that will limit the search to a specific set of feasible solutions instead of exploring all possible solutions.

This paper is organized as follows: a review of the related literature is presented in Section 2. ASRLT in S-BOM and G-BOM is introduced in Section 3. Section 4 is dedicated to the problem description and model formulation of the new SIP problem in DDMRP. The process of building the genetic algorithm to solve the problem and the computational results using the proposed GA are given in Section 5. Conclusion and future works are given in Section 6.

2. Literature Review

The inventory positioning problem has been studied by many researchers for past several years. In some previous studies, authors solved the problem of determining the optimal inventory position to increase the service fill rate and minimize the total inventory holding cost against a facility or full supply chains. In order to arrive the best service level, Whybark and Yang [3] proposed a controlled simulation experiment to decide where to place inventory. However, we determine the optimal inventory position depending on the replenishment model presented in DDMRP. Our paper has same targets including increasing service fill rate and minimizing total inventory holding cost with those in the previous literatures, but the method of calculating the lead time of the end product that the customer requires and the total inventory holding cost for all stations is different with those in previous research.

Many researchers apply Simpson’s model to solve the inventory positioning problems. Simpson Jr. [4] used the “all or nothing” policy to decide whether to place the inventory or not in some points in serial line system. He indicated that the service time is independent variable, and the total inventory holding cost is the objective function. The safety stock for time period is considered as the average inventory when he calculated the total inventory cost. He defined the safety stock as multiplying standard deviation of the demand for time period by the safety factor.

Graves and Willems [57] extended Simpson Jr. [4] model to assembly, distribution, and spanning tree network structures. They employed the periodic review base-stock-system model assuming no capacity constraints. In their paper, they considered the service time as a decision variable and minimizing the total inventory holding cost as the objective. They calculated the total inventory holding cost just as what Simpson Jr. [4] did but proposed DP algorithm.

Lesnaia [8] applied the Graves and Willems [6] model to the supply chain in a manufacturing firm. The difference is that Lesnaia considered the service time as a stochastic model instead of a deterministic model assumed by Graves and Willems. Lesnaia presented a general-network algorithm to manage the safety stock placement problem and determine the optimal service time for the path.

Magnanti et al. [9] solved the inventory positioning problem at production/assembly stages of components in an acyclic supply chain network structure which is not a spanning tree network structure. The service delivery time and inbound replenishment lead time are considered as decision variables and minimizing the total inventory cost is considered as the objective in all stages of the network. They built an efficient MIP formulation and employed a successive piecewise linear approximation approach.

Kaminsky and Kaya [10] developed an effective heuristic algorithm using linear programming modes to handle problems in supply chain including inventory positioning, short and reliable due-date quotation, and order sequencing. They defined two expected times as decision variables. One is that the safety stock of components stored at facility that are received from upstream facility for the production of type will last; another one is that the inventory of finished goods at facility required for the production of product type will last. The objective function consists of inventory costs, average lead time costs, and tardiness cost in a facility.

Inderfurth [11] and Inderfurth and Minner [12] applied the Simpson model to divergent and convergent systems in the network structure. The paper used a periodic review base-stock control policy and considered minimizing the inventory cost as the objective function and the service time as a decision variable.

Some researchers adopted the postponement concept in determining the inventory position in the supply chain such as Ioannou et al. [13]. This idea was initially introduced by Alderson [14] and then has been studied by Bucklin [15], Lee et al. [16], Davis [17], Feitzinger and Lee [18], and Ernst and Kamrad [19].

Rim et al. [2] stated the total inventory cost as the objective function and the strategic inventory position as decision variable that is different from the service time presented in previous studies. The paper used the new and important concept of ASRLT introduced by Ptak and Smith [1]. A mathematical model is built and GA method is proposed to determine the inventory position in BOM to minimize the total inventory cost while meeting the lead time required by customers. In this paper, we continue using the basic theory bought by our last paper (Rim et al. [2]). However, the BOM structure considered in last paper is simplified, so we generalize the BOM structure in this paper so that the calculation of ASRLT becomes more complicated. The new mathematical model and GA method are proposed in this paper so that the SIP problem in DDMRP is studied completely in aspect of BOM structure.

3. ASRLT in S-BOM and G-BOM

In this chapter, we introduce Actively Synchronized Replenishment (ASR) lead time. ASR lead time is defined as the longest unprotected or unbuffered sequence in the BOM for a particular parent (Ptak and Smith [1]). Figure 1 from Rim et al. [2] shows that the manufacturing lead time of part 101, the end product of the BOM, is 2 days, given that three components of level 2 are all available. The longest path to calculate the cumulative lead time is 26 days from 101 through 201, 301,402, to 501P.

Suppose the five shaded nodes 202, 301, 304P, 401, and 402 of the BOM have stocks. Then the ASRLT of the end product is 2 + 6 = 8 days. If part 203 additionally holds stock, the ASRLT of product 101 becomes 2 + 5 = 7 days. The more the nodes we allocate the stock at, the shorter (or equal) the ASRLT will be, obviously with higher total inventory cost.

In the general BOM (G-BOM), some of parts have multiple immediate parents so that the processing of the part will be divided into several sections to meet the requirement of different parent parts. Since the ASRLT will change depending on which immediate parent part to supply first, the sequence is critical.

For example, part 402 in Figure 2 has two immediate parent parts 301 and 203. The processing of part 402 is divided into 2 sections: 402(a) for 301 and 402(b) for 203. Suppose that 2 units of 402(a) and 6 units of 402(b) are required to produce one unit of 101. Then the total processing times of 402(a) and 402(b) are 1 day (=) and 3 days (=), respectively. Depending on which to produce first, two different BOMs appear, as shown in Figures 3 and 4.

In Figure 3, 402(a) is processed first for 301 and the manufacturing lead time (MLT) of 402(a) is 1. However, 402(b) has to wait processing for 203 until 402(a) is finished, so the MLT of 402(b) is 1 + 3 = 4 days. The ASRLT of 301 and 203 is 4 + 1 + 11 = 16 and 6 + 4 + 11 = 21, respectively. The common parent part of 301 and 203 is 101. (In this example, the common part happens to be the end product.) In order to just compare the ASRLT of the common parent part under these two sequences, suppose 202 has inventory. Then the ASRLT of 101 is days. On the other hand, in Figure 4, 402(b) is processed first for 203, so the MLT of 402(b) and 402(a) is 3 and 4, respectively. So ASRLT of 101 is 26 days. By comparing the ASRLT of 101 using the two processing sequences, we find that processing 402(a) first shortens the MLT of 101.

4. The Strategic Inventory Positioning Problem

In the earlier study (Rim et al. [2]) we presented a mathematical model to determine the locations among the stations in the S-BOM which hold WIP inventory in order to meet the customer lead time, while minimizing the total inventory holding cost. We named this problem as the strategic inventory positioning (SIP) problem. In this section, we consider the SIP problem for the G-BOM, in which some nodes of the BOM may have two immediate parent nodes. To distinguish the two SIP problems, the earlier SIP problem, in which all nodes have only one immediate parent node, is named as SIP-S (SIP in S-BOM) problem. The problem we address in this paper, in which some nodes may have multiple immediate parent nodes, is names as SIP-G (SIP in G-BOM) problem. For convenience, we will use the same notations that we used in our earlier paper and add a few new notations to describe the SIP-G problem as follows:: th part counted from the left in the th level of the BOM;: set of the two immediate parent parts of ;: required quantity of to make its immediate parent , ;: total processing time of to make its all immediate parent parts;: processing time of for parent , ;: ASRLT of ;: required quantity of to make unit end product from the leg that its immediate parent is on, ;: required quantity of to make unit end product;: set of immediate children of ;: average daily usage of ;: annual inventory cost of ;: average inventory quantity of ;: unit price of ;: percentage usage of over ASR lead time;: percentage of red zone base of that red zone safety accounts for;: annual inventory holding cost rate;: average daily usage of the end product;: lead time for the end product requested by the customer (i.e., service time);: yellow zone quantity of ;: green zone quantity of ;: red zone base quantity of ;: red zone safety quantity of ;

4.1. Formulation of SIP Problem

The objective function in (2) is to minimize the total inventory cost for all items in the BOM. Equation (3) defines the annual inventory cost of by unit price multiplied by average inventory quantity and annual inventory holding cost rate:

Ptak and Smith [1] propose that in the demand-driven materials requirement planning (DDMRP), the average inventory quantity can be calculated to be the total quantity of the red zone plus half quantity of the green zone, as (4). The yellow zone is usually set to be equal to the usage over ASR lead time, as (5):

The lead time category of the buffer profile has a direct impact on the size of the green and red zone bases. Table 1 describes the recommended ranges of impact for the green and red zone bases for each lead time category. The minimum order quantity (MOQ) in this paper is not considered so that we can express the red zone base and green zone as (6) and (7), respectively:

The variability category will size the red zone safety portion of the total red zone. Red zone safety is an expression of a percentage of the red zone base. Table 2 shows how red zone safety is sized. Red zone safety can be shown as

In G-BOM, one part may have multiple immediate parent parts and is divided into several sections. We can get the number of requirements of each section to process one end product as (9). Then the number of requirements for one end product is the summation of all sections as (10). Average daily usage of is the average daily usage of end product multiplied by required quantity of for one end product, as in (11):

We know the processing time of a part to produce all its immediate parent part. The processing time for each section can be obtained according to the number of requirements to process one end product of each section portion of all requirements of this part to process one end product, as follows:

4.2. ASRLT for G-BOM

Turning the S-BOM to the G-BOM in this paper leads to more complicated expression of ASRLT. There can be more than one parent for one part in the G-BOM. In other words, at least two parts in the BOM have the common child part. In Figure 2, part 402 has two parent parts, 301 and 203. If there is no inventory holding for part 402, it has to be processed for its parent parts in make-to-order manufacturing. Among multiple parent parts of a common child part, which parent part to feed first becomes a scheduling problem. For parents case, ! different sequences need to be compared. In this paper, to make the logic easy to understand and simplify the expression of ASRLT, we consider the problem where each part has at most two parent parts.

Suppose that and are the two parent parts of ; that is, . We assume that the common part is located in one lower level than that of the two parent parts whichever is lower; that is, .

Sequence 1. is processed first to feed . The ASRLT of is maximum of two branches, as in the following:According to (12), we get

We trace the indirect parent (which is the “ancestors” above the immediate parent in the BOM) of and along the “leg” in the BOM. Suppose that and are the parents of and , respectively. Since and are at the same level , the two parts will have a common immediate parent part. Define .

Sequence 2. is processed first to feed : Now we get the new and . Define .

We choose Sequence 2 if , or Sequence 1, otherwise. We want to find a set which minimizes the total cost, given that the ASR lead time of the end product is no longer than the given service time . Now the SIP-G problem is mathematically formulated as follows:

The formulation given in (17) is a nonlinear programming (NLP) model so that it is hard to solve. Therefore, a heuristic search algorithm is required to solve the model. In this study, a genetic algorithm is developed to solve the model. In the next section, we will present the genetic algorithm to efficiently solve the problem given in (17).

5. GA-Based Solution Procedure

5.1. Modeling

We solved the SIP-S problem using GA algorithm (Rim et al., 2014). In this paper, GA that combines the genetic operators such as selection, crossover, and mutation is also used to find the near-optimal solution for SIP-G problem. For the SIP-G problem, if there are parts in the G-BOM, more parts will exist in the new S-BOM after the BOM is split. So we define a chromosome as an array of binary genes instead of in the previous paper, each of which also represents whether the corresponding part in the G-BOM holds inventory or not. To form the initial population of size , chromosomes are randomly generated until valid chromosomes are collected. A chromosome is a valid one only if two conditions are satisfied. One is that the ASR lead time of the end product is no longer than the lead time requested by the customer, as is in the SIP-S problem; the other is that all common parts in the S-BOM after splitting have to keep the same inventory status. For example, there are 13 parts in the G-BOM shown in Figure 2 and 15 parts in the S-BOMs in Figures 3 and 4. Parts 402(1) and 402(2) must have the same inventory status; and all of part 501P must have the same inventory status.

To evaluate the chromosome, we set the fitness value of the chromosome equal to a large enough number minus the total cost generated from the chromosome. Roulette wheel selection and elitism methods are applied to carry on the selection operation. In order to explore better solution, we use single-point crossover and mutation operations. In order to terminate the process of GA, termination conditions must be given. In this paper, we stop the algorithm after several generations, depending on the number of all solutions.

5.2. Computational Results

In order to demonstrate the application of the proposed methodology and to evaluate its performance, in this section two numerical examples are given.

5.2.1. Small Instances

There are 13 parts in the G-BOM given in Figure 2. Table 3 shows the parameter values. We assume that every part has the same variability percentage set as 50% for all parts. The lead time percentage is 30%, 50%, and 80% for long, medium, and short lead time, respectively. The setting of the variability percentage and lead time percentage depends on the company environment. In much previous GA work, the crossover rate is typically set from 0.7 to 0.9; and the mutation rate is set from 0.05 to 0.15. We set the crossover rate as 0.8; and the mutation rate as 0.1. For small instances, we set the population size as 80. The iteration will stop after 30 generations. As we mentioned, Figure 2 can be split to Figures 3 and 4, they have the same data of , , shown in Table 4. In the S-BOM, denote the number counted from the left of the immediate parent part, just the same as in the previous paper (e.g., , in Figures 3 and 4).

In order to show the effectiveness of the proposed GA to solve G-SIP problem, the result obtained from GA is compared with those from enumeration. We accept the solutions only when 402(1) and 402(2) are in the same inventory status, and two 501P have the same inventory status. The number of all solutions is = 49,152. However, 10 days of the ASRLT of the end product has to be met, so 2,400 feasible solutions are ranked from the minimum to maximum total cost generated from those feasible solutions. The number of feasible solutions increases as the lead time of the end product requested by the customer is longer. When we get the solution from GA, the ranking in the enumeration will be matched in the list of feasible solutions. Here, we run GA for 20 times and get the average ranking. The results are summarized in Table 5. Error is not employed to test the effectiveness of the GA here ( represents the total cost of the best solution obtained from GA; is the minimum cost from enumeration), because the total cost of the feasible solution is not consecutive. Table 6 shows the top 10 minimum total costs obtained from the enumeration.

5.2.2. Medium-Sized Problem

In the previous section, we evaluate the effectiveness of GA to solve the G-SIP problem with small instances compared with the result obtained from enumeration. However, since GA is used to solve the larger problem, enumeration will take excessively long time. In this example, we have 50 parts in G-BOM, which means that there are at least 250 solutions, the exact number of solutions depends on the number of parts that has multiple immediate parents and the number of parents. Table 7 lists the parameters for medium-sized problem in which the number of , , , , and is the same as those in small instances and is set as 16. For the same size problem, the computation time depends on the population size and running generations that we set as the termination condition. We set to stop when the programming runs 200 generations. When is selected as 80, the running time of the process is around 10 seconds in the PC. The computation time increases almost proportionally with the increased population size. When is selected to be 80, 120, and 200, the convergence path graphs are presented in Figure 5, where the total cost converges faster and earlier for larger .

In general, the larger population size () we use in genetic algorithm, the better solution we will obtain since more alternatives are considered in each generation as the result of crossover and mutation. That is why the solution quality is higher as increases for the same number of generations in Figure 5. However, even if keeps increasing, the solution quality will not keep improving as it converges to the optimum so fast. If is set over 200, we will get the better solution with less number of generations; however, the total running time will increase. In the GA, we need to balance the solution quality with the computation time. So we need to take the appropriate ; that is to say, it is not true that the bigger the is, the better it is.

6. Conclusion

In this paper, we present a solution procedure for the strategic inventory positioning problem with general BOM, in which certain part has not only multiple child parts but also multiple parent parts. The paper concentrates on finding the optimal strategic inventory position to minimize the total inventory cost in the G-BOM, while satisfying a given ASR lead time of the end product. The G-BOM makes the calculation of ASRLT more complicated. The G-BOM can be split into S-BOMs as many as the factorial of the number of parents of an item in G-BOM, so computing the ASRLT becomes more complicated.

The practical contribution of this paper is that we provide a systematic solution procedure to determine the optimal location of stocks in the make-to-order manufacturing process in order to minimize the total inventory cost, while meeting the manufacturing lead time of the end product requested by the customer, when certain part has multiple parent parts in the BOM.

For further research related with current study, the limited capacity of the processes needs to be considered since common machines can be used in some stations; and due to the limited capacity of the machine, some processing needs to be in a queue until other processing is finished.

Conflict of Interests

The authors declare that there is no conflict of interest regarding the publication of this paper.