# Dynamic MRR Function Optimization Using Calculus of Variations

Authors: Tian-Syung Lan
Journal: Mathematical Problems in Engineering
 (2010)

---

## 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.


---

## Body

## 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 CNCmachines. 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:a: average volume of material machined per unit part,B: upper limit of material removal rate,c: overall holding cost for unit chip per unit time in the machine, where c=h1+(h2/a),bM′(t): marginal operation cost at the material removal rate M'(t); where b is a constant,bM'2(t): operational cost at time t,h1: chip holding cost for unit chip per unit time,h2: part holding cost for unit finished part per unit time,P: contribution per unit part machined,[0,T]: controlling time interval of the tool life.


### 2.3. Decision Functions

One has the followingM(t): cumulated volume of material machined during time interval [0,t],M′(t): material removal rate at time t.


## 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:a: average volume of material machined per unit part,B: upper limit of material removal rate,c: overall holding cost for unit chip per unit time in the machine, where c=h1+(h2/a),bM′(t): marginal operation cost at the material removal rate M'(t); where b is a constant,bM'2(t): operational cost at time t,h1: chip holding cost for unit chip per unit time,h2: part holding cost for unit finished part per unit time,P: contribution per unit part machined,[0,T]: controlling time interval of the tool life.


## 2.3. Decision Functions

One has the followingM(t): cumulated volume of material machined during time interval [0,t],M′(t): 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 dynamicMRR control is formulated. In the mathematical model; P represents the revenue of each machined product, T is the selected fixed tool life of a cutting tool, M(t) describes the cumulated amount of material machined during [0,t], M′(t) illustrates the MRR at time t, M(T) denotes the cumulated volume of material machined before the fixed tool life T, and a indicates the average volume of material to be machined per unit part. Therefore, P(M(T)/a) 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 bM′(t), where b 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, ∫0TbM'2(t)dt signifies the operational cost during time interval [0,T].As all chips from cutting and finished products are usually held and stored at the machine until a tool replacement,∫0Th1M(t)dt denotes the chip holding cost during time interval [0,T], where h1 is the chip holding cost per unit chip per unit time. And, ∫0Th2(M(t)/a)dt denotes the finished part holding cost during time interval [0,T], where h2 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:(3.1)Max{PM(T)a-∫0T[bM'2(t)+h1M(t)+h2M(t)a]dt}.Letc=h1+(h2/a) be the overall holding cost for unit chip per unit time. Thus, the dynamic machining model is then developed as shown below:  (3.2)Max{PM(T)a-∫0T[bM'2(t)+cM(t)]dt},s.t.M(0)=0,M(T)isfree,0≤M′(t)≤B,∀t∈[0,T].


## 4. Optimal Solution

LetM*′(t) and M*(t) be the optimal solution of the dynamic machining model, and assume that time interval (0,t̃) is the maximal subinterval of [0,T] satisfying the Euler Equation [12, 13] FM(t,M*(t),M*′(t))=dFM'(t,M*(t),M*′(t))/dt.There are two probable situations to be discussed in this study.


### 4.1. Situation 1

M*′(t)will  not touchBbefore T.The optimal solution forSituation 1 is shown as follows:(4.1)M*′(t)=c2bt+12b(Pa-cT),M*(t)=c4bt2+12b(Pa-cT)t.The detail is developed in AppendixA.From (4.1), it is found that the optimal material removal rate, M*′(t), is a strictly increasing linear function of t. In addition, by considering the constraint M′(t)≥0 into (4.1) at t=0, it is obtained that (4.2)cT≤Pa.This denotes that the overall holding cost of unit chip for the periodT is smaller than the contribution per unit material machined. Otherwise, the optimal material removal rate will touch zero.Applying the constraintM′(t)≤B into (4.1) at t=T, we have(4.3)Pa≤2bB.This represents that twice of the marginal cost at upper limitB 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 forSituation 2, one Property is proposed and illustrated as follows.Property 1.
If the liney=M*′(t) touches the line y=B, these two lines should overlap to be y=B from the touch point t̃ to the end point T.The proof of the Property is illustrated in AppendixB.


### 4.2. Situation 2

M*′(t)will  touchBat  timet̃: t̃∈[0,T).The optimal solution forSituation 2 is shown as follows:(4.4)t̃=T-1c(Pa-2bB),M*′(t)={c2bt+12b(Pa-cT),ift∈[0,t̃],B,ift∈(t̃,T],(4.5)M*(t)={c4bt2+12b(Pa-cT)t,ift∈[0,t̃],M(t̃)+B(t-t̃),ift∈(t̃,T].The detail is developed in AppendixC.


### 4.3. Decision Criteria

From (4.4) and (4.5), the maximum value of M*(t) is found at t̃=0 when t=T. Therefore, the following criteria are then made.WhenP/a≤2bB, it means that t̃≥T. 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.WhenP/a>2bB, it means that t̃<T. 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.


## 4.1. Situation 1

M*′(t)will  not touchBbefore T.The optimal solution forSituation 1 is shown as follows:(4.1)M*′(t)=c2bt+12b(Pa-cT),M*(t)=c4bt2+12b(Pa-cT)t.The detail is developed in AppendixA.From (4.1), it is found that the optimal material removal rate, M*′(t), is a strictly increasing linear function of t. In addition, by considering the constraint M′(t)≥0 into (4.1) at t=0, it is obtained that (4.2)cT≤Pa.This denotes that the overall holding cost of unit chip for the periodT is smaller than the contribution per unit material machined. Otherwise, the optimal material removal rate will touch zero.Applying the constraintM′(t)≤B into (4.1) at t=T, we have(4.3)Pa≤2bB.This represents that twice of the marginal cost at upper limitB 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 forSituation 2, one Property is proposed and illustrated as follows.Property 1.
If the liney=M*′(t) touches the line y=B, these two lines should overlap to be y=B from the touch point t̃ to the end point T.The proof of the Property is illustrated in AppendixB.


## 4.2. Situation 2

M*′(t)will  touchBat  timet̃: t̃∈[0,T).The optimal solution forSituation 2 is shown as follows:(4.4)t̃=T-1c(Pa-2bB),M*′(t)={c2bt+12b(Pa-cT),ift∈[0,t̃],B,ift∈(t̃,T],(4.5)M*(t)={c4bt2+12b(Pa-cT)t,ift∈[0,t̃],M(t̃)+B(t-t̃),ift∈(t̃,T].The detail is developed in AppendixC.


## 4.3. Decision Criteria

From (4.4) and (4.5), the maximum value of M*(t) is found at t̃=0 when t=T. Therefore, the following criteria are then made.WhenP/a≤2bB, it means that t̃≥T. 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.WhenP/a>2bB, it means that t̃<T. 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 rateUand a minimum material removal rate Lderived 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 Bfor the dynamic machining model. From the well-known Taylor’s expression of the tool life [8], it is then modified as B×Tn=k, when the cutting speed and depth of cut are selected fixed. Therefore, with this reconstructed equation, each fixed tool life T is then obtained with respect to every feasible upper limit Bselected. 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, kshould be given before the following algorithm. Initialize B=L.Step 1.
ComputeT for both models, then compute the profit
(5.1)Ob=P(BT/a)-(bB2T+(cB/2)T2)
and the production quantity Q=BT/a for the traditional model.
Go toStep2.Step 2.
PlotC(B,Ob) and F(B,Q), then go to Step 3.Step 3.
IfP/a>2bB, go to Step 5.
Otherwise, go to Step4.Step 4.
Compute the profit(5.2)O̅b=(PM(T)/a)-∫0T[bM'2(t)+cM(t)]dt
and the production quantity Q̅=M(T)/a for Situation 1 of the dynamic model.
Go to Step6.Step 5.
Computet̃, the profit
(5.3)O̅b=(PM(T)/a)-∫0t̃[bM'2(t)+cM(t)]dt-∫t̃T[bB2+c(M(t̃)+∫t̃tBdt)]dt
and the production quantity Q̅=M(T)/a for Situation 2of the dynamic model.
Go to Step6.Step 6.
PlotC̅(B,O̅b) and F̅(B,Q̅), then go to Step 7.Step 7.
IfB≥U, stop the program.
Otherwise, setB=B+0.05 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 aCNC lathe with FANUC controller. All data collected are transformed and listed as follows: (6.1)P=13.00dollars(US),a=1.00in3,b=1.05dollars(US)-minin3,c=0.03dollars(US)in3,U=12.00in3min,L=10.00in3min,n=0.2,k=16.0.The simulated results of the MATLAB program are presented as the profit analysis and production quantity analysis in Figures2 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 Bexceeds approximately 10.45(in3/min). 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 Figure3, 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 rateM*′(t) is a strictly increasing linear function of t before reaching the upper limit B. Second, by PROPERTY, if the optimal material removal rate M*′(t) 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 theMRR 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.



---

*Source: 671062-2010-06-24.xml*