Next Article in Journal
Markov Moment Problem and Sandwich Conditions on Bounded Linear Operators in Terms of Quadratic Forms
Next Article in Special Issue
Relaxed Variable Metric Primal-Dual Fixed-Point Algorithm with Applications
Previous Article in Journal
Some New Mathematical Integral Inequalities Pertaining to Generalized Harmonic Convexity with Applications
Previous Article in Special Issue
A High-Precision Surrogate Modeling Method Based on Parallel Multipoint Expected Improvement Point Infill Criteria for Complex Simulation Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effective Optimization Based on Equilibrium Optimizer for Dynamic Cutting Force Coefficients of the End-Milling Process

1
Department of Mechanical Engineering, Thai Nguyen University of Technology, 3/2 Street, Tich Luong Ward, Thai Nguyen 250000, Vietnam
2
Department of Electrical Engineering, National Kaohsiung University of Science and Technology, Kaohsiung 807618, Taiwan
3
Department of Electrical Engineering, Faculty of Engineering at Shoubra, Benha University, Cairo 11629, Egypt
4
Faculty of International Training, Thai Nguyen University of Technology, 3/2 Street, Tich Luong Ward, Thai Nguyen 250000, Vietnam
5
Department of Electrical Engineering, College of Engineering, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Mathematics 2022, 10(18), 3287; https://doi.org/10.3390/math10183287
Submission received: 10 August 2022 / Revised: 6 September 2022 / Accepted: 7 September 2022 / Published: 10 September 2022
(This article belongs to the Special Issue Numerical Analysis and Optimization: Methods and Applications)

Abstract

:
This study aims to develop an accurate dynamic cutting force model in the milling process. In the proposed model, the estimated cutting force tackles the effect of the self-excited vibration that causes machining instability during the cutting process. In particular, the square root of the residual cutting force between the prediction and the actual cutting force is considered as an objective function for optimizing the cutting force coefficients using the equilibrium optimizer (EO) approach instead of the trial-and-error approach. The results confirm that the proposed model can provide higher prediction accuracy when the EO is applied. In addition, the proposed EO has a minimum integral square error (ISE) of around 1.12, while the genetic algorithm (GA) has an ISE of around 1.14 and the trial-and-error method has an ISE of around 2.4. Moreover, the proposed method can help to investigate the cutting stability and to suspend the chatter phenomenon by selecting an optimal set of cutting parameters.

1. Introduction

Milling technology plays an important role in the manufacturing sector, with many applications in the automotive, aerospace, and mobility industries, among others [1]. Milling is the process of removing material with a rotary cutter by advancing a cutter into a workpiece. Thus, improvements in the product quality and productivity of the milling process are key contributors to enhancing efficiency and decreasing the cost when making parts and components [2]. Generally, modeling of dynamic cutting forces during the cutting is required to optimize the performance of the machining process [3]. Moreover, an accurate prediction of cutting force is very useful for tool condition monitoring and investigating cutting stability [4]. However, one of the major issues during the cutting process is self-excited vibration under aggressive cutting conditions, generating a phenomenon called chatter [5,6,7]. These severe vibrations not only cause reduced quality of the surface finish, poor dimensional accuracy, and noise, but also damage the machine tool itself. Undoubtedly, chatter investigation—specifically, detecting the early state of chatter—is important to suppress chatter and maintain stable cutting conditions [8,9,10,11]. Many researchers have utilized the cutting signal, vibration, and sound signals, as well as stability lobe diagrams, for chatter detection and the prediction of milling stability [12,13,14,15,16]. The sound signal can become unreliable, based on some previous work where the vibration signal only focused on spontaneous high frequency during the cutting process. In this context, the cutting force is utilized as a physical quantity to demonstrate the machining process performance and quality. It is critical to determine the precise cutting force in order to investigate chatter characteristics, overload, and tool wear issues [17]. Taner et al. [3] introduced a common cutting force approach that predicts cutting force by detecting cutting force factors. Sun et al. [18] presented a dynamic milling parameter prediction model that takes into account the deformation problem. Zhang et al. [19] established an effective cutting force prediction approach for a five-axis milling machine of a sculptured surface using an adequate calibration strategy for cutter runout factors and particular cutting force characteristics. However, this cutting force model has disadvantages in that it is designed for a specific situation and the model coefficient must be calibrated when the cutting conditions change. In order to be more convenient for the analysis of machine spindles, Ji et al. [20] suggested a technique for predicting the dynamic behaviors of an advanced tool tip by taking into account the contact dynamic features within the tool and its holder. This modeling utilizes the current tool hammer experiment and finite element approach. Zhang et al. [21] introduced a substructure response analysis approach that utilizes the machine tool as a coupling of discrete pieces, allowing the dynamic characteristics of the tool tip to be discovered to investigate chatter throughout the cutting process. In milling processes, the size of the cutting force is related to the previous cutting behavior. By extracting the cutting force, Liu et al. [22] exploited a method of a dynamic cutting force model, in which the cutting force is only generated when there is contact between the tool and the workpiece. It forms an intermittent input with a time delay term in the cutting dynamics equation, which is a nonlinear system with progressive chaos characteristics (route-to-chaos). In this case, the Fourier spectrum cannot be used for general analysis of the chatter due to the nonlinear cutting behavior. Therefore, a method based on Tony L. Schmitz’s model [23] to discuss the cause and frequency range of chatter was proposed, and it was compared with actual milling experiments. On the other hand, the Lyapunov exponent index used in chaotic dynamics was used to analyze the time-domain changes in cutting force and establish an objective basis for judging chatter phenomena [24]. Hilbert–Huang Transform (HHT) with respect to chatter frequency and energy through the decomposing signal to empirical modes can be used to deal with nonlinear cutting behavior [25,26].
The accuracy of the cutting force model in cutting processes depends largely on the modeling and identification of cutting force coefficients [27,28]. There are two main methods to identify cutting force coefficients: (1) using cutting mechanics and tool geometry, and (2) using specific coefficients derived from actual experimental results. For the cutting mechanics and tool geometry, one of the most common approaches is called “orthogonal to oblique transformation”, proposed by Budak et al. [29]. This is a generic method for identifying cutting force coefficients for various cutting tools and procedures using data derived from orthogonal cutting tests. Cutting mechanics coefficients are more adaptable since they may be applied to any tool geometry due to the orthogonal to oblique transformation. On the other hand, under the conditions of the same tool–material combination used in an experimental test, the calculation of specific coefficients from actual experiments normally gives a higher-accuracy result. The estimation of the cutting force during the process completely relies on the precision of the empirical cutting force factors. In order to adjust the factors, extensive experiments have been performed in various cutting conditions [30]. The experiment is usually conducted in low-speed experiments in order to control the dynamic problem generated by cutting force measurement devices. However, the primary disadvantage of this technique is that the determined coefficients are used in the simulation of a generic machining process at various spindle speeds. This might be a problem because the cutting process and chip formation mechanics change with cutting speed, implying a change in coefficients as well. Moreover, the shear angle oscillations generated by vibration, cutting speed, and tool flank–wavy surface contact mechanism all influence dynamic cutting forces [31].
Recently, enormous optimization strategies have been developed for the purpose of tuning the adjustable factors [32,33,34]. A hybrid tuning approach was created based on the genetic algorithm (GA), particle swarm optimization (PSO), and neural network algorithm (NN) for coefficient gear fault detection [35]. In [36,37,38], another probabilistic neural network was utilized to solve the tuning issue based on the NN procedure in [39]. Entrapment at local optima represents a big issue of these optimization approaches [40,41]. Different optimization algorithms have been applied to different applications to tackle this issue [42,43,44]. In [45], parametric tuning of a CNC-drilling system was performed based on the Taguchi–whale optimization algorithm. The Divide-and-Conquer Bat Algorithm was applied to tune cutting parameters in CNC Turnings in [46]. In [47], a quadratic interpolation approach was introduced with a whale optimization algorithm to solve high-dimensional global optimization problems. A machine learning technique was combined with a constrained coral reef optimization algorithm to tune multi-reservoir processing [48]. In [49], a local escaping operator and orthogonal learning were utilized to improve the Archimedes optimization algorithm for PEM fuel cell parameter tuning. Among these algorithms, the equilibrium optimizer (EO) approach has proved to be an effective solution for tuning issues in various applications [50,51,52,53,54]. This algorithm can tackle the local optimum trapping issue and demonstrates effective tuning with a fast convergence rate and few parameters. This paper introduces the EO for the optimization of the cutting force factors of a milling machine in place of the traditional approaches that depend on the trial-and-error method of the designer. The proposed EOA utilizes a few adjustable factors to improve the characteristics of the cutting process, performs tuning procedures with a high-speed convergence rate, and tackles the local optimum trapping issue instead of other approaches. Furthermore, the performance of the milling machine based on the proposed algorithm is compared with the traditional method for the cutting force factors. Various experiments were conducted to confirm the effectiveness of the developed approaches under different cutting conditions. The accomplishments of this research work are summarized as follows:
  • A new tuning approach based on EO is introduced to improve the cutting characteristics of milling machines;
  • The tuning issue of the cutting force coefficients is tackled based on the developed EO instead of the conventional approaches that depend on the trial-and-error method of the designer;
  • The introduced EO approach can improve the cutting conditions with few adjustable factors and overcome the local optimum trapping issue;
  • The integral square error (ISE) index is utilized to evaluate the performance of the milling machine based on the proposed algorithm as compared with the traditional method of cutting force factors and genetic algorithm (GA);
  • The proposed EO has a minimum ISE of around 1.12, while the genetic algorithm (GA) has an ISE of around 1.14 and the trial-and-error method has an ISE of around 2.4;
  • The experimental tests confirm the effectiveness of the developed approaches under different cutting conditions.
The rest of the paper is organized as follows: Section 2 characterizes the chatter vibration issue in the milling process. In Section 3, the equilibrium optimizer approach is derived. Section 4 presents the experimental results and related discussion. Finally, the main conclusions are visualized.

2. Chatter Vibration Phenomenon in the Milling Process

Figure 1 presents a model of the end mill including a milling cutter around two orthogonal degrees of freedom in the feed direction (x) and the normal direction (y). The tool vibration in the x and y directions is projected into the surface normal direction to evaluate the instantaneous chip thickness. The tool vibration in the normal direction is denoted by n.
n = x . sin θ j y . cos θ j
The dynamic displacements x and y are caused when the cutting force excites the structure in those two directions. Therefore, the instantaneous chip thickness for the ith disk element, the jth tooth, and the kth angular position is rewritten in Equation (2) [2]:
h θ j = f t sin θ j + n t τ n t g θ j
where f t is the feed per tooth. The tooth period is τ = 60 / Ω N t . For the end mill, the instantaneous immersion angle, θ j , is determined by considering the lag angle δ .
θ j = ϕ j + 2 π k 1 N j δ
δ = z r tan γ
ϕ j = 6 Ω t
The term [n(tτ) − n(t)] represents the dynamic chip thickness produced owing to vibrations at the present time t and one spindle revolution period prior. The switching function g θ j describes the engagement between each tooth and the workpiece in one revolution, as shown in Equation (6).
g θ i , j , k = 1 ,   θ s θ j θ e   0 ,                     o t h e r w i s e
The instantaneous chip thickness from Equation (2) can also be rewritten using the x- and y-dynamic displacements as:
h θ j = Δ x sin θ j + Δ y cos θ j g θ j
where Δ x = x x 0   and   Δ y = y y 0 are the dynamic displacements of the cutter in the x and y directions, respectively. The inner and outer modulations are represented by x , y and x 0 , y 0 .
Assuming that the milling cutter is approximated as a system with two orthogonal degrees of freedom in the feed direction (x) and the normal direction (y), the cutting force excites the structure in the x and y directions and causes the dynamic displacements x and y. The equations of motion in the x and y directions, including process damping, are expressed in Equation (8). Here, total damping in the model is assumed, including structural damping and process damping:
m x ¨ + c + c d , x x ˙ + k x = F x m y ¨ + c + c d , y y ˙ + k y = F y
where c d , x c d , y = cos ϕ j sin ϕ j sin ϕ j cos ϕ j c d t c d n .
The damping coefficient in the normal direction c d n and that in the tangential direction c d t are determined as follows:
c d n = K d b π A   π / 2 ω c 3 π / 2 ω c S d cos ( ω c t ) d t   and   c d t = μ c d n
in which μ presents the coefficient of contact friction between the flank face and the machined surface. The contact friction coefficients can be determined based on the workpiece and tool material.
Furthermore, the equations of dynamic forces can be converted to (x, y) coordinates in matrix form:
F x F y = 1 2 b K t a x x a x y a y x a y y Δ x Δ y = 1 2 b K t A Δ x Δ y
where b is the axial depth of the cut, Kt is the tangential cutting force coefficient, and matrix [A] is the directional coefficients relating the dynamic displacements to the dynamic cutting forces. These directional dynamic force coefficients are calculated via Equations (11)–(14).
α x x = j = 1 N t g ( θ j ) sin ( 2 θ j ) + K n 1 cos ( 2 θ j )
α x y = j = 1 N t g ( θ j ) 1 + cos ( 2 θ j ) + K n sin ( 2 θ j )
α y x = j = 1 N t g ( θ j ) 1 cos ( 2 θ j ) K n sin ( 2 θ j )
α y y = j = 1 N t g ( θ j ) sin ( 2 θ j ) K n 1 + cos ( 2 θ j )
Because the matrix [ A ( θ ) ] can be expressed as a function of time [ A ( t ) ] , the term [ A ( t ) ] can then be expanded into a Fourier series with Fourier coefficients calculated via Equation (15).
A ( t ) = A r e j r ω t o o t h t
A r = 1 N t τ 0 N t τ A ( t ) e j r ω t o o t h t d t
The higher-order terms of the Fourier series are neglected. Then, the dynamic cutting force becomes:
F ( t ) = F x F y = 1 2 b K t A 0 Δ x Δ y
A 0 = N t 2 π α x x α x y α y x α y y
The vibration response in the frequency domain can be represented as Equation (19).
x ˜ y ˜ = Φ x x ( j ω ) 0 0 Φ y y ( j ω ) F ˜ x F ˜ y  
where Φ x x ( j ω ) and Φ y y ( j ω ) are direct transfer functions in the x and y directions. The characteristic equation can be written as in Equation (20).
det I + Λ Φ o r = 0
Λ = N t 4 π b K t ( 1 e j ω c τ )
Φ o r = α x x Φ x x α x y Φ y y α y x Φ x x α y y Φ y y
This results in a quadratic equation, as shown in Equation (23).
Φ x x Φ y y α x x α y y α x y α y x Λ 2 + α x x Φ x x + α y y Φ y y Λ + 1 = 0
The limiting chip width or critical depth of the cut is evaluated from the real and imaginary parts of the eigenvalue Λ = Λ Re + j Λ Im , and its relationship with Ω is expressed by Equation (24).
b lim = 2 π Λ Re N t K t 1 + Λ Im Λ Re 2
The phase shift between subsequent tooth passages, tooth passing periods, and spindle speeds vs. blim is determined as shown in Equations (25)–(27).
ε p h = π 2 μ ;   μ = a tan Λ Im Λ Re
τ = 1 ω c ε p h + j 2 π
Ω = 60 N t τ
The stability lobe diagram is generated by plotting spindle speeds vs. two limiting chip width values for each chatter frequency.

3. Equilibrium Optimizer Approach

This paper proposes the use of the EO as a recent optimization algorithm that requires few adjustable parameters. Furthermore, cooperation between multiple agents is the focus in EO to improve the exploration manner, supporting global search and avoiding the issue of entrapment in a particular local optimum [52]. Equilibrium Optimizer is a physics-based algorithm that imitates the dynamic mass stabilization behavior within the control volume. The concentration of nonreactive constituents is described as a function of various source and sink dynamics based on the mass balance formulation. The generic mass stabilization is formulated based on a first-order differential equation to describe the variation in the dynamic system due to the mass entering and mass leaving as follows:
V d C d t = F C e q F + G
where V is the control volume, C is the concentration, F is the flow rate, Ceq is the concentration of the equilibrium state, and G is the generation rate of mass.
The EO algorithm consists of different stages to find the best solutions, and it is defined as follows.
Stage 1: Initialization
The EO starts the initial populations randomly, similarly to other population-based metaheuristics. In this stage, the initial concentrations of each particle are determined randomly within the limits of the system gains as follows:
C i = C min + r ( C min C max )
where Cmin is the minimum concentrations, Cmax is the maximum concentrations, r is random vector within [0, 1], and i is the population index.
After the determination of the concentration of each particle, the fitness function of the optimization problem is evaluated for each particle. Then, the fitness values are sorted to clarify the equilibrium candidates. In the EO algorithm, the four best particles and the mean of these candidates are chosen to represent the equilibrium pool in order to enhance exploration and exploitation behaviors.
Stage 2: Equilibrium pool construction
In this stage, the equilibrium pool is constructed based on the best four particles and their average from the above stage as follows:
C p o o l = C e q ( 1 ) , C e q ( 2 ) , C e q ( 3 ) , C e q ( 4 ) , C e q ( a v g )
Each particle within the population updates its concentration during the iteration’s scope based on each candidate in the selected equilibrium pool with the same probability. In addition, each particle utilizes the same number of updates within the updating process from all of the candidate solutions until the end of the optimization process.
Stage 3: Balancing exploration and exploitation
Obtaining a balance between exploration and exploitation is an essential procedure in optimization problems. The EO algorithm utilizes the exponential term ‘E’ to provide a proper balance between exploration and exploitation. This exponential term ‘E’ is defined as follows:
E = exp ( δ ( t t 0 ) )
where δ is a random factor within [0, 1], while the factors t and t0 are defined as follows:
t = 1 I t e r I t e r max a 2 × I t e r I t e r max
t 0 = 1 δ ln a 1 s i g n ( r 0.5 [ 1 e δ t ] ) + t
We substitute t and t0 into Equation (43); then, the exponential term ‘E’ can be described as:
E = a 1 s i g n ( r 0.5 ) [ e δ t 1 ]
The factors a1 and a2 control the behavior of the exponential term, leading to balance in the exploration and exploitation of the EO algorithm. The term s i g n ( r 0.5 ) controls the direction of exploration and exploitation within the optimization process.
Stage 4: Convergence to the optimal global solution
In this stage, convergence to the optimal global solution is carried out by utilizing an exponential generation rate factor that is described as follows:
G = G 0 × exp ( δ ( t t 0 ) ) = G 0 × E
where G0 is the initial value of the generation rate factor and is defined as:
G 0 = β ( C e q δ × C )
Here, β is a controlling factor for the generation rate and is utilized to adjust the exploitation and exploration of each particle as follows:
β = 0.5 r 1       r 2 > ρ 0       r 1 < ρ
where r1 and r2 are uniformly distributed random numbers in [0, 1], while ρ is a probability factor. The above stages are performed every iteration to find the best solutions until the stopping criteria are achieved and the best solution for the optimization problem is identified.

4. Results and Discussion

A series of experiments were conducted on a three-axis CNC milling machine (with a Heidenhain TNC620 controller) to determine the specific cutting force coefficients and cutting force angle for the milling force model. The workpiece was a block of Al6061-T6, which is commonly utilized in the automobile and aerospace industries because of its high strength-to-weight ratio. An end mill cutter with a diameter of 12 mm, helix angle of 26°, and two flutes was used. Besides this, the tooth-to-tooth radius error was 6 lm, and the cutting force signals were measured by a Kistler dynamometer mounted between the workpiece and workbench. In this study, lot milling tests were conducted with start cutting angle ϕs = 0 and existing cutting angle ϕe = 180°. The average cutting force in the x and y directions, in this case, is described in Equations (38) and (39).
F ¯ x = N t b K n 4 f t + N t b K n e π
F ¯ y = N t b K t 4 f t + N t b K t e π
Six slot milling tests corresponding to different feed rates were conducted to achieve average cutting force in the x and y directions. The experimental parameters are provided in Table 1. Figure 2 shows the measured cutting forced signal at a spindle speed of 2500 rpm, cut axial depth of 1 mm, and feed rate of 0.03 mm/tooth in the (a) x-direction and (b) y-direction.
The linear regression method summarized by Schmitz [23] was then utilized to determine the four unknown cutting force coefficients in Equations (38) and (39). Those values can be obtained by Equations (40) and (41):
K n = 4 a 1 x N t b ;         K n e = π a 0 x N t b
K t = 4 a 1 y N t b ;         K t e = π a 0 y N t b
where K n , K t are the radial and tangential cutting coefficients, respectively, and K n e , K t e are the radial and tangential edge coefficients, respectively.
a 1 x = i = 1 n f t , i F ¯ x , i i = 1 n f t , i i = 1 n F ¯ x , i n i = 1 n f t , i 2 i = 1 n f t , i 2
a 0 x = 1 n i = 1 n F ¯ x , i a 1 x 1 n i = 1 n f t , i
a 1 y = i = 1 n f t , i F ¯ y , i i = 1 n f t , i i = 1 n F ¯ y , i n i = 1 n f t , i 2 i = 1 n f t , i 2
a 0 y = 1 n i = 1 n F ¯ y , i a 1 y 1 n i = 1 n f t , i
The results of linear regression are shown in Figure 3. The four unknown cutting force coefficients are presented in Table 2. The specific cutting force coefficient K s and cutting force angle β were then calculated following Equation (46).
K s = K t 1 + K n β = a tan 1 K n
Once the FRFs in the x and y directions were measured, a model was defined by performing a modal fit to the measured data. To identify the modal parameters, the fitting approach was a peak-picking method wherein we used the real and imaginary parts of the system FRFs. This work was done on TXFTM software, and the model fit results are shown in Figure 4a,b in which six modes were selected in both the x and y directions.
The peak values of real/imaginary parts were selected, and the corresponding values of frequencies in the x and y directions were applied. Consequently, the model parameters were calculated by using Equations (47)–(49).
ζ q i = ω r e a l min i ω r e a l max i 2 ω n i
k q i = 1 min Im F R F i 2 ζ q i
m q i = k q i ω n i 2
Figure 5 illustrates the procedure for optimization of dynamic cutting force coefficients for the end-milling process based on the EO algorithm. The first step was to determine the modal parameters via modal analysis, then the instantaneous chip thickness was calculated using Equation (2). The dynamic cutting force coefficients were generated by the linear regression method. The proposed dynamic cutting force model of milling was used to estimate the cutting forces. Finally, the residual cutting force between the prediction and the actual cutting force was considered as an objective function for optimizing the cutting force coefficients using the EO approach. The boundary of the parameters was adjusted as lower bound = [600 × 106 50 10 × 103 5 × 103 1 × 103] and upper bound = [1000 × 106 90 50 × 103 30 × 103 20 × 103]. As a result, the optimum set of specific cutting force coefficients by the EO method was obtained, as presented in Table 3.
The simulation cutting forces under various cutting conditions are shown in Figure 6, in which the stable cutting condition in Figure 6a indicates that the cutting force behavior is repetitive from one to the next cutting revolution. On the other hand, in Figure 6b, the repetitive behavior of the cutting force does not remain from one to the next cutting revolution of the cutter; instead, the magnitude of the cutting force eventually increases during the cutting.
Figure 7 shows the performance of the cutting force model when experimentally validated by the measured cutting force. A comparison between the simulated and experimental results for both stable and unstable cutting conditions is presented in Figure 7a,b, respectively. The simulated results coincide with the measured cutting forces. Moreover, this study provides a new solution for the tuning of the force factors of the milling process based on an intelligent algorithm named EO instead of trial-and-error methods. The EO algorithm finds the best force factors based on the minimization of the integral square error (ISE), defined as follows:
I S E = e r 2 d t
where
e r = F x F x
As the result, the errors of the cutting forces of the proposed dynamic cutting force model with and without the EO method are illustrated in Figure 8. The figure shows that the performance of the cutting force model is improved when applying the EO algorithm to correct the specific cutting force coefficients. The errors between simulated and measured cutting forces are significantly reduced with the EO method in both stable and unstable cutting conditions. Moreover, the integral square error values based on trial and error, GA, and the proposed EO method are shown in Figure 9 and Table 4. It is clear that the proposed EO has a minimum ISE of around 1.12, while the GA has an ISE of around 1.14 and the trial-and-error method has an ISE of around 2.4. The ISE values of the cutting force model are significantly decreased when correcting the cutting force coefficients by the EO algorithm.

5. Conclusions

During the milling process, the engagement between the tool and the workpiece relatively oscillates and affects the cutting depth. This might include both cutting force vibration and cutting chatter vibration and creates difficulty in modeling the cutting force. In this study, we proposed a dynamic cutting force model of the milling process, in which the modal parameters were determined via impact testing and modal analysis. The dynamic cutting force coefficients were generated by the linear regression method. The proposed dynamic cutting force model of milling was used to estimate the cutting forces. Finally, the square error of the residual cutting force between the prediction and the actual cutting force was considered as an objective function to optimize the cutting force coefficients using the EO approach. The linear regression method was firstly used to determine the specific cutting force coefficients, and then the EO algorithm was adopted to optimize the set of specific cutting force coefficients based on minimization of the integral square error between the simulated and measured cutting forces. The performance of the proposed force model was experimentally validated under various cutting conditions. The results indicate that the cutting force model was significantly improved by applying the EO algorithm to correct the specific cutting force coefficients. The proposed EO had a minimum ISE of around 1.12, while the GA had an ISE of around 1.14 and the trial-and-error method had an ISE of around 2.4. The errors between the simulated and measured cutting forces were also significantly reduced by the EO method under both stable and unstable cutting conditions.

Author Contributions

Conceptualization, M.E. and M.-Q.T.; Data curation, M.E. and M.-Q.T.; Formal analysis, M.E., M.-Q.T. and V.Q.V.; Funding acquisition, F.A. and S.S.M.G.; Investigation, M.-Q.T., V.Q.V., F.A. and S.S.M.G.; Methodology, M.E. and M.-Q.T.; Project administration, F.A.; Resources, M.E. and F.A.; Software, M.-Q.T.; Supervision, M.E. and S.S.M.G.; Validation, V.Q.V., F.A. and S.S.M.G.; Visualization, V.Q.V., F.A. and S.S.M.G.; Writing—original draft, M.E. and M.-Q.T.; Writing—review and editing, V.Q.V., F.A. and S.S.M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Taif University Researchers Supporting Project number (TURSP-2020/97), Taif University, Taif, Saudi Arabia and in part by the Ministry of Science and Technology (MOST) of Taiwan (grant numbers: MOST 110-2222-E-011-002- and MOST 110-2222-E-011-013-).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Taif University Researchers Supporting Project number (TURSP-2020/97), Taif University, Taif, Saudi Arabia for supporting this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pereira, R.B.D.; Brandão, L.C.; de Paiva, A.P.; Ferreira, J.R.; Davim, J.P. A review of helical milling process. Int. J. Mach. Tools Manuf. 2017, 120, 27–48. [Google Scholar] [CrossRef]
  2. Altintas, Y. Manufacturing Automation: Metal Cutting Mechanics, Machine Tool Vibrations, and CNC Design, 2nd ed.; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  3. Taner, T.L.; Ömer, Ö.; Erhan, B. Generalized cutting force model in multi-Axis milling using a new engagement boundary determination approach. Int. J. Adv. Manuf. Technol. 2015, 77, 341–355. [Google Scholar] [CrossRef]
  4. Chung, C.; Tran, M.-Q.; Liu, M.-K. Estimation of Process Damping Coefficient Using Dynamic Cutting Force Model. Int. J. Precis. Eng. Manuf. 2020, 21, 623–632. [Google Scholar] [CrossRef]
  5. Yue, C.; Gao, H.; Liu, X.; Liang, S.Y.; Wang, L. A review of chatter vibration research in milling. Chin. J. Aeronaut. 2019, 32, 215–242. [Google Scholar] [CrossRef]
  6. Moon, F.C.; Kalmár-Nagy, T. Nonlinear Models for Complex Dynamics in Cutting Materials. Philos. Trans. Math. Phys. Eng. Sci. 2001, 359, 695–711. [Google Scholar] [CrossRef]
  7. Kasahara, N.; Sato, H.; Tani, Y. Phase Characteristics of Self-Excited Chatter in Cutting. J. Eng. Ind. 1992, 114, 393–399. [Google Scholar] [CrossRef]
  8. Munoa, J.; Beudaert, X.; Dombovari, Z.; Altintas, Y.; Budak, E.; Brecher, C.; Stepan, G. Chatter suppression techniques in metal cutting. CIRP Ann. 2016, 65, 785–808. [Google Scholar] [CrossRef]
  9. Zengxi, P.; Hui, Z. Analysis and suppression of chatter in robotic machining process. In Proceedings of the 2007 International Conference on Control, Automation and Systems, Seoul, Korea, 17–20 October 2007; pp. 595–600. [Google Scholar] [CrossRef]
  10. Ding, L.; Sun, Y.; Xiong, Z. Active Chatter Suppression in Turning by Simultaneous Adjustment of Amplitude and Frequency of Spindle Speed Variation. J. Manuf. Sci. Eng. 2019, 142, 021004. [Google Scholar] [CrossRef]
  11. Altintas, Y.; Stepan, G.; Budak, E.; Schmitz, T.; Kilic, Z.M. Chatter Stability of Machining Operations. J. Manuf. Sci. Eng. 2020, 142, 110801. [Google Scholar] [CrossRef]
  12. Tran, M.-Q.; Liu, M.-K.; Elsisi, M. Effective multi-Sensor data fusion for chatter detection in milling process. ISA Trans. 2021, 125, 514–527. [Google Scholar] [CrossRef]
  13. Mohammadi, Y.; Ahmadi, K. Frequency domain analysis of regenerative chatter in machine tools with Linear Time Periodic dynamics. Mech. Syst. Signal Process. 2019, 120, 378–391. [Google Scholar] [CrossRef]
  14. Mou, W.; Zhu, S.; Jiang, Z.; Song, G. Vibration signal-Based chatter identification for milling of thin-Walled structure. Chin. J. Aeronaut. 2022, 35, 204–214. [Google Scholar] [CrossRef]
  15. Perrelli, M.; Cosco, F.; Gagliardi, F.; Mundo, D. In-Process Chatter Detection Using Signal Analysis in Frequency and Time-Frequency Domain. Machines 2022, 10, 24. [Google Scholar] [CrossRef]
  16. Tran, M.-Q.; Elsisi, M.; Liu, M.-K. Effective feature selection with fuzzy entropy and similarity classifier for chatter vibration diagnosis. Measurement 2021, 184, 109962. [Google Scholar] [CrossRef]
  17. Yoon, M.C.; Chin, D.H. Cutting force monitoring in the endmilling operation for chatter detection. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2005, 219, 455–465. [Google Scholar] [CrossRef]
  18. Sun, Y.; Jiang, S. Predictive modeling of chatter stability considering force-Induced deformation effect in milling thin-Walled parts. Int. J. Mach. Tools Manuf. 2018, 135, 38–52. [Google Scholar] [CrossRef]
  19. Zhang, X.; Zhang, J.; Pang, B.; Zhao, W. An accurate prediction method of cutting forces in 5-Axis flank milling of sculptured surface. Int. J. Mach. Tools Manuf. 2016, 104, 26–36. [Google Scholar] [CrossRef]
  20. Ji, Y.; Bi, Q.; Zhang, S.; Wang, Y. A new receptance coupling substructure analysis methodology to predict tool tip dynamics. Int. J. Mach. Tools Manuf. 2018, 126, 18–26. [Google Scholar] [CrossRef]
  21. Zhang, J.; Li, J.; Xie, Z.; Du, C.; Gui, L.; Zhao, W. Rapid dynamics prediction of tool point for bi-Rotary head five-Axis machine tool. Precis. Eng. 2017, 48, 203–215. [Google Scholar] [CrossRef]
  22. Liu, M.-K.; Tran, M.-Q.; Chung, C.; Qui, Y.-W. Hybrid model- and signal-Based chatter detection in the milling process. J. Mech. Sci. Technol. 2020, 34, 1–10. [Google Scholar] [CrossRef]
  23. Smith, K.S.; Schmitz, T.L. Machining Dynamics: Frequency Response to Improved Productivity; Springer International Publishing: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  24. Tran, M.-Q.; Liu, M.-K.; Tran, Q.-V. Analysis of Milling Chatter Vibration Based on Force Signal in Time Domain. In Advances in Engineering Research and Application; Sattler, K.-U., Nguyen, D.C., Vu, N.P., Long, B.T., Puta, H., Eds.; Springer International Publishing: Berlin/Heidelberg, Germany, 2021; pp. 192–199. [Google Scholar]
  25. Liu, M.-K.; Tran, Q.M.; Qui, Y.-W.; Chung, C.-H. Chatter Detection in Milling Process Based on Time-Frequency Analysis. In Proceedings of the ASME 2017 12th International Manufacturing Science and Engineering Conference Collocated with the JSME/ASME 2017 6th International Conference on Materials and Processing, Los Angeles, CA, USA, 4–8 June 2017; Volume 1, p. V001T02A025. [Google Scholar] [CrossRef]
  26. Liu, C.; Zhu, L.; Ni, C. The chatter identification in end milling based on combining EMD and WPD. Int. J. Adv. Manuf. Technol. 2017, 91, 3339–3348. [Google Scholar] [CrossRef]
  27. Gao, G.; Wu, B.; Zhang, D.; Luo, M. Mechanistic identification of cutting force coefficients in bull-Nose milling process. Chin. J. Aeronaut. 2013, 26, 823–830. [Google Scholar] [CrossRef]
  28. Adem, K.A.M.; Fales, R.; El-Gizawy, A.S. Identification of cutting force coefficients for the linear and nonlinear force models in end milling process using average forces and optimization technique methods. Int. J. Adv. Manuf. Technol. 2015, 79, 1671–1687. [Google Scholar] [CrossRef]
  29. Budak, E.; Altintas, Y.; Armarego, E.J.A. Prediction of Milling Force Coefficients from Orthogonal Cutting Data. J. Manuf. Sci. Eng. 1996, 118, 216–224. [Google Scholar] [CrossRef]
  30. Chen, X.; Zhang, D.; Wang, Q. Cutting Force Transition Model Considering the Influence of Tool System by Using Standard Test Table. Sensors 2021, 21, 1340. [Google Scholar] [CrossRef] [PubMed]
  31. Lamikiz, A.; de Lacalle, L.N.L.; Sanchez, J.A.; Bravo, U. Calculation of the specific cutting coefficients and geometrical aspects in sculptured surface machining. Mach. Sci. Technol. 2005, 9, 411–436. [Google Scholar] [CrossRef]
  32. Elsisi, M.; Ebrahim, M.A. Optimal design of low computational burden model predictive control based on SSDA towards autonomous vehicle under vision dynamics. Int. J. Intell. Syst. 2021, 36, 6968–6987. [Google Scholar] [CrossRef]
  33. Halim, A.H.; Ismail, I.; Das, S. Performance assessment of the metaheuristic optimization algorithms: An exhaustive review. Artif. Intell. Rev. 2021, 54, 2323–2409. [Google Scholar] [CrossRef]
  34. Elsisi, M.; Tran, M.Q.; Hasanien, H.M.; Turky, R.A.; Albalawi, F.; Ghoneim, S.S. Robust model predictive control paradigm for automatic voltage regulators against uncertainty based on optimization algorithms. Mathematics 2021, 9, 2885. [Google Scholar] [CrossRef]
  35. Ding, J.; Xiao, D.; Li, X. Gear fault diagnosis based on genetic mutation particle swarm optimization VMD and probabilistic neural network algorithm. IEEE Access 2020, 8, 18456–18474. [Google Scholar] [CrossRef]
  36. Seyfipour, N.; Menhaj, M.; Nik, R. A New Optimization Method by Ring Probabilistic Logic Neural Networks. AMIRKABIR 2003, 14, 43–57. [Google Scholar]
  37. Azizi, A.; Vatankhah Barenji, A.; Hashmipour, M. Optimizing radio frequency identification network planning through ring probabilistic logic neurons. Adv. Mech. Eng. 2016, 8, 1687814016663476. [Google Scholar] [CrossRef]
  38. Menhaj, M.B.; Seifipour, N. Function optimization by RPLNN. In Proceedings of the 2002 International Joint Conference on Neural Networks, IJCNN’02 (Cat. No. 02CH37290), Honolulu, HI, USA, 12–17 May 2002; Volume 2, pp. 1522–1527. [Google Scholar]
  39. Qu, M.; Tang, J. Probabilistic logic neural networks for reasoning. arXiv 2019, arXiv:1906.08495. [Google Scholar]
  40. Guo, X.; Albalawi, F.; Laleg-Kirati, T.-M. Observer-based economic model predictive control for direct contact membrane distillation. Chem. Eng. Res. Des. 2020, 156, 86–99. [Google Scholar] [CrossRef]
  41. Cao, Y.; Zhang, H.; Li, W.; Zhou, M.; Zhang, Y.; Chaovalitwongse, W.A. Comprehensive learning particle swarm optimization algorithm with local search for multimodal functions. IEEE Trans. Evol. Comput. 2018, 23, 718–731. [Google Scholar] [CrossRef]
  42. Alanqar, A.; Durand, H.; Albalawi, F.; Christofides, P.D. An economic model predictive control approach to integrated production management and process operation. AIChE J. 2017, 63, 1892–1906. [Google Scholar] [CrossRef]
  43. Li, W.; Wang, G.G.; Gandomi, A.H. A survey of learning-Based intelligent optimization algorithms. Arch. Comput. Methods Eng. 2021, 28, 3781–3799. [Google Scholar] [CrossRef]
  44. Albalawi, F.; Durand, H.; Christofides, P.D. Distributed economic model predictive control with Safeness-Index based constraints for nonlinear systems. Syst. Control Lett. 2017, 110, 21–28. [Google Scholar] [CrossRef]
  45. Sahoo, A.K.; Jeet, S.; Bagal, D.K.; Barua, A.; Pattanaik, A.K.; Behera, N. Parametric optimization of CNC-drilling of Inconel 718 with cryogenically treated drill-Bit using Taguchi-Whale optimization algorithm. Mater. Today Proc. 2022, 50, 1591–1598. [Google Scholar] [CrossRef]
  46. Huang, X.; He, Z.; Chen, Y.; Xie, S. A Divide-And-Conquer Bat Algorithm with Direction of Mean Best Position for Optimization of Cutting Parameters in CNC Turnings. Comput. Intell. Neurosci. 2022, 2022, 4719266. [Google Scholar] [CrossRef]
  47. Sun, Y.; Yang, T.; Liu, Z. A whale optimization algorithm based on quadratic interpolation for high-Dimensional global optimization problems. Appl. Soft Comput. 2019, 85, 105744. [Google Scholar] [CrossRef]
  48. Emami, M.; Nazif, S.; Mousavi, S.F.; Karami, H.; Daccache, A. A hybrid constrained coral reefs optimization algorithm with machine learning for optimizing multi-Reservoir systems operation. J. Environ. Manag. 2021, 286, 112250. [Google Scholar] [CrossRef] [PubMed]
  49. Houssein, E.H.; Helmy, B.E.D.; Rezk, H.; Nassef, A.M. An enhanced Archimedes optimization algorithm based on Local escaping operator and Orthogonal learning for PEM fuel cell parameter identification. Eng. Appl. Artif. Intell. 2021, 103, 104309. [Google Scholar] [CrossRef]
  50. Elsheikh, A.H.; Shehabeldeen, T.A.; Zhou, J.; Showaib, E.; Abd Elaziz, M. Prediction of laser cutting parameters for polymethylmethacrylate sheets using random vector functional link network integrated with equilibrium optimizer. J. Intell. Manuf. 2020, 32, 1377–1388. [Google Scholar] [CrossRef]
  51. Yu, H.; Yang, B.; Wang, S.; Wang, Y.; Wang, S.; Wang, Z.; Wang, Z. An effective multi-Part dedicatedflow-Line reconfiguration model considering the optimal selection of machining process path and machines. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2022, 09544054221100873. [Google Scholar] [CrossRef]
  52. Faramarzi, A.; Heidarinejad, M.; Stephens, B.; Mirjalili, S. Equilibrium optimizer: A novel optimization algorithm. Knowl.-Based Syst. 2020, 191, 105190. [Google Scholar] [CrossRef]
  53. Gao, Y.; Zhou, Y.; Luo, Q. An efficient binary equilibrium optimizer algorithm for feature selection. IEEE Access 2020, 8, 140936–140963. [Google Scholar] [CrossRef]
  54. Houssein, E.H.; Dirar, M.; Abualigah, L.; Mohamed, W.M. An efficient equilibrium optimizer with support vector regression for stock market prediction. Neural Comput. Appl. 2022, 34, 3165–3200. [Google Scholar] [CrossRef]
Figure 1. Dynamic cutting force model of the milling process.
Figure 1. Dynamic cutting force model of the milling process.
Mathematics 10 03287 g001
Figure 2. Measured cutting force signal at a feed rate of 0.03 mm/tooth in the (a) x-direction and (b) y-direction.
Figure 2. Measured cutting force signal at a feed rate of 0.03 mm/tooth in the (a) x-direction and (b) y-direction.
Mathematics 10 03287 g002
Figure 3. Relationship between feed per tooth and cutting force.
Figure 3. Relationship between feed per tooth and cutting force.
Mathematics 10 03287 g003
Figure 4. Relative frequency responses in the (a) x-direction and (b) y-direction.
Figure 4. Relative frequency responses in the (a) x-direction and (b) y-direction.
Mathematics 10 03287 g004
Figure 5. Procedure for optimization of dynamic cutting force coefficients for the end-milling process.
Figure 5. Procedure for optimization of dynamic cutting force coefficients for the end-milling process.
Mathematics 10 03287 g005
Figure 6. Simulation of cutting force under different cutting conditions: (a) stable cutting: 3750 rpm spindle and 0.4 mm cut depth; (b) unstable cutting: 3750 rpm spindle and 1.4 mm cut depth.
Figure 6. Simulation of cutting force under different cutting conditions: (a) stable cutting: 3750 rpm spindle and 0.4 mm cut depth; (b) unstable cutting: 3750 rpm spindle and 1.4 mm cut depth.
Mathematics 10 03287 g006
Figure 7. Experimental validation of the cutting force model: (a) stable cutting: 2250 rpm spindle and 1.4 mm cut depth; (b) unstable cutting: 3750 rpm spindle and 1.4 mm cut depth.
Figure 7. Experimental validation of the cutting force model: (a) stable cutting: 2250 rpm spindle and 1.4 mm cut depth; (b) unstable cutting: 3750 rpm spindle and 1.4 mm cut depth.
Mathematics 10 03287 g007
Figure 8. Residual values of the proposed dynamic cutting force model with and without the EO method: (a) stable cutting: 1500 rpm spindle and 1.4 mm cut depth; (b) stable cutting: 2250 rpm spindle and 1.4 mm cut depth; (c) unstable cutting: 3750 rpm spindle and 1.4 mm cut depth.
Figure 8. Residual values of the proposed dynamic cutting force model with and without the EO method: (a) stable cutting: 1500 rpm spindle and 1.4 mm cut depth; (b) stable cutting: 2250 rpm spindle and 1.4 mm cut depth; (c) unstable cutting: 3750 rpm spindle and 1.4 mm cut depth.
Mathematics 10 03287 g008
Figure 9. Integral square error values from simulation with and without the EO method.
Figure 9. Integral square error values from simulation with and without the EO method.
Mathematics 10 03287 g009
Table 1. Cutting parameters for cutting force coefficients.
Table 1. Cutting parameters for cutting force coefficients.
Spindle Speed
(rpm)
Axial Depth of Cut
(mm)
Radial Depth of Cut
(mm)
Feed Rate
(mm/min)
25001.012100, 150, 175, 200, 225, 250
Table 2. Specific cutting force coefficients.
Table 2. Specific cutting force coefficients.
KsβKneKte
N/mmdegN/mmN/mm
78065.3212.3224.9
Table 3. The optimum set of specific cutting force coefficients by EO.
Table 3. The optimum set of specific cutting force coefficients by EO.
KsβKneKte
N/mmdegN/mmN/mm
60050.6828.8518.01
Table 4. ISE values in the case of the proposed EO and different techniques in literature.
Table 4. ISE values in the case of the proposed EO and different techniques in literature.
Trial-and-Error MethodGAProposed EO
ISE2.41.141.12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tran, M.-Q.; Elsisi, M.; Vu, V.Q.; Albalawi, F.; Ghoneim, S.S.M. Effective Optimization Based on Equilibrium Optimizer for Dynamic Cutting Force Coefficients of the End-Milling Process. Mathematics 2022, 10, 3287. https://doi.org/10.3390/math10183287

AMA Style

Tran M-Q, Elsisi M, Vu VQ, Albalawi F, Ghoneim SSM. Effective Optimization Based on Equilibrium Optimizer for Dynamic Cutting Force Coefficients of the End-Milling Process. Mathematics. 2022; 10(18):3287. https://doi.org/10.3390/math10183287

Chicago/Turabian Style

Tran, Minh-Quang, Mahmoud Elsisi, Viet Q. Vu, Fahad Albalawi, and Sherif S. M. Ghoneim. 2022. "Effective Optimization Based on Equilibrium Optimizer for Dynamic Cutting Force Coefficients of the End-Milling Process" Mathematics 10, no. 18: 3287. https://doi.org/10.3390/math10183287

APA Style

Tran, M. -Q., Elsisi, M., Vu, V. Q., Albalawi, F., & Ghoneim, S. S. M. (2022). Effective Optimization Based on Equilibrium Optimizer for Dynamic Cutting Force Coefficients of the End-Milling Process. Mathematics, 10(18), 3287. https://doi.org/10.3390/math10183287

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