# Design, Dynamics, and Workspace of a Hybrid-Driven-Based Cable Parallel Manipulator

Authors: Bin Zi; Jianbin Cao; Zhencai Zhu; Peter Mitrouchev
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

The design, dynamics, and workspace of a hybrid-driven-based cable parallel manipulator (HDCPM) are presented. The HDCPM is able to perform high efficiency, heavy load, and high-performance motion due to the advantages of both the cable parallel manipulator and the hybrid-driven planar five-bar mechanism. The design is performed according to theories of mechanism structure synthesis for cable parallel manipulators. The dynamic formulation of the HDCPM is established on the basis of Newton-Euler method. The workspace of the manipulator is analyzed additionally. As an example, a completely restrained HDCPM with 3 degrees of freedom is studied in simulation in order to verify the validity of the proposed design, workspace, and dynamic analysis. The simulation results, compared with the theoretical analysis, and the case study previously performed show that the manipulator design is reasonable and the mathematical models are correct, which provides the theoretical basis for future physical prototype and control system design.


---

## Body

## 1. Introduction

In manipulation systems, it is always a goal to continuously improve the range of motion, accuracy, and load capacities. Compared with the traditional mechanical structures, parallel mechanisms have the advantages of high dexterity and accuracy. Cable parallel manipulators (CPMs) are categorized as a type of parallel manipulators. They are a variant of the Gough-Stewart platform in which rigid extensible legs are replaced by cables stored spools. CPMs have several advantages over rigid-link mechanisms, such as (i) remote location of motors and controls, (ii) rapid deployability, (iii) potentially large workspaces and (iv) high load capacity; (v) reliability [1, 2]. For the preceding characteristics, the CPMs have received a lot of attention since 1980s [3]. These manipulators, referred to as the overhead crane and rotary crane, are useful for many applications in hazardous environments in order to enable the manipulation of heavy objects [4, 5].The first written information on use of hoist mechanism in ancient Greece appeared around 530 BC, which is concerned with the construction of the first temple of Artemis in Ephesus, and temples at Selinous, temples of Apollo at Syracuse and at Corinth. After 515 BC, cranes were in common use [6, 7]. The CPMs consist of a fixed base and a centrally located end-effector, attached to moving payload, connected by cables whose tension is maintained along the tracked trajectory [8, 9]. However, due to the unilateral driving property of the cables, maintaining positive cable tension is essential to maneuver the end-effector. As a result, the number of driving cables must be more than the number of degrees of freedom (DOF) of the end-effector in order to obtain a completely restrained manipulation [10]. This type of technology is now well known, and several studies of the CPMs have been undertaken. For instance, Gosselin and Barrette [11] presented a fundamental and systematic analysis of a planar CPM. Design and performances of a 4-cable driven parallel manipulator have been presented by Ottaviano et al. [12]. Flatness-based trajectory tracking control for a three-cable suspension manipulator has been investigated by Heyden and Woernle [13]. Due to cable slippage on the drum, Baser and Konukseven developed an analytical method for predicting the transmission error [14]. Different aspects like kinematics [15], workspace characteristics [16], dimensional optimization [17], and statics [18] were extensively investigated.As we are aware, there are relatively fewer published works on the drive system of the CPMs, which is the main contribution of this work. The existing research on the drive system of the CPMs usually uses the low power controllable motor. In the recent years, as the CPMs have been investigated widely for various applications in modern manufacturing, they are required not only for operations with high accuracy and high payload, but also for output with greater flexibility, which can change the law of output motion quickly and conveniently [19]. Nevertheless, the low power controllable motor cannot directly drive the high-loading CPM due to the restrictions of power and torque. Hence, it is necessary to carry out for a new-type driven system for the CPMs.The hybrid-driven planar five-bar mechanism (HDPM) is a kind of machine whose drive system consists of a constant velocity (CV) motor and a servomotor [20]. The CV motor provides main power and motion required; however, it lacks high controllability. On the other hand, the servomotor acts as a motion regulation device which suffers from the limited power capacity. Therefore, the HDPM can take the advantage of the complementary characteristics of both motors to generate a programmable range of highly nonlinear output motions with high power capacities [21]. In fact, the HDPM can be seen as a type of servomechanism [22]. In this investigation the design of a hybrid-driven-based cable parallel manipulator (HDCPM) which combines four groups of HDPMs with a 4-cable parallel manipulator is presented. The HDCPM can take the advantages of the characteristics of both mechanisms to generate a high-performance output motion with high efficiency and heavy payload [23].It is well known that dynamic performance has a great impact on the operation of the manipulator, and this is on the basis of its dynamic control. As demonstrated in [24], servomechanism dynamics constitute an important component of the complete robot dynamics. Therefore, for the HDCPM, the dynamics should not only consider the CPM, but also include the HDPM (servomechanism). However, the literature on the dynamics of the CPM system including the linkage (servomechanism) dynamics is sparse. Aiming at the HDCPM, this paper is intended to study the completely dynamic modeling including the HDPM actuator dynamics on the basis of Newton-Euler method. The latter is used to solve the inverse kinematic and dynamic problems of the HDCPM on condition that the operation path of the end-effector has been planned.The workspace is also one of the most useful measures for the evaluation of the manipulator performances [25–27]. Pusey et al. [28] defined the workspace of a CPM as the set of points where the centre of gravity of the end-effector can be reached with admissible tensions in all cables. Apart from the tension condition, workspace of the HDCPM is also limited by geometrical conditions, which include the workspace of the HDPM actuator, length of cables, and motion range of end-effector for the HDCPM. This paper discusses the process of workspace formation of the HDCPM in detail.In the rest of the paper, Section2 describes the design model of the HDCPM. In the following, dynamic modeling of the HDCPM is performed based on Newton-Euler method in Section 3. Workspace analysis of the HDCPM is provided in Section 4. In Section 5, illustrative simulation studies highlight its performances. Finally, some concluding remarks and future studies are summarized in Section 6.


## 2. Mechanism Description

In this paper, for the purpose of analytical modeling and numerical analysis, the three-dimensional design model of the completely restrained HDCPM with three translational motions is taken as an example (see Figure1). The HDCPM suspends an end-effector by four cables and restrains all motion degrees of freedom for the object using the cable force when the end-effector moves within the workspace. For each cable, one end is connected to the end-effector, the other one rolls through a pulley fixed on the top of the relative cable tower rack and then is fed into the HDPM. The HDCPM comprises of two modules: (1) the CPM consisting of four-cable tower racks, four cables, pulley struts, pulleys, girder, and cargo (i.e., end-effector); (2) four groups of HDPMs containing three-phase asynchronous motors, servomotors, reducers, and double crank five-bar linkage; the asynchronous motors are connected using the pulley transmission mechanisms, while the servomotors and the reducers are linked by couplings. At the same time, the pulley transmission mechanisms and reducers are joined to the double crank planar five-bar linkage. The double crank planar five-bar linkage is made up of a large crank disk, a long connecting rod, a small crank disk, and a short connecting rod. Both long and short connecting rods are articulated with each other and attached to the cable. Moreover, each tower rack of the CPM is equipped with a cable guide pulley. The mechanical configuration of the HDCPM is illustrated in Figure 2.Figure 1
Three-dimensional model of the HDCPM.Mechanical configuration of the HDCPM: (a) front view, (b) top view.
(a)
(b)The working process of the HDCPM is as follows. The CV three-phase asynchronous motors and servomotors are power sources. On the one hand, the servomotor and the reducer are linked by the coupling, then the servomotors through a small crank disk are connected with a short rod connection thus adjusting the output motion of the HDPM. On the other hand, each CV three-phase asynchronous motor is connected with the pulley transmission mechanisms, then it links to a long connecting rod by a large crank disk, thus providing the main power for the HDCPM. These two types of input motions are hinged through the short and long connecting rod, so that the power distribution and other characteristics of the HDPM are improved while ensuring the output motions. Four groups of the HDPM with the same structure produce rotary motion, thus the four cables can be driven in order to realize the output motion trajectory of the end-effector.


## 3. Dynamic Modeling

A simple schematic sketch of the HDCPM structure model with the associated coordinate systems is depicted in Figure3. At the bottom of one cable tower rack a global coordinate system O(XYZ) is established. The end-effector has location coordinates G(x,y,z). The distance between each cable tower rack top Mi(xi,yi,zi) and the end-effector is li(i=1,2,3,4). The four cable tower racks have same heighth and are arrayed in a rectangle on the ground, where deformation is ignored. In order to simplify the model, the cables are treated as a massless body with no deformation, with length s. The bottom of the HDPM corresponds to the joint A as the origin, where a local coordinate system A(x′y′z′) is established.Structure model (a) of the HDPM and schematic sketch (b) of the HDCPM.
(a)
(b)LinkAB driven by a CV motor and link ED driven by a servomotor are the two driving links (Figure 2(a)), the joint C is the output. There are two independent constraints in the HDPM, which can be used to derive the kinematic equation of the HDPM:
(1)l1′cosθ1+l2′cosθ2-l3′cosθ3-l4′cosθ4-l5′=0,l1′sinθ1+l2′sinθ2-l3′sinθ3-l4′sinθ4=0In (1) θ1 and θ4 are known. Consequently θ2 and θ3 can be calculated as
(2)θ3(θ1,θ4)=arccos(CA2+B2)+arctan(μλ),θ2(θ1,θ4)=arctan(μ+l3′sinθ3λ+l3′cosθ3),
where A=2l3′λ, B=2l3′μ, C=l2′2-l3′2-λ2-μ2, λ=l4′cosθ4-l1′cosθ1+l5′, μ=l4′sinθ4-l1′sinθ1.Let(3)xC′=l1′cosθ1+l2′cosθ2=l5′+l3′cosθ3+l4′cosθ4,zC′=l1′sinθ1+l2′sinθ2=l3′sinθ3+l4′sinθ4
be the coordinates of joint C in the local coordinate system. According to the geometric relationship of the HDCPM, l1 can be derived as
(4)l1=S-(p+l5′-xC′)2+(h-zC′)2,
where p represents the distance EO¯ and s is the total length of cable. Similarly, l2, l3, and l4 can be derived, respectively.The relationships between the cable lengthli and the end-effector location G(x,y,z) can be easily obtained as follows:
(5)l1=x2+y2+(z-h)2,l2=(x-a)2+y2+(z-h)2,l3=(x-a)2+(y-b)2+(z-h)2,l4=x2+(y-b)2+(z-h)2
resulting in
(6)x=l12-l22+a22a,y=l12-l42+b22b,z=h-l12-x2-y2.Newton’s law for the end-effector of the HDCPM results in(7)m[x¨y¨z¨]=∑i=14Ti→+[00-mg],
where Ti(i=1,2,3,4) is a tension cable force, that is, the driving force of cable i exerted on the end-effector for the HDCPM, and Ti→ is force vector.Thus the dynamic model can be expressed as follows:(8)N3×4T4×1=M3×1,
where
(9)N3×4=[-xl1a-xl2a-xl3-xl4-yl1-yl2b-yl3b-yl4h-zl1h-zl2h-zl3h-zl4]
is the matrix of coordinate transformer, T4×1=[T1T2T3T4]T is the vector of cable tension, M3×1=[mx¨my¨mg+mz¨]T is the vector of inertia force, m is mass of the end-effector, and g the gravitational acceleration.The dynamic model of the HDCPM is presented in two parts. The first is directed to the structural model (CPM) above, the other one to the actuator dynamics (servomechanism). Applying the Newton-Euler equation establishes the dynamic model of the HDPM. The force analysis of each link of the HDPM is shown in Figure4. τ1 and τ4 represent the driving torque of the CV motor on link AB and the servomotor on link ED, respectively; β= arctan  ((h-l1′sinθ1-l2′sinθ2)/(p+l5′-l1′cosθ1-l2′cosθ2)) is the angle between cable tension and x′-coordinate.Free body diagrams of the HDPM: (a) linkAB, (b) link BC, (c) link CD, and (d) link ED.
(a)
(b)
(c)
(d)By the force analysis mentioned previously, the Newton-Euler formulation of the links can be listed as follows.LinkAB:
(10)FAx′-FBx′=m1x¨1′,FAz′-m1g-FBz′=m1z¨1′,τ1+FAx′r1sinθ1-FAz′r1cosθ1+FBx′(l1′-r1)sinθ1-FBz′(l1′-r1)cosθ1=J1θ¨1.LinkBC:
(11)T1cosβ+FBx′-FCx′=m2x¨2′,T1sinβ+FBz′-FCz′-m2g=m2z¨2′,FBx′r2sinθ2-FBz′r2cosθ2+(FCx′-T1cosβ)(l2′-r2)sinθ2+(T1sinβ-FCz′)(l2′-r2)cosθ2=J2θ¨2.LinkCD:
(12)T1cosβ+FCx′-FDx′=m3x¨3′,T1sinβ+FCz′-m3g-FDz′=m3z¨3′,(FCz′+T1sinβ)r3cosθ3-(FCx′+T1cosβ)r3sinθ3-FDx′(l3′-r3)sinθ3+FDz′(l3′-r3)cosθ3=J3θ¨3.LinkED:
(13)FDx′-FEx′=m4x¨4′,FDz′-m4g-FEz′=m4z¨4′,τ4-FDx′r4sinθ4+FDz′r4cosθ4-FEx′(l4′-r4)sinθ4+FEz′(l4′-r4)cosθ4=J4θ¨4,
where mj,  lj′, Jj, rj, (xj′,zj′)(j=1,2,3,4), represent the mass, length, inertia, distance between the mass centre and the joint, and coordinates of the mass centre of each link, respectively.The equation can be summarized as(14)A12×14H14×1=B12×1,where (15)A12×14=[10-100000000000010-10000000000r1s1-r1c1(l1′-r1)s1(r1-l1′)c100000000100010-101000000000010-10100000000r2s2-r2c2(l2′-r2)s2(r2-l2′)c2(r2-l2′)s2(l2′-r2)c200000000001010-100000000001010-100000000-r3s3r3c3-r3s3r3c3(r3-l3′)s3(l3′-r3)c300000000000010-100000000000010-10000000000-r4s4r4c4(r4-l4′)s4(l4′-r4)c401]is the matrix of force transformer, si=sinθi, ci=cosθi.H
14
×
1
=
[
F
A
x
′
,
F
A
z
′
,
F
B
x
′, FBz′,Fcx′,Fcz′, T1cosβ, T1sinβ, FDx′, FDz′,FEx′, FEz′,τ1,τ4]T is the vector of resultant and active force.B
12
×
1
=
[
m
1
x
¨
1
′
,
m
1
g
+
m
1
z
¨
1
′, J1θ¨1,m2x¨2′,m2g+m2z¨2′, J2θ¨2,m3x¨3′,m3g+m3z¨3′, J3θ¨3,m4x¨4′, m4g+m4z¨4′,J4θ¨4]T is the vector of inertia force.Taking into account (8) and (14), the dynamic model of the HDCPM is derived from the dynamic equations of the CPM and four groups of the HDPM as follows:
(16)N3×4T4×1=M3×1,A12×14iH14×1i=B12×1i(i=1,2,3,4).


## 4. Workspace Analysis

The workspace of the HDCPM is characterized as the set of points where the end-effector can be positioned while all cables are in tension (Ti>0). At each point within the possible workspace, (8) describing the force in each cable is used to see if tension is obtainable. However, it is not sufficient to obtain the actual workspace of the HDCPM if it only depends on the tension condition. Accordingly, when analyzing the workspace of the HDCPM, one should also consider the following constraints such as the workspace of the HDPM and the motion range of the end-effector.(1)
The workspace of the HDPM. According to (6), the set of points (x,y,z) where the end-effector can be positioned are influenced by cable length li. Cable length li is calculated by the coordinates joint C(xC′,zC′) from (4). Thus, the workspace of the joint C in x′-z′ plane is one of the constraints of the workspace of the HDCPM, and it can be expressed as (only the first group of the HDPM)
(17)|l1′-l2′|<LAC<l1′+l2′,|l3′-l4′|<LCE<l3′+l4′,
where LAC=xC′2+zC′2, LCE=(xC′-l5′)2+zC′2.(2)
The motion range of the end-effector for the HDCPM can be written in the following form:(18)0<x=l12-l22+a22a<a,0<y=l12-l42+b22b<b,0<z=h-l12-x2-y2<h.A general numerical workspace generation approach is employed here. The possible motion range of jointC in the HDPM is first discretized into a number of points. According to (17), a series of the cable length li can be calculated. Then, from (18), various combinations of li can calculate the end-effector coordinates (x,y,z). Thus, a series of points in Cartesian coordinate system are determined. At the end, the points are checked and those satisfying the tension condition (Ti>0) will generate the workspace. The flowchart of generating the workspace of the HDCPM is shown in Algorithm 1.Algorithm 1:Flowchart of generating the workspace of the HDCPM.
xmin′=-(l1′+l2′),xmax′=l3′+l4′+l5′
z
min
′
=
-
max
{
l
1
′
+
l
2
′
,
l
3
′
+
l
4
′
}
,
z
max
′
=
max
{
l
1
′
+
l
2
′
,
l
3
′
+
l
4
′
}
x
′
=
x
min
′
,
z
′
=
z
min
′
Do While  z′≤zmax′
Do While  x′≤xmax′
If  |l1′-l2′|<LAC<l1′+l2′and|l3′-l4′|<LCE<l3′+l4′  Then
Froml1=s-(p+l5′-xC′)2+(h-zC′)2,l1 is determined.
Similarly,l2,l3, and l4 can be calculated, respectively.
Recordli
End If
x′=x′+Δx′
Loop
z′=z′+Δz′
Loop
Various combinations ofli for calculating coordinates (x,y,z) by
x
=
(
l
1
2
-
l
2
2
+
a
2
)
/
2
a
,
y
=
(
l
1
2
-
l
4
2
+
b
2
)
/
2
b
,
z
=
h
-
l
1
2
-
x
2
-
y
2
and those satisfying motion range of the end-effector(0<x<a,0<y<b,0<z<h)
and tension condition(Ti>0) will generate the workspace


## 5. Results and Discussions

Simulation studies were performed with the software named MATLAB 2010. A three-dimensional simulation model of the HDCPM has been established, where the four groups of HDPMs have the same structure parameters and symmetry in the three-dimensional space. The parameters of the HDCPM are listed in Table1. All computations were carried out on a computer with a 2.99-GHz CPU and 2.00 GB memory; the circle trajectory took about 9.67 s, while the line segment trajectory requires only 2.93 s.Table 1
Parameters of the HDCPM.
System parameters
Value
Mass of the end-effectorm
20 kg
Acceleration due to gravityg
9.81 m/s2
Height of the cable tower rackh
1 m
Side length of four towers distributed in rectanglea
1 m
Side length of four towers distributed in rectangleb
1 m
Length of the cables
2.3 m
Length of theEOp
1 m
Length of the LinkABl1′
0.2 m
Length of the LinkBCl2′
0.5 m
Length of the LinkCDl3′
0.5 m
Length of the LinkEDl4′
0.28 m
Length of the LinkAEl5′
0.51 m
Density of the Links
7.85
×
10
-
3 kg/m3
Cross-sectional area of the Links
2
×
10
-
4 m2


### 5.1. Dynamic Simulation

This section presents two motion cases of the end-effector for dynamic simulation.Case 1 (circle trajectory).
The equation of the circle with a radius of 0.25 m parallel to theX-Y plane in global coordinate system is
(19)x=0.25cost+0.5,y=0.25sint+0.5,z=0.5,
where 0≤t≤6.2(s). The end-effector moves along the trajectory with the constant velocity v=0.25 m/s as shown in Figure 5.Figure 5
Following trajectory of the circle motion.Figure6 shows cable motion in the process of simulation. The cable lengths li are generated by (5) and vary symmetrically. Figure 7 shows the curves of cable tension exerted on the end-effector related to the circle trajectory. Since the end-effector moves on a horizontal circle trajectory, the changes of cable tensions in Figure 7 are reasonable and the transition is smooth.Figure 6
Time behavior of cable length for Case1.Figure 7
Curves of cable tension for Case1.The kinematic parameters of the HDPM can be obtained by the inverse kinematics of the HDCPM, for given end-effector position. The angle of jointA is driven by a CV motor and assumed to be
(20)θ1=ωt,
where ω=1 rad/s. According to Figure 3(b) the following geometrical relation is verified:
(21)(l5′+p-l1′cosθ1-l2′cosθ2)2+(h-l1′sinθ1-l2′sinθ2)2=(s-l1)2,
where l1 is known from (5). Combining (20), and (21), θ2 can be calculated as
(22)θ2=2arctan(E1-F12+E12-G12F1-G1),
where E1=2l2′(l2′sinθ1-h), F1=2l2′(l1cosθ1-l5′-p), G1 = (l5′+p)2- 2l1′(hsinθ1+(l5′+p)cosθ1)+l1′2+l2′2+h2-(s-l1)2.From (3), (20) and (22), the coordinates of joint C(xC′,zC′) can be determined as
(23)xC′=l1′cosθ1+l2′cosθ2,zC′=l1′sinθ1+l2′sinθ2.According to (3), the coordinates of joint C(xC′,zC′) can be rewritten as
(24)xC′=l5′+l3′cosθ3+l4′cosθ4,zC′=l3′sinθ3+l4′sinθ4.Solving (24), one gets
(25)θ4=2arctan(E2-F22+E22-G22F2-G2),θ3=arcsin(zC′-l4′sinθ4l3′),
where E2=-2zC′l4′,F2=2l4′l5′-2xC′l4′ and G2=(xC′-l5′)2+zC′2+l4′2-l3′2.According to the inverse kinematics analysis of the HDCPM above, one can calculate the angular variables for all links of the HDPM as the end-effector moves along the circle. Figure8 shows the angle curves of the links of the HDPM. It can be noted that the angular variable law for link AB is a line, while for the other links are curves. This is in line with the motion law of the HDPM. Moreover, angular curves of all the links maintain good continuity, and the numerical values are reasonable. The angular velocity curves and angular acceleration curves of links can be generated by taking the first and second time derivatives of the angular change (dθi/dt, d2θi/dt2). The coordinates of the joint C in local x′-z′ plane can be computed by (23). The running trajectory of joint C is depicted in Figure 9.Curves of the links motion parameters of the HDPM for Case1: (a) angle, (b) angular velocity, and (c) angular acceleration.
(a)
(b)
(c)Figure 9
Trajectory of jointC in x′-z′ plane for Case 1.Figure10 depicts the resulting force curves of the joints of the HDPM. From it one can see that the resultant force curves of the joints are smooth, which is consistent with the desired motion law.Figure 10
Resultant force curves of the joints of the HDPM for Case1.If the cables are driven only by coils and servomotors, which is a classical approach for the conventional CPM, the range and speed of obtainable output motions may be limited by certain servomotor power capacity and high cost. This problem can be solved by using hybrid-driven machines instead of coils and servomotors. For the hybrid-driven machines, the CV motor provides the majority of power supply to drive the mechanism’s motion. The servomotor is real-time controllable and offline programmable [20]. Therefore, the hybrid-driven mechanisms take the advantage of the complementary characteristics of two types of motors to generate programmable output motions with high power capacities in low cost.The CV motor power and servomotor power of the HDPM can be calculated by(26)PCV1=τ1×θ˙1,PSV1=τ4×θ˙4,
where PCV1 is the CV motor power of the first group of the HDPM and PSV1 is the servomotor power of the first group of the HDPM. Similarly, the other motor power, PCV2, PSV2, PCV3, PSV3, PCV4, PSV4, of the HDPM can be derived, respectively.Figure11 shows motor absolute value power curves of four groups of the HDPMs of the HDCPM. From Figure 11, for each of the HDPM, the range of the CV motor power (PCVi) is larger than that of the servomotor power (PSVi). It is illustrated also that the CV motors undertake the majority of power for the HDCPM system.Figure 11
Motor absolute value power curves of the four groups of HDPMs of the HDCPM.As the motion circular trajectory of the end-effector is the same, and as the structural parameters of the classical approach for the conventional CPM are the same (cables are driven only by coils and servomotors), the servomotor power of the conventional CPM can be obtained by(27)Pi=Ti×l˙i(i=1,2,3,4),
where Pi is the ith servomotor power, Ti represents the ith cable tension, l˙i denotes the ith cable velocity.Figure12 shows the servomotor absolute value power curves of the classical approach for the conventional CPM. For comparison, absolute value power curves of the first group of HDPM of the HDCPM and the first servomotor of the conventional CPM are given in Figure 13. As it is seen, large range of servomotor power |P1| of the conventional CPM may be decreased by using the HDPM.Figure 12
Servomotor absolute value power curves of the conventional CPM.Figure 13
Comparison of absolute value power curves between the HDCPM and conventional CPM.In order to investigate the power property of the HDCPM further, the following simulation is performed as the CV motors of the four groups of the HDPMs of the HDCPM run with 1.1 rad/s (ω=1.1 rad/s). Figure 14 displays curves of the link angle and velocity of the HDPM. Angular curves of all links maintain good continuity, and the numerical values are reasonable. Figure 15 shows absolute value power curves of four groups of the HDPMs of the HDCPM. Figure 16 illustrates absolute value power curves of the first group of HDPM of the HDCPM and the first servomotor of the conventional CPM. From the simulations, we can see that the CV motors provide the majority of power supply and the servomotors mainly act as a motion regulation device as depicted in Figures 11, 13, 15, and 16.Curves of the links motion parameters of the HDPM asω=1.1 rad/s: (a) angle and (b) angular velocity.
(a)
(b)Figure 15
Motor absolute value power curves of the four groups of HDPMs of the HDCPM asω=1.1 rad/s.Figure 16
Comparison of absolute value power curves between the HDCPM asω=1.1 rad/s and conventional CPM.Case 2 (line segment trajectory).
The equation of the line in global coordinate system is(28)x=t6+0.25,y=t12+13,z=t6+0.25,
where 0<t<2(s). From (28), it can be noted that the center of gravity of the end-effector moves the geometric path that is line segment with the constant velocity v=0.25 m/s, the initial point is (0.25, 1/3, 0.25) (m). The trajectory tracking of the end-effector related to Case 2 is shown in Figure 17. From this figure the formulation tracks the planned trajectory relatively well.
Figure18 shows the changes of cables lengths. As can be seen, the variations of the curves of cables lengths are in accordance with the reference path (28) of the end-effector.
Figure19 depicts the curves of cable tensions for Case 2. One can see that the curves are smooth, which is the ideal data for the HDCPM system which operates steadily.
The angular variables of links of the HDPM for Case2 can be obtained by the inverse kinematics analysis of the HDCPM. The angle curves of the HDPM links are shown in Figure 20(a). The corresponding angular velocity and angular acceleration are shown in Figures 20(b) and 20(c). Figure 21 shows the running trajectory of the joint C in local x′o′z′ plane.
The resulting force curves of the joints of the HDPM for Case2 can be solved by the dynamics modeling of the HDCPM in Section 3, as shown in Figure 22.  It can be noted that the resulting force curves of the joints have the same trend and change smoothly. Comparison of absolute value power curves of the HDPMs in the HDCPM and the servomotors of the conventional CPM for Case 2 are shown in Figure 23. Referring to Figure 23, it can be seen that the CV motors provide the majority of power supply and the servomotors mainly act as a motion regulation device.
From the results and discussions, illustrative simulation studies highlight its performances, which lay a foundation for the further research on optimization and real-time control.Figure 17
Following trajectory of the line segment motion.Figure 18
Curves of cable length for Case2.Figure 19
Curves of cable tension for Case2.Curves of the links motion parameters of the HDPM for Case2: (a) angle, (b) angular velocity and (c) angular acceleration.
(a)
(b)
(c)Figure 21
Running trajectory of the jointC in x′-z′ plane for Case 2.Figure 22
Resultant force curves of the joints of the HDPM for Case2.Figure 23
Comparison of absolute value power curves between the HDCPM and conventional CPM for Case2.


### 5.2. Workspace Simulation

The workspace of jointC is determined by geometric constraints (17) given in Figure 24. The simulation results of the reachable 3D workspace for the HDCPM and the projections of workspace onto the X-Y plane and X-Z plane are shown in Figure 25. The whole workspace is approximately an upside-down cone. It should be noted that most of the workspace volume is concentrated in the upper part of the tower cube, which is important for trajectory planning of the end-effector.Figure 24
Workspace of jointC in x′-z′ plane.Workspace of the HDCPM (a) and projections of workspace ontoX-Y plane (b) and X-Z plane (c).
(a)
(b)
(c)


## 5.1. Dynamic Simulation

This section presents two motion cases of the end-effector for dynamic simulation.Case 1 (circle trajectory).
The equation of the circle with a radius of 0.25 m parallel to theX-Y plane in global coordinate system is
(19)x=0.25cost+0.5,y=0.25sint+0.5,z=0.5,
where 0≤t≤6.2(s). The end-effector moves along the trajectory with the constant velocity v=0.25 m/s as shown in Figure 5.Figure 5
Following trajectory of the circle motion.Figure6 shows cable motion in the process of simulation. The cable lengths li are generated by (5) and vary symmetrically. Figure 7 shows the curves of cable tension exerted on the end-effector related to the circle trajectory. Since the end-effector moves on a horizontal circle trajectory, the changes of cable tensions in Figure 7 are reasonable and the transition is smooth.Figure 6
Time behavior of cable length for Case1.Figure 7
Curves of cable tension for Case1.The kinematic parameters of the HDPM can be obtained by the inverse kinematics of the HDCPM, for given end-effector position. The angle of jointA is driven by a CV motor and assumed to be
(20)θ1=ωt,
where ω=1 rad/s. According to Figure 3(b) the following geometrical relation is verified:
(21)(l5′+p-l1′cosθ1-l2′cosθ2)2+(h-l1′sinθ1-l2′sinθ2)2=(s-l1)2,
where l1 is known from (5). Combining (20), and (21), θ2 can be calculated as
(22)θ2=2arctan(E1-F12+E12-G12F1-G1),
where E1=2l2′(l2′sinθ1-h), F1=2l2′(l1cosθ1-l5′-p), G1 = (l5′+p)2- 2l1′(hsinθ1+(l5′+p)cosθ1)+l1′2+l2′2+h2-(s-l1)2.From (3), (20) and (22), the coordinates of joint C(xC′,zC′) can be determined as
(23)xC′=l1′cosθ1+l2′cosθ2,zC′=l1′sinθ1+l2′sinθ2.According to (3), the coordinates of joint C(xC′,zC′) can be rewritten as
(24)xC′=l5′+l3′cosθ3+l4′cosθ4,zC′=l3′sinθ3+l4′sinθ4.Solving (24), one gets
(25)θ4=2arctan(E2-F22+E22-G22F2-G2),θ3=arcsin(zC′-l4′sinθ4l3′),
where E2=-2zC′l4′,F2=2l4′l5′-2xC′l4′ and G2=(xC′-l5′)2+zC′2+l4′2-l3′2.According to the inverse kinematics analysis of the HDCPM above, one can calculate the angular variables for all links of the HDPM as the end-effector moves along the circle. Figure8 shows the angle curves of the links of the HDPM. It can be noted that the angular variable law for link AB is a line, while for the other links are curves. This is in line with the motion law of the HDPM. Moreover, angular curves of all the links maintain good continuity, and the numerical values are reasonable. The angular velocity curves and angular acceleration curves of links can be generated by taking the first and second time derivatives of the angular change (dθi/dt, d2θi/dt2). The coordinates of the joint C in local x′-z′ plane can be computed by (23). The running trajectory of joint C is depicted in Figure 9.Curves of the links motion parameters of the HDPM for Case1: (a) angle, (b) angular velocity, and (c) angular acceleration.
(a)
(b)
(c)Figure 9
Trajectory of jointC in x′-z′ plane for Case 1.Figure10 depicts the resulting force curves of the joints of the HDPM. From it one can see that the resultant force curves of the joints are smooth, which is consistent with the desired motion law.Figure 10
Resultant force curves of the joints of the HDPM for Case1.If the cables are driven only by coils and servomotors, which is a classical approach for the conventional CPM, the range and speed of obtainable output motions may be limited by certain servomotor power capacity and high cost. This problem can be solved by using hybrid-driven machines instead of coils and servomotors. For the hybrid-driven machines, the CV motor provides the majority of power supply to drive the mechanism’s motion. The servomotor is real-time controllable and offline programmable [20]. Therefore, the hybrid-driven mechanisms take the advantage of the complementary characteristics of two types of motors to generate programmable output motions with high power capacities in low cost.The CV motor power and servomotor power of the HDPM can be calculated by(26)PCV1=τ1×θ˙1,PSV1=τ4×θ˙4,
where PCV1 is the CV motor power of the first group of the HDPM and PSV1 is the servomotor power of the first group of the HDPM. Similarly, the other motor power, PCV2, PSV2, PCV3, PSV3, PCV4, PSV4, of the HDPM can be derived, respectively.Figure11 shows motor absolute value power curves of four groups of the HDPMs of the HDCPM. From Figure 11, for each of the HDPM, the range of the CV motor power (PCVi) is larger than that of the servomotor power (PSVi). It is illustrated also that the CV motors undertake the majority of power for the HDCPM system.Figure 11
Motor absolute value power curves of the four groups of HDPMs of the HDCPM.As the motion circular trajectory of the end-effector is the same, and as the structural parameters of the classical approach for the conventional CPM are the same (cables are driven only by coils and servomotors), the servomotor power of the conventional CPM can be obtained by(27)Pi=Ti×l˙i(i=1,2,3,4),
where Pi is the ith servomotor power, Ti represents the ith cable tension, l˙i denotes the ith cable velocity.Figure12 shows the servomotor absolute value power curves of the classical approach for the conventional CPM. For comparison, absolute value power curves of the first group of HDPM of the HDCPM and the first servomotor of the conventional CPM are given in Figure 13. As it is seen, large range of servomotor power |P1| of the conventional CPM may be decreased by using the HDPM.Figure 12
Servomotor absolute value power curves of the conventional CPM.Figure 13
Comparison of absolute value power curves between the HDCPM and conventional CPM.In order to investigate the power property of the HDCPM further, the following simulation is performed as the CV motors of the four groups of the HDPMs of the HDCPM run with 1.1 rad/s (ω=1.1 rad/s). Figure 14 displays curves of the link angle and velocity of the HDPM. Angular curves of all links maintain good continuity, and the numerical values are reasonable. Figure 15 shows absolute value power curves of four groups of the HDPMs of the HDCPM. Figure 16 illustrates absolute value power curves of the first group of HDPM of the HDCPM and the first servomotor of the conventional CPM. From the simulations, we can see that the CV motors provide the majority of power supply and the servomotors mainly act as a motion regulation device as depicted in Figures 11, 13, 15, and 16.Curves of the links motion parameters of the HDPM asω=1.1 rad/s: (a) angle and (b) angular velocity.
(a)
(b)Figure 15
Motor absolute value power curves of the four groups of HDPMs of the HDCPM asω=1.1 rad/s.Figure 16
Comparison of absolute value power curves between the HDCPM asω=1.1 rad/s and conventional CPM.Case 2 (line segment trajectory).
The equation of the line in global coordinate system is(28)x=t6+0.25,y=t12+13,z=t6+0.25,
where 0<t<2(s). From (28), it can be noted that the center of gravity of the end-effector moves the geometric path that is line segment with the constant velocity v=0.25 m/s, the initial point is (0.25, 1/3, 0.25) (m). The trajectory tracking of the end-effector related to Case 2 is shown in Figure 17. From this figure the formulation tracks the planned trajectory relatively well.
Figure18 shows the changes of cables lengths. As can be seen, the variations of the curves of cables lengths are in accordance with the reference path (28) of the end-effector.
Figure19 depicts the curves of cable tensions for Case 2. One can see that the curves are smooth, which is the ideal data for the HDCPM system which operates steadily.
The angular variables of links of the HDPM for Case2 can be obtained by the inverse kinematics analysis of the HDCPM. The angle curves of the HDPM links are shown in Figure 20(a). The corresponding angular velocity and angular acceleration are shown in Figures 20(b) and 20(c). Figure 21 shows the running trajectory of the joint C in local x′o′z′ plane.
The resulting force curves of the joints of the HDPM for Case2 can be solved by the dynamics modeling of the HDCPM in Section 3, as shown in Figure 22.  It can be noted that the resulting force curves of the joints have the same trend and change smoothly. Comparison of absolute value power curves of the HDPMs in the HDCPM and the servomotors of the conventional CPM for Case 2 are shown in Figure 23. Referring to Figure 23, it can be seen that the CV motors provide the majority of power supply and the servomotors mainly act as a motion regulation device.
From the results and discussions, illustrative simulation studies highlight its performances, which lay a foundation for the further research on optimization and real-time control.Figure 17
Following trajectory of the line segment motion.Figure 18
Curves of cable length for Case2.Figure 19
Curves of cable tension for Case2.Curves of the links motion parameters of the HDPM for Case2: (a) angle, (b) angular velocity and (c) angular acceleration.
(a)
(b)
(c)Figure 21
Running trajectory of the jointC in x′-z′ plane for Case 2.Figure 22
Resultant force curves of the joints of the HDPM for Case2.Figure 23
Comparison of absolute value power curves between the HDCPM and conventional CPM for Case2.


## 5.2. Workspace Simulation

The workspace of jointC is determined by geometric constraints (17) given in Figure 24. The simulation results of the reachable 3D workspace for the HDCPM and the projections of workspace onto the X-Y plane and X-Z plane are shown in Figure 25. The whole workspace is approximately an upside-down cone. It should be noted that most of the workspace volume is concentrated in the upper part of the tower cube, which is important for trajectory planning of the end-effector.Figure 24
Workspace of jointC in x′-z′ plane.Workspace of the HDCPM (a) and projections of workspace ontoX-Y plane (b) and X-Z plane (c).
(a)
(b)
(c)


## 6. Concluding Remarks

This paper describes the theoretical investigations of a class of the completely restrained hybrid-driven based cable parallel manipulators (HDCPMs) with three translational motions. The HDCPM combines the advantages of four groups of HDPMs and a 4-cable parallel manipulator in a way to provide a solution for moving heavy objects with high efficiency and high performance. Detailed design and implementation of a virtual prototype model of the HDCPM is presented as well. Dynamic modeling of the full system including both the mechanical system and the driving system is addressed. The inverse kinematic and dynamic problems of the HDCPM system are resolved on condition that an operation path of the end-effector has been planned. Then, the main factors, which influence the workspace of the HDCPM such as workspace of the HDPM and the length of cables, are analyzed based on the structure constraints of the HDCPM. Finally, illustrative examples for trajectory tracking and workspace demonstrate that the manipulator design is reasonable and the associated models are correct.The prototype verification of the designed HDCPM is currently under investigation. Furthermore, more work is in progress to implement control system design, design modification, and optimum synthesis for the HDCPM. In addition, the methodology presented in this work could be used as a guide for the design and analysis of other manipulators.



---

*Source: 914653-2013-03-14.xml*