Computer-aided design tool for typical flexible
mechanisms synthesis by means of evolutionary
algorithms

Mohamed Amine Ben Abdallah, Imed Khemili, Nizar Aifaoui, Med Amine Laribi

► To cite this version:

Mohamed Amine Ben Abdallah, Imed Khemili, Nizar Aifaoui, Med Amine Laribi. Computer-aided design tool for typical flexible mechanisms synthesis by means of evolutionary algorithms. Robotica, 2024, 42 (4), pp.1172-1211. 10.1017/s0263574724000171 . hal-04717134

HAL Id: hal-04717134

https://hal.science/hal-04717134v1

Submitted on 1 Oct 2024

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

# Computer-aided design tool for typical flexible mechanisms synthesis by means of evolutionary algorithms

Mohamed Amine Ben Abdallah¹, Imed Khemili², Nizar Aifaoui³ and Med Amine Laribi⁴

¹Department of Engineering and Mathematics, College of Business and Technology, Sheffield Hallam University, Sheffield, UK, ²Laboratory of Mechanics of Sousse, National Engineering School of Sousse, University of Sousse, Sousse, Tunisia, ³Laboratory of Mechanical Engineering, National Engineering School of Monastir, University of Monastir, Monastir, Tunisia, and ⁴Département Génie Mécanique et systèmes complexes, Institut Pprime, University of Poitiers, Poitiers, France

**Corresponding author:** Mohamed Amine Ben Abdallah; Email: M.Ben-abdallah@shu.ac.uk

**Received:** 21 October 2022; **Revised:** 29 December 2023; **Accepted:** 3 January 2024; **First published online:** 4 March 2024

**Keywords:** dynamic synthesis; metaheuristics; flexible mechanism; computer-aided design for multibody system synthesis tool
"CADMS" optimization

## Abstract

Accurate prediction for mechanisms' dynamic responses has always been a challenging task for designers. For modeling easiness purposes, mechanisms' synthesis and optimization have been mostly limited to rigid systems, making consequently the designer unable to vow that the manufactured mechanism satisfies the target responses. To address this limitation, flexible mechanism synthesis is aimed in this work. Two benchmark mechanisms being the core of myriad mechanical devices are of scope, mainly, the flexible slider-crank and the four-bar. In addition to the mechanism dimensions, materials properties have been embedded in the synthesis problem. Two responses are of interest for the slider-crank mechanism, the slider velocity, and the midpoint axial displacement for the flexible connecting rod. Whereas five responses have been compiled for the four-bar mechanism synthesis. A comparative analysis of seven optimization techniques to solve the synthesis problem for both mechanisms has been performed. Subsequently, an executable computer-aided design tool for mechanisms synthesis has been developed under MATLAB®. Numerical outcomes emphasize the limits of a single-response-based synthesis for a flexible mechanism. It has been proven that combining different responses alleviates possible error and fulfill high-accuracy requirement.

## 1. Introduction

Nowadays, accurate trajectories and dynamic responses of multibody systems play a key role in their performance. Moreover, mechanical systems are invading a broad range of applications, that is, medical, welding, and manufacturing robots. Thus, they became a crucial part of production systems. Consequently, ceaseless improvement and development of these systems to deliver more accurate responses emerge as an emergency. Nevertheless, the development is facing many challenges and limits deterring to meet the level of expectations mainly because of modeling limitations. In order to satisfy an already defined trajectory, designers proceed with path, motion, or function generation synthesis. Despite that mechanisms are not modeled as flexible in most synthesis approach, realistically this assumption results in errors between the real and predicted responses. Once a mechanism is manufactured, its prescribed trajectory and dynamics responses are different to the predicted ones using the synthesis approach. Therefore, compliant modeling should be deployed for the mechanism's synthesis. Although synthesis analysis is limited to the mechanism end-effector trajectory, other responses should be considered for flexible mechanisms. Three different paradigms for formulating synthesis problems
have been presented in the literature, mainly, path generation, motion generation, and function generation synthesis.

Rigid mechanisms synthesis based on path generation has been the scope of several works. Shaoping Bai *et al.* [1] proposed a synthesis approach for a four-bar mechanism based on its coupler path. Sang Min Han *et al.* [2] have dealt with a topology optimization based on Fourier descriptors for path generation synthesis. The proposed approach has been tested through a set of illustrative examples. Yixin Shao *et al.* [3] have proposed a robotized gain rehabilitation system optimization, subsuming slider-crank and seven-bar mechanisms. The synthesis has been subject to dual-objective paths. The genetic algorithm (GA) has been deployed to optimize the system’s design parameters subject to constraints. A comparative analysis of the GA, the differential evolution (DE), and the imperialist competitive algorithm (ICA) for a rigid four-bar mechanism synthesis has been investigated in Saeed Ibrahimi *et al.* [4]. Volkan Parlaktaş *et al.* [5] have optimized the design parameters for a geared four-bar mechanism subsuming collinear input and output shafts. Sahand Hadizadeh Kafash *et al.* [6] have treated a four-bar mechanism synthesis based on its generated path. By means of the circular proximity function, the error separating the obtained and target paths for each of its constitutive points has been evaluated. G. Ganesan *et al.* [7] proposed a synthesis approach to overcome the filleted rectangular path generation challenges for an adjustable four-bar mechanism.

Another possible paradigm of the mechanisms synthesis relies on the motion generation. Jianwei Sun *et al.* [8] treated a planar four-bar mechanism synthesis based on the motion generation. Jean-François Collard *et al.* [9] have dealt with extensible links mechanisms synthesis. The extensible links have been modeled by means of springs. Two optimization strategies have been deployed to solve the synthesis problem. The first minimizes each variable deviation during one cycle of the mechanism motion. Whereas the parameters have been simultaneously involved in the second strategy. The proposed approach efficiency has been investigated for two benchmark mechanisms synthesis, a four-bar, and a six-bar with extensible links. A multiphase motion generation synthesis for the four-bar mechanism has been the focus of Venkatesh Venkataramanujam *et al.* [10].

In addition to path and motion generation synthesis, mechanisms synthesis can be carried out by means of function generation. A spatial spherical four-bar mechanism has been the interest of Rasim Alizade *et al.* [11]. By means of six independent design variables, the Quaternion algebra has been deployed with the Chebychev approximation to provide optimal design variables. It has been proven that the Chebychev approximation error owns the best performance besides interpolation approximation, least square approximation, and fitting error extremums.

Although the major interest has been devoted to rigid mechanism, some works investigated compliant mechanism synthesis. Zhen Luo *et al.* [12] have treated a topology optimization for a compliant micro-inverter. To this end, a finite element approach, avoiding the Courant-Friedrichs-Lewy issues, has been developed. Three performance criteria have been considered for the topology optimization, mainly the geometrical advantage (GA), mechanical efficiency (ME), and mechanical advantage (MA). An enhanced synthesis approach combining both the dimensional and pose errors for a robot arm has been discussed in ref. [13]. Numerical outcomes confirm the approach robustness to solve the synthesis problem for several scenarios.

Deepak *et al.* [14] dealt with a compliant mechanism topology synthesis. A comparative analysis for a set of compliant mechanisms formulation, that is, inverter and gripper, has been presented. The mechanism’s efficiency, characteristic stiffness, and artificial spring formulations have been the key performance indicators for the numerical outcomes. A continuum 16-degree-of-freedom snake robot has been the interest of Abdelkhalick Mohammad *et al.* [15]. An online algorithm has been developed allowing to follow the snake path and adjust it whenever needed. To this end, the coiled and uncoiled configurations of the robot have been considered. Another application of continuum robot has been treated in Anzhu work [16]. Different configurations for a contact-aided compliant mechanisms have been targeted in the numerical simulations. Experimental tests have been carried out to approve the numerical simulation outcomes. Chikhaoui *et al.* [17] have focused on a design optimization of continuum robot. Both the design and control aspects of the dual-arm continuum robot have been optimized.
Soft medical robot for cardiovascular application has been the scope of Liu Wang *et al.* [18]. The GA has been deployed to solve the design and optimization problem. Numerical outcomes corroborate well with the experimental data. The particle swarm optimization (PSO) has been deployed to find out optimal parameters subject to selected guides. Compliant mechanisms, taking into account their flexibility, have been treated in Nishiwaki *et al.* work [19]. Based on the mutual energy concept, the flexible body has been modeled as linear elastic. Two cost function have been considered for multi-objective optimization problem. The first cost function has taken into account the kinematic aspect of the problem, whereas the second modeled the structural aspect. An illustrative example for a compliant gripper has been used to confirm the proposed approach efficiency.

Thanks to their outstanding performance to solve complex optimization problems, metaheuristics techniques have been deployed in complex mechanical systems optimization such as two fingers gripper [20], belt mechanism [21], unmanned vehicle [22], self-reconfigurable robot [23], parallel manipulator [24], delta parallel robot [25], biped robot [26], chamfered rectangle [27], multi-point haptic device [28], designing a bracket of a vehicle [29–32], robot gripper optimization [33], design of gear train [29], car crashworthiness [34–35], clutch diaphragm [36], suspension arm [37–38], and parallelogram synthesis [39]. All the references are compiled and classified based on the synthesis paradigm in Appendix Table AI.

In this work, the focus is concentrated on flexible mechanism synthesis. For the sake of more compliant modeling, the end-effector dynamic responses have been considered as target responses instead of its path. Moreover, the axial displacement of the flexible body has been involved in the synthesis problem to reduce emerging errors from considering such characteristic. Two benchmark flexible mechanisms, slider-crank and four-bar have been used for validation. The synthesis involves in first instance separately the slider velocity, and the axial displacement taking place in the flexible connecting rod. Subsequently, the previous responses, in addition to the slider acceleration, have been simultaneously combined into the same cost function. For the four-bar mechanism synthesis, five responses have been gathered, namely the crank, the flexible coupler, and the flexible follower midpoint paths in addition to the axial displacements of the flexible coupler and follower. A straightforward tool (Computer-Aided Design for Multibody systems Synthesis “CADMS”) has been created to perform the synthesis of the mechanism of interest using one or several optimization techniques among the GA, simulating annealing (SA), DE, PSO, ICA, Artificial Bee Colony (ABC), and the Ant Colony (AC). The main contributions of this work are as follows: (1) propose a synthesis approach for mechanisms synthesis considering the flexible parts, (2) consider both the dimensional and material characteristics as parameters in the synthesis approach, and (3) provide a straightforward executable tool “CADMS” under MATLAB software for flexible mechanism synthesis offering a wide range of selection among a set of metaheuristic techniques.

The paper is structured as follows: Section 2 deals with the dynamic modeling and the synthesis problem formulation of the two flexible mechanisms of scope. The different optimization algorithms used in this paper are discussed in Section 3. Section 4 denotes an insight into the proposed CADMS tool. Section 5 focuses on the main obtained results. The most important conclusions have been summarized in Section 6.

## 2. Mathematical formulation

To accomplish the objective of this work, an approach based on two modules is designed. The first is responsible on the dynamic modeling of flexible mechanism, whereas the second deals with the optimization process as depicted in Appendix Figure A2. They communicate together to ensure an exchange of data, that is, dynamic responses or optimal parameters.

### 2.1. Dynamics of a flexible mechanism

Dynamics modeling of flexible multibody systems has been of interest in tremendous works. Equation of motion are often a combination of algebraic equations defining the system’s holonomic constraints and differential equations. They involve the time derivative for the mechanism-independent parameters.
<figure><img src="image_2.png" /><figcaption>Figure 1. The general coordinates for a flexible beam.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Figure 2. The flexible connecting rod in both rigid and deformed configurations.</figcaption></figure>

Different approaches for solving differential-algebraic equations of flexible and rigid mechanisms have been discussed in refs. [40–44]. In this work, the dyad finite elements method has been used [42] for the dynamic analysis of flexible slider-crank and four-bar mechanisms. It is suitable for flexible mechanism modeling with affordable computational burden. The flexible bodies of the mechanism are considered to be homogeneous and isotropic. It is noteworthy that this study is limited only to purely linear elastic modeling of all flexible bodies with small deformation. Hence, all the assumptions of infinitesimal strain theory are fulfilled. Among the possible existing models in the literature for flexible beams, that is, Euler, Kirchhoff, and Cosset, the authors used Euler formulation with the dyad finite elements methods discussed in ref. [42]. The ability to provide accurate results of the chosen method has been approved through experimental investigation presented in the aforementioned reference.

The synthesis process is governed by an interaction between two parts. This process will be implemented in a graphical user interface tool discussed later in Section 4. While the optimization technique has been selected by the user (arrow 1), an initial random population is generated. Thereby, as depicted in Appendix Figure A2, a connection will be established with the dynamic model as pointed out by arrow number 2. Subsequently, the dynamic model will simulate their responses (arrow 3). In order to assess the accuracy of the dynamic responses, they are sent back to the optimization module (narrow 5). If the maximum iteration number is not yet reached, the algorithm will go back to arrow 1. Else, it will follow arrows 6 and 7 to stop the process. More details will be provided later in Section 4 for the implementation of this process into a graphical user interface. In what follows, the mathematical modeling for the Euler beam based on the Dyad finite elements [42] is detailed.

The six general coordinates describing the flexible beam are detailed in Fig. 1. Mainly, two dyads constitute the flexible beam. Figure 2 depicts the six general coordinates $u_1, u_2, u_3, u_4, u_5, u_6$, for the flexible connecting rod. The absolute acceleration for a beam element can be written as shown in Eq. (1):

$$ \ddot{u}_a = \ddot{u}_r + \ddot{u} + a_n + a_c + a_t \quad (1) $$

where $\ddot{u}_a$ is the absolute acceleration, $\ddot{u}_r$ is the rigid element acceleration, $\ddot{u}$ is the generalized relative acceleration, $a_n$ is the normal acceleration, $a_c$ is Coriolis, and $a_t$ is tangential accelerations.
For this work, due to the small strain formulation assumption, only two components of the Eq. (1) are considered, the rigid element acceleration and the generalized relative acceleration. The new equation for the beam element acceleration yields

$$ \ddot{u}_a = \ddot{u}_r + \ddot{u} \quad (2) $$

where $\ddot{u}_a$ is the absolute acceleration, $\ddot{u}_r$ is the rigid element acceleration, and $\ddot{u}$ is the generalized relative acceleration.

Since the connecting rod is considered flexible, an axial displacement should take place along the x-axis and expressed as follows:

$$ v(\bar{x}, t) = \phi_1(\bar{x})u_1(t) + \phi_4(\bar{x})u_4(t) \quad (3) $$

where $u_1(t)$ and $u_4(t)$ are the general coordinates detailed in Fig. 1. $\phi_1(\bar{x})$ and $\phi_4(\bar{x})$ are the shape functions. The y-axis transversal displacement is written in Eq. (4) as:

$$ w(\bar{x}, t) = \phi_2(\bar{x})u_2(t) + \phi_3(\bar{x})u_3(t) + \phi_5(\bar{x})u_5(t) + \phi_6(\bar{x})u_6(t) \quad (4) $$

where $u_2(t)$, $u_3(t)$, $u_5(t)$, and $u_6(t)$ are the general coordinates detailed in Fig. 1. $\phi_2(\bar{x})$, $\phi_3(\bar{x})$, $\phi_5(\bar{x})$, and $\phi_6(\bar{x})$ are the shape functions.

Several conditions should be fulfilled by the shape function used to model the flexible bodies of the mechanism. The shape functions should ensure the continuity within the element as well as in the boundaries related to other elements. In addition, the shape function should consider the stress and strain during the link deformation. Moreover, the boundaries compatibility for the displacement must be satisfied. To meet these requirements, the shape functions have been chosen as:

$$ \phi_1(\bar{x}) = 1 - \frac{\bar{x}}{L}; \quad \phi_4(\bar{x}) = \frac{\bar{x}}{L}; \quad \phi_2(\bar{x}) = 3\left(\frac{L-\bar{x}}{L}\right)^2 - 2\left(\frac{L-\bar{x}}{L}\right)^3; \quad \phi_3(\bar{x}) = \bar{x}\left(\frac{L-\bar{x}}{L}\right)^2; \\ \phi_5(\bar{x}) = 3\left(\frac{\bar{x}}{L}\right)^2 - 2\left(\frac{\bar{x}}{L}\right)^3; \quad \phi_6(\bar{x}) = -L(L-\bar{x})\left(\frac{\bar{x}}{L}\right)^2 $$

where L is the total length of the flexible beam.

The motion equation, based on the Lagrange principle for a flexible beam, yields

$$ \frac{d}{dt} \left( \frac{\partial T}{\partial \dot{u}_i} \right) - \frac{\partial T}{\partial u_i} + \frac{\partial U}{\partial u_i} = \bar{Q}_i \quad (5) $$

where T, U and $\bar{Q}_i$ are, respectively, the kinetic energy, the deformation energy, and the external forces applied to the beam element. $u_i$ and $\dot{u}_i$ are respectively the general coordinate and the velocity of the i-th coordinate.

Consequently, the kinetic energy matrix form can be expressed as:

$$ T = \frac{1}{2} \dot{u}^T \bar{m} \dot{u} \quad (6) $$

where $\dot{u}$, $\dot{u}^T$ are the first derivative with respect of time for the local coordinates vector and its transpose, respectively. $\bar{m}$ is the elementary inertia matrix in the local coordinates of the beam.

The reduced matrix formulation for the deformation energy yields

$$ U = \frac{1}{2} u^T \bar{k} u \quad (7) $$

where u and $u^T$ are the general coordinates vector and its transpose. $\bar{k}$ is the elementary stiffness matrix in the local coordinates of the beam.

Following [42], a single dyad motion equation, in the flexible beam, regarding its local frame is written as:

$$ \bar{m} \ddot{u}_a + \bar{k} u = \bar{Q} \quad (8) $$

where $\bar{Q}$ is the general forces vector expressed in the local frame.
The local and global coordinates are related by means of the following equation:

$$u = RU \tag{9}$$

The rotation matrix is expressed as shown in Eq. 10:

$$R = \begin{bmatrix} \lambda & \mu & 0 & 0 & 0 & 0 \\ -\mu & \lambda & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & \lambda & \mu & 0 \\ 0 & 0 & 0 & -\mu & \lambda & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 \end{bmatrix} \tag{10}$$

where $\lambda = \cos \theta_3$, $\mu = \sin \theta_3$ and $U$ is the general coordinates vector expressed in the global frame.

The global kinetic energy regarding the global frame is written as shown in Eq. 11:

$$T = \frac{1}{2} \dot{U}^T m \dot{U} \tag{11}$$

where

$$m = R^T \bar{m} R \tag{12}$$

$m$ is the mass matrix in the global frame.

In a similar way, the flexible beam deformation energy yields

$$V = \frac{1}{2} U^T k U \tag{13}$$

wherein

$$k = R^T \bar{k} R \tag{14}$$

$k$ is the stiffness matrix in the global coordinates.

Basically, two dyads constitute the flexible beam instilled in the mechanism. Hereby, both kinetic and deformation energies for the whole mechanism yield

$$T = \frac{1}{2} \{\dot{U}_i\}^T M \{\dot{U}_i\} \tag{15}$$

where $M$ is the assembling mass matrix.

$$U = \frac{1}{2} \{U_i\}^T K \{U_i\} \tag{16}$$

where $K$ is the assembling stiffness matrix.

The global motion equation relative to the flexible mechanism relying on the dyad finite element method, considering the beam structural damping, yields

$$M \ddot{U} + C \dot{U} + K U = -M \ddot{U}_r \tag{17}$$

where $\ddot{U}_r$ is the acceleration vector.

The damping matrix has been defined in ref. [45] as:

$$C_k = 2 \zeta_k \omega_k \theta_k \theta_k^T \tag{18}$$

where $\theta_k = M \phi_k$, $\zeta_k$ is the damping coefficient, $\omega_k$ is the $k^{th}$ eigenfrequency, and $\phi_k$ is the $k^{th}$ eigenmode.

Based on the selected eigenmodes numbers, $n$, the damping matrix yields

$$C = \sum_{k=1}^{n} C_k \tag{19}$$

The following Eq. (20) allows to solve the free vibration systems as first step. Consequently, the eigenvalues for the system can be defined. The next step for the resolution will be to solve the motion equation
<table><caption>**Table I.** The inputs parameters.</caption><thead><tr><th></th><th>Flexible slider-crank</th><th>Flexible four-bar</th></tr></thead><tbody><tr><td>The crank angular velocity (θ<sub>2</sub>)</td><td>26.18 rad/s</td><td>26.18 rad/s</td></tr><tr><td>The crank length (l<sub>2</sub>)</td><td>50 mm</td><td>100 mm</td></tr><tr><td>The flexible connecting rod/coupler length (l<sub>3</sub>)</td><td>350 mm</td><td>250 mm</td></tr><tr><td>The flexible follower length (l<sub>4</sub>)</td><td>—</td><td>250 mm</td></tr><tr><td>The frame length (l<sub>1</sub>)</td><td>—</td><td>285 mm</td></tr><tr><td>The flexible connecting rod density</td><td>2660 kg/m<sup>3</sup></td><td>2660 kg/m<sup>3</sup></td></tr><tr><td>The flexible connecting rod Young's modulus</td><td>71,000 MPa</td><td>71,000 MPa</td></tr><tr><td>The flexible connecting rod cross section</td><td>0.36 10<sup>-4</sup> m<sup>2</sup></td><td>0.36 10<sup>-4</sup> m<sup>2</sup></td></tr><tr><td>The crank cross section</td><td>4.68 10<sup>-4</sup> m<sup>2</sup></td><td>4.68 10<sup>-4</sup> m<sup>2</sup></td></tr><tr><td>The slider mass</td><td>0.245 kg</td><td>—</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** The flexible mechanisms: (a) slider crank and (b) four-bar.</figcaption></figure>

for the whole system. It should be highlighted that resolution of the motion equation is impossible without the first step:

$$K\phi_n = \omega_n^2 M\phi_k \quad (20)$$

where $\phi_n$ and $\omega_n$ are, respectively, the eigenmode and the eigenfrequency.

To be able to decouple Eq. (17), modal coordinates should be used as written in Eq. (21):

$$U = \Phi\eta, \dot{U} = \Phi\dot{\eta}, \ddot{U} = \Phi\ddot{\eta} \quad (21)$$

where $\Phi$ and $\eta$ are, respectively, the modal matrix and amplitude.

Finally, Eq. (17) yields

$$\Phi^T M \Phi \ddot{\eta} + \Phi^T C \Phi \dot{\eta} + \Phi^T K \Phi \eta = -\Phi^T M \ddot{U}_r \quad (22)$$

As detailed previously, the target responses to satisfy are defined following the designer convenience (arrow 6 Appendix Figure A2). Given responses for the two mechanisms of scope, the flexible slider-crank and four-bar have been considered as a case of study to investigate the proposed synthesis approach efficiency (Fig. 3).

The characteristic parameters for the target mechanism responses are detailed in Table I.

## 2.2. The synthesis problem formulation

As pointed out in Section 1, the main novelty of this paper is flexible mechanism synthesis based on a single or multiple dynamic responses, that is, velocity, acceleration of the end-effector, or axial displacement for a flexible body. The difference between the proposed approach, compared to classic ones, is that it does not rely on the end-effector path but the dynamic responses instead. Moreover, the flexible body’s axial displacement is involved in the synthesis analysis. Hereby, the cost function ought to evaluate different response types with different scales. The synthesis problem is extended from its most commonly treated case dealing with dimensions, to involve additional specifications such as material density and young modulus. Mechanism synthesis has been often formulated as an optimization problem. Thus, acting on a set of parameters, the error between the target and obtained responses is minimized.
<figure><img src="image_5.png" /><figcaption>**Figure 4.** The error evaluation.</figcaption></figure>

The cost function should fit and accommodate all the responses. To do so, the error in each time increment is divided by the module of the difference between the extremum obtained in the response curve as witnessed in Fig. 4. Consequently, the error will be dimensionless and scaled independently from the response types allowing to compile different responses as linear combination. Based on the same core cost function, the synthesis problem has been formulated for the two mechanisms of interest in this work.

### 2.2.1 The mathematical formulation of the flexible slider-crank mechanism synthesis

Six design variables have been selected for the mechanism synthesis combining both dimensional and characteristic parameters. The idea behind picking similar parameters is to optimize and control flexible bodies stiffness in order to satisfy a desired axial displacement threshold. The design variables are $l_2$ for the crank length and $l_3$, $\rho$, $E$, $h$, $a$, respectively, the flexible connecting rod length's, material density's, Young modulus's, high's as well as width's.

The core objective function (Eq. 23) is formulated to be generic and could apply to any response type, whether velocity, acceleration, or axial displacement. It provides normalized and dimensionless values. Thereby, combining different response types as linear combination into the same cost function could be henceforth, possible:

$$ fitness = \text{Minimize}(f(R(X))) = \text{Minimize} \left( f(R(l_2, l_3, \rho, E, h, a)) \right) \quad (23) $$

where

$$ f(R(l_2, l_3, \rho, E, h, a)) = \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left[ \frac{R_i(l_2, l_3, \rho, E, h, a) - R_i^*(l_2^*, l_3^*, \rho^*, E^*, h^*, a^*)}{\text{Max}(R_i(l_2, l_3, \rho, E, h, a)) - \text{Min}(R_i(l_2, l_3, \rho, E, h, a))} \right]^2} \quad (24) $$

Subject to:

$$ l_2 \cos \theta_2 + l_3 \cos \theta_3 - x_3 = 0 $$

$$ l_2 \sin \theta_2 + l_3 \sin \theta_3 = 0 $$

$$ L_{b2} \le l_2 \le U_{b2} $$

$$ L_{b3} \le l_3 \le U_{b3} $$

$$ l_2 < l_3 $$

$$ L_{b4} \le \rho \le U_{b4} $$

$$ L_{b5} \le E \le U_{b5} $$

$$ L_{b6} \le a \le U_{b6} $$

$$ L_{b7} \le h \le U_{b7} $$

(25)
where R is the considered response and $X = (l_2, l_3, \rho, E, h, a)$ representing the design variables. The parameter $R_i(l_2, l_3, \rho, E, h, a)$ represents the mechanism given response using the design variable vector $X$, whereas $R_i^*(l_2^*, l_3^*, \rho^*, E^*, h^*, a^*)$ represents the target response defined by the user. The subscript $i$ for both the generated and target responses is related to the number of prescribed points ($N$) equal to 121 for an average of one acquisition point each 0.004 s along 0.48 s. The simulation time is chosen to cover two crank's revolution.

Substituting the letter R in the core function (Eq. 24) yields to equations Eq. (26), (27), and (28) dealing respectively with the acceleration error $f(A(l_2, l_3, \rho, E, h, a))$, velocity error $f(V(l_2, l_3, \rho, E, h, a))$, and axial displacement error $f(D(l_2, l_3, \rho, E, h, a))$. The slider velocity, acceleration, as well as the connecting rod midpoint axial displacement are strongly dependent on the set of following parameters $(l_2, l_3, \rho, E, h, a)$:

$$f(A(l_2, l_3, \rho, E, h, a)) = \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left[ \frac{A_i(l_2, l_3, \rho, E, h, a) - A_i^*(l_2^*, l_3^*, \rho^*, E^*, h^*, a^*)}{\mathrm{Max}(A_i(l_2, l_3, \rho, E, h, a)) - \mathrm{Min}(A_i(l_2, l_3, \rho, E, h, a))} \right]^2} \quad (26)$$

$$f(V(l_2, l_3, \rho, E, h, a)) = \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left[ \frac{V_i(l_2, l_3, \rho, E, h, a) - V_i^*(l_2^*, l_3^*, \rho^*, E^*, h^*, a^*)}{\mathrm{Max}(V_i(l_2, l_3, \rho, E, h, a)) - \mathrm{Min}(V_i(l_2, l_3, \rho, E, h, a))} \right]^2} \quad (27)$$

$$f(D(l_2, l_3, \rho, E, h, a)) = \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left[ \frac{D_i(l_2, l_3, \rho, E, h, a) - D_i^*(l_2^*, l_3^*, \rho^*, E^*, h^*, a^*)}{\mathrm{Max}(D_i(l_2, l_3, \rho, E, h, a)) - \mathrm{Min}(D_i(l_2, l_3, \rho, E, h, a))} \right]^2} \quad (28)$$

The velocity (Eq. 27)) and acceleration (Eq. (26)) considered in this work are relative to the end-effector component. However, in Eq. (28), the considered axial displacement $D_i$ is taking place at the connecting rod midpoint. The combined optimization compiles the three mechanism responses into the same cost function. A weight coefficient has been allocated to each response as recommended in ref. [46]. The combined error is governed by Eq. (29) as follows:

$$\begin{aligned} f(D, V, A) = & \alpha f(D(l_2, l_3, \rho, E, h, a)) + \beta f(V(l_2, l_3, \rho, E, h, a)) + \gamma f(A(l_2, l_3, \rho, E, h, a)) \\ = & \alpha \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left[ \frac{D_i(l_2, l_3, \rho, E, h, a) - D_i^*(l_2^*, l_3^*, \rho^*, E^*, h^*, a^*)}{\mathrm{Max}(D(l_2, l_3, \rho, E, h, a)) - \mathrm{Min}(D(l_2, l_3, \rho, E, h, a))} \right]^2} \\ & + \beta \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left[ \frac{V_i(l_2, l_3, \rho, E, h, a) - V_i^*(l_2^*, l_3^*, \rho^*, E^*, h^*, a^*)}{\mathrm{Max}(V(l_2, l_3, \rho, E, h, a)) - \mathrm{Min}(V(l_2, l_3, \rho, E, h, a))} \right]^2} \\ & + \gamma \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left[ \frac{A_i(l_2, l_3, \rho, E, h, a) - A_i^*(l_2^*, l_3^*, \rho^*, E^*, h^*, a^*)}{\mathrm{Max}(A(l_2, l_3, \rho, E, h, a)) - \mathrm{Min}(A(l_2, l_3, \rho, E, h, a))} \right]^2} \end{aligned} \quad (29)$$

where $\alpha = 0.4$, $\beta = \gamma = 0.3$ based on the parametric study carried out in ref. [46].

### 2.2.2 The mathematical formulation for the flexible four-bar mechanism synthesis

Since the four-bar mechanism subsumes two flexible links, its synthesis should be treated differently. Usually, for rigid four-bar mechanisms, the focus is devoted to the coupler path. However, when flexible parts are considered, this approach could not be accurate. Consequently, the crank, flexible coupler, as well as the flexible follower midpoint drawn paths, are considered. In addition, to keep the axial displacement during the synthesis process in acceptable thresholds, both axial displacements of the coupler and follower have been considered toward the mechanism synthesis. By means of the proposed
core cost function in Eq. (24), combining all the five responses is possible providing a scalable and dimensionless error:

$$f(l_1, l_2, l_3, l_4) = \sqrt{ \frac{1}{N} \sum_{i=1}^{N} \left[ \left( \frac{x_{cri}(l_1, l_2, l_3, l_4) - x_{cri}^*(l_1^*, l_2^*, l_3^*, l_4^*)}{\mathrm{Max}(x_{cri}(l_1, l_2, l_3, l_4)) - \mathrm{Min}(x_{cri}(l_1, l_2, l_3, l_4))} \right)^2 + \left( \frac{y_{cri}(l_1, l_2, l_3, l_4) - y_{cri}^*(l_1^*, l_2^*, l_3^*, l_4^*)}{\mathrm{Max}(y_{cri}(l_1, l_2, l_3, l_4)) - \mathrm{Min}(y_{cri}(l_1, l_2, l_3, l_4))} \right)^2 } + \sqrt{ \frac{1}{N} \sum_{i=1}^{N} \left[ \left( \frac{x_{coi}(l_1, l_2, l_3, l_4) - x_{coi}^*(l_1^*, l_2^*, l_3^*, l_4^*)}{\mathrm{Max}(x_{coi}(l_1, l_2, l_3, l_4)) - \mathrm{Min}(x_{coi}(l_1, l_2, l_3, l_4))} \right)^2 + \left( \frac{y_{coi}(l_1, l_2, l_3, l_4) - y_{coi}^*(l_1^*, l_2^*, l_3^*, l_4^*)}{\mathrm{Max}(y_{coi}(l_1, l_2, l_3, l_4)) - \mathrm{Min}(y_{coi}(l_1, l_2, l_3, l_4))} \right)^2 } + \sqrt{ \frac{1}{N} \sum_{i=1}^{N} \left[ \left( \frac{x_{foi}(l_1, l_2, l_3, l_4) - x_{foi}^*(l_1^*, l_2^*, l_3^*, l_4^*)}{\mathrm{Max}(x_{foi}(l_1, l_2, l_3, l_4)) - \mathrm{Min}(x_{foi}(l_1, l_2, l_3, l_4))} \right)^2 + \left( \frac{y_{foi}(l_1, l_2, l_3, l_4) - y_{foi}^*(l_1^*, l_2^*, l_3^*, l_4^*)}{\mathrm{Max}(y_{foi}(l_1, l_2, l_3, l_4)) - \mathrm{Min}(y_{foi}(l_1, l_2, l_3, l_4))} \right)^2 } + \sqrt{ \frac{1}{N} \sum_{i=1}^{N} \left[ \frac{D_{coxi}(l_1, l_2, l_3, l_4) - D_{coxi}^*(l_1^*, l_2^*, l_3^*, l_4^*)}{\mathrm{Max}(D_{coxi}(l_1, l_2, l_3, l_4)) - \mathrm{Min}(D_{coxi}(l_1, l_2, l_3, l_4))} \right]^2 + \sqrt{ \frac{1}{N} \sum_{i=1}^{N} \left[ \frac{D_{foxi}(l_1, l_2, l_3, l_4) - D_{foxi}^*(l_1^*, l_2^*, l_3^*, l_4^*)}{\mathrm{Max}(D_{foxi}(l_1, l_2, l_3, l_4)) - \mathrm{Min}(D_{foxi}(l_1, l_2, l_3, l_4))} \right]^2 } } \quad (30)$$

$$\begin{array}{l}
\text{Subject to:} \\
l_2 + l_4 < l_1 + l_3 \\
l_2 + l_3 < l_1 + l_4 \\
l_1 + l_2 < l_3 + l_4 \\
L_{b1} \le l_1 \le U_{b1} \\
L_{b2} \le l_2 \le U_{b2} \\
L_{b3} \le l_3 \le U_{b3} \\
L_{b4} \le l_4 \le U_{b4}
\end{array}
\qquad (31)$$

Errors following the X and Y axis for the prescribed paths of the crank, coupler, and follower beside the target ones as well as axial displacements for the coupler and follower have been assessed. Errors have been measured among 121 points of the drawn responses. The cost function is formulated Eq. (30), where $l_1, l_2, l_3, l_4$ and $l_1^*, l_2^*, l_3^*, l_4^*$ are respectively the generated and the reference design variables, where $l_1$ is the frame length, $l_2$ is the crank length, $l_3$ is the flexible coupler length, and $l_4$ is the flexible follower length. The parameters $x_{foi}y_{foi}x_{coi}y_{coi}x_{cri}y_{cri}$ represent the coordinates along x and y axis for the follower, coupler, and crank midpoints, respectively. All variables labeled with the subscript (*) such as $x^*_{foi}, y^*_{foi}, x^*_{coi}, y^*_{coi}, x^*_{cri}$, and $y^*_{cri}$ characterize the target response coordinates, respectively, for the follower, coupler, and crank. $D_{cox}$ and $D_{fox}$ are respectively the coupler and follower axial displacements; meanwhile, $D^*_{cox}$ and $D^*_{fox}$ are the ones relative to the target response. $L_b$ and $U_b$ represent the lower and upper bounds for the search intervals of the design variables, respectively.

### 3. Resolution algorithms

Metaheuristics optimization techniques seem the most prominent to solve similar synthesis problems owing to their complexity. Therefore, several optimization techniques have been deployed in this work, that is, the GA [47], SA [48], DE [49], PSO [50], Imperialist Colony Algorithm [51], Artificial Bee Colony [52], and the Ant Colony [53]. The different algorithms pseudo-codes and how they have been adopted to the actual synthesis problem have been discussed detailed in ref. [54]. Search intervals for the different parameters of the flexible slider-crank mechanism, namely, the crank length ($l_c$), the flexible connecting rod length ($l_s$), material density ($\rho$), Young Modulus ($E$) as well as the connecting rod high
<table><caption>**Table II.** The different parameter search intervals for the slider-crank mechanism.</caption><thead><tr><th></th><th>l<sub>2</sub> (mm)</th><th>l<sub>3</sub> (mm)</th><th>E (kPa)</th><th>ρ (kg/m<sup>3</sup>)</th><th>h (mm)</th><th>a (mm)</th></tr></thead><tbody><tr><td>Search interval [L<sub>b</sub>, U<sub>b</sub>]</td><td>[20;80]</td><td>[200;500]</td><td>[6.5 10<sup>10</sup>;7.5 10<sup>10</sup>]</td><td>[2000;3000]</td><td>[1;2]</td><td>[20;30]</td></tr></tbody></table>

<table><caption>**Table III.** The different parameter search intervals for the four-bar mechanism.</caption><thead><tr><th></th><th>l<sub>1</sub> (mm)</th><th>l<sub>2</sub> (mm)</th><th>l<sub>3</sub> (mm)</th><th>l<sub>4</sub> (mm)</th></tr></thead><tbody><tr><td>Search interval [L<sub>b</sub>, U<sub>b</sub>]</td><td>[100;300]</td><td>[70;200]</td><td>[200;300]</td><td>[200;300]</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 5.** The Computer-Aided Design for Multibody Synthesis (CADMS) graphical interface.</figcaption></figure>

(h) and width (a), are summarized in Table II. L<sub>b</sub>, U<sub>b</sub> represent the lower and upper bounds for the search intervals of the design variables, respectively.

Four parameters have been of scope for the flexible four-bar mechanism, namely the frame length (l<sub>1</sub>), the crank length (l<sub>2</sub>), the coupler length (l<sub>3</sub>), and the follower length (l<sub>4</sub>). The different parameters search intervals have been summarized in Table III.

## 4. Computer-Aided Design for Multibody systems Synthesis

The graphical interface is created to provide the user with a straightforward tool subsuming all necessary commands. The welcoming window gives the user the possibility to select the mechanism of interest as illustrated in Figure 5. The second window contains the run button, in addition to four main sections as shown in Figure 6 (response type synthesis, design variables, optimization techniques, and target responses). The slider acceleration or velocity, the connecting rod axial displacement, and the combined synthesis constitute the cornerstone of the slider-crank mechanism menu. On the other hand, the crank, the flexible follower, and the flexible coupler paths or a combination of all of them constitute the synthesis menu for the flexible four-bar mechanism.

Six design variables are of scope for the flexible slider-crank mechanism, mainly the crank (l<sub>2</sub>), the flexible connecting rod (l<sub>3</sub>) lengths, flexible connecting rod material density (ρ), and Young Modulus
<figure><img src="image_7.png" /><figcaption>**Figure 6.** The flexible slider-crank mechanism synthesis menu.</figcaption></figure>

(E) as well as the connecting rod height (h) and width (a) as exhibited in Fig. 5. However, regarding the four-bar, the frame (l1), crank (l2), flexible coupler (l3), and flexible follower (l4) lengths are optimized (Fig. 7).

Possible search intervals avoiding any constraint violation are the constituents of the design variables menu. The tool can automatically evaluate any possible incoherence of the selected search intervals. If so, an error message will be displayed to help the user adjusting them. The choice of the optimization technique is possible through the third part of this tool titled optimization techniques. The last step before running the simulation is to introduce the target responses. Depending on the selected or available response data, the auspicious fitness function, expressed in Eqs. (26), (27), (28), (29), and (30), will be selected.

Finally, thanks to the run button the user can start the synthesis process. Subsequently, a post processing window pops up allowing to pick the desired response to plot as shown in Fig. 8. The illustrated case in Fig. 8, is for a combined synthesis selecting all the optimization techniques.

The next section compiles numerical outcomes for the proposed synthesis approach implemented in the Computer-Aided Design for Multibody systems Synthesis (CADMS) tool.

## 5. Results and discussions

Towards this section, the proposed approach for flexible mechanisms synthesis based on dynamic responses or axial displacement of a flexible part is investigated. Two responses, mainly the slider velocity and flexible connecting rod midpoint axial displacement, have been involved separately for the flexible slider-crank mechanism synthesis (Eqs. (27), (28)). Subsequently, the previous responses, in addition to the slider acceleration (Eq. (26)), have been involved in a combined synthesis (Eq. (29)). Five responses are of interest for the flexible four-bar mechanism, namely the crank, the coupler, and the follower midpoint paths' as well as the coupler and follower axial displacements. The difference between
<figure><img src="image_8.png" /><figcaption>**Figure 7.** The four-bar mechanism synthesis menu.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 8.** The postprocessing window for the flexible slider-crank mechanism.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 9.** *The slider velocity based on the slider velocity synthesis.*</figcaption></figure>

the chosen responses for both mechanisms is strictly due to their industrial applications. Usually, for the slider-crank mechanism, the designer is interested in the end-effector dynamic responses. However, four-bar mechanisms are commonly used in complicated devices to satisfy a given path. It is noteworthy that considering the different component flexibility in the four-bar mechanism, the resulting path is pulled away from the target one. The obtained optimal solutions have been chosen based on statistical analysis for the obtained results over 100 runs of all the optimization techniques. A fitting for the outcome histograms has been made to approximate the distribution law and the key statistical metrics governing these results. The numerical simulations have been performed using an Intel I7 9th generation CPU computer.

**5.1. Flexible slider-crank mechanism synthesis**

**5.1.1. Synthesis based on the slider velocity**

As a first investigation for the proposed synthesis approach, the slider velocity has been chosen as a target response among the synthesis process. Statistical metrics confirm that the synthesis results corroborate well with a normal law distribution, wherein the mean $\mu$ and the standard deviation $\sigma$ are detailed in Appendix Table AII. The authors opted for the solution having the highest probability and the same or nearest mean $\mu$ of the normal law distribution fitting as depicted in Appendix Figure A1. The deployed optimization techniques' performances are compared according to this selection standard. It could be inferred that the Ant Colony (AC) technique outperforms the other ones proposing a set of design variables settling at a nil error. The DE method is less accurate with an error estimated of $1.8503 \times 10^{-17}$. Accordingly, this is reflected in results of Figs. 9 and 10, respectively, for the slider velocity and acceleration through an overlapping for target and optimal responses. Nevertheless, a proportional discrepancy to the errors in Appendix Table AII is witnessed for the other optimization techniques. Although its low accuracy, the GA presents an interesting CPU time to conclude the required number of iterations. The results in Fig. 11 reveal difficulty to satisfy the axial displacement threshold based on the velocity synthesis. In perfect agreement with Kim's work [55] for the first- and second-order derivative curves for mechanisms synthesis, satisfying the mechanism velocity results consequently in satisfying its acceleration.

It can be witnessed in Fig. 12 that the error evolution as function of iteration number contains some flat area. Within these areas, the error remains constant. Optimization algorithms have a set of steps allowing
<figure><img src="image_11.png" /><figcaption>**Figure 10.** The slider acceleration based on the slider velocity synthesis.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 11.** The axial displacement based on the slider velocity synthesis.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 12.** The error evolution based on the slider velocity synthesis.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>Figure 13. Graphical summary for the performance of the optimization techniques based on the slider velocity, (a) for the slider velocity, (b) for the slider acceleration, and (c) for the axial displacement.</figcaption></figure>

them to test different random combination of parameters from the search interval. If the algorithm could not manage to retrieve a better solution, it will keep the one from the previous iteration. Hence, this justifies the constant error for some iterations.

To emphasize the gap between the target and optimal responses retrieved by the optimization algorithms, a zoom focusing on specific area has been made as per Figs. 9 and 10. The responses of scope in this work are obtained by means of numerical integration of the motion equations. The minimum time step chosen in this work is 0.004 s. Consequently, any two successive points, separated by the time step, are linked with a straight line. For clarity purposes, it is useful to zoom beyond the time step interval. Therefore, linear areas can be visible in the zoom area.

It can be witnessed in Fig. 12 that the error evolution as function of the iteration number differs from one technique to another. As each of the used techniques is based on different algorithm, it should be expected that their performances will be different. This is in perfect agreement with previous works where different techniques have been used in optimization problems [4, 56–60]. The GA common process consists of randomly generating initial population from the search intervals, selection, crossover, and mutation. This same process is governing the DE where the mutation is taking place before the crossover. For the SA technique, the process starts with generating the initial population to create subsequently new solution based on the mutation and evaluation steps. The mutation process in this case helps the algorithm to assess each solution with its neighbourhood in order to retrieve the best one. This will result consequently in building new population and updating the ever-best solution. The ant colony, particle swarm, and artificial bee colony are swarm-based algorithms. They are based on the swarm collaboration in their search for food. By analogy, the best food site is considered as optimal solution for the optimization problem. The search is based on the swarm communication and collaboration to find the best food site. Therefore, the algorithms manage to find the optimal solution. As the communication differs from birds to bees or ants, the optimal solution differs from one technique to another explaining the different levels of error in Fig. 12.

Fig. 13 depicts the performance of the optimization techniques. Ideally, the heptagon should be tangent to the green circle where the performance of the optimization technique is considered excellent. However, as in Table IV, some of the techniques performances are limited or moderated, and the heptagon is tangent to either the red or amber circles for some techniques. Fig. 13 is a visual illustration of Table IV. An in-depth analysis for the results obtained in Appendix Table AII reveals that some of techniques can be more advantageous compared to others for the mechanism synthesis based on the slider velocity. Two main criteria are considered to assess the advantages and downsides of these techniques, mainly the calculation time to achieve the required number of iterations and the accuracy. The AC settle at a zero error for a CPU time of 3248 s. The DE manages to provide accurate results as well. However, the CPU time required is around 97% higher than the AC. It is worth considering the results of the GA
<table><caption>**Table IV.** Summary of the synthesis results based on the velocity synthesis.</caption><thead><tr><th rowspan="2">Optimization techniques</th><th colspan="3">Performance based on the slider velocity synthesis</th></tr><tr><th>Slider velocity</th><th>Slider acceleration</th><th>Axial displacement</th></tr></thead><tbody><tr><td>Genetic algorithm (GA)</td><td>Limited</td><td>Limited</td><td>Excellent</td></tr><tr><td>Simulating annealing (SA)</td><td>Moderated</td><td>Moderated</td><td>Limited</td></tr><tr><td>Differential evolution (DE)</td><td>Excellent</td><td>Excellent</td><td>Limited</td></tr><tr><td>Imperialist competitive algorithm (ICA)</td><td>Moderated</td><td>Moderated</td><td>Limited</td></tr><tr><td>Artificial bee colony (ABC)</td><td>Moderated</td><td>Moderated</td><td>Limited</td></tr><tr><td>Ant colony (AC)</td><td>Excellent</td><td>Excellent</td><td>Moderated</td></tr><tr><td>Particle warm optimization (PSO)</td><td>Excellent</td><td>Excellent</td><td>Moderated</td></tr></tbody></table>

<figure><img src="image_15.png" /><figcaption>**Figure 14.** The error evolution based on the flexible connecting rod midpoint axial displacement.</figcaption></figure>

as it requires as little as 30% of the CPU time for the AC. The low CPU time is the major advantage of the GA, making it very prominent for preliminary design stage synthesis.

To make the results more understandable and readable, the performance of each optimization technique, based on the responses involved in the synthesis process, is summarized in Table IV. The performance of the optimization techniques has been assessed with, limited, moderated, or excellent based on how accurate the responses of the mechanism. In the next section, the mechanism synthesis is performed based on the axial displacement taking place in the connecting rod midpoint.

### 5.1.2 The axial displacement synthesis

This section is devoted to the synthesis results based on axial displacement threshold for the flexible connecting rod midpoint. It can be witnessed in Fig. 14 that this synthesis-based response is a burdensome task to accomplish. The errors evolution plotted in Fig. 14 have the mean $\mu$ parameter value of the normal distribution shown in Appendix Figure A3 and Appendix Table AIII. Figure 15 for the slider velocity depicts a considerable discrepancy beside the one taking place for the slider velocity synthesis. As witnessed in Fig. 16, the AC solution responses almost overlap the target ones. It is interesting to highlight that the ABC technique exhibits a good performance, as illustrated in Fig. 17, wherein error worth $1.25 \times 10^3$ despite a CPU time of 12,723 s.

The proposed design variables fulfil only the axial displacement illustrated in Fig. 17. These results corroborate well with the slider velocity synthesis results discussed in Section 5.1.1 proving that except
<figure><img src="image_16.png" /><figcaption>**Figure 15.** The slider velocity based on the axial displacement synthesis.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 16.** The slider acceleration based on the axial displacement synthesis.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 17.** The axial displacement based on the axial displacement synthesis.</figcaption></figure>
<table><caption>**Table V.** Summary of the synthesis results based on the axial displacement synthesis.</caption><thead><tr><th rowspan="2">Optimization techniques</th><th colspan="3">Performance based on the axial displacement synthesis</th></tr><tr><th>Slider velocity</th><th>Slider acceleration</th><th>Axial displacement</th></tr></thead><tbody><tr><td>Genetic algorithm (GA)</td><td>Limited</td><td>Limited</td><td>Limited</td></tr><tr><td>Simulating annealing (SA)</td><td>Moderated</td><td>Moderated</td><td>Moderated</td></tr><tr><td>Differential evolution (DE)</td><td>Moderated</td><td>Moderated</td><td>Moderated</td></tr><tr><td>Imperialist competitive algorithm (ICA)</td><td>Moderated</td><td>Moderated</td><td>Moderated</td></tr><tr><td>Artificial bee colony (ABC)</td><td>Moderated</td><td>Moderated</td><td>Moderated</td></tr><tr><td>Ant colony (AC)</td><td>Limited</td><td>Limited</td><td>Moderated</td></tr><tr><td>Particle swarm optimization (PSO)</td><td>Limited</td><td>Limited</td><td>Moderated</td></tr></tbody></table>

<figure><img src="image_19.png" /><figcaption>**Figure 18.** Graphical summary for the performance of the optimization techniques based on the axial displacement, (a) for the slider velocity, (b) for the slider acceleration, and (c) for the axial displacement.</figcaption></figure>

for the involved response the other ones are dismissed. The performance of the different optimization techniques based on the axial displacement synthesis is assessed in Table V. Fig. 18 summarizes the results presented in Table V. The comparison of the different optimization techniques reveals that four out of the seven implemented techniques in this work have moderated performance. This is clearly seen through tangency of the heptagon to the amber circle for both the slider's velocity and acceleration. The three other techniques have limited performances for the slider's velocity and acceleration. For the axial displacement, the figure depicts almost a perfect heptagon included inside the moderated performance circle except from the GA edge where it is tangent to the red circle indicating a limited performance.

### 5.1.3 The combined error synthesis

Based on the outcomes of Sections 5.1.1 and 5.1.2, the designer is subject to several constraints. Whether to satisfy the target dynamic responses namely the slider velocity and acceleration or to perfectly fit the defined threshold for the midpoint axial displacement of the flexible connecting rod. In order to simultaneously satisfy the three responses of interest, they have been involved in a combined synthesis. A weight coefficient is allocated to each response error, $α = 0.4$, $β = 0.3$, and $γ = 0.3$, respectively, for midpoint axial displacement, the slider velocity, and slider acceleration errors, respectively.

Numerical outcomes rely on statistical analysis. To this end, histograms, for 100 runs of the different optimization techniques, have been plotted. As it can be seen in Appendix Figure A4, histograms dispersion could be appropriately modeled as a normal law distribution. Consequently, all the approved
<figure><img src="image_20.png" /><figcaption>**Figure 19.** *The error evolution based on the combined synthesis.*</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Figure 20.** *The slider velocity responses based on the combined synthesis.*</figcaption></figure>

solutions for the set of optimization techniques of interest represent the mean values of the normal distribution as shown in Appendix Table AIV. The error evolution plot for the set of implemented optimization techniques is depicted in Fig. 19. It could be inferred that in sharp contrast with the single response synthesis, the PSO wrests the lead from the AC. The relative error evolution curve for the PSO technique converges in about 240 iterations for a combined error worth $7.9712 \times 10^{-5}$.

As witnessed in Fig. 20 for the slider velocity responses, the discrepancy is considerably reduced compared to the midpoint axial displacement synthesis. Moreover, all the optimization techniques' responses seem almost overlapping the target ones for the slider acceleration as per Fig. 21. Figure 22 proves acceptable given responses' accuracy beside the target ones. The combined synthesis provides design variable parameters' dealing simultaneously with the three responses of scope concluding a perfect trade-off for them. Although their different levels of accuracy as shown in Table VI, the overall performance for the optimization techniques is acceptable.

The different optimization techniques compared in this work are based on different algorithms. Therefore, they perform with slightly different levels of accuracy. The two main criteria considered to assess the advantages and downsides of these techniques are the calculation time required and the accuracy. This will help the designer to select the most auspicious technique that could be used based on the design stage, that is, preliminary or final design stage, high or low accuracy application, etc. The analysis of results confirms that the DE, ICA, Ant Colony, and PSO have the advantage of providing
<table><caption>**Table VI.** Summary of the synthesis results based on the combined synthesis.</caption><thead><tr><th rowspan="2">Optimization techniques</th><th colspan="3">Performance based on the combined synthesis</th></tr><tr><th>Slider velocity</th><th>Slider acceleration</th><th>Axial displacement</th></tr></thead><tbody><tr><td>Genetic algorithm (GA)</td><td>Limited</td><td>Limited</td><td>Limited</td></tr><tr><td>Simulating annealing (SA)</td><td>Moderated</td><td>Moderated</td><td>Moderated</td></tr><tr><td>Differential evolution (DE)</td><td>Excellent</td><td>Excellent</td><td>Excellent</td></tr><tr><td>Imperialist competitive algorithm (ICA)</td><td>Excellent</td><td>Excellent</td><td>Excellent</td></tr><tr><td>Artificial bee colony (ABC)</td><td>Moderated</td><td>Moderated</td><td>Moderated</td></tr><tr><td>Ant colony (AC)</td><td>Excellent</td><td>Excellent</td><td>Excellent</td></tr><tr><td>Particle swarm optimization (PSO)</td><td>Excellent</td><td>Excellent</td><td>Excellent</td></tr></tbody></table>

<figure><img src="image_22.png" /><figcaption>**Figure 21.** The slider acceleration responses based on the combined synthesis.</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>**Figure 22.** The axial displacement responses based on the combined synthesis.</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>**Figure 23.** Graphical summary for the performance of the optimization techniques based on the combined synthesis, (a) for the slider velocity, (b) for the slider acceleration, and (c) for the axial displacement.</figcaption></figure>

outstandingly accurate results. The PSO is taking the lead being the most advantageous methods for the accuracy. However, its main downside is the calculation time needed to conclude the 250 iterations. Among the four aforementioned techniques, the ant colony perfectly compromise the accuracy and the calculation time.

The ABC and SA could be classified as second rank techniques for the combined synthesis. These two methods are incurring the drawback of a substantial calculation time. Despite its large error compared to other techniques, the GA manages to run the 250 iterations in shorter time. Thus, it could be interesting to perform preliminary design calculation using the GA to have an idea about the overall parameters of the mechanism. However, for applications requiring high accuracy, the authors would recommend the use of the PSO or the Aunt Colony. Fig. 23 is a graphical summary for the results discussed in Table VI.

## 5.2. Four-bar mechanism synthesis

The responses involved in the flexible four-bar mechanism synthesis are slightly different to the flexible slider-crank mechanism. Three paths are considered for the crank, the flexible coupler, and the flexible follower midpoints. In addition, the axial displacements of the flexible coupler and follower have been made part of the mechanism synthesis. It is important to incorporate the midpoint displacement of a flexible parts to overcome the classic synthesis approaches limits. Commonly, the synthesis is only limited to midpoints' paths. Using the approach proposed in this paper, the synthesis remains possible based on a single response. However, based on the results of Section 5.1.3, only the combined synthesis is of scope to satisfy all the target responses simultaneously.

As exhibited in Appendix Figure A5, the chosen solutions in Appendix Table AV perfectly fit with the normal distribution law mean $\mu$ obtained over 100 runs. Even though the obtained errors are spread over several solutions, it still in a very tiny interval as confirmed by the standard deviation $\sigma$ for all the optimization techniques. The simulation is set for two revolutions of the crank. Therefore, 121 acquisition points have been considered over a period of time equal to 0.48 s. Thus, the time step is 0.004 s. The analysis of Fig. 24 for the error evolution confirms that some of the deployed techniques are more burdened than others to converge. The DE proposed design variables responses for the crank midpoint paths are nearly overlapping the target response (Fig. 25). A more in-depth analysis of Fig. 26 reveals that the AC, SA, and DE proposed design variables gap to the target response is narrow. The DE proposed design variables' responses almost overlap the flexible coupler midpoint target path. The coupler represents an important component of the four-bar mechanism. Satisfying as maximum as possible the target response will have a considerable impact to reduce any eventual gap in case the coupler is connected to another mechanism.
<figure><img src="image_25.png" /><figcaption>**Figure 24.** *The error evolution.*</figcaption></figure>

<figure><img src="image_26.png" /><figcaption>**Figure 25.** *The crank midpoint path.*</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>**Figure 26.** *The flexible coupler midpoint path.*</figcaption></figure>

The third response involved in the synthesis process is the generated path of the flexible follower midpoint. The analysis of Fig. 27 confirms that the generated path of the follower midpoint is the most affected response beside the crank and the coupler midpoint prescribed trajectories. This is mainly due to the previous gaps for the crank and coupler responses resulting to enlarge the gap between the target
<figure><img src="image_28.png" /><figcaption>**Figure 27.** *The flexible follower midpoint path.*</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>**Figure 28.** *The flexible coupler midpoint axial displacement.*</figcaption></figure>

and obtained responses for the follower. Thereby, the witnessed discrepancy for the follower response is substantially larger. However, the SA technique exhibits a far response from the target one.

In addition to the three aforementioned responses, the mechanism is subject to given axial displacement responses. Furthermore, controlling the axial displacement contributes to alleviate the flexible bodies' impact on the mechanism response. Contrary to solid multibody systems, flexible bodies deviate the mechanism from its target response. Figures 28 and 29 represent the obtained results for the set of optimization techniques tested in this work. It can be observed that, apart from the ABC and GA, performances of other techniques remain acceptable. For the DE and AC techniques, the axial displacement responses almost overlap the target ones. Hence, the designer can consider the possibility of making the coupler or the follower midpoints' part of another mechanism to build more complicated mechanical devices. The performance of each optimization technique based on the different responses has been assessed in Table VII and graphically summarized in Fig. 30.

Comparing the different optimization techniques for the four-bar mechanism synthesis reveals that some techniques could be appealing. It is worth mentioning that being more advantageous depends imperatively on the application type and design stage. For applications wherein the accuracy is a major concern, the trio SA, DE, and Ant Colony optimization are taking the lead. Despite being the technique managing to deplete the error to the lowest values, the SA requires more resources for the calculation time. The Ant Colony requires only 15% of the CPU time, whereas the DE requires 47% of the CPU time required by the SA. The increase of accuracy for the SA is 6% compared to the AC and 3% to the
<figure><img src="image_30.png" /><figcaption>**Figure 29.** The flexible follower midpoint axial displacement.</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>**Figure 30.** Graphical summary for the performance of the optimization techniques based on the combined synthesis, (a) for the crank midpoint path, (b) for the flexible coupler midpoint path, (c) for the flexible follower midpoint path, (d) for the flexible coupler midpoint axial displacement, and (e) for the flexible follower midpoint axial displacement.</figcaption></figure>

DE. The PSO, ICA, and ABC performances' are relatively similar. Among the three aforementioned techniques, the ICA seems to be the most convenient as it offers the highest accuracy with the minimum CPU time. Even though the ABC technique has the advantage of providing design parameters with decent accuracy, its abundantly high CPU time is a major downside.
<table><caption>**Table VII.** Summary of the synthesis results for the four-bar mechanism.</caption><thead><tr><th rowspan="2">Optimization techniques</th><th colspan="5">Performance for the four-bar mechanism synthesis</th></tr><tr><th>The crank midpoint path</th><th>The flexible coupler midpoint path</th><th>The flexible follower midpoint path</th><th>The flexible coupler midpoint axial displacement</th><th>The flexible follower midpoint axial displacement</th></tr></thead><tbody><tr><td>Genetic algorithm (GA)</td><td>Excellent</td><td>Limited</td><td>Limited</td><td>Limited</td><td>Limited</td></tr><tr><td>Simulating annealing (SA)</td><td>Limited</td><td>Excellent</td><td>Excellent</td><td>Excellent</td><td>Excellent</td></tr><tr><td>Differential evolution (DE)</td><td>Excellent</td><td>Excellent</td><td>Excellent</td><td>Excellent</td><td>Excellent</td></tr><tr><td>Imperialist competitive algorithm (ICA)</td><td>Moderated</td><td>Excellent</td><td>Excellent</td><td>Excellent</td><td>Excellent</td></tr><tr><td>Artificial bee colony (ABC)</td><td>Limited</td><td>Limited</td><td>Limited</td><td>Moderated</td><td>Moderated</td></tr><tr><td>Ant colony (AC)</td><td>Moderated</td><td>Excellent</td><td>Excellent</td><td>Excellent</td><td>Excellent</td></tr><tr><td>Particle swarm optimization (PSO)</td><td>Moderated</td><td>Moderated</td><td>Excellent</td><td>Moderated</td><td>Moderated</td></tr></tbody></table>

It can be confirmed from Appendix Table AV that the error values are within the same range of values as in refs. [4, 58, 60, 61]. This confirms the suitability of the proposed cost function by the authors to provide significant and meaningful results for the mechanism synthesis based on dynamic responses.

## 6. Conclusion

This work denotes an insight into flexible mechanisms synthesis. A flexible slider-crank and four-bar mechanisms have been used as illustrative examples. Six design variables have been of scope for the flexible slider-crank mechanism encompassing dimensional and material characteristics parameters. Whereas four dimensional parameters have been of interest for the flexible four-bar mechanism. Seven optimization techniques have been used for the mechanism synthesis based on either a dynamic response or axial displacement for flexible parts. Some important conclusions about the performance of the optimization techniques can be drawn. The ant colony and the DE techniques outstandingly perform compared to the other techniques for a single synthesis-based response. It should be highlighted that for flexible mechanisms, a single synthesis-based response satisfies only the response involved in the synthesis process. For the combined synthesis, all the responses involved are satisfied simultaneously concluding a perfect trade-off between them. The proposed approach offers to the designer a set of possible solutions for the synthesis problem based on noncommon responses, namely dynamic responses or a flexible link axial displacement. Moreover, the material characteristics such as Young modulus and material density are parts of the parameters to optimize. It is interesting to extend this work to other mechanisms such as robot arms and exoskeleton. Including clearance into the mechanical joints alongside flexible bodies is also of scope in future works.

**Author contributions.** Mohamed Amine Ben Abdallah and Imed Khemili conceived and designed the study. Mohamed Amine Ben Abdallah conducted data gathering and analyses of the results and wrote the article. Imed Khemili, Nizar Aifaoui, and Med Amine Laribi reviewed the first draft of the manuscript.
**Financial support.** This research received no specific grant from any funding agency, commercial, or not-for-profit sectors.

**Competing interests.** The authors declare that there is no competing interests.

**References**

[1] S. Bai and J. Angeles, “Coupler-curve synthesis of four-bar linkages via a novel formulation,” *Mech. Mach. Theory* **94**, 177–187 (2015).

[2] S. M. Han, S. In Kim and Y. Y. Kim, “Topology optimization of planar linkage mechanisms for path generation without prescribed timing,” *Struct. Multidiscip. Optimiz.* **56**(3), 501–517 (2017).

[3] Y. Shao, Z. Xiang, H. Liu and L. Li, “Conceptual design and dimensional synthesis of cam-linkage mechanisms for gait rehabilitation,” *Mech. Mach. Theory* **104**, 31–42 (2016).

[4] S. Ebrahimi and P. Payvandy, “Efficient constrained synthesis of path generating four-bar mechanisms based on the heuristic optimization algorithms,” *Mech. Mach. Theory* **85**, 189–204 (2015).

[5] V. Parlaktaş, E. Söylemez and E. Tanik, “On the synthesis of a geared four-bar mechanism,” *Mech. Mach. Theory* **45**(8), 1142–1152 (2010).

[6] S. Hadizadeh Kafash and A. Nahvi, “Optimal synthesis of four-bar path generator linkages using Circular Proximity Function,” *Mech. Mach. Theory* **115**, 18–34 (2017).

[7] G. Ganesan and M. Sekar, “Optimal synthesis and kinematic analysis of adjustable four-bar linkages to generate filleted rectangular paths,” *Mech. Based Des. Struct. Mach.* **45**(3), 363–379 (2017).

[8] J. Sun, P. Wang, W. Liu and J. Chu, Non-integer-period motion generation of a planar four-bar mechanism using wavelet series,” *Mech. Mach. Theory* **121**,28–41 (2018).

[9] J. F. Collard, P. Duysinx and P. Fisette, “Optimal synthesis of planar mechanisms via an extensible-link approach,” *Struct. Multidiscip. Optimiz.* **42**(3), 403–415 (2010).

[10] V. Venkataramanujam and P. M. Larochelle, “Design and development of planar reconfigurable motion generators,” *Mech. Based Des. Struct. Mach.* **44**(4), 426–439 (2016).

[11] R. Alizade and E. Gezgin, “Synthesis of function generating spherical four bar mechanism for the six independent parameters,” *Mech. Mach. Theory* **46**(9), 1316–1326 (2011).

[12] Z. Luo, L. Tong, M. Y. Wang and S. Wang, “Shape and topology optimization of compliant mechanisms using a parameterization level set method,” *J. Comput. Phys.* **227**(1), 680–705 (2007).

[13] M. Russo, L. Raimondi, X. Dong, D. Axinte and J. Kell, Task-oriented optimal dimensional synthesis of robotic manipulators with limited mobility,” *Robot Comput. Integr. Manuf.* **69**, 102096 (2021).

[14] S. R. Deepak, M. Dinesh, D. K. Sahu and G. K. Ananthasuresh, “A comparative study of the formulations and benchmark problems for the topology optimization of compliant,” *J. Mech. Robot.* **1**(1), 011003 (2009).

[15] A. Mohammad, M. Russo, Y. Fang, X. Dong, D. Axinte and J. Kell, “An efficient follow-the-leader strategy for continuum robot navigation and coiling,” *IEEE Robot. Autom. Lett.* **6**(4), 7493–7500 (2021).

[16] A. Gao, J. Li, Y. Zhou, Z. Wang and H. Liu, “Modeling and task-oriented optimization of contact-aided continuum robots,” *IEEE/ASME Trans. Mechatron.* **25**(3), 1444–1455 (2020).

[17] M. T. Chikhaoui, J. Granna, J. Starke and J. Burgner-Kahrs, “Toward motion coordination control and design optimization for dual-arm concentric tube continuum robots,” *IEEE Robot. Autom. Lett.* **3**(3), 1793–1800 (2018).

[18] L. Wang, D. Zheng, P. Harker, A. B. Patel, C. F. Guo and X. Zhao, “Evolutionary design of magnetic soft continuum robots,” *Proc. Natl. Acad. Sci.* **118**(21), e2021922118 (2021).doi:10.1073/pnas.2021922118.

[19] S. Nishiwaki, M. I. Frecker, S. Min and N. Kikuchi, “Topology optimization of compliant mechanisms using the homogenization method,” *Int. J. Numer. Methods Eng.* **42**(3), 535–559 (1998).

[20] I. Gholaminezhad and A. Jamali, “A multi-objective differential evolution approach based on ε-elimination uniform-diversity for mechanism design,” *Struct. Multidiscip. Optimiz.* **52**(5), 861–877 (2015).

[21] W. Y. Lin and K. M. Hsiao, “A new differential evolution algorithm with a combined mutation strategy for optimum synthesis of path-generating four-bar mechanisms,” *Proc. Inst. Mech. Eng. C J. Mech. Eng. Sci.* **231**(14), 2690–2705 (2017).

[22] Y. Liu and R. Bucknall, “A survey of formation control and motion planning of multiple unmanned vehicles,” *Robotica* **36**(7), 1019–1047 (2018).

[23] T. Zhang, W. Zhang and M. M. Gupta, “An underactuated self-reconfigurable robot and the reconfiguration evolution,” *Mech. Mach. Theory* **124**, 248–258 (2018).

[24] S. Botello-Aceves, S. I. Valdez, H. M. Becerra and E. Hernandez, “Evaluating concurrent design approaches for a Delta parallel manipulator,” *Robotica* **36**(5), 697–714 (2018).

[25] J. Brinker, B. Corves and Y. Takeda, “Parallel robots based on motion/force transmission indices,” *Mech. Mach. Theory* **125**, 111–125 (2018).

[26] G. Gupta and A. Dutta, “Trajectory generation and step planning of a 12 DoF biped robot on uneven surface,” *Robotica* **36**(7), 945–970 (2018).

[27] G. Ajith Kumar, G. Ganesan and M. Sekar, “Near perfect path generation of corners chamfered rectangle and single synthesis cam-link mechanism to generate special-slot path,” *Mech. Based Des. Struct. Mach.* **46**(4), 483–498 (2018).

[28] R. Roberts and E. Rodriguez-Leal, “Kinematics and workspace-based dimensional optimization of a novel haptic device for assisted navigation,” *Mech. Based Des. Struct. Mach.* **44**(1-2), 43–57 (2016).

[29] A. R. Yildiz and M. U. Erdaş, “A new hybrid taguchi-salp swarm optimization algorithm for the robust design of real-world engineering problems,” *Mat. Test.* **63**(2), 157–162 (2021).
[30] B. S. Yıldız, V. Patel, N. Pholdee, S. M. Sait, S. Bureerat and A. R. Yıldız, "Conceptual comparison of the ecogeography-based algorithm, equilibrium algorithm, marine predators algorithm and slime mold algorithm for optimal product design," *Mat. Test.* **63**(4), 336-340 (2021).

[31] B. S. Yildiz, "Robust design of electric vehicle components using a new hybrid salp swarm algorithm and radial basis function-based approach," *Int. J. Vehicle Des.* **83**(1), 38-53 (2020).

[32] N. Panagant, N. Pholdee, S. Bureerat, K. Kaen, A. R. Yıldız and S. M. Sait, "Seagull optimization algorithm for solving real-world design optimization problems," *Mat. Test.* **62**(6), 640-644 (2020). doi: 10.3139/120.111529.

[33] B. S. Yıldız, N. Pholdee, S. Bureerat, A. R. Yıldız and S. M. Sait, "Robust design of a robot gripper mechanism using new hybrid grasshopper optimization algorithm," *Expert Syst.* **38**(3), e12666 (2021).

[34] S. Bureerat, S. M. Sait, C. M. Aye, N. Pholdee and A. R. Yıldız, "Multi-surrogate-assisted metaheuristics for crashworthiness optimisation," *Int. J. Vehicle Des.* **80**(2-4), 223-240 (2019).

[35] B. S. Yıldız, "Slime mould algorithm and kriging surrogate model-based approach for enhanced crashworthiness of electric vehicles," *Int. J. Vehicle Des.* **83**(1), 54-68 (2020). Available from: https://www.inderscienceonline.com/doi/abs/10.1504/IJVD.2020.114786

[36] B. S. Yıldız, N. Pholdee, N. Panagant, S. Bureerat, A. R. Yıldız and S. M. Sait, "A novel chaotic Henry gas solubility optimization algorithm for solving real-world engineering problems," *Eng. Comput.* **1**(2), 871-883 (2022).

[37] B. S. Yıldız, N. Pholdee, S. Bureerat, A. R. Yıldız and S. M. Sait, "Enhanced grasshopper optimization algorithm using elite opposition-based learning for solving real-world engineering problems," *Eng. Comput.* **38**(5), 4207-4219 (2022).

[38] A. B. S. Yıldız, N. Pholdee, S. Bureerat, A. R. Yıldız and S. M. Sait, "Sine-cosine optimization algorithm for the conceptual design of automobile components," *Mat. Test.* **62**(7), 744-748 (2020). doi: 10.3139/120.111541.

[39] S. I. Valdez Peña, E. Chávez-Conde, E. E. Hernandez and M. Ceccarelli, "Structure-control design of a mechatronic system with parallelogram mechanism using an estimation of distribution algorithm," *Mech. Based Des. Struc. Mach.* **44**(1-2), 58-71 (2016).

[40] A. A. Shabana, "Computational dynamics," *Encycl. Aerosp. Eng.* **1**, 1-545 (2010).

[41] A. A. Shabana, *Dynamics of Multibody Systems*. 4th edition (Cambridge University Press, England, 2013) pp. 233.

[42] L. Romdhane, H. Dhuibi and H. B. H. Salah, "Dynamic analysis of planar elastic mechanisms using the dyad method," *Proc. Inst. Mech. Eng. Part K J. Mult. Dyn.* **217**(1), 1-14 (2003).

[43] P. Flores and J. Ambrósio, "Revolute joints with clearance in multibody systems," *Comput. Struct.* **82**(17-19), 1359-1369 (2004).

[44] P. E. Nikravesh, *Computer-aided analysis of mechanical systems* (Prentice Hall, New Jersey, USA, 1988). pp. 1-388.

[45] G. N. Sandor and A. E. Erdman, *Advanced Mechanism Design: Analysis and Synthesis*, vol. 2 (Prentice Hall Inc, Englewood Cliffs, NJ, 1988) pp. 1-388.

[46] I. Khemili, M. A. Ben Abdallah and N. Aifaoui, "Multi-objective optimization of a flexible slider-crank mechanism synthesis, based on dynamic responses," *Eng. Optimiz.* **51**(6), 978-999 (2019).

[47] D. E. Golderberg, *Genetic Algorithms in Search, Optimization and Machine Learning* (Wesley, Boston, USA, 1953).

[48] S. Kirkpatrick, C. D. Gelatt and M. P. Vecchi, "Optimization by simulated annealing," *Science* **220**(4598), 671-680 (1983).

[49] R. Storn and K. Price, "Differential evolution – A simple and efficient heuristic for global optimization over continuous spaces," *J. Global Optim.* **11**(4), 341-359 (1997).

[50] M. Clerc, "Particle Swarm Optimization," *In: Introduction to Mathematical Optimization* vol. 21, (2006) pp. 107-117.

[51] E. Atashpaz-Gargari and C. Lucas, "Imperialist Competitive Algorithm: An Algorithm for Optimization Inspired by Imperialistic Competition," *In: 2007 IEEE Congress on Evolutionary Computation*, Singapore (IEEE, 2007) pp. 4661-4667.

[52] D. Karaboga and B. Basturk, On the performance of artificial bee colony (ABC) algorithm," *J. Theor. Appl. Inf. Technol.* **47**(2), 434-459 (2013).

[53] M. E. El-Hawary, *Modern Heuristic Optimization Techniques, Theory and Applications to Power Systems* (John Wiley, New Jersey, 2008) pp. 71-83.

[54] M. A. Ben Abdallah, I. Khemili and N. Aifaoui, "Flexible slider crank mechanism synthesis using meta-heuristic optimization techniques: A new designer tool assistance for a compliant mechanism synthesis," *Artif. Intell. Rev.* **53**(4), 2809-2840 (2020). doi: 10.1007/s10462-019-09747-y.

[55] J.-W. Kim, T. W. Seo and J. Kim, "A new design methodology for four-bar linkage mechanisms based on derivations of coupler curve," *Mech. Mach. Theory* **100**, 138-154 (2016).

[56] Q. Zhang, X. Ning, Y. Li, L. Pan, R. Gao and L. Zhang, "Path planning of patrol robot based on modified grey wolf optimizer," *Robotica* **41**(7), 1947-1975 (2023).

[57] Y.-J. Pak, Y.-S. Kong and J.-S. Ri, "Robust PID optimal tuning of a Delta parallel robot based on a hybrid optimization algorithm of particle swarm optimization and differential evolution," *Robotica* **41**(4), 1159-1178 (2023).

[58] K. Zhang, M. Yang, Y. Zhang and Q. Huang, "Error feedback method (EFM) based dimension synthesis optimisation for four-bar linkage mechanism," *Appl. Soft Comput.* **144**, 110424 (2023).

[59] J. Liang, H. Guo, K. Chen, K. Yu, C. Yue and X. Li, "An improved Kalman particle swarm optimization for modeling and optimizing of boiler combustion characteristics," *Robotica* **41**(4), 1087-1097 (2023).

[60] B. K. Rout and R. K. Mittal, "Optimal design of manipulator parameter using evolutionary optimization techniques," *Robotica* **28**(3), 381-395 (2010).

[61] J. L. Torres-Moreno, N. C. Cruz, J. D. Álvarez, J. L. Redondo and A. Giménez-Fernandez, "An open-source tool for path synthesis of four-bar mechanisms," *Mech. Mach. Theory* **169**, 104604 (2022).
Appendix A

<figure><img src="image_32.png" /><figcaption>**Figure A1.** Histograms of the objective function computed in 100 runs of the slider velocity synthesis.</figcaption></figure>
<table><caption>**Table AI. Classification of the bibliographic references.**</caption><thead><tr><td>Reference</td><td>Mechanism</td><td>Synthesis-based response</td><td>Flexible or rigid multibody systems</td><td>Synthesis technique</td></tr></thead><tbody><tr><td>Shaoping Bai et al., [1]</td><td>Four-bar mechanism</td><td>Path generation</td><td>Rigid</td><td>Algebraic formulation</td></tr><tr><td>Sang Min Han et al., [2]</td><td>Four-bar mechanism</td><td>Non-prescribed time for the generated path</td><td>Rigid</td><td>Fourier descriptors</td></tr><tr><td>Yixin Shao et al., [3]</td><td>Slider-crank mechanism and a seven-bar mechanism</td><td>Path generation</td><td>Rigid</td><td>The genetic algorithm</td></tr><tr><td>Saeed Ibrahimi et al., [4]</td><td>Four-bar mechanism</td><td>Path generation</td><td>Rigid</td><td>Genetic algorithm, differential evolution, and imperialist Competitive algorithm</td></tr><tr><td>Volkan Parlaktaş et al., [5]</td><td>Geared four-bar mechanism</td><td>Path generation</td><td>Rigid</td><td>Derived transmission angle</td></tr><tr><td>Sahand Hadizadeh Kafash et al., [6]</td><td>Four-bar mechanism</td><td>Path generation</td><td>Rigid</td><td>Differential Evolution (DE)</td></tr><tr><td>G. Ganesan et al., [7]</td><td>Adjustable four-bar mechanism</td><td>Path generation</td><td>Rigid</td><td>Genetic Algorithm and Pattern Search</td></tr><tr><td>Jianwei Sun et al., [8]</td><td>Planar four-bar mechanism</td><td>Motion generation</td><td>Rigid</td><td>Haar wavelet series theory</td></tr><tr><td>Jean-François Collard et al., [9]</td><td>Four-bar and a six-bar mechanism with extensible</td><td>Motion generation</td><td>Rigid</td><td>The nucleation method</td></tr><tr><td>Venkatar amanujam et al., [10].</td><td>Four-bar mechanism</td><td>Motion generation</td><td>Rigid</td><td>Algebraic method</td></tr><tr><td>Rasim Alizade et al., [11]</td><td>Spatial spherical four bar mechanism</td><td>Function generation</td><td>Rigid</td><td> Quaternion algebra and Chebychev approximation</td></tr></tbody></table>
<table><caption>Table AI. Continued.</caption><thead><tr><td>Reference</td><td>Mechanism</td><td>Synthesis-based response</td><td>Flexible or rigid multibody systems</td><td>Synthesis technique</td></tr></thead><tbody><tr><td>Gholaminezhad et al., [20]</td><td>Two fingers gripper</td><td>Griping forces and acting transmission forces.</td><td>Rigid</td><td>DE</td></tr><tr><td>Lin et al., [21]</td><td>Belt mechanism</td><td>Target path</td><td>Rigid</td><td>Transformation of pole maps</td></tr><tr><td>Yuanchang Liu et al., [22]</td><td>Unmanned vehicle</td><td>Path planning</td><td>Rigid</td><td>PSO/ GA</td></tr><tr><td>Tan Zhang [23]</td><td>Self-reconfigurable robot</td><td>Motion planning</td><td>Rigid</td><td>Genetic algorithm</td></tr><tr><td>Salvador Botello-Aceves et al., [24]</td><td>Parallel manipulator</td><td>End-effector position and orientation</td><td>Rigid</td><td>Evolutionary algorithms and quasi-Newton method</td></tr><tr><td>J. Brinker et al., [25].</td><td>Delta parallel robot</td><td>Motion transmission indices</td><td>Rigid</td><td>Nondominated Sorting Genetic Algorithm II</td></tr><tr><td>Gaurav Gupta et al., [26]</td><td>biped robot</td><td>Trajectory generation</td><td>Rigid</td><td>The genetic algorithm coupled to the artificial neural network</td></tr><tr><td>Yaldiz et al., [29], [30], [31], [32] Panagant et al., [32]</td><td>Vehicle's seat bracket</td><td>Stress and mass</td><td>Rigid</td><td>New hybrid Taguchi salp swarm algorithm (HTSSA) Ecogeography-based algorithm Equilibrium algorithm Marine predators algorithm Slime mold algorithm</td></tr><tr><td>Yaldiz et al., [33]</td><td>Robot gripper</td><td>Grip force</td><td>Rigid</td><td>HGOANM</td></tr><tr><td>Yaldiz et al., [56]</td><td>Design of gear train</td><td>Design parameters</td><td>Rigid</td><td>Chaotic Lévy flight distribution (CLFD) algorithm</td></tr><tr><td>Bureerat et al., [34] Yaldiz [35]</td><td>Crashworthiness of a car</td><td>Topology optimization</td><td>Rigid</td><td>Multi-surrogate-assisted optimization</td></tr></tbody></table>
<table><caption>Table AI. Continued.</caption><thead><tr><th>Reference</th><th>Mechanism</th><th>Synthesis-based response</th><th>Flexible or rigid multibody systems</th><th>Synthesis technique</th></tr></thead><tbody><tr><td>Zhen Luo et al., [12]</td><td>Compliant micro-inverter</td><td>The geometrical advantage (GA), mechanical efficiency (ME), and mechanical advantage (MA)</td><td>Flexible</td><td>Topology optimization</td></tr><tr><td>Matteo Russo et al., [13]</td><td>Six degrees of freedom Robot arm</td><td>The dimensional and pose error</td><td>Flexible</td><td>MATLAB optimization toolbox fmincon</td></tr><tr><td>Sangamesh R. Deepak et al., [14]</td><td>Compliant inverter, gripper, and crimper.</td><td>Mechanism advantage, efficiency, characteristic stiffness, and artificial spring formulations</td><td>Flexible</td><td>MATLAB optimization toolbox</td></tr><tr><td>Abdelkhalick Mohammad et al., [15]</td><td>A 16-degree-of-freedom snake robots</td><td>Path optimization</td><td>Flexible</td><td>Online optimization algorithm</td></tr><tr><td>Anzhu Gao et al., [16]</td><td>Continuum robot</td><td>Path optimization</td><td>Flexible</td><td>fmincon function in MATLAB</td></tr><tr><td>Shinjishinjishinji Nishiwaki et al., [19]</td><td>Compliant gripper and clamp</td><td>Kinematic and structural aspects</td><td>Flexible</td><td>DSPLP from the SLATEC library of simplex method</td></tr><tr><td>Mohamed Taha Chikhaoui et al., [17]</td><td>Continuum robot</td><td>End-effector distance control and trajectory tracking</td><td>Flexible</td><td>Particle swarm optimization (PSO) algorithm</td></tr><tr><td>Liu Wang et al., [18]. This work</td><td>Soft medical robot<br>Slider-crank<br>Four-bar</td><td>Path optimization<br>- Motion generation<br>- Path generation<br>- Axial deformation of the flexible body</td><td>Flexible<br>Flexible</td><td>The genetic algorithm<br>GA, PSO, SA, DE, AC, ABC, and ICA</td></tr></tbody></table>
<figure><img src="image_33.png" /><figcaption>**Figure A2.** The automatic interaction between optimization module and dynamic analysis module for the Computer-Aided Design for Multibody Systems.</figcaption></figure>
<table><caption>**Table AII.** The obtained mechanism design variables based on the slider velocity synthesis.</caption><thead><tr><th></th><th>l<sub>2</sub> (mm)</th><th>l<sub>3</sub> (mm)</th><th>ρ (kg/m<sup>3</sup>)</th><th>E (kPa)</th><th>h (mm)</th><th>a (mm)</th><th>Error</th><th>CPU time (sec)</th><th>μ</th><th>σ</th></tr></thead><tbody><tr><td>GA</td><td>50.44</td><td>363.41</td><td>2980.9</td><td>7.491 10<sup>10</sup></td><td>1.6945</td><td>27.931</td><td>3.0858 10<sup>-3</sup></td><td>1061.1</td><td>3.0938 10<sup>-3</sup></td><td>7.55559 10<sup>-5</sup></td></tr><tr><td>SA</td><td>50.0001</td><td>350.0000</td><td>2.1493 10<sup>3</sup></td><td>6.7575 10<sup>10</sup></td><td>1</td><td>28.3475</td><td>7.2296 10<sup>-7</sup></td><td>26009.2</td><td>7.21652 10<sup>-7</sup></td><td>1.80909 10<sup>-8</sup></td></tr><tr><td>DE</td><td>50.0000</td><td>350.0000</td><td>2.5909 10<sup>3</sup></td><td>6.5000 10<sup>10</sup></td><td>1.2273</td><td>25.6861</td><td>1.8503 10<sup>-17</sup></td><td>9050.9</td><td>2.26308 10<sup>-17</sup></td><td>5.2379 10<sup>-17</sup></td></tr><tr><td>PSO</td><td>50.0000</td><td>349.9999</td><td>2.468 10<sup>3</sup></td><td>6.5844 10<sup>10</sup></td><td>2</td><td>30</td><td>1.0156 10<sup>-9</sup></td><td>6390.9</td><td>9.95411 10<sup>-10</sup></td><td>4.95851 10<sup>-10</sup></td></tr><tr><td>ICA</td><td>50.0000</td><td>350.0001</td><td>3.000 10<sup>3</sup></td><td>7.5 10<sup>10</sup></td><td>1</td><td>30</td><td>1.2491 10<sup>-8</sup></td><td>9325.2</td><td>2.64607 10<sup>-8</sup></td><td>8.00423 10<sup>-8</sup></td></tr><tr><td>ABC</td><td>50.0006</td><td>349.9535</td><td>2.3041 10<sup>3</sup></td><td>7.0635 10<sup>10</sup></td><td>1.1502</td><td>29.3891</td><td>6.3428 10<sup>-6</sup></td><td>12567.9</td><td>6.51672 10<sup>-6</sup></td><td>5.29079 10<sup>-7</sup></td></tr><tr><td>AC</td><td>50.0000</td><td>350.0000</td><td>2.6282 10<sup>3</sup></td><td>7.1090 10<sup>10</sup></td><td>1.5753</td><td>29.4323</td><td>0</td><td>3248.6</td><td>4.10942 10<sup>-19</sup></td><td>3.8388 10<sup>-17</sup></td></tr></tbody></table>
<figure><img src="image_34.png" /><figcaption>**Figure A3.** Histograms of the objective function computed in 100 runs of the axial displacement synthesis.</figcaption></figure>
<table><caption>**Table AIII.** The obtained mechanism design variables based on the flexible connecting rod midpoint axial displacement.</caption><thead><tr><th></th><th>l<sub>1</sub>(mm)</th><th>l<sub>2</sub>(mm)</th><th>Rho(kg/m<sup>3</sup>)</th><th>Young's modulus E (kPa)</th><th>h (mm)</th><th>a (mm)</th><th>Error</th><th>CPU time (sec)</th><th>μ</th><th>σ</th></tr></thead><tbody><tr><td>GA</td><td>39.3</td><td>396.1</td><td>2.2711 10<sup>3</sup></td><td>7.262 10<sup>10</sup></td><td>1.3835</td><td>20.052</td><td>2.6953 10<sup>-2</sup></td><td>1060.3</td><td>2.6953 10<sup>-2</sup></td><td>3.85768 10<sup>-10</sup></td></tr><tr><td>SA</td><td>61.4654</td><td>388.3530</td><td>2.8002 10<sup>3</sup></td><td>6.6418 10<sup>10</sup></td><td>1.9983</td><td>27.3342</td><td>2.3108 10<sup>-3</sup></td><td>24428.294</td><td>2.3108 10<sup>-3</sup></td><td>2.97019 10<sup>-10</sup></td></tr><tr><td>DE</td><td>56.0655</td><td>369.4410</td><td>2.2691 10<sup>3</sup></td><td>7.3051 10<sup>10</sup></td><td>1.5556</td><td>30</td><td>1.7181 10<sup>-3</sup></td><td>6108.869</td><td>1.7181 10<sup>-3</sup></td><td>7.74508 10<sup>-10</sup></td></tr><tr><td>PSO</td><td>70</td><td>414.3858</td><td>2.3525 10<sup>3</sup></td><td>6.9302 10<sup>10</sup></td><td>2</td><td>30</td><td>3.8959 10<sup>-3</sup></td><td>6327</td><td>3.8959 10<sup>-3</sup></td><td>4.0529 10<sup>-10</sup></td></tr><tr><td>ICA</td><td>59.4769</td><td>381.9328</td><td>2.7804 10<sup>3</sup></td><td>7.5 10<sup>10</sup></td><td>1.7761</td><td>30.0000</td><td>1.9122 10<sup>-3</sup></td><td>8037.070s</td><td>1.9122 10<sup>-3</sup></td><td>5.55546 10<sup>-10</sup></td></tr><tr><td>ABC</td><td>55.9045</td><td>370.0060</td><td>2.5070 10<sup>3</sup></td><td>6.9703 10<sup>10</sup></td><td>1.6420</td><td>29.4908</td><td>1.2559 10<sup>-3</sup></td><td>12723.194</td><td>1.2559 10<sup>-3</sup></td><td>1.04487 10<sup>-9</sup></td></tr><tr><td>AC</td><td>52.0020</td><td>357.4133</td><td>2.623810<sup>3</sup></td><td>6.873610<sup>10</sup></td><td>1.6086</td><td>22.0406</td><td>7.2535 10<sup>-4</sup></td><td>2206.752</td><td>7.2535 10<sup>-4</sup></td><td>6.72066 10<sup>-10</sup></td></tr></tbody></table>
<figure><img src="image_35.png" /><figcaption>**Figure A4.** Histograms of the objective function computed in 100 runs of the combined synthesis.</figcaption></figure>
<table><caption>Table AIV. The obtained mechanism design variables based on the combined synthesis.</caption><thead><tr><th></th><th>l<sub>1</sub>(mm)</th><th>l<sub>2</sub>(mm)</th><th>ρ (kg/m<sup>3</sup>)</th><th>Young's modulus (kPa)</th><th>h (mm)</th><th>a (mm)</th><th>Error</th><th>CPU time (sec)</th><th>μ</th><th>σ</th></tr></thead><tbody><tr><td>GA</td><td>52.074</td><td>369.06</td><td>2.6627 10<sup>3</sup></td><td>7.167 10<sup>10</sup></td><td>1.6448</td><td>20.535</td><td>1.5546 10<sup>-2</sup></td><td>1999.1</td><td>1.5546 10<sup>-2</sup></td><td>5.89379 10<sup>-10</sup></td></tr><tr><td>SA</td><td>49.9977</td><td>349.1044</td><td>2.3516 10<sup>3</sup></td><td>7.3308 10<sup>10</sup></td><td>1.4114</td><td>21.1616</td><td>4.0028 10<sup>-4</sup></td><td>26364.162</td><td>4.0028 10<sup>-4</sup></td><td>3.7539 10<sup>-10</sup></td></tr><tr><td>DE</td><td>50.0028</td><td>350.3338</td><td>2.000 10<sup>3</sup></td><td>6.5663 10<sup>10</sup></td><td>1.3812</td><td>22.8525</td><td>1.9583 10<sup>-4</sup></td><td>6405.587</td><td>1.9583 10<sup>-4</sup></td><td>5.42008 10<sup>-10</sup></td></tr><tr><td>PSO</td><td>49.9865</td><td>349.9833</td><td>2.7070 10<sup>3</sup></td><td>7.3394 10<sup>10</sup></td><td>1.4881</td><td>30</td><td>7.9712 10<sup>-5</sup></td><td>4208.386</td><td>7.9712 10<sup>-5</sup></td><td>6.3814 10<sup>-10</sup></td></tr><tr><td>ICA</td><td>49.9955</td><td>349.5904</td><td>2.4743 10<sup>3</sup></td><td>6.5 10<sup>10</sup></td><td>1.5093</td><td>30</td><td>2.4199 10<sup>-4</sup></td><td>7820.585 s</td><td>2.4199 10<sup>-4</sup></td><td>2.74033 10<sup>-10</sup></td></tr><tr><td>ABC</td><td>50.3554</td><td>350.5856</td><td>2.5092 10<sup>3</sup></td><td>7.2491 10<sup>10</sup></td><td>1.4500</td><td>25.8143</td><td>2.1640 10<sup>-3</sup></td><td>13009.927</td><td>2.164 10<sup>-3</sup></td><td>2.74033 10<sup>-10</sup></td></tr><tr><td>AC</td><td>49.9993</td><td>350.3910</td><td>2.4902 10<sup>3</sup></td><td>6.7035 10<sup>10</sup></td><td>1.5263</td><td>29.4434</td><td>1.6045 10<sup>-4</sup></td><td>3359.949</td><td>1.6045 10<sup>-4</sup></td><td>1.5847 10<sup>-10</sup></td></tr></tbody></table>
<table><caption>**Table AV.** The obtained results for the flexible four-bar mechanism based on the combined synthesis.</caption><thead><tr><th></th><th>l<sub>1</sub>(mm)</th><th>l<sub>2</sub>(mm)</th><th>l<sub>3</sub>(mm)</th><th>l<sub>4</sub>(mm)</th><th>Error</th><th>μ</th><th>σ</th></tr></thead><tbody><tr><td>GA</td><td>274.7782</td><td>100.1322</td><td>248.2180</td><td>248.1663</td><td>2.30 10<sup>-1</sup></td><td>2.3 10<sup>-1</sup></td><td>4.18861 10<sup>-8</sup></td></tr><tr><td>SA</td><td>289.4258</td><td>100.2761</td><td>253.2160</td><td>253.2160</td><td>8.09 10<sup>-2</sup></td><td>8.09 10<sup>-2</sup></td><td>4.88032 10<sup>-9</sup></td></tr><tr><td>DE</td><td>287.7952</td><td>100.2856</td><td>252.1081</td><td>252.1011</td><td>8.31 10<sup>-2</sup></td><td>8.31 10<sup>-2</sup></td><td>7.52815 10<sup>-9</sup></td></tr><tr><td>PSO</td><td>280.3504</td><td>100.3652</td><td>251.1854</td><td>251.1537</td><td>1.35 10<sup>-1</sup></td><td>1.35 10<sup>-1</sup></td><td>3.55105 10<sup>-9</sup></td></tr><tr><td>ICA</td><td>284.7969</td><td>100.2849</td><td>251.3476</td><td>251.3331</td><td>1.0563 10<sup>-1</sup></td><td>1.0563 10<sup>-1</sup></td><td>4.90483 10<sup>-9</sup></td></tr><tr><td>ABC</td><td>285.0343</td><td>100.0724</td><td>249.5127</td><td>249.4702</td><td>1.2795 10<sup>-1</sup></td><td>1.2795 10<sup>-1</sup></td><td>1.76551 10<sup>-9</sup></td></tr><tr><td>AC</td><td>287.4285</td><td>100.4849</td><td>252.1730</td><td>252.1738</td><td>8.55424 10<sup>-2</sup></td><td>8.55424 10<sup>-2</sup></td><td>6.34207 10<sup>-9</sup></td></tr></tbody></table>

<figure><img src="image_36.png" /><figcaption>**Figure A5.** Histograms of the objective function computed in 100 runs for the flexible four-bar mechanism based on the combined synthesis.</figcaption></figure>
Appendix B. Glossary.

<table><thead><tr><th>Variable</th><th>Designation</th></tr></thead><tbody><tr><td>u<sub>1</sub>, u<sub>2</sub>, u<sub>3</sub>, u<sub>4</sub>, u<sub>5</sub>, u<sub>6</sub></td><td>The general coordinates of the flexible beam</td></tr><tr><td>ü<sub>a</sub></td><td>The beam's absolute acceleration</td></tr><tr><td>ü<sub>r</sub></td><td>The rigid element acceleration,</td></tr><tr><td>ü</td><td>The generalized relative acceleration,</td></tr><tr><td>a<sub>n</sub></td><td>Normal acceleration</td></tr><tr><td>a<sub>c</sub></td><td>Coriolis acceleration</td></tr><tr><td>a<sub>t</sub></td><td>Tangential accelerations</td></tr><tr><td>v(¯x, t)</td><td>The axial displacement taking place along the x-axis</td></tr><tr><td>x&#773;</td><td>The x-axis</td></tr><tr><td>t</td><td>Time</td></tr><tr><td>w(¯x, t)</td><td>The axial displacement taking place along the y-axis</td></tr><tr><td>φ<sub>1</sub>(¯x), φ<sub>2</sub>(¯x), φ<sub>3</sub>(¯x), φ<sub>4</sub>(¯x), φ<sub>5</sub>(¯x), φ<sub>6</sub>(¯x)</td><td>The shape functions</td></tr><tr><td>T</td><td>The kinetic energy</td></tr><tr><td>U</td><td>The deformation energy</td></tr><tr><td>Q&#773;<sub>i</sub></td><td>The external forces applied to the beam element</td></tr><tr><td>ü&#773;, ü<sup>T</sup></td><td>The first derivative with respect of time for the local coordinates and its transpose</td></tr><tr><td>m&#773;</td><td>The elementary inertia matrix in the local coordinates of the beam</td></tr><tr><td>k&#773;</td><td>The elementary stiffness matrix in the local coordinates of the beam</td></tr><tr><td>U</td><td>The global coordinates</td></tr><tr><td>U&#773;</td><td>The first derivative with respect to time for the global coordinates</td></tr><tr><td>m</td><td>The global mass matrix</td></tr><tr><td>k</td><td>The global stiffness matrix</td></tr><tr><td>R</td><td>The rotation matrix</td></tr><tr><td>M</td><td>The assembling mass matrix for two elements</td></tr><tr><td>K</td><td>The assembling stiffness matrix for two elements</td></tr><tr><td>C</td><td>The assembling damping matrix for two elements</td></tr></tbody></table>

Appendix C. List of abbreviations

<table><thead><tr><th>Acronym</th><th>Significance</th></tr></thead><tbody><tr><td>PSO</td><td>Particle swarm optimization</td></tr><tr><td>GA</td><td>Genetic algorithm</td></tr><tr><td>ICA</td><td>Imperialist competitive algorithm</td></tr><tr><td>ABC</td><td>Artificial bee colony</td></tr><tr><td>DE</td><td>Differential evolution</td></tr><tr><td>SA</td><td>Simulating annealing</td></tr><tr><td>AC</td><td>Ant colony</td></tr><tr><td>CPU</td><td>Computational processing unit</td></tr></tbody></table>

Cite this article: M.A. Ben Abdallah, I. Khemili, N. Aifaoui and M.A. Laribi (2024). "Computer-aided design tool for typical flexible mechanisms synthesis by means of evolutionary algorithms", Robotica **42**, 1172–1211. https://doi.org/10.1017/S0263574724000171