Equilibrium Conditions of Class 1 Tensegrity Structures

Transcription

1 Revue francaise de genie civil, July, 2003 Equilibrium Conditions of Class 1 Tensegrity Structures Darrell Williamson, Robert E. Skelton, JeongHeon Han University of Wollongong, Faculty of Informatics, Wollongong, NSW 2522 Australia, darrel@titr.uow.edu.au University California San Diego Dept. of Mechanical and Aerospace Engineering 9500 Gilman Dr, La Jolla, CA 92093, USA, bobskelton@ucsd.edu University California San Diego, Dept. of Mechanical and Aerospace Engineering, 9500 Gilman Dr., La Jolla CA 92093, USA, jehan@ucsd.edu ABSTRACT. Static models of tensegrity structures are reduced to linear algebra problems, after first characterizing the problem in a vector space where direction cosines are not needed. That is, we describe the components of all member vectors as opposed to the usual practice of characterizing the statics problem in terms of the magnitude of tension vectors. While our approach enlarges (by a factor of 3) the vector space required to describe the problem, the computational space is not increased. The advantage of enlarging the vector space makes the mathematical structure of the problem amenable to linear algebra treatment. Using the linear algebraic techniques, many variables are eliminated from the final existence equations. This paper characterizes the existence conditions for all class 1 tensegrity equilibria KEYWORDS: tensegrity, statics, linear algebra.

2 292 RFGC - 7/2003. Tenségrité 1. Introduction The Tensegrity structures represent a wonderful blend of geometry and mechanics. In addition, they have engineering appeal in problems requiring large changes is structural shape. We define class 1 tensegrity as a stable connection of bars and strings with only one bar connected to any given node without any constraints. Nodes connecting more than one bar or other constraints form class 2 tensegrity. Thus, class 1 structures are unconstrained, whereas, class 2 may have boundary constraints or bar to bar connections. Most existing smart structure methods are limited to small displacements. Since class 1 tensegrity structures have no bar-to-bar connections, the control of tendons allows very large shape changes. Therefore, an efficient set of analytical tools could be the enabler to a hoist of new engineering concepts for deployable and shape controllable structures. Connelly, Motro, Pelligrino and others have described many details of the static equilibria of tensegrity structures. However, the dynamics and control of these structures is relatively new. This paper characterizes the static equilibria of class 1 tensegrity structures. Furthermore, we use vectors to describe each element (bars and tendons), eliminating the need to use direction cosines and the subsequent transcendental functions that follow their use. It is well known in a variety of mathematical problems that enlarging the domain in which the problem is posed can often simplify the mathematical treatment. Many nonlinear problems admit solutions by linear techniques by enlarging the domain of the problem. The purpose of this paper is to show that by enlarging the vector space in which we characterize the tensegrity statics problem, the mathematical structure of the equations admit treatment by linear algebra methods, for the most part. Our results characterize the equilibria conditions of class 1 tensegrity structures in terms of a very small number of variables, since the necessary and sufficient conditions of the linear algebra treatment allows the elimination of several of the original variables. These results can be programmed into object-oriented software to design and simulate a large class of tensegrity structures. This paper provides the enabling technology for efficient algorithms to design tensegrity structures, which have been around for fifty years without efficient design procedures. Section 2 introduces the network representations of tensegrity structures as an oriented graph in real three dimensional space. Geometric connectivity, equilibrium, and a coordinate transformation will be introduced. Section 3 introduces the algebraic equilibrium conditions for a class 1 tensegrity structure. After we derive necessary and sufficient conditions for the existence of an unloaded tensegrity structure in equilibrium, we write the necessary and sufficient conditions for the externally loaded structure in equilibrium. A couple of examples will show how to construct a tensegrity structure that concludes the paper. For notation, we let I n define the n n identity matrix, and 0 define an n m matrix of zeros. (The dimensions of 0 will be clear from the context.) We also let ρ(a) define the rank of the matrix A.

3 Equilibrium Conditions Network Representation of Structures In this paper, we choose to represent a tensegrity structure as an oriented graph in real three dimensional space R 3 defined in terms of n p nodes and n s + n b directed branches which are all represented as vectors in R 3. A loop is any closed path in the graph. As we shall see, the advantage of this approach is that the both the magnitude and the direction cosines of the forces are contained in vectors which can be solved using linear algebra. Thus linear algebra plays a larger role in this approach compared to the usual approach in mechanics and finite element methods using direction cosines. In this oriented graph, the nodes consist of the ends of the bars as represented by the n p nodes (or vectors) {p k }. Hence if there are n b bars, then there are n p = 2n b nodes. We choose to identify two different types of directed branches ; the n s string branches (or vectors) {s n } and the n b bar branches (or vectors) {b m }. Geometric Connectivity Each directed branch can undergo a displacement in reaching its equilibrium state. String vectors can change both their length and orientation while bar vectors can only change their orientation. Node vectors can change both their length and orientation but subject to a Law of Geometric Connectivity which we state as follows : The vector sum of all branch vectors in any loop is zero. [1 These equations are in the form of a set of linear algebraic equations in the branch vectors. Force Equilibrium In our study of tensegrity structures, we are concerned with structures in which bars sustain compressive forces. We therefore choose to distinguish between the string (or tensile) forces {t n } and the bar (or compressive) forces {f m } which are defined in terms of the string and bar vectors respectively as follows. Definition 1 Given the tensile force t n in the string characterized by the string vector s n and the compressive force f n in the bar characterized by the bar vector b n, the tensile force coefficient γ m 0 and the compressive force coefficient λ n 0 are defined by t n = γ n s n ; f m = λ m b m [2 The force of the tensegrity structure is defined by the external force vector w R 3np, the tensegrity compression vector f R 3n b, and the tensegrity tension vector t R 3ns where w T = [w T 1,wT 2,...,wT n p ; f T = [f T 1,fT 2,...,fT n b ; t T = [t T 1,tT 2,...,tT n s [3

4 294 RFGC - 7/2003. Tenségrité Force Convention Suppose each node p k is subjected to compressive vector forces {f mk }, tensile vector forces {t nk } and external force w k. Then the Law for Static Equilibrium may be stated as follows : t nk f mk w k = 0 [4 n m where a positive sign is assigned to a (tensile, compressive, external) force leaving a node, and a negative sign is assigned to a (tensile, compressive, external) force entering a node. The negative sign in (4) is a consequence of the fact that we choose to define nonnegative compressive force coefficients λ n. Consider a class 1 tensegrity structure consisting of n p nodes, n b bars and n s strings. Suppose the positions of the positions of the nodes are described by the n p vectors {p 1,p 2,...,p np }, the positions of the bars are described by the n b vectors {b 1,b 2,...,b nb }, and the positions of the strings are described by the n s vectors {s 1,s 2,...,s ns }. Definition 2 The geometry of the tensegrity structure is defined by the tensegrity node vector p R 3np, the tensegrity bar vector b R 3n b, and the tensegrity string vector s R 3ns where p T = [p T 1,p T 2,...,p T n p ; b T = [b T 1,b T 2,...,b T n b ; s T = [s T 1,s T 2,...,s T n s [5 From network analysis and the law for static equilibrium, we have the following result. Lemma 1 There exists an 3n p 3n s matrix S and an 3n p 3n b matrix B such that or equivalently [ t A f = w ; A = [S,B [6 St = Bf + w. [7 In particular, if we consecutively number the n s + n b branches {s 1,s 2,...,s ns, b 1,b 2,...,b nb } as {1, 2,..., n s, n s + 1,..., n s + n b }, then the 3n p (3n s + 3n b ) matrix A = [A ij is defined by I 3 if tension (compression) j leaves (enters) node i A ij = I 3 if tension (compression) j enters (leaves) node i 0 if tension (compression) j is not incident with node i Also : (i) each column of A has exactly one block I 3 and one block I 3 with all other column blocks 0, and (ii) for any row i there exists a column j such that A ij = ±I 3. [8

5 Equilibrium Conditions 295 Figure 1. 2-Bar 4-String Class 1 Tensegrity In network analysis, the matrix A is known as the incidence matrix. This matrix is not the reduced incidence matrix since we have included the datum node which means that one block row of equations in (6) are dependent on the other rows. This fact does not cause any difficulties in subsequent developments. On the contrary, some symmetry is preserved in the algebraic equations. Network analysis also gives us the following result. Lemma 2 The string vector t and the bar vector b are related to the node vector p by the equation A T p = [ s b or S T p = s B T p = b [9 The proof is constructive. Firstly, from the network, it follows that components of the string vector s and the bar vector b can be written as a linear combination of components of the node vector p. Also if branch k is a bar which leaves node i and enters node j, then b k = p j p i, whereas if branch k is a string which leaves node i and enters node j, then s k = p j p i. Hence we have Cp = [s T,b T where the matrix C consists only of block matrices of the form {0, ±I 3 }. A closer examination then reveals that C = A T. Example 1 Consider the 2-bar 4-string class 1 tensegrity structure illustrated in Fig. 1 with tensile force vectors {t 1,t 2,t 3,t 4 } and compressive force vectors {f 1,f 2 }. The Geometric Connectivity conditions are : p 1 + p 3 = s 1 ; p 2 p 3 = s 2 p 2 + p 4 = s 3 ; p 1 p 4 = s 4 [10

6 296 RFGC - 7/2003. Tenségrité Also, in terms of the stated force convention, the conditions for Static Equilibrium at nodes 1-4 are : t 1 t 4 f 1 w 1 = 0 ; t 2 + t 3 + f 1 w 2 = 0 t 1 + t 2 f 2 w 3 = 0 ; t 3 + t 4 + f 2 w 4 = 0 [11 The static equilibrium conditions and the geometric conditions can be written in the form (6), (8), (9) where with A = I I 3 I I 3 I 3 0 I 3 0 I 3 I I I 3 I 3 0 I 3 S = I I 3 0 I 3 I 3 0 I 3 I I 3 I 3 ; B = I 3 0 I I 3 0 I 3 We now derive a canonical network representation in which the bar vectors are a subset of the transformed network nodes. Lemma 3 Given the tensegrity node vector p, define the nodal transformation for some nonsingular matrix P R 3np 3np where q = P 1 p [12 P = [P T 1,PT 2 ; P 1,P 2 R 3n b 3n p [13 Then : (i) In terms of the transformed tensegrity node q, the tensegrity geometry is given by where S T q q = s ; BT q q = b [14 S T q = [ST 1,ST 2 ; BT q = [BT 1,BT 2 [15

7 Equilibrium Conditions 297 with S 1 = P 1 S ; S 2 = P 2 S ; B 1 = P 1 B ; B 2 = P 2 B [16 (ii) The tensegrity force equilibrium is given by S 1 t = B 1 f + P 1 w ; S 2 t = B 2 f + P 2 w [17 Proof : Part (i) follows directly from the definition of P. Part (ii) follows from the expansion of the equilibrium force equation S q t = B q f + P T w. Lemma 4 Given the tensegrity node vector p, bar vector b and string vector s, the geometry of the class 1 tensegrity structure can be described by the algebraic equations where the 3n p 3n b matrix B is given by B T = B T p = b ; S T p = s [18 I 3 I I 3 I I 3 I 3 [19 for some 3n p 3n s matrix S which consist only of block matrices {I 3, I 3,0 3 } with In particular : (i) P 1,P 2 R 3n b 3n p in (13) are given by n p = 2n b [20 P 1 = P 2 = I I I I I 3 I I 3 I I 3 I I 3 I 3 [21

8 298 RFGC - 7/2003. Tenségrité (ii) The inverse transformation P 1 = [P T 1,P T 2 1 = [ B T J T [22 where B is given by (19) and J R 3np 3n b is given by J T = 0 3 I I I I 3 [23 (iii) The transformed tensegrity node q is given by where q b,q e R 3n b are given by q T = [q T b,qt e [24 with q b = b ; q T e = [qt 2,qT 4,...,qT n p [25 S T q = [ST 1,ST 2 ; BT q = [I n b,0 nb [26 where S 1 R 3n b 3n s ; S 2 R 3n b 3n s are given by (16). (iv) The transformed force equilibrium conditions (17) are given by S 1 t = f + P 1 w ; S 2 t = P 2 w [27 Proof : Since each bar has two end points, n p = 2n b. Without any loss of generality, we can assume in a class 1 tensegrity that b m = p 2m 1 + p 2m, m = 1, 2,..., n b which then gives B T p = b with B as in (19). The block structure of the S matrix comes from the vector equations which relate the {s n,p r }. It follows from (19) and (21) that P T B = [ P1 P 2 B = which then gives (27) and (16). [ I3nb 0 3nb ; [ S1 = S 2 [ P1 P 2 S = P T S

9 Equilibrium Conditions 299 Example 2 The equations (10), (11) of the 2-bar 4-string tensegrity introduced in Example 1 can be written in the form (6) and (9) where B T = [ I3 I ; S 0 0 I 3 I T = 3 Then from Lemma 3, we have (27) where I 3 0 I I 3 I I 3 0 I 3 I I 3 S 1 = P 1 = [ [ I3 0 0 I 3 I3 I ; S I 3 I = 3 I 3 I 3 I 3 I 3 I 3 I 3 [ [ I I3 I ; P 0 0 I = I 3 I 3 [28 3. Algebraic Conditions for Equilibrium Definition 3 A tensegrity structure with tensile force coefficients {γ n 0}, compressive force coefficients {λ n 0}, node vector p, string vector s and bar vector b is said to be in equilibrium if the element relationships (2), the force equations (7) and the geometric equations (9) are all satisfied. From lemma 1 and lemma 2, we have the that t T S T p f T B T p = w T p, or equivalently, t T s f T b = w T p. Hence from (2), we deduce the mechanical equivalent of Tellegen s Theorem for electrical networks. Corollary 1 For a tensegrity structure in equilibrium γ n s n 2 n m λ m b m 2 = k w T k p k [29 Connelly defined a tensegrity framework G(p) as a graph on p where each edge is either a cable or a bar. A stress state {γ n, λ m } for G(p) is said to be a self-stress if in the absence of an external load (i.e. w = 0), force equilibrium is satisfied at each node. The force t n in a linear spring of length s n having spring constant k n and rest length l n0 is given by t n = k n ( s n l n0 )

10 300 RFGC - 7/2003. Tenségrité Hence in terms of the tensile force coefficient γ n in (2), we have The potential energy PE n in the spring is also given by γ n = k n (1 l n0 s n ) [30 so that when the cable has zero rest length : PE n = 1 2 k n( s n l n0 ) 2 [31 γ n = k n is constant ; PE n = 1 2 γ n s n 2 [32 In analogy with total potential energy, Connelly defined the energy form associated with the stress state {γ n, λ m } as : E(p) = 1 2 { n γ n s n 2 m λ m b m 2 } [33 where all members are assumed to behave as elastic strings in which the cables have zero rest length and the bars have infinite rest length. The idea is that when the end points of a bar are displaced, the energy builds up as the square of the extension. The function E(p) when w = 0 is said to have an absolute minimum corresponding to the rest length of the element [Connelly In this context, the more general statement in Corollary 1 then says then in the presence of external forces {w k }, the energy associated with the stressed state must be balanced by the sum of the external work done w T k p k at each node. Requirements for Equilibrium : Given an external force vector w, the problem of determining the geometric and force configuration of a tensegrity structure consisting of n s strings and n b bars in equilibrium is therefore equivalent to finding a solution b,q e R 3n b of the equations : s = S T 1 b + S T 2 q e [34 t = Γs ; Γ = diag{γ 1 I 3, γ 2 I 3,..., γ ns I 3 } [35 S 2 t = P 2 w [36 f = S 1 t P 1 w [37 f = Λb ; Λ = diag{λ 1 I 3, λ 2 I 3,..., λ nb I 3 } [38 for given matrices

11 Equilibrium Conditions 301 where S 1,S 2 R 3n b 3n s ; P 1,P 2 R 3n b 3n p [39 n p = 2n b ; n s > n b [40 The diagonal matrices {Γ, Λ} shall be referred to as the tensile force matrix and compressive force matrix respectively. Materials Properties : Conditions for equilibrium will depend on both the tensile force coefficients of the strings and the compressive force coefficients of the bars. A tensile force coefficient γ n of a (linear) string may be designed by selection of the spring constant k n and the ratio l m0 / s n of the rest to stretched length of the string subject only to the requirement that the yield force (i.e. the force that causes the string to reach its elastic limit) is not exceeded. In the (mathematically ideal) case of zero rest length, k n = γ n. On the other hand, the compressive force coefficient λ n of a bar cannot be designed by an adjustment of the materials properties. However the bar force must not exceed the buckling force where an ideal bar of radius r 0 and length L 0 buckles at a compressive force F 0 given by F 0 = E 0πr 4 0 L 2 0 ; E 0 = Young s Modulus of bar Using an earlier derivation, the next result provides an expression for both the tensile and compressive force coefficients in an equilibrium structure. Lemma 5 Suppose the transformed tensegrity structure is in equilibrium. Then : (i) Assuming a linear force relationship for the strings, the tensile force coefficient {γ n } in the string s n is given by (30) where k n > 0 is the spring constant, l n0 > 0 is the rest length of the nth string, and s n is the stretched length of the nth string (ii) The compressive force coefficient {λ m } in the bar b m is given by λ m = f m L m ; L m = bm [41 where L m is the (constant) length of the bar, and f m = S m1 Γs P m1 w ; s = S T q q Γ = diag{γ 1 I 3, γ 2 I 3,.., γ ns I 3 } [42 where S m1,p m1 R 3 3np for 1 m n b are the mth blocks of {S 1,P 1 } respectively in (21), (16) where

12 302 RFGC - 7/2003. Tenségrité S 1 = S 11 S 21 S nb,1 ; P 1 = P 11 P 21 P nb,1 [43 Proof : From (2), (27) and (43), we have λ m b m = f m. Then since λ m 0, we get (41) Prestressed Equilibrium Structure We now proceed to derive necessary and sufficient conditions for the existence of a structure in equilibrium that is prestressed in the absence of any external load (i.e. w = 0). The derivation relies heavily on the singular value decomposition svd(a) = {U A,Σ A,V A } of a matrix A. Lemma 6 (i) Suppose an n m matrix A has rank r A, then there exists an n n unitary matrix U A, an m m unitary matrix V A and a positive definite r A r A diagonal matrix Σ 1A such that A = U A Σ A V T A ; Σ A = [ Σ1A [44 (ii) If {U A,V A } are partitioned such that with U A = [U 1A,U 2A ; V = [V 1A,V 2A [45 U 1A R n ra ; U 2A R n (n ra) ; V 1A R m ra ; V 2A R m (m ra) [46 then U T 1A.U 1A = I ra ; U T 1A.U 2A = 0 ; U T 2A.U 2A = I n ra V T 1A.V 1A = I ra ; V T 1A.V 2A = 0 ; V T 2A.V 2A = I m ra U T 2A A = 0 ; AV 2A = 0 [47 (iii) The algebraic equation

13 Equilibrium Conditions 303 Ax = b has a solution if and only if U T 2Ab = 0. When this condition is satisfied, then all solutions x are of the form where z x R n ra is arbitrary. x = V 1A Σ 1 1A UT 1Ab + V 2A z x ; U T 2Ab = 0 We now establish necessary and sufficient conditions for a solution of equations (34)- (40) in the absence of external forces (i.e. w = 0) by examining each of these equations in turn beginning with the solution of (34), (35). The next result follows from lemma 6. Lemma 7 Suppose r = ρ(s 2 ) min{3n p, 3n s } [48 and let S 2 have the singular value decomposition {U,Σ,V} given by where S 2 = UΣV T R 3n b 3n s [49 U = [U 1,U 2 ; Σ = [ Σ ; V = [V 1,V 2 U 1,V 1 R 3n b r ; U 2,V 2 R 3n b (3n b r) [50 Then a necessary and sufficient condition for (34) to have a solution q e R 3n b is given by V2 T (s ST 1 b) = 0 [51 Furthermore, when (51) is satisfied, all solutions q e are of the form where z e R 3n b r is arbitrary. q e = U 1 Σ 1 11 VT 1 (s S T 1 b) + U 2 z e [52 We now consider the solution of (36) when w = 0.

14 304 RFGC - 7/2003. Tenségrité Lemma 8 When w = 0, all solutions t of (36) which guarantee (51) are of the form t = V 2 M 1 V T 2 S T 1 b ; M = V T 2 Γ 1 V 2 [53 Proof : From (49), (50) and lemma 6, we have t = V 2 z t where z t is the free solution of (36). Then from (35) V T 2 (s ST 1 b) = VT 2 (Γ 1 V 2 z t S T 1 b) Since V 2 has full column rank, the matrix M = V2 T Γ 1 V 2 is invertible. Hence (51) is satisfied when z t = M 1 V2 TST 1 b, and this gives (53). We now consider the solution of (37), (38) when w = 0. Lemma 9 When w = 0, a necessary and sufficient condition for (37), (38) to have a solution b R 3n b is given by where (X Λ)b = 0 [54 In particular, define X = (S 1 V 2 )M 1 (S 1 V 2 ) T ; M = V T 2 Γ 1 V 2 [55 r b = ρ(λ X) [56 Then : (i) When r b = 3n b, then b = 0 is the only solution of (54), (ii) When r b = 0, any b R 3n b is a solution of (54), and (iii) When 0 < r b < 3n b, all solutions b satisfy the equation [I V X1 Σ 1 X11 UT X1 Λb = U X2z b [57 where z b R 3n b r b is free, and where {U X,Σ X,V X } is the singular value decomposition of the matrix X R 3n b 3n b ; that is U X = [U X1,U X2 ; Σ b = [ ΣX ; V X = [V X1,V X2 ; r X = rank(x) with U X1 R 3n b r X, V X2 R 3n b (3n b r X).

15 Equilibrium Conditions 305 Proof : From (53), (38) and (36) S 1 t f = S 1 V 2 M 1 V T 2 S T 1 b Λb = (X Λ)b The result then follows from the singular value decomposition of X after writing (54) in the form Xb = f, f = Λb. Now Γ positive definite implies M in (55) is positive definite. Then M positive definite implies X in (55) is positive semidefinite with at least one non-zero eigenvalue. Since (54) has a nontrivial solution b if and only if det(x Γ) = 0, the following lemma is important. Lemma 10 Given a symmetric positive semidefinite matrix X. Then there exists a positive definite diagonal matrix Λ such that if and only if X 0. det(x Λ) = 0 Proof : Since X = X T 0, X cannot have negative eigenvalues. Suppose λ 0 is a positive eigenvalue of X. Then Λ = λ 0 I implies det(λ X) = det(λ 0 I X) = 0. Remarks : For 0 < r b < 3n b, (Λ X)b = 0 always has at least one solution b 0 of the form (57) where the freedom in the choice of b is available in the choice of a free vector z b. However for n b 3, an additional necessary condition for the existence of a nontrivial tensegrity structure is that the n b bar vectors {b m } span R 3. This condition places additional requirements on the component matrices of the singular value decomposition of the matrix X. In the 2-bar 4-string planar tensegrity and the (3; 9; 3) 3-bar 9-string single stage shell tensegrity structure, no nontrivial tensegrity structures have been found in the case when r b 0 ; only in the case when r b = 0 do equilibrium solutions exist. We therefore specialize our results to this case. Theorem 1 Consider a class 1 tensegrity structure as defined by the geometry and force equations in the absence of external load as described by the geometric conditions B T p = b ; S T p = s ; p R 3np, b R 3n b, s R 3ns with n p = 2n b, and the equilibrium force equations where for {γ m > 0, λ n > 0}, St = Bf ; t = Γs, f = Λb Γ = diag{γ 1 I 3, γ 2 I 3,..., γ ns I 3 } ; Λ = diag{λ 1 I 3, λ 2 I 3,..., λ nb I 3 }

16 306 RFGC - 7/2003. Tenségrité Then given any tensile force coefficients {γ m > 0, 1 m n s }, there exits compressive force coefficients {λ n > 0, 1 m n b } which define an equilibrium structure if {Λ,Γ} satisfy the condition : Λ = X ; X = (S 1 V 2 )M 1 (S 1 V 2 ) T, M = V 2 Γ 1 V T 2 [58 Moreover, given condition (58), any vector p of the form where p = PQ[b T,z T e [59 Q = [ I3nb 0 ; L = U L U 1 Σ 1 11 VT 1 [Γ 1 V 2 M 1 V2 T I 3n s S T 1 [60 2 for arbitrary vectors {b R 3n b, z e R 3n b r } defines a corresponding tensegrity node. Equivalently, the tensegrity node is constrained according to the equation U T Q2 P 1 p = 0 [61 where U Q Σ Q VQ T is the singular value decomposition of Q in (60) with U Q = [U Q1,U Q2 ; U Q2 R 3np (3np rq), r Q = ρ(q) [62 The corresponding tensegrity tension vector t, string vector s, and compression vector f are given in terms of the bar vector b by t = V 2 M 1 V T 2 ST 1 b ; s = Γ 1 t ; f = Λb [63 Construction Procedure One procedure for construction of a class 1 tensegrity structure in equilibrium is provided as follows : 1) Given the matrices S and B from the network topology, find a nonsingular matrix P = [P T 1,P T 2 such that B T q = B T P = [I nb,0 nb, and calculate {S 1 = P 1 S,S 2 = P 2 S}. (Note that if B is defined as in (19), then {P 1,P 2 } are given by (21).) 2) Choose {γ m > 0} such that X 0 in (55) is diagonal, and then select {λ n > 0} such that Λ = X. 3) Select any bar vector b such that {b 1,b 2,...,b nb } span R 3, and calculate {t,s,f} from (63). 4) Choose the free vector z e in (59), (60) to give a suitable node vector p.

17 Equilibrium Conditions 307 In the design of tensegrity structures, an alternative approach for the selection of the node vector p and bar vector b may be useful. That is, we can first select the node vector p as in (59), (60) (which satisfies the constraint (61)) subject (if necessary) to other design constraints. The corresponding bar vector b is then determined from the first 3n b components of the vector P 1 p Computation of Equilibrium Configurations An inspection of all matrices in all algebraic expressions reveals that each matrix A = [a ij is a tensor product of the form A = à I = [a ij I where à is m n and the identity matrix I is of either dimension 3n 3m, or 2n 2m for planar tensegrity structures. Hence all matrix manipulations with the exception of the final determination of the vectors {p, b, s, t, f} can be carried out in lower dimension thereby reducing the computational burden. In particular, the sufficient condition (58) for the existence of a prestress equilibrium can be written in the form where Λ = X ; λ = Λ I, X = X I X = ( S 1 Ṽ 2 ) M 1 ( S 1 Ṽ 2 ) T ; M = Ṽ T 2 Λ 1 Ṽ 2 We now illustrate the construction procedure Example 3 For the 2-bar 4-string planar tensegrity structure, we can replace {S,B,S 1, S 2,P 1,P 2 } in examples 1 and 2 by { S, B, S 1, S 2, P 1, P 2 } where S = S 1 = P 1 = ; B = [ [ ; S = [ [ ; P2 = Then the singular value decomposition {Ũ, Σ,Ṽ} of the (reduced dimension) matrix S 2 of (reduced) rank r = ρ( S 2 ) = 1 is given by Σ 11 = 2 2 and

18 308 RFGC - 7/2003. Tenségrité and Ũ 1 = 1 2 [ 1 1 ; Ũ 2 = 1 2 [ 1 1 ; Ṽ 1 = Ṽ 2 = S 1 Ṽ 2 = [ We also have that ( S 1 Ṽ 2 )( S 1 Ṽ 2 ) T = I 2. (i) Suppose {γ k = 1; k = 1, 2, 3, 4} so that the (reduced) matrix Γ = I 4. It then follows that the (reduced) matrix X = I 4 which satisfies condition (58) in Theorem 1. Hence when there is no external load (i.e. w = 0), (54) has a (nonzero) solution b T = [b T 1,bT 2 when λ 1 = λ 2 = 1. For example, when b 1 = [4, 0,b 2 = [0, 4, the nodes {p 1 = [ 3, 1,p 2 = [1, 1,p 3 = [ 1, 3,p 4 = [ 1, 1} define an equilibrium solution. (ii) We hypothesize that there is no equilibrium solution unless : for some γ, {γ n = γ} for 1 n 4 and {λ m = γ} for 1 m 2. Such tensegrity structures may be said to have force symmetry (but not necessarily geometric symmetry). No equilibrium solutions have been found for other choices of {γ n, λ n }. An analytical investigation was undertaken for the case {γ 1 = γ 3 = 1, γ 2 = γ 4 = a}. In this case, it may be shown that X = 0.5 [ a + 1 a 1 a 1 a + 1 with singular values (and eigenvalues) given by {1, a}. In all choices for {λ n } that led to the rank of Λ X having rank 2, the 4 2 matrix Ṽb2 is of the form ṼT b2 = [ÃT 1, ±ÃT 1 ; that is, in (57), the two bar vectors {b 1,b 2 } are always parallel, so the equilibrium structure is trivial with {p 1 = p 3,p 2 = p 4 }. Example 4 A (3; 9; 3) class 1 tensegrity consisting of one stage of n b = 3 bars (so n p = 6 nodes) with n s = 9 strings has topology as illustrated in Fig. 2 The corresponding (reduced) geometric matrices { S, B} are given by

19 Equilibrium Conditions 309 Figure 2. 3-Bar 9-String (3; 9; 3) Class 1 Tensegrity S = ; B = A symmetrical force configuration was investigated with equal base string force coefficients {γ 1 = γ 2 = γ 1 = γbase }, equal top string coefficients {γ 5 = γ 6 = γ 7 = γ top }, and equal vertical string coefficients {γ 4 = γ 8 = γ 9 = γvert }. The following conclusions were reached : (i) No equilibrium solution is possible when γ base = γ top = γ vert. In this case X is not diagonal, and once again, all choices of Λ such that ρ( Λ X) 0 led to parallel bar vectors. (ii) For all choices of {γ base, γ top }, an equilibrium solution with Λ = X is provided by λ 1 = λ 2 = λ 3 = λ where γ vert = λ. For example, when {γ base = 1, γ top = 1}, then γ vert = λ = , and when {γ base = 1, γ top = 5}, then γ vert = λ = In all cases, X is a scaled identity matrix. Recall that once the condition Λ = X is satisfied, the choice for the bar vectors was completely arbitrary subject only to {b 1,b 2,b 3 } spanning R 3. Alternatively, as required in (61), the construction may be achieved subject to the node vector p being in the null space of the matrix U T Q2 P 1. It then follows that parameters α 1 = γ base γ top ; α 2 = λ γ top

20 310 RFGC - 7/2003. Tenségrité may be used to optimize the node positions with respect to some performance measure (such as being restricted to lie on a given surface). 4. Conclusion This paper characterizes the static equilibria of a class 1 tensegrity structure. Analytical expressions are derived for the equilibrium condition of a tensegrity structure in terms of member force coefficients and string and bar connectivity information. We use vectors to describe each element (bars and tendons), eliminating the need to use direction cosines and the subsequent transcendental functions that follow their use. By enlarging the vector space in which we characterize the problem, the mathematical structure of the equations admit treatment by linear algebra methods, for the most part. This reduces the study of a significant portion of the tensegrity equilibria to a series of linear algebra problems. Our results characterize the equilibria conditions of tensegrity structures in terms of a very small number of variables since the necessary and sufficient conditions of the linear algebra treatment has eliminated several of the original variables. This formulation offers insight and identifies the free parameters that may be used to achieve desired structural shapes. Since all conditions are necessary and sufficient, these results can be used in the design of any tensegrity structure. 5. References R. Connelly, Rigidity and Energy, Inventiones Mathematicae, 66 (1), pp , R. Connelly, Rigidity in : P.M. Gruber and J.M. Wills, ed., Handbook of Convex Geometry, Elsevier Publishers Ltd, pp , R.E. Skelton and C. Sultan, Controllable Tensegrity, Proc. SPIE 4th Ann. Symp. Smart Str. & Mats., San Diego, March R.E. Skelton, W.J. Helton and R. Adhikari, Mechanics of Tensegrity Beams, UCSD, Structural Systems & Contr. Lab., Rep. No , R.E. Skelton, J.P. Pinaud, and D.L. Mingori, Dynamics of the shell class of tensegrity structures, Journal of the Franklin Institute, 338/2-3, pp , D. Williamson and R.E. Skelton, A General Class of Tensegrity Systems : Geometric Definition, Engineering Mechanics for the 21st Century, ASCE Conference, 1998 R. Motro, Tensegrity Systems : The State of the Art, Int. J. Space Structures, 7(2) : pp , D. E. Ingber, Architecture of Life, Scientific American, Jan., pp , Tibert and S. Pellegrino, Review of Form-Finding Methods for Tensegrity Strucutres, submitted to Int. J. Space Structures