Sustainability-Based Optimization of the Rotary Turning of the Hardened Steel

: The rotary turning is an e ﬀ ective manufacturing method to machine hardened metals due to longer tool life, higher production rate, and acceptable quality. However, sustainability-based optimization of the rotary turning has not been thoroughly considered because of the huge e ﬀ orts. This study presents an optimization to enhance the energy e ﬃ ciency ( EF R ), turning cost ( CT ), average roughness ( R a ), and the operational safety ( POS ) for the rotary turning of the hardened steel. Four key process parameters considered are the inclined angle ( α ), depth of cut ( a p ), feed rate ( f ), and cutting speed ( v c ). The improved Kriging ( IK ) models were used to construct the relations between the parameters and performances. The optimum varied factors were obtained utilizing the neighborhood cultivation genetic algorithm (NCGA). The ﬁndings revealed that the performance models are primarily a ﬀ ected by the feed rate, depth of cut, speed, and inclined angle, respectively. The optimal values of the α , a p , f , and v c are 26 ◦ , 0.44 mm, 0.37 mm / rev, and 200 mm / min, respectively. The improvements in energy e ﬃ ciency, average roughness, and cost are 8.91%, 20.00%, and 14.75%, as compared to the initial values. Moreover, the NCGA may perform an e ﬃ cient operation to obtain the optimal outcomes, as compared to conventional algorithms. conducted, power components are automatically recorded. The experimental data is stored on a ﬂashcard and analyzed with the aid of the KEW6305 software. An interval of 0.1 sec. is used to improve the accuracy of the measured data. A roughness tester labeled Mitutoyo SJ-301 is used to measure the roughness in ﬁve investigated points. The sampling length of 4 mm is used to ensure the measuring precision.


Introduction
The rotary turning is an alternative solution to machine hardened metals, in which the round insert rotates continuously about its axis in the cutting period. The rotary turning can be divided into two primary methods, including the actively driven (ADRT) and self-propelled rotary tool (SPRT) cutting operations. In the ADRT machining, the rotational speed of the insert is controlled using an external motor, while the rotary motion is produced by the interaction with the workpiece in the SPRT turning. The rotary motion of the insert provides the rest period, which results in minimizing the cutting temperatures, longer tool life, and higher production rate [1]. Moreover, the fresh insert in the turning time has a significant contribution to increased surface quality. The rotary turning operation possesses various advantages, including cost-effectiveness, less machining time, and ensuring machined quality, as compared to the traditional approaches, such as grinding and polishing operations. the tool wear, heat generation, surface roughness, and cutting force were explored. The experiments were conducted to validate the simulation results. The authors stated that the roughness and turning force were decreased by 50%, as compared to the conventional method. The heat concentration in the contact zone was suppressed with the aid of the rotational motion. The UART operation had a significant contribution to a reduction in the tool wear. The impacts of the oil mist on the turning temperature and tool wear in the ADRT turning of the stainless steel at various cutting speeds were explored by Shasahara et al. [12]. The findings indicated that the oil mist helped to decrease the cutting temperature and flank wear. The cutting temperature and tool wear of the ADRT turning were lower than that of the fixed tool. The range of the cutting speed could be increased in the rotary turning.
A novel method was developed to predict the machining position and feed direction-dependent tool point frequency response functions (FRFs) based on the modal theory, matrix transformation, the backpropagation neural network, and particle swarm optimization algorithm. The feasibility of the proposed approach was checked by case studies on a vertical machining center. The authors stated that the developed approach can contribute to select the optimal process planning based on actual milling conditions at the design stage [13]. Additionally, a reacceptance coupling-based approach was developed to describe the dynamic properties of the machine-tool-workpiece system [14]. The experimental validation was used to check the fidelity of the proposed procedure. The outcomes revealed that the knowledge could be employed to select the technological parameter, such as the spindle speed for the computerized numerical control lathe. The author stated that the proposed method could be applied in both old and new machine tools to avoid expensive hardware implementations. Lin et al. developed the predictive models of the surface roughness and machining vibration in terms of the cutting parameters using the multiple regression analysis (MRA) and artificial neural network (ANN) methods [15]. The authors concluded that the ANN model could provide higher prediction accuracy, as compared to the MRA one. The developed correlations could be applied in an online roughness monitoring system for the surface milling test. The proposed models considering the machining stability could be used to enhance the machining accuracy in milling processes.
As a result, the cutting speed, depth, feed, inclined angle, and lubrication characteristics are considered as the optimized parameters. The machining performances are turning forces, roughness, tool wear, cutting temperature, and the formation of the chip were analyzed. Different optimization approaches, such as RSM, neural network, and simulation method were employed to propose the performance models and select the optimal factors. However, the selection of appropriate parameters to enhance energy efficiency and reduce the roughness, as well as costs, has not been presented in the aforementioned works. To achieve sustainable production, the optimal varied factors have to be selected to satisfy the environmental, economic, and social aspects.
To overcome the analyzed gaps, a parameter-based optimization of the rotary turning using SPRT has been investigated and presented to improve the economic, social, and environmental performances. The energy efficiency (EF R ) is an important indicator of the environmental impact. The average roughness (R a ) and turning cost (CT) are the primary indicators of social and economic metrics. In this work, the self-propelled rotary tool has been designed and fabricated. The improved Kriging models are applied to render nonlinear relationships between the process parameters and turning responses [16,17]. The neighborhood cultivation genetic algorithm (NCGA) is employed to select the optimal values of the inputs and objectives [18].

Objective Functions
In this work, four turning responses, including energy efficiency, average roughness, turning cost, and the index of the personal health and operational safety are optimized simultaneously. The energy efficiency of the rotary turning (EF R ) is defined as the ratio of the turning energy consumed (TEC) to the machined energy consumed (MEC) during turning. The value of the EF R is calculated using Equation (1): where: ATP-the active turning power; MPC-the machine power consumed; t c -the turning time.
The average roughness (R a ) is calculated as: where: R ai -the value of average roughness at the ith position; n-the number of measured points. The turning cost model under dry condition (CT) is expressed as follow: where: C TM -the turning cost; C E -the energy cost; C CH -the tool changing cost; C T -the tool costs. The total rotary turning time is calculated as: where: t o -the start-up time; t st -the set-up time; t air -the air-turning; t ch -the tool change time.
The turning cost (C TM ), the energy cost (C E ), the tool changing cost (C CH ), and the tool cost (C T ) are calculated using Equations (5)-(9): C E = x e (P o t o + P st t st + P air t air + P st t ch t T T + P T t T ) = x e (P o t o + P st t st + P air t air + P st t ch × 3.9 × 10 −6 V + [P air + kMRR] π×D×L×60 1000×v c × f ) where: x-the total cost for labor, machine, and management (VND/s); x e -the energy cost (VND/kWh); y c -the tooling cost (VND/piece); MRR-the material removal rate (mm 3 /min); k-the specific energy requirement in machining operations (kJ/cm 3 ); P o -the consumed power in the start-up time; P st -the consumed power in the set-up time; P air -the consumed power in the air-turning; P ch -the consumed power in the tool change time; T-the tool life; c-the coefficient related to the turning conditions; α and β-the positive constant parameters depending on tool material and workpieces. The index of the personal health and operational safety (POS) can be divided into five primary components, which are expressed as follows: 1.
The operational safety index (POS tc ) is used to exhibit the toxic chemicals exposure during cutting fluid preparation or cutting process itself. The value of the POS tc can be 1, 2, or 3 depending on its harmful impacts due to the type and concentration. 2. The operational safety index (POS hts ) is used to exhibit the high-temperature surface exposure during the cutting processes. The value of the POS hts is expressed as: where T o denotes the cutting temperature. 3.
The operational safety index (POS hss ) is used to exhibit the high-speed surface exposure during the cutting processes. The value of the POS hss is expressed as: where n max presents the maximum spindle speed. 4.
The personal health index (POS wes ) is used to exhibit the working environment noise level during the cutting processes. The value of the POS wes is expressed as: where NI presents the measured noise level. 5.
The personal health index (POS wea ) is used to exhibit the working environment atmospheric conditions during the cutting processes. The value of the POS wea is expressed as: where WBGT presents the wet-bulb globe temperature. 6.
The personal health index (POS wei ) is used to exhibit the working environment illumination level during the cutting processes. The value of the POS wei is expressed as: where IL presents the measured illumination level.
Therefore, the total index of personal health and operational safety is computed as:

Varied Parameters of the Rotary Turning
For the rotary turning, the inputs considered are the process parameters (depth of cut, inclined angle, speed, and feed rate), the insert characteristics (geometry, dimensions, and materials), lubrication properties (type and lubrication method), and the workpiece materials.
In the current work, the characteristics of the turning insert and machining specimen are considered as constants. All experiments are performed under the dry condition; hence, the influences of the lubrication are neglected. Therefore, the inclined angle, cutting speed, depth of cut, feed rate, and nose radius are listed as the primary variables that intensively influence the machining performances of the rotary turning.
Four key factors having the ranges, including the inclined angle, cutting speed, depth of cut, and feed rate are exhibited in Table 1. The parameter levels are identified based on the recommendations of the insert manufacturer, the characteristics of the fabricated rotary tool, and material properties. These values are validated using the mechanical handbooks and aforementioned publications.

Optimization Framework
Sustainable production has become the top priority for industrial manufacturers due to an increase in energy costs and environmental legislation. Three primary aspects of sustainable production in terms of the environment, economy, and society are depicted in Figure 1. The environmental aspects of the turning process are lower power consumed (maximizing energy efficiency or minimizing carbon emissions), higher power factor, reduced tool as well as workpiece consumption, and elimination of the hazardous materials. The economic aspects are achieved through the reduction of production costs, higher productivity, and reduced lead time. The social aspects include higher surface quality, lower-dimensional deviation, and operational safety. Different cooling-lubrication techniques [19,20] were reported to exhibit superior performance and cleaner machining system with enhanced productivity, product quality and temperature. However, that was not enough to promote complete sustainability which inspired the optimization of the machining system [21].

Optimization Framework
Sustainable production has become the top priority for industrial manufacturers due to an increase in energy costs and environmental legislation. Three primary aspects of sustainable production in terms of the environment, economy, and society are depicted in Figure 1. The environmental aspects of the turning process are lower power consumed (maximizing energy efficiency or minimizing carbon emissions), higher power factor, reduced tool as well as workpiece consumption, and elimination of the hazardous materials. The economic aspects are achieved through the reduction of production costs, higher productivity, and reduced lead time. The social aspects include higher surface quality, lower-dimensional deviation, and operational safety. Different cooling-lubrication techniques [19,20] were reported to exhibit superior performance and cleaner machining system with enhanced productivity, product quality and temperature. However, that was not enough to promote complete sustainability which inspired the optimization of the machining system [21]. The research approach depicted in Figure 2 includes the following steps.
Step 2: The predictive models of the EFR and Ra are then established in terms of the inputs using the improved Kriging models. The Kriging model is an interpolation-based correlation using the Gaussian function, which is expressed as a function of the global model u and local deviation Z(x) [25]: The distance between the x i and x j is calculated as: where θk is the correlative factor of the kth element. The correlation between the Z(x i ) and Z(x j ) is calculated as: The Kriging predictor at a specific value of x is calculated by: The research approach depicted in Figure 2 includes the following steps.
Step 2: The predictive models of the EF R and R a are then established in terms of the inputs using the improved Kriging models. The Kriging model is an interpolation-based correlation using the Gaussian function, which is expressed as a function of the global model u and local deviation Z(x) [25]: The distance between the x i and x j is calculated as: where θ k is the correlative factor of the kth element.
The correlation between the Z(x i ) and Z(x j ) is calculated as: The Kriging predictor at a specific value of x is calculated by: r is the vector of the ith element.
The correlation factor is estimated as: The values of theμ andσ 2 at the specific θ k are defined as: In order to improve the predictive accuracy, the additional points distributed in the interesting region (e.g., Pareto set) are used to learn and rebuild the Kriging models. The additional data for generating the Kriging models are selected by the possibility of existing in this region. The convergence tolerance ε without considering the magnitudes of outputs is determined using the following constraint: where MAE r is the relatively max absolute error. y i andŷ i present the experimental and predicted values. The adequacy of the developed models can be evaluated using the normalized root mean squared error (NRSME). The value of NRSME is calculated as: Step 3: An evolutionary algorithm labeled neighborhood cultivation genetic algorithm (NCGA) is employed to identify the optimal values of process parameters and responses. The NCGA is an efficient technique to select the optimal solution in the highly non-linear design space. The primary principle of the NCGA is to add a neighborhood crossover operation into the working mechanism. The together individuals are sorted into groups generated by the population. The genetic operations are performed in a group that is consisted of the two individuals. The operating algorithm of the NCGA is listed as bellows ( Figure 3  -Termination: If the terminal condition is satisfied, the simulation is terminated. If it is not satisfied, the simulation returns to step 2.

Experimental Details
The steel bars labeled 40XC (diameter 40 mm and length 260 mm) are used as specimens in this work. The hardness of the hardened workpiece is around 45 HRC. The material is chosen due to extensive applications in the automotive, molding, and marine industries. The workpiece is divided into different machining segments. Each workpiece is machined in a length of 26 mm. Its chemical composition is shown in Table 2. Table 2. Chemical compositions of hardened 40XC [26].

Experimental Details
The steel bars labeled 40XC (diameter 40 mm and length 260 mm) are used as specimens in this work. The hardness of the hardened workpiece is around 45 HRC. The material is chosen due to extensive applications in the automotive, molding, and marine industries. The workpiece is divided into different machining segments. Each workpiece is machined in a length of 26 mm. Its chemical composition is shown in Table 2.
The turning trials are conducted under dry conditions using the TiAN coated rotary insert (RHHW 1605MOTN IN 2005) with 15° rake angle provided by INGERSOLL (Figure 4). The outside diameter, inside diameter, and thickness of the round insert are 16, 5.2, and 5.2 mm, respectively. The hardness of the employed insert is around 78 HRC. The mold steel labeled 2311 is used to fabricate the tool shank.
All the turning trials are performed on a CNC lathe labeled CTX 400E manufactured by GILDEMEISTER (Bielefeld, Germany). The controller namely HEIDENHAIN CNC-PILOT is employed to perform all machining operations. The machine used had a power spindle of 15 kW and a maximum of 4500 rpm. The turning experiments are repeated two times to ensure the data repeatability and a new insert is employed in each run. A power meter labeled KEW6305 is employed to capture power components during the rotary turning. Three clamp sensors are connected to the three-phase power lines of CNC lathe with correct direction. The power meter is turned on after the appropriate connection. When the rotary turning is conducted, power components are automatically recorded. The experimental data is stored on a flashcard and analyzed with the aid of the KEW6305 software. An interval of 0.1 sec. is used to improve the accuracy of the measured data. A roughness tester labeled Mitutoyo SJ-301 is used to measure the roughness in five investigated points. The sampling length of 4 mm is used to ensure the measuring precision.  The obtained values of total power consumption and average roughness are exhibited in Figure 5. All the turning trials are performed on a CNC lathe labeled CTX 400E manufactured by GILDEMEISTER (Bielefeld, Germany). The controller namely HEIDENHAIN CNC-PILOT is employed to perform all machining operations. The machine used had a power spindle of 15 kW and a maximum of 4500 rpm. The turning experiments are repeated two times to ensure the data repeatability and a new insert is employed in each run.
A power meter labeled KEW6305 is employed to capture power components during the rotary turning. Three clamp sensors are connected to the three-phase power lines of CNC lathe with correct direction. The power meter is turned on after the appropriate connection. When the rotary turning is conducted, power components are automatically recorded. The experimental data is stored on a flashcard and analyzed with the aid of the KEW6305 software. An interval of 0.1 sec. is used to improve the accuracy of the measured data. A roughness tester labeled Mitutoyo SJ-301 is used to measure the roughness in five investigated points. The sampling length of 4 mm is used to ensure the measuring precision.
The obtained values of total power consumption and average roughness are exhibited in Figure 5.

Development of Improved Kriging Models
The coefficients for the cost models are exhibited in Table 3. Experimental results are shown in Table 4. The experimental data from 1 to 27 were used to construct improved Kriging models. The trials from 28 to 33 were adopted in order to explore the precision of the models proposed.

Development of Improved Kriging Models
The coefficients for the cost models are exhibited in Table 3. Experimental results are shown in Table 4. The experimental data from 1 to 27 were used to construct improved Kriging models. The trials from 28 to 33 were adopted in order to explore the precision of the models proposed.  The normalized root mean squared error (NRMSE) is used to evaluate the soundness of the improved Kriging models. The NRMSE values of the EF R and R a are 0.0018 and 0.0011, respectively ( Figure 6). It can be stated that the fidelity of the improved Kriging models is acceptable. ANOVA analysis having 95% confidence is employed to check the adequacy of the developed models [27]. The factors having p-value less than 0.05 are significant. ANOVA results of the energy efficiency model are shown in Table 5. The R 2 value of 0.9812 indicated that the energy efficiency model presented 98.12% of the total variety of the experimental data. The adjusted R 2 value of 0.9632 indicated that the proposed model exhibited a variation of 96.32% with significant factors. Moreover, the predicted R 2 value of 0.9526 indicated that the proposed model exhibited a variation of 9526% with any new data. The significant factors are the single terms (α, ap, f, and vc), interaction term (αap and apvc), and quadratic terms (ap 2 , f 2 , α 2 , and vc 2 ). Other terms having p-value higher than 0.05 are insignificant factors. ANOVA results of the average roughness model are shown in Table 6. The R 2 value of 0.9834 indicated that the average roughness model presented 98.34% of the total variety of the ANOVA analysis having 95% confidence is employed to check the adequacy of the developed models [27]. The factors having p-value less than 0.05 are significant. ANOVA results of the energy efficiency model are shown in Table 5. The R 2 value of 0.9812 indicated that the energy efficiency model presented 98.12% of the total variety of the experimental data. The adjusted R 2 value of 0.9632 indicated that the proposed model exhibited a variation of 96.32% with significant factors. Moreover, the predicted R 2 value of 0.9526 indicated that the proposed model exhibited a variation of 9526% with any new data. The significant factors are the single terms (α, a p , f, and v c ), interaction term (αa p and a p v c ), and quadratic terms (a p 2 , f 2 , α 2 , and v c 2 ). Other terms having p-value higher than 0.05 are insignificant factors. ANOVA results of the average roughness model are shown in Table 6. The R 2 value of 0.9834 indicated that the average roughness model presented 98.34% of the total variety of the experimental data. The adjusted R 2 value of 0.9664 indicated that the proposed model exhibited a variation of 96.64% with significant factors. Moreover, the predicted R 2 value of 0.9552 indicated that the proposed model exhibited a variation of 95.52% with any new data. The significant factors are the single terms (α, a p , f, and v c ), interaction term (a p v c ), and quadratic terms (a p 2 , f 2 , α 2 , and v c 2 ). Other terms having p-value higher than 0.05 are insignificant factors. To confirm the accuracy of the developed models, the experimental outcomes are compared with the predictive values. Figure 7 depicts the comparisons between the predicted and experimental values. The small errors indicated that the models have high consistency and the models can be employed to predict the response values.
Metals 2020, 10, x FOR PEER REVIEW 13 of 22 experimental data. The adjusted R 2 value of 0.9664 indicated that the proposed model exhibited a variation of 96.64% with significant factors. Moreover, the predicted R 2 value of 0.9552 indicated that the proposed model exhibited a variation of 95.52% with any new data. The significant factors are the single terms (α, ap, f, and vc), interaction term (apvc), and quadratic terms (ap 2 , f 2 , α 2 , and vc 2 ). Other terms having p-value higher than 0.05 are insignificant factors. To confirm the accuracy of the developed models, the experimental outcomes are compared with the predictive values. Figure 7 depicts the comparisons between the predicted and experimental values. The small errors indicated that the models have high consistency and the models can be employed to predict the response values.

Parametric Effects
Energy efficiency is a significant indicator to evaluate the environmental impact of the dry turning process and high energy efficiency is desirable. The impacts of the varied factors on the energy efficiency of the rotary turning process are shown in Figure 8.

Parametric Effects
Energy efficiency is a significant indicator to evaluate the environmental impact of the dry turning process and high energy efficiency is desirable. The impacts of the varied factors on the energy efficiency of the rotary turning process are shown in Figure 8.  Figure 8a exhibits the effects of the depth of cut and inclined angle on energy efficiency. It can be stated that higher efficiency is associated with an increment in the depth and/or inclined angle. Additionally, a further value of the inclined angle leads to a reduction in energy efficiency. The phenomena can be explained as follows.
The thickness of the chip increases with an increment in the depth of cut, which causes higher contact region between the rotary insert and workpiece. A higher amount of the material is turned, which leads to an increase in the degree of plastic deformation. Higher frictional resistance is produced. As a result, a higher active power is required to overcome friction, and energy efficiency increases [28]. An increase in the inclined angle leads to a lower contact area between the workpiece and the turning insert. A lower value of the total power consumption is employed to turn the material, hence, energy efficiency is improved. In contrast, a further value of the inclined angle causes a higher turning area due to the perpendicular between the rotary insert and the workpiece. A higher amount of the material is processed and a higher degree of plastic deformation is obtained. Greater resistance causes higher the total power consumed in the CNC lathe, hence, energy efficiency decreases. Consequently, the highest depth of cut can be employed to enhance energy efficiency, while the middle angle is recommended to use. Figure 8b exhibits the effects of the feed rate and speed on energy efficiency. It can be stated that increased speed and/or feed leads to an increment in energy efficiency. A higher power of the spindle motor is required with an increment in the speed. The turning momentum also increases to  Figure 8a exhibits the effects of the depth of cut and inclined angle on energy efficiency. It can be stated that higher efficiency is associated with an increment in the depth and/or inclined angle. Additionally, a further value of the inclined angle leads to a reduction in energy efficiency. The phenomena can be explained as follows.
The thickness of the chip increases with an increment in the depth of cut, which causes higher contact region between the rotary insert and workpiece. A higher amount of the material is turned, which leads to an increase in the degree of plastic deformation. Higher frictional resistance is produced. As a result, a higher active power is required to overcome friction, and energy efficiency increases [28]. An increase in the inclined angle leads to a lower contact area between the workpiece and the turning insert. A lower value of the total power consumption is employed to turn the material, hence, energy efficiency is improved. In contrast, a further value of the inclined angle causes a higher turning area due to the perpendicular between the rotary insert and the workpiece. A higher amount of the material is processed and a higher degree of plastic deformation is obtained. Greater resistance causes higher the total power consumed in the CNC lathe, hence, energy efficiency decreases. Consequently, the highest depth of cut can be employed to enhance energy efficiency, while the middle angle is recommended to use. Figure 8b exhibits the effects of the feed rate and speed on energy efficiency. It can be stated that increased speed and/or feed leads to an increment in energy efficiency. A higher power of the spindle motor is required with an increment in the speed. The turning momentum also increases to satisfy a higher setting of the cutting speed. As a result, higher active turning power is obtained, hence, energy efficiency improves. Additionally, the strength and hardness of the machined specimen may decrease due to an increment in the temperature at the turning region, which causes a softer material. A lower value of the total power consumed is then obtained and energy efficiency increases. According to Figure 8b, when the feed rate increases, the reaction force on the spindle system increases. Therefore, the active cutting power of the spindle system increases. Additionally, a higher value of the feed rate leads to an increment in the machining temperature and may cause the work-hardening on the workpiece surface. This results in an increase in the turning force due to a greater resistance, hence, higher active cutting power is consumed. Moreover, an increment in the feed causes higher contact region between the tool and workpiece. Obviously, higher active cutting power is required to process material and remove the chip. Therefore, increased feed causes higher energy efficiency.
The parametric contributions of the varied factors on the energy efficiency model are exhibited in Figure 8c. The blue and red bars denote the negative and positive impacts, respectively. The feed rate is the most effective factor (19.87%), followed by the depth (19.57%), speed (18.05%), and inclined angle (3.93%), respectively. The contribution of the α 2 , a p 2 , f 2 , and v c 2 are 24.80%, 7.90%, 2.25%, and 0.88%, respectively. The surface roughness is an important evaluator to explore the social impact of the dry turning process and a low roughness is preferred. Figure 9 depicted the influences of the varied parameters on the average roughness.
Metals 2020, 10, x FOR PEER REVIEW 15 of 22 satisfy a higher setting of the cutting speed. As a result, higher active turning power is obtained, hence, energy efficiency improves. Additionally, the strength and hardness of the machined specimen may decrease due to an increment in the temperature at the turning region, which causes a softer material. A lower value of the total power consumed is then obtained and energy efficiency increases. According to Figure 8b, when the feed rate increases, the reaction force on the spindle system increases. Therefore, the active cutting power of the spindle system increases. Additionally, a higher value of the feed rate leads to an increment in the machining temperature and may cause the work-hardening on the workpiece surface. This results in an increase in the turning force due to a greater resistance, hence, higher active cutting power is consumed. Moreover, an increment in the feed causes higher contact region between the tool and workpiece. Obviously, higher active cutting power is required to process material and remove the chip. Therefore, increased feed causes higher energy efficiency.
The parametric contributions of the varied factors on the energy efficiency model are exhibited in Figure 8c. The blue and red bars denote the negative and positive impacts, respectively. The feed rate is the most effective factor (19.87%), followed by the depth (19.57%), speed (18.05%), and inclined angle (3.93%), respectively. The contribution of the α 2 , ap 2 , f 2 , and vc 2 are 24.80%, 7.90%, 2.25%, and 0.88%, respectively.
The surface roughness is an important evaluator to explore the social impact of the dry turning process and a low roughness is preferred. Figure 9 depicted the influences of the varied parameters on the average roughness.   Figure 9a exhibits the effects of the depth of cut and inclined angle on the roughness. It can be stated that higher roughness is associated with an increment in the depth. A reduction in the roughness is associated with an increased inclined angle. Unfortunately, a further value of the inclined angle leads to an increase in the roughness. The turning area between the tool and the specimen is increased with a higher depth of cut. The tuned chip is heavily removed from the machined surface, which causes an increase in the roughness. The turned material volume is decreased with a larger inclined angle and material is easily processed. Lower roughness is obtained [29]. When the inclined angle increases from the middle to the highest levels, the turned material volume is increased due to the perpendicular between the cutting tool and the workpiece. A higher amount of material is processed and higher roughness is produced.
The impacts of the feed rate and cutting speed on the average roughness are shown in Figure 9b. The roughness is significantly decreased at the highest value of the cutting speed. Additionally, a higher feed rate leads to an increase in the roughness. The phenomena can be explained as follows. Higher speed results in increased temperature at the turning region, which causes a reduction in the hardness and strength of the workpiece. The chip is smoothly removed from the turned surface, hence, roughness is decreased. Low feed rate leads to a small distance between the turned peaks, which causes a low roughness. Higher distance between the turned peaks is produced at higher feed, hence, higher feed mark is generated on the machined surface.
The parametric contributions of the varied factors on the roughness model are exhibited in Figure 9c. The feed rate is the most effective factor (19.66%), followed by the speed (15.97%), depth (14.70%), and inclined angle (9.44%), respectively. The contribution of the α 2 , a p 2 , v c 2 , and f 2 are 18.12%, 7.41%, 5.50%, and 3.41%, respectively. The turned images at different process inputs are depicted in Figure 10. As a result, a low roughness is obtained at low values of the depth of cut as well as feed and the middle value of the inclined angle, in which the smooth cut, small grooves, and small pits are produced (Figure 10a). A higher roughness is produced at high values of the depth of cut as well as feed rate and the lowest inclined angle, in which the bigger grooves and larger pits as well as voids are observed (Figure 10b). roughness is associated with an increased inclined angle. Unfortunately, a further value of the inclined angle leads to an increase in the roughness. The turning area between the tool and the specimen is increased with a higher depth of cut. The tuned chip is heavily removed from the machined surface, which causes an increase in the roughness. The turned material volume is decreased with a larger inclined angle and material is easily processed. Lower roughness is obtained [29]. When the inclined angle increases from the middle to the highest levels, the turned material volume is increased due to the perpendicular between the cutting tool and the workpiece. A higher amount of material is processed and higher roughness is produced.
The impacts of the feed rate and cutting speed on the average roughness are shown in Figure  9b. The roughness is significantly decreased at the highest value of the cutting speed. Additionally, a higher feed rate leads to an increase in the roughness. The phenomena can be explained as follows. Higher speed results in increased temperature at the turning region, which causes a reduction in the hardness and strength of the workpiece. The chip is smoothly removed from the turned surface, hence, roughness is decreased. Low feed rate leads to a small distance between the turned peaks, which causes a low roughness. Higher distance between the turned peaks is produced at higher feed, hence, higher feed mark is generated on the machined surface.
The turned images at different process inputs are depicted in Figure 10. As a result, a low roughness is obtained at low values of the depth of cut as well as feed and the middle value of the inclined angle, in which the smooth cut, small grooves, and small pits are produced (Figure 10a). A higher roughness is produced at high values of the depth of cut as well as feed rate and the lowest inclined angle, in which the bigger grooves and larger pits as well as voids are observed ( Figure  10b).
The SEM images of the round insert at various machining conditions are depicted in Figure 11. As depicted in Figure 11a, the wear is not much prominent without grooves or fracture at the low values of the depth of cut and feed rate. An increment in the depth and/or feed leads to bigger grooves and fractures on the rake face. Higher machining conditions may cause excessive pressure and stress, which decreases the hardness of the rotary insert (Figure 11b).  The SEM images of the round insert at various machining conditions are depicted in Figure 11. As depicted in Figure 11a, the wear is not much prominent without grooves or fracture at the low values of the depth of cut and feed rate. An increment in the depth and/or feed leads to bigger grooves and fractures on the rake face. Higher machining conditions may cause excessive pressure and stress, which decreases the hardness of the rotary insert (Figure 11b).

Optimization Results
In the current work, the proposed models can be adopted to obtain the optimal inputs for improvements in energy efficiency, roughness, and turning cost. Table 7 presents the setting parameter of the NCGA algorithm. The 3D plots generated by NCGA are depicted in Figure 12. When higher operator safety is selected, the POS of 6 is chosen. On the contrary, the POS of 5 is selected with a lower impact on operational safety. In this study case, a CNC late having a protecting cover is used to separate the operator from the turning process, hence, the lower curve (POS = 5) should be selected.
The Pareto fronts are shown in Figure 13. As a result, it is difficult to select the optimal solution, which can satisfy the optimization requirement. An increment in the energy efficiency leads to an increase in the average roughness. To decrease the average roughness, higher machining cost is required. Various optimal solutions are exhibited on the Pareto fronts, which can be used to satisfy different turning purposes. Three typical points, including points 1, 2, and 3 are selected to evaluate the turning performances and to compare with the initial values. As shown in Table 8, the average roughness increases at the first solution, which does not satisfy the optimization requirement. For the second point, the energy efficiency decreases, which cannot consider as a technical solution. For the third solution, the roughness and cost decrease while the energy efficiency improves. Consequently, the third point can be selected as a potential solution to enhance turning performances.

Optimization Results
In the current work, the proposed models can be adopted to obtain the optimal inputs for improvements in energy efficiency, roughness, and turning cost. Table 7 presents the setting parameter of the NCGA algorithm. The 3D plots generated by NCGA are depicted in Figure 12. When higher operator safety is selected, the POS of 6 is chosen. On the contrary, the POS of 5 is selected with a lower impact on operational safety. In this study case, a CNC late having a protecting cover is used to separate the operator from the turning process, hence, the lower curve (POS = 5) should be selected.
The Pareto fronts are shown in Figure 13. As a result, it is difficult to select the optimal solution, which can satisfy the optimization requirement. An increment in the energy efficiency leads to an increase in the average roughness. To decrease the average roughness, higher machining cost is required. Various optimal solutions are exhibited on the Pareto fronts, which can be used to satisfy different turning purposes. Three typical points, including points 1, 2, and 3 are selected to evaluate the turning performances and to compare with the initial values. As shown in Table 8, the average roughness increases at the first solution, which does not satisfy the optimization requirement. For the second point, the energy efficiency decreases, which cannot consider as a technical solution. For the third solution, the roughness and cost decrease while the energy efficiency improves. Consequently, the third point can be selected as a potential solution to enhance turning performances. (a) Average roughness and energy efficiency (a) Average roughness and energy efficiency  To evaluate the effectiveness of the NCGA, common optimization algorithms, including the Non dominated sorting genetic algorithm II (NSGA II) and Multiple objective particle swarm optimization (MOPSO) are adopted to select the optimal varied factors. As shown in Table 9, the optimal results generated by the NCGA are better than the MOPOSO and NSGA II. It can be stated that the NCGA can provide reliable outcomes when dealing with multi-objective optimization, as compared to the traditional ones. The optimal values of the inclined angle, depth, feed, and speed are 26°, 0.44 mm, 0.37 mm/rev, and 200 m/min, respectively. Consequently, the improvements in energy efficiency, average roughness, and turning cots are 16.11% and 24.35%, and 17.45%, respectively.  To evaluate the effectiveness of the NCGA, common optimization algorithms, including the Non dominated sorting genetic algorithm II (NSGA II) and Multiple objective particle swarm optimization (MOPSO) are adopted to select the optimal varied factors. As shown in Table 9, the optimal results generated by the NCGA are better than the MOPOSO and NSGA II. It can be stated that the NCGA can provide reliable outcomes when dealing with multi-objective optimization, as compared to the traditional ones. The optimal values of the inclined angle, depth, feed, and speed are 26 • , 0.44 mm, 0.37 mm/rev, and 200 m/min, respectively. Consequently, the improvements in energy efficiency, average roughness, and turning cots are 16.11% and 24.35%, and 17.45%, respectively.

Conclusions
This paper presented a machining parameter-based optimization for the rotary turning operation to achieve the improvements in energy efficiency, average roughness, turning costs, and operational safety. The predictive correlations were proposed in terms of the depth of cut, inclined angle, feed rate, and cutting speed with the aid of the improved Kriging models. ANOVA results were employed to identify the significances of the proposed models and parameters. The NCGA was used to obtain the optimal factors. The main conclusions are listed as follows: