Chapter 5 the Orientation and Stress Tensors

Total Page:16

File Type:pdf, Size:1020Kb

Chapter 5 the Orientation and Stress Tensors Chapter 5 The Orientation and Stress Tensors Introduction The topic of tensors typically produces significant anxiety for students of structural geology. That is due, at least in part, to the fact that the term is studiously avoided until it is sprung on the student when introducing the topic of stress, thus conflating two difficult to grasp concepts. In our case, the previous chapters have already introduced the concept of vectors as a type of first order tensor so the X3 X′2 X′3 X2 X′1 X1 fundamental definition — an entity that can be transformed from one coordinate system to another, changing its components in predictable ways such that its fun- damental nature doesn’t change — is already understood on an intuitive level. In the case of our vector, it has the same magnitude and orientation on the page re- gardless of the orientation of the coordinate axes (Fig. 4.3). We are now fully pre- pared to tackle second order tensors, which will allow us to quantify the relations between different families of vectors. These second order tensors will take us to the very core of structural geology. We will still be doing just multiplication and addi- tions (mostly), but first, we’ll need to introduce some concepts about matrices. CHAPTER 5 ORIENTATION & STRESS TENSORS Matrices and Indicial Notation In the last chapter, we saw that the transformation of a vector could be de- scribed with three simple equations (4.4), repeated here: v1′ = a11v1 + a12v2 + a13v3 " v2′ = a21v1 + a22v2 + a23v3 (5.1) v3′ = a31v1 + a32v2 + a33v3 As was pointed out at the time, the subscripts in these equations, all of which refer to specific coordinate axes, vary in an extremely systematic and precise way. That is because these equations are an alternative way of writing down a matrix multi- plication: ⎡ ⎤ ⎛ ⎞ ⎡ ⎤ v1′ a11 a12 a13 v1 ⎢ ⎥ ⎜ ⎟ ⎢ ⎥ " v′ = av = ⎢ v2′ ⎥ = ⎜ a21 a22 a23 ⎟ ⎢ v2 ⎥ (5.2) ⎢ v′ ⎥ ⎜ a a a ⎟ ⎢ v ⎥ ⎣ 3 ⎦ ⎝ 31 32 33 ⎠ ⎣ 3 ⎦ Note that the subscripts of the matrix a vary in exactly the same way as the sub- scripts, or indices, of aij in (5.1). Equation (5.2) works because there are the same number of columns of a as there are rows of v. For this reason, you cannot write the left side of the equation as v′ = va because v has only one column but a has three rows. That is, matrix multiplication is non-commutative: the order of the multiplication matters. This will be an extremely important insight when we come to strain in a few more chapters; strain is mathematically represented as a matrix multiplication so the order in which strain is superimposed determines the final outcome. Matrix multiplication gives us another way to write the dot product of two vectors, where we use the transpose of one of the vectors (i.e., one of the vectors is “flipped”): ⎡ v ⎤ ⎢ 1 ⎥ " u i v = uvT = ⎡ u u u ⎤ v = u v + u v + u v (5.3) ⎣ 1 2 3 ⎦⎢ 2 ⎥ 1 1 2 2 3 3 ⎢ v ⎥ ⎣ 3 ⎦ MODERN STRUCTURAL PRACTICE "88 R. W. ALLMENDINGER © 2015-16 CHAPTER 5 ORIENTATION & STRESS TENSORS You can imagine that it gets pretty tedious to keep writing out the equations in (5.1) or the right side of (5.2) and the bold text on the left side of (5.2), known as matrix notation, is one type of shorthand. There is a second type of shorthand, known as the summation convention, that is more convenient, especially when it comes to implementing these equations in a programing language. Using the sum- mation convention, we would write Equations (5.1) and (5.2) as: " vi′ = aijvj (5.4) Let’s break this down and see how it works. As we have already seen i and j can vary from 1 to 3 in value because there are three axes to our Cartesian coordinate system. There is one i one each side of the equation and it is referred to as the free suffix; that means there will be three equations and, in each, i will have a constant value. The index j is know as the dummy suffix and appears twice on the right hand side, only, of Equation (5.4). Thus, the summation occurs with respect to j in each of the three equations as shown in the following equation: 3 " vi′ = ∑ aijvj = ai1v1 + ai2v2 + ai3v3 (5.5) j=1 The “recipe” for the summation convention takes a little practice but it turns out to be quite powerful and so is worth learning. We will encounter various properties of matrices in the subsequent chapters but for right now, we’ll wrap up with three important terms: and matrix is sym- metric if there are six independent values and the three off-diagonal components above the principal diagonal are the same as the three below the principal diagonal (Fig. 5.1). That is, T12 = T21, T13 = T31, and T23 = T . If a matrix is asymmetric then there are nine 32 T T T 11 12 13 off- independent components and T12 ≠ T21, T13 ≠ T31, diagonal T21 T22 T23 and T23 ≠ T32. Finally, a matrix can be antisym- metric with three independent components where T31 T32 T33 ( ) principal the values along the principal diagonal are all zeros off- diagonal diagonal and T12 = –T21, T13 = –T31, and T23 = –T32. Figure 5.1 — Anatomy of a 3 × 3 matrix. MODERN STRUCTURAL PRACTICE "89 R. W. ALLMENDINGER © 2015-16 CHAPTER 5 ORIENTATION & STRESS TENSORS Tensors In this Chapter, when we use the term “tensors”, we are specifically referring to second order or second rank tensors. All tensors can be expressed as matrices (a 3 × 3 matrix in the case of second order tensors) but not all matrices are ten- sors. In order to be a tensor, the matrix must represent an entity or quantity that transforms like a tensor; that is, the components must change in a logical and sys- tematic way during a coordinate transformation. Note that the transformation ma- trix, a, is itself not a tensor. A common convention is that tensors are written with brackets, [ ], whereas matrices that are not tensors are written with parentheses, ( ). We used this convention, for example, in Equation (5.2). Because all tensors can be written as matrices, they can be symmetric, asymmetric, etc. and any operation that can be performed on a matrix can likewise be performed on a tensor. Tensors as Linear Vector Operators The best way to think of a second order tensor is that it relates two families of vectors. In the case of stress, discussed later in this chapter, the stress tensor re- lates the stress vector (or traction) on a plane to the orientation of the plane. In the case of strain, as we shall see in a subsequent chapter, the displacement gradient tensor (one of many tensors related to strain!) relates the position vector of a point to the displacement of the point during a deformation. If we know the stress tensor, we can calculate the stress vector on a plane of any orientation within a body. You can imagine that this is hugely important in any study of earthquakes, induced seismicity, faulting, etc. Likewise, with the displacement gradient tensor, we can cal- culate how all points within a body are displaced as a function of position during a deformation. The basic way that this relationship is written is as follows: "ui = Tijvj (5.6a) Where u and v are generic vectors and T a generic tensor that relates the two vec- tors. This equation expands following the rules of summation or tensor notation: MODERN STRUCTURAL PRACTICE "90 R. W. ALLMENDINGER © 2015-16 CHAPTER 5 ORIENTATION & STRESS TENSORS ⎡ u ⎤ ⎡ T T T ⎤⎡ v ⎤ ⎢ 1 ⎥ ⎢ 11 12 13 ⎥⎢ 1 ⎥ " ⎢ u2 ⎥ = ⎢ T21 T22 T23 ⎥⎢ v2 ⎥ (5.6b) ⎢ u ⎥ ⎢ T T T ⎥⎢ v ⎥ ⎣ 3 ⎦ ⎣ 31 32 33 ⎦⎣ 3 ⎦ u1 = T11v1 + T12v2 + T13v3 " u2 = T21v1 + T22v2 + T23v3 (5.6c) u3 = T31v1 + T32v2 + T33v3 You already know how to do this expansion because it looks a lot like Equations (5.1-5.3). All of these equations represent the multiplication of a 3 × 3 and a 3 × 1 matrix. Beyond this mathematical similarity, there are entirely different. First, in Equation (5.1) it is the same vector, v, on both sides of the equation whereas u and v are two entirely different vectors in Equation (5.6). Secondly, T is a tensor as shown by the brackets in (5.6b), where as a in Equation (5.2) is not. There is a second way to make a second order tensor out of two vectors by taking the dyad (tensor) product of those two vectors. In indicial notation, we can write: T " T = u ⊗ v = u v or Tij = uivj (5.7a) The right hand side of (5.7) does not involve any summation because there is no dummy (i.e., repeated) suffix. So T simply works out to: ⎡ u v u v u v ⎤ ⎢ 1 1 1 2 1 3 ⎥ " Tij = ⎢ u2v1 u2v2 u2v3 ⎥ (5.7b) ⎢ u v u v u v ⎥ ⎣ 3 1 3 2 3 3 ⎦ The dyad product has important applications in earthquake and faulting studies and also in constructing the orientation tensor as we shall see below. Principal Axes of a Tensor Because our second order tensor is, well, a tensor, we can rotate the axes of the coordinate system and the values of the tensor will change. It turns out that, for MODERN STRUCTURAL PRACTICE "91 R. W. ALLMENDINGER © 2015-16 CHAPTER 5 ORIENTATION & STRESS TENSORS X3 X′2 X′3 Figure 5.2 — The magnitude ellip- soid of a tensor and two different coordinate axes through the ellip- soid.
Recommended publications
  • 21. Orthonormal Bases
    21. Orthonormal Bases The canonical/standard basis 011 001 001 B C B C B C B0C B1C B0C e1 = B.C ; e2 = B.C ; : : : ; en = B.C B.C B.C B.C @.A @.A @.A 0 0 1 has many useful properties. • Each of the standard basis vectors has unit length: q p T jjeijj = ei ei = ei ei = 1: • The standard basis vectors are orthogonal (in other words, at right angles or perpendicular). T ei ej = ei ej = 0 when i 6= j This is summarized by ( 1 i = j eT e = δ = ; i j ij 0 i 6= j where δij is the Kronecker delta. Notice that the Kronecker delta gives the entries of the identity matrix. Given column vectors v and w, we have seen that the dot product v w is the same as the matrix multiplication vT w. This is the inner product on n T R . We can also form the outer product vw , which gives a square matrix. 1 The outer product on the standard basis vectors is interesting. Set T Π1 = e1e1 011 B C B0C = B.C 1 0 ::: 0 B.C @.A 0 01 0 ::: 01 B C B0 0 ::: 0C = B. .C B. .C @. .A 0 0 ::: 0 . T Πn = enen 001 B C B0C = B.C 0 0 ::: 1 B.C @.A 1 00 0 ::: 01 B C B0 0 ::: 0C = B. .C B. .C @. .A 0 0 ::: 1 In short, Πi is the diagonal square matrix with a 1 in the ith diagonal position and zeros everywhere else.
    [Show full text]
  • Vectors, Matrices and Coordinate Transformations
    S. Widnall 16.07 Dynamics Fall 2009 Lecture notes based on J. Peraire Version 2.0 Lecture L3 - Vectors, Matrices and Coordinate Transformations By using vectors and defining appropriate operations between them, physical laws can often be written in a simple form. Since we will making extensive use of vectors in Dynamics, we will summarize some of their important properties. Vectors For our purposes we will think of a vector as a mathematical representation of a physical entity which has both magnitude and direction in a 3D space. Examples of physical vectors are forces, moments, and velocities. Geometrically, a vector can be represented as arrows. The length of the arrow represents its magnitude. Unless indicated otherwise, we shall assume that parallel translation does not change a vector, and we shall call the vectors satisfying this property, free vectors. Thus, two vectors are equal if and only if they are parallel, point in the same direction, and have equal length. Vectors are usually typed in boldface and scalar quantities appear in lightface italic type, e.g. the vector quantity A has magnitude, or modulus, A = |A|. In handwritten text, vectors are often expressed using the −→ arrow, or underbar notation, e.g. A , A. Vector Algebra Here, we introduce a few useful operations which are defined for free vectors. Multiplication by a scalar If we multiply a vector A by a scalar α, the result is a vector B = αA, which has magnitude B = |α|A. The vector B, is parallel to A and points in the same direction if α > 0.
    [Show full text]
  • Glossary of Linear Algebra Terms
    INNER PRODUCT SPACES AND THE GRAM-SCHMIDT PROCESS A. HAVENS 1. The Dot Product and Orthogonality 1.1. Review of the Dot Product. We first recall the notion of the dot product, which gives us a familiar example of an inner product structure on the real vector spaces Rn. This product is connected to the Euclidean geometry of Rn, via lengths and angles measured in Rn. Later, we will introduce inner product spaces in general, and use their structure to define general notions of length and angle on other vector spaces. Definition 1.1. The dot product of real n-vectors in the Euclidean vector space Rn is the scalar product · : Rn × Rn ! R given by the rule n n ! n X X X (u; v) = uiei; viei 7! uivi : i=1 i=1 i n Here BS := (e1;:::; en) is the standard basis of R . With respect to our conventions on basis and matrix multiplication, we may also express the dot product as the matrix-vector product 2 3 v1 6 7 t î ó 6 . 7 u v = u1 : : : un 6 . 7 : 4 5 vn It is a good exercise to verify the following proposition. Proposition 1.1. Let u; v; w 2 Rn be any real n-vectors, and s; t 2 R be any scalars. The Euclidean dot product (u; v) 7! u · v satisfies the following properties. (i:) The dot product is symmetric: u · v = v · u. (ii:) The dot product is bilinear: • (su) · v = s(u · v) = u · (sv), • (u + v) · w = u · w + v · w.
    [Show full text]
  • Chapter 5 ANGULAR MOMENTUM and ROTATIONS
    Chapter 5 ANGULAR MOMENTUM AND ROTATIONS In classical mechanics the total angular momentum L~ of an isolated system about any …xed point is conserved. The existence of a conserved vector L~ associated with such a system is itself a consequence of the fact that the associated Hamiltonian (or Lagrangian) is invariant under rotations, i.e., if the coordinates and momenta of the entire system are rotated “rigidly” about some point, the energy of the system is unchanged and, more importantly, is the same function of the dynamical variables as it was before the rotation. Such a circumstance would not apply, e.g., to a system lying in an externally imposed gravitational …eld pointing in some speci…c direction. Thus, the invariance of an isolated system under rotations ultimately arises from the fact that, in the absence of external …elds of this sort, space is isotropic; it behaves the same way in all directions. Not surprisingly, therefore, in quantum mechanics the individual Cartesian com- ponents Li of the total angular momentum operator L~ of an isolated system are also constants of the motion. The di¤erent components of L~ are not, however, compatible quantum observables. Indeed, as we will see the operators representing the components of angular momentum along di¤erent directions do not generally commute with one an- other. Thus, the vector operator L~ is not, strictly speaking, an observable, since it does not have a complete basis of eigenstates (which would have to be simultaneous eigenstates of all of its non-commuting components). This lack of commutivity often seems, at …rst encounter, as somewhat of a nuisance but, in fact, it intimately re‡ects the underlying structure of the three dimensional space in which we are immersed, and has its source in the fact that rotations in three dimensions about di¤erent axes do not commute with one another.
    [Show full text]
  • A Some Basic Rules of Tensor Calculus
    A Some Basic Rules of Tensor Calculus The tensor calculus is a powerful tool for the description of the fundamentals in con- tinuum mechanics and the derivation of the governing equations for applied prob- lems. In general, there are two possibilities for the representation of the tensors and the tensorial equations: – the direct (symbolic) notation and – the index (component) notation The direct notation operates with scalars, vectors and tensors as physical objects defined in the three dimensional space. A vector (first rank tensor) a is considered as a directed line segment rather than a triple of numbers (coordinates). A second rank tensor A is any finite sum of ordered vector pairs A = a b + ... +c d. The scalars, vectors and tensors are handled as invariant (independent⊗ from the choice⊗ of the coordinate system) objects. This is the reason for the use of the direct notation in the modern literature of mechanics and rheology, e.g. [29, 32, 49, 123, 131, 199, 246, 313, 334] among others. The index notation deals with components or coordinates of vectors and tensors. For a selected basis, e.g. gi, i = 1, 2, 3 one can write a = aig , A = aibj + ... + cidj g g i i ⊗ j Here the Einstein’s summation convention is used: in one expression the twice re- peated indices are summed up from 1 to 3, e.g. 3 3 k k ik ik a gk ∑ a gk, A bk ∑ A bk ≡ k=1 ≡ k=1 In the above examples k is a so-called dummy index. Within the index notation the basic operations with tensors are defined with respect to their coordinates, e.
    [Show full text]
  • Solving the Geodesic Equation
    Solving the Geodesic Equation Jeremy Atkins December 12, 2018 Abstract We find the general form of the geodesic equation and discuss the closed form relation to find Christoffel symbols. We then show how to use metric independence to find Killing vector fields, which allow us to solve the geodesic equation when there are helpful symmetries. We also discuss a more general way to find Killing vector fields, and some of their properties as a Lie algebra. 1 The Variational Method We will exploit the following variational principle to characterize motion in general relativity: The world line of a free test particle between two timelike separated points extremizes the proper time between them. where a test particle is one that is not a significant source of spacetime cur- vature, and a free particles is one that is only under the influence of curved spacetime. Similarly to classical Lagrangian mechanics, we can use this to de- duce the equations of motion for a metric. The proper time along a timeline worldline between point A and point B for the metric gµν is given by Z B Z B µ ν 1=2 τAB = dτ = (−gµν (x)dx dx ) (1) A A using the Einstein summation notation, and µ, ν = 0; 1; 2; 3. We can parame- terize the four coordinates with the parameter σ where σ = 0 at A and σ = 1 at B. This gives us the following equation for the proper time: Z 1 dxµ dxν 1=2 τAB = dσ −gµν (x) (2) 0 dσ dσ We can treat the integrand as a Lagrangian, dxµ dxν 1=2 L = −gµν (x) (3) dσ dσ and it's clear that the world lines extremizing proper time are those that satisfy the Euler-Lagrange equation: @L d @L − = 0 (4) @xµ dσ @(dxµ/dσ) 1 These four equations together give the equation for the worldline extremizing the proper time.
    [Show full text]
  • Geodetic Position Computations
    GEODETIC POSITION COMPUTATIONS E. J. KRAKIWSKY D. B. THOMSON February 1974 TECHNICALLECTURE NOTES REPORT NO.NO. 21739 PREFACE In order to make our extensive series of lecture notes more readily available, we have scanned the old master copies and produced electronic versions in Portable Document Format. The quality of the images varies depending on the quality of the originals. The images have not been converted to searchable text. GEODETIC POSITION COMPUTATIONS E.J. Krakiwsky D.B. Thomson Department of Geodesy and Geomatics Engineering University of New Brunswick P.O. Box 4400 Fredericton. N .B. Canada E3B5A3 February 197 4 Latest Reprinting December 1995 PREFACE The purpose of these notes is to give the theory and use of some methods of computing the geodetic positions of points on a reference ellipsoid and on the terrain. Justification for the first three sections o{ these lecture notes, which are concerned with the classical problem of "cCDputation of geodetic positions on the surface of an ellipsoid" is not easy to come by. It can onl.y be stated that the attempt has been to produce a self contained package , cont8.i.ning the complete development of same representative methods that exist in the literature. The last section is an introduction to three dimensional computation methods , and is offered as an alternative to the classical approach. Several problems, and their respective solutions, are presented. The approach t~en herein is to perform complete derivations, thus stqing awrq f'rcm the practice of giving a list of for11111lae to use in the solution of' a problem.
    [Show full text]
  • The Dot Product
    The Dot Product In this section, we will now concentrate on the vector operation called the dot product. The dot product of two vectors will produce a scalar instead of a vector as in the other operations that we examined in the previous section. The dot product is equal to the sum of the product of the horizontal components and the product of the vertical components. If v = a1 i + b1 j and w = a2 i + b2 j are vectors then their dot product is given by: v · w = a1 a2 + b1 b2 Properties of the Dot Product If u, v, and w are vectors and c is a scalar then: u · v = v · u u · (v + w) = u · v + u · w 0 · v = 0 v · v = || v || 2 (cu) · v = c(u · v) = u · (cv) Example 1: If v = 5i + 2j and w = 3i – 7j then find v · w. Solution: v · w = a1 a2 + b1 b2 v · w = (5)(3) + (2)(-7) v · w = 15 – 14 v · w = 1 Example 2: If u = –i + 3j, v = 7i – 4j and w = 2i + j then find (3u) · (v + w). Solution: Find 3u 3u = 3(–i + 3j) 3u = –3i + 9j Find v + w v + w = (7i – 4j) + (2i + j) v + w = (7 + 2) i + (–4 + 1) j v + w = 9i – 3j Example 2 (Continued): Find the dot product between (3u) and (v + w) (3u) · (v + w) = (–3i + 9j) · (9i – 3j) (3u) · (v + w) = (–3)(9) + (9)(-3) (3u) · (v + w) = –27 – 27 (3u) · (v + w) = –54 An alternate formula for the dot product is available by using the angle between the two vectors.
    [Show full text]
  • Concept of a Dyad and Dyadic: Consider Two Vectors a and B Dyad: It Consists of a Pair of Vectors a B for Two Vectors a a N D B
    1/11/2010 CHAPTER 1 Introductory Concepts • Elements of Vector Analysis • Newton’s Laws • Units • The basis of Newtonian Mechanics • D’Alembert’s Principle 1 Science of Mechanics: It is concerned with the motion of material bodies. • Bodies have different scales: Microscropic, macroscopic and astronomic scales. In mechanics - mostly macroscopic bodies are considered. • Speed of motion - serves as another important variable - small and high (approaching speed of light). 2 1 1/11/2010 • In Newtonian mechanics - study motion of bodies much bigger than particles at atomic scale, and moving at relative motions (speeds) much smaller than the speed of light. • Two general approaches: – Vectorial dynamics: uses Newton’s laws to write the equations of motion of a system, motion is described in physical coordinates and their derivatives; – Analytical dynamics: uses energy like quantities to define the equations of motion, uses the generalized coordinates to describe motion. 3 1.1 Vector Analysis: • Scalars, vectors, tensors: – Scalar: It is a quantity expressible by a single real number. Examples include: mass, time, temperature, energy, etc. – Vector: It is a quantity which needs both direction and magnitude for complete specification. – Actually (mathematically), it must also have certain transformation properties. 4 2 1/11/2010 These properties are: vector magnitude remains unchanged under rotation of axes. ex: force, moment of a force, velocity, acceleration, etc. – geometrically, vectors are shown or depicted as directed line segments of proper magnitude and direction. 5 e (unit vector) A A = A e – if we use a coordinate system, we define a basis set (iˆ , ˆj , k ˆ ): we can write A = Axi + Ay j + Azk Z or, we can also use the A three components and Y define X T {A}={Ax,Ay,Az} 6 3 1/11/2010 – The three components Ax , Ay , Az can be used as 3-dimensional vector elements to specify the vector.
    [Show full text]
  • Derivation of the Two-Dimensional Dot Product
    Derivation of the two-dimensional dot product Content 1. Motivation ....................................................................................................................................... 1 2. Derivation of the dot product in R2 ................................................................................................. 1 2.1. Area of a rectangle ...................................................................................................................... 2 2.2. Area of a right-angled triangle .................................................................................................... 2 2.3. Pythagorean Theorem ................................................................................................................. 2 2.4. Area of a general triangle using the law of cosines ..................................................................... 3 2.5. Derivation of the dot product from the law of cosines ............................................................... 5 3. Geometric Interpretation ................................................................................................................ 5 3.1. Basics ........................................................................................................................................... 5 3.2. Projection .................................................................................................................................... 6 4. Summary.........................................................................................................................................
    [Show full text]
  • Rotation Matrix - Wikipedia, the Free Encyclopedia Page 1 of 22
    Rotation matrix - Wikipedia, the free encyclopedia Page 1 of 22 Rotation matrix From Wikipedia, the free encyclopedia In linear algebra, a rotation matrix is a matrix that is used to perform a rotation in Euclidean space. For example the matrix rotates points in the xy -Cartesian plane counterclockwise through an angle θ about the origin of the Cartesian coordinate system. To perform the rotation, the position of each point must be represented by a column vector v, containing the coordinates of the point. A rotated vector is obtained by using the matrix multiplication Rv (see below for details). In two and three dimensions, rotation matrices are among the simplest algebraic descriptions of rotations, and are used extensively for computations in geometry, physics, and computer graphics. Though most applications involve rotations in two or three dimensions, rotation matrices can be defined for n-dimensional space. Rotation matrices are always square, with real entries. Algebraically, a rotation matrix in n-dimensions is a n × n special orthogonal matrix, i.e. an orthogonal matrix whose determinant is 1: . The set of all rotation matrices forms a group, known as the rotation group or the special orthogonal group. It is a subset of the orthogonal group, which includes reflections and consists of all orthogonal matrices with determinant 1 or -1, and of the special linear group, which includes all volume-preserving transformations and consists of matrices with determinant 1. Contents 1 Rotations in two dimensions 1.1 Non-standard orientation
    [Show full text]
  • Coordinate Transformation
    Coordinate Transformation Coordinate Transformations In this chapter, we explore mappings – where a mapping is a function that "maps" one set to another, usually in a way that preserves at least some of the underlyign geometry of the sets. For example, a 2-dimensional coordinate transformation is a mapping of the form T (u; v) = x (u; v) ; y (u; v) h i The functions x (u; v) and y (u; v) are called the components of the transforma- tion. Moreover, the transformation T maps a set S in the uv-plane to a set T (S) in the xy-plane: If S is a region, then we use the components x = f (u; v) and y = g (u; v) to …nd the image of S under T (u; v) : EXAMPLE 1 Find T (S) when T (u; v) = uv; u2 v2 and S is the unit square in the uv-plane (i.e., S = [0; 1] [0; 1]). Solution: To do so, let’s determine the boundary of T (S) in the xy-plane. We use x = uv and y = u2 v2 to …nd the image of the lines bounding the unit square: Side of Square Result of T (u; v) Image in xy-plane v = 0; u in [0; 1] x = 0; y = u2; u in [0; 1] y-axis for 0 y 1 u = 1; v in [0; 1] x = v; y = 1 v2; v in [0; 1] y = 1 x2; x in[0; 1] v = 1; u in [0; 1] x = u; y = u2 1; u in [0; 1] y = x2 1; x in [0; 1] u = 0; u in [0; 1] x = 0; y = v2; v in [0; 1] y-axis for 1 y 0 1 As a result, T (S) is the region in the xy-plane bounded by x = 0; y = x2 1; and y = 1 x2: Linear transformations are coordinate transformations of the form T (u; v) = au + bv; cu + dv h i where a; b; c; and d are constants.
    [Show full text]