A Jacobian-based topology optimization method is recently proposed for compliant parallel mechanisms (CPMs), in which the CPMs' Jacobian matrix and characteristic stiffness are optimized simultaneously to achieve kinematic and stiffness requirement, respectively. Lately, it is found that the characteristic stiffness fails to ensure a valid topology result in some particular cases. To solve this problem, an improved stiffness evaluation based on the definition of stiffness is adopted in this paper. This new stiffness evaluation is verified and compared with the characteristic stiffness by using several design examples. In addition, several typical benchmark problems (e.g., displacement inverter, amplifier, and redirector) are solved by using the Jacobian-based topology optimization method to show its general applicability.

## Introduction

Compliant mechanisms are elastic structures that can transmit force or motion from input to output. Due to the combined characteristics of mechanism and structure, the analysis and design of compliant mechanisms are more challenging than that of rigid-body mechanisms. There are two main synthesis approaches for compliant mechanisms, i.e., rigid-body replacement approach [1–3] and topology optimization approach [4–7].

The rigid-body replacement approach has a wide design scope including all kinds of degree-of-freedom (DOF) compliant mechanisms. This approach synthesizes most of compliant parallel mechanisms (CPMs), whose synthesis method is the focus of this paper, by replacing the kinematic joints of existing rigid-body mechanisms with flexure hinges. Thus, the compliant mechanisms designed by this approach rely on the topologies of rigid-body mechanisms. While this approach is successful in designing multi-DOF CPMs for precision applications, it is limited by the fact that a compliant mechanism may still be unable to fully reproduce the motion of its rigid-body counterpart even using rigorous analysis and optimization techniques [8]. Moreover, this approach cannot select the best topology for a specific problem currently, which is quite important for developing mechanisms with high performance.

The topology optimization approach regards the synthesis of compliant mechanism as finding out the optimal material distribution within a given design domain, by maximizing the motion or force transmission between specific input and output ports. Due to this advantage, the topology optimization approach has been successfully applied to the synthesis problems of multiple input and output ports [5,9–13], multiple materials or physics [9–11], three-dimensional simple compliant mechanisms [14–16], etc. The compliant mechanisms designed by this approach possess structural type of topology, i.e., no flexure hinges, which is different from the flexure hinge-based mechanisms obtained by the rigid-body replacement approach.

To introduce the idea of topology optimization into the design scope of the rigid-body replacement approach, our previous works [17–19] tried to synthesize the flexure-based compliant mechanisms with simple motion based on the idea of topology optimization. Lum et al. [20–22] presented a hybrid topological and structural optimization method. This method first synthesizes the compliant joints with the optimal stiffness characteristics by topology optimization. The resulting compliant joints are then assembled into a CPM based on existing rigid-body mechanism topology.

Recently, we proposed a Jacobian-based topology optimization method [23] for the optimal synthesis of planar CPMs. Traditional topology optimization methods realize multiple outputs by predefining specific output displacements at output ports. The premise of this realization is that the position and direction of output ports are known. However, the output motion of multi-DOF CPMs is unknown. The traditional way of predefining specific output ports cannot be applied to multi-DOF CPMs topology optimization directly. To solve this problem, Jacobian matrix [24,25] is introduced into the field of topology optimization by the proposed method. The Jacobian matrix describes all the freedoms of CPM's mobile platform in a unified and concise form, and contains the information of CPM's DOF and direct kinematics simultaneously. By optimizing the Jacobian matrix, one can synthesize a compliant mechanism with desired DOF (mechanism's function) and optimized direct kinematics (mechanism's performance). In addition to the above kinematic realization, the mechanism's input and output characteristic stiffness [26,27] (C-stiffness for short) are also optimized to achieve enough stiffness to bear the external loads.

Lately, we found that the C-stiffness fails to ensure a valid topology result in some particular problems. In this paper, an improved stiffness evaluation based on its definition is incorporated into the problem formulation, and will be compared with the C-stiffness formulation by using several design problems.

The rest of this paper is organized as follows: Section 2 describes the problem formulations of Jacobian-based topology optimization method. Section 3 illustrates the topology analysis of CPMs. The sensitivity analysis and optimization algorithm are described in Sec. 4. Section 5 gives numerical examples to verify the stiffness formulation. The conclusions are presented in Sec. 6.

## Jacobian-Based Topology Optimization Method

### Problem Description.

The topology of a CPM is determined by the number, arrangement, and topology structures of its constituent compliant limbs. The proposed method regards the problem of CPMs topology optimization as finding the best topology of compliant limbs within several given design domains.

A CPM composed of *n* compliant limbs is shown in Fig. 1. The mechanism's compliant limbs are assumed to be synthesized within *n* predefined design domains, $\Gamma ={\Omega 1,\Omega 2,\u2026,\Omega n}$. The corresponding design variables to these design domains are illustrated by $X={X1,X2,\u2026,Xn}$. The number of compliant limbs, *n*, is determined by the number of CPM's DOF. The position of input points *a _{i}* ($i=1,2,\u2026,n$) and output point

*o*is defined by designers.

### Properties of Jacobian Matrix.

**J**describes transmission relation between the displacements of input and output freedoms, i.e., forward kinematics of the CPM

where $Uo=[ux,uy,\theta z]T$ contains the three freedoms' displacement of output point *o* at the mobile platform, and $Ua=[u1,u2,\u2026,un]T$ is the vector of input displacements at input points *a _{i}*. Since the compliant mechanisms are usually driven by linear motion actuators, we limit the displacements in

**U**

*to the translational freedoms of input points.*

_{a}The element *J _{ji}* in

**J**represents the geometry advantage (GA) between the

*j*th freedom of

**U**

*and the*

_{o}*i*th input. If all the elements of the

*j*th row in

**J**are equal or close to zero, all the inputs in

**U**

*will not produce displacement in the*

_{a}*j*th freedom of mobile platform, i.e., this freedom is suppressed. While the other freedoms corresponding to nonzero row vectors are considered as the CPM's DOF. By maximizing the absolute value of the elements in the nonzero row vectors, the forward kinematics and motion transmission performance of CPM can be optimized. Thus, the Jacobian matrix contains the information of CPM's function and performance simultaneously. By optimizing the Jacobian matrix, we can synthesize a compliant mechanism with desired DOF to realize the mechanism's function and optimized direct kinematics to achieve higher performance.

### Problem Formulation.

#### Differences Between C-Stiffness and Stiffness.

*f*

_{in}and

*f*

_{out}) and displacements (

*u*

_{in}and

*u*

_{out}) at input and output ports can be described by a mechanism stiffness matrix

**K**

*[28] shown in the following equation:*

_{m}**K**

*[28] are the C-stiffness [26]. As shown in Eq. (3), they incorporate the C-stiffness at input and output ports (*

_{m}*k*

_{11}and

*k*

_{22}) into the formulation to achieve topology optimization of hinge-free compliant mechanisms

where GA is the geometry advantage of the mechanism, and GA* is the desired geometry advantage, *f* and *S* represent the kinematic and stiffness requirement, respectively.

According to Eq. (2), the physical meanings of *k*_{11} and *k*_{21} are the forces that should be acted on input and output ports, if *u*_{in} = 1 and *u*_{out} = 0 are expected. The input C-stiffness *k*_{11} just describes part of the stiffness relationship between *f*_{in} and *u*_{in}, and so does the output C-stiffness *k*_{22}.

*f*

_{in}= 1 and

*f*

_{out}= 0, and solve Eq. (2). The resulting input displacement

*u*

_{in}is the input compliance

*f*

_{in}= 0 and

*f*

_{out}= 1, and solve Eq. (2). The resulting output displacement

*u*

_{out}is the output compliance

By comparing *k*_{11} and *k*_{22} with *k*_{in} and *k*_{out}, respectively, C-stiffness is part of the stiffness. While the C-stiffness has successfully evaluated the mechanism's stiffness property in many design problems, some problems show that the stiffness is more reliable than C-stiffness. Thus, this paper will modify the formulations of Jacobian-based topology optimization method by replacing the C-stiffness with stiffness.

#### Formulations Using New Stiffness Evaluation.

*f*and stiffness requirement

*S*are two conflicting subobjectives. On one hand, the mechanism should be soft enough to deform and deliver motion. On the other hand, it should be stiff enough to transmit forces to the mobile platform and bear the external force. A general problem formulation for CPMs topology optimization can be written as follows:

where *ω* (0 < *ω* < 1) is the weight indicating the relative significance of kinematic requirement *f*. *V*(**X**) is the volume fraction of topology candidate, and *V _{o}* is the allowed volume fraction.

*f*has two different forms according to the design problem. The first form of kinematic requirement

*f*

_{1}is suitable for the CPMs whose kinematics is simple enough to be predefined by designers. This form is to force the Jacobian matrix

**J**of CPM to be close to a desired Jacobian matrix

**J**

^{*}by minimizing the differences between

*J*and $Jji*$, which are the elements of

_{ji}**J**and

**J**

^{*}, respectively. As a result, the desired DOF and kinematic properties of CPM can be expressed in the desired Jacobian matrix

**J**

^{*}

*J*evaluates the workspace in the

_{j}*j*th freedom, which is the quadratic sum of the elements in corresponding row vector shown in Eq. (10). $Jjd$ is the workspace of a desired freedom, while $Jjc$ is the workspace of a constrained freedom. The natural exponential function forces each $Jjc$ to be close to zero

*S*is achieved by maximizing the input and output C-stiffness of the CPM. This paper uses the stiffness calculated by its definition as a new stiffness evaluation instead of the C-stiffness. Mathematically, the stiffness requirement is formulated as the product of the input stiffness $kini$ and output stiffness $kouti$

### Unification of the Units in Rotational and Translational Freedoms.

*l*is introduced in this paper. Based on the characteristic length

_{c}*l*, we can define equivalent moment $M\u0303$ and equivalent rotational displacement $\theta \u0303z$ as follows:

_{c}*J*is the ratio of the

_{ji}*j*th freedom of

**U**

*and the*

_{o}*i*th input. Since all the inputs are translational freedoms, only the rotational displacement

*θ*in

_{z}**U**

*should be transformed into equivalent rotation displacement $\theta \u0303z$. Thus, the elements in the third row vector of*

_{o}**J**are multiplied by

*l*to obtain equivalent Jacobian matrix $J\u0303$

_{c}## Topology Analysis

This section shows how the Jacobian matrix, and input and output stiffnesses of CPMs can be obtained by using the finite element analysis and matrix methods [31].

### Discretization and Parameterization.

*and part of the mobile platform $\Omega pi$. The compliant limb ($\Omega i+\Omega pi$) is discretized by using the quadrilateral elements that possess three freedoms (*

_{i}*u*,

_{x}*u*,

_{y}*θ*) at each node, and parameterized by the simplified isotropic material with penalization scheme [32]. The stiffness matrix of the

_{z}*i*th compliant limb can be obtained by the following equation:

where *N _{i}* is the number of the elements in design domain Ω

*, $Npi$ is the number of the elements in $\Omega pi$,*

_{i}**K**

*is the element stiffness matrix in the global level,*

_{e}*x*is the material density (design variable) of each element in Ω

^{e}*with value between the lower limit $xmine$ (void) and 1 (solid),*

_{i}*ρ*is the penalty factor, the elements in $\Omega pi$ are solid.

### Stiffness Modeling of CPMs.

*a*and the three freedoms of endpoint

_{i}*o*are considered in its stiffness modeling. First, a compliance matrix $Caoi$ that characterizes the compliance relationship between these four concerned freedoms is calculated according to physical meaning of compliance. Four load cases $Fj\u2009(j=1,2,\u2026,4)$, in which a unit dummy load is applied to each concerned freedom of

_{i}*a*and

_{i}*o*in sequence, are used to calculate corresponding displacements by solving the following equation:

_{i}*C*in a compliance matrix

_{jk}**C**is the displacement of the

*j*th freedom due to a unit load that only acts on the

*k*th freedom. According to this physical meaning, the element of $Caoi$ is obtained by Eq. (19). The displacements of the four concerned freedoms in displacement vector

**U**

*form the*

_{k}*k*th column in $Caoi$

*o*is transferred into the coordinate of the output point

_{i}*o*by a transformation matrix $Tio$, whereas the compliance of the input freedom at

*a*remains in its local coordinate. The inverse of the resulting compliance matrix, i.e., stiffness matrix $Kaoi$, is shown in Eq. (20). For more information about the transformation matrix $Tio$, please refer to Refs. [31,33].

_{i}*o*in $Kaoi$ of all the compliant limbs are superimposed to form stiffness of the output point

*o*at mobile platform, while the stiffness of all the input freedoms at points

*a*($i=1,2,\u2026,n$) remain in their local coordinates. The transformation is shown in Eq. (22). The resulting

_{i}**K**

*is the mechanism stiffness matrix that characterizes the stiffness relationship between the*

_{m}*n*input freedoms and the three freedoms of output point

*o*

### Kinematic Analysis of CPMs.

**U**

*and*

_{a}**U**

*, respectively, for displacements of the input and output freedoms. As shown in Eq. (23), the mechanism loads are also partitioned into two sets as*

_{o}**F**

*and*

_{a}**F**

*, accordingly. This would in effect partition the mechanism stiffness matrix*

_{o}**K**

*into the following form:*

_{m}**F**

*=*

_{o}**0,**by solving the second equation in Eq. (23), the relationship between the input displacement

**U**

*and output displacement*

_{a}**U**

*, i.e., the Jacobian matrix*

_{o}**J**, can be obtained

### Input and Output Stiffness of CPMs.

**K**

*, we can calculate the input and output stiffness of CPMs according to the physical meaning of compliance. As shown in Eq. (25),*

_{m}*n*+ 3 unit load vectors ($Fini$ and $Foutj$) are applied to the freedoms of

**K**

*in sequence to obtain corresponding displacement vectors $Uini$ and $Uoutj$*

_{m}## Sensitivity Analysis

*J*, $kini$, and $koutj$. The sensitivity of

_{ji}*J*can be extracted from sensitivity of

_{ji}**J,**which is calculated as follows:

**K**

_{22}and

**K**

_{21}are parts of

**K**

*, the sensitivity of*

_{m}**J**is determined by the sensitivity of

**K**

*. According to Eq. (26), the sensitivity of $kini$ is calculated by Eq. (28), while the sensitivity of $koutj$ can be obtained by the same way*

_{m}On the basis of sensitivity analysis, the topology optimization problem is solved by modifying the 99 line matlab code proposed by Sigmund [34]. The optimality criteria-based optimizer and filtering technique [7] of the 99 line Matlab code are used to update the design variables of each domain and ensure existence of solutions, respectively. For each numerical example in Sec. 5, the initial design **X** is defined by setting the material density of each element to be the value of the allowed volume fraction, i.e., *x ^{e}* =

*V*. The convergence criterion is the change in design variables, which is set to 0.005 in this paper. The move limit in the heuristic updating scheme is 0.1. The filter radius

_{o}*r*

_{min}is set to 1.2, i.e., the filter length scale

*d*

_{min}(

*d*

_{min}= 2

*r*

_{min}) is 2.4. It should be pointed out that the volume constraint is active during the whole optimization process. For more detail about the optimality criteria-based optimizer and filtering technique, readers can refer to Ref. [34].

## Numerical Studies

This section will compare several topology optimization results obtained by using the C-stiffness and stiffness formulations, respectively. The artificial material properties for these examples are described as: Young's modulus is *E* = 1 GPa and Poisson's ratio is *υ* = 0.3. The characteristic length *l _{c}* is set to 10 mm in this study. All the numerical examples are carried out on a computer with Intel Core

*i*7 – 6700 (3.40 GHz) CPU, 8.00GB RAM, and Matlab

*R*2009

*a*. Note that the filter length scale

*d*

_{min}= 2.4 is represented by a red bar in each figure showing the final topology.

### Design of 2DOF CPMs.

In this section, a 2DOF CPM will be synthesized using two asymmetrical compliant limbs within the design domain shown in Fig. 4. As can be seen, the two compliant limbs have the same size and boundary conditions. Each compliant limb is discretized by 100 × 100 finite elements for elastic analysis. The allowable amount of material is 20%.

#### Solved by Using f_{1}.

**J**

^{*}, this design problem can be solved using the first form of objective function in Eq. (8). The

**J**

^{*}for this example is given by Eq. (33), in which

*a*

_{1}

*and*

_{x}*a*

_{2}

*are the input freedoms. The zero vector in the third row of the*

_{y}**J**

^{*}means that the freedom

*θ*should be suppressed, i.e., the CPM is expected to have only two translational freedoms (

_{z}*u*and

_{x}*u*). Moreover, the two translational freedoms are expected to be decoupled. For example, the input

_{y}*a*

_{1}

*only induces the translational displacement*

_{x}*u*with desired geometry advantage GA

_{x}^{*}and has no impact on the other two freedoms

The GA^{*} is first set to −3 and *ω* is set to 0.5. The corresponding topology optimization problem is solved by using the C-stiffness and stiffness formulations, respectively. The optimizations were run for 200 iterations. The resulting topologies of the two formulations are shown in Fig. 5, which shows that both of C-stiffness and stiffness formulations can obtain valid topologies in this case. Figures 6 and 7 show the iteration history of objective value, kinematic requirement *f,* and stiffness requirement *S* in the optimization process of the two formulations, respectively. It can be seen that oscillations exist in the iteration curves. The oscillations may be caused by the material distribution at some specific elements. Fortunately, the topologies at later period of iteration are stable, which can be regarded as the optimal topology. The corresponding **J** of the two final topologies are listed in the first two rows of Table 1. It should be noted that only the elements in the first column vector of **J** are displayed for brevity, since *J*_{11} ≃ *J*_{22}, *J*_{21} ≃ *J*_{12}, and *J*_{31} ≃ *J*_{32}. One can see that both of C-stiffness and stiffness formulations can force **J** of CPMs to be close to **J**^{*}, i.e., the kinematics requirement *f* is realized. In addition, the C-stiffness (*Ck _{a}*

_{1}and

*Ck*) and stiffness (

_{ox}*k*

_{a}_{1}and

*k*) of the input and output ports in

_{ox}*x*-axis of the two final topologies are given in the first two rows of Table 2. The results show that both of C-stiffness and stiffness formulations realize the stiffness requirement

*S*effectively. The resulting output stiffness is smaller than the input stiffness to achieve GA, whereas the value of C-stiffness is larger than stiffness for the same mechanism.

However, when GA^{*} is set to be positive, e.g., GA^{*} = 2.5, it is found that the C-stiffness formulation fails to ensure the stiffness requirement *S* and results in invalid topologies. Figure 8(a) shows the final topology obtained by using the C-stiffness in the case of GA^{*} = 2.5. Although the input and output points of the CPM are connected by solid material successfully, there is no material connection between the compliant limbs and fixed ports. Consequently, the displacements of the input and output ports in one direction are equal, i.e., *J*_{11} = 1 (in the third row of Table 1). As shown in the third row of Table 2, the C-stiffness of the input and output ports in *x*-axis is 101.4 N/mm, whereas the corresponding stiffness *k _{a}*

_{1}and

*k*of this invalid topology are approximate to zero. Obviously, the C-stiffness fails to evaluate the stiffness of mechanism in this case.

_{ox}Then, this problem is solved by using the stiffness formulation. The final topology is shown in Fig. 8(b), where *ω* is set to 0.7. One can see that there are valid material connections between the input, output and fixed ports of the CPM. As shown in the last rows of Tables 1 and 2, the resulting **J** is close to **J**^{*}, e.g., *J*_{11} = 2.27. The values of its C-stiffness and stiffness are reasonable. Thus, the problem of C-stiffness is avoided and a valid final topology can be ensured by the stiffness formulation.

#### Solved by Using f_{2}.

*f*

_{2}shown in Eq. (9). According to this design problem, the kinematic requirement

*f*

_{2}is to maximize the workspace of the two translational freedoms and minimize the workspace of the rotational freedom, which is formulated as follows:

*f*

_{2}is combined with the stiffness requirement

*S*based on C-stiffness, the resulting final topology is similar to the topology in Fig. 8(a). The final topology obtained by using

*f*

_{2}and stiffness based

*S*is shown in Fig. 9.

*ω*is set to 0.7. The corresponding Jacobian matrix is described by Eq. (35). One can see that the workspace of the two translational freedoms (

*u*and

_{x}*u*) is much larger than that of rotational freedom

_{y}*θ*. The resulting values of

_{z}*J*

_{11}and

*J*

_{22}are positive, i.e.,

*f*

_{2}cannot control the sign of

*J*like

_{ji}*f*

_{1}does

### Design of 3DOF CPMs.

The second design problem of CPMs is to synthesize a 3DOF CPM with three symmetrically arranged compliant limbs. The positions of the CPM's input, output, and fixed points are shown in Fig. 10. Each design domain Ω* _{i}* is discretized by 50 × 80 finite elements for elastic analysis under the same boundary condition. The allowable amount of material is 20%. Since the kinematics of 3DOF CPM is complex, it is hard to predefine a desired Jacobian matrix for the optimization. The second form of kinematic requirement

*f*

_{2}shown in Eq. (9) will be used in the objective function. For the planar 3DOF CPMs, no freedom should be suppressed, i.e., the design objective is to maximize the workspace of these three freedoms.

*ω*= 0.7. Its corresponding Jacobian matrix

**J**is shown in Eq. (36). While the final topology obtained by using the stiffness formulation (

*ω*= 0.9) and its resulting Jacobian matrix

**J**is shown in Fig. 12 and Eq. (37), respectively. One can see that both of the C-stiffness and stiffness formulations are able to achieve valid topology in this example, whereas the topology obtained by the stiffness formulation is stiffer.

### Solving Benchmark Problems.

Although the Jacobian-based topology optimization method is developed for the CPMs, this method is applicable to the typical compliant mechanisms designed by the current topology optimization methods, e.g., displacement inverter, amplifier, and redirector.

#### Displacement Inverter and Amplifier.

*a*and output point

*o*are in the middle of the left and right sides, respectively. The whole design domain is discretized using 100 × 100 finite elements for elastic analysis. The material usage is restricted to 20%. The Jacobian matrix of the 1DOF compliant mechanism is a 3 × 1 vector, whose desired form is shown in Eq. (38). One can see that the desired output motion of this mechanism is in the direction of

*x*-axis, whereas the freedoms

*u*and

_{y}*θ*should be suppressed. Obviously, when only the element

_{z}*J*

_{11}of Jacobian matrix is considered, our objective function is equal to the formulation proposed by Chen and Wang (Eq. (3))

When GA^{*} is negative, e.g., GA^{*} = −3, the design problem is to synthesize a displacement inverter. The corresponding topology optimization problem is also solved by using the C-stiffness and stiffness formulations, respectively. Figure 14(a) shows the resulting topology obtained by using the C-stiffness formulation and setting the weight as *ω* = 0.5. Its corresponding Jacobian matrix is **J** = [−2.7, 0, 0]^{T}, i.e., the optimized GA of displacement inverter is −2.7. The C-stiffness at the input and output ports are 40.9 and 4.7 N/mm, while the stiffness at input and output ports are 6.4 and 0.7 N/mm, respectively. Figure 14(b) shows the resulting topology obtained by using the stiffness formulation and setting *ω* to 0.5. Its corresponding Jacobian matrix is **J** = [−2.8, 0, 0]^{T}. The C-stiffness at the input and output ports are 57.2 and 2.7 N/mm, while the stiffness at input and output ports are 41.3 and 2.0 N/mm, respectively. One can see that both of the C-stiffness and stiffness formulations can obtain valid topologies in this case.

When GA^{*} is positive, e.g., GA^{*} = 3, the design problem is to synthesize a displacement amplifier. It is found that the C-stiffness formulation results in the invalid topology shown in Fig. 15(a). Without material connected to the fixed ports, the Jacobian matrix of the resulting displacement amplifier is **J** = [1, 0, 0]^{T}. Both of the input and output C-stiffness are 116.1 N/mm, whereas both of the corresponding input and output stiffness are 4.3 × 10^{−7} N/mm. Figure 15(b) shows the resulting topology obtained by using the stiffness formulation and setting *ω* to 0.6. Its corresponding Jacobian matrix is **J** = [2.3, 0, 0]^{T}. The C-stiffness at the input and output ports are 79.5 and 6.2 N/mm, while the stiffness at input and output ports are 47.8 and 3.8 N/mm, respectively. This case shows again that the C-stiffness formulation fails to ensure the stiffness requirement.

#### Solved by Using the Artificial I/O Spring and MSE/SE Formulations.

^{−10}to 10

^{10}N/mm by ten times. On the contrary, the MSE/SE formulation (Eq. (40)) ensures the valid topology result shown in Fig. 16.

The computational expense of the proposed method is compared with the artificial I/O spring formulation by solving the design problem of displacement inverter in Sec. 5.3.1. It takes the artificial I/O spring formulation 47.59 s and 220 iterations to find out the solution, i.e., 216 ms per iteration. The C-stiffness formulation spends 181.5 s and 300 iterations to obtain the topology in Fig. 14(a), i.e., 605 ms per iteration. The stiffness formulation spends 209.37 s and 300 iterations to obtain the topology in Fig. 14(b), i.e., 698 ms per iteration. Obviously, the proposed method is more expensive than the artificial I/O spring formulation in computation. One reason for this is that the finite element with 12 nodal freedoms increases the computational expense of finite element analysis. The other reason is that all the freedoms of the output point *o* are considered by the proposed method, whereas the spring formulation only concerns the freedom *u _{x}*.

#### Displacement Redirector.

*a*is at the middle of the left side and causes two output displacement at ports

*b*and

*c*, respectively. The whole design domain is descretized using 100 × 100 finite elements for elastic analysis. The material usage is restricted to 20%. In this case, only the freedoms of the output ports

*b*and

*c*are considered, so the Jacobian matrix

**J**of the displacement redirector is a 2 × 1 vector, whose desired form is shown in the following equation:

The corresponding topology optimization problem is solved by using the C-stiffness and stiffness formulations, respectively. The optimization process was run for 300 iterations. Figure 18(a) shows the resulting topology obtained by using the C-stiffness formulation and setting *ω* to 0.5. Its corresponding Jacobian matrix is **J** = [1.7, −1.7]^{T}. The C-stiffness at the input and output ports are 41.2, 5.9, and 5.9 N/mm, while the stiffness at input and output ports are 10.7, 2.3, and 2.3 N/mm, respectively. Figure 18(b) shows the resulting topology obtained by using the stiffness formulation and setting *ω* to 0.5. Its corresponding Jacobian matrix is **J** = [1.5, −1.5]^{T}. The C-stiffness at the input and output ports are 58.1, 4.1, and 4.1 N/mm, while the stiffness at input and output ports are 40.1, 3.3, and 3.3 N/mm, respectively. This case shows that the proposed method is applicable to the compliant mechanisms with multiple output ports, and both of the C-stiffness and stiffness formulations can obtain valid topologies.

### Analysis of Mesh Independency.

This section is devoted to analyzing the mesh independency of the proposed method. The problem of 2DOF CPM, whose parameter settings except discretization are the same as the case in Fig. 5, is used to illustrate the mesh independency. By using three different element discretizations of 40 × 40, 200 × 200 and 300 × 300, we obtain the corresponding results shown in Fig. 19. The three topologies on the left side of Fig. 19 are solved by the C-stiffness formulation, whereas the other side contains the results of the stiffness formulation. One can see that the results are almost stable under mesh refinement or mesh coarsening. In other words, the proposed method is mesh independent, which is ensured by the filtering technique of the 99 line Matlab code.

## Conclusions

This paper presents a new stiffness evaluation based on the definition of stiffness for the Jacobian-based topology optimization method. The proposed stiffness formulation is compared with the C-stiffness formulation by using two synthesis problems of CPMs and three traditional benchmark design problems. The results show that both of the two formulations can achieve valid topologies in most of design cases. In some cases like displacement amplifier, the C-stiffness formulation even the artificial I/O spring formulation cannot obtain valid result, while the stiffness formulation gives an improved stiffness evaluation. Besides, the Jacobian-based topology optimization method shows a general applicability in multi-DOF CPMs and benchmark design problems.

According to the results, the topologies produced by the proposed method are easy to exhibit hinges, especially when the kinematic requirement is much higher than stiffness requirement. Relatively speaking, the stiffness formulation has a better performance than the C-stiffness formulation in avoiding hinges, e.g., the two topologies in Fig. 18. The strategies for alleviating these hinges will be addressed in our future works.

## Funding Data

National Natural Science Foundation of China (Grant Nos. 51275174, 51605166, U1609206, and 51675189).

Natural Science Foundation of Guangdong Province (Grant No. 2014A030313460).

Fundamental Research Funds for the Central Universities.