Next Article in Journal
The Study on Single-Event Effects and Hardening Analysis of Frequency Divider Circuits Based on InP HBT Process
Next Article in Special Issue
Micro-Scale Particle Tracking: From Conventional to Data-Driven Methods
Previous Article in Journal
Impact of Ceramic Micropillar Array and Fiber Layer Composite Structure on Kinematic and Heat Transfer Characteristics of Single Droplet Impacting a Wall
Previous Article in Special Issue
A Survey of Recent Developments in Magnetic Microrobots for Micro-/Nano-Manipulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compliance and Kinetostatics of a Novel 2PRS-2PSS Compliant Parallel Micromanipulator: Modeling and Analysis

Hubei Key Laboratory of Modern Manufacturing Quantity Engineering, School of Mechanical Engineering, Hubei University of Technology, Wuhan 430068, China
*
Author to whom correspondence should be addressed.
Micromachines 2024, 15(4), 526; https://doi.org/10.3390/mi15040526
Submission received: 14 March 2024 / Revised: 8 April 2024 / Accepted: 12 April 2024 / Published: 14 April 2024
(This article belongs to the Special Issue Advanced Micro-/Nano-Manipulation and Positioning Techniques)

Abstract

:
A novel 2PRS-2PSS (P represents the prismatic pair, R represents the revolute hinge, S represents the spherical hinge) compliant parallel micromanipulator with two translational DOFs and two rotational DOFs is presented, and its compliance model and kinetostatic model are sequentially developed and analyzed. Initially, an analytical model used to describe the compliance of this micromanipulator was developed using the compliance matrix method (CMM). Through a comparison with finite element analysis, the accuracy of this analytical model is confirmed, and the influence of various dimensional and structural parameters on the compliance behavior is investigated. Subsequently, the micromanipulator is treated as an equivalent spring system, allowing for the derivation of its governing equation based on the established compliance model. From this equation, a kinetostatic model relating input forces to output displacements is derived. Validation of this model is performed by comparing analytical results with finite element simulations under specific motion trajectories, revealing a maximum relative error of 6.18%. This close agreement verifies the accuracy of the kinetostatic model. Finally, the impact of the parameters of the flexure hinge on the mapping matrix is examined to offer insights into minimizing undesired displacements, providing valuable guidance for optimizing the micromanipulator’s performance.

1. Introduction

In recent years, fields such as biology, medicine, and industrial assembly have been moving toward smaller operating objects and higher precision requirements. Examples include cell injection, microsurgery, and fiber optic docking. There is an increasing demand for micromanipulators with excellent performance [1,2]. Traditional rigid manipulators, due to issues like clearance, friction, and backlash in their mechanisms, have made it difficult to meet the high-precision requirements of micromanipulation systems. In contrast, micromanipulators based on compliant parallel mechanisms utilize the deformation of flexure hinges to transmit motion, effectively addressing these issues. They also offer advantages such as high stiffness and compact structure, enabling high-precision manipulations [3,4,5,6,7,8]. Therefore, it has attracted a large number of scholars to conduct research on micromanipulators based on compliant parallel mechanisms with various configurations.
Compliance plays a crucial role in the design and analysis of compliant parallel mechanisms, as it serves as a key performance metric. The study of the kinematics and dynamics of compliance mechanisms relies on the foundation provided by the compliance model [9,10,11]. Zhang et al. [12] conducted research on an integrated compliant redundant parallel mechanism for the XY micro-adjustment of photolithography projection lens optical components. They employed the compliance matrix method to calculate the compliance of individual limbs within the mechanism, including output compliance, input compliance, and the input–output ratio of the system. Li et al. [13] proposed a novel piezoelectric-driven XY micro-motion platform. They used the compliance matrix method to establish the compliant model of the mechanism and conducted static and dynamic analyses. Mishra et al. [9] examined a 6-DOF serial-parallel compliant micromanipulator built on a Stewart platform featuring a displacement amplification structure. By applying the compliance matrix method, they developed a compliance model for the mechanism and investigated the impact of various parameters on its compliance. Ren et al. [14] proposed the design of a novel n-4R-compliant parallel micro-pointing mechanism. They utilized the compliance matrix method to develop a compliance model for the mechanism and investigated how different parameters affected the compliance model.
In addition to compliance, kinetostatics is another indispensable area of research for compliant parallel mechanisms. The kinetostatic model of these mechanisms serves as the foundation for motion control and significantly aids in studying the reduction of parasitic displacement within the mechanism. It is widely acknowledged that compliant parallel mechanisms offer advantages such as high precision, the absence of friction, the absence of backlash, and a compact structure due to the utilization of flexure hinges. However, analyzing the kinetostatics in a compliant parallel mechanism presents greater challenges compared to a rigid one, primarily due to the inherent relationship between kinematics and the elasticity in flexure hinges. This coupling necessitates a comprehensive modeling approach that extends beyond traditional kinematics or statics analysis [15]. In the past few decades, researchers have devised numerous techniques and methodologies to model the kinetostatics model of compliant mechanisms, such as the pseudo-rigid-body model, Castigliano’s second theorem, compliance matrix method, elastic beam theory, etc. [16]. Midha et al. [17] analyzed the fixed-oriented compliant beam with an inflection point by utilizing the concept of a pseudo-rigid-body model and improved the calculation accuracy by setting the inflection point. Ni et al. [18] established the kinematic, static, and dynamic models of the compliant piezoelectric micro-cracker based on the three-stage amplification by using the pseudo-rigid-body model method. Wadikhaye et al. [19] analyzed the input and output stiffness of the serial motion XYZ scanner by utilizing Castigliano’s second theorem and estimated the first natural frequency and the travel range of the scanner using the input–output stiffness. Dong et al. [20] developed a bridge amplifier and established a motion-statics model based on the compliance matrix method and a dynamics model based on the Lagrange method for the mechanism. Wang et al. [21] developed a kinetostatic model utilizing the Cosserat rod theory. The proposed model aimed to minimize positioning errors resulting from external loads. Ren et al. conducted a series of studies on compliant parallel mechanisms. In one study [14], they designed an n-4R-compliant micro-pointing mechanism and developed its kinetostatic model using the compliance matrix method. They investigated the impact of parameters on the kinetostatic mapping matrix. In a subsequent study [22], Ren et al. designed a novel 3PSS/S-compliant micro-turntable. They established the kinetostatic model of the mechanism using both the compliance matrix method and the pseudo-rigid-body model method. The two models were then compared in terms of accuracy. More recently, Ren et al. [23] introduced a generalized 3-PSS-compliant parallel micro-motion platform. They focused on investigating the compliance and kinetostatic models of this platform. Ling et al. [24] proposed a semi-analytical matrix displacement method for modeling the kinetostatics of compliant mechanisms. This method is particularly useful for complex compliant mechanisms that comprise serial-parallel substructures. Based on the work of Ling et al., Arredondo-Soto et al. [25,26] combined various methods from other literature to derive a more general approach for kinetostatic modeling. They used this approach to establish the relationship between input force/displacement and the displacement of the mobile platform for a 3-RRR spherical-compliant parallel mechanism.
Currently, there is limited research on 4-DOF-compliant parallel mechanisms. The 4-DOF parallel mechanism offers broader application possibilities compared to mechanisms with fewer degrees of freedom (less than 4-DOF). Simultaneously, the 4-DOF parallel mechanism features a simpler structure and lower control complexity compared to its 6-DOF counterpart. Therefore, this paper presents a novel 4-DOF 2PRS-2PSS compliant parallel micromanipulator. The compliance model and kinetostatic model of the micromanipulators are sequentially developed and analyzed. Firstly, the compliance matrix of a single flexure hinge is computed, followed by deriving the compliance of the PSS branch and PRS branch separately. Subsequently, the overall compliance model of the mechanism is constructed using the compliance matrix method, and its accuracy is verified through finite element simulation. The study then investigates the effects of variations in both the dimensional parameters of the mechanism and the structural parameters of the flexure hinges on the overall compliance. Secondly, the mechanism is simplified as an equivalent spring system. The governing equation of this equivalent spring system is derived from the previously established compliant model. By utilizing this governing equation, the kinetostatic model of the mechanism is formulated. The accuracy of the kinetostatic model is validated by comparing analytical calculations with finite element simulations of the specified motion trajectory of the mechanism. Finally, the impact of the structural parameters of the flexure hinges on the mapping matrix of the kinetostatic model is analyzed.

2. Structure of 2PRS-2PSS Compliant Parallel Micromanipulator

The 2PRS-2PSS compliant parallel micromanipulator is a 4-DOF micromanipulation platform with two translational and two rotational degrees of freedom. As shown in Figure 1a, it consists of a mobile platform, a fixed platform, two symmetrically arranged PRS branches, and two symmetrically arranged PSS branches. The links of each branch are connected to the mobile platform through flexure spherical hinges. The links of each PRS branch are connected to the prismatic pair through flexure revolute hinges, while the links of each PSS branch are connected to the prismatic pair through flexure spherical hinges. To enhance the compliance and adaptability of the mechanism, an equivalent compliant prismatic pair consisting of eight flexure revolute hinges is employed. The flexure revolute hinges that connect to the link have the same structure as the eight flexure revolute hinges, but their structural parameters are different.
The structural parameters and coordinate system of the mechanism are set as shown in Figure 1a. The length of the links is denoted as l. The radius of the circle formed by the rotation centers of the four hinges connected to the equivalent compliant prismatic pair is defined as the fixed platform radius, Ra, with the center of the circle denoted as O’. The radius of the circle formed by the rotation centers of the four hinges connected to the mobile platform is defined as the mobile platform radius, Rb, with the center of the circle denoted as O. The assumption is made that the radius of a fixed platform is always larger than that of the mobile platform. The mechanism utilizes four piezoelectric ceramics for actuation. By controlling the inputs of these actuators, the mobile platform can generate movement along the x and z axes and rotation around the x and y axes. Therefore, the micro-manipulator can achieve position adjustments in the x and z directions during operation, while also having the capability to adjust the orientation angles around the x and y axes simultaneously.

3. Modeling and Analysis of the Compliance of the 2PRS-2PSS Compliant Parallel Micromanipulator

The organization of this section is as follows. Firstly, the compliance model of the micromanipulator is developed using the compliance matrix method. Secondly, the compliance model is validated through finite element simulation to ensure its accuracy. Finally, the impact of some structural parameters on the overall compliance of the micromanipulator is analyzed.

3.1. Compliance Matrix Method

The compliance matrix method considers the flexure hinge as a compliant element while treating other parts as rigid bodies. During the process of compliance modeling, the flexure hinge is regarded as a multi-dimensional hinge, and the overall compliance model of the mechanism is obtained through coordinate transformations. Therefore, the compliance matrix method offers higher accuracy and wider applicability.
Under the condition that one end of the flexure hinge is fixed, deformation will be generated by exerting the forces and moments at the free end of the flexure hinge. Under the assumption of small elastic deformations, the principle of linear superposition holds, meaning that the deformation of the flexure hinge in each direction can be calculated by summing the individual deformations caused by each applied force or moment. Assuming the forces and moments are represented as F = [fx, fy, fz, mx, my, mz]T, and the corresponding linear displacements and angular displacements are represented as X = [δx, δy, δz, θx, θy, θz]T, the compliance matrix of a single flexure hinge is defined as C. The relationship between the input forces and displacements at the end of the hinge is defined as follows [27]:
X = C F

3.2. Compliance Matrix of Single Flexure Hinge or Compliant Element

Before calculating the overall compliance of the micromanipulator, the compliance of a single flexure hinge used in the mechanism should be calculated first. The 2PRS-2PSS compliant parallel micromanipulator employs flexure spherical hinges and flexure revolute hinges, so the compliance of these two types of hinges needs to be calculated separately. Additionally, the equivalent prismatic pair is treated as a compliant element, and its compliance is calculated through the compliance matrix method.
Figure 2a illustrates the structural parameters and coordinate frame configuration of the flexure revolute hinge. In this figure, rR, w, and tR represent the cutting radius, width, and minimum thickness of the flexure hinge, respectively. Based on its structural characteristics, it can be observed that this hinge exhibits significant compliance in the rotation around the x-axis and can be used as a single-degree-of-freedom revolute hinge. The compliance matrix of the flexure revolute hinge can be expressed as Equation (2) [25], and the formulas for calculating the compliance matrix CR can be found in Appendix A.
The structure parameters and coordinate frame setting of the flexure spherical hinge are shown in Figure 2b, where rS and tS represent the cutting radius and minimum thickness of the flexure hinge, respectively. It has three degrees of freedom for rotation around the x, y, and z axes and can be used as a universal joint. The compliance matrix of this hinge is denoted as CS and its representation method is consistent with Equation (2), and the formulas for calculating the compliance matrix CS can be found in Appendix A.
C R = C δ x , f x 0 0 0 C δ x , m y 0 0 C δ y , f y 0 C δ y , m x 0 0 0 0 C δ z , f z 0 0 0 0 C θ x , f y 0 C θ x , m x 0 0 C θ y , f x 0 0 0 C θ y , m y 0 0 0 0 0 0 C θ z , m z
The structure of the equivalent compliant prismatic pair is shown in Figure 3. The prismatic pair consists of four identical branches, each of which includes two flexure revolute hinges. These branches are named P-branch i. The flexure revolute hinges used here have the same structure as described in Figure 2a but with different dimensional parameters. The cutting radius, width, and minimum thickness of the flexure revolute hinge used here are expressed as rP, wP, and tP, respectively. The structural parameters a, b and c in the equivalent compliant prismatic pair represent the distance of the symmetrical branch, the distance of the adjacent branch and the length of the connecting rod connecting the two flexure revolute hinges on a single branch, respectively. To facilitate the subsequent modeling of overall compliance, the compliance of the equivalent compliant prismatic pair is modeled at coordinate system P-xyz to form local compliant elements. Thus, the compliance of P-branch 1 and 2 can be calculated as follows:
C P b m P = n = 1 2 T P m n P C R P m n P T P m n P T   ,   m = 1 , 2
T P m n P = R P m n P R P m n P P P m n P 0 R P m n P
where T P m n P is an adjoint transformation matrix that represents the transformation of the compliance matrix from the local coordinate system, Pmn-xyz, to the reference coordinate system, P-xyz. m represents the branch number, and n represents the position of the hinge on the branch. Defining the rotation transformation matrix, R P m n P , and the position transformation matrix, P P m n P , their specific expressions are as follows:
R P m n P = c γ c β c γ s β s α s γ c α c γ s β c α + s γ s α s γ c β s γ s β s α + c γ c α s γ s β c α c γ s α s γ c γ s α c γ c α ,   P P m n P = 0 h 3 h 2 h 3 0 h 1 h 2 h 1 0
where α, β, and γ represent the angle of rotation around the x, y, and z axes, respectively. s and c represent sin and cos, respectively. The elements h1, h2, and h3 in P P m n P represent the positions of local coordinates in the reference coordinate system; their coordinate form is h = (h1, h2, h3). The parameters of the adjoint matrix of the coordinate transformation in the prismatic pair are constructed as shown in Table 1.
According to the principles of stiffness and compliance superposition in series-parallel hybrid mechanisms [27], it is known that the overall compliance of the equivalent prismatic pair at the coordinate P-xyz is as follows:
C P = m = 1 4 C P B m P 1 1
The structure and structural parameters of the local branch are similar, so the compliance of P-branch 3 and P-branch 4 in the reference coordinate system P-xyz can be obtained by rotating the compliance of P-branch 1 and P-branch 2. The compliance matrix is as follows:
C P B 3 P = T 1 3 C P B 1 P T 1 3 T C P B 4 P = T 2 4 C P B 2 P T 2 4 T ,   T 1 3 = T 2 4 = R z , π R z , π

3.3. Compliance Modeling of Single Branch

In general, it is meaningful to discuss the compliance of a compliant mechanism only when it is uniformly discussed in the same coordinate system [27]. Since the 2PRS-2PSS compliant parallel micromanipulator consists of PRS and PSS branches, the compliance of a single branch should be calculated first. The structures of branch 1 and branch 2 are shown in Figure 4. The compliance of the flexure hinge needs to be unified into the reference coordinate system first. Taking branch 1 (PSS branch) and branch 2 (PRS branch) as an example, the conversion Equation (8) describes the process of unifying the compliance of the flexure hinge on the branch to the reference coordinate system O-xyz.
C P 11 O = T O 11 O C P 11 O 11 T O 11 O T C S 12 O = T O 12 O C S 12 O 12 T O 12 O T C S 13 O = T O 13 O C P 13 O 13 T O 13 O T ,   C P 21 O = T O 21 O C P 21 O 21 T O 21 O T C R 22 O = T O 22 O C R 22 O 22 T O 22 O T C S 23 O = T O 23 O C S 23 O 23 T O 23 O T  
where T O i j O is the transformation matrix of 6 × 6, which shows that the compliance matrix is transformed from the local coordinate system, Oij-xijyijzij, to the global coordinate system, O-xyz. The rotation transformation matrix is denoted as R O i j O , and the translation transformation matrix is denoted as P O i j O . Subscript i represents the number of branches in the mechanism, and j represents the position index of the flexure hinge or compliance element on the branch, following the numbering rules shown in Figure 4.
When the flexure hinges or compliant elements are connected in series, the compliance is superimposed [27]. Therefore, the compliance of branch 1 and branch 2 in the reference coordinate system can be obtained as follows:
C B 1 O = C P 11 O + C S 12 O + C S 13 O C B 2 O = C P 21 O + C S 22 O + C S 23 O

3.4. Compliance Model of the Micromanipulator

According to the principles of stiffness and compliance superposition in series-parallel hybrid mechanisms [27], the compliance of the 2PRS-2PSS compliant parallel micromanipulator can be determined as follows:
C 2 PRS - 2 PSS = i = 1 4 C B i O 1 1
where C2PRS-2PSS represents the total compliance of the 2PRS-2PSS compliant parallel micromanipulator, and C B i O denotes the compliance of each individual branch within the reference coordinate system.
Due to the symmetry of the mechanism, the compliance matrix of branch 3 and branch 4 can be easily obtained by performing coordinate transformations on the compliance matrix of branch 1 and branch 2, as shown in Equation (11).
C P B 3 P = T 1 3 C P B 1 P T 1 3 T C P B 4 P = T 2 4 C P B 2 P T 2 4 T ,   T 1 3 = T 2 4 = R z , π R z , π
where T 1 3 and T 2 4 are the adjoint transformation matrix and Rz,π is the rotation transformation matrix, indicating that the compliance matrix rotates 180° around the z-axis of the reference coordinate system, O-xyz. The related parameters of the adjoint matrix of the overall compliance coordinate transformation of the mechanism are listed in Table 2. In this table, θ represents the angle between the links and the z-axis of the reference coordinate system, and its value can be expressed as follows:
θ = arcsin R a R b l + 2 r S

3.5. Verification of Effectiveness of the Compliance Model

To validate the accuracy of the compliance model, finite element analysis is conducted in this section. Table 3 outlines the dimensional parameters of the mechanism, as well as the structural parameters of the flexure hinges utilized. The total compliance can be expressed by Equation (13). Substituting the parameters (listed in Table 3) into Equation (10), one can obtain the analytical results of the compliance of the 2PRS-2PSS compliant parallel micromanipulator.
C 2 PRS - 2 PSS = C δ x , f x 0 0 0 C δ x , m y 0 0 C δ y , f y 0 C δ y , m x 0 0 0 0 C δ z , f z 0 0 0 0 C θ x , f y 0 C θ x , m x 0 0 C θ y , f x 0 0 0 C θ y , m y 0 0 0 0 0 0 C θ z , m z
The compliance model of the 2PRS-2PSS compliant parallel micromanipulator underwent validation using the commercial software ANSYS 19.2, affirming the precision of the analytical findings. To optimize computational efficiency and accuracy, the flexure hinge was discretized into tetrahedral meshes with a size of 0.5 mm, while the remaining components were segmented into tetrahedral meshes with a size of 3 mm. The material parameters pertaining to the flexure hinge are listed in Table 4.
By comparing the compliance elements on the main diagonal between the analytical results and the finite element results, the relative errors are calculated and summarized in Table 5. The results show that the relative errors for the six main compliances of the compliance matrix are all below 7.2%. This high level of consistency indicates the accuracy of the established theoretical compliance model. The discrepancies between the analytical and finite element results can be attributed to the following factors: (1) errors in the compliance model of individual flexure hinges; and (2) theoretical modeling is carried out based on the assumption that all components except for the flexure hinges are rigid. However, in the finite element simulation, these ‘rigid’ components simulated their rigid behavior by being set with larger stiffness values. Thus, even if their stiffness is set to be sufficiently large, they may still undergo slight deformations.

3.6. Analysis of Compliance Performance of 2PRS-2PSS compliant Parallel Micromanipulator

When designing compliant parallel mechanisms, it is necessary to consider various performance indicators such as compliance, stiffness, and accuracy, etc. And different application scenarios have different requirements for these performances. For applications that require strong adaptability and sensitivity to external disturbances, greater compliance is beneficial. However, for applications that demand high precision and stiffness, it is necessary to control compliance to ensure the stability and accuracy of the mechanism. Therefore, analyzing the influence of parameters in the mechanism on compliance is essential [28].
In this section, the compliance model established in Section 3.4 is used to analyze the compliance of the micromanipulator. Here, Cδ and Cθ are defined as translational compliance and rotational compliance, respectively. The parameters related to compliance and their variation ranges are listed in Table 6, and parameters rS and rR are always the same, so we unified their symbols as r. The variation in overall compliance with these parameters is shown in Figure 5.
From Figure 5, we can draw the following conclusions: (1) in translational compliance Cδ, Cδx,fx is always the largest and most affected by parameter changes, while Cδz,fz is always the smallest. At the same time, in rotational compliance Cθ, Cθx,mx is always smaller than Cθy,my. (2) The structural parameters of the flexure hinges have a more significant influence on Cδ and Cθ compared to the dimensional parameters of the mechanism. Consequently, when designing the mechanism, it is suggested that the structural parameters of the flexure hinges be adjusted with priority to achieve the desired compliance in the micromanipulator.
Apart from the influence of these mentioned parameters, it is sometimes necessary to understand the way the micromanipulator compliance scales with the defining geometric diameters. Given the utilization of flexure hinges in the mechanism, there are two possible scaling scenarios: scaling only the dimensional parameters of the mechanism (represented by the scaling coefficient e1) or simultaneously scaling both the dimensional parameters of the mechanism and the flexure hinge parameters (represented by the scaling coefficient e2). The variations of the mechanism’s compliance with the scaling coefficients e1 and e2 are shown in Figure 6 and Figure 7, respectively. It can be observed from Figure 6 that Cδ and Cθ are directly and inversely correlated with the scaling coefficient e1, respectively. Comparatively, Cδx,fx is greatly affected, Cδy,fy is less affected, while Cδz,fz is almost unaffected, as shown in Figure 6a. It is also noted from Figure 6b that the overall decrease of Cθx,mx, Cθy,my, and Cθz,mz is generally gradual. It can be seen from Figure 7 that both Cδ and Cθ are inversely correlated with the scaling coefficient e2. And Figure 7a exhibits the same feature as Figure 6a in that Cδx,fx is greatly affected, Cδy,fy is less affected, while Cδz,fz is almost unaffected. Compared to Figure 6b, Cθ in Figure 7b decreases more sharply in the range with smaller proportion coefficients.

4. Kinetostatic Modeling of 2PRS-2PSS compliant Parallel Micromanipulator

The force–displacement relationship at the center point of the end-effector platform can be easily obtained by the compliance model. However, the input force and output displacement of 2PRS-2PSS compliant parallel micromanipulation exist in different coordinate systems. Therefore, it is of vital necessity to establish the kinetostatic model of this micromanipulation.

4.1. Kinetostatic Modeling under Single Input Force

To achieve the output displacements of the mobile platform, four input forces are exerted at the center point of the bottom surface of the equivalent prismatic pairs, as shown in Figure 8a. The local coordinate systems are established at this center point, allowing us to obtain the local generalized input forces as follows: Fi = [fi,x, fi,y, fi,z, mi,x, mi,y, mi,z], i = 1, 2, 3, 4. The value of i denotes system Fi-xFiyFizFi. The displacement of the center point of the mobile platform relative to the global coordinate system O-xyz is expressed as follows: U2PRS-2PSS = [δx, δy, δz, θx, θy, θz]. Under the assumption of linear deformations, the force–displacement mapping relationship of the mechanism can be determined using the principle of superposition.
To facilitate the kinetostatic modeling, we simplify the 2PRS-2PSS compliant parallel micromanipulator as an equivalent spring system, as presented in Figure 8b,c. Note that the concept of equivalent stiffness can also be found in [27,29]. Since the same branches of the mechanism are symmetrically arranged, only branch 1 and branch 2 are taken for analysis. The input forces applied to these two branches are denoted as F1 and F2, respectively. Correspondingly, the resulting displacements are U1 and U2, respectively. Then, the elastic deformation of the system can be described by the governing equation (Equation (14)) based on Hooke’s law. The stiffness matrices in Equation (14) are computed by Equation (15). It should be noted that the stiffness matrix in Equation (15) can be calculated by using Equations (8) and (10), as shown in Equations (18) and (19).
K O O F i K O F i K F i O K F i F i U i U F i = F O F i
K O O F i = K B i B O + C 2 PRS - 2 PSS 1 K B i K F i F i = K B i A F i + K B i B F i K O F i = T O O T K B i B O T O F i 1 K F i O = T O F i T K B i B O T O O 1
where
T O O = I 6 × 6 ,   T O F i = R O F i P O F i T O F i 0 3 × 3 R O F i
where
P O F i = 0 z i y i z i 0 x i y i x i 0 ,   R O F i = I 3 × 3 ,   i = 1 , 2
K B 1 A O = C P 11 O 1 ,   K B 1 B O = C S 12 O + C S 13 O 1 K B 2 A O = C P 21 O 1 ,   K B 2 B O = C R 22 O + C S 23 O 1
K B i A F i = T O F i K B i A O 1 T O F i T 1 K B i B F i = T O F i K B i B O 1 T O F i T 1 ,   i = 1 , 2
In Equation (15), the superscripts O and Fi associated with each stiffness matrix denote that the stiffness matrix is defined with respect to the coordinate frames O-xyz and Fi-xFiyFizFi, respectively. In Equations (16) and (17), the symbol I represents the unit matrix, and xi, yi, and zi denote the position coordinates of the local coordinate system Fi-xFiyFizFi in the reference coordinate system.
Since there is no force applied to the mobile platform, FO in Equation (14) can be set to 0, which yields the following:
U i = C T O F i F i ,   i = 1 , 2
where
C T O F i = K O O F i K O F i K F i F i 1 K F i O 1 K O F i K F i F i 1 ,   i = 1 , 2
Then, the kinetostatic model of this micromanipulator under a single input force is obtained by Equation (20).

4.2. Kinetostatic Modeling of 2PRS-2PSS compliant Parallel Micromanipulator

Define the displacement of the center of the mobile platform as U2PRS-2PSS. Based on the principle of superposition, displacement of U2PRS-2PSS resulting from the combined exertions of forces F1, F2, F3, and F4 can be considered as the superposition of the displacements U1, U2, U3, and U4 generated by F1, F2, F3, and F4 acting alone. Therefore, the displacement of U2PRS-2PSS can be expressed as follows:
U 2 PRS - PSS = i = 1 4 U i
According to Equations (20) and (22), one can then obtain the kinetostatic model of a 2PRS-2PSS compliant parallel micromanipulator, as follows:
U 2 PRS - PSS = i = 1 4 C T O F i F i
where the matrix, C T O F i , represents the mapping relationship between the force, Fi, and the displacement, Ui (i = 1, 2, 3, 4).
Since branches 3 and 4 are symmetrically arranged with branches 1 and 2, respectively, C T O F 3 and C T O F 4 can be easily obtained through the rotational transformation of C T O F 1 and C T O F 2 , respectively.
C T O F i + 2 = T π T C T O F i ,   T π = R z , π 0 3 × 3 0 3 × 3 T z , π
where the matrix C T O F i + 2 (i = 1, 2) describes the relationship between the input forces Fi+2 (i = 1, 2) and the corresponding output displacements Ui+2 (i = 1, 2). And Rz,π represents the coordinate transformation matrix with 180° rotation around the z-axis of the coordinate system O-xyz.

5. Verification and Analysis of the Kinetostatic Model

In Section 4, the kinetostatic model of the 2PRS-2PSS compliant parallel micromanipulator has been established. Here, we validate the accuracy of the established model through a comparative analysis between theoretical calculations and finite element simulations, utilizing an example. Additionally, we investigate the impact of variations in the structural parameters of the flexure hinges on the kinetostatic model. It is important to note that the parameters in this example are consistent with those presented earlier in Section 3.5.

5.1. Numerical Calculation and Simulation Analysis of the Kinetostatic Model

The 2PRS-2PSS compliant parallel micro-manipulator is a 4-DOF mechanism. Therefore, the elements related to the movements of four functional directions of the mechanism in the model are extracted from Equation (23), and the simplified kinetostatic model is obtained as in Equation (25).
δ x δ z θ x θ y = C δ x , f 1 , z C δ x , f 2 , z C δ x , f 3 , z C δ x , f 4 , z C δ z , f 1 , z C δ z , f 2 , z C δ z , f 3 , z C δ z , f 4 , z C θ x , f 1 , z C θ x , f 2 , z C θ x , f 3 , z C θ x , f 4 , z C θ y , f 1 , z C θ y , f 2 , z C θ y , f 3 , z C θ y , f 4 , z f 1 , z f 2 , z f 3 , z f 4 , z
The mapping matrix in Equation (25) is composed of the elements of the 1st, 3rd, 4th, and 5th rows, and the 3rd column of the matrix C T O F i , and it is defined as C T O F r c . And f1,z, f2,z, f3,z, and f4,z, respectively, denote the components of forces F1, F2, F3, and F4 along the z-axis of the local coordinate system. Substituting the parameters in Table 3 into C T O F r c yields the following:
C T O F r c A n = 0 1.61 × 10 6 0 1.61 × 10 6 8.87 × 10 6 8.93 × 10 6 8.87 × 10 6 8.93 × 10 6 6.73 × 10 5 0 6.73 × 10 5 0 0 5.36 × 10 5 0 5.36 × 10 5
It can be observed from the matrix that the elements of the first row and the fourth row are proportional, which means that the output displacement δx and θy are always proportional in value (but their units are different). To facilitate the calculation of input forces from given output displacements, we extract the first three rows of the matrix, C T O F r c A n , to form a new matrix, C T O F r c ( r 1 , 2 , 3 ) A n . Obviously, Matrix C T O F r c ( r 1 , 2 , 3 ) A n has full column rank, so the input forces can be determined by taking the right inverse of C T O F r c ( r 1 , 2 , 3 ) A n as shown in Equation (27). The trajectory equation is given by Equations (28) and (29), and the motion trajectory is shown in Figure 9. The obtained input forces are shown in Figure 10.
f 1 , z f 2 , z f 3 , z f 4 , z = C T O F r c r 1 , 2 , 3 A n T C T O F r c r 1 , 2 , 3 A n C T O F r c r 1 , 2 , 3 A n T 1 δ x δ z θ x
δ x = L cos 10000 L π δ z = L sin 10000 L π ,   0 L 5 × 10 5
θ x = 2 e t / 3 sin 3 t × 10 4 ,   0 t 2 π θ y = C θ y , f 2 , z / C δ x , f 2 , z δ x
Substituting the obtained input forces into the finite element software for simulation, the corresponding trajectory can be obtained, as shown in Figure 11a. Figure 11b,c presents a comparison between the analytical trajectories and the finite element trajectories. Good consistency of the trajectories demonstrates the accuracy of the established kinetostatic model. The absolute error between the theoretical analysis and the finite element analysis for the motion trajectory is illustrated in Figure 12a,b. It is shown that the absolute errors of the moving trajectory and rotation trajectory are directly correlated with the radius of the spiral trajectory and the rotating angle, respectively.
It can be observed from Figure 12c that the maximum relative error occurs in the x-axis movement direction, stabilizing between 6.176% and 6.180%. The minimum error occurs in the z-axis movement direction, stabilizing between 3.57450% and 3.57454%. The relative error for rotation around the x-axis ranges from 3.6582% to 3.6594%, while for rotation around the y-axis, it ranges from 4.23708% to 4.23730%. Good consistency indicates the effectiveness of the established kinetostatic model.

5.2. Analysis of the Influence of Parameters on the Kinetostatic Model

In Section 4, the kinetostatic model of the 2PRS-2PSS compliant parallel micromanipulator was established, and its correctness was verified in Section 5.1. Since the micromanipulator generates parasitic displacement during its motion, it affects the positioning accuracy of the micromanipulator. Therefore, in this section, an analysis of the influence of the structural parameters of the flexure hinges on the mapping matrix is further conducted, which provides a reference on how to reduce parasitic displacement.
Taking into account the movements in non-functional directions (δy and θz), the relationship between the input forces and output displacements of the mechanism can be expressed as Equation (30), and the mapping matrix in the equation is defined as CA. Substitute the parameters in Table 3 into the mapping matrix CA to obtain the results as shown in Equation (31). The element in matrix CA is defined as the compliance coefficient. It can be seen from the mapping matrix that the mechanism will produce parasitic displacement moving along the y-axis direction. For convenience, the first three and last three rows of the mapping matrix are defined as the translation-related mapping matrix CT and rotation-related mapping matrix CR, respectively. Meanwhile, since the structural symmetry of the mechanism, some of the compliance coefficients in the mapping matrix are the same in absolute values (with the positive or negative sign only indicating the direction of motion). Therefore, we only need to analyze the compliance coefficients C δ x , f 4 , z , C δ y , f 3 , z , C δ z , f 3 , z , C δ z , f 4 , z , C θ x , f 1 , z and C θ y , f 3 , z , where C δ y , f 3 , z is the compliance coefficient related to parasitic displacement. Figure 13 shows the variations of the compliance coefficient CT in terms of the parameters r, w, tR, and tS. Variations of the compliance coefficient CR, in terms of these parameters, are shown in Figure 14. The range of variation for these parameters is listed in Table 6.
δ x δ y δ z θ x θ y θ z = C δ x , f 1 , z C δ x , f 2 , z C δ x , f 3 , z C δ x , f 4 , z C δ y , f 1 , z C δ y , f 2 , z C δ y , f 3 , z C δ y , f 4 , z C δ z , f 1 , z C δ z , f 2 , z C δ z , f 3 , z C δ z , f 4 , z C θ x , f 1 , z C θ x , f 2 , z C θ x , f 3 , z C θ x , f 4 , z C θ y , f 1 , z C θ y , f 2 , z C θ y , f 3 , z C θ y , f 4 , z C θ z , f 1 , z C θ z , f 2 , z C θ z , f 3 , z C θ z , f 4 , z f 1 , z f 2 , z f 3 , z f 4 , z
C A An = 0 1.61 × 10 6 0 1.61 × 10 6 1.34 × 10 7 0 1.34 × 10 7 0 8.87 × 10 7 8.93 × 10 7 8.87 × 10 7 8.93 × 10 7 6.73 × 10 5 0 6.73 × 10 5 0 0 5.36 × 10 5 0 5.36 × 10 5 0 0 0 0
From Figure 13, we can draw the following conclusions: (1) The parasitic displacement-related compliance coefficient C δ y , f 3 , z is always much smaller than the mobile functional-direction-related compliance coefficients ( C δ x , f 4 , z , C δ z , f 3 , z and C δ z , f 4 , z ), indicating that the mechanism has a higher motion accuracy. (2) The parasitic displacement-related compliance coefficient C δ y , f 3 , z is directly correlated to the parameter tS, and inversely correlated to the parameter w, while it is not sensitive to the variations in parameters r and tR. Therefore, reducing parameter tS and increasing parameter w are preferable measures to reduce parasitic displacement. However, such changes to these parameters will also decrease the movement in the functional direction, so designers need to make reasonable trade-offs based on actual requirements. It is worth noting that changing the parameters to reduce parasitic displacement will also affect the rotation functional-direction-related compliance coefficients ( C θ x , f 1 , z and C θ y , f 3 , z ), as shown in Figure 14. From Figure 14, we can draw the following conclusions: (1) The rotation functional-direction-related compliance coefficients, C θ x , f 1 , z , are inversely correlated to the parameter tS, while not significantly affected by the changes of parameters r, w, and tR. (2) Another rotation functional-direction-related compliance coefficient C θ y , f 3 , z is directly correlated to the parameters w and tR, and inversely correlated to the parameter tS, while not sensitive to the variation of parameter r.

6. Conclusions

A novel 2PRS-2PSS compliant parallel micromanipulator is presented, and its compliance model and kinetostatic model are successively established and analyzed. The correctness of the two models is confirmed through FE simulation. The conclusions are as follows:
(1)
In the verification of the compliance model, comparison results show that the maximum relative errors of the elements on the main diagonal of the compliance model between the theoretical calculation and the finite element analysis do not exceed 7.2%, indicating the correctness of the compliance model. Analysis of the effect of the parameters on the compliance model indicates that the structural parameters of the flexure hinges have a more significant influence on the compliance of the micromanipulator compared to the dimensional parameters of the mechanism.
(2)
By equating the 2PRS-2PSS micromanipulator to a spring system, the kinetostatic model of the micromanipulator is established based on the previously established compliance model according to Hooke’s law. In the validation of the kinetostatic model, comparison results show that the maximum relative errors between the theoretical analysis and the finite element analysis are within 6.18%, confirming the accuracy of the kinetostatic model.
(3)
In the analysis of the effect of the parameters on the kinetostatic model, it is found that compared to other parameters, the width (w) of the flexure revolute hinge and the cutting thickness (tS) of the flexure spherical hinge has the most significant influence on the parasitic displacement. Adjusting these parameters by increasing w and reducing tS can effectively minimize the parasitic displacement of the mechanism, thereby improving motion accuracy.
It is worth noting that when adjusting parameters to reduce the parasitic displacement of the micromanipulator, the overall compliance of the micromanipulator also changes. On the one hand, increasing the width (w) of flexure revolute hinges is advantageous for reducing parasitic displacement, but it also decreases overall compliance. On the other hand, reducing the minimum cutting thickness (tS) of flexure spherical hinges can also reduce parasitic displacement but increase overall compliance. All parameter adjustment schemes aimed at reducing parasitic displacement will to some extent decrease the working range of the micromanipulator because the impact of parameters on the compliance coefficient is universal and consistent. For micromanipulators, reducing parasitic displacement is beneficial for improving accuracy, and when external disturbances (i.e., external forces acting on the mobile platform) are significant, reducing overall compliance also helps improve accuracy. When external disturbances during micromanipulator operation are significant, lower overall compliance (i.e., higher overall stiffness) is beneficial. As long as the mechanism’s working range is within the allowed range, increasing the width of flexure revolute hinges can reduce parasitic displacement and increase the overall compliance of the mechanism, thereby improving operational precision. However, when external disturbances during the micromanipulator operation are minimal or negligible, it is possible to simultaneously increase the width of flexure revolute hinges and decrease the minimum cutting thickness of flexure spherical hinges to reduce parasitic displacement while meeting the working range requirements of the micromanipulator.

Author Contributions

J.R. designed and executed this study. J.R. and H.J. completed the data analysis, paper writing, and revision. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key Project of the Hubei Provincial Department of Education Research Program (grant number D20211401).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The formulas for calculating the compliance matrix of the right-circular flexure revolute hinge are as follows (w, rR, tR):
C δ x , f x = 12 E w 3 0 2 r R x 2 t 2 ( x ) d x + 1 G w 0 2 r R 1 t ( x ) d x C δ y , f y = 12 E w 0 2 r R x 2 t 3 ( x ) d x C δ z , f z = 1 E w 0 2 r R 1 t ( x ) d x C θ x , m x = 12 E w 0 2 r R 1 t 3 ( x ) d x ,   C θ y , m y = 12 E w 3 0 2 r R 1 t ( x ) d x C δ x , m y = C θ y , f x = 12 E w 3 0 2 r R x t ( x ) d x C δ y , m x = C θ x , f y = 12 E w 0 2 r R x t 3 ( x ) d x
where t ( x ) = t R + 2 r R 2 x ( 2 r R x ) ,   ( x 0 , 2 r R ) , E, and G are the Young’s modulus and shear modulus of the material, respectively.
The formulas for calculating the compliance matrix of the right-circular flexure spherical hinge are as follows [30] (rS, tS):
C δ x , f x = 64 0 2 r S x 2 t 4 ( x ) d x + κ ( 1 + υ ) 0 2 r S 1 t 2 ( x ) d x π E C δ y , f y = C δ x , f x C δ z , f z = 4 π E 0 2 r S 1 t 2 ( x ) d x C θ x , m x = 64 π E 0 2 r S 1 t 4 ( x ) d x C θ y , m y = C θ x , m x C θ z , m z = ( 1 + υ ) C θ x , m x C δ x , m y = C θ y , f x = r S C θ x , m x C δ y , m x = C θ x , f y = C δ x , m y
where t ( x ) = t S + 2 r S 2 x ( 2 r S x ) ,   ( x 0 , 2 r S ) , κ is the correction coefficient; here, we take 10/9, υ is the Poisson’s ratio. E is Young’s modulus of the material.

References

  1. Li, H.; Hao, G.; Yu, J.; Deng, Z.; Song, Z. Systematic Approach to the Design of Spatial Translational Compliant Parallel Mechanisms. J. Mech. Eng. 2018, 54, 57–65. [Google Scholar] [CrossRef]
  2. Yang, M.; Du, Z.; Sun, L.; Wei, D. Optimal design, modeling and control of a long stroke 3-PRR compliant parallel manipulator with variable thickness flexure pivots. Robot. Comput. -Integr. Manuf. 2019, 60, 23–33. [Google Scholar] [CrossRef]
  3. Ren, J.; Cao, Q. Dynamic Modeling and Frequency Characteristic Analysis of a Novel 3-PSS Flexible Parallel Micro-Manipulator. Micromachines 2021, 12, 678. [Google Scholar] [CrossRef] [PubMed]
  4. Ren, J.; Li, Q.; Wu, H.; Cao, Q. Optimal Design for 3-PSS Flexible Parallel Micromanipulator Based on Kinematic and Dynamic Characteristics. Micromachines 2022, 13, 1457. [Google Scholar] [CrossRef] [PubMed]
  5. Zhang, S.; Liu, J.; Ding, H. Design and Analysis of a Novel Compliant Mechanism with RPR Degrees of Freedom. Chin. J. Mech. Eng. 2022, 35, 124. [Google Scholar] [CrossRef]
  6. Lin, J.; Qi, C.; Gao, F.; Yue, Y.; Hu, Y.; Wei, B. Modeling and Verification for a Three-Degree-of-Freedom Flexure-Based Planar Parallel Micro Manipulator. J. Mech. Robot. 2023, 15, 041006. [Google Scholar] [CrossRef]
  7. Li, R.; Zhu, Y.; Zhou, S.; He, Z.; Sun, J.; Liu, S. Modeling and control system design of 6-DOF bionic parallel mechanism with compliant modules. Trans. Inst. Meas. Control 2023, 46, 167–182. [Google Scholar] [CrossRef]
  8. Li, H.; Guo, F.; Wang, Y.; Wang, Z.; Li, C.; Ling, M.; Hao, G. Design and modeling of a compact compliant stroke amplification mechanism with completely distributed compliance for ground-mounted actuators. Mech. Mach. Theory 2022, 167, 104566. [Google Scholar] [CrossRef]
  9. Mishra, S.K.; Kumar, C.S. Compliance modeling of a full 6-DOF series–parallel flexure-based Stewart platform-like micromanipulator. Robotica 2022, 40, 3435–3462. [Google Scholar] [CrossRef]
  10. Tian, H.-b.; Ma, H.-w.; Xia, J.; Ma, K.; Li, Z.-z. Stiffness analysis of a metamorphic parallel mechanism with three configurations. Mech. Mach. Theory 2019, 142, 103595. [Google Scholar] [CrossRef]
  11. Wang, N.; Zhang, Z.; Zhang, X. Stiffness analysis of corrugated flexure beam using stiffness matrix method. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2018, 233, 1818–1827. [Google Scholar] [CrossRef]
  12. Zhang, D.; Zhao, L. Research of monolithic XY micro-adjustment mechanism based on compilance matrix. Chin. J. Sci. Instrum. 2014, 35, 269–275. [Google Scholar] [CrossRef]
  13. Li, Y.; Huang, J.; Tang, H. A Compliant Parallel XY Micromotion Stage With Complete Kinematic Decoupling. IEEE Trans. Autom. Sci. Eng. 2012, 9, 538–553. [Google Scholar] [CrossRef]
  14. Ren, J.; Li, Q. Analysis of Compliance and Kinetostatic of a Novel Class of n-4R Compliant Parallel Micro Pointing Mechanism. Micromachines 2022, 13, 1014. [Google Scholar] [CrossRef] [PubMed]
  15. Howell, L.L. Compliant Mechanisms. In 21st Century Kinematics; Springer: London, UK, 2013; pp. 189–216. [Google Scholar]
  16. Ling, M.; Howell, L.L.; Cao, J.; Chen, G. Kinetostatic and Dynamic Modeling of Flexure-Based Compliant Mechanisms: A Survey. Appl. Mech. Rev. 2020, 72, 030802. [Google Scholar] [CrossRef]
  17. Midha, A.; Bapat, S.G.; Mavanthoor, A.; Chinta, V. Analysis of a Fixed-Guided Compliant Beam With an Inflection Point Using the Pseudo-Rigid-Body Model Concept. J. Mech. Robot. 2015, 7, 031007. [Google Scholar] [CrossRef]
  18. Ni, L.; Chen, G.; Hong, K.; Wang, G. Design and development of a compliant piezoelectric microgripper based on three-stage amplification. Microsyst. Technol. 2023, 29, 939–952. [Google Scholar] [CrossRef]
  19. Wadikhaye, S.P.; Yong, Y.K.; Moheimani, S.O.R. Design of a compact serial-kinematic scanner for high-speed atomic force microscopy: An analytical approach. Micro Nano Lett. 2012, 7, 309–313. [Google Scholar] [CrossRef]
  20. Dong, W.; Chen, F.; Gao, F.; Yang, M.; Sun, L.; Du, Z.; Tang, J.; Zhang, D. Development and analysis of a bridge-lever-type displacement amplifier based on hybrid flexure hinges. Precis. Eng. 2018, 54, 171–181. [Google Scholar] [CrossRef]
  21. Wang, X.; Ding, Y.; Zeng, L.; Zhu, C.; Wu, B.; Xu, K. Kinetostatic Modeling of Continuum Delta Robot With Variable Curvature Continuum Joints. J. Mech. Robot. 2023, 15, 031005. [Google Scholar] [CrossRef]
  22. Ren, J.; Wu, H. Kinematic Modeling of a Novel 3-PSS/S Flexible Parallel Micro-turntable: Accuracy Comparison Between Pseudo-Rigid-Body Model-Based and Compliance Matrix Method-Based Models. Iran. J. Sci. Technol. Trans. Mech. Eng. 2023, 47, 2071–2088. [Google Scholar] [CrossRef]
  23. Ren, J.; Lan, A. Compliance Modeling and Kinetostatic Analysis of a Generalized 3-PSS Compliant Parallel Micro-Motion Platform. Micromachines 2024, 15, 354. [Google Scholar] [CrossRef] [PubMed]
  24. Ling, M.; Cao, J.; Howell, L.L.; Zeng, M. Kinetostatic modeling of complex compliant mechanisms with serial-parallel substructures: A semi-analytical matrix displacement method. Mech. Mach. Theory 2018, 125, 169–184. [Google Scholar] [CrossRef]
  25. Arredondo-Soto, M.; Cuan-Urquizo, E.; Gómez-Espinosa, A. The compliance matrix method for the kinetostatic analysis of flexure-based compliant parallel mechanisms: Conventions and general force–displacement cases. Mech. Mach. Theory 2022, 168, 104583. [Google Scholar] [CrossRef]
  26. Arredondo-Soto, M.; Cuan-Urquizo, E.; Gómez-Espinosa, A.; Roman-Flores, A.; Coronado, P.D.U.; Jimenez-Martinez, M. The compliant version of the 3-RRR spherical parallel mechanism known as “Agile-Eye”: Kinetostatic analysis and parasitic displacement evaluation. Mech. Mach. Theory 2023, 180, 105160. [Google Scholar] [CrossRef]
  27. Yu, J.; Bi, S.; Pei, X.; Zhao, H.; Zong, G. Flexure Design: Analysis and Synthesis of Compliance Mechanism Beijing; High Education Press: Beijing, China, 2018; p. 621. [Google Scholar]
  28. Gong, J.; Hu, G.; Zhang, Y. Closed Loop Design Method of Micro-driving Displacement Amplifier Module Targeting for Stiffness. J. Mech. Eng. 2012, 48, 58–64. [Google Scholar] [CrossRef]
  29. Zhu, Z.; To, S.; Zhu, W.-L.; Li, Y.; Huang, P. Optimum Design of a Piezo-Actuated Triaxial Compliant Mechanism for Nanocutting. IEEE Trans. Ind. Electron. 2018, 65, 6362–6371. [Google Scholar] [CrossRef]
  30. Yao, J.; Li, L.; Yang, W.; Zhao, Y. Analytical calculation of compliance matrix for right-circular flexure spherical hinge. Opt. Precis. Eng. 2014, 22, 1857–1863. [Google Scholar] [CrossRef]
Figure 1. (a) Structure of the micromanipulator; (b) structural parameters and coordinate systems of the micromanipulator.
Figure 1. (a) Structure of the micromanipulator; (b) structural parameters and coordinate systems of the micromanipulator.
Micromachines 15 00526 g001
Figure 2. (a) Structure parameters of the right-circular flexure revolute hinge; (b) structure parameters of the right-circular flexure spherical hinge.
Figure 2. (a) Structure parameters of the right-circular flexure revolute hinge; (b) structure parameters of the right-circular flexure spherical hinge.
Micromachines 15 00526 g002
Figure 3. Structure parameters of the equivalent compliant prismatic pair.
Figure 3. Structure parameters of the equivalent compliant prismatic pair.
Micromachines 15 00526 g003
Figure 4. The structure of the branch and the setting of the local coordinate system (a) branch 1; (b) branch 2.
Figure 4. The structure of the branch and the setting of the local coordinate system (a) branch 1; (b) branch 2.
Micromachines 15 00526 g004
Figure 5. Variation of compliance with the parameters: (a) Cδ with l and Ra; (b) Cθ with l and Ra; (c) Cδ with r and w; (d) Cθ with r and w; (e) Cδ with tR and tS; (f) Cθ with tR and tS.
Figure 5. Variation of compliance with the parameters: (a) Cδ with l and Ra; (b) Cθ with l and Ra; (c) Cδ with r and w; (d) Cθ with r and w; (e) Cδ with tR and tS; (f) Cθ with tR and tS.
Micromachines 15 00526 g005aMicromachines 15 00526 g005b
Figure 6. Compliance variations in terms of the scale coefficient e1: (a) Cδ; (b) Cθ.
Figure 6. Compliance variations in terms of the scale coefficient e1: (a) Cδ; (b) Cθ.
Micromachines 15 00526 g006
Figure 7. Compliance variations in terms of the scale coefficient e2: (a) Cδ; (b) Cθ.
Figure 7. Compliance variations in terms of the scale coefficient e2: (a) Cδ; (b) Cθ.
Micromachines 15 00526 g007
Figure 8. (a) Branch number and force loading position; (b) simplification of equivalent stiffness with force F1 acting alone; (c) equivalent spring system.
Figure 8. (a) Branch number and force loading position; (b) simplification of equivalent stiffness with force F1 acting alone; (c) equivalent spring system.
Micromachines 15 00526 g008
Figure 9. (a) The given moving trajectory; (b) the given rotational trajectory.
Figure 9. (a) The given moving trajectory; (b) the given rotational trajectory.
Micromachines 15 00526 g009
Figure 10. Curves of input forces f1,z, f2,z, f3,z and f4,z.
Figure 10. Curves of input forces f1,z, f2,z, f3,z and f4,z.
Micromachines 15 00526 g010
Figure 11. (a) FE model; (b) analytical and FE moving trajectories; (c) analytical and FE rotational trajectory.
Figure 11. (a) FE model; (b) analytical and FE moving trajectories; (c) analytical and FE rotational trajectory.
Micromachines 15 00526 g011
Figure 12. (a) The absolute error of δx and δz in the moving trajectory; (b) the absolute error of θx and θy in the rotational trajectory; (c) the relative error of the motion trajectory.
Figure 12. (a) The absolute error of δx and δz in the moving trajectory; (b) the absolute error of θx and θy in the rotational trajectory; (c) the relative error of the motion trajectory.
Micromachines 15 00526 g012
Figure 13. Variation of CT in terms of the structure parameters: (a) r; (b) w; (c) tR; (d) tS.
Figure 13. Variation of CT in terms of the structure parameters: (a) r; (b) w; (c) tR; (d) tS.
Micromachines 15 00526 g013
Figure 14. Variation of CR in terms of the structure parameters: (a) r; (b) w; (c) tR; (d) tS.
Figure 14. Variation of CR in terms of the structure parameters: (a) r; (b) w; (c) tR; (d) tS.
Micromachines 15 00526 g014
Table 1. Parameters of the transformation matrix in the equivalent prismatic pair.
Table 1. Parameters of the transformation matrix in the equivalent prismatic pair.
Transformation Matrixh1h2h3αβγ
T P 11 P 0.5 a 0 r P + 0.5 t P π / 2 0 π / 2
T P 12 P 0.5 a + b + 2 r P 0 r P + 0.5 t P π / 2 0 π / 2
T P 21 P 0.5 a 0 3 r P + 1.5 t P + c π / 2 0 π / 2
T P 22 P 0.5 a + b + 2 r P 0 3 r P + 1.5 t P + c π / 2 0 π / 2
Table 2. Parameters of the transformation matrix in the mechanism.
Table 2. Parameters of the transformation matrix in the mechanism.
Transformation Matrixxyzαβγ
T P 11 P 0 R a + r S sin θ l + 3 r S cos θ 00π
T P 12 P 0 R a + r S sin θ l + 3 r S cos θ θ0π
T P 21 P 0 R a + r S sin θ r S cos θ θ0π
T P 22 P R a + r S sin θ 0 l + 3 r S cos θ 00π/2
T P 22 P R a r S sin θ 0 l + 3 r S cos θ θ0π/2
T P 22 P R b r S sin θ 0 r S cos θ θ0π/2
Table 3. Structural parameters of 2PRS-2PSS compliant parallel micromanipulator.
Table 3. Structural parameters of 2PRS-2PSS compliant parallel micromanipulator.
ItemValues (mm)ItemValues (mm)ItemValues (mm)
Ra40rR3.75c1
Rb25tR0.5rP3.75
l52.5w5tP0.5
rS3.75a20wP8
tS1b10
Table 4. Material parameters of the flexure hinge.
Table 4. Material parameters of the flexure hinge.
Type of HingeMaterialDensity (kg/m3)Young’s Modulus (GPa)Poisson Ratio
Spherical hingeCuBe280001280.3
Revolute hinge65 Mn80002060.3
Table 5. Error comparison of overall compliance.
Table 5. Error comparison of overall compliance.
ComplianceAnFERelative Error
C δ x , f x (m/N)0.000064830.000066682.77%
C δ y , f y (m/N)0.000005450.000005877.13%
C δ z , f z (m/N)0.000000720.000000693.43%
C θ x , m x (rad/N)0.002720240.002634533.25%
C θ y , m y (rad/N)0.008066180.008056650.12%
C θ z , m z (rad/N)0.008229980.008644324.79%
Table 6. Variation range in the parameters of the micromanipulator.
Table 6. Variation range in the parameters of the micromanipulator.
ItemVariation Range (mm)ItemVariation Range (mm)
l30~60w4~10
Ra35~55tR0.3~0.7
r (rS, rR)2.5~5tS0.7~1.3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ren, J.; Jiang, H. Compliance and Kinetostatics of a Novel 2PRS-2PSS Compliant Parallel Micromanipulator: Modeling and Analysis. Micromachines 2024, 15, 526. https://doi.org/10.3390/mi15040526

AMA Style

Ren J, Jiang H. Compliance and Kinetostatics of a Novel 2PRS-2PSS Compliant Parallel Micromanipulator: Modeling and Analysis. Micromachines. 2024; 15(4):526. https://doi.org/10.3390/mi15040526

Chicago/Turabian Style

Ren, Jun, and Hui Jiang. 2024. "Compliance and Kinetostatics of a Novel 2PRS-2PSS Compliant Parallel Micromanipulator: Modeling and Analysis" Micromachines 15, no. 4: 526. https://doi.org/10.3390/mi15040526

APA Style

Ren, J., & Jiang, H. (2024). Compliance and Kinetostatics of a Novel 2PRS-2PSS Compliant Parallel Micromanipulator: Modeling and Analysis. Micromachines, 15(4), 526. https://doi.org/10.3390/mi15040526

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop