Next Article in Journal
A Distributed Automatic Video Annotation Platform
Previous Article in Journal
Multiblock Analysis Applied to TD-NMR of Butters and Related Products
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Empirical Quasi-Static and Inverse Kinematics of Cable-Driven Parallel Manipulators Including Presence of Sagging

Department of Mechatronics, HCMC University of Technology and Education; Ho Chi Minh 700000, Vietnam
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(15), 5318; https://doi.org/10.3390/app10155318
Submission received: 27 June 2020 / Revised: 20 July 2020 / Accepted: 24 July 2020 / Published: 31 July 2020
(This article belongs to the Section Mechanical Engineering)

Abstract

:
Cable-driven parallel manipulators (CDPMs) have been of great interest to researchers in recent years because they have many advantages compared to the traditional parallel robot. However, in many studies they lack the cable’s elasticity that leads to flexible cables just being considered as extendable rigid links. Furthermore, an external force acts on the extremities of cable and the self-weight is relevant to the length of it. Experimentally, a small change in length produces a huge change in tension act on the entire cable. By this property, the adjusting length of cable is often added to the traditional inverse kinematic solution in order to reduce the tension force exerted on the cable. This means that the load on the actuator is also reduced. Because of the relationship between tension that acts on the cable and its length, the kinematic problem itself does not make sense without concerning the static or dynamic problems. There is often interest in planning forces for actuators and the length of cables based on a given quasi-static trajectory of the moving platform. The mentioned problem is combined with the quasi-static problem with the inverse kinematic problem of CDPM. In this study, we introduce a novel procedure to produce the quasi-static model and inverse kinematic model for CDPM with the presence of sagging by using both an analytic approach and empirical approach. The produced model is time-efficient and is generalized for spatial CDPM. To illustrate the performance of the proposed model, the numerical and experimental approaches are employed to determine particular solutions in the feasible solutions set produced by our model in order to control the two redundant actuators’ CDPM tracking on a certain desired trajectory. Its results are clearly described in the experimental section.

1. Introduction

1.1. Motivation

Cable-driven parallel manipulators (CDPMs) are very similar to the traditional parallel robot since they also have a static frame, an end-effector, and many elements that connect the static frame and the end-effector. However, CDPM can easily variate its workspace, which is significantly greater than the other, because the flexible cables are used instead of rigid links. Except for those, they also have other great advantages such as economic structure, transportability, and easy assembling or disassembling [1]. With these outstanding advantages, applications of CDPM can be found in many fields: warehousing in logistics [2], contour crafting system in the construction industry [3], camera system in the stadium, etc. Many CDPMs are enhanced its load capacity also the acceleration by increasing the number of cables. Thus, lots of them have redundant actuators and this kind of structure produces an overconstrained system [4,5,6]. Furthermore, the cables that are used to link the end-effector and the robot fixed frame have considerable mass and elasticity, which is mentioned. These kinds of CDPM do not have a unique solution for the statistic or dynamic problem, and the kinematic problem cannot be solved alone. The major problem comes from the deformation of the cables: it increases the length of cables. In practice, these elongations cannot be measured and adjusted, resulting in deviations in the position and orientation of the end-effector. To the best of our knowledge in this field, it is merely impossible to directly extract exact solutions in terms of analysis. Therefore, in the kinematic or dynamic problem, with the given pose of the moving platform, the feasible length of cables and the corresponding tensions are obtained via numerical methods. However, by the complexity of the objective functions, these methods are almost impossible to use in real-time. In order to reduce the complexity of the dynamic problem while preserving the usability of the study, quasi-static problem will be considered in this paper. Similar to the dynamic problem, the quasi-static problem is time-dependent. However, the movement of the moving platform is very slow, or the load is very light, therefore the inertia effect is ignorable. That is very useful in many applications such as contour crafting systems and camera systems, since it simplifies the dynamic problem. In this study, we propose a method that produces the approximated model to solve the quasi-static problem. From that, the length of each cable in which the elongations and the adjusting length are also taken account can be determined in an instant.

1.2. Related Work

Because of its high load capacity and stability, the overconstrained CDPM is an interesting topic for many researchers. As a consequence of the redundant actuators, formulating the model in order to directly obtain the feasible distribution of force in real-time is a challenge. The interval analysis and Levenberg–Marquardt algorithm for real-time forward kinematics are deployed by Pott and Andreas in [7,8]. The same problem in planar CDPM has been studied by Reza Oftadeh et al. in [5] by using extra force sensors. Côté et al. employed quadratic programming method in [9] to handle the tension distribution in the spatial CDPM. Linear programming method is used in [6] to distribute the tension by Lamaury et al. A versatile tension distribution algorithm for two redundant actuators’ CDPM was proposed by Gouttefarde, Marc et al. in [10]. That algorithm is based on the properties of a feasible region related to the null space of wrench matrix and unidirectional searching technique. Borgstrom, Per Henrik et al. employed a linear programming approach [11] to determine minimum L1-norm solution in feasible polyhedron (associated with the null space of wrench matrix). These numerical methods are confirmed as usable in real-time and their results are acceptable. However, not mentioning the elasticity and sagging can cause the degeneration in motion of the large-scale CDPMs. The sagging of the cable can be categorized by the source causing it: elasticity and adjusting length. Elastodynamic problem was analyzed by Han Yuan et al. in [12] by considering the dynamic stiffness of the sagging cable. Nguyen, Dinh Quan et al. are honored to establish a simplified model for cable [13] in static analysis of large-dimension CDPM. The mass and elasticity of cable are taken into account in dual problems: static and kinematic of CDPMs mentioned in [14]. Riehl, Nicolas et al. also used numerical methods to obtain a feasible solution. The combination of interval analysis and the Newton–Raphson scheme, performed to solve the forward kinematics with the presence of sagging cable, is presented in [15]. Although the method proposed by Merlet et al. produces a good solution, it cannot be used in real-time. Another article that uses interval analysis to real-time solve forward kinematics for both elastic and non-elastic cables is shown in [16]. Meanwhile, the methods used in [16] do not deal with saggable cables and ignore the mass of them.

1.3. Contributions

In this study, the quasi-static model and inverse kinematic model of spatial CDPM are constructed by using both an analytic approach and empirical approach. The empirical approach is based on the popular sagging cable model introduced by Irvine [17]. Our Method is valid for the low speed and low acceleration motion of the moving platform for any kind of spatial CDPM. This is evidenced by results from experiments and simulations detailed in Section 7.

2. Cable-Driven Parallel Manipulators Model

2.1. Sagging Cable Model

In this study, the realistic sagging cable model introduced by Irvine [17] will be used as the standard model in order to derive others. Irvine’s sagging cable model is widely used from small-scale to macro applications such as cable-stayed bridge, improving the analysis of CDPM [18], ship anchors analysis, etc. because of the validation for elastic and deformable cable whose mass is considerable. Let consider the cable that has an unstrained length L and one extremity is fixed in the origin of frame { A } , illustrated in Figure 1. The other extremity B was influenced by external force T = [ T x T z ] T , and we can treat this force as the tension of the cable at free extremity. The reaction force at A is denoted as T = [ T x T z ] . Let ρ , E, A be the linear density, Young’s modulus, and unstrained cross-section area of the considering cable, respectively. The Cartesian coordinates of a certain point on the cable with respect to { A } at the Lagrangian coordinates s is given by Equations (1) and (2).
x A ( s ) = T x A s E A + | T x A | ρ g [ sinh 1 ( T z A ρ g ( L s ) T x A ) sinh 1 ( T z A ρ g L T x A ) ]
z A ( s ) = T z A s E A + ρ g E A ( s 2 2 L s ) + 1 ρ g [ ( T x A ) 2 + ( T z A ρ g ( L s ) ) 2 ( T x A ) 2 + ( T z A ρ g L ) 2 ]
The Cartesian coordinates of extremities of cable respect to { A } can be expressed as in (3) and (4).
A A = [ x A ( 0 ) z A ( 0 ) ] T
B A = [ x A ( L ) z A ( L ) ] T

2.2. Quasi-Static and Inverse Kinematic Problem of Cable-Driven Parallel Manipulators

The kinematic problem of CDPM that just describes the geometric relationship between a redefined pose of end-effector and the cables’ length. However, each cable has its length depend on the tension which acts on it, i.e., the kinematic model alone no longer makes sense. Many works apply the static or dynamic problem to it and solve them simultaneously [3,14,15,16,19]. Meanwhile, the quasi-static problem of CDPM concerns finding the proper forces of the moving platform and forces of its constraints given a certain motion of it where the inertial effects are negligible. Once again, with the relationship of tension of each cable and its length, the inverse kinematic problem is solved. Let { W } be the reference frame attached to the fixed of CDPM and { R } = ( P ( t ) , θ ( t ) ) , illustrated in Figure 2, be the moving platform frame attached on its body, where P = [ x P y P z P ] T , θ = [ α P β P γ P ] T are the coordinate of center of mass and the vector that represents Euler angle of the moving platform with respect to { W } . From now on, any vector without a superscript related to any frame is implied to respect to { W } . We also note that x = [ P θ ] T represents the generalized coordinate of the moving platform. If the moving platform is isolated from the system, according to the Newton–Euler’s equation, the resultant force F and torque τ exert on its center of mass can be obtained via Equation (5)
[ M 0 0 Γ ] x ¨ + [ M W 0 0 W Γ ] x ˙ + [ M g 0 ] = [ F τ ]
where M is simplified into the scalar matrix of mass due to the homogeneous density and symmetric geometry of the end-effector, Γ is the inertia tensor respect to center of mass of the end-effector expressed in { W } , g is the gravitational acceleration vector and W is the skew symmetric matrix of θ ˙ . On the other hand, the inertial forces of cables are trivial since the weights of these cables are negligible compared to the end-effector. Thus, if tensions of cables are treated as the external forces of the end-effector, then the resultant force F and torque τ also can be expressed as in (6) and (7)
F x = i T x i F y = i T y i F z = i T z i m g
τ x = i ( b z i T y i b y i T z i ) τ y = i ( b x i T z i b z i T x i ) τ z = i ( b y i T x i b x i T y i )
where b i = B i P is the vector demonstrates the position where the reaction force of tension T i = [ T x i T y i T z i ] is exerted on the end-effector with respect to { R } . b i , given by (8), depends on the rotation matrix R R W that denotes the rotation from { R } to { W } defined in (9).
b i = ( R R W ) T b i R
R R W = [ C α P C β P C α P S β P S γ P S α P C γ P C α P S β P C γ P + S α P S γ P S α P C β P S α P S β P S γ P + C α P C γ P S α P S β P C γ P C α P S γ P S β P C β P S γ P C β P C γ P ]
Naturally, in the ideal statistic condition, or mobile platform moving with speed and acceleration which are not too high, the ith cable profile that is fixed at its extremities A i and B i completely lies in the plane that contains these points. Similar to the sagging cable model section, let { A i } be the 2D-frame that defines the considered plane and illustrated in Figure 2. We can define a i and T i respect to { A i } as (10) and (11), respectively
a i A i = R A i W a i
T i A i = R A i W T i
where
R A i W = [ C θ i S θ i 0 0 0 1 ]
In the matrix (12), θ i can be given by (13). Vector a i in Equation (10) can be obtained by deploying Equation (14).
θ i = a tan 2 ( a y i , a x i )
a i = P + b i A i
According to (11), the equations in (6) and (7) can be expressed as (15) and (16).
F x = i T x i A i cos ( θ i ) F y = i T x i A i sin ( θ i ) F z = i T z i A i m g
τ x = i ( b y i T z i A i + b z i T x i A i sin ( θ i ) ) τ y = i ( b z i T x i A i cos ( θ i ) b x i T z i A i ) τ z = i ( b x i T x i A i sin ( θ i ) + b y i T x i A i cos ( θ i ) )
By changing the frame which the tension relates to, the dimension of distribution forces problem is reduced. From the relationship between the length of cable and a force exerted on it, the solution of both the quasi-static and inverse kinematic problem can be expressed as in (17).
[ { T x i A i T z i A i L i } i = 1 n ] T

3. Quasi-Static and Inverse Kinematic Characteristics of Single Cable

Since everything will be discussed in general frame { A } in this section, if any vector without a denotation of frame, then this vector respect to { A } . For the convenience, T = [ T x T z ] T ϒ denotes the tension acts on free extremity B = [ x B z B ] T , ς = [ ρ E A ] T denotes the properties of materials that are used in cable, S = [ T x T z L ] T Δ denotes the solution of our study. According to (1), (2) and (4), the sagging cable model of Irvine in [17] can be used to archive the position of the free extremity B of cable and can be expressed in (18).
I : Δ × , [ S ς ] T B
The main study of this section is try to figure out the feasible tensions and feasible unstrained lengths based on the desired position of the free extremity of cable B d = [ x B d z B d ] T with given ς that can be obtained by deploying (19).
f : × Δ , [ B d ς ] T S
However, the exact (19) cannot be directly derived via (18) by its natural complexity. In previous research, based on given B d and the constant material of cable, numerical approaches have been deployed to find out feasible solution S . However, most of these methods cannot be used in real-time because of time-consumption. The approximated models were also introduced in research [13,16,20]. In this work, the approximated model of (19) will be constructed based on (18) by applying both the empirical statistical analysis approach and static analysis approach. Generally, the cost function of our approximated model can be express as Equation (20).
ε = | x B d x B | + | z B d z B |
The optimum model is obtained if ε = 0 ,   B d . However, it is impossible for any approximated model that converges to the actual model, i.e., we always have ε > 0 . However, we restrict feasible solutions with the maximum error ε < ε max . By using Irvine model with fixed ς and L , a set of feasible tension can be found to satisfy a given desired position B d . The set of feasible tension of an instance is illustrated in Figure 3.
Statement 1.
Let  ψ be the feasible tension set with given ε max , L , B d and ς . According to Irvine model, with given L , T and ς there exists only one B , so (18) is the injective function. Similarly, (19) is also the injective function. Hence, if  ε max 0 then | ψ | 1 , i.e., the set has unique tension  T ψ .
Collecting all of ψ with L R + , we can get the feasible solution set Ψ of (19) with given ς and considered B d as the input parameter. Figure 4 illustrates a specific case of Ψ with d = L D and D is the L 2 -norm of B d .
Statement 2.
Let Ψ Π be the feasible solution set with given B d and ς . If ε max 0 then Π Δ is the affine plane that contains Ψ and its associated vector space that is the orthogonal complement to ϒ . Thus, the projection of Ψ onto ϒ is the affine line π ϒ . Moreover, for each given L value, if ε max 0 then there exists only one feasible tension. Therefore, Ψ becomes the spatial curve in Δ . π can be expressed as the relationship between T z and T x in (21)
( π ) : T z = a π T x + b π
where a π and b π denote the parameters that depend on B d and ς .
The quasi-static model is established when these parameters are formed in explicit function of B d . A certain case of π is obtained by using linear regression to fit the projected solution on ϒ and it is illustrated in Figure 5.
On the other hand, Li, Hui, et al. pointed out the nonlinear relationship between T z and T x in [21]. In case of CDPM, the conditions of ρ g L < w max and T x > T x min must be satisfied in order to maintain the stability of the system where w max , T x min are the maximum self-weight of cable and minimum tension along the X-axis. Moreover, the nonlinear segment of the relationship between T z and T x is located at the low part of T x . Hence, we can ignore the nonlinear part of it and linearize the relationship between T z and T x when ρ g L < w max , by employing (21). Moreover, on the simplification of cable model in [13], Nguyen, Dinh Quan et al. proved that if T E A , then the elongation of cable caused by elasticity has an ignorable impact on the relationship between T z and T x . This is also proved in [20]. In fact, the above assumption is always satisfied since cables used in CDPM often have very high stiffness to withstand the load of the mobile platform. Therefore, a π and b π now just depend on B d and ρ . In Section 4, we will figure out the quasi-static model of a single cable by combining an analytical approach and empirical approach.

4. Empirical Quasi-Static Model of Single cable

4.1. Analyze the Natural Limit of Tension Act on Free Extremity of Cable

Let m , α be the magnitude and direction of T respect to the X-axis, respectively. Let β be the angle between A B and the X-axis. The direction of T can be expressed as (22).
α = tan 1 ( T z T x )
Since the relationship between T z and T x can be considered as (21), Equation (22) can be rewritten as (23).
α = tan 1 ( a π + b π T x )
In reality, we know that T T where T is the reaction force at A illustrated in Figure 1 due to the self-weight of the cable. Since T and T tangent to the cable at extremities A and B , respectively, there is always the sagging in the cable. This explains the phenomenon of not being able to straightly stress the cable. However, if | T z | ρ g L , then we can ignore the mass of cable and T + T 0 that leads to α β . This can be expressed as (24).
lim | T z | α = β
According to (21) and (23), if | T z | and a π 0 , then | T x | that leads to α tan 1 ( a π ) . This process can be given by (25).
lim | T z | α = lim | T x | ( tan 1 ( a π + b π T x ) ) = tan 1 ( a π )
From (24) and (25), we can conclude that if (21) is considered as the relationship between T x and T z , then the slope a π just depends on β and be expressed as (26).
a π = tan ( β )

4.2. Empirical Relationship between b π and [ B d ρ ]

4.2.1. Relationship between b π and D

In order to comfortably obtain the complete relationship between b π and [ D β ρ ] , the relationship between b π and one of the mentioned parameters are considered first and that results in new parameters. These parameters now just depend on the rest and so we recursively repeat the procedure until the last parameter. In our case, the order of execution will be D , ρ and β . By variating D , fixing β and ρ , and using linear regression to establish (21), the relationship between b π and D can be guessed by forming the scatter plot in Figure 6.
According to the results illustrated in Figure 6, the relationship between b π and D can be expressed as in Equation (27). The parameter of this equation changes if ρ changes.
b π = h b ( ρ ) D

4.2.2. Relationship between h b and ρ

In Figure 7, the values of hb are plotted with given discrete ρ ( 0 , 2 ) and fixed β . It is easy to state that h b and ρ are linearly decreasing in the considered open set. Furthermore, when β is variated, the relationship between h b and ρ does not change, i.e., h b does not depend on β . Thus, the relationship between h b and ρ can be expressed as h b = k h ρ , where k h is a constant.
Based on our taken data and proposed procedure, we have k h 4.845 . This parameter can be approximately equal to negative half of gravitational acceleration. Following this assumption, the relationship between h b and ρ can be expressed as (28).
h b = ρ g 2
The obtained parameter h b based on our procedure coincides with the simplified cable model introduced in [20].

4.3. Relationship between T z and T x

In Statement 2, we state that T z = a π T x + b π with a π and b π depend on B d and ρ . According to Equations (26) and (27), the relationship between T z and T x can be expressed as Equation (29).
T z = tan ( β ) T x + h b D
Substituting the parameter (28) to (29), the affine line π can be expressed as the set (30).
π = { ( T x , T z ) : T z = tan ( β ) T x ρ g D 2 }
The second statement pointed out that the solution of the quasi-static problem of single cable is not unique and the particular solution just can be resolved when the cable unstrained length L is determined. In Section 5, we will try to establish the relationship between cable tension act on the extremity B and the unstrained length of cable by using both analytical approach and empirical approach.

5. Relationship between Cable Unstrained Length and Cable Tension

5.1. Simplify the Problem

In order to formulate the relationship between L and T , Ψ that is the set of all feasible solutions S of (19) with given desired position B d and ς , must be established. However, the construction of Ψ directly is not simple. Meanwhile, in Statement 2, we assumed that Ψ Π , as long as ρ g L < w max and T x > T x min . By successfully constructing the projection of Π onto the ϒ , π , the relationship between T x and T z was obtained. Hence, the relationship between L and T can be formed by constructing the relationship between L and any component of tension such as T x , T z , T , α instead of directly constructing Ψ .
Since the relationship between T z and T x is described as (21), if b π 0 then the minimum magnitude of tension T min that satisfies a certain point B d with 0 < β < 90 ° and D > 0 , is always greater than 0, i.e., the function that expresses the relationship between L and T has its vertical asymptote that can be estimated as (31)
T min = | b π | a π 2 + b π 2
where a π and b π are the parameters that define π . On the other hand, when T x or T z is considered instead of T , the vertical asymptotes of the function that expresses the relationship between T x or T z , and d are T x = 0 or T z = b π respectively. Hence, to ease the calculation and effortlessly identify the model, T x is chosen. Moreover, we have L = D + d and D is determined so the relationship between L and T x can be replaced by d and T x . Let Ψ be the new solution set that can be expressed as (32).
Ψ = { ( T x , d ) : d = d ( T x ) }
Theoretically, if the cable is unstretchable then the adjusting length d cannot be lesser than 0 and vice versa. Ψ of unstretchable cable and stretch cable with given certain B d and ς are illustrated in the Figure 8.
The negative part of d of the stretchable cable illustrated in Figure 8 is caused by the elasticity of the cable and it has the boundary that obeys Hooke’s law. In general, modelling Ψ in the case of unstretchable cable is easier, since it has d = 0 as its horizontal asymptote. Meanwhile, Ψ in the case of stretchable cable has an inclined asymptote.
To avoid misunderstanding, from now on, L denotes the unstrained length of stretchable cable. Let s and Δ l ( s ) be the Lagrangian coordinates of an unstrained stretchable cable and elongation of the cable under external force at s, respectively. Δ l ( s ) can be formulated as Equation (33)
Δ l ( s ) = Δ l x ( s ) + Δ l z ( s )
where Δ l x , Δ l z are the strained length of cable at s along the X and Z axis, respectively. They can be formed by Equations (34) and (35), derived from Hooke’s law.
Δ l x ( s ) = T x s E A
Δ l z ( s ) = T z s E A ρ g ( L s ) s E A
Substituting (34), (35) to (33), the relationship between the elongation of overall cable and tension can be formulated by (36)
Δ l = T L E A
where Δ l = Δ l ( L ) . In terms of T x , (36) can be written as (37).
Δ l = T x 2 + ( a π T x + b π ) 2 L E A
Considering a stretchable cable, let L e = D + d e + Δ l be the strained length of it where d e is the adjusted length of stretchable cable and satisfies a certain point, B d , with given tension T e . According to statement 1, we can state that there exists the unique tension, T , that satisfies the same desired position, B d , when unstretchable cables are concerned, with a given cable length that equals the strained length of stretchable cable, i.e., L e = D + d . Following that, we can obtain Equation (38).
d = d e + Δ l
By applying Equation (38) to Ψ associated with stretchable cable, the approximated relationship between d and T x can be archived. In Figure 9, the blue dot illustrated Ψ taken from an Irvine model with given certain B d when the elongation part is removed. The red dot illustrated Ψ associated with stretchable cable when the elongation part is removed by (38).
According to the result that is illustrated in Figure 9, the feasible solution set Ψ can be expressed as (39) instead of (32)
Ψ = { ( T x , d ) : d ( T x , B d ) = d e ( T x , B d ) Δ l ( T x , B d ) }
where ( T x , d ) satisfies a certain point B d when considering the unstretchable cable.

5.2. Relationship between d and T x

5.2.1. Select Model

Experimentally, the rational model (40) was chosen in order to estimate Ψ . The rational model (40) that has a constant numerator and 2nd order demonstrator, has an adjustable vertical asymptote based on a d , and horizontal asymptote d = 0. Furthermore, this model is the fittest and simplest that we tried our best to figure out.
d = a d T x 2 + b d T x + c d
Since the vertical asymptote of Ψ is T x = 0, by the properties of (40), we have b d = 0 and c d = 0 . Hence, the model is simplified to (41)
d = a d T x 2
where a d depends on the desired position B d and ρ that can be expressed as (42). The only constraint of a d is a d 0 since d 0 .
a d = a d ( D , β , ρ )
Similar to constructing π , Equation (42) will be constructed based on the model (41) when one variable is variated and the others are fixed. Hence, the recently constructed model may have parameters that depend on variables that have not been investigated. In our case, the order of investigating parameters will be D , β and ρ .

5.2.2. Relationship between a d and D

By fixing B d and ρ , the adjusted length d that corresponds to a certain T x is obtained by optimizing (20). The relationship of d and T x is modeled by using nonlinear regression with respect to (41). The result is illustrated in Figure 10 with diversity of D .
In reality, if D 0 , then B d A , which leads to T x 0 and d R + . According to (41), the consequence of D 0 is a d 0 . That is also proven in Figure 11. Hence, the model estimated that the relationship between a d and d cannot have terms that are independent of D . Furthermore, if D and d is a finite number, then L D , which leads to T x . According to (41), if D then a d . By its monotonic property, the function that expresses the relationship between a d and D cannot have extrema.
According to these analyzes, either the power model or exponential model can be used to estimate the relationship between D and a d . However, the 3rd order power model (43) is the fittest model to the collected data. In the next step, β will be variated with fixed ρ .
a d = q a ( β ) D 3
The relationship (43) is always valid with constraint a d 0 since D 0 .

5.2.3. Relationship between β and q a

According to (43), each certain β gives us one value of q a and each value of q a gives us one curve expressed as the relationship between a d and D . Those curves are illustrated in Figure 12. By variating β , we can plot the relationship between q a and β , which is illustrated in Figure 13.
According to the result in Figure 13, the periodic function is valid with the relationship between β and q a . To satisfy the constraint a d 0 , the equation that expresses this relationship should have absolute value. Experimentally, the model (44) was chosen in order to estimate the relationship between β and q a . Once again, p q is the parameter that depends on ρ .
q a = p q ( ρ ) cos 4 ( β )

5.2.4. Relationship between p q and ρ

Similar to the procedure in relationship between h b and ρ , the relationship between p q and ρ will be investigated discretely in ρ ( 0 , 2 ) . The Figure 14 and Figure 15 show the relationship (44) with different ρ , and the relationship between p q and ρ , respectively.
According to these results, the model (45) can be used to estimate this relationship, since it is the simplest and fittest model for our gathered data. Furthermore, it is also valid to the constraint a d 0 .
p q = k p ρ 2
However, we cannot exactly formulate k p due to the errors in the established model. Thus, k p will be treated as the constant. Based on the non-linear regression method, we have k p 4.385 .

5.2.5. Relationship between L and T x

Substituting (37) and (41) to (38), the adjusted length in stretchable cable de can be expressed as (46).
d e = a d T x 2 T x 2 ( a π T x + b π ) 2 L E A
Based on the relationship between a d and D, β , Equation (46) turns to (47).
d e = k p ρ 2 cos 4 ( β ) D 3 T x 2 T x 2 ( a π T x + b π ) 2 L E A
Since L = D + d e , the explicit function L in terms of T x is given by (48).
L = k p E A ρ 2 cos 4 ( β ) D 3 T x 2 E A D T x 2 ( E A + T x 2 ( a π T x + b π ) 2 )
The equation of (48) is the explicit function of the unstrained length L of stretchable cable and the tension along the X-axis, T x . We also need to pay attention to parameters D and β , presented in (48). These parameters represent the position of extremity B d . Hence, (48) can be treated as the inverse kinematic model of a single cable.

6. From Quasi-Statics to Inverse Kinematics of Cable-Driven Parallel Manipulators

The equations in (15), (16) describe the total force and total torque that act on the end-effector respect to { W } . By substituting Equations (21) to (15), (16), these equations are fully expressed by the geometric parameters and properties of material of the cable with variables in terms of T x A . It can be given as the system of Equation (49)
H t x + b = [ F τ ]
where H is given as the matrix (50), t x = [ T x 1 A T x 2 A T x n A ] T is the vector of tensions along X-axis with respect to { A } and b , expressed as (51), is the vector of free coefficients related to T x A .
H = [ cos ( θ 1 ) cos ( θ 2 ) cos ( θ n ) sin ( θ 1 ) sin ( θ 2 ) sin ( θ n ) a π 1 a π 2 a π n b y 1 a π 1 + b z 1 sin ( θ 1 ) b y 2 a π 2 + b z 2 sin ( θ 2 ) b y n a π n + b z n sin ( θ n ) b z 1 cos ( θ 1 ) b x 1 a π 1 b z 2 cos ( θ 2 ) b x 2 a π 2 b z n cos ( θ n ) b x n a π n b x 1 sin ( θ 1 ) + b y 1 cos ( θ 1 ) b x 2 sin ( θ 2 ) + b y 2 cos ( θ 2 ) b x n sin ( θ n ) + b y n cos ( θ n ) ]
b = [ 0 0 i = 1 n b π i m g i = 1 n b π i b y i i = 1 n b π i b x i 0 ] T
Equation (5) gives us the wrench of end-effector respect to the redefined motion of the end-effector. Meanwhile, the wrench of end-effector, which is calculated based on a set of tensions, is given by Equation (49). Since we only investigated the active segment of cable, the wrench involved in friction on the pulley can be ignored. Therefore, Equation (49) can be substituted to (5) and expressed as the quasi-static model of the CDPM (52).
( H T H ) 1 H T { [ M 0 0 Γ ] x ¨ + [ M W 0 0 W Γ ] x ˙ + [ M g 0 ] b } = t x
In CDPM, the inverse kinematic problem relates to figuring out the length of each cable according to the given motion of the moving platform. The traditional inverse kinematic model is only suitable for robots that have rigid links. Meanwhile, the CDPM links between end-effector and fixed frame are flexible and elastic, i.e., the length of its link is respect to the load distribution. This is evidence of having independently done the quasi-static problem before the solution of the inverse kinematic problem was determined by Equation (48). With the given motion of end-effector x ( t ) , the solution of this hybrid problem can be determined by applying the overall procedure illustrated in Figure 16.
The trajectory is fed to these models: the quasi-static models of CDPM are considered first. The final result is obtained by passing the results of the quasi-static model of CDPM block together with the trajectory of end-effector to the relationship between cable length and tension block. That block can be considered as a pseudo-inverse-kinematic model of CDPM. In next section, we will employ this model into a series of experiment in order to verify its capacities.

7. Experiments and Results

In this section, the series of experiment and its results will be divided into two parts: simulation and experiment. All experiments are applied to over-constrained system with two redundant actuators because of their high popularity and application. Other parameters of the robot are mentioned in Table 1. In order to confirm the availability of the proposed model for any spatial CDPM, we considered a simulation test with two CDPM. They have the same number of cables and related parameters as summarized in Table 2, but they have different a configuration, illustrated in Figure 17. From that, we can compare the performance of the proposed model in these configurations.

7.1. Simulation

In reality, there are many means of 8-cable CDPM wiring, so we tried to offer the two most common configurations. Both configurations mentioned in this section have symmetric geometry. Figure 17a depicts the configuration where all eight wires are attached to top of the workspace and the other depicts the configuration where four wires are attached to top of the workspace and the remaining are attached to the bottom of the workspace. These configurations can be found in many prototypes. The first one is famous thanks to the introduction of CoGiRo, and also considers the optimal geometry of two redundant actuators suspended CDPM [22]. IPAnema1, presented by Pott, Andreas et al. in [23], is the well-known prototype of the second configuration illustrated in Figure 17b. In this experiment, the moving platform of simulated CDPM is expected to track on the desired trajectory illustrated in Figure 18 with the given constraint [ z P α P β P γ P ] = [ 1.5 0 0 0 ] . The 5th order polynomial control equations are applied in order to generate this trajectory. Time-based velocities of the end-effector along the X-axis and Y-axis are illustrated in Figure 19.
Because the mentioned CDPMs in this section are overconstrained systems, there is a presence of the non-unique pseudo inverse of (50). To illustrate the performance of our model during tracking on the desired trajectory, a simple numerical approach is employed to solve that issue. The goal of this procedure is to find out a smooth series of minimum tensions which also satisfies the desired trajectory and constraints of our proposed model. Firstly, the particular solution t x 0 = [ T x 1 0 A T x 2 0 A T x 8 0 A ] T of (49) is a minimum L2 norm solution. The general solution of (49) can be expressed as (53)
t x = t x 0 + N λ
where N is the 8 × 2-matrix that represents the orthonormal basis for the null space of H and the particular solution t x of (49) depends on the value of λ = [ λ 1 λ 2 ] T . The constraint of (53) is introduced, and can be written as (54)
t x > t x min , T x i min = b π i a π i + T x i a d j , T x i a d j > 0
where t x min = [ T x 1 min A T x 2 min A T x 8 min A ] . In fact, the addition of T x i a d j A to maintain the accuracy of our proposed model as we discussed in Section 3. Applying the constraint (54) to Equation (53), the feasible solution set Θ λ is obtained and it is surrounded by a planar polygon. The polygonal boundary of Θ λ is defined by the vertex set Ω λ = { λ i } i = 1 h .However, the obtained feasible solution set Θ λ is presented in the domain of N . Since we cannot directly determine suitable t x in that space, Θ λ is transformed to new set Θ t presented in the null space of H by executing the linear affine transformation (53) to set Ω λ , resulting in a new vertex set Ω t = { t x i } i = 1 h that defines the R 8 -polygon. The kth particular solution t x k is selected based on t x k 1 . To produce the smooth series of the minimum L2-norm solution, we will select t x k by optimizing the objective function (55).
min t x k Θ t k t x k t x k 1
The results of this experiment, obtained by employing our proposed model and the numerical method mentioned above, are illustrated in Figure 20 and Figure 21.
In the results of both configurations illustrated in Figure 20 and Figure 21, our proposed model produces a continuous and smooth set of tension. In the second configuration, the T x of cables remaining in bottom of the workspace is small in order to minimize the load of set of cables staying in the top, but still valid, since T x > b π / a π . This produces uneven load distribution and explains the good load capacity of the first configuration, as it has the same load but the highest load among cables of the second configuration, approximated as twice the first configuration. However, the solution is still applicable in both cases. In this test, we also make a simple controller that takes the tension forces as the input and responses by unstrained length L to simulate the motion of CDPMs. Thereby, the positioning accuracy of the moving platform is considered and illustrated in Figure 22.
The error of center of mass of moving platform during tracking on the desired trajectory is formulated in (56).
ξ ( t ) = i B i ( t ) P ( t )
It can be seen that these errors are varying depending on the trajectory of moving platform and it does not tend to decrease over time. However, theoretically, our proposed model produces the new trajectory maintain positioning error around ( 6 , 6 ) mm.

7.2. Experiment

Figure 23 illustrates our prebuilt spatial CDPM that contains two redundant actuators. Its configuration is similar to the first configuration of the simulated system and its main parameters are shown in Table 1. Figure 24 describes the control system of our prototype CDPM which consists of four similar groups. Each group contains two sets of actuator/sensors controlled by 1 PLC Siemens S7-1200. Each cable housed inside a winch is driven by a 220 V Delta AC Servo Motor ECMA-C20604RS with capacity of 400 W. The unstrained length and the cable tension force of cable is measured via an encoder and a load cell fixed on the winch, respectively. The value measured on the encoder will provide feedback to the driver, while the other, measured on the loadcell, is amplified and transmitted to PC via PLC. The PC is powered by Intel® core™ I3-2100 processor 3.10 MHz. It is used to generate and control signal for four groups through our application. One more response of PC is to synchronize signals from four PLCs, hence, the wireless communication is employed.
Due to the limit of technique in measurement positioning accuracy, the experimental evaluation is done in the static state. By employing the proposed method in static state, the valid minimum tension sets in the whole size of robot is obtained. The set of the testing position is randomly taken from this. Thereby, the proposed model is applied and the robot is controlled to these points. Our robot is ordered to move point-to-point with the linear trajectory generated by the two-loop control system. The outer loop takes the cable tension forces and current pose of moving platform as its input, compares these values to the values of reference trajectory and aims at the response of the inner loop reference input. The inner loop is mainly based on the PID controller. It directly controls an actuator and feedback by the encoder. The linear trajectory is generated based on a 5th order trajectory generator with low speed. At each point, the robot is ordered to stop while measurement and error calculation is done.
Table 2 summarizes the test in this section, that includes testing desired positions, actual positions, and positioning errors obtained by applying Equation (56). Amplified load-cells and difference of cable lengths taken from encoders signals are also plotted in Figure 25 and Figure 26 while the control system is manipulating the test. The test was done in about 14 min with 100 samples/s. However, for the convenience in tracking the signals while the moving platform is ordered to move to difference points, the intervals for stopping and calculating were replaced by vertical dash-lines represented in Figure 25 and Figure 26. The test begins and finishes at home position P 0 = [ 0 0 1.5 ] T . The torque fluctuation at the time the moving platform starts to move always has the largest amplitude and decreases gradually until it is stabilized, resulting in a jerk at the beginning. This is due to the overshooting of the PID controller. However, this jerk does not last long (about 2 s from the beginning) and the moving platform moves reliably for the rest of given trajectory.
Figure 27 shows the position of the desired point of the moving platform in 3D. The size of each point indicates the positioning error. This error is produced by many sources: the accuracy of proposed model, the measurement method, the accuracy of measurers, non-ideal characteristics of cable material, latency sync, transmission error and the cumulative error. According to the result of the experiment and our assumptions, the relationship between the desired position and positioning accuracy of the moving platform can be interpreted easily. Cable tensions at low altitude are smaller than those at high altitude. Furthermore, if the moving platform approaches a pole of robot fixed frame, then the force distribution is very unbalanced among cables. These reasons are also illustrated in Figure 25. On the other hand, the larger the cable tensions, the smaller the positioning error. In Figure 27, P 2 , P 8 and P 6 have a low altitude, resulting in the increase in positioning error. Meanwhile, points which have a high altitude (e.g., P 7 , P 5 and P 4 ), have a better accuracy.

8. Conclusions

In this study, we introduced the model that use to solve the quasi-static problem of general spatial CDPM based on our simplified quasi-static model of a single cable. The explicit function that expresses the relationship between unstrained length L and tension exerted on free extremity B along the X-axis T x is also concerned and applied to solve the inverse kinematic problem of CDPM. The main motivation of employing our model is that it is applicable in real-time and has a low mean positioning error. However, we cannot control the accuracy of the position of the moving platform due to an undefined feasible wrench workspace. Furthermore, the relation (48) is valid for as long as T x > T x min . Thus, for a better version, we propose some future work:
  • Analyzing and determining the boundaries of the wrench-feasible workspace in order to determine the set of all wrenches that CDPMs can apply without violating the boundaries of the tensions and load capacities of actuators;
  • Due to the presence of unstrained cables and large inertia of the moving platform, the vibration is easily identifiable in CDPM. Vibration has a great impact on the dynamic behavior of the positioning accuracy and deviation on force distribution. Therefore, having a good analysis of vibration is a premise to developing a good dynamic model.

Author Contributions

Conceptualization, N.T.T.; methodology, N.T.T.; software, P.G.L.; validation, P.G.L.; formal analysis, P.G.L.; resources, N.T.T.; data curation, P.G.L.; writing—original draft preparation, P.G.L.; writing—review and editing, P.G.L.; visualization, P.G.L.; supervision, N.T.T.; project administration, N.T.T.; funding acquisition, N.T.T. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This research was supported by Ho Chi Minh City University of Technology and Education, Vietnam.

Conflicts of Interest

The authors declare no conflict of interest

References

  1. Vafaei, A.; Khosravi, M.A.; Taghirad, H.D. Modeling and control of cable driven parallel manipulators with elastic cables: Singular perturbation theory. In International Conference on Intelligent Robotics and Applications; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  2. Alias, C.; Nikolaev, I.; Magallanes, E.G.C.; Noche, B. An Overview of Warehousing Applications based on Cable Robot Technology in Logistics. In Proceedings of the 2018 IEEE International Conference on Service Operations and Logistics, and Informatics (SOLI), Singapore, 31 July–2 August 2018. [Google Scholar]
  3. Bosscher, P.; Williams, R.L., II; Bryson, L.S.; Castro-Lacouture, D. Cable-suspended robotic contour crafting system. Autom. Constr. 2007, 17, 45–55. [Google Scholar] [CrossRef]
  4. El-Ghazaly, G.; Gouttefarde, M.; Creuze, V. Adaptive terminal sliding mode control of a redundantly-actuated cable-driven parallel manipulator: CoGiRo. In Cable-Driven Parallel Robots; Springer: Cham, Switzerland, 2015; pp. 179–200. [Google Scholar]
  5. Oftadeh, R.; Aref, M.M.; Taghirad, H.D. Forward kinematic analysis of a planar cable driven redundant parallel manipulator using force sensors. In Proceedings of the 2010 IEEE/RSJ International Conference on Intelligent Robots and Systems, Taipei, Taiwan, 18–22 October 2010. [Google Scholar]
  6. Lamaury, J.; Gouttefarde, M.; Chemori, A.; Hervé, P.-E. Dual-space adaptive control of redundantly actuated cable-driven parallel robots. In Proceedings of the 2013 IEEE/RSJ International Conference on Intelligent Robots and Systems, Tokyo, Japan, 3–7 November 2013. [Google Scholar]
  7. Pott, A. An algorithm for real-time forward kinematics of cable-driven parallel robots. In Advances in Robot Kinematics: Motion in Man and Machine; Springer: Dordrecht, The Netherlands, 2010; pp. 529–538. [Google Scholar]
  8. Pott, A.; Schmidt, V. On the forward kinematics of cable-driven parallel robots. In Proceedings of the 2015 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Hamburg, Germany, 28 September–2 October 2015. [Google Scholar]
  9. Côté, A.F.; Cardou, P.; Gosselin, C. A tension distribution algorithm for cable-driven parallel robots operating beyond their wrench-feasible workspace. In Proceedings of the 2016 16th International Conference on Control, Automation and Systems (ICCAS), Gyeongju, Korea, 16–19 October 2016. [Google Scholar]
  10. Gouttefarde, M.; Lamaury, J.; Reichert, C.; Bruckmann, T. A versatile tension distribution algorithm for $ n $-DOF parallel robots driven by $ n+ 2$ cables. IEEE Trans. Robot. 2015, 31, 1444–1457. [Google Scholar] [CrossRef] [Green Version]
  11. Borgstrom, P.H.; Jordan, B.L.; Sukhatme, G.S.; Batalin, M.A.; Kaiser, W.J. Rapid computation of optimally safe tension distributions for parallel cable-driven robots. IEEE Trans. Robot. 2009, 25, 1271–1281. [Google Scholar] [CrossRef] [Green Version]
  12. Yuan, H.; Courteille, E.; Deblaise, D. Elastodynamic analysis of cable-driven parallel manipulators considering dynamic stiffness of sagging cables. In Proceedings of the 2014 IEEE International Conference on Robotics and Automation (ICRA), Hong Kong, China, 31 May–7 June 2014. [Google Scholar]
  13. Nguyen, D.Q.; Gouttefarde, M.; Company, O.; Pierrot, F. On the simplifications of cable model in static analysis of large-dimension cable-driven parallel robots. In Proceedings of the 2013 IEEE/RSJ International Conference on Intelligent Robots and Systems, Tokyo, Japan, 3–7 November 2013. [Google Scholar]
  14. Riehl, N.; Gouttefarde, M.; Krut, S.; Baradat, C.; Pierrot, F. Effects of non-negligible cable mass on the static behavior of large workspace cable-driven parallel mechanisms. In Proceedings of the 2009 IEEE International Conference on Robotics and Automation, Kobe, Japan, 12–17 May 2009. [Google Scholar]
  15. Merlet, J.-P.; Alexandre-dit-Sandretto, J. The forward kinematics of cable-driven parallel robots with sagging cables. In Cable-Driven Parallel Robots; Springer: Cham, Switzerland, 2015; pp. 3–15. [Google Scholar]
  16. Kozak, K.; Zhou, Q.; Wang, J. Static analysis of cable-driven manipulators with non-negligible cable mass. IEEE Trans. Robot. 2006, 22, 425–433. [Google Scholar] [CrossRef]
  17. Irvine, H.M. Cable Structures; The MIT Press: Cambridge, MA, USA, 1981; pp. 15–24. [Google Scholar]
  18. Merlet, J.-P. Some properties of the Irvine cable model and their use for the kinematic analysis of cable-driven parallel robots. Mech. Mach. Theory 2019, 135, 271–280. [Google Scholar] [CrossRef] [Green Version]
  19. Williams, R.L.; Gallina, P.; Rossi, A. Planar cable-direct-driven robots, part i: Kinematics and statics. In Proceedings of the 2001 ASME Design Technical Conference, 27th Design Automation Conference, Pittsburgh, PA, USA, 9–12 September 2001. [Google Scholar]
  20. Gouttefarde, M.; Collard, J.-F.; Riehl, N.; Baradat, C. Simplified static analysis of large-dimension parallel cable-driven robots. In Proceedings of the 2012 IEEE International Conference on Robotics and Automation, Saint Paul, MN, USA, 14–18 May 2012. [Google Scholar]
  21. Li, H.; Zhang, X.; Yao, R.; Sun, J.; Pan, G.; Zhu, W. Optimal force distribution based on slack rope model in the incompletely constrained cable-driven parallel mechanism of FAST telescope. In Cable-Driven Parallel Robots; Springer: Berlin/Heidelberg, Germany, 2013; pp. 87–102. [Google Scholar]
  22. Gouttefarde, M.; Collard, J.-F.; Riehl, N.; Baradat, C. Geometry selection of a redundantly actuated cable-suspended parallel robot. IEEE Trans. Robot. 2015, 31, 501–510. [Google Scholar] [CrossRef] [Green Version]
  23. Pott, A.; Mütherich, H.; Kraus, W.; Schmidt, V.; Miermeister, P.; Dietz, T.; Verl, A. IPAnema: A family of cable-driven parallel robots for industrial applications. In Cable-Driven Parallel Robots; Springer: Berlin/Heidelberg, Germany, 2013; pp. 119–134. [Google Scholar]
Figure 1. Cable sag model in frame { A } .
Figure 1. Cable sag model in frame { A } .
Applsci 10 05318 g001
Figure 2. General geometric model of cable-driven parallel manipulators (CDPM).
Figure 2. General geometric model of cable-driven parallel manipulators (CDPM).
Applsci 10 05318 g002
Figure 3. The feasible tensions T with given L , B d and ς with ε max > 0 .
Figure 3. The feasible tensions T with given L , B d and ς with ε max > 0 .
Applsci 10 05318 g003
Figure 4. The feasible solutions S Ψ and ψ is the blue segment line with ε max > 0 .
Figure 4. The feasible solutions S Ψ and ψ is the blue segment line with ε max > 0 .
Applsci 10 05318 g004
Figure 5. The projection of Ψ onto ϒ is the affine line π with given B d and ς .
Figure 5. The projection of Ψ onto ϒ is the affine line π with given B d and ς .
Applsci 10 05318 g005
Figure 6. Relationship between b π and D with different ρ .
Figure 6. Relationship between b π and D with different ρ .
Applsci 10 05318 g006
Figure 7. The relationship between β and h b with given ς .
Figure 7. The relationship between β and h b with given ς .
Applsci 10 05318 g007
Figure 8. The relationship between d and T x of both of unstretchable and stretchable cable in the instance.
Figure 8. The relationship between d and T x of both of unstretchable and stretchable cable in the instance.
Applsci 10 05318 g008
Figure 9. The relationship between d and m that concerns to the unstretchable cable and the approximated one.
Figure 9. The relationship between d and m that concerns to the unstretchable cable and the approximated one.
Applsci 10 05318 g009
Figure 10. The relationship between Tx and d with different D and fixed β = 0.754   rad .
Figure 10. The relationship between Tx and d with different D and fixed β = 0.754   rad .
Applsci 10 05318 g010
Figure 11. The relationship between D and p.
Figure 11. The relationship between D and p.
Applsci 10 05318 g011
Figure 12. The relationship between D and a d with different β .
Figure 12. The relationship between D and a d with different β .
Applsci 10 05318 g012
Figure 13. The relationship between β and q a .
Figure 13. The relationship between β and q a .
Applsci 10 05318 g013
Figure 14. The relationship between β and q a with different ρ .
Figure 14. The relationship between β and q a with different ρ .
Applsci 10 05318 g014
Figure 15. The relationship between p q and ρ .
Figure 15. The relationship between p q and ρ .
Applsci 10 05318 g015
Figure 16. The proposed method used to obtain (17) based on given motion of end-effector.
Figure 16. The proposed method used to obtain (17) based on given motion of end-effector.
Applsci 10 05318 g016
Figure 17. (a) The first configuration; (b) the second configuration used in simulation test.
Figure 17. (a) The first configuration; (b) the second configuration used in simulation test.
Applsci 10 05318 g017
Figure 18. The trajectory of moving platform for simulation test.
Figure 18. The trajectory of moving platform for simulation test.
Applsci 10 05318 g018
Figure 19. The velocity of moving platform in desired trajectory.
Figure 19. The velocity of moving platform in desired trajectory.
Applsci 10 05318 g019
Figure 20. Unstrained length L , adjusting length d , tension along X-axis T x and tension along Z-axis T z respect to { A } of CDPM that has the first configuration during the simulation test.
Figure 20. Unstrained length L , adjusting length d , tension along X-axis T x and tension along Z-axis T z respect to { A } of CDPM that has the first configuration during the simulation test.
Applsci 10 05318 g020
Figure 21. Unstrained length L , adjusting length d , tension along X-axis T x and tension along Z-axis T z respect to { A } of CDPM that has the second configuration during the simulation test.
Figure 21. Unstrained length L , adjusting length d , tension along X-axis T x and tension along Z-axis T z respect to { A } of CDPM that has the second configuration during the simulation test.
Applsci 10 05318 g021
Figure 22. The positioning error of moving platform during the simulation test of both configuration (left: first configuration, right: second configuration).
Figure 22. The positioning error of moving platform during the simulation test of both configuration (left: first configuration, right: second configuration).
Applsci 10 05318 g022
Figure 23. Our prebuilt spatial CDPM with 2 redundant actuators used in experiment (a) perspective view; (b) end-effector.
Figure 23. Our prebuilt spatial CDPM with 2 redundant actuators used in experiment (a) perspective view; (b) end-effector.
Applsci 10 05318 g023
Figure 24. Diagram of the control system of our proposed prototype.
Figure 24. Diagram of the control system of our proposed prototype.
Applsci 10 05318 g024
Figure 25. Amplified signals are taken from load-cells while the moving platform moves point-to-point with the linear trajectory.
Figure 25. Amplified signals are taken from load-cells while the moving platform moves point-to-point with the linear trajectory.
Applsci 10 05318 g025
Figure 26. Different cable lengths are obtained by integrating the encoder signal while the moving platform moves point-to-point with the linear trajectory.
Figure 26. Different cable lengths are obtained by integrating the encoder signal while the moving platform moves point-to-point with the linear trajectory.
Applsci 10 05318 g026
Figure 27. Set of static poses and the positioning error of the moving platform illustrated by the size of each point.
Figure 27. Set of static poses and the positioning error of the moving platform illustrated by the size of each point.
Applsci 10 05318 g027
Table 1. Parameters of CDPMs used in experiment and simulation.
Table 1. Parameters of CDPMs used in experiment and simulation.
ParametersExperimentSimulation
Fixed frame size3 m × 6 m × 3 m15 m × 11 m × 6 m
End-effector size0.3 m × 0.3 m × 0.5 m0.3 m × 0.3 m × 0.5 m
End-effector load50 kg50 kg
ρ 0.09 kg/m0.067 kg/m
E 3.5 GPa2 GPa
A 4 mm22 mm2
Table 2. Performance of model in our CDPM based on set of static poses.
Table 2. Performance of model in our CDPM based on set of static poses.
Desired Positions (m)Actual Position (m)Positioning Error (mm)
P1(0.22, −0.45, 1.43)(0.2220, −0.4510, 1.4470)17.1464
P2(−0.13, 1.42, 0.68)(−0.1290, 1.4380, 0.7210)44.7884
P3(−0.29, −0.25, 1.52)(−0.2920, −0.2510, 1.5340)14.1774
P4(0.30, 1.26, 2.31)(0.3000, 1.2610, 2.3150)5.0990
P5(−0.47, 0.80, 2.17)(−0.4710, 0.8000, 2.1760)6.0828
P6(0.48, −0.85, 0.72)(0.4860, −0.8550, 0.7460)27.1477
P7(−1.08, −0.75, 2.54)(−1.0810, −0.7510, 2.5410)1.7321
P8(0.26, 1.06, 0.44)(0.2640, 1.0760, 0.4790)42.3438
P9(0.45, −1.73, 1.71)(0.4500, −1.7340, 1.7270)17.4642
P10(0.67, 0.52, 2.02)(0.6710, 0.5200, 2.0280)8.0623

Share and Cite

MDPI and ACS Style

Gia Luan, P.; Truong Thinh, N. Empirical Quasi-Static and Inverse Kinematics of Cable-Driven Parallel Manipulators Including Presence of Sagging. Appl. Sci. 2020, 10, 5318. https://doi.org/10.3390/app10155318

AMA Style

Gia Luan P, Truong Thinh N. Empirical Quasi-Static and Inverse Kinematics of Cable-Driven Parallel Manipulators Including Presence of Sagging. Applied Sciences. 2020; 10(15):5318. https://doi.org/10.3390/app10155318

Chicago/Turabian Style

Gia Luan, Phan, and Nguyen Truong Thinh. 2020. "Empirical Quasi-Static and Inverse Kinematics of Cable-Driven Parallel Manipulators Including Presence of Sagging" Applied Sciences 10, no. 15: 5318. https://doi.org/10.3390/app10155318

APA Style

Gia Luan, P., & Truong Thinh, N. (2020). Empirical Quasi-Static and Inverse Kinematics of Cable-Driven Parallel Manipulators Including Presence of Sagging. Applied Sciences, 10(15), 5318. https://doi.org/10.3390/app10155318

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop