Kinematics and singularity analysis of a novel hybrid industrial manipulator

Abstract This paper proposes a new type of hybrid manipulator that can be of extensive use in industries where translational motion is required while maintaining an arbitrary end-effector orientation. It consists of two serially connected parallel mechanisms, each having three degrees of freedom, of which the upper platform performs a pure translational motion with respect to the mid-platform. Closed-form forward and inverse kinematic analysis of the proposed manipulator has been carried out. It is followed by the determination of all of its singular configurations. The theoretical results have been verified numerically, and the 3D modeling and simulation of the manipulator have also been performed. A simple optimal design is presented based on optimizing the kinematic manipulability, which further demonstrates the potential of the proposed hybrid manipulator.


Introduction
Serial and parallel robot manipulators have had significant attention from academia.Serial manipulators are based on open kinematic chains, with the links being connected serially via different types of joints.They have been used in various industrial applications and elsewhere where there is a demand for large workspaces.On the other hand, parallel manipulators are based on closed-loop kinematic chains with a standard topological structure of two platforms (fixed base and moving) connected via multiple legs, where each leg is yet another serial chain.Owing to their parallel structure, they offer the advantages of higher stiffness, better operating speeds, and higher payload.However, they suffer from reduced as well as more complex overall workspace.A significant amount of research has been done on the design and kinematic analysis [1][2][3][4][5], as well as on the determination of workspace and singularities [1,3,[5][6][7][8][9][10][11][12][13] of parallel manipulators.Based on different applications, numerous lower mobility parallel manipulators have been designed [1-5, 12, 13], in which the moving platform either purely translates [1][2][3]12], or purely rotates [4,5,13], or else performs a combined translational and rotational motion [14].A purely translational 3-DOF parallel manipulator named 3-UPU was first proposed by Tsai [2].Di Gregorio et al. [1] further studied the geometric and assembly conditions necessary for the translational nature of a more general 3-RRPRR parallel manipulator in detail, with the 3-UPU being a special case of 3-RRPRR.The 3-UPU manipulator has three identical legs, each with two universal joints and a prismatic pair in U-P-U order.The 3-UPU wrist, having a pure spherical rotation of the platform with respect to the base, was proposed by Karouia et al. [4], which was further analyzed by Di Gregorio [5,13].
A hybrid manipulator combines the advantages of serial and parallel manipulators in terms of having higher rigidity and a significantly larger workspace than a parallel manipulator, thereby offering a way to perform complex tasks.Due to this, there has been a growing interest in hybrid manipulators over the past few decades, even though the available literature is currently limited compared to serial and parallel manipulators.A hybrid manipulator is generally formed by combining serial and parallel structures.Quite a few serial-parallel-type hybrid manipulators have been proposed in the literature, which are made up of parallel mechanisms connected in series.However, there are other type of hybrid manipulators as well consisting of multiple serial chains linked to a single parallel "base" (commonly seen in robotics).Shahinpoor [15] presented the kinematic analysis of a novel hybrid manipulator made by serially connecting two parallel manipulators.The inverse kinematics was completed by obtaining a set of nonlinear equations that must be solved numerically.Singh et al. [16] developed the complete kinematic model and determined the jacobian matrix of a 7-DOF spatial hybrid robotic arm for use in assisting a surgeon during medical surgery.Verification of the model was performed via a physical prototype.Maier and Woernle [17] presented a generalized inverse kinematic analysis of a cable-suspension hybrid manipulator comprising a 3-DOF cable-driven parallel mechanism and a 3-DOF serial mechanism.However, these hybrid manipulators require numerical solutions.A 6-DOF spatial manipulator with a novel architecture, designated as H6A (i.e., Hybrid 6-Axis), is discussed in ref. [18].The manipulator consists of two arm-like branches attached to a rigid waist at one end and coupled through a wrist at the other end.
Romdhane [19] presented the kinematic analysis of a hybrid manipulator, again made up of two serially connected 3-DOF parallel mechanisms.The three successive prismatic joints in the lower mechanism restricted the motion of the mid-platform only to three translations.At the same time, the top platform performed spherical rotation about the mid-platform due to the existence of a spherical joint between the two platforms.Eight inverse kinematic solutions in closed form were obtained for the upper platform.Tanio Tanev [20] proposed a novel 6-DOF hybrid manipulator made up of two novel tripod mechanisms involving prismatic, revolute, and spherical joints.Closed-form forward and inverse kinematics were presented.Tanev obtained four forward kinematic solutions each for both the tripod mechanisms.Zheng et al. [21] utilized the two special 3-UPU mechanisms, which are the pure translational 3-UPU [2] and the pure rotational 3-UPU [4], and combined them serially such that the translational 3-UPU becomes the lower mechanism and the rotational 3-UPU becomes the upper mechanism.The forward position analysis of this new hybrid manipulator was carried out using quaternions, and two and eight solutions for the lower and upper mechanisms were obtained, respectively.
Gallardo-Alvarado et al. [22,23] presented a novel topology for the 2(3-RPS) hybrid manipulator and approached its kinematic and dynamic analysis using screw theory.The simplicity and compactness of the expressions derived in the study were noted.Hu [24] presented a general approach to derive the unified forward and inverse Jacobian matrices for hybrid manipulators composed of lower mobility parallel manipulators connected serially.Huang et al. [25] proposed using a 3-DOF parallel mechanism as the main body of the 5-DOF "TriVariant" hybrid robot and discussed its design and synthesis.Its kinematic analysis was done, and observations about its performance were made.Lu et al. [26] proposed a novel hybrid manipulator composed of a 5-DOF parallel mechanism with two composite equivalent universal joints and three 1-DOF finger mechanisms.Its kinematics and dynamics are derived, and a 3D prototype of the proposed hybrid manipulator is also constructed and analyzed.Sangveraphunsiri and Chooprasird [27] proposed a 5-DOF hybrid manipulator having three translational and one rotational movement, with a single-axis rotating table.Its kinematics and dynamics are derived, where the latter is derived using the Lagrange formulation.Lingmin Xu et al. [28] presented a one-of-its-kind 3T1R 2-(PRR)-2-RH hybrid manipulator that uses two helical joints, having an equal but opposite pitch, connected to the moving platform, which imparts it full rotational capability.Peng Xu et al. [29] presented the kinematic analysis of a novel hybrid manipulator composed of mainly three parts: a 3-DOF translational parallel mechanism, a 2-DOF serial mechanism, and a turntable providing a redundant DOF.They noted its essential application in computer-controlled ultra-precision freeform polishing.A novel 5-DOF hybrid manipulator was proposed by Guo et al. [30], and its closed-form kinematic solutions were presented.More recently, He et al. [31] proposed a foldable 7-DOF hybrid manipulator to capture a noncooperative target in space.It has two parts, the upper part being a 4R serial mechanism and the lower part being a three-legged parallel manipulator with one UP leg and 2 UPS legs.In another work [32], a hybrid parallel-serial manipulator named CaHyMan (Cassino Hybrid Manipulator) is analyzed in terms of stiffness characteristics.From all these hybrid manipulators, it is clear that some hybrid manipulators [19,21] make it possible for the upper platform to have a pure 3-DOF rotational motion while maintaining some desired translation with respect to the base (due to the lower parallel mechanism being of pure translatory nature), but not the other way round.
However, in industries and manufacturing facilities, there is a frequent requirement of having a pure translational motion of the top platform with respect to an oblique plane at some arbitrary desired orientation with respect to the fixed base, such as in polishing, finishing, painting, or cleaning operations, especially where large structures are involved.One specific case could be of polyhedral domes or surfaces such as the geodesic dome or even the polyhedral structures used in deep space technology, including telescopes and satellites.Based on this motivation, this paper proposes a novel 6-DOF hybrid manipulator in which the top platform can perform 3-DOF pure translational motion with respect to any desired oblique plane (i.e., the mid-platform).This hybrid manipulator utilizes two 3-DOF parallel mechanisms: the lower one being the 1-RRR 2-SPS mechanism that resembles the tripod mechanism proposed by Tanev [20], and the upper one being the 3-UPU translational parallel manipulator (TPM) [2].Although, the upper platform H of the proposed hybrid manipulator can perform a general 3-DOF pure translational motion with respect to the mid-platform while maintaining a constant orientation with respect to the fixed base as shown in Fig. 1(a), however, of particular interest and usefulness could be the special case of planar translation along a specified direction in a plane parallel to the desired oblique plane (i.e., mid-platform M) which is at some specified offset and orientation with respect to the fixed base.The schematic of this case is shown in Fig. 1(b).This desired motion can be decoupled and expressed in four parameters: the desired orientation + offset of the mid-platform with respect to the base, the perpendicular distance between the middle and the upper platform, and the desired translation direction.So, first, the lower parallel mechanism orients and translates simultaneously to achieve the desired pose (orientation and offset).After that, the upper parallel mechanism first attains the required perpendicular distance and then performs the desired translational motion in the specified direction in a plane parallel to the desired oblique plane.Note that the maximum possible offset value is constrained by the value of the revolute limb length of the lower parallel mechanism.The structure and configuration of the proposed hybrid manipulator allow for a convenient decoupling in its inverse kinematic analysis, which facilitates the attainment of closed-form solutions that, in turn, would make the control of this manipulator computationally efficient.
The remainder of this paper is organized as follows.The geometry and structure of the hybrid manipulator are described in Section 2. It is followed by its forward and inverse kinematic analysis in Sections 3 and 4, respectively.The 3D modeling, simulation, and numerical verification are presented in Section 5. Section 6 presents the singularities, while Section 7 lays the ground for more rigorous studies to get better optimal designs for the proposed hybrid manipulator.

Description of novel 6-DOF hybrid manipulator 2.1. Geometry and structure
The proposed hybrid manipulator, as shown in Fig. 2(a), consists of two 3-DOF parallel mechanisms connected serially, making it an overall 6-DOF hybrid manipulator.All the prismatic pairs (P 2 , P 3 , P 4 , P 5 , P 6 ) and the R 2 revolute joint are actuated (active).The two platforms, M and H, are taken to be equilateral, with circumradii of h 1 and h 2 ( = h 1 )1 , respectively.The lower parallel mechanism resembles the upper parallel mechanism used by Tanev [20] in his hybrid manipulator.The lower parallel mechanism connects the base to the mid-platform M and comprises one RRR and two SPS legs.Here, axis(R 1 )2 ⊥ axis(R 2 ) ⊥ axis(R 3 ) and also axis(R 1 ) ⊥ B 1 B 2 edge of the lower triangular base.Further, axis(R 3 ) M 1 M 2 edge of the mid-platform.The upper parallel mechanism between the mid-platform M and the upper platform H (end-effector) is the standard 3-UPU TPM (also called the Tsai [2] manipulator) having three identical UPU legs.Tsai [2] first proposed the conditions required for the 3-UPU manipulator to perform pure translation, and later Gregorio et al. [1] derived the general TPM conditions using the proximal (modified) D-H convention [33], starting from the generic RRPRR3 leg.This paper uses the distal (classical) D-H convention [34].The distal form of the TPM conditions [1,2] are imposed on the upper parallel mechanism in the following manner: , and U 3,1 U 6,1 , U 3,2 U 6,2 , where U i,j denotes the jth revolute joint axis (while traversing a limb upward, i.e., from bottom to the top) of the ith universal joint • For each of the three UPU legs, θ 3 = −θ 2 and θ 4 = −θ 1 , where θ i is the rotation about the ith revolute joint of that particular limb while traversing upward • U 1,2 , U 2,2 , and U 3,2 point toward the circumcenter of platform M • U 4,1 , U 5,1 , and U 6,1 point toward the circumcenter of platform H • U 1,2 and U 4,1 are ⊥ axis(P 4 ), U 2,2 and U 5,1 are ⊥ axis(P 5 ) leg, and U 3,2 and U 6,1 are ⊥ axis(P 6 ) leg The TPM conditions ensure that platform H does pure translation with respect to platform M. Thus, this structure of the hybrid manipulator allows the upper mechanism to perform pure translation motion while being in any desired orientation defined by the mid-platform M, with respect to the base.Furthermore, since the lower parallel mechanism can also translate, platform H can perform pure translation motion even in an oblique plane having some offset with respect to the base.
It is important to note that the choice of the lower mechanism as a 1-RRR 2-SPS mechanism provides the benefit of a straightforward and convenient interpretation for the Euler angles of the mid-platform M in terms of R 1 , R 2 , and R 3 revolute joint angles.From the point of view of controlling the manipulator, this can be pretty useful as encoder feedback from the three revolute joints would directly give the orientation.Furthermore, choosing the upper mechanism as a translational 3-UPU mechanism decouples the orientation and translation of the top platform H, as now the orientation part is entirely governed by the lower mechanism actuators.This helps it achieve a fully closed-form solution structure.In addition to the various applications mentioned in Section 1, this proposed hybrid manipulator offers much flexibility as it can be converted to a standard 3-UPU manipulator just by fixing the lower mechanism, further diversifying its use cases.Finally, as seen in Section 7, even a basic optimization leads to a pretty good design obtaining relatively high manipulability for large parts of the motion range.

D-H analysis
The B 1 M 1 H 1 limb of the hybrid manipulator with all the attached D-H coordinate frames is shown in Fig. 2(b).The frame 1 (X 1 Y 1 Z 1 ) is attached to the fixed base (link #0), while the frame 8 is attached to the top moving platform H (link #7).Two separate frames 4 l and 4 u fixed to the middle platform M (link #3) are used for convenience in the analysis, which will be discussed in subsequent sections.Since the axes of the two revolute joints in a universal joint intersect, the lengths of links #4 and #6 are taken to be zero as they are the imaginary links connecting the two intersecting revolute joints of a universal joint. 4ow, Table I of the proposed hybrid manipulator.
Now, T 1 = f (θ 1 , θ 2 , θ 3 ) corresponds to the D-H Table I, while T 2 = f (θ 4 , θ 5 , L 4 ) corresponds to the D-H Table II.Therefore, the overall transformation matrix is written using the D-H Tables I and II as (details in Section A.1): Now, any position vector of a point expressed 5 in frame 8 can be transformed to the base frame using T from Eq. ( 2) as: It is to be noticed that the orientation part 6 of T depends only on θ 1 , θ 2 , and θ 3 , owing to the pure translational nature of the upper mechanism, which in turn facilitates the decoupling as shown later in Section 4.

Closed-form forward position analysis
In the direct kinematics problem, the values of all the active joints (here θ 2 , L 2 , L 3 , L 4 , L 5 , L 6 ) are specified, and the pose (position + orientation) of the end-effector (upper platform H) needs to be determined.The pose of the end-effector can be entirely described by the overall transformation matrix T because once T is determined, the position vector7 of the vertex H 1 ≡ O 8 as well as the orientation of the upper platform can be extracted from it.Using these two quantities, the position vector of the center of the upper platform can then easily be calculated, which effectively completes the description of the end-effector pose.Now, it is known from Section 2.2 that T contains both active (known) and passive (unknown) joint variables.The passive joint variables can be determined using the leg constraints for the other two limbs, that is, B 2 M 2 H 2 and B 3 M 3 H 3 .A two-stage procedure is adopted as discussed below -first determining the passive joint variables θ 1 and θ 3 of the lower mechanism in terms of the active joint variables (the values of which are given) and then using the determined values along with the active joint variables to solve for the passive joint variables θ 4 and θ 5 of the upper mechanism.The overall transformation matrix T can now be fully determined using these values.

Lower parallel (1-RRR 2-SPS) mechanism
The values of θ 2 , L 2 , L 3 are given, and the values of the passive joint variables θ 1 and θ 3 are required.Now, from the geometry of the manipulator, the following vectors8 can be written: Now the vectors in Eqs. ( 6) and ( 7) are transformed to frame 1 as: The leg vectors 1 L 2 and 1 L 3 are written as: Using the given leg lengths: Using Eqs. ( 4)- (11), the above two equations are written in the following form: where, Since Eq. ( 14) is only a function of θ 1 , it can be solved by using half-angle substitution to get two solutions (corresponding to ±) for θ 1 as: Here, tan −1 gives a solution in the range − π 2 , π 2 .Substituting the obtained value of θ 1 in Eq. ( 15) gives the following equation: where, The solution for θ 3 using Eq. ( 18) is given as: In total, 2 × 2 = 4 solution sets of θ 1 and θ 3 are obtained.For each solution set, the transformation matrix T 1 can be fully determined, which completes the forward kinematics of the lower parallel mechanism.

Upper parallel 3-UPU TPM mechanism
For analyzing the 3-UPU mechanism, a similar procedure is adopted as done for the lower parallel mechanism.Here, the values of L 4 , L 5 , L 6 are given, and the values of the passive joint variables θ 4 and θ 5 are determined.The components of the relevant vectors obtained from the geometry of the manipulator are given below: Now, the vectors in Eqs. ( 23) and ( 24) can be transformed to frame 4 u as: The leg vectors 4u L 5 and 4u L 6 can be written as: The leg length constraint equations are: Using Eqs. ( 21)- (28) in Eqs. ( 29) and ( 30), then adding Eqs. ( 29) and (30), and subtracting Eqs. ( 30) from (29), yields the following two equations: where Since Eq. ( 31) is only a function of θ 5 , the two solutions are written as: Now substituting the value of θ 5 obtained in Eq. ( 32), the following equation is obtained: where, The two solutions for θ 4 from the above equation are written as: In total, 2 × 2 = 4 solution sets of θ 4 and θ 5 are obtained.For each solution set, the transformation matrix T 2 can be fully determined, which completes the forward kinematics of the upper parallel mechanism.However, these four distinct joint space solution sets for the upper parallel mechanism result in only two distinct poses, that is, only two distinct task space configurations.The reason for this lies in the interdependence of the four joint space solution sets, that is, say the first solution set is {θ 4 = β 4 , θ 5 = β 5 }, then the other three solution sets are given as {(β 4 − π ), (π − β 5 )}, {−β 4 , β 5 }, {(π − β 4 ), (π − β 5 )}.So rotating about the two revolute axes of the universal joint U 1 following the first solution set values for θ 4 and θ 5 , and then again following the second solution set values leads to the same final configuration.
Similarly, the third and the fourth solution sets lead to the same final configuration, but it is different from the one attained by the first two solution sets.Therefore, only two distinct task space configurations are obtained, as would be numerically illustrated later in Section 5.2.1 as well.This fact can also be directly deduced easily by substituting the solution sets into the symbolic form of T 2 and observing that there are indeed only two distinct T 2 matrices, and hence only two distinct poses.Moreover, the two obtained task space configurations are related to each other -if the upper mechanism of one of the configurations is reflected in the mid-plane, the other configuration is obtained.In other words, both task space configurations are just the reflection of the upper mechanism in the mid-plane, with the lower mechanism configuration being the same.

Determining the overall transformation matrix
Overall, 4 × 4 = 16 distinct solution sets of (θ 1 , θ 3 , θ 4 , θ 5 ) quadruplet of passive variables are obtained that correspond to a total of 16 possible direct kinematic solutions in joint space for the complete hybrid manipulator.However, for the above-mentioned reason in Section 3.2, there are only 4 × 2 = 8 direct kinematic solutions in task space.
Having determined all the passive variables and consequently evaluating the matrices T 1 and T 2 , the overall transformation matrix T 2 can be determined using Eq. ( 2).Further, the position vector of vertex H 1 , as well as the orientation of the upper platform, can be extracted from it, which can be used to get the position vector of any other point (including the center) of the upper platform.Therefore, this completes the forward kinematic analysis of the proposed hybrid manipulator.

Closed-form inverse position analysis
The most crucial part in the control of any manipulator is the efficient and accurate calculation of the inverse kinematics of the manipulator, in which the required joint motions are to be found from the given pose of the end-effector.In the context of the present hybrid manipulator, given the pose of the upper platform H [i.e., given the overall transformation matrix T of Eq. ( 2)], the active joint variables θ 2 , L 2 , L 3 , L 4 , L 5 , L 6 need to be determined.Owing to the particular structure of the hybrid manipulator, the orientation part of T contains the active and the passive variables of only the lower parallel mechanism; this facilitates the decoupling of the inverse kinematic analysis into two separate parts: one for the lower mechanism and the other for the upper mechanism, as discussed below.Now, let the given numerical overall transformation matrix be T num .

Lower parallel (1-RRR 2-SPS) mechanism
Let the orientation part of T num be: From the components of T described in symbolic form in Section 2.2, the components of the orientation part are extracted to get the orientation matrix Q in symbolic form, which can then be compared with the given numerical orientation matrix Q num to solve for the active and the passive variables of the lower mechanism.So using the fact that Q ≡ Q num , and comparing their (3,1), (3,2), and (3,3) elements, the following three equations are obtained: where the RHS 9 is known, while the LHS is unknown.Multiplying Eq. (39) by √ 3 and adding it to Eq. (41), followed by suitable manipulation gives Now, dividing the RHS of Eq. ( 40) by LHS of Eq. (43) (provided s 2 = 0) gives Therefore, two solutions10 (having a difference of π ) are obtained for θ 3 .In order to get θ 2 , the obtained θ 3 value is substituted in Eq. ( 40) to get: From Eqs. ( 42) and (45), the value of θ 2 can be obtained: Now, in order to get θ 1 , the (1,2) and (2,2) elements of Q and Q num are compared and the calculated values of θ 2 and θ 3 are substituted to get the following two equations in s 1 and c 1 : Eliminating c 1 using Eqs.( 47) and ( 48), an equation in s 1 of the following form is obtained: where the coefficients are in terms of the known quantities, Substituting back s 1 in Eq. ( 47), an equation in c 1 is obtained: where the coefficients are in terms of the known quantities, From Eqs. (49) and (51), the value of θ 1 can be obtained: Finally, overall, two solutions are obtained, since for each of the two possible θ 3 values, unique solutions for θ 1 and θ 2 exist.Also, now the values of θ 1 , θ 2 , and θ 3 can be substituted into Eqs.( 10) and (11) to obtain the leg length vectors, and then L 2 and L 3 can be calculated using Eqs.( 12) and ( 13), respectively, as: It is to be noted that, for the reasons stated earlier 6 on Page 6, the limb lengths L 2 and L 3 are only taken to be positive. 11his completes the inverse kinematics for the lower parallel mechanism since there are now the values for θ 1 , θ 2 , θ 3 , L 2 and L 3 .Further, the numeric T 1num matrix can be obtained by substituting all these values in the symbolic form of T 1 .

Upper parallel 3-UPU TPM mechanism
Now, since the inverse kinematics is solved for the lower parallel mechanism, we can decouple the overall numeric transformation matrix T num using Eq. ( 2) and get the transformation matrix T 2num having the following form as: Moreover, the symbolic form of T 2 matrix found out using the D-H Table II is Thus, exploiting the fact that T 2 ≡ T 2num , and comparing their (1,4), (2,4) and (3,4) elements, the following three equations are seen: Squaring and adding Eqs. ( 58), ( 59) and (60), L 4 is obtained as: Also, dividing Eqs. ( 59) and (58), two solutions for θ 4 are obtained as: Now, using the values of L 4 and θ 4 in Eqs. ( 59) and (60), s 5 and c 5 , and therefore, θ 5 are found out as: So, again, there are overall two solutions, since for each of the two possible θ 4 values, there exist unique solutions for L 4 and θ 5 .However, the two solution sets of θ 4 and θ 5 are of the form {β 4 , β 5 } and {(β 4 − π ), (π − β 5 )}, which again lead to same task space configuration, for reasons mentioned earlier in Section 3.2.Thus, only one distinct task space configuration exists corresponding to the two inverse kinematic solutions of the upper parallel mechanism.Now, the values of θ 4 , θ 5 , and L 4 can be substituted into Eqs.( 27) and ( 28) to obtain the leg length vectors, and then L 5 and L 6 can be calculated using Eqs.( 29) and (30), respectively, as: (66) It is to be noted that, here also, for the reasons stated earlier 6 on Page 6, the limb lengths L 4 , L 5 , and L 6 are only taken to be positive.12Finally, 2 × 2 = 4 possible inverse kinematic solutions are obtained for the complete hybrid manipulator in joint space.However, only 2 × 1 = 2 inverse kinematic solutions are seen with distinct task space configurations.This fact will also be numerically illustrated later in Section 5.2.2.
Thus, this completes the inverse kinematic analysis of the hybrid manipulator as all the actuated joint variables θ 2 , L 2 , L 3 , L 4 , L 5 , L 6 have been obtained from the given pose of the end-effector.

Model
In the literature, different simulation environments have been developed and used, some of which, such as the RoboSim and RoboAnalyzer , are specialized in robotics.In contrast, others are generic 3D simulation software, such as Adams TM , MapleSim TM , and Simscape TM Multibody TM (earlier SimMechanics TM ).In this paper, the proposed hybrid manipulator is modeled and simulated using Simscape Multibody, which is based on Simulink and MATLAB .The advantage of using Simscape Multibody is that sensing various signals, including joint movements and displacements, torques, and forces, becomes much more accessible.
The 3D model of the hybrid manipulator is shown in Fig. 3, while Fig. 4(a) shows the complete hybrid manipulator block diagram built in Simscape Multibody.The whole model is made up of five parts: base, lower parallel mechanism, platform M, upper parallel mechanism, and platform H, where the lower mechanism and upper mechanism models, as shown in Fig. 4(b) and (c), respectively, have been built in accordance with the geometry and structure presented in Section 2.1.

Forward kinematics
The following active joint variable values are noted:   III.

Simulation example: translation in a specific direction in a plane parallel to an oblique plane
Consider the following example in which a planar translation of the upper platform is desired along a specified translational direction while maintaining some constant orientation with respect to the base.As described in Section 1, four parameters can be specified, of which the offset distance and the required orientation of the mid-platform are embedded in the given T 1num matrix:  IV.
The desired perpendicular distance is 50 cm.The translational motion is expected to be along the X 8 ≡ X 4u axis, such that it first translates 50 cm along −X 8 axis, and then from that point onward it translates 100 cm along the +X 8 axis, effectively ending the motion at a displacement of 50 cm along +X 8 axis.Furthermore, the same design parameters were used as mentioned in Eq. ( 68) in Section 5.2.
For the simulation, inverse kinematics was performed on the desired motion to get the actuated joint values, which were then fed to the model.The resulting simulation is presented in Fig. 7, showing the three stages of motions.Also, the simulated translational motion measured using the Transform Sensor Block in Simscape Multibody is shown in Fig. 8, and it clearly matches the desired motion.
During the simulation, the joint limits were added appropriately to only allow for feasible motions.

Singularity analysis
Singularity analysis of serial and parallel manipulators has been quite extensively studied in the literature, and different types of singularities have been identified [6].In the following subsections, the singularity analysis of the proposed hybrid manipulator is presented, and all of its singularities in the joint space of active joints are determined.Consider the following vectors defined as: ] T : vector of all active joint variables = [θ 1 (t), θ 3 (t), θ 4 (t), θ 5 (t)] T : vector of all passive joint variables in the formulation, i.e., along limb B 1 M 1 H 1 q = [θ 1 (t), θ 2 (t), θ 3 (t), θ 4 (t), θ 5 (t), L 4 (t)] T : vector of all the joint variables (active + passive) in the formulation, i.e., along limb B 1 M 1 H 1 $ = ω x , ω y , ω z , v x , v y , v z T : the twist vector of the end-effector (platform H) First, the singularities occurring due to the parallel structure of the hybrid manipulator are determined, which we loosely call the "parallel singularities."Differentiating the closed-loop constraint equations ( 12), ( 13), (29), and ( 30) with respect to time, and also introducing two dummy equations: θ2 (t) = θ2 (t) https://doi.org/10.1017/S0263574723001662Published online by Cambridge University Press  and L4 (t) = L4 (t), which are trivially satisfied, a convenient and compact matrix form for the six equations is as follows: Here, D = D ( (t)) is a diagonal matrix with the diagonal as [1, L 2 (t), L 3 (t), 1, L 5 (t), L 6 (t)], and J p = J p (q(t)) is a 6 × 6 jacobian matrix having the following form (details in Section A.2): https://doi.org/10.1017/S0263574723001662Published online by Cambridge University Press It is to be noted that the J p matrix as defined here is closely related to the K * matrix definition given by Ghosal [35] in their analysis of the gain singularities.More specifically, det(J p ) = det(K * ), and thus in accordance with [35], det(J p ) = 0 gives Type A (gain) singularity.

Type A singularity: DOF-gain singularity
For the gain singularity to occur, det(J p ) = 0 should hold.From Eq. ( 74), it is known that: det(J p ) = j p 21 j p 33 j p 64 j p 55 − j p 54 j p 65 = 0 (75) So, from Eq. ( 75), the following different cases are seen for the gain singularity to occur: Here, b 2 = 0 is a design singularity, which, however, seems practically infeasible.Now, from Eqs. ( 76) and ( 14), there are two equations in θ 1 , θ 2 , and L 2 : where the coefficients are in terms of the known quantities, Solving Eq. (77) to get two solutions (with a difference of π ) for θ 1 : Substituting each value of θ 1 in Eq. ( 78), two separate implicit equations are obtained for the two branches of the singularity curve in θ 2 -L 2 plane.For the numerical example presented earlier in Section 5.2, the two branches of the singularity curve are depicted in Fig. 9.
Case 2: j p 33 = 0 Here, again, h 1 = 0 is a design singularity, which is not feasible practically.Again Eq. ( 14) ≡ f (θ 1 , θ 2 , L 2 ) is used to get two solutions for θ 1 ≡ f (θ 2 , L 2 ).Substituting θ 1 back in Eq. (81) ≡ f (θ 1 , θ 2 , θ 3 ), there are now two solutions for θ 3 ≡ f (θ 2 , L 2 ).Finally, putting the values of θ 1 and θ 3 into Eq.( 15) ≡ f (θ 1 , θ 2 , θ 3 , L 3 ), an equation for a branch of the singularity surface in actuated joint space of θ 2 , L 2 , and L 3 is obtained.Since 2 × 2 = 4 overall solutions for θ 1 and θ 3 exist, a total of four branches are seen of the singularity surface in θ 2 -L 2 -L 3 actuated joint space.For the numerical example presented earlier in Clearly, h 1 = h 2 gives a feasible design singularity.Additionally, the following singularities are also obtained: (1) L 4 = 0 (2) θ 5 = n + 1 2 π where n = 0, ±1, ±2, . .., occurs when M 1 H 1 limb becomes parallel or antiparallel to the U 1 1 axis (3) θ 4 = nπ where n = 0, ±1, ±2, . .., occurs when M 1 H 1 limb comes in the plane of mid-platform M Here, (b) yields the same singularity surface in L 4 -L 5 -L 6 actuated joint space as shown in Fig. 12. From (c), |cos(θ 4 )| = 1.So, using Eqs.(31) and (32): Thus, Eq. ( 85) gives the singularity surface in L 4 -L 5 -L 6 actuated joint space.For the numerical example presented earlier in Section 5.2, the singularity surface is shown in Fig. 11.Now, the "serial singularities" are analyzed occurring because the two platforms are serially connected by a serial kinematic chain.Consider the end-effector point to be at the centroid of platform H.The right-invariant angular velocity matrix 1  e is written as:   T since: Now, the position vector of the end-effector is written as 8 p e = [h 2 , 0, 0, 1] T and is transformed into base frame as 1 p e = T 8 p e .Differentiating 1 p e with respect to time, the velocity 1 v e of end-effector 1 v e is finally obtained as the first three entries (i.e., ignoring the fourth homogeneous coordinate) of 1 p e .Thus, now the twist $ e of the platform H about the global coordinate system at the end-effector point can be found out by concatenating 1 ω e and 1 v e .Further simplification and rearranging gives the following equation: where J s = J s (q(t)) is another 6 × 6 jacobian matrix with the following form (details in Section A.2): Combining Eqs. ( 73) and (88), while assuming J p as nonsingular: Again, the matrix J s J p −1 D ≡ J eq matrix definition given by Ghosal [35].Thus, in accordance with the condition described in ref. [35], since det(J eq ) = 0 if any one of J s or D is singular, det(J s ) = 0 and det(D) = 0 give Type B (loss) singularity.

Type B singularity: DOF-loss singularity
From D matrix definition, det(D) = L 2 L 3 L 5 L 6 = 0 gives the trivial conditions L i = 0 for i = 2, 3, 5, 6 having obvious geometric interpretation (i.e., corresponding limb length becomes zero).Now, from Eq. ( 89), the determinant of J s matrix is det(J s ) = j s 44 j s 55 − j s 45 j s 54 j s 66 + j s 54 j s 65 − j s 55 j s 64 j s 46 + j s 45 j s 64 − j s 44 j s 65 j s 56 (91) • j s 12 j s 23 − j s 13 j s 22 (92) Substituting the expressions presented earlier for j sxy and further simplifying, the condition for singularity is written as: Again, L 4 = 0 is the trivial case.For the case sin(θ 2 ) = 0, the condition on θ 2 is θ 2 = nπ where n = 0, ±1, ±2, . . .
This has a straightforward geometric interpretation that whenever the R 2 revolute joint (i.e., θ 2 ) is actuated such that the axes of R 1 and R 3 (i.e., Z 1 and Z 3 ) become parallel or antiparallel, this type of singularity occurs.
For the case cos(θ 5 ) = 0, the condition on θ 5 is which implies that sin(θ 5 ) = ±1.Substituting in Eq. ( 31), two branches of the singularity surface in L 4 -L 5 -L 6 actuated joint space are seen: since f 1 and f 2 , as defined in Eq. ( 33), are f (L 4 , L 5 , L 6 ).For the numerical example presented earlier in Section 5.2, the two branches of the singularity surface are depicted in Fig. 12. Geometrically, this type of singularity is seen when M 1 H 1 limb becomes parallel or antiparallel to the U 1 1 axis.

Type C singularity: special case
As can be seen from Sections 6.1 and 6.2, singularities L 4 = 0 and θ 5 = n + 1 2 π satisfy the DOF-gain and DOF-loss conditions simultaneously.Hence, they are termed as the "special case" of singularities.This concludes the determination of the singularity conditions for the proposed hybrid manipulator.

Kinematic manipulability and optimal design
As seen in Section 6, the manipulator is in a singularity when the determinant of the Jacobian approaches zero.However, it is not desirable to even go close to being at a singularity when the manipulator is in motion.We use performance indicators such as dexterity and manipulability to quantify this "closeness" to a singularity.Dexterity (denoted as w d ) is defined as the determinant of the manipulator Jacobian J eq , while manipulability (denoted as w m ) is the square root of the determinant of the product of the manipulator Jacobian J eq by its transpose.However, since J eq is square, the dexterity and the manipulability are equal in magnitude:  The manipulability14 index gets its importance from the fact that it represents the volume of the manipulability ellipsoid formed by the end-effector.A higher volume indicates it is more toward isotropy, while a vanishing volume signifies a singularity.Consider the same design parameters and the active joint variables for the upper mechanism as in Section 5.2.The manipulability index as a function of the active joint variables of the lower mechanism for three different values of θ 2 has been plotted in Fig. 13.
Similarly, taking the same active joint variables for the lower mechanism as in Section 5.2, the manipulability index as a function of the active joint variables of the upper mechanism for three different values of L 4 has been plotted in Fig. 14.
While a more rigorous optimization analysis could be done as in ref. [3] taking into account the workspace and condition number as well, we present, in this paper, a preliminary analysis to obtain an optimal design for the proposed hybrid manipulator considering just the kinematic manipulability as the determining factor.As can be seen from Figs. 13 and 14, for a given set of parameters, we have the manipulability varying in order of magnitude throughout the motion range.We define "optimal" design as the set of parameters with the highest minimum value of manipulability throughout the motion range, that is, which is farthest away from singularity for the entirety of its motion range.Therefore, we formulate the problem of optimal design as a bilevel nested continuous maximization-minimization problem: max This is clearly a high-dimensional complex problem.To make the analysis a bit simpler at the cost of losing full optimality, we split it into two lower-dimensional optimization problems as follows: Furthermore, we perform optimization using the MATLAB Global Optimization Toolbox's Genetic Algorithm solver.To further reduce computational complexity, we discretize the optimization with a small enough step size, and the parameters and the active variables are taken to be in a reasonably wide parameter range and motion range, respectively.For the optimization in Eq. (100), the parameters in params 2 are assigned heuristically selected suitable values, and then for the optimization in Eq. ( 101), params 1 are assigned their optimal values from Eq. (100).The optimization was carried out, resulting in the following solution:

Conclusion
This paper introduces a new 6-DOF hybrid manipulator, which uniquely combines a 1-RRR 2-SPS mechanism with a 3-UPU mechanism.This novel architecture holds promise for industrial settings that demand both translational motion and adjustable end-effector orientation.We have provided closed-form solutions for its forward and inverse kinematics, facilitating easier control and real-time applications.Performance and singularity analyses have also been conducted, leading us to a preliminary, near-optimal manipulator design.Looking ahead, the manipulator's potential can be further explored through real-time control algorithms and trajectory planning.A comprehensive workspace analysis could offer further insights into its practical applications.Lastly, a comparative study with other hybrid manipulators in the field would help benchmark its performance and identify improvement areas.(A4)

A.2. Singularity matrices
The J p = J p (q(t)) is a 6 × 6 jacobian matrix having the following form:

Figure 1 .
Figure 1.Conceptual representation of the motivation for the proposed hybrid manipulator.(a) 3-DOF pure translational capability of the upper platform H (it is a representative (i.e., conceptual) image since the platforms are, in reality, triangular and not rectangular).(b) Desired translational motion parallel to the oblique plane.

Figure 2 .
Figure 2. Proposed hybrid manipulator: (a) schematic (the centers of upper joints of the lower mechanism (i.e., R 3 , S 4 , and S 3 ) coincide, respectively, with the centers of the lower joints of the upper mechanism (i.e., U 1 , U 2 , and U 3 ).However, they have been depicted as separate for clarity), and (b) B 1 M 1 H 1 limb with D-H frames.

Figure 3 .
Figure 3. Simulated 3D model of the hybrid manipulator in Simscape Multibody (in Simscape Multibody models, joints are not displayed in the final simulated 3D model).

Figure 4 .
Figure 4. Model of the proposed hybrid manipulator built in Simscape Multibody.

Figure 5 .
Figure 5. Eight distinct task space solutions, with each representing a pair of joint space solutions in TableIII.

Figure 6 .
Figure 6.Two distinct configurations, with each representing a pair of solutions in TableIV.

Figure 7 .
Figure 7. Simulation of the numerical example, showing the three different stages of motion.

Figure 12 .
Figure 12.The two branches of the singularity surface in L 4 -L 5 -L 6 actuated joint space.

Figure 13 .
Figure 13.Manipulability index as a function of lower mechanism's active joint variables for three different values of θ 2 .

Figure 14 .
Figure 14.Manipulability index as a function of upper mechanism's active joint variables for three different values of L 4 .

L 1 ,Figure 15 .
Figure 15.Manipulability index as a function of lower mechanism's active joint variables for three different values of θ 2 , for optimal parameters.

Figure 16 .
Figure 16.Manipulability index as a function of upper mechanism's active joint variables for three different values of L 4 , for optimal parameters.

Figures 15 and 16
Figures 15 and 16  show the manipulability index variation for the motion range for optimal parameters.It shows that the manipulability is pretty high in an appreciably large region of the motion range and how superior it is compared to Figs. 13 and 14 for the nonoptimal parameters.

Table I .
represents the D-H table for the B 1 M 1 leg of the lower mechanism, and it transforms frame 4 l to frame 1. Table II represents the D-H table for the M 1 H 1 leg of the upper mechanism after imposing the TPM conditions as listed in Section 2.1, and it transforms frame 8 to frame 4 u .The intermediate (constant) rigid transformation matrix, which transforms frame 4 u to 4 l is D-H table for lower mechanism (along B 1 M 1 limb) of the proposed hybrid manipulator.

Table III .
All 16solutions for the DK of the hybrid manipulator, with each pair of it corresponding to a single configuration.

Table IV .
All four solutions to the inverse kinematics of the hybrid manipulator, with each pair of it corresponding to a single configuration.