# An Optimized PID Controller Using Enhanced Bat Algorithm in Drilling Processes

<figure><img src="image_1.png" /></figure>

Check for
updates

Rash M. Naji*<sup>1,2</sup>, Hussien Dulaimi<sup>1,2</sup>, Huthaifa Al-Khazraji<sup>1,2</sup>

Control and System Engineering Department, University of Technology- Iraq, Baghdad 10066, Iraq

Corresponding Author Email: 60049@uotechnology.edu.iq

Copyright: ©2024 The authors. This article is published by IIETA and is licensed under the CC BY 4.0 license (http://creativecommons.org/licenses/by/4.0/).

https://doi.org/10.18280/jesa.570314

**Received:** 20 March 2024

**Revised:** 20 May 2024

**Accepted:** 7 June 2024

**Available online:** 25 June 2024

**Keywords:**

drilling process, production, PID controller,
swam optimization, enhance bat algorithm

## ABSTRACT

Drilling operation has a direct impact on the quality of the production. Insufficiently controlling the cutting force during the drilling process leads to the risk of early drill failure. Typically, the selection of the drilling parameters is determined based on machining-data handbook where the experience and skill of the operator are required. This paper presents an optimal framework to control the cutting force of the drilling process. A mathematical model that captures complex drilling dynamics between cutting force and feed rate based on system identification is used. Then, a Proportional-Integral-Derivative (PID) controller is proposed to control the cutting force. Taking advantage of up-to-date swam-based optimization technique, an Enhance Bat Algorithm (EBA) approach is used to tune the design variables of the PID controller based on the Integral Absolute Error (IAE) criterion. The results are compared with another two swam optimization, the Particle Swarm Optimization (PSO) and the Whale Optimization Algorithm (WOA). The comparison reveals that EBA can give better results in terms of improving time domain specifications and reducing error performance indices.

## 1. INTRODUCTION

Drilling operation is one of the common metal-removal processes. It is widely used in manufacturing and has a direct impact on the quality of production [1, 2]. Holes are necessary for the assembly of different items to be combined and finalize the product [3]. Therefore, improving the drilling process leads to significant economic benefits such as improving productivity and reducing costs. The friction and the temperature between the drilling operation and workpiece increase as the depth of the drilling increase. This leads to tool vibration and tool wear [4, 5]. Moreover, sometimes workpiece has a multi-layered (i.e. composite structure) that causes changes in the cutting surface stiffness [6]. Delamination of composite materials during the drilling process often happened due to insufficiently controlling the cutting force [7]. Feed rate of drilling operation is the important variable that needs to be controlled based on the cutting condition to improve the drill operation. Usually, the operator changes the feed rate during the drilling process based on his experience to work on a safe operation condition [8].

In order to control the cutting force, it is important to have an accurate mathematical model of the drilling process. Researchers have been developed a different data-driven mathematical model of the drilling process with different operating conditions. For example, Roukema and Altintas [9] proposed a mathematical model of the drilling process considering cutting speed, type of tool used and workpiece geometries as inputs and cutting forces and vibrations as outputs. The model has been formulated from a linear least squares regression. del Toro et al. [10] established a

mathematical model between resultant force and feed rate in the drilling process using conventional identification methods. In the same way, Singh and Sharma [7] developed a mathematical model of a drilling process of a composite material. A transfer function formulation and state-space representation of the process was obtained by converting between thrust force and feed rate. On another side, different controller strategies have been investigated to control the drilling process. The majority of the papers utilized PID controller as the controller strategy [11-16].

Despite the existing extensive research on the tuning the PID controller, there is still a significant gap in understanding the best tuning strategies. In practice, finding optimal tuning for PID controllers is a challenging problem [16]. Nowadays, there is an increasing trend of using optimization in solving industrial problems and engineering research [17-21]. Taking advantage of recent optimization techniques, several tuning methods based on swam based optimization has been proposed for finding the best setting of numerous controllers [22, 23]. In the context of finding the optimal tuning of the PID controller in drilling process, different swarm optimization techniques have been proposed such as Simulated Annealing (SA) [11], PSO [13], WOA [15] and Enhanced Flower Pollination (EFP) [16].

In this paper, an optimal framework to control the cutting force of the drilling process is present. A mathematical model that captures complex drilling dynamics between thrust force and feed rate based on system identification is used. Then, a PID controller is proposed to control the cutting force. An EBA approach is used to tune the design variables of the PID controller based on IAE criterion. The rest of the paper is
organized as follows: Section 2 presents the mathematical model of the drilling process. In Section 3 a brief explanation of the PID controller is given. The proposed optimization algorithm is explained in Section 4. Section 5 shows the simulation and comparative study. Section 6 summarizes the conclusion.

## 2. DRILLING PROCESS MODELLING

Mathematical model of a system provides a way to represent the behavior of that system using mathematical equations. By developing a model of a system, it can be used to gain deep insights into how the system works, predict its behavior under different work conditions, and design control strategies to achieve desired outcomes. Drilling is a machining process that involves complex interactions between the work piece and the cutting tool [12]. The mathematical model that is considered in this work of the drilling process is obtained experimentally using system identification. The specifications of the considered drilling process during the experiments are reported in Table 1 as given in literature [12].

<table><caption>**Table 1.** Drilling specification [12]</caption><thead><tr><td>Parameter</td><td>Specification</td></tr></thead><tbody><tr><td>Power</td><td>26 KW</td></tr><tr><td>Nominal feed rate</td><td>100 mm/min</td></tr><tr><td>Nominal spindle speed</td><td>870 rpm</td></tr><tr><td>Tool diameter</td><td>10 mm</td></tr><tr><td>Maximum hole depth</td><td>30 mm</td></tr><tr><td>Drilling tool coating</td><td>TiN/TiAlN</td></tr></tbody></table>

The modeling of the drilling process consists of modeling the feed drive, the spindle and the cutting process [13]. The input/output data from all the experiments are then processed. An open loop single-input single-output (SISO) mathematical model that best approximates the actual drilling process was obtained as a 3rd order transfer function. The feed rate is used as an input and the cutting force is the output. The 3rd order transfer function between cutting force and feed rate is given by [11, 12]:

$$G(s) = \frac{F(s)}{f(s)} = \frac{1958}{s^3 + 17.89s + 103.3s + 190.8} \quad (1)$$

where,

s: Laplace operator

G: Drilling-process transfer function

F: Cutting force

f: Feed rate

## 3. PID CONTROLLER

PID controller is used for stabilization and disturbance rejection of dynamical systems. It provides effective solutions to a majority of the control engineering problems [24]. PID consists of three terms of the proportional (P), integral (I) and derivative (D) as shown in Figure 1. The proportional term changes the control signal proportionally to the error. The proportional term makes the system faster and reduces the error steady-state. The integral term changes the control signal proportionally to the integration of the error. The integral term eliminates the error steady-state. The derivative term changes

the control signal proportionally to the derivative of the error. This term reduces the settling time of the system [25].

<figure><img src="image_2.png" /><figcaption>**Figure 1.** The block diagram of PID controller</figcaption></figure>

The overall action of the controller is the sum of the contributions from these three terms. The final control signal $u(t)$ of the PID controller is the sum of the three terms as given [26, 27]:

$$u(t) = K_p e(t) + K_i \int_{0}^{t} e(t) dt + K_d \frac{de}{dt} \quad (2)$$

The overall transfer function of the PID controller can be represented as follows:

$$G_{PID}(s) = \frac{U(s)}{E(s)} = K_p + \frac{K_i}{s} + K_d s \quad (3)$$

where,

G<sub>PID</sub>: PID transfer Function

U: Control signal

E: Error

K<sub>p</sub>: Proportional gain

K<sub>i</sub>: Integral gain

K<sub>d</sub>: Derivative gain

## 4. ENHANCE BAT ALGORITHM

Bat Algorithm (BA) is a population-based swarm optimization algorithm. It was developed by Yang [28]. BA mimics the echolocation behavior of bats. Echolocation is a strategy used by bats to estimate distance and distinguish between prey and obstacles [29]. By using echolocation, the bat sends a series of loud ultra-sound waves and then receives back echoes with various sound levels. This process enables bats to identify between prey and obstacles [30]. Yang [28] assumed three approximations as follows:

The echolocation is used by all bats to discover the distance, and they know the difference between prey and obstacles when an echo is received back with various sound levels [31].

All bats have a velocity v, position x and frequency f. For searching for prey, bats fly randomly with fixed frequency and varying wavelength λ and loudness A [32]. Depending how they far from the target, they repeatedly adjust the wavelength of their emitted pulses and adjust the rate of pulse emission R ∈ (1,0].

Loudness A reduces from the maximum value A₀ to its minimum value A_min.

In the bat algorithm, the frequency fᵢ, velocity vᵢ and position xᵢ are updated for each bat in as follows [28]:

$$f_i = f_{max} \times r_1 s \quad (4)$$
$$v_i(itr + 1) = v_i(itr) + (x_i(itr) - x_g) \times f_i \quad (5)$$

$$x_i(itr + 1) = x_i(itr) + v_i(itr + 1) \quad (6)$$

where,

* i: Counter for population ($i=1,2,3,...,N$)

* itr: Counter for iteration ($itr=1,2,3,...,T$)

* fᵢ: Frequency for bat i

* fmax: Maximum frequency

* r₁: Random value between [0,1]

* vᵢ: Velocity for bat i

* xᵢ: Position for bat i

* xg: Position of the bat that has the best fitness

The next search in bat algorithm is based on two parameters: loudness A and the pulse emission rate R as follows: Select a solution among best solutions. Then, select a round number (Rand) if Rand < Rᵢ, generate new local solution around the best solution as given in Eq. (7).

$$x_i(itr + 1) = x_i(itr) + r_2 A_p(itr) \quad (7)$$

where,

* r₂: Random value between [0,1]

* Aₚ(itr): Average loudness value of all bats in the current iteration

Furthermore, a new solution randomly is generated based on Eq. (8).

$$x_i(itr + 1) = x_i(itr) + r_2 A_p(itr) \quad (8)$$

If Rand > Aᵢ and the new solution fₓ₁ have better objective value than the best one found by the algorithm fₓg, the algorithm will accept the new solution and loudness and pulse emission rate will be updated based on Eqs. (9) and (10).

$$A_i(itr + 1) = \alpha A_i(itr) \quad (9)$$

$$R_i(itr + 1) = R_i(1 - e^{-\gammaitr}) \quad (10)$$

where, α and γ are constant and for simplicity could be considered equal.

It can be noticed that BA compromise between PSO and SA. However, in this work the BA enhanced by improving the update processes of velocity. The improvement is made by adding a new term to the Eq. (5) to let the bat guided by its best location found xᵢₚ in addition to the global guidance xg. This is similar to the PSO. The Eq. (5) after the modification becomes [33]:

$$v_i(itr + 1) = v_i(itr) + (x_i(itr) - x_{ip}) \times f_i + (x_i(itr) - x_g) \times f_i \quad (11)$$

The pseudo code of EBA is illustrated in Algorithm 1.

### Algorithm 1. Pseudo code of EBA

#### 1. Input

* Objective function, Population size (N), Maximum frequency ($f_{max}$), Coefficient value (α), Number of iteration (T)

#### 2. Initialization

* Initialize velocity v, position x, loudness A and the pulse emission rate R of the population

* Evaluate objective function and assign xₚᵢ and xg

#### 3. Loop:

* while (itr < T)

* For each bat in the population (N)

* Update population based on Eq. (4), Eq. (11) and Eq. (6)

* End for

* Select a random number (Rand)

* If rand < Rᵢ

Select a solution randomly

Generate new local solution based on Eq. (7)

* If rand > Aᵢ

Generate new solution randomly based on Eq. (8)

Update loudness based on Eq. (9)

Update pulse emission rate based on Eq. (10)

* Update xₚᵢ

* End for

* Update xg

* $itr =itr + 1$

* End while

#### 4. Print the Optimal Solution

### 5. SIMULATION STUDY

In this section, the simulation results of controlling the cutting force of the drilling process that is given by the transfer function in Eq. (1) using the PID controller as given in Eq. (3) are reported. EBA is employed to find the optimal value of the gain parameters of the PID controller. IAE index as given in Eq. (12) [34] is used as a cost function of the optimization.

$$IAE = \int_{t=0}^{t=T} |e(t)| dt \quad (12)$$

The parameters of the EBA are listed in Table 2.

<table><caption>**Table 2.** EBA algorithm parameters</caption><thead><tr><td>Parameters</td><td>Value</td></tr></thead><tbody><tr><td>Population size (N)</td><td>50</td></tr><tr><td>Number of iterations (T)</td><td>100</td></tr><tr><td>Maximum Frequency (f<sub>max</sub>)</td><td>0.8</td></tr><tr><td>Loudness</td><td>0.9</td></tr><tr><td>Pulse emission rate</td><td>0.5</td></tr></tbody></table>

A step change in the input force is considered and the goal of the PID controller is to follow set-point changes. The simulation was executed in MATLAB and the values of the parameters of the PID controller are found as given:

$$K_p = 0.8828, K_i = 1.8188, K_d = 0.2227$$

Figure 2 shows the response of the system to unit step input. Comparatives analysis was performed with another two swam optimizations, PSO and WOA. The value of the designed gains Kₚ, Kᵢ and K₉ of the PID controller that were obtained by the proposed EAB, PSO and WOA is given in Table 3.

For evaluation the proposed EBA, time domain specifications comparison as given in Table 3 is made with the results obtain by PSO and WOA.

It can be noticed from Figure 2 and also from the results that reported in Table 4 that the proposed EBA-PID controller is a compromise solution between the PSO-PID and WOA-PID. For instance, EBA-PID has less overshoot than WOA-PID but
need longer time to reach steady state. In the same way, EBA-PID has less settling time than PSO-PID but it has more overshoot.

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Unit step response of the PSO, WAO and EBA based PID controller</figcaption></figure>

<table><caption>**Table 3.** PID controller values using PSO, WAO and EBA based PID controller</caption><thead><tr><td>Tuning Algorithm</td><td>k<sub>p</sub></td><td>k<sub>i</sub></td><td>k<sub>d</sub></td></tr></thead><tbody><tr><td>PSO-PID [13]</td><td>0.7164</td><td>1.4665</td><td>0.0984</td></tr><tr><td>WOA-PID [15]</td><td>1.13</td><td>2</td><td>0.31488</td></tr><tr><td>EBA-PID (Proposed)</td><td>0.8828</td><td>1.8188</td><td>0.2227</td></tr></tbody></table>

<table><caption>**Table 4.** Time domain specifications of the PSO, WAO and EBA based PID controller</caption><thead><tr><td>Time Domain Specifications</td><td>PSO-PID</td><td>WOA-PID</td><td>EBA-PID</td></tr></thead><tbody><tr><td>Overshoot (%)</td><td>18</td><td>31</td><td>24</td></tr><tr><td>Peak time (s)</td><td>0.2</td><td>0.13</td><td>0.15</td></tr><tr><td>Rising time (s)</td><td>0.17</td><td>0.054</td><td>0.067</td></tr><tr><td>Settling time (s)</td><td>1.4</td><td>0.588</td><td>0.64</td></tr></tbody></table>

In addition, a number of error performance indices namely IAE as given in Eq. (12) [34], Integral Time of Absolute Error (ITAE) as given in Eq. (13) [35], Integral Square of Error (ISE) as given in Eq. (14) [36] has been employed to compare the performance of the proposed EBA algorithm as revealed in Table 5.

$$ITAE = \int_{t=0}^{t=t_{sim}} t|e(t)|dt \qquad (13)$$

$$ISE = \int_{t=0}^{t=t_{sim}} e^2(t)dt \qquad (14)$$

<table><caption>**Table 5.** Error performance indices of the PSO, WAO and EBA based PID controller</caption><thead><tr><td>Performance Indices</td><td>PSO-PID</td><td>WOA-PID</td><td>EBA-PID</td></tr></thead><tbody><tr><td>IAE</td><td>1.9047</td><td>1.5195</td><td>1.3933</td></tr><tr><td>ITAE</td><td>0.2157</td><td>0.1954</td><td>0.1848</td></tr><tr><td>ISE</td><td>1.2291</td><td>1.0626</td><td>1.0406</td></tr></tbody></table>

It can be clearly seen from Table 5 that the PID controller

designed by EBA reduces error performance indices in comparison with the PID controller designed by PSO and WOA. For example, EBA-PID reduces the IAE of the system by 26.85% and 8.3% in comparison with PSO-PID and WOA-PID respectively. Besides, the ITAE of the system with EBA-PID is reduced by 14.3% and 5.4% in comparison with PSO-PID and WOA-PID respectively. In the same way, the ISE is reduced by 15.3% and 2% in comparison with PSO-PID and WOA-PID respectively. These results show superiority of EBA over other two methods to tuning the designed gains $k_p$, $k_i$ and $k_d$ of the PID controller.

The implication of the current study is to improve the control algorithm of the drilling process using the compensation of the PID controller with the swarm optimization techniques. The economic advantages from this improvement of the closed-loop drilling control by keeping the value of the feed rate close to the set point without overshoot prevents damage to the workpiece surface and drilling tool.

## 6. CONCLUSIONS

Drilling process is the most commonly manufacturing operation in today's industries. Feed rate of drilling process is the important variable that needs to be controlled based on the cutting condition to improve the drill operation. Due to the simplicity and the low cost compared to other controllers, PID controller has been implemented in this paper to control the cutting force during the drilling process. To find the optimal value of the designed gains $k_p$, $k_i$ and $k_d$ of the PID controller, an Enhance Bat Algorithm (EBA) approach is used based on the Integral Absolute Error (IAE) criterion. A step change in the input force is considered and the goal of the PID controller is to follow set-point changes. To evaluate the proposed tuning method, the results is compared with two swarm optimization, the Particle Swarm Optimization (PSO) and the Whale Optimization Algorithm (WOA). In terms of time domain specification, it is revealed that there is an improvement such as reducing settling time and reduce overshoot. In terms of the performance induces that are based on error, the proposed tuning algorithm achieved better results than the other swarm optimization that are considered in this paper. The results show that the EBA-PID reduces the IAE of the system by 26.85% and 8.3% in comparison with PSO-PID and WOA-PID respectively. Moreover, the ITAE of the system with EBA-PID is reduced by 14.3% and 5.4% in comparison with PSO-PID and WOA-PID respectively. Besides, the ISE is reduced by 15.3% and 2% in comparison with PSO-PID and WOA-PID respectively.

Future research for further improving in the system performance can be performed in different directions. One possible improvement is by investigation a new novel optimization technique to tune the design parameters of the PID controller. Other direction for future work, the proposed tuning method for the PID can be used in other applications beyond the context of drilling process.

## REFERENCES

[1] Furness, R.J., Tsao, T.C., Rankin, J.S., Muth, M.J. Manes, K.W. (1999). Torque control for a form tool drilling operation. IEEE Transactions on Control Systems Technology, 7(1): 22-30.
[2] Sharif, A., Hussain, A., Habib, N., Alam, W., Hanif, M.I., Noon, A.A., Khan, M.I. (2023). Experimental investigation of hole quality and chip analysis during the dry drilling process of Al6061-T6. Journal of Materials and Manufacturing, 2(1): 21-30.
https://doi.org/10.5281/zenodo.8020538

[3] Kyratsis, P., Markopoulos, A.P., Efkolidis, N., Maliagkas, V., Kakoulis, K. (2018). Prediction of thrust force and cutting torque in drilling based on the response surface methodology. Machines, 6(2): 24.
https://doi.org/10.3390/machines6020024

[4] Oh, Y.T., Kim, G.D., Chu, C.N. (2003). Design of a drilling torque controller for a machining center. The International Journal of Advanced Manufacturing Technology, 22: 329-335.
https://doi.org/10.1007/s00170-002-1503-z

[5] Roukema, J.C., Altintas, Y. (2006). Time domain simulation of torsional-axial vibrations in drilling. International Journal of Machine Tools and Manufacture, 46(15): 2073-2085.
https://doi.org/10.1016/j.ijmachtools.2005.12.010

[6] Seif, A., Fathy, A., Megahed, A.A. (2023). Effect of drilling process parameters on bearing strength of glass fiber/aluminum mesh reinforced epoxy composites. Scientific Reports, 13(1): 12143.
https://doi.org/10.1038/s41598-023-39097-3

[7] Singh, A.P., Sharma, M. (2013). Modelling of thrust force during drilling of fibre reinforced plastic composites. Procedia Engineering, 51: 630-636.
https://doi.org/10.1016/j.proeng.2013.01.089

[8] Kim, J.B., Lee, S.J., Park, Y.P. (1994). Stable and efficient drilling process by active control of the thrust force. Mechanical Systems and Signal Processing, 8(5): 585-595. https://doi.org/10.1006/mssp.1994.1041

[9] Roukema, J.C., Altintas, Y. (2004). Kinematic model of dynamic drilling process. ASME International Mechanical Engineering Congress and Exposition, 47136: 955-963. https://doi.org/10.1115/IMECE2004-59340

[10] del Toro, R.L.M., Schmittdiel, M.C., Haber-Guerra, R.E., Haber-Haber, R. (2007). System identification of the high performance drilling process for network-based control. In International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Las Vegas, Nevada, USA, pp. 827-834. https://doi.org/10.1115/DETC2007-34307

[11] Haber, R.E., Haber-Haber, R., Del Toro, R.M., Alique, J.R. (2007). Using simulated annealing for optimal tuning of a PID controller for time-delay systems. An application to a high-performance drilling process. In Computational and Ambient Intelligence: 9th International Work-Conference on Artificial Neural Networks, IWANN 2007, San Sebastián, Spain, pp. 1155-1162. https://doi.org/10.1007/978-3-540-73007-1_140

[12] Haber-Haber, R., Haber, R., Schmittdiel, M., del Toro, R.M. (2007). A classic solution for the control of a high-performance drilling process. International Journal of Machine Tools and Manufacture, 47(15): 2290-2297.
https://doi.org/10.1016/j.ijmachtools.2007.06.007

[13] Giriraj Kumar, S.M., Jayaraj, D., Kishan, A.R. (2010). PSO based tuning of a PID controller for a high performance drilling machine. International Journal of

[14] Singh, A.P., Sharma, M. (2014). Modeling and PID control of thrust force during drilling in composite laminates. In 2014 Recent Advances in Engineering and Computational Sciences (RAECS), Chandigarh, India, pp. 1-5. https://doi.org/10.1109/RAECS.2014.6799513

[15] Kumar, A.A., Kumar, S.G. (2018). Application of whale optimization algorithm for tuning of a PID controller for a drilling machine. In ICAARS 2018, Coimbatore, India, No. fhlal-02314429.

[16] Mahmood, Z.N., Al-Khazraji, H., Mahdi, S.M. (2023). PID-based enhanced flower pollination algorithm controller for drilling process in a composite material. Annales de Chimie Science des Matériaux, 47(2): 91-96.
https://doi.org/10.18280/acsm.470205

[17] Al-Khazraji, H., Nasser, A.R., Khlil, S. (2022). An intelligent demand forecasting model using a hybrid of metaheuristic optimization and deep learning algorithm for predicting concrete block production. IAES International Journal of Artificial Intelligence, 11(2): 649. http://doi.org/10.11591/ijai.v11.i2.pp649-657

[18] Al-Khazraji, H. (2022). Comparative study of whale optimization algorithm and flower pollination algorithm to solve workers assignment problem. International Journal of Production Management and Engineering, 10(1): 91-98. https://doi.org/10.4995/ijpme.2022.16736

[19] Al-Khazraji, H., Khlil, S. and Alabacy, Z. (2021). Cuckoo search optimization for solving product mix problem. IOP Conference Series: Materials Science and Engineering, 1105(1): 012016.
https://doi.org/10.1088/1757-899X/1105/1/012016

[20] Al-Khazraji, H., Khlil, S., Alabacy, Z. (2020). Industrial picking and packing problem: Logistic management for products expedition. Journal of Mechanical Engineering Research and Developments, 43(2): 74-80.

[21] Al-Khazraji, H., Cole, C., Guo, W. (2018). Multi-objective particle swarm optimisation approach for production-inventory control systems. Journal of Modelling in Management, 13(4): 1037-1056. https://doi.org/10.1108/JM2-02-2018-0027

[22] Al-Khazraji, H., Naji, R.M., Khashan, M.K. (2024). Optimization of sliding mode and back-stepping controllers for AMB systems using gorilla troops algorithm. Journal Européen des Systèmes Automatisés, 57(2): 417-424. https://doi.org/10.18280/jesa.570211

[23] Ahmed, A.K., Al-Khazraji, H., Raafat, S.M. (2024). Optimized PI-PD control for varying time delay systems based on modified smith predictor. International Journal of Intelligent Engineering & Systems, 17(1): 331-342. https://doi.org/10.22266/ijies2024.0229.30

[24] Nagaraj, B., Murugananth, N. (2010). A comparative study of PID controller tuning using GA, EP, PSO and ACO. In 2010 International Conference on Communication Control and Computing Technologies, Nagercoil, India, pp. 305-313. https://doi.org/10.1109/ICCCCT.2010.5670571

[25] Ahamed, S.R., Parumasivam, P., Lipu, M.H., Hannan, M.A., Ker, P.J. (2020). A comparative evaluation of PID-based optimisation controller algorithms for DC motor. International Journal of Automation and Control, 14(5-6): 634-655. https://doi.org/10.1504/IJAAC.2020.110076

[26] Blondin, M.J., Sanchis Sáez, J., Pardalos, P.M. (2019). Control engineering from classical to intelligent control
theory - An overview. Computational Intelligence and Optimization Methods for Control Engineering, pp. 1-30. https://doi.org/10.1007/978-3-030-25446-9_1

[27] Gopi, P., Reddy, S.V., Bajaj, M., Zaitsev, I., Prokop, L. (2024). Performance and robustness analysis of V-Tiger PID controller for automatic voltage regulator. Scientific Reports, 14(1): 7867. https://doi.org/10.1038/s41598-024-58481-1

[28] Yang, X.S. (2010). A new metaheuristic bat-inspired algorithm. In Nature Inspired Cooperative Strategies for Optimization (NICSO 2010), pp. 65-74. https://doi.org/10.1007/978-3-642-12538-6_6

[29] Singh, K., Vasant, P., Elamvazuthi, I., Kannan, R. (2015). PID tuning of servo motor using bat algorithm. Procedia Computer Science, 60: 1798-1808. https://doi.org/10.1016/j.procs.2015.08.290

[30] Katal, N., Kumar, P., Narayan, S. (2014). Optimal PID controller for coupled-tank liquid-level control system using bat algorithm. In 2014 International Conference on Power, Control and Embedded Systems (ICPCES), Allahabad, India, pp. 1-4. https://doi.org/10.1109/ICPCES.2014.7062818

[31] Yuan, X., Yuan, X., Wang, X. (2021). Path planning for mobile robot based on improved bat algorithm. Sensors, 21(13): 4389. https://doi.org/10.3390/s21134389

[32] Alharbi, A., Alosaimi, W., Alyami, H., Rauf, H.T., Damaševičius, R. (2021). Botnet attack detection using local global best bat algorithm for industrial internet of things. Electronics, 10(11): 1341. https://doi.org/10.3390/electronics10111341

[33] Yılmaz, S., Küçüksille, E.U. (2015). A new modification approach on bat algorithm for solving optimization problems. Applied Soft Computing, 28: 259-275. https://doi.org/10.1016/j.asoc.2014.11.029

[34] AL-Khazraji, H., Cole, C., Guo, W. (2021). Optimization and simulation of dynamic performance of production-inventory systems with multivariable controls. Mathematics, 9(5): 568. https://doi.org/10.3390/math9050568

[35] Ahmed, A.K., Al-Khazraji, H. (2023). Optimal control design for propeller pendulum systems using gorilla troops optimization. Journal Européen des Systèmes Automatisés, 56(4): 575-582. https://doi.org/10.18280/jesa.560407

[36] Sule, A.H., Mokhtar, A.S., Jamian, J.J.B., Khidrani, A., Larik, R.M. (2020). Optimal tuning of proportional integral controller for fixed-speed wind turbine using grey wolf optimizer. International Journal of Electrical

and Computer Engineering, 10(5): 5251-5261. https://doi.org/10.11591/ijece.v10i5.pp5251-5261

## NOMENCLATURE

$A_p(itr)$ Average loudness value of all bats in the current iteration

$E$ Error

$F$ Cutting force

$f$ Feed rate

$f_i$ Frequency for bat i

$f_{max}$ Maximum frequency

$G(s)$ Drilling-process transfer function

$G_{PID}(s)$ PID transfer Function

$i$ Counter for population (i = 1,2,3,...N)

$itr$ Counter for iteration (itr = 1,2,3,...T)

$K_d$ Derivative gain

$K_i$ Integral gain

$K_p$ Proportional gain

$r_1$ Random value between [0,1]

$r_2$ Random value between [0,1]

$s$ Laplace operator

$U$ Control signal

$v_i$ Velocity for bat i

$x_i$ Position for bat i

$x_g$ Position of the bat that has the best fitness

## Greek symbols

$\alpha$ Coefficient in EBA

$\gamma$ Coefficient in EBA

## Acronym

BA Bat Algorithm

EBA Enhance Bat Algorithm

IAE Integral Absolute Error

ITAE Integral Time of Absolute Error

ISE Integral Square of Error

PSO Particle Swarm Optimization

PID Proportional-Integral-Derivative

SA Simulated Annealing

SISO Single-Input Single-Output

WOA Whale Optimization Algorithm