A novel optimized design of a piezoelectric-driven 4-stage amplified compliant microgripper using a 2-step multi-objective algorithm

Advancements in microscale technologies have prompted a demand for high precision micro-manipulation. Microgrippers are the primary means of conducting micro-scale operations, and they significantly affect the procedure's performance. This paper presents a novel optimized design for compliant microgrippers, intending to enhance functionality and durability. The mainframe of the proposed microgripper is based on a compact flexure-based compliant structure with four stages of movement amplification. Experiments were designed based on the L25 Taguchi orthogonal arrays. The experiments were conducted using the finite element method in Abaqus 6.14 workbench. Range of motion and maximum created mechanical stress are selected as the two fundamental goals of the optimization. A variety of designs are achieved using the proposed algorithm. The use of Analytical Hierarchy Process has led to the presentation of an efficient and well-defined algorithm to perform decisions. The decision process can be performed with regard to specific requirements of various applications. The presented design process of microgrippers has the potential for customized manufacturing for specific applications. Finding correlations between design parameters and outputs (Amplification factor & Stress), using Taguchi's method in design of experiments (DOE). Optimization of dimensional inputs using a multi-objective genetic algorithm process to achieve an optimal Pareto-front instead of a single design point. Selecting the desirable point on the optimal Pareto-front for specific applications using Analytic Hierarchy Process (AHP) to prevent possible decision-making errors. Finding correlations between design parameters and outputs (Amplification factor & Stress), using Taguchi's method in design of experiments (DOE). Optimization of dimensional inputs using a multi-objective genetic algorithm process to achieve an optimal Pareto-front instead of a single design point. Selecting the desirable point on the optimal Pareto-front for specific applications using Analytic Hierarchy Process (AHP) to prevent possible decision-making errors.


Introduction
In the last decades, the rapid development of Micro-electromechanical systems (MEMS) has raised global interest in micro-operations, which are used in a wide variety of applications such as bioengineering [1], drug delivery [2], optics [3], and aerospace [4]. As the importance and applications of manipulating small parts and objects increase, instruments, and tools for carrying out such operations are becoming of great favor to scientists.
Microgrippers are used to carry out micro-scale operations with the desired precision. Using a flexible mechanism in the microgripper lowers the manufacturing costs and omits friction, clearance, and the need for lubrication [5].
Comparing to other types of actuation, piezoelectric stack actuators have advantages such as fast response, high resolution, stable displacement and large output force [17,18].
The output displacement of piezo-electric stack actuators is usually much smaller than the size of objects we aim to manipulate. Displacement amplification mechanisms have been utilized to amplify the input displacement of microgrippers, resulting in a greater range of motion in the gripping jaws [19].
Due to various features expected of Microgrippers, some studies operated a multi-objective optimization on its design. Dao et al. [20] performed a multi-objective optimization based on displacement and frequency. They achieved a set of optimal dimensions using the differential evolution algorithm. Ho et al. [21] performed an optimization based on hybrid Taguchi-teaching learningbased optimization algorithm (HTLBO) choosing displacement, frequency and gripping effort as objectives. Xiao et al. [22] derived an optimal design of microgripper, based on input stiffness, safety factor and amplifier ratio. Using Radius basis functional network (RBFN) multiobjective GA (Genetic Algorithm) optimization method, they derived a set of optimal designs. Ho et al. [23] carried out an optimization using hybridization algorithm between adaptive neuro-fuzzy inference system (ANFIS) and Jaya algorithm in order to achieve an optimal displacement and resonance frequency. Grossard et al. [24] performed a flexible building block method optimization on microgripper to achieve a higher stroke amplification and force amplification. Nguyen et al. [25] optimized a sand crab-inspired compliant microgripper based on the amplification factor and the maximum stress created by a specific displacement. Das et al. [26] used computational methods to optimize a piezoelectric microgripper design to minimize parasitic motion and increase output displacement. Their experimental results demonstrated low parasitic motion and high precision motion resolution of their microgripper.
Based on what has been said, an overview of the literature shows three major weaknesses: First, due to the use of the full factorial method, the number of experiments to achieve the desired result will be very large. This is a complicated matter in cases involving theoretical analysis and costly and time consuming when conducting practical experiments.
Second, most research has done the process of optimizing the micro-gripper to reach a design point. But due to the wide range of applications of microgrippers, we will need different outputs. For example, sometimes we need to withstand high stress and sometimes a large amplification factor. Therefore, instead of achieving a single design point, it is better for the user to have a set of designs so that he/she can choose the desired design point according to his/her needs.
Third, in cases where the choice of the desired design point is left to the user, unsystematically doing of the selection process leads to decision errors due to incorrect feeling of the demands. Therefore, using a systematic multi-criteria decision-making process is necessary to solve this problem. Of course, this in itself can complicate the process. So, the decision-making method should be effective and simple at the same time.
To address the first weakness, experiments based on the full factorial method are replaced by the Taguchi orthogonal arrays, which, significantly reduce the volume of records while covering the entire data. In addition, to determine the correlation of the parameters, low effect values are eliminated and lead to a concise, agile, and effective regression equation.
For the second, instead of using single-objective optimization, a multi-objective meta-heuristic algorithm based on the genetic algorithm is used, and therefore a single optimal point is replaced by an optimal Pareto-front that provides output flexibility for different applications.
To solve the third problem, the "Analytical Hierarchy process" (AHP), which is a multi-criteria decision-making method, allows the systematic selection of the desired design point according to the real expectations of the user based on the matrix of pairwise comparisons.
It is important to note that the focus of this research is not on a specific design or specific output parameters, but is to provide a fast, flexible and efficient algorithm, yet simple and low cost to optimize the parameters of a microgripper. Therefore, the selected outputs such as amplification factor, stress, vibration modes, dynamic or thermal considerations etc., the method of identification of correlations and validation is not important and done only using the finite element method. Of course, this algorithm can be generalized to the experimental data.
In this research, we present two important factors of a microgripper design, which can help us choose a design for our specific kind of application. These two factors are the amplification factor (AF) and minimal stress. For a set of designs that have the same overall structure but vary in some features, there is usually a tradeoff between these two factors. When we increase the amplification factor, we observe an increase in the maximum stress created in the design and vice versa. Therefore, in this study, a novel method of obtaining the right design for one's needs was developed, utilizing a multi-objective optimization procedure and an importance criterion. Using this criterion, one can decide the relative importance of amplification factor and minimal stress, which will help choose the right design. This method of having minor variations in designs and choosing from them becomes especially important since for every design, achieving a well-regulated control system needs a lot of experiment and design knowledge. However, if we use a universal design with minor regulations to our work, we already know how to best control it and where to put our sensors. Flexibility and efficiency are expected from a universal design. As a result, having a great range of motion and easy means of producing movement were two important factors in our design. To ensure a great range of motion, we used four stages of amplification, and we achieved a remarkable magnitude of amplification factor.
This study aims to investigate the optimization of a four-stage microgripper based on the range of motion and stress responses. Following the introduction, the design process and the computational procedure used to perform the analyses are presented. Afterward, the design of experiments and the mathematical modelings are explained. Results are subsequently presented to demonstrate possible optimized designs. In addition, the observations on the experimental results are discussed further. Finally, conclusions are drawn based on the results and discussions.

Design of the microgripper and the finite element analysis
A main view diagram of the proposed microgripper is presented in Fig. 1. The microgripper is comprised of two stack piezoelectric ceramic actuators (SPCA), fixing holes, preload bolts, bridge amplifiers, lever amplifiers, a number of circular and rectangular flexure hinges, and finally, a pair of gripping jaws. By going through elastic deformation, flexure hinges transmit the produced displacements of the SPCAs to create the desired gripping motion. The symmetricity of the structure improves accuracy, helps balance the structure's internal stress, and doubles the displacement amplification. Preload bolts are used to adjust the preload force on each SPCA. One end of each SPCA is fixed with the displacement transmission mechanism (DTM) using the preload bolt. Several fixing holes on the microgripper are used to keep the microgripper in place. Figure 2 presents the main design parameters of the proposed microgripper design, which are also specified in Table 1.
The piezoelectric actuators produce a small output displacement; hence, four stages of displacement amplification were incorporated to ensure a great gripping stroke in the microgripper, leading to a large tip displacement. A bridge amplifier and three lever amplifiers were deployed to increase the output displacement.
Compliant mechanisms move solely by the deformation and flexure of hinges. Therefore, friction and unwanted motions are eliminated in compliant mechanisms. As a result, the microgripper's controlling system can possess a greater movement precision in such designs.
As for the material, the Aluminum alloy was used due to its favorable mechanical properties such as relatively low density, high flexibility, low cost, and being easy to be machined using electric discharge machining (EDM). Therefore, 7075 aluminum alloy was chosen, which has a modulus of elasticity E = 71 GPa , a Poisson's ratio = 0.33 , a yield strength s = 455 MPa , and a density = 2810 kg/m 3 . In order to reduce the computational costs, a halfmodel was used in simulations, and an x-symmetric constraint was applied to the symmetry plane. A uniform distributed displacement in the y-direction was input to the  surfaces, which contact the upper and lower surfaces of the piezoelectric (5 µm in the "+y" direction for the upper surface and 5 µm in the "−y" direction for the lower surface). The surfaces of fixing holes were constrained to be fully encased. The displacement of gripping tips (which yields the amplification factor when divided by the input displacement of the piezoelectric) and the maximum stress generated in the microgripper structure were the outputs of the simulation. In order to mesh the samples, linear hexahedral elements of type C3D8R with a size of 0.04 mm were used, and a static-general solver with non-linear geometry was employed. A sensitivity analysis showed that a 0.04-mm mesh was fine enough to eliminate the effect of mesh size on the results. The contours of displacement and stress generated in the microgripper can be seen in Figs. 3 and 4.

Design of experiments using the Taguchi method
The Taguchi method is one of the most powerful and stable methods of designing experiments. Comparing to the Full-Factorial method, the Taguchi method decreases the number of experiments in a drastic way, using orthogonal matrices [27]. The main advantage of the Taguchi method is to reduce price and time by decreasing the number of experiments, eliminating the effects of uncontrollable and noisy factors, and determining the main influencing factors [28,29]. Furthermore, by this method, the interactions of factors on each other are identified. An essential step in the Taguchi algorithm is defining the loss function, then defining the N/A relation function regarding the loss function [30]. The higher the signalto-noise ratio, the higher the ratio of useful information to false information and noise [31]. Based on the nature of the dependent variable, the S/N ratio function can be defined in the following three ways [21-25, 27, 32-34]: 1. Smaller is better: 2. Larger is better: 3. Nominal is better: In these equations, Y i is the i-th observation and n is the number of observations [32].
Note that in all the above steps, a higher signal-to-noise ratio indicates that the variable has a stronger significant effect on the dependent variable than the misinformation and noise at the desired level. Also, the difference between the maximum and minimum signal-to-noise ratio at different levels for an independent variable represented by the "δ" parameter indicates that the variable has a greater effect on the output or the dependent variable [32].

Mathematical Modeling (Regression)
The collected data is analyzed by Minitab version 11 software. Different methods of regression can be used to model, predict and analyze various systems [25]. In this study, in order to find the relationship between the dependent variable and independent variables, first from Simple Linear Regression (SLR) and then Stepwise Linear Regression (SWLR) and then from Full Quadratic Linear Regression (FQMLR) and then Backward Full Quadratic Linear Regression (BFQMLR) is used to analyze the data. Then, using Analysis of Variance (ANOVA) parameters such as R 2 , AdjustedR 2 and Predict R 2 , also from P Value and finally RMSE, modeling methods are compared and the most accurate method is selected [35,36]. The general equation for SLR is as follows: where x will be the predictor and 0 and 1 will be the model parameters (coefficients). Epsilon is a random component of the model that follows an independent normal distribution. The estimating equation of SLR method is as follows: where b 0 and b 1 are the estimated parameters of the model and ŷ is the approximate value of the dependent variable y.
In SWLR method, effective parameters are added one after another while the ineffective parameter are removed on by one. The general and estimated equations of SWLR method are respectively as follows: where, x i , … , x j are the predictors and 0 , … , j are the model parameters. Also, b 0 , … , b j are the estimated model parameters and ŷ represents the approximate value of y.
Then FQMLR was used, the general equation of which is as follows [36,37]: The full-quadratic type of estimation equation MLR is as follows: mations of model parameters and ŷ is the approximate predicted value of dependent variable(output). Then, to implement BFQMLR, in each step, each term of the FQMLR equation that has a P value > 0.05 is deemed as non-significant and is removed from Equation. Therefore, the final equation will have fewer parameters than FQMLR and at the same time will be more consistent with the data. Therefore, in the final comparison to choose the best method, it is enough to compare the ANOVA and Root Mean Square Error (RSME) analysis between SWLR and BFQMLR, and choose one of these two methods [35,36].

Multi-objective Optimization Using NSGA-II:
The NSGA algorithm was first introduced in 1994 by Srinivas and Deb using the ranking and non-dominated sorting process [38]. This algorithm had problems in selecting higher quality individuals among all the individuals raised, problems in elitism and complexity of calculations. Therefore, in order to fix these problems, edited version was introduced in 2002 by Deb under the title NSGA-II [39]. Therefore, NSGA-II is an algorithm with high computational efficiency, fast, Non-elitism preventing, Non-Dominated Sorting, and with less control over shared members to maintain their diversity [34,[40][41][42].
The implementation steps of NSGA-II can be summarized as follows: Generation of the initial population with size N and Gen = 0 2. Calculation of all components of the objective function for all individuals of the initial population 3. Non-dominated sorting based on crowding distance and ranking 4. Selecting a portion of the initial population which have a higher ranking, generating crossover and mutation for making offspring 5. Making a new population with size 2 N by a combination of parents and offspring populations 6. Rankin based on crowding distance and selecting N individuals with the higher rank 7. Checking end condition (it can be meet the maximum allowed generations, achieving the desired accuracy, reaching max. time limit and no improvement after in Maximum pre-defined number of generations).
The solution will consist of a set of points whose input are independent variables and whose output are all components of the objective function (dependent variables). Each point on optimal Pareto-front represents an optimal answer to the problem. It is clear that by changing one of the optimal output points, the other outputs will change. Therefore, it depends entirely on the selector will have to choose which point.

Optimum point Selection; An AHP Approach
There are several ways to select one of the several points on Optimal Pareto-front, all of which depend on the opinion of the selector. To ensure that the real and accurate opinion of the selector in determining the weight of each component of the objective function (dependent variables) is used in constructing the final output, a systematic method based on pairwise comparisons, called analytic hierarchical process, is used as the following [43]: Suppose the number of dependent variables (output) is n, denoted by f 1 to fn , in other words: In this case, we form a pairwise comparison matrix as follows: Where a ij represents the comparative significance of output f i over output f j , determined on a predetermined scale (usually from 1 to 5). It should be known that a ii = 1(i ∶ 1 to n) and we have: To calculate the weight of each output component we have: where W i will be the weight of f i . For normalized mode we have: where, w N i is the normalized weight of i-th output (or f i ). Now for each row of optimal Pareto-front, the final combined output can be expressed as y c and defined as below: Now we just have to find the maximum value of calculated y c values.
Instead of multi-objective optimization mode, it is also possible to use a single-objective optimization mode with the following objective function from the beginning of the optimization process, and just find the optimum point:

The Proposed Algorithm
The proposed algorithm can be divided into two basic steps: Modeling stage, Includes geometric designing of the microgripper, defining independent variables (inputs), defining dependent variables (outputs), design and execution of experiments based on proper definition and selection of Taguchi orthogonal matrices, recording and analyzing results and thus ranking and separating variables with specific effects against variables with uncertain effect or noise (and eliminating them in the process of making a mathematical model), deriving a mathematical model using regression in two different ways, ANOVA analysis and deciding on the appropriate model based on this analysis.
Optimization stage, includes defining the multi-component objective function, modifying the objective function based on the nature of the outputs (so that all need to be minimized), performing the NSGA-II algorithm and (14)

Uncertainties and Limitations
Regarding this fact that we often have to use the Taguchi orthogonal arrays instead of the full factorial method due to cost and time considerations on the one hand and the complexity of the calculations and (sometimes) qualitative variables, on the other hand, the weaknesses of this method may affect our research. This means that the accuracy of the constructed mathematical models may be reduced due to the omission or misinterpretation of the interaction of the variables or the small number of experiments and the discrete levels.
In the FEM, the choice of elements and mesh sizes, geometric parameters of the model, and mechanical properties of the model material can be the primary sources of uncertainty in this study.
To eliminate these shortcomings, the results should always be validated within an acceptable range by performing a sufficient number of tests.

Results
According to the proposed algorithm, the first step is the geometric design of the microgripper. According to the review of the articles mentioned in the introduction section, the proposed plan is as shown in Fig. 1.
Also, the geometric design variable parameters are X1 to X6, which are specified in Fig. 2. The minimum and maximum values as well as the levels of each variable are according to Table 3.
Also, the outputs (independent variables) as mentioned in Sect. 4.1 are amplification factor and stress, which are hereinafter referred to as AF and ST in this article. The ultimate goal is to reach the maximum AF and the minimum ST.
According to the independent variables (inputs), dependent variables (outputs) and their nature, the Taguchi orthogonal matrix will be 6 factors and each will have 5 levels of type L25 and according to Table 3.  Table 4 for both AF and ST. The results are according to Table 5.
At this stage, it is necessary to implement the Taguchi algorithm separately for each of the dependent variables AF and ST and obtain the results. Due to the nature of the dependent variables, the signal to noise ratio (S/N Ratio) function is considered larger is better for AF and smaller is better for ST. The S/N values for AF are shown in Table 6 and its diagram in Fig. 7, as well as the S/N values for ST in Table 6 and its diagram in Fig. 8.
According to Table 6 and DELTA values for the levels of independent variables, at this stage we cannot assume any variable from AF or ST completely noisy and completely remove it, so we leave this to the regression steps. The next step is to obtain a mathematical model of the problem using regression. As explained in the proposed algorithm, for each of the dependent variables AF and ST, we first use the SLR and then the SWLR, then FQMLR and therefore BFQMLR. Since SWLR and BFQMLR are complementary regression functions and are better than SLR and FQMLR, respectively, we provide the results for SWLR and BFQMLR only. The SWLR equation for AF is as (1), the BFQMLR for AF is as (2), and the BFQMLR, SWLR for ST are both identical and consistent with 3 (Table 7).   As previously stated in the description of regression methods, noisy and low-effect nerves that had a P value > 0.05 in ANOVA analysis were removed from the regression equation. Also, the Pareto chart of standardized effects for both AF and ST variables and each of the two modes SWLR and BFQMLR are shown in Figs. 9, 10 and 11. As can be seen, for the ST variable, due to the P value threshold value > 0.05, based on a = 0.15, all nonlinear nerves (and some linear nerves) are eliminated and the SWLR and BFQMLR equations are the same. Based on ANOVA analysis, the values of R 2 , ADJUST R 2 , Prediction R 2 for each of the above modes are listed in Table 8.
Then, according to Eqs. 20, 21 and 22, the validity of each of the mathematical models was checked with the values obtained from FEM, and the relative error and RSME of each model were obtained according to Tables 9, 10, and 11. The results of these tables show that for the AF variable, the BFQMLR model has a higher R 2 , which indicates a better fit of the data. In addition, for these variables, the value of R 2 (adj), which is used to compare models with different number of terms, indicates the considerable superiority of the BFQMLR model. Comparison of superiority in predicting results (excluding input points for Regression) is done by R 2 (pred), in which case there is a significant advantage in the BFQMLR model, as well. On the other hand, the average relative error percentage in this model is lower than SWLR. Comparison of RSME values in this case also shows the superiority of the BFMQLR model. Therefore, in general, the BFQMLR model is selected for the dependent variable AF. Regarding the ST variable, due to the similarity of the SWLR and BFQMLR models, there is no choice between the two modes, but the sum of the above shows that the AF model is more valid than the ST model. In the next step, in order to perform the optimization operation, we define a targeting function with two output components according to the system outputs as follows: According to the nature of AF (large is better) and also the nature of ST (where smaller is better) and that in the optimization process we want to minimize the objective function, it is necessary to multiply the AF variable  Pareto front diagram and the optimal points on it are according to Table 13 and Fig. 12.
Obviously, each of the points on the optimal Pareto front can be used as an optimal point, but if we want to have only one optimal point as the final answer, as stated in the section on optimization description (4.3), we use AHP method and create a pairwise comparison table and then ask an expert user to compare the AF and ST parameters in terms of importance (for a specific application) and fill in a table according to Table 14. In Table 14, for example, the importance of AF relative to ST is three to one. It is also noteworthy that in order to combine the   output components according to the obtained weights, it is necessary for both of these components to be normalized before multiplying by the weights. For example, to normalize AF with respect to the optimal Pareto-front, we do the following: where AF N is the normalized value of AF and AF max and AF min are the maximum and minimum values for AF in the optimal Pareto-front, respectively. After performing these (25) AF N = AF AF max − AF min steps, the table of values of the optimization parameter will be as in Table 15. The values in Y c = combined output are obtained by the following formula: As can be seen, the maximum value of the combined output in Table 15  To show the flexibility of the above method, the above point (p3) along with the optimal points obtained for the AF to ST ratios of 4 to 1 (p4) and 8 to 1 (p8) are shown in Fig. 13.

Discussion
As shown in Table 6 and Fig. 8, the input variables X3 and X1 have a very strong effect on the AF output parameter, while the effect of other variables is less in this case. In addition, the effect of these two variables has a definite upward trend, while the other 4 variables fluctuate. Regarding ST output, it can be concluded from Table 7 and the extent that the X3 and X5 regression models are expected to be more important than the other inputs. In the mathematical model of the AF variable, according to the SWLR method, all neurons that had a p value > 0.05 were removed, and as a result, only 4 variables X1, X2, X3, X5 were present, and the effect of X4, X6 noise was identified and removed. Also, in the AF mathematical model based on BFQMLR, in the linear part similar to the SWLR model, only 4 variables X1, X2, X3, X5 are present. Out of 6 square terms ( x 2 i ), only x 2 2 , x 2 3 , x 2 4 and x 2 5 are present and other cases have a noise effect. Also, out of 30 possible terms for the interactional component, only x 3 x 5 , x 2 x 4 , x 2 x 3 and x 1 x 4 are present, and the other 26 terms have a noise effect and have been removed from the model. On the other hand, in the ST mathematical model with BFQMLR method, there are only 2 linear terms that are due to the variables X3, X5 and other linear terms with all interactional square terms have a noise effect and have been removed, so SWLR regression models, BFQMLR are the same for ST. Comparison of the effects of ANOVA analysis between the two regression models expressed for AF, as mentioned in the Results section, shows the absolute superiority of the BFQMLR model. Validity check performed on the presented regression models and the error values entered for each in Table 12 and 13 indicate the validity of all 3 regression models. Of course, the validity of this model is fully verifiable if it is valid for points other than the points used to Fig. 11 Pareto chart of standardized effects for both SWLR and BFQMLR of St (only for terms which P value < 0.05. Here 1.50 is threshold) build the model. In terms of analysis of variance, such validation is performed using the parameter R 2 (pred.) , Which has an acceptable value for all 3 models according to Table 10. The optimal Pareto front shown in Fig. 12 indicates that achieving higher AF magnitudes (which is desirable) has the unintended effect of increasing stress. For this purpose, it is necessary to compromise between the values of AF and ST, which is a kind of technical contradiction, in a reliable and systematic way. For this purpose, and instead of manual weighting, which is a common method in these cases, the matrix of pairwise comparisons and determination and normalization of the weight of each output using the Analytic Hierarchy Process (AHP) has been used. The results shown in Fig. 13 confirm this claim and its conformity with the logic of the problem.s Finally, Table 16 compares the proposed algorithm with the previous ones (mentioned in the literature review). This comparison is made in 4 contexts. In the fourth column, comparisons are made from the perspective of data collection methods, and as can be seen, some studies have used estimates based on analytical relationships or arbitrary finite points, which certainly do not provide reliable coverage of the range of inputs and their optimal point may be a local extreme. Of course, this is not the case when well-known design of experiments methods such as Taguchi orthogonal arrays or Response Surface Method (RSM) are used. In the fifth column, the optimization methods used are compared, all of which are multi-objective and each has its advantages and disadvantages. Therefore, in this case, neither can be superior to the other. In the sixth column, the methods for selecting the optimal point    from the set of answers are examined, and as can be seen, except for the proposed algorithm of this research, none of the previous cases provides such an answer and limited to a maximum of 5 optimal points (Whose choice was also unsystematic). This can be a measure of the flexibility of the algorithm for different applications of microgrippers and presented in the sixth column.
For a better understanding, if each algorithm has an advantage in a specific aspect, it is highlighted in the corresponding column, and as can be seen, the only algorithm that has an advantage in all four areas at the same time is the algorithm proposed in this paper.

Conclusions
Microgrippers are the end effectors in micro-operations and micro-manipulations, making them a sensitive instrument in terms of accuracy and efficiency. As an actuation mechanism for microgrippers, piezo-electric actuation has several benefits, but there is a major drawback; the range of motion in this type of actuation is minimal. Therefore, amplification mechanisms are of great interest, especially in piezo-electric microgrippers. Amplification Factor is an important parameter to ensure a great range of motion for the gripper, thus resulting in a wider application use. Using four stages of amplification, we achieved a desirable range of motion for the microgripper. This paper proposed a novel optimization process to adjust a symmetric compliant piezoelectric actuated microgripper for specific applications. A trade-off was made between the displacement amplification factor and maximum generated mechanical stress using multi-objective optimization. Based on the optimization, a set of designs was proposed instead of a single optimum design. A selection method based on the Analytical Hierarchy Process (AHP) was suggested for the selection process of designs for specific applications. The performance of the proposed microgripper design was inspected using the finite element method. Taguchi's method of designing experiments was used to identify effective variables and obtain relations between design parameters and intended responses. Using NSGA-II instead of the traditional Multi-Objective Genetic Algorithm (MOGA) reduces the complexity of calculations, speeds up the algorithm, leads to easier elitism, and preserves the GA population's diversity, leading to more accurate and faster results. Also, using the AHP algorithm prevents possible errors in the final decision of the designer.  Author contribution All of the authors of this manuscript confirm that this work is original and has not been published elsewhere nor is it currently under consideration for publication elsewhere. They confirm that the manuscript has been read and approved by all named authors and that there are no other persons who satisfied the criteria for authorship but are not listed. The authors further confirm that the order of authors listed in the manuscript has been approved by all of them.
Funding They state that no funds, grants or other financial support were received for conducting this study and have no relevant financial or non-financial interests to disclose.

Conflict of interest
The Authors declare that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.