Abstract
The specific energy consumption of machine tools and surface roughness are important indicators for evaluating energy consumption and surface quality in processing. Accurate prediction of them is the basis for realizing processing optimization. Although tool wear is inevitable, the effect of tool wear was seldom considered in the previous prediction models for specific energy consumption of machine tools and surface roughness. In this paper, the prediction models for specific energy consumption of machine tools and surface roughness considering tool wear evolution were developed. The cutting depth, feed rate, spindle speed, and tool flank wear were featured as input variables, and the orthogonal experimental results were used as training points to establish the prediction models based on support vector regression (SVR) algorithm. The proposed models were verified with wet turning AISI 1045 steel experiments. The experimental results indicated that the improved models based on cutting parameters and tool wear have higher prediction accuracy than the prediction models only considering cutting parameters. As such, the proposed models can be significant supplements to the existing specific energy consumption of machine tools and surface roughness modeling, and may provide useful guides on the formulation of cutting parameters.
Keywords
Introduction
As the basic devices in machining system, machine tools have large quantity with high energy consumption. However, their energy efficiency is very low, with an average of less than 30%. 1 In addition, the greenhouse gas emissions caused by the use of machine tools cannot be ignored. 2 Gutowski 3 pointed out that the CO2 emissions of a computer numerical control (CNC) machine tool with spindle power in 22 kW operating 1 year is equivalent to the emissions of 61 sport utility vehicles (SUVs). Consequently, the urgent task of manufacturing industry is how to reduce energy consumption and carbon emission to achieve sustainable development on the premise of ensuring the processing quality.4,5
Reducing the energy consumption in machining can significantly improve the environmental performance of manufacturing industry. 6 In this regard, Rao 7 presented some schemes for reducing the energy consumed in grinding, which paves the way for improved energy-efficient manufacturing processes. The basic energy-saving approach is to monitor and evaluate energy consumption of machine tools. 8 Therefore, many scholars have investigated energy consumption characteristics of machine tools and made valuable contributions toward reliable and accurate energy consumption and power demand modeling. For example, Draganescu et al. 9 and Zhang et al. 10 presented energy consumption models of machine tools from the viewpoint of cutting force, respectively. Gutowski et al. 11 proposed a total power model of machine tools as a linear function of the material removal rate (MRR). Based on the relationship between energy consumption of machine tools and MRR, Kara and Li obtained an energy consumption prediction model suitable for turning and milling.12,13 It is worth noticing that MRR is the only independent variable in aforementioned models. However, the cutting experiments indicate that the energy consumption of machine tools under the same MRR and different cutting parameters is not always equal. In this regard, Zhou et al. 14 compared the prediction accuracy between the model based on MRR and the model based on cutting parameters. Their conclusion is that the latter model with cutting parameters as independent variables obtains higher prediction accuracy. Li et al. 15 developed a model as a function of MRR and spindle speed, which can realize the prediction of energy consumption in face milling. Besides, He et al. 16 presented a method to estimate the energy consumption of machine tools, which can be used to optimize the numerical control (NC) codes to improve energy efficiency. Most of the preceding studies characterize the relationship between the energy consumption of machine tools and cutting parameters based on empirical modeling. However, these traditional energy consumption models fail to take the effect of tool wear into consideration.
Actually, the cutting experiments with the same workpiece material and cutting parameters show that the cutting power is different as the tool wears gradually. Axinte and Gindy assessed the sensitivity of spindle power for tool wear in turning, milling, and drilling. 17 For the purpose of tool condition monitoring, Cuppini et al. 18 analyzed the relationship between cutting power of machine tools and tool wear. In addition, Li et al. 19 established a regression model of cutting power based on cutting parameters and tool wear using the response surface method (RSM) for tool change decision. Although these studies are not intended to estimate the energy consumption of machine tools, they indicate that tool wear can affect the energy consumption of machine tools in machining.
Recently, several scholars have taken tool wear into consideration and made further research on the energy consumption of machine tools. For instance, Yoon et al. 20 investigated the energy consumption of machine tools under slight, medium, and severe tool wear conditions. Their experimental results showed that the increase of tool wear leads to more energy consumption in machining. Additionally, Liu et al. 21 explored the effect of MRR and tool wear on energy consumption at the machine tool, spindle, and process levels in hard milling AISI H13. Their conclusion is that the tool wear has the significant influence on the specific energy consumption of the process level. However, these models were established with few conditions of tool wear, which cannot fully reflect the effect of tool wear evolution on energy consumption of machine tools. In this regard, Zhao et al. 22 developed an empirical model to predict tool tip cutting specific energy consumption based on cutting parameters and tool flank wear in milling. In fact, the condition of tool wear changes constantly because of the friction between the tool and the workpiece in cutting stage.
Although sustainable manufacturing requires manufacturers to save energy and reduce emissions as much as possible, the premise is to guarantee the processing quality. The workpiece surface quality usually has a great influence on mechanical performance of machined parts. In general, the improvement of energy efficiency in machining will lead to the deterioration of surface roughness. In order to solve this problem, techniques such as hybrid decision analysis, Taguchi method, and gray correlation analysis are used for multi-objective optimization of cutting parameters.23–25 However, accurate prediction of surface roughness is the basis for processing optimization. Thus, as an important indicator for evaluating the processing quality, surface roughness also needs to be accurately predicted. In theory, surface roughness is closely related to cutting conditions (e.g. cutting parameters, cutting fluid, tool geometric parameters, and tool-workpieces properties). Various studies have made important contributions to the prediction and control of surface roughness. For instance, Arapoğlu et al. 26 and Imani et al. 27 established surface roughness prediction models based on artificial neural networks (ANN) in turning and milling, respectively. Debnath et al. 28 studied the influence of various cutting fluid levels and cutting parameters on surface roughness in turning. Patel and Gandhi built a mathematical model for surface roughness in turning AISI D2, which is exponentially related to feed rate, cutting speed, and tool nose radius. 29 In addition, RSM is often used as a powerful design-of-experiment (DoE) technique when dealing with surface roughness. Using RSM, Kosaraju and Anne 30 proposed a model for surface roughness based on cutting speed, feed rate, cutting depth, and back rake angle to select optimal cutting parameters in turning Ti-6Al-4V. Similarly, Parida and Maity 31 developed a model for surface roughness in hot turning Monel 400, which is based on cutting speed, feed rate, cutting depth, and workpiece temperature. However, these models were based on a common assumption that the condition of tool wear remains unchanged in machining, without considering the effect of tool wear evolution on surface roughness. In truth, worn tools harmfully affect surface roughness and the performance of machined parts.
In sum, the energy consumption of machine tools and surface roughness modeling with tool wear evolution have not been solved to a satisfactory degree, and further effort is needed for a more comprehensive solution. The support vector regression (SVR) is one of machine learning (ML) algorithms, 32 which can effectively solve the problem of nonlinear regression. Compared with ANN and random forest (RF), SVR has better prediction accuracy and shorter running time when the amount of training set is small. 33 Accordingly, it has been widely used in the field of machining, such as predicting surface roughness and tool wear.34,35
The objectives of this paper are to: (1) consider the effect of tool wear evolution on energy consumption of machine tools and surface roughness; (2) introduce SVR algorithm to develop the prediction models for specific energy consumption of machine tools and surface roughness based on cutting parameters and tool wear; and (3) verify the proposed models with wet turning experiments of AISI 1045 steel.
Experimental and SVR modeling details
Experimental design
CNC turning operations
As a kind of high-quality carbon structural steel, AISI 1045 steel is widely used in the production of bolts, shafts, and guide pins. It was selected as the workpiece material, and the chemical composition of AISI 1045 is shown in Table 1. The diameter of workpieces was 55 mm, with processing length 120 mm. Wet turning of AISI 1045 steel was performed at a CNC lathe (Yishui CKJ6163, Figure 1) with carbide cutting inserts (SKYWALKER CNMG120408-PG SC4025). The workpiece surface was turned prior to any experiments in order to remove surface defects.
The chemical composition of AISI 1045.

Yishui CKJ6163 lathe.
The ranges of cutting parameters (cutting depth ap, feed rate f, and cutting speed v) were customized within the ranges recommended by the cutting insert manufacturer. In order to reduce the cost of tests and improve the efficiency, the L25 (53) orthogonal experiments with 25 sets of three factors and five levels were designed with Taguchi method. Table 2 shows cutting parameters and their levels. The spindle speed n was used in CNC programming, it can be calculated by equation (1).
where n is spindle speed in r/min; v is cutting speed in m/min; and d is workpiece diameter in mm. Experimental combinations of the cutting parameters are presented in Table 3.
Cutting parameters and their levels.
Experiment design using L25 (53) orthogonal array.
Tool wear and surface roughness measurement
Tool flank wear is always considered to be the main subject of tool wear and the standard of tool bluntness. 36 Therefore, the flank wear of the cutting inserts was measured in this research. The tool wear measurement system (as shown in Figure 2) consisted of a microscope (ZEISS Axio Lab.A1 Mat) and its image analysis software. The flank wear of cutting inserts was measured twice before and after each turning, and the average value was taken as the tool wear VB. In addition, the maximum allowable tool flank wear has been set to be VBmax = 0.27 mm, which means that the cutting inserts will be replaced if the wear length reaches 0.27 mm. Since tool wear measurement required cutting inserts to be removed from the tool holder, the tool placement error was eliminated by re-setting the tool before the next experiment.

Tool wear measurement system.
The workpiece surface roughness was measured (as shown in Figure 3) with the surface roughness tester RTP120. The surface roughness tester was placed on the sliding table after each turning. Measurements from three different locations were carried out, and the average value was taken as the surface roughness Ra in each experiment.

Surface roughness measurements.
Machine tool power measurement and specific energy consumption calculation
The complete CNC lathe working process includes start-up, standby, spindle start, no-load operation, cutting stage, no-load operation, spindle stop, standby, and shutdown. In general, cutting stage is the most energy-consuming working process. In this stage, the energy-consuming components of machine tools include machine control unit (MCU), spindle motor, feed-axis motors, cooling pump motor, and lighting device. Therefore, the energy consumption of machine tools in cutting stage is the focus of this study.
The cutting power of machine tools was measured with a power analyzer (Yokogawa WT500) and shunts. The voltage signals and current signals of lathe were connected to the power analyzer, and a personal computer (PC) was connected to the power analyzer for recording and processing power signals. Figure 4 shows the hardware configuration used for the power measurement system.

Hardware configuration for the power measurement system.
The specific energy consumption was employed to evaluate the energy efficiency of machine tools, which is defined as the energy consumption to remove a unit volume of material. 13 Once the specific energy of machine tools is known, the total energy consumption of machine tools can be accurately predicted. The specific energy consumption of machine tools can be calculated by equation (2).
where SEC is the specific energy consumption of machine tools in J/mm3; E is the energy consumption of machine tools in J; Q is the total material removal volume in mm3; P is the total power of machine tools in cutting stage in W; MRR is material removal rate in mm3/s; and t is the cutting time in s. Moreover, the material removal rate in turning can be obtained by equation (3).
where d is workpiece diameter in mm; n is spindle speed in r/min; f is feed rate in mm/r; and ap is cutting depth in mm.
SVR modeling
In this section, SVR algorithm was introduced to establish the prediction models for SEC and Ra. The essence of SVR modeling is to find an optimal classification surface, so that the distance error from all training points to the classification surface is minimized. 32 For comparative analysis, two kinds of prediction models were established using the 25 sets of orthogonal experiment data. The first kind of prediction models (SEC1 and Ra1) was based on cutting parameters, the other proposed (SEC2 and Ra2) was based on cutting parameters and tool wear.
Running of SVR
The flowchart of SVR modeling is shown in Figure 5, including the following five steps.
Step 1: transform the experiment data to the format of an SVR package.
Step 2: conduct scaling on the original experiment data. The input variables are scaled to [−1, 1] to eliminate the dimension differences and improve the training efficiency of models.
Step 3: find the optimal parameters. The radial basis function (RBF, see equation (4)) is selected as the kernel function, and the key parameters of SVR algorithm, i.e., constant coefficient C, kernel function parameter g, and loss function parameter p are searched by cross-validation.
Step 4: use the optimal parameters to train the whole training set to obtain the prediction model.
Step 5: test.

The flowchart of SVR modeling.
Prediction models based on cutting parameters
According to metal cutting theory, the cutting power and surface roughness are closely related to cutting parameters. Hence, the prediction models (SEC1 and Ra1) based on cutting parameters were presented. There were three inputs and one output in SEC1 and Ra1. Specifically, ap, f, and n were featured as input variables, with SEC and Ra were taken as target values, respectively, as shown in Figure 6.

Prediction models based on cutting parameters.
Prediction models based on cutting parameters and tool wear
Using the same workpiece and cutting parameters, the cutting power and surface roughness will increase as the cutting insert gradually wears. Therefore, specific energy consumption of machine tools and surface roughness are not only related to cutting parameters, but also tool wear condition. The prediction models (SEC2 and Ra2) based on cutting parameters and tool wear were developed. There were four inputs and one output in SEC2 and Ra2. To be specific, ap, f, n, and VB were featured as input variables, with SEC and Ra were taken as target values, respectively, as shown in Figure 7.

Prediction models based on cutting parameters and tool wear.
Results and discussion
Experimental results
About 25 sets of experimental measurement and calculation results (VB, Ra, and SEC) are presented in Table 4.
Experimental measurement results.
Numerical results
Training results of SVR models
With the help of LIBSVM 3.23 toolbox in windows cmd.exe (Appendix A), 37 both kinds of prediction models were established. The selection of parameters and training results of SVR models are as follows.
For SEC1, the optimal parameters found by cross-validation are C = 64, g = 0.5, and p = 0.0625. Additionally, the mean square error MSE of the model is 0.177892, and the determination coefficient R2 is 0.985971. The comparison of measured-predicted values from SEC1 is depicted in Figure 8. With regard to Ra1, the optimal parameters found by cross-validation are C = 4, g = 0.25, and p = 0.5. The MSE of the model is 0.125451, and R2 is 0.660473. The comparison of measured-predicted values from Ra1 is depicted in Figure 9.

Comparison of measured-predicted values from SEC1.

Comparison of measured-predicted values from Ra1.
For SEC2, the optimal parameters found by cross-validation are C = 16, g = 2.2974, and p = 0.01. The MSE of the model is 1.00584e−4, and R2 is 0.99995. Figure 10 shows the comparison of measured-predicted values from SEC2. With regard to Ra2, the optimal parameters found by cross-validation are C = 32, g = 0.5, and p = 0.0078125. The MSE of the model is 6.0567e−5, and R2 is 0.9998450. Figure 11 shows the comparison of measured-predicted values from Ra2.

Comparison of measured-predicted values from SEC2.

Comparison of measured-predicted values from Ra2.
The mean square error MSE and the determination coefficient R2 are widely used to evaluate regression effect in SVR modeling. The smaller the value of MSE and the larger the value of R2, the better the model fitting effect and the higher the prediction accuracy. The regression evaluation indices of two kinds of prediction models are presented in Table 5. From the training results, the MSE of SEC2 and Ra2 are smaller than SEC1 and Ra1, and the R2 of SEC2 and Ra2 are larger than SEC1 and Ra1, which indicates that SEC2 and Ra2 achieve higher fitting accuracy.
Regression evaluation indices of two kinds of prediction models.
Test results of SVR models
New combinations of cutting parameters were designed to test the above two kinds of prediction models. The experimental cutting parameters and measurement results are presented in Table 6, and the test results when adopting two kinds of models are compared in Figures 12 and 13.
Cutting parameters of test experiments and measurement results.

Comparison of prediction accuracy of SEC models.

Comparison of prediction accuracy of Ra models.
The SEC measurement results are 3.6144 J/mm3 and 3.2064 J/mm3, respectively. For the test No. 1, the result using SEC1 is 4.8899 J/mm3 and the prediction accuracy is 64.71%, while the result using SEC2 is 3.5526 J/mm3 and the prediction accuracy is 98.29%. For the test No. 2, the result using SEC1 is 3.0401 J/mm3 and the prediction accuracy is 94.81%, while the result using SEC2 is 3.1419 J/mm3 and the prediction accuracy is 97.99%.
The Ra measurement results are 2.0380 µm and 2.1105 µm, respectively. For the test No. 1, the result using Ra1 is 2.3178 µm and the prediction accuracy is 86.27%, while the result using Ra2 is 1.9613 µm and the prediction accuracy is 96.24%. For the test No. 2, the result using Ra1 is 1.5978 µm and the prediction accuracy is 75.71%, while the result using Ra2 is 2.292 µm and the prediction accuracy is 91.40%. According to the test results, the developed models SEC2 and Ra2 achieve higher prediction accuracy than SEC1 and Ra1.
Since the contact between tool and workpiece will certainly initiate tool wear when sufficient cutting time is given, tool wear is unavoidable. Actually, it is an accumulating process including initial, normal, and sharp wear stage. Taking the first cutting insert in this study as an example, the tool flank wear evolution is presented in Figure 14, and the representative SEM images of tool flank wear can be seen in Figure 15. In cutting stage, as the tool flank wear increases, the contact area between tool flank and the machined surface adds. This will lead to an increase in the friction force and cutting temperature between tool and workpiece. More energy is expected to be consumed in cutting stage, and increased surface roughness on the workpiece will be generated.

Tool flank wear evolution (taking the first cutting insert in this study as an example).

Representative SEM images of tool flank wear: (a) initial wear stage, (b) normal wear stage, and (c) sharp wear stage.
Hence, the effect of tool wear evolution should be taken into account in SEC and Ra modeling. Featuring cutting parameters (ap, f, and n) and tool flank wear (VB) to establish SEC and Ra models based on SVR algorithm is a more comprehensive and reliable method. This is the reason why the developed models achieve higher fitting accuracy and prediction accuracy.
Conclusion
In this paper, the tool flank wear evolution was observed and recorded in the cutting stage. The SVR models for SEC and Ra based on cutting parameters and tool wear were developed and verified in wet turning AISI 1045 steel experiments. The following conclusions can be drawn.
The tool wear has an important effect on energy consumption of machine tools and workpiece surface quality. Therefore, the wear condition of cutting inserts is a significant factor that cannot be neglected in the modeling of SEC and Ra.
The prediction accuracy of developed models SEC2 and Ra2 is 97.99% and 91.40%, respectively. The improved models based on cutting parameters and tool wear are simple with high prediction accuracy, which are useful to formulate cutting parameters that can reduce energy consumption and improve surface quality simultaneously in turning.
Only the tool flank wear is considered in the modeling of SEC and Ra. In addition, the tool rake wear is serious in cutting stage. The effect of tool rake wear on SEC and Ra will be investigated in future.
Footnotes
Appendix A 37
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Project of Shandong Province Natural Science Foundation of China (No. ZR2016EEM29) and the Project of Shandong Province key research development of China (No. 2017GGX30114).
