Next Article in Journal
A New Parameter Identification Method for Industrial Robots with Friction
Previous Article in Journal
Fault Detection for Interval Type-2 T-S Fuzzy Networked Systems via Event-Triggered Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Energy Efficiency Tool Path Optimization Method Using a Discrete Energy Consumption Path Model

1
State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou 310027, China
2
Hangzhou Innovation Institute, Beihang University, Hangzhou 310052, China
3
Wuxi Lead Intelligent Equipment Co., Ltd., Wuxi 214125, China
4
School of Mechanical Engineering, Shenyang University of Technology, Shenyang 110023, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(5), 348; https://doi.org/10.3390/machines10050348
Submission received: 18 April 2022 / Revised: 5 May 2022 / Accepted: 6 May 2022 / Published: 8 May 2022
(This article belongs to the Section Machine Design and Theory)

Abstract

:
As the energy cost accounts for about one-third of the total manufacturing cost, there is great significance in evaluating and managing energy consumption in manufacturing processes. The energy consumption during multi-axis end milling, which represents a large part of the industrial energy costs, is usually extraordinarily large, especially for complex free-form surfaces requiring multi-finish-machining. To obtain the most efficient tool path, the tool orientation is adjusted to obtain the largest cutting stripe width at each cutter contact point. However, the use of excessive driving energy consumption and cutting energy to obtain the largest cutting stripe width may reduce the energy efficiency of the tool path. To solve this problem, the geometry features of the tool path are analyzed firstly, and the global energy consumption analysis, which includes a cutting energy analysis and driving energy analysis, is conducted. The discrete energy consumption path model is constructed to find the most energy-efficient tool orientation sequence for a tool path. Finally, contrast experiments are carried out to validate the proposed method.

1. Introduction

Due to the vast energy costs involved in manufacturing, improving the energy efficiency is a good approach to improving energy savings. Multi-axis end milling is widely used in machining free-form surfaces, for which high levels of machining accuracy and efficiency are required. In aviation, aerospace, shipbuilding, and nuclear industries, complex parts such as blisks and blades are machined using high speed multi-axis end milling, especially in semi-finish and finish machining operations. The high hardness and large size of the machined parts usually lead to long machining times. Considering the high power of multi-axis machine tools, a significant amount of electric energy is consumed as a result. For example, according to Edem [1], more than 20,000 euros worth of electricity is consumed in five-axis machining, and the energy costs account for one-third of the total machining costs (including operation, maintenance, and depreciation). Molds and dies are changed according to the customers’ requirements, but blisks and blades will generally remain the same for a long period. Thus, even a small energy reduction in such a machining process would mean a huge energy-efficiency promotion in the long term.
Multi-axis end milling is a complicated process involving three independent phases. At the computer-aided manufacturing (CAM) phase, tool paths are generated according to the properties of the machining materials, tool cutter, and the machine tool configuration. Specific machining strategies are used in work piece coordinate systems (WCSs). Afterward, at the data conversion and numerical control (CNC) phase, the geometric data for the tool path is converted into a G-code, which can be executed directly by the machine tool according to inverse kinematics properties of the machine tool. Finally, at the execution stage, the G-code program is executed and the energy consumption can be estimated. This complete machining process contains a series of inter-related machining parameters. Finding the optimal parameter combination and optimizing energy consumption in the machining process draw lots of attention. A comprehensive review of energy consumption models in five-axis machining was made by Zhou et al. [2]. It was noticed that the material removal rate (MRR) is the dominant factor that determines the energy efficiency of the machining process, and the MRR can be impacted by a series of machining parameters, including the cutting force, feed rate, and scallop height requirement. Energy consumption data were studied by Aramcharoen et al. [3], and it was noticed that there was a linear correlation between the feed rate and driving power. Thus, the spiral tool path, especially the changing of the orientation of the tool, will improve the energy efficiency due to the elimination of unnecessary air cuts. Additionally, the feed rate and other machining parameters are highly coupled and determine the energy cost of the machine tool. The characteristics of energy consumption were studied using the systematic analysis method considering the material properties of the workpiece. S/N analysis was applied to investigate the relations between process parameters and energy consumption [4]. Detailed explanations have been given to illustrate the different optimization trends for two kinds of specific energy consumption. Some equipment for CNC systems has been developed to help to improve the energy efficiency of the machining process. For example, the conventional numerical control strategy of an energy-efficient stepper drive was demonstrated by Arzhanov et al. [5]. Considering the configuration properties of a feed drive system, including the masses, sizes, and moments of inertia, an analytical model optimizing the driving energy consumption machining process was constructed by Ryuta et al. [6]. From the perspective of data mining, ANN has been utilized to construct a predictive model of energy consumption, where the spindle speed, feed rate, cutting depth, and cutting stripe width are optimized simultaneously [7]. Additionally, the energy consumption of the whole manufacturing process, which involves milling, turning, and grinding, can be optimized using the modified heuristic algorithms [8,9,10,11,12].
The aforementioned studies on machining energy consumption tend to optimize the tool path at the second or third phase of the CNC machining process. For multi-axis end milling, for which the machining quality is the dominant requirement, the adjustable parameters of the last two stages are the spindle speed and feed rate. However, there is no room to adjust these parameters due to a series of constraints involved, such as machine tool kinematics and cutting condition requirements. Therefore, this paper focuses on the most dominant factor of a multi-axis end milling process, i.e., the tool path. Plenty of optimization methods have been proposed in the last two decades. As the machined stripe width can change along with the equivalent curvature of the cutter at each CC point, Marciniak et al. [12] proposed a tool path generation method that could obtain the maximum cutting stripe width by matching the principal curvatures of the work piece and the cutter. However, unforeseen overcuts and undercuts may occur when the tool orientation changes abruptly, meaning the machining quality will be reduced due to the scratching on the machined surface. In addition to obtaining the maximum cutting stripe width, the elimination of sudden changes in neighboring cutter postures also deserves researchers’ attention. Therefore, smoothing methods for the tool path, which involves a series of CC points, have been studied by many researchers. The C-space model [13], which was first used for robot joint movement planning, could significantly help minimize the movements of the feed drive systems of machine tools. The same movements of the cutter, including translational movements and rotational movements, can be executed in different ways due to the specific kinematic configuration of each machine tool [14]. Cong et al. [15] employed a modified interpolation algorithm to simultaneously optimize the tool path of the tool cutter and for interpolation of tool orientation vectors. As the energy is consumed by the machine tool, it can be concluded that optimizing the movements in the MCS is more efficient than optimizing the tool orientation in the WCS. Algorithm studies were conducted by Wang et al. [16], Yu et al. [17], and Tian et al. [18] to improve the efficiency of tool path generation methods.
The above methods optimize the tool path, aiming to improve the machining efficiency while guaranteeing the machining quality. They set the time efficiency as the optimization object but ignore the consideration of the energy consumption of the machining process. However, energy savings represent a dominant factor that could enhance the economic benefits. Thus, this paper focuses on the tool path’s energy consumption, i.e., optimizing the tool path to improve the energy efficiency during machining. It has to be noted that this paper focuses on point milling due to its high adaptivity, despite the fact that flank milling has higher machining efficiency and material removal rates. In this field, relevant research has also been carried out. Aramcharoen et al. [3] and Kara et al. [19] explored the relationships between energy consumption and different tool paths. Xu et al. [20] proposed a tool path generation method to reduce the total amount of energy consumption by constructing an algebraic model called a machine-dependent energy potential field on the surfaces of the part. Li et al. [21] presented a methodology to optimize the tool path for high efficiency, low energy consumption, and a small carbon footprint in the milling process. Hu et al. [22] introduced a model for the single objective optimization problem that minimizes the energy consumption of machining tools during the feature transitions, which include the tool path and tool change operations. However, all of these studies focused on tool path selection optimization. The tool orientation of the cutter at each cutter contact (CC) point is derived based on the largest cutting width. Hence, the tool orientation of the cutter is uniquely defined as long as the CC curve is known. However, the total energy consumption is closely related to the tool orientation of the cutter and the adjustment of the tool cutter’s posture. When the largest cutting width, there will be an admissible domain of the tool orientation of the cutter at each CC point, meaning the over-cut and global collision can be avoided in a tool path. Through more refined optimization of the tool path, for example by using the optimal tool orientation sequence, the total amount of energy consumption in a single cutting path can be further reduced. Degradation of the cutting width can be compensated for and optimized in the configuration of adjacent cutting paths. This is a feasible scheme for further reducing the total energy consumption during curved surface machining. This paper makes a contribution to this field.
The rest of this paper is organized as follows. In Section 2, a geometry analysis at the CC point is conducted and the equivalent cutting stripe width is estimated. In Section 3, the global energy consumption, including the cutting energy consumption, driving energy consumption, and idle energy consumption, is analyzed in detail. Additionally, an energy efficiency index (SEC) is defined as the optimization object. In Section 4, the discrete energy consumption path (DECP) model is constructed and the most energy-efficient tool orientation sequence is found based on this. The simulation cutting experiments are carried out in Section 5, and the energy consumption data for the proposed method and iso-scallop height method are compared, showing that the energy consumption after implementing the optimization method results in significant savings. Finally, conclusions are summarized in Section 6.

2. Geometry Analysis at the CC Point

There are actually two kinds of end cutters used in the milling process, i.e., ball-end cutters and cylinder cutters (including flat-end and toroidal cutters). A flat-end cutter is used in this paper, as shown in Figure 1. In the figure, f is the feed direction, s is the nominal surface, and s is the tolerance surface.
The local coordinate system CXCYCZC is defined at the cutter contact point (CC point) C. The cutter feeds along the direction of axis XC, ZC is the normal vector of the work piece surface at CC point, and YC is the cross-product of XC and ZC. The cutter coordinate system is defined with the center point of the cutter bottom circle O as the origin. Zt is the instantaneous cutter orientation, Xt is the tangent direction of the cutter generatrix, and Yt is the cross-product of Xt and Zt. The cutter posture can be defined with β (the angle between ZC axis and Zt axis) and σ (the angle between XC axis and Xt axis).
The following research is based on the assumption that global and local collisions are avoided. The rotary angle σ is initially assumed to be 0 to make the tilt angle the only factor that impacts the cutter orientation. The cutter’s projection along the feed direction is an ellipse. With the given cutter radius R, the equivalent radius at the CC point can be theoretically estimated as [23]:
R e = R sin β
Based on the estimation of the cutter’s equivalent cutting radius, the cutting width can be calculated according to whether the machined surface is convex or concave, as shown in Figure 1. Referring to [24], the cutting stripe width can be estimated as:
D i = 2 R sin ϕ min ( 2 R , 8 r k h )         ( C o n c a v e     c a s e ) ( C o n v e x     c a s e ) ,
where h is the constant cutting remnant depth and rk is the surface radius of curvature perpendicular to the feed direction.
To better estimate the energy efficiency of the tool path, the “specific energy consumption” index is proposed [25], which is utilized to describe the energy consumption for removing a unit volume of material. As the thickness of remnant material before the finishing machining is constant, the swept area of the tool path is needed to calculate the specific energy consumption. The swept area from CC point Pi to Pi+1, which is denoted as Areai, can be estimated as
A r e a i = 1 2 ( D i + D i + 1 ) P i + 1 P i
Additionally, the chord error e of neighboring CC points has to be within the given requirements. Correspondingly, an excessively large step size of the tool path Lf can lead to dissatisfaction of the chord error requirements. According to the geometric analysis of the machined surface, it can be realized that Lf is highly dependent on the local curvature of the machined surface. Thus, the step size can be estimated by adopting a simplistic circular arc approximation [26]:
L f = min ( 8 e 1 κ f 4 e 2 , L f f l o o r ) ,
where κf is the local curvature of the machined surface along the feed direction and Lf-floor is the set to be the floor limit if Lf is infinite. The chord error e is shown in Figure 2 below.

3. The Energy Consumption Model of a Tool Path

According to the dynamic properties of the machine tool, the energy consumption can be divided into two parts, i.e., the cutting energy consumption and the driving energy consumption. The cutting energy consumption is used to strip extraneous materials away from the workpiece and the driving energy consumption is used to adjust the position and posture of the cutter. In this section, Xu’s work [20] and our previous work [27,28,29] are extended to establish an energy consumption model of a tool path.

3.1. The Cutting Energy Consumption Analysis of the Tool Path

The profile of the flat-end mill cutter at the high rotation speeds is actually a cylinder. Theoretically, the cutting edge on a tool cutter can be modeled as:
p = ( X t , Y t , Z t ) = ( R sin φ , R cos φ , z ) ,
where φ is the rotating angle of a specific reference point on the cutting edge. As the helix angle i0 of the cutter is constant, the rotating angle φ can be estimated as:
φ = φ 0 z R tan i 0
By analyzing a single cutting edge, the rest of the cutting edges can be modeled by rotating the cutter with an angular phase [30]. Thus, the cutting force on the cutting edge can be resolved into three components, i.e., the tangent, radial, and sub-tangent differential cutting forces, as shown in Equation (7):
d F t = K t c T ( φ ) d z + K t e d z d F r = K r c T ( φ ) d z + K r e d z d F a = K a c T ( φ ) d z + K a e d z ,
where T ( φ ) = max ( 0 , f t p c ) is the chip thickness; Ktc, Krc, and Kac represent the shear force factors; and Kte, Kre, and Kae are the edge factors along the tangent, radial, and sub-tangent directions, respectively.
To calculate the cutting force on the cutting edges by integrating the differential cutting force elements, three components are expressed in the tool coordinate system as:
d F t x d F t y d F t z = sin φ cos φ 0 cos φ sin φ 0 0 0 1 d F r d F t d F a
After the transformation, the instantaneous cutting force can be obtained by integrating the differential terms, and the cutting force during one period of rotation can be estimated by integrating the differential term of time. It should be noted that when the point on the cutting edge is no longer merged into the remnant material, the differential term [dFtx,dFtx,dFtx]T should be set to be [0,0,0]T because of the disappearance of the cutting force. In this case, the tangential force Ft is the only component of the cutting force that contributes to the cutting torque on the flat-end cutter. Thus, the energy consumption is derived as [31]:
E T = P T p i + 1 p i f = 0 Δ t 0 l 0 F t ( h , l ) R S d l d t ,
where Ft (h, t) is the tangential force intensity at a certain height l during the cutting time t, R is the tool radius, and S is the spindle speed.

3.2. The Driving Energy Consumption Analysis of the Tool Path

Regarding the driving energy consumption ED, the energy consumption of each axis is included and should be investigated in the machine coordinate system (MCS), where the driving energy is needed to overcome internal frictional force torques. The frictional force of each axis changes with its velocity and acceleration. With the given chord error e, three adjacent CC points along a tool path can be easily calculated, as shown in Figure 2. A cutter posture can be defined as:
p i = ( x i , y i , z i , a i , b i , c i ) ,
where (ai, bi, ci) is the given tool orientation vector. Then inverse kinematics transformation is applied to the tool posture, and the tool cutter’s posture in the machine coordinate system can be obtained:
A x i s = ( a x i , 1 , a x i , 2 , a x i , 3 , a x i , 4 , a x i , 5 )                   = I K T ( x i , y i , z i , a i , b i , c i )                    
The detailed transformation process can be found in [32]. Thus, the velocity and acceleration of each axis can be estimated as:
v k = ( a x i + 1 , k a x i , k ) / Δ t a k = ( a x i + 1 , k 2 a x i , k + a x i 1 , k ) / Δ t 2         k = 1 , 2 , 3 , 4 , 5
According to the aforementioned analysis, the driving energy consumption of the machining process can be estimated as:
E D = k = 1 5 ( μ k J k Δ t + F k v k 2 Δ t + ζ C R ) ζ C R = 1 2 J k ( v k + a k Δ t ) 2 v k 2 , i f         a k > 0                                                           0 ,                           i f         a k 0 ,
where μ k J k Δ t + F k v k 2 Δ t represents the energy consumption when moving a driving axis at a constant speed and ζ C R represents the energy consumption needed to change the rotation speed of a driving axis.
It is worth noting that the last term in Equation (13) represents the energy consumption caused by the kinetic properties of each axis, which decreases to 0 when the corresponding axis is decelerating, i.e., ak ≤ 0.

3.3. Global Energy Consumption of a Tool Path

With ED and ET defined, the global energy consumption requirement of a machining process can be estimated as:
E = P i d l e t + E T + E D ,
where Pidle is the idle running power of the machine tool and t is the time requirement of the machining process. In general optimization methods, the width of the cutting stripe is set as the only optimization target. To obtain the widest cutting stripe, the cutter’s postures at each CC point are adjusted to make the curvatures of the tool cutter and machined surface match. According to the cutting energy consumption analysis, it can be realized that the cutting energy consumption rises with the cutting stripe width. The driving energy consumption of each axis meanwhile rises because of the adjustment of the tool cutter’s posture. However, since the cutting stripes are enlarged, the machining time can be shortened and the variation in total energy can be determined by the aforementioned analysis. This means that it is useless to obtain the largest energy-efficient cutting stripe by adjusting the tool cutter’s posture. Thus, with the pre-determined machining parameters in Equations (1)–(4), the “specific energy consumption (SEC)” index should be set as the criterion to evaluate the energy efficiency of a machining process, as shown in Equation (15).
S E C = E / A r e a = P i d l e t + E T + E D 1 2 ( D i + D i + 1 ) P i + 1 P i

4. The Energy Efficiency Optimization of the Tool Path Using the Discrete Energy Consumption Path Model

To find the optimal solution to the global energy consumption model proposed in Section 3, the following three steps are carried out: (1) discretizing the domain of admissible orientations at a CC point to transform an NP hard problem into an exhaustion problem; (2) constructing the discrete energy consumption path model to describe the exhaustion problem; (3) finding the shortest path in the discrete energy consumption path model to obtain the energy-efficient orientation sequence of a single tool path.

4.1. Discretization of the Domain of Admissible Orientations at a CC Point

As the tilt angle is the only variable expressing the tool orientation, the admissible domain of a CC point can be visualized as an arc that lies between the corresponding minimal and maximal admissible tilt angles. A minimal inclination angle is the inclination angle limit to avoid over-cutting. A maximal inclination angle is the inclination angle limit to avoid a global collision.
According to the inverse kinematic transformation analysis of the AC-type machine tool, a specific tool orientation is determined by the angle position of two rotary axes in AC-type machine tools. The mapping relation of rotary axes’ angle position and the tool orientation vector can be expressed by:
sin α sin γ = x sin α cos γ = y cos α = z ,
where (α, γ) defines the rotary axes’ angle position and (x, y, z) defines the tool orientation vector in WCS.
Inspired by Plakhotnik [33], we realized that the discretization of admissible domains of tool orientations is efficient to find an optimal solution for a tool path. The aim is minimizing the SEC of a tool path, and the tool orientation can be expressed by its admissible domain [τmin,τmax], such that the energy-efficient solution can be obtained by finding a sequence of τs at each CC point.
Each discrete domain of admissible orientations is composed of a group of nodes that represent a specific posture of the tool cutter. In order to describe the tool orientation sequence of the tool path, one node from each volume is selected and connected as shown in Figure 3. It has to be noted that the dots in Figure 3 are physically virtual, and each of them represents a potential choice of tool orientation. Regarding the edges, they represent the rotary movements transforming from one posture to another. The optimal tool orientation sequence can be solved by finding the shortest path from the start point to the end point.

4.2. Construction of Discrete Energy Consumption Path Model

Based on the discretization of admissible domains of tool orientations, the DECP model can be constructed in a similar way to that shown in Figure 4. As the cutting energy consumption ET is only relevant to the posture of the tool cutter, the cutting energy of each posture at a CC point is a form of static data that can be stored in the volumes of nodes.
Regarding the driving energy consumption ED, this is determined by how the tool orientation changes, meaning it is a form of dynamic data that can be stored in edges connecting two neighboring volumes of nodes. The swept area per feed can be stored in the same way as the driving energy consumption ED.
Based on the aforementioned profile of the DECP model, it can be detailed as follows:
Definition of nodes: According to Equations (7)–(9), the cutting energy consumption is only relevant to the posture of the tool cutter, meaning the jth node at the ith CC point corresponds to a cutting energy consumption value ET(i,j). The cutting energy consumption sequence of a tool path can be expressed as:
S T = { E T ( 1 , j 1 ) , E T ( 2 , j 2 ) , , E T ( n , j m ) } , m = 1 , 2 , , n       ,
where n is the CC point number of a tool path. Additionally, the positions of two rotary axes in the machine coordinate system are stored in volumes of nodes as:
P i j = A ( i , j ) , C ( i , j )
Thus, a specific node that contains information regarding the cutting energy and positions on rotary axes can be expressed as:
N d i j = { P i j , E T ( i , j ) }
Definition of edges: According to Equation (13), the driving energy consumption is relevant to the angle changes of two rotary axes, meaning the movements from one posture to the next should be defined first as:
Δ A x i s i = A x i s ( i + 1 , j ) A x i s ( i , j ) ,
where it should be noted that ∆Axis is a vector that can indicate whether the rotary axes are in acceleration or deceleration.
The velocity of two rotary axes can be estimated according to Equation (21):
v A i = ( A ( i + 1 , j ) A ( i , j ) ) / Δ t i v C i = ( C ( i + 1 j ) C ( i , j ) ) / Δ t i
The corresponding acceleration can be estimated similarly:
a A i = A ( i + 2 , j ) A ( i + 1 , j ) Δ t i + 1 A ( i + 1 , j ) A ( i , j ) Δ t i 1 2 ( Δ t i + Δ t i + 1 ) a C i = C ( i + 2 , j ) C ( i + 1 , j ) Δ t i + 1 C ( i + 1 , j ) C ( i , j ) Δ t i 1 2 ( Δ t i + Δ t i + 1 )
To avoid confusion regarding the variable naming in Equations (20)–(22), there should be a special illustration; if the jth node in the ith column is selected, when it comes to the (i + 1)th column, the number of selected nodes will be j’. In the same way, the naming will be conducted in the (i + 1)th column, i.e., the number of the selected node will be j’’. Thus, the driving energy consumed between two adjacent postures can be estimated as:
E D ( i , j ) ( i , j ) = B o t h       R o t a r y     A x e s ( μ J Δ t + F v i 2 Δ t + Δ E       ) Δ E   = 1 2 J ( v i + a i Δ t ) 2 v i 2 , i f         a i > 0                                                           0                       ,   i f         a i 0 ,
where μ, F, and J are respectively the friction coefficient, velocity friction coefficient, and inertia of each rotary axis. Thus, a specific edge that connects two neighboring nodes can be expressed as:
E d ( i , j ) ( i , j ) = { E D ( i , j ) ( i , j ) , A r e a i }
In order to optimize the energy efficiency of a tool path, the specific energy consumption is set as the cost function of the optimization process as:
min T h r o u g h T h e P a t h ( E T + E D ) T h r o u g h T h e P a t h A r e a

4.3. Finding the Most Energy-Efficient Path in the DECP Model

To find the most energy-efficient path in the DECP model, an energy efficiency solution graph (EESG) is first constructed. Columns of nodes and edges connecting adjacent columns are contained in the EESG. The frame of the EESG inherits that of the DECP, i.e., each column of nodes represents the specific orientations at each CC point, as shown in Figure 5. As the energy consumption Eij-i’j’ and swept area Areaij-i’j’ from one specific tool orientation to another can be calculated, the edges connecting two arbitrary nodes of neighboring columns record the energy consumption and the swept area from one tool posture to another, i.e.,
E d g e ( i j i j ) = { E i j i j , A r e a i j i j }
The value of each node is assigned as the optimal solution from the start point. The assignment can be concluded as:
N o d e i j = min ( E i 1 , 1 + E ( i 1 , 1 ) ( i , j ) A r e a i 1 , 1 + A r e a ( i 1 , 1 ) ( i , j ) ,   E i 1 , 2 + E ( i 1 , 2 ) ( i , j ) A r e a i 1 , 2 + A r e a ( i 1 , 2 ) ( i , j ) ,   , E i 1 , N + E ( i 1 , N ) ( i , j ) A r e a i 1 , N + A r e a ( i 1 , N ) ( i , j ) ) ,
where N is the total number of nodes in the (i − 1)th column.
The nodes are assigned column by column and the node with the lowest assignment represents the selected node of the optimal orientation sequence. Then, the most energy-efficient path can be obtained by following the selected edges from the end node to the start node. The validity of the assignment can be illustrated as follows:
(1)
Illustration of the initial assignments.
When I is initially set as 1, the graph shown in Figure 6 represents the first two columns of nodes. The node with the lowest assignment in (i + 1)th (that is, the 2nd) column indicates the selected node that is involved in the optimal orientation sequence.
(2)
Illustration of the transitivity.
From the aforementioned illustration, it can be concluded that each node of the (i − 1)th column represents the most energy-efficient path, meaning the nodes of ith column can inherit the superiority from the (i − 1)th column, and the minimization process of Equation (27) can select the most energy-efficient path. AS the value of i increases, each column of nodes is assigned with the lowest SEC.

5. Experiment and Discussion

Firstly, a simulation experiment is carried out. The whole optimization process is implemented in MATLAB R2017b and is run on a computer (Intel/Core i5 CPU, 3.10 GHz, and 6 GB RAM with a Windows 7 operating system). A five-axis machine tool with a double pendulum working table is chosen to experiment on. Without loss of generality, the configuration of the machine tool is chosen as the orthogonal AC type, i.e., the motion axes are composed of three orthogonal translational axes and two orthogonal rotary axes.
The machining parameters used in the simulation, including configuration parameters and graphical parameters, are shown in Table 1. Moreover, the accuracy of the optimized orientation sequence will be higher if the tool orientations are sampled at a higher density. However, the marginal unity when increasing the sampling efficiency is degressive. Therefore, the sampling density is set to 0.1° based on our previous work [29,34] in this paper. In addition, it should be pointed out that the determination of other related authentic intrinsic parameters of a particular five-axis machine, e.g., coefficients involved in the calculation of the cutting force, are in another important paper. The values of these coefficients are determined by referring to the relevant references [20,21,29] in this paper.
A free-form surface is chosen to prove the validity of the proposed method. The offset distance is given as 3 mm, i.e., the thickness of the remnant material on the surface is 3 mm. To obtain the largest cutting stripe width, the initial tool paths are generated by the iso-scallop height algorithm on the free-form surface and a single path is chosen to verify the optimization method.
The chosen path is firstly discretized into a series of points on the path according to Equation (4) to satisfy the requirement of the chord error. Thus, 94 points are extracted from the chosen path to constitute the optimization domain of the proposed energy-saving optimization method. With the length of the tool path L = 203.4 mm and the feed rate FR = 600 mm/min substituted, the idle energy consumption during the machining process along the chosen path can be easily estimated: Eidle = Pidle·L/FR = 4068 J.
Next, for comparative analysis, the cutting energy and driving energy respectively consumed by the iso-scallop height algorithm and the method proposed in this paper are calculated. For the iso-scallop height algorithm, the cutting energy and driving energy consumed are ET = 10916.34 J and ED = 2015.76 J, respectively, and the accumulated machining area is Areainit = 175.48 mm2; thus, the specific energy of the tool path is SECinit = 96.88 J/mm2 considering the idle energy consumption Eidle = 4068 J. Correspondingly, for the method proposed in this paper, the machining energy and driving energy consumed are respectively ET-Opti = 8718.37 J and ED-Opti = 2045.23 J, and the accumulated machining area is Areainit = 170.11 mm2; thus, the specific energy of the tool path is SECOpti = 87.19 J/mm2 considering the idle energy consumption Eidle = 4068 J.
As shown in Table 2, it can be found that the machining energy decreases with the deterioration of the cutting stripe width, while the driving energy stays relevantly stable. However, the specific energy consumption of the tool path is optimized, i.e., the specific energy consumption decreases from 96.88 J/mm2 to 87.19 J/mm2, which implies a 10.01% energy saving when removing a specific volume of material. Among the three energy consumption components, the idle running energy is unchanged because only a single tool path is studied, the driving energy stays relevantly stable, and the change of cutting energy is relatively large. This is due to the shearing energy no longer being the same under different tool orientation sequences, even with the same volume of material removed, because the cutter–workpiece engagement during the cutting process has changed. Based on the further comparative analysis of some segments of the tool path, it can be found that larger cutting stripe widths are obtained by changing the tool orientation acutely, although the cutting energy increases. Sometimes it is not economically feasible to drive rotary axes of the machine tool to obtain the largest cutting stripe because a nonlinear relationship exists with the driving energy, meaning the marginal benefits will be reduced in certain segments of the tool path. Therefore, it is not smart to acquire a large MMR blindly. The proposed optimization method makes the rotary motions of the machine tool smarter, meaning the machine tool can adjust its rotary axes to obtain the largest marginal energy benefit. Then, the SEC value when cutting a single tool path is reduced and also lays the foundation for further reductions in global energy consumption during the machining of complete curved surfaces.

6. Conclusions

This paper proposed an optimization method that can generate an energy-efficient tool orientation sequence for a tool path. This paper is innovative in the following aspects:
(1)
Based on the geometry analysis at the CC point, the global energy consumption model of a tool path is constructed. The global energy consumption of a tool path is divided into 3 constituents, namely the cutting energy, driving energy, and idle energy, and each aspect is evaluated;
(2)
The computation cost from searching for an optimal path in a series of feasible arcs is unacceptably large. The DECP model is constructed to discretize the tool orientation domains at each CC point within a specific accuracy requirement. Thus, the searching problem is converted to a shortest-path searching problem in a directed graph;
(3)
The most energy-efficient tool orientation sequence of the tool path can be found by finding the shortest path in the EESG. It can be found that the cutting energy decreases with the deterioration of the cutting stripe width; the driving energy stays relatively stable. According to the data on the energy consumption and machined area, the specific energy consumption of the tool path is optimized to prove the method’s validity.
Even though the energy saving of 10% is not remarkable, the quantity of the energy cost savings could be significantly large. This could be an effective approach for lean machining that could make the production process cleaner. As the proposed method can only optimize the energy efficiency of a single tool path, only simulation experiments based on a five-axis machine tool with a double pendulum working table and flat-end cutter were carried out in this paper. In future work, the overlap between adjacent tool paths will be considered, an integrated optimization method that could improve the energy efficiency of tool paths on free-form surfaces will be paid more attention, and empirical research will be carried out.

Author Contributions

Conceptualization, Y.G. and S.M.; methodology, S.M.; validation, S.M. and H.Z.; formal analysis, H.Z.; investigation, Q.W.; resources, Z.W.; data curation, S.M.; writing—original draft preparation, S.M. and Q.W.; writing—review and editing, Y.G. and H.Z.; visualization, Q.W.; supervision, Z.W.; project administration, Y.G.; funding acquisition, Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of Zhejiang Province (No. LR22E050006) and the National Natural Science Foundation of China (No. 51975386).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Edem, I.F.; Balogun, V.A.; Mativenga, P. An investigation on the impact of toolpath strategies and machine tool axes configurations on electrical energy demand in mechanical machining. Int. J. Adv. Manuf. Technol. 2017, 92, 2503–2509. [Google Scholar] [CrossRef]
  2. Zhou, L.; Li, J.; Li, F.; Meng, Q.; Li, J.; Xu, X. Energy consumption model and energy efficiency of machine tools: A comprehensive literature review. J. Clean. Prod. 2016, 112, 3721–3734. [Google Scholar] [CrossRef]
  3. Aramcharoen, A.; Mativenga, P. Critical factors in energy demand modelling for CNC milling and impact of toolpath strategy. J. Clean. Prod. 2014, 78, 63–74. [Google Scholar] [CrossRef]
  4. Xiao, Q.; Li, C.; Chen, X.; Tang, Y. An investigation into the dependence of energy efficiency on CNC process parameters with a sustainable consideration of electricity and materials. In Proceedings of the 2017 IEEE International Conference on Systems, Man, and Cybernetics (SMC), Banff, AB, Canada, 5–8 October 2017. [Google Scholar] [CrossRef]
  5. Arzhanov, K.V.; Shurygin, Y.A. Energy saving control algorithms of stepper motor drives for metal-working machines equipped with CNC unit. In Proceedings of the 2015 International Conference on Mechanical Engineering, Automation and Control Systems (MEACS), Tomsk, Russia, 1–4 December 2015. [Google Scholar] [CrossRef]
  6. Sato, R.; Shirase, K.; Hayashi, A. Energy Consumption of Feed Drive Systems Based on Workpiece Setting Position in Five-Axis Machining Center. J. Manuf. Sci. Eng. 2018, 140, 021008. [Google Scholar] [CrossRef]
  7. Kant, G.; Sangwan, K.S. Predictive Modelling for Energy Consumption in Machining Using Artificial Neural Network. Procedia CIRP 2015, 37, 205–210. [Google Scholar] [CrossRef]
  8. Mouzon, G.C.; Yildirim, M.B.; Twomey, J.M. Operational methods for minimization of energy consumption of manufacturing equipment. Int. J. Prod. Res. 2007, 45, 4247–4271. [Google Scholar] [CrossRef] [Green Version]
  9. Tian, G.; Zhou, M.; Li, P. Disassembly Sequence Planning Considering Fuzzy Component Quality and Varying Operational Cost. IEEE Trans. Autom. Sci. Eng. 2017, 15, 748–760. [Google Scholar] [CrossRef]
  10. Gao, Y.; Feng, Y.; Wang, Q.; Zheng, H.; Tan, J. A multi-objective decision making approach for dealing with uncertainty in EOL product recovery. J. Clean. Prod. 2018, 204, 712–725. [Google Scholar] [CrossRef]
  11. Gao, Y.; Feng, Y.; Zhang, Z.; Tan, J. An optimal dynamic interval preventive maintenance scheduling for series systems. Reliab. Eng. Syst. Saf. 2015, 142, 19–30. [Google Scholar] [CrossRef]
  12. Marciniak, K. Influence of surface shape on admissible tool positions in 5-axis face milling. Comput. Des. 1987, 19, 233–236. [Google Scholar] [CrossRef]
  13. Mi, Z.; Yuan, C.-M.; Ma, X.; Shen, L.-Y. Tool orientation optimization for 5-axis machining with C-space method. Int. J. Adv. Manuf. Technol. 2017, 88, 1243–1255. [Google Scholar] [CrossRef]
  14. Kim, T.; Sarma, S.E. Toolpath generation along directions of maximum kinematic performance; a first cut at machine-optimal paths. Comput. Des. 2002, 34, 453–468. [Google Scholar] [CrossRef]
  15. Cong, G.; Li, M.; Chunyu, W. A tool path control algorithm based on tool orientation vector interpolation in five-axis CNC machining. In Proceedings of the 2016 Sixth International Conference on Instrumentation & Measurement, Computer, Communication and Control (IMCCC), Harbin, China, 21–23 July 2016. [Google Scholar] [CrossRef]
  16. Wang, J.; Liu, Q.; Pi, S.; Liu, Q.; Li, Y. Research on algorithm of tool-path planning for 5-axis NC machining based on double quaternion. In Proceedings of the 2017 IEEE International Conference on Mechatronics and Automation (ICMA), Takamatsu, Japan, 6–9 August 2017. [Google Scholar] [CrossRef]
  17. Yu, W.; Bi, Y.; Li, Z. Research on Tool Path Planning Method of NURBS Surface Based on CPU-GPU Parallel Computing. In Proceedings of the 2017 International Conference on Computer Network, Electronic and Automation (ICCNEA), Xi’an, China, 23–25 September 2017. [Google Scholar] [CrossRef]
  18. Tian, G.; Ren, Y.; Feng, Y.; Zhou, M.; Zhang, H.; Tan, J. Modeling and Planning for Dual-Objective Selective Disassembly Using and/or Graph and Discrete Artificial Bee Colony. IEEE Trans. Ind. Inform. 2018, 15, 2456–2468. [Google Scholar] [CrossRef]
  19. Kara, S.; Li, W. Unit process energy consumption models for material removal processes. CIRP Ann. 2011, 60, 37–40. [Google Scholar] [CrossRef]
  20. Xu, K.; Luo, M.; Tang, K. Machine based energy-saving tool path generation for five-axis end milling of freeform surfaces. J. Clean. Prod. 2016, 139, 1207–1223. [Google Scholar] [CrossRef]
  21. Li, L.; Deng, X.; Zhao, J.; Zhao, F.; Sutherland, J.W. Multi-objective optimization of tool path considering efficiency, energy-saving and carbon-emission for free-form surface milling. J. Clean. Prod. 2018, 172, 3311–3322. [Google Scholar] [CrossRef]
  22. Hu, L.; Liu, Y.; Peng, C.; Tang, W.; Tang, R.; Tiwari, A. Minimising the energy consumption of tool change and tool path of machining by sequencing the features. Energy 2018, 147, 390–402. [Google Scholar] [CrossRef] [Green Version]
  23. Rao, N.; Ismail, F.; Bedi, S. Tool path planning for five-axis machining using the principal axis method. Int. J. Mach. Tools Manuf. 1997, 37, 1025–1040. [Google Scholar] [CrossRef]
  24. Lo, C.-C. Efficient cutter-path planning for five-axis surface machining with a flat-end cutter. Comput. Des. 1999, 31, 557–566. [Google Scholar] [CrossRef]
  25. Velchev, S.; Kolev, I.; Ivanov, K.; Gechevski, S. Empirical models for specific energy consumption and optimization of cutting parameters for minimizing energy consumption during turning. J. Clean. Prod. 2014, 80, 139–149. [Google Scholar] [CrossRef]
  26. Li, H.; Feng, H.-Y. Efficient five-axis machining of free-form surfaces with constant scallop height tool paths. Int. J. Prod. Res. 2004, 42, 2403–2417. [Google Scholar] [CrossRef]
  27. Gao, Y.; Wang, Q.; Feng, Y.; Zheng, H.; Zheng, B.; Tan, J. An Energy-Saving Optimization Method of Dynamic Scheduling for Disassembly Line. Energies 2018, 11, 1261. [Google Scholar] [CrossRef] [Green Version]
  28. Feng, Y.; Wang, Q.; Gao, Y.; Cheng, J.; Tan, J. Energy-Efficient Job-Shop Dynamic Scheduling System Based on the Cyber-Physical Energy-Monitoring System. IEEE Access 2018, 6, 52238–52247. [Google Scholar] [CrossRef]
  29. Wang, Q.; Feng, Y.; Gao, Y.; Li, Z.; Tan, J. Smooth fillet-end cutter tool path generation method on triangular-mesh surface based on Modified Butterfly subdivision. Int. J. Adv. Manuf. Technol. 2018, 98, 2831–2847. [Google Scholar] [CrossRef]
  30. Lin, B.; Wang, L.; Guo, Y.; Yao, J. Modeling of cutting forces in end milling based on oblique cutting analysis. Int. J. Adv. Manuf. Technol. 2016, 84, 727–736. [Google Scholar] [CrossRef]
  31. Xu, K.; Tang, K. An Energy Saving Approach for Rough Milling Tool Path Planning. Comput.-Aided Des. Appl. 2016, 13, 253–264. [Google Scholar] [CrossRef] [Green Version]
  32. Lin, Z.; Fu, J.; Shen, H.; Gan, W. Non-singular tool path planning by translating tool orientations in C-space. Int. J. Adv. Manuf. Technol. 2014, 71, 1835–1848. [Google Scholar] [CrossRef]
  33. Plakhotnik, D.; Lauwers, B. Graph-based optimization of five-axis machine tool movements by varying tool orientation. Int. J. Adv. Manuf. Technol. 2014, 74, 307–318. [Google Scholar] [CrossRef]
  34. Wang, Q.; Feng, Y.; Zhang, Z.; Tan, J. Tool orientation sequence smoothing method based on the discrete domain of feasible orientations. Int. J. Adv. Manuf. Technol. 2017, 92, 4501–4510. [Google Scholar] [CrossRef]
Figure 1. Cutting stripe width of flat-end milling for concave and convex cases. (a) Cutting stripe width diagram for concave case, (b) cutting stripe width diagram for concave case, and (c) five axis cutting diagram of the flat-end cutter.
Figure 1. Cutting stripe width of flat-end milling for concave and convex cases. (a) Cutting stripe width diagram for concave case, (b) cutting stripe width diagram for concave case, and (c) five axis cutting diagram of the flat-end cutter.
Machines 10 00348 g001
Figure 2. The tool cutter postures of neighboring CC points along a tool path.
Figure 2. The tool cutter postures of neighboring CC points along a tool path.
Machines 10 00348 g002
Figure 3. Discretization of admissible domains of CC points.
Figure 3. Discretization of admissible domains of CC points.
Machines 10 00348 g003
Figure 4. Discrete energy consumption path.
Figure 4. Discrete energy consumption path.
Machines 10 00348 g004
Figure 5. Energy efficiency solution graph.
Figure 5. Energy efficiency solution graph.
Machines 10 00348 g005
Figure 6. Illustration of an EESG.
Figure 6. Illustration of an EESG.
Machines 10 00348 g006
Table 1. Configuration and machining parameters in the simulation.
Table 1. Configuration and machining parameters in the simulation.
Parameter NamesParameter Values
Radius of the flat-end tool cutterR = 6 mm
Number of cutting edges on the cuttersNc = 2
Helix angle of the cutteri0 = 30°
Spindle speedS = 3000 r/min
Idle power of the machine centerPidle = 0.25 kW
Mass and moments of inertia of the five axesMX = 800 Kg, MY = 800 Kg, MZ = 100 Kg, JA = 15 Kg·m2, JC = 10 Kg·m2
Viscosity friction coefficients of five axesFX = 1000 N·m, FY = 1000 N·s/m, FZ = 1000 N·s/m, FA = 100 N·s/m, FC = 100 N·s/m
Friction coefficients of five axesμ = 1
Coefficients in calculation of cutting forceKte = 25 N/mm, Kre = 42 N/mm, Kae = 6 N/mm, Ktc = 1600 N/mm2, Krc = 440 N/mm2, Kac = 591 N/mm2
Feed rateFR = 600 mm/min
Scallop height toleranceh = 0.5 mm
Chord errore = 0.3 mm
Sampling density on admissible arc0.1°
Table 2. Comparison of the specific energy values of the tool paths.
Table 2. Comparison of the specific energy values of the tool paths.
The Iso-Scallop Height AlgorithmThe Method Proposed in This Paper
L = 203.4 mmL = 203.4 mm
FR = 600 mm/minFR = 600 mm/min
Eidle = Pidle·L/FR = 4068 JEidle = Pidle·L/FR = 4068 J
ET = 10916.34 JET-Opti = 8718.37 J
ED = 2015.76 JED-Opti = 2045.23 J
Areainit = 175.48 mm2Areainit = 170.11 mm2
SECinit = 96.88 J/mm2SECOpti = 87.19 J/mm2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gao, Y.; Mi, S.; Zheng, H.; Wang, Q.; Wei, Z. An Energy Efficiency Tool Path Optimization Method Using a Discrete Energy Consumption Path Model. Machines 2022, 10, 348. https://doi.org/10.3390/machines10050348

AMA Style

Gao Y, Mi S, Zheng H, Wang Q, Wei Z. An Energy Efficiency Tool Path Optimization Method Using a Discrete Energy Consumption Path Model. Machines. 2022; 10(5):348. https://doi.org/10.3390/machines10050348

Chicago/Turabian Style

Gao, Yicong, Shanghua Mi, Hao Zheng, Qirui Wang, and Zhe Wei. 2022. "An Energy Efficiency Tool Path Optimization Method Using a Discrete Energy Consumption Path Model" Machines 10, no. 5: 348. https://doi.org/10.3390/machines10050348

APA Style

Gao, Y., Mi, S., Zheng, H., Wang, Q., & Wei, Z. (2022). An Energy Efficiency Tool Path Optimization Method Using a Discrete Energy Consumption Path Model. Machines, 10(5), 348. https://doi.org/10.3390/machines10050348

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop