Optimization of fixture locating layout design using comprehensive optimized machine learning

Fixtures are commonly employed in production as work holding devices that keep the workpiece immobilized while machined. The workpiece’s deformation, which affects machining precision, is greatly influenced by the positioning of fixture elements around the workpiece. By positioning the locators and clamps appropriately, the workpiece’s deformation might be decreased. Therefore, it is required to model the fixture–workpiece system’s complicated behavioral relationship. In this study, long short-term memory (LSTM), multilayer perception (MLP), and adaptive neuro-fuzzy inference system (ANFIS) are three machine-learning approaches employed to model the connection between locator and clamp positions and maximum workpiece deformation throughout end milling. The hyperparameters of the developed ANFIS, MLP, and LSTM are chosen using the evolutionary algorithms, including genetic algorithm (GA), particle swarm optimization (PSO), butterfly optimization algorithm (BOA), grey wolf optimization (GWO), and wolf optimization algorithm (WOA). Among developed methods, MLP optimized using BOA (BOA-MLP) reached the highest accuracy among all developed models in predicting the response surface. The developed model had a lower computational load than the final element model in calculating the response surface during the machining process. At the final step, the prementioned five evolutionary algorithms were implemented in the developed BOA-MLP to extract the optimal parameters of the fixture to decrease the deflection of the workpiece throughout the machining. The proposed method was modeled in MATLAB. The outcomes showed that the mentioned model was efficient enough compared with the previous method, such as optimized response surface methodology in the point view of 0.0441 μm lower workpiece deflection.

Acronyms a 1 , a 2 , a 3 , b 1 , b 2 , and b 3 Linear gains of fuzzy Rules 1 and 2 A i and B i Fuzzy sets c, k, a , and I The sensory modality, the fragrance's perceived magnitude and observed power exponent, magnitude reliant on stimulus intensity, and modality in BOA c j t The memory cell of the j th neuron at time t in LSTM c j t The new memory value of the j th neuron at time t in LSTM c t The vector representations of c The forget gate of the j th neuron at time t in LSTM F c1 First clamp force F c2 Second clamp force g Centered logistic sigmoid function h j (n) The j th output node in the n th data point of MLP J and K Coefficients of WOA L The distance between the prey and the best agent m j t The output of the j th neuron in LSTM m t− 1 The vector representations of m The i th raw dataset (combination of input and output) x = x 1 , x 2 , ⋯ , x T Input sequence of LSTM y Actual outputs of MLP y j The previous neuron's output z 1 and z 2 Outputs of fuzzy Rules 1 and 2 Δw ij (n) The variation of weights in the i th input and j th output node in the n th data point of MLP The i th standardized input data Standard logistic sigmoid function ∅ � The derivative of the activation function

Introduction
Machining fixtures are an inescapable part of the work holding devices since they precisely locate and restrict the workpiece throughout the machining process [1]. Optimized machining fixture layout refers to properly selecting the position and orientation of locators and clamps around the workpiece to reduce its deformation throughout machining. Designing and optimizing fixture layouts have a distinctive role in production since they increase production rates and reach a higher quality standard [2]. However, factors including fixture component position, clamping force magnitude, and cutting force are inevitable in determining the optimum machining fixture layout. On the other hand, clamping and cutting forces must be implemented for immobilization and machining of the workpiece accordingly. They cannot be reduced any farther than a particular limit. As a result, workpiece deformation is unavoidable throughout machining, which may result in a defective workpiece. Identifying the proper location for the fixture elements, on the other hand, will reduce workpiece deformation, which would assist in machining accuracy.
While machining fixtures play a significant role in production, designers struggle to optimize the layout of the fixture [3]. Traditionally, fixture layout optimization has depended on the designer's knowledge and ability, which is expensive and laborious [4,5]. The machining inaccuracy and the fixture arrangement do not have a direct analytical connection [6,7]. The complex fixture-workpiece interactions were investigated using numerical approaches, including the finite-element method (FEM) [8,9]. Edward [10] created a quick support layout optimization model using FEM to determine the stiffness of the workpiece and a set of tests to validate the model. Li and Melkote [11] suggested an optimization approach for fixture layout. The resultant nonlinear program was solved using Zoutendjik's technique of possible directions, and the ultimate enhanced solution depends on the initial viable solution. Hamedi [12] proposed the hybrid FFNN and GA to extract the optimal fixture layout. He used the dataset generated via FEM for training the FFNN. Padmanaban et al. [13] demonstrated an ant colony algorithm-based discrete and continuous technique for optimizing the layout of the fixture. Selvakumar et al. [14] used FFNN to estimate the deformation of the workpiece during milling process based on the regenerated data using FEM and Taguchi design experiment. Sundararaman et al. [15] modeled the relations between locator and clamp position to calculate the maximum deformation of the workpiece using response surface methodology (RSM). Afterwards, optimization techniques such as sequential approximation and LINGO were used. Xing et al. [16] optimized the fixture design for complex auto-body components using a nondomination sorting social radiation method. Lu and Zhao [17] used feed-forward neural network (FFNN) and GA to optimize the layout of the sheet metal workpiece fixture utilizing the 4-2-1 locating scheme. The FFNN is employed to estimate the deformation of sheet metal workpieces under various layouts of fixtures and their induced locating error. Lu utilizes a GA and Zhao [17] to identify the optimal location for the fourth fixture locator relying on the FFNN estimation model.
The usage of machine-learning methods and evolutionary algorithms to extract the optimal fixture layout has been focused by lots of researchers in recent years to extract the optimal fixture layout. Low et al. [18] used reinforcement learning method to design the fixture layout fully automatic. Higher computational time and lack of accuracy were the disadvantages of their proposed method. Wu et al. [19] used genetic algorithm to extract the optimal position of the clamping points in blade fixture. As they predict the deformation of the blade during the milling process via FEM, the higher computational load and limited time of evaluation were their proposed method disadvantages. Rezaei Aderiani et al. [20] extract the optimal fixture layout for compliant sheet metal assemblies using GA and FEM. The higher computational load and limited times of evaluation were their model limitations; same as Wu's model [19]. Wu et al. [21] used the hybrid model of FEM and GA to extract the optimal layout of the fixture in end-milling process of the fan blade. The same methodology has been followed by Butt et al. [22] to extract the optimal fixture layout in milling operation. Du et al. [23] used simulated annealing algorithm (SAA) as a global optimization method instead of evolutionary algorithm in order to decrease the number of trials in extracting the optimal solutions. However, they still used FEM for calculation of the workpiece deformation during the milling process, which is very time consuming. Vinosh et al. [24] followed the same methodology of previous studies, which is based on the combination of evolutionary algorithm and FEM. However, they [24] used Taguchi methods for designing of experiments in order to reduce the number of trials. Michael Thomas Rex et al. [25] used the FEM and mixed discrete-integer GA to decrease the deformation of the prismatic workpiece in the pocket milling operation. Alshameri et al. [26] used SAA; same as Du et al. [23] for calculation of the set of points for designing the robust fixture. Hu [27] used hidden Markov model for optimization of digital twin-driven reconfigurable fixturing in aviation industry. Yu and Wang [28] used machine-learning method to extract the model that can imitate the same process using RSM and mathematical optimization. Li et al. [29] used whale optimization algorithm to reduce the clamping deformation of the curved thin-walled parts. Feng et al. [30] proposed the recent machine learning-based prediction model of the workpiece deformation calculator using extreme gradient boosting (XGBoost) method.
Based on reviewed articles, the previous works can be divided into two groups. The first group of the work is established based on the combination of optimization method (evolutionary, global, mathematic) with FEM [10,11,13,16,[19][20][21][22][23][24][25][26][27]29]. These studies are accurate in calculation of the workpiece deformation during the machining operation. However, the FEM is computationally expensive, and each trial takes a couple of hours. Then, optimizer should reach the optimal results using the limited number of trials. On the other hand, the second group of studies are based on the combination of machine learning developed model and optimization technique [15,17,18,28,30]. Using these developed methods based on machine learning technique, the recalculation of the workpiece deformation can be achieved in couple of seconds. It should be noted that the idea of combining machine learning and optimization (second group) is limited in using some basic machine learning techniques including RSM, FFNN, and XGBoost. However, the comprehensive study with evaluations of many machine-learning methods and evolutionary algorithms is not investigated in this field.
The main contribution of this study is to fill this gap in the optimal designing of the optimal machining fixture. Cost-effective modeling must determine the connection between fixture components and the maximum workpiece deformation. Furthermore, the evolutionary technique's implementation will be streamlined to minimize the complexity and computing time of the machining fixture layout without compromising solution quality. As a result, multilayer perceptron (MLP), long short-term memory (LSTM), and adaptive neuro-fuzzy inference system (ANFIS) are employed to predict the deformation of the workpiece based on the different 3-2-1 locator's position of the machining fixture. The hyperparameters of the developed ANFIS, MLP, and LSTM models are chosen via the evolutionary algorithms, including particle swarm optimization (PSO), GA, grey wolf optimization (GWO), butterfly optimization algorithm (BOA), and whale optimization algorithm (WOA). The proposed optimized machine-learning methods are compared with each other, and the best candidate is selected as a proper representer for calculating the workpiece elastic deformation during the machining process. In addition, the selected model is employed in evolutionary algorithms to extract the optimal 3-2-1 locator's position of the machining fixture to achieve the minimum workpiece deformation.
In the next section, the design of the fixture layout with the addition of FEM is explained in detail. Section 3 explains the proposed methodology based on the regenerated data using the FEM model. It includes the explanation of ANFIS, MLP, LSTM, and evolutionary algorithms. Section 4 describes the implementation of the proposed method in MATLAB software to prove the efficiency of the proposed method compared to the previous method. The conclusion is remarked in Sect. 5.

Fixture layout and finite element simulation
The position of fixture elements has a substantial impact on workpiece deformation throughout machining and clamp actuation. As a result, it is critical to identify the proper position for the locators and clamps to reduce the workpiece's maximum deformation. The model proposed by Sundararaman et al. [31] was used in the present study to evaluate the capabilities of the developed optimization methods in calculating the most proper positions of the locating and clamping agents in the fixture. In the mentioned research, a 2D simulation was performed using finite element analysis software to calculate the machining forces acting on a deformable part. The workpiece in the simulation was a 300 × 90 mm rectangular sheet with a thickness of 10 mm. The fixture in the present study is a 2D machining fixture with a 2-1 locating layout. This specific layout needs two clamping forces to push the workpiece toward the locators. The direction of the tool axis is perpendicular to the plane of the workpiece directed inward. In other words, the represented tool in Fig. 1 is the cross-section of the 3D endmill tool. It rotates clockwise to apply the two machining force components to the workpiece. The workpiece was steel with the modulus of elasticity and Poisson's ratio equal to 206 GPa and 0.3, respectively. At least two clamping forces must also act on the opposite edges of the locators so that the workpiece stays stable during machining. The chip removal effect is accounted for through the element death method [31]. The displacement of the clamping and locating points L 1 , L 2 , and C 2 along the y-axis is zero, whereas L 3 and C 1 along the x-axis are zero. The machining forces were calculated in two dimensions through performing simulations. The mentioned method suffers from several shortcomings. First, it does not utilize the experimentally validated machining force models to calculate the three-dimensional machining loads and torque. Second, even in the 2D simulation, the machining torque was not considered in the model, which strongly affects the locating layout. Finally, it would have been better if the machining model had been considered in three dimensions closer to the actual machining processes. The authors have identified these deficiencies, and comprehensive research is underway to address these deficiencies, which will be reported in future articles. The optimization parameters include the position of the locators and clamps' positions concerning the workpiece vertices. The objective function is to minimize the maximum amount of elastic deformation of the workpiece under the clamping and machining forces. The clamping forces were considered as F c1 = 200 and F c2 = 350 for the end-milling process. The machining forces were calculated as 100 N and 286 N in the x-and y-directions, respectively. The clamp C 2 and locator L 1 are horizontally measured from the left top and left bottom corner of the workpiece. In contrast, clamp C 1 and locator L 3 are vertically measured from the workpiece's right bottom and left bottom corners. In addition, as illustrated in Fig. 1, the locator's position L 2 is determined by horizontally measuring from the right bottom corner. A horizontal step milling operation is conducted on the part. An endmill is used to perform this machining process. To convert 3D to 2D machining, it was assumed that the cutting flutes of the tool are parallel to the tool axis such that no axial force was applied to the workpiece. Therefore, only two cutting force components in radial and tangential directions were applied to the workpiece (Fig. 2). Depending on convergence analysis, the machining route is discredited into 11 load steps, and the maximum deformation for each load step is obtained. The required result for a single tool path at the final milling operation is the biggest eleven. The geometry of the 2D workpiece was being discretized into 42 quadrilateral plane elements, containing four nodes, each of which has two degrees of freedom. Figure 3 depicts the FEM model of locator, clamp, and machining forces at various load stages. Locators with positive reaction forces keep the workpiece connected with all the locators throughout machining, while those with opposing reaction forces induce the workpiece to be decoupled from the locators. External forces are imposed via clamps, and the limitation of positive reaction forces at the locators during machining is guaranteed.
All experiments in the study were conducted using ANSYS and MATLAB. The random location for the clamps and locators in the required range was created in a MATLAB m-file to conduct exploratory tests. The maximum deformation of the workpiece is computed while considering the limitations and saved as an output text file. The design variables and the maximum deformation of the workpiece are subsequently coded under MATLAB.

Methodology
This paper aims to estimate the workpiece's elastic deformation throughout machining inside the planer fixture with a 2-1 locating scheme. Figure 4 shows the schematic proposed method in the current research as the main contribution, which is the combination of machine-learning methods including ANFIS, MLP, and LSTM with evolutionary algorithms including GA, PSO, GWO, BOA, and WOA. According to Fig. 4, the process starts with selecting different locating and clamping points. Then, the FEM using ANSYS software is designed and simulated based on these predefined parameters. After simulation, the maximum deformation of the workpiece is captured in the FEM environment to be used in the proposed machine-learning methods. Then, three machine-learning methods, as the most common methods in this field (ANFIS, MLP and LSTM), are chosen to be trained based on the calculated data on FEM software. It should be noted that the evolutionary Fig. 2 The 3D and projected 2D end-milling machining process algorithms are used to extract the optimal hyperparameters of the developed machine-learning methods and extract the optimal fixture locating and clamping positions using developed machine-learning methods.

Data recording
The design parameter range should be as narrow as possible. It is due to the optimization point of the method to increase   x y the convergence speed of the evolutionary algorithm. In addition, narrowing the search area of evolutionary optimization algorithms can decrease the chance of local optimal point problem. The areas with a significant impact on the response of maximum deformation of the workpiece should be identified reasonably. The initial tests are designed and executed using a random and systematic search strategy to discover the best area with the least amount of workpiece deformation and determine the possible range for the design indices. Random trials consisting of ten sets with experiments of ten, twenty, thirty, and one hundred were prepared and performed in reference [31] for this configuration. The most significant workpiece deformation was calculated for each trial. Table 1 shows the maximum workpiece deformations as a function of the design parameters. The connection between fixture element position and the maximum deformation of the workpiece is simulated in the area with the most significant potential to discover the optimum solution. Table 1 indicates that the associated locator and clamp position values for the two successive minimums of maximum workpiece deformation offer the ideal potential range for the design parameters. In the design input parameter L 1 , the range is set initially at 5 to 148 mm. This generates the optimum outcome and goes narrower because of the initial trials, ranging from 105.1 to 119.4 mm. Likewise, ranges are generated for additional design input parameters. The experiment plan used is a second-order rotatable design with 10 star points, 10 center points, and 32 factorial points (52 experimental situations). As mentioned before, the maximum deformation of the workpiece is determined for the appropriate number of tests using FEM simulation, and the related variables of design are saved in MATLAB. The recorded variables of design and the APDL input batch file are transmitted to ANSYS to estimate the maximum deformation of the workpiece.

Data pre-processing
Three tasks should be implemented in the data before developing the model and applying the dataset. Initially, the data out of range should be removed to increase the robustness of the system. The second step is related to the normalization or standardization of the system to decrease the data complexity for a system before training. In this work, both methods are employed to decrease the complexity of the network's input data and increase the system's accuracy. The following formula is used to calculate the data standardization: where x i and x i are respectively the i th raw and standardized dataset. Dataset is defined as the combination of input data x (fixturing and clamping position) and output data (elastic deformation of the workpiece during the end-milling process). x and x are also the functions for extracting the standard and average deviation of the data.
The normalized data can also be obtained as where n x i is the i th normalized dataset. x _ and x are also the functions for extracting the minimum and maximum values of the data. At the final network pre-tuning process, the data is divided into 86.54% and 13.46% for the testing and training process of the network. The testing data is not shown to the system until the testing stage of the network in order to reach realistic results via the proposed models. Three models, including ANFIS, MLP, and LSTM, were utilized in this paper to estimate the elastic workpiece deformation during machining based on the locating and clamping positions of the fixture.

Adaptive neuro-fuzzy inference system
Being a Takagi-Sugeno-type artificial intelligence model used to solve nonlinear and complicated problems, ANFIS is a hybrid system composed of an artificial NN and a fuzzy inference system (FIS) [32]. This implies that ANFIS is capable of self-learning and reasoning operation. ANFIS's fuzzy rules are as below: In Rules 1 and 2, the variables v and u are inputs. The outputs z 1 and z 2 are specified by Rules 1 and 2. The parameters a 1 , a 2 , a 3 , b 1 , b 2 , and b 3 are acquired during the learning process. As seen in Fig. 5, the ANFIS structure comprises five layers, each with one output and two inputs.
Each layer has the following meaning: The first layer is a layer of fuzzification that employs the functions membership in computing fuzzy clusters according to the training input data. The indices a i and b i are utilized to establish the membership function and calculate the degree of membership as in which membership degrees are u and y . Fuzzy sets are A i and B i , and x is the input, whereas Q j i is the i th node in the j th layer output.
(2) Table 1 The maximum deformation of the workpiece and experimental design matrix is estimated by FEM [31] Exp. no The rule is layer 2, which is utilized to double the firing power. The output should be calculated as Layer 3 is for normalizing the last layer's firing intensity.
The normalized values are generated by dividing the i th rule's firing strength by the summation of all firing strengths: Layer 4 is the layer of defuzzification. The result is computed as the normalized firing strength and the parameter set dot product ( p i , q i , and r i ) as The final result of layer 5 is created by adding the defuzzification outputs of each rule: The error signals from the output to the input layers are computed using the gradient descent approach's backpropagation algorithm. The algorithm attempts to minimize the training error by adjusting the customizable parameters.

Multilayer perceptron
Originated by Rosenblatt [33], the MLP is an FFNN, and a single hidden layer MLP is known as vanilla [34]. The simplest MLP comprises three layers: output, hidden, and input. The neurons use nonlinear activation functions except the input nodes.
MLP can be used for classification and regression problems via supervised training. The inputs and outputs are used during the training process to adjust the network parameters, including weights and biases, to minimize the error between target and predicted outputs. The backpropagation algorithm is employed to derive the optimal weights and biases regarding the root mean square error (RMSE), mean square error (MSE), or other measures. The stochastic gradient descent technique is used in backpropagation through a backward pass of weights and biases through the MLP. During the training process, the j th output node in the n th data point can be represented as where d and y are the targets and actual outputs. The weights can be calculated by minimizing the error using the entire network output, as follows: The variation of weights using the gradient descent method is where y i and s i are the last neuron's output and local induced field, respectively. In addition, is the learning rate that affects the convergence of the MLP. The weight derivative concerning the local induced field can be simplified for an output node as where ∅ � is the derivative of the activation function, which is constant against weight and local induced field. The backpropagation of the activation function method is used to vary the output layer weights by the hidden layer weights [35].

Long short-term memory
Hochreiter et al. [36] had demonstrated the effectiveness of the LSTM to address these disadvantages of the RNN. The gradient information eliminates or destroys if a highly long sequence is employed for learning. As an improved version of the RNN, the LSTM model can handle a long sequence of data samples using its memory blocks. The output is generated by a fully connected layer and a regression layer. It consists of three gates: input, output, and forget gates. Each gate is composed of the dot product and sigmoid function to protect the gradient information against distortion or elimination and control the information flow. Figure 6 presents a memory block depicting f multiplicative gating units to determine the information flow and memory cells to connect and memories the secular state.   ). Then, the outdated information flow resets, and the CEC weight is substituted with the multiplicative forget gate activation. After updating the memory cell based on the new memory value, the forget gate, f j t , is computed as where W f , U f , and V f are the diagonal weight matrices. The same methodology is adopted in the input gate, which determines the reserved new features as follows: where W i , U i , and V i are the diagonal weight matrices. It should be noted that the value of the three gates is between 0 and 1. The LSTM output is formulated as where g is a centered logistic sigmoid function within the range [−2, 2] , i.e., Training of the LSTM is based on a modified real-time recurrent learning (RTRL) and a truncated backpropagation through time (BPTT) along with the gradient descent optimization method. The loss function is explained as the sum of square errors. The memory cell shortens the errors by exploiting the linear CEC of the memory cell. Inside the CEC, the error recedes and is discharged from the cell in a degraded exponential manner. This is the primary capability of the LSTM in dealing with a long prediction horizon compared with the RNN.
The hyperparameters of the three investigated machinelearning methods, namely, ANFIS, MLP, and LSTM, are chosen based on five evolutionary algorithms in the last five subsections of this section.

Genetic algorithm
The GA belongs to the group of evolutionary-based algorithms based on the evolutionary algorithms' natural selection. Recently, GA has been utilized to solve a wide range of optimization problems. It is a bio-inspired operator that uses mutation, crossover, and selection. Mitchell [37] devised this approach according to the natural selection process. Genetic algorithms often start with generating a randomly produced population of chromosomes. Following that, the created chromosomes are "assessed" using an objective function, and those that best approximate the "ideal solution" to the issue have a higher probability of reproducing.

Forget Gate
Memory Cell Fig. 6 A memory block depicting the data flow at time step t GA parameters must be carefully chosen to maximize the accuracy and speed of the final result's convergence. The crossover and mutation parameters must be adjusted correctly to ensure the success of the GA and the acquired output. The high mutation rate raises the danger of missing a solution close to the present state. Setting the mutation parameter to a small value raises the likelihood of being trapped in the local optimum. In addition, crossover prevents kids from becoming identical replicas of the parents from the prior generations. The mutation parameter had to be set to a fair value and the crossover parameter to a significant value, as suggested by Mitchell [37]. In the present study, the population size is 80. The small population size imposes a constraint on the GA search capabilities. In addition, the large population increases optimization length without significantly improving the outcomes. The GA results will reveal that the mutation, crossover, generation, and population parameters have all been selected appropriately in this paper.

Particle swarm optimization
Kennedy et al. [38] developed PSO as a swarm-based algorithm category of the evolutionary algorithms. PSO is a stochastic population-based optimization method via the collaborative behavior of a different colony of species known as populations [39]. It is similar to evolutionary computation methods such as GA [37]. It starts by defining the initial particles to evaluate the objective function at every particle's location. It then decides a minor objective function as the best location. It calculates the new velocities according to the best position of the practices and their neighbors and the current velocity. The mentioned procedure continues iteratively until the algorithm fulfills the stopping criteria. PSO is capable of improving the candidate solution of an optimization problem iteratively.

Grey wolf optimization
In comparison to a GA [37] or PSO [38], Mirjalili et al. [40] proposed a GWO method that requires fewer parameters to be predefined. The adjustment of controller parameters is straightforward, with a balance between exploitation and exploration. GWO imitates grey wolves' social hunting behaviors. A grey wolf method is divided into four distinct groups: alphas (the search agent/solution of the highest rank α), betas (the second rank β), deltas (the third rank δ), and omegas (the fourth-ranked search agent ω). Omegas should lead all three ranking agents, while all other solutions are omegas. GWO's mathematical relationship is dependent on these four wolves' groups in prey hunting. A fitness function is necessary for GWO to estimate the optimal hyperparameters of the presented methods.

Butterfly optimization algorithm
Arora and Singh [41] were the first to suggest BOA, which relies on butterfly behavior in nature. Compared to wellestablished evolutionary algorithms, including GA, which retains qualifying solutions, BOA does not exclude any solutions from the search space, implying that every solution seems to have an equal opportunity to improve with the novel solution. Another notable distinction between BOA and typical evolutionary algorithms is that BOA computes the value of the fitness function for creating original solutions. In BOA, the butterflies act as search agents, moving randomly toward each other and generating further novel solutions. BOA is divided into three broad steps: initialization, iteration, and last phases. The approach defines the solution space, cost function, and boundary variables of the issue during the initialization stage by generating a population. In the second stage, BOA uses an iteration technique to identify the position of butterflies at random, computing and storing fitness values. Ultimately, BOA reaches the stop in the final stage, discovering the optimal solution to the optimization issue. From a mathematical standpoint, the butterflies act as search agents for the BOA and generate fragrance at their positions utilizing k = cI a , where c , k , a , and I represent the sensory modality, the fragrance's perceived magnitude and observed power exponent, magnitude reliant on stimulus intensity, and modality, respectively. The algorithm comprises two critical stages: local and universe searches.
Butterflies use techniques mentioned in the local and universe phases to discover activities, including food and partners. The necessities for food and partner have a noticeable fraction chance of p in such activities due to the physical proximity and several other elements, including rain and windy weather. The probability is used in this algorithm to find common approaches between the local and universe investigations in the problem's search space. BOA determines the best fitness employing the most outstanding solution discovered in the issue by satisfying the stopping criteria. It is worth noting that BOA is launched with the following parameters: c = 0.01 , a = [0.1, 0.3] , and p = 0.5.

Whale optimization algorithm
Mirjalili et al. [36] suggested WOA based on the behavior of whales. The three operators in this algorithm are bubble-net foraging behaviors, prey search, and encircling prey. Whales detect the location of their prey and encircle it. The best hunting location is not shown right at the start of the hunt. The remaining agent whales will then gravitate toward the most effective search agent whale. The whales' reported behavior was mathematically modeled as follows: where J and K denote coefficients, t is the current iteration, and p i and P best represent the position matrix for the i th whale with the number of search agents and the best search agent's position vector, respectively. Moreover, L denotes the distance between the prey and the best agent. P best should be replaced with a better option in each iteration. The coefficients J and K may also be computed using the following formula: where 0 < r < 1 is a random number that can be acquired using Eq. (26). The two basic mechanisms for bubble-net behavior used by humpback whales are the shrinking encircling mechanism and spiral updating position.

Results and discussions
Three models optimized by five evolutionary methods have been investigated in this study to estimate the elastic deformation of the workpiece throughout machining based on the different positions of the locating and clamping points. Moreover, the mentioned method is compared and validated with the previous work published in this topic using RSM [1] with the combined GA and PSO methods. The investigated methods are shown by GA, PSO, GWO, BOA, and WOA for optimized method and ANFIS, MLP, and LSTM for machine-learning methods. For instance, PSO-MLP refers to the MLP machine-learning method, which is optimized by the PSO algorithm. The models are developed in MAT-LAB software to be evaluated and prove the efficiency of the proposed methods. This section consists of two subsections regarding the verification and validation of this study.

Verification
The considered models in Sect. 3 with five optimization methods are developed in MATLAB. ANFIS is proposed according to the described model in Sect. 3.3. MLP is developed according to the described model in Sect. 3.4. LSTM is designed according to the described model in Sect. 3.5. Then, the models are optimized based on the five optimization methods in Sect. 3.6-3.10.
The ANFIS function of MATLAB is used to design the ANFIS model. Also, the feedforwardnet and trainNetwork function of MATLAB is used to design the MLP and LSTM model. The proposed algorithms in Sect. 3 are developed under the MATLAB coding environment, and all the results are extracted using the plot function. In addition, the GA and PSO optimization methods are developed using the ga and particleswarm functions of MATLAB. Also, GWO, BOA, and WOA are employed via the provided code in FileExchange [31][32][33]. Five optimization methods are used to extract the optimal hyperparameters of the proposed ANFIS, MLP, and LSTM. Initially, the evolutionary algorithms are implemented in ANFIS, then the extracted hyperparameters using GA, PSO, GWO, BOA, and WOA  Table 2. In the second step of the research, the investigated evolutionary-based optimization algorithms are implemented in the MLP model, and the extracted optimal hyperparameters are reported in Table 3 using GA, PSO, GWO, BOA, and WOA, respectively. Lastly, five investigated evolutionary-based optimization algorithms are implemented in the LSTM model. The extracted three optimal hyperparameters, including the number of layers, number of units, and learning rate, are represented in Table 4. In the next step, the extracted ANFIS, MLP, and LSTM using the extracted optimal hyperparameters are trained and tested to find the most appropriate model in estimating the workpiece deformation during the machining process.

Validation
After extracting the hyperparameters shown in Tables 2,  3 and 4, the models are trained using 86.54% of the represented dataset in Table 2. Then, the models are tested using 13.46% of the dataset, and the results are reported in Table 5. Table 5 reports the results regarding MSE, RMSE, normalized root mean square error (NRMSE), mean, and Std between the calculated workpiece deformation during end milling via FEM and proposed methods. In addition, it represents the correlation coefficient (CC) and R 2 between the calculated workpiece deformation during end milling via FEM and proposed methods.
Based on the represented data in Table 5, among the ANFIS models, BOA is the best evolutionary-based optimization method in the extraction of the optimal hyperparameters with the lowest deformation of the workpiece during the machining process. In addition, BOA is the best evolutionary-based optimization method among MLP models in extracting the optimal hyperparameters with the lowest deformation of the workpiece during the machining process. Lastly, WOA is the best evolutionary-based optimization method in calculating the workpiece deformation during the machining process. After evaluating the models, the BOA-MLP is selected as the most suitable method for calculating the workpiece deformation inside the 3-2-1 located workpiece during the machining. Figure 7a-b represents the results of the extracted BOA-MLP during the testing process as the best-extracted model for calculation of the workpiece elastic deformation during the machining based on the different locating and clamping positions. Figure 7a shows the actual and estimated workpiece elastic deformation during the machining based on the different locating and clamping positions using BOA-MLP during the testing stage. Figure 7b shows the error between the actual and estimated workpiece elastic deformation during the machining based on the different locating and clamping positions using BOA-MLP during the testing stage. Figure 7c presents the error histogram of the BOA-MLP for calculating the workpiece elastic deformation during the machining based on the different locating and clamping positions during the testing stage. Also, Fig. 7d presents the regression between the actual and estimated workpiece elastic deformation during the machining using BOA-MLP based on the different locating and clamping positions during the testing stage.
The extracted BOA-MLP is a candidate as the best model for representing the workpiece elastic deformation during the machining based on the different locating and clamping positions. Then, the five explained evolutionary-based optimization methods used this extracted model to calculate the best locating and clamping position of the fixture layout to reach the lowest elastic deformation of the workpiece. In addition, the previously proposed method using the GAand PSO-RSM model [1] is investigated and compared with our recently proposed model to prove the efficiency of our proposed method.
The convergence plot of the five investigated evolutionary algorithms, including GA, PSO, GWO, BOA, and WOA, is shown in Fig. 8a-e. Based on the represented results, PSO, GWO, and WOA can reach the results more accurately and faster than GA and BOA. Also, the BOA has the slowest convergence speed between all investigated evolutionarybased optimization methods. GWO and WOA can extract the optimal solution after three iterations. Table 6 shows the mathematical representation of the results in Fig. 8a-e with the addition of the previous method results using the GAand PSO-RSM method [1]. Based on the represented results Table 6 The results of five investigated evolutionary algorithms via the BOA-MLP model and the previously proposed model using RSM [31] Method L 1 (mm) L 2 (mm) L 3 (mm) C 1 (mm) C 2 (mm) Elastic deformation (μm) in Table 6, the elastic deformation of the workpiece during the machine of the workpiece using the extracted results of our proposed method is 0.0441 μm lower than the previous GA-and PSO-RSM method.

Conclusion
Fixtures are used in machining operation in order to locate the workpiece in precious location. Extracting the optimal fixture layout is not an easy task in order to reach the highest possible accuracy during the machining of workpieces. Previously, researchers used optimization methods, machine learning, and FEM models to extract the optimal fixture layout. However, there is a lack of comprehensive studies in evaluation of the fixture layout design using machine-learning and optimization methods. This study presented a comprehensive evaluation for optimal designing of the machining fixture using three standard machinelearning techniques, including ANFIS, MLP, and LSTM, which were optimized using five evolutionary algorithms including GA, PSO, GWO, BOA, and WOA. Also, the outcomes of the previously presented model using RSM, GA, and PSO were mentioned in this study for validation purposes. The investigated models were trained and tested using 86.54% and 13.46% of the FEM simulation software dataset. The prediction accuracy of the mentioned models was investigated using various parameters, including CC, MSE, RMSE, NRMSE, and R 2 . Among the investigated methods, BOA-MLP reached the most accurate result based on the regenerated dataset in ANSYS software for the workpiece elastic deformation desiring the machining inside the fixture. The extracted accurate model was selected as the best accurate model to be employed inside the evolutionary algorithms to extract the optimal locating and clamping position of the 3-2-1 fixture during the machining. The results showed the 0.18% lower elastic deformation of the workpiece inside the 3-2-1 fixture during the machining using our proposed method compared to the previous GA-and PSO-RSM method.