Development of a Tool Cost Optimization Model for Stochastic Demand of Machined Products

Francisco G. Pantoja^{1},
Victor Songmene^{1},
Jean-Pierre Kenné^{1},
Oluwole A. Olufayo^{1},
Michael Ayomoh^{2}

Show more

1. Introduction

Research in the field of cost optimization modelling for stochastic inventory management and control has significantly intensified over the past few decades [1] . The ability to minimize the overall production cost of a machining process while adapting to the ever growing levels of stochasticity in the open market can be considered a major step towards its competitiveness. This capability must, however occur without adversely affecting throughput and production quality amidst multiplicity of inherent market conditions. Machining operations generally constitute a large segment in the manufacturing sector with costing effectiveness serving as a critical success factor. In a machining process, the cutting tool represents an essential element to the entire machining system. This is so because its effect is directly linked to a multiplicity of factors such as: the quality of finished product, cycle time of operation and ultimately the cost of machining amongst others. The significance of the influence of machining parameters on the cutting tool, workpiece and other machining components during an operation is quite crucial to deciding the overall production cost. Some notable components that make up an ideal machining operation are as presented in Figure 1. The optimum behavior of a machining operation can be influenced by any one of these member elements.

Figure 1. An illustration of factors influencing tool performance optimization.

1.1. Predictive Challenges to Manufacturing

Companies that focus on industrial metal works are often faced with predictive challenges even though there are several quantitative techniques that can be used to improve production processes such as forecasting, probability distributions and optimization related methods amongst others [2] . Very often, applicability within this research space can be limited hence, resulting in a significant production planning problem. Some of these problems are traceable to inefficiencies in the administrative and production processes derived from a poor strategic planning.

In most companies, the behavior of product sales is random, since it is difficult to absolutely predict demand to determine sales. Due to this challenge, the focus is often placed on alternate areas where the lack of a sound process design can be mitigated and controlled. Some of these areas include:

• Production: The lack of a future programming and being left with the expectation of the needs of the clients, manufacturing companies are forced to fulfill the demand of their clients by incurring additional time in the productive processes.

• Inventories. To the poor management of inventory levels, such as the case of high inventory level, this causes high costs and prevents the expansion of the business.

• Finances: Given the high costs of inventories, companies limit their growth in the purchase of new equipment to make production processes more efficient.

To mitigate these and other problems caused by poor planning of the production, it is necessary to design a methodology that provides accurate predictions of the sales of industrial metal mechanical products or mitigate alternate costs generated elsewhere in production [3] . Actual modeling of industrial systems is complicated because of the variability of factors, in many cases it is necessary to make complex combinations of several distribution functions to resolve such problems.

1.2. Optimization Modelling in Manufacturing

Existing machining models more often than not, are seen to be analyzing the combination of different cutting materials by focusing on their conditions while conducting analytical, numerical, empirical and artificial intelligence based observations. Robust predictive models are often required to accommodate the complex interaction that exists amongst the trio of a workpiece, the cutting tool and the machine tool [4] [5] [6] [7] .

Modeling a cutting process often constitute integrating a system’s planning process to improve productivity and enhance product quality. However, to have a significant progress in the modelling of machining processes, especially on metallic work pieces, some notable variables capable of improving the quality of output should be taken into cognizance. Some of these include the effect of stress, strain rate and temperature amongst others. From the point of view of demand, when a product with intricate properties such as high strength (e.g. aerospace superalloys) is seasonally or randomly produced, it often becomes difficult for industries machining these materials to adequately estimate the quantity of cutting tools needed to meet the demand. On the other hand, working with high strength materials increases the tendency of tool penalty cost due to short tool life. These effects are central to the increase in the cost of products and a reduction in the overall profit margin. These costs could be associated with storage costs, ordering costs, stock-out or shortage costs.

1.3. Recent Studies on Cost Optimization in Machining

Traditional procurement policies are subjective and often premised on periodic supply decisions. These policies are usually based on simplified and idealistic assumptions, as well as on the expected cost criterion, without considering machining factors as well as the finished goods throughput. Some research works [8] [9] [10] [11] [12] have considered more inclusive views into production cost in manufacturing by considering machining conditions. However, the need for a more generalized and adaptable model based on a stochastic demand is needed to optimize production cost in the manufacturing sector.

The rate of tool depletion was found to be a function of tool life and machining conditions. Finding the optimal machining conditions related to tool life and obtaining the optimal order quantity could stabilize cycle length and ordering cycle. A study into production cost in the milling of titanium alloys by considering machining conditions, was conducted by Conradie, Dimitrov and Oosthuizen [12] . In their paper, they designed a cost model which considered pre-cost of manufacturing and auxiliarycost incurred based on machining parameters. Their approach generated means to estimate the cost of titanium milling based on actual machining and pre-manufacturing conditions. Furthermore, due to tools supply unpredictability and market demand variability, numerous researchers have sought to design models which would be used to maximize the efficiency in tool management to mitigate the adverse costs incurred in shortage conditions [9] [10] [11] .

Li, Sarker and Yi [9] in their research into an optimal stocking policy for machining tools, designed a model for stochastically distributed demand and tool lifespan. Their model considered the economic penalty for tool breakdown and its influence on lifespan. They also investigated these effects under a certainty distributed demand scenario [10] . Kouedeu, Kenne, Dejax, Songmene and Polotski [8] in their research on stochastic optimal control of manufacturing systems under production-dependent failure rates, proposed a stochastic dynamic programming formulation and derived numerical optimum policies to satisfy the challenge. Their research focuses on a part inventory production policy, however it does not take into consideration in its estimation the influence of the machining parameters.

In their research work on cost modelling in milling operations, Parent, Songmene and Kenné [11] used an optimization technique premised on operations research to minimize machining time and evaluate the cost of milling operation hence, demonstrating the practicality of the algorithm in solving cost related problems and as a decision-making aid for a central controller flexible manufacturing cell. Their approach yet only considers an optimization of the machining cost in the overall production cost and ignores the inventory management. This approach can be deployed for analysis of the impact of setup time while computing for the capacity of a system emanating from a new design and in determining the manufacturing price of parts.

The above mentioned works, as most resources from the literature, do not consider the attributes that influences the productive length of a tool-life while solving the general optimization problem. It is however important to recognize that the tool-life productive length has a great impact on tool procurement and management of its inventory systems as well as indirectly influence total operating cost.

1.4. Problem Statement

Based on the reality of market variability, which has got so much impact on production processes, modern day manufacturing firms are faced with a rising level of supply uncertainty and demand variability. These fluctuations are characterized with challenges capable of impacting negatively on the optimum performance level and general sustainability of the manufacturing sector. Industrial firms involved in the production of high-valued parts such as the automotive or aerospace industry amongst others are often posed with these challenges [13] . These firms often seek to identify optimum machining parameters to increase tool life while maintaining a high quality production level. Optimization of machining time as well as tool inventory cost management are key factors for attaining high production efficiency for these high-value parts. A functional tool optimization policy capable of reducing costs and sustaining the overall tool optimum conditions has become a major priority in production systems functionality.

Other factors of significant importance in manufacturing processes include machine unit costs, machine uptime costs, ordering costs and lead-time due to their impact on procurement, tool inventory policy and influences on tool lifespan. These considerations if well managed, can reduce the costs associated with production activities, decrease the levels of inventory of tools with low demand and increase the level of service for products considering their seasonality and cyclicality.

This paper presents a combined tool acquisition and cost optimization model capable of predicting and reducing the degree of uncertainty in demand. With this collective information, this will assist in determining the lot-sizing and cycle-time of cutting tools for a better manipulation of the supply chain system. This paper develops a nonlinear cost optimization model based on an inventory policy for cutting tools at optimum machining parameters for the production of high-value mechanical parts. Finally, the study inclusively offers information to a production manager regarding the life of his cutting tools in a production process and permits the flexibility to adapt a production process to suit demands consumer requirements.

Figure 2 shows a flowchart that depicts the overall production cost optimization process developed in this study. The optimization process carried out using LINGO software utilized a combined modelling effect of cost minimization and optimum quantity computation as presented in Figure 2. This results in a minimized production cost and optimized ordering quantity. However, these are achievable only when the ordering quantity is optimized.

Figure 2. Flowchart of the proposed cost optimization process.

2. Model Formulation

In this section, the model development is introduced with the mathematical formulations, assumptions and constraints. Prior to the model formulation, an illustration of the interaction between main cost components and effect of cost modelling is presented as seen in Figure 3. Herein, the dynamical relationship amongst cost, productivity and basic machining parameters are presented. The optimization process of the model in this research seeks to find the lowest point on the total cost curve amidst diverse production costs. This point marks the beginning of the high efficiency zone in choice of machining parameters where efficient machine and tool life utilization is achieved.

From Figure 3, it could also be seen that the drop in total cost crosses the rising productivity of operations hence, creating an efficient range in production. A subset of this range was found between the lowest cost and highest productivity for the highest machining efficiency. This range is based on a set of machining parameters which indicates the most favorable machining conditions. Beyond the efficient performance zone, the occurrence of severe tool damage creates a proportionate rise in production costs and drop in productivity. On the other hand, below the optimal machine utilization point, the tool utilization is not maximized hence leading to a reduction in tool costs yet resulting in a relatively high machining costs due to poor efficiency.

2.1. Hypothesis

This sub-section presents the assumptions considered in the model development process:

1) Tool orders are made on request and not in predefined cyclic batches.

Figure 3. Schematic illustration of interaction between main cost components and effect of cost modelling (adapted from [12] [14] ).

2) The cutting operations tools are mainly divided into two groups, namely: tools for roughing and finishing operations. However, this research will focus on the use a single cutting operation.

3) Once a cutting tool has reached its life span, it is no longer used to avoid tool breakage and the occurrence of a defective process. This also eliminates the penalty cost for tool break.

4) The life span of a tool is considered based on empirical models with the experimentation of similar conditions.

5) Tool vendors are located nearby resulting in a fixed lead time. Hence, all tool demands are satisfied in real-time without any need for back orders.

6) Labor cost is included within machine usage cost and considered to be a constant.

7) Inventory is considered as the average between the initial and final inventory.

8) The demand of part product is random, however, it presents patterns of seasonality, cyclic tendencies that can be followed through forecasting methods such as triple exponential smoothening as premised in (Holt-Winter’s Method). Some examples of these patterns are festive yearly seasons, economicfactors, recession and inflation.

9) Machine time used is inclusive of both operational, setup and installation time.

10) Price per unit product is constant

11) Ordering costs are constant

12) The holding cost is a constant and includes both the warehouse and preservation costs.

13) Initial purchases are below demand quantity.

Objective Function

The focus of this research is to optimize production by obtaining the minimum total cost of operations supported by the inventory system as a complete manufacturing process. This can be represented by an objective function which consists of the cost optimizations elements, these are; the cost of operations, cost of holding, cost of shortage, annual ordering and raw material costs.

$\begin{array}{c}Total\text{\hspace{0.17em}}cost=operational\text{\hspace{0.17em}}cost+shortage\text{\hspace{0.17em}}cost+holding\text{\hspace{0.17em}}cost\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+raw\text{\hspace{0.17em}}material+ordering\text{\hspace{0.17em}}cost\end{array}$

$\text{TC}\left(\text{\$}\right)={C}_{op}\left(\$\right)+{C}_{sh}\left(\$\right)+{C}_{h}\left(\$\right)+{C}_{B}\left(\$\right)+{C}_{o}(\$)$

*TC = total cost of production in over a production time period i.

The next sub-sections present the details of the elements of the objective function:

2.2. Cost Optimization Elements

2.2.1. Operational Cost

Machining time is vital in determining the factor of tool life. The cutting time during production is defined by the material volume to be cut (Length × Width × Height) and the depths of cut (radial and axial) as presented in Figure 4. A representation of Taylor’s extended tool life model was applied in this research. The machining workpiece in this research is a rectangular prism which requires different operational schedules ranging from rough cuts through smoothened finish cuts. This research would however focus on the rough cut hence a single cutting operation. The decision variables associated with the material removal process include: width of cut ( ${w}_{r}$ and ${w}_{f}$ ), depth of cut ( ${d}_{r}$ and ${d}_{f}$ ) and feeds ( ${v}_{fr}$ and ${v}_{ff}$ ).

To estimate the life span of a tool using Taylor’s empirical model as presented in Equation (2), it is necessary to determine the tool cost. The following equations are a derivation of the original Taylor’s model:

${t}_{0}=\frac{k}{{v}^{\alpha}{f}^{\beta}{d}^{r}{w}^{\epsilon}}$

In application for milling operations, the tool life expression ( ${t}_{0}$ ) was extended by adding the radial width of cut with α, β and ε serving as parametric constants. Equations based on this type of operation could be devised for either roughing or finishing operations i.e. ${t}_{0r}$ or ${t}_{0f}$ . However, this research is focused on the roughing operation as presented in Equation (3) below:

${t}_{0r}=\frac{k}{{v}_{r}^{\alpha}{f}_{r}^{\beta}{d}_{r}^{\gamma}{w}_{r}^{\lambda}}$

The parametric constants defined in this model need to be obtained through experimentation due to the fact that the machining conditions such as the tool material-make, workpiece, cutting fluid, rigidity of the assembly and the system vibration amongst other factors, vary significantly from one machine to the other. These can be determined using a statistical design of experiment (DOE) [11] .

The cost of tools utilization per hour is defined here as the ratio of the current market price to the tool life as presented in Equation (5). The tool life is usually expressed in terms of the operations the tool is subjected to over time:

Figure 4. Schematic illustration representation of the material volume to be machined [11] .

${c}_{hr}=\frac{{c}_{pr}}{{t}_{0r}}$

In this research, the primary focus is the cutting operation. Apart from identifying and analyzing the cutting operations, a special consideration is also given to the respective quantities of cutting tools ( ${Q}_{p}$ ) deployed during the overall production task. A combination of all these results in the models below:

$\text{operationalcost}={C}_{op}={c}_{\left(r,f\right)}\cdot {Q}_{p}$

where,

${c}_{r}=\left({k}_{0}\cdot {t}_{opr}+{c}_{hr}\cdot {t}_{opr}\right)$

${c}_{f}=\left({k}_{0}\cdot {t}_{opf}+{c}_{hf}\cdot {t}_{opf}\right)$

${Q}_{p}={Q}_{r}\left(\text{roughing}\right)+{Q}_{f}\left(\text{finishing}\right)$

where ${Q}_{f}=0$ in roughing only operations.

Thus:

${C}_{op}\left(\$\right)=\left({k}_{0}\cdot {t}_{op}+{c}_{hr}\cdot {t}_{op}\right)\cdot {Q}_{p}$

2.2.2. Holding Cost

The holding cost represents all the costs associated with the storage of the inventory until its depletion. These costs usually include tied-up capital, space, insurance, protection, taxes attributed to storage amongst others. The holding cost can be assessed either continuously or on a period-by-period basis. A common consideration for the holding cost is that an initial inventory level exists at the beginning of every period. Based on this, the inventory cost at the beginning of a period I can be obtained from:

${C}_{h}\left(\text{\$}\right)={c}_{h}\left\{{I}_{i}+{Q}_{s}\right\}$

where ${I}_{i}$ is considered as the average inventory level required to avoid shortage of cutting tools. This is computed based on the summation of the initial and final inventory levels divided by two.

${I}_{i}=\frac{\left({I}_{\text{initial}}+{I}_{\text{final}}\right)}{2}$

Furthermore, $\stackrel{\xaf}{{I}_{i}}$ as presented below, represents the total amount of tools required to meet the demand for a given product.

$\stackrel{\xaf}{{I}_{i}}={Q}_{s}+{I}_{i}$

In addition, considering the likely physical limitations of a warehouse in terms of its holding capacity, the capacity of storage can be represented as:

$\stackrel{\xaf}{{I}_{i}}\le {A}_{i}$

Considering that the average inventory, we can add it directly to the objective function as follows:

${C}_{h}\left(\text{\$}\right)={c}_{h}\cdot \stackrel{\xaf}{{I}_{i}}$ (6)

2.2.3. Shortage Cost

The shortage cost ( ${C}_{sh}$ ) represents the costs incurred during the time production has a depleted inventory. This is the cost incurred between this time and the period of a re-order. It includes costs of not using the machine and equipment involved in production.

${C}_{sh}={c}_{sh}\cdot {F}_{t}$ (7)

${F}_{t}\text{}=D\frac{{y}_{i}}{{Q}_{s}}$

where ${y}_{i}$ service as an indicator for when shortage is present in the equation. It can be activated or switched off during computation and holds two statuses which are 0 or 1. Whereas ${F}_{t}$ stands for the amount of missing tools in system.

2.2.4. Material Cost

The material cost of purchasing the workpiece used in production is included for a more inclusive and comprehensive review of the total costs involved in production. It is a factor of the quantity of material and the cost per unit material purchase.

${C}_{B}\left(\text{\$}\right)={Q}_{B}\cdot {c}_{B}$

2.2.5. Ordering Cost

The ordering cost is described as a function of the unit cost of ordering tools, cycle time between orders and a variableterm. The variable term ${z}_{i}$ is used to consider the alternate conditions which are assumed constant in this present model. Such conditions for variability in unit ordering, shipping, delays and alternate reasons which could influence this cost could be expanded through this variable in future research modifications of this model. However, in this research work, the ordering cost is considered as constant of value 1.

${C}_{O}\left(\$\right)={C}_{f}\frac{D}{Q}{z}_{i}$

2.3. Basic Objective Function for Single Operations

The objective function formed to minimize the production cost that involves the operation cost, shortage cost, holding cost, material cost and ordering cost is shown in the Equation (11). The cost of labor is not included in this function. The focal point of this function lies in the operational costs which is derived from the machining parameters and tool conditions.

$\text{MinTC}\left({Q}_{p}{T}_{c}\right)i={C}_{op}+{C}_{sh}+{C}_{h}+{C}_{B}+{C}_{O}$

$\begin{array}{l}\text{Minimised}\text{\hspace{0.17em}}\text{Total}\text{\hspace{0.17em}}\text{Cost}\left({Q}_{p},{T}_{c}\right)i\\ =\left({k}_{0}{t}_{op}+{c}_{hr}{t}_{op}\right){Q}_{p}+{c}_{sh}{\displaystyle {\sum}_{i=0}^{n}D\frac{{y}_{i}}{{Q}_{s}}}+{\displaystyle {\sum}_{i=0}^{n}{c}_{hi}\stackrel{\xaf}{{I}_{i}}}+{Q}_{B}\cdot {c}_{B}+{C}_{f}\frac{D}{{Q}_{o}}{z}_{i}\end{array}$

2.4. Constraints

Tool restrictions

Some decision variables found in the model are the radial width of cut for the roughing operation ( ${w}_{r}$ ) and the axial depth of cut ( ${d}_{r}$ ). This radial width ( ${w}_{r}$ ) is connected to the number of passes and the width of the volume to machine (W):

${w}_{r}\cdot {N}_{w}=W$

Also, the axial depth of cut ( ${d}_{r}$ ) which is connected to the number of axial passes ( ${N}_{p}$ ) and the height (H) of the volume to machine is a decision variable:

${d}_{r}\cdot {N}_{p}=H$

Limits are given to the width of radial cuts ( ${w}_{r}$ ) by the diameter of the tool ( ${D}_{t}$ ):

${w}_{r}\le {D}_{t}$

The depths of axial cuts ( ${d}_{r}$ ) are also limited by the length of the tool flutes ( ${h}_{f\mathrm{max}}$ ):

${d}_{r}\le {h}_{f\mathrm{max}}$

Machine-tool constraints

The power needed for the operation relies on the material removal rate (M_{RR}) and the specific power rating of the workpiece material as seen in Equation (12). It can thus be inferred that at maximum machine power, the maximum metal removal rate (M_{RR}) can be determined. This rate depends on the width of cut (
${w}_{r}$ ), the feed (
${f}_{r}$ ) to the depth (
${d}_{r}$ ), as shown in the equation below:

${w}_{r}\ast {f}_{r}\ast {d}_{r}\le {M}_{RR}$

The material removal rate ( ${M}_{RR}$ ) can be denoted as:

${M}_{RR}\le \text{energycoefficent}\left(\eta \right)\cdot \frac{\text{maximummachinecuttingpower}}{\text{Unitpowerfortheworkpiecematerial}}$ (12)

In this case,

• The cutting speed rate is constrained by the machine’s maximum spindle speed.

$\frac{1000\cdot {v}_{r}}{\text{\pi}\cdot {D}_{t}}\le {N}_{\mathrm{max}};$

• The achievable power is constrained by the maximum safe operational force. $\frac{60\cdot H{P}_{r}}{{v}_{r}}\le {F}_{\mathrm{max}};$

• The cutting feed rate is constrained by the machine’s maximum feed rate. ${f}_{r}\le {f}_{m\mathrm{max}}$

where ${f}_{m\mathrm{max}}$ is the maximum feed of the machine, ${F}_{\mathrm{max}}$ are the maximum forces obtained from the machine and ${N}_{\mathrm{max}}$ is the maximum spindle speed of the machine.

Inventory and shortage restrictions

With restrictions of inventory as follows:

$\left({I}_{i-1}+{Q}_{p}\right)-{I}_{i}={Q}_{t}i$

${I}_{i}\le {A}_{i}$

Restrictions of shortage:

${Q}_{t}-{Q}_{p}={F}_{t}$

${Q}_{p}-\stackrel{\xaf}{{I}_{i}}\left({y}_{i}\right)\le 0$

Finally, non-negativity is assumed for all the decision variables and binary variables:

${F}_{t},\stackrel{\xaf}{{I}_{i}},{A}_{i}\ge 0$ and ${y}_{i}=0,1$

This project focuses on optimizing the process thought for quantity of tool in an inventory. Therefore the study remains within basic machine limitations or overlooks tool immersion and surface finish as constraints.

3. Solution Methodology and Illustrative Example

The model was implemented in “LINGO” by LINDO® systems. “LINGO” is a commercially available optimization modelling software used for building and solving mathematical optimization models. Its package build-up provides a language tool needed to build and design models as well as include needed solving tools within a compact integrated environment. It could be applied to linear, nonlinear, quadratic, integer and stochastic optimization problems. A seven core computer with 8 gigabytes of memory was used for the LINDO analysis to minimize processing time. Note that for each result the LINGO solver provided, it declared the solution as being only a local optimum, i.e. it could not fully guarantee a global optimum.

3.1. Graphical Illustration

The illustration in below depicts the policy of the model for the annual optimization of production in a manufacturing company whose primary task and operations is focused at metal works. This new tool policy includes the forecast of sales and the annual tool planning for production. The main objectives of this illustration is to show how the model can meet a high level of service in sales and a low inventory level of production. As the proposed model is founded on the principle of the Economic order quantity model (EOQ) [15] , a generic inventory system adapted to include the influence of machining parameters is produced (Figure 5).

From the figure, the tool inventory system shows two different ordering quantity cycles which changes based on the variation of tools inventory depletion. By evaluating an accurate estimate of the production demand and correlating it with the productive tool life for machining tools, it is possible to determine the optimal order quantity ( ${Q}_{t}$ ) and optimal-order cycle ( ${T}_{c}$ ). A steep quantity slope due to reduced tool life or frequent tool damage would reduce the order cycle length. This can be prevented by accurately estimating the reorder point from the lead time in obtaining new tools [16] . A minimum safety threshold is also needed to protect against the impact of back-orders and offer enough time

Figure 5. A generic review inventory system (adapted from Li and Cheng [10] , Shang, Tadikamalla, J. Kirsch and Brown [17] ) To the other sections. It is a result with more explanantion.

to react to depleting stock. Beyond this threshold the risks of reaching shortage is high.

It is worth noting that a machining tool has an expected productive life around its mean value; therefore, there exists more risk of failure cost for prolonged use of machining tool over this value. However, a reduction in running time of these tools will potentially increase the replacement quantity and thus re-ordering cost. Therefore, production inventory estimation is contingent on requirement needs based on trends and controlling the utilization time for the tools could be used to regulate the inventory system without influencing the integrity and quality of the product.

3.2. Numerical Sample Case 1

The example shown below was adapted from the work of Wang, Kim, Katayama and Hsueh [18] . In the study, the cost analysis of face milling amedium carbon steel (150 NHB) block, with a three-flute cemented carbide mill is performed. The tool diameter used is 25 mm. The Taylor’s tool life parameters utilized in this study were also obtained from literature [19] . Table 1 and Table 2 lists all the machine tool, cost and production parameters. The solution represented for this initial example, considers only a fixed demand of 2500 produced parts and determines ideal parameters needed to optimize inventory use and reduce cost. All constraints are based on realistic machine parameters.

Table 1. Sample data utilized in numerical example 1; machine constants.

Table 2. Sample data utilized in numerical example 1; cutting parameters, workpiece dimension and cost values.

Results and Discussion

Solving the developed model using the LINGO® solver in Equation (11) gave the following cost breakdown as presented in Figure 6. From the figure, the main contributors to the total production costs were identified as the ordering costs, machine operations costs and holding costs. Ordering costs are known to vary based on the lead-time and unit cost of ordering a tool. As the lead-time to delivery increased, the ordering costs simultaneously increased. The operational cost on the other hand, is greatly susceptible to changes in tool life. An increment in the amount of tools following a poor performing tool life, was seen to significantly increase the incurred costs. The model also indicates the significant cost incurred from inventory of tools. Tool quantity stored during the production of the product bares an influence on its selling price. However, due to the successful minimization of the model, the separate cost of holding the safety stock as well as the shortage costs has become insignificant in the summation.

This elevated holding cost is due to the short cycle period of 5.5 days for reordering. A set of approximately 65 orders per year was obtained. The average forecast quantity of 7 units per order cycle was obtained from Equation (5). This resulted in an optimal order quantity ${Q}_{0}$ of about 38 units to minimize total costs of production. This can also be seen from Figure 7 and Figure 8. The value obtained in contrary to an ideal EOQ model contains the inclusions of the sum of operational, shortage, holding, and material costs. An optimal ordering

Figure 6. Bar chart result of the influence of various costs on total operating cost.

Figure 7. Schematic representation of the order optimization model with added machining op-erational cost.

Figure 8. Schematic of inventory control system with a fixed demand.

quantity was obtained at the intersection of the ordering, sum of operations and safety holding costs. The results also indicated that a safety stock level of 28 units indicating the least amount in the inventory is needed to prevent a shortage situation. Furthermore, a value for re-ordering of stock to match lead time was deduced to be 52 units.

3.3. Numerical Sample Case 2

Table 3 summarizes the data used in this numerical sample case 2. This example was adapted from the work of Parent, Songmene and Kenné [11] and modified to suit machining setup proposed in this study. A typical end milling operation of an alloyed steel part (4140, 4340), with a four-flute uncoated high speed steel (HSS) end mill is proposed for evaluation. The diameter of the tool is 19.05 mm (0.375”). The chip load per tooth is estimated to be 0.05 mm/tooth and the speed between 30 - 45 m/min with a hardness between 30 - 38 HRC. The milling machine setup and initialization values, cost parameters as well as parameters for production are listed in the Table 1 & Table 3. The solution represented in this initial part, considers only a fixed demand of 5000 produced parts and determines ideal parameters needed to optimize inventory use and reduce cost. All constraints are established on realistic machine parameters.

Results and Discussion

Figure 9 presents results acquired using the LINGO® solver. The following cost breakdown presented in the Figure 10 was obtained solving the developed model

Figure 9. Bar chart result of the influence of various costs on total operating cost.

Table 3. Sample data utilized in numerical example 2; cutting parameters, workpiece dimension and cost values.

Figure 10. Schematic representation of tool quantity per order cycle needed at minimum total cost.

in Equation (11). From the figure, the main contributors to the total production costs were identified as the ordering and machine operations cost. However, this example presents a lower inventory holding costs. This is due to the lower cost of purchasing HSS tools and its influence in the combined optimization function. Due to the successful minimization of the model, the holding costs incurred for safety stock as well as the shortage costs remain insignificant in the overall summation.

For the numerical example, a set of 22 orders per year over an approximate 16 days cycle was obtained. The average forecast quantity of 13 units per order cycle was obtained from Equation (5). This resulted in an optimal order quantity

3.4. Sensitivity Analysis

An in-depth sensitivity analysis of numerical sample case 2 is performed to establish influencing factors relevant to changes in operational costs, determine the influence of demand on inventory optimization and ascertain the key cost drivers of the model using actual industrial stochastic demand data.

3.4.1. Cutting Parameters

From an overview of the factors influencing the total cost of consequence, operational machining costs are selected to assess the effect of cutting factors. The

Figure 11. Schematic representation of the order optimization model with added machining op-erational cost.

Figure 12. Schematic of inventory control system with a fixed demand.

most important variables for minimizing this cost area are the cutting speed and the feed rate. These significantly affect the tool life criterion and could create a rising cost avalanche during production. By fixing other factors involved to optimal levels, a plot of the relationship between speed and feed to total cost is shown in the Figure 13.

From Figure 13, it can be seen that as the speed increases a corresponding increase in machining cost is observed. This is expected based on the increased wear at such cutting parameters, however, with increased tool life a corresponding increase in machining cost was also observed. This situation can be explained by the factors controlling an optimum inventory ordering system. The improvement in the tool life changes the tool usage frequency and in consequence increases the holding costs of tools. Based on these changes and the consideration of lead time for ordering, a new optimal order quantity is calculated which is more than the previous to minimize the overall total cost of production.

(a)(b)

Figure 13. Trend of the machining cost based on cutting speed and tool life (a) 3D graph (b) Contour plot.

Similar observations are seen from the feed in Figure 14. An increase in feed to a maximum value of 0.09 mm/tooth also increases the machining cost. However, this increase and effect on the machining cost is more pronounced than that experienced by an increase in speed. Over the whole range of feed rates, a maximum costs increase of $4000 was seen as compared to $2400 with speed increase. The tool life changes based on changes in feed pose a significant effect of the optimal ordering quantity and thus the total cost. The most important variable to consider in estimating machining costs is therefore the feed rate.

3.4.2. Demand Analysis and Cost Effects

The estimated demand ideally influences total production cost since inclusive components costs do increase (Figure 15). However, a cost intensive phase is identified where optimal ordering quantity and safety stock levels are within

(a)(b)

Figure 14. Trend of the machining cost based on feed and tool life (a) 3D graph (b) Contour plot..

proximity (Figure 16). It can be seen that there are two points in these graphs which should be identified when establishing cost expectations in production. The phase before point 1, indicates a non-profitable area of production with a steep increase in costs over a minimal demand. Within this stage fixed costs such as holding costs and inventory costs form a bulk of the total costs in production. As the demand increases, an increased effect of ordering and machining costs is experienced.

From Figure 16, we can also see an interplay between the safety stock and the optimum ordering quantity for minimized costs, which eventually steady out by increasing the safety stock to obtain optimize production costs. Overall costs from point 2 infer that an increased amount of safety stock is to higher demand is preferable. It is postulated that improved demands, reduces the costs impact of fixed expenses and thus increases profits. At such levels, an optimal safety stock quantity higher than the optimal ordering quantity is needed.

Figure 15. Trend showing the influence of change in demand to total cost.

Figure 16. Trend showing the influence of change in demand tool quantity.

3.4.3. Application of Demand Stochasticity

A sensitivity analysis was carried out to ascertain the key cost drivers of the model using actual industrial stochastic demand data. The application of Holt-Winters (HW) forecasting was used to establish the trend for future demands. The current Holt-Winters method of exponential smoothing displays trend and seasonality and is characterized by three smoothing equations: the smoothing for level, equation for the trend and equation for seasonality [16] [20] . As the seasonal component is variable in unit production demand, the multiplicative method is used [21] . This applies when the size of the seasonal component is proportional to the trend level [22] . Figure 17 shows the results of the triple exponential smoothing HW method on the demand for produced units.

The basic equations for the Holt-Winters additive method are:

Equation for level (Overall smoothing):

${L}_{t}=\alpha \left(\frac{{Y}_{t}}{{S}_{t-s}}\right)+\left(1-\alpha \right)\left({L}_{t-1}+{b}_{t-1}\right)$ (12)

Equation for trend (Trend Smoothing):

Figure 17. Trend showing Holt-Winters demand forecasting model for finished alloy steel parts.

${b}_{t}=\beta \left({L}_{t}-{L}_{t-1}\right)+\left(1-\beta \right){b}_{t-1}$ (13)

Equation for seasonality (seasonal index):

${S}_{t}=\gamma \frac{{Y}_{t}}{{L}_{t}}+\left(1-\gamma \right){S}_{t-s}$ (14)

Forecast for m period equals:

${F}_{t+m}={L}_{t}+{b}_{t}m+{S}_{t-sm}$ (15)

where ${Y}_{t}$ are the observed value, ${L}_{t}$ represent the smoothing of variable in time t, ${b}_{t}$ the trend estimation and ${S}_{t}$ is the estimation of seasonality. The smoothing constants α, β, γ are in the interval [0, 1], m is the number of forecast periods and s stands for the duration of seasonality (e.g. in months or quarters in a year).

Manufacturing industrial data figures utilized to establish the adequacy of the model can be found in Table 3. The sensitivity analysis of these results was conceivable by applying maximum and minimum ranges to the machine and user inputs considering what values may realistically be found in industry and accounting for price fluctuations. Industry indexes were used to collect average values for variables and rates. To identify the most important variables of the model, a One-Factor-at-a-Time (OFAT) sensitivity analysis approach was taken [23] . This would be a step in verifying the results of the cost comparison as well as a start in identifying the variables where variable uncertainty would have the highest cost impact. Sensitivity indices were calculated using an equation:

$\text{Sensitivityindex}={D}_{\mathrm{max}}-{D}_{\mathrm{min}}/{D}_{\mathrm{max}}$ (16)

where ${D}_{\mathrm{min}}$ and ${D}_{\mathrm{max}}$ are the minimum and maximum output costs resulting from changing dependent variables to the model.

The machine data and user inputs selected for sensitivity analysis and the applied maximum and minimum values are presented in Table 4.

The results of the sensitivity analysis for the variables in the model is shown in Table 4. Sensitivity analysis of the cost for the daily production period showed

Table 4. Industrial varied data ranges for SA utilized in industrial application sample.

that Machining uptime has a significant influence on the total cost. An increase in cost of approximately ±10% per additional 1000 units demand was also identified as significant. Machine tool life and cutting feed rate have been identified as the most influential cutting variables to total costs. The cutting speed rate had marginal effect on both costs and tool life. Other costs variable such as shortage costs per unit and tool costs had low sensitivity values as their effects were mitigated from the minimization process. Optimization trend of ordered quantity over changing demand displayed a correlation with stochastic changes premised on a smoothing factor. This is shown in Figure 18 with the corresponding calculated production costs.

3.5. Discussion of Influence of Production Variables

From the results in Table 4 we can see that the time of operation of the machine as well as the demand are fundamental to deciding the total cost. The ordering cost in the numerical example shows some significant influence, however, this cost is not expected to exceed the machining costs during full production cycles. The core parameter of influence on machining is the cutting speed, which also affects the tool life. A scaled chart is shown in Figure 19 to illustrate the trends between factors and the influence of on cost. From the Figure 19, an increase in machining parameters (speed, feed rate) will cause a fall in tool-life and a corresponding increase in machining costs.

The implications of the research to industrial practice, shows an adequate method to estimate tools needed in a production system based on machining and clearly define the correlations between costs to process machining conditions. An enterprise could thus estimate create a production process which optimizes machining conditions while at the same time finding the optimal settings to tool allocation for storage.

Figure 18. Optimal tool ordering quantity and cost for stochastic demand data.

Figure 19. Graphical illustration of cost components against machining parameters and tool-life.

4. Conclusions

This study has presented an inventory control analysis premised on stochastic demand of machining tools. The model developed in this research is an optimum cost model governed by experimented machining conditions. Furthermore, the model is proposed to serve as a tool management policy for manufacturing industries. The fundamental considerations of the developed model includes the cutting tool, machining conditions and applied constraints. All of these are related to realistic conditions adapted to a practical production environment.

Optimization of the model solution was performed using Lingo, a commercially available optimization software. The simulation results obtained indicates that the feed rate is the most significant cutting parameter following its effect on machining costs. This in turn has a significant incremental effect on the overall production cost. Machine uptime and demand also showed significant effect on the total production cost. The process of optimizing the demand quantity, resulted in a cost-efficient optimum demand range characterized with a low safety stock and optimum ordering quantity. Additionally, an increase in the level of demand beyond the optimum range indicated that higher safety stock would be required. The application of Holt-Winters exponential smoothing forecasting technique further validates the outputs from this research.

Some limitations apply to the results of this research. Some of which are that the study focuses mainly on the use of a single cutting operation. In addition, an in-depth breakdown of cost influences from personnel, tool refurbishing and penalty costs are not considered. Fluctuations in ordering, holding, and unexpected inventory conditions are also not assessed.

In conclusion, several practical applications of this model can be obtained. This beginning with a single production operation concept as illustrated in this research to a more robust multi-operation manufacturing assembly. Future considerations of this research involve the development of a multiple-operational production approach with the inclusion of labor and penalty costs.

Appendix

1) Cost function components and intermediate variables

${C}_{op}$ : Total cost of operation ($)

${C}_{sh}$ : Total shortage cost ($)

${C}_{B}$ : Total workpiece material cost ($)

${C}_{h}$ : Total holding cost ($)

${C}_{o}$ : Total ordering cost ($)

TC: Total production cost ($)

2) Decision variables

${Q}_{p}$ : Quantity of cutting tools required for production (tool/order)

${Q}_{r}$ : Quantity of cutting tools required for roughing operation (tool/order) in period i

${Q}_{f}$ : Quantity of cutting tools for finishing operation (tool/order)

${T}_{c}$ : Tool ordering cycle length i.e. (length of time between placement of replenishment of tool orders) (time units)

3) Cost function components and intermediate variables

${c}_{hr},{c}_{hf}$ : Tool utilization cost per hour ($/hr)

${c}_{h}$ : Expected unit holding cost ($/tool)

${c}_{B}$ : Cost per unit of workpiece ($/unit)

${c}_{pr},{c}_{pf}$ : Tool Market price ($/tool)

${k}_{0}$ : Cost of machine use per hour ($/hr)

${c}_{sh}$ : Expected total annual shortage cost ($/tool)

L: Length per unit of workpiece to be cut (mm)

W: Width per unit of workpiece to be cut (mm)

H: Height per unit of workpiece to be cut (mm)

${w}_{r}$ : Width (radial) of cut (mm)

${d}_{r}$ : Depth (axial) of cut (mm)

${D}_{t}$ : Diameter of the tool (mm)

${t}_{opr},{t}_{opf}$ : Cutting time per unit of workpiece, roughing and finishing operations, (min/unit)

${t}_{0r},{t}_{0f}$ : Tool life, roughing and finishing operations (min)

k: Taylor’s equation constant

${v}_{r},{v}_{f}$ : Cutting speed, roughing and finishing operations (m/min)

${f}_{f}{f}_{r}$ : Feed per tooth, roughing and finishing operations (mm/tooth)

$\alpha ,\beta ,\gamma ,\lambda $ : Tool life parameters

${Q}_{t}$ : Demand for tool required per year (in units)

${Q}_{s}$ : Quantity of tool shortage (units)

${Q}_{B}$ : Quantity of material ordered (units)

D: Demand rate of products (products/year)

${F}_{t}$ : Amount of missing tool in system (units)

${I}_{i}$ : Average inventory at the end of the current period (units)

${I}_{i-1}$ : Previous inventory (units)

${A}_{i}$ : Total capacity of inventory (units)

$\stackrel{\xaf}{{I}_{i}}$ : Possible holding inventory (units)

${z}_{i}$ : Factor indicating variable ordering cost (1,0)

References

[1] Tomotani, J.V. and de Mesquita, M.A. (2018) Lot Sizing and Scheduling: A Survey of Practices in Brazilian Companies. Production Planning & Control, 29, 236-246.

https://doi.org/10.1080/09537287.2017.1409370

[2] Gershwin, S.B. (2018) The Future of Manufacturing Systems Engineering. International Journal of Production Research, 56, 224-237.

https://doi.org/10.1080/00207543.2017.1395491

[3] Conceicao, S.V., et al. (2015) A Demand Classification Scheme for Spare Part Inventory Model Subject to Stochastic Demand and Lead Time. Production Planning & Control, 26, 1318-1331.

https://doi.org/10.1080/09537287.2015.1033497

[4] Dimla Snr, D.E. (2000) Sensor Signals for Tool-Wear Monitoring in Metal Cutting Operations—A Review of Methods. International Journal of Machine Tools and Manufacture, 40, 1073-1098.

https://doi.org/10.1016/S0890-6955(99)00122-4

[5] Teti, R., et al. (2010) Advanced Monitoring of Machining Operations. CIRP Annals-Manufacturing Technology, 59, 717-739.

https://doi.org/10.1016/j.cirp.2010.05.010

[6] Arrazola, P.J., et al. (2013) Recent Advances in Modelling of Metal Machining Processes. CIRP Annals-Manufacturing Technology, 62, 695-718.

https://doi.org/10.1016/j.cirp.2013.05.006

[7] Byrne, G., et al. (1995) Tool Condition Monitoring (TCM)—The Status of Research and Industrial Application. CIRP Annals-Manufacturing Technology, 44, 541-567.

https://doi.org/10.1016/S0007-8506(07)60503-4

[8] Kouedeu, A.F., et al. (2014) Stochastic Optimal Control of Manufacturing Systems under Production-Dependent Failure Rates. International Journal of Production Economics, 150, 174-187.

https://doi.org/10.1016/j.ijpe.2013.12.032

[9] Li, C.R., Sarker, B.R. and Yi, H.Z. (2014) An Optimal Stocking Policy for Machining Tools with Stochastically Distributed Lifespan and Demand. International Journal of Production Research, 52, 6175-6191.

https://doi.org/10.1080/00207543.2014.918291

[10] Li, C.R. and Cheng, J.D. (2015) An Optimal Inventory Policy under Certainty Distributed Demand for Cutting Tools with Stochastically Distributed Lifespan. Cogent Engineering, 2.

[11] Parent, L., Songmene, V. and Kenné, J.P. (2007) A Generalised Model for Optimising an End Milling Operation. Production Planning & Control, 18, 319-337.

https://doi.org/10.1080/09537280701292501

[12] Conradie, P., Dimitrov, D. and Oosthuizen, G. (2016) A Cost Modelling Approach for Milling Titanium Alloys. 7th Hpc 2016-Cirp Conference on High Performance Cutting, 46, 412-415.

[13] Braziotis, C., Tannock, J.D.T. and Bourlakis, M. (2017) Strategic and Operational Considerations for the Extended Enterprise: Insights from the Aerospace Industry. Production Planning & Control, 28, 267-280.

https://doi.org/10.1080/09537287.2016.1268274

[14] Neugebauer, R., et al. (2012) Resource and Energy Efficiency in Machining Using High-Performance and Hybrid Processes. Fifth Cirp Conference on High Performance Cutting, 1, 3-16.

https://doi.org/10.1016/j.procir.2012.04.002

[15] Andrew-Munot, M. and Ibrahim, R.N. (2013) Development and Analysis of Mathematical and Simulation Models of Decision-Making Tools for Remanufacturing. Production Planning & Control, 24, 1081-1100.

https://doi.org/10.1080/09537287.2012.654667

[16] Syntetos, A.A., Boylan, J.E. and Disney, S.M. (2009) Forecasting for Inventory Planning: A 50-Year Review. Journal of the Operational Research Society, 60, S149-S160.

https://doi.org/10.1057/jors.2008.173

[17] Shang, J., et al. (2008) A Decision Support System for Managing Inventory at GlaxoSmithKline. Decision Support Systems, 46, 1-13.

[18] Wang, Y.-C., et al. (2018) Optimization of Machining Economics and Energy Consumption in Face Milling Operations. The International Journal of Advanced Manufacturing Technology, 99, 2093-2100.

https://doi.org/10.1007/s00170-018-1848-6

[19] Wang, J., et al. (2002) Optimization of Cutting Conditions for Single Pass Turning Operations Using a Deterministic Approach. International Journal of Machine Tools and Manufacture, 42, 1023-1033.

https://doi.org/10.1016/S0890-6955(02)00037-8

[20] Mascle, C. and Gosse, J. (2014) Inventory Management Maximization Based on Sales Forecast: Case Study. Production Planning & Control, 25, 1039-1057.

https://doi.org/10.1080/09537287.2013.805343

[21] Taylor, J.W. (2003) Short-Term Electricity Demand Forecasting Using Double Seasonal Exponential Smoothing. Journal of the Operational Research Society, 54, 799-805.

https://doi.org/10.1057/palgrave.jors.2601589

[22] Ferbar, L. and Vehove, A. (2013) The Improvement of the Holt-Winters Method for Intermittent Demand: A Case of Overnight Stays of Turists for Some Community in Republic of Slovenia. 45-52.

[23] Cunningham, C.R., et al. (2017) Cost Modelling and Sensitivity Analysis of Wire and Arc Additive Manufacturing. Procedia Manufacturing, 11, 650-657.

https://doi.org/10.1016/j.promfg.2017.07.163