Research Article  Open Access
Ngoc Le Chau, ThanhPhong Dao, Van Thanh Tien Nguyen, "An Efficient Hybrid Approach of Finite Element Method, Artificial Neural NetworkBased Multiobjective Genetic Algorithm for Computational Optimization of a Linear Compliant Mechanism of Nanoindentation Tester", Mathematical Problems in Engineering, vol. 2018, Article ID 7070868, 19 pages, 2018. https://doi.org/10.1155/2018/7070868
An Efficient Hybrid Approach of Finite Element Method, Artificial Neural NetworkBased Multiobjective Genetic Algorithm for Computational Optimization of a Linear Compliant Mechanism of Nanoindentation Tester
Abstract
This paper proposes a new evolutionary multiobjective optimization technique for a linear compliant mechanism of nanoindentation tester. The mechanism design is inspired by the elastic deformation of flexure hinge. To improve overall static performances, a multiobjective optimization design was carried out. An efficient hybrid optimization approach of central composite design (CDD), finite element method (FEM), artificial neural network (ANN), and multiobjective genetic algorithm (MOGA) is developed to solve the optimization problem. In this procedure, the CDD is used to lay out the experimental data. The FEM is developed to retrieve the quality performances. And then, the ANN is developed as black box to call the pseudoobjective functions. Unlike previous studies on multiobjective evolutionary algorithms, most of which generating only one Paretooptimal solution, this proposed approach can generate more than three Paretooptimal solutions. Based on the user’s realwork problem, one of the best optimal solutions is chosen. The results showed that the optimal results were found at the displacement of 330.68 μm, stress of 140.65 MPa, and safety factor of 3.6. The statistical analysis is conducted to investigate the behavior of the MOGA. The sensitivity analysis was carried out to determine the significant contribution of each factor. The results revealed that the lengths and thickness almost significantly affect both responses. It confirms that the proposed hybrid optimization approach gains high robustness and effectiveness with flexible decision maker rules to solve complex optimization engineering problems.
1. Introduction
Nanoindentation tester has attracted great interest from researchers from academics and industry. This device is now widely used in areas of material science. It is desired to probe the mechanical properties such as harness, creep, elasticplastic modulus, and roughness surface. Materials can be tested, including hard and soft types from tissue, biological cell, nanomaterial, optics, material science, semiconductor, biomechanics, microelectromechanical systems, and electronics [1].
During the indentation process, multiple microscopes are used to record the image of sample before and after indenting test to characterize the curve of displacement versus load. High precision positioning and stability of the system are dependent not only on the microscope, controller, and imaging technique, but also a mechanicaldriven platform. The platform is utilized to bring the material in front of microscope and then a picture is taken. And then, the platform takes the material to the location of indenter so as to do an indentation. At last, the platform brings the material to return back the location of microscope to characterize the curve of displacement versus load. It can be concluded that the mechanicaldriven platform is a critical mechanism for the nanoindentation tester. In commercialization, the current instruments are relatively effective but so costly. The reason is that, in the current nanoindentation tester system, a ballscrew driven system by servo motors or linear motors is used to ensure the linear precision. However, such mechanicaldriven systems possess some common defects such as assembly gap, wear and aging, and friction error, which may restrict their further applications in the nanoindentation industries where micro or nanoscale positioning accuracy is needed. Hence, it has been a difficult challenge for researchers to make a cheaper system with better accuracy.
To eliminate the cost and enhance the positioning accuracy, a linear compliant mechanism (LCM) is proposed in this study to bring the observed material. The proposed LCM is operated based on ideal of elastic springs and has no backlash, free friction, high precision, and monolithic structure [2–5]. Meanwhile, the servodriven system causes large backlash and friction between kinematic joints. To create more precise and light positioning system, the LCM is recommended as a potential candidate. There have not been studies on the LCM for positioning the material in nanoindentation tester yet. LCM is now used to handle displacement from a few micrometers to hundreds of micrometers. The proposed LCM can meet a good precision positioning but its workspace and life time are limited because its operation relies on elastic links. In addition, the common piezoelectric actuatorguided nanopositioning stages have the problems of small travel range because the piezoelectric actuator has a limited travel range of less than 200 μm.
To meet the practice requirements of advanced material sciences, the capability of a large positioning space and long working life of the LCM should be improved further. To overcome the limitation of performances, common procedures are, for example, topology and size optimization. Topology represents the connectivity of the domain. Size optimization denotes design parameters such as thickness, crosssectional area, and length [6–8]. Shape optimization represents a suitable configuration. In general, mathematical equations are often established, and then an evolutionary optimization algorithm is used to seek a set of optimal design parameters [9]. If the established mathematical equations are not right, the optimized results will not be accurate.
For these reasons, this study introduces a new evolutionary multiobjective optimization technique for optimizing the performances of LCM to decrease the modeling errors. The proposed optimization approach is a hybrid intelligent integration of central composite design (CDD), finite element method (FEM), artificial neural network (ANN), and multiobjective genetic algorithm (MOGA). Each of the mentioned methods has been used widely but individually. For example, CDD was investigated in building and optimizing [10–15], FEM was used in various engineering fields [9, 16–18], ANN was employed for thermal or other engineering disciplines [19–23], and MOGA was utilized for optimizing multiple responses [24–26]. Although these methods were well known, hybridization of them has not been developed for optimization of the LCM. Therefore, an integration of hybrid optimization approach is a main purpose of this study. In addition, most of previous multiobjective evolutionary algorithms (MOEAs) often give a Paretooptimal front. Unlike previous MOEAs, the proposed hybrid algorithm can generate more number of optimal solutions. This helps users choose the best solution for realwork purpose.
The computationally evolutionary optimization approach is developed to satisfy a tradeoff between performance characteristics of LCM. In this hybrid optimization approach, a model is created via using finite element method. Subsequently, based on central composite design (CCD), a number of simulated experiments are generated and the results of performances are automatically generated via FEM. And then, relationship between design parameters and output performances will be generated well by using the ANN technique. ANN is as a black box to map well the inputs and outputs. Finally, MOGA is integrated to seek an optimal solution. This optimization process is totally automatic solution to gain accurate results.
The purposes of this paper are to propose an optimal design strategy for the LCM in terms of good static characteristics. Elastic joints of the LCM are designed based on elastic springs. A leverdisplacement amplifier is used to amplify the workspace. To improve the static performances, a hybrid optimization algorithm is developed. A validation is performed to evaluate predicted results, efficiency, and robustness of the proposed approach.
2. Design of Linear Compliant Mechanism
2.1. Compliant Displacement Amplification
Although the piezoelectric stack actuator (PSA) was a common actuator for precision positioning platform and manipulator, it is still limited by its working travel. To overcome this limitation, the leverdisplacement amplification is often used to amplify the travel of PSA.
A traditional leverdisplacement amplifier is depicted, as seen in Figure 1(a). The length of lever is divided into two parts with input distance L_{1} and output distance L_{o}. Figure 1(b) illustrates the kinematic diagram for operation principle. When the input load is exerted at the input port, point A moves to A’. At the same time, the output port is moved from point B to B’. The vertical link serves as elastic joint or flexure hinge.
(a)
(b)
The traditional amplification ratio can be calculated as follows:
in which, δ_{1} and δ_{2} are the input displacement and output displacement, respectively.
In the field of precision positioning system, compliant mechanisms (CMs) are operated relying on the elastic deformation of thin cross section links [27]. CMs are largely dependent on the lump or distributed elastic elements. These elastic parts work under the yield strength while the deformation of flexure hinges is eliminated. Hence, a compliant displacement amplification (CDA) mechanism with two levers was proposed in this study. The CDA included a lever in the left and another lever in the right. The CDA was intended to amplify the work travel of linear compliant mechanism (LCM). As seen in Figures 2(a) and 2(b), when a force was applied at the input port of the CDA, the movement of the output port of CDA was amplified.
(a)
(b)
The basic mathematical model for the CDA was computed as follows:
in which, n is the number of used levers.
As seen in (2), the amplification ratio of CDA was n times the traditional amplification ratio in (1).
As depicted in Figure 2, the LCM was only desired to move along the yaxis. The links with thickness t_{1} and length l_{1} served as flexure hinge and helped to eliminate the parasitic motion along the xaxis. Other links with thicknesses t_{2} and t_{3} and lengths l_{2} and l_{3} were geometrical parameters of CDA. The LCM had width of . Overall design parameters were the most important parameters affecting the performances of the precision positioning system. The Al 7075 was selected as material for the LCM. The parameters of material and design parameters were given, as in Table 1.

2.2. Primary Application for Nanoindentation Tester
In the material field, inspecting mechanical properties of a soft and thin material sample during nanoindentation test is critical work. Figure 3 illustrates a potential application of the LCM for nanoindentation tester. It includes 5 important parts with different functions. The top plate was connected with the base to serve as fixed stationary. A microscope was fixed on the top plate so as to monitor the material sample. The LCM was used to bring the material sample. At the beginning, the sample was pictured by the microscope. The microscope was used to drive a microscope/lens so as to observe the priorandlater surface of material. And then, the imaged pictures are analyzed using image technique. As a result, the mechanical properties can be extracted from the curve of loaddisplacement and this type of curve should ensure accurately. And then, the sample was moved to a location of the indenter. At last, the indenter would indent inside the sample and at the same time the microscope would defect the indentation depth and specifications of the deformed sample.
3. Formulation of Multiobjective Optimization Problem
Existing nanoindentation tester has attracted a great interest in material engineering. It is used for various materials such as biological tissues and soft gels. Commercialized instruments can monitor easily the mechanical properties of a sample but its operating costs are rather high. A high cost may come from positioning system where motors and actuators are utilized simultaneously. To fulfill higher precise probe requirements in several tens of micrometers to hundreds of millimeters, the proposed LCM was used to decrease the motors and actuators. Hence, the considerable challenges have been faced when developing the nanoindentation devices; the instruments should have the following specifications: (i) able to give a large displacement; (ii) able to speed up the indentation process; (iii) able to ensure a long working time; (iv) able to gain an excellent stability; and (v) able to record the displacement and force signal itself before and after indentation process.
Specifically, the LCM could give a large displacement and a long working time for nanoindentation tester. In this system, the LCM should fulfill the following important requirements: to satisfy faced challenges as(a)a large working stroke to allow suitability for various materials;(b)a high safety factor to enhance the fatigue life.
In general, these requirements conflicted together. A tradeoff multiple responses could be gained by developing a new hybrid optimization approach.
3.1. Design Variables
It was noted that the proposed LCM was very sensitive to geometric dimensions such as thicknesses, lengths, and width. Therefore, these parameters were considered as the design variables. The vector of design variables was set as .
The lower and upper bounds for the design variables were assigned as follows:where t represents the thickness of flexure hinges, l_{i} represents length of ith flexure hinge with , and is the width of LCM.
3.2. Objective Functions
The multiple quality performances of LCM were as follows: (i) the displacement, , was desired as large as possible; (ii) the safety factor, , is required as high as possible. To summarize, the optimization problem was briefly stated as follows:
3.3. Constraints
The resulting stress of the LCM must be under the yield strength of the material, which was described as follows:where is the resulting stress and is the yield strength of proposed material s.t.
4. Hybrid Optimization Algorithm
Before conducting an engineering optimization problem, mathematical models are traditionally established, and then a suitable optimization algorithm is utilized. The fact is that there would be an error between mathematical models and simulations or experiments. This is because those models may not be right and depend mainly on capability and knowledge of researchers about engineering and mathematical theory. Therefore, an evolution optimization algorithm can then generate inaccurate results. To overcome this limitation, a hybrid approach of RSM, FEM, ANN, and MOGA was proposed in this study so as to improve the quality performances of the LCM. The robustness and efficiency of proposed approach could be validated through simulated experiments. Figure 4 illustrates a systematic flowchart for the multiobjective optimization procedure. The optimization process was carried out by three large phases and substeps as follows.
Phase 1 (computer aid engineering application and design). Nowadays, with a fast development of computer, a high performance computer (HPC) can perform a thousands of calculations per a second, which is considered as a workstation. Based on these advantages, machine learning and artificial intelligent algorithms can be programed well to meet a research and development in engineering. It can be called computer aid engineering application and design (CAD A&D). The first stage to optimize the CDA, a CAD A&D process for the CDA, was conducted by following steps.
Step 1 (define problem). The CDA was deigned to serve as a displacement amplifier for applications in precision engineering micro or nanopositioning systems. In these devices, a piezo stack actuator (PSA) serves as an actuator but its travel is limited. Therefore, the CDA was intended to amplify the work travel for PSA, a product from Physik Instrumente Co., Ltd. Multiobjective optimization problem for the CDA was conducted to improve its quality performances.
Step 2 (mechanical structure). A mechanical structure was proposed based on designer’s experience and perception. Initially, many drafts were drawn to illustrate the highlights of CDA. A final model was chosen but its specifications cannot meet the requirements of positioning system.
Step 3 (define design variables and objective functions). The lengths and thicknesses of flexure hinge and the width of the CDA were determined as design variables. The reason was that these parameters affect the performances of the CDA. These variables can be seen as Figure 1.
To meet the practical requirements of customers, a precision positioning system must have a large workspace, a high speed, a high safety factor, and a minimum stress. All these performances can be fulfilled by optimizing the CDA because the mechanical structure and geometries of CDA contribute to the system. Hence, the mentioned quality characteristics were taken into account as four objective functions.
Step 4 (build 3D model). A 3D model was drawn by using FEM. The FEM is a computational technique used to obtain approximate solutions of problems in engineering. Postprocessor contains sophisticated routines used for sorting, printing, and plotting selected results from a finite element solution. Because the CDA was an elastic structure, during the analysis some of the following relationships were used.
According to the linear Hooke’s law, the equation describing a relation between the force and displacement was computed aswhere is the global stiffness matrix that is computed according to finite element principle. Correspondingly, the global displacement vector and force vector for all nodes can be calculated aswhere is the total number of nodes and and are the displacement and force vectors of node i, which can be determined bywhere and are the displacement along the x and yaxes of node ith. is the rotational angle of node ith. and are applied forces and is moment.
Also by Hooke’s law, the relationship between stress and strain was calculated as follows:where , , are the stress and Young’s modulus, sequentially.
Step 5 (evaluate initial performances). The 3D model was drawn in Step 4 and considered as a draft design. It had four quality performances, as given in Step 3; namely, a precision positioning system must have a large workspace, a high speed, a high safety factor, and a minimum stress. Such performances should be evaluated by implementing FEM simulations. This step ended the CAD A&D process. If the CDA’s specifications were not satisfied according to the designer’s requirements, the process would return to Step 2. If it was ok, it would move to Phase 2.
Phase 2 (response surface and regression model). Prior to optimization process, virtually pseudomathematical equations could be treated as fitness functions. To do this, a number of experiments were generated and simulated data were collected. At last, the regression models were established to map design variables and outcome performances.
Step 6 (design of experiment). By using the RSMintegrated FEM, a numbers of simulation experiments were planned via using the central composite design integrated with RSM. The number of necessary experiments was determined by the following equation:where was the total of the design points, was the number of design variables, was the factorial number (f = 0), and n_{c} = 1 was the number of replicates at the center point of the design space.
Step 7 (generate data). The simulated experimentations were collected based on an integration of RSM and FEM. First, they were based on the 3D model designed in the FEM in the Step 4. Subsequently, the estimated results for the quality responses were retrieved based on finite element analysis (FEA).
Step 8 (establish regression models). There were various regression models such as full 2ndorder polynomials, artificial neural network, and nonparametric regression. In this study, ANN was a suitable candidate for the retrieved data. It was used as a regression approach for the estimated database. It was considered as a black box to approximate the complex nonlinear relationship between design parameters and the qualities. ANN model can find a pseudoobjective functions. These objective functions were used for the MOGA algorithm.
Phase 3 (operation of MOGA). After the pseudoobjective functions were determined in Step 8, the optimization process was implemented by setup controllable parameter and programing MOGA algorithm. This algorithm was proposed for multiobjective optimization problem in this study because it can converge to the global Pareto solutions. This algorithm was used to seek a Paretooptimal set for multiobjective optimization problem. A flowchart of proposed MOGA procedure was illustrated in Figure 5.
Step 9 (initialize controllable parameters of MOGA). The multiobjective genetic algorithm (MOGA) was a variant of the nondominated sorted genetic algorithmII (NSGAII) that relied on controlled elitism concepts. It can solve multiple objectives and constraints. At last, the MOGA can find the global optimum solution. The controllable parameters of MOGA in this study were given in Table 2. This algorithm includes some important parameters such as population size, crossover probability, mutationprobability, and number of generations. To achieve an accurate optimal solution, correct selection of the parameters (crossover, mutation, size of population, and number of generations) would be conducted later. Table 2 gives the range of tuning parameters of the MOGA.

Step 10 (computing the fitness value). This step was critical to change the fitness so as to gain the fitness values of the individuals.
Step 11 (selection operator). The usual selection operator, such as Monte Carlo, was used and combined with the elitist model.
Step 12 (crossover operator). Crossover operation was utilized to produce parent generation.
Step 13 (mutation operation). Mutation operation was used to generate new generation.
Step 14 (termination criteria). If the gross generation was gained and the best individual was found through predetermined iterations, the MOGA was stopped. Otherwise, Step 10 was repeated.
Phase 4 (evaluate optimum candidates). The best candidates were found by using MOGA and then the final phase was overall evaluations of these candidates.
Step 15 (find the Paretooptimal set). Multiobjective optimization problem was solved by using MOGA, and then the tradeoff the objective optimizations was treated as Paretooptimal set. At last, the optimal results were found.
Step 16 (evaluate optimal candidates). If the optimal results were satisfying, the extra validations were conducted to evaluate the robustness and efficiency of the proposed hybrid optimization approach. The optimization process was ended herein. If they were not well refined, the optimization process would be further enhanced by Step 17 or Step 18.
Step 17 (refine objective functions). In performing the optimization process, a refinement was the most important phase to seek the optimal solutions. After optimal candidates were generated and based on the initial requirements of quality characteristics, the researchers would evaluate the candidates. If there was no candidate that was satisfying, the ranges of quality characteristics or the range of design variables must be controlled or refined again from Step 15. This step was repeated until the best candidate was found.
Step 18 (refine range of design variables). If Step 13 was not well turned, the optimization process would act on the range of design variables. The ranges could be increased or decreased. Steps 2–12 were repeated. Otherwise, Step 13 would be repeated.
5. Results and Discussion
5.1. Test Initial Performances
A 3D model was created in FEM design modeler, and then FEA simulations were conducted to test the initial performance of LCM. As shown in Table 3. The results showed that the displacement is about 295.59 μm and safety factor is approximately 2.733. The two values are not satisfying for the desired requirements by constraints in (7)(8). Therefore, to meet the used actual requirement of nanoindentation tester, the two characteristics should be further optimized.

5.2. Experimental Data
As proposed hybrid optimization approach, a 3D model of the LCM was designed via using FEM. With six design variables, the numbers of experiments were calculated using the central composite design in (15).
And then, based on Kriging regression model, each pseudoobjective function was made for the displacement, resulting stress, and safety factor. Finally, by integrating FEM, RSM, and ANN regression model, the real values of three quality responses were automatically computed, as given in Table 4.

5.3. Regression Model
Compared with other regression models such as full 2ndorder polynomials, nonparametric regression, and Kriging model, artificial neural network gave relatively correct results. The coefficient of determination of ANN model was approximately 1, as given in Table 5. This value was relatively well for a regression model. It was considered as a performance criterion. To sum up, the built ANN can be regarded as the pseudoobjective functions that described the relation between design parameters and the displacement and safety factor. The data from the table was randomly divided into training of 70% of data, validation of 15% of data, and testing of 15% of data.

The structure of an ANN was determined by the number of layers, the number of nodes in each layer, and the nature of the transfer functions. This paper used threelayered feed forward back propagation neural network (7:10:1), as shown in Figure 6.
In this study, the Bayesian Regularization was chosen as the training method for the safety factor data because it permitted a coefficient of determination (R^{2}) of approximately 1. Meanwhile, to gain the R^{2} of 1, the LevenbergMarquardt was suitable for the displacement data. Training, validation, and testing all had a mean square error of about zero, as shown in Figure 7 and Table 5. The results showed that the performances of the proposed ANN models are good for the data.
(a)
(b)
The contribution of each parameter and interaction affecting the displacement was analyzed in Table 6. The results showed that the length l_{2} has a largest contribution of 97.78% with respect to the displacement and statistical significance with pvalue of 0.000 (less than 0.05%).

As shown in Table 7, the length l_{2} also had a largest contribution on the safety factor with 50.14% and a good statistical significance with pvalue of 0.000 (less than 0.05%).

5.4. Optimal Results
By using Phase 3 and begining from Step 5, MOGA was integrated with FEM, RSM, and ANN to seek a tradeoff between the displacement and safety factor. During the optimization process, the optimal results were generated automatically. Many of potential candidates would be retrieved but some might not good and could not find any value. Therefore, this study proposed two ways to gain the candidates. The first way was an adjustment on the range of design parameters. The second way was a change in the range of objective functions. After implementing many optimization processes to conduct and prove the two suggestions, the results showed that this hybrid optimization could be effective only if the range of two quality responses was limited well. It was mainly dependent on the designer’s experiences. To suppress this dependence, there was two common rules as follows:(i)The largerthe better objective function: the range of this function should not exceed an allowable value. For example, the maximum displacement of LCM was desired to be about 400 μm, and therefore the range of this response was controlled in the range of 400 μm, as seen in Table 8.(ii)A minimum high safety was also required, and this function was constrained to be 3, as seen in Table 8.

The history charts of displacement, stress, and safety factor were retrieved from the results of the proposed hybrid optimization algorithm. As seen in Figure 8, the displacement was convergent in the range from 390 μm to 420 μm. Meanwhile, the safety factor was convergent in the range from 3 to 3.5, as shown in Figure 9. It means that the optimal results can be found in the range of desired space. As can be seen in these figures, the curves may be so noisy. It is true that there were so many data points occurring at the optimum space and the optimal cost functions were found in the ranges. Hence, they could be very noisy. This was a new approach compared with previous studies where there was only a convergent solution. The convergence in a range allowed the choice the best solution suitable for a realwork problem.
The pseudo or virtual mathematical models were found by using ANN approach. To achieve the optimal results for the LCM, the constraints for two objective functions were set up, as given in Table 5. Three potential candidates were generated. As given in Table 9, the potentially optimal design parameters were found with similar values for each factor. How to choose the best candidate was dependent on user’s requirements.

Before selecting the best optimal solution for the LCM, the best tuning parameters of MOGA algorithm were selected. As given in Table 2, the tuning parameters for MOGA were set in this study as follows: the population size is within the range of 20100, the crossover probability is within the range of 0.20.90, the mutation probability is within the range of 0.010.02, and the number of generation is within the range of 50150. To find the best parameters for the MOGA, after thirty times of computational simulations, the best parameters of MOGA were chosen as follows: population size of 100, crossover probability of 0.7, mutation probability of 0.01, and number of generation of 100.
A shown in Table 10, candidate 1 was chosen as the best optimal design because it fully satisfies the mentioned design objectives and the resulting stress of 140.6496 MPa was lowest and highest safety factor of 3.5952. The minimum stress could guarantee a fatigue life and long working time. The predicted results indicated that the optimal results were found at the displacement of 398.5 μm, a minimal stress of 140.65 MPa, and safety factor of 3.6.

The Paretooptimal fronts and results obtained from the MOGA algorithm were than compared to those obtained using the Elitist Nondominated Sorting Genetic Algorithm (NSGAII) [28]. Common parameters for NSGAII are included: the population size is within the range of 20100, the crossover probability is within the range of 0.20.9, the mutation probability is within the range of 0.10.2, and the number of generation is within the range of 50150.
A basic way to find good parameters for the algorithm is preliminary runs. After thirty times of computational simulations, the best parameters of NSGAII were chosen as follows: population size of 75, crossover probability of 0.8, mutation probability of 0.2, and number of generation of 110. The comparative results showed that the optimal performances from the MOGA are better than those obtained from the NSGAII. Specifically, the stress from MOGA was smaller than that from NSGAII, and the displacement and safety factor from MOGA were higher than those from NSGAII. Besides, the computational optimization time required for MOGA was less than that for NSGAII, as given in Table 11.

5.5. Statistical Analysis
In order to evaluate the behavior of the evolutionary algorithms, a statistical analysis is often used. In this study, the Wilcoxon’s rank signed test was applied to describe the behavior of the MOGA. It was a statistically nonparametric technique applied for various fields [29, 30]. For details about this procedure, the readers can refer to [29, 30]. To investigate the behavior of the MOGA, the NSGAII was used simultaneously to make the two algorithm pairs. A comparison of two algorithms aimed to discover the significant difference between their behaviors. The computational simulations were conducted with 30 runs for each algorithm. The Wilcoxon’s rank signed test was performed at 5% significant level. The results of Wilcoxon’s rank signed test were given in Table 12. It showed three parameters, namely, R, R+, and pvalue. The MOGA showed significantly better results than the NSGAII with pvalue of 0.006.

5.6. Analysis of Sensitivity
Along with optimization process, a sensitivity analysis is also a necessary step to determine an influence of each design parameter on each quality response. For the proposed LCM, the sensitivity analysis was mainly focused on three characteristics. Generally, there were many methods that can be applied for calculating the sensitivity such as Nelson method, modal method, matrix perturbation method, differential method, RSM, or statistical analysis [31, 32]. For example, the direct differential method takes more time for analysis because it needs to construct physical models. Therefore, RSM was chosen for this study. The sensitivity could be calculated by following formula: where f_{i}, are the quality response and the design parameter ith, respectively.
In this study, seven variables and two quality characteristics were considered. As seen in Figure 10, regarding the displacement’s sensitivity, thickness t_{3} had a highest influence or significant contribution on the displacement while the length l_{3} had a lowest contribution. It was noted that a change in t_{3} would adjust the displacement as desired. Other parameters were in relatively middle effects. Regarding the safety factor, the lengths l_{2} and l_{3} had the smallest contributions while the length l_{1} had the largest effect on the safety factor. To adjust the safety factor, the l_{1} should be changed firstly.
Particularly, Figure 11(a) illustrates the effects of the width and length l_{2} on the displacement. The results indicated that the displacement has a linear change corresponding to the width while it has a nonlinear influence with respect to the length. Figure 11(b) shows the effects of the width and length l_{2} on the safety factor. The results indicated that the safety factor has a nonlinear effect corresponding to the width and it also has a nonlinear influence with respect to the length. In particular, the length had a more significant contribution to the safety factor compared with the width because this response was changed sharply. An increase in the length made a decrease in the safety factor.
(a)
(b)
The contribution diagrams of thicknesses t_{2} and t_{3} were plotted in Figures 12(a) and 12(b), respectively. The results revealed that both parameters almost have a nonlinear influence on the displacement. When the thickness was increased, this response was lowered, as seen in Figure 10(a). Meanwhile, both design variables had a relatively linear contribution on this response. It means that the safety factor was linearly raised with respect to these parameters, as given in Figure 12(b).
(a)
(b)
As depicted in Figure 13(a), the lengths l_{1} and l_{3} had a nearly linear influence with the displacement but they had a nonlinear contribution affecting the safety factor.
(a)
(b)
As seen in Figure 14(a), the thickness t_{1} had a nonlinear influence on the displacement as well as the safety factor.
(a)
(b)
Summary, almost the mentioned design parameters had significant contributions to the displacement and safety factor. This would help designers and researchers to make a decision and meet the requirements of a specific system.
6. Verifications
To evaluate and validate the optimal performances of the proposed LCM, simulation tests were carried out. Using the optimal design variables of candidate 1 in Table 6, the prototype of LCM was created. Extra experimental validations were conducted and the average value was retrieved.
As given in Table 13, the results indicated that the optimal results were found at the displacement of 398.5 μm and safety factor of 3.5952. The results found that the predicted and experimental results are in a good agreement. It means that the proposed hybrid optimization method is actually robust approach to solve multiobjective optimization problem for the LCM. It can be applied to solve complex optimization problems. Compared with the initial design, the displacement and safety factor were improved by about 2.5% and 4.7%, respectively, as shown in Table 14.


7. Conclusions
This paper presented a new intelligent evolutionary multiobjective optimization approach for a linear compliant mechanism. The mechanism was designed based on connecting series of the leaf springs. These springs were located in symmetric configuration not only to guarantee a motion linearity but also to increase the working travel. To improve overall static performances, including the large working travel and high safety factor, a hybrid optimization approach was developed. This approach was an integration of FEM, RSM, ANN method, and MOGA. Three optimal candidates were retrieved and then candidate 1 was chosen as the finally optimum solution.
The sensitivity analysis was carried out to determine the significant contribution of each factor. The results revealed that the lengths and thickness are main influencers. The results revealed that the lengths and thickness almost significantly affect both responses. The results showed that the optimal results were found at the displacement of 330.68 μm, stress of 140.6496 MPa, and safety factor of 3.5952. The predicted results were highly consistent with the experimental results. It was confirmed that the proposed hybrid optimization approach is highly robust to solve complex optimization engineering problems.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
This research is funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 107.012016.20.
References
 W. O'Brien, “Longrange motion with nanometer precision,” Photonics Spectra, vol. 39, no. 6, pp. 8081, 2005. View at: Google Scholar
 S. Polit and J. Dong, “Development of a highbandwidth XY nanopositioning stage for highrate micro/nanomanufacturing,” IEEE/ASME Transactions on Mechatronics, vol. 16, no. 4, pp. 724–733, 2011. View at: Publisher Site  Google Scholar
 M.G. Song, H.W. Baek, N.C. Park et al., “Development of small sized actuator with compliant mechanism for optical image stabilization,” IEEE Transactions on Magnetics, vol. 46, no. 6, pp. 2369–2372, 2010. View at: Publisher Site  Google Scholar
 T.P. Dao and S.C. Huang, “Compliant thinwalled joint based on zygoptera nonlinear geometry,” Journal of Mechanical Science and Technology, vol. 31, no. 3, pp. 1293–1303, 2017. View at: Publisher Site  Google Scholar
 R. F. Fung, Y. L. Hsu, and M. S. Huang, “System identification of a dualstage XY precision positioning table,” Precision Engineering, vol. 33, no. 1, pp. 71–80, 2009. View at: Publisher Site  Google Scholar
 T.P. Dao and S.C. Huang, “Design and multiobjective optimization for a broad selfamplified 2DOF monolithic mechanism,” Sādhanā, vol. 42, no. 9, pp. 1527–1542, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 S.C. Huang and T.P. Dao, “Design and computational optimization of a flexurebased XY positioning platform using FEAbased response surface methodology,” International Journal of Precision Engineering and Manufacturing, vol. 17, no. 8, pp. 1035–1048, 2016. View at: Publisher Site  Google Scholar
 S.C. Huang and T.P. Dao, “Multiobjective optimal design of a 2DOF flexurebased mechanism using hybrid approach of GreyTaguchi coupled response surface methodology and entropy measurement,” Arabian Journal for Science and Engineering, vol. 41, no. 12, pp. 5215–5231, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 E. Burman, D. Elfverson, P. Hansbo, M. G. Larson, and K. Larsson, “Shape optimization using the cut finite element method,” Computer Methods Applied Mechanics and Engineering, vol. 328, pp. 242–261, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 S. Baugreet, J. P. Kerry, A. Brodkorb et al., “Optimisation of plant protein and transglutaminase content in novel beef restructured steaks for older adults by central composite design,” Meat Science, vol. 142, pp. 65–77, 2018. View at: Publisher Site  Google Scholar
 A. R. Fares, A. N. Elmeshad, and M. A. A. Kassem, “Enhancement of dissolution and oral bioavailability of lacidipine via pluronic P123/F127 mixed polymeric micelles: Formulation, optimization using central composite design and in vivo bioavailability study,” Drug Delivery, vol. 25, no. 1, pp. 132–142, 2018. View at: Publisher Site  Google Scholar
 S. Mosleh, M. R. Rahimi, M. Ghaedi, K. Dashtian, and S. Hajati, “Sonochemicalassisted synthesis of CuO/Cu2O/Cu nanoparticles as efficient photocatalyst for simultaneous degradation of pollutant dyes in rotating packed bed reactor: LED illumination and central composite design optimization,” Ultrasonics Sonochemistry, vol. 40, pp. 601–610, 2018. View at: Publisher Site  Google Scholar
 F. Nemati, R. ZareDorabei, M. Hosseini, and M. R. Ganjali, “Fluorescence turnon sensing of thiamine based on Arginine – functionalized graphene quantum dots (ArgGQDs): Central composite design for process optimization,” Sensors and Actuators B: Chemical, vol. 255, pp. 2078–2085, 2018. View at: Publisher Site  Google Scholar
 A. Raheem, G. Ji, A. Memon et al., “Catalytic gasification of algal biomass for hydrogenrich gas production: Parametric optimization via central composite design,” Energy Conversion and Management, vol. 158, pp. 235–245, 2018. View at: Publisher Site  Google Scholar
 S. Shariati, E. Bozorgzadeh, F. Shariati, and F. Safa, “Ionic Liquid Based UltrasoundAssisted Emulsification Microextraction for Preconcentration of Phenol Using Central Composite Design,” Journal of Analytical Chemistry, vol. 73, no. 1, pp. 36–41, 2018. View at: Publisher Site  Google Scholar
 A. Albanesi, F. Bre, V. Fachinotti, and C. Gebhardt, “Simultaneous plyorder, plynumber and plydrop optimization of laminate wind turbine blades using the inverse finite element method,” Composite Structures, vol. 184, pp. 894–903, 2018. View at: Publisher Site  Google Scholar
 A. Karamanlı and T. P. Vo, “Size dependent bending analysis of two directional functionally graded microbeams via a quasi3D theory and finite element method,” Composites Part B: Engineering, vol. 144, pp. 171–183, 2018. View at: Publisher Site  Google Scholar
 N. Le Chau, V. A. Dang, H. G. Le, and T.P. Dao, “Robust Parameter Design and Analysis of a Leaf Compliant Joint for Micropositioning Systems,” Arabian Journal for Science and Engineering, vol. 42, no. 11, pp. 4811–4823, 2017. View at: Publisher Site  Google Scholar
 T. Zhang, Z. Chen, J. Liu, and X. Li, “The artificial neural networks based on scalarization method for a class of bilevel biobjective programming problem,” Comput Intell Neurosci, vol. 2017, Article ID 1853131, 14 pages, 2017. View at: Google Scholar
 A. Zamaniyan, F. Joda, A. Behroozsarand, and H. Ebrahimi, “Application of artificial neural networks (ANN) for modeling of industrial hydrogen plant,” International Journal of Hydrogen Energy, vol. 38, no. 15, pp. 6289–6297, 2013. View at: Publisher Site  Google Scholar
 Vijay G. S., Kumar H. S., Srinivasa Pai P., Sriram N. S., and R. B. Rao, “Evaluation of Effectiveness of Wavelet Based Denoising Schemes Using ANN and SVM for Bearing Condition Classification,” Computational Intelligence and Neuroscience, vol. 2012, Article ID 582453, 12 pages, 2012. View at: Publisher Site  Google Scholar
 B. A. Garro and R. A. Vázquez, “Designing artificial neural networks using particle swarm optimization algorithms,” Computational Intelligence and Neuroscience, vol. 2015, Article ID 369298, 20 pages, 2015. View at: Publisher Site  Google Scholar
 M. Hemmat Esfe, S. Esfandeh, S. Saedodin, and H. Rostamian, “Experimental evaluation, sensitivity analyzation and ANN modeling of thermal conductivity of ZnOMWCNT/EGwater hybrid nanofluid for engineering applications,” Applied Thermal Engineering, vol. 125, pp. 673–685, 2017. View at: Publisher Site  Google Scholar
 S. H. MousaviAvval, S. Rafiee, M. Sharifi et al., “Application of multiobjective genetic algorithms for optimization of energy, economics and environmental life cycle assessment in oilseed production,” Journal of Cleaner Production, vol. 140, pp. 804–815, 2017. View at: Publisher Site  Google Scholar
 S. Wang, G. Jian, J. Wang, L. Sun, and J. Wen, “Application of entransydissipationbased thermal resistance for performance optimization of spiralwound heat exchanger,” International Journal of Heat and Mass Transfer, vol. 116, pp. 743–750, 2018. View at: Publisher Site  Google Scholar
 J. Wen, X. Gu, M. Wang, S. Wang, and J. Tu, “Numerical investigation on the multiobjective optimization of a shellandtube heat exchanger with helical baffles,” International Communications in Heat and Mass Transfer, vol. 89, pp. 91–97, 2017. View at: Publisher Site  Google Scholar
 T.P. Dao, “Multiresponse optimization of a compliant guiding mechanism using hybrid taguchigrey based fuzzy logic approach,” Mathematical Problems in Engineering, vol. 2016, 17 pages, 2016. View at: Google Scholar
 K. Deb, S. Agrawal, A. Pratap, and T. Meyarivan, “A fast elitist nondominated sorting genetic algorithm for multiobjective optimization: NSGAII,” Lecture Notes in Computer Science, vol. 1917, pp. 849–858, 2000. View at: Publisher Site  Google Scholar
 E. Cuevas, J. Gálvez, S. Hinojosa, O. Avalos, D. Zaldívar, and M. PérezCisneros, “A comparison of evolutionary computation techniques for IIR model identification,” Journal of Applied Mathematics, vol. 2014, 2014. View at: Google Scholar
 S. García, D. Molina, M. Lozano, and F. Herrera, “A study on the use of nonparametric tests for analyzing the evolutionary algorithms' behaviour: a case study on the CEC'2005 Special Session on Real Parameter Optimization,” Journal of Heuristics, vol. 15, no. 6, pp. 617–644, 2009. View at: Publisher Site  Google Scholar
 N. L. Chau, T. Dao, and V. T. Nguyen, “Optimal Design of a DragonflyInspired Compliant Joint for Camera Positioning System of Nanoindentation Tester Based on a Hybrid Integration of JayaANFIS,” Mathematical Problems in Engineering, vol. 2018, pp. 1–16, 2018. View at: Publisher Site  Google Scholar
 N. L. Ho, T. Dao, H. G. Le, and N. L. Chau, “Optimal Design of a Compliant Microgripper for Assemble System of Cell Phone Vibration Motor Using a Hybrid Approach of ANFIS and Jaya,” Arabian Journal for Science and Engineering, pp. 1–16, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Ngoc Le Chau et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.