Abstract

Flexure-based compliant mechanisms are increasingly promising in precision engineering, robotics, and other applications, thanks to the excellent advantages of no friction, no backlash, no wear, and minimal assembly. However, their design and analysis are still challenging due to the coupling of kinematic-mechanical behaviors with large nonlinear deflections in comparison to their rigid-body counterparts. Optimal design is an important aspect in the field of compliant mechanisms and has attracted much attention during the last decades. Especially, when considering a multiobjective optimization design for compliant mechanisms, the problem is becoming more complicated. Thus, this article presents a new efficient hybrid computational method to resolve multiobjective optimization design of compliant mechanisms. A Scott Russell compliant mechanism is employed as the design example and to show the advantages of the proposed optimizing method. The proposed method is developed by hybridizing the desirability function approach, fuzzy logic system, adaptive neuro-fuzzy inference system (ANFIS), and lightning attachment procedure optimization (LAPO). First of all, a 3D finite element model is created and central composite design is employed to build a numerical matrix. First, design variables are refined by using analysis of variance and Taguchi approach in terms of considerably eliminating space of design variables. Subsequently, desirability values of two objective functions are determined and transferred into the fuzzy logic system. The output of the fuzzy logic system is considered as a single combined objective function. Next, modeling for fuzzy output is established via developing the ANFIS model. At last, the LAPO algorithm is adopted for solving the multiobjective optimization problem for the mechanism. Three numerical examples are investigated to validate the feasibility and the effectiveness of performance efficiency of the proposed methodology. The results find that the proposed methodology is more efficient than traditional Taguchi-based fuzzy logic. Besides, the performance efficiency of the proposed approach outperforms the Jaya algorithm and TLBO algorithm through the Wilcoxon signed rank test and Friedman test. The results of this article give a useful approach for complex optimization problems.

1. Introduction

Scott Russell compliant mechanism (SRCM) is utilized as a shuttle to transform two input motions into output one. It is exerted by a piezoelectric actuator in applications of high precision positioning systems. The mechanism can achieve its performances through elastic deformation of flexure hinges [14]. It benefits a monolithic structure, lightweight, and free friction. Nonlinear behaviors of the mechanism are facing challenges for researchers in compliant mechanism. Over the past several decades, researchers have tried to make approaches and methods to synthesize compliant mechanisms, e.g., kinematic-based methods [58]. They showed that physical behaviors and performances of compliant mechanism should analyze displacement, stress, frequency, and parasitic error. Kinematic-based methods are still valid, but some are difficult to resolve the performances of SRCM because of their complicated applications.

In the past decades, a lot of analytical methods and numerical methods have been proposed for designing and analyzing compliant mechanisms, such as the pseudorigid-body model (PRBM) [9], Castigliano’s second theorem [10], compliance matrix method [11], empirical method [12], beam constraint model [13], Euler–Bernoulli beam theory [14], and finite element method [15]. However, these methods are useful in advancing simple compliant mechanisms with existing configurations. But they are critically challenging in analyzing the static and dynamic behaviors for complex compliant mechanisms, such as irregular shapes of flexure hinges and complex configurations and large deflections. It is known that the PRBM is still a valuable technique for modeling, but it is limited in highly nonlinear deformation of flexure hinges. For complex loads and complicated configurations of flexure hinges, the PRBM is not suitable. Especially, the PRBM is strongly dependent on the location of torsional springs, the number of torsional springs, and designer’s experiences in the field of compliant mechanisms. Castigliano’s second theorem can analyze the strain energies of flexure hinges such as tensile, shear, and bending strains, but this theorem is critically complicated for complex structures. Compliance matrix method is actually difficult for analyzing when compliant mechanisms subject multiple complicated actuation forces. Elastic beam theory is related to the deformation energy of compliant mechanisms, but it is still complicated for complex structures. Based on the two-port dynamic stiffness model, statics and dynamics can be effectively solved, but it becomes more complicated for modeling irregular structures. Compliance matrix method is used to transfer the compliance of flexure hinges from the local coordinate frame into a reference coordinate system. This method is still complicated for a more complex structure of compliant mechanisms. Empirical modeling can analyze compliant mechanisms with high accuracy, but it is very costly and time-consuming. It is difficult to analyze a large deformation of flexure hinges using the beam constraint model. Finite element method (FEM) analyzes complex geometric shapes of compliant mechanisms by discretizing compliant mechanisms into elements and nodes.

Although the previous methods are still valid, they are still challenging to analyze the complex compliant mechanisms. The estimating accuracy of modeling are based on a complexity degree of kinematic and mechanical coupling of compliant mechanisms. Especially for a large nonlinear deflections and irregular shapes, the abovementioned methods are limited. For optimal design of multiple performances of compliant mechanisms, the problem is becoming more complex. Therefore, data-driven approach, computing intelligence, and machine learning have attracted much attention to model complex systems [1618]. Data-driven approaches-coupled computing intelligence is an effective tool that has not extended to the field of compliant mechanisms yet. From reviewing the literature review, the motivation of this article is to develop a new hybrid computational method to solve the optimal design for compliant mechanisms. The proposed method is extended to solve for most compliant mechanisms from a basic structure to a complex one. Specifically, it is an efficient approach for the SRCM. In order meet wide applications, the SRCM requires both a large displacement, a high frequency, and a good enough working strength. In order to solve three quality performances simultaneously, size, shape, or topology of the proposed mechanism can be changed. However, a such random change spends a lot of computational time. At the same time, a parameter optimization is conducted by researchers to enhance the physical performances. Optimization problem includes a single response or more than two ones. In the past, the structural optimization problem of a compliant mechanism separates into three main types, e.g., topology optimization [19], shape optimization, and size optimization [14, 20, 21]. In this work, size optimization is chosen to find the best geometrical parameters of SRCM.

In order to solve multiple objective responses, a multiobjective optimization (MOO) problem is investigated. The aim of the MOO problem is to reach a trade-off between characteristics. Nowadays, the MOO problem has received a great interest for compliant mechanisms [2224]. However, a MOO optimal design for the proposed SRCM has not been studied in depth. So, parameter optimization is proposed in the present article. In a real structure of SRCM, it includes possible design parameters affecting desired performances. Besides, it is well-known that physical behaviors of the SRCM are very sensitive to a change in shape, material, or geometrical parameters. Shape and material may be directly varied according to customer’s demands. Meanwhile, structural parameters should be optimized to improve the performances. The most critical parameters strongly contribute on the performances, while some nonsignificant parameters lead high computational cost. Hence, before implementing the MOO process, this study determines the most significant design variables and nonsignificant parameters.

Basically, before dealing with a MOO problem for a specific compliant mechanism, a mathematical model is expected to be established through kinematic-based analytical approaches such as pseudorigid-body model [25], semianalytical model [26], and compliance matrix [27]. Although mentioned approaches are still used, they are complicated for modeling the SRCM. In other words, analytical methods are limited to feature its behaviors of more complex structures. In order to avoid an imprecise optimization solution since approximate approaches may be complicated, data-driven approaches are more advantageous to adopt for the MOO problem. Data-driven approach is utilized to directly map the inputs and outputs data; therefore, this approach has a good predicting accuracy in comparison with the analytical methods. During the last decades, a lot of methods have been developed to solve the MOO problems, such as desirability function approach (DFA) [28], grey relational analysis (GRA) [29], Taguchi method-based fuzzy logic (TMFL) [30]. Both DFA and GRA need a weight factor for each objective function, but weighting factor is largely dependent on experience or users. Meanwhile, TMFL is a statistical-based approach without the weight factor. Although TMFL has more benefits, it searches optimal values at discrete points from the Taguchi method. In order to avoid a local result, surrogate models are employed to connect inputs and outputs, e.g., response surface method (RSM) [31], Kriging technique [32], artificial neural network (ANN) [3335], adaptive neuro-fuzzy inference system (ANFIS) [3638], RSM-ANN [39], ANN machine learning [40], and fuzzy logic theory [41]. Among them, ANFIS is proved as an effective technique to modeling and predicting, which is applied for different engineering fields. For example, ANFIS was developed to predict for various fields, e.g., the hydropower generation [17], the scour hole shape [18], and mercury emissions in combustion flue gases [42]. However, ANFIS has not extended for modeling the Scott Russell compliant mechanism yet.

Up to now, nature-based metaheuristic algorithms have been developing very fast, e.g., genetic algorithm [43], particle swarm optimization (PSO) [44], differential evolution [45], cuckoo search algorithm [46], and other algorithms [4750]. Some of them have successfully been applied in various fields, such as grey wolf optimization for temperature-based oil-water relative permeability [16] and PSO algorithm for predicting the scour depth around bridge piers [51]. More recently, several algorithms with less or nontuned parameters have been suggested such as teaching learning-based algorithm (TLBO) [52, 53], Jaya algorithm [54], and lightning attachment procedure optimization (LAPO) algorithm [55]. Among three less-parameter algorithms, LAPO algorithm is an effective approach for a lot of engineering areas, but it has not been applied for the SRCM mechanism yet. Therefore, this article chooses the LAPO algorithm to extend the MOO design. In summary, highlight contributions of this study are as follows: (i) limit nonsignificant parameters and refine space of design variables; (ii) transfer the desirability value of objective functions into the fuzzy logic reasoning; (iii) develop the ANFIS structure to model the output of FIS system; (iv) LAPO algorithm is proposed for implementing the MOO problem to find global solution; (v) performances of proposed methodology are compared with other algorithms using statistical techniques.

The goal of this article is to develop new hybrid computational methodology to conduct MOO optimal design for the SRCM mechanism. As mentioned above, the proposed method is considered as statistical-based intelligent computation. A 3D model is created by the finite element method, and central composite design is employed to make design of experiments. Initial space of design variables is determined and then refined via analysis of variance and Taguchi technique. New spaces of design parameters are generated through sensitivity investigation. Next, desirability values of quality responses are calculated and transferred into the fuzzy logic system. ANFIS technique is developed to connect design variables and output of the FIS. MOO problem is then solved by using the LAPO algorithm. Numerical examples are investigated to evaluate the computational efficiency and performance of proposed method. In addition, performances of the proposed approach are compared with the ANFIS-based Jaya and ANFIS-based TLBO by implementing the Wilcoxon signed rank test and Friedman test.

2. Proposed Hybrid Methodology

Figure 1 describes a new computational method which is hybridized by the integrating desirability approach, fuzzy logic reasoning, ANFIS, and LAPO algorithm. The proposed method is applied to solve a MOO design for a Scott Russell compliant mechanism. Practically, this mechanism is combined with a piezoelectric actuator (PA) for use in ultrahigh precision positioning devices. Because the PA limited its travel and responding speed, the mechanism needs to improve its displacement and a high natural frequency. In addition, a small equivalent stress is considered as an additional constraint to ensure a good working strength. According to background of compliant mechanism, these specifications are very sensitive to geometrical parameters of flexure hinges. In other words, a change in any geometrical factor may lead to a change in performances. In order to achieve an improvement in performances, design parameters are optimized simultaneously. Proposed hybrid methodology consists of five main phases as follows: (i) design and analysis, (ii) desirability calculation, (iii) fuzzy logic system analysis, (iv) ANFIS modeling, and (v) conduct MOO problem by LAPO algorithm.

2.1. Phase 1: Design and Analysis

This phase concerns about a mechanical design and analysis problem. In this study, a MOO design for a SRCM is investigated. It is an example to show effective application to validate the performance effectiveness of proposed methodology. This phase experiences the following main steps.

Step 1. Mechanical structure
As discussed above, the mechanism needs a large displacement and a high frequency but must guarantee a small stress for a long operation without plastic failure.

Step 2. Define initial design variables
Structure of the mechanism has a lot of possible geometrical parameters affecting the quality responses. Such influences must be considered before implementing further the optimization process. Some factors with actual influences as well as no influences are evaluated through sensitivity investigation in this study.

Step 3. Define quality responses
After determining the initial design variables, a large displacement and a high frequency are assigned as two desired outputs. In order to satisfy in terms of enough working strength, equivalent stress is an extra constraint.

Step 4. Create 3D finite element model
A 3D finite element model (FEM) is created, and then finite element analysis (FEA) is implemented to retrieve numerical data. To complete this work, design variables and quality responses are assigned to be parametric factors during FEA simulations.

Step 5. Central composite design
While a 3D model is created, an experimental matrix is also initialized by using central composite design (CCD). Face centered is used in this work.

Step 6. Run numerical data
Boundary conditions, load, and material are set up for the 3D model. Subsequently, a set of numerical data is collected by running simulations.

Step 7. Sensitivity analysis
This is a critically important step. It determines several parameters directly affecting the quality characteristics to be design variables. Besides, some parameters with no influence are ignored. This process is conducted through analysis of variance (ANOVA) and the Taguchi method.

Step 8. Eliminate design variables
The results of sensitivity analysis in Step 7 identify some nonsignificant parameters. And then, such parameters are ignored or eliminated while other important factors are defined as key design variables. The results of refining parameters perhaps create several new population spaces for population-based optimization metaheuristic algorithms. This is a preparing step for further optimization process. The goal of Step 8 is to determine a suitable searching space so as to reduce the computational cost and make more reliable solution.

Step 9. Repeat central composite design
Several new spaces of design variables are created in Step 8. Each space of refined variables is considered as a case study. And then, Step 9 will rebuild a number of experimental matrixes by using CCD.

Step 10. Refine 3D finite element model
In order to conduct the numerical simulations, a 3D FEM model is created with respect to each case study. Subsequently, simulations are implemented to get the numerical data for each case study.

2.2. Phase 2: Desirability Calculation

Step 11. Recollect numerical data
Numerical data for each of case study are collected through FEA simulations, individually.

Step 12. Compute desirability value for performances
According to theory of the desirability function approach [56], the exponential type is transformed to each response into a desirability function. Principle of this technique is to compute the ith quality response () to become the ith desirability (di), and then all values of desirability are combined into a single desirability function (D). Desirability value can reach unity when the response’s target is achieved, and otherwise. Each individual desirability function is divided into three types based on user’s specific demand. It includes three target types, i.e., smaller-the-best, higher-the-best, and normal-the-best. Desirability for two specifications of the mechanism is computed. Displacement is millimeter, while natural frequency is in Hertz. Based on calculating the desirability, influence and difference of units are neglected during the optimization process. So, such difference no longer affects later optimal solution.

2.2.1. Smaller-the-Best

The response’s value is expected lower than an upper limit. The desirability function is identified bywhere di is the desirability value, is the ith response, and Lb and Ub are the lower and upper limits of the response.

2.2.2. Normal-the-Best

The response’s value is expected toward a target value (T). The desirability function can be determined bywhere and q are specific parameters defined by users (, q > 0) which determine the shape of di.

2.2.3. Higher-the-Best

The response’s value is expected higher than a lower limit. The desirability function can be described bywhere r is the desirability function index.

All individual desirability’s values of responses are transformed into overall desirability, D, which is considered as a single quality index. This index is determined by assigning the corresponding weight factor for different quality responses. However, weight factor of each response is dependent on priority or customer’s demands. A combined desirability index is computed aswhere D is the overall desirability index and is the weight of the ith response. D is equal to one as each di is also equal to one. Otherwise, at least, one of di is zero and D is zero.

In this paper, a larger type is used for both desired targets (displacement and frequency). Previously, in order to solve a MOO problem, a set of optimal parameters may be found by maximizing the single quality index D. Although this approach is still an effective tool, an optimal result is varied when the weight factor of each response is changed. It may result in an imprecisely optimal value. In order to overcome this uncertainty, a fuzzy logic system is then developed to deal with all individual values of desirability of both responses since this system does not require any weight factor.

2.3. Phase 3: Fuzzy Logic Reasoning System

A fuzzy logic reasoning system includes knowledge base, fuzzifier, inference engine, and defuzzifier. Details can be briefly described as follows [57].

2.3.1. Fuzzifier

Fuzzifier plays a role to transform real value into a fuzzy system. Inputs of the fuzzy logic system are called as crisp values which contain precise information of a real world. Through the fuzzifier, real value is transformed into linguistic variable.Knowledge base: knowledge base consists of rule base that forms a number of fuzzy rules (if-then). Moreover, the knowledge base contains a database which defines membership function (MF) of the fuzzy sets.Inference engine system: based on fuzzy rules, fuzzy inference system (FIS) is considered as decision-making. It handles how the rules are combined.

2.3.2. Defuzzifier

It transfers the output of the FIS system into a crisp value. In the defuzzification method, the centroid method is used to the transformation. The output of the FIS system, a nonfuzzy value, is called as a multicharacteristic performance index (MPCI).

In order to implement the FIS system, the Mamdani method is employed in the present paper. Subsequently, trapezoidal MFs are adopted for the inputs and outputs of the FIS so as to form fuzzy sets. MFs are in the range from zero to one, and MFs can describe the way a variable matches a fuzzy set. Inputs and outputs of fuzzification system are then transformed into linguistic variables. The trapezoidal MFs are defined aswhere μA denotes the MFs while k, l, m, and are parameters and x is a variable.

In this study, desirability of frequency and desirability of displacement are calculated, and then they are considered as two input variables for the FIS. These linguistic input variables will be combined into an output for the FIS. The trapezoidal MFs are employed for both the fuzzification and defuzzification. The following fuzzy rules are briefly described.

Fuzzy rule: if x1i is A1 and x2i is B1i, then yi is Ci else (i = 1, 2, 3, …, n), where x1i and x2i are the two ith inputs and yi is the ith output. Ai, Bi, and Ci are defined by corresponding MFs (μAi, μBi, and μCi), and these parameters are regarded as fuzzy subsets.

In order to compute the fuzzy logic reasoning, max-min operation of Mamdani is adopted. Subsequently, the FIS output is retrieved. The MFs of the FIS output can be described by

At last, the FIS output is transformed into real values through the defuzzification. Subsequently, a nonfuzzy value y0, which is called as MPCI, is defined by

Based on the theory of FIS system, the best solution for overall responses and a set of optimal design variables can be found by maximizing the MPCI index through the Taguchi method. This technique is divided into three types such as (i) higher-the-better, (ii) normal-the-best, and (iii) smaller-the-better. The larger-the-better type is chosen for maximizing the MPCI, which is described aswhere MPCIi is the ith MPCI index and n is the number ith experiment’s repetition.

A Taguchi-based fuzzy logic can find optimal candidates in a MOO problem, but this approach may be a local optimum solution. The reason is because the Taguchi technique is employed to minimize or maximize a single fitness function in terms of discrete values. Meanwhile, a real problem desires to search a global optimum solution. In order to overcome such situation, ANFIS is then extended to model the MPCI, and the MOO design for the SRCM mechanism can be effectively solved by using the LAPO algorithm.

2.4. Phase 4: Adaptive Neuro-Fuzzy Inference System

As aforementioned, the FIS system is a modeling technique in terms of linguistic variables where the Mamdani method is used for the FIS. Meanwhile, ANFIS is an artificial technique which is developed by integrating neural network and the FIS. Nowadays, ANFIS is considered as an intelligent modeling which creates a connection between inputs and outputs. In theory of ANFIS, the Sugeno model is employed to create fuzzy rules [58]. The fuzzy rules for the ANFIS model can be defined aswhere x1 and x2 are inputs with respect to A1 and A2 term sets, y is the output, and a, b, and c are constant values.

Figure 2 illustrates the ANFIS structure which consists of five-layer feed-forward neural network.

Layer 1 is the fuzzification layer that assigns the membership degrees for input factors based on the given MFs. The output of layer 1 can be described aswhere x is the input with respect to node ith and Ai is a linguistic label. is the MF of Ai.

Layer 2 includes the fuzzy rule and a rule node which gets inputs and identifies firing strength of the rule. Each node is labeled as a circle node, П. Each node output can be defined as

Layer 3 is a normalized layer which is used to estimate a ratio of firing strength of a given rule to total of firing strengths of all rules. Every node is labeled as a circle node, N. In this layer, is named for normalized firing strength of rules and defined as

Layer 4 is the defuzzification process, and the ith node is labeled as square node by

Layer 5 is an overall output and sum of all signals, which is defined as

In this article, the trapezoidal MFs are adopted for the ANFIS system as well.

2.5. Phase 5: Lightning Attachment Procedure Optimization Algorithm

LAPO is a new algorithm based on a physically lightning phenomena [59]. It is one of the most nontuned parameter algorithms. LAPO algorithm can solve a lot of different engineering problems, and it is superior to other metaheuristic algorithms [60]. Therefore, this study adopts the LAPO algorithm to expand its application capacity for the MOO design for the bridge amplification compliant mechanism. Readers can find details in literatures [59, 60]. It consists of five substeps as follows.

Step 1. Initialization
The algorithm needs an initial population in a search space of design variables. Unlike previous studies, the present work evaluates and refines the search spaces in which a set of decision variables are defined based on statistical analysis. In the LAPO, some test points are assigned in the cloud and ground surface, which can be defined aswhere and denote the upper and lower bounds of design variables. Rank is a randomly given parameter in range from zero to one. Fitness of the solution is defined as

Step 2. Next jump determination
The average value of the overall test points and the value of fitness of the points are computed as in equations (17) and (18), respectively:where and denote the average value and value of fitness of the points.
If the electric field, fitness function, of potential point h is greater than the average value of the test points, it is calculated asOtherwise, it is calculated as

Step 3. Branch fading
If the fitness function of new test point is better than the previous point, the branch sustains; otherwise, the branch fades asThis step is carried out for overall the test points so as to rest points moving down.

Step 4. Upward leader movement
Upward leader movement calculates an exponent factor based on the charge of downward leader aswhere noi and noimax denote a number of iterations and maximum number of iterations, respectively. The upward leaders are defined for the next strategy of test point aswhere and denote the best value and the worst solution of the population.

Step 5. Final jump
The process is ended when the upward leader meets the downward leader to each other. It is noted that the optimization process is ended if stop criterion is achieved. In addition, the maximum iteration is met, and the process is ended herein.

3. Case Studies: MOO Design for Scott Russell Compliant Mechanism

MOO design of the Scott Russell compliant mechanism is studied. This article considers three numerical examples, and the obtained results are to evaluate application’s capability and effectiveness of the proposed hybrid computational method.

3.1. Structural Description and Problem Statement

Figure 3 describes a structural design of the Scott Russell compliant mechanism (SRCM). Bottom end of the mechanism is fixed by screws, and top end is output. The mechanism is a type of bridge amplifier, but four bridge legs are alternated by using four Scott Russell mechanisms. The purpose of alternative is to enhance stiffness of traditional bridge mechanism. Each Scott Russell mechanism includes four flexure hinges and rigid links, as shown in Figure 3. The output, end effector, is at the middle rigid link. Based on symmetry, the end effector can move along the y-axis by exciting two input loads from a piezoelectric actuator. Based on this design, the mechanism can amplify a large working stroke along the y-direction, but its natural frequency is decreased. Another problem is how to ensure an enough working strength, i.e., a long fatigue working life without plastic failure. In a such situation, this study proposes a new optimal design methodology to resolve a balance between displacement and frequency. The main parameters of SRCM are illustrated, as given in Figures 3 and 4. In this article, AL T73-7075 material is chosen for the mechanism. In this structure, some parameters are design variables and others are constant values.

Table 1 gives geometrical factors and material properties of the proposed mechanism. According to the designer’s experiences and expert’s knowledge in the field, spaces of design variables are determined based on a capacity of fabrication device.

The purpose of this study is to develop a new hybrid computational approach to deal with the MOO problem for compliant mechanisms. The SRCM is an application example. Numerical investigations are carried out to validate the proposed hybrid computational method. To address about the MOO problem for the SRCM, the design variables, objective functions, and constraints are firstly identified. According to the theory and investigations on compliant mechanisms [10, 6164], geometry and locations of flexure hinges mainly contribute on the performances. As seen in Figures 3 and 4, the dimensional parameters are represented as a vector X = [θ, L0, L1, L2, L3, L4, L5, L6, D, T, R, H, L, W]T in which the key design variables are affecting the performances, including θ, L1, L2, L3, L4, and D. Regarding objective functions and constraints, the SRCM is expected to achieve a fast speed, a large working space, and a good strength. In order to reach a desired speed, the first natural frequency of the SRCM must be improved accordingly. Meanwhile, a large displacement of the SRCM can help toward a large working stroke. At last, a good strength can be gained if the resulting equivalent stress of SRCM is under the yield strength of material [10, 6164]. In addition, based on the experiences in designing and analyzing kinematics of compliant mechanisms, the mentioned parameters (θ, L1, L2, L3, L4, and D) are very sensitive to the performances. They can be chosen properly to improve the performances. The formulation of the optimization problem of the SRCM is briefly stated as follows.

Find vector of design variable: X = [θ, L1, L2, L3, L4, D]T

Subject to constraints:

Initial space of design variables:where F1(X), F2(X), and F3(X) denote the frequency, displacement, and equivalent stress, respectively. is the yield stress of proposed Al T73-7075. Specifically, F1(X) and F2(X) are two objective functions, while F3(X) is regarded as the constraint. SF is the safety factor that is chosen as 1.5.

3.2. Numerical Simulation

A 3D FEM model is created, and numerical experiments are formulated via the CCD. Numerical data are then collected by implementing FEA simulations. The boundary conditions and loads are given from Figure 3. A force, F, of 70 N is applied toward an end of inputs, and the output response is measured. Automatic technique is utilized for meshing. Meshing of flexure hinges is refined again in order to achieve accurate analysis, as shown in Figure 5. The meshing results give a number of elements of 1544 and a number of nodes of 9897. Skewness criteria, metric performance, are used to evaluate the quality of meshing. Its value is in the range from zero to one. The results found that the average value of skewness is approximately 0.3745, as given in Figure 6. On the contrary, this value guarantees the quality of meshing is a relatively good result.

4. Results and Discussion

4.1. Sensitivity Investigation and Refinement of Design Variables

The CCD is utilized to build the numerical matrix by using space of all initial design variables. And then, the most critical important parameters and suppress nonsignificant parameters are identified. Secondly, this analysis refines space of design variables prior to implement the MOO problem of proposed mechanism. In this article, six factors are divided into three levels, as given in Table 2. Each numerical experiment retrieves two quality specifications and one constraint, including natural frequency (F1(X)—Hz), displacement (F2(X)—mm), and equivalent stress (F3(X)—MPa). Forty-five experiments are generated, and the numerical results are given, as in Table 3.

As mentioned above, this article considers displacement, frequency, and stress for MOO design of the Scott Russell compliant mechanism. Geometrical parameters directly affecting those specifications are identified and may be suppressed for further optimization process. Sensitivity investigation of design variables with respect to the three specifications is analyzed, individually. In this study, analysis of variance and Taguchi approach are used and implemented by using Minitab software 18.

Case 1, considering the displacement, ANOVA is utilized to analysis of sensitivity at 95% confidence interval. The result of ANOVA in Table 4 found that the contribution of parameters θ (the-ta), L1, and L2 on this quality response is about 0.03%, 0.02%, and 0.01%, respectively. These contributions are too small compared to remaining parameters. It means that the contributions of θ, L1, and L2 can be ignored in further modeling and optimization process. In addition, the Taguchi method with the larger-the-better type is used for the displacement, as shown in Figure 7. The results also found that factors θ, L1, and L2 are also two nonsignificant factors because the displacement is a little change when those parameters are varied. Meanwhile, other parameters are very sensitive to the displacement which are taken as key variables for case study 1.

Considering the frequency, the ANOVA results in Table 5 reveal that parameters L3 and L4 have a lowest contribution on the displacement with respect to 0.09% and 0.22%, respectively. As a result, these two factors are nonsignificant variables that can be neglected in modeling and optimization process for case study 2.

Figure 8 describes the sensitivity of full variables on the frequency via using the Taguchi method with the larger-the-better type. It also notes that the parameters L3 and L4 have a very small influence on this response compared to other parameters. On the contrary, if these parameters are varied, the frequency is somehow changed. Meanwhile, a change in the other parameters has a sharp change.

Case 3 considers the equivalent stress, and the ANOVA results in Table 6 indicate that three parameters θ, L2, and L4 have the lowest contribution on the stress with 0% contribution. These three factors are nonsignificant and can be suppressed in further modeling and the MOO process. A smaller-the-better type is utilized for the stress. The results of the Taguchi analysis show that three parameters also have a very small contribution on the stress response, as shown in Figure 9.

4.2. Optimization Implementation
4.2.1. Selection of Membership Functions

Each real value of physical response is transformed into linguistic variables which are assigned for the fuzzy reasoning inference process. The linguistic variables are determined based on expert’s knowledge and designer’s institution. Table 7 gives the linguistic variables, including tiny, very small, small, small-medium, medium, medium-large, large, very large, and huge.

In this article, two inputs of FIS structure consist of desirability of frequency and desirability of displacement. The MFs for desirability of frequency are divided into three levels such as small, medium, and large, as illustrated in Figure 10(a). Meanwhile, the MFs of the desirability of displacement are separated into seven levels, including tiny, very small, small-medium, medium, medium-large, large, very large, and huge, as depicted in Figure 10(b). At last, the MFs of the output MPCI also include seven levels, as shown in Figure 11. All the MFs are defined by trapezoidal shape because this is a popular MF type for both the FIS and ANFIS models.

The results of sensitivity analysis are used to decrease space of design variables. Through this analysis, computational time is lowered and reliable solution is increased. Particularly, some parameters may be redundant, while others are critically important factors. Based on the sensitivity results, initial design variables are eliminated. From that, three following numerical examples are considered to evaluate the performance effectiveness of proposed computational methodology.

If full initial space is used, it may take a long time to search for optimum solutions. Therefore, this study limits design space and following case studies are separated to investigate and evaluate the efficiency of the proposed hybrid approach.

4.2.2. Case Study 1

By using the results of sensitivity analysis in Table 4 and Figure 7, some nonsignificant factors are suppressed since their influences are very small on the displacement. Real influencing factors with large contributions are taken in account as key design variables. Case study 1 takes two objective functions to be optimized simultaneously with the spaces of design variables being relatively reduced. In order to solve the MOO problem, the desirability of displacement and the desirability of frequency are transformed into a single combined objective function, which is called as the MPCI of the FIS system. Subsequently, the MPCI is maximized by using the proposed hybrid approach. The optimization problem of case study 1 is stated as follows.

Find vector of design variable: X = [D, L3, L4]T:s.t.

Case study 1 evaluates three design variables (D, L2, and L1) and two objective functions and one constraint. Such design variables and both responses are set as parametric variables in the FE model. Based on the number of design variables and their levels, CCD is used again to build 15 experiments. The results of frequency, displacement, and equivalent stress are collected, simultaneously.

Next, desirability of frequency and desirability of displacement are calculated, as given in Table 8. Real values of two responses have different units, and it is a main cause resulting in unprecise solutions. Therefore, desirability function approach is proposed to overcome this limitation because desirability has no unit and its value is from zero to one. The results found all stress values are satisfied in terms of allowable stress of material.

In order to conduct the MOO problem for case study 1, the MFs for desirability of frequency include three levels and the desirability of displacement is divided into seven levels, while the MFs of fuzzy MPCI output are separated into seven levels, as given in Table 9. Subsequently, matrix of fuzzy rules is built in this table.

In this study, center of gravity method is employed for defuzzification and Mamdani method is utilized for fuzzy operation. The FIS is implemented in Matlab R2015a. Figure 12 illustrates a relation of the fuzzy MPCI with the desirability of displacement and desirability of frequency. The value of inputs and outputs of the FIS is in range from zero to one. In this study, the weight of each rule is assumed to be unity.

Based on 27 fuzzy rules in Table 9, the output of FIS system is computed, as shown in Figure 13.

The results of MPCI output are determined with respect to each value of desirability of frequency and desirability of displacement, as given in Table 10.

By taking the values of design variables in Table 8 and the MPCI value in Table 10, a proper ANFIS model is developed for modeling the MPCI. Modeling is coded in Matlab. The results of ANFIS parameters are given in Table 11. In this study, 75% of data and the remaining 25% of data are utilized for the training and testing, respectively.

In order to evaluate the built ANFIS model, some performances metrics are employed, such as mean absolute percentage error (MAPE), root mean square error (RMSE), and correlation coefficient (R2) [18, 65]. Mathematical equations of MAPE, RMSE, and R2 are described bywhere yi is the ith of the measured result, is the ith of the predicted result, is the mean value of measured results, and n is the number of observations.

As the discussion in Section 3.1, a high frequency, a large displacement, and a minimal stress are desired simultaneously. According to the fuzzy logic theory [30], the largest value of the MPCI is expected to achieve the optimal solution. In other words, the MPCI must be maximized so as to achieve a global optimal solution for the SRCM. The results found that the MAPE and RMSE are close to zero, while R2 reaches one. It means that the built ANFIS model is relatively good. In order to evaluate the effectiveness and robustness of the proposed hybrid approach, the optimal results are compared with those obtained from a traditional combination of the Taguchi method and fuzzy logic (TMFL), and the optimal parameters from the TMFL are determined at D = 2.25 mm, L3 = 8.5 mm, and L4 = 18.7 mm, as given in Table 12. The optimal results of TMFL find that the frequency and displacement are approximately 24.4574 Hz and 0.1691 mm with a small equivalent stress of 251 MPa, respectively. It also ensures a high safety factor. However, these optimal values are found at discrete points according to theory of the Taguchi approach. Such points can lead to local optimum solutions. To overcome this limitation, LAPO algorithm is used to search a global optimum solution. Based on the established ANFIS modeling, a pseudoobjective function of MPCI is well established. And then, a LAPO programming is implemented in Matlab R2019b, and the results of the proposed approach found that the optimal design parameters are at D = 2.30 mm, L3 = 8.5 mm, and L4 = 18.7 mm. By using proposed methodology, the optimal solutions find the frequency and displacement are about 26.592 Hz and 0.1539 mm, respectively. The equivalent stress is also lower than the yield strength of Al material. It permits a good enough working condition with a safety factor of 1.98. Moreover, the MPCI predicted from the TMFL is smaller than that predicted from the proposed hybrid methodology. According to a larger-the-better value is best for the MPCI, the proposed approach is efficient, robust, and better than the TMFL. The optimal results are satisfied with initial requirements.

4.2.3. Case Study 2

Based on the results of sensitivity analysis in Table 5 and Figure 8, some unimportant factors are neglected in modeling and optimization problem. The factors with actual contributions are considered as key design variables. The space of design parameters is decreased, accordingly. Similarly, case study 2 also maximizes F1(X) and maximizes F2(X), simultaneously. Two objective functions are transferred into desirability values and then transformed to the MPCI value. The optimization problem is briefly expressed as below.

Find design variable X = [θ, L1, L2, D]T:s.t.

Case study 2 considers two design variables and two objective functions and three constraints. CCD is also used to build 25 numerical experiments. The results of frequency, displacement, and equivalent stress are collected by simulations. And then, desirability of displacement and desirability of frequency are calculated, as given in Table 13. The resulting stress values are under the allowable stress of Al T73-7075.

Next, the MFs for desirability of frequency also are divided into three levels, and the desirability of displacement is seven levels, while the MFs of fuzzy MPCI output are separated into seven levels, as given in Table 9. Subsequently, matrix of fuzzy rules is built in Table 14.

The results of MPCI output are calculated corresponding to the desirability two responses, as given in Table 15.

By taking the values of design variables in Table 13 and the MPCI value of Table 15, ANFIS model is formed for modeling the MPCI. The ANFIS parameters for case study 2 are developed, as given in Table 16. The predicting accuracy of the ANFIS model is well evaluated through equations (29) and (31) before conducting the MOO problem for the SRCM.

By using the TMFL, the optimal parameters are determined at θ = 16.5°, L1 = 7.7 mm, L2 = 18.7 mm, and D = 2.25 mm, as given in Table 17. Based on the TMFL, the optimal frequency and displacement are found at 28.94823 Hz and 0.1527 mm, respectively. Subsequently, LAPO algorithm is used to search a global optimum value. The results of the proposed hybrid approach determined that the optimal design parameters are θ = 16.5°, L1 = 8.0 mm, L2 = 14.5 mm, and D = 2.25 mm. Using the proposed methodology, the optimal results found that the frequency and displacement are about 651 Hz and 0.1981 mm, respectively. Finally, the predicted MPCI from the proposed approach is greater than that predicted from the TMFL. It means that the estimated solution from the proposed methodology is robust and better than that from the TMFL. Moreover, it finds that the equivalent stress is under the allowable stress of Al T73-7075. The optimal results are satisfied with design constraints.

4.2.4. Case Study 3

Based on the results of sensitivity analysis in Table 6 and Figure 9, the space of design variables is relatively shortened. Case study 3 also simultaneously maximizes F1(X) and maximizes F2(X). The computational principle is similar to previous cases. The statement of optimization is stated as follows.

Find design variable X = [L1, L3, D]T:s.t.

Case study 3 regards three design variables and two objective functions and one constraint. CCD is used to build 15 numerical experiments. The results of frequency, displacement, and equivalent stress are collected. And then, desirability of two responses is calculated, as given in Table 18. The results found all stress values are smaller than the allowable stress of the proposed material.

The MFs for desirability of frequency include three levels, and the desirability of displacement is consisted of seven levels, while the MFs of fuzzy MPCI output are divided into seven levels. Subsequently, a matrix of fuzzy rules is built in Table 19. The value of MPCI is determined according to each value of desirability for objective functions, as given in Table 20.

Subsequently, using the values of design variables in Table 18 and the MPCI value of Table 20, the ANFIS model is formed for mapping the design variables and the MPCI. The ANFIS parameters for case study 2 are given in Table 21. The predicting accuracy of the ANFIS model is well confirmed through equations (29) and (31) prior to implement the MOO problem for the SRCM.

By using the TMFL, the optimal parameters are determined at L1 = 7.7 mm, L3 = 10.1625 mm, and D = 2.25 mm, as shown in Table 22. Based on the TMFL, the optimal frequency and displacement are found at 158.970 Hz and 0.0957 mm, respectively. Such optimal points are local optimum solutions. Subsequently, the results of proposed hybrid methodology found that the optimal design parameters are at L1 = 6.3793 mm, L3 = 9.73 mm, and D = 2.25 mm. Using these, the proposed approach, the optimal frequency, and displacement are found at 27.006 Hz and 0.2203 mm, respectively. At last, the predicted MPCI from the proposed methodology is effective and better than that from the TMFL. In addition, the equivalent stress is also lower than the allowable stress. The optimal results are satisfied with designer’s requirements.

5. Results and Discussion

As discussed above, three cases studies are taken as numerical examples to describe application capacity and effectiveness of the proposed method. The performances of proposed methodology are compared with the TMFL. Table 23 summarizes all the optimal solutions for this comparison. In each case study, the optimal parameters are used to create a 3D model in SolidWorks 2018 software and then imported into ANSYS 2018 software to evaluate the robustness of the proposed hybrid approach by calculating the relative error. This error is calculated as follows:where ε, Rp, and Ra represent the relative error, predicted result from proposed method, and actual result, respectively.

Case study 1: using the TMFL, the results found the optimal frequency and displacement are 24.4574 Hz and 0.1691 mm, respectively. Meanwhile, the proposed hybrid approach predicts the frequency and displacement are about 26.592 Hz and 0.1679 mm, respectively. By using the proposed method, the relative errors between predicted values and FEA values are very small (around 5%), while the error value is around 7% via using TMFL. Both the proposed methods are reliable tools in this case. According to the Taguchi method, a larger-the-better value of MPCI is a desired value for the MOO problem. On the contrary, the proposed hybrid approach has a better performance than the TMFL because the MPCI from the proposed method is higher than that from the TMFL.

Case study 2: through the TMFL, the results found the optimal frequency and displacement are 28.9482 Hz and 0.1527 mm, respectively. By using the proposed hybrid integration, the optimal frequency and displacement are found at 651.08 Hz and 0.1981 mm, respectively. By using the proposed hybrid approach, the relative errors between predicted values and simulated values are around 4.65%, while the error between the predicted value by the TMFL and FEA value is around 9.5%. In addition, the MPCI predicted by the proposed hybrid approach is greater than that from TMFL. In other words, the performance and prediction accuracy of proposed methodology outperforms the TMFL.

Case study 3: the results found the optimal frequency and displacement are 158.970 Hz and 0.0957 mm, respectively, by using TMFL. Meanwhile, the proposed hybrid approach finds the optimal frequency and displacement approximately 651.08 Hz and 0.1981 mm, respectively. The relative errors between predicted values and simulated values are around 6%, while the error between the TMFL and FEA is about 484%. Moreover, the MPCI predicted by the proposed hybrid approach is higher than that from TMFL. It can conclude that the performance and precision of the proposed approach outperforms the TMFL.

The results of three numerical examples show that the proposed optimization scheme is better than the traditional TMFL. Among three numerical examples, case study 3 may be chosen as an optimal solution for MOO design of the SRCM because the MPCI is the highest value. The optimal design variables are at L1 = 6.3793 mm, L3 = 9.73 mm, and D = 2.25 mm. Other remain parameters are constant values. It concludes that the proposed hybrid approach is a robust optimization tool and effectiveness for solving the MOO design for the mechanism and related complex optimization fields.

6. Comparison with Metaheuristic Optimization Algorithms

As discussed above, the proposed method outperformed the traditional TMFL to solve MOO design for the mechanism. Continuously, the proposed algorithm is compared with recent less-parameter optimization algorithms such as Jaya algorithm and TLBO algorithm. Case 3 is chosen as the final optimal candidate for the proposed mechanism. Therefore, the ANFIS structure of case 3 is then coupled with Jaya and TLBO algorithms for this comparison. The results indicated that the optimal frequency and displacement predicted by the proposed algorithm are better than those estimated by the Jaya and TLBO algorithms, as given in Table 24. On the contrary, the proposed methodology has a prediction accuracy better than other metaheuristic optimization algorithms.

In order to evaluate and validate the effectiveness of three abovementioned algorithms, the nonparametric tests are used for comparing the performances of the proposed hybrid approach with other algorithms [66, 67]. The Wilcoxon’s rank signed test and Friedman test, two performance metrics, are utilized as the nonparametric tests. Both these statistical techniques consider that both different algorithms behave unsimilarly. Recently, the nonparametric tests have been widely applied in comparing the performances between algorithms, e.g., [68, 69]. The comparing results are almost based on the value. In the present article, the computational simulations are conducted with 60 runs for each algorithm. The Wilcoxon’s rank signed test is performed for frequency and displacement at 5% significant level and 95% confidence intervals. The statistical analysis is implemented in Minitab 18 software. The results in Tables 25 and 26 show that the value is less than 0.05 (5% significance level). According to nonparametric tests [66, 67], it reveals that there is a statistical difference between the proposed approach and other algorithms. Besides, the results confirm that the proposed hybrid algorithm has a performance effectiveness greater than ANFIS-based Jaya and ANFIS-based TLBO algorithms.

Another way, the Friedman test, is used for the behavior comparison of the abovementioned algorithms. The results find that the value is less than 0.05. Therefore, it confirms that there the proposed hybrid approach differs from two other algorithms, as given in Tables 27 and 28. Moreover, the effectiveness of the proposed hybrid algorithm is better than ANFIS-based Jaya algorithm and ANFIS-based TLBO algorithm. Besides, the median responses for the proposed algorithm are greater than those of two other algorithms. So, the proposed approach is a more effective tool than the other algorithms.

Although the proposed hybrid computational method has specific advantages in terms of solving MOO design of the SRCM, it has still some drawbacks: the running time is also long because the datasets are transferred into different stages. This may take a long time and high computational cost for more complex structures. It is known that the time and space complexity of an algorithm is very important. For a more complicated problem with large datasets, it needs more better memory space of a computer. These disadvantages of the proposed method can be solved as follows: (i) FEM analysis is time consuming, but this problem can be solved by using parallelism computation by using multiple processors or a high-performance computer. (ii) Automatically adaptive combination of statistical method, FEM, and metaheuristic algorithms is a feasible approach to reduce the computational cost. In other words, the proposed optimization scheme should automatically update new data if there is an adjustment of space or range of design variables. The aim of the present work is to recommend a new hybrid computational approach in solving the MOO design for compliant mechanisms. In future work, the time, space, and algorithm complexity need to be investigated to improve the proposed approach.

7. Conclusions

A new hybrid computational method is developed to solve a MOO design for a Scott Russell compliant mechanism in this article. The methodology can be considered as a combination of statistical techniques, finite element method, and intelligent computation. It is an integration of the desirability function technique, fuzzy logic reasoning, ANFIS, and LAPO algorithm. A 3D FEM model is created, and CCD is used to build the experimental matrix. Next, numerical data are collected through FEM simulations. Based on the numerical results, a sensitivity investigation is analyzed by using ANOVA and Taguchi method. The refinement process for design variables is to determine really searching spaces of design parameters. As a result, three case studies are assigned as application example of the proposed approach. Next, the desirability of two responses is computed and transferred into the FIS system. ANFIS technique is then utilized to model the MPCI, which is considered as a combined objective function. LAPO algorithm is adopted for optimizing the MPCI. The results find that the predicting accuracy of proposed hybrid methodology is better than that of the traditional Taguchi-based fuzzy logic. Case study 3 is chosen to be optimal candidate. Finally, the Wilcoxon signed rank test and Friedman test are analyzed to compare the performances of the proposed approach with the ANFIS-based Jaya algorithm and ANFIS-based TLBO algorithm. The results found the performance efficiency of proposed approach is better than those of other approaches. The proposed methodology can be an efficient tool for more complex optimization problems, including multiple objective functions, design variables, and constraints.

The results of the present study will be extended to more complex compliant mechanisms. Future research will concentrate on fabricating physical prototypes, and then an extra number of testing experiments is performed to validate the theoretical models and simulations.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant no. 107.01-2019.14.