Modeling and Analysis of a Novel 3R Parallel Compliant Mechanism

: This paper presents and investigates a new three-rotation (3R) parallel compliant mechanism that uses compliant rods to achieve three rotations. The mechanism is designed for use in pointing devices or as a spatial parallel manipulator. The mobility analysis is based on the Cosserat rod model and Lagrangian dynamics equations. The dynamics equations are then effectively solved using the back-propagation neural network and chaos-enhanced accelerated particle swarm optimization. After studying the mobility of the moving platform, a simpliﬁed model is proposed and used for kinematic analysis. The analysis of motion includes discussions on forward kinematics, inverse kinematics, singularities


Introduction
A compliant mechanism (CM) is a device that employs the deformation of its compliant elements to transform or transmit loads and motions.One of the significant benefits of CMs is their ability to reduce costs and improve performance by minimizing the number of parts required, shortening assembly times, simplifying manufacturing processes, and decreasing wear, friction, and noise.These advantages make CMs widely used in different applications, such as compliant assembly systems, vibratory bowl feeders, and highprecision manipulators.
Recently, CMs have drawn much attention from researchers, and have been used in precision position systems, metrology instruments, MEMS/NEMS devices and other fields.
For precision position systems, Schitter [1] presented a novel design of a scanning unit for atomic force microscopy (AFM).Kim [2] designed a new AFM with a 2D plane CM for the establishment of a standard technique of nano-length measurement.Minh [3] designed a decoupled 6-DOF compliant parallel mechanism.
For metrology instrumentation, Khalid [4] applied a CM in 3D coordinate measurement.Jin [5] and Hao [6] used CMs as force/moment sensors.Hansen [7] and Gao [8] designed displacement/acceleration sensors with CMs.For bio-medical/health devices, Sung [9] designed an ankle rehabilitation with CMs.Chen [10] applied CMs to a bodygravity compensation device.Awtar [11] presented a new minimally invasive surgical tool design paradigm that enables enhanced dexterity.
For MEMS/NEMS devices, Liew [12] designed a bulk-micromachined CMOS micromirror.In [13,14], Olfatnia applied electrostatic comb-drive actuators to drive a large-range dual-axis micro-stage.Aten [15] presented a self-reconfiguring metamorphic nanoinjector for injection into mouse zygotes.For compliant space mechanisms, Fowler [16] studied compliant space mechanisms applied to a frontier and a space-pointing mechanism.Throughout all the CMs in these applications, most of them use plate-shape compliant elements which mainly provide deformation in a plane.In this paper, compliant elements that can provide deformation in three-dimensional space will be used.Thalman [17] proposed an approach to design Flexure Pivot time bases.
The main challenge in CM application and analysis is the modeling of compliant element deformation.Accurate modeling of CMs is consistently desired to provide quick insight into the effect of material properties, displacements, geometrical parameters, and loads on the performance characteristics of CMs.The key issue in accurate modeling is how to describe the coordinates of any point on the compliant rod.There are many emerging modeling methods [18][19][20][21][22][23][24][25][26][27][28][29][30][31] for compliant elements.
For energy methods [22,23], Sen [22] used a closed-form nonlinear model when analyzing the constraint characteristics of symmetric spatial beams.In [23], Awter proposed the nonlinear strain energy formulation for two-dimensional beam flexures.
For numerical solutions [24][25][26][27][28][29][30][31], the elliptical integration method is commonly used to address large geometrical nonlinearity [24].Saxena and Kramer [25] solved these equations with numerical integration using Gauss-Chebyshev quadrature formula.Howell et al. [26] first proposed the pseudo-rigid body method (PRBM), which can help simplify the design and analysis of compliant mechanisms.Pietro [27] used the pseudo-rigid body method to synthesize a sub-optimal lumped compliance solution.More recently, other PRBM with prismatic (P) joints [28] was proposed to incorporate the elastic extension effect and even the effects from shearing deflection and cross-section changes for Timoshenko beams.With the Cosserat rod model widely used in computer graphics, some researchers [29,30] have tried to apply it to robotics.Weak-form Cosserat rod models can be formulated and solved using a finite-element [31], finite-difference, or discrete-differential geometry approach, and will be used in this study.
In the above studies, the exact solution of the compliant rod requires a huge amount of computation.This is not conducive to real-time precise control of CMs.How to balance the contradiction between accuracy and calculation amount is another challenge in compliant rod modeling.Therefore, this paper will combine neural networks to simplify the mapping relationship between the input and output variables of CMs, and discuss the modeling and kinematics analysis of the novel compliant mechanism.
The rest of this paper is organized as follows.Section 2 introduces the novel compliant mechanism and its structure.Section 3 introduces the theoretical models used in analyzing compliant rods, and discusses the mobility of the novel mechanism.In Section 4, the mechanism is simplified to a rigid model, and kinematic variables are calculated based on the simplified model.Next, the singularity is discussed with force analysis, and the workspace without singularity positions.Finally, a prototype is manufactured and used to verify the mobility analysis and simplified model.

Design of the 3R Parallel Compliant Mechanism
In this section, the construction of the novel 3R parallel compliant mechanism prototype structure and its analysis model will be proposed.

CAD Model
The basic design concept of the novel mechanism is illustrated in Figure 1.The prototype comprises a moving platform, a base, three flexible bars, and three drive cables.Both the moving platform and the base are designed as equilateral circles.The ends of each compliant rod are fixed at the edge of the moving platform and the base, respectively.The tangential direction of the rod's end is parallel to the moving platform.The three compliant rods are connected at their midpoint by a restraining ring, and the tangential directions at the rods' midpoint are identical to each other.One end of the drive cable is fixed to the edge of the moving platform, and the other end is fastened to the winch through a guide hole at the base's edge.A servo motor drives the winch to rotate, and changes in the length of the drive cables result in the movement of the moving platform.
totype comprises a moving platform, a base, three flexible bars, and thre Both the moving platform and the base are designed as equilateral circle each compliant rod are fixed at the edge of the moving platform and the bas The tangential direction of the rod's end is parallel to the moving platf compliant rods are connected at their midpoint by a restraining ring, and directions at the rods' midpoint are identical to each other.One end of th fixed to the edge of the moving platform, and the other end is fastene through a guide hole at the base's edge.A servo motor drives the winch changes in the length of the drive cables result in the movement of the mo

Mechanism Analysis Model
To analyze the novel mechanism's movement characteristics, the po ture of the mechanism will be calculated for any input parameter (three length li i = 1, 2, 3).As shown in Figure 2, the base coordinate is established origin point is set at the restraining ring's center, and X-axis is in the dire origin point to the middle point of line segment A1B1.Z-axis is in the dire base platform's center to the moving platform's center, and Y-axis is determ right-hand rule.Three vertices of the moving platform are defined as A1, same as the base platform' vertices to B1, B2 and B3.The length of driven c A1 is defined as l1, the same as l2 and l3.
Next, force analysis will be discussed, which is the hardest part in cal pliant mechanism's kinematics.To analyze the force in the mechanism, th separated into the upper part and lower part.The upper part (shown in tains the moving platform and compliant rods' segments from the constrai moving platform, and the remaining parts belong to the lower part (show In analyzing a compliant rod in the upper part, one rod's end is fixed at th and the another end exerts forces form the moving platform pfAi and drive 2, 3), as in the lower part.As mentioned in Section 2.1, the moving platfor platform have the same structure dimension.When the mechanism is rot axis in X-Y plane at 180 degrees, the mechanism will keep the same shap This means that the upper part and lower part can be replaced by each o force pFA1 has equal numerical value and opposite direction compared to pF pFA2 and pFA3.

Mechanism Analysis Model
To analyze the novel mechanism's movement characteristics, the position and posture of the mechanism will be calculated for any input parameter (three driven cables' length l i i = 1, 2, 3).As shown in Figure 2, the base coordinate is established as follows: the origin point is set at the restraining ring's center, and X-axis is in the direction from the origin point to the middle point of line segment A 1 B 1 .Z-axis is in the direction from the base platform's center to the moving platform's center, and Y-axis is determined using the right-hand rule.Three vertices of the moving platform are defined as A 1 , A 2 and A 3 , the same as the base platform' vertices to B 1 , B 2 and B 3 .The length of driven cable from B 1 to A 1 is defined as l 1 , the same as l 2 and l 3 .
Machines 2023, 11, x FOR PEER REVIEW 4 of 18 0( 1,2,3) Since the cables can only bear tension, the force Fl1 is coincident with cable direction, the same as Fl2 and Fl3.Next, force analysis will be discussed, which is the hardest part in calculating a compliant mechanism's kinematics.To analyze the force in the mechanism, the mechanism is separated into the upper part and lower part.The upper part (shown in Figure 2b) contains the moving platform and compliant rods' segments from the constraining ring to the moving platform, and the remaining parts belong to the lower part (shown in Figure 2d).In analyzing a compliant rod in the upper part, one rod's end is fixed at the origin point, and the another end exerts forces form the moving platform p f Ai and driven cable f i (i = 1, 2, 3), as in the lower part.As mentioned in Section 2.1, the moving platform and the base platform have the same structure dimension.When the mechanism is rotated along any axis in X-Y plane at 180 degrees, the mechanism will keep the same shape and position.This means that the upper part and lower part can be replaced by each other.Thus, the force p F A1 has equal numerical value and opposite direction compared to p F B1 , the same as p F A2 and p F A3 .
Since the cables can only bear tension, the force F l1 is coincident with cable direction, the same as F l2 and F l3 .
For the moving platform, this just bears the moment of the force from three compliant rods at the rods' ends.When the system reaches equilibrium, the amount of external force in the moving platform is zero.
According to the geometric constraint, the direction of the rods' ends is collinear to the platform's triangle midline, and the direction of the other end is collinear to the Z-axis: where d 3ij (i = A, B; j = 1, 2, 3) is the normal of the rods' cross-section, as explained in Section 3.2.
According to the structure, the length between the rods' ends is equal to the triangle's side length: Assuming there is a function between forces and position at the compliant rod's end: {x, y, z, q 1 , q 2 , q 3 , According to the lengths of driven cables: By solving Equations ( 1)-( 7) above, the moving platform's position and pose can be obtained.However, the relationship between the compliant rod's deformation and force is unknown.On the other hand, the equations are nonlinear and huge, and almost cannot be solved with analytical solutions.To solve these problems, three methods will be applied here and introduced in next section.

Theoretical Model
In this section, three methods will be introduced to analyze the mechanism's movement characteristics.At first, the Cosserat rod model will be applied to the analysis of the compliant rod.The relationship between force and deformation of a compliant rod will be investigated with the Cosserat rod model and Lagrangian dynamics equation.Next, the relationship will be trained as a neural network, and the neural network will be used to calculate the mechanism's kinematics.To solve a large number of nonlinear equations in the mechanism's kinematics, a chaos-enhanced accelerated particle swarm optimization method will be used.

Method 1: Cosserat Rod Model
This section provides a brief introduction to the Cosserat rod model of elastic rods.An elastic rod can be visualized as a long thin deformable body.If a rod's length is much greater than its radius, the rod's continuous configuration can be characterized by the centerline r(w) = {r x (w), r y (w), r z (w)} T , where r(w) assigns a position in space to each line parameter value w ∈ [0, 1].An orthonormal frame with basis vectors {d 1 (w), d 2 (w), d 3 (w)} is attached to every point on the curve, such that d 1 and d 2 span the plane of the rod's cross-section, and d 3 is the normal to the cross-section (as shown in Figure 3).The dynamic equilibrium configuration of an elastic rod is charac point of the Lagrangian.Details can be found in Goldstein [32].
where the gi∈{rx, ry, rz, q1, q2, q3, q4} are coordinates and Fe are external f T and V = vs.+ Vb + Vt are the kinetic energy and the potential energy the total dissipation energy; γt is the translational internal friction coef tational internal friction coefficient; Cp and Cq are holonomic constrain grangian multipliers; Bk and Bk 0 are a constant skew-symmetric matri vector at point r; r To relate the directions d i to the reference frame, quaternions are chosen as a representation of rotation.Quaternions can be present as a quadruple q = {q 1 , q 2 , q 3 , q 4 } T with q i ∈ R.Only unit quaternions represent pure rotations, thus the q i is not independent but coupled by the constraint q = 1.The directions d i in terms of the quaternion q are given by: 1 − q 2 2 − q 2 3 + q 2 4 , 2(q 1 q 2 + q 3 q 4 ), 2(q 1 q 3 − q 2 q 4 ) T d 2 = 2(q 1 q 2 − q 3 q 4 ), −q 2 1 + q 2 2 − q 2 3 + q 2 4 , 2(q 2 q 3 + q 1 q 4 ) For elastic rods, the potential energy consists of three parts, namely the energy versus the stretch deformation, energy V b of the bending deformation, and energy V t of the torsional deformation.However, the bending deformation and torsional deformation are similar in expression, and can be integrated to V b,t .The three energies are given as: where K s = E s πr r 2 is the stretching stiffness; E s is the stretch modulus; the û k conform to the intrinsic bending and torsion of the rod; K k is the stiffness tensor (K 1 = K 2 = 0.25 Eπr r 2 , K 3 = 0.5 Gπr r 2 ); E is the Young's modulus governing the bending resistance; G is the shear modulus governing the torsional resistance; and r r is the radius of the rod's cross-section.The dynamic equilibrium configuration of an elastic rod is characterized as a critical point of the Lagrangian.Details can be found in Goldstein [32].
where the g i ∈{r x , r y , r z , q 1 , q 2 , q 3 , q 4 } are coordinates and F e are external forces and torques; T and V = vs.+ V b + V t are the kinetic energy and the potential energy of the system; D is the total dissipation energy; γ t is the translational internal friction coefficient; γ r is the rotational internal friction coefficient; C p and C q are holonomic constraints; λ and µ are Lagrangian multipliers; B k and B k 0 are a constant skew-symmetric matrix; r' is the tangent vector at point r; . r is the derivative of r .

Method 2: Back-Propagation Neural Network
The back-propagation (BP) network is a widely used multi-layer feed-forward network with error inverse propagation.It is generated using the BP algorithm, which is a type of supervised learning algorithm.The BP algorithm continuously calculates the network's weights and the error's trend in the direction that minimizes the error.By doing so, it gradually reduces the error until the sum of squares of errors is minimal.The weights and the error's trend are positively correlated with the network's errors and are transmitted in reverse to every layer.The BP network can learn and store numerous input-output mapping relationships without a mathematical expression.It can be used to map the relationship between forces and the position and direction at a compliant rod's end.
A BP network consists of an input layer, an output layer, and multiple hidden layers.In this study, positions and quaternions of a compliant rod are taken as the input layer, and the forces at the rod's end are taken as the output layer.This BP network's structure is shown in Figure 4, and can replace the function F f in Equation (6).To train the neural network, 6.4 × 10 7 sets of data are calculated and used as the samples.
Machines 2023, 11, x FOR PEER REVIEW 7 of 18 A BP network consists of an input layer, an output layer, and multiple hidden layers.In this study, positions and quaternions of a compliant rod are taken as the input layer, and the forces at the rod's end are taken as the output layer.This BP network's structure is shown in Figure 4, and can replace the function Ff in Equation (6).To train the neural network, 6.4 × 10 7 sets of data are calculated and used as the samples.

Method 3: Chaos-Enhanced Accelerated Particle Swarm Optimization Method
Kennedy and Eberhart developed the standard Particle Swarm Optimization (PSO) in 1995.Extensive studies have demonstrated that PSO is highly efficient for optimizing many problems.However, for highly multimodal problems, it may suffer from premature convergence.The standard PSO uses g (the current global best) and x (the individual best) to increase the diversity in quality solutions.To simplify the algorithm and accelerate convergence, the global best can be used exclusively.In 2008, accelerated particle swarm optimization (APSO) [33] was proposed, which generates the velocity vector using a simpler formula.
where r random is drawn from N (0,1).
During each iteration, the particle's position is updated by The randomization term αr allows the system to escape local optima, and r can be solved by a probability distribution.Typically, a value of β ranging from 0.2 to 0.7 is adequate for most situations.Another enhancement to APSO involves reducing randomness as iterations progress.This can be achieved using the follow function: where 0 < δ < 1 is an annealing-like parameter whose value can be taken as 0.7 to 0.9 for most applications; t is the t-th iteration.
There is no requirement to maintain a constant value of β in standard APSO.However, a varying β may provide advantages and lead to faster convergence of the algorithm.Since all chaotic maps are normalized, a chaos-enhanced APSO (known as chaotic APSO) uses the chaotic map to tune the parameter β.In [34], after comparing various chaotic APSO variants, the sinusoidal map is deemed the best chaotic APSO, and is expressed as follows:

Numerical Computation
In this section, all these three methods will be integrated.There are two steps in this integration.The first is to analyze a single compliant rod with the Cosserat rod model and BP network.The second is to build the compliant mechanism's kinematical model.
The purpose of the first step is to establish the relationship between forces, positions and quaternion of a compliant rod.First, the rod will be simplified to a Cosserat rod model.By solving the Lagrange dynamic equations, the rod's position and quaternion can be calculated under forces.In addition, all the material parameters for the compliant rod are shown in Table 1.To facilitate the solution of the motion's equations, the rod is discretized into elements.First, the centerline of the rod is expressed as a chain of The centerline elements r i+1 − r i may differ in size.The orientations of the centerline elements are represented by quaternions q j , j ∈ [1, N − 1], as illustrated in Figure 2. The discrete spatial derivative r i of the centerline and q j are obtained as Next, a large number of the rod's positions and quaternions can be calculated as samples in the training BP network.Taking the positions and quaternions as the input parameters and the forces and moments as the output parameters, the BP network is trained with 10 hidden neurons and the Levenberg-Marquardt algorithm.The second step is to build the compliant mechanism's kinematical model.The model is introduced in Section 2.2.However, as it is hard to solve these equations with an analytical solution, the chaos-enhanced accelerated particle swarm optimization method is used here.As shown in Figure 5, the numerical computation can be separated into serval steps: (a) taking the mechanism as a rigid mechanism, the position and quaternions can be calculated.The result can be taken as the initial value.Next, all particle swarms can be initialed. (

Calculate fitness value
Renew speed and location of particle

Results and Summary
For numerical computation, the radii of the moving platform and base platform is set as 65 mm.The range of the driven cable's length is set to be 65 mm-185 mm.To reduce the calculation amount and ensure the result's reliability, 124 sets inputs are chosen ran-

Results and Summary
For numerical computation, the radii of the moving platform and base platform is set as 65 mm.The range of the driven cable's length is set to be 65 mm-185 mm.To reduce the calculation amount and ensure the result's reliability, 124 sets inputs are chosen randomly from the range of the driven cable's length.
Next, all these positions will be analyzed.The relative motion between a moving platform and base platform will be investigated.According to the moving platform and base platform positions, the relative rotating center (RRC) and the relative rotating radius (RRR) will be defined here.As shown in Figure 6, cml 1 is the normal line of the moving platform via the moving platform's center point O A , the same as cml 2 to base platform.Next, the common vertical line of cml 1 and cml 2 could be found, and the middle point of common vertical line is set as O e .Point O e is approximately seen as the relative rotational center.Next, the distance between O e and the moving platform's center point O A is approximately seen as its relative rotational radii l e1 , the same as l e2 to the base platform.The mean of l e1 and l e2 are seen as the mechanism's RRR.approximately seen as its relative rotational radii le1, the same as le2 to the base platform.
The mean of le1 and le2 are seen as the mechanism's RRR.According to the position results, a set of RRC and RRR can be calculated.The results are shown in Figure 7. From Figure 7c, the RRC is close to the origin point O, and the mean distance is 2.7 mm. On the other hand, the RRR ranges from 54 mm to 57 mm, and the mean is 56 mm.It is obvious that 2.7 is far less than 56 mm, which means that origin point O could be taken as the mechanism's RRC.In addition, the RRR can be seen as a constant.A constant rotate radius and a fixed rotational center means that the moving platform and base platform move along a spherical surface.In summary, the moving platform has three spatial rotational degrees of freedom.According to the position results, a set of RRC and RRR can be calculated.The results are shown in Figure 7. From Figure 7c, the RRC is close to the origin point O, and the mean distance is 2.7 mm. On the other hand, the RRR ranges from 54 mm to 57 mm, and the mean is 56 mm.It is obvious that 2.7 is far less than 56 mm, which means that origin point O could be taken as the mechanism's RRC.In addition, the RRR can be seen as a constant.A constant rotate radius and a fixed rotational center means that the moving platform and base platform move along a spherical surface.In summary, the moving platform has three spatial rotational degrees of freedom.approximately seen as its relative rotational radii le1, the same as le2 to the base platform.

(a)
The mean of le1 and le2 are seen as the mechanism's RRR.According to the position results, a set of RRC and RRR can be calculated.The results are shown in Figure 7. From Figure 7c, the RRC is close to the origin point O, and the mean distance is 2.7 mm. On the other hand, the RRR ranges from 54 mm to 57 mm, and the mean is 56 mm.It is obvious that 2.7 is far less than 56 mm, which means that origin point O could be taken as the mechanism's RRC.In addition, the RRR can be seen as a constant.A constant rotate radius and a fixed rotational center means that the moving platform and base platform move along a spherical surface.In summary, the moving platform has three spatial rotational degrees of freedom. (a)

Kinematics Analysis of 3R Parallel Compliant Mechanism
From the conclusion in Section 3.5, the novel compliant mechanism is a device whose moving platform can rotate around the center point in any direction.Thus, it can be considered to be a rigid parallel mechanism with three rotational degrees of freedom.As shown in Fig. 8, three flexural limbs are replaced by a spatial ball hinge at a center point.On the other hand, the driven cable can only bear the tension.That means the cable can just drive the moving platform rotate counterclockwise along the Z-axis.To drive the moving platform to rotate clockwise at the initial position, there must be a prerotation at the initial position.As shown in Figure 8c, the position OA 1 OA 2 OA 3 is the moving platform in assembly-complete state without pre-rotation.Next, a pre-rotation is applied to the mechanism.Then, the mechanism reaches the initial position A 1 A 2 A 3 with a counterclockwise degree θ 0 .

Kinematics Analysis of 3R Parallel Compliant Mechanism
From the conclusion in Section 3.5, the novel compliant mechanism is a device whose moving platform can rotate around the center point in any direction.Thus, it can be considered to be a rigid parallel mechanism with three rotational degrees of freedom.As shown in Fig. 8, three flexural limbs are replaced by a spatial ball hinge at a center point.On the other hand, the driven cable can only bear the tension.That means the cable can just drive the moving platform rotate counterclockwise along the Z-axis.To drive the moving platform to rotate clockwise at the initial position, there must be a pre-rotation at the initial position.As shown in Figure 8c, the position OA1OA2OA3 is the moving platform in assembly-complete state without pre-rotation.Next, a pre-rotation is applied to the mechanism.Then, the mechanism reaches the initial position A1A2A3 with a counterclockwise degree θ0.

Direct Kinematics
The direct kinematics analysis of a parallel mechanism aims to calculate moving platform attitude angles and positions according to given joint variables.For the novel mechanism, the three driven cables' lengths li (I = 1, 2, 3) are joint variables, and the three attitude angles are the moving platform's pose.In this study, Z-Y-X Euler angles are applied to express the moving platform's pose, as the sensor in the experiment used it.The Z-Y-X Euler angles mean that the moving platform rotates θz degrees along the Z-axis of the fixed frame first, acquiring a new frame X′Y′Z′.Then, the moving platform rotates θy degrees along the Y′-axis of the new frame X′Y′Z′, acquiring another new frame X′′Y′′Z′′.After that, the moving platform rotates θx degrees along X′′-axis of the new frame X′′Y′′Z′′.Next, the rotation matrix of the moving platform can be expressed as follows.
where Ci and Si represent cos(θi) and sin(θi), and i = x, y, z.

Direct Kinematics
The direct kinematics analysis of a parallel mechanism aims to calculate moving platform attitude angles and positions according to given joint variables.For the novel mechanism, the three driven cables' lengths l i (I = 1, 2, 3) are joint variables, and the three attitude angles are the moving platform's pose.In this study, Z-Y-X Euler angles are applied to express the moving platform's pose, as the sensor in the experiment used it.The Z-Y-X Euler angles mean that the moving platform rotates θ z degrees along the Z-axis of the fixed frame first, acquiring a new frame X Y Z .Then, the moving platform rotates θ y degrees along the Y -axis of the new frame X Y Z , acquiring another new frame X"Y"Z".After that, the moving platform rotates θ x degrees along X"-axis of the new frame X"Y"Z".Next, the rotation matrix of the moving platform can be expressed as follows.
where C i and S i represent cos(θ i ) and sin(θ i ), and i = x, y, z.At the initial position, the vertexes' coordinates of the moving platform and base platform can be expressed as follows.
Then, the vertexes' coordinate of the moving platform at any position can be written as follows.
According to the given driven cables' lengths l i , forward kinematics equations can be obtained.
By solving Equation ( 19), three Euler angles can be obtained.

Inverse Kinematics
The inverse kinematics of the novel mechanism aims to solve the three driven cable's lengths according to a given moving platform pose.According to Equations ( 16)-( 19), the inverse kinematics can be solved.

Results and Summary
A mechanism is driven by actuators to move.When a mechanism experiences changes in its kinematical or dynamic performance at certain positions, it may become stuck in dead-point positions, lose balance stability, or alter its degree of freedom.Consequently, the mechanism cannot transfer movements or be controlled properly.These positions are known as singularities.In the case of a novel mechanism, the movement of the moving platform depends on compliant rods and driven cables.Thus, the mechanism will only encounter singularity positions when the deformation of the compliant rods disappears or the tension in the driven cables disappears.

Compliant Rods' Deformation Disappearance
Compliant rods' deformation includes bending, stretching, and twisting.In this mechanism, the compliant rods mainly bear radial forces, so stretching will not be discussed here.The key points are bending and twisting deformation.When rod bending disappears, rods will recover to become straight rods.However, whatever position the moving platform is in, the compliant rod will always bend.As a result, the rod bending will not disappear.
For rod twisting deformation, when the mechanism is at initial position, there is a prerotation along the Z-axis.When the twisting deformation disappeared, the mechanism's pre-rotation will disappear.Thus, when the moving platform's rotation angle along the Z-axis clockwise is bigger than θ 0 , the mechanism is in singularity positions.

Driven Cables' Tension Disappearance
Considering that cables only can bear tension, when the driven cables' tension disappears, the moving platform will lose the control along with the power from driven cables.Thus, the mechanism is in a singularity position.
When the three driven cables' length are given, the position and quaternion can be calculated according to Section 4.1.Next, the forces at the rods' ends can be calculated from Section 3.1, and the driven cable's force F li can be obtained too.
When the directions of cables' force were reverse to vector A i B i , the mechanism is in singularity positions.

Workspace
In this section, the parasitic motion workspace and the reachability of the novel mechanism are obtained based on forward position analysis and singularity analysis.The key points in analyzing the workspace are the mechanism's kinematical constraints and singularity.The kinematical constraints include action range and physical joint range.
The actuator's range limits the mechanism's workspace.Because of the actuator's physical structure, its range is usually between a minimum value and a maximum value.The corresponding actuator's constraint condition can be expressed as follows: Physical joints also limit the mechanism's range of motion.As mentioned in Section 3.5, the compliant rods are equivalently replaced by a ball hinge joint.However, a ball hinge joint's range is limited by its physical structure.For this mechanism, the moving platform and base platform have the same structure size, so the maximum angle between these two platforms is 90 degrees.Thus, the angle between the moving platform's normal and Z-axis should be less then 90 degrees.According to the rotation matrix, this constraint can be present as According to Section 4.3, the mechanism has two kinds of singularity position.The first depends on the pre-rotation angle.When the prototype is manufactured, the prerotation angle θ 0 is set as 60 degrees.The second singularity position is related to the mechanism's forward kinematic and force analysis in Section 2.2.
Figure 9a shows the mechanism's workspace with singularity positions.The blue plane is the border of the pre-rotation angle, and the green parts are the positions of that mechanism's driven cable without tension.After excluding singularities, the mechanism's final workspace is shown in Figure 9b.

Validation and Results
To verify the simplified model and improve the algorithm's accuracy and sta two group experiments are designed in this section.An angle sensor is fixed on the ing platform to measure its Z-Y-X Euler angles.

Validation and Results
To verify the simplified model and improve the algorithm's accuracy and stability, two group experiments are designed in this section.An angle sensor is fixed on the moving platform to measure its Z-Y-X Euler angles.

Experimental Design
These two group experiments are designed to test the mechanism's forward and inverse kinematics model.In the forward kinematics experiment, a set of input parameters are designed and input to the prototype and simulation.Next, the output Euler angles from the prototype and simulation will be recorded.In the inverse kinematics experiment, a set of moving platform movement trajectories is designed, and the input parameters are calculated according to the simplified model.Next, the input parameters are input to the prototype, and the output angles of the moving platform are recorded.

Experimental Setup
The experimental devices (shown in Figure 10) include a novel mechanism's prototype, an angle sensor and an analysis system.The prototype applies three high-performance actuators (dynamixel mx64t).The actuators whose maximum output torque is 5 Nm can achieve 360 • rotation.The angle sensor (jy901) weighs about 50 g and integrates a nine-axes accelerometer and dynamic Kalman filter, and can output angles, angular velocities and angular accelerations in three axes.The analysis system is a high-performance workstation with Intel i7-6850 CPU, Amd w5100 graphic card.Communication with the prototype and the sensor is carried out using Matlab installed on a workstation.Two serial ports are used in Matlab.One sends input parameters to the prototype with RS485 protocol.The other one receives angle data from the sensor with the RS232 protocol.The moving platform of the prototype weighs about 95 g.In the forward kinematics experiment, three sets of input parameters (shown in Fig ure 11a) are designed.The first set is designed with only one actuator driven along a straight line.The second set is driven by two actuators along the same straight line.In the third set, three actuators move along a straight line, but one of them moves in the opposite direction.The workstation sends these parameters into the prototype and drives the pro totype movement.Next, the Euler angle data from the angle sensor is recorded.
In an inverse kinematics experiment, three sets of moving platform movement tra jectories (shown in Figure 11b) are designed.In these three trajectories, the moving plat form rotates along three axes.Next, the input parameters are calculated with the simpli fied model, and be sent to the prototype.Then, the moving platform's angles are recorded In the forward kinematics experiment, three sets of input parameters (shown in Figure 11a) are designed.The first set is designed with only one actuator driven along a straight line.The second set is driven by two actuators along the same straight line.In the third set, three actuators move along a straight line, but one of them moves in the opposite direction.The workstation sends these parameters into the prototype and drives the prototype movement.Next, the Euler angle data from the angle sensor is recorded.

Conclusions
In this paper, a novel compliant mechanism with three rotational degrees of freedom is proposed.The compliant rods are simulated with Cosserat rod models, and the mechanism is analyzed with a back-propagation neural network and chaos-enhanced accelerated particle swarm optimization method.Next, the mechanism's mobility is verified by numerical analysis, and the mechanism is simplified into a rigid mechanism whose compliant rods are replaced by a spatial ball hinge.With the simplified model, the mechanism's forward and inverse kinematics can be calculated more effectively.In addition, with force analysis, the mechanism's singularity is calculated, and the workspace without singularity positions is obtained.Finally, experiments with a prototype verify the mobility analysis and simplified model.In Figure 12a, the experiment data almost coincides with the simulation data in the Xand Y-axes, showing the same curve trend.However, there are still some differences, as the experimental data's amplitude is lower than the simulation data.Two main reasons explain the difference.First, plastic deformation reduces the rotation angle, and deviations occur mainly in positions with large angles.Second, deviations also occur at inflection points, especially at junctions between two segments.The actuator's acceleration is not continuous at inflection points, causing a shiver in the moving platform and magnifying the sensor's error.In the Z-axis, the experimental data shows the same trend as the simulation data, but its deviation is more obvious than in the X-and Y-axes.This is mainly due to the sensor, which measures the Z-axis angle differently from the X-and Y-axes.The angle in the Z-axis is obtained by measuring the Earth's magnetic induction line, which is easily influenced by external factors, such as the electrical motor's magnetic field.
In Figure 12b, the experiment data in the inverse kinematic experiment almost coincides with the simulation data for the three axes, and there are still errors because of the plastic deformation, inflection points, and external interferences.However, the errors are lower than in the forward kinematic experiment.This means the simplified model performs better in inverse kinematic, and is good for controlling the mechanism.
In summary, the experiment results tally with the simulation results in both forward and inverse kinematics.The experiment results not only verify the simplified model's accuracy, but also prove that the novel mechanism has three rotational degrees of freedom.

Conclusions
In this paper, a novel compliant mechanism with three rotational degrees of freedom is proposed.The compliant rods are simulated with Cosserat rod models, and the mechanism is analyzed with a back-propagation neural network and chaos-enhanced accelerated particle swarm optimization method.Next, the mechanism's mobility is verified by numerical analysis, and the mechanism is simplified into a rigid mechanism whose compliant rods are replaced by a spatial ball hinge.With the simplified model, the mechanism's forward and inverse kinematics can be calculated more effectively.In addition, with force analysis, the mechanism's singularity is calculated, and the workspace without singularity positions is obtained.Finally, experiments with a prototype verify the mobility analysis and simplified model.

Figure 1 .
Figure 1.Three-dimensional model of the novel mechanism.

Figure 1 .
Figure 1.Three-dimensional model of the novel mechanism.

Figure 2 .
Figure 2. The structure diagram of the mechanism (a), moving platform (b), the base platform (d); the force analysis of the moving platform (c) and the base platform (e).

Figure 2 .
Figure 2. The structure diagram of the mechanism (a), moving platform (b), the base platform (d); the force analysis of the moving platform (c) and the base platform (e).

Figure 6 .
Figure 6.Relative rotational center and relative rotation radii.

Figure 6 .
Figure 6.Relative rotational center and relative rotation radii.

Figure 6 .
Figure 6.Relative rotational center and relative rotation radii.

Figure 7 .
Figure 7. Simulation of RRC (a), RRR (b) and the distance between the RRC and origin point (c).

Machines 2023 , 18 Figure 7 .
Figure 7. Simulation of RRC (a), RRR (b) and the distance between the RRC and origin point (c).

Figure 8 .
Figure 8.The overall simplified model (a), separation model (b), and top view (c) of the novel mechanism.

Figure 8 .
Figure 8.The overall simplified model (a), separation model (b), and top view (c) of the novel mechanism.

Figure 9 .
Figure 9. Workspace with the singularity positions (a) and workspace without the singularity positions (b).

Figure 12 .
Figure 12.Experiment and simulation results of the forward kinematics (a) and inverse kinematics (b).

Figure 12 .
Figure 12.Experiment and simulation results of the forward kinematics (a) and inverse kinematics (b).

Table 1 .
Material parameters for compliant rods.
b) calculating the fitness value.(c) renewing the P best , G best and particles' speed and location.(d) judging the termination condition.If the G best meets the termination condition, output the G best ; if not, come back to step (a).