Open access

# Wire Robots Part II: Dynamics, Control & Application

Written By

Tobias Bruckmann, Lars Mikelsons, Thorsten Brandt, Manfred Hiller and Dieter Schramm

Published: April 1st, 2008

DOI: 10.5772/5366

From the Edited Volume

## Parallel Manipulators

Edited by Jee-Hwan Ryu

Chapter metrics overview

View Full Metrics

## 1. Introduction

In (Bruckmann et al., 2008) the kinematics, analysis and design of wire robots were presented. This chapter focuses on control and applications of wire robots. Wire robots are a very recent area of research. Nevertheless, the are well studied and already in application (see section 5). Due to their possible lightweight structure, wire robots can operate at very high velocities. Hence, as can be seen by experiment, only positioning control using the inverse kinematics is not sufficient. In particular, slackness in the wires can be observed at highly dynamic motions. To overcome this problem, force control can be employed. In section 4 different control schemes are proposed. The required dynamical model is obtained in section 2, while for the calculation of feasible wire force distributions are proposed in section 3. Since wire robots are kinematically redundant the latter is not straightforward, but requires advanced approaches. The same holds for the control schemes, since a CRPM as well as a RRPM is a non-linear, coupled, redundant system (Ming and Higuchi, 1994).

## 2. Dynamics

According to figure 1 a wire robot can be considered as a multibody system with m unilateral constraints. In contrast to the generally complicated forward kinematics (Bruckmann et al., 2008) the dynamical equations of motion are comparably easy to formulate with respect to the base frame . The wrench w wire of the wires acting on the platform can be written as (see 2)

E1

Since the forces act along the wires

E4

holds. It follows

E3

E4
E5

with

m p : the mass of platform,

I 3×3 : inertia tensor defined with respect to the inertial system which is an expression of rotation angles,

Ω = [ φ ϑ ψ ]T : orientation of the platform in ,

f E : vector of external forces,

τ E : vector of external torques.

The equations eqn. 4 can be rewritten by

E6

with

M p : mass matrix of platform,

E : identity matrix,

g C n×1 : Cartesian space vector of Coriolis and centrifugal forces and torques,

g E n×1 : vector of the generalized applied forces and torques, not including the resultants of wire tensions.

Taking wire force limits f min and f max (see (Bruckmann et al., 2008)) into account it follows

E7
E8

## 3. Wire force calculation

In section 2 a description of the force equilibrium was presented. Here methods for the calculation of a feasible force distribution f, i.e. a force distribution f which satisfies eqn. 7 and the constraints in eqn. 8, are presented. Obviously eqn. 7 represents an underdetermined system of linear equations. Its solution space is r-dimensional. Hence isolating the force distribution f leads to

E9

where A +T denotes the Moore-Penrose Pseudo-Inverse of A T. Thus the task of finding a feasible wire force distribution has been transformed to the task of finding λ r such that f> 0 holds. Note that H is the nullspace or kernel of A T defined as

E10

where

E11

In other words, a linear combination of the columns of H describes force distributions creating an inner tension in the system without applying wrenches w wire onto the end effector. In case of an homogenious problem, i.e. w = 0, it describes the possible solutions of eqn. 7 for f. Now the problem of satisfying the constraints of eqn. 8 arises, i.e. the force limits also have to be considered. Thus plugging eqn. 9 into eqn. 8 leads to

E12

Therefore the task of identifying a feasible force distribution is equivalent to the problem of identifying λ r such that eqn. 12 holds. In other words, the boundaries of the wire forces form a m-dimensional hypercube C m. All force distributions satisfying eqn. 9 obviously form a r-dimensional subspace S m spanned by the kernel of the structure matrix (see fig. 3). Hence, if the intersection F of the hypercube C and the subspace S is non-empty, feasible solutions f exist, i.e. F = C S (= ), where F is a r-dimensional manifold in the m. A more detailed introduction is given in (Oh and Agrawal, 2005) and (Mikelsons et al., 2008). Noteworthy, the r-dimensional solution space generally allows to compute force distributions with different characteristics: While for fast motion, smallest possible forces are demanded, for applications requiring a high stiffness, high forces are advantageous (Kawamura et al., 2000), (Fang, 2005).

### 3.1. Linear optimization

Looking at the geometric interpretation of finding feasible force distributions, the most intuitive way is to search for a convenient characterization of the manifold F. Since F is completely determined by its vertices, the computation of those seems to be a promising way. In this work, two approaches following this idea are shown: In section 3.3, a method using the kernel as a transformation is presented. This leads to (m/r) r -dimensional linear systems of equations. Alternatively, the approach presented in this section presumes no knowledge of the kernel but solves (m/r) n dimensional linear systems of equations. Hence, the method to be applied has to be chosen depending on m and n.

Examining eqn. 7, one needs to set r forces in the wire force distribution to get a quadratic system. Obviously the desired points are located on the faces of the cube C. It can be shown that a point belongs to the workspace if and only if a valid wire force distribution f that satisfies

For a set A, | A | denotes the cardinal number of A

E13

exists

Using the kernel as a transformation from the r into the m (see section 3.3), the feasible force distribution form a polyhedron bounded by the force limits. r force limits determine a vertex. This finishes the proof.

. Therefore, r wire forces can be set to their minimum or maximum value, respectively. It is unknown in advance which wire forces have to be preset to get a feasible distribution. Thus, in the worst case all combinations of r wires have to be tested, leaving m×m systems of linear equations to be solved for every combination. For sure every vertex represents a valid wire force distribution. Choosing the vertex, which minimizes the 1-norm could be an appropriate procedure. The resulting procedure can formally be expressed as a Linear Optimization Problem

In (Oh & Agrawal, 2005) a Linear Programming approach is presented to solve the problem in the r . Note that for control purposes, the Linear Optimization approach may deliver inadequate results since along a trajectory through the workspace, the result may be discontinuous.

### 3.2. Nonlinear optimization

Due to the formulation of the cost function, the Linear Programming method may deliver discontinuous solutions along a continuous trajectory. This leads to jumps in the time history of the wire forces, causing stability problems and additional mechanical wear. In (Verhoeven, 2004) it is proven that cost functions using a p-norm (1 < p <), lead to guaranteed continuous wire forces along a continuous trajectory. The resulting formulation of the optimization problem is as follows:

In (Verhoeven, 2004), also an effective algorithm is presented which solves the problem employing the knowledge of the solution structure, based on an iterative approximation of the optimal solution. However, this algorithm has the drawback to fail in specific configurations, i.e. solutions might be not found although they exist. To obtain the lowest possible force distribution (according to a p-norm), the unbounded polyhedron Plow is introduced, which is limited by the lower wire force limits:

E14

Furthermore, the wire force distribution flow is introduced, which has minimal p-norm:

E15

It should be mentioned that for 1 < p < flow is unique, which is essential for the continuity of flow. The algorithm works as follows

1. Compute an initial guess f ˜ low for flow.

2. If low is not contained in Plow, move low towards Plow until it is placed on the polyhedron.

3. Minimize the p-norm of low.

The initial guess is obtained by the orthogonal projection low of fmin onto the manifold of feasible force distributions F. Note that low is not always contained in Plow. The second step of the algorithm is performed by moving along the negative gradient of the distance between the polyhedron Plow and low. The distance is measured in the squared 2-norm. Finally, the minimization of low is done using a gradient based method again. Analogously, a vector fhigh representing the highest possible solution in the chosen p−norm can be obtained. Hence, choosing a wire force distribution on the line between flow and fhigh allows either fast motions due to low wire forces or high stiffness due to high wire forces. This approach is very effective in terms of computation time since the initial guess is often already a feasible solution, but suffers from the fact that a solution is not always found.

### 3.3. Barycentric force calculation

The shown approaches require the usage of an optimizer to deliver continuous results as shown in ((Verhoeven, 2004),(Nahon and Angeles, 1991), (Bruckmann et al., 2006), (Voglewede and Ebert-Uphoff, 2004) and (Bosscher and Ebert-Uphoff, 2004)). Standard optimizer implementations as LAPACK or the NAG ® library require iterative computations, which may not be used within a realtime control system due to their normally non-predictable worst-case runtime. In this section, a non-iterative algorithm is shown, which provides continuous force distributions furthermost from the force limits. The algorithm provides a force distribution, which lies in the center of gravity (CoG or barycenter) of the intersection manifold F.

The structure matrix AT has the dimension n × m. Hence, within the workspace, the kernel can be computed as H = (h1... hr ) m×r. Here, the kernel is used to define a map from the r to S m, i.e. for all λ , eqn. 12 must hold, where Λ is the (convex) polyhedron-shaped preimage of the manifold F under the mapping : r m, λ −A+Tw+H λ . In other words, since maps the r onto the solution subspace S, it maps the polyhedron r onto the solution manifold F. Since there is no explicit expression for , a convenient representation is sought. As mentioned above is a polyhedron. Thus, its vertices determine completely. Componentwise evaluation of both sides of eqn. 12 gives 2m hyperplanes in r. The vertices of are intersection points of r hyperplanes. Hence, all those intersection points are calculated and examined with respect to their compatibility with all inequalities. Obviously a vertex of the polyhedron has to satisfy all inequalities of eqn. 12. In order to compute the center of gravity of the obtained polyhedron, " is triangulated, i.e. splitted into r-simplexes. In the case of r = 2 this just means dividing into triangles. Advanced techniques as shown in (Cignoni et al., 1998) are required in the case of higher dimensions. Triangulation delivers a list of ns simplexes Pk with each having r + 1 vertices vkj with k = 1... ns and j = 1... r + 1. The volumes Vk of the simplexes can be determined by integration (Hammer et al., 1956). Furthermore their CoG λ sk are computed by the equation

E16

which is used to calculate the CoG λ s of the polyhedron via

E17

Finally, the solution is transformed back using the mapping

E18

where fs is the center of gravity of the manifold F.

#### 3.3.1. Proof-of-Concept

In this section we prove that the CoG of the manifold F can be computed by calculating the CoG of the convex polyhedron. Without loss of generality w = 0 is assumed. The CoG of the manifold F can be computed componentwise as

E19

The theorem for integration on manifolds states

E20

where H* : F, λ H is a map from to F and (DH)* is the Jacobian of H* which is equal to H itself since it is linear. Furthermore, d e t ( H T H )   is independent from and can hence be canceled in the next step. Additionally splitting into the simplexes gives:

E21

Since H is independent from , it can be moved out of the integral:

E22

Using eqn. 19 and eqn. 17 this can be rewritten as

E23

Therefore fs = H s holds where s denotes the CoG of in r.

#### 3.3.2. Continuity of solution

In this section the continuity of of the solution of the developed algorithm in the pnorm p ( p 1 , ) is proven, i.e. the function : m n n, which maps a matrix A m×n (considered as a vector in m n) onto the center of gravity as described before, is continuous on the set of points of the workspace.

Proof

Again without loss of generality w = 0 is assumed. First is splitted into two mappings Ker : m n nr and GravC : n r n. The latter maps a vector p from n r onto the center of gravity of the manifold F spanned by the r n-dimensional downwards listed vectors in p. Ker : m n nr maps a matrix A on its kernel H represented as a vector p in n r. In calculations the kernel is still denoted with H for simplicity. Continuity of Ker and GravC implies continuity of , since = GravC Ker.

First the continuity of GravC will be proven. Therefore 0 is assumed (i.e. the intersection of hypercube C and subspace S is non-empty and thus also the CoG exists), since continuity inside of C is to be proven. The CoG s is considered:

E24

Let λ ˜ s be the CoG of , where Λ ˜ is the preimage of F ˜ , which is obtained from H ˜ = H +E. The matrices H = [h1... hr]T n×r and E = [e1... er]T n×r are considered as vectors in n r. Then the p-norm of H is . It follows

E25
E26

Since the vertices of the polyhedron λ ˜ are obtained from the inequality

E27
E28

and the vertices of the polyhedron are obtained from (12), it is obvious that

E29
E30

Hence

E31

holds, because and are bounded. This yields together with eqn. (18)

E32
E33

This implies the continuity of GravC.

The continuity of Ker follows from the fact that the solution of a full ranked linear system of equations depends continuously on the coefficient matrix.

## 4. Control

Wire robots allow for very high velocities and accelerations when handling lightweight goods. In this case, wire robots benefit from their lightweight structure and low moved masses. Contrariwise, wire-based mechanisms like cranes, winches or lifting blocks are used widely to move extremely heavy loads. Thus, the wide range of application demands for a robust and responsive control. To move the platform along a trajectory precisely, position control is mandatory. On the other hand, the usage of wires claims for a careful observation and control of the applied tensions to guarantee a safe and accurate operation. Pure force control suffers from the drawbacks of model based control, e.g. model mismatch and parameter uncertainties. Thus force control is not sufficient and a combined force and position control is advised. Beside this, the relatively high elasticity of the wires may demand for a compensation by control. (Fang, 2005) shows more details of the shown concepts.

### 4.1. Elastic wire compensation

Compared to a conventional parallel kinematic machine (e.g. Stewart platform), a wire robot has generally a higher elasticity in the kinematic chains connecting the base and the platform. This is both due to the stiffness of the wire material as well as due to the wire construction (e.g. laid/twisted, braided or plaited)(Feyrer, 2000). Approximating the dynamical characteristics of the wires by a linear spring-damper model and considering the unilateral constraint, the wire model can be described as

E34

with 1 < i < m, ci and di denoting the stiffness and damping coefficients, respectively and Δ li denoting the length change due to elasticity. Assuming the untensed wire length is li,0, Δ li can be computed as Δ li = li − li, 0. The stiffness coefficient ci depends on the actual wire length. Using the wire cross section A and Young’s modulus E, ci can be calculated as

E35

with

E36

Note that this is only a linear approach. Taking into account long and heavy wires, a specific wire composition and applied tensions close to the admittible work load, advanced non-linear models have to be utilized. Especially the damping coefficient di may be hard to estimate (Wehking et al., 1999) and thus, experiments have to be carried out (Vogel & Götzelmann, 2002).

### 4.2. Motion control in joint space

The idea of motion control in joint space is to use a feedback position control and a feedforward force controller. The feedforward control employs an inverse dynamics model to calculate the winch torques necessary for the accelerations belonging to the desired trajectory. Since the used dynamic model usually will not cover all mechanical influences (e.g. friction), the remaining position errors can be compensated by the position control which employs the inverse kinematics. Noteworthy, the inverse dynamics is calculated for the desired platform position. Optionally, one may think of tracking control to guide the platform along the desired trajectory for the price of additional calculations. Referring to eqn. 6, the inverse system dynamics (i.e. the wire force distribution) can be computed by methods shown in section 3 (where the loads w include the inertia and gravity loads). Assuming the winch drives are adressable by desired torques (which is normally the case for DC/EC motors, preferably with digital current control), the motor dynamics can be modeled as

E37

where MM m×m is the inertial matrix of the drive units, η is the radius of the drums and D m×m depends on the structure of the motors. Combining the feedforward force control and the feedback position control leads to the following controller output:

E38

denoting the feedback gain matrices Kp m×m and Kd m×m and the actual and desired motor angles Θ and Θ d, respectively. Due to the decoupled position controllers, these may be designed as decentralized, simple and high control rate devices. To compensate for elastic tendons, the following correction may be applied:

E39

where Θ ^ d,i corresponds to the uncompensated drum angle (1 i m).

### 4.3. Motion control in operational space

Observing the sections above, independent linear PD controllers are applied. Practical experiences show that this is possible even though the system dynamics are described by a nonlinear, coupled system of equations due to the parallel topology of the robot, represented by the pose dependent structure matrix. Nevertheless, it is difficult to determine stable or even optimal controller parameters since the usual tools of the linear control theory may only be applied for locally linearized configurations of the robot. For predefined trajectories, this may be possible (e.g. by defining a cost function accumulating the control errors in simulation and applying a nonlinear optimizer to obain values for Kp and Kd), but is is desirable to have a globally linear system to avoid this only locally valid approach. From literature (Schwarz, 1991) (Woernle, 1995), exact linearization approaches are known which eliminate the nonlinear system characteristics by feedback. Using this as an inner loop, an outer linear controller may now be applied to the resulting linear system. Eqns. 37 and 6 deliver

E40

Since the final control law is formulated in the operational space, this equation is transformed into cartesian coordinates using the inverse kinematics relations

E41
E42

In cartesian coordinates the dynamical equations are then given by

E43

Instead of using the motor torques u as the system input, the resulting forces and torques acting onto the platform F ν are chosen to represent the actuator torques. Now a global linearization is desired. Setting F ν = Meq ν +N delivers

E44

and is therefore a proper choice. This linear system is now controlled by a PD controller for the position. Thus, the new system input is extended by

E45

Substituting eqn. 45 into eqn. 43, F ν can be found as

E46

which describes the required wrench onto the platform w which allows to calculate the desired wire forces by the methods shown in section 3. Optionally, the desired forces can be controlled by an outer feedback loop to enhance the control precision.

## 5. Applications

As already mentioned before, wire-based manipulation and construction is used since millenia, mostly taking advantage from the principle of the lifting block. In ancient civilisations like the Egypt of the Pharaos, probably wires and winches were applied to build the pyramids - wether using ramps or lifting mechanisms (see fig. 7(a)). Crane technology was only possible due to the usage of wires and especially the old Romans deleloped this technology to a remarkable state - they already lifted loads around 7 tons with cranes driven by 4 workers. With industrialisation, the transport and manipulation of heavy goods became very important, and hence, cranes using steel cables completed the transport chain for cargo handling. In the last few years, the automatisation of crane technology was subject to extensive research, e.g. in the project RoboCrane ® by the National Institute of Standards and Technology (NIST) (Bostelman et al., 2000). At the University of Rostock, the prototype CABLEV (Cable Levitation) (Maier, 2004),(Heyden, 2006) was build up, see fig. 8. It uses a gantry crane and three wires to guide the load along a trajectory. Thew load is stabilized by a tracking control for IRPM systems which eliminates

oscillations. In Japan, the Tadokoro Laboratory of the Tohoku University in Japan proposes the application of wires for rescue robots (Takemura et al., 2005) (Maeda et al., 1999). A problem solved very smart by usage of wires is the positioning of a large telescope. Several projects, e.g. the world’s largest telescope at Arecibo (fig. 7(b)), deal with the usage of wires to place the receiver module. The Arecibo project (900t receiver, approximately 300m satellite dish diameter) uses three wires guided by three mast heads while other project use an inverse configuration, lifting the receiver by balloons (see (Su et al., 2001), (Taghirad and Nahon, 2007a), (Taghirad and Nahon, 2007b)). Another popular application of wire robots is the usage as a manipulator for aerodynamical models in wind tunnels as proposed in (Lafourcade et al., 2002), (Zheng, 2006) and (Yaqing et al., 2007). Here, the experiments take advantage from the very thin wires since undisturbed air flow is mandatory. On the other hand, the wire robot can perform high dynamical motion as for example the FALCON (Fast Load Conveyance) robot (Kawamura et al., 1995). In the past few years at the Chair of Mechatronics the testbed for wire robots SEGESTA Seilgetriebene Stewart-Plattformen in Theorie und Anwendung) (Hiller et al., 2005b) has been developed. It is currently operated with seven (see fig. 9) wires in an CRPM configuration or eight wires for a RRPM setup. One focus of research is the development of fast and reliable methods for workspace calculation (Verhoeven and Hiller, 2000) and robot design. Another focus is the development of robust and realtime-capable control concepts(Mikelsons et al., 2008). Since the teststand is available, the theoretical results can be tested and verified (Hiller et al., 2005a). The system performs accelerations up to 10g and velocities arount 10m/s.

Another very recent application area has been created by Visual Act AB®. As pictured in fig. 10 a snowboard simulator was built up. The snowboarder is connected to four wires leading to three translational d.o.f.. Hence, the snowboarder can be guided along a trajectory in a setting consisting of ramps to grind on while moving freely in the air. (Visualact AB, 2006). A completely different field is the application of wire robots for rehabilitation which was demonstrated by the system String Man by the Fraunhofer-Institut für Produktionsanlagen und Konstruktionstechnik (IPK) in Berlin, Germany (Surdilovic et al., 2007) and the Lokomat protoype by the Eidgenössische Technische Hochschule in Zürich, Switzerland (Duschau-Wicke et al., 2007),(Wellner et al., 2007), (von Zitzewitz et al., 2006). Another prototype for rehabilitation is described in (Frey et al., 2006). The application of wire robots as a tracking device was proposed in (Ottaviano & Ceccarelli, 2006), (Thomas et al., 2003) and (Ottaviano et al., 2005). Here, the wire robot is not actively supporting a load but attached to an object which is tracked by the robot.

## Acknowledgments

This work is supported by the German Research Council (Deutsche Forschungsgemeinschaft) under HI370/24-1, HI370/19-3 and SCHR1176/1-2. The authors would like to thank Martin Langhammer for contributing the figure design.

## References

1. 1. Bosscher P. Ebert-Uphoff I. 2004 Wrench-based analysis of cable-driven robots. Proceedings of the 2004 IEEE International Conference on Robotics & Automation, 4950 4955 .
2. 2. Bostelman R. Jacoff A. Proctor F. 2000 Cable-based reconfigurable machines for large scale manufacturing. In Japan/USA Flexible Automation Conference Proceedings , University of Michigan, Ann Arbor, MI.
3. 3. Bruckmann T. Mikelsons L. Brandt T. Hiller M. Schramm D. 2008 Wire robots part I- kinematics, analysis & design. In Lazinica, A., editor, Parallel Manipulators, ARS Robotic Books. I-Tech Education and Publishing, Vienna, Austria. 978-3-90261-320-2
4. 4. Bruckmann T. Pott A. Hiller M. 2006 Calculating force distributions for redundantly actuated tendon-based Stewart platforms. In Lenarcic, J. and Roth, B., editors, Advances in Robot Kinematics- Mechanisms and Motion, 403 413 , Ljubljana, Slowenien. Advances in Robotics and Kinematics 2006, Springer Verlag, Dordrecht, The Netherlands.
5. 5. Cignoni P. Montani C. Scopigno R. 1998 Dewall: A fast divide and conquer delaunay triangulation algorithm in ed. Computer-Aided Design, 30 5 333 341 .
6. 6. Duschau-Wicke A. v. Zitzewitz J. Banz R. Riener R. . June 1215 2007). Iterative learning synchronization of robotic rehabilitation tasks. Proceedings of the 2007 IEEE 10th International Conference on Rehabilitation Robotics, 335 340 . Noordwijk, The Netherlands.
7. 7. Fang S. 2005 Design, Modeling and Motion Control of Tendon-based Parallel Manipulators. Ph. D. dissertation, Gerhard-Mercator-University, Duisburg, Germany. Fortschritt-Berichte VDI, Reihe 8, Nr. 1076, Düsseldorf.
8. 8. Feyrer K. 2000 Drahtseile. Springer Verlag Berlin.
9. 9. Frey M. Colombo G. Vaglio M. Bucher R. Jörg M. Riener R. 2006 A novel mechatronic body weight support system. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 14 3 311 321 .
10. 10. Hammer P. C. Marlowe O. P. Stroud A. H. 1956 Numerical integration over simplexes and cones. Math. Tables Aids Comp., 10 55 130 137 .
11. 11. Heyden T. 2006 Bahnregelung eines seilgeführten Handhabungssystems mit kinematisch unbestimmter Lastführung. PhD thesis, Universität Rostock. 3-18510-008-5Berichte VDI, Reihe 8, Nr. 1100, Düsseldorf.
12. 12. Hiller M. Fang S. Hass C. Bruckmann T. 2005a Analysis, realization and application of the tendon-based parallel robot segesta. In Last, P., Budde, C., and Wahl, F., editors, Robotic Systems for Handling and Assembly, 2 of International Colloquium of the Collaborative Research Center SFB 562, 185 202 , Braunschweig, Germany. Aachen, Shaker Verlag.
13. 13. Hiller M. Fang S. Mielczarek S. Verhoeven R. Franitza D. 2005b Design, analysis and realization of tendon-based parallel manipulators. Mechanism and Machine Theory, 40.
14. 14. Kawamura S. Choe W. Tanaka S. Pandian S. R. 1995 Development of an ultrahigh speed robot falcon using wire drive system. IEEE International Conference on Robotics and Automation, 215 220 .
15. 15. Kawamura S. Kino H. Won C. 2000 High-speed manipulation by using parallel wire-driven robots. Robotica, 18 1 13 21 .
16. 16. Lafourcade P. Llibre M. Reboulet C. . October 34 2002). Design of a parallel wire-driven manipulator for wind tunnels. In Gosselin, C. M. and Ebert-Uphoff, I., editors, Workshop on Fundamental Issues and Future Research Directions for Parallel Mechanisms and Manipulators.
17. 17. Maeda K. Tadokoro S. Takamori T. Hattori M. Hiller M. Verhoeven R. 1999 On design of a redundant wire-driven parallel robot warp manipulator. Proceedings of IEEE International Conference on Robotics and Automation, 895 900 .
18. 18. Maier T. 2004 Bahnsteuerung eines seilgeführten Handhabungssystems- Modellbildung, Simulation und Experiment. PhD thesis, Universität Rostock, Brandenburg. Fortschritt-Berichte VDI, Reihe 8, Nr. 1047, Düsseldorf.
19. 19. Mikelsons L. Bruckmann T. Hiller M. Schramm D. 2008 A real-time capable force calculation algorithm for redundant tendon-based parallel manipulators. appears in Proceedings on IEEE International Conference on Robotics and Automation.
20. 20. Ming A. Higuchi T. 1994 Study on multiple degree of freedom positioning mechanisms using wires, part 1- concept, design and control. International Journal of the Japan Society for Precision Engineering, 28 131 138 .
21. 21. Nahon M. Angeles J. 1991 Real-time force optimization in parallel kinematics chains under inequality constraints. In IEEE International Conference on Robotics and Automation, 2198 2203 , Sacramento.
22. 22. Oh S. R. Agrawal S. K. 2005 Cable suspended planar robots with redundant cables: Controllers with positive tensions. In IEEE Transactions on Robotics.
23. 23. Ottaviano E. Ceccarelli M. 2006 Numerical and experimental characterization of singularities of a six-wire parallel architecture. Robotica, 25 3 315 324 .
24. 24. Ottaviano E. Ceccarelli M. Paone A. Carbone G. (April 18-22 2005). A low-cost easy operation 4-cable driven parallel manipulator. In Proceedings of the 2005 IEEE International Conference on Robotics and Automation, pages 4008-4013, Barcelona, Spain.
25. 25. Schwarz H. 1991 Nichtlineare Regelungssysteme. Oldenbourg, München. ISBN- 13 978 3486218336 .
26. 26. Su Y. X. Duan B. Y. Nan R. D. Peng B. 2001 Development of a large parallel-cable manipulator for the feed-supporting system of a next-generation large radio telescope. In Journal of Robotic Systems, 18 633 643 .
27. 27. Surdilovic D. Zhang J. Bernhardt R. 1315 June 2007). String-man: Wirerobot technology for safe, flexible and human-friendly gait rehabilitation. In Proccedings of IEEE 10th International Conference on Rehabilitation Robotics, 2007, 446 453 , Noordwijk, Netherlands. 978-1-42441-320-1
28. 28. Taghirad H. Nahon M. 2007a Forward kinematics of a macro-micro parallel manipulator. In Proceedings of the 2007 IEEE/ASME International Conference on Advanced Intelligent Mechatronics (AIM2007), Zurich, Switzerland.
29. 29. Taghirad H. Nahon M. 2007b Jacobian analysis of a macro-micro parallel manipulator. In Proceedings of the 2007 IEEE/ASME International Conference on Advanced Intelligent Mechatronics (AIM2007), Zurich, Switzerland.
30. 30. Takemura F. Enomoto M. Tanaka T. Denou K. Kobayashi Y. Tadokoro S. 2005 Development of the balloon-cable driven robot for information collection from sky and proposal of the search strategy at a major disaster. In Proceedings on IEEE/ASME International Conference on Advanced Intelligent Mechatronics, 658 663 , Monterey.
31. 31. Thomas F. Ottaviano E. Ros L. Ceccarelli M. . September 1419 2003). Coordinate-free formulation of a 3-2-1 wire-based tracking device using cayleymenger determinants. In Proceedings of the 2003 IEEE International Conference on Robotics and Automation, 355 361 , Taipei, Taiwan.
32. 32. Verhoeven R. 2004 Analysis of the Workspace of Tendon-based Stewart Platforms. PhD thesis, University of Duisburg-Essen.
33. 33. Verhoeven R. Hiller M. 2000 Estimating the controllable workspace of tendonbased Stewart platforms. In Proceedings of the ARK’00: 7th. International Symposium on Advances in Robot Kinematics, 277 284 , Portoroz, Slovenia.
34. 34. Visualact A. B. 2006 Visual act 3d. http://www.visualact.net/.
35. 35. Vogel W. Götzelmann B. 2002 Kraft in Faserseilen bei ausgewählten stossartigen Beanspruchungen. EUROSEIL, 121(3):44/45.
36. 36. Voglewede P. Ebert-Uphoff I. 2004 On the connections between cable-driven robots, parallel manipulators and grasping. In IEEE International Conference on Robotics and Automation, 5 4521 4526 , New Orleans. IEEE.
37. 37. von Zitzewitz. J. Duschau Wicke. A. Wellner M. Lünenburger L. Riener R. 2006 Path control: A new approach in patient-cooperative gait training with the rehabilitation robot lokomat. Gemeinsame Jahrestagung der Deutschen, Österreichischenund Schweizerischen Gesellschaften für Biomedizinische Technik. Zürich, Schweiz.
38. 38. Wehking K. H. Vogel W. Schulz R. 1999 Dämpfungsverhalten von Drahtseilen. F+H Fördern und Heben, 49(1-2):60-61.
39. 39. Wellner M. Guidali M. Zitzewitz J. Riener R. . (June 12-15, 2007). Using a robotic gait orthosis as haptic display- a perception-based optimization approach. Proceedings of the 2007 IEEE 10th International Conference on Rehabilitation Robotics, 81 88 . Noordwijk, The Netherlands.
40. 40. Woernle C. 1995 Regelung von Mehrkörpersystemen durch externe Linearisierung. 517 in Fortschrittberichte VDI, Reihe 8. VDI Verlag, Düsseldorf.
41. 41. Yaqing Z. Qi L. Xiongwei L. . (June 18-2, 2007) Initial test of a wiredriven parallel suspension system for low speed wind tunnels. In Proceedings on 12thIFToMM World Congress, Besançon, France.
42. 42. Zheng Y.-Q. 2006 Feedback linearization control of a wire-driven parallel support system in wind tunnels. Sixth International Conference on Intelligent Systems Design and Applications, 3 9 13 .

## Notes

• For a set A, | A | denotes the cardinal number of A
• Using the kernel as a transformation from the ℝ r into the ℝ m (see section 3.3), the feasible force distribution form a polyhedron bounded by the force limits. r force limits determine a vertex. This finishes the proof.

Written By

Tobias Bruckmann, Lars Mikelsons, Thorsten Brandt, Manfred Hiller and Dieter Schramm

Published: April 1st, 2008