Topology optimization of a 6-DOF spatial compliant mechanism based on Stewart propotype platform

To improve the global stiffness and conveniently build a model of a compliant mechanism with spatial multiple degrees of freedom (DOF), the topology optimization method, combined with the isomorphic mapping matrix, is proposed in this paper for structure synthesis of a 6-DOF spatial compliant mechanism. By using the differential approximation method, the Jacobian matrix of the Stewart prototype platform is calculated as the isomorphic mapping matrix, and its eigenvalues and eigenvectors are considered. Combining the isomorphic mapping matrix with the solid isotropic material with the penalization topology optimization method, the topological model of the 6-DOF spatial compliant mechanism is constructed, and a topological structure of the 6-DOF spatial compliant mechanism is derived which has the same differential kinematic characteristics as the Gough–Stewart prototype platform. Piezoelectric actuators are mounted inside the topological structure during the three-dimensional printing manufacturing process, and its driver directions are in accordance with the driver configuration directions of the Gough–Stewart prototype platform. The effectiveness of the proposed method for topological structure synthesis of the 6-DOF spatial compliant mechanism is demonstrated through several numerical examples and experimental studies.


Introduction
Compared with traditional rigid mechanisms, compliant mechanisms have the characteristics of less parts, light weight and fewer joints, which result in reduced friction, wear, impact, vibration and noise in the movement of mechanisms, and the utility mechanism can improve the precision of the mechanisms, increase the reliability, and reduce the maintenance cost [1]. In the past two decades, the important components of a compliant mechanism, such as a flexible beam, flexure hinge and structural synthesis composited by these components, have become the focus of research [2][3][4][5]. In order to establish a mathematic model of a compliant mechanism composed of these components, an equivalent stiffness method called pseudo-rigid-body method had been proposed. For example, Yu and Zhu [6] proposed a 5R pseudo-rigid-body model for inflection beams in compliant mechanisms, and Zhu and Yu [7] proposed two types of pseudo-rigid-body models to simulate the large deflection of a flexible beam with an inflection point in different configurations. Based on finite elements analysis and on the principle of minimum potential energy, Jin et al. [8] presented a numerical method for analyzing the pseudo-rigid-body model of compliant mechanisms consisting of the finite elements and springs. The pseudo-rigid-body method has been widely used to design and analyze the process of compliant mechanisms [9][10][11][12][13].
Due to large deformation, a large-deflection slender beam is usually used for the structure synthesis of compliant mechanisms. However, it is difficult to construct the spatial multiple degrees of freedom (multi-DOF) compliant mechanism with complex geometric constraint conditions. How to realize multi-DOF for micro-displacement movement in micro-/nano-manufacturing fields is a problem that needs to be solved urgently. To solve this problem, the topology optimization method is proposed and applied to the structure synthesis of compliant mechanisms. The main advantage of topology optimization is that the optimal design is automatically suggested for prescribed design domains, boundary conditions and functional specifications. There is no need to pre-determine the number of links or the location of the flexural joints in the device [14]. In the past decades, many different techniques, such as the homogenization method [15], solid isotropic material with penalization (SIMP) method [16][17][18], level set method [19,20], evolutionary structural optimization (ESO) method [21,22] and others [23][24][25], have been developed. To a great extent, structure synthesis methods of spatial compliant mechanisms are faced with many complicated factors, especially the following two main issues.
(1) As the topology optimization methods are adapted to the structure synthesis of multi-DOF spatial compliant mechanisms, the topological results show that the topological structure is a hingeless configuration. The pseudo-rigid-body modeling method for a traditional flexible beam is not applicable to the spatial compliant mechanism derived by the topology optimization method. (2) The kinematic characteristics of spatial compliant mechanisms are very complex because we need to consider the three-dimensional (3D) constraints, forces and displacements simultaneously.
The structure synthesis of parallel mechanisms brings inspiration to the solution to the topology optimization problems for spatial compliant mechanisms. This approach can achieve the topology optimization theory, shape and dimension of spatial compliant mechanisms from the design specification directly, without using a known rigidbody mechanism. In other words, this approach adapts the Jacobian matrix of a parallel mechanism which has the mapping relationship between the operation space and the joint space, and it may be applied during the topology optimization process.
The Gough-Stewart platform with 6-DOF has been widely used in different mechatronic devices since 1965 [26]. Due to the square matrix with 6 9 6 dimensions, the forward/backward matrices of the Gough-Stewart platform can be calculated conveniently. Therefore, compared with the spatial structure of the Gough-Stewart platform, researchers have tried to use the flexible hinge replacement method to reconstruct a spatial compliant mechanism with the same kinematic characteristics as the Gough-Stewart platform [27,28]. In furtherance of this research, Zhao et al. [29] proposed a parallel compliant mechanism with six universal-prismatic-universal-revolute (UPUR) flexible joints to achieve a long-range, six-axis force sensor; the force Jacobian matrix is obtained by using screw theory in two cases of the ideal state, and the state of flexibility of each flexible joint is considered. However, the structure design method with hinge replacement cannot fundamentally solve the problem of overall structural stiffness reduction caused by the accumulation of flexible hinges. Furthermore, it is difficult to establish the equivalent mapping relationship by using the pseudo-rigid-body method between the rigid-elastic model composed of flexible hinges and the rigid-body model.
To overcome these shortcomings, a novel structure synthesis method considering the topology performance for 6-DOF compliant mechanisms is presented in this paper. Firstly, a new extension of the topology optimization technique with a six-inputs-six-outputs conditions is proposed using the SIMP method. Secondly, as a supplement to the traditional topology optimization method, the Jacobian matrix of the Gough-Stewart platform [30][31][32][33] is used as the isomorphic mapping matrix between the 6-DOF compliant mechanism and conventional Gough-Stewart platform, and then, this mapping matrix is applied to the topological process. By using this method, the kinematic characteristics of 6-DOF compliant mechanism are consistent with the conventional Gough-Stewart platform.
In the present paper, the main aim is to introduce the topology optimization theory combining the isomorphic mapping matrix with the SIMP method as an effectiveness design method for structure synthesis of spatial compliant mechanisms. The rest of this paper is organized as follows. The basic idea of the proposed method is firstly introduced in Sect. 2. Section 3 illustrates the Jacobian matrix of the conventional Gough-Stewart platform. Section 4 describes how to formulate, build and solve the model of topology optimization with the Jacobian matrix. Section 5 gives the experiment results and simulations to verify the proposed method. In addition, the vibration modal analysis with three orders is conducted. Finally, discussions and conclusions are presented in Sect. 6.

Basic idea of the proposed method
The general Stewart-Gough parallel manipulator consists of a base platform, a moving platform, and six limbs connected at six distinct points on the base platform and the moving platform respectively, as shown in Fig. 1a, and it is used to implement ultra-high-precision tasks under large external loads (e.g., flight simulator) mainly due to their high stiffness and high load carrying capacity. However, the backlash or joint clearance of the rigid hinge may reduce the accuracy. To overcome this shortcoming, the flexure hinge replacement method is used to reconstruct the structure of the Gough-Stewart platform, as shown in Fig. 1b. The stiffness of the mechanism obtained by the flexure hinge replacement method is obviously lower than that of the conventional parallel mechanism with rigid hinges, which is related to the vibration of the mechanism during movement. On the other hand, topology optimization is a kind of optimization method based on the minimization of mechanism flexibility, and the configuration of a mechanism derived by this method can improve the overall stiffness of the mechanism. This idea has been confirmed by our previous research [34]. Based on the relationship between the constraints and the motion characteristics of a parallel mechanism given by Ref. [35], there are no constraints on the moving platform. In order to avoid a driver configuration singularity, the multi-inputs in topology optimization should be presented according to the conventional driver configuration of the Gough-Stewart parallel mechanism. A design domain of a 6-DOF spatial compliant mechanism is shown in Fig. 2.
A cube structure with 200 mm Â 200 mm Â 200 mm is selected as the design domain of topology optimization, and the multi-inputs are set according to the driver configuration form of conventional Gough-Stewart platform.
The objective of the minimum compliance problem is to find a material element density distribution q e that maximizes the structural deformation under the prescribed support and loading conditions. The driver configurations of the Gough-Stewart platform are placed inside the design domain as the piezoelectric PbZrTiO 3 (PZT) actuator installation locations. According to the screw theory, the number of constraints on the outputs of the topological optimization domain equals zero, so we give three ranges of translation and three ranges of rotations according to the actual expected output values of the 6-DOF compliant mechanism, respectively. It is shown in Fig. 2.
Another problem we have to solve is how to build the kinematic model of a 6-DOF compliant mechanism. The premise of a kinematic model established by using the   pseudo-rigid-body method can be described by flexible joints. However, the hinge (or joint) and beam form does not exist in the topological structure, so the pseudo-rigidbody method is not suitable for the topological kinematic modeling. On the other hand, as we know, the Jacobian matrix represents the mapping relationship between the operation space and the joint space. Therefore, if we can use the Jacobian matrix of the corresponding parallel prototype mechanism in the design process of topology optimization, then the kinematics model of the parallel prototype mechanism can be used as the kinematics model of the compliant mechanism. This novel method can also be extended to the topology optimization design for a compliant mechanism with an arbitrary DOF.

Kinematics Jacobian matrix of the Gough-Stewart platform
Consider a 6-DOF conventional Gough-Stewart platform, as shown in Fig. 2a. Each link has a prismatic actuator (driver). The links are joined at the ends to the base platform (BP) and moving platform (MP) by a universal and a spherical joint, respectively. Let the connection points of the i-th link on the BP and MP be denoted by A i and C i , respectively. Also, let the two coordinate frames o À xyz f g and o 0 À uvw f gbe attached to the BP and MP at the points o and o 0 , respectively. The kinematics design problem involves determination of coordinates of the connection point for each link. Hence, six design parameters are associated with each link. Thus, a 6-links, non-redundant, parallel mechanism has 36 design parameters, which are required to be evaluated so that the design criterion is satisfied. But, since the number of these design parameters is quite large and governing equations are nonlinear, their values must be determined numerically.

Traditional solution method for a kinematics Jacobian matrix of the Gough-Stewart platform
The Gough-Stewart platform is a fully parallel kinematic linkage system that has major mechanical differences from typical serial link robots. Its closed kinematic chain and parallel linkage structure give it great rigidity and a high force-to-weight ratio. However, due to the lack of efficient algorithms for solving the kinematics equations, its potential application as a robotic manipulator is difficult to realize. In this paper, as the motion characteristics with a micro-/nanoscale of the 6-DOF compliant mechanism, we focus on the Jacobian matrix with displacement. A simplified algorithm was proposed by Liu et al. [36] to solve the forward kinematics of the Gough-Stewart platform by solving only three nonlinear simultaneous equations. Assuming that the output position vector of the Gough-Stewart platform is defined as  Applying the dot product of l i to both sides, Eq. (3) can be arranged with respect to l i as below where p is the displacement vector from the center of the base platform to the center of the moving platform.
Assuming that the initial frame parameters are represented by superscript 0 yields Following the same procedure for the other five kinematics links, the constraint function for the entire system can be derived as Topology optimization of a 6-DOF spatial compliant mechanism based on Stewart propotype platform  : ð6Þ To solve these equations, Tang et al. [37] proposed a numerical method based on an equivalent model for the Gough-Stewart platform forward kinematics problem by using the iteration method. In this method, the critical values of the allowable errors e i of a; b; c; d x ; d y ; d z are given, respectively. The iterations terminate when the posed difference Dl i is smaller than e i for all i ¼ 1; 2; . . .; 6.

Differential equivalent analysis method for a kinematics Jacobian matrix of the Gough-Stewart platform
The conditions of displacement output are set with the rotation value of the angle being [0, 1 9 10 -7 ] rad and the translation value of the displacement being [0, 1 9 10 -8 ] mm. Taking into account the differential kinematics characteristics of the conventional Gough-Stewart platform with rigid hinges, without loss of generality, the coordinate transformation homogeneous Eq. (1) can be rewritten by Substituting the initial coordinate parameters of C i into Eq. (7) yields As ; the differential equivalent method is adapted as follows where l i are the initial lengths of links; according to Fig. 1, they are equal to each other. Substituting Eq. (8) with the initial vector coordinate parameters into Eq. (9) yields Simplifying Eq. (10), we obtain Equation (11) can be rewritten as a matrix form Based on Fig. 2, the initial coordinates are shown in Table 1.
By substituting the parameters given by Table 1, the Jacobian matrix can be derived as  To ensure the homogeneity and decoupling of the Jacobian matrix, the orthogonal and right matrix deformation method (QR) is used to derive the orthogonal matrix, eigenvalues and eigenvectors.
The orthogonal matrix can be calculated by Matlab software, which yields And the eigenvalues E value and eigenvectors E vector are also calculated as : 4 Topology optimization model of the 6-DOF spatial compliant mechanism In this paper, a cube optimization design domain is selected as Fig. 2, and the Gough-Stewart prototype platform is placed inside this cube in order to determine the placement position of the six PZT actuators. During the topology optimization process, there are two key problems that need to be solved, that is, one is the how to integrate the Jacobian matrix with the topology optimization method, and the other is how to put the PZT actuators inside the topological structure when 3D printing technology is used to manufacture the topological structure. Our solution to the second problem is to put the PZT actuators in the slot when the reserved actuator slot is about to be closed, and the SIMP method to integrate the Jacobian matrix of the 6-DOF Gough-Stewart prototype platform will be used to solve the first problem.

Finite element method based on transfer characteristic between force/macro and strain/micro
The isomorphic Jacobian mapping matrix proposed as previously can be used to express the mapping relationship between macro inputs and outputs, that is, it can be ensured that the topological configuration after topology optimization can be consistent with the Gough-Stewart prototype platform, and then the model of the 6-DOF spatial compliant mechanism can be established. To realize macro/micro unification during the topology optimization process, three steps need to be considered as follows. Kinematic characteristics of the topological structure of 6-DOF complaint mechanisms can be consistent with the Stewart prototype platform. To realize macro/micro unification during the topology optimization process, three steps need to be considered as follows: Step 1: The stress and strain generated by the six driving forces provided by the PZT actuators need to be solved, and they are the initial conditions for solving the stress and strain between the element meshes.
Step 2: The finite element method is adapted during the topology optimization process, and the SIMP method will be used in this paper.
Step 3: The final stress/strain results of element meshes should correspond to the output displacements to achieve unified micro and macro motion characteristics, and the convergence condition is satisfied by the given Jacobian mapping matrix between the input forces and output displacements. Topology optimization of a 6-DOF spatial compliant mechanism based on Stewart propotype platform 1049

Transfer from the driving forces to stress/strain of element meshes
Under the action of external force, the object will produce deformation and internal force at the same time, and the deformation will be stopped when the internal force and the external force are kept in balance. Without loss of generality, according to the mesh grid with a hexahedral element in the topology optimization method, we assume that the driving force acts on a micro element with a hexahedral. The i-th driving force can be denoted as . . .; 6 and the macro expression between the driving force and displacement is satisfied as follows where u v w ð Þ T is the function of displacement with 5 is the whole elastic matrix.
And then, the micro-deformation in element meshes will be studied in order to connect with macro-displacement caused by driving forces. The coordinate transformation can deduce to locate a point into the evenly shaped triangle network. The projection plane and micro-deformation are shown in Fig. 3a, b, respectively.
By using the Taylor series expansion formula and neglecting the higher-order small quantity, the deformation coordinates values for the points projected on the X-Y plane are given by So, the length of line segment AB projected on the x axis is calculated, yielding Similarly, we have d dy ð Þ ¼ ov oy

Stress/strain transfer between element meshes
The basic idea of the finite element method is to divide the continuum into finite elements connected at the nodes, and then approximate the unknown field functions on the whole solution domain by approximating the assumed functions in each element. Therefore, the stress and strain transfer between the elements becomes an important factor during the topology optimization process. The equilibrium equation expressed by displacement is given by Compared with Eq. (17), the initial displacement caused by external force F i can be transferred by stress and strain inside of element meshes. In the initial stage of topology optimization, a reference domain X in R 3 should be chosen to allow by definition the applied loads and boundary conditions. In the topological design domain, the material distribution problem is formulated with optimal shape, design formula, and solved according to the load or support constraints. The reference domain X is discredited into n elements, and the relative density q e of material corresponding to the e-th element is a design variable. The objective function of these design variables is used for modifying the stiffness matrix so that it can be interpreted as a density of material. The penalized proportional stiffness model is given by where q x ð Þ is the density of the design function, E 0 ijkl represents the material properties of a given isotropic material, and E ijkl x ð Þ is the stiffness tensor which is a variable over all the design domain and is modified by the optimal algorithm; it is equivalent to the left side of Eq. (20) in the projective plane. V denotes the design domain which is satisfied with V 2 X, the design volume (not including the non-design area or fixed domain) is evaluated as R V q x ð ÞdV, and E ijkl x ð Þ is given by the density interpolates with constraint conditions, yielding According to the isotopic material properties and 3D topology optimization admissible, the power of p satisfies that p ! max 15 1 À m 0 7 À 5m 0 ; where m 0 is the Poisson ratio with the given stiffness tensor of the isotopic material. Generally, it is chosen to be greater than 3. The important problem is that the 6-DOF compliant mechanism includes the multiple force inputs and displacement outputs. The couple between each force input and displacement output should be considered at the initial condition of the topology optimization process as Eq. (15). The model of the 6-DOF spatial compliant mechanism can be given by two cases as following.
Case one. The Jacobian matrix is not considered during the topology optimization process q p e k e u ¼ KU, e ¼ 1; 2; . . .; n; where the design domain is divided into n elements by using a hexahedral structure, and K is the global stiffness matrix and satisfies K ¼ P 6 i¼1 P n e¼1 k e . q min is a lower bound of the density in order to prevent the possible singularity; generally, it is set q min ¼ 1 Â 10 À3 . By using the projective plane and the adjoint method, the function c q e ð Þ is rewritten by adding the zero function Because of the symmetry and consistency of the design structure, we do not use a subscript to represent the different force inputs and displacement outputs, andũ is any arbitrary fixed real vector.
Derivative c q e ð Þ with to respect q e , we obtain that It can be rearranged as whereũ satisfies the adjoint equation f T Àũ T k À Á ¼ 0. Assuming thatũ is equal to u, the differentiating of compliance c q e ð Þ is simple as Obviously, the sensitivity of density is negative for all elements.
Case two. Multiple-force inputs-displacement outputs with a Jacobian matrix are considered.
In case two, maximization of the fundamental eigenvalue is used as design functions for dynamic loads, yielding ; Topology optimization of a 6-DOF spatial compliant mechanism based on Stewart propotype platform 1051 wherex is the design variable denoting the density filter function and given byx i ¼ , H ij is a weight factor, K is the global stiffness matrix and k e is the element stiffness matrix. v ¼ v 1 ; v 2 . . .v n ½ T is element volume, v is the iterative convergence objective respect to the design volume V, g i is the eigenvalue given by Eq. (16) according to the Stewart prototype platform and u i is the eigenvector given by Eq. (15); M is the mass matrix.
With derivative g i with respect to q e , the sensitivity of the eigenvalue can be derived The eigenvector has been normalized with respect to the kinetic energy and satisfies u T i Mu i ¼ 1, and the sensitivities of eigenvalues opposed to the sensitivities of the compliance objective may take negative as well as positive values.
The sensitivity of the volume constraint vx ð Þ and compliance are given by derivation with respect to the design variable x e , yielding where F T À Ux ð Þ T K ¼ 0 satisfies the adjoint equation, and And then Eq. (32) can be rewritten by where u i is the nodal displacement of the i-th element, and k 0 i is definitely positive.
The sensitivity of the objective function is given by According to the Jacobian matrix J D derived by QR decomposition of the Stewart prototype platform, without loss of generality, the force inputs are assumed independent from the design variables, which yields   where m is the constant value with moving constraint, and m ¼ 0:1À0:3. 1 is the damping factor, and 1 ¼ 0:4À0:5.
D ðkÞ e is given by the expression where K k ð Þ is the Lagrange multiplier with the volume constraint at the k-th step iteration, and V e is the element volume of k-th step iteration.
The iteration procedure is shown in Fig. 4.
5 Topological structure of the 6-DOF spatial compliant mechanism In this section, the configuration of the actuators in 3D space will be considered for manufacture, and then the way that the six actuators are embedded in the topological structure will be determined. The prototype structure of the Gough-Stewart platform is adapted and shown as Fig. 2, and the six directions according to the prototype structure are accepted as six force inputs. Six holes along the direction of the installation of piezoelectric actuators are set to the non-design area. When the height of the 3D printing process is parallel to the piezoelectric actuator, the 3D printing process will be terminated, and the PZT actuator is packaged in a topological structure simultaneously. The design domain is a cube of dimension 200 mm 9 200 mm 9 200 mm discredited with 2000 eight-node elements with a hexahedral structure. The composition material of the 6-DOF spatial compliant mechanism is defined with a Young's modulus of 2.06 9 10 -5 F/mm 2 and Poisson's ratio of 0.3. The 6 force inputs are assumed as F in1 ¼ F in2 ¼ Á Á Á ¼ F in6 ¼ 180 N, the volume restriction is 20% of the given design domain, and for the elastic workpiece, a value of 2 Â 10 6 mN/m is assigned to the spring constant.

Topological structure of case study one
In case one, we only consider the multiple force inputs and displacement outputs using the SIMP topology optimization method, and the constraints of the displacement outputs are limited in 1 Â 10 À5 mm or 1 Â 10 À6 rad scope. The six PZT actuators are installed along the same directions with the drivers of the Stewart prototype platform in order to avoid singularity.
The optimization results for the 6-DOF spatial compliant mechanism are presented below. Figure 5 includes the optimal topologies. Figure 6 displays the convergence process, and the stress nephogram is shown in Fig. 7.
Extracting displacement data from each iterative step, the displacement values of the moving platform caused by the overall elastic deformation of the material are shown as Table 2.

Topological structure of case two
In case two, we also only consider the multiple force inputs and displacement outputs using the SIMP topology optimization method, and the constraints of the displacement outputs are limited in 1 Â 10 À5 mm or 1 Â 10 À6 rad scope. The six PZT actuators are installed along the same directions with the drivers of the Stewart prototype platform in order to avoid singularity. The optimization results for the 6-DOF spatial compliant mechanism of case two are presented below. Figure 8 includes the optimal topologies. Figure 9 displays the convergence process, and the stress nephogram is shown as Fig. 10.
By extracting displacement data from each iterative step, the displacement values of the moving platform caused by the overall elastic deformation of the material with case two are shown as Table 3.
Compared the displacement and rotational values between case one and two, we can find that the displacements and rotational angles order of magnitude are close to each other in two cases, and it is decided by the initial Rotation around the z axis (rad) 5.1076 9 10 -8 1.5998 9 10 -7 9.6931 9 10 -4 -1.4937 9 10 -9 4.6365 9 10 -10 1.0958 9 10 -12 Fig. 8 Topologies for the 6-DOF spatial compliant mechanism under the different iterative steps and the topological structure (case two). a Topological structure at the iteration 1 step. b Topological structure at the iteration 8 step. c Topological structure at the iteration 12 step. d Topological structure of the 6-DOF spatial compliant mechanism convergence conditions of the topology optimization process, and the values of the displacements and rotational angles are consistent with the given parallel prototype manipulator with Jacobian mapping matrix.

Vibration modal analysis with the first four orders
The vibration modal is the most important factor in precision manufacturing, and it can provide accurate information for the location of sensors. By using the Optistruct software, the first four orders of the vibration modal are proposed with the same load input conditions and boundary constraint conditions. The first four orders of the vibration modal of case one and case two are shown as Figs. 11 and 12, respectively. The results shows that the natural frequency between case one and case two is almost close to each other. This highly natural frequency can suppress the vibration of the proposed spatial compliant mechanism and reduce the uncertainty of the precision positioning.

Discussions
In this paper, an isomorphic Jacobain mapping matrix derived by the differential kinematic method with the Gough-Stewart prototype platform is proposed. Compared with the Gough-Stewart prototype platform, a differential approximation method is adapted to calculate and construct the Jacobian matrix of the Gough-Stewart prototype platform. Its eigenvalues and eigenvectors are derived to build the topological model of 6-DOF spatial compliant mechanisms. By using this isomorphic mapping matrix combined with the topology optimization method, we can easily establish the kinematic/dynamic model of 6-DOF spatial compliant mechanisms with topological structure which cannot be realized by Fig. 9 Typical convergence history for 6-DOF spatial compliant mechanism (case two) Fig. 10 Stress nephogram of the topological structure of 6-DOF spatial compliant mechanism with the maximum stress 0.3739, and minimum stress 4.780 9 10 -5 , respectively (case two) the pseudo-rigid-body method. Piezoelectric actuators are mounted inside the topological structure during the 3D printing manufacture process, and its driving directions are in accordance with the driver directions of the Gough-Stewart prototype platform. The SIMP method is used for the solving topological model which combined the multiple force inputs and multiple displacement outputs with the isomorphic mapping Jacobian matrix. The sensitivity of volume, penalty and compliance with respect to the design function are proposed. In summary, we can draw the following conclusions.
(1) The topological model of a 6-DOF spatial compliant mechanism is constructed by the combination of a SIMP optimization model and the Jacobian isomorphic mapping matrix of the Gough-Stewart prototype platform, and the topological structures in two cases are derived in this paper. (2) The differential kinematics, including displacements and rotational angles of the moving platform, are calculated with the two cases, and simulation results are given to analyze the isomorphic kinematic characteristic between the topological structures of a 6-DOF spatial compliant mechanism with the Gough-Stewart prototype platform. Simulation results show that the topology optimization method proposed in this paper is effective. (3) The topological structure of the 6-DOF spatial compliant mechanism with isotropic material is repaired and imported by using Hypermesh @ software for executing finite element static analysis, and then the sign and order of magnitude are consistent in all directions. It is validated that the integer part of the experimental results are sufficiently accurate to represent the design method of topology optimization with a Jacobian isomorphic mapping matrix derived from the prototype parallel mechanism for the structure synthesis of spatial compliant mechanisms.