Appendix A: Matrices and Tensors
Total Page:16
File Type:pdf, Size:1020Kb
Appendix A: Matrices and Tensors A.1 Introduction and Rationale The purpose of this appendix is to present the notation and most of the mathematical techniques that will be used in the body of the text. The audience is assumed to have been through several years of college level mathematics that included the differential and integral calculus, differential equations, functions of several variables, partial derivatives, and an introduction to linear algebra. Matrices are reviewed briefly and determinants, vectors, and tensors of order two are described. The application of this linear algebra to material that appears in undergraduate engineering courses on mechanics is illustrated by discussions of concepts like the area and mass moments of inertia, Mohr’s circles and the vector cross and triple scalar products. The solutions to ordinary differential equations are reviewed in the last two sections. The notation, as far as possible, will be a matrix notation that is easily entered into existing symbolic computational programs like Maple, Mathematica, Matlab, and Mathcad etc. The desire to represent the components of three-dimensional fourth order tensors that appear in anisotropic elasticity as the components of six-dimensional second order tensors and thus represent these components in matrices of tensor components in six dimensions leads to the nontraditional part of this appendix. This is also one of the nontraditional aspects in the text of the book, but a minor one. This is described in Sect. A.11, along with the rationale for this approach. A.2 Definition of Square, Column, and Row Matrices An r by c matrix M is a rectangular array of numbers consisting of r rows and c columns, S.C. Cowin, Continuum Mechanics of Anisotropic Materials, 335 DOI 10.1007/978-1-4614-5025-2, # Springer Science+Business Media New York 2013 336 Appendix A: Matrices and Tensors 2 3 M11 M12 :::M1c 6 ::: 7 6 M21 M22 M2c 7: M ¼ 4 ::::::5 (A.1) Mr1 ::::Mrc The typical element of the array, Mij, is the ith element in the jth column; in this text the elements Mij will be real numbers or functions whose values are real numbers. The transpose of the matrix M is denoted by MT and is obtained from M by interchanging the rows and columns 2 3 M11 M21 :::Mr1 6 ::: 7 T 6 M12 M22 Mr2 7: M ¼ 4 ::::::5 (A.2) M1c ::::Mrc The operation of obtaining MT from M is called transposition. In this text we are interested in special cases of the r by c matrix M. These special cases are those of the square matrix, r ¼ c ¼ n, the case of the row matrix, r ¼ 1, c ¼ n, and the case of column matrix, r ¼ n, c ¼ 1. Further, the special sub-cases of interest are n ¼ 2, n ¼ 3, and n ¼ 6; the sub-case n ¼ 1 reduces all three special cases to the trivial situation of a single number or scalar. A square matrix A has the form 2 3 A11 A12 :::A1n 6 ::: 7 6 A21 A22 A2n 7; A ¼ 4 ::::::5 (A.3) An1 ::::Ann while row and column matrices, r and c, have the forms 2 3 c1 6 7 6 c2 7 6 7 6 : 7 r ¼½r1 r2 :::rn ; c ¼ 6 7; (A.4) 6 : 7 4 : 5 cn respectively. The transpose of a column matrix is a row matrix, thus T c ¼½c1 c2 :::cn : (A.5) To save space in books and papers the form of c in (A.5) is used more frequently than the form in the second of (A.4). Wherever possible, square matrices will be denoted by upper case boldface Latin letters, while row and column matrices will be denoted by lower case boldface Latin letters as is the case in equations (A.3) and (A.4). A.3 The Types and Algebra of Square Matrices 337 A.3 The Types and Algebra of Square Matrices The elements of the square matrix A given by (A.3) for which the row and column indices are equal, namely the elements A11, A22,..., Ann, are called diagonal elements. The sum of the diagonal elements of a matrix is a scalar called the trace of the matrix and, for a matrix A, it is denoted by trA, trA ¼ A11 þ A22 þÁÁÁþAnn: (A.6) If the trace of a matrix is zero, the matrix is said to be traceless. Note also that trA ¼ trAT. A matrix with only diagonal elements is called a diagonal matrix, 2 3 A11 0 ::: 0 6 ::: 7 6 0 A22 0 7: A ¼ 4 ::::::5 (A.7) 0 ::::Ann The zero and the unit matrix, 0 and 1, respectively, constitute the null element, the 0, and the unit element, the 1, in the algebra of square matrices. The zero matrix is a matrix whose every element is zero and the unit matrix is a diagonal matrix whose diagonal elements are all one: 2 3 2 3 00:::0 10:::0 6 ::: 7 6 ::: 7 6 00 0 7; 6 01 0 7: 0 ¼ 4 : :::::5 1 ¼ 4 ::::::5 (A.8) 0 ::::0 0 ::::1 A special symbol, the Kronecker delta dij, is introduced to represent the components of the unit matrix. When the indices are equal, i ¼ j, the value of the Kronecker delta is one, d11 ¼ d22 ¼ ... ¼ dnn ¼ 1 and when they are unequal, i 6¼ j, the value of the Kronecker delta is zero, d12 ¼ d21 ¼ ... ¼ dn1 ¼ d1n ¼ 0. The multiplication of a matrix A by a scalar is defined as the multiplication of every element of the matrix A by the scalar a, thus 2 3 aA11 aA12 :::aA1n 6 ::: 7 6 aA21 aA22 aA2n 7: aA ¼ 4 ::::::5 (A.9) aAn1 ::::aAnn It is then easy to show that 1A ¼ A, À1A ¼A,0A ¼ 0, and a0 ¼ 0. The addition of matrices is defined only for matrices with the same number of rows and columns. The sum of two matrices, A and B, is denoted by A þ B, where 338 Appendix A: Matrices and Tensors 2 3 A11 þ B11 A12 þ B12 :::A1n þ B1n 6 ::: 7 6 A21 þ B21 A22 þ B22 A2n þ B2n 7: A þ B ¼ 4 ::::::5 (A.10) An1 þ Bn1 ::::Ann þ Bnn Matrix addition is commutative and associative, A þ B ¼ B þ A and A þðB þ CÞ¼ðA þ BÞþC; (A.11) respectively. The following distributive laws connect matrix addition and matrix multiplication by scalars: aðA þ BÞ¼aA þ aB and ða þ bÞA ¼ aA þ bA; (A.12) where a and b are scalars. Negative square matrices may be created by employing the definition of matrix multiplication by a scalar (A.8) in the special case when a ¼1. In this case the definition of addition of square matrices (A.10) can be extended to include the subtraction of square matrices, AÀB. A matrix for which B ¼ BT is said to be a symmetric matrix and a matrix for which C ¼CT is said to be a skew-symmetric or anti-symmetric matrix. The symmetric and anti-symmetric parts of a matrix, say AS and AA, are constructed from A as follows: 1 symmetric part of AS ¼ ðA þ ATÞ; and (A.13) 2 1 anti À symmetric part of AA ¼ ðA À ATÞ: (A.14) 2 It is easy to verify that the symmetric part of A is a symmetric matrix and that the skew-symmetric part of A is a skew-symmetric matrix. The sum of the symmetric part of A and the skew-symmetric part of A, AS þ AA,isA: 1 1 A ¼ AS þ AA ¼ ðA þ ATÞþ ðA À ATÞ: (A.15) 2 2 This result shows that any square matrix can be decomposed into the sum of a symmetric and a skew-symmetric matrix or anti-symmetric matrix. Using the trace operation introduced above, the representation (A.15) can be extended to three-way decomposition of the matrix A, trA 1 trA 1 A ¼ 1 þ A þ AT À 2 1 þ ðA À ATÞ: (A.16) n 2 n 2 A.3 The Types and Algebra of Square Matrices 339 The last term in this decomposition is still the skew-symmetric part of the matrix. The second term is the traceless symmetric part of the matrix and the first term is simply the trace of the matrix multiplied by the unit matrix. Example A.3.1 Construct the three-way decomposition of the matrix A given by: 2 3 123 A ¼ 4 4565: 789 Solution: The symmetric and skew-symmetric parts of A, AS, and AA, as well as the trace of A are calculated, 2 3 2 3 135 0 À1 À2 1 6 7 1 6 7 AS ¼ ðA þ ATÞ¼4 3575; AA ¼ ðA À ATÞ¼4 10À1 5; trA ¼ 15; 2 2 579 21 0 then, since n ¼ 3, it follows from (A.16) that 2 3 2 3 2 3 500 À435 0 À1 À2 A ¼ 4 0505 þ 4 3075 þ 4 10À1 5: 005 574 21 0 Introducing the notation for the deviatoric part of an n by n square matrix A, trA devA ¼ A À 1; (A.17) n the representation for the matrix A given by (A.16) may be rewritten as the sum of three terms trA A ¼ 1 þ devAS þ AA; (A.18) n where the first term is called the isotropic,orspherical (or hydrostatic as in hydraulics) part of A. Note that 1 devAS ¼ ðdevA þ devATÞ: (A.19) 2 Example A.3.2 Show that tr (devA) ¼ 0.