Modeling and Analysis of a 2-DOF Spherical Parallel Manipulator

The kinematics of a two rotational degrees-of-freedom (DOF) spherical parallel manipulator (SPM) is developed based on the coordinate transformation approach and the cosine rule of a trihedral angle. The angular displacement, angular velocity, and angular acceleration between the actuators and end-effector are thus determined. Moreover, the dynamic model of the 2-DOF SPM is established by using the virtual work principle and the first-order influence coefficient matrix of the manipulator. Eventually, a typical motion plan and simulations are carried out, and the actuating torque needed for these motions are worked out by employing the derived inverse dynamic equations. In addition, an analysis of the mechanical characteristics of the parallel manipulator is made. This study lays a solid base for the control of the 2-DOF SPM, and also provides the possibility of using this kind of spherical manipulator as a 2-DOF orientation, angular velocity, or even torque sensor.


Introduction
Compared with traditional serial manipulators, the parallel manipulator is a more promising branch of robotic mechanism which has the distinct advantages of high stiffness, high precision, large load-weight ratio, high speed, and high acceleration. These advantages lead to its accordance with the trends of modern electrical and mechanical equipment. Therefore, the parallel manipulator has been a popular topic in manufacturing for decades [1][2][3][4]. With the development of robotics and automation engineering, spherical parallel mechanisms of various forms have attracted wide attention, and have found successful applications in industrial fields. In addition to the general characteristics of the parallel manipulator, the spherical parallel manipulator has the special advantages of simple structure, large workspace, interference-free of limbs, straightforward kinematics computation, reliable control, and so on. As early as 1931, Gwinnett proposed an amusement instrument based on a spherical parallel mechanism [5]. Cox developed a 3-RRR (R indicates the revolute joint) spherical 3-degrees-of-freedom (DOF) parallel mechanism [6] which is characterized by the free rotation of the mobile platform around the intersections of the revolute joint axes. Cox and Tesar proposed a 3-DOF parallel robotic shoulder module [7] which is considered as a typical spherical parallel manipulator with triangular platforms and three identical chains, each of which consisting of three revolute joints. The center of rotation of this mechanism is located at the intersection of the revolute joint axes. Asada and Granito proposed a spherical parallel manipulator with three collinear actuators [8] which improves the symmetry and convenience of the assembly of the mechanism. Gosselin and Angeles proposed a spherical 3-DOF parallel manipulator with three coplanar actuators [9] wherein the revolute joint axes of the base and mobile platforms intersect at one point. This special structure results in its preferable symmetry in The degrees of freedom of this manipulator can be verified according to the revised Kutzbach-Grubler formula [3] is the number of common constraints for the spherical mechanism, n is the number of parts, g is the number of joints, fi is the degrees of freedom of the ith joint, v is the number of redundant constraints, and ζ is the isolated degree of freedom. In this research, n = 6, and the number of joints is g = 7. For the universal joint, 1  system. In branch 1, the axis X1j (j = 1, 2, 3) of the local coordinate system constantly points to the center of the revolute joint from origin O, and Y11 coincides with Y21, Y21, and 1 1  T B are in the same direction, Y13 is in the same direction with the Y axis of the global coordinate system, and Z1j (j = 1, 2, 3) can be determined by the right hand rule. In branch 2, the Y2j (j = 1, 2, 3) of the local coordinate system constantly points to the center of the revolute joint from the origin O, X21 coincides with X11, X22 and 2 2  B T are in the same direction, X23 is in the same direction with the X axis of the global coordinate system, and Z1j (j = 1, 2, 3) can be determined by the right hand rule. O-X11Y11Z11 coincides with O-X21Y21Z21. The position vectors of the six revolute joints in their local coordinate systems are presented, respectively, as follows: The degrees of freedom of this manipulator can be verified according to the revised is the number of common constraints for the spherical mechanism, n is the number of parts, g is the number of joints, f i is the degrees of freedom of the ith joint, v is the number of redundant constraints, and ζ is the isolated degree of freedom. In this research, n = 6, and the number of joints is g = 7. For the universal joint, f 1 = 2, for the other six revolute joints, f i = 1. There is no redundant constraint and isolated degree of freedom, so v = 0, ζ = 0 in this research. Thus, the degrees of freedom of this mechanism are M = 2.
In order to describe this 2-DOF SPM, a global coordinate system O-XYZ and the local coordinate systems O-X ij Y ij Z ij were established, as shown in Figure 2. The subscripts i and j represent the jth revolute joint of the ith branch, respectively, where i = 1, 2, and j = 1, 2, 3. The coordinate conventions are as follows. The center point O of the universal joint is the origin of each coordinate system. The projection of X axes on the base are same direction coincident with vector − − DB 1 , and those of the Y axes are with − − DB 2 ; thus, Z axes can be determined according to the right hand rule. P 1 , T 1 , B 1 , P 2 , T 2 and B 2 represent the position vectors of center point of each revolute joint in the global coordinate system. In branch 1, the axis X 1j (j = 1, 2, 3) of the local coordinate system constantly points to the center of the revolute joint from origin O, and Y 11 coincides with Y 21 , Y 21 , and T 1 × B 1 are in the same direction, Y 13 is in the same direction with the Y axis of the global coordinate system, and Z 1j (j = 1, 2, 3) can be determined by the right hand rule. In branch 2, the Y 2j (j = 1, 2, 3) of the local coordinate system constantly points to the center of the revolute joint from the origin O, X 21 coincides with X 11 , X 22 and B 2 × T 2 are in the same direction, X 23 is in the same direction with the X axis of the global coordinate system, and Z 1j (j = 1, 2, 3) can be determined by the right hand rule. O-X 11 Y 11 Z 11 coincides with O-X 21 Y 21 Z 21 . The position vectors of the six revolute joints in their local coordinate systems are presented, respectively, as follows: Figure 2. Schematic diagram of the coordinate system.

Representation of the Joints
Due to the mechanical constraint of the universal joint, the mobile platform practically has only two degrees of freedom, including roll angle x  and pitch angle y  . Based on the knowledge of robotics [3,21], the coordinate transformation matrices of the mobile platform in Euler angle form are: 1 0 0 ( ) 0 cos sin 0 sin cos where x  and y  are roll and pitch angles of the mobile platform, respectively.
In order to work out the kinematics of the 2-DOF SPM, this paper conducts a straightforward deduction. On one hand, once the mobile platform is assumed to be in a given orientation, the representation of the position of related points such as P1, P2, T1 and T2 can be obtained. On the other hand, the angular displacement of actuators will also lead to the representation of T1 and T2. So, the kinematic equations can be formulated according to the equivalence of T1 and T2. The details are as follows.
Assume that the original state of the 2-DOF SPM is that the Since the lower links are directly connected to the actuators, their position vectors can be represented with the input angular displacement θ1, θ2.

Representation of the Joints
Due to the mechanical constraint of the universal joint, the mobile platform practically has only two degrees of freedom, including roll angle φ x and pitch angle φ y . Based on the knowledge of robotics [3,21], the coordinate transformation matrices of the mobile platform in Euler angle form are: where φ x and φ y are roll and pitch angles of the mobile platform, respectively. In order to work out the kinematics of the 2-DOF SPM, this paper conducts a straightforward deduction. On one hand, once the mobile platform is assumed to be in a given orientation, the representation of the position of related points such as P 1 , P 2 , T 1 and T 2 can be obtained. On the other hand, the angular displacement of actuators will also lead to the representation of T 1 and T 2 . So, the kinematic equations can be formulated according to the equivalence of T 1 and T 2 . The details are as follows.
Assume that the original state of the 2-DOF SPM is that the local coordinate system O-X 11 Y 11 Z 11 coincides with the global coordinate system O-XYZ. For the first operation, let the local coordinate rotate around the Y axis by φ y to get a new coordinate system. Then, for the second operation, let the new coordinate rotate around the X axis by φ x to get the final coordinate system. So, the resulting vectors are: Since the lower links are directly connected to the actuators, their position vectors can be represented with the input angular displacement θ 1 , θ 2 .
Besides, the actuator locations can be described in global coordinates as, where α 1 , β are structural parameters as shown in Figure 1. R is the radius of the sphere of the 2-DOF SPM.

Inverse Kinematics
The inverse kinematics of a manipulator calculate the motion of actuators given the orientation of the mobile platform, which is the theoretical base for the control of the manipulator.
In branch 1, the angle between P 1 and B 1 is denoted as α 31 , thus The planes related to angles α 31 , α 2 , and α 1 intersect at point O. The angle between the planes of α 31 and α 1 is denoted as θ 11 . According to cosine rule of a trihedral angle, the geometric Equation (11) is obtained.
The angle of the plane of α 31 and O-XZ is set to θ 12 . According to Equation (4), P 1 and B 1 are in the plane of O-XZ. Thus θ 12 = 0 holds. The input angle of the actuator in branch 1 is θ 1 = θ 11 + θ 12 = θ 11 .

Forward Kinematics
The forward kinematics of this 2-DOF SPM calculate the orientation of the mobile platform, provided the angular displacement of the two actuators. This procedure is of great importance for the orientation sensor application of this SPM. More generally, the semi-closed loop control of the parallel mechanism also needs its forward kinematics so that the orientation can be obtained by reading out the feedback values of the actuators and a further computation.
The angle between vectors T 1 and P 1 is constantly equal to α 2 , which can be expressed by By solving Equation (15), the following solution of φ y can be obtained where Similarly, the angle between T 2 and P 2 is constantly equal to α 2 ; thus, one can obtain By calculating Equation (17), the following parameter can be obtained From the solutions to the forward kinematics mentioned above, one can infer that the rotation of the mobile platform around the Y axis is only dependent on the motion of actuator 1. Compared with this conclusion, the rotation of the mobile platform around X axis is dependent on the motion of both actuators.

Analysis of Velocity and Acceleration
The velocity and acceleration analysis plays a key role in the mechanical design, controller design, and hardware configuration.
be the orientation, rotation velocity, and rotation acceleration vectors of the mobile platform, respectively.
be the angular displacement, angular velocity, and angular acceleration vectors of the actuators, respectively. This section will deal with the input-output relation of the velocity and acceleration of the 2-DOF SPM.
where A and B are matrices as shown in Equation (19). From this result, one can likewise notice that the angular acceleration of the mobile platform relates to the input angular velocity, input angular acceleration, and the orientation and structural parameters of the manipulator.

The Dynamics Model Based on the Virtual Work Principle
The investigation of dynamics is dealing with the quantitative relation between the motion of the mobile platform and the input torques provided by the actuating joints. Furthermore, calculating the force or torque of each actuator on the active joint for the desired trajectory of the mobile platform is called inverse dynamics. In this research, since the mobile platform has thirty times greater mass than that of the four linkages, the mass of the four linkages is thus negligible for simplicity.
The inertial matrix tenor of the mobile platform is where I xx = I yy = 1 /4mr 2 , m is the mass of the mobile platform, and r is the radius of the mobile platform. It is to be noted that for this 2-DOF SPM prototype, the center of gravity lies in the point above the center of the universal joint by a distance of h, as shown in Figure 3. Therefore, the position vector of the center of gravity of the mobile platform is where φ x and φ y are roll and pitch angles of the mobile platform, respectively. The differentiation of pcog with respect to time results in: cos cos sin sin cos 0 The differentiation of p cog with respect to time results in: Equation (25) can be rewritten in the following form: where [G v ] is the first-order translation influence coefficient matrix on the motion of the center of gravity of the mobile platform. The translation here especially indicates the translation of the center of gravity of the mobile platform. By differentiating Equation (26) with respect to time, one can obtain: Based on the virtual work principle, the dynamic equation where g = (0, 0, −g 0 ) T is the gravitational acceleration vector. G φ and [G v ] are the two influence matrices of the input on the rotation and motion of the center of gravity of the mobile platform, respectively. m and [I c ] are the mass and inertial tensor of the mobile platform, respectively. a is the intermediate motion parameter as defined in Equation (27). This dynamic equation describes the relation of input torque and the resulting motion of the mobile platform.

Numerical Simulation
To validate the theoretical results derived above, the motion planning, inverse, forward kinematic simulation, and inverse dynamics simulation were conducted. The parameters of the 2-DOF SPM for this simulation are listed in Table 1. For the first typical trajectory motion following simulation, a linear orientation variation for 2 s is designed in which the orientation of the mobile platform can be described by Equation (29) during time 0 ≤ t ≤ 2 s. The angular displacement of the mobile platform is shown in Figure 4. The resultant rotational motion of the mobile platform displays strict linearity. By using the inverse kinematics developed in Section 3.2, the angular displacement curves of the two actuators can be worked out, as shown in Figure 5. One can observe that the actuator displays parabolic curves as a result of the quadratic rotational motion function defined in Equation (29). Then, the angular velocity and angular acceleration of the mobile platform can be derived as shown in Figures 6 and 7, respectively. The angular velocity and angular acceleration of the two actuators can be obtained as shown in Figures 8 and 9, respectively. The differential relation between the velocity and acceleration and also the symmetry in two rotational degrees of freedom can be observed. Based on the dynamic equation shown in Equation (28), the driving torque curves of the actuators are shown in Figures 10 and 11. In each of them, the dashed line (T i_1g ) and solid line (T i_0g ) represent the driving torque with and without gravity, respectively, from which one can notice that the linear angular acceleration of the mobile platform leads to its relatively small constant inertial force supplied by the actuators, while the gravity of the mobile platform accounts for a large proportion of the driving torque in this type of uniformly accelerated motion.
For the first typical trajectory motion following simulation, a linear orientation variation for 2 s is designed in which the orientation of the mobile platform can be described by Equation (29) during time 0 2 t   s.
The angular displacement of the mobile platform is shown in Figure 4. The resultant rotational motion of the mobile platform displays strict linearity. By using the inverse kinematics developed in Section 3.2, the angular displacement curves of the two actuators can be worked out, as shown in Figure 5. One can observe that the actuator displays parabolic curves as a result of the quadratic rotational motion function defined in Equation (29). Then, the angular velocity and angular acceleration of the mobile platform can be derived as shown in Figures 6 and 7, respectively. The angular velocity and angular acceleration of the two actuators can be obtained as shown in Figures 8  and 9, respectively. The differential relation between the velocity and acceleration and also the symmetry in two rotational degrees of freedom can be observed. Based on the dynamic equation shown in Equation (28), the driving torque curves of the actuators are shown in Figures 10 and 11. In each of them, the dashed line (Ti_1g) and solid line (Ti_0g) represent the driving torque with and without gravity, respectively, from which one can notice that the linear angular acceleration of the mobile platform leads to its relatively small constant inertial force supplied by the actuators, while the gravity of the mobile platform accounts for a large proportion of the driving torque in this type of uniformly accelerated motion.   is designed in which the orientation of the mobile platform can be described by Equation (29) during time 0 2 t   s.
The angular displacement of the mobile platform is shown in Figure 4. The resultant rotational motion of the mobile platform displays strict linearity. By using the inverse kinematics developed in Section 3.2, the angular displacement curves of the two actuators can be worked out, as shown in Figure 5. One can observe that the actuator displays parabolic curves as a result of the quadratic rotational motion function defined in Equation (29). Then, the angular velocity and angular acceleration of the mobile platform can be derived as shown in Figures 6 and 7, respectively. The angular velocity and angular acceleration of the two actuators can be obtained as shown in Figures 8  and 9, respectively. The differential relation between the velocity and acceleration and also the symmetry in two rotational degrees of freedom can be observed. Based on the dynamic equation shown in Equation (28), the driving torque curves of the actuators are shown in Figures 10 and 11. In each of them, the dashed line (Ti_1g) and solid line (Ti_0g) represent the driving torque with and without gravity, respectively, from which one can notice that the linear angular acceleration of the mobile platform leads to its relatively small constant inertial force supplied by the actuators, while the gravity of the mobile platform accounts for a large proportion of the driving torque in this type of uniformly accelerated motion.                     In addition, for the typical sinusoidal trajectory motion following simulation, the orientation of the mobile platform can be described by Equation (30) during time 0 2 t   s.
By using the inverse kinematics developed in Section 3.2, the angular displacement curves of the two actuators are shown in Figure 12. One can observe that the actuator displays periodical curves slightly different from canonical sinusoidal or cosine curves. Moreover, Figure 13 illustrates the comparison between the desired orientation curve and the forward kinematic result derived from curves in Figure 12. It is clear that the inverse and forward kinematics obtained in Section 3.3 agree with each other, which validates the forward kinematics deduced in this research. As this structure contains two serial chains, the forward kinematics has the analytical form instead of iterative solutions usually appearing in parallel mechanisms of all types. It is efficient to calculate the orientation of the mobile platform once to get the angular information of the active joints. This fact is one of the distinct characteristics of the possibility of using this 2-DOF SPM as a 2-DOF orientation sensor.  In addition, for the typical sinusoidal trajectory motion following simulation, the orientation of the mobile platform can be described by Equation (30) during time 0 ≤ t ≤ 2s.
By using the inverse kinematics developed in Section 3.2, the angular displacement curves of the two actuators are shown in Figure 12. One can observe that the actuator displays periodical curves slightly different from canonical sinusoidal or cosine curves. Moreover, Figure 13 illustrates the comparison between the desired orientation curve and the forward kinematic result derived from curves in Figure 12. It is clear that the inverse and forward kinematics obtained in Section 3.3 agree with each other, which validates the forward kinematics deduced in this research. As this structure contains two serial chains, the forward kinematics has the analytical form instead of iterative solutions usually appearing in parallel mechanisms of all types. It is efficient to calculate the orientation of the mobile platform once to get the angular information of the active joints. This fact is one of the distinct characteristics of the possibility of using this 2-DOF SPM as a 2-DOF orientation sensor.  In addition, for the typical sinusoidal trajectory motion following simulation, the orientation of the mobile platform can be described by Equation (30) during time 0 2 t   s.
By using the inverse kinematics developed in Section 3.2, the angular displacement curves of the two actuators are shown in Figure 12. One can observe that the actuator displays periodical curves slightly different from canonical sinusoidal or cosine curves. Moreover, Figure 13 illustrates the comparison between the desired orientation curve and the forward kinematic result derived from curves in Figure 12. It is clear that the inverse and forward kinematics obtained in Section 3.3 agree with each other, which validates the forward kinematics deduced in this research. As this structure contains two serial chains, the forward kinematics has the analytical form instead of iterative solutions usually appearing in parallel mechanisms of all types. It is efficient to calculate the orientation of the mobile platform once to get the angular information of the active joints. This fact is one of the distinct characteristics of the possibility of using this 2-DOF SPM as a 2-DOF orientation sensor.  As for the derivatives of the angular displacement, the angular velocity and angular acceleration of the two actuators can be derived as shown in Figures 14 and 15, respectively. Because of the sharp variation of the angular displacement of actuator 2 at the moment of 0 t  , its angular velocity is quite considerable, as shown in Figure 14. Especially for the angular acceleration, the peak value gets As for the derivatives of the angular displacement, the angular velocity and angular acceleration of the two actuators can be derived as shown in Figures 14 and 15, respectively. Because of the sharp variation of the angular displacement of actuator 2 at the moment of t = 0, its angular velocity is quite considerable, as shown in Figure 14. Especially for the angular acceleration, the peak value gets as much as 248 • /s 2 . This kind of rigid impulse is very harmful for heavy load applications. Therefore, in a practical motion of the SPM, this sharp jump of acceleration should be avoided or decreased by trapezoidal or high-order polynomial velocity motion planning. As for the derivatives of the angular displacement, the angular velocity and angular acceleration of the two actuators can be derived as shown in Figures 14 and 15, respectively. Because of the sharp variation of the angular displacement of actuator 2 at the moment of 0 t  , its angular velocity is quite considerable, as shown in Figure 14. Especially for the angular acceleration, the peak value gets as much as 248 °/s 2 . This kind of rigid impulse is very harmful for heavy load applications. Therefore, in a practical motion of the SPM, this sharp jump of acceleration should be avoided or decreased by trapezoidal or high-order polynomial velocity motion planning.  Similarly, based on the dynamic equation expressed in Equation (28), the driving torque curve of the actuators in this sinusoidal trajectory tracking simulation are shown in Figures 16 and 17, in which Ti_1g and Ti_0g represent the driving torque with and without gravity, respectively. From these figures, one can observe that the maximum torque of both actuators is less than 2.5 Nm, which satisfies the provided performance of general commercial middle and small inertia AC servo motors. As for the comparison of the driving torque under conditions with and without gravity, the  As for the derivatives of the angular displacement, the angular velocity and angular acceleration of the two actuators can be derived as shown in Figures 14 and 15, respectively. Because of the sharp variation of the angular displacement of actuator 2 at the moment of 0 t  , its angular velocity is quite considerable, as shown in Figure 14. Especially for the angular acceleration, the peak value gets as much as 248 °/s 2 . This kind of rigid impulse is very harmful for heavy load applications. Therefore, in a practical motion of the SPM, this sharp jump of acceleration should be avoided or decreased by trapezoidal or high-order polynomial velocity motion planning.  Similarly, based on the dynamic equation expressed in Equation (28), the driving torque curve of the actuators in this sinusoidal trajectory tracking simulation are shown in Figures 16 and 17, in which Ti_1g and Ti_0g represent the driving torque with and without gravity, respectively. From these figures, one can observe that the maximum torque of both actuators is less than 2.5 Nm, which satisfies the provided performance of general commercial middle and small inertia AC servo motors. As for the comparison of the driving torque under conditions with and without gravity, the Similarly, based on the dynamic equation expressed in Equation (28), the driving torque curve of the actuators in this sinusoidal trajectory tracking simulation are shown in Figures 16 and 17, in which T i_1g and T i_0g represent the driving torque with and without gravity, respectively. From these figures, one can observe that the maximum torque of both actuators is less than 2.5 Nm, which satisfies the provided performance of general commercial middle and small inertia AC servo motors. As for the comparison of the driving torque under conditions with and without gravity, the conclusion for the linear acceleration case does not apply in this simulation due to the fact that the accelerations are highly time-varying and in larger magnitude.    To investigate the singularity of the 2-DOF SPM, its first-order influence coefficient matrix is considered as the Jacobian matrix J in this research. The mapping of condition number when the orientation of the mobile platform is an arbitrary value within the workspace is shown in Figure 18. Since the condition number lies in the rational finite interval for the orientation of the mobile platform, one can reach the conclusion that this mechanism is singularity-free in the workspace. Additionally, the condition number varies with the combination of the two rotational angles rather than keeping constant, which implies that this mechanism fails to achieve homogenous performance, especially for the error and static generalized force transmission from the active joints to the mobile platform. This is also the reason why the driving force of the actuator displays a high value for certain orientations. As a tentative prototype, the 2-DOF SPM as shown in Figure 19 was fabricated, the diameter of the mobile platform was 600 mm, the height of the vertical column was 378.6 mm, and the total mass was 73.67 kg. Two YASKAWA MSMA AC servo motors (SGMGV-09A) were employed as the actuators. The basic motion was achieved with the Googol CPAC motion controller as the master To investigate the singularity of the 2-DOF SPM, its first-order influence coefficient matrix is considered as the Jacobian matrix J in this research. The mapping of condition number when the orientation of the mobile platform is an arbitrary value within the workspace is shown in Figure 18. Since the condition number lies in the rational finite interval for the orientation of the mobile platform, one can reach the conclusion that this mechanism is singularity-free in the workspace. Additionally, the condition number varies with the combination of the two rotational angles rather than keeping constant, which implies that this mechanism fails to achieve homogenous performance, especially for the error and static generalized force transmission from the active joints to the mobile platform. This is also the reason why the driving force of the actuator displays a high value for certain orientations.  To investigate the singularity of the 2-DOF SPM, its first-order influence coefficient matrix is considered as the Jacobian matrix J in this research. The mapping of condition number when the orientation of the mobile platform is an arbitrary value within the workspace is shown in Figure 18. Since the condition number lies in the rational finite interval for the orientation of the mobile platform, one can reach the conclusion that this mechanism is singularity-free in the workspace. Additionally, the condition number varies with the combination of the two rotational angles rather than keeping constant, which implies that this mechanism fails to achieve homogenous performance, especially for the error and static generalized force transmission from the active joints to the mobile platform. This is also the reason why the driving force of the actuator displays a high value for certain orientations. As a tentative prototype, the 2-DOF SPM as shown in Figure 19 was fabricated, the diameter of the mobile platform was 600 mm, the height of the vertical column was 378.6 mm, and the total mass was 73.67 kg. Two YASKAWA MSMA AC servo motors (SGMGV-09A) were employed as the actuators. The basic motion was achieved with the Googol CPAC motion controller as the master As a tentative prototype, the 2-DOF SPM as shown in Figure 19 was fabricated, the diameter of the mobile platform was 600 mm, the height of the vertical column was 378.6 mm, and the total mass was 73.67 kg. Two YASKAWA MSMA AC servo motors (SGMGV-09A) were employed as the actuators. The basic motion was achieved with the Googol CPAC motion controller as the master control unit. An elementary motion control containing two degrees of freedom was firstly carried out, validating the structural design and kinematic model. As future work, the thorough verification of the precision, workspace, and dynamic performance should be conducted. control unit. An elementary motion control containing two degrees of freedom was firstly carried out, validating the structural design and kinematic model. As future work, the thorough verification of the precision, workspace, and dynamic performance should be conducted. Figure 19. Field photo of the 2-DOF SPM prototype.

Conclusions
The kinematics of a 2-DOF SPM was elaborated based on the coordinate transformation approach and cosine rule of a trihedral angle. Moreover, the angular displacement, angular velocity, and angular acceleration relation between the actuators and the end-effector were determined. The dynamic model of the parallel manipulator were established by employing the virtual work principle

Conclusions
The kinematics of a 2-DOF SPM was elaborated based on the coordinate transformation approach and cosine rule of a trihedral angle. Moreover, the angular displacement, angular velocity, and angular acceleration relation between the actuators and the end-effector were determined. The dynamic model of the parallel manipulator were established by employing the virtual work principle and the first-order influence coefficient matrix of the manipulator. Finally, typical numerical simulations of the 2-DOF SPM reveal the quantitative relation between the input and output in terms of angular displacement, angular velocity, and acceleration. This preliminary result lays the foundation for better planning motion and control of the 2-DOF SPM, and also provides the fundamental verification of the possibility of using this kind of spherical mechanism as 2-DOF orientation or torque sensors.