Dynamic Modeling and Analysis of 5-PSS/UPU Parallel Mechanism with Elastically Active Branched Chains

To study the characteristics of the 5-prismatic–spherical–spherical (PSS)/universal–prismatic–universal (UPU) parallel mechanism with elastically active branched chains, the dynamics modeling and solutions of the parallel mechanism were investigated. First, the active branched chains and screw sliders were considered as spatial beam elements and plane beam element models, respectively, and the dynamic equations of each element model were derived using the Lagrange method. Second, the equations of the 5-PSS/UPU parallel mechanism were obtained according to the kinematic coupling relationship between the active branched chains and moving platform. Finally, based on the parallel mechanism dynamic equations, the natural frequency distribution of the 5-PSS/UPU parallel mechanism in the working space and elastic displacement of the moving platform were obtained. The results show that the natural frequency of the 5-PSS/UPU parallel mechanism under a given motion situation is greater than its operating frequency. The maximum position error is − 0.096 mm in direction Y, and the maximum orientation error is − 0.29° around the X-axis. The study provides important information for analyzing the dynamic performance, dynamic optimization design, and dynamic control of the 5-PSS/UPU parallel mechanism with elastically active branched chains.


Introduction
The sea surface recovery platform is constantly moving and swaying, thereby leading to the tilting of the recovery spacecraft and causing the spacecraft recovery mission to fail. Currently, only NASA has successfully developed a surface recovery platform for spacecraft; the first stage of the "Falcon 9" carrier rocket landed vertically on the sea platform in April 2016 [1]. For some dynamic balancing devices, only the condition of the vertically launched instantaneous missile is suitable [2,3]. However, during the spacecraft recovery process, the dynamic balancing device needs to be constantly balanced under a large variable load impact. Therefore, the existing dynamic balancing device of the ship-borne missile launch is not suitable for the surface recovery platform of spacecraft. Parallel mechanisms are increasingly used in applications where precision is of great importance [4][5][6]. To take advantage of the parallel mechanism by applying it to the surface recovery platform, a 5-prismatical-spherical-spherical (PSS)/ universal-prismatical-universal (UPU) parallel mechanism is proposed as a dynamic balancing device [7] in this paper. To reduce the load on the device and improve its driving ability, all the drives are placed on the frame. The position adjustment of five degrees of freedom of the moving platform can be realized through the real-time drive control of each active branch, and the dynamic balance can be maintained.
Under the condition of high speed and heavy load, the elastic deformation of each component will have a certain influence on the motion accuracy of the moving platform and the positioning accuracy of the load [8][9][10][11][12][13]. Therefore, to reduce the influence of elastic deformation on the moving platform and improve its output accuracy, the research on the elastic dynamics modeling and dynamic characteristics of the 5-PSS/UPU mechanism are absolutely necessary [14][15][16]. Liu et al. [17] established an elastic dynamics modeling of a flexible 3-RRS parallel robot using the simplified KED (kineto-elastodynamics) method and analyzed its dynamic characteristics in detail. Fattah et al. [18] obtained the whole mechanism dynamic equation of the 3-RRS parallel mechanism using the natural orthogonal complement method, and the influence of the output precision of the moving platform with flexible links was studied. Zhao et al. [19] established the elastic dynamic equations of each moving member based on the idea of substructure, and the dynamic equations of threedegrees-of-freedom translational parallel mechanism was obtained using the displacement coordination. Xie et al. [20] considered the parallel moving platform and active branched chains as a spatial beam element modeling, and an overall dynamic equation was established based on the motion constraints. Furthermore, the Newmark numerical method was used, and the elastic dynamic model was solved discretely. Zhao et al. [21] investigated the elastic dynamic characteristics of the 6-PSS and 8-PSS parallel robots, and the result shows that the redundant 8-PSS parallel robots have higher natural frequencies and better dynamic characteristics. Shan et al. [22] established the elastic dynamics modeling of a novel 2 (3HUS + U) parallel hip joint simulator, and the natural frequency and stiffness of the mechanism in its working space were analyzed. As mentioned above, the investigations on the parallel mechanism are mainly limited to the simple planar or three-degrees-of-freedom parallel mechanism, and only a few studies have been reported on the dynamic characteristics of parallel mechanism with complex structure, especially the five-degrees-of-freedom parallel mechanism.
In this study, a novel five-degrees-of-freedom parallel mechanism-5-PSS/UPU parallel mechanism-is proposed [23,24], considering the moving platform, and the active branching chains and screw slider as the spatial beam element and plane beam element modeling, respectively. Next, according to the kinematic coupling relationship of each component, the dynamic equation of the parallel mechanism is constructed. Finally, the natural frequency distribution of the 5-PSS/UPU parallel mechanism in the working space and the elastic displacement of the moving platform are obtained.

Dynamic Equations
The basic structure of the parallel mechanism is shown in Figure 1. It consists of a fixed base and a moving platform connected using six branched chains. The six branches include five PSS joint branches and one UPU joint branch. The two ends of five PSS joint branch chain are respectively connected to the moving platform and linear module through two spherical joints, and the linear module is fixed on the base. The UPU joint branch chain is connected to the geometric center of the moving platform and fixed base through a universal joint. The power input is the moving pair of five PSS joint branches,and the UPU joint branch chains only provide constraints to institutions.
As shown in Figure 2, the rigid moving pentagon platform is a 1 a 2 a 3 a 4 a 5 , and the radius of the circumscribed circle is r. The fixed pentagon base B 1 B 2 B 3 B 4 B 5 is assumed rigid, and the radius of the circumscribed circle is R. Each active branched chain contains a link and a linear motor. The active branched chains are connected to moving platform a 1 a 2 a 3 a 4 a 5 through spherical joints and coupled to the slider through spherical joints at b i , where i = 1, 2, 3, 4, 5. The UPU branched chain is connected to moving platform a 1 a 2 a 3 a 4 a 5 through the universal joint and coupled to the fixed base through the universal joint. The length of link a i b i is L.
The basal Cartesian coordinate frame, designated as the o−xyz frame, is fixed at the center of the base platform, with the z-axis pointing vertically upward and x-axis pointing towards joint B 1 . Similarly, a coordinate frame p−x p y p z p is assigned to the center of the moving platform, with the z p -axis normal to the platform and x p -axis pointing towards jointa 1 . The local frame w i −x i y i z i (i = 1, 2, 3, 4, 5) is fixed on the base platform joint B i , with y i -axis pointing from B i to b i and x i -axis pointing vertically to x-axis. The parallel mechanism has five degrees of freedom and five linear motors to drive the actuated joints.
To analyze the characteristics of the 5-PSS/UPU with elastically active branched chains, the dynamics modeling of the parallel mechanism was investigated based on the finite element theory. To ensure the accuracy of the analysis and reduce the complexity of modeling, the following assumptions are made: 1. The deformation of the flexible components is very small and can be regarded as a small elastic deformation. Thus, the actual movement of the member can be regarded as a linear superposition of the elastic and rigid displacement. 2. Considering that the UPU chain is composed of an electric cylinder and universal joints, its elastic deformation is relatively small compared to the active branch chains. Thus, considering the moving platform, peripheral bracket, and UPU branched chains as rigid, the torsional deformation of active branched chains (including flexible link and linear motor) is mainly considered.

Flexible Link Dynamic Equation
As shown in Figure 3, the flexible link is considered a spatial flexible beam element for elastic dynamic modeling. δ li = [δ 1 , δ 2 , . . . , δ 18 ] T represents the vector of generalized coordinates of beam elements, where δ 1 −δ 3 and δ 10 −δ 12 , δ 4 −δ 6 and δ 13 −δ 15 , and δ 7 −δ 9 and δ 16 −δ 18 represent the axial or transverse displacements, rotary angles, and curvatures at nodes b i and a i , respectively. It is supported that a spatial flexible beam element is subjected to axial, lateral, and torsional deformational. A point in the element has elastic displacement in the direction of u i v i w i -axes. According to the deformation characteristics and requirement of the flexible components, the lateral elastic, axial elastic, and elastic angular displacements around w i -axis of beam elements are expressed using the quintic Hermite, linear, and cubic interpolation functions, respectively. Next, the functions can be obtained based on the set of boundary conditions of the flexible beam element as follows: where Φ ui , Φ vi , Φ wi , and Φ ϕi are interpolation vectors, and the functions of w. The expressions are specified as follows [25]: where n i (i = 1, 2, . . . , 10) is the type function of displacement of elements and w is the axial displacement of elements.
The absolute acceleration at a random point on the beam element is considered to be the sum of the acceleration of the movement of rigid body and the acceleration of elastic deformation. Hence, the velocity of the random point with the coordinates of w on beam element is shown as follows: (1)  where Φ aui , Φ avi , and Φ awi are the absolute velocities of a given point on beam element along the u-, v-, and w-axes, respectively; Φ rui , Φ rvi , and Φ rwi are the velocities of the moving rigid body along u-, v-, and w-axes, respectively; Φ ui , Φ vi , and Φ wi are the velocities of elastic deformation of a given point on beam element along the u-, v-, and w-axes, respectively; Φ aϕi , Φ rϕi , and Φ ϕi are the absolute angular velocity, angular velocity of the rigid body, and angular velocity of elastic deformation around w-axis of a given point on beam element, respectively; and δ rl represents the vector of generalized coordinates of the rigid body. The kinetic energy of the spatial flexible beam element is expressed as follows: where L and ρ ab are the length and mass density of beam element, respectively; I p is the polar moment of inertia of cross sections of beam element about w-axis; S ab is the cross-sectional area of beam element; and δ w i is the w iaxial displacement. Moreover, M abi is the function of the mass distribution of beam element: Ignoring the shear deformation of the beam element and coupling between the axial displacement and lateral displacement, the potential energy of the branched component is expressed as follows: where E ab and G ab are the elastic modulus of tension/ compression and shearing modulus of elasticity of the material, respectively. I abiu and I abiv are the principal moments of inertia of cross sections of beam element to u-and v-axes, respectively. Moreover, K abi is the function of rigidity of beam element: where ∂δ wi . By substituting Eqs. (4) and (6)  where F abi is the array of the generalized force of external load of beam element, F ebai is the array of forces on the studied beam element exerted by other beam elements connected with the studied one, and δ ′′ rli is the acceleration of rigid bodies.
Because two end nodes of the branched link are spherical joints, the curvature in three directions is zero. Consequently, the generalized coordinates ( Figure 4) of the branch can be expressed in the global coordinates as follows: then where R li is the transformation matrix of the element coordinate system to the global coordinate system.
The element equations are expressed in terms of the O-XYZ frame as follows: where

Linear Motor Dynamic Equation
Because the linear module only performs linear motion in a fixed plane, the linear motor is considered a plane beam element model for elastic dynamic modeling. The linear motor is divided into two units. The joint of the coupling and the screw at the position slider b i is the first unit, and the position of the slider at the end of the screw bearing is the second unit. Each node has then three degrees of freedom, and generalized coordinates ( Figure 5) of the linear motor can be expressed as follows: (12) With the boundary constraints, the transverse displacements and rotary angles of the linear motor in the element coordinates can be expressed as follows: Considering the dynamic accuracy design requirements of the moving platform, the lateral elastic displacement and axial elastic displacement of elements are expressed using the quintic Hermite interpolation function and linear interpolation function, respectively.
where N Lx and N Ly are the vectors of interpolation polynomials. The expressions are expressed as follows [26]: where y i and L sj are axial displacement and axial length of the element, respectively.
According to the kinetic energy of the linear motor element unit, the mass matrix can be obtained as follows:  where ρ L , S L , and L sj are the mass density, cross-sectional area, and slider displacement of linear motor element, respectively.
With the boundary conditions of the element unit, the overall mass matrix M Li9×9 of the linear motor is assembled.
The shear deformation of the beam element, and the coupling between the axial displacement and lateral displacement are ignored. Consequently, the function of rigidity of beam element K Lj6×6 can be expressed as follows: where E L is the elastic modulus of tension/compression of the material, and I Lx is the principal moments of inertia of the cross sections of linear motor element along x i -axes.
With the boundary conditions of the element unit, the overall stiffness matrix K Li9×9 of the linear motor is assembled.
The linear motor dynamic equation in the terms of local coordinate can be obtained as follows: where F Li is the array of generalized force of external load of linear motor element, F Lei is the array of forces on the studied linear motor element exerted by other elements connected with the studied one, and ü Lri is the acceleration of rigid bodies.
In the global coordinates, the generalized coordinates of the linear motor can be expressed as follows: , and R Li is the rotation matrix of the base coordinate system to the global coordinate system of the linear motor element.
Next, the element equations are expressed in terms of the O-XYZ frame as follows:

Kinematic Constraint Equations
As shown in Figure 6, the elastic deformation of the moving platform is much smaller than that of the branched chains; hence, it can be regarded as a rigid body [27]. is (α, β, γ ) , and the coordinate of point p in coordinate system O-xyz is (p x , p y , p z ) T . The transformation matrix of coordinate system p−xyz to the global coordinate system O−XYZ is R o p , the transformation matrix of coordinate system p ′ − x ′ y ′ z ′ to the global coordinate system O−XYZ is R o p ′ , and the transformation matrix of coordinate system p ′ − x ′ y ′ z ′ to the coordinate system p−xyz is R p p ′ . Based on small deformation, the approximate expression of transformation matrix R p p ′ is defined as follows: Next, the transformation matrix R o p can be expressed as follows: In this case, the elastic line displacement at the end of the branched chains a i can be obtained as follows: where J pi = I T pi 0 I .
The elastic displacement of joint b i between the linear motor and slider can be expressed as follows: (25) a i = J pi p, (26) where p is the lead for the screw and R si = 1 0 p 2π . The elastic displacement of the slider is the same as the elastic displacement of the end of the branched chains in y i axial direction of w i -x i y i z i , �S i = δ 3i .

Dynamic Model of 5-PSS/UPU Parallel Mechanism
Based on the above analysis, the generalized coordinates of the 5-PSS/UPU parallel mechanism system are as follows: In the generalized coordinates of the system, the generalized coordinates of the branched chains can be expressed as follows: By substituting Eq. (28) into Eq. (11), the mass matrix and function of the rigidity of beam element in the generalized coordinate can be obtained as follows: In the generalized coordinates of the system, the generalized coordinates of the linear motor can be expressed as follows: where R Li * = E 12×12 0 12×12 . By substituting Eq. (30) into Eq. (19), the mass matrix and function of linear motor in the generalized coordinate can be obtained as follows: The generalized coordinates of the slider can be expressed as follows: where R Si * = 0 1×3 R si R Li 0 1×18 . In this case, the mass matrix of slider in the generalized coordinate can be obtained as follows: Therefore, based on the above analysis, the total mass matrix and stiffness matrix of the single branch of the 5-PSS/UPU parallel mechanism in the generalized system coordinates are expressed as follows: (27) The moving platform is considered a rigid body. Therefore, in coordinate system O-xyz, the kinematic equation of the moving platform can be obtained as follows: where M p is the mass matrix of the moving platform, T p is the array of forces exerted on the studied moving platform by branched chains with the studied one, T ep is the array of generalized force of external load of the moving platform, and ü rp is the acceleration of rigid bodies.
It is supposed that the U is the array of generalized coordinates of the 5-PSS/UPU parallel mechanism system, which is expressed as follows: The dynamic equation of the parallel mechanism with rigid and flexible couplings can be obtained by assembling differential equations of motion of kinematic chains, and combining the kinematic and dynamics constraint equations: where U is the array of generalized coordinates of the 5-PSS/UPU parallel mechanism system. M and K are the mass matrix and rigidity matrix, respectively. C = 1 M + 2 K is the damping matrix of the system, and 1 and 2 are the Rayleigh damping coefficients. F is the array of generalized forces. Ü r is the array of the accelerations of rigid bodies, which can be obtained through the kinematic analysis of the rigid body of the 5-PSS/UPU parallel mechanism. (37)

Analysis of Natural Frequency Characteristics
To prevent the resonance of parallel mechanism caused by impact load [28] and reduce its effects on the dynamic accuracy of the moving platform, it is necessary to analyze the natural frequency characteristics of the 5-PSS/ UPU parallel mechanism. Table 1 shows the parameters of the 5-PSS/UPU parallel mechanism. Neglecting the damp term from Eq. (37), we obtain the main vibration equation of the system and its characteristic equation as follows: where ω n denotes the natural frequency. The natural frequency with respect to the 5-PSS/UPU parallel mechanism can be achieved using Eq. (38) during the simulation, which is very significant to realize the dynamic decoupling of the parallel mechanism [29].
Next, we can obtain the distribution of the natural frequency of the 5-PSS/UPU parallel mechanism in the position space at z p = 0 and the distribution under the attitude working space. Figure 7 shows the natural frequency of the 5-PSS/UPU parallel mechanism positional space with consideration of elastically active branched chains candidates. It can be observed that the variations in each natural frequency in the position working space are [2.19, 20.89] Hz, [23.49, 64.64] Hz, and [65.31, 100.33] Hz. The first-order natural frequency and the third-order natural frequency are relatively large at the center position (X = 0; Y = 0), and the natural frequency is relatively small at the working edge position. The second-order natural frequency is small at the center position (X = 0; Y = 0), and the natural frequency is large at the working edge position [30]. Figure 8 shows the natural frequency of the 5-PSS/UPU parallel mechanism attitude space with consideration of elastically active branched chains candidates. It can be observed that the variations in each natural frequency in the attitude working space are [3.89, 20.89] Hz, [23.49, 69.48] Hz, and [37.49, 100.33] Hz. The first-order natural frequency and the third-order natural frequency are relatively large at the center position ( α = 0 o ; β = 0 o ), and the natural frequency is relatively small at the working edge position. The second-order natural frequency is small at the center position ( α = 0 o ; β = 0 o ), and the natural frequency is large at the working edge position.
In addition, it can be observed through the above natural frequency numerical analysis that the natural frequency of the 5-PSS/UPU parallel mechanism in the moving working space is greater than its operating frequency of 0.5 Hz; hence, the system does not cause resonance.

Elastic Displacement of Moving Platform
During simulation, the moving platform is set to move on a trajectory given as follows: Using the Newmark method [31] to solve Eq. (36), the displacement errors of the moving platform in the global frame O-XYZ are obtained. Figure 9 illustrates the displacement errors of the moving platform in the global frame O-XYZ, where ε X , ε Y , and ε Z refer to the position errors in directions X, Y, and Z, respectively. ε γ and ε β refer to the orientation errors around X-axis and Y-axis, respectively. With these preconditions, it is found that the maximum position error is − 0.096 mm in direction Y, and the maximum orientation error is − 0.29° around X-axis. Hence, the influence of elastic links should not be ignored, especially when the moving platform moves at a relatively high speed. It is necessary to analyze the displacement errors of the flexible 5-PSS/UPU parallel mechanism.

Conclusions
A novel five-degree-of-freedom parallel mechanism is proposed as the dynamic balance device in the spacecraft recovery research filed because of its simple structure, better dynamic, and stiffness qualities.
According to the method of the finite element theory and Lagrange, the dynamics modeling of flexible branched links and linear motors are investigated, and the 5-PSS/UPU parallel mechanism according to the kinematic coupling relationship between each components of parallel mechanism is obtained.
The natural frequency and accuracy of the moving platform of 5-PSS/UPU parallel mechanism by considering elastically active branched chains are investigated through numerical simulation. The results show that the natural frequency of the system in the moving working space is greater than its operating frequency, the maximum position errors in direction of X-axis, Y-axis, and Z-axis are − 0.084 mm, − 0.096 mm, and 0.018 mm, respectively, and the maximum orientation errors around X-axis and Y-axis are − 0.29° and 0.26°, respectively. It is of great significance to improve the kinematics performance of the parallel mechanism.