Research Article | Open Access

Tian-Syung Lan, "Dynamic MRR Function Optimization Using Calculus of Variations", *Mathematical Problems in Engineering*, vol. 2010, Article ID 671062, 12 pages, 2010. https://doi.org/10.1155/2010/671062

# Dynamic MRR Function Optimization Using Calculus of Variations

**Academic Editor:**Jyh Horng Chou

#### Abstract

A dynamic machining model to optimize the control of material removal rate (MRR) for a cutting tool undergoing the considerations of fixed tool life and maximum machining rate is established in this paper. This study not only applies material removal rate mathematically into the objective function, but also implements Calculus of Variations to comprehensively optimize the control of material removal rate. In addition, the optimal solution for the dynamic machining model to gain the maximum profit is provided, and the decision criteria for selecting the optimal solution of the dynamic machining model are then recommended. Moreover, the computerized analyses to simulate both the dynamic and traditional machining models for a numerical example are also promoted. This study definitely contributes an applicable approach to the dynamic function of material removal rate and provides the efficient tool to concretely optimize the profit of a cutting tool for operation planning and control in modern machining industry with profound insight.

#### 1. Introduction

The cutting conditions of a cutting tool have been the most critical variables in the machining process. Cutting speed, feed rate, and depth of cut were considered as three factors of input cutting parameters. Liew and Ding [1] have expressed several methods to be used under stepwise constant variation in feed, speed, or depth of cut, but none is practically applicable when two or more cutting conditions are changed. Hence, the research of controlling cutting conditions with fixed material removal rate has been introduced [2, 3]. For most studies with this viewpoint, the material removal rate is fixed because of the expensive observation of control. However, through the computer-integrated interface to program the feed rate on modern computer numerical controlled (*CNC*) machines [2] with fixed cutting speed and depth of cut, the material removal rate is capable of being dynamically controlled.

In addition, Cámara et al. [3] described that tool life is also a critical parameter of the cutting process. Iqbal et al. [4] proposed a suitable implementation to predict tool life. Ferrer et al. [5] proposed a method of optimal control to ensure maximum tool life. Kim et al. [6] also provided a modified Taylor tool life equation to minimize tool cost. However, the maximum tool life or the minimum tool cost will not guarantee the maximum profit of a cutting tool. Besides, the various tool inspecting periods for a tool change from different machine tool operators will decrease the productivity and increase the cost during machining significantly. In order to well manage the productivity and the consumption of cutting tools, a fixed tool life is practically considered into the cutting process.

Moreover, the cost to machine each part is a function of the machining time, and the marginal cost of production is a linear function of production rate. Therefore, the marginal cost of operation for the machining process is also proposed to be a linear function of the material removal rate in this study. This explicates that more machining rate causes more operational cost such as machine maintenance, loading-unloading, and machine depreciation costs.

According to Zong et al. [7], the dependency of a reliability model on the cutting conditions is the aim to optimize the manufacturing system [8]. Although several time series modeling on the control of machining process are mentioned for decision making [9], none is capable to achieve the maximum profit. They are all emphasizing on the maximum tool usage or minimum tool cost. Practically, the profit and the productivity of machining operation are the mostly concerned problems confronting the manufacturing industry. Besides, the need of operating *CNC* machines efficiently to obtain the required payback is even greater in the case of rough machining, since a greater amount of material is removed thus increasing possible savings [10]. With the reasons above, there is an economic demand to control the material removal rate of rough machining operation on *CNC *machines. Thus, the essentiality of finding the optimal solution to achieve the maximum profit of a cutting tool with fixed tool life is absolutely arising.

#### 2. Assumptions and Notations

Before formulating the problem, several assumptions and notations are to be constructed. They are described as follows.

##### 2.1. Assumptions

(1)The cutting process is a continuous rough turning operation with one type of tool.(2)Each tool performs the same fixed length of cutting time (tool life). (3)The upper limit of material removal rate is generated from the maximum cutting conditions suggested in the handbook, and the fixed tool life is derived from the Taylor’s tool life equation [11] with these maximal conditions. Thus, no tool will break before this fixed tool life even with the upper limit of material removal rate.(4)The marginal cost of operation is a linear function of the material removal rate [12].(5)There is no chattering or scrapping of parts occuring during the machining process.(6)All chips from cutting and finished parts are held in the machine until a tool change.(7)All machined parts are paid at a given price immediately at the tool change.##### 2.2. Notations

One has the following:: average volume of material machined per unit part,: upper limit of material removal rate,: overall holding cost for unit chip per unit time in the machine, where ,: marginal operation cost at the material removal rate ; where *b* is a constant, : operational cost at time *t*,: chip holding cost for unit chip per unit time,: part holding cost for unit finished part per unit time,: contribution per unit part machined,: controlling time interval of the tool life.

##### 2.3. Decision Functions

One has the following: cumulated volume of material machined during time interval ,: material removal rate at time *t*.

#### 3. Model Formulation

When considering the machining profit of an individual cutting tool under fixed tool life, the mathematical model representing the machining profit under dynamic *MRR* control is formulated. In the mathematical model; represents the revenue of each machined product, is the selected fixed tool life of a cutting tool, describes the cumulated amount of material machined during , illustrates the *MRR* at time , denotes the cumulated volume of material machined before the fixed tool life , and indicates the average volume of material to be machined per unit part. Therefore, explains the contribution of an individual cutting tool under the machining operation with a fixed life *T*.

Since the marginal operation cost of production is a linear increasing function of production rate [12], the marginal operation cost of a machine is proposed in direct proportion to the *MRR* that is described as , where is the operational constant that is obtained through the equipment management strategy of the *CNC* machine. It can often be evaluated by dividing the initial purchasing cost of the *CNC* machine from the value of multiplying the expected life of the machine with the suggested average *MRR* from the equipment manufacturer.

Whilst the marginal operation cost of production is a linear increasing function of production rate and the operational cost is directly proportional to the square of the production rate [12], the operational cost of a machine is also proposed in direct proportion to the square of the material removal rate in this study. Thus, signifies the operational cost during time interval .

As all chips from cutting and finished products are usually held and stored at the machine until a tool replacement, denotes the chip holding cost during time interval , where is the chip holding cost per unit chip per unit time. And, denotes the finished part holding cost during time interval , where is the finished part holding cost per unit finished part per unit time. Therefore, the objective function and its constraints are listed as shown below:

Let be the overall holding cost for unit chip per unit time. Thus, the dynamic machining model is then developed as shown below:

#### 4. Optimal Solution

Let and be the optimal solution of the dynamic machining model, and assume that time interval is the maximal subinterval of [0,*T*] satisfying the Euler Equation [12, 13] .

There are two probable situations to be discussed in this study.

##### 4.1. Situation 1

*will not touch ** before T*.

The optimal solution for *Situation 1* is shown as follows:

The detail is developed in Appendix A.

From (4.1), it is found that the optimal material removal rate, , is a strictly increasing linear function of . In addition, by considering the constraint into (4.1) at , it is obtained that

This denotes that the overall holding cost of unit chip for the period *T* is smaller than the contribution per unit material machined. Otherwise, the optimal material removal rate will touch zero.

Applying the constraint into (4.1) at , we have

This represents that twice of the marginal cost at upper limit *B* must be larger than the contribution per unit material machined. Otherwise, the optimal material removal rate will touch the upper limit *B*.

Before solving the optimal solution for *Situation 2*, one Property is proposed and illustrated as follows.

*Property 1. *If the line touches the line , these two lines should overlap to be from the touch point to the end point .

The proof of the Property is illustrated in Appendix B.

##### 4.2. Situation 2

*will touch **at time *: .

The optimal solution for *Situation 2* is shown as follows:

The detail is developed in Appendix C.

##### 4.3. Decision Criteria

From (4.4) and (4.5), the maximum value of is found at when . Therefore, the following criteria are then made.

When , it means that . This denies the assumption of *Situation 2*. It is that the optimal material removal rate will not reach the upper limit *B* within the tool life. The optimal solution is *Situation 1*.

When , it means that . This is that the optimal material removal rate will reach the upper limit *B* within the tool life. Therefore, the optimal solution is *Situation 2*.

#### 5. Computerized Analyses

For a specific turning operation, there are bounds for cutting conditions suggested in the machining handbook. Therefore, there must exist a maximum material removal rate *U *and a minimum material removal rate *L *derived from the maximum and the minimum cutting conditions, respectively. Each material removal rate between *U* and *L* can feasibly be selected as the upper limit *B *for the dynamic machining model. From the well-known Taylor’s expression of the tool life [8], it is then modified as , when the cutting speed and depth of cut are selected fixed. Therefore, with this reconstructed equation, each fixed tool life is then obtained with respect to every feasible upper limit *B *selected. Thus, for each feasible upper limit, there exist an optimal profit for the dynamic machining model. In addition, the traditional machining model with a fixed *MRR* between *U* and *L* is also potential to be optimal. To locate the maximum profit solution among all these possibilities, a computer program written in MATLAB to analyze the problem is then constructed. The concept of the flow chart is shown as follows: *P, a, b, c, U, L, n, k *should be given before the following algorithm. Initialize *. *

*Step 1. *Compute *T* for both models, then compute the profit
and the production quantity for the traditional model.

Go to *Step 2*.

*Step 2. *Plot and , then go to Step 3.

*Step 3. *If , go to Step 5.

Otherwise, go to Step 4.

*Step 4. *Compute the profit
and the production quantity for *Situation 1* of the dynamic model.

Go to Step 6.

*Step 5. *Compute , the profit
and the production quantity for *Situation 2 *of the dynamic model.

Go to Step 6.

*Step 6. *Plot and , then go to Step 7.

*Step 7. *If , stop the program.

Otherwise, set as initialize, and return to Step 1.

#### 6. Simulated Results

With a machining project from AirTAC Corporation in Taipei, Taiwan, the numerical example referenced to a rough turning operation of specific cylindrical parts with step diameters is studied. The machining operation is assigned to a *CNC* lathe with FANUC controller. All data collected are transformed and listed as follows:

The simulated results of the MATLAB program are presented as the profit analysis and production quantity analysis in Figures 2 and 3, respectively. From Figure 2, it is noticed that the dynamic model is superior to the traditional model on the aspect of profit, when the selected upper limit *B *exceeds approximately . In addition, the profit of dynamic model is much more steady within the whole range of allowable machining rate, since the profit of traditional model is greatly decreasing as the *MRR* increases. Therefore, it is educed that the dynamic machining model is much more profitable than the traditional machining model for the higher range of machining rate, which is usually used in the real-world machining industry for higher productivity.

From Figure 3, it is also observed that the production quantity of a cutting tool under dynamic function is always less than what is under traditional control. However, as the machining rate increases, the difference of production quantity between two models is discovered to be decreasing. Thus, it is derived that the dynamic function will turn more competitive to the traditional control on the aspect of production quantity and productivity for the higher range of machining rate.

With the analyses of the dynamic model above, when there is a need for the cutting tool to acquire other productivity than the maximum profit solution, it is always feasible to slightly compromise on the profit for the required productivity. This will fractionally change the profit, but satisfactorily increases the productivity for the dynamic machining model.

#### 7. Conclusions

The tool life, operational cost, holding costs, contribution per unit part machined, average volume of material machined per unit part, and upper speed limit are considered simultaneously to dynamically optimize the control of material removal rate. It is an extremely hard-solving and complicated issue. However, through the dynamic machining model, the problem becomes concrete and solvable.

In addition, four characteristics of the dynamic machining model are illustrated as follows. First, the optimal material removal rate is a strictly increasing linear function of before reaching the upper limit *B*. Second, by PROPERTY, if the optimal material removal rate touches the upper limit *B*, it will stay to be *B* for the rest of tool life *T*. Third, by (4.2), the overall holding cost of unit chip for the tool life period *T* must be smaller than the contribution per unit material machined. Otherwise, the optimal material removal rate will touch zero. Fourth, with (4.3), two times of the marginal cost at upper limit *B* must be larger than the contribution per unit material machined. Otherwise, the optimal material removal rate will touch the upper limit *B*.

Moreover, from the computerized analyses and numerical simulation with the MATLAB program, the four remarks are then provided. First, the dynamic model receives much more balanced profit than the traditional model within the whole range of allowable machining rate. Second, the dynamic machining model is much more profitable than the traditional machining model for the higher range of machining rate. Third, when the machining rate increases, the difference of production quantity between two models will decrease. Fourth, when there is a need to obtain higher productivity than the optimal solution, it is always feasible for the dynamic machining model to slightly compromise on the profit in achieving the productivity. With these remarks above, the applicability and flexibility of the dynamic machining model are significantly expended.

The material removal rate is an important control factor of machining operation, and the control of machining rate is also critical for production planning. This study not only delivers the idea of controlling the material removal rate to the machining technology, but also leads a cutting tool towards to achieve maximum profit. This study focuses on the modeling and optimization of the *MRR* in *CNC* turning. Future researches with the cutting experiments to prove the proposed model as well as the modeling of dynamic optimization on multitool machining process and machining project control are absolutely encouraged. In summary, this study surely generates a reliable and applicable concept of machining control to the techniques, and also provides a better and practical solution to this field.

#### Appendices

#### A. Appendix A

The optimal solution for *Situation 1*.

Suppose that the material removal ratewill not reach the upper limit* B *before tool life *T*.

From the Euler Equation [12, 13], , it is derived that

There is a constant satisfying

With the transversality condition of salvage value for free end value [12, 13], it is obtained that
Integrating (A.2) with* t*, then
Using the boundary condition and (A.4), we have
Comparing (A.3) and (A.2) at , hence

Substituting (A.5) and (A.6) into (A.2) and (A.4), the optimal solution for *Situation 1* is then derived.

#### B. Appendix B

The proof of Property.

*Proof. *From (4.1), the optimal material removal rate is a strictly increasing linear function of . It holds for any subinterval of subject to . Therefore, the straight line in the time interval (shown in Figure 1) cannot exist, because it contradicts the Euler Equation [12, 13] to be a decreasing linear function of *t*, the PROPERTY is then verified.

#### C. Appendix C

The optimal solution for *Situation 2*.

It is assumed that the material removal rate will touch the upper limit *B *before tool life *T*. Before reaching the upper limit, the optimal solution for *Situation 1* is satisfied for this situation either.

The objective function of the dynamic machining model is then modified as

With the transversality condition of salvage value for free end point [12, 13] and Property; the time to reach the upper limit is derived.

Thus, with and the Property; the optimal solution for *Situation 2 *is then obtained.

#### Acknowledgment

The author would like to thank the anonymous referees who kindly provided the suggestions and comments to improve this paper.

#### References

- W. Y. H. Liew and X. Ding, “Wear progression of carbide tool in low-speed end milling of stainless steel,”
*Wear*, vol. 265, no. 1-2, pp. 155–166, 2008. View at: Publisher Site | Google Scholar - W. Grzesik, “Influence of tool wear on surface roughness in hard turning using differently shaped ceramic tools,”
*Wear*, vol. 265, no. 3-4, pp. 327–335, 2008. View at: Publisher Site | Google Scholar - M. Cámara, J. Ortega, and F. de Toro, “A single front genetic algorithm for parallel multi-objective optimization in dynamic environments,”
*Neurocomputing*, vol. 72, no. 16–18, pp. 3570–3579, 2009. View at: Publisher Site | Google Scholar - A. Iqbal, H. Ning, I. Khan, L. Liang, and N. U. Dar, “Modeling the effects of cutting parameters in MQL-employed finish hard-milling process using D-optimal method,”
*Journal of Materials Processing Technology*, vol. 199, no. 1–3, pp. 379–390, 2008. View at: Publisher Site | Google Scholar - I. Ferrer, J. Rios, and J. Ciurana, “An approach to integrate manufacturing process information in part design phases,”
*Journal of Materials Processing Technology*, vol. 209, no. 4, pp. 2085–2091, 2009. View at: Publisher Site | Google Scholar - D. W. Kim, Y. S. Lee, M. S. Park, and C. N. Chu, “Tool life improvement by peck drilling and thrust force monitoring during deep-micro-hole drilling of steel,”
*International Journal of Machine Tools and Manufacture*, vol. 49, no. 3-4, pp. 246–255, 2009. View at: Publisher Site | Google Scholar - W. J. Zong, Z. Q. Li, T. Sun, K. Cheng, D. Li, and S. Dong, “The basic issues in design and fabrication of diamond-cutting tools for ultra-precision and nanometric machining,”
*International Journal of Machine Tools and Manufacture*, vol. 50, no. 4, pp. 411–419, 2010. View at: Publisher Site | Google Scholar - J. V. Kadam, M. Schlegel, B. Srinivasan, D. Bonvin, and W. Marquardt, “Dynamic optimization in the presence of uncertainty: from off-line nominal solution to measurement-based implementation,”
*Journal of Process Control*, vol. 17, no. 5, pp. 389–398, 2007. View at: Publisher Site | Google Scholar - J. Chae and S. S. Park, “High frequency bandwidth measurements of micro cutting forces,”
*International Journal of Machine Tools and Manufacture*, vol. 47, no. 9, pp. 1433–1441, 2007. View at: Publisher Site | Google Scholar - S. Binsaeid, S. Asfour, S. Cho, and A. Onar, “Machine ensemble approach for simultaneous detection of transient and gradual abnormalities in end milling using multisensor fusion,”
*Journal of Materials Processing Technology*, vol. 209, no. 10, pp. 4728–4738, 2009. View at: Publisher Site | Google Scholar - I. Deiab, K. Assaleh, and F. Hammad, “On modeling of tool wear using sensor fusion and polynomial classifiers,”
*Mechanical Systems and Signal Processing*, vol. 23, no. 5, pp. 1719–1729, 2009. View at: Publisher Site | Google Scholar - A. Chiang,
*Dynamic Optimization*, McGraw-Hill, New York, NY, USA, 1992. - Q. Meng, J. A. Arsecularatne, and P. Mathew, “Calculation of optimum cutting conditions for turning operations using a machining theory,”
*International Journal of Machine Tools & Manufacture*, vol. 40, no. 12, pp. 1709–1733, 2000. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2010 Tian-Syung Lan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.