Wrench-Closure Condition of Cable-Driven Parallel Manipulators

: Cable-driven parallel manipulators (CDPM) have parallel structures that consist of moving platforms connected to the ﬁxed platform through many ﬂexible cables. The moving platform is driven by many winches, and because of the unidirectional property of cables (cables can only pull and not push in the moving platform), some speciﬁc workspaces of CDPM are limited and often do not exist (Gouttefarde and Gosselin 2006). Therefore, determining workspaces for CDPM become an important task, in order to easily plan the trajectory or control the robot. In this paper, we are interested in a set of poses of moving platforms, in which CDPM is always able to generate wrenches that balance any given external wrench exerted on the moving platform. This set of poses is also called wrench-closure workspace (WCW). In this work, we propose a novel procedure used to determine whether a pose of moving platform belongs to WCW. The condition used to check the feasibility of a certain pose in WCW is also called wrench closure condition (WCC). The proposed WCC is able to apply to completely or redundantly restrained CDPM. By analyzing the geometric properties of CDPM and applying the method used to check the feasibility of a system of inequalities, the algorithm used to check the presence of a given pose in WCW is established. To verify the performance of the proposed WCC, the test has been done in two different CDPMs that are clearly described in the introduction section.


Motivation
The cable-driven parallel manipulators (CDPMs) have a similar structure to the traditional parallel kinematic machines that consist of a moving platform (end-effector), the robot's fixed platform, and many links used to connect the moving platform to the fixed platform. The main difference between them is that, in the CDPMs, the rigid links are replaced by the flexible light-weight cables, and their lengths are adjusted by winches that are usually mounted on the fixed platform. By changing the length of these cables, the position and orientation of the moving platform are controlled. In recent decades, CDPMs are of great interest to researchers, because they overcome many disadvantages of conventional parallel mechanisms, such as easy assembling and disassembling, simple configuration, high-load capacity, etc. By these dominant advantages, CDPMs are widely applicable in many fields: contour crafting system in construction [1,2], air vehicle simulation [3], etc. The giant CDPM developed by Salafian and Iman was applied in agriculture [4]. Aerial robotic cameras used for broadcasting are studied in [5]. Besides the advantages that it brings, the use of cables instead of rigid links also presents many issues for CDPM. In order to deal with a trajectory planning task for the moving platform while balancing an external wrench exerted at the moving platform, the cable must be in stretching, i.e., there must always be cable tension forces exerted in them. The relationship between the elongation of cable and cable tension forces was also coordinated in our previous work [6]. Due to the unilateral property of the cable tension force, the cable can only be pulled but not be pushed, resulting in the wrench matrix of CDPM perhaps being in full rank, but not providing any solution. This leads to the limitation of a workspace of this kind of mechanism. In practice, many redundant actuators are added, not only for increasing load capacity, but also for reducing the limitation of workspace. The cable failure in CDPM is also researched in many articles [7,8], since this affects the limitation of workspace. In this study, we are interested in a set of poses of the moving platform in which CDPM can balance any external wrench by stretching the cables through winches. This set of poses is also called wrench-closure workspace (WCW). Details of other important workspaces of CDPM are clearly introduced in [9]. According to Gouttefarde and Gosselin (2006), WCW depends on the geometry of the mechanism only [10]. Hence, WCW is the prerequisite for an optimal design configuration and checking the validation for a solution of the trajectory planning task of CDPM. Therefore, WCW is the most important workspace and has been widely studied. Currently, the determination of WCW can be classified into two groups: the determination of the boundary of WCW and determination based on wrench closure condition (WCC). The first group directly sketches the WCW by using an analytical approach, while the second must be combined with some search algorithms to complete the same task. However, in this paper, the method that determines whether a pose of the moving platform belongs to WCW is of interest.

Related Works
Due to its importance, WCW has been studied and analyzed in many articles to determine the active area of the CDPM as well as its configuration to achieve the maximum operational space. WCW of planar parallel cable-driven mechanisms was analyzed by Gouttefarde and Gosselin (2006). They proposed many theorems to bound the wrench closure workspace in the reachable workspace depended on number of cables used in given mechanisms. According to Gouttefarde and Gosselin (2005), the procedure to disclose the WCW of spatial CDPM with seven cables in use is introduced [11]. Gouttefarde, Merlet and Daney (2006) once again introduce the theorem to determine the boundary of wrench-closure workspace of general spatial overconstraint CDPM [12]. A hybrid analytical-numerical approach to determine WCW of CDPM was established by Lau, Oetomo and Halgamuge (2011) in [13]. However, this method can be applied to a single redundant actuator CDPM. Li, etc. proposed the method of algebraic expression of WCW boundary for planar CDPM based on antipodal theorem [14]. This method gives more exact results but it is difficult to apply for spatial overconstrained CDPM. A close form-based approach to determine WCW for CDPM is introduced by Pott and Kraus in [15]. This approach is verified in the least time-consuming way, but it is only tested in a single redundant actuator CDPM. A recursive reducing space dimension algorithm in order to check the feasibility of a given pose of moving platform was proposed in [16]. However, this algorithm is highly time-consuming and was not given mathematic proof. One more article related to determining WCW of single redundant actuator CDPM is introduced by Diao and Ma in [17]. Their approach is also derived based on convex analysis and has a good performance.

Contributions
In this study, we introduce the procedure to check whether a certain pose of moving platform belongs to WCW. This procedure is based on the theory of convex analysis and feasibility determination approach for an inequalities system proposed by Dines in [18]. Diao and Ma also introduced the wrench closure condition to check the presence of a pose of moving platform in WCW [17]. However, that method is only applicable to one redundant actuator CDPM. Meanwhile, our method can apply for a completely or redundantly restrained CDPM. In the following sections, in order to verify the performance of the proposed method and give examples, the details of used CDPM configuration in this paper were illustrated in Figure 1. Determining WCW is a part of the design phase, when the accuracy of the model is unnecessary, so that the cable mass and cable elasticity are neglected in our study. redundantly restrained CDPM. In the following sections, in order to verify the performance of the proposed method and give examples, the details of used CDPM configuration in this paper were illustrated in Figure 1. Determining WCW is a part of the design phase, when the accuracy of the model is unnecessary, so that the cable mass and cable elasticity are neglected in our study.

Cable-driven Parallel Manipulators Kinematics
Let's consider the general structure of CDPM illustrated in Figure 2 that contains a moving platform connected with a fixed platform through several flexible cables. The reference frame {W} is attached to the fixed platform and the moving frame {R} is attached to the center of mass of moving platform. For convenience in notation, the rightsuperscript denotes the frame where the mentioned object reference to, and any vector without a right-superscript is implied to reference to {W}.
where W R R is the rotational matrix that represents the orientational of {R} respect to {W} and it can be defined by (2).
P P P P P P P P P P P P P P P P P P P P P P P P P P P P P R (2) where Cx and Sx represent for cos(x) and sin(x), respectively.

Cable-Driven Parallel Manipulators Kinematics
Let's consider the general structure of CDPM illustrated in Figure 2 that contains a moving platform connected with a fixed platform through several flexible cables. The reference frame {W} is attached to the fixed platform and the moving frame {R} is attached to the center of mass of moving platform. For convenience in notation, the right-superscript denotes the frame where the mentioned object reference to, and any vector without a right-superscript is implied to reference to {W}. The general coordinate of {R} with respect to {W} is defined by the vector x = [P θ] T where P = [x p y p z p ] T and θ = [α p β p γ p ] T are the positioning vector of center of mass and rotational vector that represents Euler angle (z-y-x) of the moving platform, respectively. Let A i and B i be the positioning vectors that represent the end of the ith cable at the pole and at the moving platform, respectively. According to Figure 2, we have a i = A i B i and b i = PB i . If the general coordinate of the moving platform is given, then the direction of tension exerted on any cable can be obtained through Equation (1) where W R R is the rotational matrix that represents the orientational of {R} respect to {W} and it can be defined by (2).
where C x and S x represent for cos(x) and sin(x), respectively. In order to balance the external wrench ψ e that exerted on P, the actuators mounted in the fixed frame are controlled to create cable tension forces. Moreover, the cable tension forces also affect cable lengths, while cable lengths are desired results in the kinematic problem of CDPM. Hence, the derivation of the relationship between ψ e and cable tension forces is an interesting issue. Let m and n be the degree of freedom of moving platform and the number of cables used in system, respectively. Considering the static equilibrium system of the moving platform, the relationship between the cable tension force and external wrench is established in Equation (3).
where W ∈ R mxn is the wrench matrix that has ith column vector is defined in (4), and t = [t 1 . . . t 2 ] is the vector that represents the set of magnitude of cable tension forces with t i ≥ 0, ∀i = {0, . . . , n}. In order to balance the external wrench ψe that exerted on P, the actuators mounted in the fixed frame are controlled to create cable tension forces. Moreover, the cable tension forces also affect cable lengths, while cable lengths are desired results in the kinematic problem of CDPM. Hence, the derivation of the relationship between ψe and cable tension forces is an interesting issue. Let m and n be the degree of freedom of moving platform and the number of cables used in system, respectively. Considering the static equilibrium system of the moving platform, the relationship between the cable tension force and external wrench is established in Equation (3). e = Wt ψ (3) where W∈R mxn is the wrench matrix that has ith column vector is defined in (4), and t = [t1 … t2] is the vector that represents the set of magnitude of cable tension forces with ti ≥

The Properties of Wrench-Closure Workspace
In the following section, for any vector v ∈ R n and scalar a, we denote that v = a means that vi = a, ∀i = {0…,n} and the other binary conditional operators are analogous to that. According to Equation (3), let ω = {w1,…,wn} be the set of column vectors of given wrench matrix; the external wrench ψe can be considered as a non-negative combination of ω since t ≥ 0. Let Sψ be the set of all nonnegative combinations of ω and can be expressed as (5).
where pos(ω) denotes the nonnegative linear span of ω. Again, on the problem of kinematics, the wrench matrix W depends on the pose x of the moving platform. Hence, Sψ also depends on x. On the other hand, in this paper, we

The Properties of Wrench-Closure Workspace
In the following section, for any vector v ∈ R n and scalar a, we denote that v = a means that v i = a, ∀i = {0 . . . , n} and the other binary conditional operators are analogous to that.
According to Equation (3), let ω = {w 1 , . . . , w n } be the set of column vectors of given wrench matrix; the external wrench ψ e can be considered as a non-negative combination of ω since t ≥ 0. Let S ψ be the set of all nonnegative combinations of ω and can be expressed as (5).
where pos(ω) denotes the nonnegative linear span of ω. Again, on the problem of kinematics, the wrench matrix W depends on the pose x of the moving platform. Hence, S ψ also depends on x. On the other hand, in this paper, we are interested in the set of poses, in which there always exists at least one set of non-negative tension forces balance arbitrary, given external wrench ψ e . The set of these poses is also known as the wrench-closure workspace, and this can be expressed as (6).
We can express (6) as the condition for pose x to be in WCW by (7).
In convexity analysis, if ω non-negatively spans R m , then the spanning set ω must hold at least m + 1 vectors, and if ω just has m + 1 vectors, then they should be affinely independent.
According to the results from Theorem 3.6 of Davis [19] and Theorem 2.3 of Conn, Scheinberg and Vicente [20], we can extract the condition (7) into two simpler conditions, which is stated in Theorem 1.

Theorem 1.
A general coordinate x of moving platform belongs to the wrench-closure workspace of a certain CDPM that must satisfy these following conditions: The Equation (8) is the necessary condition and it is equivalent to the wrench matrix W has full-rank. If dim(W) < m then the system (3) has no solution, even if t is arbitrarily chosen in R n . On the contrary, if the condition (8) is satisfied, then ω(x) non-negatively linear span the set which always contains a simplex cone. The condition (9) is sufficient for the set ω(x) to nonnegatively linear span R m . The Equation (9) is equivalent to Equation (10), in which the geometrical interpretation is easily archived.
where m-first column vectors in W are linearly independent. For convenience, let ω A = {w 1 , . . . , w m } be the set of m-first column vectors in W and ω B = {w m+1 , . . . , w n } be the set of remaining vectors. Geometrically, the condition (10) can be interpreted that the intersection between interior of cones spanned by ω B and -ω A must exist.
For instance, let us consider a 3 × 5 wrench matrix W = [w 1 w 2 . . . w 5 ] whose column vectors are divided into two sets, ω A = {w 1 , w 2 , w 3 } and ω B = {w 4 , w 5 }. These sets and their positive spanning sets are illustrated in Figure 3. Since ω A is a linearly independent set (w 1 , w 2 and w 3 do not lie in a common plane), the Equation (8) is satisfied. Secondly, the intersection between pos(−ω A ) and pos(ω B ) is not empty, hence, the Equation (9) is also satisfied. Therefore, the given wrench matrix that has column vectors illustrated in Figure 3 has a general coordinate presented in the wrench-closure workspace.  The special case of this problem is CDPM contains only one redundant actuator, i.e., n − m = 1. In that case, the condition (10) is simplified to (11).
where λ = [λ1 λ2 … λm] T and λi = ti/tn. Indeed, vector wn non-negatively spans a ray and this must belong to a cone spanned by −ωA to satisfy the Equation (10).
where λ = [λ 1 λ 2 . . . λ m ] T and λ i = t i /t n . Indeed, vector wn non-negatively spans a ray and this must belong to a cone spanned by −ω A to satisfy the Equation (10).
For instance, let us consider a 3 × 4 wrench matrix W = [w 1 w 2 w 3 w 4 ] which has column vectors illustrated in Figure 4. Similar to the last example, since w 1 , w 2 and w 3 are linearly independent, we have that span({w 1 , w 2 , w 3 , w 4 }) = R 3 . Moreover, w 4 belongs to the negative spanning set of w 1 , w 2 and w 3 then pos({w 1 , w 2 , w 3 , w 4 }) = R 3 .
Indeed, vector wn non-negatively spans a ray and this must belong to a cone spanned by −ωA to satisfy the Equation (10).

Simplify Problem by Using Linear Equivalence
Let M ∈ R mxm be the square matrix established by taking m-column vectors in W and H ∈ R mx(n-m) consist of the remaining column vectors in W. According to Theorem 1, if x ∈ WCW ⊂ R m then ω must contain m linearly independent vectors, i.e., there exist a certain set of m-column vectors in W that form the nonsingular matrix M.

Simplify Problem by Using Linear Equivalence
Let M ∈ R mxm be the square matrix established by taking m-column vectors in W and H ∈ R mx(n−m) consist of the remaining column vectors in W. According to Theorem 1, if x ∈ WCW ⊂ R m then ω must contain m linearly independent vectors, i.e., there exist a certain set of m-column vectors in W that form the nonsingular matrix M. Theorem 2. Given that N ∈ R mxm is the non-singular matrix, if ω non-negatively spans R m then N.ω = {Nw 1 , Nw 2 , . . . , Nw n } also non-negatively spans R m . If ω does not nonnegatively span R m , then N.ω also does not non-negatively span R m .
According to Theorem 2, we can state that if ω does or does not non-negatively span R m than M −1 ω also does or does not nonnegatively span R m . Thus, Equation (3) can be replaced by any of the Equations (12)- (14), but the condition represented by (9) is not influent.
where t A = [t 1 t 2 . . .  Figure 5 illustrated the specific case of 2 redundant actuators-planar CDPM when a general coordinate is presented in wrench-closure workspace and linear equivalence was applied to the wrench matrix generated by the considered general coordinate.
Based on Equation (12), we can write M −1 ω = {e1, e2, …, em, M −1 wm+1, …, M −1 wn}, where  {e1, e2, …, em} is the standard basis of R m . Following the Equation (10), the set of m-first column vectors in W is now represented as ωA = {e1, e2, …, em} and the set of remaining vectors is ωB = {M −1 wm+1, …, M −1 wn}. Geometrically, ωA non-negatively spans closed nonnegative hyperoctant of R m , R m +. Hence, if x belongs to WCW, then the intersection between the interior of cone ωB and open negative hyperoctant R mshould not be empty. For instance, Figure 5 illustrated the specific case of 2 redundant actuators-planar CDPM when a general coordinate is presented in wrench-closure workspace and linear equivalence was applied to the wrench matrix generated by the considered general coordinate.  In the special case, when n − m = 1, the sufficient condition for the general coordinate belonging to WCW is simplified as (15).
In the general case, when n − m > 1, we cannot determine whether x does or does not belong to WCW only based on the signs of elements in ω B . However, the linear equivalence nonnegative span set of ω simplifies the condition (10). Since ω A already non-negatively spans R m + , in Equation (14), it is equivalent to It A > 0, ∀ t A > 0, the sufficient condition for the general coordinate x belongs to WCW is defined by (16).
For a more compact system of inequalities, we can rewrite (16) as (17).
Theorem 3. A general coordinate x belongs to the wrench-closure workspace of a certain CDPM if (8) and (17) are satisfied.
In the next section, we will introduce the procedure in order to determine whether (17) does or does not have solution.

Checking the Feasibility of a System of Inequalities
Based on the matrix theory for solving the system of linear equations, L. Dines developed his own theory on determining the general solution of system of linear inequations. He also proposed the method used to check the feasibility of them. In this section, we will introduce L. Dines' theory, make our conclusions and establish our algorithm to solve the system (17).
Let's consider the system of inequations in form of Ax > 0 where A is a coefficient matrix. L. Dines proposed the concepts I-positive, I-negative and I-definite to establish his own work on checking the feasibility of Ax > 0. These concepts can be summarized by the Definition 1. However, conclusion 1 is not sufficient to check the feasibility of system Ax > 0, i.e., If the matrix A is not I-definite, then we cannot conclude that Ax > 0 is feasible or not. In this case, L. Dines introduced the concepts of I-complement and I-minors. The I-complement of the rth column of matrix A, A 1 r , is the matrix that can be established as follows.
Then, we orderly top-down establish each row of A 1 r by orderly combining a ir and a jr , and used them to compute each element of the corresponding row by the secondorder determinant (19), where value of k corresponding to the order of elements in the considered row.
a ir a ik a jr a jk , k = {1, 2 . . . , m}\r The last rows of A 1 r are derived from the rows of A that contain a kr and formed by taking these rows without elements a kr .
By that procedure, the matrix A 1 r which has P = N + Z rows and n − 1 columns is established. As we can see, if the matrix A has n columns, then there were n I-complements {A 1 1 , A 1 2 , . . . , A 1 n }, and the set of these matrices is called I-minors of n − 1 columns of the matrix A. By employing the above procedure to each I-minor of n − 1 columns of the matrix A, I-minors of n − 2 columns of the matrix A can be archived. Recursively, we can obtain I-minors of n − h columns of the matrix A where h < n.
Based on Theorem 1 stated by Dines in [18], we can derive the sufficient condition in order to check the feasibility of system of inequalities in form of Ax > 0. Theorem 4. The system of inequalities in form of Ax > 0 is feasible if there exists at least one I-minor of n-h column which is I-definite where h < n.
Since the system (17) in form Ax > 0, so we can apply Theorem 4 to it and make our conclusion. (17) is feasible if there exists at least one I-minor of h columns which is I-positive.

Conclusion 2. The system of inequalities
After the whole process of checking the feasibility of given system (17), we have a sufficient condition to determine whether the pose x which generates the wrench matrix W Appl. Sci. 2021, 11, 4228 9 of 12 belongs to WCW or does not. The overview of the whole process to check the feasibility of certain pose is illustrated in Figure 6.

Case Studies and Simulation
In this section, the proposed wrench closure condition is applied to search for the wrench-closure workspace of two configurations CDPM defined in Figure 1. For convenience, we define three types of WCW: • Total orientation wrench closure workspace (TOWCW): Other than the definition in [22], in this article, TOWCW presents the whole WCW when position and orientation of moving platform are considered as variables. Notice that, because of the limitation in visualization, we can plot the TOWCW of 3-dof CDPM. As shown in Figure 7, WCW of the first configuration CDPM shown in Figure 1 is generated by point-wisely employing proposed WWC in the grid map, which has a dimension of 1 m × 1 m × 2π and resolution of 0.02 m × 0.02 m × 0.02 rad. Total orientation WCW of planar CDPM is a subset of 3-dimensional space (xP-yP-γP); it is illustrated in Figure 7a. Moreover, 2-dimensional sections of total orientation WCW of planar CDPM obtained at γP = 0 o , 0.052 rad and −0.052 rad are also plotted in Figure 7b-d. The orange, yellow and purple closed curves in Figure 7a represent the boundaries of slices that present in Figure 7c, Figure 7b and Figure 7d, respectively. However, this configuration is simple, since it has 3-dof and only one redundant actuator. Hence, it cannot verify the feature "applicable for a completely or redundantly restrained CDPM" of proposed WWC. The second configuration CDPM shown in Figure 1 is a 2-redundantly actuated spatial CDPM, and it can be used to reinforce the flexibility of proposed WWC.

Case Studies and Simulation
In this section, the proposed wrench closure condition is applied to search for the wrench-closure workspace of two configurations CDPM defined in Figure 1. For convenience, we define three types of WCW: • Total orientation wrench closure workspace (TOWCW): Other than the definition in [22], in this article, TOWCW presents the whole WCW when position and orientation of moving platform are considered as variables. Notice that, because of the limitation in visualization, we can plot the TOWCW of 3-dof CDPM. As shown in Figure 7, WCW of the first configuration CDPM shown in Figure 1 is generated by point-wisely employing proposed WWC in the grid map, which has a dimension of 1 m × 1 m × 2π and resolution of 0.02 m × 0.02 m × 0.02 rad. Total orientation WCW of planar CDPM is a subset of 3-dimensional space (x P -y P -γ P ); it is illustrated in Figure 7a. Moreover, 2-dimensional sections of total orientation WCW of planar CDPM obtained at γ P = 0 • , 0.052 rad and −0.052 rad are also plotted in Figure 7b-d. The orange, yellow and purple closed curves in Figure 7a represent the boundaries of slices that present in Figure 7b-d, respectively. However, this configuration is simple, since it has 3-dof and only one redundant actuator. Hence, it cannot verify the feature "applicable for a completely or redundantly restrained CDPM" of proposed WWC. The second configuration CDPM shown in Figure 1 is a 2-redundantly actuated spatial CDPM, and it can be used to reinforce the flexibility of proposed WWC.
As shown in Figure 8, WCW of the second configuration CDPM shown in Figure 1 is generated by applying proposed WWC. Constant orientation WCW with α P = 0, β P = 0 and γ P = 0.052 rad, and constant position WCW with x P = 0, y P = 0 and z P = 0 are shown in Figure 8a,b, respectively. As shown in Figures 7c and 8a, the cross-section perpendicular to the z-axis of COWCW of the spatial CDPM is trickly equivalent to the COWCW of the planar CDPM due to their similar dimensions. This comment is also verified by Pham et al. in [16]. In Figure 8b, the CPWCW is not continuous and in practice, this type of CDPM just can work on a small set of orientational spaces.  As shown in Figure 8, WCW of the second configuration CDPM shown in Figure 1 is generated by applying proposed WWC. Constant orientation WCW with αP = 0, βP = 0 and γP = 0.052 rad, and constant position WCW with xP = 0, yP = 0 and zP = 0 are shown in Figure 8a and Figure 8b, respectively. As shown in Figures 7c and 8a, the cross-section perpendicular to the z-axis of COWCW of the spatial CDPM is trickly equivalent to the COWCW of the planar CDPM due to their similar dimensions. This comment is also verified by Pham et al. in [16]. In Figure 8b, the CPWCW is not continuous and in practice, this type of CDPM just can work on a small set of orientational spaces.
In simulations, the time-consumption of proposed WCC used to determine whether a pose presents in WCW is also considered. Time-consumption is not the same in each performance due to two factors: size of wrench matrix and the interruption of algorithm. The first factor depends on configuration of CDPM, it is related to the number of used cables and DoF of moving platform. In the worst case, when a pose belongs to WCW and the algorithm for checking the feasibility of system of inequations must evaluate all I-minor of given wrench matrix. The second factor is the interruption of the algorithm. Our algorithm will orderly evaluate I-minor of the given wrench matrix until one of them is Idefinite. In the best case, the wrench matrix is not full-rank, the pose absolutely does not belong to WCW and the algorithm does not take place.

Conclusions
In this paper, the wrench closure condition that can be applied for a completely or redundantly restrained CDPM is our main contribution. Based on the properties of WCW, theorem 1 is given to determine the feasibility of a pose in WCW. By classifying the type In simulations, the time-consumption of proposed WCC used to determine whether a pose presents in WCW is also considered. Time-consumption is not the same in each performance due to two factors: size of wrench matrix and the interruption of algorithm. The first factor depends on configuration of CDPM, it is related to the number of used cables and DoF of moving platform. In the worst case, when a pose belongs to WCW and the algorithm for checking the feasibility of system of inequations must evaluate all I-minor of given wrench matrix. The second factor is the interruption of the algorithm. Our algorithm will orderly evaluate I-minor of the given wrench matrix until one of them is I-definite. In the best case, the wrench matrix is not full-rank, the pose absolutely does not belong to WCW and the algorithm does not take place.

Conclusions
In this paper, the wrench closure condition that can be applied for a completely or redundantly restrained CDPM is our main contribution. Based on the properties of WCW, theorem 1 is given to determine the feasibility of a pose in WCW. By classifying the type of CDPM, n = m + 1 and n > m + 1, and using linear equivalence theory, the sufficient condition of WCC is simplified. Finally, applying the algorithm to check the feasibility of system of inequalities, given that the pose of the moving platform determined belongs to WCW or does not.
Our proposed WCC proved to be synthetic and to be an algorithm to rapidly determine whether a pose belongs to the WCW or not. The overview of our proposed algorithm is illustrated in Figure 6. This algorithm can be combined with any search algorithm to generate the WCW. Two different cases study were also given, the first one is one redundant actuator planar CDPM and the other is two redundant actuators spatial CDPM. By employing the proposed algorithm with point-wise search method, the WCW or spatial WCW of these CDPM are determined and evaluated.