1. Introduction
The creative design of mechanisms is widely regarded as one of the most important branches of the theory of mechanisms [Reference Yang, Ding and Kecskeméthy1]. As a new type of mechanism, parallel driving mechanism (PDM) [Reference Shi, Zhang, Wang, Li and Huang2] represents a subset of hybrid manipulator [Reference Li, Ye, Shen, Wang and Geng3–Reference Gupta, Gupta and Ekta5]. By definition, the PDM consists of an operating mechanism, multiple driving mechanisms, and connecting joints. Based on this, this paper proposed a 2T2R PDM. The 2T2R manipulator organically combines spatial translation and rotation, and is widely used in industry. For example, this manipulator can be combined with linear guide rail to form a five-axis linkage CNC machine tool [Reference Zhou, Li, Li, Meng, Li and Xu6, Reference Tang and Zhang7] and can also be used as motion simulator [Reference Wei, Fang and Liu8], vibration isolation platform [Reference Lu, Wu, Ding and Chen9], and so on. Among which, CNC machine tools have requirements for high accuracy of the manipulator, while the others have certain requirements for high-frequency and high-acceleration characteristics. In general, the PDM involves more sub-closed loops [Reference To and Webb10, Reference Nguyen, Zhou, Kang and Le11] and over-constraints [Reference Chen, Zhang, Qin and Sun12, Reference Xu, Liu, Yao and Zhao13]. Excluding actuators, its operating mechanism is independent of the driving mechanism, which makes the mechanism perform well in terms of structural stiffness [Reference Zhao, Gao, Li, Liu and Zhao14, Reference Wu, Li, Wang and Wang15] and high acceleration [Reference Shan, Li, Liu and Huang16, Reference Li-xin and Yong-gang17]. Similar to parallel mechanisms, the PDM has drawbacks such as small workspace, difficulty in solving forward kinematics, and complex control. In addition, there are new challenges in dynamic modeling due to over-constraints and the closed-loop structure.
In contrast to the open-chain serial manipulators, the dynamic analysis of parallel manipulators is inherently complex due to their closed-loop structure. Nevertheless, dynamic modeling plays an essential role in structural design and control scheme design for parallel manipulators. However, the requirements for dynamic models are different, which depends on the exact application. For the control scheme design, the dynamic model must be efficient enough for the real-time calculation; for the structural design, it is necessary to calculate the driving forces/torque and the constrained forces/moments simultaneously [Reference Chen, Xu, Zhang and Li18]. In general, dynamic modeling methods [Reference Chen, Wu, Deng and Wang19–Reference Tian, Flores and Lankarani21] can be divided into four categories: Lagrangian formula [Reference Flores, Ambrósio, Claro, Lankarani and Koshy22, Reference McGrath, Howard and Baker23], Newton–Euler method [Reference Chen, Xu, Zhang and Li18, Reference Marques, Roupa, Silva, Flores and Lankarani24, Reference Dasgupta and Mruthyunjaya25], virtual work principle [Reference Chen, Zhang, Qin and Sun26, Reference Zhao and Gao27], and Kane equation [Reference Pedrammehr, Nahavandi and Abdi28]. Among them, the Lagrangian formula is advantageous in solving the driving force of manipulator, with its explicit dynamic form conducive to dynamic analysis and control. Newton–Euler method can be used not only to solve the driving force but also to solve all the constraints/torques of joints effectively. However, the number of equations is generally smaller than the unknown number under the over-constrained mechanism, which makes it challenging to solve the dynamic model. At present, there have been some studies conducted on the dynamic modeling of the over-constrained parallel manipulators. Xu et al. [Reference Xu, Liu, Yao and Zhao29] proposed a new approach to the force analysis of over-constrained parallel mechanisms with lower mobility parallel mechanism, with consideration given to the spatial compound elastic deformation of the rod. Chen et al. [Reference Chen, Xu, Zhang and Li18] used the Newton–Euler formula and the natural orthogonal complement method to establish dynamic models with and without constraint force/moment simultaneously, with the deformation caused by over-constraint force/moment deduced. Bi et al. [Reference Bi and Kang30] constructed the dynamic model of the over-constrained mechanism Exechon by combining Newton–Euler formula with compliance conditions. The established dynamic model can be combined with the control system to improve the precision and dynamic performance of real-time control. For a mechanism subjected to the active and passive redundant constraints without parasitic motion, Chen et al. [Reference Chen, Yu, Li and Wang31] established an efficient dynamic model in line with the principle of virtual work by introducing classical optimization criteria and deformation coordination conditions. To sum up, the above-mentioned studies have contributed great ideas to solving the over-constrained mechanisms.
As can be seen from the above, introducing link flexibility is a common method for analyzing over-constrained mechanisms. The links’ elastic deformation can not only supplement the deformation coordination equation [Reference Elgolli, Houidi, Mlika and Romdhane32, Reference Yi, Li and Wang33] in the over-constrained direction but also cause pose error of the end-effector. The former provides a basis for the lightweight design. It focuses on using the links’ elastic deformation to solve the constraint force/moment [Reference Wang, Liu and Huang34] and joint friction [Reference Sharifzadeh, Arian, Salimi, Masouleh and Kalhor35] of the manipulator. As for the latter, it focuses on the impact of the rigid–flexible coupling [Reference Chen, Jiang and Wang36–Reference Xiao, Yan, Chen, Zhang and Jiang38] characteristics, large-range rigid body motion, and the elastic deformation of the rod, on the errors in position and orientation of the end-effector. However, the hybrid manipulator is free from large-range rigid body motion in the over-constrained direction, and the link flexibility conforms to the assumption of linear elasticity and small deformation. Therefore, it is proposed in this paper to supplement the deformation compatibility equation in the over-constrained direction by introducing the flexibility matrix and the superposition principle into the process of theoretical calculation. To facilitate simulation, the assumed mode method is used for the link flexibility in ABAQUS. Then, dynamic verification is carried out using the rigid–flexible coupling dynamic model of the co-simulation of ADAMS and ABAQUS.
To sum up, there have been some methods developed to solve the over-constraint issue. However, the PDM is still more challenging than the traditional parallel mechanism in terms of structural complexity and static indeterminate times. In addition, sub-closed loops are another characteristic of PDM, and the over-constrained mechanisms containing sub-closed loops [Reference To and Webb39, Reference Nguyen, Zhou, Kang and Le40] are usually coupled in a nonlinear manner, which makes it difficult to conduct kinematic and dynamic analysis of them [Reference Chen, Wu, Deng and Wang19, Reference Li and Staicu41, Reference Shao, Wang, Yang and Liu42]. At present, there are still few reports on the dynamic modeling of the PDMs. In this paper, with a new 4-DOF PDM as the research object, the over-constrained direction is first analyzed through the screw theory. Then, the flexibility matrix and superposition principle are introduced to supplement the deformation compatibility equation in the over-constrained direction. Finally, a complete dynamic model involving all active forces/moments and restraint forces/moments is constructed. This modeling method is also applicable to other over-constrained PDMs.
The remainder of this paper is organized as follows. In Section 2, the system conceptual design is described. Then, the kinematics modeling of the PDM is detailed in Section 3, where the positions, velocities, and accelerations of all bodies in the manipulator are determined. In Section 4, dynamics modeling is performed using the Newton–Euler method. Moreover, compatibility equations are introduced in Section 5 to solve the over-constraint forces/moments. Then, numerical simulations and experiments are conducted in Section 6. Lastly, the conclusions of this study are drawn in Section 7.
2. System conceptual design
2.1. Mechanism configuration design
In this paper, a new 4-DOF PDM with two rotational DOFs and two translational DOFs is taken as the research object. The virtual prototype for the mechanism is represented in Fig. 1. $\boldsymbol{\$ }_{T1}$ and $\boldsymbol{\$ }_{T2}$ denote the translation along x-axis and y-axis, respectively. $\boldsymbol{\$ }_{R3}$ and $\boldsymbol{\$ }_{R4}$ represent rotation around z-axis and x-axis, respectively. The 4-DOF PDM is composed of a fixed platform, a moving platform (MP), an operating mechanism and driving mechanisms. The driving mechanisms consist of four distributed limbs, with two identical RPR limbs and two identical SPU limbs. Each limb is formed by a cylinder and a piston and the prismatic joints of limbs are actuated, which are the inputs of the system, the position and orientation of the MP is the output. The operating mechanism consists of a six-bar mechanism and a revolute joint in series. As shown in Fig. 2, $R_{1i},R_{2i},R_{3i},R_{4i}$ , and $R_{4}(i=1,2)$ represent the revolute joint; $U_{2j}(j=3,4)$ represents the universal joint; $S_{2j}$ denotes the spherical joint; $P_{1i}$ and $P_{2j}$ represent the prismatic joint. The reference frame P -xyz and the moving frame A -xyz are attached to the MP and the base, respectively. $l_{1i}$ and $l_{2j}$ denote the location of the prismatic joint in each limb. The geometric parameters are given in Appendix A.
2.2. Workspace analysis
The boundary of the reachable workspace is determined with the help of boundary search method as shown in Fig. 3(a), where the z-axis is parasitic motion. Figure 3(b) shows its projection in the x–y direction. The link extensions, $l_{11},l_{12},l_{23}$ , and $l_{24}$ , are limited between −50 (mm) and 50 (mm) since the left and right grounded prismatic links are assumed to be physically constrained by the ground at those boundary angles. In addition, other link angles are also subject to certain constraints based on physical properties.
2.3. Over-constraint analysis
The PDM perform well in structural stiffness due to its closed-loop structure. Meanwhile, the closed-loop structure brings more over-constraints, which is necessary for the mobility and force analysis. Next, the over-constrained analysis of the 4-DOF PDM will be derived in detail. The joint type and number are shown in Table I.
The number of over-constraints on the 4-DOF PDM is calculated as follows:
3. Kinematic analysis
Kinematic analysis is an indispensable step in the evaluation and application. Therefore, position, velocity, and acceleration analysis of the 4-DOF PDM are presented in this section. In the following, the kinematics of the hybrid manipulator is solved in modules. As shown in Fig. 1, the module-1 represents the operating mechanism, and the module-2 denotes the driving mechanisms, including the 1st, 2nd RPR limbs and the 3rd, 4th SPU limbs. Meanwhile, the module-1 is solved prior to the module-2.
3.1. Inverse position analysis
3.1.1. Inverse position of module-1
Inverse position analysis determines the positions of each limb when the position and orientation of the MP are given. Any vector in the frame P -xyz can map to the frame A -xyz.
in which
where T and P denote the orientation matrix and the position vector of the MP, respectively.
For the 1st and 2nd driving limbs, the position vector $R_{3i}(i=1,2)$ with respect to the frame A -xyz can be obtained based on the homogeneous transformation. In addition, the position vector $R_{2i}$ with respect to the frame A -xyz can be expressed by $\theta _{\mathrm{1}i}$ . We can derive the closed-loop motion equation as
Eq. (5) can be simplified as follows
in which
where atan2 stands for four-quadrant inverse tangent function. $a_{1i},b_{1i}$ , and $d_{1i}$ represent the simplified values of Eq. (5), and their expressions are detailed in Appendix B.
3.1.2. Inverse position of module-2
According to the constraint relationship, the length of $l_{1i}$ can be obtained
Similar to Eq. (8), the position vector of $U_{2j}(j=3,4)$ with respect to the frame O -xyz can be obtained. We can derive the length of $l_{2j}$ as
3.2. Velocity and acceleration analysis
The displacement of the prismatic joint is given by Eqs. (8) and (9). However, to derive the inverse dynamic formulation for the 4-DOF PDM with Newton–Euler method, the positions, velocities, and accelerations of all bodies in the hybrid manipulator are required to be derived. In the following, we will derive the inverse kinematics in detail.
3.2.1. Velocity and acceleration of module-1
The angular velocity, acceleration, and linear acceleration of the MP can be obtained
in which
where $\boldsymbol{\omega }_{mp}$ and $\boldsymbol{\varepsilon }_{mp}$ represent the angular velocity and acceleration of the MP, respectively. $\boldsymbol{a}_{mp}^{n}$ and $\boldsymbol{a}_{mp}^{\tau }$ denote the normal and tangential acceleration of the MP, respectively.
For the operating mechanism, $\theta _{2i}$ can be written as follows:
Furthermore, differentiating Eq. (12) with respect to time yields
Since the operating mechanism is a planar mechanism, the angular velocity and acceleration of $l_{3i}$ can be written as follows:
Similarly, the acceleration of the operating mechanism can be written as
in which
3.2.2. Velocity and acceleration of module-2
The Euler angle of the 1st and 2nd driving limbs can be written as follows:
The angular velocity and angular acceleration of $l_{i}$ can be obtained by differentiating Eq. (17). Similar to Eq. (14), the angular velocity and acceleration vector can be expressed as
Furthermore, the acceleration at the center of mass of the cylinder can be written as follows:
in which
Similarly, the acceleration at the center of mass of the piston can be written as
in which
where $\boldsymbol{a}_{li2}^{k}$ denote the Coriolis acceleration of the piston in 1st and 2nd driving limbs.
For the U joint in the 3rd and 4th driving limbs, many papers assume that the angular velocity and acceleration are perpendicular to the direction along SPU limb for the convenience of derivation. In this paper, the actual angular velocity and acceleration values are obtained through geometric constraint relationship. The local coordinate frame attached to the U joint is shown in Fig. 4.
The coordinate systems of the 3rd and 4th driving limbs are shown in Fig. 4. First, we assume that the orientation matrix of the 3rd and 4th driving limbs frame Q -xyz with respect to the reference frame A -xyz can be derived in terms of three rotational angles $\alpha,\beta$ , and $\gamma$ satisfying the z-x-z convention:
where s and c represent the abbreviations of sine and cosine, respectively. In addition, $\boldsymbol{Q}_{j}$ can be obtained by the following equation:
Hence, taking Eqs. (23) and (24) into account, the Euler angle of $l_{2j}$ can be written as follows:
Differentiating Eq. (25) with respect to time, $\dot{\alpha }, \ddot{\alpha }, \dot{\beta }$ , and $\ddot{\beta }$ can be obtained, respectively. To derive the Euler angle $\gamma$ , it is necessary to supplement equations based on the geometric constraint relationship. Hence, we can derive $\boldsymbol{U}_{j2}$ as
where vector $\boldsymbol{U}_{j1}$ denotes the axis along with the MP at the U joint, vector $\boldsymbol{U}_{j2}$ denotes the axis along with the 3rd and 4th SPU limb. From Eqs. (23) to (26), $\gamma$ can be expressed as
Differentiating Eq. (27) with respect to time, $\dot{\gamma }, \ddot{\gamma }$ can be obtained, respectively. In addition, $\boldsymbol{U}_{j3}$ can be obtained by vector cross product, which is the direction of the constrained torque at the U joint.
Furthermore, to obtain the angular velocity of each link, we need to map the Euler angles rate to the angular velocity. The mapping principle is shown in Fig. 5.
According to the above transformation, Eq. (29) can be obtained
By sorting out Eq. (29), the following equation can be obtained
in which
Furthermore, differentiating Eq. (30) with respect to time, the following equation can be obtained
For the 1st and 2nd driving limbs, L 1i represents the vector between the two joints $R_{2i}$ and $R_{4i}$
where $\boldsymbol{r}_{2i}$ and $\boldsymbol{r}_{4i}$ represent the vectors of joint $R_{2i}$ and $R_{4i}$ with respect to frame { A }, respectively. Moreover, use Q 1i to denote the unit vector along L 1i
The velocity of any point in the moving frame can be obtained by vector cross-product
where $\boldsymbol{\omega }$ and V represent the angular velocity and the linear velocity of the MP, respectively.
Differentiating Eq. (7) with respect to time yields:
where $\dot{\boldsymbol{P}}$ denotes the velocity matrix of the MP.
Similarly, linear velocity of l i by differentiating Eq. (8) can be obtained
The linear acceleration of the MP can be obtained by the following formula
where $\boldsymbol{\varepsilon }$ denotes the angular acceleration of the MP.
Furthermore, differentiating Eq. (36) with respect to time yields:
where $\ddot{\boldsymbol{P}}$ denotes the acceleration matrix of the MP.
Similarly, linear acceleration of $l_{i}$ can be obtained as time derivative of Eq. (37)
For the 3rd and 4th driving limbs, the linear velocity can be calculated by Eq. (35)
Differentiating Eq. (42) with respect to time, the following formula can be obtained
Hence, the velocity of the $l_{2j}$ can be expressed as follows:
in which
According to the Eq. (38), the following formula can be obtained
Furthermore, linear acceleration of $l_{2j}$ can be obtained by differentiating Eq. (44)
Hence, the acceleration of the $l_{2j}$ can be written as follows:
So far, all the kinematic parameters required for the inverse dynamics are determined. In the next section, attention will be paid on the deduction of dynamic equations of the 4-DOF PDM.
4. Dynamic analysis
For the structural design, the driving forces and the constrained forces/moments should be calculated simultaneously to evaluate the force condition of the system. Therefore, this section presents the dynamics modeling of the 4-DOF PDM. First, the force analysis of individual bodies, limbs are derived with Newton–Euler method. Then, based on the flexibility matrix and compatible deformation are employed to supplementary equations. Similar to the kinematic analysis, the dynamic derivation of the 4-DOF PDM is also solved in modules.
4.1. Forces and moments analysis for the module-1
The forces and moments analysis diagrams of link $r_{21}$ are given in Fig. 6. The reaction forces and moments at fixed points $R_{1i}$ are denoted by $\boldsymbol{f}_{e1}$ and $\boldsymbol{M}_{e1}$ . The forces and moments at moving points $R_{2i}$ are denoted by $\boldsymbol{f}_{32}$ and $\boldsymbol{M}_{32}$ , respectively.
The Newton–Euler equations for link $r_{21}$ are given as
Similar to link $r_{21}$ , other links’ forces and moments of the operating mechanism can be obtained.
Furthermore, the forces and moments analysis for the MP are shown in Fig. 7.
The Newton–Euler equations for MP are given as
4.2. Forces and moments analysis of the module-2
Assume that each driving limb consists of two parts, the cylinder and the piston, where the velocities and accelerations of their centers of masses are determined in Section 3. Assume that the centers of masses of the cylinder and the piston are located at a distance of $c_{i1}$ and $c_{i3}$ , and their masses are denoted by $m_{i1}$ and $m_{i2}$ , respectively. Moreover, consider that the cylinder and the pistons are symmetric about their axes, and their centers of masses lie at their midlengths.
The forces and moments analysis diagrams of 1st and 2nd driving limb are given in Fig. 8. The reaction forces and moments at fixed points $R_{4i}$ are denoted by $\boldsymbol{f}_{fi}$ and $\boldsymbol{M}_{fi}$ . The forces and moments at moving points $R_{2i}$ are denoted by $\boldsymbol{f}_{ci}$ and $\boldsymbol{M}_{ci}$ , and the internal forces and moments between cylinders and pistons are denoted by $\boldsymbol{f}_{mi}$ and $\boldsymbol{M}_{mi}$ , respectively.
The Newton–Euler equation for 1st and 2nd driving limbs are given as
where ${}^{A}{\boldsymbol{I}}{_{ci_{1}}^{}}$ denotes the inertia matrix of the cylinder evaluated at a frame parallel to { A } located at point $c_{i1}$ . Since the cylinder is moving with respect to this frame, this matrix can be evaluated by rotation rule from the inertia matrix easily obtained about { A i } attached to the cylinder at point $c_{i1}$ . This can be done through the following transformation:
Note the formula of vector mixed product is as follows:
Multiply $\boldsymbol{t}_{i}$ on both sides of Eq. (51), and the following equation can be obtained
Forces and moments analysis for piston are shown in Fig. 9.
The Newton–Euler equations for piston are given as
By multiplying $\boldsymbol{s}_{i}$ on both sides of Eq. (55), formulas that only include driving force and kinematic parameters can be derived as follows:
The forces and moments analysis diagrams of 3rd and 4th driving limbs are given in Fig. 10. The reaction forces and moments at fixed points $S_{2j}$ are denoted by $\boldsymbol{f}_{sj}$ and $\boldsymbol{M}_{sj}$ . The forces and moments at moving points $U_{2j}$ are denoted by $\boldsymbol{f}_{uj}$ and $\boldsymbol{M}_{uj}$ , and the internal forces and moments between cylinders and pistons are denoted by $\boldsymbol{f}_{mj}$ and $\boldsymbol{M}_{mj}$ , respectively.
The Newton–Euler equation for 3rd and 4th driving limbs are given as
According to Eq. (56), the following equation can be obtained
Vector cross-product formula is as follows:
Furthermore, by multiplying $\boldsymbol{s}$ on both sides of Eq. (57) and combing Eq. (59), the following equation can be obtained
The forces and moments analysis of piston are shown in Fig. 11.
The Newton–Euler equations for piston segments are given as
Multiplying $\boldsymbol{s}$ on both sides of Eq. (61) and substituting Eq. (60) into Eq. (61), we have
To sum up, there are 84 equations based on the Newton–Euler method, but 93 unknown variables need to be solved due to the over-constraints of the 4-DOF PDM. In the next section, nine compatibility equations are introduced considering the flexibility matrix and compatible deformation.
5. Deformation compatibility analysis
Based on the redundant constraint analysis in Section 2.2, we know that the 4-DOF PDM possesses nine passive redundant constraints, although the four actuator driving forces can be obtained with the classical method, the internal force analysis of this manipulator is still a statically indeterminate problem, and therefore the joint reaction forces cannot be uniquely determined. In this section, to obtain the unique solution, the over-constrained direction is analyzed by the screw theory, and the flexibility of the limbs is considered and the deformation compatibility condition is introduced.
5.1. Link’s elastic deformation
In this paper, the link is regarded as a homogeneous rod. As shown in Fig. 12, when the end link is subjected to external forces $\boldsymbol{F}\boldsymbol{=}[\begin{array}{cccccc} F_{x} & F_{y} & F_{z} & M_{x} & M_{y} & M_{z} \end{array}]^{T}$ , corresponding six-dimensional deformation will be formed, including tension, torsion, and bending.
The relation between the external force F and the corresponding small deflection X along the applied force axis can be written as follows
in which
where C represents the flexibility matrix with respect to frame $\{\boldsymbol{O}\}$ . $I_{y},I_{z}$ , and $I_{p}$ denote cross-section inertia parameters of the rod, respectively. $l$ and $A$ are rod length and cross-sectional area of rod. E and G represent elastic modulus and shear modulus of materials, respectively.
For the two-link mechanism depicted in Fig. 13, according to the force balance condition of the rod, the Jacobian matrix J as the projection of forces $\boldsymbol{F}_{i}$ to external forces $\boldsymbol{F}$ are given
in which
where ${}_{ok}^{oi}{\boldsymbol{R}}{}$ represents rotational transformation matrix between frame $\{\boldsymbol{O}_{\boldsymbol{k}}\}$ and frame $\{\boldsymbol{O}_{\boldsymbol{i}}\}$ . ${}^{oi}{\boldsymbol{P}}{_{ok}^{}}$ is the position vector with respect to the frame $\{\boldsymbol{O}_{\boldsymbol{k}}\}$ . $\boldsymbol{S}({}^{oi}{\boldsymbol{P}}{_{ok}^{}})$ denotes antisymmetric matrix of position vector ${}^{oi}{\boldsymbol{P}}{_{ok}^{}}$ .
Referring to Eq. (63), the relation between the force $\boldsymbol{F}_{i}$ and the corresponding small deflection $\boldsymbol{X}_{i}$ along the applied force axis can be written as follows
Hence, the deformation of end link $\Delta \boldsymbol{X}_{i}$ caused by the deformation of the rod i is as follows
Therefore, the deformation of the end link under the external force F can be derived by the sum of linear deformation of all rod i
5.2. Rigid-flexible coupling simulation
As indicated above, the 4-DOF PDM possesses nine passive redundant constraints. Since the SPU limb own six degrees-of-freedom, the over-constraints are caused by other links besides the SPU limbs. In this section, based on the screw theory, the over-constrained direction is analyzed. Moreover, the over-constraints were verified via the co-simulation with ABAQUS and ADAMS software. First, the flexible links are meshed by hexahedron elements in ABAQUS, and the mnf files of the meshed links are obtained via ABAQUS. Then, loading the link files, rigid links are replaced with flexible bodies in ADAMS. The simulation model is shown as Fig. 14.
From the above analysis, we know that over-constraints are caused by operating mechanism and two RPR limbs. In order to supplement the deformation coordination equations, we need to further determine the direction of over-constraints. Based on the screw theory, the constraint screw of the limb is as follows:
As shown in Fig. 15, the operating mechanism and two RPR limbs possess the constrained force along z -axis and constrained couples around x -axis, y -axis, which means each link in the operating mechanism and two RPR limb own three over-constrained motions.
5.3. Deformation compatibility equation
In this section, on the condition of the numbers and direction of over-constraints, we derive the deformation compatibility equation in detail. First, the operating mechanism and two RPR limbs are divided into three modules. Then derive the deformation of each module according to the link’s elastic deformation in Section 5.1. Finally, the deformation compatibility equations are derived at joints $R_{21},R_{22}$ , and $R_{32}$ .
The elastic deformation of limbs 1 and 2 at joint $R_{32}$ as shown in Fig. 16.
Referring to Eq. (65), we have
where ${}^{B1}{\boldsymbol{F}}{_{r32,r4}^{}}$ is the force applied to link $r_{4}$ by link $r_{32}$ with respect to the frame { A }. $\boldsymbol{J}_{B2,A}$ represents the transformation matrix between frame { B 2} and frame { A }.
Based on the Eq. (69), the elastic deformations of the links $r_{21}$ and $r_{31}$ at the joint $R_{31}$ with respect to the frame { B 3} are
Hence, the elastic deformation of the links $r_{21},r_{31}$ , and $r_{4}$ at the joint $R_{32}$ with respect to the frame { A } are
Similarly, the elastic deformation of limb2 at $R_{32}$ is given as
The deformation caused by limbs 3 and 4 at joint $R_{32}$ as shown in Fig. 17.
Referring to Eq. (69), the elastic deformation of limb 3 at joint $R_{21}$ is given as
Similarly, the elastic deformation of limb 4 at joint $R_{21}$ is given as
The deformation caused by limbs 5 and 6 at joint $R_{32}$ as shown in Fig. 18.
Based on the Eq. (69), the elastic deformation of limb 5 at joint $R_{22}$ are given as
Similarly, the elastic deformation of limb 6 at joint $R_{22}$ are given as
The deformations of the operating mechanism and two RPR limbs have to be compatible with each other to satisfy the geometric constraints. Hence, translational deformation along z and rotational deformation about x and y axis at joints $R_{21},R_{22},R_{32}$ should be equal, namely
So far, the complete dynamics models with constrained forces/moments of the over-constrained 4-DOF PDM with sub-closed chains are established. The actuation forces and constrained forces/moments can be obtained simultaneously, which is essential for structure design. Also, clearance and friction can be easily integrated into the system equation.
6. Simulation and experiment
We set the trajectory of the end-effector of the manipulator with respect to the inertial reference frame P -xyz as follows:
where $S_{x}$ and $S_{y}$ denote the displacement along the x-axis and y-axis with respect to the fixed reference coordinate system, $S_{\alpha }$ and $S_{\gamma }$ denote the angle along the x-axis and z-axis with respect to the fixed reference coordinate system, respectively.
6.1. Simulation validation
According to the derivation of the above formula, the displacements, velocities, accelerations, driving forces, and all constraint forces/moments of the 4-DOF PDM are obtained. The numerical simulation results are qualitatively consistent with the theoretical results. As shown in Fig. 19, $l_{1,a}, v_{1,a}$ , and $a_{1,a}$ (i = 1, 2, 3, 4) denote the theoretical displacements, velocities, and accelerations of the slider respectively. $l_{i,n}, v_{i,n}$ , and $a_{i,n}$ (i = 1, 2, 3, 4) denote the numerical displacements, velocities, and accelerations of the slider, respectively.
The results of co-simulation with ADAMS and ABAQUS are shown in Fig. 14. Figure 20 shows the active forces computed by the proposed method and ADAMS; Figure 21 shows the active forces computed by the proposed method and Lagrangian method; Figure 22(a) and (b) show the constraint forces/moments in the over-constrained direction at joints $R_{32}$ and $R_{21}$ , respectively. $f_{i\_ z}^{a},m_{i\_ x}^{a}$ and $m_{i\_ y}^{a}(i=1,2)$ represent the theoretical value and $f_{i\_ z}^{n},m_{i\_ x}^{n}$ and $m_{i\_ y}^{n}(i=1,2)$ represent the simulation value.
It is assumed that the cylinders and pistons are symmetrically distributed in this paper. That is, the influence of the eccentricity of the electric cylinder on the constraint torque at the U-joint is not considered. In order to verify the influence of acceleration on the constraint torque and driving force at the U-joint, the constraint torques at the U-joint are calculated when the maximum translational acceleration along y-axis is up to 50, 500, and 1000 mm/s2. Moreover, the driving forces are solved with and without constraint torque at the U-joint, and the driving forces error percentage are plotted as shown in Fig. 23(b).
6.2. Experiment validation
According to the CAD model displayed in Fig. 1, a physical prototype has been built for the motion control experiment. As shown in Fig. 24, actuators are driven by the electric piston-cylinder, which are connected to an embedded controller (CX5130) made by the company Beckhoff through an Ethercat bus. Based on the position loop mode PID controller to achieve the target trajectory, it is one of the commonly used controllers in academic research and industrial products. The motion control experiment based on PID controller is shown in Figure 25. The measured values of displacements, velocities, and accelerations of the four driving limbs are shown in Fig. 26, which are consistent with the theoretical calculation.
In this paper, the control cycle of the embedded controller (CX5130) is 2 ms. As shown in Fig. 26, the measured position and velocity of the four driving limbs are consistent with the theoretical value, but the lag of phase is 20 ms. Due to the short control cycle, the position signal obtained by the encoder undergoes two differential and filtering processes, resulting in significant signal distortion. As shown in Fig. 26, the measured acceleration not only have phase lag but also have attenuation in amplitude.
7. Conclusion
In this paper, a 4-DOF over-constrained PDM with sub-closed chains is proposed and its prototype is developed. By considering the flexibility of the operating mechanism and two RPR limbs, the deformation compatibility condition is introduced to solve a statically indeterminate problem. The actuation forces and constrained forces/moments are obtained simultaneously, which is conducive to the design of the PDM.
In the symmetrical structure of the electric piston-cylinder, the value of constraint moments at the U-joint under different accelerations is analyzed by taking into account its impact on driving forces. According to the results of comparative experiment conducted on the end-effector under different accelerations in a certain direction, the impact on the active force is about 1%, regardless of the constraint moment at the U-joint. The effectiveness of the kinematic model is demonstrated by ADAMS and the experiment based on classical PID control, while that of the dynamic model is confirmed by the co-simulation with MATLAB, ADAMS and ABAQUS. In addition, the actuation forces are also verified by using the Lagrangian method. The present study provides a valuable reference for future study on structural optimization. In future studies, the rigid–flexible coupling dynamic model will also be established to obtain rigid–flexible coupling dynamic characteristics and an accurate constraint forces/moments.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Author contribution
Yangyang Huang, Jinzhu Zhang, and Xiaoyan Xiong conceived and designed the study. Shengxiang Liu conducted data gathering. Yangyang Huang wrote the article.
Financial support
The authors would like to acknowledge (1) Project (51905367) supported by The National Nature Science Foundation of China, (2) Project (201901D211011) supported by The Foundation of Applied Basic Research General Youth Program of Shanxi.
Competing interests
The authors declare no competing interests exist.
Ethical approval
Not applicable.
Appendix A
The kinematic and dynamic parameters of the PDM virtual prototype used in the simulation are as follows (default in mm unit):
Appendix B