Elastostatic Stiffness Analysis for the US/UPS Parallel Manipulators

The virtual joint method (VJM) cannot calculate the strain energy stored in each rod. In order to solve the problem, a modeling method of the elastostatic stiﬀness was investigated for the UP/UPS parallel manipulators (PMs), taking the example of the 6-SPS PM. The modeling method was based on screw theory, Castigliano’s theorem, and strain energy (where U, P, and S, respectively, denote universal, prismatic, and spherical joints). First, the actuator and constraint wrenches of the mechanism were obtained by screw theory. Second, compact limb stiﬀness matrices were obtained in terms of strain energy and Castigliano’s second theorem. Finally, analytic expressions for the overall stiﬀness matrix of the mechanism and the amplitudes of the actuator force were obtained by adopting the virtual work principle and the balance equation for the mobile platform. All relative errors between the results of the analytical model and the ﬁnite element model are below 2%, which validates the eﬀectiveness of the elastostatic stiﬀness model. The virtual work index was adopted to evaluate the stiﬀness performance of the mechanism, and the results show that the stiﬀness is not only related to position and orientation but also closely related to the directions of external loads. It is also demonstrated that the method has general adaptability for the stiﬀness analysis for the US/UPS PMs, laying the foundation for further reasonable dynamic design and optimization of such manipulators.

As an essential class of the UP/UPS PMs, the 6-SPS PM is an example of successful engineering applications [9], including parallel robotic manipulators [10], positioning devices [11], remote center-of-compliance devices [12], minimally invasive surgeries [13], and haptic interface mechanisms [14]. e 6 degrees of freedom (DOFs) of the Stewart PM based on the 6-SPS PM [15] are well known as a typical example to provide enough orientation capability and high stiffness, and it has been widely used in many industrial applications, for example, the flight simulators [16].
To meet the demand for high precision from the industry, the stiffness analysis of PMs is also essential in the pre-design stage, which should consider the influences of joints and links, and is also the basis of the multi-objective optimization and structural design of PMs [17,18]. Considering the influence and importance of stiffness on the precision of PMs, many progresses have been presented by the academic community for the stiffness characteristics of PMs, which are focused on the various methods of the stiffness modeling and performance evaluation index [19][20][21][22][23].
Although the stiffness of the 6-SPS PM can be numerically simulated by using the finite element analysis (FEA) software, from which the results with high accuracy can be obtained, the process using the FEA software is very timeconsuming, and it has to rebuild the model under different configurations [24,25]. erefore, developing the analytical or semi-analytical models in the pre-design stage is a better choice for estimating the stiffness characteristics of the Stewart PMs in the whole workspace [26]. Some analytical studies about the stiffness analysis of the Stewart PM have been proposed. Using the kinematics error model of the Stewart robot, Ding et al. analyzed the stiffness performance of the Stewart PM equipped with decoupled sensors and actuators [27]. Additionally, the FEA method under the action of (nanoscale) ultrahigh accuracy positioning and the large external load was used to verify the stiffness results. Adam et al. analyzed and evaluated the stiffness matrix of a double Stewart PM truss adopting an analytical method [28]. Based on the principle of linear superposition, Li et al. studied and established the stiffness matrix of the Stewart platform by considering the deformation of the mechanism support [29]. Svinin et al. analyzed the stiffness and stability of the Stewart PM under the condition of internal forces preloading [30]. Huang et al. considered the flexibility factor of the Stewart PM using the law of the conservation of energy and proposed a new conservative stiffness mapping [31]. ey captured the compliance characteristics of the mechanism using an additional stiffness matrix. Moreover, since the rod of two-force only subjected to tension and compression, the virtual joint method is used to analyze the elastostatic stiffness performance of PMs [32,33].
Besides the stiffness modeling, the evaluation of stiffness performance is also essential for ensuring that the performance meets the requirement [34,35]. e stiffness matrix of the PMs is a 6×6 mapping matrix, which can evaluate the relationships between the external wrenched and infinitesimal deformations. To evaluate the stiffness performance of the PMs, the stiffness matrix must be transformed into a quantifiable index [36,37]. e most used stiffness indices at present include the determinant of the stiffness matrix [38], the average value of the eigenvalues [39], the maximum and minimum eigenvalues [40], and the ratio of the minimum and maximum eigenvalues [41]. However, owing to the different units used for position and orientation, there is a common shortcoming of all the above indices; i.e., the dimensions are not uniform, and the index values cannot be well illustrated. e virtual work index (VWI), used in this paper, is reciprocal of the virtual work done by the external load [35]. In contrast to the case for other indices, the units of position and orientation are uniformly used as virtual work within the VWI. Additionally, the VWI can be used to measure the deformation of PMs resisting external load in a specific direction.
In order to solve the problem of the fact that the virtual joint method cannot calculate the strain energy stored in each rod, the main content of this paper is to investigate a modeling method for elastostatic stiffness analysis for the UP/UPS PMs, taking the example of the 6-SPS PM, which is based on the screw theory, Castigliano's theorem, and strain energy method. rough calculating the strain energy stored in each rod, the stiffness matrices of each rod and PMs can be obtained, which have a clear physical meaning and concise expression. e remainder of the paper is organized as follows. Section 2 adopts screw theory, strain energy, and Castigliano's second theorem to derive a general analytical expression for the stiffness matrix of the limb of the 6-SPS PM. e elastostatic stiffness of the 6-SPS PM using the proposed method was analyzed. Section 3 uses the finite element commercial software ANSYS to develop FEA models in different configurations with external wrenches and makes a verification to show the correctness of the theoretical method. Based on the virtual work, Section 4 introduces and evaluates the stiffness performance of the 6-SPS PM by using VWI. Section 5 discusses the results of the simulation. Section 6 presents the conclusions of the study and suggestions for future work.

Process of Elastostatic Stiffness Modeling for the 6-SPS PM.
Based on the strain energy method, Castigliano's second theorem, and the screw system of the limb stiffness matrix for the Stewart PM and the deformation coordination equation (DCE), the global stiffness matrix of parallel mechanism is derived [42]. e calculation process is shown in Figure 1. At the same time, the constrained reaction force exerted on the mobile platform by the supporting limb is obtained. e calculation steps are described in detail as follows: Step 1: using reciprocal screw theory, the constraint wrenches exerted on the mobile platform were defined by limbs $ i1 � ($ i1 $ i2 . . .) Step 2: based on the assumption of small deformations and the mechanics of the materials, we calculate the strain energy are the magnitudes of the limb constraint wrenches and U i a function of f i Step 3: based on Castigliano's second theorem, the elastic deformation at the end of each limb in the direction of the screw axis of the constraint wrenches is calculated as follows: Step 4: we map elastic deformation δ i to constraint wrenches f i , namely, δ i � C i f i ; the stiffness matrix of the limb can be obtained easily, namely, Step 5: the equilibrium equation of the mobile platform in a specific coordinate system was formulated, where W � G F f f Step 6: using the results from steps 4 and 5 along with the virtual work principle, the DCEs are obtained as where D is the infinitesimal twist of the geometric center point of the mobile platform and J i is the mapping matrix between δ i and D Step 7: using the DCEs, the overall stiffness matrix K is

Mathematical Problems in Engineering
Step 8: from matrix K, the magnitudes of the constraint wrenches are obtained, where

Description for the UP/UPS PMs.
In the following, the application of the stiffness modeling method for the UP/UPS PMs is studied. e UP/UPS PMs feature a static platform and mobile platform connected through both links (legs) of type US having constant length and legs of type UPS having controllable distance, as shown in Figure 2.
A typical characteristic of the US/UPS PMS is the same topological structure within the workspace. In the process of parameterization, generalized matrices with the same independent variables and similar structures are analyzed. A class of three-degree-of-freedom (3-DOF) US/UPS PMs, in which the relative motion between the base and platform is spherical, has been proposed, where mechanisms belonging to this class have at least three legs of type US and only three legs of type UPS [43]. A class of 6-DOF US/UPS PMs, in which the relative motion between the base and platform is general, has been proposed, where mechanisms belonging to this class have legs of type UPS only [44].
e diagram of the model and the coordinate system of the 6-SPS PM are shown in Figure 3. Six identical linkages of the 6-SPS PM connect the mobile platform and static platform at B i and A i through a spherical joint. Both the mobile platform and static platform are semi-regular hexagons. e circumferential radius of the mobile platform and static platform is, respectively, denoted by R b and R a .    Mathematical Problems in Engineering e center angles of the long sides A 1 A 2 and B 2 B 3 are β 0 . As is well known, there are six SPS limbs for 6-SPS PM.

Stiffness Matrix of Each Limb for the 6-SPS PM.
A force analysis diagram of the mobile platform of the 6-SPS PM is shown in Figure 4. e motion pairs of the limb from the base to the mobile platform are S pairs, P pairs, and S pairs. Each limb of the 6-SPS is subjected only to an axial force, resulting in the limb rod generating a tensile or compressive elastic deformation. According to the screw theory, there is a constrained wrench exerted on the SPS limb of the mobile platform along the center for the two S pairs. erefore, the SPS rod can regard each limb as a two-force rod only subject to tension and compression. Force analyses of the SPS limb and mobile platform are shown in Figure 5: e strain energy of the SPS limb is given as where f Ni is the axial force of the i th limb, L i is the length of the i th limb, E i is the elastic modulus of the i th limb, and A i is the cross-sectional area of the i th limb. e total strain energy of the mechanism can be obtained as e elastic deformation of the limb end along the direction of the constraint wrench screw system is obtained from Castigliano's second theorem and the relation e screw flexibility matrix of the constraint system of the supporting limb is obtained as And the screw stiffness matrix of the constraint system is further obtained as According to the giveCcn position and attitude, the compliance matrix of the SPS limb is C i � L i /E i A i , as expressed by equation (5). e stiffness matrix is is the unit screw of the constraint wrench vector direction applied by the ith limb to the mobile platform. S i is the unit vector along the path of the A i B i rod, b i � pB i , and J i � ($ r i ). S i and b i are expressed in the definite coordinate system o−xyz. e overall stiffness matrix K i can thus be obtained from equation (6).

Overall Stiffness Matrix.
e mobile platform equilibrium equation is where W � [F M] is an extra wrench acting on the mobile platform and $ r ij is the unit screw of the j th constraint force that the i th limb exerts on the mobile platform.
Under the action of the constraint wrench screw, elastic deformation at the end of each limb of the mechanism is generated along the axis of the constraint wrench screw. e deformation coordination between the infinitesimal displacements of the mobile platform under the action of external forces is expressed as where J i � ($ r i1 . . . $ r ij ) T and D is the tiny twist of the geometric center point of the mobile platform. According to the knowledge of material mechanics, we have en, by substituting equations (8) and (9) into equation (7), we have where n is the number of limbs of the Stewart PM. Equation (11) is the general analytical expression of the overall Mathematical Problems in Engineering stiffness matrix and is more concise. e analytical model of the stiffness matrix is given from the perspective of strain energy and Castigliano's second theorem. (10), we further obtain

Amplitude of the Driving Force. From equation
en, by substituting equations (8) and (12) into equation (9), we obtain the inverse constrained solution equation as

Comparison with the Virtual Joint Method.
According to the calculation requirements of the VJM, the rod is regarded as a rigid rod, supported by a virtual spring at both ends of the rod. e limb of the SPS can be treated as a two-force rod, and the virtual spring is a one-way spring along the axial direction of the rod. e stiffness coefficient of the spring is e overall stiffness matrix K is Formally, equation (15) is the same as equation (11). Due to the fact that the VJM method cannot calculate the strain energy stored in each rod, the method proposed in this paper can calculate the strain energy stored in each rod and overall parallel manipulators. A new approach to solving the problem of each component for the stiffness performance of the PM is proposed.

Comparison of the Proposed Model and
FEA Model e finite-element commercial software ANSYS is used to develop a unified FEA model and verify the correctness and feasibility of the method investigated in this paper. Analytical expression theory and the FEA simulation model are implemented.
Physical parameters of the 6-SPS PM are given in Table 1, where R is the radius of the circumferential circle of the fixed platform, r is the radius of the circumferential circle of the mobile platform, d is the diameter of the cross section of the limb rod, and E and G are, respectively, the elastic modulus and shear modulus of the material.
At the same time, the five poses and five groups of external forces and their couples acting on the mobile platform of the 6-SPS PM are given in Table 2.
Five groups of data to make a theoretical analysis calculation and FEA simulation calculation for the established analytical model are taken. e results obtained using the two methods are compared one by one. e theoretical calculation values, FEA simulation values, and relative errors of each group are listed in Table 3. It shows that the maximum relative error is 1.581% for the fifth group. All relative errors between the results of the analytical model and the workbench finite element model are below 2%, which validates the effectiveness of the proposed elastostatic stiffness model of the 6-SPS PM. Table 3: Deformation of the geometric center point of the mobile platform of the 6-SPS PM (comparison of the results of the FEA and theoretical models).

Pose
Method

Stiffness Index of the 6-SPS PM
e VWI is used to evaluate the stiffness of PMs. In contrast with other indices, this index unifies the units of displacement and angle relating to virtual work. By correlating the value of the stiffness index with the direction of the external load, the ability of the mechanism to resist the external load in one direction is measured.
Supposing that a unit wrench W is imposed on the mobile platform, the virtual work done by the unit wrench can be expressed as e infinitesimal twist is δξ � K −1 W , and the VWI of the stiffness for the 6-SPS can thus be expressed as e VWI of the stiffness for the 6-SPS PM is a measure of the ability of the mechanism to resist deformation under a given external wrench. It allows the stiffness matrix to be transformed into a single index value. To calculate the value of the VWI, the virtual work completed by the force and moment must have the same units.  different directions. Figure 7(a) shows the results for a vertical load. e stiffness of the 6-SPS PM is best for the central position, and the stiffness performance of the mechanism decreases rapidly with deviation from this position.

Results and Discussion of the Simulation
Similarly, we can see from Figure 7(c) that the stiffness performance of the 6-SPS PM is best when the position is not deflected, and the stiffness performance of the mechanism decreases rapidly with deviation from this position. Figures 7(b) and 7(d) show that the optimal position for the 6-SPS PM stiffness under random loading is not at the central position of the mechanism. e results show that the stiffness performance of the 6-SPS PM is closely related to the direction of the external load and different under different external loads. It is shown that the elastic deformation of the 6-SPS PM is closely associated with the stiffness, configuration, and external loads.

Conclusion
In this work, a modeling method is investigated for elastostatic stiffness analysis for the UP/UPS PMs, taking the example of the 6-SPS PM, which is based on the screw theory, Castigliano's theorem, and strain energy method. e global stiffness matrix was derived. e virtual work index was adopted to evaluate the stiffness of the 6-SPS PM, and the VWI was used to measure the deformation of the 6-SPS PM in resisting external loads in specific directions. e main conclusions of the study are as follows: (1) Although the theoretical method used in this paper is much simplified and has many limitations, the results show that the relative errors of the theoretical model and finite element model are less than 2%, which verifies the correctness of the analytical model. During the predesign stage, the stiffness analysis model used in this paper can be regarded as an alternative to the finite element model in evaluating the stiffness performance of the 6-SPS PM. (2) e VWI was used to evaluate the stiffness performance of the 6-SPS PM. e results show that the stiffness of the mechanism is related to the direction of the action of the external wrench. (3) In contrast, VJM cannot calculate the strain energy stored in each rod, and the described method can be applied to not only the 6-SPS PM but also other US/UPS PMs, such as two-DOF PMs, three-DOF PMs, and other six-DOF PMs. e modeling method for the elastostatic stiffness of the US/UPS PMs was derived in this paper. e limitation of the model is that the flexibility of a joint is ignored. e study did not establish an experimental model for verification and did not consider the problem of dynamics. ese issues need further study in the future.
Data Availability e MATLAB and ANSYS files data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.