Equilibrium Conditions of Class 1 Tensegrity Structures
|
|
- Sydney James
- 6 years ago
- Views:
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
Equilibrium Conditions of a Tensegrity Structure
Proceedings of the Third World Conference on Structural Control (3WCSC) Como, Italy, April 7-12, 2002 Equilibrium Conditions of a Tensegrity Structure Darrell Williamson, Robert E. Skelton, and JeongHeon
More informationVolume 112, Advances in the Astronautical Sciences Spaceflight Mechanics 2002
Volume 112, Advances in the Astronautical Sciences Spaceflight Mechanics 2002 AAS 02-177 EQUILIBRIUM CONDITIONS OF A CLASS I TENSEGRITY STRUCTURE Robert E. Skelton 1, Darrell Williamson 2, JeongHeon Han
More informationPrestress stability. Lecture VI. Session on Granular Matter Institut Henri Poincaré. R. Connelly Cornell University Department of Mathematics
Prestress stability Lecture VI Session on Granular Matter Institut Henri Poincaré R. Connelly Cornell University Department of Mathematics 1 Potential functions How is the stability of a structure determined
More informationEconomic sensor/actuator selection and its application to flexible structure control
To appear in Proc. SPIE Int. Soc. Opt. Eng. 2004 Economic sensor/actuator selection and its application to flexible structure control Robert E. Skelton a and Faming Li a a Department of Mechanical and
More informationOPTIMUM PRE-STRESS DESIGN FOR FREQUENCY REQUIREMENT OF TENSEGRITY STRUCTURES
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm OPTIMUM PRE-STRESS DESIGN FOR FREQUENCY REQUIREMENT OF TENSEGRITY STRUCTURES Seif Dalil
More information(a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? Solution: dim N(A) 1, since rank(a) 3. Ax =
. (5 points) (a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? dim N(A), since rank(a) 3. (b) If we also know that Ax = has no solution, what do we know about the rank of A? C(A)
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationLecture notes: Applied linear algebra Part 1. Version 2
Lecture notes: Applied linear algebra Part 1. Version 2 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 1 Notation, basic notions and facts 1.1 Subspaces, range and
More informationDynamics and control of tensegrity systems
Dynamics and control of tensegrity systems Robert Skelton Feb 24, 24 Abstract This paper introduces rigid body dynamics in a new form, as a matrix differential equation, rather than the traditional vector
More informationSimon D. GUEST Reader University of Cambridge Cambridge, UK.
Multi-stable Star-shaped Tensegrity Structures Jingyao ZHANG Lecturer Ritsumeikan University Kusatsu, Shiga, JAPAN zhang@fc.ritsumei.ac.jp Robert CONNELLY Professor Cornell University Ithaca, NY, USA rc46@cornell.edu
More informationA note on an unusual type of polar decomposition
A note on an unusual type of polar decomposition H. Faßbender a, Kh. D. Ikramov b,,1 a Institut Computational Mathematics, TU Braunschweig, Pockelsstr. 14, D-38023 Braunschweig, Germany. b Faculty of Computational
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationSemidefinite Programming
Semidefinite Programming Notes by Bernd Sturmfels for the lecture on June 26, 208, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra The transition from linear algebra to nonlinear algebra has
More informationChapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors
Chapter 7 Canonical Forms 7.1 Eigenvalues and Eigenvectors Definition 7.1.1. Let V be a vector space over the field F and let T be a linear operator on V. An eigenvalue of T is a scalar λ F such that there
More informationTopology Optimization of Tensegrity Structures Based on Nonsmooth Mechanics. Yoshihiro Kanno. November 14, 2011 ACOMEN 2011
Topology Optimization of Tensegrity Structures Based on Nonsmooth Mechanics Yoshihiro Kanno November 14, 2011 ACOMEN 2011 tensegrity definition tension + integrity [Fuller 75] [Emmerich], [Snelson] pin-jointed
More information7. Symmetric Matrices and Quadratic Forms
Linear Algebra 7. Symmetric Matrices and Quadratic Forms CSIE NCU 1 7. Symmetric Matrices and Quadratic Forms 7.1 Diagonalization of symmetric matrices 2 7.2 Quadratic forms.. 9 7.4 The singular value
More informationChapter 1: Systems of Linear Equations
Chapter : Systems of Linear Equations February, 9 Systems of linear equations Linear systems Lecture A linear equation in variables x, x,, x n is an equation of the form a x + a x + + a n x n = b, where
More informationEIGENVALUES AND EIGENVECTORS 3
EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices
More informationMTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education
MTH 3 Linear Algebra Study Guide Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education June 3, ii Contents Table of Contents iii Matrix Algebra. Real Life
More informationSHAPE CONTROL OF A MULTI-AGENT SYSTEM USING TENSEGRITY STRUCTURES
SHAPE CONTROL OF A MULTI-AGENT SYSTEM USING TENSEGRITY STRUCTURES Benjamin Nabet,1 Naomi Ehrich Leonard,1 Mechanical and Aerospace Engineering, Princeton University, Princeton, NJ 08544 USA, {bnabet, naomi}@princeton.edu
More informationBistable Regimes in an Elastic Tensegrity System
Bistable Regimes in an Elastic Tensegrity System Andrea Micheletti Department of Civil Engineering and Computer Science Engineering University of Rome Tor Vergata Via Politecnico 1, 00133, Rome, Italy
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2017 LECTURE 5
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 17 LECTURE 5 1 existence of svd Theorem 1 (Existence of SVD) Every matrix has a singular value decomposition (condensed version) Proof Let A C m n and for simplicity
More informationA MARCHING PROCEDURE FOR FORM-FINDING FOR TENSEGRITY STRUCTURES
JOURNAL OF MECHANICS OF MATERIALS AND STRUCTURES Vol., No., A MARCHING PROCEDURE FOR FORM-FINDING FOR TENSEGRITY STRUCTURES ANDREA MICHELETTI AND WILLIAM O. WILLIAMS We give an algorithm for solving the
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More informationSingular Value Decomposition
Singular Value Decomposition CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Singular Value Decomposition 1 / 35 Understanding
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
More informationLecture Summaries for Linear Algebra M51A
These lecture summaries may also be viewed online by clicking the L icon at the top right of any lecture screen. Lecture Summaries for Linear Algebra M51A refers to the section in the textbook. Lecture
More informationChapter 12. Static Equilibrium and Elasticity
Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial
More informationMath 408 Advanced Linear Algebra
Math 408 Advanced Linear Algebra Chi-Kwong Li Chapter 4 Hermitian and symmetric matrices Basic properties Theorem Let A M n. The following are equivalent. Remark (a) A is Hermitian, i.e., A = A. (b) x
More informationChapter 2: Matrix Algebra
Chapter 2: Matrix Algebra (Last Updated: October 12, 2016) These notes are derived primarily from Linear Algebra and its applications by David Lay (4ed). Write A = 1. Matrix operations [a 1 a n. Then entry
More informationLinear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations.
POLI 7 - Mathematical and Statistical Foundations Prof S Saiegh Fall Lecture Notes - Class 4 October 4, Linear Algebra The analysis of many models in the social sciences reduces to the study of systems
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationDiscrete Optimization 23
Discrete Optimization 23 2 Total Unimodularity (TU) and Its Applications In this section we will discuss the total unimodularity theory and its applications to flows in networks. 2.1 Total Unimodularity:
More informationKinestatic Analyses of Mechanisms with Compliant Elements
Kinestatic Analyses of Mechanisms with Compliant Elements Carl Crane How can such a simple mechanism have such a high order solution? Tensegrity structures comprised of struts in compression and ties in
More informationThe Singular Value Decomposition (SVD) and Principal Component Analysis (PCA)
Chapter 5 The Singular Value Decomposition (SVD) and Principal Component Analysis (PCA) 5.1 Basics of SVD 5.1.1 Review of Key Concepts We review some key definitions and results about matrices that will
More informationSolvability of Linear Matrix Equations in a Symmetric Matrix Variable
Solvability of Linear Matrix Equations in a Symmetric Matrix Variable Maurcio C. de Oliveira J. William Helton Abstract We study the solvability of generalized linear matrix equations of the Lyapunov type
More informationMatrices and Vectors. Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A =
30 MATHEMATICS REVIEW G A.1.1 Matrices and Vectors Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A = a 11 a 12... a 1N a 21 a 22... a 2N...... a M1 a M2... a MN A matrix can
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 2
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 2 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory April 5, 2012 Andre Tkacenko
More informationSingular Value Decomposition
Singular Value Decomposition Motivatation The diagonalization theorem play a part in many interesting applications. Unfortunately not all matrices can be factored as A = PDP However a factorization A =
More informationMA 527 first midterm review problems Hopefully final version as of October 2nd
MA 57 first midterm review problems Hopefully final version as of October nd The first midterm will be on Wednesday, October 4th, from 8 to 9 pm, in MTHW 0. It will cover all the material from the classes
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationMATH 583A REVIEW SESSION #1
MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationEIGENVALUE PROBLEMS. Background on eigenvalues/ eigenvectors / decompositions. Perturbation analysis, condition numbers..
EIGENVALUE PROBLEMS Background on eigenvalues/ eigenvectors / decompositions Perturbation analysis, condition numbers.. Power method The QR algorithm Practical QR algorithms: use of Hessenberg form and
More informationLinear Systems. Carlo Tomasi. June 12, r = rank(a) b range(a) n r solutions
Linear Systems Carlo Tomasi June, 08 Section characterizes the existence and multiplicity of the solutions of a linear system in terms of the four fundamental spaces associated with the system s matrix
More informationSymmetric Prismatic Tensegrity Structures: Part II. Symmetry-adapted Formulations
Symmetric Prismatic Tensegrity Structures: Part II. Symmetry-adapted Formulations J.Y. Zhang a S.D. Guest b, M. Ohsaki a a Dept. of Architecture & Architectural Engineering, Kyoto University, Japan b Dept.
More informationAdvanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One
Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to
More informationLecture Note 12: The Eigenvalue Problem
MATH 5330: Computational Methods of Linear Algebra Lecture Note 12: The Eigenvalue Problem 1 Theoretical Background Xianyi Zeng Department of Mathematical Sciences, UTEP The eigenvalue problem is a classical
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationNotes on singular value decomposition for Math 54. Recall that if A is a symmetric n n matrix, then A has real eigenvalues A = P DP 1 A = P DP T.
Notes on singular value decomposition for Math 54 Recall that if A is a symmetric n n matrix, then A has real eigenvalues λ 1,, λ n (possibly repeated), and R n has an orthonormal basis v 1,, v n, where
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More information18.06 Problem Set 10 - Solutions Due Thursday, 29 November 2007 at 4 pm in
86 Problem Set - Solutions Due Thursday, 29 November 27 at 4 pm in 2-6 Problem : (5=5+5+5) Take any matrix A of the form A = B H CB, where B has full column rank and C is Hermitian and positive-definite
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationMathematical Properties of Stiffness Matrices
Mathematical Properties of Stiffness Matrices CEE 4L. Matrix Structural Analysis Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall, 0 These notes describe some of the
More informationTruss Structures: The Direct Stiffness Method
. Truss Structures: The Companies, CHAPTER Truss Structures: The Direct Stiffness Method. INTRODUCTION The simple line elements discussed in Chapter introduced the concepts of nodes, nodal displacements,
More informationMATH 2331 Linear Algebra. Section 2.1 Matrix Operations. Definition: A : m n, B : n p. Example: Compute AB, if possible.
MATH 2331 Linear Algebra Section 2.1 Matrix Operations Definition: A : m n, B : n p ( 1 2 p ) ( 1 2 p ) AB = A b b b = Ab Ab Ab Example: Compute AB, if possible. 1 Row-column rule: i-j-th entry of AB:
More informationLinear Systems. Carlo Tomasi
Linear Systems Carlo Tomasi Section 1 characterizes the existence and multiplicity of the solutions of a linear system in terms of the four fundamental spaces associated with the system s matrix and of
More informationChapter 7: Symmetric Matrices and Quadratic Forms
Chapter 7: Symmetric Matrices and Quadratic Forms (Last Updated: December, 06) These notes are derived primarily from Linear Algebra and its applications by David Lay (4ed). A few theorems have been moved
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationRigidity of Skeletal Structures
CHAPTER 3 Rigidity of Skeletal Structures 3. INTRODUCTION The rigidity of structures has been studied by pioneering structural engineers such as Henneberg [79] and Müller-Breslau [76]. The methods they
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationThe Singular Value Decomposition
The Singular Value Decomposition Philippe B. Laval KSU Fall 2015 Philippe B. Laval (KSU) SVD Fall 2015 1 / 13 Review of Key Concepts We review some key definitions and results about matrices that will
More informationSolving Homogeneous Systems with Sub-matrices
Pure Mathematical Sciences, Vol 7, 218, no 1, 11-18 HIKARI Ltd, wwwm-hikaricom https://doiorg/112988/pms218843 Solving Homogeneous Systems with Sub-matrices Massoud Malek Mathematics, California State
More informationDefinition 2.3. We define addition and multiplication of matrices as follows.
14 Chapter 2 Matrices In this chapter, we review matrix algebra from Linear Algebra I, consider row and column operations on matrices, and define the rank of a matrix. Along the way prove that the row
More informationLinear Algebra. Linear Equations and Matrices. Copyright 2005, W.R. Winfrey
Copyright 2005, W.R. Winfrey Topics Preliminaries Systems of Linear Equations Matrices Algebraic Properties of Matrix Operations Special Types of Matrices and Partitioned Matrices Matrix Transformations
More information8. Diagonalization.
8. Diagonalization 8.1. Matrix Representations of Linear Transformations Matrix of A Linear Operator with Respect to A Basis We know that every linear transformation T: R n R m has an associated standard
More informationLecture 7: Positive Semidefinite Matrices
Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.
More informationLinear Algebra- Final Exam Review
Linear Algebra- Final Exam Review. Let A be invertible. Show that, if v, v, v 3 are linearly independent vectors, so are Av, Av, Av 3. NOTE: It should be clear from your answer that you know the definition.
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationStudy Guide for Linear Algebra Exam 2
Study Guide for Linear Algebra Exam 2 Term Vector Space Definition A Vector Space is a nonempty set V of objects, on which are defined two operations, called addition and multiplication by scalars (real
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationNotes on Eigenvalues, Singular Values and QR
Notes on Eigenvalues, Singular Values and QR Michael Overton, Numerical Computing, Spring 2017 March 30, 2017 1 Eigenvalues Everyone who has studied linear algebra knows the definition: given a square
More informationON THE SINGULAR DECOMPOSITION OF MATRICES
An. Şt. Univ. Ovidius Constanţa Vol. 8, 00, 55 6 ON THE SINGULAR DECOMPOSITION OF MATRICES Alina PETRESCU-NIŢǍ Abstract This paper is an original presentation of the algorithm of the singular decomposition
More informationCHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS. W. Erwin Diewert January 31, 2008.
1 ECONOMICS 594: LECTURE NOTES CHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS W. Erwin Diewert January 31, 2008. 1. Introduction Many economic problems have the following structure: (i) a linear function
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5
Semidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5 Instructor: Farid Alizadeh Scribe: Anton Riabov 10/08/2001 1 Overview We continue studying the maximum eigenvalue SDP, and generalize
More informationChapter 1. Vectors, Matrices, and Linear Spaces
1.6 Homogeneous Systems, Subspaces and Bases 1 Chapter 1. Vectors, Matrices, and Linear Spaces 1.6. Homogeneous Systems, Subspaces and Bases Note. In this section we explore the structure of the solution
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationLeast Sparsity of p-norm based Optimization Problems with p > 1
Least Sparsity of p-norm based Optimization Problems with p > Jinglai Shen and Seyedahmad Mousavi Original version: July, 07; Revision: February, 08 Abstract Motivated by l p -optimization arising from
More informationMATH36001 Generalized Inverses and the SVD 2015
MATH36001 Generalized Inverses and the SVD 201 1 Generalized Inverses of Matrices A matrix has an inverse only if it is square and nonsingular. However there are theoretical and practical applications
More informationGeometric Mapping Properties of Semipositive Matrices
Geometric Mapping Properties of Semipositive Matrices M. J. Tsatsomeros Mathematics Department Washington State University Pullman, WA 99164 (tsat@wsu.edu) July 14, 2015 Abstract Semipositive matrices
More information. = V c = V [x]v (5.1) c 1. c k
Chapter 5 Linear Algebra It can be argued that all of linear algebra can be understood using the four fundamental subspaces associated with a matrix Because they form the foundation on which we later work,
More informationOnline Exercises for Linear Algebra XM511
This document lists the online exercises for XM511. The section ( ) numbers refer to the textbook. TYPE I are True/False. Lecture 02 ( 1.1) Online Exercises for Linear Algebra XM511 1) The matrix [3 2
More information2. Linear algebra. matrices and vectors. linear equations. range and nullspace of matrices. function of vectors, gradient and Hessian
FE661 - Statistical Methods for Financial Engineering 2. Linear algebra Jitkomut Songsiri matrices and vectors linear equations range and nullspace of matrices function of vectors, gradient and Hessian
More informationsublinear time low-rank approximation of positive semidefinite matrices Cameron Musco (MIT) and David P. Woodru (CMU)
sublinear time low-rank approximation of positive semidefinite matrices Cameron Musco (MIT) and David P. Woodru (CMU) 0 overview Our Contributions: 1 overview Our Contributions: A near optimal low-rank
More informationProperties of Linear Transformations from R n to R m
Properties of Linear Transformations from R n to R m MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Topic Overview Relationship between the properties of a matrix transformation
More informationLEC 2: Principal Component Analysis (PCA) A First Dimensionality Reduction Approach
LEC 2: Principal Component Analysis (PCA) A First Dimensionality Reduction Approach Dr. Guangliang Chen February 9, 2016 Outline Introduction Review of linear algebra Matrix SVD PCA Motivation The digits
More informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More informationEXERCISE SET 5.1. = (kx + kx + k, ky + ky + k ) = (kx + kx + 1, ky + ky + 1) = ((k + )x + 1, (k + )y + 1)
EXERCISE SET 5. 6. The pair (, 2) is in the set but the pair ( )(, 2) = (, 2) is not because the first component is negative; hence Axiom 6 fails. Axiom 5 also fails. 8. Axioms, 2, 3, 6, 9, and are easily
More informationLecture 2 INF-MAT : , LU, symmetric LU, Positve (semi)definite, Cholesky, Semi-Cholesky
Lecture 2 INF-MAT 4350 2009: 7.1-7.6, LU, symmetric LU, Positve (semi)definite, Cholesky, Semi-Cholesky Tom Lyche and Michael Floater Centre of Mathematics for Applications, Department of Informatics,
More information1. Let m 1 and n 1 be two natural numbers such that m > n. Which of the following is/are true?
. Let m and n be two natural numbers such that m > n. Which of the following is/are true? (i) A linear system of m equations in n variables is always consistent. (ii) A linear system of n equations in
More informationMath 240 Calculus III
The Calculus III Summer 2015, Session II Wednesday, July 8, 2015 Agenda 1. of the determinant 2. determinants 3. of determinants What is the determinant? Yesterday: Ax = b has a unique solution when A
More informationLinear Algebra Exam 1 Spring 2007
Linear Algebra Exam 1 Spring 2007 March 15, 2007 Name: SOLUTION KEY (Total 55 points, plus 5 more for Pledged Assignment.) Honor Code Statement: Directions: Complete all problems. Justify all answers/solutions.
More informationStiffness Matrices, Spring and Bar Elements
CHAPTER Stiffness Matrices, Spring and Bar Elements. INTRODUCTION The primary characteristics of a finite element are embodied in the element stiffness matrix. For a structural finite element, the stiffness
More informationLinear Algebra M1 - FIB. Contents: 5. Matrices, systems of linear equations and determinants 6. Vector space 7. Linear maps 8.
Linear Algebra M1 - FIB Contents: 5 Matrices, systems of linear equations and determinants 6 Vector space 7 Linear maps 8 Diagonalization Anna de Mier Montserrat Maureso Dept Matemàtica Aplicada II Translation:
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More information3.1 Limit Analysis and Design of Structures Formulated as LP Problems
Linear Programming 3 Mathematical programming is concerned with the extremization of a function f defined over an n-dimensional design space R n and bounded by a set S in the design space. The set S may
More informationMaths for Signals and Systems Linear Algebra in Engineering
Maths for Signals and Systems Linear Algebra in Engineering Lectures 13 15, Tuesday 8 th and Friday 11 th November 016 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE
More informationNotes taken by Graham Taylor. January 22, 2005
CSC4 - Linear Programming and Combinatorial Optimization Lecture : Different forms of LP. The algebraic objects behind LP. Basic Feasible Solutions Notes taken by Graham Taylor January, 5 Summary: We first
More informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationChapter 7. Linear Algebra: Matrices, Vectors,
Chapter 7. Linear Algebra: Matrices, Vectors, Determinants. Linear Systems Linear algebra includes the theory and application of linear systems of equations, linear transformations, and eigenvalue problems.
More information