Research Article  Open Access
Research on Preference Polyhedron Model Based Evolutionary Multiobjective Optimization Method for Multilink Transmission Mechanism Conceptual Design
Abstract
To make the optimal design of the multilink transmission mechanism applied in mechanical press, the intelligent optimization techniques are explored in this paper. A preference polyhedron model and new domination relationships evaluation methodology are proposed for the purpose of reaching balance among kinematic performance, dynamic performance, and other performances of the multilink transmission mechanism during the conceptual design phase. Based on the traditional evaluation index of single target of multicriteria design optimization, the robust metrics of the mechanism system and preference metrics of decisionmaker are taken into consideration in this preference polyhedron model and reflected by geometrical characteristic of the model. At last, two optimized multilink transmission mechanisms are designed based on the proposed preference polyhedron model with different evolutionary algorithms, and the result verifies the validity of the proposed optimization method.
1. Introduction
To improve work efficiency of mechanical press and acquire specific kinematic and dynamic output of slider, multilink transmission mechanisms are applied to replace the traditional cranklink mechanisms. The determination of structure parameters is the most important link in the design process of the mechanism, which determines the ultimate characteristics and performances of the product to a great extent, especially in the conceptual design phase. The conceptual development of the product accomplishes almost 70 percent of the performances and characteristics while occupying only 1 percent of the whole lifecycle cost as shown in Figure 1. For this reason, this research focuses on the optimal design problems of the multilink transmission mechanism in the conceptual design phase.
The design optimization of multilink transmission mechanism is competitive not only for kinematic and dynamic performance, but also for manufacturability, serviceability, and overall lifecycle cost. Considering the conflicting objectives, as well as the highly complex search space and constraints, a rigorous and quantitative multidisciplinary design methodology and evaluation standard of design scheme are needed for solving such multiobjective optimization problem (MOOP). The evolutionary algorithms (EAs) could provide efficient solutions to the above problems [1, 2]. In such optimization problems, the design objectives describe all the product functions as well as the constraints under which these functions should be realized. Instead of a single optimal objective, multitudinous and competing objectives of actual engineering problems often provide a set of compromise solutions, which is generally called Pareto front. As the basis of decisionmaking, the Pareto front can offer valuable tradeoff information between objectives.
The feasibility of applying evolutionary algorithms to the solutions of multiobjective engineering optimization problems has been explored by many previous researches.
Marler and Arora [3] presented a survey of current continuous nonlinear multiobjective optimization (MDO) concepts and method. It was found that no single method was superior, while the selection of a specific method is depending on the type of information provided in the problem, the user’s preference, the solution requirements, and the availability of software.
Saravanan et al. [4] realized the multicriteria design optimization of robot grippers by using evolutionary algorithms. Three different robot gripper configurations with five competing objectives, nine constraints, and seven design variables were optimized by using intelligent techniques. What is more, the paper compared the computing time and the quality of the solution obtained from MOGA, NSGAII, MODE, and the proposed optimization procedures.
Castillo et al. [5] described the use of a genetic algorithm (GA) for the problem of offline pointtopoint autonomous mobile robot path planning based on two criteria, length and difficulty, and showed that GAs were effective tools for solving such pathplanning problems.
However, the researches mentioned above are based on traditional ranking and domination method. This research aims to extend the evaluation process of Pareto front solutions with the incorporation of robust design and concept of multiobjective risk decision. A preference polyhedron model is proposed to provide further dominant relationship information [6]. Thus, the Pareto front can be further improved and more valuable optimization results can be obtained. The rest of this paper consists of the following—Section 2 reviews some literatures. Section 3 introduces the proposed methodology and the key technologies for multilink transmission mechanism design. Section 4 describes the establishment of optimization model. And the experimental results are depicted in Section 5. Finally, the conclusions and further work are summarized in Section 6.
2. Basic Theory
This section will introduce the main concepts and previous research works related to this project, covering topics like basic evolutionary algorithm theory, multiobjective optimization, existing Pareto optimal front, and design optimization of multilink transmission mechanism. Moreover, the basic conception of preference polyhedron model will also be introduced briefly in this section.
2.1. Evolutionary Algorithms
Evolutionary algorithms are random exploring optimization algorithm based on the idea of the biological evolutionary. Different from traditional optimization algorithms, evolutionary algorithms have no strict requirements for the problem. After a repetitive loop or a series of generations, it can find the fittest individual or individuals to solve the particular problem. Thus, evolutionary algorithms are widely used to solve complex optimization problems [7, 8].
Genetic algorithm is a typical example of evolutionary algorithms, which regards each solution to the particular optimization problem as an individual in the evolution process of large population. The fitness of each individual is determined by a given fitness function which evaluates the level of aptitude that particular individual solves the given optimization problem, while each generation in the evolution process will create a new set of individuals through genetic operators: crossover and mutation operation. New child individuals produced by the above operations will then be selected by a selection method and finally reinserted into the population by a replacement method. Such process will be repeated generationally at a user defined number. While solving real engineering design problems, for example, to observe the optimal structure parameters of a complex mechanical system, the individual characteristics always can be coded into a finite set of such design parameters. These parameters make up the chromosome with genetic information that represents the realworld structure of the individual, which in this case is a solution to the particular optimization problem. The process can be concluded as shown in Figure 2.
2.2. Multiobjective Optimization Problems
In a real complex engineering optimization problem, lots of various performance factors or characteristic factors should be taken into consideration while such factors are always difficult to be quantified and competing with each other [9–11]. As a result, in such multiobjective optimization problems, it is difficult to reach an optimal solution with respect to all of the objectives evaluated individually. Creating a singleobjective function or converting the objectives into restrictions imposed on the optimization problem is one of the conventional ways [12–15]. For example, the evaluation function of a specific problem can be expressed as minimizing or maximizing the fitness evaluation function :where is weight factor, reflecting the contribution of each subobjective to the final evaluation function, and is the value of th subobjective.
However, the above method always overlooks the stability of the optimization result in the optimum iterative procedure. Figure 3 shows a typical example about deterministic optimal design point and robust optimal design point. As a result of the existence of error due to uncertainty factors (such as manufacturing and fixing error, external disturbance, and fatigue wear), the real product performances should be considered as an interval. Although the deterministic optimal design point provides the best performance, the performance fluctuation of the deterministic optimal design point may be larger than the other design points when the product performance is considered within the same and specific interval range. Sometimes the performance fluctuation due to error caused by uncertainty factors may be beyond the limit of design objective, so it is necessary to sacrifice some performances to get a more reliable and steady designing scheme.
2.3. Preference Polyhedron Model
Considering the interaction of different objectives in the multiobjective optimization problems, it is necessary to evaluate the importance of multiple objectives and determine the degree to which objectives need to be modified [16–18]. Human elements must be regarded as a key component in any successful system design methodology. As mentioned above, many existing methods for solving such multiobjective optimization problems were developed, and the human elements always took part in the iterative optimization process based on the concept of fuzzy mathematics and theory of probability [19, 20]. To better evaluate the optimal Pareto front, extended evaluation standards used to make further reflection to the quality of the alternative have been defined, for example, graphical illustration. In order to select the most preferred alternative, Miettinen researched different graphical visualization tools used to support and help the decisionmaker in understanding similarities and differences between the alternatives and graphical illustration [21]. Moreover, Trinkaus and Hanne represented a novel multicriteria decision support system (MCDSS) to solve multicriteria decision problems [22]. Visualization methods are both used in the above two researches for better description and analysis of the alternatives. In this paper, a preference polyhedron model is proposed and used to describe the alternatives, as shown in Figure 4. Under the basis of the description of the alternatives, the alternatives could be filtered by using the proposed preference polyhedron model. And the most preferred alternative satisfying the preference of decisionmakers will be automatically obtained.
The geometric size and shape parameters of the polygon models in Figure 4 reflect different evaluation information including quantized value of different design index and the stability of whole performance. For example, the length parameter of the full line in Figure 4 can reflect the value of the th indicator , and then it is easy to get the global evaluation function . However, the evaluation function of a specific problem should not include only evaluation function ; the error function reflecting the stability of performance should also be taken into consideration and can be defined as follows:
Moreover, the above polygon model can offer more geometric size and shape information which can be used to build extra evaluation function to improve the Pareto front of optimal results when the evaluation function . In contrast with the proposed methods in [23–25] all centered around the concept of Pareto optimality and the Pareto optimal set, in this paper, the preference of decisionmaker and the risk assessment of design scheme based on the area attribute of polygon model are used to define related dominant relationship.
2.4. Design Optimization Problem of Multilink Transmission Mechanism
The design optimization problem of multilink transmission mechanism is a typical multiobjective optimization problem, where various factors such as the kinematic performance, the dynamic performance, the lifecycle cost, manufacturability, and serviceability should be considered. This means that the evolutionary should optimize possible structure parameters based on at least such following criteria: kinematic performance, dynamic performance, manufacturing costs, and structural stability. In addition, the design parameters of multilink transmission mechanism include the length of each link and the angle parameters of component placement, which will ultimately influence the performance of the mechanism. The design parameters are usually defined as a design vector which contains all variables that form the design space. The design vector of this transmission mechanism can be expressed as follows:
The objective of the design of multilink transmission mechanism can be represented as a vector . It can be divided into different types, and each of them includes several subobjectives:
It is the most important thing to determine an optimal parameter set which can make the objective vector the balance of state where all disciplines can complement each other. However, interactions and conflictions among the above elements increase the difficulty of optimal solution of this complex engineering problem. The multiobjective optimization algorithm offered an effective solution for such problem. The typical framework of MOOD has been shown in Figure 5.
3. Optimization Model for Multilink Transmission Mechanism Design
This section briefly introduces the input and output parameters in the conceptual design process of multilink transmission mechanism. Moreover, the preference polyhedron model and related dominant relationship will also be introduced in detail in this section. By the iteration of the proposed optimization model, new multiobjective optimal solution based on risk assessment and robust design will ultimately be adopted. Figure 6 shows a comparison between the traditional cranklink mechanism and multilink transmission mechanism, which are both driven by a motor on the point and ultimately turn the rotational motion of motor to the straight reciprocating motion of the slider by the connecting rods.
3.1. Design Parameters and Constraint Equations
This simplified model of multilink transmission mechanism was derived from a real multilink mechanical press, so the design variables depend upon the original value of real design parameters. The design vector and objective vector can be expressed as follows:
The basic kinematic and dynamic performances and other performances concerned in this optimization are all defined by the elements of the vector and closely associated with the value of the design parameter in vector . In order to meet some geometric dimensions and motion requirements, some constraint equations must be satisfied during the optimization process. In addition, some elements in objective vector should be transformed into constraints to simplify the optimization model. These constraint equations are as follows:where are ideal design index and are reasonable deviation.
So, the optimization model is simplified to determine the optimal set of vectors to maximize or minimize the objective under the above constraints, where the vector is simplified as follows:
As a result of the error of the manufacturing and assembling, the design parameters always deviate from the ideal optimal design value and bring about the fluctuation of the output performance. The evaluation function of such performance fluctuation is defined in the following way:where in (9) reflects the performance fluctuation and reflects the impact of such fluctuation on the whole evaluation.
3.2. Preference Polyhedron Model and Related Dominant Relationship
The above multiobjective optimization problems were ultimately expressed by the design vector expressed by (5), the objective vector expressed by (6), the evaluation function expressed by (9), and some constraints. And multiple objectives are combined into scalar objective via weight vector, which can be determined randomly or adaptively through direct assignment, eigenvector method, empty method, and minimal information method. However, if the objective function is simply weighted and added to produce a single fitness, the value of the ultimate evaluation function may be dominated by the objective with large range and small weight. To avoid such condition, the first step to establish the optimization model is to normalize each objective. The normalizing condition can be defined in the following way: where in (10) is the feasible domain of design vector and is the lower limit of the value of objective function. In this way, all objective functions can be evaluated in the same metric. The preference polyhedron model can be created as shown in Figure 7.
The practical engineering problems always have many upper and lower limits, and the inside and outside red lines of the polyhedron in Figure 7 reflect such limits. The length of radial lines sent out from the center of the polyhedron reflects the value of object function. The parameter is determined by the angle between adjacent objective functions. Based on the definition of the confidence interval, the overall and individual confidence interval is expressed by the area of polyhedrons and . The confidence level of each individual in the iterative optimization process can be defined as follows:Equation (11) can be described by the following:where the above equation should meet the following constraint equations:
When the value of is out of limit just as shown in Figure 7, a penalty function is brought into the confidence function to reduce or increase the confidence level of the individual. The penalty functions are divided into two kinds of situations.
Definition 1. As shown in Figure 8, when and , the penalty function is defined as follows:In this definition, is penalty factor which is related to the importance level of the subobjective, the parameter is the angle between and and the parameter is the angle between and . What is more, based on the geometric relationship, the above parameters meet the following constraint equations:
Definition 2. When and , the penalty function is defined as follows:The confidence function can be ultimately revised as the following form:
In this way, the individual with high confidence level will be easier to be adopted to compare with the lower ones. Based on the definitions above, two individuals in the Pareto front (as shown in Figure 9) can be further compared by the confidence level in the following ways.
It is assumed that and , respectively, indicate the value of the th subobjective for the design vector , while and indicate the th subobjective. Based on the definition above, the domination relationship can be evaluated by the confidence level (the reliability of dominating ) and defined as follows:while in (18), if and , (18) can be corrected as follows:
Apparently, it can be easily proved that , so that the fact that dominate when can be true. However, it is necessary to increase the value of the righthand side of the inequality function to decrease the risk of decisionmaking.
4. The Establishment of Optimization Model
Different metal materials show different malleability under different working conditions including working speed and working temperature. In order to obtain high forming quality products with less energy consumption, the end effector of the transmission mechanism must satisfy some specific kinematic and dynamic performance requirements. More specifically, the slider of press should maintain constant working velocity and acceleration during the working process. In addition, the transmission mechanism should also have a good effect on the reinforcement to consume less energy. Based on the original structure design parameters and performance parameters, this section will discuss the establishment of the multiobjective optimization model.
4.1. The Initial Design Information
The main design objectives in this research are demonstrated in Table 1.

The initial values of input and output parameters for the multilink transmission mechanisms in Figure 6 are proposed in Tables 2 and 3. The initial values of structural design parameters are demonstrated in Table 2. Under the premise of meeting the basic design objectives, the real corresponding output performance parameters are demonstrated in Table 3.


4.2. The Design Information of Optimization Model
In order to confirm the value range of design parameters, correlation analysis is implemented between input and output parameters. Figure 10 is the correlation graph, which conveys the sensitivity information of output performance parameters to different design parameters.
From Figure 10, where the blue bars mean that the design variables are positively related to the output performance while the red bars indicate the negative correlations, it is obvious that the design parameter is the most sensitive design parameter for the output parameters , , and , which should be set at a large range of value under specific constraints. Thus, the preference of design parameters is determined as shown in Table 4.

4.3. The Information of Performance Stability
It is supposed that the error of design parameters in the processing and installation is ±5 mm and ±2° and follows the standard distribution. The effect caused by such errors is shown in Table 5 by random combination tests.

4.4. The Optimization Flow of the Model
The basic process of this optimization is just as shown in Figure 11; the most outstanding characteristic is the selection mechanism in contrast with the traditional optimization process which is based on the genetic algorithm. In the selection mechanism, individuals based on the multiobjectives weighted value created by the genetic operators will be ranked firstly. And, then, the first individuals will be selected to go into the next phase. After that, the above individuals will be ranked based on the level of confidence once again and finally the first individuals will be selected to be the ultimate generation.
The optimization process based on the evolutionary algorithms can be deterministic as you want by either changing the number of iterations or other termination criteria. This experiment adopts NSGAII multiobjective algorithms as the basic framework and takes the number of iterations as the terminal condition. Relative evolutionary parameters are set as Table 6 indicated.

It must be said that the initial values of output parameters are taken as the base value. And new evaluation parameters are obtained by the iteration where the comparison of dominant relationships among individuals is implemented by the preference polyhedron model based on (11) to (17) in Section 3.2.
5. Experimental Results
The information about new evaluation parameters can be obtained, after rationalizing and quantifying the product design problem by selecting appropriate objective functions, design variables, parameters, and constraints. Figures 12, 13, 14, and 15 indicate the history data of the parameters , , , and separately. From these figures, it can be shown that, with the increase of the number of iterations, each history data of parameters presents convergent tendency. What is more, the convergent tendency of Figures 12, 14, and 15 is more apparent than that of Figure 13. The reason for this is because the parameter is more sensitive in the range of design parameter selection. But it can be seen from this data that the output performance is improved apparently, which can be also indicated from Table 7. The output performance improvement can be explained by the percentage in Table 7, where the initial value was defined as 1.

6. Conclusion
This paper explored the use of intelligent optimization techniques to obtain optimum design of a multilink transmission mechanism. In contrast with traditional evolutionary algorithms, this paper adopted a polygon model into the iterative optimization process to describe the domination relationships of the individuals on the Pareto front. In addition, robust design optimization was used in the optimization process. From the experiment results, the output performance of the multilink transmission mechanism was improved apparently. Interesting future work will include a more detailed study of the polygon model on the description of the individuals on the Pareto front.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (no. 51505218) and the Fundamental Research Funds for the Central Universities (no. NS2015050). The supports are gratefully acknowledged.
References
 K. Deb, MultiObjective Optimization Using Evolutionary Algorithms, WileyInterscience, New York, NY, USA, 2001. View at: MathSciNet
 C. Coello, D. Veldhuizen, D. Van, and G. Lamont, Evolutionary Algorithms for Solving MultiObjective Problems, Genetic Algorithms and Evolutionary Computation, Kluwer Academic, Norwell, Mass, USA, 2002.
 R. T. Marler and J. S. Arora, “Survey of multiobjective optimization methods for engineering,” Structural & Multidisciplinary Optimization, vol. 26, no. 6, pp. 369–395, 2004. View at: Publisher Site  Google Scholar
 R. Saravanan, S. Ramabalan, N. G. R. Ebenezer, and C. Dharmaraja, “Evolutionary multi criteria design optimization of robot grippers,” Applied Soft Computing Journal, vol. 9, no. 1, pp. 159–172, 2009. View at: Publisher Site  Google Scholar
 O. Castillo, L. Trujillo, and P. Melin, “Multiple objective genetic algorithms for pathplanning optimization in autonomous mobile robots,” Soft Computing, vol. 11, no. 3, pp. 269–279, 2007. View at: Publisher Site  Google Scholar
 D. W. Gong, N. N. Qin, and X. Y. Sun, “Evolutionary optimization algorithm for multiobjective optimization problems with interval parameters,” in Proceedings of 5th IEEE International Conference on BioInspired Computing: Theories and Applications, pp. 411–420, Changsha, China, September 2010. View at: Google Scholar
 M. A. Abido, “Multiobjective evolutionary algorithms for electric power dispatch problem,” IEEE Transactions on Evolutionary Computation, vol. 10, no. 3, pp. 315–329, 2006. View at: Publisher Site  Google Scholar
 K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGAII,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182–197, 2002. View at: Publisher Site  Google Scholar
 M. Emmerich, N. Beume, and B. Naujoks, “An EMO algorithm using the hypervolume measure as selection criterion,” in Evolutionary MultiCriterion Optimization: Third International Conference, EMO 2005, Guanajuato, Mexico, March 9–11, 2005. Proceedings, vol. 3410 of Lecture Notes in Computer Science, pp. 62–76, Springer, Berlin, Germany, 2005. View at: Publisher Site  Google Scholar
 S. Gunawan, A. FarhangMehr, and S. Azarm, “Multilevel multiobjective genetic algorithm using entropy to preserve diversity,” in Evolutionary MultiCriterion Optimization: Second International Conference, EMO 2003, Faro, Portugal, April 8–11, 2003. Proceedings, vol. 2632 of Lecture Notes in Computer Science, pp. 148–161, Springer, Berlin, Germany, 2003. View at: Publisher Site  Google Scholar
 R. Battiti and A. Passerini, “Braincomputer evolutionary multiobjective optimization: a genetic algorithm adapting to the decision maker,” IEEE Transactions on Evolutionary Computation, vol. 14, no. 5, pp. 671–687, 2010. View at: Publisher Site  Google Scholar
 J. Branke, K. Deb, K. Miettinen, and R. Slowinski, Multiobjective OptimizationInteractive and Evolutionary Approaches, Springer, Heidelberg, Germany, 2008.
 K. Deb, A. Sinha, P. J. Korhonen, and J. Wallenius, “An interactive evolutionary multiobjective optimization method based on progressively approximated value functions,” IEEE Transactions on Evolutionary Computation, vol. 14, no. 5, pp. 723–739, 2010. View at: Publisher Site  Google Scholar
 H. Eskandari, C. D. Geiger, and R. Bird, “Handling uncertainty in evolutionary multiobjective optimization: SPGA,” in Proceedings of the IEEE Congress on Evolutionary Computation (CEC '07), pp. 4130–4137, IEEE, Singapore, September 2007. View at: Publisher Site  Google Scholar
 J. W. Fowler, E. S. Gel, M. M. Köksalan, P. Korhonen, J. L. Marquis, and J. Wallenius, “Interactive evolutionary multiobjective optimization for quasiconcave preference functions,” European Journal of Operational Research, vol. 206, no. 2, pp. 417–425, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. Rachmawati and D. Srinivasan, “Incorporating the notion of relative importance of objectives in evolutionary multiobjective optimization,” IEEE Transactions on Evolutionary Computation, vol. 14, no. 4, pp. 530–546, 2010. View at: Publisher Site  Google Scholar
 L. B. Said, S. Bechikh, and K. Ghédira, “The rdominance: a new dominance relation for interactive evolutionary multicriteria decision making,” IEEE Transactions on Evolutionary Computation, vol. 14, no. 5, pp. 801–818, 2010. View at: Publisher Site  Google Scholar
 J. D. Schaffer, “Multiple objective optimization with vector evaluated genetic algorithms,” in Proceedings of the 1st International Conference on Genetic Algorithms, pp. 93–100, Hillsdale, NJ, USA, 1985. View at: Google Scholar
 A. Sinha, P. Korhonen, J. Wallenius, and K. Deb, “An interactive evolutionary multiobjective optimization method based on polyhedral cones,” in Learning and Intelligent Optimization, C. Blum and R. Battiti, Eds., vol. 6073 of Lecture Notes in Computer Science, pp. 318–332, 2010. View at: Publisher Site  Google Scholar
 J. Sun, D. Gong, and X. Sun, “Solving interval multiobjective optimization problems using evolutionary algorithms with preference polyhedron,” in Proceedings of the 13th Annual Conference on Genetic and Evolutionary Computation (GECCO '11), pp. 729–736, ACM, July 2011. View at: Publisher Site  Google Scholar
 K. Miettinen, “Survey of methods to visualize alternatives in multiple criteria decision making problems,” OR Spectrum, vol. 36, no. 1, pp. 3–37, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 H. L. Trinkaus and T. Hanne, “KnowCube: a visual and interactive support for multicriteria decision making,” Computers and Operations Research, vol. 32, no. 5, pp. 1289–1309, 2005. View at: Publisher Site  Google Scholar
 C. M. Fonseca and P. J. Fleming, “Genetic algorithms for multiobjective optimization: formulation discussion and generalization,” in Proceedings of the 5th International Conference on Genetic Algorithms, pp. 416–423, 1993. View at: Google Scholar
 D. E. Goldberg, Genetic Algorithms in Search, Optimization, and Machine Learning, Compact MRI for TB of the Distal Radius, 1989.
 N. Srinivas and K. Deb, “Multiobjective optimization using nondominated sorting in genetic algorithms,” Evolutionary Computation, vol. 2, no. 3, pp. 221–248, 1994. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Haihua Zhu 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.