Prediction of Thrust Force and Torque for High-Speed Drilling of AL6061 with TMPTO-Based Bio-Lubricants Using Machine Learning

: This study was designed to examine the effects of a trimethylolpropane trioleate (TMPTO)- based lubricant on thrust force and torque under the high-speed drilling of Al-6061 as an effective environmentally friendly cutting ﬂuid. The tribological performance of three lubricant blends was evaluated based on ASTM standards. TMPTO base oil, notably enhances load-carrying capacity under extreme pressure conditions, with a seizer load of 7848 N. The best-performing oil was further optimized using a Taguchi-based design experiment to investigate the effect of different additive concentrations on thrust force and torque under actual contact conditions. Experiments were conducted using three critical machining parameters: additive concentration, spindle speed, and feed rate. The results of the ANOVA analysis reveal that spindle speed contributes most substantially (62.99%) to torque, with feed rate (23.72%) and additive concentration (7.74%) also showing signiﬁcant impacts. On the other hand, thrust force is primarily inﬂuenced by feed rate (73.52%), followed by spindle speed (16.82%), and additive concentration (6.28%). Furthermore, a machine learning model was developed to predict and compare a few signiﬁcant aspects of high-speed drilling machinability, including thrust force and torque. Three different error metrics were utilized in order to assess the performance of the predicted values, namely the coefﬁcient of determination (R2), mean absolute percentage error (MAPE) and mean square error (MSE), which are all based on the coefﬁcient of determination. Compared to other models, decision tree produces more accurate prediction values for cutting forces. The present study provides a novel approach for evaluating the most promising biodegradable lube oils and predicting cutting forces by formulating a perfect blend.


Introduction
A collaborative effort between researchers and industry experts is actively involved in improving the characteristics of vegetable-based lubricants to make them sustainable cutting fluids. Additive technology and the chemical modification of vegetable oils are some methods being explored to enhance their tribological properties and overall performance. With ongoing research and development, vegetable-based lubricants will likely improve performance and eventually replace mineral-based cutting fluids in various applications [1,2]. Cutting fluid dissipates the heat generated in the metal cutting process and In the context of drilling fluid applications, trimethylolpropane ester (TMPE) was investigated as an additive, specifically as a nonionic surfactant. It was studied alongside other potential polyol esters, including neopentyl glycol ester (NPGE) and pentaerythritol ester. The purpose of this examination was likely to assess the suitability and effectiveness of TMPE, NPGE, and pentaerythritol ester as additives in drilling fluids, potentially evaluating their impact on various properties such as viscosity, lubricity, stability, or other relevant performance factors [7,20,21].
Machine learning (ML) has found numerous applications in machining processes across various contexts, such as roughness quality prediction, tool condition monitoring, process optimization, quality control, etc. [21,22]. The versatility and adaptability of machine learning make it a powerful tool for optimizing and enhancing various aspects of machining operations and cutting tools [23,24]. Various techniques are used for the prediction of cutting forces and torque in drilling, such as artificial neural networks [25,26], the response surface methodology [27,28], and the use of different machine learning algorithms [29][30][31] under dry as well as lubricated conditions.
The predictions developed using advanced machine learning techniques are more accurate in optimizing machining operations and reducing production costs. It is worth noting that publicly accessible research on machine learning models specifically developed for machining with cutting fluids may be relatively limited compared to that for dry machining. Moreover, there is very limited information available regarding the prediction of thrust force and torque under biodegradable lubrication conditions in order to be able to make accurate predictions. Cutting fluids can affect temperature distribution, chip disposal, friction, and tool life, influencing the forces involved in cutting operations and torque. Including cutting fluid-related parameters in the model can help capture these effects and improve the accuracy of the predictions.
Investigating the impact of TMPTO-based lubricant on thrust force and torque during the high-speed drilling of Al-6061 is essential to enhancing machining performance and prolong tool life, improving surface quality, boosting energy efficiency, and ensuring overall process sustainability. The objective of this study was to investigate the drilling process of Al 6061 using a TMPTO-based bio-lubricant as the drilling fluid. The main focus was on understanding and predicting the cutting forces and torque involved in high-speed drilling operations. In addition, the current study presents four regression-based machine learning techniques (linear regression, decision Trees, random forest and support vector machine) to construct predictive models that determine how machining forces develop during turning. An optimization problem was developed to find the exact combination of parameters so that maximum accuracy could be achieved with the fewest resources possible. This would ensure that the model was efficient and effective in its predictions.

Lubricants
TMPTO (trimethylolpropane trioleate) was chosen as the base oil due to its significant rheological properties and long-chain fatty acid content [6,7,32,33]. Three additives were obtained from the local market to formulate the blend. The blends were thoroughly mixed using a magnetic stirrer to ensure uniformity. Among the three formulated samples, sample A contains chemically modified sulfurized vegetable oil additive, which has a sulfur content of 15%, while sample B contains a sulfurized vegetable oil-based olefin, which has a sulfur content of 26.6%, and sample C has a sulfur content of 39.7%. The rheological properties of additives are presented in Table 1 as influenced by a lube oil formulator.
Infrared spectroscopy was used to characterize the molecular structure of the formulated TMPTO base oil, as depicted in Figure 1. Olefinic (alkane) C-H stretching and bending vibration peaks were observed at 3004.38 cm −1 and 1461.37 cm −1 to 1354.68 cm −1 , respectively. Peaks for the carbonyl (C=O) stretching vibration and methylene (C-H) bending vibration were observed at 1740.90 cm −1 and 1461.37 cm −1 , respectively. Trimethylolpropane trioleate (TMPTO) peaked at 721.33 cm −1 due to its long alkyl chain. Infrared spectroscopy was used to characterize the molecular structure of the formulated TMPTO base oil, as depicted in Figure 1. Olefinic (alkane) C-H stretching and bending vibration peaks were observed at 3004.38 cm −1 and 1461.37 cm −1 to 1354.68 cm −1 , respectively. Peaks for the carbonyl (C=O) stretching vibration and methylene (C-H) bending vibration were observed at 1740.90 cm −1 and 1461.37 cm −1 , respectively. Trimethylolpropane trioleate (TMPTO) peaked at 721.33 cm −1 due to its long alkyl chain.

Four-Ball Test
Experiments were carried out on the DUCOMTR-30L four-ball tester to investigate the effect of additives in TMPTO base oil on wear. The test setup consists of four steel balls (made up of AISI E52100 chrome alloy) with 62-64 HRC arranged in a triangular configuration. Three balls are fixed at the base, forming a lower stationary ball pot, while the fourth ball, known as the top ball, is placed on top and loaded against the stationary balls. The ball holder assembly securely holds the four balls in place during the test. The holder usually has the ability to rotate the top ball at a specific speed. The ball holder assembly allows for the lubricant under study to be applied to the balls during the experiment. It also houses the temperature and load control mechanisms. An arrangement is in place to apply a specific load on the top ball to create contact pressure between the balls. This load is critical in determining the wear and extreme pressure properties of the lubricant. The temperature of the lubricant in the ball pot, the rotation speed of the top ball, and the applied load were all continuously monitored using an electronic controller.

Four-Ball Test
Experiments were carried out on the DUCOMTR-30L four-ball tester to investigate the effect of additives in TMPTO base oil on wear. The test setup consists of four steel balls (made up of AISI E52100 chrome alloy) with 62-64 HRC arranged in a triangular configuration. Three balls are fixed at the base, forming a lower stationary ball pot, while the fourth ball, known as the top ball, is placed on top and loaded against the stationary balls. The ball holder assembly securely holds the four balls in place during the test. The holder usually has the ability to rotate the top ball at a specific speed. The ball holder assembly allows for the lubricant under study to be applied to the balls during the experiment. It also houses the temperature and load control mechanisms. An arrangement is in place to apply a specific load on the top ball to create contact pressure between the balls. This load is critical in determining the wear and extreme pressure properties of the lubricant. The temperature of the lubricant in the ball pot, the rotation speed of the top ball, and the applied load were all continuously monitored using an electronic controller.

Antiwear and Extreme Pressure Tests
Ducom's four-ball tester is employed for conducting AW/EP (Anti-Wear/Extreme Pressure) tests per ASTM D4172-94 and ASTM D2783-88 standards [34,35]. The lubricant evaluation process involves assembling three 12.7 mm diameter steel balls tightly together and covering them with the lubricant under examination. A 12.7 mm diameter steel ball is pressed with a force of either of 392 N into the space formed by the three clamped balls, creating a three-point contact configuration. The temperature of the test lubricant is precisely controlled at 75 • C, and then the top ball is set in motion, rotating at 1200 rpm) for a duration of 60 min. The lubricants are compared based on the average size of the wear scars formed on the three lower clamped balls [34].
After the test, a solvent organic compound such as acetone removes the oil from the wear scars on the balls. Subsequently, an optical microscope is utilized to measure the diameter of the wear scars.
EP tests determine the lubricant's ability to prevent wear and indicate the load at which it fails. ASTM D2783 is used to measure EP properties. The parameters used to measure the device's properties are time = 10 s, temperature = 27 • C, RPM = 1760, and load = 6-800 kgf, as defined in the standard [35]. The test procedure involves conducting a series of tests, each lasting 10 s. The tests are performed at increasing loads until welding occurs. Welding refers to the phenomenon where the load-carrying ability of the lubricant is exceeded, resulting in the welding of the steel balls due to extreme pressure. The load at which welding takes place provides a measure of the extreme pressure properties of the lubricant. Higher loads before welding indicate better EP characteristics of the lubricating fluid. The results and discussion section of the study includes a comprehensive analysis and discussion of the AW (anti-wear) and EP (extreme pressure) test outcomes.

Drilling Condition
We intended in this work to conduct a comparative analysis of a variety of different approaches for modeling the drilling process of Al 6061 in order to compare the results. Aluminum 6061 (Al-6061) is a widely used alloy in various industries due to its excellent combination of mechanical properties, good formability, and high corrosion resistance. It belongs to the 6xxx series of aluminum alloys and is commonly used for structural applications, such as in aerospace, automotive, marine, and general engineering. The tensile strength of Aluminum 6061 typically ranges from about 140 to 290 MPa. Experiments on an Al 6061 workpiece were conducted to obtain the datasets for the models. A 150 mm × 150 mm × 20 mm plate was used for this investigation, the chemical composition and major mechanical properties of which are depicted in Tables 2 and 3, respectively. Drilling experiments were performed using HSS drill tools on a BFW Chakra BMV 60 CNC machining center. Thrust force and torque measurements were captured using a wireless Spike tool holder dynamometer. Experimental cutting force values were monitored and recorded using a data acquisition system, as depicted in Figure 2.
an Al 6061 workpiece were conducted to obtain the datasets for the models. A 150 mm × 150 mm × 20 mm plate was used for this investigation, the chemical composition and major mechanical properties of which are depicted in Tables 2 and 3, respectively. Drilling experiments were performed using HSS drill tools on a BFW Chakra BMV 60 CNC machining center. Thrust force and torque measurements were captured using a wireless Spike tool holder dynamometer. Experimental cutting force values were monitored and recorded using a data acquisition system, as depicted in Figure 2.

Experimental Design
This study aimed to investigate the influence of additive concentration, spindle speed, and feed rate on thrust force and torque. The chosen parameters directly impact the lubrication conditions between contacting surfaces, affecting film thickness, friction, and tool wear. Understanding the relationship between these parameters and thrust force and torque is essential to optimizing lubrication performance and enhancing the efficiency and reliability of the tribological system under actual contact conditions. A Taguchi L9 orthogonal array that systematically varies these parameters while measuring thrust force and torque can provide valuable insights into their interaction and effects in a practical setting. The data obtained from such experiments can be used to develop predictive models and optimize lubrication strategies for specific applications.
A preliminary screening result of AW/EP was used to select the blend for this investigation. Three levels, as well as three key control factors, were considered in this investigation. For the most effective type of lube oil, the control factors and corresponding levels were additive concentrations, with a 4%, 6% and 8%, spindle speed; 3053 RPM, 4580 RPM and 6107 RPM; and a feed rate of 0.1 mm; 0.2 mm, and 0.3 mm. The level of control factor determines values taken for the experiment. The experiment was conducted using the L27 orthogonal array. Under Taguchi's orthogonal array, the proposed investigation entails the conduction of 27 experiments. In order to analyze the variance and single-to-noise ratios of the results, each response parameter was replicated twice, leading to a total of 54 observations. Upon completion of the experiment, the response parameters are transformed into a signal-to-noise (S/N) ratio. Calculated S/N ratio values are compared to the obtained results to see if they deviate from or are close to them. In the present study, the objective function follows a "smaller is better" approach, as indicated by Equation (1). This is because all the response parameters are anticipated to have smaller values for improved performance.

Multiple Linear Regression
The multiple linear regression algorithm (MLR) is a linear regression algorithm used in datasets where the dependent variable is a single one, and the independent variables are multivariate [23]. Linear regression is a type of regression analysis that demonstrates a linear relationship between one or more independent variables and a dependent vari-able. By using multiple linear regression, the dependent variable (y) is modeled via the independent variables (x 1 , x 2 , x 3 . . . x n ) as shown in Equation (2), where, ., x n are the independent variables; y is the dependent variable; n is a positive integer; b is the bias or the intercept; w 1 , w 2 , w 3 , . . ., w n are the weights.

Random Forest Regression
The random forest regression algorithm is a class of machine learning algorithms that combines multiple random decision trees, each of which is trained on a subset of the data set to generate a model [36]. As a result of using multiple trees, the algorithm becomes more stable, and the variance is reduced. Random forest regression is a popular model in the field of regression analysis due to its ability to work well with large and diverse data sets.
The algorithm creates each tree based on a different sample of input data, which are used as input information. The features are randomly selected for splitting at each node, and the trees are run in parallel at each node without any interaction. In order to produce a single prediction, the predictions from all of the trees are then averaged to produce a single result that can be referred to as the prediction of the random forest. A visual representation of a random forest is shown in Figure 3. where, x1, x2, x3, …, xn are the independent variables; y is the dependent variable; n is a positive integer; b is the bias or the intercept; w1, w2, w3, …, wn are the weights.

Random Forest Regression
The random forest regression algorithm is a class of machine learning algorithms that combines multiple random decision trees, each of which is trained on a subset of the data set to generate a model [36]. As a result of using multiple trees, the algorithm becomes more stable, and the variance is reduced. Random forest regression is a popular model in the field of regression analysis due to its ability to work well with large and diverse data sets.
The algorithm creates each tree based on a different sample of input data, which are used as input information. The features are randomly selected for splitting at each node, and the trees are run in parallel at each node without any interaction. In order to produce a single prediction, the predictions from all of the trees are then averaged to produce a single result that can be referred to as the prediction of the random forest. A visual representation of a random forest is shown in Figure 3.

Support Vector Machines
Support vector machines are used to classify data points based on a hyperplane in Ndimensional space (N-number of features), as shown in Figure 4. As a result, several possible hyperplanes can be used to separate the two data points [36]. The objective is to find a plane with the maximum margin, representing the maximum distance between two data points. By increasing the margin's distance, future data points can be classified with greater confidence because they will be reinforced with some reinforcement. Data points are classified using hyperplanes, which are decision boundaries. Depending on which side of the hyperplane the data points fall on, they can be classified into different categories.

Support Vector Machines
Support vector machines are used to classify data points based on a hyperplane in N-dimensional space (N-number of features), as shown in Figure 4. As a result, several possible hyperplanes can be used to separate the two data points [36]. The objective is to find a plane with the maximum margin, representing the maximum distance between two data points. By increasing the margin's distance, future data points can be classified with greater confidence because they will be reinforced with some reinforcement. Data points are classified using hyperplanes, which are decision boundaries. Depending on which side of the hyperplane the data points fall on, they can be classified into different categories.  Furthermore, in the context of feature space, the dimensionality of the hyperplane is determined by the number of input features. When two input features exist, the hyperplane reduces to a line in the feature space. Three input features convert the hyperplane into a two-dimensional plane. As the number of features increases, visualizing the data and the relationships among the features becomes more challenging.

Decision Tree
In its most basic form, a decision tree can be explained as a flowchart showing the clear path that leads to a decision. Figure 5 shows a graphical representation of the decision tree. In data analytics, decision tree is an algorithm used to determine whether data can be classified based on a condition. Initially, a decision tree starts at a point called a node, then branches away from it in two or more directions. Generally, each branch has many possible outcomes, and it is a combination of various decisions and chance events that lead to the outcome. Furthermore, in the context of feature space, the dimensionality of the hyperplane is determined by the number of input features. When two input features exist, the hyperplane reduces to a line in the feature space. Three input features convert the hyperplane into a two-dimensional plane. As the number of features increases, visualizing the data and the relationships among the features becomes more challenging.

Decision Tree
In its most basic form, a decision tree can be explained as a flowchart showing the clear path that leads to a decision. Figure 5 shows a graphical representation of the decision tree.  Furthermore, in the context of feature space, the dimensionality of the hyperplane is determined by the number of input features. When two input features exist, the hyperplane reduces to a line in the feature space. Three input features convert the hyperplane into a two-dimensional plane. As the number of features increases, visualizing the data and the relationships among the features becomes more challenging.

Decision Tree
In its most basic form, a decision tree can be explained as a flowchart showing the clear path that leads to a decision. Figure 5 shows a graphical representation of the decision tree. In data analytics, decision tree is an algorithm used to determine whether data can be classified based on a condition. Initially, a decision tree starts at a point called a node, then branches away from it in two or more directions. Generally, each branch has many possible outcomes, and it is a combination of various decisions and chance events that lead to the outcome. In data analytics, decision tree is an algorithm used to determine whether data can be classified based on a condition. Initially, a decision tree starts at a point called a node, then branches away from it in two or more directions. Generally, each branch has many possible outcomes, and it is a combination of various decisions and chance events that lead to the outcome.

Model Performance Evaluation
Training and testing subsets were divided in 80:20 ratios, widely accepted and wellestablished. Three statistical indices were used to assess the accuracy of the proposed model: mean square error (MSE), mean absolute per cent error (MAPE), and coefficient of determination (R2). These indices provide quantitative measures to assess the performance and predictive capability of the model. The MSE quantifies the average squared difference between predicted and actual values, the MAPD measures the average percentage deviation between predicted and actual values, and the R2 indicates the proportion of variance in the dependent variable explained by the model. The accuracy of any data-driven modeling technique is commonly measured by using all of these indices [31]. The coefficient of determination (R2) value is calculated by using Equation (3).
To determine the performance of a model, the mean squared error (MSE) is usually used per Equation (4).
y P represents the predicted values, y A represents the actual values, and n is the number of samples. It is also possible to evaluate the reliability of the neural network model's result using the mean absolute percentage error (MAPE). MAPE is calculated by using Equation (5): However, MSE and MAPE are both measures of the accuracy of a model, but they differ in how they handle outliers. R-squared is a measure of how well the model accounts for the variance in the data and is not affected by outliers. MSE is more affected by outliers because it squares all errors, while MAPE is more robust since it takes the absolute value of the error. R-squared is not affected by outliers since it is a measure of the variance explained by the model, not the accuracy of the model itself.
In studies involving animals or humans and other research requiring ethical approval, it is crucial to provide details regarding the authority that granted approval and the corresponding ethical approval code. This information ensures transparency and demonstrates compliance with ethical standards.

Tribological Performance
A microscope was used to measure the diameter of the wear scars on the surface of the test ball. In comparison to lube oil A (0.44 mm) and lube oil B (0.40 mm), lube oil C showed the smallest AWSD (0.37 mm) and best AW properties, as shown in Figure 6.

Model Performance Evaluation
Training and testing subsets were divided in 80:20 ratios, widely accepted and wellestablished. Three statistical indices were used to assess the accuracy of the proposed model: mean square error (MSE), mean absolute per cent error (MAPE), and coefficient of determination (R2). These indices provide quantitative measures to assess the performance and predictive capability of the model. The MSE quantifies the average squared difference between predicted and actual values, the MAPD measures the average percentage deviation between predicted and actual values, and the R2 indicates the proportion of variance in the dependent variable explained by the model. The accuracy of any datadriven modeling technique is commonly measured by using all of these indices [31]. The coefficient of determination (R2) value is calculated by using Equation (3).
To determine the performance of a model, the mean squared error (MSE) is usually used per Equation (4).
Yp represents the predicted values, YA represents the actual values, and n is the number of samples. It is also possible to evaluate the reliability of the neural network model's result using the mean absolute percentage error (MAPE). MAPE is calculated by using Equation (5): However, MSE and MAPE are both measures of the accuracy of a model, but they differ in how they handle outliers. R-squared is a measure of how well the model accounts for the variance in the data and is not affected by outliers. MSE is more affected by outliers because it squares all errors, while MAPE is more robust since it takes the absolute value of the error. R-squared is not affected by outliers since it is a measure of the variance explained by the model, not the accuracy of the model itself.
In studies involving animals or humans and other research requiring ethical approval, it is crucial to provide details regarding the authority that granted approval and the corresponding ethical approval code. This information ensures transparency and demonstrates compliance with ethical standards.

Tribological Performance
A microscope was used to measure the diameter of the wear scars on the surface of the test ball. In comparison to lube oil A (0.44 mm) and lube oil B (0.40 mm), lube oil C showed the smallest AWSD (0.37 mm) and best AW properties, as shown in Figure 6.  A combination of TMPTO-based vegetable oils with long fatty acid chains and a proportion of sulfur has a synergistic effect on the performance of lubricants, suggesting that they enhance the performance of the oil itself. A series of extreme pressure tests were carried out until the point of welding was reached. Based on the observations, it has been noted that lube oil A exhibits poor performance under extreme pressure conditions, specifically ceasing at a load of 2000 N. In comparison, lube oil B ceases at a higher load of 6082 N, while lube oil C exhibits even better performance, ceasing at 7848 N. In conclusion, the experimental findings indicate that lube oil C, which incorporates a chemically modified sulfuric additive in conjunction with TMPTO base oil, notably enhances load-carrying capacity under extreme pressure and temperature conditions. This highlights the additive's effectiveness in improving the lubricant's performance in demanding operational scenarios. This modification delivers superior performance and long-term stability, making lube oil C the best option for industrial applications requiring heavy-duty performance.

Preliminary Results of Experiment
According to the previous section, sample C showed the best anti-wear and extreme pressure performance results, followed by samples B and A. An in-depth analysis of sample C in actual contact during the drilling operation was conducted using a Taguchi technique to determine the thrust force, torque, surface roughness, and power consumption.
A summary of the experimental values is shown in Figure 7. It is evident that cutting parameters and cutting forces are correlated. Despite the increased additive concentration, both cutting forces remained minimal and stable. Increasing feed rate increased cutting forces. Based on the results of Wei et al. [37], the change in thrust force concerning spindle speed is minimal or negligible. Therefore, this study's result is consistent with the results we found for drilling aluminum alloys. Furthermore, twenty-seven (27) drilling experiments were conducted utilizing a Taguchi-based orthogonal array design. The thrust force and torque were individually modeled using a polynomial mathematical model, allowing for the independent analysis and optimization of these parameters. A combination of TMPTO-based vegetable oils with long fatty acid chains and a proportion of sulfur has a synergistic effect on the performance of lubricants, suggesting that they enhance the performance of the oil itself. A series of extreme pressure tests were carried out until the point of welding was reached. Based on the observations, it has been noted that lube oil A exhibits poor performance under extreme pressure conditions, specifically ceasing at a load of 2000 N. In comparison, lube oil B ceases at a higher load of 6082 N, while lube oil C exhibits even better performance, ceasing at 7848 N. In conclusion, the experimental findings indicate that lube oil C, which incorporates a chemically modified sulfuric additive in conjunction with TMPTO base oil, notably enhances load-carrying capacity under extreme pressure and temperature conditions. This highlights the additive's effectiveness in improving the lubricant's performance in demanding operational scenarios. This modification delivers superior performance and long-term stability, making lube oil C the best option for industrial applications requiring heavy-duty performance.

Preliminary Results of Experiment
According to the previous section, sample C showed the best anti-wear and extreme pressure performance results, followed by samples B and A. An in-depth analysis of sample C in actual contact during the drilling operation was conducted using a Taguchi technique to determine the thrust force, torque, surface roughness, and power consumption.
A summary of the experimental values is shown in Figure 7. It is evident that cutting parameters and cutting forces are correlated. Despite the increased additive concentration, both cutting forces remained minimal and stable. Increasing feed rate increased cutting forces. Based on the results of Wei et al. [37], the change in thrust force concerning spindle speed is minimal or negligible. Therefore, this study's result is consistent with the results we found for drilling aluminum alloys. Furthermore, twenty-seven (27) drilling experiments were conducted utilizing a Taguchi-based orthogonal array design. The thrust force and torque were individually modeled using a polynomial mathematical model, allowing for the independent analysis and optimization of these parameters.
(a)  Figure 8 shows sample of acquisition for the thrust and torque force obtained from the wireless Spike tool holder dynamometer. Additionally, at the highest spindle speed chosen for this study, a slight increase in thrust force and torque was observed when employing a high feed rate. This can be attributed to the increased number of holes drilled, leading to heightened tool wear due to the elevated friction between the tool and the workpiece. Consequently, higher forces are generated as a result. Adding a chemically sulfurized vegetable oil additive (GS 440L) to   Figure 8 shows sample of acquisition for the thrust and torque force obtained from the wireless Spike tool holder dynamometer. Additionally, at the highest spindle speed chosen for this study, a slight increase in thrust force and torque was observed when employing a high feed rate. This can be attributed to the increased number of holes drilled, leading to heightened tool wear due to the elevated friction between the tool and the workpiece. Consequently, higher forces are generated as a result. Adding a chemically sulfurized vegetable oil additive (GS 440L) to Additionally, at the highest spindle speed chosen for this study, a slight increase in thrust force and torque was observed when employing a high feed rate. This can be attributed to the increased number of holes drilled, leading to heightened tool wear due to the elevated friction between the tool and the workpiece. Consequently, higher forces are generated as a result. Adding a chemically sulfurized vegetable oil additive (GS 440L) to TMPTO base oil significantly increases layer formation between the workpiece and tool.

Analysis of Thrust and Torque Force
As a result of this composition, the trio film remains stable and hydrodynamic lubrication is maintained during drilling. On the other hand, slight effects of additive concentrations in the range of 4-8% seem to be observed on thrust force and torque force. Table 4 depicts the measured response value and single-to-noise ratios (S/N) of thrust force and torque. To determine the proportionate change in the signal-to-noise (S/N) ratio for thrust force, parameter ranking and an analysis of variance (ANOVA) were conducted using MINITAB software (version 21.4.1). These analyses provide valuable insights into the relative importance of the parameters and the significance of their effects on thrust force. Tables 5 and 6 illustrate the results of these calculations.   Figure 9 shows the main effects plots of the S/N ratios for the various parameters of the model that correspond to the various S/N ratios.  Figure 9 shows the main effects plots of the S/N ratios for the various parameters of the model that correspond to the various S/N ratios. The analysis of Table 6 indicates that the most significant factor controlling thrust force is feed rate (73.52%), followed by speed (16.82%) and additive concentration (6.28%) as the least significant factor.  Tables 7 and 8 illustrate the ranking of the parameters and analysis of variance for torque. Figure 10 shows the main effects plots of the S/N ratios for the various parameters  The analysis of Table 6 indicates that the most significant factor controlling thrust force is feed rate (73.52%), followed by speed (16.82%) and additive concentration (6.28%) as the least significant factor. Tables 7 and 8 illustrate the ranking of the parameters and analysis of variance for torque. Figure 10 shows the main effects plots of the S/N ratios for the various parameters   Table 8 indicates that the most significant factor controlling torque force is speed (62.99%), followed by feed rate (23.72%) and additive concentration (7.74%) as the least significant factor.

Prediction of Responses Using Different Machine Learning Models
Four different regression-based machine learning models were used to predict the response parameters of high-speed drilling operations. Figure 11 shows the proposed scheme for predicting response parameters using machine learning techniques. In this session, 54 data points derived from the experimental data were used to create and test the model. The 80% input data were randomly selected to train the model, and the 20% input data were used for testing and model evaluation.  Table 8 indicates that the most significant factor controlling torque force is speed (62.99%), followed by feed rate (23.72%) and additive concentration (7.74%) as the least significant factor.

Prediction of Responses Using Different Machine Learning Models
Four different regression-based machine learning models were used to predict the response parameters of high-speed drilling operations. Figure 11 shows the proposed scheme for predicting response parameters using machine learning techniques. In this session, 54 data points derived from the experimental data were used to create and test the model. The 80% input data were randomly selected to train the model, and the 20% input data were used for testing and model evaluation. Three input variables are used to predict responses by different models, including additive concentration, spindle speed, and feed rate. A standard scalar function was used to scale the training and testing data to minimize errors caused by unit differences in input parameters. The model is cross-validated using GridSearch CV to ensure that the param- Three input variables are used to predict responses by different models, including additive concentration, spindle speed, and feed rate. A standard scalar function was used to scale the training and testing data to minimize errors caused by unit differences in input parameters. The model is cross-validated using GridSearch CV to ensure that the parameters are the best and the model is not overfit or underfit. The model is evaluated through training and test errors to determine whether it is underfit or overfit. Authors compare the four models' predictions for thrust force depicted in Table 9 and Figure 12. Three performance indicators are used to determine the accuracy of the predicted models. Table 10 shows the errors associated with testing and training for thrust force. It was found that decision tree gave the most accurate prediction for thrust force based on the error metrics used, followed by the Tree Support vector machine (SVM), linear regression, and random forest.  Table 11 and Figure 13 represent a comparison between the experimental torque values and the predicted values. According to the performance indicator for torque, as depicted in Table 12, it was found that the Tree Support vector machine (SVM) gave the most accurate prediction for torque based on the error metrics used, followed by Linear regression, decision tree and random forest.      Support vector machines (SVMs) have been established as highly efficient and accurate models. This study's results validate this assertion, confirming SVMs as the optimal choice for regression tasks. However, the selection of the most suitable model depends on the nature of the problem at hand. For linear problems, SVMs tend to yield the highest accuracy. Conversely, when dealing with nonlinear problems, neural networks or random  Support vector machines (SVMs) have been established as highly efficient and accurate models. This study's results validate this assertion, confirming SVMs as the optimal choice for regression tasks. However, the selection of the most suitable model depends on the nature of the problem at hand. For linear problems, SVMs tend to yield the highest accuracy. Conversely, when dealing with nonlinear problems, neural networks or random forests may be more appropriate choices [36,37].
The future scope of work involves conducting a comprehensive investigation into the impact of the TMPTO-based lubricant on thrust force and torque during the highspeed drilling of Al-6061 in order to enhance machining performance, extending tool longevity, enhancing surface quality, increasing energy efficiency, and upholding the overall sustainability of the process.

Conclusions
This experimental study investigated four responses using three key point parameters for Al 6061 material during drilling operations. From this investigation, the following conclusions may be drawn: (1) Incorporating the chemically sulfurized vegetable oil-based additive GS 440L with TMPTO base oil yields substantial enhancements in high-speed metal cutting applications, particularly when working with non-ferrous materials. This additive combination improves the lubricating properties and performance of the cutting process, leading to improved cutting efficiency. (2) The investigated properties reveal an additive concentration range that is optimized when combined with TMPTO base oil. TMPTO base oil, notably enhances loadcarrying capacity under extreme pressure conditions, with a seizer load at 7848 N. (3) Based on the SN ratio analysis and ANOVA ranking of the test parameters, selected test parameters were statistically significant for the response variables. The results of the ANOVA analysis reveal that spindle speed contributes most substantially (62.99%) to torque, with feed rate (23.72%) and additive concentration (7.74%) also showing significant impacts. On the other hand, thrust force is primarily influenced by feed rate (73.52%), followed by spindle speed (16.82%), and additive concentration (6.28%). (4) Experimental and machine learning predictions confirmed strong correlations, and the proposed model offered acceptable results with relatively high prediction accuracy. (5) The forthcoming research will offer valuable insights into the applicability of TMPTObased lubricants in high-speed drilling operations in order to enhance machining performance, prolong tool life, improve surface quality, boost energy efficiency, and ensure the overall sustainability of metalworking processes.