Open Access
Issue
Mechanics & Industry
Volume 25, 2024
Article Number 5
Number of page(s) 19
DOI https://doi.org/10.1051/meca/2023044
Published online 14 February 2024

© Y. Sheng and X. Chen Published by EDP Sciences 2024

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Spatial parallel mechanism [1] has some competitive advantages compared with other mechanisms, such as more compact mechanical structure, smaller kinematic inertia, higher motion accuracy and good isotropy [2], it has been successfully applied in the field of coordinate measuring mechanism. The use of spatial parallel mechanism improves the measurement accuracy and mechanism stability of coordinate measuring mechanism [35]. Spherical joints are extensively used linking devices with the advantages of higher transmission efficiency, but existence of clearance between joint elements is inevitable because of assembly requirements. Spherical clearance joints exert a great influence on the dynamic behavior of spatial parallel coordinate measuring mechanism, reducing measurement accuracy and causing unpredictable vibration [6]. Studying the influence of clearance joints on spatial parallel coordinate measuring mechanism can be demonstrated by the dynamic behavior under different research parameters [7]. Based on this, it is critical to establish a dynamics modeling and analysis method applied for spatial parallel coordinate measuring mechanism, and the method can also enable the research of the effect of spherical clearance joints on dynamic behavior of spatial parallel mechanisms.

In the past, research about dynamic of spatial parallel mechanism is imperfect, and research direction is mainly focused on the plane multi-link mechanism [811]. Li et al. [12] developed a comprehensive methodology for predicting wear characteristic in crank slider mechanism with multiple clearance joints, and found a result that changing boundary conditions have large effect on wear characteristics of the mechanism. Bai et al [13] researched the influence of joints with axial and radial clearances on the dynamic response of planar mechanical systems, and discussed dynamic response of hybrid clearance joints to crank-slider mechanism with different examples. Chen et al. [14] compared the experimental data with numerical simulations to investigate the effect of clearance joints on the dynamical characteristics of planar multibody systems, and verified the correctness of modeling and analysis methods. Wu et al. [15] analyzed the system dynamics of planar crank-slider mechanism with multiple clearance joints by studying the correlation dimension and bifurcation action, and found that the response of system has different sensitivities to different clearance joints. Wang et al. [16] explored dynamic behavior of planar four-bar mechanism considering different clearance joint mounting positions and number of clearance joints, and inferred that the larger the number of clearance joints, the worse the stability of mechanism.

The research on chaotic characteristics of spatial parallel mechanisms is immature, and the current research mainly focuses on the dynamic response of spatial parallel mechanisms [1720]. Erkaya et al. [21] explored influence of vibration generated by clearance joint on chaos of mechanism, concluded that vibration reduces dynamic performance and making rod flexible can reduce the negative impact which clearance bring to mechanism. Hou et al. [22] explored dynamic response of four-bar mechanism with rigid rod, flexible rod and flexible socket, and inferred that flexible rod and flexible socket can weaken adverse impact on the system produced by clearance joint. Wang et al. [23] built dynamic model of parallel mechanism with spherical clearance and flexible input rod, found a result that making the rod flexible is very necessary to study the influence of clearance joint on mechanism. Zhang et al. [24] analyzed the dynamics of a new parallel mechanism using implicit integration algorithm and Baumgart stabilization technology, and verified the validity of numerical results of the dynamic response using ADAMS simulation software. Chen et al. [25] used the principle of virtual work to analyze dynamics of a new type of over constrained parallel mechanism, and verified the validity of dynamic simulation results through the physical model simulation in Simscape and rigid flexible coupling model simulation in Adams.

In this paper, because it is immature and imperfect in the field of researching dynamic response and chaotic characteristics of spatial parallel mechanisms at the same time, a universal dynamics modeling and analysis method applied for spatial parallel coordinate with multiple spherical clearance joints was proposed. In Section 2, motion collision process of spherical clearance joints is described in detail, kinematic model of spherical clearance joint is established, and different force models are introduced to analyze contact force at the contact collision point. In Section 3, kinematic model of spatial parallel coordinate measuring mechanism is established, and accompanying motion of mechanism is explained, and clearance dynamic model based on Lagrange multiplier technology is established. In Section 4, the correctness of clearance dynamics modeling method is verified by comparison with Adams virtual simulation, and taking clearance value and friction coefficient as single input variables, effect of spherical clearance joints on response and chaotic characteristics of mechanism is studied.

2 Motion description of spherical clearance joint

2.1 Kinematic model

Spatial parallel coordinate measuring mechanism (SPCMM) is made up of three identical kinematic chains connected to a mobile platform, SPCMM can realize data measurement of three degrees of freedom (one translation measurement and two rotation measurement). Every branch consists of a driving rod and a driven rod, driving rod is fixed connection with the frame through the translational joint, the driving rod is fixed connection with the driven rod through spherical joint, and driven rod is fixed connection with the terminal platform through rotating joint.

Figure 1 is the three-dimensional virtual model of SPCMM. Figure 2 is structural sketch of SPCMM, Pi represents translational joint, Si represents spherical joint, Ri represents rotating joint, OA-XYZ is world coordinate system fixed connection with the frame which is expressed as {A}, its vertical distance from point P is a1, oi is object coordinate system fixed connection with other moving components which is expressed as {oi}, and {oi}(i = 1 ~ 6) is fixed at the centroid position of each component, and o7 is established at the center position of R1 and R2 on the terminal platform. Length hi(i = 1 ~ 3) is the initial installation position of the translational joint from point P.

Spherical joints S1, S2  and S3 are clearance joints, with the ball socket of the ball joint fixedly connected to the driving rod, and the ball joint ball fixedly connected to the driven rod. Existing clearance makes the position constraint of spherical joint invalid, and the ball can get a freely motion in six degrees of freedom, therefore, the modeling of spherical clearance joint can be seen in Figure 3. Coordinate system {A} is world coordinate system, Point Zi and Zj are center points of the socket and ball, R and r is radius of socket and ball. Vector Ann and Att are normal vector and tangential vector of plane where the contact point is located. Vector eij is eccentric vector of center points in spherical joint. Vector FN and FT are force at the point of impact.

The position coordinate about origin of object coordinate system fixed connection with driving rod in world coordinate system is

(1)

The Euler matrix of object coordinate system in world coordinate system is written as

(2)

The position coordinate about center point of ball and socket in world coordinate system is written as follow

(3)

Vector opzi and opZj are position coordinate of center point of ball and socket in object coordinate system and they can be expressed as

(4)

Symbol L1 and L2 are the length of driving rod and driven rod.

The eccentric vector Aeij of spherical joint is written as

(5)

The amplitude of eccentric vector Aeij at clearance joint is expressed as

(6)

The relative impact velocity of center points between ball and socket is

(7)

The normal vector of the plane where the impact point Ci and Dj are located is

(8)

The clearance value of the spherical joint is

(9)

Pseudo-penetration depth between spherical joint elements is expressed as

(10)

According to the value of pseudo-penetration depth, the specific motion state of spherical joint elements can be judged as follow

(11)

when spherical joint elements are in contact and collision state, the position coordinate of potential collision points Ci and Dj in world coordinate system are

(12)

Velocity vector of collision point in world coordinate system is

(13)where, AVZi and AVZj are velocity vector of center point in world coordinate system, they can be written as

(14)

where, AWoi is the Euler velocity vector, ArZi is position vector of center point Zi in object coordinate system which transformed from object coordinate system to world coordinate system, AVoi is velocity vector of the origin of object coordinate system in world coordinate system.

Normal component of relative collision velocity is

(15)

Tangential component of relative collision velocity is

(16)

Amplitude of tangential component of relative collision velocity can be expressed as

(17)

So tangential vector of plane where impact point Ci and Dj are located is

(18)

thumbnail Fig. 1

Characteristic map of SPCMM.

thumbnail Fig. 2

The kinematic scheme of SPCMM.

thumbnail Fig. 3

Spherical clearance joint model.

2.2 Normal contact force and tangential friction force model

Lankarani-Nikravesh (L-N) contact force model [26] has been confirmed by many experiments, and this model has been widely used [2729]. L-N contact force model performs well in dealing with contact problems with low energy dissipation and low collision velocities at high coefficient of restitution [30, 31]. Therefore, L-N contact force model was chosen in this paper to deal with normal contact force at collision point in spherical clearance joins of spatial parallel mechanism. Normal contact force model based on L-N contact force model can be written as

(19)

Where, KijKij is the stiffness coefficient, n is the power exponent of material and δij is pseudo-penetration depth, is pseudo-penetration velocity, it comes from the derivative of δij with respect to time, is initial value of pseudo-penetration velocity, ce is coefficient of restitution, υi and υj are Poisson's ratio of socket and ball, Ei and Ej are elastic modulus, R is radius of socket and r is radius of ball.

In equation (19), and in extended formula of Kij, δi and δj are fixed values related to material properties, and, .

In order to describe tangential force at collision point, the tangential friction model need to be established. Among the tangential force models for multi-body system dynamics, the Coulomb friction model is more common [32]. However, when the velocity approaches zero, this model does not describe the tangential friction at collision point well. Ambrosio [33] proposed an improved Coulomb friction model by introducing dynamic correction coefficients, improved Coulomb friction model compensates for the reduced stability of numerical integration due to the direction changed of friction force when collision velocity is zero and also the improved Coulomb friction model can well describe the dry friction and the viscous phenomena in spherical clearance joints at relatively low speeds. The improved Coulomb friction model can be expressed as

(20)

Where cf is sliding friction coefficient, cg is dynamic correction coefficient, Att is scalar of tangential velocity.

(21)

Symbols υa and υθ are threshold of given limiting velocity.

According to equations (19) and (20), the contact force ball exerts on socket can be expressed as

(22)

Reaction force of socket to ball is

(23)

Under the action of contact force, the moment at the centroid of components i and j can be written as

(24)

3 Dynamic model of SPCMM with spherical clearance joints

3.1 Kinematic model

Kinematic modeling of SPCMM is the basis for solving clearance dynamics, therefore, it is critical to build an appropriate kinematic model.

Position coordinate of origin of object coordinate system {07} in world coordinate system is OAo7 = [x y z]T.

Position coordinates of the spherical joints Si in world coordinate system are

(25)

Position coordinates of rotating joints Ri in world coordinate system are

(26)

Vectors of driven rods in world coordinate system are

(27)

Position coordinates of the rotating joints Ri in object coordinate system{o7} are

(28)

The transformational vector of the above vector changed from object coordinate system {o7} to world coordinate system {A} is

(29)

Vector εi is axis vector of rotation joints, it can be written as

(30)

According to the constraint relationship at the rotating joints Ri in the parallel mechanism, SiRi is always perpendicular to εi, constraint equation can be written as

(31)

By solving equation (31), the accompanying motion of the mechanism can be obtained. The Euler matrix of object coordinate system under world coordinate system can be written as

(32)

The Euler transformation process from point OA to point Ri can be written as

(33)

Simplified equation (33), and simplified results can be written as

(34)

Solved equation (34), the Euler angle γi and βi can be calculated as

(35)

(36)

Symbol Fi, Gi and Ni are simplification parameters.

The kinematics closed-loop vector equation of mechanical system is established as

(37)

According to the coordinate correspondence, the following equations can be obtained

(38)

where, symbol f(Pi) is driving function, by simplified equation (38), f(Pi) can be expressed as

(39)

The symbol A, B, C and D is simplification parameters.

The position coordinate OAo7 = [x y z]T of terminal platform is a function of variable x, β and γ. So linear velocity vp and angular velocity wp of terminal platform can be expressed by independent variable parameter, it can be expressed as

(40)

Matrix J1A is linear velocity transformation matrix.

(41)

The matrix J2A is angular velocity transformation matrix.

Velocity vector of rotating joints in world coordinate system is written as

(42)

Also, the acceleration vector of rotating joints in world coordinate system can be written as

(43)

3.2 Dynamic model with spherical clearance joints

The generalized coordinates of each component are selected to accurately describe the pose changes of the mechanism, it is expressed as

(44)

where ,Oi = (xi yi zi)T, φi = (αi βi γi)T(i = 1, 2, ..., 7), Oi is position coordinate of the origin of object coordinate system in world coordinate system, ϕi is Euler angle vector of object coordinate system in world coordinate system.

Spherical joints restrict three-dimensional movement and provide three-dimensional rotation, so its constraint equation is written as

(45)

where OAoi and OAoj are position coordinate of origin of object coordinate system {oi} and {oj} in world coordinate system{A}, oici is position coordinate about center point of socket in object coordinate system {oi}, ojcj is position coordinate about center point of ball in object coordinate system {oj},ARi  is Euler transformation matrix.

Sliding joints only provide translation in one direction and limits the pose in the other five directions, its constraint equation can be written as

(46)

where ARi is Euler matrix, Pi and Pj are center of two components that make up sliding joint, dxi, dyi, dzi are mutually perpendicular unit vectors fixed at dxj, dyj, dzj are mutually perpendicular unit vectors fixed at Pj, dxi and dxj are parallel to the translational direction of sliding joint. The first two equations restrict the movement of the sliding joint in the other two directions, and the last three equations restrict the rotation in three directions.

The rotating pair provides a rotational pose around the axis, so the movement in three directions and rotation in the other two directions should be restricted, and the constraint equation can be written as

(47)

where ARi is Euler matrices, ui and uj are rotational axis vectors of shaft and bearing, OAoi and OAoj are position coordinate of origin of object coordinate system in world coordinate system, ri and rj are center points of shaft and bearing

The rigid constraint equation of SPCMM including all joints can be written as

(48)

where , f(Pi) is driving function of translation joints.

Formula (48) derivation of time t, the velocity constraint equation is written as

(49)

where ΦQ is Jacobean matrix, is generalized velocity vector, Φt means that constraint equation Φ(Q) derivate of time t, σ  is the other side of the equation.

The equation (49) derivate of time t, acceleration constraint equation is written as

(50)

where represents generalized acceleration vector of mechanism which means derivate of time t, means derivate of the generalized coordinate Q, ΦQt means ΦQ derivate of time t, Φtt means Φt derivate of time t, γ is other side of the equation.

Mass matrix of SPCMM is followed

(51)

where Nk = diag[mk mk mk] (k = 1, 2...7), the rotational inertia Jk = diag[Ixk Iyk Izk] (k = 1, 2...7).

The above results are substituted into dynamic equation with Lagrange multiplier, equation may write as follow

(52)

where λ is Lagrange multiplier, it only related to the internal force and internal moment of connected rods, FQ is generalized force vector and it included external force and external torque.

When SPCMM carried out with no load, the generalized force vector FQ is expressed as

(53)

where Mk = [mkg 0 0 0 0 0]T (k = 1, 2...7), g is the acceleration of gravity.

Expression of formula (50) and (52) in form of differential algebra, the union matrix is shown as

(54)

Because displacement constraint and velocity constraint of multibody system are not used in equation (54), so there may be a violated constraint. Baumgarte [34] proposed a method to solve this problem using feedback control theory to enhance the stability of numerical calculation, the new union matrix is written as

(55)

where ac and βc are correction parameter bigger than zero, is that the rigid constraint equation Φ(Q) derivate of t.

In clearance dynamics model, the ball moves freely in the socket, the force constraint in spherical joint will replace the displacement constraint in rigid dynamics model, and the constraint equation and generalized force vector considering spherical clearance joints S1, S2 and S3 are as follows:

(56)

(57)

In formula (57),

where the force is component of contact force of spherical clearance joint in world coordinate system, is the moment of ball and socket from collision point to component centroid.

The union matrix with spherical clearance joints is written as

(58)

4 Dynamics behavior analysis of SPCMM

4.1 Solving process and applied parameters

Numerical technique is used to solve dynamic model of SPCMM with spherical clearance joints, and detailed structural parameters and simulation parameters is written in Tables 1 and 2. Simulation trajectory is given as following

(59)

According to Figure 4, the detailed numerical simulation process can be summarized as

  • Solve initial displacement and velocity of each component at time=0 by inverse kinematics, these data will be loaded into dynamics program with multiple spherical clearance joints.

  • Solve displacement constraint equations including all prismatic joints, spherical joints and rotating joints. The velocity and acceleration constraint equations can be obtained by taking time derivatives of the displacement constraint equations.

  • Establish kinematics model of spherical clearance joint, and solve normal contact force and tangential contact force at the impact point.

  • The force constraint in spherical clearance joint is introduced into Baumgarte stability equation instead of position constraint of component to enhance the stability of numerical simulation.

  • Judge the contact state at the contact point and calculate collision depth, if collision depth is greater than 0, calculate normal contact force and tangential contact force, and generalized displacement, velocity and acceleration can be got at the same time.

  • Perform iterative calculation until the end of motion time.

Table 1

design parameters of SPCMM.

Table 2

Clearance dynamics simulation parameters of SPCMM.

thumbnail Fig. 4

Flow chart of numerical simulation.

4.2 Dynamic response of SPCMM with multiple spherical clearance joints

4.2.1 Dynamic response considering clearance values

It is inevitable that clearance will cause an unpredictable impact and reduce the performance of mechanism. Concreting dynamic response and chaotic characteristics of parallel mechanism under different clearances can provide theoretical guidance for selecting the appropriate clearance in production and assembly, and set the clearance value as the single variable, set other variables to fixed values and the friction coefficient is fixed at 0.2.

From Figure 5a, the displacement curves of the terminal platform in X direction has no obvious deviation under different clearance values, and three curves coincided together. In the Figure 5b, curves have obvious fluctuations within the 0–0.3 s, and the amplitude of fluctuation increases while clearance value becomes bigger. Maximum deviation is 0.012 rad when clearance value is 0.7 mm, maximum deviation is 0.008 rad when clearance value is 0.5 mm, and maximum deviation is 0.006 rad when clearance value is 0.3 mm. After 0.3 s, the curves basically coincide and mechanism start a stable operation state. It is a conclusion change of clearance values have a great influence on the displacement response in Beta in the initial stage of movement, curve deviation becomes bigger when matched a bigger clearance value.

From Figure 6a, It is obvious clearance values caused a big different fluctuations within 0–0.4 s. Maximum fluctuation peak is 0.186 m/s while matched a fixed clearance value 0.7 mm, maximum fluctuation peak is 0.167 m/s while matched a fixed clearance value 0.5 mm, maximum fluctuation peak is 0.147 m/s while matched a fixed clearance value 0.3 mm. In Figure 6b, and the fluctuation is greater when matched a bigger clearance value, maximum fluctuation peak is 1.511 rad/s while matched a fixed clearance value 0.7 mm, maximum fluctuation peak is 1.104 rad /s while matched a fixed clearance value 0.5 mm, maximum fluctuation peak is 0.906 rad /s while matched a fixed clearance value 0.3 mm, but the fluctuations weaken very quickly, the mechanical system becomes stable and the curves basically coincide when it is time 0.4 s. It can be concluded that clearance has a great effect on the velocity in X and Beta of mechanism, when the impact phenomenon in joints disappear and the spherical joints elements are in contact, the mechanism can work in a stable motion state.

From Figure 7a, acceleration fluctuation becomes larger while matched a greater clearance value, the peak acceleration in X direction is 667.5 m/s2 when clearance value set a fixed value 0.7 mm, the peak acceleration in X direction is 446.2 m/s2 when clearance value set a fixed value 0.5 mm, the peak acceleration is 430.7 m/s2 when clearance value set a fixed value 0.3 mm. In Figure 7b, the acceleration in Beta is bigger than X, this means the clearance value does a great effect on acceleration in Beta, the initial fluctuation range of acceleration is very large, and then it weakened rapidly, until the time comes to 0.4 s the motion state becomes stable. The reason for this fluctuation rule of acceleration in X and Beta is that coupling motion and mutual interference occurs between the three spherical joints.

In Figure 8, contact forces fluctuate greatly within 0–0.21 s in different joints, violent collision occurred between joints elements at start of motion, and fluctuation peaks is different. Fluctuation peak is 809.7 N in S1 joint when clearance value matched a fixed value 0.3 mm, fluctuation peak is 463.3 N in S2 joint and fluctuation peak is 755.8 N in S3 joint, after 0.21 s, collision disappeared and mechanism are working steadily, here is a conclusion impact phenomenon in S1 is more obvious than impact phenomenon in S2 and S3. In Figure 9, the contact forces fluctuate greatly within 0–0.33 s in different joints S1, S2 and S3, after 0.33 s, collision disappeared and mechanism are working smoothly. When clearance value matched a fixed value 0.5 mm, fluctuation peak is 840 N in S1 joint, fluctuation peak is 1061 N in S2 joint and fluctuation peak is 982.7 N in S3 joint, it is concluded impact phenomenon in S2 is more obvious than impact phenomenon in S1 and S3. In Figure 10, contact forces fluctuate greatly within 0–0.39 s in different joints S1, S2 and S3, after 0.39 s, collision disappeared and mechanism are working smoothly. When clearance value matched a fixed value 0.7 mm, fluctuation peak is 1094 N in S1 joint, fluctuation peak is 917.4 N in S2 joint and fluctuation peak is 1176 N in S3 joint, it is concluded impact phenomenon in S3 is more obvious than impact phenomenon in S1 and S2. Compare contact forces considering different clearance values in joints S1, S2 and S3 and find that contact force increased while clearance values increased, it takes longer time for mechanism to reach stable motion state when matched a bigger clearance value.

thumbnail Fig. 5

Displacement response in X and Beta under different clearance values.

thumbnail Fig. 6

Velocity response in X and Beta under different clearance values.

thumbnail Fig. 7

Acceleration response in X and Beta under different clearance values.

thumbnail Fig. 8

Contact forces of S1, S2 and S3 with c is 0.3 mm.

thumbnail Fig. 9

Contact forces of S1, S2 and S3 with c is 0.5 mm.

thumbnail Fig. 10

Contact forces of S1, S2 and S3 with c is 0.7 mm.

4.2.2 Dynamic response considering friction coefficients

Oversized contact force existing in spherical clearance joint caused by violent collision may lead to an unbalanced motion state of the mechanism, and the main influence factor of the contact force is the friction coefficient, so it is particularly critical to research the dynamic response and chaotic characteristics of mechanism under different friction coefficients, here the clearance value takes a fixed 0.3 mm.

In Figure 11a, the displacement curves of different friction coefficients coincide, it can be concluded that different friction coefficients have little effect on the displacement curves in X direction. In Figure 11b, the curves fluctuate significantly before 0.3 s, and the fluctuation amplitude increases with increase of friction coefficient. Maximum deviation is 6.61 × 10−3 rad when friction coefficient is 0.15, maximum deviation is 6.77 × 10−3 rad when friction coefficient is 0.20, and maximum deviation is 6.90 × 10−3 rad when friction coefficient is 0.25. Fluctuation amplitude decreases gradually with movement of the mechanism, after 0.3 s, the curve does not fluctuate and the curves coincide. It can be concluded that friction coefficient has much effect on displacement curve at beginning of movement, the greater the friction coefficient, the greater the fluctuation of the curve.

From Figure 12a, the velocity curves in X fluctuate obviously within 0–0.3 s, and the peak of fluctuation reached 0.161 m/s, 0.156 m/s and 0.149 m/s separately while cf is 0.25, 0.20 and 0.15, and the curves fluctuate significantly before 0.3 s, the fluctuation amplitude decreases gradually within 0–0.3 s, and after 0.3 s, the curve has no obvious fluctuation. Here is a concluded that friction coefficient does a great effect on the velocity curve in X within 0–0.3 s, after 0.3 s, the friction coefficients have little effect on the velocity curve in X. In Figure 12b, the fluctuation of velocity curve in Beta is very obvious within 0–0.3 s, different friction coefficients lead to different fluctuation peaks, the fluctuation peak corresponding to friction coefficient 0.25 is 1.279 rad/s, the fluctuation peak corresponding to friction coefficient 0.20 is 1.277 rad/s, the fluctuation peak corresponding to friction coefficient 0.15 is 1.265 rad/s, after 0.3 s the fluctuation of curve is not obvious and curves basically coincide. It can be concluded that different friction coefficients have much effect on the velocity curve in Beta within 0–0.3 s, after 0.3 s, the friction coefficient has little effect on the curve, the mechanism reaches a stable motion state.

In Figure 13, the acceleration curves in X and Beta fluctuate violently within 0–0.3 s, after 0.3 s, the curves do not fluctuate and the curves present a stable periodic working state. In Figure 13a, biggest peak value is 533.2 m/s2 while friction coefficient matched a fixed value 0.15, peak value is 430.7 m/s2 while friction coefficient matched a fixed value 0.20, the peak value is 426 m/s2 while friction coefficient matched a fixed value 0.25, Here is a concluded that the greater the friction coefficient, the greater the fluctuation peak. In Figure 13b, the biggest peak value is 2193 m/s2 while friction coefficient is 0.20, the peak value is 1582 m/s2 while friction coefficient matched a fixed value 0.25 and peak value is 993.1 m/s2 while friction coefficient matched a fixed value 0.15. Here is a concluded that due to existence of coupled motion, influence of friction coefficient on acceleration in Beta is not following the rule the greater the friction coefficient, the greater the peak value of the curve, the acceleration curve in Beta has higher sensitivity than the acceleration curve in X when the friction coefficient changes.

From Figure 14, the contact force curves fluctuate greatly within 0–0.22 s, after 0.22 s, mechanism are working steadily. Peak value is 1057 N in joint S2, peak value is 966.4 N in joint S1 and the peak value is 866.8 N in joint S3, it is a conclusion collisions in S2 is more severe. From Figure 15, the contact force curves fluctuate greatly within 0–0.23 s, after 0.23 s, mechanism are working steadily. Peak value is 922.3 N in joint S2, peak value is 552.6 N in joint S1 and peak value is 444.7 N in joint S3, it is a conclusion collisions in S2 is more severe. In Figure 16, the contact force curves fluctuate greatly within 0–0.22 s, after 0.22 s, mechanism are working steadily. Peak value is 719.4 N in joint S1, peak value is 560.8 N in joint S2 and peak value is 500.2 N in joint S3, it is a conclusion collisions in S1 is more severe.

thumbnail Fig. 11

Displacement response in X and Beta under friction coefficients.

thumbnail Fig. 12

Velocity response in X and Beta under different friction coefficients.

thumbnail Fig. 13

Acceleration response in X and Beta under different friction coefficients.

thumbnail Fig. 14

Contact forces of S1 S2 and S3 with cf is 0.15

thumbnail Fig. 15

Contact forces of S1 S2 and S3 with cf is 0.20.

thumbnail Fig. 16

Contact forces of S1 S2 and S3 with cf is 0.25.

4.2.3 Reliability verification of dynamic model with multiple spherical clearance joints

To test the reliability of dynamic model for SPCMM with multiple spherical clearance joints, Matlab calculation results, Adams simulation results and ideal motion results without clearance joints were compared, the comparison results are shown in Figure 17.

Figure 17 shows simulation comparison curves between Adams and Matlab. Figure 17a and 17b shows the displacement curves in X and Beta, and Adams results, Matlab result and the no clearance results basically consistent. Figure 17c and 17d shows that velocity curves fluctuate greatly at the beginning of operation, the results solved by Matlab is bigger than the results solved by Adams, the curves are basically consistent when the mechanism reaches the stable state. Figure 17(e) and (f) shows that the acceleration curves fluctuate greatly at the beginning of operation, the peak fluctuation of Adams curve in X is 291 m/s2, the peak fluctuation of Matlab curve in X is 319.2 m/s2, the peak fluctuation of Adams curve in Beta is 1571 rad/s2, the peak fluctuation of Matlab curve in Beta is 1754 rad/s2, when the mechanism reaches the stable state, curves are basically consistent. It can be concluded that reliability of dynamics model of SPCMM with multiple spherical clearance joints is verified through comparison results.

thumbnail Fig. 17

Simulation comparison curves of Beta between Adams and Matlab.

4.3 Chaotic characteristics of SPCMM with multiple spherical clearance joints

Through Phase diagram, Poincare map and Bifurcation diagram of terminal platform in X and Beta direction, effects of different clearance values and friction coefficients on chaotic characteristics of mechanism can be researched and analyzed intuitively.

From Figure 18a, the Phase diagrams in X coincide basically and the images are smooth ellipses. It can be concluded different clearance values have little effect on Poincare map in X, the mechanism is in periodic motion in X analyzed by Poincare map. In Figure 18b, the picture is presented as ellipses of different widths, and the width of the ellipse becomes bigger when increases clearance value, it can be concluded that with increase of clearance value, influence of clearance on Poincare map in Beta becomes more obvious.

In Figure 19a, the Poincare maps in Beta corresponding to different clearance values are presented as three points at different positions, and the larger the clearance value, the more right the point is. It can be concluded that mechanism is in quasi periodic motion at this moment. From Figure 19b, the image is represented as three circular areas, and circular area becomes bigger when increases clearance value, it can be concluded that clearance value has much influence on the Poincare map in Beta, the mechanism is in quasi periodic motion state in Beta.

In Figure 20a, the curve has no divergence and presents as a stable straight line, it is resulted that mechanism is in periodic motion state all time while clearance value increases from 0.01 mm to 1 mm. In Figure 20b, the image gradually diverges while increases clearance value, and it can be concluded the clearance value does a great effect on the chaotic characteristics in Beta direction, and mechanism varied from periodic motion state to chaotic motion state with increase of clearance value.

In Figure 21a, the image appears as three smooth ellipses which coincide together, it can be concluded that friction coefficient has little effect on the phase diagram in X, and the mechanism is in periodic motion under different friction coefficients. From Figure 21b, there are three ellipses with wider contours which coincide together, the contours of the ellipse do not become wider when the friction coefficient becomes bigger, here is a concluded that friction coefficient has less effect on the phase diagram in Beta, and mechanism is in quasi periodic motion.

From Figures 22a and (b), the Poincare maps corresponding to different friction coefficients are located in the same position, the closer the mapping points are to the center, the denser the distribution is. It can be concluded that changing friction coefficient has little effect on Poincare mapping in X and Beta, and the mechanism is in quasi periodic motion.

From Figure 23, the Bifurcation diagrams of terminal platform in X and Beta present two convergent straight lines when friction coefficient increases from 0.01 to 0.4, it can be concluded that the mechanism is always in periodic working state, friction coefficient has little effect on the Bifurcation diagram of terminal platform in X and Beta.

thumbnail Fig. 18

Phase diagram of terminal platform.

thumbnail Fig. 19

Poincare map of terminal platform.

thumbnail Fig. 20

Bifurcation diagram of terminal platform.

thumbnail Fig. 21

Phase diagram of terminal platform.

thumbnail Fig. 22

Poincare map of terminal platform.

thumbnail Fig. 23

Bifurcation diagram of terminal platform.

5 Conclusions

As an indispensable resource for studying dynamic behavior of parallel mechanism, dynamic response and chaos characteristics of spatial parallel coordinate measuring mechanism with multi spherical clearance joints are studied.

  • A dynamic modeling method with multiple spherical clearance joints for spatial parallel coordinate measuring mechanism based on Lagrange multiplier technology is developed. The reliability of dynamics modeling method is verified by comparing the results from Adams and Matlab, also a detailed contact collision model of spherical clearance joints under different collision states are established.

  • Dynamic response of spatial parallel coordinate measuring mechanism is studied by comparing displacement, velocity, acceleration curves and contact force curves. The clearance values are respectively taken as fixed values of 0.3 mm, 0.5 mm and 0.7 mm and friction coefficient are respectively taken as fixed values of 0.15, 0.20 and 0.25. Comparing the analysis results, it is found that dynamic response of mechanism becomes more violent with increase of clearance value and friction coefficient.

  • Chaos characteristics of spatial parallel coordinate measuring mechanism are researched by comparing Phase diagram, Poincare map and bifurcation diagram. Comparing the analysis results, get a conclusion that clearance value has a great influence on chaotic characteristics of Beta direction, chaotic characteristics became more obvious when matched a bigger clearance value, and different friction coefficient has little effect on Chaos characteristics.

Funding

This research is supported by Shandong Provincial Natural Science Foundation (Grant no. ZR2022ME040).

Conflicts of interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. The authors have no conflict to disclose.

Data availability statement

All data has already been reported in the manuscript.

Author contribution statement

Xiulong Chen came up with the idea and Yongchao Sheng wrote the article.

References

  1. Y.H. Jia, X.L. Chen, Application of a new conformal contact force model to nonlinear dynamic behavior analysis of parallel robot with spherical clearance joints, Nonlinear Dyn. 108, 2161–2191 (2022) [CrossRef] [Google Scholar]
  2. X. Liang, Y. Takeda, An iterative method for the inverse kinematics of lower-mobility parallel mechanism with three RS or SR chains based on kinematically equivalent mechanism, Mech. Mach. Theory 141, 40–51 (2019) [CrossRef] [Google Scholar]
  3. M.X. Wang, Q.S. Chen, H.T. Liu, et al., Evaluation of the kinematic performance of a 3-RRS parallel mechanism, Robotica 39, 606–617 (2021) [CrossRef] [Google Scholar]
  4. Y.F. Zheng, H.P. Wan, J.Y. Zhang, et al., Local-coordinate representation for spatial revolute clearance joints based on a vector-form particle-element method, Int. J. Struct. Stabil. Dyn. 21, 2150093 (2021) [CrossRef] [Google Scholar]
  5. X.L. Chen, J.Y. Guo, Effects of spherical clearance joint on dynamics of redundant driving spatial parallel mechanism, Robotic 39, 1064–1080 (2021) [CrossRef] [Google Scholar]
  6. Y.L. Hou, Y.J. Deng, D.X. Zeng, et al., Dynamic modelling and properties analysis of 3RSR parallel mechanism considering spherical joint clearance and wear, J. Central South Univ. 28, 712–727 (2021) [CrossRef] [Google Scholar]
  7. C.Y. Rao, L.M. Xu, Q.H. Chen, et al., Dynamic modeling and performance evaluation of a 2UPR-PRU parallel kinematic machine based on screw theory, J. Mech. Sci. Technol. 35, 2369–2381 (2021) [CrossRef] [Google Scholar]
  8. S. Jiang, T. Wang, L.J. Xiao, Experiment research and dynamic behavior analysis of multi-link mechanism with wearing clearance joint, Nonlinear Dyn. 109, 1325–1340 (2022) [CrossRef] [Google Scholar]
  9. S. Fan, S.W. Fan, Difference between the ideal and combined spherical joints and its effects on parallel manipulators, Proc. Inst. Mech. Eng. C 234, 1112–1129 (2020) [CrossRef] [Google Scholar]
  10. C.U. Dogruer, B. Yildirim, Simulation of planar mechanism dynamics with index. I Formulation and projection method solution of differential algebraic equations, Arab. J. Sci. Eng. 47, 9305–9317 (2022) [Google Scholar]
  11. Z.S. Chen, L.M. Xu, W.Z. Zhang, et al., Closed-form dynamic modeling and performance analysis of an over-constrained 2PUR-PSR parallel manipulator with parasitic motions, Nonlinear Dyn. 96, 517–534 (2019) [CrossRef] [Google Scholar]
  12. B. Li, M.S. Wang, C.J. Gantes, et al., Modeling and simulation for wear prediction in planar mechanical systems with multiple clearance joints, Nonlinear Dyn. 108, 887–910 (2022) [CrossRef] [Google Scholar]
  13. Z.F. Bai, T.X. Liu, J.Y. Li, J. Zhao, Numerical and experimental study on dynamic characteristics of planar mechanism with mixed clearances, Mech. Based Des. Struct. Mach. 51, 6142–6165 (2022) [Google Scholar]
  14. Y. Chen, J. Feng, X. Peng, et al., An approach for dynamic analysis of planar multibody systems with revolute clearance joints, Eng. Comput. 37, 2159–2172 (2021) [CrossRef] [Google Scholar]
  15. X.Z. Wu, Y. Sun, Y. Wang, et al., Correlation dimension and bifurcation analysis for the planar slider-crank mechanism with multiple clearance joints, Multibody Syst. Dyn. 52, 95–116 (2021) [CrossRef] [MathSciNet] [Google Scholar]
  16. X.P. Wang, W.Z. Lin, X.M. Ji, et al., Dynamic analysis of a planar multibody system with multiple revolute clearance joints, Proc. Inst. Mech. Eng. C 233, 3429–3443 (2019) [CrossRef] [Google Scholar]
  17. Z.C. Song, S.R. Jiang, B. Chen, et al., Analysis of dynamic modeling and solution of 3-RPS parallel mechanism based on conformal geometric algebra, Meccanica 57, 1443–1455 (2022) [CrossRef] [MathSciNet] [Google Scholar]
  18. R.J. Zhang, J. Wu, C.D. Wu, et al., Two-step optimum design of a six-axis loading device with parallel mechanism considering dynamic coupling, Int. J. Robot. Autom. 37, 498–511 (2022) [Google Scholar]
  19. Y. Li, H. Zheng, B. Chen, et al., Dynamic modeling and analysis of 5PSS/UPU parallel mechanism with elastically active branched chains, Chinese J. Mech. Eng. 33, 33–44 (2020) [CrossRef] [Google Scholar]
  20. M. Zhang, C. Wang, C. Yan, et al., Design and dynamic analysis of a four-degree-of-freedom chaotic vibrating screen, Shock Vibr. 2021, 1–10 (2021) [Google Scholar]
  21. S. Erkaya, Experimental investigation of flexible connection and clearance joint effects on the vibration responses of mechanisms, Mech. Mach. Theory 121, 515–529 (2018) [CrossRef] [Google Scholar]
  22. X.L. Chen, J.Y. Guo, Effects of spherical clearance joint on dynamics of redundant driving spatial parallel mechanism, Robotica 39, 1064–1080 (2021) [CrossRef] [Google Scholar]
  23. G. Wang, L. Wang, Dynamics investigation of spatial parallel mechanism considering rod flexibility and spherical joint clearance, Mech. Mach. Theory 137, 83–107 (2019) [CrossRef] [Google Scholar]
  24. J.T. Zhang, Y.M. Song, D. Liang, Mathematical modeling and dynamic characteristic analysis of a novel parallel tracking mechanism for inter-satellite link antenna, Appl. Math. Model. 93, 618–643 (2021) [CrossRef] [MathSciNet] [Google Scholar]
  25. M. Chen, Q. Zhang, Y.F. Ge, et al., Dynamic analysis of an over-constrained parallel mechanism with the principle of virtual work, Math. Comput. Model. Dyn. Syst. 27, 347–372 (2021) [CrossRef] [MathSciNet] [Google Scholar]
  26. G.X. Wang, C.S. Liu, Further investigation on improved viscoelastic contact force model extended based on hertz's law in multibody system, Mech. Mach. Theory 153, 103986 (2020) [CrossRef] [Google Scholar]
  27. P. Flores, C.S. Koshy, H.M. Lankarani, et al., Numerical and experimental investigation on multibody systems with revolute clearance joints, Nonlinear Dyn. 65, 383–398 (2011) [CrossRef] [Google Scholar]
  28. C.S. Koshy, P. Flores, H.M. Lankarani, Study of the effect of contact force model on the dynamic response of mechanical systems with dry clearance joints: computational and experimental approaches, Nonlinear Dyn. 73, 325–338 (2013) [CrossRef] [Google Scholar]
  29. S.B. Farahan, M.R. Ghazavi, S. Rahmanian, Bifurcation in a planar four-bar mechanism with revolute clearance joint, Nonlinear Dyn. 87, 1–19 (2016) [Google Scholar]
  30. M. Machado, P. Moreira, P. Flores, et al., Compliant contact force models in multibody dynamics: evolution of the Hertz contact theory, Mech. Mach. Theory 53, 99–121 (2012) [CrossRef] [Google Scholar]
  31. A. Gummer, B. Sauer, Modeling planar slider-crank mechanisms with clearance joints in RecurDyn, Multibody Syst. Dyn. 31, 127–145 (2014) [CrossRef] [Google Scholar]
  32. Q. Tian, Y. Zhang, L. Chen, et al., Dynamics of spatial flexible multibody systems with clearance and lubricated spherical joints, Comput. Struct. 87, 913–929 (2009) [CrossRef] [Google Scholar]
  33. P. Flores, J. Ambrosio, On the contact detection for contact-impact analysis in multibody systems, Multibody Syst. Dyn. 24, 103–122 (2010) [CrossRef] [MathSciNet] [Google Scholar]
  34. P. Flores, M. Machado, E. Seabra, et al., A parametric study on the baumgarte stabilization method for forward dynamics of constrained mulibody systems, J. Comput. Nonlinear Dyn. 6, 73–82 (2011) [Google Scholar]
  35. Y. Wang, J. Angeles, et al., Kinematic analysis and optimum design of a novel 2PUR-2RPU parallel robot, Mech. Mach. Theory 139, 407–423 (2019) [CrossRef] [Google Scholar]

Cite this article as: Y. Sheng, X. Chen Dynamics behavior analysis of spatial parallel coordinate measuring mechanism with spherical clearance joints, Mechanics & Industry 25, 5 (2024)

All Tables

Table 1

design parameters of SPCMM.

Table 2

Clearance dynamics simulation parameters of SPCMM.

All Figures

thumbnail Fig. 1

Characteristic map of SPCMM.

In the text
thumbnail Fig. 2

The kinematic scheme of SPCMM.

In the text
thumbnail Fig. 3

Spherical clearance joint model.

In the text
thumbnail Fig. 4

Flow chart of numerical simulation.

In the text
thumbnail Fig. 5

Displacement response in X and Beta under different clearance values.

In the text
thumbnail Fig. 6

Velocity response in X and Beta under different clearance values.

In the text
thumbnail Fig. 7

Acceleration response in X and Beta under different clearance values.

In the text
thumbnail Fig. 8

Contact forces of S1, S2 and S3 with c is 0.3 mm.

In the text
thumbnail Fig. 9

Contact forces of S1, S2 and S3 with c is 0.5 mm.

In the text
thumbnail Fig. 10

Contact forces of S1, S2 and S3 with c is 0.7 mm.

In the text
thumbnail Fig. 11

Displacement response in X and Beta under friction coefficients.

In the text
thumbnail Fig. 12

Velocity response in X and Beta under different friction coefficients.

In the text
thumbnail Fig. 13

Acceleration response in X and Beta under different friction coefficients.

In the text
thumbnail Fig. 14

Contact forces of S1 S2 and S3 with cf is 0.15

In the text
thumbnail Fig. 15

Contact forces of S1 S2 and S3 with cf is 0.20.

In the text
thumbnail Fig. 16

Contact forces of S1 S2 and S3 with cf is 0.25.

In the text
thumbnail Fig. 17

Simulation comparison curves of Beta between Adams and Matlab.

In the text
thumbnail Fig. 18

Phase diagram of terminal platform.

In the text
thumbnail Fig. 19

Poincare map of terminal platform.

In the text
thumbnail Fig. 20

Bifurcation diagram of terminal platform.

In the text
thumbnail Fig. 21

Phase diagram of terminal platform.

In the text
thumbnail Fig. 22

Poincare map of terminal platform.

In the text
thumbnail Fig. 23

Bifurcation diagram of terminal platform.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.